Daphnia magna transcriptome by RNA-Seq across 12 environmental stressors

A Novel Method for the Assessment of Feeding Rate as a Phenotypic Endpoint for the Impact of Pollutants in Daphnids

Traditional approaches for monitoring aquatic pollution primarily rely on chemical analysis and the detection of pollutants in the aqueous environments. However, these methods lack realism and mechanistic insight and, thus, are increasingly supported by effect-based methods, which offer sensitive endpoints. In this context, daphnids, a freshwater species used extensively in molecular ecotoxicology, offer fast and noninvasive approaches to assess the impact of pollutants. Among the phenotypic endpoints used, feeding rate is a highly sensitive approach because it provides evidence of physiological alterations even in sublethal concentrations. However, there has been no standardized method for measuring feeding rate in daphnids, and several approaches follow different protocols. There is a diversity among tests employing large volumes, extensive incubation times, and high animal densities, which in turn utilize measurements of algae via fluorescence, radiolabeling, or counting ingested cells. These tests are challenging and laborious and sometimes require cumbersome instrumentation. In the present study, we optimized the conditions of a miniaturized fast, sensitive, and high-throughput assay to assess the feeding rate based on the ingestion of fluorescent microparticles. The protocol was optimized in neonates in relation to the concentration of microplastic and the number of animals to increase reproducibility. Daphnids, following exposures to nonlethal concentrations, were incubated with microplastics; and, as filter feeders, they ingest microparticles. The new approach revealed differences in the physiology of daphnids in concentrations below the toxicity limits for a range of pollutants of different modes of action, thus proving feeding to be a more sensitive and noninvasive endpoint in pollution assessment. Environ Toxicol Chem 2024;00:1-11. © 2024 The Author(s). Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

RNA-seq analysis reveals changes in mRNA expression during development in Daphnia mitsukuri

Temporal transcriptional variation is a major contributor to functional evolution and the developmental process. Parthenogenetic water fleas of the genus Daphnia (Cladocera) provide an ideal model to characterize gene expression patterns across distinct developmental stages. Herein, we report RNA-seq data for female Daphnia mitsukuri at three developmental stages: the embryo, juvenile (three timepoints) and adult. Comparisons of gene expression patterns among these three developmental stages and weighted gene co-expression network analysis based on expression data across developmental stages identified sets of genes underpinning each of the developmental stages of D. mitsukuri. Specifically, highly expressed genes (HEGs) at the embryonic developmental stage were associated with cell proliferation, ensuring the necessary foundation for subsequent development; HEGs at the juvenile stages were associated with chemosensory perception, visual perception and neurotransmission, allowing individuals to enhance detection of potential environmental risks; HEGs at the adult stage were associated with antioxidative defensive systems, enabling adults to mount an efficient response to perceived environmental risks. Additionally, we found a significant overlap between expanded gene families of Daphnia species and HEGs at the juvenile stages, and these genes were associated with visual perception and neurotransmission. Our work provides a resource of developmental transcriptomes, and comparative analyses that characterize gene expression dynamics throughout development of Daphnia.

Dietary exposure of the water flea Daphnia galeata to microcystin-LR

Harmful substances like the cyanotoxin microcystin-leucine-arginine (MC-LR) are commonly found in eutrophic freshwater environments, posing risks to aquatic organisms. The water flea, Daphnia, is a well-established model organism for environmental toxicology research. Nevertheless, there is currently insufficient research on the genes that respond to MC-LR in Daphnia galeata. This study aimed to gain insights into the notable genes that react significantly to MC-LR. In this study, we generated an extensive RNA-Seq sequences isolated from the D. galeata HK strain, Han River in Korea. This strain was nourished with a diet of the green microalga Chlorella vulgaris and treated with pure MC-LR at a concentration of 36 ug/L. The transcriptome profile in response to the MC-LR treatment was obtained and 336 differentially expressed genes were subjected to Gene Ontology (GO) and euKaryotic Orthologous Groups of proteins analyses. GO enrichment analysis showed that chemical stimulus, amino sugar metabolic and catabolic process, oxidative stress, and detoxification were highly enriched, in reverse, proteolysis and fucosylation were underpresented. Detoxification process related genes such as peroxidase-like, chorion, and thyroid peroxidase-like were enriched for eliminating or neutralizing MC_LR from an organism's body. Furthermore, functional protein classification revealed an upregulation of lipid and inorganic ion transport processes, while amino acid and carbohydrate transport processes were found to be downregulated. These findings offer insights into how organisms respond to ecotoxic stimuli, providing valuable information for understanding adaptation or defense pathways.

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.

Responses to predation pressure involve similar sets of genes in two divergent species of Daphnia

Species that are not closely related can express similar inducible traits, but molecular mechanisms underlying the observed responses are often unknown, nor is it known if these mechanisms are shared between such species. Here, we compared transcriptional profiles of two Daphnia species (D. mitsukuri and D. sinensis) from different subgenera, at both juvenile and adult developmental stages. Both species were exposed to the same predation threat (fish kairomones), and both showed similar induced morphological changes (reduced body length). At the early developmental stage, response to predation risk resulted in similar changes in expression levels of 23 orthologues in both species. These orthologues, involved in 107 GO categories, changed in the same direction in both species (over- or underexpressed), in comparison to non-exposed controls. Several of these orthologues were associated with DNA replication, structural constituents of cuticle or innate immune response. In both species, the differentially expressed (DE) genes on average had higher ω (dN /dS ) values than non-DE genes, suggesting that these genes had experienced greater positive selection or lower purifying selection than non-DE genes. Overall, our results suggest that similar suites of genes, responding in similar ways to predation pressure, have been retained in Daphnia for many millions of years.

Characterizing nutritional phenotypes using experimental nutrigenomics: Is there nutrient-specificity to different types of dietary stress?

The ability to directly measure and monitor poor nutrition in individual animals and ecological communities is hampered by methodological limitations. In this study, we use nutrigenomics to identify nutritional biomarkers in a freshwater zooplankter, Daphnia pulex, a ubiquitous primary consumer in lakes and a sentinel of environmental change. We grew animals in six ecologically relevant nutritional treatments: nutrient replete, low carbon (food), low phosphorus, low nitrogen, low calcium and high Cyanobacteria. We extracted RNA for transcriptome sequencing to identify genes that were nutrient responsive and capable of predicting nutritional status with a high degree of accuracy. We selected a list of 125 candidate genes, which were subsequently pruned to 13 predictive potential biomarkers. Using a nearest-neighbour classification algorithm, we demonstrate that these potential biomarkers are capable of classifying our samples into the correct nutritional group with 100% accuracy. The functional annotation of the selected biomarkers revealed some specific nutritional pathways and supported our hypothesis that animal responses to poor nutrition are nutrient specific and not simply different presentations of slow growth or energy limitation. This is a key step in uncovering the causes and consequences of nutritional limitation in animal consumers and their responses to small- and large-scale changes in biogeochemical cycles.

QTL study reveals candidate genes underlying host resistance in a Red Queen model system

