Social structure and Escherichia coli sharing in a group-living wild primate, Verreaux’s sifaka
© Springer et al. 2016
Received: 14 July 2015
Accepted: 28 January 2016
Published: 12 February 2016
Epidemiological models often use information on host social contacts to predict the potential impact of infectious diseases on host populations and the efficiency of control measures. It can be difficult, however, to determine whether social contacts are actually meaningful predictors of transmission. We investigated the role of host social structure in the transmission of Escherichia coli in a wild population of primates, Verreaux’s sifakas (Propithecus verreauxi). Using multilocus sequence typing (MLST), we compared genetic similarities between E. coli isolates from different individuals and groups to infer transmission pathways.
Correlation of social and transmission networks revealed that membership to the same group significantly predicted sharing of E. coli MLST sequence types (ST). Intergroup encounter rate and a measure of space-use sharing provided equally potent explanations for type sharing between social groups when closely related STs were taken into account, whereas animal age, sex and dispersal history had no influence. No antibiotic resistance was found, suggesting low rates of E. coli spillover from humans into this arboreal species.
We show that patterns of E. coli transmission reflect the social structure of this group-living lemur species. We discuss our results in the light of the species’ ecology and propose scent-marking, a type of social contact not considered in previous epidemiological studies, as a likely route of transmission between groups. However, further studies are needed to explicitly test this hypothesis and to further elucidate the relative roles of direct contact and environmental transmission in pathogen transfer.
KeywordsFecal-oral disease transmission Propithecus verreauxi Escherichia coli Multilocus sequence typing Scent-marking Social contacts Social network analysis Spillover
Pathogens and parasites—broadly defined as organisms that live on and draw nutrients from a host —may impact their hosts’ survival and fitness considerably [1–3], often leading to declines in wildlife populations [e.g. 4–6]. Epidemiological models predict the potential impact of infectious diseases on host populations and the efficiency of control measures. However, whereas most traditional models have assumed homogeneous mixing of individuals within a host population , it is now recognized that heterogeneity in contact patterns arising from both the social behavior of individual hosts and the spatial structure of host populations affect pathogen transmission [8–11]. Epidemiological models accounting for social contact heterogeneities in wildlife have been introduced in the 2000s [e.g. 12–14] and may yield predictions that differ dramatically from mean-field models assuming homogeneous mixing [15, 16], especially with regard to threshold population sizes for disease invasion  as well as transmission  and mortality rates .
Therefore, increasing efforts are being made to accurately assess transmission-relevant social contact patterns in humans  and wildlife  through direct observations , live-trapping  or the use of proximity loggers [20, 24]. It can be difficult, however, to establish whether the contact patterns that have been measured are meaningful for actual transmission . Moreover, it remains a major goal in disease ecology to determine the relative importance of social pathogen transmission compared to uptake from environmental reservoirs [25–27].
Escherichia coli is an ideal model organism to trace transmission of fecal-orally transmitted microorganisms . It is the main facultative anaerobic commensal colonizing the gastrointestinal tract of mammals and exhibits a clonal population structure that is little affected by genetic transfer and mutation . Healthy humans and dogs usually carry one predominant, resident strain of E. coli, which is present for months to years, and one to several transient strains [30–32]. Thus, if two individuals share the same or genetically very similar strains, it can be assumed that either direct transmission has occurred between them, or that they have been exposed to a common source of E. coli .
Among wild mammals, concordance between contact networks and E. coli strain sharing networks has been demonstrated in the fission–fusion society of giraffes , in a multihost system of African ungulates  and in solitary mountain brushtail possums . A similar study in elephants has been hampered by the fact that E. coli survive in water sources, which can act as infection reservoirs, masking social transmission . None of these studies however, examined the relationship between social networks and bacterial transmission across multiple discrete, adjacent social groups, which is a common setting in epidemiological models [35–38].
We studied social determinants of E. coli transmission in 10 groups of Verreaux’s sifakas (Propithecus verreauxi), an endemic Malagasy primate. Verreaux’s sifakas are diurnal, folivorous lemurs that are strictly arboreal. They live in small stable, multi-male multi-female groups and defend core areas of their home range by scent-marking in overlap areas of regular intergroup encounters [39, 40]. Males disperse from their natal group at around 3–5 years of age . In contrast to most other primates, the small group and home range size of sifakas permit simultaneous study of multiple neighboring groups [e.g. 39, 40], and thus a rare opportunity to study transmission within and between social groups of primates in a natural setting.
Microorganisms can be taken up directly from the environment or through contact with members of other species, potentially masking social transmission between conspecifics. However, because Verreaux’s sifakas do not drink from waterholes, but rely on licking dew from trees as well as on the water content of their diet, E. coli transmission via environmental reservoirs is unlikely. Furthermore, environmental conditions in the habitat are considered unfavorable for E. coli survival during most of the year [42, 43], as a long dry season is accompanied by high-amplitude daily temperature fluctuations, high levels of light exposure and low humidity .
We investigated whether E. coli transmission, measured as multilocus sequence type (MLST ST) sharing, is influenced by group membership, GPS-derived intergroup encounter rates and a measure of space-use sharing between groups, the Utilization Distribution Overlap Index (UDOI), while controlling for host sex, age and dispersal history. If direct social contacts were important for E. coli transmission, we predicted that ST sharing would be more prevalent within than between groups, and we expected intergroup encounter rates to better explain strain-sharing than joint space-use.
While E. coli transmission from environmental reservoirs may be unlikely, contact with humans represents a potential source of pathogens, as the study population has been exposed to visits by several thousand tourists per year. It has been shown that other lemurs living in habitats frequented by tourists harbor pathogenic Enterobacteria, which do not occur in undisturbed populations . E. coli serves as an important indicator of potential microbial spillover from humans or livestock into wildlife populations [45, 46]. Although E. coli is most often a commensal, many pathogenic strains exist, causing severe intestinal and extra-intestinal disease [47, 48], also in lemurs . Among these, multi-resistant, extended-spectrum beta-lactamase (ESBL)-producing strains have recently emerged worldwide as important causes of disease both in humans and animals [50, 51]. By testing for antibiotic resistance of E. coli isolates, we assess the risk of microbial spillover from humans to a species belonging to the most endangered group of mammals, the lemurs of Madagascar .
Study area and host population
The study was carried out in Kirindy Forest, Western Madagascar, located at 44°39′E, 20°03′S. The 90-ha study area is part of a field site operated by the German Primate Center, where lemurs are habituated and individually marked with unique collars, including radio (Holohil Systems, Carp, Ontario, Canada) and GPS (e-obs, Grünwald, Germany) units, as part of an ongoing long-term study , which has been approved by the Ethics Committee of the German Primate Center and the Ministère des Eaux et Forêts of Madagascar. All necessary research permits were obtained from the respective Malagasy and German authorities (Ministère des Eaux et Forêts of Madagascar; Commission ad hoc Flore et Faune (CAFF) of Madagascar; Centre National de Formation, d’Etudes et de Recherche en Environnement et Foresterie (CNFEREF); The Federal Agency for Nature Conservation of Germany). Regarding animal welfare, we followed the “Code of Best Practices for Field Primatology 2014” of the International Primatological Society.
We studied 10 social groups of Verreaux’s sifakas, 8 of which were adjacent. Samples from two groups living about 2 km away from the principal study population were also included. Groups ranged in size from 3 to 7 individuals, comprising 36–38 animals in total. Censuses of group membership are carried out two to three times a week .
Estimation of home range size, home range overlap and intergroup encounter rates
One adult animal in each of the 8 adjacent groups was equipped with a GPS collar. Scan sampling of group cohesion during animal observations showed that all other group members were in a radius of less than 10 m from collared individuals in more than 50 % of the time (range 53–90 %, mean 69.4 %), thus we considered ranging data from those individuals to be representative for the group as a whole. Animals considered likely to roam (i.e. young males) were not chosen to be collared. Collars were set to simultaneously record GPS coordinates every 15 min, from 04:00 to 20:00 h local time. GPS data were collected from August to December 2013 and from March to July 2014. GPS data were not available for the two groups living outside of the principal study area.
We calculated 95 % kernel home ranges and their overlaps for bi-weekly intervals using the adehabitatHR package  in R version 3.0.2. To quantify space-use sharing between the different groups, i.e. how much the animals actually use the overlap area, we calculated the Utilization Distribution Overlap Index (UDOI) . This index usually ranges from 0 to 1 but can be >1 if overlap is high and space-use is non-uniformly distributed. To derive intergroup encounter rates from the GPS data, we used the linear movement model contained in the R package movement Analysis , assuming linear movement between location measurements. An encounter was inferred if two groups were in ≤42 m distance based on the interpolated trajectories. The 42 m distance threshold was derived by calculating the distance between the groups’ GPS locations during directly observed intergroup encounters, based on an extended dataset of observations and GPS data collection over the course of one year . Using this threshold best predicted observed encounters. A new encounter was recorded if the two GPS-bearing individuals from different groups were at a distance >42 m for at least 30 min, until this threshold was crossed again.
Encounter rates were calculated as encounters per day for bi-weekly intervals. We tested for correlation between UDOI and intergroup encounter rate using Spearman rank correlation.
In order to collect data on animal behavior during encounters, direct observations were conducted in the 8 adjacent groups during two periods, once during the dry, lean season from August to October 2013 and once during the wet season from February to May 2014. Focal animal observations of 1 h per individual were carried out in an alternating order for 3 h in the morning and 3 h in the afternoon, resulting in 4 statistical days per individual, which amounted to a total of 860 h of focal animal observations. During intergroup encounters, the identities and proportion of participating animals, all close contacts (i.e. grooming, body contact or proximity of <1 m) between members of different groups and their durations as well as the total duration of the encounter were recorded ad libitum.
Sample collection, E. coli isolation and genotyping
Rectal swabs (Transwab® Amies, Medical Wire and Equipment, Corsham, Wiltshire, UK) were taken in the course of routine biomedical examinations of immobilized animals. Fecal samples were collected within 2 min of defecation from those individuals from whom rectal swabs could not be obtained. Samples were obtained during 3 periods, each spanning at most 25 days: In March and April 2013 and in August 2013, we took rectal swabs from 14 and 11 individuals, respectively, and in March and April 2014, 24 rectal swabs and 15 fecal samples were taken, resulting in a total of 66 samples from 48 individuals. Individuals from all groups were sampled in random order in each period (Additional file 1), to exclude the possibility that MLST sharing within groups could arise as an artifact of the time point of sampling. Four individuals were repeatedly sampled. Pre-cultivation of bacteria was undertaken in the field laboratory to ensure that a sufficient number would survive storage and transport and thus maximize E. coli recovery: Rectal swabs were streaked within 48 h onto MacConkey and Columbia blood agar and used to inoculate glucose-containing nutrient broth (agar and broth: Oxoid GmbH, Wesel, Germany). Broth and agar plates were incubated for 24–32 h at 37 °C. After this first incubation period, broth was streaked onto both MacConkey and Columbia blood agar and incubated for another 24 h. To maximize recovery of E. coli, we randomly picked several colonies from all four agar plates, representing all colony morphology types present on the plate, dissolved them in sterile 0.9 % sodium chloride solution with an addition of 20 % glycerol and froze them at −20 °C until shipment and further processing. Fecal samples were treated slightly differently: They were used to inoculate glucose-containing nutrient broth, and after incubation of 24–32 h, an aliquot of the broth was frozen at −20 °C with an addition of 20 % glycerol.
In the laboratory in Germany, samples were streaked out onto MacConkey and Columbia blood agar. After an incubation period of 24–48 h, up to 4 colonies typical for E. coli were isolated per sample and subjected to species identification using matrix-assisted laser desorption ionization time-of-flight mass-spectrometry (MALDI-ToF MS; Bruker GmbH, Bremen, Germany). Antimicrobial susceptibility testing was done by agar disc diffusion as recommended by the European Committee on Antimicrobial Susceptibility Testing (EUCAST) applying EUCAST clinical breakpoints for categorization of susceptible, intermediate and resistant isolates. All isolates were tested for multidrug resistance due to the production of extended-spectrum beta-lactamases (ESBL) using a chromogenic agar plate (chromID™ESBL; Bio Mérieux, Marcy l’Etoile, France). Isolates belonging to ST131 were additionally tested for resistance against the following antibiotics: Ampicillin, Piperacillin, Cefuroxim, Cefotaxim, Cefpodoxim, Ceftazidim, Cefepime, Piperacin/Tazobactam, Imipenem, Meropenem, Ertapenem, Trimethoprim/Sulfamethoxazol, Tigecyclin, Gentamicin, Amikacin, Ciprofloxacin, Fosfomycin and Nitrofurantoin. These include, but are not limited to, the most prevalent resistance-causing antibiotics used in Madagascar .
For molecular subtyping, each isolate was characterized using MLST . This typing method relies on determination of the sequence of internal fragments of seven housekeeping genes  and STs were assigned according to the E. coli MLST website (http://mlst.warwick.ac.uk/mlst/dbs/Ecoli). Sequences were further analyzed using the SeqSphere+ software version 1 (Ridom GmbH, Münster, Germany). The minimum spanning tree based on the MLST was generated also using the SeqSphere+ software.
Construction of the transmission network
We constructed an E. coli transmission network by assigning a link to a dyad if the two animals shared the same ST. Links were unweighted, i.e. we did not make a distinction between the number of types that animals shared, and undirected as we did not make any assumptions about the direction of transmission. We also constructed a second network, taking strains that only differed in one of the seven loci (single locus variants, SLV) into account as evidence for less recent transmission. Animals from which we did not obtain any isolate were not included.
We tested for significant differences in E. coli recovery rate between the two sampling methods, rectal swab and fecal sample, using a χ2-test. We then tested whether the number of STs isolated from an individual was significantly correlated with the animals’ age, the size of its’ social group and sampling effort, i.e. the number of samples that were taken from this individual, using Spearman rank correlations. Because adult males are expected to roam more than adult females, we tested for differences in the number of sequence types between adult males and adult females using a Wilcoxon rank sum test. We used nonparametric tests, because the data violated assumptions of generalized linear mixed models.
To determine which factors influence the likelihood of a link occurring in the transmission network, we analysed pair-wise strain-sharing using Bayesian regression modelling. Because of the non-independence of network data, statistical methods that assume data independence are not appropriate . We controlled for non-independence by modelling the identities of animals as random effects with a multi-membership structure in generalized linear mixed models (logit link) using the R package MCMCglmm .
belonging to the same sex,
being born in the same year,
being born in the same social group,
being a member of the same social group at the point of sampling,
having been sampled in the same month and
having been sampled by the same method (rectal swab or fecal sample)
Including only the 8 neighboring groups, we further tested whether UDOIs and intergroup encounter rates were significantly correlated with E. coli type sharing between dyads belonging to different groups. For UDOIs as well as encounter rates, we used the arithmetic mean of all bi-weekly measurements to best represent the relation between the two groups throughout the study period, and we z-transformed those values to a mean of zero and a standard deviation of one for better comparability. We constructed separate models for these 2 predictors because of their strong correlation, and controlled for birth cohort (i.e. same age), sex, natal group, sampling month and sampling method. In a second set of models, we included sharing of single locus variants as evidence for less recent transmission.
We used non-informative priors and fixed the observation-level variance (R-structure) at 1. We ran the models for 1,000,000 iterations, with a thinning parameter of 500 and a burn in phase of 10,000 to achieve low autocorrelation between recorded iterations. This resulted in a sample size of 1980 for estimating the posterior distribution. All statistical analyses were performed in R v. 3.2.2.
Home range overlap and intergroup encounters
A total of 581 encounters were inferred based on GPS data, and mean bi-weekly intergroup encounter rates varied from 0.005 to 0.7 encounters per day. Fourty-nine intergroup encounters (8.4 % of all encounters) were directly observed. Mean duration (±SD) of observed intergroup encounters was 12.5 ± 11.5 min (range 2–60 min). Scent-marking occurred in 75.5 % of those encounters, while proximity <1 m between members of different groups occurred in 32.7 % and body contact and grooming only occurred during 16.3 and 12.2 % of encounters, respectively. Mean UDOIs and mean intergroup encounter rates of the 14 neighboring group-dyads were strongly correlated (Spearman rank correlation, N = 14, S = 12, rho = 0.97, P < 0.001).
Recovery of E. coli types
Predictors of individual ST richness and sharing
There was no significant correlation between the number of STs isolated from an individual and its age, the size of its social group, or the number of samples collected from an individual (Spearman rank correlation, rho = 0.05, P = 0.76; rho = −0.12, P = 0.42; rho = 0.25 and P = 0.08, respectively). Adult males did not harbor more STs than adult females (Wilcoxon rank sum test, W = 183.5, P = 0.54).
Bayesian general linear mixed models testing the influence of several predictor variables on the probability of sharing the same E. coli MLST
Posterior mean coefficient
95 % Confidence interval
MCMCglmm P value
Same birth year
Same natal group
Same sampling month
Same sampling method
Same birth year
Same natal group
Space-use sharing (UDOI)
Same sampling month
Same sampling method
Same birth year
Same natal group
Intergroup encounter rate
Same sampling month
Same sampling method
Bayesian general linear mixed models testing the influence of several predictor variables on the probability of sharing the same E. coli MLST or a single locus variant
Posterior mean coefficient
95 % Confidence interval
MCMCglmm P value
Same birth year
Same natal group
Same sampling month
Same sampling method
Same birth year
Same natal group
Space-use sharing (UDOI)
Same sampling month
Same sampling method
Same birth year
Same natal group
Intergroup encounter rate
Same sampling month
Same sampling method
Using E. coli as a model for other, potentially more harmful fecal-orally transmitted microorganisms, this study revealed that the social structure of a group-living primate shapes pathogen transmission. Group membership was the single best explanatory variable for E. coli strain sharing, while both intergroup encounter rates and space-use sharing qualified as predictors for inter-group transmission when single locus variants of E. coli STs were taken into account as evidence for less recent transmission. These measures of sociality should therefore be considered in epidemiological modeling of fecal-orally transmitted infectious agents. In contrast, there was no evidence for E. coli spillover from humans, despite exposure to anthropogenic activities, including human defecation, in the study area. The strictly arboreal lifestyle of sifakas, in combination with reliance on ephemeral water sources, may therefore hamper environmental uptake of E. coli and other similarly transmitted microorganisms, although comparison with a sympatric lemur species that spends more time on the ground (Eulemur rufifrons) will be needed to confirm this hypothesis [cf. 63].
Extensive E. coli strain sharing has been documented in members of the same household, including pets [64, 65], suggesting that households function as a microbiological unit. Our results show that wild groups mirror these microbiological units in Verreaux’s sifakas. Belonging to the same social group was the only significant predictor of E. coli ST sharing as compared to possible non-social influences such as host sex and age, which were important determinants of E. coli community composition in elephants . Being born in the same group also had no significant effect, indicating that adult animals acquire the E. coli types of their new group after dispersal, and that their E. coli community is not static throughout life, but responds to changing social conditions, as in humans . Belonging to the same birth cohort even exhibited a negative correlation with E. coli type sharing. This makes sense in light of the fact that in most years, one or two juveniles are born into each social group. Thus, each birth cohort consists of members of different groups, which have a higher likelihood of harboring the same E. coli type as their group members rather than that of conspecifics of the same age.
We recovered a maximum of 4 different STs from an individual and this number adequately represented within-host E. coli diversity found in humans  and other animals [26, 31]. One limitation of the study, however, is that we were unable to assess strain-turnover in individual hosts, as only 4 animals could be sampled more than once due to practical constraints of studying wild primates. In those animals, we only once recovered the same ST twice. Future research will have to reveal whether this was due to strain-turnover, variation in strain abundance within the host, or the fact that sometimes subdominant strains can be more easily cultured than the dominant, resident strain . However, E. coli strains have been shown to persist for 1–3 years in humans and their pets .
Furthermore, if strain-turnover on the population level was indeed high, sampling during different periods should obscure any social signature rather than strengthen it, especially since sampling was equally distributed over all of the groups during each sampling period. To investigate whether our results were influenced by the different sampling periods or by the fact that we used both rectal swabs and fecal samples, we tested whether the probability of ST sharing was influenced by the animals being sampled in the same month or by the same method, but none of these terms was significant in any of the models. While there might be temporal turnover of strains within hosts, the same strains could circulate within social groups, as has been shown for mountain brushtail possums .
It is possible that animals within a group only shared the same strains because they live in the same habitat. However, transmission from water sources on the ground can be excluded because sifakas do not drink from water holes. This is in contrast to a wetland elephant population, where infection from this environmental reservoir masked social factors . Unfortunately, there is no empirical information on how long E. coli bacteria are able to survive on leaves and substrates used by sifakas, such as tree bark. However, light exposure in the dry deciduous Kirindy Forest is high, especially during the long dry season from April to October, and is accompanied by high-amplitude temperature fluctuations and little rainfall [41, 43], conditions which are regarded hostile for E. coli survival . Conditions may be more favorable during the short wet season, however, which is characterized by more stable temperatures and regular rainfall.
Regarding transmission between different social groups, we found no significant influences on sharing the same ST. However, if sharing of single locus variants was included, both intergroup encounter rates and space-use sharing were significantly correlated with the transmission network. Thus, the effect of social structure on inter-group transmission became only apparent when we included evidence for less recent transmission events, indicating that the frequency of transmission between groups may be considerably lower than within groups. Unfortunately, we were not able to test the influence of specific social behaviors on transmission, due to disappearance of animals during the study period and a correspondingly low sample size of behavioral observations.
Because space sharing is a prerequisite for encounters to occur, encounter rate was expected and found to strongly correlate with our measure of range overlap, which takes into account the intensity of space use in the overlap area (UDOI). Consequently, we evaluated these parameters in two different Bayesian GLMMS. Both models had very similar deviance information criterion (DIC) values, indicating that neither model should be preferred over the other. Furthermore, the standardized effect sizes of intergroup encounter rates and UDOI regarding E. coli transmission were almost identical, making it difficult to determine whether transmission takes place predominantly through direct contact or through the subsequent use of the same substrates. We can thus not rule out environmental transmission. In contrast, association of individuals correlated with E. coli type sharing in giraffes, whereas home range overlap did not , but the intensity of use of the overlap area was not controlled for.
Most observed sifaka intergroup encounters did not involve physical contact. We recorded members of different groups to be in a distance of <1 m in only 33 % of observed encounters, and body contact occurred in only 16 % of them. Alternatively, we propose that transmission might be mediated by scent-marking. Sifakas scent-mark by rubbing their anogenital region and, in the case of males, chest glands on trees, and subsequent olfactory inspection of scent-marks and overmarking, i.e. placing a scent-mark directly on the scent-mark of another individual, are frequent behaviors in both male and female sifakas  and other lemurs . Scent-marking has been shown to occur at higher rates in zones of home range overlap than in core areas  and occurred during 75.5 % of observed encounters. Recently, it has been shown that meerkats belonging to the same social group harbor similar bacterial communities in their anal scent glands, indicating that sociality plays a role in the acquisition of these bacterial assemblages . Moreover, olfactory inspection of conspecific cues has been suggested as an explanation of Salmonella transmission in sleepy lizards . Thus, olfactory communication involving signal inspection and overmarking, which are also common in other mammals [73, 74], and potentially other types of indirect species-specific social contacts are worth considering as mechanisms of social pathogen transmission in future studies of many other mammals.
Spillover of antimicrobial-resistant bacteria from humans into wildlife populations, which do not themselves come into contact with antimicrobial agents, has become an increasing concern, as these wild animal populations might constitute reservoirs for human infections [51, 75]. ESBL-producing E. coli have been found in wildlife populations ranging from seagulls and birds of prey to rodents and wild boar [reviewed in 51]. More recently, multidrug-resistant E. coli have been identified at high prevalence in banded mongoose living close to human settlements  and in a non-habituated gorilla living in a protected area .
Most E. coli isolates recovered in this study have not been described before in humans or domestic animals and were non-ESBL-producing. Thus, while we certainly did not sample the entire E. coli diversity harbored by sifakas, the data suggest that spillover of E. coli from humans into this population might be low, even though the study area is frequented regularly by tourists, researchers and occasionally villagers and livestock. We did not investigate E. coli carriage in humans in this study, but a survey including human stool samples from all over Madagascar revealed that more than 80 % of E. coli isolates were resistant to the most widely used antibiotics, including those tested in our study . Bublitz et al.  recently found a higher prevalence of pathogenic Enterobacteria in lemurs living in habitats with frequent exposure to humans than in undisturbed forests. Lemurs found to be positive did not include Propithecus sp., however.
Nonetheless, we did recover ST131 isolates in 3 individuals. The spread of a single clonal group with an ST131 profile has largely been responsible for the sudden increase in ESBL-producing E. coli during the last decades and today constitutes the predominant lineage among extraintestinal pathogenic E. coli worldwide [77, 78]. Varying prevalences in healthy human populations suggest that ST131 may be a dominant strain also in healthy individuals . ESBL-producing bacteria, including E. coli ST131, occur with high prevalence in the region of Antananarivo, the capital of Madagascar [79, 80]. The ST131 isolates which we recovered in this study were non-ESBL-producing, however. Additional susceptibility testing on all isolates belonging to this ST revealed full susceptibility to all antimicrobial agents tested. Along with some studies on humans, our results confirm that non-ESBL-producing, antibiotic-susceptible ST131 isolates exist [81, 82]. Another type known for ESBL-production, ST101, was also isolated, but was also non-ESBL-producing.
In conclusion, our results provide suggestive evidence for socially-structured transmission of potentially pathogenic microorganisms within and between wild groups of primates. We showed that E. coli strain-sharing is most prevalent within groups of wild lemurs, while intergroup relationships affect population-wide spread. We propose that species-specific patterns of scent-marking and overmarking might constitute a likely route of social transmission, especially between groups. Further studies are needed to explicitly test this hypothesis and to further elucidate the relative roles of direct contact and environmental transmission. The absence of antibiotic resistance indicates that E. coli spillover into the population might be low, despite relatively high human impact, although spillover of other pathogens cannot be excluded. More detailed genetic analyses will be needed to clarify the extent to which the ST131 isolates discovered in this study are related to disease-causing human ST131 isolates.
deviance information criterion
- E. coli :
European Committee on Antimicrobial Susceptibility Testing
general linear mixed model
global positioning system
- MALDI-ToF MS:
matrix-assisted laser desorption ionization time-of-flight mass-spectrometry
multi-locus sequence typing
utilization distribution overlap index
AS, AM, CF and PMK designed research. AS conducted the field work, performed the statistical analyses and drafted the manuscript. AM performed E. coli isolation, genotyping and analyses. PMK and CF analysed data and contributed to manuscript writing. All authors read and approved the final manuscript.
We thank the team of the Kirindy field station and the Malagasy Ministère de l’Environnement et des Eaux et Forêts, the Département Biologie Animale de l’Université d’Antananarivo, and the Centre National de Formation, d’Etudes et de Recherche en Environnement et Foresterie for supporting and authorizing our long-term research in Kirindy. We are grateful to H. Sennhenn-Reulen for statistical advice and to C. L. Nunn and F. Leendertz as well as three anonymous reviewers for valuable comments on an earlier version of this manuscript. We thank all members of the research group FOR 2136 “Sociality and health in primates” for fruitful discussions. This research was funded by the Deutsche Forschungsgemeinschaft (Ka 1082/29-1) and the German Primate Center. The publication of this article was funded by the Open Access fund of the Leibniz Association.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Nunn CL, Altizer S. Infectious diseases in primates: behaviour, ecology and evolution. Oxford: Oxford University Press; 2006.View ArticleGoogle Scholar
- Hillegass MA, Waterman JM, Roth JD. Parasite removal increases reproductive success in a social African ground squirrel. Behav Ecol. 2010;21:696–700.View ArticleGoogle Scholar
- Milton K. Effects of bot fly (Alouattamyia baeri) parasitism on a free-ranging howler monkey (Alouatta palliata) population in Panama. J Zool. 1996;239:39–63.View ArticleGoogle Scholar
- Bermejo M, Rodríguez-Teijeiro JD, Illera G, Barroso A, Vilà C, Walsh PD. Ebola outbreak killed 5000 gorillas. Science. 2006;314:1564.View ArticlePubMedGoogle Scholar
- Packer C, Altizer S, Appel M, Brown E, Martenson J, O’Brien SJ, et al. Viruses of the Serengeti: patterns of infection and mortality in African lions. J Anim Ecol. 1999;68:1161–78.View ArticleGoogle Scholar
- LaDeau SL, Kilpatrick AM, Marra PP. West Nile virus emergence and large-scale declines of North American bird populations. Nature. 2007;447:710–3.View ArticlePubMedGoogle Scholar
- McCallum H, Barlow N, Hone J. How should pathogen transmission be modelled? Trends Ecol Evol. 2001;16:295–300.View ArticlePubMedGoogle Scholar
- Griffin R, Nunn C. Community structure and the spread of infectious disease in primate social networks. Evol Ecol. 2012;26:779–800.View ArticleGoogle Scholar
- Craft ME, Volz E, Packer C, Meyers LA. Disease transmission in territorial populations: the small-world network of Serengeti lions. J R Soc Interface. 2011;8:776–86.PubMed CentralView ArticlePubMedGoogle Scholar
- Ezenwa VO. Host social behavior and parasitic infection: a multifactorial approach. Behav Ecol. 2004;15:446–54.View ArticleGoogle Scholar
- Altizer S, Nunn CL, Thrall PH, Gittleman JL, Antonovics J, Cunningham AA, et al. Social organization and parasite risk in mammals: integrating theory and empirical studies. Annu Rev Ecol Evol Syst. 2003;34:517–47.View ArticleGoogle Scholar
- Cross P, Lloyd-Smith J, Bowers J, Hay C, Hofmeyr M, Getz W. Integrating association data and disease dynamics in a social ungulate: bovine tuberculosis in African buffalo in the Kruger National Park. Ann Zool Fenn. 2004;41:879–92.Google Scholar
- Keeling MJ, Eames KTD. Networks and epidemic models. J R Soc Interface. 2005;2:295–307.PubMed CentralView ArticlePubMedGoogle Scholar
- Perkins SE, Ferrari MF, Hudson PJ. The effects of social structure and sex-biased transmission on macroparasite infection. Parasitology. 2008;135:1561–9.View ArticlePubMedGoogle Scholar
- Manlove KR, Cassirer EF, Cross PC, Plowright RK, Hudson PJ. Costs and benefits of group living with disease: a case study of pneumonia in bighorn lambs (Ovis canadensis). Proc R Soc B Biol Sci. 2014;. doi:10.1098/rspb.2014.2331.Google Scholar
- Carne C, Semple S, Morrogh-Bernard H, Zuberbühler K, Lehmann J. The risk of disease to great apes: simulating disease spread in orang-utan (Pongo pygmaeus wurmbii) and chimpanzee (Pan troglodytes schweinfurthii) association networks. PLoS ONE. 2014;. doi:10.1371/journal.pone.0095039.Google Scholar
- Lloyd-Smith JO, Cross PC, Briggs CJ, Daugherty M, Getz WM, Latto J, et al. Should we expect population thresholds for wildlife disease? Trends Ecol Evol. 2005;20:511–9.View ArticlePubMedGoogle Scholar
- Molina C, Stone L. Modelling the spread of diseases in clustered networks. J Theor Biol. 2012;315:110–8.View ArticlePubMedGoogle Scholar
- Hamede R, Bashford J, Jones M, McCallum H. Simulating devil facial tumour disease outbreaks across empirically derived contact networks. J Appl Ecol. 2012;49:447–56.View ArticleGoogle Scholar
- Obadia T, Silhol R, Opatowski L, Temime L, Legrand J, Thiébaut ACM, et al. Detailed contact data and the dissemination of Staphylococcus aureus in hospitals. PLoS Comput Biol. 2015;. doi:10.1371/journal.pcbi.1004170.PubMed CentralPubMedGoogle Scholar
- Craft ME. Infectious disease transmission and contact networks in wildlife and livestock. Philos Trans R Soc B Biol Sci. 2015;. doi:10.1098/rstb.2014.0107.Google Scholar
- MacIntosh AJJ, Jacobs A, Garcia C, Shimizu K, Mouri K, Huffman MA, et al. Monkeys in the middle: parasite transmission through the social network of a wild primate. PLoS ONE. 2012;. doi:10.1371/journal.pone.0051144.Google Scholar
- VanderWaal K, Atwill E, Hooper S, Buckle K, McCowan B. Network structure and prevalence of Cryptosporidium in Belding’s ground squirrels. Behav Ecol Sociobiol. 2013;67:1951–9.View ArticleGoogle Scholar
- Vander Wal E, Laforge MP, McLoughlin PD. Density dependence in social behaviour: home range overlap and density interacts to affect conspecific encounter rates in a gregarious ungulate. Behav Ecol Sociobiol. 2013;68:383–90.View ArticleGoogle Scholar
- Chiyo PI, Grieneisen LE, Wittemyer G, Moss CJ, Lee PC, Douglas-Hamilton I, et al. The influence of social structure, habitat, and host traits on the transmission of Escherichia coli in wild elephants. PLoS ONE. 2014;. doi:10.1371/journal.pone.0093408.Google Scholar
- VanderWaal KL, Atwill ER, Isbell LA, McCowan B. Linking social and pathogen transmission networks using microbial genetics in giraffe (Giraffa camelopardalis). J Anim Ecol. 2013;83:406–14.View ArticlePubMedGoogle Scholar
- Kappeler PM, Cremer S, Nunn CL. Sociality and health: impacts of sociality on disease susceptibility and transmission in animal and human societies. Philos Trans R Soc Lond B Biol Sci. 2015;. doi:10.1098/rstb.2014.0116.PubMed CentralGoogle Scholar
- Archie EA, Luikart G, Ezenwa VO. Infecting epidemiology with genetics: a new frontier in disease ecology. Trends Ecol Evol. 2009;24:21–30.View ArticlePubMedGoogle Scholar
- Tenaillon O, Skurnik D, Picard B, Denamur E. The population genetics of commensal Escherichia coli. Nat Rev Microbiol. 2010;8:207–17.View ArticlePubMedGoogle Scholar
- Sears HJ, Brownlee I, Uchiyama JK. Persistence of individual strains of Escherichia coli in the intestinal tract of man. J Bacteriol. 1950;59:293–301.PubMed CentralPubMedGoogle Scholar
- Sears HJ, Janes H, Saloum R, Brownlee I, Lamoreaux LF. Persistence of individual strains of Escherichia coli in man and dog under varying conditions. J Bacteriol. 1956;71:370–2.PubMed CentralPubMedGoogle Scholar
- Caugant DA, Levin BR, Selander RK. Genetic diversity and temporal variation in the E. coli population of a human host. Genetics. 1981;98:467–90.PubMed CentralPubMedGoogle Scholar
- VanderWaal KL, Atwill ER, Isbell LA, McCowan B. Quantifying microbe transmission networks for wild and domestic ungulates in Kenya. Biol Conserv. 2014;169:136–46.View ArticleGoogle Scholar
- Blyton MDJ, Banks SC, Peakall R, Lindenmayer DB, Gordon DM. Not all types of host contacts are equal when it comes to E. coli transmission. Ecol Lett. 2014;17:970–8.View ArticlePubMedGoogle Scholar
- Nunn CL, Thrall PH, Leendertz FH, Boesch C. The spread of fecally transmitted parasites in socially-structured populations. PLoS ONE. 2011;. doi:10.1371/journal.pone.0021677.Google Scholar
- Bonnell TR, Sengupta RR, Chapman CA, Goldberg TL. An agent-based model of red colobus resources and disease dynamics implicates key resource sites as hot spots of disease transmission. Ecol Modell. 2010;221:2491–500.View ArticleGoogle Scholar
- Caillaud D, Craft ME, Meyers LA. Epidemiological effects of group size variation in social species. J R Soc Interface. 2013;. doi:10.1098/rsif.2013.0206.PubMed CentralPubMedGoogle Scholar
- Leslie E, Cowled B, Garner MG, Toribio J, Ward MP. Effective surveillance strategies following a potential classical swine fever incursion in a remote wild pig population in north-western Australia. Transbound Emerg Dis. 2014;61:432–42.View ArticlePubMedGoogle Scholar
- Richard AF. Social boundaries in a Malagasy prosimian, the sifaka (Propithecus verreauxi). Int J Primatol. 1985;6:553–68.View ArticleGoogle Scholar
- Benadi G, Fichtel C, Kappeler P. Intergroup relations and home range use in Verreaux’s sifaka (Propithecus verreauxi). Am J Primatol. 2008;70:956–65.View ArticlePubMedGoogle Scholar
- Kappeler PM, Fichtel C. A 15-year perspective on the social organization and life history of sifaka in Kirindy Forest. In: Kappeler PM, Watts DP, editors. Long-term field studies of primates. Berlin: Springer; 2012. p. 101–21.View ArticleGoogle Scholar
- van Elsas JD, Semenov AV, Costa R, Trevors JT. Survival of Escherichia coli in the environment: fundamental and public health aspects. ISME J. 2011;5:173–83.PubMed CentralView ArticlePubMedGoogle Scholar
- Ottoson JR, Nyberg K, Lindqvist R, Albihn A. Quantitative microbial risk assessment for Escherichia coli O157 on lettuce, based on survival data from controlled studies in a climate chamber. J Food Protect. 2011;74:2000–7.View ArticleGoogle Scholar
- Bublitz DC, Wright PC, Rasambainarivo FT, Arrigo-Nelson SJ, Bodager JR, Gillespie TR. Pathogenic enterobacteria in lemurs associated with anthropogenic disturbance. Am J Primatol. 2014;77:330–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Pesapane R, Ponder M, Alexander KA. Tracking pathogen transmission at the human–wildlife interface: banded mongoose and Escherichia coli. EcoHealth. 2013;10:115–28.View ArticlePubMedGoogle Scholar
- Rwego IB, Isabirye-Basuta G, Gillespie TR, Goldberg TL. Gastrointestinal bacterial transmission among humans, mountain gorillas, and livestock in Bwindi Impenetrable National Park, Uganda. Conserv Biol. 2008;22:1600–7.View ArticlePubMedGoogle Scholar
- Nataro JP, Kaper JB. Diarrheagenic Escherichia coli. Clin Microbiol Rev. 1998;11:142–201.PubMed CentralPubMedGoogle Scholar
- Smith JL, Fratamico PM, Gunther NW. Extraintestinal pathogenic Escherichia coli. Foodborne Pathog Dis. 2007;4:134–63.View ArticlePubMedGoogle Scholar
- Williams CV. Chapter 36 - Prosimians. In: Fowler REME, editor. Fowler’s zoo and wild animal medicine, vol. 8. St. Louis: W.B. Saunders; 2015. p. 291–301.Google Scholar
- Oteo J, Perez-Vazquez M, Campos J. Extended-spectrum beta-lactamase producing Escherichia coli: changing epidemiology and clinical impact. Curr Opin Infect Dis. 2010;23:320–6.View ArticlePubMedGoogle Scholar
- Guenther S, Ewers C, Wieler LH. Extended-spectrum beta-lactamases producing E. coli in wildlife, yet another form of environmental pollution? Front Microbiol. 2011. doi:10.3389/fmicb.2011.00246.
- Schwitzer C, Mittermeier RA, Johnson SE, Donati G, Irwin M, Peacock H, et al. Averting lemur extinctions amid Madagascar’s political crisis. Science. 2014;343:842–3.View ArticlePubMedGoogle Scholar
- Calenge C. The package “adehabitat” for the R software: a tool for the analysis of space and habitat use by animals. Ecol Modell. 2006;197:516–9.View ArticleGoogle Scholar
- Fieberg J, Kochanny CO. Quantifying home-range overlap: the importance of the utilization distribution. J Wildl Manage. 2005;69:1346–59.View ArticleGoogle Scholar
- Sijben S. Movement analysis: analysis of trajectory data using linear or Brownian motion model. R package. 0.6 ed2013.
- Koch de Vasconcellos F. Intergroup relationships in Verreaux’s sifakas (Propithecus verreauxi). Göttingen: Georg-August-University; 2015.
- Randrianirina F, Ratsima EH, Ramparany L, Randremanana R, Rakotonirina HC, Andriamanantena T, et al. Antimicrobial resistance of bacterial enteropathogens isolated from stools in Madagascar. BMC Infect Dis. 2014;. doi:10.1186/1471-2334-14-104.Google Scholar
- Maiden MCJ, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, et al. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998;95:3140–5.PubMed CentralView ArticlePubMedGoogle Scholar
- Wirth T, Falush D, Lan R, Colles F, Mensa P, Wieler LH, et al. Sex and virulence in Escherichia coli: an evolutionary perspective. Mol Microbiol. 2006;60:1136–51.PubMed CentralView ArticlePubMedGoogle Scholar
- Croft DP, Madden JR, Franks DW, James R. Hypothesis testing in animal social networks. Trends Ecol Evol. 2011;26:502–7.View ArticlePubMedGoogle Scholar
- Hadfield JD. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J Stat Softw. 2010;33:1–22.View ArticleGoogle Scholar
- Nicolas-Chanoine MH, Blanco J, Leflon-Guibout V, Demarty R, Alonso MP, Canica MM, et al. Intercontinental emergence of Escherichia coli clone O25:H4-ST131 producing CTX-M-15. J Antimicrob Chemother. 2008;61:273–81.View ArticlePubMedGoogle Scholar
- 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:441–52.View ArticlePubMedGoogle Scholar
- Johnson JR, Miller S, Johnston B, Clabots C, DebRoy C. Sharing of Escherichia coli sequence type ST131 and other multidrug-resistant and urovirulent E. coli strains among dogs and cats within a household. J Clin Microbiol. 2009;47:3721–5.PubMed CentralView ArticlePubMedGoogle Scholar
- Johnson JR, Clabots C, Kuskowski MA. Multiple-host sharing, long-term persistence, and virulence of Escherichia coli clones from human and animal household members. J Clin Microbiol. 2008;46:4078–82.PubMed CentralView ArticlePubMedGoogle Scholar
- Skurnik D, Bonnet D, Bernède-Bauduin C, Michel R, Guette C, Becker J-M, et al. Characteristics of human intestinal Escherichia coli with changing environments. Environ Microbiol. 2008;10:2132–7.View ArticlePubMedGoogle Scholar
- de Muinck EJ, Øien T, Storrø O, Johnsen R, Stenseth NC, Rønningen KS, et al. Diversity, transmission and persistence of Escherichia coli in a cohort of mothers and their infants. Environ Microbiol Rep. 2011;3:352–9.View ArticlePubMedGoogle Scholar
- Blyton MDJ, Banks SC, Peakall R, Gordon DM. High temporal variability in commensal Escherichia coli strain communities of a herbivorous marsupial. Environ Microbiol. 2013;15:2162–72.View ArticlePubMedGoogle Scholar
- Lewis RJ. Sex differences in scent-marking in sifaka: mating conflict or male services? Am J Phys Anthropol. 2005;128:389–98.View ArticlePubMedGoogle Scholar
- Kappeler PM. To whom it may concern: the transmission and function of chemical signals in Lemur catta. Behav Ecol Sociobiol. 1998;42:411–21.View ArticleGoogle Scholar
- Leclaire S, Nielsen JF, Drea CM. Bacterial communities in meerkat anal scent secretions vary with host sex, age, and group membership. Behav Ecol. 2014;. doi:10.1093/beheco/aru074.Google Scholar
- Bull CM, Godfrey SS, Gordon DM. Social networks and the spread of Salmonella in a sleepy lizard population. Mol Ecol. 2012;21:4386–92.View ArticlePubMedGoogle Scholar
- Ferkin M. Scent over-marking and adjacent-marking as competitive tactics used during chemical communication in voles. In: Johnston R, Müller-Schwarze D, Sorensen P, editors. Advances in chemical signals in vertebrates. Springer: US; 1999. p. 239–46.View ArticleGoogle Scholar
- Jordan NR, Manser MB, Mwanguhya F, Kyabulima S, Rüedi P, Cant MA. Scent marking in wild banded mongooses: 1. Sex-specific scents and overmarking. Anim Behav. 2011;81:31–42.View ArticleGoogle Scholar
- Radhouani H, Silva N, Poeta P, Torres C, Correia S, Igrejas G. Potential impact of antimicrobial resistance in wildlife, environment and human health. Front Microbiol. 2014;. doi:10.3389/fmicb.2014.00023.PubMed CentralPubMedGoogle Scholar
- Janatova M, Albrechtova K, Petrzelkova KJ, Dolejska M, Papousek I, Masarikova M, et al. Antimicrobial-resistant Enterobacteriaceae from humans and wildlife in Dzanga–Sangha Protected Area, Central African Republic. Vet Microbiol. 2014;171:422–31.View ArticlePubMedGoogle Scholar
- Rogers BA, Sidjabat HE, Paterson DL. Escherichia coli O25b-ST131: a pandemic, multiresistant, community-associated strain. J Antimicrob Chemother. 2010;66:1–14.View ArticlePubMedGoogle Scholar
- Nicolas-Chanoine M-H, Bertrand X, Madec J-Y. Escherichia coli ST131, an intriguing clonal group. Clin Microbiol Rev. 2014;27:543–74.PubMed CentralView ArticlePubMedGoogle Scholar
- Andriatahina T, Randrianirina F, Hariniana ER, Talarmin A, Raobijaona H, Buisson Y, et al. High prevalence of fecal carriage of extended-spectrum β-lactamase-producing Escherichia coli and Klebsiella pneumoniae in a pediatric unit in Madagascar. BMC Infect Dis. 2010;. doi:10.1186/1471-2334-10-204.PubMed CentralPubMedGoogle Scholar
- Herindrainy P, Randrianirina F, Ratovoson R, Hariniana ER, Buisson Y, Genel N, et al. Rectal carriage of extended-spectrum beta-lactamase-producing gram-negative bacilli in community settings in Madagascar. PLoS ONE. 2011;. doi:10.1371/journal.pone.0022738.PubMed CentralPubMedGoogle Scholar
- Leflon-Guibout V, Blanco J, Amaqdouf K, Mora A, Guize L, Nicolas-Chanoine M-H. Absence of CTX-M enzymes but high prevalence of clones, including clone ST131, among fecal Escherichia coli isolates from healthy subjects living in the area of Paris, France. J Clin Microbiol. 2008;46:3900–5.PubMed CentralView ArticlePubMedGoogle Scholar
- Brisse S, Diancourt L, Laouénan C, Vigan M, Caro V, Arlet G, et al. Phylogenetic distribution of CTX-M- and non-extended-spectrum-Β-lactamase-producing Escherichia coli isolates: Group B2 isolates, except clone ST131, rarely produce CTX-M enzymes. J Clin Microbiol. 2012;50:2974–81.PubMed CentralView ArticlePubMedGoogle Scholar