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.     

Hybridization of Bombina bombina and B. variegata (Anura, Discoglossidae) at a sharp ecotone in western Ukraine: comparisons across transects and over time

Abstract Bombina bombina and B. variegata are two anciently diverged toad taxa that have adapted to different breeding habitats yet hybridize freely in zones of overlap where their parapatric distributions meet. Here, we report on a joint genetic and ecological analysis of a hybrid zone in the vicinity of Stryi in western Ukraine. We used five unlinked allozyme loci, two nuclear single nucleotide polymorphisms and a mitochondrial DNA haplotype as genetic markers. Parallel allele frequency clines with a sharp central step occur across a sharp ecotone, where transitions in aquatic habitat, elevation, and terrestrial vegetation coincide. The width of the hybrid zone, estimated as the inverse of the maximum gradient in allele frequency, is 2.3 km. This is the smallest of four estimates derived from different clinal transects across Europe. We argue that the narrow cline near Stryi is mainly due to a combination of habitat distribution and habitat preference. Adult toads show a preference for either ponds ( B. bombina ) or puddles ( B. variegata ), which is known to affect the distribution of genotypes within the hybrid zones. At Stryi, it should cause a reduction of the dispersal rate across the ecotone and thus narrow the cline. A detailed comparison of all five intensively studied Bombina transects lends support to the hypothesis that habitat distribution plus habitat preference can jointly affect the structure of hybrid zones and, ultimately, the resulting barriers to gene flow between differentiated gene pools. This study also represents a resampling of an area that was last studied more than 70 years ago. Our allele-frequency clines largely coincide with those that were described then on the basis of morphological variation. However, we found asymmetrical introgression of B. variegata genes into B. bombina territory along the bank of a river.     

The maintenance of reproductive isolation in a mosaic hybrid zone between the fire-bellied toads Bombina bombina and B. variegata

Abstract. Mosaic hybrid zones arise when ecologically differentiated taxa hybridize across a network of habitat patches. Frequent interbreeding across a small-scale patchwork can erode species differences that might have been preserved in a clinal hybrid zone. In particular, the rapid breakdown of neutral divergence sets an upper limit to the time for which differences at marker loci can persist. We present here a case study of a mosaic hybrid zone between the fire-bellied toads Bombina bombina and B. variegata (Anura: Discoglossidae) near Apahida in Romania. In our 20 × 20 km study area, we detected no evidence of a clinal transition but found a strong association between aquatic habitat and mean allele frequencies at four molecular markers. In particular, pure populations of B. bombina in ponds appear to cause massive introgression into the surrounding B. variegata gene pool found in temporary aquatic sites. Nevertheless, the genetic structure of these hybrid populations was remarkably similar to those of a previously studied transect near Pescenica (Croatia), which had both clinal and mosaic features: estimates of heterozygote deficit and linkage disequilibrium in each country are similar. In Apahida, the observed strong linkage disequilibria should stem from an imperfect habitat preference that guides most (but not all) adults into the habitats to which they are adapted. In the absence of a clinal structure, the inferred migration rate between habitats implies that associations between selected loci and neutral markers should break down rapidly. Although plausible selection strengths can maintain differentiation at those loci adapting the toads to either permanent or temporary breeding sites, the divergence at neutral markers must be transient. The hybrid zone may be approaching a state in which the gene pools are homogenized at all but the selected loci, not dissimilar from an early stage of sympatric divergence.     

Genetic and phenotypic variation across a hybrid zone between ecologically divergent tree squirrels (Tamiasciurus)

A hybrid zone along an environmental gradient should contain a clinal pattern of genetic and phenotypic variation. This occurs because divergent selection in the two parental habitats is typically strong enough to overcome the homogenizing effects of gene flow across the environmental transition. We studied hybridization between two parapatric tree squirrels (Tamiasciurus spp.) across a forest gradient over which the two species vary in coloration, cranial morphology and body size. We sampled 397 individuals at 29 locations across a 600‐km transect to seek genetic evidence for hybridization; upon confirming hybridization, we examined levels of genetic admixture in relation to maintenance of phenotypic divergence despite potentially homogenizing gene flow. Applying population assignment analyses to microsatellite data, we found that Tamiasciurus douglasii and T. hudsonicus form two distinct genetic clusters but also hybridize, mostly within transitional forest habitat. Overall, based on this nuclear analysis, 48% of the specimens were characterized as T. douglasii, 9% as hybrids and 43% as T. hudsonicus. Hybrids appeared to be reproductively viable, as evidenced by the presence of later‐generation hybrid genotypes. Observed clines in ecologically important phenotypic traits—fur coloration and cranial morphology—were sharper than the cline of putatively neutral mtDNA, which suggests that divergent selection may maintain phenotypic distinctiveness. The relatively recent divergence of these two species (probably late Pleistocene), apparent lack of prezygotic isolating mechanisms and geographic coincidence of cline centres for both genetic and phenotypic variation suggest that environmental factors play a large role in maintaining the distinctiveness of these two species across the hybrid zone.     

