Fine scale diversity in the lava: genetic and phenotypic diversity in small populations of Arctic charr Salvelinus alpinus

BACKGROUND

A major goal in evolutionary biology is to understand the processes underlying phenotypic variation in nature. Commonly, studies have focused on large interconnected populations or populations found along strong environmental gradients. However, studies on small fragmented populations can give strong insight into evolutionary processes in relation to discrete ecological factors. Evolution in small populations is believed to be dominated by stochastic processes, but recent work shows that small populations can also display adaptive phenotypic variation, through for example plasticity and rapid adaptive evolution. Such evolution takes place even though there are strong signs of historical bottlenecks and genetic drift. Here we studied 24 small populations of the freshwater fish Arctic charr (Salvelinus alpinus) found in groundwater filled lava caves. Those populations were found within a few km2-area with no apparent water connections between them. We studied the relative contribution of neutral versus non-neutral evolutionary processes in shaping phenotypic divergence, by contrasting patterns of phenotypic and neutral genetic divergence across populations in relation to environmental measurements. This allowed us to model the proportion of phenotypic variance explained by the environment, taking in to account the observed neutral genetic structure.

RESULTS

These populations originated from the nearby Lake Mývatn, and showed small population sizes with low genetic diversity. Phenotypic variation was mostly correlated with neutral genetic diversity with only a small environmental effect.

CONCLUSIONS

Phenotypic diversity in these cave populations appears to be largely the product of neutral processes, fitting the classical evolutionary expectations. However, the fact that neutral processes did not explain fully the phenotypic patterns suggests that further studies can increase our understanding on how neutral evolutionary processes can interact with other forces of selection at early stages of divergence. The accessibility of these populations has provided the opportunity for long-term monitoring of individual fish, allowing tracking how the environment can influence phenotypic and genetic divergence for shaping and maintaining diversity in small populations. Such studies are important, especially in freshwater, as habitat alteration is commonly breaking populations into smaller units, which may or may not be viable.

Environmental heterogeneity promotes individual specialisation in habitat selection in a widely distributed seabird

Individual specialisations in behaviour are predicted to arise where divergence benefits fitness. Such specialisations are more likely in heterogeneous environments where there is both greater ecological opportunity and competition-driven frequency dependent selection. Such an effect could explain observed differences in rates of individual specialisation in habitat selection, as it offers individuals an opportunity to select for habitat types that maximise resource gain while minimising competition; however, this mechanism has not been tested before. Here, we use habitat selection functions to quantify individual specialisations while foraging by black-legged kittiwakes Rissa tridactyla, a marine top predator, at 15 colonies around the United Kingdom and Ireland, along a gradient of environmental heterogeneity. We find support for the hypothesis that individual specialisations in habitat selection while foraging are more prevalent in heterogeneous environments. This trend was significant across multiple dynamic habitat variables that change over short time-scales and did not arise through site fidelity, which highlights the importance of environmental processes in facilitating behavioural adaptation by predators. Individual differences may drive evolutionary processes, and therefore these results suggest that there is broad scope for the degree of environmental heterogeneity to determine current and future population, species and community dynamics.

Index of host habitat preference explored by movement-based simulations and trap captures

Animal species likely have different strengths of host habitat preference (HHP) that might be characterized by a standardized index ranging from 0 (no preference) to 1 (maximum preference). We hypothesized that in some species, HHP may result from individuals dispersing out of the host habitat having a probability of turning back at the boundary, or after entering host habitat by reducing speed or increasing size of turning angles. Computer simulations of individuals moving between various sized patches of host and nonhost habitat were conducted based on these three behaviours hypothesized to affect HHP. In the rebounding model, simulations resulted in equilibria of animal numbers inside and outside of host habitat that depend on sizes of these areas, initial number and the rebounding probability. Curvilinear regression of simulation results suggested an equation that predicted numbers in the host habitat and was solved for rebounding probability. A modified equation that sampled population densities (e.g., insect pheromone trap catches) inside and outside host habitat areas gave the rebounding probability, an index of HHP, without requiring the sizes of the areas. Simulations with traps and moving animals verified that the modified equation could predict the index correctly. The modified equation also estimates an index of HHP from sampled densities due to speed reductions and a combination of this and rebounding. Changes in angular turning size upon entering host habitat, however, did not affect habitat preference. Using pheromone trap captures, we found that the lesser date moth Batrachedra amydraula has a HHP for date Phoenix dactylifera plantations of 0.96. Host habitat preference indexes also were calculated from sampled insect densities reported in the literature. The new index of HHP is useful to characterize habitat patches of many organisms and aid understanding of animal spatial distributions and speciation processes. In addition, the index can be applied in studies of invasive species, trap crops of pest insects and conservation management.

