Skip to main content
Download PDF

Eco-epidemiology of Rickettsia amblyommatis and Rickettsia parkeri in naturally infected ticks (Acari: Ixodida) from South Carolina

Parasites & Vectors volume 17, Article number: 33 (2024) Cite this article

Abstract

Background

Spotted fever group Rickettsia (SFGR) is the largest group of Rickettsia species of clinical and veterinary importance emerging worldwide. Historically, SFGR cases were linked to Rickettsia rickettsii, the causal agent of Rocky Mountain spotted fever; however, recently discovered species Rickettsia parkeri and Rickettsia amblyommatis have been shown to cause a wide range of clinical symptoms. The role of R. amblyommatis in SFGR eco-epidemiology and the possible public health implications remain unknown.

Methods

This study evaluated statewide tick surveillance and land-use classification data to define the eco-epidemiological relationships between R. amblyommatis and R. parkeri among questing and feeding ticks collected across South Carolina between 2021 and 2022. Questing ticks from state parks and feeding ticks from animal shelters were evaluated for R. parkeri and R. amblyommatis using reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) on pooled samples. A Bayesian multivariable logistic regression model for pool testing data was used to assess associations between R. parkeri or R. amblyommatis infection and land-use classification variables among questing ticks. The Spearman correlation was used to evaluate the relationship between the two tested pathogens.

Results

The infection prevalence for R. amblyommatis was 24.8% (23.4–26.3%) among questing ticks, and 39.5% (37.4–42.0%) among feeding ticks; conversely, for R. parkeri it was 19.0% (17.6–20.5%) among questing ticks and 22.4% (20.3–24.5%) among feeding ticks. A negative, refractory correlation was found between the species, with ticks significantly more likely to contain one or the other pathogen, but not both simultaneously. The Bayesian analysis revealed that R. amblyommatis infection was positively associated with deciduous, evergreen, and mixed forests, and negatively associated with hay and pasture fields, and emergent herbaceous wetlands. Rickettsia parkeri infection was positively associated with deciduous, mixed, and evergreen forests, herbaceous vegetation, cultivated cropland, woody wetlands, and emergent herbaceous wetlands, and negatively associated with hay and pasture fields.

Conclusions

This is the first study to evaluate the eco-epidemiological factors driving tick pathogenicity in South Carolina. The negative interactions between SFGR species suggest the possible inhibition between the two pathogens tested, which could have important public health implications. Moreover, land-use classification factors revealed environments associated with tick pathogenicity, highlighting the need for tick vector control in these areas.

Graphical Abstract

Background

Rickettsia spp. within the spotted fever group (SFGR) are a large group of tick-borne intracellular bacteria comprising several clinical and veterinary important species, and pose a major re-emerging public health concern worldwide [1]. SFGR pathogenicity varies from severe to mild disease: Rickettsia rickettsii, the causal agent of Rocky Mountain spotted fever (RMSF), is considered the most pathogenic species, while other species such as Rickettsia parkeri result in a similar, often less severe disease [2]. The main R. rickettsii vectors in the United States are Dermacentor variabilis, Dermacentor andersoni, and Rhipicephalus sanguineus [3, 4]. Rickettsia parkeri is also an important disease-causing agent, first described in 2004 among Amblyomma maculatum ticks [2]. The clinical symptoms associated with this species are considered different from those caused by R. rickettsii; however, unless molecular diagnostic testing is performed, the serological tests cannot be used to confirm species, and therefore clinical confirmation is based entirely on symptomology. One species in particular, Rickettsia amblyommatis (previously Candidatus R. amblyommii), has an undefined pathogenicity that has been highly debated by the clinical community. Some propose this species to be pathogenic, while others argue that it is non-pathogenic [5, 6]. Clinical diagnostic tests traditionally do not distinguish between SFGR species, and thus the ability to study veritable R. parkeri or other Rickettsia spp. clinical cases is challenging, contributing to the ongoing lack of knowledge regarding the pathogenic potential of other species.

Laboratory mice and guinea pig studies have demonstrated that R. amblyommatis can cause clinical symptoms under experimental conditions [7, 8], suggesting that mild human cases could be related to this species. Serological evidence in patients with a history of disease suggests that humans can become infected and produce high antibody titers against R. amblyommatis [9]. In South Carolina, USA, a pediatric hospitalization case was reported in the late 2010s (Chris Evans, personal communication, May 4, 2022) further suggesting the pathogenic potential of this species. We hypothesize that a portion of SFGR clinical cases linked to R. amblyommatis are misdiagnosed due to Rickettsia spp. antibody cross-reactivity, given the large number of RMSF cases in areas with high R. amblyommatis-infected tick prevalence [10]. For instance, R. rickettsii has been found in approximately 0.1% of ticks, while R. amblyommatis can be found present in more than 50% of ticks, indicating a greater R. amblyommatis exposure risk [11, 12]. Given the potential for R. amblyommatis to cause human disease, the need for targeted tick pathogen surveillance is warranted.