Evidence for nonallopatric speciation among closely related sympatric Heliotropium species in the Atacama Desert

The genetic structure of populations of closely related, sympatric species may hold the signature of the geographical mode of the speciation process. In fully allopatric speciation, it is expected that genetic differentiation between species is homogeneously distributed across the genome. In nonallopatric speciation, the genomes may remain undifferentiated to a large extent. In this article, we analyzed the genetic structure of five sympatric species from the plant genus Heliotropium in the Atacama Desert. We used amplified fragment length polymorphisms (AFLPs) to characterize the genetic structure of these species and evaluate their genetic differentiation as well as the number of loci subject to positive selection using divergence outlier analysis (DOA). The five species form distinguishable groups in the genetic space, with zones of overlap, indicating that they are possibly not completely isolated. Among‐species differentiation accounts for 35% of the total genetic differentiation (F_ST = 0.35), and F_ST between species pairs is positively correlated with phylogenetic distance. DOA suggests that few loci are subject to positive selection, which is in line with a scenario of nonallopatric speciation. These results support the idea that sympatric species of Heliotropium sect. Cochranea are under an ongoing speciation process, characterized by a fluctuation of population ranges in response to pulses of arid and humid periods during Quaternary times.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

Interspecific competition promotes habitat and morphological divergence in a secondary contact zone between two hybridizing songbirds

Interspecific competition is assumed to play an important role in the ecological differentiation of species and speciation. However, empirical evidence for competition's role in speciation remains surprisingly scarce. Here, we studied the role of interspecific competition in the ecological differentiation and speciation of two closely related songbird species, the Common Nightingale (Luscinia megarhynchos) and the Thrush Nightingale (Luscinia luscinia). Both species are insectivorous and ecologically very similar. They hybridize in a secondary contact zone, which is a mosaic of sites where both species co‐occur (syntopy) and sites where only one species is present (allotopy). We analysed fine‐scale habitat data for both species in both syntopic and allotopic sites and looked for associations between habitat use and bill morphology, which have been previously shown to be more divergent in sympatry than in allopatry. We found that the two nightingale species differ in habitat use in allotopic sites, where L. megarhynchos occurred in drier habitats and at slightly higher elevations, but not in syntopic sites. Birds from allotopic sites also showed higher interspecific divergence in relative bill size compared to birds from syntopic sites. Finally, we found an association between bill morphology and elevation. Our results are consistent with the view that interspecific competition in nightingales has resulted in partial habitat segregation in sympatry and that the habitat‐specific food supply has in turn very likely led to bill size divergence. Such ecological divergence may enhance prezygotic as well as extrinsic postzygotic isolation and thus accelerate the completion of the speciation process.     

Genomic analysis of a migratory divide reveals candidate genes for migration and implicates selective sweeps in generating islands of differentiation

Differential gene flow, reductions in diversity following linked selection and/or features of the genome can structure patterns of genomic differentiation during the process of speciation. Possible sources of reproductive isolation are well studied between coastal and inland subspecies groups of Swainson's thrushes, with differences in seasonal migratory behaviour likely playing a key role in reducing hybrid fitness. We assembled and annotated a draft reference genome for this species and generated whole‐genome shotgun sequence data for populations adjacent to the hybrid zone between these groups. We documented substantial genomewide heterogeneity in relative estimates of genetic differentiation between the groups. Within population diversity was lower in areas of high relative differentiation, supporting a role for selective sweeps in generating this pattern. Absolute genetic differentiation was reduced in these areas, further suggesting that recurrent selective sweeps in the ancestral population and/or between divergent populations following secondary contact likely occurred. Relative genetic differentiation was also higher near centromeres and on the Z chromosome, suggesting that features of the genome also contribute to genomewide heterogeneity. Genes linked to migratory traits were concentrated in islands of differentiation, supporting previous suggestions that seasonal migration is under divergent selection between Swainson's thrushes. Differences in migratory behaviour likely play a central role in the speciation of many taxa; we developed the infrastructure here to permit future investigations into the role several candidate genes play in reducing gene flow between not only Swainson's thrushes but other species as well.     

