TOXIC HYDROGEN SULFIDE AND DARK CAVES: PHENOTYPIC AND GENETIC DIVERGENCE ACROSS TWO ABIOTIC ENVIRONMENTAL GRADIENTS INPOECILIA MEXICANA

Sex-specific life-history trait expression in hybrids of a cave- and surface-dwelling fish (Poecilia mexicana, Poeciliidae)

Evaluating the fitness of hybrids can provide important insights into genetic differences between species or diverging populations. We focused on surface- and cave-ecotypes of the widespread Atlantic molly Poecilia mexicana and raised F1 hybrids of reciprocal crosses to sexual maturity in a common-garden experiment. Hybrids were reared in a fully factorial 2 × 2 design consisting of lighting (light vs. darkness) and resource availability (high vs. low food). We quantified survival, ability to realize their full reproductive potential (i.e., completed maturation for males and 3 consecutive births for females) and essential life-history traits. Compared to the performance of pure cave and surface fish from a previous experiment, F1s had the highest death rate and the lowest proportion of fish that reached their full reproductive potential. We also uncovered an intriguing pattern of sex-specific phenotype expression, because male hybrids expressed cave molly life histories, while female hybrids expressed surface molly life histories. Our results provide evidence for strong selection against hybrids in the cave molly system, but also suggest a complex pattern of sex-specific (opposing) dominance, with certain surface molly genes being dominant in female hybrids and certain cave molly genes being dominant in male hybrids.

Evolutionary Rate Shifts in Coding and Regulatory Regions Underpin Repeated Adaptation to Sulfidic Streams in Poeciliid Fishes

Adaptation to extreme environments often involves the evolution of dramatic physiological changes. To better understand how organisms evolve these complex phenotypic changes, the repeatability and predictability of evolution, and possible constraints on adapting to an extreme environment, it is important to understand how adaptive variation has evolved. Poeciliid fishes represent a particularly fruitful study system for investigations of adaptation to extreme environments due to their repeated colonization of toxic hydrogen sulfide-rich springs across multiple species within the clade. Previous investigations have highlighted changes in the physiology and gene expression in specific species that are thought to facilitate adaptation to hydrogen sulfide-rich springs. However, the presence of adaptive nucleotide variation in coding and regulatory regions and the degree to which convergent evolution has shaped the genomic regions underpinning sulfide tolerance across taxa are unknown. By sampling across seven independent lineages in which nonsulfidic lineages have colonized and adapted to sulfide springs, we reveal signatures of shared evolutionary rate shifts across the genome. We found evidence of genes, promoters, and putative enhancer regions associated with both increased and decreased convergent evolutionary rate shifts in hydrogen sulfide-adapted lineages. Our analysis highlights convergent evolutionary rate shifts in sulfidic lineages associated with the modulation of endogenous hydrogen sulfide production and hydrogen sulfide detoxification. We also found that regions with shifted evolutionary rates in sulfide spring fishes more often exhibited convergent shifts in either the coding region or the regulatory sequence of a given gene, rather than both.

Battle of the Bites: The Effect of Sewage Effluent Exposure on Mosquitofish Biocontrol of Mosquitoes in Residential Louisiana

Mosquitofish, Gambusia affinis, are eponymous larval mosquito predators. Their ability to colonize and survive in habitats that are uninhabitable by other potential predators allows them to naturally manage larval mosquito populations in most ground pools they are present in. However, effluent from residential onsite wastewater treatment systems (OWTSs) appears to limit the presence of fish predators. This is especially problematic in Louisiana, where regulations allow the discharge of OWTS effluent into open drainage conveyances. To determine the effect of effluent on the capacity of mosquitofish for biocontrol in contaminated areas, we assessed the body condition metrics of populations from two effluent-exposed sites and two sites not exposed to effluent, determined the lethal effect of effluent-contaminated drainage water on fish, and measured the prey consumption rates in the presence of effluent. Female fish collected from effluent-impacted sites had a reduced somatic body condition and most females examined displayed masculinized anal fins resembling the male gonopodium structure. This trait was not seen in fish collected from the control sites and has not yet been documented in association with OWTSs or in the state of Louisiana. Fish from the control sites survived at effluent-contaminated water levels < 70%, and the prey clearance rates increased with dilution. Onsite wastewater treatment system effluent has significant effects on both the short- and long-term persistence of mosquitofish, their body composition, reproductive health, and larval mosquito consumption. These effects likely release mosquito larvae from suppression and may increase the threat of mosquito-transmitted pathogens in effluent-contaminated locations.

Wnt11 acts on dermomyotome cells to guide epaxial myotome morphogenesis

The dorsal axial muscles, or epaxial muscles, are a fundamental structure covering the spinal cord and vertebrae, as well as mobilizing the vertebrate trunk. To date, mechanisms underlying the morphogenetic process shaping the epaxial myotome are largely unknown. To address this, we used the medaka zic1/zic4-enhancer mutant Double anal fin (Da), which exhibits ventralized dorsal trunk structures resulting in impaired epaxial myotome morphology and incomplete coverage over the neural tube. In wild type, dorsal dermomyotome (DM) cells reduce their proliferative activity after somitogenesis. Subsequently, a subset of DM cells, which does not differentiate into the myotome population, begins to form unique large protrusions extending dorsally to guide the epaxial myotome dorsally. In Da, by contrast, DM cells maintain the high proliferative activity and mainly form small protrusions. By combining RNA- and ChIP-sequencing analyses, we revealed direct targets of Zic1, which are specifically expressed in dorsal somites and involved in various aspects of development, such as cell migration, extracellular matrix organization, and cell-cell communication. Among these, we identified wnt11 as a crucial factor regulating both cell proliferation and protrusive activity of DM cells. We propose that dorsal extension of the epaxial myotome is guided by a non-myogenic subpopulation of DM cells and that wnt11 empowers the DM cells to drive the coverage of the neural tube by the epaxial myotome.

Population genetics and independently replicated evolution of predator-associated burst speed ecophenotypy in mosquitofish

Many species show replicated ecophenotypy due to recurring patterns of natural selection. Based on the presence or absence of pursuit predators, at least 17 species of fish repeatedly differentiated in body shape in a manner that increases burst swimming speed and the likelihood of predator escape. The predator-associated burst speed (PABS) ecophenotype is characterized by a small head and trunk and enlarged caudal region. Mechanisms promoting replicated phenotype-environment association include selection (without evolution), a single instance of adaptive evolution followed by biased habitat occupation, repeated instances of local adaptation, or adaptive phenotypic plasticity. Common garden rearing of mosquitofish, Gambusia affinis, demonstrated a likely heritable basis for PABS phenotypy, but it is unknown whether populations are otherwise genetically distinct or whether replicated ecophenotypy represents a single or replicated instances of adaptation. To genetically characterize the populations and test hypotheses of single or multiple adaptations, we characterized variation in 12 polymorphic DNA microsatellites in the previously studied G. affinis populations. Populations were genetically distinct by multilocus analysis, exhibited high allelic diversity, and were heterozygote deficient, which effects were attributed to G. affinis's shoaling nature and habitat patchiness. Genetic and phenotypic distances among populations were correlated for non-PABS but not PABS morphology. Multilocus analysis demonstrated ecophenotype polyphyly and scattered multivariate genetic structure which support only the replicated-adaptation model. As all of the diverse tests performed demonstrated lack of congruence between patterns of molecular genetic and PABS differentiation, it is likely that divergent natural selection drove multiple instances of adaptive evolution.

