Speciation in the Artemia genus: mitochondrial DNA analysis of bisexual and parthenogenetic brine shrimps

Mitogenomic phylogeny and divergence time estimation of Artemia Leach, 1819 (Branchiopoda: Anostraca) with emphasis on parthenogenetic lineages

The brine shrimp Artemia, a crustacean adapted to the extreme conditions of hypersaline environments, comprises nine regionally distributed sexual species scattered (island-like) over heterogeneous environments and asexual (parthenogenetic) lineages with different ploidies. Such sexual and asexual interaction within the genus raises questions regarding the origin and time of divergence of both sexual species and asexual lineages, including the persistence of the latter over time, a problem not yet clarified using single mitochondrial and nuclear markers. Based on the complete mitochondrial genome of all species and parthenogenetic lineages, this article first describes the mitogenomic characteristics (nucleotide compositions, genome mapping, codon usage, and tRNA secondary structure) of sexual species and asexual types and, secondly, it provides a comprehensive updated phylogenetic analysis. Molecular dating and geographical evidence suggest that the ancestral Artemia taxon originated in ca. 33.97 Mya during the Paleogene Period. The mitogenomic comparisons suggest that the common ancestor of diploid and triploid parthenogenetic lineages (ca. 0.07 Mya) originated from a historical ancestor (ca. 0.61 Mya) in the Late Pleistocene. Additionally, the common ancestor of tetraploid and pentaploid parthenogenetic lineages (ca. 0.05 Mya) diverged from a historical maternal ancestor with A. sinica (ca. 0.96 Mya) in the early Pleistocene. The parthenogenetic lineages do not share a direct ancestor with any sexual species. The Asian clade ancestor diverged more recently (ca. 14.27 Mya, Middle Miocene). The mitogenomic characteristics, maternal phylogenetic tree, and especially divergence time prove that A. monica and A. franciscana are two biological species.

Conflicts in Mitochondrial Phylogenomics of Branchiopoda, with the First Complete Mitogenome of Laevicaudata (Crustacea: Branchiopoda)

Conflicting phylogenetic signals are pervasive across genomes. The potential impact of such systematic biases may be reduced by phylogenetic approaches accommodating for heterogeneity or by the exclusive use of homoplastic sites in the datasets. Here, we present the complete mitogenome of Lynceus grossipedia as the first representative of the suborder Laevicaudata. We employed a phylogenomic approach on the mitogenomic datasets representing all major branchiopod groups to identify the presence of conflicts and concordance across the phylogeny. We found pervasive phylogenetic conflicts at the base of Diplostraca. The homogeneity of the substitution pattern tests and posterior predictive tests revealed a high degree of compositional heterogeneity among branchiopod mitogenomes at both the nucleotide and amino acid levels, which biased the phylogenetic inference. Our results suggest that Laevicaudata as the basal clade of Phyllopoda was most likely an artifact caused by compositional heterogeneity and conflicting phylogenetic signal. We demonstrated that the exclusive use of homoplastic site methods combining the application of site-heterogeneous models produced correct phylogenetic estimates of the higher-level relationships among branchiopods.

The complete mitochondrial genome of the Arctic fairy shrimp Branchinectapaludosa (Müller, 1788) (Anostraca, Branchinectidae) from Sirius Passet, North Greenland

Here we report the complete mitochondrial genome of the Arctic fairy shrimp, Branchinectapaludosa (Müller, 1788) (Anostraca, Branchinectidae), which was collected in the High Arctic of North Greenland. A complete 16,059 bp mitochondrion of B.paludosa was sequenced and assembled with the Illumina next generation sequencing platform. The B.paludosa mitogenome contains 13 PCGs, 22 tRNAs and 2 rRNA genes that are commonly observed in most metazoans and shows the conserved gene arrangement pattern of Anostraca. Our results of the phylogenomic analysis are consistent with the previous phylogenetic relationship, based on nuclear 18S ribosomal DNA. The B.paludosa mitogenome will be useful for understanding the geographical distribution and phylogenetic relationship of anostracans.

The complete mitochondrial genome of the Antarctic fairy shrimp Branchinectagaini Daday, 1910 (Branchiopoda, Anostraca, Branchinectidae)

The complete mitochondrial genome of Antarctic fairy shrimp Branchinectagaini Daday, 1910 was sequenced, assembled and annotated using next-generation sequencing technology. The mitogenome of B.gaini is circular at 15,536 bp in length, consisting of 13 protein-coding genes, 23 tRNAs, two rRNAs and two major non-coding regions. In particular, there are two tRNA^Gly genes and one non-coding region between these two tRNA^Gly genes. A phylogenetic tree was constructed using concatenated amino acid sequences of 13 protein-coding genes. It reveals that B.gaini is clustered with the Anostraca group within the Branchiopoda clade. This study helps us understand the evolution of Anostraca.

