Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Fecal parasite risk in the endangered proboscis monkey is higher in an anthropogenically managed forest environment compared to a riparian rain forest in Sabah, Borneo

  • Annette Klaus,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    Current address: Veterinary Medicine & Parasitology Department, Springer Heidelberg, Heidelberg, Baden-Wuerttemberg, Germany

    Affiliations Institute of Zoology, University of Veterinary Medicine Hannover, Hannover, Lower Saxony, Germany, Institute for Parasitology, Centre for Infection Medicine, University of Veterinary Medicine Hannover, Hannover, Lower Saxony, Germany

  • Christina Strube,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – review & editing

    Affiliation Institute for Parasitology, Centre for Infection Medicine, University of Veterinary Medicine Hannover, Hannover, Lower Saxony, Germany

  • Kathrin Monika Röper,

    Roles Data curation, Writing – review & editing

    Current address: Zoology Department, Zoo Hannover, Hannover, Lower Saxony, Germany

    Affiliation Institute of Zoology, University of Veterinary Medicine Hannover, Hannover, Lower Saxony, Germany

  • Ute Radespiel,

    Roles Conceptualization, Formal analysis, Methodology, Validation, Writing – review & editing

    Affiliation Institute of Zoology, University of Veterinary Medicine Hannover, Hannover, Lower Saxony, Germany

  • Frank Schaarschmidt,

    Roles Formal analysis, Methodology, Software, Writing – review & editing

    Affiliation Institute for Biostatistics, Leibniz University, Hannover, Lower Saxony, Germany

  • Senthilvel Nathan,

    Roles Project administration, Writing – review & editing

    Affiliation Sabah Wildlife Department, Kota Kinabalu, Sabah, Malaysia

  • Benoit Goossens,

    Roles Conceptualization, Project administration, Writing – review & editing

    Affiliations Sabah Wildlife Department, Kota Kinabalu, Sabah, Malaysia, Organisms and Environment Division, School of Biosciences, Cardiff University, Cardiff, Wales, United Kingdom, Danau Girang Field Centre, c/o Sabah Wildlife Department, Kota Kinabalu, Sabah, Malaysia, Sustainable Places Research Institute, Cardiff University, Cardiff, Wales, United Kingdom

  • Elke Zimmermann

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – review & editing

    elke.zimmermann@tiho-hannover.de

    Affiliation Institute of Zoology, University of Veterinary Medicine Hannover, Hannover, Lower Saxony, Germany

Abstract

Understanding determinants shaping infection risk of endangered wildlife is a major topic in conservation medicine. The proboscis monkey, Nasalis larvatus, an endemic primate flagship species for conservation in Borneo, is endangered through habitat loss, but can still be found in riparian lowland and mangrove forests, and in some protected areas. To assess socioecological and anthropogenic influence on intestinal helminth infections in N. larvatus, 724 fecal samples of harem and bachelor groups, varying in size and the number of juveniles, were collected between June and October 2012 from two study sites in Malaysian Borneo: 634 samples were obtained from groups inhabiting the Lower Kinabatangan Wildlife Sanctuary (LKWS), 90 samples were collected from groups of the Labuk Bay Proboscis Monkey Sanctuary (LBPMS), where monkeys are fed on stationary feeding platforms. Parasite risk was quantified by intestinal helminth prevalence, host parasite species richness (PSR), and eggs per gram feces (epg). Generalized linear mixed effect models were applied to explore whether study site, group type, group size, the number of juveniles per group, and sampling month predict parasite risk. At the LBPMS, prevalence and epg of Trichuris spp., strongylids, and Strongyloides spp. but not Ascaris spp., as well as host PSR were significantly elevated. Only for Strongyloides spp., prevalence showed significant changes between months; at both sites, the beginning rainy season with increased precipitation was linked to higher prevalence, suggesting the external life cycle of Strongyloides spp. to benefit from humidity. Higher prevalence, epgs, and PSR within the LBPMS suggest that anthropogenic factors shape host infection risk more than socioecological factors, most likely via higher re-infection rates and chronic stress. Noninvasive measurement of fecal parasite stages is an important tool for assessing transmission dynamics and infection risks for endangered tropical wildlife. Findings will contribute to healthcare management in nature and in anthropogenically managed environments.

Introduction

Parasite infection patterns are shaped by multifaceted variables. Understanding the drivers of parasitism is a main goal of conservation medicine [1] and crucial to assess the impact of parasites and diseases on endangered species [2, 3]. Especially nonhuman primate populations are increasingly threatened by massive anthropogenic changes of natural habitats. These changes have recently been shown to also induce alterations of parasite risks on nonhuman primate reproduction and survival [4, 5, 6].

Anthropogenic influence was described to promote helminth infection measures indirectly via different modes of intervention, such as provisioning of food resources leading to increased host contact and parasite transmission rates [7, 8]. In natural forests affected by fragmentation, host population density was documented to correlate positively with parasite richness and prevalence [9]. In different forest types, it was shown that the risk of infection may rise with the presence of more social group partners within a group [10], but that larger groups may also disrupt direct transmission dynamics through fission as a counter-strategy [11, 12]. Additionally, host age and developing immunity [13] as well as co-infection status [14] may affect parasite risk. Furthermore, hormonally mediated sex differences may lead to divergent infection patterns in different types of social groups like harem (one male-multi female) and bachelor (all male) units, as immunosuppressive steroid hormone levels were stated to cause a male bias in parasitism [15]. Moreover, it was shown that climatic changes and seasonal challenges in moist and dry tropical forests shape parasite risk; increased precipitation may promote survival of specific parasite intermediates [16] whereas drier habitat conditions may influence host susceptibility towards parasite infections [17].

