Open Access

Why do snails have hairs? A Bayesian inference of character evolution

  • Markus Pfenninger1Email author,
  • Magda Hrabáková2,
  • Dirk Steinke3 and
  • Aline Dèpraz4
BMC Evolutionary Biology20055:59

DOI: 10.1186/1471-2148-5-59

Received: 14 July 2005

Accepted: 04 November 2005

Published: 04 November 2005

Abstract

Background

Costly structures need to represent an adaptive advantage in order to be maintained over evolutionary times. Contrary to many other conspicuous shell ornamentations of gastropods, the haired shells of several Stylommatophoran land snails still lack a convincing adaptive explanation. In the present study, we analysed the correlation between the presence/absence of hairs and habitat conditions in the genus Trochulus in a Bayesian framework of character evolution.

Results

Haired shells appeared to be the ancestral character state, a feature most probably lost three times independently. These losses were correlated with a shift from humid to dry habitats, indicating an adaptive function of hairs in moist environments. It had been previously hypothesised that these costly protein structures of the outer shell layer facilitate the locomotion in moist habitats. Our experiments, on the contrary, showed an increased adherence of haired shells to wet surfaces.

Conclusion

We propose the hypothesis that the possession of hairs facilitates the adherence of the snails to their herbaceous food plants during foraging when humidity levels are high. The absence of hairs in some Trochulus species could thus be explained as a loss of the potential adaptive function linked to habitat shifts.

Background

Evolutionary theory predicts that costly structures must convey a fitness advantage to their bearers in order to be maintained over evolutionary time [1]. Flightlessness in birds and insects, limblessness in lizards and sightlessness in cave-dwelling organisms are some prominent examples of phenotypic regression due to the loss of adaptive function (reviewed in [2]). Molluscs in general and gastropods in particular display a fascinating diversity of elaborate shell structures [3, 4] and have attracted considerable research efforts to explain them in adaptive terms [57]. The proposed roles invoked mechanical stability [8], defence against predators [9], sexual selection [10] and climatic selection [11]. However, the potential selective advantage of hair-like shell ornamentation of certain land snail species remains unknown.

These hairs can reach varying densities (up to 20 per squaremilimetre) and lengths (up to three millimetres). In some cases hardly visible, they confer an almost furry impression to the shell in others. These semi-rigid structures are part of the periostracum, a thin protein layer (conchiolin) secreted by the snail to cover the calcareous shell [12]. Building hairs requires the snail to have specialised glandular tissue and complex strategies to form them. Consequently, this trait can be assumed to be costly and should thus present a selective advantage to its bearers in order to be conserved.

Haired shells occur in several species of the Stylommatophoran families Polygyridae, Helicidae and Hygromiidae. These families are only distantly related [13], suggesting that this features has evolved several times independently. Haired shells are almost exclusively observed in species living in moist microhabitats, like layers of fallen leaves, broad-leaved vegetation, damp meadows or wet scree [14]. Such a correlation suggests an adaptive significance of the trait in such a habitat [1]; it was thus speculated that the hygrophobic hairs facilitate the movement in wet environments by relieving surface tension [14, 15]. A correlation between haired shells and humid habitats is thus expected. In order to test this, we employed the recent Bayesian extensions of the comparative method, allowing to take mapping and phylogenetic uncertainty simultaneously into account [16]. With a diversity hotspot in South Germany, Eastern France and Switzerland, the land snail genus Trochulus s. str. (common name: Hairy snails) is particularly suited to address our question: its species exhibit variability in both hairiness and ecology. This study present the first comprehensive molecular phylogeny for the genus Trochulus Chemnitz, 1786 (until recently Trichia, Hartmann 1840) based on mitochondrial and nuclear loci. Finally, we tested experimentally whether the possession of haired shells indeed facilitates locomotion.

Results

Lineage identification and phylogenetic relations

The initial phylogenetic analysis on a COI data set of the presumed Trochulus species resolved 18 terminal clades, each with 0.99 posterior probabilities or higher (Figure 1). The uncorrected sequence divergence among those clades ranged from 0.029 to 0.173 (Table 2). Out of these lineages, nine could be assigned to existing taxa, because the species were sampled from the type locality and/or were morphologically unmistakable. The nine remaining clades, however, could not be unequivocally attributed to a taxonomic name. All eighteen identified lineages were used as molecularly defined operational taxonomic units in the subsequent analyses [17].
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-5-59/MediaObjects/12862_2005_Article_171_Fig1_HTML.jpg
Figure 1

Unrooted consensus tree of 90,000 trees sampled by the Markov-chain in Bayesian analysis for the COI-fragment. Numbers on nodes indicate the Bayesian posterior probability.

Table 2

Pairwise uncorrected COI sequence divergence among lineages and species (mean ± s.d.).

 

A

B

C

villosulus

montanus

clandest.

caelatus

D

E

F

G

