Quantitative proteogenomics of human pathogens using DIA-MS

Proteome profiling of Campylobacter jejuni 81-176 at 37 °C and 42 °C by label-free mass spectrometry

BACKGROUND

The main natural reservoir for Campylobacter jejuni is the avian intestinal tract. There, C. jejuni multiplies optimally at 42 °C - the avian body temperature. After infecting humans through oral intake, the bacterium encounters the lower temperature of 37 °C in the human intestinal tract. Proteome profiling by label-free mass spectrometry (DIA-MS) was performed to examine the processes which enable C. jejuni 81-176 to thrive at 37 °C in comparison to 42 °C. In total, four states were compared with each other: incubation for 12 h at 37 °C, for 24 h at 37 °C, for 12 h at 42 °C and 24 h at 42 °C.

RESULTS

It was shown that the proteomic changes not only according to the different incubation temperature but also to the length of the incubation period were evident when comparing 37 °C and 42 °C as well as 12 h and 24 h of incubation. Altogether, the expression of 957 proteins was quantifiable. 37.1 - 47.3% of the proteins analyzed showed significant differential regulation, with at least a 1.5-fold change in either direction (i.e. log2 FC ≥ 0.585 or log2 FC ≤ -0.585) and an FDR-adjusted p-value of less than 0.05. The significantly differentially expressed proteins could be arranged in 4 different clusters and 16 functional categories.

CONCLUSIONS

The C. jejuni proteome at 42 °C is better adapted to high replication rates than that at 37 °C, which was in particular indicated by the up-regulation of proteins belonging to the functional categories "replication" (e.g. Obg, ParABS, and NapL), "DNA synthesis and repair factors" (e.g. DNA-polymerase III, DnaB, and DnaE), "lipid and carbohydrate biosynthesis" (e.g. capsular biosynthesis sugar kinase, PrsA, AccA, and AccP) and "vitamin synthesis, metabolism, cofactor biosynthesis" (e.g. MobB, BioA, and ThiE). The relative up-regulation of proteins with chaperone function (GroL, DnaK, ClpB, HslU, GroS, DnaJ, DnaJ-1, and NapD) at 37 °C in comparison to 42 °C after 12 h incubation indicates a temporary lower-temperature proteomic response. Additionally the up-regulation of factors for DNA uptake (ComEA and RecA) at 37 °C compared to 42 °C indicate a higher competence for the acquisition of extraneous DNA at human body temperature.

The Investigation of Protein Profile and Meat Quality in Bovine Longissimus thoracic Frozen under Different Temperatures by Data-Independent Acquisition (DIA) Strategy

The influence of freezing on the protein profile and quality traits in bovine Longissimus thoracic (LT) muscle was investigated by the data-independent acquisition (DIA) technique. Compared to fresh meat, a total of 262 proteins were identified as differential abundance proteins (DAPs) in four frozen groups (−12 °C, −18 °C, −38 °C, and −80 °C). According to the bioinformatics analysis, most of the DAPs in the significant Go terms and the KEGG pathway were structure proteins and enzymes. Proteome changes in the frozen bovine muscle at −12 °C and −18 °C were more significant than those at −38 °C and −80 °C. The result was consistent with the deterioration trend of the meat quality. The correlation analysis revealed that 17 proteins were correlated closely with the color, shear force, thawing loss, and cooking loss of the frozen meat, which could be used as putative biomarkers for frozen meat quality. MYO18A and ME3 are newly discovered proteins that are associated with frozen beef quality. In addition, CTTN and SERPINB6 were identified in frozen groups, which exhibited a significant inverse correlation with thawing loss (p < 0.01). These findings reveal the quality changes induced by freezing at the protein molecular level and provide new insights into the control of quality deterioration.

Identification of Protein Isoforms Using Reference Databases Built from Long and Short Read RNA-Sequencing

