Mapping the expression of the sex determining factor Doublesex1 in Daphnia magna using a knock-in reporter

Zooplankton-based adverse outcome pathways: A tool for assessing endocrine disrupting compounds in aquatic environments

Endocrine disrupting compounds (EDCs) pose a significant ecological risk, particularly in aquatic ecosystems. EDCs have become a focal point in ecotoxicology, and their identification and regulation have become a priority. Zooplankton have gained global recognition as bioindicators, benefiting from rigorous standardization and regulatory validation processes. This review aims to provide a comprehensive summary of zooplankton-based adverse outcome pathways (AOPs) with a focus on EDCs as toxicants and the utilisation of freshwater zooplankton as bioindicators in ecotoxicological assessments. This review presents case studies in which zooplankton have been used in the development of AOPs, emphasizing the identification of molecular initiating events (MIEs) and key events (KEs) specific to zooplankton exposed to EDCs. Zooplankton-based AOPs may become an important resource for understanding the intricate processes by which EDCs impair the endocrine system. Furthermore, the data sources, experimental approaches, advantages, and challenges associated with zooplankton-based AOPs are discussed. Zooplankton-based AOPs framework can provide vital tools for consolidating toxicological knowledge into a structured toxicity pathway of EDCs, offering a transformative platform for facilitating enhanced risk assessment and chemical regulation.

Identification of gene isoforms and their switching events between male and female embryos of the parthenogenetic crustacean Daphnia magna

The cladoceran crustacean Daphnia exhibits phenotypic plasticity, a phenomenon that leads to diverse phenotypes from one genome. Alternative usage of gene isoforms has been considered a key gene regulation mechanism for controlling different phenotypes. However, to understand the phenotypic plasticity of Daphnia, gene isoforms have not been comprehensively analyzed. Here we identified 25,654 transcripts derived from the 9710 genes expressed during environmental sex determination of Daphnia magna using the long-read RNA-Seq with PacBio Iso-Seq. We found that 14,924 transcripts were previously unidentified and 5713 genes produced two or more isoforms. By a combination of Illumina short-read RNA-Seq, we detected 824 genes that implemented switching of the highest expressed isoform between females and males. Among the 824 genes, we found isoform switching of an ortholog of CREB-regulated transcription coactivator, a major regulator of carbohydrate metabolism in animals, and a correlation of this switching event with the sexually dimorphic expression of carbohydrate metabolic genes. These results suggest that a comprehensive catalog of isoforms may lead to understanding the molecular basis for environmental sex determination of Daphnia. We also infer the applicability of the full-length isoform analyses to the elucidation of phenotypic plasticity in Daphnia.

Novel peptide identified from viable-cell based phage display technique regulates growth cycle of Daphnia magna

Phage display is a widely used technique for selecting specific binding peptides, but presenting antigens in their natural form can be challenging, as protein coating may induce structural changes. In this study, we employed a whole cell-based phage display technique without a coating step to select peptides that bind specifically to Daphnia magna eggs. Boiled eggs were used as a control to ensure that antigens were presented in their natural forms. We identified a peptide, DEP1 (LYALPLSHLKSHGGG), with the highest binding affinity to D. magna eggs. DEP1 did not affect zebrafish eggs, but it inhibited normal hatching and reproductive ability in D. magna eggs, and hindered growth in neonates before their first ecdysis. Morphological analysis revealed that DEP1 caused intestinal damage and tissue abnormalities. Our findings demonstrate that the whole cell-based phage display technique is successful in presenting antigens in their natural form, and that the DEP1 peptide can be applied to regulate the growth cycle of D. magna. These results have implications for the use of phage display in environmental research and the potential use of DEP1 for hazardous organisms in aquatic systems.

Risk assessment of bisphenol analogues towards mortality, heart rate and stress-mediated gene expression in cladocerans Moina micrura

