Living specimens were collected from several localities in Tunisia during two periods: spring 2014, and winter 2015/2016. Geographic coordinates were recorded using GPS (see Table 1). For subsequent molecular analysis, specimens were preserved and stored in 80% ethanol until dissection and DNA extraction.

Empty shells were also collected (see section material under the species description) in order to complete the distributional record of the species. Specimens used in this study (both shells and preserved animals) are housed in the voucher collections of the High Institute of Biotechnology of Monastir and the Natural History Museum Bern; all sequenced specimens are deposited in the museum’s collection.

First assessments of the shell morphological characters were done by using simple magnifying glasses. Preserved animals were dissected under LEICA M212 stereo-microscope using thin tweezers. The genital organs of the specimens were removed from the body, the genital situs (i.e. the outer morphology of the complete hermaphroditic genital organ) and further morphological details were investigated. After that, shells, genital situs, and details of the genital organs were photographed with a LEICA DFC 425 camera combined with a LEICA M205 C. The multifocal images were processed by using an imaging software (Imagic Switzerland).

MVHN Museu Valencià d’Historia Natural;

MHNG-MOLLMuseum d’Histoire Naturelle de Genève, malacological collection;

NMBENaturhistorisches Museum der Burgergemeinde Bern;

ZMHZoological Museum of the University of Hamburg.

D: shell diameter; H: shell height; PD: peristome diameter; PH: peristome height; W: number of whorls.

Fourteen specimens of Xerocrassa from southern Tunisia could be used in this study, originating from five localities. Sequenced specimens are housed in the voucher collection of the NMBE (Table 2). In the analysis, sequences of four Cretan Xerocrassa species were also included (Sauer and Hausdorf 2009), and eleven Spanish and Balearic Xerocrassa species from the work recently published by Chueca et al. (2017).

As outgroup species Cernuella virgata, Trochoidea elegans, and Trochoidea pyramidata were used. All three species are each represented by one specimen from Tunisian localities, and complemented by one specimen of Hygromia limbata, one Xerosecta adolfi, and one T. elegans (Razkin et al. 2014). All specimens used to produce phylogenetic trees are listed in Table 2. Specimens where nuclear markers are not available were excluded from the analysis of the concatenated mitochondrial // nuclear dataset. Thus, all Cretan Xerocrassa specimens, except two specimens of X. cretica (recently collected by Neubert), and the Tunisian Trochoidea and Cernuella species were not used in this type of analysis.

Total genomic DNA was extracted from the foot muscle tissue using the standard phenol chloroform method (Estoup et al. 1996). Two mitochondrial gene fragments and one rDNA region were chosen to be analysed in the current study. Mitochondrial markers were consisting of Cytochrome c oxidase subunit I (COI) and the 16S ribosomal RNA subunit (16S) gene. The nuclear marker was formed by the 3’ end of the 5.8s ribosomal RNA, the complete ITS2 region and the 5’end of the large subunit of the 28S rRNA. Polymerase chain reactions (PCR) were performed in a reaction mixture, containing 15 ng of DNA template, 1×1.5 mM buffer reaction, 0.1 mM of each selected couple primers, 0.2 mMdNTPs, Taq polymerase (1.25U) and adjusted till a total volume of 25 µl with DNAase free water/sterilized water (UNIMED) (H2O). PCR reactions were run under following conditions: 3 min at 95°C, followed by 35cycles of 1 min at 95°C, 1 min at 40°C and 1 min at 72°C and finally, 5 min at 72°C for COI. For 16S the amplification conditions were: 3 min at 95°C, followed 35 cycles of 1 min at 95°C, 1 min at 50°C and 1 min at 72°C. To amplify the ribosomal cluster, two pairs of primers were used to get a sequence of 1300 bp: the standard LSU1/LSU3 and the 28SF/28SR (see Table 3). PCR reactions were run under the following conditions: 3 min at 96°C, followed by 35 cycles of 1 min at 94°C, 1 min at 50°C and 1 min at 72°C and finally, 5 min at 72°C for LSU1/LSU3 and 5 min at 95°C, followed by 35cycles of 1 min at 95°C, 30 s at 62°C and 1 min at 72°C and finally, 10 min at 72°C for 28SF/28SR. PCR products were sequenced using automated and standardised ABI 3730 XL sequencing run with a read length up to 1100 bp (PHRED20 quality) and using the same primers as for the PCR (Table 3).

