Proteomics Goes to Court: A Statistical Foundation for Forensic Toxin/Organism Identification Using Bottom-Up Proteomics.

Bottom-up proteomics is increasingly being used to characterize unknown environmental, clinical, and forensic samples. Proteomics-based bacterial identification typically proceeds by tabulating peptide "hits" (i.e., confidently identified peptides) associated with the organisms in a database; those organisms with enough hits are declared present in the sample. This approach has proven to be successful in laboratory studies; however, important research gaps remain. First, the common-practice reliance on unique peptides for identification is susceptible to a phenomenon known as signal erosion. Second, no general guidelines are available for determining how many hits are needed to make a confident identification. These gaps inhibit the transition of this approach to real-world forensic samples where conditions vary and large databases may be needed. In this work, we propose statistical criteria that overcome the problem of signal erosion and can be applied regardless of the sample quality or data analysis pipeline. These criteria are straightforward, producing a p-value on the result of an organism or toxin identification. We test the proposed criteria on 919 LC-MS/MS data sets originating from 2 toxins and 32 bacterial strains acquired using multiple data collection platforms. Results reveal a > 95% correct species-level identification rate, demonstrating the effectiveness and robustness of proteomics-based organism/toxin identification.

Production and Characterization of Desmalonichrome Relative Binding Affinity for Uranyl Ions in Relation to Other Siderophores.

Siderophores are iron (Fe)-binding secondary metabolites that have been investigated for their uranium-binding properties. Previous work has focused on characterizing hydroxamate types of siderophores, such as desferrioxamine B, for their uranyl (UO2)-binding affinity. Carboxylate forms of these metabolites hold potential to be more efficient chelators of UO2, yet they have not been widely studied. Desmalonichrome is a carboxylate siderophore that is not commercially available and so was obtained from the fungus Fusarium oxysporum cultivated under Fe-depleted conditions. The relative affinity for UO2 binding of desmalonichrome was investigated using a competitive analysis of binding affinities between UO2 acetate and different concentrations of Fe(III) chloride using electrospray ionization mass spectrometry. In addition to desmalonichrome, three other siderophores, including two hydroxamates (desferrioxamine B and desferrichrome) and one carboxylate (desferrichrome A), were studied to understand their relative affinities for the UO2(2+) ion at two pH values. The binding affinities of hydroxamate siderophores to UO2(2+) ions were observed to decrease with increasing Fe(III)Cl3 concentration at the lower pH. On the other hand, decreasing the pH has a smaller impact on the binding affinities between carboxylate siderophores and the UO2(2+) ion. Desmalonichrome in particular was shown to have the greatest relative affinity for UO2 at all pH and Fe(III) concentrations examined. These results suggest that acidic functional groups in the ligands are important for strong chelation with UO2 at lower pH.

Assessment of marker proteins identified in whole cell extracts for bacterial speciation using liquid chromatography electrospray ionization tandem mass spectrometry.

Staphylococcal strains (CoNS) were speciated in this study. Digests of proteins released from whole cells were converted to tryptic peptides for analysis. Liquid chromatography electrospray ionization tandem mass spectrometry (LC-ESI MS/MS, Orbitrap) was employed for peptide analysis. Data analysis was performed employing the open-source software X!Tandem which uses sequenced genomes to generate a virtual peptide database for comparison to experimental data. The search database was modified to include the genomes of the 11 Staphylococcus species most commonly isolated from man. The number of total peptides matching each protein along with the number of peptides specifically matching to the homologue (or homologues) for strains of the same species were assessed. Any peptides not matching to the species examined were considered conflict peptides. The proteins typically identified with the largest percentage of sequence coverage, number of matched peptides and number of peptides corresponding to only the correct species were elongation factor Tu (EF Tu) and enolase (Enol). Additional proteins with consistently observed peptides as well as peptides matching only homologues from the same species were citrate synthase (CS) and 1-pyrroline-5-carboxylate dehydrogenase (1P5CD). Protein markers, previously identified from gel slices, (aconitate hydratase and oxoglutarate dehydrogenase) were found to provide low confidence scores when employing whole cell digests. The methodological approach described here provides a simple yet elegant way of identification of staphylococci. However, perhaps more importantly the technology should be applicable universally for identification of any bacterial species.

