Since the publication of the Age-Related Eye Disease Study (AREDS2) in 2013, the macular pigment carotenoids lutein and zeaxanthin have become well known to both the eye care community and the public. It is a fascinating aspect of evolution that primates have repurposed photoprotective pigments and binding proteins from plants and insects to protect and enhance visual acuity. Moreover, utilization of these plant-derived nutrients has been widely embraced for preventing vision loss from age-related macular degeneration (AMD). More recently, there has been growing awareness that these nutrients can also play a role in improving visual performance in adults. On the other hand, the potential benefits of lutein and zeaxanthin supplementation at very young ages have been underappreciated. In this review, we examine the biochemical mechanisms and supportive data for lutein and zeaxanthin supplementation throughout the lifespan, with particular emphasis on prenatal supplementation. We propose that prenatal nutritional recommendations may aim at improving maternal and infant carotenoid status. Prenatal supplementation with lutein and zeaxanthin might enhance infant visual development and performance and may even prevent retinopathy of prematurity, possibilities that should be examined in future clinical studies.

The LDL receptor-related protein-1 (LRP1) is highly expressed in numerous cell types, and its impairment is associated with obesity, diabetes, and fatty liver disease. However, the mechanisms linking LRP1 to metabolic disease are not completely understood. Here, we compared the metabolic phenotype of C57BL/6J wild type and LRP1 knock-in mice carrying an inactivating mutation in the distal NPxY motif after feeding a low fat (LF) diet or high fat diets with (HFHC) or without (HF) cholesterol supplementation. In response to HF feeding, both groups developed hyperglycemia, hyperinsulinemia, and hyperlipidemia, as well as increased adiposity with adipose tissue inflammation and liver steatosis. However, when animals were fed the HF diet supplemented with cholesterol, the LRP1 NPxY mutation prevents hypercholesterolemia, reduces adipose tissue and brain inflammation, and limits liver progression to steatohepatitis. Nevertheless, insulin signaling is impaired in LRP1 NPxY mutant hepatocytes and this mutation does not protect against HFHC-induced insulin resistance. The selective metabolic improvement observed in HFHC-fed LRP1 NPxY mutant mice is due to an apparent increase of hepatic LDL receptor levels, leading to an elevated rate of plasma lipoprotein clearance and lowering of plasma and hepatic cholesterol levels. The unique metabolic phenotypes displayed by LRP1 NPxY mutant mice in response to HF or HFHC diet feeding indicate an LRP1-cholesterol axis in modulating tissue inflammation. The LRP1 NPxY mutant mouse phenotype differs from phenotypes observed in mice with tissue-specific LRP1 inactivation, thus highlighting the importance of an integrative approach to evaluate how global LRP1 dysfunction contributes to metabolic disease development.

The acyltransferase LCAT mediates FA esterification of plasma cholesterol. In vitro studies have shown that LCAT also FA-esterifies several oxysterols, but in vivo evidence is lacking. Here, we measured both free and FA-esterified forms of sterols in 206 healthy volunteers and 8 individuals with genetic LCAT deficiency, including familial LCAT deficiency (FLD) and fish-eye disease (FED). In the healthy volunteers, the mean values of the ester-to-total molar ratios of the following sterols varied: 4β-hydroxycholesterol (4βHC), 0.38; 5,6α-epoxycholesterol (5,6αEC), 0.46; 5,6β-epoxycholesterol (5,6βEC), 0.51; cholesterol, 0.70; cholestane-3β,5α,6β-triol (CT), 0.70; 7-ketocholesterol (7KC), 0.75; 24S-hydroxycholesterol (24SHC), 0.80; 25-hydroxycholesterol (25HC), 0.81; 27-hydroxycholesterol (27HC), 0.86; and 7α-hydroxycholesterol (7αHC), 0.89. In the individuals with LCAT deficiency, the plasma levels of the FA-esterified forms of cholesterol, 5,6αEC, 5,6βEC, CT, 7αHC, 7KC, 24SHC, 25HC, and 27HC, were significantly lower than those in the healthy volunteers. The individuals with FLD had significantly lower FA-esterified forms of 7αHC, 24SHC, and 27HC than those with FED. It is of note that, even in the three FLD individuals with negligible plasma cholesteryl ester, substantial amounts of the FA-esterified forms of 4βHC, 5,6αEC, 7αHC, 7KC, and 27HC were present. We conclude that LCAT has a major role in the FA esterification of many plasma oxysterols but contributes little to the FA esterification of 4βHC. Substantial FA esterification of 4βHC, 5,6αEC, 7αHC, 7KC, and 27HC is independent of LCAT.

The backbone of all sphingolipids (SLs) is a sphingoid long-chain base (LCB) to which a fatty acid is N-acylated. Considerable variability exists in the chain length and degree of saturation of both of these hydrophobic chains, and recent work has implicated ceramides with different LCBs and N-acyl chains in distinct biological processes; moreover, they may play different roles in disease states and possibly even act as prognostic markers. We now demonstrate that the half-life, or turnover rate, of ceramides containing diverse N-acyl chains is different. By means of a pulse-labeling protocol using stable-isotope, deuterated free fatty acids, and following their incorporation into ceramide and downstream SLs, we show that very-long-chain (VLC) ceramides containing C24:0 or C24:1 fatty acids turn over much more rapidly than long-chain (LC) ceramides containing C16:0 or C18:0 fatty acids due to the more rapid metabolism of the former into VLC sphingomyelin and VLC hexosylceramide. In contrast, d16:1 and d18:1 ceramides show similar rates of turnover, indicating that the length of the sphingoid LCB does not influence the flux of ceramides through the biosynthetic pathway. Together, these data demonstrate that the N-acyl chain length of SLs may not only affect membrane biophysical properties but also influence the rate of metabolism of SLs so as to regulate their levels and perhaps their biological functions.

