Maintenance of clinal variation for shell colour phenotype in the flat periwinkle Littorina obtusata

Clines can signal spatially varying selection and therefore have long been used to investigate the role of environmental heterogeneity in maintaining genetic variation. However, clinal patterns alone are not sufficient to reject neutrality or to establish the mechanism of selection. Indirect, inferential methods can be used to address neutrality and mechanism, but fully understanding the adaptive significance of clinal variation ultimately requires a direct approach. Ecological model systems such as the rocky intertidal provide a useful context for direct experimentation and can serve as a complement to studies in more traditional genetic model systems. In this study, we use indirect and direct approaches to investigate the role of environmental heterogeneity in the maintenance of shell colour polymorphism in the flat periwinkle snail, Littorina obtusata. We document replicated clines in shell colour morph frequencies over thermal gradients at two spatial scales, contrasting with patterns at previously reported microsatellite loci. In addition, experimental results demonstrate that that shell colour has predictable effects on shell temperature and that these differences in temperature, in turn, coincide with patterns of survivorship under episodic thermal stress. Direct manipulation of shell colour revealed that shell colour, and not a correlated character, was the target of selection. Our study provides evidence that spatially varying selection via thermal regime contributes to the maintenance of shell colour phenotype variation in L. obtusata in the sampled areas of the Gulf of Maine.     

Change and stability in a steep morph‐frequency cline in the snail Cepaea nemoralis (L.) over 43 years

Populations of the polymorphic land snail Cepaea nemoralis (L.) from Deepdale, Derbyshire, UK, sampled in 1965–67, showed a pattern of area effects, with steep clines among groups of populations differing in shell colour and banding morph frequencies. In 2010, most of these populations were resampled. In particular, a continuous transect made in 1967 of 42 quadrats (18.34 × 18.34 m) across a steep cline in several morph frequencies was completely resampled. In the dale as a whole, yellow shells had increased in frequency. In the transect, the frequencies of banding morphs showed no significant changes, although colour morphs showed some changes. Pink shells had increased in frequency in a section in which scrub had developed, and brown shells had increased in frequency in the area in which they had originally been at the highest frequency. In each case, the selection coefficients were of the order of 4%. Yellow had increased elsewhere. Nevertheless, both in the dale as a whole and in the transect, the pattern of geographical change in morph frequencies had remained essentially the same. Estimates of migration based on previous studies of marked snails and on modelling of the effect of drift and migration suggest that, regardless of whether the cline is a product of differential selection or of the gradual merging of previously separate founding populations, it has been in existence for a long time, and that migration occurs over greater distances than estimated from direct observation on marked snails. Although we can demonstrate that selection has occurred, the origin and maintenance of the cline and others similar to it remain in doubt; the development and maintenance of polymorphism in this species may require consideration of several processes operating on different time scales. © 2012 The Linnean Society of London     

Colour polymorphism in the intertidal snail Littorina rudis Maton

The tictors affecting the variation of shell colour morph frequencies of the intertidal snail Littorina rudis Maton have been examined using ordination and correlation techniques. The major gradient underlying variation in shell colour is related to habitat type and structure. Dark, patterned shells are characteristic of cliff and salt marsh habitats whilst light, unpatterned shells reach high frequencies on boulder shores. Colour morphs do not differ significantly in shell thickness and between-shore variation in colour morph frequencies cannot be directly attributed to habitat dilterences in the risk of shell injury, but to other factors associated with habitat type and structure, e.g. parasitism. It is suggested that the maintenance of shell colour polymorphism in L. rudis is probably mediated by selection on pleiotropic characters rather than shell colour alone.     

Irish Cepaea nemoralis Land Snails Have a Cryptic Franco-Iberian Origin That Is Most Easily Explained by the Movements of Mesolithic Humans