Bisphenol A (BPA) is a well-known endocrine-disrupting compound that causes several toxic effects on human and aquatic organisms. The restriction of BPA in several applications has increased the substituted toxic chemicals such as bisphenol F (BPF) and bisphenol S (BPS). A native tropical freshwater cladoceran, Moina micrura, was used as a bioindicator to assess the adverse effects of bisphenol analogues at molecular, organ, individual and population levels. Bisphenol analogues significantly upregulated the expressions of stress-related genes, which are the haemoglobin and glutathione S-transferase genes, but the sex determination genes such as doublesex and juvenile hormone analogue genes were not significantly different. The results show that bisphenol analogues affect the heart rate and mortality rate of M. micrura. The 48-h lethal concentration (LC₅₀) values based on acute toxicity for BPA, BPF and BPS were 611.6 µg L^-1, 632.0 µg L^-1 and 819.1 µg L^-1, respectively. The order of toxicity based on the LC₅₀ and predictive non-effect concentration values were as follows: BPA > BPF > BPS. Furthermore, the incorporated method combining the responses throughout the organisation levels can comprehensively interpret the toxic effects of bisphenol analogues, thus providing further understanding of the toxicity mechanisms. Moreover, the output of this study produces a comprehensive ecotoxicity assessment, which provides insights for the legislators regarding exposure management and mitigation of bisphenol analogues in riverine ecosystems.

CELF1 represses Doublesex1 expression via its 5' UTR in the crustacean Daphnia magna

In sex determination of the crustacean Daphnia magna, male-specific expression of DM-domain transcription factor Doublesex1 (Dsx1) orchestrates the male developmental program triggered by environmental stimuli. We previously identified the CELF1 ortholog as a candidate of proteins associated with the 5' UTR of the Dsx1α isoform. Here we report the CELF1-dependent suppression of Dsx1 expression in D. magna. During embryogenesis, CELF1 expression was not sexually dimorphic. Silencing of CELF1 led to the activation of Dsx1 expression both in female and male embryos. Overexpression of CELF1 in male embryos resulted in a reduction of Dsx1 expression. By these manipulations of CELF1 expression, the Dsx1 transcript level was not significantly changed. To investigate whether the CELF1 controls Dsx1 expression via its 5' UTR, we injected the GFP reporter mRNA having intact Dsx1α 5' UTR or mutated one lacking the GU-rich element (GRE) that is known as a binding site of the CELF1 ortholog. We found that deletion of the GRE significantly increased the reporter gene expression. These results indicate that CELF1 suppresses Dsx1 expression both in females and males, possibly at the post-transcriptional level. We speculate that CELF1 may avoid unintended Dsx1 expression and generation of sexual ambiguity by setting a threshold of Dsx1 expression.

Regulation of Doublesex1 Expression for Environmental Sex Determination in the Cladoceran Crustacean Daphnia

The cladoceran crustacean Daphnia produces only females by parthenogenesis in a healthy population. However, in response to environmental declines such as crowding and lack of foods, it produces eggs destined to become males that are genetically identical to females. During the development of the sexually committed eggs, DM domain-containing transcription factor Doublesex1 (Dsx1) orchestrates male trait formation globally both in somatic and gonadal tissues. Recent studies have revealed that Dsx1 expression is tightly controlled at transcriptional, post-transcriptional, and epigenetic levels to avoid sexual ambiguity. In this review, together with basic information on Dsx1 structure and expression, we introduce the multi-layered Dsx1 regulation and discuss how each regulation is interconnected for controlling male development in environmental sex-determining Daphnia.

Juvenile hormone synthesis and signaling disruption triggering male offspring induction and population decline in cladocerans (water flea): Review and adverse outcome pathway development