Reprint of "Assessment of marker proteins identified in whole cell extracts for bacterial speciation using liquid chromatography electrospray ionization tandem mass spectrometry".

Staphylococcal strains (CoNS) were speciated in this study. Digests of proteins released from whole cells were converted to tryptic peptides for analysis. Liquid chromatography electrospray ionization tandem mass spectrometry (LC-ESI MS/MS, Orbitrap) was employed for peptide analysis. Data analysis was performed employing the open-source software X!Tandem which uses sequenced genomes to generate a virtual peptide database for comparison to experimental data. The search database was modified to include the genomes of the 11 Staphylococcus species most commonly isolated from man. The number of total peptides matching each protein along with the number of peptides specifically matching to the homologue (or homologues) for strains of the same species were assessed. Any peptides not matching to the species examined were considered conflict peptides. The proteins typically identified with the largest percentage of sequence coverage, number of matched peptides and number of peptides corresponding to only the correct species were elongation factor Tu (EF Tu) and enolase (Enol). Additional proteins with consistently observed peptides as well as peptides matching only homologues from the same species were citrate synthase (CS) and 1-pyrroline-5-carboxylate dehydrogenase (1P5CD). Protein markers, previously identified from gel slices, (aconitate hydratase and oxoglutarate dehydrogenase) were found to provide low confidence scores when employing whole cell digests. The methodological approach described here provides a simple yet elegant way of identification of staphylococci. However, perhaps more importantly the technology should be applicable universally for identification of any bacterial species.

Identification of staphylococcal species based on variations in protein sequences (mass spectrometry) and DNA sequence (sodA microarray).

This report is among the first using sequence variation in newly discovered protein markers for staphylococcal (or indeed any other bacterial) speciation. Variation, at the DNA sequence level, in the sodA gene (commonly used for staphylococcal speciation) provided excellent correlation. Relatedness among strains was also assessed using protein profiling using microcapillary electrophoresis and pulsed field electrophoresis. A total of 64 strains were analyzed including reference strains representing the 11 staphylococcal species most commonly isolated from man (Staphylococcus aureus and 10 coagulase negative species [CoNS]). Matrix assisted time of flight ionization/ionization mass spectrometry (MALDI TOF MS) and liquid chromatography-electrospray ionization tandem mass spectrometry (LC ESI MS/MS) were used for peptide analysis of proteins isolated from gel bands. Comparison of experimental spectra of unknowns versus spectra of peptides derived from reference strains allowed bacterial identification after MALDI TOF MS analysis. After LC-MS/MS analysis of gel bands bacterial speciation was performed by comparing experimental spectra versus virtual spectra using the software X!Tandem. Finally LC-MS/MS was performed on whole proteomes and data analysis also employing X!tandem. Aconitate hydratase and oxoglutarate dehydrogenase served as marker proteins on focused analysis after gel separation. Alternatively on full proteomics analysis elongation factor Tu generally provided the highest confidence in staphylococcal speciation.

Reprint of "Identification of staphylococcal species based on variations in protein sequences (mass spectrometry) and DNA sequence (sodA microarray)".

This report is among the first using sequence variation in newly discovered protein markers for staphylococcal (or indeed any other bacterial) speciation. Variation, at the DNA sequence level, in the sodA gene (commonly used for staphylococcal speciation) provided excellent correlation. Relatedness among strains was also assessed using protein profiling using microcapillary electrophoresis and pulsed field electrophoresis. A total of 64 strains were analyzed including reference strains representing the 11 staphylococcal species most commonly isolated from man (Staphylococcus aureus and 10 coagulase negative species [CoNS]). Matrix assisted time of flight ionization/ionization mass spectrometry (MALDI TOF MS) and liquid chromatography-electrospray ionization tandem mass spectrometry (LC ESI MS/MS) were used for peptide analysis of proteins isolated from gel bands. Comparison of experimental spectra of unknowns versus spectra of peptides derived from reference strains allowed bacterial identification after MALDI TOF MS analysis. After LC-MS/MS analysis of gel bands bacterial speciation was performed by comparing experimental spectra versus virtual spectra using the software X!Tandem. Finally LC-MS/MS was performed on whole proteomes and data analysis also employing X!tandem. Aconitate hydratase and oxoglutarate dehydrogenase served as marker proteins on focused analysis after gel separation. Alternatively on full proteomics analysis elongation factor Tu generally provided the highest confidence in staphylococcal speciation.