The origins of flora and fauna that are only found in Ireland and Iberia, but which are absent from intervening countries, is one of the enduring questions of biogeography. As Southern French, Iberian and Irish populations of the land snail Cepaea nemoralis sometimes have a similar shell character, we used mitochondrial phylogenies to begin to understand if there is a shared “Lusitanian” history. Although much of Europe contains snails with A and D lineages, by far the majority of Irish individuals have a lineage, C, that in mainland Europe was only found in a restricted region of the Eastern Pyrenees. A past extinction of lineage C in the rest of Europe cannot be ruled out, but as there is a more than 8000 year continuous record of Cepaea fossils in Ireland, the species has long been a food source in the Pyrenees, and the Garonne river that flanks the Pyrenees is an ancient human route to the Atlantic, then we suggest that the unusual distribution of the C lineage is most easily explained by the movements of Mesolithic humans. If other Irish species have a similarly cryptic Lusitanian element, then this raises the possibility of a more widespread and significant pattern.     

Quantitative measures and 3D shell models reveal interactions between bands and their position on growing snail shells

The nature of shell growth in gastropods is useful because it preserves the ontogeny of shape, colour, and banding patterns, making them an ideal system for understanding how inherited variation develops, is established and maintained within a population. However, qualitative scoring of inherited shell characters means there is a lack of knowledge regarding the mechanisms that control fine variation. Here, we combine empirical measures of quantitative variation and 3D modeling of shells to understand how bands are placed and interact. By comparing five‐banded Cepaea individuals to shells lacking individual bands, we show that individual band absence has minor but significant impacts upon the position of remaining bands, implying that the locus controlling band presence/absence mainly acts after position is established. Then, we show that the shell grows at a similar rate, except for the region below the lowermost band. This demonstrates that wider bands of Cepaea are not an artifact of greater shell growth on the lower shell; they begin wider and grow at the same rate as other bands. Finally, we show that 3D models of shell shape and banding pattern, inferred from 2D photos using ShellShaper software, are congruent with empirical measures. This work therefore establishes a method that may be used for comparative studies of quantitative banding variation in snail shells, extraction of growth parameters, and morphometrics. In the future, studies that link the banding phenotype to the network of shell matrix proteins involved in biomineralization and patterning may ultimately aid in understanding the diversity of shell forms found in molluscs.     

SHELL STRENGTH OF COLOUR MORPHS OF THE MANGROVE SNAIL LITTORARIA PALLESCENS

The load required to break shells of the snail Littor-aria pallescens(Philippi) (Gastropoda: Littorinidae) has been established.Mass increases more than proportionally with size. Size forsize, relatively heavy shells are thicker than relatively lightones, and the thicker shells are stronger. Three samples fromdifferent localities differ significantly in strength. Littorariapallescens is polymorphic for shell colour. It is shown thatin two of the three samples, both from Papua New Guinea, individualsof the dark morph are significantly heavier and stronger thanthose of the yellow morph, while there is no difference betweenmorphs in the third sample, from Thailand. In the sample withthe most robust shells the dark individuals are stronger, massfor mass, than yellows. Thus, colour morph is sometimes, butnot always, associated with differences in shell characteristicswhich are likely to affect fitness. Where there is a difference,the dark morph is estimated to be on average about 22% strongerthan the yellow morph. The difference between morphs is smallcompared with the large variation in shell characters associatedwith size. (Received 20 May 1992; accepted 22 June 1992)     

Geographic structuring of molecular and morphological polymorphism in Pyrenean populations of the snail Cepaea nemoralis

To determine whether geographic patterns of variation in the frequency of shell polymorphisms in the land snail Cepaea nemoralis mark regions of extensive genetic differentiation, allele frequencies at six polymorphic enzyme-encoding loci were analyzed electrophoretically in 74 samples from the central part of the Spanish Pyrenees. Within a large area effect for unbanded shells on the southern slope of the Pyrenees, there is no enzyme locus at which allele frequencies are homogeneous. Because this morphological area effect does not mark a region of pervasive genetic change, the genetic structure of populations in the region studied does not conform to recent models of area effect speciation. Factor analysis of allele frequencies at three loci controlling shell polymorphisms and six enzyme loci separated populations on the southern slope into two groups, one of which is similar to the group of populations on the northern slope. There is therefore considerable geographic structuring of both molecular and morphological polymorphism in this species.     

