Previous Article | Next Article ![]()
Journal of Bacteriology, December 2006, p. 8385-8394, Vol. 188, No. 24
0021-9193/06/$08.00+0 doi:10.1128/JB.01081-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
,
,
C. A. Cummings,3,4,
H. G. J. van der Heide,1
M. van Gent,1
S. Liew,3,4
D. A. Relman,3,4,5 and
F. R. Mooi1*
Laboratory for Vaccine Preventable Diseases, National Institute of Public Health and the Environment, Bilthoven, The Netherlands,1 Eijkman Winkler Institute, University Medical Center, Utrecht, The Netherlands,2 Department of Microbiology and Immunology, Stanford University School of Medicine, Stanford, California,3 VA Palo Alto Health Care System, Palo Alto, California,4 Department of Medicine, Stanford University School of Medicine, Stanford, California5
Received 21 July 2006/ Accepted 2 October 2006
|
|
|---|
|
|
|---|
In addition to B. holmesii, a number of other pathogenic species are found in the genus Bordetella. Best studied are the respiratory pathogens known as the classical or mammalian bordetellae: B. pertussis, human B. parapertussis, ovine B. parapertussis, and B. bronchiseptica. B. pertussis and human B. parapertussis are the causative agents of whooping cough, or pertussis, a disease which has reemerged despite widespread vaccination (8, 42). Ovine B. parapertussis infects the respiratory tracts of sheep, and B. bronchiseptica infects a wide range of mammalian species. The nonclassical Bordetella species include B. avium, B. hinzii, B. petrii, "B. ansorpii" (proposed name), and B. trematum, and these form a genetically diverse group that is clearly distinct from the mammalian bordetellae (16, 21).
Comparative analysis of their 16S rRNA sequences suggested that B. holmesii is very closely related to B. pertussis (48), a hypothesis supported by the discovery of the B. pertussis-specific insertion sequence element IS481 in B. holmesii (34). However, subsequent sequencing of housekeeping genes, analysis of cellular fatty acid composition, and characterization of the bvgAS locus suggested that B. holmesii may not be as closely related to B. pertussis as was first assumed (15, 48).
To identify genetic factors that may have contributed to the emergence of B. holmesii as a human pathogen and to clarify its relation to other members of the Bordetella genus, we used microarray-based comparative genome hybridization (CGH) and sequencing of housekeeping genes. Although gene sequencing indicated that B. holmesii is a uniform group that is more closely related to the nonclassical B. hinzii and B. avium than to the mammalian bordetellae, CGH detected a genomic region of approximately 66 kb that is highly conserved between B. pertussis and B. holmesii. This genomic island, encoding genes that could promote pathogenesis, may have been transferred from B. pertussis to B. holmesii, possibly contributing to the emergence of B. holmesii as a human pathogen.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. B. holmesii strain characteristics
|
For each spot on the array, the fluorescence intensity of each of the two dyes was calculated by subtracting the mean intensity of the background pixels from the median intensity of the spot pixels. The logarithm of the ratio of background-subtracted Cy5 intensity over background-subtracted Cy3 intensity [log2(Cy5/Cy3)] is near 0 when both samples hybridize equally well to the probe and significantly below 0 when the B. holmesii DNA hybridizes poorly compared to the reference DNA. Data normalization was achieved by calculating the mean log2(Cy5/Cy3) for 16 probes within the iron uptake island (IUI) that were verified by sequencing to be at least 99% identical to the B. holmesii genome (see Table S1 in the supplemental material) and then subtracting this value from the log2(Cy5/Cy3) for each probe on the array. Data from replicate spots for the same probe were averaged. A probe was considered to have a positive B. holmesii hybridization signal when its hybridization ratio exceeded the fifth percentile of ratios from a control hybridization of B. pertussis Tohama I. Further details are in Methods in the supplemental material.
Organization of the IUI. The IUI was analyzed by amplifying overlapping 5- to 10-kb DNA fragments across the region. PCRs were analyzed by agarose gel electrophoresis to compare the fragment sizes from B. holmesii and B. pertussis Tohama I (see Table S2 in the supplemental material).
A total of 21 genes in the IUI were partially sequenced from three B. holmesii isolates. The nucleotide sequence of the 4.8-kb B. holmesii insert, absent from the classical bordetellae, was determined by genome walking inwards from the adjacent genes. In order to assess the organization of the BB3888 homolog, one of the open reading frames (ORFs) contained in the IUI, in B. pertussis and B. holmesii, we designed PCR primers in BB3888 that span the IS481 insertion and rearrangement point from B. pertussis Tohama I. Nucleotide sequence data adjacent to the left breakpoint of the IUI was obtained by PCR, using one primer inside the IUI, in BB0794, and one primer outside the IUI, in BB0795, which was not detected in B. holmesii by CGH. The sequence adjacent to the right breakpoint of the IUI was obtained using the TOPO walker kit (Invitrogen, Carlsbad, CA). Further details are in Methods in the supplemental material, and primers are listed in Table S3 in the supplemental material.
Phylogenetic analysis and broad-range PCR. PCR primers for phylogenetic analysis were designed to hybridize to conserved regions of the atpD, rpoB, tuf, and rnpB genes (see Table S3 in the supplemental material). PCR products were directly sequenced. Concatenated sequences of these four genes were aligned for the construction of neighbor-joining trees using PHYLIP software (J. Felenstein. 2005. Phylogeny Inference Package, version 3.6 [distributed by the author]. Department of Genomic Sciences, University of Washington, Seattle). Details are in Methods in the supplemental material.
B. holmesii 16S rRNA genes were amplified by broad-range PCR using the 8F and 806R primers (35), and 23S rRNA genes were amplified using the MS37 and MS38 primers (22). PCR products were cloned using the Topo TA cloning kit (Invitrogen). Twelve clones were picked for each strain, and plasmids were purified using the Wizard miniprep kit (Promega) and sequenced using M13 forward and reverse primers.
Southern blot hybridization. The genomic DNA of B. holmesii, B. pertussis, and B. avium isolates was digested with the restriction enzymes ClaI and NcoI. Southern blotting and hybridization with biotin-labeled oligonucleotide probes were performed essentially as described by Schouls et al. (37). The same blot was used for both 16S rRNA gene probe hybridizations, with a stripping step between experiments. Table S3 in the supplemental material lists the probes that were hybridized to the genomes of the Bordetella isolates.
Transcription of alcaligin genes. B. holmesii strains were grown to mid-logarithmic phase in iron-depleted or iron-replete medium. Total RNA was isolated using the QIAGEN RNeasy kit (QIAGEN Inc, Valencia, CA). Real-time PCR experiments were performed with a LightCycler (Roche Diagnostics, Indianapolis, IN) using the LightCycler DNA master hybridization probe kit (Roche Diagnostics). alcA and fauA transcript levels were compared to that of rpoB, which served as the internal reference. Further details are in Methods in the supplemental material.
MS. Culture supernatant samples were analyzed by nanoscale reversed-phase liquid chromatography (LC) on a model HP1100 LC system (Hewlett Packard GmbH, Waldbronn, Germany) coupled to an electrospray mass spectrometer (MS) (liquid chromatography quadrupole classic ion trap), essentially as described by Meiring et al. (25). Technical details are in Methods in the supplemental material.
Microarray data accession number. Raw microarray data have been deposited at ArrayExpress under accession number E-TABM-55.
Nucleotide sequence accession numbers. GenBank accession numbers DQ402382 to DQ402420 and DQ420056 to DQ420164 are published for the first time.
|
|
|---|
![]() View larger version (47K): [in a new window] |
FIG. 1. Comparative genome hybridization of B. holmesii strains to a classical Bordetella DNA microarray. (A) CGH of 12 B. holmesii isolates to a microarray comprising the genomes of B. pertussis Tohama I, B. parapertussis 12822, and B. bronchiseptica RB50. The moving average (with a sliding window of three adjacent values) of the mean log2(Cy5/Cy3) of 12 B. holmesii genomes is plotted on the x axis. Microarray probes that are represented in the B. pertussis Tohama I genome are arranged on the y axis in B. pertussis Tohama I genome order. The dashed line indicates the fifth percentile of B. pertussis Tohama I hybridization ratios; ratios below the fifth percentile are considered to indicate the lack of B. holmesii hybridization. (B) Probes that hybridized to the B. holmesii genome with a strength comparable to that of the reference and adjacent nonhybridizing probes are shown in detail for individual B. holmesii strains and for B. pertussis Tohama I. ORF and gene designations for selected loci are provided as landmarks. Probes representing insertion sequence elements are also shown. BB2492 is a putative IS3 family transposase found in B. bronchiseptica RB50 but not in B. pertussis. Strain numbers are indicated above the columns. Each row represents one probe in B. pertussis Tohama I gene order. The logarithm of the hybridization ratio [log2(Cy5/Cy3)] is indicated by the yellow-black-blue color scale. Missing data are represented in gray. The data represented in this figure are available as Table S4 in the supplemental material.
|
Molecular characterization of the B. holmesii IUI. Partial nucleotide sequences of 21 genes in the IUI, indicated to be present by CGH, were obtained from B. holmesii isolates B0436, B1850, and B1852 after PCR amplification using B. pertussis primer pairs (GenBank accession numbers DQ402382 to DQ402418). All of these sequences were at least 99.4% identical to the B. pertussis sequence, with the exception of BP0794 (94.1%), located at the putative right end of the IUI (see below; see also Table S1 in the supplemental material).
The GC content of the IUI (66.0%) is substantially higher than the average GC content of B. holmesii genes (61.5 to 62.3%) (48; see below) but comparable to that of the B. pertussis genome (67.7%) (30).
In B. bronchiseptica and B. parapertussis, but not in B. pertussis, the two genomic regions identified by CGH constitute a single continuous chromosomal fragment. The presence of IS481 sequences at the putative rearrangement breakpoint in BP0787 and BP2450, which are the two halves of a single ORF in B. bronchiseptica and B. parapertussis, suggests that this region has been rearranged in the B. pertussis lineage by recombination between transposon copies (Fig. 2). PCR was performed using a primer pair that spans the breakpoint (BB3888-480f and BB3888-720r; see also Table S3 in the supplemental material). This reaction yielded a 240-bp fragment from B. bronchiseptica RB50 but failed to amplify a product from B. pertussis Tohama I due to the chromosomal rearrangement between these primers.
![]() View larger version (17K): [in a new window] |
FIG. 2. Comparison of IUI genomic organizations in B. holmesii, B. pertussis, B. bronchiseptica, and B. avium. Deletions and insertions between species are indicated by gray surfaces. The dashed lines connect the ORFs at the borders of the orthologous sequences. The ORFs are color coded as described in the key. The ORF composition of B. holmesii was deduced from PCR and CGH data, while the representations of B. avium 197N, B. pertussis Tohama I, and B. bronchiseptica RB50 ORFs were based on published genome sequences.
|
Fourteen PCR primer pairs that amplify overlapping fragments from the B. pertussis Tohama I IUI sequence were used to map the IUIs in 12 B. holmesii strains by PCR amplification (see Table S2 in the supplemental material). As shown by agarose gel electrophoresis, most B. holmesii PCR amplicons were identical in size to those from B. pertussis, suggesting conserved genomic organization, but small insertions or deletions or point mutations that could result in pseudogene formation would not have been detected by electrophoresis of PCR amplicons or CGH. Two primer pairs yielded PCR products from B. holmesii that were approximately 1 kb shorter than the corresponding Tohama I products; sequencing of these products indicated that the IS481 elements BP2485 and BP2492 were missing in B. holmesii. Like B. holmesii, B. pertussis 18323 lacked these two IS481 elements. A third IS481 element, BP2477, was variably present among the B. holmesii strains.
PCR amplification of one genomic region in the IUI, downstream of the alc operon, revealed a 4.8-kb insertion in B. holmesii which is not present in the classical bordetellae (Fig. 2). Nucleotide sequencing identified the left terminus of this insert 237 nucleotides downstream of the BP2464 stop codon and its right terminus in BP2465, which has been partially deleted (GenBank accession number DQ402419).
BLASTX comparison of the 4.8-kb B. holmesii-specific sequence against GenBank identified five putative genes, designated bhoABCDE, with highly significant hits (E value < 1 x 1017) to known proteins (Table 2). bhoA is predicted to encode an IS3 family transposase with two out-of-frame overlapping ORFs (orfA and orfB), as previously described for IS3 transposases (39). However, a frameshift mutation detected in the 3' ORF may have inactivated the transposase. Southern blot hybridization with a bhoA-specific probe identified 40 to 50 copies in the genome of B. holmesii but no copies in the genome of B. pertussis (data not shown). bhoB and bhoC are homologous to adjacent genes from the Ralstonia metallidurans CH34 genome. These encode a probable extracytoplasmic solute receptor and a D-isomer-specific 2-hydroxyacid dehydrogenase. bhoD encodes a putative extracytoplasmic function sigma factor, and bhoE encodes a homolog of the FecR proteins (Table 2). The GC content of the insert is 62%, similar to the average GC content of B. holmesii genes (61.5 to 62.3%) (48) but lower than the shared component of the IUI (66.0%).
|
View this table: [in a new window] |
TABLE 2. Annotation of the B. holmesii-specific insertion in IUI
|
The IUI sequence was also compared to the B. avium 197N genome sequence (38). Using the BLAST programs, a syntenic region with homology extending beyond the ends of the IUI was identified (Fig. 2). The average nucleotide sequence identity of the B. avium region to the B. holmesii IUI was 73%, indicating that these sequences are not highly conserved across all members of the genus, but sequences near the ends of the IUI were more highly conserved (at least 88%) than sequences in the middle. Although gene order is otherwise largely conserved, the B. avium genome shows no detectable homology to the alcABCDERS, the fauA, the bhoABCDE, and the IS481 loci (Fig. 2; see also Table S1 in the supplemental material).
Alcaligin expression and detection. To determine whether the IUI alcaligin biosynthesis, export, and uptake locus is expressed in B. holmesii, the transcript levels of alcA and fauA were measured by quantitative reverse transcription-PCR of B0437 grown under iron-depleted and iron-replete conditions. The transcript abundances of alcA and fauA were 3.3- and 4.8-fold higher, respectively, in the iron-depleted sample than in the nondepleted sample, indicating that the locus is expressed and iron regulated in vitro.
LC-MS of culture supernatant from iron-starved B. holmesii detected a compound in the iron-depleted supernatant but not in the iron-replete supernatant with m/z 405 (M + H)+ (Fig. 3), corresponding to the mass of desferri-alcaligin (27). MS fragmentation analysis detected two fragment ions with m/z 283 and m/z 387, corresponding to previously described alcaligin fragment ions (27).
![]() View larger version (15K): [in a new window] |
FIG. 3. LC-MS spectra of B. holmesii supernatants from iron-depleted and iron-replete cultures. (A) LC-MS spectra of iron-depleted (Fe) and iron-replete (+Fe) B. holmesii supernatants, depicting ions with m/z 405 (corresponding to alcaligin) and with retention times between 0 and 45 min. (B) Fragment ions detected after collisional activation (35% energy) of the peak from (A), with a retention time of 36.99 min. After measurement of the reference standard, the calibration deviated by approximately 1.5 mass units, explaining the difference in m/z values compared to those of the fragment ions as reported in reference 22. The nomenclature of fragment ions is according to that in reference 27.
|
![]() View larger version (8K): [in a new window] |
FIG. 4. Neighbor-joining phylogeny of Bordetella spp. and related ß-proteobacteria based on four genes outside the rRNA operon. The unrooted tree is based on 3,559 fully informative nucleotides from an alignment of concatenated sequences of atpD, rpoB, tuf, and rnpB from seven B. holmesii strains, B. pertussis Tohama I, B. parapertussis 12822, B. bronchiseptica RB50, B. avium 197N, B. hinzii BC304, B. trematum DSM11334, and B. petrii SE1111R, and from A. xylosoxidans ATCC 15173 and B. pseudomallei K96243. Bootstrapping values greater than 50% (in 1,000 resamplings) are indicated at branches.
|
![]() View larger version (40K): [in a new window] |
FIG. 5. Southern blot hybridizations of B. pertussis, B. holmesii, and B. avium with 16S rRNA gene probes. Genomic DNA from each of the three strains was restriction digested with ClaI or NcoI as indicated and hybridized to a B. pertussis-specific probe (left panel) and a broad-range 16S rRNA gene probe (B-16S8F [right panel]) (15). Biotinylated DNA markers are in the first lane, and numbers at the left are fragment sizes (in base pairs).
|
Among the sequenced Bordetella genomes, the intergenic transcribed spacer region (ITS) between the 16S and 23S rRNA genes is more varied than the sequences of the structural rRNA genes (Table 3). However, ITS sequences from seven B. holmesii strains, all of which were identical to each other, were 99.6% identical to that from B. pertussis (GenBank accession numbers DQ420075 to DQ420095) and more similar than even the ITS sequence from B. bronchiseptica. A partial sequence of the 23S rRNA gene (B. pertussis bases 1043 to 1919) was also determined following broad-range PCR of three B. holmesii strains. Nine sequences (three from each strain; GenBank accession numbers DQ420059 to DQ420061) were identical to each other but only 97.2% identical to B. pertussis and more divergent than the B. avium sequence (Table 3). For comparison, B. pertussis and B. bronchiseptica vary at only a single nucleotide in this region. Therefore, the region of unexpectedly high identity between the B. pertussis and B. holmesii rRNA gene operons includes the 16S rRNA gene and the ITS but does not include the 3' end of the 23S rRNA gene. The boundary between nearly identical (>99.5%) and divergent sequences lies in the first kilobase of the 23S gene, which was not sequenced in this study.
|
View this table: [in a new window] |
TABLE 3. Percentages of nucleotide identity in the rRNA operon between B. holmesii and selected Bordetella species
|
|
|
|---|
Only 2 of 3,666 nucleotide positions in the housekeeping genes were polymorphic among seven geographically distinct B. holmesii strains. This level of sequence conservation is comparable to the three polymorphisms of 3,666 nucleotides observed between B. bronchiseptica and B. pertussis, which diverged and expanded clonally from the B. bronchiseptica lineage 2 to 5 million years ago (12, 30). Likewise, B. holmesii may have diverged and expanded within a similar time frame.
Comparative genomic hybridization identified a 66-kb DNA region, designated the IUI, which, unlike most of the B. holmesii genome, was highly conserved in B. holmesii and the classical bordetellae. No conserved chromosomal regions were identified when B. avium genomic DNA was hybridized to the classical Bordetella microarray (data not shown). Detailed molecular characterization verified that the shared sequences are more than 99% identical, much higher even than the expectation based on strongly conserved phylogenetic marker genes (under the assumption that the nucleotide substitution rates are similar across the genome). Cloning of the IUI ends also demonstrated sharp boundaries between the conserved IUI sequence and the weakly conserved flanking sequences. These data are consistent with the recent transfer of the IUI from a classical Bordetella genome to the B. holmesii genome. Transfer from a B. bronchiseptica or B. parapertussis genome cannot be ruled out, but the presence in the IUI of IS481 elements, which have been found in all B. pertussis and B. holmesii strains analyzed but in only two B. bronchiseptica strains (12, 45), argues strongly that a B. pertussis strain was the donor of this DNA. Although the IUI region is rearranged in B. pertussis Tohama I relative to this region in B. holmesii and the other mammalian bordetellae, another B. pertussis strain, 18323, had a genomic architecture and an IS481 element distribution very similar to those observed in the B. holmesii IUI, suggesting that the B. holmesii IUI may have been transferred from an 18323-like B. pertussis strain. B. pertussis 18323, isolated from a pertussis patient in 1947, is an atypical B. pertussis strain by many criteria, including pulsed-field gel electrophoresis (17), multilocus sequence typing (12), multilocus enzyme electrophoresis (45), and CGH (11) results, but a genetically related strain, CZ, was isolated from a pertussis case in 1993 (4, 17), suggesting that this is a rare, but circulating, variant.
The data presented here do not rule out the possibility that the IUI was transferred in the other direction, from B. holmesii to the classical bordetellae. However, our interpretation is supported by two lines of evidence. First, the GC content of the IUI is similar to that of the classical Bordetella chromosomes but significantly higher than that of B. holmesii protein-coding genes outside this region. Because atypical nucleotide composition is a hallmark of recently acquired DNA (reviewed in reference 13), this suggests that the IUI was recently transferred into B. holmesii. Second, the presence of multiple IS481 elements in the B. holmesii and B. pertussis IUIs, but not in the other classical Bordetella genomes, is more easily explained by our model, which does not require transfer of an IS481-containing DNA fragment to the ancestor of the classical Bordetella lineage, followed by precise excision of those transposons in two of three derived species. B. pertussis and B. holmesii have both been isolated from the respiratory tracts of humans, making the human airway a likely environment in which the transfer of the IUI could have occurred. At the resolution of CGH and PCR analysis, IUIs were identical in all examined B. holmesii isolates, with the exception of a variably present IS481 element, suggesting that the IUI was acquired recently or is under selective pressure to maintain its genomic organization.
A region of synteny between B. avium, B. holmesii, and the classical bordetellae extended through the IUI and into the flanking sequences on either side, indicating that the gene order of this chromosomal region is conserved across distantly related Bordetella species (Fig. 2). This finding, together with the failure to identify signatures of mobile elements (e.g., phage or conjugative transposon genes) at the ends of the IUI, suggests that the most likely mechanism for the proposed integration of the IUI into the B. holmesii genome is homologous recombination between the laterally transferred IUI and the ancestral B. holmesii genome, which may have diverged from the B. pertussis sequence to an extent similar to the extent of the divergence of the B. avium sequence. The degree of nucleotide sequence identity between the B. holmesii IUI and the orthologous sequences in B. avium is highest in the vicinities of the left and right insertion breakpoints, indicating that these genes may be more highly conserved across the Bordetella genus than the average gene. These sequences are proposed to have served as the substrates for the putative double-recombination event that replaced the B. holmesii sequence with a transferred B. pertussis fragment.
Given the evidence presented above for lateral transfer of the IUI from B. pertussis to B. holmesii, lateral transfer is also a plausible explanation for the unexpectedly high sequence similarity of their 16S rRNA genes. The failure to identify a 16S rRNA gene within, or in close proximity to, the IUI suggests that it was independently transferred from B. pertussis or that a subsequent rearrangement separated the 16S rRNA gene from the IUI. The proposed exchange includes the 16S rRNA gene and the ITS and may include a 5' fragment of the 23S rRNA gene.
Immediately following the proposed 16S rRNA gene transfer event, the recipient B. holmesii chromosome would have carried a single copy of the B. pertussis 16S rRNA gene and, unless there was only a single rRNA locus that was replaced by homologous recombination with the foreign DNA, one or more copies of its ancestral 16S rRNA gene. Other microbial genomes have been found to harbor one or more foreign copies of the small-subunit rRNA locus in addition to their presumed native copies (1, 37, 49). In order to detect putative ancestral B. holmesii 16S rRNA loci, independent hybridization- and PCR-based methods were used to search for copies of the 16S rRNA gene that diverged from the B. pertussis-like sequence. However, only B. pertussis-like 16S rRNA genes were found, suggesting that all copies of the ancestral recipient strain's 16S rRNA gene have been either lost or converted by recombination since acquisition of the B. pertussis-like 16S rRNA locus. The putative replacement of the ancestral B. holmesii 16S rRNA, which, by analogy to B. avium, might have been more than 98% identical to the B. pertussis sequence, may not have had serious functional consequences for the organism, because heterologous rRNA genes have been shown to function in some cases. For example, the Salmonella enterica serovar Typhimurium rRNA operon, although only 97% identical to that of Escherichia coli, can functionally replace the native E. coli loci (2).
The complement of genes carried by the IUI is of particular interest because it may include virulence factors or host specificity determinants that have contributed to the emergence of B. holmesii as a human pathogen. Of the 57 putative and confirmed genes (excluding transposons) carried on the IUI, the alcaligin locus, consisting of 8 genes, is the most recognizable candidate for a role in virulence. Alcaligin is an important iron-scavenging siderophore in B. pertussis and B. bronchiseptica (3, 6, 7, 19, 31) and is required for the maximal virulence of B. bronchiseptica in a piglet infection model (33). The demonstration of iron-regulated alcA and fauA transcription and the detection of alcaligin in the supernatants of iron-limited B. holmesii cultures indicated that the alcaligin biosynthesis and export locus was functional in B. holmesii. Although B. avium is capable of acquiring iron through other mechanisms (e.g., heme uptake [28]), alcaligin production has not been detected in B. avium, and its genome does not appear to encode an alcaligin biosynthesis locus, indicating that not all Bordetella species possess the ability to produce this siderophore. Therefore, prior to the IUI acquisition, the hypothetical progenitor of B. holmesii may not have been competent to produce alcaligin. Wholesale acquisition of this function by lateral transfer from B. pertussis could have provided B. holmesii with a new, highly efficient iron uptake system, potentially leading to an enhancement of its ability to colonize the human host, in which free iron is sequestered. Additionally, BhoD and BhoE, encoded by ORFs in the B. holmesii-specific IUI fragment, are homologous to E. coli FecI and FecR, which regulate ferric citrate transport (5, 14, 46, 47), and to related Bordetella heme uptake regulators (20, 32, 43, 44), raising the possibility that these may regulate iron uptake systems in B. holmesii.
Other genes in the IUI may also contribute to B. holmesii host specificity or virulence. For example, mar, adjacent to the alcaligin locus, encodes a member of the MarC family of multiple antibiotic resistance efflux pumps. Multidrug resistance transporters have been shown to confer resistance to specific host antimicrobial peptides in diverse bacterial pathogens (e.g., Neisseria meningitidis [40] and Staphylococcus aureus [23]). Perhaps the mar gene product enhances colonization of the host airway by B. holmesii by increasing resistance to host defensins and other endogenous antimicrobial peptides. Some ORFs in the IUI encode putative cell surface or exported proteins, including a putative lipoprotein (BP0793) and a putative exported protein (BP2486), that could influence pathogenesis by modifying interactions with the host epithelium and immune system. The presence of the B. pertussis Bvg-repressed gene vrg6 in the IUI is interesting, but this gene is not required for virulence in B. bronchiseptica (24), making its role in B. holmesii virulence unlikely. The numerous IUI genes encoding so-called housekeeping functions (e.g., metabolic enzymes, cytochromes, potassium transporters, cell division proteins, and molecular chaperones) might contribute to pathogenicity by enhancing fitness and so cannot be ruled out as virulence factors without further functional characterization. The potential roles of the hypothetical and conserved hypothetical ORFs, which account for about 20% of the genes in IUIs, remain to be explored.
Only recently has B. holmesii been identified as a human pathogen and recognized to cause pertussis-like disease. However, in spite of its similarities in host range and clinical course to those of B. pertussis, B. holmesii is not as closely related to B. pertussis as are the other mammalian bordetellae. The results presented here suggest that lateral transfer of a genomic island from B. pertussis may have contributed to the emergence of B. holmesii and its adaptation to the human host. These results represent a significant advance in the characterization of B. holmesii and in our understanding of the evolution of virulence in the important group of pathogens comprising the genus Bordetella.
This work was supported by a travel grant from The Netherlands Organization for Scientific Research (NWO). C.A.C. was supported by an American Lung Association research training fellowship. D.A.R. received grant support from the NIH (grants AI054970 and AI057188).
Published ahead of print on 13 October 2006. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
D.A.D. and C.A.C. contributed equally to this work. ![]()
Present address: Department of Microbiology and Immunology, University of Melbourne, Parkville, Victoria, Australia. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»