Specific interactions of host and parasite genotypes can lead to balancing selection, maintaining genetic diversity within populations. In order to understand the drivers of such specific coevolution, it is necessary to identify the molecular underpinnings of these genotypic interactions. Here, we investigate the genetic basis of resistance in the crustacean host, Daphnia magna, to attachment and subsequent infection by the bacterial parasite, Pasteuria ramosa. We discover a single locus with Mendelian segregation (3:1 ratio) with resistance being dominant, which we call the F locus. We use QTL analysis and fine mapping to localize the F locus to a 28.8-kb region in the host genome, adjacent to a known resistance supergene. We compare the 28.8-kb region in the two QTL parents to identify differences between host genotypes that are resistant versus susceptible to attachment and infection by the parasite. We identify 13 genes in the region, from which we highlight eight biological candidates for the F locus, based on presence/absence polymorphisms and differential gene expression. The top candidates include a fucosyltransferase gene that is only present in one of the two QTL parents, as well as several Cladoceran-specific genes belonging to a large family that is represented in multiple locations of the host genome. Fucosyltransferases have been linked to resistance in previous studies of Daphnia-Pasteuria and other host-parasite systems, suggesting that P. ramosa spore attachment could be mediated by changes in glycan structures on D. magna cuticle proteins. The Cladoceran-specific candidate genes suggest a resistance strategy that relies on gene duplication. Our results add a new locus to a growing genetic model of resistance in the D. magna-P. ramosa system. The identified candidate genes will be used in future functional genetic studies, with the ultimate aim to test for cycles of allele frequencies in natural populations.

High-quality genome of Diaphanosoma dubium provides insights into molecular basis of its broad ecological adaptation

Diaphanosoma dubium Manuilova, 1964, is a widespread planktonic water flea in Asian freshwater. Although sharing similar ecological roles with species of Daphnia, studies on D. dubium and its congeners are still few and lacking a genome for the further studies. Here, we assembled a high quality and chromosome level genome of D. dubium by combining long reads sequencing and Hi-C technologies. The total length of assembled genome was 101.8 Mb, with 98.92 Mb (97.2%) anchored into 22 chromosomes. Through comparative genomic analysis, we found the genes, involved in anti-ROS, detoxification, protein digestion, germ cells regulation and protection, underwent expansion in D. dubium. These genes and their expansion helpfully explain its widespread geographical distribution and dominance in eutrophic waters. This study provides insight into the adaptive evolution of D. dubium at genomic perspectives, and the present high quality genomic resource will be a footstone for future omics studies of the species and its congeners.

Daphnia as a model organism to probe biological responses to nanomaterials-from individual to population effects via adverse outcome pathways

The importance of the cladoceran Daphnia as a model organism for ecotoxicity testing has been well-established since the 1980s. Daphnia have been increasingly used in standardised testing of chemicals as they are well characterised and show sensitivity to pollutants, making them an essential indicator species for environmental stress. The mapping of the genomes of D. pulex in 2012 and D. magna in 2017 further consolidated their utility for ecotoxicity testing, including demonstrating the responsiveness of the Daphnia genome to environmental stressors. The short lifecycle and parthenogenetic reproduction make Daphnia useful for assessment of developmental toxicity and adaption to stress. The emergence of nanomaterials (NMs) and their safety assessment has introduced some challenges to the use of standard toxicity tests which were developed for soluble chemicals. NMs have enormous reactive surface areas resulting in dynamic interactions with dissolved organic carbon, proteins and other biomolecules in their surroundings leading to a myriad of physical, chemical, biological, and macromolecular transformations of the NMs and thus changes in their bioavailability to, and impacts on, daphnids. However, NM safety assessments are also driving innovations in our approaches to toxicity testing, for both chemicals and other emerging contaminants such as microplastics (MPs). These advances include establishing more realistic environmental exposures via medium composition tuning including pre-conditioning by the organisms to provide relevant biomolecules as background, development of microfluidics approaches to mimic environmental flow conditions typical in streams, utilisation of field daphnids cultured in the lab to assess adaption and impacts of pre-exposure to pollution gradients, and of course development of mechanistic insights to connect the first encounter with NMs or MPs to an adverse outcome, via the key events in an adverse outcome pathway. Insights into these developments are presented below to inspire further advances and utilisation of these important organisms as part of an overall environmental risk assessment of NMs and MPs impacts, including in mixture exposure scenarios.

Multigenerational DNA methylation responses to copper exposure in Daphnia: Potential targets for epigenetic biomarkers?

Epigenetic mechanisms are moving to the forefront of environmental sciences, as environmentally induced epigenetic changes shape biological responses to chemical contamination. This work focused on Daphnia as a representative of potentially threatened freshwater biota, aiming to gain an insight into the involvement of epigenetic mechanisms in their response and eventual adaptation to metal contamination. Copper-induced DNA methylation changes, their potential transgenerational inheritance, and life-history traits were assessed. Organisms with different histories of past exposure to copper were exposed to toxic levels of the element for one generation (F0) and then monitored for three subsequent unexposed generations (F1, F2, and F3). Overall, methylation changes targeted important genes for counteracting the effects of metals and oxidative stress, including dynein light chain, ribosomal kinase and nuclear fragile X mental retardation-interacting protein. Also, contrasting overall and gene-specific methylation responses were observed in organisms differing in their history of exposure to copper, with different transgenerational methylation responses being also identified among the two groups, without apparent life-history costs. Taken together, these results demonstrate the capacity of copper to promote epigenetic transgenerational inheritance in a manner related explicitly to history of exposure, thereby supporting the development and incorporation of epigenetic biomarkers in risk assessment frameworks.

Negative Effects of Cyanotoxins and Adaptative Responses of Daphnia

The plethora of cyanobacterial toxins are an enormous threat to whole ecosystems and humans. Due to eutrophication and increases in lake temperatures from global warming, changes in the distribution of cyanobacterial toxins and selection of few highly toxic species/ strains are likely. Globally, one of the most important grazers that controls cyanobacterial blooms is Daphnia, a freshwater model organism in ecology and (eco)toxicology. Daphnia-cyanobacteria interactions have been studied extensively, often focusing on the interference of filamentous cyanobacteria with Daphnia's filtering apparatus, or on different nutritional constraints (the lack of essential amino acids or lipids) and grazer toxicity. For a long time, this toxicity only referred to microcystins. Currently, the focus shifts toward other deleterious cyanotoxins. Still, less than 10% of the total scientific output deals with cyanotoxins that are not microcystins; although these other cyanotoxins can occur just as frequently and at similar concentrations as microcystins in surface water. This review discusses the effects of different cyanobacterial toxins (hepatotoxins, digestive inhibitors, neurotoxins, and cytotoxins) on Daphnia and provides an elaborate and up-to-date overview of specific responses and adaptations of Daphnia. Furthermore, scenarios of what we can expect for the future of Daphnia-cyanobacteria interactions are described by comprising anthropogenic threats that might further increase toxin stress in Daphnia.

Stochastic variation in gut bacterial community affects reproductive rates in the water flea Daphnia magna

It is well-documented that perturbation of the gut bacterial community can influence the reproductive rates of the host. Less is known about how natural ecological processes can change the bacterial composition in the gut and how such changes influence the reproductive rate of the host. Here, we provide novel experimental insights into such processes using the clonally reproducing water flea, Daphnia magna. A total of 20 replicate cultures were reared for 5 weeks (Phase 1) to allow for divergence of bacterial communities through stochastic processes (i.e. drift, founder effects, and/or colonization). Duplicate cultures created from each of these were reared for 21 days (Phase 2) while recording reproductive rates. There was a significant repeatability in reproductive rates between these duplicates, suggesting that divergence of the bacterial communities during Phase 1 translated into reproductive rate effects during Phase 2. This was further supported by significant differences in the relative abundance of gut bacteria (investigated by amplicon sequencing of a part of the 16S rRNA gene) between cultures with high and low reproductive rate in Phase 2. These results are consistent with the hypothesis that stochastic processes can cause natural variation in the bacterial composition in the gut, which in turn affect host reproductive rates.

Thermal stress and mutation accumulation increase heat shock protein expression in Daphnia

