Gene response profiles for Daphnia pulex exposed to the environmental stressor cadmium reveals novel crustacean metallothioneins

History of exposure to copper influences transgenerational gene expression responses in Daphnia magna

The establishment of transgenerational effects following chemical exposure is a powerful phenomenon, capable of modulating ecosystem health beyond exposure periods. This study assessed the transgenerational effects occurring due to copper exposure in the invertebrate D. magna at the transcriptional level, while evaluating the role of exposure history on such responses. Thus, daphnids acclimated for several generations in a copper vs. clean medium were then exposed for one generation (F0) to this metal, and monitored for the following non-exposed generations (F1, F2 and F3). Organisms differing in exposure histories showed remarkably different transcriptional profiles at the F0, with naïve organisms being more profoundly affected. These trends were confirmed for F3 treatments, which presented different transcriptional patterns for genes involved in detoxification, oxidative stress, DNA damage repair, circadian clock functioning and epigenetic regulation. Furthermore, regardless of exposure history, a great number of histone modifier genes were always found transcriptionally altered, thus suggesting the involvement of histone modifications in the response of Daphnia to metal exposure. Lastly, remarkably distinct transgenerational transcriptional responses were found between naïve and non-naïve organisms, thereby highlighting the influence of exposure history on gene expression and confirming the capacity of metals to determine transgenerational transcriptional effects across non-exposed generations.

Genome-wide annotation and comparative analysis revealed conserved cuticular protein evolution among non-biting midges with varied environmental adaptability

Chironomidae, non-biting midges, a diverse and abundant insect group in global aquatic ecosystems, represent an exceptional model for investigating genetic adaptability mechanisms in aquatic insects due to their extensive species diversity and resilience to various environmental conditions. The cuticle in insects acts as the primary defense against ecological pressures. Cuticular Proteins (CPs) determine cuticle characteristics, facilitating adaptation to diverse challenges. However, systematic annotation of CP genes has only been conducted for one Chironomidae species, Propsilocerus akamusi, by our team. In this study, we expanded this annotation by identifying CP genes in eight additional Chironomidae species, covering all Chironomidae species with available genome data. We identified a total of 889 CP genes, neatly categorized into nine CP families: 215 CPR RR1 genes, 272 CPR RR2 genes, 23 CPR RR3 genes, 21 CPF genes, 16 CPLCA genes, 19 CPLCG genes, 28 CPLCP genes, 77 CPAP genes, and 37 Tweedle genes. Subsequently, we conducted a comprehensive phylogenetic analysis of CPs within the Chironomidae family. This expanded annotation of CP genes across diverse Chironomidae species significantly contributes to our understanding of their remarkable adaptability.

Efficient CRISPR genome editing and integrative genomic analyses reveal the mosaicism of Cas-induced mutations and pleiotropic effects of scarlet gene in an emerging model system

Despite the revolutionary impacts of CRISPR-Cas gene editing systems, the effective and widespread use of CRISPR technologies in emerging model organisms still faces significant challenges. These include the inefficiency in generating heritable mutations at the organismal level, limited knowledge about the genomic consequences of gene editing, and an inadequate understanding of the inheritance patterns of CRISPR-Cas-induced mutations. This study addresses these issues by 1) developing an efficient microinjection delivery method for CRISPR editing in the microcrustacean Daphnia pulex; 2) assessing the editing efficiency of Cas9 and Cas12a nucleases, examining mutation inheritance patterns, and analyzing the local and global mutation spectrum in the scarlet mutants; and 3) investigating the transcriptomes of scarlet mutants to understand the pleiotropic effects of scarlet underlying their swimming behavior changes. Our reengineered CRISPR microinjection method results in efficient biallelic editing with both nucleases. While indels are dominant in Cas-induced mutations, a few on-site large deletions (>1kb) are observed, most likely caused by microhomology-mediated end joining repair. Knock-in of a stop codon cassette to the scarlet locus was successful, despite complex induced mutations surrounding the target site. Moreover, extensive germline mosaicism exists in some mutants, which unexpectedly produce different phenotypes/genotypes in their asexual progenies. Lastly, our transcriptomic analyses unveil significant gene expression changes associated with scarlet knock-out and altered swimming behavior in mutants, including several genes (e.g., NMDA1, ABAT, CNTNAP2) involved in human neurodegenerative diseases. This study expands our understanding of the dynamics of gene editing in the tractable model organism Daphnia and highlights its promising potential as a neurological disease model.

Genome-wide annotation of cuticular protein genes in non-biting midge Propsilocerus akamusi and transcriptome analysis of their response to heavy metal pollution

The insect cuticle is a sophisticated chitin-protein extracellular structure for mutable functions. The cuticles varied their structures and properties in different species, and the same species but in different regions or at different stages, to fill the requirements of different functions. The alteration of cuticle structures may also be induced due to challenges by some environmental crises, such as pollution exposures. The physical properties of the cuticle were determined by the cuticle proteins (CPs) they contain. The cuticle proteins are large protein groups in all insects, which are commonly divided into different families according to their conserved protein sequence motifs. Although Chironomidae is an abundant and universal insect in global aquatic ecosystems and a popular model for aquatic toxicology, no systematic annotation of CPs was done for any species in Chironomidae before. In this work, we annotated the CP genes of Propsilocerus akamusi, the most abundant Chironomidae species in Asia. A total of 160 CP genes were identified, and 97 of them could be well classified into eight CP families: 76 CPR genes can be subdivided into three groups (further divided into three subgroups: 36 RR1 genes, 37 RR2 genes, and 3 RR3 genes), 2 CPF genes, 3 CPLCA genes, 1 CPLCG gene, 8 CPAP genes, and 3 Tweedle genes. Additionally, we analyzed the response of P. akamusi CP genes at expression level to Cu exposure, which is related to the high heavy metal tolerance and the earlier onset of pupariation in heavy metal polluted water.

Characterization of cuticular proteins in CPR family in the wolf spider, Pardosa pseudoannulata, and the response of one subfamily genes to environmental stresses

Cuticular protein (CP) plays an essential role in the construction and function of exoskeleton in arthropods. CPR family, CP with Rebers and Riddiford (R&R) Consensus, is the largest CP family in insects, but it lacks systematic research in non-insect arthropods. In this study, we explored CPRs in the wolf spider, Pardosa pseudoannulata, a predator to many insect pests. We totally identified 152 CPRs in P. pseudoannulata genome, which were divided into two subgroups based on R&R Consensus sequences, with 12 CPRs in RR-1 and 140 in RR-2. All RR-2 members presented a novel Consensus with 34 amino acids, G-x(8)-G-x(6)-Y-x-A-x(3)-G-x(7)-N-E-x-G, which was a common characteristic for RR-2 CPRs in chelicerates. Transcriptome data was used to document the expression patterns of CPR genes in different tissues and ecdysis processes. The specific expressions were found for part CPR genes, such as five RR-2 genes that were specifically expressed in male genital bulbs and eleven RR-1 genes that were highly expressed in the integument. Due to the limited number and integument-specific expression of RR-1 genes, we further analyzed their responses to different environmental stresses at the transcriptional level. Except for PapsCPR11, ten RR-1 genes responded to at least one environmental stress, among with the expression of PapsCPR12 was significantly changed by three stresses (dryness, low temperature and imidacloprid treatments). Silencing PapsCPR12 increased the tolerance of P. pseudoannulata to imidacloprid. Overall, the results presented novel Consensus characteristics of CPRs in P. pseudoannulata, which was helpful for the identification and evolution analysis of CPRs in non-insect arthropods.

Emerging trends in nanoparticle toxicity and the significance of using Daphnia as a model organism

Nanoparticle production is on the rise due to its many uses in the burgeoning nanotechnology industry. Although nanoparticles have growing applications, there is great concern over their environmental impact due to their inevitable release into the environment. With uncertainty of environmental concentration and risk to aquatic organisms, the microcrustacean Daphnia spp. has emerged as an important freshwater model organism for risk assessment of nanoparticles because of its biological properties, including parthenogenetic reproduction; small size and short generation time; wide range of endpoints for ecotoxicological studies; known genome, useful for providing mechanistic information; and high sensitivity to environmental contaminants and other stressors. In this review, we (1) highlight the advantages of using Daphnia as an experimental model organism for nanotoxicity studies, (2) summarize the impacts of nanoparticle physicochemical characteristics on toxicity in relation to Daphnia, and (3) summarize the effects of nanoparticles (including nanoplastics) on Daphnia as well as mechanisms of toxicity, and (4) highlight research uncertainties and recommend future directions necessary to develop a deeper understanding of the fate and toxicity of nanoparticles and for the development of safer and more sustainable nanotechnology.

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.

