Skip to main content

Epidemiology of Bartonella henselae infection in pet and stray cats in Croatia with risk factors analysis

Abstract

Background

Cats are the primary reservoirs of the bacterium Bartonella henselae, the main cause of cat-scratch disease in humans. The main vector of the bacterium is the cat flea, Ctenocephalides felis. In southeastern Europe, data are lacking on the prevalence of B. henselae infection in cats, the strains of B. henselae involved and the risk factors associated with the infection.

Methods

Blood samples collected in ethylenediaminetetraacetic acid-containing tubes from 189 domestic cats (156 pet cats and 33 stray cats) from Zagreb, the capital city of Croatia, and 10 counties throughout Croatia were cultured for Bartonella spp. Following culture, bacterial isolates were genotyped at eight loci after using PCR to amplify 16S ribosomal RNA (rRNA) and the internal transcribed spacer region between the 16S and 23S rRNA sequences. Univariate and multivariate logistic regression were used to identify risk factors for B. henselae infection in cats.

Results

Bartonella spp. was detected in 31 cats (16.4%), and subsequent genotyping at the eight loci revealed B. henselae in all cases. Thirty complete multilocus sequence typing profiles were obtained, and the strains were identified as four sequence types that had been previously reported, namely ST5 (56.7%), ST6 (23.3%), ST1 (13.3%) and ST24 (3.3%), as well as a novel sequence type, ST33 (3.3%). The univariate analysis revealed a significantly higher risk of B. henselae infection in cats residing in coastal areas of Croatia (odds ratio [OR] 2.592, 95% confidence interval [CI] 1.150–5.838; P = 0.0191) and in cats with intestinal parasites (OR 3.207, 95% CI 1.088–9.457; P = 0.0279); a significantly lower risk was identified in cats aged > 1 year (OR 0.356, 95% CI 0.161–0.787; P = 0.0247) and in cats sampled between April and September (OR 0.325, 95% CI 0.147–0.715; P = 0.005). The multivariate analysis that controlled for age showed a positive association with the presence of intestinal parasites (OR 4.241, 95% CI 1.243–14.470; P = 0.0119) and coastal residence (OR 2.567, 95% CI 1.114–5.915; P = 0.0216) implying increased risk of infection, and a negative association with sampling between April and September (OR 0.379, 95% CI 0.169–0.848; P = 0.018) implying a decreased risk of infection. After controlling for the season, an increased risk of infection remained for the coastal region (OR 2.725, 95% CI 1.200–6.186; P = 0.012).

Conclusions

Bartonella henselae is prevalent throughout Croatia and is a public health threat. Environmental and host factors can significantly affect the risk of infection, and these should be explored in more detail. The presence of intestinal parasites highlights the need to eliminate the flea vector, Ctenocephalides felis, as the most effective approach to control infections in cats and humans.

Graphical Abstract

Background

Bartonella spp. are slow-growing alpha-proteobacteria [1], with more than 40 species and 20 Candidatus [2,3,4], of which at least 13 are zoonotic [5, 6]. These bacteria are well adapted to their mammalian animal reservoirs, in which they usually cause long-term intraerythrocytic bacteremia that is mostly asymptomatic [3,4,5, 7, 8]. Several species of Bartonella have been detected in the blood of domestic cats (Felis catus), of which Bartonella henselae is the dominant infective species. However, cats also serve as natural reservoirs for less common Bartonella spp., such as B. clarridgeiae and B. koehlerae, and as accidental hosts for B. bovis, B. quintana and B. vinsonii subsp berkhoffii [3, 5, 9, 10, 15, 61]. Bartonella spp. are mainly transmitted through infected blood-sucking arthropod vectors, although other modes of transmission, such as feline blood transfusions, have also been observed [5]. The major competent vector of B. henselae is the cat flea, Ctenocephalides felis [10,11,12,13,14]. In both cats and humans, infection primarily occurs following a cat scratch, since multiplied bacteria from flea feces contaminate the cat's claws and then enter the injured skin [2, 14, 17]. Humans exposed to infected cats can develop cat-scratch disease [3, 4, 10, 15, 16]. Ticks are also recognized as possible vectors of Bartonella spp. [3, 5, 10], but there is currently no consensus on the ability of Ixodes ricinus ticks to transmit B. henselae to cats [8, 12,13,14,15, 70].

The prevalence of Bartonella spp. infection in cats appears to vary widely around the world [5, 10, 18], ranging from 1.6% among pet cats in Canada [19] to 62% among stray cats in France [20], as determined by culture, or from 4.7% in stray cats in Greece [64] to 39.9% in sheltered cats in Brazil [60], as determined by PCR. Prevalence therefore depends not only on the diagnostic methods used, which vary in sensitivity [15, 35, 57], but also on the geographic area, characteristics of the cat population and living conditions, all of which appear to strongly influence the risk of feline infection [5, 10, 18]. Frequently published risk factors for Bartonella spp. infection in cats include young age and flea infestation [21, 22], living outdoors or as a stray [19,20,21, 23, 24], warm climate [22] and lack of prophylaxis against ectoparasites [24,25,26]. In the USA, the risk of cats having bacteria in the bloodstream was 80% lower in cats aged ≥ 13 months than in cats aged up to 6 months (odd ratio [OR] 0.18, 95% confidence interval [CI] 0.05–0.61) [22]; the same lower risk was found in cats that lived strictly indoors compared to outdoor cats (OR 0.18, 95% CI 0.06–0.54) [21]. In Poland, cats from the warmest region of the country were found to have a 2.01-fold higher risk of Bartonella spp. infection than cats from cooler regions (OR 2.01, CI 1.16–3.47). Similarly, cats that were not consistently treated against ectoparasites had a 2.02-fold higher risk of bacteremia, with a 95% CI of 1.41–2.92 [26]. The risk for B. henselae infection in stray cats in China was reported to be more than double that of client-owned cats (OR 2.283, 95% CI 1.093–4.772) [24]. Additionally, a study in the USA reported that the risk of Bartonella spp. infection was threefold higher in cats with fleas than in uninfested cats (OR 2.82, 95% CI 1.1–7.3) [22]. The broad geographic variation in prevalence of Bartonella spp. infection in cats and the potentially widespread distribution of risk factors for this infection highlight the need for each country or region to assess these parameters locally in order to guide control and monitoring efforts.

Data are lacking on the prevalence and genetic diversity of B. henselae in cats in Croatia and most of the surrounding countries of Central and Southeastern Europe [5, 27]. The bacterium is known to infect cats in Croatia [28,29,30] and to circulate in humans [31,32,33]. Therefore, the aim of the present study was to assess the prevalence of Bartonella spp. in cats based on multilocus sequence typing (MLST), and to examine the risk factors for such infection.

Methods

Sampling of cats

For this study, we analysed blood samples collected in ethylenediaminetetraacetic acid-containing tubes (EDTA-blood samples) between 2014 and 2017 from 189 domestic cats (156 pet cats, 33 stray cats) that had been examined or undergone minor surgical procedures at 20 veterinary clinics at 13 locations around Croatia (nine inland locations, four locations along the coast). The majority of blood samples included in the analysis were collected from apparently healthy cats (N = 149), while 40 (21%) cats showed clinical signs.

