Skip to main content
Advertisement
  • Loading metrics

Drivers of Echinococcus multilocularis Transmission in China: Small Mammal Diversity, Landscape or Climate?

  • Patrick Giraudoux ,

    patrick.giraudoux@univ-fcomte.fr

    Affiliations Chrono-environment, University of Franche-Comte (CRNS), Besançon, France, Institut Universitaire de France, Paris, France

  • Francis Raoul,

    Affiliation Chrono-environment, University of Franche-Comte (CRNS), Besançon, France

  • David Pleydell,

    Affiliation Institut National de la Recherche Agronomique (INRA), Unité Mixte de Recherche (UMR), Contrle des Maladies Animales Exotiques et Émergentes (CMAEE), Petit-Bourg, Guadeloupe, France

  • Tiaoying Li,

    Affiliation Sichuan Center for Disease Control, Chengdu, China

  • Xiuming Han,

    Affiliation Qinghai Center for Disease Control, Xining, China

  • Jiamin Qiu,

    Affiliation Sichuan Center for Disease Control, Chengdu, China

  • Yan Xie,

    Affiliation Institute of Zoology, Chinese Academy of Sciences, Beijing, China

  • Hu Wang,

    Affiliation Qinghai Center for Disease Control, Xining, China

  • Akira Ito,

    Affiliation Department of Parasitology, Asahikawa Medical University, Asahikawa, Japan

  • Philip S. Craig

    Affiliation School of Environment and Life Sciences, University of Salford, Salford, United Kingdom

Abstract

Background

Human alveolar echinococcocosis (AE) is a highly pathogenic zoonotic disease caused by the larval stage of the cestode E. multilocularis. Its life-cycle includes more than 40 species of small mammal intermediate hosts. Therefore, host biodiversity losses could be expected to alter transmission. Climate may also have possible impacts on E. multilocularis egg survival. We examined the distribution of human AE across two spatial scales, (i) for continental China and (ii) over the eastern edge of the Tibetan plateau. We tested the hypotheses that human disease distribution can be explained by either the biodiversity of small mammal intermediate host species, or by environmental factors such as climate or landscape characteristics.

Methodology/findings

The distributions of 274 small mammal species were mapped to 967 point locations on a grid covering continental China. Land cover, elevation, monthly rainfall and temperature were mapped using remotely sensed imagery and compared to the distribution of human AE disease at continental scale and over the eastern Tibetan plateau. Infection status of 17,589 people screened by abdominal ultrasound in 2002–2008 in 94 villages of Tibetan areas of western Sichuan and Qinghai provinces was analyzed using generalized additive mixed models and related to epidemiological and environmental covariates. We found that human AE was not directly correlated with small mammal reservoir host species richness, but rather was spatially correlated with landscape features and climate which could confirm and predict human disease hotspots over a 200,000 km2 region.

Conclusions/Significance

E. multilocularis transmission and resultant human disease risk was better predicted from landscape features that could support increases of small mammal host species prone to population outbreaks, rather than host species richness. We anticipate that our study may be a starting point for further research wherein landscape management could be used to predict human disease risk and for controlling this zoonotic helminthic.

Author Summary

The loss or gain of certain host species may either dilute or amplify the risk of pathogen and parasite infection through direct or indirect effects. The relative contribution of host communities combined with climate and landscape characteristics on non-vector-borne parasite transmission to humans has been a relatively neglected area of investigation. Here we show that zoonotic transmission of the cestode E. multilocularis in China was not directly correlated with small mammal host species richness, but spatially correlated with alpine meadows, forest characteristics and rainfall that confirmed hotspots of human disease in a 200,000 km2 region of the eastern Tibetan plateau. Our findings indicate more intensive transmission in homogeneous landscapes with larger areas of optimal habitats for one or some host species in low diversity small mammal communities, making multi-annual population outbreaks more likely. Landscape features that could support large population outbreaks of small mammal host species were better predictors of E. multilocularis transmission to humans than indices of host species richness per se. Our results support the notion that landscape, small mammal host biodiversity and their population dynamics may protect humans from zoonotic parasite transmission where they prevent population outbreaks of a few specific small mammal host species.

Introduction

Ecologic systems are nested within one another. This well-known fundamental hierarchical organization [1] is easy to detect in nature but has been generally undervalued as a source of influence on the structure and development of pathogen transmission patterns, and also as a means of understanding the crucial connections between local processes and large-scale distribution patterns. At a community level, Guernier et al. [2] explored the worldwide distribution of human parasitic and infectious diseases (PID) and found that, after correcting for cofactors, PID richness (as for free-living species), was strongly correlated with latitude: PID species diversity decreased as one moved from the equator, and the strongly nested pattern of their global distribution was confirmed. They also pointed out how, along such gradient, the maximum range of precipitation and monthly temperature might be intimately connected in generating the observed pattern of disease diversity. This similarity in the diversity patterns of free living organisms and PIDs suggests that common processes are at work which might be explained at large by the climatically based energy hypotheses (energy availability generates and maintains species richness gradients) [2]. However, for a given PID species distribution patterns at regional scales may be more complex. Transmission may depend not only on species richness but also on host assemblage composition [3]. For instance, reducing host diversity can increase disease transmission when the lost species are either not hosts or suboptimal for the pathogen, especially when population size of optimal hosts are inversely correlated to species richness. On the other hand, a large number of competent host species may provide a much more stable transmission system for transmission, that is robust against environmental disturbances, anthropogenic or natural, that temporarily decrease the density of some host populations. However, empirical examples of the relationships between host biodiversity and parasite transmission are still relatively rare, in part because suitable datasets that may allow comparisons are deficient [4].

Further complexities are encountered when the distribution range of a pathogen covers a large number of host communities and climatic zones. Such a wide-extent distribution-range consists of a nested hierarchy of transmission systems that can be inter-connected in space and time via dispersion [5]. At these scales, datasets with high resolution and precision, describing hosts, disease distribution and environmental factors (such as climate and land cover) are most often heterogeneous.

