Fine-scale genetic structure, phylogeny and systematics of threatened crayfish species complex

Fine-scale genetic structure, phylogeny and systematics of threatened crayfish species complex

Molecular Phylogenetics and Evolution 61 (2011) 1–11 Contents lists available at ScienceDirect Molecular Phylogenetics and Evolution journal homepag...

835KB Sizes 3 Downloads 80 Views

Molecular Phylogenetics and Evolution 61 (2011) 1–11

Contents lists available at ScienceDirect

Molecular Phylogenetics and Evolution journal homepage: www.elsevier.com/locate/ympev

Fine-scale genetic structure, phylogeny and systematics of threatened crayfish species complex S. Chiesa a,b, M. Scalici a, R. Negrini c, G. Gibertini a, F. Nonnis Marzano b,⇑ a

Department of Biology, University of ‘‘Roma Tre’’, V. le G. Marconi 446, 00146 Rome, Italy Department of Evolutionary and Functional Biology, University of Parma, V. le Usberti 11, 43100 Parma, Italy c Institute of Zootechnical Sciences, Catholic University of Sacred Heart, Via Emilia Parmense 84, 29100 Piacenza, Italy b

a r t i c l e

i n f o

Article history: Received 25 June 2010 Revised 15 March 2011 Accepted 30 March 2011 Available online 19 April 2011 Keywords: Molecular evolution Phylogenetics Austropotamobius pallipes AFLP COI

a b s t r a c t Systematic uncertainties in the crayfish Austropotamobius pallipes are well grounded by the number of species and subspecies described using different approaches, causing scientists to define this taxon as ‘‘complex’’. However, a key task that conservation programmes are facing regarding the recent and drastic decline of European populations, is the coherent systematic classification of this threatened species. Here we present results obtained by coupling mtDNA and genome analysis suggestive of a novel evolutionary framework to explain the relationships among phylogenetic lineages of A. pallipes. The direct sequencing of mtDNA COI gene fragment revealed a strong geographic structure with four distinct haplogroups separated by a range of 5–25 mutations. However, mitochondrial data were not supported by genomic fingerprinting based on 535 AFLP polymorphisms. Nuclear markers showed an unexpected moderate level of genetic differentiation and the absence of any geographic structure. Consequently, this study proposes that the taxonomic hypothesis of a single species of A. pallipes settling the Italian continental waters, is affected by complex evolutionary events. To solve the paradox, we hypothesized an evolutive scenario in which the separation of ancient mtDNA lineages likely occurred before the latest glacial periods. However, the speciation process remained incomplete due to secondary intensive postglacial contacts that forced the mingling of the genomes, and confounds the phylogeographic signature still detectable within mtDNA. Postglacial dispersion and the following demographic events, such as founder effects, drift and bottlenecks, abruptly depleted the local mtDNA variation, and shaped the current genetic population structure of white-clawed crayfish. Ó 2011 Elsevier Inc. All rights reserved.

1. Introduction In the last decade a number of studies investigated the taxonomy and evolutionary relationships among phylogenetic lineages of crayfish species (Fezner and Crandall, 2003; Buhay and Crandall, 2005; Song et al., 2008). Although based on several different molecular approaches, the proposed taxonomic and systematic models do not appropriately fit with the observed genetic variability. Controversial data induced the scientists to describe a considerably high number of different ‘groups’ or lineages and to temporarily define crayfish species as a ‘‘complex’’ (Austin, 1996; Santucci et al., 1997; Mathews et al., 2008) or cryptic species. This means that the evolutionary and biogeographic histories of these decapods are far from being completely understood (Crandall

⇑ Corresponding author. Fax: +39 0521 905657. E-mail addresses: [email protected] (S. Chiesa), [email protected] (M. Scalici), [email protected] (R. Negrini), [email protected] (G. Gibertini), [email protected] (F. Nonnis Marzano). 1055-7903/$ - see front matter Ó 2011 Elsevier Inc. All rights reserved. doi:10.1016/j.ympev.2011.03.031

et al., 2000; Scholtz, 2002; Fezner and Crandall, 2003; Buhay and Crandall, 2005). Phylogeny and adaptive radiation of the family Astacidae in European regions have been interpreted as the result of past limited speciation events. Only five native European species were described up to now (Holdich, 2002; Machino and Holdich, 2006), as a result of relict taxa. The systematic uncertainty also suggests that elucidating phylogeographic relationships in this group is a pressing issue, especially regarding the implications for its conservation (Mathews et al., 2008). In this context, the white-clawed crayfish Austropotamobius pallipes (Lereboullet, 1858) represents a great example of such problems. It is the most common native river dwelling crayfish in Italy, once widely distributed in Western Europe (Souty Grosset et al., 2006), however nowadays it is suffering a progressive and drastic decline of both biomass and number of populations due to ecological and environmental factors (Gherardi and Holdich, 1999; Holdich et al., 1999; Souty Grosset et al., 2006). It is considered ‘‘vulnerable’’ in the IUCN Red List, ‘‘critically endangered’’ in the IUCN Italy, included in Annex II and V of EU Habitats Directive

2

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11

(92/43/EEC) and in Annex III and IV of Bern’s Convention. With more than 50 years of research, from Bott (1950) to Fratini et al. (2005), description of A. pallipes comprised of a maximum of four species and six subspecies (Starobogatov, 1995) to a minimum of one species (Trontelj et al., 2005) (Table 1). Accordingly, the nomenclature has been revised several times increasing the uncertainties (see Manganelli et al., 2006) and negatively impacting its conservation. Until now, the sequencing of mitochondrial genes has been the most common approach to investigate A. pallipes phylogeny, phylogeography and systematics (Grandjean et al., 2002; Fratini et al., 2005; Trontelj et al., 2005; Zaccara et al., 2005; Diéguez-Uribeondo et al., 2008; Stefani et al., 2011). However, studies based exclusively on mitochondrial data reflect the maternal pattern of inheritance, only providing a partial view of the complex phylogeographic patterns of white-clawed crayfish. Many recent studies show the contrast between mitochondrial and nuclear signal, in different zoological groups (Shaw, 2002; Koblmüller et al., 2007; Babbucci et al., 2010; Yang and Kenagy, 2009; Leaché, 2010; Rato et al., 2010): discrepancy between mitochondrial and

nuclear markers emerged when reconstructing the history of a species, extending the conflicts also in strongly support mtDNA clades. Concerns include polymorphisms in common ancestral populations, male-biased gene flow, selection on mtDNA, gene flow following hybridization and transfer of mtDNA to the nucleus (Moritz and Cicero, 2004). This make usefulness of single mitochondrial analysis controversial (Dasmahapatra and Mallet, 2006). Here we utilized for the first time an integrated approach, combining the effectiveness of COI sequences in phylogeographic analysis of genus Austropotamobius (Trontelj et al., 2005; Zaccara et al., 2005; Diéguez-Uribeondo et al., 2008; Stefani et al., 2011), with the well documented potential of the high throughput AFLP technique to assess genetic diversity within and among closely related taxa (Mendelson and Shaw, 2005; Maldini et al., 2006; Koblmüller et al., 2007). Present results based on the integration of mitochondrial data with nuclear markers challenge previously proposed hypothesis on the taxonomy of white-clawed crayfish, and suggest a novel evolutive and phylogenetic scenario for A. pallipes populations of Southern Europe.

Table 1 Systematic and taxonomic revisions of Austropotamobius pallipes complex (from Grandjean et al. (2000) modified). Authors

Genus

Species

Subspecies

Data

Lereboullet (1858) Bott (1950, 1972)

Austropotamobius Austropotamobius