Rickettsia amblyommatis was first isolated from an Amblyomma americanum tick collected in Tennessee in 1974, and repeatedly detected in naturally infected A. americanum ticks throughout the USA [13, 14]. North Carolina is one of the states with the highest RMSF incidence [15]; however, since R. rickettsii is rarely found in ticks, investigators have postulated that some cases could be due to R. amblyommatis infection [9, 16, 17]. A tick surveillance initiative in North Carolina detected R. amblyommatis in 37% of A. americanum tick pools and 8.7% of D. variabilis tick pools but found R. rickettsii in 0% of pools, indicating that those bitten by ticks in North Carolina are far more likely to encounter R. amblyommatis when bitten [9]. Similarly, in another study, almost 55% of A. americanum ticks were R. amblyommatis-positive across all tested sites in Chatham County, NC [16]. An Oklahoma study detected R. amblyommatis in dogs exposed to environmental ticks (> 90% A. americanum, 3% A. maculatum, and 6% D. variabilis) [18]. Therefore, R. amblyommatis has been the most widely detected SFGR species across several studies investigating ticks attached to humans, although limited clinical signs have ever been described related to this species [17, 19]. In one case, a R. amblyommatis-positive tick was removed from the skin of a patient who developed a macular rash, although no other clinical signs were detected once the patient began doxycycline treatment [20]. Nevertheless, eco-epidemiological studies suggest that R. amblyommatis could be related to milder cases, since most SFGR cases happen in A. americanum high-prevalence areas (this tick is not considered a primary vector for R. rickettsii), and most patients with SFGR antibodies show higher titers against R. amblyommatis than any other Rickettsia spp. [11].

Despite the circumstantial evidence, little is known about the role of R. amblyommatis in SFGR cases and the eco-epidemiological factors associated with exposure. Under laboratory conditions, there have been several findings evaluating Rickettsia species interactions. One study found that R. amblyommatis-positive ticks co-feeding with R. parkeri-infected ticks are less likely to acquire the former, suggesting a possible inhibitory effect [21]. Other investigations have found that previous exposure to other Rickettsia spp. could reduce the severity of disease caused by R. rickettsii; however, no clear epidemiological evidence has been described in humans yet [22]. Thus, there is a need for epidemiological studies to unravel these pathogen interactions.

The South Carolina landscape promotes outdoor recreational activities at multiple state parks, creating tick exposure opportunities and subsequently SFGR exposure opportunities [23]. In 2020, a tick surveillance initiative was implemented in South Carolina where questing ticks were collected from state parks and public submissions, and host-attached ticks were collected from South Carolina animal shelters [24]. Given that South Carolina reports a moderate number of SFGR cases, and A. americanum is the most commonly found questing tick, R. amblyommatis is likely to be present in ticks collected throughout the state. These ticks should be evaluated to better understand the epidemiological characteristics of this infection [15, 24]. To evaluate the eco-epidemiological risk factors associated with R. amblyommatis in metastriate ticks, this study aims to investigate R. amblyommatis and R. parkeri distribution and prevalence utilizing reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) and land-use classification environmental factors using a Bayesian analysis across the state parks of South Carolina.

Methods

Study area and sampling locations

South Carolina is a coastal state in the Southeastern USA, with a total area of 82,933 km2, and is bordered by North Carolina to the north and northeast, Tennessee to the northwest, Georgia to the southwest, and the Atlantic Ocean to the southeast. Slightly over 5 million residents live in the state, with 45% living within the counties encompassing the major metroplexes: Charleston, Columbia, Greenville, and Myrtle Beach [25]. The state can be divided into three main areas, the Blue Ridge Mountain in the northwestern corner (rising to 1085 m elevation), the Piedmont from the mountains to the Sandhills in the Southeastern Plains (90 to 365 m elevation), and the Coastal Plain comprising Southern and Middle Atlantic regions (from sea level to 90 m elevation) delimited by the Sandhills and Coastal zone. The state has a subtropical climate, with hot, humid summers (average 23–33 °C, 15 days precipitation in July) and mild winters (average 3–15 °C, 9 days precipitation in January) [26]. Vegetation varies from woodlands in the Blue Ridge region, row crops in the Piedmont region with some loblolly pine forest, and gum, live oaks, cypresses, and magnolias across the coastal plain. Most wildlife is spread across the state, but some species follow geographical boundaries between the Coastal Plain and Piedmont [23]. Wildlife is best represented by white-tailed deer in the Piedmont and Coastal Plain, woodchucks and red squirrels found in the Blue Ridge, and American beaver, wild turkey, red foxes, and European wild pigs across the state.

For this study, we used metastriate ticks (Acari: Ixodidae) collected from January 2021 to December 2022 from a statewide tick surveillance effort recently implemented in South Carolina [24]. The ticks originated from either animal shelter submissions (host-attached ticks removed at the time of admission to one of 20 participating humane shelters statewide) or state park-based surveillance (collection of questing ticks using 10 CO2-baited traps, body-found ticks from the staff performing collections, and 50 m2 density dragging at each of 39 state parks) [27]. This program aims to elucidate tick-borne pathogen distribution across the state of South Carolina by collecting state park questing ticks, and tallying animal shelter feeding ticks.

Pathogen testing

Ticks were morphologically identified by species, sex, and life stage using dichotomous taxonomical keys by Keirans and Litwak [28], Egizi [29], and Clifford and Anastos [30]. Only metastriate ticks were selected and included, and non-metastriate ticks were kept in a 70% ethanol solution for future usage. After identification, host-attached ticks and questing ticks were bisected longitudinally and used for DNA extraction. Questing ticks were pooled by species, sex, life stage, collection method, collection date, and location in the following pool sizes: three for adults, five for nymphs, and a single pool for all the larvae. Host-attached ticks were tested individually. Each testing pool or individual tick was homogenized using 180 µl Qiagen Buffer ATL and two 5-mm stainless-steel beads using a bead homogenizer (TissueLyser, Qiagen, Germantown, MD, USA). DNA was extracted from the homogenized pools using QIAmp 96 DNA QIAcube Mini Kit (Qiagen, Germantown, MD, USA) on the QIAcube HT workstation (Qiagen, Germantown, MD, USA), following the manufacturer’s instructions. Gene amplification was performed to evaluate R. amblyommatis and R. parkeri DNA presence. DNA detection was performed by RT-qPCR using validated primers and probes for the outer membrane B gene ompB specific to R. amblyommatis and R. parkeri (Additional file 1: Table S1). The amplifications were performed using the QuantStudio 5 Real-Time PCR system (Applied Biosystems, Foster City, CA, USA) following the manufacturer’s instructions. Cycling conditions were 95 °C for 3 min and a two-step cycling for 40 cycles. Samples were considered positive for R. amblyommatis or R. parkeri when the cycle threshold (Ct) value was ≤ 40.

