The scale of divergence: A phylogenetic appraisal of intercontinental allopatric speciation in a passively dispersed freshwater zooplankton genus

Balancing selection and the functional effects of shared polymorphism in cryptic Daphnia species

The patterns of genetic variation within and between related taxa represent the genetic history of a species. Shared polymorphisms, loci with identical alleles across species, are of unique interest as they may represent cases of ancient selection maintaining functional variation post-speciation. In this study, we investigate the abundance of shared polymorphism in the Daphnia pulex species complex. We test whether shared mutations are consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting, or convergent evolution. We analyzed over 2,000 genomes from North American and European D. pulex and several outgroup species to examine the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that while North American and European D. pulex diverged over ten million years ago, they retained tens of thousands of shared alleles. We found that the number of shared polymorphisms between North American and European D. pulex cannot be explained by hybridization or incomplete lineage sorting alone. Instead, we show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting young and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild. Ultimately, our findings provide insights into the genetic basis of adaptation and the maintenance of genetic diversity between species.

Across the highest mountain on earth: discordant phylogeographic patterns and recent dispersal of Tibetan stone loaches (Triplophysa) in the Himalayas

Phylogeographic congruence among co-distributed taxa is regarded as an inherent inference to vicariance events. Nonetheless, incongruent patterns of contemporary lineage divergence among taxa indicated that species differ in their response to common past events. To investigate the role of past events, ecological traits and lineage diversification time in shaping the contemporary phylogeographic patterns, comparative analyses were conducted for Tibetan stone loaches in the Himalayas using three gene markers and two ecological traits (depth of caudal peduncle in their length and presence/absence of posterior chamber of the air bladder). By a thorough sampling in two flanks of the Himalayas, the authors detected that phylogenetic breaks were spatially discordant and divergences of populations were also temporally asynchronous in co-distributed loaches. Estimated divergence time using fossil-calibrated node dating indicated that the Tibetan stone loaches colonised into the south flank of the Himalayas until the Pleistocene. The demographic expansions were also disconcerted between populations in north and south flanks, or east and west Himalayas. Ongoing gene flows between populations in north and south sides implied that the Himalayas do not strictly impede dispersal of cold-adapted species. The results highlight that the quaternary climatic oscillation, in conjunction with ecological traits and lineage diversification time, shaped contemporary phylogenetic patterns of stone loaches in the Himalayas and provide new insights into the biodiversity and composition of species in the Himalayas and surrounding region.

The hologenome of Daphnia magna reveals possible DNA methylation and microbiome-mediated evolution of the host genome

Properties that make organisms ideal laboratory models in developmental and medical research are often the ones that also make them less representative of wild relatives. The waterflea Daphnia magna is an exception, by both sharing many properties with established laboratory models and being a keystone species, a sentinel species for assessing water quality, an indicator of environmental change and an established ecotoxicology model. Yet, Daphnia's full potential has not been fully exploited because of the challenges associated with assembling and annotating its gene-rich genome. Here, we present the first hologenome of Daphnia magna, consisting of a chromosomal-level assembly of the D. magna genome and the draft assembly of its metagenome. By sequencing and mapping transcriptomes from exposures to environmental conditions and from developmental morphological landmarks, we expand the previously annotates gene set for this species. We also provide evidence for the potential role of gene-body DNA-methylation as a mutagen mediating genome evolution. For the first time, our study shows that the gut microbes provide resistance to commonly used antibiotics and virulence factors, potentially mediating Daphnia's environmental-driven rapid evolution. Key findings in this study improve our understanding of the contribution of DNA methylation and gut microbiota to genome evolution in response to rapidly changing environments.

Microcrustaceans (Cladocera and Copepoda) of the Boreal/Tropical Transition Zone in the Russian Far East: A Case Study of Species Associations in Three Large Lakes

The Far East of Russia is a region where boreal and tropical faunas mix; it is also a zone of cladoceran endemism. The present study aimed to compare a set of microcrustacean (Cladocera and Copepoda) associations in three large lakes of the Russian Far East: Khanka, Bolon, and Chukchagir. The associations of the microcrustaceans were identified based on the function of the discrete hypergeometric distribution. Many of the 108 taxa found here were unaffiliated with an association. Interestingly, the portion of taxa involved and “not involved” in species associations differed among geographic faunistic complexes. The rate of endemism was significantly higher among the taxa incorporated into the associations as compared to the “not involved” taxa. In all the lakes, there were large clusters of phytophilous species characteristic of the macrophyte zone (and its margins) and clusters characteristic of pelagic and sublittoral plankton. We found that in the three lakes, the microcrustaceans formed a set of functionally similar associations, but the taxonomic composition of each functional association was specific to each lake. We hypothesize that the composition of functional clusters reflects the history of colonization for each water body. That is, the founder effects and subsequent “monopolization” of habitats have affected species associations.

Daphnia japonica sp. nov. (Crustacea: Cladocera) an eastern Palearctic montane species with mitochondrial discordance

The Daphnia longispina complex (Crustacea: Cladocera) contains several keystone freshwater species such as D. longispina O.F. Müller (D. rosea Sars is a junior synonym), D. galeata Sars, D. cucullata Sars, and D. dentifera Forbes. The complex is common throughout the Holarctic, but there are several geographic regions where local forms have been assigned to European species names based on a superficial morphological resemblance. Here we examine the species status of a form that was previously assigned to D. rosea from a montane bog pond on Honshu, Japan. We used two nuclear non-coding loci (nDNA), mitochondrial sequences (the ND2 protein-coding region) and morphology for evidence. The mitochondrial gene evidence supported the existence of a divergent lineage that is more closely related to D. galeata than to D. dentifera. However, morphology and the nuclear DNA data indicated a lineage that is most closely related to D. dentifera. As our evidence supported the existence of a cohesive divergent lineage, we described a new species, Daphnia japonica sp. nov. Recognition of local and subalpine diversity in this group is critical as ongoing anthropogenic disturbance has been associated with introductions, local extirpations, and hybridization.

