- Open Access
Patterns of sex-biased dispersal are consistent with social and ecological constraints in a group-living cichlid fish
BMC Ecology and Evolution volume 22, Article number: 21 (2022)
Sex-biased dispersal is a common and widespread phenomenon that can fundamentally shape the genetic structure of the social environments in which animals live. For animals that live in and move between social groups, sex-biased dispersal can result in an asymmetry in the degree of relatedness among cohabiting males and females, which can have strong implications for their social evolution. In this study, we measured the relatedness structure within and across groups of a wild population of Neolamprologus multifasciatus, a highly-social, shell-dwelling cichlid fish endemic to Lake Tanganyika, East Africa. In total, we genotyped 812 fish from 128 social groups at 20 microsatellite loci. Neolamprologus multifasciatus live at high densities, and also experience strong ecological constraints on free movement throughout their habitat. At the same time, they exhibit sex differences in the degree of reproductive competition within their groups and this makes them an excellent model system for studying the factors associated with sex-biased dispersal.
Social groups of N. multifasciatus consist of multiple males and females living together. We found that cohabiting females were unrelated to one another (Lynch-Ritland estimates of relatedness = 0.045 ± 0.15, average ± SD), while males shared much higher, albeit variable, levels of relatedness to other males in their groups (0.23 ± 0.27). We uncovered a pronounced decline in relatedness between males living in separate groups as the spatial separation between them increased, a pattern that was not evident in females. Female dispersal was also markedly constrained by the distribution and availability of nearby territories to which they could emigrate.
Our results indicate female-biased dispersal in N. multifasciatus. Our study also highlights how the spatial distribution of suitable dispersal destinations can influence the movement decisions of animals. We also emphasize how sex-biased dispersal can influence the relatedness structure of the social environment in which individuals interact and compete with one another.
Dispersal is a ubiquitous life history trait that, through the mixing of individuals, influences the genetic structure of animal groups and populations, setting the social context in which selection operates [1, 2]. Dispersal decisions are thought to be part of a strategy that improves the breeding conditions of individuals, and can generally be categorized into natal dispersal, i.e., the movement of individuals from their birthplaces to the locations of their first breeding, and breeding dispersal, i.e., the movement of individuals between consecutive reproductive attempts [1, 3, 4]. For many organisms, dispersal is not uninhibited, but is rather influenced by their dispersal capacity and the availability of suitable patches on which to settle that are within reachable distances . For example, in Mauritius kestrels, Falco punctatus, females typically disperse within 2 km of their natal patch, but can be forced to disperse over much farther distances when local breeding sites are limiting . Also, in the spotted hyena, Crocuta crocuta, male dispersal patterns are dictated by the distribution of available clans to join and their relative qualities . In general, dispersal is thought to be variable and plastic, with decisions about when to disperse, how far to disperse, and where to settle differing among individuals and across different environmental circumstances [2, 5]. Therefore, in order to better understand the forces that shape the genetic structure of populations, groups, and their social environments, it is important to determine which individuals disperse, what factors determine their dispersal behaviour, and whether these factors affect certain individuals disproportionately to others.
Dispersal frequently differs between males and females of a species with one sex dispersing more readily and over greater distances than the other . For social animals that live in discrete groups, sex biases in the propensity for individuals to emigrate will shape the relatedness structure of their groups, and therefore also the social arena in which individuals interact and compete with one another. Sex-biased dispersal can cause one sex to be more closely related to their same-sex group mates than the other, resulting in asymmetric distributions of relatedness between interacting males and females and sex differences in selection for behavioral strategies. Sex-biased dispersal can therefore have far-reaching consequences for the evolution and expression of cooperative and social behaviours [9,10,11,12] and is fundamental to our understanding of social evolution .
Male-biased dispersal is particularly common in mammals (e.g., common voles, Microtus arvalis, ), while female-biased dispersal is typically more prevalent in birds (e.g., Seychelles warblers, Acrocephalus sechellensis, ). The current body of work on sex-biased dispersal suggests that selection can result in drastically different dispersal strategies for each sex [8, 16], and a myriad of ecological and social factors have been put forward in an attempt to explain these sex-specific patterns in nature. For example, sex biases in dispersal can arise when one sex experiences more intense local mate or resource competition [17,18,19], greater risk of inbreeding [19,20,21,22], or more pronounced constraints on emigration and settlement than the other [16, 23]. Many of the existing hypotheses for the evolution of sex-biased dispersal stem from a wealth of early work that focused on mammals and birds [8, 17, 24], but examinations of sex biases in dispersal have been increasing for other taxa as well. Fishes offer intriguing systems in which to study sex-biased dispersal because they display tremendous variation in mating systems, sociality, territorial behaviour, sex determination mechanisms, and parental care strategies [25,26,27,28], all factors that could affect sex-biased dispersal by influencing the costs and benefits of philopatric versus dispersive behaviours [8, 16]. Studies investigating sex-biased dispersal in fishes have so far produced examples of both male-biased (e.g., brook trout, Salvelinus fontinalis, ; Neolamprologus pulcher, ; three-spined sticklebacks, Gasterosteus aculeatus, ) and female-biased dispersal (e.g., Asian Seabass, Lates calcarifer, ). The variation that exists across species in terms of whether males or females are the more dispersing sex is particularly well-represented within fish family Cichlidae, which is also well-regarded for its extreme cross-species variation in social traits [33, 34]. Here, male-biased dispersal has been detected in rock- and sand-dwelling cichlids of Lake Malawi Pseudotropheus spp.  and Copadichromis spp.  as well as in N. pulcher from Lake Tanganyika . No evidence of sex biased dispersal has been found in the Central American species, Amphilophus astorquii , and female-biased dispersal has been uncovered in some Tanganyikan species, such as Eretmodus cyanostictus , N. caudopunctatus , and N. obscurus . It is interesting to note that dispersal biases in favour of one sex can be observed in cichlid species that also exhibit vastly different mating systems and breeding strategies. For example, E. cyanostictus is a biparental mouth brooder , N. caudopunctatus is a shoaling, but also biparental substrate spawner , and N. obscurus is a group-living, cooperative breeder , but all species display female-biased dispersal. The finding so far that neither male-biased nor female-biased dispersal is noticeably more prevalent within the cichlids stands in contrast to the more consistent patterns seen in other taxa, such as male-biased dispersal in mammals and female-biased dispersal in birds (see ), which makes them an intriguing system in which study sex-biased dispersal and its links to social evolution.
In this study, we tested for sex-biased dispersal as well as for environmental constraints on dispersal in a wild population of Neolamprologus multifasciatus, a group-living, highly social Lamprologine cichlid fish. Neolamprologus multifasciatus is among the smallest cichlid species endemic to Lake Tanganyika, East Africa, with males and females rarely exceeding 3.0 cm and 2.1 cm in standard length respectively in the wild [41, 42]. Neolamprologus multifasciatus can be found living in groups of up to ~ 20 individuals , and these groups consist of a single dominant male, several adult females and subordinate males, as well as immature, juvenile individuals [41,42,43,44]. Social groups of N. multifasciatus hold territories on vast regions of the lake floor called ‘shell beds’ where accumulations of empty gastropod shells cover the ground and the fish excavate these shells from the sandy substrata. The fish use these shells as shelters for evading predators and as brood chambers for females to raise their offspring. Over hundreds of hours of behavioural studies in the field, we have never observed clear instances of alloparental care in N. multifasciatus (AB, LK, AJ personal observations). Dispersal in N. multifasciatus appears to be delayed until individuals near or reach sexual maturity and adults can sometimes switch their group memberships between broods . Thus, it is possible that N. multifasciatus movement patterns are comprised of both natal and breeding dispersal. While males and females may move between groups over their lifetimes, large males can also establish new territories on their own when space is available to them [42, 44]. Neolamprologus multifasciatus territories can be extremely densely distributed across the shell bed, sometimes with only ~ 30 cm separating nearest neighbouring territories [42, 44]. Therefore, in certain areas of the shell bed there may be ample alternative groups nearby to serve as potential dispersal destinations, whereas in other areas of the shell bed such options can be scarcer. Furthermore, dispersal distances are typically very short—often no more than 200 cm—as gauged by the spatial separation between parents and their offspring when living apart in the wild . Such a limited scale of movement is likely due to a high risk of predation from piscivorous fishes that share the same habitat space .
Neolamprologus multfasciatus experience a number of ecological and social conditions that led us to predict dispersal would be sex-biased in this species, and that females might be the more-dispersing sex. To begin, we expected both sexes to be similarly constrained in their ability to disperse between territories on the shell bed as both males and females are susceptible to predation when moving away from their territories , though body size differences between the sexes might lead one sex to be more susceptible to predation than the other. Most reproduction is secured by the dominant male in their groups, which means that small, subordinate males achieve very little reproduction, both at home and in neighbouring territories . Subordinate males are also unable to establish independent territories of their own until they have grown large enough and sufficient space is available nearby . Thus, while subordinate males suffer low reproductive success at home, most of them also have few prospects for better success elsewhere. This might result in males remaining philopatric where they can queue for a breeding position. Relatedness between males can then promote cooperation and tolerance between cohabiting dominant and subordinate males. Unlike males, however, cohabiting females can reproduce concurrently , suggesting that young females have better prospects to reproduce soon after they emigrate than young males. Females should therefore be less likely to form queues in their natal groups, because the do have reproductive options elsewhere. Females are also highly competitive and will engage in agonistic interactions over limited resources within their groups, especially when offspring are being cared for . This suggests little scope for cooperation among females, and implies low or negligible inclusive fitness benefits for female relatives that reside near to one another. Lastly, while both sexes experience aggressive resistance from resident fish when attempting to join a new group, potential male joiners receive more aggression overall than potential female joiners , which suggests that group switching may be costlier for males than for females.
In this study, we evaluated the hypothesis that dispersal is female-biased in N. multifasciatus by using pairwise relatedness estimates from extensive microsatellite genotyping of a wild population along with high-resolution spatial data. In total, we sampled 812 fish from 128 social groups. We assessed the relatedness structure of social groups and also examined how relatedness among individuals changes with increasing spatial separation across the neighbourhood; in doing so, we tested for sex differences in dispersal probabilities and dispersal distances . Furthermore, if the dispersal potential of individuals is influenced by the distribution and availability of suitable territories to which to emigrate, then we also predicted to see higher within-group relatedness in groups from sparse regions of the shell bed relative to groups living in denser regions.
All 128 N. multifasciatus territories that we sampled from our study quadrat in the wild (Fig. 1) contained a dominant male, though we were unable to capture seven of them for measurement and fin clipping. Each group contained on average (± SD) 6.3 ± 4.4 fish, which comprised 1.9 ± 1.4 females (range = 0–6, Fig. 2A), 1.5 ± 0.9 males (range = 1–5, Fig. 2B), and 3.0 ± 3.0 juveniles (range = 0–14, Fig. 2C). The most common sex composition was a single dominant male with one adult female (Fig. 2D). Groups containing more males also contained more females (Spearman correlation coefficient = 0.38).
High within-group relatedness among males, but not among females
We calculated pairwise Lynch-Ritland estimators (rLR, ) to assess relatedness among the fish in our study. Females had low relatedness to other females in their groups, with an average (± SD) rLR of 0.045 ± 0.15 (Fig. 3). Males were more related to each other within their groups (average (± SD) rLR of 0.23 ± 0.27) than females were to each other (LMM, Est. ± SE = 0.16 ± 0.027, z = 6.08, P < 0.001). Adult males and females living in the same groups had average relatedness estimates of 0.092 ± 0.20 to each other, while juveniles were related to other juveniles in their groups at a level of 0.24 ± 0.24. Table 1 summarizes all Tukey contrasts between cohabiting fish types.
Relatedness among males between groups falls more steeply with spatial separation than for females
Female-female relatedness between groups declined with geographic separation (GAM, edf = 2.43, F = 8.49, P < 0.0001), though the slope was very shallow (Fig. 4). Male-male relatedness, on the other hand, declined more steeply with geographic separation than female-female relatedness (GAM, edf = 3.91, F = 36.6, P < 0.0001), and this sex difference was primarily seen within a two-meter radius around focal individuals (Fig. 4).
Female dispersal, but not male dispersal, is influenced by neighbourhood density
We calculated neighbourhood densities at multiple spatial scales by counting the number of other N. multifasciatus territories within 50 cm, 100 cm, 150 cm, and 200 cm of each focal group (Fig. 1). Then, using permutation tests, we found that within-group female-female relatedness was higher when there were few neighbouring territories, and was lower when there were many neighbouring territories. This pattern was evident at all spatial scales tested: 50 cm (P = 0.0019), 100 cm (P = 0.007), 150 cm (P = 0.018), and 200 cm (P = 0.043, Fig. 5A). In contrast, within-group male-male relatedness was unrelated to the number of neighbouring territories at all spatial scales: 50 cm (P = 0.51), 100 cm (P = 0.16), 150 cm (P = 0.25), and 200 cm (P = 0.30, Fig. 5B).
In this study, we uncovered evidence of female-biased dispersal in Neolamprologus multifasciatus, a highly social, group-living cichlid fish. Estimates of relatedness among cohabiting females were low, especially when compared to males, who shared a higher level of relatedness to other males in their groups (rLR = 0.23 ± 0.27, average ± SD, for males compared to 0.045 ± 0.15 for females). Furthermore, average relatedness coefficients among females from different territories were low, and this was independent of whether females resided in closely neighbouring or more distant territories. In contrast, males displayed elevated relatedness to their male neighbours relative to males living farther away. There was a distinct pattern of decreasing relatedness with geographic separation among males within our study quadrat, which was discernable over a spatial scale of approximately two meters. Taken together, these data indicate female-biased dispersal within our study system both in terms of probability for females to emigrate from their groups, as revealed by within-group relatedness patterns, and the distances over which females travel, as revealed by between-group relatedness patterns.
Our results indicate that the social environments within N. multifasciatus groups are characterized by a marked asymmetry in how related males and females are to their same-sex group mates. While males tended to be related to other males in their group, females were on average unrelated to other females (here, relatedness was computed with reference to the whole study quadrat). The question now, however, is whether such relatedness structure within groups coincides with patterns of sex-specific behaviour observed in our study species. Low average relatedness among females is consistent with the findings of previous studies that have documented high female-female conflict in this species [43, 45]; however, competition for resources and mates in N. multifasciatus is also thought to occur on extremely local spatial scales, which can lead to competition among individuals that is less affected, or even unaffected, by their relatedness . In N. multifasciatus, females divide the territory space into discrete sub-territories that they defend from each other, a behavioural pattern that is thought to be the outcome of intrasexual competition [45, 46]. Interactions among cohabiting males have been relatively unstudied to date, however, dominant males can clearly tolerate subordinate males to live on their territories. Given that dominant males are larger and more competitive than subordinates and that subordinates are often limited in their ability to establish their own independent breeding groups , this means that subordinates are ill-equipped to compete for reproduction at home or to move elsewhere, and may therefore be forced to make the best-of-a-bad-job by queuing to inherit an eventual breeding position. Tolerance of subordinate males forming such queues can be facilitated by genetic relatedness and/or by subordinates offering direct fitness benefits to the dominant male (i.e., pay-to-stay, [50, 51]). Indeed, we found that groups with additional males also supported more females (Fig. 2D), suggesting that subordinate males might provide direct fitness benefits to the dominant male by helping to maintain a larger and more productive territory. Field observations indicate that subordinate males, along with females, will contribute to territory defense, especially when intruders approach their sub-territories [41, 46]. We therefore suggest that while females disperse more often and move further afield, males will generally remain on their natal territories as subordinates and either inherit a future dominant breeding position, or move to a nearby vacancy when one becomes accessible.
Numerous proximate factors are commonly invoked to explain sex-biases in dispersal, and they include inbreeding avoidance (19–22; but see [16, 52]), asymmetric handicaps [8, 16, 23] and local mate or resource competition [17,18,19]. Female-biased dispersal in N. multifasciatus supports two of these non-mutually exclusive hypotheses. Since N. multifasciatus groups represent polygynous harems in which females have a single mating partner (i.e., the dominant male), while the dominant male may have multiple partners, the inbreeding avoidance hypothesis would predict that females should be more likely to disperse in order to avoid incestuous matings comprising a large proportion of their lifetime reproduction (but see [16, 52, 53]). Males also experience a significant handicap relative to females when attempting to switch groups because they are met with intense aggression from resident fish when acting as prospective joiners . Neolamprologus multifasciatus also express sexual dimorphism with males being larger than females . This raises the possibility that large body sizes impose higher dispersal handicaps than small body sizes in this system. For example, larger bodied fish may find it more difficult to seek shelter when venturing away from their own shells, perhaps due to the relative scarcity of large gastropod shells  into which they fit. This remains an open question for future research in N. multifasciatus and other shell-dwelling cichlids. Males experiencing higher dispersal costs than females have been detected in other taxa as well including the Seychelles warbler, Acrocephalus sechellensis, where males suffer higher mortality rates during extra-territorial forays than females . Thus, multiple factors may be working jointly to promote female-biased dispersal in N. multifasciatus, and future work will be needed to tease apart their relative contributions. Interestingly, the local mate or resource competition hypothesis [17,18,19] predicts male-biased dispersal in N. multifasciatus, because while both sexes engage in within-group competition over reproduction, reproductive skew is much higher among same-group males than females . However, as mentioned above, most subordinate males have poor options for establishing new territories of their own or being allowed to settle on a neighbouring territory, and so they cannot readily escape their local competitive environment.
Though females are the more-dispersing sex in N. multifasciatus, we also detected a marked constraint on their dispersal, namely the availability of nearby groups that might allow them to join. The shell-beds of Lake Tanganyika can be extremely densely packed with shell-dwelling cichlids , and our study population displays an average of ~ 30 cm separating nearest neighbouring N. multifasciatus territories [42, 44]. While certain regions of the shell bed offer a dispersing fish numerous destination options nearby, other regions can be more sparsely populated. We found that the more neighbouring territories there were within a reachable distance (i.e., tested up to two meters), the less related females were to their other female group-mates—a pattern that was not detected among males. Such a sex difference might arise because males are more philopatric to begin with, but also because males can establish their own territories, or forcefully take over neighbouring territories if they have grown large enough . Another, non-mutually exclusive, explanation is that in sparser regions of the shell bed, females may be discouraged from dispersing due to wider average distances of open sand separating territories acting as a movement barrier. Our results highlight the importance of accounting for the spatial distribution of suitable dispersal destinations in shaping the movement patterns of individuals. This is in agreement with previous studies that have similarly emphasized the role of dispersal site availability (e.g., great tits, Parus major, [55, 56]; lesser kestrels, Falco naumanni, ).
In conclusion, we present evidence of female-biased dispersal in the group-living and highly social cichlid N. multifasciatus, and also suggest that the spatial distribution and availability of nearby territories on which to settle represents an important ecological constraint on their movement. The genetic structure of the social environment of N. multifasciatus is highly asymmetrical, with most cohabiting females being unrelated, and most males sharing a significantly higher degree of relatedness to each other. This asymmetry in relatedness among males and females means that their social behaviour will likely generate complex fitness costs and benefits, which can in turn influence selection on decisions of when and where to disperse in an evolutionary feedback loop . Cichlid fishes showcase a great diversity of mating systems, parental care strategies , and social structures [34, 58] and we suggest that they can also be a highly valuable system for gaining new insights into the evolution and expression of sex-biased dispersal and its relationship with social evolution.
Study site and field sampling
Our study site was located in a large shell-bed on the lake floor at the southern tip of Lake Tanganyika, Zambia (8°42′49.0″ S 31°07′22.9″ E). Between September and October 2019, we identified all N. multifasciatus territories in a study quadrat on the shell bed (quadrat size measured approximately 10 × 10 m and was located at a depth of 10–11 m). Our study quadrat enclosed a cluster of territories that were separated from the rest of the shell bed by a border of open sand that was at least one meter wide. While on SCUBA, we took downward facing video footage of the entire study quadrat (GoPro Hero 7 camera set to 1080p resolution, 30 fps, and a ‘linear’ field of view) and used Structure-from-Motion photogrammetry [59, 60] to recreate the spatial layout of the territories in the study area (Fig. 1). From this spatial layout, we calculated the pairwise distances between all N. multifasciatus territories based on their Cartesian coordinates as placed in ImageJ (v 1.53e). No territories of any other shell-dwelling Lamprologine cichlids (see ) were present within this study quadrat. A buddy pair of two divers systematically sampled all individuals that could be captured from each territory. Capturing N. multifasciatus individuals can be done by hand as the fish hide within shells when approached by predators or SCUBA divers. While underwater, the shells containing hiding fish were picked up and the fish were extracted from their shells and sedated with clove oil. We then sexed each fish by inspecting their urogenital papillae, measured them for standard length (cm, SL), and recorded them as either an adult or a juvenile based on the presence of dark banding patterns along the sides of their bodies, which denote sexual maturity . Social groups were defined as the fish living in close proximity to one another, often all within a single crater excavated from the sand on the lake floor containing a collection of empty gastropod shells. These groups are also easily identified as they are physically separated from neighbouring craters constituting other N. multifasciatus territories. Fish within a social group would also share abutting or overlapping home ranges and frequently interact with one another . We recorded the largest male in each group as the dominant male and any further males within the group as subordinate males [41, 42]. We fin-clipped large fish (> 1.7 cm in SL) on their anal fins (taking at most 2 × 2 mm of tissue). When the fish had fully recovered from sedation, we returned them and their shells to their original territories. Any fish smaller than 1.7 cm (~ 20% of our whole sample) were euthanized with an overdose of clove oil and sampled whole because of the relatively large amount of tissue that fin-clipping would have removed. All tissue samples and whole fish were stored in 99% ethanol, and transported back to the lab for later microsatellite genotyping. To ensure that we had captured as many of the fish from each group as possible, we returned to each territory on at least two further occasions to check for unclipped fish. Any missed fish were similarly captured, sexed, measured, and fin-clipped. After two consecutive visits to each territory, unclipped fish were exceptionally rare, lending high confidence that we had sampled all or nearly all of the individuals living in the study quadrat. In total, we sampled 812 fish (239 adult females, 191 adult males, and 382 juveniles) from 128 territories, which constituted all territories in the study area.
Microsatellite genotyping and marker polymorphism
The microsatellite data obtained from the fish in this study quadrat were also used in other studies (e.g., 42). In the lab, DNA was extracted from the tissues using a standard Chelex protocol . All individuals were genotyped at 20 microsatellite loci divided into three multiplexes (Table 2). We used 3 µL of Qiagen Type-it Multiplex PCR Master Mix for the multiplex PCRs, along with 1 µL of template DNA, and 0.5 µL of primer mix (see Table 2 for concentrations). Total PCR volume was 5.5 µL, and each forward primer was labeled with one of the fluorescent dyes HEX, FAM, NED, ATTO550 and ATTO565. We used the following PCR program settings: denaturation at 95 °C for 5 min, followed by 30 cycles at 95 °C for 30 s, annealing at 55 °C (for multiplex1), 54 °C (for multiplex 2), or 53 °C (for multiplex 3) for 90 s, extension at 72 °C for 30 s, and a final extension at 60 °C for 30 min. We scored allele sizes against an internal standard (GeneScan 500 LIZ, Applied Biosystems) in an automatic sequencer (3130xL Genetic Analyzer, Applied Biosystems) and GeneMapper software (v 3.7, Applied Biosystems).
We estimated population allele frequencies in CERVUS (v 3.0.7; ), using a subset of fish sampled from the quadrat. To reduce the influence of within-group kinship structure, we chose the dominant male and up to one random female from each territory (N = 233 fish). The markers were highly polymorphic with an average of 16.4 alleles per locus, a mean heterozygosity of 0.75, and all markers adhered to Hardy–Weinberg equilibrium (Table 2).
All analyses in this study were carried in out in R . We used the R package “Demerelate”  to calculate symmetric Lynch-Ritland pairwise relatedness estimators (rLR, ) for each pair of fish in our dataset that had at least 10 microsatellite loci successfully genotyped (772 out of 812 or 95.1% of all individuals sampled, which was 97.9% of all adult females, 99.0% of all adult males, and 91.4% of all juveniles). We chose rLR as an estimator because it performs well for weakly related or unrelated individuals [76,77,78], which we expected to constitute the majority of cases in our dataset.
We used the R package “glmmTMB”  to fit all linear mixed effects models (LMMs), and inspected model diagnostics using the R packages “performance”  and “Dharma” . We used the R package “mgcv”  to fit our generalized additive models (GAMs).
We first examined the sex and life stage compositions of our N. multifasciatus groups, by generating frequency histograms reflecting the number of males, females, and juveniles comprising each group. We also related the number of males residing on each territory to the number of females using a linear regression model (LM).
Is sex-biased dispersal revealed by the relatedness structure of N. multifasciatus groups?
We now asked whether pairwise relatedness estimates (rLR) differed among cohabiting females, males, and juveniles. We assembled all rLR corresponding to pairs of fish living in the same groups as each other and examined the female-female, male-male, male–female, and juvenile-juvenile pairings (N = 1514 fish pairings living in 109 groups). Note that some groups were not able to contribute relatedness estimates for all pairings; for example, if a group contained only one female then we could not calculate a within-group estimate of female-female relatedness. Furthermore, if a territory was occupied by a solitary male, then this male could only contribute to between-group relatedness estimates (see below). We fit a LMM to the rLR estimates and included ‘fish pairing’ as a predictor variable (4-level categorical). Note that each rLR value corresponds to a pair of fish—Fish 1 and Fish 2—yet we opted not to include the identities of Fish 1 and Fish 2 as random effects because it was not possible to structure the random effects such that the intra- and inter-individual variation attributable to fish identities could be accounted for. Instead, we included a random intercept of ‘Group ID’. All pairwise contrasts between the levels of ‘fish pairing’ were then tested for statistical significance with the “multcomp” R package  using the Tukey method.
How does relatedness change with geographic separation between individuals?
We examined how male-male and female-female relatedness changed with geographic distance between groups. Here, we assembled all rLR corresponding to pairs of same-sex adult fish living in different groups (N = 44,246 fish pairings across all 128 territories). We used a generalized additive model (GAM) assuming a Gaussian error distribution to accommodate any potential non-linear relationships between rLR and the geographic distance between groups. We first applied a Yeo-Johnson power transformation to the rLR values to improve the normality and symmetry of the model residuals . We included ‘fish pairing’ as a predictor variable (2-level categorical: male-male vs female-female), as well as geographic distance between groups (in cm) as a cubic regression spline composed of five knots. We also looked for a sex difference in the potentially non-linear relationship between geographic distance and rLR.
Does neighbourhood density influence dispersal patterns?
We tested whether within-group relatedness among same-sex individuals varied with the number of territories in the nearby environment. For each N. multifasciatus territory in our quadrat, we calculated the number of other territories within radii of 50 cm, 100 cm, 150 cm, and 200 cm (Fig. 1). When not living in the same groups as their offspring, the vast majority of N. multifasciatus parents can be found living within 200 cm of their progeny  and so this represents an ecologically relevant spatial scale over which to look for dispersal patterns. We fit a series of linear regression models, one for each radius tested. We included average within-group female-female rLR as the response variable, and the number of neighbouring territories that each group had within their radius as a predictor variable. We then repeated this using average within-group male-male rLR values. Note that the female-female analyses only considered multi-female groups (N = 65), and the male-male analyses only used multi-male groups (N = 46). We applied a permutation-based approach to calculate P-values in which we ran 10,000 permutations of each model randomizing the response variable without replacement. P-values were calculated as the proportion of randomizations yielding slope values for the relationship between average within-group rLR and neighbourhood density that were more extreme (i.e., more negative) than our observed slopes.
Clobert J, Danchin E, Dhondt AA, Nichols JD, editors. Dispersal. Oxford: Oxford University Press; 2001.
Bowler DE, Benton TG. Causes and consequences of animal dispersal strategies: relating individual behaviour to spatial dynamics. Biol Rev. 2005;2005:205–25.
Greenwood PJ, Harvey PH. The natal and breeding dispersal of birds. Annu Rev Ecol Syst. 1982;13:1–21.
Clobert J, Ims RA, Rousset F. Causes, mechanisms and consequences of dispersal. Ecol Genet Evol Metapopulations. 2004;8:307–35.
Serrano D, Carrete M, Tella JL. Describing dispersal under habitat constraints: A randomization approach in lesser kestrels. Basic Appl Ecol. 2008;9(6):771–8.
Nevoux M, Arlt D, Nicoll M, Jones C, Norris K. The short- and long-term fitness consequences of natal dispersal in a wild bird population. Ecol Lett. 2013;16:438–45.
Davidian E, Courtiol A, Wachter B, Hofer H, Höner OP. Why do some males choose to breed at home when most other males disperse? Sci Adv. 2016;2(3):e1501236.
Trochet A, Stevens VM, Baguette M. Evolution of sex-biased dispersal. Q Rev Biol. 2016;91(3):297–320.
Hamilton WD. The genetical evolution of social behaviour. I J Theor Biol. 1964;7:1–16.
Hamilton WD. The genetical evolution of social behaviour. II J Theor Biol. 1964;7:17–52.
West SA, Pen I, Griffin AS. Cooperation and competition between relatives. Science. 2002;296:72–5.
Kay T, Lehmann L, Keller L. Kin selection and altruism. Curr Biol. 2019;29(11):R438-442.
Gardner A. Sex-biased dispersal of adults mediates the evolution of altruism among juveniles. J Theor Biol. 2010;262(2):339–45.
Schweizer M, Excoffier L, Heckel G. Fine-scale genetic structure and dispersal in the common vole (Microtus arvalis). Mol Ecol. 2007;16(12):2463–73.
Kingma SA, Komdeur J, Burke T, Richardson DS. Differential dispersal costs and sex-biased dispersal distance in a cooperatively breeding bird. Behav Ecol. 2017;28(4):1113–21.
Li XY, Kokko H. Sex-biased dispersal: a review of the theory. Biol Rev. 2019;94(2):721–36.
Greenwood PJ. Mating systems, philopatry and dispersal in birds and mammals. Anim Behav. 1980;28(4):1140–62.
Stephen DF. Competition for mates and predominant juvenile male dispersal in mammals. Anim Behav. 1982;30(4):1183–92.
Perrin N, Mazalov V. Local competition, inbreeding, and the evolution of sex-biased dispersal. Am Nat. 2000;155(1):116–27.
Hamiltion WD, May RM. Dispersal in stable habitats. Nature. 1977;269(5629):578–81.
Pusey AE. Sex-biased dispersal and inbreeding avoidance in birds and mammals. Trends Ecol Evol. 1987;2(10):295–9.
Perrin N, Goudet J. Inbreeding, kinship, and the evolution of natal dispersal. In: Clobert J, Danchin E, Dhondt AA, Nichols JD, editors. Dispersal. Oxford: Oxford University Press; 2001. p. 123–42.
Bonte D, Van Dyck H, Bullock JM, Coulon A, Delgado M, Gibbs M, et al. Costs of dispersal. Biol Rev. 2012;87(2):290–312.
Mabry KE, Shelley EL, Davis KE, Blumstein DT, van Vuren DH. Social mating system and sex-biased dispersal in mammals and birds: A phylogenetic analysis. PLoS ONE. 2013;8(3):1–9.
Avise JC, Jones AG, Walker D, DeWoody JA. Genetic mating systems and reproductive natural histories of fishes: lessons for ecology and evolution. Annu Rev Genet. 2002;36(1):19–45.
Coleman SW, Jones AG. Patterns of multiple paternity and maternity in fishes. Biol J Linn Soc. 2011;103(4):735–60.
Balshine S, Sloman KA. Parental Care in Fishes. In: Encyclopedia of Fish Physiology: From Genome to Environment. 2011. p. 670–677.
Renn SCP, Hurd PL. Epigenetic regulation and environmental sex determination in cichlid fishes. Sex Dev. 2021;15(1–3):93–107.
Hutchings JA, Gerber L. Sex-biased dispersal in a salmonid fish. Proc R Soc B Biol Sci. 2002;269:2487–93.
Stiver KA, Desjardins JK, Fitzpatrick JL, Neff B, Quinn JS, Balshine S. Evidence for size and sex- specific dispersal in a cooperatively breeding cichlid fish. Mol Ecol. 2007;16(14):2974–84.
Cano JM, Mäkinen HS, Merilä J. Genetic evidence for male-biased dispersal in the three-spined stickleback (Gasterosteus aculeatus). Mol Ecol. 2008;17(14):3234–42.
Yue GH, Xia JH, Liu F, Lin G. Evidence for female-biased dispersal in the protandrous hermaphroditic Asian seabass. Lates calcarifer PLoS One. 2012;7:6.
Jordan A, Taborsky B, Taborsky M. Cichlids as a model system for studying social behaviour and evolution. In: Abate ME, Noakes DLG, editors. The Behavior, Ecology and Evolution of Cichlid Fishes. Fish & Fisheries Series 40. Springer, Dordrecht; 2021. p. 587–635.
Lein E, Jordan A. Studying the evolution of social behaviour in Darwin’s Dreampond – a case for the Lamprologine shell cichlids. Hydrobiologia. 2021;848:3699–726.
Knight ME, Van Oppen MJH, Smith HL, Rico C, Hewitt GM, Turner GF. Evidence for male-biased dispersal in Lake Malawi cichlids from microsatellites. Mol Ecol. 1999;8(9):1521–7.
Anseeuw D, Maes GE, Busselen P, Knapen D, Snoeks J, Verheyen E. Subtle population structure and male-biased dispersal in two Copadichromis species (Teleostei, Cichlidae) from Lake Malawi. East Africa Hydrobiologia. 2008;615(1):69–79.
Lehtonen TK, Elmer KR, Lappalainen M, Meyer A. Genetic evidence for panmixia in a colony- breeding crater lake cichlid fish. Sci Rep. 2018;8(1):1–8.
Taylor MI, Morley JI, Rico C, Balshine S. Evidence for genetic monogamy and female-biased dispersal in the biparental mouthbrooding cichlid Eretmodus cyanostictus from Lake Tanganyika. Mol Ecol. 2003;12(11):3173–7.
van Dongen WFD, Wagner RH, Moodley Y, Schaedelin FC. Sex biases in kin shoaling and dispersal in a cichlid fish. 2014. p. 965–974.
Tanaka H, Heg D, Takeshima H, Takeyama T, Awata S, Nishida M, et al. Group composition, relatedness, and dispersal in the cooperatively breeding cichlid Neolamprologus obscurus. Behav Ecol Sociobiol. 2015;69(2):169–81.
Kohler U. Zur Struktur und Evolution des Sozialsystems von Neolamprologus multifasciatus (Cichlidae, Pisces), dem kleinsten Schneckenbuntbarsch des Tanganjikasees. Max-Planck Institut für Verhaltensphysiologie in Seewiesen. 1998.
Bose APH, Dabernig-Heinz J, Koch L, Grimm J, Lang S, Hegedűs B, et al. Parentage analysis across age cohorts reveals sex differences in reproductive skew in a group-living cichlid fish. Neolamprologus multifasciatus Mol Ecol. 2022;8:67.
Schradin C, Lamprecht J. Female-biased immigration and male peace-keeping in groups of the shell-dwelling cichlid fish Neolamprologus multifasciatus. Behav Ecol Sociobiol. 2000;48(3):236–42.
Jordan LA, Maguire SM, Hofmann HA, Kohda M. The social and ecological costs of an “over- extended” phenotype. Proc R Soc B. 2016;283:20152359.
Bose APH, Nührenberg P, Jordan A. Female-female conflict is higher during periods of parental care in a group-living cichlid fish. Anim Behav. 2021;182:91–105.
Gübel J, Bose APH, Jordan A. Social and spatial conflict drive resident aggression toward outsiders in a group-living fish. Behav Ecol. 2021;9:1–9.
Prugnolle F, de Meeus T. Inferring sex-biased dispersal from population genetic tools: a review. Heredity (Edinb). 2002;88:161–5.
Lynch M, Ritland K. Estimation of pairwise relatedness with molecular markers. Genetics. 1999;152(4):1753–66.
Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biometrical J. 2008;50(3):346–63.
Kokko H, Johnstone RA, Wright J. The evolution of parental and alloparental effort in cooperatively breeding groups: when should helpers pay to stay? Behav Ecol. 2002;13(3):291–300.
Hamilton IM, Taborsky M. Unrelated helpers will not fully compensate for costs imposed on breeders when they pay to stay. 2005;272:445–54.
Lehmann L, Perrin N. Inbreeding avoidance through kin recognition: Choosy females boost male Dispersal. Am Nat. 2003;162(5):638–52. https://doi.org/10.1086/378823.
De BRA, Vega-trejo R, Kotrschal A, Fitzpatrick JL. Meta-analytic evidence that animals rarely avoid inbreeding. Nat Ecol Evol. 2021;5:949–64.
Bose APH, Windorfer J, Böhm A, Ronco F, Indermaur A, Salzburger W, et al. Structural manipulations of a shelter resource reveal underlying preference functions in a shell-dwelling cichlid fish. Proc R Soc B. 2020;287:20200127.
van Noordwijk AJ. Quantitative Genetics in Natural Populations of Birds Illustrated with Examples from the Great Tit, Parus major. In: Wöhrmann K, Loeschcke V, editors. Population Biology and Evolution Proceedings in Life Sciences. Berlin: Springer; 1984. p. 67–79.
Vardakis M, Goos P, Adriaensen F, Matthysen E. Discrete choice modelling of natal dispersal: ‘Choosing’ where to breed from a finite set of available areas. Methods Ecol Evol. 2015;6:997–1006.
Sefc KM. Mating and Parental Care in Lake Tanganyika’s Cichlids. Int J Evol Biol. 2011;2011:e470875.
Fryer G, Iles TD. The cichlid fishes of the Great Lakes of Africa. Edinburgh: Oliver and Boyd; 1972.
Westoby MJ, Brasington J, Glasser NF, Hambrey MJ, Reynolds JM. ‘Structure-from-Motion’ photogrammetry: A low-cost, effective tool for geoscience applications. Geomorphology. 2012;179:300–14.
Jungwirth A, Nührenberg P, Jordan A. On the importance of defendable resources for social evolution: Applying new techniques to a long-standing question. Ethology. 2021;127:872–85.
Walsh PS, Metzger DA, Higuchi R. Chelex 100 as a medium for simple extraction of DNA for PCR- based typing from forensic material. Biotechniques. 1991;10(4):506–13.
Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16:1099–106.
Crispo E, Hagen C, Glenn T, Geneau G, Chapman LJ. Isolation and characterization of tetranucleotide microsatellite markers in a mouth-brooding haplochromine cichlid fish (Pseudocrenilabrus multicolor victoriae) from Uganda. Mol Ecol Notes. 2007;7(6):1293–5.
Carleton KL, Streelman JT, Lee BY, Garnhart N, Kidd M, Kocher TD. Rapid isolation of CA microsatellites from the tilapia genome. Anim Genet. 2002;33(2):140–4.
Lee BY, Lee WJ, Streelman JT, Carleton KL, Howe AE, Hulata G, et al. A second-generation genetic linkage map of tilapia (Oreochromis spp). Genetics. 2005;170(1):237–44.
Taylor MI, Meardon F, Turner G, Seehausen O, Mrosso HDJ, Rico C. Characterization of tetranucleotide microsatellite loci in a Lake Victorian, haplochromine cichlid fish: A Pundamilia pundamilia x Pundamilia nyererei hybrid. Mol Ecol Notes. 2002;2(4):443–5.
Maeda K, Takeshima H, Mizoiri S, Okada N, Nishida M, Tachida H. Isolation and characterization of microsatellite loci in the cichlid fish in Lake Victoria Haplochromis chilotes. Mol Ecol Resour. 2008;8(2):428–30.
Lee WJ, Kocher TD. Microsatellite DNA markers for genetic mapping in Oreochromis niloticus. J Fish Biol. 1996;49(1):169–71.
Parker A, Kornfield I. Polygynandry in Pseudotropheus zebra, a cichlid fish from Lake Malawi. Environ Biol Fishes. 1996;47(4):345–52.
Zardoya R, Vollmer DM, Craddock C, Streelman JT, Karl S, Meyer A. Evolutionary conservation of microsatellite flanking regions and their use in resolving the phylogeny of cichlid fishes (Pisces: Perciformes). Proc R Soc B Biol Sci. 1996;263:1589–98.
Albertson RC, Streelman JT, Kocher TD. Directional selection has shaped the oral jaws of Lake Malawi cichlid fishes. Proc Natl Acad Sci U S A. 2003;100(9):5252–7.
Brandtmann G, Scandura M, Trillmich F. Female-female conflict in the harem of a snail cichlid (Lamprologus ocellatus): Behavioural interactions and fitness consequences. Behaviour. 1999;136:1123–44.
Van Oppen MJH, Rico C, Deutsch JC, Turner GF, Hewitt GM. Isolation and characterization of microsatellite loci in the cichlid fish Pseudotropheus zebra. Mol Ecol. 1997;6(4):387–8.
R Core Team. R: a language and environment for statistical computing. Vienna (Austria): R Foundation for Statistical Computing. Vienna (Austria): R Foundation for Statistical Computing. http://www.R-project.org/. 2020.
Kraemer P, Gerlach G. Demerelate: calculating interindividual relatedness for kinship analysis based on codominant diploid genetic markers using R. Mol Ecol Resour. 2017;17(6):1371–7.
Bose APH, Henshaw JM, Zimmermann H, Fritzsche K, Sefc KM. Inclusive fitness benefits mitigate costs of cuckoldry to socially paired males. BMC Biol. 2019;17(2):1–16.
Van De Casteele T, Galbusera P, Matthysen E. A comparison of microsatellite-based pairwise relatedness estimators. Mol Ecol. 2001;10(6):1539–49.
Csilléry K, Johnson T, Beraldi D, Clutton-Brock T, Coltman D, Hansson B, et al. Performance of marker-based relatedness estimators in natural populations of outbred vertebrates. Genetics. 2006;173(4):2091–101.
Brooks ME, Kristensen K, van Benthem KJ, Magnusson A, Berg CW, Nielsen A, et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 2017;9(2):378–400.
Lüdecke D, Makowski D, Waggoner P, Patil I. Assessment of regression models performance. CRAN. 2020. https://easystats.github.io/performance/
Hartig F. DHARMa: Residual diagnostics for hierarchical (multi-level / mixed) regression models. CRAN. 2020. https://cran.r-project.org/package=DHARMa
Wood SN. Generalized Additive Models: An Introduction with R. New York: CRC Press; 2017.
Yeo I-K, Johnson RA. A new family of power transformations to improve normality or symmetry. Biometrika Trust. 2000;87(4):954–9.
We are grateful to the Department of Fisheries in Mpulungu, Zambia for supporting our research at Lake Tanganyika. We also thank the staff at the Tanganyika Science Lodge for their hospitality and our guides and boat drivers for their incredibly hard work.
Open Access funding enabled and organized by Projekt DEAL. This research was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—EXC 2117—422037984 and by the Department of Collective Behaviour, Max Planck Institute of Animal Behavior, as well as by the Austrian Science Fund (FWF, Grant Number P 27605-B25 to KMS) and the University of Graz. AB was further supported by an Alexander von Humboldt Research Fellowship and a Young Scholars Fund grant from the University of Konstanz.
Ethics approval and consent to participate
All methods adhered to the ASAB/ABS Guidelines for the Use of Animals in Research. All work is reported in accordance with ARRIVE guidelines. Field work was carried out with approval from the Fisheries Department at the Ministry of Fisheries and Livestock Zambia, under study permits issued by the government of Zambia (Nos. G7067690 and C3195368) and in conjunction with a memorandum of understanding with the University of Zambia (MOU 101/14/11). The study species is listed as ‘Least Concern’ on the IUCN Red List of Threatened Species.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Between-group data. Pairwise Lynch-Ritland relatedness estimates for each pairing of fish that reside in different social groups from one another. Dataset also includes the spatial distances separating social groups (measured in cm).
R script. R script to replicate all analyses in this paper.
Group density data. Dataset containing the number of social groups found near each focal group (neighbourhood density as measured over 50 cm, 100 cm, 150 cm, and 200 cm). Dataset also contains the average within-group relatedness estimates (pairwise Lynch-Ritland) for both male-male and female-female comparisons.
Within-group data. Pairwise Lynch-Ritland relatedness estimates for each pairing of fish that reside together in the same social group.
About this article
Cite this article
Bose, A.P.H., Koch, L., Dabernig-Heinz, J. et al. Patterns of sex-biased dispersal are consistent with social and ecological constraints in a group-living cichlid fish. BMC Ecol Evo 22, 21 (2022). https://doi.org/10.1186/s12862-022-01980-4
- Relatedness estimators
- Microsatellite genotyping
- African cichlid
- Habitat saturation
- Social conflict
- Population structure
- Ecological constraints