Skip to main content
Pathogens and Global Health logoLink to Pathogens and Global Health
. 2021 Nov 17;116(4):220–235. doi: 10.1080/20477724.2021.1989185

Global meta-analysis on Babesia infections in human population: prevalence, distribution and species diversity

Solomon Ngutor Karshima a,, Magdalene Nguvan Karshima b, Musa Isiyaku Ahmed c
PMCID: PMC9132453  PMID: 34788196

ABSTRACT

Human babesiosis is an emerging tick-borne protozoan zoonosis caused by parasites of the genus Babesia and transmitted by ixodid ticks. It was thought to be a public health problem mainly for the immunocompromised, however the increasing numbers of documented cases among immunocompetent individuals is a call for concern. In this systematic review and meta-analysis, we reported from 22 countries and 69 studies, an overall pooled estimate (PE) of 2.23% (95% CI: 1.46–3.39) for Babesia infections in humans. PEs for all sub-groups varied significantly (p < 0.05) with a continental range of 1.54% (95% CI: 0.89–2.65) in North America to 4.17% (95% CI: 2.11–8.06) in Europe. PEs for country income levels, methods of diagnosis, study period, sample sizes, Babesia species and targeted population ranged between 0.43% (95% CI: 0.41–0.44) and 7.41% (95% CI: 0.53–54.48). Babesia microti recorded the widest geographic distribution and was the predominant specie reported in North America while B. divergens was predominantly reported in Europe. Eight Babesia species; B. bigemina, B. bovis, B. crassa-like, B. divergens, B. duncani, B. microti, B. odocoilei and B. venatorum were reported in humans from different parts of the world with the highest prevalence in Europe, lower middle income countries and among individuals with history of tick bite and other tick-borne diseases. To control the increasing trend of this emerging public health threat, tick control in human settlements, the use of protective clothing by occupationally exposed people and the screening of transfusion blood in endemic countries are recommended.

Abbreviations AJOL: African Journals OnLine, CI: Confidence interval, CIL: Country income level, df: Degree of freedom, HIC: Higher-income countries, HQ: High quality, I2: Inverse variance index, IFAT: Indirect fluorescent antibody test, ITBTBD: Individuals with tick-bite and tick-borne diseases, JBI: Joanna Briggs Institute, LIC: Lower-income countries, LMIC: Lower middle-income countries, MQ: Medium quality, NA: Not applicable, N/America: North America, OEI: Occupational exposed individuals, OR: Odds ratio, PE: Pooled estimates, PCR: Polymerase chain reaction, Prev: Prevalence, PRISMA: Preferred Reporting System for Systematic Reviews and Meta-Analyses, Q: Cochran’s heterogeneity statistic, QA: Quality assessment, Q-p: Cochran’s p-value, qPCR: Quantitative polymerase chain reaction, S/America: South America, Seq: Sequencing, UMIC: Upper middle-income countries, USA: United States of America

KEYWORDS: Geographic distribution, human babesiosis, prevalence, species diversity, zoonotic Babesia species

Introduction

Zoonotic babesiosis is an emerging tick-borne protozoan disease caused by intraerythrocytic parasites of the genus Babesia, family Babesiidae, order Piroplasmorida and phylum Apicomplexa. The disease is predominantly caused by three parasites namely; Babesia divergens, B. microti and B. venatorum and transmitted by tick vectors. Although these zoonotic Babesia species are capable of infecting a wide range of vertebrate hosts across the globe, they require competent vertebrates and arthropod vectors; particularly ticks to complete their life cycle [1]. The sexual stage of the life cycle occurs in the tick vectors while the asexual stage occurs within the erythrocytes of the vertebrate host [2,3].

Babesia species have distinct geographic distributions. Transmission primarily occurs through bites of ticks of the family Ixodidae particularly those of the genus Ixodes. However, transfusion-associated [4,5] and vertical transmission from mothers to infants particularly for B. microti [6,7] are well documented. Quarter to half of infections are usually asymptomatic; usually confused with malaria due to nonspecific symptoms like fever and hemolytic anemia [8]. In immunocompromised population, these parasites infect host’s erythrocytes causing hemolytic anemia, intravascular coagulopathy, hepatomegaly and splenomegaly; leading to complications such as respiratory distress syndrome, heart failure, inflammation of the central nervous system and death [8].

Babesia divergens was the cause of the first human case of babesiosis in Croatia [9]. It is the most common cause of zoonotic babesiosis in Europe [10] and to a lesser extent in North America where the tick vectors are distributed [11]. The parasite is transmitted principally by Ixodes ricinus and cattle are the natural host [12]. Clinical symptoms due to B. divergens in humans appear within 1 to 3 weeks post infected tick bite and may include hemoglobinuria, jaundice resulting from severe hemolysis as well as renal failure and pulmonary edema in severe cases [8].

Babesia microti was first reported in 1966 [13] and is the principal cause of zoonotic babesiosis in the North America [14], however, an expanding geographic range to all New England states, north into Maine, west into the upper Hudson River valley, eastern Pennsylvania and New Jersey and as far south as Maryland has also been reported [15–18]. The principal vector of B. microti in the North America is Ixodes scapularis [19] and rodents are the reservoirs of infection [12]. Clinically, B. microti causes asymptomatic to severe disease with multi-organ failures and up to 20% fatality rates among the elderly, immunocompromised and asplenic patients [20,21].

Babesia venatorum is associated with zoonotic babesiosis in Asia and Europe [22] and is transmitted by the vectors Ixodes persulcatus in China [23] and Ixodes ricinus in Austria and Italy [24]. The roe deer serve as the primary vertebrate host of the parasite [22,25]. The first three reported cases of Babesia venatorum (formally known as Babesia sp. EU1) were in splenectomized occupational hunters from Austria, Italy [24] and Germany [26] who were above 50 years of age.

It is speculated that zoonotic babesiosis has suffered neglect and underestimation of burden as a result of misdiagnosis especially in resource limited countries with limited diagnostic capabilities. Though this disease is asymptomatic in immunocompetent population, the risk of fatality associated with old age [27], splenectomy [28], HIV infection [29] and malignancy [28] suggest the need for more attention to this emerging public health threat. Furthermore, the possibility of transmission to non-endemic countries as a result of asymptomatic blood donors may also be an important epidemiologic risk. For a better understanding of the worldwide burden of zoonotic babesiosis in humans, we undertook a systematic review and meta-analysis to determine its global prevalence, distribution and Babesia species diversity in humans.

Materials and methods

Study protocol and eligibility criteria

This systematic review and meta-analysis was performed and documented in accordance with the recommendations published by Moher et al. [30]. Studies were incorporated based on the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) checklist available from https://journals.plos.org/plosmedicine/article?id=10.1371/journal.pmed.1000097. The review protocol was registered on PROSPERO international prospective register of systematic reviews with registration number CRD42020153018 and available from https://www.crd.york.ac.uk/prospero/display_record.php?ID=CRD42020153018.

Data sources and search strategy

We systematically searched African Journals OnLine (AJOL), Google Scholars, PubMed, Science Direct, Scopus and Web of Science for a period of 18 months (1 March 2020–31 August 2021) for articles published on human Babesia infections worldwide from inception to 31 August 2021. References of retrieved studies were manually searched for additional articles. We requested full texts of relevant articles with only visible online titles and abstracts through texting and e-mailing of authors and editors of the publishing journals. We employed the MeSH search string ‘prevalence’ OR ‘sero-prevalence’ OR ‘epidemiology’ OR ‘sero-epidemiology’ OR ‘detection’ OR ‘molecular detection’ AND ‘zoonotic babesiosis’ OR ‘human babesiosis’ OR ‘human babesiasis’ OR ‘zoonotic Babesia species’ OR ‘Babesia (B.) divergens’ OR ‘B. microti’ OR ‘B. venatorum’ OR ‘B. duncani’ OR ‘B. crassa’ AND ‘humans’ OR ‘man’ OR ‘people’.

Study selection

All the authors (SNK, MNK and MIA) screened the title and abstract of each article and removed irrelevant and duplicate articles. Apparently relevant articles were further subjected to full text review for the identification and extraction of relevant information. A study was considered eligible only if it fulfill the following conditions: (i) it investigated human infections with Babesia species, (ii) was published in English, (iii) it stated numbers of positive cases and sample size, (iv) it stated study location, (v) was published until 31 August 2021, (vi) it stated the method of diagnosis used and (vii) it identified Babesia parasites to species level.

Data extraction

Data extraction was conducted by all authors (SNK, MNK and MIA) independently using a database developed in Microsoft Excel (Supplementary file 1). In cases of discrepancies, all authors double checked data simultaneously and discuss issues until unanimity was reached. From each of the articles, we extracted data including author’s name, year of conduct and publication, number of positive cases, sample sizes, study location, study design, method of diagnosis, Babesia species associated with the infection and where stated, the targeted population.

