Mitogenomics of Cladocera (Branchiopoda): Marked gene order rearrangements and independent predation roots

Cladocera (Crustacea: Branchiopoda) is a key group of invertebrates. Despite a long history of phylogenetic research, relationships within this group remain disputed. We here provide new insights based on 15 new mitochondrial genomes obtained from high-throughput sequencing (HTS) and 40 mitogenomes extracted from published HTS datasets. Together with 25 mitogenomes from GenBank, we generated a matrix of 80 mitogenomes, 44 of them belonging to Cladocera. We also obtained a matrix with 168 nuclear orthologous genes to further assess the phylogenetic result from mitogenomes based on published data and one new HTS data ofLeptodora. Maximum likelihood and Bayesian phylogenetic analyses recovered all Branchiopoda orders as monophyletic and supported a sister-group relationship between Anomopoda and Onychopoda, making the taxon Gymnomera paraphyletic and supporting an independent origin of predatory Haplopoda and Onychopoda. The nuclear phylogeny and topological tests also support Gymnomera as paraphyletic, and the nuclear phylogeny strongly supports a sister-group relationship between Ctenopoda and Haplopoda. We provide a fossil-calibrated time tree, congruent with a Carboniferous origin for Cladocera and a subsequent diversification of the crown group of Anomopoda, Onychopoda, and Ctenopoda, at least in the Triassic. Despite their long evolutionary history, non-Cladoceran Branchiopoda exhibited high mitogenome structural stability. On the other hand, 21 out of 24 gene rearrangements occurred within the relatively younger Cladocera. We found the differential base compositional skewness patterns between Daphnia s.s. and Ctenodaphnia, which might be related to the divergence between these taxa. We also provide evidence to support the recent finding that Spinicaudata possesses mitogenomes with inversed compositional skewness without gene rearrangement. Such a pattern has only been reported in Spinicaudata.

Three nuclear protein-coding genes corroborate a recent phylogenomic model of the Branchiopoda (Crustacea) and provide estimates of the divergence times of the major branchiopodan taxa

The class Branchiopoda (Crustacea) shows great diversity in morphology and lifestyle among its constituent higher-level taxa: Anostraca, Notostraca, Laevicaudata, Spinicaudata, Cyclestherida and Cladocera. The phylogenetic relationships among these taxa have long been controversial. We sequenced three orthologous nuclear genes that encode the catalytic subunit of DNA polymerase delta and the largest and second-largest subunits of RNA polymerase II in the expectation that the amino acid sequences encoded by these genes might be effective in clarifying branchiopod phylogeny and estimating the times of divergence of the major branchiopodan taxa. The results of phylogenetic analyses based on these amino acid sequences support the monophyly of Branchiopoda and provide strong molecular evidence in support of the following phylogenetic relationships: (Anostraca, (Notostraca, (Laevicaudata, (Spinicaudata, (Cyclestherida, Cladocera))))). Within Cladocera, comparison of the nucleotide sequences of these same genes shows Ctenopoda to be the sister group of Haplopoda + Anomopoda. Three statistical tests based on the present amino acid sequence data-the approximately unbiased test, Kishino-Hasegawa test and weighted Shimodaira-Hasegawa test-tend to refute most of the previous molecular phylogenetic studies on Branchiopoda, which have placed Notostraca differently than here; however, our results corroborate those of one recent phylogenomic study, thus confirming the effectiveness of these three genes to investigate relationships among branchiopod higher taxa. Divergence time estimates calibrated on the basis of fossil evidence suggest that the first divergence of extant branchiopods occurred about 534 Ma during the early Cambrian period and that diversification within the extant branchiopod lineages started in or after the late Permian.

Tetraconatan phylogeny with special focus on Malacostraca and Branchiopoda: highlighting the strength of taxon-specific matrices in phylogenomics

