Classification of Bartonella Strains Associated with Straw-Colored Fruit Bats (Eidolon helvum) across Africa Using a Multi-locus Sequence Typing Platform

Ying Bai; David T. S. Hayman; Clifton D. McKee; Michael Y. Kosoy

doi:10.1371/journal.pntd.0003478

Abstract

Bartonellae are facultative intracellular bacteria and are highly adapted to their mammalian host cell niches. Straw-colored fruit bats (Eidolon helvum) are commonly infected with several bartonella strains. To elucidate the genetic diversity of these bartonella strains, we analyzed 79 bartonella isolates from straw-colored fruit bats in seven countries across Africa (Cameroon, Annobon island of Equatorial Guinea, Ghana, Kenya, Nigeria, Tanzania, and Uganda) using a multi-locus sequencing typing (MLST) approach based on nucleotide sequences of eight loci (ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA). The analysis of each locus but ribC demonstrated clustering of the isolates into six genogroups (E1 – E5 and Ew), while ribC was absent in the isolates belonging to the genogroup Ew. In general, grouping of all isolates by each locus was mutually supportive; however, nuoG, gltA, and rpoB showed some incongruity with other loci in several strains, suggesting a possibility of recombination events, which were confirmed by network analyses and recombination/mutation rate ratio (r/m) estimations. The MLST scheme revealed 45 unique sequence types (ST1 – 45) among the analyzed bartonella isolates. Phylogenetic analysis of concatenated sequences supported the discrimination of six phylogenetic lineages (E1 – E5 and Ew) corresponding to separate and unique Bartonella species. One of the defined lineages, Ew, consisted of only two STs (ST1 and ST2), and comprised more than one-quarter of the analyzed isolates, while other lineages contained higher numbers of STs with a smaller number of isolates belonging to each lineage. The low number of allelic polymorphisms of isolates belonging to Ew suggests a more recent origin for this species. Our findings suggest that at least six Bartonella species are associated with straw-colored fruit bats, and that distinct STs can be found across the distribution of this bat species, including in populations of bats which are genetically distinct.

Author Summary

Bats, with over 1000 recognized species, represent about 20% of all classified mammalian species worldwide. These mammals have a wide range of ecologies and life-history traits, and are now widely recognized as important reservoirs of many pathogens. Bartonella species have been found distributed in a wide range of mammalian species, including bats. About half of recognized Bartonella species, including one bat-associated species, have been associated with human illness. Previous studies have shown that Bartonella species are extremely diverse, with or without evident specificity to their mammalian hosts. Possessing many unique aspects, bartonellae can serve as a useful biological marker to study how microorganisms have evolved and diversified along with their animal hosts in evolutionary history. In this study, we applied multi-locus sequence typing, or MLST, to study the genetic differences of straw-colored fruit bat (Eidolon helvum)-associated Bartonella species. Our studies suggest Bartonella species have both exchanged genetic materials among species through recombination events and lost genes that are perhaps superfluous to their life cycles, which includes an intracellular stage in mammals.

Citation: Bai Y, Hayman DTS, McKee CD, Kosoy MY (2015) Classification of Bartonella Strains Associated with Straw-Colored Fruit Bats (Eidolon helvum) across Africa Using a Multi-locus Sequence Typing Platform. PLoS Negl Trop Dis 9(1): e0003478. https://doi.org/10.1371/journal.pntd.0003478

Editor: David H. Walker, University of Texas Medical Branch, UNITED STATES

Received: October 16, 2014; Accepted: December 13, 2014; Published: January 30, 2015

This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: The authors received no specific funding for this work.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Bartonellae are both Gram-negative alpha-proteobacteria and hemotropic bacteria highly adapted to facultative intracellular lifestyle in a wide variety of mammals, such as rodents, bats, insectivores, carnivores, ungulates, and other vertebrates. During the last two decades, progressively more bacterial species belonging to the genus Bartonella have been recognized with over 30 species described from different mammalian hosts. A number of Bartonella species were found to be associated with human illnesses and are associated with a growing spectrum of emerging diseases, including life-threatening endocarditis [1–9]. Animal reservoirs have been identified for some of the human pathogens, while remain unknown for others. Knowledge of the transmission of bartonella bacteria between mammalian hosts is incomplete. However, hematophagous arthropods such as fleas, flies, lice, mites, and ticks have been found naturally infected and are frequently implicated in transmitting Bartonella species [10–15].

Increasing recognition of bats as natural reservoirs of many emerging pathogens has drawn considerable attentions to study these mammals [16]. Multiple investigations of bartonella in bats have been conducted in different regions of the world [17–23]. These studies reported that bartonella infections are highly prevalent in many bat species and bartonella communities associated with bats are extremely diverse with co-circulation of numerous Bartonella species in the same bat populations. The data on relationships between Bartonella species and bats are quite contradictive depending on the investigated geographic region. Specifically, in Central and South America, multiple bat species may share the same Bartonella species without an evident host-specificity [19,20], while investigations from Asia and Africa demonstrated that a bat population typically harbors one or few Bartonella species specific for a particular bat species [17–18,21–22].

The straw-colored fruit bat (Eidolon helvum) is widely distributed in Africa. Colonies of these bats can be large and often found near human populations, with some roosts containing millions of individuals. Like other bats, straw-colored fruit bats have long life spans, with some individuals reaching at least 14 years of age [24]. Local human residents have a close contact with the bats, as they are frequently hunted for “bush meat”. A previous study in Kenya showed that straw-colored fruit bats are infected with bartonellae with high prevalence reaching 26% [18]. The same study demonstrated that Bartonella associated with the bats were genetically distant and belonged to four distinct genogroups based on sequence variation in the citrate synthase gene (gltA) [18]. Similar observations were reported in a very recent study conducted in Nigeria [22]. Considering the broad distribution and other ecological characteristics of straw-colored fruit bats, it was expected that additional Bartonella genogroups associated with this bat species can be identified in Africa.

In the present study, we aim to expand our knowledge of bartonella infections in straw-colored fruit bats and to better understand how multiple Bartonella species can co-habit populations of one bat species. We compared genetic differences of bartonella isolates obtained from straw-colored fruit bats captured in seven countries across Africa using a multi-locus sequence typing (MLST) approach. Based on comparison of nucleotide sequences derived from multiple loci, MLST has been shown to provide high discriminatory power in epidemiological and genetic analysis of bacterial strain populations, while retaining signatures of longer-term evolutionary relationships or clonal stability [25–28]. In addition, sequencing of multiple loci can detect evidences of micro-evolutionary events, particularly homologous recombination, among identified sequence types [29,30]. The main objective of the current study was to analyze bartonella isolates obtained from naturally infected straw-colored fruit bats from different parts of Africa to determine whether well-defined phylogenetic lineages correspond to currently accepted criteria for discrimination of bacterial species [31]. This, in turn, can help to enhance our understanding of population structure in the bacteria and the relationships between allelic profiles and the animal host. To achieve this goal, we developed a MLST scheme that incorporates eight loci which have been previously proposed for characterization of Bartonella species. In this study, clusters defined by the phylogenetic analysis were called either a “genogroup” or a “lineage”. A genogroup is defined when a single locus sequence was applied for the phylogenetic analysis, while a lineage is defined from concatenated sequences of all loci.