Ceriodaphnia (Cladocera: Daphniidae) in China: Lineage diversity, phylogeography and possible interspecific hybridization

The distribution and species/lineage diversity of freshwater invertebrate zooplankton remains understudied in China. Here, we explored the species/lineage diversity and phylogeography of Ceriodaphnia species across China. The taxonomy of this genus is under-explored. Seven morphospecies of Ceriodaphnia (C. cornuta, C. laticaudata, C. megops, C. pulchella, C. quadrangula, C. rotunda and C. spinata) were identified across 45 of 422 water bodies examined. Rather little morphological variation was observed within any single morphospecies regardless of country of origin. Nevertheless, we recognized that some or all of these morphospecies might represent species complexes. To investigate this, phylogenetic relationships within and among these morphospecies were investigated based on mitochondrial (partial cytochrome c oxidase subunit I gene) and nuclear (partial 28S rRNA gene) markers. The mitochondrial marker placed these populations in nine lineages corresponding to the morphospecies: C. laticaudata and C. pulchella were each represented by two lineages, suggesting that both are species complexes. The remaining five morphospecies were each represented by a single mtDNA lineage. Three of the nine mitochondrial lineages (belonging to C. pulchella, C. rotunda and C. megops) are newly reported and exhibited a restricted distribution within China. The nuclear-DNA phylogeny also recognized seven Ceriodaphnia taxa within China. We detected occasional mito-nuclear discordances in Ceriodaphnia taxa across China, suggesting interspecific introgression and hybridization. Our study contributes to an understanding of the species/lineage diversity of Ceriodaphnia, a genus with understudied taxonomy.

Mitochondrial Lineage Diversity and Phylogeography of Daphnia (Daphnia) (Crustacea: Cladocera) in North-East Russia

The variability of the 12S gene fragment of the mtDNA for taxa belonging to subgenus Daphnia (Daphnia) O.F. Müller, 1776 (Crustacea: Cladocera) in NE Russia is studied, and their phylogenetic analysis performed. We identified (based both on morphological and molecular data) nine species belonging to four species complexes, namely: (A) D. longispina s.l.: (1) D. longispina O.F. Müller, 1776; (2) D. dentifera Forbes, 1893; (3) D. galeata Sars, 1864; (4) D. umbra Taylor, Hebert et Colbourne, 1996; (B) D. cristata s.l.: (5) D. cristata Sars, 1862; (6) D. longiremis Sars, 1862; (C) D. curvirostris s.l.: (7) D. curvirostris Eylmann, 1887; (D) D. pulex s.l.: (8) D. pulex Leydig, 1860; (9) D. middendorffiana Fischer, 1851. Rare arcto-mountainous taxon D. umbra was found in the mountains of the Sakha (Yakutia) Republic for the first time. Species diversity in NE Asia is relatively low, and the most revealed taxa are trans-Beringian. We also performed a phylogeographic analysis of D. dentifera and D. pulex s.l., the two most common species in NE Russia. Our new data allow us to assume that the daphniids of NE Asia have undergone various evolutionary scenarios during the Pleistocene period: survival is within some local refugia, and re-colonization from these areas and from North America through the Beringian land bridge, etc. We agree with previous authors who revealed that the patterns in the studied species groups are relatively recent (of Late Pleistocene or even Holocene age), although the main phylogenetic daphniid lineages (mainly congruent with the biological species) are very old. Our results provide convincing evidence for the hypothesis that NE Russia is a very important source of modern haplotypic diversity for the cladocerans.

Multiple Recent Colonizations of the Australian Region by the Chydorus sphaericus Group (Crustacea: Cladocera)

Biotic introductions are an ongoing disruption for many ecosystems. For passively dispersed freshwater zooplankton, transcontinental introductions have been common but are poorly studied in the southern hemisphere. Here we assess the hypothesis of recent introduction for populations of the Chydorus sphaericus group (Crustacea: Cladocera) in Australia. We analyzed 254 sequences (63 original sequences) from the cytochrome oxidase I region of mitochondrial DNA of Chydorus sp., which included global representation. Three Australian populations were connected with separate clades in the northern hemisphere, suggesting multiple colonization events for Australia. The timescale of the divergences was consistent with recent (Quaternary) dispersal. As Australian populations are exposed to migrating birds from the northern hemisphere, both avian and anthropogenic sources are candidates for dispersal vectors. We concluded that recent cross-hemisphere dispersal in the Chydorus sphaericus group is more common than previously believed.

Daily feeding rhythm linked to microbiome composition in two zooplankton species

Host-associated microbial communities are impacted by external and within-host factors, i.e., diet and feeding behavior. For organisms known to have a circadian rhythm in feeding behavior, microbiome composition is likely impacted by the different rates of microbe introduction and removal across a daily cycle, in addition to any diet-induced changes in microbial interactions. Here, we measured feeding behavior and used 16S rRNA sequencing to compare the microbial community across a diel cycle in two distantly related species of Daphnia, that differ in their life history traits, to assess how daily feeding patterns impact microbiome composition. We find that Daphnia species reared under similar laboratory conditions have significantly different microbial communities. Additionally, we reveal that Daphnia have daily differences in their microbial composition that correspond with feeding behavior, such that there is greater microbiome diversity at night during the host’s active feeding phase. These results highlight that zooplankton microbiomes are relatively distinct and are likely influenced by host phylogeny.

