https://orcid.org/0000-0001-6795-3541
Figures
Abstract
In this study, we estimated exposure for Scrub typhus (STG), Typhus (TG) and Spotted fever groups (SFG) rickettsia using serology at a fine scale (a whole sub-district administration level) of local communities in Nan Province, Thailand. Geographical characteristics of the sub-district were divided into two landscape types: lowland agricultural area in an urbanized setting (lowland-urbanized area) and upland agricultural area located close to a protected area of National Park (upland-forested area). This provided an ideal contrast between the two landscapes with low and high levels of human-altered habitats to study in differences in disease ecology. In total, 824 serum samples of participants residing in the eight villages were tested by screening IgG ELISA, and subsequently confirmed by the gold standard IgG Immunofluorescent Assay (IFA). STG and TG IgG positivity were highest with seroprevalence of 9.8% and 9.0%, respectively; whereas SFG positivity was lower at 6.9%. Inhabitants from the villages located in upland-forested area demonstrated significantly higher STG exposure, compared to those villages in the lowland-urbanized area (chi-square = 51.97, p < 0.0001). In contrast, TG exposure was significantly higher in those villagers living in lowland-urbanized area (chi-square = 28.26, p < 0.0001). In addition to the effect of landscape types, generalized linear model (GLM) analysis identified socio-demographic parameters, i.e., gender, occupation, age, educational level, domestic animal ownership (dog, cattle and poultry) as influential factors to explain the level of rickettsial exposure (antibody titers) in the communities. Our findings raise the public health awareness of rickettsiosis as a cause of undiagnosed febrile illness in the communities.
Author summary
Evidence of human exposures to rickettsial pathogens were reported from a cross-sectional study at a whole sub-district scale of local communities in Nan Province, Thailand. Seroprevalence and level of rickettsial exposures demonstrated differences between the habitat types, ecological aspects and socio-demographic factors. In addition, abundance of domestic animals in the community appeared to be one of significant factors influencing levels of human exposure to rickettsial pathogens. Our findings will benefit the local public health by raising awareness of rickettsial infections as one of potential health concerns in the community. Inclusion of rickettsioses in routine laboratory diagnosis would help to differentiate unknown febrile illness and guide appropriate treatment. Further studies are required, particularly in the fields of disease ecology as well as medical and veterinary entomology, in order to better understand epidemiology and potential zoonotic transmission of these neglected rickettsioses in endemic areas.
Citation: Chaisiri K, Tanganuchitcharnchai A, Kritiyakan A, Thinphovong C, Tanita M, Morand S, et al. (2022) Risk factors analysis for neglected human rickettsioses in rural communities in Nan province, Thailand: A community-based observational study along a landscape gradient. PLoS Negl Trop Dis 16(3): e0010256. https://doi.org/10.1371/journal.pntd.0010256
Editor: Wen-Ping Guo, Chengde Medical University, CHINA
Received: September 3, 2021; Accepted: February 12, 2022; Published: March 23, 2022
Copyright: © 2022 Chaisiri et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. For the purpose of Open Access, the author has applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: SM received the French ANR Project FutureHealthSEA (grant number: ANR-17-CE35-0003-02) “Predictive scenarios of health in Southeast Asia: linking land use and climate changes to infectious diseases”. SDB is funded by the Wellcome Trust of the United Kingdom [220211]. This research was funded in whole, or in part, by the Wellcome Trust [220211]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Human rickettsioses are neglected emerging and re-emerging diseases caused by a group of obligate intracellular bacteria in the Order Rickettsiales (e.g., Rickettsia and Orientia). Rickettsial pathogens causing disease in humans can be divided into the three main groups, known as scrub typhus (STG), typhus (TG) and spotted fever groups (SFG) [1]. These zoonotic bacteria are transmitted to humans and animals after bites by arthropod vectors, including ticks, mites, fleas or lice. Co-occurrences of rickettsial species among arthropod vectors, domestic animals and humans as well as increasing opportunities for interactions at the human-animal interface are of importance in zoonotic transmission cycles [2–4]. These vector-borne diseases have increasingly been recognized as one of the common causes of febrile illness in Southeast Asia, in addition to dengue, leptospirosis and malaria [1,5,6].
Scrub typhus is a chigger-borne rickettsiosis caused by intracellular bacteria in the genus Orientia, including Orientia tsutsugamushi (reported primarily in Asia-Pacific region), Orientia chuto (in the Arabian Peninsula and African countries) and Orientia sp. (Candidatus Orientia chiloensis, a potential new species from Chile) [7–11]. Larval stage of mites in the family Trombiculidae, known as “chiggers”, are the main vector of the disease. Previous studies have suggested that human scrub typhus is more common in rural or sub-urban areas [12,13] with the majority of disease in farmers, military soldiers and jungle trekkers [7,14,15]. Typhus group (TG) mainly involves murine typhus (also known as endemic typhus) infections caused by Rickettsia typhi, a flea-borne rickettsiosis with worldwide distribution, notably in tropical coastal regions [16,17]; and Rickettsia prowazekii which causes epidemic typhus, a louse-borne rickettsiosis in the areas of poor socioeconomic status and high prevalence of human body louse infestation [18]. Transmission of typhus group rickettsioses generally occur when the bacteria in the feces of arthropod vectors are rubbed into the bite wound or the vectors themselves are crushed and broke in to the skin [19]. For the spotted fever group (SFG), there are more than 25 rickettsial species and related strains distributed globally, and known to cause disease in human [20], including Rickettsia africae, Rickettsia asiatica, Rickettsia australis, Rickettsia conorii, Rickettsia heilongjiangensis, Rickettsia helvetica, Rickettsia honei, Rickettsia japonica, Rickettsia parkeri, Rickettsia raoultii, Rickettsia rickettsia, Rickettsia sibirica and Rickettsia tamurae. These SFG rickettsia are mainly harbored by ticks, and transmitted to human or other animal hosts via their bites. Other ectoparasites such as fleas and mites are also involved as the vectors of several SFG rickettsia including the flea-borne Rickettsia felis and the mite-borne Rickettsia akari [2,21,22].
Diagnosis of these pathogens relies on molecular techniques such as PCR-based amplification to demonstrate the presence of the infectious agent, as well as serological methods as a surrogate including enzyme-linked immunosorbent assay (ELISA) and immunofluorescence assays (IFA) [23]. Serological investigations allow epidemiologists to investigate the history of rickettsial exposure and estimation of seroprevalence by detecting the occurrence of preexisting levels of antibodies in a population [24–26]. IFA detection of IgM and/or IgG is recognized as a gold standard for rickettsial serological diagnosis, particularly for STG, however there remains uncertainty as to the most suitable diagnostic cut-off for each geographic location [23,27]. ELISA-based diagnostics give acceptable levels of sensitivity and specificity for STG, SFG and TG antibody detection [23,28,29] as well as having increased diagnostic throughput and less subjectivity [29,30]. However, serological results are normally confined to the serogroup level. Although generally, there are cross-reactions within the serogroups, they occur less so between serogroups with the exception that cross-reactions between SFG and TG antibodies are recognised to a limited extent [1,31].
In Thailand, previous studies have indicated that all three rickettsial groups exist and are a cause for public health concern [2,5,32–34]. In terms of human febrile illness, scrub typhus is the most prevalent rickettsial pathogen which has markedly increased during the last two decades although this may be due to increased awareness and access to improved diagnostics. During 2003–2018, there were approximately 100,000 cases reported to the National Disease Surveillance system (Ministry of Public Health, Thailand). The disease is reported nation-wide, with the majority of the cases from the northern region accounting approximately half (53%) of the total reported cases per year [34]. Murine typhus is also recognized a rickettsiosis in Thailand. Case reports of R. typhi infection have been documented in several provinces, either in municipalities, rural communities or border areas [5,35–38]. Epidemiological data regarding SFG rickettiosis is rather more scarce compared to scrub typhus and murine typhus. In 1994, three human SFG rickettsiosis cases from Chiang Mai Province were confirmed by serology, and documented as the first report in Thailand [39]. The first human case of R. felis infection in Thailand (and in Asia) was reported at the Thai-Myanmar border in Kanchanaburi Province (32). Thai tick typhus Rickettsia TT-118 (also known as R. honei) was identified from a patient in urban area of Bangkok who had history of camping at a National Park prior to the onset of fever [40]. In addition, Rickettsia spp. closely related to R. japonica and R. helvetica have also reported to infect humans in Thailand [41–43].
Most of the case reports and epidemiological studies in Thailand were recognised in the context of passive surveillance or retrospective studies from stored clinical samples [44], while active surveillance data of human rickettsiosis is still limited. Here, within the research framework of FutureHealthSEA project (Predictive scenarios of health in Southeast Asia: linking land use and climate changes to infectious diseases, ANR-17-CE35-0003), we performed screening for human rickettsiosis at the sub-district administration level of local communities in Nan Province, northern Thailand. The objectives of the study were: (1) to estimate level of past exposure to rickettsial pathogens using serological diagnosis approach; and (2) to determine potential factors, such as socio-demographic, zoonotic potential and environmental parameters influencing rickettsial exposures in a local scale setting.
Materials and methods
Ethics statement
All procedures involving human samples and data collection were reviewed and approved by the Ethics Committee of the Faculty of Tropical Medicine, Mahidol University (document no. MUTM 2018-035-01). Signed formal consent form was obtained from all the participants and also from the parent/guardian of the child participants.
Study population and individual demographic information
In 2018, we conducted a cross-sectional study at the local communities, all of the 8 villages in Saen Thong Sub-district, Tha Wang Pha District, Nan Province (latitude 19.133 and longitude 100.7680). Geographical distribution of the 8 villages is shown in Fig 1. In terms of geographical characteristics, the sub-district can be roughly divided into the two landscape types: (1) lowland agricultural area in an urbanized setting (Lowland-urbanized area), closed to Tha Wang Pha city on the east; and (2) upland agricultural area closed to the forested area of Nantaburi National Park (Upland-forested area) on the west.
Distribution and level of (A) STG, (B) TG and (C) SFG individual exposures from each village. Green circles indicate negativity. Red circles show positive exposures, and size of the circles indicate level of IgG IFA positivity at a diagnostic titer cut-off of 1:100.
The target population of this study was randomly selected with the following inclusion criteria: both males and females; age over 9 years old; recently and permanently residing in the study area; able to adequately respond to the questions and voluntary participation in the research project. Individual data regarding participants’ personal information was collected including age group (children at 9–17 years old; young adult at 18–35 years old; adult at 36–60 years old; and elderly at >60 years old), gender (male and female) and occupation (farmer or non-farmer). Information on companion pets and livestock ownership, which included numbers of dogs, cats, poultry, pigs and cattle, were also obtained to determine potential zoonotic association. Every household that the participants’ resided were geo-tagged to collect GPS coordinates information for further spatial case mapping.
Blood sample collection and serological diagnosis of rickettsial past exposure
Staff from the local primary health care unit (PCU) made an appointment with the village chiefs for preparation of participant recruitment for blood collection. Blood samples were collected via venipuncture (5ml) in non-anticoagulant tubes (VACUETTE TUBE Serum Clot Activator, Greiner Bio-One, Austria), and then centrifuged at 2,795 xg for 5 min within the same day of collection to separate serum samples in sterile 1.5 ml collection tubes (Axygen, US). The serum samples were frozen in -80°C prior to laboratory investigations.
Due to the large number of samples, sera were screened for the presence of IgG antibodies at 1:100 dilution using STG, TG or SFG IgG ELISA using diagnostic cut off with a high sensitivity of 0.50 nett OD. ELISA-positive samples were titrated in the “gold standard” IFA to determine the endpoint antibody titer. The MORU in-house STG IgG ELISA used specific antigens of O. tsutsugamushi Karp, Kato, Gilliam, and TA716 strains to detect scrub typhus IgG antibodies [45]. For the TG ELISA, R. typhi Wilmington strain and for SFG ELISA, R. conorii and R. honei antigens were applied using the same methodology. Consequently, we performed IgG IFA to produce quantitative results, with the full methodology described elsewhere [46]. Briefly, sera were serially 2-fold diluted from 1:100 to 1:25,600 and the endpoint was determined as the highest titer displaying specific fluorescence. Titer >1:100 was considered positive for the presence of IgG antibodies.
Statistical analysis
Estimation of STG, TG and SFG by village, landscape type and personal attributes were performed using descriptive statistics, i.e., seroprevalence with 95% confidence interval. A Venn diagram was created using “ggplot2” package [47] in R freeware [48] to illustrate number and proportion of participants exposed with STG, TG, SFG and co-exposures. Seropositive cases were geographically mapped to reveal distribution patterns of STG, TG and SFG using “tmap” [49] and “raster” [50] packages in R freeware [48].
Generalized linear model (GLM) was performed using R freeware [48] to determine socio-demographic and environmental factors explaining the level of STG, TG and SFG exposures (titers). Similar initial models were applied to the three rickettsial groups, comprising with the following variables: Landscape type, Gender, Occupation, Age, Educational level, as well as discrete variables for population numbers of dogs, cats, pigs, poultry and cattle. Poisson distribution was selected for model fitting after checking data distribution using “fitdistrplus” package [51]. GLM with Akaike’s information criterion corrected for sample size (AICc) was performed with negative binomial distribution using “glmulti” package [52]. The best models for each rickettsial group were selected for further discussion based on low AICc and high Akaike’s weight (Wr) scores [53]. To identify the degree of multicollinearity among explicative variables, the variance inflation factor (VIF) was computed using “car” package [54]. In addition, quality of the selected models was assessed with goodness-of-fit and power analysis tests using “sjstats” [55] and “lmSupport” [56] packages, respectively.
Results and discussion
Socio-demographic information of Saen Thong subdistrict
In total, 824 sera from the eight villages were subjected to serological testing yielding a participation rate of 19.8% (total registered population of 4,145 in the subdistrict). Of the total population registered in the system, it is likely that a proportion of the population (particularly those working age) did not live in the area because they moved to work in other cities during the study which may have artificially lowered the participation rate. Serum samples were collected from participants who lived in 524 households covering almost a half (48.3%) of the total number of eligible households (1,083 registered households in total). Adults were the dominant age group (477 samples, 57.8%), followed by the elderly (250 samples, 30.3%), children (45 samples, 5.4%), young adults (37 samples, 4.5%), however 15 samples (1.8%) were not able to assigned an age group category as this information was not recorded during the interview. The majority of the participants were farmers (601 samples, 72.9%), with non-farmers accounting for 21.9% (181 samples) and unknown occupation for 42 samples (5.1%). We were unable to assign their occupation as this information was not recorded during the interview (more details in S1 Table).
Over a half of the participants (483 individuals, 58.6%) kept or raised animals including companion pets and farm livestock at home, whereas 341 individuals (41.3%) responded as having no animals. Among these, poultry (mainly backyard chickens) were the most popular (375 people, 45.5%) followed by 185 people (22.4%) responded to have dogs, 57 people (6.9%) cats, 14 people (1.7%) cattle and 6 people (0.7%) raised pigs at home.
Rickettsiosis past exposure in Saen Thong subdistrict: Associations of environmental and socio-demographic factors
Rickettsial exposure was noted in the Saen Thong subdistrict with seroprevalence of STG, TG and SFG of 9.8%, 9.0% and 6.9%, respectively (Table 1). Comparing the results of STG seroprevalence presented here to previous publications, the study population in Nan Province demonstrated a similar rate of exposure. This agreed with the review data of Bonell et al. (2017) that Orientia tsutsugamushi infection endemic in several Asian countries (i.e., Bangladesh, Indonesia, Lao PDR, Malaysia, Papua New Guinea and Sri Lanka), with seroprevalence ranging from 9.3%-27.9% [44]. Compared with previous studies in Thailand, the STG seroprevalence in Nan province was slightly higher than a previous reported among the populations in the three previously studied provinces of Thailand (4.2%), i.e., Chiang Rai, Khon Kaen and Nakhon Phanom provinces [33], however it should be noted that the study used a different antibody cut-off titer (≥ 1:128 compared with ≥1:100 in this study to detect IgG using IFA) to define positivity. In addition, our result demonstrated much lower STG positivity rate than the population in the western provinces, i.e., 62.5% Ratchaburi, 64.6% Petchaburi and 49.1% Kanchanaburi province [57], however it is important to recognise that a lower diagnostic cut-off value was used for interpretation (IgG using IFA, ≥1:50 compared with ≥1:100 in our study). In terms of murine typhus, seroprevalence of TG in the present study was slightly higher (9.0%) than the study in Chiang Rai, Khon Kaen and Nakhon Phanom provinces (4.2%), [33]; and similar rate of exposure with a previous study (8.0%) in suburban Bangkok which used the indirect immunoperoxidase assay with a diagnostic cut-off of ≥1:50 [58]. For SFG, our result in Nan province demonstrated higher seroprevalence (6.9%) than that from a previous study in Chiangrai, Khon Kaen and Nakhon Phanom provinces (0.8%), [33]. In contrast, another seropositivity study reported much higher SFG exposures in Chiang Rai (33%) and Tak provinces (27.3%), but the study used different assays (SFGR Enzyme Immunoassay IgG Antibody Kit), [59]. Therefore, it is important to consider the serological diagnostic methods and selection of the diagnostic cut-off values when comparing exposure rates, otherwise this has the potential to cause analysis bias and misinterpretation.
There was evidence of multiple positivity with the three rickettsial serogroups in a small proportion of participants. However, no individual demonstrated multiple positivity with all three rickettsial groups. Multiple positivity was occurred in 10 (5.2%), 6 (3.1%) and 3 (1.6%) cases between STG-SFG, TG-SFG and STG-TG, respectively (Fig 2). However, it is recognised that serological cross-reactivity among SFG and TG [60] may occur, which could account for the dual-positivity between these antigenic groups.
Relative proportions of the rickettsiosis exposures are provided in parentheses.
GLM analysis revealed that environmental factors, including landscape type and socio-demographic parameters, such as gender, occupation, age, educational level, and number of domestic animals in possession (dog, cattle and poultry) was significantly associated with higher rickettsial antibody levels in the communities (Table 2). Inhabitants from the villages located in upland-forested area (Village 4, 5, 6 and 7) demonstrated significantly higher STG seroprevalence (9.2% to 35.2%) and higher level of STG exposure (titer), compared to those villages (Village 1, 2, 3 and 8) in the lowland-urbanized area (none to 4.1%), (chi-square = 51.97, p < 0.0001), (Figs 1A and 3). Individuals located in upland-forested areas (increased habitat complexity and biodiversity) would have higher rates of exposures to chigger mites resulting in higher STG seroprevalence than those who live in urban environments (7,34). The reverse trend was observed for TG exposure, where the villages in the lowland-urbanized area demonstrated significantly higher seroprevalence (7.1% to 24.2%) and higher TG titer level compared to those villages in the upland-forested area (1.6% to 4.1%), (chi-square = 28.26, p < 0.0001), (Figs 1B and 4). Individuals located in urban environments may have potentially higher rates of exposure to fleas (the vector for murine typhus) parasitizing on rats and companion pets [61,62] resulting in higher TG seroprevalence compared with those individuals that are located in more forested or non-urban environments. These trends are in accordance with the study of scrub typhus and murine typhus cases in Lao PDR [12] and Malaysia [63]. We found no significant effect of landscape type neither on SFG seroprevalence nor antibody titer level of SFG exposure. Distribution maps of STG-, TG- and SFG-positive individuals from each village are provided in S1–S3 Figs.
Similar initial models are applied comprising with the following variables: Landscape type, Gender, Occupation, Age, Educational level, Dog number, Cat number, Pig number, Poultry number and Cattle number (with Poisson distribution after checking data distribution using “fitdistrplus” package in R freeware). Analysis of deviance (type II test) significant level for each variable is indicated as following: *p = 0.05–0.01; **p = 0.01–0.001 and ***p < 0.001. VIF = Variance inflation factor; Wr = Akaike’s weights.
Male participants demonstrated significantly higher STG seroprevalence (15.9%) when compared with females (5.3%), (chi-square = 19.47, p < 0.0001), (see Table 1). This may be explained by the need for males to enter into the forest for hunting or animal tending purposes. For TG and SFG positivity, no clear difference was observed between the two genders. Surprisingly, there was no discernible difference in rickettsial seroprevalence between farmer and non-farmer participants. One would normally associate farming with higher levels of scrub typhus seroprevalence, however it would depend where the farming was performed for example, whether the farming was performed close or remote from forested areas. In terms of educational level, the rickettsial exposures were found to be significantly more common among those with lower levels of education compared to those with higher education which had been previously identified in Lao PDR [12].
Domestic animals as influencing factor of rickettsiosis exposure in Saen Thong subdistrict
In terms of the human-animal interface, our results demonstrate a significant link between number of domestic animals owned and the level of human rickettsial exposures. Higher levels of STG, SFG and TG exposures were positively correlated with number of dogs, poultry and both dogs and cattle, respectively (Figs 3–5). Domestic animals can act as natural vertebrate reservoir hosts for Rickettsia spp. in Southeast Asia [2], together with the close contact between people, domestic animals and associated ectoparasites (e.g., ticks, mites, lice and fleas). These bring emerging or remerging zoonoses closer to humans.
Marks on the X-axis represent individual observations.
Marks on the X-axis represent individual observations.
Marks on the X-axis represent individual observations.
Dogs are one of the prime examples of this phenomenon as they have long been associated with humans for approximately 23,000 years ago, being the first domesticated species on earth [64]. Using a data science approach and network analysis, Morand et al. (2020), [3] revealed that domestic dogs (Canis lupus familiaris) were the most important species sharing a great number of potential zoonotic rickettsial species with human (e.g., R. conorii, R. felis, R. rickettsii, R. typhi, O. tsutsugamushi, Ehrlichia canis and Neorickettsia sennetsu) when compared with other animals. Recent publications have confirmed the role of domestic dogs as important reservoirs for Rickettsia felis [4,65]. The study, Ng-Nguyen et al. (2020), [4] demonstrated that dogs have ability to sustain R. felis for a long period being asymptomatic during the course of infection; and they can also act as biological vehicles passing the rickettsial bacteria to uninfected fleas (Ctenocephalides felis) through horizontal transmission. In Sri Lanka, all the three rickettsiosis groups were reported in dogs with a high exposure rates (seroprevalence) at 49% [66]. The information presented here, demonstrates that dogs may play a very important roles in epidemiology and transmission cycle of human rickettsioses. For surveillance and disease ecology aspects, dogs could be used as a sentinel animal to identify high disease burden foci to mitigate risk of rickettsiosis outbreaks in the future.
In addition to dogs, the numbers of cattle and poultry were also positively correlated with levels of human TG and SFG exposures, respectively. Normally, evidence of R. typhi exposure in companion pets (e.g., dogs and cats) or in wildlife (e.g., rodents and opossums) translates to increased risk of human infections [22]. To our knowledge, cattle has never been linked directly with transmission cycle of R. typhi. However, previous publications revealed potential occurrences of flea infestation on cattle; e.g., Ctenocephalides felis felis was found on buffalo calves in India [67] and on dairy calves in Australia [68]; Pulex irritans infested on cow and sheep in Iran [69]. In terms of poultry, several ectoparasite species such as biting lice (e.g., Lipeurus caponis, Goniodes dissimilis, G. gallinae, Menacanthus stramineus and Menopon gallinae), fleas (e.g., Echidnophaga gallinacea), ticks (e.g., Haemaphysalis wellingtoni) and mites (e.g., Ornithonyssus bursa, Dermanyssus gallinae, Megninia cubitalis and Pterolichus obtusus) parasitize domestic chickens in Thailand [70,71], and might be play a vectorial role for SFG rickettsial transmission. However, the study of Rickettsia-associated poultry is very limited in Southeast Asia. In the other parts of the world, information regarding rickettsial infections either in poultry or in ectoparasites associated with avian species have been documented but only in a limited manner. Antibodies against Rickettsia spp. (i.e., R. rickettsiii, R. parkeri and/or R. bellii) were detected in domestic chickens from Brazil [72]. Rickettsia sp. (closely related to R. africae) were detected in flea specimens, E. gallinacea collected from domestic chickens in Madagascar [73]. SFG rickettsiae (e.g., R. aeschlimannii, R. africae, R. helvetica, R. massiliae, R. monacensis and R. slovaca) were isolated from ticks of migratory birds in European countries [74–76]. Furthermore, R. helvetica and R. monacensis were detected in ticks removed from birds in Taiwan, alongside with other tick-borne pathogens such as Borrelia turdi, Anaplasma sp. and Ehrlichia sp [77]. These results highlight a prominent relationship among avian hosts and their ectoparasites in carrying and maintaining SFG rickettsiae in nature.
Further studies on ectoparasite burden in dogs, cattle, poultry and wild rodents as potential reservoirs of several ectoparasites of public health and veterinary importance. These studies should be combined with investigations of rickettsial infections in the animal hosts and their associated arthropod vectors in contrasting ecologies to better determine the role of these animals in zoonotic cycle of rickettsial transmission.
Supporting information
S1 Table. Demographic data of each village in Saen Thong Sub-district, Tha Wang Pha District, Nan Province, Thailand.
https://doi.org/10.1371/journal.pntd.0010256.s001
(DOCX)
S1 Fig. Mapping of scrub typhus group (STG) seropositive cases by villages in Saenthong subdistrict, Nan Province.
Green and red circles indicate sero-negativity and sero-positivity of STG, respectively. Size of the red circles indicate level of STG positivity (IFA titers).
https://doi.org/10.1371/journal.pntd.0010256.s002
(TIF)
S2 Fig. Mapping of typhus group (TG) seropositive cases by villages in Saenthong subdistrict, Nan Province.
Green and red circles indicate sero-negativity and sero-positivity of TG, respectively. Size of the red circles indicate level of TG positivity (IFA titers).
https://doi.org/10.1371/journal.pntd.0010256.s003
(TIF)
S3 Fig. Mapping of spotted fever group (SFG) seropositive cases by villages in Saenthong subdistrict, Nan Province.
Green and red circles indicate sero-negativity and sero-positivity of SFG, respectively. Size of the red circles indicate level of SFG positivity (IFA titers).
https://doi.org/10.1371/journal.pntd.0010256.s004
(TIF)
Acknowledgments
We would like to express our gratitude to the local administration in Nan Province, including Nan Provincial Public Health Office, Tha Wang Pha District Public Health Office, Tha Wang Pha Hospital and Saen Thong Subdistrict Health Promoting Hospital for supporting and facilitating research in the field. Also, we would like to extend our sincere thanks to the community leaders and the villagers for their involvement in this study. CT, AK and SM thank the Thailand International Cooperation Agency (TICA) with the supporting project “Innovative Animal Health”.
References
- 1. Robinson MT, Satjanadumrong J, Hughes T, Stenos J, Blacksell SD. Diagnosis of spotted fever group Rickettsia infections: the Asian perspective. Epidemiol Infect. 2019 Oct;147:e286. pmid:31587667
- 2. Low VL, Tan TK, Khoo JJ, Lim FS, AbuBakar S. An overview of rickettsiae in Southeast Asia: Vector-animal-human interface. Acta Trop. 2020 Feb;202:105282. pmid:31778642
- 3. Morand S, Chaisiri K, Kritiyakan A, Kumlert R. Disease Ecology of Rickettsial Species: A Data Science Approach. Trop Med Infect Dis. 2020 Apr;5(2). pmid:32349270
- 4. Ng-Nguyen D, Hii S-F, Hoang M-TT, Nguyen V-AT, Rees R, Stenos J, et al. Domestic dogs are mammalian reservoirs for the emerging zoonosis flea-borne spotted fever, caused by Rickettsia felis. Sci Rep [Internet]. 2020;10(1):4151. Available from: pmid:32139802
- 5. Suputtamongkol Y, Suttinont C, Niwatayakul K, Hoontrakul S, Limpaiboon R, Chierakul W, et al. Epidemiology and clinical aspects of rickettsioses in Thailand. Ann N Y Acad Sci. 2009 May;1166:172–9. pmid:19538278
- 6. Acestor N, Cooksey R, Newton PN, Ménard D, Guerin PJ, Nakagawa J, et al. Mapping the aetiology of non-malarial febrile illness in Southeast Asia through a systematic review—terra incognita impairing treatment policies. PLoS One. 2012;7(9):e44269. pmid:22970193
- 7. Elliott I, Pearson I, Dahal P, Thomas N V, Roberts T, Newton PN. Scrub typhus ecology: a systematic review of Orientia in vectors and hosts. Parasit Vectors. 2019 Nov;12(1):513. pmid:31685019
- 8. Horton KC, Jiang J, Maina A, Dueger E, Zayed A, Ahmed AA, et al. Evidence of Rickettsia and Orientia Infections Among Abattoir Workers in Djibouti. Am J Trop Med Hyg. 2016 Aug;95(2):462–5. pmid:27273647
- 9. Izzard L, Fuller A, Blacksell SD, Paris DH, Richards AL, Aukkanit N, et al. Isolation of a novel Orientia species (O. chuto sp. nov.) from a patient infected in Dubai. J Clin Microbiol. 2010 Dec;48(12):4404–9. pmid:20926708
- 10. Weitzel T, Martínez-Valdebenito C, Acosta-Jamett G, Jiang J, Richards AL, Abarca K. Scrub Typhus in Continental Chile, 2016-2018(1). Emerg Infect Dis. 2019 Jun;25(6):1214–7. pmid:30835200
- 11. Abarca K, Martínez-Valdebenito C, Angulo J, Jiang J, Farris CM, Richards AL, et al. Molecular Description of a Novel Orientia Species Causing Scrub Typhus in Chile. Emerg Infect Dis. 2020 Sep;26(9):2148–56. pmid:32818400
- 12. Vallée J, Thaojaikong T, Moore CE, Phetsouvanh R, Richards AL, Souris M, et al. Contrasting spatial distribution and risk factors for past infection with scrub typhus and murine typhus in Vientiane City, Lao PDR. PLoS Negl Trop Dis. 2010 Dec;4(12):e909. pmid:21151880
- 13. Cadigan FCJ, Andre RG, Bolton M, Gan E, Walker JS. The effect of habitat on the prevalence of human scrub typhus in Malaysia. Trans R Soc Trop Med Hyg. 1972;66(4):582–7. pmid:4561007
- 14. Harris PNA, Oltvolgyi C, Islam A, Hussain-Yusuf H, Loewenthal MR, Vincent G, et al. An outbreak of scrub typhus in military personnel despite protocols for antibiotic prophylaxis: doxycycline resistance excluded by a quantitative PCR-based susceptibility assay. Microbes Infect. 2016 Jun;18(6):406–11. pmid:27005452
- 15. Weitzel T, Aylwin M, Martínez-Valdebenito C, Acosta-Jamett G, Abarca K. Scrub typhus in Tierra del Fuego: a tropical rickettsiosis in a subantarctic region. Clinical microbiology and infection : the official publication of the European Society of Clinical Microbiology and Infectious Diseases. England; 2020.
- 16. Peniche Lara G, Dzul-Rosado KR, Zavala Velázquez JE, Zavala-Castro J. Murine Typhus: Clinical and epidemiological aspects. Colomb medica (Cali, Colomb. 2012 Apr;43(2):175–80. pmid:24893060
- 17. Rauch J, Eisermann P, Noack B, Mehlhoop U, Muntau B, Schäfer J, et al. Typhus Group Rickettsiosis, Germany, 2010-2017(1). Emerg Infect Dis. 2018 Jul;24(7):1213–20. pmid:29912688
- 18. Bechah Y, Capo C, Mege J-L, Raoult D. Epidemic typhus. Lancet Infect Dis. 2008 Jul;8(7):417–26. pmid:18582834
- 19.
Walker DH, Ismail N, Olano JP, Valbuena G, McBride J. Pathogenesis, immunity, pathology, and pathophysiology in rickettsial diseases. In: Rickettsial Diseases. New York: Informa Healthcare; 2007. p. 15–26.
- 20. Parola P, Paddock CD, Socolovschi C, Labruna MB, Mediannikov O, Kernif T, et al. Update on tick-borne rickettsioses around the world: a geographic approach. Clin Microbiol Rev. 2013 Oct;26(4):657–702. pmid:24092850
- 21. Akram SM, Jamil RT, Gossman WG. Rickettsia Akari. In Treasure Island (FL); 2021.
- 22. Brown LD, Macaluso KR. Rickettsia felis, an Emerging Flea-Borne Rickettsiosis. Curr Trop Med reports. 2016;3:27–39. pmid:27340613
- 23. Abdad MY, Abou Abdallah R, Fournier P-E, Stenos J, Vasoo S. A Concise Review of the Epidemiology and Diagnostics of Rickettsioses: Rickettsia and Orientia spp. J Clin Microbiol. 2018 Aug;56(8). pmid:29769278
- 24. Paddock CD. Perspectives on the laboratory diagnosis of rickettsial diseases in the 21st century. Vol. 55, Acta Médica Costarricense. scielo; 2013. p. 13–24.
- 25. Um J, Nam Y, Lim JN, Kim M, An Y, Hwang SH, et al. Seroprevalence of scrub typhus, murine typhus and spotted fever groups in North Korean refugees. Int J Infect Dis [Internet]. 2021;106:23–8. Available from: https://www.sciencedirect.com/science/article/pii/S120197122100196X pmid:33676004
- 26. Elders PND, Swe MMM, Phyo AP, McLean ARD, Lin HN, Soe K, et al. Serological evidence indicates widespread distribution of rickettsioses in Myanmar. Int J Infect Dis [Internet]. 2021;103:494–501. Available from: pmid:33310022
- 27. Kim Y-J, Yeo S-J, Park S-J, Woo Y-J, Kim M-W, Kim S-H, et al. Improvement of the diagnostic sensitivity of scrub typhus using a mixture of recombinant antigens derived from Orientia tsutsugamushi serotypes. J Korean Med Sci. 2013 May;28(5):672–9. pmid:23678257
- 28. Paris DH, Dumler JS. State of the art of diagnosis of rickettsial diseases: the use of blood specimens for diagnosis of scrub typhus, spotted fever group rickettsiosis, and murine typhus. Curr Opin Infect Dis. 2016 Oct;29(5):433–9. pmid:27429138
- 29. Elders PND, Dhawan S, Tanganuchitcharnchai A, Phommasone K, Chansamouth V, Day NPJ, et al. Diagnostic accuracy of an in-house Scrub Typhus enzyme linked immunoassay for the detection of IgM and IgG antibodies in Laos. PLoS Negl Trop Dis. 2020 Dec;14(12):e0008858. pmid:33284807
- 30. Blacksell SD, Bryant NJ, Paris DH, Doust JA, Sakoda Y, Day NPJ. Scrub typhus serologic testing with the indirect immunofluorescence method as a diagnostic gold standard: a lack of consensus leads to a lot of confusion. Clin Infect Dis. 2007 Feb;44(3):391–401. pmid:17205447
- 31. Luce-Fedrow A, Mullins K, Kostik AP, St John HK, Jiang J, Richards AL. Strategies for detecting rickettsiae and diagnosing rickettsial diseases. Future Microbiol. 2015;10(4):537–64. pmid:25865193
- 32. Parola P, Miller RS, McDaniel P, Telford SR 3rd, Rolain J-M, Wongsrichanalai C, et al. Emerging rickettsioses of the Thai-Myanmar border. Emerg Infect Dis. 2003 May;9(5):592–5. pmid:12737744
- 33. Bhengsri S, Baggett HC, Edouard S, Dowell SF, Dasch GA, Fisk TL, et al. Sennetsu Neorickettsiosis, Spotted Fever Group, and Typhus Group Rickettsioses in Three Provinces in Thailand. Am J Trop Med Hyg. 2016 Jul;95(1):43–9. pmid:27139448
- 34. Wangrangsimakul T, Elliott I, Nedsuwan S, Kumlert R, Hinjoy S, Chaisiri K, et al. The estimated burden of scrub typhus in Thailand from national surveillance data (2003–2018). PLoS Negl Trop Dis. 2020 Apr;14(4):e0008233. pmid:32287307
- 35. Takada N, Fujita H, Yano Y, Huang WH, Khamboonruang C. Serosurveys of spotted fever and murine typhus in local residents of Taiwan and Thailand compared with Japan. Southeast Asian J Trop Med Public Health. 1993 Jun;24(2):354–6. pmid:8266242
- 36. Sankasuwan V, Pongpradit P, Bodhidatta P, Thonglongya K, Winter PE. Murine typhus in Thailand. Trans R Soc Trop Med Hyg [Internet]. 1969 Jan 1;63(5):639–43. Available from: pmid:5824285
- 37. Duffy PE, Guillouzic H Le, Gass RF, Innis BL. Murine Typhus Identified as a Major cause of Febrile Illness in a Camp for Displaced Khmers in Thailand. Am J Trop Med Hyg [Internet]. 43(5):520–6. Available from: https://www.ajtmh.org/view/journals/tpmd/43/5/article-p520.xml
- 38. Silpapojakul K, Woodtayagone J, Lekakula A, Vimuktalaba A, Krisanapan S. Murine typhus in southern Thailand. J Med Assoc Thai. 1987 Feb;70(2):55–62. pmid:3585210
- 39. Sirisanthana T, Pinyopornpanit V, Sirisanthana V, Strickman D, Kelly DJ, Dasch GA. First cases of spotted fever group rickettsiosis in Thailand. Am J Trop Med Hyg. 1994 Jun;50(6):682–6. pmid:8024059
- 40. Sangkasuwan V, Chatyingmongkol T, Sukwit S, Eamsila C, Chuenchitra T, Rodkvamtook W, et al. Description of the first reported human case of spotted fever group rickettsiosis in urban Bangkok. Am J Trop Med Hyg. 2007 Nov;77(5):891–2. pmid:17984348
- 41. Fournier P-E, Allombert C, Supputamongkol Y, Caruso G, Brouqui P, Raoult D. Aneruptive fever associated with antibodies to Rickettsia helvetica in Europe and Thailand. J Clin Microbiol. 2004 Feb;42(2):816–8. pmid:14766859
- 42. Gaywee J, Sunyakumthorn P, Rodkvamtook W, Ruang-areerate T, Mason CJ, Sirisopana N. Human infection with Rickettsia sp. related to R. japonica, Thailand. Vol. 13, Emerging infectious diseases. 2007. p. 657–9.
- 43. Takada N, Fujita H, Kawabata H, Ando S, Sakata A, Takano A, et al. Spotted fever group Rickettsia sp. closely related to Rickettsia japonica, Thailand. Vol. 15, Emerging infectious diseases. 2009. p. 610–1. pmid:19331747
- 44. Bonell A, Lubell Y, Newton PN, Crump JA, Paris DH. Estimating the burden of scrub typhus: A systematic review. PLoS Negl Trop Dis. 2017 Sep;11(9):e0005838. pmid:28945755
- 45. Phanichkrivalkosil M, Tanganuchitcharnchai A, Jintaworn S, Kantipong P, Laongnualpanich A, Chierakul W, et al. Determination of Optimal Diagnostic Cut-Offs for the Naval Medical Research Center Scrub Typhus IgM ELISA in Chiang Rai, Thailand. Am J Trop Med Hyg. 2019 May;100(5):1134–40. pmid:30860022
- 46. Paris DH, Blacksell SD, Nawtaisong P, Jenjaroen K, Teeraratkul A, Chierakul W, et al. Diagnostic accuracy of a loop-mediated isothermal PCR assay for detection of Orientia tsutsugamushi during acute Scrub Typhus infection. PLoS Negl Trop Dis. 2011 Sep;5(9):e1307. pmid:21931873
- 47.
Wickham H. ggplot2: Elegant Graphics for Data Analysis [Internet]. New York: Springer-Verlag; 2016. Available from: https://ggplot2.tidyverse.org
- 48.
R Core Team. R: A Language and Environment for Statistical Computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2020. Available from: https://www.r-project.org/
- 49. Tennekes M. tmap: Thematic Maps in R. J Stat Software; Vol 1, Issue 6 [Internet]. 2018 Apr 20; Available from: https://www.jstatsoft.org/v084/i06
- 50. Hijmans RJ, van Etten J. raster: Geographic analysis and modeling with raster data. R package; 2012.
- 51. Delignette-Muller ML, Dutang C. fitdistrplus: An R Package for Fitting Distributions. J Stat Software; Vol 1, Issue 4 [Internet]. 2015 Mar 20; Available from: https://www.jstatsoft.org/v064/i04
- 52. Calcagno V, de Mazancourt C. glmulti: An R Package for Easy Automated Model Selection with (Generalized) Linear Models. J Stat Software; Vol 1, Issue 12 [Internet]. 2010 May 31; Available from: https://www.jstatsoft.org/v034/i12
- 53.
Burnham KP, Anderson DR. Model Selection and Multimodel Inference: a Practical Information-Theoretic Approach. 2nd ed. New York: Springer International Publishing; 2002. 488 p.
- 54.
Fox J, Weisberg S. An {R} Companion to Applied Regression [Internet]. 2nd ed. California; 2011. Available from: http://socserv.socsci.mcmaster.ca/jfox/Books/Companion
- 55. Lüdecke D. sjstats: Statistical Functions for Regression Models (Version 0.18.1) [Internet]. 2021. Available from: https://cran.r-project.org/package=sjstats
- 56. Curtin J. lmSupport: Support for Linear Models [Internet]. R package version 2.9.13; 2018. Available from: https://cran.r-project.org/package=lmSupport
- 57. Chanyasanha C, Kaeburong K, Chenchittikul M, Sujirarat D. Seroprevalence of scrub typhus infection in patients with pyrexia at some malaria clinics in three western provinces of Thailand. Asian Pacific J allergy Immunol. 1998;16(2–3):119–25. pmid:9876950
- 58. Strickman D, Tanskul P, Eamsila C, Kelly DJ. Prevalence of antibodies to rickettsiae in the human population of suburban Bangkok. Am J Trop Med Hyg. 1994 Aug;51(2):149–53. pmid:8074248
- 59. Blacksell SD, Kantipong P, Watthanaworawit W, Turner C, Tanganuchitcharnchai A, Jintawon S, et al. Underrecognized arthropod-borne and zoonotic pathogens in northern and northwestern Thailand: serological evidence and opportunities for awareness. Vector Borne Zoonotic Dis. 2015 May;15(5):285–90. pmid:25988437
- 60. Wächter M, Wölfel S, Pfeffer M, Dobler G, Kohn B, Moritz A, et al. Serological differentiation of antibodies against Rickettsia helvetica, R. raoultii, R. slovaca, R. monacensis and R. felis in dogs from Germany by a micro-immunofluorescent antibody test. Parasit Vectors. 2015 Mar;8:126. pmid:25889200
- 61. Blanton LS, Walker DH. Flea-Borne Rickettsioses and Rickettsiae. Am J Trop Med Hyg. 2017 Jan;96(1):53–6. pmid:27799640
- 62. Noden BH, Davidson S, Smith JL, Williams F. First Detection of Rickettsia typhi and Rickettsia felis in Fleas Collected From Client-Owned Companion Animals in the Southern Great Plains. J Med Entomol. 2017 Jul;54(4):1093–7. pmid:28399210
- 63. Tay ST, Ho TM, Rohani MY, Devi S. Antibodies to Orientia tsutsugamushi, Rickettsia typhi and spotted fever group rickettsiae among febrile patients in rural areas of Malaysia. Trans R Soc Trop Med Hyg. 2000;94(3):280–4. pmid:10974999
- 64. Perri AR, Feuerborn TR, Frantz LAF, Larson G, Malhi RS, Meltzer DJ, et al. Dog domestication and the dual dispersal of people and dogs into the Americas. Proc Natl Acad Sci U S A. 2021 Feb;118(6). pmid:33495362
- 65. Hii SF, Kopp SR, Abdad MY, Thompson MF, O’Leary CA, Rees RL, et al. Molecular evidence supports the role of dogs as potential reservoirs for Rickettsia felis. Vector Borne Zoonotic Dis. 2011 Aug;11(8):1007–12. pmid:21612534
- 66. Nanayakkara DM, Rajapakse RPVJ, Wickramasinghe S, Kularatne SAM. Serological evidence for exposure of dogs to Rickettsia conorii, Rickettsia typhi, and Orientia tsutsugamushi in Sri Lanka. Vector Borne Zoonotic Dis. 2013 Aug;13(8):545–9. pmid:23930973
- 67. Singh NK, Haque M, Jyoti , Rath SS, Ghosh S. First report of Ctenocephalides felis felis infestation of buffalo calves in Punjab, India. J Parasit Dis. 2011 Oct;35(2):235–6. pmid:23024513
- 68. Tsoi FMB, Šlapeta J, Reynolds M. Ctenocephalides felis (cat flea) infestation in neonatal dairy calves managed with deltamethrin pour-on in Australia. Vet Parasitol. 2020 Mar;279:109039. pmid:32044501
- 69. Ramezani Awal Riabi H, Atarodi A. A Survey on Fauna of Fleas (Order: Siphonaptera) of Cow and Sheep in Southern Khorasan-e-Razavi Province, Iran. Zahedan J Res Med Sci [Internet]. 2015;17(6). Available from: https://sites.kowsarpub.com/zjrms/articles/998.html
- 70. Changbunjong T, Buddhirongawatr R, Suwanpakdee S, Siengsanan J, Yongyuttawichai P, Cheewajorn K, et al. A survey of ectoparasitic arthropods on domestic animals in Tak Province, Thailand. Southeast Asian J Trop Med Public Health. 2009 May;40(3):435–42. pmid:19842427
- 71. Sangvaranond A. Ectoparasites of Domesticated Chickens in Thailand (in Thai). Kasetsart Vet [Internet]. 2003;13(2):42–54. Available from: http://jkv.vet.ku.ac.th/new/images/journal/13-2/05.pdf
- 72. Maciel JF, da Silva Krawczak F, Sobotyk de Oliveira C, Moraes-Filho J, Labruna MB, Botton S de A, et al. Detection of anti: Rickettsia spp. antibodies in domestic chickens of extensive breeding in an endemic area for spotted fever in the state of Rio Grande do Sul, Brazil. Ciência Rural. 2013;43(11):2037–41.
- 73. Ehlers J, Krüger A, Rakotondranary SJ, Ratovonamana RY, Poppert S, Ganzhorn JU, et al. Molecular detection of Rickettsia spp., Borrelia spp., Bartonella spp. and Yersinia pestis in ectoparasites of endemic and domestic animals in southwest Madagascar. Acta Trop [Internet]. 2020;205:105339. Available from: https://www.sciencedirect.com/science/article/pii/S0001706X19315700
- 74. Mărcuţan I-D, Kalmár Z, Ionică AM, D’Amico G, Mihalca AD, Vasile C, et al. Spotted fever group rickettsiae in ticks of migratory birds in Romania. Parasit Vectors. 2016 May;9(1):294. pmid:27207258
- 75. Elfving K, Olsen B, Bergström S, Waldenström J, Lundkvist A, Sjöstedt A, et al. Dissemination of spotted fever rickettsia agents in Europe by migrating birds. PLoS One. 2010 Jan;5(1):e8572. pmid:20052286
- 76. Wallménius K, Barboutis C, Fransson T, Jaenson TGT, Lindgren P-E, Nyström F, et al. Spotted fever Rickettsia species in Hyalomma and Ixodes ticks infesting migratory birds in the European Mediterranean area. Parasit Vectors [Internet]. 2014;7(1):318. Available from: pmid:25011617
- 77. Kuo C-C, Lin Y-F, Yao C-T, Shih H-C, Chung L-H, Liao H-C, et al. Tick-borne pathogens in ticks collected from birds in Taiwan. Parasit Vectors. 2017 Nov;10(1):587. pmid:29178908