- Research article
- Open Access
Recent evolution of alternative reproductive modes in the 'living fossil' Triops cancriformis
BMC Evolutionary Biology volume 7, Article number: 161 (2007)
The Notostraca is a small but ancient crustacean order with a contrasting combination of a conservative morphology and a wide range of reproductive modes. The tadpole shrimp Triops cancriformis, includes bisexual – the putatively ancestral state -, androdioecious and hermaphrodite populations. As hermaphroditism and androdioecy confer a colonisation advantage, we expect the postglacial colonisation of northern Europe to have been effected by lineages with such reproductive modes. Therefore, N European populations should be composed of closely related lineages reflecting a recent range expansion. In contrast, glacial refugia in the south should contain bisexual populations with high haplotype diversity and more population structuring. To test these hypotheses, we analysed the geographic distribution of reproductive modes based on new and published sex ratio data. In addition, we investigated the European phylogeography of T. cancriformis by sequencing over a 1000 bp of mitochondrial DNA (mtDNA) in individuals from a large sample of populations of the three recognised subspecies.
Bisexual populations were only found in the Iberian Peninsula, with the rest of European populations showing low male proportions or no males. Androdioecious populations were found in Central and Eastern Europe. Regarding mtDNA diversity, Spanish and Moroccan populations of T. c. mauritanicus were highly divergent, and showed strong population structure. In contrast, Triops c. cancriformis and T. c. simplex formed a single mtDNA lineage with low haplotype diversity. This diversity was structured into two phylogenetic clades (A, B), coexisting in E Germany. Basal haplotypes of both lineages were found in the Iberian Peninsula. Most of the populations in clade A and B are either hermaphroditic or androdioecious, with the only bisexual population in these clades found in the Iberian Peninsula. The genetic divergence between these two clades suggests a split in the Late Pleistocene and their geographic distribution reflects a complex evolutionary history of European Triops populations, with possibly two episodes of range expansions – one of them by clade A – involving androdioecious and hermaphroditic populations.
As we predicted, N European populations of T. cancriformis are closely related, with few widely distributed haplotypes and indications of a recent range expansion involving hermaphroditic/androdioecious lineages. A possible second range expansion or long distance colonisation may have created the secondary contact zone between T. c. cancriformis/simplex clades A and B. The large haplotype diversity and strong genetic subdivision in the Iberian Peninsula, which is known to contain only bisexual populations, strongly suggest that this area was a Pleistocene refugium for T. cancriformis, although the occurrence of additional eastern refugia cannot be ruled out. Our data support the status of T. c. mauritanicus as a separate species and the colonisation of N Africa from the Iberian Peninsula. We suggest that hermaphroditism/androdioecy has evolved recently in T. cancriformis and has facilitated the postglacial colonisation of northern Europe.
The widespread occurrence of outcrossing sexual reproduction has puzzled evolutionary biologists since Darwin, as reproduction either through self-fertilisation or parthenogenesis provides immediate advantages [1, 2]. One of such key advantages is evident during the colonisation process [3, 4], as a single selfing hermaphrodite or parthenogenetic individual can found a new population. Organisms with mixed reproductive strategies allow a unique direct comparison of the relative advantages of different reproductive modes [5, 6]. Mixed strategies often have a geographic component (e.g. 'geographical parthenogenesis') which reflects the interplay between historical and selective factors . The process of recolonisation of northern areas during interglacial periods of the Pleistocene might have favoured passively dispersing parthenogens or selfing hermaphrodites. As a consequence, in temperate organisms with mixed reproductive strategies, hermaphroditism is often found in areas previously covered by glaciers or permafrost and therefore, recolonisation must have taken place after the last glacial maximum. Despite this, ecological correlates have been intensively sought to explain such geographic patterns adaptively [5, 7], whereas the phylogeographic context of mixed reproductive patterns has only recently begun to be investigated [8–10].
The Eurasian tadpole shrimp Triops cancriformis (Crustacea: Branchiopoda: Notostraca) has a mixed reproductive strategy which includes hermaphroditism, androdioecy (consisting of hermaphrodites and a small proportion of males) and bisexuality [11, 12], the latter being considered to be the ancestral condition . The currently recognized subspecies T. c. simplex (N Africa and NE Iberian Peninsula) and T. c. mauritanicus (NW Morocco and SW Iberian Peninsula) are considered to be bisexual, with equal or male biased sex ratios [14, 15]. The nominal subspecies T. c. cancriformis, occurring in central and northern Europe and a disjunct area in Japan, includes hermaphrodite and androdioecious populations [16, 17]. The reproductive strategy in T. c. cancriformis populations is controversial, although there is evidence to suggest that 'females' are in fact hermaphrodites as they can reproduce in isolation  and ovotestes are present . However, females from a German T. c. cancriformis population apparently lacked testicular tissue  suggesting that parthenogenesis – or a mixture of bisexual and hermaphroditic reproduction in some populations – cannot be ruled out. Our preliminary microsatellite data  show significant heterozygote deficiencies in many populations, and genotyping of individuals reared in isolation and their offspring supports hermaphroditism. Thus, and following Sassaman  we consider some T. c. cancriformis populations to be androdioecious, as populations consist of variable proportions of hermaphrodites and males, with some populations being made solely of selfing hermaphrodites. As in other androdioecious branchiopods, T. c. cancriformis hermaphrodites can reproduce either by selfing or by outcrossing with males.
T. cancriformis inhabits temporary freshwater ponds and all three reproductive modes lead to the production of resistant diapausing cysts, which survive in the pond sediments during drought periods. These cysts are also the means of passive dispersal by wind, water currents or downstream floods, and birds [21–23], amphibians or hoofed animals [24–26]. An indication of the vagility of such cysts is that populations occur on remote islands, and are apparently found wherever suitable habitats are available . Evidence for wide notostracan distribution before the Pleistocene can be derived from its abundant fossil record dating back to the Carboniferous or possibly up to the Devonian period . In fact, the striking morphological similarity of some Upper Triassic Triops sp. fossils from Germany and extant T. cancriformis  makes this notostracan one of the best examples of morphological stasis and 'living fossils' . Thus, the potentially high dispersal abilities of their diapausing cysts and the possibility of hermaphrodite reproduction could account for the wide distribution of T. cancriformis.
Here we present a phylogeographic analysis of European T. cancriformis in the context of its reproductive mode variation. So far, attempts at investigating the genetic variability of T. cancriformis have been hampered by its low genetic variability [31, 32], therefore our phylogeographic survey of T. cancriformis is based on sequence variation on over 1000 bp of mitochondrial DNA containing highly informative mitochondrial genes. We hypothesize, in agreement with Longhurst , that hermaphrodite or androdioecious populations should occur in areas unsuited for the species during Pleistocene glacial maxima (either covered with ice sheets or permafrost), and that therefore must have been recently colonised. Further, N European populations should be composed of one or few closely related lineages reflecting a recent range expansion of hermaphroditism and/or androdioecy. In contrast, glacial refugia in S Europe should contain bisexual populations with high haplotype diversity and high population structuring. To test these hypotheses, we analysed the geographic distribution of inferred reproductive modes based on reviewed and re-analysed data concerning sex-ratio in this species and new critical data of our own. In addition, we screened a large sample of populations of the three recognised subspecies for nucleotide sequence variation in two mtDNA genes (COI and ATPase). The phylogeographic pattern is interpreted in the context of the inferred population reproductive mode. Our results provide insights into the evolution of reproductive mode and population diversification in Triops.
mtDNA haplotype diversity
A total of 84 T. cancriformis individuals from 29 ponds were sequenced for two mitochondrial gene fragments: a 568 bp sequence of the cytochrome c oxidase 1 gene (COI) and a 506 bp fragment comprising partial sequences of ATP synthase FO subunits 8 (ATP8) and 6 (ATP6) genes (Table 1, Figure 1). Given that the partition-homogeneity test did not find significant differences between both fragments (p-value of 0.07), we combined them into a single alignment with a total length of 1074 bp for phylogenetic analysis.
Overall 213 variable nucleotide sites, 158 of them parsimony informative, were identified in the combined alignment. Excluding the Iberian (Extremadura, EXT; Ares del Maestre, ARE) and Moroccan (Youssofia, YOI) T. c. mauritanicus specimens, 38 sites remained variable and 26 sites were parsimony informative. Comparing the COI and ATPase (including ATP8 and 6 genes) fragments separately we found within the COI fragment a total of 20 substitutions on 1st codon position, 3 on 2nd and 91 on 3rd positions, which resulted in nine amino acid changes. The ATPase fragment is characterised by 46 changes on 1st codon position, 22 on 2nd and 99 at 3rd codon position which result in 33 amino acid changes. The combined alignment was moderately A+T rich (mean AT content = 65.2%).
Despite sequencing over 1000 bp of mtDNA, overall haplotype diversity was low, with only 19 mtDNA haplotypes identified in the combined alignment. Eleven of those haplotypes were found in the 14 T. c. cancriformis populations sampled. Fourteen out of the nineteen identified haplotypes were found in single populations (Table 1, Figure 2), while one haplotype (H8) was found in two nearby populations in Germany (DAX, HAG), H13 was found in a set of four nearby ponds in the area of Königswartha, and H1, H2 and H3 were common and widespread, found in six, nine and three ponds respectively. H1 was found in the Iberian individuals from Espolla (ESP) and El Puig (PUI), which had been previously identified as T. c. simplex , as well as in individuals of the nominal subspecies from Döberitzer Heide (DOE), Königswartha pond 12 and 21 (KOE12, KOE21) and Rhine (DAX).
The number of haplotypes per population ranged from 1–4 (average 1.35), with the maximum found in Königswartha (KOE21) (Table 1). Given the low level of genetic diversity found in our preliminary analysis, we optimised our resources by increasing the number of bp sequenced per individual and the number of populations surveyed, instead of the number of individuals sequenced per population.
Both maximum likelihood and Bayesian analyses produced identical topologies showing two highly divergent and strongly supported monophyletic lineages (Figure 3). The first lineage is represented by T. c. mauritanicu s haplotypes (nucleotide diversity π = 0.05905) (clade C). The second lineage includes all haplotypes from the subspecies T. c. cancriformis and T. c. simplex (π = 0.00920). Despite the low genetic diversity within T. c. cancriformis/simplex, two divergent clades (A and B) were identified. The closely related haplotypes H1 to H9 within clade A occurred from E Spain to N and C Europe and the Japanese population (π = 0.00114) (see network in Figure 2). This group includes the three most common and widespread haplotypes (H1, H2, H3) and a cluster of rare, geographically restricted haplotypes differing from them in single substitutions. In contrast, haplotypes from the southern populations within clade A (H10, Sicily and H11, UBA in Spain) are more divergent and differ in 9 and 6 respectively substitutions from H3. The two morphologically described subspecies T. c. cancriformis and T. c. simplex, which show different reproductive modes, are not represented by distinct monophyletic lineages, as T. c. simplex individuals from some East Iberian populations (PUI, ESP) share haplotype H1 with five T. c. cancriformis individuals from four different populations in Germany. Clade B was only found in three disjunct locations (KOE, DAN, UBA, see Figure 3). Haplotypes from clades A and B were present in Königswartha populations. Population UBA from eastern Spain contained the basal haplotype of clade A (H11) – the most internal haplotype in the network – and a haplotype of unresolved position (H12) according to the phylogenetic analyses, but which appears as most closely related to clade B from our network analysis under the limits of parsimony (Figure 2), and therefore was included in clade B. The root of the T. c. cancriformis/simplex lineage lies in between both haplotypes in this population, which seems to have retained ancestral polymorphisms in the species.
Sequence divergence and approximation of divergence times
Based on the combined alignment, pair-wise nucleotide differences among all T. cancriformis haplotypes ranged from 1 to 137 bp which translate into a maximum distance (HKY + G corrected) of 14.38% (H5 vs. H17). The mean corrected distances between clades for the combined alignment and for the COI and ATPase fragments separately are shown in Table 2. The genetic distance between clades A and B, including the two recognized subspecies T. c. cancriformis and T. c. simplex, is less than 3.5%. The mean genetic distance was considerable higher (>10%) between haplotypes from clade A and B together (T. c. cancriformis and T. c. simplex) and clade C (T. c. mauritanicus).
Calculating the corrected distances separately for the COI and ATPase fragment resulted in higher values for the ATPase alignment. This illustrates, that ATPase is evolving faster and provides more information regarding intraspecific relationships than COI (Table 2).
The crustacean COI calibration applied to obtain an approximate time of divergence between the T. c. cancriformis/simplex group (clades A and B) and T. c. mauritanicus (clade C), yielded 4.71–7.85 mya (around the Miocene-Pliocene boundary) (Table 2). The divergence time between clades A and B was estimated to have happened during the Early Pleistocene (0.68–1.13 mya).
Distribution of inferred reproductive modes in Triops cancriformis
Bisexuality was supported in two Portuguese populations of T. c. mauritanicus (Algarve) and one Spanish population of T. c. simplex (Espolla) (Table 3, Figure 4). Androdioecy was inferred for populations in Central and E Europe (Table 3). Populations with male proportions above 0.100 were identified mainly from eastern European populations (Poland and Hungary). As listed in Table 3, male-less populations occur in Germany, France, Poland and United Kingdom, and no males have been reported from Italian, Polish and Japanese populations (with no reported sample sizes) [17, 34–36]. As androdioecious populations can have low proportion of males it can not be excluded that those population are actually androdioecious. In populations where males have been recorded male proportions ranged from 0.012 to 0.281 (Table 3). As shown in Figure 4 both androdioecious and putatively hermaphrodite populations are distributed where no suitable habitats for T. cancriformis were found during the Pleistocene ice ages.
As we predicted, T. cancriformis haplotypes from most N European populations and the disjunct Japanese population – clustered in clade A – are closely related, with a few, widely distributed haplotypes. The reduced genetic diversity in T. cancriformis is a consistent and remarkable result of our data. This can be supported by the large mitochondrial fragment we investigated, the selection of the most rapid evolving gene in Triops and extensive sampling. A low mitochondrial diversity in T. cancriformis was noted by Mantovani et al.  and Korn et al.  based on 16S and 12S sequences, which are slowly evolving compared to our fragments. Furthermore, preliminary data on diversity of nuclear microsatellite loci [20, 37] suggested low genome wide polymorphism in T. c. cancriformis/simplex . This pattern supports postglacial range expansion of T. c. cancriformis/simplex lineages into areas which were unsuitable for this species during the last glacial maximum. When we take the phylogeographic and reproductive mode results together, the predominance of androdioecious and potentially hermaphroditic populations in formerly glacial or permafrost regions strongly suggests that lineages with alternative reproductive modes were responsible for the postglacial recolonization of N Europe, as hypothesized by Longhurst . Our results with T. cancriformis closely resemble the pattern found in geographic parthenogenesis. The term geographical parthenogenesis is used for the tendency of parthenogenetic populations to be distributed in high latitudes reflecting an association between parthenogenesis and environments that were strongly affected by the Pleistocene glacial cycles [6, 38]. Several organisms have been described having parthenogenetic populations in N and C Europe and sexual populations in restricted southern areas [39, 40]. Although T. cancriformis populations, being hermaphroditic, do not fulfil the geographical parthenogenesis concept, they would enjoy the colonisation advantage in a comparable way. Therefore we suggest broadening the concept of geographical parthenogenesis by introducing the term 'geographical hermaphroditism'.
In contrast to the pattern found in N Europe with populations of low genetic diversity and alternative reproductive modes, southern T. cancriformis populations, including all the identified bisexual populations, show high and geographically structured haplotype diversity. These results support the hypothesis that bisexual populations are characterised by reduced gene flow and high interpopulation divergence (regional differentiation as found between Iberian T. c. mauritanicus populations), possibly because the areas inhabited by these populations have been suitable for longer. Similar strong geographic structuring of genetic variation with very isolated populations has been reported from other sexual European large branchiopods (Anostracans) [41, 42]. These strong population structures despite high potential for gene flow through their diapausing propagules suggests that Triops populations and other large branchiopods, in a similar way to continental zooplankton undergo a process of 'Monopolisation' of their habitats . The presence of ancestral haplotypes and strong genetic subdivision in the Iberian Peninsula, strongly points to this area as a Pleistocene refugium for T. cancriformis. This inference is further supported by the presence of bisexual populations for T. c. mauritanicus and T. c. simplex in the Iberian Peninsula. However, the genetic divergence between clades A and B and their geographic distribution indicates that the phylogeographic history of European T. cancriformis is more complex as these clades are likely to have diverged earlier than the last glacial cycle. Clade B is especially puzzling as it is present only in three locations including northern and southern populations. Different scenarios could explain this pattern. One of these is the possibility that N European populations from clade B are derived from an early postglacial expansion which left relict populations in an Eastern refugium – Balkans, Middle East – which we did not sample, a pattern that has previously been reported for the European fire salamander . Another explanation could be long distance colonisation events facilitated by the passive dispersal of Triops diapausing cysts as reported for other aquatic invertebrates [45, 46]. Either way, this possible range expansion/long distance colonisation event – comparable to clade A – was linked to androdioecious lineages, which is the reproductive mode, present in clade B. In summary we suggest a split in the Late Pleistocene that originated clades A and B in T. c. cancriformis/simplex, with an early refugium in the Iberian Peninsula and possibly a second one. Further sampling from eastern populations and further Mediterranean islands where the species is known to occur should shed light into this issue. Subsequently these lineages underwent possibly two episodes of range expansions/long distance colonisation, in both cases involving hermaphroditic and/or androdioecious populations. Interestingly, those populations we infer to be androdioecious, occur in the area where both clades coexist in Germany. The occurrence of hybridisation between these clades is suggested by our microsatellite data  as two of the populations where both lineages coexist (KOE) and have been analysed with microsatellite loci display Hardy Weinberg equilibrium .
We suggest that geographic isolation within the Iberian Peninsula, possibly dating back to the Miocene-Pliocene boundary, led to the initial divergence between T. c. mauritanicus and T. c. cancriformis/simplex lineages and also to genetic isolation among T. c. mauritanicus populations in the Northwest and Southeast Iberian Peninsula. Such regional persistence suggests that long-term suitable habitats were available in separate areas within the Iberian Peninsula throughout the Pliocene and Pleistocene. Deep phylogeographic splits within the Iberian Peninsula have been observed in a wide range of taxa, indicating either the existence of several independent Pleistocene refugia within this area  or colonisation of the southern Iberian Peninsula from North Africa [48, 49]. Our analyses, containing Moroccan samples, indicate that the Iberian populations of T. c. mauritanicus are basal, which supports a scenario in which T. c. mauritanicus would have colonised North Africa from the Iberian Peninsula, possibly during the early Pliocene. A similar pattern of colonisation from the Iberian Peninsula into North Africa has been documented for lizards and amphibians [50–52]. Therefore, we suggest that the Iberian Peninsula provided a long-term subdivided glacial refugium for both T. c. mauritanicus and T. c. cancriformis/simplex.
The combined alignment of COI and ATPase sequences showed that T. c. mauritanicus and T. c. cancriformis/simplex form two highly divergent monophyletic lineages, as also suggested by a 16S and 12S analysis . The estimated 13.55% sequence divergence between both lineages is within the range separating other notostracan  and crustacean species . Our data therefore provide strong support to the proposal of Korn et al.  to reinstate T. c. mauritanicus to species status, as originally described by Ghigi . In contrast, the traditional subdivision into the subspecies T. c. cancriformis and T. c. simplex sensu Longhurst  could not be supported by our data as the populations representing T. c. simplex (ESP, PUI, UBA) shared either a genetic haplotype with T. c. cancriformis populations or were very closely related to them. We propose that the validity of T. cancriformis simplex should be reconsidered.
Our analysis of T. cancriformis reproductive mode in the context of its phylogeography allowed us to gain insights into the role of alternative reproductive modes during Pleistocene range shifts. We suggest that hermaphroditism/androdioecy has evolved recently, possibly in the Late Pleistocene and has facilitated the postglacial colonisation of northern Europe from one or two refugia, one of them likely placed in the Iberian Peninsula. Despite generally low genetic variability in T. cancriformis we were able to identify mitochondrial regions which provided sufficient resolution to reveal phylogeographic patterns across Europe. Our data support a taxonomic revision of T. cancriformis. The occurrence of three different reproductive modes within closely related genetic lineages of a single species make T. cancriformis a good model to investigate the evolution of mixed reproductive strategies.
A total of 29 water bodies in inundated floodplains, isolated ponds, puddles and fish nursery pools containing T. cancriformis populations were sampled in Europe, Japan and Morocco between 2000 and 2006 (Table 1, Figure 1). Overall 23 populations belong to the nominal subspecies, T. c. cancriformis, 3 populations were identified as the subspecies T. c. simplex and 3 populations belong to T. c. mauritanicus (Table 1). Samples included specimens collected using a dip net (5 mm pore diameter), and sediment containing diapausing cysts from dry ponds where the species was known to occur. Once specimens were needed for genetic analyses, and to be able to obtain sex ratio data, sediment was placed at the bottom of glass tanks and distilled water added to induce hatching of cysts. Individuals were grown until maturity and fed daily with fish flakes ad libitum.
Amplification and Sequencing
Total genomic DNA was isolated from ethanol-preserved tissue using commercial DNA extraction kits (Invisorb Spin Forensic Kit, Invitek; PureGene Kit, Gentra Systems). For the populations DAN and DAX, where attempts to hatch specimens from the sediment were unsuccessful, DNA was extracted from single cysts isolated from the sediment following procedures described in Gómez & Carvalho  (Table 1).
Given that our preliminary studies based on sequence variation at mitochondrial 16S and COI genes provided little resolution for intraspecific relationships within T. cancriformis  we decided to explore the published T. cancriformis and T. longicaudatus mitochondrial genomes in order to identify the most variable regions for a phylogeographic survey. We designed primers within highly conserved regions identified from an alignment of four branchiopoda sequences, Triops cancriformis (GenBank accession number NC_004465), Triops longicaudatus (NC_006079), Daphnia pulex, (DQ340834) and Artemia franciscana (X69067). The set of primers are available from the authors under request. We then carried out a preliminary survey of genetic variation for nine mitochondrial genes (ATP6, ATP8, cytb, ND2, ND3, ND5, ND4, ND4L and the D-loop) using eight T. cancriformis individuals from different populations. As a result ATP8 and ATP6 in combination with COI proved to be the most diverse fragments in terms haplotype and nucleotide diversity (data available from the authors upon request). Consequently these loci were chosen to complement our preliminary COI data.
A 568 bp fragment of COI was amplified and sequenced using primers LCO1490 and HCO2198 . A second fragment of 506 bp, covering parts of the ATP8 and parts of the ATP6 gene, was amplified using newly designed Triops specific primers. The 5' primer, ATP-Tc-3800U (5'-TCTAATCCCCCAAATAGCCCCTAGT-3'), is located in the tRNAAsp (corresponding to sites 3800 to 3825 in the reference sequence NC_004465), and the 3' primer, ATP-Tc-4405L (5'-ACAGCAAGAGTTCCGGGACGAATA-3') is located in the ATP6 mtDNA (sites 4381 to 4405 in the reference sequence NC_004465).
Amplifications were performed in 20 μl final volume containing 2 μl template DNA, 1.5 mM MgCl2, 200 μM of each nucleotide, 100 μM of each primer, 0.01 U of Taq DNA polymerase and 1× NH4-PCR Buffer (Bioline). The following cycling conditions were used: 3 min denaturating at 93°C, (45 s at 94°C, 45 s at 50°C, 1 min at 72°C) × 35, 5 min extension at 72°C. PCR products were sequenced directly using the PCR primers and the Beckman DTCS Quick Start Sequencing kit in a Beckman CEQ8000 capillary sequencer. The sequences were checked by eye with the CEQ8000 data analyser and initially aligned with Clustal W  and finally adjusted by hand. Sequences for the COI, ATP6 and ATP8 fragments were obtained for all individuals. All the sequences were deposited in GenBank (accession numbers EF675826–EF675991).
Phylogenetic relationships were reconstructed using maximum-likelihood and Bayesian approaches. The best fitting, least-parameter rich model of sequence evolution was based on hierarchical likelihood-ratio test in the program ModelGenerator0.6 . This method identified the HKY+G model of sequence evolution and four gamma distributed rate categories as the optimal model. This model of sequence evolution and its parameter estimates were used to perform a maximum likelihood (ML) algorithm in PHYML 2.4.4 . Branch support of the ML tree was assessed by 1000 bootstrap pseudo-replications. MRBAYES 3.1.1  was used to generate a Bayesian reconstruction by running a Markov chain for 2,000,000 generations and using a burn-in of 2500. The 50% consensus tree was obtained from the remaining trees from both runs sampled after the initial burn-in period.
Nodes in the ML-tree were considered to be well supported if they showed at least 70% bootstrap support . Similarly, strong branch support was inferred in the Bayesian analyses when posterior probabilities were over 80% . Trees were displayed with NJPlot . To examine the geographic and reproductive mode associations of haplotypes a haplotype network was constructed using TCS1.21 .
Molecular clock estimates
Unfortunately, the abundant fossil record of Triopsidae is not informative for molecular calibration purposes mainly due to the morphological stasis of the group since the Carboniferous period (about 302 Ma) [28, 68–70]. Therefore, currently, fossil data provide no information on splits within Triops. In consequence, no direct calibration points were available within the Triops phylogeny precluding the application of methods to estimate divergence times such as suggested by Sanderson . In order to provide a rough approximation of divergence times among T. cancriformis lineages we used two crustacean clock calibrations of Knowlton and Weigt  (1.4% sequence divergence per million years) and Schubart et al.  (1.66–2.33% sequence divergence per million years). Both calibrations are based on COI data and consequently we used only the COI fragment from our data set to estimate divergence times. Corrected average pairwise genetic sequence (using the best fit model for the COI dataset, TVM+G) between haplotypes were used in all divergence time estimates.
Sex ratio and reproductive mode
Sex ratio was determined for all T. cancriformis populations sampled in this study either from wild-caught specimens or from individuals hatched in the laboratory from sediment samples. Sexing was carried out on fixed specimens under a dissecting microscope. We considered an individual was a male when no egg sac was found in the 11th pair of thoracopods. No external diagnostic character allows discriminating between females and hermaphrodites; therefore, we calculated the sex ratio as the number of males divided by the number of female/hermaphrodites. Furthermore we included published data for 21 European populations in order to understand the geographic distribution of sex ratios. We only included literature data that provided either sample sizes and the sex ratio, or the raw numbers of males and females/hermaphrodites per sample. Sex ratio was used to distinguish between three potential reproductive modes for each population: bisexual reproduction was inferred when the sex ratio did not significantly differ from 1 based on a binominal Chi-square test; androdioecy was inferred when males where present but the sex ratio deviated significantly from 1; populations in which no males where found were treated as potentially hermaphroditic populations. Note, however, that such populations could also represent androdioecy with low frequency of males.
Otto SP, Lenormand T: Resolving the paradox of sex and recombination. Nat Rev Genet. 2002, 3: 252-261. 10.1038/nrg761.
Bell GD: The Masterpiece of Nature: The Evolution and Genetics of Sexuality. 1982, Berkley and Los Angeles: University of California Press
Baker HG: Self-compatibility and establishment after "long-distance" dispersal. Evolution. 1955, 9: 347-348. 10.2307/2405656.
Baker HG: Support for Baker's Law-as a rule. Evolution. 1967, 21: 853-856. 10.2307/2406780.
Peck JR, Yearsley JM, Waxman D: Explaining the geographic distributions of sexual and asexual populations. Nature. 1998, 391: 889-892. 10.1038/36099.
Kearney MR: Hybridization, glaciation and geographical parthenogenesis. Trends Ecol Evol. 2005, 20: 495-502. 10.1016/j.tree.2005.06.005.
Kearney MR: Why is sex so unpopular in the Australian desert?. Trends Ecol Evol. 2003, 18: 605-607. 10.1016/j.tree.2003.09.021.
Schon I, Gandolfi A, Di Masso E, Rossi V, Griffiths HI, Martens K, Butlin RK: Persistence of asexuality through mixed reproduction in Eucypris virens (Crustacea, Ostracoda). Heredity. 2000, 84: 161-169. 10.1046/j.1365-2540.2000.00647.x.
Kearney MR, Blacket MJ, Strasburg JL, Moritz C: Waves of parthenogenesis in the desert: evidence for the parallel loss of sex in a grasshopper and a gecko from Australia. Mol Ecol. 2006, 15: 1743-1748. 10.1111/j.1365-294X.2006.02898.x.
Law JH, Crespi BJ: The evolution of geographic parthenogenesis in Timema walking-sticks. Mol Ecol. 2002, 11: 1471-1489. 10.1046/j.1365-294X.2002.01547.x.
Sassaman C: Sex ratio variation in female-biased populations of Notostracans. Hydrobiologia. 1991, 212: 169-179. 10.1007/BF00025998.
Zaffagnini F, Trentini M: The distribution and reproduction of Triops cancriformis (Bosc) in Europe (Crustacea Notostraca). Monitore Zool Ital (NS). 1980, 14: 1-8.
Pannel JR: The evolution and maintenance of androdioecy. Ann Rev Ecol Syst. 2002, 33: 397-425. 10.1146/annurev.ecolsys.33.010802.150419.
Machado M, Cristo M, da Fonseca LC: Non-cladoceran branchiopod crustaceans from southwest Portugal. I Occurrence notes. Crustaceana. 1999, 72: 591-602. 10.1163/156854099503636.
Boix D, Sala J, Moreno-Amich R: Population dynamics of Triops cancriformis (Crustacea: Branchiopoda: Notostraca) of the Espolla temporary pond in the northeastern Iberian peninsula. Hydrobiologia. 2002, 486: 175-183. 10.1023/A:1021350802144.
Engelmann M, Hahn T: Vorkommen von Lepidurus apus, Triops cancriformis, Eubranchipus (Siphonophanes) grubei, Tanymastix stagnalis und Branchipus schaefferi in Deutschland und Österreich (Crustacea: Notostraca und Anostraca). Faun Abh. 2005, 25: 3-67.
Scanabissi Sabelli F, Eder E, Cesari M: Male occurrence in Austrian populations of Triops cancriformis and ultrastructural observations of the male gonad. Invertebr Biol. 2005, 124: 57-65.
Longhurst AR: Evolution in the Notostraca. Evolution. 1955, 9: 84-86. 10.2307/2405360.
Engelmann M, Hahn T, Hoheisel G: Ultrastructural characterization of the gonads of Triops cancriformis (Crustacea, Notostraca) from populations containing both females and males: no evidence for hermaphroditic reproduction. Zoomorphologie. 1997, 117: 175-180. 10.1007/s004350050042.
Zierold T: Morphological variation and genetic diversity of Triops cancriformis (Crustacea: Notostraca) and their potential for understanding the influence of postglacial distribution and habitat fragmentation. 2006, Freiberg: Technische Universität Bergakademie Freiberg, PhD thesis
Càceres CE, Soluk DA: Blowing in the wind: a field test of overland dispersal and colonisation by aquatic invertebrates. Oecologia. 2002, 131: 402-408. 10.1007/s00442-002-0897-5.
Green AJ, Figuerola J: Recent advances in the study of long-distance dispersal of aquatic invertebrates via birds. Diversity and Distributions. 2005, 11: 149-156. 10.1111/j.1366-9516.2005.00147.x.
Brendonck L, Riddoch BJ: Wind-borne short-range egg dispersal in anostracans (Crustacea: Branchiopoda). Biol J Linn Soc. 1999, 67: 87-95. 10.1006/bijl.1998.0293.
Frank PW: Conchostraca. Palaeogeogr Palaeoclimatol Palaeoecol. 1986, 62: 399-403. 10.1016/0031-0182(88)90064-8.
Thiery A: Multispecies coexistence of branchiopods (Anostraca, Notostraca & Spinicaudata) in temporary ponds in Chaouia plain (western Morocco): sympatry or syntopy between usually allopatric species. Hydrobiologia. 1991, 212: 117-136. 10.1007/BF00025994.
Bohonak AJ, Whiteman HH: Dispersal of the fairy shrimp Branchinecta coloradensis (Anostraca): Effects of hydroperiod and salamanders. Limnol Oceanogr. 1999, 44: 487-493.
Longhurst AR: A Review of the Notostraca. Bull Br Mus (Nat Hist) Zool. 1955, 3: 1-57.
Tasch P: Evolution of the Branchiopoda. Phylogeny and evolution of Crustacea. Edited by: Whittington HB, Rolfe WD. 1963, Harvard University, 145-162.
Trusheim F: Triopsiden (Crust. Phyll.) aus dem Keuper Frankens. Palaeontol Z. 1937, 19: 198-216.
Futuyma DJ: Evolutionsbiologie. 1990, Basel, Boston, Berlin: Birkhäuser
Korn M, Marrone F, Pérez-Bote JL, Machado M, Cristo M, da Fonseca C, Hundsdoerfer AK: Sister species within Triops cancriformis lineage (Crustacea, Notostraca). Zool Scripta. 2006, 35: 301-322. 10.1111/j.1463-6409.2006.00230.x.
Mantovani B, Cesari M, Scanabissi Sabelli F: Molecular taxonomy and phylogeny of the 'living fossil' lineages Triops and Lepidurus (Branchiopoda: Notostraca). Zool Scripta. 2004, 44: 367-374. 10.1111/j.0300-3256.2004.00155.x.
Boix D: Estructura i dinàmica de la comunitat animal aquàtica de l'estanyol temporani d'Espolla. PhD Dissertation. Institut d'Ecologia Aquàtica. 2002, Girona: Universitat de Girona
Simon E: Etude sur les Crustacés du sous-ordre des Phyllopodes. Ann Soc Entomol France. 1886, 6: 393-460.
Mathias P: Biologie des crustaces phyllopodes. 1937, Paris: Hermann & Cie
Akita M: Classification of Japanese tadpole shrimps. Zool Magazine. 1976, 85: 237-247.
Cesari M, Mularoni L, Scanabissi Sabelli F, Mantovani B: Characterization of dinucleotide microsatellite loci in the living fossil shrimp Triops cancriformis (Crustacea: Branchiopoda: Notostraca). Mol Ecol Notes. 2004, 4: 733-735. 10.1111/j.1471-8286.2004.00807.x.
Vandel A: La parthénogénèse géographique. Contribution à l'ètude biologiquie et cytologique de la parthénogénèse naturelle. Bull Biol France, Belgique. 1928, 12: 164-281.
Schoen I: Did Pleistocene glaciations shape genetic patterns of Europe ostracodes? A phylogeographic analysis of two species with asexual reproduction. Hydrobiologia. 2007, 575: 33-50. 10.1007/s10750-006-0276-z.
Knebelsberger T, Bohn H: Geographic parthenogenesis in the subaptera-group of Phyllodromica (Blattoptera, Blattellidae, Ectobiinae). Insect Syst Evol. 2003, 34: 427-452.
Ketmaier V, Mandatori R, De Matthaeis E, Mura G: Molecular systematics and phylogeography in the fairy shrimp (Tanymastix stagnalis) based on mitochondrial DNA. J Zool. 2005, 266: 1-10. 10.1017/S0952836904006247.
Ketmaier V, Zarattini P, De Matthaeis E, Cobolli M, Mura G: Intra-and inter-specific relationships in the six Italian species of the fairy shrimp genus Chirocephalus: combining allozyme and mtDNA data. J Zool Syst Evol Research. 2003, 41: 276-285. 10.1046/j.1439-0469.2003.00222.x.
de Meester L, Gómez A, Okamura B, Schwenk K: The monopolization hypothesis and the dispersal-gene flow paradox in aquatic organisms. Acta Oecol. 2002, 23: 121-135. 10.1016/S1146-609X(02)01145-1.
Steinfartz S, Veith M, Tautz D: Mitochondrial sequence analysis of Salamandra taxa suggests old splits of major lineages and postglacial recolonizations of Central Europe from distinct source populations of Salamandra salamandra. Mol Ecol. 2000, 9: 397-410. 10.1046/j.1365-294x.2000.00870.x.
Gómez A, Carvalho GR, Lunt DH: Phylogeography and regional endemism of a passively dispersing zooplankter: mtDNA variation of rotifer resting egg banks. Proc R Soc Lond B. 2000, 267: 2189-2197. 10.1098/rspb.2000.1268.
De Gelas K, De Meester L: Phylogeography of Daphnia magna in Europe. Mol Ecol. 2005, 14: 753-764. 10.1111/j.1365-294X.2004.02434.x.
Gómez A, Lunt DH: Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. Phylogeography in Southern European Refugia: Evolutionary Perspectives on the origins and conservation of European Biodiversity Chapter IV. Edited by: Weiss S, Ferrad D. 2007, Dordrecht, The Netherlands: Springer, 155-188.
Martínez-Solano I: Phylogeography of Iberian Discoglossus (Anura: Discoglossidae). J Zool Syst Evol Research. 2004, 42: 298-305. 10.1111/j.1439-0469.2004.00257.x.
Veith M, Mayer C, Samraoui B, Barroso DD, Bogaerts S: From Europe to Africa and vice versa: evidence for multiple intercontinental dispersal in ribbe salamanders (Genus Pleurodeles). J Biogeogr. 2004, 31: 159-171. 10.1111/j.1365-2699.2004.00957.x.
Harris DJ, Carranza S, Arnold EN, Pinho C, Ferrand N: Complex biogeographical distribution of genetic variation within Podarcis Wall lizards across the Strait of Gibraltar. J Biogeogr. 2002, 29: 1257-1262. 10.1046/j.1365-2699.2002.00744.x.
Carranza S, Harris DJ, Arnold EN, Batista V, Gonzalez de la Vega JP: Phylogeography of the lacertid lizard, Psammodromus algirus, in Iberia and across the Strait of Gibraltar. J Biogeogr. 2006, 33: 1279-1288. 10.1111/j.1365-2699.2006.01491.x.
Vargas JM, Real R, Guerrero JC: Biogeographical regions of the Iberian peninsula based on freshwater fish and amphibian distributions. Ecography. 1998, 21: 371-382. 10.1111/j.1600-0587.1998.tb00402.x.
King JL, Hanner R: Cryptic species in a "living fossil" lineage: Taxonomic and phylogenetic relationships within the genus Lepidurus (Crustacea: Notostraca) in North America. Mol Phylogenet Evol. 1998, 10: 23-36. 10.1006/mpev.1997.0470.
Hebert PDN, Ratnasingham S, deWaard JR: Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc R Soc Lond B (Suppl). 2003, 270: S96-S99. 10.1098/rsbl.2003.0025.
Ghigi A: Ricerche sui notostraci di Cirenaicae di altri paesi del mediterraneo. Atti Soc Ital Sci nat. 1921, 60: 161-188.
Gómez A, Carvalho GR: Sex, parthenogenesis and genetic structure of rotifers: microsatellite analysis of contemporary and resting egg bank populations. Mol Ecol. 2000, 9: 203-214. 10.1046/j.1365-294x.2000.00849.x.
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 Biotechnol. 1994, 3: 294-299.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiplesequence alignment through sequence weighting, position-specific gap penalties and weight matrixchoice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Farris JS, Källersjo M, Kluge A, Bult C: Testing significance of incongruence. Cladistics. 1994, 10: 315-319. 10.1111/j.1096-0031.1994.tb00181.x.
Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony. 1998
Keane TM, Naughton TJ, McInerney JO: ModelGenerator: amino acid and nucleotide substitution model selection, Version 0.6. 2004
Guindon S, Gascuel O: PHYML, version 2.4.4. 2003
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Hillis DM, Bull JJ: An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst Biol. 1993, 42: 182-192. 10.2307/2992540.
Wittingham LA, Silkas B, Winkler DW, Sheldon F: Phylogeny of tree swallow, Tachycineta (Aves: Hirundinidae), Bayesian analysis of mitochondrial DNA sequences. Mol Phylogenet Evol. 2002, 22: 430-441. 10.1006/mpev.2001.1073.
Pierrière G, Gouy M: WWW-Query: An online retrieval system for biological sequence banks. Biochemie. 1996, 78: 364-369. 10.1016/0300-9084(96)84768-7.
Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.
Fryer G: A new classification of the branchiopod Crustacea. Zool J Linnean Soc. 1987, 91: 357-383.
Wallossek D: The Upper Cambrian Rehbachiella and the phylogeny of Branchiopoda and Crustacea. Fossils and Strata. 1993, 32: 1-202.
Wallossek D: The Upper Cambrian Rehbachiella, its larval development, morphology and significance for the phylogeny of Branchiopoda and Crustacea. Hydrobiologia. 1995, 298: 1-13. 10.1007/BF00033797.
Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood appoach. Mol Biol Evo. 2002, 19: 101-109.
Knowlton N, Weigt LA: New dates and new rates for divergence across the Isthmus of Panama. Proc R Soc Lond B. 1998, 265: 2257-2263. 10.1098/rspb.1998.0568.
Schubart CD, Diesel R, Hedges SB: Rapid evolution to terrestrial life in Jamaican crabs. Nature. 1998, 393: 363-365. 10.1038/30724.
Finlayson C, Carrión JS: Rapid ecological turnover and its impact on Neanderthal and other human populations. Trends Ecol Evol. 2007, 22: 213-222. 10.1016/j.tree.2007.02.001.
Hempel-Zawitkowska J: Some remarks on morphological variability in Triops cancriformis (Bosc). Zool Poloniae. 1968, 18: 301-316.
Heidecke D, Neumann V: Zur Verbreitung und Ökologie von Triops cancriformis BOSC und Lepidurus apus L. in der DDR. Hercynia N F. 1987, 24: 166-173.
Gaschott O: Beobachtungen und Versuche an Triops cancriformis (Bosc). Zool Anz. 1928, 75: 267-280.
Kozubowsky A: Über den männlichen Apus cancriformis (L.). Z Naturgesch. 1857, 23: 312-318.
Braem F: Bemerkungen über die Gattung Apus. Z wissenschaftl Zool. 1893, 56: 165-187.
von Siebold CTE: Beiträge zur Parthenogenesis der Arthropoden. 1871, Leipzig: Verlag W. Engelmann
Abonyi S: The males of Apus cancriformis (Schäffer), described on the base of the specimens collected in the region of the lake Balaton biological station of Révfülöp. Arch Balaton. 1926, 1: 71-90.
Machado M, Cristo M, Reis J, da Fonseca C: Biological data on Triops cancriformis mauritanicus (Ghigi, 1921) and Cyzicus grubei (Simon, 1886) – Crustacea, Branchiopoda – in SW Portugal temporary ponds. Limnetica. 1999, 16: 1-7.
We are grateful to J. deWaard, who sent us unpublished T. australiensis COI and 16S sequences. We thank Drs. J. L. Pérez-Bote, J. Rueda Sevilla, F. Mezquita, J. K. Vanderkerkhove, D. Boix Sala, L. Griffin, A. Petrusek, L. Simon, M. Kawata, L. Forró, F. Suhling, L. Simon and M. Korn for kindly providing samples. D. Boix Sala provided us with unpublished data from his PhD thesis. D. H. Lunt, R. L. Hammond, M. Engelmann and two anonymous referees provided constructive comments and discussion in earlier versions of the manuscript. This research was supported by grants from NERC (NE/D007496/1), HERI and the German Environmental Foundation (DBU 20002/243). We further wish to thank English Nature, Scottish Heritage and the Sieland Foundation (Germany) for sample permissions.
TZ participated in the design of the study, carried out the molecular genetic studies, sequence alignment and sex-ratio analysis and drafted the manuscript. BH participated in the design of the study and the statistical analysis and helped to draft the manuscript. AG conceived the study, and participated in its design, the analysis and interpretation of data and helped to draft and revise the manuscript it critically for important intellectual content. All authors read and approved the final manuscript.