Alternative splicing can lead to distinct protein isoforms. These can have different functions in specific cells and tissues or in different developmental stages. In this study, we explored whether transcripts assembled from long read, nanopore-based, direct RNA-sequencing (RNA-seq) could improve the identification of protein isoforms in human K562 cells. By comparing with Illumina-based short read RNA-seq, we showed that a large proportion of Ensembl transcripts (5949/14,326) and genes expressing alternatively spliced transcripts (486/2981) identified with long direct reads were missed by short paired-end reads. By co-analyzing proteomic and transcriptomic data, we also showed that some peptides (826/35,976), proteins (262/3215), and protein isoforms arising from distinct transcript variants (574/1212) identified with isoform-specific peptides via custom long-read-based databases were missed in Illumina-derived databases. Finally, we generated unequivocal peptide evidence for a set of protein isoforms and showed that long read, direct RNA-seq allows the discovery of novel protein isoforms not already in reference databases or custom databases built from short read RNA-seq data. Our analysis highlights the benefits of long read RNA-seq data in the generation of reference databases to increase tandem mass spectrometry (MS/MS) identification of protein isoforms.

Proteogenomic Workflow Reveals Molecular Phenotypes Related to Breast Cancer Mammographic Appearance

Proteogenomic approaches have enabled the generat̲ion of novel information levels when compared to single omics studies although burdened by extensive experimental efforts. Here, we improved a data-independent acquisition mass spectrometry proteogenomic workflow to reveal distinct molecular features related to mammographic appearances in breast cancer. Our results reveal splicing processes detectable at the protein level and highlight quantitation and pathway complementarity between RNA and protein data. Furthermore, we confirm previously detected enrichments of molecular pathways associated with estrogen receptor-dependent activity and provide novel evidence of epithelial-to-mesenchymal activity in mammography-detected spiculated tumors. Several transcript-protein pairs displayed radically different abundances depending on the overall clinical properties of the tumor. These results demonstrate that there are differentially regulated protein networks in clinically relevant tumor subgroups, which in turn alter both cancer biology and the abundance of biomarker candidates and drug targets.

Cerebrospinal fluid proteome maps detect pathogen-specific host response patterns in meningitis

Meningitis is a potentially life-threatening infection characterized by the inflammation of the leptomeningeal membranes. Many different viral and bacterial pathogens can cause meningitis, with differences in mortality rates, risk of developing neurological sequelae, and treatment options. Here, we constructed a compendium of digital cerebrospinal fluid (CSF) proteome maps to define pathogen-specific host response patterns in meningitis. The results revealed a drastic and pathogen-type specific influx of tissue-, cell-, and plasma proteins in the CSF, where, in particular, a large increase of neutrophil-derived proteins in the CSF correlated with acute bacterial meningitis. Additionally, both acute bacterial and viral meningitis result in marked reduction of brain-enriched proteins. Generation of a multiprotein LASSO regression model resulted in an 18-protein panel of cell- and tissue-associated proteins capable of classifying acute bacterial meningitis and viral meningitis. The same protein panel also enabled classification of tick-borne encephalitis, a subgroup of viral meningitis, with high sensitivity and specificity. The work provides insights into pathogen-specific host response patterns in CSF from different disease etiologies to support future classification of pathogen type based on host response patterns in meningitis.

Automated Workflow for Peptide-Level Quantitation from DIA/SWATH-MS Data

Data-independent acquisition (DIA) is a powerful method to acquire spectra from all ionized precursors of a sample. Considering the complexity of the highly multiplexed spectral data, sophisticated workflows have been developed to obtain peptides quantification. Here we describe an open-source and easy-to-use workflow to obtain a quantitative matrix from multiple DIA runs. This workflow requires as prior information an "assay library," which contains the MS coordinates of peptides. It consists of OpenSWATH, pyProphet, and DIAlignR software. For the ease of installation and to isolate operating system-related dependency, docker-based containerization is utilized in this workflow.

A Proteogenomic Resource Enabling Integrated Analysis of Listeria Genotype-Proteotype-Phenotype Relationships

