Scholarly article on topic 'High lability of sexual system over 250 million years of evolution in morphologically conservative tadpole shrimps'

High lability of sexual system over 250 million years of evolution in morphologically conservative tadpole shrimps Academic research paper on "Biological sciences"

Share paper
Academic journal
BMC Evol Biol
OECD Field of science

Academic research paper on topic "High lability of sexual system over 250 million years of evolution in morphologically conservative tadpole shrimps"

Evolutionary Biology

High lability of sexual system over 250 million years of evolution in morphologically conservative tadpole shrimps

Mathers et al.

Bio Med Central

Mathers et al. BMC Evolutionary Biology 2013,13:30

Evolutionary Biology


High lability of sexual system over 250 million years of evolution in morphologically conservative tadpole shrimps

Thomas C Mathers1*, Robert L Hammond2, Ronald A Jenner3, Thorid Zierold1,4, Bernd Hänfling1 and Africa Gomez1


Background: Sexual system is a key factor affecting the genetic diversity, population structure, genome structure and the evolutionary potential of species. The sexual system androdioecy - where males and hermaphrodites coexist in populations - is extremely rare, yet is found in three crustacean groups, barnacles, a genus of clam shrimps Eulimnadia, and in the order Notostraca, the tadpole shrimps. In the ancient crustacean order Notostraca, high morphological conservatism contrasts with a wide diversity of sexual systems, including androdioecy. An understanding of the evolution of sexual systems in this group has been hampered by poor phylogenetic resolution and confounded by the widespread occurrence of cryptic species. Here we use a multigene supermatrix for 30 taxa to produce a comprehensive phylogenetic reconstruction of Notostraca. Based on this phylogenetic reconstruction we use character mapping techniques to investigate the evolution of sexual systems. We also tested the hypothesis that reproductive assurance has driven the evolution of androdioecy in Notostraca. Results: Character mapping analysis showed that sexual system is an extremely flexible trait within Notostraca, with repeated shifts between gonochorism and androdioecy, the latter having evolved a minimum of five times. In agreement with the reproductive assurance hypothesis androdioecious notostracans are found at significantly higher latitudes than gonochoric ones indicating that post glacial re-colonisation may have selected for the higher colonisation ability conferred by androdioecy.

Conclusions: In contrast to their conserved morphology, sexual system in Notostraca is highly labile and the rare reproductive mode androdioecy has evolved repeatedly within the order. Furthermore, we conclude that this lability of sexual system has been maintained for at least 250 million years and may have contributed to the long term evolutionary persistence of Notostraca. Our results further our understanding of the evolution of androdioecy and indicate that reproductive assurance is a recurrent theme involved in the evolution of this sexual system.

Keywords: Androdioecy, Character evolution, Sexual system, Notostraca, Phylogeny


Plants and animals have evolved a great diversity of sexual systems that range from the extremes of obligatory self-fertilisation to complete outcrossing. Transitions between these sexual systems have long fascinated biologists due to the impacts they have on key parameters such as inbreeding depression, genetic diversity, population structure, genome structure and the evolutionary potential of species [1-7]. Transitions between sexual

* Correspondence:

1Schoolof Biological, Biomedicaland EnvironmentalSciences, University of Hull, HullHU6 7RX, UK

Fulllist of author information is available at the end of the article

Bio Med Central

systems often present tradeoffs between short and long term selective advantages and can have significant connotations for the long-term viability of species. For example, selection for reproductive assurance and colonisation advantage due to mate limitation during range expansions, or as a result of high population turnover in metapopulations, can drive transitions to self-fertilisation strategies [8-10]. These transitions occur despite the deleterious effects of self-fertilisation, which include inbreeding depression, reduction in effective recombination rates and reduction in effective population size [11].

Transitions to androdioecy (AD) - a sexual system where males and hermaphrodites co-occur in varying

© 2013 Mathers et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, anc reproduction in any medium, provided the original work is properly cited.

frequencies in populations, with different levels of self-fertilisation and outcrossing - are extremely rare in plants and animals [12-16]. In animals, AD has only been described in five groups, rhabditid nematodes, the killi-fish Kryptolebias marmoratus and three crustacean groups; barnacles, a genus of clam shrimps Eulimnadia, and in the order Notostraca, the tadpole shrimps [17,18]. AD can evolve either through the invasion of males into hermaphrodite only populations, as in barnacles [18,19], or through the replacement of females with hermaphrodites in gonochoric populations (where males and females are found in approximate equality), as in the plants Mercurialis annua [20,21] and Datisca glomerata [22]. As models to describe the evolution and maintenance of AD only predict its evolution under stringent conditions, AD has historically been considered an unstable, transitional sexual system between gonochorism and hermaphroditism (or vice versa) [12-14,23,24]. This view is borne out by the scarcity of AD in nature [12,17], although recent research in the branchiopod Eulimnadia has revealed an unexpected stability of androdioecy [25].

Notostraca, or tadpole shrimps, is a small order of branchiopod crustaceans characterised by a high level of morphological stasis. Fossils dating back as far as the Triassic are almost indistinguishable from contemporary species leading them to be referred to as 'living fossils' [26-30]. In contrast, Notostraca has diverse sexual systems, including gonochorism, self-fertile hermaphrodit-ism and AD, with variation occurring on both an interspecific and intraspecific level [31,32]. Remarkably, AD is found in species from both notostracan genera, Triops and Lepidurus, suggesting that transitions in reproductive system might have evolved repeatedly in the order. Despite this, the evolutionary history of reproductive systems in Notostraca is unknown due to the lack of a resolved phylogeny [33,34], and the poor knowledge of the diversity of the group, partly due to the widespread presence of cryptic species [32,35-38]. Gono-chorism has been hypothesized to be the ancestral state in the group, and the evolution of self-fertile hermaphro-ditism and AD has been linked to reproductive assurance in the context of range expansions, possibly after glacial retreat [8,31,39,40], although this has never been explicitly tested.

Here we combine newly generated and GenBank sequence data to assess Notostraca taxonomic diversity, identifying considerable cryptic diversity, and employ a multigene phylogenetic approach to create a well-supported, global phylogeny of Notostraca. Information on sexual system was compiled and Maximum Parsimony (MP) and model-based Maximum Likelihood (ML) character mapping approaches were used on the phylogeny to investigate sexual system evolution across the order. We also tested the hypothesis that reproductive

assurance has driven the evolution of self-fertilisation across Notostraca [8,39]. Taxa found at higher latitudes are likely to have experienced bouts of colonisation during post glacial range expansions, which would select for AD/hermaphroditism. We therefore compared the latitudes that AD/hermaphroditic and gonochoric taxa are found using a phylogenetic t-test. Our analyses reveal high levels of reproductive lability with frequent transitions occurring to and from androdioecy. Furthermore, this flexibility is conserved across Notostraca, and may have been maintained for at least 250 million years. Additionally, AD/hermaphroditic taxa are found at significantly higher latitudes than gonochoric ones suggesting that colonisation advantage through reproductive assurance is likely to be involved in transitions between sexual systems in Notostraca.


Delimitation of significant taxonomic units

