Daphnia comes of age: an ecological model in the genomic era

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

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

Supplementary Information

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

Intraspecific Variation in the Rates of Mutations Causing Structural Variation in Daphnia magna

Mutations that cause structural variation are important sources of genetic variation upon which other evolutionary forces can act, however, they are difficult to observe and therefore few direct estimates of their rate and spectrum are available. Understanding mutation rate evolution, however, requires adding to the limited number of species for which direct estimates are available, quantifying levels of intraspecific variation in mutation rates, and assessing whether rate estimates co-vary across types of mutation. Here, we report structural variation-causing mutation rates (svcMRs) for six categories of mutations (short insertions and deletions, long deletions and duplications, and deletions and duplications at copy number variable sites) from nine genotypes of Daphnia magna collected from three populations in Finland, Germany, and Israel using a mutation accumulation approach. Based on whole-genome sequence data and validated using simulations, we find svcMRs are high (two orders of magnitude higher than base substitution mutation rates measured in the same lineages), highly variable among populations, and uncorrelated across categories of mutation. Furthermore, to assess the impact of scvMRs on the genome, we calculated rates while adjusting for the lengths of events and ran simulations to determine if the mutations occur in genic regions more or less frequently than expected by chance. Our results pose a challenge to most prevailing theories aimed at explaining the evolution of the mutation rate, underscoring the importance of obtaining additional mutation rate estimates in more genotypes, for more types of mutation, in more species, in order to improve our future understanding of mutation rates, their variation, and their evolution.

Engines of change: Transposable element mutation rates are high and variable within Daphnia magna

Transposable elements (TEs) represent a major portion of most eukaryotic genomes, yet little is known about their mutation rates or how their activity is shaped by other evolutionary forces. Here, we compare short- and long-term patterns of genome-wide mutation accumulation (MA) of TEs among 9 genotypes from three populations of Daphnia magna from across a latitudinal gradient. While the overall proportion of the genome comprised of TEs is highly similar among genotypes from Finland, Germany, and Israel, populations are distinguishable based on patterns of insertion site polymorphism. Our direct rate estimates indicate TE movement is highly variable (net rates ranging from -11.98 to 12.79 x 10-5 per copy per generation among genotypes), differing both among populations and TE families. Although gains outnumber losses when selection is minimized, both types of events appear to be highly deleterious based on their low frequency in control lines where propagation is not limited to random, single-progeny descent. With rate estimates 4 orders of magnitude higher than base substitutions, TEs clearly represent a highly mutagenic force in the genome. Quantifying patterns of intra- and interspecific variation in TE mobility with and without selection provides insight into a powerful mechanism generating genetic variation in the genome.

Calmodulin inhibition as a mode of action of antifungal imidazole pharmaceuticals in non-target organisms

To improve assessment of risks associated with pharmaceutical contamination of the environment, it is crucial to understand effects and mode of action of drugs in non-target species. The evidence is accumulating that species with well-conserved drug targets are prone to be at risk when exposed to pharmaceuticals. An interesting group of pharmaceuticals released into the environment is imidazoles, antifungal agents with inhibition of ergosterol synthesis as a primary mode of action in fungi. However, imidazoles have also been identified as competitive antagonists of calmodulin (CaM), a calcium-binding protein with phylogenetically conserved structure and function. Therefore, imidazoles would act as CaM inhibitors in various organisms, including those with limited capacity to synthesize sterols, such as arthropods. We hypothesized that effects observed in crustaceans exposed to imidazoles are related to the CaM inhibition and CaM-dependent nitric oxide (NO) synthesis. To test this hypothesis, we measured (i) CaM levels and its gene expression, (ii) NO accumulation and (iii) gene expression of NO synthase (NOS1 and NOS2), in the cladoceran Daphnia magna exposed to miconazole, a model imidazole drug. Whereas significantly increased CaM gene expression and its cellular allocation were observed, supporting the hypothesized mode of action, no changes occurred in either NO synthase expression or NO levels in the exposed animals. These findings suggest that CaM inhibition by miconazole leads to protein overexpression that compensates for the loss in the protein activity, with no measurable downstream effects on NO pathways. The inhibition of CaM in D. magna may have implications for effect assessment of exposure to mixtures of imidazoles in aquatic non-target species.

Variation in the Microbiota Associated with Daphnia magna Across Genotypes, Populations, and Temperature