Forward and reverse sequences were assembled, checked for ambiguities and aligned using default settings of “Clustal W” implemented in Bioedit V 7.2.5 (Hall 1999). Aligned sequences of Tunisian Xerocrassa species were analysed using DnaSP v5.10.01 software (Librado and Rozas 2009) to estimate number of informative sites and nucleotide diversity for each marker used. The p-distance values within Tunisian samples were estimated using Mega v.6 (Tamura et al. 2013). The relationships of inferred haplotypes of mitochondrial nuclear and concatenated dataset of Tunisian Xerocrassa species were estimated using the TCS method (Clement et al. 2002) implemented with Popart software v1.7 (Leigh et al. 2015).

Our data consist of two mitochondrial markers and one nuclear ribosomal cluster. The data was partitioned used the PartitionFinder software v1.1.1 (Lanfear et al. 2012), in six partitions: three codon positions of the COI, the 16S the rRNA 5.8S and 28S were considered as a single partition and finally the ITS2.

For the mitochondrial dataset as well as for the concatenated data, we produced two phylogenetic trees within the Mediterranean Xerocrassa species using the Maximum Likelihood (ML) and the Bayesian inference (BI). The ML analyses were conducted using RAxML v7.2.6 (Boc et al. 2012, Stamatakis 2006) under the GTRGAMMA model, with 1000 nonparametric bootstrap replicates to estimate node support. For the Bayesian Inference, we used Mr Bayes v3.2.2 (Ronquist and Huelsenbeck 2003) using partition scheme and substitutions models suggested by PartitionFinder v1.1.1 (Lanfear et al. 2012). Four independent runs were conducted for 10⁶ generations, sampling every 1000. The first 25% trees were discarded as default burn-in and a majority rule consensus tree was calculated from the remaining trees. The topology obtained, and the posterior probabilities of each node were displayed on Figtree V1.4.0 (Rambaut 2012).

Both, the results of our morphological research on the genital organs as well as the molecular study, prove the affiliation of Helix latastei and Helix latasteopsis to the genus Xerocrassa Monterosato, 1892. For the subgeneric placement refer to the chapter “Discussion,”

Xerocrassa (Xerocrassa) latastei (Letourneux in Letourneux & Bourguignat, 1887)

Figs 1, 2, 3

1887 Helix latastei Letourneux in Letourneux & Bourguignat, Prodrome de la malacologie terrestre et fluviatile de la Tunisie: 63 [Ketenna et dans le vallon de l’Oued El-Ftour, ainsi qu’à l’oasis du Hammam de Gabès. Plaine entre Ras-el-Aïn et Sidi-Salem-Bouguerara. Bir-el-Ahmar. Bords de l’Oued Medzesar et de l’Oued Taferma entre Aïn-Magroun et Fratis. Ras-ed-Djerf, vis-à-vis de Djerba; Zarzis, etc. (Let.). — En Algérie: Ouled Naïl près de Biskraoù, à Aïn-Gussera, à Bou-Ghezoul sur les hauts plateaux, entre Boghar et Laghouat et entre cette ville et Djelfa].

Type specimens