Veterinarians who agreed to participate in the study were sent a questionnaire that had been specially designed to provide information that could shed light on risk factors for Bartonella spp. infections in cats. Cat owners were informed of the purpose of sampling and agreed to provide the requested information on the animals and on the presence/absence of cat-scratch disease in their homes. Age of the cat at the time of blood sampling was recorded by veterinarians based on information provided by the cat owners or based on the veterinarian’s estimate (see Table 2 for data on animal health, treatments, living and environmental conditions, and other data collected by means of the questionnaires and selected for the risk factor analysis).

Venepuncture was performed aseptically, and blood was collected directly into commercial Vacuette tubes (Greiner Bio-One, Kremsmünster, Austria) containing ethylenediaminetetraacetic acid dipotassium salt dihydrate (K2-EDTA) as anticoagulant. Samples were stored frozen until their analysis at the Croatian Veterinary Institute (Zagreb, Croatia) for the present study.

Culture of Bartonella spp. in cat blood

Prior to feline blood inoculation, we conducted initial tests on the culture media using two reference strains of cat-related species (B. henselae ATCC 49882 and B. clarridgeiae ATCC 700095), purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA). These tests aimed to assess laboratory conditions, growth quality and colony recognition. Since solid media were our first choice, aliquots of each of the 189 samples of thawed and well-mixed cat’s blood were inoculated directly onto the surface of the specific agar medium and the plates rotated to ensure an even distribution. All agar media used included defibrinated sheep or rabbit blood: tryptic soy agar [36], brain–heart agar [37], Columbia agar [20, 38], chocolate agar [22, 39] and previously unexplored esculin-blood agar (No. 2 blood agar base, 5% sheep blood and 0.275 g/l esculin powder, prepared as Koch-sterilized 1% aqueous solution and added to the medium at 27.5 ml per liter). The volume of blood in the sample, which ranged from 0.5 to 3.0 ml, determined how many agar plates per sample could be used; thus, each blood sample was inoculated onto two to five different types of agar media to increase the chances of obtaining isolates.

For a small number of blood samples there was enough blood for additional inoculations, using one or two types of biphasic media consisting of a liquid medium layered on an agar slant in a test tube: tryptic soy agar with tryptic soy broth [1] and brain–heart agar [37] with Brucella broth [38]. All media were freshly prepared from dehydrated bases manufactured by Merck (Darmstadt, Germany), with the exception of the Brucella broth (BD; Sparks, MD, USA).

Inoculated media were incubated for at least 4–8 weeks at 37 °C in a humid atmosphere containing 8% CO2 and checked for bacterial growth twice a week. Agar surfaces were inspected for typical Bartonella spp. colonies, which appear as single, hard colonies that are firmly embedded in agar [9, 35]. The isolation rate for each medium was calculated by dividing the number of B. henselae isolates obtained on that type of medium with the total number of that medium inoculated with the blood of infected cats.

Multilocus sequence typing

Primary or subcultured colonies from solid and biphasic media were subjected to multilocus sequencing as follows. Bacterial colonies were mixed with DNA extraction solution containing 180 µl of lysis buffer (“Buffer AL”) and 20 µl of proteinase K from the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany). The mixtures were incubated at 56 °C for 20 min until lysis was complete, then centrifuged at 6000 g for 1 min and the supernatant processed using the QIAcube automated DNA extraction system (Qiagen) according to the manufacturer's protocol for blood and body fluids. As an isolation control, 200 μl of phosphate-buffered saline (PBS) was processed in parallel. The isolated DNA from each sample (200 µl) was stored at − 20 °C until PCR testing.

Conventional PCR testing was conducted on aliquots of extracted DNA, including the isolation control sample, using the Pro Flex PCR System (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA) and previously described primers [69] to amplify regions of the 16S ribosomal RNA (rRNA) and the internal transcribed spacer between the 16S and 23S rRNA. Analysis of the 16S rRNA allows differentiation of genotype I (“Houston-1”) and genotype II (“Marseille”) of B. henselae [20, 41]. Each 20-µl PCR reaction volume contained 10 µl HotStarTaq Master Mix solution (Qiagen), 7 µl of water free of DNase and RNase (Qiagen), 0.5 µM of each primer (Macrogen, Amsterdam, The Netherlands) and template DNA or PBS (2 µl). During each PCR run, there was both an amplification negative control containing 2 µl DNase/RNase-free sterile water and a positive control containing 2 µl of DNA purified from the reference strain B. henselae Huston-1 (ATCC 49882). The thermal cycling conditions were: one cycle at 95 °C for 15 min; followed by 40 cycles of 95 °C for 45 s, 55 °C for 45 s and 72 °C for 45 s; with a final extension at 72 °C for 10 min, followed by a hold at 4 °C. PCR amplicons were visualized on a capillary electrophoresis device (QIAxcel System; Qiagen) and stored at − 20 °C for 1 to 24 months, until sequence typing.

To prevent contamination, we performed the DNA extraction, preparation of PCR reagents, addition of DNA template to PCR reactions, running of the PCR reactions and capillary electrophoresis in separate, dedicated rooms using disposable pipette tips, gloves and aprons.

The DNA amplified from isolates was subjected to sequence typing at eight gene loci, 16S rRNA and genes involved in protein production: batR (coding for two-component regulator), ftsZ (coding for cell division protein), gltA (coding for citrate synthase), groEL (coding for heat shock protein 60 chaperon), nlpD (coding for cell surface glycoprotein), ribC (coding for riboflavin synthase) and rpoB (coding for RNA polymerase subunit), as previously described [42, 43]. This widely used typing method has been used to define 37 sequence types of B. henselae differing in virulence and zoonotic potential [18, 27, 44, 45]. The DNA was purified using ExoSAP-IT® (USB, Cleveland, OH, USA) according to the manufacturer's instructions, then sequenced by Macrogen. The sequences obtained were used to define alleles and sequence types based on comparisons with the PubMLST database containing all reported B. henselae sequence types (https://pubmlst.org/organisms/bartonella-henselae) [27].

Phylogenetic analysis of sequence types

Phylogenetic analysis was performed by comparing sequences at the eight target genes using BioNumercis 7.6 (Applied Maths, Gent, Belgium) and data in the PubMLST database. The unweighted pair group method with arithmetic mean (UPGMA) was used to cluster the concatenated allelic sequences of all sequence types identified in this study, together with 36 previously reported sequence types, and to construct a “minimum-spanning” phylogenetic tree (MST) showing their relatedness and evolutionary relationships.

Statistical analysis

Data were processed statistically using STATA 13.1 (StataCorp, College Station, TX, USA). The normality of data distribution was assessed based on plots and the non-parametric Shapiro–Wilk test. Results of statistical tests were considered significant if associated with P < 0.05.

Univariate and multivariate logistic regressions were performed to identify risk factors for B. henselae infection in cats. Since the dependent variable in our study is binary (B. henselae-positive or -negative cat), we used logistic regression also known as binominal logistic regression. Variables that emerged as significant from the univariate analysis were entered step-wise into the multivariate regression. Where appropriate, results were expressed as ORs with 95% CIs.

Ethics approval

This study was approved by the Committee for Ethics in Veterinary Medicine (640-01/16-17/63, 251-61-01/139-16-2). Cat owners who provided information to assist veterinarians to complete the study questionnaire were informed of the purposes of this research, and they consented for the anonymized information to be analysed and published.

Results

