Jana Asselman

Identification of Pathways, Gene Networks, and Paralogous Gene Families in Daphnia pulex Responding to Exposure to the Toxic Cyanobacterium Microcystis aeruginosa

Although cyanobacteria produce a wide range of natural toxins that impact aquatic organisms, food webs, and water quality, the mechanisms of toxicity are still insufficiently understood. Here, we implemented a whole-genome expression microarray to identify pathways, gene networks, and paralogous gene families responsive to Microcystis stress in Daphnia pulex . Therefore, neonates of a sensitive isolate were given a diet contaminated with Microcystis to contrast with those given a control diet for 16 days. The microarray revealed 2247 differentially expressed (DE) genes (7.6% of the array) in response to Microcystis , of which 17% are lineage-specific (i.e., these genes have no detectable homology to any other gene in currently available databases) and 49% are gene duplicates (paralogues). We identified four pathways/gene networks and eight paralogous gene families affected by Microcystis . Differential regulation of the ribosome, including three paralogous gene families encoding 40S, 60S, and mitochondrial ribosomal proteins, suggests an impact of Microcystis on protein synthesis of D. pulex . In addition, differential regulation of the oxidative phosphorylation pathway (including the NADH:ubquinone oxidoreductase gene family) and the trypsin paralogous gene family (a major component of the digestive system in D. pulex ) could explain why fitness is reduced based on energy budget considerations.

Daphnia magna transcriptome by RNA-Seq across 12 environmental stressors

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

Genome-wide transcription profiles reveal genotype-dependent responses of biological pathways and gene-families in Daphnia exposed to single and mixed stressors.

The present study investigated the possibilities and limitations of implementing a genome-wide transcription-based approach that takes into account genetic and environmental variation to better understand the response of natural populations to stressors. When exposing two different Daphnia pulex genotypes (a cadmium-sensitive and a cadmium-tolerant one) to cadmium, the toxic cyanobacteria Microcystis aeruginosa, and their mixture, we found that observations at the transcriptomic level do not always explain observations at a higher level (growth, reproduction). For example, although cadmium elicited an adverse effect at the organismal level, almost no genes were differentially expressed after cadmium exposure. In addition, we identified oxidative stress and polyunsaturated fatty acid metabolism-related pathways, as well as trypsin and neurexin IV gene-families as candidates for the underlying causes of genotypic differences in tolerance to Microcystis. Furthermore, the whole-genome transcriptomic data of a stressor mixture allowed a better understanding of mixture responses by evaluating interactions between two stressors at the gene-expression level against the independent action baseline model. This approach has indicated that ubiquinone pathway and the MAPK serine-threonine protein kinase and collagens gene-families were enriched with genes showing an interactive effect in expression response to exposure to the mixture of the stressors, while transcription and translation-related pathways and gene-families were mostly related with genotypic differences in interactive responses to this mixture. Collectively, our results indicate that the methods we employed may improve further characterization of the possibilities and limitations of transcriptomics approaches in the adverse outcome pathway framework and in predictions of multistressor effects on natural populations.

A new reference genome assembly for the microcrustacean Daphnia pulex

Comparing genomes of closely related genotypes from populations with distinct demographic histories can help reveal the impact of effective population size on genome evolution. For this purpose, we present a high quality genome assembly of Daphnia pulex (PA42), and compare this with the first sequenced genome of this species (TCO), which was derived from an isolate from a population with >90% reduction in nucleotide diversity. PA42 has numerous similarities to TCO at the gene level, with an average amino acid sequence identity of 98.8 and >60% of orthologous proteins identical. Nonetheless, there is a highly elevated number of genes in the TCO genome annotation, with ∼7000 excess genes appearing to be false positives. This view is supported by the high GC content, lack of introns, and short length of these suspicious gene annotations. Consistent with the view that reduced effective population size can facilitate the accumulation of slightly deleterious genomic features, we observe more proliferation of transposable elements (TEs) and a higher frequency of gained introns in the TCO genome.

Transcription patterns of genes encoding four metallothionein homologs in Daphnia pulex exposed to copper and cadmium are time- and homolog-dependent

