A Simulation Study of the Ecological Speciation Conditions in the Galician Marine Snail Littorina saxatilis.

In the last years, the interest in evolutionary divergence at small spatial scales has increased and so did the study of speciation caused by ecologically based divergent natural selection. The evolutionary interplay between gene flow and local adaptation can lead to low-dispersal locally adapted specialists. When this occurs, the evolutionary interplay between gene flow and local adaptation could eventually lead to speciation. The L. saxatilis system consists of two ecotypes displaying a microhabitat-associated intraspecific dimorphism along the wave-exposed rocky shores of Galicia. Despite being a well-known system, the dynamics of the ecotype formation remain unclear and cannot be studied from empirical evidence alone. In this study, individual-based simulations were used to incorporate relevant ecological, spatial, and genetic information, to check different evolutionary scenarios that could evolve non-random mating preferences and finally may facilitate speciation. As main results, we observed the evolution of intermediate values of choice which matches the estimates from empirical data of L. saxatilis in Galician shores and coincides with previous theoretical outcomes. Also, the use of the mating correlation as a proxy for assortative mating led to spuriously inferring greater reproductive isolation in the middle habitat than in the others, which does not happen when directly considering the choice values from the simulations. We also corroborate the well-known fact that the occurrence of speciation is influenced by the strength of selection. Taken together, this means, also according to other L. saxatilis systems, that speciation is not an immediate consequence of local divergent selection and mating preferences, but a fine tuning among several factors including the ecological conditions in the shore levels, the selection strength, the mate choice stringency, and cost to choosiness. The L. saxatilis system could correspond to a case of incomplete reproductive isolation, where the choice intensity is intermediate and local adaptation within the habitat is strong. These results support previous interpretations of the L. saxatilis model system and indicate that further empirical studies would be interesting to test whether the mate choice mechanism functions as a similarity-like mechanism as has been shown in other littorinids.

Incipient post-zygotic barrier in a model system of ecological speciation with gene flow.

The role of post-zygotic isolation in nonallopatric ecological speciation is still mostly unknown and information on the nature and strength of these barriers in well-known speciation models is essential for a deeper understanding of such processes. The Galician ecotypes of the marine snail Littorina saxatilis represent one of the best studied cases of nonallopatric ecological speciation. Here, we test the existence of incipient post-zygotic isolation by comparing the fertility of male hybrids with that of both pure forms [ridged and banded (RB) and smooth and unbanded (SU) ecotypes]. We analysed the degree of sperm DNA fragmentation (SDF) of individuals morphologically classified as RB, SU and hybrids, sampled from two locations. SDF analyses were chosen to study sperm quality because, in other animal species, SDF rates correlate with important parameters for speciation research, such as fertilization and abortion rates and viability of adult progeny. In the present work, hybrids showed significantly higher SDF rates than RB and SU males in one location and significantly higher variances in both locations. These results suggest the existence of an incipient post-zygotic barrier, the strength of which may vary across the Galician shore, and highlight the potential of SDF analyses for speciation research.

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail.

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged-Banded (RB) and Smooth-Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid-shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome-wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.

Proteomic evidence of a paedomorphic evolutionary process within a marine snail species: a strategy for adapting to extreme ecological conditions?

The exposed and sheltered ecotypes of the marine snail Littorina saxatilis from European rocky shores are considered a key model system to study adaptation and ecological speciation. Previous studies showed that two ecotypes (RB and SU) of this species in NW Spain have adapted differently to different shore levels and microhabitats. In order to understand how this divergent adaptive process has been accomplished, we followed a quantitative proteomic approach to investigate the proteome variation in a number of different biological factors, that is, ecotype, ontogeny and their interactions. This approach allowed testing the hypothesis that one of the ecotypes has evolved by paedomorphosis, and also whether or not the molecular mechanisms related to ecotype differentiation are set up in early developmental stages. Additionally, the identification of some candidate proteins using mass spectrometry provides some functional insights into these evolutionary processes. Results from this study provided evidence of higher ontogenetic differentiation at proteome level in the RB (metamorphic) than in SU (paedomorphic) ecotype that point to the possibility of juvenile stage retention in this latter ecotype. The level of protein expression (proteome) differences between ecotypes maintained nearly constant from late embryonic stages to adulthood, although some proteins involved in these changes considerably differed in embryonic compared to other ontogenetic stages. Paedomorphosis may be the evolutionary response of the SU ecotype of solving the trade-off during sexually immaturity that is caused by the evolution of small size arising from adaptation to the wave-exposed habitat. Some potential candidate genes of adaptation related to energetic metabolism have been identified, providing a promising baseline for future functional analyses.

Comparing three different methods to detect selective loci using dominant markers.