For the purpose of our analysis, studies which used diagnostic methods targeting antibodies or antigens were grouped as studies diagnosed by serology, while those that utilized polymerase chain reaction (PCR) either alone or in combination with sequencing were grouped as diagnosed by molecular methods. Individuals with tick bite history and/or tick-borne diseases (ITBTBD) referred to people that had history of tick bite and those who tested positive for human granulocytic anaplasmosis and/or Lyme borreliosis as well as those with tick-borne disease symptoms like erythema migrans. Occupational exposed group referred to people whose occupations put them at high risk of infection with tick-borne diseases like foresters, hunters, farmers, livestock owners and veterinarians. The category others included all normal individuals, hospitalized individuals without any specified disease and individuals with conditions other than tick-borne diseases targeted by the individual studies. More so, zoonotic Babesia species refer to species of Babesia like B. divergens, B. microti and B. venatorum which are infective to man and animals under natural conditions. Analysis for country-based pooled estimates were only performed for countries that reported at least two studies. We grouped countries income levels based on the World Bank year 2020 classifications [31].

Quality assessment for individual studies

All authors (SNK, MNK and MIA) independently screened articles explicitly for quality based on the checklist of the Joanna Briggs Institute (JBI) critical appraisal instrument for studies reporting prevalence data published by Munn et al. [32]. This was carried out on a database developed in Microsoft Excel (Supplementary file 2), and discrepancies in data collected were simultaneously doubled checked by all authors and discussed until consensus was reached. Each article was examined for appropriateness of sample frame and the way the study participants were sampled, completeness of the description of study subjects and settings as well as the sufficiency of data analysis of the identified sample. The articles were also examined for, validity of the methods employed for the detection of Babesia infections, reliability of the measurement of the condition in all participants, appropriateness of the statistical analysis used and adequacy of the response rate and its management. We assigned scores of 0 or 1 for no or yes answers and NA when the item was not applicable for any given article. We slightly modified the grading by grouping studies based on total scores earned as follows; 0–3 (low quality), 4–6 (medium quality) and 7–9 (high quality).

Pooling and heterogeneity analysis

Data were managed in Microsoft Excel and transferred to Graph-Pad Prism version 4.0 and Comprehensive Meta-Analysis version 3.0 for statistical and meta-analysis respectively. Prevalence of each included study was determined by expressing the proportion of positive cases and sample size as a percentage. We used the random-effects model to determine pooled estimates of human Babesia infections and their 95% Confidence interval [33]. Furthermore, we used the Cochran’s Q-test to determine variations across studies and quantified the percentage variations among studies due to heterogeneity using the inverse variance index (I2); which is expressed as I2 = 100(Q-df)/Q; where Q represents Cochran heterogeneity statistic and df, the degree of freedom; which is equal to number of studies analyzed minus 1. I2 values of 0, 25, 50 and 75% represented absence of, low, moderate and substantive heterogeneities respectively [34,35].

Publication bias, sensitivity, sub-group and meta-regression analyses

Across study bias was investigated by the funnel plots, while its statistical significance was examined using Egger’s regression asymmetry test [36]. We also employed the Duval and Tweedie non-parametric ‘fill and trim’ linear random method to test for unbiased estimates [37]. We used one study deletion analysis at a time to test the effect of each eligible study on the pooled estimate of human Babesia infections, and a study was considered to have no influence if the pooled estimate without it was within the 95% confidence limits of the overall PE [38].

We performed sub-groupings for continents (Africa, Asia, Australia, Europe, North America and South America), country income levels (high-income countries; HIC, low-income countries; LIC, lower middle-income countries; LMIC, and upper middle-income countries; UMIC), methods of diagnosis (serology and molecular), study period (1975–1990, 1991–2005 and 2006–2021), sample sizes (≤500, 501–1000 and >1000), Babesia species (B. bigemina, B. bovis, B. crassa-like, B. divergens, B. duncani, B. microti, B. odocoilei, and B. venatorum) as well as targeted population (blood donors, individuals with tick-bite and/or tick-borne diseases; ITBTBD, occupationally exposed individuals; OEI and others).

To identify possible sources of heterogeneity in our analysis, we performed meta-regression for different strata by comparing South America with other continents (Africa, Asia, Australia, Europe and North America) where studies were reported, UMIC with other country income level groups (HIC, LIC and LMIC), serological methods with molecular methods, 2006–2021 with other study periods (1975–1990 and 1991–2005), sample sizes >1000 with other sample sizes (≤500 and 501–1000), Babesia venatorum with other species (B. bigemina, B. bovis, B. crassa-like, B. divergens, B. duncani, B. microti and B. odocoilei) and others with other target populations (blood donors, ITBTBD and OEI).

Results

Study selection

Overall, our literature search revealed a total of 1027 studies (1003 form databases and 24 from manual search of references). A total of 792 duplicates were excluded after title screening and another 142 animal-based studies were excluded after scrutinizing abstracts. We finally subjected 93 studies to full text review and 24 studies were excluded for the following reasons; case reports (n = 16) and insufficient data on prevalence (n = 8). We included in this quantitative synthesis a total of 69 studies (Figure 1).

Figure 1.

Figure 1.

Flow chart for the selection of eligible studies

Study characteristics, quality assessment and distribution of Babesia species

Table 1 shows the baseline characteristics of eligible studies and is ordered by continent/country. Studies were reported from six continents: Africa and Australia reported one study each, South America (2/69) Asia (9/69), Europe (22/69) and North America (34/69). Individual countries with the highest concentration of studies across continents were China in Asia (7/69; 10.14%), Poland in Europe (10/69; 14.49%) and the USA in North America (33/69; 47.83%). The 2 studies in South America were from Bolivia (n = 1) and Colombia (n = 1), while the only study in Africa was from Tanzania. Prevalence of zoonotic Babesia in individual studies analyzed ranged between 0.00 and 81.91%. Studies were reported from LIC (1/69), LMIC (3/69), UMIC (8/69) and HIC (57/69). Forty studies diagnosed babesiosis using serological methods, while the remaining 29 utilized molecular techniques Figure 2. The majority of studies were conducted during 2006–2021 (49/69; 71.01%) with the remainder during 1991–2005 (14/69; 20.29%) and 1975–1990 (6/69; 8.70%). Furthermore, 37 studies had sample sizes of ≤500, 9 had sample sizes of 501–1000 and 23 had sample sizes of >1000. Babesia bigemina, B. bovis, B. crassa-like and B. odocoilei were reported by one study each, B. duncani (2 studies), B. venatorum (3 studies), B. divergens (11 studies) and B. microti (57 studies). Overall, B. microti had the widest geographic distribution globally, while continental distribution revealed predominance of B. divergens and B. microti in Europe and North America respectively. In addition, 2 of the 3 reports on B. venatorum infections in humans were from Asia. The population targeted by the 69 eligible studies were blood donors (n = 23), ITBTBD (n = 26), OEI (n = 11) and other individuals (n = 11) as shown in Table 1. No study was excluded for lack of merit based on assessment of quality. The majority of studies (62/69; 89.86%) were studies of medium quality (4–6 points) and the remainder (7/69; 10.14%) were of high quality (7–9 points) as shown in Table 1.

Table 1.

List and baseline characteristics of studies on Babesia infections in humans

Continent/Country CIL Study year Targeted population MOD Babesia species Sample size Cases Prev. (%) 95% CI QA Study reference
Africa                      
Tanzania LIC 2017 Normal children and adults IFAT B. microti 299 4 1.34 0.37–3.39 MQ [39] Bloch et al., 2019b
Asia                      
China UMIC 2011–14 Individuals with recent tick-bite Npcr B. venatorum 2912 48 1.65 1.22–2.18 MQ [123] Jiang et al., 2015
China UMIC 2016 Blood donors IFAT B. microti 1000 13 1.30 0.69–2.21 MQ [40] Bloch et al., 2018
China UMIC 2015/16 Individuals with recent tick-bite Mic & PCR/Seq B. crassa-like 1125 58 5.16 3.94–6.61 MQ [124] Jia et al., 2018
China UMIC 2010/11 Individuals bitten by ticks PCR/Seq B. venatorum 180 2 1.11 0.13–3.96 MQ [41] Liu et al., 2019
China UMIC 2009 Anemic patients PCR/Seq B. divergens 377 2 0.53 0.06–1.90 MQ [113] Qi et al., 2011
China UMIC 2018 Individuals with history of tick bite PCR/Seq B. divergens 754 10 1.33 0.64–2.43 HQ [114] Wang et al., 2019
China UMIC 2013 Humans with malaria-like symptoms PCR/Seq B. microti 449 10 2.23 1.07–4.06 MQ [112] Zhou et al., 2013
Japan HIC 1985 Individuals from a tick-borne disease endemic area IFAT B. microti 1335 18 1.35 0.80–2.12 MQ [42] Arai et al., 2003
Mongolia LMIC 2011 Stock farmers PCR B. microti 100 3 3.00 0.62–8.52 MQ [117] Hong et al., 2014
Australia                      
Australia HIC 2009–11 Self-referred patients for LD FISH, IFA & PCR All species 41 13 31.71 18.08–48.09 MQ [43] Mayne, 2011
          B. duncani 41 10 24.39 12.36–40.30    
          B. microti 41 7 17.07 7.15–32.06    
Europe                      
Austria HIC 2009 Blood donors IFAT All species 988 26 2.63 1.73–3.83 MQ [44] Sonnleitner et al., 2014
          B. divergens 988 21 2.13 1.32–3.23    
          B. microti 988 5 0.51 0.16–1.18    
Belgium HIC 2010 Humans with tick bite and tickborne disease-like symptoms IFAT All species 199 163 81.91 75.85–87.00 MQ [95] Lempereur et al., 2015
          B. divergens 199 66 33.17 26.67–40.17    
          B. microti 199 18 9.05 5.45–13.92    
          B. venatorum 199 79 39.70 32.85–46.86    