pallipes pallipes

pallipes pallipes italicus lusitanicus

Morphology Morphology

Karaman (1963)

Austropotamobius

berndhauseri pallipes italicus

Albrecht (1982)

Austropotamobius

pallipes

Brodsky (1983)

Austropotamobius

pallipes italicus

Pretzmann (1987)

Austropotamobius

Starobogatov (1995)

Atlantoastacus

pallipes pallipes pallipes fulcisianus

orientalis lusitanicus pallipes italicus pallipes

Santucci et al. (1997)

Austropotamobius

Grandjean et al. (1998)

Austropotamobius

Grandjean et al. (2000)

Austropotamobius

Grandjean et al. (2002)

Austropotamobius

Fratini et al. (2005)

Austropotamobius

pallipes italicus

Trontelj et al. (2005)

Austropotamobius

pallipes

pallipes italicus

berndhauseri? pallipes italicus

Morphology italicus lusitanicus carsicus pallipes lombardicus carinthiacus trentinicus pallipes bispinosus italicus lusitanicus carsicus pallipes fulcisianus = italicus pallipes rhodanicus fulcisianus = berndhauseri + A. p. lombardicus + A. p. trentinicus italicus orientalis = carsicus carinthiacus

Morphology

Morphology

Morphology Morphology

Allozyme pallipes italicus lusitanicus pallipes italicus carsicus carinthiacus Not sampled pallipes italicus carsicus carinthiacus pallipes italicus meridionalis carsicus carinthiacus pallipes

mtDNA RFLP Morphology 16S rDNA

16S rDNA

16S rDNA

COI mtDNA

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11

2. Materials and methods 2.1. Collecting samples In total, 356 white-clawed crayfish were collected by hand and nets in 38 different sampling sites distributed from the Alps to the most southern peninsular areas of Italy, a ‘‘hotspot’’ for the diversity of the species (Souty Grosset et al., 2006). From each individual, a portion of muscle tissue was collected from fragments of pleiopods and chelae taken from live samples, avoiding the death of the animals. Crayfish were released alive in the same collection site after tissue sampling. Genomic DNA was extracted from all samples but to ensure high reliability of the data, only 199 samples from 35 populations having high quality DNA were retained for molecular analysis (Fig. 1, Table 2). In the next sections the sampled populations will be labeled with alpha numerical codes based on the 16S mtDNA haplotype (see Grandjean et al., 2002; Fratini et al., 2005) and sampling location following a northern–southern gradient. Populations previously characterized by 16S mtDNA (Grandjean et al., 2002; Fratini et al., 2005) were now submitted to COI haplotyping and AFLP fingerprinting. COI haplogroups acronyms used reflect those proposed by Trontelj et al. (2005). In addition, three individuals of Austropotamobius torrentium (Schrank, 1803) from the Magyar Termeszettudomanyi Muzeum of Budapest (Hungary) were included as out-group for both the mitochondrial and AFLP analysis. All tissue samples are actually

3

conserved in absolute ethanol at 20 °C, at the Laboratory of Molecular Zoology, Department of Evolutionary and Functional Biology, University of Parma, Parma (Italy). 2.2. DNA extraction and purification High molecular weight genomic DNAs were extracted and purified from ethanol-fixed muscle tissue stored at 20 °C. DNA was extracted either according to the classical SDS-proteinase K and phenol–chloroform technique described by Moore (1999) or alternatively by Aquapure genomic DNA kit (Biorad). DNA quality was visually inspected by 1% agarose gel electrophoresis in TAE buffer. All DNA samples are actually conserved in rehydration solution at 20 °C, at the Laboratory of Molecular Zoology, Department of Evolutionary and Functional Biology, University of Parma, Parma (Italy). 2.3. Sequencing of mitochondrial DNA COI gene fragment was amplified using the primers HCO 2198 and LCO 1490 (Folmer et al., 1994; Trontelj et al., 2005; Dawnay et al., 2007), which presumably should avoid the problem of pseudogenes detection, since they have been extensively used in many invertebrate and vertebrate species. A reaction volume of 25 ll containing 1U di GoTaq (Promega), Mg2+ 2.5 mM and dNTPs 0.2 mM was used. PCR profile was set as follows: 40 cycles of

Fig. 1. Sampling sites of the A. pallipes populations analysed. The different A. pallipes populations labeled with alpha numerical codes based on the 16S mtDNA haplotype (see Grandjean et al., 2002; Fratini et al., 2005), and sampling location following a northern–southern gradient.

4

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11 Table 2 Population identification number, sampling sites, 16S haplotype and number of analysed specimen of the A. pallipes and A. torrentium, submitted to COI characterization (bold) and AFLP fingerprinting. CS = 16S carsicus haplotype; CR = 16S carinthiacus haplotype, P = 16S pallipes haplotype, I = 16S italicus haplotype; M = 16S meridionalis haplotype (haplotypes according to Grandjean et al. (2002) and Fratini et al. (2005)). No.

Sampling sites

16S Haplotype

No. individuals AFLP

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36

Udine Rio Valturcana (Belluno) Roggia di Terlago (Trento) Rio Gamberi (Verbania) Risorgive del Ticino (Novara) Fontanile (Cremona) Bragamonti (Cuneo) Valle Peso (Cuneo) Arandolo (Cuneo) Torrente Rossana (Cuneo) Torrente Oxentina (Imperia) Rio Rialasco (Genova) Torrente Nenno (Genova) Torrente Chiusola (La Spezia) Torrente Stirone (Parma) Lago Pranda (Reggio Emilia) Torrente Zena (Bologna) Pistoia Rio Canvella (Prato) Il Mulino (Firenze) Castagno d’Andrea (Firenze) Rio Ceppeta (Arezzo) Fosso Doglio (Perugia) Fosso Argentina (Terni) Fosso Olpeta (Viterbo) Fosso Ariana (Rieti) Fosso dei Pantani (Rieti) Torrente Licenza (Roma) Fosso Duranna (Roma) Chieti Torrente Verde (Chieti) Fosso Vetrina (Chieti) Bussento (Salerno) Bussentino (Salerno) Villa D’Agri (Potenza) Magyar Termeszettudomanyi Muzeum of Budapest

CS1 CS2 CS3 CR4 CR5 CR6 P7 P8 P9 P10 P11 P12 P13 CR14 CR15 CR16 I17 I18 I19 I20 I21 I22 I23 M24 M25 M26 M27 M28 M29 M30 M31 M32 M33 M34 M35 A. torrentium

7 20 5 4 4 5 3 3 5 7 3 3 10 2 4 5 2 10 7 1 11 14 5 2 3 4 4 12 3 2 3 3 4 7 12 3

30 s at 95 °C, 45 s at 45 °C, and 1 min at 72 °C, after an initial 10 min denaturation step at 95 °C and a final extension at 72 °C for 10 min. PCR products were purified by elution from a 2.5% agarose gel, then precipitated with three volumes of 100% ethanol, and washed with 70% ethanol. Fragments sequencing was performed by CEQ™ DTCS-Quick Start Kit (Beckman Coulter) and analyzed on ‘‘CEQ™ 8000 DNA Analysis System’’ (Beckman Coulter). The amplified fragments consistently reached about 650 bp length. They were compared with sequences available in genomic databases using Blast and multiple alignments of sense and antisense sequences were conducted using such softwares as Clustal X (Thompson et al., 1997) and Sequencer 4.2 (Gene Code Corporation), verifying the correctness of the alignment also at the amino acid level. 2.4. Production and scoring of fluorescent AFLP markers AFLP markers were produced following the protocol already validated by Papa et al. (2005) and Maldini et al. (2006). Compared to the original protocol by Vos et al. (1995), herein the extract and purified genomic DNAs were cut by TaqI/EcoRI endonucleases, which showed higher polymorphism in animal genomes. Five primer combinations were used (in brackets the selective nucleotides): E32-T33 (eAAC-tAAG); E33-T32 (eAAG-tAAC); E32-T32 (eAAC-tAAC); E40-T37 (eAGC-tACG); E33-T37 (eAAG-tACG). An amount of 1.7 ll of PCR product and 0.3 ll of DNA internal size standards (CEQ DNA Size standard-600 Beckman-Coulter, Fullerton, CA) were added to 40 ll of deionized formamide (J.T. Baker, Phillipsburg, NJ). Samples were then loaded into the ‘‘CEQ™ 8000

