Full Native timsTOF PASEF-Enabled Quantitative Proteomics with the i2MassChroQ Software Package.

Ion mobility mass spectrometry has become popular in proteomics lately, in particular because the Bruker timsTOF instruments have found significant adoption in proteomics facilities. The Bruker's implementation of the ion mobility dimension generates massive amounts of mass spectrometric data that require carefully designed software both to extract meaningful information and to perform processing tasks at reasonable speed. In a historical move, the Bruker company decided to harness the skills of the scientific software development community by releasing to the public the timsTOF data file format specification. As a proteomics facility that has been developing Free Open Source Software (FOSS) solutions since decades, we took advantage of this opportunity to implement the very first FOSS proteomics complete solution to natively read the timsTOF data, low-level process them, and explore them in an integrated quantitative proteomics software environment. We dubbed our software i2MassChroQ because it implements a (peptide)identification-(protein)inference-mass-chromatogram-quantification processing workflow. The software benchmarking results reported in this paper show that i2MassChroQ performed better than competing software on two critical characteristics: (1) feature extraction capability and (2) protein quantitative dynamic range. Altogether, i2MassChroQ yielded better quantified protein numbers, both in a technical replicate MS runs setting and in a differential protein abundance analysis setting.

Fast alignment of mass spectra in large proteomics datasets, capturing dissimilarities arising from multiple complex modifications of peptides.

BACKGROUND: In proteomics, the interpretation of mass spectra representing peptides carrying multiple complex modifications remains challenging, as it is difficult to strike a balance between reasonable execution time, a limited number of false positives, and a huge search space allowing any number of modifications without a priori. The scientific community needs new developments in this area to aid in the discovery of novel post-translational modifications that may play important roles in disease. RESULTS: To make progress on this issue, we implemented SpecGlobX (SpecGlob eXTended to eXperimental spectra), a standalone Java application that quickly determines the best spectral alignments of a (possibly very large) list of Peptide-to-Spectrum Matches (PSMs) provided by any open modification search method, or generated by the user. As input, SpecGlobX reads a file containing spectra in MGF or mzML format and a semicolon-delimited spreadsheet describing the PSMs. SpecGlobX returns the best alignment for each PSM as output, splitting the mass difference between the spectrum and the peptide into one or more shifts while considering the possibility of non-aligned masses (a phenomenon resulting from many situations including neutral losses). SpecGlobX is fast, able to align one million PSMs in about 1.5 min on a standard desktop. Firstly, we remind the foundations of the algorithm and detail how we adapted SpecGlob (the method we previously developed following the same aim, but limited to the interpretation of perfect simulated spectra) to the interpretation of imperfect experimental spectra. Then, we highlight the interest of SpecGlobX as a complementary tool downstream to three open modification search methods on a large simulated spectra dataset. Finally, we ran SpecGlobX on a proteome-wide dataset downloaded from PRIDE to demonstrate that SpecGlobX functions just as well on simulated and experimental spectra. We then carefully analyzed a limited set of interpretations. CONCLUSIONS: SpecGlobX is helpful as a decision support tool, providing keys to interpret peptides carrying complex modifications still poorly considered by current open modification search software. Better alignment of PSMs enhances confidence in the identification of spectra provided by open modification search methods and should improve the interpretation rate of spectra.

Benefits of Iterative Searches of Large Databases to Interpret Large Human Gut Metaproteomic Data Sets.

The gut microbiota are increasingly considered as a main partner of human health. Metaproteomics enables us to move from the functional potential revealed by metagenomics to the functions actually operating in the microbiome. However, metaproteome deciphering remains challenging. In particular, confident interpretation of a myriad of MS/MS spectra can only be pursued with smart database searches. Here, we compare the interpretation of MS/MS data sets from 48 individual human gut microbiomes using three interrogation strategies of the dedicated Integrated nonredundant Gene Catalog (IGC 9.9 million genes from 1267 individual fecal samples) together with the Homo sapiens database: the classical single-step interrogation strategy and two iterative strategies (in either two or three steps) aimed at preselecting a reduced-sized, more targeted search space for the final peptide spectrum matching. Both iterative searches outperformed the single-step classical search in terms of the number of peptides and protein clusters identified and the depth of taxonomic and functional knowledge, and this was the most convincing with the three-step approach. However, iterative searches do not help in reducing variability of repeated analyses, which is inherent to the traditional data-dependent acquisition mode, but this variability did not affect the hierarchical relationship between replicates and all other samples.