Notostraca is known to contain cryptic species complexes e.g. [35,37] so we first used a cytochrome oxidase subunit one (COI)-based barcoding approach to identify significant taxonomic units (STUs) for inclusion in our phylogeny prior to the multigene analysis. Including available GenBank data and 12 newly generated sequences for this study, 243 Notostraca COI sequences were aligned. We applied a generalized mixed Yule co-alescent (GMYC) model to identify independently evolving clusters in our COI dataset, which correspond to STUs. The GMYC model identified 26 STUs (Additional file 1: Figure S1). Uncorrected mean genetic distances in COI between STUs ranged from 2.3% to 24.3%. Four Notostraca lineages did not have COI data available, but are represented by other genes used in our multigene phylogenetic analysis; T. gadensis, T. cf. granarius (Tunisia) L. bilobatus and L. cryptus. As the species status of these lineages has been confirmed in regional studies of cryptic diversity in Notostraca [35-37,41] they were included as additional STUs for the multigene phylogenetic analysis. In total we recognise 30 STUs within Notostraca.

Notostracan phylogeny

Phylogeny, based on a concatenated supermatrix of 110 sequences (54 of which were newly generated for this study) from three mitochondrial genes and four nuclear genes for 30 STUs, was inferred by ML and Bayesian Markov-chain Monte Carlo methods. Both methods of phylogenetic reconstruction gave congruent topologies, with most branches having high levels of support (Figure 1). The two recognised notostracan genera, Triops and Lepidurus, formed highly supported clades. Within Triops, four main monophyletic lineages with a strong geographic association (Australian, N American, Palearctic and African/Asian respectively)

100 87

-T. cf. australiensis sp. 3

-T. cf. australiensis sp. 2

— australiensis sp. 1

-T. cf. australiensis sp. A

-T. cf. australiensis sp. A

961-T. cf. longicaudatus sp

LÍ7Í-T. newberryi

1 '-T. cf longicaudatus sp. 1

- T. cf granarius (Russia)

1-T. cf granarius (Namibia) granarius (Japan) -T. cf granarius (Russia)

- T. vicentinus

- T. emeritensis

- T. mauritanicus

- T. simplex

- T. cf mauritanicus (E. Spain) -T. gadensis

T. baeticus T. cancriformis L. lemmoni L. cryptus

L. bilobatus


-L. cf. couesii (Apulia) ~L. couesii couesii (Sardinia) -L. arcticus

-L. apus — L. lubbocki

Figure 1 Phylogenetic relationships in Notostraca based on a multigene supermatrix. Tree topology shows the best scoring ML tree from the RAxML analysis. Values above nodes show Bayesian posterior probabilities and values below nodes, in bold, show bootstrap support (1000 replicates). Branches with bootstrap support lower than 50 are collapsed. The Leptestheria outgroup was removed after rooting.

were highly supported. The analysis supported a close relationship between the T. australiensis complex (Australia) and the T. longicaudatus complex (North America), and a sister relationship of these with the T. granarius complex (Asia/Africa). The T. cancriformis/mauritanicus complex (Palearctic) appears as the sister group to the rest of Triops. Within Lepidurus, L. lubbocki (Mediterranean), and L. apus sensu stricto (N European) have long branches and are sister species to the rest of Lepidurus. Four North American species [35,41] L. packardi, L. cryptus, L. biloba-tus and L. lemmoni, with narrowly endemic, mostly allo-patric distributions in western North America, form a well supported group. L. arcticus, a circumpolar species from Arctic and Subarctic regions, forms a sister relationship with a clade containing L. couesii, which forms a widely distributed species complex.

Sexual system assignment

Our literature review identified sexual system data for 22 STUs; 18 of these have at least one barcoded population (either COI, 12S or 16S). In the few cases where an STU does not have a barcode sequence for the same population for which sexual system data is derived (L. arcticus, L. couesii, L. lemmoni, L. packardi) these species have been well studied, no further cryptic diversity has been identified and the populations used for sexual system inference fall inside known species ranges (see [35,41,42]). We found two polymorphic STUs; T. cancriformis and T. cf. longicaudatus sp.2,

which include both androdioecious and gonochoric populations, and four androdioecious STUs; T. newberryi, T. cf. australiensis sp. B, L. apus (sensu stricto) and L. arcticus (Table 1 and Additional file 1, Table S6). Sex ratio of populations assigned AD ranged from 0% to 27% males. In addition, either histological data or reproduction in isolation data, or both, confirmed the presence of anatomical hermaphrodites. Fourteen taxa were considered gonochoric on the basis of histology, lack of reproduction in isolation and sex ratio. The quality and quantity of information on sexual system was heterogeneous amongst gonochoric taxa, with actual sex ratios tending to be underreported and histo-logical studies lacking for many taxa. Eight taxa could not be assigned to a sexual system due to an absence of data or equivocal sex ratio.

Sexual system evolution

MP reconstruction of ancestral character states infers that gonochorism is the ancestral state of Notostraca (Figure 2). Furthermore, AD appears to have multiple origins in Notostraca having evolved three times in Triops and twice in Lepidurus. Sexual system is highly flexible across Notostraca and varies even between closely related species (T. cf. australiensis A vs. T. cf. australiensis B; T. newberryi vs. T. cf. longicaudatus sp.1) or shows intraspecific variation (T. cancriformis; T. cf. longi-caudatus sp.2). Model based ML methods showed that a two-parameter model, which allows distinct transition

Table 1 Notostraca sexual system information

STU Sex ratio Reproduction in isolation Ovotestis Sexual system References

T. baeticus > 45 Gonochoric [36]

T. cancriformis 0 - 53* Yes** Yes** Gonochoric/AD [26,31,40,43,44]

T. cf. australiensis sp. 1 No data

T. cf. australiensis sp. 2 No data

T. cf. australiensis sp. 3 No data

T. cf. australiensis sp. A > 45 No Gonochoric [45]

T. cf. australiensis sp. B < 30 Yes AD [45]

T. emeritensis >45 Gonochoric [36]

T. gadensis 36 Equivocal [36]

T. cf. granarius (Japan) > 45 No Gonochoric [26,44,46,47]

T. cf. granarius (Namibia) Even*** No Gonochoric [26,37,44]

T. cf. granarius (Tunisia) Even*** No Gonochoric [26,37,44]

T. cf. granarius (Russia) No data

T. cf. longicaudatus sp. 1 > 45 Gonochoric [32]

T. cf. longicaudatus sp.2 0 - 68* Yes** Gonochoric/AD [32,48,49]

T. cf. mauritanicus (E Spain) No data

T. mauritanicus No data

T. newberryi < 30 Yes AD [32,48,50]

T. simplex > 45 Gonochoric [38]

T. vicentinus > 45 Gonochoric [36,51]

L. apus < 45 Yes AD [26,40,44,52-55]

L. arcticus <45 Yes AD [26,40,44,54,56-58]

L. bilobatus 35 Equivocal

L. cf. couesii (Apulia) No data

L. cf. couesii (Sardinia) >45 Gonochoric [59]

L. couesii (Canada) >45 Gonochoric [60]

L. cryptus No data

L. lemmoni >45 Gonochoric [56,61]

