Skip to main content
Download PDF

Patterns of helminth infection in Kenyan elephant populations

Parasites & Vectors volume 13, Article number: 145 (2020) Cite this article

Abstract

Background

The dynamics of helminth infection in African elephant populations are poorly known. We examined the effects of age, sex, social structure and the normalized difference vegetation index (NDVI) as primary drivers of infection patterns within and between elephant populations.

Methods

Coprological methods were used to identify helminths and determine infection patterns in distinct elephant populations in Maasai Mara National Reserve, Tsavo East National Park, Amboseli National Park and Laikipia-Samburu Ecosystem. Gaussian finite mixture cluster analyses of egg dimensions were used to classify helminth eggs according to genera. Generalized linear models (GLM) and Chi-square analyses were used to test for variation in helminth infection patterns and to identify drivers in elephant populations.

Results

Helminth prevalence varied significantly between the studied populations. Nematode prevalence (96.3%) was over twice as high as that of trematodes (39.1%) in elephants. Trematode prevalence but not nematode prevalence varied between populations. Although we found no associations between helminth infection and elephant social groups (male vs family groups), the median helminth egg output (eggs per gram, epg) did vary between social groups: family groups had significantly higher median epg than solitary males or males in bachelor groups. Young males in mixed sex family groups had lower epg than females when controlling for population and age; these differences, however, were not statistically significant. The average NDVI over a three-month period varied between study locations. Cluster analyses based on egg measurements revealed the presence of Protofasciola sp., Brumptia sp., Murshidia sp., Quilonia sp. and Mammomonogamus sp. GLM analyses showed that the mean epg was positively influenced by a three-month cumulative mean NDVI and by social group; female social groups had higher epg than male groups. GLM analyses also revealed that epg varied between elephant populations: Samburu-Laikipia elephants had a higher and Tsavo elephants a lower epg than Amboseli elephants.

Conclusions

Elephants had infection patterns characterized by within- and between-population variation in prevalence and worm burden. Sociality and NDVI were the major drivers of epg but not of helminth prevalence. Gastrointestinal parasites can have a negative impact on the health of wild elephants, especially during resource scarcity. Thus, our results will be important when deciding intervention strategies.

Background

Most studies on the helminths parasitizing African elephants have in the past focused on helminth taxonomy and more recently on within population infection dynamics [1,2,3,4], but no studies have simultaneously examined inter-population and intra-population infection dynamics and their drivers. The most common helminths infecting African elephants are nematodes and trematodes; two groups of helminths that have environmentally dependent transmission mechanisms. Species of Murshidia, Quilonia and Khalilia are the most common nematodes infecting African elephants [5]. The free-living environmental stages of gastrointestinal nematodes are strongly affected by climate, e.g. extreme temperatures are detrimental to their development and survival. Moisture is needed for the development and transition of larvae from soil to pasture and so rainfall and vegetation may be limiting factors on transmission and may influence patterns of inter-population variability in infection patterns [6]. Some studies in Africa have found a significant positive correlation between mean annual precipitation (rainfall and relative humidity) and nematode infection rates [6,7,8]. Furthermore, some studies have also found associations between precipitation and certain qualitative measurements of egg burden (mean nematode species richness, mean number of nematode worms and infection intensity per individual host) [9,10,11].

Several species of trematodes are known to infect elephants [12] and some are associated with pathological lesions in starving animals [12, 13]. Trematodes of the family Fasciolidae usually have a complex life-cycle that involves a vertebrate host, in which they reproduce sexually, and an intermediate snail host, in which asexual reproduction occurs. The transmission of trematodes is largely driven by the presence of water and snails, suggesting that water availability and precipitation are important factors in their life histories [14]. Given the association between climatic factors and the propagation and transmission of both nematodes and trematodes, we used a normalized difference vegetation index (NDVI) since it is generally strongly correlated with climatic parameters (precipitation and temperature) and soil moisture content. These factors directly or indirectly influence host-parasite relationships and the propagation of environmentally transmitted helminths and may be important driver of between population variation in infection patterns.

The African elephant utilizes a wide range of habitats and lives in socially structured contact networks. Within elephant populations, individual animals live in structured groups that exhibit fission-fusion dynamics varying between sex-age groups [15,16,17,18]. Females and their offspring form fusion-fission matriarchal social groups where adult females and their calves live in stable units that coalesce with other similar cow-calf groups to create family and bond groups, thereby allowing adult females to form a nested or hierarchical social structure [15]. Males, on the other hand, form fluid social groups of mixed or similarly aged males in bachelor groups that have periodical contact with matriarchal groups when searching for mating opportunities [16]. Due to these more fluid social dynamics and greater mobility, males rove more widely than females [19] and adult males have larger home ranges than the immature males that still form part of the family groups [20,21,22].

The aims of this study were to examine helminth infection patterns between and within the most important elephant populations in Kenya found in a number of different agro-ecological zones, and to test the importance of influence of age, social structure and NDVI as drivers of these infection patterns.

Methods

Study area

The study was carried out in Tsavo East National Park, Laikipia-Samburu Ecosystem, Maasai Mara National Reserve and Amboseli National Park (Fig. 1), the four conservation areas that hold the largest elephant populations in Kenya. These populations remain separate and do not mix.

Fig. 1
figure 1

Map of Kenya showing the locations of the four major elephant populations in Kenya

Full size image