X!TandemPipeline: A Tool to Manage Sequence Redundancy for Protein Inference and Phosphosite Identification.

X!TandemPipeline is a software designed to perform protein inference and to manage redundancy in the results of phosphosite identification by database search. It provides the minimal list of proteins or phosphosites that are present in a set of samples using grouping algorithms based on the principle of parsimony. Regarding proteins, a two-level classification is performed, where groups gather proteins sharing at least one peptide and subgroups gather proteins that are not distinguishable according to the identified peptides. Regarding phosphosites, an innovative approach based on the concept of phosphoisland is used to gather overlapping phosphopeptides. The graphical interface of X!TandemPipeline allows the users to launch X!tandem identification, to inspect spectra and to manually validate their assignment to peptides, to launch the grouping program, and to visualize elementary data as well as grouping and redundancy information. Identification results obtained from other search engines can also be processed. X!TandemPipeline results can be exported as ready-to-use tabulated files or as XML files that can be directly used by the PROTICdb database or by the MassChroQ quantification software. X!TandemPipeline runs fast, is easy to use, and can process hundreds of samples simultaneously. It is freely available under the GNU General Public License v3.0 at http://pappso.inra.fr/bioinfo/xtandempipeline/ .

Interactions between genotype and environment drive the metabolic phenotype within Escherichia coli isolates.

To gain insights into the adaptation of the Escherichia coli species to different environments, we monitored protein abundances using quantitative proteomics and measurements of enzymatic activities of central metabolism in a set of five representative strains grown in four contrasted culture media including human urine. Two hundred and thirty seven proteins representative of the genome-scale metabolic network were identified and classified into pathway categories. We found that nutrient resources shape the general orientation of metabolism through coordinated changes in the average abundances of proteins and in enzymatic activities that all belong to the same pathway category. For example, each culture medium induces a specific oxidative response whatever the strain. On the contrary, differences between strains concern isolated proteins and enzymes within pathway categories in single environments. Our study confirms the predominance of genotype by environment interactions at the proteomic and enzyme activity levels. The buffering of genetic variation when considering life-history traits suggests a multiplicity of evolutionary strategies. For instance, the uropathogenic isolate CFT073 shows a deregulation of iron demand and increased oxidative stress response.

O-Glycosylation of the N-terminal region of the serine-rich adhesin Srr1 of Streptococcus agalactiae explored by mass spectrometry.

Serine-rich (Srr) proteins exposed at the surface of Gram-positive bacteria are a family of adhesins that contribute to the virulence of pathogenic staphylococci and streptococci. Lectin-binding experiments have previously shown that Srr proteins are heavily glycosylated. We report here the first mass-spectrometry analysis of the glycosylation of Streptococcus agalactiae Srr1. After Srr1 enrichment and trypsin digestion, potential glycopeptides were identified in collision induced dissociation spectra using X! Tandem. The approach was then refined using higher energy collisional dissociation fragmentation which led to the simultaneous loss of sugar residues, production of diagnostic oxonium ions and backbone fragmentation for glycopeptides. This feature was exploited in a new open source software tool (SpectrumFinder) developed for this work. By combining these approaches, 27 glycopeptides corresponding to six different segments of the N-terminal region of Srr1 [93-639] were identified. Our data unambiguously indicate that the same protein residue can be modified with different glycan combinations including N-acetylhexosamine, hexose, and a novel modification that was identified as O-acetylated-N-acetylhexosamine. Lectin binding and monosaccharide composition analysis strongly suggested that HexNAc and Hex correspond to N-acetylglucosamine and glucose, respectively. The same protein segment can be modified with a variety of glycans generating a wide structural diversity of Srr1. Electron transfer dissociation was used to assign glycosylation sites leading to the unambiguous identification of six serines and one threonine residues. Analysis of purified Srr1 produced in mutant strains lacking accessory glycosyltransferase encoding genes demonstrates that O-GlcNAcylation is an initial step in Srr1 glycosylation that is likely required for subsequent decoration with Hex. In summary, our data obtained by a combination of fragmentation mass spectrometry techniques associated to a new software tool, demonstrate glycosylation heterogeneity of Srr1, characterize a new protein modification, and identify six glycosylation sites located in the N-terminal region of the protein.