Hybridization Dynamics and Extensive Introgression in the Daphnia longispina Species Complex: New Insights from a High-Quality Daphnia galeata Reference Genome

Hybridization and introgression are recognized as an important source of variation that influence adaptive processes; both phenomena are frequent in the genus Daphnia, a keystone zooplankton taxon in freshwater ecosystems that comprises several species complexes. To investigate genome-wide consequences of introgression between species, we provide here the first high-quality genome assembly for a member of the Daphnia longispina species complex, Daphnia galeata. We further resequenced 49 whole genomes of three species of the complex and their interspecific hybrids both from genotypes sampled in the water column and from single resting eggs extracted from sediment cores. Populations from habitats with diverse ecological conditions offered an opportunity to study the dynamics of hybridization linked to ecological changes and revealed a high prevalence of hybrids. Using phylogenetic and population genomic approaches, we provide first insights into the intra- and interspecific genome-wide variability in this species complex and identify regions of high divergence. Finally, we assess the length of ancestry tracts in hybrids to characterize introgression patterns across the genome. Our analyses uncover a complex history of hybridization and introgression reflecting multiple generations of hybridization and backcrossing in the Daphnia longispina species complex. Overall, this study and the new resources presented here pave the way for a better understanding of ancient and contemporary gene flow in the species complex and facilitate future studies on resting egg banks accumulating in lake sediment.

Phylogenetic analysis of eight species of Anomopoda based on transcriptomic and mitochondrial DNA sequences

Anomopoda is the widespread planktonic microcrustacean, which plays a crucial role in aquatic ecosystem. There are few studies about the evolutionary relationships among various Anomopoda basing on molecular data. In the present study, phylogenetic analysis of eight Anomopoda was carried out. Firstly, the culture system was developed to breed cladocerans. By using this system, eight species (Daphnia magna, D. pulex, D. sinensis, Ceriodaphnia reticulata, Moina micrura, Scapholeberis kingi, Simocephalus vetulus and Eurycercus lamellatus) were purified and cultured stably in the laboratory. Then, transcriptomic sequences and partial mitochondrial DNA sequences were both used to reconstruct the phylogenetic tree among 8 species. Transcriptomic sequences were sequenced on Illumina Hiseq 2500 platform. After assembly and annotation, transcriptomic sequences were spliced together and aligned for phylogenetic analysis. Basing on the orthologous genes derived from transcriptomic sequences, the phylogenetic analysis showed that 4 genera of Daphniidae were clustered into one group, and among the 4 genera, Ceriodaphnia was closer to Daphnia than Simocephalus, while Scapholeberis was farthest from other species. In addition, Eurycercidae was closer to Daphniidae than Moinidae. The phylogenetic trees based on both 12S rRNA and 16S rRNA sequences were similar with that based on transcriptomic sequences. Meanwhile, the phylogenetic tree based on 16S rRNA sequences was more suitable than that based on 12S rRNA sequences. These results suggested that the phylogenetic analysis basing on the transcriptomic sequences was available in cladocerans, which will help us to effectively understand the phylogenetic relationships among various cladocerans.

Speciation in Daphnia

The microcrustacean Daphnia is arguably one of the most studied zooplankton species, having a well understood ecology, life history, and a relatively well studied evolutionary history. Despite this wealth of knowledge, species boundaries within closely related species in this genus often remain elusive and the major evolutionary forces driving the diversity of daphniids remain controversial. This genus contains more than 80 species with multiple cryptic species complexes, with many closely related species able to hybridize. Here, we review speciation research in Daphnia within the framework of current speciation theory. We evaluate the role of geography, ecology, and biology in restricting gene flow and promoting diversification. Of the 253 speciation studies on Daphnia, the majority of studies examine geographic barriers (55%). While evidence shows that geographic barriers play a role in species divergence, ecological barriers are also probably prominent in Daphnia speciation. We assess the contribution of ecological and nonecological reproductive isolating barriers between closely related species of Daphnia and found that none of the reproductive isolating barriers are restricting gene flow completely. Research on reproductive isolating barriers has disproportionally focused on two species complexes, the Daphnia pulex and Daphnia longispina species complexes. Finally, we identify areas of research that remain relatively unexplored and discuss future research directions that build our understanding of speciation in daphniids.

Next-generation sequencing of DNA from resting eggs: signatures of eutrophication in a lake's sediment

Hatching resting stages of ecologically important organisms such as Daphnia from lake sediments, referred to as resurrection ecology, is a powerful approach to assess changes in alleles and traits over time. However, the utility of the approach is constrained by a few obstacles, including low and/or biased hatching among genotypes. Here, we eliminated such bottlenecks by investigating DNA sequences isolated directly (i.e. without hatching) from resting eggs found in the sediments of Lake Constance spanning pre-, peri-, and post-eutrophication. While we expected genome-wide changes, we specifically expected changes in alleles related to pathways involved in mitigating effects of cyanobacterial toxins. We used pairwise F_ST-analyses to identify transcripts that showed strongest divergence among the four different populations and a clustering analysis to identify correlations between allele frequency shifts and changes in abiotic and biotic lake parameters. In a cluster that correlated with the increased abundance of cyanobacteria in Lake Constance we find genes that have been reported earlier to be differentially expressed in response to the cyanobacterial toxin microcystin and to microcystin-free cyanobacteria. We further reveal the enrichment of gene ontology terms that have been shown to be involved in microcystin-related responses in other organisms but not yet in Daphnia and as such are candidate loci for adaptation of natural Daphnia populations to increased cyanobacterial abundances. In conclusion this approach of investigating DNA extracted from Daphnia resting stages allowed to determine frequency changes of loci in a natural population over time.

