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

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

Predator-induced shape plasticity in Daphnia pulex.

All animals and plants respond to changes in the environment during their life cycle. This flexibility is known as phenotypic plasticity and allows organisms to cope with variable environments. A common source of environmental variation is predation risk, which describes the likelihood of being attacked and killed by a predator. Some species can respond to the level of predation risk by producing morphological defences against predation. A classic example is the production of so-called 'neckteeth' in the water flea, Daphnia pulex, which defend against predation from Chaoborus midge larvae. Previous studies of this defence have focussed on changes in pedestal size and the number of spikes along a gradient of predation risk. Although these studies have provided a model for continuous phenotypic plasticity, they do not capture the whole-organism shape response to predation risk. In contrast, studies in fish and amphibians focus on shape as a complex, multi-faceted trait made up of different variables. In this study, we analyse how multiple aspects of shape change in D. pulex along a gradient of predation risk from Chaoborus flavicans. These changes are dominated by the neckteeth defence, but there are also changes in the size and shape of the head and the body. We detected change in specific modules of the body plan and a level of integration among modules. These results are indicative of a complex, multi-faceted response to predation and provide insight into how predation risk drives variation in shape and size at the level of the whole organism.

Adaptive phenotypic plasticity is under stabilizing selection in Daphnia.

The adaptive nature of phenotypic plasticity is widely documented. However, little is known about the evolutionary forces that shape genetic variation of plasticity within populations. Whether genetic variation in plasticity is driven by stabilizing or diversifying selection and whether the strength of such forces remains constant through time, remain open questions. Here, we address this issue by assessing the evolutionary forces that shape genetic variation in antipredator developmental plasticity of Daphnia pulex. Antipredator plasticity in D. pulex is characterized by the growth of a pedestal and spikes in the dorsal head region upon exposure to predator cue. We characterized genetic variation in plasticity using a method that describes the entire dorsal shape amongst >100 D. pulex strains recently derived from the wild. We observed the strongest reduction in genetic variation in dorsal areas where plastic responses were greatest, consistent with stabilizing selection. We compared mutational variation (Vm) to standing variation (Vg) and found that Vg/Vm is lowest in areas of greatest plasticity, again consistent with stabilizing selection. Our results suggest that stabilizing selection operates directly on phenotypic plasticity in Daphnia and provide a rare glimpse into the evolution of fitness-related traits in natural populations.

Genetic Variation in Reproductive Investment Across an Ephemerality Gradient in Daphnia pulex.

Species across the tree of life can switch between asexual and sexual reproduction. In facultatively sexual species, the ability to switch between reproductive modes is often environmentally dependent and subject to local adaptation. However, the ecological and evolutionary factors that influence the maintenance and turnover of polymorphism associated with facultative sex remain unclear. We studied the ecological and evolutionary dynamics of reproductive investment in the facultatively sexual model species, Daphnia pulex. We found that patterns of clonal diversity, but not genetic diversity varied among ponds consistent with the predicted relationship between ephemerality and clonal structure. Reconstruction of a multi-year pedigree demonstrated the coexistence of clones that differ in their investment into male production. Mapping of quantitative variation in male production using lab-generated and field-collected individuals identified multiple putative quantitative trait loci (QTL) underlying this trait, and we identified a plausible candidate gene. The evolutionary history of these QTL suggests that they are relatively young, and male limitation in this system is a rapidly evolving trait. Our work highlights the dynamic nature of the genetic structure and composition of facultative sex across space and time and suggests that quantitative genetic variation in reproductive strategy can undergo rapid evolutionary turnover.

Dissecting the Transcriptomic Basis of Phenotypic Evolution in an Aquatic Keystone Grazer.