Variation in female mate preference across a grasshopper hybrid zone

A method for studying variation in female preference among populations free from the confounding effects of variation in male characters is described. This method is applied to sixteen populations in a transect across a hybrid zone in the grasshopper, Chorthippus parallelus. Significant variation in preference is revealed among populations. The pattern of this variation is analysed in the context of the reinforcement model of speciation. While the data do indicate an increased homogamic preference on one side of the zone relative to more distant populations, this pattern cannot be distinguished statistically from an abrupt transition in preference over a distance of less than 1 km. We argue that either of these patterns implies that selection operates on female preferences in the hybrid zone.     

Nonallopatric and parallel origin of local reproductive barriers between two snail ecotypes

Theory suggests that speciation is possible without physical isolation of populations (hereafter, nonallopatric speciation), but recent nonallopatric models need the support of irrefutable empirical examples. We collected snails (Littorina saxatilis) from three areas on the NW coast of Spain to investigate the population genetic structure of two ecotypes. Earlier studies suggest that these ecotypes may represent incipient species: a large, thick-shelled 'RB' ecotype living among the barnacles in the upper intertidal zone and a small, thin-shelled 'SU' ecotype living among the mussels in the lower intertidal zone only 10-30 m away. The two ecotypes overlap and hybridize in a midshore zone only 1-3 m wide. Three different types of molecular markers [allozymes, mitochondrial DNA (mtDNA) and microsatellites] consistently indicated partial reproductive isolation between the RB and the SU ecotypes at a particular site. However, each ecotype was related more closely to the other ecotype from the same site than to the same ecotype from another site further along the Galician coast (25-77 km away). These findings supported earlier results based solely on allozyme variation and we could now reject the possibility that selection produced these patterns. The patterns of genetic variation supported a nonallopatric model in which the ecotypes are formed independently at each site by parallel evolution and where the reproductive barriers are a byproduct of divergent selection for body size. We argue that neither our laboratory hybridization experiments nor our molecular data are compatible with a model based on allopatric ecotype formation, secondary overlap and introgression.     

Differential barrier strength and allele frequencies in hybrid zones maintained by sex-biased hybrid incompatibilities

In animals, hybrid sterility and inviability between closely related species often affect only the heterogametic sex (XY). This widespread phenomenon, known as Haldane's rule, is an early speciation event found across broad taxa, but the role of heterogametic hybrid incompatibilities, as opposed to homogametic ones, as a barrier in a speciation process remains obscure. It has been hypothesized that heterogametic incompatibility may be a more efficient mechanism in driving speciation. The population dynamics after (rather than before) the occurrence of sex-biased incompatibilities may account for Haldane's rule. In this study, a recursion model of hybrid zones was developed to investigate the differences between heterogametic and homogametic incompatibilities. The selection strengths and selection patterns of sex chromosome-linked, two-locus Bateson-Dobzhansky-Muller (BDM) genetic incompatibilities were examined. It is noted that a sex-biased hybrid incompatibility in a hybrid zone confers asymmetric and uneven impedance to gene flow. The clines of different loci in such a hybrid zone displayed diverse differentiation in their width, steepness and asymmetry. Alleles involved in the incompatibility face much stronger resistance to cross a hybrid zone. Different sex-biased BDM incompatibilities also affect the flow of neutral alleles differently. Compared to a homogametic one, heterogametic incompatibility is a weaker but more asymmetric barrier. These unique patterns of gene flow may explain uneven divergence among different genomic regions during speciation between some closely related species.     

Ecological differentiation and habitat unsuitability maintaining a ground beetle hybrid zone

Exogenous selection via interactions between organisms and environments may influence the dynamics of hybrid zones between species in multiple ways. Two major models of a hybrid zone allowed us to hypothesize that environmental conditions influence hybrid zone dynamics in two ways. In the first model, an environmental gradient determines the mosaic distribution at the boundary between ecologically differentiated species (mosaic hybrid zone model). In the second model, a patch of unsuitable habitat traps a hybrid zone between species whose hybrids are unfit (tension zone model). To test these, we examined the environmental factors influencing the spatial structure of a hybrid zone between the ground beetles Carabus maiyasanus and C. iwawakianus using GIS‐based quantification of environmental factors and a statistical comparison of species distribution models (SDMs). We determined that both of the hypothetical processes can be important in the hybrid zone. We detected interspecific differences in the environmental factors in presence localities and their relative contribution in SDMs. SDMs were not identical between species even within contact areas, but tended to be similar within the range of each species. These results suggest an association between environments and species, and provide evidence that ecological differentiation between species plays a role in the maintenance of the hybrid zone. Contact areas were characterized by a relatively high temperature, low precipitation, and high topological wetness. Thus, the contact areas were regarded as being located in an unsuitable habitat with a drier climate, where those populations are likely to occur in patches with limited precipitation concentrated. A comparison of spatial scales suggests that exogenous selection via environmental factors may be weaker than endogenous selection via genitalic incompatibility.     

