International Journal of Evolutionary BiologyVolume 2012 (2012), Article ID 148745, 14 pagesdoi:10.1155/2012/148745

Research Article

Divergent Evolution of Male Aggressive Behaviour: Another Reproductive Isolation Barrier in Extremophile Poeciliid Fishes?

David Bierbach,¹ Moritz Klein,¹ Vanessa Sassmannshausen,¹ Ingo Schlupp,² Rüdiger Riesch,³ Jakob Parzefall,⁴ and Martin Plath¹

¹Department of Ecology and Evolution, J.W. Goethe University Frankfurt, Siesmayerstrasse 70-72, 60054 Frankfurt am Main, Germany²Department of Zoology, University of Oklahoma, 730 Van Vleet Oval, Norman, OK 73019, USA³Department of Biology & W. M. Keck Center for Behavioral Biology, North Carolina State University, 127 David Clark Labs, Raleigh, NC 27695-7617, USA⁴Department of Ethology, University of Hamburg, Martin-Luther-King-Platz 3, 20146 Hamburg, Germany

Received 15 July 2011; Revised 20 September 2011; Accepted 14 October 2011

Academic Editor: Kyoichi Sawamura

Copyright © 2012 David Bierbach et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Reproductive isolation among locally adapted populations may arise when immigrants from foreign habitats are selected against via natural or (inter-)sexual selection (female mate choice). We asked whether also intrasexual selection through male-male competition could promote reproductive isolation among populations of poeciliid fishes that are locally adapted to extreme environmental conditions [i.e., darkness in caves and/or toxic hydrogen sulphide (H₂S)]. We found strongly reduced aggressiveness in extremophile Poecilia mexicana, and darkness was the best predictor for the evolutionary reduction of aggressiveness, especially when combined with presence of H₂S. We demonstrate that reduced aggression directly translates into migrant males being inferior when paired with males from nonsulphidic surface habitats. By contrast, the phylogenetically old sulphur-endemic P. sulphuraria from another sulphide spring area showed no overall reduced aggressiveness, possibly indicating evolved mechanisms to better cope with H₂S.

1. Introduction

1.1. Ecological Speciation

Divergent natural selection has the potential to drive adaptive trait divergence along environmental gradients [1], but can also lead to the evolution of reproductive isolating barriers [2, 3]. During ecological speciation, reproductive isolation results from ecologically based divergent selection, and prezygotic isolation may arise as a byproduct of local adaptation if immigrants from ecologically divergent habitats are selected against [3]. This can be owing to natural selection, if immigrants show reduced viability [4–6], or sexual selection, if poorly adapted individuals have a disadvantage in mate competition [5, 7, 8]. Furthermore, ecological speciation may also be driven by selection against hybrids with intermediate phenotypes [9], behavioural isolation based on a “magic trait” [10–14], and sensory drive [15].

Our present paper briefly collates our current knowledge regarding trait divergence and especially mechanisms of reproductive isolation among different locally adapted populations of livebearing fishes (Poeciliidae), currently undergoing ecological speciation processes in response to “extreme” conditions (see below). Using both lab-reared as well as wild-caught fish we then demonstrate that divergent evolution of male competitive abilities (aggressive behaviour) in extremophile fishes may play yet another role in maintaining reproductive isolation among different locally adapted populations: adaptation to extreme habitat conditions appears to have selected for reduced aggressiveness, and we show that this renders potential migrant males from extreme habitats less competitive in intrasexual combat when the resident males inhabiting benign habitats show “normal” aggressive behaviour. As the mating system of our study species is based on male dominance hierarchies, with dominant males aggressively defending small shoals of females from intruders [16–18], we argue that this pattern directly translates into reproductive inferiority of such migrant males.

1.2. Life in Extreme Habitats

Habitats can be considered extreme if certain characteristics of the environment are outside of the range normally experienced by a species and if organisms colonizing this particular habitat type experience an initial reduction in fitness [19, 20]. For example, some extreme environments are characterized by exceptionally high concentrations of hydrogen sulphide (H₂S): deep-sea hydrothermal vents, hydrocarbon seeps, as well as intertidal zones, salt marshes, mudflats, and sewage outfalls, where H₂S is usually of biogenic origin [21–24]. H₂S inhibits aerobic respiration due to its interference with mitochondrial respiration and blood oxygen transport, but also leads to extreme hypoxia in the water [21, 22]. This makes H₂S acutely toxic to most metazoans even in micromolar amounts, and accordingly, pulses of H₂S discharge have been reported to be the source of mass mortalities [22].

An environmental toxicant like hydrogen sulphide that requires energetically costly behavioural (i.e., actively avoiding microhabitats with high levels of toxicity) and physiological adaptations (various forms of detoxification) by animals exposed to it will certainly have a profound influence on the evolutionary trajectories of populations experiencing the stressor [6]. For instance, when exposed to H₂S and hypoxia, livebearing fishes resort to aquatic surface respiration (ASR) and, thus, exploit the more oxygen-rich air-water interface [25]. Under experimental conditions, Atlantic mollies (Poecilia mexicana) have been shown to spend more than 60% of their time performing ASR when exposed to sulphidic water [25], and in natural populations P. mexicana have been observed to spend up to 84% of their time performing ASR [26]. However, while access to the water surface (i.e., the possibility to perform ASR) is a strong predictor of short-term survival in fish exposed to H₂S-containing water [25], time spent at the water surface clearly trades off against the time fish can spend foraging. Hence, fish from H₂S-containing habitats tend to have less food in their guts and lower body condition than conspecifics from nonsulphidic habitats [26–29].

Beside toxicants, perpetual darkness, like in cave ecosystems, can represent an extreme condition for typical surface-dwelling organisms like P. mexicana [30]. Darkness renders visual orientation and navigation an impossible task, and cave organisms need to develop specific adaptations to cope with this situation [31–34]. Cave animals (especially crustaceans and fishes) are widely used model organisms to study the evolutionary effects of permanent darkness on various traits, including improved nonvisual sensory systems and increased starvation tolerance (e.g., [34–39]).

1.3. Ecological Speciation in Extremophile Poeciliid Fishes

Notwithstanding all the adverse effects of H₂S, several species of livebearing fishes (Poeciliidae) have been documented to thrive (and speciate) in waters containing exorbitant concentrations of H₂S. Among them are sulphur endemics like the sulphur molly (Poecilia sulphuraria) and widemouth gambusia (Gambusia eurystoma) [8, 27, 40], as well as species that are currently undergoing ecological speciation, like certain populations of P. mexicana [5, 6, 28, 41–43].

Of particular interest are different locally adapted P. mexicana populations in the Cueva del Azufre system (Tabasco, Mexico), a system that is characterized by the simultaneous action of two strong selective forces: permanent darkness in subterranean parts of streams and toxic H₂S [6, 30, 44, 45] of volcanic origin [46–48]. Within a small geographic range of only few kilometres, reproductively isolated populations of P. mexicana inhabit environments characterized by all possible combinations of these two factors: a toxic cave (Cueva del Azufre, CA), a nontoxic cave (Cueva Luna Azufre, LA), and toxic surface waters; however, a small cascade separates all extreme habitats from nonsulphidic, normoxic sites (for discussion see [42]).

Another system considered in our present study is the sulphur molly system situated at the Baños del Azufre near Teapa (Tabasco, Mexico). This system is characterized by even higher H₂S concentrations (around 230 μM [8, 43]). Just like in the Cueva del Azufre system, no barriers other than presence of environmental stressors prevent movement of fish among different habitat types in this system [43]. P. sulphuraria forms a monophyletic sister clade with phylogenetic affinity to a northern clade of P. mexicana rather than P. mexicana inhabiting the clear-water habitats in the vicinity of the Baños del Azufre [6]. Thus, sulphur mollies appear to represent a phylogenetically old sulphur-adapted lineage and have been considered a potential “endpoint” of H₂S adaptation [27].

Extremophile P. mexicana in the Cueva del Azufre system are characterized by site-specific local adaptations in several behavioural (e.g., [25, 26, 49–51]), dietary [52], female and male life-history [27–29, 53, 54], morphological [17, 43, 45, 55, 56], and physiological traits [17, 57], and there is strong evidence for convergent patterns of H₂S adaptations across both aforementioned sulphur systems [6, 27].

1.4. Reproductive Isolating Barriers in Extremophile Poeciliids

Gene flow between populations with different ecological backgrounds in the Cueva del Azufre system is virtually absent with the exception of some degree of genetically detectible migrants from CA found outside of that cave (inside the El Azufre River, EA; [41, 42]). This may be due, in part, to the release of Barbasco, a fish toxicant containing rotenone, during an annual fertility ceremony (La Pesca) of the indigenous Zoque people inside the CA. As Barbasco does not lead to 100% mortality rates, it was suggested that gene flow between the two habitat types may actually be mediated by a certain degree of downstream drift of sedated individuals [58].