The role of habitat choice in micro-evolutionary dynamics: An experimental study on the Mediterranean killifish Aphanius fasciatus (Cyprinodontidae)

Habitat choice is defined as a nonrandom distribution of genotypes in different microhabitats. Therefore, it could exert a great impact on the genetic variance of natural populations by promoting genetic divergence, local adaptation, and may even lead to sympatric speciation. Despite this potential role in micro‐ and macro‐evolutionary processes, there is little empirical evidence that the various genotypes within a population may differ in habitat choice‐related behaviors. Here, we tested whether habitat choice may have contributed to genetic divergence within a local population of the Mediterranean killifish Aphanius fasciatus, which emerged between groups inhabiting microhabitats with different oxygen concentrations during previous field studies. In a first experiment, we studied the distribution of individuals in conditions of hypoxia and normoxia to test whether they had a different ability to shy away from a hypoxic environment; in a second experiment, we analyzed the individual behavior of fish separately in the two conditions, to verify whether they showed peculiar behavioral responses linked to a possible differential distribution. We then analyzed the six allozyme loci, whose allelic and genotypic frequencies were significantly divergent in the previous studies. In the first test, we found that the distribution of the two homozygote genotypes of the glucose‐6‐phosphate isomerase‐1 locus (GPI‐1) was significantly different between the hypoxic and the normoxic conditions. During the second test, all individuals were more active in hypoxic conditions, but the two GPI‐1 homozygotes showed a significant difference in time spent performing surface breathing, which was consistent with their distribution observed in the first experiment. These results provide evidence that individual behavioral traits, related to genetic features, may lead to a nonrandom distribution of genotypes in heterogeneous although contiguous microhabitats and, consequently, that habitat choice can play a significant role in driving the micro‐evolutionary dynamics of this species.

Sexual isolation promotes divergence between parapatric lake and stream stickleback

Speciation can be initiated by adaptive divergence between populations in ecologically different habitats, but how sexually based reproductive barriers contribute to this process is less well understood. We here test for sexual isolation between ecotypes of threespine stickleback fish residing in adjacent lake and stream habitats in the Lake Constance basin, Central Europe. Mating trials exposing females to pairings of territorial lake and stream males in outdoor mesocosms allowing for natural reproductive behaviour reveal that mating occurs preferentially between partners of the same ecotype. Compared to random mating, this sexual barrier reduces gene flow between the ecotypes by some 36%. This relatively modest strength of sexual isolation is surprising because comparing the males between the two ecotypes shows striking differentiation in traits generally considered relevant to reproductive behaviour (body size, breeding coloration, nest size). Analysing size differences among the individuals in the mating trials further indicates that assortative mating is not related to ecotype differences in body size. Overall, we demonstrate that sexually based reproductive isolation promotes divergence in lake-stream stickleback along with other known reproductive barriers, but we also caution against inferring strong sexual isolation from the observation of strong population divergence in sexually relevant traits.

Phenotypic plasticity of mate recognition systems prevents sexual interference between two sympatric leaf beetle species