Comparative Mitogenomic Analyses and New Insights into the Phylogeny of Thamnocephalidae (Branchiopoda: Anostraca)

Thamnocephalidae, a family of Anostraca which is widely distributed on all continents of the world except Antarctica, currently consists of six genera and approximately 63 recognized species. The relationships among genera in Thamnocephalidae and the monophyly of Thamnocephalidae, determined using morphological characteristics or gene markers, remain controversial. In order to address the relationships within Thamnocephalidae, we sequenced Branchinella kugenumaensis mitogenomes and conducted a comparative analysis to reveal the divergence across mitogenomes of B. kugenumaensis. Using newly obtained mitogenomes together with available Anostracan genomic sequences, we present the most complete phylogenomic understanding of Anostraca to date. We observed high divergence across mitogenomes of B. kugenumaensis. Meanwhile, phylogenetic analyses based on both amino acids and nucleotides of the protein-coding genes (PCG) provide significant support for a non-monophyletic Thamnocephalidae within Anostraca, with Asian Branchinella more closely related to Streptocephalidae than Australian Branchinella. The phylogenetic relationships within Anostraca were recovered as follows: Branchinectidae + Chirocephalidae as the basal group of Anostraca and halophilic Artemiidae as a sister to the clade Thamnocephalidae + Streptocephalidae. Both Bayesian inference (BI)- and maximum likelihood (ML)-based analyses produced identical topologies.

Complete mitochondrial genomes of three fairy shrimps from snowmelt pools in Japan

Background

Fairy shrimps belong to order Anostraca, class Branchiopoda, subphylum Crustacea, and phylum Arthropoda. Three fairy shrimp species (Eubranchipus uchidai, E. asanumai, and E. hatanakai) that inhabit snowmelt pools are currently known in Japan. Whole mitochondrial genomes are useful genetic information for conducting phylogenetic analyses. Mitochondrial genome sequences for Branchiopoda members are gradually being collated.

Results

Six whole mitochondrial genomes from the three Eubranchipus species are presented here. Eubranchipus species share the anostracan pattern of gene arrangement in their mitochondrial genomes. The mitochondrial genomes of the Eubranchipus species have a higher GC content than those of other anostracans. Accelerated substitution rates in the lineage of Eubranchipus species were observed.

Conclusion

This study is the first to obtain whole mitochondrial genomes for Far Eastern Eubranchipus species. We show that the nucleotide sequences of cytochrome oxidase subunit I and the 16S ribosomal RNA of E. asanumai presented in a previous study were nuclear mitochondrial DNA segments. Higher GC contents and accelerated substitution rates are specific characteristics of the mitochondrial genomes of Far Eastern Eubranchipus. The results will be useful for further investigations of the evolution of Anostraca as well as Branchiopoda.

Settling taxonomic and nomenclatural problems in brine shrimps, Artemia (Crustacea: Branchiopoda: Anostraca), by integrating mitogenomics, marker discordances and nomenclature rules