Borneo is known to harbor a vast number of endemic plant and animal species threatened by deforestation [18]. The proboscis monkey (Nasalis larvatus) is a primate flagship species for conservation in Sabah [19, 20]. Linked to lowland riverine and coastal mangrove forests, proboscis monkeys inhabit a habitat type most affected by logging and agricultural land use [21]. Social groups consist of up to 30 individuals and are organized as either harems (one male-multi female-groups with their offspring) or bachelor (all male) units [22]. Since the 1990s, proboscis monkeys became a major focus for ecotourism [23]. The Lower Kinabatangan Wildlife Sanctuary (LKWS) is a protected area open for tourism and some of the well-known wildlife tourism projects include the Labuk Bay Proboscis Monkey Sanctuary (LBPMS) [24]. The LKWS is situated along the floodplain of the Kinabatangan River in Sabah, Malaysian Borneo, and consists of 10 blocks of remaining forest, so called ‘Lots’, along the river (Fig 1). Most of the area surrounding the LKWS was converted to oil palm plantations during the last decades. Local proboscis monkey groups (Table 1) can range between adjacent Lots or cross the Kinabatangan [25]. The LBPMS comprises private land near Samawang Village at Labuk Bay, Sandakan, in the northeast of Malaysian Sabah. It consists of a stretch of mangrove forest and coastal swamps along the coastline [26] and is further surrounded by agriculture (Fig 1). The anthropogenically managed sanctuary shelters a population of proboscis monkeys organized in six harems and three bachelor groups (personal communication of rangers in 2012) (Table 1).

thumbnail
Fig 1. Sample collection sites in Sabah, Malaysian Borneo, Southeast Asia.

(a) Map with Borneo, Southeast Asia, and the position of sampling sites framed. (b) Enlarged map with Sabah, Malaysian Borneo, and the position of the two sampling sites: the LKWS with Lots 1–10 (dark grey) along the Kinabatangan River (black line) with DGFC (white house) in Lot 6, and the LBPMS near the coastline (grey house). LKWS = Lower Kinabatangan Wildlife Sanctuary, LBPMS = Labuk Bay Proboscis Monkey Sanctuary, DGFC = Danau Girang Field Centre.

https://doi.org/10.1371/journal.pone.0195584.g001

thumbnail
Table 1. Local proboscis monkey density across two sample collection sites in Sabah, Malaysian Borneo.

https://doi.org/10.1371/journal.pone.0195584.t001

Noninvasive rapid assessment tools for estimating parasite risks are of imminent importance for proboscis monkeys, since parasites may be a risk for survival of this endangered monkey species. Here, we explore whether the noninvasive collection of fecal samples can be used to assess the relative impact of anthropogenic (site of origin) and ecological parameters (season = time of sampling/month of the year; group size; group type = harem or bachelor, and group demography = proportion of juveniles per group) on parasite risk in the LKWS and the LBPMS. Host parasite risk was quantified by parasite prevalence, parasite species richness, and egg shedding intensity. The following hypotheses will be explored: (1) study site has a significant effect on parasite risk: parasite prevalence, species richness, and egg shedding intensity are higher in the forest area with feeding platforms and higher host density (LBPMS); (2) group size, group type or demography will shape parasite risk: larger groups or bachelor groups have a generally higher parasite prevalence, species richness, and epgs, whereas groups with more juveniles show a higher prevalence of typical juvenile parasitoses; (3) season will affect parasite risk specifically according to parasite transmission cycles: prevalence of parasite species which benefit from humidity increases with rainfall in particular months of the year.

Materials and methods

Sampling sites

In Sabah, Malaysian Borneo, mean temperature variation over the year varies monthly between 21° and 34°C [29]. The annual precipitation inland is 2,800 mm on average [30] and is influenced by two yearly raining seasons; the stronger raining season usually lasts from October to March and the weaker raining season from April to September (Fig 2) [31].

thumbnail
Fig 2. Climate diagram of 2012 based on data from the Danau Girang Field Centre in Lot 6 of the Lower Kinabatangan Wildlife Sanctuary (LKWS) in Sabah.

https://doi.org/10.1371/journal.pone.0195584.g002

Lower Kinabatangan Wildlife Sanctuary (LKWS).

