Skip to main content

Potential distribution of Leptotrombidium scutellare in Yunnan and Sichuan Provinces, China, and its association with mite-borne disease transmission



Leptotrombidium scutellare is one of the six main vectors of scrub typhus in China and is a putative vector of hemorrhagic fever with renal syndrome (HFRS). This mite constitutes a large portion of the chigger mite community in southwest China. Although empirical data on its distribution are available for several investigated sites, knowledge of the species’ association with human well-being and involvement in the prevalence of mite-borne diseases remains scarce.


Occurrence data on the chigger mite were obtained from 21 years (2001–2021) of field sampling. Using boosted regression tree (BRT) ecological models based on climate, land cover and elevation variables, we predicted the environmental suitability for L. scutellare in Yunnan and Sichuan Provinces. The potential distribution range and shifts in the study area for near-current and future scenarios were mapped and the scale of L. scutellare interacting with human activities was evaluated. We tested the explanatory power of the occurrence probability of L. scutellare on incidences of mite-borne diseases.


Elevation and climate factors were the most important factors contributing to the prediction of the occurrence pattern of L. scutellare. The most suitable habitats for this mite species were mainly concentrated around high-elevation areas, with predictions for the future showing a trend towards a reduction. Human activity was negatively correlated with the environmental suitability of L. scutellare. The occurrence probability of L. scutellare in Yunnan Province had a strong explanatory power on the epidemic pattern of HFRS but not scrub typhus.


Our results emphasize the exposure risks introduced by L. scutellare in the high-elevation areas of southwest China. Climate change may lead to a range contraction of this species towards areas of higher elevation and lessen the associated exposure risk. A comprehensive understanding of the transmission risk requires more surveillance efforts.

Graphical Abstract


Mite-borne diseases represent a significant proportion of global infectious diseases [1,2,3]. Reports in recent decades have highlighted an increasing number and range expansion of infectious diseases, such as scrub typhus and hemorrhagic fever with renal syndrome (HFRS) [4,5,6]. Scrub typhus (also known as tsutsugamushi disease) is caused by the agent Orientia tsutsugamushi, a typical mite-borne bacterial pathogen causing zoonotic disease. It is well established that the larvae of some chigger mites (Actinotrichida: Trombiculidae) serve as the exclusive vector of scrub typhus, which can be transmitted from rodents and other wild animals to humans during the parasitic phase through the feeding canal, i.e. a stylostome [7, 8]. HFRS is a viral zoonosis caused by members of the genus Hantavirus, and some chigger species are suspected to be possible vectors of the disease [9].

Six chigger species belonging to the genus Leptotrombidium, including L. scutellare (Nagayo et al., 1921), L. deliense (Walch, 1922), L. rubellum Wang and Liao, 1984, L. wenense Wu, Wen, Yang and Wu, 1982, L. insulare Wei, Wang and Tong, 1989 and L. sialkotense Vercammen-Grandjean and Langston, 1976, have been confirmed as the main vectors of scrub typhus in mainland China [10]. These mite species are known to occur in Yunnan Province and the adjacent areas, with L. scutellare and L. deliense being the dominant species [11,12,13,14]. Leptotrobidium scutellare is considered to be one of the possible vectors of Hantavirus, an emerging zoonotic disease [3, 15,16,17]. No longer than a decade ago, L. scutellare was mostly reported in northern China where it was considered to be the main vector of scrub typhus, occurring at latitudes higher than 31° N [8, 18]. However, field observations in the last decade demonstrated an extensive distribution of L. scutellare in the provinces of southwest China [19, 20], effectively challenging the notion that L. deliense was the main vector of scrub typhus in the region [18]. The outcomes of these observations prompted more efforts in studying the distribution and epidemic relevance of L. scutellare. A better understanding of the distribution of L. scutellare may provide knowledge contribution to the prevention of associated diseases and benefit public health interventions.

Recently, considerable efforts have been made to better understand the distribution of L. scutellare and its role in disease transmission. It has been reported that the presence of L. scutellare and the epidemic incidence of scrub typhus were identical in South Korea [21]. These results suggest that changes in the distribution of the vector may affect the corresponding epidemic patterns. A recent study outlined the potential distribution range under simulated future environments of the two most widely distributed vector chiggers L. deliense and L. scutellare in mainland China [22]. Another study mapping the distributions of blood-sucking mites and mite-borne agents also indicated that L. scutellare harbors the highest variety of pathogens [23], which is turn increases the exposure risk to relevant diseases for the human population where the vector species occurs.

Although these earlier studies contributed significant data towards improving our knowledge of the spatial and temporal dynamics of chigger vectors, they relied on the presence data of the mites with coordinates extracted from center points of the corresponding administrative areas, i.e., the counties [24]. Thus, these projections may underestimate the heterogeneity of environmental suitability for the species within the areas, especially when field observations have shown that mites such as L. scutellare showed specialized habitat preference, i.e. uncultivated habitats at mid to high elevations [12, 25]. Taking environmental differences within the sampled areas into account when evaluating the distribution may improve the performance of ecological models. Such an attempt is realistic only when the collection site is accurately documented. So far, the occurrence of L. scutellare in southwest China has been reported at the county level [19, 24]. The suitable area and potential distribution of L. scutellare in southwest China have not been exhaustively discussed, and the interactions between L. scutellare and anthropic activities, such as cropping, logging or tourism, are unexplored. Also, the epidemic role of L. scutellare in mite-borne disease prevalence in southwest China is unclear.

Based on the occurrence data of chigger mites collected continuously in 21 years of field investigations, we explored three main questions: (i) what is the potential distribution range of L. scutellare in the study area under near current and future climate scenarios; (ii) how much does the distribution range of L. scutellare overlap with human activities; and (iii) does the potential distribution range of L. scutellare explain the prevalence of scrub typhus and HFRS in the study region (see the workflow in Fig. 1). Through these efforts, we aimed to evaluate the exposure risk to scrub typhus and HFRS from L. scutellare in Yunnan and Sichuan Provinces on both spatial and temporal scales.

Fig. 1
figure 1