Protein markers for identification of Yersinia pestis and their variation related to culture.

The detection of high consequence pathogens, such as Yersinia pestis, is well established in biodefense laboratories for bioterror situations. Laboratory protocols are well established using specified culture media and a growth temperature of 37 °C for expression of specific antigens. Direct detection of Y. pestis protein markers, without prior culture, depends on their expression. Unfortunately protein expression can be impacted by the culture medium which cannot be predicted ahead of time. Furthermore, higher biomass yields are obtained at the optimal growth temperature (i.e. 28 °C-30 °C) and therefore are more likely to be used for bulk production. Analysis of Y. pestis grown on several types of media at 30 °C showed that several protein markers were found to be differentially detected in different media. Analysis of the identified proteins against a comprehensive database provided an additional level of organism identification. Peptides corresponding to variable regions of some proteins could separate large groups of strains and aid in organism identification. This work illustrates the need to understand variability of protein expression for detection targets. The potential for relating expression changes of known proteins to specific media factors, even in nutrient rich and chemically complex culture medium, may provide the opportunity to draw forensic information from protein profiles.

Characterization of residual medium peptides from Yersinia pestis cultures.

Here we demonstrate that when Yersinia pesitis is grown in laboratory media, peptides from the medium remain associated with cellular biomass even after washing and inactivation of the bacteria by different methods. These peptides are characteristic of the type of growth medium and of the manufacturer of the medium, reflecting the specific composition of the medium. We analyzed biomass-associated peptides from cultures of two attenuated strains of Yersinia pestis [KIM D27 (pgm-) and KIM D1 (lcr-)] grown in several formulations of 4 different media (tryptic soy broth (TSB), brain-heart infusion (BHI), Luria-Bertani broth (LB), and glucose (G) medium) made from components purchased from different suppliers. Despite the range of growth medium sources and the associated manufacturing processes used in their production, a high degree of peptide similarity was observed for a given medium recipe; however, notable differences in the termination points of select peptides were observed in media formulated using products from some suppliers, presumably reflecting the process by which a manufacturer performed protein hydrolysis for use in culture media. These results may help explain the presence of peptides not explicitly associated with target organisms during proteomic analysis of microbes and other biological systems that require culturing. While the primary aim of this work is to outline the range and type of medium peptides associated with Yersinia pestis biomass and improve the quality of proteomic measurements, these peptides may also represent a potentially useful forensic signature that could provide information about microbial culturing conditions.

Forensic proteomics of poxvirus production.

The field of microbial forensics has recently sought to develop methods to discern biological signatures to indicate production methods for biological agents. Viral agents have received less attention to date. Their obligate propagation in living cells makes purification from cellular material a challenge. This leads to potential carryover of protein-rich signatures of their production system. Here we have explored a proteomic analysis of vaccinia virus as a model poxvirus system in which to compare samples of virus propagated in different cell lines and subjected to different purification schemes. The proteomic data sets indicated viral, host cell and culture medium proteins. Several layers of data analysis were applied to build confidence in the peptide identification and capture information on the taxonomic utility of each. The analysis showed clear shifts in protein profiles with virus purification, with successive gradient purification steps showing different levels of viral protein enrichment. Peptides from cellular proteins, including those present in purified virus preparations, provided signatures which enabled discrimination of cell line substrates, including distinguishing between cells derived from different primate species. The ability to discern multiple aspects of viral production demonstrates the potential value of proteomic analysis as tool for microbial forensics.

Proteomic analysis of bronchoalveolar lavage fluid proteins from mice infected with Francisella tularensis ssp. novicida.