High morphological plasticity in populations of brine shrimp subjected to different environmental conditions, mainly salinity, hindered for centuries the identification of the taxonomic entities encompassed within Artemia. In addition, the mismatch between molecular and morphological evolution rates complicates the characterization of evolutionary lineages, generating taxonomic problems. Here, we propose a phylogenetic hypothesis for Artemia based on two new complete mitogenomes, and determine levels of congruence in the definition of evolutionary units using nuclear and mtDNA data. We used a fossil of Artemia to calibrate the molecular clock and discuss divergence times within the genus. The hypothesis proposed herein suggests a more recent time frame for lineage splitting than previously considered. Phylogeographic analyses were performed using GenBank available mitochondrial and nuclear markers. Evidence of gen e flow, identified through discordances between nuclear and mtDNA markers, was used to reconsider the specific status of some taxa. As a result, we consider Artemia to be represented by five evolutionary units: Southern Cone, Mediterranean-South African, New World, Western Asian, and Eastern Asian Lineages. After an exhaustive bibliographical revision, unavailable names for nomenclatural purposes were discarded. The remaining available names have been assigned to their respective evolutionary lineage. The proper names for the evolutionary units in which brine shrimps are structured remain as follows: Artemia persimilis Piccinelli & Prosdocimi, 1968 for the Southern Cone Lineage, Artemia salina (Linnaeus, 1758) for the Mediterranean-SouthAfrican Lineage, Artemia urmiana Günther, 1899 for the Western Asian Lineage, and Artemia sinica Cai, 1989 for the Eastern Asian Lineage. The name Artemia monica Verrill, 1869 has nomenclatural priority over A. franciscana Kellogg, 1906 for naming the New World Lineage. New synonymies are proposed for A. salina (= C. dybowskii Grochowski, 1896 n. syn., and A. tunisiana Bowen & Sterling, 1978 n. syn.), A. monica (= A. franciscana Kellogg, 1906 n. syn., and A. salina var. pacifica Sars, 1904 n. syn.); A. urmiana (= B. milhausenii Fischer de Waldheim, 1834 n. syn., A. koeppeniana Fischer, 1851 n. syn., A. proxima King, 1855 n. syn., A. s. var. biloba Entz, 1886 n. syn., A. s. var. furcata Entz, 1886 n. syn., A. asiatica Walter, 1887 n. syn., A. parthenogenetica Bowen & Sterling, 1978 n. syn., A. ebinurica Qian & Wang, 1992 n. syn., A. murae Naganawa, 2017 n. syn., and A. frameshifta Naganawa & Mura, 2017 n. syn.). Internal deep nuclear structuring within the A. monica and A. salina clades, might suggest the existence of additional evolutionary units within these taxa.

Mitochondrial genome diversity and evolution in Branchiopoda (Crustacea)

BACKGROUND

The crustacean class Branchiopoda includes fairy shrimps, clam shrimps, tadpole shrimps, and water fleas. Branchiopods, which are well known for their great variety of reproductive strategies, date back to the Cambrian and extant taxa can be mainly found in freshwater habitats, also including ephemeral ponds. Mitochondrial genomes of the notostracan taxa Lepidurus apus lubbocki (Italy), L. arcticus (Iceland) and Triops cancriformis (an Italian and a Spanish population) are here characterized for the first time and analyzed together with available branchiopod mitogenomes.

RESULTS

Overall, branchiopod mitogenomes share the basic structure congruent with the ancestral Pancrustacea model. On the other hand, rearrangements involving tRNAs and the control region are observed among analyzed taxa. Remarkably, an unassigned region in the L. apus lubbocki mitogenome showed a chimeric structure, likely resulting from a non-homologous recombination event between the two flanking trnC and trnY genes. Notably, Anostraca and Onychocaudata mitogenomes showed increased GC content compared to both Notostraca and the common ancestor, and a significantly higher substitution rate, which does not correlate with selective pressures, as suggested by dN/dS values.

CONCLUSIONS

Branchiopod mitogenomes appear rather well-conserved, although gene rearrangements have occurred. For the first time, it is reported a putative non-homologous recombination event involving a mitogenome, which produced a pseudogenic tRNA sequence. In addition, in line with data in the literature, we explain the higher substitution rate of Anostraca and Onychocaudata with the inferred GC substitution bias that occurred during their evolution.

Long-term prevalence data reveals spillover dynamics in a multi-host (Artemia), multi-parasite (Microsporidia) community

In the study of multi-host parasites, it is often found that host species contribute asymmetrically to parasite transmission. Yet in natural populations, identifying which hosts contribute to parasite transmission and maintenance is a recurring challenge. Here, we approach this issue by taking advantage of natural variation in the composition of a host community. We studied the brine shrimps Artemia franciscana and Artemia parthenogenetica and their microsporidian parasites Anostracospora rigaudi and Enterocytospora artemiae. Previous laboratory experiments had shown that each host can transmit both parasites, but could not predict their actual contributions to the parasites' maintenance in the field. To resolve this, we gathered long-term prevalence data from a metacommunity of these species. Metacommunity patches could contain either or both of the Artemia host species, so that the presence of the hosts could be linked directly to the persistence of the parasites. First, we show that the microsporidian A. rigaudi is a spillover parasite: it was unable to persist in the absence of its maintenance host A. parthenogenetica. This result was particularly striking, as A. rigaudi displayed both high prevalence (in the field) and high infectivity (when tested in the laboratory) in both hosts. Moreover, the seasonal presence of A. parthenogenetica imposed seasonality on the rate of spillover, causing cyclical pseudo-endemics in the spillover host A. franciscana. Second, while our prevalence data was sufficient to identify E. artemiae as either a spillover or a facultative multi-host parasite, we could not distinguish between the two possibilities. This study supports the importance of studying the community context of multi-host parasites, and demonstrates that in appropriate multi-host systems, sampling across a range of conditions and host communities can lead to clear conclusions about the drivers of parasite persistence.

