Scalable Data Analysis in Proteomics and Metabolomics Using BioContainers and Workflows Engines

Innovative analytical methodologies for characterizing chemical exposure with a view to next-generation risk assessment

The chemical burden on the environment and human population is increasing. Consequently, regulatory risk assessment must keep pace to manage, reduce, and prevent adverse impacts on human and environmental health associated with hazardous chemicals. Surveillance of chemicals of known, emerging, or potential future concern, entering the environment-food-human continuum is needed to document the reality of risks posed by chemicals on ecosystem and human health from a one health perspective, feed into early warning systems and support public policies for exposure mitigation provisions and safe and sustainable by design strategies. The use of less-conventional sampling strategies and integration of full-scan, high-resolution mass spectrometry and effect-directed analysis in environmental and human monitoring programmes have the potential to enhance the screening and identification of a wider range of chemicals of known, emerging or potential future concern. Here, we outline the key needs and recommendations identified within the European Partnership for Assessment of Risks from Chemicals (PARC) project for leveraging these innovative methodologies to support the development of next-generation chemical risk assessment.

TopDownApp: An open and modular platform for analysis and visualisation of top-down proteomics data

Although Top-down (TD) proteomics techniques, aimed at the analysis of intact proteins and proteoforms, are becoming increasingly popular, efforts are needed at different levels to generalise their adoption. In this context, there are numerous improvements that are possible in the area of open science practices, including a greater application of the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles. These include, for example, increased data sharing practices and readily available open data standards. Additionally, the field would benefit from the development of open data analysis workflows that can enable data reuse of public datasets, something that is increasingly common in other proteomics fields.

Tissue-based absolute quantification using large-scale TMT and LFQ experiments

Relative and absolute intensity-based protein quantification across cell lines, tissue atlases and tumour datasets is increasingly available in public datasets. These atlases enable researchers to explore fundamental biological questions, such as protein existence, expression location, quantity and correlation with RNA expression. Most studies provide MS1 feature-based label-free quantitative (LFQ) datasets; however, growing numbers of isobaric tandem mass tags (TMT) datasets remain unexplored. Here, we compare traditional intensity-based absolute quantification (iBAQ) proteome abundance ranking to an analogous method using reporter ion proteome abundance ranking with data from an experiment where LFQ and TMT were measured on the same samples. This new TMT method substitutes reporter ion intensities for MS1 feature intensities in the iBAQ framework. Additionally, we compared LFQ-iBAQ values to TMT-iBAQ values from two independent large-scale tissue atlas datasets (one LFQ and one TMT) using robust bottom-up proteomic identification, normalisation and quantitation workflows.

UmetaFlow: an untargeted metabolomics workflow for high-throughput data processing and analysis

Metabolomics experiments generate highly complex datasets, which are time and work-intensive, sometimes even error-prone if inspected manually. Therefore, new methods for automated, fast, reproducible, and accurate data processing and dereplication are required. Here, we present UmetaFlow, a computational workflow for untargeted metabolomics that combines algorithms for data pre-processing, spectral matching, molecular formula and structural predictions, and an integration to the GNPS workflows Feature-Based Molecular Networking and Ion Identity Molecular Networking for downstream analysis. UmetaFlow is implemented as a Snakemake workflow, making it easy to use, scalable, and reproducible. For more interactive computing, visualization, as well as development, the workflow is also implemented in Jupyter notebooks using the Python programming language and a set of Python bindings to the OpenMS algorithms (pyOpenMS). Finally, UmetaFlow is also offered as a web-based Graphical User Interface for parameter optimization and processing of smaller-sized datasets. UmetaFlow was validated with in-house LC-MS/MS datasets of actinomycetes producing known secondary metabolites, as well as commercial standards, and it detected all expected features and accurately annotated 76% of the molecular formulas and 65% of the structures. As a more generic validation, the publicly available MTBLS733 and MTBLS736 datasets were used for benchmarking, and UmetaFlow detected more than 90% of all ground truth features and performed exceptionally well in quantification and discriminating marker selection. We anticipate that UmetaFlow will provide a useful platform for the interpretation of large metabolomics datasets.

Deoxyhypusine hydroxylase: A novel therapeutic target differentially expressed in short-term vs long-term survivors of glioblastoma