Human alveolar echinococcocosis (AE) is a highly pathogenic parasitic disease caused by the larval stage of the cestode Echinococcus. multilocularis, which usually results in a slow-growing multivesicualted tumor-like lesion in the liver of cases. The parasite's life-cycle can exploit a large number of small mammal intermediate hosts (>40 species known to date) and several definitive host species (e.g. foxes, coyotes, wolf, dog, etc.). Human infection arises from accidental ingestion of E. multilocularis eggs from direct contact, or via food contaminated by carnivore definitive host faeces. Although patchy, its distribution range covers the Northern hemisphere from the Arctic to the 28th parallel on the Tibetan plateau. In Eurasia, although other carnivores can theoretically sustain the transmission of E. multilocularis, the whole range of the parasite is actually included in the range of its main definitive host, the red fox, Vulpes vulpes, except on the Tibetan plateau where the cestode circulates through the Tibetan fox, Vulpes ferrilata [6]. In contrast, throughout Eurasia, E. multilocularis transmission is sustained by a large variety of different small mammal communities [7]. Within Eurasia, continental China stretches from Siberia in North Xinjiang, Inner Mongolia and Heilongjiang to tropical rain forest in Yunnan, and includes high altitude areas such as the Tibetan plateau, several deserts (some of them below sea level), coastal and agricultural areas (Figure 1). Based on mammal and plant species comparisons a total of 25 biogeographical regions and 77 sub-regions are defined [8]. Transmission of E. multilocularis is sustained in diverse small mammal communities of western and northern China. While some of these communities have been investigated [5], [9] many remain unexplored.

thumbnail
Figure 1. Map of China and main locations cited in the text.

Dotted line, area where regional predictions were computed (see Figure 7), TAR, Tibetan Autonomous Region, LPS, Liupan Shan (Liupan Mountains). Coordinates reference system in degrees (WGS84). The background is Global Land One-kilometer Base Elevation model, provided by the US National Oceanic and Atmospheric Administration.

https://doi.org/10.1371/journal.pntd.0002045.g001

Here we examine the distribution of human AE disease in China on two spatial scales (continental and regional) and test the hypotheses that human disease distribution can be explained by the distribution of intermediate host species richness (i.e. host biodiversity) or by environmental factors such as climate and land cover or an interaction between them. Our central hypothesis is that human AE distribution is best explained by considering the impact of landscape on low species-rich communities prone to population outbreaks, rather than host diversity per se. In that case, for endemic areas, basic information on land cover and climate could be used as a proxy to predict high risk transmission systems for this pathogenic helminthic zoonosis.

Materials and Methods

Continental China

Small mammal data.

Small mammal distributions across China were obtained from a digital version of maps published by Smith et al. [10]. Species nomenclature were standardized according to Wilson and Reeder [11]. Since the main intermediate hosts of E. multilocularis are small mammals, the following taxa were selected: Order Rodentia (family Sciuridae -excluding flying squirrel-, Gliridae, Castoridae, Dipodidae, Platacanthomyidae, Spalacidae, Cricetidae and Muridae), Order Lagomorpha (family Ochotonidae and Leporidae), Erinaceomorpha (family Erinaceidae), 0rder Soricomorpha (family Soricidae and Talpidae), totaling 274 species. A total of 967 point locations were defined on a regular grid (100 km mesh, coordinates reference system: UTM47) covering continental China and a buffer of 100 km radius centered on each node. The ecology and population dynamics of small mammal species is still little known in China. However, some species whose population outbreaks impact agriculture or forestry have been identified as pest species in the literature [5], [10], [12][15]. Evidence of interactions between land cover, human activities and small mammal population dynamics has been provided in Europe and Japan for a number of grassland and forest species [16][18]. Therefore, in the present study, the category ‘pest species’ was considered as a surrogate for those small mammal species whose populations are prone to outbreak and reach large densities locally. Intermediate host species were those previously reported for E. multilocularis [5], [7], [19][22]. Thus the number of species (species richness) belonging to the following categories were computed within each buffer: ‘all species’, ‘not flying small mammals’, ‘species known to be intermediate hosts’, ‘species known to be pests’, ‘species known to be intermediate hosts and pests’.

Environmental and climate data.

Monthly and yearly average rainfall and temperature for China were extracted from the 0.5 degree [23] gridded dataset of the FAO Global maps (http://www.fao.org/geonetwork). Altitude maps were derived from the Global Land One-km Base Elevation (GLOBE) Project [24]. Using the same grid and buffers as above we computed the mean, median, minimum and maximum altitude. Information on land cover was obtained from the Global Land Cover 2000 project [25]. It was based on Spot Vegetation data in the period 1.1 to 31.12.2000 and consisted in 24 land-cover classes mapped at a pixel resolution of 1 km and was designed to cover the full diversity of land-cover found across China, using the FAO Land Cover Classification System.

Human data.

Zhou et al. [26], made an extensive literature search that identified counties (administrative subdivisions of provinces in China) where one or more human AE cases were recorded and derived a distribution map. This map, following Danson et al. [27] was used in the present study.

Eastern Tibetan plateau

Human data.

Mass screening programs for human abdominal echinococcosis were performed from 2002 to 2008 using portable ultrasound in communities that comprised a total of 94 villages, located in Western Sichuan and Qinghai provinces. Persons with confirmative or suspected AE or cystic echinococcosis (CE, caused by E. granulosus) lesions or with other space-occupying lesions in the liver were asked to give a five ml venous blood sample for detection of Echinococcus antibodies using ELISA with recombinant AgB as antigen for CE or AE, and ELISA with recombinant Em18 antigen for AE [28]. In cooperation with County level Centers for Disease Control (CDC) and local health administrators, information about the purpose of the screening program was spread to the villagers. People screened were self-selected volunteers and agreed by informed consent. For each participant a questionnaire was completed with the help of Tibetan registration auxiliaries, to obtain information on gender, age, occupation and ethnicity (see [28] for further details).

Ethics statement.

Approval for the surveys was given by the relevant ethics Committee of the Centers of Disease Control in Sichuan and Qinghai, China. Written consent was obtained from all adult participants and from parents of minors who agreed to participate [28]. Persons with confirmative or suspected AE or CE were assured free diagnosis and medical treatment with long-term albendazole drug therapy for echinococcosis if diagnosis was positive [28]. Recommendation was also provided for possible surgical intervention (cyst/lesion removal or drainage) if appropriate. Persons with other infections or medical conditions were examined and referred to local health clinics for further investigation or treatment.

Data handling and models.

Land cover was estimated as the proportion of pixels belonging to a Global Land Cover 2000 category in question within a circular buffer centered at a given pixel. Elevation, yearly rainfall and yearly average temperature, were extracted from a 0.1 degree gridded dataset of the FAO Global maps, and assigned the median pixel values within the buffer. Such metrics vary with respect not only to landscape composition but also to the buffer radius R. Therefore radii of 10 km, 20 km, 50 km and 100 km were considered. This was repeated for every pixel in a 50×32 grid with a 20×20 km resolution. To avoid resolution differences between human data collected in each village and environmental data available for every pixel of the 20×20 km resolution grid, human data were assigned to the nearest pixel centroid. Generalized additive mixed models (GAMMs) with a logistic link function [29] were used to investigate non-linearities between age, environmental variables and AE status (presence/absence of hepatic AE) of subjects. A within-pixel (or within-village for age) random effect was added to account for pixel or village specific peculiarities in prevalence arising from unobserved pixel or village specific factors. Land cover effects (within buffer percentage of alpine meadows, alpine steppe – termed ‘alpine and sub-alpine plain grass’ in the global land cover 2000 nomenclature –, forest), altitude and rainfall, were modeled similarly for all buffer sizes. The ratio of farmland was nominal within this range of altitude and was not kept for this analysis. The most pertinent value for buffer radius R was identified via Akaike Information Criterion (AIC) based model selection. Selected GAMMs were further subjected to Bayesian analysis using the software BayesX [30], [31] which enabled the inclusion of a spatial random effect to account for otherwise unexplained between-pixel spatial autocorrelation. The spatial random effect was modeled as a Markov random field [32], [33] on the 20 km resolution grid. Non-linearity was modeled using a degree 3 P-splines on 20 equidistant knots and a second order random walk was used to penalize against over fitting, [34], [35]. Non-informative inverse gamma priors were assumed for variance components of the P-spline and random effects with shape and scale hyper-parameters initialized at 0.001). The parameter space of the Bayesian GAMMs was sampled using Markov Chain Monte Carlo (MCMC) techniques. For each model, MCMC was run for 500,000 iterations, discarding an initial burn-in period of 50,000 iterations and sub-sampling every 400th sample thereafter. Model comparison was performed using the deviance information criterion (DIC) [36] with model selection based upon the model returning the lowest DIC. Coefficients and corresponding credibility intervals of selected models were examined. The probability for each linear coefficient to be either strictly positive or strictly negative was computed from their empirical density distribution, considering that the direction of the effect could be expected based on earlier works and preliminary statistics.

