Conserved Transcription Factors Steer Growth-Related Genomic Programs in Daphnia

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

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

Environmental RNA: A Revolution in Ecological Resolution?

Current advancements in environmental RNA (eRNA) exploit its relatively fast turnover rate relative to environmental DNA (eDNA) to assess 'metabolically active' or temporally/spatially recent community diversity. However, this focus significantly underutilizes the trove of potential ecological information encrypted in eRNA. Here, we argue for pushing beyond current species-level eDNA detection capabilities by using eRNA to detect any organisms with unique eRNA profiles, potentially including different life-history stages, sexes, or even specific phenotypes within a species. We also discuss the future of eRNA as a means of assessing the physiological status of organisms and the ecological health of populations and communities, reflecting ecosystem-level conditions. We posit that eRNA has the potential to significantly improve the resolution of organism detection, biological monitoring, and biomonitoring applications in ecology.

A complex gene regulatory architecture underlies the development and evolution of cuticle morphology in Drosophila

The cuticle of insects is decorated with non-sensory hairs called trichomes. A few Drosophila species independently lost most of the dorso-lateral trichomes on first instar larvae. Genetic experiments revealed that this naked cuticle phenotype was caused by the evolution of enhancer function at the ovo/shavenbaby (ovo/svb) locus. Here we explore how this discovery catalyzed major new insights into morphological evolution in different developmental contexts, enhancer pleiotropy in gene regulation and the functionality and evolution of the Svb gene regulatory network (GRN). Taken together this highlights the importance of understanding the architecture and evolution of gene regulatory networks in detail and the great potential for further study of the Svb GRN.

Reduced stress defence responses contribute to the higher toxicity of a pesticide under warming

There is a pressing need to identify the molecular mechanisms underlying the, often magnifying, interactive effects between contaminants and natural stressors. Here we test our hypothesis that lower general stress defence responses contribute to synergistic interactions between stressors. We focus on the widespread pattern that many contaminants are more toxic at higher temperatures. Specifically, we tested the effects of an environmentally realistic low-effect and high-effect concentration of the pesticide chlorpyrifos under warming at the gene expression level in the northern house mosquito Culex pipiens molestus (Forskal, 1775). By applying the independent action model for combined stressors on RNA-sequencing data, we identified interactive gene expression patterns under combined exposure to chlorpyrifos and warming for general stress defence responses: protection of macromolecules, antioxidant processes, detoxification and energy metabolism/allocation. Most of these general stress defence response genes showed upregulated antagonistic interactions (i.e., were less upregulated than expected under the independent action model). This indicates that when pesticide exposure was combined with warming, the general stress defence responses were no longer buffering increased stress levels, which may contribute to a higher sensitivity to toxicants under warming. These upregulated antagonistic interactions were stronger for the high-effect chlorpyrifos concentration, indicating that exposure to this concentration under warming was most stressful. Our results highlight that quantitative analysis of the frequency and strength of the interaction types of general stress defence response genes, specifically focusing on antagonistic upregulations and synergistic downregulations, may advance our understanding of how natural stressors modify the toxicity of contaminants.

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

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

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.

Core transcriptional signatures of phase change in the migratory locust

Phenotypic plasticity plays fundamental roles in successful adaptation of animals in response to environmental variations. Here, to reveal the transcriptome reprogramming in locust phase change, a typical phenotypic plasticity, we conducted a comprehensive analysis of multiple phase-related transcriptomic datasets of the migratory locust. We defined PhaseCore genes according to their contribution to phase differentiation by the adjustment for confounding principal components analysis algorithm (AC-PCA). Compared with other genes, PhaseCore genes predicted phase status with over 87.5% accuracy and displayed more unique gene attributes including the faster evolution rate, higher CpG content and higher specific expression level. Then, we identified 20 transcription factors (TFs) named PhaseCoreTF genes that are associated with the regulation of PhaseCore genes. Finally, we experimentally verified the regulatory roles of three representative TFs (Hr4, Hr46, and grh) in phase change by RNAi. Our findings revealed that core transcriptional signatures are involved in the global regulation of locust phase changes, suggesting a potential common mechanism underlying phenotypic plasticity in insects. The expression and network data are accessible in an online resource called LocustMine (http://www.locustmine.org:8080/locustmine).

The mlpt/Ubr3/Svb module comprises an ancient developmental switch for embryonic patterning

Small open reading frames (smORFs) encoding 'micropeptides' exhibit remarkable evolutionary complexity. Conserved peptides encoded by mille-pattes (mlpt)/polished rice (pri)/tarsal less (tal) are essential for embryo segmentation in Tribolium but, in Drosophila, function in terminal epidermal differentiation and patterning of adult legs. Here, we show that a molecular complex identified in Drosophila epidermal differentiation, comprising Mlpt peptides, ubiquitin-ligase Ubr3 and transcription factor Shavenbaby (Svb), represents an ancient developmental module required for early insect embryo patterning. We find that loss of segmentation function for this module in flies evolved concomitantly with restriction of Svb expression in early Drosophila embryos. Consistent with this observation, artificially restoring early Svb expression in flies causes segmentation defects that depend on mlpt function, demonstrating enduring potency of an ancestral developmental switch despite evolving embryonic patterning modes. These results highlight the evolutionary plasticity of conserved molecular complexes under the constraints of essential genetic networks.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Towards a Dynamic Interaction Network of Life to unify and expand the evolutionary theory

The classic Darwinian theory and the Synthetic evolutionary theory and their linear models, while invaluable to study the origins and evolution of species, are not primarily designed to model the evolution of organisations, typically that of ecosystems, nor that of processes. How could evolutionary theory better explain the evolution of biological complexity and diversity? Inclusive network-based analyses of dynamic systems could retrace interactions between (related or unrelated) components. This theoretical shift from a Tree of Life to a Dynamic Interaction Network of Life, which is supported by diverse molecular, cellular, microbiological, organismal, ecological and evolutionary studies, would further unify evolutionary biology.

Identification of General Patterns of Sex-Biased Expression in Daphnia, a Genus with Environmental Sex Determination

Daphnia reproduce by cyclic-parthenogenesis, where phases of asexual reproduction are intermitted by sexual production of diapause stages. This life cycle, together with environmental sex determination, allow the comparison of gene expression between genetically identical males and females. We investigated gene expression differences between males and females in four genotypes of Daphnia magna and compared the results with published data on sex-biased gene expression in two other Daphnia species, each representing one of the major phylogenetic clades within the genus. We found that 42% of all annotated genes showed sex-biased expression in D. magna This proportion is similar both to estimates from other Daphnia species as well as from species with genetic sex determination, suggesting that sex-biased expression is not reduced under environmental sex determination. Among 7453 single copy, one-to-one orthologs in the three Daphnia species, 707 consistently showed sex-biased expression and 675 were biased in the same direction in all three species. Hence these genes represent a core-set of genes with consistent sex-differential expression in the genus. A functional analysis identified that several of them are involved in known sex determination pathways. Moreover, 75% were overexpressed in females rather than males, a pattern that appears to be a general feature of sex-biased gene expression in Daphnia.