Juvenile hormone (JH) are a family of multifunctional hormones regulating larval development, molting, metamorphosis, reproduction, and phenotypic plasticity in arthropods. Based on its importance in arthropod life histories, many insect growth regulators (IGRs) mimicking JH have been designed to control harmful insects in agriculture and aquaculture. These JH analogs (JHAs) may also pose hazards to nontarget species by causing unexpected endocrine-disrupting (ED) effects such as molting and metamorphosis defects, larval lethality, and disruption of the sexual identity. This critical review summarizes the current knowledge of the JH-mediated effects in the freshwater cladoceran crustaceans such as Daphnia species on JHA-triggered endocrine disruptive outputs to establish a systematic understanding of JHA effects. Based on the current knowledge, adverse outcome pathways (AOPs) addressing the JHA-mediated ED effects in cladoceran leading to male offspring production and subsequent population decline were developed. The weight of evidence (WoE) of AOPs was assessed according to established guidelines. The review and AOP development aim to present the current scientific understanding of the JH pathway and provide a robust reference for the development of tiered testing strategies and new risk assessment approaches for JHAs in future ecotoxicological research and regulatory processes.

The freshwater water flea Daphnia magna NIES strain genome as a resource for CRISPR/Cas9 gene targeting: The glutathione S-transferase omega 2 gene

The water flea Daphnia magna is a small freshwater planktonic animal in the Cladocera. In this study, we assembled the genome of the D. magna NIES strain, which is widely used for gene targeting but has no reported genome. We used the long-read sequenced data of the Oxford nanopore sequencing tool for assembly. Using 3,231 genetic markers, the draft genome of the D. magna NIES strain was built into ten linkage groups (LGs) with 483 unanchored contigs, comprising a genome size of 173.47 Mb. The N50 value of the genome was 12.54 Mb and the benchmarking universal single-copy ortholog value was 98.8%. Repeat elements in the D. magna NIES genome were 40.8%, which was larger than other Daphnia spp. In the D. magna NIES genome, 15,684 genes were functionally annotated. To assess the genome of the D. magna NIES strain for CRISPR/Cas9 gene targeting, we selected glutathione S-transferase omega 2 (GST-O2), which is an important gene for the biotransformation of arsenic in aquatic organisms, and targeted it with an efficient make-up (25.0%) of mutant lines. In addition, we measured reactive oxygen species and antioxidant enzymatic activity between wild type and a mutant of the GST-O2 targeted D. magna NIES strain in response to arsenic. In this study, we present the genome of the D. magna NIES strain using GST-O2 as an example of gene targeting, which will contribute to the construction of deletion mutants by CRISPR/Cas9 technology.

Sense-overlapping lncRNA as a decoy of translational repressor protein for dimorphic gene expression

Long noncoding RNAs (lncRNAs) are vastly transcribed and extensively studied but lncRNAs overlapping with the sense orientation of mRNA have been poorly studied. We analyzed the lncRNA DAPALR overlapping with the 5´ UTR of the Doublesex1 (Dsx1), the male determining gene in Daphnia magna. By affinity purification, we identified an RNA binding protein, Shep as a DAPALR binding protein. Shep also binds to Dsx1 5´ UTR by recognizing the overlapping sequence and suppresses translation of the mRNA. In vitro and in vivo analyses indicated that DAPALR increased Dsx1 translation efficiency by sequestration of Shep. This regulation was impaired when the Shep binding site in DAPALR was deleted. These results suggest that Shep suppresses the unintentional translation of Dsx1 by setting a threshold; and when the sense lncRNA DAPALR is expressed, DAPALR cancels the suppression caused by Shep. This mechanism may be important to show dimorphic gene expressions such as sex determination and it may account for the binary expression in various developmental processes.

Long noncoding RNAs are vastly transcribed throughout the genome. Among them, RNAs overlapping the protein-coding RNA in sense orientation have been poorly studied because of the difficulty in differentiating their sequences from their overlapping coding RNAs although this class of RNAs has been reported to comprise the majority of the long noncoding RNAs. In the crustacean Daphnia magna, a long noncoding RNA, called DAPALR, is transcribed from the male determining gene, Doublesex1, and overlaps with the Doublesex1 5´ UTR. DAPALR activates Doublesex1 but this regulatory mechanism remains unknown. We found the RNA binding protein Shep bound to the Doublesex1 5´ UTR. In vitro and in vivo experiments indicated that Shep suppresses translation of the mRNA and DAPALR increases Doublesex1 translation efficiency by sequestration of Shep. Since male-specific expression of Doublesex1 is also regulated at the transcriptional level, we propose that Shep cancels the unexpected expression of Doublesex1 and maintains the feminized state for sexual dimorphism but DAPALR suppresses this repression by sequestration of Shep. We infer that this mechanism is not only for binary sex regulation but could function in the binary regulation of other genes in various biological processes.