Insights into the genetic basis of predator-induced response in Daphnia galeata

Phenotypic plastic responses allow organisms to rapidly adjust when facing environmental challenges-these responses comprise morphological, behavioral but also life-history changes. Alteration of life-history traits when exposed to predation risk have been reported often in the ecological and genomic model organism Daphnia. However, the molecular basis of this response is not well understood, especially in the context of fish predation. Here, we characterized the transcriptional profiles of two Daphnia galeata clonal lines with opposed life histories when exposed to fish kairomones. First, we conducted a differential gene expression, identifying a total of 125 candidate transcripts involved in the predator-induced response, uncovering substantial intraspecific variation. Second, we applied a gene coexpression network analysis to find clusters of tightly linked transcripts revealing the functional relations of transcripts underlying the predator-induced response. Our results showed that transcripts involved in remodeling of the cuticle, growth, and digestion correlated with the response to environmental change in D. galeata. Furthermore, we used an orthology-based approach to gain functional information for transcripts lacking gene ontology (GO) information, as well as insights into the evolutionary conservation of transcripts. We could show that our candidate transcripts have orthologs in other Daphnia species but almost none in other arthropods. The unique combination of methods allowed us to identify candidate transcripts, their putative functions, and evolutionary history associated with predator-induced responses in Daphnia. Our study opens up to the question as to whether the same molecular signature is associated with fish kairomones-mediated life-history changes in other Daphnia species.

Partial revision of the neustonic genus Scapholeberis Schoedler, 1858 (Crustacea: Cladocera): decoding of the barcoding results

Water fleas (Crustacea: Cladocera) are among the most intensively studied freshwater invertebrates. However, ecologically important daphniids that live on the surface layer (neuston) remain taxonomically confused. Here we attempt to reconcile genetic and morphological information for the neustonic genus Scapholeberis Schoedler, 1858 (Cladocera: Daphniidae) and present the first revision of the Scapholeberis kingii species group. We analyzed new and existing mitochondrial DNA sequences (сytochrome C oxidase subunit I gene region) together with morphology for all but one of the known species of the neustonic daphniids. Morphological comparisons of available populations, belonging to the Scapholeberis kingii species group from several Australian, Asian and African localities, revealed, that they are almost identical according to parthenogenetic females. However, Australian populations can be reliably distinguished from Asian ones based on the morphology of gamogenetic females. Mitochondrial DNA data analyses revealed divergent lineages (>17% for the DNA barcoding COI region) for the three different species (Australia, Asia and Africa). Based on this set of data, we redescribed S. kingii Sars, 1888 from Australia, its terra typica, and described a new species, S. smirnovi sp.nov. from the Russian Far East, Korea and Japan. The status of populations from Ethiopia and the Republic of South Africa remained unclear, because in the African material and the putative type material, we found only parthenogenetic females. Our results provide an integrative revision of the S. kingii species group and improve the taxonomic scaffold used for barcoding and genomics for the remaining species groups in the daphniid genus Scapholeberis.

Optimized and affordable high-throughput sequencing workflow for preserved and nonpreserved small zooplankton specimens

Genomic analysis of hundreds of individuals is increasingly becoming standard in evolutionary and ecological research. Individual-based sequencing generates large amounts of valuable data from experimental and field studies, while using preserved samples is an invaluable resource for studying biodiversity in remote areas or across time. Yet, small-bodied individuals or specimens from collections are often of limited use for genomic analyses due to a lack of suitable extraction and library preparation protocols for preserved or small amounts of tissues. Currently, high-throughput sequencing in zooplankton is mostly restricted to clonal species, that can be maintained in live cultures to obtain sufficient amounts of tissue, or relies on a whole-genome amplification step that comes with several biases and high costs. Here, we present a workflow for high-throughput sequencing of single small individuals omitting the need for prior whole-genome amplification or live cultures. We establish and demonstrate this method using 27 species of the genus Daphnia, aquatic keystone organisms, and validate it with small-bodied ostracods. Our workflow is applicable to both live and preserved samples at low costs per sample. We first show that a silica-column based DNA extraction method resulted in the highest DNA yields for nonpreserved samples while a precipitation-based technique gave the highest yield for ethanol-preserved samples and provided the longest DNA fragments. We then successfully performed short-read whole genome sequencing from single Daphnia specimens and ostracods. Moreover, we assembled a draft reference genome from a single Daphnia individual (>50× coverage) highlighting the value of the workflow for non-model organisms.

A new species group from the Daphnia curvirostris species complex (Cladocera: Anomopoda) from the eastern Palaearctic: taxonomy, phylogeny and phylogeography

The eastern Palaearctic is a centre of diversity for freshwater cladocerans (Crustacea), but little is known about the evolution and taxonomy of this diversity. Daphnia curvirostris is a Holarctic species complex that has most of its diversity in the eastern Palaearctic. We examined the phylogeography, rates of evolution and taxonomic status for each clade of the D. curvirostris complex using morphological and genetic evidence from four genes. The cybertaxonomical and morphological evidence supported an eastern Palaearctic clade, with at least four species (described here as the Daphnia korovchinskyi sp. nov. group) having diagnostic morphological characters. We also detected convergent morphological characters in the D. curvirostris complex that provided information about species boundaries. Two of the new species (Daphnia koreana sp. nov. and Daphnia ishidai sp. nov.) are known from single ponds and are threatened by human activity. Divergence time estimates suggested an ancient origin (12–28 Mya) for the D. korovchinskyi group, but these estimates are complicated by the small number of calibration points.