Francisella tularensis causes the zoonosis tularemia in humans and is one of the most virulent bacterial pathogens. We utilized a global proteomic approach to characterize protein changes in bronchoalveolar lavage fluid from mice exposed to one of three organisms, F. tularensis ssp. novicida, an avirulent mutant of F. tularensis ssp. novicida (F.t. novicida-ΔmglA), and Pseudomonas aeruginosa. The composition of bronchoalveolar lavage fluid (BALF) proteins was altered following infection, including proteins involved in neutrophil activation, oxidative stress, and inflammatory responses. Components of the innate immune response were induced including the acute phase response and the complement system; however, the timing of their induction varied. F. tularensis ssp. novicida infected mice do not appear to have an effective innate immune response in the first hours of infection; however, within 24 h, they show an upregulation of innate immune response proteins. This delayed response is in contrast to P. aeruginosa infected animals which show an early innate immune response. Likewise, F.t. novicida-ΔmglA infection initiates an early innate immune response; however, this response is diminished by 24 h. Finally, this study identifies several candidate biomarkers, including Chitinase 3-like-1 (CHI3L1 or YKL-40) and peroxiredoxin 1, that are associated with F. tularensis ssp. novicida but not P. aeruginosa infection.

Integration of gas chromatography mass spectrometry methods for differentiating ricin preparation methods.

The investigation of crimes involving chemical or biological agents is infrequent, but presents unique analytical challenges. The protein toxin ricin is encountered more frequently than other agents and is found in the seeds of Ricinus communis, commonly known as the castor plant. Typically, the toxin is extracted from castor seeds utilizing a variety of different recipes that result in varying purity of the toxin. Moreover, these various purification steps can also leave or differentially remove a variety of exogenous and endogenous residual components with the toxin that may indicate the type and number of purification steps involved. We have applied three gas chromatography-mass spectrometry (GC-MS) based analytical methods to measure the variation in seed carbohydrates and castor oil ricinoleic acid, as well as the presence of solvents used for purification. These methods were applied to the same samples prepared using four previously identified toxin preparation methods, starting from four varieties of castor seeds. The individual data sets for seed carbohydrate profiles, ricinoleic acid, or acetone amount each provided information capable of differentiating different types of toxin preparations across seed types. However, the integration of the data sets using multivariate factor analysis provided a clear distinction of all samples based on the preparation method, independent of the seed source. In particular, the abundance of mannose, arabinose, fucose, ricinoleic acid, and acetone were shown to be important differentiating factors. These complementary tools provide a more confident determination of the method of toxin preparation than would be possible using a single analytical method.

Proteomic Signatures of Antimicrobial Resistance in Yersinia pestis and Francisella tularensis.

Antimicrobial resistance (AMR) is a well-recognized, widespread, and growing issue of concern. With increasing incidence of AMR, the ability to respond quickly to infection with or exposure to an AMR pathogen is critical. Approaches that could accurately and more quickly identify whether a pathogen is AMR also are needed to more rapidly respond to existing and emerging biological threats. We examined proteins associated with paired AMR and antimicrobial susceptible (AMS) strains of Yersinia pestis and Francisella tularensis, causative agents of the diseases plague and tularemia, respectively, to identify whether potential existed to use proteins as signatures of AMR. We found that protein expression was significantly impacted by AMR status. Antimicrobial resistance-conferring proteins were expressed even in the absence of antibiotics in growth media, and the abundance of 10-20% of cellular proteins beyond those that directly confer AMR also were significantly changed in both Y. pestis and F. tularensis. Most strikingly, the abundance of proteins involved in specific metabolic pathways and biological functions was altered in all AMR strains examined, independent of species, resistance mechanism, and affected cellular antimicrobial target. We have identified features that distinguish between AMR and AMS strains, including a subset of features shared across species with different resistance mechanisms, which suggest shared biological signatures of resistance. These features could form the basis of novel approaches to identify AMR phenotypes in unknown strains.

Analysis of carbohydrate and fatty acid marker abundance in ricin toxin preparations for forensic information.