Statistical analysis

Descriptive statistics were performed to determine the pathogen distribution across tick species between questing and host-attached ticks, life stage, feeding status, and sampling method. Rickettsia amblyommatis and R. parkeri prevalence and 95% confidence intervals (CIs) were calculated for questing and attached ticks from PCR-positive pools using the Epitools® prevalence calculator for pooled samples (Ausvet, Fremantle, WA, Australia). This method accounts for different pool sizes to estimate prevalence [31]. A Spearman correlation was performed to evaluate the co-presence of R. amblyommatis and R. parkeri within pools. Additionally, the presence of R. amblyommatis and R. parkeri was mapped to evaluate the overlapping geographical distribution between the two pathogen species. For visualization, a layer including the ecological regions from South Carolina was obtained from the US Environmental Protection Agency [32].

Bayesian logistic regression models for group testing data [33] were used to evaluate the association between R. amblyommatis and R. parkeri presence and tick species, collection method, life stage, and land-use classification parameters. The Bayesian regression models were fitted for questing ticks collected in South Carolina state parks only due to the unknown geographical origin of the animals hosting attached ticks. The land-use classification variables were obtained from raster data from the National Land Cover Database from 2019 [34]. Values for percent land-use classification variables at each park were obtained to create a 10-km radius buffer within the park extension and summarized as a categorical raster to provide numerical values of the percentage area covered. The description of each land-use classification variable can be found in Additional file 1: Table S2. The values were exported to a data file and incorporated into the final dataset. All variables were evaluated using an unadjusted model and were considered to be included in the final adjusted model. Variables included in the final model were chosen using step-wise variable selection and model convergence. The models were run for 10,000 Markov chain Monte Carlo (MCMC) iterations with the first 5999 iterations discarded as burn-in. Convergence was assessed with trace plots (Additional file 2). A variable was considered to have a statistically important relationship with infection status if the 95% equal-tail credible interval for the variable’s regression coefficient did not contain 0. Statistical analysis and data visualization were performed using Epitools® (Ausvet, Fremantle, WA, Australia), RStudio, R version 4.1.3 (Free Software Foundation, Boston, MA, USA), and ArcGIS Pro 2.8.3 (Esri, Redlands, CA, USA).

Results

Between January 2021 and December 2022, a total of 4412 metastriate ticks were collected or submitted from South Carolina state parks (n = 2909) and animal shelters or humane societies (n = 1503), respectively. Of the total ticks collected, 2556 were adults, 1611 were nymphs, and 245 were larvae. There were five metastriate species collected: 3364 (76.2%) A. americanum, 799 (18.1%) D. variabilis, 195 (4.4%), A. maculatum, 44 (1.0%) R. sanguineus, and 10 (0.2%) Haemaphysalis longicornis. As shown in Table 1, only two species were obtained from parks—2188 (49.36%) were collected using CO2 traps, followed by 619 (14.0%) collected from density dragging. Animal shelters submitted ticks representing all five species. Pathogen point prevalence estimates among questing ticks—irrespective of species—were 24.8% (95% CI 23.4–26.3%) for R. amblyommatis and 19.0% (95% CI 17.6–20.5%) for R. parkeri. Among all feeding ticks, the estimated point prevalence was 39.5% (95% CI 37.4–42.0%) for R. amblyommatis and 22.4% (95% CI 20.3–24.5%) for R. parkeri.

Table 1 Classification of collected ticks by species, life stage, sampling method, and feeding status
Full size table

A significant moderate negative correlation between R. amblyommatis and R. parkeri positivity was found for shelter submissions (ρ = −0.142), for park collections (ρ = −0.422), for A. americanum (ρ = −0.301), and for R. sanguineus (ρ = −0.38) pools. Among life stages, nymphs showed the strongest negative correlation between R. amblyommatis and R. parkeri positivity (ρ = −0.388). Questing ticks showed the strongest negative correlation (ρ = −0.422) (Table 2).

Table 2 Rickettsia amblyommatis and R. parkeri positivity correlation by tick species, life stage, and location
Full size table

Ticks were collected from every region of the state, although the Blue Ridge region yielded the fewest ticks. Two state parks, Edisto State Park located in the Southern Coastal Plains and Sesquicentennial State Park located in the Sandhills, produced the largest volume of metastriate ticks. Similarly, two animal shelter locations, both located in the Piedmont region, produced the highest volume of host-attached ticks (Fig. 1). Some geographical overlapping existed between R. parkeri- and R. amblyommatis-infected ticks (Fig. 2); however, one or the other was typically predominant, and these were not present simultaneously in 50% or more of the pools from the same region. Rickettsia amblyommatis was generally collected throughout the southern part of the state, whereas R. parkeri-positive ticks were collected in the more northern part of the state.

Fig. 1
figure 1

Distribution and number of ticks collected statewide between January 2021 and December 2022

Full size image
Fig. 2
figure 2

Bivariate distribution of R. parkeri and R. amblyommatis metastriate tick positivity by host-seeking status. *Note: Questing ticks were collected at state parks, and host-attached ticks were collected from animal shelters and humane societies

Full size image

