Population density and group size effects on reproductive behavior in a simultaneous hermaphrodite
© Sprenger et al; licensee BioMed Central Ltd. 2011
Received: 26 October 2010
Accepted: 18 April 2011
Published: 18 April 2011
Despite growing evidence that population dynamic processes can have substantial effects on mating system evolution, little is known about their effect on mating rates in simultaneous hermaphrodites. According to theory, mating rate is expected to increase with mate availability because mating activity is primarily controlled by the male sexual function. A different scenario appears plausible in the hermaphroditic opisthobranch Chelidonura sandrana. Here, field mating rates are close to the female fitness optimum, suggesting that mating activity remains unresponsive to variation in mate availability.
Applying an experimental design that aims at independent experimental manipulation of density and social group size, we find substantial increases in mate encounter rate with both factors, but no statistically detectable effects on mating rate in C. sandrana. Instead, mating rate remained close to the earlier determined female fitness optimum.
We demonstrate that mating rate in C. sandrana is largely unresponsive to variation in mate availability and is maintained close to the female fitness optimum. These findings challenge the prevailing notion of male driven mating rates in simultaneous hermaphrodites and call for complementary investigations of mating rate effects on fitness through the male sexual function.
Keywordspopulation density mating group size optimal mating rate sexual conflict simultaneous hermaphrodite
There is increasing awareness that ecological processes such as changes in population demography can substantially affect the evolution of mating systems (reviewed in ). For example, theoretical and empirical work shows that female selectivity generally increases at higher population densities, resulting in stronger sexual selection on males [2–8]. Yet, few studies to date have investigated the effect of variation in population density on realized mating rates . In water striders and dung flies, two systems where males are capable of enforcing female matings, empirical work showed that mating rate increases with population density as expected [9, 10]. Alternatively, in systems where females have full control over mating probability and mating partners are chosen randomly, theoretical work predicts that mating rates should remain unresponsive to variation in population density . This is because, in the absence of male harassment, females can maintain mating rates close to their reproductive optimum by modulating the number of matings with males accordingly. Yet, empirical investigations of this scenario are lacking to date.
In simultaneous hermaphrodites, the situation is less straightforward than in separate sex species. Mating rates of both sexual functions are often directly linked via reciprocal mating, where individuals donate and receive sperm during a single copulatory bout [12, 13]. Theory suggests that, in simultaneous hermaphrodites, mating rates are largely driven by the male sex function for two, reasons. First, similar to many separate sex species and in accordance with Bateman's Principle , the male function is thought to show a steeper fitness gain with increasing mating success than the female function ([15, 16] but see ). Second, individual costs associated with matings in the female function can be compensated by a gain in paternity , rendering hermaphrodites likely to accept even costly female matings. If mating rate is indeed controlled by the male function, one would expect that mating rate increases at higher population densities because higher densities should increase mate encounter rates and thus mating opportunities. In support of this prediction a recent study on the hermaphroditic flatworm Macrostomum lignano showed that mating activity increased with partner availability . Similarly, the freshwater snail Lymnaea stagnalis was found to receive more inseminations when kept in groups of eight compared to groups of four or pairs .
A different scenario being concordant with the predictions by Härdling & Kaitala  appears plausible for the simultaneously hermaphroditic opisthobranch Chelidonura sandrana. In this non-selfing internal fertilizer, multiple matings are frequent in both sexual functions . Experimental work showed that maternal per-egg investment depends on the number of different mating partners , being maximized at an intermediate mating rate of 2-3 matings per day . At this mating rate, fecundity costs due to multiple matings appear to be offset through reduced embryo mortality [23, 24]. Field surveys further revealed that this intermediate daily mating rate established in the laboratory closely matches realized field mating rates . All these findings suggest that the detected mating rate is close to the female fitness optimum. Moreover, because mate encounter rates typically exceed the number of copulations in this species (personal observations) the natural mating rate may be maintained close to the female fitness optimum independent of mate availability.
Against this background our main objective was to experimentally test for mate availability effects on mating rate in C. sandrana. Given that frequent matings in the female function result in reduced fecundity  and that comparable studies are rare in simultaneous hermaphrodites, we further tested for mate availability effects on female fecundity. Because mate availability can be both a function of social group size (the total number of interacting individuals) and density (the number of individuals per area) we applied an experimental approach that, in theory, is capable of differentiating between both effects. We hypothesized that (i) variation in mate availability has little effect on mating activity and that (ii) the realized mating rate will be close to the earlier determined fitness optimum whenever mate encounters are sufficiently high to maintain this optimal mating rate. Because our predictions match those of the theoretical work by Härdling and Kaitala  developed for separate sex species (see above), we discuss our findings accordingly.
General Linear Mixed Model results for the effects of group size (nominal factor) and density (continuous factor) on mating behaviour and fitness components.
Type III SS
Mate encounter rate (per individual & day)
Full model (linear)
log Mating rate (per individual & day)
Full model (linear)
Average egg mass weight
Full model (quadratic)
Total egg mass weight per individual
Full model (linear)
N egg masses per individual
Full model (linear)
We found that mate encounter rates and thus mating opportunities increased with both group size and density. Despite this, the average individual mating rate remained largely constant and close to the earlier detected fitness optimum of the female function in C. sandrana . Importantly, the observed increase in mate encounters across densities by ~100% did not transform in a statistically detectable increase in mating rate, where the maximum change rendered undetectable due to noise in the data does not exceed a 15% increase across densities (upper 95% CI bound of regression slope). This at best weak effect can largely be attributed to reduced mating activity at the lowest mate encounter rates, indicating that under such conditions optimal mating rates are difficult to maintain, whereas mating rates remain stable at higher densities where mate encounters are sufficiently frequent.
Although the existing variation in mating rate implies strong between-individual differences in mating activity, the documented absence of significant group size and density effects clearly contradicts the general notion of a tight positive relationship between mate availability and average mating rate as previously reported from other simultaneously hermaphroditic [20, 19] and separate sex [9, 10] species. Instead, our findings indicate that the average daily mating rate in C. sandrana is largely independent of mate availability and close to the female fitness optimum, as long as the latter can be realized. These findings challenge the general notion of male driven mating rates in simultaneous hermaphrodites and suggest three alternative scenarios. First, conflict over mating rate may be minute or even absent if both sexual functions share similar mating rate optima, i.e. when the here detected average mating rate is close to the female and male mating optimum. In this case, no sex function would ultimately control matings, but realized mating rates result from mutual interest. Such low male mating rate optima are plausible if remating in the male sexual function generates accelerating costs, e.g. via sperm digestion by the sperm recipient . However, previous work in C. sandrana contradicts this scenario of accelerating male costs because animals are easily capable of donating sperm up to eight times within 9 h . Second, realized mating rates may be intermediate between divergent male and female optima while providing roughly similar fitness returns for both sexual functions (representing inherent sexual antagonism;). Under this scenario, mating rate represents the balanced result of opposing sex-specific interests. Our previous finding that realized field mating rates are close to the female fitness optimum  do not lend support to this scenario. Third, assuming that the male mating rate optimum exceeds that of the female function (see  for conforming data in freshwater planarians), then mating rate should be largely controlled by the female function as it is close to the female fitness optimum of C. sandrana . This last scenario is currently most concordant with the available data but detailed information on mating rate effects on male fitness are clearly needed to substantiate this conclusion.
Interestingly, our findings conform to recent theoretical predictions for separate sex species by Härdling & Kaitala . Their model predicts that females evolve constant mating rates that are largely independent of mate availability if the following three central assumptions are met. First, female fitness solely depends on the number of different male partners, not on the individual quality of certain males. Thereby the model explicitly excludes female choosiness for high quality males as an adaptive trait. Second, mating probability is under full female control, excluding systems in which females accept matings due to sexual harassment. Third, female fitness is maximized at an intermediate mating rate at which the fecundity benefits and the mortality costs of multiple mating balance. All these three key model assumptions appear to be met in our study system C. sandrana. First, polyandry-mediated benefits primarily depend on the number of different mating partners, not on the identity or quality of mates . The production of multiply sired egg masses seems to represent a diversifying genetic bet-hedging strategy that increases the probability of offspring survival under fluctuating environmental conditions . Female fitness is thus a direct function of mating frequency irrespective of male quality. Second, male acting individuals do not show any harassment or otherwise manipulative behavior of female mating activity [24, 28]. This is further supported by the observation that realized mating rates both under laboratory and field conditions  are close to the female fitness optimum. Third, female fitness is maximized at an intermediate mating rate, where fecundity costs of multiple matings appear to be offset through increased offspring viability .
We further found that fitness measures for the female function were independent of variation in group size and solely a function of density, with fecundity being maximized at an intermediate density. Because mating rate remained constant across all these treatments, density-dependent differences in fecundity are unlikely to represent mere differences in allosperm availability. Instead, reduced fecundity at low densities might represent costs due to increased mate searching. Interestingly, previous work in C. sandrana showed that repeated matings with the same male result in decreased fecundity relative to repeated matings with different males . Because the probability of mating with the same partner is higher at lower densities, the same currently unidentified mechanism might contribute to reduced offspring production at low densities. In this context it remains puzzling, however, why group size, i.e. the actual number of available mating partners, had no comparable effect on fecundity in the present study, because the likelihood of copulating repeatedly with the same partner equals 1 in our pair treatment. Perhaps, differences in social group size do not reflect differences in mating group size, with the latter being primarily affected by the actual distance between animals (i.e. density).
Decreasing fecundity at the highest densities may be the result of a shift in sex allocation towards the male function. Although our data provide no direct measure of sex allocation, various studies confirmed rapid strategic reallocation of resources towards a hermaphrodite's male function with increasing social group size (reviewed in ). For example, in Macrostomum lignano, testis size tends to increase with social group size while ovary size significantly decreases . Although constancy in average daily mating rate may suggest little variation in mean sex allocation in C. sandrana, resource allocation may have been adjusted in response to mate availability rather than to mating rate. Unlike in M. lignano, however, sex allocation cannot be measured in vivo in C. sandrana, rendering definite conclusions on sex allocation adjustment difficult. Possible components of resource re-allocation that should be explicitly addressed in follow-up studies include the amount of sperm or seminal fluid transferred, the composition of these seminal fluids (e.g. with respect to manipulative substances, ), or energetic investment in other male components such as precopulatory interactions, all of which may contribute to compromised female fecundity at high densities. Finally, declining fecundity at high densities may have been caused by waste products or certain metabolites, which have been suggested to act as egg laying inhibitory substances in pulmonate gastropods . Although water was regularly exchanged and experimental containers were frequently cleaned to reduce such effects, it is possible that inhibitory substances accumulated at higher densities or were actively produced by animals in response to higher mate encounter rates.
The here documented differential effects of group size and density on reproductive behaviour imply that studies on mate availability effects need to carefully disentangle both factors by applying an appropriate experimental design. However, our findings also indicate that a non-random distribution of animals over the available space may make a clear differentiation between the relative contributions of group size and density on mating opportunities difficult, even though the applied experimental design allows, in theory, to do so. For example, while mate encounter rate would be expected to increase at higher densities as the distance between individuals is reduced , mate encounter rates in our study also increased with group size independent of density. Given that mating aggregations are prevalent in many animal systems, future experimental work that directly manipulates spatial distribution is needed in order to shed light on the detailed relationship between aggregation behavior and reproductive behavior.
Contrary to classic predictions by theory, we here show for the first time that mating rate is largely unresponsive to variation in mate availability in a simultaneous hermaphrodite. With mating rates being close to the female fitness optimum, our findings challenge the prevailing notion of male driven mating rates in simultaneous hermaphrodites and call for complementary investigations of mating rate effects on fitness through the male sexual function.
Chelidonura sandrana Rudman 1973 (Cephalaspidea, Aglajidae) inhabits tropical marine shallow water sand flats. Copulations are unilateral, but sexual roles are usually alternated during a copulatory bout, resulting in reciprocal sperm exchange (details in ). Following mating, egg masses containing hundreds of singly encapsulated eggs are deposited on the substrate every few days . Embryos develop into feeding veliger larvae with a planktonic period of unknown length. Animals are sexually active throughout the Australian summer (personal observation) and occur in patchy aggregations with up to dozens of individuals/m2.
Sampling and maintenance
Individuals were collected in 0.5-8 m depth around Lizard Island, Queensland, Australia, in December 2007 and January 2008 (permit G05/15308·1 of the Great Barrier Reef Marine Park Authority). On the day of collection, body wet weight was measured to the nearest mg. Animals were then randomly distributed among the experimental groups and kept at 26°C water temperature and a natural diurnal cycle. Water was stirred twice per day and experimental containers were thoroughly cleaned and refilled with fresh seawater every fourth day to minimize accumulation of waste products. Every other day animals were fed with their natural food source, the flatworm Wulguru cuspidata Winsor, 1988 ad libitum.
Behavioral and fitness measures
Each container was observed 3 times per day (8:00, 12:00, 16:00) for 30 minutes and the following parameters were recorded: (i) Mate encounter rate: Because behavioral activity ceases at night (personal observation) we assumed a daily activity period of 12 h, i.e. 8 times the daily observation period of 1.5 h. Daily mate encounter rate per individual thus calculates as 8*((N observed mate encounters per container * 2 individuals/7 observation days)/group size). (ii) Mating rate: an estimate of the average number of copulations for each individual per day. Given that this species typically alternates sexual roles within a copulatory bout  the obtained mating rates can be easily transformed into estimated numbers of matings in both sex roles. Considering a daily activity period of 12 h (see above), individual daily mating rate calculates as 8*((N observed matings per container * 2 individuals/7 observation days)/group size). Matings that were already in progress at the start of observations or that continued beyond the observation period were weighted by 0.5 to avoid artificial inflation of mating rate estimates. Assuming that an individual would rarely be involved in more than one copulation per 30' observation interval, with average copulation duration ~10 min ([28, 34]), our observation scheme allowed detecting individual daily mating rates ranging between zero and 24.
For fitness measurements, egg masses were collected daily from each container and the following parameters were recorded: (i) Mean egg mass weight per container: Total egg mass weight per container divided by the number of egg masses; (ii) Total egg mass weight per capita: Total egg mass weight per container divided by the number of individuals; (iii) N egg masses per individual: egg mass count per container divided by the number of individuals.
We used a General Linear Mixed Model to test for the effects of group size (nominal fixed factor) and density (continuous fixed factor) on female reproductive parameters as dependent variables. All dependent variables were entered as means per experimental container, and their values were log-transformed if necessary to meet the assumptions of normality in data and residual distributions. To account for temporal effects, we included the term experimental run as a random effect. In order to assess whether density effects varied between the three group sizes, we further included the interaction group size * density as a fixed factor. The interaction term was omitted from the final models whenever P exceeded 0.25 . Larger groups typically yield less variable estimates of traits than smaller groups, therefore violating the assumption of equality in variances. To meet requirements of parametric testing, we thus weighted each dependent variable differentially for each group size . Weighting factors were 1 for group size 8, 2 for group size 4, and 4 or multiples of 4 for group size 2, until no significant differences in variance (Levene's test) were detectable between groups . In cases where visual inspection indicated a non-linear effect of density on mean egg mass weight, we added density as a second-order polynomial term and tested whether the reduced or the full model provided a better fit using partial F statistics . For near significant parameters (0.05 ≤ P ≤ 0.1) we calculated the upper bound of the 95% confidence limits for the regression slope over the complete density range (representing a 300% increase from the lowest to the highest density) to quantify the maximum undetected change across the corresponding predictor variable range. All statistical analysis were performed with JMP IN version 8·0·1 (SAS Institute Inc., Cary, NC, USA).
We are grateful to Hanna Kokko, Roger Härdling, and two anonymous referees for valuable comments on a previous draft of this manuscript, Kate Kunigelis for support during field work, and the staff at Lizard Island Research Station, Queensland, Australia, for excellent logistic support. This work was funded by a DFG grant to Nico K. Michiels and NA (Mi 482/7-4).
- Kokko H, Rankin DJ: Lonely hearts or sex in the city? Density-dependent effects in mating systems. Phil Trans R Soc Lond B. 2006, 361: 319-334. 10.1098/rstb.2005.1784.View ArticleGoogle Scholar
- Gwynne DT: Sexual selection and sexual differences in Mormon crickets (Orthoptera: Tettigoniidae, Anabrus simplex). Evolution. 1984, 38: 1011-1022. 10.2307/2408435.View ArticleGoogle Scholar
- Hubbell SP, Johnson LK: Environmental variance in lifetime mating success, mate choice, and sexual selection. Am Nat. 1987, 130: 91-112. 10.1086/284700.View ArticleGoogle Scholar
- Crowley PH, Travers SE, Linton MC, Cohn SL, Sih A, Sargent RC: Mate density, predation risk, and the seasonal sequence of mate choices: a dynamic game. Am Nat. 1991, 137: 567-596. 10.1086/285184.View ArticleGoogle Scholar
- Palokangas P, Alatalo RV, Korpimäki E: Female choice in the kestrel under different availability of mating options. Anim Behav. 1992, 43: 659-665. 10.1016/S0003-3472(05)81024-3.View ArticleGoogle Scholar
- Souroukis K, Murray AM: Female mating behvaior in the field cricket, Gryllus pennsylvanicus (Orthoptera, Gryllidae) at different operational sex ratios. J Insect Behav. 1995, 8: 269-279. 10.1007/BF01988910.View ArticleGoogle Scholar
- Kokko H, Monaghan P: Predicting the direction of sexual selection. Ecol Lett. 2001, 4: 159-165. 10.1046/j.1461-0248.2001.00212.x.View ArticleGoogle Scholar
- Kokko H, Mappes J: Sexual selection when fertilization is not guaranteed. Evolution. 2005, 59: 1876-1885.View ArticlePubMedGoogle Scholar
- Lauer MJ, Sih A, Krupa JJ: Male density, female density and inter-sexual conflict in a stream-dwelling insect. Anim Behav. 1996, 52: 929-939. 10.1006/anbe.1996.0241.View ArticleGoogle Scholar
- Martin OY, Hosken DJ: The evolution of reproductive isolation through sexual conflict. Nature. 2003, 423: 979-982. 10.1038/nature01752.View ArticlePubMedGoogle Scholar
- Härdling R, Kaitala A: The evolution of repeated mating under sexual conflict. J Evol Biol. 2005, 18: 106-115.View ArticlePubMedGoogle Scholar
- Michiels NK: Mating conflicts and sperm competition in simultaneous hermaphrodites. Sperm competition and sexual selection. Edited by: Birkhead TR, Møller AP. 1998, London: Academic Press, 219-254. full_text.View ArticleGoogle Scholar
- Anthes N: Mate choice and reproductive conflict in simultaneous hermaphrodites. Animal Behaviour: Evolution and Mechanisms. Edited by: Kappeler P. 2010, Heidelberg: Springer, 329-358. full_text.View ArticleGoogle Scholar
- Bateman AJ: Intra-sexual selection in Drosophila. Heredity. 1948, 2: 349-368. 10.1038/hdy.1948.21.View ArticlePubMedGoogle Scholar
- Charnov EL: Simultaneous hermaphroditism and sexual selection. Proc Natl Acad Sci USA. 1979, 76: 2480-2484. 10.1073/pnas.76.5.2480.View ArticlePubMedPubMed CentralGoogle Scholar
- Anthes N, David P, Auld JR, Hoffer JNA, Jarne P, Koene JM, Kokko H, Lorenzi MC, Péllisié B, Sprenger D, Staikou A, Schärer L: Bateman Gradients in Hermaphrodites: An Extended Approach to Quantify Sexual Selection. Am Nat. 2010, 176: 249-263. 10.1086/655218.View ArticlePubMedGoogle Scholar
- Leonard JL: Bateman's principle and simultaneous hermaphrodites: a paradox. Integr Compar Biol. 2005, 45: 10.1093/icb/45.5.856.Google Scholar
- Michiels NK, Koene JM: Sexual selection favors harmful mating in hermaphrodites more than in gonochorists. Integr Compar Biol. 2006, 46: 473-480. 10.1093/icb/icj043.View ArticleGoogle Scholar
- Janicke T, Schärer L: Determinants of mating and sperm-transfer success in a simultaneous hermaphrodite. J Evol Biol. 2009, 22: 405-415. 10.1111/j.1420-9101.2008.01660.x.View ArticlePubMedGoogle Scholar
- Koene JM, Ter Maat A: Coolidge effect in pond snails: male motivation in a simultaneous hermaphrodite. BMC Evolutionary Biology. 2007, 7: 212-10.1186/1471-2148-7-212.View ArticlePubMedPubMed CentralGoogle Scholar
- Anthes N, Putz A, Michiels NK: Hermaphrodite sex role preferences: the role of partner body size, mating history and female fitness in the sea slug Chelidonura sandrana. Behav Ecol Sociobiol. 2006, 60: 359-367. 10.1007/s00265-006-0173-5.View ArticleGoogle Scholar
- Sprenger D, Anthes N, Michiels NK: Multiple mating affects offspring size in the opisthobranch Chelidonura sandrana. Mar Biol. 2008, 153: 891-897. 10.1007/s00227-007-0861-3.View ArticleGoogle Scholar
- Sprenger D, Faber J, Michiels NK, Anthes N: Natural female mating rate maximizes hatchling size in a marine invertebrate. J Anim Ecol. 2008, 77: 696-701. 10.1111/j.1365-2656.2008.01376.x.View ArticlePubMedGoogle Scholar
- Sprenger D, Lange R, Michiels NK, Anthes N: Sources of phenotypic variation in offspring size and performance in an internally fertilizing marine hermaphrodite. Evol Ecol. 2010, 24: 185-194. 10.1007/s10682-009-9300-x.View ArticleGoogle Scholar
- Greeff JM, Michiels NK: Sperm digestion and reciprocal sperm transfer can drive hermaphrodite sex allocation to equality. Am Nat. 1999, 153: 421-430. 10.1086/303184.View ArticleGoogle Scholar
- Bedhomme S, Bernasconi G, Koene JM, Lankinen A, Arathi HS, Michiels NK, Anthes N: How does breeding system variation modulate sexual antagonism?. Biol Lett. 2009, 5: 717-720. 10.1098/rsbl.2009.0401.View ArticlePubMedPubMed CentralGoogle Scholar
- Pongratz N, Michiels NK: High multiple paternity and low last-male sperm precedence in a hermaphroditic planarian flatworm: consequences for reciprocity patterns. Mol Ecol. 2003, 12: 1425-1433. 10.1046/j.1365-294X.2003.01844.x.View ArticlePubMedGoogle Scholar
- Anthes N, Michiels NK: Do "sperm trading" simultaneous hermaphrodites always trade sperm?. Behav Ecol. 2005, 16: 188-195. 10.1093/beheco/arh150.View ArticleGoogle Scholar
- Schärer L: Tests of sex allocation theory in simultaneously hermaphroditic animals. Evolution. 2009, 63: 1377-1405.View ArticlePubMedGoogle Scholar
- Janicke T, Schärer L: Sex allocation predicts mating rate in a simultaneous hermaphrodite. Proc R Soc Lond B. 2009, 276: 4247-4253. 10.1098/rspb.2009.1336.View ArticleGoogle Scholar
- Hoffer J, Ellers J, Koene J: Costs of receipt and donation of ejaculates in a simultaneous hermaphrodite. BMC Evolutionary Biology. 2010, 10: 393-10.1186/1471-2148-10-393.View ArticlePubMedPubMed CentralGoogle Scholar
- Jordaens K, Dillen L, Backeljau T: Effects of mating, breeding system and parasites on reproduction in hermaphrodites: pulmonate gastropods (Mollusca). Anim Biol. 2007, 57: 137-195. 10.1163/157075607780377965.View ArticleGoogle Scholar
- Anthes N, Michiels NK: Reproductive morphology, mating behavior, and spawning ecology of cephalaspid sea slugs (Aglajidae and Gastropteridae). Invert Biol. 2007, 126: 335-365. 10.1111/j.1744-7410.2007.00096.x.View ArticleGoogle Scholar
- Sprenger D, Lange R, Michiels NK, Anthes N: The role of body size in early mating behavior in a simultaneous hermaphrodite, Chelidonura sandrana. Behav Ecol Sociobiol. 2009, 63: 953-958. 10.1007/s00265-009-0738-1.View ArticleGoogle Scholar
- Engqvist L: The mistreatment of covariate interaction terms in linear model analyses of behavioural and evolutionary ecology studies. Anim Behav. 2005, 70: 967-971. 10.1016/j.anbehav.2005.01.016.View ArticleGoogle Scholar
- Quinn GP, Keough MJ: Experimental design and data analysis for biologists. 2002, Cambridge: Cambridge University PressView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.