Social structure and Escherichia coli sharing in a group-living wild primate, Verreaux’s sifaka
BMC Ecology volume 16, Article number: 6 (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.
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 .
We investigated in a multivariate model whether the following factors
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)
influenced a dyad’s log-odds of having the same E. coli ST.
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
Mean bi-weekly 95 % kernel home ranges varied in size from 5.67 to 14.15 ha among groups, and the corresponding pair-wise home range overlap of the 14 overlapping dyads varied between 0.64 and 3.88 ha (Additional file 2). Mean UDOIs of neighboring groups varied from 0.00017 to 0.098. Home ranges and overlaps from one exemplary bi-weekly period are shown in Fig. 1.
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
We characterized 83 E. coli isolates, belonging to 39 individuals from 10 social groups (mean 1.7 isolates per individual, range 0–6). Twenty-nine individuals were natal to their social group, while 10 animals (one female, nine males) were immigrants. Isolation success did not differ significantly between rectal swabs and fecal samples (76.5 % of Nrectal = 51 vs. 60 % of Nfecal = 15, χ2 = 0.78, P = 0.37). We isolated 24 distinct MLST STs, 13 of which (54.2 %) occurred in multiple hosts (Fig. 2). Previously unknown STs were deposited in the E. coli MLST database (http://mlst.warwick.ac.uk/mlst/dbs/Ecoli). Up to 2 different STs were obtained per sample and up to 4 different STs per individual host (mean: 1.25 STs). Of the 4 animals that were sampled repeatedly, we only once found the same type twice in subsequent samples. All isolates were non-ESBL-producing. Further susceptibility tests on isolates of type ST131 (N = 5), an ST that is shared with the major antibiotic resistant E. coli lineage , revealed full susceptibility to all antimicrobial agents tested.
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).
Density of the E. coli transmission network was 0.116 (Fig. 3b), indicating that 11.6 % of all possible connections were present. Dyads that belonged to the same social group were more likely to share the same or closely related STs as compared to dyads from different social groups with home range overlap (Fig. 4). We examined whether this pattern might be driven by mother-offspring ST-sharing, but none of the 10 mother-offspring dyads in the dataset shared the same ST. The likelihood of sharing STs was smallest for animals belonging to non-neighboring groups. We tested whether several dyad-level attributes influenced the likelihood of a link occurring between two individuals, but the full E. coli transmission network was only significantly correlated with the group membership network; all other predictor variables were not statistically significant (Table 1, model A). In fact, belonging to the same social group improved the odds of sharing the same ST 3.1 times. Regarding links between social groups, none of the variables tested were statistically significant predictors of sharing the same ST (Table 1, models B and C). However, in the second set of models which included sharing of closely related E. coli STs, UDOIs and intergroup encounter rates were strongly correlated with the intergroup transmission network (Table 2, models B and C), qualifying both as conduits of social transmission between different groups. Furthermore, the effect sizes of these parameters were almost identical. An increase of one standard deviation in home range overlap increased the odds of sharing the same or a closely related ST 1.54 times, whereas an increase of one standard deviation in intergroup encounter rates resulted in 1.58 times higher odds. Belonging to the same birth cohort was also significant, but the correlation coefficient was negative and the corresponding odds ratio <1 in both models, indicating that in this data subset, animals born in the same year had a smaller chance of harboring the same E. coli strain than by chance.
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
Nunn CL, Altizer S. Infectious diseases in primates: behaviour, ecology and evolution. Oxford: Oxford University Press; 2006.
Hillegass MA, Waterman JM, Roth JD. Parasite removal increases reproductive success in a social African ground squirrel. Behav Ecol. 2010;21:696–700.
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.
Bermejo M, Rodríguez-Teijeiro JD, Illera G, Barroso A, Vilà C, Walsh PD. Ebola outbreak killed 5000 gorillas. Science. 2006;314:1564.
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.
LaDeau SL, Kilpatrick AM, Marra PP. West Nile virus emergence and large-scale declines of North American bird populations. Nature. 2007;447:710–3.
McCallum H, Barlow N, Hone J. How should pathogen transmission be modelled? Trends Ecol Evol. 2001;16:295–300.
Griffin R, Nunn C. Community structure and the spread of infectious disease in primate social networks. Evol Ecol. 2012;26:779–800.
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.
Ezenwa VO. Host social behavior and parasitic infection: a multifactorial approach. Behav Ecol. 2004;15:446–54.
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.
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.
Keeling MJ, Eames KTD. Networks and epidemic models. J R Soc Interface. 2005;2:295–307.
Perkins SE, Ferrari MF, Hudson PJ. The effects of social structure and sex-biased transmission on macroparasite infection. Parasitology. 2008;135:1561–9.
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.
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.
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.
Molina C, Stone L. Modelling the spread of diseases in clustered networks. J Theor Biol. 2012;315:110–8.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Archie EA, Luikart G, Ezenwa VO. Infecting epidemiology with genetics: a new frontier in disease ecology. Trends Ecol Evol. 2009;24:21–30.
Tenaillon O, Skurnik D, Picard B, Denamur E. The population genetics of commensal Escherichia coli. Nat Rev Microbiol. 2010;8:207–17.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Richard AF. Social boundaries in a Malagasy prosimian, the sifaka (Propithecus verreauxi). Int J Primatol. 1985;6:553–68.
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.
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.
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.
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.
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.
Pesapane R, Ponder M, Alexander KA. Tracking pathogen transmission at the human–wildlife interface: banded mongoose and Escherichia coli. EcoHealth. 2013;10:115–28.
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.
Nataro JP, Kaper JB. Diarrheagenic Escherichia coli. Clin Microbiol Rev. 1998;11:142–201.
Smith JL, Fratamico PM, Gunther NW. Extraintestinal pathogenic Escherichia coli. Foodborne Pathog Dis. 2007;4:134–63.
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.
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.
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.
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.
Fieberg J, Kochanny CO. Quantifying home-range overlap: the importance of the utilization distribution. J Wildl Manage. 2005;69:1346–59.
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.
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.
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.
Croft DP, Madden JR, Franks DW, James R. Hypothesis testing in animal social networks. Trends Ecol Evol. 2011;26:502–7.
Hadfield JD. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J Stat Softw. 2010;33:1–22.
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.
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.
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.
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.
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.
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.
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.
Lewis RJ. Sex differences in scent-marking in sifaka: mating conflict or male services? Am J Phys Anthropol. 2005;128:389–98.
Kappeler PM. To whom it may concern: the transmission and function of chemical signals in Lemur catta. Behav Ecol Sociobiol. 1998;42:411–21.
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.
Bull CM, Godfrey SS, Gordon DM. Social networks and the spread of Salmonella in a sleepy lizard population. Mol Ecol. 2012;21:4386–92.
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.
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.
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.
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.
Rogers BA, Sidjabat HE, Paterson DL. Escherichia coli O25b-ST131: a pandemic, multiresistant, community-associated strain. J Antimicrob Chemother. 2010;66:1–14.
Nicolas-Chanoine M-H, Bertrand X, Madec J-Y. Escherichia coli ST131, an intriguing clonal group. Clin Microbiol Rev. 2014;27:543–74.
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.
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.
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.
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.
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.
Additional file 1. Isolates and host characteristics. A table including all of the E. coli isolates recovered in this study, together with information on their hosts, sampling date, and MLST profile.
Additional file 2. Dyadic home range overlap and intergroup encounter rates. A table of home range overlap and intergroup encounter rates per group-dyad.
About this article
Cite this article
Springer, A., Mellmann, A., Fichtel, C. et al. Social structure and Escherichia coli sharing in a group-living wild primate, Verreaux’s sifaka. BMC Ecol 16, 6 (2016). https://doi.org/10.1186/s12898-016-0059-y