No. individuals COI mtDNA 2 3 2 4

10 1

8

3

1 1

DNA Analysis System’’ (Beckman Coulter). Running conditions for capillary electrophoresis, raw data elaboration and reproducibility was tested according to Papa et al. (2005). A minimum fluorescence threshold value of 2000 was chosen on the basis of several comparative analyses, but lower peaks have sometimes been considered, especially for high molecular weight fragments when their resolution was comparable to those of ladder fragments having similar size. Fragments useful to assess genetic differentiation within species were identified as those reaching the threshold peak and being well distinct from nearby peaks (Papa et al., 2005). Analysed data was imported in Genographer software (Vers.1.6.0, Benham J.J., Montana State University 2001), for band scoring. Genographer allowed the construction of a virtual gel with bands shaped on the base of peak height, resolution and mobility and a thorough analysis of single fragments. Fragments in each individual were binary scored as 1 or 0, following the ‘‘band-based’’ approach (Bonin et al., 2007). To avoid genotyping errors in the AFLP procedure all the recommendations by Pompanon et al. (2005) and Bonin et al. (2007) were accomplished: (1) high quality DNA (2) standardized and validated procedures (3) replication of 30% of the samples, coupled with systematic duplication of samples (4) reliability of the results checked among different laboratories (5) double independent reading of the profiles carried out by researchers. 2.5. Phylogenetic and statistical analysis The Median Joining Network of COI gene sequences was built using Network 4.02-05 (http://www.fluxus-technology.com)

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11

according to Bandelt et al. (1999). It has been chosen to include missing haplotypes by generating median vectors. Moreover, it has been tested by Cassens et al. (2005) and Mardulyn et al. (2009) to be one of the most appropriate methods to infer genealogical relationships among intraspecific sequences, as one of the methods generating the smallest number of errors. The evolutionary tree of DNA sequences has been estimated by the maximum likelihood optimal criterion developed by Felsenstein and Churchill (1996). Over the traditional parsimony algorithms, which can give misleading results if rates of evolution differ in different lineages, the ML model assumes that each lineages and each site evolves independently and it allows the expected frequencies of the four bases to be unequal. The confidence level of the phylogenetic tree has been assessed by the bootstrapping technique proposed by Felsenstein (1985). The confidence for each clade of the observed tree has been estimated on the proportion of 1000 bootstrap trees showing that same clade. Finally, Factorial Correspondence Analysis of AFLP markers was undertaken by Genetix software (Belkhir et al., 1996–2002) where correspondences are adapted to the diploid genotypes and can be depicted graphically in 3D if desired. The FCA analyzed objects (group of individuals) are seen like a group of dots in a hyperspace which has as many dimensions as terms (alleles-bands) for all the variables (alleles-bands at different loci). The algorithm seeks the independent directions (orthogonal) in this hyperspace the length whose inertia – size which, by analogy with physics, represents the integral of the mass (here e.g. the number of individuals in a point of the hyperspace) multiplied by the square of the distance to the

5

center of the co-ordinates (here called center of gravity) – is maximized. These directions, which are defined by the appropriate vectors of the matrix (C–D from linear combinations of the original vector variables), determine a series of factorial axes. By convention, the first axis is that which has the strongest total contribution to the inertia. With AFLP bands, to use the individual genotypic data each individual is represented by his score for each term of each variable, that is to say 0 for the absence, 1 for the presence of the AFLP band. For each axis determined at the time of the analysis, a set of coefficients for each individual and of the locus are also calculated, as a matter of: 1. Of the absolute contributions: expressing the share taken by an element given (individual or alleles) in the inertia explained by a factor. 2. Of the relative contributions: expresses the share taken by the axis in the total contribution of the individual or the allele to the inertia (representing the dispersion of the cloud of the points). 3. Of the co-ordinates of all the points individuals and alleles with the various axes, safeguarded in a file. This last can be used in another graphic utility to draw groups of dots. These co-ordinates are also used to draw groups of dots, into two or three dimensions, to which it is possible to make undergo rotations and zooms, to visualize them under different angles. 3. Results 3.1. COI mitochondrial DNA Thirty-three samples belonging to nine populations representative of the geographic sampling area were selected for COI gene

Fig. 2. Median Joining Network of COI gene sequences. Median Joining Network obtained from the COI gene sequences with mutational steps numbered. For example, P13_1 means P haplotype group, population 13, specimen 1. See Table 1 for details. WEURO: Western Europe haplotype group (Trontelj et al., 2005); NWITAL: North Western Italy haplotype group (Trontelj et al., 2005); APPEN: Apennine haplotype group (Trontelj et al., 2005); SEA-WB: South Eastern Alps–Western Balkans haplotype group (Trontelj et al., 2005); ISTRIA: Istria haplotype group (Trontelj et al., 2005); CASER: Caserta sequences from Southern Apennines (Aceto et al., 2008).

6

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11

Fig. 3. Maximum Likelihood tree inferred by the COI gene sequences analysis. For example, P13 1 means P haplotype, population 13, specimen 1. See Table 1 for details.

analysis (Table 2). We amplified and sequenced 650 bp of the gene, while the alignments were based on 399 bp GenBank reference sequences (Supplementary Table, Supporting Information for online publication only). A range of 5–25 mutational steps were detected among the aligned sequences of different haplotype groups of white-clawed crayfish, with limited or null transversions

(Table 3). The total number of nucleotide substitutions between the A. pallipes and the outgroup A. torrentium ranged between 45 and 52 (considering both transitions and transversions) (Table 3). A detailed description of nucleotide diversity is reported in Table 3. Median Joining Network revealed a strong geographic structure with at least four distinct haplogroups separated by a variable

Table 3 Total number of nucleotide substitutions between A. pallipes and the outgroup A. torrentium, and among A. pallipes lineages (Tr = Transitions; Tv = Transversions). WEurope = Western Europe haplotype group; NWItaly = North Western Italy haplotype group; SEAlps-WBalkans = South Eastern Alps–Western Balkans haplotype group; Apennine = Apennine haplotype group (geographical abbreviations according to Trontelj et al. (2005)). Sequences of A. pallipes were grouped to detect haplogroups on the basis of highly diagnostic mutations.

A. torrentium WEurope NWItaly SEAlps-WBalkans Apennine

A. torrentium

WEurope

NWItaly

SEAlps-WBalkans

Apennine

– 52 48 45 48

52 – 23 25 23

48 23 – 11 12

45 (30 Tr-15 Tv) 25 (21 Tr-4 Tv) 11 (11 Tr-0 Tv) – 5 (5 Tr-0 Tv)

48 (30 Tr18 Tv) 23 (18 Tr-5 Tv) 12 (12 Tr-0 Tv) 5 (5 Tr-0 Tv) –

(33 (30 (30 (30

Tr-19 Tr-18 Tr-15 Tr-18

Tv) Tv) Tv) Tv)

(33 Tr-19 Tv) (18 Tr-5 Tv) (21 Tr4 Tv) (18 Tr-5 Tv)

(30 Tr-18 Tv) (18-5 Tv) (11 Tr-0 Tv) (12 Tr-0 Tv)

7

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11

number of mutations (Fig. 2). Accordingly, Maximum Likelihood tree revealed four groups supported by bootstrap values ranging between 650 and 1000 (Fig. 3) and consistent with the haplogroups previously described by Trontelj et al. (2005). No match with the additional Istria reference sequence was identified. The Istrian haplotype identified by Trontelj et al. (2005) was not found and therefore omitted in Table 3. In detail, the populations from Northern and Central Italy (CR4, CR5, CR15 and I18) grouped with the ‘‘North Western Italy’’ haplotype group, far apart from sequences of the population of Western Italy P13, which were grouped in the haplotype group of ‘‘Western Europe’’ (seven polymorphic sites). The sequences from populations M28 and M35 of Central and Southern Italy, clustered together with the ‘‘Apennine’’ haplotype group; while the specimen from North Eastern Italy, belonging to the populations CS1 and CS2, were assigned to the ‘‘South Eastern Alps–Western Balkans’’ haplotype group. Genetic variability within haplogroups ranged between 7 and 32 point mutations considering both our samples and reference sequences while within the populations analysed in this work the number was comprised between 5 and 10.

3.2. AFLP markers A total of 1010 AFLP reliable fingerprinting profiles were obtained genotyping 202 samples (199 A. pallipes and 3 A. torrentium) with 5 EcoRI/TaqI primer combinations carrying three selective nucleotides (Fig. 4). The repeatability was greater than 99%. A total of 535 AFLP bands in a range 50–500 bp, 479 of which polymorphic (89%) and only 56 monomorphic across the dataset, were binary scored. The mean number of bands per primer combination was 107; the monomorphic bands ranged from 8 to 15, while the polymorphic ones ranged from 86 to 107 depending on the primer pairs (Table 4).

Fig. 4. Genographer virtual gel with AFLP profiles. Fingerprinting profiles of single crayfish (columns) and AFLP loci (rows) detected by Genographer software.

The expected heterozygosis calculated within populations ranged from 0.07 (S.d. = 0.02) in M24 to 0.31 (S.d. = 0.01) in M35, with a mean value of 0.17 (S.d. = 0.01). The Gst indexes were 0.48 among the whole dataset (outgroup included) and 0.46 among A. pallipes populations indicating an overall genetic differentiation between different populations. The Nei’s distances between A. pallipes populations and A. torrentium ranged from 0.27 (CR15 vs A. torrentium) to 0.59 (P13 vs torrentium) while within the pallipes populations the distances ranged between 0.06 (CR5 vs I17) and 0.32 (P12 vs M33). Focusing on samples belonging to the same COI haplotype group, Nei’s genetic distances based on AFLP polymorphisms were in the range 0.15–0.38 for N Western Italy (respectively I18 vs CR4 and I18 vs CR15), 0.16 for Apennine (M28 vs M35), 0.14 for S Eastern Alps–Western Balkans (CS1 vs CS2) . The Factorial Correspondence Analysis (FCA) at individual level (Fig. 5) showed a clear separation between A. torrentium and A. pallipes only. Indeed, no clear clusters gathering populations of A. pallipes were detected among the 35 analyzed groups. The inertia explained by the first three axes were 15.53%; 10.51% and 8.68% respectively. In spite of the mitochondrial divergence, the nuclear genome showed a weak differentiation. For example, individuals belonging to I18 and P13 described as different species by Grandjean et al. (2002) on the basis of mitochondrial 16S lineages, showed COI mtDNA separated by 23 nucleotide mutational steps, but they shared 144 monomorphic AFLP bands out of 535 (27%). Only 42 bands out of 535 (9%) were diagnostic (fixed in a group and absent in the other).

4. Discussion Until now, sequencing of mitochondrial gene fragments has been the most common method to study crayfish phylogeny, either for the genus Austropotamobius, Euastacus (Clark, 1936) (Shull et al., 2005), Orconectes (Cope, 1872) (Mathews et al., 2008), Cambaroides (Faxon, 1884) (Ahn et al., 2006; Braband et al., 2006), and Cherax (Erichson, 1846) (Munasinghe et al., 2004; Nguyen et al., 2004). Analyses of COI gene fragment of A. pallipes from widely distributed Italian populations identified four major haplogroups, overlapping those previously described by Trontelj et al. (2005) and Stefani et al. (2009). However, results were only partially matching the mitochondrial lineages described by Grandjean et al. (2002) and Fratini et al. (2005) with a different marker. More precisely, CR (carinthiacus) and I (italicus) lineages, considered distinct subspecies by the 16S mtDNA (Grandjean et al., 2002; Fratini et al., 2005), in this work resulted in the same COI haplogroup defined ‘‘North Western Italy’’ (assimilated to previously defined italicus lineage). As a consequence, the hypothesis of two different subspecies, already considered doubtful (Machino, 1997; Trontelj et al., 2005), has also been heavily challenged by our results.

Table 4 Number of monomorphic, polymorphic, total loci, and percentage of polymorphism scored by AFLP, totally and for each primer combination. Primer combination

Monomorphic loci

Polymorphic loci

Total

% Polymorphism

eAGC tACG eAAG tACG eAAG tAAC eAAC tAAC eAAC tAAG Total 5

12 9 8 15 12 56

86 104 86 96 107 479

98 113 94 111 119 535

88 92 91 86 90 90

8

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11

Fig. 5. Two-dimensional FCA of the 202 analysed individuals generated by Genetix software. Individuals were grouped to detect putative species A. pallipes–A. italicus (sensu Grandjean et al., 2002; Fratini et al., 2005) genotyped by AFLP. Squares: A. pallipes; Triangles: A. italicus; Circles: A. torrentium.

The median-joining network showed both a differentiation between lineages and a clear geographic pattern reflecting local genetic variation. Each population is characterized by a single haplotype, even if consistent mutational steps could be identified among sequences of the same population. However, incongruence related to the effectiveness in species recognition by mitochondrial sequences have been demonstrated in the recent years (Song et al., 2008; Galtier et al., 2009), mainly due to the identification of nuclear mitochondrial pseudogenes in crustacean genomes (Schneider-Broussard and Neigel, 1997; Bucklin et al., 1999; Bensasson et al., 2001; Williams and Knowlton, 2001; Nguyen et al., 2002; Song et al., 2008; Galtier et al., 2009). Similarly, Galtier et al. (2009) showed that in such particular groups as arthropods, characterization of mtDNA for taxonomy and systematics is not completely reliable due to the common hybrid introgression. Therefore, an integrated approach is potentially useful to investigate the composite relationships among species/subspecies and populations, and AFLPs revealed to be a valid class of nuclear markers to be coupled to mitochondrial analysis. As explained by Vernesi and Bruford (2009), conservation genetics has largely remained a field were a small number of molecular markers are applied to few populations of a single species, some of which may be threatened. The usefulness of AFLP markers for taxonomy and definition of closely related species has been largely proved in many zoological groups, as arthropods (Mendelson and Shaw, 2005; Ming and Wang, 2006). The AFLP markers possess features suitable for animal conservation genetics (Vernesi and Bruford, 2009): for surveys of intrapopulation variability and genetic structures among populations, AFLPs have proved to be effective as standards markers (Vernesi and Bruford, 2009). Furthermore, AFLPs can be successfully used for phylogenetics inferences (Sullivan et al., 2004). Flexibility, no need for

extensive technique development, and low cost make AFLPs potentially effective for conservation genetics studies (Lucchini, 2003). The limited use of AFLPs in animal conservation genetics seems even more paradoxical (Vernesi and Bruford, 2009). The percentage of polymorphic AFLPs (90%) in crayfish agreed with those reported in several studies on crustacean and insects, ranging from 65% to 95% (Fezner and Crandall, 1999; Mendelson and Shaw, 2002; Salvato et al., 2002; Cannas et al., 2003; Gili et al., 2004; Mendelson et al., 2004; Sun et al., 1999; Shuster et al., 2005), up to 99% in the Ortoptera Locusta migratoria manilensis (Meyen, 1835) (Zheng et al., 2006). Only few AFLP markers among five hundred fragments were highly diagnostic, and produced an unexpected limited overall differentiation between individuals and lineages. Only 9% of polymorphic bands were highly diagnostic between PAL (pallipes) and ITA (italicus) lineages, while in sibling species like Helicoverpa armigera (Hübner, 1805) and Helicoverpa assulta (Guenèe, 1852) (Lepidoptera) the percentage of species-specific bands was at least threefold higher (about 30% Ming and Wang, 2006). Despite the differentiated mitochondrial haplotype groups, AFLP fingerprinting indicated that a consistent part of the genome is shared among distant populations of A. pallipes. In particular, individuals belonging to the two putative species A. pallipes (P) and A. italicus (I) were clearly overlapping. AFLPs on A. pallipes therefore open a new perspective on its evolutionary history, suggesting the existence of a unique species in Italy. In similar cases AFLPs solved disputes regarding the taxonomy of putative species and subspecies of crustaceans, such as the ones belonging to genus Eriocheir (De Haan, 1835) (Lu et al., 2000) and Artemia (Leach, 1819) (Tryantaphyllidis et al., 1997), and revealed that the common pearl oysters Pinctada fucata (Gould, 1850) is monospecific (Yu et al., 2006).

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11

Incongruence between mitochondrial and nuclear signals were already reported in the phylogeny of Arthropods (Gomez-Zurita and Vogler, 2003), and in other taxa when comparing mtDNA sequencing with polymorphic markers as AFLPs (Egger et al., 2007; Koblmüller et al., 2007). The controversy is not only limited to our work, but it reflects the contradiction of DNA evolution in such groups, as the freshwater decapods, where mitochondrial and nuclear DNA could show different evolutionary histories. So how can we explain this complex scenario? The distribution of mitochondrial haplotype groups allowed Trontelj et al. (2005) to pinpoint the possible glacial refugia, such as microrefugia in Southern France (Alpes-Maritimes), NorthWestern Italy (Liguria), Apennine and Istrian region, the latter considered the primary center of radiation of A. pallipes. After the split between A. torrentium and A. pallipes in the Miocene (8–12 MYA); the separation of the basal mitochondrial lineages of A. pallipes (Western Europe, North Western Italy, South Eastern Alps– Western Balkans) occurred during the Pliocene (3–6 MYA), while the split of the Apennine lineage occurred during the Pleistocenic glaciations (0.1–1.5 MYA) (Trontelj et al., 2005). Despite different methodological approaches, our paper is in agreement with Santucci et al. (1997) that reached the following conclusion: ‘‘a wide genetic heterogeneity that appears mainly related to range fragmentation and subsequent recolonizations from multiple refugia during the glacial events’’. Such a pattern indicates that, beyond a general expansionspeciation model (Hewitt, 1996), the white-clawed crayfish experienced a more complex evolutionary history. This consideration was driven by the existence of several microrefugia for other freshwater invertebrates (Pongratz et al., 2003; Verovnik et al., 2004) and fish such as Thymallus thymallus (Linnaeaus,1758) (Weiss et al., 2002), Cottus gobio (Linnaeaus, 1758) (Englbrecht et al., 2000), and Leuciscus (Cuvier, 1816) (sin. Squalius Berg, 1949) species complex (Costedoat et al., 2006). The Balkan microrefugia generated the ‘‘South Eastern Alps–Western Balkans’’ (carsicus) lineage while the Southern France (Costedoat et al., 2006), and Alps Maritimes (Trontelj et al., 2005; Stefani et al., 2011) microrefugia differentiated the ‘‘Western Europe’’ (pallipes) and the ‘‘North Western Italy’’ (italicus + carinthiacus) lineages respectively. Also Central and Southern Apennines have been considered important refugial areas for Italian fauna (Trontelj et al., 2005), from which the ‘‘Apennine’’ (meridionalis) lineage emerged. Interestingly, Canestrelli et al. (2006) demonstrated that similar evolutionary events in amphibians were driven by fragmentation of the Apennine chain as separated islands during plio-pleistocene periods. To summarize, the white-clawed crayfish likely originated in the headwaters of the Eastern Alps–Western Balkans, migrated westward onto the Apuane Alps, and then proceeded southward through the Apennine in the Central and Southern Italy. During the Pleistocene, the alteration of glacial and interglacial phases caused the sea level variations that resulted in a minimum Mediterranean level 100–130 m below the present day level (Pielou, 1979). This process produced the emergence of the upper and middle Adriatic and Tyrrhenian continental shelves, favoring the river confluence of downstream courses (Bianco, 1990). In particular, the ancient location of the Po river delta in the Central-Southern Adriatic area (Southern Italy) permitted the direct connection between Italian and Croatian rivers. This event might have enabled several crayfish settling between the West and East Adriatic coasts, supporting the Eastern Adriatic regions/ Balkans as the center of original radiation (Trontelj et al., 2005). Bottleneck and genetic drift have then played a major role in shaping the present-day distribution of the mitochondrial haplotype groups. In fact, the number of median vectors connecting

9

the network, suggests the existence of several intermediate unsampled lineages, likely already extinct. The distribution of the haplotypes might be also affected by the dispersal capability and gene flow of A. pallipes – a k-selected species – nowadays limited by habitat fragmentation, and severe demographic fluctuations, while in ancient times were enhanced by consistent river connection and better environmental quality. Finally, artificial propagations such as the well known translocation to the Iberian Peninsula from Northwestern Italy (Grandjean et al., 2002; Trontelj et al., 2005) and the introduction of crayfish from the Apennine group to Southern France and Northwestern Italy (Zaccara et al., 2004) probably had a massive influence on the distribution of A. pallipes in recent years. In conclusion, the analysis of COI gene sequences revealed a zoogeographic pattern not detectable in the nuclear genome of A. pallipes. In our view, AFLP fingerprinting witnesses a secondary intensive contact among the different phylogeographic lineages occurred during the postglacial events, due to the changes of habitats and freshwater drainages. This scenario was proposed also in other groups where nuclear genomic variation was not correlated withmtDNA differentiation (Koblmüller et al., 2007; Yang and Kenagy, 2009). In particular, Egger et al. (2007) found significant incongruence between phylogenetic reconstructions from mitochondrial DNA and AFLP data, suggesting incomplete sorting of mitochondrial haplotypes as well as frequent introgression between differentiated lineages of cichlid fishes. Hence, in absence of complete reproductive barriers, hybridization of nuclear genomes has been shown to be possible, and sister taxa still hybridize despite millions of years of divergence demonstrated by molecular data (Hewitt, 2000). Therefore, considering experimental molecular data and the above cited ecological considerations, the previously proposed haplotypes of A. pallipes are clearly nested within a unique taxon and do not support the existence of two sister species A. pallipes and A. italicus. Thus, we hereby named all the white-clawed crayfish populations as A. pallipes being the recognized species name due to priority. Acknowledgments This research was supported by the Italian Ministry of University and Research, as a Ph.D. Project at the University of ‘‘Roma Tre’’, Rome (Italy). The authors would like to thank people and Authorities for their help with collecting samples of A. pallipes: G. Caricato (Basilicata Region), M. De Biaggi (Piedmont), G. De Luise (Friuli Venezia-Giulia), F. Gherardi and S. Fratini (Tuscany), T. Pagliani (Abruzzo and Marche), E. Marconato (Veneto and Lombardia), S. Salvidio (Liguria), M. Zanetti (Veneto) and Lazlo Forrò from the Magyar Termeszettudomanyi Muzeum of Budapest (Hungary) for A. torrentium samples. Authors would like to thank anonymous referees for improving suggestions and also thank Dr. Kristin Brabender for English language corrections. Appendix A. Supplementary material Supplementary data associated with this article can be found, in the online version, at doi:10.1016/j.ympev.2011.03.031. References Aceto, S., De Vico, G., Sica, M., Maio, N., Carella, F., Cataldi, M., Canestrini, F., Froglia, C., 2008. Mitochondrial COI gene sequences of Austropotamobius pallipes complex from Reggia di Caserta (Italy). Direct Submission. Available from: www.ncbi.nlm.nih.gov.

10

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11

Ahn, D.A., Kawai, T., Kim, S.J., Rho, H.S., Jung, J.W., Kim, W.L., Byung, J., Kim, M.S., Min, G.S., 2006. Phylogeny of northern hemisphere freshwater crayfishes based on 16S rRNA gene analysis. Korean J. Genet. 28, 185–192. Albrecht, H., 1982. Das system der europäischen flusskrebse (Decapoda, Astacidae): vorschlag und begrundung. Mitt. Hamburg. Zool. Museum Inst. 79, 187–210. Austin, C.M., 1996. Systematics of the freshwater crayfish genus Cherax Erichson (Decapoda: Parastacidae) in northern and eastern Australia: electrophoretic and morphological variation. Aust. J. Zool. 44, 259–296. Babbucci, M., Buccoli, S., Cau, A., Cannas, R., Goñi, R., Díaz, D., Marcato, S., Zane, L., Patarnello, T., 2010. Population structure, demographic history, and selective processes: contrasting evidences from mitochondrial and nuclear markers in the European spiny lobster Palinurus elephas (Fabricius, 1787). Mol. Phylogenet. Evol. 56, 1040–1050. Bandelt, H.J., Forster, P., Rohl, A., 1999. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48. Belkhir, K., Borsa, P., Chikhi, L., Rafauste, N., Catch, F., 1996–2002. GENETIX4.04, software under Windows TM for the genetics of the populations. Laboratory Genome, Populations, Interactions, CNRS UMR 5000, University of Montpellier II, Montpellier, France. Bensasson, D., Zang, D.X., Hartl, D.L., Hewitt, G.M., 2001. Mitochondrial pseudo genes: evolution’s misplaced witnesses. Trends Ecol. Evol. 16, 314–321. Bianco, P.G., 1990. Potential role of the paleohistory of the Mediterranean and Paratethys basin on the early dispersal of Europe-Mediterranean freshwater fishes. Ichthyol. Explor. Freshwaters 1, 167–184. Bonin, A., Ehrich, D., Manel, S., 2007. Statistical analysis of amplified fragment length polymorphism data: a toolbox for molecular ecologists and evolutionists. Mol. Ecol. 16, 3737–3758. Bott, R., 1950. Die flüssekrebse europas (Decapoda, Astacidae). Proc. Senckenberg Nat. Soc. 483, 1–36. Bott, R., 1972. Besiedlungsgeschichte und systematik der Astacieden West-Europas unter besonderer beruchksichtigung der shweizt. Rev. Suisse Zool. 79, 387– 408. Braband, A., Kawai, T., Scholtz, G., 2006. The phylogenetic position of the east Asian freshwater crayfish Cambaroides within the Northern Hemisphere Astacoidea (Crustacea, Decapoda, Astacida) based on molecular data. J. Zool. Syst. Evol. Res. 44, 17–24. Brodsky, S.Y., 1983. On the systematics of paleartic crayfishes (Crustacea, Astacidae). Freshwater Crayfish 5, 464–470. Bucklin, A., Guarnieri, M., Hill, R.S., Bentley, A.M., Kaartvedt, S., 1999. Taxonomic and systematic assessment of planktonic copepods using mitochondrial COI sequence variation and competitive, species-specific PCR. Hydrobiologia 401, 239–254. Buhay, J.E., Crandall, K.A., 2005. Subterranean phylogeography of freshwater crayfishes shows extensive gene flow and surprisingly large population sizes. Mol. Ecol. 14, 4259–4273. Canestrelli, D., Cimmaruta, R., Costantini, V., Nascetti, G., 2006. Genetic diversity of yellow-bellied toad Bombina pachypus, with implications for conservation. Mol. Ecol. 15, 3741–3754. Cannas, A., Deiana, A.M., Tagliavini, J., 2003. AFLP-fingerprinting delle aragoste mediterranee Palinurus elephas e P. mauritanicus. Dati preliminari. Biol. Mar. Medit. 10, 379–382. Cassens, I., Mardulyn, P., Milinkovitch, M., 2005. Evaluating intraspecific ‘‘network’’ construction methods using simulated sequence data: do existing algorithms outperform the global maximum parsimony approach? Syst. Biol. 54, 363–372. Costedoat, C., Chappaz, R., Barascud, B., Guillard, O., Gilles, A., 2006. Heterogeneous colonization pattern of European Cyprinids, as highlighted by the dace complex (Teleostei: Cyprinidae). Mol. Phylogenet. Evol. 41, 127–148. Crandall, K.A., Harris, D.J., Fetzner Jr., J.W., 2000. The monophyletic origin of freshwater crayfish estimated from nuclear and mitochondrial DNA sequences. Proc. Roy. Soc. Lond. B 267, 1679–1686. Dasmahapatra, K.K., Mallet, J., 2006. DNA barcodes: recent success and future prospects. Heredity 97, 254–255. Dawnay, N., Ogden, R., McEwing, R., Carvalho, G.R., Thorpe, R.S., 2007. Validation of the barcoding gene COI for use in forensic genetic species identification. Forensic Sci. Int. 173, 1–6. Diéguez-Uribeondo, J., Royo, F., Souty-Grosset, C., Ropiquet, A., Grandjean, F., 2008. Low genetic variability of the white clawed crayfish in the Iberian Peninsula: its origin and management implications. Aquat. Conserv. 18, 19–31. Egger, B., Koblmüller, S., Sturmbauer, C., Sefc, K.M., 2007. Nuclear and mitochondrial data reveal different evolutionary processes in the Lake Tanganyika cichlid genus Tropheus. BMC Evol. Biol. 7, 137. Englbrecht, C.C., Freyhof, J., Nolte, A., Rassman, K., Schliewen, U., Tautz, D., 2000. Phylogeography of the bullhead Cottus gobio Pisces: Teleostei: Cottidae suggests a pre-Pleistocene origin of the major central European populations. Mol. Ecol. 9, 709–722. Felsenstein, J., 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39, 783–791. Felsenstein, J., Churchill, G.A., 1996. A hidden Markov model approach to variation among sites in rate of evolution. Mol. Biol. Evol. 13, 93–104. Fezner Jr., J.W., Crandall, K.A., 1999. Genetic variability within and among populations of the golden crayfish (Orconectes luteus): a comparison using amplified fragments length polymorphisms (AFLPs) and mitochondrial 16S gene sequences. Freshwater Crayfish 12, 396–412. Fezner Jr., J.W., Crandall, K.A., 2003. Linear habitats and the nested clade analysis: an empirical evaluation of geographic versus river distances using an Ozark crayfish (Decapoda: Cambaridae). Evolution 57, 2101–2118.

Folmer, O., Black, M., Hoeh, W., Lutz, R., Vrijenhoek, R., 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 5, 294–299. Fratini, S., Zaccara, S., Barbaresi, S., Grandjean, F., Souty-Grosset, C., Crosa, G., Gherardi, F., 2005. Phylogeography of the threatened crayfish (genus Austropotamobius) in Italy: implications for its taxonomy and conservation. Heredity 94, 108–118. Galtier, N., Nabholz, B., Glemin, S., Hurst, G.D.D., 2009. Mitochondrial DNA as a marker of molecular diversity: a reappraisal. Mol. Ecol. 18, 4541–4550. Gherardi, F., Holdich, D.M., 1999. Crayfish in Europe as Alien Species. How to Make the Best of a Bad Situation? A.A. Balkema, Rotterdam. Gili, M., Monaghan, M.T., Spaak, P., 2004. Amplified fragments length polymorphisms (AFLP) reveals species-specific markers in the Daphnia galeata-hyalina species complex. Hydrobiologia 526, 63–71. Gomez-Zurita, J., Vogler, A.P., 2003. Incongruent nuclear and mitochondrial phylogeographic patters in the Timarcha gottingensis species complex Coleptera, Crisomelidae. J. Evol. Biol. 16, 833–843. Grandjean, F., Gouin, N., Frelon, M., Souty-Grosset, C., 1998. Genetic and morphological systematic studies on the crayfish Austropotamobius pallipes (Decapoda: Astacidae). J. Crustacean Biol. 18, 549–555. Grandjean, F., Harris, F.D.J., Souty-Grosset, C., Crandall, K.A., 2000. Systematic of the European endangered crayfish species, Austropotamobius pallipes (Decapoda: Astacidae). J. Crustacean Biol. 20, 522–529. Grandjean, F., Frelon-Raimond, M., Souty-Grosset, C., 2002. Compilation of molecular data for the phylogeny of the genus Austropotamobius: one species or several? Bull. Fr. Pêche Pisc. 367, 671–680. Hewitt, G.M., 1996. Some genetic consequences of ice ages and their role in divergence and speciation. Biol. J. Linn. Soc. 58, 247–276. Hewitt, G.M., 2000. The genetic legacy of the Quaternary ice aces. Nature 405, 907– 913. Holdich, D.M., 2002. Biology of Freshwater Crayfish. Blackwell Science, Oxford. Holdich, D.M., Rogers, W.D., Reynolds, J.D., 1999. Crayfish in the British Isles. In: Gherardi, F., Holdich, D.M. (Eds.), Crayfish in Europe as Alien Species. How to Make the Best of a Bad Situation? A.A. Balkema, Rotterdam, pp. 221–235. Karaman, M.S., 1963. Studie der Astacidae (Crustacea, Decapoda). Hydrobiologia 22, 111–132. Koblmüller, S., Duftner, N., Sefc, K.M., Aibara, M., Stipacek, M., Blanc, M., Egger, B., Sturmbauer, C., 2007. Reticulate phylogeny of gastropod-shell-breeding cichlids from Lake Tanganyika – the result of repeated introgressive hybridization. BMC Evol. Biol. 7, 7. Leaché, A.D., 2010. Species trees for spiny lizards (genus Sceloporus): identifying points of concordance and conflict between nuclear and mitochondrial data. Mol. Phylogenet. Evol. 54, 162–171. Lereboullet, A., 1858. Description de deux nouvelles espèces d’écrevisses de nos rivières. Mémories Soc. Sci. Nat. 5, 1–11. Lu, R., Qui, T., Xiang, C., Chaomei, X.H., Zhang, J., 2000. RAPD and AFLP techniques for the analysis of genetic relationships in two genera of Decapoda. Crustaceana 73, 1027–1036. Lucchini, V., 2003. AFLP: a useful approach for biodiversity, conservation and management. CR Soc. Biol. 326, S43–S48. Machino, Y., 1997. New white-clawed crayfish Austropotamobius pallipes (Lereboullet, 1858) occurrences in Carinthia, Austria. Bull. Fr. Peche Pisc. 347, 713–720. Machino, Y., Holdich, D.M., 2006. Distribution of crayfish in Europe and adjacent countries: updates and comments. Freshwater Crayfish 15, 292–323. Maldini, M., Nonnis Marzano, F., González Fortes, G., Papa, R., Gandolfi, G., 2006. Fish and seafood traceability based on AFLP markers: elaboration of a species database. Aquaculture 261, 487–494. Manganelli, G., Favilli, L., Fiorentino, V., 2006. Taxonomy and nomenclature of Italian white-clawed crayfish. Crustaceana 79, 633–640. Mardulyn, P., Cassens, I., Milinkovitch, M., 2009. A comparison of methods for constructing evolutionary networks from intraspecific DNA sequences. In: Bertorelle, G., Bruford, M.W., Hauffe, H.C., et al. (Eds.), Population Genetics for Animal Conservation. Cambridge University Press, Cambridge, pp. 104–120. Mathews, L.M., Adams, L., Anderson, E., Basile, M., Gottardi, E., Buckholt, M.A., 2008. Genetic and morphological evidence for substantial hidden biodiversity in a freshwater crayfish species complex. Mol. Phylogenet. Evol. 48, 126–135. Mendelson, T.C., Shaw, K.L., 2002. Genetic and behavioral components of the cryptic species boundary between Laupala cerasina and L. kohalensis (Orthoptera: Grillidae). Genetica 116, 301–310. Mendelson, T.C., Shaw, K.L., 2005. Use of AFLP markers in surveys of arthropod diversity. Methods Enzymol. 395, 161–177. Mendelson, T.C., Siegel, A.M., Shaw, K.L., 2004. Testing geographic pathways of speciation in a recent island radiation. Mol. Ecol. 13, 3787–3796. Ming, Q.L., Wang, C.Z., 2006. Genetic differentiation of Helicoverpa armigera (Hübner) and H. Assulta (Guenée) (Lepidoptera: Noctuidae) based on AFLP markers. Insect Sci. 13, 437–444. Moore, D., 1999. Manipulation of DNA. In: Ausubel, F.M., Brent, R., Kingston, R.E., et al. (Eds.), Current Protocols in Molecular Biology. John Wiley & Sons Inc., New York, pp. 2.1.1–2.2.3. Moritz, C., Cicero, C., 2004. DNA barcoding: promises and pitfall. Plos Biol. 2, e354. Munasinghe, D.H.N., Burridge, C.P., Austin, C.M., 2004. Molecular phylogeny and zoogeography of the freshwater crayfish genus Cherax Erichson (Decapoda: Parastacidae) in Australia. Biol. J. Linn. Soc. 81, 553–563. Nguyen, T.T.T., Murphy, N.P., Austin, C.M., 2002. Amplification of multiple copies of mitochondrial Cytochrome b gene fragments in the Australian freshwater

S. Chiesa et al. / Molecular Phylogenetics and Evolution 61 (2011) 1–11 crayfish, Cherax destructor Clark (Parastacidae: Decapoda). Anim. Genet. 33, 304–308. Nguyen, T.T.T., Austin, C.M., Meewan, M.M., Schultz, M.B., Jerry, D.R., 2004. Phylogeography of the freshwater crayfish Cherax destructor Clark (Parastacidae) in inland Australia: historical fragmentation and recent range expansion. Biol. J. Linn. Soc. 83, 539–550. Papa, R., Troggio, M., Ajmone-Marsan, P., Nonnis Marzano, F., 2005. An improved protocol for the production of AFLP markers in complex genomes by means of capillary electrophoresis. J. Anim. Breed. Genet. 122, 62–68. Pielou, E.C., 1979. Biogeography. Wiley, New York. Pompanon, F., Bonin, A., Bellemain, E., Taberlet, P., 2005. Genotyping errors: causes, consequences and solutions. Nat. Rev. Genet. 6, 1–13. Pongratz, N., Storhas, M., Carranza, S., Miciels, N.K., 2003. Phylogeography of competing sexual and parthenogenetic forms of a freshwater flatworm: patterns and explanations. BMC Evol. Biol. 3, 23. Pretzmann, G., 1987. A contribution to a historic analysis of Mediterranean freshwater decapods chorology. Invest. Pesq. 51, 17–25. Rato, C., Carranza, S., Perera, A., Carretero, M.A., Harris, D.J., 2010. Conflicting patterns of nucleotide diversity between mtDNA and nDNA in the Moorish gecko, Tarentola mauritanica. Mol. Phylogenet. Evol. 56, 962–971. Salvato, P., Battisti, A., Concato, S., Masutti, L., Patarnello, T., Zane, L., 2002. Genetic differentiation in the winter pine processionary moth (Thaumetopoea pityocampa-wilkinsoni complex) inferred by AFLP and mitochondrial DNA markers. Mol. Ecol. 11, 2435–2444. Santucci, F., Iaconelli, M., Andreani, R., 1997. Allozyme diversity of European freshwater crayfish of the genus Austropotamobius. Bull. Fr. Peche Pisc. 347, 663–676. Schneider-Broussard, R., Neigel, J.E., 1997. A large-subunit mitochondrial ribosomal DNA sequence translocated to the nuclear genome of two stone crabs (Menippe). Mol. Biol. Evol. 14, 156–165. Scholtz, G., 2002. Phylogeny and evolution. In: Holdich, D.M. (Ed.), Biology of Freshwater Crayfish. Blackwell Science, Oxford. Shaw, K.L., 2002. Conflict between nuclear and mitochondrial DNA phylogenies of a recent species radiation: what mtDNA reveals and conceals about modes of speciation in Hawaiian crickets. Proc. Natl. Acad. Sci. USA 99, 16122– 16127. Shull, H.C., Pérez-Losada, M., Blair, D., Sewell, K., Sinclair, A.E., Lawler, S., Ponniah, M., Crandall, K.A., 2005. Phylogeny and biogeography of the freshwater crayfish Euastacus (Decapoda: Parastacidae) based on nuclear and mitochondrial DNA. Mol. Phylogenet. Evol. 37, 249–263. Shuster, M.S., Miller, M.P., Lang, B.K., Zorich, N., Huynh, L., Keim, P., 2005. The effects of controlled propagation on an endangered species: genetic differentiation and divergence in body size among native and captive populations of the Socorro Isopod (Crustacea: Flabellifera). Conserv. Genet. 6, 355–368. Song, H., Buhay, J.E., Whiting, M.F., Crandall, K.A., 2008. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudo genes are co amplified. Proc. Natl. Acad. Sci. USA 105, 13486–13491. Souty Grosset, C., Holdich, D.M., Noël, P., Reynolds, J.D., Haffner, P., 2006. Atlas of Crayfish in Europe. Muséum de National d’Histoire Naturelle, Paris. Starobogatov, Y.I., 1995. Taxonomy and geographical distribution of crayfishes of Asia and East Europe (Crustacea: Decapoda: Astacidae). Arthropoda Selecta 4, 3–25.

11

Stefani, F., Zaccara, S., Delmastro, G.B., Buscarino, M., 2011. The endangered whiteclawed crayfish Austropotamobius pallipes (Decapoda, Astacidae) east and west of the Maritime Alps: a result of human translocation? Conserv. Genet. 12, 51– 60. Sullivan, J.P., Lavoue, S., Arnegard, M.E., Hopkins, C.D., 2004. AFLPs resolve phylogeny and reveal mitochondrial introgression within a species flock of African electric fish (Mormyroidea: Teleostei). Evolution 58, 825–841. Sun, Y., Song, W.Q., Zhong, Y.C., Zhang, R.S., Abatzopoulos, T.J., Chen, R.Y., 1999. Diversity and genetic differentiation in Artemia species and populations detected by AFLP markers. Int. J. Salt Lake Res. 8, 341–350. Thompson, J.D., Gibson, T.J., Plewniak, F., Jeanmougin, F., Higgins, D.G., 1997. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl. Acids Res. 25, 4876–4882. Trontelj, P., Machino, Y., Sket, B., 2005. Phylogenetic and phylogeographic relationships in the crayfish genus Austropotamobius inferred from mitochondrial COI gene sequences. Mol. Phylogenet. Evol. 34, 212–226. Tryantaphyllidis, G.V., Criel, G.R.L., Abatzopoulos, T.J., Thomas, K.M., Peleman, J., Beardmore, J.A., Sorgeloos, P., 1997. International study on Artemia. LVII. Morphological and molecular characters suggest conspecificy of all bisexual European and North African Artemia populations. Mar. Biol. 129, 477–487. Vernesi, C., Bruford, M.W., 2009. Recent developments in molecular tools for conservation. In: Bertorelle, G., Bruford, M.W., Hauffe, H.C., et al. (Eds.), Population Genetics for Animal Conservation. Cambridge University Press, Cambridge, pp. 321–344. Verovnik, R., Sket, B., Trontelj, P., 2004. Phylogeography of subterranean and surface populations of water lice Asellus aquaticus (Crustacea: Isopoda). Mol. Ecol. 13, 1519–1532. Vos, P., Hogers, R., Bleeker, M., 1995. AFLP: a new technique for DNA fingerprinting. Nucl. Acids Res. 23, 4407–4414. Weiss, S., Persat, H., Eppe, R., Scholotter, C., Uiblein, F., 2002. Complex patterns of colonization and refugia revealed for European grayling Thymallus thymallus, based on complete sequencing of the mitochondrial DNA control region. Mol. Ecol. 11, 1393–1407. Williams, S.T., Knowlton, N., 2001. Mitochondrial pseudo genes are pervasive and often insidious in the snapping shrimp genus Alpheus. Mol. Biol. Evol. 18, 1484– 1493. Yang, D., Kenagy, G.J., 2009. Nuclear and mitochondrial DNA reveal contrasting evolutionary processes in populations of deer mice (Peromyscus maniculatus). Mol. Ecol. 18, 5115–5125. Yu, D.H., Jia, X., Chu, K.H., 2006. Common pearl oysters in China, Japan, and Australia are conspecific: evidence from ITS sequences and AFLP. Fisheries Sci. 72, 1183– 1190. Zaccara, S., Stefani, F., Galli, P., Nardi, P.A., Crosa, G., 2004. Taxonomic implications in conservation management of the white-clawed crayfish (Austropotamobius pallipes) (Decapoda, Astacidae) in Northern Italy. Biol. Conserv. 10, 1–10. Zaccara, S., Stefani, F., Crosa, G., 2005. Diversity of mitochondrial DNA of the endangered white-clawed crayfish (Austropotamobius italicus) in the Po River catchment. Freshwater Biol. 50, 1262–1272. Zheng, X., Zong, Y., Duan, Y., LI, C., Dang, L., Guo, Y., Ma, E., 2006. Genetic variation and population structure of oriental migratory locust, Locusta migratoria manilensis, in China by Allozyme, SSRP-PCR, and AFLP Markers. Biochem. Genet. 44, 333–347.