Croatia HIC 1999 Individuals with history of tick bite IFAT B. microti 102 1 0.98 0.02–5.34 MQ [96] Topolovec et al., 2003
Finland, Norway, Sweden HIC 2007 High risk tick bitten individuals IFAT B. microti 86 0 0.00 0.00–4.20 MQ [97] Ocias et al., 2020
Finland, Sweden HIC 2019 Humans bitten by Babesia species positive ticks IFAT B. microti 159 7 4.40 1.79–8.86 MQ [98] Wilhelmsson et al., 2020
France HIC 2003 Foresters IFAT All species 810 21 2.59 1.61–3.94 MQ [99] Rigaud et al., 2016
          B. divergens 810 1 0.12 0.00–0.69    
          B. microti 810 20 2.47 1.51–3.79    
Germany HIC 1999 Blood donors IFAT All species 467 42 8.99 6.56–11.96 MQ [45] Hunfeld et al., 2002
          B. divergens 467 17 3.64 2.13–5.76    
          B. microti 467 25 5.35 3.49–7.80    
Italy HIC 2009 Farmers, livestock keepers, veterinary practitioners and hunters IFAT All species 432 37 8.56 6.10–11.61 MQ [100] Gabrielli et al., 2014
          B. divergens 432 17 3.94 2.31–6.23    
          B. microti 432 20 4.63 2.85–7.06    
Netherlands HIC 2008 Individuals with tick bites history and erythema migrans PCR/Seq B. divergens 291 1 0.34 0.01–1.90 MQ [101] Jahfari et al., 2016
Norway HIC 2016 Healthy adults IFAT B. microti 1537 33 2.15 1.48–3.00 MQ [102] Thortveit et al., 2020
Poland HIC 2008 Farmers, foresters, blood donors IFAT B. microti 190 10 5.26 2.55–9.74 MQ [103] Chmielewska-Badora et al., 2012
Poland HIC 2012 Tick-borne encephalitis patients PCR/Seq B. microti 110 1 0.91 0.02–4.96 MQ [104] Moniuszko et al., 2014
Poland HIC 2013 Individual with history of tick bite PCR/Seq B. microti 548 6 1.09 0.40–2.37 MQ [105] Moniuszko-Malinowska et al., 2016
Poland HIC 2010 Foresters IFAT B. microti 114 5 4.39 1.4–-9.94 MQ [106] Pancewicz et al., 2011
Poland HIC 2013 Farmers and foresters IFAT B. microti 93 1 1.08 0.03–5.85 MQ [107] Panczuk et al., 2016
Poland HIC 2018 HIV-patients and blood donors IFAT B. microti 426 23 5.40 3.45–7.99 MQ [46] Pawelczyk et al., 2019
Poland HIC 2008 Lyme disease patients PCR/Seq B. divergens 24 1 4.17 0.11–21.12 MQ [108] Welc-Faleciak et al., 2010
Poland HIC 2011 Foresters PCR/Seq B. microti 58 2 3.45 0.42–11.91 MQ [109] Welc-Faleciak et al., 2015
Poland HIC 2013 Foresters IFAT B. microti 216 50 23.15 17.70–29.35 MQ [110] Zukiewicz-Sobczak et al., 2014
Sweden HIC 2015 Lyme disease patients IFAT All species 283 19 6.71 4.09–10.29 HQ [47] Svensson et al., 2019
          B. divergens 283 8 2.83 1.23–5.49    
          B. microti 283 11 3.89 2.96–6.85    
Switzerland HIC 1998 Blood donors IFAT B. microti 396 6 1.52 0.56–3.27 MQ [48] Foppa et al., 2002
Ukraine LMIC 2020 Blood donors as well as HIV and LD patients IFA All species 145 15 10.34 5.91–16.49 MQ [49] Bondarenko et al., 2021
          B. divergens 145 10 6.90 3.36–12.32    
          B. microti 145 5 3.45 1.13–7.86    
N/America                      
Canada HIC 2018 Blood donors TMA B. microti 50,752 1 0.00 0.00, 0.00 MQ [50] Tonnetti et al., 2019a
USA HIC 1999–07 Patient records Mic/PCR Unspecified 326,081 353 0.11 0.10–0.12 MQ [51] Asad et al., 2009
USA HIC 2013 Blood donors PCR B. microti 26,702 134 0.50 0.42–0.59 MQ [52] Bloch et al., 2016
USA HIC 1985 Children on an Indian reservation in North Carolina IFA B. microti 185 6 3.24 1.20–6.92 MQ [53] Chisholm et al., 1986
USA HIC 2015 Lyme disease patients IFAT B. microti 130 35 26.92 19.52–35.40 MQ [54] Curcio et al., 2016
USA HIC 1978 Healthy individuals IFAT B. microti 136 6 4.41 1.64–9.36 MQ [111] Filstein et al., 1980
USA HIC 1995 High risk adult population IFAT B. microti 671 7 1.04 0.42–2.14 MQ [55] Hilton et al., 1999
USA HIC 2007 Blood donors IFAT B. microti 17,465 267 1.53 1.35–1.72 MQ [56] Johnson et al., 2009
USA HIC 2005 Blood donors PCR B. microti 17,422 208 1.19 1.04–1.37 MQ [57] Johnson et al., 2012
USA HIC 2009 Blood donors PCR B. microti 1002 3 0.30 0.06–0.87 MQ [118] Johnson et al., 2013
USA HIC 2004–13 Patient records Mic & PCR B. microti 3138 22 0.70 0.44–1.06 MQ [58] Kowalski et al., 2015
USA HIC 1989 Lyme disease patients, randomly selected outpatients, college students and healthy individuals IFAT B. microti 1285 81 6.30 5.04–7.77 MQ [59] Krause et al., 1991
USA HIC 1991 Children and adults IFAT B. microti 574 52 9.06 6.84–11.71 HQ [60] Krause et al., 1992
USA HIC 1990–94 Long-term residents of island community and LD patients Serology Unspecified 1396 40 2.87 2.05–3.88 MQ [61] Krause et al., 1996
USA HIC 1998 Blood donors IFAT B. microti 2006 9 0.45 0.21–0.85 MQ [62] Leiby et al., 2002
USA HIC 1999 Blood donors IFAT B. microti 3490 30 0.86 0.58–1.22 MQ [63] Leiby et al., 2005
USA HIC 2012 Blood donors IFAT B.microti 1272 14 1.10 0.60–1.84 HQ [64] Levin et al., 2014
USA HIC 2013 Blood donors EIA, IFA, Mic & PCR B. microti 26,703 56 0.21 0.16–0.27 MQ [65] Levin et al., 2016
USA HIC 1997 Healthy blood donors from a highly endemic area IFA/PCR B. microti 150 5 3.33 1.09–7.61 MQ [66] Linden et al., 2000
USA HIC 1990 Lyme disease patients and normal individuals IFAT B. microti 172 3 1.74 0.36–5.01 MQ [67] Magnarelli et al., 1995
USA HIC 1996 Human granulocytic ehrlichiosis patients IFAT B. microti 494 47 9.51 7.0–12.45 MQ [68] Magnerelli et al., 1998
USA HIC 1995 LD and HGE patients IFA B. microti 115 6 5.22 1.94–11.01 MQ [69] Mitchell et al., 1996
USA HIC 2014 Blood donors PCR B. microti 84,209 331 0.39 0.35–0.44 MQ [70] Moritz & Stramer, 2015
USA HIC 2011 Blood donors AFIA B. microti 13,269 44 0.33 0.24–0.45 HQ [71] Moritz et al., 2014
USA HIC 2014 Blood donors PCR B. microti 89,152 67 0.08 0.06–0.10 MQ [72] Moritz et al., 2016
USA HIC 1984 Blood donors IFAT B, microti 779 29 3.72 2.51–5.30 MQ [73] Popovsky et al., 1988
USA HIC 2017 Individual with history of tick bite Qpcr B. microti 192 92 47.92 40.67–55.23 MQ [74] Primus et al., 2018
USA HIC 2008/09 Blood donors and individuals requesting for B. duncani test IFA All species 900 22 2.44 1.54–3.68 MQ [75] Prince et al., 2010
          B. duncani 900 18 2.00 1.19–3.14    
          B. microti 900 4 0.44 0.12–1.13    
USA HIC 2000–16 Patients who received medical care within 2000–2016 Mic & PCR B. microti 2956 213 7.21 6.30–8.20 HQ [76] Rau et al., 2020
USA HIC 2020 Humans with babesiosis symptoms PCR/Seq B. odocoilei 19 2 10.53 1.30–33.14 MQ [77] Scott et al., 2021
USA HIC 2001 Patients with erythema migrans PCR B. microti 93 2 2.15 0.26–7.55 MQ [78] Steere et al., 2003
USA HIC 2010/11 Blood donors IFA/RT-PCR B. microti 2150 42 1.95 1.41–2.63 MQ [79] Tonnetti et al., 2013
USA HIC 2018 Blood donors AFIA B. microti 506,540 1299 2.96 2.43–2.71 MQ [80] Tonnetti et al., 2019b
USA HIC 2015 Lyme disease patients IFAT B.microti 52 4 7.69 2.14–18.54 MQ [81] Wormser et al., 2019
S/America                      
Bolivia LMIC 2013 Healthy rural residents PCR B. microti 271 9 3.32 1.53–6.21 MQ [82] Gabrielli et al., 2016
Colombia UMIC 2014/15 Herders ELISA, IFA & nPCR All species 300 6 2.00 0.74, 4.30 HQ [83] Gonzalez et al., 2018
          B. bigemina 300 2 0.67 0.08–2.39    
          B. bovis 300 4 1.33 0.36–3.38    

