Genomic resources for a model in adaptation and speciation research: characterization of the Poecilia mexicana transcriptome

Developmental Regulation of Transcription in Touriga Nacional Berries under Deficit Irrigation

Grapevine (Vitis vinifera L.) is one of the most economically important crops worldwide, especially due to the economic relevance of wine production. Abiotic stress, such as drought, may contribute to low yield, shifts in quality, and important economic loss. The predicted climate change phenomena point to warmer and dryer Mediterranean environmental conditions; as such, it is paramount to study the effects of abiotic stress on grapevine performance. Deficit irrigation systems are applied to optimize water use efficiency without compromising berry quality. In this research, the effect of two deficit irrigation strategies, sustained deficit irrigation (SDI) and regulated deficit irrigation (RDI), in the grape berry were assessed. The effects of different levels of drought were monitored in Touriga Nacional at key stages of berry development (pea size, véraison, and full maturation) through RNA-Seq transcriptome analysis and by specific differentially expressed genes (DEGs) monitoring through RT-qPCR. Handy datasets were obtained by bioinformatics analysis of raw RNA-Seq results. The dominant proportion of transcripts was mostly regulated by development, with véraison showing more upregulated transcripts. Results showed that primary metabolism is the functional category more severely affected under water stress. Almost all DEGs selected for RT-qPCR were significantly upregulated in full maturation and showed the highest variability at véraison and the lowest gene expression values in the pea size stage.

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.

Embryonic origin and genetic basis of cave associated phenotypes in the isopod crustacean Asellus aquaticus

Characteristics common to animals living in subterranean environments include the reduction or absence of eyes, lessened pigmentation and enhanced sensory systems. How these characteristics have evolved is poorly understood for the majority of cave dwelling species. In order to understand the evolution of these changes, this study uses an invertebrate model system, the freshwater isopod crustacean, Asellus aquaticus, to examine whether adult differences between cave and surface dwelling individuals first appear during embryonic development. We hypothesized that antennal elaboration, as well as eye reduction and pigment loss, would be apparent during embryonic development. We found that differences in pigmentation, eye formation, and number of segments of antenna II were all present by the end of embryonic development. In addition, we found that cave and surface hatchlings do not significantly differ in the relative size of antenna II and the duration of embryonic development. To investigate whether the regions responsible for eye and pigment differences could be genetically linked to differences in article number, we genotyped F2 hybrids for the four previously mapped genomic regions associated with eye and pigment differences and phenotyped these F2 hybrids for antenna II article number. We found that the region previously known to be responsible for both presence versus absence of pigment and eye size also was significantly associated with article number. Future experiments will address whether pleiotropy and/or genetic linkage play a role in the evolution of cave characteristics in Asellus aquaticus.

Consistent individual differences in associative learning speed are not linked to boldness in female Atlantic mollies

Recent studies on consistent individual differences in behavioural tendencies (animal personality) raised the question of whether individual differences in cognitive abilities can be linked to certain personality types. We tested female Atlantic mollies (Poecilia mexicana) in two different classical conditioning experiments. For the first time, we provide evidence for highly consistent individual differences in associative learning speed in fish. We characterized the same individuals for boldness in two experimental situations (latency to emerge from shelter and freezing time after a simulated predator attack) and found high behavioural repeatability. When we tested for a potential correlation between associative learning speed and boldness, however, there was no evidence for a link between them. Our study design included several steps to avoid typical pitfalls of disadvantaging shy individuals during learning tests. We caution that other experimental studies may have suffered from erroneous interpretations due to a more cautious coping style of shy individuals in the respective setup used to assess learning.

The aquatic animals' transcriptome resource for comparative functional analysis

Background

Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome.

Results

To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis.

Conclusion

In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4463-x) contains supplementary material, which is available to authorized users.

Molecular evolution and expression of oxygen transport genes in livebearing fishes (Poeciliidae) from hydrogen sulfide rich springs

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S—along with related lineages from non-sulfidic environments—to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

Comparative transcriptomics of cyprinid minnows and carp in a common wild setting: a resource for ecological genomics in freshwater communities

Comparative transcriptomics can now be conducted on organisms in natural settings, which has greatly enhanced understanding of genome–environment interactions. Here, we demonstrate the utility and potential pitfalls of comparative transcriptomics of wild organisms, with an example from three cyprinid fish species (Teleostei:Cypriniformes). We present extensively filtered and annotated transcriptome assemblies that provide a valuable resource for studies of genome evolution (e.g. polyploidy), ecological and morphological diversification, speciation, and shared and unique responses to environmental variation in cyprinid fishes. Our results and analyses address the following points: (i) ‘essential developmental genes’ are shown to be ubiquitously expressed in a diverse suite of tissues across later ontogenetic stages (i.e. juveniles and adults), making these genes are useful for assessing the quality of transcriptome assemblies, (ii) the influence of microbiomes and other exogenous DNA, (iii) potentially novel, species-specific genes, and (iv) genomic rearrangements (e.g. whole genome duplication). The data we present provide a resource for future comparative work in cypriniform fishes and other taxa across a variety of sub-disciplines, including stress response, morphological diversification, community ecology, ecotoxicology, and climate change.