Blood from 189 cats ranging in age from 5 months to 16 years were sampled at veterinary clinics around Croatia (Table 1). Most of the cats in our sample (N = 156, 82.5%) were pets, and the remaining 33 cats were street cats (4.2%) or from shelters (13.2%) (Table 2). Bartonella spp. was detected in 31 animals (16.4%, 95% CI 11.1–21.7) from nine of 13 locations, and prevalence at individual locations ranged from 8.8% to 66.7%. Infections in animals were not detected at four distant locations, which were spread throughout the country (Table 1).

Table 1 Prevalence of Bartonella henselae in 189 cats sampled in Croatia by sampling locations and geographical distribution of sequence types
Table 2 Univariate analysis to identify risk factors for B. henselae infection in cats (N = 189)

Bartonella henselae isolates grew more frequently on solid culture media (75.5%; 77/102) and less frequently on biphasic culture media (60.0%; 6/10). The best isolation efficiency was achieved on brain–heart agar (87.5%) and Columbia agar (82.4%), followed by esculin-blood agar (80.0%), biphasic tryptic soy agar and broth (80.0%) and chocolate agar (77.3%). The least effective culture media were tryptic soy agar (54.2%) and biphasic brain–heart agar with Brucella broth (40.0%). The first colonies appeared in primary isolation after 4 to 56 days of incubation.

All isolates were found to be B. henselae, and 30 complete multilocus sequencing profiles were obtained that belonged to five sequence types (STs): ST5 (accounting for 17 of the 30 profiles; 56.7%); ST6 (7 profiles; 23.3%); ST1 (4 profiles; 13.3%); ST24 and ST33 (each 1 profile; 3.3%). It was not possible to identify one isolate (“DUB5”) because we failed to obtain usable sequence data for the gltA gene, which meant we could not determine whether it belonged to ST5 or ST24. While four of the sequence types that we detected have previously been reported, ST33 appears to be a new sequence type and was assigned the numerical code 2-3-3-1-2-1-1-2 based on the gene order 16S, batR, ftsZ, gltA, groEl, nlpD, ribC and rpoB. This new B. henselae isolate was entered in PubMLST under the name CRO_MICIKA and accession ID 447.

The five sequence types were distributed across eight locations in Croatia (Table 1). The three most prevalent sequence types (ST1, ST5 and ST6) were detected in samples collected at both inland and coastal locations. ST5 occurred more often in samples collected at inland locations (64.7%) than in those collected at coastal locations (46.2%). Conversely, ST6 occurred more often in samples collected at coastal locations (38.5% vs 11.8%). Both ST24 and ST33 occurred in samples collected at the inland locations of Bjelovar and Jastrebarsko, respectively. The unidentified isolate DUB5 was detected in a sample collected at a coastal location (Dubrovnik).

We were able to assign all the sequence types in our study to one of three previously identified clonal complexes of B. henselae [43, 44]. The MST (Fig. 1) showed ST33 to be closely related to an isolate of ST26 in Germany [46] and of ST27 in the UK [44]. Adding ST33 to clonal complex 2 increased the cluster from eight to nine members.

Fig. 1
figure 1

“Minimum-spanning” phylogenetic tree showing putative relationships among B. henselae sequence types, based on the unweighted pair group method with arithmetic mean. Isolates are represented as circles, whose size reflects their frequency in the data pooled from the present study and the public database PubMLST, and whose color reflects the country of origin. Numbers refer to sequence types. Thick lines indicate minor allelic differences (in only 1 of 8 genes analysed) and therefore closer relatedness. Thin lines indicate sequence variations in 2 or 3 genes. The sequence type ST33, reported in the present study for the first time, is marked in the middle of the tree

In the univariate analysis age > 1 year was associated with significantly lower risk of infection than age ≤ 1 year (OR 0.356, 95% CI 0.161–0.787; P = 0.0247), with 19 of 76 cats (25%) in the younger age group showing infection, compared to only 12 of 113 (10.6%) older cats. Risk of infection in coastal regions was 2.6-fold higher than the risk in inland regions (OR 2.592, 95% CI 1.150–5.838; P = 0.0191), and risk among cats with intestinal parasites was more than threefold higher than the risk among those without parasites (OR 3.207, 95% CI 1.088–9.457; P = 0.0279). The risk of infection was almost 70% lower in cats sampled between April and September (OR 0.325, 95% CI 0.147–0.715; P = 0.005) than in cats sampled between October and March. The difference in prevalence between sampling sites (OR 1.129, 95% CI 1.021–1.249; P = 0.007), and between cats by age (OR 0.985, 95% CI 0.974–0.997; P = 0.0047) was also statistically significant, although weakly associated (Table 2).

We then stepwise included the five factors into the multivariate models and controlled for each factor in turn (Table 3). The association between infection and presence of intestinal parasites remained significant and became even stronger after controlling for age (OR 4.241, 95% CI 1.243–14.470; P = 0.0119), but it lost significance after controlling for other factors. The association between infection and coastal residence remained similarly significant after controlling for age (OR 2.567, 95% CI 1.114–5.915; P = 0.0216) and season (OR 2.725, 95% CI 1.200–6.186; P = 0.012). The association between less frequent infection and sampling between April and September remained significant after controlling for age (OR 0.379, 95% CI 0.169–0.848; P = 0.018). The significant association between infection and sampling locations remained weak after controlling for age (OR 1.125, 95% CI 1.014–1.249; P = 0.0257).

Table 3 Multivariate analysis to identify independent predictors of B. henselae infection in cats (N = 189)

Discussion

We detected a prevalence of Bartonella spp. of 16.4% in 189 domestic cats in Croatia using the culture method, similar to the 17% reported for domestic cats in Auckland, New Zealand and to the 16.5% reported for domestic cats in Paris, France [23, 47], but less than the 39.5% reported for pet and stray cats in northern California State, USA [21]. Similar to the studies in New Zealand and the USA, we detected only B. henselae in our study animals, while the study in France also detected a small prevalence of B. clarridgeiae. Detection of B. henselae as the sole species in present study is not surprising since other Bartonella spp. associated with cats, such as B. clarridgeiae and B. koehlerae, are less common and challenging to isolate by culture [3, 5, 9, 10, 20, 22].

The relatively low prevalence of Bartonella spp. in our study may reflect, in part, that > 80% of the cats in our sample were pets living with owners and, therefore, more likely to receive better nutrition and flea treatment [24, 61]. In comparison, a study in Spain involving a sample similar to that in our study, with 87% of pet cats living with owners in the former study, reported a low prevalence of Bartonella spp. of 7% [38]. Other European studies that sampled cats brought to veterinary clinics reported an even lower prevalence, such as 2.2% in Germany and 5.8% in the UK [44, 46]. In contrast, studies involving higher proportions of stray cats than that in our study reported a higher prevalence, such as 22.6% in Denmark [48] or 28.1% in Turkey [49], while studies of only stray cats in France reported a prevalence exceeding 50% [20, 50].

Although access to the outside world intensifies contacts with other cats and ectoparasites, neither status as a pet versus stray nor as indoor versus outdoor living appeared to significantly affect the risk of infection in our sample. In contrast, studies in several other countries have identified one or both of these as risk factors [18, 21, 25, 26, 53, 62]. We also failed to detect an association between infection and flea infestation [22, 23, 49, 60], or lack of prophylactic flea/tick treatment [24,25,26, 61], which could also be the result of more appropriate care when cats have owners.

