Micro- and macrogeographic allozyme variation in Littorina fabalis; do sheltered and exposed forms hybridize?

Earlier studies of Swedish populations of the marine snail Littorina fabalis show that snails from different microhabitats (with a greater and lesser exposure to wave action) have almost diagnostic differences in one allozyme locus (arginine kinase, Ark ), and differ in adult size. Snails with 'sheltered' and 'exposed' Ark genotypes occur in sympatry in intermediary exposed sites and here adult sizes remain distinct. Approaching the microgeographic differentiation we studied the parts of two populations where the frequency of Ark changes dramatically over zones 50–120 m wide. The aim was to test if the transitional zones are best described as areas of mixing of two genetically separate populations, or if hybridization between the exposed and sheltered groups occurs. Heterozygotes were in deficiency along both clines but were still roughly twice as common as expected from a pure mixing of 'sheltered' and 'exposed' groups suggesting hybridization. Hybridization was also supported by the observation that snails homozygous for sheltered and exposed alleles mated at random with each other in both populations. On the macrogeographic scale, we found populations from exposed and sheltered sites in France and Wales being fixed for the same exposed and sheltered Ark alleles as found in Sweden. However, variation in three other highly polymorphic loci indicated geographic affinity rather than habitat similarity being the main factor of genetic coherence. These observations support a hypothesis of gene flow between exposed and sheltered populations of L. fabalis . Two Spanish populations were remarkably different with unique alleles at high frequencies in three of four strongly polymorphic loci.     

Habitat-related genetic substructuring in a marine snail (Littorina fabalis) involving a tight link between an allozyme and a DNA locus

Metapopulation structure and genetic differentiation among subpopulations will be tightly related to patterns and processes of local adaptation and microevolution. Understanding the mechanisms behind genetic substructuring will aid in the interpretation of species' ecological performances and strategies. The marine gastropod Littorina fabalis occurs in two size morphs – a small and a large – found in microhabitats of different wave exposure, but overlapping in distribution where wave exposure is intermediate. Earlier studies have found substantial genetic differentiation linked to morph in one allozyme locus (arginine kinase), while 29 other allozyme loci reveal no or minute differences between morphs. Here we add new results showing DNA variation in a RAPD marker being tightly linked to the allozyme variation. Indeed, 97% of the snails homozygotic for one of the Ark alleles had a unique DNA band, while 89% of the snails homozygotic for the other Ark allele lacked the marker. We discuss alternative hypotheses explaining the genetic substructure and suggest that the linkage of size, allozyme and DNA traits might be due to a paracentric chromosomal inversion involving loci coding for these traits. A genetic linkage of traits might promote microhabitat specialization of this species, and such a chromosomal transformation may therefore be adaptive.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 81 , 301–306.     

Habitat related allozyme variation on a microgeographic scale in the marine snail Littorina mariae (Prosobranchia: Littorinacea)

The marine snail Littorina mariae Sacchi & Rastelli occurs in high numbers in the littoral zone on fucoid macro-algae. The eggs are laid on the seaweed and d e velopment to miniature snails takes place without any pelagic larval stage. We have mapped the genetic variation of 30 enzyme loci in populations from eight small islands within 15 km of each other on th swedish west coast. The original intention was to investigate the magnitude of gene flow within and between islands. However, we soon realized that our basic assumption of neutral genetic variation was forcefully violated in at least one locus, arginine kinase ( Ask ). Allele frequencies of Ark were strongly associated with type of habitat. Therefore the main part of this study focused on allele frequency distribution in different habitats. Eight of the 30 loci screened were polymmorphic but we mainly considered the four most polymorphic ones (total heterozygosity between 0.32 and 0.57). All four showed significant heterogeneity between subpopulations, Ask especially so; 42% of the total variation of Ark was explained by differentiation between samples, and 90% of this variation was attributed to differences between different types of habitats (more or less exposed to wave action). In contrast, peptidase ( Pep-1 ) and phosphoglucomutase ( Pgm-2 ) varied in a way predicted by neutral theory; between sample variation being mainly attributed to differentiation bvetween islands. The variation in phosphoglucose isomerase ( Pgi ) was less consistent. In some islands there w as an obvious difference between different habitats, but on other islands we found no significant difference.Not taking into account the extreme, presumably sclected, variation in Ark , we concluded that the metapopulation of Littoriana mariae that we studied was divided into semi-isolated populations, between which the average rate of migration was in the range of a fe w individuals per generation.     