Gudgeon fish with and without genetically determined countershading coexist in heterogeneous littoral environments of an ancient lake

Countershading, characterized by a darker dorsal surface and lighter ventral surface, is common among many animals. This dorsoventral pigment polarity is often thought to be adaptive coloration for camouflage. By contrast, noncountershaded (melanistic) morphs often occur within a species due to genetic color polymorphism in terrestrial animals. However, the polymorphism with either countershaded or melanistic morphs is poorly known in wild aquatic animals. This study explored the genetic nature of diverged color morphs of a lineage of gudgeon fish (genus Sarcocheilichthys) in the ancient Lake Biwa and propose this system as a novel model for testing hypotheses of functional aspects of countershading and its loss in aquatic environments. This system harbors two color morphs that have been treated taxonomically as separate species; Sarcocheilichthys variegatus microoculus which occurs throughout the littoral zone and Sarcocheilichthys biwaensis which occurs in and around rocky areas. First, we confirmed that the divergence of dorsoventral color patterns between the two morphs is under strict genetic control at the levels of chromatophore distribution and melanin‐related gene expression under common garden rearing. The former morph displayed sharp countershading coloration, whereas the latter morph exhibited a strong tendency toward its loss. The crossing results indicated that this divergence was likely controlled by a single locus in a two‐allele Mendelian inheritance pattern. Furthermore, our population genomic and genome‐wide association study analyses detected no genome‐wide divergence between the two morphs, except for one region near a locus that may be associated with the color divergence. Thus, these morphs are either in a state of intraspecific color polymorphism or two incipient species. Evolutionary forces underlying this polymorphism appear to be associated with heterogeneous littoral environments in this lake. Future ecological genomic research will provide insight into adaptive functions of this widespread coloration, including the eco‐evolutionary drivers of its loss, in the aquatic world.     

Differences in predatory pressure on terrestrial snails by birds and mammals

The evolution of shell polymorphism in terrestrial snails is a classic textbook example of the effect of natural selection in which avian and mammalian predation represents an important selective force on gene frequency. However, many questions about predation remain unclear, especially in the case of mammals. We collected 2000 specimens from eight terrestrial gastropod species to investigate the predation pressure exerted by birds and mice on snails. We found evidence of avian and mammalian predation in 26.5% and 36.8% of the shells. Both birds and mammals were selective with respect to snail species, size and morphs. Birds preferred the brown-lipped banded snail Cepaea nemoralis (L.) and mice preferred the burgundy snail Helix pomatia L. Mice avoided pink mid-banded C. nemoralis and preferred brown mid-banded morphs, which were neglected by birds. In contrast to mice, birds chose larger individuals. Significant differences in their predatory pressure can influence the evolution and maintenance of shell size and polymorphism of shell colouration in snails.     

Patterns of morphological and molecular polymorphism in the land snail Cepaea nemoralis

In the land snail Cepaea nemoralis , allele frequencies at loci controlling shell polymorphisms often show large areas of remarkable constancy which are separated by steep clines from neighbouring areas with strikingly different allele frequencies. It has recently been claimed that these 'area effects' exemplify a general tendency for population differentiation without geographic isolation in a variety of organisms of relatively low mobility. As such they could represent an early phase in the process of speciation. If this is true, population differentiation of shell polymorphisms in Cepaea would be expected to be accompanied by parallel differentiation at other gene loci, such as those detected by gel electrophoresis.
We have studied populations of C. nemoralis in North Wales and in the Valle de Aran of the Pyrenees. Levels of molecular heterogeneity are comparable to those found in related animals which show much less visible polymorphism. In spite of some statistical problems inherent in the analysis ol overlapping geographic patterns, there is no clear association between the patterns of geographic variation at the visible and molecular levels. Claims that genetic differentiation in the visible polymorphisms between C. nemoralis populations are a special case of the formation of geographic races are therefore probably not justified.     