Studies of how the microbiome varies among individuals, populations, and abiotic conditions are critical for understanding this key component of an organism's biology and ecology. In the case of Daphnia, aquatic microcrustaceans widely used in population/community ecology and environmental science studies, understanding factors that influence microbiome shifts among individuals is useful for both basic and applied research contexts. In this study, we assess differences in the microbiome among genotypes of D. magna collected from three regions along a large latitudinal gradient (Finland, Germany, and Israel). After being reared in the lab for many years, we sought to characterize any differences in genotype- or population-specific microbial communities, and to assess whether the microbiota varied among temperatures. Our study is similar to a recent comparison of the microbial communities among D. magna genotypes raised in different temperatures published by Sullam et al. (Microb Ecol 76(2):506-517, 2017), and as such represents one of the first examples of a reproducible result in microbiome research. Like the previous study, we find evidence for a strong effect of temperature on the microbiome of D. magna, although across a much smaller temperature range representing potential near-future climates. In addition, we find evidence that the microbiomes of D. magna genotypes from different regions are distinct, even years after being brought into the laboratory. Finally, our results highlight a potentially common finding in the expanding area of microbiome research-differences among treatments are not necessarily observed in the most abundant taxonomic groups. This highlights the importance of considering sampling scheme and depth of coverage when characterizing the microbiome, as different experimental designs can significantly impact taxon-specific results, even when large-scale effects are reproduced.

Intraspecific Variation in Microsatellite Mutation Profiles in Daphnia magna

Microsatellite loci (tandem repeats of short nucleotide motifs) are highly abundant in eukaryotic genomes and often used as genetic markers because they can exhibit variation both within and between populations. Although widely recognized for their mutability and utility, the mutation rates of microsatellites have only been empirically estimated in a few species, and have rarely been compared across genotypes and populations within a species. Here, we investigate the dynamics of microsatellite mutation over long- and short-time periods by quantifying the starting abundance and mutation rates for microsatellites for six different genotypes of Daphnia magna, an aquatic microcrustacean, collected from three populations (Finland, Germany, and Israel). Using whole-genome sequences of these six starting genotypes, descendent mutation accumulation (MA) lines, and large population controls (non-MA lines), we find each genotype exhibits a distinctive initial microsatellite profile which clusters according to the population-of-origin. During the period of MA, we observe motif-specific, highly variable, and rapid microsatellite mutation rates across genotypes of D. magna, the average of which is order of magnitude greater than the recently reported rate observed in a single genotype of the congener, Daphnia pulex. In our experiment, genotypes with more microsatellites starting out exhibit greater losses and those with fewer microsatellites starting out exhibit greater gains—a context-dependent mutation bias that has not been reported previously. We discuss how genotype-specific mutation rates and spectra, in conjunction with evolutionary forces, can shape both the differential accumulation of repeat content in the genome and the evolution of mutation rates.

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.

Transcriptional profiling of predator-induced phenotypic plasticity in Daphnia pulex

Background

Predator-induced defences are a prominent example of phenotypic plasticity found from single-celled organisms to vertebrates. The water flea Daphnia pulex is a very convenient ecological genomic model for studying predator-induced defences as it exhibits substantial morphological changes under predation risk. Most importantly, however, genetically identical clones can be transcriptionally profiled under both control and predation risk conditions and be compared due to the availability of the sequenced reference genome. Earlier gene expression analyses of candidate genes as well as a tiled genomic microarray expression experiment have provided insights into some genes involved in predator-induced phenotypic plasticity. Here we performed the first RNA-Seq analysis to identify genes that were differentially expressed in defended vs. undefended D. pulex specimens in order to explore the genetic mechanisms underlying predator-induced defences at a qualitatively novel level.

Results

We report 230 differentially expressed genes (158 up- and 72 down-regulated) identified in at least two of three different assembly approaches. Several of the differentially regulated genes belong to families of paralogous genes. The most prominent classes amongst the up-regulated genes include cuticle genes, zinc-metalloproteinases and vitellogenin genes. Furthermore, several genes from this group code for proteins recruited in chromatin-reorganization or regulation of the cell cycle (cyclins). Down-regulated gene classes include C-type lectins, proteins involved in lipogenesis, and other families, some of which encode proteins with no known molecular function.

Conclusions

The RNA-Seq transcriptome data presented in this study provide important insights into gene regulatory patterns underlying predator-induced defences. In particular, we characterized different effector genes and gene families found to be regulated in Daphnia in response to the presence of an invertebrate predator. These effector genes are mostly in agreement with expectations based on observed phenotypic changes including morphological alterations, i.e., expression of proteins involved in formation of protective structures and in cuticle strengthening, as well as proteins required for resource re-allocation. Our findings identify key genetic pathways associated with anti-predator defences.