Refuge function of marine algae complicates selection in an intertidal snail

Species with restricted gene flow often show trait-shifts from one type of environment to another. In those rock-dwelling marine gastropods that lack larval dispersal, size generally decreases in wave-exposed habitats reducing risk of dislodgement, while increases in less exposed habitats to resist crab-crushing. In Littorina fabalis, however, snails of moderately exposed shores are generally much larger (11–14 mm) than snails of sheltered shores (5–8 mm). Observations from the White Sea (where crabs are not present) indicate that in the absence of crabs snails are small (6–7 mm) in both habitats. We assumed that the optimal size for L. fabalis in the absence of crabs is less than 8 mm, and thus that increased size in moderately exposed habitats in areas with crabs might be a response to crab predation. In a crab-rich area (Sweden) we showed that crab predation is an important mortality factor for this snail species in both sheltered and moderately exposed habitats. In sheltered habitats, snails were relatively more protected from crab-predation when dwelling on their habitual substrate, fucoid algae, than if experimentally tethered to rocks below the algae. This showed that algae function as snail refuges. Snail dislodgement increased, however, with wave exposure but tethering snails in moderately exposed habitats showed that large snails survived equally well on rocks under the algae as in the canopy of the algae. Thus in sheltered habitats a small snail size is favored, probably due to life-history reasons, while increased risk of being dislodged from the algae refuges promotes a large size in moderately exposed habitats. This study shows an example of selection of a trait depends on complex interactions of different factors (life-history optimization, crab predation, wave induced dislodgement and algal refuges).     

Distribution and biomass of microphytes measured by benthic chlorophyll a in a tropical lagoon (New Caledonia,South Pacific)

Starch-gel electrophoretic techniques were applied to the investigation of molecular genetic variation in populations of the rough periwinkle Littorina saxatilis. The investigation comprised two phases: a) technique development to resolve as many as possible of the allozyme loci reported in the literature as having been screened in the genus Littorina and in Melarhaphe neritoides; b) the use of these loci to assess levels of genetic variation in and patterns of genetic differentiation among populations of L. saxatilis from a relatively isolated group of populations from Galway Bay, Ireland. More than 43 allozyme loci (of which four were screened for the first time here in this species), coding for 37 enzymes, were investigated and thirteen of these loci (including two loci screened for the first time here) were found to be variable and reliably scorable. Samples from five pairs of transects were collected from Inismór, Aran Islands, from sites with known exposure levels; one transect within each pair was collected from an exposed site and the other from a nearby, but relatively sheltered site. UPGMA for eleven loci, (ARK and PGDH were excluded from cluster and FST analysis as they were unscorable in a few samples), showed that the samples cluster mostly by pair, reflecting their geographic origin and is indicative of little gene flow between populations. Levels of population differentiation were high among samples from the top of the shore, but unusually so at AAT-1 which showed nearly three times the mean FST value for the eleven loci. There was also a significant regression of frequency of AAT-1100 against level of exposure. In addition, among midshore samples, there was a consistently higher frequency of AAT-1100 in sheltered habitats. These results support the findings of others, indicating that this locus may be subject to natural selection.     

Microgeographic variation in allozyme and shell characters in Littorina saxatilis Olivi (Prosobranchia: Littorinidae)