Knowledge of the molecular basis of phenotypic responses to environmental cues is key to understanding the process of adaptation. Insights to adaptation at an evolutionary time scale can be gained by observing organismal responses before and after a shift in environmental conditions, but such observations can rarely be made. Using the ecological and genomic model Daphnia, we linked transcriptomic responses and phosphorus (P)-related phenotypic traits under high and low P availability. We mapped weighted gene coexpression networks to traits previously assessed in resurrected ancient (600 years old) and modern Daphnia pulicaria from a lake with a historic shift in P-enrichment. Subsequently, we assessed evolutionary conservation or divergence in transcriptional networks of the same isolates. We discovered highly preserved gene networks shared between ancient genotypes and their modern descendants, but also detected clear evidence of transcriptional divergence between these evolutionarily separated genotypes. Our study highlights that phenotypic evolution is a result of molecular fine-tuning on different layers ranging from basic cellular responses to higher order phenotypes. In a broader context, these findings advance our understanding how populations are able to persist throughout major environmental shifts.

The transcriptomic and proteomic responses of Daphnia pulex to changes in temperature and food supply comprise environment-specific and clone-specific elements.

BACKGROUND: Regulatory adjustments to acute and chronic temperature changes are highly important for aquatic ectotherms because temperature affects their metabolic rate as well as the already low oxygen concentration in water, which can upset their energy balance. This also applies to severe changes in food supply. Thus, we studied on a molecular level (transcriptomics and/or proteomics) the immediate responses to heat stress and starvation and the acclimation to different temperatures in two clonal isolates of the model microcrustacean Daphnia pulex from more or less stressful environments, which showed a higher (clone M) or lower (clone G) tolerance to heat and starvation. RESULTS: The transcriptomic responses of clone G to acute heat stress (from 20 °C to 30 °C) and temperature acclimation (10 °C, 20 °C, and 24 °C) and the proteomic responses of both clones to acute heat, starvation, and heat-and-starvation stress comprised environment-specific and clone-specific elements. Acute stress (in particular heat stress) led to an early upregulation of stress genes and proteins (e.g., molecular chaperones) and a downregulation of metabolic genes and proteins (e.g., hydrolases). The transcriptomic responses to temperature acclimation differed clearly. They also varied depending on the temperature level. Acclimation to higher temperatures comprised an upregulation of metabolic genes and, in case of 24 °C acclimation, a downregulation of genes for translational processes and collagens. The proteomic responses of the clones M and G differed at any type of stress. Clone M showed markedly stronger and less stress-specific proteomic responses than clone G, which included the consistent expression of a specific heat shock protein (HSP60) and vitellogenin (VTG-SOD). CONCLUSIONS: The expression changes under acute stress can be interpreted as a switch from standard products of gene expression to stress-specific products. The expression changes under temperature acclimation probably served for an increase in energy intake (via digestion) and, if necessary, a decrease in energy expenditures (e.g, for translational processes). The stronger and less stress-specific proteomic responses of clone M indicate a lower degree of cell damage and an active preservation of the energy balance, which allowed adequate proteomic responses under stress, including the initiation of resting egg production (VTG-SOD expression) as an emergency reaction.

Conserved Transcription Factors Steer Growth-Related Genomic Programs in Daphnia.

Ecological genomics aims to understand the functional association between environmental gradients and the genes underlying adaptive traits. Many genes that are identified by genome-wide screening in ecologically relevant species lack functional annotations. Although gene functions can be inferred from sequence homology, such approaches have limited power. Here, we introduce ecological regulatory genomics by presenting an ontology-free gene prioritization method. Specifically, our method combines transcriptome profiling with high-throughput cis-regulatory sequence analysis in the water fleas Daphnia pulex and Daphnia magna. It screens coexpressed genes for overrepresented DNA motifs that serve as transcription factor binding sites, thereby providing insight into conserved transcription factors and gene regulatory networks shaping the expression profile. We first validated our method, called Daphnia-cisTarget, on a D. pulex heat shock data set, which revealed a network driven by the heat shock factor. Next, we performed RNA-Seq in D. magna exposed to the cyanobacterium Microcystis aeruginosa. Daphnia-cisTarget identified coregulated gene networks that associate with the moulting cycle and potentially regulate life history changes in growth rate and age at maturity. These networks are predicted to be regulated by evolutionary conserved transcription factors such as the homologues of Drosophila Shavenbaby and Grainyhead, nuclear receptors, and a GATA family member. In conclusion, our approach allows prioritising candidate genes in Daphnia without bias towards prior knowledge about functional gene annotation and represents an important step towards exploring the molecular mechanisms of ecological responses in organisms with poorly annotated genomes.