Understanding the short- and long-term consequences of climate change is a major challenge in biology. For aquatic organisms, temperature changes and drought can lead to thermal stress and habitat loss, both of which can ultimately lead to higher mutation rates. Here, we examine the effect of high temperature and mutation accumulation on gene expression at two loci from the heat shock protein (HSP) gene family, HSP60 and HSP90. HSPs have been posited to serve as 'mutational capacitors' given their role as molecular chaperones involved in protein folding and degradation, thus buffering against a wide range of cellular stress and destabilization. We assayed changes in HSP expression across 5 genotypes of Daphnia magna, a sentinel species in ecology and environmental biology, with and without acute exposure to thermal stress and accumulated mutations. Across genotypes, HSP expression increased ~ 6× in response to heat and ~ 4× with mutation accumulation, individually. Both factors simultaneously (lineages with high mutation loads exposed to high heat) increased gene expression ~ 23×-much more than that predicted by an additive model. Our results corroborate suggestions that HSPs can buffer against not only the effects of heat, but also mutations-a combination of factors both likely to increase in a warming world.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10682-022-10209-1.

Using the Daphnia magna Transcriptome to Distinguish Water Source: Wetland and Stormwater Case Studies

A major challenge in ecotoxicology is accurately and sufficiently measuring chemical exposures and biological effects given the presence of complex and dynamic contaminant mixtures in surface waters. It is impractical to quantify all chemicals in such matrices over space and time, and even if it were practical, concomitant biological effects would not be elucidated. Our study examined the performance of the Daphnia magna transcriptome to detect distinct responses across three water sources in Minnesota: laboratory (well) waters, wetland waters, and storm waters. Pyriproxyfen was included as a gene expression and male neonate production positive control to examine whether gene expression resulting from exposure to this well-studied juvenoid hormone analog can be detected in complex matrices. Laboratory-reared (<24 h) D. magna were exposed to a water source and/or pyriproxyfen for 16 days to monitor phenotypic changes or 96 h to examine gene expression responses using Illumina HiSeq 2500 (10 million reads per library, 50-bp paired end [2 × 50]). The results indicated that a unique gene expression profile was produced for each water source. At 119 ng/L pyriproxyfen (~25% effect concentration) for male neonate production, as expected, the Doublesex1 gene was up-regulated. In descending order, gene expression patterns were most discernable with respect to pyriproxyfen exposure status, season of stormwater sample collection, and wetland quality, as indicated by the index of biological integrity. However, the biological implications of the affected genes were not broadly clear given limited genome resources for invertebrates. Our study provides support for the utility of short-term whole-organism transcriptomic testing in D. magna to discern sample type, but highlights the need for further work on invertebrate genomics. Environ Toxicol Chem 2022;41:2107-2123. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Rapid communication: effects of cadmium exposure on the growth-related genes of Daphnia magna

Few data are available regarding the effects of gene expression on growth in Daphnia magna. The aim of this study was to determine the influence of cadmium (Cd) exposure on global gene transcription and growth-related genes in D. magna using RNASeq generated data. Our results demonstrated that Cd exposure decreased gene expression, but did not adversely affect the expression of growth-related genes, suggesting differential allocation of resources to growth avoids the deleterious effect of the toxicant on this trait.

Coping with antagonistic predation risks: Predator-dependent unique responses are dominant in Ceriodaphnia cornuta

Inducible defences of prey are evolved under diverse and variable predation risks. However, during the co-evolution of prey and multiple predators, the responses of prey to antagonistic predation risks, which may put the prey into a dilemma of responding to predators, remain unclear. Based on antagonistic predation pressure from an invertebrate (Chaoborus larvae) and a vertebrate (Rhodeus ocellatus) predator, we studied the responses of multiple traits and transcriptomes of the freshwater crustacean Ceriodaphnia cornuta under multiple predation risks. Chaoborus predation risk altered the expression of genes encoding cuticle proteins and modulated the biosynthesis of steroid hormones, cutin, suberine, and wax, leading to the development of horns and increase in size at the late developmental stage. Meanwhile, fish predation risk primarily triggered genes encoding ribosomes and those involved in unsaturated fatty acid biosynthesis and cysteine and methionine metabolism, resulting in smaller individual size and earlier reproduction. Inducible responses of both transcriptome and individual traits revealed that predator-dependent unique responses were dominant and the dilemma of antagonistic responses was relatively limited. However, the unique individual traits in response to invertebrate predation could be significantly impaired by vertebrate predation risk, even though the unique responses to different predators were extremely weakly correlated and could be elicited simultaneously. These results indicate that diverse predator-dependent unique responses are favoured by Ceriodaphnia during its co-evolution with multiple predators. Nonetheless, Ceriodaphnia is not a generalist that can fully adopt all predator-dependent unique responses simultaneously under multiple predation risks.

Degradation of Sulfoxaflor in Water and Soil: Kinetics, Degradation Pathways, Transformation Product Identification, and Toxicity

Hydrolysis, photolysis, and soil degradation are important degradation pathways of pesticides and might generate toxic chemicals and pose threats to the environment. Sulfoxaflor is a widely used neonicotinoid pesticide, but few studies have been conducted to research its environmental behaviors and residues. Herein, the hydrolysis and photolysis of sulfoxaflor in water and degradation in four typical Chinese soils were systematically studied. In addition, degradation products, pathways, and toxicity to Daphnia magna were also investigated. Sulfoxaflor can undergo photolysis and soil degradation, and the degradation percentage was greater than 90% after 96 h or 96 days, respectively. However, sulfoxaflor was not degraded or only slightly degraded during in hydrolysis and was not photodegraded in acidic water or sterilized soil. Four degradation products were screened by UHPLC-Q-Orbitrap-HRMS, three candidates (X11719474, X11721061, and X11718922) were synthesized, and the photolysis and soil degradation kinetics were explored. The possible pathways were elucidated. Sulfoxaflor, X11718922, and X11721061 had a low toxicity, and X11719474 (48 h EC₅₀ 0.74 mg/L) had a high toxicity to Daphnia magna. Thus, sulfoxaflor and its degradation products could induce tissue damage in Daphnia magna. This work offers a theoretical basis for the application and ecological risk assessment of sulfoxaflor.

Population-genomic analysis identifies a low rate of global adaptive fixation in the proteins of the cyclical parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared with Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction.

Draft genomes for one Microcystis-resistant and one Microcystis-sensitive strain of the water flea, Daphnia pulicaria

Daphnia species are well-suited for studying local adaptation and evolutionary responses to stress(ors) including those caused by algal blooms. Algal blooms, characterized by an overgrowth (bloom) of cyanobacteria, are detrimental to the health of aquatic and terrestrial members of freshwater ecosystems. Some strains of Daphnia pulicaria have demonstrated resistance to toxic algae and the ability to mitigate toxic algal blooms. Understanding the genetic mechanism associated with this toxin resistance requires adequate genomic resources. Using whole-genome sequence data mapped to the Daphnia pulex reference genome (PA42), we present reference-guided draft assemblies from one tolerant and one sensitive strain of D. pulicaria, Wintergreen-6 (WI-6), and Bassett-411 (BA-411), respectively. Assessment of the draft assemblies reveal low contamination levels, and high levels (95%) of genic content. Reference scaffolds had coverage breadths of 98.9-99.4%, and average depths of 33X and 29X for BA-411 and WI-6, respectively. Within, we discuss caveats and suggestions for improving these draft assemblies. These genomic resources are presented with a goal of contributing to the resources necessary to understand the genetic mechanisms and associations of toxic prey resistance observed in this species.

Pharmacological Modulation of Behaviour, Serotonin and Dopamine Levels in Daphnia magna Exposed to the Monoamine Oxidase Inhibitor Deprenyl