Bile acids (BAs) have been established as ubiquitous regulatory molecules implicated in a large variety of healthy and pathological processes. However, the scope of BA heterogeneity is often underrepresented in current literature. This is due in part to inadequate detection methods, which fail to distinguish the individual constituents of the BA pool. Thus, the primary aim of this study was to develop a method that would allow the simultaneous analysis of specific C24 BA species, and to apply that method to biological systems of interest. Herein, we describe the generation and validation of an LC-MS/MS assay for quantification of numerous BAs in a variety of cell systems and relevant biofluids and tissue. These studies included the first baseline level assessment for planar BAs, including allocholic acid, in cell lines, biofluids, and tissue in a nonhuman primate (NHP) laboratory animal, Macaca mulatta, in healthy conditions. These results indicate that immortalized cell lines make poor models for the study of BA synthesis and metabolism, whereas human primary hepatocytes represent a promising alternative model system. We also characterized the BA pool of M. mulatta in detail. Our results support the use of NHP models for the study of BA metabolism and pathology in lieu of murine models. Moreover, the method developed here can be applied to the study of common and planar C24 BA species in other systems.

Oxylipins are potent lipid mediators involved in a variety of physiological processes. Their profiling has the potential to provide a wealth of information regarding human health and disease and is a promising technology for translation into clinical applications. However, results generated by independent groups are rarely comparable, which increases the need for the implementation of internationally agreed upon protocols. We performed an interlaboratory comparison for the MS-based quantitative analysis of total oxylipins. Five independent laboratories assessed the technical variability and comparability of 133 oxylipins using a harmonized and standardized protocol, common biological materials (i.e., seven quality control plasmas), standard calibration series, and analytical methods. The quantitative analysis was based on a standard calibration series with isotopically labeled internal standards. Using the standardized protocol, the technical variance was within ±15% for 73% of oxylipins; however, most epoxy fatty acids were identified as critical analytes due to high variabilities in concentrations. The comparability of concentrations determined by the laboratories was examined using consensus value estimates and unsupervised/supervised multivariate analysis (i.e., principal component analysis and partial least squares discriminant analysis). Interlaboratory variability was limited and did not interfere with our ability to distinguish the different plasmas. Moreover, all laboratories were able to identify similar differences between plasmas. In summary, we show that by using a standardized protocol for sample preparation, low technical variability can be achieved. Harmonization of all oxylipin extraction and analysis steps led to reliable, reproducible, and comparable oxylipin concentrations in independent laboratories, allowing the generation of biologically meaningful oxylipin patterns.

Lipopolysaccharide (LPS) is a key player for innate immunity activation. It is therefore a prime target for sepsis treatment, as antibiotics are not sufficient to improve outcome during septic shock. An extracorporeal removal method by polymyxin (PMX) B direct hemoperfusion (PMX-DHP) is used in Japan, but recent trials failed to show a significant lowering of circulating LPS levels after PMX-DHP therapy. PMX-DHP has a direct effect on LPS molecules. However, LPS is not present in a free form in the circulation, as it is mainly carried by lipoproteins, including LDLs. Lipoproteins are critical for physiological LPS clearance, as LPSs are carried by LDLs to the liver for elimination. We hypothesized that LDL apheresis could be an alternate method for LPS removal. First, we demonstrated in vitro that LDL apheresis microbeads are almost as efficient as PMX beads to reduce LPS concentration in LPS-spiked human plasma, whereas it is not active in PBS. We found that PMX was also adsorbing lipoproteins, although less specifically. Then, we found that endogenous LPS of patients treated by LDL apheresis for familial hypercholesterolemia is also removed during their LDL apheresis sessions, with both electrostatic-based devices and filtration devices. Finally, LPS circulating in the plasma of septic shock and severe sepsis patients with gram-negative bacteremia was also removed in vitro by LDL adsorption. Overall, these results underline the importance of lipoproteins for LPS clearance, making them a prime target to study and treat endotoxemia-related conditions.

Lipins are eukaryotic proteins with functions in lipid synthesis and the homeostatic control of energy balance. They execute these functions by acting as phosphatidate phosphatase enzymes in the cytoplasm and by changing gene expression after translocation into the cell nucleus, in particular under fasting conditions. Here, we asked whether nuclear translocation and the enzymatic activity of Drosophila Lipin serve essential functions and how gene expression changes, under both fed and fasting conditions, when nuclear translocation is impaired. To address these questions, we created a Lipin null mutant, a mutant expressing Lipin lacking a nuclear localization signal (Lipin^NLS), and a mutant expressing enzymatically dead Lipin. Our data support the conclusion that the enzymatic but not nuclear gene regulatory activity of Lipin is essential for survival. Notably, adult Lipin^NLS flies were not only viable but also exhibited improved life expectancy. In contrast, they were highly susceptible to starvation. Both the improved life expectancy in the fed state and the decreased survival in the fasting state correlated with changes in metabolic gene expression. Moreover, increased life expectancy of fed flies was associated with a decreased metabolic rate. Interestingly, in addition to metabolic genes, genes involved in feeding behavior and the immune response were misregulated in Lipin^NLS flies. Altogether, our data suggest that the nuclear activity of Lipin influences the genomic response to nutrient availability with effects on life expectancy and starvation resistance. Thus, nutritional or therapeutic approaches that aim at lowering nuclear translocation of lipins in humans may be worth exploring.