Genome-Wide Analysis of Cadmium-Induced, Germline Mutations in a Long-Term Daphnia pulex Mutation-Accumulation Experiment

BACKGROUND

Germline mutations provide the raw material for all evolutionary processes and contribute to the occurrence of spontaneous human diseases and disorders. Yet despite the daily interaction of humans and other organisms with an increasing number of chemicals that are potentially mutagenic, precise measurements of chemically induced changes to the genome-wide rate and spectrum of germline mutation are lacking.

OBJECTIVES

A large-scale Daphnia pulex mutation-accumulation experiment was propagated in the presence and absence of an environmentally relevant cadmium concentration to quantify the influence of cadmium on germline mutation rates and spectra.

RESULTS

Cadmium exposure dramatically changed the genome-wide rates and regional spectra of germline mutations. In comparison with those in control conditions, Daphnia exposed to cadmium had a higher overall A : T → G : C mutation rates and a lower overall C : G → G : C mutation rate. Daphnia exposed to cadmium had a higher intergenic mutation rate and a lower exonic mutation rate. The higher intergenic mutation rate under cadmium exposure was the result of an elevated intergenic A : T → G : C rate, whereas the lower exon mutation rate in cadmium was the result of a complete loss of exonic C : G → G : C mutations-mutations that are known to be enriched at 5-hydroxymethylcytosine. We experimentally show that cadmium exposure significantly reduced 5-hydroxymethylcytosine levels.

DISCUSSION

These results provide evidence that cadmium changes regional mutation rates and can influence regional rates by interfering with an epigenetic process in the Daphnia pulex germline. We further suggest these observed cadmium-induced changes to the Daphnia germline mutation rate may be explained by cadmium's inhibition of zinc-containing domains. The cadmium-induced changes to germline mutation rates and spectra we report provide a comprehensive view of the mutagenic perils of cadmium and give insight into its potential impact on human population health. https://doi.org/10.1289/EHP8932.

Daphnia populations are similar but not identical in tolerance to different protease inhibitors

Cyanobacterial blooms often produce different classes and chemical variants of toxins such as dietary protease inhibitors (PIs) that affect the keystone grazer Daphnia. However, it has been shown that Daphnia populations are able to locally adapt to frequently occurring dietary PIs by modulating their digestive proteases. Up until now, local adaptation has exclusively been tested by making use of single cyanobacterial strains and by measuring average population tolerance. In contrast, we measured juvenile somatic growth rates and egg numbers of several individual clones per each of three different D. magna populations that have previously been found to be either tolerant or sensitive to the Microcystis strain BM25. Clones from the three D. magna populations were either treated with BM25 that produces three different protease inhibitor variants of the class of Ahp-cyclodepsipeptides or another Microcystis strain that produces two other Ahp-cyclodepsipeptide variants. Subsequently, the population growth was calculated as mean of the single-clone growth rates. Both tolerant populations (which originate from ponds with a cyanobacterial history) proved to be similarly tolerant to both Microcystis strains. However, single genotypes of the populations differed in their response to the different strains. Both the tolerant and the sensitive populations contained both sensitive and tolerant genotypes but in different proportions. Furthermore, the genotypes from the sensitive population showed a higher variance in response to one or both strains. Trade-offs between somatic growth rate and clutch size were found in one of the tolerant populations that originated from a pond where cyanobacteria were frequent in the past but completely absent since the pond's restoration. Because of those intra-population difference, we conclude that the tolerant populations were putatively selected by different Ahp-cyclodepsipeptide variants in the past and that all populations still possess the potential to adapt to other environmental conditions by genotype frequency shifts.

Reversible and irreversible transgenerational effects of metal exposure on nine generations of a tropical micro-crustacean

Micro-crustaceans are important grazers that control the algal blooms in eutrophic lakes. However, we know little about how these key species may be affected by long-term exposure to contaminants and when the transgenerational effects are reversible and irreversible. To address this, we investigated the effects of lead (Pb, 100 μg L^-1) exposure on morphology and reproduction of Moina dubia for nine consecutive generations (F1-F9) in three treatments: control, Pb, and pPb (M. dubia from Pb-exposed parents returned to the control condition). In F1-F2, Pb did not affect morphology, and reproduction of M. dubia. In all later generations, Pb-exposed M. dubia had a smaller body and shorter antennae than those in control. In F3-F6, pPb-exposed animals showed no differences in body size and antennae compared to the control, suggesting recoverable effects. In F7-F9, the body size and antennae of pPb-exposed animals did not differ compared to Pb-exposed ones, and both were smaller than the control animals, suggesting irreversible effects. Pb exposure reduced the brood size, number of broods and total neonates per female in F3-F9, yet the reproduction could recover in pPb treatment until F7. No recovery of the brood size and number of broods per female was observed in pPb-exposed animals in the F8-F9. Our study suggests that long-term exposure to metals, here Pb, may cause irreversible impairments in morphology and reproduction of tropical urban micro-crustaceans that may lower the top-down control on algal blooms and functioning of eutrophic urban lakes.

A comprehensive system for detection of behavioral change of D. magna exposed to various chemicals

The purpose of the present study was to develop a sensitive and comprehensive method, based on D. magna swimming behavior, for toxicity assessment of environmental chemicals. Firstly, D. magna swimming in several chambers with different diameters were compared to determine the most suitable container, and then baseline behaviors during light/dark periods as well as reactions to light/dark switching and vibration stimulation were determined. Secondly, after exposure to sub-lethal concentrations of the selected 42 typical chemicals, which were classified into heavy metals, pesticides, fungicides and flame retardants, the alterations in the swimming parameters were evaluated. Our results indicated the 48-well plate was the most suitable chamber for behavioral monitoring of D. magna, and specific responsive patterns of D. magna neonates to light/dark switching and vibration stimulation were observed. The results of the behavioral assays of chemicals suggested that D. magna was the most sensitive to methylmercury-chloride and then to abamectin and chlorpyrifos. The three chemicals at several to dozens of ng/L significantly changed swimming behaviors of D. magna. Furthermore, the alteration in the behavioral parameters (average swimming speed, etc.) induced by the selected chemicals could be ascribed to various modes of actions, confirming the reliability and practicability of the monitoring method.

Nitrate-N-mediated toxicological responses of Scenedesmus acutus and Daphnia pulex to cadmium, arsenic and their binary mixture (Cd/Asmix) at environmentally relevant concentrations

Several biomarkers used for ecological risk assessment have been established for single contaminant toxicity, many of which are less predictive of the influence of media and/or dietary nutrients on toxicity outcomes of contaminant mixtures. In this study, we investigate toxicological responses and life traits of Scenedesmus acutus and Daphnia pulex to heavy metals (cadmium-Cd, arsenic-As, binary mixture-Cd/As_mix) in media and diets with varied nutrient (nitrate-N) conditions (low-LN, median-MN, optimum-COMBO). Results showed that nitrate-N-mediated metal inhibitory effects on growth and productivity of primary producer (S. acutus) were significantly interactive (p < 0.05; effect size, ƞ²≤56 %). Cadmium toxicities (Cd-IC₅₀s) in S. acutus were 1.2×, 5.3×, and 4.3× As-IC₅₀s in LN, MN and COMBO media, respectively, while mixture (Cd/As_mix) toxicities were synergistic in MN medium and partial additivity in COMBO and LN media. Nitrate-N and metal exposure effects on S. acutus nutrient stoichiometry, metal uptake and bioaccumulation were significantly interactive (p < 0.05, ƞ²≤100 %). Moreover, survival of primary consumer (D. pulex) was significantly impaired by single and mixed dietary-metal exposures with greater effect under LN condition coupled with significant interactive effects on reproductive capacity (p < 0.05, ƞ²≤21.2 %) but not on swimming activity. We recommend that nitrate-N-mediated metal exposure effects/toxicity in bioindicator species should be considered during ecological risk assessments.

Differential gene expression profile of male and female copepods in response to cadmium exposure