Twenty-three enzymes and five shell parameters were screened in 11 subpopulations of Littorina saxatilis Olivi (= L. rudis Maton) occupying different habitats over a 1 km stretch of coastline. Shell morphology varied considerably and consistently with respect to degree of exposure, and since there is evidence that such morphology is at least partly under genetic control, it is likely that natural selection selects particular genotypes at particular locations. There was significant allozyme heterogeneity between neighbouring subpopulations, sometimes only metres apart, but little of the allozyme variability could be related directly to environmental pressures. Thus, with the exception of the Odh locus, the considerable morphological differentiation between snails from exposed and sheltered sites was not reflected in differentiation of those genes coding for electrophoretically assayed enzymes. At the Odh locus, virtually all the genetic differentiation between subpopulations was attributable to differentiation between habitat types. Two loci, Sod-1 and Aat-1, showed highly significant genetic disequilibrium, and possible reasons for this are explored. The population structure as assessed electrophoretically accords well with the stepping-stone model which permits greater differentiation of neighbouring populations than the island model, and which seems realistic in the ovoviviparously reproducing L. saxatilis, where the greater part of gene flow is likely to occur through the occasional migration of adults between contiguous populations.     

Genetic heterogeneity among intertidal habitats in the flat periwinkle, Littorina obtusata

Comparisons among patterns exhibited by functionally distinct genetic markers have been widely used to infer the impacts of demography and selection in structuring genetic variation in natural populations. However, such multilocus comparisons remain an indirect evaluation of selection at particular candidate loci; ideally, the identification of a candidate gene by comparative genetic methodologies should be complemented by functional analyses and experimental manipulations of genotypes in the laboratory or field. We examined genotype frequency variation among replicated intertidal habitats at two spatial scales in the grazing snail Littorina obtusata. Both of the candidate allozyme markers varied predictably with environment, and these patterns were consistent at both spatial scales. Three of four reference loci were spatially homogeneous, but one microsatellite exhibited significant structure at both geographical and mesoscales. To initiate a direct examination of whether the observed genotype frequency variation at one of the candidate markers, mannose-6-phosphate isomerase (MPI), was impacted by differential survivorship of genotypes, we conducted a series of laboratory-based thermal stress assays using snails from two geographically disparate source populations. When snails were exposed to bouts of thermal/desiccation stress, patterns of mortality were nonrandom with respect to MPI genotype. Furthermore, patterns of mortality in the laboratory manipulation coincided with the observed distribution of genotypes in the field. The data suggest the operation of selection at the Mpi or a linked locus, but functional studies and further experimentation are required to establish the relationship between MPI genotype and fitness across heterogeneous intertidal environments.     

Character displacement, frequency-dependent selection, and divergence of shell colour in land snails Mandarina (Pulmonata)

Endemic land snails of the genus Mandarina of the oceanic Bonin Islands offer an exceptional example of habitat and character divergence among closely related species. In this study, microhabitat differences between sympatric ground-dwelling species were studied by distinguishing habitats on the basis of vegetation and types of litter. In all sites where two ground species coexisted, segregation occurred with each species showing preference for the microhabitat in which they were found. When they were in sympatry, one species was predominant in relatively wet and sheltered sites and the other in relatively dry and exposed sites. Although most species can live in both types of habitat, occupation by one species is inhibited by occupation by another. This suggests that competitive interaction between sympatric species caused segregation. Except for populations that have undergone interspecific hybridization, no examples were found of sympatric populations of two ground species sharing a similar shell colour. Species that were predominant in relatively wet and sheltered sites possessed shells with dark coloration and their colour patterns were mostly of one type. Species that were predominant in relatively dry and exposed sites possessed shells with bright coloration and their color patterns were polymorphic. Most populations from areas in which single species were distributed had shells with medium coloration. Microhabitat differentiation between sympatric species possibly caused diversification of shell colour, because bright shells are advantageous in sites where snails are largely exposed, and dark shells are advantageous in sites in where they are mostly sheltered from sunlight. In addition, frequency-dependent selection by predators hunting by sight may have operated to maintain colour polymorphism in the populations which are restricted to exposed habitats by competition with other sympatric species. This reveals the importance of interaction among closely related species as a cause of diversification in ecological and morphological traits.     

Temporal variation in parasite loading in relation to life history patterns of Littorina obtusata and L. fabalis