Phospholipids, including ether phospholipids, are composed of numerous isomeric and isobaric species that have the same backbone and acyl chains. This structural resemblance results in similar fragmentation patterns by collision-induced dissociation of phospholipids regardless of class, yielding complicated MS/MS spectra when isobaric species are analyzed together. Furthermore, the presence of isobaric species can lead to misassignment of species when made solely based on their molecular weights. In this study, we used normal-phase HPLC for ESI-MS/MS analysis of phospholipids from bovine heart mitochondria. Class separation by HPLC eliminates chances for misidentification of isobaric species from different classes of phospholipids. Chromatography yields simple MS/MS spectra without interference from isobaric species, allowing clear identification of peaks corresponding to fragmented ions containing monoacylglycerol backbone derived from losing one acyl chain. Using these fragmented ions, we characterized individual and isomeric species in each class of mitochondrial phospholipids, including unusual species, such as PS, containing an ether linkage and species containing odd-numbered acyl chains in cardiolipin, PS, PI, and PG. We also characterized monolysocardiolipin and dilysocardiolipin, the least abundant but nevertheless important mitochondrial phospholipids. The results clearly show the power of HPLC-MS/MS for identification and characterization of phospholipids, including minor species.

HMG-CoA reductase (Hmgcr) is the rate-limiting enzyme in the mevalonate pathway and is inhibited by statins. In addition to cholesterol, Hmgcr activity is also required for synthesizing nonsterol isoprenoids, such as dolichol, ubiquinone, and farnesylated and geranylgeranylated proteins. Here, we investigated the effects of Hmgcr inhibition on nonsterol isoprenoids in the liver. We have generated new genetic models to acutely delete genes in the mevalonate pathway in the liver using AAV-mediated delivery of Cre-recombinase (AAV-Cre) or CRISPR/Cas9 (AAV-CRISPR). The genetic deletion of Hmgcr by AAV-Cre resulted in extensive hepatocyte apoptosis and compensatory liver regeneration. At the biochemical level, we observed decreased levels of sterols and depletion of the nonsterol isoprenoids, dolichol and ubiquinone. At the cellular level, Hmgcr-null hepatocytes showed ER stress and impaired N-glycosylation. We further hypothesized that the depletion of dolichol, essential for N-glycosylation, could be responsible for ER stress. Using AAV-CRISPR, we somatically disrupted dehydrodolichyl diphosphate synthase subunit (Dhdds), encoding a branch point enzyme required for dolichol biosynthesis. Dhdds-null livers showed ER stress and impaired N-glycosylation, along with apoptosis and regeneration. Finally, the combined deletion of Hmgcr and Dhdds synergistically exacerbated hepatocyte ER stress. Our data show a critical role for mevalonate-derived dolichol in the liver and suggest that dolichol depletion is at least partially responsible for ER stress and apoptosis upon potent Hmgcr inhibition.

A comprehensive and standardized system to report lipid structures analyzed by MS is essential for the communication and storage of lipidomics data. Herein, an update on both the LIPID MAPS classification system and shorthand notation of lipid structures is presented for lipid categories Fatty Acyls (FA), Glycerolipids (GL), Glycerophospholipids (GP), Sphingolipids (SP), and Sterols (ST). With its major changes, i.e., annotation of ring double bond equivalents and number of oxygens, the updated shorthand notation facilitates reporting of newly delineated oxygenated lipid species as well. For standardized reporting in lipidomics, the hierarchical architecture of shorthand notation reflects the diverse structural resolution powers provided by mass spectrometric assays. Moreover, shorthand notation is expanded beyond mammalian phyla to lipids from plant and yeast phyla. Finally, annotation of atoms is included for the use of stable isotope-labeled compounds in metabolic labeling experiments or as internal standards. This update on lipid classification, nomenclature, and shorthand annotation for lipid mass spectra is considered a standard for lipid data presentation.

S-Acylation, a reversible post-translational lipid modification of proteins, controls the properties and function of various proteins, including ion channels. Large conductance Ca2+-activated potassium (BK) channels are S-acylated at two sites that impart distinct functional effects. Whereas the enzymes that attach lipid groups are known, the enzymes mediating lipid removal (i.e. deacylation) are largely unknown. Here, McClafferty et al. identify two enzymes, ABHD17a and ABHD17c, that excise BK channel lipid groups with remarkable precision. These findings lend insights into mechanisms that orchestrate the (de)acylation that fine-tunes ion channel function in physiology and disease.

S-Acylation, the reversible post-translational lipid modification of proteins, is an important mechanism to control the properties and function of ion channels and other polytopic transmembrane proteins. However, although increasing evidence reveals the role of diverse acyl protein transferases (zDHHC) in controlling ion channel S-acylation, the acyl protein thioesterases that control ion channel deacylation are very poorly defined. Here we show that ABHD17a (α/β-hydrolase domain-containing protein 17a) deacylates the stress-regulated exon domain of large conductance voltage- and calcium-activated potassium (BK) channels inhibiting channel activity independently of effects on channel surface expression. Importantly, ABHD17a deacylates BK channels in a site-specific manner because it has no effect on the S-acylated S0–S1 domain conserved in all BK channels that controls membrane trafficking and is deacylated by the acyl protein thioesterase Lypla1. Thus, distinct S-acylated domains in the same polytopic transmembrane protein can be regulated by different acyl protein thioesterases revealing mechanisms for generating both specificity and diversity for these important enzymes to control the properties and functions of ion channels.

The divalent anion sodium symporter (DASS) family (SLC13) plays critical roles in metabolic homeostasis, influencing many processes, including fatty acid synthesis, insulin resistance, and adiposity. DASS transporters catalyze the Na+-driven concentrative uptake of Krebs cycle intermediates and sulfate into cells; disrupting their function can protect against age-related metabolic diseases and can extend lifespan. An inward-facing crystal structure and an outward-facing model of a bacterial DASS family member, VcINDY from Vibrio cholerae, predict an elevator-like transport mechanism involving a large rigid body movement of the substrate-binding site. How substrate binding influences the conformational state of VcINDY is currently unknown. Here, we probe the interaction between substrate binding and protein conformation by monitoring substrate-induced solvent accessibility changes of broadly distributed positions in VcINDY using a site-specific alkylation strategy. Our findings reveal that accessibility to all positions tested is modulated by the presence of substrates, with the majority becoming less accessible in the presence of saturating concentrations of both Na+ and succinate. We also observe separable effects of Na+ and succinate binding at several positions suggesting distinct effects of the two substrates. Furthermore, accessibility changes to a solely succinate-sensitive position suggests that substrate binding is a low-affinity, ordered process. Mapping these accessibility changes onto the structures of VcINDY suggests that Na+ binding drives the transporter into an as-yet-unidentified conformational state, involving rearrangement of the substrate-binding site–associated re-entrant hairpin loops. These findings provide insight into the mechanism of VcINDY, which is currently the only structurally characterized representative of the entire DASS family.

