Hostname: page-component-cd9895bd7-dk4vv Total loading time: 0 Render date: 2024-12-23T11:58:43.489Z Has data issue: false hasContentIssue false
  • English
  • Français

Multiple piroplasm parasites (Apicomplexa: Piroplasmida) in northeastern populations of the invasive Asian longhorned tick, Haemaphysalis longicornis Neumann (Ixodida: Ixodidae), in the United States

Published online by Cambridge University Press:  04 October 2023

Heidi Herb ,
Julia González ,
Francisco C. Ferreira Open the ORCID record for Francisco C. Ferreira [Opens in a new window]  and
Dina M. Fonseca Open the ORCID record for Dina M. Fonseca [Opens in a new window]
Heidi Herb
Affiliation:
Center for Vector Biology, Rutgers University, New Brunswick, NJ 08901, USA Department of Ecology and Evolution, Rutgers University, New Brunswick, NJ 08901, USA
Julia González
Affiliation:
Center for Vector Biology, Rutgers University, New Brunswick, NJ 08901, USA
Francisco C. Ferreira
Affiliation:
Center for Vector Biology, Rutgers University, New Brunswick, NJ 08901, USA
Dina M. Fonseca*
Affiliation:
Center for Vector Biology, Rutgers University, New Brunswick, NJ 08901, USA Department of Ecology and Evolution, Rutgers University, New Brunswick, NJ 08901, USA Department of Entomology, Rutgers University, New Brunswick, NJ 08901, USA
*
Corresponding author: Dina M. Fonseca; Email: dina.fonseca@rutgers.edu
Rights & Permissions [Opens in a new window]

Abstract

Piroplasms, which include the agents of cattle fever and human and dog babesiosis, are a diverse group of blood parasites of significant veterinary and medical importance. The invasive Asian longhorned tick, Haemaphysalis longicornis, is a known vector of piroplasms in its native range in East Asia and invasive range in Australasia. In the USA, H. longicornis has been associated with Theileria orientalis Ikeda outbreaks that caused cattle mortality. To survey invasive populations of H. longicornis for a broad range of piroplasms, 667 questing H. longicornis collected in 2021 from 3 sites in New Jersey, USA, were tested with generalist piroplasm primers targeting the 18S small subunit rRNA (395–515 bp, depending on species) and the cytochrome b oxidase loci (1009 bp). Sequences matching Theileria cervi type F (1 adult, 5 nymphs), an unidentified Theileria species (in 1 nymph), an undescribed Babesia sensu stricto (‘true’ Babesia, 2 adults, 2 nymphs), a Babesia sp. Coco (also a ‘true Babesia’, 1 adult, 1 nymph), as well as Babesia microti S837 (1 adult, 4 nymphs) were recovered. Babesia microti S837 is closely related to the human pathogen B. microti US-type. Additionally, a 132 bp sequence matching the cytochrome b locus of deer, Odocoileus virginanus, was obtained from 2 partially engorged H. longicornis. The diverse assemblage of piroplasms now associated with H. longicornis in the USA spans 3 clades in the piroplasm phylogeny and raises concerns of transmission amplification of veterinary pathogens as well as spillover of pathogens from wildlife to humans.

Keywords

agricultureinvasive spilloverlivestockOne Healthpublic healthwildlife
Type
Research Article
Information
Parasitology , Volume 150 , Issue 11 , September 2023 , pp. 1063 - 1069
Check for updates
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press

Introduction

Invasive ticks, mosquitoes and other blood-feeding arthropods may introduce and transmit (i.e. vector) exotic pathogens for which local populations have little or no immunity. Resulting disease can range from mild to severe to fatal and can have a significant impact on human health (e.g. Zika fever, West Nile virus encephalitis), animal health (e.g. redwater fever, blue tongue virus), hasten extinction (e.g. bird malaria) and cause economic damage to agriculture, tourism and other industries (Athni et al., Reference Athni, Shocket, Couper, Nova, Caldwell, Caldwell, Childress, Childs, De Leo, Kirk, MacDonald, Olivarius, Pickel, Roberts, Winokur, Young, Cheng, Grant, Kurzner, Kyaw, Lin, Lopez, Massihpour, Olsen, Roache, Ruiz, Schultz, Shafat, Spencer, Bharti and Mordecai2021). Invasive vectors can also potentially spread existing wildlife pathogens by creating new transmission pathways, which can have significant ecological and public health implications, particularly in the context of One Health (sensu Lerner and Berg, Reference Lerner and Berg2015).

The phylum Apicomplexa includes well-known blood-borne protozoa such as Plasmodium falciparum and P. vivax, the primary agents of human malaria (Votýpka et al., Reference Votýpka, Modrý, Oborník, Šlapeta and Lukeš2016). The Apicomplexa class Piroplasmida includes Babesia, Theileria and Cytauxzoon that are primarily transmitted by hard ticks (Ixodida: Ixodidae) and can affect a wide range of hosts (Onyiche et al., Reference Onyiche, Răileanu, Fischer and Silaghi2021; Almazán et al., Reference Almazán, Scimeca, Reichard and Mosqueda2022). While piroplasms were once classified based on morphology and host associations alone, the advent of molecular methods has greatly advanced the overall understanding of the diversity of Piroplasmida (Garrett et al., Reference Garrett, Hernandez, Balsamo, Barron, Beasley, Brown, Cloherty, Farid, Gabriel, Groves, Hamer, Hill, Lewis, McManners, Nemeth, Oesterle, Ortiz, Peshock, Schnellbacher, Schott, Straif-Bourgeois and Yabsley2019). A recent analysis by Jalovecka et al. (Reference Jalovecka, Sojka, Ascencio and Schnittger2019) indicates that there are at least 10 distinct clades within Piroplasmida, with both Babesia spp. and Theileria spp. comprising polyphyletic groups in need of taxonomic revision.