Understanding the evolution of Tetraconata or Pancrustacea-the clade that includes crustaceans and insects-requires a well-resolved hypothesis regarding the relationships within and among its constituent taxa. Here, we assembled a taxon-rich phylogenomic dataset focusing on crustacean lineages based solely on genomes and new-generation Illumina-generated transcriptomes, including 89 representatives of Tetraconata. This constitutes, to our knowledge, the first phylogenomic study specifically addressing internal relationships of Malacostraca (with 26 species included) and Branchiopoda (36 species). Seven matrices comprising 81-684 orthogroups and 17 690-242 530 amino acid positions were assembled and analysed under five different analytical approaches. To maximize gene occupancy and to improve resolution, taxon-specific matrices were designed for Malacostraca and Branchiopoda. Key tetraconatan taxa (i.e. Oligostraca, Multicrustacea, Branchiopoda, Malacostraca, Thecostraca, Copepoda and Hexapoda) were monophyletic and well supported. Within Branchiopoda, Phyllopoda, Diplostraca, Cladoceromorpha and Cladocera were monophyletic. Within Malacostraca, the clades Eumalacostraca, Decapoda and Reptantia were well supported. Recovery of Caridoida or Peracarida was highly dependent on the analysis for the complete matrix, but it was consistently monophyletic in the malacostracan-specific matrices. From such examples, we demonstrate that taxon-specific matrices and particular evolutionary models and analytical methods, namely CAT-GTR and Dayhoff recoding, outperform other approaches in resolving certain recalcitrant nodes in phylogenomic analyses.

The doublesex gene integrates multi-locus complementary sex determination signals in the Japanese ant, Vollenhovia emeryi

A female diploid, male haploid sex determination system (haplodiploidy) is found in hymenopteran taxa, such as ants, wasps, bees and sawflies. In this system, a single, complementary sex-determination (sl-CSD) locus functions as the primary sex-determination signal. In the taxa that has evolved this system, females and males are heterozygous and hemi/homozygous at the CSD locus, respectively. While the sl-CSD system enables females to alter sex ratios in the nest, it carries a high cost in terms of inbreeding, as individuals that are homozygous at the CSD locus become sterile diploid males. To counter this risk, some of hymenopteran species have evolved a multi-locus CSD (ml-CSD) system, which effectively reduces the proportion of sterile males. However, the mechanism by which these multiple primary signals are integrated and how they affect the terminal sex-differentiation signal of the molecular cascade have not yet been clarified. To resolve these questions, we examined the molecular cascade in the Japanese ant Vollenhovia emeryi, which we previously confirmed has two CSD loci. Here, we showed that the sex-determination gene, doublesex (dsx), which is highly conserved among phylogenetically distant taxa, is responsible for integrating two CSD signals in V. emeryi. After identifying and characterizing dsx, genotypes containing two CSD loci and splicing patterns of dsx were found to correspond to the sexual phenotype, suggesting that two primary signals are integrated into dsx. These findings will facilitate future molecular and functional studies of the sex determination cascade in V. emeryi, and shed light on the evolution and diversification of sex determination systems in insects.

Molecular cloning of doublesex genes of four cladocera (water flea) species

Background

The gene doublesex (dsx) is known as a key factor regulating genetic sex determination in many organisms. We previously identified two dsx genes (DapmaDsx1 and DapmaDsx2) from a freshwater branchiopod crustacean, Daphnia magna, which are expressed in males but not in females. D. magna produces males by parthenogenesis in response to environmental cues (environmental sex determination) and we showed that DapmaDsx1 expression during embryonic stages is responsible for the male trait development. The D. magna dsx genes are thought to have arisen by a cladoceran-specific duplication; therefore, to investigate evolutionary conservation of sex specific expression of dsx genes and to further assess their functions in the environmental sex determination, we searched for dsx homologs in four closely related cladoceran species.

Results

