Low reproductive skew despite high male-biased operational sex ratio in a glass frog with paternal care
© Alexandra et al. 2015
Received: 15 June 2015
Accepted: 27 August 2015
Published: 3 September 2015
Reproductive skew, the uneven distribution of reproductive success among individuals, is a common feature of many animal populations. Several scenarios have been proposed to favour either high or low levels of reproductive skew. Particularly a male-biased operational sex ratio and the asynchronous arrival of females is expected to cause high variation in reproductive success among males. Recently it has been suggested that the type of benefits provided by males (fixed vs. dilutable) could also strongly impact individual mating patterns, and thereby affecting reproductive skew. We tested this hypothesis in Hyalinobatrachium valerioi, a Neotropical glass frog with prolonged breeding and paternal care.
We monitored and genetically sampled a natural population in southwestern Costa Rica during the breeding season in 2012 and performed parentage analysis of adult frogs and tadpoles to investigate individual mating frequencies, possible mating preferences, and estimate reproductive skew in males and females. We identified a polygamous mating system, where high proportions of males (69 %) and females (94 %) reproduced successfully. The variance in male mating success could largely be attributed to differences in time spent calling at the reproductive site, but not to body size or relatedness. Female H. valerioi were not choosy and mated indiscriminately with available males.
Our findings support the hypothesis that dilutable male benefits - such as parental care - can favour female polyandry and maintain low levels of reproductive skew among males within a population, even in the presence of direct male-male competition and a highly male-biased operational sex ratio. We hypothesize that low male reproductive skew might be a general characteristic in prolonged breeders with paternal care.
A wide range of mating systems has evolved in sexually reproducing animals [1, 2]. Mating systems are shaped by the two components of sexual selection: direct competition among individuals for access to mates (intrasexual selection) and the actual mating decision (intersexual selection). These two processes act in the context of general traits of a species and the current environmental conditions, and are shaped by specific characteristics of males and females in a certain population [1, 3]. The prevailing mating system in a group of breeding individuals thus dynamically derives from the optimal, or at least available, mating strategies of all the involved individuals at a given time [4, 5].
Reproductive skew, the uneven distribution of reproductive success among individuals, is a common feature of many animal populations . High male reproductive skew occurs if most females show a preference for a single male (or a few males) in the population, or, in turn, if single males can monopolize the majority of receptive females. Several scenarios have been proposed that would predict high or low levels of reproductive skew, respectively. For example, a highly male-biased operational sex ratio (OSR), which will increase the intensity of male competition for females, is expected to increase the variance in reproductive success among males [1, 7]. The OSR of a given population is, in turn, influenced by parameters such as the population-wide sex ratio, potential reproductive rates in males and females, their mating frequencies, and the asynchronous arrival of females in the breeding population [1, 3, 8].
Amphibian mating patterns are typically divided into two broad categories: explosive breeders, with individuals breeding over a few days to a few weeks, and prolonged breeders, with individuals breeding longer than a month. However, these categories should rather be seen as the two ends of a continuum . High mating pressure on females, the short time-frame available for reproduction, scramble mating with associated sperm competition when females get amplexed by multiple males simultaneously, and stray sperm in aquatic environments presumably preclude the possibility of active female choice in explosively breeding species [10; but see 11]. In species with prolonged breeding, females might be able to selectively choose among available males due to the relaxed temporal constraints. Indeed, discriminatory female choice for certain male traits that indicate physical or genetic quality has been found in field studies of various prolonged breeders (e.g. Engystomops (Physalaemus) pustulosus, ; Uperoleia rugosa, ; Scinax ruber (Ololygon rubra), ; Dendropsophus, (Hyla elegans), ). Several studies have also shown female preferences for certain call characteristics in two-speaker phonotaxis experiments in the laboratory, although under more realistic conditions these preferences sometimes diminished or vanished altogether (reviewed in [16, 17]).
Extended breeding periods are expected to lead to high reproductive skew in males . The rationale behind this hypothesis is that the longer time-frame for mating, and the associated asynchronous arrival of females at the reproductive site, theoretically enables single males to monopolize large shares of the total number of approaching females, compared to explosive breeders. In stark contrast to these theoretical predictions, high male mating skew has been reported for many explosively breeding anuran species (cf. [18, 19]), while comparatively low reproductive skew has recently been reported for a prolonged breeding poison frog .
A complementary conceptual framework on this issue has been proposed in a recent study , suggesting that higher mating skew can be expected in situations where males provide fixed benefits for females (i.e. benefits are not shared or divided between other females under male polygyny, for example ‘good genes’) compared to situations where the provided benefits are dilutable (i.e. benefits are shared between females mating with the same male). One prominent and common example for such a dilutable benefit is male parental performance, as quality and quantity of paternal care can be expected to decrease with an increasing number of eggs/clutches. In amphibians more eggs/clutches are usually associated with higher levels of polygyny . We speculate that low male reproductive skew can be expected in species with paternal care, even in the presence of strong intra-sexual competition and a male-skewed operational sex ratio.
We performed parentage analyses of adults and larvae in a H. valerioi population in order to gain insights into individual mating frequencies, possible mating preferences, and reproductive skew in males and females. We expected the mating system to be at least polygynous (if not polygamous) because males guard several clutches simultaneously. For the same reason, we assumed a higher mating frequency for successful males than for females despite paternal care. And as male reproductive success is probably limited by their access to females, we expected strong intra-sexual competition among H. valerioi males. We chose to combine observational field data and molecular parentage analyses for our study to overcome the logistic constraints on gaining unbiased estimates of mating and reproductive success from field observations alone. Molecular parentage analyses are a powerful tool for reconstructing genealogies within populations and to elucidate individual mating patterns and sexual selection .
During the study period of 101 days (Aug-Nov 2012) we monitored and sampled 142 individual H. valerioi (93 males, 47 females, 2 sex undefined) and 374 of their larvae from 193 clutches (two larvae extracted from 181 clutches, one larva from 12 clutches). For males, we obtained 2419 focal observations in total, 1350 at night and 1069 during the day; for females we obtained 95 observations in total, 64 and 31 at night and day, respectively. Individual males were observed on average 21 times and up to 99 times (median = 21, iqr = 4–41) in the course of the study period, while most females were encountered only once (median = 1, iqr = 1–2, max = 6). The genotype of one male (m138) was not included as a candidate father for the pedigree reconstruction due to low PCR amplification success. However, since no parental genotypes were simulated by COLONY, we assumed zero mating and reproductive success for this male. The MMMeans estimator predicted a total number of 93 males at the study site, which equals the total number of males that we sampled. Since females were only recorded during mating activity, we could not estimate female sampling coverage due to the low number of recaptures.
For 124 clutches (240 larvae) both assigned parents were known from the study area, and for 69 clutches (134 larvae) a known male was assigned together with a simulated female. In total, 30 female genotypes were simulated by the parentage analysis. The genotypes of the unsexed individuals were included in both the putative mothers as well as in the putative fathers for the COLONY analysis. However, both individuals were never assigned as matching parental genotypes to any of the larvae. Thus these individuals were not included in any of the further analyses. For 190 (98.45 %) out of 193 clutches the father that was observed closest to the clutch was also assigned as the genetic father by COLONY. One single clutch was attended by a different male than the one assigned in the parentage analysis. No attending male was observed at two clutches, but paternity of both clutches could be assigned to two sampled males via the parentage analysis. Most clutches (N = 150, 77.72 %) were found after the female had already left the breeding site, thus no predictions about maternity could be made in these cases. For the remaining 43 clutches (22.28 %) where a female was observed in amplexus with a male or close to a recent clutch, the genetic data confirmed the field observations.
The genetic mating system and opportunity for sexual selection
There were significant differences in the number of clutches between males and females, while differences in the number of mates were only marginally significant (Mann–Whitney U-test, Nm/f = 93/77, clutches: U = 2919, P = 0.035; mates: U = 3000, P = 0.064). However, these differences were only due to a larger number of non-reproducers in males, as no significant differences were found for successfully reproducing males and females (Mann–Whitney U-test, Nm/f = 64/72, clutches: U = 2081.5, P = 0.321, mates: U = 2000.5, P = 0.173, Fig. 3).
Opportunity for sexual selection. Standardized variances in mating and reproductive success for male and female H. valerioi, calculated for all individuals and successful reproducers
Number of clutches
I = 0.86
I = 0.32
ΔI clutches = 0.54
I = 0.28
I = 0.23
ΔI clutches = 0.04
Number of mates
I = 0.78
I = 0.33
ΔI mates = 0.45
I = 0.22
I = 0.24
ΔI mates = −0.02
Mating and reproductive success
During our daily night-time surveys we recorded 0–34 (mean ± SD = 17.76 ± 7.79) males present along the transect per night. During individual nights on average only 4.84 % of present males were found in amplexus, while the remaining 95.16 % were encountered alone. Of those males present at the reproductive site 67.82 % were also calling. The parameters ‘observation period’ (i.e. timespan between first and last day encountered), ‘#nights recorded’, and ‘#nights calling’ were all significantly correlated with each other (all rho > 0.819, all P < 0.001), but neither of the parameters correlated significantly with body size (all rho < 0.164, all P > 0.128). All three temporal parameters showed a high variation among males (‘observation period’: mean ± SD = 51.04 ± 35.97; ‘#nights recorded’: mean ± SD = 14.52 ± 12.28; ‘#nights calling: mean ± SD = 10.53 ± 9.26), while the variation in body size among males was rather low (SUL: mean ± SD = 20.74 ± 1.00). We decided to only include ‘#nights calling’ as a surrogate for all temporal parameters, because we considered this variable to have the highest biological relevance with regard to the potential impact on mating success in males. We thus performed a stepwise linear regression to test whether the parameters ‘SUL’ and/or ‘#nights calling’ were significant predictors for the number of clutches or mates per male. The parameter ‘#nights calling’ but not ‘SUL’ proved to be a significant predictor for the number of clutches (#nights calling: t = 5.298, P < 0.001; SUL: t = 0.329, P = 0.743; Model fit: R2 = 0.248, F 1,85 = 28.071, P < 0.001) or mating partners per male (#nights calling: t = 5.003, P < 0.001; SUL: t = 0.652, P = 0.516; Model fit: R2 = 0.228, F 1,85 = 25.033, P < 0.001). The same effect was found when the parameter ‘#nights calling’ was replaced by ‘observation period’ or ‘#nights present’.
For females, body size was not a significant predictor for the total number of clutches, the average number of eggs per clutch, or their rate of clutch production (clutches: t = −0.434, F 1,39 = 0.188, P = 0.667; eggs/clutch: t = 0.318, F 1,34 = 0.101, P = 0.753; mating frequency: t = 1.423, F 1,29 = 2.024, P = 0.166).
Patterns of relatedness
Among all 7084 possible male–female dyads, we observed 138 full-sibling pairs, 1216 half-sibling pairs, and 5646 unrelated pairs, corresponding to probabilities of 1.95 % for full-sib matings, 17.17 % for half-sib matings, and 79.70 % for matings between not closely related individuals when assuming fully random mating. In our actual offspring sample we detected 0.52 % of all matings to have occurred between full-siblings, 20.73 % between half-siblings and 78.76 % between not closely related individuals. The matings we observed did not significantly differ from proportions expected under random mating (2 proportions Z-test; full siblings: Z = 1.342, P = 0.153, half siblings: Z = −1.291, P = 0.197, unrelated: Z = 0.322, P = 0.749). The pairwise relatedness values of females and their actual mating partners did not significantly differ from the relatedness of females to the rest of males (paired t-test, N = 72, t = −0.797, P = 0.426).
The genetic mating system and opportunity for sexual selection
For males, the MMMeans estimator equalled the number of actually sampled individuals. This indicates that for our study population the sampling of H. valerioi males was complete during the observation period. Females were almost exclusively recorded when engaged in mating activities, while otherwise their detection rate, and as a result their recapture rate, was very low. As the calculation of demographic estimators depends on sufficiently high recapture rates, it was thus not possible to reliably estimate female sampling coverage. However, given our long study period, we assume that all females that were in principle able to reproduce during that time were also present at the reproductive site at least once during our sampling. Since females approach males and were never observed to be rejected by males (pers. obs. KT and AM), and since male availability did not seem to be a limiting factor (only about 5 % of the males present at night were found in courtship or already amplexus), we presume that females face no general external restrictions on their opportunity to breed. Therefore we assume that the integration of information from the parentage analysis (i.e. simulated individuals) with the field observations resulted in a reliable approximation of the actual number of females in our study population, thereby yielding an unbiased and representative assessment of the reproductive activity in 2012.
Our results show a high degree of sequential polygamy in our population, where most individuals from both sexes mated more than once and with multiple partners during one single breeding season (cf. ; see also  for similar mating patterns in another frog species). For frogs, this type of mating system is generally associated with prolonged breeding seasons , which are particularly common in tropical environments  as a result of the prevailing climatic conditions. The mating network shows that reproductive activity is evenly distributed across the whole population, and that the majority of individuals form one large reproductive cluster. We observed that almost 70 % of all males recorded at the reproductive site sired at least one clutch, which is very high, compared to estimates of male mating success in other glass frog species (H. fleischmanni: 10 %, C. prosoblepon: 9 %, both based only on behavioural observations, ) and other anurans (cf. ). For the remaining 30 % of males that were unsuccessful we presume that failure to occupy suitable calling sites led to limited access to females. As the ability to establish and defend a calling site is presumably highly energetically expensive, male territoriality likely serves as honest signal for male quality to females. We did not detect alternative reproductive tactics, such as sneaking  or clutch piracy , neither during behavioural observations, nor from genetic parentage assignments. Thus, those alternative strategies might not have evolved in H. valerioi at all, or at least do not play a major role in male reproductive behaviour.
In our study population, we only observed five (6 %) females without reproductive success. All of these females were observed during courtship, but without apparent mating success (i.e. no clutch that could be attributed to those females was found). All other females (94 %) could be assigned to at least one clutch. Given the highly biased operational sex ratio in our study population, we propose that, in general, all females that attempt to mate will also be able to do so. Predation soon after oviposition (i.e. before we were able to genetically sample the clutch) might have been a reason for the lack of reproductive success in these females.
While males' mating success is presumably limited by their ability to advertise and defend a calling site as well as by the number of available females at the stream, females are mainly restricted physiologically by their rate of egg production (cf. ). Mating frequencies were significantly higher in males than in females, but this was likely due to the physiological latency (i.e. recovery time) after ovulation in females, leading to relatively few females with short mating frequencies compared to males (Fig. 4). While several studies have found a positive relationship between the number of eggs per clutch and female body size in anurans [36–38], we did not find such an association in our study, nor were the parameters ‘eggs/clutch’ or ‘mating frequency’ related to body size. Thus we conclude that body size is not a reliable predictor for female fecundity in H. valerioi.
Being choosy or not?
Given the prolonged breeding period, high male site fidelity, and the highly male-biased operational sex ratio in the H. valerioi population under study, females in principle have the possibility to choose their ‘favourite’ male. However, we did not find any indication for females selecting males based on body size or relatedness. No mating preferences in H. valerioi females were also reported in previous studies. The clutch removal experiment by  showed that females do not preferentially choose males that attend a high number of clutches, although previous mating success could potentially serve as an indicator of male quality (e.g. female copying ). In turn, a high number of clutches might even limit the calling activity in males. Likewise,  did not find a correlation between male mating success and the surface area of the leaf or the height of oviposition sites.
Female preference for large males has been found in several anuran species . In H. valerioi, bigger males might be more attractive to females as they are potentially better in defending high quality reproductive resources from other males or protecting egg clutches from approaching predators. Our findings however show that females did not preferentially mate with bigger males.
Inbreeding avoidance is known from other amphibian species (e.g. American toads, Bufo americanus, ). However, when looking at the relatedness of mating partners, our results show no preference for males of low relatedness in H. valerioi females, i.e. the relatedness between mating partners was as high/low as would be expected by chance. In our study population, we found a rather low risk of inbreeding (mating probability between full siblings was 2 %), and due to the lack of this selective pressure, assortative mating by relatedness might not have evolved in H. valerioi.
Under female choice, two general strategies might increase female reproductive success. On the one hand, if there is high variation in the quality of single males, females are expected to evolve preferences for certain male traits that serve as ‘honest signals’ for male quality . On the other hand, females might also increase their reproductive success by allocating successive matings to multiple mating partners [42, 43]. Being choosy or not may also depend on the trade-off between finding the best suitable mate, and time and energy spent searching for him . Moreover, benefits provided by males might decrease with increasing level of male polygyny, as is the case for parental care (cf. ; see also ). Females of H. valerioi have up to five mating partners during one reproductive season, one of the highest levels of polyandry reported in anurans so far [cf. 20, 45]. The fact that H. valerioi females rarely mate twice with the same male could be because females actively try to maximise the number of mating partners. Polyandry has repeatedly been discussed as a powerful strategy to increase fitness by spreading the risk of low quality matings, as well as to increase the genetic diversity among offspring [42, 46]. Through polyandrous mating females may also distribute clutches across many males, which might lower the risk of insufficient parental care by a single male (cf. [20, 45]; see  for partitioning of larvae). However, as the number of clutches and mating partners of successful males and females did not differ from a Poisson distribution, we instead interpret the observed polyandry as a result of females roaming through the breeding habitat and mating with the closest, or most prominent, calling male present at the time of ovulation (cf. [20, 48]).
Male reproductive success
We found that only a fraction of males from the entire study population was present at the reproductive site at any individual night. The observation period, the total number of nights present at the stream, and the number of nights spent calling were significantly correlated with each other and also highly variable across males. The parameter ‘nights calling’ was taken as a proxy for all temporal parameters, because we considered this variable to be most biologically relevant concerning its potential impact on mating success in males. Indeed, the number of nights a male spent calling at the stream was a significant predictor for his mating and reproductive success. This indicates that the variance in success across males can largely be attributed to differences in time spent calling at the reproductive site. Similar results were found for the glass frog species H. fleischmanni and Espadarana prosoblepon: mating success was correlated with the number of nights a male was present at the site, but not with body size, call pitch, or call duration [33, 49]. A significant positive relationship between male mating success and the number of nights present at the reproductive site has also been found in other anuran species (e.g. Bufo woodhousei, ; Physalaemus pustolosus, ; Bufo calamita, ; Hyla chrysoscelis, ; Phyllomedusa rohdei, ; all studies based only on behavioural observation). As calling, mating, and caring can be energetically expensive, the difference in the presence of males at the reproductive site probably reflects the variance in the males’ physiological constitutions (cf. [54, 55]). Although prolonged calling activity might in principal serve as an ‘honest signal’ for male quality to females, in practice, males that call more often will at the same time automatically increase their likelihood of being chosen as a mating partner simply because they are more conspicuous than silent males. If this is the case, variation in mating success is not due to selective female choice but rather due to direct male-male competition. Nonetheless, if presence and calling at the reproductive site actually relates to male quality, females that mate with any closely calling male are more likely to mate with high-quality males without incurring costs of extended mate searching (cf. ).
Low reproductive skew - a general pattern in prolonged breeders with paternal care?
In prolonged breeders, receptive females enter the breeding site successively and/or cyclically (as in the case of iteroparity), rather than simultaneously . In such a highly male-skewed situation they might more easily choose their ‘best mate’ from the available males than explosive breeders. Following the same logic, but from the male perspective, this means that single ‘outstanding’ males might be able to monopolize most of the approaching females. Such a scenario is particularly likely if the ‘quality’ of a specific male – if assessable – is equally attractive to all females, and the benefits of a good male will not decrease as the number of matings increases (i.e. fixed benefits, such as ‘good genes’). However, our data did not support the hypothesis that extended breeding periods and a highly male-biased OSR lead to high reproductive skew across males. In turn, we found high levels of reproductive success in males and females. In our system, it is presumably less costly for females to be polyandrous than to achieve monandry. On the one hand, males are not present and do no call every night (, this study), thus a female’s previous mating partner might not always be available when she is ready to mate. When male territoriality, advertised by calling, as such is an honest signal, females who approach any nearby calling male can be sure to mate with a male of sufficient quality. On the other hand, polyandry in our study population presumably helps to further reduce the risk of repeatedly choosing a male who is genetically inferior or incompatible, or who provides insufficient parental care (cf. ). Fitness benefits associated with parental care distributed across several males and locations might have outbalanced any alleged benefits of single males, and thereby precluded the evolution of mate preferences in H. valerioi. These findings are likely to generalize to many other prolonged breeding species, particularly other vertebrates with male parental care (see  for a summary on parental care in vertebrates). When social and environmental conditions associated with reproduction and offspring survival are highly unpredictable, strong fitness benefits can be expected for females that spread matings across time, space, and mating partners. Such partitioning strategies might be particularly beneficial for females in species with paternal care. Typically, male parenting quality and quantity are more unpredictable than female care, as generally males have gametes more readily available than females and can more easily increase their reproductive success with additional matings than is the case for females. We propose that low male reproductive skew and high rates of polyandry might be common features in prolonged breeding species with paternal care, even in the presence of direct male-male competition, and male-biased sex ratios. Supporting evidence for this hypothesis has also been found in a recent study on mating patterns in a natural population of a hyper-dispersed, territorial frog with paternal care . In this study 35.5 % of the males produced offspring that survived until adulthood, and, due to the commonly high mortality at the larval stage , significantly higher male mating success rates can be assumed. In this species, selection was found to act only on the males’ opportunity to breed (i.e. whether they were able to establish a territory), while females mated non-discriminatorily with the advertising males. To better understand how differential selective pressures act on mating systems and patterns of reproductive success in natural populations, comparative studies with other glass frog species that have less intense (e.g. H. fleischmanni, ) or absent (e.g. Cochranella pulverata, ) parental care are needed.
For our study we combined observational field data with molecular parentage analyses to achieve maximal information on individual mating activities. Our results show that based on field observations alone paternity could be assigned for more than 99 % of all clutches, while for one out of 191 cases the genetic assignments did not match the field observations. This single case was due to a territory takeover after the clutch has been produced. Females were present at only 21 % of the clutches, and their maternity was confirmed in all cases by the genetic analyses. This emphasizes the advantage of molecular tools when investigating reproductive activity in species where individuals are highly cryptic, mobile and/or do not remain close to their offspring.
In the present study we monitored and genetically sampled a natural population of H. valerioi over one breeding season and used microsatellite markers and parentage analysis to characterize the genetic mating system. Our results suggest that females are not choosy, and mate indiscriminately with available males. The variance in male mating success could largely be attributed to the variation in the time that individual males spent at the reproductive site, which in turn could reflect differences in their physiological constitutions. We found low reproductive skew in female and male H. valerioi and high levels of polyandry, which were amongst the highest levels of polyandry reported in anurans so far. Our findings support the hypothesis that dilutable male benefits - such as parental care - can favour female polyandry and maintain low levels of reproductive skew within a population. We show that such low levels of skew can occur even in the presence of direct male-male competition and a highly male-biased operational sex ratio. We hypothesize that low male reproductive skew is a general characteristic in prolonged breeders with paternal care.
The H. valerioi population under study is located in tropical lowland rainforest along the ‘Quebrada Negra’, a small lowland stream close to the tropical research station ‘La Gamba’, situated on the Pacific side of Costa Rica (N 8°42’61”, W 83°12’97”). Detailed characteristics of the drainage system were published by . The area receives on average 5923 mm of rainfall annually with a pronounced wet season between May and November. The highest precipitation occurs from August to November, and the mean annual temperature at the site is 28.2° . The study area is situated within a secondary forest located at the edge of the ‘Piedras Blancas National Park’ at 70 m above sea level. In a preliminary survey, H. valerioi was only detected along a 425 m section of the stream, which then was chosen as the transect for monitoring and sampling adult frogs and clutches. The stream has its source in the hilly slopes at the margins of primary forest. After passing several small waterfalls, the stream then (still very narrow at that point) flows through a small patch of primary forest in a flat area upstream of the transect. The streambed then widens to 3–5 m along the transect where H. valerioi is found. The stream enters open agricultural grassland downstream of the transect, habitat that is not suitable for H. valerioi. During our study we verified at irregular intervals that no H. valerioi were present for at least 150 m both up- and downstream of the transect.
Monitoring and DNA sampling
We performed fieldwork from 15 August to 23 November 2012, covering the period with the highest breeding activity of H. valerioi . We divided the transect into intervals of 5 m and used flagging tape as distance markers. We performed two surveys a day throughout the whole study period, one during the day (between 8 a.m. and 1 p.m.) and one at night (between 6 p.m. and 3 a.m.), with the exception of five days when high water levels due to heavy rainfall impeded field work. For each survey we slowly walked back and forth along the transect. At night we approached every calling male frog and scanned the vegetation for clutches and further frogs. We determined the exact location of frogs and clutches to the nearest meter. During the day we revisited already known clutches and calling sites to check for the presence of frogs and to sample clutches. We used a ladder to reach frogs and clutches up to a maximum height of 5 m above the streambed and bank respectively, depending on the characteristics of ground and vegetation. This dense sampling regime allowed us to gather reliable daily presence-absence data for males. The cryptic female behaviour and the accordingly much lower encounter rate only yielded presence data for females.
We determined the sampling coverage for males by calculating an asymptotic population size estimator (MMMeans, [66, 67]), based on individual capture histories using EstimateS 8.2.0. . Females were almost exclusively recorded when engaged in courtship, leading to a lower detection probability and much lower recapture rate than for males, thereby precluding the reliable estimation of female sampling coverage.
DNA extraction and microsatellite genotyping
For DNA extraction we used a Proteinase K digestion and a standard phenol-chloroform protocol  to isolate total genomic DNA. We amplified DNA samples at eight polymorphic microsatellite loci (Hval04, Hval10, Hval16, Hval17, Hval20, Hval21, Hval22 and Hval24) using PCR primers and protocols described by . After amplification, we diluted the products with water, mixed them with the internal size standard ROX350 (Applied Biosystems) and ran them on an ABI 3130xl Genetic Analyser. We identified all loci visually using PeakScanner 1.0 (Applied Biosystems) and determined final allele sizes using the binning software TANDEM v1.08 . All microsatellite genotypes are given in [Additional file 1].
We carried out parentage assignment using COLONY 2.0 , a program that uses a group-likelihood approach for sibship reconstruction. We used the full likelihood model with medium precision and allowed for polygamous mating in both sexes. We treated both tadpoles that had been sampled from the same clutch as ‘known maternal sibs’. Given the reproductive biology of H. valerioi (e.g. amplexus, male territoriality), we excluded the possibility of multiple maternity within single clutches but did not exclude the possibility of multiple paternity (e.g. through clutch piracy or delayed fertilization by sneakers). When the program could not assign a genotype from the sampled parents to a given offspring-pair, it provided a simulated parental genotype. For the following analyses of reproductive contributions we only used ‘Best (ML) Configuration’ assignments with the maximum likelihood obtained at the end of the computation (cf. COLONY user guide) [see Additional file 2].
The genetic mating system and opportunity for sexual selection
We constructed a network graph showing all inferred matings with the program Cytoscape . We tested for differences in the number of mating partners and clutches between successfully reproducing males and females with Mann–Whitney U-tests. We also tested for differences in the mating frequency (i.e. number of days between two successful matings) in males and females for individuals that had at least two clutches.
To identify sex-specific differences in the opportunity for selection, we calculated the standardized variance in mate acquisitions (I mates ) and clutch production (I clutches ) [2, 74]. This analysis was conducted for all males and females in our study population and again separately for just the successful reproducers. We also determined the Bateman gradient (the relative number of mating partners in relation to the relative number of clutches) for successfully reproducing males and females (cf. ).
Male and female mating success
We investigated whether male body size and/or previous male reproductive activity were significant predictors of male reproductive success [see Additional file 3]. As the parameters ‘observation period’ (timespan between first and last day encountered), ‘#nights recorded’ (number of nights recorded), and ‘#nights calling’ (number of nights calling) were all significantly correlated with each other (see results), ‘#nights calling’ was taken as a proxy for all temporal parameters, since we considered this variable to be most biologically relevant with regard to its potential impact on mating success in males. We then performed a stepwise linear regression to test whether, for each male, body size (SUL) or ‘#nights calling’ were significantly correlated with either the number of clutches or the number of mates.
To investigate whether female body size is related to fecundity, we tested whether the parameter SUL was a significant predictor for the average number of eggs per clutch (‘eggs/clutch’), the total number of clutches per female (‘clutches’), or their mating frequency using stepwise linear regression [see Additional file 4].
We used the program KINGROUP v2  to determine the pairwise relatedness coefficients r  for all male–female pairs. For the simulated individuals, we used the genotype with the highest probability as calculated by COLONY. We used the implemented simulation function to obtain reference intervals (first to third quartile) for expected pairwise relatedness values for 100 full siblings, 100 half siblings, and 100 unrelated individuals, which were [0.35;061],[0.14;0.35], and [−0.11;0.11], respectively. Accordingly, we defined full-siblings as all individual pairs with r values > 0.351, half siblings with values between 0.123 and 0.351, and unrelated individuals with values below 0.123. To test whether mating is random with respect to relatedness, we calculated expected probabilities for full and half siblings to mate in our study population under the assumption of random mating, and compared them with the observed numbers of matings between full and half siblings. Furthermore, we tested whether pairwise relatedness values of females and their mating partners (‘r_mates’) differed from their relatedness to all other males (‘r_nonmates’) using a paired t-test.
All statistical tests were performed in IBM SPSS Statistics 22.0.0. Normality of variables was tested with the Kolmogorov–Smirnov test. Medium values, interquartile ranges (iqr), and non-parametric tests were applied in those cases where normality was rejected.
Availability of supporting data
The datasets supporting the results in this article are included in the additional files.
We want to thank Werner Huber for logistic support in Costa Rica, the staff members of the Biological Research Station La Gamba for support, and Conagebio for issuing the research permits. Work in the field was funded by the University of Vienna (Förderstipendium) and the state Vorarlberg (Vorarlberg Stipendium). Work in the laboratory was supported by the Division of Tropical Ecology and Animal Biodiversity at the Department of Botany and Biodiversity Research (University of Vienna) and by the Austrian Science Fund (FWF): P24788-B22 (PI: Eva Ringler). The Biological Research Station La Gamba is supported by the University of Vienna. Thanks to Gesche Westphal-Fitch for language editing.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Shuster SM, Wade MJ. Mating systems and strategies. Princeton: Princeton University Press; 2003.Google Scholar
- Shuster SM. Sexual selection and mating systems. Proc Natl Acad Sci U S A. 2009;106:10009–16.PubMed CentralView ArticlePubMedGoogle Scholar
- Emlen ST, Oring LW. Ecology, sexual selection, and the evolution of mating systems. Science. 1977;197:215–23.View ArticlePubMedGoogle Scholar
- Kokko H, Jennions MD. Parental investment, sexual selection and sex ratios. J Evol Biol. 2008;21:919–48.View ArticlePubMedGoogle Scholar
- Alonzo SH. Social and coevolutionary feedbacks between mating and parental investment. TREE. 2009;25:99–108.PubMedGoogle Scholar
- Nonacs P, Hager R. The past, present and future of reproductive skew theory and experiments. Biol Rev. 2011;86:271–98.View ArticlePubMedGoogle Scholar
- Klug H, Heuschele J, Jennions MD, Kokko H. The mismeasurement of sexual selection. J Evol Biol. 2010;23:447–62.View ArticlePubMedGoogle Scholar
- Kokko H, Klug H, Jennions MD. Unifying cornerstones of sexual selection: operational sex ratio, Bateman gradient and the scope for competitive investment. Ecol Lett. 2012;15:1340–51.View ArticlePubMedGoogle Scholar
- Wells KD. The social behaviour of anuran amphibians. Anim Behav. 1977;25:666–93.View ArticleGoogle Scholar
- Wells KD. The ecology and behavior of amphibians. Chicago: The University of Chicago Press; 2007.View ArticleGoogle Scholar
- Hettyey A, Baksay S, Vági B, Hoi H. Counterstrategies by female frogs to sexual coercion by heterospecifics. Anim Behav. 2009;78:1365–72.View ArticleGoogle Scholar
- Ryan MJ. The Tungara Frog. Chicago: University of Chicago Press; 1985.Google Scholar
- Robertson JGM. Female choice, male strategies and the role of vocalizations in the Australian frog Uperoleia rugosa. Anim Behav. 1986;34:773–84.View ArticleGoogle Scholar
- Bourne GR. Lekking behavior in the neotropical frog Ololygon rubra. Behav Ecol Sociobiol. 1992;31:173–80.View ArticleGoogle Scholar
- Bastos RP, Haddad CFB. Breeding activity of the neotropical treefrog Hyla elegans (Anura, Hylidae). J Herpetol. 1996;30:355–60.View ArticleGoogle Scholar
- Jennions MD, Petrie M. Variation in mate choice and mating preferences: a review of causes and consequences. Biol Rev. 1997;72:283–327.View ArticlePubMedGoogle Scholar
- Gerhardt HC, Huber F. Acoustic communication in insects and anurans: common problems and diverse solutions. Chicago: University of Chicago Press; 2002.Google Scholar
- Schmeller DS, Merilä J. Demographic and genetic estimates of effective population and breeding size in the amphibian Rana temporaria. Conserv Biol. 2007;21:142–51.View ArticlePubMedGoogle Scholar
- Ficetola GF, Padoa‐Schioppa E, Wang J, Garner TW. Polygyny, census and effective population size in the threatened frog. Rana latastei Anim conserv. 2010;13:82–9.View ArticleGoogle Scholar
- Ursprung E, Ringler M, Jehle R, Hödl W. Strong male/male competition allows for nonchoosy females: high levels of polygynandry in a territorial frog with paternal care. Mol Ecol. 2011;20:1759–71.View ArticlePubMedGoogle Scholar
- Tazzyman SJ, Seymour RM, Pomiankowski A. Fixed and dilutable benefits: female choice for good genes or fertility. P Roy Soc Lond B Bio. 2012;279:334–40.
- Summers K. Paternal care and the cost of polygyny in the green dart-poison frog. Behav Ecol Sociobiol. 1990;27:307–13.View ArticleGoogle Scholar
- Savage JM. The amphibians and reptiles of Costa Rica: a herpetofauna between two continents, between two seas. Chicago: University of Chicago Press; 2002.Google Scholar
- Cisneros-Heredia DF, McDiarmid RW. Revision of the characters of Centrolenidae (Amphibia: Anura: Athesphatanura), with comments on its taxonomy and the description of new taxa of glassfrogs. Zootaxa. 2007;1572:1–82.Google Scholar
- Kubicki B. Ranas de vidrio de Costa Rica. Santo Domingo de Heredia: Instituto Nacional de Biodiversidad INBio; 2007.Google Scholar
- McDiarmid RW, Adler K. Notes on territorial and vocal behavior of neotropical frogs of the genus Centrolenella. Herpetologica. 1974;30:75–8.Google Scholar
- Vockenhuber EA, Hödl W, Karpfen U. Reproductive behaviour of the glass frog Hyalinobatrachium valerioi (Anura: Centrolenidae) at the tropical stream Quebrada Negra (La Gamba, Costa Rica). Stapfia. 2008;88:335–48.Google Scholar
- Vockenhuber EA, Hödl W, Amézquita A. Glassy fathers do matter: egg attendance enhances embryonic survivorship in the glass frog Hyalinobatrachium valerioi. J Herpetol. 2009;43:340–4.View ArticleGoogle Scholar
- Jones AG, Arden WR. Methods of parentage analysis in natural populations. Mol Ecol. 2003;12:2511–23.View ArticlePubMedGoogle Scholar
- Byrne PG, Roberts JD. Evolutionary causes and consequences of sequential polyandry in anuran amphibians. Biol Rev. 2012;87:209–28.View ArticlePubMedGoogle Scholar
- Arak A. Male-male competition and mate choice in anuran amphibians. In: Bateson P, editor. Mate choice. Cambridge: Cambridge University Press; 1983. p. 181–210.Google Scholar
- Crump ML. Reproductive strategies in a tropical anuran community. Univ Kans Mus Nat Hist Misc Publ. 1974;61:1–68.Google Scholar
- Jacobson SK. Reproductive behavior and male mating success in two species of glass frogs (Centrolenidae). Herpetologica. 1985;41:396–404.Google Scholar
- Zamudio KR, Chan LM. Alternative reproductive tactics in amphibians. In: Oliveira R, Taborsky M, Brockmann J, editors. Alternative Reproductive Tactics: An Integrative Approach. Cambridge: Cambridge University Press; 2008. p. 300–31.View ArticleGoogle Scholar
- Vieites DR, Nieto-Román S, Barluenga M, Palanca A, Vences M, Meyer A. Post-mating clutch piracy in an amphibian. Nature. 2004;431:305–8.View ArticlePubMedGoogle Scholar
- Howard RD, Kluge AG. Proximate mechanisms of sexual selection in wood frogs. Evolution. 1985;39:260–77.View ArticleGoogle Scholar
- Cummins CP. Temporal and spatial variation in egg size and fecundity in Rana temporaria. J Anim Ecol. 1986;55:303–16.View ArticleGoogle Scholar
- Tejedo M. Effects of body size and timing of reproduction on reproductive success in female natterjack toads (Bufo calamita). J Zool. 1992;228:545–55.View ArticleGoogle Scholar
- Jones AG. On the opportunity for sexual selection, the Bateman gradient and the maximum intensity of sexual selection. Evolution. 2009;63:1673–84.View ArticlePubMedGoogle Scholar
- Unger LM, Sargent RC. Allopaternal care in the fathead minnow, Pimephales promelas: females prefer males with eggs. Behav Ecol Sociobiol. 1988;23(1):27–32. doi:10.1007/BF00303054.View ArticleGoogle Scholar
- Waldman B, Rice JE, Honeycutt RL. Kin recognition and incest avoidance in toads. Am Zool. 1992;32:18–30.Google Scholar
- Yasui Y. The genetic benefits' of female multiple mating reconsidered. TREE. 1998;13:246–50.PubMedGoogle Scholar
- Jennions MD, Petrie M. Why do females mate multiply A review of the genetic benefits. Biol Rev. 2000;75:21–64.View ArticlePubMedGoogle Scholar
- Johnstone RA, Reynolds JD, Deutsch JC. Mutual mate choice and sex differences in choosiness. Evolution. 1996;50:1382–91.View ArticleGoogle Scholar
- Byrne PG, Keogh JS. Extreme sequential polyandry insures against nest failure in a frog. P Roy Soc Lond B Bio. 2009;276:115–20.View ArticleGoogle Scholar
- Evans JP, Magurran AE. Multiple benefits of multiple mating in guppies. Proc Natl Acad Sci U S A. 2000;97:10074–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Erich M, Ringler M, Hödl W. Ringler E Brood-partitioning behaviour in unpredictable environments: hedging the bets? Behav Ecol Sociobiol. 2015;69:1011–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Meuche I, Brusa O, Linsenmair KE, Keller A, Pröhl H. Only distance matters – non-choosy females in a poison frog population. Front Zool. 2013;10:29.PubMed CentralView ArticlePubMedGoogle Scholar
- Greer BJ, Wells KD. Territorial and reproductive behavior of the tropical American frog Centrolenella fleischmanni. Herpetologica. 1980;36:318–26.Google Scholar
- Woodward B. Male persistence and mating success in Woodhouse's toad (Bufo woodhousei). Ecology. 1982;63:583–5.View ArticleGoogle Scholar
- Arak A. Female mate selection in the natterjack toad: active choice or passive attraction? Behav Ecol Sociobiol. 1988;22:317–27.Google Scholar
- Morris MR. Female choice of large males in the treefrog Hyla chrysoscelis: the importance of identifying the scale of choice. Behav Ecol Sociobiol. 1989;25:275–81.View ArticleGoogle Scholar
- Wogel H, Abrunhosa PA, Pombal Jr JP. Breeding behaviour and mating success of Phyllomedusa rohdei (Anura, Hylidae) in south-eastern Brazil. J Nat Hist. 2005;39:2035–45.View ArticleGoogle Scholar
- Murphy CG. Determinants of chorus tenure in barking treefrogs (Hyla gratiosa). Behav Ecol Sociobiol. 1994;34:285–94.View ArticleGoogle Scholar
- Halliday T, Tejedo M. Intrasexual selection and alternative mating behaviour. In: Heatwole H, Sullivan BK, editors. Amphibian Biology, vol. 2 Social behaviour. Chipping Norton: Surrey Beatty; 1995. p. 419–68.Google Scholar
- Friedl TW, Klump GM. Sexual selection in the lek-breeding European treefrog: body size, chorus attendance, random mating and good genes. Anim Behav. 2005;70:1141–54.View ArticleGoogle Scholar
- Fox CW, Rauter CM. Bet-hedging and the evolution of multiple mating. Evol Ecol Res. 2003;5:273–86.Google Scholar
- Balshine S. Patterns of parental care in vertebrates. In: Royle NJ, Smiseth PT, Kölliker M, editors. The evolution of parental care. Oxford: Oxford University Press; 2012. p. 62–80.View ArticleGoogle Scholar
- Vonesh JR, De la Cruz O. Complex life cycles and density dependence: assessing the contribution of egg mortality to amphibian declines. Oecologia. 2002;133:325–33.View ArticleGoogle Scholar
- Delia J, Cisneros-Heredia DF, Whitney J, Murrieta-Galindo R. Observations on the reproductive behavior of a Neotropical glassfrog, Hyalinobatrachium fleischmanni (Anura: Centrolenidae). S Am J Herpetol. 2010;5:1–12.View ArticleGoogle Scholar
- Schiemer F, Huber W, Weissenhofer A. Stream ecosystems of Costa Rica. Vienna: Verein zur Förderung der Tropenstation La Gamba; 2010.Google Scholar
- Rakosy D, Speckmaier M, Weber A, Huber W, Weissenhofer A. Orchids: Botanical jewels of the Golfo Dulce Region, Costa Rica. Vienna: Verein zur Förderung der Tropenstation La Gamba; 2013.Google Scholar
- Abramoff MD, Magalhaes PJ, Ram SJ. Image processing with imageJ. Biophoton Int. 2004;11:36–42.Google Scholar
- Gosner KL. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica. 1960;16:183–90.Google Scholar
- ASAB. Guidelines for the treatment of animals in behavioural research and teaching. Anim Behav. 2006;71:245–53.View ArticleGoogle Scholar
- Colwell RK, Mao CX, Chang J. Interpolating, extrapolating, and comparing incidence-based species accumulation curves. Ecology. 2004;85:2717–27.View ArticleGoogle Scholar
- Brose U, Martinez ND. Estimating the richness of species with variable mobility. Oikos. 2004;105:292–300.View ArticleGoogle Scholar
- Colwell RK, Estimate S. Statistical Estimation of Species Richness and Shared Species from Samples. Release 8.2.0. 2006. Available from: http://purl.oclc.org/estimates.Google Scholar
- Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual. 2nd ed. New York: Cold spring harbor laboratory press; 1989.Google Scholar
- Ringler E, Mangold A, Trenkwalder K, Ringler M, Hödl W. Characterisation of nine new polymorphic microsatellite loci in the reticulated glass frog Hyalinobatrachium valerioi (Centrolenidae). Amphibia-Reptilia. 2014;35:243–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Matschiner M, Salzburger W. TANDEM: integrating automated allele binning into genetics and genomics workflows. Bioinformatics. 2009;25:1982–3.View ArticlePubMedGoogle Scholar
- Jones OR, Wang J. COLONY: a program for parentage and sibship inference from multilocus genotype data. Mol Ecol Res. 2010;10:551–5.View ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. CYTOSCAPE: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.PubMed CentralView ArticlePubMedGoogle Scholar
- Wade MJ, Arnold SJ. The intensity of sexual selection in relation to male sexual behaviour, female choice, and sperm precedence. Anim Behav. 1980;28:446–61.View ArticleGoogle Scholar
- Konovalov DA, Manning C, Henshaw MT. Kingroup: a program for pedigree relationship reconstruction and kin group assignments using genetic markers. Mol Ecol Notes. 2004;4:779–82.View ArticleGoogle Scholar
- Queller DC, Goodnight KF. Estimating relatedness using genetic markers. Evolution. 1989;43:258–75.View ArticleGoogle Scholar