H

I

biconicus

nov. spec.

villosus

B

0.111 ± 0.015

 

striolatus/plebeius

             

C

0.118 ± 0.015

0.091 ± 0.013

clade

             

villosulus

0.104 ± 0.014

0.073 ± 0.013

0.091 ±

0.014

             

montanus

0.143 ± 0.017

0.144 ± 0.017

0.138 ±

0.017

0.117 ± 0.016

            

clandestinus

0.135 ± 0.017

0.129 ± 0.016

0.140 ±

0.017

0.125 ± 0.016

0.083 ± 0.014

Jura clade

          

caelatus

0.144 ± 0.016

0.152 ± 0.018

0.141 ±

0.017

0.130 ± 0.016

0.092 ± 0.014

0.093 ± 0.013

          

D

0.146 ± 0.017

0.105 ± 0.016

0.125 ±

0.016

0.127 ± 0.017

0.142 ± 0.017

0.126 ± 0.016

0.138 ± 0.017

         

E

0.148 ± 0.017

0.109 ± 0.016

0.135 ±

0.016

0.117 ± 0.016

0.134 ± 0.016

0.104 ± 0.014

0.107 ± 0.015

0.079 ± 0.013

     

sericeus/hispidus

  

F

0.152 ± 0.017

0.112 ± 0.015

0.126 ±

0.015

0.119 ± 0.016

0.146 ± 0.016

0.113 ± 0.015

0.114 ± 0.015

0.072 ± 0.013

0.029 ± 0.008

    

clade

  

G

0.142 ± 0.017

0.139 ± 0.017

0.136 ±

0.018

0.109 ± 0.016

0.132 ± 0.017

0.117 ± 0.016

0.147 ± 0.016

0.120 ± 0.017

0.111 ± 0.016

0.105 ± 0.015

      

H

0.164 ± 0.017

0.153 ± 0.018

0.170 ±

0.017

0.169 ± 0.018

0.162 ± 0.017

0.126 ± 0.015

0.147 ± 0.017

0.116 ± 0.015

0.100 ± 0.014

0.101 ± 0.014

0.128 ± 0.016

     

I

0.150 ± 0.017

0.151 ± 0.017

0.155 ±

0.017

0.144 ± 0.017

0.148 ± 0.016

0.124 ± 0.015

0.133 ± 0.016

0.106 ± 0.014

0.102 ± 0.014

0.099 ± 0.013

0.106 ± 0.014

0.061 ± 0.010

    

biconicus

0.158 ± 0.018

0.136 ± 0.017

0.141 ±

0.018

0.128 ± 0.018

0.164 ± 0.020

0.154 ± 0.018

0.167 ± 0.019

0.145 ± 0.017

0.145 ± 0.018

0.145 ± 0.018

0.143 ± 0.018

0.179 ± 0.018

0.167 ± 0.017

   

nov. spec.

0.153 ± 0.017

0.137 ± 0.017

0.151 ±

0.017

0.142 ± 0.017

0.151 ± 0.016

0.128 ± 0.015

0.137 ± 0.016

0.100 ± 0.015

0.082 ± 0.014

0.081 ± 0.013

0.106 ± 0.015

0.080 ± 0.013

0.057 ± 0.011

0.161 ±

0.018

  

villosus

0.144 ± 0.018

0.139 ± 0.018

0.154 ±

0.018

0.152 ± 0.018

0.160 ± 0.018

0.153 ± 0.018

0.160 ± 0.018

0.168 ± 0.019

0.149 ± 0.018

0.151 ± 0.018

0.173 ± 0.019

0.171 ± 0.019

0.165 ± 0.018

0.163 ±

0.019

0.168 ± 0.019

villosa/alpicola

alpicolus

0.142 ± 0.016

0.141 ± 0.015

0.151 ±

0.018

0.153 ± 0.017

0.159 ± 0.014

0.151 ± 0.019

0.162 ± 0.019

0.167 ± 0.018

0.147 ± 0.016

0.150 ± 0.019

0.173 ± 0.019

0.173 ± 0.020

0.162 ± 0.016

0.161 ±

0.020

0.170 ± 0.020

0.006 ±

0.004

The Bayesian phylogenetic analysis of the entire data set (COI, 16S and ITS-1) showed the monophyly of the genus Trochulus within the Hygromiinae with high posterior probability, except for T. lubomirskii, which seems to be only distantly related to this genus (Figure 2). In addition to the early branching T. villosus/alpicolus clade, the genus is composed of three well supported subclades: first, a clade containing the T. striolatus/plebeius-like lineages together with T. villosulus, a second clade with ecologically divers species confined to the Jura mountains and the neighbouring Mittelland-plain and finally, a T. hispidus/sericeus-like clade, containing also T. biconicus and a new species.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-5-59/MediaObjects/12862_2005_Article_171_Fig2_HTML.jpg
Figure 2