Materials and Methods

Bartonella isolates selection

A total of 79 isolates (Table 1) originating from straw-colored fruit bats from seven countries in different regions across Africa were selected for the MLST analysis, including Annobón island of Equatorial Guinea (n = 9), Cameroon (n = 3). Ghana (n = 20), Kenya (n = 16), Nigeria (n = 11), Tanzania (n = 15), and Uganda (n = 5). These isolates were obtained from previous studies completed by Bartonella Laboratory, CDC at Fort Collins, Colorado, and culturing procedures were published elsewhere [18,20].

Download:

Table 1. Allelic profiles, sequence types (ST), and lineage identification for the 79 Bartonella isolates from different geographic locations.

https://doi.org/10.1371/journal.pntd.0003478.t001

Multi-locus sequence typing (MLST)

Eight loci (ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA) that have been previously used for bartonella description [27] were selected for MLST characterization of the bartonella isolates. Information on primers is provided in Table 2. A specific fragment of each locus was amplified by PCR for each of the 79 isolates. PCR products of each locus were purified with the QIAquick PCR Purification Kit (Qiagen, Germantown, MD) and sequenced in both directions using an Applied Biosystems Model 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA). Using the Lasergene software package (DNASTAR, Madison, WI), obtained sequences were aligned by each locus and compared between the isolates and with other bat-originated Bartonella strains (S1 Table) and known Bartonella species (S2 Table). Based on the allelic profile (Table 1), each unique variant was designated as a sequence type (ST) and sequences for the eight loci were concatenated.

Download:

Table 2. Characteristics of the eight loci evaluated for the straw-colored fruit bat (Eidolon helvum)-associated Bartonella MLST scheme.

https://doi.org/10.1371/journal.pntd.0003478.t002

Analysis of nucleotide polymorphism and diversity

Sequence polymorphisms of the eight MLST loci were examined in MEGA v6.06 separately for each Bartonella lineage and among all 45 STs. For protein-coding loci (ftsZ, gltA, nuoG, ribC, and rpoB), the appropriate open reading frame that contained no stop codons was selected for each locus by checking all possible starting positions. To test for evidence of selection in these loci, d_N/d_S ratios were calculated using the Nei-Gojobori method for each of the loci considering each lineage separately and for all 45 STs together. Nucleotide diversity (π), determined by the number of nucleotide differences per site, was calculated as p-distance for all eight loci individually and for concatenated sequences considering each lineage separately and for all 45 STs together. For concatenated sequences, the total length included gap regions and missing genes, thus the pairwise deletion option was used.

Phylogenetic analysis

A neighbor-joining tree based on the concatenated MLST alleles alone was constructed using the Clustal W program within the Lasergene 11 package of DNASTAR (version 8). An additional phylogeny was constructed, in which known Bartonella species and Bartonella strains isolated from bats in Old World were included in addition to representative strains from straw-colored fruit bats (E1—E5, Ew), and Brucella abortus was included as an outgroup, This phylogeny was inferred using sequences of ftsZ, gltA, nuoG, ribC, rpoB, ssrA, and 16S rRNA; ITS sequences were not included due to the large number of gaps among the strains that could not be resolved. Sequences from each locus were aligned using Clustal X v2.1, trimmed to equal lengths, and concatenated. The best model of nucleotide substitution was determined using MEGA. Based on this model, a maximum-likelihood tree was generated in MEGA with 1000 bootstrap replicates. Due to missing genes among the strains as well as some alignment gaps, we used the pairwise deletion option when inferring the tree.

Recombination tests

To visualize potential recombination events among the sequence types, a phylogenetic network was inferred from concatenated sequences from the 45 STs using the Neighbor-Net algorithm in SplitsTree v4.13.1 with 1000 bootstrap replicates. Gaps in aligned sequences and from missing genes were included. The pairwise homoplasy index (PHI) was implemented in SplitsTree to test for significant recombination among the isolates.

ClonalFrame v1.1 was used to estimate the relative contribution of recombination and mutation in generating polymorphisms among the 45 Bartonella STs. Based on recommendations by Vos and Didelot [32], ITS and genes coding for RNA (ssrA and 16S rRNA) were not included because of potential confounding effects of selection on the detection of homologous recombination rates. Therefore, we used only the five protein-coding loci (ftsZ, gltA, nuoG, ribC, and rpoB) in the study. Two independent runs were performed in ClonalFrame using 200,000 MCMC iterations. The initial mutation rate (θ) was set using Watterson's moment estimator while the remaining initial parameters used the default settings in ClonalFrame. We assessed convergence and mixing properties of the dataset through visual inspection of the traces for the likelihood and model parameters.

Nucleotide sequence accession numbers

Newly-identified alleles from the current study were submitted to GenBank with the following accession numbers: KJ999677 to KJ999694 (ftsZ), KM030503 to KM030526 (gltA), KM030527 to KM030544 (nuoG), KM215178 to KM215188 (ribC), KM215189 to KM215208 (rpoB), KM233456 to KM233470 (ssrA), KM233471 to KM233471 (ITS), and KJ944271 to KJ944279 (16S rRNA).

Results

Individual sequence analysis of gltA and other loci

All 79 isolates were first compared based on the 751bp gltA fragment for the initial identification. The gltA sequence alignment revealed 247 variable sites with 24 variants, delineating six unique genogroups with divergence of 7.3–23% among the genogroups and with similarities of ≥ 96.8% within a genogroup. Four genogroups were the same as previously identified in straw-colored fruit bats from Kenya, which were named E1, E2, E3, and Ew [18]; whereas the remaining two genogroups were distinct from these four and from any other previously reported Bartonella genotypes. Continuing the same naming system proposed for discrimination of Bartonella strains discovered in Kenya [18], these two new genogroups were designated as E4 and E5.

Following the initial identification, additional analyses were performed with the other seven loci. Each analyzed locus, except ribC, showed that all bartonella isolates obtained from straw-colored fruit bats also fell into the same six genogroups identified by gltA. Unexpectedly, analysis of the ribC locus revealed a large variation in length of the examined fragment among the isolates. Depending on the group identified by other markers, the examined ribC fragment was either fully presented (535bp), or reduced (382bp), or non-amplifiable (PCR negative). The phylogenetic analysis based on ribC sequences of the 79 isolates revealed five genogroups related to E1—E5 clusters that correspond to the grouping by other loci. Absence of the ribC locus was indicative for the genogroup Ew. All isolates with a reduced ribC fragment (partially missing) belonged to the genogroup E5 (Table 1).

Allelic profiles, sequence types (ST), and phylogenetic analysis

The size of sequenced fragments ranged between 282bp and 1172bp at different loci. The investigated loci showed different degrees of variation, with 38–249 variable sites and 9–24 alleles (Table 2). The length of concatenated sequences ranged from 4,622bp to 5,160bp as a result of the variation in the fragment length of ribC (from 0 to 535bp), ITS (from 315 to 352bp), and ssrA (from 282 to 289bp). Based on the allelic profiles, the MLST analysis distinguished 45 sequence types (ST) among the 79 isolates, showing high heterogeneity. Most STs were represented by a single isolate, while some were represented by 2–4 isolates, and ST1 was found in 17 isolates (Table 1). Phylogenetic analysis demonstrated that all of the 45 STs resolved into six lineages, namely, E1, E2, E3, E4, E5, and Ew, to match the names proposed for grouping the strains based on sequences of one locus (Fig. 1). The divergence was 5.1–14% among lineages and ≤ 4.2% within a lineage.