Our analysis among cats in Croatia detected several factors as being associated with higher risk of infection that have also been reported in several other countries. We found a higher infection risk among cats aged up to 1 year, which is consistent with a greater risk among younger animals reported in previous studies [21, 22, 36, 49, 50]. The exact reasons behind the increased susceptibility to infections of kittens and young cats remain unclear. In addition to a relatively weaker immune systems in young animals, there is evidence that fleas have a preference for young cats [71]; in addition, there have been recent suggestions of transmission from mothers to offspring [72]. Furthermore, recent studies have shown that antibodies provide sufficient protection against Bartonella spp. by preventing their attachment to red blood cells [73], but also by long-term prevention of reinfection with the same strain [74]. This supports earlier hypotheses that younger animals could lack the protection due to their undeveloped antibodies [15, 21, 50], although more verification of this possibility is required.

Our finding that cats from coastal regions have a greater risk of infection than those from inland areas is consistent with the higher prevalence of Bartonella spp. in warmer areas. For example, in the USA, prevalence of Bartonella spp. in cats was significantly higher in California and Florida than in Chicago, Washington DC or Michigan [22, 56]. In eastern Poland, prevalence of Bartonella spp. in cats was found to be substantially higher in the warmer Subcarpathian region than in cooler areas [26]. These observations may reflect the fact of warmer habitats being more hospitable for the fleas that may carry the bacteria.

We detected a higher risk of Bartonella spp. infection in cats between October and March than between April and September, which mostly overlaps with the known seasonality of cat-scratch disease, which in Croatia usually occurs between August and March or in France between September and April [32]. In the present study, for example, we isolated B. henselae from three cats in January, March and October, when members of the same households were also diagnosed with cat-scratch disease [28,29,30], which fully matches aforementioned seasonal pattern. Since cats are more often indoors when the weather turns cold [32], and because they can remain bacteremic long after infection [32], the chances of infection might be increased during the fall and winter.

Interestingly, none of the 53 cats that were sampled in June, July and August during three consecutive years in our study showed Bartonella spp. infection, although fleas are generally most abundant during the summer [71]; this result is consistent with the lower prevalence of infection in the summer months seen in a previous study from the USA [57]. Low infection rates in the summer may reflect an increased use of flea protection products, which can prevent B. henselae transmission [58, 59], or perhaps low and currently undetectable levels of bacteremia due to the cyclical circulation pattern in cats [4, 5, 29, 57]. In the southern hemisphere, in Brazil, a significantly higher risk for infection was also found in cats in autumn, after the end of the hot summer, possibly reflecting a persistent activity of fleas or a flea treatment has been missed [60]. Therefore, knowledge of the seasonal nature of human and feline infections may be helpful to clinicians when cat-scratch disease is suspected.

In contrast to these risk factors reported above, which are common to Croatia and several other countries, our study detected one risk factor that has yet to be reported: the presence of intestinal parasites. Such an association is new but not surprising because cats with intestinal parasites are more likely to have the cat flea C. felis [63], which also serves as an intermediate host for the tapeworm Dipylidium caninum [63, 64]. Indeed, four of the six cats with fleas and B. henselae infection in our study also had the intestinal parasites D. caninum or Toxocara cati. Therefore, the practice of simultaneously treating of cats for external and internal parasites is strongly recommended. The association between intestinal parasites and B. henselae infection has not yet been discussed, as we are aware of only two studies that have analysed intestinal parasites as potential risk factors [38, 60]. Endoparasites, particularly in young cats, exacerbate an already weaker immune system and increase the likelihood of Bartonella spp. infections. As five of these six cats with fleas and B. henselae infection in our study were pets, the worrying implications for public health should be addressed; both T. cati and D. caninum are parasites of zoonotic importance, with a real possibility of transmission to humans [63, 64]. Since intestinal parasites in our study were detected through visual observation of the perianal area or feces, our findings should be verified and extended in studies based on coprology. These considerations highlight the need to include intestinal parasites when analysing risk factors of Bartonella spp. infection in cats.

The most reliable method for identifying active Bartonella spp. infection is inoculation of cat blood onto solid culture media, followed by molecular typing of isolated bacteria [3, 15, 34, 35]. The use of agar plates for cultivation of Bartonella spp. is generally the preferred diagnostic method due to its ability to yield a sufficient quantity and higher quality of DNA for complex genetic investigations [35, 40], such as MLST. We detected five sequence types of B. henselae in the Croatian cats included in our study, consistent with the detection of diverse strains of B. henselae in cats in several other European countries [18, 44, 46] and elsewhere in the world [42, 43, 51]. The three sequence types that together explained 93% of infections in our sample (ST5, ST6 and ST1) have also been frequently detected in cats from northwestern and Mediterranean Europe [43, 44].

The most frequent sequence types in our study, ST5 (56.7%) and ST6 (23.3%), were also reported to be the most prevalent in Spain (61.5% and 15.4%, respectively), but ST1 was not detected in cats in the Spanish study despite its detection in human isolates [18]. ST5 has also been reported to be quite prevalent in cats from Italy, France, Greece, Germany and the USA [43, 46, 51], and has also been identified in cats in the UK [44], Australia [43] and Algeria (Africa) [37]. In contrast, ST6 was found to be more prevalent (40.0%) in the UK and Australia [42, 44], and also identified in Italy and the USA [43], while it was not detected in Germany [46].

Although ST1 was the third most prevalent sequence type in our study (13.3%), this subtype is actually the most common feline and human sequence type worldwide, having been detected in Europe, Australia, North and South America, Asia and Africa [37, 39, 42,43,44, 51,52,53,54]. ST1 is quite common among cats from the Mediterranean area of Europe (Italy, Greece) [43, 51] and Asia (Israel) [43], but less prevalent in Germany and the UK [43, 44]. Thus, our detection of ST1 and ST5 seems to fit with a potentially widespread presence of these sequence types throughout the Mediterranean region, which remains to be verified in the future.

The ST24 that we detected was previously reported only in Germany [46] and Algeria [37], while ST33 appears to be novel. Since sequence types are thought to develop locally through clonal evolution [46], it is unclear whether ST33 arose in the inland area of Jastrebarsko, where the infected cat lived during sampling, or in the coastal town of Primošten, where it was previously adopted.

The sequence types in our study did not show obvious geographical clustering within Croatia, with the exception that ST6 was observed to occur more frequently at coastal locations (38.5%) than at inland ones (11.8%). Geographically biased distribution was observed in Spain, with ST6 isolates localized to only one of two regions with a similar climate, making it unclear whether such clustering reflects climate or other causes [18]. In contrast, locational clustering of sequence types was not observed in the UK [44].

We also identified diverse sequence types within relatively small areas (approx. 50 km2) of cities, such as areas of the cities of Jastrebarsko (ST1, ST5, ST6 and ST33) and Pula (ST1, ST5 and ST6). Diversity was observed even between cats within the same households in the cities of Jastrebarsko and Rijeka (ST5 and ST6), but also within the same shelter in the city of Pula (ST1 and ST6). Such local diversity may reflect the strains carried by fleas. Since cats by nature tend to move only over a limited area, it is assumed that the transmission of B. henselae is also quite local [46, 49]. This possibility may also help explain the detection of nine sequence types in Hannover, Germany [46], which appears to be an example of the greatest variability of sequence types in one city. In contrast, infected cats and humans in Asia contain a smaller diversity of strains, dominated by ST1 [39, 52, 53]. Future studies should examine whether these observations reflect differences in local flea populations that presumably infect the cats. However, we are aware of only one study involving fleas conducted in Spain which has demonstrated fleas as carriers of ST5 [18].