Consensus tree of 90,000 trees sampled by the Markov-chain in Bayesian analysis for the total data set (1383 bp of COI, 16S and ITS1). Numbers on nodes indicate the Bayesian posterior probabilities.

Correlation of shell hairiness with habitat

The PCA on habitat humidity describing variables resulted in two meaningful axes, representing 79.7% and 13.4% of the total variation. The first component opposed sampling sites in shady woods and sites in sun exposed, open areas. This axis can therefore be interpreted as an evaporation gradient. The second axis is a gradient of the summer precipitation on one hand and the humidity demand of the vegetation on the other (Figure 3). It can thus be considered as a humidity gradient. The sampling sites appear as two distinct clusters that could be classified as either moist or dry (Table 1). The outlier (TA) was also considered to be humid, according to its high humidity levels. For each population, at least ten adult individuals were scored for the presence or absence of hairs (mixed populations were not found). Non-haired populations exclusively corresponded to species described in the literature as having smooth shells (Table 1). When plotting the hairiness of each population on the PCA, a complete congruence between humidity and hairiness became apparent: haired shells tended to occur at sites with low evaporation and/or high precipitation while smooth shells were found at places with high evaporation and/or lower precipitation (Figure 3).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-5-59/MediaObjects/12862_2005_Article_171_Fig3_HTML.jpg
Figure 3

Two first components of the PCA of the sampled localities over 9 environmental variables. Black dots: populations with haired individuals. Open squares: populations with hair-less individuals. Sampling sites above the dotted line are considered moist whereas those under it are dry.

Table 1

Table of sampling sites, presumed taxon, habitat characterisation and presence or absence of hairs.

Sampling site

Abbreviation

Geographical position

Presumed taxon

Habitat

Humidity

Hairs

Burgsinn, Bayern, Germany

D3

50°09' 31"N 09°40'34"E

T. striolatus/plebeius

wood

moist

yes

Habichtstal, Bayern, Germany

D4

50°02'54"N 09°25'41"E

T. striolatus/plebeius

wood

moist

yes

Dommershausen, Rheinland-Pfalz, Germany

D5

50°07'45"N 07°23'47"E

T. striolatus/plebeius

wood

moist

yes

Bingen, Rheinland-Pfalz, Germany

D6

49°55'56"N 07°58'57"E

T. striolatus/plebeius

wood

moist

yes

Eltville, Hessen, Germany

D7

50°0059"N 08°04'28"E

T. striolatus/plebeius T. sericeus/hispidus

wood

moist

yes

Büchsenberg, Baden-Württemberg, Germany

D8

48°04'55"N 07°37'23"E

T. striolatus/plebeius

wood

moist

yes

St. Seine l'Abbaye, Côte d'Or, France

MdO

47°26'04"N 04°46'55"E

T. sericeus/hispidus

wood

-

yes

La Neirigue, Fribourg, Switzerland

CH3

46°42'16"N 06°55'12"E

T. clandestinus

riverbank vegetation

moist

no

Barrage des Rossens, Fribourg, Switzerland

CH6

46°43'33"N 07°06'55"E

T. sericeus/hispidus

wood

moist

yes

Gorges de la Jogne, Fribourg, Switzerland

CH8

46°36'45"N 07°07'10"E

T. sericeus/hispidus

gorge

moist

yes

Ste Croix, Vaud, Switzerland

CH10

46°50'44"N 06°32'02"E

T. montanus

open wood

dry

no

La Côte aux Fées, Neuchâtel, Switzerland

CH11+12

46°50'74"N 06°32'42"E

T. montanus

grassland

dry

no

Col des Mosses, Vaud, Switzerland

CH13

46°25'49"N 07°08'22"E

T. villosus

wood

moist

yes

Vallée du Rhône, Vaud, Switzerland

CH15

46°19'43"N 06°13'39"E

T. sericeus/hispidus

wood

moist

yes

Sensetal, Bern, Switzerland

CH18

46°49'46"N 07°19'19"E

T. clandestinus

riverbank vegetation

dry

no

Birseschlucht, Bern, Switzerland

CH21

47°17'54"N 07°23'00"E

T. caelatus

cliff

dry

no

Birseschlucht, Bern, Switzerland

CH24

47°16'56"N 07°23'13"E

T. sericeus/hispidus

wood

moist

yes

Château d'Oex, Vaud, Switzerland

CHAT

46°16'23"N 07°21'42"E

T. nov. spec.

alpine meadow

dry

no

Bannalppass, Nidwalden, Switzerland

TA

46°53'40"N 08°27'15"E

T. alpicolus

alpine meadow

moist

yes

Bannalppass, Nidwalden, Switzerland

TB

46°53'43"N 08°27'21"E

T. biconicus

alpine meadow

dry

no

Velká Javořina, Velká nad Veličkou, Czech Republic

TVJ

48°51'26"N 17°39'11"E

