Genome dynamics of Bartonella grahamii in micro-populations of woodland rodents
BMC Genomics volume 11, Article number: 152 (2010)
Rodents represent a high-risk reservoir for the emergence of new human pathogens. The recent completion of the 2.3 Mb genome of Bartonella grahamii, one of the most prevalent blood-borne bacteria in wild rodents, revealed a higher abundance of genes for host-cell interaction systems than in the genomes of closely related human pathogens. The sequence variability within the global B. grahamii population was recently investigated by multi locus sequence typing, but no study on the variability of putative host-cell interaction systems has been performed.
To study the population dynamics of B. grahamii, we analyzed the genomic diversity on a whole-genome scale of 27 B. grahamii strains isolated from four different species of wild rodents in three geographic locations separated by less than 30 km. Even using highly variable spacer regions, only 3 sequence types were identified. This low sequence diversity contrasted with a high variability in genome content. Microarray comparative genome hybridizations identified genes for outer surface proteins, including a repeated region containing the fha gene for filamentous hemaggluttinin and a plasmid that encodes a type IV secretion system, as the most variable. The estimated generation times in liquid culture medium for a subset of strains ranged from 5 to 22 hours, but did not correlate with sequence type or presence/absence patterns of the fha gene or the plasmid.
Our study has revealed a geographic microstructure of B. grahamii in wild rodents. Despite near-identity in nucleotide sequence, major differences were observed in gene presence/absence patterns that did not segregate with host species. This suggests that genetically similar strains can infect a range of different hosts.
Emerging infectious diseases have increased significantly during recent decades, with major implications for human health and global economy . A majority has been caused by bacteria, and many of these are zoonotic, i.e. they are accidentally transmitted to humans from other animal reservoirs in which they do not cause disease . Bacteria of the genus Bartonella represent an excellent model system for studies of host adaptation patterns in zoonotic agents. Bartonella infect the red blood cells of a wide variety of mammals, and are transmitted among hosts by blood-sucking arthropods. The genus contains two recognized human-specific pathogens, Bartonella bacilliformis and Bartonella quintana, the causative agents of Carrion's disease and trench fever, respectively. Bartonella henselae is naturally adapted to felines but can incidentally infect humans, manifested as cat-scratch disease. Recently, an increasing number of Bartonella species with natural animal reservoirs have been associated with disease manifestations in humans .
Rodent-associated Bartonella are particularly attractive for studies of emerging infectious diseases, given the abundance of rodents in nature, their importance as a reservoir for infectious diseases  and the observation that they carry several pathogenic Bartonella species [4–9]. Field studies in many different countries have shown that Bartonella is present in most rodent populations, with conflicting results regarding host specificity [10–18]. These studies have revealed that Bartonella grahamii is one of the most prevalent Bartonella species in wild rodents. B. grahamii infects many species of mice and voles, is transmitted by the rodent flea Ctenophthalmus nobili and has been associated with two reported cases of human disease [5, 6].
The genome of B. grahamii strain as4aup, isolated from a wood mouse in Sweden, was recently sequenced and found to consist of a 2.3 Mb circular chromosome and a 28 kb plasmid, pBGR3 . This genome is slightly smaller than the 2.6 Mb genome of the closely related rat-associated species Bartonella tribocorum, but larger than the genomes of the three major human pathogens  (TIGR, unpublished). Sixteen genomic islands (BgGI 1-16), containing many genes for type IV and V secretion systems and phage genes, were identified in B. grahamii by comparison to other published Bartonella genomes . Many of these genes are also present in genomic islands in B. henselae, albeit in lower copy numbers. Comparative microarray hybridizations of a global collection of B. henselae strains and one strain of its close relative Bartonella koehlerae, showed that genomic islands have been lost independently from different lineages [23, 24]. Pulsed-field gel electrophoresis of the B. henselae strains revealed numerous rearrangements across the terminus of replication with breakpoints in regions encompassing the genomic islands . In addition to these genome-scale analyses, the intra-species diversity has been studied by multi locus sequence typing (MLST) or sequencing of variable spacer regions (MST) in B. henselae, B. quintana and B. grahamii[23, 25–29].
MLST schemes have been developed for all major human pathogens, and it has been the standard technique for classification of bacterial isolates [30, 31]. However, while being useful for inferring the population structure, MLST only rarely correlates with virulence properties . In contrast, microarray-based studies offer the opportunity to analyze the complete genome of closely related strains, and have been applied to study the diversity and the difference in virulence properties of many bacterial pathogens . For both types of analysis, the standard has been to include a diverse set of strains, often collected at different time points in different parts of the world. While this strategy does well in grasping as much variability as possible, it gives a very limited understanding of the mechanisms and rates whereby hyper-variable sequences evolve, and of the environmental factors influencing these processes.
In this study, we have used a collection of B. grahamii strains isolated from wild rodents in three adjacent but non-contiguous geographic areas to analyze the genome content and the geographic microstructure of natural bacterial populations. The results revealed a geographic pattern of these strains, with few SNPs contrasting with dramatic variability in genome structure and repertoire of secretion systems.
We have analyzed 27 B. grahamii isolates collected from three geographic sites, separated by circa 30 km, in central Sweden (Table 1) . Fourteen isolates were sampled in Håtunaholm, nine in Kumla and four in Ålbo. Sixteen were obtained from yellow-necked mice (Apodemus flavicollis), five from each of wood mice (Apodemus sylvaticus) and bank voles (Myodes glareolus), and one from a house mouse (Mus musculus). The prevalence of B. grahamii in the 236 captured rodents ranged between 8-16% among sites, and 6-20% among host species .
Few sequence polymorphisms in geographically adjacent B. grahamii populations
Partial sequencing of four housekeeping genes (gltA, ftsZ, batR and cycK) revealed only a single polymorphism, located in ftsZ (Table 1, Table 2). To further differentiate the strains, we sequenced 13 loci expected to be more variable (Table 2), most of which correspond to intergenic regions that were successfully used for typing of B. henselae strains [28, 29]. Polymorphisms were observed in only three of these regions (Table 1). In total, there were 5 SNPs in the 12 kb alignment, with three distinct sequence types (STs). We observed no sign of host specificity, for example the ST1-strains were sampled from all four host species (Table 1). In contrast, the results indicate some degree of geographic pattern, with all strains from Håtunaholm belonging to ST1, together with strains af115up and cg147up from Kumla. ST2 and ST3 were present in both Kumla and Ålbo.
Differences in gene repertoires for type IV and type V secretion systems
We screened for gene content variations using comparative genome hybridizations (CGH) to a microarray covering 96% of the genes in B. grahamii as4aup . The results suggested that all ST1-strains have a virtually identical gene content (Figure 1), with the limitation that possible extra genes cannot be detected with CGH. All ST2 and ST3-strains displayed decreased microarray signal in comparison to the reference strain in the so-called fha- repeat, which is present in 5 genomic islands in the sequenced genome . The fha/hec operon codes for a two-partner type V secretion system (T5SS), where the hec gene product probably assists in the secretion of filamentous hemagglutinin (FHA). In addition, the fha- repeat contains a number of putative virulence genes and phage genes .
The lower microarray signal in the fha-repeat can be due to a complete loss of these genes, a lower copy number or a high sequence divergence compared to the reference strain. Given the few identified SNPs and the observation that several other rapidly evolving genes, e.g. badA, did not show any indications of low signal, it is unlikely that the weak signal in the fha-repeat is due to sequence divergence. Attempts to amplify different fha genes by PCR using several sets of primers, all of which were tested successfully in the sequenced strain, resulted in one single PCR product in strain af140up and no products in strains af165up or as211up (see Additional file 1). Thus, all ST2 and ST3-strains are likely to have few or no copies of these genes.
In contrast, several strains displayed a relative increase in hybridization signal in the region containing the multicopy badA T5SS genes (Figure 1). BadA is a trimeric autotransporter that is essential for adhesion to host cells and bloodstream infection in Bartonella species [21, 34–37]. Since the badA region is already much longer in B. grahamii than in any other sequenced Bartonella genome , the finding that some strains seem to have additional copies of this gene was unexpected. Variability was also observed in one of the type IV secretion systems (T4SS), namely the plasmid-encoded vbh T4SS (Figure 2). The entire plasmid appeared to be absent from the strains af206up (ST2) and as211up (ST3) from Ålbo (Figure 2).
The ST2 and ST3-strains also displayed reduced hybridization signals in genes for two families of hypothetical proteins that are repeated in the genomes of B. grahamii and B. tribocorum, but have not been identified in any other sequenced Bartonella genome . The first family (e.g. Bgr_15790) is present in a total of 17 copies in the sequenced strain, located adjacent to some of the fha genes. Homologs have mainly been identified in Gamma-proteobacteria, e.g. in a genomic island in Photorhabdus luminescens, in Pseudomonas syringae and Listeria monocytogenes, in one single copy in each genome. The second family (e.g. Bgr_18020) is present in 6 copies in the sequenced strain, located between vbl 5 and vbl 6 in the chromosomal vbh gene cluster in BgGI 16 and on pBGR3, as well as shortly upstream of the chromosomal virB system. No homologs to this gene have been found in other genera. The remaining variability is mostly due to phages, and includes loss/divergence of single phage genes, replication of complete prophages in BgGI 2 and BgGI 4 and putatively phage-derived bi-directional run-off replication around BgGI 15, in line with previous observations  (Figure 1).
Genome structure diversity
Pulsed-field gel electrophoresis (PFGE) with the NotI enzyme revealed four restriction patterns in 20 representative strains and a relatively homogenous genome size of the isolates (Figure 3). One pattern was shared by all ST1-strains, two were unique to three strains from Kumla and Ålbo, respectively, and the last was shared between the remaining strains from Kumla and Ålbo. This result suggests some degree of correlation between genome structure and geographic origin. However, one of the restriction sites for the largest NotI band is located very close to the terminus of replication (Figure 4), and rearrangements around the terminus may thus have escaped detection. Therefore, we also employed the SgfI enzyme, which has more restriction sites than NotI in the sequenced genome of B. grahamii strain as4aup. With this enzyme, we observed a much greater diversity with 15 distinct patterns that showed no obvious correlation with geographic origin (Figure 3). This suggests that genome rearrangements occur frequently, in line with previous observations in B. henselae and B. quintana[23, 27].
The 490 and 610 kb fragments obtained after cleavage of the sequenced strain with the SgfI enzyme were found to be the least conserved in the other strains (Figure 3). Both these bands contain one copy of the main prophage (Figure 4). The 490 kb band also contains several copies of the fha-repeat in BgGI 4 and some additional phage genes. Rearrangements within and across duplicated genes in these islands as well as across the terminus of replication could explain the variability observed after digestion with the SgfI enzyme. The other copies of the fha-repeat are located in BgGI 7, BgGI 10, BgGI 11 and BgGI 14. The three latter are encompassed in the 760 kb fragment of the SgfI digest and in the 120 kb and 200 kb bands of the NotI digest of strain as4aup. It is interesting to note that the 120 kb and 200 kb bands are identical in size in all ST1 strains that contain the fha-repeat but differ in the ST2 and ST3-strains, which have no or fewer copies of these genes.
Genomotypes do not have a random geographical distribution
We defined the genomotype based on the ST and the two major gene content variations: presence/absence of the fha-repeat and presence/absence of the plasmid, where absence also includes cases of lower copy number. Figure 5 summarizes the distribution of the genomotypes across geographical sites and the rodent hosts associated with each type. The most striking result is that only one genomotype was recovered from Håtunaholm (type I, which is of ST1 and harbors both the fha-repeat and the plasmid), although these isolates were obtained from all four different rodent host species. Genomotype I was also present in two isolates from Kumla, while absent in Ålbo. A more diverse set of isolates was identified in Kumla and Ålbo with three and four genomotypes, respectively. The fha-repeat displayed low signal in all ST2 and ST3-strains, suggesting that it was lost or gained only once. In contrast, the plasmid was variably present in both ST2 and ST3 and is likely to have experienced a minimum of two gain/loss events, unless one of the variants is the result of a recombination event.
B. grahamii strains exhibit large differences in growth rate in liquid culture
In order to compare genomic differences to phenotypic traits, we selected five representative strains and studied their growth in liquid culture for ten days (Figure 6). We observed large differences in growth rate, with the only strain isolated from Mus musculus (mm3up) growing significantly slower than the others and reaching a lower maximal optical density (Table 3). The maximal generation times, ranging from approximately 5 (as206up) to 22 (mm3up) hours, did not correlate with ST, presence/absence of the fha-repeat or presence/absence of the plasmid. We also performed aggregation assays for three strains (as4aup, af165up and af206up), but observed no significant differences between them (data not shown).
In this study, we have analyzed the genomic diversity of B. grahamii from wild rodents in three adjacent geographic areas in central Sweden. To our knowledge, this is the first genome-scale analysis of bacterial populations from wild animals in a geographically restricted area. Since these rodents typically have a home area of less than 1 hectare, and seldom migrate distances longer than 30 km [38, 39], genetic exchange between the three sites is expected to be limited. Indeed, we observed a geographic pattern in the distribution of genomotypes, showing that microstructures of bacterial populations inhabiting hosts with limited migration are visible also in the comparison of closely located sites. One of the sites, Håtunaholm, displayed a lower degree of genomic diversity than the others, with all isolates being of the same genomotype. In contrast, the presence of strains of ST2 and ST3 in both Kumla and Ålbo suggests that exchange has occurred at least twice between these regions. A possible reason for the clonal population in Håtunaholm is that this site is surrounded by water and only connected to Kumla and Ålbo by a narrow piece of land, where the small town Bålsta is situated (Figure 5). This environmental barrier is likely to restrict the genetic flow by isolation of hosts and vectors. The presence of two ST1-strains in Kumla suggests that the Håtunaholm population is a recent expansion of an immigrant from the Kumla region, or that the clonal variants in Håtunaholm occasionally escape the inland seawater barriers.
Despite sequencing regions expected to be among the most variable in the genome, we observed only five SNPs in 27 strains, corresponding to only three STs. As expected, this is much less diversity than what was observed in a previous MLST analysis of a global collection of 31 B. grahamii strains, where 10-16 variants were identified for each gene . However, previous MLST analysis of 7 house-keeping genes in the likewise intracellular Borrelia lucitaiae collected from ticks on lizards in two nearby geographic areas identified 13 types in 16 isolates , suggesting that the low sequence diversity observed in our population is not necessarily a general feature of vector-borne bacterial populations in restricted geographic areas. Future studies of similar bacterial populations will reveal the degree of diversity that they typically display, and help elucidate the factors that govern these dynamics.
In light of the limited sequence diversity, the gene content of the studied populations was very dynamic, with the major differences being associated with the plasmid pBGR3 and the fha-repeat. Thus, gene gain and loss occur so frequently in nature that even strains sampled from the same geographic site with essentially no nucleotide sequence differences in a subset of their housekeeping genes can have major differences in gene content. This result shows that traditional MLST may not always be appropriate for inferring phenotypic or virulence properties, and illustrates the importance of a whole-genome approach when analyzing bacterial diversity. Furthermore, our results show that comparisons of genomic diversity in different geographic regions should be based on multiple strains from each location, since selecting only one or a few strains from each region for sequencing might easily underestimate the true genomic diversity. The real diversity could also be underestimated if there is considerable within-host heterogeneity and each host is represented by only one or a few strains. However, for B. grahamii we consider such within-host genomic diversity unlikely since a random sampling of a heterogeneous micro-population would recover more than a single genomotype from Håtunaholm.
The function of FHA is not known in Bartonella, but by homology inference to Bordetella, we assume that it acts as an adhesin and/or immunomodulator. In pathogenic Bordetella species, the fhaB/fhaC genes encode a two-partner type V secretion system, where FhaB is translocated to the cell surface with the aid of FhaC, which is homologous to Hec in Bartonella[42, 43]. The mature FHA protein is required for adhesion to epithelial cells and has a role in apoptosis and immunosuppression [44–48]. The fhaB/fhaC genes are located within a well conserved single copy gene cluster in all pathogenic Bordetella species, however, several paralogs to fhaB have been identified . Comparative genome hybridizations have shown no losses of fhaB/fhaC, but occasional losses of the fhaB paralogs [50–52]. In Bordetella bronchiseptica, it was suggested that one of the FhaB paralogs (FhaS) might be involved in host-specificity . Although FhaS is not required for infection, a ∆fhaS strain was out-competed by the wild type during coinfection .
One could think of a similar situation in Bartonella, namely that some of the various fha gene copies confer host specificity and/or are advantageous only under certain conditions. Such a scenario would be consistent with their location within genomic islands that have been integrated/excised multiple times independently in B. grahamii, B. henselae and B. tribocorum[20, 22], and with the many cases of pseudogenization of individual fha genes. Some degree of variability was observed among strains in B. henselae although none had completely lost the fha genes . In contrast, the fha genes have been lost from all strains of B. quintana investigated (data not shown) and are also absent from B. bacilliformis[20, 22]. The absence of the fha gene in several strains and species suggests that the gene may be lost from individual Bartonella populations, both temporarily and permanently, and is therefore not essential for growth and infection per se. Indeed, we observed no significant differences in either growth rate or level of autoaggregation between strains with different copy numbers of the fha gene.
Since rodents have an adaptive immune system, shedding of surface-exposed proteins could potentially enable super-infection, i.e. a second infection of a previously infected animal. Immunity against super-infection of the same Bartonella species has indeed been demonstrated in rodent hosts . Under this scenario, loss of the fha gene might be the only way to escape acquired immunity in the rodent host population. If selection for functionality is balanced by counter-selection from the host immune system, this could trigger a periodic cycle in which strains with and without fha repeatedly out-compete each other. The possibility for super-infection is likely to depend on the transmission rates, and thus the abundance of hosts and vectors. A high density of hosts and vectors will increase transmission rates, and might favor loss of fha genes from many individuals in a population. Vice versa, if hosts and vectors are low in abundance, super-infections will be rare and host-acquired immunity less of a problem, which might set the stage for the return of strains with fha.
The only putative host-interaction genes located on the plasmid are those encoding a vbh T4SS that is closely related to the chromosomal vbh system, and more distantly related to the virB system, which has been shown to be essential for infection in B. tribocorum and secrets Bartonella effector proteins into the host cell [55–57]. The plasmid also contains a set of tra genes, including the DNA binding traA, which are usually involved in conjugation. It is not known whether the vbh gene cluster encodes proteins involved in host interactions or if it is only part of a conjugation system. The presence of a putative addiction module, consisting of a killer and an antidote protein, on pBGR3 (Bgr_p00250-Bgr_p00260) may ensure the survival of the plasmid in the population even in the absence of selection.
The observation that the Håtunaholm strains were sampled from four different rodent species shows that B. grahamii isolates with similar gene contents are capable of infecting diverse hosts. Such a lack of host-specificity contrasts with experimental infections of the cotton rat (Sigmodon hispidus) and the white-footed mouse (Peromyscus leucopus), which indicates that Bartonella strains are only able to infect hosts that are closely related to the host from which it was originally isolated . Although S. hispidus and P. leucopus belong to the same family of rodents (Cricetidae), cross-infections were not observed . The rodents sampled in Håtunaholm belong to two different families, Cricetidae and Muridae, which are genetically more diverged , yet harbor similar B. grahamii strains. This suggests a lack of host-specificity, although minor sequence differences among strains may have gone undetected in our survey. Since the only parameter showing some correlation with the generation times was the host of isolation, it is possible that the growth rate of each strain has been optimized for its particular host. It remains to be shown, however, whether a specific isolate from Håtunaholm is actually capable of infecting both mice and voles, or if there are nucleotide sequence or expression differences in outer surface proteins that account for host specificity but were not detected in this analysis.
In the future, results on the natural genomic diversity of zoonotic bacterial populations should be analyzed in an ecological context, including data about e.g. host and vector abundance and climate, to allow correlations between genomic diversity and bacterial prevalence or seasonal variation. Our results also urge for functional characterization of the role of FHA in Bartonella, and offer an opportunity to study the infectious process of strains with observed differences in genes for surface-exposed proteins in rodent hosts.
Based on studies of SNPs and gene contents by comparative microarray genome hybridization, we present an illustration of how different genomotypes of B. grahamii are distributed across nearby geographic locations. Despite a very limited number of SNPs, we observed dramatic changes in terms of genome structure and content, with the major gains/losses being a plasmid encoding a vbh T4SS and a repeated region containing the fha/hec genes for a T5SS. Whereas a mixture of different genomotypes was recovered from rodents in Ålbo and Kumla, all strains isolated from rodents in Håtunaholm were identical in both ST and gene content. This suggests that environmental barriers in the form of water can lead to isolation and loss of genomic variability in host-associated bacterial populations. Neither the STs nor the gene content variability correlated with host of isolation, indicating a lack of host-specificity in these populations. However, the generation time was significantly longer for a strain isolated from M. musculus (more than 20 hours), than for strains isolated from Apodemus mice, only 5 to 12 hours. This might indicate growth rate optimization in response to the host cell environment. Given the variability in gene content and generation times for strains isolated from within such a small geographic region, we emphasize the importance of selecting more than a single strain for genomic and functional genomics analyses.
Bacterial strains and DNA isolation
All B. grahamii strains used in this study were isolated from wild rodents in September 1999 . For production of bacterial stock suspensions, bacteria were harvested from hematin agar plates after four days and stored at -80°C in a buffer containing 150 mM NaCl, 50 mM Tris (pH 7) and 11.2% glycerol. The number of viable bacteria in the frozen stocks was determined by plating serial dilutions of the suspension and calculating the number of colony forming units (CFU). DNA was extracted from bacteria grown for 5-10 days on hematin agar plates, as described previously .
Sequence data was collected for seventeen genomic regions, using primers listed in Table 2. Primers used for PCR of fha are listed in Additional file 1. For gltA and ftsZ, PCR amplifications were performed with a Perkin-Elmer GeneAmp 9600 thermocycler using a PCR Master kit from Boehringer Mannheim Scandinavia AB (Bromma, Sweden). The 50 μl reaction mixture consisted of the template, forward and reverse primers (10 pmol/primer), 25 μl PCR master mix, and distilled water of PCR-grade up to 50 μl. The reaction conditions were 5 min at 95°C, 10 cycles at 94°C for 20 s, 50°C for 1 min, 72°C for 90 s, annealing temperature lowered 1°C per cycle until reaching 40°C. This was followed by 40 cycles at 94°C for 20 s, 40°C for 1 min, 72°C for 90 s and a 7 min extension period at 72°C. PCR products were purified with the QIAquick purification kit from QIAGEN, Inc., Chatsworth, Calif., U.S., according to the manufacturer's instructions. PCR products were sequenced with an ABI 310 Genetic Analyzer (Perkin-Elmer Corp., Norwalk, Conn., US). The gltA product was sequenced with primers CS443f  and BhCS.781p  in addition to the primers used for PCR. Sequencing reactions were performed using a DNA Sequencing Kit with AmpliTaq® DNA Polymerase, FS, for the BigDye™ Terminator Cycle Sequencing Ready Reaction protocol (Perkin-Elmer, Applied Biosystems, Warrington, GB). Remaining PCR reactions were performed as described previously .
Analysis of sequence data
Sequences for gltA and ftsZ were analyzed with ABI Prism™ DNA Sequencing Analysis Software Version 3.0 (PE Applied Biosystems, Foster City, CA, US) and Sequencher™ (Gene Codes Corporation, Ann Arbor, MI, US). Remaining sequences were assembled and edited with Phred, Phrap and Consed [62–64]. The allele of the sequenced strain was named allele 1 and the combination of alleles was used to define the sequence types.
Microarray comparative genome hybridizations
Microarrays were designed and manufactured as reported previously  and cross-linked at 250 mJ/cm2 (slides printed at the Royal Institute of Technology) or 800 mJ (slides printed at Uppsala University). Genomic DNA from the sequenced strain was used as reference in all hybridizations and two or three hybridizations were performed for each strain. Prehybridization, DNA labeling, hybridization, scanning and image analysis were performed as described previously . The channel used for the reference strain is referred to as Ch1 and the test strain is referred to as Ch2, and spots meeting the following criteria were removed from further analysis: spots flagged as bad or not found during quantification, spots with more than 10% saturated pixels in either channel, Ch1 spot median below 2 times the Ch1 background median, less than 90% of pixels having a Ch1 intensity higher than background intensity plus one standard deviation, less than 75% of pixels having a Ch1 intensity higher than background intensity plus two standard deviations or spots having less than 70 pixels. M-values were computed as log2(Ch2/Ch1) and normalization was performed as described previously . Median M-values of all replicate spots and arrays were computed for each strain. A gene was defined as absent if the median M-value of all probes in that gene was ≤ -2.
Pulsed-field gel electrophoresis
Bacteria grown for 10 days were digested with 10 U of NotI (New England Biolabs) or SgfI (Promega) restriction endonucleases. Experimental procedures for digestion were as reported previously . The DNA fragments were separated in 0.9% PFGE-grade agarose (SeaKem® Gold; Cambrex Bio Science) in 0.5× TBE buffer in GenNavigator® System apparatus (Amersham Biosciences) at 14°C and 5.6 V/cm, for a total of 65 hours. The total run was separated in six phases; switch times ramped from 5 to 150 s. The sizes of the fragments were estimated using PFGE λ-ladder and Yeast Chromosome PFG marker (New England Biolabs).
Growth curves in liquid culture
For growth in liquid culture, frozen bacterial stocks were thawed on ice and inoculated at a concentration of approximately 106 CFU per ml into Schneider's insect medium (Sigma) supplemented with 10% FBS (heat inactivated at 56°C for 30 min) and 5% sucrose, as described by Riess et al. . In addition, the medium was supplemented with 0.1 M HEPES, to keep the pH near 7. Bacterial growth was determined by measuring the optical densities at 600 nm (OD600) in triplicates at 24-h intervals. For strain as4aup, the number of CFU was determined at each time point by plating 10-fold serial dilutions. We used the formula described by Reiss et al.  to estimate the maximum growth rate, expressed as generation time, for each strain, assuming a similar amount of bacteria at a certain optical density as for as4aup.
Bacteria grown in liquid cultures as described above were sampled in the mid-exponential phase. The bacterial cells were centrifuged, and pellets were suspended in phosphate buffered saline (PBS) to an OD600 of 0.6. After gentle vortexing, standing tubes with 5 ml of bacterial suspension were incubated at 35°C, and the OD600 of the upper part of the suspensions was measured throughout time (0 min, 5 min, 15 min, 1 h, 4 h and 22 h).
Nucleotide sequence and microarray accession numbers
The novel sequences of ftsZ and three spacers have been submitted to GenBank under the accession numbers FN597604. FN597607. The microarray data have been deposited in the ArrayExpress database of the European Bioinformatics Institute under the accession number the accession numbers A-MEXP-1741 and A-MEXP-1576 for the array designs and E-TABM-850 for the experimental data.
multi locus sequence typing
multi spacer typing
comparative genome hybridization
pulsed-field gel electrophoresis
single nucleotide polymorphism
Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, Daszak P: Global trends in emerging infectious diseases. Nature. 2008, 451: 990-993. 10.1038/nature06536.
Boulouis HJ, Chang CC, Henn JB, Kasten RW, Chomel BB: Factors associated with the rapid emergence of zoonotic Bartonella infections. Vet Res. 2005, 36: 383-410. 10.1051/vetres:2005009.
Woolhouse ME, Gowtage-Sequeria S: Host range and emerging and reemerging pathogens. Emerg Infect Dis. 2005, 11: 1842-1847.
Daly JS, Worthington MG, Brenner DJ, Moss CW, Hollis DG, Weyant RS, Steigerwalt AG, Weaver RE, Daneshvar MI, O'Connor SP: Rochalimaea elizabethae sp. nov. isolated from a patient with endocarditis. J Clin Microbiol. 1993, 31: 872-881.
Kerkhoff FT, Bergmans AM, Zee van Der A, Rothova A: Demonstration of Bartonella grahamii DNA in ocular fluids of a patient with neuroretinitis. J Clin Microbiol. 1999, 37: 4034-4038.
Serratrice J, Rolain JM, Granel B, Ene N, Conrath J, Avierinos JF, Disdier P, Raoult D, Weiller PJ: [Bilateral retinal artery branch occlusions revealing Bartonella grahamii infection]. Rev Med Interne. 2003, 24: 629-630. 10.1016/S0248-8663(03)00224-8.
Welch DF, Carroll KC, Hofmeister EK, Persing DH, Robison DA, Steigerwalt AG, Brenner DJ: Isolation of a new subspecies, Bartonella vinsonii subsp. arupensis, from a cattle rancher: identity with isolates found in conjunction with Borrelia burgdorferi and Babesia microti among naturally infected mice. J Clin Microbiol. 1999, 37: 2598-2601.
Fenollar F, Sire S, Raoult D: Bartonella vinsonii subsp. arupensis as an agent of blood culture-negative endocarditis in a human. J Clin Microbiol. 2005, 43: 945-947. 10.1128/JCM.43.2.945-947.2005.
Kosoy M, Murray M, Gilmore RD, Bai Y, Gage KL: Bartonella strains from ground squirrels are identical to Bartonella washoensis isolated from a human patient. J Clin Microbiol. 2003, 41: 645-650. 10.1128/JCM.41.2.645-650.2003.
Ellis BA, Regnery RL, Beati L, Bacellar F, Rood M, Glass GG, Marston E, Ksiazek TG, Jones D, Childs JE: Rats of the genus Rattus are reservoir hosts for pathogenic Bartonella species: an Old World origin for a New World disease?. J Infect Dis. 1999, 180: 220-224. 10.1086/314824.
Ying B, Kosoy MY, Maupin GO, Tsuchiya KR, Gage KL: Genetic and ecologic characteristics of Bartonella communities in rodents in southern China. Am J Trop Med Hyg. 2002, 66: 622-627.
Castle KT, Kosoy M, Lerdthusnee K, Phelan L, Bai Y, Gage KL, Leepitakrat W, Monkanna T, Khlaimanee N, Chandranoi K: Prevalence and diversity of Bartonella in rodents of northern Thailand: a comparison with Bartonella in rodents from southern China. Am J Trop Med Hyg. 2004, 70: 429-433.
Pretorius AM, Beati L, Birtles RJ: Diversity of bartonellae associated with small mammals inhabiting Free State province, South Africa. Int J Syst Evol Microbiol. 2004, 54: 1959-1967. 10.1099/ijs.0.03033-0.
Jardine C, Appleyard G, Kosoy MY, McColl D, Chirino-Trejo M, Wobeser G, Leighton FA: Rodent-associated Bartonella in Saskatchewan, Canada. Vector Borne Zoonotic Dis. 2005, 5: 402-409. 10.1089/vbz.2005.5.402.
Inoue K, Maruyama S, Kabeya H, Yamada N, Ohashi N, Sato Y, Yukawa M, Masuzawa T, Kawamori F, Kadosaka T: Prevalence and genetic diversity of Bartonella species isolated from wild rodents in Japan. Appl Environ Microbiol. 2008, 74: 5086-5092. 10.1128/AEM.00071-08.
Birtles RJ, Harrison TG, Molyneux DH: Grahamella in small woodland mammals in the U.K.: isolation, prevalence and host specificity. Ann Trop Med Parasitol. 1994, 88: 317-327.
Holmberg M, Mills JN, McGill S, Benjamin G, Ellis BA: Bartonella infection in sylvatic small mammals of central Sweden. Epidemiol Infect. 2003, 130: 149-157. 10.1017/S0950268802008075.
Engbaek K, Lawson PA: Identification of Bartonella species in rodents, shrews and cats in Denmark: detection of two B. henselae variants, one in cats and the other in the long-tailed field mouse. Apmis. 2004, 112: 336-341. 10.1111/j.1600-0463.2004.apm1120603.x.
Bown KJ, Bennet M, Begon M: Flea-borne Bartonella grahamii and Bartonella taylorii in bank voles. Emerg Infect Dis. 2004, 10: 684-687.
Berglund EC, Frank AC, Calteau A, Vinnere Pettersson O, Granberg F, Eriksson AS, Naslund K, Holmberg M, Lindroos H, Andersson SG: Run-off replication of host-adaptability genes is associated with gene transfer agents in the genome of mouse-infecting Bartonella grahamii. PLoS Genet. 2009, 5: e1000546-10.1371/journal.pgen.1000546.
Saenz HL, Engel P, Stoeckli MC, Lanz C, Raddatz G, Vayssier-Taussat M, Birtles R, Schuster SC, Dehio C: Genomic analysis of Bartonella identifies type IV secretion systems as host adaptability factors. Nat Genet. 2007, 39: 1469-1476. 10.1038/ng.2007.38.
Alsmark CM, Frank AC, Karlberg EO, Legault BA, Ardell DH, Canback B, Eriksson AS, Naslund AK, Handley SA, Huvet M: The louse-borne human pathogen Bartonella quintana is a genomic derivative of the zoonotic agent Bartonella henselae. Proc Natl Acad Sci USA. 2004, 101: 9716-9721. 10.1073/pnas.0305659101.
Lindroos H, Vinnere O, Mira A, Repsilber D, Naslund K, Andersson SG: Genome rearrangements, deletions, and amplifications in the natural population of Bartonella henselae. J Bacteriol. 2006, 188: 7426-7439. 10.1128/JB.00472-06.
Lindroos HL, Mira A, Repsilber D, Vinnere O, Naslund K, Dehio M, Dehio C, Andersson SG: Characterization of the genome composition of Bartonella koehlerae by microarray comparative genomic hybridization profiling. J Bacteriol. 2005, 187: 6155-6165. 10.1128/JB.187.17.6155-6165.2005.
Iredell J, Blanckenberg D, Arvand M, Grauling S, Feil EJ, Birtles RJ: Characterization of the natural population of Bartonella henselae by multilocus sequence typing. J Clin Microbiol. 2003, 41: 5071-5079. 10.1128/JCM.41.11.5071-5079.2003.
Arvand M, Feil EJ, Giladi M, Boulouis HJ, Viezens J: Multi-locus sequence typing of Bartonella henselae isolates from three continents reveals hypervirulent and feline-associated clones. PLoS ONE. 2007, 2: e1346-10.1371/journal.pone.0001346.
Foucault C, La Scola B, Lindroos H, Andersson SG, Raoult D: Multispacer typing technique for sequence-based typing of Bartonella quintana. J Clin Microbiol. 2005, 43: 41-48. 10.1128/JCM.43.1.41-48.2005.
Li W, Chomel BB, Maruyama S, Guptil L, Sander A, Raoult D, Fournier PE: Multispacer typing to study the genotypic distribution of Bartonella henselae populations. J Clin Microbiol. 2006, 44: 2499-2506. 10.1128/JCM.00498-06.
Li W, Raoult D, Fournier PE: Genetic diversity of Bartonella henselae in human infection detected with multispacer typing. Emerg Infect Dis. 2007, 13: 1178-1183.
Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA: Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998, 95: 3140-3145. 10.1073/pnas.95.6.3140.
Maiden MC: Multilocus sequence typing of bacteria. Annu Rev Microbiol. 2006, 60: 561-588. 10.1146/annurev.micro.59.030804.121325.
Turner KM, Feil EJ: The secret life of the multilocus sequence type. Int J Antimicrob Agents. 2007, 29: 129-135. 10.1016/j.ijantimicag.2006.11.002.
Dorrell N, Hinchliffe SJ, Wren BW: Comparative phylogenomics of pathogenic bacteria by microarray analysis. Curr Opin Microbiol. 2005, 8: 620-626. 10.1016/j.mib.2005.08.012.
Riess T, Andersson SG, Lupas A, Schaller M, Schafer A, Kyme P, Martin J, Walzlein JH, Ehehalt U, Lindroos H: Bartonella adhesin a mediates a proangiogenic host cell response. J Exp Med. 2004, 200: 1267-1278. 10.1084/jem.20040500.
Riess T, Raddatz G, Linke D, Schafer A, Kempf VA: Analysis of Bartonella adhesin A expression reveals differences between various B. henselae strains. Infect Immun. 2007, 75: 35-43. 10.1128/IAI.00963-06.
Kaiser PO, Riess T, Wagner CL, Linke D, Lupas AN, Schwarz H, Raddatz G, Schafer A, Kempf VA: The head of Bartonella adhesin A is crucial for host cell interaction of Bartonella henselae. Cell Microbiol. 2008, 10: 2223-2234. 10.1111/j.1462-5822.2008.01201.x.
Zhang P, Chomel BB, Schau MK, Goo JS, Droz S, Kelminson KL, George SS, Lerche NW, Koehler JE: A family of variably expressed outer-membrane proteins (Vomp) mediates adhesion and autoaggregation in Bartonella quintana. Proc Natl Acad Sci USA. 2004, 101: 13630-13635. 10.1073/pnas.0405284101.
Szacki J, Liro A: Movements of small mammals in the heterogenous landscape. Landscape Ecology. 1991, 5: 219-224. 10.1007/BF00141436.
Corbet GB, Harris S: The Handbook of British Mammals. 1991, Blackwell Scientific Publications, Oxford, Third
Inoue K, Kabeya H, Kosoy MY, Bai Y, Smirnov G, McColl D, Artsob H, Maruyama S: Evolutional and geographical relationships of Bartonella grahamii isolates from wild rodents by multi-locus sequencing analysis. Microb Ecol. 2009, 57: 534-541. 10.1007/s00248-009-9488-x.
Vitorino LR, Margos G, Feil EJ, Collares-Pereira M, Ze-Ze L, Kurtenbach K: Fine-scale phylogeographic structure of Borrelia lusitaniae revealed by multilocus sequence typing. PLoS One. 2008, 3: e4002-10.1371/journal.pone.0004002.
Jacob-Dubuisson F, El-Hamel C, Saint N, Guedin S, Willery E, Molle G, Locht C: Channel formation by FhaC, the outer membrane protein involved in the secretion of the Bordetella pertussis filamentous hemagglutinin. J Biol Chem. 1999, 274: 37731-37735. 10.1074/jbc.274.53.37731.
Jacob-Dubuisson F, Kehoe B, Willery E, Reveneau N, Locht C, Relman DA: Molecular characterization of Bordetella bronchiseptica filamentous haemagglutinin and its secretion machinery. Microbiology. 2000, 146 (Pt 5): 1211-1221.
Kimura A, Mountzouros KT, Relman DA, Falkow S, Cowell JL: Bordetella pertussis filamentous hemagglutinin: evaluation as a protective antigen and colonization factor in a mouse respiratory infection model. Infect Immun. 1990, 58: 7-16.
Relman DA, Domenighini M, Tuomanen E, Rappuoli R, Falkow S: Filamentous hemagglutinin of Bordetella pertussis: nucleotide sequence and crucial role in adherence. Proc Natl Acad Sci USA. 1989, 86: 2637-2641. 10.1073/pnas.86.8.2637.
Cotter PA, Yuk MH, Mattoo S, Akerley BJ, Boschwitz J, Relman DA, Miller JF: Filamentous hemagglutinin of Bordetella bronchiseptica is required for efficient establishment of tracheal colonization. Infect Immun. 1998, 66: 5921-5929.
Abramson T, Kedem H, Relman DA: Proinflammatory and proapoptotic activities associated with Bordetella pertussis filamentous hemagglutinin. Infect Immun. 2001, 69: 2650-2658. 10.1128/IAI.69.4.2650-2658.2001.
Inatsuka CS, Julio SM, Cotter PA: Bordetella filamentous hemagglutinin plays a critical role in immunomodulation, suggesting a mechanism for host specificity. Proc Natl Acad Sci USA. 2005, 102: 18578-18583. 10.1073/pnas.0507910102.
Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, Holden MT, Churcher CM, Bentley SD, Mungall KL: Comparative analysis of the genome sequences of Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica. Nat Genet. 2003, 35: 32-40. 10.1038/ng1227.
Cummings CA, Brinig MM, Lepp PW, Pas van de S, Relman DA: Bordetella species are distinguished by patterns of substantial gene loss and host adaptation. J Bacteriol. 2004, 186: 1484-1492. 10.1128/JB.186.5.1484-1492.2004.
King AJ, van Gorkom T, Pennings JL, Heide van der HG, He Q, Diavatopoulos D, Heuvelman K, van Gent M, van Leeuwen K, Mooi FR: Comparative genomic profiling of Dutch clinical Bordetella pertussis isolates using DNA microarrays: identification of genes absent from epidemic strains. BMC Genomics. 2008, 9: 311-10.1186/1471-2164-9-311.
Caro V, Hot D, Guigon G, Hubans C, Arrive M, Soubigou G, Renauld-Mongenie G, Antoine R, Locht C, Lemoine Y, Guiso N: Temporal analysis of French Bordetella pertussis isolates by comparative whole-genome hybridization. Microbes Infect. 2006, 8: 2228-2235. 10.1016/j.micinf.2006.04.014.
Julio SM, Cotter PA: Characterization of the filamentous hemagglutinin-like protein FhaS in Bordetella bronchiseptica. Infect Immun. 2005, 73: 4960-4971. 10.1128/IAI.73.8.4960-4971.2005.
Kosoy MY, Regnery RL, Kosaya OI, Childs JE: Experimental infection of cotton rats with three naturally occurring Bartonella species. J Wildl Dis. 1999, 35: 275-284.
Schulein R, Dehio C: The VirB/VirD4 type IV secretion system of Bartonella is essential for establishing intraerythrocytic infection. Mol Microbiol. 2002, 46: 1053-1067. 10.1046/j.1365-2958.2002.03208.x.
Schulein R, Guye P, Rhomberg TA, Schmid MC, Schroder G, Vergunst AC, Carena I, Dehio C: A bipartite signal mediates the transfer of type IV secretion substrates of Bartonella henselae into human cells. Proc Natl Acad Sci USA. 2005, 102: 856-861. 10.1073/pnas.0406796102.
Schmid MC, Scheidegger F, Dehio M, Balmelle-Devaux N, Schulein R, Guye P, Chennakesava CS, Biedermann B, Dehio C: A translocated bacterial protein protects vascular endothelial cells from apoptosis. PLoS Pathog. 2006, 2: e115-10.1371/journal.ppat.0020115.
Kosoy MY, Saito EK, Green D, Marston EL, Jones DC, Childs JE: Experimental evidence of host specificity of Bartonella infection in rodents. Comp Immunol Microbiol Infect Dis. 2000, 23: 221-238. 10.1016/S0147-9571(99)00075-2.
Steppan S, Adkins R, Anderson J: Phylogeny and divergence-date estimates of rapid radiations in muroid rodents based on multiple nuclear genes. Syst Biol. 2004, 53: 533-553. 10.1080/10635150490468701.
Birtles RJ, Raoult D: Comparison of partial citrate synthase gene (gltA) sequences for phylogenetic analysis of Bartonella species. Int J Syst Bacteriol. 1996, 46: 891-897.
Norman AF, Regnery R, Jameson P, Greene C, Krause DC: Differentiation of Bartonella-like isolates at the species level by PCR-restriction fragment length polymorphism in the citrate synthase gene. J Clin Microbiol. 1995, 33: 1797-1803.
Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.
Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.
Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8: 195-202.
Riess T, Dietrich F, Schmidt KV, Kaiser PO, Schwarz H, Schafer A, Kempf VA: Analysis of a novel insect cell culture medium-based growth medium for Bartonella species. Appl Environ Microbiol. 2008, 74: 5224-5227. 10.1128/AEM.00621-08.
Ehrenborg C, Wesslen L, Jakobson A, Friman G, Holmberg M: Sequence variation in the ftsZ gene of Bartonella henselae isolates and clinical samples. J Clin Microbiol. 2000, 38: 682-687.
We thank Kirsten Ellegaard, Anna Åsman and Marcel Verhaart for technical assistance, Magnus Lundgren for graphical work and Björn Nystedt for helpful discussions. We thank Annelie Waldén at the Royal Institute of Technology, Stockholm, Sweden, and Niclas Olsson at Uppsala University, Uppsala, Sweden, for microarray printing. This work was supported by grants to S.G.E.A. from the Swedish Research Council, the Göran Gustafsson Foundation, the Swedish Foundation for Strategic Research and the Knut and Alice Wallenberg Foundation.
ECB performed the microarray CGH experiments, carried out the sequence and microarray data analyses and prepared the manuscript. CE participated in the sequence analysis and drafted part of the manuscript. OVP performed the PFGE experiments. FG performed the cultivation and aggregation experiments. KN performed the cultivation and sequencing experiments. MH contributed strains, analyzed experimental data and helped in the design of the study. SGEA designed the study, analyzed experimental data, prepared the manuscript and coordinated the project. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
About this article
Cite this article
Berglund, E.C., Ehrenborg, C., Vinnere Pettersson, O. et al. Genome dynamics of Bartonella grahamii in micro-populations of woodland rodents. BMC Genomics 11, 152 (2010). https://doi.org/10.1186/1471-2164-11-152