In this study, the whole transcriptome and sex-specific differential gene expression of the copepod Pseudodiaptomus annandalei exposed to cadmium (Cd) were investigated. P. annandalei were exposed to 40 μg/L Cd from the naupliar stage to male and female adults. High-throughput transcriptome sequencing (RNA-seq) was performed with copepod samples using an Illumina Hiseq™ 2000 platform. TransDecoder analysis found 32,625 putative open reading frame contigs. At p-values of <0.001, a total of 4756 differentially expressed genes (DEGs) (2216 up-regulated and 2540 down-regulated genes) were found in male copepods. Whereas a total of 2879 DEGs (2007 up-regulated and 872 down-regulated genes) were found in female copepods. A few selected cellular stress response genes, involved in xenobiotic metabolism, energy metabolism, growth, and development as a result of Cd exposure in the copepods were discussed. The study showed that most of these processes were changed in a sex-specific manner, accounting for the different sensitivities of male and female copepods. Results suggest and reinforce that sex is an important factor to be considered in ecotoxicogenomics.

Protective effect assessment of Moringa oleifera against cadmium-induced toxicity in HCT116 and HEK293 cell lines

The cadmium (Cd) is considered one of the widespread toxic metals in the aquatic and terrestrial environments, which is due to its long half-life, non-degradable characteristic, and toxicity. Aqueous extract of freeze-dried Moringa oleifera (Moringaceae family) leaves was examined for protective effect and antioxidant power against Cd toxicity. The results revealed that Moringa aqueous extract (MAE) has contents of total polyphenols and flavonoids about 30.14 mg GAE/g and 18.35 mg QE/g respectively. Furthermore, phenolic compounds in leaves of Moringa were studied using a high-performance liquid chromatography coupled with mass spectrometry (HPLC-MS). Results showed that the largest number of phenolic compounds determined in leaves of Moringa belongs to flavonoids. Moreover, biological properties were determined by radical scavenging capacity (DPPH) and ferric-reducing power (FRAP). Cytoprotective effect and antioxidant power of Moringa extract were assessed using the mitochondrial activity testing method (MTT test), malondialdehyde (MDA), and reactive oxygen species (ROS) production. Results indicate that Moringa aqueous extract have a significant (i) proliferative, (ii) antioxidant, and (iii) cytoprotective effect on HCT116 and HEK293 cells against metal toxicity.

Alternative splicing is highly variable among Daphnia pulex lineages in response to acute copper exposure

BACKGROUND

Despite being one of the primary mechanisms of gene expression regulation in eukaryotes, alternative splicing is often overlooked in ecotoxicogenomic studies. The process of alternative splicing facilitates the production of multiple mRNA isoforms from a single gene thereby greatly increasing the diversity of the transcriptome and proteome. This process can be important in enabling the organism to cope with stressful conditions. Accurate identification of splice sites using RNA sequencing requires alignment to independent exonic positions within the genome, presenting bioinformatic challenges, particularly when using short read data. Although technological advances allow for the detection of splicing patterns on a genome-wide scale, very little is known about the extent of intraspecies variation in splicing patterns, particularly in response to environmental stressors. In this study, we used RNA-sequencing to study the molecular responses to acute copper exposure in three lineages of Daphnia pulex by focusing on the contribution of alternative splicing in addition to gene expression responses.

RESULTS

By comparing the overall gene expression and splicing patterns among all 15 copper-exposed samples and 6 controls, we identified 588 differentially expressed (DE) genes and 16 differentially spliced (DS) genes. Most of the DS genes (13) were not found to be DE, suggesting unique transcriptional regulation in response to copper that went unnoticed with conventional DE analysis. To understand the influence of genetic background on gene expression and alternative splicing responses to Cu, each of the three lineages was analyzed separately. In contrast to the overall analysis, each lineage had a higher proportion of unique DS genes than DE genes suggesting that genetic background has a larger influence on DS than on DE. Gene Ontology analysis revealed that some pathways involved in stress response were jointly regulated by DS and DE genes while others were regulated by only transcription or only splicing.

CONCLUSIONS

Our findings suggest an important role for alternative splicing in shaping transcriptome diversity in response to metal exposure in Daphnia, highlighting the importance of integrating splicing analyses with gene expression surveys to characterize molecular pathways in evolutionary and environmental studies.

Trophic-Level Interactive Effects of Phosphorus Availability on the Toxicities of Cadmium, Arsenic, and Their Binary Mixture in Media-Exposed Scenedesmus acutus and Media and Dietary-Exposed Daphnia pulex

Various anthropogenic activities simultaneously alter essential mineral nutrients and contaminant content in the environment. Depending on essential nutrient conditions, the uptake and effects of contaminants in exposed organisms may be altered. The addressing of ecological risk assessment (ERA) of contaminant mixtures has proven difficult. Furthermore, most assessments involving single contaminant exposures do not consider the interaction of essential nutrients on toxicological end points. Hypotheses for toxicological effects of cadmium (Cd), arsenic (As), and their binary mixture (Cd/As_mix) include alteration under varying dietary and media phosphorus (P) conditions. However, interactive effects and effect size (η²) are largely unknown. Here, we investigated the toxicities of Cd-, As-, and Cd/As_mix-treated media and diets on Scenedesmus acutus (a primary producer) and Daphnia pulex (a primary consumer), under varied media and dietary P conditions [low (LP), median (MP), and optimum (COMBO)]. Our results showed significant (p < 0.05) interactive effects and concentration dependent growth inhibition of S. acutus. The toxicity (at day 7) of Cd against S. acutus was 2×, 11×, and 4× that of As in LP, MP, and COMBO conditions, respectively, while the joint toxicity effects of Cd/As_mix were partially additive in LP and COMBO, and synergistic in MP media. Furthermore, acute lethal toxicity (96 h) of Cd in D. pulex was ∼60× that of As, while Cd/As_mix joint toxicity was synergistic. Chronic toxicity (14 d) in D. pulex showed significant (p < 0.05) interaction of As and P-availability on survival, reproduction, and behavior (distance moved, velocity, acceleration and mobility), while Cd and P availability showed significant interactive effect on rotational behavior. Dose response effects of Cd, As, and Cd/As_mix in S. acutus and D. pulex were either monophasic or biphasic under varying nutrient conditions. This study provides empirical evidence of the interactive effects of media/dietary P and toxic metals (Cd, As, and Cd/As_mix) at environmentally relevant concentrations, emphasizing the need for consideration of such interactions during ERA.

Development of the transcriptome for a sediment ecotoxicological model species, Chironomus dilutus

Chironomus dilutus is a prominent model species in conventional sediment toxicity testing and sediment contamination diagnosis. However, lack of genomic data significantly limited its application in identifying toxicological mode of action (MOA) and molecular biomarkers of toxicants. Here the transcriptome of C. dilutus in full life span and both sexes (1st, 2nd, 3rd and 4th instar larvae, pupae, and adults) were developed and temporal gene expression across adjacent life stages were investigated to understand the regulation of development. Furthermore, transcriptional response of Midges (the 4th instar larvae) exposed to chemicals of different MOAs (CdCl₂, nonylphenol and triclosan) were profiled based on the reference transcriptome. Consequently, a complete transcriptome of 31132 unigenes with N50 of 3117bp, covering 98.8% of the arthropod single-copy orthologs were assembled. While 364 genes were differentially expressed among adjacent larval stages, 7142 and 2127 of transcripts were significantly changed for the transition of larvae-pupae and pupae-adults, respectively. Finally, chemical-specific gene expression profile were identified in the midges, showed its potential in classifying distinct contaminants. Overall, the comprehensive transcriptome of C. dilutus developed here could not only facilitate the mechanistic understanding of environmental toxicants during critical life stage of aquatic insects, but also provide molecular diagnostic tools in sediment ecotoxicology.

Dynamics of Cadmium Acclimation in Daphnia pulex: Linking Fitness Costs, Cross-Tolerance, and Hyper-Induction of Metallothionein