Computing and graphical displays were performed using BayesX 2.0.1, R 2.11.1 [37], the R packages BayesX 0.2–4 [38], coda [39], fields 6.3 [40], gstat 0.9–69 [41], mgcv 1.6–2 [29], maptools 0.7–34 [42], rgdal 0.6–2 [43], sp 0.9–65 [44], splancs 2.01–27 [45] and pgirmess 1.4.7 [46].

Results

Continental China

Continental maps of China for distribution of rainfall, altitude and temperature can be found in Supplementary Material. In summary these abiotics factors did not correlate well to human AE disease distribution. Furthermore, the continental distribution of human AE was not clearly correlated to small mammal species richness, even considering the distribution and number of species that are known to be potential intermediate hosts for E. multilocularis or pest species (Figure 2). However, the overlap between the spatial distribution of human AE and the spatial distribution of a combination of two Global Land Cover 2000 categories, i.e. ‘meadows’ and ‘alpine and subalpine meadow’ was very clear. The main endemic areas correspond to the central and eastern Tibetan meadows, those of the Tien Shan and Altai Mountains and northern Inner Mongolia (Figure 3a). On the Tibetan plateau, the human AE endemic foci corresponds to the Global land cover category ‘alpine and subalpine meadow’, characterized by alpine meadows densely covered with thick perennial sedges (Kobresia spp) and various forbs, lying generally below 4500 m as described by Schaller [47]. In the western part of the human AE focus, alpine meadow becomes largely riparian with streams, seepages, swamps and lakes. The area is characteristically composed of a 10–40 cm thick moisture retaining sod layer providing a long growing season. Livestock tend to concentrate on grassy habitat (28–70 animals/km2) and continual grazing and trampling coupled with solifluxion, cause extensive erosion. This anthropogenic activity helps to maintain extensive open habitats favorable to the population surges of the plateau pika (Ochotona curzoniae), a known intermediate hosts of E. multilocularis [48]. By contrast, human AE distribution did not correlate with the Global Land Cover 2000 category ‘alpine and sub-alpine plain grass’ (Figure 3b). This category, lying generally between 4500 and 5000 m, corresponds to an Alpine steppe, dry, cold and windy. Plant coverage is sparse (Stipa spp, Festuca spp, Poa spp, Carex moorcroftii) seldom more than 30%, and soil is poor without a sod layer [47].

thumbnail
Figure 2. Biodiversity of non-flying small mammals.

a, all species; b, pest species; c, intermediate hosts species; d, pest and intermediate host species; red line, limits of human AE areas; circle, Zhang-Puma counties, Gansu; triangle, Xiji-Guyan counties, Ningxia. Vertical bars represent the number of species.

https://doi.org/10.1371/journal.pntd.0002045.g002

thumbnail
Figure 3. Grasslands in China (ratio of total land in a 100 km radius buffer).

a, ‘meadows’ and ‘subalpine meadow’; b, ‘alpine and sub-alpine plain grass’ (Global land cover class terminology) actually corresponding to high altitude steppes; red line, limits of human AE areas; circle, Zhang-Puma counties, Gansu; triangle, Xiji-Guyuan counties, Ningxia.

https://doi.org/10.1371/journal.pntd.0002045.g003

Despite its great size (several 100,000 km2), the Eastern Tibetan plateau focus of human AE disease correlates well with Alpine meadows. The apparent absence of human AE cases on the western and north western Tibetan plateau may simply be due to the fact that human population density decreases from east to west to low population areas of high altitude desert (i.e. Chang Tang area). From east to west on the Tibetan plateau, climatic conditions get drier and colder, which impacts grass growth and yak/goat/sheep survival on which human populations depend. This aridity can also be expected to affect E. multilocularis egg survival. On the other hand, AE cases may possibly just be undetected in small discrete nomadic Tibetan populations of the western Tibetan plateau, since they are typically isolated from public health facilities and no mass screening has occurred to date west of Naqu in Tibet Autonomous Region. We return to this issue in the following section.