Parallel shifts of visual sensitivity and body coloration in replicate populations of extremophile fish

Visual sensitivity and body pigmentation are often shaped by both natural selection from the environment and sexual selection from mate choice. One way of quantifying the impact of the environment is by measuring how traits have changed after colonization of a novel habitat. To do this, we studied Poecilia mexicana populations that have repeatedly adapted to extreme sulphidic (H₂ S-containing) environments. We measured visual sensitivity using opsin gene expression, as well as body pigmentation, for populations in four independent drainages. Both visual sensitivity and body pigmentation showed significant parallel shifts towards greater medium-wavelength sensitivity and reflectance in sulphidic populations. Altogether we found that sulphidic habitats select for differences in visual sensitivity and pigmentation. Shifts between habitats may be due to both differences in the water's spectral properties and correlated ecological changes.

Population divergence in aggregation and sheltering behaviour in surface- versus cave-adapted Asellus aquaticus (Crustacea: Isopoda)

Aggregation (gathering together) and sheltering (hiding in cover) are basic behaviours that can reduce the risk of predation. However, both behaviours have costs, such as increased competition over resources and high prevalence of contact-spread parasites (aggregation) or lost opportunities for foraging and mating (sheltering). Therefore, variation in these behaviours is expected between populations with varying levels of predation risk. We compared aggregation and sheltering in surface- (various predators) and cave-adapted (no predator) populations of the isopod Asellus aquaticus in a common garden experiment. Given that the cave environment is constantly dark, we also tested for population variation in light-induced behavioural plasticity. Variation in sheltering was explained by habitat type: cave individuals sheltered less than surface individuals. We found high between-population variation in aggregation with or without shelters and their light-induced plasticity, which was not explained by habitat type. Cave individuals decreased (habituation) whereas surface individuals increased sheltering with time (sensitization). We suggest that population variation in sheltering is driven by predation, whereas variation in aggregation must be driven by other, unaccounted environmental factors, in a similar manner to light-induced behavioural plasticity. Based on habituation/sensitization patterns, we suggest that predation-adapted populations are more sensitive to disturbance related to routine laboratory procedures.

Epigenetic inheritance of DNA methylation changes in fish living in hydrogen sulfide-rich springs

Environmental factors can promote phenotypic variation through alterations in the epigenome and facilitate adaptation of an organism to the environment. Although hydrogen sulfide is toxic to most organisms, the fish Poecilia mexicana has adapted to survive in environments with high levels that exceed toxicity thresholds by orders of magnitude. Epigenetic changes in response to this environmental stressor were examined by assessing DNA methylation alterations in red blood cells, which are nucleated in fish. Males and females were sampled from sulfidic and nonsulfidic natural environments; individuals were also propagated for two generations in a nonsulfidic laboratory environment. We compared epimutations between the sexes as well as field and laboratory populations. For both the wild-caught (F0) and the laboratory-reared (F2) fish, comparing the sulfidic and nonsulfidic populations revealed evidence for significant differential DNA methylation regions (DMRs). More importantly, there was over 80% overlap in DMRs across generations, suggesting that the DMRs have stable generational inheritance in the absence of the sulfidic environment. This is an example of epigenetic generational stability after the removal of an environmental stressor. The DMR-associated genes were related to sulfur toxicity and metabolic processes. These findings suggest that adaptation of P. mexicana to sulfidic environments in southern Mexico may, in part, be promoted through epigenetic DNA methylation alterations that become stable and are inherited by subsequent generations independent of the environment.

Mitochondrial physiology and responses to elevated hydrogen sulphide in two isogenic lineages of an amphibious mangrove fish

Hydrogen sulphide (H2S) is toxic and can act as a selective pressure on aquatic organisms, facilitating a wide range of adaptations for life in sulphidic environments. Mangrove rivulus (Kryptolebias marmoratus) inhabit mangrove swamps and have developed high tolerance to environmental H2S. They are hermaphroditic and can self-fertilize, producing distinct isogenic lineages with different sensitivity to H2S. Here, we tested the hypothesis that observed differences in responses to H2S are the result of differences in mitochondrial functions. For this purpose, we performed two experimental series, testing (1) the overall mitochondrial oxidizing capacities and (2) the kinetics of apparent H2S mitochondrial oxidation and inhibition in two distinct lineages of mangrove rivulus, originally collected from Belize and Honduras. We used permeabilized livers from both lineages, measured mitochondrial oxidation, and monitored changes during gradual increases of sulphide. Ultimately, we determined that each lineage has a distinct strategy for coping with elevated H2S, indicating divergences in mitochondrial function and metabolism. The Honduras lineage has higher anaerobic capacity substantiated by higher lactate dehydrogenase activity and higher apparent H2S oxidation rates, likely enabling them to tolerate H2S by escaping aquatic H2S in a terrestrial environment. However, Belize fish have increased cytochrome c oxidase and citrate synthase activities as well as increased succinate contribution to mitochondrial respiration, allowing them to tolerate higher levels of aquatic H2S without inhibition of mitochondrial oxygen consumption. Our study reveals distinct physiological strategies in genetic lineages of a single species, indicating possible genetic and/or functional adaptations to sulphidic environments at the mitochondrial level.

A Novel Body Plan Alters Diversification of Body Shape and Genitalia in Live-Bearing Fish

Major evolutionary innovations can greatly influence subsequent evolution. While many major transitions occurred in the deep past, male live-bearing fishes (family Poeciliidae) more recently evolved a novel body plan. This group possesses a three-region axial skeleton, with one region—the ano-urogenital region—representing a unique body region accommodating male genitalic structures (gonopodial complex). Here we evaluate several hypotheses for the evolution of diversity in this region and examine its role in the evolution of male body shape. Examining Gambusia fishes, we tested a priori predictions for (1) joint influence of gonopodial-complex traits on mating performance, (2) correlated evolution of gonopodial-complex traits at macro- and microevolutionary scales, and (3) predator-driven evolution of gonopodial-complex traits in a post-Pleistocene radiation of Bahamas mosquitofish. We found the length of the sperm-transfer organ (gonopodium) and its placement along the body (gonopodial anterior transposition) jointly influenced mating success, with correlational selection favoring particular trait combinations. Despite these two traits functionally interacting during mating, we found no evidence for their correlated evolution at macro- or microevolutionary scales. In contrast, we did uncover correlated evolution of modified vertebral hemal spines (part of the novel body region) and gonopodial anterior transposition at both evolutionary scales, matching predictions of developmental connections between these components. Developmental linkages in the ano-urogenital region apparently play key roles in evolutionary trajectories, but multiple selective agents likely act on gonopodium length and cause less predictable evolution. Within Bahamas mosquitofish, evolution of hemal-spine morphology, and gonopodial anterior transposition across predation regimes was quite predictable, with populations evolving under high predation risk showing more modified hemal spines with greater modifications and a more anteriorly positioned gonopodium. These changes in the ano-urogenital vertebral region have facilitated adaptive divergence in swimming abilities and body shape between predation regimes. Gonopodium surface area, but not length, evolved as predicted in Bahamas mosquitofish, consistent with a previously suggested tradeoff between natural and sexual selection on gonopodium size. These results provide insight into how restructured body plans offer novel evolutionary solutions. Here, a novel body region—originally evolved to aid sperm transfer—was apparently co-opted to alter whole-organism performance, facilitating phenotypic diversification.