We identified homologs of both dsx genes from, D. pulex, D. galeata, and Ceriodaphnia dubia, yet only a single dsx gene was found from Moina macrocopa. The deduced amino acid sequences of all 9 dsx homologs contained the DM and oligomerization domains, which are characteristic for all arthropod DSX family members. Molecular phylogenetic analysis suggested that the dsx gene duplication likely occurred prior to the divergence of these cladoceran species, because that of the giant tiger prawn Penaeus monodon is rooted ancestrally to both DSX1 and DSX2 of cladocerans. Therefore, this result also suggested that M. macrocopa lost dsx2 gene secondarily. Furthermore, all dsx genes identified in this study showed male-biased expression levels, yet only half of the putative 5’ upstream regulatory elements are preserved in D. magna and D. pulex.

Conclusions

The all dsx genes of five cladoceran species examined had similar amino acid structure containing highly conserved DM and oligomerization domains, and exhibited sexually dimorphic expression patterns, suggesting that these genes may have similar functions for environmental sex determination in cladocerans.

The ecoresponsive genome of Daphnia pulex

We describe the draft genome of the microcrustacean Daphnia pulex, which is only 200 megabases and contains at least 30,907 genes. The high gene count is a consequence of an elevated rate of gene duplication resulting in tandem gene clusters. More than a third of Daphnia's genes have no detectable homologs in any other available proteome, and the most amplified gene families are specific to the Daphnia lineage. The coexpansion of gene families interacting within metabolic pathways suggests that the maintenance of duplicated genes is not random, and the analysis of gene expression under different environmental conditions reveals that numerous paralogs acquire divergent expression patterns soon after duplication. Daphnia-specific genes, including many additional loci within sequenced regions that are otherwise devoid of annotations, are the most responsive genes to ecological challenges.

Coexisting cyclic parthenogens comprise a holocene species flock in Eubosmina

BACKGROUND

Mixed breeding systems with extended clonal phases and weak sexual recruitment are widespread in nature but often thought to impede the formation of discrete evolutionary clusters. Thus, cyclic parthenogens, such as cladocerans and rotifers, could be predisposed to "species problems" and a lack of discrete species. However, species flocks have been proposed for one cladoceran group, Eubosmina, where putative species are sympatric, and there is a detailed paleolimnological record indicating a Holocene age. These factors make the Eubosmina system suitable for testing the hypotheses that extended clonal phases and weak sexual recruitment inhibit speciation. Although common garden experiments have revealed a genetic component to the morphotypic variation, the evolutionary significance of the morphotypes remains controversial.

METHODOLOGY/PRINCIPAL FINDINGS

In the present study, we tested the hypothesis of a single polymorphic species (i.e., mixing occurs but selection maintains genes for morphology) in four northern European lakes where the morphotypes coexist. Our evidence is based on nuclear DNA sequence, mitochondrial DNA sequence, and morphometric analysis of coexisting morphotypes. We found significant genetic differentiation, genealogical exclusivity, and morphometric differentiation for coexisting morphotypes.

CONCLUSIONS

We conclude that the studied morphotypes represent a group of young species undergoing speciation with apparent reproductive barriers despite coexistence in the freshwater pelagic zone.

Meiosis genes in Daphnia pulex and the role of parthenogenesis in genome evolution

Background

Thousands of parthenogenetic animal species have been described and cytogenetic manifestations of this reproductive mode are well known. However, little is understood about the molecular determinants of parthenogenesis. The Daphnia pulex genome must contain the molecular machinery for different reproductive modes: sexual (both male and female meiosis) and parthenogenetic (which is either cyclical or obligate). This feature makes D. pulex an ideal model to investigate the genetic basis of parthenogenesis and its consequences for gene and genome evolution. Here we describe the inventory of meiotic genes and their expression patterns during meiotic and parthenogenetic reproduction to help address whether parthenogenesis uses existing meiotic and mitotic machinery, or whether novel processes may be involved.

Results

