A novel clinical metaproteomics workflow enables bioinformatic analysis of host-microbe dynamics in disease

Clinical metaproteomics has the potential to offer insights into the host-microbiome interactions underlying diseases. However, the field faces challenges in characterizing microbial proteins found in clinical samples, which are usually present at low abundance relative to the host proteins. As a solution, we have developed an integrated workflow coupling mass spectrometry-based analysis with customized bioinformatic identification, quantification and prioritization of microbial and host proteins, enabling targeted assay development to investigate host-microbe dynamics in disease. The bioinformatics tools are implemented in the Galaxy ecosystem, offering the development and dissemination of complex bioinformatic workflows. The modular workflow integrates MetaNovo (to generate a reduced protein database), SearchGUI/PeptideShaker and MaxQuant (to generate peptide-spectral matches (PSMs) and quantification), PepQuery2 (to verify the quality of PSMs), and Unipept and MSstatsTMT (for taxonomy and functional annotation). We have utilized this workflow in diverse clinical samples, from the characterization of nasopharyngeal swab samples to bronchoalveolar lavage fluid. Here, we demonstrate its effectiveness via analysis of residual fluid from cervical swabs. The complete workflow, including training data and documentation, is available via the Galaxy Training Network, empowering non-expert researchers to utilize these powerful tools in their clinical studies.

Mass Spectrometry-Based Strategies for Assessing Human Exposure Using Hemoglobin Adductomics

Hemoglobin (Hb) adducts are widely used in human biomonitoring due to the high abundance of hemoglobin in human blood, its reactivity toward electrophiles, and adducted protein stability for up to 120 days. In the present paper, we compared three methods of analysis of hemoglobin adducts: mass spectrometry of derivatized N-terminal Val adducts, mass spectrometry of N-terminal adducted hemoglobin peptides, and limited proteolysis mass spectrometry . Blood from human donors was incubated with a selection of contact allergens and other electrophiles, after which hemoglobin was isolated and subjected to three analysis methods. We found that the FIRE method was able to detect and reliably quantify N-terminal adducts of acrylamide, acrylic acid, glycidic acid, and 2,3-epoxypropyl phenyl ether (PGE), but it was less efficient for 2-methyleneglutaronitrile (2-MGN) and failed to detect 1-chloro-2,4-dinitrobenzene (DNCB). By contrast, bottom-up proteomics was able to determine the presence of adducts from all six electrophiles at both the N-terminus and reactive hemoglobin side chains. Limited proteolysis mass spectrometry, studied for four contact allergens (three electrophiles and a metal salt), was able to determine the presence of covalent hemoglobin adducts with one of the three electrophiles (DNCB) and coordination complexation with the nickel salt. Together, these approaches represent complementary tools in the study of the hemoglobin adductome.

Discovery of Modified Metabolites, Secondary Metabolites, and Xenobiotics by Structure-Oriented LC-MS/MS

Exogenous compounds and metabolites derived from therapeutics, microbiota, or environmental exposures directly interact with endogenous metabolic pathways, influencing disease pathogenesis and modulating outcomes of clinical interventions. With few spectral library references, the identification of covalently modified biomolecules, secondary metabolites, and xenobiotics is a challenging task using global metabolomics profiling approaches. Numerous liquid chromatography-coupled mass spectrometry (LC-MS) small molecule analytical workflows have been developed to curate global profiling experiments for specific compound groups of interest. These workflows exploit shared structural moiety, functional groups, or elemental composition to discover novel and undescribed compounds through nontargeted small molecule discovery pipelines. This Review introduces the concept of structure-oriented LC-MS discovery methodology and aims to highlight common approaches employed for the detection and characterization of covalently modified biomolecules, secondary metabolites, and xenobiotics. These approaches represent a combination of instrument-dependent and computational techniques to rapidly curate global profiling experiments to detect putative ions of interest based on fragmentation patterns, predictable phase I or phase II metabolic transformations, or rare elemental composition. Application of these methods is explored for the detection and identification of novel and undescribed biomolecules relevant to the fields of toxicology, pharmacology, and drug discovery. Continued advances in these methods expand the capacity for selective compound discovery and characterization that promise remarkable insights into the molecular interactions of exogenous chemicals with host biochemical pathways.

Metaproteomic Analysis of Nasopharyngeal Swab Samples to Identify Microbial Peptides in COVID-19 Patients

During the COVID-19 pandemic, impaired immunity and medical interventions resulted in cases of secondary infections. The clinical difficulties and dangers associated with secondary infections in patients necessitate the exploration of their microbiome. Metaproteomics is a powerful approach to study the taxonomic composition and functional status of the microbiome under study. In this study, the mass spectrometry (MS)-based data of nasopharyngeal swab samples from COVID-19 patients was used to investigate the metaproteome. We have established a robust bioinformatics workflow within the Galaxy platform, which includes (a) generation of a tailored database of the common respiratory tract pathogens, (b) database search using multiple search algorithms, and (c) verification of the detected microbial peptides. The microbial peptides detected in this study, belong to several opportunistic pathogens such as Streptococcus pneumoniae, Klebsiella pneumoniae, Rhizopus microsporus, and Syncephalastrum racemosum. Microbial proteins with a role in stress response, gene expression, and DNA repair were found to be upregulated in severe patients compared to negative patients. Using parallel reaction monitoring (PRM), we confirmed some of the microbial peptides in fresh clinical samples. MS-based clinical metaproteomics can serve as a powerful tool for detection and characterization of potential pathogens, which can significantly impact the diagnosis and treatment of patients.

Integrative meta-omics in Galaxy and beyond

BACKGROUND

'Omics methods have empowered scientists to tackle the complexity of microbial communities on a scale not attainable before. Individually, omics analyses can provide great insight; while combined as "meta-omics", they enhance the understanding of which organisms occupy specific metabolic niches, how they interact, and how they utilize environmental nutrients. Here we present three integrative meta-omics workflows, developed in Galaxy, for enhanced analysis and integration of metagenomics, metatranscriptomics, and metaproteomics, combined with our newly developed web-application, ViMO (Visualizer for Meta-Omics) to analyse metabolisms in complex microbial communities.

RESULTS

In this study, we applied the workflows on a highly efficient cellulose-degrading minimal consortium enriched from a biogas reactor to analyse the key roles of uncultured microorganisms in complex biomass degradation processes. Metagenomic analysis recovered metagenome-assembled genomes (MAGs) for several constituent populations including Hungateiclostridium thermocellum, Thermoclostridium stercorarium and multiple heterogenic strains affiliated to Coprothermobacter proteolyticus. The metagenomics workflow was developed as two modules, one standard, and one optimized for improving the MAG quality in complex samples by implementing a combination of single- and co-assembly, and dereplication after binning. The exploration of the active pathways within the recovered MAGs can be visualized in ViMO, which also provides an overview of the MAG taxonomy and quality (contamination and completeness), and information about carbohydrate-active enzymes (CAZymes), as well as KEGG annotations and pathways, with counts and abundances at both mRNA and protein level. To achieve this, the metatranscriptomic reads and metaproteomic mass-spectrometry spectra are mapped onto predicted genes from the metagenome to analyse the functional potential of MAGs, as well as the actual expressed proteins and functions of the microbiome, all visualized in ViMO.

CONCLUSION

Our three workflows for integrative meta-omics in combination with ViMO presents a progression in the analysis of 'omics data, particularly within Galaxy, but also beyond. The optimized metagenomics workflow allows for detailed reconstruction of microbial community consisting of MAGs with high quality, and thus improves analyses of the metabolism of the microbiome, using the metatranscriptomics and metaproteomics workflows.

A Galaxy of informatics resources for MS-based proteomics

INTRODUCTION

Continuous advances in mass spectrometry (MS) technologies have enabled deeper and more reproducible proteome characterization and a better understanding of biological systems when integrated with other 'omics data. Bioinformatic resources meeting the analysis requirements of increasingly complex MS-based proteomic data and associated multi-omic data are critically needed. These requirements included availability of software that would span diverse types of analyses, scalability for large-scale, compute-intensive applications, and mechanisms to ease adoption of the software.

AREAS COVERED

The Galaxy ecosystem meets these requirements by offering a multitude of open-source tools for MS-based proteomics analyses and applications, all in an adaptable, scalable, and accessible computing environment. A thriving global community maintains these software and associated training resources to empower researcher-driven analyses.