Listeria monocytogenes is an opportunistic foodborne pathogen responsible for listeriosis, a potentially fatal foodborne disease. Many different Listeria strains and serotypes exist, but a proteogenomic resource that bridges the gap in our molecular understanding of the relationships between the Listeria genotypes and phenotypes via proteotypes is still missing. Here, we devised a next-generation proteogenomics strategy that enables the community to rapidly proteotype Listeria strains and relate this information back to the genotype. Based on sequencing and de novo assembly of the two most commonly used Listeria model strains, EGD-e and ScottA, we established two comprehensive Listeria proteogenomic databases. A genome comparison established core- and strain-specific genes potentially responsible for virulence differences. Next, we established a DIA/SWATH-based proteotyping strategy, including a new and robust sample preparation workflow, that enables the reproducible, sensitive, and relative quantitative measurement of Listeria proteotypes. This reusable and publicly available DIA/SWATH library covers 70% of open reading frames of Listeria and represents the most extensive spectral library for Listeria proteotype analysis to date. We used these two new resources to investigate the Listeria proteotype in states mimicking the upper gastrointestinal passage. Exposure of Listeria to bile salts at 37 °C, which simulates conditions encountered in the duodenum, showed significant proteotype perturbations including an increase of FlaA, the structural protein of flagella. Given that Listeria is known to lose its flagella above 30 °C, this was an unexpected finding. The formation of flagella, which might have implications on infectivity, was validated by parallel reaction monitoring and light and scanning electron microscopy. flaA transcript levels did not change significantly upon exposure to bile salts at 37 °C, suggesting regulation at the post-transcriptional level. Together, these analyses provide a comprehensive proteogenomic resource and toolbox for the Listeria community enabling the analysis of Listeria genotype-proteotype-phenotype relationships.

The Core Proteome of Biofilm-Grown Clinical Pseudomonas aeruginosa Isolates

Comparative genomics has greatly facilitated the identification of shared as well as unique features among individual cells or tissues, and thus offers the potential to find disease markers. While proteomics is recognized for its potential to generate quantitative maps of protein expression, comparative proteomics in bacteria has been largely restricted to the comparison of single cell lines or mutant strains. In this study, we used a data independent acquisition (DIA) technique, which enables global protein quantification of large sample cohorts, to record the proteome profiles of overall 27 whole genome sequenced and transcriptionally profiled clinical isolates of the opportunistic pathogen Pseudomonas aeruginosa. Analysis of the proteome profiles across the 27 clinical isolates grown under planktonic and biofilm growth conditions led to the identification of a core biofilm-associated protein profile. Furthermore, we found that protein-to-mRNA ratios between different P. aeruginosa strains are well correlated, indicating conserved patterns of post-transcriptional regulation. Uncovering core regulatory pathways, which drive biofilm formation and associated antibiotic tolerance in bacterial pathogens, promise to give clues to interactions between bacterial species and their environment and could provide useful targets for new clinical interventions to combat biofilm-associated infections.

A quantitative Streptococcus pyogenes-human protein-protein interaction map reveals localization of opsonizing antibodies

A fundamental challenge in medical microbiology is to characterize the dynamic protein–protein interaction networks formed at the host–pathogen interface. Here, we generate a quantitative interaction map between the significant human pathogen, Streptococcus pyogenes, and proteins from human saliva and plasma obtained via complementary affinity-purification and bacterial-surface centered enrichment strategies and quantitative mass spectrometry. Perturbation of the network using immunoglobulin protease cleavage, mixtures of different concentrations of saliva and plasma, and different S. pyogenes serotypes and their isogenic mutants, reveals how changing microenvironments alter the interconnectivity of the interaction map. The importance of host immunoglobulins for the interaction with human complement proteins is demonstrated and potential protective epitopes of importance for phagocytosis of S. pyogenes cells are localized. The interaction map confirms several previously described protein–protein interactions; however, it also reveals a multitude of additional interactions, with possible implications for host–pathogen interactions involving other bacterial species.

Characterizing host-pathogen protein interactions can help elucidate the molecular basis of bacterial infections. Here, the authors use an integrative proteomics approach to generate a quantitative map of protein interactions between Streptococcus pyogenes and human saliva and plasma.

Computational Health Engineering Applied to Model Infectious Diseases and Antimicrobial Resistance Spread

