Dispersal and grouping patterns form the foundations of social interactions
in group-living mammals and are the outcomes of a complex interplay between
inbreeding avoidance, kin cooperation and competition, predation pressure and
food resource distribution. In species where both sexes disperse, the
potential for kin-biased associations would seem limited. In one such
species, the western lowland gorilla (WLG), short-term data suggest that
female kin associations may be present due to directed local dispersal
decisions, but monitoring of groups over longer timescales is needed to
better elucidate this pattern. Using autosomal genotyping of 419 faecal
samples representing 85 unhabituated gorillas collected non-invasively over
5 years in a 132 km
Group living in many mammal species is generally explained as an arrangement
that maximizes the potential benefits of kin cooperation and predator
avoidance (Clutton-Brock and Lukas, 2012; Handley and Perrin, 2007). The
potential costs of inbreeding and competition, particularly with kin, are
usually invoked to explain the prevalence of sex-biased natal dispersal in
mammals (Handley and Perrin, 2007; Thierry, 2008). However, the dispersing
sex would seem to lose the potential for fitness-enhancing interactions with
life-long associates who may also be close kin (Langergraber et al., 2007,
2009; Silk, 2009; Silk et al., 2006). Unlike most mammal species which
display predominantly male dispersal (Greenwood, 1980), great apes such as
eastern gorillas, western lowland gorillas (WLGs,
It is currently thought that female gorillas gain little from residing with kin but obtain benefits from dispersing by avoiding potential inbreeding in their natal community and by having the ability to choose a new male in a new group (Breuer et al., 2012; Caillaud et al., 2008; Guschanski et al., 2008; Harcourt and Stewart, 2007). However, genetic analysis of multiple adjacent WLG groups at one site suggested that, despite routine dispersal by both sexes, female WLGs may maintain kin associations post-emigration (Bradley et al., 2007). Furthermore, in mountain gorillas, a WLG sister taxon displaying non-obligate female and male dispersal, non-dispersing related females are more tolerant of each other than of non-relatives (Watts, 1994a, b). On the other hand, Douadi et al. (2007) suggested that the proportion of related WLG female dyads (including both adults and immatures) within groups did not differ from the proportion between groups, while another study found few WLG female kin dyads in their study population (Inoue et al., 2013). However, the latter two studies relied upon estimation of dyadic relatedness values, which is not considered an accurate way to identify related pairs in a population (Csillery et al., 2006), thus leaving the question of female kin associations in WLGs unresolved.
Because female gorillas only transfer between groups during intergroup encounters or after group dissolution and are never solitary (Harcourt and Stewart, 2007), the timing and distance of single dispersal events are limited. Furthermore, because the risk of infanticide would appear to inhibit transfer by females with dependent young, opportunities for female dispersal are constrained in comparison to males (Robbins et al., 2013; Robbins, 2010). Nevertheless, females may move across multiple groups over a lifetime as the result of multiple transfers which are believed to reflect female choice (Breuer, 2010; Jeffery et al., 2007; Robbins et al., 2004). While females typically disperse singly, the transfer of pairs of females has been inferred twice, with two females observed in one group, and shortly thereafter in a new group after the death of their group's silverback (Parnell, 2002; Stokes et al., 2003). Such transfers subsequent to group dissolution are generally considered involuntary in nature. Based on these considerations and the observation of higher than expected average relatedness among females within groups at one site, it is hypothesized in WLGs that females may be able to maintain kin relationships despite natal and secondary dispersal (Bradley et al., 2007), which could in turn confer fitness benefits (Chapais, 2001; Clutton-Brock and Huchard, 2013; Hatchwell, 2010).
The dynamics of WLG group composition have been observed at several study sites (Bermejo, 2004; Cipolletta, 2004; Doran-Sheehy et al., 2004; Gatti et al., 2004; Jeffery et al., 2007; Magliocca et al., 1999; Parnell, 2002; Robbins et al., 2004; Stokes et al., 2003; Tutin, 1996), and these studies have elucidated some of the apparent contrasts between WLGs and mountain gorillas. For example, all-male bachelor groups, mixed-sex non-reproductive groups and single-male breeding groups have all been documented in WLGs, but in contrast to mountain gorillas, multi-male breeding groups have not been observed except for a few transient cases (Gatti et al., 2004; Jeffery et al., 2007; Magliocca et al., 1999; Parnell, 2002; Robbins et al., 2004; Stokes et al., 2003; Tutin, 1996). Furthermore, group fissions (the splitting of one group into two or more) have also not been observed in WLGs, and the only mode of group formation observed has been by the acquisition of females by lone silverback males (Gatti et al., 2004; Robbins et al., 2004). Group dissolutions (the cessation of existence of a group, wherein the group silverback no longer retains any females and females transfer to several different new groups) have been observed at one WLG site at a rate approximately 5 times higher than in mountain gorillas (Robbins et al., 2004). In general, however, the difficulty of simultaneous observational monitoring of multiple groups of a long-lived species has limited the direct study of patterns of dispersal and group formation in WLGs.
Genetic analyses provide an opportunity for understanding the dispersal
choices of individuals in such elusive species. In this study we use samples
collected within a 132 km
As detailed previously, we opportunistically collected a total of 396 gorilla
faecal samples across a 132 km
We stored samples for up to 1 year after collection and extracted DNA from
them using the QIAmp Stool kit (QIAGEN) with slight modifications (Nsubuga et
al., 2004). We determined the sex of the gorillas by amplifying a portion of
the amelogenin locus (Bradley et al., 2000). We also amplified the samples at
16 autosomal microsatellite loci using a two-step multiplex PCR method
described in detail elsewhere (Arandjelovic et al., 2009, 2010). We combined
up to four different PCR products and electrophoresed them on an ABI PRISM
3100 genetic analyser. We sized alleles relative to an internal size standard
(ROX-labelled HD400) using GeneMapper software version 3.7 (Applied
Biosystems). For accuracy, we validated heterozygous genotypes by observing
each allele in two or more independent reactions (Arandjelovic et al., 2009).
To ensure with
We used CERVUS 3.0 to compare genotypes to detect cases of different samples
originating from the same individual. For every case in which genotypes from
two samples had no mismatches and were identical at a minimum of seven loci, the
probability that two genotypes were exact matches despite originating from
two different individuals (
We evaluated minimum group membership as detailed in Arandjelovic et
al. (2010). In brief, we first assumed that when samples from two or more
individuals were collected on the same day at the same GPS location (same
nest site or multiple fresh faecal remains found within 20 m of each other),
these samples belonged to the same group. Once individuals were assigned to a
group, when samples from those individuals were then identified on other
sampling days, even if they were collected alone, then those individuals were
still considered part of the initial group assignment. Some identified groups
may thus actually be part of a larger group that was not well sampled over
time due to the heterogeneous and opportunistic sampling protocol used
(Arandjelovic et al., 2010). However, for the purposes of these analyses, for
lack of contrary evidence, and supported by the identification of at least
eight breeding groups at Loango by remote camera traps (Head et al., 2013), we
treat each identified group as its own social unit. When we detected an
individual in a particular group and then later on in the study period in
another group, we identified this as a dispersal event. We can only assume
that we have
To calculate the rate of dispersal and rate of group dissolutions, we estimated the number of recorded gorilla social unit years (sum of years each social unit was observed) by taking the total number of months between the initial and final sampling of any individual in any given gorilla social unit (group or lone silverback). Because we did not sample all individuals consistently over the study period and large time gaps occur between records of the same group, this is a very approximate measure of total detection time. With these caveats in mind, we have a total of 19 gorilla social unit years recorded.
We reconstructed partial pedigree relationships for individuals in the Loango gorilla study population using parentage analysis to identify possible offspring in the data set. We first did a simulation of parentage analysis in CERVUS 3.0 to determine the power of our microsatellite loci to assign mother–father–offspring trios with 85, 95 and 99 % confidence in our data set. We simulated genotypes for 10 000 offspring using the allele frequencies from our data set. We then ran a series of assignment simulations, assuming that we sampled from 30 % of mothers and fathers. Analyses were repeated assuming 20 and 10 % of parents were sampled to ensure results were robust (data not shown). The proportion of genotyped loci was 88 % and we conservatively allowed for a 1 % genotyping error rate.
In the parentage analysis, we used all males as potential fathers, all
females as potential mothers and all individuals as potential offspring and
allowed up to one mismatch for parent–offspring trios to account for
genotyping errors and mutational events. Firstly, we assigned paternity and
maternity when the trio-confidence score was
We evaluated dyadic relatedness (
First, we tested whether females in the study area are more related than
expected by chance. To do so, we estimated average
The inferred composition and dynamics of the Loango gorilla groups over the 5-year study period. Crosses designate that the gorilla was detected at least once in the corresponding month. Black arrows depict transfers between social units. Hatched lines represent the time at which the last nest group was found from the Mandondo or Achilles groups. LSB: lone silverback male; UKF: unknown female; UKM: unknown male. Unknown individuals were captured alone and once only and thus have no social unit affiliation. To calculate the number of recorded gorilla social unit years, we took the total number of months between the initial and final samplings of any individual in any given gorilla social unit (group or LSB) which corresponds to the shaded area for each group.
Second, we tested whether females occur in groups containing kin more often
than expected by chance (Bradley et al., 2007). In the Loango research area,
like at all WLG research sites, there is a high degree of home range overlap
for all of the study groups (Arandjelovic et al., 2010); thus even if females
are limited to transferring to neighbouring groups, we assume in this test
that all groups in the study area are equally possible transfer options. We
calculated the average dyadic relatedness of females within groups by
averaging across all pairwise comparisons and comparing this to the
distribution of average relatedness values obtained by permutation analysis
under a simplified model of random dispersal (Bradley et al., 2007; Lukas et
al., 2005; Manly, 1997). For the permutation, we pooled all females found in
groups from 2006 to 2007 (
Third, because we cannot be sure that all groups are equally likely transfer
options as assumed in the second analyses above, we tested whether females
occur in groups containing kin more often than expected by a scenario where
females preferentially disperse to neighbouring groups. To account for the
possibility that dispersal to closer groups is more likely than dispersal to
more distant groups, we permuted the females such that they had higher
likelihoods of being redistributed to geographically closer groups.
Specifically, we assigned groups with probabilities equaling
Using genetic analysis, we identified a total of 85 different gorillas during
2005–2009 in a 132 km
Individual gorillas are shown as circles (females) or rectangles
(males). Social units are
During the study, two of the groups apparently dissolved (Mandondo and Achilles) as we never resampled the individuals together nor were the groups directly observed after the suggested dissolution date (Figs. 1, 2a). Furthermore, when individuals from these groups were resampled they were either outside of their original range or found in association with new groups. Two groups formed over the course of the study (group J, a potential all-male group, and Atananga), both mostly comprising individuals previously found in association with other social units. Sample collection in 2009 was motivated by the suspected formation of a new group from direct observations. Even though only 23 samples were collected in 2009, these samples originated from five nest groups (allowing for repeated sampling, and therefore linking, of individuals found at the same site) and represent individuals from only two groups: Atananga (resampled from four nest group sites) and group B (sampled at one nest group site and containing three individuals detected from 2006 to 2007).
We confirmed the formation of the new Atananga group and deduced that its members consisted of a previously identified lone silverback male, four previously identified females from three different groups (two of which had dissolved) and two females of unknown origin. The two females that originated from the dissolved Achilles group were not a mother–daughter pair. We also inferred eight cases of female group transfer and three cases of male transfer between social units and found one dead adult male (G72) (Fig. 2a).
Samples from 2 of the 12 groups contained only males, and so these are possible all-male non-breeding groups (Figs. 1, 2a). Five groups contained both males and females, but we did not identify any pedigree relationships between these individuals using our criteria, and in the absence of age data we cannot be sure whether these are breeding or non-breeding groups (Levréro et al., 2006). From the pedigree analysis, we identified 15 parent–offspring trios in the remaining five groups (Fig. 2b), so we can be confident that this is the minimum number of breeding groups existing over the course of the study period. Although not set as a prior, parent–offspring trio assignments always identified candidate parents and offspring within the same group, except in two cases: (1) we found that a lone silverback (G19) is the probable father of a breeding female (G01) found in the Achilles group (with female G12 from the Mandondo group as the assigned mother), and (2) we identified the LayonA group silverback (G29) along with LayonA female G43 as the probable parents of group H female G44 (Fig. 2b).
In the first permutation analysis where average female relatedness was
compared to the permuted relatedness values of all adults across the study
area, we found that females were on average more related to each other
(
Specific examples of potential kin associations between females can also be extracted from the data. In the Indegho and Achilles groups, we identified one and two cases, respectively, where pairs of females shared an allele at all loci and thus represent potential mother–offspring, but mismatched the breeding male at two loci or more (Fig. 2b). Thus, these three females are not daughters of the silverback but either immigrated into groups which contained their kin; transferred into groups with their daughters; or, following group dissolution due to the death of a silverback, joined the same new silverback. Furthermore, putative half-sisters G06 and G01 (sharing Mandondo female G12 as their mother) both apparently immigrated into the Achilles group. In the Achilles group we also identify females G91, G97 and G01 as first-order relatives (potential mother–daughters or sisters) who were also unrelated to the group silverback.
By using genetic analysis to monitor multiple groups over several years, we indirectly observed the dynamics of group formation, group dissolution and individual movements in WLGs. As we have observed not only harem groups but also probable multi-male and mixed-sex non-reproductive groups at Loango, our results add to the growing body of evidence suggesting that the structure of WLG society contains various types of social units (Gatti et al., 2004; Levréro et al., 2006; Robbins et al., 2004). Furthermore, not all males in the reproductive groups are sons of the silverback, atypical for classical harem social units (Fig. 2b). Over a 5-year period, while genetically monitoring 18 social units, we inferred the dissolution of two groups, the formation of two groups, the death of one male and the dispersal of 13 individuals between social units. Taking into account the limits of our approximated number of recorded gorilla social unit years as detailed in the methods, the rate of group dissolution (2 in 19 social unit years: 11 %) is similar to that observed at Mbeli Bai (5 dissolutions in 63 gorilla group years: 8 %) (Robbins et al., 2004), as is the rate of female transfer: 8 cases in 19 social unit years (42 %) at Loango and 27 incidences in 63 group years (43 %) at Mbeli Bai (Stokes et al., 2003). It is important to note that these results represent the minimum number of dispersal events, social unit formations and dissolutions and that with increased sampling it is likely that more such events would be revealed. Similarly, this suggests that our rates of group dissolution and female transfer underestimate the true rates at Loango.
In many species, there are clear benefits for females to reside with female kin (Silk, 2009), so much so that most mammals are female philopatric (Greenwood, 1980; Pusey and Packer, 1987). Despite low levels of feeding competition and no evidence of between-group competition for resources, related mountain gorilla females exhibit more affiliative behaviours towards kin than non-kin when they live in the same group (Watts, 1994a) and female kin are more likely to form coalitions than non-kin in competition over food (Harcourt and Stewart, 1989; Watts, 2001). However, no studies have yet tested to see whether residing with and supporting kin improves their fitness (as shown in other species; meerkats: Clutton-Brock et al., 2010; white-faced capuchins: Fedigan et al., 2008; howler monkeys: Pope, 2000; savannah baboons: Van Horn et al., 2007).
As suggested in previous studies, our analyses imply that female WLGs do not disperse far from their natal range as the average relatedness of females across all groups at the study site was higher than expected by chance (Bradley et al., 2007; Douadi et al., 2007; Harcourt and Stewart, 2007; Jeffery et al., 2007). Furthermore, we found that, assuming that all groups at Loango are equally viable dispersal options, females co-occur in groups containing their kin more often than would be expected than if they dispersed randomly. When we attempted to correct for group proximity, we obtained mixed results depending on the relatedness estimator used. Therefore, we cannot distinguish between whether the patterns we observe are due to a process of related females dispersing to geographically proximate groups or whether there may be other factors that lead to related females occurring in the same post-dispersal groups. It is important to note, however, that due to limitations in the test procedure our geographic proximity permutation test is biased towards permuting co-residing females together. As we found that females within groups are more likely to be related than expected by chance, this means that our geographic proximity permutation test is biased towards finding related individuals occurring in the same group and therefore finding a non-significant result. This suggests that the significant result obtained with one of the two relatedness estimators may indeed be the more accurate of the two results and that females may not be basing their dispersal decisions on just geographic proximity. Female dispersal decisions are likely influenced by male quality (Breuer et al., 2010, 2012; Caillaud et al., 2008) and probably by habitat/vegetation composition familiarity as seen in mountain gorillas (Ganas et al., 2004; Guschanski et al., 2009). Our tests cannot determine the primary driver(s) of female dispersal, nor rule out that co-residency of female kin may be a secondary effect of other drivers of female dispersal decisions. Further long-term studies taking these factors into account will better elucidate how female relatives end up together post-dispersal.
Finding related post-dispersal females in the same group is consistent with a previous study on a smaller geographic scale (Bradley et al., 2007), but not with two others (Douadi et al., 2007; Inoue et al., 2013). The studies which failed to find a non-random distribution of female kin relied on assessment of dyadic kin relationships, an approach known to be less reliable than average group relatedness measures due to high levels of noise in the pairwise estimators (Csillery et al., 2006). Furthermore, as most samples came from a one-time sampling of unhabituated animals, some group members may have been omitted in the sampling, dampening the signal of within-group female relatedness. Finally, the different conclusions may be a result of smaller sample sizes in previous studies as it is well established that small sample sizes lead to reductions in statistical power. Inoue et al. (2013) included 12 adult females from 5 groups in their analyses and Douadi et al. (2007) included 14 females from 6 groups, while Bradley et al. (2007) analysed 22 adult females from 8 groups and we evaluated 35 adult females from 8 groups. We also inferred several instances of apparent natal and secondary transfer where females moved between groups either separately or together, which has also been observed in other populations (Robbins et al., 2004).
Bradley et al. (2007) suggest that post-dispersal female kin associations could occur in a scenario where female half-siblings (daughters of their group's silverback) reach dispersal age at similar times and transfer into the same group. We have evidence here that, in at least one case, the co-occurrence of half-siblings G01 and G06 in Achilles group occurred over a longer time period as G01 had already produced at least one offspring with the Achilles silverback when G06 immigrated into the group. Thus, females may not need to be in the same age cohort (and possibly may have never even co-resided in their natal group) but may still transfer into the same group over a longer time period; further long-term monitoring is needed to determine the kin knowledge of females prior to dispersal.
It has also been suggested that related females may co-reside in a new group if they co-transfer after the disintegration of their natal group (Bradley et al., 2007). In our data set we have one case of two Achilles group females (not a mother–daughter pair) who transferred to the Atananga group. We are unable, however, to specify the degree of relatedness of these females other than to say that they are not daughters of the Achilles silverback and thus that their transfer was probably a secondary dispersal event. Further sampling across the Loango study site to determine into which groups the Achilles females emigrated after the group dissolution should better elucidate the pattern of female kin associations in WLGs.
Our multi-year genetic monitoring study reinforces previous findings that WLGs exhibit dynamic group compositions and social structures as observed across their range. We also show that females, despite natal and secondary dispersal, appeared to retain associations with their same-sex kin post-emigration as females co-reside in groups containing female kin. The lack of long-term behavioural data on multiple western gorilla groups limits the inferences we can make on the function of these kin relations; however, dispersal need not be treated as an impediment to the possibility of kin-biased behaviours in the species.
M. Arandjelovic, L. Vigilant, M. M. Robbins and C. Boesch conceived and designed the study. J. Head, C. Boesch and M. M. Robbins collected samples, provided field information and field-based infrastructure. L. Vigilant contributed reagents and materials. M. Arandjelovic performed the experiments and analyses. M. Arandjelovic prepared the manuscript with contributions from all co-authors.
We thank the Agence Nationale des Parcs Nationaux (ANPN) and the Centre National de la Recherche Scientifique et Technique (CENAREST) of Gabon for permission to conduct our research in Loango National Park. We also thank L. Rabanal, L. Makaga, E. R. Guizard, E. Fairet, M. Gregoire, L. Rankin, E. Wright and the other field assistants of the Loango Ape project for their help collecting the samples. We are very grateful to D. Lukas, D. Caillaud, D. Morgan, C. Sanz, G. Schubert, O. Thalmann, B. J. Bradley, J. Ganas, K. Langergraber, T. Breuer, G. Campbell, J. Junker, A. Kalan, M. Whitten and C. Zacharias for helpful discussions, A. Abraham and C. Lang for laboratory assistance and R. Mundry and D. Lukas for statistical assistance and programming. We are also immensely grateful to the Editor and two anonymous reviewers for their constructive feedback and insightful suggestions which greatly improved this manuscript. This collaborative project is financially supported by, the Société pour la Conservation et le Développement (SCD), Wildlife Conservation Society (WCS), and the Max Planck Society. Edited by: J. Tung Reviewed by: two anonymous referees