Electronic supplementary material

The online version of this article (doi:10.1186/s12983-015-0109-x) contains supplementary material, which is available to authorized users.

The effect of spontaneous mutations on competitive ability

Understanding the impact of spontaneous mutations on fitness has many theoretical and practical applications in biology. Although mutational effects on individual morphological or life-history characters have been measured in several classic genetic model systems, there are few estimates of the rate of decline due to mutation for complex fitness traits. Here, we estimate the effects of mutation on competitive ability, an important complex fitness trait, in a model system for ecological and evolutionary genomics, Daphnia. Competition assays were performed to compare fitness between mutation-accumulation (MA) lines and control lines from eight different genotypes from two populations of Daphnia pulicaria after 30 and 65 generations of mutation accumulation. Our results show a fitness decline among MA lines relative to controls as expected, but highlight the influence of genomic background on this effect. In addition, in some assays, MA lines outperform controls providing insight into the frequency of beneficial mutations.

Identifying multiple origins of polyploid taxa: a multilocus study of the hybrid cloak fern (Astrolepis integerrima; Pteridaceae)

PREMISE OF THE STUDY

Molecular studies have shown that multiple origins of polyploid taxa are the rule rather than the exception. To understand the distribution and ecology of polyploid species and the evolutionary significance of polyploidy in general, it is important to delineate these independently derived lineages as accurately as possible. Although gene flow among polyploid lineages and backcrossing to their diploid parents often confound this process, such post origin gene flow is very infrequent in asexual polyploids. In this study, we estimate the number of independent origins of the apomictic allopolyploid fern Astrolepis integerrima, a morphologically heterogeneous species most common in the southwestern United States and Mexico, with outlying populations in the southeastern United States and the Caribbean.

METHODS

Plastid DNA sequence and AFLP data were obtained from 33 A. integerrima individuals. Phylogenetic analysis of the sequence data and multidimensional clustering of the AFLP data were used to identify independently derived lineages.

KEY RESULTS

Analysis of the two datasets identified 10 genetic groups within the 33 analyzed samples. These groups suggest a minimum of 10 origins of A. integerrima in the northern portion of its range, with both putative parents functioning as maternal donors, both supplying unreduced gametes, and both contributing a significant portion of their genetic diversity to the hybrids.

CONCLUSIONS

Our results highlight the extreme cryptic genetic diversity and systematic complexity that can underlie a single polyploid taxon.

Transcontinental phylogeography of the Daphnia pulex species complex

Daphnia pulex is quickly becoming an attractive model species in the field of ecological genomics due to the recent release of its complete genome sequence, a wide variety of new genomic resources, and a rich history of ecological data. Sequences of the mitochondrial NADH dehydrogenase subunit 5 and cytochrome c oxidase subunit 1 genes were used to assess the global phylogeography of this species, and to further elucidate its phylogenetic relationship to other members of the Daphnia pulex species complex. Using both newly acquired and previously published data, we analyzed 398 individuals from collections spanning five continents. Eleven strongly supported lineages were found within the D. pulex complex, and one lineage in particular, panarctic D. pulex, has very little phylogeographical structure and a near worldwide distribution. Mismatch distribution, haplotype network, and population genetic analyses are compatible with a North American origin for this lineage and subsequent spatial expansion in the Late Pleistocene. In addition, our analyses suggest that dispersal between North and South America of this and other species in the D. pulex complex has occurred multiple times, and is predominantly from north to south. Our results provide additional support for the evolutionary relationships of the eleven main mitochondrial lineages of the D. pulex complex. We found that the well-studied panarctic D. pulex is present on every continent except Australia and Antarctica. Despite being geographically very widespread, there is a lack of strong regionalism in the mitochondrial genomes of panarctic D. pulex--a pattern that differs from that of most studied cladocerans. Moreover, our analyses suggest recent expansion of the panarctic D. pulex lineage, with some continents sharing haplotypes. The hypothesis that hybrid asexuality has contributed to the recent and unusual geographic success of the panarctic D. pulex lineage warrants further study.

Histaminergic signaling in the central nervous system of Daphnia and a role for it in the control of phototactic behavior