The plasma membrane of a cell is characterized by an asymmetric distribution of lipid species across the exofacial and cytofacial aspects of the bilayer. Regulation of membrane asymmetry is a fundamental characteristic of membrane biology and is crucial for signal transduction, vesicle transport, and cell division. The type IV family of P-ATPases, or P4-ATPases, establishes membrane asymmetry by selection and transfer of a subset of membrane lipids from the lumenal or exofacial leaflet to the cytofacial aspect of the bilayer. It is unclear how P4-ATPases sort through the spectrum of membrane lipids to identify their desired substrate(s) and how the membrane environment modulates this activity. Therefore, we tested how the yeast plasma membrane P4-ATPase, Dnf2, responds to changes in membrane composition induced by perturbation of endogenous lipid biosynthetic pathways or exogenous application of lipid. The primary substrates of Dnf2 are glucosylceramide (GlcCer) and phosphatidylcholine (PC, or their lyso-lipid derivatives), and we find that these substrates compete with each other for transport. Acutely inhibiting sphingolipid synthesis using myriocin attenuates transport of exogenously applied GlcCer without perturbing PC transport. Deletion of genes controlling later steps of glycosphingolipid production also perturb GlcCer transport to a greater extent than PC transport. In contrast, perturbation of ergosterol biosynthesis reduces PC and GlcCer transport equivalently. Surprisingly, application of lipids that are poor transport substrates differentially affects PC and GlcCer transport by Dnf2, thus altering substrate preference. Our data indicate that Dnf2 exhibits exquisite sensitivity to the membrane composition, thus providing feedback onto the function of the P4-ATPases.

Mass spectrometry has become an indispensable tool for the characterization of glycosylation across biological systems. Our ability to generate rich fragmentation of glycopeptides has dramatically improved over the last decade yet our informatic approaches still lag behind. Although glycoproteomic informatics approaches using glycan databases have attracted considerable attention, database independent approaches have not. This has significantly limited high throughput studies of unusual or atypical glycosylation events such as those observed in bacteria. As such, computational approaches to examine bacterial glycosylation and identify chemically diverse glycans are desperately needed. Here we describe the use of wide-tolerance (up to 2000 Da) open searching as a means to rapidly examine bacterial glycoproteomes. We benchmarked this approach using N-linked glycopeptides of Campylobacter fetus subsp. fetus as well as O-linked glycopeptides of Acinetobacter baumannii and Burkholderia cenocepacia revealing glycopeptides modified with a range of glycans can be readily identified without defining the glycan masses before database searching. Using this approach, we demonstrate how wide tolerance searching can be used to compare glycan use across bacterial species by examining the glycoproteomes of eight Burkholderia species (B. pseudomallei; B. multivorans; B. dolosa; B. humptydooensis; B. ubonensis, B. anthina; B. diffusa; B. pseudomultivorans). Finally, we demonstrate how open searching enables the identification of low frequency glycoforms based on shared modified peptides sequences. Combined, these results show that open searching is a robust computational approach for the determination of glycan diversity within bacterial proteomes.

Influenza A virus (IAV) mutates rapidly, resulting in antigenic drift and poor year-to-year vaccine effectiveness. One challenge in designing effective vaccines is that genetic mutations frequently cause amino acid variations in IAV envelope protein hemagglutinin (HA) that create new N-glycosylation sequons; resulting N-glycans cause antigenic shielding, allowing viral escape from adaptive immune responses. Vaccine candidate strain selection currently involves correlating antigenicity with HA protein sequence among circulating strains, but quantitative comparison of site-specific glycosylation information may likely improve the ability to design vaccines with broader effectiveness against evolving strains. However, there is poor understanding of the influence of glycosylation on immunodominance, antigenicity, and immunogenicity of HA, and there are no well-tested methods for comparing glycosylation similarity among virus samples. Here, we present a method for statistically rigorous quantification of similarity between two related virus strains that considers the presence and abundance of glycopeptide glycoforms. We demonstrate the strength of our approach by determining that there was a quantifiable difference in glycosylation at the protein level between WT IAV HA from A/Switzerland/9715293/2013 (SWZ13) and a mutant strain of SWZ13, even though no N-glycosylation sequons were changed. We determined site-specifically that WT and mutant HA have varying similarity at the glycosylation sites of the head domain, reflecting competing pressures to evade host immune response while retaining viral fitness. To our knowledge, our results are the first to quantify changes in glycosylation state that occur in related proteins of considerable glycan heterogeneity. Our results provide a method for understanding how changes in glycosylation state are correlated with variations in protein sequence, which is necessary for improving IAV vaccine strain selection. Understanding glycosylation will be especially important as we find new expression vectors for vaccine production, as glycosylation state depends greatly on the host species.