A Crab Is Not a Fish: Unique Aspects of the Crustacean Endocrine System and Considerations for Endocrine Toxicology

Crustaceans-and arthropods in general-exhibit many unique aspects to their physiology. These include the requirement to moult (ecdysis) in order to grow and reproduce, the ability to change color, and multiple strategies for sexual differentiation. Accordingly, the endocrine regulation of these processes involves hormones, receptors, and enzymes that differ from those utilized by vertebrates and other non-arthropod invertebrates. As a result, environmental chemicals known to disrupt endocrine processes in vertebrates are often not endocrine disruptors in crustaceans; while, chemicals that disrupt endocrine processes in crustaceans are often not endocrine disruptors in vertebrates. In this review, we present an overview of the evolution of the endocrine system of crustaceans, highlight endocrine endpoints known to be a target of disruption by chemicals, and identify other components of endocrine signaling that may prove to be targets of disruption. This review highlights that crustaceans need to be evaluated for endocrine disruption with consideration of their unique endocrine system and not with consideration of the endocrine system of vertebrates.

Two Doublesex1 mutants revealed a tunable gene network underlying intersexuality in Daphnia magna

In recent years, the binary definition of sex is being challenged by repetitive reports about individuals with ambiguous sexual identity from various animal groups. This has created an urge to decode the molecular mechanism underlying sexual development. However, sexual ambiguities are extremely uncommon in nature, limiting their experimental value. Here, we report the establishment of a genetically modified clone of Daphnia magna from which intersex daphniids can be readily generated. By mutating the conserved central sex determining factor Doublesex1, body-wide feminization of male daphniid could be achieved. Comparative transcriptomic analysis also revealed a genetic network correlated with Doublesex1 activity which may account for the establishment of sexual identity in D. magna. We found that Dsx1 repressed genes related to growth and promoted genes related to signaling. We infer that different intersex phenotypes are the results of fluctuation in activity of these Dsx1 downstream factors. Our results demonstrated that the D. magna genome is capable of expressing sex in a continuous array, supporting the idea that sex is actually a spectrum.

A comprehensive epigenomic analysis of phenotypically distinguishable, genetically identical female and male Daphnia pulex

BACKGROUND

Daphnia species reproduce by cyclic parthenogenesis involving both sexual and asexual reproduction. The sex of the offspring is environmentally determined and mediated via endocrine signalling by the mother. Interestingly, male and female Daphnia can be genetically identical, yet display large differences in behaviour, morphology, lifespan and metabolic activity. Our goal was to integrate multiple omics datasets, including gene expression, splicing, histone modification and DNA methylation data generated from genetically identical female and male Daphnia pulex under controlled laboratory settings with the aim of achieving a better understanding of the underlying epigenetic factors that may contribute to the phenotypic differences observed between the two genders.

RESULTS

In this study we demonstrate that gene expression level is positively correlated with increased DNA methylation, and histone H3 trimethylation at lysine 4 (H3K4me3) at predicted promoter regions. Conversely, elevated histone H3 trimethylation at lysine 27 (H3K27me3), distributed across the entire transcript length, is negatively correlated with gene expression level. Interestingly, male Daphnia are dominated with epigenetic modifications that globally promote elevated gene expression, while female Daphnia are dominated with epigenetic modifications that reduce gene expression globally. For examples, CpG methylation (positively correlated with gene expression level) is significantly higher in almost all differentially methylated sites in male compared to female Daphnia. Furthermore, H3K4me3 modifications are higher in male compared to female Daphnia in more than 3/4 of the differentially regulated promoters. On the other hand, H3K27me3 is higher in female compared to male Daphnia in more than 5/6 of differentially modified sites. However, both sexes demonstrate roughly equal number of genes that are up-regulated in one gender compared to the other sex. Since, gene expression analyses typically assume that most genes are expressed at equal level among samples and different conditions, and thus cannot detect global changes affecting most genes.