We report an inventory of 130 homologs representing over 40 genes encoding proteins with diverse roles in meiotic processes in the genome of D. pulex. Many genes involved in cell cycle regulation and sister chromatid cohesion are characterized by expansions in copy number. In contrast, most genes involved in DNA replication and homologous recombination are present as single copies. Notably, RECQ2 (which suppresses homologous recombination) is present in multiple copies while DMC1 is the only gene in our inventory that is absent in the Daphnia genome. Expression patterns for 44 gene copies were similar during meiosis versus parthenogenesis, although several genes displayed marked differences in expression level in germline and somatic tissues.

Conclusion

We propose that expansions in meiotic gene families in D. pulex may be associated with parthenogenesis. Taking into account our findings, we provide a mechanistic model of parthenogenesis, highlighting steps that must differ from meiosis including sister chromatid cohesion and kinetochore attachment.

Nucleotide polymorphism and within-gene recombination in Daphnia magna and D. pulex, two cyclical parthenogens

Theory predicts that partially asexual organisms may make the "best of both worlds": for the most part, they avoid the costs of sexual reproduction, while still benefiting from an enhanced efficiency of selection compared to obligately asexual organisms. There is, however, little empirical data on partially asexual organisms to test this prediction. Here we examine patterns of nucleotide diversity at eight nuclear loci in continentwide samples of two species of cyclically parthenogenetic Daphnia to assess the effect of partial asexual reproduction on effective population size and amount of recombination. Both species have high nucleotide diversities and show abundant evidence for recombination, yielding large estimates of effective population sizes (300,000-600,000). This suggests that selection will act efficiently even on mutations with small selection coefficients. Divergence between the two species is less than one-tenth of previous estimates, which were derived using a mitochondrial molecular clock. As the two species investigated are among the most distantly related species of the genus, this suggests that the genus Daphnia may be considerably younger than previously thought. Daphnia has recently received increased attention because it is being developed as a model organism for ecological and evolutionary genomics. Our results confirm the attractiveness of Daphnia as a model organism, because the high nucleotide diversity and low linkage disequilibrium suggest that fine-scale mapping of genes affecting phenotypes through association studies should be feasible.

The nauplius eye complex in 'conchostracans'(Crustacea, Branchiopoda: Laevicaudata, Spinicaudata, Cyclestherida) and its phylogenetic implications

The nauplius eye in Cyclestherida, Laevicaudata and Spinicaudata (previously collectively termed Conchostraca) consists of four cups of inverse sensory cells separated by a pigment layer and a tapetum layer. There are two lateral and two medial cups, a ventral medial cup and a posterior medial cup. The pigment and tapetum layers contain two different kinds of pigment granules, the inner pigment layer relatively large, dark (and electron dense) granules, and the outer tapetum layer light, reflective pigment granules. The presence of four cups and two different kinds of pigment granules are interpreted as autapomorphies of Phyllopoda. The position and shape of the nauplius eye in Spinicaudata is very distinct and herein interpreted as an autapomorphy of this taxon. Additional frontal eyes might be present dorsally or ventrally in varying proximity to the nauplius eye, but they have separate nerves from their sensory cells to the nauplius eye centre in the protocerebrum. Rhabdomeric structures are present in all these frontal eyes, evidencing their light sensitivity. In Lynceus biformis and L. tatei (Laevicaudata), two pairs of frontal eyes were found. In Cyclestheria hislopi (Cyclestherida), an unpaired ventral frontal eye is present. We did not find additional frontal eyes in Limnadopsis parvispinus and Caenestheriella sp. (Spinicaudata).

Probing the relationships of the branchiopod crustaceans

The Branchiopoda display extraordinary variation in body form, even within the morphologically diverse crustaceans. To fully understand the origin and evolution of these morphological reconfigurations, a robust phylogeny of the group is essential. To infer the affinities among branchiopods, we employed two approaches to taxon and gene sampling, presented new sequence data from three genes, incorporated previously published sequence data from three additional genes, and utilized comprehensive techniques of phylogeny reconstruction. The results provided support for a number of longstanding hypotheses concerning the relationships among the orders. For example, we obtained support for the Cladoceromorpha and Gymnomera, and favoured a unique arrangement of the cladoceran orders. A few affinities remain to be resolved, particularly at the base of the Phyllopoda and within the Anomopoda. However, the results suggest that increased gene sampling is recommended for future investigations of branchiopod systematics.