Despite the continued analysis of HDAC inhibitors in clinical trials, the heterogeneous nature of the protein complexes they target limits our understanding of the beneficial and off-target effects associated with their application. Among the many HDAC protein complexes found within the cell, Sin3 complexes are conserved from yeast to humans and likely play important roles as regulators of transcriptional activity. The presence of two Sin3 paralogs in humans, SIN3A and SIN3B, may result in a heterogeneous population of Sin3 complexes and contributes to our poor understanding of the functional attributes of these complexes. Here, we profile the interaction networks of SIN3A and SIN3B to gain insight into complex composition and organization. In accordance with existing data, we show that Sin3 paralog identity influences complex composition. Additionally, chemical cross-linking MS identifies domains that mediate interactions between Sin3 proteins and binding partners. The characterization of rare SIN3B proteoforms provides additional evidence for the existence of conserved and divergent elements within human Sin3 proteins. Together, these findings shed light on both the shared and divergent properties of human Sin3 proteins and highlight the heterogeneous nature of the complexes they organize.

Kir2.1, a strong inward rectifier potassium channel encoded by the KCNJ2 gene, is a key regulator of the resting membrane potential of the cardiomyocyte and plays an important role in controlling ventricular excitation and action potential duration in the human heart. Mutations in KCNJ2 result in inheritable cardiac diseases in humans, e.g. the type-1 Andersen-Tawil syndrome (ATS1). Understanding the molecular mechanisms that govern the regulation of inward rectifier potassium currents by Kir2.1 in both normal and disease contexts should help uncover novel targets for therapeutic intervention in ATS1 and other Kir2.1-associated channelopathies. The information available to date on protein-protein interactions involving Kir2.1 channels remains limited. Additional efforts are necessary to provide a comprehensive map of the Kir2.1 interactome. Here we describe the generation of a comprehensive map of the Kir2.1 interactome using the proximity-labeling approach BioID. Most of the 218 high-confidence Kir2.1 channel interactions we identified are novel and encompass various molecular mechanisms of Kir2.1 function, ranging from intracellular trafficking to cross-talk with the insulin-like growth factor receptor signaling pathway, as well as lysosomal degradation. Our map also explores the variations in the interactome profiles of Kir2.1^WT versus Kir2.1^314-315, a trafficking deficient ATS1 mutant, thus uncovering molecular mechanisms whose malfunctions may underlie ATS1 disease. Finally, using patch-clamp analysis, we validate the functional relevance of PKP4, one of our top BioID interactors, to the modulation of Kir2.1-controlled inward rectifier potassium currents. Our results validate the power of our BioID approach in identifying functionally relevant Kir2.1 interactors and underline the value of our Kir2.1 interactome as a repository for numerous novel biological hypotheses on Kir2.1 and Kir2.1-associated diseases.

Synaptic transmission leading to release of neurotransmitters in the nervous system is a fast and highly dynamic process. Previously, protein interaction and phosphorylation have been thought to be the main regulators of synaptic transmission. Here we show that sialylation of N-linked glycosylation is a novel potential modulator of neurotransmitter release mechanisms by investigating depolarization-dependent changes of formerly sialylated N-linked glycopeptides. We suggest that negatively charged sialic acids can be modulated, similarly to phosphorylation, by the action of sialyltransferases and sialidases thereby changing local structure and function of membrane glycoproteins. We characterized site-specific alteration in sialylation on N-linked glycoproteins in isolated rat nerve terminals after brief depolarization using quantitative sialiomics. We identified 1965 formerly sialylated N-linked glycosites in synaptic proteins and found that the abundances of 430 glycosites changed after 5 s depolarization. We observed changes on essential synaptic proteins such as synaptic vesicle proteins, ion channels and transporters, neurotransmitter receptors and cell adhesion molecules. This study is to our knowledge the first to describe ultra-fast site-specific modulation of the sialiome after brief stimulation of a biological system.

Tandem mass tag (TMT) is a multiplexing technology widely-used in proteomic research. It enables relative quantification of proteins from multiple biological samples in a single MS run with high efficiency and high throughput. However, experiments often require more biological replicates or conditions than can be accommodated by a single run, and involve multiple TMT mixtures and multiple runs. Such larger-scale experiments combine sources of biological and technical variation in patterns that are complex, unique to TMT-based workflows, and challenging for the downstream statistical analysis. These patterns cannot be adequately characterized by statistical methods designed for other technologies, such as label-free proteomics or transcriptomics. This manuscript proposes a general statistical approach for relative protein quantification in MS- based experiments with TMT labeling. It is applicable to experiments with multiple conditions, multiple biological replicate runs and multiple technical replicate runs, and unbalanced designs. It is based on a flexible family of linear mixed-effects models that handle complex patterns of technical artifacts and missing values. The approach is implemented in MSstatsTMT, a freely available open-source R/Bioconductor package compatible with data processing tools such as Proteome Discoverer, MaxQuant, OpenMS, and SpectroMine. Evaluation on a controlled mixture, simulated datasets, and three biological investigations with diverse designs demonstrated that MSstatsTMT balanced the sensitivity and the specificity of detecting differentially abundant proteins, in large-scale experiments with multiple biological mixtures.

Natural compounds that can stimulate salivary secretion are of interest in developing treatments for xerostomia, the perception of a dry mouth, that affects between 10 and 30% of the adult and elderly population. Chemesthetic transient receptor potential (TRP) channels are expressed in the surface of the oral mucosa. The TRPV1 agonists capsaicin and piperine have been shown to increase salivary flow when introduced into the oral cavity but the sialogogic properties of other TRP channel agonists have not been investigated. In this study we have determined the influence of different TRP channel agonists on the flow and protein composition of saliva. Mouth rinsing with the TRPV1 agonist nonivamide or menthol, a TRPM8 agonist, increased whole mouth saliva (WMS) flow and total protein secretion compared with unstimulated saliva, the vehicle control mouth rinse or cinnamaldehyde, a TRPA1 agonist. Nonivamide also increased the flow of labial minor gland saliva but parotid saliva flow rate was not increased. The influence of TRP channel agonists on the composition and function of the salivary proteome was investigated using a multi-batch quantitative MS method novel to salivary proteomics. Inter-personal and inter-mouth rinse variation was observed in the secreted proteomes and, using a novel bioinformatics method, inter-day variation was identified with some of the mouth rinses. Significant changes in specific salivary proteins were identified after all mouth rinses. In the case of nonivamide, these changes were attributed to functional shifts in the WMS secreted, primarily the over representation of salivary and nonsalivary cystatins which was confirmed by immunoassay. This study provides new evidence of the impact of TRP channel agonists on the salivary proteome and the stimulation of salivary secretion by a TRPM8 channel agonist, which suggests that TRP channel agonists are potential candidates for developing treatments for sufferers of xerostomia.