Download:

Figure 1. Phylogenetic relationships of the 45 sequence types from 79 bartonella isolates obtained from straw-colored fruit bats (Eidolon helvum) in seven African countries/regions.

The number of isolates belonging to each sequence type is given in parentheses. The phylogenetic tree was constructed from concatenated sequences (4,622bp—5,160bp) of eight loci (ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA) using the neighbor-joining method. Bootstrap values were calculated with 1000 replicates. The sequence types are grouped into six phylogenetic lineages (boxed clades) named as E1—E5 and Ew, with each lineage presumably representing a separate and unique Bartonella species.

https://doi.org/10.1371/journal.pntd.0003478.g001

Lineage Ew contained 21 isolates of two similar sequence types (ST1 and ST2) with only one nucleotide difference between the two STs. The 17 isolates presenting ST1 were discovered in each of the seven studied countries. ST2 was identified in four isolates from Ghana, Kenya, Nigeria, and Tanzania (one isolate per country) (Table 1). None of the isolates belonging to this lineage possessed the fragment for ribC. The lineage Ew was very common among the six lineages identified and accounted for 26.6% (21/79) of all analyzed isolates.

Lineage E1 contained 13 isolates of 10 sequence types (ST12—ST21). The distance among the STs was less than 2.1%. Eight isolates were from Kenya, and the others were from Ghana, Nigeria, and Tanzania (Table 1).

Lineage E2 included seven isolates of six sequence types (ST22—ST27). The distance among the STs was less than 2.1%. Isolates belonging to this lineage were from Kenya, Ghana, Tanzania, and Uganda (Table 1).

Lineage E3 contained 24 isolates of 15 sequence types (ST28—ST42). The distance among the STs was less than 4.2%. ST30, ST36, ST37, and ST40 were represented in four, three, two, and three isolates, respectively. The remaining STs were each represented by a single strain.

Lineage E3 was the most common among all bartonella lineages detected in straw-colored fruit bats, and accounted for 30.4% (26/79) of all isolates analyzed. Notably, mismatches in assignment of isolates to specific ST lineages were observed for a few loci on several occasions. Specifically, isolate B32119 (ST38) from a Kenyan bat and isolate B23797 (ST39) from a Nigerian bat were assigned to lineage E3 by ST classification, but separate analyses of the gltA and nuoG sequences of the two isolates indicated their closeness to genogroup Ew. Similarly, isolate B40005 (ST28) from a Cameroonian bat was identified as E3 by ST classification, but was closer to genogroup E1 by the gltA and nuoG sequences.

Lineage E4 contained four isolates of three sequence type (ST43—ST45). These STs are very similar among themselves with distance ≤ 0.1%. The isolates were from Ghana, Tanzania, and Uganda (Table 1).

Lineage E5 contained 10 isolates of eight sequence types (ST3—ST11). Compared to other lineages, the STs within this lineage were more distant (0.2–3.5%). All isolates/STs in this lineage presented shorter concatenated sequences due to the partial fragment missing in ribC. One isolate (B39286) from this lineage also showed a mismatch in lineage assignment when analyzed by individual loci. The isolate was from a Ghanaian bat and belonged to lineage E5 by ST classification, but assigned to genogroup E2 by rpoB and nuoG sequences.

Based on MEGA results, the best nucleotide substitution model for the concatenated sequences of gltA, ftsZ, nuoG, ribC, rpoB, ssrA, and 16S rRNA was determined to be GTR+G+I [33]. Comparison with known Bartonella species and other bat-associated Bartonella strains demonstrated that lineages E1, E2, E3, and E5 belong to a clade that includes other Bartonella strains found in Old World bats, while E4 and Ew appear to have arisen independently (Fig. 2).

Download:

Figure 2. Phylogenetic relationship of straw-colored fruit bats (E. helvum) with known Bartonella species and other bat-associated Bartonella strains.

The maximum-likelihood tree was inferred using concatenated sequence of seven loci (ftsZ, gltA, nuoG, ribC, rpoB, ssrA, and 16S rRNA) from 31 Bartonella species, other bat-associated Bartonella strains, and 6 sequence types from E. helvum representing the lineages E1—E5 and Ew (bold text). Groups of bat-associated Bartonella strains are indicated with a bat silhouette. The numbers at the nodes correspond to bootstrap values greater than 60% based on 1000 replicates.

https://doi.org/10.1371/journal.pntd.0003478.g002

Patterns of selection and diversity in nucleotide sequences

The d_N/d_S ratios calculated for protein-coding sequences ranged from 0.022 for ftsZ to 0.242 for ribC when comparing all 45 STs. The values for each locus differed when each Bartonella lineage was analyzed separately. Nucleotide diversity among STs varied by locus, ranging from 1.02% in 16S rRNA to 14.38% in ribC. Among the Bartonella lineages, E2 and E5 had the highest nucleotide diversity when all loci were considered together, with 1.56% and 1.24%, respectively, while Ew had the lowest (0.02%).

Recombination test

A network phylogeny based on the concatenated sequences for the 45 STs in SplitsTree (v4.13.1) shows that the 45 STs fall into the same six identified lineages identified by single loci. However, isolates B32119 from Nigeria, B23797 from Kenya, B40005 from Cameroon, and B39286 from Ghana are split from the network, reflecting their mixed ancestry (Fig. 3). The pairwise homoplasy index (PHI) [34] found significant evidence of recombination in the network (mean = 0.2, variance = 4.5x10^–6, p < 0.01). Using ClonalFrame, the r/m value (the ratio of probabilities that a site is altered by recombination or mutation) was estimated to be 1.14 (95% CI 0.63, 1.82) and 1.06 (95% CI 0.53, 1.67) based on two independent runs (average = 1.1), while the ρ/θ value (the ratio of rates at which recombination and mutation occur) was estimated to be 0.049 (95% CI 0.023, 0.087) and 0.037 (95% CI 0.016, 0.065).

Download:

Figure 3. Network phylogeny of the 45 bartonella sequence types obtained from Eidolon helvum.

The network was constructed in SplitsTree using the NeighborNet algorithm based on concatenated sequences of eight loci (ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA). Clusters of sequence types were named according to phylogenetic lineages (E1—E5, Ew). Individual isolate labels indicate samples with mixed ancestry due to possible recombination.

https://doi.org/10.1371/journal.pntd.0003478.g003

Discussion

Straw-colored fruit bats have been reported to be commonly infected with multiple Bartonella strains based on gltA sequences [18,22]. In this study, we developed a MLST scheme using eight loci which have been commonly used for Bartonella species description to characterize 79 bartonella isolates obtained from straw-colored fruit bats that were captured in different regions across Africa to better elucidate the genetic diversity among these isolates.