Non-toxic and toxic Microcystis aeruginosa reduce the tolerance of Daphnia pulex to low calcium in different degrees: Based on the changes in the key life-history traits

Calcium decline and cyanobacterial blooms pose a serious threat to the crustacean zooplankton Daphnia, which has a high demand for calcium. In the present study, we exposed two different clones of Daphnia pulex to different combinations of calcium concentrations (0.1, 0.5, 1.0, 5.0, 10.0 mg L^-1) and food types (100% Chlorella; 80% Chlorella and 20% non-toxic Microcystis; 80% Chlorella and 20% toxic Microcystis) for 16 days, recorded the key life-history traits, and then used an exponential rise function to fit the traits. Results showed toxic Microcystis and low calcium together negatively affected the survival, development, and reproduction of Daphnia. The negative effect of non-toxic Microcystis and low calcium only affected the development and reproduction. The survival time and reproductive performance increased exponentially with increasing calcium concentration and then approached an asymptotic maximum. Both non-toxic and toxic Microcystis reduced the asymptotic maximum of the reproductive performance. The rising rate at which they reached the asymptotes differed significantly among the three food types; i.e., the reproductive performance of Daphnia was affected in a wider range of calcium concentrations under bad food quality. The findings indicated that Microcystis impaired the tolerance of Daphnia to low calcium, which may cause serious consequences in freshwater ecosystems.

Predator-specific inducible morphological defenses of a water flea against two freshwater predators

The expression of inducible morphological defenses in Daphnia in response to a single predator is a well-known phenomenon. However, predator-specific modifications of the same defensive traits as an adaption to different predator regimes is so far only described for Daphnia barbata. It is unknown if this accounts only for this species or if it is a more widespread, general adaptive response in the genus Daphnia. In the present study, we therefore investigated whether a clone of the pond-dwelling species Daphnia similis responds to different predatory invertebrates (Triops cancriformis; Notonecta maculata) with the expression of predator-specific modifications of the same defensive traits. We showed that Triops-exposed individuals express a significantly longer tail-spine, while body width decreased in comparison to control individuals. Additionally, they also expressed inconspicuous defenses, that is, significantly longer spinules on the dorsal ridge. The Notonecta-exposed D. similis showed a significantly longer tail-spine, longer spinules and a larger spinules bearing area on the dorsal ridge than control individuals as well. However, a geometric morphometric analysis of the head shape revealed significant, predator-specific changes. Triops-exposed individuals expressed a flattened head shape with a pronounced dorsal edge, while Notonecta-exposed individuals developed a high and strongly rounded head. Our study describes so far unrecognized inducible defenses of D. similis against two predators in temporary waters. Furthermore, the predator-dependent change in head shape is in concordance with the 'concept of modality', which highlights the qualitative aspect of natural selection caused by predators.

Unexpected endemism in the Daphnia longispina complex (Crustacea: Cladocera) in Southern Siberia

The biological significance of regional cladoceran morphotypes in the montane regions of the central Palearctic remains poorly understood. In the Holarctic Daphnia longispina complex (Cladocera: Daphniidae), several variants, lineages and species have been proposed as endemic for Southern Siberia. Daphnia turbinata Sars, for example, named after its unusual head shape, is known only from Southern Siberia. Here we sequence DNA of Daphnia from three mitochondrial genes (12S rRNA, 16S rRNA, and NADH dehydrogenase subunit 2, ND2) from 57 localities in Russia and Mongolia (the majority being from Southern Siberia) and place them in evolutionary context with existing data. Our aim was to examine regional endemism of the Daphnia longispina complex in Southern Siberian; to improve the phylogenetic understanding with improved taxonomic and regional sampling, and to better understand the influence of Pleistocene glaciation on the biogeography of these lineages. At least three lineages showed genetic evidence for endemism in Southern Siberia. There was strong support for D. turbinata as a sister lineage to to D. longispina/D. dentifera. Another endemic, Siberian D. cf. longispina, is a sister group to the longispina group in general. Within D. longispina s. str. there was an endemic Siberian clade with a western range boundary near the Yenisei River Basin. Gene flow estimates among populations (based on F_ST values) were very low for clades of D. longispina on a regional (the original 12S dataset), and on a pan-Eurasian (the extended 12S dataset) scale. Negative values of Fu’s F_S and Tajima’s D tests prevailed for the species examined with significant values found for two D. longispina clades, D. dentifera, D. galeata and D. cristata. Our results support the notion that Southern Siberia is an important biogeographic region for cladocerans as it contained unexpected diversity of endemics (such as D. turbinata, D. cf. longispina and lineages of D. umbra and D. longsipina s.str.) and from being the geographic meeting place of expanding postglacial lineages from eastern and western refugia.

Shedding light on the effects of hydraulic fracturing flowback and produced water on phototactic behavior in Daphnia magna