Tsavo East National Park (TENP) is situated in south-east of Kenya and enjoys a semi-arid savannah climate with a bimodal annual rainfall pattern. Heavy rains occur in April-May, while light rains fall in November-December. Overall, rainfall is erratic and low, with an annual average of 300–600 mm. This area holds 7727 elephants according to the 2018 large mammal census conducted by the Kenya Wildlife Service. Laikipia-Samburu Ecosystem (LSE) is in central Kenya and is covered by arid savannah grassland with annual rainfall of 300–700 mm. Rainfall is bimodal and falls in April-May and November-December. The Laikipia-Samburu Ecosystem hosts an elephant population estimated at 7166 during the 2017 large mammal census carried out by Kenya Wildlife Service. The Maasai Mara National Reserve (MMNR) is located in southern Kenya close to the border between Kenya and Tanzania, where it is contiguous with the Serengeti. Overall, this area consists of a large expanse of savannah grassland with annual rainfall ranging from 650 mm in the south-east to about 1300 mm in the north-west. Most rain falls in March-May, although some also falls in October-November. The 2017 large mammal census reported 2493 elephants in this national reserve. Lastly, Amboseli National Park (ANP) is located at the base of Mount Kilimanjaro in southern Kenya. It holds 2127 elephants according to the 2018 large mammal census report by the Kenya Wildlife Service. It consists mostly of arid dry savannah open grassland land, mixed with patches of scrub and Acacia xanthophloea woodland. Average annual rainfall is 340 mm with an annual range of 141–757 mm (https://amboselibaboons.nd.edu/downloads/). A network of marshes fed by underground water originating as snow melt from Mount Kilimanjaro provides a permanent water supply.

Faecal sampling

Faecal sampling was carried out in February-November 2017 using a cross-sectional study design. Individuals in a social family herd, male bachelor herds and lone bulls were tracked until they defecated. From each animal defecation, a fresh dung bolus was carefully opened and approximately 20 g of the dung were scooped out and preserved in 10% formalin. The following information was recorded for each sample: age of the individual animal (adult, subadult or juvenile), sex, date of collection, GPS coordinates at the time of sighting, and type of social group (whether part of a female or male social group). A female social group was defined as a group consisting of females and their offspring and occasional males, whereas a male social group was taken to be as a solitary male or a group of two or more males seen in proximity at the time of observation. A total of 243 faecal samples, 71 from independent male groups or solitary males and 172 from family social groups, were collected in the four study areas. Totals of 62 family groups were sampled in MMNR, 37 in TENP, 27 in ANP and 19 in LSE, while 19 male social groups were sampled in MMNR, 22 in TENP, 16 in ANP and 14 in LSE.

Coprological analyses

Sedimentation technique

A method described by VanderWaal et al. [23] with a slightly modified procedure was applied. Approximately 4 g of dung was weighed, mixed with 45 ml of tap water in a 50 ml centrifuge tube, and stirred until the mixture became a slurry. The dung slurry was then sieved and left to stand for 30 min. Decanting and re-suspension of the sediment was repeated 2–3 times until the suspension cleared. A dropping pipette was used to place ~ 0.05 ml of the sediment on a glass slide for examination under a Leica DM500 microscope. The presence of nematode and trematode eggs was assessed and micrographs of at least 10 eggs of different morphotypes were taken. Measurements of eggs (length and width) were taken from the photomicrographs using Leica LAS EZ software (Leica Microsystems GmbH, Wetzlar, Germany).

Flotation technique

Faecal samples were thoroughly homogenized with a stirring stick so that parasite eggs would be uniformly distributed throughout the sample. Initially, a faecal floatation fluid with specific gravity of 1.27 was prepared. Briefly, 454 g of table sugar was weighed and mixed with 355 ml of distilled water. The mixture was heated over low heat whilst being stirred until all the sugar had dissolved. The sugar solution was left to cool before use as the floatation fluid. Faecal samples were homogenized (as in the sedimentation technique) and prepared by weighing approximately 4 g of the elephant dung. The sample was mixed with 12 ml tap water, stirred and sieved through a tea strainer, before being transferred to a 15 ml plastic centrifuge tube. If the filtrate was less than 15 ml, it was topped up with tap water and the tube capped. It was then centrifuged at 1500× rpm for 10 min. The supernatant was decanted out and the sediment re-suspended using the flotation fluid to fill up half the test tube. The sediment was mixed thoroughly with the flotation fluid using a stirring stick. The tube was then filled to the top with more flotation fluid until a slight bulging meniscus formed. A coverslip was gently placed on the centre of the top of the tube. The tubes were then centrifuged for 10 min at 1500× rpm. After centrifugation, the coverslip was gently removed and placed directly onto a clean glass slide for examination under the microscope. Helminth eggs were qualitatively assessed. Photomicrographs of at least 10 eggs of different morphotypes were taken and processed as described in the sedimentation section.

McMaster technique

The helminth eggs were counted using a quantitative technique based on a calibrated McMaster chamber. Egg counts give an estimate of the number of eggs per gram (epg) in the faecal sample. The faecal sample was prepared as described for the floatation technique. A pipette was used to transfer the mixture to each of the two chambers of the McMaster slide. The preparation on the slide was left to settle for at least 5 min and then examined under the microscope. The eggs present in each chamber were counted. The total count for the slide was multiplied by a constant (50) to give the number of eggs per gram.

Normalized difference vegetation index (NDVI) analysis

NDVI is a measure of reflectance and absorbance of the light spectra by vegetation and depends on the phenology and density of the vegetation being tested. Green vegetation reflects mostly green and near-infrared light spectra and absorbs red and blue light spectra and is often used as an index of productivity as it is correlated to plant phenology and nitrogen content. Overall, NDVI is strongly influenced by climatic parameters such as precipitation and temperature, and by soil moisture. These factors directly or indirectly influence host-parasite relationships and the transmission of environmentally transmitted helminths.

NDVI data were taken from satellite images obtained from a Landsat 8 satellite using the Operational Land Imager and Thermal Infrared Sensors for data capture. Satellite images of 30 m resolution were retrieved from the Libra development seed website (https://developmentseed.org/projects/libra/). Shape files of the four study areas were obtained from the Kenya Wildlife Service and were used to resize the relevant satellite images. The Amboseli satellite images were downloaded from 168 path and 062 row in April-November 2017; images for Maasai Mara were downloaded from 169 path and 061 row in December 2016 and January-March 2017; images for Laikipia-Samburu were downloaded from 168 path and 060 row in April-August 2017; and images for Tsavo East were downloaded from 167 path 062 row and from 163 path and 062 row in October-December 2016 and January-March 2017. These periods coincided with the sampling months and the three months prior to sampling. The Tsavo East satellites images were mosaicked into a single image using Q GIS software (Creative Commons, Mountain View, USA).

Selected satellite images were pre-processed using Q GIS to remove both radiometric and geometrical errors. The corrected images were used to generate the normalized difference vegetation index (NDVI) at a resolution of 250 × 250 m from 100 randomly selected points in each protected area. The final NDVI was calculated using the equation: NDVI = (NIR − RED)/(NIR + RED), where NIR represents the near-infrared electromagnetic ray and RED the visible red ray. Given that NDVI is a standardized method used to evaluate the health status of vegetation by quantifying the difference between the near-infrared electromagnetic ray (which vegetation strongly reflects) and the visible red light (which vegetation absorbs), the formula generates values between − 1 and + 1: negative values indicate the presence of water, while values close to 0 indicate bare soils; values between 0.1 and 0.5 indicate low-to-medium vegetation density cover, while values between 0.5 to + 1 indicate high vegetation density. To generate the NDVI raster, the calculator tool in Q GIS was used. As well, random points were generated within the study area, for which the NDVI values were extracted. This was done because there was an assumption that elephants move within their habitats. The NDVI values generated for the random points were compared to the actual sampled elephant locations for both the families and males recorded in the field.

Statistical analyses

Assigning eggs of strongylid nematodes into taxonomic classes using measurements is a challenge as there is some degree of overlap in the dimensions of the eggs of these taxa. Moreover, when the measurements are multidimensional, discordance in measurements taken in a one dimension can lead to biases when assigning eggs to genera. However, by using model-based clustering, these problems can be overcome as information on the variation in the densities of the measurements across the taxa and the covariance of the different measurements is used to minimize the assignment bias. We employed an unsupervised multivariate cluster model using Gaussian finite mixture analysis to group nematode and trematode eggs into operational taxonomic units (OTUs) using egg measurements. The Gaussian finite mixture model (GMM), assumes that the measurements of helminth eggs taken from each taxon (species or genus) will follow a normal distribution resulting in a (multivariate) Gaussian distribution; each taxonomic component will form a cluster of unique density, centred at the mean vector, and with other geometric features such as volume, shape and orientation of the measurements determined by the covariance matrix. The volume, shape and orientation of the covariance’s can be constrained to be equal or variable across groups, giving rise to 14 possible models characterized by unique geometric characteristics [24]. The most parsimonious parameterisation of the covariance matrix is obtained using eigen-decomposition. The Gaussian finite mixture clustering process provides a model estimate for the data that allows for overlapping clusters and produces a probabilistic clustering that quantifies the uncertainty of observations belonging to the components of the mixture. The unsupervised GMM was performed using the mclust package [24] of the R statistical software [25]. The OTUs of the nematode and trematode eggs were assigned to taxonomic classes of helminth eggs based on mean length and width measurements taken from published records (Additional file 1: Table S1). We present here the data at the generic rather than the species level due to the variability in helminth egg sizes reported in previous studies.

To test our hypotheses, we conducted both bivariate and multivariate analyses to evaluate whether the variation in the dependent covariates such as helminth prevalence and helminth epg are influenced by independent covariates such as NDVI and sociality in the presence of unbalanced data. Any discordance in bivariate and multivariate models suggests that imbalance in data is causing spurious partial covariate effects. Using bivariate analyses and Friedman and Kruskal-Wallis tests, we tested for differences in epg between populations. We conducted multivariate analyses using Poisson and negative binomial generalized linear models (GLM) including hurdle and simple count models. The best model was selected based on parsimony criteria using Akaike information criteria (AIC). To examine the influence of social group type, NDVI and age on epg, we used the negative binomial hurdle GLM with the glmmTMB package [26] in the R statistical software [25].

Results

The best Gaussian finite mixture cluster model for trematodes was one with two components characterized by ellipsoidal, equal volume, shape and orientation (EE2) (Fig. 2). This model revealed that elephant populations in Kenya were infected by trematodes that can be characterized by two OTUs. Based on published egg dimensions of trematodes infecting African elephants, we found that trematode OTU1 had mean egg lengths and widths that were similar to those for Protofasciola robusta, while OTU2 had egg dimensions that were similar to Brumptia bicaudata (Figs. 2, 3; Table 1, Additional file 1: Table S1). For nematodes, the best supported model based on BIC consisted of five diagonal components with equal shape [27] indicating the presence of five OTUs (Fig. 2, Table 1, Additional file 1: Table S1): one group putatively belongs to the genus Murshidia (OTU1) with similar egg measurements to Murshidia dawoodi; three groups belong to the genus Quilonia (OTU2, OTU3 and OTU4) (Fig. 3), with egg measurements very similar to those of Quilonia apiensis (OTU2), Q. africana (OTU3) and Q. magna (OTU4). Finally, OTU5 had large egg measurements went beyond the range for the genus Quilonia but were similar to those recorded for Mammomonogamus loxodontis (Table 1, Additional file 1: Table S1).

Fig. 2
figure 2

Model-based classification into operational taxonomic groups of elephants’ a trematode and b Strongylidae (nematode) eggs. Trematodes are classified into two (OTU1, blue; OTU2, red) and nematodes into five (OTU1, green; OTU2, orange; OTU3, purple; OTU4, red; and OTU5, blue) operational taxonomic units

Full size image
Fig. 3
figure 3

Photomicrographs of eggs of nematodes and trematodes of five genera. aMurshidia (69 × 37 µm). b, c Different egg sizes of Quilonia: 96 × 57 µm (b) and 84 × 52 µm (c). dMammomonogamus (101 × 59 µm). eProtofasciola (90 × 49 µm). fBrumptia (115 × 59 µm). Scale-bars: 50 µm

Full size image
Table 1 Results of unsupervised classification of trematode and nematode eggs into operational taxonomic units (OTUs)
Full size table

The prevalence of infection determined from sedimentation was 97.5%, whereas the prevalence obtained from floatation was 92.6%; this difference, however, was not statistically significant (χ 2(1, n = 243)  = 0.769, P = 0.366; Table 2). Therefore, all analyses of the prevalence were based on results obtained using the sedimentation technique. The prevalence of helminth infection determined using sedimentation varied between populations and was statistically significant (χ 2(3, n = 243)  = 8.972, P = 0.030); however, there was no association between prevalence and elephant social group (male social groups vs female social groups, χ 2(1, n = 243)  = 0.461, P = 0.497). The prevalence of nematodes was 96.3% (95% CI: 93.09–98.29%) and was significantly higher than that of trematodes, which was 39.1% (95% CI: 32.92–45.54%; χ 2(1, n = 243)  = 179.18, P < 0.001). There was no significant influence of either social group (χ 2(1, n = 243)  = 1.952, P = 0.162) or sampling location (χ 2(3, n = 243)  = 5.956, P = 0.114) on nematode prevalence (Table 3). By contrast, trematode prevalence was significantly influenced by the location and elephant population (χ 2(3, n = 243)  = 53.13, P < 0.001, Table 2) but not by social group (χ 2(1, n = 243)  = 0.254, P = 0.614; Table 3).

Table 2 Variation in the prevalence of helminths in elephant populations and social groups in Kenya estimated using sedimentation and floatation methods
Full size table
Table 3 Prevalence of nematodes and trematodes in male and family social groups in different populations estimated using the faecal sedimentation method
Full size table

The quantitative analysis using the McMaster technique revealed that the mean epg varied within elephant populations and between elephant social groups. Bivariate analyses revealed that elephants sampled in family groups had significantly higher median epg than solitary males and/or males in bachelor groups when controlling for epg variation across sampling locations or elephant populations (Friedman χ 2(1, n = 243)  = 4, P = 0.046; Fig. 4). Elephant family groups in the four elephant populations showed significant differences in mean epg (Kruskal-Wallis χ 2(3, n = 243)  = 40.942, P < 0.001; Fig. 4). Similarly, male groups from various elephant populations differed in their mean epg (Kruskal-Wallis χ 2(3, n = 243)  = 9.38, P = 0.025; Fig. 4). Young males in mixed sex family groups had a lower epg than that of females (Table 4) when controlling for population or age; however, the differences were not statistically significant (population: Friedman χ 2(1, n = 243)  = 1, P = 0.317; age: Friedman χ 2(1, n = 243)  = 2, P = 0.157).

Fig. 4
figure 4

Mean egg burden (epg faeces) of helminths for each social group in the studied Kenyan elephant populations

Full size image
Table 4 Mean helminth burden (epg faeces) for each sex and type of social group in Kenyan elephant populations
Full size table

The normalized difference vegetation index (NDVI) was generally very low in the areas occupied by the study populations. The average NDVI over a three-month period varied between the four study locations and these differences were statistically significant (Kruskal-Wallis; χ 2(3, n = 243)  = 9.18, P = 0.027). The lowest three-month mean NDVIs were recorded in Amboseli (mean ± standard deviation, 0.091 ± 0.002) and Tsavo East (0.118 ± 0.006) but were relatively higher in Laikipia-Samburu (0.16 ± 0.013) and Maasai Mara (0.239).

The most parsimonious multivariate model for variation in helminth epg was a hurdle GLM with a negative binomial distribution. This model indicated that the variation in non-zero positive counts of epg were driven by three-month cumulative mean NDVIs, social group type and elephant population, and sampling location (Table 5). We observed a positive association between mean epg and the three-month cumulative mean NDVI (Fig. 5). Among elephant social groups, female social groups had a higher mean epg than male social groups (Fig. 5). Among elephant populations or protected areas, the elephant populations in Samburu-Laikipia had a significantly higher epg than in Amboseli, while in Tsavo elephants had a significantly lower epg than in Amboseli (Fig. 6). However, Maasai Mara elephants showed no differences in their epg from Amboseli. In the binomial part of the model, which shows the presence or absence of a detectable epg, location and age and, to a lesser extent, NDVI all had a significant influence on detectable helminth infection. Laikipia-Samburu elephants had higher helminth prevalence than Amboseli elephants, whereas the elephants from Maasai Mara and Tsavo East had prevalence that were similar to Amboseli (Table 5). Adult elephants had a higher detectable epg than sub-adults and juveniles combined.

Table 5 A multivariate hurdle GLM showing important factors explaining variations in epg between Kenyan elephant populations
Full size table
Fig. 5
figure 5

Scatterplot showing the relationship between NDVI and egg burden (epg faeces) for each social group for all elephant populations combined

Full size image
Fig. 6
figure 6

Scatterplot showing the relationship between NDVI and egg burden (epg faeces) for each social group for elephant populations treated separately

Full size image

Discussion

The African elephant is a mega herbivore and a keystone species in conservation whose ecological impact on the diversity and survival of habitats and other species is enormous [28]. Its numbers have continued to dwindle in many parts of its African range where populations are split into separate subpopulations [29]. As such, subpopulations may suffer different rates of parasite infestation. Here we present the first study to have examined helminth infection patterns in distinct African elephant populations, and the first to have evaluated factors associated with intra- and inter-population variability in prevalence and egg burden (epg). This study used only egg dimensions for helminth identification and the inference of prevalence and load since the obtaining of worms is often invasive or opportunistic. We identified two genera of trematodes, Protofasciola and Brumptia, and three of nematodes, Murshidia, Mammomonogamus and Quilonia (tribe Quiloninea, subfamily Cyathostominae). Using unsupervised classification of strongylid nematodes, we recovered OTUs that had egg dimensions corresponding to species that are known to occur in East Africa, in particular Kenya. Specifically, we recovered eggs from species including P. robusta, B. bicaudata, Murshidia africana, Quilonia africana, Q. magna and Mammomonogamus loxodontis, which demonstrates that, using information on egg measurements, it is possible to employ model-based clustering to group eggs into taxonomic units matching actual species. This type of model-based classification for nematode eggs has previously been tested with discriminations within the range of 80–95% for sheep nematodes [30, 31].

Although cluster-based modelling of the dimensions of strongylid-type eggs was used to identify genera and, potentially, species of the Strongylidae, the utility of OTUs detected by model-based clustering depends on the knowledge of egg measurements from species known to infect elephant populations. Such data will help identify any variation within species from different host populations, thereby providing useful information for egg identification through model clustering. A potential handicap with this method is that egg measurements for a single species taken from different host populations seem to vary greatly. The cause of this variation is not clear but could be due to the misidentification of larval nematodes from hatched eggs or to high inherent variance in nematode egg size that varies between host populations. The factors causing variance in egg measurements within species across host populations could hamper matching model based OTUs with known species. Thus, our species identification remains tentative and should not be taken as conclusive; nevertheless, it reveals the potential diversity in helminths that exists both within and between host populations.

The overall prevalence of helminths in elephants was 97.5%. However, this prevalence was characterized by a significantly higher proportion of nematodes (96.3%) than of trematodes (39.1%). In our study, we modified the sedimentation method commonly used for examining trematodes (which involves examining all of the sediment in a Petri dish under a dissecting microscope) [3]. However, our results were comparable with infection in elephants elsewhere, which suggests that our modifications did not significantly affect the sensitivity of the methodology. For instance, several studies have reported infection patterns in elephants in which nematode prevalence is 2–3 times greater than trematode prevalence. For instance, nematodes in elephant populations in Burkina Faso, West Africa, had an overall prevalence of 97.7% compared to 30.9% for trematodes [32]. A similar pattern has been observed in Botswana, South Africa, where elephants had a 100% prevalence of nematodes and 26% of trematodes [3]. This pattern is not restricted to the African elephant as a comparable prevalence of nematodes (92–96%) has been recorded in Asian elephants [33]. Furthermore, in our analysis, we observed significant inter-population variation in helminth prevalence. This was probably mainly due to trematode prevalence since, unlike nematodes, trematodes require the presence of an intermediate host (a snail) that depends on the presence of a permanent water source. Specifically, elephants from the Laikipia-Samburu ecosystem had the highest trematode prevalence, while elephants from Tsavo East had the lowest. By contrast, Amboseli elephants, that are known to be exposed to permanent water sources [34], had more moderate trematode prevalence than expected when compared to the Laikipia-Samburu population. The factors that determine trematode prevalence may be linked to the environmental variables influencing the abundance of snails, the intermediate hosts of trematodes [35]. As expected, the prevalence of trematodes in Amboseli elephants exposed to permanent water was nearly double that of elephants using the seasonally water-logged Okavango delta, where the prevalence of trematodes was 23% [3].

Mean egg burden was higher in family groups than in male social groups, which contradicts male-biased parasitism known to be associated with both hormonal and behavioural differences often seen in other animals. In cattle, younger animals and males tend to have higher levels of gastrointestinal parasite infection than older and female animals [36]. In most mammals, males exhibit higher infection rates than females (i.e. humans, ungulates, rodents, bats and birds) [27, 37,38,39,40]. The hormone testosterone is associated with immunosuppression in males, leading to greater parasite infection. However, Thurber et al. [1] found no effect of testosterone on parasite burden in male elephants in musth, individuals expected to have the highest level of testosterone. It is unlikely that testosterone will have more immunosuppressive effects on elephants than on other mammal species [3]. The higher parasite infection observed in female than in male elephants suggests that elephant social structures have a significant influence on mean egg burden since group-living exposes group members to higher parasite infection risks than individuals with solitary lifestyles or who form transient associations. This social dichotomy in infection patterns in elephants may be related to habitat use and ranging patterns, which drive the exposure and transmission of parasites such as helminths. In elephants, the ranging patterns of the female-led (matriarch) family groups are predictable as they often remain within reach of water, and long-distance movement is avoided due to the presence of juveniles. Solitary males or bachelor groups, on the other hand, have no such constraints and range over greater distances [19,20,21,22]. Moreover, infectious helminth propagules build up in frequently used habitats, hence family social groups suffer a higher risk of infection [41]. The influence of sociality on egg burden or egg shedding could also be attributed to the sex composition of the group, especially if the effects of male-biased parasitism are taken into account. However, we did not find any significant difference in egg burden between male and female individuals in the family social groups.

Although the influence of social structure on egg burden has been observed in other elephant populations, the egg burden detected by our study was much lower. We detected a mean egg burden of 172 in female groups and 89 in male groups, figures that are much lower than in elephant populations in the Okavango Delta, Botswana, where female groups had a mean egg burden of 1116 and males 529 [3]. A study of male elephants in Etosha National Park in Namibia revealed that in one year, the average strongylid egg burden varied between 1409–2204 in two different years [42]. In addition, previous studies on Rhodesian elephants have recorded much higher mean egg burden reaching 2072 with a range of 512–4382 [5]. Faecal egg counts or egg burden are often used to assess parasite burdens but have inherent pitfalls as they are subject to numerous variables that confound cause-effect relationships [43]. The few studies that have ever correlated egg burden to worm burdens have had variable outcomes [44,45,46,47]. Therefore, we believe that care should be taken with egg burden values we recorded from elephants since they may not correspond to the total worm burdens. A previous study has reported an average of around 30,000 worms per elephant with a range of 3837–105,294 [48]. The factors that influence variations in egg burden are not clear but may include factors intrinsic to these parasites including variations in the life histories of infecting worm species, the number of immature stages, worm sex imbalance and host-environmental factors [49].

Our results show that NDVI, a measure of vegetation productivity, biomass and habitat structure, were variable but generally low in all four studied habitats. However, NDVI was positively correlated with egg burden. Given that NDVI is correlated with environmental variables such as rainfall, soil moisture and habitat structure [50,51,52,53,54], it can both directly and indirectly determine the survival and transmission of infective stages, the maturation of immature worms in hosts, and the shedding rates of eggs by definitive hosts. Moreover, since NDVI is strongly correlated to primary production and the nutritional content of forage [55,56,57,58,59], it can also influence the distribution and abundance of the susceptible hosts and can enhance the heterogeneity of host spatial distribution [60,61,62]. Evidence of increased transmission of nematodes during the rainy season has been reported in a study of African elephants [3]. A study of human helminth infection also found a positive relationship between the prevalence of helminth infection and NDVI [63].

Conclusions

Overall, our study shows that helminth infection in elephants is characterized by statistically significant inter-population variation in prevalence and egg burden. Sociality in elephants did not influence helminth prevalence but did have an influence on egg burden. Given that NDVI significantly varied between the four habitats and was positively correlated with mean egg burden it is likely that NDVI is an important driver of variation in egg burden in elephant populations.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its additional file. Raw data used and/or analysed during the current study are available from the first and corresponding author upon request.

Abbreviations

AIC:

Akaikeʼs information criterion

ANP:

Amboseli National Park

epg:

eggs per gram

GLM:

generalized linear model

GMM:

Gaussian finite mixture model

LSE:

Laikipia-Samburu ecosystem

MMNR:

Maasai Mara National Reserve

OTU:

operational taxonomic unit

TENP:

Tsavo East National Park

NDVI:

normalized difference vegetation index

References

  1. Thurber MI, O’Connell-Rodwell CE, Turner WC, Nambandi K, Kinzley C, Rodwell TC, et al. Effects of rainfall, host demography, and musth on strongyle fecal egg counts in african elephants (Loxodonta africana) in Namibia. J Wildl Dis. 2011;47:172–81.

    Article  CAS  PubMed  Google Scholar 

  2. Kinsella JM, Deem SL, Blake S, Freeman AS. Endoparasites of African forest elephants (Loxodonta africana cyclotis) from the Republic of Congo and Central African Republic. Comp Parasitol. 2004;71:104–10.

    Article  Google Scholar 

  3. Baines L, Morgan ER, Ofthile M, Evans K. Occurrence and seasonality of internal parasite infection in elephants, Loxodonta africana, in the Okavango delta, Botswana. Int J Parasitol Parasites Wildl. 2015;4:43–8.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Condy JB. Observations on internal parasites in Rhodesian elephant, Loxodonta africana Blumenbach, 1797. Proc Trans Rhodesia Sci Assoc. 1974;55:67–99.

    Google Scholar 

  5. Van Der Westhuysen OP. A monograph of the helminth parasites of the elephant. Onderstepoort J Vet Sc Anim Ind. 1938;10:49–190.

    Google Scholar 

  6. Morgan ER, van Dijk J. Climate and the epidemiology of gastrointestinal nematode infections of sheep in Europe. Vet Parasitol. 2012;189:8–14.

    Article  CAS  PubMed  Google Scholar 

  7. Pandey V, Chitate F, Nyanzunda T. Epidemiological observations on gastro-intestinal nematodes in communal land cattle from the highveld of Zimbabwe. Vet Parasitol. 1993;51:99–106.

    Article  CAS  PubMed  Google Scholar 

  8. Pfukenyi DM, Mukaratirwa S. A review of the epidemiology and control of gastrointestinal nematode infections in cattle in Zimbabwe. Onderstepoort J Vet Res. 2013;80:612.

    Article  PubMed  Google Scholar 

  9. Froeschke G, Harf R, Sommer S, Matthee S. Effects of precipitation on parasite burden along a natural climatic gradient in southern Africa—implications for possible shifts in infestation patterns due to global changes. Oikos. 2010;119:1029–39.

    Article  Google Scholar 

  10. Apio A, Plath M, Wronski T. Patterns of gastrointestinal parasitic infections in the bushbuck Tragelaphus scriptus from the Queen Elizabeth National Park, Uganda. J Helminthol. 2006;80:213–8.

    CAS  PubMed  Google Scholar 

  11. Turner WC, Getz WM. Seasonal and demographic factors influencing gastrointestinal parasitism in ungulates of Etosha National Park. J Wildl Dis. 2010;46:1108–19.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Fowler ME, Mikota SK. Biology, medicine, and surgery of elephants. Oxford: Blackwell Publishing Ltd; 2006.

    Book  Google Scholar 

  13. Obanda V, Iwaki T, Mutinda NM, Gakuya F. Gastrointestinal parasites and associated pathological lesions in starving free-ranging African elephants. S Afr J Wildl Res. 2011;41:167–72.

    Article  Google Scholar 

  14. Sherrard-Smith E, Chadwick EA, Cable J. Climatic variables are associated with the prevalence of biliary trematodes in otters. Int J Parasitol. 2013;43:729–37.

    Article  PubMed  Google Scholar 

  15. Wittemyer G, Douglas-Hamilton I, Getz WM. The socioecology of elephants: analysis of the processes creating multitiered social structures. Anim Behav. 2005;69:1357–71.

    Article  Google Scholar 

  16. Chiyo PI, Archie EA, Hollister-Smith JA, Lee PC, Poole JH, Moss CJ, et al. Association patterns of African elephants in all-male groups: the role of age and genetic relatedness. Anim Behav. 2011;81:1093–9.

    Article  Google Scholar 

  17. Moss CJ, Poole JH. Relationships and social structure of african elephants. In: Hinde RA, editor. Primate social relationships: an integrated approach. Oxford: Blackwell Scientific; 1983. p. 315–25.

    Google Scholar 

  18. Archie EA, Moss CJ, Alberts SC. The ties that bind: genetic relatedness predicts the fission and fusion of social groups in wild African elephants. Proc R Soc B. 2006;273:513–22.

    Article  CAS  PubMed  Google Scholar 

  19. Barnes RFW. Mate searching behaviour of elephant bulls in a semi-arid environment. Anim Behav. 1982;30:1217–23.

    Article  Google Scholar 

  20. Shannon G, Page B, Slotow R, Duffy K. African elephant home range and habitat selection in Pongola Game Reserve, South Africa. Afr Zool. 2006;41:37–44.

    Article  Google Scholar 

  21. Leggett KEA. Home range and seasonal movement of elephants in the Kunene region, northwestern Namibia. Afr Zool. 2006;41:17–36.

    Article  Google Scholar 

  22. Mills EC, Poulsen JR, Fay JM, Morkel P, Clark CJ, Meier A, et al. Forest elephant movement and habitat use in a tropical forest-grassland mosaic in Gabon. PLoS ONE. 2018;13:e0199387.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. VanderWaal K, Omondi GP, Obanda V. Mixed-host aggregations and helminth parasite sharing in an East African wildlife-livestock system. Vet Parasitol. 2014;205:224–32.

    Article  PubMed  Google Scholar 

  24. Scrucca L, Fop M, Murphy TB, Raftery AE. mclust 5: clustering, classification and density estimation using Gaussian finite mixture models. R J. 2016;8:289–317.

    Article  PubMed  PubMed Central  Google Scholar 

  25. R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.

  26. Magnusson A, Skaug H, Nielsen A, Berg C, Kristensen K, Maechler M, et al. glmmTMB: generalized linear mixed models using a template model builder. R package version 01. 2017;3.

  27. Guerra-Silveira F, Abad-Franch F. Sex bias in infectious disease epidemiology: patterns and processes. PLoS ONE. 2013;8:e62390.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Western D. The ecological role of elephants in Africa. Pachyderm. 1989;12:42–5.

    Google Scholar 

  29. Chase MJ, Schlossberg S, Griffin CR, Bouché PJ, Djene SW, Elkan PW, et al. Continent-wide survey reveals massive decline in African savannah elephants. PeerJ. 2016;4:e2354.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Sommer C. Digital image analysis and identification of eggs from bovine parasitic nematodes. J Helminthol. 1996;70:143–51.

    Article  CAS  PubMed  Google Scholar 

  31. Christie M, Jackson F. Specific identification of strongyle eggs in small samples of sheep faeces. Res Vet Sci. 1982;32:113–7.

    Article  CAS  PubMed  Google Scholar 

  32. Nakandé A, Belem AMG, Nianogo AJ, Jost C. Parasites gastro-intestinaux des éléphants dans la Réserve Partielle de Pama, Burkina Faso. Pachyderm. 2007;42:22–32.

    Google Scholar 

  33. Vidya TNC, Sukumar R. The effect of some ecological factors on the intestinal parasite loads of the Asian elephant (Elephas maximus) in southern India. J Biosci. 2002;27:521–8.

    Article  CAS  PubMed  Google Scholar 

  34. Kioko J, Kiringe J, Omondi P. Human-elephant conflict outlook in the Tsavo-Amboseli ecosystem, Kenya. Pachyderm. 2006;41:53–60.

    Google Scholar 

  35. Lodge DM, Brown KM, Klosiewski SP, Stein R, Covich A, Leathers B, et al. Distribution of freshwater snails: spatial scale and the relative importance of physicochemical and biotic factors. Am Malacol Bull. 1987;5:73–84.

    Google Scholar 

  36. Gunathilaka N, Niroshana D, Amarasinghe D, Udayanga L. Prevalence of gastrointestinal parasitic infections and assessment of deworming program among cattle and buffaloes in Gampaha District, Sri Lanka. BioMed Res Int. 2018;2018:3048373.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Zuk M, McKean KA. Sex differences in parasite infections: patterns and processes. Int J Parasitol. 1996;26:1009–24.

    Article  CAS  PubMed  Google Scholar 

  38. Martínez-Guijosa J, Martínez-Carrasco C, López-Olvera JR, Fernández-Aguilar X, Colom-Cadena A, Cabezón O, et al. Male-biased gastrointestinal parasitism in a nearly monomorphic mountain ungulate. Parasites Vectors. 2015;8:165.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Connors V, Nickol B. Effects of Plagiorhynchus cylindraceus (Acanthocephala) on the energy metabolism of adult starlings, Sturnus vulgaris. Parasitology. 1991;103:395–402.

    Article  PubMed  Google Scholar 

  40. Krasnov BR, Bordes F, Khokhlova IS, Morand S. Gender-biased parasitism in small mammals: patterns, mechanisms, consequences. Mammalia. 2012;76:1–13.

    Article  Google Scholar 

  41. Hausfater G, Meade BJ. Alternation of sleeping groves by yellow baboons (Papio cynocephalus) as a strategy for parasite avoidance. Primates. 1982;23:287–97.

    Article  Google Scholar 

  42. Brumfitt KN. Assessment of body condition of African elephants (Loxodonta africana) in north east of Etosha National Park, Namibia: how it relates to strongyle parasite egg counts and nutrition value of feed. MSc. thesis, The University of Namibia, Windhoek, Namibia; 2015.

  43. Coulson G, Cripps JK, Garnick S, Bristow V, Beveridge I. Parasite insight: assessing fitness costs, infection risks and foraging benefits relating to gastrointestinal nematodes in wild mammalian herbivores. Philos Trans R Soc B. 2018;373:20170197.

    Article  CAS  Google Scholar 

  44. Cripps J, Beveridge I, Martin JK, Borland D, Coulson G. Temporal dynamics of helminth infections in eastern grey kangaroos (Macropus giganteus) in Victoria. Aust J Zool. 2015;63:163–74.

    Article  Google Scholar 

  45. Seivwright L, Redpath S, Mougeot F, Watt L, Hudson P. Faecal egg counts provide a reliable measure of Trichostrongylus tenuis intensities in free-living red grouse Lagopus lagopus scoticus. J Helminthol. 2004;78:69–76.

    Article  CAS  PubMed  Google Scholar 

  46. Gassó D, Feliu C, Ferrer D, Mentaberre G, Casas-Díaz E, Velarde R, et al. Uses and limitations of faecal egg count for assessing worm burden in wild boars. Vet Parasitol. 2015;209:133–7.

    Article  PubMed  Google Scholar 

  47. Kim J-H, Choi M-H, Bae YM, Oh J-K, Lim MK, Hong S-T. Correlation between discharged worms and fecal egg counts in human clonorchiasis. PLoS Negl Trop Dis. 2011;5:e1339.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Condy JB. Observations on internal parasites in Rhodesian elephants Loxodonta africana Blumenbach 1797. Proc Trans Rhod Sci Assoc. 1974;55:67–99.

    Google Scholar 

  49. McKenna P. The diagnostic value and interpretation of faecal egg counts in sheep. N Z Vet J. 1981;29:129–32.

    Article  CAS  PubMed  Google Scholar 

  50. Gamon JA, Field CB, Goulden ML, Griffin KL, Hartley AE, Joel G, et al. Relationships between NDVI, canopy structure, and photosynthesis in three Californian vegetation types. Ecol Appl. 1995;5:28–41.

    Article  Google Scholar 

  51. Nicholson SE, Farrar TJ. The influence of soil type on the relationships between NDVI, rainfall, and soil moisture in semiarid Botswana. I. NDVI response to rainfall. Remote Sens Environ. 1994;50:107–20.

    Article  Google Scholar 

  52. Farrar TJ, Nicholson SE, Lare AR. The influence of soil type on the relationships between NDVI, rainfall, and soil moisture in semiarid Botswana. II. NDVI response to soil moisture. Remote Sens Environ. 1994;50:121–33.

    Article  Google Scholar 

  53. Davenport ML, Nicholson SE. On the relation between rainfall and the normalized difference vegetation index for diverse vegetation types in East Africa. Int J Remote Sens. 1993;14:2369–89.

    Article  Google Scholar 

  54. Li J, Lewis J, Rowland J, Tappan G, Tieszen LL. Evaluation of land performance in Senegal using multi-temporal NDVI and rainfall series. J Arid Environ. 2004;59:463–80.

    Article  Google Scholar 

  55. Paruelo JM, Epstein HE, Lauenroth WK, Burke IC. ANPP estimates from NDVI for the central grassland region of the United States. Ecology. 1997;78:953–8.

    Article  Google Scholar 

  56. Prince SD. Satellite remote sensing of primary production: comparison of results for Sahelian grasslands 1981–1988. Int J Remote Sens. 1991;12:1301–11.

    Article  Google Scholar 

  57. Tucker CJ, Vanpraet CL, Sharman MJ, Van Ittersum G. Satellite remote sensing of total herbaceous biomass production in the Senegalese Sahel: 1980–1984. Remote Sens Environ. 1985;17:233–49.

    Article  Google Scholar 

  58. Albayrak S. Use of reflectance measurements for the detection of N, P, K, ADF and NDF contents in sainfoin pasture. Sensors. 2008;8:7275–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Ryan SJ, Cross PC, Winnie J, Hay C, Bowers J, Getz WM. The utility of normalized difference vegetation index for predicting African buffalo forage quality. J Wildl Manag. 2012;76:1499–508.

    Article  Google Scholar 

  60. Sousa WP, Grosholz ED. The influence of habitat structure on the transmission of parasites. In: Bell SS, McCoy ED, Mushinsky HR, editors. Habitat structure; the physical arrangement of objects in space. Suffolk: Springer Science & Business Media; 1991. p. 300–24.

    Chapter  Google Scholar 

  61. Mueller T, Olson KA, Fuller TK, Schaller GB, Murray MG, Leimgruber P. In search of forage: predicting dynamic habitats of Mongolian gazelles using satellite-based estimates of vegetation productivity. J Appl Ecol. 2008;45:649–58.

    Article  Google Scholar 

  62. Leimgruber P, McShea WJ, Brookes CJ, Bolor-Erdene L, Wemmer C, Larson C. Spatial patterns in relative primary productivity and gazelle migration in the Eastern Steppes of Mongolia. Biol Conserv. 2001;102:205–12.

    Article  Google Scholar 

  63. Lai Y-S, Zhou X-N, Utzinger J, Vounatsou P. Bayesian geostatistical modelling of soil-transmitted helminth survey data in the People’s Republic of China. Parasites Vectors. 2013;6:359.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to the Director of the Kenya Wildlife Service for permission to carry out this study and to the Department of Veterinary Service for support with the infrastructure. We are specifically grateful to Nelson Mwangi, Amos Muthiuru, Felix Micheni, Vasco Nyaga, Humphrey Kinyua, Stephen Nyaga, Joseph Kyalo, Lilian Apollo, Christine Mwende and Katana for helping with the sampling in the various protected areas, and to Dr Ndambiri for his support in Amboseli National Park.

Funding

This research was funded by the Kenya Wildlife Service and the Department of Animal Pathology, Veterinary Faculty, University of Santiago de Compostela, Lugo, Spain.

Author information

Authors and Affiliations

  1. Department of Animal Pathology (INVESAGA Group), Veterinary Faculty, University of Santiago de Compostela, Lugo, Spain

    Edward King’ori & Patrocinio Morrondo

  2. Veterinary Department, Kenya Wildlife Service, Nairobi, Kenya

    Edward King’ori & Vincent Obanda

  3. Institute of Primate Research, National Museums of Kenya, Nairobi, Kenya

    Patrick I. Chiyo

  4. Estación Biológica de Doñana, Consejo Superior de Investigaciones Científicas (CSIC), Sevilla, Spain

    Ramon C. Soriguer & Samer Angelone

  5. Institute of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, Switzerland

    Samer Angelone

Authors
  1. Edward King’ori
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Vincent Obanda
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Patrick I. Chiyo
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ramon C. Soriguer
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Patrocinio Morrondo
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Samer Angelone
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

EK, VO and PM developed the research concept and design. EK and VO carried out the fieldwork and laboratory analysis. EK, VO, PIC, PM, RS and SA performed the data analysis, interpretation, drafting and critical reviews in the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Samer Angelone.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Research and Ethics Committee of the Kenya Wildlife Service (permit number: KWS/BRM/5001), which is mandated to protect and conserve wildlife in Kenya. The study employed non-invasive methods to study helminth infection patterns and observational methods to identify the sex and age of subjects of the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Table S1.

Egg measurements (in µm) of gastrointestinal nematodes and trematodes infecting African elephants (data compiled from the literature).

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

King’ori, E., Obanda, V., Chiyo, P.I. et al. Patterns of helminth infection in Kenyan elephant populations. Parasites Vectors 13, 145 (2020). https://doi.org/10.1186/s13071-020-04017-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-020-04017-1

Keywords

Parasites & Vectors

ISSN: 1756-3305

Contact us