Acclimation increases tolerance to stress in individuals but is assumed to contribute fitness costs when the stressor is absent, though data supporting this widely held claim are sparse. Therefore, using clonal (i.e., genetically identical) cultures of Daphnia pulex, we isolated the contributions of acclimation to the regulation of the metal response gene, metallothionein 1 (MT1), and defined the reproductive benefits and costs of cadmium (Cd)-acclimation. Daphnia pulex were exposed for 50 parthenogenetic generations to environmentally realistic levels (1 μg Cd/L), and tolerance to Cd and other metals assessed during this period via standard toxicity tests. These tests revealed (1) increased tolerance to Cd compared to genetically identical nonacclimated cultures, (2) fitness costs in Cd-acclimated Daphnia when Cd was removed, and (3) cross-tolerance of Cd-acclimated Daphnia to zinc and silver, but not arsenic, thereby defining a functional role for metallothionein. Indeed, Cd-acclimated clones had significantly higher expression of MT1 mRNA than nonacclimated clones, when Cd exposed. Both the enhanced induction of MT1 and tolerant phenotype were rapidly lost when Cd was removed (1-2 generations), which is further evidence of acclimation costs. These findings provide evidence for the widely held view that acclimation is costly and are important for investigating evolutionary principles of genetic assimilation and the survival mechanisms of natural populations that face changing environments.

Comparative Proteomics on Deep-Sea Amphipods after in Situ Copper Exposure

The interest in deep-sea mining increased along with the environmental concerns of these activities to the deep-sea fauna. The discovery of optimal biomarkers of deep-sea mining activities in deep-sea species is a crucial step toward the supply of important ecological information for environmental impact assessment. In this study, an in situ copper exposure experiment was performed on deep-sea scavenging amphipods. Abyssorchomene distinctus individuals were selected among all the exposed amphipods for molecular characterization. Copper concentration within the gut was assessed, followed by a tandem mass tag-based coupled with two-dimensional liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) applied to identify and quantify the protein expression changes after 48 h of exposure. 2937 proteins were identified and annotated, and 1918 proteins among all identified proteins were assigned by at least two nonambiguous peptides. The screening process was performed based on the differences in protein abundance and the specific correlation between the proteins and copper in previous studies. These differentially produced proteins include Na⁺/K⁺ ATPase, cuticle, chitinase, and proteins with unknown function. Their abundances showed correlation with copper and had high sensitivity to indicate the copper level, being here proposed as biomarker candidates for deep-sea mining activities in the future. This is a key step in the development of environmental impact assessment of deep-sea mining activities integrating ecotoxicological data.

Comparison between transcriptomic responses to short-term stress exposures of a common Holarctic and endemic Lake Baikal amphipods

Background

Lake Baikal is one of the oldest freshwater lakes and has constituted a stable environment for millions of years, in stark contrast to small, transient bodies of water in its immediate vicinity. A highly diverse endemic endemic amphipod fauna is found in one, but not the other habitat. We ask here whether differences in stress response can explain the immiscibility barrier between Lake Baikal and non-Baikal faunas. To this end, we conducted exposure experiments to increased temperature and the toxic heavy metal cadmium as stressors.

Results

Here we obtained high-quality de novo transcriptome assemblies, covering mutiple conditions, of three amphipod species, and compared their transcriptomic stress responses. Two of these species, Eulimnogammarus verrucosus and E. cyaneus, are endemic to Lake Baikal, while the Holarctic Gammarus lacustris is a potential invader.

Conclusions

Both Baikal species possess intact stress response systems and respond to elevated temperature with relatively similar changes in their expression profiles. G. lacustris reacts less strongly to the same stressors, possibly because its transcriptome is already perturbed by acclimation conditions.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6024-3) contains supplementary material, which is available to authorized users.

Rapid evolution in response to warming does not affect the toxicity of a pollutant: Insights from experimental evolution in heated mesocosms

While human-induced stressors such as warming and pollutants may co-occur and interact, evolutionary studies typically focus on single stressors. Rapid thermal evolution may help organisms better deal with warming, yet it remains an open question whether thermal evolution changes the toxicity of pollutants under warming. We investigated the effects of exposure to a novel pollutant (zinc oxide nanoparticles, nZnO) and 4°C warming (20°C vs. 24°C) on key life history and physiological traits of the water flea Daphnia magna, a keystone species in aquatic ecosystems. To address the role of thermal evolution, we compared these effects between clones from an experimental evolution trial where animals were kept for two years in outdoor mesocosms at ambient temperatures or ambient +4°C. The nZnO was more toxic at 20°C than at 24°C: only at 20°C, it caused reductions in early fecundity, intrinsic growth rate and metabolic activity. This was due to a higher accumulated zinc burden at 20°C than at 24°C, which was associated with an upregulation of a metallothionein gene at 20°C but not at 24°C. Clones from the heated mesocosms better dealt with warming than clones from the ambient mesocosms, indicating rapid thermal evolution. Notably, rapid thermal evolution did not change the toxicity of nZnO, neither at 20°C nor at 24°C, suggesting no pleiotropy or metabolic trade-offs were at work under the current experimental design. Evaluating whether thermal evolution influences the toxicity of pollutants is important for ecological risk assessment. It provides key information to extrapolate laboratory-derived toxicity estimates of pollutants both in space to warmer regions and in time under future global warming scenarios. In general, studying how the evolution of tolerance to one anthropogenic stressor influence tolerance to other anthropogenic stressors should get more attention in a rapidly changing world where animals increasingly face combinations of stressors.

Variation in transcriptional responses to copper exposure across Daphnia pulex lineages

Copper pollution is pervasive in aquatic habitats and is particularly harmful to invertebrates sensitive to environmental changes such as Daphnia pulex. Mechanisms of toxicity and tolerance to copper are not well understood. We used RNA-sequencing to investigate these mechanisms in three genetically distinct D. pulex clonal lineages with different histories of copper exposure. Upregulated genes after copper exposure were enriched with Gene Ontology (GO) categories involved in digestion, molting and growth, whereas downregulated genes after copper exposure were enriched in the metal-regulatory system, immune response and epigenetic modifications. The three D. pulex clones in our study show largely similar transcriptional patterns in response to copper, with only a total of twenty genes differentially expressed in a single clonal lineages. We also detected lower relative expression of some genes known to be important for copper tolerance, metallothionein and glutathione-S-transferase, in a sensitive lineage sampled from an uncontaminated habitat. Daphnia-specific genes (without orthologs outside the genus) and Daphnia-specific duplications (genes duplicated in the Daphnia lineage) were overrepresented in differentially expressed genes, highlighting an important role for newly emerged genes in tolerating environmental stressors. The results indicate that the D. pulex lineages tested in this study generally respond to copper stress using the same major pathways, but that the more resistant clone with previous copper exposure might be better able to regulate key genes. This finding highlights the important nuances in gene expression among clones, shaped by historical exposure and influencing copper tolerance.

Aberrant gene expression profiles in Mediterranean sea urchin reproductive tissues after metal exposures

Marine organisms are simultaneously exposed to numerous pollutants, among which metals probably represent the most abundant in marine environments. In order to evaluate the effects of metal exposure at molecular level in reproductive tissues, we profiled the sea urchin transcriptional response after non-lethal exposures using pathway-focused mRNA expression analyses. Herein, we show that exposures to relatively high concentrations of both essential and toxic metals hugely affected the gonadic expression of several genes involved in stress-response, detoxification, transcriptional and post-transcriptional regulation, without significant changes in gonadosomatic indices. Even though treatments did not result in reproductive tissues visible alterations, metal exposures negatively affected the main mechanisms of stress-response, detoxification and survival of adult P. lividus. Additionally, transcriptional changes observed in P. lividus gonads may cause altered gametogenesis and maintenance of heritable aberrant epigenetic effects. This study leads to the conclusion that exposures to metals, as usually occurs in polluted coastal areas, may affect sea urchin gametogenesis, thus supporting the hypothesis that parental exposure to environmental stressors affects the phenotype of the offspring.

Age-dependent survival, stress defense, and AMPK in Daphnia pulex after short-term exposure to a polystyrene nanoplastic