CONCLUSIONS

The epigenetic differences between male and female in Daphnia pulex are vast and dominated by changes that promote elevated gene expression in male Daphnia. Furthermore, the differences observed in both gene expression changes and epigenetic modifications between the genders relate to pathways that are physiologically relevant to the observed phenotypic differences.

Phylogenetic analysis and embryonic expression of panarthropod Dmrt genes

Background

One set of the developmentally important Doublesex and Male-abnormal-3 Related Transcription factors (Dmrt) is subject of intense research, because of their role in sex-determination and sexual differentiation. This likely non-monophyletic group of Dmrt genes is represented by the Drosophila melanogaster gene Doublesex (Dsx), the Caenorhabditis elegans Male-abnormal-3 (Mab-3) gene, and vertebrate Dmrt1 genes. However, other members of the Dmrt family are much less well studied, and in arthropods, including the model organism Drosophila melanogaster, data on these genes are virtually absent with respect to their embryonic expression and function.

Results

Here we investigate the complete set of Dmrt genes in members of all main groups of Arthropoda and a member of Onychophora, extending our data to Panarthropoda as a whole. We confirm the presence of at least four families of Dmrt genes (including Dsx-like genes) in Panarthropoda and study their expression profiles during embryogenesis. Our work shows that the expression patterns of Dmrt11E, Dmrt93B, and Dmrt99B orthologs are highly conserved among panarthropods. Embryonic expression of Dsx-like genes, however, is more derived, likely as a result of neo-functionalization after duplication.

Conclusions

Our data suggest deep homology of most of the panarthropod Dmrt genes with respect to their function that likely dates back to their last common ancestor. The function of Dsx and Dsx-like genes which are critical for sexual differentiation in animals, however, appears to be much less conserved.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0322-0) contains supplementary material, which is available to authorized users.

Evolution of Gene Expression during a Transition from Environmental to Genetic Sex Determination

Genetic sex determination (GSD) can evolve from environmental sex determination (ESD) via an intermediate state in which both coexist in the same population. Such mixed populations are found in the crustacean Daphnia magna, where non-male producers (NMP, genetically determined females) coexist with male producers (MP), in which male production is environmentally inducible and can also artificially be triggered by exposure to juvenile hormone. This makes Daphnia magna a rare model species for the study of evolutionary transitions from ESD to GSD. Although the chromosomal location of the NMP-determining mutation has been mapped, the actual genes and pathways involved in the evolution of GSD from ESD remain unknown. Here, we present a transcriptomic analysis of MP and NMP females under control (female producing) and under hormone exposure conditions. We found ∼100 differentially expressed genes between MP and NMP under control conditions. Genes in the NMP-determining chromosome region were especially likely to show such constitutive expression differences. Hormone exposure led to expression changes of an additional ∼100 (MP) to ∼600 (NMP) genes, with an almost systematic upregulation of those genes in NMP. These observations suggest that the NMP phenotype is not determined by a simple "loss-of-function" mutation. Rather, homeostasis of female offspring production under hormone exposure appears to be an active state, tightly regulated by complex mechanisms involving many genes. In a broader view, this illustrates that the evolution of GSD, while potentially initiated by a single mutation, likely leads to secondary integration involving many genes and pathways.

Sensory Ecology of Predator-Induced Phenotypic Plasticity