One challenge in the forensic analysis of ricin samples is determining the method and extent of sample preparation. Ricin purification from the source castor seeds is essentially a protein purification through removal of the nonprotein fractions of the seed. Two major, nonprotein constituents in the seed are the castor oil and carbohydrates. We used derivatization of carbohydrate and fatty acid markers followed by identification and quantification using gas chromatography/mass spectrometry (GC/MS) to assess compositional changes in ricin samples purified by different methods. The loss of ricinoleic acid indicated steps for oil removal had occurred, and a large decrease of ricinoleic acid was observed between unextracted mash and solvent extracted and protein precipitate preparations. Changes to the carbohydrate content of the sample were also observed following protein precipitation. The differential loss of arabinose relative to mannose was observed indicating the removal of the major carbohydrate fraction of the seed and enrichment of the protein content. When the data is combined and multivariate principle component analysis is applied, these changes in fatty acid and carbohydrate abundance are discriminating enough to be indicative of the preparation method used for each sample.

Residual agar determination in bacterial spores by electrospray ionization mass spectrometry.

Presented here is an analytical method to detect residual agar from a bacterial spore sample as an indication of culturing on an agar plate. This method is based on the resolubilization of agar polysaccharide from a bacterial spore sample, enzymatic digestion, followed by electrospray ionization tandem mass spectrometry (ESI-MS(n)) analysis for detection of a specific agar fragment ion. A range of Bacillus species and strains were selected to demonstrate the effectiveness of this approach. The characteristic agar fragment ion was detected in the spores grown on agar that were washed from 1 to 5 times, irradiated or nonirradiated, and not in the spores grown in broth. A sample containing approximately 10(8) spores is currently needed for confident detection of residual agar from culture on agar plates in the presence of bacterial spores with a limit of detection of approximately 1 ppm agar spiked into a broth-grown spore sample. The results of a proficiency test with 42 blinded samples are presented demonstrating the utility of this method with no false positives and only three false negatives for samples that were below the detection level of the method as documented.

Applications and challenges of forensic proteomics.

Mass spectrometry-based proteomics has been a useful tool for addressing numerous questions in basic biology research for many years. This success, combined with the maturity of mass spectrometric instrumentation, the ever-increasing availability of protein sequence databases derived from genome sequencing, and the growing sophistication of data analysis methods, places proteomics in a position to have an important role in biological forensics. Because proteins contain information about genotype (sequence) and phenotype (expression levels), proteomics methods can both identify biological samples and characterize the conditions that produced them. In addition to serving as a valuable orthogonal method to genomic analyses, proteomics can be used in cases where nucleic acids are absent, degraded, or uninformative. Mass spectrometry provides both broad applicability and exquisite specificity, often without customized detection reagents like primers or antibodies. This review briefly introduces proteomics methods, and surveys a variety of forensic applications (including criminal justice, historical, archaeological, and national security areas). Finally, challenges and crucial areas for further research are addressed.

Bayesian-integrated microbial forensics.

In the aftermath of the 2001 anthrax letters, researchers have been exploring ways to predict the production environment of unknown-source microorganisms. Culture medium, presence of agar, culturing temperature, and drying method are just some of the broad spectrum of characteristics an investigator might like to infer. The effects of many of these factors on microorganisms are not well understood, but the complex way in which microbes interact with their environments suggests that numerous analytical techniques measuring different properties will eventually be needed for complete characterization. In this work, we present a Bayesian statistical framework for integrating disparate analytical measurements. We illustrate its application to the problem of characterizing the culture medium of Bacillus spores using three different mass spectral techniques. The results of our study suggest that integrating data in this way significantly improves the accuracy and robustness of the analyses.

Detection of agar, by analysis of sugar markers, associated with Bacillus anthracis spores, after culture.

Detection of small quantities of agar associated with spores of Bacillus anthracis could provide key information regarding its source or growth characteristics. Agar, widely used in growth of bacteria on solid surfaces, consists primarily of repeating polysaccharide units of 3,6-anhydro-l-galactose (AGal) and galactose (Gal) with sulfated and O-methylated galactoses present as minor constituents. Two variants of the alditol acetate procedure were evaluated for detection of potential agar markers associated with spores. The first method employed a reductive hydrolysis step, to stabilize labile anhydrogalactose, by converting to anhydrogalactitol. The second eliminated the reductive hydrolysis step simplifying the procedure. Anhydrogalactitol, derived from agar, was detected using both derivatization methods followed by gas chromatography-mass spectrometry (GC-MS) analysis. However, challenges with artifactual background (reductive hydrolysis) or marker destruction (hydrolysis) respectively lead to the use of an alternative agar marker. A minor agar component, 6-O-methyl galactose (6-O-M gal), was readily detected in agar-grown but not broth-grown bacteria. Detection was optimized by the use of gas chromatography-tandem mass spectrometry (GC-MS-MS). With appropriate choice of sugar marker and analytical procedure, detection of sugar markers for agar has considerable potential in microbial forensics.

