Parallel evolution of similar phenotypes provides strong evidence for the operation of natural selection. Where these phenotypes contribute to reproductive isolation, they further support a role for divergent, habitat-associated selection in speciation. However, the observation of pairs of divergent ecotypes currently occupying contrasting habitats in distinct geographical regions is not sufficient to infer parallel origins. Here we show striking parallel phenotypic divergence between populations of the rocky-shore gastropod, Littorina saxatilis, occupying contrasting habitats exposed to either wave action or crab predation. This divergence is associated with barriers to gene exchange but, nevertheless, genetic variation is more strongly structured by geography than by ecotype. Using approximate Bayesian analysis of sequence data and amplified fragment length polymorphism markers, we show that the ecotypes are likely to have arisen in the face of continuous gene flow and that the demographic separation of ecotypes has occurred in parallel at both regional and local scales. Parameter estimates suggest a long delay between colonization of a locality and ecotype formation, perhaps because the postglacial spread of crab populations was slower than the spread of snails. Adaptive differentiation may not be fully genetically independent despite being demographically parallel. These results provide new insight into a major model of ecologically driven speciation.

Distinct ecotypes of the snail Littorina saxatilis, each linked to a specific shore microhabitat, form a mosaic-like pattern with narrow hybrid zones in between, over which gene flow is 10 -30% of within-ecotype gene flow. Multi-locus comparisons cluster populations by geographic affinity independent of ecotype, while loci under selection group populations by ecotype. The repeated occurrence of partially reproductively isolated ecotypes and the conflicting patterns in neutral and selected genes can either be explained by separation in allopatry followed by secondary overlap and extensive introgression that homogenizes neutral differences evolved under allopatry, or by repeated evolution in parapatry, or in sympatry, with the same ecotypes appearing in each local site. Data from Spain, the UK and Sweden give stronger support for a non-allopatric model of ecotype formation than for an allopatric model. Several different non-allopatric mechanisms can, however, explain the repeated evolution of the ecotypes: (i) parallel evolution by new mutations in different populations; (ii) evolution from standing genetic variation; and (iii) evolution in concert with rapid spread of new positive mutations among populations inhabiting similar environments. These models make different predictions that can be tested using comprehensive phylogenetic information combined with candidate loci sequencing.

The marine gastropod Littorina saxatilis has different ecotypes in shores only a few meters apart. This has both taxonomic and evolutionary implications. Here we report on an extreme type of within-shore dimorphism in shell characters. In the wave-exposed rocky shores in northwestern Spain, we found one form of L. saxatilis in the upper-level barnacle zone. It had a white, ridged shell, with black bands in the grooves. Another form confined to the lower-shore mussel belt had a smooth shell that was either white and tessellated or darkly colored. These two forms cooccured in a narrow midshore zone together with individuals that had combined characters, but were present in low frequencies (II %-29%). We used principal-component analysis of metric shell characters to study variation in shell size and shape. We found that the upper-shore form was larger than the lower-shore form. We also found small but significant differences in shell shape. Experiments in a common laboratory environment suggested the differences in shell ornamentation and color are inherited, but the individuals did not develop the morph-specific characters until a shell height of about 3 mm. The occurrence of mainly two distinct forms may suggest the presence of two species that hybridize. An analysis of five polymorphic enzyme loci in populations of snails from three geographically separated sites indicated, however, that there was no positive correlation between morphological distances and genetic distances among populations on a geographic scale (tens ofkilometers). Thus, we rejected the hypothesis oftwo species. However, on a microgeographic scale (meters), genetic differentiation between groups with the same form was less than differentiation between forms. This indicated a partial barrier to gene flow between the two forms, and preliminary mate choice data suggested this was caused by nonrandom mating in the midshore zone of overlap.

Abstract. 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 ϫ 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.

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.

No abstract

The study of speciation in recent populations is essentially a study of the evolution of reproductive isolation mechanisms between sub-groups of a species. Prezygotic isolation can be of central importance to models of speciation, either being a consequence of reinforcement of assortative mating in hybrid zones, or a pleiotropic effect of morphological or behavioral adaptation to different environments. To suggest speciation by reinforcement between incipient species one must at least know that gene flow occurs, or have recently occurred, and that assortative mating has been established in the hybrid zone. In Galician populations of the marine snail Littorina saxatilis, two main morphs appear on the same shores, one on the upper-shore barnacle belt and the other in the lower-shore mussel belt. The two morphs overlap in distribution in the midshore where hybrids are found together with pure forms. Allozyme variation indicates that the two parental morphs share a common gene pool, although within shores, gene flow between morphs is less than gene flow within morphs. In this study, we observed mating behavior in the field, and we found that mating was not random in midshore sites, with a deficiency of heterotypic pairs. Habitat selection, assortative mating, and possibly sexual selection among females contributed to the partial reproductive isolation between the pure morphs. Sizes of mates were often positively correlated, in particular, in the upper shore, indicating size-assortative mating too. However, this seemed to be a consequence of nonrandom microdistributions of snails of different sizes. Because we also argue that the hybrid zone is of primary rather than secondary origin, this seems to be an example of sympatric reproductive isolation, either established by means of reinforcement or as a by-product to divergent selection acting on other characters.

AFLP markers are becoming one of the most popular tools for genetic analysis in the fields of evolutionary genetics and ecology and conservation of genetic resources. The technique combines a high-information content and fidelity with the possibility of carrying out genomewide scans. However, a potential problem with this technique is the lack of homology of bands with the same electrophoretic mobility, what is known as fragment-size homoplasy. We carried out a theoretical analysis aimed at quantifying the impact of AFLP homoplasy on the estimation of within-and between-neutral population genetic diversity in a model of a structured finite population with migration among subpopulations. We also investigated the performance of a currently used method (DFDIST software) to detect selective loci from the comparison between genetic differentiation and heterozygosis of dominant molecular markers, as well as the impact of AFLP homoplasy on its effectiveness. The results indicate that the biases produced by homoplasy are: (1) an overestimation of the frequency of the allele determining the presence of the band, (2) an underestimation of the degree of differentiation between subpopulations, and (3) an overestimation or underestimation of the heterozygosis, depending on the allele frequency of the markers. The impact of homoplasy is quickly diminished by reducing the number of fragments analyzed per primer combination. However, substantial biases on the expected heterozygosity (up to 15-25%) may occur with $50-100 fragments per primer combination. The performance of the DFDIST software to detect selective loci from dominant markers is highly dependent on the number of selective loci in the genome and their average effects, the estimate of genetic differentiation chosen to be used in the analysis, and the critical bound probability used to detect outliers. Overall, the results indicate that the software should be used with caution. AFLP homoplasy can produce a reduction of up to 15% in the power to detect selective loci. et al. 1995) is becoming one of the most popular methods in the fields of conservation and evolutionary genetics and ecology (Mueller and Wolfenbarger 1999;Bensch and Akesson 2005;Bonin et al. 2007;Meudt and Clarke 2007), as it combines a high reproducibility and information content with the possibility of making genomewide screenings. Because of the anonymous nature of the fragments generated by the AFLP technique, however, one major concern is the incidence of size homoplasy due to the lack of homology of comigrating fragments. This implies that fragments of a given size migrating in a band may involve more than one locus of the genome and, therefore, the inferences obtained from the band can produce misleading conclusions. T HE amplified fragment length polymorphism (AFLP) technique (VosSeveral empirical approaches have been used to estimate levels of homoplasy in AFLP data sets. A few studies have demonstrated the presence of homoplasy in AFLP data by sequencing comigrating fragments within t...