Based on this China-wide continental scale analysis, some interesting contrasting results were apparent. Regional areas of south Gansu province (i.e. Zhang-Puma counties), and south Ningxia region (i.e. Xiji-Guyuan counties, close to the North Liu Pan Mountains) were not expected to have large number of human AE cases. However, local foci over an area of 400–2500 km2 were confirmed to have significant AE prevalence in humans reaching an average of 4.1% and 3.0% respectively. These largely montane agricultural zones are considered to be significantly influenced by anthropogenic landscape disturbance caused by deforestation [9], [49], [50]. This process led to the regional opening of forest areas to agriculture, with transitional stages of grassland and pastures triggering small mammal population outbreaks, thus potentially fostering transmission of E. multilocularis locally over a relatively short time-span (10–20 years) when optimal small mammal habitats become temporarily available. Thus, in both the provincial Gansu and Ningxia endemic zones, risk of human AE prevalence increased during the 1970s–90s in response to anthropogenic local landscape changes. However these landscape changes cannot be seen on the coarse resolution land use map.

Eastern Tibetan plateau

A total of 17,589 people were screened for hepatic echinococcosis in 94 villages and the overall prevalence of alveolar echinococcosis was 3.28% (3.02–3.56, 95%CI). Despite substantial sample sizes, no cases were detected to the west of Qinghai Lake (Figure 4). People screened in this area (n = 1,975) were 95.5% Tibetan (71.7% herdsmen, 25.6% children). Geographically isolated from the southern hotspot by the lower Qaidam Basin and the Chaka Yanhu depression, this area was not considered for further analysis (see study area subset Figure 4). Thus, the resulting filtered data set included 15,614 people and 81 villages, with a raw prevalence of 3.7% (3.41–4.01, 95%CI) (Table 1). Univariate analysis detected significantly higher AE prevalence among females than among males (X2 = 11.2, df = 1, p = 0.0008) and among Tibetan than among Han ethnic groups (X2 = 6.05, df = 1, p = 0.01). Evidence of occupation differences was found (X2 = 129, df = 5, p<0.000001), with AE prevalence greater in herdsmen compared to a pool of other categories (X2 = 29.5, df = 1, p<0.000001) and lower in employee (X2 = 5.57, df = 1, p = 0.02), student and children (X2 = 41.9, df = 1, p<0.000001), semi-herdsman (X2 = 7.02, df = 1, p = 0.008) and farmers (X2 = 32.1, df = 1, p<0.000001). Table 2 shows how different occupations were represented between ethnic groups. The link between age categories and some occupations (e.g. ‘student and children’) is obvious. A non-linear relationship of human AE prevalence to age was found (Figure 5).

thumbnail
Figure 4. Raw prevalence and sample size of the 94 villages.

https://doi.org/10.1371/journal.pntd.0002045.g004

thumbnail
Figure 5. Relationships between age, environmental factors and AE status.

Predicted values (solid line) and 95% CI (dotted line) were computed using a general additive mixed model with a binomial link, a cubic regression spline smoothed term on age and a random effect on villages [29]. The rug on the abscissa represents the observations (some at the same position may overlap).

https://doi.org/10.1371/journal.pntd.0002045.g005

thumbnail
Table 1. Prevalence for the current subset excluding Qinghai Lake area (see text).

https://doi.org/10.1371/journal.pntd.0002045.t001

thumbnail
Table 2. Percentage of occupations in each ethnicity category.

https://doi.org/10.1371/journal.pntd.0002045.t002

The buffer radius maximizing the likelihood of the GAMM was found to be 100 km for all land cover variables, altitude, rainfall and temperature. Figure 5 shows the relationships between AE prevalence in humans and environmental variables. AE prevalence was found to increase exponentially with the ratio of Alpine meadows and to decrease with the ratio of forest (corresponding to a linear relationship on the linear predictor – data not shown). Human AE prevalence increased with the percentage of Alpine steppe to a maximum at 15.6% of total land cover, and then decreased to low prevalence. Similar patterns were obtained for rainfall, altitude, and temperature with maximum at 564 mm, 4385 m and −5.62°C respectively. This suggests an increased risk for human AE in areas with a larger percentage of alpine grassland at altitudes ranging between 4200–4600 m. In comparison, landscapes dominated by alpine steppe at higher altitude, mostly observed in the western part of the study area, appeared to be clearly sub-optimal.

We included in models either altitude as a proxy for other ecological parameters (such as rainfall, temperature and other unknown factors) or rainfall and temperature as potential ecological factors impacting directly e.g. Echinococcus egg survival. Models containing non-linear effects on all three variables frequently gave rise to computational difficulties and are not reported here. Models including the spatial random effect consistently returned lower DICs than their non-spatial counterpart (see Table S1 in Supplementary material for the list of models). Evidence for the importance of a non-spatial (pixel) random effect was less overwhelming. DIC comparisons show that 5 models among the 33 fitted were quasi-equivalent with a DIC difference lower than 2. All of them included all variables related to human populations (ethnicity, gender, age, occupation), all landscape variables and at least one variable related to climate or altitude. The mean and 95% credibility intervals of posterior samples for regression coefficients, and the probability that the coefficient be ≤0 (Tibetan, female, herders, meadows) or ≥0 (forest), are shown in Table 3. The proportion of positive coefficient samples were > = 0.98 for Tibetan, > = 0.99 for female and < = 0.03 for forest in all models. 95%CIs for meadow coefficients did not exclude zero or negative values when forest was included, however, >95% of posterior samples for that coefficient were positive. Although the coefficient for herdsmen was greater than zero in less than 90% of samples, dropping this term from models largely increased DIC.

thumbnail
Table 3. Summary statistics of posterior samples for regression coefficients of the 5 best models.

https://doi.org/10.1371/journal.pntd.0002045.t003

For prediction of human AE disease risk, model 5 was selected in order to avoid altitude (a proxy variable whose importance may change with latitude outside of the study area), and in order to maintain the ecologically meaningful variables such as rainfall and temperature which appear not to lead to excessive over-parameterization, even keeping a pixel random effect. Posterior means and standard deviations of the variance parameters for the P-splines and Markov random field were 0.0147 (sd 0.02306) for age, 0.3831 (sd 0.6105) for rainfall, 0.1468 (sd 0.5690) for temperature, 1.2036 (sd 0.5164) for pixel spatial random effect, and 0.0586 (sd 0.0963) for pixel random effect. Figure 6a presents the predicted prevalence for a hypothetical 31.6 year old male, non-Tibetan, non-herdsman (corresponding to the mean age of the sampled population and to zero values for co-factors). The range of the 95% credibility intervals for each pixel of the spatial random effect is presented in Figure 6b. As a result of this analysis, one large hotspot of human AE disease was indicated in an area at the south-east border of Qinghai Province and at the north-west of Sichuan Province. Over a total area of 290,400 km2, an endemic area of 193,600 km2 was calculated to have a human AE prevalence higher than 1 per thousand; for human AE prevalence greater than 1% the predicted transmission zone was 67,200 km2.

