- Open Access
Only distance matters – non-choosy females in a poison frog population
Frontiers in Zoologyvolume 10, Article number: 29 (2013)
Females have often been shown to exhibit preferences for certain male traits. However, little is known about behavioural rules females use when searching for mates in their natural habitat. We investigated mate sampling tactics and related costs in the territorial strawberry poison frog (Oophaga pumilio) possessing a lek-like mating system, where both sequential and simultaneous sampling might occur. We continuously monitored the sampling pattern and behaviour of females during the complete period between two successive matings.
We found no evidence that females compared males by visiting them. Instead females mated with the closest calling male irrespective of his acoustic and physical traits, and territory size. Playback experiments in the natural home ranges of receptive females revealed that tested females preferred the nearest speaker and did not discriminate between low and high call rates or dominant frequencies.
Our results suggest that females of O. pumilio prefer the closest calling male in the studied population. We hypothesize that the sampling tactic in this population is affected by 1) a strongly female biased sex ratio and 2) a low variance in traits of available males due to strong male-male competition, preventing low quality males from defending a territory and mating.
The function of extravagant male characters and females preferences for these traits has been a major focus of research in behavioural and evolutionary biology . In recent years, particular attention has been concentrated on determining the origin and adaptiveness of female mating preferences. It has been demonstrated that females can gain direct  and indirect benefits [3–6] from choosing among males and from their preferences for certain male traits.
Despite a steady increase in empirical and theoretical investigations about sexually selected traits , relatively few studies have been devoted to information gathering by females about male traits and decision making based on that information [8, 9]. Detailed knowledge about female mate sampling behaviour is important because it directly affects the evolutionary trajectory of male characteristics [10–12]. Such knowledge is also important to understand how selection acts on female mating preferences, since the expression of mating preferences depends on the mate sampling tactic. Several theoretical models as to how females sample potential mates have been proposed [10, 11, 13–18]. The simplest and least costly female mate sampling tactic (in terms of time and energy expenditure) is to mate with the first male encountered . This random mating implies that accepted males can be of any quality . Alternatively, females may use a search tactic, in which they compare a number of males and then choose the male with the highest quality within their sample . Under this scenario, females may compare only the last encountered males [21, 22], a certain number of males ["best-of-N" tactic: 10] [11, 17] or sample males for a certain amount of time  before returning to the highest quality male. The quality of accepted males should be higher for the respective female than the quality of rejected males . Females may also use a threshold tactic whereby the first male that exceeds an established threshold criterion is accepted [10, 11, 22]. Females adopting a fixed threshold criterion use the same threshold regardless of the quality of males that are available . Females adopting an adjustable threshold criterion adjust their threshold to the quality of available males [10, 11, 17]. Both tactics predict that females always mate with the last encountered male, whereas the number of sampled males can vary. This implies that some females accept the first male encountered and that the selected males should be of higher quality than males rejected on first visits .
Mate sampling has been investigated in several taxonomic groups [9, 24–26]. Consistent with predictions of theoretical models, empirical studies have revealed that the female mate sampling tactic is influenced by the costs of searching [21, 27–31] as well as the variation in quality among potential males [32–35]. Females should become choosier with an increase in benefits and a decrease in costs [36, 37]. In contrast to Janetos’s  conclusions, the best-of-N tactic does not inevitably lead to the highest fitness for females. When search costs are considered, an adjustable threshold tactic (one-step decision tactic) generates higher fitness gain . If search costs are very high, females should accept the first male they encounter.
Theoretical models of mate choice [10, 18, 22] assume that females assess males sequentially. However, especially in lekking and chorusing species of different taxonomic groups, females have opportunities to assess several males simultaneously without closely approaching them [24, 38–40] and can thereby avoid costs associated with a close approach . In these species, sampling is more likely to be limited by the perceptual abilities of females in distinguishing males, and females must weigh the potential for improvements in detection of male quality against the costs of reaching males .
Observations of anuran mate sampling suggest that both sequential  and simultaneous [24, 43, 44], assessment can occur, making them an outstanding system for studies on mate sampling tactics. Several studies have revealed that preferences of female anurans are based on male advertisement calls . Males that produce energetically very costly calls are often more attractive for females . Different call properties require differing energy investments and females were shown to prefer calls with more components, higher intensity, longer duration or higher rate of repetition . Some properties of acoustic signals (e.g. dominant frequency of calls) can contain information about male body size  which often influences female mate choice decisions [40, 47].
In the present study we investigated female preferences for certain male characteristics (acoustic and physical traits, territory size) as well as mate sampling tactics in the strawberry dart-poison frog (Oophaga pumilio). Within the last years, the ecology, reproductive and territorial behaviour, sexual selection as well as geographic, phenotypic and genetic variation in O. pumilio has received a great deal of attention [48–64]. O. pumilio provides an excellent model system for examining mate choice tactics for several reasons. These diurnal aposematically coloured frogs can easily be monitored over several months of the prolonged breeding season when females exhibit, contrary to most other anuran, very high mating rates (average interval between two matings of a female approx. 4–5 days in our monomorphic, red coloured study population ). Due to the lack of an amplexus males have no possibility of forcing a female to mate, thus females are free to choose among potential mates. Mating takes place in the territories of the males and after oviposition males moisten the terrestrial eggs in the leaf litter. Ten days after oviposition females start transporting their tadpoles to small water-filled leaf axils of bromeliads or bananas (only one tadpole per axil) and feed them with nutritive eggs [65–68]. O. pumilio has a lek-like mating system with long-term territoriality. Males use advertisement calls to attract females, and advertisement calls, aggressive calls and physical combat to defend their territories against intruding males. Males compete for access to mating partners and their distribution depends on the distribution of the females [49, 55] while females defend feeding territories  and are limited in their distribution to areas with sites suited for tadpole-rearing . In our study population females possess large overlapping home ranges (~36 m2; covering several male territories) but smaller territories (~2 m2) than males (~15 m2) [49, 55, 69].
In lek-like mating systems females have opportunities to compare several males sequentially and simultaneously. If female O. pumilio sample males sequentially we expect to find movements among potential mates and acceptance of mates according to the predictions of theoretical models [10, 11, 16, 22]. However, it is possible that females sample several males simultaneously from a greater distance. Therefore we also investigated whether females preferred and chose individuals with certain characteristics (see above) regardless of whether females had visited these males before or not.
Because most data on mate choice in amphibians come from experiments, Richards-Zawacki et al.  emphasised the need and importance of field studies for understanding the patterns of mate choice exhibited in natural populations. However, due to the complexity of female mate choice, observations of mate sampling behaviour alone may also be insufficient and experiments are required to distinguish among different mate choice tactics [17, 43, 70, 71]. Thus, in addition to field observations on female behaviour we used playback experiments to verify female preferences for certain male call traits. Since in a previous study male mating success of strawberry poison frogs was positively correlated with dominant frequency and call rate , we tested female preferences for these two call properties.
In many taxonomic groups observations on female mate sampling are restricted to a certain period before mating [e.g. [24, 28, 39, 72, 73]]. To our knowledge no information is available on the behaviour of females for the complete period between two matings. In the present study, we continuously observed individual females of the strawberry poison frog during the whole period between two successive ovipositions and simultaneously monitored the behaviour and spatial distribution of the surrounding males.
We observed the mating behaviour of 20 females (N = 5 in 2004; N = 15 in 2005). The mean intervals between two successive ovipositions were 5.8 ± 2.44 days. Ovipositions were initiated between 06:12 and 11:30 am. During courtship males led females to oviposition places (typically dry leaves) inside their territory. Males left the oviposition site 18 ± 4 min (median ± interquartile range; N = 19) min after entering the oviposition leaf. Thus they stayed a significantly shorter time at the oviposition site (Wilcoxon matched pairs test: T = 0, Z = 3.82, N = 19, p = 0.0001) than females (median ± interquartile range: 45 ± 13 min; N = 19). Clutch size varied between 3 and 8 eggs (mean: 4.89 ± 1.09 eggs (N = 20)). In only 4 different occasions we observed males leading more than one courting female successively to exactly the same oviposition place. These females deposited their clutches closely next to each other and it was not possible to distinguish between both clutches.
Hatching success was determined 7 days after oviposition because females (N = 4) revisited their clutches 7.5 ± 0.6 days after oviposition. In both years hatching success was very low. Only 20.38 ± 25.17% of the eggs per clutch (N = 19 females) survived until the seventh day after oviposition. Main reasons for egg mortality were mould, desiccation as well as predation (e.g. by ants). Two clutches did not develop.
Males stopped calling during and after oviposition of their mate. The duration of these calling interruptions was on average 32.57 ± 13.1 min (N = 14). Three males did not begin to call again after oviposition in the course of the data collection on the particular day. These males were omitted from the calculations.
Female mate choice behaviour
Number of contact males
We observed the mate choice behaviour of 5 (2004) and 15 (2005) females. In the year 2005 for each female 6.4 ± 1.92 and 11.0 ± 2.62 males were available as mates in a radius of 5 and 10 m, respectively, measured from the border of their territory. However, females showed restricted mate search tactics (Table 1). Twenty-five percent of the focal females (N = 20) were in contact (approximation closer than 50 cm to a male) with only one male during the observation period (time between two successive ovipositions) and chose this male as mate. Additionally, 10% (N = 20) of females were in contact with males only during the time of mating. Twenty-five percent of females (N = 20) accepted a male with which they had no previous contact during the observation period (time between two successive ovipositions).
On the oviposition day females were in contact with significantly more males than on other days (Friedman ANOVA: F = 7.89; N = 20, df = 2; P = 0.02). Nonetheless, on the day of oviposition 85% of females (N = 20) were in contact with only one male and chose this male as a mate. Thus, they accepted the first male they encountered (Figure 1). The other 15% of females (N = 20) were in contact with 1.33 ± 0.58 other males prior to mating on the day of oviposition and accepted the second or third contact male.
All initiated courtships resulted in oviposition except two cases when females laid unfertilised eggs without a male. In both cases courtship was interrupted by external disturbances (e.g. territorial fighting between males, biting ants, or several females trying to mate with the same male simultaneously). The females left the first male and moved fast and directly toward the next male. Also courtship with these males was disturbed by one or more of the mentioned factors. Abruptly both females showed no further courtship behaviour although the males were still calling and oviposition took place in the absence of a male in the late afternoon.
Comparing first and second ovipositions 10 of 20 focal females chose the same male again: In 93% of the cases (N = 15 females, year 2005) the male selected during the first oviposition was calling 0–60 min before initiation of the second mating.
Based on the systematic observation of males, information about 2.47 ± 1.81 (range: 0–6) additional matings of focal females (N = 14) after the regular observation period (from one oviposition to the next oviposition) was available in 2005. Thus, females were observed to choose 2 ± 0.65 (range: 1–3) different males as mates. Only 3 females oviposited repeatedly with the same male. Territories of all mates overlapped with the home ranges of females and were located in a radius of 5 m measured from the border of females’ territory.
Eighty-two percent of all females (N = 11) spent the night before oviposition inside their territories. On the day of oviposition four of 15 females did not leave their territory for mating initiation and approached and accepted the closest calling male as mate (Figure 1). Eleven females initiated mating outside their territory. On the day of oviposition nine of these 11 females left their territory within the last 45 min (22.78 ± 9.91 min) before mating initiation and moved directly towards the (at that moment) closest calling male (N = 7 females) or towards one of the two or three closest calling males that were located at the same distance from the female (N = 2 females) and accepted this male as mate.
Quality of males
The average snout-vent length (SVL) of males (N = 16) in the study area was 23.82 ± 0.80 mm. Their mean weight was 1.04 ± 0.08 g. The average male condition (condition index [74–76]) amounted to −0.02 ± 0.09 g. Descriptive statistics and coefficient of variation (CV) for all call properties are given in Table 2. The ratio CVb/CVw (CV at the between-male level/CV at the within-male level) was always smaller than 2. Thus, differences in CVs between and within males were small. CVs for the number of pulses were higher and CVs for dominant frequency were lower than CVs for other properties.
Although females were in contact with only a few males they may be able to assess calls of males simultaneously without closely approaching them. Because focal females did not mate with a male more distant than the second closest to them (distance calculation on the day of oviposition at the moment when females left their territory; see above), we compared the qualities of the two closest males irrespective of prior contact between a female and a male. Using univariate analysis we did not find female preferences for certain call properties or physical properties (SVL, weight, body condition; Table 3).
Additionally, chosen males did not differ from the average of all other contact males (Additional file 1) or the average of all other males whose territory overlapped with the home range of the females (Additional file 2) regarding the measured traits.
Our PCA analyses yielded 10 dimensions in total, however 74% of male trait variability was representable by the first three dimensions (Table 4, Additional files 3, 4, 5). These components did not correlate with success rate, thus, mate choice by females is not significantly explainable by these factors. This is in accordance with our results of the univariate analyses (see above). We did not find differences in PCA scores between chosen and rejected males (Table 5).
According to our power analyses for the t-tests (with n = 15, p < 0.05, and power = 0.8, paired design), the sample size is good for our parameters with high effect size (territory size and calling activity) and marginally adequate for medium effect sizes (number of pulses and pulse rate). The effect size is low for the remaining parameters (frequency, call duration, call rate, condition, SVL, weight).
Thirty-six of 45 females tested in playback experiments responded by approaching a speaker. Females showed no preference for the higher or lower call rate or higher or lower dominant frequency (Figure 2). However, they showed a clear preference for the closest speaker (Figure 2).
During each experiment there was no difference between females choosing one speaker and females choosing the other speaker with regard to their body condition, SVL, weight, time and air temperature (t-test or U-test; p > 0.05).
Our results suggest that the distance to potential mates is the crucial criterion for mate choice in female strawberry poison frogs in our study population. The tested acoustic and physical traits, and territory size did not seem to influence female mate choice decision.
In many taxa female preferences are based on male acoustic properties [reviewed by ]. Pröhl  found that in O. pumilio males, mating success is positively correlated with calling activity and call rate and negatively with dominant frequency in some study years. However, female preferences for call traits are not the only possible explanation for this observation. Instead, males with certain characteristics (call properties, fighting abilities) could have the ability to defend preferred areas with high female density [49, 55] and thus gain a mating advantage. In many anuran species [78–80] including O. pumilio males reduce their dominant call frequency during social interactions with other males. Possibly, frequency is involved in the assessment of fighting ability between rival males [80, 81]. Since male strawberry poison frogs that defend territories with high female density also have a higher mating success , the correlation between certain call properties of males and their mating success  could be due to intrasexual interactions.
In most anurans, male vocalizations are the primary signal modality. However, many studies suggest that different modalities are involved in information transfer and sexual selection in anurans [82–87] and other taxonomic groups [reviewed by ]. Recent experiments with polymorphic and monomorphic Panamanian populations of O. pumilio revealed that females prefer males of their own phenotype to males with other colours [59, 60, 89–91]. Beside colouration, other parameters such as brightness or conspicuousness may also be important cues for female mate choice . However, the extreme colour variability in this species seems to be restricted to the islands of the Bocas del Toro archipelago and the adjacent mainland of Panama [92, 93]. Whether brightness, conspicuousness and colouration influences female mate choice decisions in the monomorphic red or reddish orange populations of the remaining distribution range in Costa Rica and Nicaragua remains to be addressed in future studies.
Meuche et al.  showed that males possess smaller territories in areas with high female density and strong intrasexual competition for these places. In the present study accepted males tend to occupy smaller territories than their neighbours. This is consistent with a situation where males compete for places with high female density because the best competitors would be closest to more females than males that defend territories with low female density.
Females might prefer other parameters such as resources defended by males [9, 77, 94]. The quality of oviposition sites inside the territory of each male could influence female preferences if the quality of the defended area has an impact on hatching success [95, 96]. Additionally, male parental care can affect hatching success of clutches and female preferences . With one exception (Hyla rosenbergi) there are no indications that female anurans prefer better fathers as mates. In O. pumilio overall hatching success is very low [this study, ] and probably depending on many ecological factors that cannot be controlled by males (e.g. predation, drought). Oviposition sites (mainly rolled-up leaves of bananas, Heliconia and caoutchouc trees) are very abundant and do not seem to be a limited resource . Furthermore, during our study no female rejected a male after visiting the oviposition site. Therefore, we assume that female mate decisions are not influenced by male parental care and the quality of oviposition sites. Additionally, there is no evidence that mate choice in O. pumilio is based on tadpole rearing sites since 1) males do not control the access to larval habitats and 2) females usually transport tadpoles to water bodies outside the territory of the fathers .
We are aware that sample sizes are low for a part of our statistical analyses and their non-significance should be considered with caution. However, as some of the parameters with low effect-size are connected to the high or medium effect-size parameters (I. Meuche, unpublished data), we expect them to be similar regarding their effect on mating success. Strong positive or negative effects should be identifiable at least in some of the parameters, if present. We thus consider low effect-size parameters also as not significant and assume this effect to be consistent even when studying larger sample sizes.
Mate sampling tactic
Our data suggest that in our study population, female strawberry poison frogs use a mate sampling tactic that could be defined as “accept the closest calling male”. The probability that a female mated with a male was not influenced by the analysed characteristics (Table 3) and 85% of all focal females (N = 20) accepted the first male they encountered on the day of oviposition. Non-choosy females or the acceptance of the closest male as mate sampling tactic has also been observed in fishes [101, 102] and some anuran species with a lek-like mating system [24, 103, 104]. For example in the dendrobatid frog Allobates femoralis spatial and genetic pattern of reproductive success suggest that females are likely to mate with any spatially proximate territorial male and only territory possession determines if males are considered as mating partner [105, 106]. Additionally, recent studies suggest that females of Hyla gratiosa use a simultaneous sampling tactic [24, 43]. Observations showed that females moved to a location where they were able to compare several males simultaneously, selecting and approaching the preferred male directly. Such a tactic could also be applied by females of O. pumilio. However in that case we would expect females to choose their mates according to certain characteristics. Our data did not reveal any female preference for male traits that might be possible indicators of fitness (Tables 3, 4, 5).
The choice of the closest caller as a mate was not the result of preliminary movements towards the preferred male. Most females spent the night and morning before oviposition inside their territories before approaching the closest caller.
Experience from previous sampling bouts early in the mating season has been shown to reduce mate sampling in later time periods especially in species were male display sites are stable [107, 108]. In O. pumilio both males  and females (I. Meuche pers. obs.) are defending their territories for several months or years. Although, long-term experience could explain restricted mate sampling found in females of O. pumilio, it is not consistent with the clear preference for the closest caller. Additionally, it is unlikely that the choice of the closest calling male is simply the result of establishing a territory near the preferred males because female distribution does not depend on male distribution  and seems to be associated with feeding sites . However, further research is necessary to determine the influence of long-term experience on mate sampling in O. pumilio.
Selective phonotactic movements toward the closest male could be a cost-reducing strategy in terms of search time and competition for mates. Computer simulations have shown that the number of sampled males decreases as sampling costs increase . Especially time constraints seem to affect the optimal sampling tactic in different taxonomic groups [9, 29, 111–114]. Even in mating systems where females can simultaneously evaluate a large number of mates, the time-related costs would force them to assess only a small number of males  for a relatively short time . In O. pumilio, time constraints seemed to increase the risk of failure to mate. Females that were not able to obtain a mate within a certain period of time laid unfertilised eggs in the absence of a male.
Additionally, indirect competition between females for access to males can cause time constraints if the probability of finding an available male decreases with increasing competition [103, 116]. The intensity of this indirect competition depends on the sex ratio of adults. There is an increased risk of failure to mate when the sex ratio is more biased towards the choosing sex . For example in fishes [117–120] and birds  it has been shown that choosiness is reduced to minimise the risk of not reproducing when the sex ratio is more biased towards the choosing sex. Female strawberry poison frogs are only able to localise males, which call. We found that males did not call for an average of 31 min during and after oviposition of their mate and thus were not localisable from greater distances by other females. Additionally, many males called only if females were in sight. This observation, the strongly female biased sex ratio , and the fact that most matings concentrate on days with optimal climatic conditions , reduced the probability of finding an available male over greater distances in our study population. In contrast to our results, Pröhl & Hödl  reported mate sampling and courtship interruptions by females from a population with an 1:1 sex ratio. We conclude that the mate sampling tactic might be a flexible tactic depending on the actual sex ratio, thereby offering more opportunities for females to be choosy in male biased populations. In contrast, mate sampling and choosiness decrease in more female biased populations . The female biased sex ratio in our study area seems to rely on the much higher abundance of potential tadpole rearing sites (i.e. female reproductive resources) in comparison to areas where the sex ratio is equal [55, 109]. Since female home ranges largely overlap  more females might fit into an area of a certain size. However, the number of males in the same area might be restricted due to their pronounced territoriality and bioacoustic spacing.
Energetic limitations or predation can also influence optimal sampling tactics. Sampling several males can reduce energy reserves that could be allocated otherwise into the offspring . In some studies energetic constraints influenced mate choice [28, 123, 124]. For example, female sticklebacks are highly selective between dull and bright males when sampling costs are low. However, females reduced their selectivity when they had expended energy swimming against a current . No information is available on energetic costs of female mate sampling in the strawberry poison frog; however, it can be assumed that energetic expenditures increase when females are covering greater distances. As females of O. pumilio defend their territory against other females  and the sex ratio is strongly female biased, intensive mate sampling could lead to increasing energetic expenditures [see ], stress and risk of losing the territory to another female.
In different taxonomic groups it has been shown that energetic constraints imposed by parasites are reducing female selectivity [[125–127], but see ]. In O. pumilio ecto- and endoparasites have been found [129–131], and parasite infections seem to vary among individuals  and populations (S. Hagemann, pers. comm.). In our study population we never found ectoparasites in the skin of any of the females or males. Additionally, the load of intestinal parasites load was very low in both sexes [I. Meuche unpublished data, ] in comparison to other anuran species from the same geographic region . Because condition-dependent characteristics of males are not associated with parasite load  we assume that parasites did not impose energetic constraints on the females in our study population.
Predation risk has also been shown to influence female choosiness in other taxa [133, 134]. In contrast to other anurans [104, 135, 136] predation risk for the aposematically coloured populations of strawberry poison frogs seems to be low . However, populations differ in their colouration, toxicity, and local predator communities [63, 93, 138–141]. Recent studies showed that the red colour morph posseses a more effective warning signal and a lower predation rate than other colour morphs [137, 138]. These findings are supported by the fact that in our study area the mortality rate of all territorial males in the course of 12 months was 0%. Therefore, we assume that predation risk and its influence on female mate choice was very low in our red coloured study population. However, further studies are necessary to determine whether the mate sampling tactic used by female strawberry poison frogs is a flexible tactic depending on colouration, availability of defensive alkaloids, and/or predation pressure.
In addition to the costs of intensive mate sampling, benefits also affect the intensity and precision of a preference [33, 142]. When the benefits of intensive mate sampling do not exceed the costs of being choosy, females should be less selective and sample fewer males [11, 36]. Theoretical models show that selectivity increases as variation among potential mates increases because benefits of mating with high quality males increases . If the variance in male quality is low, and the costs of comparing males are high , selection should favour females mating with the first male encountered . In accordance with Pröhl , the variation in bioacoustic traits among males of O. pumilio was low in our study population. Additionally, males defending territories with high female density over several years  are expected to be high quality males. In many populations, very low-quality males are eliminated by male-male competition, so all individuals owning a territory are acceptable [47, 105]. Consequently, females get high-quality mates even if they accept the closest calling male [see ].
High egg mortality as well as the risk of losing the whole clutch by laying unfertilised eggs and the probably low benefits of intensive mate sampling support our assumption that acceptance of the closest calling male represents an optimal mate sampling tactic in female strawberry poison frogs in the investigated population. However, given the large variation in traits (e.g. sex ratio, colouration, toxicity, predation) among different populations of O. pumilio, further studies in a natural situation as well as experimental set-up are needed to evaluate the potential within-species differences in sampling tactics.
The present study was conducted at the Hitoy Cerere Biological Reserve on the Caribbean side of Costa Rica (9°40′ N, 83°05′ W) between August 2004 and May 2006. We established two study areas (located at 100 m above sea level) near the Hitoy Cerere river. The study areas were chosen with regard to high frog density (> 2 males per 100 m2) and accessibility for the researchers. The density was confirmed by capturing and re-capturing frogs at the beginning of the study.
The vegetation consisted mostly of bananas (Musa spp., Musaceae) and Heliconia spp. (Heliconiaceae) as well as Carludovica rotundifolia (Cyclanthaceae) and caoutchouc trees (Castilla elastica). The two areas measured 2400 m2 (2004) and 505 m2 (2005) and were divided by nylon strings into 1 m2 quadrats for the recording of spatial distribution of all frogs.
The study was carried out in accordance with the legal and ethical standards of Costa Rica and was approved by the corresponding government authorities (MINAE/SINAC permit No. 137-2004-OFAU, 038-2005-OFAU, 004 ASP AC/AC). In the following we describe the behavioural observations of males and females. In general, we minimised the impact of our presence on animal behaviour by moving slowly. As a result, no individuals jumped away or seemed to hide when we approached.
From September to December 2004 as well as from April to August 2005, we observed the mate choice behaviour of 20 females for an average of 70 h per female. We successively monitored the behaviour of each focal female from one oviposition to the next oviposition (approximately 5 days) using “focal animal sampling” . The position of the female in the grid system was recorded every 5 min during daylight hours (05:00 am – 06:00 pm). It was not necessary to observe individual frogs during the night because they do not leave their sleeping sites in the leaf litter until dawn (I. Meuche, pers. obs.). Individuals were easily distinguishable by their natural marks such as black spots.
Behavioural observations were made every minute using “one–zero-sampling” . We recorded contacts with males (yes/no), the identity of these males, courtship (yes/no) and mating behaviour (yes/no). In O. pumilio certain male traits (e.g. weight, condition) are not correlated with male advertisement call properties  and, therefore, cannot be assessed by females over large distances. As the dense secondary forest floor creates a visual barrier, it is unlikely that females can assess these properties from a distance of more than 50 cm. Thus, contact with males was defined as an approach closer than 50 cm to a male (or vice versa) . Courtship was an interaction in which a female showed a directional movement toward a calling male and followed this male at least for a short distance . Tactile, visual, and acoustic stimuli were also possible . A mating was defined as a courtship followed by an oviposition event. Because both sexes were visually undetectable after entering the oviposition leaf, the start of an oviposition was defined when both sexes entered the leaf. The oviposition was completed when the leaf was abandoned after egg laying. After each oviposition the number of eggs was determined and the clutches were inspected daily determining the rates and probable causes for mortality.
Data collection was stopped when the focal female transported tadpoles during the observation period (time between two successive ovipositions) as tadpole rearing females do not continue to mate .
To minimise the effect of our presence on female behaviour during data collection, focal females were not captured until the end of each observation period. We measured their weight to the nearest 0.01 g with a Voltcraft PS-250 scale and SVL to the nearest 0.1 mm with a manual calliper. Each focal female was marked by toe-clipping for individual identification – a conventional marking method  that does not affect the behaviour of O. pumilio. We used toe-clipping in this study only after a careful consideration of welfare implications. We only marked two toes or fingers per frog and only one per limb. Furthermore, the first finger was never marked to avoid impacting their climbing abilities. Additionally, each female was photographed for re-identification via its individual patterns of black spots and short lines on the red dorsal ground colour. However, due to changes of individual patterns over time , photos alone were insufficient for unequivocal re-identification in our study population.
At the beginning of the behavioural observations, male frogs (N = 16) were captured for measuring their weight and SVL. Individuals were recognised by distinctive natural marks such as black spots (using photographs) and their toe-clipping pattern. Males were observed over a period of 65 days (from April 2005 to August 2005) by 1–2 other persons during the same time when female behaviour was observed. All observers were trained prior to the observation period ensuring thereby high consistency in data recording.
As calling activity was most intense from 07:00 to 10:00 a.m. and mating as well as oviposition only occurs in the morning , male behaviour was studied from 06:00 to 10:00 a.m. Males were successively located at five minute intervals (“scan sampling” ) and their position in the grid system was recorded. Observation time per male was approximately 15 seconds during each interval. The sampling order of males was maintained during each five minute interval. We recorded calling activity (yes/no) using “one–zero sampling“ . Furthermore, we recorded every mating activity (yes/no) with our focal females to obtain information about additional matings of focal females after the regular observation period (from one oviposition to the next oviposition).
Vocalisations of males were recorded during the morning using a Sony Professional Walkman (WM DC6) and a Sennheiser Directional Microphone (MZA 14 P48) placed approximately 50 cm from the males. We recorded at least six advertisement call sequences for each male on different days. Subsequently, air temperature was measured using a digital thermometer (Voltcraft DT-8820 Environment Meter) with an accuracy of 0.1°C.
We analysed 190 advertisement call sequences of 17 males (mean: 11.2 sequences/male) using the Digital Sound Analysing System Avisoft-SASLab Pro (Fa. Specht, Berlin, Germany) at a sampling rate of 44100 Hz. Power spectrograms (hamming window, FFT-length: 256) were used to measure dominant frequency. Temporal properties were measured using an oscillogram. A call sequence is a series of single calls produced in succession during several seconds or minutes. Between call sequences males pause for irregular time intervals. A call is defined as a complete sequence of pulses separated from the next call by an intercall interval [see ]. Ten consecutive calls of each call sequence were analysed, quantifying the following call properties: call duration (s), dominant frequency (kHz), number of pulses per call and pulse rate (pulses/ms). Additionally, we determined the call rate (calls/s) by counting the number of calls during three periods of two seconds. The mean of the above mentioned call properties was calculated for every call sequence. The duty cycle was determined as an indicator of the percentage of time an acoustic signal was produced within a call sequence. It was calculated as the call rate multiplied by average call duration per call sequence [see [146, 147]]. To avoid variability caused by different thermal conditions, all temperature-dependent call parameters  were adjusted to a temperature of 25.6°C (average temperature during call recordings) using standard regression methods. For every male and each call parameter the average value of the call sequences was used. Additionally, we determined a coefficient of variation (CV) of the call properties at the within-male level (CVw) and between-male level (CVB), and calculated the ratio CVB/CVw.
Territories and home ranges
For each male in the study area as well as each focal female, we calculated the position and size of territories as well as home ranges, respectively. The position data from the grid system were transformed in x- and y-coordinates, in which each coordinate corresponds to the centre of the grid cell. Home range and territory sizes were calculated with the adaptive-kernel method  using the software ArcView GIS 3.3 (Fa. Environmental Systems Research Institute (ESRI)) as well as the supplementary programs Spatial Analyst 2.0a (Fa. ESRI), Cad2Shape 3.0, and Home Range Extension .
Using the kernel method we calculated the density function for 95% of the observation points allowing us to determine the interpolated home ranges (focal females) or territories (males) [49, 69]. Additionally, we calculated the density function for 50% of the observation points of focal females and defined this area as the females’ territory .
Since individual conditions can influence female mate choice  a condition index [74–76] was calculated. The condition of each individual was determined by regression between the initial body mass of all males and females and their SVL (linear regression in males: F = 21.5; R = 0.61; R2 = 0.37; T = 4.6; N = 38; P < 0.001; y = −0.795351992 + 0.0778729712*x), (linear regression in females: R = 0.68; R2 = 0.46; N = 38; P < 0.001; y = −0.4712 + 0.0664*x). The condition of an individual was defined as the deviation of its body mass from the body mass predicted by the above regression (residuals).
We quantified the total and daily number of different contact males to determine mate search pattern of focal females in 2004 and 2005. Additionally, we compared the number of contact males on different days (oviposition day, one day before oviposition, and ≥ 2 days before oviposition) using a Friedman-ANOVA. We averaged the number of contact males for the days 2–11 before oviposition because depending on the female the time-interval between two successive ovipositions was 2–11 days.
In 2005 we also localised the sleeping places of females, allowing us to eliminate the possibility that females shift their sleeping places toward preferred males in the night before mating. Additionally, we determined whether females initiated mating inside or outside their territory. For females that initiated mating outside their territory we determined the moment when females left their territory and approached a calling male on the oviposition day. Then, we calculated the distance between the focal females and all calling males in the study area based on their actual location (position coordinates corresponded to the center of the grid cell) at that moment. For females that mated inside their territory, we calculated the distance between the current position (center of the grid cell) of the focal females and all calling males at the moment when females started approaching the male. We did this to find out whether the females mated with the closest male. It is generally unknown when the females of our study species make their decision, which in turn might depend on the moment of ovulation. However, at the day of oviposition most females stayed within their territory before directly approaching and choosing a mate (see results). Therefore, we assume that females make their decision before leaving their territories or prior to approaching a male.
With regard to male qualities we determined the condition, SVL, body weight, territory size as well as male call properties (call duration, dominant frequency, the number of pulses per call, pulse rate, call rate, duty cycle), and calling activity (the number of five minute periods of calling activity).
The focal females did not mate with a calling male more distant than the second closest to them (distance calculation at the moment when females left their territory; see above). Consequently, we compared the above mentioned qualities of these two closest males, regardless of whether females were in contact with them before or not, using a paired t-test. However, females could assess males over a longer period than just the day of oviposition. Therefore, a paired t-test demonstrated whether there was a difference between the qualities of the chosen male and the average qualities of all other contact males or the average qualities of all other males whose territories overlapped with the home range of each focal female. Females that were in contact only with the chosen male or whose home range overlapped with the territory of only one male were omitted from these analyses. A power analysis was conducted to determine if the sample size was adequate to show statistically significant differences.
Furthermore, in multivariate analyses we compared the success-rate (number of successes/total number of comparisons) of (1) all contact males and (2) all males within the home range of the female according to female choices for male traits. For this, we computed a principal component analysis (PCA) for male traits prior to correlation with success-rate because (1) many of these traits were shown to correlate with one another , (2) different studies have highlighted the importance of multiple traits in female mate choice decisions [46, 150, 151] and (3) PCA reduces the probability of inflated correlations due to excessive number of traits. Our PCA was computed with R (R Development Core Team, 2012). The parameters “territory size“, “clutch survival rate “and “calling activity“, included missing values for some males. We accounted for that by imputing missing entries using an iterative PCA algorithm on bootstrapped PCA models as described in Josse et al.  and implemented in the R package “missMDA” . Related diagnostics plots can be found in the Additional files 3, 4, 5. Success-rate was correlated with the first three dimensions of the PCA by using the “vegan” package . P-values were based on 999 permutations. The R-code for this analysis is available upon request from the authors. Finally, by using a paired t-test we compared the PCA scores between (1) the two closest calling males, (2) the chosen male and the average of all other contact males, and (3) the chosen male and the average of all other within the home range of the female. A power analysis was conducted using the R package “pwr” . If not stated otherwise, means ± SD (standard deviation) are given.
Acoustic stimuli employed in playback experiments were synthesised using the software Avisoft-SASLab Pro (Fa. Sprecht, Berlin, Germany) and were based on a typical call according to average characteristics found in the same population by Pröhl  (Figure 3). The call sequences were recorded non-stop onto high quality Sony UX-S Chrome Class Premium audio tapes using the Software Goldwave (Fa. Goldwave Inc.).
From November 2005 to May 2006 between 08:00 am and 01:00 pm hours phonotactic responses of females were tested in two-choice discrimination experiments under natural conditions in the home range of each female. Because preliminary playback tests demonstrated that non-courting females did not show positive phonotactic movements, we only tested receptive i.e. courting females. Females were considered to be courting when they approached and followed a calling male . Before each test we captured the present male as well as all other calling males within a radius of 5 m. To avoid disturbances we left the female untouched at the site. Some females were tested more than once but never for the same stimulus pair to avoid pseudo replication.
After the experiments females were captured, identified or marked and photographed, measured (SVL), and weighed. Additionally, air temperature as well as the time and duration of each test were recorded.
Call rate experiment
To ensure equal distances between both speakers as well as between the speaker and the female in each test, we located two water-resistant broadband loudspeakers (VISATON FRS 10 WP) in a triangle of 45°. The loudspeakers were placed on a 14 cm × 212 cm long plank with a distance of two metres between each other and were removable at both ends (Figure 4). The distance between both speakers was comparable to the natural distances (1–4 m) between neighbouring territorial males [52, 65, 157]. Each speaker was connected to a Sony Professional Walkman WM-D6C, an amplifier (Roadstar AM-311) and a 12V battery. Two nails and a 50 cm long stick were placed in the middle of the plank in a 90° triangle to allow adjustments of the set up.
Before each test we measured the sound pressure level (SPL) of the two speakers with the RadioShack sound level meter (measurement error: ± 2db at 114 db) and set the volume to 70 db at a distance of 50 cm . The set-up had been adjusted prior to the experimental start. The 50 cm long stick was used to obtain the correct distance of 1m between the center of the plank and the female. Both nails acted as a makeshift to ensure the correct angle of 90° between the female and both loudspeakers. Subsequently, synthetic calls were presented simultaneously from both speakers at different call rates (5 calls/sec vs. 8 calls/sec (mean of the population ± 2SD)). In accordance with other studies on dendrobatid frogs [e.g. ] we scored a positive phonotactic response when females approached a speaker within a 20 cm radius or less. Between each experiment, the side used for one stimulus was switched. A two-sided binominal test was used to determine significant preferences.
Dominant frequency experiment
The same procedure was used for the recording of dominant frequency as in the call rate experiment. However, the SPL was adjusted to 67 db at a distance of 140 cm (female start point). We presented females with two synthetic calls that differed in dominant frequency (3.6 kHz vs. 4.4 kHz (mean of the population ± 2SD); natural range of the study population: 3.2–4.5 kHz). All other call properties corresponded to the average values found in the natural population .
This experiment was conducted to test whether female mate choice is affected by the distance to potential mates. One speaker was closer (1 m) to the female than the other speaker (2.2 m) (Figure 5). The distance between both speakers was 2.3 m. The female and the speakers formed an angle of 80°. Two nails and a 66 cm long stick were placed in the middle of the plank for adjusting the set-up. The SPL of each speaker was adjusted to 70 db at a distance of 50 cm. Both speakers broadcasted identical calls with standard characteristics and a frequency of 4.4 kHz. The SPL of the closer speaker was 55 dB and the SPL of the distant speaker was 63 dB at the starting point of the female.
Andersson MB: Sexual Selection. 1994, Princeton: Princeton University Press
Hill GE: Plumage coloration is a sexually selected indicator of male quality. Nature. 1991, 350: 337-339. 10.1038/350337a0.
Doty GV, Welch AM: Advertisement call duration indicates good genes for offspring feeding rate in gray tree frogs (Hyla versicolor). Behav Ecol Sociobiol. 2001, 49: 150-156. 10.1007/s002650000291.
Hasselquist D, Bensch S, Schantz Von T: Correlation between male song repertoire, extra-pair paternity and offspring survival in the great reed warbler. Nature. 1996, 381: 229-232. 10.1038/381229a0.
Welch AM, Semlitsch D, Gerhardt HC: Call duration as an indicator of genetic quality in male gray tree frogs. Science. 1998, 280: 1928-1930. 10.1126/science.280.5371.1928.
Wilkinson GS, Presgraves DC, Crymes L: Male eye span in stalk-eyed flies indicates genetic quality by meiotic drive suppression. Nature. 1998, 391: 276-279. 10.1038/34640.
Shuster SM, Wade MJ: Mating Systems and Strategies. 2003, Princeton: Princeton University Press
Gibson RM, Langen TA: How do animals choose their mates?. Trends Ecol Evol. 1996, 11: 468-470. 10.1016/0169-5347(96)10050-1.
Jennions MD, Petrie M: Variation in mate choice and mating preferences: a review of causes and consequences. Biol Rev. 1997, 72: 283-327. 10.1017/S0006323196005014.
Janetos AC: Strategies of female mate choice: a theoretical analysis. Behav Ecol Sociobiol. 1980, 7: 107-112. 10.1007/BF00299515.
Real L: Search theory and mate choice. I. Models of single-sex discrimination. Am Nat. 1990, 136: 376-405. 10.1086/285103.
Seger J: Unifying genetic models for the evolution of female choice. Evolution. 1985, 39: 1185-1193. 10.2307/2408776.
Dombrovsky Y, Perrin N: On adaptive search and optimal stopping in sequential mate choice. Am Nat. 1994, 144: 355-361. 10.1086/285680.
Luttbeg B: A Comparative Bayes tactic for mate assessment and choice. Behav Ecol. 1996, 7: 451-460. 10.1093/beheco/7.4.451.
Luttbeg B: Assessing the robustness and optimality of alternative decision rules with varying assumptions. Anim Behav. 2002, 63: 805-814. 10.1006/anbe.2001.1979.
Wiegmann DD, Mukhopadhyay K, Real LA: Sequential search and the influence of male quality on female mating decisions. J Math Biol. 1999, 39: 193-216. 10.1007/s002850050168.
Wiegmann DD, Real LA, Capone TA, Ellner S: Some distinguishing features of models of search behavior and mate choice. Am Nat. 1996, 147: 188-204. 10.1086/285846.
Wiegmann DD, Angeloni LM, Seubert SM, Wade JG: Mate choice decisions by searchers. Current Zoology. 2013, 59: 184-199.
Reid ML, Stamps JA: Female mate choice tactics in a resource-based mating system: field tests of alternative models. Am Nat. 1997, 150: 98-121. 10.1086/286058.
Downhower JF, Lank DB: Effect of previous experience on mate choice by female mottled sculpins. Anim Behav. 1994, 47: 369-372. 10.1006/anbe.1994.1050.
Choudhury S, Black JM: Mate-selection behaviour and sampling strategies in geese. Anim Behav. 1993, 46: 747-757. 10.1006/anbe.1993.1252.
Wittenberger JF: Tactics of mate choice. Mate choice. Edited by: Bateson P. 1983, Cambridge: Cambridge University Press
Moore AJ, Moore PJ: Female strategy during mate choice: threshold assessment. Evolution. 1988, 42: 387-391. 10.2307/2409241.
Murphy CG, Gerhardt HC: Mate sampling by female barking treefrogs (Hyla gratiosa). Behav Ecol. 2002, 13: 472-480. 10.1093/beheco/13.4.472.
Beckers OM, Wagner WE: Mate sampling strategy in a field cricket: evidence for a fixed threshold strategy with last chance option. Anim Behav. 2011, 81: 519-527. 10.1016/j.anbehav.2010.11.022.
Ivy TM, Sakaluk SK: Sequential mate choice in decorated crickets: females use a fixed internal threshold in pre- and postcopulatory choice. Anim Behav. 2007, 74: 1065-1072. 10.1016/j.anbehav.2007.01.017.
Dale S, Rinden H, Slagsvold T: Competition for a mate restricts mate search of female pied flycatchers. Behav Ecol Sociobiol. 1992, 30: 165-176. 10.1007/BF00166699.
Byers JA, Byers AA, Dunn SJ: A dry summer diminishes mate search effort by pronghorn females: evidence for a significant cost of mate search. Ethology. 2006, 112: 74-80. 10.1111/j.1439-0310.2006.01127.x.
Backwell PRY, Passmore NI: Time constraints and multiple choice criteria in the sampling behaviour and mate choice of the fiddler crab, Uca annulipes. Behav Ecol Sociobiol. 1996, 38: 407-416. 10.1007/s002650050258.
Johnson JB, Basolo AL: Predator exposure alters female mate choice in the green swordtail. Behav Ecol. 2003, 14: 619-625. 10.1093/beheco/arg046.
Willis PMP, Rosenthal GGG, Ryan MJM: An indirect cue of predation risk counteracts female preference for conspecifics in a naturally hybridizing fish Xiphophorus birchmanni. PLoS ONE. 2011, 7: e34802-e34802.
Collins SA: The effect of recent experience on female choice in zebra finches. Anim Behav. 1995, 49: 479-486. 10.1006/anbe.1995.0062.
Kvarnemo C, Forsgren E: The influence of potential reproductive rate and variation in mate quality on male and female choosiness in the sand goby, Pomatoschistus minutus. Behav Ecol Sociobiol. 2000, 48: 378-384. 10.1007/s002650000246.
Kvarnemo C, Simmons LW: Variance in female quality, operational sex ratio and male mate choice in a bushcricket. Behav Ecol Sociobiol. 1999, 45: 245-252. 10.1007/s002650050559.
Rosenqvist G, Houde A: Prior exposure to male phenotypes influences mate choice in the guppy, Poecilia reticulata. Behav Ecol. 1997, 8: 194-198. 10.1093/beheco/8.2.194.
Parker GA: Mate quality and mating decisions. Mate choice. Edited by: Bateson P. 1983, Cambridge: Cambridge University Press, 141-166.
Fawcett TW, Johnstone RA: Mate choice in the face of costly competition. Behav Ecol. 2003, 14: 771-779. 10.1093/beheco/arg075.
Gibson RM: Female choice in sage grouse: the roles of attraction and active comparison. Behav Ecol Sociobiol. 1996, 39: 55-59. 10.1007/s002650050266.
Rintamäki PT, Alatalo RV, Höglund J, Lundberg A: Mate sampling behaviour of black grouse females (Tetrao tetrix). Behav Ecol Sociobiol. 1995, 37: 209-215. 10.1007/BF00176719.
Gerhardt HC, Huber F: Acoustic Communication in Insects and Anurans. 2002, Chicago: University of Chicago Press
Otter KA, Ratcliffe L: Enlightened decisions: female assessment and communication networks. 2005, Cambridge: Cambridge University Press: In Animal communication networks. Edited by McGregor PK, 133-151. 2005
Ryan MJ: The Túngara Frog, a Study in Sexual Selection and Communication. 1985, Chicago: University of Chicago Press
Murphy CG: Simultaneous mate-sampling by female barking treefrogs (Hyla gratiosa). Behav Ecol. 2012, 10.1093/beheco/ars093. Advanced access
Schwartz JJ, Huth K, Hutchin T: How long do females really listen? Assessment time for female mate choice in the grey treefrog, Hyla versicolor. Anim Behav. 2004, 68: 533-540. 10.1016/j.anbehav.2003.09.016.
Gerhardt HC, Klump GM: Masking of acoustic signals by the chorus background noise in the green tree frog: A limitation on mate choice. Anim Behav. 1988, 36: 1247-1249. 10.1016/S0003-3472(88)80090-3.
Poole K, Murphy CG: Preferences of female barking treefrogs, Hyla gratiosa, for larger males: univariate and composite tests. Anim Behav. 2007, 73: 513-524. 10.1016/j.anbehav.2006.09.008.
Wells KD: The Ecology and Behavior of Amphibians. 2007, Chicago: University of Chicago Press
Meuche I, Pröhl H: Alternative mating tactics in the strawberry poison frog (Oophaga pumilio). Herpetol J. 2011, 21: 275-277.
Meuche I, Linsenmair KE, Pröhl H: Intrasexual competition, territoriality and acoustic communication in male strawberry poison frogs (Oophaga pumilio). Behav Ecol Sociobiol. 2012, 66: 613-621. 10.1007/s00265-011-1309-9.
Baugh JR, Forester DC: Prior residence effect in the dart-poison frog, Dendrobates pumilio. Behaviour. 1994, 131: 207-224. 10.1163/156853994X00442.
Pröhl H: Variation in male calling behaviour and relation to male mating success in the strawberry poison frog (Dendrobates pumilio). Ethology. 2003, 109: 273-290. 10.1046/j.1439-0310.2003.00863.x.
Bee MA: A test of the “dear enemy effect” in the strawberry dart-poison frog (Dendrobates pumilio). Behav Ecol Sociobiol. 2003, 54: 601-610. 10.1007/s00265-003-0657-5.
Gardner EA, Graves BM: Responses of resident male Dendrobates pumilio to territory intruders. J Herpetol. 2005, 39: 248-253. 10.1670/95-04A.
Pröhl H: Territorial behaviour of the strawberry poison-dart frog, Dendrobates pumilio. Amphibia-Reptilia. 1997, 18: 437-442. 10.1163/156853897X00495.
Pröhl H, Berke O: Spatial distributions of male and female strawberry poison frogs and their relation to female reproductive resources. Oecologia. 2001, 129: 534-542.
Stynoski JL: Discrimination of offspring by indirect recognition in an egg-feeding dendrobatid frog, Oophaga pumilio. Anim Behav. 2009, 78: 1351-1356. 10.1016/j.anbehav.2009.09.002.
Stynoski JL, Noble VR: To beg or to freeze: multimodal sensory integration directs behavior in a tadpole. Behav Ecol Sociobiol. 2012, 66: 191-199. 10.1007/s00265-011-1266-3.
Maan ME, Cummings ME: Sexual dimorphism and directional sexual selection on aposematic signals in a poison frog. Proc Natl Acad Sci USA. 2009, 106: 19072-19077. 10.1073/pnas.0903327106.
Richards-Zawacki CL, Cummings ME: Intraspecific reproductive character displacement in a polymorphic poison dart frog, Dendrobates pumilio. Evolution. 2011, 65: 259-267. 10.1111/j.1558-5646.2010.01124.x.
Richards-Zawacki CL, Wang IJ, Summers K: Mate choice and the genetic basis for colour variation in a polymorphic dart frog: inferences from a wild pedigree. Mol Ecol. 2012, 21: 3879-3892. 10.1111/j.1365-294X.2012.05644.x.
Brown JL, Maan ME, Cummings ME, Summers K: Evidence for selection on coloration in a Panamanian poison frog: a coalescent-based approach. J Biogeogr. 2010, 37: 891-901. 10.1111/j.1365-2699.2009.02260.x.
Wang IJ, Summers K: Genetic structure is correlated with phenotypic divergence rather than geographic isolation in the highly polymorphic strawberry poison-dart frog. Mol Ecol. 2010, 19: 447-458. 10.1111/j.1365-294X.2009.04465.x.
Rudh AA, Rogell BB, Håstad OO, Qvarnström AA: Rapid population divergence linked with co-variation between coloration and sexual display in strawberry poison frogs. Evolution. 2011, 65: 1271-1282. 10.1111/j.1558-5646.2010.01210.x.
Crothers LR, Cummings ME: Warning signal brightness variation: Sexual selection may work under the radar of natural selection in populations of a polytypic poison frog. Am Nat. 2013, 181: E000-
Pröhl H, Hödl W: Parental investment, potential reproductive rates, and mating system in the strawberry dart-poison frog, Dendrobates pumilio. Behav Ecol Sociobiol y. 1999, 46: 215-220. 10.1007/s002650050612.
Limerick S: Courtship behavior and oviposition of the poison-arrow frog Dendrobates pumilio. Herpetologica. 1980, 36: 69-71.
Weygoldt P: Complex brood care and reproductive behavior in captive poison-arrow frogs, Dendrobates pumulio O. Schmidt. Behav Ecol Sociobiol. 1980, 7: 329-332. 10.1007/BF00300674.
Weygoldt P: Evolution of parental care in dart poison frogs (Amphibia: Anura: Dendrobatidae). J Zool Syst Evol Res. 1987, 25: 51-67.
Meuche I, Linsenmair KE, Pröhl H: Female territoriality in the strawberry poison frog (Oophaga pumilio). Copeia. 2011, 2011: 351-356. 10.1643/CE-08-135.
Bakker TCM, Milinski M: Sequential female choice and the previous male effect in sticklebacks. Behav Ecol Sociobiol. 1991, 29: 205-210. 10.1007/BF00166402.
Shuker DM, Day TH: Mate sampling and the sexual conflict over mating in seaweed flies. Behav Ecol. 2002, 13: 83-86. 10.1093/beheco/13.1.83.
Forsgren E: Mate sampling in a population of sand gobies. Anim Behav. 1997, 53: 267-276. 10.1006/anbe.1996.0374.
Dale S, Amundsen T, Lifjeld JT, Slagsvold T: Mate sampling behavior of female pied flycatchers: evidence for active mate choice. Behav Ecol Sociobiol. 1990, 27: 87-91. 10.1007/BF00168450.
Jakob EM, Marshall SD, Uetz GW: Estimating fitness: a comparison of body condition indices. Oikos. 1996, 77: 61-67. 10.2307/3545585.
Murphy CG: Determinants of chorus tenure in barking treefrogs (Hyla gratiosa). Behav Ecol Sociobiol. 1994, 34: 285-294. 10.1007/BF00183479.
Meuche I, Grafe TU: Supplementary feeding affects the breeding behaviour of male European treefrogs (Hyla arborea). BMC Ecol. 2009, 9: 1-10.1186/1472-6785-9-1.
Ryan MJ, Keddy-Hector A: Directional patterns of female mate choice and the role of sensory biases. Am Nat. 1992, 139: S4-S35. 10.1086/285303.
Howard RD, Young JR: Individual variation in male vocal traits and female mating preferences in Bufo americanus. Anim Behav. 1998, 55: 1165-1179. 10.1006/anbe.1997.0683.
Lopez PT, Narins PM, Lewis ER, Moore SW: Acoustically induced call modification in the white-lipped frog, Leptodactylus albilabris. Anim Behav. 1988, 36: 1295-1308. 10.1016/S0003-3472(88)80198-2.
Wagner WE: Graded aggressive signals in Blanchard's cricket frog: vocal responses to opponent proximity and size. Anim Behav. 1989, 38: 1025-1038. 10.1016/S0003-3472(89)80141-1.
Davies NB, Halliday TR: Deep croaks and fighting assessment in toads Bufo bufo. Nature. 1978, 274: 683-685. 10.1038/274683a0.
Grafe TUT, Preininger DD, Sztatecsny MM, Kasah RR, Dehling JMJ, Proksch SS, Hödl WW: Multimodal communication in a noisy environment: a case study of the Bornean Rock Frog Staurois parvus. PLoS ONE. 2012, 7: e37965-e37965. 10.1371/journal.pone.0037965.
Grafe TU, Wanger TC: Multimodal signaling in male and female foot-flagging frogs Staurois guttatus (Ranidae): An alerting function of calling. Ethology. 2007, 113: 772-781. 10.1111/j.1439-0310.2007.01378.x.
Hirschmann W, Hödl W: Visual signaling in Phrynobatrachus krefftii Boulenger, 1909 (Anura: Ranidae). Herpetologica. 2006, 62: 18-27. 10.1655/04-06.1.
Amézquita A, Hödl W: How, when, and where to perform visual displays: the case of the Amazonian frog Hyla parviceps. Herpetologica. 2004, 60: 420-429. 10.1655/02-51.
de Luna AG, Hoedl W, Amézquita A: Colour, size and movement as visual subcomponents in multimodal communication by the frog Allobates femoralis. Anim Behav. 2010, 79: 739-745. 10.1016/j.anbehav.2009.12.031.
Gomez D, Richardson C, Lengagne T, Plenet S, Joly P, Lena J-P, Thery M: The role of nocturnal vision in mate choice: females prefer conspicuous males in the European tree frog (Hyla arborea). Proc R Soc Lond B. 2009, 276: 2351-2358. 10.1098/rspb.2009.0168.
Ptacek MB: The role of mating preferences in shaping interspecific divergence in mating signals in vertebrates. Behav Process. 2000, 51: 111-134. 10.1016/S0376-6357(00)00123-6.
Summers K, Symula R, Clough M, Cronin T: Visual mate choice in poison frogs. Proc R Soc Lond B. 1999, 266: 2141-2145. 10.1098/rspb.1999.0900.
Reynolds RGR, Fitzpatrick BMB: Assortative mating in poison-dart frogs based on an ecologically important trait. Evolution. 2007, 61: 2253-2259. 10.1111/j.1558-5646.2007.00174.x.
Maan ME, Cummings ME: Female preferences for aposematic signal components in a polymorphic poison frog. Evolution. 2008, 62: 2334-2345. 10.1111/j.1558-5646.2008.00454.x.
Daly JW, Myers CW: Toxicity of Panamanian poison frogs (Dendrobates): some biological and chemical aspects. Science. 1967, 156: 970-973. 10.1126/science.156.3777.970.
Summers K, Cronin TW, Kennedy T: Variation in spectral reflectance among populations of Dendrobates pumilio, the strawberry poison frog, in the Bocas del Toro Archipelago, Panama. J Biogeogr. 2003, 30: 35-53. 10.1046/j.1365-2699.2003.00795.x.
Sullivan BK, Ryan MJ, Verrell PA: Female choice and mating system structure. Amphibian biology. Volume 2. Social Behavior. Edited by: Heatwole H, Sullivan BK. 1995, Chipping Norton: Surrey Beatty and Sons, 469-517.
Howard RD: The evolution of mating strategies in bullfrogs, Rana catesbeiana. Evolution. 1978, 32: 850-871. 10.2307/2407499.
Howard RD: The Influence of Male Defended Oviposition Sites on Embryo Mortality in Bullfrogs. Ecology. 1978, 59: 789-798. 10.2307/1938783.
Forsgren E: Female sand gobies prefer good fathers over dominant males. Proc R Soc Lond B, Volume 264. 1997, 1283-1286.
Kluge AG: The Life History, Social Organization, and Parental Behavior of Hyla Rosenbergi Boulenger, a Nest-Building Gladiator Frog. 1981, University of Michigan: Museum of Zoology
Pröhl H: Clutch loss affects the operational sex ratio in the strawberry poison frog Dendrobates pumilio. Behav Ecol Sociobiol. 2005, 58: 310-315. 10.1007/s00265-005-0915-9.
Pröhl H: Territorial behavior in dendrobatid frogs. J Herpetol. 2005, 39: 354-365. 10.1670/162-04A.1.
Alonzo SH, Heckman KL: The unexpected but understandable dynamics of mating, paternity and paternal care in the ocellated wrasse. Proc R Soc Lond B. 2010, 277: 115-122. 10.1098/rspb.2009.1425.
Alonzo SH: Uncertainty in territory quality affects the benefits of usurpation in a Mediterranean wrasse. Behav Ecol. 2004, 15: 278-285. 10.1093/beheco/arh007.
Arak A: Female mate selection in the natterjack toad: active choice or passive atraction?. Behav Ecol Sociobiol y. 1988, 22: 317-327.
Grafe TU: Costs and benefits of mate choice in the lek-breeding reed frog, Hyperolius marmoratus. Anim Behav. 1997, 53: 1103-1117. 10.1006/anbe.1996.0427.
Ursprung E, Ringler M, Jehle R, Hoedl 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-1771. 10.1111/j.1365-294X.2011.05056.x.
Ringler EE, Ringler MM, Jehle RR, Hödl WW: The female perspective of mating in A. femoralis, a territorial frog with paternal care - a spatial and genetic analysis. PLoS ONE. 2012, 7: e40237-e40237. 10.1371/journal.pone.0040237.
Uy JA, Patricelli GL, Borgia G: Complex mate searching in the satin bowerbird Ptilonorhynchus violaceus. Am Nat. 2001, 158: 530-542. 10.1086/323118.
Uy JA, Patricelli GL, Borgia G: Dynamic mate-searching tactic allows female satin bowerbirds Ptilonorhynchus violaceus to reduce searching. Proc R Soc Lond B. 2000, 267: 251-256. 10.1098/rspb.2000.0994.
Pröhl H: Population differences in female resource abundance, adult sex ratio, and male mating success in Dendrobates pumilio. Behav Ecol. 2002, 13: 175-181. 10.1093/beheco/13.2.175.
Slagsvold T, Dale S: Mate choice models: can cost of searching and cost of courtship explain mating patterns of female Pied Flycatchers?. Ornis Scand. 1991, 22: 319-326. 10.2307/3676503.
Alatalo RV, Carlson A, Lundberg A: The search cost in mate choice of the pied flycatcher. Anim Behav. 1988, 36: 289-291. 10.1016/S0003-3472(88)80272-0.
Hedrick AV, Dill LM: Assessing female choice for mates when the males' characters vary during the sampling period. Anim Behav. 1990, 40: 780-782. 10.1016/S0003-3472(05)80709-2.
Sullivan MS: Mate choice as an information gathering process under time constraint: implications for behaviour and signal design. Anim Behav. 1994, 47: 141-151. 10.1006/anbe.1994.1016.
Lynch K, Rand AS, Ryan MJ, Wilczynski W: Plasticity in female mate choice associated with changing reproductive states. Anim Behav. 2005, 69: 689-699. 10.1016/j.anbehav.2004.05.016.
Castellano S, Cermelli P: Sampling and assessment accuracy in mate choice: a random-walk model of information processing in mating decision. J Theor Biol. 2011, 274: 161-169. 10.1016/j.jtbi.2011.01.001.
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.
Berglund A: The operational sex ratio influences choosiness in a pipefish. Behav Ecol. 1994, 5: 254-258. 10.1093/beheco/5.3.254.
Berglund A: Many mates make male pipefish choosy. Behaviour. 1995, 132: 213-218. 10.1163/156853995X00702.
Jirotkul M: Operational sex ratio influences female preference and male-male competition in guppies. Anim Behav. 1999, 58: 287-294. 10.1006/anbe.1999.1149.
Vitousek MN, Mitchell MA, Woakes AJ, Niemack MD, Wikelski M: Does competition allow male mate choosiness in threespine sticklebacks?. Am Nat. 2009, 173: e567-e567.
Pröhl H: Patrón reproductivo en Dendrobates pumilio (Anura: Dendrobatidae). Rev Biol Trop. 1997, 45: 1669-1676.
Vitousek MN, Mitchell MA, Woakes AJ, Niemack MD: High costs of female choice in a lekking lizard. PLoS ONE. 2007, 2 (6): e567-10.1371/journal.pone.0000567.
Milinski M, Bakker TCM: Costs influence sequential mate choice in sticklebacks, Gasterosteus aculeatus. Proc R Soc Lond B. 1992, 250: 229-233. 10.1098/rspb.1992.0153.
Slagsvold T, Lifjeld JT, Stenmark G, Breiehagen T: On the cost of searching for a mate in female pied flycatchers Ficedula hypoleuca. Anim Behav. 1988, 36: 433-442. 10.1016/S0003-3472(88)80013-7.
Poulin R: Mate choice decisions by parasitized female upland bullies, Gobiomorphus breviceps. Proc R Soc Lond B. 1994, 256: 183-187. 10.1098/rspb.1994.0068.
Pfennig KS, Tinsley RC: Different mate preferences by parasitized and unparasitized females potentially reduces sexual selection. J Evol Biol. 2002, 15: 399-406. 10.1046/j.1420-9101.2002.00406.x.
López S: Parasitized female guppies do not prefer showy males. Anim Behav. 1999, 57: 1129-1134. 10.1006/anbe.1998.1064.
Buchholz R: Effects of parasitic infection on mate sampling by female wild turkeys (Meleagris gallopavo): should infected females be more or less choosy?. Behav Ecol. 2004, 15: 687-694. 10.1093/beheco/arh066.
Martinez SA, Maggenti AR: Cosmocerca panamensis sp. N. (Nematoda: Cosmocercidae) from the Panamanian Poison-arrow Frog, Dendrobates pumilio Schmidt, 1857, with a discussion of prodelphy, the type species and family authorship in Cosmocerca Diesing, 1861. Proc Helminthol Soc Wash, Volume 56. 1989, 97-103.
Bursey CR, Brooks DR: Nematode parasites of 41 anuran species from the area de conservación Guanacaste, Costa Rica. Comp Parasitol. 2010, 77: 221-231. 10.1654/4418.1.
Pröhl H, Eulenburg J, Meuche I, Bolaños F: Parasite infection has little effect on sexual signals and reproductive behaviour in strawberry poison frogs. Evol Ecol. in press
Goldberg SR, Bursey CR: Helminth biodiversity of Costa Rican Anurans (Amphibia: Anura). J Nat Hist. 2010, 44: 1755-1787. 10.1080/00222931003764931.
Forsgren E: Predation risk affects mate choice in a gobiid fish. Am Nat. 1992, 140: 1041-1049. 10.1086/285455.
Hedrick AV, Dill LM: Mate choice by female crickets is influenced by predation risk. Anim Behav. 1993, 46: 193-196. 10.1006/anbe.1993.1176.
Friedl TWP, Klump GM: Sexual selection in the lek-breeding European treefrog (Hyla arborea): body size, chorus attendance, random mating and good genes. Anim Behav. 2005, 70: 1141-1154. 10.1016/j.anbehav.2005.01.017.
Ryan MJ, Tuttle MD, Taft LK: The costs and benefits of frog chorusing behavior. Behav Ecol Sociobiol. 1981, 8: 273-278. 10.1007/BF00299526.
Saporito RA, Zuercher R, Roberts M, Gerow KG, Donnelly MA: Experimental evidence for aposematism in the dendrobatid poison frog Oophaga pumilio. Copeia. 2007, 2007: 1006-1011. 10.1643/0045-8511(2007)7[1006:EEFAIT]2.0.CO;2.
Hegna RH, Saporito RA, Donnelly MA: Not all colors are equal: predation and color polytypism in the aposematic poison frog Oophaga pumilio. Evol Ecol. 2012, 10.1007/s10682-012-9605-z.
Maan ME, Cummings ME: Poison frog colors are honest signals of toxicity, particularly for bird predators. Am Nat. 2012, 179: E1-14. 10.1086/663197.
Saporito RA, Donnelly MA, Garraffo HM, Spande TF, Daly JW: Geographic and seasonal variation in alkaloid-based chemical defenses of Dendrobates pumilio from Bocas del Toro, Panama. J Chem Ecol. 2006, 32: 795-814. 10.1007/s10886-006-9034-y.
Saporito RA, Donnelly MA, Jain P, Martin Garraffo H, Spande TF, Daly JW: Spatial and temporal patterns of alkaloid variation in the poison frog Oophaga pumilio in Costa Rica and Panama over 30 years. Toxicon. 2007, 50: 757-778. 10.1016/j.toxicon.2007.06.022.
Klump GM, Gerhardt HC: Use of non-arbitrary acoustic criteria in mate choice by female gray tree frogs. Nature. 1987, 326: 286-288. 10.1038/326286a0.
Martin P, Bateson P: Measuring Behaviour: an Introductory Guide. 2007, Cambridge: Cambridge University Press, 3
Ursprung E, Ringler M, Jehle R, Hoedl W: Toe regeneration in the neotropical frog Allobates femoralis. Herpetol J. 2011, 21: 83-86.
Meuche I: Changes of individual colour patterns in the Central American strawberry poison frog, Oophaga pumilio (Amphibia: Dendrobatidae). Salamandra. 2009, 45: 177-179.
Gerhardt HC: Female mate choice in treefrogs: static and dynamic acoustic criteria. Anim Behav. 1991, 42: 615-635. 10.1016/S0003-3472(05)80245-3.
Taigen TL, Wells KD: Energetics of vocalization by an anuran amphibian (Hyla versicolor). J Comp Physiol B. 1985, 155: 163-170. 10.1007/BF00685209.
Worton BJ: Kernel methods for estimating the utilization distribution in home-range studies. Ecology. 1989, 70: 164-168. 10.2307/1938423.
Rodgers AR, Carr AP: HRE: the home range extension for ArcView. User’s manual. 1998, Thunder Bay, Ontario, Canada: Centre for Northern Forest Ecosystem Research, Ontario Ministry of Natural Resources
Candolin U: The use of multiple cues in mate choice. Biol Rev. 2003, 78: 575-595. 10.1017/S1464793103006158.
Hamilton PS, Sullivan BK: Female mate attraction in ornate tree lizards, Urosaurus ornatus: a multivariate analysis. Anim Behav. 2005, 69: 219-224. 10.1016/j.anbehav.2004.03.011.
Josse J, Husson F, Pagès J: Gestion des données manquantes en analyse en composantes principales. Journal de la SFdS. 2009, 150: 28-51.
missMDA: Handling missing values with/in multivariate data analysis (principal component methods). 2012, R package version 1.4. URL http://www.agrocampus-ouest.fr/math/husson
Oksanen J, Blanchet FG, Kindt R, Legendre P, O'Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H: vegan: Community Ecology Package. 2013, R package version 2.0-7. http://cran.r-project.org/web/packages/vegan/index.html
Champely S: pwr: Basic functions for power analysis. R package version 1.1.1. 2009,http://cran.r-project.org/web/packages/pwr/index.html,
Donnelly M: Demographic effects of reproductive ressource supplementation in a territorial frog - Dendrobates pumilio. Ecol Monogr. 1989, 59: 207-221. 10.2307/1942599.
Bunnell P: Vocalizations in the territorial behavior of the frog Dendrobates pumilio. Copeia. 1973, 1973: 277-284. 10.2307/1442966.
Ursprung E, Ringler M, Hödl W: Phonotactic approach pattern in the neotropical frog Allobates femoralis: a spatial and temporal analysis. Behaviour. 2009, 146: 153-170. 10.1163/156853909X410711.
We wish to thank the government of Costa Rica (MINAE; Servicio de Parques Nacionales) for permission to conduct the research in the Hitoy Cerere Biological Reserve, and especially the staff at the Biological Reserve for their support. For cooperation and/or discussion, we wish to thank Federico Bolaños, Ulmar Grafe, Marine Joly, Astrid Krug, Frank Scherbarth, Sabine Schmidt, Aline Horwath and Sönke von den Berg. For invaluable assistance in the field we would like to thank Susanne Becker, Ulrike Freisen, Franziska Grözinger, Ulrike Hartmann, Angelika Kronhardt, Björn Marschall, Karin Moll, Esther Müller, Claudia Ramenda, Isabelle Schulz and Jairo Reyes Dominguez. The German Academic Exchange Service (DAAD) and the German Research Foundation (DFG) supported this work.
The authors declare that they have no competing interests.
IM, OB, AK, KEL and HP participated in design of the study and/or statistical analysis. IM and AK performed the data collection. IM wrote the manuscript. OB, AK, KEL and HP helped to draft the manuscript. All authors read and approved the final manuscript.