The analyses of selection pattern showed none of the d_N/d_S ratios exceeded 1. This suggested that all of the protein-coding loci in the study are undergoing purifying selection. The analyzed bartonella isolates exhibited a high level of heterogeneity with 9–24 alleles identified by different loci. The MLST scheme resolved 45 STs among the 79 isolates and these isolates/STs clustered into six distinct lineages (E1—E5, and Ew) with genetic distances of 5.1–14% among the lineages. MLST has been regarded as a highly discriminatory tool by a number of previous studies. According to the widely accepted criteria for a new bacterial species based on nucleotide divergence greater than 5% in housekeeping genes [31], the MLST analyses in the present study demonstrated that each lineage identified in straw-colored fruit bats in this study likely represents a new Bartonella species. These Bartonella species have never been found in other bat species, suggesting their specific co-adaptation to Eidolon bats, though future studies are necessary to confirm this. The performed analyses allow us to conclude that there are at least six Bartonella species associated with straw-colored fruit bats. Phylogenetic analysis based on multiple loci demonstrates that four of these lineages, E1, E2, E3, and E5, belong to a phylogenetic lineage that includes other Bartonella strains isolated from Old World bats [18,21–23], while E4 and Ew appear to have evolved independently.

Lineages E3 and Ew, containing 24 and 21 isolates, respectively, and accounting for 57% of all isolates, are likely to be the more common Bartonella species circulating among this bat species. However, because the original bat samples from which the analyzed bartonella isolates obtained were collected in different settings and during different seasons by multiple investigators, there could be some bias in temporal or seasonal variations. Furthermore, the culturing procedure that was necessary for performing such bacterial classification may also contribute to misrepresentation of comparative prevalence of Bartonella species as some of them could hardly be cultured [35]. Due to the limited number of isolates selected from each part of Africa, we did not attempt to investigate any patterns in the geographic distribution of Bartonella species in this study. The presented classification, however, can be used for future studies aimed at describing geographic patterns of bartonella distribution among Eidolon bats.

Most lineages consisted of multiple STs (from 4 to 14). Nevertheless, lineage Ew consisted of only two STs with a large number of isolates, and in fact there is only one nucleotide difference between the two STs. This indicates a lower diversity of allelic polymorphisms in the isolates from this lineage and allows a speculation that Ew has diverged more recently compared to the other Bartonella lineages. Further, our observations of partially missing or non-amplifiable sequences of ribC suggests gene loss, a common evolutionary process in bacteria [36], possibly has occurred and resulted in the appearance of two new Bartonella species (E5 and Ew). As a housekeeping gene, ribC participates in the regulation of riboflavin biosynthesis. The non-amplifiability of the gene might not be true but suggest the gene is in some undetected form. If this gene is truly missing, its loss may be compensated by the fact that bartonellae are intracellular parasites which can perhaps utilize host riboflavin or they have lost the need of it. If this were true, fitness costs may lead to selection against this superfluous genetic material [37]. Overall, little is known about the mechanisms of regulation of bacterial riboflavin genes [38]. Although we cannot explain biological consequences of the loss of ribC, the 'presence/absence' characteristics make ribC a marker for the differentiation and classification of Eidolon-associated Bartonella strains. In addition, the dramatic variation in ribC among the straw-colored fruit bat-associated Bartonella species presents an interesting model for understanding evolutionary trends in developing host-bacterial parasite relations [36].

Another intriguing observation in this study was the set of mismatches between MLST lineages and some genogroups comparing phylograms built based either on entire concatenated sequences or individual locus sequences for four isolates. Particularly, three isolates belonging to E3 based on the ST classification were either Ew (two isolates) or E1 (one isolate) by nuoG or gltA alone and one isolate of E5 by ST classification would be E2 by nuoG or rpoB alone. The loci rpoB and gltA have been shown to be the most potent for species discrimination in other studies [31]. According to our observations in the present study, an analysis limited to two loci may still lead to a misclassification. Thus characterization of multiple loci should be preferably used for Bartonella species identification. The mismatches between lineages and genogroups identified by different loci indicate that recombination and/or lateral gene transfer events between different Bartonella species are ongoing processes.

Our recombination test results indicate that while recombination happens infrequently, an individual recombination event introduces a large number of polymorphisms. Recombination in rpoB, gltA, and other housekeeping loci has been observed in other Bartonella species [29,30,39]. Estimated r/m values for B. grahamii alone (1.7 [38] and 6.81 [30]), B. taylorii alone (3.77 [27]), or B. grahamii and B. taylorii considered together (4.06 [30]) compare favorably with our estimate for recombination among six species of Bartonella in E. helvum (r/m = 1.1). These estimates contrast with results from B. henselae (r/m = 0.1) [40] and the perception that intracellular bacteria have low rates of recombination and horizontal gene transfer [32]. Intermediate to high rates of recombination in these studies suggest that multiple Bartonella species may coexist at some point in the infection cycle, potentially in their mammalian hosts and/or arthropod vectors, thus facilitating gene exchange. Moreover, gene exchange appears to play an important role in generating sequence diversity among co-circulating Bartonella species. Further studies that could measure co-infection of Bartonella species and resulting rates of recombination would help to explain these interesting dynamics.

Our study provides information on the genetic diversity of straw-colored fruit bat-associated Bartonella species that can be used for the discrimination of the lineages at the level corresponding to separate species. Additional information is needed to understand the interaction between these Bartonella species and their bat hosts under natural conditions, as well as their putative ectoparasite vectors [41]. A biological approach to define Bartonella species has been discussed [42] with emphasis on host-association as a phenomenon that promotes biological isolation of Bartonella species. However, such an approach has limitations for situations when several Bartonella species are associated with one animal host. A similar situation was also described among Bartonella strains in grasshopper mice [43] and cotton rats [44]. Nevertheless, among all mammalian species described as a host of one or multiple Bartonella species, the number of presumptive species associated with E. helvum bats is clearly the highest. Although there is no evidence suggesting that Bartonella species associated with African fruit bats may cause illnesses in humans or in bats themselves, a recent study reported the presence of Bartonella mayotimonensis, a reported etiologic agent of endocarditis in humans, in the Daubenton's bat (Myotis daubentonii) and the Northern bat (Eptesicus nilssonii) in Europe [23], suggesting a potential role of bats as reservoirs for human bacterial pathogens. It is important to understand the role of straw-colored fruit bats in biotic communities and their importance as reservoir hosts contributing to the maintenance and transmission of bartonellae to other animals and humans.

The great diversity of Bartonella species associated with E. helvum, however, clearly presents interesting questions from the perspective of microbial evolution. What are the mechanisms that generate and maintain such diversity? How can these bacteria share the same ecological niche? Analyzing such relationships, Chan and Kosoy [45] hypothesized that co-circulation of independent Bartonella species in populations of one host species may represent an escape mechanism to circumvent the host immune responses. We ultimately believe this system will be very informative for understanding the evolution of bacteria in a host-vector-pathogen system.

Supporting Information

S1 Table. GenBank accession numbers for ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA sequences for bat-associated Bartonella strains.

https://doi.org/10.1371/journal.pntd.0003478.s001

(DOCX)

S2 Table. GenBank accession numbers for ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA sequences of reference Bartonella species.

https://doi.org/10.1371/journal.pntd.0003478.s002

(DOCX)

Acknowledgments