Influence of shell colour on the mortality rate of the land snail Arianta arbustorum under different microclimatic regimes

The mortality of phenotypic shell colour morphs and age classes of the snail Arianta arbustorum in several microclimatic conditions was recorded. An analysis of variance was performed on five factors: adaptation temperature, relative humidity, shell colour, age class and test temperature. There were no significant differences in the mortality between different adaptation temperatures or relative humidities, but the interaction of these two factors was highly significant. There were significant differences in mortality rate between test temperature and age class. The mortality of the brown morph was higher than the yellow one at all adaptation temperatures (overall, P < 0.05). There were no significant differences between the mortality rates of the two morphs at different relative humidities. The mortality of the brown morph was higher than that of the yellow at six out of seven test temperatures. Juvenile snails survived significantly better than adults at all test temperatures.     

Evolutionary change in Cepaea nemoralis shell colour over 43 years

We compared shell colour forms in the land snail Cepaea nemoralis at 16 sites in a 7 × 8 km section of the Province of Groningen, the Netherlands, between 1967 and 2010. To do so, we used stored samples in a natural history collection and resampled the exact collection localities. We found that almost all populations had experienced considerable evolutionary change in various phenotypes, possibly due to population bottlenecks and habitat change after repeated land consolidation schemes in the area. More importantly, we found a consistent increase in yellow effectively unbanded snails at the expense of brown snails. This is one of the expected adaptations to climate change (this area of the Netherlands has warmed by 1.5–2.0 °C over the time period spanned by the two sampling years), and the first clear demonstration of this in C. nemoralis.     

Gene frequencies, genetic background and environment in Pyrenean populations of Cepaea nemoralis (L.)

Microgeographical genetic differentiation of populations of the polymorphic land snail Cepaea nemoralis (L.) has been the subject of intensive investigation. The importance of environmental selection acting at the loci which control the shell polymorphism in determining the distribution of genes in this species remains a matter of considerable controversy. Some populations of C. nemoralis in the central Pyrenees show striking local associations of morph frequency (and particularly of the frequency of unhanded shells) with topography which have been ascribed to the action of natural selection.
In the Ter valley (eastern Spanish Pyrenees) populations described here, however, the pattern of association of the frequencies of the shell morphs with topography is quite different from that found in the central Pyrenees. In particular, the distribution of unhanded shells is almost the inverse of that found in some previously studied populations. The general picture of variation in relation to topography in Pyrenean C. nemoralis populations is now seen to be one of great local consistency within valley systems, but of considerable differences between valley systems. This variation in population structure may be due either to cryptic and undetected environmental differences between valley systems or to evolutionary divergence between their C. nemoralis populations. The difficulty of deciding between these two conflicting hypotheses emphasizes the limitations of the technique of simply examining the distribution of individual genes in natural populations as a method of analysing their evolutionary history.     

A countrywide survey of colour morph frequencies in the spider Enoplognatha ovata (Glerck) (Araneae: Theridiidae): evidence for natural selection

Data on the frequencies of the main colour morphs ( lineata, redimita and ovata ) of Enoplognatha ovata have been collected from a total of 454 Ordnance Survey 10 km squares distributed throughout Great Britain. Only c. 0.5% of samples were monomorphic, for lineata in each case. Multiple regression analysis has been used to assess possible associations between morph frequencies and principal components derived from nine environmental variables. The distribution of morph frequencies is not random but shows weak clines associated with certain climatic factors. These large scale clines indicate the action of natural selection although very local variations in morph frequencies may result from selection and/or drift.     

Rapid,habitat-related evolution of land snail colour morphs on reclaimed land

