Abstract

Brucellosis is a prevalent zoonotic disease and is endemic in the Middle East, South America, and other areas of the world. In this study, complete inventories of putative functional ABC systems of five Brucella species have been compiled and compared. ABC systems of Brucella melitensis 16M, Brucella abortus 9-941, Brucella canis RM6/66, Brucella suis 1330, and Brucella ovis 63/290 were identified and aligned. High numbers of ABC systems, particularly nutrient importers, were found in all Brucella species. However, differences in the total numbers of ABC systems were identified (B. melitensis, 79; B. suis, 72; B. abortus 64; B. canis, 74; B. ovis, 59) as well as specific differences in the functional ABC systems of the Brucella species. Since B. ovis is not known to cause human brucellosis, functional ABC systems absent in the B. ovis genome may represent virulence factors in human brucellosis.

1. Introduction

Brucella species are the causative agents of brucellosis, the world’s most prevalent zoonotic disease, with high occurrences in endemic areas including the Middle East, Asia, Mexico, and the Mediterranean [1]. The bacteria are small nonmotile, Gram-negative, nonspore-forming coccobacilli that reside within the subphylum α-proteobacteria, which also includes nitrogen-fixing bacteria of the genus Nitrobacter, Rhizobium, Agrobacterium, and Rickettsia [2]. They are considered facultative intracellular pathogens.

There are six traditionally recognised Brucella species that have different host preferences: Brucella melitensis (which usually infects sheep and goats), Brucella abortus (cattle), Brucella suis (pigs), Brucella ovis (sheep), Brucella canis (dogs), and Brucella neotomae (desert wood rats). Furthermore, there are three newly identified Brucella species isolated from marine mammals: Brucella pinnipedialis (seals) [3], Brucella ceti (dolphins and porpoises) [3], and Brucella microti (voles) [4]. Although Brucella are primarily animal pathogens causing infectious abortions in females and orchitis in males [5], four of the nine species may infect humans (B. melitensis, B. abortus, B. suis, and occasionally B. canis, in order of disease severity) causing a range of flu-like symptoms including fever, sweats, malaise, and nausea [6]. Transmission to humans takes place via three recognised channels: (i) the consumption of infected animal products, (ii) direct contact with infected animal birth products, and (iii) the inhalation of aerosolised Brucella. Due to the nature of the human disease and the ability to be infectious via aerosol, Brucella species have been classified as category B threat agents by the US Centre for Disease Control and Prevention (CDC) [7].

Genome sequence analysis of B. melitensis 16M [8], B. suis 1330 [9], B. abortus 9-941 [10], B. canis RM6/66 (NCBI: NC_009504 and NC_009505, unpublished), and B. ovis 63/290 (NCBI: NC_010103 and NC_0010104, unpublished) has demonstrated the close relatedness of these organisms [11, 12]. The genomic DNA of each strain comprises two chromosomes of approximately 2.1 Mb and 1.2 Mb. DNA-DNA hybridisations between the species had previously revealed over similarity between the species, leading to the suggestion that all Brucella species should be classified as B. melitensis [13, 14]. However, it is widely believed that the differences in host specificity and pathogenicity are related to Brucella genetics; although there is currently little experimental evidence to support this, a few studies have found differences between the Brucella species genomes that may support this hypothesis [10, 15, 16]. A significant proportion of the Brucella genomes appear to code for ATP-binding cassette (ABC) systems.

ABC transporters are responsible for the import and export of many different substances across cellular membranes [17]. Although ABC transporters are extremely versatile, they all contain one defining feature, the ability to hydrolyse ATP to ADP, providing the energy needed for active transport. ABCs have three main conserved motifs known as Walker A (G-X-X-G-X-G-K-S/T, where X represents any amino acid residue), Walker B (ø-ø-ø-ø-D, where ø designates a hydrophobic residue), and a signature sequence (LSGGQ) [18]. The Walker A and Walker B motifs form tertiary structure enabling ATP-binding and can be found in all ATP-binding molecules. The signature sequence is well conserved in all ABC proteins and is also known as the linker peptide or C motif [19]. Although the configuration of ABC systems varies, the majority of ABC systems comprise of two hydrophilic ABC domains associated with two hydrophobic membrane-spanning domains (inner membrane (IM) proteins). Import systems are only found in prokaryotic organisms and contain both ABC domains and IM domains, along with extra-cytoplasmic binding proteins (BPs) designed to bind the specific allocrite of that ABC system. In Gram-negative bacteria the BPs are located in the periplasm whereas, in Gram-positive bacteria, they are anchored to the outer membrane of the cell via N-terminal lipid groups [20]. ABC systems import a diverse range of substrates into the bacterial cell including peptides [21], polyamines [22], metal ions [23], amino acids [24], iron [25], and sulphates [26]. In comparison, ABC systems involved in export functions usually contain only IM and ABC domains fused together via either the N-terminus (IM-ABC) or the C-terminus (ABC-IM), which homodimerise to create a functional system [27]. Substances exported by ABC transporters include antibiotics in both producing and resistant bacteria [28, 29], fatty acids in Gram-negative bacteria [27], and toxins [30]. In addition to transporters, many ABC proteins have roles in house-keeping functions, such as regulation of gene expression [31] and DNA repair [27, 32]. These proteins do not contain IM domains but are constituted of two fused ABC domains (ABC2) [27]. There is now increasing evidence that ABC systems can play roles in bacterial virulence [3336] and can be used as targets for vaccine development [37].

