Previous Article | Next Article ![]()
Infection and Immunity, May 2005, p. 3128-3136, Vol. 73, No. 5
0019-9567/05/$08.00+0 doi:10.1128/IAI.73.5.3128-3136.2005
Steven L. Sviat,1 and
Michael Frace2
Bacterial Zoonoses Branch, Division of Vector-Borne Infectious Diseases, National Center for Infectious Diseases, Centers for Disease Control and Prevention, Fort Collins, Colorado,1 Genome Sequencing Laboratory Biotechnology Core Facility, National Center for Infectious Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia2
Received 27 December 2004/ Accepted 11 January 2005
|
|
|---|
|
|
|---|
Bartonella infections are associated with arthropod vector transmission; e.g., B. bacilliformis is transmitted by sand flies, B. henselae by cat fleas, and B. quintana by the human body louse. There is recent evidence of Bartonella association with ticks, fleas, and flies (10, 14, 15, 46). Furthermore, several distinct Bartonella strains have been isolated from various rodent (5, 7, 8, 22, 23, 29) and ruminant (19, 6, 13) species throughout the world. Because of the ubiquitous nature of Bartonella and their association with arthropod vectors, this genus has been implicated as a potential causative agent in diverse disease manifestations seen in humans for which there have been no definitive clinical diagnoses. Our laboratory was interested in exploring whether the prevalence of Bartonella residing in wild-mammal populations may have had an association with unexplained human febrile illnesses in the North American Southwest, as there was some serological evidence suggesting such a link (M. Y. Kosoy et al., Abstr. Am. Soc. Rickettsiol. Bartonella Emerg. Pathog. Group 2001 Joint Conf., abstr. 108, 2001; F. Koster et al., Abstr. Am. Soc. Rickettsiol. Bartonella Emerg. Pathog. Group 2001 Joint Conf., abstr. 133, 2001). Moreover, a Bartonella (B. washoensis) strain isolated from a Nevada patient with myocarditis was shown to be identical to strains cultured from ground squirrels in the same area (27). We wished to further investigate the putative role of rodent-associated Bartonella in human infections and, as an initial step, to identify Bartonella species-specific antigens that elicit antibody responses during infection. The goal was to initiate a primary analysis of immunodominant antigens related to rodent-borne Bartonella and to further develop a more detailed comparative analysis of immunogenic proteins from human Bartonella pathogens to define genus-specific antigens versus species-specific antigens. As a result of these aims, this report describes a novel immunodominant, surface-associated protein from B. vinsonii subsp. arupensis encoded by one of the largest bacterial genes yet identified. Additionally, at least two more related genes are present upstream, forming a multigene family. The genes in this family are characterized by multiple internal repetitive regions that are conserved within and among the individual genes, which we have therefore termed brp (bartonella repeat protein).
Identification of the brp gene complex. Mice were infected with live Bartonella strains isolated from rodent species to generate antibodies that could be used to clone and characterize genes encoding proteins associated with infection. Strains na19103nm (isolated from Neotoma albigula [wood rat] in New Mexico), Pm15590co and Pm136co (isolated from Peromyscus maniculatus [deer mouse] in Colorado), and Sb944nv (isolated from Spermophilus beecheyi [California ground squirrel] in Nevada) were obtained from Michael Kosoy, Division of Vector-Borne Infectious Diseases, Centers for Disease Control and Prevention, Fort Collins, Colorado. B. henselae, B. quintana, B. elizabethae, B. vinsonii subsp. berkhoffii, B. vinsonii subsp. arupensis, and B. bacilliformis were purchased from the American Type Culture Collection, Manassas, Va. Bartonella strains were cultivated on Trypticase soy agar supplemented with 5% sheep blood. Plates were incubated at 37°C in 5% CO2 until growth was visible, usually between 5 and 10 days. B. bacilliformis was incubated at 26°C under ambient atmosphere. Antiserum was generated from mice infected with Bartonella as follows. Female BALB/cJ mice (The Jackson Laboratory, Bar Harbor, Maine) were inoculated intraperitoneally with live Bartonella strain Pm15590co resuspended in 0.5 ml of phosphate-buffered saline (PBS) following harvesting from a plate culture. A booster injection with live Pm15590co was given 5 weeks after the primary injection, followed by a third booster 6 weeks later with strain Sb944nv. The mice were bled 3 weeks later and provided a source of anti-Bartonella polyclonal antiserum. Polyclonal antibodies against each of strains Sb944nv and Pm136co were generated in a similar fashion. The rodent-associated strains of Bartonella were used in this initial experiment because there was an interest in identifying specific antigens related to these types of Bartonella infections.
The progression from initially finding an immunoreactive clone from a rodent Bartonella strain genomic library to identifying the brp gene complex in B. vinsonii subsp. arupensis is summarized in Fig. 1. A positive expression clone, termed 6.1, was identified from the Bartonella strain na19203nm genomic library using the mouse antiserum described in the previous paragraph. Clone 6.1 contained a 1.5-kb insert and expressed a recombinant protein product of 58 kDa that reacted in an immunoblot with the same antibody used in the screening process (Fig. 1). DNA sequencing of the 6.1 insert revealed a truncated gene with no start or stop codon at the 5' or 3' end, respectively. To extend and obtain the DNA sequence of the 5' region of the 6.1 clone, probes specific to that area were used for hybridization to the genomic library, which yielded several new clones that were sequenced. It was found that these clones, although identical at the hybridizing ends, were entirely different from one another as the sequences extended. This was the first evidence that this gene was composed of repeat regions and/or was part of a paralogous gene family.
|
View larger version (14K): [in a new window] |
FIG. 1. Flow chart of steps describing the identification of the initial 6.1 clone through the DNA sequencing of the 26-kb Bva-Brp fosmid clone. (1) Bartonella strain na19103nm genomic lambda library screened with antiserum from a mouse infected with Bartonella strains isolated from Peromyscus and Spermophilus reveals several positive plaques. A positive clone, 6.1, was selected for further characterization. (2) Western blot of the 58-kDa recombinant protein expressed by clone 6.1 reacted with the antiserum used in the library screening in step 1. Lane A, Escherichia coli with expression vector with no insert; lane B, E. coli with the 6.1 cloned insert. (3) The 1.5-kb insert of 6.1 is sequenced. Repeat regions are recognized. (4) PCR primers are synthesized from the 6.1 termini and used for amplification from Bartonella species genomic DNA. A similar-sized amplicon is observed from B. vinsonii subsp. arupensis (lane d). Lanes: a, B. quintana; b, B. henselae; c, B. vinsonii subsp. berkhoffii; d, B. vinsonii subsp. arupensis; e, B. elizabethae; f, B. bacilliformis; g, strain na19103 nm; h, clone 6.1. (5) B. vinsonii subsp. arupensis fosmid library is hybridized with the clone 6.1 probe to identify large inserts harboring analogous region. Positive clones are selected for analysis. (6) A B. vinsonii subsp. arupensis fosmid clone containing a 26-kb insert, Bva-Brp, was sequenced and analyzed.
|
Genomic DNA was purified from Bartonella strains for Southern blot analysis by phenol-chloroform extraction according to standard methods. Two µg of genomic DNA was digested to completion with BamHI or XbaI restriction enzymes (New England Biolabs, Waverly, Mass.) and fractionated on a 0.8% Tris-acetate-EDTA agarose gel, whereby the gel was depurinated, denatured, and neutralized by respective soaking in 0.25 M HCl; 0.5 M NaOH, 1.5 M NaCl; and 0.5 M Tris, pH 7.5, 1.5 M NaCl. Hybridization was performed utilizing the AlkPhos direct labeling reagents (Amersham Biosciences, Little Chalfont, Buckinghamshire, England), and the hybridization probe was the 1.5-kb PCR amplicon amplified from B. vinsonii subsp. arupensis genomic DNA using primers designed from the 6.1 clone of strain na19103nm (Fig. 1). The forward primer, 6.1F, was 5' ATTGCCATTGGTAAGGAAACAGGC 3'. The reverse primer, 6.1R, was 5' TTGTAAGGGCACTGACCTCTTCTG 3'. Southern blot analysis revealed orthologs of this gene present in various Bartonella species known to infect humans (Fig. 2). Hybridizing bands were present in B. henselae, B. elizabethae, B. vinsonii subsp. arupensis, and B. vinsonii subsp. berkhoffii, with less prominent hybridizing bands in B. quintana and B. bacilliformis. A single distinct B. vinsonii subsp. arupensis BamHI fragment of approximately 25 kb was observed in the Southern blot; therefore, a fosmid library was constructed with that genomic DNA to ensure a cloned insert large enough to harbor the entire gene plus the flanking regions. A B. vinsonii subsp. arupensis fosmid library was created with the CopyControl fosmid library production kit (Epicentre, Madison, Wis.). Purified genomic DNA was randomly sheared, and fragments were size selected from 25 kb to 40 kb. The sheared DNA was blunt ended and ligated into the fosmid vector, pCC1FOS, and packaged using MaxPlax lambda packaging extracts (Epicentre). Packaged phage particles were allowed to infect E. coli EPI300-T1 (Epicentre), with subsequent colony plating on chloramphenicol selection plates. Three positive fosmid clones were selected by colony hybridization, analyzed by BamHI restriction digestion, and rehybridized with the 1.5-kb probe to ensure the presence of the band seen in the Southern blot. Once confirmed, one of the clones containing a 26-kb insert, Bva-Brp (Fig. 1), was sequenced in its entirety. Purified Bva-Brp fosmid DNA was sequenced by primer walking on both strands to an average fourfold redundancy. Electrophoresis was performed using ABI (PE Biosystems, Foster City, Calif.) Big-Dye 3.1 dye chemistry and an ABI 3730XL automated DNA sequencer. Chromatogram data were analyzed using the program Phred (University of Washington) for base calling and quality assessment. Consensus sequence assembly was performed using the program Phrap, with manual editing being performed using the program Consed. Additional sequencing analyses were performed using SeqManII and MegAlign (DNASTAR, Madison, Wis.) for assembly, DNA and protein alignments, and dot plot matrix analyses of repetitive regions. Protein and DNA GenBank database search and alignments were performed using BLAST (http://www.ncbi.nlm.nih.gov/BLAST/).
![]() View larger version (79K): [in a new window] |
FIG. 2. Southern blot of Bartonella spp. genomic DNA. Genomic DNA samples were digested with either BamHI or XbaI as indicated and hybridized to the clone Bva-6.1 probe (for the location within the gene, see Fig. 3A). Bh = B. henselae, Bq = B. quintana, Be = B. elizabethae, Bvb = B. vinsonii subsp. berkhoffii, Bva = B. vinsonii subsp. arupensis, Bb = B. bacilliformis. Molecular weight markers in kb are noted on the left. A 23- to 25-kb weakly hybridizing band for B. quintana is present.
|
![]() View larger version (25K): [in a new window] |
FIG. 3. Organization of genes on the 26-kb Bva-Brp clone of B. vinsonii subsp. arupensis. (A) brpA, brpB, and brpC coding sequences are represented by rectangular bars, with the arrows showing the direction of transcription. The hatched bars within the coding sequences denote portions of the genes composed of repetitive sequences. Genes downstream of brpA are abbreviated as HP (hypothetical protein), EPP (exopolyphosphatase), ftsJ (cell division protein), and GBP (GTP-binding protein). The 409-bp intervening sequence between the stop codon of brpC and the start codon of brpB is shown. The bidirectional arrows indicate inverted repeats characteristic of transcriptional termination sites for brpC. Two putative brpB promoter sites based on consensus sequences for the ribosome binding site (RBS) (underlined), the 10 site (underlined with a dot-and-dash line), and the 35 site (heavily underlined line) are denoted. The region immediately downstream of the stop codon for brpB is shown with inverted repeats illustrated with the bidirectional arrows. The immediate upstream noncoding region is shown for brpA with putative RBS, 10, and 35 sites. The location of the original truncated Bva-6.1 clone that generated the recombinant protein used in this study is shown above brpA. A scale in kbp is above the diagram. (B) Signal peptides deduced from the amino acid sequences for brpA and brpB. Charged residues are italicized, the hydrophobic central region is in bold, and putative SPase I cleavage sites are underlined. Cleavage sites predicted by the SignalP program follow the underlined VNA and VVA of BrpA and BrpB, respectively.
|
A search of the GenBank protein database with BrpA found the most significant match with portions of a B. henselae protein encoded by a gene termed badA1, submitted as part of the genome sequencing project recently reported by Alsmark et al. (2). The gene product, BadA, has recently been described as a 340-kDa surface protein adhesin of B. henselae that mediates proangiogenic responses in infected cells (40). The Clustal W alignment program of DNASTAR indicated there was a 35.9% amino acid identity between BadA and BrpA. Portions of BrpA also shared significant homology with a recently described B. quintana family of four variably expressed outer membrane proteins (Vomps), A to D (51). The B. quintana Vomps, like BadA, also function as adhesins and mediate autoaggregation. The Vomp gene family encodes proteins of approximately 100 kDa, smaller than the Brp family. Both the B. quintana Vomps and B. henselae BadA contain repeated regions and were shown to be immunogenic, consistent with properties observed with the B. vinsonii subsp. arupensis Brp proteins. The B. henselae badA1 and the B. quintana Vomp gene family members are likely to be orthologous companions of the B. vinsonii subsp. arupensis brp gene family, and it is not surprising that the comparably smaller sized Vomp proteins may be a result of the genome reduction of B. quintana as noted by the genome descriptions of Alsmark et al. (2).
There were four additional genes downstream of brpA, each showing similarity to database proteins by BLAST analysis. With the exception of the gene immediately downstream from brpA, the other genes were conserved and identified from analogous genes present in several other prokaryotic organisms. The gene downstream of brpA encoded a hypothetical protein of 22.3 kDa and has highly conserved orthologous partners in B. henselae and B. quintana, with amino acid identities of 82 and 81%, respectively. Next downstream was a gene encoding a 49-kDa protein described as an exopolyphosphatase, with 79% amino acid identity to the corresponding B. henselae gene. There was no BLAST match to this protein in B. quintana. Following the exopolyphosphatase gene were genes encoding the cell division protein ftsJ (also termed a 23S rRNA methyltransferase) and the GTP-binding protein. Both of these genes were present in B. henselae and B. quintana.
Repeat regions in brp genes. brpA, brpB, and brpC were each found to contain multiple extensive regions of internal repeats sometimes encoding long stretches of repetitive amino acid residues. These regions were not necessarily tandemly oriented but were separated by other repeats in no particular order or arrangement. Additionally, the repeats, although quite similar, were not identical. Figure 4 graphically illustrates the complexity of the repeat regions by dot plot matrix analysis of the brp nucleotide sequences. The dot plot of the entire 26.4-kb insert plotted against itself clearly shows overlapping homologous areas within each gene and also the homology among genes. Dot plot analysis of the amino acid sequence of BrpA reveals the extensive repetitive nature of various lengths of sequence throughout the deduced protein sequence (Fig. 5A). For example, the longest repeats were seen at amino acid positions 1279 to 1504, 1800 to 2022, 2504 to 2738, and 2825 to 3046 (Fig. 5A). These particular repeats range from 221 to 234 amino acids and are aligned in Fig. 5B. The percentages of identity between each repeat ranged from 62.8% (amino acids 2504 to 2738 and 2825 to 3046) to 75.8% (amino acids 1800 to 2022 and 2504 to 2738). Several more smaller repetitive regions are interspersed throughout the BrpA coding region, as illustrated in the dot plot matrix. The apparent randomness of the repetitive arrangements, with some longer than 200 amino acids, are a unique characteristic of these genes. The larger brpA repeats are not tandemly oriented but are interspaced at various lengths, as are the many shorter repetitive regions. The repetitive regions of BrpA are limited to the internal majority of the protein, with the amino and carboxy ends of the protein having few or no repeats. The recent publication describing the genomes of B. henselae and B. quintana reported that these organisms contain repeated or partially repeated gene families, with B. henselae containing the higher proportion, with 78 repeat families of 2 to 14 members (2). The precise number of brp genes in B. vinsonii subsp. arupensis remains to be counted, and comparison of their organization and arrangements to the badA1 and Vomp gene families within Bartonella species remains to be performed.
![]() View larger version (37K): [in a new window] |
FIG. 4. Dot plot matrix analysis at the nucleotide level of the 26-kb Bva-Brp insert plotted against itself. (A) Analysis performed at an 80% homology level. Each line represents an area with a corresponding homologous area with an identity level of at least 80%. The bold diagonal in the center represents the 100% match of the sequence with itself. Areas illustrating the locations of brpA, brpB, and brpC are denoted along the diagonal. The boxed regions of the matrix indicate homologous segments to brpA within brpA, brpB, and brpC. (B) Same as panel A, except the analysis was performed at the 70% identity level, which illustrates the areas of homology within and among the genes more clearly. Numbers on the x and y axes denote the nucleotide positions of the 26-kb Bva-Brp cloned insert.
|
![]() View larger version (35K): [in a new window] |
FIG. 5. (A) Dot plot analysis at the amino acid level of BrpA plotted against itself at the 70% identity level. The diagonal in the center indicates the 100% match of the sequence with itself. Circles surround the largest repeated regions corresponding to amino acid positions 1279 to 1504, 1800 to 2022, 2504 to 2738, and 2825 to 3046. Numbers on the x and y axes denote the amino acid positions for BrpA. (B) Amino acid alignment of repeats 1279 to 1504, 1800 to 2022, 2504 to 2738, and 2825 to 3046. Residues in bold italics represent divergence from the consensus.
|
![]() View larger version (50K): [in a new window] |
FIG. 6. (A) Immunoblot of whole-cell Bartonella lysates with anti-Bva 6.1. Anti-Bva-6.1 was used at a 1:1,000 dilution. Proteins were fractionated on a 7.5% polyacrylamide gel. Molecular size markers (in kilobases) are noted at the left. (B) Recombinant Bva-6.1 protein blotted against antisera from mice infected with rodent Bartonella strains. Sources of antibodies are indicated. Antibodies were used at a 1:500 dilutin. The lower reacting bands in the anti-Bva-6.1 blot are breakdown products of recombinant BrpA. (C) Immunoblot of recombinant Bva-6.1 protein blotted with polyclonal antisera generated against Bartonella sp. whole-cell lysates. Antibodies were used at a 1:500 dilution. Bh = B. henselae, Bq = B. quintana, Bb = B. bacilliformis, Be = B. elizabethae, Bvv = B. vinsonii subsp. vinsonii.
|
B. vinsonii subsp. arupensis and B. quintana cells were subjected to immunofluorescent staining using anti-Bva-6.1 protein (Fig. 7). B. vinsonii subsp. arupensis and B. quintana cells from a culture plate were separately mixed in 10 µl sterile PBS on a microscope slide and allowed to air dry. The mouse antibody generated against the recombinant Bva-6.1 was incubated with the cells at a 1:50 dilution at 37°C for 1 h. Following five washes in PBS, the cells were incubated with a 1:50 dilution of fluorescein isothiocyanate-conjugated goat anti-mouse immunoglobulin G at 37°C for 1 h. After another series of PBS washes, the slide was air dried and viewed by epifluorescence microscopy. The staining showed strong fluorescent reactivity with B. vinsonii subsp. arupensis cells, providing evidence for an outer surface localization for at least the part of the protein to which the antibody was directed. The antibody reactivity was concentrated in discrete areas of the cell, usually at the ends, indicating a polarized location for BrpA following secretion or export to the outer surface (Fig. 7A). B. quintana cells stained with the same antibody did not demonstrate fluorescent activity (Fig. 7B), reflecting the result seen in the immunoblot. Epitopes in the Bva-6.1 segment may also be present in other Brp proteins; therefore, these could be detected on the surface as well.
![]() View larger version (28K): [in a new window] |
FIG. 7. Immunofluorescent staining of Bartonella cells. (A) B. vinsonii subsp. arupensis stained with anti-Bva 6.1; (B) B. quintana stained with anti-Bva 6.1; (C) B. vinsonii subsp. arupensis stained with preimmunized mouse serum. Arrows in panel A point to the concentrated staining in the polar regions of the cells. Magnification, x1,000.
|
BrpA is encoded by one of the largest prokaryotic genes thus far identified, and to our knowledge, the brp family is novel in having genes encoding such large proteins with the extensive and complex mosaic of repeated domains. The definition of the numbers of brp gene family members in Bartonella spp., outer membrane protein locations, kinetics of expression, secretion mechanisms, the role of the repeat regions, and the biological function of the Brps are some of the many issues awaiting future investigations.
Nucleotide sequence accession number. The GenBank accession number for the Bva-Brp 26-kb clone is AY730759.
(Portions of this work were presented at the 18th meeting of the American Society for Rickettsiology, 28 September to 1 October 2003, Cumberland, Maryland, and at the 104th General Meeting of the American Society for Microbiology, 23 to 27 May 2004, New Orleans, Louisiana.)
Present address: Microbiology and Molecular Genetics Dept., University of Vermont, Burlington, VT 05405. ![]()
|
|
|---|
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»