We are indebted to many people who collected blood samples from bats working in the field. We wish specifically to thank Alison Peel, Richard Suu-Ire, and Ivan Kuzmin for their assistance.

Author Contributions

Conceived and designed the experiments: YB MYK DTSH. Performed the experiments: YB. Analyzed the data: YB CDM. Contributed reagents/materials/analysis tools: MYK. Wrote the paper: YB CDM MYK DTSH.

References

1. Welch DF, Pickett DA, Slater LN, Steigerwalt AG, Brenner DJ (1992) Rochalimaea henselae sp. nov., a cause of septicemia, bacillary angiomatosis, and parenchymal bacillary peliosis. J Clin Microbiol 30: 275–280. pmid:1537892
- View Article
- PubMed/NCBI
- Google Scholar
2. Daly JS, Worthington MG, Brenner DJ, Moss CW, Hollis DG, et al (1993) Rochalimaea elizabethae sp. nov. isolated from a patient with endocarditis. J Clin Microbiol 31: 872–881. pmid:7681847
- View Article
- PubMed/NCBI
- Google Scholar
3. Bass JW, Vincent JM, Person DA (1997) The expanding spectrum of bartonella infections II. Cat scratch disease. Pediatr Infect Dis J 16: 163–179. pmid:9041596
- View Article
- PubMed/NCBI
- Google Scholar
4. Kerkhoff FT, Bergmans AM, van Der Zee A, Rothova A (1999) Demonstration of Bartonella grahamii DNA in ocular fluids of patient with neuroretinitis. J Clin Microbiol 37: 4034–4038. pmid:10565926
- View Article
- PubMed/NCBI
- Google Scholar
5. Roux V, Eykyn SJ, Wyllie S, Raoult D (2000) Bartonella vinsonii subsp. berkhoffii as an agent of a febrile blood culture-negative endocarditis in a human. J Clin Microbiol 38: 1698–1700. pmid:10747175
- View Article
- PubMed/NCBI
- Google Scholar
6. Kosoy MY, Murray M, Gilmore RD, Bai Y, Gage K (2003) Bartonella strains from ground squirrels are identical to Bartonella washoensis isolated from a human patient. J Clin Microbiol 41: 645–650. pmid:12574261
- View Article
- PubMed/NCBI
- Google Scholar
7. Eremeeva ME, Gerns HL, Lydy SL, Goo JS, Ryan ET, et al (2007) Bacteremia, fever, and splenomegaly caused by a newly recognized Bartonella species. N Engl J Med 256: 2381–2387.
- View Article
- Google Scholar
8. Kosoy M, Morway C, Sheff K, Bai Y, Colborn J, et al (2008) Bartonella tamiae sp. nov., a newly recognized pathogen isolated from human patients from Thailand. J Clin Microbiol 46: 772–775. pmid:18077632
- View Article
- PubMed/NCBI
- Google Scholar
9. Welch DF, Carroll KC, Hofmeister EK, Persing DH, Robison DA, et al (1999) 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 37: 2598–2601. pmid:10405408
- View Article
- PubMed/NCBI
- Google Scholar
10. Garcia-Caceres U, Garcia FU (1991) Bartonellosis: an immunodepressive disease and the life of Daniel Alcides Carrion. Am J Clin Pathol 95: s58–s66. pmid:2008885
- View Article
- PubMed/NCBI
- Google Scholar
11. Chomel BB, Kasten RW, Floyd-Hawkins KA, Chi B, Yamamoto K, et al (1996). Experimental transmission of Bartonella henselae by the cat flea. J Clin Microbiol 34: 1952–1956. pmid:8818889
- View Article
- PubMed/NCBI
- Google Scholar
12. Higgins JA, Radulovic S, Jaworski DC, Azad AF (1996) Acquisition of the cat scratch disease agent Bartonella henselae by cat fleas (Siphonaptera: Pulicidae). J Med Entomol 33: 490–495. pmid:8667399
- View Article
- PubMed/NCBI
- Google Scholar
13. Pappalardo BL, Correa MT, York CC, Peat CY, Breitschwerdt EB (1997) Epidemiologic evaluation of the risk factors associated with exposure and seroreactivity to Bartonella vinsonii in dogs. Am J Vet Res 58: 467–471. pmid:9140552
- View Article
- PubMed/NCBI
- Google Scholar
14. Roux V, Raoult D (1999) Body lice as tools for diagnosis and surveillance of reemerging diseases. J Clin Microbiol 37: 596–599. pmid:9986818
- View Article
- PubMed/NCBI
- Google Scholar
15. Billeter SA, Sangmaneedet S, Kosakewich RC, Kosoy MY (2012) Bartonella species in dogs and their ectoparasites from Khon Kaen Province, Thailand. Southeast Asian J Trop Med Public Health 43: 1186–1192. pmid:23431825
- View Article
- PubMed/NCBI
- Google Scholar
16. Luis AD, Hayman DTS, O’Shea TJ, Cryan PM, Gilbert AT, et al (2013) A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special? Proc Biol Sci 280: 20122753. pmid:23378666
- View Article
- PubMed/NCBI
- Google Scholar
17. Concannon R, Wynn-Owen K, Simpson VR, Birtles RJ (2005) Molecular characterization of haemoparasites infecting bats (Microchiroptera) in Cornwall, UK. Parasitology131: 489–496. pmid:16174413
- View Article
- PubMed/NCBI
- Google Scholar
18. Kosoy M, Bai Y, Lynch T, Kuzmin I, Niezgoda M, et al. (2010) Bartonella spp. in bats, Kenya. Emerg Infect Dis 16: 1875–1881. pmid:21122216
- View Article
- PubMed/NCBI
- Google Scholar
19. Bai Y, Kosoy M, Recuenco S, Alvarez D, Moran D, et al (2011) Bartonella spp. in bats, Guatemala. Emerg Infects Dis 17: 1269–1271. pmid:21762584
- View Article
- PubMed/NCBI
- Google Scholar
20. Bai Y, Recuenco S, Turmelle A, Osikowicz LM, Gomez J, et al (2012) Prevalence and diversity of Bartonella spp. in bats in Peru. Am J Trop Med & Hyg 87: 518–523. pmid:25548383
- View Article
- PubMed/NCBI
- Google Scholar
21. Lin JW, Hsu YM, Chomel BB, Lin LK, Pei JC, et al (2012) Identification of novel Bartonella spp. in bats and evidence of Asian gray shrew as a new potential reservoir of Bartonella. Vet Microbiol 156: 119–126. pmid:22005177
- View Article
- PubMed/NCBI
- Google Scholar
22. Kamani J, Baneth G, Mitchell M, Mumcuoglu KY, Gutiérrez R, et al (2014) Bartonella species in bats (Chiroptera) and bat flies (Nycteribiidae) from Nigeria, West Africa. VBZD 14: 625–632.
- View Article
- Google Scholar
23. Veikkolainen V, Vesterinen EJ, Lilley TM, Pulliainen AT (2014) Bats as reservoir hosts of human bacterial pathogen, Bartonella mayotimonensis. Emerg Infect Dis 20: 960–967. pmid:24856523
- View Article
- PubMed/NCBI
- Google Scholar
24. Hayman DTS, McCrea R, Suu-Ire R, Wood JLN, Cunningham AA, et al (2012) Straw-colored fruit bat demography in Ghana. J Mammal 93: 1393–1404 pmid:23525358
- View Article
- PubMed/NCBI
- Google Scholar
25. Arvand M, Raoult D, Feil EJ (2010) Multi-locus sequence typing of a geographically and temporally diverse sample of the highly clonal human pathogen Bartonella quintana. PLoS One 5: e9765. pmid:20333257
- View Article
- PubMed/NCBI
- Google Scholar
26. Chaloner GL, Ventosilla P, Birtles RJ (2011) Multi-locus sequence analysis reveals profound genetic diversity among isolates of the human pathogen Bartonella bacilliformis. PLoS Negl Trop Dis 7: e1248. pmid:21811647
- View Article
- PubMed/NCBI
- Google Scholar
27. Mietze A, Morick D, Köhler H, Harrus S, Dehio C, et al (2011) Combined MLST and AFLP typing of Bartonella henselae isolated from cats reveals new sequence types and suggests clonal evolution. Vet Microbiol 148: 238–245. pmid:20863631
- View Article
- PubMed/NCBI
- Google Scholar
28. Bai Y, Malania L, Castillo DA, Moran D, Boonmar S, et al (2013) Global distribution of bartonella infections in domestic bovine and characterization of Bartonella bovis strains using multi-locus sequence typing. PloS One.
29. Paziewska A, Harris PD, Zwolińska L, Bajer A, Siński E (2011) Recombination within and between species of the alpha proteobacterium Bartonella infecting rodents. Microb Ecol 61: 134–145. pmid:20740281
- View Article
- PubMed/NCBI
- Google Scholar
30. Buffet JP, Pisanu B, Brisse S, Roussel S, Félix B, et al (2013) Deciphering Bartonella diversity, recombination, and host specificity in a rodent community. PLoS One 8: e68956. pmid:23894381
- View Article
- PubMed/NCBI
- Google Scholar
31. La Scola B, Zeaiter Z, Khamis A, Raoult D (2003) Gene-sequence-based criteria for species definition in bacteriology: the Bartonella paradigm. Trends Microbiol 11: 318–321. pmid:12875815
- View Article
- PubMed/NCBI
- Google Scholar
32. Vos M, Didelot X (2009) A comparison of homologous recombination rates in bacteria and archaea. ISME J 3: 199–208. pmid:18830278
- View Article
- PubMed/NCBI
- Google Scholar
33. Tavaré S (1986) Some probabilistic and statistical problems in the analysis of DNA sequences. Lectures Math Life Sci 17: 57–86.
- View Article
- Google Scholar
34. Bruen TC, Phillippe H, Bryant D (2006) A simple and robust statistical test for detecting the presence of recombination. Genetics 172: 2665–2681. pmid:16489234
- View Article
- PubMed/NCBI
- Google Scholar
35. Gundi VA, Billeter SA, Rood MP, Kosoy MY (2012) Bartonella spp. in rats and zoonoses, Los Angeles, California, USA. Emerg Infect Dis 18: 631–633. pmid:22469313
- View Article
- PubMed/NCBI
- Google Scholar
36. Zhu Q, Kosoy M, Olival KJ, Dittmar K (2014) Horizontal Transfers and Gene Losses in the Phospholipid Pathway of Bartonella Reveal Clues about Early Ecological Niches. Genome Biol Evol 6: 2156–2169. pmid:25106622
- View Article
- PubMed/NCBI
- Google Scholar
37. Koskiniemi S, Sun S, Berg OG, Andersson DI (2012) Selection-driven gene loss in bacteria. PLoS Genet 8: e1002787. pmid:22761588
- View Article
- PubMed/NCBI
- Google Scholar
38. Vitreschak AG, Rodionov DA, Mironov AA, Gelfand MS (2002) Regulation of riboflavin biosynthesis and transport genes in bacteria by transcriptional and translational attenuation. Nucleic Acids Res. 30: 3141–351. pmid:12136096
- View Article
- PubMed/NCBI
- Google Scholar
39. Berglund EC, Ellegaard K, Granberg F, Xie Z, Maruyama S, et al (2010) Rapid diversification by recombination in Bartonella grahamii from wild rodents in Asia contrasts with low levels of genomic divergence in Northern Europe and America. Mol Ecol 19: 2241–2255. pmid:20465583
- View Article
- PubMed/NCBI
- Google Scholar
40. Arvand M, Feil EJ, Giladi M, Boulouis HJ, Viezens J (2007) Multi-locus sequence typing of Bartonella henselae isolates from three continents reveals hypervirulent and feline-associated clones. PLoS One 2: e1346. pmid:18094753
- View Article
- PubMed/NCBI
- Google Scholar
41. Billeter SA, Hayman DTS, Peel A, Baker KS, Cunningham AA, et al (2012) Bartonella species in bat flies (Diptera: Nycteribiidae) from western Africa. Parasitology 139: 324–329. pmid:22309510
- View Article
- PubMed/NCBI
- Google Scholar
42. Kosoy M, Hayman DT, Chan KS (2012) Bartonella bacteria in nature: where does population variability end and a species start? Infect Genet Evol 12: 894–904. pmid:22449771
- View Article
- PubMed/NCBI
- Google Scholar
43. Bai Y, Kosoy MY, Cully JF, Bala T, Ray C, et al (2007) Acquisition of nonspecific Bartonella strains by the northern grasshopper mouse (Onychomys leucogaster). FEMS Microbiol Ecol 61: 438–48. pmid:17672850
- View Article
- PubMed/NCBI
- Google Scholar
44. Kosoy M, Mandel E, Green D, Marston E, Jones D, et al (2004) Prospective studies of Bartonella of rodents. Part II. Diverse infections in a single rodent community. Vector Borne Zoonotic Dis 4: 296–305. pmid:15671736
- View Article
- PubMed/NCBI
- Google Scholar
45. Chan KS, Kosoy M (2010) Analysis of multi-strain Bartonella pathogens in natural host population—do they behave as species or minor genetic variants? Epidemics 2: 165–172. pmid:21352787
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Welch DF, Pickett DA, Slater LN, Steigerwalt AG, Brenner DJ (1992) Rochalimaea henselae sp. nov., a cause of septicemia, bacillary angiomatosis, and parenchymal bacillary peliosis. J Clin Microbiol 30: 275–280. pmid:1537892
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Daly JS, Worthington MG, Brenner DJ, Moss CW, Hollis DG, et al (1993) Rochalimaea elizabethae sp. nov. isolated from a patient with endocarditis. J Clin Microbiol 31: 872–881. pmid:7681847
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Bass JW, Vincent JM, Person DA (1997) The expanding spectrum of bartonella infections II. Cat scratch disease. Pediatr Infect Dis J 16: 163–179. pmid:9041596
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Kerkhoff FT, Bergmans AM, van Der Zee A, Rothova A (1999) Demonstration of Bartonella grahamii DNA in ocular fluids of patient with neuroretinitis. J Clin Microbiol 37: 4034–4038. pmid:10565926
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Roux V, Eykyn SJ, Wyllie S, Raoult D (2000) Bartonella vinsonii subsp. berkhoffii as an agent of a febrile blood culture-negative endocarditis in a human. J Clin Microbiol 38: 1698–1700. pmid:10747175
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Kosoy MY, Murray M, Gilmore RD, Bai Y, Gage K (2003) Bartonella strains from ground squirrels are identical to Bartonella washoensis isolated from a human patient. J Clin Microbiol 41: 645–650. pmid:12574261
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Eremeeva ME, Gerns HL, Lydy SL, Goo JS, Ryan ET, et al (2007) Bacteremia, fever, and splenomegaly caused by a newly recognized Bartonella species. N Engl J Med 256: 2381–2387.
View Article
Google Scholar