Brandt (1959: 113) considered four taxa of hygromiid species described by Letourneux and Bourguignat, 1887 to constitute the species H. latastei. Our investigation of the type specimens of these taxa revealed that the species Helix fratisiana and Helix tafermica, which had been listed by him in the synonymy of H. latastei, belong to species of the Hygromiidae living in Tunisia. In order to stabilize nomenclature, we herewith select MHNG-MOLL 115121 as lectotype for Helix latastei [hic!]. Thus, the type locality of this species is herewith restricted to Ketenna [= Kettana]: mouth of Oued El Ferd, Gouv. Gabés, at 33.7575 10.2047; paralectotypes MHNG-MOLL 115121b/4, MHNG-MOLL 115128/2.

Additional specimens examined

Bou Hedma, 29.3.1997, leg. J. Gugel, 34.4958°N 9.488°E, NMBE 516753/1; Boughrara, Medenine, 6.12.2015, leg. Ezzine, NMBE 541952/3, ditto, NMBE 547176/3; ditto, NMBE 541955/7 (preserved); Jorf (El Djorf), Mednine, 6.12.2015, leg. Ezzine, NMBE 541956/1 (preserved), NMBE 549907/1 (anat.); “plaine entre Ras-el-Aïn et Sidi-Salem-Bouguerara”, MHNG-MOLL 115118/3, MHNG-MOLL 115126/6, MHNG-MOLL 115127/6, MHNG-MOLL 115129/4; Bir-el-Ahmar MHNG-MOLL 115119/1; Zarzis MHNG-MOLL 115120/2; “Oued el Ftour près de Gabès” MHNG-MOLL 115124/6; “Ras-ed-Djerf, vis-à-vis de Djerba” MHNG-MOLL 115125/1. — Specimens recorded from literature: ruins of Gighti close Djorf (Djerba) (Jaeckel 1963).

Diagnosis

Shell small to medium sized, thick, basic colour white; protoconch brownish to blackish; three first whorls with granulations; whorls ribbed; suture moderately deep; umbilicus very small, conical.

Description

Shell small to medium sized, depressed globular, thick, basic colour creamy white; protoconch very small, brownish to blackish, smooth, consisting of 1½ whorls; teleoconch consisting of 5½ slightly flattened whorls, sculptured by moderately sized axial ribs; three first whorls brown with whitish granules; lower teleoconch whorls with up to 5 brown spherical bands; suture moderately deep; underside often white; aperture sub-spherical, slightly descending; columellar peristome thick; umbilicus moderately small, conical.

Genital anatomy

The description of the genital organs is taken from an adult and mature specimen collected in El Djorf. Figure 1B shows the lumen of the atrium with its internal structures.

Figure 1. 

Anatomy of genital organs of Xerocrassa latastei; Jorf, 6.12.2015, leg. Ezzine, NMBE 549907/1; A situs B atrium. Abbreviations: a = atrium; al = atrial lumen; ds = dart sac(s); ep = epiphallus; fl = flagellum; gm = glandulae mucosae; itc = internal tissue cone; mrp = penial retractor muscle; p = penis; ped = pedunculus; po = pore of penial papilla; pp = penial papilla; sod = spermoviduct; sti = stimulator (?); vd = vas deferens.

Male part. Penis club-shaped, thick; epiphallus longer than penis; penial retractor muscle inserting at the boundary between penis and epiphallus, with a strong fascia enveloping the genitals; flagellum short; penial papilla subdivided in a simple basal shaft and a subsequent part characterised by deep perpendicular grooves, terminal part of the penial papilla strongly kinked, with central pore at its tip.

Genital atrium. Considerably thickened, lumen filled by two structures: 1) a strong crest of fleshy tissue (here called stimulator), auricle-shaped, the interior side (i.e. opposite to the penial papilla) with zigzag-shaped longitudinal pilasters becoming smooth when entering the interior wall of the atrium, and 2) a longitudinal spoon- or tongue-shaped tissue plate (here called internal tissue cone), with the outer rims bent upwards forming a hollow structure.

Female part. Two very small, almost spherical dart sacs in opposite position; glandulae mucosae simple, tubes randomly attached on the vaginal wall between dart sacs and pedunculus; vagina moderately long, pedunculus formed by a quite strong tube.