Linking post-translational modifications and variation of phenotypic traits.

Enzymes can be post-translationally modified, leading to isoforms with different properties. The phenotypic consequences of the quantitative variability of isoforms have never been studied. We used quantitative proteomics to dissect the relationships between the abundances of the enzymes and isoforms of alcoholic fermentation, metabolic traits, and growth-related traits in Saccharomyces cerevisiae. Although the enzymatic pool allocated to the fermentation proteome was constant over the culture media and the strains considered, there was variation in abundance of individual enzymes and sometimes much more of their isoforms, which suggests the existence of selective constraints on total protein abundance and trade-offs between isoforms. Variations in abundance of some isoforms were significantly associated to metabolic traits and growth-related traits. In particular, cell size and maximum population size were highly correlated to the degree of N-terminal acetylation of the alcohol dehydrogenase. The fermentation proteome was found to be shaped by human selection, through the differential targeting of a few isoforms for each food-processing origin of strains. These results highlight the importance of post-translational modifications in the diversity of metabolic and life-history traits.

(p)ppGpp modifies RNAP function to confer ß-lactam resistance in a peptidoglycan-independent manner.

(p)ppGpp is a nucleotide alarmone that controls bacterial response to nutrient deprivation. Since elevated (p)ppGpp levels confer mecillinam resistance and are essential for broad-spectrum ß-lactam resistance as mediated by the ß-lactam-insensitive transpeptidase YcbB (LdtD), we hypothesized that (p)ppGpp might affect cell wall peptidoglycan metabolism. Here we report that (p)ppGpp-dependent ß-lactam resistance does not rely on any modification of peptidoglycan metabolism, as established by analysis of Escherichia coli peptidoglycan structure using high-resolution mass spectrometry. Amino acid substitutions in the ß or ß' RNA polymerase (RNAP) subunits, alone or in combination with the CRISPR interference-mediated downregulation of three of seven ribosomal RNA operons, were sufficient for resistance, although ß-lactams have no known impact on the RNAP or ribosomes. This implies that modifications of RNAP and ribosome functions are critical to prevent downstream effects of the inactivation of peptidoglycan transpeptidases by ß-lactams.

Management and dissemination of MS proteomic data with PROTICdb: example of a quantitative comparison between methods of protein extraction.

High throughput MS-based proteomic experiments generate large volumes of complex data and necessitate bioinformatics tools to facilitate their handling. Needs include means to archive data, to disseminate them to the scientific communities, and to organize and annotate them to facilitate their interpretation. We present here an evolution of PROTICdb, a database software that now handles MS data, including quantification. PROTICdb has been developed to be as independent as possible from tools used to produce the data. Biological samples and proteomics data are described using ontology terms. A Taverna workflow is embedded, thus permitting to automatically retrieve information related to identified proteins by querying external databases. Stored data can be displayed graphically and a "Query Builder" allows users to make sophisticated queries without knowledge on the underlying database structure. All resources can be accessed programmatically using a Java client API or RESTful web services, allowing the integration of PROTICdb in any portal. An example of application is presented, where proteins extracted from a maize leaf sample by four different methods were compared using a label-free shotgun method. Data are available at http://moulon.inra.fr/protic/public. PROTICdb thus provides means for data storage, enrichment, and dissemination of proteomics data.