microRNA expression variation as a potential molecular mechanism contributing to adaptation to hydrogen sulphide

microRNAs (miRNAs) are post-transcriptional regulators of gene expression and can play an important role in modulating organismal development and physiology in response to environmental stress. However, the role of miRNAs in mediating adaptation to diverse environments in natural study systems remains largely unexplored. Here, we characterized miRNAs and their expression in Poecilia mexicana, a species of small fish that inhabits both normal streams and extreme environments in the form of springs rich in toxic hydrogen sulphide (H₂ S). We found that P. mexicana has a similar number of miRNA genes as other teleosts. In addition, we identified a large population of mature miRNAs that were differentially expressed between locally adapted populations in contrasting habitats, indicating that miRNAs may contribute to P. mexicana adaptation to sulphidic environments. In silico identification of differentially expressed miRNA-mRNA pairs revealed, in the sulphidic environment, the downregulation of miRNAs predicted to target mRNAs involved in sulphide detoxification and cellular homeostasis, which are pathways essential for life in H₂ S-rich springs. In addition, we found that predicted targets of upregulated miRNAs act in the mitochondria (16.6% of predicted annotated targets), which is the main site of H₂ S toxicity and detoxification, possibly modulating mitochondrial function. Together, the differential regulation of miRNAs between these natural populations suggests that miRNAs may be involved in H₂ S adaptation by promoting functions needed for survival and reducing functions affected by H₂ S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H₂ S. Overall, this study provides a critical stepping-stone towards a comprehensive understanding of the regulatory mechanisms underlying the adaptive variation in gene expression in a natural system.

Hydrogen sulphide sensitivity and tolerance in genetically distinct lineages of a selfing mangrove fish (Kryptolebias marmoratus)

Mangroves are critical marine habitats. High hydrogen sulphide (H2S) is a feature of these important ecosystems and its toxicity creates a challenge for mangrove inhabitants. The mangrove rivulus (Kryptolebias marmoratus) is a selfing, hermaphroditic, amphibious fish that can survive exposure to 1116 μM H2S in the wild. These fish rely on cutaneous respiration for gas and ion exchange when emerged. We hypothesized that the skin surface is fundamentally important in H2S tolerance in these mangrove fish by limiting H2S permeability. To test our hypothesis, we first disrupted the skin surface in one isogenic lineage and measured H2S tolerance and sensitivity. We increased water H2S concentration until emersion as a measure of the ability to sense and react to H2S, which we refer to as sensitivity. We then determined H2S tolerance by preventing emersion and increasing H2S until loss of equilibrium (LOE). The H2S concentration at emersion and LOE were significantly affected by disrupting the skin surface, providing support that the skin is involved in limiting H2S permeability. Capitalizing on their unique reproductive strategy, we used three distinct isogenic lineages to test the hypothesis that there would be genetic differences in H2S sensitivity and tolerance. We found significant differences in emersion concentration only among lineages, suggesting a genetic component to H2S sensitivity but not tolerance. Our study also demonstrated that external skin modifications and avoidance behaviours are two distinct strategies used to tolerate ecologically relevant H2S concentrations and likely facilitate survival in challenging mangrove habitats.

Expression analyses of cave mollies (Poecilia mexicana) reveal key genes involved in the early evolution of eye regression

Eye regression occurs across cave-dwelling populations of many species and is often coupled with a decrease or loss in eye function. Teleost fishes are among the few vertebrates to undergo widespread colonization of caves and often exhibit eye regression with blindness. Cave populations of the poeciliid fish Poecilia mexicana (cave molly) exhibit reduced-albeit functional-eyes, offering the opportunity to investigate partial eye regression. We sequenced eye transcriptomes of cave and surface populations of P. mexicana to identify differentially expressed genes that potentially underlie eye regression in cave mollies. We identified 28 significantly differentially expressed genes, 20 of which were directly related to light sensitivity, eye structure and visual signaling. Twenty-six of these genes were downregulated in cave compared to surface populations. Functional enrichment analysis revealed eye-related gene ontologies that were under-represented in cave mollies. In addition, a set of co-expressed genes related to vision and circadian rhythm was correlated with habitat type (cave versus surface). Our study suggests that differential gene expression plays a key role in the beginning evolutionary stages of eye regression in P. mexicana, shedding further light on regressive evolution in cavefish.

Hydrogen sulphide toxicity and the importance of amphibious behaviour in a mangrove fish inhabiting sulphide-rich habitats

We investigated amphibious behaviour, hydrogen sulphide (H₂S) tolerance, and the mechanism of H₂S toxicity in the amphibious mangrove rivulus (Kryptolebias marmoratus). We found that fish emersed (left water) in response to acutely elevated [H₂S] (~ 130-200 µmol l^-1). The emersion response to H₂S may be influenced by prior acclimation history due to acclimation-induced alterations in gill morphology and/or the density and size of neuroepithelial cells (NECs) on the gills and skin. Thus, we acclimated fish to water (control), H₂S-rich water, or air and tested the hypotheses that acclimation history influences H₂S sensitivity due to acclimation-induced changes in (i) gill surface area and/or (ii) NEC density and/or size. Air-acclimated fish emersed at significantly lower [H₂S] relative to fish acclimated to control or H₂S-rich water, but exhibited no change in gill surface area or in NEC density or size in the gills or skin. Despite possessing exceptional H₂S tolerance, all fish lost equilibrium when unable to emerse from environments containing extremely elevated [H₂S] (2272 ± 46 µmol l^-1). Consequently, we tested the hypothesis that impaired blood oxygen transport (i.e., sulphemoglobin formation) causes H₂S toxicity in amphibious fishes. In vitro exposure of red blood cells to physiologically relevant [H₂S] did not cause a substantial increase in sulphemoglobin formation. We found evidence, however, for an alternative hypothesis that H₂S toxicity is caused by impaired oxidative phosphorylation (i.e., cytochrome c oxidase inhibition). Collectively, our results show that amphibious behaviour is critical for the survival of K. marmoratus in H₂S-rich environments as fish experience impaired oxidative phosphorylation when unable to emerse.

Dissecting a potential spandrel of adaptive radiation: Body depth and pectoral fin ecomorphology coevolve in Lake Malawi cichlid fishes