Measurements. Lectotype latastei: D: 15.9 mm; H: 12.39 mm; PD: 8.58 mm; PH: 6.72; W: 6.25.

Distribution

(Fig. 2). This species is currently known from the coastal and neighbouring inland areas of central and southern Tunisia. It occurs almost in sympatry with H. latasteopsis in some areas of the province Medenine and Sidi Bouzid.

The Senckenberg Museum houses a considerable number of dry shells under the name H. latatstei from Libya, based on the collections of Brandt (1959: 112 ff.). They were examined by Neubert during the last years, and they in fact are very similar to X. latastei from Tunisia. However, all these shells were collected in the Cyrenaica and its hinterland with the westernmost locality being Marsa Brega (ca. 200 km SSW of Bengasi). So far we have not seen any shells from the Sirte nor the Tripolitanian area towards Tunisia, which embraces almost half of the coastal stripe of Libya. The gap to the Tunisian populations is more than 800 km as the crow flies. This area was visited several times by Kaltenbach (Kaltenbach 1950a; 1950b), but there are no records for X. latastei from this area in his rich collection, which is also housed in SMF. As long as no preserved specimens from the Cyrenaica are available, we consider these populations as not conspecific.

Figure 2. 

Distribution of Xerocrassa latastei, Xerocrassa latasteopsis and Xerocrassa lacertara.

Remark

Specimens of this species are characterized by a globose shell with a quite small umbilicus if compared to the large Cernuella species, which live sympatrically in southern Tunisia.

The internal structures in the genital atrium are poorly understood. However, when dissecting the atrium, the internal tissue cone is always found to almost completely envelop the penial papilla; the situation shown in Fig. 1B is the result of pulling the penial papilla out of the internal tissue cone. Spreading the opened atrium then leads to a position of this organ on the right side.

Figure 3. 

Anatomy of genital organs of Xerocrassa latasteopsis; A X. latasteopsis, NMBE 551301, Henchir el Zitouna: situs, penial lumen and atrial lumen B X. latasteopsisNMBE 549906, Sidi Aich 2, situs, penial lumen; and atrial lumen.

This dataset consists of two mitochondrial markers (COI and 16S) and one nuclear cluster (5.8S-ITS2-28S). The mitochondrial data was analysed first and afterwards the nuclear marker was added to confirm the results.

The results of the anatomical and morphological studies of the Tunisian samples show that there are two Xerocrassa species existing: X. latasteopsis and X. latastei. The nucleotide divergence of these two morphological groups is studied, and a haplotype network is produced. Among fourteen sequences of 1090 bp (655 bp of COI and 435 bp of 16S) of Tunisian Xerocrassa, six haplotypes were identified using both markers, suggesting a high haplotype diversity (Hd=0.8022). The haplotypes obtained cluster in two divergent haplo-groups: the first is formed by samples collected from SidiAich1, SidiAich2, and Henchir El Zitouna, and the second was formed by samples collected from El Djorf and Boughrara (Fig. 8A). A high number of variable sites could be found in-between the groups (190 sites: COI: 118 and 16S: 72), but only a low number within the groups (maximum of 60 sites within the group of Boughrara_El Djorf). Thus, the nucleotide divergence was of 17% between haplo-groups, and 5.5% within the haplo-group of Boughrara-El Djorf and a low divergence within Gafsa-Henchir El Zitouna (<1%).

The analyses of each mitochondrial separately showed some differences between COI and 16S. The nucleotide divergence of COI sequences reached 18% between groups, and varied between 0.4% and 6% within haplo-groups. Additionally, the amino acid composition of the partial COI sequence (218 amino acids) displayed eight different amino acids between haplo-groups of which two are of different polarity.

The ribosomal gene 16S showed a high nucleotide divergence between haplo-groups (16%) and a low divergence with a maximum of 2% of nucleotide divergence within the haplo-group of Boughrara-El Djorf and a monomorphic haplo-group formed by the sequences of Gafsa and Henchir el Zitouna.