The widespread occurrence and accumulation of micro- and nanoplastics in aquatic environments has become a growing global concern. Generally, natural aquatic populations are characterized by a variety of multi-structured age groups, for which physiological and biochemical responses typically differ. The freshwater cladoceran, Daphnia pulex, is a model species used extensively in environmental monitoring studies and ecotoxicology testing. Here, the effects of a polystyrene nanoplastic on the physiological changes (i.e., survival) and expression levels of stress defense genes (i.e., those encoding antioxidant-mediated and heat shock proteins) in this freshwater flea were measured. Results from acute bioassays were used to determine the respective nanoplastic LC50 values for five age groups (1-, 4-, 7-, 14- and 21-day-old individuals): the obtained values for the 1- and 21-day-old D. pulex groups were similar (i.e., not significantly different). The expression levels of genes encoding key stress defense enzymes and proteins-SOD, CAT, GST, GPx, HSP70, and HSP90-were influenced by the nanoplastic in all the age groups, but not in the same way for each. Significant differences were observed among all age groups in their expression of the gene encoding the energy-sensing enzyme AMPK (adenosine monophosphate-activated protein kinase) α, β, and γ following exposure to the nanoplastic. Moreover, the expression of AMPK α was significantly increased in the 1-, 7-, and 21-day-old individuals exposed to nanoplastic relative to the control group. Together, these results indicate that age in D. pulex affects the sensitivity of its individuals to pollution from this nanoplastic, primarily via alterations to vital physiological and biochemical processes, such as cellular energy homeostasis and oxidation, which were demonstrated in vivo. We speculate that such age-related effects may extend to other nanoplastics and forms of pollution in D. pulex and perhaps similar marine organisms.

Identification of reference genes for RT-qPCR data normalization in Gammarus fossarum (Crustacea Amphipoda)

Gene expression profiling via RT-qPCR is a robust technique increasingly used in ecotoxicology. Determination and validation of optimal reference genes is a requirement for initiating RT-qPCR experiments. To our best knowledge, this study is the first attempt of identifying a set of reference genes for the freshwater crustacean Gammarus fossarum. Six candidate genes (Actin, TUB, UB, SDH, Clathrin and GAPDH) were tested in order to determine the most stable ones in different stress conditions and to increase the robustness of RT-qPCR data. SDH and Clathrin appeared as the most stable ones. A validation was performed using G. fossarum samples exposed for 15 days to AgNO₃, silver nanoparticles (AgNPs) 40 nm and gold nanoparticles (AuNPs) 40 nm. Effects on HSP90 were evaluated and data normalized using Clathrin and SDH. A down-regulation of HSP90 was observed when G. fossarum were exposed to AuNPs 40 nm whereas no effects were observed when G. fossarum were exposed to AgNPs 40 nm. This study highlights the importance of the preliminary determination of suitable reference genes for RT-qPCR experiments. Additionally, this study allowed, for the first time, the determination of a set of valuable genes that can be used in other RT-qPCR studies using G. fossarum as model organism.

Thermal evolution offsets the elevated toxicity of a contaminant under warming: A resurrection study in Daphnia magna

Synergistic interactions between temperature and contaminants are a major challenge for ecological risk assessment, especially under global warming. While thermal evolution may increase the ability to deal with warming, it is unknown whether it may also affect the ability to deal with the many contaminants that are more toxic at higher temperatures. We investigated how evolution of genetic adaptation to warming affected the interactions between warming and a novel stressor: zinc oxide nanoparticles (nZnO) in a natural population of Daphnia magna using resurrection ecology. We hatched resting eggs from two D. magna subpopulations (old: 1955-1965, recent: 1995-2005) from the sediment of a lake that experienced an increase in average temperature and in recurrence of heat waves but was never exposed to industrial waste. In the old "ancestral" subpopulation, exposure to a sublethal concentration of nZnO decreased the intrinsic growth rate, metabolic activity, and energy reserves at 24°C but not at 20°C, indicating a synergism between warming and nZnO. In contrast, these synergistic effects disappeared in the recent "derived" subpopulation that evolved a lower sensitivity to nZnO at 24°C, which indicates that thermal evolution could offset the elevated toxicity of nZnO under warming. This evolution of reduced sensitivity to nZnO under warming could not be explained by changes in the total internal zinc accumulation but was partially associated with the evolution of the expression of a key metal detoxification gene under warming. Our results suggest that the increased sensitivity to the sublethal concentration of nZnO under the predicted 4°C warming by the end of this century may be counteracted by thermal evolution in this D. magna population. Our results illustrate the importance of evolution to warming in shaping the responses to another anthropogenic stressor, here a contaminant. More general, genetic adaptation to an environmental stressor may ensure that synergistic effects between contaminants and this environmental stressor will not be present anymore.

Comparative Transcriptome Analysis for Understanding Predator-Induced Polyphenism in the Water Flea Daphnia pulex

The crustacean Daphnia pulex is one of the best model organisms for studying inducible defense mechanisms due to their inducible morphology in response to the predator Chaoborus larvae. In this study, multiple developmental stages of D. pulex were exposed to C. flavicans larvae and transcriptome profiles of samples from late embryo to fifth instar were sequenced by the RNA-seq technique to investigate the genetic background underlying inducible defenses. In comparison, differentially expressed genes between defensive and normal morphs were identified, including 908 genes in late embryo, 1383 genes in the first-third (1–3) instar, and 1042 genes in fourth-fifth (4–5) instar. Gene ontology enrichment analysis showed that structural constituents of the cuticle and structural molecule activity genes were prominent up-regulated genes in late embryos. Down-regulated genes in late embryos and 1–3 instar comprised metabolic process, hydrolase activity, and peptidase activity gene classes. Pathway analysis indicated that small molecule neurotransmitter pathways were potentially involved in the development of inducible defenses. The characterization of genes and pathways in multiple developmental stages can improve our understanding of inducible defense responses of D. pulex to predation at the molecular level.

Gene Coexpression Networks Drive and Predict Reproductive Effects in Daphnia in Response to Environmental Disturbances

Increasing effects of anthropogenic stressors and those of natural origin on aquatic ecosystems have intensified the need for predictive and functional models of their effects. Here, we use gene expression patterns in combination with weighted gene coexpression networks and generalized additive models to predict effects on reproduction in the aquatic microcrustacean Daphnia. We developed models to predict effects on reproduction upon exposure to different cyanobacteria, different insecticides and binary mixtures of cyanobacteria and insecticides. Models developed specifically for groups of stressors (e.g., either cyanobacteria or insecticides) performed better than general models developed on all data. Furthermore, models developed using in silico generated mixture gene expression profiles from single stressor data were able to better predict effects on reproduction compared to models derived from the mixture exposures themselves. Our results highlight the potential of gene expression data to quantify effects of complex exposures at higher level organismal effects without prior mechanistic knowledge or complex exposure data.

Genetic Diversity in Insect Metal Tolerance

Insects encounter a variety of metals in their environment, many of which are required at some concentration for normal organismal homeostasis, but essentially all of which are toxic at higher concentrations. Insects have evolved a variety of genetic, and likely epigenetic, mechanisms to deal with metal stress. A recurring theme in all these systems is complexity and diversity; even simple, single gene, cases are complex. Of the known gene families, the metallothioneins are perhaps the best understood and provide good examples of how diverse metal response is. Interestingly, there is considerable diversity across taxa in these metal-responsive systems, including duplications to form small gene families and complex expression of single loci. Strikingly, different species have evolved different mechanisms to cope with the same, or similar, stress suggesting both independent derivation of, and plasticity in, the pathways involved. It is likely that some metal-response systems evolved early in evolutionary time and have been conserved, while others have diverged, and still others evolved more recently and convergently. In addition to conventional genetics, insects likely respond to environmental metal through a variety of epigenetic systems, but direct tests are lacking. Ultimately, it is likely that classical genetic and epigenetic factors interact in regulating insect metal responses. In light of this diversity across species, future studies including a broad-based examination of gene expression in non-model species in complex environments will likely uncover additional genes and genetic and epigenetic mechanisms.

Four Transgenerational Demographic Performance of Moina macrocopa Exposed to Chronic Levels of Cadmium

In this study, we quantified intergenerational, demographic variability of Moina macrocopa subjected to cadmium stress. Exposure of M macrocopa to cadmium (0.2, 0.3, and 0.4 mg/L as CdCl₂) through 4 consecutive generations revealed changes in demographic responses not only in survivorship variables but also in reproductive parameters. Long-term demographic responses varied differently, depending on the demographic trait and the concentration of heavy metal in the medium. With the exception of generation time, all life history traits were significantly and adversely influenced due to increase in Cd concentrations. The average life span of M macrocopa varied up to 40% depending on Cd level and the generation of exposure. The highest gross reproductive rates were recorded in controls, while the lowest (∼30% less) were recorded at the highest Cd level. Survival-weighted net reproductive rates were reduced by nearly 50% due to Cd toxicity. The rate of population increase per day of M macrocopa was also significantly affected (∼25%) by Cd as compared to controls. This cladoceran showed a dose–response to Cd toxicity with a significance in both magnitude and frequency of offspring production.