The recent sequencing of the genomes of B. melitensis 16M [8], B. abortus 9-941 [10], B. suis 1330 [9], B. ovis 63/290 (NCBI: NC_009504 and NC_009505, unpublished), and B. canis RM6/66 (NCBI: NC_010103 and NC_0010104, unpublished) has enabled the genomic comparison of different Brucella species. We report the creation and comparison of reannotated inventories of the functional ABC systems in Brucella. This improved annotation has assisted in understanding Brucella lifestyles and the identification of ABC systems that may be involved in virulence.

2. Methods

The prediction of ABC systems in sequenced bacterial genomes is based on annotation- and similarity-based homology assessment of identified or predicted ABC proteins from heterologous bacterial systems. The Artemis viewer (available from http://www.sanger.ac.uk) was used to visualise the sequenced genomes of B. melitensis 16 M, B. suis 1330, B. abortus 9-941, B. canis RM6/66, and B. ovis 63/290 [810]. Using the annotated genomes, ABC proteins were searched for using an array of related words, specifically “ATP-binding cassettes,” “binding protein”, or “outer membrane protein.” For completeness all proteins that were labelled as hypothetical or conserved hypothetical proteins were also checked. Hits from this search were compiled and then genes upstream and downstream were also checked to ensure that all genes from one system were found. After the genome searches were completed, protein sequences were aligned using the basic local alignment search tool (BlastP) against other ABC proteins using the ABC systems: Information on Sequence Structure and Evolution (ABCISSE) database [27, 38]. The ABCISSE database comprises 24000 proteins from 9500 annotated systems over 795 different organisms. Proteins searched against ABCISSE that scored a threshold -value of 10-6 were assigned to an ABC family and subfamily based on the hits from the ABCISSE database. Where searches on ABCISSE were unclear or hits for multiple families were produced, proteins were aligned using BlastP searches against the Genbank protein database. Use of this larger database increased the number of positive hits and functions that could be assigned. An ABC system was defined as a series of contiguous ORFs that shared the same family, subfamily, and substrate. A complete signal sequence (LSGGQ) was identified in the majority of the ABC proteins identified, and all of the other ABC proteins contained remnants of a complete signal sequence. Walker A and Walker B sequences were not sought during these searches.

The ABC system inventories compiled in this study include systems that contain genes with predicted frame shift mutations and premature stop codons. For example, the B. melitensis 16M gene BMEII0099 is a known pseudogene with multiple premature stop codons. However, this gene is part of an ABC system that is encoded by another four genes (BMEII0098, BMEII00101, BMEII102, and BMEII0103), all of which are predicted to be functional; the mutation in BMEII0099 might render the whole system nonfunctional or it is possible that the other four genes could create a partially functional system. Due to the inability to determine the functionality of ABC systems using bioinformatic techniques, the ABC systems where one or more components were predicted to be nonfunctional were excluded from the total ABC system numbers and functions of the ABC systems. Within the genomes of all Brucella species single components of ABC systems (mainly BP) not attached to individual systems were located. These were included in ABC system inventories and termed lone components but were not included in total functional ABC system counts.

3. Results and Discussion

The genome structures of Brucella species are very similar [10–12], and although it is widely believed that the differences in Brucella species virulence and host preferences are related to their genetic composition, there is little experimental evidence to support this belief. However, there are a few studies that have uncovered differences between the genomes [10, 15, 16]. In this study we have compared the presence of putative functional ABC systems in the genomes of B. melitensis 16M (BM), B. suis 1330 (BS), B. abortus 9-941(BS), B. canis RM6/66 (BC), and B. ovis 63/290 (BO). In the original annotations of these genomes, a uniform nomenclature was not used and functional assignment of the systems varied considerably. Here we describe a reannotation of the ABC systems of these bacterial strains, leading to new predicted functions of the systems and predictions about how the individual components combine to form functional systems. Complete inventories of the ABC systems of BM, BS, BA, BC, and BO are shown in Table 1.

The Brucella strains investigated in this study all have approximately 3.3 Mb genomes comprising two chromosomes of approximately 2.1 Mb and 1.2 Mb. The total number of predicted functional ABC systems encoded by the genomes of the Brucella strains is similar but does show some variability (BM = 79, BS = 72, BA = 64, BC = 74, BO = 59). Our evaluation of the Brucella genomes confirms that these species encode a relatively high proportion of ABC system genes when compared to other bacteria [39], with an average of 8.8% of their genomes dedicated to predicted functional ABC system genes (if lone components and mutated genes are included this figure increases to 9.3%). This may reflect their relatedness to environmental -proteobacteria such as Nitrobacter and Agrobacterium which also encode high numbers of ABC systems [39] that may assist in their survival in diverse conditions.

This work reports the first full inventories of ABC systems within five genome-sequenced Brucella strains. There are a number of specific ABC systems/genes that have previously been identified in the published literature. For example, Paulsen et al. describe two ABC systems that are present in B. suis and absent in B. melitensis. The first of these is an ABC importer encoded by BR0952 (IM), BR0953 (IM), and BR0955 (BP) [9]. Although this particular system is listed in the inventory, the ABC protein component of the system was not located in the BS genome and so this system was deemed incomplete and unlikely to be functional. The system was almost completely missing from the BM genome which is consistent with the findings of Paulsen et al. [9]. The second reported system is encoded by BRA0630, BRA0631, BRA0632, BRA0633, BRA0634, and BRA0635. However, when these genes were assessed using ABCISSE, only two of the five genes were predicted to be ABC transporter binding proteins (BRA0631 and BRA0632) and no other ABC components were located. Thus we deem this system also likely to be nonfunctional. Other genes that have been identified in the literature are BRA1080 (a dipeptide ABC transporter protein indentified in BS), BMEI1742 (a mitochondrial export ABC transporter identified in BM), and BRA0749-BRA0750 (involved in oligopeptide import) [10], all of which are present in our inventories.

4. ABC System Functions

In this study, we have classified the ABC systems of BM, BS, BA, BC, and BO into classes, families, and subfamilies according to the functional classification system described by Dassa and Bouige [27] (Table 2). The Brucella strains encode 8–12 class 1 systems, characterised by an ABC-IM domain fusion and comprising predicted export systems, and 5 class 2 systems, characterised by a duplicated fused ABC and with predicted functions in antibiotic resistance and house-keeping functions. However, we have found that most of the ABC systems of Brucella species belong to class 3 with roles predicted in import processes. The further classification of Brucella ABC systems into families and subfamilies shows that there are a high number of ABC systems of specific importer families, particularly the MOI (minerals and organic ions), MOS (monosaccharide), OPN (oligopeptides and nickel), OSP (oligosaccharides and polyols), and OTCN (osmoprotectants taurine cyanate and nitrate) families, all of which primarily function to acquire nutrients.

The predicted functionality of the ABC systems within the Brucella genomes is dominated by ABC systems involved in the import of nutrients (Figure 1), and although this is not uncommon amongst bacteria, it is probable that Brucella species utilise ABC transporters to provide most of the nutrients they require [8, 39]. In support of the findings of Paulsen et al. [9], the 2.1 Mb chromosome encodes a large proportion of the ABC systems involved in molecular export and cellular process whereas the ABC systems located on the smaller chromosome are largely biased toward nutrient acquisition, leading to the idea that this second chromosome is important in the acquisition and processing of nutrients in Brucella.

Since the ABC systems were identified by homology searches, it is possible to assign each ABC importer with a predicted substrate that it imports, providing an overview of the ABC system-based import ability of the Brucella species. Table 3 shows the range of predicted substrates imported via ABC transporters within the Brucella genomes. Overall, our results show that there is little difference in the import ability between strains of the four species of Brucella that are pathogenic to humans (BM, BS, BA, and BC). However, BO lacks the ability to import 8 of the 26 listed nutrients via ABC systems. In fact, all of the 29 pseudogenes that are present within the BO ABC system inventory occur within nutrient importers. The nutrients that BO appears to be unable to import using ABC systems include polyamines (specifically spermidine and putrescine), nickel, thiamine, glycine betaine, erythritol, xylose, and molybdenum. It is possible that the defective uptake of one or more of these substrates by B. ovis may contribute to its likely lack of virulence in humans. For example, polyamines have recently been associated with bacterial virulence and pathogenicity in human pathogens [40] and polyamine transporters have therefore been targeted as novel vaccine candidate targets for human pathogens [41, 42].

Two predicted erythritol transport systems have been reported that have yet to be confirmed by experimental data [8, 43]. Although the erythritol transporter identified in this study has also been identified by Crasta et al. [43], it should be noted that B. abortus S19 has this transport system inactivated by pseudogenes and yet it is still able to incorporate erythritol [43], indicating that this ABC system might not be wholly responsible for erythritol transport. Another study has demonstrated that B. ovis does not utilise erythritol as readily as other sugars [44].

In this study we have identified one ABC system in BM that we have categorised within a new ABC system family (currently labelled NEW1; See Table 1). This system includes BP and IM proteins related to those of the MOS family and ABC proteins that are different to those from the MOS family. We previously identified a similar ABC system in the genomes of Burkholderia pseudomallei and Burkholderia mallei strains [45]. Clearly, experimental data is required to define the function of this system.

5. Differences between Brucella Species

Although there is similarity between the ABC system inventories of the Brucella strains studied in this work, we have identified systems that are absent in one or several Brucella species (Table 4). The systems that are absent from species are not critical for bacterial survival but could contribute to differences in the lifestyles or virulence of the Brucella species. Our data shows that there are ABC systems absent from all of the Brucella strains studied. In particular, BO (5 systems), BC (4 systems), and BA (4 systems) lack systems that are present in BM and/or BS. The absence of the ISB (formally ABCX) system from BO and BC is an interesting observation since the ISB systems are soluble complexes involved in labile [Fe-S] biogenesis, which is important in resistance to oxidative stresses. This could indicate that B. ovis and B. canis reside in environments that are low in oxygen or high in oxygen reducatants, or that they lack enzymes that need labile [Fe-S] centres [46, 47]. Furthermore, this difference may be a factor contributing to the reduced virulence for humans of B. ovis and B. canis when compared to B. melitensis, B. suis, and B. abortus. The CDI system absent from B.ovis is comprised of two proteins, FtsE (ABC protein) and FtsX (IM protein) [48], and has been studied in E. coli and other bacteria including Bacillus subtilis [49] and Mycobacterium tuberculosis [50]. This CDI system is involved in cell division. E. coli mutants of ftsE show a reduced growth capacity [51]. The MKL system absent from BC may play a role in toluene tolerance, since Tn5 insertions within the ttgA2 gene coding for the MKL ABC protein in Pseudomonas putida elicited a toluene-sensitive phenotype [52].

6. Conclusions

In this study the ABC systems of B. melitensis strain 16 M, B. suis strain 1330, B. abortus 9-941, B. canis strain RM6/66, and B. ovis strain 63/290 have been reannotated using the ABCISSE database in order to provide a new and improved set of annotated Brucella ABC systems for the strains studied. The information obtained and the uniform annotation and classification of ABC systems in these closely related species has enabled a more detailed analysis of the roles of ABC systems in Brucella species, contributing to an improved understanding of Brucella lifestyle and pathogenicity. Previous analysis of the Brucella genomes has shown that there is over 90% genome similarity between the Brucella species [13, 14]. Similarly, the ABC system inventory compiled in this work reflects the close similarities of the Brucella species. However, despite the high genetic homology of Brucella, this work highlighted differences in the predicted numbers and functions of the ABC systems encoded by each Brucella species. It is widely accepted that the three species that may cause the most human brucellosis are B. melitensis, B. suis, and B. abortus (and occasionally B. canis). This study has shown that these four species of Brucella have a larger set of ABC systems encoded within their genomes than B. ovis, which is not known to cause human disease. Although it is difficult to ascertain the exact effect of the loss of these ABC systems on B. ovis, it is possible to hypothesise that, along with other genetic differences observed [15], they contribute to its overall reduced virulence in humans. It should also be noted there that four further Brucella strains have been genome sequenced since this work was completed: B. melitensis 63/9, B. abortus 2308, B. abortus S19, and B. suis Thomsen. Compiling ABC systems inventories of these strains may identify further differences between strains that may have biological relevance. Among the newly sequenced strains are B. suis Thomsen, a strain which is not known to cause disease in humans, and B. abortus S19, a vaccine strain. ABC system inventories of these strains would be of particular interest since they are considered less pathogenic than the wild-type strains and yet the reasons for this lack of pathogenicity are currently unknown. Overall, the identified differences observed in the ABC system inventory of the Brucella strains studied should contribute to a greater understanding of differences in the lifestyles of the Brucella species.