CI (Confidence interval), CIL (Country income level), HIC (Higher income countries), HQ (High quality), IFAT (Indirect fluorescent antibody test), LIC (Lower income countries), LMIC (Lower middle income countries), MOD (Method of diagnosis), MQ (Medium quality), N/America (North America), PCR (Polymerase chain reaction), Prev. (Prevalence), QA (Quality assessment), qPCR (Quantitative polymerase chain reaction), S/America (South America), Seq (Sequencing), UMIC (Upper middle income countries), USA (United States of America)

Figure 2.

Figure 2.

Global distribution of studies on human babesiosis

Pooled estimates, sub-group and heterogeneity analyses

Pooled estimates of human Babesia infections and heterogeneities among studies are presented in Tables 2 and 3 as well as Figures 3 and 4. Exactly 4198 of 1,198,469 individuals examined during the period under review were positive for Babesia infections, yielding an overall PE of 2.23% (95% CI: 1.46–3.39). PEs of human Babesia infections were 1.54% (95% CI: 0.89–2.65) in North America, 1.76% (95% CI: 1.07–2.87) in Asia, 2.71% (95% CI: 1.64–4.45) in South America and 4.17% (95% CI: 2.11–8.06) in Europe. Single studies reported prevalence of 1.34% (95% CI: 0.50–3.51) in Africa and 31.71% (95% CI: 19.39–47.27) in Australia. Based on country income level, the highest PE was in LMIC (5.03%; 95% CI: 2.04–11.89), followed by HIC (2.26%; 95% CI: 1.41–3.60) and then UMIC (1.76%; 95% CI: 1.03–2.98). A single study from LIC reported a prevalence of 1.34% (95% CI: 0.50–3.51). PE by molecular techniques was 1.30% (95% CI: 0.65–2.57), while that by serological methods was 3.24% (95% CI: 1.78–5.81). The prevalence of human Babesia infections decreased by 2.74% during the 45 years under review. The highest occurred during 1975–1990 (3.17%; 95% CI: 1.79–5.56), followed by 1991–2005 (2.44%; 95% CI: 1.56–3.18) and then 2006–2021 (0.43%; 95% CI: 0.42–0.44).

Table 2.

Sub-group analysis for the pooled estimates of Babesia infections in humans

Variables No. of
Studies
Pooled Estimates
(95% CI)  
Heterogeneity
Meta-regression
Sample size Cases Prev. (%) P-value Q-value I2 (%) Q-p P- value OR (95% CI)
Continent                      
Africa 1 299 4 1.34 0.50–3.51 <0.001 0.00 0.00 1.000 0.1091 −0.67 (−4.62, 3.28)
Asia 9 8232 164 1.76 1.07–2.87   62.27 87.15 <0.001   −0.45 (−2.94, 2.04)
Australia 1 41 13 31.71 19.39–47.27   0.00 0.00 1.000   2.86 (−1.02, 6.74)
Europe 22 7674 470 4.17 2.11–8.06   729.37 97.12 <0.001   0.49 (−1.86, 2.85)
North America 34 1,181,652 3532 1.54 0.89–2.65   7310.19 99.55 <0.001   −0.53 (−2.84, 1.79)
South America 2 571 15 2.71 1.64–4.45   0.95 0.00 0.329   Reference
Income level                      
HIC 57 1,190,557 4018 2.26 1.41–3.60 <0.001 10,280.30 99.46 <0.001 0.8216 0.34 (−1.00, 1.68)
LIC 1 299 4 1.34 0.50–3.51   0.00 0.00 1.000   0.09 (−3.89, 4.06)
LMIC 3 516 27 5.03 2.04–11.89   9.57 79.11 0.008   1.12 (−1.29, 3.52)
UMIC 8 7097 149 1.76 1.03–2.98   56.09 87.52 <0.001   Reference
MOD                      
Molecular 29 636,111 1661 1.30 0.65–2.57 <0.001 4661.85 99.40 <0.001 0.0491 0.94 (0.00, 1.87)
Serology 40 562,358 2537 3.24 1.78–5.81   5974.95 99.35 <0.001   Reference
Study period                      
1975–1990 6 3892 143 3.17 1.79–5.56 <0.001 41.18 87.86 <0.001 0.8595 0.39 (−1.11, 1.89)
1991–2005 14 28,186 476 2.44 1.56–3.18   458.52 96.95 <0.001   0.15 (−0.91, 1.21)
2006–2021 49 1,166,391 3579 0.43 0.42–0.44   9845.20 99.50 <0.001   Reference
Sample size                      
≤500 37 7546 641 4.90 3.00–7.91 <0.001 905.48 96.02 <0.001 <0.001 1.97 (1.25, 2.70)
501–1000 9 7024 186 2.28 1.35–3.82   96.67 91.72 <0.001   1.14 (0.08, 2.19)
>1000 23 1,183,899 3371 0.73 0.42–1.26   4926.18 99.55 <0.001   Reference
Overall 69 1,198,469 4198 2.23 1.46–3.39   10,678.92 99.36 <0.001    

CI (Confidence interval), HIC (High income countries), I2 (Inverse variance index), IFAT (Indirect fluorescent antibody technique), LIC (Low income countries), LMIC (Lower middle income countries), OR (Odds ratio), PCR (Polymerase chain reaction), Prev. (Prevalence), Q-p (Cochran’s p-value), UMIC (Upper middle income countries)

Table 3.

Pooled estimates of human Babesia infections in relation to target population and species

Variables No. of
Studies
Pooled Estimates
95% CI  
Heterogeneity
Meta-regression
Sample size Cases Prev. (%) P-value Q-value I2 (%) Q-p P- value OR (95% CI)
Target population                    
Blood donors 23 847,300 2672 0.87 0.55–1.36 <0.001 2237.12 99.01 <0.001 <0.001 −1.03 (−2.13, 0.07)
ITBTBD 26 12,127 665 5.00 2.62–9.35   1166.74 97.86 <0.001   0.82 (−0.27, 1.91)
OEI 11 3169 149 3.80 1.88–7.53   147.29 93.21 <0.001   0.49 (−0.81, 1.78)
Others 11 335,848 712 2.31 0.51–9.75   2894.13 99.65 <0.001   Reference
Babesia species                      
B. bigemina 1 300 2 0.67 0.17–2.63 <0.001 0.00 0.00 1.000 0.7823 0.18 (−4.14, 4.50)
B. bovis 1 300 4 1.33 0.50–3.50   0.00 0.00 1.000   0.88 (−3.33, 5.08)
B. crassa-like 1 1125 58 5.16 4.01–6.61   0.00 0.00 1.00   2.27 (−1.83, 6.37)
B. divergens 11 4770 154 2.39 0.89–6.27   279.36 96.42 <0.001   1.46 (−1.13, 4.05)
B. duncani 2 941 28 7.41 0.53–54.48   40.33 97.52 <0.001   2.65 (−0.71, 6.02)
B. microti 57 864,895 3460 2.10 1.36–3.24   7638.00 99.27 <0.001   1.34 (−1.07, 3.75)
B. odocoilei 1 19 2 10.53 2.65–33.74   0.00 0.00 1.000   3.04 (−1.30, 7.39)
B. venatorum 3 3291 129 4.91 0.28–49.02   329.24 99.39 <0.001   2.26 (−0.83, 5.34)
Unspecified spp. 2 327,477 393 0.56 0.02–12.55   382.03 99.74 <0.001   Reference

CI (Confidence interval), I2 (Inverse variance index), ITBTBD (Individuals with tick-bite and tick-borne diseases), OEI (Occupational exposed individuals), OR (Odds ratio), Prev. (Prevalence), Q-p (Cochran’s p-value)

Figure 3.

Figure 3.

Country-based pooled estimates of Babesia infections in humans

Figure 4.

Figure 4.

Forest plot for the global pooled estimates of Babesia infections in humans

Prevalence of human Babesia infections decreased with increasing sample size, with the highest among studies with sample size ≤500 (4.90%; 95% CI: 3.00–7.91), followed by studies with sample size 501–1000 (2.28%; 95% CI: 1.35–3.82) and then those with sample size >1000 (0.73%; 95% CI: 0.42–1.26). Babesia infections were most prevalent among ITBTBD (5.00%; 95% CI: 2.62–9.35), with OEI in second place (3.80%; 95% CI: 1.88–7.53). PE among the others category was 2.31% (95% CI: 0.51–9.75), with the lowest observed among blood donors (0.87%; 95% CI: 0.55–1.36). Eight different Babesia species were reported in humans during the period under review. The highest PE was recorded by B. duncani (7.41%; 95% CI: 0.53–54.48), followed by B. venatorum (4.91%; 95% CI: 0.28–49.01), then B. divergens (2.39%; 95% CI: 0.89–6.27) and finally B. microti (2.10%; 95% CI: 1.36–3.24). Single studies reported the prevalence of 0.67% (95% CI: 0.17–2.63) for B. bigemina, 1.33% (95% CI: 0.50–3.50) for B. bovis, 5.16% (95% CI: 4.01–6.61) for B. crassa-like and 10.53% (95% CI: 2.65–33.74) for B. odocoilei. Country-based PEs ranged from 1.72% (95% CI: 0.95–3.08) in China to 4.95% (95% CI: 2.55–9.46) in Sweden. Overall heterogeneity was 99.36%, with a range of 0.00–99.74% across sub-groups.