[26] View Article

[27] Google Scholar

[ref8] 8. Kosoy M, Morway C, Sheff K, Bai Y, Colborn J, et al (2008) Bartonella tamiae sp. nov., a newly recognized pathogen isolated from human patients from Thailand. J Clin Microbiol 46: 772–775. pmid:18077632
View Article
PubMed/NCBI
Google Scholar

[29] View Article

[30] PubMed/NCBI

[31] Google Scholar

[ref9] 9. Welch DF, Carroll KC, Hofmeister EK, Persing DH, Robison DA, et al (1999) 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 37: 2598–2601. pmid:10405408
View Article
PubMed/NCBI
Google Scholar

[33] View Article

[34] PubMed/NCBI

[35] Google Scholar

[ref10] 10. Garcia-Caceres U, Garcia FU (1991) Bartonellosis: an immunodepressive disease and the life of Daniel Alcides Carrion. Am J Clin Pathol 95: s58–s66. pmid:2008885
View Article
PubMed/NCBI
Google Scholar

[37] View Article

[38] PubMed/NCBI

[39] Google Scholar

[ref11] 11. Chomel BB, Kasten RW, Floyd-Hawkins KA, Chi B, Yamamoto K, et al (1996). Experimental transmission of Bartonella henselae by the cat flea. J Clin Microbiol 34: 1952–1956. pmid:8818889
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref12] 12. Higgins JA, Radulovic S, Jaworski DC, Azad AF (1996) Acquisition of the cat scratch disease agent Bartonella henselae by cat fleas (Siphonaptera: Pulicidae). J Med Entomol 33: 490–495. pmid:8667399
View Article
PubMed/NCBI
Google Scholar