This study assessed the effects of the monoamine oxidase (MAO) inhibitor deprenyl in Daphnia magna locomotor activity. The mechanisms of action of deprenyl were also determined by studying the relationship between behaviour, MAO activity and neurotransmitter levels. Modulation of the D. magna monoamine system was accomplished by 24 h exposure to two model psychotropic pharmaceuticals with antagonistic and agonistic serotonin signalling properties: 10 mg/L of 4-chloro-DL-phenylalanine (PCPA) and 1 mg/L of deprenyl, respectively. Contrasting behavioural outcomes were observed for deprenyl and PCPA reflected in decreased basal locomotor activity and enhanced habituation for the former compound and delayed habituation for the latter one. Deprenyl exposure inhibited monoamine oxidase (MAO) activity and increased the concentrations of serotonin, dopamine and the dopamine metabolite 3-methoxytyramine in whole D. magna extracts. Our findings indicate that D. magna is a sensitive and useful nonvertebrate model for assessing the effects of short-term exposure to chemicals that alter monoamine signalling changes.

Transcriptomic response of Daphnia magna to nitrogen- or phosphorus-limited diet

Effects of nutrient-imbalanced diet on the growth and fitness of zooplankton were widely reported as key issues to aquatic ecology. However, little is known about the molecular mechanisms driving the physiological changes of zooplankton under nutrient stress.In this study, we investigated the physiological fitness and transcriptomic response of Daphnia magna when exposed to nitrogen (N)-limited or phosphorus (P)-limited algal diet (Chlamydomonas reinhardtii) compared to regular algae (N and P saturated).D. magna showed higher ingestion rates and overexpression of genes encoding digestive enzymes when fed with either N-limited or P-limited algae, reflecting the compensatory feeding. Under P-limitation, both growth rate and reproduction rate of D. magna were greatly reduced, which could be attributed to the downregulated genes within the pathways of cell cycle and DNA replication. Growth rate of D. magna under N-limitation was similar to normal group, which could be explained by the high methylation level (by degradation of methionine) supporting the body development.Phenotypic changes of D. magna under nutrient stress were explained by gene and pathway regulations from transcriptome data. Generally, D. magna invested more on growth under N-limitation but kept maintenance (e.g., cell structure and defense to external stress) in priority under P-limitation. Post-translational modifications (e.g., methylation and protein folding) were important for D. magna to deal with nutrient constrains.This study reveals the fundamental mechanisms of zooplankton in dealing with elemental imbalanced diet and sheds light on the transfer of energy and nutrient in aquatic ecosystems.

Extensive standing genetic variation from a small number of founders enables rapid adaptation in Daphnia

We lack a thorough understanding of the origin and maintenance of standing genetic variation that enables rapid evolutionary responses of natural populations. Whole genome sequencing of a resurrected Daphnia population shows that standing genetic variation in over 500 genes follows an evolutionary trajectory that parallels the pronounced and rapid adaptive evolution of multiple traits in response to predator-driven natural selection and its subsequent relaxation. Genetic variation carried by only five founding individuals from the regional genotype pool is shown to suffice at enabling the observed evolution. Our results provide insight on how natural populations can acquire the genomic variation, through colonization by a few regional genotypes, that fuels rapid evolution in response to strong selection pressures. While these evolutionary responses in our study population involved hundreds of genes, we observed no evidence of genetic erosion.

New insights on the effects of ionic liquid structural changes at the gene expression level: Molecular mechanisms of toxicity in Daphnia magna

Knowledge on the molecular basis of ionic liquids' (ILs) ecotoxicity is critical for the development of these designer solvents as their structure can be engineered to simultaneously meet functionality performance and environmental safety. The molecular effects of ILs were investigated by using RNA-sequencing following Daphnia magna exposure to imidazolium- and cholinium-based ILs: 1-ethyl-3-methylimidazolium chloride ([C₂mim]Cl), 1-dodecyl-3-methylimidazolium chloride ([C₁₂mim]Cl) and cholinium chloride ([Chol]Cl)-; the selection allowing to compare different families and cation alkyl chains. ILs shared mechanisms of toxicity focusing e.g. cellular membrane and cytoskeleton, oxidative stress, energy production, protein biosynthesis, DNA damage, disease initiation. [C₂mim]Cl and [C₁₂mim]Cl were the least and the most toxic ILs at the transcriptional level, denoting the role of the alkyl chain as a driver of ILs toxicity. Also, it was reinforced that [Chol]Cl is not devoid of environmental hazardous potential regardless of its argued biological compatibility. Unique gene expression signatures could also be identified for each IL, enlightening specific mechanisms of toxicity.

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.

DNMT3.1 controls trade-offs between growth, reproduction, and life span under starved conditions in Daphnia magna

The cladoceran crustacean Daphnia has long been a model of energy allocation studies due to its important position in the trophic cascade of freshwater ecosystems. However, the loci for controlling energy allocation between life history traits still remain unknown. Here, we report CRISPR/Cas-mediated target mutagenesis of DNA methyltransferase 3.1 (DNMT3.1) that is upregulated in response to caloric restriction in Daphnia magna. The resulting biallelic mutant is viable and did not show any change in growth rate, reproduction, and longevity under nutrient rich conditions. In contrast, under starved conditions, the growth rate of this DNMT3.1 mutant was increased but its reproduction was reciprocally reduced compared to the wild type when the growth and reproduction activities competed during a period from instar 4 to 8. The life span of this mutant was significantly shorter than that of the wild type. We also compared transcriptomes between DNMT3.1 mutant and wild type under nutrient-rich and starved conditions. Consistent with the DNMT3.1 mutant phenotypes, the starved condition led to changes in the transcriptomes of the mutant including differential expression of vitellogenin genes. In addition, we found upregulation of the I am not dead yet (INDY) ortholog, which has been known to shorten the life span in Drosophila, explaining the shorter life span of the DNMT3.1 mutant. These results establish DNMT3.1 as a key regulator for life span and energy allocation between growth and reproduction during caloric restriction. Our findings reveal how energy allocation is implemented by selective expression of a DNMT3 ortholog that is widely distributed among animals. We also infer a previously unidentified adaptation of Daphnia that invests more energy for reproduction than growth under starved conditions.

Nutrient-imbalanced conditions shift the interplay between zooplankton and gut microbiota

Background

Nutrient stoichiometry of phytoplankton frequently changes with aquatic ambient nutrient concentrations, which is mainly influenced by anthropogenic water treatment and the ecosystem dynamics. Consequently, the stoichiometry of phytoplankton can markedly alter the metabolism and growth of zooplankton. However, the effects of nutrient-imbalanced prey on the interplay between zooplankton and their gut microbiota remain unknown. Using metatranscriptome, a 16 s rRNA amplicon-based neutral community model (NCM) and experimental validation, we investigated the interactions between Daphnia magna and its gut microbiota in a nutrient-imbalanced algal diet.

Results

Our results showed that in nutrient-depleted water, the nutrient-enriched zooplankton gut stimulated the accumulation of microbial polyphosphate in fecal pellets under phosphorus limitation and the microbial assimilation of ammonia under nitrogen limitation. Compared with the nutrient replete group, both N and P limitation markedly promoted the gene expression of the gut microbiome for organic matter degradation but repressed that for anaerobic metabolisms. In the nutrient limited diet, the gut microbial community exhibited a higher fit to NCM (R² = 0.624 and 0.781, for N- and P-limitation, respectively) when compared with the Control group (R² = 0.542), suggesting increased ambient-gut exchange process favored by compensatory feeding. Further, an additional axenic grazing experiment revealed that the growth of D. magna can still benefit from gut microbiota under a nutrient-imbalanced diet.

Conclusions

Together, these results demonstrated that under a nutrient-imbalanced diet, the microbes not only benefit themselves by absorbing excess nutrients inside the zooplankton gut but also help zooplankton to survive during nutrient limitation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07333-z.

Roundup causes embryonic development failure and alters metabolic pathways and gut microbiota functionality in non-target species

BACKGROUND

