The 1960s and 1970s saw the emergence of two research topics within evolutionary biology that could be described as revolutions. Kin selection provided an explanation for the evolution of altruistic traits and more generally for the evolution of social behavior (Hamilton, 1963, 1964a,b). Meanwhile, post-copulatory sexual selection extended the scope of sexual selection from Darwin’s view of pre-copulatory sexual selection (Darwin, 1871) to selection taking place after copulation, or more generally, after ejaculation, with consequences that rival those of pre-copulatory sexual selection (Parker, 2020). A key form of post-ejaculatory sexual selection is sperm competition (Parker, 1970), with far-reaching consequences for the evolution of both gametic and organismal traits (Parker, 2020). Kin selection and sperm competition have both since become central research areas within evolutionary biology (Gardner et al., 2011; Lehtonen, 2020a; Parker, 2020). Whilst the intersection of sexual selection, sexual conflict, and kin selection has been the focus of considerable theoretical attention in recent years (e.g., Rankin, 2011; Wild et al., 2011; Faria et al., 2015, 2020), there has been relatively little theoretical work focusing explicitly on the intersection of sperm competition and kin selection. In this article, we discuss ways in which the two fields intersect, and how sperm competition, as well as the related topics of local sperm competition, sperm limitation, and sex allocation can be viewed from a kin selection perspective. We focus on an overview of these broad themes, taking a theoretical evolutionary ecology perspective (for a detailed discussion of ‘sperm sociality’ (see, e.g., Pizzari and Foster, 2008; Higginson and Pitnick, 2011). We then derive a novel model of sex allocation in simultaneous hermaphrodites that makes use of and illuminates all these topics.

In kin selection theory, fitness effects on genetic relatives are explicitly accounted for: individuals acting on their relatives indirectly affect the propagation of copies of their own genes (Hamilton, 1964a). An individual can increase its inclusive fitness by increasing either its direct fitness effect by its own reproductive success or its indirect fitness by increasing the reproductive success of its genetic relatives (Hamilton, 1964a). A central aspect of kin selection theory is the relatedness coefficient that measuress the relative quality of individuals as propagators of the focal individual’s genes. When related individuals interact, the evolution of altruistic traits can be promoted, and the evolution of selfish traits can be restricted to the extent that it maximises the overall propagation of the focal individual’s genes into future generations, via its own offspring (direct fitness) as well as those of other, related individuals (indirect fitness).

Considering relatedness from the gametic perspective, Trivers expressed the question: ‘If close degrees of relatedness predispose individuals to altruism, why not search for such examples in sex cells, specifically sperm?’ (Trivers, 1985), justifying the question with the fact that sperm from a single male are related to each other by a relatedness coefficient of ½ (under outbreeding). This quote illustrates the reason why kin selection is relevant and useful in thinking about sperm competition. Kin selection is important in interactions between organisms when related individuals, such as siblings interact. By similar logic, kin selection is relevant in interactions between male gametes even when the sperm descend from the same male. Therefore, as the quote by Trivers suggests, in competition between gametes such kin interactions are almost inevitable. From such considerations follows our central message: if we shift our perspective to gametes as analogues of interacting individuals, sperm competition models can typically be interpreted as kin selection models. Furthermore, a kin selection perspective can bring with it a unified structure across models that is readily extendable to situations where competing gametes are related beyond a simple sibling relationship. An example of this kind of scenario is the sex allocation model we present in this paper, where relatedness between competing gametes depends on the rate and cost of dispersal in the population.

Trivers’ quote above seems to attribute control of the sperm phenotype to the haploid genotype. The extent of haploid control in gametes is an active area of ongoing research, and the view that haploid transcription is insignificant in animals has in recent years been challenged and at least partially overturned (Immler, 2019). However, current evidence still shows that haploid expression is much more extensive in plants than it is in animals (Immler, 2019) and the prevailing view seems to be that sperm form and function is mainly under diploid control (e.g., Higginson and Pitnick, 2011). Thus, whilst its importance remains debated, the evidence is clear that haploid gene expression does exist and that there is a continuum of expression levels.

The theoretical predictions of sperm competition theory depend crucially on whether the evolution of sperm behavior is controlled completely by the male (diploid control: Parker, 1993) or by the male and the sperm, the sperm having sufficiently strong haploid gene expression. Haploid control (sperm phenotype driven by the spermatozoon itself) has been modelled by Parker and Begon (1993). Under haploid control, there is potential for conflict between the evolutionary interests of the parent and gamete, because there is competition both within and between ejaculates, whilst under diploid (parental) control there is competition only between ejaculates from different males. The magnitude of this conflict decreases as the magnitude of inter-ejaculate competition increases (Parker and Begon, 1993). Typically, the term ‘sperm competition’ refers to inter-ejaculate competition, and this default meaning is advocated by Parker and Begon (1993), whilst ‘inter-ejaculate’ and ‘intra-ejaculate’ sperm competition can be used where necessary to differentiate evolution under haploid expressed sperm competition.

This can also be seen from a kin selection perspective: the appropriate kin selection relatedness coefficient depends on the ‘control individual’ (Taylor, 1988a). If control lies with the gamete (haploid control), then the coefficient of ½ between sibling gametes in Trivers’ quote is correct: from a focal gamete’s perspective it measures how likely a sibling gamete is to transmit a copy of the same allele as the focal does, and under outbreeding this probability is ½. It is essentially a measure of the similarity of the sibling gamete to the focal gamete. However, if the parent has control, the relevant relatedness coefficient is instead a measure of the relative similarity of the two gametes to the control individual (Taylor, 1988a), and this relatedness coefficient equals 1: from the parent’s perspective, one of its own gametes is not any better than its sibling as a propagator of genes. Therefore, sperm cells should be predisposed to altruism as Trivers suggests, but under the typical assumption of diploid control this is due to a relatedness coefficient of 1, rather than a relatedness coefficient ½. We should therefore be careful in directly comparing the evolution of cooperation between gametes to, e.g., that between siblings. Pizzari and Foster (2008) remark that if only one male has mated with a female, from the perspective of the diploid genome of the male parent all sperm competing for fertilisation are equally valuable, whilst different sperm haplotypes from the same male are in evolutionary conflict with each other. As noted above, under haploid control, the conflict between the evolutionary interests of the parent and gamete decreases as inter-ejaculate sperm competition increases (Parker and Begon, 1993). A kin selection interpretation of this effect is that as the number of individuals contributing to sperm competition increases, the relatedness of a focal spermatozoon to a random competitor spermatozoon decreases, regardless of whether control lies with the parent or gamete: the fraction of unrelated competitors increases in both cases. However, conflict does not necessarily disappear entirely even under extremely high sperm competition (Parker and Begon, 1993). The kin selection explanation for this is that whilst relatedness between sperm competing for fertilisations may approach 0 under intense sperm competition, we must additionally consider relatedness between spermatozoa competing for resources, and these are still sibling gametes regardless of the level of sperm competition. Regardless of the intensity of sperm competition, it may pay off for mutant sperm under haploid control to sap resources from non-mutant sperm of the same parent (with relatedness ½ between sibling gametes), whilst under diploid control (with relatedness 1 from the parent’s perspective) this does not pay off.

The above applies to situations where the parent individuals are not necessarily related to each other. However, the parent individuals can of course also be related to each other, and modelling suggests that sperm expenditure decreases with increasing relatedness between males in sperm competition (Parker, 2000). This is a very intuitive outcome: it does not pay off (or at least is not as valuable) to invest into competition with individuals who are likely to pass on copies of the focal individual’s genes anyway. Empirical tests of this theory have, however, so far not found clear evidence for this prediction (e.g., Thomas and Simmons, 2008). On the theoretical side, Parker and Pizzari (2010) note that Parker’s (2000) models may have implicitly assumed that the population had low viscosity (i.e., high dispersal). It may be worth pursuing these relatively understudied questions further in future work, accounting explicitly for factors such as dispersal rate.

Sperm competition refers to competition between sperm for eggs, typically assuming that the total number of fertilisations is fixed and limited only by the number of eggs. However, the fertilisation success of eggs may depend on the number or density of spermatozoa, and this problem may be particularly acute in marine broadcast spawner settings (Levitan and Petersen, 1995; Yund, 2000). This concept of sperm limitation (or more generally, gamete limitation) can be described mathematically using fertilisation functions (an early example can be found in Rothschild and Swann, 1951, and for reviews see Lehtonen and Dardare, 2019, and Crimaldi and Zimmer, 2014 for more complex physical settings).

Sperm limitation can interact with sperm competition, together driving sperm and ejaculate evolution (Ball and Parker, 1996). Similarly, it has been argued that gamete limitation can drive the evolution of anisogamy (i.e., the ancestral divergence of the sexes), either by itself (Iyer and Roughgarden, 2008; Roughgarden and Iyer, 2011) or in conjunction with gamete competition (Lehtonen and Kokko, 2011; Lehtonen and Parker, 2014). The latter models unify the proposed evolutionary pressures of gamete competition and gamete limitation in a common framework, and in these models local spawning group size modulates the relative evolutionary pressure via gamete competition and gamete limitation. This can again be viewed through a kin selection lens: group size influences the relatedness of gametes within the group, and the average relatedness of a parent to those gametes. Consider an extreme case of a group of one adult spawning individual. Then, all gametes within the group will be siblings, and the parent will be maximally related to the gametes in the group. Then again, in very large groups, or in panmictic populations with no group structure, the vast majority of gamete interactions occur between unrelated individuals. The relative importance of gamete limitation as a selective pressure becomes emphasised in small groups whilst that of gamete competition becomes reduced (Parker and Lehtonen, 2014). Taking a step further, it has been suggested that more generally and beyond the group size effect, relatedness within groups modulates the selective effects of gamete competition and gamete limitation (Lehtonen, 2016).

However, this last argument comes with a caveat not mentioned in the article (Lehtonen, 2016). Beyond group size, the simplest and most likely source of increased relatedness in a mating group is limited dispersal between groups: if some offspring remain in their natal patch to breed in the next generation, relatedness within the patch increases. However, classic studies in the kin selection literature show that under such circumstances, there is limited scope for increased relatedness to promote the evolution of altruistic traits. The reason for this is the secondary effect of increased competition between relatives: if one individual helps a related individual make more offspring, and some of those offspring remain within the natal patch, they increase competition for reproductive opportunities within the patch. Under many circumstances, these effects cancel each other exactly (Taylor, 1992a,b). Similarly, in models of sex ratio evolution in organisms with separate sexes, the effect of increased relatedness and increased competition (both due to limited dispersal) tend to cancel each other out (Bulmer, 1986; Frank, 1986; Taylor, 1988b), so that relatedness has no effect on the stable sex ratio. It therefore seems plausible that a similar effect might arise, at least in some circumstances, with gamete limitation. An example of this will be seen in the model of hermaphrodite sex ratio evolution we present later.

A model of sex allocation in simultaneously hermaphroditic organisms (Charnov, 1980) predicted that a larger fraction of reproductive resources should be invested in egg production (and conversely, less in sperm production) when the mating group size is small. Conversely, with increasing mating group size sex allocation becomes less female gamete biassed, approaching a 1:1 sex allocation ratio in large groups. Fischer (1981) derived a similar conclusion, apparently independent from (Charnov, 1980; see also West, 2009, p. 83 and references therein).

At its core this model is built on gamete interactions, including sperm competition. In accordance with the central message stated above, the logic of these models is analogous to models of sex ratio evolution under local mate competition in gonochoristic species (Hamilton, 1967; West, 2009). The main difference is that the ‘local competition’ for fertilisations happens on the level of gametes in simultaneous hermaphrodites, instead of the level of adult organisms in organisms with separate sexes. Due to the analogy, Schärer (2009) has proposed the term ‘local sperm competition’ (LSC) to describe such situations in hermaphrodites.

From a broader, kin selection perspective, in simultaneous hermaphrodites female-biassed sex allocation arises when competition occurs between related sperm for the fertilisation of ova of the mating partner (Schärer, 2009), emphasising the fact that the sperm need not be siblings and instead it suffices that they are related (Schärer, 2009). Relatedness between sperm in an ejaculate is thus analogous to the relatedness between gonochoristic brothers, but other factors may also increase relatedness – in particular limited dispersal of zygotes or adult individuals after fertilisation. However, it is again uncertain how this effect would interact with the secondary effects of increased local competition for breeding spots that must be taken into account (see previous section). The model on hermaphrodite sex ratio evolution that follows shows that in a manner similar to separate sexes, primary and secondary effects tend to cancel out, so that local competition between related gametes need not influence the sex ratio, beyond the effect driven by group size alone.

In some sense, local sperm competition can be thought to be the inverse of sperm competition (Singh and Schärer, 2022): under selfing (or in the absence of selfing, in mating groups of two individuals), there is no sperm competition and maximally strong LSC; in slightly larger mating groups, there is weak sperm competition and weak LSC; and in large mating group size, there is strong sperm competition and no LSC (Schärer, 2009). Thus, increasing mating group size, that is, increasing the number of unrelated sperm competing to fertilise the eggs, increases sperm competition between unrelated ejaculates and increases the investment to sperm (Charnov, 1980).

We now present a novel model combining kin selection theory with several classical features of the theory of gametic interactions: gamete competition, local gamete competition, gamete limitation, and allocation between male and female gametes.

Consider simultaneous hermaphrodites reproducing in numerous patches of n individuals. This could correspond to either internal fertilisers where one individual transfers sperm to another (e.g., many species of barnacles which informed the model of Charnov, 1980) or broadcast spawners or spermcasters, where either one type or both types of gametes are released into the water column (common modes of fertilisation in hermaphrodites: Strathmann et al., 1984; Henshaw et al., 2014; Parker et al., 2018). Reproduction via obligate outcrossing and possible self-fertilisation are considered separately. A fraction f of eggs in a patch are fertilised, and this fraction may depend on the density of spermatozoa available for fertilisation (in other words, f corresponds to a ‘fertilisation function’ Lehtonen and Dardare, 2019). After fertilisation but prior to reproducing themselves, a fraction d of zygotes disperse to a different, random patch, whilst a fraction (1-d) remain in their natal patch. Dispersal incurs a cost, so that a fraction c of the dispersers perish. After this partial dispersal event, the mature individuals in each patch compete for n breeding spots, so that each patch is again culled to n individuals. Fertilisation amongst these n individuals takes place and the life cycle restarts.

This scenario differs from and expands upon the classic model of Charnov (1980) in numerous ways, but contains it as a special case where fertilisation is internal, there is no selfing, all gametes are fertilised (f = 1), and all individuals disperse after fertilisation (d = 1). Internal vs. external fertilisation does not significantly change the model and altering the possibility of selfing is relatively straightforward (see Appendix). However, the combined additions of sperm limitation (with f < 1 and dependent on sperm density) and limited, costly dispersal bring with them new challenges which we address in the model that follows.

Taking a kin selection approach, the model could be constructed in two equivalent ways: direct fitness or inclusive fitness (Taylor et al., 2007). Here we present a direct fitness derivation of the random fertilisation scenario where both outbreeding and selfing are possible, and an outcrossing-only model is presented in the Appendix. With the direct fitness approach, we consider fitness via male and female gametes of a focal parent, accounting for all the effects on those gametes due to the sex ratio strategy of the focal parent itself, due to those of everyone else in their local group, and to those of the population as a whole. For clarity of exposition, we construct the model in sequential steps and eventually pull all those pieces together in a single direct fitness function. The notation is presented in Table 1.

We begin with the reproductive output of a focal individual in its focal patch:

For transparency, we have included variables which cancel out in the final expression. The logic of equation (1) is as follows: Each individual has a total gametic resource budget of M units which it can allocate into female and male gametes. Variables x and y denote the sex allocation trait (fraction of resources allocated into male gametes) of the focal individual, and its mean value in the focal patch, respectively. The focal individual thus allocates ( 1 − x ) M uwnits of resources into female gametes, making a total of ( 1 − x ) M / m e gametes of size m e . Total egg production in the patch is n ( 1 − y ) M / m e , and a fraction f(y) of all eggs are fertilised, potentially dependent on mean sperm production in the patch (y). The focal individual makes x M / m s sperm of size m_s, and under fair raffle sperm competition (Parker and Pizzari, 2010) they gain a fraction x M / m s n y M / m s of all fertilised eggs in the patch.

The reproductive output in its focal patch (w_p) only accounts for the success of the focal individual’s gametes in its patch prior to dispersal. However, this by itself is not sufficient. Dispersing offspring are additionally affected by the cost of dispersal. Furthermore, both non-dispersers and surviving dispersers must additionally make it through population regulation in the patch they settle in, where the probability of survival for a non-disperser is

Here z indicates the mean value of the sex allocation trait over the whole population. The symbol n in the denominator indicates the number of breeding spots being competed over, whilst the denominator indicates the total number of individuals competing in the natal patch, including both non-dispersers and surviving dispersers from other patches. For a disperser who has already made it through the dispersal stage, the probability of surviving through population regulation is

Combining these fitness components and simplifying, we find that many constants cancel out, and we are left with

At this point we can perform a brief check by calculating fitness in a neutral population. Setting (x = y = z) and simplifying, we find w = 2, as it should be: in a population regulated at stable size and where all individuals are identical, in each generation a parent must on average achieve one unit of fitness via its female gametes and one via its male gametes.

To find the stable sex ratio based on this fitness function, we then need to infer how the fitness of the focal parent, achieved via its own male and female gametes, is influenced by a small change in sex allocation of the focal individual, and by a small change in the average allocation strategy in the focal group, weighted by the appropriate relatedness coefficients. The fitness effects are the partial derivatives ∂ w ∂ x | x = y = z and ∂ w ∂ y | x = y = z (Taylor and Frank, 1996). Then, if the relatedness of the focal individual to its own gametes is r_o and the expected relatedness of the focal individual to a randomly picked gamete in its patch is r_p, we can find the candidate ESS sex ratio by solving the equation

or equivalently,

where R = r_p/r_o is the relatedness of a random gamete relative to the relatedness to own gametes, from the perspective of the focal parent (who controls the phenotype, but gains fitness only via the gametes). This relatedness coefficient in turn is equal to the consanguinity (Bulmer, 1994) of the parent with a random gamete relative to the consanguinity of the parent with one of its own gametes, p_p/p_o. A random haploid gamete represents an allele picked randomly from the parental population, so that p_p must be equal to the consanguinity p_rg between two random gametes in the same generation. Similarly, a haploid gamete of the parent itself represents an allele picked randomly from the parent, so that p_o must be equal to the consanguinity p_sg between two random sibling gametes (i.e., gametes from the same parent).

Therefore, to find R it suffices to find the ratio R = p_rg/p_sg, which can be computed using the recursion

where 1 n is the probability that two random gametes have the same parent, and k is the probability that the parents of two randomly picked gametes originate from the focal patch, i.e., k = 1 − d 1 − d + ( 1 − c ) d = 1 − d 1 − c d . Dividing through by p_sg we have

from which we can solve

Similar equations have been derived for various biological scenarios. For example, equation (13) in Taylor (1988a) for the relatedness between two mated diploid females on the same patch is identical in form to equation (9) above.

The partial derivatives ∂ w ∂ x | x = y = z and ∂ w ∂ y | x = y = z are relatively complex, but nevertheless straightforward to calculate using standard rules of differentiation (e.g., with computer software such as Mathematica). Putting everything together and simplifying, we find that the necessary ESS condition is

where q = ( c − 1 ) d ( d + c d − 2 ) ( 1 − 2 d ( 1 + ( c − 1 ) n ) + d 2 ( 1 + ( c 2 − 1 ) n ) ) ( z − 1 ) z f ( z ) . Clearly, any solutions to the equation with respect to z are determined solely by the factor in square brackets, so it is sufficient to solve

or equivalently,

where f ′ ( z ) f ( z ) could alternatively be written as d d z ln ( f ( z ) ) .

When fertilisation is limited by sperm abundance, the derivative on the right side of equation (12) is positive (i.e. an increase in sperm concentration increases fertilisation success). Sperm limitation f ′ ( z ) > 0 promotes less female-biased sex allocation (Figure 1), thus increasing investment into sperm and improving the fertilisation prospects of eggs. This result is analogous to that of Lehtonen and Schwanz (2018) for sex allocation under mate limitation and local mate competition, but with full dispersal (i.e. d = 1) in organisms with separate sexes. This has biological relevance particularly for broadcast spawning and spermcasting simultaneous hermaphrodites, given the relative prevalence of gamete limitation in the sea (Levitan and Petersen, 1995; Yund, 2000).

Whilst the effect of gamete limitation is interesting, it is equally interesting to note that c and d are both absent from equations (11, 12). Limited dispersal and increased cost of dispersal increase the relatedness of gametes within a patch, and one might expect relatedness in turn to influence the stable sex ratio. Yet no such effect is seen in the result. This is an example of the cancellation effect of increased relatedness vs. increased competition for reproductive spots amongst related females, noted in studies of sex ratio evolution in organisms with separate sexes and no sperm limitation or mate limitation (Bulmer, 1986; Frank, 1986; Taylor, 1988b). Our results show that the cancellation effect extends to simultaneous hermaphrodites, even under sperm limited conditions. Lehtonen (2016) argued that relatedness determines the relative importance of gamete limitation and gamete competition as selective forces. However, our analysis shows that the question is more nuanced when relatedness is modulated by dispersal. Whilst it is true that the effect of sperm limitation is modulated by relatedness, in the present model this effect is cancelled by secondary, competitive effects in the manner described above.

In the absence of sperm limitation (i.e. f’(z)=0) equation (12) reduces to 1 + n(2z − 1) = 0, with the solution

the classic result of sex allocation in hermaphrodites under local sperm competition when both selfing and outcrossing are allowed. We have shown that the result is not altered by increased relatedness under limited dispersal, analogous to the classic result in sex ratio evolution in organisms with separate sexes (e.g., Taylor, 1988b).

Under strict outcrossing (Appendix) we find the solution

in the absence of sperm limitation, and similarly find that it is independent of dispersal and of mortality under dispersal, and that sperm limitation drives a less female-biassed sex allocation (see Appendix). Note that equation (14) is equivalent to Charnov’s (1980) result: where we denote the total number of individuals in the patch with n, Charnov instead uses k to denote the number of competitors in the patch, aside from the focal individual itself. Therefore, k = n − 1 and substituting n = k + 1 into equation (14) we recover k − 1 2 k − 1 as was found by Charnov.

Kin selection is a cornerstone of modern evolutionary theory and provides arguably the most powerful set of tools we have for the theoretical investigation of adaptation under long-term evolution when genetically related individuals interact. Whilst it is most commonly used in models on interactions between organisms, we have argued that kin selection can be an equally powerful tool in understanding the evolutionary consequences of interactions at the gametic level. Interactions between related gametes are ubiquitous and do not require spatial structure in the adult population if fertilisation is internal or gametes otherwise do not fully disperse after release.

We have discussed how kin selection relates to some central themes in gamete evolution and interactions: sperm competition, sperm limitation, local sperm competition, and sex allocation. All of these themes converge in the model we have presented. We draw two main conclusions from the model.

First, neither limited dispersal nor costly dispersal influence the outcome of sex ratio evolution in simultaneous hermaphrodites. Whilst the evolution of sex allocation in simultaneous hermaphrodites can be characterised in terms of ‘local sperm competition’ (LSC), i.e., sperm competition between related sperm, our results show that the effect of relatedness should be taken with caution, just as it should in models of sex allocation under ‘local mate competition’ (LMC) with separate sexes: limited dispersal is perhaps the simplest way for relatedness to arise between interacting individuals, and likely very common in nature, yet dispersal rate and dispersal cost do not affect sex allocation in our models.

Second, sperm limitation does influence the evolution of sex allocation in simultaneous hermaphrodites. This may be particularly important for understanding sex allocation evolution in broadcast spawning or spermcasting marine hermaphrodites, where sperm limitation is a real risk. We show examples of sex allocation outcomes with a specific, well-known mathematical formulation of sperm limitation in Figure 1. Under sperm-limited conditions the influence of mating/spawning group size on the sex ratio is significantly reduced, and the group size effect disappears entirely under extremely sperm limited conditions. We should therefore not necessarily expect simultaneous hermaphrodites under sperm-limited conditions to conform to the predictions of classic sex ratio models.

Whilst there are novel aspects to the model, there are also several clear parallels to earlier models. First, classic theoretical results for sex allocation in simultaneous hermaphrodites (e.g., Charnov, 1980; see Charnov, 1982; West, 2009 for reviews) are contained as a special case. Second, we have shown that the cancellation of increased relatedness and increased competition due to limited dispersal – widely studied in sex allocation theory for separate sexes (Bulmer, 1986; Frank, 1986; Taylor, 1988b) – applies equally to simultaneous hermaphrodites, with implications for the ‘local sperm competition’ effect as described above. Third, analogous to theoretical studies of sex allocation under mate limitation in organisms with separate sexes (e.g., Nagelkerke and Hardy, 1994; Lehtonen and Schwanz, 2018), we have shown that sperm limitation has a significant effect on sex allocation: strong sperm limitation is compensated by increased allocation into male function. The novel aspect of this model is the combination of phenomena it addresses. For example, whilst previous models have explored the analogous effect of mate limitation on sex allocation, we are not aware of any model analysing the effect of mate limitation or sperm limitation under limited dispersal, and therefore under increased relatedness. Without a model examining this combination, the possibility of an interaction between these factors remains. Superficially, the importance of sperm limitation as an evolutionary pressure is expected to increase with increasing relatedness (Lehtonen, 2016). However, our model shows that whilst sperm limitation does strongly influence sex allocation, here too the cancellation effect holds.

What do our results say about the importance of considering kin selection in such models? Playing the devil’s advocate, one could argue that we have shown kin selection to be unimportant after all, at least in this model. However, our view is the opposite: the scenario we have modelled would be extremely difficult to analyse and understand without carefully considering kin selection effects. Whilst the final result might suggest that relatedness plays no role in the evolution of sex allocation, we would not know this without a formal model, and might even be misled to predict that limited, costly dispersal should (via increased relatedness) modulate sex allocation. Kin selection does play a central role in the model: its various effects at different stages of the life cycle happen to cancel each other exactly, but seeing this requires computing those effects, as well as the relatedness coefficient that arises due to the combined effect of limited patch size, limited dispersal, and mortality during dispersal. All three factors appear in the relatedness coefficient (equation 9), but only the first appears in the equilibrium sex ratio.

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

LM and JL contributed to the literature review, to planning the article, and writing the manuscript. JL contributed the mathematical model. All authors contributed to the article and approved the submitted version.

This research was funded by the Academy of Finland (grant number 340130, awarded to JL).

We would like to thank Jeremy Van Cleve and John Fitzpatrick for very helpful and constructive comments.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.