The Role of Omics in the Application of Adverse Outcome Pathways for Chemical Risk Assessment.

In conjunction with the second International Environmental Omics Symposium (iEOS) conference, held at the University of Liverpool (United Kingdom) in September 2014, a workshop was held to bring together experts in toxicology and regulatory science from academia, government and industry. The purpose of the workshop was to review the specific roles that high-content omics datasets (eg, transcriptomics, metabolomics, lipidomics, and proteomics) can hold within the adverse outcome pathway (AOP) framework for supporting ecological and human health risk assessments. In light of the growing number of examples of the application of omics data in the context of ecological risk assessment, we considered how omics datasets might continue to support the AOP framework. In particular, the role of omics in identifying potential AOP molecular initiating events and providing supportive evidence of key events at different levels of biological organization and across taxonomic groups was discussed. Areas with potential for short and medium-term breakthroughs were also discussed, such as providing mechanistic evidence to support chemical read-across, providing weight of evidence information for mode of action assignment, understanding biological networks, and developing robust extrapolations of species-sensitivity. Key challenges that need to be addressed were considered, including the need for a cohesive approach towards experimental design, the lack of a mutually agreed framework to quantitatively link genes and pathways to key events, and the need for better interpretation of chemically induced changes at the molecular level. This article was developed to provide an overview of ecological risk assessment process and a perspective on how high content molecular-level datasets can support the future of assessment procedures through the AOP framework.

Adjustments of serine proteases of Daphnia pulex in response to temperature changes.

Elevated temperatures considerably challenge aquatic invertebrates, and enhanced energy metabolism and protein turnover require adjustments of digestion. In Daphnia, the serine proteases chymotrypsin and trypsin represent the major proteolytic enzymes. Daphnia pulex acclimated to different temperature conditions or subjected to acute heat stress showed increased expression level of serine proteases with rising temperatures. Transcripts of trypsin isoforms were always present in higher amounts than observed for chymotrypsin. Additionally, trypsin isoform transcripts were induced by elevated temperatures to a larger extent. Correspondingly, trypsin activity dominated in cold-acclimated animals. However, the enzymatic activity of chymotrypsin increased at elevated temperatures, whereas trypsin activity slightly decreased, resulting in a shift to dominating chymotrypsin activity in warm-acclimated animals. Zymograms revealed eight bands with proteolytic activity in the range of 20 to 86 kDa. The single bands were assigned to trypsin or chymotrypsin activity applying specific inhibitors or from casein cleavage products identified by mass spectrometric analysis. The total amount of proteolytic activity was elevated with acclimation temperature increase and showed a transient decrease under acute heat stress. The contribution of the different isoforms to protein digestion indicated induction of chymotrypsin with increasing acclimation temperature. For trypsin, the share of one isoform decreased with elevated temperature, while another isoform was enhanced. Thus differential expression of serine proteases was observed in response to chronic and acute temperature changes. The observed phenotypic plasticity adjusts the set of active proteases to the altered needs of protein metabolism optimizing protein digestion for the temperature conditions experienced in the habitat.

Differential transcriptomic responses of ancient and modern Daphnia genotypes to phosphorus supply.