Research around the weedkiller Roundup is among the most contentious of the twenty-first century. Scientists have provided inconclusive evidence that the weedkiller causes cancer and other life-threatening diseases, while industry-paid research reports that the weedkiller has no adverse effect on humans or animals. Much of the controversial evidence on Roundup is rooted in the approach used to determine safe use of chemicals, defined by outdated toxicity tests. We apply a system biology approach to the biomedical and ecological model species Daphnia to quantify the impact of glyphosate and of its commercial formula, Roundup, on fitness, genome-wide transcription and gut microbiota, taking full advantage of clonal reproduction in Daphnia. We then apply machine learning-based statistical analysis to identify and prioritize correlations between genome-wide transcriptional and microbiota changes.

RESULTS

We demonstrate that chronic exposure to ecologically relevant concentrations of glyphosate and Roundup at the approved regulatory threshold for drinking water in the US induce embryonic developmental failure, induce significant DNA damage (genotoxicity), and interfere with signaling. Furthermore, chronic exposure to the weedkiller alters the gut microbiota functionality and composition interfering with carbon and fat metabolism, as well as homeostasis. Using the "Reactome," we identify conserved pathways across the Tree of Life, which are potential targets for Roundup in other species, including liver metabolism, inflammation pathways, and collagen degradation, responsible for the repair of wounds and tissue remodeling.

CONCLUSIONS

Our results show that chronic exposure to concentrations of Roundup and glyphosate at the approved regulatory threshold for drinking water causes embryonic development failure and alteration of key metabolic functions via direct effect on the host molecular processes and indirect effect on the gut microbiota. The ecological model species Daphnia occupies a central position in the food web of aquatic ecosystems, being the preferred food of small vertebrates and invertebrates as well as a grazer of algae and bacteria. The impact of the weedkiller on this keystone species has cascading effects on aquatic food webs, affecting their ability to deliver critical ecosystem services. Video Abstract.

Roles of and cross-talk between ecdysteroid and sesquiterpenoid pathways in embryogenesis of branchiopod crustacean Daphnia magna

The ecdysteroid and sesquiterpenoid pathways control growth, developmental transition, and embryogenesis in insects. However, the function of orthologous genes and the cross-talk between both pathways remain largely uncharacterized in non-insect arthropods. Spook (Spo) and Juvenile hormone acid o-methyltransferase (Jhamt) have been suggested to function as rate-limiting factors in ecdysteroid and sesquiterpenoid biosynthesis, respectively, in insects. In this study, we report on the functions of Spo and Jhamt and the cross-talk between them in embryos of the branchiopod crustacean Daphnia magna. Spo expression was activated at the onset of gastrulation, with the depletion of Spo transcript by RNAi resulting in developmental arrest at this stage. This phenotype could be partially rescued by supplementation with 20-hydroxyecdysone, indicating that Spo may play the same role in ecdysteroid biosynthesis in early embryos, as reported in insects. After hatching, Spo expression was repressed, while Jhamt expression was activated transiently, despite its silencing during other embryonic stages. Jhamt RNAi showed little effect on survival, but shortened the embryonic period. Exposure to the sesquiterpenoid analog Fenoxycarb extended the embryonic period and rescued the Jhamt RNAi phenotype, demonstrating a previously unidentified role of sesquiterpenoid in the repression of precocious embryogenesis. Interestingly, the knockdown of Jhamt resulted in the derepression of ecdysteroid biosynthesis genes, including Spo, similar to regulation during insect hormonal biosynthesis. Sesquiterpenoid signaling via the Methoprene-tolerant gene was found to be responsible for the repression of ecdysteroid biosynthesis genes. It upregulated an ortholog of CYP18a1 that degrades ecdysteroid in insects. These results illuminate the conserved and specific functions of the ecdysteroid and sesquiterpenoid pathways in Daphnia embryos. We also infer that the common ancestor of branchiopod crustaceans and insects exhibited antagonism between the two endocrine hormones before their divergence 400 million years ago.

Spatio-temporal patterns in the gene expression of the calanoid copepod Temora longicornis in the Belgian part of the North Sea

Marine zooplankton are increasingly being affected by recent environmental changes, such as climate change, and respond with profound spatial relocations and shifts in phenology and physiology. In order to predict whether populations are able to persist or adapt to such new conditions, it is essential to understand the molecular basis of such adaptations, which ultimately get translated into these physiological responses. To explore variation in population gene expression across time and space, we investigated transcriptome-level profiles of the calanoid copepod Temora longicornis, that were collected at four different locations in the Belgian Part of the North Sea (BPNS) on three different time points (April, June, October) in 2018. RNA-seq analysis of field collected adults identified large seasonal differences in gene expression, mainly between spring-summer and autumn samples. The largest log-fold changes occurred in a set of genes encoding for ribosomal and myosin (heavy chain) transcripts. Enrichment analysis revealed a strong seasonal pattern in vitellogenin, cuticle and glycolytic gene expression as well. We also found a positive correlation between vitellogenin expression and densities of T. longicornis. No clear spatial variation in expression patterns was found in the BPNS. This study underlines the potential of field gene expression studies for biomonitoring purposes and the significance of considering seasonal variation in future studies.

Two Doublesex1 mutants revealed a tunable gene network underlying intersexuality in Daphnia magna

In recent years, the binary definition of sex is being challenged by repetitive reports about individuals with ambiguous sexual identity from various animal groups. This has created an urge to decode the molecular mechanism underlying sexual development. However, sexual ambiguities are extremely uncommon in nature, limiting their experimental value. Here, we report the establishment of a genetically modified clone of Daphnia magna from which intersex daphniids can be readily generated. By mutating the conserved central sex determining factor Doublesex1, body-wide feminization of male daphniid could be achieved. Comparative transcriptomic analysis also revealed a genetic network correlated with Doublesex1 activity which may account for the establishment of sexual identity in D. magna. We found that Dsx1 repressed genes related to growth and promoted genes related to signaling. We infer that different intersex phenotypes are the results of fluctuation in activity of these Dsx1 downstream factors. Our results demonstrated that the D. magna genome is capable of expressing sex in a continuous array, supporting the idea that sex is actually a spectrum.

Assessment of Transcriptomic and Apical Responses of Daphnia magna Exposed to a Polyethylene Microplastic in a 21-d Chronic Study

There is global concern regarding the fate and effects of microplastics in the environment, particularly in aquatic systems. In the present study, ethylene acrylic acid copolymer particles were evaluated in a chronic toxicity study with the aquatic invertebrate Daphnia magna. The study design included a natural particle control treatment (silica) to differentiate any potential physical effects of a particle from the intrinsic toxicity of the test material. In addition to the standard endpoints of survival, growth, and reproduction, the transcriptomic profiles of control and ethylene acrylic acid copolymer-exposed D. magna were evaluated at the termination of the 21-d toxicity study. No significant effects on D. magna growth, survival, or reproduction were observed in comparison with both particle and untreated control groups. Significant transcriptomic alterations were induced at the highest treatment level of 2.3 × 10¹² particles of the ethylene acrylic acid copolymer/L in key pathways linked to central metabolism and energy reserves, oxidative stress, and ovulation and molting, indicating a global transcriptomic response pattern. To put the results in perspective is challenging at this time, because, to date, microplastic environmental monitoring approaches have not been equipped to detect particles in the nanosize range. However, our results indicate that ethylene acrylic acid copolymer microplastics in the upper nanosize range are not expected to adversely affect D. magna growth, survival, or reproductive outcomes at concentrations of up to 10¹² particles/L. Environ Toxicol Chem 2020;39:1578-1589. © 2020 SETAC.

Daphnia magna modifies its gene expression extensively in response to caloric restriction revealing a novel effect on haemoglobin isoform preference