Workflow of the present study


Study area

The study was conducted in Yunnan and Sichuan Provinces, two adjacent provinces in southwest China. The uplift of the Qinghai-Tibet Plateau resulted in the formation of various mountain ranges in this region and, consequently, an extremely high heterogeneity in its topography and climate (Fig. 2). Elevation generally increases from the southeast to northwest, with the mean annual temperature decreasing in the same direction. Moisture delivered by the seasonal monsoon causes a decline in precipitation from west to east. Based on these variations, the study area covers almost all types of existing climates worldwide, from wet or dry tropical to alpine tundra [26]. The western part of the study region is dominated by a spectacular mountain-valley system with vast elevational differences (between 76 and 7556 m a.s.l.), which is known as the Longitudinal Range and Gorge Region of southwest China [27]. Various vegetation types are distributed among this mountain-valley system, including humid evergreen broad-leaved forest and semi-arid bushes along the valleys, mixed broadleaf-conifer forest and dark coniferous forest on mid- and high-elevation slopes, and meadows and tundra in alpine areas. Eastern Sichuan is marked by a low-elevation basin (250–700 m a.s.l.) surrounded by mountains (1000–3000 m a.s.l.). This basin occupies approximately 33% of the area of Sichuan Province and is mostly covered by subtropical evergreen broad-leaved forest except for the areas with human modifications. The western Sichuan area is formed by a vast plateau with an average elevation of > 4000 m a.s.l., which is predominantly covered by oak bushes, meadows and scree. The dark coniferous forest can also be seen in a few valleys in this area. Eastern and middle Yunnan Province is generally dominated by mountainous highlands (approx. 2000 m a.s.l.). Mixed broadleaf-conifer forests are widely distributed in these highlands. Low-elevation plains and hills predominantly constitute the topography in southern Yunnan (< 1000 m a.s.l.), which is generally covered by tropical rainforest and humid evergreen broad-leaved forest. Human settlements and croplands are spread over the valleys, plateaus and low-elevation plains in the study region.

Fig. 2
figure 2

Map showing the topography of the study area and sampling points investigated for the presence/absence of Leptotrombidium scutellare

Field sampling and occurrence data of L. scutellare

Field sampling was conducted in 149 sites across 69 counties in Yunnan and Sichuan Provinces from 2001 to 2021. Animal hosts were collected in cage traps (18 × 12 × 9 cm; Guixi Mousetrap Apparatus Factory, Guixi, Jiangxi, China) baited with fresh peanuts, corn or other standard baits [28]. Chiggers were collected from the body surface of each host separately, following a standardized protocol [29]. The collected chiggers were mounted on glass slides in Hoyer’s medium and identified under an optical microscope. The specimen vouchers are deposited in the specimen repository of the Institute of Pathogens and Vectors, Dali University, China. A total of 178,942 chiggers were collected and subsequently identified as 302 species, with 1053 individuals of L. scutellare determined to be present in 40 of the 149 investigated sampling sites (Fig. 2).

Environmental variables