EXPERT OPINION

The community-supported Galaxy ecosystem remains a crucial contributor to basic biological and clinical studies using MS-based proteomics. In addition to the current status of Galaxy-based resources, we describe ongoing developments for meeting emerging challenges in MS-based proteomic informatics. We hope this review will catalyze increased use of Galaxy by researchers employing MS-based proteomics and inspire software developers to join the community and implement new tools, workflows, and associated training content that will add further value to this already rich ecosystem.

An optimized workflow for MS-based quantitative proteomics of challenging clinical bronchoalveolar lavage fluid (BALF) samples

BACKGROUND

Clinical bronchoalveolar lavage fluid (BALF) samples are rich in biomolecules, including proteins, and useful for molecular studies of lung health and disease. However, mass spectrometry (MS)-based proteomic analysis of BALF is challenged by the dynamic range of protein abundance, and potential for interfering contaminants. A robust, MS-based proteomics compatible sample preparation workflow for BALF samples, including those of small and large volume, would be useful for many researchers.

RESULTS

We have developed a workflow that combines high abundance protein depletion, protein trapping, clean-up, and in-situ tryptic digestion, that is compatible with either qualitative or quantitative MS-based proteomic analysis. The workflow includes a value-added collection of endogenous peptides for peptidomic analysis of BALF samples, if desired, as well as amenability to offline semi-preparative or microscale fractionation of complex peptide mixtures prior to LC-MS/MS analysis, for increased depth of analysis. We demonstrate the effectiveness of this workflow on BALF samples collected from COPD patients, including for smaller sample volumes of 1-5 mL that are commonly available from the clinic. We also demonstrate the repeatability of the workflow as an indicator of its utility for quantitative proteomic studies.

CONCLUSIONS

Overall, our described workflow consistently provided high quality proteins and tryptic peptides for MS analysis. It should enable researchers to apply MS-based proteomics to a wide-variety of studies focused on BALF clinical specimens.

Lung proteome and metabolome endotype in HIV-associated obstructive lung disease

Purpose

Obstructive lung disease is increasingly common among persons with HIV, both smokers and nonsmokers. We used aptamer proteomics to identify proteins and associated pathways in HIV-associated obstructive lung disease.

Methods

Bronchoalveolar lavage fluid (BALF) samples from 26 persons living with HIV with obstructive lung disease were matched to persons living with HIV without obstructive lung disease based on age, smoking status and antiretroviral treatment. 6414 proteins were measured using SomaScan® aptamer-based assay. We used sparse distance-weighted discrimination (sDWD) to test for a difference in protein expression and permutation tests to identify univariate associations between proteins and forced expiratory volume in 1 s % predicted (FEV₁ % pred). Significant proteins were entered into a pathway over-representation analysis. We also constructed protein-driven endotypes using K-means clustering and performed over-representation analysis on the proteins that were significantly different between clusters. We compared protein-associated clusters to those obtained from BALF and plasma metabolomics data on the same patient cohort.

Results

After filtering, we retained 3872 proteins for further analysis. Based on sDWD, protein expression was able to separate cases and controls. We found 575 proteins that were significantly correlated with FEV₁ % pred after multiple comparisons adjustment. We identified two protein-driven endotypes, one of which was associated with poor lung function, and found that insulin and apoptosis pathways were differentially represented. We found similar clusters driven by metabolomics in BALF but not plasma.

Conclusion

Protein expression differs in persons living with HIV with and without obstructive lung disease. We were not able to identify specific pathways differentially expressed among patients based on FEV₁ % pred; however, we identified a unique protein endotype associated with insulin and apoptotic pathways.

Lung and Plasma Metabolome in HIV-Associated Obstructive Lung Disease

BACKGROUND

HIV is a risk factor for obstructive lung disease (OLD), independent of smoking. We used mass spectrometry (MS) approaches to identify metabolomic biomarkers that inform mechanistic pathogenesis of OLD in persons with HIV (PWH).

METHODS

We obtained bronchoalveolar lavage fluid (BALF) samples from 52 PWH, in case:control (+OLD/-OLD) pairs matched on age, smoking status, and antiretroviral treatment. Four hundred nine metabolites from 8 families were measured on BALF and plasma samples using a MS-based Biocrates platform. After filtering metabolites with a high proportion of missing values and values below the level of detection, we performed univariate testing using paired t tests followed by false discovery rate corrections. We used distance-weighted discrimination (DWD) to test for an overall difference in the metabolite profile between cases and controls.

RESULTS

After filtering, there were 252 BALF metabolites for analysis from 8 metabolite families. DWD testing found that collectively, BALF metabolites differentiated cases from controls, whereas plasma metabolites did not. In BALF samples, we identified 3 metabolites that correlated with OLD at the false discovery rate of 10%; all were in the phosphatidylcholine family. We identified additional BALF metabolites when analyzing lung function as a continuous variable, and these included acylcarnitines, triglycerides, and a cholesterol ester.

CONCLUSIONS

Collectively, BALF metabolites differentiate PWH with and without OLD. These included several BALF lipid metabolites. These findings were limited to BALF and were not found in plasma from the same individuals. Phosphatidylcholine, the most common lipid component of surfactant, was the predominant lipid metabolite differentially expressed.

Catching the Wave: Detecting Strain-Specific SARS-CoV-2 Peptides in Clinical Samples Collected during Infection Waves from Diverse Geographical Locations

The Coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resulted in a major health crisis worldwide with its continuously emerging new strains, resulting in new viral variants that drive "waves" of infection. PCR or antigen detection assays have been routinely used to detect clinical infections; however, the emergence of these newer strains has presented challenges in detection. One of the alternatives has been to detect and characterize variant-specific peptide sequences from viral proteins using mass spectrometry (MS)-based methods. MS methods can potentially help in both diagnostics and vaccine development by understanding the dynamic changes in the viral proteome associated with specific strains and infection waves. In this study, we developed an accessible, flexible, and shareable bioinformatics workflow that was implemented in the Galaxy Platform to detect variant-specific peptide sequences from MS data derived from the clinical samples. We demonstrated the utility of the workflow by characterizing published clinical data from across the world during various pandemic waves. Our analysis identified six SARS-CoV-2 variant-specific peptides suitable for confident detection by MS in commonly collected clinical samples.

The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2022 update

Galaxy is a mature, browser accessible workbench for scientific computing. It enables scientists to share, analyze and visualize their own data, with minimal technical impediments. A thriving global community continues to use, maintain and contribute to the project, with support from multiple national infrastructure providers that enable freely accessible analysis and training services. The Galaxy Training Network supports free, self-directed, virtual training with >230 integrated tutorials. Project engagement metrics have continued to grow over the last 2 years, including source code contributions, publications, software packages wrapped as tools, registered users and their daily analysis jobs, and new independent specialized servers. Key Galaxy technical developments include an improved user interface for launching large-scale analyses with many files, interactive tools for exploratory data analysis, and a complete suite of machine learning tools. Important scientific developments enabled by Galaxy include Vertebrate Genome Project (VGP) assembly workflows and global SARS-CoV-2 collaborations.

Quantitative Proteogenomic Characterization of Inflamed Murine Colon Tissue Using an Integrated Discovery, Verification, and Validation Proteogenomic Workflow

Chronic inflammation of the colon causes genomic and/or transcriptomic events, which can lead to expression of non-canonical protein sequences contributing to oncogenesis. To better understand these mechanisms, Rag2^−/−Il10^−/− mice were infected with Helicobacter hepaticus to induce chronic inflammation of the cecum and the colon. Transcriptomic data from harvested proximal colon samples were used to generate a customized FASTA database containing non-canonical protein sequences. Using a proteogenomic approach, mass spectrometry data for proximal colon proteins were searched against this custom FASTA database using the Galaxy for Proteomics (Galaxy-P) platform. In addition to the increased abundance in inflammatory response proteins, we also discovered several non-canonical peptide sequences derived from unique proteoforms. We confirmed the veracity of these novel sequences using an automated bioinformatics verification workflow with targeted MS-based assays for peptide validation. Our bioinformatics discovery workflow identified 235 putative non-canonical peptide sequences, of which 58 were verified with high confidence and 39 were validated in targeted proteomics assays. This study provides insights into challenges faced when identifying non-canonical peptides using a proteogenomics approach and demonstrates an integrated workflow addressing these challenges. Our bioinformatic discovery and verification workflow is publicly available and accessible via the Galaxy platform and should be valuable in non-canonical peptide identification using proteogenomics.