L. lubbocki >45 No Gonochoric [62,63]

L. packardi >45 Gonochoric [64]

Sex ratio (percent male), ability of ovisac bearing individuals to reproduce in isolation and the presence or absence of ovotestis in ovisac bearing individuals is shown. See Table S5 for detailed information.

* Depending on population. ** Only in populations where sex ratio is < 30%. *** Reported gonochoric with even sex ratio, exact numbers not given.

rates for AD to gonochorism and from gonochorism to AD, was a significantly better fit for the data than a model where both transitions have an equal rate (Table 2), or models where transitions were restricted to one direction, either from gonochorism to AD or AD to gonochorism, indicating that in Notostraca changes in sexual system could be bidirectional. Overall, the ML model suggests that transition rates between sexual systems were high and in particular, transitions from AD to gonochorism were more than three times higher than transitions from gonochorism to AD. This result is in striking contrast to the MP results which suggested repeated evolution in the

opposite direction, to AD from gonochorism. This indicates that, once evolved, AD may be unstable and likely to revert back to gonochorism. The high rates of change across the order meant that, unlike for the MP analysis where a minimum number of transitions is inferred, ancestral sexual systems for all nodes were equally likely to be either gonochoric or AD. Virtually identical results were achieved using an ultrametric phylogeny constructed in BEAST v1.7.4 [65] with a lognormal relaxed molecular clock both from the full dataset and with a reduced dataset containing only the mitochondrial genes COI, 12S and 16S (see Additional file 1).

T. cf. australiensis sp. 3 T. cf. australiensis sp. 2 T. cf. australiensis sp. 1 T. cf. australiensis sp. B T. cf. australiensis sp. A T. cf. longicaudatus sp. 2 T. newberryi T. cf. longicaudatus sp. 1 T. cf. granarius (Tunisia) T. cf. granarius (Namibia) T. cf. granarius (Japan) T. cf. granarius (Russia) T. vicentinus T. emeritensis T. mauritanicus T. simplex

T. cf. mauritanicus (E. Spain) T. gadensis T. baeticus T. cancriformis L. lemmoni L. cryptus L. bilobatus L. packardi L. cf. couesii (Apulia) L. couesii

L. cf. couesii (Sardinia) L. arcticus L. apus L. lubbocki

Figure 2 Maximum parsimony ancestral character state reconstruction of sexual systems in Notostraca. Sexualsystem is mapped onto the best scoring ML tree, and indicated by the colour of the square that precedes the taxon names. Blue = gonochoric, Orange = Androdioecy. STUs without squares preceding the taxon name do not have sexualsystem information or sexualsystem is equivocal. Bi-coloured squares indicate that both sexualsystems are found. Branches are coloured according to MP reconstruction of ancestralsexualsystems based on ar unordered modelwith grey branches denoting an equivocalancestralstate assignment. The Leptestheria outgroup was removed after rooting and left uncoded for sexualsystem.

Testing the reproductive assurance hypothesis

The mean latitude of STUs where AD is found was significantly higher than for gonochoric ones (phylogenetic t-test, t = 2.922, p = 0.010, d.f. = 17) with an average latitude of 49.329 compared to 37.256. ML value for \ was estimated to be 0.508 and its inclusion significantly improved the

Table 2 Comparison of models of sexual system evolution in Notostraca

Model lnL q ga qAG P

Unrestricted -7.229 23.729 110.640 -

Equalrates -10.262 5.730 5.730 0.0138

AD to gonochorism only -16.855 0.000 6.673 < 0.0001

Gonochorism to AD only -12.0864 1.755 0.000 0.0018

Models of transitions in sexual system across the best scoring ML estimation of Notostraca phylogeny compared using the ML implementation of BayesMultistate. lnL = log-likelihood of model, q^ = transition rate from gonochorism to androdioecy, qAG = transition rate from androdioecy to gonochorism, p = p-value for D with 1 degree of freedom comparing the restricted models to the unrestricted model.

likelihood of the model (D, p = 0.028, d.f. = 1) indicating that the latitude of STUs has moderate phylogenetic signal.


Sexual system evolution

Our analyses reveal that sexual system is a highly labile trait within Notostraca. The rare sexual system AD has evolved multiple times in both notostracan genera, with MP indicating at least five independent origins across the whole group (Figure 2). This lability is also supported by the model based ML analysis which infers high transition rates between sexual systems. Unlike the MP analysis, the ML analysis did not resolve the ancestral sexual system for any nodes. This suggests that MP may give an oversimplified reconstruction of the evolutionary history of sexual system in Notostraca and highlights the importance of incorporating branch lengths in ancestral state reconstructions [66]. The ML analysis also contrasts with MP by inferring that transitions

between gonochorism and AD can occur in both directions as is the case for pedunculate barnacles [18], rather than in a unidirectional manner as implied by MP. As there are no known biological constraints to transitions between sexual systems in either direction there is no reason to discredit the model of sexual system evolution inferred by either method. For example, a transition from AD to gonochorism could occur as a result of selection for outcrossing over selfing, followed by the loss of testis lobes in hermaphrodites. In conclusion, although it is not possible to infer the history of sexual system change in Notostraca with current data or methods of ancestral state reconstruction, our results do display a consistent pattern of high lability of sexual system across Notostraca.

Although flexibility in sexual system has previously been shown in T. cancriformis, where a recent transition from gonochorism to androdioecy has occurred [31,43], our results demonstrate that reproductive flexibility is a general feature of Notostraca as a whole. Given the age of the order - well-preserved notostracan fossils date back to the Carboniferous [67] and Triops and Lepidurus are known from the Permian and Triassic [28,29,68], dating the split in the two genera to at least ~250 million years ago - it can be inferred that reproductive lability has been maintained for at least 250 million years within tadpole shrimps and may have aided their long term evolutionary persistence. Such lability of sexual system in Notostraca contrasts strongly with the dynamics of sexual system evolution in the clam shrimps of the genus Eulimnadia, the other branchiopod crustacean where AD is found. In Eulimnadia, AD evolved once and has persisted for at least 24 million years, passing through many speciation events [25,69-71]. The contrast in sexual system evolution between both taxa is striking considering that tadpole shrimps and clam shrimps inhabit similar ecological niches [26,72], in some cases occupying the same pool [73,74], and have similar life histories, producing long-lived dormant cysts that survive during adverse periods and also provide the dispersal stage [75-77]. They also have a similar genetic mechanism of sex determination where males are recessive to hermaphrodites in AD populations [78-80]. Finally, in both groups, hermaphrodites are morphologically derived from obligatory out-crossing females, and can either self-fertilise or outcross with males but, unlike in plants or barnacles, are unable to fertilise each other [26,70]. Within Branchiopoda therefore, superficially similar ecological backgrounds and genetic mechanisms of sex determination have resulted in two very different, but equally unusual, evolutionary outcomes for sexual system; stable and conserved in Eulimnadia and labile and dynamic in Notostraca. In contrast, the flexibility of sexual system in notostracans resembles that of barnacles, where AD has evolved on