T. villosulus

wood

moist

yes

Bohuslavice u Zlína, Czech Republic

PLB

49°09'19"N 17°37'29"E

T. lubomirskii

meadow

-

yes

Character state evolution

As the occurrence in moist habitats was systematically linked to the presence of hairs in Trochulus s.str., only a single analysis was necessary for both characters. The Bayesian analysis of character evolution suggested with high posterior probability that the most recent common ancestor of the genus Trochulus most likely possessed hairs and lived in a moist habitat (Figure 4). The analysis also revealed considerable mapping- and/or phylogenetic uncertainty in the reconstruction of crucial ancestral nodes (nodes 1–3 in Figure 4). The average Bayesian parameter estimate for the character change ratio was 2.50 ± 0.11 (mean ± s.d.), indicating that a loss of hairs associated with a transition from wet to dry habitats occurred more frequently than vice versa. This was in concordance with the parsimony reconstruction of character state changes on all different topologies of the 99% credibility set of trees. A minimum number of three independent losses of hairs / habitat transitions had a higher probability (0.59) than the only other observed pattern of two losses/one gain or three losses/no gain (0.41).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-5-59/MediaObjects/12862_2005_Article_171_Fig4_HTML.jpg
Figure 4

Bayesian reconstruction of ancestral states (hairs/no hairs, moist/dry habitat, respectively) on the topology of the Bayesian consensus tree (restricted to the Trochulus-clade).

Functional analysis

The analysis of variance showed that on a water-covered leaf surface, hairy shells required a significantly higher minimum force to overcome the adhesion (F = 720, d.f. = 2, p < 0.00001). There was no difference on a dry surface (F = 0.47, d.f. = 2, p = 0.37; Figure 5).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-5-59/MediaObjects/12862_2005_Article_171_Fig5_HTML.jpg
Figure 5

Mean (+/- s.d.) minimum necessary force to move haired and smooth shells on wet and dry leaf surfaces.

Discussion

Considering the limited number of sites sampled, we found a relatively large number of lineages, most of which could not be attributed to described species. This suggests that many other more or less morphologically similar entities may exist throughout the range of the genus. The existence of cryptic lineages could explain at least in part the current taxonomic uncertainty in Trochulus [1820]. For example, several subspecies have been described for T. striolatus [21], which may well represent distinct evolutionary lineages such as described here. Given that the sequence divergence among the nine unidentified lineages is of the same magnitude as among described, morphologically and ecologically distinct species (Figure 1, Table 2), it can be reasonably assumed that the cryptic lineages within the striolatus/plebeius and hispidus/sericeus clades correspond to good species. Even under the assumption of an exceptionally fast molecular clock in land snails of up to 5% sequence divergence per one million years [22], the lineages in the striolatus/plebeius clade, for example, persisted for at least two million years as independent evolutionary entities. The existence of more or less cryptic lineages or species is not an unusual finding in land snails [2325]. In contrast to the high divergence of the unidentified lineages, the comparatively small genetic distance between T. villosus and T. alpicolus indicated a questionable specific distinction between these two taxa. Detailed phylogeographic analyses in addition to morphometric and ecological studies will be necessary to disentangle the species limits of these cryptic Trochulus complexes, clarify the taxonomy and reveal their evolutionary history. In addition, the species T. lubomirskii, which was placed by Schileyko [26] into the subgenus Plicuteria, may not belong to the genus Trochulus at all.

A haired shell appears as the ancestral state in the genus Trochulus. This inference is strengthened by the observation that some of the hair-less species do possess some as juveniles. During the evolutionary history of the genus Trochulus, hairs appear to have been lost several times independently (Fig. 3, Table 1) and this was always correlated with a shift in habitat (i.e. hairs are only present in moist habitats, mostly woodlands). This suggests that hairs potentially have an adaptive function in humid habitats and once the presumed selective pressure for the maintenance of these costly protein structures is relieved, they are lost. Such a correlation makes certain potential adaptive explanations for hairiness unlikely: defence against predators or mechanical stability have no obvious reasons to co-vary with the humidity characteristics of a habitat.