An overview of technologies for MS-based proteomics-centric multi-omics

INTRODUCTION

Mass spectrometry-based proteomics reveals dynamic molecular signatures underlying phenotypes reflecting normal and perturbed conditions in living systems. Although valuable on its own, the proteome has only one level of moleclar information, with the genome, epigenome, transcriptome, and metabolome, all providing complementary information. Multi-omic analysis integrating information from one or more of these other domains with proteomic information provides a more complete picture of molecular contributors to dynamic biological systems.

AREAS COVERED

Here, we discuss the improvements to mass spectrometry-based technologies, focused on peptide-based, bottom-up approaches that have enabled deep, quantitative characterization of complex proteomes. These advances are facilitating the integration of proteomics data with other 'omic information, providing a more complete picture of living systems. We also describe the current state of bioinformatics software and approaches for integrating proteomics and other 'omics data, critical for enabling new discoveries driven by multi-omics.

EXPERT COMMENTARY

Multi-omics, centered on the integration of proteomics information with other 'omic information, has tremendous promise for biological and biomedical studies. Continued advances in approaches for generating deep, reliable proteomic data and bioinformatics tools aimed at integrating data across 'omic domains will ensure the discoveries offered by these multi-omic studies continue to increase.

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.

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.

Perspectives on automated composition of workflows in the life sciences

Scientific data analyses often combine several computational tools in automated pipelines, or workflows. Thousands of such workflows have been used in the life sciences, though their composition has remained a cumbersome manual process due to a lack of standards for annotation, assembly, and implementation. Recent technological advances have returned the long-standing vision of automated workflow composition into focus. This article summarizes a recent Lorentz Center workshop dedicated to automated composition of workflows in the life sciences. We survey previous initiatives to automate the composition process, and discuss the current state of the art and future perspectives. We start by drawing the "big picture" of the scientific workflow development life cycle, before surveying and discussing current methods, technologies and practices for semantic domain modelling, automation in workflow development, and workflow assessment. Finally, we derive a roadmap of individual and community-based actions to work toward the vision of automated workflow development in the forthcoming years. A central outcome of the workshop is a general description of the workflow life cycle in six stages: 1) scientific question or hypothesis, 2) conceptual workflow, 3) abstract workflow, 4) concrete workflow, 5) production workflow, and 6) scientific results. The transitions between stages are facilitated by diverse tools and methods, usually incorporating domain knowledge in some form. Formal semantic domain modelling is hard and often a bottleneck for the application of semantic technologies. However, life science communities have made considerable progress here in recent years and are continuously improving, renewing interest in the application of semantic technologies for workflow exploration, composition and instantiation. Combined with systematic benchmarking with reference data and large-scale deployment of production-stage workflows, such technologies enable a more systematic process of workflow development than we know today. We believe that this can lead to more robust, reusable, and sustainable workflows in the future.

Novel 4-Hydroxybenzyl Adducts in Human Hemoglobin: Structures and Mechanisms of Formation

Humans are exposed to large numbers of electrophiles from their diet, the environment, and endogenous physiological processes. Adducts formed at the N-terminal valine of hemoglobin are often used as biomarkers of human exposure to electrophilic compounds. We previously reported the formation of hemoglobin N-terminal valine adducts (added mass, 106.042 Da) in the blood of human smokers and nonsmokers and identified their structure as 4-hydroxybenzyl-Val. In the present work, mass spectrometry-based proteomics was utilized to identify additional sites for 4-hydroxybenzyl adduct formation at internal nucleophilic amino acid side chains within hemoglobin. Hemoglobin isolated from human blood was treated with para-quinone methide (para-QM) followed by global nanoLC-MS/MS and targeted nanoLC-MS/MS to identify amino acid residues containing the 4-hydroxybenzyl modification. Our experiments revealed the formation of 4-hydroxybenzyl adducts at the αHis20, αTyr24, αTyr42, αHis45, βSer72, βThr84, βThr87, βSer89, βHis92, βCys93, βCys112, βThr123, and βHis143 residues (in addition to N-terminal valine) through characteristic MS/MS spectra. These amino acid side chains had variable reactivity toward para-QM with αHis45, αTyr42, βCys93, βHis92, and βSer72 forming the largest numbers of adducts upon exposure to para-QM. Two additional mechanisms for formation of 4-hydroxybenzyl adducts in humans were investigated: exposure to 4-hydroxybenzaldehyde (4-HBA) followed by reduction and UV-mediated reactions of hemoglobin with tyrosine. Exposure of hemoglobin to a 5-fold molar excess of 4-HBA followed by reduction with sodium cyanoborohydride produced 4-hydroxybenzyl adducts at several amino acid side chains of which αHis20, αTyr24, αTyr42, αHis45, βSer44, βThr84, and βHis92 were verified in targeted mass spectrometry experiments. Similarly, exposure of human blood to ultraviolet radiation produced 4-hydroxybenzyl adducts at αHis20, αTyr24, αTyr42, αHis45, βSer44, βThr84, and βSer89. Overall, our results reveal that 4-hydroxybenzyl adducts form at multiple nucleophilic sites of hemoglobin and that para-QM is the most likely source of these adducts in humans.

A rigorous evaluation of optimal peptide targets for MS-based clinical diagnostics of Coronavirus Disease 2019 (COVID-19)

BACKGROUND

The Coronavirus Disease 2019 (COVID-19) global pandemic has had a profound, lasting impact on the world's population. A key aspect to providing care for those with COVID-19 and checking its further spread is early and accurate diagnosis of infection, which has been generally done via methods for amplifying and detecting viral RNA molecules. Detection and quantitation of peptides using targeted mass spectrometry-based strategies has been proposed as an alternative diagnostic tool due to direct detection of molecular indicators from non-invasively collected samples as well as the potential for high-throughput analysis in a clinical setting; many studies have revealed the presence of viral peptides within easily accessed patient samples. However, evidence suggests that some viral peptides could serve as better indicators of COVID-19 infection status than others, due to potential misidentification of peptides derived from human host proteins, poor spectral quality, high limits of detection etc. METHODS: In this study we have compiled a list of 636 peptides identified from Sudden Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) samples, including from in vitro and clinical sources. These datasets were rigorously analyzed using automated, Galaxy-based workflows containing tools such as PepQuery, BLAST-P, and the Multi-omic Visualization Platform as well as the open-source tools MetaTryp and Proteomics Data Viewer (PDV).

RESULTS

Using PepQuery for confirming peptide spectrum matches, we were able to narrow down the 639-peptide possibilities to 87 peptides that were most robustly detected and specific to the SARS-CoV-2 virus. The specificity of these sequences to coronavirus taxa was confirmed using Unipept and BLAST-P. Through stringent p-value cutoff combined with manual verification of peptide spectrum match quality, 4 peptides derived from the nucleocapsid phosphoprotein and membrane protein were found to be most robustly detected across all cell culture and clinical samples, including those collected non-invasively.

CONCLUSION

We propose that these peptides would be of the most value for clinical proteomics applications seeking to detect COVID-19 from patient samples. We also contend that samples harvested from the upper respiratory tract and oral cavity have the highest potential for diagnosis of SARS-CoV-2 infection from easily collected patient samples using mass spectrometry-based proteomics assays.

Updates on metaQuantome Software for Quantitative Metaproteomics