Morphological and molecular clues for recording the first appearance of Artemia franciscana () in Egypt

Artemia franciscana is a native species to the New World, and became an exotic species to most parts of the world. The Egyptian hypersaline, continental Qaroun Lake (Fayoum Governorate, Middle of Egypt) is subjected to a gradually increasing salinity rates that approximate or exceed these of seawater. Artemia populations there are known to be parthenogenetic. Yet, these populations started to exhibit abnormal morphologies. Therefore, Qaroun Lake samples of Artemia were subjected to several morphological, biometric, and molecular phylogenetic analyses for accurate species identification and phylogeographic origin approximation. These analyses revealed the existence of the alien sexual species of brine shrimp A. franciscana in Qaroun Lake. The characteristics of the subspherical frontal knob with several spines on the top, ovisac lateral triangular lobe on both sides and its projection together with the biometrics confirmed this species morphotype. DNA barcoding and other molecular analyses based on PCR-based amplification and sequencing of the barcode region of the cytochrome oxidase subunit I gene (COI) exhibited that all the collected samples belong to five haplotypes. Egyptian A. franciscana COI sequences phylogeny and pairwise distances analysis exhibited closer proximity to Latin American strains than to the Northern American ones. A. franciscana presence may be ascribed to the migratory birds present in Qaroun Lake protectorate, since no marine aquaculture activity in Qaroun Lake is known. Therefore, and for the best of our knowledge, this is the first record of the invasive A. franciscana in Egypt.

PHYLOGENY AND EVOLUTION OF PARTHENOGENETIC WEEVILS OF THE ARAMIGUS TESSELLATUS SPECIES COMPLEX (COLEOPTERA: CURCULIONIDAE: NAUPACTINI): EVIDENCE FROM MITOCHONDRIAL DNA SEQUENCES

Molecular-phylogenetic studies of parthenogenetic animals have been a valuable recent addition to the literature on the evolutionary biology of sex. By illuminating the origins and ages of parthenogenetic lineages, such studies can help to define the temporal scale at which selection acts against parthenogenetic lineages, as well as provide an essential framework for further study. Although parthenogenetic weevils have played an important role in cytogenetic and protein-electrophoretic studies of parthenogenesis, they have not previously been subjects of DNA-based molecular-phylogenetic study. A mitochondrial DNA study of Aramigus tessellatus, a species complex of weevils native to South America, indentified 12 distinct (1-9% divergent) maternal lineages, of which 2 represent sexual populations, while at least 9 represent parthenogenetic lineages. These lineages partially correspond to lineages previously recognized by morphological differences. Phylogenetic analysis found 14 most parsimonious trees, according to which parthenogenesis appears to have arisen 3-7 times. There is a monophyletic group of lineages (the "brown clade"), having up to 4.5% sequence divergence within it, which may be primitively parthenogenetic and over 2 million years old.

Phylogeographic analyses of bacterial endosymbionts in fig homotomids (Hemiptera: Psylloidea) reveal codiversification of both primary and secondary endosymbionts

While obligate primary (P-) endosymbionts usually cospeciate with their insect hosts, less is known about codiversification of secondary (S-) endosymbionts that are generally considered facultative. Typically, insects of the superfamily Psylloidea harbour one P- (Carsonella) and at least one S-endosymbiont, thought to compensate for Carsonella genome reduction. Most codiversification studies have used phylogenies of psyllids and their endosymbionts across and within host families or genera, but few have explored patterns within species. We focussed on P- and S-endosymbionts of three Mycopsylla (Homotomidae) species to explore whether they have congruent phylogenies and within-species geographic structures. The P-endosymbiont Carsonella, a S-endosymbiont and Wolbachia all had 100% prevalence, while Arsenophonus was only found in one species at low prevalence. Congruent phylogenies of Mycopsylla and P-endosymbionts across populations and species support strict cospeciation. S-endosymbiont phylogenies were also congruent across host species but low genetic variation in the S-endosymbiont was not correlated with host phylogeography, possibly due to a shorter evolutionary association. Between species, Wolbachia and Mycopsylla phylogenies were incongruent, probably due to horizontal transmission events. Our study is the first to explore endosymbionts of Mycopsylla and further supports the codivergence of Psylloidea hosts and P-endosymbionts, with obligate host interactions for both P- and S-endosymbionts.

Mitochondrial genome sequences of Artemia tibetiana and Artemia urmiana: assessing molecular changes for high plateau adaptation