Caloric restriction (CR) produces clear phenotypic effects within and between generations of the model crustacean Daphnia magna. We have previously established that micro-RNAs and cytosine methylation change in response to CR in this organism, and we demonstrate here that CR has a dramatic effect on gene expression. Over 6,000 genes were differentially expressed between CR and well-fed D. magna, with a bias towards up-regulation of genes under caloric restriction. We identified a highly expressed haemoglobin gene that responds to CR by changing isoform proportions. Specifically, a transcript containing three haem-binding erythrocruorin domains was strongly down-regulated under CR in favour of transcripts containing fewer or no such domains. This change in the haemoglobin mix is similar to the response to hypoxia in Daphnia, which is mediated through the transcription factor hypoxia-inducible factor 1, and ultimately the mTOR signalling pathway. This is the first report of a role for haemoglobin in the response to CR. We also observed high absolute expression of superoxide dismutase (SOD) in normally fed individuals, which contrasts with observations of high SOD levels under CR in other taxa. However, key differentially expressed genes, like SOD, were not targeted by differentially expressed micro-RNAs. Whether the link between haemoglobin and CR occurs in other organisms, or is related to the aquatic lifestyle, remains to be tested. It suggests that one response to CR may be to simply transport less oxygen and lower respiration.

Liver transcriptome resources of four commercially exploited teleost species

The generation of omic resources is central to develop adequate management strategies for species with economic value. Here, we provide high-coverage RNA-seq datasets of liver tissue (containing between 80,2 and 88,4 million of paired-end reads) from four wildtype teleost species with high commercial value: Trachurus trachurus (TTR; Atlantic horse mackerel), Scomber scombrus (SSC; Atlantic mackerel), Trisopterus luscus (TLU; pout), and Micromesistius poutassou (MPO; blue whiting). A comprehensive assembly pipeline, using de novo single and multi-kmer assembly approaches, produced 64 single high-quality liver transcriptomes - 16 per species. The final assemblies, with N50 values ranging from 2,543-3,700 bp and BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness values between 81.8-86.5% of the Actinopterygii gene set, were subjected to open reading frame (ORF) prediction and functional annotation. Our study provides the first transcriptomic resources for these species and offers valuable tools to evaluate both neutral and selected genetic variation among populations, and to identify candidate genes for environmental adaptation assisting in the investigation of the effects of global changes in fisheries.

An alternative route of bacterial infection associated with a novel resistance locus in the Daphnia-Pasteuria host-parasite system

To understand the mechanisms of antagonistic coevolution, it is crucial to identify the genetics of parasite resistance. In the Daphnia magna–Pasteuria ramosa host–parasite system, the most important step of the infection process is the one in which P. ramosa spores attach to the host’s foregut. A matching-allele model (MAM) describes the host–parasite genetic interactions underlying attachment success. Here we describe a new P. ramosa genotype, P15, which, unlike previously studied genotypes, attaches to the host’s hindgut, not to its foregut. Host resistance to P15 attachment shows great diversity across natural populations. In contrast to P. ramosa genotypes that use foregut attachment, P15 shows some quantitative variation in attachment success and does not always lead to successful infections, suggesting that hindgut attachment represents a less-efficient infection mechanism than foregut attachment. Using a Quantitative Trait Locus (QTL) approach, we detect two significant QTLs in the host genome: one that co-localizes with the previously described D. magna PR locus of resistance to foregut attachment, and a second, major QTL located in an unlinked genomic region. We find no evidence of epistasis. Fine mapping reveals a genomic region, the D locus, of ~13 kb. The discovery of a second P. ramosa attachment site and of a novel host-resistance locus increases the complexity of this system, with implications for both for the coevolutionary dynamics (e.g., Red Queen and the role of recombination), and for the evolution and epidemiology of the infection process.

A comprehensive epigenomic analysis of phenotypically distinguishable, genetically identical female and male Daphnia pulex

BACKGROUND

Daphnia species reproduce by cyclic parthenogenesis involving both sexual and asexual reproduction. The sex of the offspring is environmentally determined and mediated via endocrine signalling by the mother. Interestingly, male and female Daphnia can be genetically identical, yet display large differences in behaviour, morphology, lifespan and metabolic activity. Our goal was to integrate multiple omics datasets, including gene expression, splicing, histone modification and DNA methylation data generated from genetically identical female and male Daphnia pulex under controlled laboratory settings with the aim of achieving a better understanding of the underlying epigenetic factors that may contribute to the phenotypic differences observed between the two genders.

RESULTS

In this study we demonstrate that gene expression level is positively correlated with increased DNA methylation, and histone H3 trimethylation at lysine 4 (H3K4me3) at predicted promoter regions. Conversely, elevated histone H3 trimethylation at lysine 27 (H3K27me3), distributed across the entire transcript length, is negatively correlated with gene expression level. Interestingly, male Daphnia are dominated with epigenetic modifications that globally promote elevated gene expression, while female Daphnia are dominated with epigenetic modifications that reduce gene expression globally. For examples, CpG methylation (positively correlated with gene expression level) is significantly higher in almost all differentially methylated sites in male compared to female Daphnia. Furthermore, H3K4me3 modifications are higher in male compared to female Daphnia in more than 3/4 of the differentially regulated promoters. On the other hand, H3K27me3 is higher in female compared to male Daphnia in more than 5/6 of differentially modified sites. However, both sexes demonstrate roughly equal number of genes that are up-regulated in one gender compared to the other sex. Since, gene expression analyses typically assume that most genes are expressed at equal level among samples and different conditions, and thus cannot detect global changes affecting most genes.

CONCLUSIONS

The epigenetic differences between male and female in Daphnia pulex are vast and dominated by changes that promote elevated gene expression in male Daphnia. Furthermore, the differences observed in both gene expression changes and epigenetic modifications between the genders relate to pathways that are physiologically relevant to the observed phenotypic differences.

Hemoglobin in Arthropods-Daphnia as a Model

Hemoglobin is the respiratory protein of many arthropods, enhancing the oxygen transport capacity of the hemolymph. One example, that has been subject of extensive studies, is the hemoglobin of the crustacean genus Daphnia. Here the characteristics of this oxygen binding protein are reviewed. The genetic structure is the result of repeated duplication events in the evolution, leading to a variety of di-domain isoforms. Adjustments to environmental changes thus result from differential expression of these paralogs. The biochemical properties, including spectral characteristics, concentration ranges, molecular mass of monomers and native oligomers, are compared. Structural differences between isoforms can be correlated to functional properties of oxygen binding characteristics. The mechanism of hemoglobin induction via hypoxia-inducible factor 1 allows the response to altered oxygen and temperature conditions. Changes of the hemoglobin suite in quantity and functional quality can be linked to their benefits for the animals' physiological performance. However, there is a large inter- and intra-specific variability of this induction potential. The consequences of altered hemoglobin characteristics for the animals' success within their habitat are discussed.

The Mercury-Tolerant Microbiota of the Zooplankton Daphnia Aids in Host Survival and Maintains Fecundity under Mercury Stress

Many aquatic organisms can thrive in polluted environments by having the genetic capability to withstand suboptimal conditions. However, the contributions of microbiomes under these stressful environments are poorly understood. We investigated whether a mercury-tolerant microbiota can extend its phenotype to its host by ameliorating host survival and fecundity under mercury-stress. We isolated microbiota members from various clones of Daphnia magna, screened for the mercury-biotransforming merA gene, and determined their mercury tolerance levels. We then introduced the mercury-tolerant microbiota, Pseudomonas-10, to axenic D. magna and quantified its merA gene expression, mercury reduction capability, and measured its impact on host survival and fecundity. The expression of the merA gene was up-regulated in Pseudomonas-10, both in isolation and in host-association with mercury exposure. Pseudomonas-10 is also capable of significantly reducing mercury concentration in the medium. Notably, mercury-exposed daphnids containing only Pseudomonas-10 exhibited higher survival and fecundity than mercury-exposed daphnids supplemented with parental microbiome. Our study showed that zooplankton, such as Daphnia, naturally harbor microbiome members that are eco-responsive and tolerant to mercury exposure and can aid in host survival and maintain host fecundity in a mercury-contaminated environment. This study further demonstrates that under stressful environmental conditions, the fitness of the host can depend on the genotype and the phenotype of its microbiome.