I made use of the known dates of reclamation (and of afforestations) in the IJsselmeerpolders in The Netherlands to assess evolutionary adaptation in Cepaea nemoralis. At 12 localities (three in each polder), I sampled a total of 4390 adult individuals in paired open and shaded habitats, on average 233 m apart, and scored these for genetic shell colour polymorphisms. The results show (highly) significant differentiation at most localities, although the genes involved differed per locality. Overall, though, populations in shaded habitats had evolved towards darker shells than those in adjacent open habitats, whereas a ‘Cain & Sheppard'' diagram (proportion yellow shells plotted against ‘effectively unbanded'' shells) failed to reveal a clear pattern. This might suggest that thermal selection is more important than visual selection in generating this pattern. Trait differentiation, regardless of whether they were plotted against polder age or habitat age, showed a linear increase of differentiation with time, corresponding to a mean rate of trait evolution of 15–31 kilodarwin. In conclusion, C. nemoralis is capable of rapid and considerable evolutionary differentiation over 1–25 snail generations, though equilibrium may be reached only at longer time scales.     

Citizen science reveals unexpected continental-scale evolutionary change in a model organism

Organisms provide some of the most sensitive indicators of climate change and evolutionary responses are becoming apparent in species with short generation times. Large datasets on genetic polymorphism that can provide an historical benchmark against which to test for recent evolutionary responses are very rare, but an exception is found in the brown-lipped banded snail (Cepaea nemoralis). This species is sensitive to its thermal environment and exhibits several polymorphisms of shell colour and banding pattern affecting shell albedo in the majority of populations within its native range in Europe. We tested for evolutionary changes in shell albedo that might have been driven by the warming of the climate in Europe over the last half century by compiling an historical dataset for 6,515 native populations of C. nemoralis and comparing this with new data on nearly 3,000 populations. The new data were sampled mainly in 2009 through the Evolution MegaLab, a citizen science project that engaged thousands of volunteers in 15 countries throughout Europe in the biggest such exercise ever undertaken. A known geographic cline in the frequency of the colour phenotype with the highest albedo (yellow) was shown to have persisted and a difference in colour frequency between woodland and more open habitats was confirmed, but there was no general increase in the frequency of yellow shells. This may have been because snails adapted to a warming climate through behavioural thermoregulation. By contrast, we detected an unexpected decrease in the frequency of Unbanded shells and an increase in the Mid-banded morph. Neither of these evolutionary changes appears to be a direct response to climate change, indicating that the influence of other selective agents, possibly related to changing predation pressure and habitat change with effects on micro-climate.     

A two-stage model for Cepaea polymorphism

The history of the study of snails in the genus Cepaea is briefly outlined. Cepaea nemoralis and C. hortensis are polymorphic for genetically controlled shell colour and banding, which has been the main interest of the work covered. Random drift, selective predation and climatic selection, both at a macro- and micro-scale, all affect gene frequency. The usual approach to understanding maintenance of the polymorphism, has been to look for centripetal effects on frequency. Possible processes include balance of mutation pressure and drift, heterozygote advantage, relational balance heterosis, frequency-dependent predation, multi-niche selective balance, or some combination of these. Mutational balance is overlaid by more substantial forces. There is some evidence for heterosis. Predation by birds may protect the polymorphism, and act apostatically to favour distinct morphs. Although not substantiated for Cepaea, many studies show that predators behave in the appropriate manner, while shell colour polymorphisms in molluscs occur most commonly in species exposed to visually searching predators. It is not known whether different thermal properties of the shells help to generate equilibria. Migration between colonies is probably greater than originally thought. The present geographical range has been occupied for less than 5000 generations. Climatic and human modification alter snail habitats relatively rapidly, which in turn changes selection pressures. A simple simulation shows that migration coupled with selection which fluctuates but is not centripetal, may retain polymorphism for sufficiently long to account for the patterns we see today. There may therefore be a two-stage basis to the polymorphism, comprising long-term but weak balancing forces coupled with fluctuating selection which does not necessarily balance but results in very slow elimination. Persistence of genetic variants in this way may provide the conditions for evolution of a balanced genome.     

Genome-wide population genetic analysis identifies evolutionary forces establishing continuous population divergence