ST5, as the most frequent sequence type in our study (56.7%), is also the one most frequently associated with cat-scratch disease patients, particularly in Europe [18, 43, 44]. A study in Spain, for example, identified this sequence type in > 50% of infected humans, mostly those with the typical form of cat-scratch disease [18]. In the USA, Lindroos et al. [51] detected ST5 in two cats and in samples from their owners with cat-scratch disease, implying transmission from pets to humans. We also confirmed the presence of ST5 in three cats from households with patients diagnosed with cat-scratch disease [28,29,30]; although the sequences types were not analysed in the affected families, cat-to-human transmission was strongly suspected. All three cat-owner dyads were typical of cat-scratch disease: the cats were young (5 months, 10 months and 2 years), asymptomatic, and had a history of flea infestation. Both cats and humans were diagnosed in the autumn and winter. Two of the three patients were minors (12 and 16 years old), consistent with the greater prevalence of cat-scratch disease among children and young people [10, 16, 31, 32, 55].

In contrast, the presence of ST1 in Spain appeared to be more associated with clinically atypical B. henselae infections in humans, such as endocarditis, fever of unknown origin and hepatic peliosis [18]. ST6 has been linked to cat-scratch disease in Australia [42] and France [43]. However, studies in the UK and Spain failed to detect this sequence type in human patients, even though it was detected in cats [18, 44]. Other members of the cluster to which ST6 belongs (clonal complex 2), including the apparently novel ST33, have only been found in cats, but not yet in humans. This coincides with the knowledge obtained from the older division based on 16S rRNA gene analysis, which confirms that all sequence types from our study (ST5, ST6, ST24, ST33) with the exception of ST1 belong to B. henselae genotype II, which has been more frequently detected in cats from Europe, USA and Australia than in humans [18, 21,22,23, 42, 44, 46, 48, 62]. Whether this suggests a lower virulence of the new ST33 for humans remains to be seen in future studies.

Our study describes, apparently for the first time, the culture of Bartonella spp. using esculin-blood agar. Widely used to culture mastitis-causing bacteria, such as staphylococci, streptococci and coliforms in diary animals [65, 66], this medium is also used routinely in our laboratories on milk samples, food of animal origin and as a general-purpose bacteriological medium, as well as in scientific studies [67, 68]. It has also proved to be suitable for primary isolation of the fastidious bacterium B. henselae through culturing of feline blood samples [30], which represents a novelty.

Conclusions

We have detected B. henselae in 31 cats (24 pet cats and 7 stray cats) from across Croatia, suggesting that the pathogen is widespread throughout the country, posing a public health threat to humans. The situation may be similar in other parts of Southeast Europe, which should be explored in future work. The association between intestinal parasites and B. henselae in cats highlights the need to eliminate the flea vector, C. felis, as the most effective approach to controlling infections in cats and humans.

Availability of data and materials

All data generated or analysed during this study are included in this article and are available from the corresponding authors upon the request.

Abbreviations

16S rRNA:

16S ribosomal RNA

batR :

Gene coding for the two-component regulator

CI:

Confidence interval

ftsZ :

Gene coding for the cell division protein

gltA :

Gene coding for the citrate synthase

groEL :

Gene coding for the heat shock protein 60 (Hsp60) chaperone

K2-EDTA:

Ethylenediaminetetraacetic acid dipotassium salt dihydrate

MLST:

Multilocus sequence typing

MST:

Minimum spanning tree

nlpD :

Gene coding for the cell surface glycoprotein

OR:

Odds ratio

PBS:

phosphate buffered saline

ribC :

Gene coding for the riboflavin synthase

rpoB :

Gene coding for the RNA polymerase beta subunit

ST:

Sequence type

UPGMA:

Unweighted pair group method with arithmetic mean

