Search

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.

The burden of tick-borne diseases continues to increase in the United States. Tick surveillance has been implemented to monitor changes in the distribution and prevalence of human disease-causing pathogens in ticks that frequently bite humans. Such efforts require accurate identification of ticks to species and highly sensitive and specific assays that can detect and differentiate pathogens from genetically similar microbes in ticks that have not been demonstrated to be pathogenic in humans. We describe a modification to a next generation sequencing pathogen detection assay that includes a target that accurately identifies Ixodes ticks to species. We show that the replacement of internal control primers used to ensure assay performance with primers that also act as an internal control and can additionally differentiate tick species, retains high sensitivity and specificity, improves efficiency, and reduces costs by eliminating the need to run separate assays to screen for pathogens and for tick identification.

Lyme disease is the most commonly reported vector-borne disease in the United States and is transmitted by Ixodes scapularis in the eastern US and I. pacificus in the west. The causative agents, Borrelia burgdorferi sensu stricto (Bbss) and B. mayonii belong to the B. burgdorferi sensu lato (Bbsl) species complex. An additional eight species of Bbsl have been identified in Ixodes species ticks in the US, but their geographic distribution, vector associations, human encounter rates and pathogenicity in humans are poorly defined. To better understand the geographic distribution and vector associations of Bbsl spirochetes in frequent and infrequent human-biting Ixodes species ticks in the US, we previously screened 29,517 host-seeking I. scapularis or I. pacificus ticks and 692 ticks belonging to eight other Ixodes species for Borrelia spirochetes using a previously described tick testing algorithm that utilizes a combination of real-time PCR and Sanger sequencing for Borrelia species identification. The assay was designed to detect known human pathogens spread by Ixodes species ticks, but it was not optimized to detect Bbsl co-infections. To determine if such co-infections were overlooked particularly in ticks infected with Bbss, we retested and analyzed a subsample of 845 Borrelia infected ticks using a next generation sequencing multiplex PCR amplicon sequencing (MPAS) assay that can identify Borrelia species and Bbsl co-infections. The assay also includes targets that can molecularly confirm identifications of Ixodes species ticks to better inform pathogen-vector associations. We show that Bbss is the most prevalent species in I. scapularis and I. pacificus; other Bbsl species were rarely detected in I. scapularis and the only Bbsl co-infections identified in I. scapularis were with Bbss and B. mayonii. We detected B. andersonii in I. dentatus in the Mid-Atlantic and Upper Midwest regions, B. kurtenbachii in I. scapularis in the Upper Midwest, B. bissettiae in I. pacificus and I. spinipalpis in the Northwest, and B. carolinensis in I. affinis in the Mid-Atlantic and Southeast, and B. lanei in I. spinipalpis in the Northwest. Twelve of 62 (19.4%) Borrelia-infected I. affinis from the Mid-Atlantic region were co-infected with Bbss and B. carolinensis. Our data support the notion that Bbsl species are maintained in largely independent enzootic cycles, with occasional spill-over resulting in multiple Bbsl species detected in Ixodes species ticks.

The genus Bartonella includes a group of species that are associated with a wide range of mammalian species, including human. It is challenging to detect all Bartonella species using a single molecular target due to its high genetic diversity. To solve this issue, we developed a quadruplex PCR amplicon sequencing assay using next-generation sequencing (NGS) technology for the detection and differentiation of Bartonella species. Our objective was to obtain the specific sequences of a minimum of two of the four target genes as confirmation of the identity of a particular Bartonella species using the assay. Four pairs of primers targeting specific regions on gltA, groEL, rpoB, and ssrA were evaluated for their capability of differentiating Bartonella species individually and collectively by performing singular PCR amplicon sequencing and quadruplex PCR amplicon sequencing. Using the quadruplex PCR amplicon sequencing, 24 Bartonella reference species were tested, all of which were successfully differentiated by at least two targets. Bartonella species were accurately identified from the artificially mixed DNA templates developed to simulate coinfections. The limit of detection was determined to be 1 fg based on testing a series of 10-fold dilutions of DNA from the Bartonella species. Testing of high DNA concentrations of 19 non-Bartonella species showed high specificity with none of the non-Bartonella species misclassified as Bartonella. Finally, the assay was evaluated by testing DNA extracts from field-collected body lice (Pediculus humanus humanus) and Norway rats (Rattus norvegicus): Bartonella quintana was detected and confirmed by three targets in the lice and Bartonella tribocorum was detected and confirmed by two targets in the rats. These results demonstrated that Bartonella species could be accurately and rapidly detected and differentiated into different tissue types using the quadruplex sequencing assay.

The Centers for Disease Control and Prevention's national tick and tick-borne pathogen surveillance program collects information to better understand the regional distribution, prevalence, and exposure risk of host-seeking medically important ticks in the United States. A recently developed next generation sequencing (NGS) targeted multiplex PCR amplicon sequencing (MPAS) assay has enhanced the detection capabilities for Ixodes-associated human pathogens found in Ixodes scapularis and Ixodes pacificus ticks compared to the routinely used real-time PCR assay. To operationalize the MPAS assay for the large number of tick surveillance submissions processed each year, a reproducible high throughput bioinformatics pipeline is needed. We describe the development and validation of the MPAS pipeline, a bioinformatics pipeline that identifies and summarizes amplicon sequences produced by the MPAS assay. This pipeline is portable and reproducible across different computing environments, and flexible by allowing modifications to input parameters, assay primer and reference sequences. The automation of the summary report, BLAST report, and phylogenetic analysis reduces the amount of time needed for downstream analysis. To validate this pipeline, we compared the analysis of a MPAS assay dataset consisting of 175 I. scapularis nymphs with the MPAS pipeline and previously published results analyzed with a CLC Genomic Workbench workflow. The MPAS pipeline identified the same number of positive ticks for Anaplasma phagocytophilum and Babesia species as the original analysis, but the MPAS pipeline provided enhanced sequencing resolution of Borrelia burgdorferi sensu lato co-infected samples. The reproducibility, flexibility, analysis automation, and improved sequence resolution of the MPAS pipeline make it well suited for a high throughput tick pathogen surveillance program.