thumbnail
Figure 6. Predictions of human AE prevalence.

Bounding box is visible in Figure 1. (a) prevalence, (b) range of the credibility interval from the posterior simulation of the spatial random effect (on the linear predictor scale). Units of horizontal and vertical axes are 1/1000 UTM tile 47, the area corresponds to the area subset of Figure 4. Trend terms were those of model 5. The model succeeds in smoothing the raw prevalence to within a plausible range. The main hotspot is identified in the area at the south-east border of Qinghai and at the north-west of Sichuan. Axes units are in kilometers (UTM tile 47).

https://doi.org/10.1371/journal.pntd.0002045.g006

Predictions were extrapolated across a larger area where rainfall and temperature data remained within the same ranges as the area used to train the model (Figure 7). Over a total area of 902,800 km2, human AE prevalence was predicted to be higher than 1 per thousand over 664,000 km2, and greater than 1% over 210,000 km2. Three other hotspots of human AE disease were predicted from the landscape model. Two were located in the Tibetan Autonomous Region (TAR), one region at the north-west of Naqu and the other to the south of Dingqing. The south Dingqing putative focus was predicted to be connected to the Qinghai-Sichuan hotspot by a crescent shape corridor of elevated prevalence. This appears to be confirmed because a recent initial screening study (n = 232 persons) in Dingqing found a very high AE prevalence of 4.7% (Feng X and Craig PS, unpublished observations). Another hotspot was predicted in Gansu at the north-east limit of the transmission area. Actually, the Gansu hotspot was the first confirmed in the late 1980s [51] (see discussion). Furthermore, the area west of Qinghai Lake was predicted as a low prevalence area, which was confirmed by the present study (this was the subset of 1975 people – 13 villages - where no AE cases could be detected, thus excluded from model fitting).

thumbnail
Figure 7. Prediction of human AE prevalence on the eastern Tibetan plateau.

The vertical gauge gives the predicted prevalence in percent. Empty circles, villages surveyed in the present or earlier studies, triangles,localities cited in the text. Axes units are in kilometers (UTM tile 47).

https://doi.org/10.1371/journal.pntd.0002045.g007

Discussion

Zoonotic pathogens are of growing concern and the role of animal host biodiversity is both of theoretical and practical consideration regarding pathogen stability and the effect of anthropogenic perturbations on the transmission ecosystem [4], [52]. Human alveolar echinococcosis (AE) is a chronic potentially fatal hepatic infection, for which the causative pathogen Echinococcus multilocularis, is transmitted in indirect life-cycles involving mammalian wildlife predator-prey interactions between fox definitive hosts and small mammal intermediate hosts [6], [7]. This cestode zoonosis is known to be highly endemic in parts of Eurasia including China, and provides an excellent model to investigate the role of reservoir host diversity, landscape characters and climate for a non-vector borne zoonotic macro-parasite.

Landscape, human AE distribution and small mammal assemblages

For continental China, we found that the spatial distribution of human AE disease was not clearly correlated to small mammal species richness, even when limiting consideration to the distribution and number of species that are known to be both potential intermediate hosts and/or agricultural/grassland pests. There was however a clear spatial overlap between the distribution of certain grassland types, and the spatial distribution of human AE. The main endemic areas corresponded to the central and eastern Tibetan meadows, those of the Tien Shan and Altai Mountains and of northern Inner Mongolia. Rainfall, altitude and temperature did not directly correlate with human AE disease distribution at continental scale.

Human AE prevalence was found to exponentially increase with the ratio of alpine meadows and decrease with the ratio of forests in our regional analysis. The resolution of small mammal atlases in China did not facilitate inclusion of biodiversity indices in regional models. However, the Qinghai-Sichuan human AE disease hotspot lies in one of the areas with the lowest regional biodiversity of small mammals in China (Figure 1). A comprehensive small mammal survey carried out in Shiqu county (see location Figure 6), in the middle of the Qinghai-Sichuan human AE hotspot, recorded 6 species only of which Microtus spp and Ochotona curzoniae were classified as pests [48]. By contrast, 15 small mammal species were found in the Rangtang/Maerkang locality of west Sichuan, in the more forested part of the study area, among which a maximum of 10 species were present in some forest habitats [53]. Some of those species potentially outbreaking were shared with the east Tibetan Shiqu area (i.e. Ochotona cansus, Microtus irene, M. limnophilus), but no indications of small mammal population surges were found from the survey and from local farmer interviews. Moreover human AE prevalence in these areas was lower (1.5% versus 6% in Shiqu county) [28].

Impact of landscape on small mammal population dynamics and fox diet

Evidence of interactions between landscape and arvicolid vole population dynamics at various spatial scales has been provided in earlier studies in Western Europe. On the regional scale (area of about 2500 km2) larger variations in vole population densities occur where permanent grassland cover exceeded 50% and 85% of the total land for M. arvalis and A. terrestris respectively [16], [17]. In complex ecosystems, the population dynamics of small mammals is regulated by both top-down (predation, parasitism) and bottom-up (resources) forces in a multivariate context [54][57]. In general, this means that population surges are less likely in biodiverse communities of small mammals, in regions characterized by heterogeneous landscapes (e.g. forest mosaics). Towards the far west and north of the Tibetan plateau, the increasing aridity of Alpine steppes and semi-deserts and the more patchy distribution of Alpine meadows decrease both primary production and connectivity of optimal habitats, thus the probability of large scale small mammal outbreaks. By contrast extensive grasslands with higher grass production (Alpine ‘meadows’) of the Eastern Tibetan plateau forms optimal conditions on a regional scale for surges of potentially cyclic small mammal species such as the plateau pika (Ochotona curzoniae), several species of Microtus (M. leucurus, M. irene, M. limnophilus), and probably hamsters (Cricetulus sp.). In such contexts, definitive hosts, foxes, take advantage of the most abundant and accessible resources and specialize on them [58][60]. Raoul et al. [61] have shown that infection of foxes with E. multilocularis responds quickly and asymptotically to small increases in the densities of favored prey species. This may cause intense environmental contamination and human exposure, directly from fox feces or more likely via dogs infected from preying upon abundant (infectious) small mammal intermediate host reservoirs [62]. Such system processes may explain why, on the Tibetan plateau, regional landscape variables were found to help to predict human AE distribution.