metaQuantome is a software suite that enables the quantitative analysis, statistical evaluation. and visualization of mass-spectrometry-based metaproteomics data. In the latest update of this software, we have provided several extensions, including a step-by-step training guide, the ability to perform statistical analysis on samples from multiple conditions, and a comparative analysis of metatranscriptomics data. The training module, accessed via the Galaxy Training Network, will help users to use the suite effectively both for functional as well as for taxonomic analysis. We extend the ability of metaQuantome to now perform multi-data-point quantitative and statistical analyses so that studies with measurements across multiple conditions, such as time-course studies, can be analyzed. With an eye on the multiomics analysis of microbial communities, we have also initiated the use of metaQuantome statistical and visualization tools on outputs from metatranscriptomics data, which complements the metagenomic and metaproteomic analyses already available. For this, we have developed a tool named MT2MQ ("metatranscriptomics to metaQuantome"), which takes in outputs from the ASaiM metatranscriptomics workflow and transforms them so that the data can be used as an input for comparative statistical analysis and visualization via metaQuantome. We believe that these improvements to metaQuantome will facilitate the use of the software for quantitative metaproteomics and metatranscriptomics and will enable multipoint data analysis. These improvements will take us a step toward integrative multiomic microbiome analysis so as to understand dynamic taxonomic and functional responses of these complex systems in a variety of biological contexts. The updated metaQuantome and MT2MQ are open-source software and are available via the Galaxy Toolshed and GitHub.

A rigorous evaluation of optimal peptide targets for MS-based clinical diagnostics of Coronavirus Disease 2019 (COVID-19)

The Coronavirus Disease 2019 (COVID-19) global pandemic has had a profound, lasting impact on the world's population. A key aspect to providing care for those with COVID-19 and checking its further spread is early and accurate diagnosis of infection, which has been generally done via methods for amplifying and detecting viral RNA molecules. Detection and quantitation of peptides using targeted mass spectrometry-based strategies has been proposed as an alternative diagnostic tool due to direct detection of molecular indicators from non-invasively collected samples as well as the potential for high-throughput analysis in a clinical setting; many studies have revealed the presence of viral peptides within easily accessed patient samples. However, evidence suggests that some viral peptides could serve as better indicators of COVID-19 infection status than others, due to potential misidentification of peptides derived from human host proteins, poor spectral quality, high limits of detection etc. In this study we have compiled a list of 639 peptides identified from Sudden Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) samples, including from in vitro and clinical sources. These datasets were rigorously analyzed using automated, Galaxy-based workflows containing tools such as PepQuery, BLAST-P, and the Multi-omic Visualization Platform as well as the open-source tools MetaTryp and Proteomics Data Viewer (PDV). Using PepQuery for confirming peptide spectrum matches, we were able to narrow down the 639 peptide possibilities to 87 peptides which were most robustly detected and specific to the SARS-CoV-2 virus. The specificity of these sequences to coronavirus taxa was confirmed using Unipept and BLAST-P. Applying stringent statistical scoring thresholds, combined with manual verification of peptide spectrum match quality, 4 peptides derived from the nucleocapsid phosphoprotein and membrane protein were found to be most robustly detected across all cell culture and clinical samples, including those collected non-invasively. We propose that these peptides would be of the most value for clinical proteomics applications seeking to detect COVID-19 from a variety of sample types. We also contend that samples taken from the upper respiratory tract and oral cavity have the highest potential for diagnosis of SARS-CoV-2 infection from easily collected patient samples using mass spectrometry-based proteomics assays.

ASaiM-MT: a validated and optimized ASaiM workflow for metatranscriptomics analysis within Galaxy framework

The Earth Microbiome Project (EMP) aided in understanding the role of microbial communities and the influence of collective genetic material (the 'microbiome') and microbial diversity patterns across the habitats of our planet. With the evolution of new sequencing technologies, researchers can now investigate the microbiome and map its influence on the environment and human health. Advances in bioinformatics methods for next-generation sequencing (NGS) data analysis have helped researchers to gain an in-depth knowledge about the taxonomic and genetic composition of microbial communities. Metagenomic-based methods have been the most commonly used approaches for microbiome analysis; however, it primarily extracts information about taxonomic composition and genetic potential of the microbiome under study, lacking quantification of the gene products (RNA and proteins). On the other hand, metatranscriptomics, the study of a microbial community's RNA expression, can reveal the dynamic gene expression of individual microbial populations and the community as a whole, ultimately providing information about the active pathways in the microbiome.  In order to address the analysis of NGS data, the ASaiM analysis framework was previously developed and made available via the Galaxy platform. Although developed for both metagenomics and metatranscriptomics, the original publication demonstrated the use of ASaiM only for metagenomics, while thorough testing for metatranscriptomics data was lacking.  In the current study, we have focused on validating and optimizing the tools within ASaiM for metatranscriptomics data. As a result, we deliver a robust workflow that will enable researchers to understand dynamic functional response of the microbiome in a wide variety of metatranscriptomics studies. This improved and optimized ASaiM-metatranscriptomics (ASaiM-MT) workflow is publicly available via the ASaiM framework, documented and supported with training material so that users can interrogate and characterize metatranscriptomic data, as part of larger meta-omic studies of microbiomes.

ASaiM-MT: a validated and optimized ASaiM workflow for metatranscriptomics analysis within Galaxy framework

The Human Microbiome Project (HMP) aided in understanding the role of microbial communities and the influence of collective genetic material (the 'microbiome') in human health and disease. With the evolution of new sequencing technologies, researchers can now investigate the microbiome and map its influence on human health. Advances in bioinformatics methods for next-generation sequencing (NGS) data analysis have helped researchers to gain an in-depth knowledge about the taxonomic and genetic composition of microbial communities. Metagenomic-based methods have been the most commonly used approaches for microbiome analysis; however, it primarily extracts information about taxonomic composition and genetic potential of the microbiome under study, lacking quantification of the gene products (RNA and proteins). Conversely, metatranscriptomics, the study of a microbial community's RNA expression, can reveal the dynamic gene expression of individual microbial populations and the community as a whole, ultimately providing information about the active pathways in the microbiome.  In order to address the analysis of NGS data, the ASaiM analysis framework was previously developed and made available via the Galaxy platform. Although developed for both metagenomics and metatranscriptomics, the original publication demonstrated the use of ASaiM only for metagenomics, while thorough testing for metatranscriptomics data was lacking.  In the current study, we have focused on validating and optimizing the tools within ASaiM for metatranscriptomics data. As a result, we deliver a robust workflow that will enable researchers to understand dynamic functional response of the microbiome in a wide variety of metatranscriptomics studies. This improved and optimized ASaiM-metatranscriptomics (ASaiM-MT) workflow is publicly available via the ASaiM framework, documented and supported with training material so that users can interrogate and characterize metatranscriptomic data, as part of larger meta-omic studies of microbiomes.

Metaproteomics Analysis of SARS-CoV-2-Infected Patient Samples Reveals Presence of Potential Coinfecting Microorganisms

In this Letter, we reanalyze published mass spectrometry data sets of clinical samples with a focus on determining the coinfection status of individuals infected with SARS-CoV-2 coronavirus. We demonstrate the use of ComPIL 2.0 software along with a metaproteomics workflow within the Galaxy platform to detect cohabitating potential pathogens in COVID-19 patients using mass spectrometry-based analysis. From a sample collected from gargling solutions, we detected Streptococcus pneumoniae (opportunistic and multidrug-resistant pathogen) and Lactobacillus rhamnosus (a probiotic component) along with SARS-Cov-2. We could also detect Pseudomonas sps. Bc-h from COVID-19 positive samples and Acinetobacter ursingii and Pseudomonas monteilii from COVID-19 negative samples collected from oro- and nasopharyngeal samples. We believe that the early detection and characterization of coinfections by using metaproteomics from COVID-19 patients will potentially impact the diagnosis and treatment of patients affected by SARS-CoV-2 infection.

Evaluation of the potential of Pap test fluid and cervical swabs to serve as clinical diagnostic biospecimens for the detection of ovarian cancer by mass spectrometry-based proteomics

BACKGROUND

The purpose of this study was to determine whether the residual fixative from a liquid-based Pap test or a swab of the cervix contained proteins that were also found in the primary tumor of a woman with high grade serous ovarian cancer. This study is the first step in determining the feasibility of using the liquid-based Pap test or a cervical swab for the detection of ovarian cancer protein biomarkers.

METHODS

Proteins were concentrated by acetone precipitation from the cell-free supernatant of the liquid-based Pap test fixative or eluted from the cervical swab. Protein was also extracted from the patient's tumor tissue. The protein samples were digested into peptides with trypsin, then the peptides were run on 2D-liquid chromatography mass spectrometry (2D-LCMS). The data was searched against a human protein database for the identification of peptides and proteins in each biospecimen. The proteins that were identified were classified for cellular localization and molecular function by bioinformatics integration.

RESULTS