Parasitic infection of Littorina obtusata (L.) and Littorina fabalis Turton by larval digeneans was studied in summer (August) and winter (February) on a sheltered shore at St Michael's Island, Isle of Man. Sampling was timed to investigate parasite loading at different stages of the littorinids' life history. Nine species of digeneans were identified and infection rates were similar between the species, with a slightly higher prevalence in females for L. obtusata and males in L. fabalis. Temporal differences in infection were small, although metacercariae and degenerate infections were greater in female L. obtusata in February. Analysis of shell parameters and size of the penis in males and the pallial oviduct in females revealed morphological differences between parasitised and uninfected individuals. Parasitic stunting of the penis was evident in both species and the pallial oviduct was also smaller in parasitised females. There was also a tendency for parasitised individuals to have larger shells; although size also varied with season for L. fabalis.     

Island isolation and habitat heterogeneity correlate with DNA variation in a marine snail (Littorina saxatilis)

The null assumption of molecular variation is that most of it is neutral to natural selection. This is in contrast to variation in morphological traits that we generally assume is maintained by selection, and therefore often by selection coupled to environmental heterogeneity in time and space. Examples of molecular variation that vary over habitat-shifts, particularly in allozymes, show that the relative impact of non-neutral variation as compared to neutral variation might be substantial in some systems. To assess the importance of habitat-generated variation in relation to variation generated by random processes in nuclear DNA markers at small spatial scales, we compared the effects of island isolation and habitat heterogeneity on genetic substructuring in a rocky shore snail ( Littorina saxatilis ). This species has a restricted migration among islands owing to the lack of free-floating larvae. Earlier studies show that allozymes vary extensively as a consequence of isolation by water barriers among islands, but also as a consequence of divergent selection among different microhabitats within islands. In the DNA markers we observed genetic differentiation owing to island isolation at three of nine loci. In addition, variation at three loci correlated with habitat type, but the correlation for two of the loci was weak. Overall, isolation contributed slightly more to the genetic variation among populations than did habitat-related factors but the difference was small. It is concluded that both island isolation, which interrupts gene flow, and a heterogeneous habitat cause genetic substructuring at the DNA level in L. saxatilis in the studied area, and thus in this species we need to be somewhat concerned about habitat heterogeneity also at DNA loci.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 377–384.     

Genetic differentiation on multiple spatial scales in an ecotype-forming marine snail with limited dispersal: Littorina saxatilis

The population genetic structure of marine species lacking free-swimming larvae is expected to be strongly affected by random genetic drift among populations, resulting in genetic isolation by geographical distance. At the same time, ecological separation over microhabitats followed by direct selection on those parts of the genome that affect adaptation might also be strong. Here, we address the question of how the relative importance of stochastic vs. selective structuring forces varies at different geographical scales. We use microsatellite DNA and allozyme data from samples of the marine rocky shore snail Littorina saxatilis over distance scales ranging from metres to 1000 km, and we show that genetic drift is the most important structuring evolutionary force at distances > 1 km. On smaller geographical scales (< 1 km), divergent selection between contrasting habitats affects population genetic structure by impeding gene flow over microhabitat borders (microsatellite structure), or by directly favouring specific alleles of selected loci (allozyme structure). The results suggest that evolutionary drivers of population genetic structure cannot a priori be assumed to be equally important at different geographical scales. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 31–40.     

Intron sequences of arginine kinase in an intertidal snail suggest an ecotype-specific selective sweep and a gene duplication

Many species with restricted gene flow repeatedly respond similarly to local selection pressures. To fully understand the genetic mechanisms behind this process, the phylogeographic history of the species (inferred from neutral markers) as well as the loci under selection need to be known. Here we sequenced an intron in the arginine kinase gene (Ark), which shows strong clinal variation between two locally adapted ecotypes of the flat periwinkle, Littorina fabalis. The ‘small-sheltered'' ecotype was almost fixed for one haplotype, H1, in populations on both sides of the North Sea, unlike the ‘large-moderately exposed ecotype'', which segregated for ten different haplotypes. This contrasts with neutral markers, where the two ecotypes are equally variable. H1 could have been driven to high frequency in an ancestral population and then repeatedly spread to sheltered habitats due to local selection pressures with the colonization of both sides of the North Sea, after the last glacial maximum (∼18 000 years ago). An alternative explanation is that a positively selected mutation, in or linked to Ark, arose after the range expansion and secondarily spread through sheltered populations throughout the distribution range, causing this ecotype to evolve in a concerted fashion. Also, we were able to sequence up to four haplotypes consistently from some individuals, suggesting a gene duplication in Ark.     

