Identification and Validation of Specific Markers of Bacillus anthracis Spores by Proteomics and Genomics Approaches

Bacillus anthracis is the causative bacteria of anthrax, an acute and often fatal disease in humans. The infectious agent, the spore, represents a real bioterrorism threat and its specific identification is crucial. However, because of the high genomic relatedness within the Bacillus cereus group, it is still a real challenge to identify B. anthracis spores confidently. Mass spectrometry-based tools represent a powerful approach to the efficient discovery and identification of such protein markers. Here we undertook comparative proteomics analyses of Bacillus anthracis, cereus and thuringiensis spores to identify proteoforms unique to B. anthracis. The marker discovery pipeline developed combined peptide- and protein-centric approaches using liquid chromatography coupled to tandem mass spectrometry experiments using a high resolution/high mass accuracy LTQ-Orbitrap instrument. By combining these data with those from complementary bioinformatics approaches, we were able to highlight a dozen novel proteins consistently observed across all the investigated B. anthracis spores while being absent in B. cereus/thuringiensis spores. To further demonstrate the relevance of these markers and their strict specificity to B. anthracis, the number of strains studied was extended to 55, by including closely related strains such as B. thuringiensis 9727, and above all the B. cereus biovar anthracis CI, CA strains that possess pXO1- and pXO2-like plasmids. Under these conditions, the combination of proteomics and genomics approaches confirms the pertinence of 11 markers. Genes encoding these 11 markers are located on the chromosome, which provides additional targets complementary to the commonly used plasmid-encoded markers. Last but not least, we also report the development of a targeted liquid chromatography coupled to tandem mass spectrometry method involving the selection reaction monitoring mode for the monitoring of the 4 most suitable protein markers. Within a proof-of-concept study, we demonstrate the value of this approach for the further high throughput and specific detection of B. anthracis spores within complex samples.Bacillus anthracis is a highly virulent bacterium, which is the etiologic agent of anthrax, an acute and often lethal disease of animals and humans (1). According to the Centers for Disease Control and Prevention, B. anthracis is classified as a category A agent, the highest rank of potential bioterrorism agents (http://www.bt.cdc.gov/agent/agentlist-category.asp). The infectious agent of anthrax, the spore, was used as a bioterrorism weapon in 2001 in the United States when mailed letters containing B. anthracis spores caused 22 cases of inhalational and/or cutaneous anthrax, five of which were lethal (2). These events have emphasized the need for rapid and accurate detection of B. anthracis spores.Bacillus anthracis is a member of the genus Bacillus, Gram-positive, rod-shaped bacteria characterized by the ability to form endospores under aerobic or facultative anaerobic conditions (3). The genus Bacillus is a widely heterogeneous group encompassing 268 validly described species to date (http://www.bacterio.net/b/bacillus.html, last accessed on August 9^th 2013). B. anthracis is part of the B. cereus group which consists of six distinct species: B. anthracis, B. cereus, B. thuringiensis, B. mycoides, B. pseudomycoides, and B. weihenstephanensis (4, 5). The latter three species are generally regarded as nonpathogenic whereas B. cereus and B. thuringiensis could be opportunistic or pathogenic to mammals or insects (5, 6). B. cereus is a ubiquitous species that lives in soil but is also found in foods of plant and animal origin, such as dairy products (7). Its occurrence has also been linked to food poisoning and it can cause diarrhea and vomiting (6, 8). B. thuringiensis is primarily an insect pathogen, also present in soil, and often used as a biopesticide (9).B. anthracis is highly monomorphic, that is, shows little genetic variation (10), and primarily exists in the environment as a highly stable, dormant spore in the soil (1). Specific identification of B. anthracis is challenging because of its high genetic similarity (sequence similarity >99%) with B. cereus and B. thuringiensis (5, 11). The fact that these closely related species are rather omnipresent in the environment further complicates identification of B. anthracis. The main difference among these three species is the presence in B. anthracis of the two virulence plasmids pXO1 and pXO2 (1), which are responsible for its pathogenicity. pXO1 encodes a tripartite toxin (protective antigen (PA), lethal factor (LF), and edema factor (EF)) which causes edema and death (1), whereas pXO2 encodes a poly-γ-d-glutamate capsule which protects the organism from phagocytosis (1). B. anthracis identification often relies on the detection of the genes encoded by these two plasmids via nucleic acid-based assays (12–14). Nevertheless, the occasionally observed loss of the pXO2 plasmid within environmental species may impair the robustness of detection (1). In addition, in recent years a series of findings has shown that the presence of pXO1 and pXO2 is not a unique feature of B. anthracis. Indeed, Hu et al. have demonstrated that ∼7% of B. cereus/B. thuringiensis species can have a pXO1-like plasmid and ∼1.5% a pXO2-like plasmid (15). This was particularly underlined for some virulent B. cereus strains (i.e. B. cereus strains G9241, B. cereus biovar anthracis strains CA and CI) (16–20).Because of these potential drawbacks, the use of chromosome-encoded genes would be preferable for the specific detection of B. anthracis. Such genes (rpoB, gyrA, gyrB, plcR, BA5345, and BA813) have been reported as potential markers (21–25), but concerns have also been raised about their ability to discriminate B. anthracis efficiently from closely related B. cereus strains (26). Ahmod et al. have recently pointed out, by in silico database analysis, that a specific sequence deletion (indel) occurs in the yeaC gene and exploited it for the specific identification of B. anthracis (27). Nevertheless, a few B. anthracis strains (e.g. B. anthracis A1055) do not have this specific deletion and so may lead to false-negative results (27).In the last few years, protein profiling by MS, essentially based on matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF MS), has emerged as an alternative (or a complement) to genotypic or phenotypic methods for the fast and efficient identification of microorganisms (28, 29). Such an approach is based on the reproducible acquisition of global bacterial protein fingerprints/patterns. The combination of MS-based protein patterns and chemometric/bioinformatic tools has been demonstrated to efficiently differentiate members of the B. cereus group from other Bacillus species (30). However, the specific discrimination of B. anthracis from the closely related B. cereus and B. thuringiensis remains difficult (30). This study of Lasch and coworkers, performed on vegetative cells, identified a few ribosomal and spore proteins as being responsible for this clustering (30). Closer inspection of the data revealed that B. anthracis identification was essentially based on one particular isoform of the small acid-soluble spore protein B (SASP-B)¹ (30–34), which is exclusively expressed in spores, as the samples were shown to contain residual spores. However, the specificity of SASP-B has recently been questioned as the published genomes of B. cereus biovar anthracis CI and B. thuringiensis BGSC 4CC1 strains have been shown to share the same SASP-B isoform as B. anthracis (35). Altogether these results underline that the quest for specific markers of B. anthracis needs to be pursued.Mass spectrometry also represents a powerful tool for the discovery and identification of protein markers (36, 37). In the case of B. anthracis, this approach has more commonly been used for the comprehensive characterization of given bacterial proteomes. For example, the proteome of vegetative cells with variable plasmid contents has been extensively studied (38–40), as the proteomes of mature spores (41, 42) and of germinating spores (43, 44). Only one recent study, based on a proteo-genomic approach, was initiated to identify protein markers of B. anthracis (45). In this work, potential markers were characterized but using a very limited number of B. cereus group strains (three B. cereus and two B. thuringiensis). Moreover, this study was done on vegetative cells, whereas the spore proteome is drastically different. To our knowledge, no study has characterized and validated relevant protein markers specific to B. anthracis spores, which constitute the dissemination form of B. anthracis and are often targeted by first-line immunodetection methods (46).Here we report comparative proteomics analyses of Bacillus anthracis/cereus/thuringiensis spores, undertaken to identify proteoforms unique to B. anthracis. Preliminary identification was performed on a limited set of Bacillus species both at the peptide (after enzymatic digestion) and protein levels by liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) using a high resolution/high mass accuracy LTQ-Orbitrap instrument. The pertinence of 11 markers was further demonstrated using proteomics and genomics approaches on a representative larger set of up to 55 different strains, including the closely related B. cereus biovar anthracis CI, CA, and B. thuringiensis 9727. Lastly, as a proof-of-concept study, we also report for four B. anthracis markers the implementation of a targeted LC-MS/MS method using selected reaction monitoring (SRM), based on the extension of a previous one focused on SASP-B (35). Preliminary results regarding method usefulness for the high throughput and accurate detection of B. anthracis spores in complex samples were also obtained and will be reported herein.