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.
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.
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.
Country-based pooled estimates of Babesia infections in humans
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.
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
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
Data Availability Statement
The data supporting the conclusion of this article are all included in the article and supplementary files 1 and 2.