Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, United Kingdom,1 Department of Molecular and Cellular Biology, University of Guelph, Guelph, Ontario, N1G 2W1, Canada,2 Centre for Veterinary Science, Department of Veterinary Medicine, University of Cambridge, Madingley Road, Cambridge, CB3 0ES, United Kingdom,3 Department of Biophysics and Biochemistry, University of Rochester, Rochester, New York,4 The Witebsky Center for Microbial Pathogenesis and Immunology, Department of Microbiology and Immunology, The University of Buffalo, The State University of New York, Buffalo, New York 14214,5 Department of Population Health and Pathobiology, College of Veterinary Medicine, North Carolina State University, Raleigh, North Carolina 27612,6 Department of Biology, Drew University, Madison, New Jersey 07940,7 Department of Biology, MSC7801, James Madison University, Harrisonburg, Virginia 228018
Received 16 December 2005/ Accepted 2 June 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Of the other bordetellae described, B. hinzii is usually found as a commensal of birds but has occasionally been seen as an opportunistic pathogen of immunocompromised humans (21, 52), and B. holmesii has been associated with septicemia (114) and, more recently, with respiratory illness possibly resembling pertussis (67, 116), while B. trematum has been described as an occasional pathogen of humans causing wound infections and otitis media (108). B. petrii was first isolated from a bioreactor and is not thought to be a pathogen (113).
B. avium is the causative agent of bordetellosis in avian species, which is similar to respiratory disease caused by other bordetellae (48, 53). B. avium infects many different wild and domesticated birds (47, 49, 84), but commercially raised turkeys are particularly susceptible (94). Similar to B. bronchiseptica (75, 76), B. avium survives for long periods of time in water and dilute salt solutions, which may provide a means of survival between hosts (84). B. avium has a specific tropism for ciliated respiratory tissue, a property that it shares with B. pertussis, B. parapertussis, and B. bronchiseptica (7).
Among them, B. pertussis, B. parapertussis, and B. bronchiseptica share a large complement of virulence factors which can be broadly categorized as adhesins or toxins. The adhesins include filamentous hemagglutinin (FHA), pertactin (Prn), tracheal colonization factor, serum resistance protein (BrkA), and fimbriae. The toxins include pertussis toxin (PTX), adenylate cyclase-hemolysin (CyaA), dermonecrotic toxin (DNT), and tracheal cytotoxin (TCT) (65). Of these, only FHA, fimbriae, DNT, and TCT have been identified in B. avium (36, 95). The expression of these virulence-activated genes and another set of virulence-repressed genes is regulated by a two-component master regulatory system encoded by the bvgAS locus (reviewed in reference 24), also previously identified in B. avium (95).
In this report, we present the complete genome sequence of B. avium strain 197N and a comparative analysis between this and the genomic sequences of B. bronchiseptica, B. pertussis, and B. parapertussis (71). We also present some experimental data related to phenomena suggested by the bioinformatic analysis of the genomic sequence data. Comparison of the genomes of these four bordetellae will provide insights into the virulence factors of B. avium, provide information about the mechanisms underlying the host specificities and evolution of members of the Bordetella genus, and suggest many testable hypotheses about the pathogenesis, ecology, and biology of the bordetellae.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Sequence analysis and annotation. The sequence was assembled, finished, and annotated as described previously (71), using Artemis (88) to collate data and facilitate annotation. The DNA and predicted protein sequences of B. avium were compared to genomic sequences of the other three Bordetella species using the Artemis Comparison Tool (15). Orthologous gene sets were calculated by reciprocal best-match FASTA comparisons. Pseudogenes had one or more mutations that would prevent translation; each of the inactivating mutations was checked against the original sequencing data.
Polymorphism sequence analysis. For identifying polymorphisms in bvgS, the following primers were used: 5'GCACGCCCATCTATGCCATCAC3' and 5'TGCAGCAGATGTTCGGCATCC3'. Amplification of the 301-base-pair fragment spanning the poly(C) region was accomplished utilizing GC-cDNA polymerase from Clontech. Amplicons were sequenced by MWG Biotechnologies.
Construction of mutants. Internal fragments of the target genes (bfrH, bfrB, and bfeA) were amplified from B. avium genomic DNA using PCR and cloned into the HinDIII and BamHI sites of pFUS2 (6). Gene disruptions were performed by conjugation using B. avium strain KO1, a bhuR mutant of wild-type strain 4169rif (69), to generate KO1bfrH, KO1bfrB, and KO1bfeA. Transconjugants were selected on brain heart infusion (BHI) agar containing gentamicin.
Outer membrane protein isolation. A modification of the protocol of Leyh and Griffith (59) was used to purify B. avium outer membrane proteins (OMPs) from Fe-replete and Fe-stressed cultures and was described elsewhere (69). Fe-stressed conditions were achieved by adding enough ethylene-di-o-hydroxyphenylacetic acid (EDDHA) to reduce the cell density (measured by the optical density at 600 nm) of an 18-h culture (cultured at 37°C and 300 rpm) to 50 to 60% of the growth of duplicate culture grown under Fe-replete conditions. OMPs were harvested from B. avium cultured at 37°C to stationary phase in 200 ml of BHI (Difco Laboratories, Detroit, MI) broth containing either 36 µM FeSO4 (Fe replete) or 50 µM EDDHA (Sigma Biochemicals, St. Louis, Mo.) (Fe stressed). The total protein concentration of the OMP suspensions was determined with a Bio-Rad protein assay using bovine serum albumin as the standard. OMPs were resolved by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) (56) and visualized by staining with colloidal Coomassie brilliant blue (70).
LPS preparation and analysis. For lipopolysaccharide (LPS) analysis and immunological staining, bacteria were prepared as described previously (81). Briefly, approximately 1 x 109 bacteria were pelleted from a 1-ml suspension and resuspended in 500 µl of phosphate-buffered saline. Next, 250 µl of LPS buffer I (0.1875 M Tris-HCl, pH 6.8, 6% [wt/vol] glycerol) was added, and the mixture was heated in a boiling water bath for 5 min. Then, 10 µl of the cell lysate was added to 35 µl of LPS buffer II (0.0624 M Tris-HCl, pH 6.8, 0.1% SDS, 10% [wt/vol] glycerol, 0.1% [wt/vol] bromophenol blue) along with 10 µl of proteinase K (25 mg/ml). The sample was incubated at 55°C for 12 to 16 h. Samples were heated in a boiling water bath for 5 min prior to SDS-PAGE. Twenty microliters of sample was loaded onto SDS-PAGE gels. LPS samples were analyzed by SDS-PAGE in the Tricine-buffered gel system originally described by Schagger and von Jagow (89), with the modifications described by Lesse and colleagues (58). Silver staining was performed according to the method of Tsai and Frasch (105). Western blot analysis was performed as previously described (3) using the monoclonal antibodies BL-2 (band A LPS specific) and BL-8 (band B LPS specific) (64).
Nucleotide sequence accession number. The genome sequence has been submitted to EMBL with the accession number AM167904.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
To distinguish CDSs orthologous between B. avium and B. bronchiseptica or unique to each (which may include paralogous genes), reciprocal FASTA analysis followed by manual curation was performed. B. avium and B. bronchiseptica share only 2,380 orthologous CDSs (Fig. 3). These shared CDSs encode primarily core/essential functions, whereas the species-specific CDSs are involved mainly in accessory functions (Fig. 4), which is likely to reflect the difference in lifestyle and pathogenicity of these two organisms (for example, infection of avian versus mammalian respiratory epithelia). The genome of B. avium contains 1,109 CDSs (31.5%) that do not have orthologues in B. bronchiseptica (Fig. 1 and 3). B. bronchiseptica, with nearly 46% more coding regions than B. avium, has 2,628 CDSs (52.5%) that are absent from B. avium (Fig. 3). These figures are larger than the variations observed within other species (e.g., Escherichia coli O157:H7 contains 1,387 genes [26%] not present in E. coli K-12, which has 528 [12%] unique genes [73]) and roughly comparable to those identified from other congeneric species comparisons (e.g., Bacteroides fragilis versus Bacteroides thetaiotaomicron, where 1,941 [44.3%] and 2,337 [48.9%] of the genes are unique, respectively [16]).
|
|
There are 68 pseudogenes in the B. avium genome, more than in B. bronchiseptica (18) but substantially fewer than in B. parapertussis (220) or B. pertussis (358). B. pertussis and B. parapertussis are clonal and of recent descent. The presence of relatively few pseudogenes suggests that B. avium has not undergone a recent evolutionary bottleneck and may be older than B. pertussis and B. parapertussis. This correlates with the fact that B. avium more closely resembles B. bronchiseptica in host range and ability to survive in multiple niches (84).
Host interaction/virulence. The overall comparison between B. avium and B. bronchiseptica shows a very large degree of variation in genes for surface structures, both proteins and polysaccharides. It is likely that these differences are central to their different specificities for avian and mammalian tracheae. These components include agglutinins/adhesins, LPSs, capsules and extracellular polymers, fimbriae and pili, autotransporters, large secreted proteins, secretion systems, and toxins and are discussed in detail in the following sections.
Agglutinins/adhesins. FHA, encoded by the fhaB gene, is the major adhesin and hemagglutinin in B. pertussis, B. parapertussis, and B. bronchiseptica (54, 87). In B. bronchiseptica and B. pertussis, FHA is produced in the form of a preprotein, FhaB, which is exported to the cell surface via a two-partner secretion system, mediated by the secretory protein FhaC (50). The C-terminal region of FhaB is removed by a specific, anchored, subtilisin autotransporter, SphB1, to yield the mature protein FHA, which is partitioned between the cell surface and the extracellular milieu (26). B. avium has orthologues of fhaB and fhaC (BAV1966 and BAV1961, respectively) in exactly the same arrangement as B. bronchiseptica, with fhaC located as the fifth gene in the fimABCD operon (see "Fimbriae" below). FhaB of B. avium exhibits a low level (29% identity) of sequence similarity to and is shorter than (2,621 amino acids [aa]) that of B. bronchiseptica (3,652 aa). Furthermore, B. avium fhaB and fhaC insertion mutants retain hemagglutination ability but are avirulent in the turkey poult model (95). B. bronchiseptica has two additional FHA-like proteins, FhaL and FhaS, orthologues of which are also encoded in the B. avium genome (BAV2159 and BAV1551). In addition, the genome of B. avium contains six novel genes encoding FHA-like proteins (BAV0121, BAV0712, BAV0937, BAV1217, BAV2819, and BAV2825), of which two (BAV0937 and BAV1217) are pseudogenes.
Among the novel genes are two, BAV2824 (hagA) and BAV2825 (hagB), that encode the hemagglutination phenotype of B. avium (102); a mutation in either gene abolishes hemagglutination (unpublished data). The paired loci are similar to the fhaB/fhaC pair of B. avium and the other sequenced bordetellae, suggesting that they may also form a two-partner secretion system. This gene pair is physically located 4 CDSs from a very similar pair of genes, BAV2818 and BAV2819, with predicted proteins of 96% (hagB and BAV2819) and 53% (hagA and BAV2818) identity to each other at the amino acid level. Apart from the extreme N and C termini, the FhaB-like predicted proteins from hagB and BAV2819 are identical at the amino acid level across most of their lengths. While hagA and hagB are divergently transcribed, BAV2818 and BAV2819 occur in tandem and may constitute an operon. Work to determine the role, if any, of the second set of genes in hemagglutination and virulence is ongoing.
Two additional CDSs unique to B. avium, BAV1656 and BAV0856, encode adhesin-like proteins. The protein encoded by BAV1656 is 29% identical to E. coli Tia, an outer membrane protein involved in adherence and invasion of epithelial cells (35); it is also similar to Hek, a hemagglutinin encoded in a pathogenicity island of uropathogenic E. coli (32). BAV0856 encodes a potential adhesin that is 31% identical to Bap, a Staphylococcus aureus surface protein involved in biofilm formation (28), and 24% identical to E. coli intimin (51).
LPS. LPS plays an important role in the virulence of the bordetellae (80). In B. avium, it is required for colonization of turkeys, serum resistance, tracheal ring binding, and resistance to antimicrobial peptide-mediated killing (96). It also serves as the attachment site for the B. avium bacteriophage Ba1 (92, 96). There is very little information regarding the structure of B. avium LPS except for the characterization of the O antigen domain (see below). Thus, the genome sequence provides the first clues as to the relatedness of the LPS of B. avium and the previously sequenced bordetellae. B. avium contains orthologues of many of the LPS loci characterized in B. bronchiseptica (Table 2), which suggests that B. avium LPS is structurally similar to B. bronchiseptica LPS, although it should be noted that single amino acid differences between some lipid A biosynthesis genes can result in an altered substrate specificity that can affect the precise lipid A structure and have profound effects on endotoxin activity, for example (100). Structural characterization of B. avium LPS to address issues such as these is under way.
|
The O antigen is the distal domain of the LPS molecule relative to the bacterial cell surface. In the B. avium type strain, it has been determined as a polymer of 1-4-linked 2-acetamidino-3-[3-hydroxybutanamido]-2,3-dideoxy-ß-D-glucopyranosyluronic acid (57). The O antigens of B. bronchiseptica and B. parapertussis are a polymer of 2,3-diacetamido-2,3-dideoxy-
-L-galactopyranosyluronic acid (2,3diNAcGalA) (31), and in B. hinzii the O antigen is composed of trisaccharide repeating units of both the gluco and galacto isomers of the same uronic acid (8, 111). Thus, all of these O antigens are comprised of similar 2,3-dideoxy-hexuronic acids. This suggests that the genetics of O antigen expression in these bordetellae are similar. In B. bronchiseptica and B. parapertussis, the locus responsible for O antigen expression, wbm, contains 24 genes, including those encoding enzymes that are presumably involved in modification of the O antigen sugars (e.g., epimerase/dehyratases) and an ABC transport system similar to other O antigen export systems (79). At an equivalent location in B. avium, a putative polysaccharide biosynthesis locus is present. However, despite the structural similarity between the B. avium and B. bronchiseptica/B. parapertussis O antigens, the B. avium locus displays no similarities to wbm. It comprises only nine genes, including those encoding an O antigen translocase (BAV0081), a possible polysaccharide repeat unit exporter/flippase (BAV0084), a glycosyl transferase (BAV0086), an amidotransferase (BAV0087), and an acyltransferase (BAV0089), but no obvious transport system. While it remains to be confirmed that this locus is involved in B. avium O antigen expression, the difference between this locus and B. bronchiseptica/B. parapertussis wbm may suggest that it is involved in modifications of 2,3diNAcGlcA (for example, the unusual 3-hydroxybutanamido and the 2-acetamidino substitutions) and that the genes responsible for expression of the O antigen polymer backbone are elsewhere in the chromosome. However, we have not identified any obvious wbm homologues in the B. avium genome, and thus it must be concluded that B. avium utilizes a novel biosynthesis pathway for expression of its O antigen.
In B. pertussis, B. parapertussis, and B. bronchiseptica, the wlb gene cluster (notated wlbA to wlbL), which is immediately adjacent to the wbm locus, is involved in the addition of a trisaccharide to the distal portion of the LPS core (4), generating a structure known as band A LPS. B. avium also possesses a wlb locus at an equivalent location that is highly similar to that of B. bronchiseptica, which suggests that the two loci perform the same function; thus, it is suggested that B. avium expresses an LPS containing a band A trisaccharide structure. However, there are two differences between the B. avium and B. bronchiseptica loci, in that wlbD and wlbI are missing from B. avium. The function of WlbI is unknown, but WlbD is an epimerase that is required for expression of full-length LPS. The phenotype of a B. pertussis wlbD mutant suggests that WlbD is involved in the synthesis of the middle sugar of the trisaccharide 2,3diNAcManA (83). The absence of wlbD in B. avium suggests that any band A-like trisaccharide expressed by B. avium would comprise the gluco rather than the manno epimer of the di-N-acetylated uronic acid. However, LPS from B. avium is antigenically related to that of B. bronchiseptica, as shown by Western analysis using anti-LPS antibodies (Fig. 5). Antibodies against band A (the core plus trisaccharide) reacted with crude LPS preparations from B. avium strain 197N. In these experiments, LPS from B. avium wlbA and wlbL mutants did not react with anti-band A antibodies. Further, B. avium wlbA and wlbL mutations were complemented with the wlb region from B. bronchiseptica, as shown by normal LPS synthesis and phage Ba1 sensitivity (96).
|
Capsular polysaccharide. B. avium is known to produce a capsule (53), and a putative capsular polysaccharide biosynthesis locus is present in the B. avium genome (BAV2635-BAV2653). This large locus contains several genes that are predicted to encode putative sugar modification functions, including a homologue of GDP-mannose pyrophosphorylase/phosphomannose isomerase, four glycosyltransferases, an oxidoreductase, and an epimerase/dehydratase. The clustering of these genes suggests that the locus is involved in polysaccharide biosynthesis, but bioinformatic analysis does not suggest a structure. The remainder of the locus comprises genes of unknown function. There is no obvious polysaccharide export/transport apparatus encoded by this locus. This locus is different from that in the B. bronchiseptica genome (BB2918-BB2934), suggesting that B. avium and B. bronchiseptica produce different types of capsular polysaccharides. In contrast, B. pertussis and B. parapertussis appear to be devoid of capsular polysaccharide, an observation supported by the presence of pseudogenes in their capsular polysaccharide loci.
The B. avium capsular polysaccharide locus is directly adjacent to a 12-gene operon, BAV2623-BAV2634, which carries CDSs similar in both sequence and organization to those required for the biosynthesis of a cellulose-like polymer (wssA-wssJ) and which contributes to biofilm formation, niche specialization, and fitness of Pseudomonas fluorescens (97). The presence of a cellulose-like polymer has not been reported for B. avium, but if it exists, it may well suggest that B. avium is able to form biofilms.
Fimbriae. In B. pertussis and B. bronchiseptica, fimbriae bind to sulfated sugars that are present on cell surfaces of the respiratory tract and therefore play an important role in tracheal colonization and development of a humoral immune response (38, 66). There are 11 putative fimbrial subunit genes in the B. avium genome (BAV0313, BAV0314, BAV0315, BAV0316, BAV0321, BAV0465, BAV1776, BAV1777, BAV1965, BAV2661, and BAV3039), whereas B. bronchiseptica, B. pertussis, and B. parapertussis have 6 intact genes, 3 intact genes and 1 pseudogene, and 4 intact genes and 2 pseudogenes, respectively (71).
Among the B. avium fimbrial genes, only one, BAV1965, is a true orthologue of a fimbrial subunit gene, fimA found in B. bronchiseptica, B. pertussis, and B. parapertussis. All the others are different in terms of both sequence similarity and chromosomal location.
In B. avium and the previously sequenced bordetellae, the fimA gene is part of an operon (fimABCD fhaC) encoding a chaperone-usher fimbrial biosynthesis system, which consists of a fimbrial subunit (FimA/BAV1965), a chaperone (FimB/BAV1964), a fimbrial usher (FimC/BAV1963), and a fimbrial adhesin (FimD/BAV1962). The location of fhaC in this fimbrial biogenesis operon is identical to the arrangement in B. pertussis and B. bronchiseptica (see "Agglutinins/adhesins" above). The fimABCD fhaC operon is functional and essential for virulence, since transposon mutants in fhaC and fimC were identified by signature tag mutagenesis (STM) using infection of turkey poults as a screening model (95). A second locus (BAV1773-BAV1777), unique to B. avium, carries two adjacent fimbrial subunits (BAV1776 and BAV1777), a chaperone (BAV1775), a fimbrial usher (BAV1774), and a fimbrial adhesin (BAV1773). Four other putative B. avium fimbrial subunit genes, BAV0313, BAV0314, BAV0315, and BAV0316, are contiguous and form a fimbrial gene cluster. Five of the B. avium fimbrial genes (BAV0313, BAV0315, BAV0316, BAV0321, and BAV2661) contain between six and nine repeats of a 7-bp sequence in their 5' ends. Slipped-strand mispairing within these repeats may provide a mechanism for phase variation among the different fimbrial structural proteins. It is interesting that under conditions of normal laboratory growth in rich medium, such as BHI broth, negatively stained cells examined by transmission electron microscopy exhibit few if any visible fimbriae (unpublished observations). This may mean that the genes encoding the putative structures are expressed under conditions that we have not yet tested, such as within the host. Work to determine which, if any, of these many potential fimbrial genes are expressed and under what conditions is ongoing.
Type IV pilus. The type IV pilus is involved in adherence and twitching motility and is important for virulence in some microorganisms (reviewed in reference 27). For example, Neisseria spp. require these pili for initial attachment to human epithelial cells, and type IV PAK pili mediate adhesion of Pseudomonas aeruginosa to mucosal surfaces. There are three potential type IV pilus gene clusters in the B. avium genome, two that appear intact (BAV2283-BAV2291 and BAV2357-BAV2364) and one that is probably nonfunctional due to the presence of several pseudogenes (BAV2525A-BAV2534). The same three loci are present in B. bronchiseptica, B. pertussis, and B. parapertussis, but all those in B. pertussis contain pseudogenes (71). Moreover, B. bronchiseptica has a fourth type IV pilus gene cluster (BB0776-BB0792) that is absent in B. pertussis, B. parapertussis, and B. avium.
Autotransporters. Autotransporters are modular proteins consisting of a conserved C-terminal pore-forming domain through which the N-terminal passenger domain is exported (46). The variable passenger domains have diverse functions. Several autotransporters (BapA, BapB, BapC, BipA, BrkA, Prn, TcfA, and Vag8) from the previously sequenced bordetellae are known to be important for host interaction and virulence (46).
B. avium has seven intact autotransporter CDSs (BAV0466, BAV0811, BAV1572, BAV1641, BAV1864, BAV2846, BAV3269), while BAV1640 is a pseudogene (containing a stop codon, a frameshift, and a premature truncation) and BAV2452 is potentially phase variable, as it contains two frameshifts, one at a nine-guanine [G(9)] tract and one at a G(10). This tally compares to 21 intact autotransporters in B. bronchiseptica, of which 2 are potentially phase variable (71). Although some of the B. avium autotransporters share similarity with those present in the other bordetellae (for example, BipA; see below), they are not true orthologues. This observation suggests that these autotransporters were independently acquired by ancestral strains of B. avium and the other sequenced bordetellae or that they may be different samples from a common gene pool. One autotransporter gene, BAV2846, was previously identified by virtue of its homology to an adhesin, AIDA-I from diarrheagenic E. coli, which is involved in host cell attachment, autoaggregation, and biofilm formation (93). This gene was subsequently cloned, sequenced, expressed in E. coli, and shown to be expressed in B. avium by reverse transcription-PCR (unpublished data).
Large novel surface proteins. The adjacent CDSs BAV1944 (4,342 aa, 447 kDa) and BAV1945 (6,460 aa, 650 kDa) are the largest in the genome. They are flanked by CDSs encoding components of a type I secretion system and potentially encode cell surface or exported proteins. Both proteins contain numerous locally repeated sequences. BAV1944 encodes a glycine- and serine-rich protein that is weakly similar to several proteins from a wide range of bacteria. A highly repetitive region occurs in the C-terminal half, with 10 occurrences of a 44-aa repeat sequence that are identical at 39 residues. These motifs are similar to hemolysin-type calcium-binding repeats (PF00353) that are conserved in proteins involved in host interactions (33). BAV1945 is most similar to a predicted bacterial immunoglobulin-like domain protein from Yersinia pseudotuberculosis (18) and to a putative cell wall-associated protein, srpA, of Streptococcus cristatus (44). These two predicted proteins, BAV1944 and BAV1945, are also weakly similar to each other and to other predicted autotransporters in the B. avium genome.
Motility. B. bronchiseptica is motile, while B. pertussis and B. parapertussis are not; this can be ascribed to the fact that several genes encoding motility- or chemotaxis-related functions are inactivated or deleted in B. pertussis and B. parapertussis (71). B. avium appears to have a complete and intact set of motility and chemotaxis genes. With regard to its role in virulence, motility appears to be completely shut off above 40°C, approaching the body temperature of turkeys (102). Kersters et al. (53) previously reported a significantly greater number of motile rods in cultures grown at room temperature than in those grown at 35°C. Motility in B. avium has an effect on bird-to-bird transmission but has no effect on the 50% infective dose in turkeys (102). This suggests a role in environmental spread that may also be the case in B. bronchiseptica (1, 2). A B. avium bvgS mutant (95) (see below) exhibits motility identical to the wild type at 35°C and 42°C (unpublished data); thus, unlike B. bronchiseptica (2), motility of B. avium may not be regulated by the BvgAS two-component system.
Toxins. As expected from previous reports (36), the two gene clusters encoding components of the pertussis toxin (PTX) and adenylate cyclase (cya) and their respective type IV and type I secretion systems are not present in the B. avium genome. Adenylate cyclase is produced by all species of Bordetella that are pathogens of mammals. On the other hand, PTX is produced by the human pathogen B. pertussis but not B. bronchiseptica and B. parapertussis, in spite of the presence of all the genes encoding PTX components. Thus, it appears that these two toxins are not required for infection of birds, and hence, upon its adaptation to these hosts, B. avium may have lost the genes encoding PTX and Cya. Alternatively, these toxin genes could have been acquired by the common ancestor of the mammalian species of Bordetella since diverging from the ancestor of B. avium.
The gene encoding another Bordetella toxin, dermonecrotic toxin (dnt), is present in the B. avium genome but at a location different from those of the previously sequenced species. The genes flanking the dnt gene in B. avium and the sequenced bordetellae are different. Although they are 41% identical, they may have been independently acquired and hence not true orthologues. DNT production is controlled by the presence of nicotinic acid or magnesium sulfate (36), implicating the bvgAS system (see below) in its expression control; however, bvgA and bvgS insertion mutants still produce DNT (unpublished observations). In contrast, a transposon insertion that abolished the production of DNT in strain WBA16 (102) was mapped in BAV0989, which encodes a LysR-family regulator (unpublished data). There are 54 LysR-type transcriptional regulators in the B. avium genome, many of which have orthologues in the previously sequenced bordetellae; interestingly, BAV0989 does not have an orthologue in those species. The DNT-negative strain is avirulent in the turkey poult model (102). These data indicate that global regulation of virulence genes in B. avium may be different from the other bordetellae.
Tracheal cytotoxin (TCT) is a glycopeptide cell wall fragment that is released in large quantities by B. pertussis in particular, to a lesser extent by B. bronchiseptica and B. parapertussis, and in small but unquantified amounts by B. avium (36). This molecule is cytotoxic in cell cultures via induction of nitric oxide biosynthesis pathways (39, 45). AmpG permease (ampG) is required in E. coli (17) and presumably other gram-negative bacteria for recycling of these cell wall anhydromuropeptides; this gene is present in the previously sequenced bordetellae. A gene similar to ampG occurs in B. avium, but in B. pertussis, this gene has an insertion sequence immediately upstream that may affect transcription (71), offering a potential explanation for the copious quantities of TCT produced in culture supernatants of B. pertussis. Whether TCT released by B. avium is significant in the pathogenesis of bordetellosis is unknown.
BvgAS system. The two-component regulatory system BvgAS is the master expression regulator of many virulence factors and is the basis of phase variation in Bordetella mammalian pathogens (25). In response to environmental signals, the inner membrane sensor histidine kinase, BvgS, initiates a complex phosphorelay involving three of its domains, the transmitter, receiver, and C terminus, before phosphorylating the response regulator, BvgA, which then acts as a transcriptional regulator of a number of genes. Translational frameshifting at a poly(C) site near the 3' end of bvgS, between codons for two of the phosphorelay domains, is responsible for phase variation via slip-strand mispairing in B. pertussis (98). Phase variation has been reported for B. avium strain 197 (37), its derivative strain 197N (the sequenced strain), and the type strain, ATCC 25086 (unpublished data). Interestingly, the B. avium phase variants have been reported to revert to the wild-type phenotype upon passage through the natural host, turkey (37). The molecular basis for phase variation in B. avium is unknown.
Unlike in the previously reported Bordetella genomes, the bvgAS genes in B. avium are not linked to the fhaB-fimABCD-fhaC gene cluster (this report and reference 95). Compared to those of B. bronchiseptica, the B. avium bvgAS genes vary in similarity in different regions of the predicted proteins (95). However, in the predicted BvgA protein, the phosphorylation site (D54) and helix-turn-helix motif are both retained, and in the predicted B. avium BvgS protein, the conserved phosphorelay sites in the putative transmitter (H729), receiver (D1023), and C terminus (H1172) domains (25) are also conserved.
The B. avium bvgS gene in the sequenced strain 197N has a frameshift within a homopolymeric tract of 8 cytosine residues in a location similar to the tract in B. pertussis described previously (98). This tract occurs between the predicted second and third phosphorylation sites described above. The C(8) tract appears in all clones of this region that were sequenced, and it predicts an earlier termination of the protein compared to other sequenced bordetellae. The presence of this bvgS variant in virulent strain 197N might lead to the conclusion that the B. avium bvgS gene is a pseudogene and that the bvgAS system is not involved in virulence in B. avium. However, Spears et al. (95) showed that bvgS was essential for virulence by using STM of strain 197N. To address this conundrum, we analyzed the poly(C) tract from chromosomal amplicons or cloned fragments and made the following observations (Table 3). (i) All PCR products from a number of different B. avium strain 197N derivatives that had been constructed and stored over the last 10 years had identical C(8) tracts, as in the sequenced strain. This group included the STM-derived strain PAS356, a nonvirulent mutant with a transposon insertion in bvgS (95), thus indicating that the shorter predicted BvgS protein in strain 197N is functional. To test the possibility that in vivo growth could select or enrich for bvgS mutants with an altered poly(C) tract (e.g., of 7 or 4 residues), we inoculated 1-week-old turkey poults with 197N (102) and sequenced bacterial DNA directly from tracheal swabs of infected turkeys taken 2 weeks after inoculation. We found that five of five amplicons tested contained C(8) tracts. Clearly, there was no in vivo selection for a frameshift in this tract of Cs, at least under the conditions we employed. However, a bvgAS clone from a chromosomal library prepared from strain 197N in 1996 was sequenced and had a C(7) tract and thus a predicted full-length protein product (95). (ii) A PCR amplicon of the bvgS gene from an older strain, WBA16 (102), constructed shortly after 197N was derived from 197 about 12 years ago, had a C(4) tract predicting the longer BvgS. (iii) Similarly, eight different B. avium isolates other than 197N derivatives, including ATCC type strain 35086, isolated at different times in different geographical regions all had C(4) tracts. It is important to note that we have not tested the virulence of those strains with C(4) or C(7) tracts under controlled in vivo conditions that we routinely employ.
|
In addition to the activation of a set of genes, the BvgAS two-component system of B. pertussis, B. parapertussis, and B. bronchiseptica represses the expression of another set of genes that are not expressed in the virulent state: virulence-repressed genes. This repression occurs via the activation of expression of a negative regulator, BvgR (68). In the previously sequenced bordetellae, the bvgR gene is located downstream of the bvgAS locus but is transcribed in the opposite orientation from it. In B. avium, no orthologue of bvgR is present. BvgR contains the EAL domain (PF00563) that is found in diverse bacterial signaling proteins. There are 10 CDSs in the B. avium genome that have this domain. One in particular, BAV1874, is immediately upstream of the bvgA orthologue. It is possible that this B. avium CDS encodes a regulatory protein that performs a function similar to BvgR of the other sequenced bordetellae.
Another group of genes expressed in the "intermediate" phase is best characterized for B. bronchiseptica and may represent a state essential for transmission between hosts (24). The best-studied member of this group is BipA, a cell surface protein with repetitive regions (99). The B. avium CDS, BAV0900, encodes a protein similar to BipA in B. bronchiseptica, but it differs from the latter in the number of internal repeats and does not appear to be a true bipA orthologue.
Secretion systems. Secretion of proteins into the extracellular environment or target cells is required for various aspects of the lifestyle of pathogenic bacteria. In gram-negative bacteria, at least six major transport machineries have been identified (type I to VI; for a recent review, see reference 55). All of these systems except type II are present in the previously sequenced bordetellae (71). In contrast, B. avium possesses all of the systems except type III and IV, which are important to virulence of other bordetellae (type III for B. bronchiseptica [117] and type IV for B. pertussis [19, 85]). Unique in B. avium compared to the previously sequenced bordetellae is a gene cluster that potentially encodes a type II secretion system (BAV0331-BAV0345). However, the proteins that may be exported through the type II secretion system in B. avium are unknown. Also, the putative type I system genes are located next to large genes encoding proteins similar to virulence factors in other bacteria mentioned above.
Iron acquisition. In animals, iron (Fe), an essential nutrient, is sequestered in high-affinity Fe-binding complexes, such as transferrin, lactoferrin, and hemoproteins (5), leaving the estimated free Fe concentration in biological fluids as low as 1018 M (43, 115). To survive within the Fe-limited environment of the host, bacterial pathogens have evolved different mechanisms to acquire Fe (90), including contact-dependent mechanisms, which involve specific outer membrane (OM) receptors that directly bind the host Fe source (22, 23, 41), and contact-independent mechanisms, which involve specific OM receptors that bind to and transport Fe-bound siderophores. Energy required to transport host-derived Fe or heme across the OM is provided by the proton motive force of the cytoplasmic membrane via the TonB-ExbB-ExbD proteins (13). TonB has been shown to be essential for virulence in diverse organisms, including B. pertussis (77, 86, 101, 104).
The genome of B. bronchiseptica encodes 16 putative (71) or experimentally described (9-11, 14, 78, 109) TonB-dependent ferric complex receptors (42). In comparison, analysis of the B. avium genome revealed only eight genes encoding putative or described (69) TonB-dependent ferric complex receptors, which are listed in Table 4.
|
|
Prophages and integrated plasmid. There are three prophages in the B. avium genome compared to four in the sequenced strain of B. bronchiseptica.
The B. avium prophage A (BAV
-A; BAV0391-BAV0443) has an internal region (BAV0416-BAV0430) that is highly similar to an internal region (Bbp12-Bbp27) of the B. bronchiseptica bacteriophage BPP1 that undergoes tropism switching (62). However, the B. avium prophage lacks the major tropism determinant and the reverse transcriptase responsible for the tropism switching.
Adjacent to BAV
-A there is a plasmid-like conjugative element (BAV0364-BAV0389). A similar element is inserted at the same location in the B. bronchiseptica genome. Both elements share highly similar conjugal transfer genes (vir/tra), but the B. avium version carries three genes encoding components of a type I restriction-modification system.
The B. avium prophage B (BAV
-B; BAV1280-BAV1342) is inserted at the 3' end of the stable RNA gene ssrA. In the B. bronchiseptica genome, a different prophage occupies the same attachment site. There is evidence for a complete life cycle of the B. avium phage B due to the presence of phage-specific DNA in culture supernatants (data not shown). This analysis was done on DNA prepared after supernatants were treated with DNase, indicating that the DNA was protected by a capsid; however, electron microscopy on the supernatants has not been performed, and to date there has not been a sensitive strain identified for this phage. Phage Ba1 (91), which infects strain 197N, is >80% identical at the DNA level over two-thirds of its genome with prophage B (data not shown); Ba1 integrates at a different chromosomal site and confers a distinctly different immunity to Ba1 lysogens (91). Ba1 lysogens are still virulent, but 50% of spontaneous Ba1-resistant mutants are attenuated, and all resistant mutants have noticeable differences in LPS structure or amount (91).
The third B. avium prophage (BAV
-C; BAV1423-BAV1482) also shares an attachment site (tRNAleu) with another prophage in the B. bronchiseptica genome. These attachment sites therefore appear to represent preferred locations for the integration of mobile genetic elements in the bordetellae genomes.
Nutrition and metabolism. B. avium, like the other sequenced bordetellae, cannot utilize glucose as a sole carbon or energy source (53). The genome reveals that the glycolytic pathway is incomplete due to the absence of genes encoding glucokinase and phosphofructokinase, similar to B. pertussis, B. parapertussis, and B. bronchiseptica. B. avium possesses genes that indicate a fully functional tricarboxylic acid cycle. Synthetic pathways for glutamate and aspartate appear to be complete, indicating that these amino acids are added to minimal medium for B. avium simply as sources of carbon and not to satisfy auxotrophic requirements (60). A nutritional requirement for cystine and methionine (53) is explained by the absence of a seven-gene cluster, present in other sequenced bordetellae, encoding components of the sulfate transport system (Sbp and CysUWA) and sulfate assimilation (CysHDN). The BAV3142-BAV3146 genes encode the components of a formate dehydrogenase, which requires selenocysteine. They are adjacent to CDSs encoding the components of the selenocysteine incorporation pathway; selenocysteinyl-tRNA(Sec) synthase (BAV3147; selA), selenocysteine-specific elongation factor (BAV3148; selB), selenide, water dikinase (BAV3149; selD), and the selenocysteil-tRNA (selC) itself. Both the formate dehydrogenase and the selenocysteine incorporation genes (except selD) are unique to B. avium. Other than these, there are few differences in the core metabolic genes, which probably reflects the similar nutritional environments of the niche inhabited by these organisms.
| Conclusion. |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
This work was supported by a grant to L. M. Temple and J. Parkhill from the USDA (539235), a grant to L. M. Temple and D. M. Miyamoto from NIH-AREA (526112), an NIH training grant (T32 DE007034) to N. D. King, and a grant to P. E. Orndorff and L. M. Temple from the USDA (2002-35204-12236).
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|