Modern Metaproteomics: A Unique Tool to Characterize the Active Microbiome in Health and Diseases, and Pave the Road towards New Biomarkers-Example of Crohn's Disease and Ulcerative Colitis Flare-Ups.

Thanks to the latest developments in mass spectrometry, software and standards, metaproteomics is emerging as the vital complement of metagenomics, to make headway in understanding the actual functioning of living and active microbial communities. Modern metaproteomics offers new possibilities in the area of clinical diagnosis. This is illustrated here, for the still highly challenging diagnosis of intestinal bowel diseases (IBDs). Using bottom-up proteomics, we analyzed the gut metaproteomes of the same twenty faecal specimens processed either fresh or after a two-month freezing period. We focused on metaproteomes of microbial cell envelopes since it is an outstanding way of capturing host and host-microbe interaction signals. The protein profiles of pairs of fresh and frozen-thawed samples were closely related, making feasible deferred analysis in a distant diagnosis centre. The taxonomic and functional landscape of microbes in diverse IBD phenotypes-active ulcerative colitis, or active Crohn's disease either with ileo-colonic or exclusive colonic localization-differed from each other and from the controls. Based on their specific peptides, we could identify proteins that were either strictly overrepresented or underrepresented in all samples of one clinical group compared to all samples of another group, paving the road for promising additional diagnostic tool for IBDs.

A combined test for feature selection on sparse metaproteomics data-an alternative to missing value imputation.

One of the difficulties encountered in the statistical analysis of metaproteomics data is the high proportion of missing values, which are usually treated by imputation. Nevertheless, imputation methods are based on restrictive assumptions regarding missingness mechanisms, namely "at random" or "not at random". To circumvent these limitations in the context of feature selection in a multi-class comparison, we propose a univariate selection method that combines a test of association between missingness and classes, and a test for difference of observed intensities between classes. This approach implicitly handles both missingness mechanisms. We performed a quantitative and qualitative comparison of our procedure with imputation-based feature selection methods on two experimental data sets, as well as simulated data with various scenarios regarding the missingness mechanisms and the nature of the difference of expression (differential intensity or differential presence). Whereas we observed similar performances in terms of prediction on the experimental data set, the feature ranking and selection from various imputation-based methods were strongly divergent. We showed that the combined test reaches a compromise by correlating reasonably with other methods, and remains efficient in all simulated scenarios unlike imputation-based feature selection methods.

MassChroQ: a versatile tool for mass spectrometry quantification.

Recently, many software tools have been developed to perform quantification in LC-MS analyses. However, most of them are specific to either a quantification strategy (e.g. label-free or isotopic labelling) or a mass-spectrometry system (e.g. high or low resolution). In this context, we have developed MassChroQ (Mass Chromatogram Quantification), a versatile software that performs LC-MS data alignment and peptide quantification by peak area integration on extracted ion chromatograms. MassChroQ is suitable for quantification with or without labelling and is not limited to high-resolution systems. Peptides of interest (for example all the identified peptides) can be determined automatically, or manually by providing targeted m/z and retention time values. It can handle large experiments that include protein or peptide fractionation (as SDS-PAGE, 2-D LC). It is fully configurable. Every processing step is traceable, the produced data are in open standard formats and its modularity allows easy integration into proteomic pipelines. The output results are ready for use in statistical analyses. Evaluation of MassChroQ on complex label-free data obtained from low and high-resolution mass spectrometers showed low CVs for technical reproducibility (1.4%) and high coefficients of correlation to protein quantity (0.98). MassChroQ is freely available under the GNU General Public Licence v3.0 at http://pappso.inra.fr/bioinfo/masschroq/.

mineXpert2: Full-Depth Visualization and Exploration of MSn Mass Spectrometry Data.

mineXpert is a mass spectrometric data visualization and exploration software supporting only MS1 data that is aimed at proteomics scientists who do rarely require manual MS/MS data visualization and exploration (Rusconi, F. J. Proteome Res. 2019, 18, 2254-2259). In order to adapt it to new use cases in our facility and to widen its user base, mineXpert was entirely rewritten with the main aim of implementing MSn data support. Other feature additions were new data visualization and exploration methods, with an overhaul of the data plotting code to allow more flexible uses of mass data integration results. Further, the whole mass spectral data set can now be explored in a table view where the user may filter the data using a number of criteria that can be logically combined to pinpoint the smallest feature of interest. Ion mobility mass spectrometry is supported with specific data exploration and plotting. With mineXpert2, we provide a software program that will be of use to all mass spectrometrists, without restrictions on the field of endeavor, from pure chemistry to proteomics and metabolomics. As staff members of a mass spectrometry facility, we want to provide all users with a mass spectrometry data visualization and exploration software solution that frees them from the need to use closed-source vendor software. After conversion of the mass data to mzML, mineXpert2 requires no proprietary software whatsoever. The reference implementation is version 7.0.0 or greater. The software, a detailed user manual, and video tutorials are available at http://www.msxpertsuite.org.

Critical Assessment of MetaProteome Investigation (CAMPI): a multi-laboratory comparison of established workflows.

Metaproteomics has matured into a powerful tool to assess functional interactions in microbial communities. While many metaproteomic workflows are available, the impact of method choice on results remains unclear. Here, we carry out a community-driven, multi-laboratory comparison in metaproteomics: the critical assessment of metaproteome investigation study (CAMPI). Based on well-established workflows, we evaluate the effect of sample preparation, mass spectrometry, and bioinformatic analysis using two samples: a simplified, laboratory-assembled human intestinal model and a human fecal sample. We observe that variability at the peptide level is predominantly due to sample processing workflows, with a smaller contribution of bioinformatic pipelines. These peptide-level differences largely disappear at the protein group level. While differences are observed for predicted community composition, similar functional profiles are obtained across workflows. CAMPI demonstrates the robustness of present-day metaproteomics research, serves as a template for multi-laboratory studies in metaproteomics, and provides publicly available data sets for benchmarking future developments.

Comparative proteomics of leaf, stem, and root tissues of synthetic Brassica napus.

Comparative proteomics was applied to three vegetative organs of Brassica napus, the leaf, stem, and root using 2-DE. Among the >1600 analyzed spots, 43% were found to be common to all three organs, suggesting the existence of a "basal" or ubiquitous proteome composed of housekeeping proteins. The green organs, leaf, and stem, were closely related (approximately 80% common spots) while the root displayed more organ-specific polypeptides (approximately 10%). Reference maps were established using MS, allowing the identification of 93, 385, and 266 proteins in leaf, stem, and root proteomes, respectively. Bioinformatic analyses were also performed; in silico functional categorization and cellular localization allow obtaining a precise picture of the cell molecular network within vegetative organs. These proteome maps can be explored using the PROTICdb software at the following address: http://bioinformatique.moulon.inra.fr/proticdb/web_view/.

A method based on bead flows for spot detection on 2-D gel images.

A method has been developed to detect spots on 2-D gel images. It is based on an analogy with beads flowing uphill on the surface of the gel image and on the analysis of their paths. This method does not need image smoothing, thereby allowing for conservation of the resolution of the original image. In addition, it allows for the detection of overlapping spots, even when there is no valley between their centres.

The PROTICdb database for 2-DE proteomics.

PROTICdb is a web-based database mainly designed to store and analyze plant proteome data obtained by 2D polyacrylamide gel electrophoresis (2D PAGE) and mass spectrometry (MS). The goals of PROTICdb are (1) to store, track, and query information related to proteomic experiments, i.e., from tissue sampling to protein identification and quantitative measurements; and (2) to integrate information from the user's own expertise and other sources into a knowledge base, used to support data interpretation (e.g., for the determination of allelic variants or products of posttranslational modifications). Data insertion into the relational database of PROTICdb is achieved either by uploading outputs from Mélanie, PDQuest, IM2d, ImageMaster(tm) 2D Platinum v5.0, Progenesis, Sequest, MS-Fit, and Mascot software, or by filling in web forms (experimental design and methods). 2D PAGE-annotated maps can be displayed, queried, and compared through the GelBrowser. Quantitative data can be easily exported in a tabulated format for statistical analyses with any third-party software. PROTICdb is based on the Oracle or the PostgreSQLDataBase Management System (DBMS) and is freely available upon request at http://cms.moulon.inra.fr/content/view/14/44/.

The Arabidopsis TOR Kinase Specifically Regulates the Expression of Nuclear Genes Coding for Plastidic Ribosomal Proteins and the Phosphorylation of the Cytosolic Ribosomal Protein S6.

Protein translation is an energy consuming process that has to be fine-tuned at both the cell and organism levels to match the availability of resources. The target of rapamycin kinase (TOR) is a key regulator of a large range of biological processes in response to environmental cues. In this study, we have investigated the effects of TOR inactivation on the expression and regulation of Arabidopsis ribosomal proteins at different levels of analysis, namely from transcriptomic to phosphoproteomic. TOR inactivation resulted in a coordinated down-regulation of the transcription and translation of nuclear-encoded mRNAs coding for plastidic ribosomal proteins, which could explain the chlorotic phenotype of the TOR silenced plants. We have identified in the 5' untranslated regions (UTRs) of this set of genes a conserved sequence related to the 5' terminal oligopyrimidine motif, which is known to confer translational regulation by the TOR kinase in other eukaryotes. Furthermore, the phosphoproteomic analysis of the ribosomal fraction following TOR inactivation revealed a lower phosphorylation of the conserved Ser240 residue in the C-terminal region of the 40S ribosomal protein S6 (RPS6). These results were confirmed by Western blot analysis using an antibody that specifically recognizes phosphorylated Ser240 in RPS6. Finally, this antibody was used to follow TOR activity in plants. Our results thus uncover a multi-level regulation of plant ribosomal genes and proteins by the TOR kinase.

A two-dimensional proteome map of maize endosperm.

We have established a proteome reference map for maize (Zea mays L.) endosperm by means of two-dimensional gel electrophoresis and protein identification with LC-MS/MS analysis. This investigation focussed on proteins in major spots in a 4-7 pI range and 10-100 kDa M(r) range. Among the 632 protein spots processed, 496 were identified by matching against the NCBInr and ZMtuc-tus databases (using the SEQUEST software). Forty-two per cent of the proteins were identified against maize sequences, 23% against rice sequences and 21% against Arabidopsis sequences. Identified proteins were not only cytoplasmic but also nuclear, mitochondrial or amyloplastic. Metabolic processes, protein destination, protein synthesis, cell rescue, defense, cell death and ageing are the most abundant functional categories, comprising almost half of the 632 proteins analyzed in our study. This proteome map constitutes a powerful tool for physiological studies and is the first step for investigating the maize endosperm development.

PROTICdb: a web-based application to store, track, query, and compare plant proteome data.

PROTICdb is a web-based application, mainly designed to store and analyze plant proteome data obtained by two-dimensional polyacrylamide gel electrophoresis (2-D PAGE) and mass spectrometry (MS). The purposes of PROTICdb are (i) to store, track, and query information related to proteomic experiments, i.e., from tissue sampling to protein identification and quantitative measurements, and (ii) to integrate information from the user's own expertise and other sources into a knowledge base, used to support data interpretation (e.g., for the determination of allelic variants or products of post-translational modifications). Data insertion into the relational database of PROTICdb is achieved either by uploading outputs of image analysis and MS identification software, or by filling web forms. 2-D PAGE annotated maps can be displayed, queried, and compared through a graphical interface. Links to external databases are also available. Quantitative data can be easily exported in a tabulated format for statistical analyses. PROTICdb is based on the Oracle or the PostgreSQL Database Management System and is freely available upon request at the following URL: http://moulon.inra.fr/ bioinfo/PROTICdb.