Habitat preferences and positive assortative mating in an avian hybrid zone

The relatively narrow hybrid zone between hooded and carrion crow is characterised by a high proportion of breeding hybrids, no known intrinsic postzygotic costs, and most likely, high gene flow between populations. It is therefore of interest to identify factors preventing a breakdown of the zone. In this study, we investigate the possibility that habitat choice could create assortative mating in the German–Danish part of the zone. In a GIS analysis, we used two approaches to investigate the role of habitat selection in the maintenance of this zone. We tested if there are 1) any correlations between habitat and the shape of the hybrid zone, and 2) any differences in habitat preferences between territorial individuals of the crow phenotypes. The most extensive habitat, non‐irrigated arable land (Habitat 1), which covers almost 50% of the area in the zone was preferred by all three phenotypes (carrion, hooded and hybrid crows), but showed no change in frequency across the zone. Pastures (Habitat 2) and agricultural areas mixed with natural vegetation (Habitat 4) correlated with the shape of the zone, but only H4 differed in preference between the crow phenotypes. H4 was preferred by hybrids and hooded crows, but not by carrion crows. This habitat is present only in 17% of the hooded crow territories, and thus, it is not likely to have a strong influence on the maintenance of the hybrid zone. There was evidence of assortative mating in pure phenotypes, whereas hybrids showed no consistent mating pattern. Because our GIS analyses suggest that species specific habitat preferences are not involved in shaping the assortativeness, we discuss other possible mechanisms.     

Hybrid songbirds employ intermediate routes in a migratory divide

Migratory divides are contact zones between populations that use different routes to navigate around unsuitable areas on seasonal migration. Hybrids in divides have been predicted to employ intermediate and potentially inferior routes. We provide the first direct test of this hypothesis, using light‐level geolocators to track birds breeding in a hybrid zone between Swainson's thrushes in western Canada. Compared to parental forms, hybrids exhibited increased variability in their migratory routes, with some using intermediate routes that crossed arid and mountainous regions, and some using the same routes as one parental group on fall migration and the other on spring migration. Hybrids also tended to use geographically intermediate wintering sites. Analysis of genetic variation across the hybrid zone suggests moderately strong selection against hybrids. These results indicate that seasonal migratory behaviour might be a source of selection against hybrids, supporting a possible role for migration in speciation.     

Genomic regions with a history of divergent selection affect fitness of hybrids between two butterfly species

Speciation is the process by which reproductively isolated lineages arise, and is one of the fundamental means by which the diversity of life increases. Whereas numerous studies have documented an association between ecological divergence and reproductive isolation, relatively little is known about the role of natural selection in genome divergence during the process of speciation. Here, we use genome-wide DNA sequences and Bayesian models to test the hypothesis that loci under divergent selection between two butterfly species (Lycaeides idas and L. melissa) also affect fitness in an admixed population. Locus-specific measures of genetic differentiation between L. idas and L. melissa and genomic introgression in hybrids varied across the genome. The most differentiated genetic regions were characterized by elevated L. idas ancestry in the admixed population, which occurs in L. idas-like habitat, consistent with the hypothesis that local adaptation contributes to speciation. Moreover, locus-specific measures of genetic differentiation (a metric of divergent selection) were positively associated with extreme genomic introgression (a metric of hybrid fitness). Interestingly, concordance of differentiation and introgression was only partial. We discuss multiple, complementary explanations for this partial concordance.     

Habitat selection, acoustic adaptation, and the evolution of reproductive isolation