several occasions in response to mate limitation, although in this system AD has evolved from a hermaphroditic ancestor rather than a gonochoric one [18,81]. Modelling has explained the persistence of AD in Eulim-nadia and highlights the importance of the presence of males for limiting inbreeding depression and that a proportion of progeny produced by selfing - those that are homogametic - have reduced fitness due to the expression of sex-linked genetic load [82,83]. In Notostraca, however, these key parameters, along with other factors that contribute to the evolution and maintenance of AD, such as male-hermaphrodite encounter rates and fecundity, have not been investigated. Further research is therefore required to ascertain whether the dynamics of sexual system evolution in Notostraca necessitate the reevaluation of current models for the evolution of AD in animals. A lack of phylogenetic signal regarding sexual systems in Notostraca suggests that switches in sexual system occur in response to ecological factors. Notos-traca share similarities in life history with many plants [8], particularly those that exist in highly dynamic metapopulations, where colonisation of new habitats is a key feature of survival causing selection for reproductive assurance [84-86]. In many cases selection for reproductive assurance has caused a breakdown in plant self-incompatibility systems and the evolution of self-fertile hermaphrodites, which make optimal pioneer genotypes [87]. If reproductive assurance drives the evolution of self-fertile hermaphroditism and hence AD in Notostraca, as has been hypothesised [39,40], we would expect AD STUs to occur in areas which have recently become available for colonisation. These predictions appear to be met in T. cancriformis, where hermaphroditic and androdioecious populations occur in previously glaciated areas whereas known gonochoric populations occur in what were previously unglaciated refugia [31]. Furthermore, our results support the reproductive assurance hypothesis for the whole of Notostraca as our phylogenet-ically informed analysis reveals that AD notostracan STUs occur at significantly higher latitudes than gonochoric ones. These tests indicate that AD STUs occur in areas where habitat disruption through glacial cycles is more likely, providing further support for the role of colonisation advantage and reproductive assurance in a metapopulation as drivers for the evolution of AD.

Highly fluctuating population densities, which result in mate limitation, could also provide an advantage to lineages containing hermaphrodites through reproductive assurance [13,86]. In the case of Branchiopoda, the role of low population densities and mate limitation in the evolution of sexual systems is still poorly understood [82]. Ephemeral ponds are often very small [88] with strong inter-annual fluctuations in branchiopod population density [25], which could provide a context in

which self-fertilising hermaphrodites would enjoy an increased fitness relative to females, driving the evolution of AD. Further research on Notostraca species distribution, genetic diversity, metapopulation dynamics and phylogeography will help to understand the underlying factors behind changes in sexual systems in this group.

Taxonomic implications of notostracan phylogeny

Our results resolve the phylogenetic relationships of Notostraca, supporting the previously disputed [33,34] monophyly of Triops. Although the four main lineages in Triops coincide with the four species recognised by Longhurst [26] in the last revision of Notostraca, our analyses support previous work revealing that each of these lineages is made up of cryptic species. Each main lineage has maintained a remarkably stable, mostly allopat-ric geographic distribution (only Japan has representatives of the T. cancriformis, T. cf. granarius and T. cf. longicau-datus complexes, and N Africa shares both T. mauritani-cus and T. cf. granarius). In Lepidurus, L. lubbocki was found to be the sister species to the rest of the genus and as suggested by Mantovani et al. [89] full species status is warranted. Further cryptic diversity was also identified in the L. couesii complex and given that the type locality for L. couesii is in N America [41] and the level of divergence with Apulian (S. Italy) and Sardinian lineages, we propose these latter lineages are new species. Overall, although highlighting the need for further research into Notostraca taxonomic diversity, our phylogeny represents an excellent framework for the study of evolutionary processes within the order.


Our analyses show that sexual systems are extremely flexible in Notostraca, with repeated switches between gono-chorism and androdioecy possibly driven by postglacial range expansions. This unconstrained, labile pattern, strongly contrasts with the single origin of androdioecy in the genus Eulimnadia despite the similarity of their habitats and life histories and resembles the pattern found in barnacles. Despite the idiosyncratic evolution in these crustacean taxa, reproductive assurance in the face of fluctuating population sizes, habitat turnover or climate changes, appears to be a recurrent theme in the evolution of androdioecy. Flexibility in sexual system evolution has been maintained throughout the evolutionary history of Notostraca (over 250 my) and given the extreme morphological and life history conservatism in the group, could have facilitated their evolutionary persistence.


Sampling and sequencing

We produced de novo sequence data from 12 taxa from the two notostracan genera, Triops and Lepidurus [26].

Samples consisted of either sediments containing resting eggs or wild caught individuals preserved in 100% etha-nol (Additional file 1: Table S1). Total genomic DNA was extracted from ethanol-preserved individuals using a Qiagen DNeasy Blood and Tissue Kit (Qiagen, Valencia, CA) or directly from individual resting eggs using a modified 'HotSHOT' procedure [90]. DNA sequences were generated for three mitochondrial gene fragments, COI, 12S rDNA and 16S rDNA, and four nuclear gene fragments; elongation factor 1 alpha, glycogen synthase, RNA polymerase II and 28S rDNA. We used primer pairs known to amplify across Notostraca for the mitochon-drial and ribosomal genes [34,37,91] and designed new primers for the nuclear protein coding genes based on alignments of available sequences from notostracans and other branchiopods using PriFi [92] (see Additional file 1: Table S2 for primer sequences and optimised reaction conditions). Reactions were carried out in a final volume of 50 ^l containing 2 ^l of template DNA, 200 ^M of each primer, 200 ^M of each nucleotide, 0.01 U of Bio-Taq DNA polymerase (Bioline), 1x NH4 buffer (Bioline) and 2-3 mM MgCl2. Amplified fragments were purified and sequenced for both forward and reverse strands by Macrogen using an ABI 3730xl DNA Analyser (Macrogen Inc, Seoul, Korea). Sequences were manually edited using CodonCode Aligner v3.5 (CodonCode Corporation, Dedham, MA) with consensus sequences produced for each forward and reverse pair.

GMYC model based species delimitation

Available Notostraca COI sequences were downloaded from GenBank (Additional file 1: Table S3) and aligned with our newly generated sequences in MEGA 5 [93] using MUSCLE [94] with default parameters. We applied a generalized mixed Yule coalescent (GMYC) model [95] to identify independently evolving clusters in our COI dataset, which correspond to STUs. First, we created an ultrametric phylogeny based on our COI alignment using BEAST v1.6.2 [96]. The phylogenetic analysis was run for 600,000 iterations with trees printed every 1,000 iterations and the first 100,000 iterations removed as burnin. A GTR + r nucleotide substitution model was used with a strict molecular clock with the rate fixed to 1. From this the ultrametric maximum clade credibility consensus tree was constructed. The GMYC analysis was performed in R v2.14.1 (R Development Core Team, 2011) with the package splits v1.0-11 ( splits/). Clusters defined by the GMYC analysis were then assigned an STU I.D. based on the geographic location and species assignment of the accessions they contained. Un-corrected mean genetic distances in COI between STUs were calculated in MEGA 5 [93] with all positions containing gaps or missing data removed.

Alignment and supermatrix construction