Infectious diseases are the primary cause of mortality worldwide. The dangers of infectious disease are compounded with antimicrobial resistance, which remains the greatest concern for human health. Although novel approaches are under investigation, the World Health Organization predicts that by 2050, septicaemia caused by antimicrobial resistant bacteria could result in 10 million deaths per year. One of the main challenges in medical microbiology is to develop novel experimental approaches, which enable a better understanding of bacterial infections and antimicrobial resistance. After the introduction of whole genome sequencing, there was a great improvement in bacterial detection and identification, which also enabled the characterization of virulence factors and antimicrobial resistance genes. Today, the use of in silico experiments jointly with computational and machine learning offer an in depth understanding of systems biology, allowing us to use this knowledge for the prevention, prediction, and control of infectious disease. Herein, the aim of this review is to discuss the latest advances in human health engineering and their applicability in the control of infectious diseases. An in-depth knowledge of host–pathogen–protein interactions, combined with a better understanding of a host’s immune response and bacterial fitness, are key determinants for halting infectious diseases and antimicrobial resistance dissemination.

Greedy de novo motif discovery to construct motif repositories for bacterial proteomes

BACKGROUND

Bacterial surfaces are complex systems, constructed from membranes, peptidoglycan and, importantly, proteins. The proteins play crucial roles as critical regulators of how the bacterium interacts with and survive in its environment. A full catalog of the motifs in protein families and their relative conservation grade is a prerequisite to target the protein-protein interaction that bacterial surface protein makes to host proteins.

RESULTS

In this paper, we propose a greedy approach to identify conserved motifs in large sequence families iteratively. Each iteration discovers a motif de novo and masks all occurrences of that motif. Remaining unmasked sequences are subjected to the next round of motif detection until no more significant motifs can be found. We demonstrate the utility of the method through the construction of a proteome-wide motif repository for Group A Streptococcus (GAS), a significant human pathogen. GAS produce numerous surface proteins that interact with over 100 human plasma proteins, helping the bacteria to evade the host immune response. We used the repository to find that proteins part of the bacterial surface has motif architectures that differ from intracellular proteins.

CONCLUSIONS

We elucidate that the M protein, a coiled-coil homodimer that extends over 500 A from the cell wall, has a motif architecture that differs between various GAS strains. As the M protein is known to bind a variety of different plasma proteins, the results indicate that the different motif architectures are responsible for the quantitative differences of plasma proteins that various strains bind. The speed and applicability of the method enable its application to all major human pathogens.

An integrative strategy to identify the entire protein coding potential of prokaryotic genomes by proteogenomics

Accurate annotation of all protein-coding sequences (CDSs) is an essential prerequisite to fully exploit the rapidly growing repertoire of completely sequenced prokaryotic genomes. However, large discrepancies among the number of CDSs annotated by different resources, missed functional short open reading frames (sORFs), and overprediction of spurious ORFs represent serious limitations. Our strategy toward accurate and complete genome annotation consolidates CDSs from multiple reference annotation resources, ab initio gene prediction algorithms and in silico ORFs (a modified six-frame translation considering alternative start codons) in an integrated proteogenomics database (iPtgxDB) that covers the entire protein-coding potential of a prokaryotic genome. By extending the PeptideClassifier concept of unambiguous peptides for prokaryotes, close to 95% of the identifiable peptides imply one distinct protein, largely simplifying downstream analysis. Searching a comprehensive Bartonella henselae proteomics data set against such an iPtgxDB allowed us to unambiguously identify novel ORFs uniquely predicted by each resource, including lipoproteins, differentially expressed and membrane-localized proteins, novel start sites and wrongly annotated pseudogenes. Most novelties were confirmed by targeted, parallel reaction monitoring mass spectrometry, including unique ORFs and single amino acid variations (SAAVs) identified in a re-sequenced laboratory strain that are not present in its reference genome. We demonstrate the general applicability of our strategy for genomes with varying GC content and distinct taxonomic origin. We release iPtgxDBs for B. henselae, Bradyrhizobium diazoefficiens and Escherichia coli and the software to generate both proteogenomics search databases and integrated annotation files that can be viewed in a genome browser for any prokaryote.

Clinical applications of quantitative proteomics using targeted and untargeted data-independent acquisition techniques

INTRODUCTION