Daphnia magna and Daphnia pulex are well-established model organisms in the fields of ecotoxicology and toxicogenomics. Among the many assays used for determining the effects of environmental and anthropogenic stressors on these animals is monitoring for changes in their phototactic behavior. In most arthropods, histamine has been shown to play a key role in the visual system. Currently, nothing is known about histaminergic signaling in either D. magna or D. pulex. Here, a combination of immunohistochemistry and genome mining was used to identify and characterize the histaminergic systems in these daphnids. In addition, a behavioral assay was used to assess the role of histamine in their phototactic response to ultraviolet (UV) light exposure. An extensive network of histaminergic somata, axons and neuropil was identified via immunohistochemistry within the central nervous system of both daphnids, including labeling of putative photoreceptors in the compound eye and projections from these cells to the brain. Mining of the D. pulex genome using known Drosophila melanogaster proteins identified a putative ortholog of histidine decarboxylase (the rate-limiting biosynthetic enzyme for histamine), as well as two putative histamine-gated chloride channels (hclA and hclB orthologs). Exposure of D. magna to cimetidine, an H2 receptor antagonist known to block both hclA and hclB in D. melanogaster, inhibited their negative phototactic response to UV exposure in a reversible, time-dependent manner. Taken collectively, our results show that an extensive histaminergic system is present in Daphnia species, including the visual system, and that this amine is involved in the control of phototaxis in these animals.

Genomic analyses of the Daphnia pulex peptidome

Genome mining has provided a valuable tool for peptide discovery in many species, yet no crustacean has undergone this analysis. Currently, the only crustacean with a sequenced genome is the cladoceran Daphnia pulex, a model organism in many fields of biology. Here, we have mined the D. pulex genome for peptide-encoding genes. For each gene identified, the encoded precursor protein was deduced, and its mature peptides predicted. Twenty-four peptide-encoding genes were identified, including ones predicted to produce members of the A-type allatostatin, B-type allatostatin, C-type allatostatin, allatotropin (ATR), bursicon α, bursicon β, calcitonin-like diuretic hormone, corazonin, crustacean cardioactive peptide, crustacean hyperglycemic hormone, ecdysis-triggering hormone, eclosion hormone (EH), insulin-like peptide (ILP), molt-inhibiting hormone, neuropeptide F, orcokinin (two genes), pigment-dispersing hormone, proctolin, red pigment concentrating hormone/adipokinetic hormone (RPCH/AKH), short neuropeptide F, SIFamide, sulfakinin, and tachykinin-related peptide (TRP) families/subfamilies. In total, 96 peptides were predicted from these genes. Our identification of isoforms of corazonin, EH, ILP, proctolin, RPCH/AKH, sulfakinin and TRP are the first for D. pulex, while our prediction of ATR from this species is the first from any crustacean. The number of peptides predicted in our study shows the power of genome mining for peptide discovery, and provides a model for future genomic analyses of the peptidomes of other crustaceans. In addition, the data presented in our study provide foundations for future molecular, biochemical, anatomical, and physiological investigation of peptidergic signaling in D. pulex and other cladoceran species.

Bioinformatic analyses of the publicly accessible crustacean expressed sequence tags (ESTs) reveal numerous novel neuropeptide-encoding precursor proteins, including ones from members of several little studied taxa

ESTs have been generated for many crustacean species, providing an invaluable resource for peptide discovery in members of this arthropod subphylum. Here, these data were mined for novel peptide-encoding transcripts, with the mature peptides encoded by them predicted using a combination of online peptide prediction programs and homology to known arthropod sequences. In total, 70 mature full-length/partial peptides representing members of 16 families/subfamilies were predicted, the vast majority being novel; the species from which the peptides were identified included members of the Branchiopoda (Daphnia carinata and Triops cancriformis), Maxillopoda (Caligus clemensi, Caligus rogercresseyi, Lepeophtheirus salmonis and Lernaeocera branchialis) and Malacostraca (Euphausia superba, Marsupenaeus japonicus, Penaeus monodon, Homarus americanus, Petrolisthes cinctipes, Callinectes sapidus and Portunus trituberculatus). Of particular note were the identifications of an intermediate between the insect adipokinetic hormones and crustacean red pigment concentrating hormone and a modified crustacean cardioactive peptide from the daphnid D. carinata; Arg(7)-corazonin was also deduced from this species, the first identification of a corazonin from a non-decapod crustacean. Our data also include the first reports of members of the calcitonin-like diuretic hormone, FMRFamide-related peptide (neuropeptide F subfamily) and orcokinin families from members of the Copepoda. Moreover, the prediction of a bursicon alpha from the euphausid E. superba represents the first peptide identified from any member of the basal eucaridean order Euphausiacea. In addition, large collections of insect eclosion hormone- and neuroparsin-like peptides were identified from a variety of species, greatly expanding the number of known members of these families in crustaceans.