We identified almost 5000 proteins total in the three matched biospecimens. More than 2000 proteins were expressed in each of the three biospecimens, including several known ovarian cancer biomarkers such as CA125, HE4, and mesothelin. By Scaffold analysis of the protein Gene Ontology categories and functional analysis using PANTHER, the proteins were classified by cellular localization and molecular function, demonstrating that the Pap test fluid and cervical swab proteins are similar to each other, and also to the tumor extract.

CONCLUSIONS

Our results suggest that Pap test fixatives and cervical swabs are a rich source of tumor-specific biomarkers for ovarian cancer, which could be developed as a test for ovarian cancer detection.

Application of Proteomics in Sarcoidosis

Sarcoidosis is a multisystem disease with heterogeneity in manifestations and outcomes. System-level studies leveraging "omics" technologies are expected to define mechanisms contributing to sarcoidosis heterogeneous manifestations and course. With improvements in mass spectrometry (MS) and bioinformatics, it is possible to study protein abundance for a large number of proteins simultaneously. Contemporary fast-scanning MS enables the acquisition of spectral data for deep coverage of the proteins with data-dependent or data-independent acquisition MS modes. Studies leveraging MS-based proteomics in sarcoidosis have characterized BAL fluid (BALF), alveolar macrophages, plasma, and exosomes. These studies identified several differentially expressed proteins, including protocadherin-2 precursor, annexin A2, pulmonary surfactant A2, complement factors C3, vitamin-D-binding protein, cystatin B, and amyloid P, comparing subjects with sarcoidosis with control subjects. Other studies identified ceruloplasmin, complement factors B, C3, and 1, and others with differential abundance in sarcoidosis compared with other interstitial lung diseases. Using quantitative proteomics, most recent studies found differences in PI3K/Akt/mTOR, MAP kinase, pluripotency-associated transcriptional factor, and hypoxia response pathways. Other studies identified increased clathrin-mediated endocytosis and Fcγ receptor-mediated phagocytosis pathways in sarcoidosis alveolar macrophages. Although studies in mixed BAL and blood cells or plasma are limited, some of the changes in lung compartment are detected in the blood cells and plasma. We review proteomics for sarcoidosis with a focus on the existing MS data acquisition strategies, bioinformatics for spectral data analysis to infer protein identity and quantity, unique aspects about biospecimen collection and processing for lung-related proteomics, and proteomics studies conducted to date in sarcoidosis.

Survey of metaproteomics software tools for functional microbiome analysis

To gain a thorough appreciation of microbiome dynamics, researchers characterize the functional relevance of expressed microbial genes or proteins. This can be accomplished through metaproteomics, which characterizes the protein expression of microbiomes. Several software tools exist for analyzing microbiomes at the functional level by measuring their combined proteome-level response to environmental perturbations. In this survey, we explore the performance of six available tools, to enable researchers to make informed decisions regarding software choice based on their research goals. Tandem mass spectrometry-based proteomic data obtained from dental caries plaque samples grown with and without sucrose in paired biofilm reactors were used as representative data for this evaluation. Microbial peptides from one sample pair were identified by the X! tandem search algorithm via SearchGUI and subjected to functional analysis using software tools including eggNOG-mapper, MEGAN5, MetaGOmics, MetaProteomeAnalyzer (MPA), ProPHAnE, and Unipept to generate functional annotation through Gene Ontology (GO) terms. Among these software tools, notable differences in functional annotation were detected after comparing differentially expressed protein functional groups. Based on the generated GO terms of these tools we performed a peptide-level comparison to evaluate the quality of their functional annotations. A BLAST analysis against the NCBI non-redundant database revealed that the sensitivity and specificity of functional annotation varied between tools. For example, eggNOG-mapper mapped to the most number of GO terms, while Unipept generated more accurate GO terms. Based on our evaluation, metaproteomics researchers can choose the software according to their analytical needs and developers can use the resulting feedback to further optimize their algorithms. To make more of these tools accessible via scalable metaproteomics workflows, eggNOG-mapper and Unipept 4.0 were incorporated into the Galaxy platform.

Novel protein pathways in development and progression of pulmonary sarcoidosis

Pulmonary involvement occurs in up to 95% of sarcoidosis cases. In this pilot study, we examine lung compartment-specific protein expression to identify pathways linked to development and progression of pulmonary sarcoidosis. We characterized bronchoalveolar lavage (BAL) cells and fluid (BALF) proteins in recently diagnosed sarcoidosis cases. We identified 4,306 proteins in BAL cells, of which 272 proteins were differentially expressed in sarcoidosis compared to controls. These proteins map to novel pathways such as integrin-linked kinase and IL-8 signaling and previously implicated pathways in sarcoidosis, including phagosome maturation, clathrin-mediated endocytic signaling and redox balance. In the BALF, the differentially expressed proteins map to several pathways identified in the BAL cells. The differentially expressed BALF proteins also map to aryl hydrocarbon signaling, communication between innate and adaptive immune response, integrin, PTEN and phospholipase C signaling, serotonin and tryptophan metabolism, autophagy, and B cell receptor signaling. Additional pathways that were different between progressive and non-progressive sarcoidosis in the BALF included CD28 signaling and PFKFB4 signaling. Our studies demonstrate the power of contemporary proteomics to reveal novel mechanisms operational in sarcoidosis. Application of our workflows in well-phenotyped large cohorts maybe beneficial to identify biomarkers for diagnosis and prognosis and therapeutically tenable molecular mechanisms.

Multi-omics Visualization Platform: An extensible Galaxy plug-in for multi-omics data visualization and exploration

BACKGROUND

Proteogenomics integrates genomics, transcriptomics, and mass spectrometry (MS)-based proteomics data to identify novel protein sequences arising from gene and transcript sequence variants. Proteogenomic data analysis requires integration of disparate 'omic software tools, as well as customized tools to view and interpret results. The flexible Galaxy platform has proven valuable for proteogenomic data analysis. Here, we describe a novel Multi-omics Visualization Platform (MVP) for organizing, visualizing, and exploring proteogenomic results, adding a critically needed tool for data exploration and interpretation.

FINDINGS

MVP is built as an HTML Galaxy plug-in, primarily based on JavaScript. Via the Galaxy API, MVP uses SQLite databases as input-a custom data type (mzSQLite) containing MS-based peptide identification information, a variant annotation table, and a coding sequence table. Users can interactively filter identified peptides based on sequence and data quality metrics, view annotated peptide MS data, and visualize protein-level information, along with genomic coordinates. Peptides that pass the user-defined thresholds can be sent back to Galaxy via the API for further analysis; processed data and visualizations can also be saved and shared. MVP leverages the Integrated Genomics Viewer JavaScript framework, enabling interactive visualization of peptides and corresponding transcript and genomic coding information within the MVP interface.

CONCLUSIONS

MVP provides a powerful, extensible platform for automated, interactive visualization of proteogenomic results within the Galaxy environment, adding a unique and critically needed tool for empowering exploration and interpretation of results. The platform is extensible, providing a basis for further development of new functionalities for proteogenomic data visualization.

The Cancer Microbiome: Distinguishing Direct and Indirect Effects Requires a Systemic View

The collection of microbes that live in and on the human body - the human microbiome - can impact on cancer initiation, progression, and response to therapy, including cancer immunotherapy. The mechanisms by which microbiomes impact on cancers can yield new diagnostics and treatments, but much remains unknown. The interactions between microbes, diet, host factors, drugs, and cell-cell interactions within the cancer itself likely involve intricate feedbacks, and no single component can explain all the behavior of the system. Understanding the role of host-associated microbial communities in cancer systems will require a multidisciplinary approach combining microbial ecology, immunology, cancer cell biology, and computational biology - a systems biology approach.

Abstract DP-002: COMPARISON OF POTENTIAL OVARIAN CANCER BIOMARKERS BY MASS SPECTROMETRY-BASED PROTEOMIC ANALYSIS OF RESIDUAL PAP TEST FLUID, CERVICAL SWABS, AND TUMOR TISSUE FROM AN OVARIAN CANCER PATIENT