We carried out a simulation study to compare the efficiency of three alternative programs (DFDIST, DETSELD and BAYESCAN) to detect loci under directional selection from genome-wide scans using dominant markers. We also evaluated the efficiency of correcting for multiple testing those methods that use a classical probability approach. Under a wide range of scenarios, we conclude that BAYESCAN appears to be more efficient than the other methods, detecting a usually high percentage of true selective loci as well as less than 1% of outliers (false positives) under a fully neutral model. In addition, the percentage of outliers detected by this software is always correlated with the true percentage of selective loci in the genome. Our results show, nevertheless, that false positives are common even with a combination of methods and multitest correction, suggesting that conclusions obtained from this approach should be taken with extreme caution.

Semi-quantitative differences in gene transcription profiles between sexes of a marine snail by a new variant of cDNA-AFLP analysis.

A variant of the cDNA-AFLP method coupled to an automated sequencer was used to quantify transcripts differentially expressed between sexes of the marine snail Littorina saxatilis. First, we conducted a validation study of the technique using known concentrations of a commercial marker. Second, we analysed six replicates of males and females from a population showing no apparent sexual dimorphism. The results confirm that the method can be properly used within the range of DNA concentrations utilized. In addition, we detected a small percentage of spots (1.8%) differentially expressed between sexes, as expected from a low to moderately sexual dimorphic species.

Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis.

The Galician sympatric ecotypes of Littorina saxatilis have been proposed as a model system for studying parallel ecological speciation. Such a model system makes a clear prediction: candidate loci (for divergent adaptation) should present a higher level of geographical differentiation than noncandidate (neutral) loci. We used 2356 amplified fragment length polymorphisms (AFLPs) and four microsatellite loci to identify candidate loci for ecological adaptation using the F(ST) outlier method. Three per cent of the studied AFLP loci were identified as candidate loci associated with adaptation, after multitest adjustments, thus contributing to ecotype differentiation (candidate loci were not detected within ecotypes). Candidate and noncandidate loci were analysed separately at four different F(ST) partitions: differences between ecotypes (overall and local), differences between localities and micro-geographical differences within ecotypes. The magnitude of F(ST) differed between candidate and noncandidate loci for all partitions except in the case of micro-geographical differentiation within ecotypes, and the microsatellites (putatively neutral) showed an identical pattern to noncandidate loci. Thus, variation in candidate loci is determined partially independent by divergent natural selection (in addition to stochastic forces) at each locality, while noncandidate loci are exclusively driven by stochastic forces. These results support the evolutionary history described for these particular populations, considered to be a clear example of incomplete sympatric ecological speciation.

Disassociation between weak sexual isolation and genetic divergence in a hermaphroditic land snail and implications about chirality.

Examination of the association between reproductive isolation and genetic divergence in a variety of organisms is essential for elucidating the mechanisms causing speciation. However, such studies are lacking for hermaphrodites. We measured premating (sexual) isolation in species pairs of the hermaphroditic land snail Albinaria and we compared it with their genetic divergence. We did not find substantial sexual isolation barriers between the species studied. The absence of strong sexual isolation between species implies its minor effect in the evolution of this genus, because distributional, population and life-history characteristics of Albinaria make mate-choice possibly redundant. Furthermore, we found disassociation between genetic divergence and sexual isolation, suggesting that they do not form necessarily a cause-effect duet. However, Albinaria voithii, the only dextral Albinaria species, shows strong sexual isolation against the other sinistral species. We discuss whether change in coiling either has triggered instantaneous speciation, or is an example of character displacement.

Genetic differentiation and estimation of effective population size and migration rates in two sympatric ecotypes of the marine snail Littorina saxatilis.

On exposed rocky shores in Galicia (northwest Spain), a striking polymorphism exists between two ecotypes (RB and SU) of Littorina saxatilis that occupy different levels of the intertidal zone and exhibit an incomplete reproductive isolation. The setting has been suggested to represent ongoing sympatric speciation by ecological adaptation of the two ecotypes to their respective habitats. In this article we address whether or not the ecotypes have developed their own population structures in response to the rigors of their corresponding environments and life histories. We analyzed four to five allozymic loci from three surveys of the same sites, spanning a 14-year period. An experimental design including three localities with two transects per locality and three shore levels allowed studying temporal and spatial population structure and estimation of effective population sizes (N(e)), neighborhood sizes (N(n)), and migration rates (m). Genetic differentiation was significantly lower in RB populations (theta(ST) = 0.067) than in SU ones (theta(ST) = 0.124). Mean estimates of N(e), N(n), and m did not differ significantly between ecotypes, but local ecotype differences in migration between the two closest localities (larger migration rates in RB than in SU populations) could explain the pattern in population differentiation.