Glioblastoma (GB) is the most aggressive neoplasm of the brain. Poor prognosis is mainly attributed to tumor heterogeneity, invasiveness and drug resistance. Only a small fraction of GB patients survives longer than 24 months from the time of diagnosis (ie, long-term survivors [LTS]). In our study, we aimed to identify molecular markers associated with favorable GB prognosis as a basis to develop therapeutic applications to improve patients' outcome. We have recently assembled a proteogenomic dataset of 87 GB clinical samples of varying survival rates. Following RNA-seq and mass spectrometry (MS)-based proteomics analysis, we identified several differentially expressed genes and proteins, including some known cancer-related pathways and some less established that showed higher expression in short-term (<6 months) survivors (STS) compared to LTS. One such target found was deoxyhypusine hydroxylase (DOHH), which is known to be involved in the biosynthesis of hypusine, an unusual amino acid essential for the function of the eukaryotic translation initiation factor 5A (eIF5A), which promotes tumor growth. We consequently validated DOHH overexpression in STS samples by quantitative polymerase chain reaction (qPCR) and immunohistochemistry. We further showed robust inhibition of proliferation, migration and invasion of GB cells following silencing of DOHH with short hairpin RNA (shRNA) or inhibition of its activity with small molecules, ciclopirox and deferiprone. Moreover, DOHH silencing led to significant inhibition of tumor progression and prolonged survival in GB mouse models. Searching for a potential mechanism by which DOHH promotes tumor aggressiveness, we found that it supports the transition of GB cells to a more invasive phenotype via epithelial-mesenchymal transition (EMT)-related pathways.

A peptide-centric approach to analyse quantitative proteomics data- an application to prostate cancer biomarker discovery

Bottom-up proteomics is a popular approach in molecular biomarker research. However, protein analysts have realized the limitations of protein-based approaches for identifying and quantifying proteins in complex samples, such as the identification of peptides sequences shared by multiple proteins and the difficulty in identifying modified peptides. Thus, there are many exciting opportunities to improve analysis methods. Here, an alternative method focused on peptide analysis is proposed as a complement to the conventional proteomics data analysis. To investigate this hypothesis, a peptide-centric approach was applied to reanalyse a urine proteome dataset of samples from prostate cancer patients and controls. The results were compared with the conventional protein-centric approach. The relevant proteins/peptides to discriminate the groups were detected based on two approaches, p-value and VIP values obtained by a PLS-DA model. A comparison of the two strategies revealed high inconsistency between protein and peptide information and greater involvement of peptides in key PCa processes. This peptide analysis unveiled discriminative features that are lost when proteins are analyzed as homogeneous entities. This type of analysis is innovative in PCa and integrated with the widely used protein-centric approach might provide a more comprehensive view of this disease and revolutionize biomarker discovery. SIGNIFICANCE: In this study, the application of a protein and peptide-centric approaches to reanalyse a urine proteome dataset from prostate cancer (PCa) patients and controls showed that many relevant proteins/peptides are missed by the conservative nature of p-value in statistical tests, therefore, the inclusion of variable selection methods in the analysis of the dataset reported in this work is fruitful. Comparison of protein- and peptide-based approaches revealed a high inconsistency between protein and peptide information and a greater involvement of peptides in key PCa processes. These results provide a new perspective to analyse proteomics data and detect relevant targets based on the integration of peptide and protein information. This data integration allows to unravel discriminative features that normally go unnoticed, to have a more comprehensive view of the disease pathophysiology and to open new avenues for the discovery of biomarkers.

psm_utils: A High-Level Python API for Parsing and Handling Peptide-Spectrum Matches and Proteomics Search Results

A plethora of proteomics search engine output file formats are in circulation. This lack of standardized output files greatly complicates generic downstream processing of peptide-spectrum matches (PSMs) and PSM files. While standards exist to solve this problem, these are far from universally supported by search engines. Moreover, software libraries are available to read a selection of PSM file formats, but a package to parse PSM files into a unified data structure has been missing. Here, we present psm_utils, a Python package to read and write various PSM file formats and to handle peptidoforms, PSMs, and PSM lists in a unified and user-friendly Python-, command line-, and web-interface. psm_utils was developed with pragmatism and maintainability in mind, adhering to community standards and relying on existing packages where possible. The Python API and command line interface greatly facilitate handling various PSM file formats. Moreover, a user-friendly web application was built using psm_utils that allows anyone to interconvert PSM files and retrieve basic PSM statistics. psm_utils is freely available under the permissive Apache2 license at https://github.com/compomics/psm_utils.