The facilitation of locomotion by decreasing the adhesion to water films in humid environments had been previously hypothesised to be the selective advantage of a haired shell [14, 15]. However, the results of our experiments have shown that the opposite is true. The presence of hairs significantly increased the minimum force necessary to move shells over wet surfaces. Having thus shown that the initial hypothesis [14] is at least in this case not applicable, we propose an alternative: haired shells may confer an selective advantage by increasing the adhesion to the water film on the unstable, moving leaves of their feeding plants during foraging (Figure 5). Indeed, snails are mostly active during phases of high ambient humidity [27] when leaves are covered with a water film due to rain, fog or dew. This water film is usually in contact with the shell during locomotion (Figure 6). Observation shows that Trochulus species in moist habitats preferentially forage on large-leaved herbaceous plants like Adenostyles, Urtica, Homogyne or Tussilago [28]. Hence, falling off the leaf and needing to crawl up again to this feeding site (that can be one meter above ground) represents a considerable effort given the exceedingly costly and ineffective locomotion of land snails [29]. In dry habitats on the contrary, snail species avoid the hard plant matter typical for this habitat and preferentially feed on dead material lying on the ground [28, 30], where a mechanism increasing shell adhesion offers no obvious advantage to its bearer. This interpretation is supported by the fact that phylogenetically distantly related haired species, such as Helicodonta obvoluta and Isognomostoma isognomostoma, are found in the same habitats and have in general similar life-styles [31]. However, as long as the positive effect of increased adherence to food plants on the individual fitness is not proven, this remains a hypothesis and does not preclude additional or even other adaptive functions of haired shells.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-5-59/MediaObjects/12862_2005_Article_171_Fig6_HTML.jpg
Figure 6

Active T. villosus foraging on a leaf. Note that the water-film on the leaf is adhering to the shell.

Conclusion

The present comparative analysis suggested that hairs on the shell confer a selective advantage in humid habitats only and are thus lost in drier habitats. In other words, the variability of hairiness within the genus Trochulus could be explained in terms the loss of its adaptive function in a selectively different environment.

Methods

Taxon sampling

Analyses were undertaken on twelve of the about 15 currently recognised species presumed to belong to the genus Trochulus s. str. Chemnitz, 1786 (Hygromiidae, Stylommatophora). However, the exact number of existing species is not known, because the species limits of the widely distributed T. hispidus and T. sericeus on the one hand and T. plebeius and T. striolatus on the other are equivocal [19, 20], the validity of several described taxa is disputed [18, 32] and newly discovered species are not yet formally described (Pfenninger, unpublished data). Since initial analyses showed the existence of cryptic lineages, several populations for each of the putative species were sampled (Table 1). Four species from other genera of the subfamily Hygromiinae and two species of the family Helicidae were used as potential outgroups [33] (GenBank accession numbers AY546263, AY546343, AY546303, AY546284, AY546364, AY546324, AY546283, AY546363, AY546323, AY546291, AY546371, AY546331).

DNA sequencing, lineage identification and phylogenetic analysis

Entire snails were crushed and vortexed in 10% w/v laundry detergent solution for storage at room temperature and tissue digestion [34]. For 78 individuals, a 512 bp segment of the cytochrome oxidase subunit I gene (COI) was amplified with PCR and sequenced. For selected individuals representing the major evolutionary lineages inferred in the previous analysis, a 362 bp fragment of the large subunit mitochondrial ribosomal gene (16S) and 509 bp of the internal transcribed spacer 1 (ITS-1) from the nuclear ribosomal cluster were additionally amplified and sequenced. An amount of 0.2 to 1 ng total DNA (quantified on a 1% agarose gel using a λ Hind III marker) were used as template in polymerase chain reaction (PCR). Specific PCRs were performed with the primers, amplification conditions and temperature profiles shown in Table 2. Primers were used for both specific PCR and subsequential automated direct sequencing. PCR products were purified using E.N.Z.A. Cycle Pure Kit (peqlab, Erlangen, Germany). Ten ng per sample were subjected to cycle sequencing using the ABI Prism Big Dye terminator kit (Perkin-Elmer, Norwalk, CT, USA). Sequencing reactions were electrophoresed on an ABI 377 automated DNA sequencer. In order to verify the results, gene products were sequenced in both directions and the two strands were aligned with SEQUENCE NAVIGATOR 1.0.1 (Perkin-Elmer, Norwalk, CT, USA). Sequences were deposited in GenBank under accession numbers DQ217794-DQ217831. The orthologous DNA sequences were initially aligned using the default settings of CLUSTAL X [35] and optimised by eye. The most likely models of sequence evolution and their parameters according to the Akaike information criterion were inferred for each DNA data partition using MODELTEST v. 3.4 [36]. In an initial analysis, we used the COI data set to identify evolutionary lineages. A 99.9% credible set of phylogenetic trees was estimated with the program MRBAYES [37] by sampling the tree space using a Metropolis coupled Monte Carlo Markov chain, implementing a TN+I+Γ model of COI sequence evolution (where TN denotes Tamura-Nei, Γ is the shape parameter of the gamma distribution and I the proportion of invariant sites). Initial runs as well as a posterior inspection of the likelihoods in the final run showed that a burn-in phase of 10,000 generations was largely sufficient for both analyses to allow the likelihood values to reach convergence. The chain was run for 10,000,000 generations and sampled every 100th generation. An unrooted majority consensus tree was computed from the sampled trees, excluding the trees sampled in the burn-in phase. The procedure was repeated for the phylogenetic data set where the Markov chain was run with separate models of sequence evolution for each data partition (GTR (general time reversible)+I+G for 16S and TVM (transversional model)+ Γ for ITS-1). Outgroup status was assigned to Helixaspersa [33].