Little is known about the role of transcriptomic changes in driving phenotypic evolution in natural populations, particularly in response to anthropogenic environmental change. Previous analyses of Daphnia genotypes separated by centuries of evolution in a lake using methods in resurrection ecology revealed striking genetic and phenotypic shifts that were highly correlated with anthropogenic environmental change, specifically phosphorus (P)-driven nutrient enrichment (i.e. eutrophication). Here, we compared the transcriptomes of two ancient (~700-year-old) and two modern (~10-year-old) genotypes in historic (low P) and contemporary (high P) environmental conditions using microarrays. We found considerable transcriptomic variation between 'ancient' and 'modern' genotypes in both treatments, with stressful (low P) conditions eliciting differential expression (DE) of a larger number of genes. Further, more genes were DE between 'ancient' and 'modern' genotypes than within these groups. Expression patterns of individual genes differed greatly among genotypes, suggesting that different transcriptomic responses can result in similar phenotypes. While this confounded patterns between 'ancient' and 'modern' genotypes at the gene level, patterns were discernible at the functional level: annotation of DE genes revealed particular enrichment of genes involved in metabolic pathways in response to P-treatments. Analyses of gene families suggested significant DE in pathways already known to be important in dealing with P-limitation in Daphnia as well as in other organisms. Such observations on genotypes of a single natural population, separated by hundreds of years of evolution in contrasting environmental conditions before and during anthropogenic environmental changes, highlight the important role of transcriptional mechanisms in the evolutionary responses of populations.

Hypoxia-inducible haemoglobins of Daphnia pulex and their role in the response to acute and chronic temperature increase.

Daphnia pulex is challenged by severe oxygen and temperature changes in its habitat. In response to hypoxia, the equipment of oxygen transport proteins is adjusted in quantity and quality by differential expression of haemoglobin isoforms. This study focuses on the response of 20°C acclimated animals to elevated temperature using transcriptomic and proteomic approaches. Acute temperature stress (30°C) induced the hypoxia-inducible Hb isoforms most strongly, resulting in an increase of the haemoglobin mRNA pool by 70% within 8h. Long-term-acclimation to moderately elevated temperature (24°C) only evoked minor changes of the Hb mRNA suite. Nevertheless, the concentration of the hemolymph pool of haemoglobin was elevated by 80%. In this case, the constitutive Hb isoforms showed the strongest increase, with Hb01 and Hb02 contributing by 64% to the total amount of respiratory protein. The regulation patterns upon acute temperature stress likely reflect temperature-induced tissue hypoxia, whereas in case of persisting exposure to moderately elevated temperature, acclimation processes enabled the successful return to oxygen homeostasis. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Acute changes in temperature or oxygen availability induce ROS fluctuations in Daphnia magna linked with fluctuations of reduced and oxidized glutathione, catalase activity and gene (haemoglobin) expression.

BACKGROUND INFORMATION: ROS (reactive oxygen species) as well as components of the antioxidant redox systems may act as signals. To link acute environmental change with gene expression, changes in ROS and GSH/GSSG (reduced/oxidized glutathione) level were measured upon acute changes in temperature or oxygen availability in the aquatic key species Daphnia magna together with HIF-1 (hypoxia-inducible factor 1)-mediated Hb (haemoglobin) expression. RESULTS: Acute exposures to 30°C or hypoxia, which induced tissue hypoxia (and possibly elevated mitochondrial ROS production), caused resembling fluctuations of ROS and GSH levels, with frequency and number of peaks increasing and their delay decreasing with the magnitude of environmental change (size of tissue hypoxia). Acute hyperoxia induced an initial decrease in ROS level. Evidence is also provided for the promoting effects of ROS on catalase activity. A signalling function of the ROS fluctuations upon acute changes in temperature was found in the case of Hb, the expression of which is known to respond to temperature changes, by detecting corresponding time courses of both transcription and protein formation. CONCLUSION: ROS-dependent signalling was affected by changes in temperature or oxygen availability. Feedback interactions between ROS and the glutathione redox system, possibly driven by elevated mitochondrial ROS production, likely contributed to the appearance of the ROS and GSH fluctuations upon acute environmental change. Fluctuating ROS levels, which reflect for the magnitude of environmental change, could be a way to transfer information on ROS production to subsequent processes (gene expression) while avoiding too-high and damaging ROS levels.