Geographic variation in the distribution of Anaplasma phagocytophilum variants in host-seeking Ixodes scapularis nymphs and adults in the eastern United States elucidated using next generation sequencing.

Human anaplasmosis cases, caused by Anaplasma phagocytophilum, are increasing in the United States. This trend is explained, in part, by expansion in the geographic range of the primary vector, Ixodes scapularis. Multiple variants of A. phagocytophilum have been identified in field collected ticks, but only a single variant (human active, or "Ap-ha," variant) has been shown to be pathogenic in humans. Until recently, laboratory methods used to differentiate variants were cumbersome and seldomly used in large scale assessments of the pathogen's geographic distribution. As a result, many surveys reported A. phagocytophilum without segregating variants. Lack of discrimination among A. phagocytophilum variants could lead to overestimation of anaplasmosis risk to humans. Next Generation Sequencing (NGS) assays were recently developed to efficiently detect multiple Ixodes scapularis-borne human pathogens including Ap-ha. In this study, we utilized NGS to detect and differentiate A. phagocytophilum variants (Ap-ha vs. non ha) in host-seeking I. scapularis nymphs and adults collected across 23 states in the eastern United States from 2012 to 2023 as part of national tick surveillance efforts and research studies. Many of the included ticks were tested previously using a TaqMan PCR assay that could detect A. phagocytophilum but could not differentiate variants. We retested A. phagocytophilum infected ticks with NGS to differentiate variants. Anaplasma phagocytophilum (any variant) was identified in 165 (35 %) of 471 counties from which ticks were tested, whereas Ap-ha was detected in 70 (15 %) of 469 counties where variants were differentiated. Both variants were identified in 32 % (n = 40) of 126 counties with either variant detected. Among states where A. phagocytophilum (any variant) was detected, prevalence ranged from 2 % to 19 % in unfed adults and from 0.2 % to 7.8 % in unfed nymphs; prevalence of Ap-ha variant ranged from 0.0 % to 16 % in adults, and 0.0 % to 4.6 % in nymphs.