Bordetella petrii is the only environmental species hitherto found among the otherwise host-restricted and pathogenic members of the genus Bordetella. Phylogenetically, it connects the pathogenic Bordetellae and environmental bacteria of the genera Achromobacter and Alcaligenes, which are opportunistic pathogens. B. petrii strains have been isolated from very different environmental niches, including river sediment, polluted soil, marine sponges and a grass root. Recently, clinical isolates associated with bone degenerative disease or cystic fibrosis have also been described.
In this manuscript we present the results of the analysis of the completely annotated genome sequence of the B. petrii strain DSMZ12804. B. petrii has a mosaic genome of 5,287,950 bp harboring numerous mobile genetic elements, including seven large genomic islands. Four of them are highly related to the clc element of Pseudomonas knackmussii B13, which encodes genes involved in the degradation of aromatics. Though being an environmental isolate, the sequenced B. petrii strain also encodes proteins related to virulence factors of the pathogenic Bordetellae, including the filamentous hemagglutinin, which is a major colonization factor of B. pertussis, and the master virulence regulator BvgAS. However, it lacks all known toxins of the pathogenic Bordetellae.
The genomic analysis suggests that B. petrii represents an evolutionary link between free-living environmental bacteria and the host-restricted obligate pathogenic Bordetellae. Its remarkable metabolic versatility may enable B. petrii to thrive in very different ecological niches.
The genus Bordetella comprises several human and animal pathogens and nine species are currently described [see Additional file 1]. Human-restricted B. pertussis and B. parapertussis cause whooping cough, whereas B. bronchiseptica and B. avium are responsible for respiratory infections in many mammals and birds, respectively [1–3]. Despite the availability of vaccines, there are still 300,000 deaths/year caused by B. pertussis and significant economic losses associated with infections in poultry and cattle. The genomes of these Bordetellae were recently sequenced and analyzed [4, 5]. The genome sequence of B. avium revealed that it is quite divergent from the mammalian pathogens. In contrast, the genomes of the human pathogens B. pertussis and B. parapertussis showed that they are independent derivatives of B. bronchiseptica-like ancestors. They both underwent significant gene loss, probably mediated by insertion sequence (IS) elements, during the process of host adaptation, as indicated by their smaller genome sizes compared to that of B. bronchiseptica and the presence of many pseudogenes. The evolution of human-adapted species is characterized by massive genome reduction, lack of horizontal acquisition of genetic material and significant differences in virulence gene expression among species [4, 6].
The type strain of the only environmental species of the genus, B. petrii, was isolated from a dechlorinating bioreactor enriched by river sediment . B. petrii strains were repeatedly found in environmental samples, e.g. microbial consortia degrading aromatic compounds [8, 9]. Recently, two B. petrii strains were isolated from patients suffering from mandibular osteomyelitis and chronic suppurative mastoiditis [10, 11]. Bacteria closely related to B. petrii were also found in cystic fibrosis patients . Moreover, B. petrii strains were also isolated from marine sponges and a grass root consortium [13, 14]. Phylogenetic analysis shows that the Bordetellae are closely related to Achromobacter and Alcaligenes [see Additional file 1] [2, 7]. These environmental organisms occasionally cause severe nosocomial infections and they are also isolated from patients with cystic fibrosis [15, 16]. In addition to their medical relevance, these microorganisms also have an intrinsic biotechnological interest, because they show a remarkable capacity to degrade aromatic compounds and detoxify heavy metals [17, 18].
Since B. petrii exhibits properties of host-associated and environmental bacteria, it is important to perform a genomic analysis of this bacterium in order to understand the genetic basis of its versatility. The B. petrii genome provides insights into the evolutionary history of the Bordetellae, acting also as a spyglass to the Alcaligenes/Achromobacter group for which, despite their biotechnological potential, no sequenced genomes are available yet.
Results and discussion
Genome features and comparative analysis
The general features of the Bordetella genomes are shown in Table 1 and Figure 1. The genome of B. petrii has a similar size (5,287,950 bp) and coding-capacity (5,034 CDSs) to that of B. bronchiseptica (5,339,179 bp; 5,019 CDSs). The genomes of B. pertussis, B. parapertussis and B. avium are much smaller, showing massive gene loss during host-adaptation. The overall GC content (65.5%) of B. petrii is intermediate between that of B. bronchiseptica (68.08%) and B. avium (61.58%). B. petrii has 125 pseudogenes, less than B. pertussis (358) and B. parapertussis (220), but more than B. bronchiseptica (18) and B. avium (68).
General features of the Bordetella genomes
GC content (%)
Average gene size (bp)
IS3-family and others
The genomes of B. petrii and B. avium show lower similarity and synteny than those of mammalian-adapted Bordetellae, suggesting that they are most distantly related. Reciprocal BLASTP analysis with other sequenced Bordetellae identified coding sequences unique to B. petrii (Figure 2) [4, 5]. B. petrii, B. bronchiseptica and B. avium share 2,049 CDSs, which represent the Bordetella core gene set, which was likely inherited from a common ancestor. B. petrii shares more CDSs with B. bronchiseptica (2,884) than with B. avium (2,229). Interestingly, B. petrii has 1,825 CDSs that are not present in either B. bronchiseptica or B. avium. A large number of these B. petrii-unique CDSs (1,157 CDSs and 188 genes encoding transposases) are located on mobile elements and encode accessory metabolic functions (Table 2).
Major features of the B. petrii genomic islands
Genomic islands (CDS)
Coordinates (size in bp)
No integrase and direct repeats (DR)
Capsular polysaccharide genes
Metabolism of phthalate and protocatechuate, urea amidohydrolase
clc-like, integrase, DR (3' end of tRNAGly)
Metabolism of an unknown aromatic compound (Bpet1116-1123)
clc-like, integrase, DR (3' end of tRNAGly)
Metabolism of benzoate, benzylalcohol, 3-hydroxybenzoate, putative monooxygenase (Bpet1330)
Putative (chloro)phenol monooxygenase (Bpet1331)
nitrile hydratase (Bpet1415-1416)
Almost identical to the clc element, integrase, DR (3' end of tRNAGly)
Metabolism of chlorocatechol, anthranilate
Short chain dehydrogenase (Bpet1512)
Tn4371-like, integrase, (T)6 DR
Integrase, features of conjugative transposons
Metabolism of chlorocatechol (2 gene clusters), tetrachlorobenzene
4 glutathione S-transferase (Bpet3724-3727)
clc-like, integrase, DR (3' end of tRNAGly)
Multidrug efflux pump
Iron transport system.
Integrase, DR (3' end of tRNAPhe), features of conjugative transposons
Several heavy metal resistance systems
Mobile genetic elements - a mosaic genome
The most prominent feature that distinguishes B. petrii from other Bordetellae is the presence of seven laterally acquired genomic islands (GI1-GI7). The locations of these GIs on the chromosome correlate with atypical GC-content and GC-bias (Figure 1). These GIs have the characteristic features of mobile elements collectively known as integrative and conjugative elements (ICE). Four of the B. petrii GIs (GI1, GI2, GI3 and GI6) are related to the clc element of Pseudomonas knackmussii B13, which codes for proteins involved in the biodegradation of chloroaromatics . The similarity between these four GIs and the clc element is mostly confined to the conjugal transfer region and varies substantially (42%–98%) between individual genes within this region [see Additional file 2]. However, among these GIs, GI3 is the most similar (78%–98%) to the clc element and also carries all the chloroaromatic degredation genes found in the clc element. Like the clc element in B13, these four GIs are integrated into the 3'-end of tRNAGly genes, are flanked by a 15–18 bp direct repeat generated by the duplication of the 3'-end of tRNA, and each encodes an integrase highly similar (88%, 69%, 100% and 52% amino acid identity) to that of the clc element (IntB13). GI1, GI2 and GI3 are located in close proximity within a 510 kb region (Figures 1, 3). The clc-type mobile elements are self transmissible and B. petrii may have obtained them from environmental Pseudomonads, with which it shares habitat . During the conjugal transfer process, ICEs are excised from the chromosome and form a circular intermediate . Interestingly, during in vitro culture the whole 510 kb region encompassing GI1, GI2 and GI3 is deleted in a fraction of the bacterial population. Sequence analysis of the deletion point revealed a regenerated 18 bp integration site (Figure 3), indicating that these GIs are still active in terms of excision from the chromosome. In fact, circular intermediates were detected for all GIs in B. petrii cells during normal growth (Lechner and Gross, unpublished results).
The remaining GIs (GI4, GI5 and GI7) encode a conjugal transfer system similar to the type IV secretion systems associated with plasmids. GI4 is highly similar to Tn4371 of Ralstonia eutropha A5/Cupriavidus oxalaticus, a 55-kb conjugative transposon that carries genes involved in the degradation of biphenyl or 4-chlorobiphenyl . The features of all GIs are summarized in Table 2.
It has been previously reported that the chromosomes of B. pertussis and B. parapertussis carry a large number of IS elements (261 and 112, respectively). These IS elements are responsible for the extensive shuffling and deletions in their genomes . In contrast, B. bronchiseptica and B. avium carry none. B. petrii carries 105 IS elements, most of which belong to the IS3-family. Interestingly, there are six copies of IS1663 which was found in B. pertussis and one truncated copy of IS1001 (Bpet1198) which was found in B. parapertussis. The vast majority (66%) of the B. petrii IS are located within GIs, whereas 34% are randomly dispersed around the chromosome. B. petrii also carries two intact prophages (phage 1, Bpet0942-1008; phage 2, Bpet4384-4432) (Figure 1), and phage-like particles were indeed observed by electron microscopy of the sequenced strain . Phage 1 is inserted at the 3'end of a gene encoding a tRNAcys, while phage 2 is related to bacteriophage Mu. In addition, there are several regions on the B. petrii genome which may encode prophage remnants (Figure 1). In contrast, B. avium carries three and B. bronchiseptica carries four prophages, however, none of them is related to the B. petrii phages.
The GIs (GI1-GI7), IS elements and prophages represent 22% of the B. petrii genome. This highly mosaic genome structure reflects extensive DNA exchange and provides first evidence for a considerable horizontal gene transfer-driven evolution in the genus Bordetella.
Metabolism of B. petrii - genomic basis for metabolic versatility
The central metabolic pathways are remarkably similar between B. petrii and the pathogenic Bordetellae. For instance, the glycolytic pathway is incomplete, lacking the determinants for glucokinase and phosphofructokinase. B petrii encodes all the enzymes needed for the TCA cycle and for the synthesis of all essential nucleotides, cofactors, amino acids and fatty acids. However, the in silico analysis of the genome sequence of B. petrii suggests the presence of a set of auxiliary pathways for the use of alternative nutrients, such as gluconate, degraded plant products, cyanate (as N-source), and various aromatic compounds (see below), which may contribute to the ability of B. petrii to thrive in different environments (Figure 4).
Gluconate is an intermediate in the degradation of sugar acids and sugar derivatives, which are components of more complex molecules, such as plant materials and secondary metabolites. The genome analysis of B. petrii suggests that this microorganism can also utilize gluconate via a variant of the Entner-Doudoroff (ED) pathway, which is only found in few bacteria, such as Rhodobacter sphaeroides, but absent in all the other sequenced Bordetellae . In this pathway, gluconate is first transformed by a gluconate dehydratase (Bpet3300) to 2-dehydro-3-deoxy-D-gluconate, which is then imported into the cytoplasm via 2-dehydro-3-deoxy-D-gluconate permease (Bpet1985) and phosphorylated to 2-dehydro-3-deoxy-D-gluconate-6-phosphate by a 2-dehydro-3-deoxy-D-gluconate kinase (KdgK, Bpet0876). This compound is then converted to glyceraldehyde 3-phosphate and pyruvate by the action of 2-dehydro-3-deoxy-D-gluconate-6-phosphate aldolase (KdgA, Bpet0875). A IclR family transcription regulator (Bpet0874) is present immediately upstream of the kdgA and kdgK genes suggesting that this variant ED pathway may be inducible by this regulator depending on the availability of the external carbon source. In this regard, it is surprising that although a periplasmic glucose dehydrogenase (gcd, Bpet4644), which converts glucose to gluconate, is encoded in the genome, B. petrii is unable to metabolize glucose .
Interestingly, B. petrii encodes a protein (Bpet0241) with significant similarity to pectate lyases, extracellular enzymes secreted by phytopathogens, such as Erwinia carotovora, which catalyze the cleavage of pectate [24–26]. The reaction products are channeled into the variant ED pathway . In agreement with this, the determinants for two other important enzymes of these pathways, 2-deoxy-D-gluconate 3-dehydrogenase (Bpet0046, Bpet1132) and altronate hydrolase (Bpet0414) are present, indicating that B. petrii uses plant cell-wall constituents as nutrient sources, which is in line with the isolation of a B. petrii strain from a plant root consortium . The presence of a cyanate transporter (Bpet3878) and a cyanate lyase (Bpet3621) catalyzing the decomposition of cyanate into ammonia and bicarbonate, suggests that B. petrii is not only able to cope with the toxicity of environmental cyanate, but like E. coli may use it as a nitrogen source . A striking difference between B. petrii and the other Bordetellae is that it has a facultative anaerobic energy metabolism and can utilize nitrate as the terminal electron acceptor during anaerobic respiration . Accordingly, on its core genome B. petrii encodes the factors (nitrate-, nitrite-, and nitrous oxide reductases) necessary to carry out complete denitrification from nitrate to N2. Interestingly, only some of the respective genes are present in the pathogenic Bordetellae. B. bronchiseptica appears to encode a functional nitrate reductase, while B. parapertussis has a truncated napD gene probably leading to a non-functional nitrate reductase, and B. pertussis does not encode genes involved in denitrification at all.
Metabolic talents of B. petrii - degradation of aromatic compounds
The in silico analysis of the B. petrii genome revealed an unusual metabolic capability regarding the degradation of aromatic compounds summarized in Figure 5. For the metabolisation of some of these compounds including benzoate, phthalate, 3-hydroxybenzoate, 3-hydroxyphenylacetate and 4-hydroxyphenylacetate we have obtained experimental evidence (data not shown).
A prominent feature is the presence of a surprisingly large number of genes coding for enzymes of the chloroaromatic metabolism. One complete set of chlorocatechol pathway genes (Bpet1533-1538) is localized on GI3 highly related to the clc element of Pseudomonas knackmussii B13. A second set of chlorocatechol pathway genes (Bpet3747-3752) is located on GI5 and shows high similarity to genes identified in the 1,2,4-trichlorobenzene degrading Pseudomonas strain P51 . A third possible chlorocatechol gene cluster encoding proteins only distantly related to those previously characterized is also located on GI5. This multiplicity of chlorocatechol gene clusters probably endows this strain with the ability to effectively metabolize a broad range of differently substituted chlorocatechols, providing thereby a valuable competitive advantage to thrive in polluted environments [30–32].
On GI5, chlorocatechol genes are preceded by a gene cluster (Bpet3738-3742) encoding a chlorobenzene dioxygenase and chlorobenzene dihydrodiol dehydrogenase, highly homologous to those of the 1,2,4-trichlorobenzene degraders Cupriavidus strain PS12 and Pseudomonas strain P51 [33, 34]. The absence of any linked extradiol dioxygenase encoding gene suggests that similarly to what is observed in the strains PS12 and P51, this cluster has specifically evolved for the transformation of chlorobenzenes [29, 30, 35].
B. petrii also harbors four different central pathways for aromatic metabolism via dihydroxylated central intermediates, the gentisate (Bpet0516-0517 and Bpet1429-1430, the latter located on GI2), the homogentisate (Bpet 4016–4017) and the catechol and protocatechuate branch of the 3-oxoadipate pathway (see below). Such versatility is similar to that observed in bacteria well known for their high metabolic capabilities, such as P. putida KT2440 and Cupriavidus necator . Only genes encoding enzymes of the homogentisate pathway were detected in all other sequenced Bordetellae, and those involved in the gentisate pathway are only present in B. avium.
Various natural aromatic compounds are known to be degraded via either the catechol or the protocatechuate branch of the 3-oxoadipate pathway, which thus plays a key role in bacterial aromatic catabolism. This pathway is widespread in Proteobacteria, however, all previously sequenced Bordetella strains are devoid of complete pathways and only a catIJCD cluster (Bpet4528-4532) is present as a remnant. In B. petrii, complete pathways are formed by genes of the protocatechuate branch (Bpet0260-0254) located in a genomic region with a quite low GC content, but lacking other typical features of a genomic island (named GI in Table 2 and in Additional file 3), and two clusters localized on GI2, encoding the enzymes of the catechol branch. We believe this to be the first example of natural horizontal transfer of this central aromatic metabolic pathway. Notably, also peripheral pathways for funneling benzoate (Bpet1396-1399) and benzylalcohol (Bpet1387-1388) to catechol are located on GI2. A peripheral pathway for channeling phthalate to protocatechuate (Bpet0248-0252) is located on GI1. The low number of peripheral routes contrasts with the situation in other Proteobacteria, such as B. xenovorans LB400, where five substrates are degraded via protocatechuate .
A strategy of aerobic degradation of aromatics alternative to the metabolism via dihydroxylated intermediates is initiated by CoA ligases [38, 39]. B. petrii carries gene clusters required for the metabolism of all three substrates, which are known to be degraded by such a strategy (Bpet0549-0554 for anthranilate metabolism and Bpet1924-1935 for phenylacetate metabolism, which are present in other Bordetellae, and Bpet3568-3575 for benzoate metabolism). As degradation via CoA derivatives needs less oxygen molecules compared to degradation via dihydroxylated intermediates, this may provide B. petrii with a selective advantage to also thrive in oxygen limited environments. For an overview about the aromatic degradation capabilities of B. petrii see Figure 5 and Additional file 3.
"Virulence factors" of B. petrii and host-restricted Bordetellae
Pathogenic Bordetellae produce many surface associated or secreted virulence factors . B. petrii does not code for the well-characterized toxins of classical Bordetellae (i.e., pertussis toxin, adenylate cyclase toxin and dermonecrotic toxin). However, it produces a low amount (22–53 nM/OD) of the non-protein tracheal cytotoxin (TCT), a spontaneously released muramylpeptide that disrupts epithelia . This expression range is comparable to the expression levels observed for the broad host range pathogen B. bronchiseptica.
Autotransporter or type V secretion proteins play a major role in Bordetella virulence . B. bronchiseptica encodes 21 autotransporter proteins, whereas B. petrii encodes only four (Bpet0190, Bpet0193, Bpet3980, Bpet4156) . Two of them (Bpet0190 and Bpet0193) are paralogous and show significant homology to the putative serine protease autotransporter BB2301 of B. bronchiseptica, whereas the other two do not show significant similarity to any database entry. Pathogenic Bordetellae also encode a two-partner secretion system which allows secretion of the colonization factor filamentous hemagglutinin (FhaB) . This protein requires the FhaC protein for translocation across the outer membrane. B. petrii encodes proteins related to FhaB and FhaC (Bpet600-601). B. petrii FhaB has 3039 amino acids and its N-terminus, which is required for export, is well-conserved with respect to pathogenic Bordetellae. However, the overall similarity with FhaB of pathogenic Bordetellae is very low, being more similar to putative secreted proteins of plant pathogens (e.g., Ralstonia, Xylella and Burkholderia). This suggests that the fha genes found in B. petrii and the pathogenic Bordetellae may not be orthologous.
Other protein secretion systems are differentially distributed among Bordetellae. In contrast to the mammalian pathogens, B. petrii and B. avium lack type I and III protein secretion systems, but encode a type II secretion system (Bpet2855-2866 and BAV0331-0343). Type IV secretion systems are either engaged in protein transport to the extracellular medium or into eukaryotic cells, or they are involved in DNA-transfer. There are three gene clusters in B. petrii, which are associated with GIs and appear to encode complete type IV secretion systems (Bpet2204-2215, Bpet3701-3719, Bpet4615-4626). None of these systems appears to be related to the pertussis toxin exporting type IV system of B. pertussis, being probably involved in conjugal DNA transfer.
The B. petrii genome carries three fimbrial gene clusters (Bpet0111-0114, Bpet4129-4132 and Bpet4464-4468), each encoding a major fimbrial subunit, a fimbrial assembly chaperone, an outer membrane usher protein and a fimbrial adhesin. The chaperone and usher proteins have significant sequence similarities with their counterparts in the pathogenic Bordetellae, with e-values ranging from 6e-23 to 3e-35 for the chaperones and e-values ranging from 3e-71 to 6e-103 for the usher proteins. In contrast, the homology of the major structural proteins and fimbrial adhesins with the proteins of the pathogenic Bordetellae is very low (e-values > 1e-5). Similarly to B. avium, B. petrii lacks the genes for the biosynthesis of siderophores. However, it specifies 12 TonB-dependent iron scavenging receptors, which may be important for iron acquisition in soil and plant-associated environments.
Most virulence factors are coordinately regulated in pathogenic Bordetellae by a single master regulator, the BvgAS two-component system [42, 43]. This system consists of two proteins, a response regulator (BvgA) and a multidomain histidine kinase (BvgS), consisting of a periplasmic solute-binding domain, and, on the cytoplasmic side, of a PAS, a transmitter, a receiver and an HPt domain. Although the B. petrii BvgA (Bpet 4471) is highly homologous to the BvgA proteins of the pathogenic Bordetellae (65% amino acid similarity), the structure of the B. petrii equivalents of BvgS is more complex. First of all, there are two histidine kinases, BvgS1 (Bpet4469) and BvgS2 (Bpet4472), with significant sequence similarities to BvgS. BvgS1 contains a PAS, a transmitter and a receiver domain, whereas BvgS2 contains a transmitter and a receiver domain. In contrast to BvgS, both B. petrii proteins are devoid of a large periplasmic domain. The HPt domain, which in pathogenic Bordetellae is part of the multidomain BvgS protein, in B. petrii is an independent protein (Bpet4470) (see Additional file 4. To see legends for Additional files 1, 2 and 4 please see Additional file 5 ). This suggests that during their evolution to pathogens, the virulence regulatory system has been streamlined by the omission of one histidine kinase and incorporation of the HPt domain into the multidomain histidine kinase BvgS.
Therefore, B. petrii is endowed with several putative virulence factors for plant and animal hosts, such as adhesins (e.g., Fha, fimbriae) and a master virulence regulator (Bvg), but lacks Bordetella toxins. This suggests that pathogenic Bordetellae might have acquired their toxin genes horizontally, while adapting to an obligate pathogenic lifestyle. Alternatively, it is possible that toxin genes have been lost from the B. petrii lineage.
The five sequenced Bordetella species, B. bronchiseptica, B. pertussis, B. parapertussis, B. avium and B. petrii, represent different states of environmental adaptation. These five organisms have followed different evolutionary paths. Evolution of the B. bronchiseptica-derived host restricted-pathogens B. pertussis and B. parapertussis was dominated by substantial gene decay and loss. This might also hold true for B. avium, which is characterized by a relatively small genome containing 68 pseudogenes. In contrast, the evolution of the environmental isolate B. petrii was dominated by horizontal acquisition of large genomic islands expanding its metabolic capacities. B. bronchiseptica, with a broad host range and the potential to survive in the environment, might represent an intermediate evolutionary state. If the Bordetellae derive from a common environmental ancestor, B. petrii appears to be more closely related to this ancestor than the pathogenic Bordetellae, since it is the only true environmental isolate within the genus, which is also found in association with host organisms. Moreover, B. petrii is the only facultative anaerobic Bordetella species . Since closely related Achromobacter species are also facultative anaerobes, the common ancestor of the Bordetellae probably had a facultative anaerobic energy metabolism. During adaptation to an obligate host association, the genes for several metabolic functions, e.g. those required for denitrification, may have been deleted from the genomes of pathogenic Bordetellae, leaving just a few remnants in the various species.
Whole genome shotgun sequencing
For the shotgun phase we produced 103,308 paired-end sequences derived from three pTZ18R genomic libraries (with insert sizes of 1 kb, 3 kb and 10 kb) yielding 11.72-fold coverage, using dye terminator chemistry on ABI377 and MegaBACE 1000 automated sequencers. As a scaffold we used 2,688 paired-end sequences pCC1BAC libraries with insert sizes of 17–30 kb and 40–70 kb (leading to a total coverage of 5.14-fold). Assembly of the paired-end sequences was performed by PHRAP (P. Green, unpublished data; http://www.phrap.org/). The complete genome sequence was finished by using the Staden package (GAP v4.8b1) . Another 5,154 sequencing reads were generated to close the remaining gaps and to polish the sequence.
Genome analysis and annotation
Curation and annotation of the genome was done using the genome annotation system GenDB . A combined gene prediction strategy was applied on the assembled sequences using GLIMMER and CRITICA . Putative ribosomal binding sites and transfer RNA genes were identified with RBSFINDER and tRNAscan-SE, respectively [47, 48]. In a first step, automatic annotation was computed based on various tools: similarity searches were performed against different databases including SWISS-PROT, TrEMBL, KEGG, Pfam, TIGRFAM and InterPro. Furthermore, SignalP, helix-turn-helix and TMHMM were applied. Additionally, each gene was functionally classified by assigning of a Cluster of Orthologous Groups (COG) number and its corresponding COG category and Gene Ontology numbers [49, 50]. Finally, the annotation of each predicted gene was manually inspected and adjusted if necessary.
For comparative analyses, the annotated genome sequences of the following bacteria were imported into the genome annotation system GenDB : Bordetella avium 197N (gb;AM167904), Bordetella pertussis Tohama I (gb;NC_002929), Bordetella parapertussis 12822 (gb;NC_002928), Bordetella bronchiseptica RB50 (gb;NC_002927). Similarity searches were conducted against the genomes on the nucleotide and amino acid sequence level by using BLASTN and BLASTP, respectively . Synteny plots were created with GenDB . For the determination of the core genome common to all sequenced Bordetellae all the CDS were compared to each other reciprocally using BLASTP with an e-value < 1e-20 as a cut-off.
Detection of regions with atypical GC content
Genomic regions with atypical GC content were identified using the 'sliding window' technique with a window size of 3,000 bp and a step size of 1,000 bp. For this purpose, the GC content was assumed to follow a Gaussian distribution and regions with at least 2.0 standard deviation differences from the mean were determined.
Use of aromatic compounds as sole carbon source by B. petrii
Growth experiments were carried out in a minimal medium containing 10% salt solution (2.5 g NaCl, 0.5 g KH2PO4, 0.2 g KCl, 0.1 g MgCl2 × H2O, 1.5 g Tris/HCl resuspended in 100 ml distilled water), 1.8 μM CaCl2, 36 μM FeSO4 × H2O, 2 mM (NH4)SO4. Aromatic carbon sources were added to final concentrations of 2–5 mM each.
Determination of TCT production
Tracheal Cytotoxin (TCT) was quantified as previously described .
The complete, annotated genome sequence has been submitted to EBI under the accession number AM902716.
This work was supported by the Pathogenomics Competence Network of the BMBF and by the priority research programme SFB479/A2 from the Deutsche Forschungsgemeinschaft.
Chair of Microbiology, Biocenter, University of Würzburg
Department of Vaccinology and Applied Microbiology, Helmholtz Center for Infection Research
Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus
Division of Molecular Biotechnology, Helmholtz Center for Infection Research
Department of Microbial Pathogenesis, Helmholtz Center for Infection Research
Institut de Recherche pour le Développement UMR 5096, CNRS-UP-IRD 911
Center for Biotechnology (CeBiTec), Bielefeld University
Chair of Genetics, Bielefeld University
Institute for Microbiology and Hygiene
Department of Molecular Microbiology, Washington University School of Medicine
Department of Microbiology and Immunology, University of North Carolina
Department of Genome Analysis, Helmholtz Center for Infection Research
Mathematics and Computer Science Division, Argonne National Laboratory
Roche Diagnostics GmbH
BASF Plant Science GmbH
Mattoo S, Cherry JD: Molecular pathogenesis, epidemiology, and clinical manifestations of respiratory infections due toBordetella pertussisand otherBordetellasubspecies.Clin Microbiol Rev 2005, 18:326–382.View ArticlePubMed
Gerlach G, von Wintzingerode F, Middendorf B, Gross R: Evolutionary trends in the genusBordetella
. Microbes Infect 2001, 3:61–72.View ArticlePubMed
Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, Holden MT, Churcher CM, Bentley SD, Mungall KL, Cerdeño-Tárraga AM, Temple L, James K, Harris B, Quail MA, Achtman M, Atkin R, Baker S, Basham D, Bason N, Cherevach I, Chillingworth T, Collins M, Cronin A, Davis P, Doggett J, Feltwell T, Goble A, Hamlin N, Hauser H, Holroyd S, Jagels K, Leather S, Moule S, Norberczak H, O'Neil S, Ormond D, Price C, Rabbinowitsch E, Rutter S, Sanders M, Saunders D, Seeger K, Sharp S, Simmonds M, Skelton J, Squares R, Squares S, Stevens K, Unwin L, Whitehead S, Barrell BG, Maskell DJ: Comparative analysis of the genome sequences ofBordetella pertussis, Bordetella parapertussisandBordetella bronchiseptica. Nat Genet 2003, 35:32–40.View ArticlePubMed
Sebaihia M, Preston A, Maskell DJ, Kuzmiak H, Connell TD, King ND, Orndorff PE, Miyamoto DM, Thomson NR, Harris D, Goble A, Lord A, Murphy L, Quail MA, Rutter S, Squares R, Squares S, Woodward J, Parkhill J, Temple LM: Comparison of the genome sequence of the poultry pathogenBordetella aviumwith those ofB. bronchiseptica
, B. pertussis, andB. parapertussisreveals extensive diversity in surface structures associated with host interaction.J Bacteriol 2006, 188:6002–6015.View ArticlePubMed
Cummings CA, Brinig MM, Lepp PW, Pas S, Relman DA:Bordetellaspecies are distinguished by patterns of substantial gene loss and host adaptation.J Bacteriol 2004, 186:1484–1492.View ArticlePubMed
von Wintzingerode F, Schattke A, Siddiqui RA, Rösick U, Göbel UB, Gross R:Bordetella petriisp. nov., isolated from an anaerobic bioreactor, and emended description of the genusBordetella
. Int J Syst Evol Microbiol 2001, 51:1257–1265.
Bianchi F, Careri M, Mustat L, Malcevschi A, Musci M: Bioremediation of toluene and naphthalene: development and validation of a GC-FID method for their monitoring.Ann Chim 2005, 95:515–524.View ArticlePubMed
Wang F, Grundmann S, Schmid M, Dörfler U, Roherer S, Charles Munch J, Hartmann A, Jiang X, Schroll R: Isolation and characterization of 1,2,4-trichlorobenzene mineralizingBordetellasp. and its bioremediation potential in soil.Chemosphere 2007, 67:896–902.View ArticlePubMed
Fry NK, Duncan J, Malnick H, Warner M, Smith AJ, Jackson MS, Ayoub A:Bordetella petriiclinical isolate.Emerg Infect Dis 2005, 11:1131–1133.PubMed
Stark D, Riley LA, Harkness J, Marriott D:Bordetella petriifrom a clinical sample in Australia: isolation and molecular identification.J Med Microbiol 2007, 56:435–437.View ArticlePubMed
Spilker T, Liwienski AA, LiPuma JJ: Identification ofBordetellaspp. in respiratory specimens from individuals with cystic fibrosis.Clin Microbiol Infect 2008, 14:504–506.View ArticlePubMed
Sfanos K, Harmody D, Dang P, Ledger A, Pomponi S, McCarthy P, Lopez J: A molecular systematic survey of cultured microbial associates of deep-water marine invertebrates.Syst Appl Microbiol 2005, 28:242–264.View ArticlePubMed
Chowdhury SP, Schmid M, Hartmann A, Tripathi AK: Identification of diazotrophs in the culturable bacterial community associated with roots ofLasiurus sindicus, a perennial grass of Thar desert, India.Microb Ecol 2007, 54:82–90.View ArticlePubMed
Gomez-Cerezo J, Suárez I, Ríos JJ, Peña P, García de Miguel MJ, de José M, Monteagudo O, Linares P, Barbado-Cano A, Vázquez JJ:Achromobacter xylosoxidansbacteremia: a 10-year analysis of 54 cases.Eur J Clin Microbiol Infect Dis 2003, 22:360–363.View ArticlePubMed
Conway SP, Brownlee KG, Denton M, Peckham DG: Antibiotic treatment of multidrug-resistant organisms in cystic fibrosis.Am J Respir Med 2005,2(4):321–32.
Rehfuss M, Urban J:Alcaligenes faecalissubsp.phenolicussubsp. nov. a phenol-degrading, denitrifying bacterium isolated from a graywater bioprocessor.Syst Appl Microbiol 2005, 28:421–429.View ArticlePubMed
Vinas M, Sabate J, Espuny MJ, Solanas AM: Bacterial community dynamics and polycyclic aromatic hydrocarbon degradation during bioremediation of heavily creosote-contaminated soil.Appl Environ Microbiol 2005, 71:7008–7018.View ArticlePubMed
Gaillard M, Vallaeys T, Vorhölter FJ, Minoia M, Werlen C, Sentchilo V, Pühler A, Meer JR: Theclcelement ofPseudomonassp. strain B13, a genomic island with various catabolic properties.J Bacteriol 2006, 188:1999–2013.View ArticlePubMed
Meer JR, Sentchilo V: Genomic islands and the evolution of catabolic pathways in bacteria.Curr Opin Biotechnol 2003, 14:248–254.View ArticlePubMed
Ravatn R, Studer S, Zehnder AJ, Meer JR: Int-B13, an unusual site-specific recombinase of the bacteriophage P4 integrase family, is responsible for chromosomal insertion of the 105-kilobaseclcelement ofPseudomonassp. strain B13.J Bacteriol 1998, 180:5505–5514.PubMed
Toussaint A, Merlin C, Monchy S, Benotmane MA, Leplae R, Mergeay M, Springael D: The biphenyl- and 4-chlorobiphenyl-catabolic transposon Tn a member of a new family of genomic islands related to IncP and Ti plasmids.Appl Environ Microbiol 4371, 69:4837–4845.View Article
Kim S, Lee SB: Identification and characterization ofSulfolobus solfataricusD-gluconate dehydratase: a key enzyme in the non-phosphorylated Entner-Doudoroff pathway.Biochem J 2005, 387:271–280.View ArticlePubMed
Collmer A, Keen NT: The role of pectic enzymes in plant pathogenesis.Ann Rev Phytopathol 1986, 24:383–409.View Article
Kotoujansky A: Molecular genetics of pathogenesis by soft rot Erwinias.Ann Rev Phytopathol 1987, 25:405–430.View Article
Bartling S, Wegener C, Olsen O: Synergism betweenErwiniapectate lyase isozymes that depolymerise both pectate and pectin.Microbiology 1995, 141:873–881.View ArticlePubMed
Condemine G, Hugouvieux-Cotte-Pattat N, Robert-Baudouy J: Isolation ofErwinia chrysanthemi kduDmutants altered in pectin degradation.J Bacteriol 1993,165(3):937–41.
Sung YC, Fuchs JA: Identification and characterization of a cyanate permease inEscherichia coliK-12.J Bacteriol 1998,171(9):4674–8.
Meer JR, Frijters AC, Leveau JH, Eggen RI, Zehnder AJ, de Vos WM: Characterization of thePseudomonassp. strain P51 genetcbR, a LysR-type transcritional activator of thetcbCDEFchlorocatechol oxidative operon, and analysis of the regulatory region.J Bacteriol 1991, 173:3700–3708.PubMed
Meer JR, Zehnder AJB, Vos WM: Identification of a novel composite transposable element, Tn carrying chlorobenzene dioxygenase genes ofPseudomonassp. strain P51.J Bacteriol 5280, 173:7077–7083.
Potrawfke T, Armengaud J, Wittich RM: Chlorocatechols at positions 4 and 5 are substrates of the broad-spectrum chlorocatechol 1,2-dioxygenase Pseudomonas chlororaphis RW71.J Bacteriol 2001, 183:997–1011.View ArticlePubMed
Liu S, Ogawa N, Senda T, Hasebe A, Miyashita K: Amino acids in positions 48, 52, and 73 differentiate the substrate specificities of the highly homologous chlorocatechol 1,2-dioxygenases CbnA and TcbC.J Bacteriol 2005, 187:5427–5436.View ArticlePubMed
Beil S, Happe B, Timmis KN, Pieper DH: Genetic and biochemical characterization of the broad-spectrum chlorobenzene dioxygenase fromBurkholderiasp. strain PS12: Dechlorination of 1,2.4,5-tetrachlorobenzene.Eur J Biochem 1997, 247:190–199.View ArticlePubMed
Raschke H, Meier M, Burken JG, Hany R, Müller MD, Meer JR, Kohler HP: Biotransformation of various substituted aromatic compounds to chiral dihydrodihydroxy derivatives.Appl Env Microbiol 2001, 67:3333–3339.View Article
Beil S, Timmis KN, Pieper DH: Genetic and biochemical analyses of thetecoperon suggest a route for evolution of chlorobenzene degradation genes.J Bacteriol 1999, 181:341–346.PubMed
Jimenez JI, Minambres B, Garcia JL, Diaz E: Genomic analysis of the aromatic catabolic pathways fromPseudomonas putidaKT2440.Env Microbiol 2002, 4:824–841.View Article
Chain PS, Denef VJ, Konstantinidis KT, Vergez LM, Agulló L, Reyes VL, Hauser L, Córdova M, Gómez L, González M, Land M, Lao V, Larimer F, LiPuma JJ, Mahenthiralingam E, Malfatti SA, Marx CJ, Parnell JJ, Ramette A, Richardson P, Seeger M, Smith D, Spilker T, Sul WJ, Tsoi TV, Ulrich LE, Zhulin IB, Tiedje JM:Burkholderia xenovoransLB400 harbors a multi-replicon, 9.73-Mbp genome shaped for versatility.Proc Natl Acad Sc USA 2006, 103:15280–15287.View Article
Mohamed ME, Fuchs G: Purification and characterization of phenylacetate-coenzyme A ligase from a denitrifyingPseudomonassp., an enzyme involved in the anaerobic degradation of phenylacetate.Arch Microbiol 1993, 159:554–562.View Article
Altenschmidt U, Fuchs G: Novel aerobic 2-aminobenzoate metabolism. Purification and characterization of 2-aminobenzoate-CoA ligase, localisation of the genes on a 8-kbp plasmid, and cloning and sequencing of the genes from a denitrifyingPseudomonassp.Eur J Biochem 1992, 205:721–727.View ArticlePubMed
Flak TA, Heiss LN, Engle JT, Goldman WE: Synergistic epithelial responses to endotoxin and a naturally occurring muramyl peptide.Infect Immun 2000, 68:1235–1242.View ArticlePubMed
Locht C, Antoine R, Jacob-Dubuisson F:Bordetella pertussis, molecular pathogenesis under multiple aspects.Curr Opin Micobiol 2001, 4:82–89.View Article
Beier D, Gross R: Regulation of bacterial virulence by two-component systems.Curr Opin Microbiol 2006, 9:143–152.View ArticlePubMed
Cummings CA, Bootsma HJ, Relman DA, Miller JF: Species- and strain-specific control of a complex, flexible regulon byBordetellaBvgAS.J Bacteriol 2006, 188:1775–1785.View ArticlePubMed
Staden R, Beal KF, Bonfield JK: The Staden package.Methods Mol Biol 2000, 132:115–130.PubMed
Meyer F, Goesmann A, McHardy AC, Bartels D, Bekel T, Clausen J, Kalinowski J, Linke B, Rupp O, Giegerich R, Pühler A: GenDB - an open source genome annotation system for prokaryote genomes.Nucleic Acids Res 2003, 31:2187–2195.View ArticlePubMed
McHardy AC, Goesmann A, Pühler A, Meyer F: Development of joint application strategies for two microbial gene finders.Bioinformatics 2004, 20:1622–1631.View ArticlePubMed
Suzek BE, Ermolaeva MD, Schreiber M, Salzberg SL: A probabilistic method for identifying start codons in bacterial genomes.Bioinformatics 2001, 17:1123–1130.View ArticlePubMed
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence.Nucleic Acids Res 1997, 25:955–964.View ArticlePubMed
Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA: The COG database: an updated version includes eukaryotes.BMC Bioinformatics 2003, 4:41.View ArticlePubMed
Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS, Sethuraman A, Theesfeld CL, Botstein D, Dolinski K, Feierbach B, Berardini T, Mundodi S, Rhee SY, Apweiler R, Barrell D, Camon E, Dimmer E, Lee V, Chisholm R, Gaudet P, Kibbe W, Kishore R, Schwarz EM, Sternberg P, Gwinn M, Hannick L, Wortman J, Berriman M, Wood V, de la Cruz N, Tonellato P, Jaiswal P, Seigfried T, White R, Gene Ontology Consortium: The Gene Ontology (GO) database and informatics resource.Nucleic Acids Res 2004, 32:D258–261.View ArticlePubMed
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.Nucleic Acids Res 1997, 25:3389–3402.View ArticlePubMed
Cookson BT, Cho HL, Herwaldt LA, Goldman WE: Biological activities and chemical composition of purified tracheal cytotoxin ofBordetella pertussis
. Infect Immun 1989, 57:2223–2229.PubMed
Antoine R, Jacob-Dubuisson F, Drobecq H, Willery E, Lesjean S, Locht C: Overrepresentation of a gene family encoding extracytoplasmic solute receptors inBordetella
. J Bacteriol 2003, 185:1470–1474.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.