PLoS OnePLoS ONEplosplosonePLoS ONE1932-6203Public Library of ScienceSan Francisco, USA242783423836770PONE-D-13-3434810.1371/journal.pone.0080894Research ArticleGlobal Distribution of Bartonella Infections in Domestic Bovine and Characterization of Bartonella bovis Strains Using Multi-Locus Sequence TypingBartonella in Bovids and MLST of B. bovisBaiYing1*MalaniaLile2Alvarez CastilloDanilo3MoranDavid3BoonmarSumalee4ChanlunAran5SuksawatFanan5MaruyamaSoichi6KnobelDarryn78KosoyMichael1Bacterial Disease Branch, Division of Vector-Borne Disease, Centers for Disease Control and Prevention, Fort Collins, Colorado, United States of AmericaNational Center for Disease Control and Public Health, Tbilisi, Republic of GeorgiaCentro de Estudios en Salud, Universidad del Valle de Guatemala, Guatemala, GuatemalaInternational Emerging Infectious Program, Thailand MOPH-US CDC Collaboration, Nonthaburi, ThailandKhon Kaen University, Khon Kaen, ThailandDepartment of Veterinary Medicine, College of Bioresource Science, Nihon University, Kanagawa, JapanDepartment of Veterinary Tropical Diseases, University of Pretoria, Pretoria, South AfricaKenya Medical Research Institute/Centers for Disease Control and Prevention Public Health Collaboration, Kisumu, KenyaZilbersteinDanEditorTechnion-Israel Institute of Technology Haifa 32000 Israel, Israel* E-mail: bby5@cdc.gov

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

Conceived and designed the experiments: YB MK. Performed the experiments: YB LM SB SM. Analyzed the data: YB. Contributed reagents/materials/analysis tools: LM SB FS DAC DM SM DK AC. Wrote the manuscript: YB MK DK.

201321112013811e80894882013171020132013This 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.

Bartonella bovis is commonly detected in cattle. One B. bovis strain was recently isolated from a cow with endocarditis in the USA, suggesting its role as an animal pathogen. In the present study, we investigated bartonella infections in 893 cattle from five countries (Kenya, Thailand, Japan, Georgia, and Guatemala) and 103 water buffaloes from Thailand to compare the prevalence of the infection among different regions and different bovid hosts. We developed a multi-locus sequence typing (MLST) scheme based on nine loci (16S rRNA, gltA, ftsZ, groEL, nuoG, ribC, rpoB, ssrA, and ITS) to compare genetic divergence of B. bovis strains, including 26 representatives from the present study and two previously described reference strains (one from French cows and another from a cow with endocarditis in the USA). Bartonella bacteria were cultured in 6.8% (7/103) of water buffaloes from Thailand; all were B. bovis. The prevalence of bartonella infections in cattle varied tremendously across the investigated regions. In Japan, Kenya, and the Mestia district of Georgia, cattle were free from the infection; in Thailand, Guatemala, and the Dusheti and Marneuli districts of Georgia, cattle were infected with prevalences of 10–90%. The Bartonella isolates from cattle belonged to three species: B. bovis (n=165), B. chomelii (n=9), and B. schoenbuchensis (n=1), with the latter two species found in Georgia only. MLST analysis suggested genetic variations among the 28 analyzed B. bovis strains, which fall into 3 lineages (I, II, and III). Lineages I and II were found in cattle while lineage III was restricted to water buffaloes. The majority of strains (17/28), together with the strain causing endocarditis in a cow in the USA, belonged to lineage I. Further investigations are needed to determine whether B. bovis causes disease in bovids.

Specimens from cattle in Kenya were collected under a project supported by the Wellcome Trust, UK (grant number 081828/B/06/Z). The funders had no role to publish, or preparation of the manuscript.
Introduction

Found throughout much of the world, cattle are the most common type of large domesticated ungulate. They comprise hundreds of breeds that are recognized worldwide. Previously described as two separate subspecies (Bos taurus subsp. taurus and B. taurus subsp. indicus), all cattle belong to the sub-tribe Bovina under tribe Bovini according to a very recent classification based on multiple autosomal introns [1]. Water buffaloes (Bubalus bubalis) resemble cattle in many characteristics but belong to Bubalina, another sub-tribe under Bovini [1]. These animals are agriculturally and economically important to humans as they are widely used for a variety of purposes, such as meat, dairy products, and draught power. However, a variety of infections in cattle, such as bovine spongiform encephalopathy, foot-and-mouth diseases, brucellosis, anthrax, infectious endocarditis, and many others, have caused either economical loss of these animals or big public health concerns as some of these infections are transmissible to humans and potentially can cause death.