Sympatric and Allopatric Divergence of MHC Genes in Threespine Stickleback

Parasites can strongly affect the evolution of their hosts, but their effects on host diversification are less clear. In theory, contrasting parasite communities in different foraging habitats could generate divergent selection on hosts and promote ecological speciation. Immune systems are costly to maintain, adaptable, and an important component of individual fitness. As a result, immune system genes, such as those of the Major Histocompatability Complex (MHC), can change rapidly in response to parasite-mediated selection. In threespine stickleback (Gasterosteus aculeatus), as well as in other vertebrates, MHC genes have been linked with female mating preference, suggesting that divergent selection acting on MHC genes might influence speciation. Here, we examined genetic variation at MHC Class II loci of sticklebacks from two lakes with a limnetic and benthic species pair, and two lakes with a single species. In both lakes with species pairs, limnetics and benthics differed in their composition of MHC alleles, and limnetics had fewer MHC alleles per individual than benthics. Similar to the limnetics, the allopatric population with a pelagic phenotype had few MHC alleles per individual, suggesting a correlation between MHC genotype and foraging habitat. Using a simulation model we show that the diversity and composition of MHC alleles in a sympatric species pair depends on the amount of assortative mating and on the strength of parasite-mediated selection in adjacent foraging habitats. Our results indicate parallel divergence in the number of MHC alleles between sympatric stickleback species, possibly resulting from the contrasting parasite communities in littoral and pelagic habitats of lakes.     

The adaptive role of Phosphoglucomutase and other allozymes in a marine snail across the vertical rocky-shore gradient

Natural selection can play an important role in the maintenance of genetic polymorphisms, despite ongoing gene flow. In the present study, we use previously analysed allozymic loci and perform an F _ST outlier-based analysis to detect the signatures of divergent selection between sympatric ecotypes of the marine snail Littorina saxatilis at different localities. The results obtained show that different allozyme polymorphisms are affected (directly or indirectly) by selection at distinct geographical regions. The Phosmogluco mutase-2 locus was the best candidate for adaptation and further biochemical analyses were performed. The kinetic properties of the three more common genotypes of Pgm-2 were studied. The results obtained are concordant with two alternative hypotheses: (1) natural selection is acting directly on this locus or, more probably, (2) selection is affecting a genomic region tightly linked to the enzyme locus. In both cases, the known existence of a parallel and partially independent origin of these ecotypes would explain why different candidate loci were detected in different localities.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 225–233.     

SYMPATRIC SPECIATION VIA HABITAT SPECIALIZATION DRIVEN BY DELETERIOUS MUTATIONS

Theoretical studies have suggested that the evolution of habitat (host) races, regarded as a prelude to sympatric speciation, requires strong trade-offs in adaptation to different habitats: alleles that improve fitness in some habitats and have deleterious effects of similar magnitude in other habitats must be segregating in the population. I argue that such trade-offs are not necessary; the evolution of habitat races can also be driven by genetic variation due to loci that affect fitness in one habitat and are neutral or nearly so in others, that is, when performance in different habitats is genetically independent. One source of such genetic variation are deleterious mutations with habitat-specific fitness effects. I use deterministic two-locus and multilocus models to show that the presence of such mutations in the gene pool results in indirect selection favoring habitat fidelity or habitat preference over acceptance of both suitable habitats. This leads to the evolution of largely genetically isolated populations that use different habitats, from a single panmictic population of individuals accepting both habitats. This study suggests that the conditions favoring habitat race formation, and thus possibly sympatric speciation, are much less stringent than previously thought.     

Sexual selection and assortative mating by size and their roles in the maintenance of a polymorphism in Swedish Littorina saxatilis populations