The Bayesian logistic regressions yielded statistically important variables that showed negative and positive associations with R. amblyommatis and R. parkeri infection (Table 3). In both models, female ticks were more likely to be positive for R. amblyommatis than males or immature ticks (odds ratio [OR]: 1.5, 95% CI 1.10–2.03; adjusted OR [aOR]: 1.43, 95% CI 1.03–1.99). Ticks collected through dragging were more likely to be positive for R. amblyommatis than ticks collected through CO2 traps or found on the body (OR: 1.53, 95% CI 1.11–2.06; aOR 1.55, 95% CI 1.11–2.15). In the unadjusted model, R. amblyommatis was more likely to be infecting ticks in the presence of deciduous forest (OR: 1.05, 95% CI 1.02–1.08), evergreen forest (OR: 1.06, 95%CI:1.03–1.08), and mixed forest (OR: 1.08, 95% CI 1.05–1.13). Both the unadjusted and adjusted models showed that R. amblyommatis infection was less likely to occur in the presence of emergent herbaceous wetlands (OR: 0.98, 95% CI 0.97–0.98; aOR: 0.90, 95% CI 0.83–0.98). In contrast, both unadjusted and adjusted models indicated that R. parkeri infection was positively associated with herbaceous vegetation (OR: 1.07, 95% CI 1.01–1.13; aOR: 1.61, 95% CI 1.22–2.15), cultivated crops (OR: 1.23, 95% CI 1.12–1.37; aOR: 19.93, 95% CI 4.28–69.12), and woody wetlands (OR: 1.03, 95% CI 1.01–1.05; aOR: 1.89, 95% CI 1.17–2.66). In the unadjusted model, R. parkeri was negatively associated with hay and pasture fields (OR: 0.86, 95% CI 0.78–0.95). In the adjusted model, only R. parkeri was positively associated with the presence of deciduous forest (aOR: 2.18, 95% CI1.51–2.95), mixed forest (aOR: 2.94, 95% CI 1.65–4.96), and emergent herbaceous wetlands (aOR: 2.01, 95% CI 1.51–2.61).

Table 3 Adjusted and unadjusted Bayesian logistic regression for pooled data evaluating R. amblyommatis and R. parkeri positivity prediction through environmental factors
Full size table

Discussion

This is the first study evaluating the ecological factors driving R. amblyommatis and R. parkeri infection among metastriate ticks in South Carolina and one of a few in the Southeastern USA. In this analysis, we estimated an overall R. amblyommatis prevalence of 24.8% among questing ticks and prevalence of 39.5% among host-attached ticks. In comparison, R. parkeri infection was slightly lower, with 19.0% prevalence among questing ticks and 22.4% among host-attached ticks. This analysis found a negative correlation between R. amblyommatis and R. parkeri presence within tested tick pools, potentially suggesting an antagonistic relationship between these two species. Ticks were less likely to be R. amblyommatis-positive in the presence of emergent herbaceous wetlands, whereas R. parkeri-positive ticks were more likely in herbaceous vegetation, cultivated crops, and woody wetlands. In the unadjusted models, both infections were less likely to be present in hay or pasture fields.

Environmental characteristics mediate the sustained transmission of tick-borne pathogens, and knowing the underlying variables associated with vector–pathogen presence can aid in vector control decision-making. It has been previously identified that increased suitable habitat presence correlates with greater species diversity, and thus enzootic cycles have more opportunity to propagate [35]. Considering that deciduous and mixed forests (environments home to greater species diversity) were both positively associated with R. amblyommatis and R. parkeri presence, sylvatic cycles, not evaluated in this analysis, likely serve as pathogen influencers in these tick populations [36]. In particular, hay pasture fields (environments with lower species diversity and potential acaricide use) [37] showed a negative association with R. amblyommatis and R. parkeri tick infection. On the contrary, cultivated cropland was positively associated with R. parkeri tick infection. Similarly, the scientific literature suggests that rural areas are associated with greater SFGR human risk [38, 39]; rural habitats typically sustain sylvatic transmission foci. Lastly, R. amblyommatis infections were less likely in emergent herbaceous wetlands. These landscapes are inhabited mostly by reptile, amphibian, and bird species, which are poorer propagators of tick-borne bacteria, potentially explaining the negative statistical association with R. amblyommatis in this study [40,41,42,43]. Future studies should aim for a One Health approach to evaluate interactions between hosts, environment, and ticks to better understand the pathogen distributions.

Rickettsia amblyommatis has been previously described in 17 countries and within 27 states in the USA [10]. First identified in A. americanum, the geographical distribution of this tick species in the USA overlaps with R. amblyommatis, further supported by high infection rates among this tick species [10, 44]. Amblyomma americanum is a known human-biting species, implicated in the transmission of several tick-borne pathogens [45]. From the tested pools in this study, R. amblyommatis was found in 61.8% of A. americanum, which is comparable to previously reported infection rates in questing ticks in SC and neighboring states ranging from 29 to 87% [46,47,48,49]. Given its wide distribution and high tick infection rates, R. amblyommatis could play an important role in SFGR epidemiology. In humans, seroconversion specific to R. amblyommatis has been observed, with some individuals expressing greater antibody titers to this than to other SFGR species, suggesting that humans might be exposed to R. amblyommatis at greater rates [50, 51]. Human pathogenicity from this species has been hypothesized and demonstrated in animal models, and anecdotal evidence suggests that R. amblyommatis infections could be implicated in undetermined febrile illness [9, 20]. Given the difficulty in determining species using an indirect fluorescent antibody (IFA) test due to cross-reactivity between Rickettsia spp., some authors suggest that probable mild RMSF cases could be caused by R. amblyommatis, and thus future studies should evaluate the pathogenicity of this species [9, 50, 51].