Babesia spp. are broadly divided into Babesia sensu stricto and Babesia sensu lato, with the former representing a monophyletic group considered as ‘true Babesia’, distinguishable from other piroplasms by their ability to infect the reproductive organs of the tick and to be transmitted to the eggs (transovarial transmission) (Schreeg et al., Reference Schreeg, Marr, Tarigo, Cohn, Bird, Scholl, Levy, Wiegmann and Birkenheuer2016; Jalovecka et al., Reference Jalovecka, Sojka, Ascencio and Schnittger2019; Schnittger et al., Reference Schnittger, Ganzinelli, Bhoora, Omondi, Nijhof and Florin-Christensen2022). Within Babesia sensu lato, one of the best-characterized clades is the Babesia microti group, which contains the piroplasms responsible for most human babesiosis cases worldwide, especially in the northern United States (Renard and Mamoun, Reference Renard and Mamoun2021). Importantly, B. microti includes at least 2 different genetic lineages pathogenic to humans and several only known from reservoir hosts such as mice, voles and skunks (Goethert, Reference Goethert2021). Moreover, although still relatively rare, human disease caused by other Babesia species such as B. divergens in Europe, B. venatorum in Asia and Europe and B. duncani in North America has been increasingly reported (Scott and Scott, Reference Scott and Scott2018; Hong et al., Reference Hong, Kim, Song, Roh, Cho, Kim, Um, Kwak, Cho and Lee2019; Kumar et al., Reference Kumar, O'Bryan and Krause2021; Scott et al., Reference Scott, Sajid, Pascoe and Foley2021).

Rates of human babesiosis have been increasing in the USA, particularly in the northeastern states (Almazán et al., Reference Almazán, Scimeca, Reichard and Mosqueda2022; Swanson et al., Reference Swanson, Pickrel, Williamson and Montgomery2023). While those infected with Babesia may experience fever, chills, headache, muscle aches, fatigue and red or brown urine, some may not have any symptoms at all, especially if their immune systems are not compromised (Almazán et al., Reference Almazán, Scimeca, Reichard and Mosqueda2022). People may remain infected for years, and even if asymptomatic, can transmit piroplasms through blood transfusions or organ transplants (Bloch et al., Reference Bloch, Kumar and Krause2019), and transmission from infected mothers to developing fetuses has been demonstrated (Horowitz and Freeman, Reference Horowitz and Freeman2020). Piroplasmid infections are typically treated with various combinations of atovaquone, azithromycin, clindamycin and quinine; however, concerns regarding side-effects, drug resistance and drug efficacy indicate the need for development of novel treatment options (Renard and Mamoun, Reference Renard and Mamoun2021).

Bovine babesiosis caused by Babesia bigemina and Babesia bovis has long been a concern to American cattle ranchers and several Theileria species can sicken horses, cervids and bovids (Almazán et al., Reference Almazán, Scimeca, Reichard and Mosqueda2022; Osbrink et al., Reference Osbrink, Thomas, Lohmeyer and Temeyer2022). Both babesiosis and theileriosis cause significant economic losses annually to North American and Australasian agricultural industries due to reduced production, death, abortions, restrictions on animal movement and costs associated with preventive measures and treatments (Dinkel et al., Reference Dinkel, Herndon, Noh, Lahmers, Todd, Ueti, Scoles, Mason and Fry2021; Almazán et al., Reference Almazán, Scimeca, Reichard and Mosqueda2022; Osbrink et al., Reference Osbrink, Thomas, Lohmeyer and Temeyer2022; Schnittger et al., Reference Schnittger, Ganzinelli, Bhoora, Omondi, Nijhof and Florin-Christensen2022). Companion animals are also at risk, with infections of Babesia vulpes, Babesia conradae, Babesia vogeli, Babesia gibsoni and Babesia sp. Coco capable of causing mild to severe disease in dogs in the United States (Dear and Birkenheuer, Reference Dear and Birkenheuer2022).

Since the initial discovery of the invasive Asian longhorned tick (Haemaphysalis longicornis) in the United States in 2017 (Rainey et al., Reference Rainey, Occi, Robbins and Egizi2018), there have been concerns regarding the potential threats this ectoparasite may pose. In North America, as in Australasia where it expanded to in the early 20th century, H. longicornis reproduces asexually (clonally) by parthenogenesis (Schappach et al., Reference Schappach, Krell, Hornbostel and Connally2020), which underlies the ability of this species to develop large populations very quickly. In its Australasian range, H. longicornis represents a major threat to domestic livestock, heavily parasitizing large ruminants, and impeding production (Heath, Reference Heath2016).

Globally, H. longicornis is a known vector of piroplasms that infect humans, livestock and companion animals, including B. ovata, B. gibsoni, B. microti, T. uilenbergi and T. orientalis (Li et al., Reference Li, Luo, Guan, Ma, Liu, Liu, Ren, Niu, Lu, Gao, Liu, Dang, Tian, Zhang, He, Bai and Yin2009; Wu et al., Reference Wu, Cao, Zhou, Zhang, Gong and Zhou2017; Gray et al., Reference Gray, Estrada-Peña and Zintl2019; Dinkel et al., Reference Dinkel, Herndon, Noh, Lahmers, Todd, Ueti, Scoles, Mason and Fry2021; Dear and Birkenheuer, Reference Dear and Birkenheuer2022) and may also be a vector of Babesia caballi, the agent of equine babesiosis (Bautista et al., Reference Bautista, Ikadai, You, Battsetseg, Igarashi, Nagasawa and Fujisaki2001). As in Australia (Marendy et al., Reference Marendy, Baker, Emery, Rolls and Stutchbury2020), in Virginia, USA, H. longicornis has already been implicated in the transmission of the virulent Theileria orientalis Ikeda genotype that resulted in multiple cattle deaths (Thompson et al., Reference Thompson, White, Shaw, Egizi, Lahmers, Ruder and Yabsley2020; Dinkel et al., Reference Dinkel, Herndon, Noh, Lahmers, Todd, Ueti, Scoles, Mason and Fry2021).

Although humans are not favoured hosts of H. longicornis, opportunistic feeding is well documented both in the native and invasive ranges (Bickerton and Toledo, Reference Bickerton and Toledo2020; Wormser et al., Reference Wormser, McKenna, Piedmonte, Vinci, Egizi, Backenson and Falco2020). In East Asia, H. longicornis vectors severe fever with thrombocytopenia syndrome virus (SFTSV), an emerging human tick-borne disease recently reclassified as Dabie bandavirus (Liu et al., Reference Liu, Zhou, Sun, Yao, Li, Wang, Mu, Li, Yang, Gray, Cui, Yin, Fang, Yu and Cao2015; Luo et al., Reference Luo, Zhao, Wen, Zhang, Liu, Fang, Xue, Ma, Zhang, Ding, Lei and Yu2015; Li et al., Reference Li, Liu, Wu, Liu, Huang, Li, Liu, Li, Wang, Li and Liang2021). Under laboratory conditions, US lineages of H. longicornis can vector the closely related Heartland virus (Raney et al., Reference Raney, Perry and Hermance2022a) as well as Powassan virus (Raney et al., Reference Raney, Herslebs, Langohr, Stone and Hermance2022b), 2 native pathogenic viruses emergent in parts of the USA. They can also vector Rickettsia rickettsii, the causative agent of Rocky Mountain spotted fever (Stanley et al., Reference Stanley, Ford, Snellgrove, Hartzer, Smith, Krapiunaya and Levin2020). In addition, Bourbon virus has been detected from a larval pool, 2 nymphs and 1 adult field collected H. longicornis in Virginia, USA (Cumbie et al., Reference Cumbie, Trimble and Eastwood2022). While H. longicornis is currently not perceived as a major public health threat in the USA, this status may change given the enormous densities it can reach in favourable habitats (Bickerton and Toledo, Reference Bickerton and Toledo2020; Schappach et al., Reference Schappach, Krell, Hornbostel and Connally2020; González et al., Reference González, Fonseca and Toledo2023; Rochlin et al., Reference Rochlin, Egizi, Narvaez, Bonilla, Gallagher, Williams, Rainey, Price and Fonseca2023).

The objective of this study was to assess the potential role of H. longicornis as a vector of piroplasms in New Jersey (NJ), the most urbanized US state that, maybe surprisingly to many, also boasts the highest density of horses (Rankins and Malinowski, Reference Rankins and Malinowski2020).

Materials and methods

Study areas

This study was conducted in 3 sites approximately 1.2 km from each other within the Rutgers University Cook Campus in New Brunswick, NJ (please refer to Ferreira et al., Reference Ferreira, González, Milholland, Tung and Fonseca2023 for a map). Surveys for H. longicornis were initiated at these sites in 2018 when the species was first detected on a grassy area next to a goat pen (Egizi et al., Reference Egizi, Occi, Price and Fonseca2019b), a site that became known as the ‘Goat Farm’ (40.47444° N, 74.43683° W). The ‘Rutgers Gardens’ site (40.47455° N, 74.42030° W) is inside a 180-acre botanical garden, consisting of designed gardens, plant collections and natural habitats. Finally, the ‘University Inn’ site (40.48413 N, 74.43051 W) is a meadow and forested park behind the Rutgers University Inn & Conference Center. At all sites, local forest is dominated by oak and maple trees and huckleberry and blueberry shrubs (Breden et al., Reference Breden, Alger, Strakosch Walz and Windisch2001), with grassy ecotones.

Tick surveillance

From June through September 2021, concomitant with surveys for ticks on mammals at the same sites, questing ticks were sampled from 50–75 m2 of the vegetation at each site (Ferreira et al., Reference Ferreira, González, Milholland, Tung and Fonseca2023). Tick sampling was performed using a white crib flannel sweep measuring 50 × 100 cm with a PVC pipe handle (Egizi et al., Reference Egizi, Occi, Price and Fonseca2019b). The sweep cloth was checked in 1–2 m intervals since H. longicornis does not attach firmly to the flannel and often drops off over longer intervals (Bickerton et al., Reference Bickerton, McSorley and Toledo2021). Ticks were collected from both sides of the sweep and identified morphologically in the laboratory to the species level using a stereomicroscope (Leica S8 APO, Leica Microsystems, Deerfield, IL, USA) following appropriate taxonomical keys (Keirans and Litwak, Reference Keirans and Litwak1989; Egizi et al., Reference Egizi, Robbins, Beati, Nava, Evans, Occi and Fonseca2019a). The larvae of H. longicornis were not stored during these surveys and were not available for pathogen testing. A few questing H. longicornis that were found partially engorged (sensu Price et al., Reference Price, Witmier, Eckert and Boyer2022) were processed separately (see section below on ‘Bloodmeal analysis of partially engorged specimens’).

DNA extraction and pathogen detection

Each tick was placed in 180 μl of Qiagen buffer ATL with 20 μl Qiagen Proteinase K (10 mg mL−1) in microfuge tubes and homogenized with a 5 mm sterile glass bead (Fisher Scientific, Waltham, MA, USA) in a TissueLyser (Qiagen Inc., Valencia, CA, USA). DNA from individual ticks was extracted using Dneasy Blood and Tissue 96-well plate kits (Qiagen Inc., Valencia, CA, USA) following the manufacturer's instructions. DNA was eluted from each column twice with 50 μl of Qiagen's elution buffer AE into separate labelled microtubes.

After reviewing the literature, primers were chosen targeting the multi copy 18S rRNA locus (Table 1) to match a broad range of piroplasm species (Casati et al., Reference Casati, Sager, Gern and Piffaretti2006) and all ticks were tested individually. To further characterize an undescribed new Babesia sp. detected, specimens positive for that Babesia were also tested using primers targeting the mitochondrial cytochrome oxidase b (cytb) locus (Table 1) shown to work across multiple Babesia species (Rajapakshage et al., Reference Rajapakshage, Yamasaki, Hwang, Sasaki, Murakami, Tamura, Lim, Nakamura, Ohta and Takiguchi2012).

Table 1. Primers used for piroplasm detection in H. longicornis

Ta = annealing temperature. The expected 18S rRNA amplicon size varies depending on the species of piroplasm from 395 to 515 bp; the expected cytochrome b amplicon size is 1009 bp.

*Instead of a specific Ta, a ‘touch-down’ approach was used starting with Ta = 54°C and decreasing the Ta by 1 degree for 4 additional cycles. This was followed by 40 additional cycles at Ta = 52°C.

The targeted loci were amplified in 20 μl reactions with Amplitaq Gold Master Mix (ThermoFisher Scientific, Waltham, MA, USA) following the manufacturer's protocol. After visualizing the amplification in a 1% agarose gel, polymerase chain reactions (PCR) were cleaned with ExoSAP-IT (ThermoFisher Scientific) and Sanger sequenced separately with both primers at Azenta Genewiz (South Plainfield, NJ, USA). The sequences were trimmed and aligned with Geneious Prime 2023.0.1 (Biomatters Inc., San Diego, CA, USA) and the consensus was used as a query in NCBI's Basic Local Alignment Search Tool, BLASTn (Altschup et al., Reference Altschup, Gish, Miller, Myers and Lipman1990).

Phylogenetic analysis

Consensus sequences were aligned to available NCBI Genbank sequences (Benson et al., Reference Benson, Cavanaugh, Clark, Karsch-Mizrachi, Lipman, Ostell and Sayers2012) representative of the major piroplasm clades (Jalovecka et al., Reference Jalovecka, Sojka, Ascencio and Schnittger2019) and the alignments were trimmed to the same size (557 bp). Maximum likelihood phylogenetic trees based on the 18S rRNA and cytb loci were constructed using IQ-TREE with 1000 ultrafast bootstrap replicates (Nguyen et al., Reference Nguyen, Schmidt, von Haeseler and Minh2015; Hoang et al., Reference Hoang, Chernomor, von Haeseler, Minh and Vinh2018). ModelFinder was used to choose the best-fitting substitution model based on Bayesian information criterion (Kalyaanamoorthy et al., Reference Kalyaanamoorthy, Minh, Wong, von Haeseler and Jermiin2017). Cardiosporidium cionae (GenBank accession number EU052685) and B. microti (GenBank accession number NC034637) were used as outgroups for the 18S rRNA and ctyb phylogenetic trees, respectively.

Statistical analysis

To determine whether there are statistically significant differences in infection rates in nymphal vs female H. longicornis, a binomial generalized linear model using the glm function in R (R Core Team, 2022) was used. To determine the correlation between sample size and infection rates, the lm function in R was used.

Bloodmeal analysis of partially engorged specimens

DNA from 2 partially engorged H. longicornis collected on 15 July 2021 (1 nymph and 1 adult) was isolated using DNeasy Blood and Tissue columns (Qiagen, Valencia, CA, USA). An extraction control was included, and all work was performed in a dedicated clean lab inside a laminar flow hood (Mystaire, Creedmor, NC, USA). Primers CytbVertR1 (5'-GGACGAGGACTATACTACGG-3' from Egizi et al., Reference Egizi, Healy and Fonseca2013) and BMF1 (5'-AAACTGCAGCCCCTCAGAATGATATTTGTCCTCA-3'), originally called H15149 (Kocher et al., Reference Kocher, Thomas, Meyer, Edwards, Paabo, Villablancatt and Wilson1989), were used to amplify a 132-nucleotide fragment in the cytochrome oxidase b locus using an annealing temperature of 55°C. A PCR product obtained from the adult H. longicornis was purified (ExoSAP-IT, Affymetrix, Santa Clara, CA, USA) then sequenced at Azenta Genewiz (South Plainfield). The sequences were trimmed and aligned with Geneious Prime 2023.0.1 (Biomatters Inc.) and the consensus was used as a query in NCBI's Basic Local Alignment Search Tool, BLASTn (Altschup et al., Reference Altschup, Gish, Miller, Myers and Lipman1990).

Results

Overall, 667 H. longicornis nymphs and adults collected from the environment were screened and evidence of piroplasm parasites was found in 18 ticks (2.7%, Table 2). Adult and nymph infection rates did not differ (infection rates of 2.4 and 2.8%, respectively, P value = 0.60). Infection rates among sampling sites reflected sample size (r 2 = 0.92, P value <0.01; using adult and nymph data separately to increase the statistical power).

Table 2. Numbers of H. longicornis nymphs and adults collected from the environment at 3 sites at Rutgers Cook campus and tested for piroplasm DNA (# Pos represents the number of ticks that were positive). Collections were started on 24 June 2021 and proceeded approximately bi-weekly until 10 September 2021 (refer to Ferreira et al., Reference Ferreira, González, Milholland, Tung and Fonseca2023 for details).

1 Babesia microti S837 (5 positive ticks)

2 Babesia sp. Coco (2 positive ticks)

3 Babesia sp. (4 positive ticks)

4 Theileria cervi (a total of 6 positive ticks)

5 Theileria sp. (1 positive ticks)

There was no evidence of piroplasm coinfections (such as double chromatogram peaks) in the positive ticks. The primers targeting the 18S rRNA gene amplified fragments ranging in size from 395 to 515 base pairs (bp) spanning the V4 hypervariable region (Cauvin et al., Reference Cauvin, Hood, Shuman, Orange, Blackburn, Sayler and Wisely2019).

The 18S rRNA fragment from a nymph collected at University Inn had a 100% pairwise identity to 5 Theileria sp. sequences (GenBank accession numbers MW008536, MW008531, MK262962, MK262963 and MK262959) that are considered ‘Type X’ or ‘divergent’ (Cauvin et al., Reference Cauvin, Hood, Shuman, Orange, Blackburn, Sayler and Wisely2019; Olafson et al., Reference Olafson, Buckmeier, May and Thomas2020). This sequence differed by 20 bp (17 mismatches and 3 deletions) from 6 other sequences obtained from H. longicornis also collected from the University Inn site. The closest match for these 6 sequences (99.8–100% pairwise identity) was a Theileria cervi type F sequence (GenBank U97054). Three of the 6 sequences, 2 from nymphs and 1 from an adult, were identical to U97054 while the remaining 3, all from nymphs, differed by 1 bp (GenBank accession numbers OR612075, OR612078, OR612080). All sequences clustered within a clade containing T. cervi and T. orientalis (Fig. 1a).

Figure 1. Phylogenetic trees with piroplasm sequences obtained from Haemaphysalis longicornis (in red). (a) Tree based on 18S rRNA locus. Constructed with TIM2 + F + I + G4 substitution model. (b) Tree based on cytochrome b locus. Constructed with K3Pu + F + I + G4 substitution model.

Two 18S rRNA sequences, one from a Goat Farm adult and the other from a University Inn nymph, were identical to a Babesia sp. isolate NYT-435 (GenBank accession number MW665118) recently sequenced from a pool of lone star ticks (Ambylomma americanum) from Staten Island, NYC (Jain et al., Reference Jain, Tagliafierro, Marques, Sanchez-Vicente, Gokden, Fallon, Mishra, Briese, Kapoor, Sameroff, Guo, Marcos, Hu, Lipkin and Tokarz2021). The sequence of the Babesia sp. isolate NYT-435 was also the closest match for 2 other sequences, 1 from another Goat Farm adult and a second University Inn nymph (99.7 and 99.5% pairwise identities with a 1 and 2 bp difference, respectively; GenBank accession numbers OR612068, OR612072, respectively). The next best match for these 4 sequences (97.9–98.2% pairwise identity) was a Babesia sp. from a white-tailed deer in Texas (GenBank accession number HQ264120). These 4 sequences clustered within the Babesia sensu stricto clade with B. bovis and Babesia ovis (Fig. 1a). From one of the specimens that were positive for the unknown Babesia, a 1009 bp DNA fragment was amplified with the cytb primers (GenBank accession number OR610156). The closest match to this sequence was Babesia motasi isolate Lintan (GenBank acc. num. MN605889) with 90.8% pairwise identity. The sequence from this H. longicornis clustered with B. motasi in the Babesia sensu stricto clade (Fig. 1b).

The 18S rRNA from 2 sequences, one from a University Inn nymph and another from a Goat Farm adult, matched Babesia sp. Coco. The sequence from the Goat Farm adult was 100% pairwise identical to Genbank EU109716, whereas the sequence from the University Inn nymph differed by 1 bp and had a 99.8% pairwise identity (GenBank accession number OR612082).

Finally, 5 18S rRNA sequences, from 2 nymphs and 1 adult from the Goat Farm and 2 nymphs from the University Inn, had a 100% pairwise identity to B. microti isolate S837 (GenBank accession number AY144698). In the phylogenetic tree, these sequences clustered with B. vulpes within the B. microti-like group (Fig. 1a; Jalovecka et al., Reference Jalovecka, Sojka, Ascencio and Schnittger2019).

A 132 bp DNA fragment was amplified and sequenced from a partially engorged adult H. longicornis and was 100% identical to a cytb fragment from white-tailed deer, Odocoileus virginianus (GenBank accession number AF535863).

Discussion

DNA sequences belonging to 3 piroplasm clades were detected in 2.7% of field-collected H. longicornis nymphs and adults in NJ, United States. Infection rates for Theileria, Babesia sensu stricto and B. microti, were 1, 0.9 and 0.8%, respectively. While these may be low piroplasm infection rates compared to studies that detected T. orientalis in 12.7% of H. longicornis in Virginia (Thompson et al., Reference Thompson, White, Shaw, Egizi, Lahmers, Ruder and Yabsley2020, Reference Thompson, White, Doub, Sharma, Frierson, Dominguez, Shaw, Weaver, Vigil, Bonilla, Ruder and Yabsley2022), those were collected from the cattle farm where the first US outbreak of T. orientalis Ikeda occurred, which would have increased the likelihood that local vectors were infected. Overall, there have been few exploratory examinations of piroplasms in H. longicornis in the United States.

This is the first report of Theileria species besides T. orientalis Ikeda in H. longicornis in the United States. Most sequences matched T. cervi type F, a piroplasm that commonly infects deer and other cervids (Cauvin et al., Reference Cauvin, Hood, Shuman, Orange, Blackburn, Sayler and Wisely2019; Olafson et al., Reference Olafson, Buckmeier, May and Thomas2020). However, 1 sequence differed in at least 20 bp from T. cervi type F and instead matched a strain denoted ‘divergent’ or ‘type X’ found in wild and farmed deer in Florida and in Anocenter nitens, the tropical horse tick, parasitizing white-tailed deer in Texas (Cauvin et al., Reference Cauvin, Hood, Shuman, Orange, Blackburn, Sayler and Wisely2019; Olafson et al., Reference Olafson, Buckmeier, May and Thomas2020). US populations of H. longicornis have often been reported feeding on white-tailed deer (Tufts et al., Reference Tufts, Goodman, Benedict, Davis, VanAcker and Diuk-Wasser2021), and finding deer DNA in a partially engorged tick supports this. These findings indicate that H. longicornis may be involved in the transmission cycle of Theileria in NJ.

This is also the first report of Babesia in questing un-engorged H. longicornis in the United States. Specifically, the phylogenetic analyses using both 18S and cytochrome b loci indicate the unknown Babesia sequences fall within the Babesia sensu stricto clade. Babesia sensu stricto (also referred to as ‘true Babesia’) are pathogens of both veterinary and medical importance capable of transovarial transmission in their tick vectors, allowing the parasite to propagate in the absence of vertebrate reservoirs (Jalovecka et al., Reference Jalovecka, Sojka, Ascencio and Schnittger2019).

Babesia sp. Coco, another ‘true Babesia’ was also detected in 2 H. longicornis. Babesia sp. Coco can be pathogenic to dogs but usually only if they are immunocompromised (Birkenheuer et al., Reference Birkenheuer, Neel, Ruslander, Levy and Breitschwerdt2004; Holman et al., Reference Holman, Backlund, Wilcox, Stone, Stricklin and Bardin2009; Sikorski et al., Reference Sikorski, Birkenheuer, Holowaychuk, McCleary-Wheeler, Davis and Littman2010; Dear and Birkenheuer, Reference Dear and Birkenheuer2022), so it is unclear whether dogs are incidental hosts or are an important reservoir species that manifest clinical disease when immunocompromised. Dogs are considered an important blood host for H. longicornis in the United States (Trout Fryxell et al., Reference Trout Fryxell, Vann, Butler, Paulsen, Chandler, Willis, Wyrosdick, Schaefer, Gerhold, Grove, Ivey, Thompson, Applegate, Sweaney, Daniels, Beaty, Balthaser, Freye, Mertins, Bonilla and Lahmers2021; Thompson et al., Reference Thompson, White, Doub, Sharma, Frierson, Dominguez, Shaw, Weaver, Vigil, Bonilla, Ruder and Yabsley2022).

Finally, this is the first report of B. microti genotype S837 in H. longicornis. This genotype is found in skunks Mephitis mephitis (Goethert, Reference Goethert2021), which are an important host for this tick species in NJ (Ferreira et al., Reference Ferreira, González, Milholland, Tung and Fonseca2023) and, critically, H. longicornis has been shown to be a competent vector of B. microti under experimental conditions (Wu et al., Reference Wu, Cao, Zhou, Zhang, Gong and Zhou2017). Of note, there is an overall lack of knowledge of the biology and epidemiology of wildlife piroplasms in the northeastern USA. In areas endemic for human babesiosis, molecular studies of wildlife piroplasms rarely employ sequencing to confirm parasite identity, although piroplasm parasites thought to only infect wildlife and/or domestic animals have recently been reported infecting humans in North America (Scott et al., Reference Scott, Sajid, Pascoe and Foley2021) and elsewhere (Hong et al., Reference Hong, Kim, Song, Roh, Cho, Kim, Um, Kwak, Cho and Lee2019).

Although the health risk of H. longicornis to US livestock was established by the outbreak of T. orientalis Ikeda in cattle in the state of Virginia (Thompson et al., Reference Thompson, White, Shaw, Egizi, Lahmers, Ruder and Yabsley2020; Dinkel et al., Reference Dinkel, Herndon, Noh, Lahmers, Todd, Ueti, Scoles, Mason and Fry2021), there is limited research regarding the veterinary and medical significance of H. longicornis in the USA, which has focused primarily on testing the ability of US specimens to transmit pathogens of known public health concern (Breuner et al., Reference Breuner, Ford, Hojgaard, Osikowicz, Parise, Rosales Rizzo, Bai, Levin, Eisen and Eisen2020; Stanley et al., Reference Stanley, Ford, Snellgrove, Hartzer, Smith, Krapiunaya and Levin2020; Raney et al., Reference Raney, Perry and Hermance2022a). The discovery of various piroplasm parasites in questing H. longicornis specimens in NJ brings forth new concerns. As a result, it becomes crucial to prioritize studies that delve into the role of H. longicornis as an actual vector for these pathogens.

Data availability statement

Nucleotide sequence data reported in this paper are available in GenBank™, EMBL and DDBJ databases under the accession numbers OR610156 and OR612065-OR612082.

Acknowledgements

We thank Nicole Wagner MS from the Rutgers Center for Vector Biology and Julia Brennan, now developing her PhD at Virginia Tech, for their technical support. We also thank Dr Glen Scoles, USDA-ARS Beltsville, for sharing his piroplasm expertise.

Author contributions

Heidi E. Herb: conceptualization, methodology, formal analysis, investigation, data curation, writing – original draft, writing – review and editing, visualization. Julia González and Francisco C. Ferreira: investigation, data curation, writing – review & editing. Dina M. Fonseca: conceptualization, methodology, resources, writing – review & editing, supervision, project administration, funding acquisition.

Financial support

This work was supported in part by a United States Department of Agriculture National Institute of Food and Agriculture Multistate Grant, NE1943, to D. M. F., a State of New Jersey FY22 Tick Research and Control-Special Purpose Funding to CVB/NJAES and by Cooperative Agreement Number 1U01CK000509 between the Centers for Disease Control and Prevention (CDC) and the Northeast Center of Excellence in Vector Borne Diseases (PI at Rutgers, D. M. F). The contents of this work are solely the responsibility of the authors and do not necessarily represent the official views of the Centers for Disease Control and Prevention or the Department of Health and Human Services.

Competing interests

None.

Ethical standards

Not applicable.

Footnotes

*

Current address: Department of Entomology, Texas A&M University, College Station, TX 77843, USA.

A similar manuscript has been uploaded to biorxiv – https://www.biorxiv.org/content/10.1101/2023.06.26.546593v1

References

Almazán, C, Scimeca, RC, Reichard, MV and Mosqueda, J (2022) Babesiosis and theileriosis in North America. Pathogens 11, 168.CrossRefGoogle ScholarPubMed
Altschup, SF, Gish, W, Miller, W, Myers, EW and Lipman, DJ (1990) Basic local alignment search tool. Journal of Molecular Biology 215, 403–410.CrossRefGoogle Scholar
Athni, TS, Shocket, MS, Couper, LI, Nova, N, Caldwell, IR, Caldwell, JM, Childress, JN, Childs, ML, De Leo, GA, Kirk, DG, MacDonald, AJ, Olivarius, K, Pickel, DG, Roberts, SO, Winokur, OC, Young, HS, Cheng, J, Grant, EA, Kurzner, PM, Kyaw, S, Lin, BJ, Lopez, RC, Massihpour, DS, Olsen, EC, Roache, M, Ruiz, A, Schultz, EA, Shafat, M, Spencer, RL, Bharti, N and Mordecai, EA (2021) The influence of vector-borne disease on human history: socio-ecological mechanisms. Ecology Letters 24, 829–846.CrossRefGoogle ScholarPubMed
Bautista, JLR, Ikadai, H, You, M, Battsetseg, B, Igarashi, I, Nagasawa, H and Fujisaki, K (2001) Molecular evidence of Babesia caballi (Nuttall and Strickland, 1910) parasite transmission from experimentally-infected SCID mice to the ixodid tick, Haemaphysalis longicornis (Neuman, 1901). Veterinary Parasitology 102, 185–191.CrossRefGoogle Scholar
Benson, DA, Cavanaugh, M, Clark, K, Karsch-Mizrachi, I, Lipman, DJ, Ostell, J and Sayers, EW (2012) GenBank. Nucleic Acids Research 41, D36–D42.CrossRefGoogle ScholarPubMed
Bickerton, M and Toledo, A (2020) Multiple pruritic tick bites by Asian longhorned tick larvae (Haemaphysalis longicornis). International Journal of Acarology 46, 373–376.CrossRefGoogle Scholar
Bickerton, M, McSorley, K and Toledo, A (2021) A life stage-targeted acaricide application approach for the control of Haemaphysalis longicornis. Ticks and Tick-borne Diseases 12, 101581.CrossRefGoogle ScholarPubMed
Birkenheuer, AJ, Neel, J, Ruslander, D, Levy, MG and Breitschwerdt, EB (2004) Detection and molecular characterization of a novel large Babesia species in a dog. Veterinary Parasitology 124, 151–160.CrossRefGoogle Scholar
Bloch, EM, Kumar, S and Krause, PJ (2019) Persistence of Babesia microti infection in humans. Pathogens 8, 102.CrossRefGoogle ScholarPubMed
Breden, TF, Alger, Y, Strakosch Walz, K and Windisch, A (2001) Classification of vegetation communities of New Jersey: second iteration. Office of Natural Lands Management, Division of Parks and Forestry, NJDEP.Google Scholar
Breuner, NE, Ford, SL, Hojgaard, A, Osikowicz, LM, Parise, CM, Rosales Rizzo, MF, Bai, Y, Levin, ML, Eisen, RJ and Eisen, L (2020) Failure of the Asian longhorned tick, Haemaphysalis longicornis, to serve as an experimental vector of the Lyme disease spirochete, Borrelia burgdorferi sensu stricto. Ticks and Tick-borne Diseases 11, 101311.CrossRefGoogle ScholarPubMed
Casati, S, Sager, H, Gern, L and Piffaretti, J-C (2006) Presence of potentially pathogenic Babesia sp. for human in Ixodes ricinus in Switzerland. Annals of Agricultural and Environmental Medicine: AAEM 13, 65–70.Google Scholar
Cauvin, A, Hood, K, Shuman, R, Orange, J, Blackburn, JK, Sayler, KA and Wisely, SM (2019) The impact of vector control on the prevalence of Theileria cervi in farmed Florida white-tailed deer, Odocoileus virginianus. Parasites & Vectors 12, 100.CrossRefGoogle ScholarPubMed
Cumbie, AN, Trimble, RN and Eastwood, G (2022) Pathogen spillover to an invasive tick species: first detection of bourbon virus in Haemaphysalis longicornis in the United States. Pathogens 11, 454.CrossRefGoogle Scholar
Dear, JD and Birkenheuer, A (2022) Babesia in North America. Veterinary Clinics of North America: Small Animal Practice 52, 1193–1209.CrossRefGoogle ScholarPubMed
Dinkel, KD, Herndon, DR, Noh, SM, Lahmers, KK, Todd, SM, Ueti, MW, Scoles, GA, Mason, KL and Fry, LM (2021) A U.S. isolate of Theileria orientalis, Ikeda genotype, is transmitted to cattle by the invasive Asian longhorned tick, Haemaphysalis longicornis. Parasites & Vectors 14, 157.CrossRefGoogle Scholar
Egizi, A, Healy, SP and Fonseca, DM (2013) Rapid blood meal scoring in anthropophilic Aedes albopictus and application of PCR blocking to avoid pseudogenes. Infection, Genetics and Evolution 16, 122–128.CrossRefGoogle ScholarPubMed
Egizi, AM, Robbins, RG, Beati, L, Nava, S, Evans, CR, Occi, JL and Fonseca, DM (2019a) A pictorial key to differentiate the recently detected exotic Haemaphysalis longicornis Neumann, 1901 (Acari, Ixodidae) from native congeners in North America. ZooKeys 818, 117–128.CrossRefGoogle Scholar
Egizi, AM, Occi, JL, Price, DC and Fonseca, DM (2019b) Leveraging the expertise of the New Jersey mosquito control community to jump start standardized tick surveillance. Insects 10, 219.CrossRefGoogle ScholarPubMed
Ferreira, FC, González, J, Milholland, MT, Tung, GA and Fonseca, DM (2023) Ticks (Acari: Ixodida) on synanthropic small and medium-sized mammals in areas of the northeastern United States infested with the Asian longhorned tick, Haemaphysalis longicornis. International Journal for Parasitology. doi: 10.1016/j.ijpara.2023.06.003. Online ahead of print.CrossRefGoogle ScholarPubMed
Garrett, KB, Hernandez, SM, Balsamo, G, Barron, H, Beasley, JC, Brown, JD, Cloherty, E, Farid, H, Gabriel, M, Groves, B, Hamer, S, Hill, J, Lewis, M, McManners, K, Nemeth, N, Oesterle, P, Ortiz, S, Peshock, L, Schnellbacher, R, Schott, R, Straif-Bourgeois, S and Yabsley, MJ (2019) Prevalence, distribution, and diversity of cryptic piroplasm infections in raccoons from selected areas of the United States and Canada. International Journal for Parasitology: Parasites and Wildlife 9, 224–233.Google ScholarPubMed
Goethert, HK (2021) What Babesia microti is now. Pathogens 10, 1168.CrossRefGoogle ScholarPubMed
González, J, Fonseca, DM and Toledo, A (2023) Seasonal dynamics of tick species in the ecotone of parks and recreational areas in Middlesex County (New Jersey, USA). Insects 14, 258.CrossRefGoogle ScholarPubMed
Gray, JS, Estrada-Peña, A and Zintl, A (2019) Vectors of babesiosis. Annual Review of Entomology 64, 149–165.CrossRefGoogle ScholarPubMed
Heath, A (2016) Biology, ecology and distribution of the tick, Haemaphysalis longicornis Neumann (Acari: Ixodidae) in New Zealand. New Zealand Veterinary Journal 64, 10–20.CrossRefGoogle ScholarPubMed
Hoang, DT, Chernomor, O, von Haeseler, A, Minh, BQ and Vinh, LS (2018) UFBoot2: improving the ultrafast bootstrap approximation. Molecular Biology and Evolution 35, 518–522.CrossRefGoogle ScholarPubMed
Holman, PJ, Backlund, BB, Wilcox, AL, Stone, R, Stricklin, AL and Bardin, KE (2009) Detection of a large unnamed Babesia piroplasm originally identified in dogs in North Carolina in a dog with no history of travel to that state. Journal of the American Veterinary Medical Association 235, 851–854.CrossRefGoogle Scholar
Hong, S-H, Kim, S-Y, Song, BG, Roh, JY, Cho, CR, Kim, C-N, Um, T-H, Kwak, YG, Cho, S-H and Lee, S-E (2019) Detection and characterization of an emerging type of Babesia sp. similar to Babesia motasi for the first case of human babesiosis and ticks in Korea. Emerging Microbes & Infections 8, 869–878.CrossRefGoogle ScholarPubMed
Horowitz, R and Freeman, PR (2020) Archives of medical case reports case report healthy fetal outcomes using a novel treatment for maternal Lyme disease and babesiosis during consecutive pregnancies: a case study and literature review. Archives of Medical Case Reports 2, 1–19.Google Scholar
Jain, K, Tagliafierro, T, Marques, A, Sanchez-Vicente, S, Gokden, A, Fallon, B, Mishra, N, Briese, T, Kapoor, V, Sameroff, S, Guo, C, Marcos, LA, Hu, L, Lipkin, WI and Tokarz, R (2021) Development of a capture sequencing assay for enhanced detection and genotyping of tick-borne pathogens. Scientific Reports 11, 12384.CrossRefGoogle ScholarPubMed
Jalovecka, M, Sojka, D, Ascencio, M and Schnittger, L (2019) Babesia life cycle – when phylogeny meets biology. Trends in Parasitology 35, 356–368.CrossRefGoogle ScholarPubMed
Kalyaanamoorthy, S, Minh, BQ, Wong, TKF, von Haeseler, A and Jermiin, LS (2017) ModelFinder: fast model selection for accurate phylogenetic estimates. Nature Methods 14, 587–589.CrossRefGoogle ScholarPubMed
Keirans, JE and Litwak, TR (1989) Pictorial key to the adults of hard ticks, family ixodidae (Ixodida: Ixodoidea), East of the Mississippi River. Journal of Medical Entomology 26, 435–448.CrossRefGoogle Scholar
Kocher, TD, Thomas, WK, Meyer, A, Edwards, SV, Paabo, S, Villablancatt, FX and Wilson, AC (1989) Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proceedings of the National Academy of Sciences 86, 6196–6200.CrossRefGoogle ScholarPubMed
Kumar, A, O'Bryan, J and Krause, PJ (2021) The global emergence of human babesiosis. Pathogens 10. doi: 10.3390/pathogens10111447CrossRefGoogle ScholarPubMed
Lerner, H and Berg, C (2015) The concept of health in One Health and some practical implications for research and education: what is One Health? Infection Ecology & Epidemiology 5, 25300.CrossRefGoogle ScholarPubMed
Li, Y, Luo, J, Guan, G, Ma, M, Liu, A, Liu, J, Ren, Q, Niu, Q, Lu, B, Gao, J, Liu, Z, Dang, Z, Tian, Z, Zhang, B, He, Z, Bai, Q and Yin, H (2009) Experimental transmission of Theileria uilenbergi infective for small ruminants by Haemaphysalis longicornis and Haemaphysalis qinghaiensis. Parasitology Research 104, 1227–1231.CrossRefGoogle ScholarPubMed
Li, A, Liu, L, Wu, W, Liu, Y, Huang, X, Li, C, Liu, D, Li, J, Wang, S, Li, D and Liang, M (2021) Molecular evolution and genetic diversity analysis of SFTS virus based on next-generation sequencing. Biosafety and Health 3, 105–115.CrossRefGoogle Scholar
Liu, K, Zhou, H, Sun, R-X, Yao, H-W, Li, Y, Wang, L-P, Mu, D, Li, X-L, Yang, Y, Gray, GC, Cui, N, Yin, W-W, Fang, L-Q, Yu, H-J and Cao, W-C (2015) A national assessment of the epidemiology of severe fever with thrombocytopenia syndrome, China. Scientific Reports 5, 9679.CrossRefGoogle ScholarPubMed
Luo, L-M, Zhao, L, Wen, H-L, Zhang, Z-T, Liu, J-W, Fang, L-Z, Xue, Z-F, Ma, D-Q, Zhang, X-S, Ding, S-J, Lei, X-Y and Yu, X (2015) Haemaphysalis longicornis ticks as reservoir and vector of severe fever with thrombocytopenia syndrome virus in China. Emerging Infectious Diseases 21, 1770–1776.CrossRefGoogle ScholarPubMed
Marendy, D, Baker, K, Emery, D, Rolls, P and Stutchbury, R (2020) Haemaphysalis longicornis: the life-cycle on dogs and cattle, with confirmation of its vector status for Theileria orientalis in Australia. Veterinary Parasitology 277, 100022.CrossRefGoogle Scholar
Nguyen, L-T, Schmidt, HA, von Haeseler, A and Minh, BQ (2015) IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Molecular Biology and Evolution 32, 268–274.CrossRefGoogle ScholarPubMed
Olafson, PU, Buckmeier, BG, May, MA and Thomas, DB (2020) Molecular screening for rickettsial bacteria and piroplasms in ixodid ticks surveyed from white-tailed deer (Odocoileus virginianus) and nilgai antelope (Boselaphus tragocamelus) in southern Texas. International Journal for Parasitology: Parasites and Wildlife 13, 252–260.Google Scholar
Onyiche, TE, Răileanu, C, Fischer, S and Silaghi, C (2021) Global distribution of Babesia species in questing ticks: a systematic review and meta-analysis based on published literature. Pathogens 10, 230.CrossRefGoogle ScholarPubMed
Osbrink, WLA, Thomas, DB, Lohmeyer, KH and Temeyer, KB (2022) Climate change and alternative hosts complicate the eradication of cattle fever ticks (Acari: Ixodidae) in the Southern United States, a review. Annals of the Entomological Society of America 115, 39–55.CrossRefGoogle Scholar
Price, KJ, Witmier, BJ, Eckert, RA and Boyer, CN (2022) Recovery of partially engorged Haemaphysalis longicornis (Acari: Ixodidae) ticks from active surveillance. Journal of Medical Entomology 59, 1842–1846.CrossRefGoogle ScholarPubMed
Rainey, T, Occi, JL, Robbins, RG and Egizi, A (2018) Discovery of Haemaphysalis longicornis (Ixodida: Ixodidae) parasitizing a sheep in New Jersey, United States. Journal of Medical Entomology 55, 757–759.CrossRefGoogle ScholarPubMed
Rajapakshage, BK, Yamasaki, M, Hwang, S-J, Sasaki, N, Murakami, M, Tamura, Y, Lim, SY, Nakamura, K, Ohta, H and Takiguchi, M (2012) Involvement of mitochondrial genes of Babesia gibsoni in resistance to diminazene aceturate. Journal of Veterinary Medical Science 74, 1139–1148.CrossRefGoogle Scholar
Raney, WR, Perry, JB and Hermance, ME (2022a) Transovarial transmission of heartland virus by invasive Asian longhorned ticks under laboratory conditions. Emerging Infectious Diseases 28, 726–729.CrossRefGoogle ScholarPubMed
Raney, WR, Herslebs, EJ, Langohr, IM, Stone, MC and Hermance, ME (2022b) Horizontal and vertical transmission of Powassan virus by the invasive Asian longhorned tick, Haemaphysalis longicornis, under laboratory conditions. Frontiers in Cellular and Infection Microbiology 12. doi: 10.3389/fcimb.2022.923914CrossRefGoogle ScholarPubMed
Rankins, EM and Malinowski, K (2020) Horse racing and veterinary practices in New Jersey. Journal of Equine Veterinary Science 85, 102879.CrossRefGoogle ScholarPubMed
R Core Team (2022) R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. Available at https://www.R-project.org/Google Scholar
Renard, I and Mamoun, CB (2021) Treatment of human babesiosis: then and now. Pathogens 10, 1120.CrossRefGoogle ScholarPubMed
Rochlin, I, Egizi, A, Narvaez, Z, Bonilla, DL, Gallagher, M, Williams, GM, Rainey, T, Price, DC and Fonseca, DM (2023) Microhabitat modeling of the invasive Asian longhorned tick (Haemaphysalis longicornis) in New Jersey, USA. Ticks and Tick-borne Diseases 14, 102126.CrossRefGoogle ScholarPubMed
Schappach, BL, Krell, RK, Hornbostel, VL and Connally, NP (2020) Exotic Haemaphysalis longicornis (Acari: Ixodidae) in the United States: biology, ecology, and strategies for management. Journal of Integrated Pest Management 11. doi: 10.1093/jipm/pmaa019CrossRefGoogle Scholar
Schnittger, L, Ganzinelli, S, Bhoora, R, Omondi, D, Nijhof, AM and Florin-Christensen, M (2022) The Piroplasmida Babesia, Cytauxzoon, and Theileria in farm and companion animals: species compilation, molecular phylogeny, and evolutionary insights. Parasitology Research 121, 1207–1245.CrossRefGoogle ScholarPubMed
Schreeg, ME, Marr, HS, Tarigo, JL, Cohn, LA, Bird, DM, Scholl, EH, Levy, MG, Wiegmann, BM and Birkenheuer, AJ (2016) Mitochondrial genome sequences and structures aid in the resolution of piroplasmida phylogeny. PLoS ONE 11, e0165702.CrossRefGoogle ScholarPubMed
Scott, JD and Scott, CM (2018) Human babesiosis caused by Babesia duncani has widespread distribution across Canada. Healthcare 6. doi: 10.3390/healthcare6020049CrossRefGoogle ScholarPubMed
Scott, JD, Sajid, MS, Pascoe, EL and Foley, JE (2021) Detection of Babesia odocoilei in humans with babesiosis symptoms. Diagnostics 11, 947.CrossRefGoogle ScholarPubMed
Sikorski, LE, Birkenheuer, AJ, Holowaychuk, MK, McCleary-Wheeler, AL, Davis, JM and Littman, MP (2010) Babesiosis caused by a large Babesia species in 7 immunocompromised dogs. Journal of Veterinary Internal Medicine 24, 127–131.CrossRefGoogle ScholarPubMed
Stanley, HM, Ford, SL, Snellgrove, AN, Hartzer, K, Smith, EB, Krapiunaya, I and Levin, ML (2020) The ability of the invasive Asian longhorned tick Haemaphysalis longicornis (Acari: Ixodidae) to acquire and transmit Rickettsia rickettsii (Rickettsiales: Rickettsiaceae), the agent of Rocky Mountain spotted fever, under laboratory conditions. Journal of Medical Entomology 57, 1635–1639.CrossRefGoogle ScholarPubMed
Swanson, M, Pickrel, A, Williamson, J and Montgomery, S (2023) Trends in reported babesiosis cases – United States, 2011–2019. Morbidity and Mortality Weekly Report 72, 273–277.CrossRefGoogle ScholarPubMed
Thompson, AT, White, S, Shaw, D, Egizi, A, Lahmers, K, Ruder, MG and Yabsley, MJ (2020) Theileria orientalis Ikeda in host-seeking Haemaphysalis longicornis in Virginia, U.S.A. Ticks and Tick-borne Diseases 11, 101450.CrossRefGoogle ScholarPubMed
Thompson, AT, White, SA, Doub, EE, Sharma, P, Frierson, K, Dominguez, K, Shaw, D, Weaver, D, Vigil, SL, Bonilla, DL, Ruder, MG and Yabsley, MJ (2022) The wild life of ticks: using passive surveillance to determine the distribution and wildlife host range of ticks and the exotic Haemaphysalis longicornis, 2010–2021. Parasites & Vectors 15, 331.CrossRefGoogle ScholarPubMed
Trout Fryxell, RT, Vann, DN, Butler, RA, Paulsen, DJ, Chandler, JG, Willis, MP, Wyrosdick, HM, Schaefer, JJ, Gerhold, RW, Grove, DM, Ivey, JZ, Thompson, KW, Applegate, RD, Sweaney, J, Daniels, S, Beaty, S, Balthaser, D, Freye, JD, Mertins, JW, Bonilla, DL and Lahmers, K (2021) Rapid discovery and detection of Haemaphysalis longicornis through the use of passive surveillance and collaboration: building a state tick-surveillance network. International Journal of Environmental Research and Public Health 18, 7980.CrossRefGoogle ScholarPubMed
Tufts, DM, Goodman, LB, Benedict, MC, Davis, AD, VanAcker, MC and Diuk-Wasser, M (2021) Association of the invasive Haemaphysalis longicornis tick with vertebrate hosts, other native tick vectors, and tick-borne pathogens in New York City, USA. International Journal for Parasitology 51, 149–157.CrossRefGoogle ScholarPubMed
Votýpka, J, Modrý, D, Oborník, M, Šlapeta, J and Lukeš, J (2016) Apicomplexa. In Archibald JM, Simpson AGB, Slamovits CH, Margulis L, Melkonian M, Chapman DJ and Corliss JO (eds), Handbook of the Protists. Cham: Springer, pp 58. https://doi.org/10.1007/978-3-319-32669-6_20-1.Google Scholar
Wormser, GP, McKenna, D, Piedmonte, N, Vinci, V, Egizi, AM, Backenson, B and Falco, RC (2020) First recognized human bite in the United States by the Asian longhorned tick, Haemaphysalis longicornis. Clinical Infectious Diseases 70, 314–316.CrossRefGoogle ScholarPubMed
Wu, J, Cao, J, Zhou, Y, Zhang, H, Gong, H and Zhou, J (2017) Evaluation on infectivity of Babesia microti to domestic animals and ticks outside the Ixodes genus. Frontiers in Microbiology 8. doi: 10.3389/fmicb.2017.01915.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Primers used for piroplasm detection in H. longicornis

Figure 1

Table 2. Numbers of H. longicornis nymphs and adults collected from the environment at 3 sites at Rutgers Cook campus and tested for piroplasm DNA (# Pos represents the number of ticks that were positive). Collections were started on 24 June 2021 and proceeded approximately bi-weekly until 10 September 2021 (refer to Ferreira et al., 2023 for details).

Figure 2

Figure 1. Phylogenetic trees with piroplasm sequences obtained from Haemaphysalis longicornis (in red). (a) Tree based on 18S rRNA locus. Constructed with TIM2 + F + I + G4 substitution model. (b) Tree based on cytochrome b locus. Constructed with K3Pu + F + I + G4 substitution model.