[45] View Article

[46] PubMed/NCBI

[47] Google Scholar

[ref13] 13. Pappalardo BL, Correa MT, York CC, Peat CY, Breitschwerdt EB (1997) Epidemiologic evaluation of the risk factors associated with exposure and seroreactivity to Bartonella vinsonii in dogs. Am J Vet Res 58: 467–471. pmid:9140552
View Article
PubMed/NCBI
Google Scholar

[49] View Article

[50] PubMed/NCBI

[51] Google Scholar

[ref14] 14. Roux V, Raoult D (1999) Body lice as tools for diagnosis and surveillance of reemerging diseases. J Clin Microbiol 37: 596–599. pmid:9986818
View Article
PubMed/NCBI
Google Scholar

[53] View Article

[54] PubMed/NCBI

[55] Google Scholar

[ref15] 15. Billeter SA, Sangmaneedet S, Kosakewich RC, Kosoy MY (2012) Bartonella species in dogs and their ectoparasites from Khon Kaen Province, Thailand. Southeast Asian J Trop Med Public Health 43: 1186–1192. pmid:23431825
View Article
PubMed/NCBI
Google Scholar

[57] View Article

[58] PubMed/NCBI

[59] Google Scholar

[ref16] 16. Luis AD, Hayman DTS, O’Shea TJ, Cryan PM, Gilbert AT, et al (2013) A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special? Proc Biol Sci 280: 20122753. pmid:23378666
View Article
PubMed/NCBI
Google Scholar

[61] View Article

[62] PubMed/NCBI

[63] Google Scholar

[ref17] 17. Concannon R, Wynn-Owen K, Simpson VR, Birtles RJ (2005) Molecular characterization of haemoparasites infecting bats (Microchiroptera) in Cornwall, UK. Parasitology131: 489–496. pmid:16174413
View Article
PubMed/NCBI
Google Scholar

[65] View Article

[66] PubMed/NCBI

[67] Google Scholar

[ref18] 18. Kosoy M, Bai Y, Lynch T, Kuzmin I, Niezgoda M, et al. (2010) Bartonella spp. in bats, Kenya. Emerg Infect Dis 16: 1875–1881. pmid:21122216
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref19] 19. Bai Y, Kosoy M, Recuenco S, Alvarez D, Moran D, et al (2011) Bartonella spp. in bats, Guatemala. Emerg Infects Dis 17: 1269–1271. pmid:21762584
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref20] 20. Bai Y, Recuenco S, Turmelle A, Osikowicz LM, Gomez J, et al (2012) Prevalence and diversity of Bartonella spp. in bats in Peru. Am J Trop Med & Hyg 87: 518–523. pmid:25548383
View Article
PubMed/NCBI
Google Scholar

[77] View Article

[78] PubMed/NCBI

[79] Google Scholar

[ref21] 21. Lin JW, Hsu YM, Chomel BB, Lin LK, Pei JC, et al (2012) Identification of novel Bartonella spp. in bats and evidence of Asian gray shrew as a new potential reservoir of Bartonella. Vet Microbiol 156: 119–126. pmid:22005177
View Article
PubMed/NCBI
Google Scholar

[81] View Article

[82] PubMed/NCBI

[83] Google Scholar

[ref22] 22. Kamani J, Baneth G, Mitchell M, Mumcuoglu KY, Gutiérrez R, et al (2014) Bartonella species in bats (Chiroptera) and bat flies (Nycteribiidae) from Nigeria, West Africa. VBZD 14: 625–632.
View Article
Google Scholar

[85] View Article

[86] Google Scholar

[ref23] 23. Veikkolainen V, Vesterinen EJ, Lilley TM, Pulliainen AT (2014) Bats as reservoir hosts of human bacterial pathogen, Bartonella mayotimonensis. Emerg Infect Dis 20: 960–967. pmid:24856523
View Article
PubMed/NCBI
Google Scholar

[88] View Article

[89] PubMed/NCBI

[90] Google Scholar

[ref24] 24. Hayman DTS, McCrea R, Suu-Ire R, Wood JLN, Cunningham AA, et al (2012) Straw-colored fruit bat demography in Ghana. J Mammal 93: 1393–1404 pmid:23525358
View Article
PubMed/NCBI
Google Scholar

[92] View Article

[93] PubMed/NCBI

[94] Google Scholar

[ref25] 25. Arvand M, Raoult D, Feil EJ (2010) Multi-locus sequence typing of a geographically and temporally diverse sample of the highly clonal human pathogen Bartonella quintana. PLoS One 5: e9765. pmid:20333257
View Article
PubMed/NCBI
Google Scholar

[96] View Article

[97] PubMed/NCBI

[98] Google Scholar

[ref26] 26. Chaloner GL, Ventosilla P, Birtles RJ (2011) Multi-locus sequence analysis reveals profound genetic diversity among isolates of the human pathogen Bartonella bacilliformis. PLoS Negl Trop Dis 7: e1248. pmid:21811647
View Article
PubMed/NCBI
Google Scholar