Questing ticks were collected from across the state; however, the majority of ticks were collected in the Southern Coastal Plains region compared to the Blue Ridge and the Piedmont regions. This could be explained by biodiversity differences across regions affecting tick populations. Moreover, CO2 traps, placed deeper in the forest compared with density dragging performed on the trails, are better suited for A. americanum, the predominant species collected [52]. Nonetheless, additional drivers such as anthropogenic activities and species diversity likely affect tick density in the Southern Coastal Plain region [53]. Conversely, greater numbers of host-attached ticks were collected in the Piedmont region. Additionally, more tick species were submitted that were host-attached than were collected in state parks, likely due to the low efficiency in capturing a diversity of species through CO2 traps and dragging. Interestingly, the collection methods were associated with R. amblyommatis positivity but not R. parkeri positivity. Ticks collected through dragging were more likely to be R. amblyommatis-positive ticks than those collected by CO2 traps; this suggests that R. amblyommatis infection impacts questing behavior, potentially affecting human SFGR exposure [54]. Despite this finding, the relationship between R. amblyommatis vector infection, questing behavior, and its epidemiological implications remains unknown.

A negative correlation was observed between R. amblyommatis and R. parkeri infections among A. americanum and R. sanguineus ticks. This negative correlation was seen for all life stages except larvae, suggesting the existence of mechanisms preventing co-infection between Rickettsia species. We hypothesize that R. amblyommatis infection acts as an inhibitory system for other Rickettsia spp. transmission, previously observed in laboratory conditions [21, 55, 56]. A study revealed that both Rickettsia rhipicephali and Rickettsia montanensis have the ability to inhibit R. parkeri during transovarial transmission in D. variabilis. It was reported that ticks infected with either species would not permit the transovarial transmission of the second species, indicating that ticks will not maintain both species equally, impacting the SFGR eco-epidemiology [55]. Although our results were not significant for negative correlation among larvae, a similar trend was observed (alpha = 0.1), which supports these findings. The current analysis supports the premise of naturally occurring inhibition of R. parkeri infection among R. amblyommatis-positive ticks during co-feeding. A laboratory evaluation confirmed that among A. americanum, R. parkeri infection was less likely to occur during co-feeding in R. amblyommatis-infected ticks than pathogen-free ticks [21]. In the analysis, this was not observed among A. maculatum, the most important R. parkeri vector. Unlike in A. americanum, the inhibitory phenomenon has not been observed in A. maculatum, and co-infections have been described, which could explain why A. maculatum is considered the predominant R. parkeri vector [56]. Interestingly, A. maculatum had the lowest R. amblyommatis positivity among host-attached ticks, and because co-feeding was not evaluated, the host-tick-pathogen interaction needs to be further evaluated. Despite R. amblyommatis being considered mildly pathogenic, these findings could reveal the mechanisms driving the inhibition of SFGR transmission and the public health implications of this species.

Some study limitations are worth mentioning. First, animal shelter-submitted ticks could not be evaluated for co-feeding as there was no information on which animals the ticks were coming from, nor which ticks were sharing a host. Therefore, the prevalence estimates, or the associations between pathogen presence and host-attached tick species could be overestimated. Second, due to the unknown location where the animal host became infested with ticks, the environmental land-use classification variables analysis was only performed for questing ticks, reducing the number of ticks included in the final analyses. Third, not all the collected ticks could be analyzed for R. amblyommatis or R. parkeri due to the low quantity of DNA extracted in some tick pools, therefore the correlation analysis was limited to high-quality DNA present specimens which are not representative of all populations. Finally, there was a disproportionate distribution of tick species despite systematic sampling methods used, with over 75% being A. americanum, which reduces the representation of some findings to other tick species. Despite these limitations, this is the first eco-epidemiological evaluation of R. amblyommatis and R. parkeri in South Carolina, and therefore these results remain meaningful ecological and epidemiological findings. Future directions should aim at including animal host densities, blood meal analysis, and possible tick species interactions within hosts, to better understand how the environmental impact on animal host and tick vector species can influence pathogen distribution.

Conclusions

In conclusion, R. amblyommatis and R. parkeri are widely distributed tick-borne pathogens across South Carolina that pose human health concerns. Ecological drivers, particularly those in deciduous and mixed forests and agriculture/livestock land, were found to play crucial roles in Rickettsia spp. distribution statewide. Despite limited clinical evidence of R. amblyommatis pathogenicity, this SFGR species has the potential to play a regulatory role through a possible inhibitory interaction with R. parkeri. Ticks carrying R. amblyommatis are suggested to be less likely to be infected with and possibly transmit fewer SFGR pathogens; thus this species has major public health interest. Given the higher rates of A. americanum across the state and the low presence of R. parkeri cases described, we hypothesized that higher prevalence of R. amblyommatis-infected A. americanum could translate to lower SFGR transmission among humans. Finally, the differences in environmental factors and Rickettsia species positivity suggest that enzootic cycles are influenced by microclimate conditions. These suggestions should be further explored using a One Health approach to better understand the possible relationships between SFGR species, tick populations, animal hosts, and the environmental factors affecting their distribution.

Availability of data and materials

The data supporting the findings of the study must be available within the article and/or its supplementary materials, or deposited in a publicly available database.

Abbreviations

aOR:

Adjusted odds ratio

CI:

Confidence interval

CO2:

Carbon dioxide

ETCI:

Equal tail credible interval

IFA:

Indirect fluorescent antibody assay

MCMC:

Markov chain Monte Carlo

NC:

North Carolina

OR:

Odds ratio

PCR:

Polymerase chain reaction

qPCR:

Quantitative polymerase chain reaction

RMSF:

Rocky Mountain spotted fever

RT-PCR:

Real time polymerase chain reaction

SC:

South Carolina

SFGR:

Spotted fever group Rickettsia sp.

References

  1. Blanton LS. The rickettsioses: a practical update. Infect Dis Clin North Am. 2019;33:213–29.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Paddock CD, Sumner JW, Comer JA, Zaki SR, Goldsmith CS, Goddard J, et al. Rickettsia parkeri: a newly recognized cause of spotted fever rickettsiosis in the United States. Clin Infect Dis. 2004;38:805–11.

    Article  PubMed  Google Scholar 

  3. Ammerman NC, Swanson KI, Anderson JM, Schwartz TR, Seaberg EC, Glass GE, et al. Spotted-fever group Rickettsia in Dermacentor variabilis. Maryland Emerg Infect Dis. 2004;10:1478–81.

    Article  PubMed  Google Scholar 

  4. Ribeiro CM, Carvalho JLB, Bastos PAS, Katagiri S, Batalha EY, Okano W, et al. Prevalence of Rickettsia rickettsii in ticks: systematic review and meta-analysis. Vector-Borne Zoonot Dis. 2021;21:557–65.

    Google Scholar 

  5. Bishop A, Borski J, Wang H-H, Donaldson TG, Michalk A, Heldman S, et al. Increasing incidence of spotted fever group rickettsioses in the United States, 2010–2018. Vector-Borne Zoonot Dis. 2022;22:491–7.

    Article  Google Scholar 

  6. Macaluso KRP, Christopher D. Tick-borne Spotted Fever Group Rickttsioses and Rickettsia species. In: Sonenshine DER, R. Michael, editor. Biology of Ticks. Second ed: Oxford University Press; 2014.

  7. Yen W-Y, Stern K, Mishra S, Helminiak L, Sanchez-Vicente S, Kim HK. Virulence potential of Rickettsia amblyommatis for spotted fever pathogenesis in mice. Pathogens Dis. 2021;79:ftab024.

    Article  CAS  Google Scholar 

  8. Snellgrove AN, Krapiunaya I, Scott P, Levin ML. Assessment of the Pathogenicity of Rickettsia amblyommatis, Rickettsia bellii, and Rickettsia montanensis in a Guinea Pig Model. Vector-Borne Zoonot Dis. 2021;21:232–41.

    Article  Google Scholar 

  9. Apperson CS, Engber B, Nicholson WL, Mead DG, Engel J, Yabsley MJ, et al. Tick-borne diseases in North Carolina: Is “Rickettsia amblyommii” a possible cause of rickettsiosis reported as Rocky Mountain spotted fever? Vector-Borne Zoonot Dis. 2008;8:597–606.

    Article  Google Scholar 

  10. Richardson EA, Roe RM, Apperson CS, Ponnusamy L. Rickettsia amblyommatis in ticks: a review of distribution, pathogenicity, and diversity. Microorganisms. 2023;11:493.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Jiang J, Yarina T, Miller MK, Stromdahl EY, Richards AL. Molecular detection of Rickettsia amblyommii in Amblyomma americanum parasitizing humans. Vector Borne Zoonot Dis. 2010;10:329–40.

    Article  Google Scholar 

  12. Hecht JA, Allerdice MEJ, Dykstra EA, Mastel L, Eisen RJ, Johnson TL, et al. Multistate survey of American dog ticks (Dermacentor variabilis) for Rickettsia Species. Vector-Borne Zoonot Dis. 2019;19:652–7.

    Article  Google Scholar 

  13. Burgdorfer W, Hayes SF, Thomas LA, Lancaster J. New spotted fever group Rickettsia from the lone star tick, Amblyomma americanum. Rickettsiae and rickettsial diseases/edited by W Burgdorfer; RL Anacker. 1981.

  14. Weller SJ, Baldridge GD, Munderloh UG, Noda H, Simser J, Kurtti TJ. Phylogenetic placement of rickettsiae from the ticks Amblyomma americanum and Ixodes scapularis. J Clin Microbiol. 1998;36:1305–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Center for Disease Control and Prevention. Rocky Mountain Spotted Fever (RMSF). Epidemiology and statistics. [website] cdc.gov2022. 2022. https://www.cdc.gov/rmsf/stats/index.html.

  16. Smith MP, Ponnusamy L, Jiang J, Ayyash LA, Richards AL, Apperson CS. Bacterial pathogens in ixodid ticks from a Piedmont County in North Carolina: prevalence of rickettsial organisms. Vector-Borne Zoonot Dis. 2010;10:939–52.

    Article  Google Scholar 

  17. Lee S, Kakumanu ML, Ponnusamy L, Vaughn M, Funkhouser S, Thornton H, et al. Prevalence of Rickettsiales in ticks removed from the skin of outdoor workers in North Carolina. Parasit Vectors. 2014;7:607.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Barrett A, Little SE, Shaw E. “Rickettsia amblyommii” and R. montanensis infection in dogs following natural exposure to ticks. Vector-Borne Zoonot Dis. 2014;14:20–5.

    Article  Google Scholar 

  19. Williamson PC, Billingsley PM, Teltow GJ, Seals JP, Turnbough MA, Atkinson SF. Borrelia, Ehrlichia, and Rickettsia spp. in ticks removed from persons, Texas, USA. Emerg Infect Dis. 2010;16:441.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Billeter SA, Blanton HL, Little SE, Levy MG, Breitschwerdt EB. Detection of “Rickettsia amblyommii” in association with a tick bite rash. Vector-Borne Zoonot Dis. 2007;7:607–10.

    Article  Google Scholar 

  21. Wright CL, Sonenshine DE, Gaff HD, Hynes WL. Rickettsia parkeri transmission to Amblyomma americanum by cofeeding with Amblyomma maculatum (Acari: ixodidae) and potential for spillover. J Med Entomol. 2015;52:1090–5.

    Article  CAS  PubMed  Google Scholar 

  22. Levin ML, Schumacher LBM, Snellgrove AN. Effects of Rickettsia amblyommatis infection on the vector competence of Amblyomma americanum ticks for Rickettsia rickettsii. Vector-Borne Zoonot Dis. 2018;18:579–87.

    Article  Google Scholar 

  23. Winberry JJ, Bushman DO. South Carolina. Encyclopedia Britannica (updated 8 November 202214 December 2022). 2022. https://www.britannica.com/place/South-Carolina.

  24. Nolan MS, Dye-Braumuller K, Evans C, Rustin L, Lachenmyer E, Neault M. The emergence of ticks and tick-borne disease in South Carolina: establishing a statewide Initiative leveraging academic-public partnerships. American Society of Tropical Medicine and Hygiene; Seattle, WA: American Journal of Tropical Medicine and Hygiene; 2022. p. 237.

  25. South Carolina 2020. https://data.census.gov/cedsci/profile?g=0400000US45. Accessed 16 Jan 2022.

  26. South carolina department of natural resources. South Carolina state climatology office dnr.sc.gov2020 https://www.dnr.sc.gov/climate/sco/ClimateData/cli_sc_climate.php. Accessed 2023.

  27. Center for disease control and prevention. Guide to the surveillance of metastriate ticks (Acari: ixodidae) and their pathogens in the United States. 2020.

  28. Keirans JE, Litwak TR. Pictorial key to the adults of hard ticks, family ixodidae (ixodida: ixodoidea), east of the Mississippi River. J Med Entomol. 1989;26:435–48.

    Article  CAS  PubMed  Google Scholar 

  29. Egizi AM, Robbins RG, Beati L, Nava S, Vans CR, Occi JL, et al. A pictorial key to differentiate the recently detected exotic Haemaphysalislongicornis Neumann, 1901 (Acari, Ixodidae) from native congeners in North America. Zookeys. 2019;818:128.

    Article  Google Scholar 

  30. Clifford CM, Anastos G. The use of chaetotaxy in the identification of larval ticks (Acarina: ixodidae). J Parasitol. 1960;46:567–78.

    Article  CAS  PubMed  Google Scholar 

  31. Williams CJ, Moffitt CM. A critique of methods of sampling and reporting pathogens in populations of fish. J Aquat Anim Health. 2001;13:300–9.

    Article  Google Scholar 

  32. Environmental Protection Agency. Ecoregion Download Files by State - Region 4. (updated September 23, 202205/14/2023). 2022. https://www.epa.gov/eco-research/level-iii-and-iv-ecoregions-state.

  33. McMahan CS, Tebbs JM, Hanson TE, Bilder CR. Bayesian regression for group testing data. Biometrics. 2017;73:1443–52.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Dewitz J. NLCD 2019 Land Cover Conterminous United States. In: (USGS) USGS, editor.2021.

  35. Diuk-Wasser MA, VanAcker MC, Fernandez MP. Impact of land use changes and habitat fragmentation on the eco-epidemiology of tick-borne diseases. J Med Entomol. 2020;58:1546–64.

    Article  Google Scholar 

  36. Nicholson WR, Sonenshine DE, Noden BH, Brown RN. Ticks (Ixodida). In: Mullen GR, Durden LA, editors. Medical and veterinary entomology. Third edition ed: Academic press; 2019.

  37. Tsao JI, Hamer SA, Han S, Sidge JL, Hickling GJ. The contribution of wildlife hosts to the rise of ticks and tick-borne diseases in North America. J Med Entomol. 2021;58:1565–87.

    Article  PubMed  Google Scholar 

  38. Satjanadumrong J, Robinson MT, Hughes T, Blacksell SD. Distribution and ecological drivers of spotted fever group Rickettsia in Asia. EcoHealth. 2019;16:611–26.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Quintero VJC, LE Paternina T, Uribe YA, Muskus C, Hidalgo M, Gil J, et al. Eco-epidemiological analysis of rickettsial seropositivity in rural areas of Colombia: a multilevel approach. PLoS Negl Trop Dis. 2017;11:e0005892.

    Article  Google Scholar 

  40. Carlson CJ, Burgio KR, Dougherty ER, Phillips AJ, Bueno VM, Clements CF, et al. Parasite biodiversity faces extinction and redistribution in a changing climate. Sci Adv. 2017;3:e1602422.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Lane RS, Fedorova N, Kleinjan JE, Maxwell M. Eco-epidemiological factors contributing to the low risk of human exposure to ixodid tick-borne borreliae in southern California, USA. Ticks Tick Borne Dis. 2013;4:377–85.

    Article  PubMed  Google Scholar 

  42. Allan BF, Goessling LS, Storch GA, Thach RE. Blood meal analysis to identify reservoir hosts for Amblyomma americanum ticks. Emerg Infect Dis. 2010;16:433–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Harmon JR, Scott MC, Baker EM, Jones CJ, Hickling GJ. Molecular identification of Ehrlichia species and host bloodmeal source in Amblyomma americanum L. from two locations in Tennessee United States. Ticks Tick Borne Dis. 2015;6:246–52.

    Article  PubMed  Google Scholar 

  44. Center for Disease Control and Prevention. Regions Where Ticks Live. cdc.gov2022. 2023. https://www.cdc.gov/ticks/geographic_distribution.html.

  45. Childs JE, Paddock CD. The ascendancy of Amblyomma americanum as a vector of pathogens affecting humans in the United States. Annu Rev Entomol. 2003;48:307–37.

    Article  CAS  PubMed  Google Scholar 

  46. Ayres BN, James AM, Wehtje ME, Nicholson WL. Rickettsial agents detected in ixodid ticks (Acari: ixodidae) collected from Sus scrofa (Artiodactyla: suidae) in Florida and South Carolina. J Entomol Sci. 2022;57:351–62.

    Google Scholar 

  47. Mixson TR, Campbell SR, Gill JS, Ginsberg HS, Reichard MV, Schulze TL, et al. Prevalence of Ehrlichia, Borrelia, and Rickettsial agents in Amblyomma americanum (Acari: ixodidae) collected from Nine States. J Med Entomol. 2006;43:1261–8.

    PubMed  Google Scholar 

  48. De Jesus CE, Ganser C, Kessler WH, White ZS, Bhosale CR, Glass GE, et al. A survey of tick-borne bacterial pathogens in Florida. Insects. 2019;10:297.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Stanley H, Rhodes DVL. Presence of Rickettsia species in ticks collected from companion animals in Northeastern Georgia United States. Veterin Sci. 2021;8:37.

    Google Scholar 

  50. Delisle J, Mendell NL, Stull-Lane A, Bloch KC, Bouyer DH, Moncayo AC. Human infections by multiple spotted fever group rickettsiae in Tennessee. Am Soc Tropic Med Hygiene. 2016;94:1212–7.

    Article  CAS  Google Scholar 

  51. Vaughn MF, Delisle J, Johnson J, Daves G, Williams C, Reber J, et al. Seroepidemiologic study of human infections with spotted fever group rickettsiae in North Carolina. J Clin Microbiol. 2014;52:3960–6.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Kensinger BJ, Allan BF. Efficacy of dry ice-baited traps for sampling Amblyomma americanum (Acari: ixodidae) varies with life stage but not habitat. J Med Entomol. 2011;48:708–11.

    Article  CAS  PubMed  Google Scholar 

  53. Raghavan RK, Goodin DG, Hazlicek GA, Zolnerowich G, Dryden MS, Anderson GA, et al. Maximum entropy-based ecological niche model and bio-climatic determinants of lone star tick (Amblyomma americanum) Niche. Vector-Borne Zoonot Dis. 2016;16:205–11.

    Article  Google Scholar 

  54. Benelli G. Pathogens manipulating tick behavior—through a glass. Darkly Pathogens. 2020;9:664.

    Article  CAS  PubMed  Google Scholar 

  55. Macaluso KR, Sonenshine DE, Ceraul SM, Azad AF. Rickettsial infection in Dermacentor variabilis (Acari: ixodidae) inhibits transovarial transmission of a second Rickettsia. J Med Entomol. 2002;39:809–13.

    Article  PubMed  Google Scholar 

  56. Budachetri K, Browning RE, Adamson SW, Dowd SE, Chao C-C, Ching W-M, et al. An insight into the microbiome of the Amblyomma maculatum (Acari: ixodidae). J Med Entomol. 2014;51:119–29.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgement

Not applicable.

Funding

This work was partially funded through the US Centers for Disease Control and Prevention’s Epidemiology and Laboratory Capacity for Prevention and Control of Emerging Infectious Diseases Cooperative Agreement #5 NU50CK000542-04–00. The conclusions, findings, and opinions expressed by authors contributing to this article do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions.

Author information

Authors and Affiliations

  1. Department of Epidemiology and Biostatistics, University of South Carolina, Columbia, SC, USA

    Lídia Gual-Gonzalez, Stella C. W. Self, Kia Zellars, Madeleine Meyer, Kyndall C. Dye-Braumuller, Myriam W. Torres & Melissa S. Nolan

  2. Vector-Borne Diseases Laboratory, South Carolina Department of Health and Environmental Control, South Carolina Public Health Laboratory, Columbia, SC, USA

    Chris L. Evans

  3. Grupo de Biología y Control de Enfermedades Infecciosas, Universidad de Antioquia, Medellín, Colombia

    Omar Cantillo-Barraza

Authors
  1. Lídia Gual-Gonzalez
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Stella C. W. Self
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Kia Zellars
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Madeleine Meyer
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Kyndall C. Dye-Braumuller
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Chris L. Evans
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Omar Cantillo-Barraza
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Myriam W. Torres
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Melissa S. Nolan
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LG-G, SCWS, CLE, OC-B, MET, and MSN participated in the study design process; KCD-B and CLE collected the data. LG-G, KCD-B, KZ, and MM developed, validated, and performed laboratory analysis. LG-G and SCWS developed and performed data analysis. LG-G, MSN, and SCWS interpreted the data. All authors reviewed the manuscript.

Corresponding author

Correspondence to Melissa S. Nolan.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for Publication

Not applicable.

Competing interests

Stella Self receives research, travel, and consulting funds from the Companion Animal Parasite Council.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Primers and probes used for R. amblyommatis and R. parkeri amplification. Table S2. Land-use classification name descriptions.

Additional file 2:

Markov chain Monte Carlo analysis technical results.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gual-Gonzalez, L., Self, S.C.W., Zellars, K. et al. Eco-epidemiology of Rickettsia amblyommatis and Rickettsia parkeri in naturally infected ticks (Acari: Ixodida) from South Carolina. Parasites Vectors 17, 33 (2024). https://doi.org/10.1186/s13071-023-06099-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-023-06099-z

Keywords

Parasites & Vectors

ISSN: 1756-3305

Contact us