Recent evolution of alternative reproductive modes in the 'living fossil' Triops cancriformis
© Zierold et al; licensee BioMed Central Ltd. 2007
Received: 21 February 2007
Accepted: 13 September 2007
Published: 13 September 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
Triops cancriformis sample locations, population code, location of inferred subspecies, source of material, number of individuals sequenced for COI and ATPase and observed haplotypes (and numbers individuals per haplotype) per population, geographic are given based on the combined COI and ATPase analysis.
54°08'27" N, 003°30'58" E
53°34'10" N, 011°25'34" E
Döberitzer Heide, Germany
52°31'43" N, 013°00'05" E
51°48'12" N, 014°24'02" E
51°48'05" N, 014°23'37" E
51°47'47" N, 014°23'31" E
51°19'50" N, 014°18'30" E
51°19'49" N, 014°18'31" E
51°19'48" N, 014°18'12" E
51°19'46" N, 014°18'14" E
51°19'43" N, 014°18'23" E
Godshill pond, UK
50°55'37" N, 001°46'30" E
49°43'41" N, 008°25'55" E
Daxlander Au, Germany
49°01'10" N, 008°17'00" E
49°00'40" N, 008°16'25" E
48°59'30" N, 008°16'25" E
Morava, Czech Republic
48°10'38" N, 021°37'18" E
47°47'40" N, 016°52'31" E
47°39'01" N, 020°43'37" E
42°09'02" N, 002°45'60" E
Ares del Maestre, Spain
40°25'19" N, 000°04'13" W
El Puig, Spain
39°34'24" N, 000°16'43' E
Laguna de la Gitanilla, Spain
39°27'00" N, 006°15'54" W
Ullal de Baldovi, Spain
39°14'55" N, 000°19'03" E
Gorgo di Baglio Cofano Sicily, Italy
38°06'11" N, 012°40'39" E
Yamagata Prefercture, Japan
38°02'58" N, 140°10'48" E
32°17'06" N, 008°19'22" W
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%).
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.
Sequence divergence and approximation of divergence times
Mean corrected percent sequence divergences between Triops cancriformis clade/population comparisons and estimated divergence times based on corrected COI distance. Distances are given for the combined COI-ATPase fragment as well as separately for COI and ATPase sequences. Because no crustacean molecular clock calibration factors are available for ATPase we used the COI corrected distance for approximations of divergence times. For population codes see Table 1, for clade definition see Figure 2.
Corrected % sequence divergence
Estimated divergence times (mya)
1.4% sequence divergence per mya
2.3% sequence divergence per mya
clade A vs. clade B
clade A and B vs. C
ARE vs. YOI
ARE&YOI vs. EXT
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
Sex ratio and inferred reproductive mode in European Triops cancriformis. The number (n) of individual males and females/hermaphrodites and the male proportion per pond are given. The results of the Chi-Square test to infer deviations from a null hypothesis of equal sex ratio and the inferred type of reproductive mode are provided (see text). In the populations investigated here, the collection year is indicated under the Reference heading.
prop. of males
p value (Chi-Square)
Inferred reproductive mode
own data (2004–2006)
Godshill pond, United Kingdom
this paper (2004–2006)
Königswartha (KOE25), Germany
this paper (2006)
Königswartha (KOE12), Germany
this paper (2004 to 2006)
Königswartha (KOE27), Germany
this paper (2006)
unknown location, Hungary
Lake Balaton, Hungary
Lake Balaton, Hungary
this paper (2006)
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.
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.
- Otto SP, Lenormand T: Resolving the paradox of sex and recombination. Nat Rev Genet. 2002, 3: 252-261. 10.1038/nrg761.View ArticlePubMedGoogle Scholar
- Bell GD: The Masterpiece of Nature: The Evolution and Genetics of Sexuality. 1982, Berkley and Los Angeles: University of California PressGoogle Scholar
- Baker HG: Self-compatibility and establishment after "long-distance" dispersal. Evolution. 1955, 9: 347-348. 10.2307/2405656.View ArticleGoogle Scholar
- Baker HG: Support for Baker's Law-as a rule. Evolution. 1967, 21: 853-856. 10.2307/2406780.View ArticleGoogle Scholar
- Peck JR, Yearsley JM, Waxman D: Explaining the geographic distributions of sexual and asexual populations. Nature. 1998, 391: 889-892. 10.1038/36099.View ArticleGoogle Scholar
- Kearney MR: Hybridization, glaciation and geographical parthenogenesis. Trends Ecol Evol. 2005, 20: 495-502. 10.1016/j.tree.2005.06.005.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Sassaman C: Sex ratio variation in female-biased populations of Notostracans. Hydrobiologia. 1991, 212: 169-179. 10.1007/BF00025998.View ArticleGoogle Scholar
- Zaffagnini F, Trentini M: The distribution and reproduction of Triops cancriformis (Bosc) in Europe (Crustacea Notostraca). Monitore Zool Ital (NS). 1980, 14: 1-8.Google Scholar
- Pannel JR: The evolution and maintenance of androdioecy. Ann Rev Ecol Syst. 2002, 33: 397-425. 10.1146/annurev.ecolsys.33.010802.150419.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- Longhurst AR: Evolution in the Notostraca. Evolution. 1955, 9: 84-86. 10.2307/2405360.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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 thesisGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Frank PW: Conchostraca. Palaeogeogr Palaeoclimatol Palaeoecol. 1986, 62: 399-403. 10.1016/0031-0182(88)90064-8.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Bohonak AJ, Whiteman HH: Dispersal of the fairy shrimp Branchinecta coloradensis (Anostraca): Effects of hydroperiod and salamanders. Limnol Oceanogr. 1999, 44: 487-493.View ArticleGoogle Scholar
- Longhurst AR: A Review of the Notostraca. Bull Br Mus (Nat Hist) Zool. 1955, 3: 1-57.View ArticleGoogle Scholar
- Tasch P: Evolution of the Branchiopoda. Phylogeny and evolution of Crustacea. Edited by: Whittington HB, Rolfe WD. 1963, Harvard University, 145-162.Google Scholar
- Trusheim F: Triopsiden (Crust. Phyll.) aus dem Keuper Frankens. Palaeontol Z. 1937, 19: 198-216.View ArticleGoogle Scholar
- Futuyma DJ: Evolutionsbiologie. 1990, Basel, Boston, Berlin: BirkhäuserView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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 GironaGoogle Scholar
- Simon E: Etude sur les Crustacés du sous-ordre des Phyllopodes. Ann Soc Entomol France. 1886, 6: 393-460.Google Scholar
- Mathias P: Biologie des crustaces phyllopodes. 1937, Paris: Hermann & CieGoogle Scholar
- Akita M: Classification of Japanese tadpole shrimps. Zool Magazine. 1976, 85: 237-247.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- Knebelsberger T, Bohn H: Geographic parthenogenesis in the subaptera-group of Phyllodromica (Blattoptera, Blattellidae, Ectobiinae). Insect Syst Evol. 2003, 34: 427-452.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Ghigi A: Ricerche sui notostraci di Cirenaicae di altri paesi del mediterraneo. Atti Soc Ital Sci nat. 1921, 60: 161-188.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony. 1998Google Scholar
- Keane TM, Naughton TJ, McInerney JO: ModelGenerator: amino acid and nucleotide substitution model selection, Version 0.6. 2004Google Scholar
- Guindon S, Gascuel O: PHYML, version 2.4.4. 2003Google Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Fryer G: A new classification of the branchiopod Crustacea. Zool J Linnean Soc. 1987, 91: 357-383.View ArticleGoogle Scholar
- Wallossek D: The Upper Cambrian Rehbachiella and the phylogeny of Branchiopoda and Crustacea. Fossils and Strata. 1993, 32: 1-202.Google Scholar
- 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.View ArticleGoogle Scholar
- Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood appoach. Mol Biol Evo. 2002, 19: 101-109.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Schubart CD, Diesel R, Hedges SB: Rapid evolution to terrestrial life in Jamaican crabs. Nature. 1998, 393: 363-365. 10.1038/30724.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Hempel-Zawitkowska J: Some remarks on morphological variability in Triops cancriformis (Bosc). Zool Poloniae. 1968, 18: 301-316.Google Scholar
- 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.Google Scholar
- Gaschott O: Beobachtungen und Versuche an Triops cancriformis (Bosc). Zool Anz. 1928, 75: 267-280.Google Scholar
- Kozubowsky A: Über den männlichen Apus cancriformis (L.). Z Naturgesch. 1857, 23: 312-318.Google Scholar
- Braem F: Bemerkungen über die Gattung Apus. Z wissenschaftl Zool. 1893, 56: 165-187.Google Scholar
- von Siebold CTE: Beiträge zur Parthenogenesis der Arthropoden. 1871, Leipzig: Verlag W. EngelmannView ArticleGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.