Correlation of habitat humidity with shell hairiness

The direct estimation of humidity levels for sampling sites is difficult without long-term observation. However, the precipitation regime, habitat structure and vegetation at a sampling site can give clues on the degree of humidity experienced by the snails. For this behalf, five variables were recorded for all but one population belonging to Trochulus s.str. species. To characterise the microhabitat conditions, the mean light- and humidity indicator values [38] of the three most abundant herbaceous plant species at each sampling site were recorded (variables LIGHTIND and HUMIND). The evaporation regime is strongly influenced locally by the exposure to sun and wind, which was accounted for by characterising each sampling site as either i) entirely shadowed (2), partially or sometimes shadowed (1) and never shadowed (0) (variable SHADOW) and either ii) situated in a closed wood (2), open wood or forest edge (1) or not in a wood (0) (variable WOOD). Ultimately, the humidity conditions of a site depend on the precipitation in the area. As Trochulus species are active mainly during summer, we have recorded the average long-term precipitation from April to September (variable SUMMERPREC). This information was extracted from the climate layers with a spatial resolution of 0.5 min implemented in the computer program DIVA-GIS version 4.2 for the spatial analysis of biodiversity [39]. The variables were summarised in a principal component analysis (PCA).
Table 3

Primers used for specific PCR and direct sequencing, amplification conditions and temperature profiles.

Primer

Sequence

amplification conditions

temperature profile

COI universal [43]

5'-GGTCAACAATCATAAAGATATTGG-3' 5'-TAAACTTCAGGGTGACCAAAAAATCA-3'

total volume 25 μl with: 0.17 mM dNTPs 3 mM MgCl2 in 1 × PCR buffer 0.13 μM of each primer 1 unit Taq polymerase (Invitrogen)

1 cycle of 2.5 min at 94°C

40 cycle 30s at 90°C

1 min at 48°C

1 min at 72°C

1 cycle of 10 min at 72°C

16S universal [44]

5'-CGGCCGCCTGTTT ATCAAAAACAT-3' 5'-GGAGCTCCGGTTTGAACTCAGATC-3'

total volume 15 μl with: 0.1 mM dNTPs 2.5 mM MgCl2 in 1 × PCR buffer 0.2 μM of each primer 0.5 unit Taq polymerase (Invitrogen)

1 cycle of 2.5 min at 90°C

10 cycles of 50s at 92°C

30s at 44°C

40s at 72°C

36 cycles of 30s at 92°C

40s at 48°C

40s at 72°C

1 cycle of 3 min at 72°C

ITS-1 mollusc specific [45]

5'-TAACAAGGTTTCCGTAGGTGAA-3' 5'GCTGCGTTCTTCATCGATGC-3'

total volume 15 μl with: 0.3 mM dNTPs 2.5 mM MgCl2 in 1 × PCR buffer 0.18 μM of each primer 0.5 unit Taq polymerase (Invitrogen)

1 cycle of 3 min at 94°C

40 cycles of 30s at 92°C

30s at 52°C

1 min at 72°C

1 cycle of 5 min at 72°C

For all Trochulus s.str. populations investigated, the presence or absence of hairs on the shell of at least 10 adult individuals was recorded. As the hairs may wear off during adulthood (although rarely completely), the lack of the typical hair pits in the fine sculpture of the shell was taken as evidence for their principal absence. The presence or absence of hairs of the respective populations was then plotted on the PCA ordination.

Bayesian estimation of ancestral character states

In a first approach, we derived the posterior probability distribution of ancestral character states and their rate of change from 3000 trees sampled at random from the 99.9% credibility set of phylogenetic trees, using the Bayesian approach as implemented in the program MULTISTATEBAYES [40]. Applying an uninformative (uniform) prior on the rate parameter distribution, a Markov chain was run for 1,000,000 generations after it reached convergence. The estimated rate parameter ratio for both directions of character change as well as the reconstructed ancestral states for each internal node of the tree investigated was sampled every 200th generation. This procedure estimates i) the probability that the ancestral node existed in the first place and ii) the probabilities of both character states at the respective node. These three probabilities sum up to 1, thus simultaneously taking phylogenetic and character mapping uncertainty into account. In a second approach, the most parsimonious number of character state changes was reconstructed for each of the 99.9% credibility set of phylogenetic trees using the ANCESTRAL STATE RECONSTRUCTION module in MESQUITE [41]. The different reconstructions were then weighted according to the posterior probability of the corresponding tree [42].

Adhesion experiments