Sampling was conducted in Lot 6 of the government-controlled conservation area of the LKWS (5° 24.231’-5° 24.985’N; 117° 59.682’-118° 04.084’E) and always started at the Danau Girang Field Centre (DGFC) (Fig 1). Using ACME Planimeter measuring the area in Google Maps, the estimated size of Lot 6 was determined as 1,870 ha (http://www.acme.com/planimeter, retrieved 3rd November 2016). According to repeated surveys close to DGFC and a nearby tributary, about 15 groups of proboscis monkeys were inhabiting Lot 6 in 2010 [28]. For the present study, boat surveys for the localization of groups and subsequent sample collection were conducted 4.7 km downstream (3.0 km linear distance) and 6.1 km upstream (3.7 km linear distance) in 2012 along the south river bank starting from DGFC.

Labuk Bay Proboscis Monkey Sanctuary (LBPMS).

Sample collection took place near feeding platforms on the private grounds of the LBPMS close to Samawang Village at Labuk Bay, Sandakan (Fig 1). Here, at two public feeding platforms A (5° 56.130’N; 117° 48.570’E) and B (5° 56.188’N; 117° 47.880’E), free-ranging groups are offered supplementary food, namely man-made cake and cucumber [32] by staff of the LBPMS twice a day. Proboscis monkey groups visited these artificial feeding platforms regularly, most often shortly before being fed.

Sample collection

During the end of the dry and beginning of the rainy season from June to September 2012, a total of 634 fecal samples (6-56/group/day) were obtained in the LKWS, originating from 16 sampling events from harem groups and four sampling events from bachelor groups (Table 2). The weather conditions, in particular storms and rainfall, strongly influenced the feasibility of boat excursions, localization of groups, and successful sample collection and are responsible for the varying number of samples per month. Before sunset (17:00 h– 19:00 h local time), groups near the river bank were localized because proboscis monkeys commonly return to riverside trees for sleeping [33]. Until group members fell asleep, data on group type, group size, the number of juvenile group members, and GPS coordinates were documented. On the next day before sunrise (5:00 h– 8:00 h local time), the same group was approached again to collect as many fresh fecal samples per group per day as possible. Immediately after the animals woke up and foraged deeper into the forest, collection of samples took place. Samples obtained from different days were treated as disparate sampling units. Home ranges of proboscis monkey groups can extensively overlap [34]. Therefore, individual identification of groups was uncertain. Nevertheless, based on our observations of group sizes and group composition, the distribution of groups along the Kinabatangan River, and published estimates of home ranges [35], we assume that we sampled at least nine different harem groups and three different bachelor groups. Therefore, we expect that most of the daily group samples stemmed from different social groups and only a minority of group samples would constitute resampling events.

thumbnail
Table 2. Distribution of the number of monthly fecal samples collected from proboscis monkeys at two sampling sites in Malaysian Borneo in 2012.

https://doi.org/10.1371/journal.pone.0195584.t002

Between 28th September and 18th October 2012, 12 days of feces collection were conducted at the LBPMS. In total, 90 fecal samples (1-18/group/day) were collected during seven sampling events from three harem groups and during four sampling events from two bachelor groups (Table 2) which could be individually identified with the help of the local rangers. Fecal samples were taken after five feeding sessions on platform A (9:30 h local time) and six feeding sessions on platform B (16:30 h local time) when the animals did leave the platforms and returned into the mangroves. During each feeding session observations, data were collected on group type, group size, and the number of juvenile group members.

Treatment of fecal samples and sample analyses

In order to impede a contamination with parasite stages from the ground, the surface of each fecal sample was removed with the help of a spatula. Collected in a 15 mL tube, each fecal sample was weighed and preserved in 10% buffered formalin in a 1:3 ratio (for 1 L: 900 mL distilled water, 100 mL formaldehyde 37–40%, 6.5 g anhydrous disodium hydrogen phosphate Na2HPO4, 4.0 g sodium dihydrogen phosphate NaH2PO4). Analysis of fecal samples was performed using a modified sodium nitrate flotation method [36]. Sodium nitrate solution was made of 600 g sodium nitrate (NaNO3) dissolved in 1 L distilled water for a specific gravity of 1.3, which was steadily verified with a hydrometer (Hydrometer g/mL Tp. 20°C 55mN/m ST. N° 0310, ALLA France). At first, the preservative formalin was largely discarded from the feces by centrifugation at 700 x g for 10 minutes (Heraeus Multifuge 3 S-R, Thermo Fisher Scientific, Waltham, MA, USA). Supernatant was removed and feces were washed through homogenization in distilled water and centrifuged a second time. Subsequently, from each sample a maximum of two gram of feces were mixed with sodium nitrate filling up the sample tube up to 15 mL and centrifuged at 860 x g for 10 minutes. With the help of an eyelet, the whole liquid surface was transferred onto a slide and covered with a coverslip. With an Axiophot microscope (Carl Zeiss MicroImaging, Oberkochen, Germany) each sample was instantly examined microscopically at x10 and x40 magnification. To determine the number of eggs per gram feces, all detected parasitic stages were counted. Using an Olympus ColorView IIIu digital camera (OLYMPUS Soft Imaging Solutions, Munster, Germany), detected parasitic stages (eggs only) were photographed and subsequently measured with the software cell^ B (version 3.1, Olympus Soft Imaging Solutions). In order to avoid bias during investigations, analyses of samples were conducted in a randomized order.

Data analyses

Three parameters of parasite risk were analyzed and tested as dependent variables; (1) parasite prevalence (number of infected individuals as percentage of the whole number of examined animals), (2) host parasite species richness (number of different parasite species morphotypes found in a single host, PSR), and (3) host parasite egg shedding intensity (number of helminth eggs per gram collected feces, epg). Due to the extremely low observed prevalence (0.2%), Trichuris morphotype T4 was excluded from statistical analyses. In total, 724 fecal samples were included in the analyses; 634 samples were from the LKWS and 90 samples from the LBPMS.

To investigate variables influencing prevalence of nematode infections and PSR, two-step generalized linear mixed models (GLMMs) were calculated for each helminth group (a) using binomial presence/absence data (for prevalence) and (b) under the assumption of a Poisson distribution (for PSR). To control for the possibility that wild groups were sampled repeatedly, sampling events (= sampled groups of proboscis monkeys) were set as a random factor in each model. Model comparison and selection were conducted via the anova function (analysis of variance) using the packages lme4, version 1.1–10 [37] and effects [38] in the R software environment (version 0.97.551, R Core Team). In preliminary analyses, all independent variables were first tested separately. In the next step of analyses, the influence of sample weight and sampling site (LKWS vs. LBPMS) on parasite risk were tested. Sample weight might affect detection rates of parasite stages shed by species with low levels of egg shedding, e.g. Strongyloides spp. [39]. Furthermore, sample weight could have a diluting effect on high numbers of egg shedding, e.g. of Trichuris spp. [40]. If sample weight had a significant influence on the dependent variable, it was retained in the following models; otherwise it was no longer included. If the sampling site had a significant effect, the subsequent models were run separately for each sampling site; otherwise one joint model was calculated. In the next step, full models were built that contained at least four fixed factors that could be of biological relevance: (1) group type (harem/bachelor), (2) group size, (3) proportion of juveniles per group (0; <0.3; >0.3), and (4) time of sampling = month of the year (June; July; August; September; October). In addition, an interaction between group size and the proportion of juveniles was also included for typical juvenile parasitoses, such as Strongyloides spp. and Ascaris spp. infections [41] in a third step. Likewise, the interaction term was included in statistical models for PSR. If the inclusion of this term improved the model, it was retained; if not, the interaction term was excluded from the model. Whenever the independent variable (4) time of sampling was a significant predictor variable in analyses for samples collected from the LKWS, all months were compared with a post hoc Tukey-test using the package multcomp [42] to identify significant pairwise differences between sampling months.

Epg data were log(1+x)-transformed and analyzed via two-step linear mixed effect models accordingly, with (1) sampling site, (2) group type, (3) group size, (4) proportion of juveniles per group as well as (5) month of the year as fixed effects. Analysis of model residuals and random effect estimates in QQ-plots showed that a normal distribution on the log-scale was tenable only for most abundant taxonomic groups (Trichuris spp. and strongyles). For the rarer taxonomic groups residuals were clearly right-skewed after log-transformation (Strongyloides spp., Ascaris spp., Enterobius spp., and Anatrichosoma spp.).

Ethics statement

The authors confirm that they did not interact with or disrupt the proboscis monkeys in any way. Fecal sample collection was performed noninvasively from free-ranging individuals without direct animal contact or disturbing the animals and adhered to the Code of Best Practices for Field Primatology of the International Primatological Society (IPS), the ethical guidelines of the German Primate Society (GfP), and the German Animal Protection Act. Sampling in this study required approval and was authorized by the Sabah Wildlife Department, the Sabah Biodiversity Council, the Danau Girang Field Centre as well as the Labuk Bay Proboscis Monkey Sanctuary and complied with Malaysia’s law on foreign research (Licence N° UPE: 40/200/19/2822).

Results

General parasite diversity across both study sites

Overall, 10 helminth egg morphotypes were identified in fecal samples collected from both sites, the LKWS and the LBPMS; namely four different morphotypes of Trichuris spp. (T1-T4) as well as eggs from Anatrichosoma spp., Trichostrongylus spp. [43], Oesophagostomum/Ternidens spp. [41, 44], an unknown strongylid, Strongyloides spp., and Ascaris lumbricoides [45]. In addition, Enterobius spp. eggs were detected in samples from the LKWS [14].

Parasite prevalence

Overall parasite prevalence was linked to parasite taxa and ranged between 1.46 and 91.11% across study sites (Table 3). For trichurids, in particular Trichuris morphotype T3 as well as for Oesophagostomum/Ternidens spp. and Strongyloides spp., but not for the other helminth taxa, prevalence was significantly higher at the LBPMS than at the LKWS (Fig 3, S1 Table). Sample weight had a significant effect on the prevalence of trichurids, namely Trichuris spp. T1-T3, the unknown strongylid, and Strongyloides spp. and was therefore included as a fixed effect in respective statistical models (S1 Table). No significant effect of the variables group type, group size, and the proportion of juvenile group members on helminth prevalence was observed.

thumbnail
Table 3. Prevalence of different helminth egg morphotypes in fecal samples from proboscis monkeys at two sampling sites in Malaysian Borneo.

https://doi.org/10.1371/journal.pone.0195584.t003

thumbnail
Fig 3. Prevalence of intestinal helminth parasites found in proboscis monkey feces at the two sampling sites.

LKWS = Lower Kinabatangan Wildlife Sanctuary, LBPMS = Labuk Bay Proboscis Monkey Sanctuary. * p≤0.05; ** p≤0.001; *** p≤0.0001.

https://doi.org/10.1371/journal.pone.0195584.g003

The time of sampling (month of the year) had no significant effect on the prevalence of helminth species at the two sampling sites, except for Strongyloides spp., in which prevalence significantly increased from June on towards September/October at both sampling sites (S1 Table).

Parasite species richness (PSR)

Parasite species richness was significantly higher in the LBPMS than in the LKWS (Fig 4). Among fecal samples analyzed for different egg morphotypes of helminths (596), 131 samples contained eggs of one morphotype (22.0%), 172 samples contained eggs of two morphotypes (28.9%), 158 samples had eggs of three morphotypes (26.5%), 91 showed eggs of four morphotypes (15.3%), 33 samples contained eggs of five morphotypes (5.5%), and eleven samples had eggs of six morphotypes (1.8%). Overall mean PSR was 2.6 (± 1.2 SD). For the variables group type, group size, and the proportion of juvenile group members no significant effect on PSR was identified (S2 Table).

thumbnail
Fig 4. Parasite species richness of proboscis monkeys at two sampling sites in Sabah.

LKWS = Lower Kinabatangan Wildlife Sanctuary, LBPMS = Labuk Bay Proboscis Monkey Sanctuary. Band in boxes is the median, boxes define the 25th and 75th percentiles, whiskers extend to maximum +/- 1.5 times the interquartile range (IQR = middle 50% of the records). * p≤0.05; ** p≤0.001; *** p≤0.0001.

https://doi.org/10.1371/journal.pone.0195584.g004

Egg shedding intensity

Egg shedding intensity (eggs per gram feces epg) was linked to helminth taxa and sampling site had a significant effect on the epg values of Trichuris sp. T3, Oesophagostomum/Ternidens spp., and Strongyloides spp. but not on the other helminth taxa (Fig 5, S3 Table). In all cases, egg shedding intensity was significantly higher in the LBPMS than in the LKWS. Neither the variables group type, group size, the proportion of juvenile group members nor the time of sampling significantly influenced epg values.

thumbnail
Fig 5. Eggs per gram feces (epg) of positive Oesophagostomum/Ternidens spp., Strongyloides spp., and Trichuris sp. T3 fecal samples from proboscis monkeys at two sampling sites in Sabah.

LKWS = Lower Kinabatangan Wildlife Sanctuary, LBPMS = Labuk Bay Proboscis Monkey Sanctuary. Band in boxes is the median, boxes define the 25th and 75th percentiles, whiskers extend to maximum +/- 1.5 times the interquartile range (IQR = middle 50% of the records). * p≤0.05; ** p≤0.001; *** p≤0.0001.

https://doi.org/10.1371/journal.pone.0195584.g005

Discussion

Our study shows for the first time that fecal collection provides an important noninvasive tool to assess parasite risk of endangered wildlife in different tropical forest settings of Borneo using the proboscis monkey.

As expected, study site was the strongest predictor variable of the intestinal parasite risk in this large bodied primate. This finding suggests that an anthropogenically managed forest area may facilitate intestinal parasite transmission between hosts. Several factors may explain this difference between study sites. Studies on various nonhuman primate species showed that host density is a key determinant for parasite infections [5, 46] and that human provisioning of wildlife with food further affects parasite spread through induced host aggregation and proximity [47]. Based on published estimates for proboscis monkeys along the LKWS [28], and own assessments at the LBPMS, the population density with about 1–1.2 individuals/ha is more than 10 times higher than in Lot 6 of the LKWS with 0.09 individuals/ha. An increased overall host density may cause an accumulation of infective soil-transmitted parasite stages in the environment, which favors the outbreak of significantly more and stronger (re-)infections [48, 49]. The concentrated feeding regime in the LBPMS which is based on a small number of feeding platforms is furthermore increasing encounter rates and the potential exchange of parasite stages via the joint use of the same substrate. The enhanced prevalence of trichurid morphotypes, Oesophagostomum/Ternidens spp., and Strongyloides spp. at the LBPMS goes hand in hand with significantly increased numbers of eggs per gram feces (epg) for Trichuris sp. T3, Oesophagostomum/Ternidens spp., and Strongyloides spp. as well as a significantly higher individual parasite species richness (PSR). An enhanced host density can also be associated with a higher competition for food, sleeping sites, and mating partners as well as frequent contact to visitors and staff, and thus may cause stress, immunosuppression [50] and a higher risk for worm infection. Our findings coincide with those on Thai long-tailed macaques (Macaca fascicularis) which showed increased Strongyloides spp. prevalence and epgs in an anthropogenically modified environment, correlated with higher portions of human-provided food and the time spend foraging on the ground [8]. Likewise, in other wildlife like wild boar (Sus scrofa), anthropogenic management with supplemental feeding was proven to drive the probability and intensity of gastrointestinal parasite infections as well as parasite species richness [51]. Future studies should additionally focus on dietary differences between the proboscis monkey populations. It is well known for mammals that a varying nutritional uptake affects the structure of the gut’s microbiome [52], which itself was recently shown to interact with the intestinal parasite fauna [53]. Ingredients of man-made cake which is supplementary fed at the LBPMS, like carbohydrates from wheat and rice flour [32], may thus be potentially related to the observed parasite infection patterns as well.

In the literature, ecological factors like social group size [54], group demography, and seasonal climatic variations [55] are discussed to predict parasite risk. Here, we considered the potential impact of several ecological parameters on parasite risk for proboscis monkeys at the LKWS and the LBPMS with the help of generalized linear mixed effect models. No effect of the variable group size was observed in any of the parasite types present in this study. Results could potentially be explained by social dynamics, as fission-fusion patterns of proboscis monkey groups could lead to temporally changing group densities or sizes [56], masking potential effects of individual group sizes on infection risk. However, our findings correspond to results on red colobus (Procolobus rufomitratus) in which gastrointestinal helminth infections did not increase with group size either [11]. Additional studies on the relationship between social group size and dynamics, spatial spread and parasite infestation in proboscis monkeys are desirable to elucidate actual transmission pathways and behavioral strategies of differently sized groups. Further, our data did not support the expectation that the number of juveniles in social groups (= proportion of juveniles) has an effect on parasite risk in those parasites known to predominantly affect juveniles (here: Strongyloides and Ascaris spp.). It is possible, though, that not all group members were equally represented in collected fecal samples, and that feces of smaller juveniles may not have been spotted and sampled as readily as those of adults due to age-related fecal size differences or lower fecal production rate in juveniles. In other nonhuman primates like olive baboons (Papio cynocephalus anubis) and Japanese macaques (Macaca fuscata), Strongyloides spp. were described to be more prevalent among juveniles [57, 58]. Furthermore, group type was no significant predictor variable of parasite risk in the present study. This result aligns with missing effects by the proportion of juvenile group members on parasite risk, as those are only found in harem groups. Overall, our results do support the increasing number of studies [10, 59, 60] disproving a generalized effect of increased susceptibility of male hosts (or bachelor groups) to parasite infection due to immunosuppressive effects of the sex steroid testosterone [61].

Moreover, we investigated whether seasonal climatic changes shape parameters of parasite risk at both study sites. This study revealed temporal prevalence dynamics only for excretion of Strongyloides eggs. A specific nature of this genus is a generational change, which comprises an external non-parasitic and an internal parasitic cycle. From eggs shed in the environment, both infective larvae as well as free-living, fertile adult threadworms can develop. The external life cycle benefits from increasing humidity by a higher cumulative survival of adult generations and larvae in the environment which may then result in higher (re-)infection rates. Strongyloides prevalence was shown to increase with enhanced amounts of precipitation in different wildlife species like eastern chimpanzees (Pan troglodytes schweinfurthii) [62] or springbok populations (Antidorcas marsupialis) [63]. Additionally, higher temperatures rather favor generation of adult threadworms in the environment, multiplying the occurrence of infective larvae. During sampling months, monthly precipitation increased and monthly mean temperature rose towards the end of the year (Fig 2). In accordance, the prevalence of Strongyloides spp. significantly increased towards the year’s end. Ongoing studies comprising complete annual seasons at different sampling sites would be desirable to link different parameters of parasite risk to climatic patterns. Summing up, our findings revealed that site of origin outweighed the effects of ecological parameters on parasite risk in proboscis monkeys. Ecological effects tested in the present study (group size, the proportion of juvenile group members per group, group type and seasonal climatic changes) did not explain differences in parasite risk among sampling sites.

No clinical symptoms (such as diarrhea or poor body condition) could be observed in proboscis monkeys inhabiting the LBPMS. However, presented results revealed that an anthropogenically managed area can influence the parasite risk in proboscis monkeys with potentially pathogenic species. Although infestation with Trichuris spp. and Oesophagostomum spp. usually remains subclinical, fatal whipworm infection with heavy worm burdens, secondary complications like bacterial infection, and stress have proven clinically important resulting in severe enteritis, anorexia or death of the host [45, 64]. Likewise, Trichuris spp. infection was described to decrease the population size of a wild nonhuman primate species [5]. Moreover, strongyloidiasis poses a serious health problem in nonhuman primate colonies [65] and was detected as a primary cause of death in rehabilitant orangutans (Pongo pygmaeus) [66]. Common symptoms are debilitation, reduced growth rates, dyspnea caused by migrating larvae, and death [45]. Therefore, appropriate strategies and prophylactic measures for proboscis monkeys and other anthropogenically managed wildlife populations infected with detected parasite species are required.

Important management tools to face helminth infections comprise a consistent, noninvasive endoparasite monitoring to identify pathogenic parasites together with a targeted antiparasitic treatment which can reduce parasite infestation and may support breeding success [67]. For different nonhuman primate hosts, effective anthelmintic regimes against detected species are available [45] and strongly suggested to improve welfare. If animals are fed on a daily basis like proboscis monkeys at the LBPMS, regular deworming would be feasible and advisable. However, successful eradication of Strongyloides and other soil-transmitted species like Trichuris and strongyles requires further strict sanitary practices. To interrupt the external development and survival of these soil-transmitted helminthes in enclosures, surfaces should be kept dry, feces needs to be removed daily from the feeding platforms and surroundings, and contamination of food and water needs to be impeded [45]. Likewise, periodic replacement of soil or rotation of platform use is recommended to prevent from recurrent reinfections [7]. Additionally, appropriate group sizes and overall host densities should be considered to limit stress that is known to increase susceptibility to parasite infections. Reference values for appropriate group sizes of proboscis monkeys and other colobines in captivity are available in specialist guidelines like the expert report on the minimum requirements for the keeping of mammals published by the German Federal Ministry of Food and Agriculture (http://www.bmel.de/EN/Animals/AnimalWelfare.html, retrieved 3rd February 2017).

Conclusions

In conclusion, our study supported the hypothesis that parasite risk is elevated in anthropogenically modified habitats [68]. The site of origin outweighed the effects of ecological parameters on parasite risk in the proboscis monkey. Especially through increased host density in public sanctuaries and on feeding platforms, endangered proboscis monkeys are at risk of parasitic infections, particularly during the rainy season. Whereas climatic and temporal dynamics are evident in both study sites, group structure and composition had no significant effect for transmission dynamics and disease risk with intestinal parasites. Thus, management strategies for this large endangered primate should focus on regular parasite assessment, consistent deworming procedures and sanitary practices. All in all, this study has provided strong empirical evidence that noninvasive fecal collection represents an important tool for conservation medicine to rapidly assess the intestinal parasite risk of endangered wildlife in tropical environments.

Supporting information

S1 Table. Best models (GLMM) for prevalence of nematode infections and influencing variables among helminth groups found in proboscis monkeys, calculated with binomial presence/absence data.

If sampling site had a significant effect, two separate models–one for each sampling site—were calculated per helminth group.

https://doi.org/10.1371/journal.pone.0195584.s001

(DOCX)

S2 Table. Best models (GLMM) for parasite species richness (PSR) of nematode infections and influencing variables among helminth groups found in proboscis monkeys, calculated under the assumption of a Poisson distribution.

If sampling site had a significant effect, two separate models—one for each sampling site—were calculated.

https://doi.org/10.1371/journal.pone.0195584.s002

(DOCX)

S3 Table. Best linear mixed effect models for log(1+x)-transformed numbers of eggs per gram feces (epg) of nematode infections and influencing variables among helminth groups found in proboscis monkeys.

If sampling site had a significant effect, two separate models for sampling sites were calculated per helminth group.

https://doi.org/10.1371/journal.pone.0195584.s003

(DOCX)

Acknowledgments

The authors gratefully acknowledge permission by the director of Labuk Bay Proboscis Monkey Sanctuary (LBPMS), Mr. Michael Lee Hing Huat, as well as Mr. Sean Lee, to perform observations and data collection. Likewise we would like to thank the Sabah Wildlife Department and the Sabah Biodiversity Centre for giving permission to carry out field work in the Lower Kinabatangan Wildlife Sanctuary. Furthermore, we thank all staff of LBPMS and the field assistants of Danau Girang Field Centre for their support in field work as well as Sönke von den Berg and Tjard Bergmann for great aid with artwork.

References

  1. 1. Stuart MD, Strier KB. Primates and parasites—a case for a multidisciplinary approach. Int J Primatol. 1995; 16(4): 577–93.
  2. 2. Pedersen AB, Jones KE, Nunn CL, Altizer S. Infectious diseases and extinction risk in wild mammals. Conserv Biol. 2007; 21(5): 1269–79. pmid:17883492
  3. 3. Gillespie TR, Chapman CA. Prediction of parasite infection dynamics in primate metapopulations based on attributes of forest fragmentation. Conserv Biol. 2006; 20(2): 441–8. pmid:16903105
  4. 4. Nguyen N, Fashing PJ, Boyd DA, Barry TS, Burke RJ, Goodale CB, et al. Fitness impacts of tapeworm parasitism on wild gelada monkeys at Guassa, Ethiopia. Am J Primatol. 2015; 77(5): 579–94. pmid:25716944
  5. 5. Chapman CA, Gillespie TR, Speirs ML. Parasite prevalence and richness in sympatric colobines: effects of host density. Am J Primatol. 2005; 67(2): 259–66. pmid:16229007
  6. 6. Gillespie TR, Chapman CA. Forest fragmentation, the decline of an endangered primate, and changes in host-parasite interactions relative to an unfragmented forest. Am J Primatol. 2008; 70(3): 222–30. pmid:17879941
  7. 7. Becker DJ, Streicker DG, Altizer S. Linking anthropogenic resources to wildlife–pathogen dynamics: a review and meta-analysis. Ecol lett. 2015; 18(5): 483–95. pmid:25808224
  8. 8. Wenz-Mücke A, Sithithaworn P, Petney TN, Taraschewski H. Human contact influences the foraging behaviour and parasite community in long-tailed macaques. Parasitology. 2013; 140(06): 709–18.
  9. 9. Mbora DN, McPeek MA. Host density and human activities mediate increased parasite prevalence and richness in primates threatened by habitat loss and fragmentation. J Anim Ecol. 2009; 78(1): 210–8. pmid:19120603
  10. 10. Clough D, Heistermann M, Kappeler PM. Host intrinsic determinants and potential consequences of parasite infection in free-ranging red-fronted lemurs (Eulemur fulvus rufus). Am J Phys Anthropol. 2010; 142(3): 441–52. pmid:20091843
  11. 11. Chapman CA, Rothman JM, Hodder SAM. Can parasite infections be a selective force influencing primate group size? A test with red colobus. In: Huffman MA, Chapman CA, editors. Primate parasite ecology: the dynamics and study of host-parasite relationships. Cambridge: Cambridge University Press; 2009. p. 423–40.
  12. 12. Snaith TV, Chapman CA, Rothman JM, Wasserman MD. Bigger groups have fewer parasites and similar cortisol levels: a multi-group analysis in red colobus monkeys. Am J Primatol. 2008; 70(11): 1072–80. pmid:18666135
  13. 13. MacIntosh AJ, Hernandez AD, Huffman MA. Host age, sex, and reproductive seasonality affect nematode parasitism in wild Japanese macaques. Primates; journal of primatology. 2010; 51(4): 353–64. pmid:20711744
  14. 14. Klaus A, Zimmermann E, Röper KM, Radespiel U, Nathan S, Goossens B, et al. Co-infection patterns of intestinal parasites in arboreal primates (proboscis monkeys, Nasalis larvatus) in Borneo. Int J Parasitol Parasites Wildl. 2017; 6(3): 320–9.
  15. 15. Zuk M, McKean KA. Sex differences in parasite infections: patterns and processes. Int J Parasitol. 1996; 26(10): 1009–24. pmid:8982783
  16. 16. Eckert KA, Hahn NE, Genz A, Kitchen DM, Stuart MD, Averbeck GA, et al. Coprological surveys of Alouatta pigra at two sites in Belize. Int J Primatol. 2006; 27(1): 227–38.
  17. 17. Radespiel U, Schaber K, Kessler S, Schaarschmidt F, Strube C. Variations in the excretion patterns of helminth eggs in two sympatric mouse lemur species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar. Parasitol Res. 2015; 114(3): 941–54. pmid:25563604
  18. 18. Sodhi NS, Posa MRC, Lee TM, Bickford D, Koh LP, Brook BW. The state and conservation of Southeast Asian biodiversity. Biodivers Conserv. 2010; 19(2): 317–28.
  19. 19. Goossens B, Ambu LN. Sabah Wildlife Department and 10 years of research: towards a better conservation of Sabah’s wildlife. J Oil Palm, Environ and Health (JOPEH). 2012; 3: 38–51.
  20. 20. Röper KM, Scheumann M, Wiechert AB, Nathan S, Goossens B, Owren MJ, et al. Vocal acoustics in the endangered proboscis monkey (Nasalis larvatus). Am J Primatol. 2014; 76(2): 192–201. pmid:24123122
  21. 21. Meijaard E, Nijman V. Distribution and conservation of the proboscis monkey (Nasalis larvatus) in Kalimantan, Indonesia. Biol Conserv. 2000; 92(1): 15–24.
  22. 22. Sha JC, Bernard H, Nathan S. Status and conservation of proboscis monkeys (Nasalis larvatus) in Sabah, East Malaysia. Primat Conserv. 2008; 23(1): 107–20.
  23. 23. Boonratana R. Fragmentation and its significance on the conservation of proboscis monkey (Nasalis larvatus) in the Lower Kinabatangan, Sabah (North Borneo). In: Marsh LK, Chapman CA, editors. Primates in Fragments Complexity and Resilience. Heidelberg, Dordrecht, London: Springer New York; 2013. p. 459–75.
  24. 24. King N, Nair V, Mohamed B, Bahauddin A, editors. Gap analysis of future studies on wildlife tourism in Lower Kinabatangan, Sabah. Proceedings of the International Conference on Tourism Development; 2013 Feb 4–5; Penang, Malaysia. Penang, Malaysia: Sustainable Tourism Research Cluster, Universiti Sains Malaysia; 2013.
  25. 25. Matsuda I, Tuuga A, Akiyama Y, Higashi S. Selection of river crossing location and sleeping site by proboscis monkeys (Nasalis larvatus) in Sabah, Malaysia. Am J Primatol. 2008; 70(11): 1097–101. pmid:18651612
  26. 26. LBPMS. Labuk Bay Proboscis Monkey Sanctuary [Internet]. Sandakan, Sabah, Malaysia, Mile 19, Jalan Labuk (Junction SPS 3): Labuk Bay Proboscis Monkey Sanctuary; c2013 [cited 2016 Dec 20]. http://www.proboscis.cc.
  27. 27. Azmi R. Natural vegetation of the Kinabatangan floodplain. Part 1. An introduction to its natural vegetation, including a preliminary plant checklist of the region: Report WWF Malaysia; 1998.
  28. 28. Stark D, Nijman V, Lhota S, Robins J, Goossens B. Modeling population viability of local proboscis monkey Nasalis larvatus populations: conservation implications. Endanger Species Res. 2012; 16(1): 31–43.
  29. 29. Ancrenaz M, Goossens B, Gimenez O, Sawang A, Lackman-Ancrenaz I. Determination of ape distribution and population size using ground and aerial surveys: a case study with orang-utans in lower Kinabatangan, Sabah, Malaysia. Anim Conserv. 2004; 7: 375–85.
  30. 30. Bidin K, Chappell NA. Characteristics of rain events at an inland locality in northeastern Borneo, Malaysia. Hydrol Process. 2006; 20(18): 3835–50.
  31. 31. Semire FA, Mohd-Mokhtar R, Ismail W, Mohamad N, Mandeep J. Ground validation of space-borne satellite rainfall products in Malaysia. Adv Space Res. 2012; 50(9): 1241–9.
  32. 32. Gahalaut VK. Borneo’s proboscis monkey—a study of its diet of mineral and phytochemical concentrations. Curr Sci India. 2005; 89(3): 454.
  33. 33. Matsuda I, Tuuga A, Higashi S. Effects of water level on sleeping-site selection and inter-group association in proboscis monkeys: why do they sleep alone inland on flooded days? Ecol Res. 2010; 25(2): 475–82.
  34. 34. Yeager CP. Feeding ecology of the proboscis monkey (Nasalis larvatus). Int J Primatol. 1989; 10(6): 497–530.
  35. 35. Stark DJ, Vaughan IP, Saldivar DAR, Nathan SK, Goossens B. Evaluating methods for estimating home ranges using GPS collars: A comparison using proboscis monkeys (Nasalis larvatus). Plos One. 2017; 12(3): e0174891. pmid:28362872
  36. 36. Gillespie TR. Noninvasive assessment of gastrointestinal parasite infections in free-ranging primates. Int J Primatol. 2006; 27(4): 1129–43.
  37. 37. Bates D, Maechler M., Bolker B., Walker S. Fitting Linear Mixed-Effects Models using lme4. J Stat Softw. 2015; 67(1): 1–48.
  38. 38. Fox J. Effect displays in R for Generalised Linear Models. J Stat Softw. 2003; 8(15): 1–27.
  39. 39. Ganesh S, Cruz RJ Jr. Strongyloidiasis: a multifaceted disease. Gastroenterol Hepatol (NY). 2011; 7(3): 194–6.
  40. 40. Bethony J, Brooker S, Albonico M, Geiger SM, Loukas A, Diemert D, et al. Soil-transmitted helminth infections: ascariasis, trichuriasis, and hookworm. Lancet. 2006; 367(9521): 1521–32. pmid:16679166
  41. 41. Cogswell F. Parasites of non-human primates. In: Baker DG, editor. Flynn’s parasites of laboratory animals. 2 ed. Iowa: Blackell Publishing; 2007. p. 693–743.
  42. 42. Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biom J. 2008; 50(3): 346–63. pmid:18481363
  43. 43. Jessee MT, Schilling PW, Stunkard J. Identification of intestinal helminth eggs in old world primates. Lab Anim Care. 1970; 20(1): 83–7. pmid:4244728
  44. 44. Ghai RR, Chapman CA, Omeja PA, Davies TJ, Goldberg TL. Nodule worm infection in humans and wild primates in Uganda: cryptic species in a newly identified region of human transmission. PLoS Negl Trop Dis. 2014; 8(1): e2641. pmid:24421915
  45. 45. Strait K, Else JG, Eberhard ML. Parasitic diseases of nonhuman primates. In: Abee CR, Mansfield K, Tardif SD, Morris T, editors. Nonhuman primates in biomedical research: Diseases. 2nd ed. London: Academic Press; 2012. p. 197–298.
  46. 46. Nunn CL, Altizer S, Jones KE, Sechrest W. Comparative tests of parasite species richness in primates. Am Nat. 2003; 162(5): 597–614. pmid:14618538
  47. 47. Lawson B, Robinson RA, Colvile KM, Peck KM, Chantrey J, Pennycott TW, et al. The emergence and spread of finch trichomonosis in the British Isles. Philos Trans R Soc Lond, B, Biol Sci. 2012; 367(1604): 2852–63. pmid:22966140
  48. 48. Dunn J. The foraging strategy of Alouatta palliata: Variations in diet and feeding behaviour associated with food availability (habitat size and seasonality) and their effects on physiological health: PhD Dissertation, University of Barcelona, Barcelona, Spain; 2009.
  49. 49. Martínez-Mota R, Valdespino C, Sánchez-Ramos M, Serio-Silva J. Effects of forest fragmentation on the physiological stress response of black howler monkeys. Anim Conserv. 2007; 10(3): 374–9.
  50. 50. Padgett DA, Glaser R. How stress influences the immune response. Trends Immunol. 2003; 24(8): 444–8. pmid:12909458
  51. 51. Navarro-Gonzalez N, Fernández-Llario P, Pérez-Martín JE, Mentaberre G, López-Martín JM, Lavín S, et al. Supplemental feeding drives endoparasite infection in wild boar in Western Spain. Vet parasitol. 2013; 196(1): 114–23.
  52. 52. Albenberg LG, Wu GD. Diet and the intestinal microbiome: associations, functions, and implications for health and disease. Gastroenterology. 2014; 146(6): 1564–72. pmid:24503132
  53. 53. Zaiss M, Harris N. Interactions between the intestinal microbiome and helminth parasites. Parasite Immunol. 2016; 38(1): 5–11. pmid:26345715
  54. 54. Freeland WJ. Primate social groups as biological islands. Ecology. 1979; 60(4): 719–28.
  55. 55. Nunn C, Altizer SM. Infectious diseases in primates: behavior, ecology and evolution. Oxford: Oxford University Press; 2006.
  56. 56. Matsuda I, Kubo T, Tuuga A, Higashi S. A Bayesian analysis of the temporal change of local density of proboscis monkeys: implications for environmental effects on a multilevel society. Am J Phys Anthropol. 2010; 142(2): 235–45. pmid:20091847
  57. 57. Müller-Graf C, Collins D, Woolhouse M. Intestinal parasite burden in five troops of olive baboons (Papio cynocephalus anubis) in Gombe Stream National Park, Tanzania. Parasitology. 1996; 112(05): 489–97.
  58. 58. Gotoh S. Regional differences in the infection of wild Japanese macaques by gastrointestinal helminth parasites. Primates; journal of primatology. 2000; 41(3): 291–8.
  59. 59. Benavides JA, Huchard E, Pettorelli N, King AJ, Brown ME, Archer CE, et al. From parasite encounter to infection: multiple-scale drivers of parasite richness in a wild social primate population. Am J Phys Anthropol. 2012; 147(1): 52–63. pmid:21989507
  60. 60. Fuxjager MJ, Foufopoulos J, Diaz-Uriarte R, Marler CA. Functionally opposing effects of testosterone on two different types of parasite: implications for the immunocompetence handicap hypothesis. Funct Ecol. 2011; 25(1): 132–8.
  61. 61. Klein S. Hormonal and immunological mechanisms mediating sex differences in parasite infection. Parasite Immunol. 2004; 26(6–7): 247–64. pmid:15541029
  62. 62. Gillespie TR, Lonsdorf EV, Canfield EP, Meyer DJ, Nadler Y, Raphael J, et al. Demographic and ecological effects on patterns of parasitism in eastern chimpanzees (Pan troglodytes schweinfurthii) in Gombe National Park, Tanzania. Am J Phys Anthropol. 2010; 143(4): 534–44. pmid:20623606
  63. 63. Turner WC, Versfeld WD, Kilian JW, Getz WM. Synergistic effects of seasonal rainfall, parasites and demography on fluctuations in springbok body condition. J Anim Ecol. 2012; 81(1): 58–69. pmid:21831195
  64. 64. Crestian J, Crespeau F. Oesophagostomum stephanostomum in the chimpanzee. Receuil de medicine Veterinaire de l’Ecole d’Alfort. 1975; 151: 13–8.
  65. 65. Brack M, Göltenboth R, Rietschel W. Primaten. In: Göltenboth R, Klös H-G, editors. Krankheiten der Zoo-und Wildtiere. Berlin: Blackwell Wissenschafts-Verlag; 1995. p. 25–66.
  66. 66. Warren KS. Orang-utan conservation: Epidemiological aspects of health management and population genetics [PhD thesis]: Murdoch University; 2001.
  67. 67. Zerbe P, Schwarz C, Deb A, Borjal R, Arif A, Taha A, et al. Long-term monitoring of endoparasites in birds-of-paradise at Al Wabra Wildlife Preservation, Doha. Int Zoo Yearbook. 2009; 43(1): 149–58.
  68. 68. Chapman CA, Gillespie TR, Goldberg TL. Primates and the ecology of their infectious diseases: How will anthropogenic change affect host-parasite interactions? Evol Anthropol. 2005; 14(4): 134–44.