Transcriptome sequencing of a keystone aquatic herbivore yields insights on the temperature-dependent metabolism of essential lipids

Background

Nutritional quality of phytoplankton is a major determinant of the trophic transfer efficiency at the plant-herbivore interface in freshwater food webs. In particular, the phytoplankton’s content of the essential polyunsaturated omega-3 fatty acid eicosapentaenoic acid (EPA) has been repeatedly shown to limit secondary production in the major zooplankton herbivore genus Daphnia. Despite extensive research efforts on the biological model organism Daphnia, and the availability of several Daphnia genomes, little is known regarding the molecular mechanisms underlying the limitations in Daphnia related to dietary EPA availability.

Results

We used RNA-seq to analyse the transcriptomic response of Daphnia magna which were fed with two different diets — each with or without supplementation of EPA — at two different temperature levels (15 and 20 °C). The transcripts were mapped to the D. magna genome assembly version 2.4, containing 26,646 translations. When D. magna fed on green alga, changing the temperature provoked a differential expression of 2001 transcripts, and in cyanobacteria-fed daphnia, 3385 transcripts were affected. The supplementation of EPA affected 1635 (on the green algal diet), or 175 transcripts (on the cyanobacterial diet), respectively. Combined effects for diet and temperature were also observed (669 for the green algal and 128 transcripts for the cyanobacterial diet). Searching for orthologous genes (COG-analysis) yielded a functional overview of the altered transcriptomes. Cross-matched transcript sets from both feed types were compiled to illuminate core responses to the factors temperature and EPA-supplementation.

Conclusions

Our highly controlled eco-physiological experiments revealed an orchestrated response of genes involved in the transformation and signalling of essential fatty acids, including eicosanoid-signalling pathways with potential immune functions. We provide an overview of downstream-regulated genes, which contribute to enhance growth and reproductive output. We also identified numerous EPA-responsive candidate genes of yet unknown function, which constitute new targets for future studies on the molecular basis of EPA-dependent effects at the freshwater plant-herbivore interface.

Effects of Single and Combined Low Concentrations of Neuroactive Drugs on Daphnia magna Reproduction and Transcriptomic Responses

Assessing the risk of neuroactive pharmaceuticals in the environment requires an understanding of their joint effects at low concentrations across species. Here, we assessed reproductive and transcriptional effects of single and ternary equi-effective mixture exposure to propranolol, diazepam, and carbamazepine on the crustacean Daphnia magna at environmentally relevant concentrations. The three compounds enhanced reproduction in adults and induced specific transcriptome changes in preadolescent individuals. Comparison of the results from single exposures to a ternary equi-effective mixture of the three compounds showed additive action. Transcriptomic analyses identified 3248 genes affected by at least one of the treatments, which were grouped into four clusters. Two clusters (1897 gene transcripts in total) behaved similarly, appearing either over- or under-represented relative to control, in all single and mixture treatments. The third and fourth clusters grouped genes differently transcribed upon exposure to diazepam and propranolol, respectively. Functional transcriptomics analysis indicated that the four clusters shared major deregulated signaling pathways implicated on energy, growth, reproduction, and neurologically related processes, which may be responsible for the observed reproductive effects. Thus, our study showed additive effects at the transcriptional and physiological level and provides a novel approach to the analysis of environmentally relevant mixtures of neuroactive compounds.

Within species expressed genetic variability and gene expression response to different temperatures in the rotifer Brachionus calyciflorus sensu stricto

Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).

Species-specific transcriptomic responses in Daphnia magna exposed to a bio-plastic production intermediate

Hydroxymethylfurfural (HMF) is a plant-based chemical building block that could potentially substitute petroleum-based equivalents, yet ecotoxicological data of this compound is currently limited. In this study, the effects of HMF on the reproduction and survival of Daphnia magna were assessed through validated ecotoxicological tests. The mechanism of toxicity was determined by analysis of transcriptomic responses induced by exposure to different concentrations of HMF using RNA sequencing. HMF exerted toxicity to D. magna with an EC₅₀ for effects on reproduction of 17.2 mg/l. HMF exposure affected molecular pathways including sugar and polysaccharide metabolism, lipid metabolism, general stress metabolism and red blood cell metabolism, although most molecular pathways affected by HMF exposure were dose specific. Hemoglobin genes, however, responded in a sensitive and dose-related manner. No induction of genes involved in the xenobiotic metabolism or oxidative stress metabolism pathway could be observed, which contrasted earlier observations on transcriptional responses of the terrestrial model Folsomia candida exposed to the same compound in a similar dose. We found 4189 orthologue genes between D. magna and F. candida, yet only twenty-one genes of those orthologues were co-regulated in both species. The contrasting transcriptional responses to the same compound exposed at a similar dose between D. magna and F. candida indicates limited overlap in stress responses among soil and aquatic invertebrates. The dose-related expression of hemoglobin provides further support for using hemoglobin expression as a biomarker for general stress responses in daphnids.

Daphnia stressor database: Taking advantage of a decade of Daphnia '-omics' data for gene annotation