Larval development ofLynceus brachyurus (Crustacea, Branchiopoda, Laevicaudata): Redescription of unusual crustacean nauplii, with special attention to the molt between last nauplius and first juvenile

The larval development of "conchostracans" has received only scattered attention. Here I present the results of a study on the larval (naupliar) development and the metamorphosis of Lynceus brachyurus, a member of the bivalved branchiopod order the Laevicaudata. Lynceus brachyurus is the only species of the "Conchostraca" in Denmark. The phylogenetic position of the Laevicaudata has traditionally been a source of controversy, and this study does not solve the question completely. This work focuses on features potentially important for phylogeny. The general appearance of the larvae of L. brachyurus has been known for more than a century and a half, and some of its unique features include a large, larval dorsal shield; a huge, plate-like labrum; and a pair of immovable, horn-like antennules. However, many details relating to limb morphology, potentially important for phylogeny, have not been studied previously. Based on size categories, five or six larval stages can be recognized. The larvae approximately double their length and width during development (length: 230-520 microm). Most morphological features stay largely unchanged during development, but the antennal coxal masticatory spines are significant exceptions: they become bifid after one of the first molts. In all larval stages only the antennae and the mandibles actively move. In late naupliar stages the trunk limbs become visible as rows of laterally placed, undeveloped, and still immovable lobes. Swimming is performed by the antennae, whereas the mandibles appear to be involved mainly in feeding, as in other branchiopod larvae. The last naupliar stage undergoes a small metamorphosis to the first juvenile stage, the details of which in part were studied by following the premolt juvenile condition through the cuticle of the last stage nauplius. Among other changes there is a characteristic change in the shape and morphology of the univalved dorsal naupliar shield to a bivalved juvenile carapace. The general morphologies of the antennae and the mandibles are very similar to those of other branchiopod larvae and fall well within the "branchiopod naupliar feeding apparatus" recognized as a branchiopod synapomorphy by Olesen (2003), but some specific features shared with the larvae of other "conchostracans" are also identified. These special "conchostracan" features include: 1) a similar antennular setation; 2) a similar comb-like setulation of the bifid antennal coxal processes; and 3) mandibular palpsetae with setules condensed. In light of recent suggestions concerning branchiopod phylogeny (Cyclestheria as a sister group to the Cladocera), these similarities probably do not support a monophyletic "Conchostraca" but rather are symplesiomorphies of this taxon. A final decision must await a phylogenetic analysis of a more complete set of characters.

On the ontogeny of Leptodora kindtii (Crustacea, Branchiopoda, Cladocera), with notes on the phylogeny of the Cladocera

Leptodora kindtii, a large predaceous cladoceran, is among the most deviant species of the Cladocera. Therefore, its phylogenetic position has traditionally proven difficult to determine. Its many peculiar features include, among others, long, stenopodous, forwardly directed trunk limbs, a posteriorly placed dorsal brood pouch, a tri-lobed lower lip, and a long, segmented abdomen. This study describes the ontogeny of L. kindtii (Haplopoda), including general body proportions, appendages, the carapace, and other external structures in an attempt to facilitate the comparison of its aberrant morphology to that of other branchiopods. In general, the early embryos are similar to the early embryos of other cladoceran taxa with respect to body shape and size and position and orientation of the early limb buds. Many of the unusual features of L. kindtii appear late in ontogeny. The carapace appears at an early stage as a pair of dorsolateral swellings in a position corresponding to the gap between the mandibles and the first pair of trunk limbs; it later becomes posteriorly transposed by a gradual fusion of its more anterior parts to the dorsal side of the thorax. The tri-lobed "lower lip," under the labrum of the late embryo and the adult, develops as a fusion of the first maxillae (lateral lobes) to an elevated sternal region behind the mouth (median lobe). The stenopodous, segmented trunk limbs in the adult develop from embryonic, elongate, subdivided limb buds, similar to those seen in early stages of other branchiopods. Two conflicting possibilities for the phylogeny of the Cladocera, involving two different positions of L. kindtii (Haplopoda), are discussed. Several characters support a sister-group relationship between the Haplopoda and Onychopoda. However, some characters support the Anomopoda and Onychopoda as sister groups, leaving the Haplopoda outside this clade. In contrast to recent suggestions, we prefer to retain the term "Cladocera" in its original sense as comprising the Haplopoda, Ctenopoda, Anomopoda, and Onychopoda.