Single sequences for each gene for each STU identified were selected (where available) for inclusion in our multigene phylogenetic analysis (Additional file 1: Table S4). Sequences generated in this study were used preferentially but, where only GenBank sequences were available, records were checked to confirm that samples were from the same or close geographic location to samples used for STU identification.

The final alignment of each nuclear protein coding gene and the mitochondrial gene COI was carried out in MEGA 5 [93] with MUSCLE [94] using default parameters. The ribosomal genes were aligned based on secondary structure information using RNAsalsa v0.8.1 [97] with Apis mellifera structural data used as a constraint. Weakened constraint settings (S1, S2 and S3 =0.51) were used to preserve structural information as described by Letsch and Kjer [98]. To confirm that the individual alignments were suitable for concatenation, phylogenetic congruence was tested with Concaterpillar v1.4 [99,100] using the GTR model and an a-level cut off of 0.05. No significant phylogenetic incongruence was identified (p = 0.55) and so all genes were concatenated using FASconCAT v1.0 [101]. The final supermatrix contained 5253 positions with 54% missing data. This number represents the overall missing data, not including indels, in the supermatrix alignment. It reflects the fact that most taxa retrieved from GenBank do not have coverage for all the genes used in this study and we could not obtain sequences for some genes for a few of our samples.

Phylogenetic analysis