The META tool optimizes metagenomic analyses across sequencing platforms and classifiers

A major challenge in the field of metagenomics is the selection of the correct combination of sequencing platform and downstream metagenomic analysis algorithm, or "classifier". Here, we present the Metagenomic Evaluation Tool Analyzer (META), which produces simulated data and facilitates platform and algorithm selection for any given metagenomic use case. META-generated in silico read data are modular, scalable, and reflect user-defined community profiles, while the downstream analysis is done using a variety of metagenomic classifiers. Reported results include information on resource utilization, time-to-answer, and performance. Real-world data can also be analyzed using selected classifiers and results benchmarked against simulations. To test the utility of the META software, simulated data was compared to real-world viral and bacterial metagenomic samples run on four different sequencers and analyzed using 12 metagenomic classifiers. Lastly, we introduce "META Score": a unified, quantitative value which rates an analytic classifier's ability to both identify and count taxa in a representative sample.

Proteomic repository data submission, dissemination, and reuse: key messages

INTRODUCTION

The creation of ProteomeXchange data workflows in 2012 transformed the field of proteomics, consisting of the standardization of data submission and dissemination and enabling the widespread reanalysis of public MS proteomics data worldwide. ProteomeXchange has triggered a growing trend toward public dissemination of proteomics data, facilitating the assessment, reuse, comparative analyses, and extraction of new findings from public datasets. By 2022, the consortium is integrated by PRIDE, PeptideAtlas, MassIVE, jPOST, iProX, and Panorama Public.

AREAS COVERED

Here, we review and discuss the current ecosystem of resources, guidelines, and file formats for proteomics data dissemination and reanalysis. Special attention is drawn to new exciting quantitative and post-translational modification-oriented resources. The challenges and future directions on data depositions including the lack of metadata and cloud-based and high-performance software solutions for fast and reproducible reanalysis of the available data are discussed.

EXPERT OPINION

The success of ProteomeXchange and the amount of proteomics data available in the public domain have triggered the creation and/or growth of other protein knowledgebase resources. Data reuse is a leading, active, and evolving field; supporting the creation of new formats, tools, and workflows to rediscover and reshape the public proteomics data.

MaxQuant and MSstats in Galaxy Enable Reproducible Cloud-Based Analysis of Quantitative Proteomics Experiments for Everyone

Quantitative mass spectrometry-based proteomics has become a high-throughput technology for the identification and quantification of thousands of proteins in complex biological samples. Two frequently used tools, MaxQuant and MSstats, allow for the analysis of raw data and finding proteins with differential abundance between conditions of interest. To enable accessible and reproducible quantitative proteomics analyses in a cloud environment, we have integrated MaxQuant (including TMTpro 16/18plex), Proteomics Quality Control (PTXQC), MSstats, and MSstatsTMT into the open-source Galaxy framework. This enables the web-based analysis of label-free and isobaric labeling proteomics experiments via Galaxy's graphical user interface on public clouds. MaxQuant and MSstats in Galaxy can be applied in conjunction with thousands of existing Galaxy tools and integrated into standardized, sharable workflows. Galaxy tracks all metadata and intermediate results in analysis histories, which can be shared privately for collaborations or publicly, allowing full reproducibility and transparency of published analysis. To further increase accessibility, we provide detailed hands-on training materials. The integration of MaxQuant and MSstats into the Galaxy framework enables their usage in a reproducible way on accessible large computational infrastructures, hence realizing the foundation for high-throughput proteomics data science for everyone.

The PRIDE database resources in 2022: a hub for mass spectrometry-based proteomics evidences