The distribution range and population persistence of chigger mites have been found to be determined by climate factors and landscape [7, 12, 29, 30]. However, there is no consensus on which climate factor is the most important. To estimate the species-environment relationship and map the potential distribution of L. scutellare, layers of 19 bioclimatic variables relating to temperature and precipitation measures spanning 1970–2000 were downloaded from the Worldclim 2.1 web server (for full description for the 19 bio-climate variables, see The ‘select07’ function in the ‘mecofun’ package was then used to identify the most important climate variables and reduce multi-collinearity [31]. A linear regression model was imposed on these variables to evaluate their explanatory power (adjusted R2) on the occurrence of L. scutellare.

Land cover types (e.g. type of vegetation) were considered a predominant determinant of chigger occurrences in previous studies [12, 32, 33]. We obtained land cover data indexing cropland, shrubland, grassland or forest from a global land use dataset (1-km resolution) [34]. Multiple ecological factors can vary dramatically along the elevational gradient [35, 36]. For L. scutellare, field observation data suggested that elevations between 2000 and 2500 m a.s.l. were the preferred range [12]. Thus, we included elevation as a synthesis variable representing any undefined factors. We used the 30-s resolution elevational layer from the Shuttle Radar Topography Mission downloaded from the Worldclim 2.1 web server.

For future scenarios, we incorporated the projections based on the Coupled Model Intercomparison Project Phase 6 (CMIP6) [37] for climate variables and the Cellular Automata models for land cover change [34]. The CMIP6 replaced the previously used Representative Concentration Pathways (RCPs) with the updated Shared Socio-economic Pathways (SSPs). Four SSPs scenarios (SSP1-2.6, SSP2-4.5, SSP4-6.0, and SSP5-8.5) have levels of radiative forcing that correspond with the RCPs (RCP2.6, RCP4.5, RCP6.0 and RCP8.5) [38]. The BCC-CSM2-MR climate model (a high-resolution climate system model released by the Beijing Climate Center [BCC]) has been found to outperform other general circulation models in certain parameters [39, 40]. Thus, SSP2-4.5 and SSP5-8.5 derived from the model were selected to simulate moderate and high greenhouse gas concentration scenarios, respectively. The layers containing future climate projections (2050, as a climatology computed from 2041–2060; 2090, as a climatology computed from 2081–2100) were also obtained from the Worldclim 2.1 web server. The 30-s resolution land cover layers of the years 2050 and 2090 under the RCP4.5 and RCP8.5 scenarios were obtained from Li et al. [34]. Ultimately, elevation and bioclimatic layers comprising two future periods (2050 and 2090) under moderate (SSP2-4.5 or RCP4.5) and high (SSP5-8.5 or RCP8.5) greenhouse scenarios were used for future distribution projection. All eco-geographic layers were cropped by the extent of Yunnan and Sichuan Provinces, respectively, and resampled to 30 s when necessary for further assessments.

Model fitting and distribution projection

Environmental variables were extracted from the selected layers of the near-current period by the presence-absence points with the ‘extract’ function provided by the ‘terra’ package [41]. We fitted a boosted regression tree (BRT) ecological model based on a case-control research strategy for species distribution modeling [42]. The data points of the investigated sites were assigned as “cases” when L. scutellare is present (1) and “controls” when the species is absent (0). Of all data points, 75% were extracted randomly as a training set, and the remaining 25% were treated as a test set [23]. The training set was used to build a BRT model and to acquire the relative contributions of the variables using the ‘gbm.step’ function provided in the ‘dismo’ package [43]. The model was run with a tree complexity of 5, a learning rate of 0.005 and a bagging fraction of 75% based on the satisfactory performance in previous practice runs [23, 44]. Tenfold cross-validation was adopted to produce the receiver-operating characteristic (ROC) curves. The model performance and predicted power were evaluated by the areas under the curve (AUC). The AUC of the fitted model was derived with the ‘evalSDM’ function available in the ‘mecofun’ package using the test set for validation [31]. The training set extraction, BRT model fitting and predictive performance evaluation were repeated 100 times. We averaged the values provided by the repetitions to increase the robustness of the modeling performance [22].

For each focused scenario and period, the occurrence probability was predicted based on the BRT model products of 100 repetitions with the ‘predict.gbm’ function available in the ‘gbm’ package [45]. The outputs were also averaged for geographical mapping. The ‘evalSDM’ function also returns a cut-off value (threshold) based on the maximum sensitivity+specificity along the ROC curve that can be used to translate the predicted probabilities of L. scutellare into binary levels [31]. Pixels with probabilities beyond the threshold were considered suitable (1), while those below the threshold were considered unsuitable (0) for L. scutellare to map the binarizing occurrence patterns.

Relating occurrence probability with human activities and mite-borne diseases

The Human Footprint index shows the gradient of human activities [46]. It is a weighted sum of maps of where people live (population density) and the accessibility to different areas (e.g. infrastructure). We downloaded the layer mapping the Human Footprint for the year 2020, released by the Wildlife Conservation Society (, and then re-sampled (bilinear interpolation) the Human Footprint layer to 30 s. We calculated the Pearson’s correlation coefficients between the Human Footprint index and the occurrence probability of L. scutellare under near-current and future projections with the ‘raster.cor.matrix’ function provided in the ‘ENMTools’ package [47].

Finally, we fitted linear regression models to evaluate the explanatory power of the occurrence probability of L. scutellare for the mite-borne disease frequencies in recent decades. Because case reports for the mite-borne disease were not publicly available for Sichuan Province, in the present study we only analyzed data on Yunnan Province. An area with more suitable habitats for L. scutellare should have a higher occurrence probability and, consequently, a higher infection rate for mite-borne diseases. We first averaged the predicted probability of L. scutellare at the city level and then derived the epidemiological patterns of mite-borne diseases from published literature. We obtained the case reports of scrub typhus in Yunnan Province during the period 2006–2017 [5]. The case reports of HFRS in Yunnan Province had been aggregated by the periods 1976–2012 [48] and 2014–2015 [49]. The number of case reports was used as a response variable, and the probability of L. scutellare was used as an explanatory factor when fitting the linear regression models. The occurrence data were primarily organized in Microsoft Excel (Microsoft Corp., Redmond, WA, USA), while the eco-geographic layers and statistical analyses were processed in the R 4.1.0 environment [50]. The output maps were assembled in QGIS 3.18 [51].


Model ensembles and predictive performance

During the field investigation, L. scutellare was collected from 34 species of small mammals, with Eothenomys miletus, Apodemus chevrieri, Tupaia belangeri and Neotetracus sinensis being the dominant hosts. The average prevalence of L. scutellare in sites where the species occurred was 12.1% (n = 1053/8728). The variable selection process retained six of the initial set of 19 bioclimatic variables while maintaining most of the explanatory power for L. scutellare occurrence after removal of highly correlated variables. These six variables were temperature mean diurnal range, temperature seasonality, mean temperature of warmest quarter, precipitation of wettest month, precipitation seasonality and precipitation of coldest quarter. Linear regression showed that these six climate variables together explained 18.9% of the presence-absence pattern of L. scutellare in the study area.

The average (± standard deviation [SD]) model AUC was 0.85 ± 0.07, indicating that the ensembled BRT model had a high predictive performance. Elevation (mean ± SD, 24.4 ± 5.7%) yielded the highest contribution to predicting environmental suitability (Table 1), followed by temperature seasonality (22.1 ± 5.5%) and precipitation seasonality (19.1 ± 7.3%). Shrubland contributed the least (0.3 ± 0.4%) of the 11 variables. Model performance evaluation returned a cut-off value of 0.35 ± 0.18 to determine the binary suitability for L. scutellare in the study area. The marginal effects of the 11 environmental variables are visualized in Additional file 1: Figure S1.

Table 1 Relative contributions of the 11 environmental variables used to predict the distribution of Leptotrombidium scutellare ensembled by the boosted regression tree model

Distribution range shifts

Geographical projection suggested that the suitable areas for L. scutellare were mostly spread over the western Yunnan plateau and southwest Yunnan mountain ridges (Fig. 3). Compared with these areas, the plateau in western Sichuan was less suitable. The Sichuan basin and the southern Yunnan lowlands were consistently recognized as unsuitable in the near-current and future periods (Fig. 4). The predicted areas suitable for L. scutellare were 214,272.40 km2 in the near-current period (Additional file 1: Figure S2). In the future simulations of moderate greenhouse gas concentrations, the suitable areas were reduced to 55,916.57 km2 in 2050 and 52,064.95 km2 in 2090. Under high greenhouse gas concentrations, the suitable areas in the study region were further reduced to 51,304.36 km2 in 2050 and 41,217.22 km2 in 2090, respectively (Additional file 1: Figure S3).

Fig. 3
figure 3

Potential suitable areas for Leptotrombidium scutellare occurrence in Yunnan and Sichuan, China in the near-current period

Fig. 4
figure 4

Projections of the suitable areas for Leptotrombidium scutellare occurrence in two future periods (2050 and 2090) under moderate (a, c) and high (b, d) greenhouse gas concentrations

Correlations between L. scutellare and exposure risks to mite-borne diseases

The correlations between the Human Footprints index and the predicted occurrence probabilities of L. scutellare in the near-current and future periods were negative (Table 2), indicating the dispreference of habitats with intensive human activities. The correlation coefficient was weakly negative in the near-current period (Pearson’s r = − 0.159). The negative correlations between the occurrence probabilities of L. scutellare and anthropic activities tended to be enhanced over time as the greenhouse effect increased (Table 2).

Table 2 Pearson’s correlation coefficients between Human Footprints and the projected distribution patterns of Leptotrombidium scutellare in different periods and greenhouse gas concentration scenarios

In the near-current period, the cities with the highest average occurrence probability of L. scutellare were Panzhihua (0.43) in Sichuan Province and Dali (0.64) in Yunnan Province, while Suining (0.05) and Zhaotong (0.15) reach the lowest average occurrence probability in Sichuan and Yunnan Provinces, respectively (Additional file 1: Table S1). The linear regression models showed that the occurrence probabilities of L. scutellare in Yunnan Province had high explanatory power for the incidence of HFRS (adjusted R2 = 0.41, F[14] = 11.32, P = 0.005) but not for the incidence of scrub typhus (adjusted R2 = − 0.07, F[14] = 0.01, P = 0.923) (Fig. 5).

Fig. 5
figure 5

Linear correlations between the potential occurrence probabilities of Leptotrombidium scutellare and mite-borne disease incidences in the cities of Yunnan Province. The blue line shows the linear regression line and the gray area is the 95% confidence interval. Black dots mark the cities used to fit the regression model. HFRS, Hemorrhagic fever with renal syndrome


Potential distribution range in near-current and future periods

In the present study, the results of the BRT models showed that elevation and climate factors explained most of the occurrence pattern of L. scutellare (> 91%), with elevation contributing the strongest to the prediction of the occurrence pattern. This result is in contrast with those of a recent study in which it was found that elevation contributed little (no more than 5%) to the species in their models [22]. This difference can be explained by the fact that our point data contains more accurate information on the presence of L. scutellare along the elevational gradients. The mountainous areas of southwest China are marked by a vast variation in elevation [27], with corresponding variations in environmental factors, such as temperature, water resources, vegetation cover and soil chemistry [36]. Such variations in elevation have been found to impact the distribution of small mammal hosts [52] and their ectoparasitic species [25]. Our results are in line with the previous finding that L. scutellare prefers mid- to high-elevation habitats [12], further highlighting the important influence of topography in the distribution of the species.

Climate variables constituted a considerable portion of the contribution to predicting the occurrence pattern. Findings from previous studies suggested a negative correlation between temperature and L. scutellare presence [53, 54]. Leptotrombidium scutellare is known to be widely distributed in northern China, which is colder in winter compared to the southern regions of the country [24, 55]. The species is also known to be more active during the cold seasons and to act as the main vector of scrub typhus in southern China during winter [18, 56]. Thus, an area with a cold winter in low latitudes may represent optimization of the environmental suitability for L. scutellare. The high contribution of precipitation seasonality to predicting the occurrence pattern of L. scutellare may be attributed to the life history traits as it has been noted that rainfall heavily influences the reproduction and growth of this chigger species, as well as that of many other species [57]. Overall, the results of the present study support the view that the occurrence of L. scutellare is highly sensitive to climatic and topographic heterogeneity [7, 12].

Land cover variables explained only a small portion of the occurrence pattern (< 9%). It has been suggested that there is no consistent effect of land covers such as woodlands, scrublands, hedgerows and grasslands on the abundance of pests and pathogens and that land cover has a small effect compared to climate [58]. Our results support the idea that climate generally outweighs land covers in determining the occurrence of L. scutellare. Previous studies supported a higher prevalence of chiggers in densely vegetated lands [2, 32]. Our results showed that croplands and forests are generally more suitable for L. scutellare than shrublands and grasslands. There are a number of possible explanations for this result. Firstly, forest and cropland constitute most of the land cover area in the study region, with shrublands and grasslands appearing only as small patches amidst them. Secondly, the local people in southwest China have developed unique cropping methods due to the mountainous terrain, and these differ from those used in flat plains; for example, people in southwest China developed the famous terrace farming to cultivate crops alongside mountainous forest habitats. This cropping method may create a transitional zone that is favored by small mammal hosts and their ectoparasitic chiggers. A similar positive correlation has been reported between scrub typhus incidence and Taiwan's mosaic landscapes of cropland and forest [59].

The geographical projection showed that under the near-current period scenario, plateau and high-elevation mountains are generally more suitable than plains and valleys. These results suggest that people living in mountainous areas likely have higher transmission risks related to L. scutellare. The occurrence probabilities of suitable areas for L. scutellare in Yunnan and Sichuan Provinces decreased with increasing greenhouse effect, possibly because this species prefers a colder environment compared to other vector mite species in Yunnan Province [60, 61]. Therefore, under future temperature increases, the range of L. scutellare may contract to high-elevation mountains. Such a tendency was also reported in other mountain-dwelling species, especially those adapted to cold temperatures and high-elevation habitats, for example, the Himalayan coldwater specialist snow trout (Schizothorax richardsonii) [62] and herbages restricted to alpine habitats [63]. However, our results are in contrast with a previous projection in the distribution tendency. A recent study found that both L. scutellare and L. deliense had broader suitability ranges in the high RCP scenario, suggesting a future northward expansion of geographic distribution [22]. This discrepancy may be attributed to their occurrence points not representing the preferable elevation of L. scutellare [24]. Leptotrombidium scutellare and L. deliense can co-exist in a given area yet occupy different elevations and landscapes [12, 19]. Thus, including accurate environmental information on the collection site in the ecological model may improve our understanding of environmental suitability for the species.

Despite the difference in directions, our results are in agreement with previous studies in regarding that the greenhouse effect may optimize habitat suitability in a colder environment, i.e. high elevation or high latitude. Combined, our results of the distribution mapping study indicate that increased temperature may drive the upward or northward shifts for L. scutellare to pursue a more suitable habitat [64, 65]. Moreover, although the potential distribution range of L. scutellare may be reduced in the future, the temperature rise may expand the suitable habitats for other vector mites, such as L. deliense and L. rubellum, which favor warmer conditions. For example, opposite responses to climate change were detected between two small mammal hosts in the mountains of southwest China [66].

Evaluation of the risks introduced by L. scutellare

The results of the present study indicate a low correlation between the potential distribution patterns of L. scutellare and the Human Footprint index in the near-current period. This result is due to the preference of this species for a mountainous uncultivated landscape and because it is primarily parasitic on hosts that live in the wilderness [12, 56]. Studies in South Korea [33] and Thailand [32] suggested that human disturbance and urbanization may mitigate the risk of rickettsioses transmission because the vector chiggers generally prefer habitats less modified by people. Nevertheless, the increasing occurrence probability of L. scutellare may increase the infection risk for people living in high-elevation areas. These people may have a higher opportunity to engage the vector due to their residences and activity patterns being close to the suitable habitats of L. scutellare. On the other hand, L. scutellare was found to peak in population size in cold seasons [7] and was noted to act as the main vector of autumn-winter-type scrub typhus [67]. The forest fire-safe periods in winter and spring, which forbid most human entry into the mountainous areas in Yunnan and Sichuan Provinces, may help to reduce the exposure risk. The interactions between L. scutellare and human activities were expected to decline gradually under future scenarios due to the contraction of suitable habitats. It is worth noting that these results were valid only in a certain context because they were built on the assumption that the intensity of anthropic activities remains constant for the following decades. However, the contraction of L. scutellare range does not necessarily correspond to a low exposure risk [44] because other vectors that coexist in the study region may expand their distribution under climate change scenarios [22, 66]. Therefore, a study on a multi-species level focusing on the study region is needed to improve the evaluation of exposure risk to mite-borne diseases.

Our results showed that the occurrence probability of L. scutellare did not explain the incidence of scrub typhus in Yunnan Province. As such, our findings are not in agreement with those of a study in South Korea that suggested a congruent pattern between the mapped distribution of L. scutellare and the incidence of scrub typhus [21]. A plausible explanation for our results is that the distribution of L. scutellare was mostly restricted to less urbanized areas in southwest China. This leads to a confined influence of the species on the majority of the human population. The incidences of scrub typhus were mostly reported in cities representative of the low-elevation areas of southern Yunnan [5]. These cities were also likely favored by L. deliense—another major O. tsustugamushi vector adapted to warmth and low elevation [61]. However, high suitability for L. scutellare in the plateaus of Yunnan Province may still drive the increasing trend of scrub typhus in these areas, such as in Kunming and Dali [5, 68].

We found strong support for a positive correlation between the occurrence probability of L. scutellare and the incidence of HFRS in Yunnan Province. This result highlights the potential for L. scutellare playing a considerable role in HFRS transmission. The hantavirus has been segregated from multiple small mammal hosts in the wild, including members of the genera Apodemus and Eothenomys [69, 70]. Leptotrombidium scutellare has been reported to be a dominant ectoparasitic species on these hosts [12]. It is worth noting that to date the epidemiological evidence of transmission of HFRS through L. scutellare remained unsatisfactory [56]. Thus, more research may be needed to establish a straightforward association between L. scutellare and HFRS.

Our analysis encountered several unavoidable limitations. Many areas in the study region have not been adequately investigated. Thus, more surveillance efforts in the future are required to increase our knowledge of the actual distribution of L. scutellare and other potential vector mites. The predictions of exposure risk were based on the human activity data available for the contemporary period. This interaction may be exacerbated or mitigated depending on the forest management policy, the progress of urbanization and the development of farming, tourism and infrastructure construction. We could not address the correlations between the occurrence of the chigger with relevant diseases such as scrub typhus and HFRS in Sichuan Province. Further study should target such an issue to better understand the role of potential vectors underlying the epidemic patterns.


Our results showed that people living in high-elevation areas are likely exposed to mite-borne diseases transmitted by L. scutellare. Future climate change is expected to reduce the suitable areas for L. scutellare in Yunnan and Sichuan Provinces. Such range shifts may mitigate the range overlap between humans and L. scutellare. We found that the occurrence probability of the species yielded contrasting associations between different mite-borne diseases, emphasizing the necessity of more study to understand the underlying epidemiological mechanisms. Overall, our findings provide vital insight into prophylactic interventions for improving public health in the targeted region.

Availability of data and materials

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.



Boosted regression tree


Coupled Model Intercomparison Project Phase 6


Hemorrhagic fever with renal syndrome


  1. Luce-Fedrow A, Lehman ML, Kelly DJ, Mullins K, Maina AN, Stewart RL, et al. A review of scrub typhus (Orientia tsutsugamushi and related organisms): then, now, and tomorrow. Trop Med Infect Dis. 2018;3:8.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Traub R, Wisseman CL. The ecology of chigger-borne rickettsiosis (scrub typhus). J Med Entomol. 1974;11:237–303.

    Article  CAS  PubMed  Google Scholar 

  3. Song G. Epidemiological progresses of hemorrhagic fever with renal syndrome in China. Chin Med J. 1999;112:472–7.

    CAS  PubMed  Google Scholar 

  4. Peng PY, Xu L, Wang GX, He WY, Yan TL, Guo XG. Epidemiological characteristics and spatiotemporal distribution of scrub typhus in mainland China in 1952–1989 and 2006–2017. Chin J Zoonoses. 2022;38:818–29 (in Chinese with English abstract).

    Google Scholar 

  5. Yuan Q, Chen M, Yang X. Epidemiological analysis of scrub typhus in Yunnan province during 2006–2017. Chin J Vector Biol Control. 2018;29:69–72 (in Chinese with English abstract).

    Google Scholar 

  6. Zou LX, Sun L. Analysis of hemorrhagic fever with renal syndrome using wavelet tools in mainland China, 2004–2019. Front Public Health. 2020;8:571984.

  7. Lv Y, Guo XG, Jin DC, Song WY, Peng PY, Lin H, et al. Infestation and seasonal fluctuation of chigger mites on the Southeast Asian house rat (Rattus brunneusculus) in southern Yunnan Province China. Int J Parasitol Parasites Wildl. 2021;14:141–9.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Elliott I, Pearson I, Dahal P, Thomas NV, Roberts T, Newton PN. Scrub typhus ecology: a systematic review of Orientia in vectors and hosts. Parasites Vectors. 2019;12:513.

  9. Santibáñez P, Palomar AM, Portillo A, Oteo JA, Santibáñez S. The role of chiggers as human pathogens. In: Sami A, editor. An overview of tropical diseases. London: IntechOpen; 2015. Chapter 9.

  10. Wu GH, Jiang ZK, Wang L, Ding LY, Mao CQ, Ma BY. Accordance and identification of vector chigger mites of tsutsugamushi disease in China. Chin J Hyg Insectic Equip. 2013;19:286–92 (in Chinese with English abstract).

    Google Scholar 

  11. Lv Y, Guo XG, Jin DC, Song WY, Fan R, Zhao CF, et al. Relative abundance of a vector of scrub typhus, Leptotrombidium sialkotense, in southern Yunnan Province China. Korean J Parasitol. 2020;58:153–9.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Peng PY, Guo XG, Jin DC, Dong WG, Qian TJ, Qin F, et al. Landscapes with different biodiversity influence distribution of small mammals and their ectoparasitic chigger mites: a comparative study from southwest China. PLoS ONE. 2018;13:e0189987.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Niu AQ, Guo XG, Men XY, Dong WG, Bao HE. Community structure of chigger mites on small mammals in the surrounding areas of Erhai lake in Dali of Yunnan. Chin J Vector Biol Control. 2005;16:434–7 (in Chinese with English abstract).

    Google Scholar 

  14. Yang GR, Yu ZZ, Gong ZD, Tao KH, Bing FG, Yang H, et al. A study on chigger mites in Xigu Flatland district besides Lancan River. Chin J Vector Biol Control. 1993;4:428–30 (in Chinese with English abstract).

    Google Scholar 

  15. Wu GH, Zhang Y, Guo HB, Jiang KJ, Zhang JJ, Gan YH. The role of Leptotrombidium scutellare in the transmission of human diseases. Chin Med J. 1996;109:670–3.

    CAS  PubMed  Google Scholar 

  16. Zhang Y, Li XF, Zhu J, Tang JQ, Li YX, Wu GH, et al. Preliminary studies on proliferation of hemorrhagic fever with renal syndrom virus in Leptotrombidium (L.) Scutellare. Chin J Exp Clin Virol. 1997;11:259–62 (in Chinese with English abstract).

    CAS  Google Scholar 

  17. Yu J, Deng XZ, Yang ZQ, Yao PP, Zhu HP, Xiong HR, et al. Study on the transmission of Hantaan virus and Orientia tsutsugamushi by naturally dual infected Leptotrombidium scutellare through stinging. Chin J Prev Med. 2010;44:324–8 (in Chinese with English abstract).

    Google Scholar 

  18. Li GC, Li XZ, Liu QY, Lu L, Mu QZ, Wang YJ, et al. Division of epidemic areas for summer- and autumn-type scrub typhus in China. Chin J Vector Biol Control. 2019;30:233–6 (in Chinese with English abstract).

    Google Scholar 

  19. Peng PY, Guo XG, Ren TG, Dong WG, Song WY. An updated distribution and hosts: trombiculid mites (Acari: Trombidiformes) associated with small mammals in Yunnan Province, southwest China. Parasitol Res. 2016;115:1923–38.

    Article  PubMed  Google Scholar 

  20. Peng PY, Guo XG, Song WY, Hou P, Zou YJ, Fan R, et al. Analysis of ectoparasites infestation on Eothenomys miletus in 39 counties of three provinces in Southwest China. Mod Prev Med. 2015;42:3016–21 (in Chinese with English abstract).

    Google Scholar 

  21. Roh JY, Song BG, Park WI, Shin EH, Park C, Park M-Y, et al. Coincidence between geographical distribution of Leptotrombidium scutellare and scrub typhus incidence in South Korea. PLoS ONE. 2014;9:e113193.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Ma T, Hao M, Chen S, Ding F. The current and future risk of spread of Leptotrombidium deliense and Leptotrombidium scutellare in mainland China. Sci Total Environ. 2022;843:156986.

    Article  CAS  PubMed  Google Scholar 

  23. Wang T, Meng F, Che T, Chen J, Zhang H, Ji Y, et al. Mapping the distributions of blood-sucking mites and mite-borne agents in China: a modeling study. Infect Dis Poverty. 2022;11:41.

  24. Meng FF, Xu Q, Chen JJ, Ji Y, Zhang WH, Fan ZW, et al. A dataset of distribution and diversity of blood-sucking mites in China. Sci Data. 2021;8:204.

  25. Peng PY, Guo XG, Ren TG, Song WY. Faunal analysis of chigger mites (Acari: Prostigmata) on small mammals in Yunnan province, southwest China. Parasitol Res. 2015;114:2815–33.

    Article  PubMed  Google Scholar 

  26. Zhang RZ. Physical geography of Hengduan mountains. Beijing: Science Press; 1997 (in Chinese).

  27. Pan T, Wu S, He D, Dai E, Liu Y. Effects of longitudinal range-gorge terrain on the eco-geographical pattern in Southwest China. J Geogr Sci. 2012;22:825–42.

    Article  Google Scholar 

  28. Mohd-Taib FS, Ishak SN. Bait preferences by different small mammal assemblages for effective cage-trapping. Malaysian J Sci. 2021;40:1–15.

    Article  Google Scholar 

  29. Lv Y, Guo XG, Jin DC, Fan R, Zhao CF, Zhang ZW, et al. Distribution and host selection of the chigger mite vector of scrub typhus, Leptotrombidium deliense, in southwest China. Int J Acarol. 2021;47:233–41.

    Article  Google Scholar 

  30. Van Peenen PFD, Lien J-C, Santana FJ, See R. Correlation of chigger abundance with temperature at a hyperendemic focus of scrub typhus. J Parasitol. 1976;62:653–4.

    Article  PubMed  Google Scholar 

  31. Zurell D: mecofun: Useful functions for macroecology and species distribution modelling. R package version 0.1.1. 2020. Accessed 21 Dec 2022.

  32. Linsuwanon P, Auysawasdi N, Wongwairot S, Leepitakrat S, Rodkhamtook W, Wanja E, et al. Assessing scrub typhus and rickettsioses transmission risks in the Chiang Rai province of northern Thailand. Travel Med Infect Dis. 2021;42:102086.

    Article  PubMed  Google Scholar 

  33. Park S-W, Ha N-Y, Ryu B, Bang JH, Song H, Kim Y, et al. Urbanization of scrub typhus disease in South Korea. PLoS Negl Trop Dis. 2015;9:e0003814.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Li XC, Yu L, Sohl T, Clinton N, Li WY, Zhu ZL, et al. A cellular automata downscaling based 1km global land use datasets (2010–2100). Sci Bull. 2016;61:1651–61.

    Article  CAS  Google Scholar 

  35. Barry RG. Mountain weather and climate. 3rd ed. Cambridge: Cambridge University Press; 2008.

    Book  Google Scholar 

  36. Ming QZ, Shi ZT, Deng YJ, Dong M. High gradient effects of mountain regions——taking the effects of nature-human landscape in the Hengduan Mountains for example. J Glaciol Geocryol. 2006;28:925–30.

    Google Scholar 

  37. Eyring V, Bony S, Meehl GA, Senior CA, Stevens B, Stouffer RJ, et al. Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geosci Model Dev. 2016;9:1937–58.

    Article  Google Scholar 

  38. Meinshausen M, Nicholls ZRJ, Lewis J, Gidden MJ, Vogel E, Freund M, et al. The shared socio-economic pathway (SSP) greenhouse gas concentrations and their extensions to 2500. Geosci Model Dev. 2020;13:3571–605.

    Article  CAS  Google Scholar 

  39. Wang D, Liu J, Wang H, Shao W, Mei C, Ding X. Performance evaluations of CMIP6 and CMIP5 models for precipitation simulation over the Hanjiang River Basin. China J Water Clim Chang. 2022;13:2089–106.

    Article  Google Scholar 

  40. Shiru MS, Chung E-S. Performance evaluation of CMIP6 global climate models for selecting models for climate projection over Nigeria. Theor Appl Climatol. 2021;146:599–615.

    Article  Google Scholar 

  41. Hijmans RJ: terra: Spatial data analysis. R package version 1.6–47. 2022. Accessed 02 Dec 2022.

  42. Elith J, Leathwick JR, Hastie T. A working guide to boosted regression trees. J Anim Ecol. 2008;77:802–13.

    Article  CAS  PubMed  Google Scholar 

  43. Hijmans RJ, Phillips S, Leathwick J, Elith J: dismo: species distribution modeling. R package version 1.3–9. 2022. Accessed 19 Sept 2022.

  44. Yao HW, Wang YX, Mi XM, Sun Y, Liu K, Li XL, et al. The scrub typhus in mainland China: spatiotemporal expansion and risk prediction underpinned by complex factors. Emerg Microbes Infect. 2019;8:909–19.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Greenwell B, Boehmke B, Cunningham J, GBM Developers: gbm: generalized boosted regression models. R package version 2022. Accessed 11 Aug 2022.

  46. Sanderson EW, Jaiteh M, Levy MA, Redford KH, Wannebo AV, Woolmer G. The Human Footprint and the Last of the Wild: the human footprint is a global map of human influence on the land surface, which suggests that human beings are stewards of nature, whether we like it or not. Bioscience. 2002;52:891–904.

    Article  Google Scholar 

  47. Warren D, Dinnage R: ENMTools: analysis of niche evolution using niche and distribution models. R package version 1.0.7. 2022. Accessed 18 Nov 2022.

  48. Jin SH, Zhang HL, Yang WH, Zhang YZ, Li JM. Epidemiological analysis of hemorrhagic fever with renal syndrome during 1976 to 2012 in Yunnan province. China Chin J Vector Biol Control. 2013;24:498–502 (in Chinese with English abstract).

    Google Scholar 

  49. Guo T. Epidemic situation of hemorrhagic fever syndrome and investigation of the host animal in some areas in Yunnan Province, 2014-2015. Dali: Dali University; 2016 (in Chinese with English abstract).

  50. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.

  51. QGIS Development Team. QGIS 3.18. Geographic Information System. Open source geospatial foundation project. 2021. Accessed 19 Jan 2023.

  52. Chen ZZ, He K, Cheng F, Khanal L, Jiang XL. Patterns and underlying mechanisms of non-volant small mammal richness along two contrasting mountain slopes in southwestern China. Sci Rep. 2017;7:13277.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Jenkins DW. Trombiculid mites affecting man: II. Control of larval behavior for disease transmission studies. Am J Epidemiol. 1948;48:36–44.

    Article  CAS  Google Scholar 

  54. Jones BM. A method for studying the distribution and bionomics of trombiculid mites (Acarina: Trombidiidae). Parasitology. 1950;40:1–13.

    Article  CAS  PubMed  Google Scholar 

  55. Lin SJ, Guo XG. Study on Leptotrombidium Scutellare and its relation with human disease. J Anhui Agric Sci. 2012;40:188–90 (in Chinese with English abstract).

    Google Scholar 

  56. Xiang R, Guo XG. Research Advances of Leptotrombidium scutellare in China. Korean J Parasitol. 2021;59:1–8.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Rose JB, Epstein PR, Lipp EK, Sherman BH, Bernard SM, Patz JA. Climate variability and change in the United States: potential impacts on water- and foodborne diseases caused by microbiologic agents. Environ Health Perspect. 2001;109:211–20.

    PubMed  PubMed Central  Google Scholar 

  58. Delaune T, Ouattara MS, Ballot R, Sausse C, Felix I, Maupas F, et al. Landscape drivers of pests and pathogens abundance in arable crops. Ecography. 2021;44:1429–42.

    Article  Google Scholar 

  59. Wardrop NA, Kuo CC, Wang HC, Clements AC, Lee PF, Atkinson PM. Bayesian spatial modelling and the significance of agricultural land use to scrub typhus infection in Taiwan. Geospat Health. 2013;8:229–39.

    Article  PubMed  Google Scholar 

  60. Suzuki T. Studies on the bionomics and chemical control of Tsutsugamushi (Scrub-Typhus Mites). Part II. Trombicula scutellaris Nagayo et al. in Southern Kanto of Japan. Jpn J Exp Med. 1954;24:181–97.

    CAS  PubMed  Google Scholar 

  61. Lv Y, Guo XG, Jin DC. Research progress on Leptotrombidium deliense. Korean J Parasitol. 2018;56:313–24.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Sharma A, Dubey VK, Johnson JA, Rawal YK, Sivakumar K. Is there always space at the top? Ensemble modeling reveals climate-driven high-altitude squeeze for the vulnerable snow trout Schizothorax richardsonii in Himalaya. Ecol Indicat. 2021;120:106900.

    Article  Google Scholar 

  63. Wershow ST, Dechaine EG. Retreat to refugia: Severe habitat contraction projected for endemic alpine plants of the Olympic Peninsula. Am J Bot. 2018;105:762–78.

    Article  Google Scholar 

  64. Diez J, Kauserud H, Andrew C, Heegaard E, Krisai-Greilhuber I, Senn-Irlet B, et al. Altitudinal upwards shifts in fungal fruiting in the Alps. Proc Biol Sci. 2020;287:20192348.

    PubMed  PubMed Central  Google Scholar 

  65. Lenoir J, Gegout JC, Marquet PA, de Ruffray P, Brisse H. A significant upward shift in plant species optimum elevation during the 20th century. Science. 2008;320:1768–71.

    Article  CAS  PubMed  Google Scholar 

  66. Hu WH, Onditi KO, Jiang XL, Wu HL, Chen ZZ. Modeling the potential distribution of two species of shrews (Chodsigoa hypsibia and Anourosorex squamipes) under climate change in China. Diversity. 2022;14:87.

    Article  Google Scholar 

  67. Jiang ZK, Wu GH, Wang CJ, Han ZJ, Ding LY, Wang L, et al. General situation of epidemiology of Tsutsugamushi disease of autumn-winter type in China. Chin J Hyg Insectic Equip. 2014;20:209–15 (in Chinese with English abstract).

    Google Scholar 

  68. Xuan Q, Liu AH, Bao FK. Current progrsss in reseach of tsutsugamushi disease in Yunnan Province in China. China Trop Med. 2010;10:1278–80 (in Chinese with English abstract).

    Google Scholar 

  69. Zhang YZ, Yuan JF, Yang XL, Zhou JH, Yang WH, Peng C, et al. A novel hantavirus detected in Yunnan red-backed vole (Eothenomys miletus) in China. J Gen Virol. 2011;92:1454–7.

    Article  CAS  PubMed  Google Scholar 

  70. Li JM, Zhang HL, Yang WH, Zhou JH, Wang JL, Zang YZ, et al. Study on host animals and their genotyping of Hantavirus in Yunnan Province of China in 2007. Endem Dis Bull. 2009;24:1–11 (in Chinese with English abstract).

    Google Scholar 

Download references


We would like to express our sincere thanks to the following people who contributed to the field investigations and laboratory work on chigger mites: Yun-Ji Zou, Qiao-Hua Wang, Rong Fan, Cheng-Fu Zhao, Zhi-Wei Zhang, Ke-Yu Mao, Yong Zhang, Cong-Hua Gao, Nan Zhao, Jian-Chang He, Xue-Song He, De-Cai Ouyang, Chang-Ji Pu, A-Si Di, Cheng-Wei He, some colleagues and college students. We thank Dr. Kenneth Otieno Onditi for the English editing.


This work was supported by the National Natural Science Foundation of China (81960380 and 32260277), the Major Science and Technique Programs in Yunnan Province (202102AA310055-X) and the Open Project from the State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences (GREKF22-03).

Author information

Authors and Affiliations



W-YS and X-GG conceived and designed the study. W-YS, YL, X-GG and P-WY performed the field work. YL and P-WY carried out species identification and organized the data. X-GG confirmed the species identification. W-YS, YL and Y-YY conducted the analyses and wrote the manuscript. All authors contributed to the article. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xian-Guo Guo.

Ethics declarations

Ethics approval and consent to participate

The animal study was reviewed and approved by Dali University, Yunnan, China (approval ID: DLDXLL2020-1104).

Consent for publication

The work is all original research carried out by the authors. All authors agree with the contents of the manuscript and its submission to the journal. No part of the research has been published in any form elsewhere, unless it is fully acknowledged in the manuscript.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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: Figure S1. Marginal effect curves of the 11 environmental variablesover 100 boosted regression treemodels. Red lines indicate the average effect curves and the gray zone marks the predicted 95% confidence interval. Figure S2. Binarizing distribution range of L. scutellare in the near-current period. Figure S3. Binarizing distribution range of L. scutellare in the future periods. a Moderate greenhouse effect in 2050; b high greenhouse effect in 2050, c moderate greenhouse effect in 2090, d high greenhouse effect in 2090. Table S1. Average occurrence probability of L. scutellare in the cities of Yunnan and Sichuan Provinces, and the numbers of case reports of HFRS and Scrub typhus in Yunnan.

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 The Creative Commons Public Domain Dedication waiver ( 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

Song, WY., Lv, Y., Yin, PW. et al. Potential distribution of Leptotrombidium scutellare in Yunnan and Sichuan Provinces, China, and its association with mite-borne disease transmission. Parasites Vectors 16, 164 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Climate change
  • CMIP6
  • Ectoparasitic mite
  • Landscape
  • Prophylactic intervention
  • Species distribution modeling
  • Zoonosis