References

  1. Diddi K, Chaudhry R, Sharma N, Dhawan B. Strategy for identification & characterization of Bartonella henselae with conventional & molecular methods. Indian J Med Res. 2013;137:380–7.

    PubMed  PubMed Central  Google Scholar 

  2. Okaro U, Addisu A, Casanas B, Anderson B. Bartonella species, an emerging cause of blood-culture-negative endocarditis. Clin Microbiol Rev. 2017;30:709–46. https://doi.org/10.1128/CMR.00013-17.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Breitschwerdt EB. Bartonellosis, One Health and all creatures great and small. Vet Dermatol. 2017;28:96-e21. https://doi.org/10.1111/vde.12413.

    Article  PubMed  Google Scholar 

  4. Deng H, Pang Q, Zhao B, Vayssier-Taussat M. Molecular mechanisms of Bartonella and mammalian erythrocyte interactions: a review. Front Cell Infect Microbiol. 2018;8:431. https://doi.org/10.3389/fcimb.2018.00431.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Álvarez-Fernández A, Breitschwerdt EB, Solano-Gallego L. Bartonella infections in cats and dogs including zoonotic aspects. Parasit Vectors. 2018;11:624. https://doi.org/10.1186/s13071-018-3152-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Gadila SKG, Embers ME. Antibiotic susceptibility of Bartonella grown in different culture conditions. Pathogens. 2021;10:718. https://doi.org/10.3390/pathogens10060718.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Jacomo V, Kelly PJ, Raoult D. Natural history of Bartonella infections (an exception to Koch’s postulate). Clin Diagn Lab Immun. 2002;9:8–18. https://doi.org/10.1128/CDLI.9.1.8-18.2002.

    Article  CAS  Google Scholar 

  8. Chomel BB, Boulouis HJ, Breitschwerdt EB, Kasten RW, Vayssier-Taussat M, Birtles RJ, et al. Ecological fitness and strategies of adaptation of Bartonella species to their hosts and vectors. Vet Res. 2009;40:29. https://doi.org/10.1051/vetres/2009011.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Boulouis HJ, Chang CC, Henn JB, Kasten RW, Chomel BB. Factors associated with the rapid emergence of zoonotic Bartonella infections. Vet Res. 2005;36:383–410. https://doi.org/10.1051/vetres:2005009.

    Article  PubMed  Google Scholar 

  10. Chomel BB, Kasten RW. Bartonellosis, an increasingly recognized zoonosis. J Appl Microbiol. 2010;109:743–50. https://doi.org/10.1111/j.1365-2672.2010.04679.x.

    Article  CAS  PubMed  Google Scholar 

  11. Chomel BB, Kasten RW, Floyd-Hawkins K, Chi B, Yamamoto K, Roberts-Wilson J, et al. Experimental transmission of Bartonella henselae by the cat flea. J Clin Microbiol. 1996;34:1952–6. https://doi.org/10.1128/jcm.34.8.1952-1956.1996.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Billeter SA, Hayman DTS, Peel AJ, Baker KS, Wood JLN, Cunningham AA, et al. Bartonella species in bat flies (Diptera: Nycteribiidae) from western Africa. Parasitology. 2012;139:324–9. https://doi.org/10.1017/S0031182011002113.

    Article  CAS  PubMed  Google Scholar 

  13. Kosoy M, Hayman DT, Chan KS. Bartonella bacteria in nature: where does population variability end and a species start? Infect Genet Evol. 2012;12:894–904. https://doi.org/10.1016/j.meegid.2012.03.005.

    Article  PubMed  Google Scholar 

  14. Bouhsira E, Franc M, Boulouis HJ, Jacquiet P, Raymond-Letron I, Liénard E. Assessment of persistence of Bartonella henselae in Ctenocephalides felis. Appl Environ Microb. 2013;79:7439–44. https://doi.org/10.1128/AEM.02598-13.

    Article  CAS  Google Scholar 

  15. Guptill L. Bartonellosis. Vet Microbiol. 2010;140:347–59. https://doi.org/10.1016/j.vetmic.2009.11.011.

    Article  PubMed  Google Scholar 

  16. Nelson CA, Saha S, Mead PS. Cat-scratch disease in the United States, 2005–2013. Emerg Infect Dis. 2016;22:1741–6. https://doi.org/10.3201/eid2210.160115.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Iannino F, Sulli N, Maitino A, Pascucci I, Pampiglione G, Salucci S. Fleas of dog and cat: species, biology and flea-borne diseases. Vet Ital. 2017;53:277–88. https://doi.org/10.12834/VetIt.109.303.3.

    Article  PubMed  Google Scholar 

  18. Gil H, Escudero R, Pons I, Rodríguez-Vargas M, García-Esteban C, Rodríguez-Moreno I, et al. Distribution of Bartonella henselae variants in patients, reservoir hosts and vectors in Spain. PLoS One. 2013;8:e68248. https://doi.org/10.1371/journal.pone.0068248.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Kamrani A, Parreira VR, Greenwood J, Prescott JF. The prevalence of Bartonella, hemoplasma, and Rickettsia felis infections in domestic cats and in cat fleas in Ontario. Can J Vet Res. 2008;72:411–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. La Scola B, Davoust B, Boni M, Raoult D. Lack of correlation between Bartonella DNA detection within fleas, serological results, and results of blood culture in a Bartonella-infected stray cat population. Clin Microbiol Infect. 2002;8:345–51. https://doi.org/10.1046/j.1469-0691.2002.00434.x.

    Article  PubMed  Google Scholar 

  21. Chomel BB, Abbott RC, Kasten RW, Floyd-Hawkins KA, Kass PA, Glaser CA, et al. Bartonella henselae prevalence in domestic cats in California: risk factors and association between bacteremia and antibody titers. J Clin Microbiol. 1995;33:2445–50. https://doi.org/10.1128/jcm.33.9.2445-2450.1995.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Guptill L, Wu CC, Hogenesch H, Slater LN, Glickman N, Dunham A, et al. Prevalence, risk factors, and genetic diversity of Bartonella henselae infections in pet cats in four regions of the United States. J Clin Microbiol. 2004;42:652–9. https://doi.org/10.1128/JCM.42.2.652-659.2004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Gurfield AN, Boulouis HJ, Chomel BB, Kasten RW, Heller R, Bouillin C, et al. Epidemiology of Bartonella infection in domestic cats in France. Vet Microbiol. 2001;80:185–98. https://doi.org/10.1016/S0378-1135(01)00304-2.

    Article  CAS  PubMed  Google Scholar 

  24. Zhang Y, Zhang Z, Lou Y, Yu Y. Prevalence of hemoplasmas and Bartonella species in client-owned cats in Beijing and Shanghai. China J Vet Med Sci. 2021;83:793–7. https://doi.org/10.1292/jvms.20-0681.

    Article  PubMed  Google Scholar 

  25. Persichetti MF, Pennisi MG, Vullo A, Masucci M, Migliazzo A, Solano-Gallego L. Clinical evaluation of outdoor cats exposed to ectoparasites and associated risk for vector-borne infections in southern Italy. Parasit Vectors. 2018;11:136. https://doi.org/10.1186/s13071-018-2725-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Mazurek L, Carbonero A, Skrzypczak M, Winiarczyk S, Adaszek Ł. Epizootic situation of feline Bartonella infection in eastern Poland. J Vet Res. 2020;64:79–83. https://doi.org/10.2478/jvetres-2020-0019.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018;3:124. https://doi.org/10.12688/wellcomeopenres.14826.1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Stepanić M, Duvnjak S, Reil I, Špičić S, Kompes G, Jurković D, et al. Prvi dokaz sekvencijskog tipa 5 Bartonella henselae u mačaka: najvjerojatniji izvor zaraze za djecu oboljelu od bolesti mačjeg ogreba [First report of Bartonella henselae sequence type 5 in cats: probable source of infection for children with cat scratch disease]. Collection of abstracts. Workshop “Emergent and neglected zoonoses in the context of One Health", Zagreb; 2018. p. 81–82.

  29. Stepanić M, Duvnjak S, Reil I, Špičić S, Kompes G, Beck R. First isolation and genotyping of Bartonella henselae from a cat living with a patient with cat scratch disease in Southeast Europe. BMC Infect Dis. 2019;19:299. https://doi.org/10.1186/s12879-019-3929-z.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Stepanić M. Uzgoj i genska tipizacija bakterija roda Bartonella izdvojenih iz mačaka na području Republike Hrvatske [Bacterial culturing and genotyping of bacteria from the genus Bartonella isolated from cats in Republic of Croatia]. Dissertation. Faculty of Veterinary Medicine, University of Zagreb, Zagreb; 2022.

  31. Pandak N, Đaković-Rode O, Čabraja I, Krištof Ž, Kotarac S. Prevalence of Bartonella henselae antibodies in children and blood donors in Croatia. Infection. 2009;37:166–7. https://doi.org/10.1007/s15010-008-8113-0.

    Article  CAS  PubMed  Google Scholar 

  32. Vilibić-Čavlek T, Karlović-Martinković D, Ljubin-Sternak S, Tabain I, Peršić Z, Mlinarić-Galinović G. High prevalence of Bartonella henselae and Bartonella quintana antibodies in Croatian patients presenting with lymphadenopathy. Pol J Microbiol. 2012;61:315–7. https://doi.org/10.33073/pjm-2012-043.

    Article  PubMed  Google Scholar 

  33. Knežević S, Rončević-Filipović M, Gorup L, Cekinović-Grbeša Đ, Pavić I, Trošelj Vukić B, et al. Mnogobrojna lica bolesti mačjega ogreba—prikaz četiri bolesnika [The many faces of cat scratch disease—a report of four cases]. Infektol Glasn. 2020;40:159–64. https://doi.org/10.37797/ig.40.4.7.

    Article  Google Scholar 

  34. Fabbi M, De Giuli L, Tranquillo M, Bragoni R, Casiraghi M, Genchi C. Prevalence of Bartonella henselae in Italian stray cats: evaluation of serology to assess the risk of transmission of Bartonella to humans. J Clin Microbiol. 2004;42:264–8. https://doi.org/10.1128/JCM.42.1.264-268.2004.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Gutiérrez R, Vayssier-Taussat M, Buffet JP, Harrus S. Guidelines for the isolation, molecular detection, and characterization of Bartonella species. Vector-Borne Zoonotic Dis. 2017;17:42–50. https://doi.org/10.1089/vbz.2016.1956.

    Article  PubMed  Google Scholar 

  36. Koehler JE, Glaser CA, Tappero JW. Rochalimaea henselae infection. A new zoonosis with the domestic cat as reservoir. J Am Med Assoc. 1994;271:531–5. https://doi.org/10.1001/jama.1994.03510310061039.

    Article  CAS  Google Scholar 

  37. Azzag N, Haddad N, Durand B, Petit E, Ammouche A, Chomel B, et al. Population structure of Bartonella henselae in Algerian urban stray cats. PLoS One. 2012;7:e43621. https://doi.org/10.1371/journal.pone.0043621.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Pons I, Sanfeliu I, Quesada M, Anton E, Sampere M, Font B, et al. Prevalence of Bartonella henselae in cats in Catalonia, Spain. Am J Trop Med Hyg. 2005;72:453–7. https://doi.org/10.4269/ajtmh.2005.72.453.

    Article  PubMed  Google Scholar 

  39. Yanagihara M, Tsuneoka H, Hoshide S, Ishido E, Umeda A, Tsukahara M, et al. Molecular typing of Bartonella henselae DNA extracted from human clinical specimens and cat isolates in Japan. FEMS Immunol Med Microbiol. 2010;60:44–8. https://doi.org/10.1111/j.1574-695X.2010.00711.x.

    Article  CAS  PubMed  Google Scholar 

  40. Kosoy M, Mckee C, Albayrak L, Fofanov Y. Genotyping of Bartonella bacteria and their animal hosts: current status and perspectives. Parasitology. 2018;145:543–62. https://doi.org/10.1017/S0031182017001263.

    Article  CAS  PubMed  Google Scholar 

  41. Viezens J, Arvand M. Simultaneous presence of two different copies of the 16S rRNA gene in Bartonella henselae. Microbiology. 2008;154:2881–6. https://doi.org/10.1099/mic.0.2008/018630-0.

    Article  CAS  PubMed  Google Scholar 

  42. Iredell J, Blanckenberg D, Arvand M, Grauling S, Feil EJ, Birtles RJ. Characterization of the natural population of Bartonella henselae by multilocus sequence typing. J Clin Microbiol. 2003;41:5071–9. https://doi.org/10.1128/JCM.41.11.5071-5079.2003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Arvand M, Feil EJ, Giladi M, Boulouis HJ, Viezens J. Multi-locus sequence typing of Bartonella henselae isolates from three continents reveals hypervirulent and feline-associated clones. PLoS One. 2007;2:e1346. https://doi.org/10.1371/journal.pone.0001346.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Chaloner GL, Harrison TG, Coyne KP, Aanensen DM, Birtles RJ. Multilocus sequence typing of Bartonella henselae in the United Kingdom indicates that only a few, uncommon sequence types are associated with zoonotic disease. J Clin Microbiol. 2011;49:2132–7. https://doi.org/10.1128/JCM.00275-11.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Furquim MEC, Do Amaral R, Dias CM, Gonçalves LR, Perles L, Lima CAP, et al. Genetic diversity and multilocus sequence typing analysis of Bartonella henselae in domestic cats from Southeastern Brazil. Acta Trop. 2021;222:106037. https://doi.org/10.1016/j.actatropica.2021.106037.

    Article  CAS  PubMed  Google Scholar 

  46. Mietze A, Morick D, Kohler H, Harrus S, Dehio C, Nolte I, et al. Combined MLST and AFLP typing of Bartonella henselae isolated from cats reveals new sequence types and suggests clonal evolution. Vet Microbiol. 2011;148:238–45. https://doi.org/10.1016/j.vetmic.2010.08.012.

    Article  CAS  PubMed  Google Scholar 

  47. Joseph AK, Wood CW, Robson JM, Paul SL, Morris AJ. Bartonella henselae bacteraemia in domestic cats from Auckland. N Z Vet J. 1997;45:185–7. https://doi.org/10.1080/00480169.1997.36023.

    Article  CAS  PubMed  Google Scholar 

  48. Chomel BB, Boulouis HJ, Petersen H, Kasten RW, Yamamoto K, Chang CC, et al. Prevalence of Bartonella infection in domestic cats in Denmark. Vet Res. 2002;33:205–13. https://doi.org/10.1051/vetres:2002008.

    Article  PubMed  Google Scholar 

  49. Siğirci BD, Ilgaz A. Detection of the presence of Bartonella henselae in cats in Istanbul. J Fac Vet Med Istanbul Univ. 2013;39:209–17.

    Google Scholar 

  50. Heller R, Artois M, Xemar V, De Briel D, Gehin H, Jaulhac B, et al. Prevalence of Bartonella henselae and Bartonella clarridgeiae in stray cats. J Clin Microbiol. 1997;35:1327–31. https://doi.org/10.1128/jcm.35.6.1327-1331.1997.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Lindroos H, Vinnere O, Mira A, Repsilber D, Näslund K, Andersson SG. Genome rearrangements, deletions, and amplifications in the natural population of Bartonella henselae. J Bacteriol. 2006;188:7426–39. https://doi.org/10.1128/JB.00472-06.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Zhao F, Song XP, Li DM, Huang RT, Li ZF, Liu QY. Multilocus sequence typing analysis for Bartonella henselae isolates in China. Chin J Zoonoses. 2011;27:592–6.

    CAS  Google Scholar 

  53. Yuan C, Zhu C, Wu Y, Pan X, Hua X. Bacteriological and molecular identification of Bartonella species in cats from different regions of China. PLoS Neglect Trop D. 2011;5:e1301. https://doi.org/10.1371/journal.pntd.0001301.

    Article  Google Scholar 

  54. Dias CM, do Amaral RB, Perles L, Muniz ALDS, Rocha TFG, Machado RZ, et al. Multi-locus sequencing typing of Bartonella henselae isolates reveals coinfection with different variants in domestic cats from Midwestern Brazil. Acta Trop. 2023;237:106742. https://doi.org/10.1016/j.actatropica.2022.106742.

    Article  CAS  PubMed  Google Scholar 

  55. Ridder-Schröter R, Marx A, Beer M, Tappe D, Kreth HW, Girschick HJ. Abscess-forming lymphadenopathy and osteomyelitis in children with Bartonella henselae infection. J Med Microbiol. 2008;57:519–24. https://doi.org/10.1099/jmm.0.47438-0.

    Article  PubMed  Google Scholar 

  56. Namekata DY, Kasten RW, Boman DA, Straub MH, Siperstein-Cook L, Couvelaire K, et al. Oral shedding of Bartonella in cats: correlation with bacteremia and seropositivity. Vet Microbiol. 2010;146:371–5. https://doi.org/10.1016/j.vetmic.2010.05.034.

    Article  PubMed  Google Scholar 

  57. Fleischman DA, Chomel BB, Kasten RW, Stuckey MJ, Scarlet J, Liu H, et al. Bartonella infection among cats adopted from a San Francisco shelter, revisited. Appl Environ Microb. 2015;81:6446–50. https://doi.org/10.1128/AEM.01864-15.

    Article  CAS  Google Scholar 

  58. Lappin MR, Davis WL, Hawley JR, Brewer M, Morris A, Stanneck D. A flea and tick collar containing 10% imidacloprid and 4.5% flumethrin prevents flea transmission of Bartonella henselae in cats. Parasites Vectors. 2013;6:26. https://doi.org/10.1186/1756-3305-6-26.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Greco G, Brianti E, Buonavoglia C, Carelli G, Pollmeier M, Schunack B, et al. Effectiveness of a 10% imidacloprid/4.5% flumethrin polymer matrix collar in reducing the risk of Bartonella spp. infection in privately owned cats. Parasit Vectors. 2019;12:69. https://doi.org/10.1186/s13071-018-3257-y.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Raimundo JM, Guimarães A, Amaro GM, Da Silva AT, Botelho CFM, Massard CL, et al. Molecular survey of Bartonella species in shelter cats in Rio De Janeiro: clinical, hematological, and risk factors. Am J Trop Med Hyg. 2019;100:1321–7. https://doi.org/10.4269/ajtmh.18-0585.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Srisanyong W, Takhampunya R, Boonmars T, Kerdsin A, Suksawat F. Prevalence of Bartonella henselae, Bartonella clarridgeiae, and Bartonella vinsonii subsp berkhoffii in pet cats from four provincial communities in Thailand. Thai J Vet Med. 2016;46:663–70.

    Article  Google Scholar 

  62. Arvand M, Klose AJ, Schwartz-Porsche D, Hahn H, Wendt C. Genetic variability and prevalence of Bartonella henselae in cats in Berlin, Germany, and analysis of its genetic relatedness to a strain from Berlin that is pathogenic for humans. J Clin Microbiol. 2001;39:743–6. https://doi.org/10.1128/JCM.39.2.743-746.2001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Beugnet F, Labuschagne M, Fourie J, Jacques G, Farkas R, Cozma V, et al. Occurrence of Dipylidium caninum in fleas from client-owned cats and dogs in Europe using a new PCR detection assay. Vet Parasitol. 2014;205:300–6. https://doi.org/10.1016/j.vetpar.2014.06.008.

    Article  CAS  PubMed  Google Scholar 

  64. Diakou A, Di Cesare A, Accettura PM, Barros L, Iorio R, Paoletti B, et al. Intestinal parasites and vector-borne pathogens in stray and free-roaming cats living in continental and insular Greece. PLoS Neglect Trop Dis. 2017;11:e0005335. https://doi.org/10.1371/journal.pntd.0005335.

    Article  Google Scholar 

  65. Hassan AA, Akineden O, Usleber E. Identification of Streptococcus canis isolated from milk of dairy cows with subclinical mastitis. J Clin Microbiol. 2005;43:1234–8. https://doi.org/10.1128/JCM.43.3.1234-1238.2005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Cvetnić L, Benić M, Habrun B, Kompes G, Stepanić M. Samardžija M Najčešći uzročnici mastitisa u krava i koza u Republici Hrvatskoj/Most common causes of mastitis in cows and goats in Republic of Croatia. Vet Stanica. 2016;47:109–16.

    Google Scholar 

  67. Stojević D. Biokemijska i molekularna analiza izolata bakterije Escherichia coli izdvojenih iz hrane životinjskog podrijetla i obrisaka klaoničkih trupova [Biochemical and molecular analysis of Escherichia coli strains isolated from food of animal origin and carcass swabs]. Dissertation. Faculty of Veterinary Medicine, University of Zagreb, Zagreb; 2017.

  68. Lukačević D, Vidić Z, Katić S, Listeš E, Galić J. Štalska epizootija giardioze goveda u Dalmaciji. [Epizootic giardiasis in barn-housed cattle in Dalmatia]. Hrvatski Vet Vjesnik. 2016;24:74–9.

    Google Scholar 

  69. Gil H, García-Esteban C, Barandika JF, Peig J, Toledo A, Escudero R, et al. Variability of Bartonella genotypes among small mammals in Spain. Appl Environ Microbiol. 2010;76:8062–70. https://doi.org/10.1128/AEM.01963-10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Król N, Militzer N, Stöbe E, Nijhof AM, Pfeffer M, Kempf VAJ, et al. Evaluating transmission paths for three different Bartonella spp. in Ixodes ricinus ticks using artificial feeding. Microorganisms. 2021;9:901. https://doi.org/10.3390/microorganisms9050901.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Farrell S, McGarry J, Noble PM, Pinchbeck GJ, Cantwell S, Radford AD, et al. Seasonality and other risk factors for fleas infestations in domestic dogs and cats. Med Vet Entomol. 2023;37:359–70. https://doi.org/10.1111/mve.12636.

    Article  CAS  PubMed  Google Scholar 

  72. Manvell C, Ferris K, Maggi R, Breitschwerdt EB, Lashnits E. Prevalence of vector-borne pathogens in reproductive and non-reproductive tissue samples from free-roaming domestic cats in the South Atlantic USA. Pathogens. 2021;10:1221. https://doi.org/10.3390/pathogens10091221.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Siewert LK, Korotaev A, Sedzicki J, Fromm K, Pinschewer DD, Dehio C. Identification of the Bartonella autotransporter CFA as a protective antigen and hypervariable target of neutralizing antibodies in mice. Proc Natl Acad Sci USA. 2022;119:e2202059119. https://doi.org/10.1073/pnas.2202059119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Rodríguez-Pastor R, Hasik AZ, Knossow N, Bar-Shira E, Shahar N, Gutiérrez R, et al. Bartonella infections are prevalent in rodents despite efficient immune responses. Parasit Vectors. 2023;16:315. https://doi.org/10.1186/s13071-023-05918-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to the veterinarians and veterinary clinic employees from the following veterinary practices who voluntarily participated in the study: Anubis Klub (Vukovar), Bobanović (Dubrovnik), Cestica (Cestica), Hajster (Pula), Kućni ljubimci /now Ferarica/ (Rijeka), Pilić (Split), Veterinarska stanica Bjelovar (Bjelovar), Veterinarska stanica Jastrebarsko (Jastrebarsko), Veterinarska stanica Osijek (Osijek), Veterinarska ambulanta Pula (Pula), Veterinarska stanica Velika Gorica (Velika Gorica) and Veterinarska stanica Zabok (Zabok); also clinics in Zagreb: Buba, Veterina Branimir, Dodo-vet, Fabela, Goldi, Hop (now Vetvision), Ljubimac, Šegota. We also thank the cat owners and guardians for providing the samples.

Funding

This research was funded by GENOTICKTRECK of the Croatian Science Foundation (project no. 1957) and by the Croatian Veterinary Institute.

Author information

Authors and Affiliations

Authors

Contributions

RB and MS designed the study and drafted the manuscript. SH and VK substantively revised the original version. MS and SŠ cultured B. henselae. SD and IR performed MLST and phylogeny analyses. ŽM was responsible for the statistical analysis. All authors contributed to data acquisition, analysis and interpretation. All authors proofread and approved manuscript to be published.

Corresponding authors

Correspondence to Maja Stepanić or Relja Beck.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the ethics committee of the Faculty of Veterinary Medicine, University of Zagreb (640-01/16-17/63, 251-61-01/139-16-2). Consent to participate was provided by each cat owner or guardian, within the context of the doctoral thesis of the first author, Maja Stepanić.

Consent for publication

Owners and guardians of cats in this study gave written informed consent for publication of their data.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Stepanić, M., Duvnjak, S., Reil, I. et al. Epidemiology of Bartonella henselae infection in pet and stray cats in Croatia with risk factors analysis. Parasites Vectors 17, 48 (2024). https://doi.org/10.1186/s13071-024-06117-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-024-06117-8

Keywords