Ecological communities are organized in trophic levels that share manifold interactions forming complex food webs. Infochemicals can further modify these interactions, e.g., by inducing defenses in prey. The micro-crustacean Daphnia is able to respond to predator-specific chemical cues indicating an increased predation risk. Daphnia shows plastic responses by adapting its morphology, behavior, and physiology, increasing organism, and population fitness. This stabilizes community structures. This review will describe the progress that has been made in understanding the high degree of plasticity observed in the model crustacean Daphnia. I summarize current knowledge on the processes of predator detection, ranging from the nature of biologically active chemical cues to the underlying neurophysiological mechanisms. With this, I aim to provide a comprehensive overview on the molecular mechanisms of ad hoc environmental phenotypic adaptation. In times of climate change and pollution understanding information transfer in aquatic systems is valuable as it will allow us to predict whether and how community structures are being affected.

Generation of white-eyed Daphnia magna mutants lacking scarlet function

The crustacean Daphnia magna is an important model in multi-disciplinary scientific fields such as genetics, evolutionary developmental biology, toxicology, and ecology. Recently, the draft genome sequence and transcriptome data became publicly available for this species. Genetic transformation has also been achieved via the introduction of plasmid DNA into the genome. The identification of a screenable marker gene and generation of mutant strains are essential to further advance D. magna functional genomics. Because crustaceans are closely related to insects, we hypothesized that, similar to Drosophila genetic studies, eye color-related genes can function as marker genes in Daphnia. We searched orthologs of Drosophila eye pigment transporters White, Scarlet, and Brown in the genome of D. magna. Amino acid sequence alignment and phylogenetic analysis suggested that D. magna has six white and one scarlet orthologs, but lacks the brown ortholog. Due to the multiplicity of white orthologs, we analyzed the function of the scarlet ortholog, DapmaSt, using RNA interference. DapmaSt RNAi embryos showed disappearance of black pigments both in the compound eye and in the ocellus, suggesting that DapmaSt is necessary for black pigmentation in Daphnia eyes. To disrupt DapmaSt using the Crispr/Cas9 system, we co-injected DapmaSt-targeting gRNAs with Cas9 mRNAs into eggs and established white-eyed DapmaSt mutant lines that lack eye pigments throughout their lifespan. Our results suggest that DapmaSt can be used as a transformation marker in D. magna and the DapmaSt mutants would be an important resource for genetic transformation of this species in the future.

Co-option of the bZIP transcription factor Vrille as the activator of Doublesex1 in environmental sex determination of the crustacean Daphnia magna

Divergence of upstream regulatory pathways of the transcription factor Doublesex (Dsx) serves as a basis for evolution of sex-determining mechanisms in animals. However, little is known about the regulation of Dsx in environmental sex determination. In the crustacean Daphnia magna, environmental sex determination is implemented by male-specific expression of the Dsx ortholog, Dsx1. Transcriptional regulation of Dsx1 comprises at least three phases during embryogenesis: non-sex-specific initiation, male-specific up-regulation, and its maintenance. Herein, we demonstrate that the male-specific up-regulation is controlled by the bZIP transcription factor, Vrille (Vri), an ortholog of the circadian clock genes-Drosophila Vri and mammalian E4BP4/NFIL3. Sequence analysis of the Dsx1 promoter/enhancer revealed a conserved element among two Daphnia species (D. magna and D. pulex), which contains a potential enhancer harboring a consensus Vri binding site overlapped with a consensus Dsx binding site. Besides non-sex-specific expression of Vri in late embryos, we found male-specific expression in early gastrula before the Dsx1 up-regulation phase begins. Knockdown of Vri in male embryos showed reduction of Dsx1 expression. In addition, transient overexpression of Vri in early female embryos up-regulated the expression of Dsx1 and induced male-specific trait. Targeted mutagenesis using CRISPR/Cas9 disrupted the enhancer on genome in males, which led to the reduction of Dsx1 expression. These results indicate that Vri was co-opted as a transcriptional activator of Dsx1 in environmental sex determination of D. magna. The data suggests the remarkably plastic nature of gene regulatory network in sex determination.