Laboratory strains of Bacillus anthracis exhibit pervasive alteration in expression of proteins related to sporulation under laboratory conditions relative to genetically related wild strains.

The spore forming pathogen Bacillus anthracis is the etiologic agent of anthrax in humans and animals. It cycles through infected hosts as vegetative cells and is eventually introduced into the environment where it generates an endospore resistant to many harsh conditions. The endospores are subsequently taken up by another host to begin the next cycle. Outbreaks of anthrax occur regularly worldwide in wildlife and livestock, and the potential for human infection exists whenever humans encounter infected animals. It is also possible to encounter intentional releases of anthrax spores, as was the case in October 2001. Consequently, it is important to be able to rapidly establish the provenance of infectious strains of B. anthracis. Here, we compare protein expression in seven low-passage wild isolates and four laboratory strains of B. anthracis grown under identical conditions using LC-MS/MS proteomic analysis. Of the 1,023 total identified proteins, 96 had significant abundance differences between wild and laboratory strains. Of those, 28 proteins directly related to sporulation were upregulated in wild isolates, with expression driven by Spo0A, CodY, and AbrB/ScoC. In addition, we observed evidence of changes in cell division and fatty acid biosynthesis between the two classes of strains, despite being grown under identical experimental conditions. These results suggest wild B. anthracis cells are more highly tuned to sporulate than their laboratory cousins, and this difference should be exploited as a method to differentiate between laboratory and low passage wild strains isolated during an anthrax outbreak. This knowledge should distinguish between intentional releases and exposure to strains in nature, providing a basis for the type of response by public health officials and investigators.

A ricin forensic profiling approach based on a complex set of biomarkers.

A forensic method for the retrospective determination of preparation methods used for illicit ricin toxin production was developed. The method was based on a complex set of biomarkers, including carbohydrates, fatty acids, seed storage proteins, in combination with data on ricin and Ricinus communis agglutinin. The analyses were performed on samples prepared from four castor bean plant (R. communis) cultivars by four different sample preparation methods (PM1-PM4) ranging from simple disintegration of the castor beans to multi-step preparation methods including different protein precipitation methods. Comprehensive analytical data was collected by use of a range of analytical methods and robust orthogonal partial least squares-discriminant analysis- models (OPLS-DA) were constructed based on the calibration set. By the use of a decision tree and two OPLS-DA models, the sample preparation methods of test set samples were determined. The model statistics of the two models were good and a 100% rate of correct predictions of the test set was achieved.

Method development for comprehensive extraction and analysis of marine toxins: Liquid-liquid extraction and tandem liquid chromatography separations coupled to electrospray tandem mass spectrometry.

A variety of toxins are produced by marine and freshwater microorganisms that present a threat to human health. These toxins have diverse chemical properties and specifically, a range of hydrophobicity. Methods for extraction and identification of these toxins are often geared toward specific classes of toxin depending on the sample type. There is a need for a general method of toxin extraction and identification for screening samples where the likely toxin content is not known a priori. We have applied a general method for metabolite extraction to toxin containing samples. This method was coupled with a simple dual liquid chromatography approach for separating a broad range of toxins. This liquid chromatography approach was coupled to triple quadrupole and quadrupole time-of-flight MS/MS platforms. The method was testing on a fish matrix for recovery of palytoxin as well as marine corals for detection of natural mixtures of palytoxin analogues. The recovery of palytoxin was found to produce a linear response (R2 of 0.95) when spiked into the fish matrix with a limit of quantitation of 2.5 ng/µL and recovery efficiency of 73% + /- 9%. The screening of corals revealed varying amount of palytoxin, and in one case, different palytoxin structural analogues. This demonstration illustrates the potential utility of this method for toxin extraction and detection.