The neuronal basis of complex social behavior is still poorly understood. In honeybees, reproductive investment decisions are made at the colony-level. Queens develop from female-destined larvae that receive alloparental care from nurse bees in the form of ad-libitum royal jelly (RJ) secretions. Typically, the number of raised new queens is limited but genetic breeding of "royal jelly bees" (RJBs) for enhanced RJ production over decades has led to a dramatic increase of reproductive investment in queens. Here, we compare RJBs to unselected Italian bees (ITBs) to investigate how their cognitive processing of larval signals in the mushroom bodies (MBs) and antennal lobes (ALs) may contribute to their behavioral differences. A cross-fostering experiment confirms that the RJB syndrome is mainly due to a shift in nurse bee alloparental care behavior. Using olfactory conditioning of the proboscis extension reflex, we show that the RJB nurses spontaneously respond more often to larval odors compared with ITB nurses but their subsequent learning occurs at similar rates. These phenotypic findings are corroborated by our demonstration that the proteome of the brain, particularly of the ALs differs between RJBs and ITBs. Notably, in the ALs of RJB newly emerged bees and nurses compared with ITBs, processes of energy and nutrient metabolism, signal transduction are up-regulated, priming the ALs for receiving and processing the brood signals from the antennae. Moreover, highly abundant major royal jelly proteins and hexamerins in RJBs compared with ITBs during early life when the nervous system still develops suggest crucial new neurobiological roles for these well-characterized proteins. Altogether, our findings reveal that RJBs have evolved a strong olfactory response to larvae, enabled by numerous neurophysiological adaptations that increase the nurse bees' alloparental care behavior.

Glutathionylation is an important posttranslational modification that protects proteins from further oxidative damage as well as influencing protein structure and activity. In the present study, we demonstrate that the cysteine-42 residue in protein arginine N-methyltransferase 5 (PRMT5) is glutathionylated in aged mice or in cells that have been exposed to oxidative stress. Deglutathionylation of this protein is catalyzed by glutaredoxin-1 (Grx1). Using mutagenesis and subsequent biochemical analyses, we show that glutathionylation decreased the binding affinity of PRMT5 with methylosome protein-50 (MEP50) and reduced the methyltransferase activity of PRMT5. Furthermore, overexpression of PRMT5-C42A mutant caused a significant increase in histone methylation in HEK293T and A549 cells and promoted cell growth, whereas overexpression of the PRMT5-C42D mutant, a mimic of glutathionylated PRMT5, inhibited cell proliferation. Taken together, our results demonstrate a new mechanism of regulation of PRMT5 methyltransferases activity and suggest that PRMT5 glutathionylation is partly responsible for reactive oxygen species-mediated cell growth inhibition.

Studies in the yeast Saccharomyces cerevisiae have helped define mechanisms underlying the activity of the ubiquitin–proteasome system (UPS), uncover the proteasome assembly pathway, and link the UPS to the maintenance of cellular homeostasis. However, the spectrum of UPS substrates is incompletely defined, even though multiple techniques—including MS—have been used. Therefore, we developed a substrate trapping proteomics workflow to identify previously unknown UPS substrates. We first generated a yeast strain with an epitope tagged proteasome subunit to which a proteasome inhibitor could be applied. Parallel experiments utilized inhibitor insensitive strains or strains lacking the tagged subunit. After affinity isolation, enriched proteins were resolved, in-gel digested, and analyzed by high resolution liquid chromatography-tandem MS. A total of 149 proteasome partners were identified, including all 33 proteasome subunits. When we next compared data between inhibitor sensitive and resistant cells, 27 proteasome partners were significantly enriched. Among these proteins were known UPS substrates and proteins that escort ubiquitinated substrates to the proteasome. We also detected Erg25 as a high-confidence partner. Erg25 is a methyl oxidase that converts dimethylzymosterol to zymosterol, a precursor of the plasma membrane sterol, ergosterol. Because Erg25 is a resident of the endoplasmic reticulum (ER) and had not previously been directly characterized as a UPS substrate, we asked whether Erg25 is a target of the ER associated degradation (ERAD) pathway, which most commonly mediates proteasome-dependent destruction of aberrant proteins. As anticipated, Erg25 was ubiquitinated and associated with stalled proteasomes. Further, Erg25 degradation depended on ERAD-associated ubiquitin ligases and was regulated by sterol synthesis. These data expand the cohort of lipid biosynthetic enzymes targeted for ERAD, highlight the role of the UPS in maintaining ER function, and provide a novel tool to uncover other UPS substrates via manipulations of our engineered strain.