Strong reproductive isolation among populations from ecologically divergent habitat types appears to be the result of a combination of natural selection (i.e., direct effects of toxicity, darkness, and predation) and sexual selection through female choice [5, 8, 26, 59], both of which are acting against immigrant individuals. Specifically, H₂S was shown to be a strong selective force in the aforementioned systems as revealed by reciprocal translocation experiments between nonsulphidic and sulphidic surface habitats [5]. Fish from nonsulphidic habitats had low survival in sulphidic habitats, whereas fish from sulphidic habitats performed poorly under nonsulphidic conditions. Those differences are underlined by tests on H₂S tolerances as fish from sulphidic habitats exhibited consistently higher tolerances than fish from nonsulphidic habitats [6]. The high mortalities of fish in translocations from sulphidic into nonsulphidic environments were hypothesized to be caused by oxidative stress, as oxygen is inherently toxic due to its biotransformation into reactive oxygen species, and organisms have evolved biochemical pathways with antioxidant activity (e.g., superoxide dismutase, catalase, and glutathione systems [60]). During hypoxia, the expression of antioxidant enzymes is often downregulated [61, 62], such that subsequent exposure to normoxic conditions causes substantial oxidative stress with profound fitness consequences [61, 63]. Oxidative stress, possibly in combination with the often poor body condition and energy limitation of fish from sulphidic habitats [25, 27–29, 52], may explain the high mortality seen in migrants from sulphidic to sulphide-free environments.

Contrary to translocations between sulphidic and nonsulphidic habitats, a transfer of fish between sulphidic cave and surface habitats had no effect on survival in either direction. This is not unexpected, as presence or absence of light is unlikely to affect survival within only 24 h. Nevertheless, common garden experiments found that while surface females fail to reproduce in darkness, cave females reared in light are not affected [59]. This is congruent with the aforementioned pattern of unidirectional gene flow from the inside of the caves towards the outside in the Cueva del Azufre system [6, 41, 42]. A further natural selection factor against immigrants was uncovered through similar translocation experiments (outside versus inside cave) that involved the presence of a predator (a giant water bug of the genus Belostoma) as heteropterans were more likely to attack cavefish in light but surface fish within the cave [64].

Beside environmental factors acting more or less directly on the viability of migrants in foreign habitats also sexual selection was found to constitute a reproductive isolation barrier. Thus far, only effects of intersexual selection (female mate choice) were assessed. For example, females from the Cueva del Azufre system, including normal surface habitats, sulphidic surface habitats (EA), and the Cueva del Azufre cave (front chambers of CA, which still receive some dim light), discriminate against males from foreign habitats and preferentially associate with males from their own habitat type [5]. Similarly, in the sulphur molly system female P. mexicana show strong assortative mating under nonsulphidic conditions, that is, associated less with males of the sulphur-endemic P. sulphuraria [8]. Immigrant males from ecologically divergent habitats are consequently at a disadvantage by sexual selection (see also [65, 66]).

In the present paper, we addressed another aspect of sexual selection, namely, intrasexual selection, and asked whether divergent evolution of male aggressive behaviour (i.e., competitive abilities) could play another role in facilitating reproductive isolation among diverging populations by selecting against (maladapted) migrant males. Cave mollies from CA are well known for their reduced aggressiveness [16, 17], and this reduction appears to increase gradually from the entrance to the innermost parts of the cave [67, 68]. When analysing aggressive behaviour with light of various intensities fights occurred at first at 5 lux [69]. When hybrids and backcrosses between cave and epigean fish were tested [68], the frequency of distribution patterns for aggressive fin erection and S-position revealed a genetically based reduction of the aggressive behaviour within CA fish. The F₁ generation had an intermediate value for the average, and the variability was practically halved in comparison to the epigean forms. It was concluded that the reduction for aggression is based and controlled by a polygenic genetic system. Furthermore, it was suggested that costly aggressive behaviours lack stabilising selection in darkness where visual perception of an opponent is prohibited; accordingly, reduced aggression was interpreted as a consequence of cave adaptation, that is, evolution under perpetual darkness [17]. Due to the young age of the CA cave molly this reduction process is thought to be still ongoing, eventually leading to the complete reduction of aggressive behaviour in this cave-dwelling population. Despite the extensive work on male aggression in fish from the CA, nothing is known about male aggressive behaviour of P. mexicana from the newly discovered sulphide-free Cueva Luna Azufre (LA) which is thought to have been colonized even more recently than the neighbouring CA cave [45]. Moreover, little is known about whether or not presence of toxic H₂S also plays a role for the reduction of aggressive behaviour and, if this was the case, whether evidence for convergent evolution in other drainages containing H₂S can be uncovered. Our hypothesis that not only darkness in caves, but also H₂S might affect the evolution of aggressive behaviour is based on the following considerations. Fish from H₂S-containing waters were found to have lower body conditions and fat stores [27–30, 52], most probably due to altered time budgets because of the amount of time being spent in ASR [26] and the physiological cost of H₂S detoxification [22]. These factors have been hypothesized to account for the observed heritable reduction of male sexual activity and sexual harassment of females found in all extremophile populations [50, 70].

In the present study, we asked the following specific questions.(a)What are the independent and interactive effects of H₂S and darkness on the evolution of aggressive behaviour in the Cueva del Azufre system? Do both stressors (H₂S and darkness) select for reduced aggression? We observed the outcome of dyadic aggressive interactions in male pairs from all divergent populations in the Cueva del Azufre system. For P. mexicana ecotypes from this system broad-sense heritability of population differences in the tendency to respond aggressively could be estimated by investigating laboratory- (i.e., common garden-) reared fish.(b)Is there evidence for convergent evolution (i.e., reduction) of aggressive behaviour in another system with high and sustained H₂S, namely, P. sulphuraria inhabiting the Baños del Azufre? While fish from the Cueva del Azufre system are easy to maintain and readily reproduce in the laboratory under nonsulphidic light conditions [17, 59], none of our attempts to breed P. sulphuraria have been successful so far, as fish would typically die within some weeks upon transfer to the lab. Therefore, for the comparisons among ecotypes in this system we had to rely on wild-caught fish and conducted our experiments on site in Southern Mexico.(c)Cave-adapted blind characids (Astyanax mexicanus) show reduced aggression [71], but were found to increase aggressiveness and to defend small feeding territories when starved [72]. Based on these findings, we asked if P. mexicana from CA and EA (i.e., populations showing reduced aggression) would also become more aggressive when starved and thus compared aggressive behaviour of male dyads that had undergone different feeding treatments (high diet versus one week starvation).(d)Does divergent evolution of aggressive behaviour in extremophile mollies translate into males being inferior in competition with more aggressive males from populations evolving under benign conditions? We simulated a potential migration scenario where the least aggressive CA males were paired with males from a nonsulphidic, normoxic surface stream and investigated their aggressive interactions as well.

2. Materials and Methods

2.1. Study System

The Atlantic molly, P. mexicana, is widespread in freshwater surface habitats along the Atlantic versant of Central America [73]. For our experiments we used both wild-caught fish (experiment 3) and lab-reared descendents of wild-caught fish (all other experiments). Laboratory-reared P. mexicana originated from the Río Oxolotán (Ox), a river with mostly clear water in the vicinity of the caves [6, 30], and from the brackish coastal waters near Tampico (Tam; Tamaulipas, eastern Mexico). As representatives from extreme habitats we used descendents from the sulfidic El Azufre (EA), a creek flowing out of the Cueva del Azufre [6, 30]. We furthermore used fish from three distinct cave chambers of the sulphidic Cueva del Azufre (chamber II (CA-II), chamber V (CA-V), and chamber X (CA-X); after [44]) and males from the newly discovered nonsulphidic Luna Azufre cave (LA, [45]). Wild-caught fish for experiment  3 were P. mexicana from the nonsulphidic Río Ixtapangajoya (IX, [74]) and from chamber II of the Cueva del Azufre (CA-II), as well as male P. sulphuraria (PS) from the Baños del Azufre [40]. GPS coordinates for all sampling localities are given in Table 1.

Table 1: Overview of populations used in this study. Given are relevant ecological habitat parameters [light absent (−) or present (+); H₂S absent (−) or present (+)], origin of test individuals [lab-reared (lr) or wild-caught (wc)], as well as coordinates of the sampling sites.

Light
Sulphide
Origin
Latitude
Longitude

Tampico (Tam)
+
−
lr
22.29632
−97.90022

Río Oxolotán (Ox)
+
−
lr
17.44444
−92.76293