The genus Bartonella contains B. bovis and many other species and/or subspecies. These fastidious Gram-negative bacteria infect and persist in mammalian erythrocytes and endothelial cells and are found in a wide range of wild and domesticated mammals, including rodents, insectivores, carnivores, ungulates, and others. A number of Bartonella species have been associated with human illnesses and are responsible for a growing spectrum of emerging diseases, including endocarditis [2-10]. 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 [11-16].

Infection with Bartonella spp. (mainly B. bovis) in cattle has been reported in a number of studies [17-20], with B. bovis originally described in cows from France [21]. For water buffaloes, there are no previous reports of infections with Bartonella spp. The prevalence of B. bovis in cattle is generally high but varies widely across studies from different countries, i.e., 50-89%, 70%, 36%, 24%, and 20%, in the USA, French Guyanna, France, Italy, and West Africa, respectively [18-23]. In a very recently study from Poland, the prevalence of B. bovis in cattle was much lower (6.8%) [24]. In addition to B. bovis, two more Bartonella spp. (B. chomelii and B. schenobuensis) are occasionally found in cattle [17,25]. Although Bartonella spp. do not usually cause clinical signs in cattle, B. bovis is associated with bovine endocarditis [26,27]. Understanding whether any specific genetic variations are responsible for the pathogenic outcome in animals and people is important and requires separate investigation.

In the present study, we first investigated the prevalence of Bartonella spp. in cattle from five countries across the world (Kenya, Thailand, Japan, Georgia, and Guatemala) and in water buffaloes from Thailand. Following that, we characterized the B. bovis strains obtained from the present study and from reference collections (France and USA), and then compared the genetic difference among these strains. We employed multi-locus sequence typing (MLST) that is based on comparison of nucleotide sequences derived from multiple loci. This approach has been applied to study genetic diversity of different agents, including Bartonella spp. [28-30]. It has been shown to provide high discriminatory power in epidemiological and genetic analysis of strain populations while retaining signatures of longer-term evolutionary relationships or clonal stability. This in turn can help to enhance our understanding of population structure of the bacteria and the relationships between sequence type and animal host. Our aims were: (1) to assess the apparent prevalence of bartonella culture-positive infections in cattle and water buffaloes from different countries; (2) to identify Bartonella spp. in positive samples using the citrate synthase gene (gltA); and (3) to develop a MLST scheme for B. bovis, and to validate this scheme against a diverse sample of strains representing each country in this study to allow comparison of their sequence diversity.

Materials and MethodsEthics Statement

Specimen collection from cattle in Kenya was approved by the Institutional Animal Care and Use Committees (IACUC) of the Kenya Medical Research Institute (protocol no. 1191) and CDC (1562BRETBDX). Sampling cattle does not require IACUC approval in Guatemala. There was no IACUC or ethics committee in Georgia, Thailand, and Japan at the time of the conducted field study.