The assessment of the nuclear ribosomal cluster 58S-ITS2-28S (1320 bp) showed that all 5.8S and 28S sequences used were identical. The sequences of ITS2 displayed only a single insertion/deletion mutation and one substitution between X. latasteopsis and X. latastei suggesting an extreme conservation of nuclear sequences in Tunisian Xerocrassa species. Adding the ribosomal cluster did not affect the haplotype diversity obtained using the mitochondrial data. In fact, we observed six haplotypes grouped in two haplo-groups with a haplotype diversity of 0.8022, a nucleotide divergence of 8% between X. latasteopsis and X. latastei (Fig. 8B).

Both topologies of the mitochondrial (mt) data from the ML and BI analyses are identical. The tree obtained is rooted by two Trochoidea species, C. virgata, Xerosecta adolfi and H. limbata. Mediterranean Xerocrassa species were divided in two groups following the geographical distribution pattern (Fig. 9A): 1) An East-Mediterranean group formed by the Tunisian and the Cretan Xerocrassa species. 2) A West-Mediterranean group composed by three clades: one clade comprising the Balearic radiation, and two continental Spanish clades. Both groups are supported by high bootstrap values (BS=100%) and posterior probability (PP=0.83). The East-Mediterranean group shows three highly supported clades: one formed by Tunisian species, one composed by the Cretan X. cretica and one formed by the rest of the Cretan Xerocrassa radiation. In Tunisia, the Xerocrassa species split into two well supported (BS=100%, PP=1) monophyletic clades. In X. latasteopsis samples, which were collected from two distinct areas, were grouped in rake proving a low mitochondrial diversity within species.

The position of “Xerocrassa meda” close to Trochoidea is quite unexpected. In case it is not a mix-up with a specimen from the highly polymorphic Trochoidea spratti-group, then the mitochondrial sequences are not informative. Inclusion of nuclear markers in the analysis will probably yield a better result.

For the concatenated tree, all Cretan Xerocrassa species except two samples of X. cretica had to be excluded because of nuclear data deficiency. Both trees based on Ml and BI analyses show identical topologies; they are rooted by T. elegans, Xerosecta adolfi and H. limbata. Again, two main Mediterranean clades appear, their node is well supported (Fig. 9B): all Spanish Xerocrassa species cluster together forming a single clade (BS=100%, PP=1), which in itself is divided in the three groups, one insular and two continental. Here, the Cretan X. cretica clusters with Tunisian Xerocrassa species composing a well-supported group (BS=100%, PP=1). The two Tunisian species are well separated and strongly supported (BS=100%, PP=1).

The type species of Xerocrassa Monterosato, 1892 is the east Mediterranean species Helix seetzeni L. Pfeiffer, 1847 (by monotypy). Xerocrassa is currently characterized by a symmetrical dart apparatus consisting of two small accessory sacs (= appendicula sensu auct.) and usually four branched glandulae mucosae around the vagina, irregular folds at the inner side of the wall of the vagina and the lack of a well-developed appendix at the atrium; the penis is innervated from the right cerebral ganglion (Hausdorf and Sauer 2009: 375). The absence of the atrial appendix is basically the only character state that separates Xerocrassa from Trochoidea Brown, 1827. Hausdorf and Sauer (2009) report the presence of an atrial “bulge-like stimulatory structure” in some Xerocrassa species such as X. cretica, X. franciscoi and X. heraclea, which can be seen as a small protuberance at the side of the atrium in the Tunisian X. latastei and X. latasteopsis (see Figs 1, 3), the atrium is much wider, and bulge considerably more pronounced. The homology of this organ with the atrial appendix seen in Trochoidea is not clear, and there is no other evidence than the similar position at the atrium. The internal tissue crest is here called a stimulator referring to the similarly shaped stimulator found in the atrium of many helicid genera. It seems to consist of two parts, a massive block of tissue, and a separate tongue- or cone-shaped stylus. A similar structure is illustrated by Giusti et al. (1995) for Trochoidea species (CAA in their nomenclature). Thus, a morphological separation of the two genera remains difficult, and our results show that the two genera are closely related.