Brine shrimps, Artemia (Crustacea, Anostraca), inhabit hypersaline environments and have a broad geographical distribution from sea level to high plateaus. Artemia therefore possess significant genetic diversity, which gives them their outstanding adaptability. To understand this remarkable plasticity, we sequenced the mitochondrial genomes of two Artemia tibetiana isolates from the Tibetan Plateau in China and one Artemia urmiana isolate from Lake Urmia in Iran and compared them with the genome of a low-altitude Artemia, A. franciscana. We compared the ratio of the rate of nonsynonymous (Ka) and synonymous (Ks) substitutions (Ka/Ks ratio) in the mitochondrial protein-coding gene sequences and found that atp8 had the highest Ka/Ks ratios in comparisons of A. franciscana with either A. tibetiana or A. urmiana and that atp6 had the highest Ka/Ks ratio between A. tibetiana and A. urmiana. Atp6 may have experienced strong selective pressure for high-altitude adaptation because although A. tibetiana and A. urmiana are closely related they live at different altitudes. We identified two extended termination-associated sequences and three conserved sequence blocks in the D-loop region of the mitochondrial genomes. We propose that sequence variations in the D-loop region and in the subunits of the respiratory chain complexes independently or collectively contribute to the adaptation of Artemia to different altitudes.

Complete mitochondrial genome of a Chinese scorpionMesobuthus martensii(Chelicerata, Scorpiones, Buthidae)

The complete mitochondrial genome (15,034 bp) of a Chinese scorpion Mesobuthus martensii (Buthidae) was sequenced and characterized in detail. The genome contains 13 protein-coding genes, 21 transfer RNA genes, two ribosomal RNA genes and a large non-coding region ( = CR). Its gene arrangement pattern is identical to that of Limulus polyphemus (Chelicerata, Xiphosura), with the exceptions of the tRNA(Glu)-tRNA(Ile)-tRNA(Met) (Q-I-M) arrangement and tRNA(Asp)-loss. Additional interesting features are found and discussed: high frequency of Leu(UUG) codon use, low A+T content of the genome (66.75%), and six repeat units (five 60-nt-long and one 58-nt-long repeats) in the 998-nt CR. Bayesian analysis based on amino acid sequences of the 12 proteincoding genes (excluding ATP8) reveals that the family Buthidae (Order Scorpiones) and the class Arachnida form strong monophyletic groups within Chelicerata, respectively. It indicated that the scorpions are the most ancestral arachnids.

First divergence time estimate of spiders, scorpions, mites and ticks (subphylum: Chelicerata) inferred from mitochondrial phylogeny

Spiders, scorpions, mites and ticks (chelicerates) form one of the most diverse groups of arthropods on land, but their origin and times of diversification are not yet established. We estimated, for the first time, the molecular divergence times for these chelicerates using complete mitochondrial sequences from 25 taxa. All mitochondrial genes were evaluated individually or after concatenation. Sequences belonging to three missing genes (ND3, 6, and tRNA-Asp) from three taxa, as well as the faster-evolving ribosomal RNAs (12S and 16S), tRNAs, and the third base of each codon from 11 protein-coding genes (PCGs) (COI-III, CYTB, ATP8, 6, ND1-2, 4L, and 4-5), were identified and removed. The remaining concatenated sequences from 11 PCGs produced a completely resolved phylogenetic tree and confirmed that all chelicerates are monophyletic. Removing the third base from each codon was essential to resolve the phylogeny, which allowed deep divergence times to be calculated using three nodes calibrated with upper and lower priors. Our estimates indicate that the orders and classes of spiders, scorpions, mites, and ticks diversified in the late Paleozoic, much earlier than previously reported from fossil date estimates. The divergence time estimated for ticks suggests that their first land hosts could have been amphibians rather than reptiles. Using molecular data, we separated the spider-scorpion clades and estimated their divergence times at 397 +/- 23 million years ago. Algae, fungi, plants, and animals, including insects, were well established on land when these chelicerates diversified. Future analyses, involving mitochondrial sequences from additional chelicerate taxa and the inclusion of nuclear genes (or entire genomes) will provide a more complete picture of the evolution of the Chelicerata, the second most abundant group of animals on earth.

The complete mitochondrial genome of Symphylella sp. (Myriapoda: Symphyla): Extensive gene order rearrangement and evidence in favor of Progoneata