While selected/multiple-reaction monitoring (SRM or MRM) is considered the gold standard for quantitative protein measurement, emerging data-independent acquisition (DIA) using high-resolution scans have opened a new dimension of high-throughput, comprehensive quantitative proteomics. These newer methodologies are particularly well suited for discovery of biomarker candidates from human disease samples, and for investigating and understanding human disease pathways. Areas covered: This article reviews the current state of targeted and untargeted DIA mass spectrometry-based proteomic workflows, including SRM, parallel-reaction monitoring (PRM) and untargeted DIA (e.g., SWATH). Corresponding bioinformatics strategies, as well as application in biological and clinical studies are presented. Expert commentary: Nascent application of highly-multiplexed untargeted DIA, such as SWATH, for accurate protein quantification from clinically relevant and disease-related samples shows great potential to comprehensively investigate biomarker candidates and understand disease.

Targeted Proteomics and Absolute Protein Quantification for the Construction of a Stoichiometric Host-Pathogen Surface Density Model

Sepsis is a systemic immune response responsible for considerable morbidity and mortality. Molecular modeling of host-pathogen interactions in the disease state represents a promising strategy to define molecular events of importance for the transition from superficial to invasive infectious diseases. Here we used the Gram-positive bacterium Streptococcus pyogenes as a model system to establish a mass spectrometry based workflow for the construction of a stoichiometric surface density model between the S. pyogenes surface, the surface virulence factor M-protein, and adhered human blood plasma proteins. The workflow relies on stable isotope labeled reference peptides and selected reaction monitoring mass spectrometry analysis of a wild-type strain and an M-protein deficient mutant strain, to generate absolutely quantified protein stoichiometry ratios between S. pyogenes and interacting plasma proteins. The stoichiometry ratios in combination with a novel targeted mass spectrometry method to measure cell numbers enabled the construction of a stoichiometric surface density model using protein structures available from the protein data bank. The model outlines the topology and density of the host-pathogen protein interaction network on the S. pyogenes bacterial surface, revealing a dense and highly organized protein interaction network. Removal of the M-protein from S. pyogenes introduces a drastic change in the network topology, validated by electron microscopy. We propose that the stoichiometric surface density model of S. pyogenes in human blood plasma represents a scalable framework that can continuously be refined with the emergence of new results. Future integration of new results will improve the understanding of protein-protein interactions and their importance for bacterial virulence. Furthermore, we anticipate that the general properties of the developed workflow will facilitate the production of stoichiometric surface density models for other types of host-pathogen interactions.

Proteogenomic studies on cancer drug resistance: towards biomarker discovery and target identification

INTRODUCTION

Chemoresistance is a major obstacle for current cancer treatment. Proteogenomics is a powerful multi-omics research field that uses customized protein sequence databases generated by genomic and transcriptomic information to identify novel genes (e.g. noncoding, mutation and fusion genes) from mass spectrometry-based proteomic data. By identifying aberrations that are differentially expressed between tumor and normal pairs, this approach can also be applied to validate protein variants in cancer, which may reveal the response to drug treatment. Areas covered: In this review, we will present recent advances in proteogenomic investigations of cancer drug resistance with an emphasis on integrative proteogenomic pipelines and the biomarker discovery which contributes to achieving the goal of using precision/personalized medicine for cancer treatment. Expert commentary: The discovery and comprehensive understanding of potential biomarkers help identify the cohort of patients who may benefit from particular treatments, and will assist real-time clinical decision-making to maximize therapeutic efficacy and minimize adverse effects. With the development of MS-based proteomics and NGS-based sequencing, a growing number of proteogenomic tools are being developed specifically to investigate cancer drug resistance.

Automated SWATH Data Analysis Using Targeted Extraction of Ion Chromatograms

Targeted mass spectrometry comprises a set of methods able to quantify protein analytes in complex mixtures with high accuracy and sensitivity. These methods, e.g., Selected Reaction Monitoring (SRM) and SWATH MS, use specific mass spectrometric coordinates (assays) for reproducible detection and quantification of proteins. In this protocol, we describe how to analyze, in a targeted manner, data from a SWATH MS experiment aimed at monitoring thousands of proteins reproducibly over many samples. We present a standard SWATH MS analysis workflow, including manual data analysis for quality control (based on Skyline) as well as automated data analysis with appropriate control of error rates (based on the OpenSWATH workflow). We also discuss considerations to ensure maximal coverage, reproducibility, and quantitative accuracy.