The minimum force necessary to move Trochulus shells (upwards oriented apex) with or without hairs over dry and wet, horizontal leaf surfaces was measured. For this behalf, we have chosen the largest species, T. villosus. It would have been desirable to use shells of other lineages as well, however, it was not possible to measure the force necessary to move smaller shells with the necessary accuracy. Twelve T. villosus shells were glued to thin nylon strings. The strings were led over a roll with a small aluminium basket fastened on the other end. Small weights were incrementally added to the basket until the shell began to slide. This was replicated five times for each shell on both water film covered and on dry surfaces. Then, the hairs were mechanically removed to obtain smooth shells and the procedure was repeated. For each condition, differences in minimum force needed to move the shells with or without hairs were tested for significance with an ANOVA design.

Declarations

Acknowledgements

We thank Holger Geupel for technical assistance. MP acknowledges the financial support of the A. Messer Stiftung. Jacques Hausser, Aris Parmakelis, Christian Albrecht and an anonymous referee gave valuable comments on the manuscript. We thank Margret Gosteli and the Museum of Natural History in Bern for samples and additional information.

Authors’ Affiliations

(1)
Abteilung Ökologie & Evolution, J.W. Goethe-Universität
(2)
Deparment of Zoology, Charles University
(3)
Department of Biology, University of Konstanz
(4)
Département d'Ecologie et Evolution, Université de Lausanne