The effluent produced during hydraulic fracturing (i.e. flowback and produced water; FPW), is a complex hyper-saline solution that is known to negatively impact the survival and the fitness of the water flea Daphnia magna, but to date effects on behavior are unstudied. In the current study, the effects of FPW on phototactic behavior of D. magna were examined. Exposure of naïve animals to FPW resulted in a dose-dependent increase in the speed of appearance of daphnids in the illuminated zone of the test apparatus (i.e. a faster positive phototaxis response). A similar dose-dependent response was observed in a test solution where the salt content of FPW was recreated in the absence of other components, suggesting that the effect was largely driven by salinity. The effect of FPW was significant when the raw FPW sample was diluted to 20% of its initial strength, while the effect of salt-matched solution was significant at a 10% dilution. A distinct effect was observed following FPW pre-exposure. After a 24 h pre-exposure to 1.5% FPW, Daphnia displayed a significantly inhibited positive phototaxis response when examined in control water, relative to control animals that were not pre-exposed to FPW. This effect was not observed in salinity pre-exposed animals, however these daphnids displayed a significantly reduced phototactic response when tested in saline waters, indicating a loss of the positive phototaxis seen in naïve organisms. These data indicate that FPW can induce perturbations in the behavior of aquatic invertebrates, an effect that may influence processes such as feeding and predation rates.

New lineages and old species: Lineage diversity and regional distribution of Moina (Crustacea: Cladocera) in China

The distribution and genetic diversity of freshwater zooplankton is understudied in the Eastern Palearctic. Here, we explored the lineage diversity and regional distribution of the genus Moina in China. Members of this genus are often keystone components of freshwater ecosystems and have been frequently subjected to toxicological and physiological studies. Four species of Moina were identified, based on morphology, in 50 of 113 Chinese water bodies examined, and their phylogenetic position was analyzed using both a mitochondrial (mitochondrial cytochrome c oxidase subunit I; COI) and a nuclear marker (the nuclear internal transcribed spacer; ITS-1). Both molecular markers identified four clades corresponding broadly to the morphological species. Mitochondrial DNA analysis showed the presence of four species complexes with eleven lineages across China, five of which were new. However, some lineages (and even individual haplotypes) were widespread in Eurasia, suggesting an ability to disperse over long distances. In contrast, a few lineages exhibited restricted distributions. The nuclear phylogeny also recognized four species of Moina within China and seven very distinct clades. Interestingly, one specimen possessing Moina cf. micrura mtDNA had ITS-1 alleles of the M. cf. brachiata clade. This discordance between mtDNA and nuclear ITS-1 phylogenies is indicative of interspecific introgression and hybridization. Additionally, our COI phylogeny showed apparent paraphyly in two Moina species groups, suggesting introgression of their mitochondrial genomes. Our data shows the regional distribution/diversity of the Moina species complex in a Eurasian context.

Contrasting phylogeographic patterns and demographic history in closely related species of Daphnia longispina group (Crustacea: Cladocera) with focus on North-Eastern Eurasia

Species with large geographic distributions present a challenge for phylogeographic studies due to the logistic difficulties of obtaining adequate samples. Daphnia O.F. Müller (Anomopoda: Daphniidae) is a model genus for evolutionary biology and ecology, but many regions such as the remote areas of Siberia, remain poorly studied. Here we examined genetic polymorphism in the ribosomal 12S and the protein-coding ND2 mitochondrial genes of three closely related taxa of the Daphnia (Daphnia) longispina complex, namely D. galeata Sars, D. longispina O.F. Müller and D. dentifera Forbes. We estimated the phylogenetic relationships among these taxa based on a concatenated alignment of these two genes. Using sequences from the present study and those available in GenBank, we investigated the geographic distributions of the mitochondrial haplotypes of these species and proposed an evolutionary scenario for each taxon. Network structures, haplotype distribution patterns, and FST values indicated significant differences in the evolutionary history of the examined species. Our analysis of D. galeata populations confirmed its recent and fast expansion, without a previous phase of a strong population disconnection. In contrast, the high haplotype diversity in D. dentifera and D. longispina could be explained by the survival of different phylogroups in several glacial refugia located in different geographic regions. For all studied species, maximum haplotype diversity was recorded in the remote regions of Siberia-lakes of the Yenisei River and Transbaikalia. Our study is an important step in our understanding of the evolutionary history of the Daphnia longispina group and provides further evidence of the biogeographic significance of Siberia for freshwater taxa.

Identification of General Patterns of Sex-Biased Expression in Daphnia, a Genus with Environmental Sex Determination

Daphnia reproduce by cyclic-parthenogenesis, where phases of asexual reproduction are intermitted by sexual production of diapause stages. This life cycle, together with environmental sex determination, allow the comparison of gene expression between genetically identical males and females. We investigated gene expression differences between males and females in four genotypes of Daphnia magna and compared the results with published data on sex-biased gene expression in two other Daphnia species, each representing one of the major phylogenetic clades within the genus. We found that 42% of all annotated genes showed sex-biased expression in D. magna This proportion is similar both to estimates from other Daphnia species as well as from species with genetic sex determination, suggesting that sex-biased expression is not reduced under environmental sex determination. Among 7453 single copy, one-to-one orthologs in the three Daphnia species, 707 consistently showed sex-biased expression and 675 were biased in the same direction in all three species. Hence these genes represent a core-set of genes with consistent sex-differential expression in the genus. A functional analysis identified that several of them are involved in known sex determination pathways. Moreover, 75% were overexpressed in females rather than males, a pattern that appears to be a general feature of sex-biased gene expression in Daphnia.

Phylogeography of Daphnia magna Straus (Crustacea: Cladocera) in Northern Eurasia: Evidence for a deep longitudinal split between mitochondrial lineages