The PRoteomics IDEntifications (PRIDE) database (https://www.ebi.ac.uk/pride/) is the world's largest data repository of mass spectrometry-based proteomics data. PRIDE is one of the founding members of the global ProteomeXchange (PX) consortium and an ELIXIR core data resource. In this manuscript, we summarize the developments in PRIDE resources and related tools since the previous update manuscript was published in Nucleic Acids Research in 2019. The number of submitted datasets to PRIDE Archive (the archival component of PRIDE) has reached on average around 500 datasets per month during 2021. In addition to continuous improvements in PRIDE Archive data pipelines and infrastructure, the PRIDE Spectra Archive has been developed to provide direct access to the submitted mass spectra using Universal Spectrum Identifiers. As a key point, the file format MAGE-TAB for proteomics has been developed to enable the improvement of sample metadata annotation. Additionally, the resource PRIDE Peptidome provides access to aggregated peptide/protein evidences across PRIDE Archive. Furthermore, we will describe how PRIDE has increased its efforts to reuse and disseminate high-quality proteomics data into other added-value resources such as UniProt, Ensembl and Expression Atlas.

ProMetIS, deep phenotyping of mouse models by combined proteomics and metabolomics analysis

Genes are pleiotropic and getting a better knowledge of their function requires a comprehensive characterization of their mutants. Here, we generated multi-level data combining phenomic, proteomic and metabolomic acquisitions from plasma and liver tissues of two C57BL/6 N mouse models lacking the Lat (linker for activation of T cells) and the Mx2 (MX dynamin-like GTPase 2) genes, respectively. Our dataset consists of 9 assays (1 preclinical, 2 proteomics and 6 metabolomics) generated with a fully non-targeted and standardized approach. The data and processing code are publicly available in the ProMetIS R package to ensure accessibility, interoperability, and reusability. The dataset thus provides unique molecular information about the physiological role of the Lat and Mx2 genes. Furthermore, the protocols described herein can be easily extended to a larger number of individuals and tissues. Finally, this resource will be of great interest to develop new bioinformatic and biostatistic methods for multi-omics data integration.

Reproducible, scalable, and shareable analysis pipelines with bioinformatics workflow managers

The rapid growth of high-throughput technologies has transformed biomedical research. With the increasing amount and complexity of data, scalability and reproducibility have become essential not just for experiments, but also for computational analysis. However, transforming data into information involves running a large number of tools, optimizing parameters, and integrating dynamically changing reference data. Workflow managers were developed in response to such challenges. They simplify pipeline development, optimize resource usage, handle software installation and versions, and run on different compute platforms, enabling workflow portability and sharing. In this Perspective, we highlight key features of workflow managers, compare commonly used approaches for bioinformatics workflows, and provide a guide for computational and noncomputational users. We outline community-curated pipeline initiatives that enable novice and experienced users to perform complex, best-practice analyses without having to manually assemble workflows. In sum, we illustrate how workflow managers contribute to making computational analysis in biomedical research shareable, scalable, and reproducible.

Twenty years of proteomics in radiation biology - a look back

PURPOSE

The aim of this article is to describe the technical development in proteomics during the last two decades with the focus on its use in radiation biology. It is written from a subjective point of view and aims not to be a scientific review of the subject.

CONCLUSION

Proteomics is a fast developing technique and it has already contributed greatly to our understanding of biological mechanisms following radiation exposure. Novel proteomics approaches can be used in adequately designed cellular and animal experiments and above all in big clinical trials to investigate effects of ionizing radiation in the future.

Genetic determinants of seed protein plasticity in response to the environment in Medicago truncatula

As the frequency of extreme environmental events is expected to increase with climate change, identifying candidate genes for stabilizing the protein composition of legume seeds or optimizing this in a given environment is increasingly important. To elucidate the genetic determinants of seed protein plasticity, major seed proteins from 200 ecotypes of Medicago truncatula grown in four contrasting environments were quantified after one-dimensional electrophoresis. The plasticity index of these proteins was recorded for each genotype as the slope of Finlay and Wilkinson's regression and then used for genome-wide association studies (GWASs), enabling the identification of candidate genes for determining this plasticity. This list was enriched in genes related to transcription, DNA repair and signal transduction, with many of them being stress responsive. Other over-represented genes were related to sulfur and aspartate family pathways leading to the synthesis of the nutritionally essential amino acids methionine and lysine. By placing these genes in metabolic pathways, and using a M. truncatula mutant impaired in regenerating methionine from S-methylmethionine, we discovered that methionine recycling pathways are major contributors to globulin composition establishment and plasticity. These data provide a unique resource of genes that can be targeted to mitigate negative impacts of environmental stresses on seed protein composition.

Big data, integrative omics and network biology

A cell integrates various signals through a network of biomolecules that crosstalk to synergistically regulate the replication, transcription, translation and other metabolic activities of a cell. These networks regulate signal perception and processing that drives biological functions. The biological complexity cannot be fully captured by a single -omics discipline. The holistic study of an organism-in health, perturbation, exposure to environment and disease, is studied under systems biology. The bottom-up molecular approaches (genes, mRNA, protein, metabolite, etc.) have laid the foundation of current biological knowledge covering the horizon from viruses, bacteria, fungi, plants and animals. Yet, these techniques provide a rather myopic view of biology at the molecular level. To understand how the interconnected molecular components are formed and rewired in disease or exposure to environmental stimuli is the holy grail of modern biology. The omics era was heralded by the genomics revolution but advanced sequencing techniques are now also ubiquitous in transcriptomics, proteomics, metabolomics and lipidomics. Multi-omics data analysis and integration techniques are driving the quest for deeper insights into how the different layers of biomolecules talk to each other in diverse contexts.

Cloudy with a Chance of Peptides: Accessibility, Scalability, and Reproducibility with Cloud-Hosted Environments

Cloud-hosted environments offer known benefits when computational needs outstrip affordable local workstations, enabling high-performance computation without a physical cluster. What has been less apparent, especially to novice users, is the transformative potential for cloud-hosted environments to bridge the digital divide that exists between poorly funded and well-resourced laboratories, and to empower modern research groups with remote personnel and trainees. Using cloud-based proteomic bioinformatic pipelines is not predicated on analyzing thousands of files, but instead can be used to improve accessibility during remote work, extreme weather, or working with under-resourced remote trainees. The general benefits of cloud-hosted environments also allow for scalability and encourage reproducibility. Since one possible hurdle to adoption is awareness, this paper is written with the nonexpert in mind. The benefits and possibilities of using a cloud-hosted environment are emphasized by describing how to setup an example workflow to analyze a previously published label-free data-dependent acquisition mass spectrometry data set of mammalian urine. Cost and time of analysis are compared using different computational tiers, and important practical considerations are described. Overall, cloud-hosted environments offer the potential to solve large computational problems, but more importantly can enable and accelerate research in smaller research groups with inadequate infrastructure and suboptimal local computational resources.

BioContainers Registry: Searching Bioinformatics and Proteomics Tools, Packages, and Containers

BioContainers is an open-source project that aims to create, store, and distribute bioinformatics software containers and packages. The BioContainers community has developed a set of guidelines to standardize software containers including the metadata, versions, licenses, and software dependencies. BioContainers supports multiple packaging and container technologies such as Conda, Docker, and Singularity. The BioContainers provide over 9000 bioinformatics tools, including more than 200 proteomics and mass spectrometry tools. Here we introduce the BioContainers Registry and Restful API to make containerized bioinformatics tools more findable, accessible, interoperable, and reusable (FAIR). The BioContainers Registry provides a fast and convenient way to find and retrieve bioinformatics tool packages and containers. By doing so, it will increase the use of bioinformatics packages and containers while promoting replicability and reproducibility in research.

Simple and Efficient Data Analysis Dissemination for Individual Laboratories

Scientific progress comes as we build upon the work of others. Implicit in this advance is that we have access to and can thoroughly examine the work of others. It is important to recognize that our scholarly work as scientists encompasses not only experimental design and data collection but also our analytical methods. Thus when communicating biology experiments, especially those that utilize molecular omics data, the analysis methods that connect raw data to scientific conclusions must be presented with sufficient clarity that others can reproduce our exact work. Although there are many resources for sharing raw data files, there is currently not a widely utilized method for sharing analysis methods. We present a semistructured pattern for sharing analysis methods that is simple and efficient and can be implemented by individual laboratories using existing software. This pattern requires three types of files in a publicly accessible repository, such as GitHub: (1) data files, (2) a universal I/O script that parses all data files, and (3) analysis scripts creating figures and metrics reported in the manuscript. We suggest additional conventions to improve the readability and provide a template repository for the pattern. Sharing our exact analysis methods as software, in addition to their narrative description in a manuscript, will ensure reproducibility and transparency. Importantly, the pattern we present does not require new infrastructure and can be achieved without advanced computing skills.