We determined the complete 14,667bp mitochondrial DNA sequence of Symphylella sp., the first representative of the Scolopendrellidae (Arthropoda: Myriapoda: Symphyla). With respect to the ancestral arthropod mitochondrial gene order, two protein-coding genes, the rRNAs and 10 of the tRNAs appear to be rearranged. This rearrangement is novel in the arthropods and genes with identical transcriptional polarity are clustered except for trnE, trnN and putative control region (CR), resembling two previously reported diplopod genomes. A duplication/loss (random and non-random)-recombination model was proposed to account for the generation of the gene order in Symphylella sp. All phylogenetic analysis yielded strong support for a clade of Symphyla plus Diplopoda, i.e., Progoneata. However, the phylogenetic position of Myriapoda within Arthropoda remains unclear. The amino acid dataset gives strong support for an affinity to Pancrustacea, while the nucleotide dataset weakly supports Myriapoda grouped with Chelicerata.

Cloning and expression analysis of p26 gene in Artemia sinica

The protein p26 is a small heat shock protein that functions as a molecular chaperone to protect embryos by preventing irreversible protein damage during embryonic development. A 542 bp fragment of the p26 gene was cloned and sequenced. The fragment encoded 174 amino acid residues and the amino acid sequence contained the alpha-crystallin domain. Phylogenetic analysis showed that eight Artemia populations were divided into four major groups. Artemia sinica (YC) belonged to the East Asia bisexual group. Expression of the p26 gene at different developmental stages of A. sinica was quantified using real-time quantitative polymerase chain reaction followed by cloning and sequencing. The relationship between the quantity of p26 gene expression and embryonic development was analyzed. The results indicated that massive amounts of p26 were expressed during the development of A. sinica. At the developmental stage of 0 h, A. sinica expressed the highest level of p26. As development proceeded, expression levels of the p26 gene reduced significantly. There was a small quantity of p26 gene expression at the developmental stages of 16 h and 24 h. We concluded that p26 might be involved in protecting the embryo from physiological stress during embryonic development.

Molecular phylogenetics and asexuality in the brine shrimp Artemia

Explaining cases of long-term persistence of parthenogenesis has proven an arduous task for evolutionary biologists. Interpreting sexual-asexual interactions though has recently advanced owing to methodological design, increased taxon sampling and choice of model organisms. We inferred the phylogeny of Artemia, a halophilic branchiopod genus of sexual and parthenogenetic forms with cosmopolitan distribution, marked geographic patterns and ecological partitioning. Joint analysis of newly derived ITS1 sequences and 16S RFLP markers from global isolates indicates significant interspecific divergence as well as pronounced diversity for parthenogens, matching that of sexual ancestors. Maximum parsimony, maximum likelihood, and Bayesian methods were largely congruent in reconstructing the phylogeny of the genus. Given the current sampling, at least four independent origins of parthenogenesis are deduced. Molecular clock calibrations based on biogeographic landmarks indicate that the lineage leading to A. persimilis diverged from the common ancestor of all Artemia species between 80 and 90 MYA at the time of separation of Africa from South America, whereas parthenogenesis first appeared at least 3 MYA. Common mitochondrial DNA haplotypes delineate A. urmiana and A. tibetiana as possible maternal parents of several clonal lineages. A novel topological placement of A. franciscana as a sister clade to all Asian Artemia and parthenogenetic forms is proposed and also supported by ITS1 length and other existing data.

Recent and ancient asexuality in Timema walkingsticks

Determining the evolutionary age of asexual lineages should help in inferring the temporal scale under which asexuality and sex evolve and assessing selective factors involved in the evolution of asexuality. We used 416 bp of the mitochondrial COI gene to infer phylogenetic relationships of virtually all known Timema walkingstick species, including extensive intraspecific sampling for all five of the asexuals and their close sexual relatives. The asexuals T. douglasi and T. shepardii were very closely related to each other and evolutionarily young (less than 0.5 million years old). For the asexuals T. monikensis and T. tahoe, evidence for antiquity was weak since only one population of each was sampled, intraspecific divergences were low, and genetic distances to related sexuals were high: maximum-likelihood molecular-clock age estimates ranged from 0.26 to 2.39 million years in T. monikensis and from 0.29-1.06 million years in T. tahoe. By contrast, T. genevieve was inferred to be an ancient asexual, with an age of 0.81 to 1.42 million years. The main correlate of the age of asexual lineages was their geographic position, with younger asexuals being found further north.

Speciation in ancient cryptic species complexes: evidence from the molecular phylogeny of Brachionus plicatilis (Rotifera)