We examined barriers to gene flow in a hybrid zone of two subspecies of the song sparrow (Melospiza melodia). We focused on how mating signals and mate choice changed along an environmental gradient and gathered data on the morphology, genetics, ecology, and behavior across the zone. Melospiza m. heermanni of the Pacific slope of California and M. m. fallax of the Sonoran Desert, each distinct in plumage, meet across a steep environmental gradient in southeastern California. Although both subspecies occur in riparian habitat, their occupied habitat differs structurally, the former subspecies occurring in areas with denser understory and greater vertical heterogeneity. Song elements varied concomitantly, as predicted by the acoustic adaptation hypothesis, with heermanni having lower-pitched, more widely spaced elements. Females of both subspecies responded more strongly to homotypic than heterotypic song, and addition of subspecific plumage cues increased response if song was homotypic but not if heterotypic. Females thus assess multiple male traits, weighing song more heavily. Males of both subspecies showed significantly greater agonistic response to homotypic song. Microsatellite variation is correlated significantly with plumage variation across the zone and suggests limited gene flow between the taxa. The association of song and plumage with the environment and in turn with assortative mating suggests a means by which reproductive isolation may evolve or be maintained in hybrid zones.     

A fine-scale spatial analysis of the mosaic hybrid zone between Gryllus firmus and Gryllus pennsylvanicus

Abstract.— The pattern of character variation within a hybrid zone, the hybrid zone structure, has been used to infer the processes that maintain hybrid zones. Unfortunately it is difficult to infer process from structure alone because many different processes can produce the same pattern of character variation. Mosaic hybrid zones may be maintained by exogenous selection in a heterogeneous environment and/or endogenous selection against hybrid individuals; habitat preference, premating isolating barriers and/or fertility selection can also contribute. The spatial scale at which a hybrid zone is sampled affects its apparent structure; a hybrid zone may appear clinal at one scale and mosaic at another. Here, we sample the mosaic hybrid zone between two field crickets, Gryllus firmus and G. pennsylvanicus , at a scale that spans the boundaries between individual soil-habitat patches. From our analysis, we find that at fine scales, the mosaic hybrid zone resolves into a set of steep clines across patch boundaries. Both morphological and molecular traits exhibit sharp and generally concordant clines. However, clines for mitochondrial DNA and one anonymous nuclear marker are clearly displaced as a result of current hybridization or past introgression (the "ghost of hybridization past"). Thus, scale is important for the structure of this and probably other hybrid zones. The extremely sharp, concordant clines across patch boundaries indicate that the cricket hybrid zone is undoubtedly structured by selection. However, the detailed mechanisms responsible for the maintenance of the hybrid zone–whether endogenous selection against hybrids, exogenous selection by the environment, and/or behavioral preferences for mates or habitats– remain to be elucidated. Determining these mechanisms will depend on closer inspection of the organisms themselves and their interactions, as is the case for all hybrid zones.     

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.     

Differential gene exchange between parapatric morphs of Littorina saxatilis detected using AFLP markers

Speciation requires the acquisition of reproductive isolation, and the circumstances under which this could evolve are of great interest. Are new species formed after the acquisition of generalized incompatibility arising between physically separated populations, or may they arise as a result of the action of disruptive selection beginning with the divergence of a rather restricted set of gene loci? Here we apply the technique of amplified fragment length polymorphism (AFLP) analysis to an intertidal snail whose populations display a cline in shell shape across vertical gradients on rocky shores. We compare the F_ST values for 306 AFLP loci with the distribution of F_ST estimated from a simulation model using values of mutation and migration derived from the data. We find that about 5% of these loci show greater differentiation than expected, providing evidence of the effects of selection across the cline, either direct or indirect through linkage. This is consistent with expectations from nonallopatric speciation models that propose an initial divergence of a small part of the genome driven by strong disruptive selection while divergence at other loci is prevented by gene flow. However, the pattern could also be the result of differential introgression after secondary contact.     

Food choice behaviour may promote habitat specificity in mixed populations of clonal and sexual Potamopyrgus antipodarum

Genetic polymorphism along an environmental gradient may be maintained if disruptive selection on habitat-specific traits leads to a correlated response in traits that reduce gene flow between habitats. We studied a short-distance cline in a population of freshwater snails Potamopyrgus antipodarum in which sexual and clonal snails coexist. Sexuals and clones show a life history cline by depth: snails reproduce at a smaller size in shallower habitats. Clones are also structured genetically across habitats and seem not to mix, even though habitats are within the dispersal distance of the snails and the opportunity for gene flow via migration must be considerable. Because habitat preference may promote divergence in both clones and sexuals along the depth gradient, we investigated whether snails show habitat-specific food choice behaviour that could reduce migration. We tested the food choice behaviour of the snails by exposing them simultaneously to food from their home and adjacent habitats. Both juvenile and adult snails from the shallow shore bank and a mid-water macrophyte habitat preferentially grazed on the vegetation of their original habitats. We suggest that the observed genetic and life history cline may be maintained by food choice behaviour that may promote a partial barrier to gene flow between the habitats. Copyright 2000 The Association for the Study of Animal Behaviour.