Species with a large geographic distributions present a challenge for phylogeographic studies due to logistic difficulties of obtaining adequate sampling. For instance, in most species with a Holarctic distribution, the majority of studies has concentrated on the European or North American part of the distribution, with the Eastern Palearctic region being notably understudied. Here, we study the phylogeography of the freshwater cladoceran Daphnia magna Straus, 1820 (Crustacea: Cladocera), based on partial mitochondrial COI sequences and using specimens from populations spread longitudinally from westernmost Europe to easternmost Asia, with many samples from previously strongly understudied regions in Siberia and Eastern Asia. The results confirm the previously suspected deep split between Eastern and Western mitochondrial haplotype super-clades. We find a narrow contact zone between these two super-clades in the eastern part of Western Siberia, with proven co-occurrence in a single lake in the Novosibirsk region. However, at present there is no evidence suggesting that the two mitochondrial super-clades represent cryptic species. Rather, they may be explained by secondary contact after expansion from different refugia. Interestingly, Central Siberia has previously been found to be an important contact zone also in other cladoceran species, and may thus be a crucial area for understanding the Eurasian phylogeography of freshwater invertebrates. Together, our study provides an unprecedented complete, while still not global, picture of the phylogeography of this important model species.

Divergent clades or cryptic species? Mito-nuclear discordance in a Daphnia species complex

BACKGROUND

Genetically divergent cryptic species are frequently detected by molecular methods. These discoveries are often a byproduct of molecular barcoding studies in which fragments of a selected marker are used for species identification. Highly divergent mitochondrial lineages and putative cryptic species are even detected in intensively studied animal taxa, such as the crustacean genus Daphnia. Recently, eleven such lineages, exhibiting genetic distances comparable to levels observed among well-defined species, were recorded in the D. longispina species complex, a group that contains several key taxa of freshwater ecosystems. We tested if three of those lineages represent indeed distinct species, by analyzing patterns of variation of ten nuclear microsatellite markers in six populations.

RESULTS

We observed a discordant pattern between mitochondrial and nuclear DNA, as all individuals carrying one of the divergent mitochondrial lineages grouped at the nuclear level with widespread, well-recognized species coexisting at the same localities (Daphnia galeata, D. longispina, and D. cucullata).

CONCLUSIONS

A likely explanation for this pattern is the introgression of the mitochondrial genome of undescribed taxa into the common species, either in the distant past or after long-distance dispersal. The occurrence of highly divergent but rare mtDNA lineages in the gene pool of widespread species would suggest that hybridization and introgression in the D. longispina species complex is frequent even across large phylogenetic distances, and that discoveries of such distinct clades must be interpreted with caution. However, maintenance of ancient polymorphisms through selection is another plausible alternative that may cause the observed discordance and cannot be entirely excluded.

Phylogeography of the Chydorus sphaericus Group (Cladocera: Chydoridae) in the Northern Palearctic

The biodiversity and the biogeography are still poorly understood for freshwater invertebrates. The crustacean Chydorus sphaericus-brevilabris complex (Cladocera: Chydoridae) is composed of species that are important components of Holarctic freshwater food webs. Recent morphological and genetic study of the complex has indicated a substantial species diversity in the northern hemisphere. However, we know little of the geographic boundaries of these novel lineages. Moreover, a large section of the Palearctic remains unexamined at the genetic level. Here we attempt to address the biodiversity knowledge gap for the Chydorus sphaericus group in the central Palearctic and assess its diversity and biogeographic boundaries. We sequenced nuclear (ITS-2) and mitochondrial (COI) gene regions of Chydorus specimens across the Palearctic and compared them with already available Holarctic sequences. We detected six main clades in the C. sphaericus group in the Palearctic, of which two of the groups are novel. Three of the more divergent clades are geographically widespread. The central portion of Eurasia (the Yenisey River basin) appears to be a narrow zone of secondary contact for phylogroups that expanded from European and Beringian refugia. As such, the previously unsampled central Palearctic represents an important region for understanding the evolutionary consequences of Pleistocene climatic oscillations on the Chydorus sphaericus group.

DNA Barcoding Reveals High Cryptic Diversity in the North Eurasian Moina Species (Crustacea: Cladocera)

Species of the genus Moina Baird (Cladocera: Moinidae) often dominate freshwater crustacean communities in temporary water bodies. Several species of Moina are used as food for fish larvae in aquaculture, as bioindicators in toxicological studies, and as common subjects for physiological studies. The aim of this paper is to estimate biodiversity of Moina in northern Eurasia using the standard DNA barcoding approach based on the cytochrome c oxidase subunit I (COI) gene. We analysed 160 newly obtained and 157 existing COI sequences, and found evidence for 21 phylogroups of Moina, some of which were detected here for the first time. Our study confirmed the opinion that the actual species diversity of cladocerans is several times higher than is presently accepted. Our results also indicated that Moina has the second richest species diversity among the cladoceran genera (with only Daphnia O. F. Mueller having a greater diversity of species). Our study strongly supports division of Moina into two faunistic groups: European-Western Siberian and Eastern Siberian-Far Eastern, with a transitional zone at the Yenisey River basin (Eastern Siberia). Here, we refrain from taxonomic descriptions of new species, as this requires a thorough morphological and taxonomic study for each putative taxon.

Do saline taxa evolve faster? Comparing relative rates of molecular evolution between freshwater and marine eukaryotes

The major branches of life diversified in the marine realm, and numerous taxa have since transitioned between marine and freshwaters. Previous studies have demonstrated higher rates of molecular evolution in crustaceans inhabiting continental saline habitats as compared with freshwaters, but it is unclear whether this trend is pervasive or whether it applies to the marine environment. We employ the phylogenetic comparative method to investigate relative molecular evolutionary rates between 148 pairs of marine or continental saline versus freshwater lineages representing disparate eukaryote groups, including bony fish, elasmobranchs, cetaceans, crustaceans, mollusks, annelids, algae, and other eukaryotes, using available protein-coding and noncoding genes. Overall, we observed no consistent pattern in nucleotide substitution rates linked to habitat across all genes and taxa. However, we observed some trends of higher evolutionary rates within protein-coding genes in freshwater taxa-the comparisons mainly involving bony fish-compared with their marine relatives. The results suggest no systematic differences in substitution rate between marine and freshwater organisms.