The evolutionary forces maintaining a wild polymorphism of Littorina saxatilis: model selection by computer simulations.

Two rocky shore ecotypes of Littorina saxatilis from north-west Spain live at different shore levels and habitats and have developed an incomplete reproductive isolation through size assortative mating. The system is regarded as an example of sympatric ecological speciation. Several experiments have indicated that different evolutionary forces (migration, assortative mating and habitat-dependent selection) play a role in maintaining the polymorphism. However, an assessment of the combined contributions of these forces supporting the observed pattern in the wild is absent. A model selection procedure using computer simulations was used to investigate the contribution of the different evolutionary forces towards the maintenance of the polymorphism. The agreement between alternative models and experimental estimates for a number of parameters was quantified by a least square method. The results of the analysis show that the fittest evolutionary model for the observed polymorphism is characterized by a high gene flow, intermediate-high reproductive isolation between ecotypes, and a moderate to strong selection against the nonresident ecotypes on each shore level. In addition, a substantial number of additive loci contributing to the selected trait and a narrow hybrid definition with respect to the phenotype are scenarios that better explain the polymorphism, whereas the ecotype fitnesses at the mid-shore, the level of phenotypic plasticity, and environmental effects are not key parameters.

Testing alternative models for sexual isolation in natural populations of Littorina saxatilis: indirect support for by-product ecological speciation?

Two ecotypes of the rough periwinkle Littorina saxatilis occur at different shore levels, showing assortative mating for size and partial reproductive isolation when they meet at the mid-shore. This system represents a putative case of incomplete speciation in sympatry. Two processes contribute to the assortative mating: morph-specific microhabitat aggregation and mate choice. The estimation of mate choice coefficients in nature and a simulation of the aggregation effects on sexual isolation were used to disentangle these processes as well as to test alternative mechanisms of mate choice. Mate choice significantly increased the frequency of within-morph pairs and significantly decreased the frequency of between-morph pairs, whereas those pairs including at least one hybrid morph mated randomly. These results allow us to reject a discriminant mate choice and support a model of evolution of sexual isolation as a side-effect of size-assortative mating in a context of divergent natural selection for size in the population. This mechanism is more compatible with a model of incomplete by-product ecological speciation, as suggested by previous evidence.

Estimating sexual selection and sexual isolation effects from mating frequencies.

Sexual selection (defined as the change in genotypic or phenotypic frequencies of mated versus total population frequencies) and sexual isolation (defined as the deviation from random mating in mated individuals) show different evolutionary consequences and partially confounded causes. Traditionally, the cross-product estimator has been used to quantify sexual selection, whereas a variety of indexes, such as Yule V, Yule Q, YA, joint I, and others have been used to quantify sexual isolation. Because the two types of estimators use different scales, the effects of both processes cannot be monitored simultaneously. We describe three new related statistics that quantify both sexual selection (PSS) and sexual isolation (PSI) effects for every mating pair combination in polymorphic traits, as well as measure their combined effects (PTI = PSI x PSS). The new statistics have the advantage of providing information on every mating pair combination, quantifying the effects of sexual selection and isolation in the same units, and detecting asymmetry in sexual isolation. The ability of the new statistics to ascertain the biological causes of sexual selection and sexual isolation are investigated under different models involving distinct marginal frequencies, mate propensity, and mate choice coefficients. We also studied the use of classical isolation indexes applied on PSI coefficients, instead of on raw data. The use of the classical indexes applied to PSI coefficients considerably reduces the statistical bias of the estimates, revealing the good estimation properties of the new statistics.

Distinct genetic subdivision in sympatric and sibling species of the genus Littorina (Gastropoda: Littorinidae).

The genetic structure of two sibling and sympatric species of the genus Littorina was compared using allozymic loci. The two species are biologically and ecologically well-known and mostly show similar life history characteristics. Three populations of L. mariae Sacchi & Rastelli and L. obstusata (L.) were studied in the Muros-Noya Ria (Galicia, NW Spain). In addition, four microgeographical subsamples taken from one of the populations were analysed for each species. Age, sex and genotypes for nine polymorphic loci were studied in 1250 snails of both species. L. mariae showed larger genetic population subdivision and lower heterozygosity levels for the loci studied than did L. obtusata. Heterozygote deficiencies were found in only a few cases in natural populations of both species, usually affecting the Lap-1 locus. No significant genetic differences among age or sex classes were found. These results may be explained by the lower effective population size in L. mariae than in L. obtusata. Known differences between these species in generation interval and population density during the winter can cause the different effective population sizes suggested. These life history characteristics appear to provide the most likely explanations for the differences in genetic differentiation and heterozygosity between the two species. A previously unknown L. mariae morph from exposed shores is tentatively suggested to be conspecific.