Continental lake-dwelling zooplanktonic organisms have long been considered cosmopolitan species with little geographic variation in spite of the isolation of their habitats. Evidence of morphological cohesiveness and high dispersal capabilities support this interpretation. However, this view has been challenged recently as many such species have been shown either to comprise cryptic species complexes or to exhibit marked population genetic differentiation and strong phylogeographic structuring at a regional scale. Here we investigate the molecular phylogeny of the cosmopolitan passively dispersing rotifer Brachionus plicatilis (Rotifera: Monogononta) species complex using nucleotide sequence variation from both nuclear (ribosomal internal transcribed spacer 1, ITS1) and mitochondrial (cytochrome c oxidase subunit I, COI) genes. Analysis of rotifer resting eggs from 27 salt lakes in the Iberian Peninsula plus lakes from four continents revealed nine genetically divergent lineages. The high level of sequence divergence, absence of hybridization, and extensive sympatry observed support the specific status of these lineages. Sequence divergence estimates indicate that the B. plicatilis complex began diversifying many millions of years ago, yet has showed relatively high levels of morphological stasis. We discuss these results in relation to the ecology and genetics of aquatic invertebrates possessing dispersive resting propagules and address the apparent contradiction between zooplanktonic population structure and their morphological stasis.

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.

Wolbachia diversity in the Porcellionides pruinosus complex of species (Crustacea: Oniscidea): evidence for host-dependent patterns of infection

Porcellionides pruinosus is a cosmopolitan woodlouse. It is known to exhibit patterns of geographical variation between populations, and has been suspected to consist of several very closely related species. This species was found to carry Wolbachia endosymbionts, alpha-proteobacteria which are known to modify the reproduction of their crustacean hosts by inducing cytoplasmic incompatibility or feminization. In the P. pruinosus complex, Wolbachia induced feminization, but two different patterns of infection were reported: high prevalence and presence of infected males, or low prevalence and absence of infected males. In this study we investigated nine populations described as Porcellionides pruinosus carrying different Wolbachia strains, with the aim of investigating the possibility of coevolution between symbionts and hosts. Molecular analyses were carried out on both Wolbachia DNA (wsp gene) and host mitochondrial DNA (mt LSU rDNA). We show that (1) the nine host populations exhibited a high genetic differentiation so that they have to be split into two sibling species; (2) three different Wolbachia strains were found in the host complex but were not species-specific; and (3) the pattern of infection by Wolbachia was more host-dependent than symbiont-dependent.

International study on Artemia. LXII. Genomic relationships between Artemia franciscana and A. persimilis, inferred from chromocentre numbers

Chromocentres, i.e. heavily stainable heterochromatic areas with highly repetitive DNA (a 130-bp repeat in the order of 6 x 105 copies per haploid genome) observed in the resting nucleus, are a reliable taxonomic trait and a good marker for speciation in Artemia. This chromosome marker was evaluated in populations of two New World sibling species: A. franciscana, from North, Central and South America, and A. persimilis from Argentina. Artemia persimilis showed the characteristically low average chromocentre numbers (<3) of the ancestral Mediterranean species, whilst higher numbers (>5) were seen in A. franciscana. The increase in chromocentre numbers from A. persimilis to A. franciscana occurred in a rather steady, continuous geographical pattern with an overlapping zone between both species. A north-south steady latitudinal decline of chomocentres was observed towards the equator in A. franciscana populations from the northern hemisphere, whilst the trait increased from this point towards southern latitudes with a population in Chilean Patagonia (below 52 degrees S) exhibiting the highest chromocentre frequency. Since chromocentres are positively correlated with repetitive DNA content, differentiation between the species through the amplification of heterochromatin or an increased DNA content in A. franciscana (although this species has a lower diploid chromosome number), might nevertheless reflect chromosomal reorganizations between the species. Chromocentres, together with other traits, reveal a probable hybrid zone between the two species in Chilean Patagonia.

Sequence and organization of the mitochondrial genome of the Chagas disease vector, Triatoma dimidiata

The 17 019 bp mitochondrial genome of Triatoma dimidiata is composed of thirteen protein coding sequences, twenty-two tRNAs, small and large ribosomal units, and a control region. The gene order and orientation are identical to that of Drosophila yakuba. The nucleotide composition is biased toward adenine and thymine (69.5% A + T). The 2.1 kb putative control region, known as the A + T rich region in most insects, has an A + T bias of 66%, but contains a 400 bp sequence that is 77.5% A + T and two other distinct regions: (1) one with a lower A + T bias (60.1%) and (2) a region of eight tandem repeat units. The identified 1.4 kb nuclear copy of mitochondrial sequences encompasses the string of Gs and the beginning of the cytochrome c oxidase 1 gene but lacks the 1.8 kb region spanning the eight tandem repeats and the 5' end of the NADH dehydrogenase subunit II gene.