Co-linearities on the Tibetan plateau

The greatest difficulty in interpreting these results arises from the fact that in such regional systems, explanatory variables which may determine transmission, directly or not, are not independent. For instance, grass productivity and forest cover are correlated with rainfall and temperature, which are correlated (locally) with altitude. This may explain why apparent discrepancies were observed between univariate and multivariate statistics. For instance, most Tibetans were also livestock herders, which may have nullified the occupation variable in the multivariate analysis. Furthermore, human AE disease was found to be correlated with Alpine meadow cover in univariate statistics, but in multivariate statistics effect size of alpine meadow was smaller when the forest term was included in models.. However perfect colinearity between forest and meadow aerial covers was ruled out by examining covariation matrices (data not shown) and because discarding one of the two variables in multivariate models led to increased DICs. The spatial random effect may also have nullified other effects suggesting spatial heterogeneity in effect size or the role of other unmeasured factors (e.g. socio-economics, etc.).

Peripheral foci, landscape and infection dynamics

Large unexplained differences between neighboring villages were observed on the Tibetan plateau. Such differences are commonly reported in most other studies in endemic AE areas, for instance in Gansu [27], [49], Ningxia [50] and France [9]. Our study shows that the Zhang county local AE focus (about 400 km2) in south Gansu Province lies in an area where regionally, rainfall and temperature were in the range of those found in the Qinghai-Plateau hotspot, although at a much lower altitude (2000–2800 m). The Zhang focus, first discovered in the late 1980s [51] was comprehensively described in studies carried out in the late 90s [9], [49], [63]. An overall human AE prevalence of 4.1% was reported there with some villages reaching more than 10%. Some authors [9], [27], [64] had previously provided indications that E. multilocularis transmission in Zhang was the consequence of a transient augmentation in shrub and grass cover in the 1980s, generated by successional growth following deforestation and triggering population outbreaks of the vole M. limnophilus and the hamster C. longicaudatus. Following 20–30 years of deforestation the parasite life-cycle is no longer being maintained in the current farmland landscape due to lack of both suitable intermediate and definitive hosts [5]. This suggest that in areas where climatic conditions are expected to be favorable, but landscape unfavorable, fine grain landscape alteration (here of anthropogenic origin) may create transitory local foci of E. multilocularis transmission that are prone to extinction. How the parasite may colonize such areas from larger and more stable regional foci remains open to question (mainland/island dynamics through fox movements and/or dog trade?) [5]. The present study however shows that even those risk areas may be predicted from climate and landscape analysis at medium (20 km×20 km) resolution.

In conclusion, our results indicate that the prevalence of the cestdode zoonosis, human alveolar echinococcosis, was higher in areas of relatively low small mammal biodiversity than in more diverse host communities within which they were nested. In low diversity small mammal communities, potential intermediate host species prone to population outbreaks (e.g. rodents and lagomorphs) can reach higher densities and enhance transmission. This can be a natural occurrence, as on the homogeneous meadows of the eastern Tibetan plateau where rainfall maintains relatively higher productivity than is possible on the western half of the Tibetan plateau. Conversely, it can also be the result of anthropogenic driven landscape alteration, as in south Gansu where deforestation under favorable climatic conditions led to a transitory decrease of small mammal local biodiversity, and an augmentation of habitats and landscapes favorable to one or several cyclic species of those low diversity communities. Our results support the notion that landscape, small mammal host biodiversity and their population dynamics may protect humans from E. multilocularis transmission by preventing population outbreaks of specific small mammal host species with subsequent consequences on prey/predator relationships and associated host/parasite transmission [65]. This mechanism is clearly different to the dilution and zooprophylaxis effects that have been described and debated elsewhere as potential utilitarian advantages of biodiversity conservation [4]

Supporting Information

Figure S1.

Altitude statistics in continental China (in a 100 km radius buffer) and human alveolar echinococcosis distribution (red lines).

https://doi.org/10.1371/journal.pntd.0002045.s001

(TIF)

Figure S2.

Rainfall statistics (logarithm) in continental China (in a 100 km radius buffer) and human alveolar echinococcosis distribution (red lines).

https://doi.org/10.1371/journal.pntd.0002045.s002

(TIF)

Figure S3.

Average temperature in continental China (in a 100 km radius buffer) and human alveolar echinococcosis distribution (red lines).

https://doi.org/10.1371/journal.pntd.0002045.s003

(TIF)

Figure S4.

Land cover in continental China (ratio of total land in a 100 km radius buffer) according to GLOBAL land cover 2000 nomenclature and human alveolar echinococcosis distribution (red lines). V1, needleleaved deciduoud forest; V2, needleleaved evergreen forest; V3, broadleaved evergreen forest; V4, broadleaved deciduoud forest; V5, bush; V6, sparse woods; V7, seaside wet lands; V8, alpine and subalpine meadow; V9, slope grassland; V10, plain grassland; V11, desert grassland; V12, meadow; V13, city; V14, river; V15, lake; V16, swamp; V17, glacier; V18, bare rocks; V19, gravels; V20, desert; V21, farmland; V22, alpine and sub-alpine plain grass; V23, Mosaic of cropping; V24, Forest Mosaic/Degraded Forest;

https://doi.org/10.1371/journal.pntd.0002045.s004

(TIF)

Table S1.

Model comparison with deviance information criterion. Spatial and random effects are depicted by spatial and pixel respectively. Tibetan, Female, Herdpeople are factors coded 0/1; ratio meadows and ratio forests are the ratio of Alpine meadows and of forest to total land within a 100 km radius buffer. P-spline term was added on age, altitude, rainfall and temperature.

https://doi.org/10.1371/journal.pntd.0002045.s005

(DOC)

Acknowledgments

The authors thank Dr. Chen Xingwang, Professors Qiu Dongchuan, Wang Qian, Yang Wen and Xiao Ning (Institute of Parasitic Diseases, Sichuan Centers for Disease Control and Prevention) and Professor Mark Danson (University of Salford) for their participation and dedication throughout our projects in Sichuan. We also thank countless local people who have regularly provided invaluable support and aid during our field expeditions.

Author Contributions

Conceived and designed the experiments: PG PSC FR TL AI. Performed the experiments: PG TL XH JQ YX HW. Analyzed the data: PG FR DP. Wrote the paper: PG FR PSC DP.