Early detection is the key to increased survival for women with ovarian cancer, yet a screening tool has yet to be developed that is adequately sensitive and specific enough for use in the general population. In contrast, screening for cervical cancer by Pap tests has been routinely performed for over 50 years. In the liquid-based Pap test, cells are collected from the cervix and placed into an alcohol-based fixative and then examined for abnormal cells. Since ovarian cancer cells have been observed in Pap tests, we reasoned that ovarian cancer peptide biomarkers may also be present. Our central hypothesis is that proteins shed by ovarian cancer cells can be detected during routine Pap tests by mass spectrometry (MS)-based proteomics. In particular, when collected at the time of cervical cancer Pap test screening, the alcohol-based Pap test fixative and cervical swabs are ideal for biomarker discovery since they are derived from a site near the ovarian cancer (i.e. proteins may be secreted or shed from the tumor and flow through the fallopian tube into the uterus and out the cervical opening). Recently, the fimbria of the fallopian tube have been suggested to be the true precursor to ovarian cancer, strengthening our hypothesis that ovarian cancer proteins will be found in the lower genital tract, perhaps even at early stages. To demonstrate the feasibility of using Pap tests as a biospecimen for proteomics, we previously examined the proteins present in residual Pap test fixative samples from women with normal cervical cytology by MS and described 152 proteins in the “Normal Pap Proteome.” The objective of this study was to identify and compare the proteins from three different sources from the same ovarian cancer patient: (i) the residual Pap test fixative, (ii) a Merocel swab of the cervix, and (iii) the primary ovarian cancer tumor tissue. Proteins were concentrated from the cell-free supernatant of the Pap test fixative or eluted from the swab, and then trypsin digested using the filter-aided sample preparation method. A total protein extract from the patient's tumor tissue was digested by standard in-solution trypsin digestion. The samples were run on 2D-liquid chromatography MS/MS, followed by bioinformatics integration. We identified over 5000 proteins total in the three samples. More than 2000 proteins were expressed in all three ovarian cancer samples, including several known ovarian cancer biomarkers such as CA125. By Scaffold analysis of the Gene Ontology nomenclature of the proteins, we classified the proteins by both cellular localization and biological processes. Additional matched samples from patients will be used to build a library of proteins and peptides that are specific to ovarian cancer for use in the development of targeted MS assays. We conclude that quantification of proteins from Pap test fixatives and cervical swabs will prove to be a rich source of biomarkers for ovarian cancer detection.

Citation Format: Kristin L.M. Boylan, Anna C. Rogers, Melissa A. Geller, Peter A. Argenta, Timothy J. Griffin, and Amy P.N. Skubitz. COMPARISON OF POTENTIAL OVARIAN CANCER BIOMARKERS BY MASS SPECTROMETRY-BASED PROTEOMIC ANALYSIS OF RESIDUAL PAP TEST FLUID, CERVICAL SWABS, AND TUMOR TISSUE FROM AN OVARIAN CANCER PATIENT [abstract]. In: Proceedings of the 12th Biennial Ovarian Cancer Research Symposium; Sep 13-15, 2018; Seattle, WA. Philadelphia (PA): AACR; Clin Cancer Res 2019;25(22 Suppl):Abstract nr DP-002.

metaQuantome: An Integrated, Quantitative Metaproteomics Approach Reveals Connections Between Taxonomy and Protein Function in Complex Microbiomes

Microbiome research offers promising insights into the impact of microorganisms on biological systems. Metaproteomics, the study of microbial proteins at the community level, integrates genomic, transcriptomic, and proteomic data to determine the taxonomic and functional state of a microbiome. However, standard metaproteomics software is subject to several limitations, commonly supporting only spectral counts, emphasizing exploratory analysis rather than hypothesis testing and rarely offering the ability to analyze the interaction of function and taxonomy - that is, which taxa are responsible for different processes.Here we present metaQuantome, a novel, multifaceted software suite that analyzes the state of a microbiome by leveraging complex taxonomic and functional hierarchies to summarize peptide-level quantitative information, emphasizing label-free intensity-based methods. For experiments with multiple experimental conditions, metaQuantome offers differential abundance analysis, principal components analysis, and clustered heat map visualizations, as well as exploratory analysis for a single sample or experimental condition. We benchmark metaQuantome analysis against standard methods, using two previously published datasets: (1) an artificially assembled microbial community dataset (taxonomy benchmarking) and (2) a dataset with a range of recombinant human proteins spiked into an Escherichia coli background (functional benchmarking). Furthermore, we demonstrate the use of metaQuantome on a previously published human oral microbiome dataset.In both the taxonomic and functional benchmarking analyses, metaQuantome quantified taxonomic and functional terms more accurately than standard summarization-based methods. We use the oral microbiome dataset to demonstrate metaQuantome's ability to produce publication-quality figures and elucidate biological processes of the oral microbiome. metaQuantome enables advanced investigation of metaproteomic datasets, which should be broadly applicable to microbiome-related research. In the interest of accessible, flexible, and reproducible analysis, metaQuantome is open source and available on the command line and in Galaxy.

Distinct Cancer-Promoting Stromal Gene Expression Depending on Lung Function

Rationale: Chronic obstructive pulmonary disease is an independent risk factor for lung cancer, but the underlying molecular mechanisms are unknown. We hypothesized that lung stromal cells activate pathological gene expression programs that support oncogenesis.Objectives: To identify molecular mechanisms operating in the lung stroma that support the development of lung cancer.Methods: The study included subjects with and without lung cancer across a spectrum of lung-function values. We conducted a multiomics analysis of nonmalignant lung tissue to quantify the transcriptome, translatome, and proteome.Measurements and Main Results: Cancer-associated gene expression changes predominantly manifested as alterations in the efficiency of mRNA translation modulating protein levels in the absence of corresponding changes in mRNA levels. The molecular mechanisms that drove these cancer-associated translation programs differed based on lung function. In subjects with normal to mildly impaired lung function, the mammalian target of rapamycin (mTOR) pathway served as an upstream driver, whereas in subjects with severe airflow obstruction, pathways downstream of pathological extracellular matrix emerged. Consistent with a role during cancer initiation, both the mTOR and extracellular matrix gene expression programs paralleled the activation of previously identified procancer secretomes. Furthermore, an in situ examination of lung tissue showed that stromal fibroblasts expressed cancer-associated proteins from two procancer secretomes: one that included IL-6 (in cases of mild or no airflow obstruction), and one that included BMP1 (in cases of severe airflow obstruction).Conclusions: Two distinct stromal gene expression programs that promote cancer initiation are activated in patients with lung cancer depending on lung function. Our work has implications both for screening strategies and for personalized approaches to cancer treatment.

Improve your Galaxy text life: The Query Tabular Tool

Galaxy provides an accessible platform where multi-step data analysis workflows integrating disparate software can be run, even by researchers with limited programming expertise. Applications of such sophisticated workflows are many, including those which integrate software from different ‘omic domains (e.g. genomics, proteomics, metabolomics). In these complex workflows, intermediate outputs are often generated as tabular text files, which must be transformed into customized formats which are compatible with the next software tools in the pipeline. Consequently, many text manipulation steps are added to an already complex workflow, overly complicating the process. In some cases, limitations to existing text manipulation are such that desired analyses can only be carried out using highly sophisticated processing steps beyond the reach of even advanced users and developers. For users with some SQL knowledge, these text operations could be combined into single, concise query on a relational database. As a solution, we have developed the Query Tabular Galaxy tool, which leverages a SQLite database generated from tabular input data. This database can be queried and manipulated to produce transformed and customized tabular outputs compatible with downstream processing steps. Regular expressions can also be utilized for even more sophisticated manipulations, such as find and replace and other filtering actions. Using several Galaxy-based multi-omic workflows as an example, we demonstrate how the Query Tabular tool dramatically streamlines and simplifies the creation of multi-step analyses, efficiently enabling complicated textual manipulations and processing. This tool should find broad utility for users of the Galaxy platform seeking to develop and use sophisticated workflows involving text manipulation on tabular outputs.

Improve your Galaxy text life: The Query Tabular Tool