El Azufre (EA)
+
+
lr
17.44225
−92.77447

Cueva del Azufre II (CA-II)
−
+
lr
17.44225
−92.77447

Cueva del Azufre V (CA-V)
−
+
lr
17.44225
−92.77447

Cueva del Azufre X (CA-X)
−
+
lr
17.44225
−92.77447

Cueva Luna Azufre (LA)
−
−
lr
17.44225
−92.77447

Río Ixtapangajoya (IX)
+
−
wc
17.49450
−92.99763

Poecilia sulphuraria (PS)
+
+
wc
17.55225
−92.99859

Cueva del Azufre II (CA-II)
−
+
wc
17.44225
−92.77447

2.2. Test Fish and Their Maintenance

Laboratory stocks were maintained in large, randomly outbred single-species tanks at the Department of Ecology and Evolution of the University of Frankfurt or at the Department of Zoology at the University of Oklahoma in Norman. At both facilities, fish were reared as mixed-sex stocks in 200-L (Frankfurt: Tam, Ox, LA) or 1,000-L tanks (Norman: EA, CA-II, CA-V, CA-X) at 25–27°C under a 12 : 12 hrs light : dark cycle (Frankfurt) or ambient light conditions (Norman) and were fed ad libitum at least once daily with commercial flake food. All lab-reared fish were kept under normoxic conditions without H₂S, and test fish were descendants of wild-caught fish of the 2nd to 4th laboratory generation.

In experiment  3 we used wild-caught fish, because P. sulphuraria could not be maintained under laboratory conditions for more than some weeks, most probably due to their high degree of adaptation to H₂S-containing water [6]. Upon capture, fish were transferred into closed and aerated (38 L, 43 × 31 × 32 cm) black Sterilite containers, and we gave them 24 h to acclimate before testing them in a field laboratory as described below.

2.3. Behavioural Tests

2.3.1. General Testing Procedure

We determined male aggressive behaviours during dyadic encounters by analysing contests staged between pairs of males in a small test tank measuring 30 × 20 ×20 cm. To avoid any confounding effects of previously established dominance and/or familiarity (see [75, 76]), males of each dyad were taken from different stock tanks. We separated both males by an opaque filter sponge while all sides of the test tank were taped with grey paper to minimize disturbances from the outside. The bottom of the tank was filled with black gravel, and water was kept at 27–29°C and aerated. All experiments were performed with normoxic, nonsulphidic water. Males could habituate to the test tank overnight, and fight observations took place the next day between 09:00 and 13:00. As even size-matched males differed slightly in their fin and general body colouration and were thus easily distinguishable, we noted individual characteristics of both males prior to the fights. At the start of the experiment, the partition separating both males was lifted, and we noted male-male interactions for a maximum of 10 minutes, starting with the first male-male interaction. We focused on three aggressive behaviours that occur frequently in Poecilia spp. (after [16]). (1) S-position: this threatening behaviour usually initiates a fight. Males swim in a parallel or antiparallel position and bend their bodies in an S-shaped manner while all unpaired fins are erected. (2) Tail beats: S-positons are often followed or superimposed by tail beats, which are fast moves of head and tail in opposing direction that either touch the opponent’s body or send shock waves towards the opponent. (3) Bites: we defined all incidences of ramming and mouth attacks into the direction of the opponent as bites, because these behaviours occur too fast and are too similar to be distinguishable by the human eye.

We also recorded the duration of the fights until dominance was established. Contest outcome was indicated by behavioural differences between the competitors. Folded fins, head-down posture, and a position at the periphery of the tank typically characterized the loser of the contest [77]. Winners, on the other hand, constantly chased and further attacked the losers with spread fins while occasionally performing S-positions. We, therefore, separated both males immediately after dominance was established to avoid serious injuries. If no dominance was established after 10 minutes of observation we terminated the fight and scored fights as “no clear winner”; those trials were discarded from the analysis of fighting durations, while fight durations were scored as “0” when no aggressive behaviour occurred at all. After a contest, body size of all males was measured as standard length (SL) to the nearest millimetre by laying the fish flat on plastic foil-covered millimetre paper, and we transferred the males back into their respective stock tanks. Despite the loss of single scales, no severe injuries and no mortality related to the experiments were observed.

2.3.2. Experiment  1: When Do Fights Escalate?

The assessment of an opponent’s resource holding potential (RHP; see [78]) is crucially connected to the opponents’ body size difference in poeciliid fishes [77], and at least in swordtail fish (Xiphophorus spp.) fighting intensity (determined as numbers of bites per minute) correlates negatively with the opponents’ size difference, but simultaneously was also found to vary greatly when size differences were small [77]. We, therefore, first examined the correlation between fight intensity (determined as numbers of bites per minute) and the opponents’ absolute body size (measured as standard length) difference using Spearman’s rank order test. We tested 17 male dyads of P. mexicana (Tam), while relative size differences within each pair ranged from 0% to 47% (mean pair size:  mm). Fight intensity was plotted against opponents’ absolute body size difference, and a logarithmic model was used to generate a reference line.

Despite the predicted large degree of variability in fight intensity (see [77]), escalating fights in swordtails (defined as both males biting each other) are more often found when body size differences are small [76, 77]. As the aim of our study was to compare maximum aggressiveness of escalating fights in different P. mexicana populations we furthermore tried to estimate the maximum relative opponents’ size difference (determined as fraction of standard length the larger male exceeded the smaller male) up to which escalating fights can be observed. To do so, a score expressing how equally both males dedicated their aggressions towards each other in a dyad was calculated as a measure of escalation as: 1 – Abs((sum aggressive behaviours male one − sum aggressive behaviours male two)/sum of aggressive behaviours of male one and male two).

Escalation scores could range between 0 and 1, with values around 0 indicating that only one male showed aggressive behaviour (uneven, no escalated fight) and values around 1 indicating that both males dedicated similar amounts of aggressive behaviours towards each other (even, escalated fight). Scores were plotted against arcsine- (squareroot-) transformed relative body-size differences, and a logistic 4-parameter curve estimation (upper constraint set to 1, lower constraint to 0) was used to determine the Evenness₅₀-score (body size difference at which the score value is 0.5). When opponents’ body size differences exceeded the body size difference at the Evenness₅₀-score we assumed fights to be less escalated.

For statistical reasons, scores and body-size differences equal 0 were substituted by 0.001, as logistic models require positive nonzero data.

2.3.3. Experiment  2: Evolution of Male Aggressiveness in Response to Environmental Stressors

It was our intention to disentangle the relative effects of sulphur and darkness on the evolution of aggressive behaviour. In a first step we tested whether populations from sites with the same combination of ecological stressors would show comparable levels of aggressive behaviour and, thus, compared the two populations from nonsulphidic surface sites (Tam and Ox) as well as the three CA populations (CA-II, CA-V, CA-X) using similar MANCOVA and ANCOVA models as outlined below. The MANCOVA models with numbers of aggressive behaviours per fight as dependent variables neither detected a significant population difference between both nonsulphidic surface populations () nor between the three CA populations (). When comparing fight durations using ANCOVA models we did not find population differences as well (surface: ; CA: ).

Based on these results we analysed numbers of aggressive behaviours per fight in seven populations of P. mexicana from different ecological backgrounds in our main analysis (see Table 1 for a detailed description of collection sites): A MANCOVA model with “number of S-positions,” “number of tail-beats,” and “number of bites” as dependent variables and “sulphur” (present/absent) as well as “light” (present/absent) as fixed factors was employed, and “mean pair size” as well as “body size difference” (arcsine (square root) transformed) were included as covariates. We initially included all levels of interaction terms between both main factors and both covariates, but removed interactions from the final model as none of them were significant (not shown). Prior to all analyses, all dependent variables were log transformed and afterwards checked for normal distribution by Kolmogorov-Smirnov tests.

Fighting durations were analysed in a separate ANCOVA model with “sulphur” (present/absent) as well as “light” (present/absent) as fixed factors and “mean pair size” as well as “body size difference” as covariates. No interaction term was significant (not shown), and thus interactions were excluded from the final model.

2.3.4. Experiment  3: Aggressiveness in the Sulphur Endemic P. sulphuraria

In this experiment we compared male aggressive interactions among wild-caught individuals of the sulphur endemic P. sulphuraria (PS; ) and two populations of P. mexicana, one from a freshwater surface habitat, the Río Ixtapangajoya (IX; ), the other one from the sulphidic Cueva del Azufre (CA-II; male dyads), in a field laboratory. We used small mice cages (23 × 15 × 16.5 cm) instead of our standard test tanks and separated males overnight by opaque plastic sheets. As described for experiment  2, we analysed aggressive behaviours among those three populations through MANCOVA with “mean pair size” as well as “body size difference” as covariates and fight durations in an ANCOVA with “mean pair size” as well as “body size difference” as covariates. In both analysis, interaction terms between the main factor “population” and the covariates were initially included, but removed from the final model as neither had a significant effect.

We used Fisher’s LSD tests for pairwise post hoc comparison of overall levels of aggression (mean sum of all aggressive interactions per fight) as well as fight durations among populations. In addition, we also compared all three different kinds of aggressive behaviours separately by use of one-way ANOVA and applied Fisher’s LSD tests to identify the source of variation when a significant population effect was detected.

2.3.5. Experiment  4: Aggressiveness and Food Limitation

As food limitation is predicted to influence the occurrence of aggressive behaviours in fish [79], we compared the intensity of male fights under normal food supply (daily ad libitum feeding regime) with fights of males that were subject to a 1-week starvation period. To do so, we separated males from the CA-X and EA populations for 6 days in 50-L tanks and deprived them of food. After this period males that had not been fighting against each other (i.e., stemmed from different tanks) were transferred into our standard test tanks, and fights were observed on the following morning (hence, males starved for 7 days altogether). We analysed numbers of aggressive behaviours per fight (including data for nonstarved males from experiment  2) in a MANCOVA model with “population” (2 levels) and “treatment” (non-starved/starved) while including “mean pair size” as well as “body size difference” as covariates. Analogously, fight durations were analysed with the same factors and covariates in an ANCOVA model. In both analysis, interaction terms of the main factors “population” and “treatment” and the covariates were initially included but removed from the final model as neither had a significant effect.

2.3.6. Experiment  5: Male Aggression as Reproductive Isolation Barrier

Reproductive isolation is crucial for speciation processes, and intrasexual competition may provide one possible mechanism to exclude immigrants from reproducing in foreign habitats. To test this idea, we staged contests between size-matched males from CA (CA-X; mean SL =  mm) versus males from a sulphide-free surface habitats (Ox population; ; paired t-test on size differences: ; ). Chi² tests were employed to compare numbers of fights won by males from either population, and numbers of aggressive behaviours shown by the two ecotypes within each male dyad were analysed using paired t-tests. We further recorded and compared all sexual behaviours (nipping and copulation attempts, so-called thrusting; see [16] for a description) between both male types, as cave mollies may answer aggressions by sexually motivated behaviours [68].

3. Results

3.1. Body Size Difference between Opponents and Male Aggressive Behaviour

In our first experiment we quantified fight intensities and durations in staged contests of P. mexicana males from the Tampico population. Body size differences between both males within a dyad varied between 0 and 8 mm. Fight intensity (measured as bites per minute) was negatively correlated with the opponents’ body size difference (Spearman rank order test; ; Figure 1(a)) meaning that males fought most intensely when both opponents were closely size matched. The body size difference below which fights escalated (i.e., below which both males displayed equal numbers of aggressive behaviours; “fight evenness”) was determined as 7.7%, with the 95% confidence interval ranging between 5.1% and 12.2% (Figure 1(b); Logistic model: ). Based on these results, we made an attempt to use closely size-matched male pairs in all subsequent experiments [mean (±SD) size difference = %] and included arcsine- (square root-) transformed relative body size difference of each dyad as a covariate in all further analyses.

(a)

(b)

Figure 1: (a) Fight intensity and (b) “fight evenness” (see main text) in relation to the two opponents’ body size difference (as standard length (SL) difference). Shown are regression lines representing the best-fit (a: logarithmic model; b: logistic model) and 95% confidence intervals ( fights).

3.2. Evolution of Male Aggressiveness in the Cueva del Azufre System

MANCOVA revealed a significant effect of the factor “light” (Table 2), indicating that cave-dwelling populations displayed significantly fewer aggressive behaviours than surface fish (Figure 2(a)). The significant “light × sulphide” interaction (Table 2) further indicates that this reduction in aggressiveness is aggravated in populations evolving under both extreme conditions, while “sulphur” per se did not lead to a significant reduction in aggressive behaviours (Table 2; Figure 2(a)). Also both covariates (“mean opponent body size” and “body size difference”) had a significant influence in our model (Table 2), and post hoc Spearman rank-order tests revealed that “mean opponent body size” was positively correlated with numbers of S-positions , tail beats , and bites per fight , suggesting that fights of larger males were more intense than those of smaller ones. In contrast, the body size difference between both opponents was negatively correlated with the number of S-positions and tail beats , but not fights , indicating that the larger the opponents’ body size difference was, the less intense fights became.

Table 2: Results from (a) MANCOVA and (b) ANCOVA models analysing attributes of dyadic male aggressive interactions in Experiment  2 (lab-reared males). -ratios were approximated using Wilk’s . Partial variance was estimated using Wilk’s partial . Significant effects are in bold typeface.

df

Partial variance explained [%]

(a) MANCOVA (number of aggressive behaviours)

    Light (absent/present)
3,62
8.97
<0.001
0.30

    Sulphide (absent/present)
3,62
2.45
0.072
0.10

    Light × sulphide
3,62
3.77
0.015
0.15

    Male body size difference
3,62
3.35
0.025
0.14

    Mean pair body size
3,62
2.81
0.044
0.12

(b) ANCOVA (fight duration)

    Light (absent/present)
1
8.44
0.005
0.12

    Sulphide (absent/present)
1
7.07
0.010
0.10

    Light × sulphide
1
2.48
0.120
0.04

    Male body size difference
1
0.25
0.622
<0.01

    Mean pair body size
1
0.71
0.403
0.01

    Error
64

(a)

(b)

Figure 2: Means (±SE) of (a) numbers of aggressive interactions and (b) fight duration in seven populations of P. mexicana. From left to right: populations from nonsulphidic surface sites [Tampico, Tam (), and Río Oxolotán, Ox ()], the sulphidic creek in the Cueva del Azufre system [El Azufre, EA ()], three cave chambers of the sulphidic Cueva del Azufre [CA-II (), CA-V (), CA-X ()], and the H₂S-free cave [Luna Azufre, LA ()].

When comparing the duration of fights we found both main factors (“light” and “sulphide”) to have significant effects (Table 2). This and the nonsignificant interaction term of “light × sulphide” suggest that both the absence of light and the presence of H₂S lead to similar reductions in fighting time (Figure 2(b)).

3.3. Aggressiveness in the Sulphur Endemic P. sulphuraria

When comparing numbers of aggressive behaviours in fights of wild-caught males from two P. mexicana populations (IX and CA-II) and P. sulphuraria males by use of MANCOVA we found a significant effect of the factor “population/species” , and post hoc pairwise comparisons (Fisher’s LSD) showed levels of aggressive behaviours of surface P. mexicana and P. sulphuraria to differ significantly from those seen in P. mexicana males from CA-II (Figure 3(a)). None of the covariates had a significant effect (male body size difference: ; mean pair body size: ). One-way ANOVAs confirmed significant differences between populations in all three aggressive behaviours (S-position: ; tail beats: ; bites: ). Post hoc tests revealed that fights between P. sulphuraria males were characterized by significantly more S-positions compared to fights of IX males (), and fights of both surface populations/species displayed significantly more tail beats compared to fights of CA-II males (IX versus CA-II: ; PS versus CA-II: ). All three populations differed significantly in numbers of bites per fight (IX versus CA-II: ; IX versus PS: ; PS versus CA-II: ).

(a)

(b)

Figure 3: Means (±SE) of (a) numbers of aggressive behaviours shown by males during dyadic fights and (b) duration of fights in wild-caught males from two P. mexicana populations [the sulphide-free Río Ixtapangajoya, IX () and cave chamber II of the sulphidic Cueva del Azufre, CA-II ()], as well as the sulphur-endemic P. sulphuraria (PS) found at the Baños del Azufre ().

When analysing the durations of fights, our ANCOVA model detected a significant effect of the factor “population/species” , and pairwise comparisons showed that both surface forms (IX and PS) fought significantly longer than P. mexicana males from CA-II (Figure 3(b)). Again, both covariates were not significant (male body size difference: ; mean pair body size: ).

3.4. Aggressiveness and Food Limitation

When comparing numbers of aggressive behaviours in fights of P. mexicana from CA-X and EA under normal food supply (data from Experiment  2) and after one week of starvation in a MANCOVA we found a significant effect of the factor “food treatment” while the factor “population” as well as the interaction term “treatment × population” was not significant . This indicates that both populations reduced their aggressive behaviour in a similar fashion when food was scarce (Figure 4(a)).

(a)

(b)

Figure 4: Means (±SE) of (a) numbers of aggressive behaviours and (b) the duration of male fights in males from two extremophile populations in the Cueva del Azufre system, one from the sulphidic surface stream [El Azufre, EA ()] and one from chamber X of the sulphidic Cueva del Azufre [CA-X ()]. Prior to the tests, males were either fed on a normal diet (left bars) or starved for one week (s: right bars).

Like in experiment  2, we found the covariate “mean opponent size” to have a significant effect in the MANCOVA , and post hoc Spearman rank-order tests revealed a significant positive correlation between “mean opponent size” and numbers of S-positions and tail beats , but not bites . The covariate “body size difference” had no significant effect .

Another ANCOVA model analysing fighting durations revealed a significant effect of the factor “food treatment” while the factor “population” and the interaction term “treatment × population” were not significant . Furthermore, both covariates had no significant effects (“body size difference”: ; “mean opponent size”: ). Starvation in general reduced the duration of fights (Figure 4(b)).

3.5. Fights between Different Locally Adapted Males

In all 13 staged contests, Río Oxolotán (Ox) males established dominance over the CA-X males after a mean fight duration of 119 ± 19 s. Ox males directed significantly more aggressive behaviours towards CA-X males (S-position: ; tail-beats: ; bites: ; Figure 5) while cave molly males directed more sexually motivated behaviours towards Ox males during the fights (nipping: ; thrusting: ; Figure 5).

Figure 5: Mean (±SE) numbers of aggressive and sexual behaviours shown by males in cross-population dyadic fights () between Río Oxolotán (Ox) males and males from chamber X of the sulphidic Cueva del Azufre (CA-X). All comparisons were statistically different between populations (see main text).

4. Discussion

An increasing body of literature documents adaptation’s potential to drive genetic differentiation and ultimately speciation (e.g., [80–83]), a phenomenon that has recently been termed “isolation by adaptation” [84]. Of particular interest in the study of ecological speciation are the proximate mechanisms leading to and maintaining genetic differentiation among populations [2, 84]. During ecological speciation, prezygotic isolation may arise when immigrants from foreign, ecologically divergent habitats are selected against [85, 86]. This may occur through natural selection, if immigrants (or hybrids) have reduced viability (extrinsic reproductive isolation; e.g., [5, 6, 87]), or by sexual selection, if maladapted individuals are discriminated against during mate choice (e.g., [88]).

In the present study we examined whether—in addition to mate choice (i.e., intersexual selection)—intrasexual selection through male-male competition could also play a role in promoting prezygotic isolation. Atlantic molly males in clear-water habitats usually establish dominance hierarchies, and dominant (typically the largest) males monopolize several females which they aggressively defend against rivals [16]. This view is supported by our present findings in that fighting intensity was positively correlated with the average body size of male dyads; in other words, larger males fought more intensely, probably driven by the prospect of monopolizing females. Smaller males, by contrast, rely on a sneak-like mating tactic [50, 70], but such “alternative” mating tactics are lost in extremophile P. mexicana [50], likely owing to very similar counterselection in energy-limited habitats that, as we will discuss, may have played a role for the evolutionary reduction of aggressive behaviour (see below).

We found fight intensity to be reduced in various extremophile P. mexicana populations, and perpetual darkness in caves was the best predictor for the evolutionary reduction of aggressiveness, especially when it was combined with presence of H₂S, as seen in the CA cave. As lab-reared fish were used for this part of our study, the observed differences seem to be largely evolved (genetic) differences among ecotypes. When we considered fight durations, also a significant main effect of the factor “H₂S” was observed; fish from sulphidic habitats engage in shorter fights as an evolutionary response to the toxicant. Finally, we demonstrate that reduced aggression directly translates into males being inferior in contests, as evidenced by the fact that Ox males always won when paired with cave molly (CA) males; CA males even responded sexually to aggressive attacks, an obviously maladaptive behaviour (see also [68]). We argue that in a system where dominance hierarchies play a vital role, reduced aggressiveness translates directly into male inferiority in mate competition upon encounter of different behavioural phenotypes. Specifically, we argue that migrant males stemming from an ecological background that has selected for reduced aggressiveness may be selected against (i.e., have low reproductive fitness) in a divergent (i.e., benign) habitat type. Together with the action of natural selection against (maladapted) migrants via H₂S-toxicity, darkness and predation [5, 6, 8, 26, 59], as well as female mate discrimination against alien male phenotypes [5, 8], divergent evolution of aggressive behaviour may thus play an important role for the maintenance of genetic differentiation in this system—at least at the interface between extreme and benign (nonsulphidic surface) habitats and, hence, could represent another mechanistic link explaining the surprising small-scale genetic structuring in the CA system [41, 42]. It remains to be determined in future studies whether the comparatively small differences in the intensity of aggressive behaviour seen in males from some habitats that are directly adjoining in the Cueva del Azufre system (e.g., CA versus EA, LA versus EA) lead to an equally clearcut picture, that is, if also in those cases it is always the more aggressive males that win a combat. For practical reasons, our present study focussed on aggressive interactions between the most extreme behavioural phenotypes: males from the most aggressive (Ox) population and the least aggressive CA-X males.

Parzefall first described reduced aggression in the CA cave population of P. mexicana (from the rearmost chamber XIII; CA-XIII) [16] and interpreted his findings as an adaptation to perpetual lightless conditions [17], as most aggressive behaviours depend on visual perception of cues from opponents, which may be more difficult to perceive in darkness. Even though theory predicts a reduction of intraspecific aggression in troglobites [89, 90], some cave dwellers may even have evolved entire novel sets of aggressive behaviours while responding to nonvisual signals. For example, aggressive behaviour is well developed in the blind catfish Uegitglanis zammaranoi [91] and the blind cave salamander Proteus anguinus [92]. Furthermore, the discovery of a highly aggressive cave-dwelling Astyanax fasciatus population [93] implies that a reduction of aggressive behaviour is not an inevitable evolutionary response to the cave environment. Those authors suggested that explanations other than simply the inability to perceive visual cues triggering aggressive behaviour should be explored in order to explain the evolutionary reduction of aggressiveness in many other Astyanax cave fish populations [94, 95].

As we have argued above, previous studies have demonstrated that P. mexicana inside the two caves (CA and LA), as well as the toxic surface habitat (EA), appear to be energy limited, as evidenced by their lower body conditions and reduced fat stores [27–30, 52]. In the nontoxic LA cave, this is probably due to low resource availability, which is typical for most caves (reviewed in [37, 96]). In contrast, CA and EA are energy-rich habitats due to high chemoautotrophic primary productivity [97, 98]; however, P. mexicana spend the majority of their time at the water surface engaged in ASR [26] and probably pay a high physiological cost in order to run ATP-expensive H₂S detoxification [22]. Not surprisingly, a recent study therefore found cave mollies from CA to have higher metabolic rates compared to surface mollies even after several generations in the laboratory [99]. Altogether, this suggests that reduced aggression is most likely an evolutionary response to continued energy limitation in the Cueva del Azufre system. This hypothesis is corroborated by the results from this study, in which we found that P. mexicana from EA and CA plastically reduce their aggression even further after being starved for one week.

It is further interesting to note that not only did overall levels of aggression diverge between extremophile and nonextremophile poeciliids, but also the relative contribution of specific aggressive behaviours to the aggressive repertoire of these species. The potentially most harmful aggressive behaviour (i.e., bites and rammings) was strongly reduced in all extremophile poeciliids, while the least harmful behaviour (i.e., S-positions) was actually increased in extremophiles. As the energetic costs of threat displays were found to be low relative to the costs of escalated fighting in an African cichlid species (Tilapia zillii) [100] we argue that again this phenomenon is a response to the energy limitation experienced in extreme habitats. Moreover, P. mexicana males appear to have higher energy expenditure than females [25, 101] and, therefore, exhibit higher mortality rates under stressful conditions and perform more ASR than females [25]. Assuming that male poeciliids in H₂S-toxic habitats live near the edge of survivability [25], any injuries obtained during fights with other males (especially during biting or ramming) could indeed lead to life-threatening infections and, ultimately, premature death—a hypothesis that is further supported by a recent study reporting on higher individual parasitization rates of P. mexicana in the CA and EA compared to Ox [102].

In stark contrast to the findings from the Cueva del Azufre system, our experiments using wild-caught males from another system with high and sustained H₂S, namely, P. sulphuraria inhabiting the Baños del Azufre, found no significant difference among P. mexicana from nonsulphidic sites and the “sulphur-endemic” P. sulphuraria. So, why did extremophile males from the Cueva del Azufre system show strongly reduced aggressiveness, but P. sulphuraria did not? We propose three mutuallynot exclusive hypotheses. First, our analysis of fish in the Cueva del Azufre system found the relative contribution of the factor “sulphide” to the evolutionary reduction of aggressive behaviour to be much lower than that of the factor “light” (see partial variance explained in Table 2), so these fish may just not experience the same selective pressure to reduce aggressiveness. Second, P. sulphuraria are clearly well adapted to high concentrations of H₂S (i.e., being sulphur endemics) and accordingly could show some kind of “rebound effect”, indicating evolved mechanisms to better cope with the toxicant (see [20] for discussion). Some support for the latter idea was also found in life-history traits, as P. sulphuraria actually had the largest (not the smallest) fat stores in a comparison of poeciliids from several benign and sulphidic habitats [27]. Third, contrary to P. mexicana from the Cueva del Azufre system, which are the only permanent piscine residents in their respective extreme habitats [43], P. sulphuraria have to share their habitat with another sulphide-adapted species, the widemouth gambusia, Gambusia eurystoma [40]. Hence, increased aggression could also be a signal of interspecific competition for resources at the Baños del Azufre.

In conclusion, sulphuric waters are characterized by reduced resource availability but increased energy expenditure, leading to low body conditions and fat stores in H₂S-inhabiting P. mexicana. We suggest that in addition to darkness in caves also resource limitation might play a crucial role in the evolutionary reduction of male aggressive behaviour. Selection against costly behaviours (such as aggression) might lead sulphur-adapted mollies to trade-off aggressive behaviour to compensate for the negative effects of H₂S—similar to the proposed processes resulting in the observed heritable reduction of male sexual activity and harassment of females found in all extremophile populations (e.g., [50]), as well as patterns of life history divergence [28, 29]. On the other hand, the phylogenetically old “sulphur specialist” P. sulphuraria, which we did not find to show reduced aggression, might cope better with H₂S and thus can afford to express costly aggressive behaviour.

Acknowledgments

The authors would like to thank Holger Geupel for help with animal care and N. Herrmann, C. T. Jung, S. Klaus, B. Kunkel, M. Schulte, S. Stadler, and M. Tobler for assistance in the field. Financial support came from the Deutsche Forschungsgemeinschaft (DFG; PL 470/1–3) as well as from the Herrmann-Willkomm-Stiftung (to MP, SS, SK).

References

J. A. Endler, Natural Selection in the Wild, Princeton University Press, Princeton, NJ, USA, 1986.

D. Schluter, “Ecology and the origin of species,” Trends in Ecology and Evolution, vol. 16, no. 7, pp. 372–380, 2001.

H. D. Rundle and P. Nosil, “Ecological speciation,” Ecology Letters, vol. 8, no. 3, pp. 336–352, 2005.

P. Nosil, “Reproductive isolation caused by visual predation on migrants between divergent environments,” Proceedings of the Royal Society B, vol. 271, no. 1547, pp. 1521–1528, 2004.

M. Tobler, R. Riesch, C. M. Tobler, T. Schulz-Mirbach, and M. Plath, “Natural and sexual selection against immigrants maintains differentiation among micro-allopatric populations,” Journal of Evolutionary Biology, vol. 22, no. 11, pp. 2298–2304, 2009.

M. Tobler, M. Palacios, L. J. Chapman et al., “Evolution in extreme environments: replicated phenotypic differentiation in livebearing fish inhabiting sulfidic springs,” Evolution, vol. 65, no. 8, pp. 2213–2228, 2011.

L. K. Snowberg and C. W. Benkman, “Mate choice based on a key ecological performance trait,” Journal of Evolutionary Biology, vol. 22, no. 4, pp. 762–769, 2009.

M. Plath, R. Riesch, A. Oranth et al., “Complementary effect of natural and sexual selection against immigrants maintains differentiation between locally adapted fish,” Naturwissenschaften, vol. 97, no. 8, pp. 769–774, 2010.

S. Via, A. C. Bouck, and S. Skillman, “Reproductive isolation between divergent races of pea aphids on two hosts. II. Selection against migrants and hybrids in the parental environments,” Evolution, vol. 54, no. 5, pp. 1626–1637, 2000.

M. R. Servedio, G. S. V. Doorn, M. Kopp, A. M. Frame, and P. Nosil, “Magic traits in speciation: “magic” but not rare?” Trends in Ecology and Evolution, vol. 26, no. 8, pp. 389–397, 2011.

C. Eizaguirre, T. L. Lenz, A. Traulsen, and M. Milinski, “Speciation accelerated and stabilized by pleiotropic major histocompatibility complex immunogenes,” Ecology Letters, vol. 12, no. 1, pp. 5–12, 2009.

P. G. D. Feulner, M. Plath, J. Engelmann, F. Kirschbaum, and R. Tiedemann, “Electrifying love: electric fish use species-specific discharge for mate recognition,” Biology Letters, vol. 5, no. 2, pp. 225–228, 2009.

P. G. D. Feulner, M. Plath, J. Engelmann, F. Kirschbaum, and R. Tiedemann, “Magic trait electric organ discharge (EOD): dual function of electric signals promotes speciation in African weakly electric fish,” Communitative and Integrative Biology, vol. 2, no. 4, pp. 329–331, 2009.

S. Fukamachi, M. Kinoshita, K. Aizawa, S. Oda, A. Meyer, and H. Mitani, “Dual control by a single gene of secondary sexual characters and mating preferences in medaka,” BMC Biology, vol. 7, article 1741, p. 64, 2009.

O. Seehausen, Y. Terai, I. S. Magalhaes et al., “Speciation through sensory drive in cichlid fish,” Nature, vol. 455, no. 7213, pp. 620–626, 2008.

J. Parzefall, “Zur vergleichenden Ethologie verschiedener Mollienesia-Arten einschließlich einer Hohlenform von M. sphenops,” Behaviour, vol. 33, no. 1, pp. 1–38, 1969.

J. Parzefall, “A review of morphological and behavioural changes in the cave molly, Poecilia mexicana, from Tabasco, Mexico,” Environmental Biology of Fishes, vol. 62, no. 1–3, pp. 263–275, 2001.

M. Plath, K. U. Heubel, and I. Schlupp, “Field observations on male mating behavior in surface—and cave—dwelling Atlantic mollies (Poecilia mexicana, Poeciliidae),” Zeitschrift für Fischkunde, vol. 7, no. 2, pp. 113–119, 2005.

P. Calow, “Proximate and ultimate responses to stress in biological systems,” Biological Journal of the Linnean Society, vol. 37, no. 1-2, pp. 173–181, 1989.

L. J. Rothschild and R. L. Mancinelli, “Life in extreme environments,” Nature, vol. 409, no. 6823, pp. 1092–1101, 2001.

T. Bagarinao, “Sulfide as an environmental factor and toxicant: tolerance and adaptations in aquatic organisms,” Aquatic Toxicology, vol. 24, no. 1-2, pp. 21–62, 1992.

M. K. Grieshaber and S. Völke, “Animal adaptations for tolerance and exploitation of poisonous sulfide,” Annual Review of Physiology, vol. 60, pp. 33–53, 1998.

E. R. McMullin, D. C. Bergquist, and C. R. Fisher, “Metazoans in extreme environments: adaptations of hydrothermal vent and hydrocarbon seep fauna,” Gravitational and Space Biology Bulletin, vol. 13, no. 2, pp. 13–23, 2000.

U. Theissen and W. Martin, “Biochemical and evolutionary aspects of eukaryotes that inhabit sulfidic environments,” in Microbial Sulfur Metabolism, C. Dahl and C. G. Friedrich, Eds., pp. 36–45, Springer, Berlin, Germany, 2008.

M. Plath, M. Tobler, R. Riesch, F. J. G. de León, O. Giere, and I. Schlupp, “Survival in an extreme habitat: the roles of behaviour and energy limitation,” Naturwissenschaften, vol. 94, no. 12, pp. 991–996, 2007.

M. Tobler, R. W. Riesch, C. M. Tobler, and M. Plath, “Compensatory behaviour in response to sulphide-induced hypoxia affects time budgets, feeding efficiency, and predation risk,” Evolutionary Ecology Research, vol. 11, no. 6, pp. 935–948, 2009.

R. Riesch, M. Plath, F. J. G. de León, and I. Schlupp, “Convergent life-history shifts: toxic environments result in big babies in two clades of poeciliids,” Naturwissenschaften, vol. 97, no. 2, pp. 133–141, 2010.

R. Riesch, M. Plath, and I. Schlupp, “Toxic hydrogen sulfide and dark caves: life-history adaptations in a livebearing fish (Poecilia mexicana, Poeciliidae),” Ecology, vol. 91, no. 5, pp. 1494–1505, 2010.

R. Riesch, M. Plath, and I. Schlupp, “Toxic hydrogen sulphide and dark caves: pronounced male life-history divergence among locally adapted Poecilia mexicana (Poeciliidae),” Journal of Evolutionary Biology, vol. 24, no. 3, pp. 596–606, 2011.

M. Tobler, I. Schlupp, K. U. Heubel et al., “Life on the edge: hydrogen sulfide and the fish communities of a Mexican cave and surrounding waters,” Extremophiles, vol. 10, no. 6, pp. 577–585, 2006.

T. B. de Perera, “Spatial parameters encoded in the spatial map of the blind Mexican cave fish, Astyanax fasciatus,” Animal Behaviour, vol. 68, no. 2, pp. 291–295, 2004.

M. Yoshizawa, Š. Gorički, D. Soares, and W. R. Jeffery, “Evolution of a behavioral shift mediated by superficial neuromasts helps cavefish find food in darkness,” Current Biology, vol. 20, no. 18, pp. 1631–1636, 2010.

M. Plath, M. Rohde, T. Schröder, A. Taebel-Hellwig, and I. Schlupp, “Female mating preferences in blind cave tetras Astyanax fasciatus (Characidae, Teleostei),” Behaviour, vol. 143, no. 1, pp. 15–32, 2006.

M. Tobler, I. Schlupp, and M. Plath, “Does divergence in female mate choice affect male size distributions in two cave fish populations?” Biology Letters, vol. 4, no. 5, pp. 452–454, 2008.

D. C. Culver and T. Pipan, The Biology of Caves and Other Subterranean Habitats, Oxford University Press, New York, NY, USA, 2008.

D. C. Culver, Cave Life—Evolution and Ecology, Harvard University Press, Cambridge, Mass, USA, 1982.

K. Hüppop, “How do cave animals cope with the food scarcity in caves?” in Ecosystems of the World 30: Subterranean Ecosystems, H. Wilkens, D. C. Culver, and W. F. Humphries, Eds., pp. 159–188, Elsevier Science, Amsterdam, The Netherlands, 2000.

W. R. Jeffery, “Cavefish as a model system in evolutionary developmental biology,” Developmental Biology, vol. 231, no. 1, pp. 1–12, 2001.

W. R. Jeffery, A. G. Strickler, and Y. Yamamoto, “To see or not to see: evolution of eye degeneration in Mexican blind cavefish,” Integrative and Comparative Biology, vol. 43, no. 4, pp. 531–541, 2003.

M. Tobler, R. Riesch, F. J. G. de León, I. Schlupp, and M. Plath, “Two endemic and endangered fishes, Poecilia sulphuraria (Alvarez, 1948) and Gambusia eurystoma Miller, 1975 (Poeciliidae, Teleostei) as only survivors in a small sulphidic habitat,” Journal of Fish Biology, vol. 72, no. 3, pp. 523–533, 2008.

M. Plath, J. S. Hauswaldt, K. Moll et al., “Local adaptation and pronounced genetic differentiation in an extremophile fish, Poecilia mexicana, inhabiting a Mexican cave with toxic hydrogen sulphide,” Molecular Ecology, vol. 16, no. 5, pp. 967–976, 2007.

M. Plath, B. Hermann, C. Schröder, et al., “Catastrophic flood does not lead to loss of small—scale genetic differentiation among locally adapted fish populations,” BMC Evolutionary Biology, vol. E10, no. 1, p. 256, 2010.

M. Tobler, T. J. DeWitt, I. Schlupp et al., “Toxic hydrogen sulfide and dark caves: phenotypic and genetic divergence across two abiotic environmental gradients in Poecilia mexicana,” Evolution, vol. 62, no. 10, pp. 2643–2659, 2008.

M. S. Gordon and D. E. Rosen, “A cavernicolous form of the poeciliid fish Poecilia sphenops from Tabasco, Mexico,” Copeia, vol. 1962, no. 2, pp. 360–368, 1962.

M. Tobler, R. Riesch, F. J. G. de León, I. Schlupp, and M. Plath, “A new and morphologically distinct population of cavernicolous Poecilia mexicana (Poeciliidae: Teleostei),” Environmental Biology of Fishes, vol. 82, no. 1, pp. 101–108, 2008.

J. Pisarowicz, “Return to tabasco,” Association for Mexican Cave Studies Activities Newsletter, vol. 28, no. 1, pp. 27–57, 2005.

L. Rosales-Lagarde, P. J. Boston, A. Campbell, and K. W. Stafford, “Possible structural connection between Chichón Volcano and the sulfur-rich springs of Villa Luz Cave (a. k. a. Cueva de las Sardinas), Southern Mexico,” Association for Mexican Cave Studies Bulletin, vol. 19, no. 7, pp. 177–184, 2006.

L. Rosales-Lagarde, A. Campbell, P. J. Boston, and K. W. Stafford, “Sulfur and oxygen isotopes: evidence of H₂S spring sources, southern Mexico,” Geochimica et Cosmochimica Acta, vol. 72, no. 12, p. A805, 2008.

M. Plath, J. Parzefall, K. E. Körner, and I. Schlupp, “Sexual selection in darkness? Female mating preferences in surface- and cave-dwelling Atlantic mollies, Poecilia mexicana (Poeciliidae, Teleostei),” Behavioral Ecology and Sociobiology, vol. 55, no. 6, pp. 596–601, 2004.

M. Plath, “Male mating behavior and costs of sexual harassment for females in cavernicolous and extremophile populations of Atlantic mollies (Poecilia mexicana),” Behaviour, vol. 145, no. 1, pp. 73–89, 2008.

R. Riesch, V. Duwe, N. Herrmann et al., “Variation along the shy-bold continuum in extremophile fishes (Poecilia mexicana, Poecilia sulphuraria),” Behavioral Ecology and Sociobiology, vol. 63, no. 10, pp. 1515–1526, 2009.

M. Tobler, “Divergence in trophic ecology characterizes colonization of extreme habitats,” Biological Journal of the Linnean Society, vol. 95, no. 3, pp. 517–528, 2008.

R. Riesch, M. Tobler, M. Plath, and I. Schlupp, “Offspring number in a livebearing fish (Poecilia mexicana, Poeciliidae): reduced fecundity and reduced plasticity in a population of cave mollies,” Environmental Biology of Fishes, vol. 84, no. 1, pp. 89–94, 2009.

R. Riesch, M. Plath, I. Schlupp, and E. Marsh-Matthews, “Matrotrophy in the cave molly: an unexpected provisioning strategy in an extreme environment,” Evolutionary Ecology, vol. 24, no. 4, pp. 789–801, 2010.

N. Peters and G. Peters, “Das auge einer höhlenform von P. sphenops (Poeciliidae, Pisces). Ein versuch, morphologische gesetzmäßigkeiten der degenerativen evolution genetisch zu interpretieren,” Zeitschrift für Morphologie der Tiere, vol. 62, no. 3, pp. 211–244, 1968.

N. Peters, G. Peters, J. Parzefall, and H. Wilkens, “Über degenerative und konstruktive Merkmale bei einer phylogenetisch jungen Höhlenform von Poecilia sphenops (Pisces, Poeciliidae),” Internationale Revue der Gesamten Hydrobiologie, vol. 58, no. 3, pp. 417–436, 1973.

M. Plath and M. Tobler, “Subterranean fishes of Mexico (Poecilia mexicana, Poeciliidae),” in Biology of Subterranean Fishes, T. Trajano, M. E. Bichuette, and B. G. Kapoor, Eds., pp. 281–330, Science Publishers, Enfield, NH, USA, 1st edition, 2010.

M. Tobler, Z. W. Culumber, M. Plath, K. O. Winemiller, and G. G. Rosenthal, “An indigenous religious ritual selects for resistance to a toxicant in a fish (Poecilia mexicana, Poeciliidae),” Biology Letters, vol. 7, no. 2, pp. 229–232, 2011.

R. Riesch, M. Plath, and I. Schlupp, “Speciation in caves: experimental evidence that permanent darkness promotes reproductive isolation,” Biology Letters, vol. 7, no. 6, pp. 909–912, 2011.

B. Halliwell and M. Gutteridge, Free Radicals in Biology and Medicine, Oxford University Press, Oxford, UK, 3rd edition, 1999.

M. Hermes-Lima and T. Zenteno-Savín, “Animal response to drastic changes in oxygen availability and physiological oxidative stress,” Comparative Biochemistry and Physiology C, vol. 133, no. 4, pp. 537–556, 2002.

P. A. Olsvik, T. Kristensen, R. Waagbø, K. E. Tollefsen, B. O. Rosseland, and H. Toften, “Effects of hypo- and hyperoxia on transcription levels of five stress genes and the glutathione system in liver of Atlantic cod Gadus morhua,” Journal of Experimental Biology, vol. 209, no. 15, pp. 2893–2901, 2006.

H. Sies, “Biochemistry of oxidative stress,” Angewandte Chemie International Edition, vol. 25, no. 12, pp. 1058–1071, 1986.

M. Tobler, “Does a predatory insect contribute to the divergence between cave- and surface-adapted fish populations?” Biology Letters, vol. 5, no. 4, pp. 506–509, 2009.

E. I. Svensson, F. Eroukhmanoff, and M. Friberg, “Effects of natural and sexual selection on adaptive population divergence and premating isolation in a damselfly,” Evolution, vol. 60, no. 6, pp. 1242–1253, 2006.

P. R. Grant and B. R. Grant, “Pedigrees, assortative mating and speciation in Darwin's finches,” Proceedings of the Royal Society B, vol. 275, no. 1635, pp. 661–668, 2008.

V. J. Parzefall, “Reduction of aggressive behavior patterns in a cave form of Poecilia sphenops (Pisces, Poeciliidae),” Zeitschrift für Tierpsychologie, vol. 35, no. 1, pp. 66–84, 1974.

J. Parzefall, “Changement of behaviour during the evolution of cave animals,” Memoirs of Biospeology, vol. 8, no. 1, pp. 55–62, 1982.

J. Parzefall, U. Gagelmann, and M. Schartl, “Aggressive behaviour and optomotor response in different populations of P. mexicana (Pisces, Poecilidae),” Memoirs of Biospeology, vol. 24, no. 1, pp. 63–69, 1997.

M. Plath, J. Parzefall, and I. Schlupp, “The role of sexual harassment in cave and surface dwelling populations of the Atlantic molly, Poecilia mexicana (Poeciliidae, Teleostei),” Behavioral Ecology and Sociobiology, vol. 54, no. 3, pp. 303–309, 2003.

T. G. Langecker, B. Neumann, C. Hausberg, and J. Parzefall, “Evolution of the optical releasers for aggressive behavior in cave-dwelling Astyanax fasciatus (Teleostei, Characidae),” Behavioural Processes, vol. 34, no. 2, pp. 161–168, 1995.

C. Hausberg, Das aggressionsverhalten von Astyanax fasciatus (Characidae, Teleostei): zur ontogenie, genetik und evolution der epigäischen und hypogäischen form, Doctoral Dissertation, University of Hamburg, Hamburg, Germany, 1995.

R. R. Miller, W. L. Minckley, and S. M. Norris, Freshwater Fishes of México, Chicago University Press, Chicago, Ill, USA, 2005.

D. Bierbach, M. Schulte, N. Herrmann, et al., “Predator-induced changes of female mating preferences: innate and experiential effects,” BMC Evolutionary Biology, vol. 11, p. 190, 2011.

J. Beaugrand, C. Goulet, and D. Payette, “Outcome of dyadic conflict in male green swordtail fish, Xiphophorus helleri: effects of body size and prior dominance,” Animal Behaviour, vol. 41, no. 3, pp. 417–424, 1991.

D. Bierbach, A. Girndt, S. Hamfler et al., “Male fish use prior knowledge about rivals to adjust their mate choice,” Biology Letters, vol. 7, no. 3, pp. 349–351, 2011.

M. R. Morris, L. Gass, and M. J. Ryan, “Assessment and individual recognition of opponents in the pygmy swordtails Xiphophorus nigrensis and X.multilineatus,” Behavioral Ecology and Sociobiology, vol. 37, no. 5, pp. 303–310, 1995.

G. A. Parker, “Assessment strategy and the evolution of fighting behaviour,” Journal of Theoretical Biology, vol. 47, no. 1, pp. 223–243, 1974.

F. A. Huntingford and A. K. Turner, Animal Conflict, Chapman & Hall, London, UK, 1987.

H. D. Rundle, L. Nagel, J. W. Boughman, and D. Schluter, “Natural selection and parallel speciation in sympatric sticklebacks,” Science, vol. 287, no. 5451, pp. 306–308, 2000.

P. Nosil, B. J. Crespi, and C. P. Sandoval, “Host-plant adaptation drives the parallel evolution of reproductive isolation,” Nature, vol. 417, no. 6887, pp. 440–443, 2002.

J. R. Dettman, C. Sirjusingh, L. M. Kohn, and J. B. Anderson, “Incipient speciation by divergent adaptation and antagonistic epistasis in yeast,” Nature, vol. 447, no. 7144, pp. 585–588, 2007.

R. B. Langerhans, M. E. Gifford, and E. O. Joseph, “Ecological speciation in Gambusia fishes,” Evolution, vol. 61, no. 9, pp. 2056–2074, 2007.

P. Nosil, D. J. Funk, and D. Ortiz-Barrientos, “Divergent selection and heterogeneous genomic divergence,” Molecular Ecology, vol. 18, no. 3, pp. 375–402, 2009.

A. P. Hendry, “Selection against migrants contributes to the rapid evolution of ecologically dependent reproductive isolation,” Evolutionary Ecology Research, vol. 6, no. 8, pp. 1219–1236, 2004.

P. Nosil, T. H. Vines, and D. J. Funk, “Perspective: reproductive isolation caused by natural selection against immigrants from divergent habitats,” Evolution, vol. 59, no. 4, pp. 705–719, 2005.

D. R. Matute, C. J. Novak, and J. A. Coyne, “Temperature-based extrinsic reproductive isolation in two species of Drosophila,” Evolution, vol. 63, no. 3, pp. 595–612, 2009.

P. D. Lorch, S. Proulx, L. Rowe, and T. Day, “Condition dependent sexual selection accelerates adaptation by natural selection,” Evolutionary Ecology Research, vol. 5, no. 6, pp. 867–881, 2003.

C. C. Culver, T. C. Kane, and D. W. Fong, Adaptation and Natural Selection in Caves, Harvard University Press, Cambridge, Mass, USA, 1995.

D. L. Bechler, “The evolution of agonistic behavior in amblyopsid fishes,” Behavioral Ecology and Sociobiology, vol. 12, no. 1, pp. 35–42, 1983.

A. Ercolini, R. Berti, A. Cianfanelli, et al., “Aggressive behaviour in Uegitglanis zammaronoi Gianferrari (Clariidae, Siluriformes) an anophtalmic phreatic fish from Somalia,” Monitore Zoologico Italiano, vol. 5, no. 1, pp. 39–56, 1981.

J. Parzefall, J. P. Durand, and B. Richard, “Aggressive behavior of the european cave salamander Proteus anguinus (Proteidae, Urodela),” in Proceedings of the 13th International Congress of Speleology, vol. 2, no. 1, pp. 415–419, 1981.

L. Espinasa, Y. Yamamoto, and W. R. Jeffery, “Non-optical releasers for aggressive behavior in blind and blinded Astyanax (Teleostei, Characidae),” Behavioural Processes, vol. 70, no. 2, pp. 144–148, 2005.

H. Burchards, A. Dölle, and J. Parzefall, “Aggressive behaviour of an epigean population of Astyanax mexicanus (Characidae, Pisces) and some observations of three subterranean populations,” Behavioural Processes, vol. 11, no. 3, pp. 225–235, 1985.

J. Parzefall, “On the hereditary of behavior patterns in cave animals and their epigean relatives,” Bulletin of the National Speleological Society, vol. 47, no. 1, pp. 128–135, 1985.

K. Hüppop, “Adaptation to low food,” in Encyclopedia of Caves, D. C. Culver and W. B. White, Eds., pp. 4–10, Academic Press, Amsterdam, The Netherlands, 2005.

T. G. Langecker, H. Wilkens, and J. Parzefall, “Studies on the trophic structure of an energy-rich Mexican cave (Cueva de las Sardinas) containing sulfurous water,” Memoirs of Biospeology, vol. 23, no. 1, pp. 121–125, 1996.

A. S. Engel, “Chemoautotrophy,” in Encyclopedia of Caves, D. C. Culver and W. B. White, Eds., pp. 90–102, Academic Press, Amsterdam, The Netherlands, 2005.

R. Riesch, I. Schlupp, E. Schleucher, P. Hildenbrand, A. Köhler, and M. Plath, “Reduced starvation resistance and increased metabolic rates in an unusual cave organism: the cave molly (Poecilia mexicana, Poeciliidae),” Bulletin of Fish Biology. In press

F. C. Neat, A. C. Taylor, and F. A. Huntingford, “Proximate costs of fighting in male cichlid fish: the role of injuries and energy metabolism,” Animal Behaviour, vol. 55, no. 4, pp. 875–882, 1998.

M. Tobler, C. M. Franssen, and M. Plath, “Male-biased predation of a cave fish by a giant water bug,” Naturwissenschaften, vol. 95, no. 8, pp. 775–779, 2008.

M. Tobler, M. Plath, R. Riesch et al., “Parasites drive the evolution of MHC diversity but do not promote MHC divergence and speciation in populations of reproductively isolated fish,” In press