References

  1. 1. Allen TE, Star TB (1982) Hierarchy perspectives for ecological complexity. Chicago: University of Chicago Press.
  2. 2. Guernier V, Hochberg ME, Guegan JFO (2004) Ecology drives the worldwide distribution of human diseases. PLoS Biol 2: 740–746.
  3. 3. Keesing F, Belden L, Daszak P, Dobson A, Drew Harvell C, et al. (2010) Impacts of biodiversity on the emergence and transmission of infectious diseases. Nature 468: 647–652.
  4. 4. Randolph SE, Dobson AD (2012) Pangloss revisited: a critique of the dilution effect and the biodiversity-buffers-disease paradigm. Parasitology 139: 847–863.
  5. 5. Giraudoux P, Pleydell D, Raoul F, Quere JP, Wang Q, et al. (2006) Transmission ecology of Echinococcus multilocularis: What are the ranges of parasite stability among various host communities in China? Parasitol Int 55: S237–S246.
  6. 6. Wang Z, Wang X, Liu X (2008) Echinococcosis in China, a Review of the Epidemiology of Echinococcus spp. Ecohealth 5: 115–126.
  7. 7. Rausch RL (1995) Life-cycle patterns and geographic distribution of Echinococcus species. In: Thomson RCA, Lymbery AJ, editors. Echinococcus and hydatid disease. Wallingford - Oxon: CAB International. pp. 89–134.
  8. 8. Xie Y, Mackinnon J, Li D (2004) Study on biogeographical divisions of China. Biodivers Conserv 13: 1391–1417.
  9. 9. Giraudoux P, Craig PS, Delattre P, Bartholomot B, Bao G, et al. (2003) Interactions between landscape changes and host communities can regulate Echinococcus multilocularis transmission. Parasitology 127: 121–131.
  10. 10. Smith A, Xie Y (2008) A Guide to the Mammals of China. Princeton: Princeton University Press. 576 p.
  11. 11. Wilson DE, Reeder DM (2005) Mammal species of the world. Baltimore: Johns Hopkins University Press. 1207 p.
  12. 12. Deng Z, Wang CX (1984). Rodent control in China. In: Proceedings of the 11th Vertebrate Pest Conference. University of Nebraska - Lincoln, USA. Available: http://digitalcommons.unl.edu/cgi/viewcontent.cgi?article=1043&context=vpc11
  13. 13. Kang L, Han XG, Zhang ZB, Sun OJ (2007) Grassland ecosystems in China: review of current knowledge and research advancement. Philos Trans R Soc Lond B Biol Sci 362: 997–1008.
  14. 14. Stenseth NC, Leirs H, Skonhoft A, Davis SA, Pech RP, et al. (2003) Mice, rats, and people: the bio-economics of agricultural rodent pests. Front Ecol Environ 1: 367–375.
  15. 15. Wang ZH (1984). Acute rodenticides in the control of rodent pest in China: a review. In: Proceedings of the 11th Vertebrate Pest Conference. University of Nebraska - Lincoln, USA. Available: http://digitalcommons.unl.edu/cgi/viewcontent.cgi?article=1042&context=vpc11
  16. 16. Delattre P, Giraudoux P, Baudry J, Truchetet D, Musard P, et al. (1992) Land use patterns and types of common vole (Microtus arvalis) population kinetics. Agric Ecosyst Environ 39: 153–169.
  17. 17. Giraudoux P, Delattre P, Habert M, Quere JP, Deblay S, et al. (1997) Population dynamics of fossorial water vole (Arvicola terrestris scherman): a land usage and landscape perpective. Agric Ecosyst Environ 66: 47–60.
  18. 18. Saitoh T, Stenseth NC, Bjornstad ON (1998) The population dynamics of the vole Clethrionomys rufocanus in Hokkaido, Japan. Res Popul Ecol 40: 61–76.
  19. 19. Craig PS (2006) Epidemiology of human alveolar echinococcosis in China. Parasitol Int 55: S221–S225.
  20. 20. Eckert J, Gemmell MA, Meslin F-X, Pawłowski ZS, editors (2001) WHO/OIE Manual on Echinococcosis in Humans and Animals: a Public Health Problem of Global Concern. Paris: World Organisation for Animal Health and World Health Organization,. 286 p.
  21. 21. Vuitton DA, Zhou H, Bresson-Hadni S, Wang Q, Piarroux M, et al. (2003) Epidemiology of alveolar echinococcosis with particular reference to China and Europe. Parasitology 127: S87–S107.
  22. 22. Wang ZH, Wang XM, Liu XQ (2008) Echinococcosis in China, a review of the epidemiology of Echinococcus spp. Ecohealth 5: 115–126.
  23. 23. New M, Hulme M, Jones PD (1999) Representing twentieth century space-time climate variability. Part 1: development of a 1961–90 mean monthly terrestrial climatology J Clim 12: 829–856.
  24. 24. National Geophysical Data Center (1999) Global Land One-km Base Elevation (GLOBE) Project. Version 1.0 ed. Boulder, CO: National Geophysical Data Center, NESDIS, NOAA, U.S. Department of Commerce.
  25. 25. Wu B, Xu W, Huang H, Yan C (2003) The Land Cover Map for China in the Year 2000. GLC2000 database: European Commision Joint Research Centre.
  26. 26. Zhou HX (2001) Epidemiologie et écologie de la transmission d'Echinococcus multilocularis au Xinjang, Republique Populaire de Chine: étude preliminaire. Besançon: University of Franche-Comté.
  27. 27. Danson FM, Graham AJ, Pleydell DRJ, Campos-Ponce M, Giraudoux P, et al. (2003) Multi-scale spatial analysis of human alveolar echinococcosis risk in China. Parasitology 127: 133–141.
  28. 28. Li TY, Chen XW, Zhen R, Qiu JM, Qiu DC, et al. (2010) Widespread co-endemicity of human cystic and alveolar echinococcosis on the eastern Tibetan Plateau, northwest Sichuan/southeast Qinghai, China. Acta Trop 113: 248–256.
  29. 29. Wood SN (2006) Generalized additive models: an introduction with R. London: Chapman & Hall/CRC. 392 p.
  30. 30. Brezger A, Kneib T, Lang S (2005) Bayesx: Analyzing Bayesian structured additive regression models. J Stat Softw 14: 1–22.
  31. 31. Belitz C, Brezger A, Kneib T, Lang S (2009) BayesX - Software for Bayesian inference in structured additive regression models. version 2.0.1 (21.10.2009) ed.
  32. 32. Besag J, York J, Mollie A (1991) Bayesian image restoration with two applications in spatial statistics (with discussion). Ann I Stat Math 43 1: 59.
  33. 33. Rue H (2001) Fast sampling of gaussian markov random fields with applications. J R Stat Soc Series B Stat Methodol 63: 325–338.
  34. 34. Lang S, Brezger A (2004) Bayesian p-splines. J Comput Graph Stat 13: 183–212.
  35. 35. Brezger A, Lang S (2006) Generalized additive regression based on Bayesian P-splines. Comput Stat Data Anal 50: 967–991.
  36. 36. Spiegelhalter DJ, Best NG, Carlin BP, van der Linde A (2002) Bayesian measures of model complexity and fit. J R Stat Soc Series B Stat Methodol 65: 583–639.
  37. 37. R Development Core Team (2010) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
  38. 38. Kneib T, Heinzl F, Brezger A, Sabanes Bove D (2009) BayesX: R Utilities Accompanying the Software Package BayesX.. version 0.2-4 ed: CRAN.
  39. 39. Plummer M, Best N, Cowles K, Karen V (2010) coda: Output analysis and diagnostics for MCMC. version 0.13-5 ed: CRAN.
  40. 40. Furrer R, Nychka D, Sain S (2010) fields: Tools for spatial data. version 6.3. ed: CRAN.
  41. 41. Pebesma EJ (2004) Multivariable geostatistics in S: the gstat package. Comput Geosci 30: 683–691.
  42. 42. Lewin-Koh N, Bivand RS (2012) maptools: Tools for reading and handling spatial objects. R package version 0.8-16.: CRAN.
  43. 43. Keitt T, Bivand R, Pebesma E, Rowlingson B (2012) rgdal: Bindings for the Geospatial Data Abstraction Library. R package version 0.7-12.: CRAN.
  44. 44. Pebesma E, Bivand R (2005) Classes and methods for spatial data in R.. R News 5: 9–13.
  45. 45. Bivand RS, Pebesma EJ, Gomes-Rubio V (2008) Applied spatial data analysis with R. Dortrecht: Springer. 374 p.
  46. 46. Giraudoux P (2010) pgirmess: miscellaneous functions for data analysis in ecology. version 1.4.7 ed: CRAN.
  47. 47. Schaller GB (1998) Wildlife of Tibetan Steppe. Chicago: The University of Chicago Press. 373 p.
  48. 48. Raoul F, Quere J-P, Rieffel D, Bernard N, Takahashi K, et al. (2006) Distribution of small mammals along a grazing gradient on the Tibetan plateau of western Sichuan, China. Mammalia 42: 214–225.
  49. 49. Craig PS, Giraudoux P, Shi D, Bartholomot B, Barnish G, et al. (2000) An epidemiological and ecological study of human alveolar echinococcosis transmission in south Gansu, China. Acta Trop 77: 167–177.
  50. 50. Pleydell D, Yang Y, Wang Q, Raoul F, Danson M, et al. (2008) Landscape composition and spatial prediction of a zoonosis: the case of alveolar echinococcosis in southern Ningxia, China. PLoS Negl Trop Dis 2: e287 .
  51. 51. Craig PS, Deshan L, Macpherson CNL, Dazhong S, Reynolds D, et al. (1992) A large focus of alveolar echinococcosis in central China. Lancet 340: 826–831.
  52. 52. Keesing F, Holt RD, Ostfeld RS (2006) Effects of species diversity on disease risk. Ecol Lett 9: 485–498.
  53. 53. Vaniscotte A, Pleydell DRJ, Raoul F, Quere JP, Qiu JM, et al. (2009) Modelling and spatial discrimination of small mammal assemblages: An example from western Sichuan (China). Ecol Modell 220: 1218–1231.
  54. 54. Lidicker WZ (2000) A food web/landscape interaction model for microtine rodent density cycles. Oikos 91: 435–445.
  55. 55. Lindström J, Ranta E, Kokko H, Lundberg P, Kaitala V (2001) From arctic lemmings to adaptative dynamics: Charles Elton's legacy in population ecology. Biol Rev Camb Philos Soc 76: 129–158.
  56. 56. Korpimäki E, Brown P, Jacob J, Pech R (2004) The puzzles of population cycles and outbreaks of small mammals solved? Bioscience 54: 1071–1079.
  57. 57. Morilhat C, Bernard N, Foltête J, Giraudoux P (2008) Neighbourhood landscape effect on population kinetics of the fossorial water vole (Arvicola terrestris scherman). Landsc Ecol 23: 569–579.
  58. 58. Ferrari N, Weber JM (1995) Influence of the abundance of food resources on the feeding habits of the red fox, Vulpes vulpes, in western Switzerland. J Zool 236: 117–129.
  59. 59. Sidorovich V, Sidorovich A, Izotova I (2006) Variations in the diet and population density of the red fox Vulpes vulpes in the mixed woodlands of northern Belarus. Mamm Biol 71: 74–89.
  60. 60. Dell'Arte G, Laaksonen T, Norrdahl K, Korpimak iE (2007) Variation in the diet composition of a generalist predator, the red fox, in relation to season and density of main prey. Acta Oecol 31: 276–281.
  61. 61. Raoul F, Deplazes P, Rieffel D, Lambert JC, Giraudoux P (2010) Predator dietary response to prey density variation and consequences for cestode transmission. Oecologia 164: 129–139.
  62. 62. Vaniscotte A, Raoul F, Poulle ML, Romig T, Dinkel A, et al. (2011) Role of dog behaviour and environmental fecal contamination in transmission of Echinococcus multilocularis in Tibetan communities. Parasitology 138: 1316–1329.
  63. 63. Giraudoux P, Quere JP, Delattre P, Bao GS, Wang XM, et al. (1998) Distribution of small mammals along a deforestation gradient in southern Gansu, central China. Acta Theriol 43: 349–362.
  64. 64. Graham A, Danson F, Giraudoux P, Craig P (2004) Ecological epidemiology: landscape metrics and human alveolar echinococossis. Acta Trop 91: 267–278.
  65. 65. Patz JA, Daszak P, Tabor GM, Aguirre AA, Pearl M, et al. (2004) Unhealthy landscapes: Policy recommendations on land use change and infectious disease emergence. Environ Health Perspect 112: 1092–1098.