Daphnia magna's sense of competition: intra-specific interactions (ISI) alter life history strategies and increase metals toxicity

This work investigates whether the scale-up to multi-animal exposures that is commonly applied in genomics studies provides equivalent toxicity outcomes to single-animal experiments of standard Daphnia magna toxicity assays. Specifically, we tested the null hypothesis that intraspecific interactions (ISI) among D. magna have neither effect on the life history strategies of this species, nor impact toxicological outcomes in exposure experiments with Cu and Pb. The results show that ISI significantly increased mortality of D. magna in both Cu and Pb exposure experiments, decreasing 14 day LC50 s and 95 % confidence intervals from 14.5 (10.9-148.3) to 8.4 (8.2-8.7) µg Cu/L and from 232 (156-4810) to 68 (63-73) µg Pb/L. Additionally, ISI potentiated Pb impacts on reproduction eliciting a nearly 10-fold decrease in the no-observed effect concentration (from 236 to 25 µg/L). As an indication of environmental relevance, the effects of ISI on both mortality and reproduction in Pb exposures were sustained at both high and low food rations. Furthermore, even with a single pair of Daphnia, ISI significantly increased (p < 0.05) neonate production in control conditions, demonstrating that ISI can affect life history strategy. Given these results we reject the null hypothesis and conclude that results from scale-up assays cannot be directly applied to observations from single-animal assessments in D. magna. We postulate that D. magna senses chemical signatures of conspecifics which elicits changes in life history strategies that ultimately increase susceptibility to metal toxicity.

Gene Family Evolution Reflects Adaptation to Soil Environmental Stressors in the Genome of the Collembolan Orchesella cincta

Collembola (springtails) are detritivorous hexapods that inhabit the soil and its litter layer. The ecology of the springtail Orchesella cincta is extensively studied in the context of adaptation to anthropogenically disturbed areas. Here, we present a draft genome of an O. cincta reference strain with an estimated size of 286.8 Mbp, containing 20,249 genes. In total, 446 gene families are expanded and 1,169 gene families evolved specific to this lineage. Besides these gene families involved in general biological processes, we observe gene clusters participating in xenobiotic biotransformation. Furthermore, we identified 253 cases of horizontal gene transfer (HGT). Although the largest percentage of them originated from bacteria (37.5%), we observe an unusually high percentage (30.4%) of such genes of fungal origin. The majority of foreign genes are involved in carbohydrate metabolism and cellulose degradation. Moreover, some foreign genes (e.g., bacillopeptidases) expanded after HGT. We hypothesize that horizontally transferred genes could be advantageous for food processing in a soil environment that is full of decaying organic material. Finally, we identified several lineage-specific genes, expanded gene families, and horizontally transferred genes, associated with altered gene expression as a consequence of genetic adaptation to metal stress. This suggests that these genome features may be preadaptations allowing natural selection to act on. In conclusion, this genome study provides a solid foundation for further analysis of evolutionary mechanisms of adaptation to environmental stressors.

De novo assembly of a tadpole shrimp (Triops newberryi) transcriptome and preliminary differential gene expression analysis

Next-generation sequencing techniques, such as RNA sequencing, have provided a wealth of genomic information for nonmodel species. Transcriptomic information can be used to quantify the patterns of gene expression, which can identify how environmental differences invoke organismal stress responses and provide a gauge in predicting species adaptability. In our study, we used RNA sequencing to characterize the first transcriptome from a naupliar tadpole shrimp (Triops newberryi) to identify the genes expressed during the early life history stages and which could be important for future genomic studies. RNA was extracted from naupliar T. newberryi that were reared in a laboratory-controlled setting and in two different water types, a native and a non-native condition. A total of six replicates, three per condition, were sequenced with the Illumina Hi-Seq 2000 achieving 365 M 50-nt reads. High-quality reads were produced and de novo assembly was used to construct a T. newberryi transcriptome that was approximately 24.8 M base pairs. More than 10 000 peptides were predicted from the assembly, and genes were sorted into gene ontology categories. The use of different water conditions allowed for a preliminary differential gene expression analysis in order to compare the changes in gene expression between conditions. There were 299 differentially expressed genes between water conditions that might serve as a focal point for future genomic studies of Triops acclimation to different environments. The Triops transcriptome could serve as vital genomic information for additional studies on Branchiopod crustaceans.

Effects of fluctuating temperature and food availability on reproduction and lifespan

Experimental studies on energetics and aging often remove two major factors that in part regulate the energy budget in a normal healthy individual: reproduction and fluctuating environmental conditions that challenge homeostasis. Here we use the cyclical parthenogenetic Daphnia pulex to evaluate the role of a fluctuating thermal environment on both reproduction and lifespan across six food concentrations. We test the hypotheses that (1) caloric restriction extends lifespan; (2) maximal reproduction will come with a cost of shortened lifespan; and (3) at a given food concentration, relative to a metabolically equivalent constant temperature environment a diel fluctuating thermal environment will alter the allocation of energy to reproduction and lifespan to maintain homeostasis. We did not identify a level of food concentration that extended lifespan in response to caloric restriction, and we found no cost of reproduction in terms of lifespan. Rather, the individuals at the highest food levels generally had the highest reproductive output and the longest lifespans, the individuals at the intermediate food level decreased reproduction and maintained lifespan, and the individuals at the three lower food concentrations had a decrease in reproduction and lifespan as would be predicted with increasing levels of starvation. Fluctuating temperature had no effect on lifespan at any food concentration, but delayed time to reproductive maturity and decreased early reproductive output at all food concentrations. This suggests that a fluctuating temperature regimen activates molecular pathways that alter energy allocation. The costs of fluctuating temperature on reproduction were not consistent across the lifespan. Statistical interactions for age of peak reproduction and lifetime fecundity suggest that senescence of the reproductive system may vary between temperature regimens at the different food concentrations.

Thermal variation and factors influencing vertical migration behavior in Daphnia populations

The antipredator behavior diel vertical migration (DVM), common in aquatic keystone species Daphnia, involves daily migration from warmer surface waters before dawn to cooler deeper waters after dusk. Plasticity in Daphnia DVM behavior optimizes fitness via trade-offs between growth, reproduction, and predator avoidance. Migration behavior is affected by co-varying biotic and abiotic factors, including light, predator cues, and anthropogenic stressors making it difficult to determine each factor's individual contribution to the variation in this behavior. This study aims to better understand this ecologically significant behavior in Daphnia by: (1) determining how Daphnia pulicaria thermal preferences vary within and among natural populations; (2) distinguishing the role of temperature verses depth in Daphnia vertical migration; and (3) defining how two anthropogenic stressors (copper and nickel) impact Daphnia migratory behavior. Simulated natural lake stratification were constructed in 8L (0.5m tall, 14.5cm wide) water columns to monitor under controlled laboratory conditions the individual effects of temperature gradients, depth, and metal stressors on Daphnia vertical migration. Three major findings are reported. First, while no difference in thermal preference was found among the four populations studied, within lake populations variability among isolates was high. Second, decoupling temperature and depth revealed that depth was a better predictor of Daphnia migratory patterns over temperature. Third, exposure to environmentally relevant concentrations of copper or nickel inhibited classic DVM behavior. These findings revealed the high variability in thermal preference found within Daphnia populations, elucidated the individual roles that depth and temperature have on migratory behavior, and showed how copper and nickel can interfere with the natural response of Daphnia to fish predator cues. Thus contributing to the body of knowledge necessary to predict how natural populations of Daphnia will be affected by climate related changes in lake temperatures and increased presence of anthropogenic stressors.

MicroRNAs are involved in cadmium tolerance in Daphnia pulex