Elucidating the mechanism shaping the spatial variations of traits has long been a central concern of evolutionary biologists. Geographic clines of allele/morph frequencies along environmental gradients are suggested to be established and maintained by the balancing of two opposing evolutionary forces, namely selection that generates spatial differentiation in morph frequencies, and selection and/or stochastic factors that lead to the coexistence of multiple morphs within a population. Thus, testing for both selection and stochastic factors is necessary for a comprehensive understanding of the mechanism underlying clinal variation in morph/allele frequency in natural populations. Here, I identified the evolutionary forces responsible for clinal variation of color morph frequency in Ischnura senegalensis by comparing the population divergence of putatively neutral loci generated by high-throughput next-generation sequencing (F _STn) with that of the putative color locus (F _STc). No strong correlation was observed between F _STn and F _STc, suggesting that stochastic factors contribute less to color-locus population divergence. F _STc was less than F _STn between populations exposed to similar environmental conditions, but greater than F _STn between populations exposed to different environmental conditions, suggesting that both balancing selection and divergent selection act on the color locus. Therefore, two antagonistic selection factors rather than stochastic and historical factors contribute to establishing the clinal variation of morph frequency in I. senegalensis.     

Disequilibrium in some Cepaea populations

Phenotypic disequilibrium is a measure available, along with morph frequency, in many published sets of sample data recording the polymorphism of shell colour and pattern in the snail Cepaea nemoralis. The relation of disequilibrium to morph frequency for the colour (pink/yellow) and banding (unbanded/banded) loci has been examined for a large and widespread set of data. The direction of disequilibrium is a function of frequency at the two loci in a way that suggests that selection favours combinations of common morphs, whichever they are. This could indicate that such combinations are common because they have been selected. The data are therefore consistent with the proposal that populations are generally subject to selection of varying directions at different places and times, acting on the phenotype. In combination with migration, such selection could lead to prolonged polymorphism.     

Multi-scale altitudinal patterns in species richness of land snail communities in south-eastern France

Aim Species richness is an important feature of communities that varies along elevational gradients. Different patterns of distribution have been described in the literature for various taxonomic groups. This study aims to distinguish between species density and species richness and to describe, for land snails in south‐eastern France, the altitudinal patterns of both at different spatial scales. Location The study was conducted on five calcareous mountains in south‐eastern France (Etoile, Sainte Baume, Sainte Victoire, Ventoux and Queyras). Methods Stratified sampling according to vegetation and altitude was undertaken on five mountains, forming a composite altitudinal gradient ranging from 100 to 3100 m. Visual searching and analysis of turf samples were undertaken to collect land snail species. Species density is defined as the number of species found within quadrats of 25 m². Species richness is defined as the number of species found within an elevation zone. Different methods involving accumulation curves are used to describe the patterns in species richness. Elevation zones of different sizes are studied. Results Eighty‐seven species of land snails were recovered from 209 samples analysed during this study. Land snail species density, which can vary between 29 and 1 species per 25 m², decreases logarithmically with increasing altitude along the full gradient. However, on each mountain separately, only a linear decrease is observable. The climatic altitudinal gradient can explain a large part of this pattern, but the great variability suggests that other factors, such as heterogeneity of ground cover, also exert an influence on species density. The altitudinal pattern of species richness varies depending on the spatial resolution of the study. At fine resolution (altitudinal zones of 100 m) land snail species richness forms a plateau at altitudes below 1000 m, before decreasing with increasing altitude. At coarse resolution (altitudinal zones of 500 and 1000 m) the relationship becomes linear. Main conclusions This study reveals that land snail species density and land snail species richness form two different altitudinal patterns. Species density exhibits strong variability between sites of comparable altitude. A large number of samples seem necessary to study altitudinal patterns of species density. Species density decreases logarithmically with increasing altitude. Above a critical altitudinal threshold, this decrease lessens below the rate seen in the first 1500 m. Different methods exist to scale‐up species density to species richness but these often produce different patterns. In this study, the use of accumulation curves has yielded a pattern of species richness showing a plateau at low altitude, whereas simple plotting of known altitudinal ranges from single mountains would have produced stronger mid‐altitudinal peaks. This study shows that not only factors such as temperatures and habitat heterogeneity, but also an ecotone effect, are responsible for the observed patterns.