Pallary (1919) based his monotypic genus Ereminella on Helix latastei Letourneux in Letourneux & Bourguignat, 1887 without delivering any discriminating characters. Brandt (1959) recognized that the species recorded by Bisacchi (1932: 361, Figs 2–4) under Helicella (Xerocrassa) pseudosimulata (Germain 1921) from El Agheila and Soluch-Agedabia (Libya) was a misinterpretation, and identified it with H. latastei, which he subsequently affiliated to Trochoidea (Ereminella). In the same publication, Brandt introduced the monotypic subgenus Trochoidea (Xerobarcana) (based on Xerobarcana huggani Brandt, 1959), and Trochoidea (Xeroregima) (based on Trochoidea (Xeroregima) regimaensis Brandt, 1959). Both subgenera show the same principal construction of their genital organs and are congruent with what is considered Xerocrassa today. Xerobarcana is defined as “differing from all other subgenera of this genus [Trochoidea] by the rudimentary wart-like flagellum and the conspicuously strong vas deferens” (translated from the original German text). Today, the relative length of the flagellum is considered a character state that encodes on species level, and is widely used in hygromiid and geomitrid taxonomy (Hausdorf 2000: 62); thus, the reduced flagellum reported by Brandt simply constitutes a character state of Xerocrassa huggani. The definition of Xeroregima is as follows: “Anatomically differentiated from Trochoidea (Trochoidea) s. str. by the lack of the vaginal appendix [sic!] and the penis, which is club-like swollen at the transition between penis and epiphallus” translated from the original German text). Apparently, Brandt confused the terms vaginal with atrial, and thus exactly described the situation as known in Xerocrassa! Even the club-shaped transition between penis and epiphallus is perfectly seen in the majority of Cretan Xerocrassa species as well as in X. latastei. Consequently, we consider Xeroregima as a junior synonym of Xerocrassa.

Figures 4–7. 

Xerocrassa species. 4 Xerocrassa latastei, lectotype MHNG-MOLL 115121, Ketenna [= Kettana], D = 15.9 mm 5 Xerocrassa latasteopsis, lectotype MHNG-MOLL 115131, Foum Hallouf, D = 18.34 mm 6 Xerocrassa latasteopsis, NMBE 549906, Sidi Aich 2, D = 18.2 mm 7 Xerocrassa lacertara, syntype MHNG-MOLL 114001, “Plaines entre Djelfa et El-Aghouat”, D = 11.8 mm.

As with the morphological and anatomical investigation, the results of our molecular approach show that, independently which maker is considered, Tunisian samples divided into two species and cluster together with the selected Xerocrassa species from Crete, the Balearic Islands and Spain, and thus our generic affiliation of the species is correct. There are several remarkable findings, which require deeper examination.

The divergence of the COI sequences between Tunisian species (18%) widely exceeded the threshold of 3% as suggested by Hebert et al. (2003) to characterize animal species in general and the threshold of 4% to identify land snails (Davison et al. 2009). In Tunisia, Chott el Jerid is widely described as a geographical barrier for many taxa (Millington et al. 1989; Ben Othmen et al. 2009; Abdallah et al. 2012; Farjallah et al. 2012). Such a barrier may restrict the gene flow between geographically isolated populations resulting in independent evolution and increase the genetic divergence within species (Funk 2003). In this case, X. latasteopsis shows a low divergence between the northern (Sidi Aich 1, Sidi Aich 2) and the southern (Henchir el Zitouna) populations, which share one haplotype. This result disproves this hypothesis for the snail species concerned, suggesting that Chott el Jerid does not restrict the gene flow. It cannot be considered as a barrier for this species. In contrast, X. latastei shows a quite high divergence within the population of Boughrara (6%) which could be interpreted as individual diversity.