Maladaptive sexual interactions among heterospecific individuals (sexual interference) can prevent the coexistence of animal species. Thus, the avoidance of sexual interference by divergence of mate recognition systems is crucial for a stable coexistence in sympatry. Mate recognition systems are thought to be under tight genetic control. However, we demonstrate that mate recognition systems of two closely related sympatric leaf beetle species show a high level of host-induced phenotypic plasticity. Mate choice in the mustard leaf beetles, Phaedon cochleariae and P. armoraciae, is mediated by cuticular hydrocarbons (CHCs). Divergent host plant use causes a divergence of CHC phenotypes, whereas similar host use leads to their convergence. Consequently, both species exhibit significant behavioral isolation when they feed on alternative host species, but mate randomly when using a common host. Thus, sexual interference between these syntopic leaf beetles is prevented by host-induced phenotypic plasticity rather than by genotypic divergence of mate recognition systems.

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations - 'outlier loci' (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false-positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated FST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated FST . Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of FST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased FST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated FST in geographically distant locations. We discuss possible explanations for this unexpected result.

Fitness differences between parapatric lake and stream stickleback revealed by a field transplant

Molecular comparisons of populations diverging into ecologically different environments often reveal strong differentiation in localized genomic regions, with the remainder of the genome being weakly differentiated. This pattern of heterogeneous genomic divergence, however, is rarely connected to direct measurements of fitness differences among populations. We here do so by performing a field enclosure experiment in threespine stickleback fish residing in a lake and in three replicate adjoining streams, and displaying weak yet heterogeneous genomic divergence between these habitats. Tracking survival over 29 weeks, we consistently find that lake genotypes transplanted into the streams suffer greatly reduced viability relative to local stream genotypes and that the performance of F1 hybrid genotypes is intermediate. This observed selection against migrants and hybrids combines to a total reduction in gene flow from the lake into streams of around 80%. Overall, our study identifies a strong reproductive barrier between parapatric stickleback populations, and cautions against inferring weak fitness differences between populations exhibiting weak overall genomic differentiation.

How mechanisms of habitat preference evolve and promote divergence with gene flow

Habitat preference may promote adaptive divergence and speciation, yet the conditions under which this is likely are insufficiently explored. We use individual-based simulations to study the evolution and consequence of habitat preference during divergence with gene flow, considering four different underlying genetically based behavioural mechanisms: natal habitat imprinting, phenotype-dependent, competition-dependent and direct genetic habitat preference. We find that the evolution of habitat preference generally requires initially high dispersal, is facilitated by asymmetry in population sizes between habitats, and is hindered by an increasing number of underlying genetic loci. Moreover, the probability of habitat preference to emerge and promote divergence differs greatly among the underlying mechanisms. Natal habitat imprinting evolves most easily and can allow full divergence in parameter ranges where no divergence is possible in the absence of habitat preference. The reason is that imprinting represents a one-allele mechanism of assortative mating linking dispersal behaviour very effectively to local selection. At the other extreme, direct genetic habitat preference, a two-allele mechanism, evolves under restricted conditions only, and even then facilitates divergence weakly. Overall, our results indicate that habitat preference can be a strong reproductive barrier promoting divergence with gene flow, but that this is highly contingent on the underlying preference mechanism.

Do the same genes underlie parallel phenotypic divergence in different Littorina saxatilis populations?

Parallel patterns of adaptive divergence and speciation are cited as powerful evidence for the role of selection driving these processes. However, it is often not clear whether parallel phenotypic divergence is underlain by parallel genetic changes. Here, we asked about the genetic basis of parallel divergence in the marine snail Littorina saxatilis, which has repeatedly evolved coexisting ecotypes adapted to either crab predation or wave action. We sequenced the transcriptome of snails of both ecotypes from three distant geographical locations (Spain, Sweden and United Kingdom) and mapped the reads to the L. saxatilis reference genome. We identified genomic regions potentially under divergent selection between ecotypes within each country, using an outlier approach based on F(ST) values calculated per locus. In line with previous studies indicating that gene reuse is generally common, we expected to find extensive sharing of outlier loci due to recent shared ancestry and gene flow between at least two of the locations in our study system. Contrary to our expectations, we found that most outliers were country specific, suggesting that much of the genetic basis of divergence is not shared among locations. However, we did find that more outliers were shared than expected by chance and that differentiation of shared outliers is often generated by the same SNPs. We discuss two mechanisms potentially explaining the limited amount of sharing we observed. First, a polygenic basis of divergent traits might allow for multiple distinct molecular mechanisms generating the same phenotypic patterns. Second, additional, location-specific axes of selection that we did not focus on in this study may produce distinct patterns of genetic divergence within each site.

Competition-driven speciation in cichlid fish

Theoretically, competition can initiate divergence in habitat use between individuals of a species, leading to restricted gene flow and eventual speciation. Evidence that sister species differ in habitat use is commonplace and consistent with this mechanism, but empirical experimental support is surprisingly scarce. Here we provide evidence that competition has taken a key role in the evolution of genetically distinct ecomorphs of the Lake Tanganyika cichlid fish Telmatochromis temporalis. Experiments show that differences in substrate use between a large-bodied rock-living ecomorph and a neighbouring small-bodied shell-living ecomorph are mediated by size-dependent competition that drives assortative mate-pair formation. Specifically, adults of the larger ecomorph outcompete adults of the smaller ecomorph on favoured rock substrate, compelling the smaller adults to use shell habitat. These results support a role for competition in maintaining reproductive isolation, and highlight the need to identify ecological processes that impose selection to improve our understanding of speciation and adaptive radiation.

Evolutionary branching in complex landscapes

Divergent adaptation to different environments can promote speciation, and it is thus important to consider spatial structure in models of speciation. Earlier theoretical work, however, has been limited to particularly simple types of spatial structure (linear environmental gradients and spatially discrete metapopulations), leaving unaddressed the effects of more realistic patterns of landscape heterogeneity, such as nonlinear gradients and spatially continuous patchiness. To elucidate the consequences of such complex landscapes, we adapt an established spatially explicit individual-based model of evolutionary branching. We show that branching is most probable at intermediate levels of various types of heterogeneity and that different types of heterogeneity have, to some extent, additive effects in promoting branching. In contrast to such additivity, we find a novel refugium effect in which refugia in hostile environments provide opportunities for colonization, thus increasing the probability of branching in patchy landscapes. Effects of patchiness depend on the scale of patches relative to dispersal. Providing a needed connection to empirical research on biodiversity and conservation policy, we introduce empirically accessible spatial environmental metrics that quantitatively predict a landscape's branching propensity.

Evolution of mate choice and the so-called magic traits in ecological speciation

Non-random mating provides multiple evolutionary benefits and can result in speciation. Biological organisms are characterised by a myriad of different traits, many of which can serve as mating cues. We consider multiple mechanisms of non-random mating simultaneously within a unified modelling framework in an attempt to understand better which are more likely to evolve in natural populations going through the process of local adaptation and ecological speciation. We show that certain traits that are under direct natural selection are more likely to be co-opted as mating cues, leading to the appearance of magic traits (i.e. phenotypic traits involved in both local adaptation and mating decisions). Multiple mechanisms of non-random mating can interact so that trait co-evolution enables the evolution of non-random mating mechanisms that would not evolve alone. The presence of magic traits may suggest that ecological selection was acting during the origin of new species.

Looking for a similar partner: host plants shape mating preferences of herbivorous insects by altering their contact pheromones

The role of phenotypical plasticity in ecological speciation and the evolution of sexual isolation remains largely unknown. We investigated whether or not divergent host plant use in an herbivorous insect causes assortative mating by phenotypically altering traits involved in mate recognition. We found that males of the mustard leaf beetle Phaedon cochleariae preferred to mate with females that were reared on the same plant species to females provided with a different plant species, based on divergent cuticular hydrocarbon profiles that serve as contact pheromones. The cuticular hydrocarbon phenotypes of the beetles were host plant specific and changed within 2 weeks after a shift to a novel host plant species. We suggest that plant-induced phenotypic divergence in mate recognition cues may act as an early barrier to gene flow between herbivorous insect populations on alternative host species, preceding genetic divergence and thus, promoting ecological speciation.