Daphnia magna transcriptome by RNA-Seq across 12 environmental stressors

The full exploration of gene-environment interactions requires model organisms with well-characterized ecological interactions in their natural environment, manipulability in the laboratory and genomic tools. The waterflea Daphnia magna is an established ecological and toxicological model species, central to the food webs of freshwater lentic habitats and sentinel for water quality. Its tractability and cyclic parthenogenetic life-cycle are ideal to investigate links between genes and the environment. Capitalizing on this unique model system, the STRESSFLEA consortium generated a comprehensive RNA-Seq data set by exposing two inbred genotypes of D. magna and a recombinant cross of these genotypes to a range of environmental perturbations. Gene models were constructed from the transcriptome data and mapped onto the draft genome of D. magna using EvidentialGene. The transcriptome data generated here, together with the available draft genome sequence of D. magna and a high-density genetic map will be a key asset for future investigations in environmental genomics.

Genetic Diversity of Daphnia pulex in the Middle and Lower Reaches of the Yangtze River

Increased human activities and environmental changes may lead to genetic diversity variations of Cladocerans in water. Daphnia pulex are distributed throughout the world and often regarded as a model organism. The 16S rDNA, cytochrome c oxidase subunit I (COI), and 18S genes were used as molecular marks. The genetic diversity and phylogeny of D. pulex obtained from 10 water bodies in the middle and lower reaches of the Yangtze River were studied. For 16S rDNA, COI gene, and 18S gene, the A+T content (65.4%, 58.4%, and 54.6%) was significantly higher than the G+C content (34.6%, 41.6% and 45.4%). This result was consistent with higher A and T contents among invertebrates. Based on the genetic distances of 16S rDNA and COI genes, the genetic differences of D. pulex from 10 water bodies located in the middle and lower reaches of the Yangtze River in China was minimal (0%-0.8% for 16S rDNA and 0%-1.5% for COI gene). However, D. pulex evolved into two branches in the phylogenetic trees, which coincided with its geographical distribution. Compared with D. pulex from other countries, the average genetic distance of D. pulex obtained from 10 water bodies in the middle and lower reaches of the Yangtze River reached 9.1%-10.5%, thereby indicating that D. pulex may have evolved into different subspecies.

Molecular diversity of some species belonging to the genus Daphnia O. F. Müller, 1785 (Crustacea: Cladocera) in Turkey

Daphnia is a freshwater zooplankton species with controversial taxonomy due to its high morphological variation linked to environmental factors and inter-specific hybridization and polyploidy in some groups. The aim of the present study is to examine molecular diversity of some Daphnia species in Turkey and to establish DNA barcodes of Turkish Daphnia species. Sequence analysis was performed using 540 bp region of cytochrome oxidase subunit I gene of mitochondrial DNA. A total of 34 haplotypes have been identified for Turkey. Daphnia pulex complex was divided into two clades with 16.1% sequence divergence according to molecular taxonomy based on Kimura 2-parameter. The clade which was molecularly diverged from Daphnia pulex with 16.1% sequence divergence was found to show 99% similarity with Daphnia cf. pulicaria (sensu Alonso 1996) instead of Daphnia pulicaria Forbes, 1893. Furthermore, this study has contributed to Turkish zoogeography by demonstrating the distribution of Daphnia species in Turkey.

Genetic structure of Daphnia galeata populations in Eastern China

This study presents the first examination of the genetic structure of Daphnia longispina complex populations in Eastern China. Only one species, D. galeata, was present across the eight investigated lakes; as identified by taxon assignment using allelic variation at 15 microsatellite loci. Three genetically differentiated D. galeata subgroups emerged independent of the type of statistical analysis applied. Thus, Bayesian clustering, discriminant analysis based on results from factorial correspondence analysis, and UPGMA clustering consistently showed that populations from two neighbouring lakes were genetically separated from a mixture of genotypes found in other lakes, which formed another two subgroups. Clonal diversity was high in all D. galeata populations, and most samples showed no deviation from Hardy-Weinberg equilibrium, indicating that clonal selection had little effect on the genetic diversity. Overall, populations did not cluster by geographical origin. Further studies will show if the observed pattern can be explained by natural colonization processes or by recent anthropogenic impact on predominantly artificial lakes.

Confirmation through genetic analysis of the existence of many local phyloclades of the genus Simocephalus (Crustacea, Cladocera) in China

Previously, a series of Simocephalus taxa (Cladocera: Daphniidae) from China were described. Most were proposed to be junior synonyms in the last revision of the genus. Using original material from China and data from GenBank, we investigate the biodiversity and phylogeny of Simocephalus using sequences of the cytochrome c oxidase subunit I (COI) and the nuclear 18S genes. In both cases, neighbor-joining, maximum likelihood and Bayesian inference analyses led to highly congruent tree topologies. The grouping of the deeper clades agrees with the inter-generic classification of Orlova-Bienkowskaja (2001). Only the populations of S. serrulatus from Eurasia and North America seem to be closely related, and there are no other shared species between the two continents. Our study unambiguously confirms the existence of many lineages from the subgenera of Simocephalus (Echinocaudus) and Simocephalus s.str. in China, but their morphology needs to be reexamined by taking a wider range of characters (e.g., of female thoracic limbs and adult males) into consideration.