Across study bias, sensitivity and meta-regression analyses

The funnel plot (Figure 5) and its respective bias coefficient for studies published worldwide on human Babesia infections (b = 7.84, 95% CI = 4.10–11.59, p = 0.00009) showed significant publication bias at α 0.05. We also observed that all the PEs following one study deletion analysis remained within the 95% confidence limits of the overall pooled estimate (Supplementary file 3). Results of meta-regression analysis (Tables 2 and 3) for sub-groups excluded study locations (Q = 9.00, df = 5, p = 0.1091), country income level (Q = 0.92, df = 3, p = 0.8216), study periods (Q = 0.30, df = 2, p = 0.8595) and Babesia species diversity (Q = 4.77, df = 8, p = 0.7823) as significant causes of heterogeneity. However, methods of diagnosis (Q = 3.87, df = 1, p = 0.0491), sample sizes of the individual studies (Q = 28.19, df = 2, p < 0.001) and targeted population (Q = 18.48, df = 3, p = 0.0004) were shown to contribute to the heterogeneity.

Figure 5.

Figure 5.

Funnel plot of standard error vs logit event rate for studies published on human Babesia infections worldwide

Discussion

Several surveillance studies reported zoonotic babesiosis in humans in different countries, however, no study harmonized these individual surveillance studies to provide the global situation of this emerging public health problem. We reported here, the global status of Babesia infections in humans in terms of prevalence, distribution and species diversity. To our knowledge, this is the first systematic review and meta-analysis to report the global prevalence of zoonotic babesiosis in humans. This study became pertinent following increased numbers of asymptomatic [84,85] and fatal [27,86–88] case reports of zoonotic Babesia from different countries in both immunocompromised [28,29] and immunocompetent [89–91] individuals as well as the risk of misdiagnosis and possible transmission by travelers to non-endemic countries [91,92].

The overall low global prevalence of 2.23% revealed by the present study may be attributable to factors including low trans-stadial transmission of zoonotic Babesia species which last for 2–3 years, very low parasitemia which can be easily missed by researchers and under-diagnosis in non-or-newly endemic areas [21,91,92]. Although the present study reported low global prevalence of these emerging pathogens, the fatality caused among immunocompromised population calls for concern, especially with the increasing numbers of people living with immunocompromised conditions.

Pooled estimates of Babesia infections in humans varied across geographic regions with the highest in Europe and the lowest in North America. These variations across regions may be attributed to the abundance of competent tick vectors [93] and reservoirs which are capable of maintaining these parasites [94]. The fact that over 76% of the studies published from Europe targeted occupationally exposed populations (foresters, farmers, livestock keepers, veterinarians) as well as individuals with history of tick bites and those with other tick-borne diseases [95–110] may be another possible reason for the higher prevalence in the region. Substantive evidence showed that the prevalence of human Babesia infections can vary even within a country. For instance, the prevalence of 1.1, 3.1 and 4.2% were reported from the North east [111], Eastern [112] and South eastern Poland [108] while in China, prevalence of 0.5, 1.3 and 2.2% were reported from the Shandong [113], Heilongjiang [114] and Yunnan [112] Provinces respectively.

We observed the highest prevalence of human Babesia infections in lower middle-income countries like Bolivia and Mongolia, which we attributed to limited resources for vector control and the poor living conditions of the lower class in these countries [115]. The high sensitivity (88–96%) and specificity (90–100%) of the immunofluorescence antibody technique [116] utilized by over 90% of the serology-based studies included in the analysis may be a possible explanation for the higher prevalence revealed by studies diagnosed using serology. The fact that serology detects both active and convalescent infections may be another factor. Similar findings were reported even at country levels like in Mongolia [117] and the United States of America [118].

The higher prevalence observed among individuals with history of tick bites and other tick-borne diseases suggests bites from infected ticks and possible transmission of multiple infections by ticks. Exposure of foresters, farmers, veterinarians and herders to tick bites at work places may explain the high prevalence among occupationally exposed population. The detection of zoonotic Babesia among blood donors suggests possible risk of transfusion-associated babesiosis which is a growing problem in areas endemic for zoonotic babesiosis [5,119].

The majority of studies we included in our analysis were published between 2006 and 2021, probably due to increased understanding of the public health significance of the pathogen. Overall, the prevalence of human Babesia infections decreased by 2.74% during the period of 45 years under review. Further analysis revealed a 0.73% initial decline between 1990 and 2005, which was followed by a 2.01% decline between 2006 and 2021. This consistent decline is suggestive of possible successes of existing control programmes targeting human babesiosis, especially with the increasing number of prevalence studies and improvement in current diagnostic techniques. Babesia microti showed the widest geographic distribution globally and was the predominant species reported across North America, while B. divergens was predominantly reported in Europe, concurring with reports from existing literature [10,120–122]. On the other hand, two of the three reports on Babesia venatorum were from Asia [123,124] and one from Belgium [95].

We found significant publication bias by the funnel plot. This publication bias might be due to gap in literature arising from non-publication of results that contradict existing knowledge as well as dependence of research publication on the nature and direction of study results [125] and its statistical significance [126]. From the results of one study deletion analysis performed, all the pooled estimates were within the 95% CI of the overall PE, suggesting that no single study influenced the PE in our analysis. We observed substantive heterogeneity across studies which persisted even after sub-group analysis. Meta-regression for sub-groups showed that the sample sizes of the eligible studies, the methods of diagnosis as well as the population targeted by the individual studies were the possible sources of heterogeneity in the present analysis.

The present study has a number of implications on public health. First, the detection of zoonotic Babesia among blood donors in different countries may influence transfusion-associated babesiosis. Second, for travelers from endemic to non-endemic areas, there is the risk of misdiagnosis by physicians and possible transmission by available ticks to susceptible individuals. Third, asymptomatic individuals may serve as reservoirs of infections thereby initiating transmission through activities of peri-domestic ticks.

A number of limitations are associated with our study. We could not include a number of studies which would have increased our understanding of the global status of human babesiosis due to insufficient prevalence data. Again, over 80% (56/69) of the studies were concentrated in two continents namely; Europe (22) and North America (34) giving an uneven distribution of studies. Studies on human Babesia infections were available from only 22 of 249 countries of the world. In addition, over 57% (40/69) of the studies included in this analysis were diagnosed using serological techniques, which are limited in their cross reactivity between species, and inability to differentiate active from convalescent infections. This suggest that the results of the present study must be applied with caution. Finally, we considered only studies published in English, resulting in language bias and the possibilities of omitting credible studies published in other languages which would have provided a better insight on human babesiosis.

Conclusion

Eight Babesia species (B. bigemina, B. bovis, B. crassa-like, B. divergens, B. duncani, B. microti B. odocoilei and B. venatorum) were reported in humans from 22 countries of the world, with the highest prevalence in Europe. Babesia microti had the widest geographic distribution across the globe and was predominantly reported from the North America. To control this emerging zoonosis, we suggest the control of tick vectors around human settlements, the use of repellants and appropriate protective clothing by occupationally exposed population. The screening of blood for this parasites will also reduce the risk of transfusion-associated babesiosis.

Supplementary Material

Supplemental Material

Acknowledgments

We are grateful to all authors whose articles are included in this study. We also appreciate Mrs Juliana Tije for assisting in literature search.

Funding Statement

The author(s) reported there is no funding associated with the work featured in this article.

Disclosure statement

The authors declare that they have no competing interests.

Availability of data and materials

The data supporting the conclusion of this article are all included in the article and supplementary files 1 and 2.

Authors’ contributions

SNK: Conceived and designed the study, SNK, MNK, MIA: conducted literature search, identified articles, screened articles and extracted data. SNK: Conducted statistical and meta-analyses and wrote the manuscript. Both authors read and approved the final manuscript.

Consent for publication

Not applicable

Ethics approval and consent to participate

Not applicable

Supplementary material

Supplemental data for this article can be accessed here.