Phylogenetic analysis of mitochondrial LSU rRNA in oniscids

The phylogenetic relationships among oniscids (Crustacea, Isodopa) remain contradictory despite numerous morphological studies. We have investigated them using molecular data. Partial sequences of the mitochondrial LSU rRNA gene were obtained from 42 species of aquatic and terrestrial crustaceans from 31 genera. This gene provided well-supported information, notwithstanding the high taxonomic level of this study, indicating a useful amount of variation despite the noise due to multiple substitutions. The phylogenetic inferences demonstrated that a) Crinocheta and Synocheta sections are monophyletic and sister-groups, b) Ligiidae and Tylidae representatives are in a basal position compared to other oniscids, c) Helleria brevicornis, the only representative of the Helleriinae subfamily, has undergone different evolution, d) the relationships between aquatic isopods and ancient groups of Oniscidea are not resolved, probably due to fast radiation not discriminated by the molecular phylogeny.

Persistence of asexuality through mixed reproduction in Eucypris virens (Crustacea, Ostracoda)

The ostracod species Eucypris virens exhibits geographical parthenogenesis, with rare sexual populations in southern Europe and widespread asexual populations elsewhere. DNA sequence data from the nuclear ITS1 and mitochondrial COI regions have been used to estimate genetic variabilities and reconstruct phylogenies. The observed divergence was exceptionally high, with intraspecific maxima of 10.3% (ITS1) and 20.9% (COI) among European lineages, levels reported for interspecific comparisons of other taxa. Phylogenetic reconstructions reveal multiple origins of asexual clones from sexual populations. However, we argue that such data can only provide a lower limit on the number of origins of asexual reproduction, and an upper limit on the age of asexual lineages. Congruence between gene trees for different loci can provide support for the inference of long-term apomictic reproduction. Nuclear and mitochondrial data differ in their placement of some asexual clones, possibly indicating that genetic exchange has taken place between sexual and asexual lineages. Such intraspecific hybridization is one route to combine the benefits of both reproductive modes, and it might explain how asexuality managed to persist in E. virens even in long, evolutionary terms.

Slow molecular evolution in an ancient asexual ostracod

Genetic variability of the non-marine ostracod species Darwinula stevensoni was estimated by sequencing part of the nuclear and the mitochondrial genome. As Darwinulidae are believed to be ancient asexuals, accumulation of mutations should have occurred, both between alleles within lineages and between lineages, during the millions of years of parthenogenetic reproduction. However, our sequence data show the opposite: no variability in the nuclear ITS1 region was observed within or among individuals of D. stevensoni , despite sampling a geographical range from Finland to South Africa. Lack of allelic divergence might be explained by concerted evolution of rDNA repeats. Homogeneity among individuals may be caused either by slow molecular evolution in ITS1 or by a recent selective sweep. Variability of mitochondrial cytochrome oxidase (COI) was similar to intraspecific levels in other invertebrates, thus weakening the latter hypothesis. Calibrating interspecific, genetic divergences among D. stevensoni and other Darwinulidae using their fossil record enabled us to estimate rates of molecular evolution. Both COI and ITS1 evolve half as fast, at most, in darwinulids as in other invertebrates, and molecular evolution has significantly slowed down in ITS1 of D. stevensoni relative to other darwinulids. A reduced ITS1 mutation rate might explain this inconsistency between nuclear and mitochondrial evolution in D. stevensoni .

Taxonomic status and phylogenetic relationships of some species of the genus Gammarus (Crustacea, Amphipoda) deduced from mitochondrial DNA sequences

Mitochondrial DNA nucleotide sequence analysis was used to determine the extent of genetic differentiation among six species of the genus Gammarus (Crustacea, Amphipoda) common in France (G. fossarum, G. pulex, G. lacustris, G. wautieri, G. roeseli, G. marinus). From 23 different populations, 104 specimens were compared to examine their taxonomic status and their phylogenetic relationships. Nucleotide sequences for a 376-base segment of the cytochrome c oxidase subunit I (COI) gene were obtained through direct sequencing of amplified DNA. The strong sequence variability found at the interspecific level, generally masked by morphological or ecological similarities, is consistent with that observed for other crustaceans and supports existing species classifications. A high level of genetic variation was observed at the intraspecific level, suggesting that the COI gene will be a powerful marker for amphipod population biology. Based on the mitochondrial population biology. Based on the mitochondrial DNA results, we propose a phylogenetic hypothesis of the relationships among the six studied species, which is discussed in biogeographic and evolutionary terms.