Renal Cell Carcinoma (RCC) is one of the most commonly diagnosed cancers worldwide with research efforts dramatically improving understanding of the biology of the disease. To investigate the role of the immune system in treatment-naïve clear cell Renal Cell Carcinoma (ccRCC), we interrogated the immune infiltrate in patient-matched ccRCC tumor samples, benign normal adjacent tissue (NAT) and peripheral blood mononuclear cells (PBMCs isolated from whole blood, focusing our attention on the myeloid cell infiltrate. Using flow cytometric, MS, and ExCYT analysis, we discovered unique myeloid populations in PBMCs across patient samples. Furthermore, normal adjacent tissues and ccRCC tissues contained numerous myeloid populations with a unique signature for both tissues. Enrichment of the immune cell (CD45⁺) fraction and subsequent gene expression analysis revealed a number of myeloid-related genes that were differentially expressed. These data provide evidence, for the first time, of an immunosuppressive and pro-tumorigenic role of myeloid cells in early, clinically localized ccRCC. The identification of a number of immune proteins for therapeutic targeting provides a rationale for investigation into the potential efficacy of earlier intervention with single-agent or combination immunotherapy for ccRCC.

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.

Cross-linking MS (XL-MS) has been recognized as an effective source of information about protein structures and interactions. In contrast to regular peptide identification, XL-MS has to deal with a quadratic search space, where peptides from every protein could potentially be cross-linked to any other protein. To cope with this search space, most tools apply different heuristics for search space reduction. We introduce a new open-source XL-MS database search algorithm, OpenPepXL, which offers increased sensitivity compared with other tools. OpenPepXL searches the full search space of an XL-MS experiment without using heuristics to reduce it. Because of efficient data structures and built-in parallelization OpenPepXL achieves excellent runtimes and can also be deployed on large compute clusters and cloud services while maintaining a slim memory footprint. We compared OpenPepXL to several other commonly used tools for identification of noncleavable labeled and label-free cross-linkers on a diverse set of XL-MS experiments. In our first comparison, we used a data set from a fraction of a cell lysate with a protein database of 128 targets and 128 decoys. At 5% FDR, OpenPepXL finds from 7% to over 50% more unique residue pairs (URPs) than other tools. On data sets with available high-resolution structures for cross-link validation OpenPepXL reports from 7% to over 40% more structurally validated URPs than other tools. Additionally, we used a synthetic peptide data set that allows objective validation of cross-links without relying on structural information and found that OpenPepXL reports at least 12% more validated URPs than other tools. It has been built as part of the OpenMS suite of tools and supports Windows, macOS, and Linux operating systems. OpenPepXL also supports the MzIdentML 1.2 format for XL-MS identification results. It is freely available under a three-clause BSD license at https://openms.org/openpepxl.

Over the past decade, modern methods of MS (MS) have emerged that allow reliable, fast and cost-effective identification of pathogenic microorganisms. Although MALDI-TOF MS has already revolutionized the way microorganisms are identified, recent years have witnessed also substantial progress in the development of liquid chromatography (LC)-MS based proteomics for microbiological applications. For example, LC-tandem MS (LC-MS²) has been proposed for microbial characterization by means of multiple discriminative peptides that enable identification at the species, or sometimes at the strain level. However, such investigations can be laborious and time-consuming, especially if the experimental LC-MS² data are tested against sequence databases covering a broad panel of different microbiological taxa. In this proof of concept study, we present an alternative bottom-up proteomics method for microbial identification. The proposed approach involves efficient extraction of proteins from cultivated microbial cells, digestion by trypsin and LC–MS measurements. Peptide masses are then extracted from MS¹ data and systematically tested against an in silico library of all possible peptide mass data compiled in-house. The library has been computed from the UniProt Knowledgebase covering Swiss-Prot and TrEMBL databases and comprises more than 12,000 strain-specific in silico profiles, each containing tens of thousands of peptide mass entries. Identification analysis involves computation of score values derived from correlation coefficients between experimental and strain-specific in silico peptide mass profiles and compilation of score ranking lists. The taxonomic positions of the microbial samples are then determined by using the best-matching database entries. The suggested method is computationally efficient – less than 2 mins per sample - and has been successfully tested by a test set of 39 LC-MS¹ peak lists obtained from 19 different microbial pathogens. The proposed method is rapid, simple and automatable and we foresee wide application potential for future microbiological applications.

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.

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.

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.

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.

O-GlcNAcylation, the addition of a single N-acetylglucosamine residue to serine and threonine residues of cytoplasmic, nuclear, or mitochondrial proteins, is a widespread regulatory post-translational modification. It is involved in response to nutritional status and stress and its dysregulation is associated with diseases ranging from Alzheimer’s to diabetes.  While the modification was first detected over thirty-five years ago, research into the function of O-GlcNAcylation has accelerated dramatically in the last ten years due to the development of new enrichment and mass spectrometry techniques that facilitate its analysis.  This article summarizes methods for O-GlcNAc enrichment, key mass spectrometry instrumentation advancements, particularly those that allow modification site localization, and software tools that allow analysis of data from O-GlcNAc modified peptides.

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.

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.

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.

We developed a simple and rapid method to enrich protein N-terminal peptides, in which the protease TrypN is first employed to generate protein N-terminal peptides without Lys or Arg and internal peptides with two positive charges at their N-termini, and then the N-terminal peptides with or without N-acetylation are separated from the internal peptides by strong cation exchange chromatography according to a retention model based on the charge/orientation of peptides. This approach was applied to 20 μg of human HEK293T cell lysate proteins to profile the N-terminal proteome. On average, 1,550 acetylated and 200 unmodified protein N-terminal peptides were successfully identified in a single LC/MS/MS run with less than 3% contamination with internal peptides, even when we accepted only canonical protein N-termini registered in the Swiss-Prot database. Since this method involves only two steps, protein digestion and chromatographic separation, without the need for tedious chemical reactions, it should be useful for comprehensive profiling of protein N-termini, including proteoforms with neo-N-termini.

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.

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.

Histone post-translational modifications (PTMs) are one of the main mechanisms of epigenetic regulation. Dysregulation of histone PTMs leads to many human diseases, such as cancer. Due to its high-throughput, accuracy, and flexibility, mass spectrometry (MS) has emerged as a powerful tool in the epigenetic histone modification field, allowing the comprehensive and unbiased analysis of histone PTMs and chromatin-associated factors. Coupled with various techniques from molecular biology, biochemistry, chemical biology and biophysics, MS has been employed to characterize distinct aspects of histone PTMs in the epigenetic regulation of chromatin functions. In this review we will describe advancements in the field of MS that have facilitated the analysis of histone PTMs and chromatin biology.  

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.

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

Open searching has proven to be an effective strategy for identifying both known and unknown modifications in shotgun proteomics experiments. Rather than being limited to a small set of user-specified modifications, open searches identify peptides with any mass shift that may correspond to a single modification or a combination of several modifications. Here we present PTM-Shepherd, a bioinformatics tool that automates characterization of PTM profiles detected in open searches based on attributes such as amino acid localization, fragmentation spectra similarity, retention time shifts, and relative modification rates. PTM-Shepherd can also perform multi-experiment comparisons for studying changes in modification profiles, e.g. in data generated in different laboratories or under different conditions. We demonstrate how PTM-Shepherd improves the analysis of data from formalin-fixed paraffin-embedded samples, detects extreme underalkylation of cysteine in some datasets, discovers an artefactual modification introduced during peptide synthesis, and uncovers site-specific biases in sample preparation artifacts in a multi-center proteomics profiling study.

Gonadal soma-derived factor (gsdf) has been demonstrated to be essential for testicular differentiation in medaka (Oryzias latipes). To understand the protein dynamics of Gsdf in spermatogenesis regulation, we used a His-tag "pull-down" assay coupled with shotgun LC-MS/MS to identify a group of potential interacting partners for Gsdf, which included cytoplasmic dynein light chain 2, eukaryotic polypeptide elongation factor 1 alpha (eEF1α), and actin filaments in mature medaka testis. As for the interaction with TGFβ-dynein being critical for spermatogonial division in Drosophila melanogaster, the physical interactions of Gsdf-dynein and Gsdf-eEF1α were identified through a yeast 2-hybrid (Y2H) screening of an adult testis cDNA library using Gsdf as bait, which were verified by a paired Y2H assay. Co-immunoprecipitation of Gsdf and eEF1α was defined in adult testes as supporting the requirement of a Gsdf and eEF1α interaction in testis development. Proteomics analysis (data are available via ProteomeXchange with identifier PXD022153) and ultrastrutural observations showed that Gsdf deficiency activated eEF1α-mediated protein synthesis and ribosomal biogenesis, which in turn led to the differentiation of undifferentiated germ cells. Thus, our results provide a framework and new insight into the coordination of a Gsdf (TGFβ and eEF1α complex in the basic processes of germ cell proliferation, transcriptional and translational control of sexual RNA which may be fundamentally conserved across phyla during sexual differentiation.

The goal of clinical proteomics is to identify, quantify, and characterize proteins in body fluids or tissue to assist diagnosis, prognosis, and treatment of patients. In this way, it is similar to more mature omics technologies, such as genomics, that are increasingly applied in biomedicine. We argue that, similar to those fields, proteomics also faces ethical issues related to the kinds of information that is inherently obtained through sample measurement, although their acquisition was not the primary purpose. Specifically, we demonstrate the potential to identify individuals both by their characteristic, individual-specific protein levels and by variant peptides reporting on coding single nucleotide polymorphisms. Furthermore, it is in the nature of blood plasma proteomics profiling that it broadly reports on the health status of an individual – beyond the disease under investigation. Finally, we show that private and potentially sensitive information, such as ethnicity and pregnancy status, can increasingly be derived from proteomics data. Although this is potentially valuable not only to the individual, but also for biomedical research, it raises ethical questions similar to the incidental findings obtained through other omics technologies. We here introduce the necessity of - and argue for the desirability for - ethical and human rights-related issues to be discussed within the proteomics community. Those thoughts are more fully developed in our accompanying manuscript. Appreciation and discussion of ethical aspects of proteomic research will allow for deeper, better-informed, more diverse, and, most importantly, wiser guidelines for clinical proteomics.

To identify novel autoantibodies of Takayasu arteritis (TAK) using HuProt array-based approach. A two-phase approach was adopted. In Phase I, serum samples collected from 40 TAK patients, 15 autoimmune disease patients, and 20 healthy subjects were screened to identify TAK-specific autoantibodies using human protein (HuProt) arrays. In Phase II, the identified candidate autoantibodies were validated with TAK-focused arrays using an additional cohort comprised of 109 TAK patients, 110 autoimmune disease patients, and 96 healthy subjects. Subsequently, the TAK-specific autoantibodies validated in Phase II were further confirmed using Western blot analysis. We identified and validated eight autoantibodies as potential TAK-specific diagnostic biomarkers, including anti-SPATA7, -QDPR, -SLC25A2, -PRH2, -DIXDC1, -IL17RB, -ZFAND4, and -NOLC1 antibodies, with AUC of 0.803, 0.801, 0.780, 0.696, 0.695, 0.678, 0.635 and 0.613, respectively. SPATA7 could distinguish TAK from healthy and disease controls with 73.4% sensitivity at 85.4% specificity, while QDPR showed 71.6% sensitivity at 86.4% specificity. SLC25A22 showed the highest sensitivity of 80.7%, but at lower specificity of 67.0%. In addition, PRH2, IL17RB and NOLC1 showed good specificities of 88.3%, 85.9% and 86.9%, respectively, but at lower sensitivities (<50%). Finally, DIXDC1 and ZFAND4 showed moderate performance as compared with the other autoantibodies. Using a decision tree model, we could reach a specificity of 94.2% with AUC of 0.843, a significantly improved performance as compared to that by each individual biomarker. The performance of three autoantibodies, namely anti-SPATA7, -QDPR and -PRH2, were successfully confirmed with Western blot analysis. Using this two-phase strategy, we identified and validated eight novel autoantibodies as TAK–specific biomarker candidates, three of which could be readily adopted in a clinical setting.