Daphnia can develop tolerance to cadmium (Cd) after multi-generational exposures. Until now, Cd tolerance in this crustacean was thought to be mainly due to its sequestration via induction of metallothioneins (MTs). Our research supports other studies showing microRNAs (miRNAs) also play a role in this enhanced tolerance. We induced Cd tolerance in Daphnia pulex after exposing them for 25 generations and examined the maintenance of enhanced Cd tolerance under a Cd-free environment for an additional three generations. Acute Cd tolerance as well as long-term effects on population dynamics were measured in selected generations via 48h LC50 tests and 21 d reproductive tests, respectively. Cd tolerance was associated with differential expression of 10 miRNAs (miR-2, miR-33, miR-92, miR-96, miR-153, miR-252, miR-279, miR-283, miR-305 and miR-615). Pathway analysis revealed these miRNAs might increase Cd tolerance by suppressing cellular growth and proliferation by GTPase and cuticle protein pathways, which switch cellular energy allocation to detoxification processes. Moreover, we found increased Cd tolerance is related with induction of MT3 and MT4 and a subsequent downregulation of MT1 and MT3 expression when animals are moved to a Cd-free environment. This is the first study linking aquatic invertebrate miRNAs with induced tolerance to environmental stressors.

Gene Body Methylation Patterns in Daphnia Are Associated with Gene Family Size

The relation between gene body methylation and gene function remains elusive. Yet, our understanding of this relationship can contribute significant knowledge on how and why organisms target specific gene bodies for methylation. Here, we studied gene body methylation patterns in two Daphnia species. We observed both highly methylated genes and genes devoid of methylation in a background of low global methylation levels. A small but highly significant number of genes was highly methylated in both species. Remarkably, functional analyses indicate that variation in methylation within and between Daphnia species is primarily targeted to small gene families whereas large gene families tend to lack variation. The degree of sequence similarity could not explain the observed pattern. Furthermore, a significant negative correlation between gene family size and the degree of methylation suggests that gene body methylation may help regulate gene family expansion and functional diversification of gene families leading to phenotypic variation.

The effect of a low iron diet and early life methylmercury exposure in Daphnia pulex

Iron (Fe) deficiency increases risk for adverse health outcomes in humans; however little is known about the potential interaction with methylmercury (MeHg) exposure. Studies testing multiple stressor hypotheses are expensive and time consuming in mammalian model systems; therefore, determining relevance of alternative models is important. Daphnia pulex were fed standard or low-Fe diets of freshwater algae, Ankistrodesmus falcatus. MeHgCl (1600 ng/L) or vehicle was added to culture media for 24 h during early life, and the combinatorial effects of a low-Fe diet and MeHg exposure on lifespan, maturation time, and reproduction were evaluated. Lipid storage effects were measured using image analysis of Oil Red O staining and triacylglyceride quantification. Our results show a dose-dependent reduction in lifespan in D. pulex fed low Fe diets. Lipid analysis suggests an interactive effect of diet and MeHg exposure, with MeHg exposure increasing lipid storage in D. pulex fed a low-Fe diet. These findings suggest the effects of dietary iron intake and early life MeHg exposure in D. pulex may be mediated by changes in energetics that result in differential lipid storage. Therefore, lipid storage in D. pulex may be a useful screen for detecting long-term effects of multiple stressors early in life.

Metallothionein and Hsp70 trade-off against one another in Daphnia magna cross-tolerance to cadmium and heat stress

The association between the insensitivity of adapted ecotypes of invertebrates to environmental stress, such as heavy metal pollution, and overall low Hsp levels characterizing these organisms has been attracting attention in various studies. The present study seeks to induce and examine this phenomenon in Daphnia magna by multigenerational acclimation to cadmium in a controlled laboratory setting. In this experiment, interclonal variation was examined: two clones of D. magna that have previously been characterized to diverge regarding their cadmium resistance and levels of the stress protein Hsp70, were continuously exposed to a sublethal concentration of Cd over four generations to study the effects of acclimation on Hsp70, metallothionein (MT), reproduction and cross-tolerance to heat stress. The two clones differed in all the measured parameters in a characteristic way, clone T displaying Cd and heat resistance, lower Hsp70 levels and offspring numbers on the one hand and higher MT expression on the other hand, clone S the opposite for all these parameters. We observed only slight acclimation-induced changes in constitutive Hsp70 levels and reproductive output. The differences in MT expression between clones as well as between acclimated organisms and controls give evidence for MT accounting for the higher Cd tolerance of clone T. Overall high Hsp70 levels of clone S did not confer cross tolerance to heat stress, contrary to common expectations. Our results suggest a trade-off between the efforts to limit the proteotoxic symptoms of Cd toxicity by Hsp70 induction and those to sequester and detoxify Cd by means of MT.

Novel Cadmium Responsive MicroRNAs in Daphnia pulex

Daphnia pulex is a widely used toxicological model and is known for its sensitivity to cadmium (Cd). Recent research suggests that microRNAs (miRNAs) play a critical role in animal responses to heavy metals. To investigate the functions of D. pulex miRNAs under Cd exposure, we analyzed the miRNA profiles of D. pulex after 48 h using miRNA microarrays and validated our findings by q-PCR. miRNA dpu-let-7 was identified as a stably expressed gene and used as a reference. We identified 22 and 21 differentially expressed miRNAs under low (20 μg/L CdCl2) and high-exposure (40 μg/L CdCl2) concentrations compared to controls, respectively. Cellular functions of predicted miRNA target Cd-responsive genes included oxidative stress, ion transport, mitochondrial damage, and DNA repair. An insulin-related network was also identified in relation to several Cd-responsive miRNAs. The expression of three predicted target genes for miR-71 and miR-210 were evaluated, and expression of two of them (SCN2A and SLC31A1) was negatively correlated with the expression of their regulator miRNAs. We show miR-210 is hypoxia-responsive in D. pulex and propose Cd and hypoxia induce miR-210 via a same HIF1α modulated pathway. Collectively, this research advances our understanding on the role of miRNAs in response to heavy-metal exposure.

Hybridization and the Origin of Contagious Asexuality in Daphnia pulex

Hybridization plays a potentially important role in the origin of obligate parthenogenesis (OP) in many organisms. However, it remains controversial whether hybridization directly triggers the transition from sexual reproduction to obligate asexuality or a hybrid genetic background enables asexual species to persist. Furthermore, we know little about the specific genetic elements from the divergent, yet still hybridizing lineages responsible for this transition and how these elements are further spread to create other OP lineages. In this study, we address these questions in Daphnia pulex, where cyclically parthenogenetic (CP) and OP lineages coexist. Ancestry estimates and whole-genome association mapping using 32 OP isolates suggest that a complex hybridization history between the parental species D. pulex and D. pulicaria is responsible for the introgression of a set of 647 D. pulicaria single nucleotide polymorphism alleles that show perfect association with OP. Crossing experiments using males of OP lineages and females of CP lineages strongly support a polygenic basis for OP. Single-sperm analyses show that although normal meiotic recombination occurs in the production of haploid sperm by males of OP lineages, a significant proportion of such sperm are polyploid, suggesting that the spread of asexual elements through these males (i.e., contagious asexuality) is much less efficient than previously envisioned. Although the current Daphnia genome annotation does not provide mechanistic insight into the nature of the asexuality-associated alleles, these alleles should be considered as candidates for future investigations on the genetic underpinnings of OP.

Impact of Predator Cues on Responses to Silver Nanoparticles in Daphnia carinata

The past decades have witnessed a boom in nanotechnology that has led to increasing production and application of silver nanoparticles (AgNPs) in the textile industry due to their antimicrobial properties. Increase in the manufacture and use of NPs inevitably has resulted in their increased release into aquatic environments resulting in the exposure of organisms living in these environments. Recently, the risk of exposure to NPs and the potential interaction with biological systems has received increasing attention. The present study investigated the potential effects of predator cues on the toxicity of environmentally relevant concentrations of AgNPs in Daphnia carinata at organismal and biochemical levels. The results of this study show that exposure to environmentally relevant concentrations of AgNPs can result in adverse effects on daphnids with 24- and 48-h LC50 values of 3.56 and 1.75 μg/L, respectively. Furthermore, significant inhibition of reproduction was observed at concentrations as low as 0.5 μg/L. Exposure to predator cues alone resulted in an increase in reproduction and inhibition of superoxide dismutase activity in daphnids. However, coexposure to predator cues interacted in an antagonistic manner with AgNPs with a 24-h LC50 value of 10.81 μg/L compared with 3.56 μg/L for AgNPs alone. In summary, AgNPs could pose risks to aquatic invertebrates at environmentally relevant concentrations. Interestingly, the presence of other factors, such as predator cues, moderated the effects of exposure to AgNPs. Therefore, there is a need to further investigate the potential interactions between NPs and biological factors that can modulate toxicity of NPs for application to the risk assessment of aquatic invertebrates.