Two independent components of mating behaviour, sexual selection and assortative mating, were studied in two allopatric morphs, one sheltered boulder shore form (S-morph) and one exposed cliff shore form (E-morph), of Littorina saxatilis from the west coast of Sweden. Sexual selection was studied by comparing the sizes of copulating and non-copulating snails in the field. Size assortative mating was studied by collecting copulating pairs in the field, while assortative mating between morphs was investigated by bringing the pure morphs together in intermediary habitats and then noting the matings. The S-morph mated randomly in relation to size in two of the studied populations and exhibited a trend towards size assortative mating in a third, while the E-morph showed size assortative mating in both studied populations. The microdistribution of sizes of snails on the shores could not explain all the size assortative mating found, and instead it is argued that a size-based mate rejection behaviour also contributes to the assortative mating in at least some of these populations. There was sexual selection on size in both males and females in the S-morph, with large individuals being favoured as mates. In contrast, copulating snails of the E-morph were smaller than non-copulating ones. The significantly different sexual selection intensities between the two morphs may help to explain the size differences between them. There was random mating between the E- and the S-morphs of L. saxatilis, which suggests no incipient reproductive isolation between morphs on Swedish rocky shores. This is in agreement with earlier studies of Swedish populations, but is in contrast to the situation found in other geographical areas.     

A GENETIC ANALYSIS OF SYMPATRIC SUBPOPULATIONS OF THE SIBLING SPECIES LITTORINA SAXATILIS (OLIVI) AND LITTORINA ARCANA HANNAFORD ELLIS

A detailed analysis of shell characteristics (four parameters)and allozyme allele frequencies (eight to nine polymorphic loci)of sympatric Littorina saxatilis and Littorina arcana subpopulationsfrom Filey Brigg on the east coast of England supports theiridentification as closely related but separate species. Oneenzyme system, esterase-2, was nearly species diagnostic. Resultswere generally comparable with an earlier study of the two speciesfrom East Lothian, Scotland. Littorina saxatilis was again foundto be more variable than L. arcana. (Received 14 December 1983;     

On the neutrality of molecular genetic markers: pedigree analysis of genetic variation in fragmented populations

Many studies employ molecular markers to infer ecological and evolutionary processes, assuming that variation found at genetic loci offers a reliable representation of stochastic events in natural populations. Increasingly, evidence emerges that molecular markers might not always be selectively neutral. However, only a few studies have analysed how deviations from neutrality could affect estimates of genetic variation, using populations with known genealogy. We monitored changes in allozyme variation over eight generations in captive metapopulations of the butterfly Bicyclus anynana. Population demography was recorded by individually marking 35 000 butterflies and constructing pedigrees. We designed a computer program that simulated the inheritance of founder allozyme alleles in butterfly pedigrees. We thus tested whether the observed transmission of allozyme alleles could be explained by random genetic drift alone, or whether there was evidence for positive or negative selection. This analysis showed that in the smallest metapopulations the loss of allozyme variation exceeded the neutral rate. Possibly, linkage disequilibria between deleterious mutations and marker alleles resulted in background selection and a faster erosion of allozyme variation. In larger metapopulations, one locus (MDH) showed a significant heterozygote excess and smaller than expected loss in heterozygosity, observations consistent with (associative) overdominance. This study demonstrates that the neutrality of molecular markers cannot always be assumed, particularly in small populations with a high mutation load.     

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

PHENOTYPIC DIFFERENTIATION IN LITTORINA SAXATILIS OLIVI (MOLLUSCA, PROSOBRANCHIA) IN A SMALL AREA ON THE SWEDISH WEST COAST

The Littorina saxatilis species-complex exhibits two distinctmorphs in the area investigated. Differences in shell morphologyand number of penial glands are statistically significant (P<0.02)between populations from exposed and sheltered habitats. However,in intermediate environments intermediate specimens as wellas typical specimens of both morphs are found. The differencesare considered intraspecific and at least partly geneticallydetermined, and possible selection pressures responsible arediscussed. (Received 8 February 1981;