Galaxy provides an accessible platform where multi-step data analysis workflows integrating disparate software can be run, even by researchers with limited programming expertise. Applications of such sophisticated workflows are many, including those which integrate software from different 'omic domains (e.g. genomics, proteomics, metabolomics). In these complex workflows, intermediate outputs are often generated as tabular text files, which must be transformed into customized formats which are compatible with the next software tools in the pipeline. Consequently, many text manipulation steps are added to an already complex workflow, overly complicating the process. In some cases, limitations to existing text manipulation are such that desired analyses can only be carried out using highly sophisticated processing steps beyond the reach of even advanced users and developers. For users with some SQL knowledge, these text operations could be combined into single, concise query on a relational database. As a solution, we have developed the Query Tabular Galaxy tool, which leverages a SQLite database generated from tabular input data. This database can be queried and manipulated to produce transformed and customized tabular outputs compatible with downstream processing steps. Regular expressions can also be utilized for even more sophisticated manipulations, such as find and replace and other filtering actions. Using several Galaxy-based multi-omic workflows as an example, we demonstrate how the Query Tabular tool dramatically streamlines and simplifies the creation of multi-step analyses, efficiently enabling complicated textual manipulations and processing. This tool should find broad utility for users of the Galaxy platform seeking to develop and use sophisticated workflows involving text manipulation on tabular outputs.

Plasma lipoproteome in Alzheimer's disease: a proof-of-concept study

Background

Although total plasma lipoproteome consists of proteins that have shown promises as biomarkers that can identify Alzheimer's disease (AD), effect sizes are modest. The objective of this study is to provide initial proof-of-concept that the plasma lipoproteome more likely differ between AD cases and controls when measured in individual plasma lipoprotein fractions than when measured as total in immunodepleted plasma.

Methods

We first developed a targeted proteomics method based on selected reaction monitoring (SRM) and liquid chromatography and tandem mass spectrometry for measurement of 120 tryptic peptides from 79 proteins that are commonly present in plasma lipoproteins. Then in a proof-of concept case-control study of 5 AD cases and 5 sex- and age-matched controls, we applied the targeted proteomic method and performed relatively quantification of 120 tryptic peptides in plasma lipoprotein fractions (fractionated by sequential gradient ultracentrifugation) and in immunodepleted plasma (of albumin and IgG). Unadjusted p values from two-sample t-tests and overall fold change was used to evaluate a peptide relative difference between AD cases and controls, with lower p values (< 0.05) or greater fold differences (> 1.05 or < 0.95) suggestive of greater peptide/protein differences.

Results

Within-day and between-days technical precisions (mean %CV [SD] of all SRM transitions) of the targeted proteomic method were 3.95% (2.65) and 9.31% (5.59), respectively. Between-days technical precisions (mean % CV [SD]) of the entire plasma lipoproteomic workflow including plasma lipoprotein fractionation was 27.90% (14.61). Ten tryptic peptides that belonged to 5 proteins in plasma lipoproteins had unadjusted p values < 0.05, compared to no peptides in immunodepleted plasma. Furthermore, 27, 32, 17, and 20 tryptic peptides in VLDL, IDL, LDL and HDL, demonstrated overall peptide fold differences > 1.05 or < 0.95, compared to only 6 tryptic peptides in immunodepleted plasma. The overall comparisons, therefore, suggested greater peptide/protein differences in plasma lipoproteome when measured in individual plasma lipoproteins than as total in immunodepleted plasma. Specifically, protein complement C3's peptide IHWESASLLR, had unadjusted p values of 0.00007, 0.00012, and 0.0006 and overall 1.25, 1.17, 1.14-fold changes in VLDL, IDL, and LDL, respectively. After positive False Discovery Rate (pFDR) adjustment, the complement C3 peptide IHWESASLLR in VLDL remained statistically different (adjusted p value < 0.05).

Discussion

The findings may warrant future studies to investigate plasma lipoproteome when measured in individual plasma lipoprotein fractions for AD diagnosis.

Quantitative proteomics of forebrain subcellular fractions in fragile X mental retardation 1 knockout mice following acute treatment with 2-Methyl-6-(phenylethynyl)pyridine: Relevance to developmental study of schizophrenia

The fragile X mental retardation 1 knockout (Fmr1 KO) mouse replicates behavioral deficits associated with autism, fragile X syndrome, and schizophrenia. Less is known whether protein expression changes are consistent with findings in subjects with schizophrenia. In the current study, we used liquid chromatography tandem mass spectrometry (LC-MS/MS) proteomics to determine the protein expression of four subcellular fractions in the forebrains of Fmr1 KO mice vs. C57BL/6 J mice and the effect of a negative allosteric modulator of mGluR5-2-Methyl-6-(phenylethynyl)pyridine (MPEP)-on protein expression. Strain- and treatment-specific differential expression of proteins was observed, many of which have previously been observed in the brains of subjects with schizophrenia. Western blotting verified the direction and magnitude of change for several proteins in different subcellular fractions as follows: neurofilament light protein (NEFL) and 2',3'-cyclic-nucleotide 3'-phosphodiesterase (CNP) in the total homogenate; heterogeneous nuclear ribonucleoproteins C1/C2 (HNRNPC) and heterogeneous nuclear ribonucleoprotein D0 (HNRNPD) in the nuclear fraction; excitatory amino acid transporter 2 (EAAT2) and ras-related protein rab 3a (RAB3A) in the synaptic fraction; and ras-related protein rab 35 (RAB35) and neuromodulin (GAP43) in the rough endoplasmic reticulum fraction. Individuals with FXS do not display symptoms of schizophrenia. However, the biomarkers that have been identified suggest that the Fmr1 KO model could potentially be useful in the study of schizophrenia.

Evaluating the potential of residual Pap test fluid as a resource for the metaproteomic analysis of the cervical-vaginal microbiome

The human cervical-vaginal area contains proteins derived from microorganisms that may prevent or predispose women to gynecological conditions. The liquid Pap test fixative is an unexplored resource for analysis of microbial communities and the microbe-host interaction. Previously, we showed that the residual cell-free fixative from discarded Pap tests of healthy women could be used for mass spectrometry (MS) based proteomic identification of cervical-vaginal proteins. In this study, we reprocessed these MS raw data files for metaproteomic analysis to characterize the microbial community composition and function of microbial proteins in the cervical-vaginal region. This was accomplished by developing a customized protein sequence database encompassing microbes likely present in the vagina. High-mass accuracy data were searched against the protein FASTA database using a two-step search method within the Galaxy for proteomics platform. Data was analyzed by MEGAN6 (MetaGenomeAnalyzer) for phylogenetic and functional characterization. We identified over 300 unique peptides from a variety of bacterial phyla and Candida. Peptides corresponding to proteins involved in carbohydrate metabolism, oxidation-reduction, and transport were identified. By identifying microbial peptides in Pap test supernatants it may be possible to acquire a functional signature of these microbes, as well as detect specific proteins associated with cervical health and disease.

Multi-omic molecular profiling of lung cancer in COPD

Chronic obstructive pulmonary disease (COPD) is a known risk factor for developing lung cancer but the underlying mechanisms remain unknown. We hypothesise that the COPD stroma contains molecular mechanisms supporting tumourigenesis.We conducted an unbiased multi-omic analysis to identify gene expression patterns that distinguish COPD stroma in patients with or without lung cancer. We obtained lung tissue from patients with COPD and lung cancer (tumour and adjacent non-malignant tissue) and those with COPD without lung cancer for profiling of proteomic and mRNA (both cytoplasmic and polyribosomal). We used the Joint and Individual Variation Explained (JIVE) method to integrate and analyse across the three datasets.JIVE identified eight latent patterns that robustly distinguished and separated the three groups of tissue samples (tumour, adjacent and control). Predictive variables that associated with the tumour, compared to adjacent stroma, were mainly represented in the transcriptomic data, whereas predictive variables associated with adjacent tissue, compared to controls, were represented at the translatomic level. Pathway analysis revealed extracellular matrix and phosphatidylinositol-4,5-bisphosphate 3-kinase-protein kinase B signalling pathways as important signals in the tumour adjacent stroma.The multi-omic approach distinguishes tumour adjacent stroma in lung cancer and reveals two stromal expression patterns associated with cancer.

ULK1 phosphorylates Ser30 of BECN1 in association with ATG14 to stimulate autophagy induction