Metallothioneins are proteins that play an essential role in metal homeostasis and detoxification in nearly all organisms studied to date. Yet discrepancies between outcomes of chronic and acute exposure experiments hamper the understanding of the regulatory mechanisms of their isoforms following metal exposure. Here, we investigated transcriptional differences among four identified homologs (mt1-mt4) in Daphnia pulex exposed across time to copper and cadmium relative to a control. Transcriptional upregulation of mt1 and mt3 was detected on day four following exposure to cadmium, whereas that of mt2 and mt4 was detected on day two and day eight following exposure to copper. These results confirm temporal and metal-specific differences in the transcriptional induction of genes encoding metallothionein homologs upon metal exposure which should be considered in ecotoxicological monitoring programs of metal-contaminated water bodies. Indeed, the mRNA expression patterns observed here illustrate the complex regulatory system associated with metallothioneins, as these patterns are not only dependent on the metal, but also on exposure time and the homolog studied. Further phylogenetic analysis and analysis of regulatory elements in upstream promoter regions revealed a high degree of similarity between metallothionein genes of Daphnia pulex and Daphnia magna, a species belonging to the same genus. These findings, combined with a limited amount of available expression data for D. magna metallothionein genes, tentatively suggest a potential generalization of the metallothionein response system between these Daphnia species.

Bisulfite Sequencing with Daphnia Highlights a Role for Epigenetics in Regulating Stress Response to Microcystis through Preferential Differential Methylation of Serine and Threonine Amino Acids

Little is known about the influence that environmental stressors may have on genome-wide methylation patterns, and to what extent epigenetics may be involved in environmental stress response. Yet, studies of methylation patterns under stress could provide crucial insights on stress response and toxicity pathways. Here, we focus on genome-wide methylation patterns in the microcrustacean Daphnia magna, a model organism in ecotoxicology and risk assessment, exposed to the toxic cyanobacterium Microcystis aeruginosa. Bisulfite sequencing of exposed and control animals highlighted differential methylation patterns in Daphnia upon exposure to Microcystis primarily in exonic regions. These patterns are enriched for serine/threonine amino acid codons and genes related to protein synthesis, transport and degradation. Furthermore, we observed that genes with differential methylation corresponded well with genes susceptible to alternative splicing in response to Microcystis stress. Overall, our results suggest a complex mechanistic response in Daphnia characterized by interactions between DNA methylation and gene regulation mechanisms. These results underscore that DNA methylation is modulated by environmental stress and can also be an integral part of the toxicity response in our study species.

Combined exposure to cyanobacteria and carbaryl results in antagonistic effects on the reproduction of daphnia pulex

In aquatic ecosystems, Daphnia are exposed to a wide variety of natural and chemical stressors that can cause interactive effects resulting in an increased impact on aquatic ecosystems. The authors therefore investigated the interactive effects of harmful cyanobacteria (cyanoHABs) with carbaryl in Daphnia pulex, because cyanobacteria have become an important concern for aquatic ecosystems. Daphnia were exposed for 21 d to 4 selected cyanobacteria (Aphanizomenon sp., Cylindrospermopsis raciborskii, Microcystis aeruginosa, and Oscillatoria sp), carbaryl, and all binary combinations of carbaryl and each individual cyanobacterium. Results were analyzed with both the independent action and the concentration addition model. The estimated median effect concentration (EC50) for carbaryl was comparable between the experiments, ranging from 2.28 µg/L to 5.94 µg/L. The EC50 for cyanobacteria ranged from 13.45% of the total diet ratio for Microcystis to 66.69% of the diet ratio for Oscillatoria. In combination with carbaryl, the response of Daphnia to 3 of the 4 cyanobacteria demonstrated antagonistic deviation patterns (p < 0.05). Exposure to combinations of carbaryl and Cylindrospermopsis did not result in statistically significant deviations from both reference models. The results may have important implications for pesticide risk assessment, underlining the impact of interactive effects on aquatic organisms. Based on these results, the authors suggest that both the independent action and the concentration addition model can serve as a protective approach in risk assessment of carbaryl in the presence of cyanobacterial blooms.

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.

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.

Will genetic adaptation of natural populations to chemical pollution result in lower or higher tolerance to future climate change

Evolved co-tolerance or cross-tolerance to chemical pollution and a future climate change stressor is observed