Structural rRNA characters support monophyly of raptorial limbs and paraphyly of limb specialization in water fleas

The evolutionary success of arthropods has been attributed partly to the diversity of their limb morphologies. Large morphological diversity and increased specialization are observed in water flea (Cladocera) limbs, but it is unclear whether the increased limb specialization in different cladoceran orders is the result of shared ancestry or parallel evolution. We inferred a robust among-order cladoceran phylogeny using small-subunit and large-subunit rRNA nuclear gene sequences, signature sequence regions, novel stem-loops and secondary structure morphometrics to assess the phylogenetic distribution of limb specialization. The sequence-based and structural rRNA morphometric phylogenies were congruent and suggested monophyly of orders with raptorial limbs, but paraphyly of orders with reduced numbers of specialized limbs. These results highlight the utility of complex molecular structural characters in resolving ancient rapid radiations.

Accelerated molecular evolution in halophilic crustaceans

In contrast to the stable ionic composition of the oceans, inland waters show striking diversity, possessing salt concentrations varying from I mM to 5 M. Although species diversity is highest in fresh water, some lineages have colonized hypersaline environments where they encounter elevated levels of both ultraviolet (UV) radiation and osmotic stress. This study compares rates of evolution in halophilic and freshwater taxa for two groups of microcrustaceans, anostracans and daphniids, from Australia and North America. The results establish that halophilic species show consistent rate acceleration, involving elevated levels of both insertion/deletion events and of nucleotide substitutions. The elevated pace of molecular evolution does not appear to be linked to selection or to other agents that are known to influence the supply rate of mutations, such as UV exposure, generation length, or shifts in metabolic rate. However, variance in ionic strength, which is known to have potent effects on DNA-protein interactions as well as on the structural properties of DNA and proteins, might account for the lowered fidelity of DNA replication in life from hypersaline settings. Regardless of its cause, the consistent rate acceleration in halophiles suggests that past efforts to employ sequence divergences to date events, such as the age of asexual lineages in Artemia, have resulted in serious overestimates. More generally, the results indicate that coordinated shifts in rates of molecular evolution may occur in lineages exposed to extreme environmental conditions.

Rate acceleration and long-branch attraction in a conserved gene of cryptic daphniid (Crustacea) species

The nuclear large subunit (LSU) rRNA gene is a rich source of phylogenetic characters because of its large size, mosaic of slowly and rapidly evolving regions, and complex secondary structure variation. Nevertheless, many studies have indicated that inconsistency, bias, and gene-specific error (e.g., within-individual gene family variation, cryptic sequence simplicity, and sequence coevolution) can complicate animal phylogenies based on LSU rDNA sequences. However, most of these studies sampled small gene fragments from expansion segments--among animals only five nonchordate complete LSU sequences are published. In this study, we sequenced near-complete nuclear LSU genes from 11 representative daphniids (Crustacea). The daphniid expansion segment V6 was larger and showed more length variation (90-351 bp) than is found in all other reported LSU V6 sequences. Daphniid LSU (without the V6 region) phylogenies generally agreed with the existing phylogenies based on morphology and mtDNA sequences. Nevertheless, a major disagreement between the LSU and the expected trees involved a positively misleading association between the two taxa with the longest branches, Daphnia laevis and D. occidentalis. Both maximum parsimony (MP) and maximum likelihood (ML) optimality criteria recovered this association, but parametric simulations indicated that MP was markedly more sensitive to this bias than ML. Examination of data partitions indicated that the inconsistency was caused by increased nucleotide substitution rates in the branches leading to D. laevis and D. occidentalis rather than among-taxon differences in base composition or distribution of sites that are free to vary. These results suggest that lineage-specific rate acceleration can lead to long-branch attraction even in the conserved genes of animal species that are almost morphologically indistinguishable.