Phylogeny, based on the concatenated supermatrix, was inferred by ML and Bayesian Markov-chain Monte Carlo methods. We estimated the ML tree with RAxML using RAxMLHPC-PTHREADS v7.0.4 [100], treating each gene as an individual partition. An initial ML search using GTR + r was performed with 100 iterations to identify the best scoring ML tree. 1000 Bootstrap replicates were then conducted using GTR + r and drawn onto this best scoring ML tree. Bayesian phylogenetic analysis was conducted using BayesPhylogenies v1.0 [102] with a reversible jump mixture model [103] using a GTR model of sequence evolution with 4 discrete r rate categories. The analysis was run for 10,000,000 iterations with trees printed and saved every 10,000 iterations. Three independent rate matrices were assigned by BayesPhyloge-nies. Following this analysis the first 500,000 iterations were removed as burnin and the remaining 950 trees were used to create a consensus tree in BayesTrees v1.3 (

Sexual system

Male notostracans are readily identified by the lack of ovisacs, subtle morphological differences in carapace shape, numbers of legless rings and active mating behaviour in live individuals [26,32,48,56,104,105]. Females and hermaphrodites, although identical in external morphology and behaviour, differ histologically by the presence in hermaphrodites of an ovotestis (testicular lobes amongst the ovarian lobes) and by their ability to reproduce in isolation [40,50].

We compiled data from the literature for sex ratio, histology (i.e. presence/absence of ovotestis) and the inability/ability of females/hermaphrodites to reproduce in isolation (Table 1). Studies showing inability to reproduce in isolation were only included if reproduction in the presence of males was confirmed, to rule out poor lab rearing conditions or lack of reproductive maturity of individuals. In addition, we estimated sex ratio from available samples for a few populations (Additional file 1: Table S6). Using these data, we assigned populations as either being andro-dioecious or gonochoric. Androdioecious populations consist of hermaphrodites and males and exhibit skewed sex ratios with hermaphrodites found in greater numbers than males [12,17,48,82]. Gonochoric populations consist of males and females and have an approximately equal sex ratio. We did not categorise any population as purely hermaphroditic because this would necessitate showing a complete absence of males. Given that males in androdioecious species can be maintained by metapopulation dynamics [12,13,106] and can be present in exceedingly low frequencies (e.g. eight males per thousand in L. apus and similar proportions in T. cancriformis [52,53,107]) large samples sizes where no males are found would be needed to establish that a population is hermaphroditic [25]. In view of the sample sizes available to us we decided to conservatively categorise STUs into two sexual systems, gono-choric and androdioecious. In the AD notostracan species T. newberryi male proportions never exceed 27% [32,48] and in populations of AD Eulimnadia male proportions were always significantly lower than 50% male with a mode of ~20% [108]. Weeks et al. [108] did, however, note that upper values for population sex ratio of AD taxa overlapped with the lower values of gonochoric taxa in the 35% - 45% range. We therefore used a conservative population sex ratio cut-off of 30% male to assign an AD sexual system in the absence of additional histological or reproduction in isolation data in order to prevent misclassification due to stochastic variation in natural population sex ratios [51,109]. Populations with a male proportion of 30% - 45% were coded as equivocal and populations with male proportions greater than 45% were coded as gonochoric.

Character mapping

Sexual system was mapped onto the best scoring ML tree as a discrete character. STUs for which sexual system could not be inferred, or for which data was lacking, were left uncoded for sexual system in our analyses. The Leptestheria outgroup used to root the tree for character mapping analyses was also left uncoded for sexual system. MP reconstruction of ancestral states was conducted using Mesquite v2.74 [110] with an unordered model. In addition, we used Bayes-Multistate [111] implemented in BayesTraits v1.0 in an ML framework to evaluate four alternative models of sexual system evolution using likelihood ratio tests (D) assuming the result approximates a chi-squared distribution with degrees of freedom equal to the difference in the number of estimated parameters between the models. The simplest model is a one-parameter model, with a single rate of transition between gonochorism and AD and vice versa. The second model is a two-parameter model, where the transition rates from AD to gonochorism and vice versa can vary. The third and fourth models allow only unidirectional changes in sexual system, one from gonochorism to AD only, as in clam shrimps [25], and the other from AD to gono-chorism only. Ancestral character states were reconstructed based on the best fit model using the AddNode function of BayesTraits.

Testing the reproductive assurance hypothesis

We used a proxy for the exposure of STUs to glacial cycles, and therefore presumed range expansions, to test whether the reproductive assurance hypothesis is responsible for sexual system evolution in Notostraca [39,40]. As STUs found at higher latitudes are more likely to have recently re-colonised following the last glacial maxima than lower latitude ones, we expect AD STUs to be found at higher latitudes than gonochoric ones. The absolute latitude values at which gonochoric and AD STUs are found were compiled using the collection location of each representative STU as an unbiased representation of the latitude at which that lineage is found (Additional file 1: Table S5). We used the program BayesTraits [111] in an ML framework to conduct a t-test which accounts for the shared ancestry as implied by our best scoring ML phylogeny (phylogenetic t-test) to determine if latitude significantly differs between gonochoric STUs and ones where AD populations are found (the presence and absence of AD was incorporated using standard contrast or 'dummy' coding). We simultaneously estimated the parameter A which detects the phylogenetic signal in the data [112,113], if A is close to 1 there is strong phylogenetic signal if it is 0 there is no phylogenetic signal and the model collapses to an ordinary t-test.

Availability of Supporting Data

Supermatrix alignment used for phylogenetic inference has been deposited in the Dryad repository under doi:10.5061/dryad.480cf.

Additional file

Additional file 1: Mathers et al. High lability of sexualsystem over 250 million years of evolution in morphologically conservative tadpole shrimps.

Competing interests

The authors declare that they have no competing interest. Authors' contributions

TCM collected and analysed data and drafted the manuscript. AG, BH, RAJ and RH conceived the study. TZ collected some samples and provided some sequences. Allauthors read, contributed to and approved the final manuscript.


This work was part of T.C.M.'s Ph.D. project 'The genetics and evolutionary dynamics of reproductive mode in tadpole shrimps' funded by a NERC CASE Studentship (NE/G012318/1). We are gratefulto Diego Fontaneto, Marco Seminara, Dani Boix, Anna Badosa, Andy Green and Don Dasis for sending us samples. We would also like to thank Chris Venditti, Dave Lunt, Steve Moss and Elze Hesse for their help and advice. A.G. was funded by an Advanced NERC fellowship (NE/B501298/1).

Author details

1Schoolof Biological, Biomedicaland EnvironmentalSciences, University of Hull, HullHU6 7RX, UK. 2Department of Biology, University of Leicester, University Road, Leicester LE1 7RH, UK. 3Department of Life Sciences, The NaturalHistory Museum, CromwellRoad, London SW7 5BD, United Kingdom. 4Museum für Naturkunde Chemnitz, Moritzstrasse 20, D-09111, Chemnitz, Germany.

Received: 20 March 2012 Accepted: 25 January 2013 Published: 5 February 2013


1. Barrett SCH: The evolution of mating strategies in flowering plants. Trends Plant Sci 1998, 3:335-341.

2. Holsinger KE: Reproductive systems and evolution in vascular plants. Proc Natl Acad Sci USA 2000, 97:7037-7042.

3. Charlesworth D, Wright SI: Breeding systems and genome evolution. Curr Opin Genet Dev 2001, 11:685-690.

4. Barrett LG, Thrall PH, Burdon JJ, Linde CC: Life history determines genetic structure and evolutionary potential of host-parasite interactions. Trends Ecol Evol 2008, 23:678-685.

5. Charnov EL: The theory of sex allocation. Princeton, NJ: Princeton University Press; 1982.

6. Charnov EL, Bull JJ, Smith JM: Why be an hermaphrodite? Nature 1976, 263:125-126.

7. Darwin CR: The effects of cross and self fertilisation in the vegetable Kingdom. London: John Murray; 1876.

8. Baker HG: Self-compatibility and establishment after 'long-distance' dispersal. Evolution 1955, 9:347-349.

9. Baker HG: Support for Baker's Law-as a rule. Evolution 1967, 21:853-856.

10. Hesse E, Pannell JR: Density dependent pollen limitation and reproductive assurance in a wind pollinated herb with contrasting sexual systems. J Ecol 2011, 99:1531 -1539.

11. Glemin S, Bazin E, Charlesworth D: Impact of mating systems on patterns of sequence polymorphism in flowering plants. Proc Royal Society B-Biol Sci 2006, 273:3011 -3019.

12. Pannell JR: The evolution and maintenance of androdioecy. Annu Rev Ecol Syst 2002, 33:397-425.

13. Pannell J: The maintenanee of gynodioeey and androdioeey in a metapopulation. Evolution 1997, 51:10-20.

14. Charlesworth D: Androdioeey and the evolution of dioeey. Biol J Linn Soc 1984, 22:333-348.

15. Darwin C: The different forms of flowers on plants of the same species. London: J Murray; 1877.

16. Pannell JR: What is funetional androdioeey? Funct Ecol 2002,1б:862-865.

17. Weeks SC, Benvenuto C, Reed SK: When males and hermaphrodites eoexist: a review of androdioeey in animals. Integr Comp Biol 2006,4б:449-464.

18. Yusa Y, Yoshikawa M, Kitaura J, Kawane M, Ozaki Y, Yamato S, H0eg JT: Adaptive evolution of sexual systems in peduneulate barnaeles. Proc Royal Soc B: Biol Sci 2011, 279:959-966.

19. Urano S, Yamaguchi S, Yamato S, Takahashi S, Yusa Y: Evolution of dwarf males and a variety of sexual modes in barnaeles: an ESS approaeh. Evol Ecol Res 2009, 11:713-729.

20. Obbard DJ, Harris SA, Buggs RJA, Pannell JR: Hybridization, polyploidy, and the evolution of sexual systems in Mercurialis (Euphorbiaeeae). Evolution 2006, б0:1801 -1815.

21. Pannell JR, Obbard DJ, Buggs RJA: Polyploidy and the sexual system: what ean we learn from Mercurialis annua? Biol J Linn Soc 2004, 82:547-560.

22. Listen A, Rieseberg LH, Elias TS: Funetional androdioeey in the flowering plant Datisca glomerata. Nature 1990, 343:641-642.

23. Wolf DE, Takebayashi N: Pollen limitation and the evolution of androdioeey from dioeey. Am Nat 2004, 1б3:122-137.

24. Lloyd DG: Maintenanee of gynodioeey and androdioeey in angiosperms. Genetica 1975, 45:325-339.

25. Weeks SC, Sanderson TF, Reed SK Zofkova M, Knott B, Balaraman U, Pereira G, Senyo DM, Hoeh WR: Aneient androdioeey in the freshwater erustaeean Eulimnadia. Proc Royal Soc B-Biol Sci 2006, 273:725-734.

26. Longhurst AR: A review of the Notostraea. Bull Br Mus Nat Hist 1955,3:1-57.

27. Gore PJW: Triassie notostraeans in the Newark Supergroup, Culpeper Basin. Northern Virginia. J Paleontol 1986, б0:1086-1096.

28. Gall J-C, Grauvogel-Stamm L: The early middle Triassie 'grès à Voltzia' formation of eastern Franee: a model of environmental refugium. Paléontologie Générale 2005, 4:569-560. 584.

29. Barnard KH: A revision of the South Afriean Branehiopoda Phyllopoda. Annals South African Museum 1929, 29:181 -272.

30. Fryer G: Studies on the funetional morphology and biology of the Notostraea (Crustaeea, Branehiopoda). Philosophical Transact Royal Soc London Series B-Biol Sci 1988, 321:27-142.

31. Zierold T, Hanfling B, Gomez A: Reeent evolution of alternative reproduetive modes in the 'living fossil' Triops cancriformis. BMC Evol Biol 2007, 7:161.

32. Macdonald KS, Sallenave R, Cowley DE: Morphologieal and genetie variation in Triops (Branehiopoda: Notostraea) from ephemeral waters of the northern Chihuahuan Desert of North Ameriea. J Crustac Biol 2011, 31:468-484.

33. Mantovani B, Cesari M, Scanabissi F: Moleeular taxonomy and phylogeny of the 'living fossil' lineages Triops and Lepidurus (Branehiopoda: Notostraea). Zoologica Scripta 2004, 33:367-374.

34. Murugan G, Maeda-MartÍnez AM, Obregón-Barboza H, Hernández-Saavedra NY: Moleeular eharaeterization of the tadpole shrimp Triops (Branehiopoda: Notostraea) from the Baja California Peninsula, Méxieo: New insights on speeies diversity and phylogeny of the genus. Hydrobiologia 2002, 48б:101 -113.

35. King JL, Hanner R: Cryptie speeies in a "living fossil" lineage: taxonomie and phylogenetie relationships within the genus Lepidurus (Crustaeea: Notostraea) in North Ameriea. Mol Phylogenetics Evol 1998, 10:23-36.

36. Korn M, Green AJ, Machado M, GarcÍa-de-Lomas J, Cristo M, da Fonseca LC, Frisch D, Pérez-Bote JL, Hundsdoerfer AK: Phylogeny, moleeular eeology and taxonomy of southern Iberian lineages of Triops mauritanicus (Crustaeea: Notostraea). Org Divers Evol 2010,10:409-440.

37. Korn M, Hundsdoerfer AK: Evidenee for eryptie speeies in the tadpole shrimp Triops granarius (lueas, 18б4) (Crustaeea: Notostraea). Zootaxa 2006, 1257:57-68.

38. Korn M, Marrone F, Pérez-Bote JL, Machado M, Cristo M, da Fonseca LC, Hundsdoerfer AK: Sister speeies within the Triops cancriformis lineage (Crustaeea, Notostraea). Zoologica Scripta 2006, 35:301-322.

39. Longhurst AR: Evolution in the Notostraea. Evolution 1955, 9:84-86.

40. Bernard H: Hermaphroditism of the Apodidae. Nature 1891, 43:343-344.

41. Rogers DC: Revision of the Nearetie Lepidurus (Notostraea). J Crustac Biol 2001, 21:991-1006.

42. Hessen DO, Rueness EK, Stabell M: Circumpolar analysis of morphologieal and genetie diversity in the notostraean Lepidurus arcticus. Hydrobiologia 2004, 519:73-84.

43. Zierold T, Montero-Pau J, Hänfling B, Gómez A: Sex ratio, reproduetive mode and genetie diversity in Triops cancriformis. Freshw Biol 2009, 54:1392-1405.

44. Longhurst AR: Reproduetion in Notostraea (Crustaeea). Nature 1954, 173:781-782.

45. Murugan G, Obregón-Barboza H, Maeda-MartÍnez AM, Timms BV: Co-oeeurrenee of two tadpole shrimp, Triops cf. australiensis (Branehiopoda: Notostraea), lineages in middle Paroo, north-western New South Wales, with the first reeord of Triops hermaphrodites for the Australian eontinent. Australian J Zool 2009, 57:77-84.

46. Mitsumoto H, Yahata K: Evidenee of eross-fertilization in a gonoehorie population of the tadpole shrimp Triops numidicus (Crustaeea: Branehiopoda: Notostraea). Zool Sci 2006, 23:1109-1113.

47. Suno-Uchi N, Sasaki F, Chiba S, Kawata M: Morphologieal stasis and phylogenetie relationships in tadpole shrimps, Triops (Crustaeea: Notostraea). Biol J Linn Soc 1997, б1:439-457.

48. Sassaman C, Simovich MA, Fugate M: Reproduetive isolation and genetie differentiation in North Ameriean speeies of Triops (Crustaeea: Branehiopoda: Notostraea). Hydrobiologia 1997, 359:125-147.

49. GarcÍa-Velazco H, Obregón-Barboza H, RodrÍguez-Jaramillo C, Maeda-MartÍnez AM: Reproduetion of the tadpole shrimp Triops (Notostraea) in Mexiean waters. Curr Sci 2009, 9б:91 -97.

50. Sassaman C: Sex-ratio variation in female-biased populations of notostraeans. Hydrobiologia 1991, 212:169-179.

51. Machado M, Cristo M, Reis J, Cancela L: Biologieal data on Triops cancriformis mauritanicus (Ghigi, 1921) and Cyzicus grubei (Simon, 188б)—Crustaeea, Branehiopoda—in SW Portugal temporary ponds. Limnetica 1999, 1б:1-7.

52. Engelmann M, Hahn T: Vorkommen von Lepidurus apus, Triops cancriformis, Eubranchipus (Siphonophanes) grubei, Tanymastix stagnalis und Branchipus schaefferi in Deutsehland und Österreieh (Crustaeea: Notostraea und Anostraea). Faunistische Abhandlungen 2005, 25:3-67.

53. Simon E: Étude sur les erustaeés du sous-ordre des phyllopodes. Annales de la Société entomologique de France 1886, б:393-460.

54. Bernard HM: XL.—Hermaphroditism among the Apodidœ. J Nat History Series 6 1896,17:296-309.

55. Stephan S: Mannehen von Lepidurus apus (Notostraea, erustaeea) in der Ruhstadter Elbtalaue. Abhandlungen und Berichte fur Naturkunde 2008, 31:47-57.

56. Linder F: Contributions to the morphology and the taxonomy of the Branehiopoda Notostraea, with speeial referenee to the North Ameriean speeies. Proc United States Nat Museum 1952, 102:1 -69.

57. Wojtasik B, Brytka-Wotk M: Reproduetion and genetie strueture of a freshwater erustaeean Lepidurus arcticus from Spitsbergen. Polish Polar Res 2010, 31:33-44.

58. Bushnell JH, Byron ER: Morphologieal variability and distribution of aquatie invertebrates (prineipally Crustaeea) from the Cumberland Peninsula and Frobisher Bay regions, Baffin Island, NWT, Canada. Arct Alp Res 1979, 1 1:159-177.

59. Margraf J, Maass B: Zur Ökologie der temporären Süßwasserflaehseen des Tafelberges "Giara di Gesturi" auf Sardinien. Spixiana 1982, 5:69-99.

60. Packard A: New phyllopod erustaeeans. Am Nat 1875, 9:311 -312.

61. Lynch JE: Lepidurus lemmoni Holmes: a redeseription with notes on variation and distribution. Trans Am MicroscSoc 1966, 85:181-192.

62. Scanabissi F, Mondini C: A survey of the reproduetive biology in Italian branehiopods. Hydrobiologia 2002, 48б:263-272.

63. Kuller Z, Gasith A: Comparison of the hatehing proeess of the tadpole shrimps Triops cancriformis and Lepidurus apus lubbocki (Notostraea) and its relation to their distribution in rain-pools in Israel. Hydrobiologia 1996, 335:147-157.

64. Ahl J: Faetors affeeting eontributions of the tadpole shrimp, Lepidurus packardi, to its oversummering egg reserves. Hydrobiologia 1991,


65. Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogeneties with BEAUti and the BEAST 1.7. Mol Biol Evol 2012, 29:1969-1973.

66. Pagel M: Inferring the historical patterns of biological evolution. Nature 1999, 401:877-884.

67. Guthörl P: Die arthropoden aus dem Carbon und Perm des Saar-NahePfalz-gebietes. Preussischen Geologischen Landesanstalt 1934, 164:1-219.

68. Gand G, Garric J, Lapeyrie J: Biocénoses à triopsidés (Crustacea, Branchiopoda) du Permien du Bassin de Lodève (France). Geobios 1997, 30:673-700.

69. Weeks SC, Chapman EG, Rogers DC, Senyo DM, Hoeh WR: Evolutionary transitions among dioecy, androdioecy and hermaphroditism in limnadiid clam shrimp (Branchiopoda: Spinicaudata). J Evol Biol 2009, 22:1781-1799.

70. Zucker N, Cunningham M, Adams HP: Anatomical evidence for androdioecy in the clam shrimp Eulimnadia texana. Hydrobiologia 1997, 359:171-175.

71. Weeks SC, Posgai RT, Cesari M, Scanabissi F: Androdioecy inferred in the clam shrimp Eulimnadia agassizii (Spinicaudata: Limnadiidae). J Crustac Biol 2005, 25:323-328.

72. Marcus V, Weeks S: The effects of pond duration on the life history traits of an ephemeral pond crustacean Eulimnadia texana. Hydrobiologia 1997, 359:213-221.

73. MacKay WP, Loring SJ, Frost TM, Whitford WG: Population dynamics of a playa community in the Chihuahuan Desert. Southwest Nat 1990, 35:393-402.

74. Simovich MA: Crustacean biodiversity and endemism in California's ephemeral wetlands. In Ecology, conservation, and management of vernal pool ecosystems - proceedings from a 1996 conference. Edited by Witham CW, Bauder ET, Belk D, Ferren WR, Ornduff R. Sacramento: California Native Plant Society; 1998:107-118.

75. Takahashi F: Pioneer life of the tadpole shrimps, Triops spp. (Notostraca: Triopsidae). Applied Entomology Zool 1977, 12:104-117.

76. Hamer M, Appleton C: Life history adaptations of phyllopods in response to predators, vegetation, and habitat duration in north-eastern Natal. Hydrobiologia 1991, 212:105-116.

77. Weeks SC: Life-history variation under varying degrees of intraspecific competition in the tadpole shrimp Triops longicaudatus (LeConte).

J Crustac Biol 1990, 10:498-503.

78. Weeks SC, Benvenuto C, Sanderson TF, Duff RJ: Sex chromosome evolution in the clam shrimp, Eulimnadia texana. J Evol Biol 2010, 23:1100-1106.

79. Sassaman C, Weeks SC: The genetic mechanism of sex determination in the chonchostracan shrimp Eulimnadia texana. Am Nat 1993,141:314-328.

80. Sassaman C: Sex-ratio and sex determination in western populations of Triops (Notostraca). Am Zool 1989, 29:A66-A66.

81. Hoeg J: Sex and the single cirripede: a phylogenetic perspective. Crustacean issues 1995, 10:195-207.

82. Pannell JR: Consequences of inbreeding depression due to sex-linked loci for the maintenance of males and outcrossing in branchiopod crustaceans. Genet Res 2008, 90:73-84.

83. Otto SP, Sassaman C, Feldman MW: Evolution of sex determination in the conchostracan shrimp Eulimnadia texana. Am Nat 1993, 141:329-337.

84. Pannell JR, Dorken ME: Colonisation as a common denominator in plant metapopulations and range expansions: effects on genetic diversity and sexual systems. Landsc Ecol 2006, 21:837-848.

85. Obbard DJ, Harris SA, Pannell JR: Sexual systems and population genetic structure in an annual plant: testing the metapopulation model. Am Nat 2006, 167:354-366.

86. Pannell JR, Barrett SCH: Baker's law revisited: reproductive assurance in a metapopulation. Evolution 1998, 52:657-668.

87. Igic B, Lande R, Kohn JR: Loss of self-incompatibility and its evolutionary consequences. Int J Plant Sci 2008, 169:93-104.

88. Blaustein L, Schwartz SS: Why study ecology in temporary pools? Israel J Zool 2001, 47:303-312.

89. Mantovani B, Cesari M, Scanabissi F: Molecular taxonomy and phylogeny of Italian Lepidurus taxa (Branchiopoda: Notostraca). Italian J Zool 2009, 76:358-365.

90. Montero-Pau J, Gómez A, Muñoz J: Application of an inexpensive and high-throughput genomic DNA extraction method for the molecular ecology of zooplanktonic diapausing eggs. Limnol Oceanography-Methods 2008, 6:218-222.

91. 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.

92. Fredslund J, Schauser L, Madsen LH, Sandal N, Stougaard J: PriFi: using a multiple alignment of related sequences to find primers for amplification of homologs. Nucleic Acids Res 2005, 33:W516-W520.

93. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 2011, 28:2731-2739.

94. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004, 32:1792-1797.

95. Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler AP: Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol 2006, 55:595-609.

96. Drummond A, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 2007, 7:214.

97. Stocsits RR, Letsch H, Hertel J, Misof B, Stadler PF: Accurate and efficient reconstruction of deep phylogenies from structured RNAs. Nucleic Acids Res 2009, 37:6184-6193.

98. Letsch HO, Kjer KM: Potential pitfalls of modelling ribosomal RNA data in phylogenetic tree reconstruction: evidence from case studies in the Metazoa. BMC Evol Biol 2011, 11:146.

99. Leigh JW, Susko E, Baumgartner M, Roger AJ: Testing congruence in phylogenomic analysis. Syst Biol 2008, 57:104-115.

100. Stamatakis A: RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006, 22:2688-2690.

101. Kuck P, Meusemann K: FASconCAT: convenient handling of data matrices. Mol Phylogenetics Evol 2010, 56:1115-1118.

102. Pagel M, Meade A: A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data. Syst Biol 2004, 53:571-581.

103. Pagel M, Meade A: Mixture models in phylogenetic inference. In

Mathematics of evolution and phylogeny. Edited by Gascuel O. New York: Oxford University Press; 2005.

104. Obregon-Barboza H, Maeda-Martinez AM, Murugan G, Timms BV, Grygier MJ, Rogers DC, Rodriguez-Almaraz G, Dumont HJ: Morphology and systematic significance of the mystax, a hitherto undescribed cephalic structure of males in certain Notostraca (Branchiopoda). J Crustac Biol 2007, 27:18-23.

105. Akita M: On the reproduction of Triops longicaudatus (LECONTE). Doubutsugaku zasshi 1971, 80:242-250.

106. Pannell JR: A hypothesis for the evolution of androdioecy: the joint influence of reproductive assurance and local mate competition in a metapopulation. Evol Ecol 2000, 14:195-211.

107. Heidecke D, Neumann V: Zur Verbreitung und ökologie von Triops cancriformis BOSC und Lepidurus apus L. In der DDR. Hercynia NF 1987, 24:166-173.

108. Weeks SC, Sanderson TF, Zofkova M, Knott B: Breeding systems in the clam shrimp family Limnadiidae (Branchiopoda, Spinicaudata). Invertebr Biol 2008, 127:336-349.

109. Petrov B, Cvetkovic DM: Community structure of branchiopods (Anostraca, Notostraca and Conchostraca) in the Banat province in Yugoslavia. Hydrobiologia 1997, 359:23-28.

110. Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis v2.75.

111. Pagel M, Meade A, Barker D: Bayesian estimation of ancestral character states on phylogenies. Syst Biol 2004, 53:673-684.

112. Freckleton R, Harvey P, Pagel M: Phylogenetic analysis and comparative data: a test and review of evidence. Am Nat 2002, 160:712-726.

113. Pagel M: The maximum likelihood approach to reconstructing ancestral character states of discrete characters on phylogenies. Syst Biol 1999, 48:612-622.


Cite this article as: Mathers et al.: High lability of sexual system over 250 million years of evolution in morphologically conservative tadpole shrimps. BMC Evolutionary Biology 2013 13:30.