Cattle samples were collected in private ranches in all sites. There was no specific permission required. Verbal or written informed consent was obtained from all animal owners before specimen collection. Sample collection was done by veterinarians, experienced field technicians, and trained animal health assistants. All animal procedures were conducted in compliance with the Animal Welfare Act and other Federal statutes and regulations relating to animals and experiments involving animals and adheres to principles stated in the Guide for the Care and Use of Laboratory Animals (National Research Council [NRC] Publication. Samples from Guatemala and Thailand were obtained specifically for this study; and samples from Kenya, Japan, and Republic of Georgia were obtained as part of routine care. Three milliliters of blood were collected from each animal.

Study sites, sample collection, and reference strains

Whole blood samples were collected from healthy cattle and water buffaloes from five countries during different time periods (Table 1). Samples collected included (1) 40 cattle from Chachoeng Sao Province of Thailand and 103 water buffaloes from Chachoeng Sao, Nakhon Pathom, and Khon Kaen provinces of Thailand; (2) 305 cattle from Chiba, Hokkaido, Kagoshima, Kanagawa, and Okinawa prefectures of Japan; (3) 159 cattle from Dusheti, Marneuli, and Mestia districts of Republic of Georgia; (4) 389 cattle from Alta Verapaz, Huehuetenango, Izabal, Quetzaltenango, Quiche, and Retalhuleu departments of Guatemala; and (5) 221 cattle from Nyanza Province of Kenya. Blood samples were kept at -20° C or lower temperature until tested. All samples were sent to the Bartonella Laboratory, CDC in Fort Collins, Colorado, USA for bartonella testing, except samples from Japan, which were tested in Laboratory of Veterinary Public Health, Nihon University, Kanagawa, Japan using analogous culturing techniques [17].

10.1371/journal.pone.0080894.t001Distribution of bartonella infection in cattle from different geographic regions.
CountryRegion# Tested# PositivePrevalence (%)Collection date
ThailandChachoeng Sao Province40410September, 2008
JapanChiba Prefecture3700July, 2000 - February, 2005
Hokkaido Prefecture5200
Kagoshima Prefecture13900
Kanagawa Prefecture5700
Okinawa Prefecture2000
total30500
GeorgiaDusheti District1007373May - June, 2010
Marneuli District201890
Mestia District3900
total1599157.2
GuatemalaAlta Verapaz Department601525April - December, 2011
Huehuetenango Department1101210.9
Izabal Department60610
Quetzaltenango Department371848.6
Quiche Department622235.5
Retalhuleu Department60711.7
total3898020.6
KenyaAsembo22100January - March, 2009

Two reference strains (B18962, and B37080) that have been identified as B. bovis were included in the MLST analysis. B18962 is the subculture of strain 91-4T (the type strain of B. bovis) described in a cow in France [21], and B37080 is an isolate from a cow with endocarditis in the USA [27].

Bartonella culturing

Animal bloods were thawed at 4° C and re-suspended 1:4 in brain heart infusion broth supplemented with 5% amphotericin B (to reduce likelihood that bacterial and fungal contaminants would overgrow the overgrow the slow-growing bartonella bacteria), then plated on heart infusion agar containing 10% rabbit blood and incubated in an aerobic atmosphere with 5% carbon dioxide at 35°C up to four weeks. Bacterial growth was monitored at the end of each week. Bacterial colonies were presumptively identified as bartonella based on colony morphology. Subcultures of bartonella colonies from the original agar plate were streaked onto secondary agar plates and incubated at the same conditions until sufficient growth was observed, usually between 5-7 days. Pure cultures were harvested in 10% glycerol.

DNA preparation, PCR verification and species differentiation by <italic>gltA</italic>

We extracted genomic DNA from pure culture of each isolate using the QIAxtractor automated DNA purification system (Qiagen, Valencia, CA). We performed PCR using primers 430F and 1210R to amplify a specific region in the citrate synthase gene (gltA) of the genus of Bartonella for verification of Bartonella strains. Positive PCR products were sequenced using the same primers as the initial PCR assay at a concentration of 1.6 µM. Using Lasergene software package (DNASTAR, Madison, WI), the gltA sequences obtained from all samples in the present study were compared with B. bovis, other ruminant-associated Bartonella, and other known Bartonella spp. for species identification.

MLST analysis

Strains identified as B. bovis based on gltA sequences were further analyzed using MLST for eight additional loci (16S rRNA, ftsZ, groEL, nuoG, ribC, rpoB, ssrA, and ITS). Primers and other related information are provided in Table 2. We randomly selected one to three strains from each region of a country where we found bartonella to be present in the animals. Including the two reference strains from France and USA, a total of 28 strains were characterized for the MLST analysis (Table 3). A neighbor-joining tree based on the concatenated MLST alleles was constructed using the Clustal W program within MegAlign of the Lasergene package.

10.1371/journal.pone.0080894.t002Characteristics of the nine loci evaluated for the <italic>B. bovis</italic> MLST scheme.
LocusForward primer Reverse primer length of analyzed sequence (bp)No. variable sitesNo. alleles
16SrRNAcaggcctaacacatgcaagtcgggcggwgtgtacaaggc117201
ftsZattaatctgcaycggccagaacvgadacacgaataacacc885207
gltAgctatgtctgcattctatcagat cyt caa tca ttt ctt tcc a753185
groELgaactngaagataagttngaaaatccattccgcccattc10812712
nuoGggcgtgattgttctcgttacacgaccacggctatcaat328184
ribCtaaccgatattggttgtgttgaagtaaagctagaaagtctggcaacataacg535188
rpoBcgcattggcttacttcgtatg gtagactgattagaacgctg852128
ssrAgctatggtaataaatggacaatgaaataagcttctgttgccaggtg28745
ITScttcagatgatgatcccaagccttctggcggaaccgacgaccccctgcttgcaaag a364-398487
10.1371/journal.pone.0080894.t003Allelic profiles, sequence types (ST), and lineage group (LG) classification for the 29 <italic>B. bovis</italic> isoaltes from different geographic locations and bovine hosts.
#Isolate16S rRNAftsZgltAgroELnuoGribCrpoBssrAITSSTLGCountryHost
1B18962111111111ST11Francecattle
2B37080111211111ST21US cattle
3B33663111311112ST31Guatemalacattle
4B38041 111311112ST31Guatemalacattle
5B33664111211213ST41Guatemalacattle
6B38038111411313ST51Guatemalacattle
7B33695122521511ST62Guatemalacattle
8B38216122522514ST72Guatemalacattle
9B38241122522514ST72Guatemalacattle
10B38035131614611ST81Guatemalacattle
11B38215131715315ST91Guatemalacattle
12B38223141131616ST101Guatemalacattle
13B38240153631311ST111Guatemalacattle
14B31166141131611ST121Georgiacattle
15B31219141131611ST121Georgiacattle
16B31170141136611ST131Georgiacattle
17B31158141836611ST141Georgiacattle
18B31167141836611ST141Georgiacattle
19B31178151131231ST151Georgiacattle
20B31182151911611ST161Georgiacattle
21B25099122512424ST172Thailandcattle
22B25093122522424ST182Thailandcattle
23B25100122523424ST192Thailandcattle
24B326741641047747ST203Thailandwater buffalo
25B327801641047747ST203Thailandwater buffalo
26B328501641047747ST203Thailandwater buffalo
27B327811641147747ST213Thailandwater buffalo
28B327301751248857ST223Thailandwater buffalo

Newly-identified alleles from the current study were submitted to GenBank with the following accession numbers: KF193407 - KF193413, KF199895 - KF199899, KF212449 - KF212460, KF218206 - KF218209, KF218210 - KF218216, KF218217 - KF218224, KF218225 - KF218229, and KF218230 - KF218236 for ftsZ, gltA, groEL, nuoG, ribC, rpoB, ssrA, and ITS, respectively.

ResultsBartonella prevalence

Seven bartonella isolates were obtained from water buffaloes from Thailand, giving an apparent prevalence of 6.8% (7/103). In cattle, a total of 175 bartonella isolates were obtained from all countries. The prevalence of bartonella varied widely across the study sites. In Thailand, bartonella was isolated from 10% (4/40) of the cattle. In Guatemala, the overall bartonella prevalence was 20.6% (80/389) with a range of 10% - 48.7% across the six departments. In Georgia, bartonella prevalence was 73% (73/100) and 90% (18/20) in cattle from Dusheti and Marneuli districts, respectively, but bacteria were absent (0/39) in cattle from Mestia district. In western Kenya and all sampled sites of Japan, bartonella was absent in cattle, despite the large number of samples investigated in both sites. Details are provided in Table 1.

<italic>Bartonella</italic> spp. identification

All 182 isolates obtained from water buffaloes and cattle from the five countries in the present study were characterized by sequencing gltA. The gltA sequences showed that all seven isolates from water buffaloes and 165 isolates from cattle from Thailand, Guatemala and Georgia, were B. bovis with 99.1-100% identity to the type strain (AF293394); nine isolates recovered in cattle from Dusheti District of Georgia was B. chomelii with 98.9-100% identity to a previously described variant (AY254308) from cattle in France; and the last isolate recovered in cattle from Marneuli District of Georgia was B. schoenbuchensis with 97.2% similarity to the B. schoenbuchensis variant (AJ567635) from China.

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

The size of sequenced alleles ranged between 287bp - 1172bp at different loci (Table 2) with the total length of concatenated sequences 5085bp - 5119bp. Noticeably, the ITS locus varied in size among the strains, ranging between 364bp - 398bp. Specifically, sequences of most strains were 386bp long while sequences of all strains from Thai cattle and two strains from Guatemala (B38216 and B38241) were shorter (364bp) with a 21-bp deletion, and on contrary, sequences of all water buffalo-originated strains were longer (398bp) with a13-bp insertion (data not shown). Sequencing analysis demonstrated that 16S rRNA was invariant and was identical to the previously described variant (NR025121) from B. bovis type strain. All other loci showed large variation, with variable sites of 4 - 48, and 4 - 12 alleles by loci (Table 2).

Due to the identity of 16S rRNA for all strains, the sequences of this locus were excluded from the MLST analysis. Based on the concatenated sequences of other eight loci, the MLST analysis distinguished 22 STs among the 28 strains. ST20 was represented by three strains; and ST3, ST7, ST12, and ST14 each were represented by two strains; all other 17 STs each were represented only by one strain. Phylogenetic analysis demonstrated that all STs resolved into three lineages (Figure 1, 16S rRNA not included) with low divergence (0.1 - 1.8%) among all STs. Lineage I contains 14 STs of 17 strains with divergence among the STs of 0.1 - 0.5%. Strains in this group have diverse geographical distribution, which were found in cattle from all studied regions except for Thailand; lineage II contains five STs of six strains with divergence of 0.1 - 0.4% among the STs. The strains within lineage II were obtained from cattle from Thailand (3 strains, of ST6 and ST7) and Guatemala (3 strains, of STs17 - 19); and lineage III contains five strains obtained exclusively from water buffaloes from Thailand. The five strains belonged to three STs (STs20 - 22) with divergence of 0.1 - 0.4%. The distances between lineage III and either lineage I (1.4 - 1.5%) or lineage II (1.7 - 1.8%) were much more distant comparing to that between lineages I and II (0.7 - 1.1%).

10.1371/journal.pone.0080894.g001Phylogenetic relationships of the 28 <italic>B. bovis</italic> isolates inferred from 5085bp - 5119bp concatenated sequences of <italic>ftsZ</italic>, <italic>gltA</italic>, groEL, <italic>nuoG</italic>, <italic>ribC</italic>, <italic>rpoB</italic>, <italic>ssrA</italic>, and ITS fragments.

Following each isolate are geographical origin, bovine host source, and sequence type (ST) classification. Phylogenetic tree was constructed by N-J method, and bootstrap values were calculated with 1,000 replicates. A total of 22 STs were identified, which fall into three lineages (square circled clades).

Host association and geographical distribution of <italic>B. bovis</italic> sequence types

Strains from the same location exhibited different sequence types, showing little evidence of concordance between the MLST data and geographical origin. However, all sequence types showed clear host association, either with cattle or water buffaloes. Furthermore, except for Guatemalan strains, other STs obtained from the same area were more associated with a specific lineage. For example, all STs from Georgia, France, and the US belong to lineage I; STs obtained from Thai cattle belong to lineage II; while all STs within lineage III were exclusively from water buffaloes. For the 11 STs obtained from cattle from Guatemala, eight of them fall into lineage I while the other three fall into lineage II.

Discussion

In this study, we investigated infections with Bartonella spp. in water buffaloes from Thailand and cattle from five countries across the world. We provided the first evidence that water buffaloes are also infected with Bartonella spp., with a prevalence of 7%. In cattle, bartonella prevalence varied widely between the studied countries and between different regions in the same country. In Japan, Kenya, and the Mestia district of Georgia, the cattle were free from the infection; in contrast, infection prevalence reached as high as 90% in the Marneuli district of Georgia, and a moderate prevalence was observed in other areas. It is of interest to note that the prevalence of infection in cattle varied widely across districts in Georgia, from apparent absence to nearly 100%. The observed variation in bartonella prevalence in cattle between different countries/regions can be the result of multiple factors. Since bartonellae are mainly vector-transmitted [11-16], we speculate that the distribution and abundance of specific arthropods play a major role in this matter. It may be that a heavy level of infestation is required for transmission. Cattle ticks (Rhipicephalus microplus) and biting flies (Diptera spp.) are implicated as potential vectors that may transmit bartonella between cattle [31,32]. Interestingly, low ectoparasite infestation on cattle in Japan was noticed (Maruyama, personal communication). Such a fact could provide a plausible explanation for the absence of bartonella infection observed in cattle from Japan. In addition to ectoparasites, environmental factors, such as geographic characters, landscape, etc., may also have influences on bartonella prevalence. Of the three districts in Georgia, Mestia district (with no bartonella infection in cattle) is located in a mountainous area in the northwest of the country, while the other two districts (with high bartonella infection in cattle) are located in lower-lying areas in the northeast and southeast.

The cultures obtained from the study belonged to three Bartonella spp. based on gltA sequences. The majority of strains (172/182), identified as B. bovis, were the most common in cattle from all investigated places. This is concordant with previously reported observations. The other two species, B. chomelii and B. schoenbuchensis, are not common, although both have been previously reported in cattle [17,25]. Interestingly, both of these species were only found in cattle from Georgia. This may suggest that the bartonella community associated with cattle in Georgia is more diverse comparing to other places. All nine B. chomelii isolates were obtained from cattle in Dusheti District located in northeastern Georgia. These observations may be associated with the composition of the local ectoparasite community. Future studies of cattle ectoparasites should test hypotheses about whether any particular arthropod species act as vectors for bartonella transmission between cattle. Finally, all isolates obtained from water buffaloes in Thailand also belong to B. bovis, indicating that water buffaloes can also serve as hosts for B. bovis.

Although B. bovis is known to be widely distributed in bovines, no study has demonstrated how strains from different geographic areas and hosts vary. As a potential pathogen for domestic animals, information about genetic diversity of B. bovis is important. In this study, we developed an MLST scheme for B. bovis based on nine loci to characterize B. bovis strains among geographically diverse populations. Resolving into 22 STs among the 28 strains from different geographical regions, our MLST data demonstrated B. bovis strains can be genetically different. Although some STs can be quite close, the same ST was never found in strains from different regions.

Multi-locus sequencing typing has been applied in a number of studies previously, in which the MLST was based on the comparison of multiple housekeeping gene sequences [28-30,33]. In this study, we included ITS. Although non-functional, sequence comparison of the ITS region is widely used in taxonomy and molecular phylogeny for its ease of amplification and high degree of variation, even between closely-related species. Our data showed either insertion or deletion has occurred in some of the tested strains, which likely are associated with host species. Such results suggested that ITS might be more powerful in differentiating genetic diversities among strain populations, compared to some regularly-used housekeeping genes.

Identification of three close but distinct lineages among the STs suggests a clonal population structure for the species. All STs of cattle-originated strains fall into either lineage I or II, while STs of water buffalo-originated strains exclusively fall in lineage III, showing the specific host relationship. Meanwhile, we suspect that lineage I and lineage II each is associated with a particular lineage of Bos taurus. By general understanding, cattle from US, Europe, and Georgia belong to the ‘taurine’ lineage (former B. taurus subsp. taurus); while those from Thailand belong to the ‘zebu’ lineage (former B. taurus subsp. indicus); and in Guatemala, cattle were of mixed breeds. Based on this information, we hypothesize that lineage I is associated with cattle of ‘taurine’ lineage and lineage II is associated with cattle of ‘zebu’ lineage. Further studies are required to confirm this hypothesis.

MLST analysis also showed B. bovis strains retain geographical particularity. This observation can be explained by the association of B. bovis with specific breed of cattle as well. For example, all cattle-originated strains from Thailand belong to lineage II; while all strains from Georgian cattle together with strains from the US and France, belong to lineage I. Contrastingly, strains from Guatemalan cattle fall into both lineages I and II. Mixed breed of cattle in Guatemala likely is the cause of crossing lineages.

Conclusion

Lineage I was the most common group of B. bovis, containing most STs (13/22) that are distributed in different regions across the world. The isolate that caused cattle endocarditis in the US (B37080, ST2) was also within this lineage. Although none of the other tested strains is of the same ST, all STs in this lineage are in fact very close. While clinical symptoms were not recorded for cattle in these studies, B. bovis strains have high potential as a pathogen, being widely distributed in cattle populations globally.

This study was supported by US CDC Global Disease Detection program. Samples from Georgia were collected by the National Center for Disease Control and Public Health field group. Special thanks to Marine Ramishvili (NCDCPH, Georgia), Ramon Medrano and Maria Fernanda Rizzo (Centro de Estudios en Salud, Universidad del Valle de Guatemala, Guatemala, Guatemala) for helping collect cattle samples in Georgia and Guatemala, respectively.

ReferencesHassaninA, AnJ, RopiquetA, NguyenTT, CoulouxA (2013) Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: Application to the tribe Bovini (Cetartiodactyla, Bovidae). Mol Phylogenet Evol66: 766-775.10.1016/j.ympev.2012.11.003 PubMed: 2315989423159894BassJW, VincentJM, PersonDA (1997) The expanding spectrum of Bartonella infections II. Cat scratch disease. Pediatr Infect Dis J16: 163-179.10.1097/00006454-199702000-00002 PubMed: 90415969041596DalyJS, WorthingtonMG, BrennerDJ, MossCW, HollisDGet al. (1993) Rochalimaea elizabethae sp. nov. isolated from a patient with endocarditis. J Clin Microbiol31: 872-881 PubMed: 76818477681847EremeevaME, GernsHL, LydySL, GooJS, RyanETet al. (2007) Bacteremia, fever, and splenomegaly caused by a newly recognized Bartonella species. N Engl J Med256: 2381-2387 PubMed: 1755411917554119KerkhoffFT, BergmansAM, van der ZeeA, RothovaAet al. (1999) Demonstration of Bartonella grahamii DNA in ocular fluids of patient with neuroretinitis. J Clin Microbiol37: 4034-4038 PubMed: 1056592610565926KosoyMY, MurrayM, GilmoreRD, BaiY, GageKL (2003) Bartonella strains from ground squirrels are identical to Bartonella washoensis isolated from a human patient. J Clin Microbiol41: 645-650.10.1128/JCM.41.2.645-650.2003 PubMed: 1257426112574261KosoyM, MorwayC, SheffKW, BaiY, ColbornJet al. (2008) Bartonella tamiae sp. nov., a newly recognized pathogen isolated from human patients from Thailand. J Clin Microbiol46: 772-775.10.1128/JCM.02120-07 PubMed: 1807763218077632RouxV, EykynSJ, WyllieS, RaoultD (2000) Bartonella vinsonii subsp. berkhoffii as an agent of a febrile blood culture-negative endocarditis in a human. J Clin Microbiol38: 1698-1700 PubMed: 1074717510747175WelchDF, PickettDA, SlaterLN, SteigerwaltAG, BrennerDJ (1992) Rochalimaea henselae sp. nov., a cause of septicemia, bacillary angiomatosis, and parenchymal bacillary peliosis. J Clin Microbiol30: 275-280 PubMed: 15378921537892WelchDF, CarrollKC, HofmeisterEK, PersingDH, RobisonDAet 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 Microbiol37: 2598-2601 PubMed: 1040540810405408ChomelBB, KastenRW, Floyd-HawkinsKA, ChiB, YamamotoKet al. (1996) Experimental transmission of Bartonella henselae by the cat flea. J Clin Microbiol34: 1952-1956 PubMed: 88188898818889Garcia-CaceresU, GarciaFU (1991) Bartonellosis: an immunodepressive disease and the life of Daniel Alcides Carrion. Am J Clin Pathol95: s58-s66 PubMed: 2008885 2008885HigginsJA, RadulovicS, JaworskiDC, AzadAF (1996) Acquisition of the cat scratch disease agent Bartonella henselae by cat fleas (Siphonaptera: Pulicidae). J Med Entomol33: 490-495 PubMed: 86673998667399PappalardoBL, CorreaMT, YorkCC, PeatCY, BreitschwerdtEB (1997) Epidemiologic evaluation of the risk factors associated with exposure and seroreactivity to Bartonella vinsonii in dogs. Am J Vet Res58: 467-471 PubMed: 91405529140552RouxV, RaoultD (1999) Body lice as tools for diagnosis and surveillance of reemerging diseases. J Clin Microbiol37: 596-599 PubMed: 99868189986818BilleterSA, SangmaneedetS, KosakewichRC, KosoyMY (2012) Bartonella species in dogs and their ectoparasites from Khon Kaen Province, Thailand. Southeast Asian J Trop Med Public Health. 43: 1186-1192 PubMed: 2343182523431825ChangCC, ChomelBB, KastenRW, HellerRM, UenoHet al. (2000) Bartonella spp. isolated from wild and domestic ruminants in North America. Emerg Infect Dis6: 306-311.10.3201/eid0603.000313 PubMed: 1082712310827123MartiniM, MenandroML, MondinA, PasottoD, MazzariolSet al. (2008) Detection of Bartonella bovis in a cattle herd in Italy. Vet Rec162: 58-59.10.1136/vr.162.2.58 PubMed: 1819266018192660RaoultD, La ScolaB, KellyPJ, DavoustB, GomezJ (2005) Bartonella bovis in cattle in Africa. Vet Microbiol105: 155-156.10.1016/j.vetmic.2004.10.013 PubMed: 1562752815627528SaisongkorhW, BarrassiL, DavoustB, de BrouckerCA, RaoultDet al. (2009) First isolation of Bartonella bovis from animals in French Guyana, South America. Clin Microbiol Infect15Suppl: 124-126.10.1111/j.1469-0691.2008.02198.x PubMed: 19456818 19456818BermondD, BoulouisHJ, HellerR, Van LaereG, MonteilHet al. (2002) Bartonella bovis sp. nov. and Bartonella capreoli sp. nov., isolated from European ruminants. Int J Syst Evol Microbiol52: 383-390 PubMed: 11931146 11931146BreitschwerdtEB, SontakkeS, CannedyA, HancockSI, BradleyJM (2001) Infection with Bartonella weissii and detection of Nanobacterium antigens in a North Carolina beef herd. J Clin Microbiol39: 879-882.10.1128/JCM.39.3.879-882.2001 PubMed: 1123039811230398CherryNA, MaggiRG, CannedyAL, BreitschwerdtEB (2009) PCR detection of Bartonella bovis and Bartonella henselae in the blood of beef cattle. Vet Microbiol135: 308-312.10.1016/j.vetmic.2008.09.063 PubMed: 1901957419019574Welc-FalęciakR, GronoK (2013) The first cases of Bartonella bovis infection in cattle from Central Europe. Vet Microbiol162: 954-956.10.1016/j.vetmic.2012.09.019 PubMed: 23063105 23063105RolainJM, RoussetE, La ScolaB, DuquesnelR, RaoultD (2003) Bartonella schoenbuchensis isolated from the blood of a French cow. Ann N Y Acad Sci990: 236-238.10.1111/j.1749-6632.2003.tb07370.x PubMed: 1286063312860633MaillardR, PetitE, ChomelB, LacrouxC, Schelcher et al (2007) Endocarditis in cattle caused by Bartonella bovis. Emerg Infect Dis13: 1383-1385.10.3201/eid1309.070236 PubMed: 1825211618252116ErolE, JacksonC, BaiY, SellsS, LockeSet al. (2013) Bartonella bovis isolated from a cow with endocarditis. J Vet Diagn Invest25: 288-290.10.1177/1040638713477408 PubMed: 23512923 23512923ArvandM, RaoultD, FeilEJ (2010) Multi-locus sequence typing of a geographically and temporally diverse sample of the highly clonal human pathogen Bartonella quintana. PLOS ONE. 5: e9765.10.1371/journal.pone.0009765 PubMed: 2033325720333257ChalonerGL, VentosillaP, BirtlesRJ (2011) Multi-locus sequence analysis reveals profound genetic diversity among isolates of the human pathogen Bartonella bacilliformis. PLoS Negl Trop. Drosophila Inf Service7: e1248.10.1371/journal.pntd.0001248 MietzeA, MorickD, KöhlerH, HarrusS, DehioCet al. (2011) Combined MLST and AFLP typing of Bartonella henselae isolated from cats reveals new sequence types and suggests clonal evolution.Vet Microbiol148: 238-245.10.1016/j.vetmic.2010.08.012 PubMed: 2086363120863631ChungCY, KastenRW, PaffSM, Van HornBA, Vayssier-TaussatMet al. (2004) Bartonella spp. DNA associated with biting flies from California. Emerg Infect Dis10: 1311-1313 PubMed: 1532455715324557TsaiYL, ChomelBB, ChangCC, KassPH, ConradPAet al. (2011) Bartonella and Babesia infections in cattle and their ticks in Taiwan. Comp Immunol Microbiol Infect Dis34: 179-187.10.1016/j.cimid.2010.11.003 PubMed: 2119475021194750MargosG, GatewoodAG, AanensenDM, HanincováK, TerekhovaDet al. (2008) MLST of housekeeping genes captures geographic population structure and suggests a European origin of Borrelia burgdorferi. Proc Natl Acad Sci U S A105: 8730-8735.10.1073/pnas.0800323105 PubMed: 18574151 18574151