References

  1. Harvey P, Pagel MD: The comparative method in evolutionary biology. 1991, Oxford: Oxford University Press
  2. Porter ML, Crandall KA: Lost along the way: the significance of evolution in reverse. Trends Ecol Evol. 2003, 18: 541-547. 10.1016/S0169-5347(03)00244-1.View Article
  3. Vermeij GJ: A natural history of shells. 1995, Princeton: Princeton University Press
  4. Godan D: Molluscs: Their significance for science, medicine, commerce and culture. 1999, Berlin: Parey
  5. Palmer AR: Quantum changes in gastropod shell morphology need not reflect speciation. Evolution. 1985, 39: 699-705.View Article
  6. Hellberg ME, Balch DP, Roy K: Climate-driven range expansion and morphological evolution in a marine gastropod. Science. 2001, 292: 1707-1710. 10.1126/science.1060102.View ArticlePubMed
  7. West KCA: Morphology and behavior of crabs and gastropods from Lake Tanganyika, Africa: Implications for lacustrine predator-prey coevolution. Evolution. 1991, 45: 589-607.View Article
  8. Britton JC: The relationship between position on shore and shell ornamentation in 2 size-dependent morphotypes of Littorina striata, with an estimate of evaporative water-loss in these morphotypes and in Melarhaphe neritoides. Hydrobiologia. 1995, 309: 129-142.View Article
  9. Wilson AB, Glaubrecht M, Meyer A: Ancient lakes as evolutionary reservoirs: evidence from the thalassoid gastropods of Lake Tanganyika. Proc R Soc Lond Ser B-Biol Sci. 2004, 271: 529-536. 10.1098/rspb.2003.2624.View Article
  10. Schilthuizen M: Sexual selection on land snail shell ornamentation: a hypothesis that may explain shell diversity. BioMedCentral Evolutionary Biology. 2003, 3: 13-10.1186/1471-2148-3-13.View Article
  11. Goodfriend GA: Variation in land-snail shell form and size and its causes – a Review. Systematic Zoology. 1986, 35: 204-223.View Article
  12. Saleuddin ASM: Ultrastructural studies on formation of periostracum in Helix aspersa (Mollusca). Calcified Tissue Research. 1976, 22: 49-65. 10.1007/BF02010346.View ArticlePubMed
  13. Wade CM, Mordan PB, Clarke B: A phylogeny of the land snails (Gastropoda: Pulmonata). Proc R Soc Lond B. 2000, 268: 413-422.View Article
  14. Suvorov AN: Some mechanisms of adaptation to the wet microhabitats in higher geophila (Mollusca, Pulmonata). Zhurnal Obshchei Biologii. 1999, 60: 177-188.
  15. Solem A: The shell makers. 1974, New York: John Wiley & Sons
  16. Ronquist F: Bayesian inference of character evolution. Trends Ecol Evol. 2004, 19: 475-481. 10.1016/j.tree.2004.07.002.View ArticlePubMed
  17. Baxter ML: The promise of a DNA taxonomy. Philosophical Transactions of the Royal Society of London – Series B: Biological Sciences. 2004, 359: 669-679. 10.1098/rstb.2003.1447.View Article
  18. Falkner G: Trichia (Trichia) graminicola n.sp. aus Südbaden (Gastropoda: Helicidae). Archiv für Molluskenkunde. 1973, 103: 209-227.
  19. Forcart L: New researches on Trichia hispida (Linnaeus) and related forms. Proceedings of the first European Malacological Congress. 1962, 79-93.
  20. Perrin N, Honsberger P, Pontet A: Approche biochimique et morphologique des espèces du genre Trichia (Helicidae, Gastropoda) de Suisse occidentale. Revue Suisse de Zoologie. 1984, 92: 483-495.View Article
  21. Falkner G: Beiträge zur Nomenklatur der europäischen Binnenmollusken, VIII: Nomenklaturnotizen zu europäischen Hygromiidae (Gastropoda: Stylommatophora). Heldia. 1995, 2: 97-107.
  22. Thomaz D, Guiller A, Clarke B: Extreme divergence of mitochondrial DNA within species of pulmonate land snails. Proc R Soc Lond Ser B-Biol Sci. 1996, 263: 363-368.View Article
  23. Pinceel J, Jordaens K, Pfenninger M, Backeljau T: Rangewide phylogeography of a terrestrial slug in Europe: evidence for Alpine refugia and rapid colonization after the Pleistocene glaciations. Mol Ecol. 2005, 14: 1133-1150. 10.1111/j.1365-294X.2005.02479.x.View ArticlePubMed
  24. Pfenninger M, Magnin F: Phenotypic evolution and hidden speciation in Candidula unifasciata ssp (Helicellinae, Gastropoda) inferred by 16S variation and quantitative shell traits. Mol Ecol. 2001, 10: 2541-2554. 10.1046/j.0962-1083.2001.01389.x.View ArticlePubMed
  25. Wilke T, Duncan N: Phylogeographical patterns in the American Pacific Northwest: lessons from the arionid slug Prophysaon coeruleum. Mol Ecol. 2004, 13: 2303-2315. 10.1111/j.1365-294X.2004.02234.x.View ArticlePubMed
  26. Schileyko AA: On the systematics of Trichia s. lat. (Pulmonata: Helicoidea: Hygromiidae). Malacologia. 1978, 17: 1-56.
  27. Cook A: Behavioural Ecology: doing the right thing, in the right place at the right time. The Biology of Terrestrial Molluscs. Edited by: Barker GM. 2001, Oxon: CABI Publishers, 447-488.View Article
  28. Frömming E: Biologie der Mitteleuropäischen Landgastropoden. 1954, Berlin: Dunker & Humbold;
  29. Denny M: Locomotion: the cost of gastropod crawling. Science. 1980, 208: 1288-1290.View ArticlePubMed
  30. Kerney MP, Cameron RAD, Jungbluth JH: Die Landschnecken Nord- und Mitteleuropas. 1983, Hamburg; Berlin: Paul Parey
  31. Speiser B: Food and feeding behaviour. The Biology of Terrestrial Molluscs. Edited by: Barker GM. 2001, Oxon: CABI Publishers, 259-288.View Article
  32. Gittenberger E, Neuteboom WH: On Trichia alpicola (Eder, 1921) from Switzerland (Mollusca: Gastropoda Pulmonata:Hygromiidae) and the spiral sculpture on its shell. Zoologische Mededelingen. 1991, 65: 247-250.
  33. Steinke D, Albrecht C, Pfenninger M: Molecular phylogeny and character evolution in the Western Palaearctic Helicidea s.l. (Gastropoda, Stylommatophora). Molecular Phylogenetics & Evolution. 2004, 33: 724-734. 10.1016/j.ympev.2004.03.004.View Article
  34. Bahl A, Pfenninger M: A rapid method of DNA isolation using laundry detergent. Nucleic Acids Res. 1996, 24: 1587-1588. 10.1093/nar/24.8.1587.PubMed CentralView ArticlePubMed
  35. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMed
  36. Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMed
  37. Huelsenbeck JP, Ronquist F: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMed
  38. Ellenberg H, Weber HE, Düll R, Wirth V, Werner W, Paulißen D: Zeigerwerte von Pflanzen in Mitteleuropa. 1992, Göttingen: Goltze KG, 2
  39. Hijmans RJ, Guarino L, Cruz M, Rojas E: Computer tools for spatial analysis of plant genetic resource data: 1. DIVA-GIS. Plant Genetic Resources Newsletter. 2001, 127: 15-19.
  40. Pagel M, Meade A, Barker DM: Bayesian estimation of ancestral character states on phylogenies. Syst Biol. 2004, 53: 673-684. 10.1080/10635150490522232.View ArticlePubMed
  41. Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 1.05. 2004
  42. Huelsenbeck JP, Rannala B, Larget B: A Bayesian framework for the analysis of cospeciation. Evolution. 2000, 54: 352-364.View ArticlePubMed
  43. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotech. 1994, 3: 294-299.
  44. Palumbi S, Martin A, Romano S, McMillan W, Stine O, Grabowski G: The simple fool's guide to PCR version 2.0. 1991, Honolulu: University of Hawaii
  45. Armbruster GFJ, van Moorsel CHM, Gittenberger E: Conserved sequence patterns in the non-coding ribosomal ITS-1 of distantly related snail taxa. J Molluscan Stud. 2000, 66: 570-573. 10.1093/mollus/66.4.570.View Article

Copyright

© Pfenninger et al; licensee BioMed Central Ltd. 2005

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement