Changes in Induced-Antipredation Defense Traits and Transcriptome Regulations of Daphnia magna in Response to 5-HT1A Receptor Antagonist.

The serotonin signaling system plays a crucial role in regulating the ontogeny of crustaceans. Here, we describe the effects of different concentrations of the 5-hydroxytryptamine 1A receptor antagonist (WAY-100635) on the induced antipredation (Rhodeus ocellatus as the predator), morphological, behavioral, and life-history defenses of Daphnia magna and use transcriptomics to analyze the underlying molecular mechanisms. Our results indicate that exposure to WAY-100635 leads to changes in the expression of different defensive traits in D. magna when faced with fish predation risks. Specifically, as the length of exposure to WAY-100635 increases, high concentrations of WAY-100635 inhibit defensive responses associated with morphological and reproductive activities but promote the immediate negative phototactic behavioral defense of D. magna. This change is related to the underlying mechanism through which WAY-100635 interferes with gene expression of G-protein-coupled GABA receptors by affecting GABBR1 but promotes serotonin receptor signaling and ecdysteroid signaling pathways. In addition, we also find for the first time that fish kairomone can significantly activate the HIF-1α signaling pathway, which may lead to an increase in the rate of immediate movement. These results can help assess the potential impacts of serotonin-disrupting psychotropic drugs on zooplankton in aquatic ecosystems.

Non-negligible Toxicity to Fish in the Early Life Stages Triggered by Aqueous Leachate of Takeaway Plastic Containers.

Ordering takeout is a growing social phenomenon and may raise public health concerns. However, the associated health risk of compounds leaching from plastic packaging is unknown due to the lack of chemical and toxicity data. In this study, 20 chemical candidates were tentatively identified in the environmentally relevant leachate from plastic containers through the nontargeted chemical analysis. Three main components with high responses and/or predicted toxicity were further verified and quantified, namely, 3,5-di-tert-butyl-4-hydroxycinnamic acid (BHC), 2,4-di-tert-butylphenol (2,4-DTBP), and 9-octadecenamide (oleamide). The toxicity to zebrafish larvae of BHC, a degradation product of a widely used antioxidant Irganox 1010, was quite similar to that of the whole plastic leachate. In the same manner, RNA-seq-based ingenuity analysis showed that the affected canonical pathways of zebrafish larvae were quite comparable between BHC and the whole plastic leachate, i.e., highly relevant to neurological disease, metabolic disease, and even behavioral disorder. Longer-term exposure (35 days) did not cause any effect on adult zebrafish but led to decreased hatching rate and obvious neurotoxicity in zebrafish offspring. Collectively, this study strongly suggests that plastic containers can leach out a suite of compounds causing non-negligible impacts on the early stages of fish via direct or parental exposure.

Identification and characterization of homeobox gene clusters in harpacticoid and calanoid copepods.

In this study, we have identified the entire complement of typical homeobox (Hox) genes (Lab, Pb, Dfd, Scr, Antp, Ubx, Abd-A, and Abd-B) in harpacticoid and calanoid copepods and compared them with the cyclopoid copepod Paracyclopina nana. The harpacticoid copepods Tigriopus japonicus and Tigriopus kingsejongensis have seven Hox genes (Lab, Dfd, Scr, Antp, Ubx, Abd-A, and Abd-B) and the Pb and Ftz genes are also present in the cyclopoid copepod P. nana. In the Hox gene cluster of the calanoid copepod Eurytemora affinis, all the Hox genes were present linearly in the genome but the Antp gene was duplicated. Of the three representative copepods, the P. nana Hox gene cluster was the most compact due to its small genome size. The Hox gene expression profile patterns in the three representative copepods were stage-specific. The Lab, Dfd, Scr, Pb, Ftz, and Hox3 genes showed a high expression in early developmental stages but Antp, Ubx, Abd-A, and Abd-B genes were mostly expressed in later developmental stages, implying that these Hox genes may be closely associated with the development of segment identity during early development.

Multigenerational Mitigating Effects of Ocean Acidification on In Vivo Endpoints, Antioxidant Defense, DNA Damage Response, and Epigenetic Modification in an Asexual Monogonont Rotifer.

Ocean acidification (OA) is caused by changes in ocean carbon chemistry due to increased atmospheric pCO2 and is predicted to have deleterious effects on marine ecosystems. While the potential impacts of OA on many marine species have been studied, the multigenerational effects on asexual organisms remain unknown. We found that low seawater pH induced oxidative stress and DNA damage, decreasing growth rates, fecundity, and lifespans in the parental generation, whereas deleterious effects on in vivo endpoints in F1 and F2 offspring were less evident. The findings suggest that multigenerational adaptive effects play a role in antioxidant abilities and other defense mechanisms. OA-induced DNA damage, including double-strand breaks (DSBs), was fully repaired in F1 offspring of parents exposed to OA for 7 days, indicating that an adaptation mechanism may be the major driving force behind multigenerational adaptive effects. Analysis of epigenetic modification in response to OA involved examination of histone modification of DNA repair genes and a chromatin immunoprecipitation assay, as Bombus koreanus has no methylation pattern for CpG in its genome. We conclude that DSBs, DNA repair, and histone modification play important roles in multigenerational plasticity in response to OA in an asexual monogonont rotifer.

Differential gene expression profile of male and female copepods in response to cadmium exposure.

In this study, the whole transcriptome and sex-specific differential gene expression of the copepod Pseudodiaptomus annandalei exposed to cadmium (Cd) were investigated. P. annandalei were exposed to 40 µg/L Cd from the naupliar stage to male and female adults. High-throughput transcriptome sequencing (RNA-seq) was performed with copepod samples using an Illumina Hiseq™ 2000 platform. TransDecoder analysis found 32,625 putative open reading frame contigs. At p-values of <0.001, a total of 4756 differentially expressed genes (DEGs) (2216 up-regulated and 2540 down-regulated genes) were found in male copepods. Whereas a total of 2879 DEGs (2007 up-regulated and 872 down-regulated genes) were found in female copepods. A few selected cellular stress response genes, involved in xenobiotic metabolism, energy metabolism, growth, and development as a result of Cd exposure in the copepods were discussed. The study showed that most of these processes were changed in a sex-specific manner, accounting for the different sensitivities of male and female copepods. Results suggest and reinforce that sex is an important factor to be considered in ecotoxicogenomics.

Nanoplastic Ingestion Enhances Toxicity of Persistent Organic Pollutants (POPs) in the Monogonont Rotifer Brachionus koreanus via Multixenobiotic Resistance (MXR) Disruption.

Among the various materials found inside microplastic pollution, nanosized microplastics are of particular concern due to difficulties in quantification and detection; moreover, they are predicted to be abundant in aquatic environments with stronger toxicity than microsized microplastics. Here, we demonstrated a stronger accumulation of nanosized microbeads in the marine rotifer Brachionus koreanus compared to microsized ones, which was associated with oxidative stress-induced damages on lipid membranes. In addition, multixenobiotic resistance conferred by P-glycoproteins and multidrug resistance proteins, as a first line of membrane defense, was inhibited by nanoplastic pre-exposure, leading to enhanced toxicity of 2,2',4,4'-tetrabromodiphenyl ether and triclosan in B. koreanus. Our study provides a molecular mechanistic insight into the toxicity of nanosized microplastics toward aquatic invertebrates and further implies the significance of synergetic effects of microplastics with other environmental persistent organic pollutants.

A review of environmental epigenetics in aquatic invertebrates.

Aquatic ecosystems face significant challenges due to increasing human-induced environmental stressors. Recent studies emphasize the role of epigenetic mechanisms in the stress responses and adaptations of organisms to those stressors. Epigenetics influences gene expression, enabling phenotypic plasticity and transgenerational effects. Therefore, understanding the epigenetic responses of aquatic invertebrates to environmental stressors is imperative for aquatic ecosystem research. In this study, we organize the mechanisms of epigenetics in aquatic invertebrates and explore their roles in the responses of aquatic invertebrates to environmental stressors. Furthermore, we discuss the inheritance of epigenetic changes and their influence across generations in aquatic invertebrates. A comprehensive understanding of epigenetic responses is crucial for long-term ecosystem management and conservation strategies in the face of irreversible climate change in aquatic environments. In this review, we synthesize existing knowledge about environmental epigenetics in aquatic invertebrates to provide insights and suggest directions for future research.

In silico identification and characterization of microRNAs from rotifers, cladocerans, and copepods.

MicroRNAs (miRNAs) are short non-coding RNA molecules that regulate post-transcription and influence various biological processes across species. Despite various studies of miRNAs in vertebrates, plants, and other organisms, miRNA data in aquatic invertebrates are insufficient. In this study, we identified miRNAs from four aquatic invertebrate species that are widely used in aquatic toxicology: the rotifer Brachionus koreanus, the water flea Daphnia magna, the cyclopoid copepod Paracyclopina nana, and the harpacticoid copepod Tigriopus japonicus, using next-generation sequencing and in silico analysis. We identified total 188, 41, 47, and 100 miRNAs from each species, and target genes were predicted based on 3'-untranslated region information. Target prediction and functional annotation results provided the biological processes of these miRNAs in various development-related mechanisms, signaling transduction, and metabolism-related pathways. Moreover, the network between the miRNAs and their targets concerning defense-related and antioxidant genes suggests the suitability of miRNAs as biomarkers in ecotoxicological studies.

Genome-wide identification of 769 G protein-coupled receptor (GPCR) genes from the marine medaka Oryzias melastigma.

The marine medaka Oryzias melastigma is a useful fish model for marine and estuarine ecotoxicology studies and can be applied to field-based population genomics because of its distribution in Asian estuaries and other coastal areas. We identified 769 full-length G protein-coupled receptor genes in the O. melastigma genome and classified them into five distinct classes. A phylogenetic comparison of GPCR genes in O. melastigma to humans and two other small fish species revealed a high-level orthological relationship. Purinergic and chemokine receptors were highly differentiated in humans whereas significant differentiation of chemosensory receptors was evident in fish species. Our results suggest that the GPCR gene families among the species used in this study exhibit evidence of sporadic evolutionary processes. These results may help improve our understanding of the advanced repertoires of GPCR and expand our knowledge of physiological mechanisms of fish in response to various environmental stimuli.

Multigenerational effects of elevated temperature on host-microbiota interactions in the marine water flea Diaphanosoma celebensis exposed to micro- and nanoplastics.

Rising ocean temperatures are driving unprecedented changes in global marine ecosystems. Meanwhile, there is growing concern about microplastic and nanoplastic (MNP) contamination, which can endanger marine organisms. Increasing ocean warming (OW) and plastic pollution inevitably cause marine organisms to interact with MNPs, but relevant studies remain sparse. Here, we investigated the interplay between ocean warming and MNP in the marine water flea Diaphanosoma celebensis. We found that combined exposure to MNPs and OW induced reproductive failure in the F2 generation. In particular, the combined effects of OW and MNPs on the F2 generation were associated with key genes related to reproduction and stress response. Moreover, populations of predatory bacteria were significantly larger under OW and MNP conditions during F2 generations, suggesting a potential link between altered microbiota and host fitness. These results were supported by a host transcriptome and microbiota interaction analysis. This research sheds light on the complex interplay between environmental stressors, their multigenerational effects on marine organisms, and the function of the microbiome.

Toxicity and speciation of inorganic arsenics and their adverse effects on in vivo endpoints and oxidative stress in the marine medaka Oryzias melastigma.

Here, we investigate the effects of acute and chronic exposure to arsenate (AsV) and arsenite (AsIII) in the marine medaka Oryzias melastigma. In vivo effects, biotransformation, and oxidative stress were studied in marine medaka exposed to the two inorganic arsenics for 4 or 28 days. An investigation of embryonic development revealed no effect on in vivo parameters, but the hatching rate increased in the group exposed to AsIII. Exposure to AsIII also caused the greatest accumulation of arsenic in medaka. For acute exposure, the ratio of AsV to AsIII was higher than that of chronic exposure, indicating that bioaccumulation of inorganic arsenic can induce oxidative stress. The largest increase in oxidative stress was observed following acute exposure to AsIII, but no significant degree of oxidative stress was induced by chronic exposure. During acute exposure to AsV, the increase in the enzymatic activity of glutathione-S-transferase (GST) was twice as high compared with exposure to AsIII, suggesting that GST plays an important role in the initial detoxification process. In addition, an RNA-seq-based ingenuity pathway analysis revealed that acute exposure to AsIII may be related to cell-cycle progression. A network analysis using differentially expressed genes also revealed a potential link between the generation of inflammatory cytokines and oxidative stress due to arsenic exposure.

Synergistic adverse effects of microfibers and freshwater acidification on host-microbiota interactions in the water flea Daphnia magna.

Microfibers are the most common type of microplastics in freshwater environments. Anthropogenic climate stressors, such as freshwater acidification (FA), can interact with plastic pollution to disrupt freshwater ecosystems. However, the underlying mechanisms responsible for the interactive effects of microfibers and FA on aquatic organisms remain poorly understood. In this study, we investigated individual Daphnia magna-microbiota interactions affected by interactions between microfibers and FA (MFA). We found that the accumulated amount of microfibers in pH-treatment groups was significantly higher than in the control groups, resulting in negative consequences on reproduction, growth, and sex ratio. We also observed that MFA interactions induced immunity- and reproduction-related biological processes. In particular, the abundance of pathogenic bacteria increased only in MFA groups, indicating that MFA interactions can cause intestinal damage. Our integrated analysis of microbiomes and host transcriptomes revealed that synergistic adverse effects of MFAs are closely related to changes in microbial communities, suggesting that D. magna fitness and the microbial community are causally linked. These finding may help elucidate the toxicity mechanisms governing the responses of D. magna to microfibers and acidification interactions, and to host-microbiome-environment interactions.

Nanoplastics induce epigenetic signatures of transgenerational impairments associated with reproduction in copepods under ocean acidification.

Ocean acidification (OA) is one of many major global climate changes that pose a variety of risks to marine ecosystems in different ways. Meanwhile, there is growing concern about how nanoplastics (NPs) affect marine ecosystems. Combined exposure of marine organisms to OA and NPs is inevitable, but their interactive effects remain poorly understood. In this study, we investigated the multi- and transgenerational toxicity of NPs on copepods under OA conditions for ten generations. The findings revealed that OA and NPs have a synergistic negative effect on copepod reproduction across generations. In particular, the transgenerational groups showed reproductive impairments in the F1 and F2 generations (F1T and F2T), even though they were never exposed to NPs. Moreover, our epigenetic examinations demonstrated that the observed intergenerational reproductive impairments are associated with differential methylation patterns of specific genes, suggesting that the interaction of OA and NPs can pose a significant threat to the sustainability of copepod populations through epigenetic modifications. Overall, our findings provide valuable insight into the intergenerational toxicity and underlying molecular mechanisms of responses to NPs under OA conditions.

Oxidative stress-mediated synergistic deleterious effects of nano- and microplastics in the hypoxia-conditioned marine rotifer Brachionus plicatilis.

While pollution due to nano- and micro-sized plastics (NMPs) and hypoxic conditions both occur in coastal areas, the deleterious potential of co-exposure to hypoxia and NMPs (hypoxia and micro-sized plastics, HMPs; hypoxia and nano-sized plastics, HNPs) is largely unclear. Here, we provide evidence for multigenerational effects of HMP and HNP in the marine rotifer Brachionus plicatilis by investigating changes in its life traits, antioxidant system, and hypoxia-inducible factor (HIF) pathway using an orthogonal experimental design, with nanoscale and microscale particles measuring 0.05 µm and 6.0 µm in diameter, respectively, and hypoxic conditions of 0.5 mg/L for six generations. Combined exposure to NMPs and hypoxia caused a significant decrease in fecundity and overproduction of reactive oxygen species (ROS). The HIF pathway and circadian clock genes were also significantly upregulated in response to HMP and HNP exposure. In particular, synergistic deleterious effects of HNP were evident, suggesting that size-dependent toxicity can be a major driver of the effects of hypoxia and NMP co-exposure. After several generations of exposure, ROS levels returned to basal levels and transcriptomic resilience was observed, although rotifer reproduction remained suppressed. These findings help eluciating the underlying molecular mechanisms involved in responses to plastic pollution in hypoxic conditions.

The Genome of the Marine Rotifer Brachionus manjavacas: Genome-Wide Identification of 310 G Protein-Coupled Receptor (GPCR) Genes.

The marine rotifer Brachionus manjavacas is widely used in ecological, ecotoxicological, and ecophysiological studies. The reference genome of B. manjavacas is a good starting point to uncover the potential molecular mechanisms of responses to various environmental stressors. In this study, we assembled the whole-genome sequence (114.1 Mb total, N50 = 6.36 Mb) of B. manjavacas, consisting of 61 contigs with 18,527 annotated genes. To elucidate the potential ligand-receptor signaling pathways in marine Brachionus rotifers in response to environmental signals, we identified 310 G protein-coupled receptor (GPCR) genes in the B. manjavacas genome after comparing them with three other species, including the minute rotifer Proales similis, Drosophila melanogaster, and humans (Homo sapiens). The 310 full-length GPCR genes were categorized into five distinct classes: A (262), B (26), C (7), F (2), and other (13). Most GPCR gene families showed sporadic evolutionary processes, but some classes were highly conserved between species as shown in the minute rotifer P. similis. Overall, these results provide potential clues for in silico analysis of GPCR-based signaling pathways in the marine rotifer B. manjavacas and will expand our knowledge of ligand-receptor signaling pathways in response to various environmental signals in rotifers.

Molecular evidence for suppression of swimming behavior and reproduction in the estuarine rotifer Brachionus koreanus in response to COVID-19 disinfectants.

The increased use of disinfectants due to the spread of the novel coronavirus infection (e.g. COVID-19) has caused burden in the environment but knowledge on its ecotoxicological impact on the estuary environment is limited. Here we report in vivo and molecular endpoints that we used to assess the effects of chloroxylenol (PCMX) and benzalkonium chloride (BAC), which are ingredients in liquid handwash, dish soap products, and sanitizers used by consumers and healthcare workers on the estuarine rotifer Brachionus koreanus. PCMX and BAC significantly affected the life table parameters of B. koreanus. These chemicals modulated the activities of antioxidant enzymes such as superoxide dismutase and catalase and increased reactive oxygen species even at low concentrations. Also, PCMX and BAC caused alterations in the swimming speed and rotation rate of B. koreanus. Furthermore, an RNA-seq-based ingenuity pathway analysis showed that PCMX affected several signaling pathways, allowing us to predict that a low concentration of PCMX will have deleterious effects on B. koreanus. The neurotoxic and mitochondrial dysfunction event scenario induced by PCMX reflects the underlying molecular mechanisms by which PCMX produces outcomes deleterious to aquatic organisms.

Comparative genome analysis of the monogonont marine rotifer Brachionus manjavacas Australian strain: Potential application for ecotoxicology and environmental genomics.

This is the first study to analyze the whole-genome sequence of B. manjavacas Australian (Aus.) strain through combination of Oxford Nanopore long-read seq, resulting in a total length of 108.1 Mb and 75 contigs. Genome-wide detoxification related gene families in B. manjavacas Aus. strain were comparatively analyzed with those previously identified in other Brachionus spp., including B. manjavacas German (Ger.) strain. Most of the subfamilies in detoxification related families (CYPs, GSTs, and ABCs) were highly conserved and confirmed orthologous relationship with Brachionus spp., and with accumulation of genome data, clear differences between genomic repertoires were demonstrated the marine and the freshwater species. Furthermore, strain-specific genetic variations were present between the Aus. and Ger. strains of B. manjavacas. This whole-genome analysis provides in-depth review on the genomic structural differences for detoxification-related gene families and further provides useful information for comparative ecotoxicological studies and evolution of detoxification mechanisms in Brachionus spp.

Local lymph node assay: 5-bromo-d-deoxyuridine-ELISA method for comparative study in assessing chemical potencies and skin sensitization in BALB/c and CBA/J strains.

Local lymph node assay (LLNA) is a predictive in vivo method to provide estimates of relative potency and to contribute to risk assessment/risk management regarding skin sensitizing potency of chemicals and formulations as a stand-alone alternative test. In addition, LLNA is relatively rapid and cost-effective compared to the Buehler method (Guinea pig test), and confers important animal welfare benefits. CBA/J and BALB/c strains are widely commercially available and have been evaluated by formal LLNA validation studies. However, the LLNA method using BrdU with ELISA, unlike other LLNA methods (OECD TG 429, 442 A, 442B), has not been previously validated. Therefore, in this study a validation method was performed to evaluate if the LLNA:BrdU-ELISA method could also be used to identify sensitizers among chemicals listed in OECD TG 429 using CBA/J and BALB/c strains. Here, we newly found that the LLNA:BrdU-ELISA validation method correctly identified 12 of 13 sensitizers in the BALB/c, 11 of 13 sensitizers in the CBA/J, and 3 of 5 non-sensitizers were identified in the two strains. Collectively, we found that the results of LLNA:BrdU-ELISA method provide a similar level of performance for accuracy and sensitivity in two mouse strains BALB/c and CBA/J.

The genome of the freshwater monogonont rotifer Brachionus rubens: Identification of phase I, II, and III detoxification genes.

Monogonont rotifers are common species in aquatic environments and make model species for ecotoxicology studies. Whole genomes of several species of the genus Brachionus have been assembled, but no information on the freshwater rotifer Brachionus rubens has been reported. In this study, the whole-genome sequence of B. rubens was successfully assembled using NextDenovo. The total length of the genome was 132.7 Mb (N50 = 2.51 Mb), including 122 contigs. The GC contents accounted for 29.96% of the genome. Aquatic organisms are always exposed to various external stresses, and a comprehensive genomic analysis is needed to better understand the adverse effects on organisms. This paper focuses on the ecotoxicological aspect and conducted genome analysis of representative gene families involved in detoxification mechanisms against environmental stressors. Specifically, we identified cytochrome P450 genes (CYPs) of phase I, glutathione S-transferase genes (GSTs) of phase II, and ATP-binding cassette transporter genes (ABCs) of phase III in the genome of B. rubens. Gene duplications were found in CYP, GST, and ABC genes, as is the case for other Brachionus rotifers. Our results suggest that these detoxification-related gene families have evolved in a species-specific and/or lineage-specific manner. This paper improves our understanding of how the freshwater Brachionus rotifers respond to environmental stressors in a molecular ecotoxicology context.

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.