The evolution of body shape reflects both the ecological factors structuring organismal diversity as well as an organism's underlying anatomy. For instance, body depth in fishes is thought to determine their susceptibility to predators, attractiveness to mates, as well as swimming performance. However, the internal anatomy influencing diversification of body depth has not been extensively examined, and changes in body depth could arise as a by-product of functional changes in other anatomical structures. Using an improved phylogenetic hypothesis for a diverse set of Lake Malawi cichlid fishes, we tested the evolutionary association between body depth and the height of the pectoral girdle. To refine the functional importance of the observed substantial correlation, we also tested the coevolution of pectoral girdle height and pectoral fin area. The extensive coevolution of these traits suggests body depth in fishes like the Lake Malawi cichlids could diverge simply as a by-product of being tightly linked to ecomorphological divergence in other functional morphological structures like the pectoral fins.

Female Choice Undermines the Emergence of Strong Sexual Isolation between Locally Adapted Populations of Atlantic Mollies (Poecilia mexicana)

Divergent selection between ecologically dissimilar habitats promotes local adaptation, which can lead to reproductive isolation (RI). Populations in the Poecilia mexicana species complex have independently adapted to toxic hydrogen sulfide and show varying degrees of RI. Here, we examined the variation in the mate choice component of prezygotic RI. Mate choice tests across drainages (with stimulus males from another drainage) suggest that specific features of the males coupled with a general female preference for yellow color patterns explain the observed variation. Analyses of male body coloration identified the intensity of yellow fin coloration as a strong candidate to explain this pattern, and common-garden rearing suggested heritable population differences. Male sexual ornamentation apparently evolved differently across sulfide-adapted populations, for example because of differences in natural counterselection via predation. The ubiquitous preference for yellow color ornaments in poeciliid females likely undermines the emergence of strong RI, as female discrimination in favor of own males becomes weaker when yellow fin coloration in the respective sulfide ecotype increases. Our study illustrates the complexity of the (partly non-parallel) pathways to divergence among replicated ecological gradients. We suggest that future work should identify the genomic loci involved in the pattern reported here, making use of the increasing genomic and transcriptomic datasets available for our study system.

Testing Wallace's intuition: water type, reproductive isolation and divergence in an Amazonian fish

Alfred Russel Wallace proposed classifying Amazon rivers based on their colour and clarity: white, black and clear water. Wallace also proposed that black waters could mediate diversification and yield distinct fish species. Here, we bring evidence of speciation mediated by water type in the sailfin tetra (Crenuchus spilurus), a fish whose range encompasses rivers of very distinct hydrochemical conditions. Distribution of the two main lineages concords with Wallace's water types: one restricted to the acidic and nutrient-poor waters of the Negro River (herein Rio Negro lineage) and a second widespread throughout the remaining of the species' distribution (herein Amazonas lineage). These lineages occur over a very broad geographical range, suggesting that despite occurring in regions separated by thousands of kilometres, individuals of the distinct lineages fail to occupy each other's habitats, hundreds of metres apart and not separated by physical barrier. Reproductive isolation was assessed in isolated pairs exposed to black-water conditions. All pairs with at least one individual of the lineage not native to black waters showed significantly lower spawning success, suggesting that the water type affected the fitness and contributed to reproductive isolation. Our results endorse Wallace's intuition and highlight the importance of ecological factors in shaping diversity of the Amazon fish fauna.

Pervasive migration across rainforest and sandy coastal plain Aechmea nudicaulis (Bromeliaceae) populations despite contrasting environmental conditions

Understanding the colonization of extreme marginal habitats and the relative roles of space and environment in maintaining peripheral populations remains challenging. Here, we leverage a system of pairs of rainforest and sandy coastal plain communities that allow us to decouple spatial and environmental effects in the population structure and migration rates of the bromeliad Aechmea nudicaulis. Structure and gene flow between populations were estimated from Bayesian clustering and coalescent-based migration models applied to chloroplast sequence and nuclear microsatellite data. Contrary to our initial expectation, the sharp environmental gradient between rainforest and sandy plains does not seem to have affected the colonization and migration dynamics in A. nudicaulis. Our analyses uncover pervasive gene flow between neighbouring habitats in both chloroplast and nuclear data despite the striking differences in environmental conditions. This result is consistent with a scenario of repeated colonization of the sandy coastal plains from forest populations through seed dispersal, as well as the maintenance of gene flow between habitats through pollination. We also recovered a broad north/south population structure that has been found in other Atlantic rainforest groups and possibly reflects older phylogeographic dynamics.

The role of isolation on contrasting phylogeographic patterns in two cave crustaceans

BACKGROUND

The underlying mechanisms and processes that prompt the colonisation of extreme environments, such as caves, constitute major research themes of evolutionary biology and biospeleology. The special adaptations required to survive in subterranean environments (low food availability, hypoxic waters, permanent darkness), and the geographical isolation of caves, nominate cave biodiversity as ideal subjects to answer long-standing questions concerning the interplay amongst adaptation, biogeography, and evolution. The present project aims to examine the phylogeographic patterns exhibited by two sympatric species of surface and cave-dwelling peracarid crustaceans (Asellus aquaticus and Niphargus hrabei), and in doing so elucidate the possible roles of isolation and exaptation in the colonisation and successful adaptation to the cave environment.

RESULTS

Specimens of both species were sampled from freshwater hypogean (cave) and epigean (surface) habitats in Hungary, and additional data from neighbouring countries were sourced from Genbank. Sequencing of mitochondrial and nuclear loci revealed, through haplotype network reconstruction (TCS) and phylogenetic inference, the genetic structure, phylogeographic patterns, and divergence-time estimates of A. aquaticus and N. hrabei surface and cave populations. Contrasting phylogeographic patterns were found between species, with A. aquaticus showing strong genetic differentiation between cave and surface populations and N. hrabei lacking any evidence of genetic structure mediated by the cave environment. Furthermore, N. hrabei populations show very low levels of genetic differentiation throughout their range, which suggests the possibility of recent expansion events over the last few thousand years.

CONCLUSIONS

Isolation by cave environment, rather than distance, is likely to drive the genetic structuring observed between immediately adjacent cave and surface populations of A. aquaticus, a predominantly surface species with only moderate exaptations to subterranean life. For N. hrabei, in which populations exhibit a fully 'cave-adapted' (troglomorphic) phenotype, the lack of genetic structure suggests that subterranean environments do not pose a dispersal barrier for this surface-cave species.

Intrasexual competition enhances reproductive isolation between locally adapted populations

During adaptation to different habitat types, both morphological and behavioral traits can undergo divergent selection. Males often fight for status in dominance hierarchies and rank positions predict reproductive success. Ecotypes with reduced fighting abilities should have low reproductive success when migrating into habitats that harbor ecotypes with superior fighting abilities. Livebearing fishes in the Poecilia mexicana-species complex inhabit not only regular freshwater environments, but also independently colonized sulfidic (H₂S-containing) habitats in three river drainages. In the current study, we found fighting intensities in staged contests to be considerably lower in some but not all sulfidic surface ecotypes and the sulfidic cave ecotype compared with populations from non-sulfidic surface sites. This is perhaps due to selection imposed by H₂S, which hampers oxygen uptake and transport, as well as cellular respiration. Furthermore, migrants from sulfidic habitats may lose fights even if they do not show overall reduced aggressiveness, as physiological performance is likely to be challenged in the non-sulfidic environment to which they are not adapted. To test this hypothesis, we simulated migration of H₂S-adapted males into H₂S-free waters, as well as H₂S-adapted cave-dwelling males into sulfidic surface waters. We found that intruders established dominance less often than resident males, independent of whether or not they showed reduced aggressiveness overall. Our study shows that divergent evolution of male aggressive behavior may also contribute to the maintenance of genetic differentiation in this system and we call for more careful evaluation of male fighting abilities in studies on ecological speciation.