Gene expression patterns help to measure and characterize the effect of environmental perturbations at the cellular and organism-level. Complicating interpretation is the presence of uncharacterized or "hypothetical" gene functions for a large percentage of genomes. This is particularly evident in Daphnia genomes, which contains many regions coding for "hypothetical proteins" and are significantly divergent from many of the available arthropod model species, but might be ecologically important. In the present study, we developed a gene expression database, the Daphnia stressor database (http://www.daphnia-stressordb.uni-hamburg.de/dsdbstart.php), built from 90 published studies on Daphnia gene expression. Using a comparative genomics approach, we used the database to annotate D. galeata transcripts. The extensive body of literature available for Daphnia species allowed to associate stressors with gene expression patterns. We believe that our stressor based annotation strategy allows for better understanding and interpretation of the functional role of the understudied hypothetical or uncharacterized Daphnia genes, thereby increasing our understanding of Daphnia's genetic and phenotypic variability.

Genetic control of male production in Daphnia pulex

Daphnia normally reproduce by cyclical parthenogenesis, with offspring sex being determined by environmental cues. However, some females have lost the ability to produce males. Our results demonstrate that this loss of male-producing ability is controlled by a dominant allele at a single locus. We identified the locus by comparing whole-genome sequences of 67 nonmale-producing (NMP) and 100 male-producing (MP) clones from 5 Daphnia pulex populations, revealing 132 NMP-linked SNPs and 59 NMP-linked indels within a single 1.1-Mb nonrecombining region on chromosome I. These markers include 7 nonsynonymous mutations, all of which are located within one unannotated protein-coding gene (gene 8960). Within this single gene, all of the marker-linked NMP haplotypes from different populations form a monophyletic clade, suggesting a single origin of the NMP phenotype, with the NMP haplotype originating by introgression from a sister species, Daphnia pulicaria Methyl farnesoate (MF) is the innate juvenile hormone in daphnids, which induces the production of males and whose inhibition results in female-only production. Gene 8960 is sensitive to treatment by MF in MP clones, but such responsiveness is greatly reduced in NMP clones. Thus, we hypothesize that gene 8960 is located downstream of the MF-signaling pathway in D. pulex, with the NMP phenotype being caused by expression change of gene 8960.

Epigenetic, transcriptional and phenotypic responses in two generations of Daphnia magna exposed to the DNA methylation inhibitor 5-azacytidine

The water flea Daphnia magna is a keystone species in freshwater ecosystems and has been widely used as a model organism in environmental ecotoxicology. This aquatic crustacean is sensitive to environmental stressors and displays considerable plasticity in adapting to changing environmental conditions. Part of this plasticity may be due to epigenetic regulation of gene expression, including changes to DNA methylation and histone modifications. Because of the generally hypomethylated genome of this species, we hypothesized that the histone code may have an essential role in the epigenetic control and that histone modifications might be an early marker for stress. This study aims to characterize the epigenetic, transcriptional and phenotypic responses and their causal linkages in directly exposed adult (F0) Daphnia and peritoneal exposed neonates (F1) after a chronic (7-day) exposure to a sublethal concentration (10 mg/l) of 5-azacytidine, a well-studied vertebrate DNA methylation inhibitor. Exposure of the F0 generation significantly reduced the cumulative fecundity, accompanied with differential expression of genes in the one-carbon-cycle metabolic pathway. In the epigenome of the F0 generation, a decrease in global DNA methylation, but no significant changes on H3K4me3 or H3K27me3, were observed. In the F1 offspring generation, changes in gene expression, a significant reduction in global DNA methylation and changes in histone modifications were identified. The results indicate that exposure during adulthood may result in more pronounced effects on early development in the offspring generation, though interpretation of the data should be carefully done since both the exposure regime and developmental period is different in the two generations examined. The obtained results improve our understanding of crustacean epigenetics and the tools developed may promote use of epigenetic markers in hazard assessment of environmental stressors.

Evolution of Gene Expression during a Transition from Environmental to Genetic Sex Determination

Genetic sex determination (GSD) can evolve from environmental sex determination (ESD) via an intermediate state in which both coexist in the same population. Such mixed populations are found in the crustacean Daphnia magna, where non-male producers (NMP, genetically determined females) coexist with male producers (MP), in which male production is environmentally inducible and can also artificially be triggered by exposure to juvenile hormone. This makes Daphnia magna a rare model species for the study of evolutionary transitions from ESD to GSD. Although the chromosomal location of the NMP-determining mutation has been mapped, the actual genes and pathways involved in the evolution of GSD from ESD remain unknown. Here, we present a transcriptomic analysis of MP and NMP females under control (female producing) and under hormone exposure conditions. We found ∼100 differentially expressed genes between MP and NMP under control conditions. Genes in the NMP-determining chromosome region were especially likely to show such constitutive expression differences. Hormone exposure led to expression changes of an additional ∼100 (MP) to ∼600 (NMP) genes, with an almost systematic upregulation of those genes in NMP. These observations suggest that the NMP phenotype is not determined by a simple "loss-of-function" mutation. Rather, homeostasis of female offspring production under hormone exposure appears to be an active state, tightly regulated by complex mechanisms involving many genes. In a broader view, this illustrates that the evolution of GSD, while potentially initiated by a single mutation, likely leads to secondary integration involving many genes and pathways.

Genome-Wide Stress Responses to Copper and Arsenic in a Field Population of Daphnia

Over the past decade, significant advances have been made to unravel molecular mechanisms of stress response in different ecotoxicological model species. Within this study, we focus on population level transcriptomic responses of a natural population of Daphnia magna Straus, (1820), to heavy metals. We aim to characterize the population level transcriptomic responses, which include standing genetic variation, and improve our understanding on how populations respond to environmental stress at a molecular level. We studied population level responses to two heavy metals, copper and arsenic, and their binary mixture across time. Transcriptomic patterns identified significantly regulated gene families and genes at the population level including cuticle proteins and resilins. Furthermore, some of these differentially regulated gene families, such as cuticle proteins, were also significantly enriched for genetic variations including SNPs and MNPs. In general, genetic variation was observed in specific gene families, many of which are known to be involved in stress response. Overall, our results indicate that molecular stress responses can be identified within natural populations and that linking molecular mechanisms with genetic variation at the population level could contribute significantly to adverse outcome frameworks.

Genome-wide methylation is modified by caloric restriction in Daphnia magna

BACKGROUND

The degradation of epigenetic control with age is associated with progressive diseases of ageing, including cancers, immunodeficiency and diabetes. Reduced caloric intake slows the effects of ageing and age-related disease in vertebrates and invertebrates, a process potentially mediated by the impact of caloric restriction on epigenetic factors such as DNA methylation. We used whole genome bisulphite sequencing to study how DNA methylation patterns change with diet in a small invertebrate, the crustacean Daphnia magna. Daphnia show the classic response of longer life under caloric restriction (CR), and they reproduce clonally, which permits the study of epigenetic changes in the absence of genetic variation.

RESULTS

Global cytosine followed by guanine (CpG) methylation was 0.7-0.9%, and there was no difference in overall methylation levels between normal and calorie restricted replicates. However, 333 differentially methylated regions (DMRs) were evident between the normally fed and CR replicates post-filtering. Of these 65% were hypomethylated in the CR group, and 35% were hypermethylated in the CR group.

CONCLUSIONS

Our results demonstrate an effect of CR on the genome-wide methylation profile. This adds to a growing body of research in Daphnia magna that demonstrate an epigenomic response to environmental stimuli. Specifically, gene Ontology (GO) term enrichment of genes associated with hyper and hypo-methylated DMRs showed significant enrichment for methylation and acyl-CoA dehydrogenase activity, which are linked to current understanding of their roles in CR in invertebrate model organisms.

Tryptophan hydroxylase (TRH) loss of function mutations in Daphnia deregulated growth, energetic, serotoninergic and arachidonic acid metabolic signalling pathways

Serotonin has a pivotal function regulating development, growth, reproduction and behavior in animals. In this paper, we studied the deregulatory effects of the deprivation of serotonin in Daphnia magna TRH CRISPR-Cas9 mutants. Bi-allelic in-del THR mutants and, to a lesser extent, mono-allelic ones grew less, reproduced later, and produced smaller clutches than wild type clones. Transcriptomic and functional gene analyses showed a down-regulation of growth/molting and energy metabolism signaling pathways in TRH mutants, while revealing marked differences between mono- and bi-allelic clones. Bi-allelic mutants, lacking serotonin, presented the serotonergic synapse and arachidonic acid metabolic pathways down-regulated while the tryptophan to kynurenine was upregulated, thus indicating a cross-talk between the serotonergic and arachidonic acid metabolic pathways. Finally, the effects on the insulin growth factor-mediated signaling pathway were marginal. These changes in functional and metabolic pathways are consistent with previously reported effects in D. magna exposed to pharmaceuticals that inhibited arachidonic metabolism or enhanced the levels of serotonin.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.

Mixtures of Aluminum and Indium Induce More than Additive Phenotypic and Toxicogenomic Responses in Daphnia magna

Aquatic systems are contaminated by many metals but their effects as mixtures on organisms are not well understood. Here, we assessed effects of aluminum with fairly well-known modes of actions and indium, an understudied emerging contaminant from electronics, followed by studying equi-effective mixtures thereof. We report acute and adverse phenotypic effects in Daphnia magna adults and global transcriptomic effects employing RNA sequencing in neonates. The mixture induced more than additive activity in mortality and in physiological effects, including growth and reproduction. Similarly, transcriptomic effects were more than additive, as indicated by a markedly higher number of 463 differentially expressed transcripts in the mixture and by distinct classes of genes assigned to several biological functions, including metabolic processes, suggesting depleted energy reserves, which may be responsible for the observed impaired reproduction and growth. A gene set enrichment analysis (GSEA) of a priori known response pathways for aluminum confirmed activation of distinct molecular pathways by indium. Our study is highlighting more than additive effects at the transcriptional and physiological level and is providing a state-of-the art approach to mixture analysis, which is important for risk assessment of these metals and metal mixtures.