This high divergence between the two Tunisian Xerocrassa species (16%) was also demonstrated by analysis of the 16S marker. High values of genetic divergence were reported for the land slug Phyllocaulis (13.1%) (Gomes et al. 2010) and between congeneric species of Ariophantidae and Dyakiidae (4.3 to 10.1%) (Abu-Bakar et al. 2014). Moreover, Liew et al. (2009) reported divergence values of 5% to 25% within Everettia spp. (Dyakiidae). The divergence of the 16S between Tunisian Xerocrassa species is higher than the divergence of 11.8% between Cretan Xerocrassa species as shown by Sauer and Hausdorf (2009). The divergence seen here is quite remarkable but not completely outstanding.

The nuclear cluster 5.8S-ITS2-28S widely confirms the results obtained from the mitochondrial markers. Both, the 5.8S and 28S sequences seem to be conserved within the Tunisian species, and the ITS2 shows only one nucleotide substitution and one insertion/deletion mutation. Thus, the genetic variability is focused on the mitochondrial markers, while the nuclear markers investigated seem to be highly conserved.

Figure 8. 

Haplotype Network showed the number of variable sites. A Haplotypes and numbers of variables sites based on mitochondrial markers B Haplotypes numbers of variables sites based on concatenated mitochondrial and nuclear data.

This is the first time that a combined phylogeny for this widespread genus has been shown. As could be expected, the clades follow the distribution pattern of a west and an east Mediterranean group. Each cluster includes an island radiation and a continental radiation. The latter fall in two groups for Europe, and one of them, which includes Helix montserratensis Hidalgo, 1870 as its type species [by monotypy] may bear the subgeneric name Amandana Fagot, 1891. The results of these combined markers proved the results obtained by mitochondrial markers and confirmed the split between geographical Xerocrassa groups. Our results suggest that Tunisian Xerocrassa species are more closely related to the Cretan species than to the Spanish and Balearic species. However, within the east Mediterranean clade, the relationship between X. cretica and the rest of the Cretan radiation is not that close with a low support of 0.75 in the mitochondrial tree. A direct comparison with the Tunisian species is problematic. The eastern Mediterranean area, especially Libya, and Egypt, is heavily undersampled, and including more species from this area and the Middle East will certainly change the relative position of Tunisian species to the Cretan species as well as the position of X. cretica on the tree.

The shell morphology of land snails is extremely affected by environmental conditions (Alonso et al. 1985; Chiba 1999; Pfenninger et al. 2006). The use of these characters in the taxonomic analysis of land snail species were severely criticized (Giusti and Manganelli 1992; Schilthuizen and Gittenberger 1996; Uit de Weerd et al. 2004; Holland and Hadfield 2007). As shown by (Elejalde et al. 2008), a comparison between shell morphological and molecular characters result in incongruent data. However, the integrative approach as used here results in a distinct network of character states enabling to interpret the morphology even of the shells, and to formulate distinctive shell traits. Here, the morphological and anatomical disparity between X. latasteopsis and X. latastei has been confirmed by phylogenetic analysis. In the reverse conclusion we now can pinpoint the significance of relative flagellum length and form of penial papilla as well as the ribbing mode of the shell, its colouration and other structural details as relevant and useful for species identification within Tunisian Xerocrassa species (Hausdorf and Sauer 2009).

Figure 9. 

Phylogenetic trees obtained by Bayesian inference (BI) and Maximum Likelihood (ML) methods. A Tree inferred based on partial mitochondrial sequences of COI and 16SB Tree inferred based on mitochondrial data, partial sequences of 5.8S, complete sequence of ITS2 and partial sequence of 28S rRNA. Posterior probability (PB) obtained from Bayesian analysis and bootstrap values obtained from Maximum likelihood (ML) were presented on each node (*: BS= 100).

Figure 9. [Continued]