Genome-scale data reveal that endemic Poecilia populations from small sulphidic springs display no evidence of inbreeding

Populations with limited ranges can be highly vulnerable to changes in their environment and are, thus, of high conservation concern. Populations that experience human-induced range reductions are often highly inbred and lack genetic diversity, but it is unknown whether this is also the case for populations with naturally small ranges. The fishes Poecilia sulphuraria (listed as critically endangered) and Poecilia thermalis, which are endemic to small hydrogen sulphide-rich springs in southern Mexico, are examples of such populations with inherently small habitats. We used geometric morphometrics and population genetics to quantify phenotypic and genetic variation within and among two populations of P. sulphuraria and one population of P. thermalis. Principal component analyses revealed phenotypic and genetic differences among the populations. Evidence for inbreeding was low compared to populations that have undergone habitat reduction. The genetic data were also used to infer the demographic history of these populations to obtain estimates for effective population sizes and migration rates. Effective population sizes were large given the small habitats of these populations. Our results imply that these three endemic extremophile populations should each be considered separately for conservation purposes. Additionally, this study suggests that populations in naturally small habitats may have lower rates of inbreeding and higher genetic diversity than expected, and therefore may be better equipped to handle environmental perturbations than anticipated. We caution, however, that the inferred lack of inbreeding and the large effective population sizes could potentially be a result of colonization by genetically diverse ancestors.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

The evolution of vertebrate eye size across an environmental gradient: phenotype does not predict genotype in a Trinidadian killifish

Vertebrates exhibit substantial variation in eye size. Eye size correlates positively with visual capacity and behaviors that enhance fitness, such as predator avoidance. This foreshadows a connection between predation and eye size evolution. Yet, the conditions that favor evolutionary shifts in eye size, besides the well-known role for light availability, are unclear. We tested the influence of predation on the evolution of eye size in Trinidadian killifish, Rivulus hartii. Rivulus are located across a series of communities where they coexist with visually oriented piscivores ("high predation" sites), and no predators ("Rivulus-only" sites). Wild-caught Rivulus from high predation sites generally exhibited a smaller relative eye size than communities that lack predators. Yet, such differences were inconsistent across rivers. Second-generation common garden reared fish revealed repeatable decreases in eye size in Rivulus from high predation sites. We performed additional experiments that tested the importance of light and resources on eye size evolution. Sites that differ in light or resource availability did not differ in eye size. Our results argue that differences in predator-induced mortality underlie genetically-based shifts in vertebrate eye size. We discuss the drivers of eye size evolution in light of the nonparallel trends between the phenotypic and common garden results.

Body shape differences in a pair of closely related Malawi cichlids and their hybrids: Effects of genetic variation, phenotypic plasticity, and transgressive segregation

Phenotypic differences may have genetic and plastic components. Here, we investigated the contributions of both for differences in body shape in two species of Lake Malawi cichlids using wild-caught specimens and a common garden experiment. We further hybridized the two species to investigate the mode of gene action influencing body shape differences and to examine the potential for transgressive segregation. We found that body shape differences between the two species observed in the field are maintained after more than 10 generations in a standardized environment. Nonetheless, both species experienced similar changes in the laboratory environment. Our hybrid cross experiment confirmed that substantial variation in body shape appears to be genetically determined. The data further suggest that the underlying mode of gene action is complex and cannot be explained by simple additive or additive-dominance models. Transgressive phenotypes were found in the hybrid generations, as hybrids occupied significantly more morphospace than both parentals combined. Further, the body shapes of transgressive individuals resemble the body shapes observed in other Lake Malawi rock-dwelling genera. Our findings indicate that body shape can respond to selection immediately, through plasticity, and over longer timescales through adaptation. In addition, our results suggest that hybridization may have played an important role in the diversification of Lake Malawi cichlids through creating new phenotypic variation.

Kin Recognition in a Clonal Fish, Poecilia formosa

Relatedness strongly influences social behaviors in a wide variety of species. For most species, the highest typical degree of relatedness is between full siblings with 50% shared genes. However, this is poorly understood in species with unusually high relatedness between individuals: clonal organisms. Although there has been some investigation into clonal invertebrates and yeast, nothing is known about kin selection in clonal vertebrates. We show that a clonal fish, the Amazon molly (Poecilia formosa), can distinguish between different clonal lineages, associating with genetically identical, sister clones, and use multiple sensory modalities. Also, they scale their aggressive behaviors according to the relatedness to other females: they are more aggressive to non-related clones. Our results demonstrate that even in species with very small genetic differences between individuals, kin recognition can be adaptive. Their discriminatory abilities and regulation of costly behaviors provides a powerful example of natural selection in species with limited genetic diversity.

Extremophile Poeciliidae: multivariate insights into the complexity of speciation along replicated ecological gradients

Background

Replicate population pairs that diverge in response to similar selective regimes allow for an investigation of (a) whether phenotypic traits diverge in a similar and predictable fashion, (b) whether there is gradual variation in phenotypic divergence reflecting variation in the strength of natural selection among populations, (c) whether the extent of this divergence is correlated between multiple character suites (i.e., concerted evolution), and (d) whether gradual variation in phenotypic divergence predicts the degree of reproductive isolation, pointing towards a role for adaptation as a driver of (ecological) speciation. Here, we use poeciliid fishes of the genera Gambusia and Poecilia that have repeatedly evolved extremophile lineages able to tolerate high and sustained levels of toxic hydrogen sulfide (H₂S) to answer these questions.

Results

We investigated evolutionary divergence in response to H₂S in Gambusia spp. (and to a lesser extent Poecilia spp.) using a multivariate approach considering the interplay of life history, body shape, and population genetics (nuclear miscrosatellites to infer population genetic differentiation as a proxy for reproductive isolation). We uncovered both shared and unique patterns of evolution: most extremophile Gambusia predictably evolved larger heads and offspring size, matching a priori predictions for adaptation to sulfidic waters, while variation in adult life histories was idiosyncratic. When investigating patterns for both genera (Gambusia and Poecilia), we found that divergence in offspring-related life histories and body shape were positively correlated across populations, but evidence for individual-level associations between the two character suites was limited, suggesting that genetic linkage, developmental interdependencies, or pleiotropic effects do not explain patterns of concerted evolution. We further found that phenotypic divergence was positively correlated with both environmental H₂S-concentration and neutral genetic differentiation (a proxy for gene flow).

Conclusions

Our results suggest that higher toxicity exerts stronger selection, and that divergent selection appears to constrain gene flow, supporting a scenario of ecological speciation. Nonetheless, progress toward ecological speciation was variable, partially reflecting variation in the strength of divergent selection, highlighting the complexity of selective regimes even in natural systems that are seemingly governed by a single, strong selective agent.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0705-1) contains supplementary material, which is available to authorized users.

Phylogeny and biogeography of Poecilia (Cyprinodontiformes: Poeciliinae) across Central and South America based on mitochondrial and nuclear DNA markers

Poeciliids are a diverse group of small Neotropical fishes, and despite considerable research attention as models in ecology and evolutionary biology, our understanding of their biogeographic and phylogenetic relationships is still limited. We investigated the phylogenetic relationships of South and Central American Poecilia, by examining 2395 base pairs of mitochondrial DNA (ATPase 8/6, COI) and nuclear DNA (S7) for 18 species across six subgenera. Fifty-eight novel sequences were acquired from newly collected specimens and 20 sequences were obtained from previously published material. Analyses of concatenated and partitioned mitochondrial DNA and nuclear DNA sets resulted in a well-supported phylogeny that resolved several monophyletic groups corresponding to previously hypothesized subgenera and species complexes. A divergence-dating analysis supported the hypothesis of the genus Poecilia dispersing into Central America in the early Pliocene (ancestors of Psychropoecilia+Allopoecilia+Mollienesia: 7.3-2.0Mya) from predominantly South America. Subsequently, one lineage (subgenus Allopoecilia: 5.1-1.3Mya) expanded deeper into South America from Lower-Central America, and one lineage expanded from Nuclear-Central America into South America (subgenus Mollienesia: 0.71-0.14Mya). The subgenus Mollienesia diverged into three monophyletic groups that can be identified by nuptial male dorsal fin morphology and inner jaw dentition. A subclade of the unicuspid short-fins (subgenus Mollienesia) was the lineage that expanded into South America during the middle Pleistocene. Species in this subclade are now distributed across northern South America, where they are partially sympatric with Allopoecilia. However the P. (A.) caucana complex was not monophyletic, with P. (A.) wandae clustering in the Mollienesia subclade that expanded into South America. It is apparent that characters (body size, scale count, pigmentation, and gonopodium morphology) used to define the P. (A.) caucana complex are homoplastic and suggestive of rapid convergence in northern South America. Our improved taxon sampling and divergence-time calibration allowed for insights into the timing and direction of dispersals, and provides an improved understanding of the biogeographic history of an enigmatic group of fishes. Furthermore, we provided strong evidence for the monophyly of the subgenus Mollienesia and further substantiated its species complexes; therefore, we advise a taxonomic re-evaluation for the P. (A.) caucana complex to maintain monophyly of both Mollienesia and Allopoecilia.

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide-toxic habitats of a neotropical fish (Poecilia mexicana)

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H2 S) and compared two population pairs of sulphide-adapted and ancestral fish by sequencing population pools of >200 individuals (Pool-Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection-mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide-adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of Jcorr  = 0.0032 at the level of SNPs, divergent genome regions (Jcorr  = 0.0061) and genes therein (Jcorr  = 0.0091). At the level of metabolic pathways, the overlap was Jcorr  = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Gould's view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.

Sex-specific local life-history adaptation in surface- and cave-dwelling Atlantic mollies (Poecilia mexicana)

Cavefishes have long been used as model organisms showcasing adaptive diversification, but does adaptation to caves also facilitate the evolution of reproductive isolation from surface ancestors? We raised offspring of wild-caught surface- and cave-dwelling ecotypes of the neotropical fish Poecilia mexicana to sexual maturity in a 12-month common garden experiment. Fish were raised under one of two food regimes (high vs. low), and this was crossed with differences in lighting conditions (permanent darkness vs. 12:12 h light:dark cycle) in a 2 × 2 factorial design, allowing us to elucidate potential patterns of local adaptation in life histories. Our results reveal a pattern of sex-specific local life-history adaptation: Surface molly females had the highest fitness in the treatment best resembling their habitat of origin (high food and a light:dark cycle), and suffered from almost complete reproductive failure in darkness, while cave molly females were not similarly affected in any treatment. Males of both ecotypes, on the other hand, showed only weak evidence for local adaptation. Nonetheless, local life-history adaptation in females likely contributes to ecological diversification in this system and other cave animals, further supporting the role of local adaptation due to strong divergent selection as a major force in ecological speciation.

Linking stream ecology with morphological variability in a native freshwater fish from semi-arid Australia

Environmental variation is a potent force affecting phenotypic expression. While freshwater fishes have provided a compelling example of the link between the environment and phenotypic diversity, few studies have been conducted with arid-zone fishes, particularly those that occur in geographically isolated regions where species typically inhabit intermittent and ephemeral creeks. We investigated morphological variation of a freshwater fish (the western rainbowfish, Melanotaenia australis) inhabiting creeks in the Pilbara region of northwest Australia to determine whether body shape variation correlated with local environmental characteristics, including water velocity, habitat complexity, predator presence, and food availability. We expected that the geographic isolation of creeks within this arid region would result in habitat-specific morphological specializations. We used landmark-based geometric morphometrics to quantify the level of morphological variability in fish captured from 14 locations within three distinct subcatchments of a major river system. Western rainbowfish exhibited a range of morphologies, with variation in body depth accounting for a significant proportion (>42%) of the total variance in shape. Sexual dimorphism was also apparent, with males displaying deeper bodies than females. While the measured local habitat characteristics explained little of the observed morphological variation, fish displayed significant morphological differentiation at the level of the subcatchment. Local adaptation may partly explain the geographic patterns of body shape variation, but fine-scale genetic studies are required to disentangle the effects of genetic differentiation from environmentally determined phenotypic plasticity in body shape. Developing a better understanding of environment-phenotype relationships in species from arid regions will provide important insights into ecological and evolutionary processes in these unique and understudied habitats.

Do stressful conditions make adaptation difficult? Guppies in the oil-polluted environments of southern Trinidad

The ability of populations to rapidly adapt to new environments will determine their future in an increasingly human-modified world. Although meta-analyses do frequently uncover signatures of local adaptation, they also reveal many exceptions. We suggest that particular constraints on local adaptation might arise when organisms are exposed to novel stressors, such as anthropogenic pollution. To inform this possibility, we studied the extent to which guppies (Poecilia reticulata) show local adaptation to oil pollution in southern Trinidad. Neutral genetic markers revealed that paired populations in oil-polluted versus not-polluted habitats diverged independently in two different watersheds. Morphometrics revealed some divergence (particularly in head shape) between these environments, some of which was parallel between rivers. Reciprocal transplant experiments in nature, however, found little evidence of local adaptation based on survival and growth. Moreover, subsequent laboratory experiments showed that the two populations from oil-polluted sites showed only weak local adaptation even when compared to guppies from oil-free northern Trinidad. We conclude that guppies show little local adaptation to oil pollution, which might result from the challenges associated with adaptation to particularly stressful environments. It might also reflect genetic drift owing to small population sizes and/or high gene flow between environments.

Naturally occurring variation in tadpole morphology and performance linked to predator regime

Divergent natural selection drives a considerable amount of the phenotypic and genetic variation observed in natural populations. For example, variation in the predator community can generate conflicting selection on behavioral, life-history, morphological, and performance traits. Differences in predator regime can subsequently increase phenotypic and genetic variations in the population and result in the evolution of reproductive barriers (ecological speciation) or phenotypic plasticity. We evaluated morphology and swimming performance in field collected Bronze Frog larvae (Lithobates clamitans) in ponds dominated by predatory fish and those dominated by invertebrate predators. Based on previous experimental findings, we hypothesized that tadpoles from fish-dominated ponds would have small bodies, long tails, and large tail muscles and that these features would facilitate fast-start speed. We also expected to see increased tail fin depth (i.e., the tail-lure morphology) in tadpoles from invertebrate-dominated ponds. Our results support our expectations with respect to morphology in affecting swimming performance of tadpoles in fish-dominated ponds. Furthermore, it is likely that divergent natural selection is playing a role in the diversification on morphology and locomotor performance in this system.

From local adaptation to ecological speciation in copepod populations from neighboring lakes

Continental copepods have been derived from several independent invasive events from the sea, but the subsequent evolutionary processes that account for the current diversity in lacustrine environments are virtually unknown. Salinity is highly variable among lakes and constitutes a source of divergent selection driving potential reproductive isolation. We studied four populations of the calanoid copepod Leptodiaptomus cf. sicilis inhabiting four neighboring lakes with a common history (since the Late Pleistocene) located in the Oriental Basin, Mexico; one lake is shallow and varies in salinity periodically (1.4-10 g L(-1)), while three are deep and permanent, with constant salinity (0.5, 1.1 and 6.5 g L(-1), respectively). We hypothesized that (1) these populations belong to a different species than L. sicilis sensu stricto and (2) are experiencing ecologically based divergence due to salinity differences. We assessed morphological and molecular (mtDNA) COI variation, as well as fitness differences and tests of reproductive isolation. Although relationships of the Mexican populations with L. sicilis s.s. could not be elucidated, we identified a clear pattern of divergent selection driven by salinity conditions. The four populations can still be considered a single biological species (sexual recognition and hybridization are still possible in laboratory conditions), but they have diverged into at least three different phenotypes: two locally adapted, specialized in the lakes of constant salinity (saline vs. freshwater), and an intermediate generalist phenotype inhabiting the temporary lake with fluctuating salinity. The specialized phenotypes are poorly suited as migrants, so prezygotic isolation due to immigrant inviability is highly probable. This implication was supported by molecular evidence that showed restricted gene flow, persistence of founder events, and a pattern of allopatric fragmentation. This study showed how ecologically based divergent selection may explain diversification patterns in lacustrine copepods.

Morphological evolution of coexisting amphipod species pairs from sulfidic caves suggests competitive interactions and character displacement, but no environmental filtering and convergence

Phenotypically similar species coexisting in extreme environments like sulfidic water are subject to two opposing eco-evolutionary processes: those favoring similarity of environment-specific traits, and those promoting differences of traits related to resource use. The former group of processes includes ecological filtering and convergent or parallel evolution, the latter competitive exclusion, character displacement and divergent evolution. We used a unique eco-evolutionary study system composed of two independent pairs of coexisting amphipod species (genus Niphargus) from the sulfidic caves Movile in Romania and Frasassi in Italy to study the relative contribution and interaction of both processes. We looked at the shape of the multifunctional ventral channel as a trait ostensibly related to oxygenation and sulfide detoxification, and at body size as a resource-related trait. Phylogenetic analysis suggests that the sulfidic caves were colonized separately by ancestors of each species. Species within pairs were more dissimilar in their morphology than expected according to a null model based on regional species pool. This might indicate competitive interactions shaping the morphology of these amphipod species. Moreover, our results suggest that the shape of the ventral channel is not subject to long-term convergent selection or to the process of environmental filtering, and as such probably does not play a role in sulfide tolerance. Nevertheless, the ancestral conditions reconstructed using the comparative method tended to be more similar than null-model expectations. This shift in patterns may reflect a temporal hierarchy of eco-evolutionary processes, in which initial environmental filtering became later on superseded by character displacement or other competition-driven divergent evolutionary processes.

Reduction of Energetic Demands through Modification of Body Size and Routine Metabolic Rates in Extremophile Fish

Variation in energy availability or maintenance costs in extreme environments can exert selection for efficient energy use, and reductions in organismal energy demand can be achieved in two ways: reducing body mass or metabolic suppression. Whether long-term exposure to extreme environmental conditions drives adaptive shifts in body mass or metabolic rates remains an open question. We studied body size variation and variation in routine metabolic rates in locally adapted populations of extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulfide-rich springs and caves. We quantified size distributions and routine metabolic rates in wild-caught individuals from four habitat types. Compared with ancestral populations in nonsulfidic surface habitats, extremophile populations were characterized by significant reductions in body size. Despite elevated metabolic rates in cave fish, the body size reduction precipitated in significantly reduced energy demands in all extremophile populations. Laboratory experiments on common garden-raised fish indicated that elevated routine metabolic rates in cave fish likely have a genetic basis. The results of this study indicate that adaptation to extreme environments directly impacts energy metabolism, with fish living in cave and sulfide spring environments expending less energy overall during routine metabolism.

Body size and allometric shape variation in the molly Poecilia vivipara along a gradient of salinity and predation

Background

Phenotypic diversity among populations may result from divergent natural selection acting directly on traits or via correlated responses to changes in other traits. One of the most frequent patterns of correlated response is the proportional change in the dimensions of anatomical traits associated with changes in growth or absolute size, known as allometry. Livebearing fishes subject to predation gradients have been shown to repeatedly evolve larger caudal peduncles and smaller cranial regions under high predation regimes. Poecilia vivipara is a livebearing fish commonly found in coastal lagoons in the north of the state of Rio de Janeiro, Brazil. Similar to what is observed in other predation gradients, lagoons inhabited by P. vivipara vary in the presence of piscivorous fishes; contrary to other poeciliid systems, populations of P. vivipara vary greatly in body size, which opens the possibility of strong allometric effects on shape variation. Here we investigated body shape diversification among six populations of P. vivipara along a predation gradient and its relationship with allometric trajectories within and among populations.

Results

We found substantial body size variation and correlated shape changes among populations. Multivariate regression analysis showed that size variation among populations accounted for 66% of shape variation in females and 38% in males, suggesting that size is the most important dimension underlying shape variation among populations of P. vivipara in this system. Changes in the relative sizes of the caudal peduncle and cranial regions were only partly in line with predictions from divergent natural selection associated with predation regime.

Conclusions

Our results suggest the possibility that adaptive shape variation among populations has been partly constrained by allometry in P. vivipara. Processes governing body size changes are therefore important in the diversification of this species. We conclude that in species characterized by substantial among-population differences in body size, ignoring allometric effects when investigating divergent natural selection’s role in phenotypic diversification might not be warranted.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0251-7) contains supplementary material, which is available to authorized users.

Patterns of Macroinvertebrate and Fish Diversity in Freshwater Sulphide Springs

Extreme environments are characterised by the presence of physicochemical stressors and provide unique study systems to address problems in evolutionary ecology research. Sulphide springs provide an example of extreme freshwater environments; because hydrogen sulphide’s adverse physiological effects induce mortality in metazoans even at micromolar concentrations. Sulphide springs occur worldwide, but while microbial communities in sulphide springs have received broad attention, little is known about macroinvertebrates and fish inhabiting these toxic environments. We reviewed qualitative occurrence records of sulphide spring faunas on a global scale and present a quantitative case study comparing diversity patterns in sulphidic and adjacent non-sulphidic habitats across replicated river drainages in Southern Mexico. While detailed studies in most regions of the world remain scarce, available data suggests that sulphide spring faunas are characterised by low species richness. Dipterans (among macroinvertebrates) and cyprinodontiforms (among fishes) appear to dominate the communities in these habitats. At least in fish, there is evidence for the presence of highly endemic species and populations exclusively inhabiting sulphide springs. We provide a detailed discussion of traits that might predispose certain taxonomic groups to colonize sulphide springs, how colonizers subsequently adapt to cope with sulphide toxicity, and how adaptation may be linked to speciation processes.

Microhabitat use, population densities, and size distributions of sulfur cave-dwelling Poecilia mexicana

The Cueva del Azufre in Tabasco, Mexico, is a nutrient-rich cave and its inhabitants need to cope with high levels of dissolved hydrogen sulfide and extreme hypoxia. One of the successful colonizers of this cave is the poeciliid fish Poecilia mexicana, which has received considerable attention as a model organism to examine evolutionary adaptations to extreme environmental conditions. Nonetheless, basic ecological data on the endemic cave molly population are still missing; here we aim to provide data on population densities, size class compositions and use of different microhabitats. We found high overall densities in the cave and highest densities at the middle part of the cave with more than 200 individuals per square meter. These sites have lower H₂S concentrations compared to the inner parts where most large sulfide sources are located, but they are annually exposed to a religious harvesting ceremony of local Zoque people called La Pesca. We found a marked shift in size/age compositions towards an overabundance of smaller, juvenile fish at those sites. We discuss these findings in relation to several environmental gradients within the cave (i.e., differences in toxicity and lighting conditions), but we also tentatively argue that the annual fish harvest during a religious ceremony (La Pesca) locally diminishes competition (and possibly, cannibalism by large adults), which is followed by a phase of overcompensation of fish densities.

AFLPs and mitochondrial haplotypes reveal local adaptation to extreme thermal environments in a freshwater gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.

Parallel evolution of cox genes in H2S-tolerant fish as key adaptation to a toxic environment

Populations that repeatedly adapt to the same environmental stressor offer a unique opportunity to study adaptation, especially if there are a priori predictions about the genetic basis underlying phenotypic evolution. Hydrogen sulphide (H2S) blocks the cytochrome-c oxidase complex (COX), predicting the evolution of decreased H2S susceptibility of the COX in three populations in the Poecilia mexicana complex that have colonized H2S-containing springs. Here, we demonstrate that decreased H2S susceptibility of COX evolved in parallel in two sulphide lineages, as evidenced by shared amino acid substitutions in cox1 and cox3 genes. One of the shared substitutions likely triggers conformational changes in COX1 blocking the access of H2S. In a third sulphide population, we detect no decreased H2S susceptibility of COX, suggesting that H2S resistance is achieved through another mechanism. Our study thus demonstrates that even closely related lineages follow both parallel and disparate molecular evolutionary paths to adaptation in response to the same selection pressure.

H2S exposure elicits differential expression of candidate genes in fish adapted to sulfidic and non-sulfidic environments

Disentangling the effects of plasticity, genetic variation, and their interactions on organismal responses to environmental stressors is a key objective in ecological physiology. We quantified the expression of five candidate genes in response to hydrogen sulfide (H2S) exposure in fish (Poecilia mexicana, Poeciliidae) from a naturally sulfide-rich environment as well as an ancestral, non-sulfidic population to test for constitutive and environmentally dependent population differences in gene expression patterns. Common garden raised individuals that had never encountered environmental H2S during their lifetime were subjected to short or long term H2S exposure treatments or respective non-sulfidic controls. The expression of genes involved in responses to H2S toxicity (cytochrome c oxidase, vascular endothelial growth factor, and cytochrome P450-2J6), H2S detoxification (sulfide:quinone oxidoreductase), and endogenous H2S production (cystathionine γ lyase) was determined in both gill and liver tissues by real time PCR. The results indicated complex changes in expression patterns that--depending on the gene--not only differed between organs and populations, but also on the type of H2S exposure. Populations differences, both constitutive and H2S exposure dependent (i.e., plastic), in gene expression were particularly evident for sulfide:quinone oxidoreductase, vascular endothelial growth factor, and to a lesser degree for cytochrome P450-2J6. Our study uncovered putatively adaptive modifications in gene regulation that parallel previously documented adaptive changes in phenotypic traits.

Selection from parasites favours immunogenetic diversity but not divergence among locally adapted host populations

The unprecedented polymorphism in the major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection from parasites. However, do parasites also drive divergence at MHC loci between host populations, or do the effects of balancing selection maintain similarities among populations? We examined MHC variation in populations of the livebearing fish Poecilia mexicana and characterized their parasite communities. Poecilia mexicana populations in the Cueva del Azufre system are locally adapted to darkness and the presence of toxic hydrogen sulphide, representing highly divergent ecotypes or incipient species. Parasite communities differed significantly across populations, and populations with higher parasite loads had higher levels of diversity at class II MHC genes. However, despite different parasite communities, marked divergence in adaptive traits and in neutral genetic markers, we found MHC alleles to be remarkably similar among host populations. Our findings indicate that balancing selection from parasites maintains immunogenetic diversity of hosts, but this process does not promote MHC divergence in this system. On the contrary, we suggest that balancing selection on immunogenetic loci may outweigh divergent selection causing divergence, thereby hindering host divergence and speciation. Our findings support the hypothesis that balancing selection maintains MHC similarities among lineages during and after speciation (trans-species evolution).

Morphological variation in vanishing Mexican desert fishes of the genus Characodon (Goodeidae)

This study investigated a rapidly vanishing group of fishes in the genus Characodon (including Characodon lateralis and Characodon audax) from the upper Rio Mezquital drainage in Mexico. Using specimens from museum collections, morphological variation was assessed to quantify body shape differentiation among historic (i.e. extirpated and extant) collection sites. In both sexes, body shape (particularly head shape, the proportion of the caudal peduncle and the position and size of dorsal and anal fins) varied significantly among populations and species. Variation among collection sites could at least partially be attributed to geography, as the presence of distinct hydrographic units and a major waterfall coincided with major body shape differences. These results are discussed in the light of previously published molecular genetic analyses, as they have direct implications for taxonomic problems and the need for conservation measures for these endangered fishes.