The gonadal transcriptome of the unisexual Amazon molly Poecilia formosa in comparison to its sexual ancestors, Poecilia mexicana and Poecilia latipinna

Background

The unisexual Amazon molly (Poecilia formosa) originated from a hybridization between two sexual species, the sailfin molly (Poecilia latipinna) and the Atlantic molly (Poecilia mexicana). The Amazon molly reproduces clonally via sperm-dependent parthenogenesis (gynogenesis), in which the sperm of closely related species triggers embryogenesis of the apomictic oocytes, but typically does not contribute genetic material to the next generation. We compare for the first time the gonadal transcriptome of the Amazon molly to those of both ancestral species, P. mexicana and P. latipinna.

Results

We sequenced the gonadal transcriptomes of the P. formosa and its parental species P. mexicana and P. latipinna using Illumina RNA-sequencing techniques (paired-end, 100 bp). De novo assembly of about 50 million raw read pairs for each species was performed using Trinity, yielding 106,922 transcripts for P. formosa, 115,175 for P. latipinna, and 133,025 for P. mexicana after eliminating contaminations. On the basis of sequence similarity comparisons to other teleost species and the UniProt databases, functional annotation, and differential expression analysis, we demonstrate the similarity of the transcriptomes among the three species. More than 40% of the transcripts for each species were functionally annotated and about 70% were assigned to orthologous genes of a closely related species. Differential expression analysis between the sexual and unisexual species uncovered 2035 up-regulated and 564 down-regulated genes in P. formosa. This was exemplary validated for six genes by qRT-PCR.

Conclusions

We identified more than 130 genes related to meiosis and reproduction within the apomictically reproducing P. formosa. Overall expression of these genes seems to be down-regulated in the P. formosa transcriptome compared to both ancestral species (i.e., 106 genes down-regulated, 29 up-regulated). A further 35 meiosis and reproduction related genes were not found in the P. formosa transcriptome, but were only expressed in the sexual species. Our data support the hypothesis of general down-regulation of meiosis-related genes in the apomictic Amazon molly. Furthermore, the obtained dataset and identified gene catalog will serve as a resource for future research on the molecular mechanisms behind the reproductive mode of this unisexual species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4382-2) contains supplementary material, which is available to authorized users.

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H₂ S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H₂ S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H₂ S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H₂ S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H₂ S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

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.

Deep analysis of wild Vitis flower transcriptome reveals unexplored genome regions associated with sex specification

RNA-seq of Vitis during early stages of bud development, in male, female and hermaphrodite flowers, identified new loci outside of annotated gene models, suggesting their involvement in sex establishment. The molecular mechanisms responsible for flower sex specification remain unclear for most plant species. In the case of V. vinifera ssp. vinifera, it is not fully understood what determines hermaphroditism in the domesticated subspecies and male or female flowers in wild dioecious relatives (Vitis vinifera ssp. sylvestris). Here, we describe a de novo assembly of the transcriptome of three flower developmental stages from the three Vitis vinifera flower types. The validation of de novo assembly showed a correlation of 0.825. The main goals of this work were the identification of V. v. sylvestris exclusive transcripts and the characterization of differential gene expression during flower development. RNA from several flower developmental stages was used previously to generate Illumina sequence reads. Through a sequential de novo assembly strategy one comprehensive transcriptome comprising 95,516 non-redundant transcripts was assembled. From this dataset 81,064 transcripts were annotated to V. v. vinifera reference transcriptome and 11,084 were annotated against V. v. vinifera reference genome. Moreover, we found 3368 transcripts that could not be mapped to Vitis reference genome. From all the non-redundant transcripts that were assembled, bioinformatics analysis identified 133 specific of V. v. sylvestris and 516 transcripts differentially expressed among the three flower types. The detection of transcription from areas of the genome not currently annotated suggests active transcription of previously unannotated genomic loci during early stages of bud development.

Transcriptome sequencing and simple sequence repeat marker development for three Macaronesian endemic plant species

PREMISE OF THE STUDY

Oceanic islands offer unparalleled opportunities to investigate evolutionary processes such as adaptation and speciation. However, few genomic resources are available for oceanic island endemics. In this study, we publish transcriptome sequences from three Macaronesian endemic plant species (Argyranthemum broussonetii [Asteraceae], Descurainia bourgaeana [Brassicaceae], and Echium wildpretii [Boraginaceae]) that are representative of lineages that have radiated in the region. In addition, the utility of transcriptome data for marker development is demonstrated.

METHODS AND RESULTS

Transcriptomes from the three plant species were sequenced, assembled, and annotated. Between 1972 and 2282 simple sequence repeats (SSRs) were identified for each taxon. Primers were designed and tested for 30 of the candidate SSRs identified in Argyranthemum, of which 12 amplified well across three species and eight were polymorphic.

CONCLUSIONS

We demonstrate here that a single transcriptome sequence is sufficient to identify hundreds of polymorphic SSR markers. The SSRs are applicable to a wide range of questions relating to the evolution of island lineages.

Mechanisms Underlying Adaptation to Life in Hydrogen Sulfide-Rich Environments

Hydrogen sulfide (H2S) is a potent toxicant interfering with oxidative phosphorylation in mitochondria and creating extreme environmental conditions in aquatic ecosystems. The mechanistic basis of adaptation to perpetual exposure to H2S remains poorly understood. We investigated evolutionarily independent lineages of livebearing fishes that have colonized and adapted to springs rich in H2S and compared their genome-wide gene expression patterns with closely related lineages from adjacent, nonsulfidic streams. Significant differences in gene expression were uncovered between all sulfidic and nonsulfidic population pairs. Variation in the number of differentially expressed genes among population pairs corresponded to differences in divergence times and rates of gene flow, which is consistent with neutral drift driving a substantial portion of gene expression variation among populations. Accordingly, there was little evidence for convergent evolution shaping large-scale gene expression patterns among independent sulfide spring populations. Nonetheless, we identified a small number of genes that was consistently differentially expressed in the same direction in all sulfidic and nonsulfidic population pairs. Functional annotation of shared differentially expressed genes indicated upregulation of genes associated with enzymatic H2S detoxification and transport of oxidized sulfur species, oxidative phosphorylation, energy metabolism, and pathways involved in responses to oxidative stress. Overall, our results suggest that modification of processes associated with H2S detoxification and toxicity likely complement each other to mediate elevated H2S tolerance in sulfide spring fishes. Our analyses allow for the development of novel hypotheses about biochemical and physiological mechanisms of adaptation to extreme environments.

Local adaptation of Gymnocypris przewalskii (Cyprinidae) on the Tibetan Plateau

Divergent selection among environments affects species distributions and can lead to speciation. In this article, we investigated the transcriptomes of two ecotypes of scaleless carp (Gymnocypris przewalskii przewalskii and G. p. ganzihonensis) from the Tibetan Plateau. We used a transcriptome sequencing approach to screen approximately 250,000 expressed sequence tags (ESTs) from the gill and kidney tissues of twelve individuals from the Ganzi River and Lake Qinghai to understand how this freshwater fish has adapted to an ecological niche shift from saline to freshwater. We identified 9,429 loci in the gill transcriptome and 12,034 loci in the kidney transcriptome with significant differences in their expression, of which 242 protein-coding genes exhibited strong positive selection (K_a/K_s > 1). Many of the genes are involved in ion channel functions (e.g., Ca²⁺-binding proteins), immune responses (e.g., nephrosin) or cellular water absorption functions (e.g., aquaporins). These results have potentially broad importance in understanding shifts from saline to freshwater habitats. Furthermore, this study provides the first transcriptome of G. przewalskii, which will facilitate future ecological genomics studies and aid in the identification of genes underlying adaptation and incipient ecological speciation.

Genomics: moving behavioural ecology beyond the phenotypic gambit

Researchers studying the adaptive significance of behaviour typically assume that genetic mechanisms will not inhibit evolutionary trajectories, an assumption commonly known as the 'phenotypic gambit'. Although the phenotypic gambit continues to be a useful heuristic for behavioural ecology, here we discuss how genomic methods provide new tools and conceptual approaches that are relevant to behavioural ecology. We first describe how the concept of a genetic toolkit for behaviour can allow behavioural ecologists to synthesize both genomic and ecological information when assessing behavioural adaptation. Then we show how gene expression profiles can be viewed as complex phenotypic measurements, used to (1) predict behaviour, (2) evaluate phenotypic plasticity and (3) devise methods to manipulate behaviour in order to test adaptive hypotheses. We propose that advances in genomics and bioinformatics may allow researchers to overcome some of the logistical obstacles that motivated the inception of the phenotypic gambit. Behavioural ecology and genomics are mutually informative, providing potential synergy that could lead to powerful advances in the field of animal behaviour.

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.

Comparative analysis of the gonadal transcriptomes of the all-female species Poecilia formosa and its maternal ancestor Poecilia mexicana

Background

The Amazon molly, Poecilia formosa (Teleostei: Poeciliinae) is an unisexual, all-female species. It evolved through the hybridisation of two closely related sexual species and exhibits clonal reproduction by sperm dependent parthenogenesis (or gynogenesis) where the sperm of a parental species is only used to activate embryogenesis of the apomictic, diploid eggs but does not contribute genetic material to the offspring.

Here we provide and describe the first de novo assembled transcriptome of the Amazon molly in comparison with its maternal ancestor, the Atlantic molly Poecilia mexicana. The transcriptome data were produced through sequencing of single end libraries (100 bp) with the Illumina sequencing technique.

Results

83,504,382 reads for the Amazon molly and 81,625,840 for the Atlantic molly were assembled into 127,283 and 78,961 contigs for the Amazon molly and the Atlantic molly, respectively. 63% resp. 57% of the contigs could be annotated with gene ontology terms after sequence similarity comparisons. Furthermore, we were able to identify genes normally involved in reproduction and especially in meiosis also in the transcriptome dataset of the apomictic reproducing Amazon molly.

Conclusions

We assembled and annotated the transcriptome of a non-model organism, the Amazon molly, without a reference genome (de novo). The obtained dataset is a fundamental resource for future research in functional and expression analysis. Also, the presence of 30 meiosis-specific genes within a species where no meiosis is known to take place is remarkable and raises new questions for future research.

Gradient evolution of body colouration in surface- and cave-dwelling Poecilia mexicana and the role of phenotype-assortative female mate choice

Ecological speciation assumes reproductive isolation to be the product of ecologically based divergent selection. Beside natural selection, sexual selection via phenotype-assortative mating is thought to promote reproductive isolation. Using the neotropical fish Poecilia mexicana from a system that has been described to undergo incipient ecological speciation in adjacent, but ecologically divergent habitats characterized by the presence or absence of toxic H2S and darkness in cave habitats, we demonstrate a gradual change in male body colouration along the gradient of light/darkness, including a reduction of ornaments that are under both inter- and intrasexual selection in surface populations. In dichotomous choice tests using video-animated stimuli, we found surface females to prefer males from their own population over the cave phenotype. However, female cave fish, observed on site via infrared techniques, preferred to associate with surface males rather than size-matched cave males, likely reflecting the female preference for better-nourished (in this case: surface) males. Hence, divergent selection on body colouration indeed translates into phenotype-assortative mating in the surface ecotype, by selecting against potential migrant males. Female cave fish, by contrast, do not have a preference for the resident male phenotype, identifying natural selection against migrants imposed by the cave environment as the major driver of the observed reproductive isolation.

Speciation in fishes

The field of speciation has seen much renewed interest in the past few years, with theoretical and empirical advances that have moved it from a descriptive field to a predictive and testable one. The goal of this review is to provide a general background on research on speciation as it pertains to fishes. Three major components to the question are first discussed: the spatial, ecological and sexual factors that influence speciation mechanisms. We then move to the latest developments in the field of speciation genomics. Affordable and rapidly available, massively parallel sequencing data allow speciation studies to converge into a single comprehensive line of investigation, where the focus has shifted to the search for speciation genes and genomic islands of speciation. We argue that fish present a very diverse array of scenarios, making them an ideal model to study speciation processes.

Genetic differentiation and selection against migrants in evolutionarily replicated extreme environments

We investigated mechanisms of reproductive isolation in livebearing fishes (genus Poecilia) inhabiting sulfidic and nonsulfidic habitats in three replicate river drainages. Although sulfide spring fish convergently evolved divergent phenotypes, it was unclear if mechanisms of reproductive isolation also evolved convergently. Using microsatellites, we found strongly reduced gene flow between adjacent populations from different habitat types, suggesting that local adaptation to sulfidic habitats repeatedly caused the emergence of reproductive isolation. Reciprocal translocation experiments indicate strong selection against immigrants into sulfidic waters, but also variation among drainages in the strength of selection against immigrants into nonsulfidic waters. Mate choice experiments revealed the evolution of assortative mating preferences in females from nonsulfidic but not from sulfidic habitats. The inferred strength of sexual selection against immigrants (RI(s)) was negatively correlated with the strength of natural selection (RI(m)), a pattern that could be attributed to reinforcement, whereby natural selection strengthens behavioral isolation due to reduced hybrid fitness. Overall, reproductive isolation and genetic differentiation appear to be replicated and direct consequences of local adaptation to sulfide spring environments, but the relative contributions of different mechanisms of reproductive isolation vary across these evolutionarily independent replicates, highlighting both convergent and nonconvergent evolutionary trajectories of populations in each drainage.