References

  • [1].Diuk-Wasser MA, Vannier E, Krause PJ.. Coinfection by Ixodes tick-borne pathogens: ecological, epidemiological, and clinical consequences. Trends Parasitol. 2016;32(1):30–42. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [2].Mehlhorn H, Schein E.. The piroplasms: life cycle and sexual stages. Adv Parasitol. 1984;23:37–103. [DOI] [PubMed] [Google Scholar]
  • [3].Becker CA, Malandrin L, Larcher T, et al. Validation of BdCCp2 as a marker for Babesia divergens sexual stages in ticks. Exp Parasitol. 2013;133(1):51–56. [DOI] [PubMed] [Google Scholar]
  • [4].Tonnetti L, Eder AF, Dy B, et al. Transfusion-transmitted Babesia microti identified through hemovigilance. Transf. 2009;49(12):2557–2563. [DOI] [PubMed] [Google Scholar]
  • [5].Herwaldt B, Linden J, Gray E, et al. Transfusion-associated babesiosis in the United States: a description of cases. Ann Int Med. 2011;155(8):509–519. [DOI] [PubMed] [Google Scholar]
  • [6].Joseph JT, Purtill K, Wong SJ, et al. Vertical transmission of Babesia microti, United States. Emerg Infect Dis. 2012;18(8):1318–1321. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [7].Tufts DM, Diuk-Wasser MA.. Transplacental transmission of tick-borne Babesia microti in its natural host peromyscus leucopus. Parasit Vectors. 2018;11(1):286. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [8].Vannier E, Krause PJ. Human babesiosis. New Engl J Med. 2012;366(25):2397–2407. [DOI] [PubMed] [Google Scholar]
  • [9].Skrabalo Z, Deanovic Z. Piroplasmosis in man; report of a case. Doc Med Geogr Trop. 1957;9:11–16. [PubMed] [Google Scholar]
  • [10].Gray J, Zintl A, Hildebrandt A, et al. Zoonotic babesiosis: overview of the disease and novel aspects of pathogen identity. Ticks Tick Borne Dis. 2010;1(1):3–10. [DOI] [PubMed] [Google Scholar]
  • [11].Zintl A, Mulcahy G, Skerrett HE, et al. Babesia divergens, a bovine blood parasite of veterinary and zoonotic importance. Clin Microbiol Rev. 2003;16(4):622–636. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [12].Dantas-Torres F, Alves LC, Uilenberg G. Babesiosis. In: Marcondes CB, editor. Arthropod borne diseases. Cham: Springer International Publishing AG; 2017. p. 347–354. [Google Scholar]
  • [13].Western KA, Benson GD, Gleason NN, et al. Babesiosis in a Massachusetts resident. N Engl J Med. 1970;283(16):854–856. [DOI] [PubMed] [Google Scholar]
  • [14].Westblade LF, Simon MS, Mathison BA, et al. Babesia microti: from mice to ticks to an increasing number of highly susceptible humans. J Clin Microbiol. 2017;55(10):2903–2912. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [15].Acosta MEP, Ender PT, Smith EM, et al. Babesia microti infection, eastern Pennsylvania, USA. Emerg Infect Dis. 2013;19(7):1105–1107. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [16].Herwaldt BL, McGovern PC, Gerwel MP, et al. Endemic babesiosis in another eastern state: new Jersey. Emerg Infect Dis. 2003;9(2):184. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [17].Stafford KC, Williams SC, Magnarelli LA, et al. Expansion of zoonotic babesiosis and reported human cases, Connecticut, 2001–2010. J Med Entomol. 2014;51(1):245–252. [DOI] [PubMed] [Google Scholar]
  • [18].Pfeiffer CD, Kazmierczak JJ, Davis JP. Epidemiologic features of human babesiosis in Wisconsin, 1996–2005. WMJ Off Publ State Med Soc Wis. 2007;106:191–195. [PubMed] [Google Scholar]
  • [19].Zamoto-Niikura A, Morikawa S, Hanaki KI, et al. Ixodes persulcatus ticks as vectors for the Babesia microti U.S. Lineage in Japan. Appl Environ Microbiol. 2016;82(22):6624–6632. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [20].Bloch EM, Kumar S, Krause PJ. Persistence of Babesia microti infection in humans. Pathogens. 2019a;8(3):102. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [21].Diuk-Wasser MA, Liu Y, Steeves T, et al. Monitoring human babesiosis emergence through vector surveillance, New England, USA. Emerg Infect Dis. 2014;20(2):225–231. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [22].Michel AO, Mathis A, Ryser-Degiorgis M. Babesia spp. in European wild ruminant species: parasite diversity and risk factors for infection. Vet Res. 2014;45(1):65. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [23].Fang LQ, Liu K, Li XL, et al. Emerging tick-borne infections in mainland China: an increasing public health threat. Lancet Infect Dis. 2015;15(12):1467–1479. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [24].Herwaldt BL, Cacciò S, Gherlinzoni F, et al. Molecular characterization of a non–babesia divergens organism causing zoonotic babesiosis in Europe. Emerg Infect Dis. 2003;9(8):942–948. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [25].Zanet S, Trisciuoglio A, Bottero E, et al. Piroplasmosis in wildlife: babesia and Theileria affecting free-ranging ungulates and carnivores in the Italian Alps. Parasit Vectors. 2014;7(1):70. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [26].Häselbarth K, Tenter AM, Brade V, et al. First case of human babesiosis in Germany - Clinical presentation and molecular characterisation of the pathogen. Int J Med Microbiol. 2007;297(3):197–204. [DOI] [PubMed] [Google Scholar]
  • [27].Chang D, Hossain M, Hossain MA. Severe babesiosis masquerading as thrombotic thrombocytopenic purpura: a case report. Cureus. 2019;11(4):e4459. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [28].Krause PJ, Gewurz BE, Hill D, et al. Persistent and relapsing babesiosis in immunocompromised patients. Clin Infect Dis. 2008;46(3):370–376. [DOI] [PubMed] [Google Scholar]
  • [29].González LM, Castro E, Lobo CA, et al. First report of Babesia divergens infection in an HIV patient. Int J Infect Dis. 2015;33:202–204. [DOI] [PubMed] [Google Scholar]
  • [30].Moher D, Liberati A, Tetzlaff J, et al. The PRISMA group. preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. PLoS Med. 2009;6(7):e1000097. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [31].World Bank . Country classification by income. 2020. https://datahelpdesk.worldbank.org/knowledgebase/articles/906519
  • [32].Munn Z, Moola S, Lisy K, et al. Methodological guidance for systematic reviews of observational epidemiological studies reporting prevalence and incidence data. Int J Evid Based Healthc. 2015;13(3):147–153. [DOI] [PubMed] [Google Scholar]
  • [33].Hedges LV, Vevea JL. Fixed- and random-effects models in meta-analysis. Psychol Meth. 1998;3(4):486–504. [Google Scholar]
  • [34].Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Meth. 2002;21:1539–1558. [DOI] [PubMed] [Google Scholar]
  • [35].Higgins JPT, Thompson SG, Deeks JJ, et al. Measuring inconsistency in meta-analyses. BMJ. 2003;327(7414):557–560. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [36].Egger M, Smith GD, Schneider M, et al. Bias in meta-analysis detected by a simple graphical test. BMJ. 1997;315(7109):629–634. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [37].Duval S, Tweedie R. Trim and fill: a simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics. 2000;56(2):455–463. [DOI] [PubMed] [Google Scholar]
  • [38].Gao L, Zhang L, Jin Q. Meta-analysis: prevalence of HIV infection and syphilis among MSM in China. Sex Transm Infect. 2009;85(5):354–358. [DOI] [PubMed] [Google Scholar]
  • [39].Bloch EM, Mrango Z, Kasubi M, et al. The Babesia observational antibody (BAOBAB) study: a cross-sectional evaluation of Babesia in two communities in Kilosa district, Tanzania. PLoS Negl Trop Dis. 2019b;13(8):e0007632. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [40].Bloch EM, Yang Y, He M, et al. A pilot serosurvey of Babesia microti in Chinese blood donors. Vox Sang. 2018;113(4):345–349. [DOI] [PubMed] [Google Scholar]
  • [41].Liu HB, Wei R, Ni XB, et al. The prevalence and clinical characteristics of tick-borne diseases at one sentinel hospital in Northeastern China. Parasitol. 2019;146(2):161–167. [DOI] [PubMed] [Google Scholar]
  • [42].Arai S, Tsuji M, Kaiho I, et al. Retrospective seroepidemiological survey for human babesiosis in an area in Japan where a tick-borne disease is endemic. J Vet Med Sci. 2003;65(3):335–340. [DOI] [PubMed] [Google Scholar]
  • [43].Mayne PJ. Emerging incidence of Lyme borreliosis, babesiosis, bartonellosis, and granulocytic ehrlichiosis in Australia. Int J Gen Med. 2011;4:845–852. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [44].Sonnleitner ST, Fritz J, Bednarska M, et al. Risk assessment of transfusion-associated babesiosis in Tyrol: appraisal by seroepidemiology and polymerase chain reaction. Transf. 2014;54(7):1725–1732. [DOI] [PubMed] [Google Scholar]
  • [45].Hunfeld KP, Lambert A, Kampen H, et al. Seroprevalence of Babesia infections in humans exposed to ticks in midwestern Germany. J Clin Microbiol. 2002;40(7):2431–2436. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [46].Pawełczyk A, Bednarska M, Kowalska JD, et al. Seroprevalence of six pathogens transmitted by the Ixodes ricinus ticks in asymptomatic individuals with HIV infection and in blood donors. Sci Rep. 2019;9(1):2117. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [47].Svensson J, Hunfeld KP, Persson KEM. High seroprevalence of Babesia antibodies among Borrelia burgdorferi-infected humans in Sweden. Ticks Tick Borne Dis. 2019;10(1):186–190. [DOI] [PubMed] [Google Scholar]
  • [48].Foppa IM, Krause PJ, Spielman A, et al. Entomologic and serologic evidence of zoonotic transmission of Babesia microti, eastern Switzerland. Emerg Infect Dis. 2002;8(7):722–726. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [49].Bondarenko AV, Torianyk II, Pokhil SI, et al. Seroprevalence of babesiosis in immunocompetent and immunocompromised individuals. Pol Merkur Lekarski. 2021;49(291):193–197. [PubMed] [Google Scholar]
  • [50].Tonnetti L, O’Brien SF, Grégoire Y, et al. Prevalence of Babesia in Canadian blood donors: June-October 2018. Transf. 2019a;59(10):3171–3176. [DOI] [PubMed] [Google Scholar]
  • [51].Asad S, Sweeney J, Mermel LA. Transfusion-transmitted babesiosis in Rhode Island. Transf. 2009;49(12):2564–2573. [DOI] [PubMed] [Google Scholar]
  • [52].Bloch EM, Levin AE, Williamson PC, et al. A prospective evaluation of chronic Babesia microti infection in seroreactive blood donors. Transf. 2016;56(7):1875–1882. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [53].Chisholm ES, Sulzer AJ, Ruebush TK 2nd. Indirect immunofluorescence test for human Babesia microti infection: antigenic specificity. Am J Trop Med Hyg. 1986;35(5):921–925. [DOI] [PubMed] [Google Scholar]
  • [54].Curcio SR, Tria LP, Gucwa AL. Seroprevalence of Babesia microti in individuals with Lyme disease. Vector Borne Zoonotic Dis. 2016;16(12):737–743. [DOI] [PubMed] [Google Scholar]
  • [55].Hilton E, DeVoti J, Benach JL, et al. Seroprevalence and seroconversion for tick-borne diseases in a high-risk population in the northeast United States. Am J Med. 1999;106(4):404–409. [DOI] [PubMed] [Google Scholar]
  • [56].Johnson ST, Cable RG, Tonnetti L, et al. Seroprevalence of Babesia microti in blood donors from Babesia-endemic areas of the northeastern United States: 2000 through 2007. Transf. 2009;49(12):2574–2582. [DOI] [PubMed] [Google Scholar]
  • [57].Johnson ST, Cable RG, Leiby DA. Lookback investigations of Babesia microti-seropositive blood donors: seven-year experience in a Babesia-endemic area. Transf. 2012;52(7):1509–1516. [DOI] [PubMed] [Google Scholar]
  • [58].Kowalski TJ, Jobe DA, Dolan EC, et al. The emergence of clinically relevant babesiosis in Southwestern Wisconsin. WMJ. 2015;114(4):152–157. [PubMed] [Google Scholar]
  • [59].Krause PJ, Telford SR, Ryan R, et al. Geographical and temporal distribution of babesial infection in connecticut. J Clin Microbiol. 1991;29(1):1–4. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [60].Krause PJ, Telford SR, Pollack RJ, et al. Babesiosis: an underdiagnosed disease of children. Pediatrics. 1992;89(6):1045–1048. [PubMed] [Google Scholar]
  • [61].Krause PJ, Telford SR, Spielman A 3rd, et al. Concurrent Lyme disease and babesiosis. Evidence for increased severity and duration of illness. JAMA. 1996;275(21):1657–1660. [PubMed] [Google Scholar]
  • [62].Leiby DA, Chung AP, Cable RG, et al. Relationship between tick bites and the seroprevalence of Babesia microti and Anaplasma phagocytophila (previously Ehrlichia sp.) in blood donors. Transf. 2002;42(12):1585–1591. [DOI] [PubMed] [Google Scholar]
  • [63].Leiby DA, Chung AP, Gill JE, et al. Demonstrable parasitemia among Connecticut blood donors with antibodies to Babesia microti. Transf. 2005;45(11):1804–1810. [DOI] [PubMed] [Google Scholar]
  • [64].Levin AE, Williamson PC, Erwin JL, et al. Determination of Babesia microti seroprevalence in blood donor populations using an investigational enzyme immunoassay. Transf. 2014;54(9):2237–2244. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [65].Levin AE, Williamson PC, Bloch EM, et al. Serologic screening of United States blood donors for Babesia microti using an investigational enzyme immunoassay. Transf. 2016;56(7):1866–1874. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [66].Linden JV, Wong SJ, Chu FK, et al. Transfusion-associated transmission of babesiosis in New York State. Transf. 2000;40(3):285–289. [DOI] [PubMed] [Google Scholar]
  • [67].Magnarelli LA, Dumler JS, Anderson JF, et al. Coexistence of antibodies to tick-borne pathogens of babesiosis, ehrlichiosis, and Lyme borreliosis in human sera. J Clin Microbiol. 1995;33(11):3054–3057. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [68].Magnarelli LA, Ijdo JW, Anderson JF, et al. Human exposure to a granulocytic Ehrlichia and other tick-borne agents in connecticut. J Clin Microbiol. 1998;36(10):2823–2827. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [69].Mitchell PD, Reed KD, Hofkes JM. Immunoserologic evidence of coinfection with Borrelia burgdorferi, Babesia microti, and human granulocytic Ehrlichia species in residents of Wisconsin and Minnesota. J Clin Microbiol. 1996;34(3):724–727. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [70].Moritz ED, Stramer SL. Blood donation screening for Babesia microti: feasibility and results. ISBT Sci Ser. 2015;10(S1):169–172. [Google Scholar]
  • [71].Moritz ED, Winton CS, Johnson ST, et al. Investigational screening for Babesia microti in a large repository of blood donor samples from non-endemic and endemic areas of the United States. Transf. 2014;54(9):2226–2236. [DOI] [PubMed] [Google Scholar]
  • [72].Moritz ED, Winton CS, Tonnetti L, et al. Screening for Babesia microti in the U.S. blood supply. N Engl J Med. 2016;375(23):2236–2245. [DOI] [PubMed] [Google Scholar]
  • [73].Popovsky MA, Lindberg LE, Syrek AL, et al. Prevalence of Babesia antibody in a selected blood donor population. Transf. 1988;28(1):59–61. [DOI] [PubMed] [Google Scholar]
  • [74].Primus S, Akoolo L, Schlachter S, et al. Efficient detection of symptomatic and asymptomatic patient samples for Babesia microti and Borrelia burgdorferi infection by multiplex qPCR. PLoS ONE. 2018;13(5):e0196748. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [75].Prince HE, Lapé-Nixon M, Patel H, et al. Comparison of the Babesia duncani (WA1) IgG detection rates among clinical sera submitted to a reference laboratory for WA1 IgG testing and blood donor specimens from diverse geographic areas of the United States. Clin Vaccine Immunol. 2010;17(11):1729–1733. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [76].Rau A, Munoz-Zanzi C, Schotthoefer AM, et al. Spatio-temporal dynamics of tick-borne diseases in North-Central Wisconsin from 2000–2016. Int J Environ Res Public Health. 2020;17(14):5105. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [77].Scott JD, Sajid MS, Pascoe EL, et al. Detection of Babesia odocoilei in humans with babesiosis symptoms. Diagnostics (Basel). 2021;11(6):947. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [78].Steere AC, McHugh G, Suarez C, et al. Prospective study of coinfection in patients with erythema migrans. Clin Infect Dis. 2003;36(8):1078–1081. [DOI] [PubMed] [Google Scholar]
  • [79].Tonnetti L, Thorp AM, Deisting B, et al. Babesia microti seroprevalence in Minnesota blood donors. Transf. 2013;53(8):1698–1705. [DOI] [PubMed] [Google Scholar]
  • [80].Tonnetti L, Townsend RL, Deisting BM, et al. The impact of Babesia microti blood donation screening. Transf. 2019b;59(2):593–600. [DOI] [PubMed] [Google Scholar]
  • [81].Wormser GP, McKenna D, Scavarda C, et al. Co-infections in persons with early lyme disease, New York, USA. Emerg Infect Dis. 2019;25(4):748–752. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [82].Gabrielli S, Totino V, Macchioni F, et al. Human babesiosis, Bolivia, 2013. Emerg Infect Dis. 2016;22(8):1445–1447. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [83].Gonzalez J, Echaide I, Pabón A, et al. Babesiosis prevalence in malaria-endemic regions of Colombia. J Vector Borne Dis. 2018;55(3):222–229. [DOI] [PubMed] [Google Scholar]
  • [84].Shih C, Liu L, Chung W, et al. Human babesiosis in Taiwan: asymptomatic infection with a Babesia microti-like organism in a Taiwanese woman. J Clin Microbiol. 1997;35(2):450–454. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [85].Dobroszycki J, Herwaldt BL, Boctor F, et al. A cluster of transfusion-associated babesiosis cases traced to a single asymptomatic donor. JAMA. 1999;281(10):927–930. [DOI] [PubMed] [Google Scholar]
  • [86].Centeno-Lima S, Do Rosario V, Parreira R, et al. A fatal case of human babesiosis in Portugal: molecular and phylogenetic analysis. Trop Med Intern Hlth. 2003;8(8):760–764. [DOI] [PubMed] [Google Scholar]
  • [87].Kukina IV, Zelya OP, Guzeeva TM, et al. Severe babesiosis caused by Babesia divergens in a host with intact spleen, Russia, 2018. Tick Tick Borne Dis. 2019;10(6):101262. [DOI] [PubMed] [Google Scholar]
  • [88].Nassar Y, Richter S. Babesiosis presenting as acute liver failure. Case Rep Gastroenterol. 2017;11(3):769–773. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [89].Rajasingham R, Williams D, Meya DB, et al. Severe babesiosis in immunocompetent man, Spain, 2011. Emerg Infect Dis. 2014;20(4):724–726. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [90].Arsuaga M, Gonzalez LM, Lobo CA, et al. First report of Babesia microti-caused babesiosis in Spain. Vector Borne Zoon Dis. 2016;16(10):677–679. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [91].Lim P, Chavatte J, Vasoo S, et al. Imported human babesiosis, Singapore, 2018. Emerg Infect Dis. 2020;26(4):826–828. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [92].Stahl P, Poinsignon Y, Pouedras P, et al. Case report of the patient source of the Babesia microti R1 reference strain and implications for travellers. J Travel Med. 2017;25(1):1–8. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [93].Vannier EG, Diuk-Wasser MA, Mamoun CB, et al. Babesiosis. Infect Dis Clin North Am. 2015;29(2):357–370. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [94].Young KM, Corrin T, Wilhelm B, et al. Zoonotic Babesia: a scoping review of the global evidence. PLoS One. 2019;14(12):e0226781. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [95].Lempereur L, Shiels B, Heyman P, et al. A retrospective serological survey on human babesiosis in Belgium. Clin Microbiol Infect. 2015;21(1):96.e1–96.e967. [DOI] [PubMed] [Google Scholar]
  • [96].Topolovec J, Puntarić D, Antolović-Pozgain A, et al. Serologically detected “new” tick-borne zoonoses in eastern Croatia. Croat Med J. 2003;44(5):626–629. [PubMed] [Google Scholar]
  • [97].Ocias LF, Wilhelmsson P, Sjowall J, et al. Emerging tick-borne pathogens in the nordic countries: a clinical and laboratory follow-up study of high-risk tick-bitten individuals. Ticks Tick Borne Dis. 2020;11(1):101303. [DOI] [PubMed] [Google Scholar]
  • [98].Wilhelmsson P, Lövmar M, Krogfelt KA, et al. Clinical/serological outcome in humans bitten by Babesia species positive Ixodes ricinus ticks in Sweden and on the Åland Islands. Ticks Tick Borne Dis. 2020;11(4):101455. [DOI] [PubMed] [Google Scholar]
  • [99].Rigaud E, Jaulhac B, Garcia-Bonnet N, et al. Seroprevalence of seven pathogens transmitted by the Ixodes ricinus tick in forestry workers in France. Clin Microbiol Infect. 2016;22(8):735.e1–735.e7359. [DOI] [PubMed] [Google Scholar]
  • [100].Gabrielli S, Calderini P, Cassini R, et al. Human exposure to piroplasms in central and Northern Italy. Vet Ital. 2014;50(1):41–47. [DOI] [PubMed] [Google Scholar]
  • [101].Jahfari S, Hofhuis A, Fonville M, et al. Molecular detection of tick-borne pathogens in humans with tick bites and erythema migrans, in the Netherlands. PLoS Negl Trop Dis. 2016;10(10):e0005042. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [102].Thortveit ET, Aase A, Petersen LB, et al. Human seroprevalence of antibodies to tick-borne microbes in southern Norway. Ticks Tick Borne Dis. 2020;11(4):101410. [DOI] [PubMed] [Google Scholar]
  • [103].Chmielewska-Badora J, Moniuszko A, Żukiewicz-Sobczak W, et al. Serological survey in persons occupationally exposed to tick-borne pathogens in cases of co-infections with Borrelia burgdorferi, Anaplasma phagocytophilum, Bartonella spp. and Babesia microti. Ann Agric Environ Med. 2012;19(2):271–274. [PubMed] [Google Scholar]
  • [104].Moniuszko A, Dunaj J, Swięcicka I, et al. Co-infections with Borrelia species, Anaplasma phagocytophilum and Babesia spp. in patients with tick-borne encephalitis. Eur J Clin Microbiol Infect Dis. 2014;33(10):1835–1841. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [105].Moniuszko-Malinowska A, Swiecicka I, Dunaj J, et al. Infection with Babesia microti in humans with non-specific symptoms in North East Poland. Infect Dis. 2016;48(7):537–543. [DOI] [PubMed] [Google Scholar]
  • [106].Pancewicz S, Moniuszko A, Bieniarz E, et al. Anti-Babesia microti antibodies in foresters highly exposed to tick bites in Poland. Scand J Infect Dis. 2011;43(3):197–201. [DOI] [PubMed] [Google Scholar]
  • [107].Pañczuk A, Tokarska-Rodak M, Kozioł-Montewka M, et al. The incidence of Borrelia burgdorferi, anaplasma phagocytophilum and Babesia microti coinfections among foresters and farmers in eastern Poland. J Vector Borne Dis. 2016;53(4):348–354. [PubMed] [Google Scholar]
  • [108].Welc-Falęciak R, Hildebrandt A, Siński E. Co-infection with Borrelia species and other tick-borne pathogens in humans: two cases from Poland. Ann Agric Environ Med. 2010;17:309–313. [PubMed] [Google Scholar]
  • [109].Welc-Falęciak R, Pawełczyk A, Radkowski M, et al. First report of two asymptomatic cases of human infection with Babesia microti (Franca, 1910) in Poland. Ann Agric Environ Med. 2015;22(1):51–54. [DOI] [PubMed] [Google Scholar]
  • [110].Żukiewicz-Sobczak W, Zwoliński J, Chmielewska-Badora J, et al. Prevalence of antibodies against selected zoonotic agents in forestry workers from eastern and southern Poland. Ann Agric Environ Med. 2014;21(4):767–770. [DOI] [PubMed] [Google Scholar]
  • [111].Filstein MR, Benach JL, White DJ, et al. Serosurvey for human babesiosis in New York. J Infect Dis. 1980;141(4):518–521. [DOI] [PubMed] [Google Scholar]
  • [112].Zhou X, Li S, Chen S, et al. Co-infections with Babesia microti and Plasmodium parasites along the China-Myanmar border. Infect Dis Pov. 2013;2(1):24. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [113].Qi C, Zhou D, Liu J, et al. Detection of Babesia divergens using molecular methods in anemic patients in Shandong Province, China. Parasitol Res. 2011;109(1):241–245. [DOI] [PubMed] [Google Scholar]
  • [114].Wang J, Zhang S, Yang J, et al. Babesia divergens in human in Gansu province, China. Emerg Microbes Infect. 2019;8(1):959–961. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [115].Gigler B. Poverty, inequality and human development of indigenous peoples in Bolivia. Working Paper Series No. 17;2009. [Google Scholar]
  • [116].Krause PJ, Telford SR, Ryan R, et al. Diagnosis of babesiosis: evaluation of a serologic test for the detection of Babesia microti antibody. J Infect Dis. 1994;169(4):923–926. [DOI] [PubMed] [Google Scholar]
  • [117].Hong SH, Anu D, Jeong YI, et al. Molecular detection and seroprevalence of Babesia microti among stock farmers in Khutul City, Selenge Province, Mongolia. Korean J Parasitol. 2014;52(4):443–447. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [118].Johnson ST, Van Tassell ER, Tonnetti L, et al. Babesia microti real-time polymerase chain reaction testing of connecticut blood donors: potential implications for screening algorithms. Transf. 2013;53(11):2644–2649. [DOI] [PubMed] [Google Scholar]
  • [119].Linden JV, Prusinski MA, Crowder LA, et al. Transfusion-transmitted and community-acquired babesiosis in New York, 2004 to 2015. Transf. 2018;58(3):660–668. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [120].Herwaldt B, Montgomery S, Woodhall D, et al. Babesiosis surveillance – 18 States, 2011. MMWR. 2012;61:505–509. [PubMed] [Google Scholar]
  • [121].Yabsley MJ, Shock BC. Natural history of zoonotic Babesia: role of wildlife reservoirs. IJP: Parasit Wildl. 2013;2:18–31. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [122].Silaghi C, Woll D, Hamel D, et al. Babesia spp. and Anaplasma phagocytophilum in questing ticks, ticks parasiting rodents and the parasitized rodents – analyzing the host-pathogen vector interface in metropolitan area. Parasit Vectors. 2012;5(1):191. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • [123].Jiang JF, Zheng YC, Jiang RR, et al. Epidemiological, clinical, and laboratory characteristics of 48 cases of “Babesia venatorum ” infection in China: a descriptive study. Lancet Infect Dis. 2015;15(2):196–203. [DOI] [PubMed] [Google Scholar]
  • [124].Jia N, Zheng YC, Jiang JF, et al. Human babesiosis caused by a Babesia crassa-like pathogen: a case series. Clin Infect Dis. 2018;67(7):1110–1119. [DOI] [PubMed] [Google Scholar]
  • [125].Royle P, Waugh N. Literature searching for clinical and cost-effectiveness studies used in health technology assessment reports carried out for the national institute for clinical excellence appraisal system. Hlth Technol Assess. 2003;7:1–51. [DOI] [PubMed] [Google Scholar]
  • [126].Sutton AJ. Publication bias. In: Cooper H, Hedges LV, Valentine JC, editors. The handbook of research synthesis and meta-analysis. 2nd ed. New York: Russell Sage Foundation; 2009. p. 435–452. [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Supplemental Material

Data Availability Statement

The data supporting the conclusion of this article are all included in the article and supplementary files 1 and 2.


Articles from Pathogens and Global Health are provided here courtesy of Taylor & Francis

RESOURCES