ULK1 (unc51-like autophagy activating kinase 1) is a serine/threonine kinase that plays a key role in regulating macroautophagy/autophagy induction in response to amino acid starvation. Despite the recent progress in understanding ULK1 functions, the molecular mechanism by which ULK1 regulates the induction of autophagy remains elusive. In this study, we determined that ULK1 phosphorylates Ser30 of BECN1 (Beclin 1) in association with ATG14 (autophagy-related 14) but not with UVRAG (UV radiation resistance associated). The Ser30 phosphorylation was induced by deprivation of amino acids or treatments with Torin 1 or rapamycin, the conditions that inhibit MTORC1 (mechanistic target of rapamycin complex 1), and requires ATG13 and RB1CC1 (RB1 inducible coiled-coil 1), proteins that interact with ULK1. Hypoxia or glutamine deprivation, which inhibit MTORC1, was also able to increase the phosphorylation in a manner dependent upon ULK1 and ULK2. Blocking the BECN1 phosphorylation by replacing Ser30 with alanine suppressed the amino acid starvation-induced activation of the ATG14-containing PIK3C3/VPS34 (phosphatidylinositol 3-kinase catalytic subunit type 3) kinase, and reduced autophagy flux and the formation of phagophores and autophagosomes. The Ser30-to-Ala mutation did not affect the ULK1-mediated phosphorylations of BECN1 Ser15 or ATG14 Ser29, indicating that the BECN1 Ser30 phosphorylation might regulate autophagy independently of those 2 sites. Taken together, these results demonstrate that BECN1 Ser30 is a ULK1 target site whose phosphorylation activates the ATG14-containing PIK3C3 complex and stimulates autophagosome formation in response to amino acid starvation, hypoxia, and MTORC1 inhibition.

Disseminating Metaproteomic Informatics Capabilities and Knowledge Using the Galaxy-P Framework

The impact of microbial communities, also known as the microbiome, on human health and the environment is receiving increased attention. Studying translated gene products (proteins) and comparing metaproteomic profiles may elucidate how microbiomes respond to specific environmental stimuli, and interact with host organisms. Characterizing proteins expressed by a complex microbiome and interpreting their functional signature requires sophisticated informatics tools and workflows tailored to metaproteomics. Additionally, there is a need to disseminate these informatics resources to researchers undertaking metaproteomic studies, who could use them to make new and important discoveries in microbiome research. The Galaxy for proteomics platform (Galaxy-P) offers an open source, web-based bioinformatics platform for disseminating metaproteomics software and workflows. Within this platform, we have developed easily-accessible and documented metaproteomic software tools and workflows aimed at training researchers in their operation and disseminating the tools for more widespread use. The modular workflows encompass the core requirements of metaproteomic informatics: (a) database generation; (b) peptide spectral matching; (c) taxonomic analysis and (d) functional analysis. Much of the software available via the Galaxy-P platform was selected, packaged and deployed through an online metaproteomics "Contribution Fest" undertaken by a unique consortium of expert software developers and users from the metaproteomics research community, who have co-authored this manuscript. These resources are documented on GitHub and freely available through the Galaxy Toolshed, as well as a publicly accessible metaproteomics gateway Galaxy instance. These documented workflows are well suited for the training of novice metaproteomics researchers, through online resources such as the Galaxy Training Network, as well as hands-on training workshops. Here, we describe the metaproteomics tools available within these Galaxy-based resources, as well as the process by which they were selected and implemented in our community-based work. We hope this description will increase access to and utilization of metaproteomics tools, as well as offer a framework for continued community-based development and dissemination of cutting edge metaproteomics software.

Salivary proteomics of healthy dogs: An in depth catalog

Objective

To provide an in-depth catalog of the salivary proteome and endogenous peptidome of healthy dogs, evaluate proteins and peptides with antimicrobial properties, and compare the most common salivary proteins and peptides between different breed phylogeny groups.

Methods

36 healthy dogs without evidence of periodontal disease representing four breed phylogeny groups, based upon single nucleotide polymorphism haplotypes (ancient, herding/sighthound, and two miscellaneous groups). Saliva collected from dogs was pooled by phylogeny group and analyzed using nanoscale liquid chromatography-tandem mass spectrometry. Resulting tandem mass spectra were compared to databases for identification of endogenous peptides and inferred proteins.

Results

2,491 proteins and endogenous peptides were found in the saliva of healthy dogs with no periodontal disease. All dog phylogeny groups’ saliva was rich in proteins and peptides with antimicrobial functions. The ancient breeds group was distinct in that it contained unique proteins and was missing many proteins and peptides present in the other groups.

Conclusions and clinical relevance

Using a sophisticated nanoscale liquid chromatography-tandem mass spectrometry, we were able to identify 10-fold more salivary proteins than previously reported in dogs. Seven of the top 10 most abundant proteins or peptides serve immune functions and many more with various antimicrobial mechanisms were found. This is the most comprehensive analysis of healthy canine saliva to date, and will provide the groundwork for future studies analyzing salivary proteins and endogenous peptides in disease states.

Normalization of mass spectrometry data (NOMAD)

iTRAQ and TMT reagent-based mass spectrometry (MS) are commonly used technologies for quantitative proteomics in biological samples. Such studies are often performed over multiple MS runs, potentially resulting in introduction of MS run bias that could affect downstream analysis. Such MS data have therefore commonly been normalized using a reference sample which is included in each MS run. We show, however, that reference normalization does not effectively remove systematic MS run bias. A linear model approach was previously proposed to improve on the reference normalization approach but does not computationally scale to larger data sets. Here we describe the NOMAD (normalization of mass spectrometry data) R package which implements a computationally efficient ANOVA normalization approach with protein assembly functionality. NOMAD provides the same advantages as the linear regression solution but is more computationally efficient which allows superior scaling to larger sample sizes. Moreover, NOMAD effectively removes bias which improves valid across MS run comparisons.

An Accessible Proteogenomics Informatics Resource for Cancer Researchers

Proteogenomics has emerged as a valuable approach in cancer research, which integrates genomic and transcriptomic data with mass spectrometry-based proteomics data to directly identify expressed, variant protein sequences that may have functional roles in cancer. This approach is computationally intensive, requiring integration of disparate software tools into sophisticated workflows, challenging its adoption by nonexpert, bench scientists. To address this need, we have developed an extensible, Galaxy-based resource aimed at providing more researchers access to, and training in, proteogenomic informatics. Our resource brings together software from several leading research groups to address two foundational aspects of proteogenomics: (i) generation of customized, annotated protein sequence databases from RNA-Seq data; and (ii) accurate matching of tandem mass spectrometry data to putative variants, followed by filtering to confirm their novelty. Directions for accessing software tools and workflows, along with instructional documentation, can be found at z.umn.edu/canresgithub. Cancer Res; 77(21); e43-46. ©2017 AACR.

Increased Indoleamine-2,3-Dioxygenase Activity Is Associated With Poor Clinical Outcome in Adults Hospitalized With Influenza in the INSIGHT FLU003Plus Study

Background

Indoleamine-2,3-dioxygenase (IDO) mediated tryptophan (TRP) depletion has antimicrobial and immuno-regulatory effects. Increased kynurenine (KYN)-to-TRP (KT) ratios, reflecting increased IDO activity, have been associated with poorer outcomes from several infections.

Methods

We performed a case-control (1:2; age and sex matched) analysis of adults hospitalized with influenza A(H1N1)pdm09 with protocol-defined disease progression (died/transferred to ICU/mechanical ventilation) after enrollment (cases) or survived without progression (controls) over 60 days of follow-up. Conditional logistic regression was used to analyze the relationship between baseline KT ratio and other metabolites and disease progression.

Results

We included 32 cases and 64 controls with a median age of 52 years; 41% were female, and the median durations of influenza symptoms prior to hospitalization were 8 and 6 days for cases and controls, respectively (P = .04). Median baseline KT ratios were 2-fold higher in cases (0.24 mM/M; IQR, 0.13-0.40) than controls (0.12; IQR, 0.09-0.17; P ≤ .001). When divided into tertiles, 59% of cases vs 20% of controls had KT ratios in the highest tertile (0.21-0.84 mM/M). When adjusted for symptom duration, the odds ratio for disease progression for those in the highest vs lowest tertiles of KT ratio was 9.94 (95% CI, 2.25-43.90).

Conclusions

High KT ratio was associated with poor outcome in adults hospitalized with influenza A(H1N1)pdm09. The clinical utility of this biomarker in this setting merits further exploration.

ClinicalTrialsgov Identifier

NCT01056185.