[100] View Article

[101] PubMed/NCBI

[102] Google Scholar

[ref27] 27. Mietze A, Morick D, Köhler H, Harrus S, Dehio C, et al (2011) Combined MLST and AFLP typing of Bartonella henselae isolated from cats reveals new sequence types and suggests clonal evolution. Vet Microbiol 148: 238–245. pmid:20863631
View Article
PubMed/NCBI
Google Scholar

[104] View Article

[105] PubMed/NCBI

[106] Google Scholar

[ref28] 28. Bai Y, Malania L, Castillo DA, Moran D, Boonmar S, et al (2013) Global distribution of bartonella infections in domestic bovine and characterization of Bartonella bovis strains using multi-locus sequence typing. PloS One.

[ref29] 29. Paziewska A, Harris PD, Zwolińska L, Bajer A, Siński E (2011) Recombination within and between species of the alpha proteobacterium Bartonella infecting rodents. Microb Ecol 61: 134–145. pmid:20740281
View Article
PubMed/NCBI
Google Scholar

[109] View Article

[110] PubMed/NCBI

[111] Google Scholar

[ref30] 30. Buffet JP, Pisanu B, Brisse S, Roussel S, Félix B, et al (2013) Deciphering Bartonella diversity, recombination, and host specificity in a rodent community. PLoS One 8: e68956. pmid:23894381
View Article
PubMed/NCBI
Google Scholar

[113] View Article

[114] PubMed/NCBI

[115] Google Scholar

[ref31] 31. La Scola B, Zeaiter Z, Khamis A, Raoult D (2003) Gene-sequence-based criteria for species definition in bacteriology: the Bartonella paradigm. Trends Microbiol 11: 318–321. pmid:12875815
View Article
PubMed/NCBI
Google Scholar

[117] View Article

[118] PubMed/NCBI

[119] Google Scholar

[ref32] 32. Vos M, Didelot X (2009) A comparison of homologous recombination rates in bacteria and archaea. ISME J 3: 199–208. pmid:18830278
View Article
PubMed/NCBI
Google Scholar

[121] View Article

[122] PubMed/NCBI

[123] Google Scholar

[ref33] 33. Tavaré S (1986) Some probabilistic and statistical problems in the analysis of DNA sequences. Lectures Math Life Sci 17: 57–86.
View Article
Google Scholar

[125] View Article

[126] Google Scholar

[ref34] 34. Bruen TC, Phillippe H, Bryant D (2006) A simple and robust statistical test for detecting the presence of recombination. Genetics 172: 2665–2681. pmid:16489234
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref35] 35. Gundi VA, Billeter SA, Rood MP, Kosoy MY (2012) Bartonella spp. in rats and zoonoses, Los Angeles, California, USA. Emerg Infect Dis 18: 631–633. pmid:22469313
View Article
PubMed/NCBI
Google Scholar

[132] View Article

[133] PubMed/NCBI

[134] Google Scholar

[ref36] 36. Zhu Q, Kosoy M, Olival KJ, Dittmar K (2014) Horizontal Transfers and Gene Losses in the Phospholipid Pathway of Bartonella Reveal Clues about Early Ecological Niches. Genome Biol Evol 6: 2156–2169. pmid:25106622
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref37] 37. Koskiniemi S, Sun S, Berg OG, Andersson DI (2012) Selection-driven gene loss in bacteria. PLoS Genet 8: e1002787. pmid:22761588
View Article
PubMed/NCBI
Google Scholar

[140] View Article

[141] PubMed/NCBI

[142] Google Scholar

[ref38] 38. Vitreschak AG, Rodionov DA, Mironov AA, Gelfand MS (2002) Regulation of riboflavin biosynthesis and transport genes in bacteria by transcriptional and translational attenuation. Nucleic Acids Res. 30: 3141–351. pmid:12136096
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

[ref39] 39. Berglund EC, Ellegaard K, Granberg F, Xie Z, Maruyama S, et al (2010) Rapid diversification by recombination in Bartonella grahamii from wild rodents in Asia contrasts with low levels of genomic divergence in Northern Europe and America. Mol Ecol 19: 2241–2255. pmid:20465583
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref40] 40. Arvand M, Feil EJ, Giladi M, Boulouis HJ, Viezens J (2007) Multi-locus sequence typing of Bartonella henselae isolates from three continents reveals hypervirulent and feline-associated clones. PLoS One 2: e1346. pmid:18094753
View Article
PubMed/NCBI
Google Scholar

[152] View Article

[153] PubMed/NCBI

[154] Google Scholar

[ref41] 41. Billeter SA, Hayman DTS, Peel A, Baker KS, Cunningham AA, et al (2012) Bartonella species in bat flies (Diptera: Nycteribiidae) from western Africa. Parasitology 139: 324–329. pmid:22309510
View Article
PubMed/NCBI
Google Scholar

[156] View Article

[157] PubMed/NCBI

[158] Google Scholar

[ref42] 42. Kosoy M, Hayman DT, Chan KS (2012) Bartonella bacteria in nature: where does population variability end and a species start? Infect Genet Evol 12: 894–904. pmid:22449771
View Article
PubMed/NCBI
Google Scholar

[160] View Article

[161] PubMed/NCBI

[162] Google Scholar

[ref43] 43. Bai Y, Kosoy MY, Cully JF, Bala T, Ray C, et al (2007) Acquisition of nonspecific Bartonella strains by the northern grasshopper mouse (Onychomys leucogaster). FEMS Microbiol Ecol 61: 438–48. pmid:17672850
View Article
PubMed/NCBI
Google Scholar

[164] View Article

[165] PubMed/NCBI

[166] Google Scholar

[ref44] 44. Kosoy M, Mandel E, Green D, Marston E, Jones D, et al (2004) Prospective studies of Bartonella of rodents. Part II. Diverse infections in a single rodent community. Vector Borne Zoonotic Dis 4: 296–305. pmid:15671736
View Article
PubMed/NCBI
Google Scholar

[168] View Article

[169] PubMed/NCBI

[170] Google Scholar

[ref45] 45. Chan KS, Kosoy M (2010) Analysis of multi-strain Bartonella pathogens in natural host population—do they behave as species or minor genetic variants? Epidemics 2: 165–172. pmid:21352787
View Article
PubMed/NCBI
Google Scholar

[172] View Article

[173] PubMed/NCBI

[174] Google Scholar

Figures

Abstract

Author Summary

Introduction

Materials and Methods

Bartonella isolates selection

Multi-locus sequence typing (MLST)

Analysis of nucleotide polymorphism and diversity

Phylogenetic analysis

Recombination tests

Nucleotide sequence accession numbers

Results

Individual sequence analysis of gltA and other loci

Allelic profiles, sequence types (ST), and phylogenetic analysis

Patterns of selection and diversity in nucleotide sequences

Recombination test

Discussion

Supporting Information

S1 Table. GenBank accession numbers for ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA sequences for bat-associated Bartonella strains.

S2 Table. GenBank accession numbers for ftsZ, gltA, nuoG, ribC, rpoB, ssrA, ITS, and 16S rRNA sequences of reference Bartonella species.

Acknowledgments

Author Contributions

References