The porcelain crab transcriptome and PCAD, the porcelain crab microarray and sequence database

Background

With the emergence of a completed genome sequence of the freshwater crustacean Daphnia pulex, construction of genomic-scale sequence databases for additional crustacean sequences are important for comparative genomics and annotation. Porcelain crabs, genus Petrolisthes, have been powerful crustacean models for environmental and evolutionary physiology with respect to thermal adaptation and understanding responses of marine organisms to climate change. Here, we present a large-scale EST sequencing and cDNA microarray database project for the porcelain crab Petrolisthes cinctipes.

Methodology/Principal Findings

A set of ∼30K unique sequences (UniSeqs) representing ∼19K clusters were generated from ∼98K high quality ESTs from a set of tissue specific non-normalized and mixed-tissue normalized cDNA libraries from the porcelain crab Petrolisthes cinctipes. Homology for each UniSeq was assessed using BLAST, InterProScan, GO and KEGG database searches. Approximately 66% of the UniSeqs had homology in at least one of the databases. All EST and UniSeq sequences along with annotation results and coordinated cDNA microarray datasets have been made publicly accessible at the Porcelain Crab Array Database (PCAD), a feature-enriched version of the Stanford and Longhorn Array Databases.

Conclusions/Significance

The EST project presented here represents the third largest sequencing effort for any crustacean, and the largest effort for any crab species. Our assembly and clustering results suggest that our porcelain crab EST data set is equally diverse to the much larger EST set generated in the Daphnia pulex genome sequencing project, and thus will be an important resource to the Daphnia research community. Our homology results support the pancrustacea hypothesis and suggest that Malacostraca may be ancestral to Branchiopoda and Hexapoda. Our results also suggest that our cDNA microarrays cover as much of the transcriptome as can reasonably be captured in EST library sequencing approaches, and thus represent a rich resource for studies of environmental genomics.

Identification of putative peptide paracrines/hormones in the water flea Daphnia pulex (Crustacea; Branchiopoda; Cladocera) using transcriptomics and immunohistochemistry

The cladoceran crustacean Daphnia pulex has emerged as a model species for many biological fields, in particular environmental toxicology and toxicogenomics. Recently, this species has been the subject of an extensive transcriptome project, resulting in the generation and public deposition of over 150,000 expressed sequence tags (ESTs). This resource makes D. pulex an excellent model for protein discovery using bioinformatics. Here, in silico searches of the D. pulex EST database were conducted to identify transcripts encoding putative peptide precursors. Moreover, the mature peptides contained within the deduced prepro-hormones were predicted using online peptide processing programs and homology to known arthropod isoforms. In total, 63 putative peptide-encoding ESTs were identified encompassing 14 distinct peptide families/subfamilies: A-type allatostatin, B-type allatostatin, C-type allatostatin, bursicon (both alpha and beta subunit peptides), crustacean cardioactive peptide (CCAP), crustacean hyperglycemic hormone (CHH)/ion transport peptide (both CHH- and moult-inhibiting hormone-like subfamilies), diuretic hormone (calcitonin-like), ecdysis-triggering hormone (ETH), FMRFamide (both neuropeptide F and short neuropeptide F subfamilies), orcokinin and pigment dispersing hormone. From these transcripts, the structures of 76 full-length/partial peptides were predicted, which included the first C-type allatostatin-like peptide identified from a crustacean, the first crustacean calcitonin-like diuretic hormone, an undescribed CCAP isoform, two hitherto unknown ETH variants, and two new orcokinins. Neuronal localization of several of the identified peptide families was confirmed using immunohistochemitry (i.e. A-type allatostatin, CCAP, FMRFamide and PDH). In addition, immunohistochemical analyses identified other putative neuropeptides for which no ESTs had been found (i.e. corazonin, insect kinin, proctolin, red pigment concentrating hormone, SIFamide, sulfakinin and tachykinin-related peptide). Collectively, the data presented here not only catalog an extensive array of putative D. pulex peptide paracrines/hormones, but also provide a strong foundation for future investigations of the effects of environmental/anthropogenic stressors on peptidergic control in this model organism.