Relating suborganismal processes to ecotoxicological and population level endpoints using a bioenergetic model

Ecological effects of environmental stressors are commonly evaluated using organismal or suborganismal data, such as standardized toxicity tests that characterize responses of individuals (e.g., mortality and reproduction) and a rapidly growing body of "omics" data. A key challenge for environmental risk assessment is relating such information to population dynamics. One approach uses dynamic energy budget (DEB) models that relate growth and reproduction of individuals to underlying flows of energy and elemental matter. We hypothesize that suborganismal information identifies DEB parameters that are most likely impacted by a particular stressor and that the DEB model can then project suborganismal effects on life history and population endpoints. We formulate and parameterize a model of growth and reproduction for the water flea Daphnia magna. Our model resembles previous generic bioenergetic models, but has explicit representation of discrete molts, an important feature of Daphnia life history. We test its ability to predict six endpoints commonly used in chronic toxicity studies in specified food environments. With just one adjustable parameter, the model successfully predicts growth and reproduction of individuals from a wide array of experiments performed in multiple laboratories using different clones of D. magna raised on different food sources. Fecundity is the most sensitive endpoint, and there is broad correlation between the sensitivities of fecundity and long-run growth rate, as is desirable for the default metric used in chronic toxicity tests. Under some assumptions, we can combine our DEB model with the Euler-Lotka equation to estimate longrun population growth rates at different food levels. A review of Daphnia gene-expression experiments on the effects of contaminant exposure reveals several connections to model parameters, in particular a general trend of increased transcript expression of genes involved in energy assimilation and utilization at concentrations affecting growth and reproduction. The sensitivity of fecundity to many model parameters was consistent with frequent generalized observations of decreased expression of genes involved in reproductive physiology, but interpretation of these observations requires further mechanistic modeling. We thus propose an approach based on generic DEB models incorporating few essential species-specific features for rapid extrapolation of ecotoxicogenomic assays for Daphnia-based population risk assessment.

Global cytosine methylation in Daphnia magna depends on genotype, environment, and their interaction

The authors characterized global cytosine methylation levels in 2 different genotypes of the ecotoxicological model organism Daphnia magna after exposure to a wide array of biotic and abiotic environmental stressors. The present study aimed to improve the authors' understanding of the role of cytosine methylation in the organism's response to environmental conditions. The authors observed a significant genotype effect, an environment effect, and a genotype × environment effect. In particular, global cytosine methylation levels were significantly altered after exposure to Triops predation cues, Microcystis, and sodium chloride compared with control conditions. Significant differences between the 2 genotypes were observed when animals were exposed to Triops predation cues, Microcystis, Cryptomonas, and sodium chloride. Despite the low global methylation rate under control conditions (0.49-0.52%), global cytosine methylation levels upon exposure to Triops demonstrated a 5-fold difference between the genotypes (0.21% vs 1.02%). No effects were found in response to arsenic, cadmium, fish, lead, pH of 5.5, pH of 8, temperature, hypoxia, and white fat cell disease. The authors' results point to the potential role of epigenetic effects under changing environmental conditions such as predation (i.e., Triops), diet (i.e., Cryptomonas and Microcystis), and salinity. The results of the present study indicate that, despite global cytosine methylation levels being low, epigenetic effects may be important in environmental studies on Daphnia.

Conserved transcriptional responses to cyanobacterial stressors are mediated by alternate regulation of paralogous genes in Daphnia

Despite a significant increase in genomic data, our knowledge of gene functions and their transcriptional responses to environmental stimuli remains limited. Here, we use the model keystone species Daphnia pulex to study environmental responses of genes in the context of their gene family history to better understand the relationship between genome structure and gene function in response to environmental stimuli. Daphnia were exposed to five different treatments, each consisting of a diet supplemented with one of five cyanobacterial species, and a control treatment consisting of a diet of only green algae. Differential gene expression profiles of Daphnia exposed to each of these five cyanobacterial species showed that genes with known functions are more likely to be shared by different expression profiles, whereas genes specific to the lineage of Daphnia are more likely to be unique to a given expression profile. Furthermore, while only a small number of nonlineage-specific genes were conserved across treatment type, there was a high degree of overlap in expression profiles at the functional level. The conservation of functional responses across the different cyanobacterial treatments can be attributed to the treatment-specific expression of different paralogous genes within the same gene family. Comparison with available gene expression data in the literature suggests differences in nutritional composition in diets with cyanobacterial species compared to diets of green algae as a primary driver for cyanobacterial effects on Daphnia. We conclude that conserved functional responses in Daphnia across different cyanobacterial treatments are mediated through alternate regulation of paralogous gene families.

Synergistic interactions of biotic and abiotic environmental stressors on gene expression

Understanding the response of organisms to multiple stressors is critical for predicting if populations can adapt to rapid environmental change. Natural and anthropogenic stressors often interact, complicating general predictions. In this study, we examined the interactive and cumulative effects of two common environmental stressors, lowered calcium concentration, an anthropogenic stressor, and predator presence, a natural stressor, on the water flea Daphnia pulex. We analyzed expression changes of five genes involved in calcium homeostasis — cuticle proteins (Cutie, Icp2), calbindin (Calb), and calcium pump and channel (Serca and Ip3R) — using real-time quantitative PCR (RT-qPCR) in a full factorial experiment. We observed strong synergistic interactions between low calcium concentration and predator presence. While the Ip3R gene was not affected by the stressors, the other four genes were affected in their transcriptional levels by the combination of the stressors. Transcriptional patterns of genes that code for cuticle proteins (Cutie and Icp2) and a sarcoplasmic calcium pump (Serca) only responded to the combination of stressors, changing their relative expression levels in a synergistic response, while a calcium-binding protein (Calb) responded to low calcium stress and the combination of both stressors. The expression pattern of these genes (Cutie, Icp2, and Serca) were nonlinear, yet they were dose dependent across the calcium gradient. Multiple stressors can have complex, often unexpected effects on ecosystems. This study demonstrates that the dominant interaction for the set of tested genes appears to be synergism. We argue that gene expression patterns can be used to understand and predict the type of interaction expected when organisms are exposed simultaneously to natural and anthropogenic stressors.

Distinct expression profiles of stress defense and DNA repair genes in Daphnia pulex exposed to cadmium, zinc, and quantum dots

The ever-increasing production and use of nanocrystaline semiconductors (Quantum dots; QDs) will inevitably result in increased appearance of these nanomaterials in the aquatic environment. However, the behavior and potential toxicity of heavy metal constituted nanoparticulates in aquatic invertebrates is largely unknown, especially with regard to molecular responses. The freshwater crustacean Daphnia pulex is a well-suited toxicological and ecological model to study molecular responses to environmental stressors. In this study, D. pulex were exposed for 48 h to sublethal doses of QDs (25% and 50% of LC50) with differing spectral properties (CdTe and CdSe/ZnS QDs) and Cd and Zn salts. Our data suggest that acute exposure to both CdSO4 and Cd-based QDs leads to Cd uptake in vivo, which was biologically supported by the observation of increased expression of metallothionein (MT-1). Furthermore, Cd, Zn, and CdSe/ZnS QDs induced different patterns of gene expression regarding stress defense and DNA repair, which furthers our knowledge regarding which response pathways are affected by nanoparticulate forms of metals versus ionic forms in aquatic crustaceans.

Effects of early life exposure to methylmercury in Daphnia pulex on standard and reduced food ration

As a well-known eco-toxicological model organism, Daphnia pulex may also offer advantages in human health research for assessing long-term effects of early life exposures to coupled stressors. Here, we examine consequences of early life exposure to methylmercury (MeHg) under standard and reduced food ration. We exposed Daphnia for 24h in early life to varying concentrations of methylmercury(II) chloride (0, 200, 400, 800 and 1600ng/L) and thereafter kept Daphnia on either a standard or a reduced food ration. The data suggests an additive effect of MeHg concentration and food ration on decreasing lifespan, although MeHg concentration does not affect survival linearly. Food ration and MeHg concentration were predictive of reduced reproduction, and there is some evidence of an interaction (p=0.048). Multi-stressor work in alternative model systems may be useful for prioritizing research, taking into account potential antagonistic, additive or synergistic effects that nutritional status may have on chemical toxicity.