Macrogeographic patterns of breeding system diversity in the Daphnia pulex group from the United States and Mexico

Studies on the biogeographical patterning of reproductive systems promise to extend understanding of the factors which modulate breeding system transitions. Two closely allied cladoceran crustaceans, Daphnia pulex and D. pulicaria, show varied modes of reproduction, with populations reproducing by either cyclic or obligate parthenogenesis. Prior studies have provided a detailed understanding of their breeding system diversity in the polar and cold temperate regions of North America. The present investigation extends this analysis, characterizing breeding systems and clonal diversity at sites throughout the United States and Mexico. Genotypic diversity in these southern areas was high, but only diploids were detected, indicating that polyploids are restricted to the north. F(1) hybrids and their two parental species were present in most areas, although their frequencies varied geographically. Hybrids invariably reproduced by obligate asexuality, but both parental taxa showed regional shifts in their breeding system. The complexity of these latter patterns suggests that they reflect the interplay of historical factors and selection.

The phylogenetic relationships of "predatory water-fleas" (Cladocera: Onychopoda, Haplopoda) inferred from 12S rDNA

Within the Cladocera, the water-fleas, four major taxa can be distinguished: Anomopoda, Ctenopoda, Haplopoda, and Onychopoda. Haplopoda and Onychopoda are called "predatory water-fleas." The Haplopoda is monotypic; its only representative, Leptodora kindtii, is common in palearctic and nearctic freshwater bodies. The Onychopoda show a remarkable geographic distribution. Most of the described species are restricted to the Caspian Sea, the Aral Sea, and peripheral areas of the Black Sea, including the Sea of Azov--all remnants of the Eastern Paratethys. The remaining onychopods are either freshwater inhabitants or marine animals, widespread in the world oceans. We present molecular evidence for a sister group relationship between Haplopoda and Onychopoda within the Cladocera. The Onychopoda and its three families are monophyletic. We suggest an independent invasion into the Ponto-Caspian basin at least three times, twice originating in the palearctic freshwater bodies and once starting from the world oceans.

Affinities among anostracan (Crustacea: Branchiopoda) families inferred from phylogenetic analyses of multiple gene sequences

To gain insights into the relationships among anostracan families, molecular phylogenetic analyses were performed on nuclear (28S D1-D3 ribosomal DNA) and mitochondrial (16S rDNA, COI) gene regions for representatives of seven families and an outgroup. Data matrices used in the analyses included 951 base pairs (bp) of aligned sequences for 28S, 465 bp for 16S, and 658 bp (219 amino acids) for COI. Maximum-parsimony and maximum-likelihood methods were used to construct phylogenetic trees, enabling the evaluation of both previous hypotheses of taxonomic relationships among families based on morphology, and of the relative merits of independent versus simultaneous analyses of multiple data sets for phylogeny construction. Data from various combinations of the gene regions produced relatively congruent patterns of phylogenetic affinity. In most analyses, two monophyletic groups were resolved: one cluster included the families Polyartemiidae, Chirocephalidae, Branchinectidae, Streptocephalidae, and Thamnocephalidae, while the other contained the Artemiidae and Branchipodidae. Comparative analyses showed that combining gene regions in a single matrix generally resulted in increased resolution and support for each cluster relative to those obtained from single-gene analyses. Statistical tests demonstrated that morphology-based hypotheses of relationships among families had poorer support than those determined from molecular data, reflecting the homoplasy in characters used to differentiate families.