Abstract
We used computer simulations of growth, mating and death of cephalopods and fishes to explore the effect of different life-history strategies on the relative prevalence of alternative male mating strategies. Specifically, we investigated the consequences of single or multiple matings per lifetime, mating strategy switching, cannibalism, resource stochasticity, and altruism towards relatives. We found that a combination of single (semelparous) matings, cannibalism and an absence of mating strategy changes in one lifetime led to a more strictly partitioned parameter space, with a reduced region where the two mating strategies co-exist in similar numbers. Explicitly including Hamilton’s rule in simulations of the social system of a Cichlid led to an increase of dominant males, at the expense of both sneakers and dwarf males (“super-sneakers”). Our predictions provide general bounds on the viable ratios of alternative male mating strategies with different life-histories, and under possibly rapidly changing ecological situations.
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig1_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig2_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig3_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig4_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig5_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig6_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig7_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig8_HTML.png)
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11538-024-01330-z/MediaObjects/11538_2024_1330_Fig9_HTML.png)
Similar content being viewed by others
References
Allen L (2010) An introduction to stochastic processes with applications to biology, 2nd edn. Chapman and Hall
Alonzo SH, Warner RR (2000a) Allocation to mate guarding or increased sperm production in a Mediterranean wrasse. Am Nat 156(3):266–275
Alonzo SH, Warner RR (2000b) Dynamic games and field experiments examining intra-and intersexual conflict: explaining counterintuitive mating behavior in a Mediterranean wrasse Symphodus ocellatus. Behav Ecol 11(1):56–70
Apostólico LH, Amoroso RodriguezMarian JE (2020) Cephalopod mating systems as models for the study of sexual selection and alternative reproductive tactics: a review. Bull Mar Sci 96(2):375–398
Bourke AF (2014) Hamilton’s rule and the causes of social evolution. Philos Trans R Soc B Biol Sci 369(1642):20130362
Brawand D, Wagner CE, Li YI, Malinsky M, Keller I, Fan S, Di Palma F (2014) The genomic substrate for adaptive radiation in African cichlid fish. Nature 513(7518):375–381
Brown C, Garwood MP, Williamson JE (2012) It pays to cheat: tactical deception in a cephalopod social signalling system. Biol Let 8(5):729–732
Drilling CC, Grober MS (2005) An initial description of alternative male reproductive phenotypes in the bluebanded goby, Lythrypnus dalli (Teleostei, Gobiidae). Environ Biol Fishes 72(4):361–372
Engqvist L, Taborsky M (2016) The evolution of genetic and conditional alternative reproductive tactics. Proc R Soc B Biol Sci 283(1825):20152945
Froese R, Pauly D (Eds) (2021) FishBase. Family Labridae—Wrasses. World Wide Web electronic publication. https://www.fishbase.se/summary/FamilySummary.php?ID=362
Hall K, Hanlon R (2002) Principal features of the mating system of a large spawning aggregation of the giant Australian cuttlefish Sepia apama (Mollusca: Cephalopoda). Mar Biol 140(3):533–545
Hamilton WD (1963) The evolution of altruistic behavior. Am Nat 97(896):354–356
Hamilton WD (1972) Altruism and related phenomena, mainly in social insects. Annu Rev Ecol Syst 3(1):193–232
Hanlon RT, Messenger JB (2018) Cephalopod behaviour. Cambridge University Press
Hirohashi N, Iwata Y (2013) The different types of sperm morphology and behavior within a single species: Why do sperm of squid sneaker males form a cluster? Commun Integr Biol 6(6):e26729
Horth NA, Dodson JJ (2004) Influence of individual body size and variable thresholds on the incidence of a sneaker male reproductive tactic in Atlantic salmon. Evolution 58(1):136–144
Iba CM (2010) Cannibalism in cephalopods. Rev Fish Biol Fish 20:123–136
Ibáñez CM, Pérez-Álvarez J, Catalán J, Carrasco SA, Pardo-Gandarillas MC, Rezende EL (2019) Sexual selection and the evolution of male reproductive traits in benthic octopuses. Front Physiol 10:1238
Marentette JR, Fitzpatrick JL, Berger RG, Balshine S (2009) Multiple male reproductive morphs in the invasive round goby (Apollonia melanostoma). J Great Lakes Res 35(2):302–308
Marian JE, Apostólico LH, Chiao CC, Hanlon RT, Hirohashi N, Iwata Y, Shaw PW (2019) Male alternative reproductive tactics and associated evolution of anatomical characteristics in loliginid squid. Front Physiol 10:1281
Naud MJ, Shaw PW, Hanlon RT, Havenhand JN (2005) Evidence for biased use of sperm sources in wild female giant cuttlefish (Sepia apama). Proc R Soc B Biol Sci 272(1567):1047–1051
Norman MD, Finn J, Tregenza T (1999) Female impersonation as an alternative reproductive strategy in giant cuttlefish. Proc R Soc Lond Ser B Biol Sci 266(1426):1347–1349
O'Dor RK (1998) Squid life-history strategies. FAO Fisheries Technical Paper, pp 233–254
Ota K, Kohda M (2006) Description of alternative male reproductive tactics in a shell-brooding cichlid, Telmatochromis vittatus, Lake Tanganyika. J Ethol 24(1):9–15
Reichard M, Le Comber SC, Smith C (2007) Sneaking from a female perspective. Anim Behav 74(4):679–688
Ros AF, Bouton N, Santos RS, Oliveira RF (2006) Alternative male reproductive tactics and the immunocompetence handicap in the Azorean rock-pool blenny, Parablennius parvicornis. Proc R Soc B Biol Sci 273(1589):901–909
Sato N (2021) A review of sperm storage methods and post-copulatory sexual selection in the Cephalopoda. Biol J Lin Soc 134(2):285–302
Sato T, Hirose M, Taborsky M, Kimura S (2004) Size-dependent male alternative reproductive tactics in the shell-brooding cichlid fish Lamprologus callipterus in Lake Tanganyika. Ethology 110(1):49–62
Schiettekatte N, Brandl S, Casey J (2019) Fishualize: color palettes based on fish species. R package version 0.1.0. https://CRAN.R-project.org/package=fishualize
Schütz D, Taborsky M (2000) Giant males or dwarf females: What determines the extreme sexual size dimorphism in Lamprologus callipterus? J Fish Biol 57(5):1254–1265
Stiefel KM (2013) Why are there no eusocial fishes? Biol Theory 7(3):204–210
Stiefel KM (2014) Dishonest signaling in vertebrate eusociality. Biol Theory 9(3):325–330
Stiver KA, Harris RM, Townsend JP, Hofmann HA, Alonzo SH (2015) Neural gene expression profiles and androgen levels underlie alternative reproductive tactics in the ocellated wrasse Symphodus ocellatus. Ethology 121(2):152–167
Suzuki S, Toguchi K, Makino Y, Kuwamura T, Nakashima Y, Karino K (2008) Group spawning results from the streaking of small males into a sneaking pair: male alternative reproductive tactics in the threespot wrasse Halichoeres trimaculatus. J Ethol 26(3):397–404
Taborsky M (1994) Sneakers, satellites, and helpers: parasitic and cooperative behavior in fish reproduction. Adv Study Behav 23(1):e100
Taborsky M (1997) Bourgeois and parasitic tactics: do we need collective, functional terms for alternative reproductive behaviours? Behav Ecol Sociobiol 41(5):361–362
Taborsky M (2001) The evolution of bourgeois, parasitic, and cooperative reproductive behaviors in fishes. J Hered 92(2):100–110
Tsuchiya K, Uzu T (1997) Sneaker male in octopus. Venus (Jpn J Malacol) 56(2):177–181
Turner GF (2007) Adaptive radiation of cichlid fish. Curr Biol 17(19):827–831
Van den Berghe EP, Wernerus F, Warner RR (1989) Female choice and the mating cost of peripheral males. Anim Behav 38(5):875–884
Wagner CE, Harmon LJ, Seehausen O (2012) Ecological opportunity and sexual selection together predict adaptive radiation. Nature 487(7407):366–369
Warner RR (1982) Mating systems, sex change and sexual demography in the rainbow wrasse, Thalassoma lucasanum. Copeia 1982:653–661
Warner RR, Hoffman SG (1980) Local population size as a determinant of mating system and sexual composition in two tropical marine fishes (Thalassoma spp.). Evolution 34:508–518
Warner RR, Lejeune P (1985) Sex change limited by paternal care: a test using four Mediterranean labrid fishes, genus Symphodus. Mar Biol 87(1):89–99
Warner RR, Swearer SE (1991) Social control of sex change in the bluehead wrasse, Thalassoma bifasciatum (Pisces: Labridae). Biol Bull 181(2):199–204
Wilbur HM, Rudolf VHW (2006) Life-history evolution in uncertain environments: bet hedging in time. Am Nat 168(3):398–411
WirtzOcana S, Meidl P, Bonfils D, Taborsky M (2014) Y-linked Mendelian inheritance of giant and dwarf male morphs in shell-brooding cichlids. Proc R Soc B Biol Sci 281(1794):20140253
Zylinski S, How MJ, Osorio D, Hanlon RT, Marshall NJ (2011) To be seen or to hide: visual characteristics of body patterns for camouflage and communication in the Australian giant cuttlefish Sepia apama. Am Nat 177(5):681–690
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix: Equations and Parameters of the Models
Appendix: Equations and Parameters of the Models
Typically models simulating the evolution of populations of animals are based on dynamical systems. In such a dynamical system a state vector, x, represent the fractions of a population. If the population is divided into n unique states, we can view x as lying in a phase space \({\varvec{x}}\in X\subset {R}^{n}\), and the population can be seen as evolving according to a map, \(\varphi\):
Such a mapping can either be discrete (i.e., \(t\in {\mathbb{Z}}\)) or continuous (i.e., \(t\in {\mathbb{R}}\)). Additionally, it may or may not be smooth over time. In other words, it is conceivable that the population will evolve in different ways depending on separate periods or “seasons”. In our model, we simulate a mating season alternating with a growth season. The evolution of the state during the growth season will depend on how the state evolved during the previous mating season.
The model we used is a discrete, piecewise-smooth dynamical systems. Each timestep represents 1 day, and at each step the population is updated according to a time dependent mapping. With \(\mu\) as a set of parameters, the master equation of our model is:
Separating the growth and mating season in our model is implemented as:
where P and Q are matrices, each depending on a set of parameters, \(\upmu\). Q also depends on the values taken on by the state vector, x, during the previous mating season. The evolution of a state vector according to a transition matrix such as this is generally referred to as a Markov chain.
The parameter range chosen for figures is intended to demonstrate the full portrait of behavior, not guarantee convergence. In general, it is a classic result from probability theory that Markhov chain converges to a unique stationary distribution (equilibrium population) if the transition matrix is irreducible (Allen 2010); in this case, that transition matrix would be P60 + Q300. However, there is a time delay in our model (to let prior mating influence future birth), so instead we demonstrate convergence with longtime simulations.
Figures including results iterate this by applying \(P\) 60 times and then Q 300 times for a total of 50 cycles. Figure 10 demonstrates that 50 years is sufficient for producing the same behavior as is present in 100 year and 1000 year simulations.
Base Model: Wrasses
To model the alternative male mating strategies for wrasses, we used a piecewise system of Markov chains (Fig. 1). On any given day, a fish can move from state x to state y with probability Px→y. A one by seven state vector x[t] is initialized as follows:
The population at t = 0 is assumed to be uniformly distributed across the high and low energy states for the two mating strategies (dominant/sneaker). The population is simulated with two transition matrices:
P is applied to x[t] for the first 60 days (during the mating season), and then Q is applied to x[t] for the following 300 days (during the growth season).
The values of x[t][i] take on real numbers between zero and one representing the proportion of the population that any given state represents. Hence the following holds true for all t:
During the mating season:
During the growth season:
The element in the location (i,j) of a transition matrix is the probability that a subject in state i will transition towards state j in 1 day. The \({i}^{th}\) state is defined as the state corresponding to the \({i}^{th}\) entry of x[t]. Observe that the placement of zeros is consistent with the diagram given by Fig. 1. Making choices for these parameter values requires a consideration of both (1) how the previous mating season played out and (2) to what extent predation and food are readily present in the environment.
To incorporate (1) into the model, we utilize \({Q}_{d\to ls}\) and \({Q}_{d\to ld}\). More accurately, we devised a deterministic method for pulling subsets of the population out of the dead state and into the two low energy states. This re-drawing of individuals from a pool of the dead is a mathematical approach reflects the birth of new fish during the growth season resulting from mating in previous months, devoid of biological interpretation. We used the following method for choosing \({Q}_{d\to ls}\) and \({Q}_{d\to ld}\) at the beginning of the growth season for the \({k}^{th}\) year of iteration:
The likelihood that a new sneaker is born depends explicitly on the relative fraction of matings by sneakers in the previous cycle relative to dominant males. \(\updelta \in \left(\text{0,1}\right)\) reflects volatility in the population, or the proportion of the dead pool distributed towards sneakers and dominants during the subsequent growth season. We choose \(\delta = 0.2\).
The remaining transition probabilities were chosen according to the following relationships:
Here z corresponds to an arbitrary state. More specifically, let \(\alpha\) correspond to food availability and \(\beta\) predation. The transition probabilities are as follows:
Low Energy Sneaker Transitions:
High Energy Sneaker Transitions:
Low Energy Dominant Transitions:
High Energy Dominant Transitions:
For the two-dimensional parameter sweeps (Figs. 4 and 5), \(\alpha\) was swept from 1.20 to 2.00, \(\beta\) from 10−5 to 0.02. This choice magnifies the curve of change in ESS. These have been re-scaled to a 1–100 scale for the purpose of figures.
The specification of these proportionality constants was done with the following principles in mind:
-
The sneaker male is more likely to increase its energy level in a given day than the dominant male.
-
Though the dominant male is slower to evolve, it is less likely to subsequently loose energy.
-
The sneaker male is more likely to die than the dominant.
-
The sneaker male is less likely to transition from the high energy state to the mating state than the dominant male.
Continuing to fill in our matrix, we must consider the values along the diagonals, which describe the probabilities that a subject does not change its state. There is no immediately obvious biological intuition for these values. The columns of each transition matrix need to sum to one. Hence, we will choose each x and y such that this property is necessarily true.
The last parameter to determine is \({Q}_{hs\to ld}\): the frequency for which a male mating strategy change occurs. For this model, we take that to be a constant. The mating strategy change is known to depend on environmental variables in some species, however this dependence is not included here.
In the simulations the updates of the values of the matrix were iterated for 50 years. Note that a simulated year does not correspond directly to a chronological year, since the simulations were set up to reach an equilibrium at the fastest possible rate and aim to reproduce the equilibrium, but not the path to equilibrium of biological evolution. In Fig. 10 we compare the state of the simulations after 10, 50 and 1000 simulated years, and observe that the sharpness of the transition between sneakers and dominant males, but not the qualitative shape of the result changes.
Cephalopods
We model cuttlefish differently from wrasses by instituting four characteristic changes. First, the sneaker to dominant transition is removed as cuttlefish are not observed to change strategies. Next, cuttlefish are observed to be semelparous, so the state of the population undergoes a mass death directly following every mating season. This is executed by resetting the state vector, X[t], to be heavily biased towards the dead state once per year. The model does not reset though, given that the information from the previous mating cycle will still be encoded in the parameters \({Q}_{d\to ls}\) and \({Q}_{d\to ld}.\)
Thirdly, a cannibalism mechanism is introduced. This is executed by letting \({Q}_{ls\to d}\) and \({Q}_{hs\to d}\) be dependent on the total population of cuttlefish playing the dominant strategy at any given time. More specifically, at each time step, t, these two transition probabilities are redefined in the following way:
Additionally, the values of \({y}_{1}\) and \({y}_{2}\) are updated appropriately at each time so that \(\left(\star \right)\) holds. Finally, the growth rates of cuttlefish are higher than those of fishes. Subsequently, we increased the values of \({Q}_{ls\to hs}\) and \({Q}_{ld\to hd}\) by a factor of two.
Figure 3 includes parameter sweeps over varied predation that result from each of these four changes applied individually to the basal wrasse model. The normalized food availability of \(\alpha = 1.52\).
Lamprologus callipterus
Lastly, we shall consider the case of L. callipterus. A third mating strategy, the dwarf male, is introduced as a genetically determined sneaker. We do this by introducing three new dimensions for our state vector and transition matrices:
And during the mating season:
During the growth season:
We also included a link between the dominant and dwarf strategies. Specifically, the dwarf males are assumed to be genetically related to the dominant males. This leads to an incentive for altruistic behavior for the dominant males.
To incorporate this, we introduced Hamilton’s rule by assuming a mean relatedness between any two male members of the breeding population. The dominant Lamprologus individuals receive a benefit for their altruism that is proportional to the total number of dwarf male in the population. This benefit will materialize as an increase to their probability of transitioning to the high energy state. At every time interval we performed the following update:
The remaining transition probabilities are given as follows.
Low Energy Dwarf Transitions:
High Energy Dwarf Transitions:
Given these values, one could fill in the transition matrices P and Q and simulate the population. Figure 7 includes the resulting ESS incorporating parameter sweeps over predator presence and food availability (with and without the altruism bonus). In 7a and 7b, predation is fixed at \(\beta = 0.001\); in 7c and 7d, food availability is fixed at \(\alpha = 1.0.\) In both cases, the other parameter was swept over the same aforementioned ranged as for the wrasses.
Stochastic Considerations
We recognize that food availability and predator presence are not strictly fixed parameters in practice, but instead quantities that fluctuate with time around some initialized state. This is realized in our model but introducing stochasticity in the low to high and high to low energy transitions.
Consider the baseline wrasse model. We can introduce multiplicative noise with a stochastic variable, \(\xi\), which acts as a multiplier on the low to high energy transitions and a divider on the high to low transitions. \(\xi\)> 1 corresponds to a relatively high concentration of food and relatively low predator presence, specifically due to environmental variability. \(\xi\)< 1 is interpreted as the natural inverse. \(\xi\) is chosen via the following mechanism: at the beginning of a new season (mating or growth), reset \(\xi\)= 1. For each following day, update the variable according to the following equation:
That is, a new value of \(\xi\) is chosen from a random normal distribution with mean \(\xi\) and variance \({\sigma }^{2}\). In other words, \(\xi\) takes a random walk during each season before being reset. During a mating season, transition probabilities are updated in the following way each day:
Similarly, during the growth season:
We can then run simulations to observe that the sneaker or dominant strategy may be present in parameter regimes despite being unstable in the deterministic sense.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Landsittel, J.A., Ermentrout, G.B. & Stiefel, K.M. A Theoretical Comparison of Alternative Male Mating Strategies in Cephalopods and Fishes. Bull Math Biol 86, 98 (2024). https://doi.org/10.1007/s11538-024-01330-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11538-024-01330-z