A global map of genetic diversity in Babesia microti reveals strong population structure and identifies variants associated with clinical relapse

Human babesiosis caused by Babesia microti is an emerging zoonosis of increasing importance due to rising incidence and expanding geographic range. To investigate the origin and genetic diversity of this parasite, we sequenced the complete genomes of 42 B. microti samples from around the world, including deep coverage of clinical infections at endemic sites in the continental United States. Samples from the continental US segregate into a Northeast lineage and a Midwest lineage, with subsequent divergence of subpopulations along geographic lines. We identify parasite variants that associate with relapsing disease, including amino acid substitutions in the atovaquone-binding region of cytochrome b (cytb) and the azithromycin-binding region of ribosomal protein subunit L4 (rpl4). Our results shed light on the origin, diversity, and evolution of B. microti, suggest possible mechanisms for clinical relapse, and create the foundation for further research on this emerging pathogen.