De novo transcriptome assembly of brackish water flea Diaphanosoma celebensis based on short-term cadmium and benzo[a]pyrene exposure experiments

Oxidative stress-mediated deleterious effects of hypoxia in the brackish water flea Diaphanosoma celebensis

In this study, we investigated the acute toxicity, in vivo effects, oxidative stress, and gene expression changes caused by hypoxia on the brackish water flea Diaphanosoma celebensis. The no-observed-effect concentration (NOEC) of 48 h of hypoxia exposure was found to be 2 mg/L O2. Chronic exposure to NOEC caused a significant decline in lifespan but had no effect on total fecundity. The induction of reactive oxygen species increased in a time-dependent manner over 48 h, whereas the content of antioxidant enzymes (superoxide dismutase and catalase) decreased. The transcription and translation levels were modulated by hypoxia exposure. In particular, a significant increase in hemoglobin level was followed by up-regulation of hypoxia-inducible factor 1α gene expression and activation of the mitogen-activated protein kinase pathway. In conclusion, our findings provide a better understanding of the molecular mechanism of the adverse effects of hypoxia in brackish water zooplankton.

Pathway-dependent toxic interaction between polystyrene microbeads and methylmercury on the brackish water flea Diaphanosoma celebensis: Based on mercury bioaccumulation, cytotoxicity, and transcriptomic analysis

Given their worldwide distribution and toxicity to aquatic organisms, methylmercury (MeHg) and microplastics (MP) are major pollutants in marine ecosystems. Although they commonly co-exist in the ocean, information on their toxicological interactions is limited. Therefore, to understand the toxicological interactions between MeHg and MP (6-μm polystyrene), we investigated the bioaccumulation of MeHg, its cytotoxicity, and transcriptomic modulation in the brackish water flea Diaphanosoma celebensis following single and combined exposure to MeHg and MP. After single exposure to MeHg for 48-h, D. celebensis showed high Hg accumulation (34.83 ± 0.40 μg/g dw biota) and cytotoxicity, which was reduced upon co-exposure to MP. After transcriptomic analysis, 2, 253, and 159 differentially expressed genes were detected in the groups exposed to MP, MeHg, and MeHg+MP, respectively. Genes related to metabolic pathways and the immune system were significantly affected after MeHg exposure, but the effect of MeHg on these pathways was alleviated by MP co-exposure. However, MeHg and MP exhibited synergistic effects on the expression of gene related to DNA replication. These findings suggest that MP can reduce the toxicity of MeHg but that their toxicological interactions differ depending on the molecular pathway.

Metabolomic analysis of combined exposure to microplastics and methylmercury in the brackish water flea Diaphanosoma celebensis

Owing to their widespread distribution and high bioaccumulation, microplastics (MPs) and mercury (Hg) are considered major threats to the ocean. MP interacts with Hg because of its high adsorption properties. However, their toxicological interactions with marine organisms, especially combined effects at the molecular level, are poorly understood. This study investigated the single and combined effects of MP and Hg on the metabolic profile of the brackish water flea Diaphanosoma celebensis. A total of 238 metabolites were significantly affected by MP, Hg, or MP + Hg. Metabolite perturbation patterns showed that toxicity of Hg and MP + Hg was similar and that of MP was not significant. Among the 223 metabolites affected by Hg, profiles of 32 unannotated metabolites were significantly different from those of MP + Hg, and combined effects of MP + Hg decreased the effect of Hg on 25 of these metabolites. Only 11 annotated metabolites were significantly affected by Hg or MP + Hg and were related to carbohydrate, lipid, vitamin, and ecdysteroid metabolism. Ten metabolites were decreased by Hg and MP + Hg and were not significantly different between the exposure groups. Enrichment analysis showed that galactose, starch, and sucrose metabolism were the most affected pathways. These findings suggest that MP has negligible toxic effect, and Hg can induce energy depletion, membrane damage, and disruption of growth, development, and reproduction. Although the impact of MP was negligible, the combined effects of MP + Hg could be metabolite specific. This study provides better understanding of the combined effects of MP and Hg on marine organisms.

Effect of heavy metals on the energy metabolism in the brackish water flea Diaphanosoma celebensis

Heavy metals such as lead (Pb), cadmium (Cd), and arsenic (As) are of great concern in aquatic ecosystems because of their global distribution, persistence, and biomagnification via the food web. They can induce the expression of cellular protective systems (e.g., detoxification enzymes and antioxidant enzymes) to protect organisms from oxidative stress, which is a high-energy-consuming process. Thus, energy reserves (e.g., glycogen, lipids, and proteins) are utilized to maintain metabolic homeostasis. Although a few studies have suggested that heavy metal stress can modulate the metabolic cycle in crustaceans, information on changes in energy metabolism under metal pollution remains lacking in planktonic crustaceans. In the present study, the activity of digestive enzymes (amylase, trypsin, and lipase) and the contents of energy storage molecules (glycogen, lipid, and protein) were examined in the brackish water flea Diaphanosoma celebensis exposed to Cd, Pb, and As for 48 h. Transcriptional modulation of the three AMP-activated protein kinase (AMPK) and metabolic pathway-related genes was further investigated. Amylase activity was highly increased in all heavy metal-exposed groups, whereas trypsin activity was reduced in Cd- and As-exposed groups. While glycogen content was increased in all exposed groups in a concentration-dependent manner, lipid content was reduced at higher concentrations of heavy metals. The expression of AMPKs and metabolic pathway-related genes was distinct among heavy metals. In particular, Cd activated the transcription of AMPK-, glucose/lipid metabolism-, and protein synthesis-related genes. Our findings indicate that Cd can disrupt energy metabolism, and may be a potent metabolic toxicant in D. celebensis. This study provides insights into the molecular mode of action of heavy metal pollution on the energy metabolism in planktonic crustaceans.

De Novo Transcriptomic and Life-History Responses of Moina Micrura Under Stress Environment Conditions

Moina micrura represents a promising model species for ecological and ecotoxicological investigations in tropical freshwater ecosystems. Illumina NovaSeq™ 6000 sequencing was employed in this study to analyze M. micrura across three distinct developmental stages: juvenile, adult, and male. Current study successfully annotated 51,547 unigenes (73.11%) derived from seven (7) different databases. A total of 554 genes were found to be significantly upregulated, while 452 genes showed significant downregulation between juvenile and male. Moreover, 1001 genes were upregulated, whereas 830 genes exhibited downregulation between the adult and male. Analysis of differentially expressed genes revealed upregulation of chitin, cuticle, myosin (MYO), mitogen-activated protein kinases (MAPK), fibrillin (FBN), cytochrome (CYP), glutathione s-transferase (GST), vitellogenin (VTG), acetylcholinesterase (AChE), and transforming growth factor beta (TGFB) under unfavorable environmental conditions (male), as compared to favorable environmental conditions (juveniles and adults). These alterations in gene expression significantly impact the phenological and life-history traits of M. micrura. Furthermore, the upregulation of hemoglobin (HMB), doublesex (DSX), juvenile hormone analogs (JHA), heat shock protein (HSP), and methyltransferase (METT) genes in males initiates the sex-switching effects observed in M. micrura. These findings hold substantial value for researchers interested in determining M. micrura sequences for future investigations of gene expression and comparative reproductive genome analysis within the Moina genus and cladoceran families.

Modulation of Cd and BaP uptake rate during acute aqueous co-exposure in adult zebrafish (Danio rerio)

Cadmium and Benzo[a]pyrene are two toxicants of great environmental importance given their frequency and ability to cause extensive toxicity in aquatic organisms including fish. There is evidence that fish can modulate their respective uptake rate during simultaneous exposures, albeit the mechanism behind this is poorly understood. The present study aimed to examine this interaction by exposing adult zebrafish to either 89.3 nM Cd, 4.25 nM BaP or a combination of the two for 72 hrs prior to examining the uptake rate of either toxicant via short-term exposures (3-6 hrs) to radiotracers (¹⁰⁹Cd and ¹⁴C-BaP). Our results showed that Cd uptake rate increased significantly in the gills when animals were pre-exposed to both toxicants simultaneously, resulting in an increased maximum uptake rate (Jmax). The increased Cd uptake rate did not correspond to increased expression of gill Cd transporters such as the epithelium calcium channel (ECaC) or the divalent metal transporter 1 (DMT1). Furthermore, BaP uptake rate increased significantly at the whole-body level when animals were exposed to both 5.03 nM ¹⁴C-BaP and 89.3 nM Cd concurrently. Additionally, we ran a time-course and observed BaP uptake rate is highest in the 6-12 hrs following the beginning of the exposure. Our results provide evidence that the increased bioaccumulation of Cd and BaP observed during co-exposures is at least in part due to an increase in uptake rate and is driven by separate mechanisms.

Modulation of ecdysteroid and juvenile hormone signaling pathways by bisphenol analogues and polystyrene beads in the brackish water flea Diaphanosoma celebensis

Owing to its high production and world-wide usage, plastic pollution is an increasing concern in marine environments. Plastic is decomposed into nano- and micro-sized debris, which negative affect reproduction and development in aquatic organisms. Bisphenol A (BPA), an additive of plastic, is released into the water column upon plastic degradation, and is known as a representative endocrine-disrupting chemical. However, the reproductive effects of plastics and bisphenols at the molecular level have not yet been explored in small marine crustaceans. In this study, we investigated the effects of polystyrene (PS) beads (0.05, 0.5, and 6 - μm) and bisphenol analogues (BPs; BPA, BPS, and BPF) on reproduction and development of small marine crustaceans. Effects on transcriptional changes in ecdysteroid and juvenile hormone (JH) signaling pathway-related genes were examined in the brackish water flea Diaphanosoma celebensis exposed to PS beads and BPs for 48 h. As results, BPs and PS beads delayed emergence time of first offspring, and increased fecundity in a concentration-dependent manner. BPs differentially modulated the expression of ecdysteroid and JH signaling pathway-related genes, indicating that BP analogs can disrupt endocrine systems via mechanisms different from those of BPA. PS beads was also changed the gene expression of both pathway, depending on their size and concentration. Our findings suggest that BP analogues and PS beads disrupt the endocrine system by modulating the hormonal pathways, affecting reproduction negatively. This study provides a better understanding of the molecular mode of action of BPs and PS beads in the reproduction of small crustaceans.

Validation of reference genes for quantitative real-time PCR in chemical exposed and at different age's brackish water flea Diaphanosoma celebensis

Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), a primary approach for evaluating gene expression, requires an appropriate normalization strategy to confirm relative gene expression levels by comparison, and rule out variations that might occur in analytical procedures. The best option is to use a reference gene whose expression level is stable across various experimental conditions to compare the mRNA levels of a target gene. However, there is limited information on how the reference gene is differentially expressed at different ages (growth) in small invertebrates with notable changes such as molting. In this study, expression profiles of nine candidate reference genes from the brackish water flea, Diaphanosoma celebensis, were evaluated under diverse exposure to toxicants and according to growth. As a result, four different algorithms showed similar stabilities of genes for chemical exposures in the case of limited conditions using the same developmental stage (H2A was stable, whereas Act was fairly unstable in adults), while the results according to age showed a significantly different pattern in suite of candidate reference genes. This affected the results of genes EcRA and GST, which are involved in development and detoxification mechanisms, respectively. Our finding is the first step towards establishing a standardized real-time qRT-PCR analysis of this environmentally important invertebrate that has potential for aquatic ecotoxicology, particularly in estuarine environments.

Genome-wide identification of 216 G protein-coupled receptor (GPCR) genes from the marine water flea Diaphanosoma celebensis

G protein-coupled receptors (GPCRs) are considered to have originated from early evolution of eukaryotic species, therefore, the genome-wide identification of GPCR genes can provide insight into the adaptive strategy and evolutionary tendency in an animal taxon. Here, we identified a total 216 full-length GPCR genes in the marine water flea Diaphanosoma celebensis genome, which were classified into five distinct classes (A, B, C, F, and other). Phylogenetic comparison of GPCRs in D. celebensis to those in humans (Homo sapiens), fruitfly (Drosophila melanogaster), and freshwater water flea (Daphnia magna) reveals a high level of orthological relationship of amine, neuropeptide, and opsin receptor repertoire, while purinergic and chemokine receptors were highly differentiated in humans. Our findings suggest sporadic evolutionary processes within the GPCR gene families identified in D. celebensis. In this study, these results may provide a better understanding on the evolution of GPCRs, and expand our knowledge of the cladoceran GPCR gene repertories which in part, mediate cell physiological mechanisms in response to various environmental stimuli.

The genome of the marine water flea Diaphanosoma celebensis: Identification of phase I, II, and III detoxification genes and potential applications in marine molecular ecotoxicology

To assemble the genome of the marine water flea Diaphanosoma celebensis, a sentinel model for marine environmental monitoring, we constructed a high-quality genome using PromethION and HiSeq 2500 platforms. The total length of the assembled genome was 100.08 Mb, with N50 = 2.56 Mb (benchmarking universal single-copy orthologs, 96.9%) and consisted of 179 scaffolds. A total of 15,427 genes were annotated, and orthologous gene clusters in D. celebensis were analyzed and compared with those of the cladocerans Daphnia magna and Daphnia pulex. In addition, phase I, II, and III detoxification gene families of cytochrome P450s, glutathione S-transferases, and ATP-binding cassette were fully identified and revealed lineage-specific gene loss and/or expansion, suggesting that the evolution of detoxification gene families likely modulates fitness and susceptibility in response to environmental stressors. The study improves our understanding of the detoxification-related gene system and should contribute to future studies of molecular ecotoxicology in cladoceran species and their responses to emerging pollutants.

Identification and response of cytochrome P450 genes in the brackish water flea Diaphanosoma celebensis after exposure to benzo[α]pyrene and heavy metals

The cytochrome P450 (CYP) enzyme family is extensive; these enzymes participate in phase I enzyme metabolism and are involved in xenobiotic detoxification in all living organisms. Despite their significance in xenobiotic detoxification, little is known about the species-specific comparison of CYPs and their molecular responses in aquatic invertebrates. We identified 31 CYPs in the brackish water flea Diaphanosoma celebensis via thorough exploration of transcriptomic databases and measured the transcript profiles of 9 CYPs (within full sequences) in response to benzo[α]pyrene (B[α]P) and two heavy metals (cadmium [Cd] and copper [Cu]). Through phylogenetic analysis, the CYPs were separated and clustered into four clans: mitochondrial, CYP2, CYP3, and CYP4. The expression of 9 CYPs were differentially modulated (up- and/or downregulated) in response to B[α]P, Cd, and Cu. In particular, CYP370A15 was significantly upregulated in response to B[α]P, Cd, and Cu, suggesting that the identified CYPs are involved in xenobiotic detoxification and are useful as biomarkers in response to B[α]P, Cd, and Cu. This study aimed to comprehensively annotate cladoceran CYPs; our results will add to the existing knowledge on the potential roles of CYPs in xenobiotic detoxification in cladocerans.

Oxidative stress responses in brackish water flea exposed to microcystin-LR and algal bloom waters from Nakdong River, Republic of Korea

Microcystis blooms and the impact of their toxins, particularly microcystin (MC), in coastal ecosystems is an emerging threat, but the species-specific effects of MC and the potential for bioconcentration are not fully understood. We exposed the brackish water flea, Diaphanosoma celebensis, to MC-LR, which showed antioxidant responses measured at the molecular to enzyme levels but no acute toxicity. We extended our experimental investigation to measure the released MC and its uptake by D. celebensis exposed to river water. In a short-term exposure (48 h) experiment, D. celebensis exposed to water from an algal bloom (approximately 2 μg L^-1 MC) assimilated more than 50 pg MC per individual. The significant increase of MCs suggests the potential for the species to accumulate MCs. The dose-dependent increase in the antioxidant response observed in the mRNA levels also showed that D. celebensis exposed to diluted algal bloom waters were affected by toxins from cyanobacteria.

Complete mitochondrial DNA of the marine water flea Diaphanosoma celebensis (Cladocera, Sididae)

The complete mitochondrial genome was sequenced from the marine water flea Diaphanosoma celebensis. The sequenced mitochondrial genome size was 17,060 bp, possessing identical gene order of 13 protein-coding genes (PGCs) to those of the congeneric freshwater species Diaphanosoma dubium in the genus Diaphanosoma. The mitochondrial genome of D. celebensis had 13 PGCs, two rRNAs, and 22 tRNAs. Of 13 PGCs, three genes (CO3, ND3, and ND4) had incomplete stop codons. Furthermore, the stop codons of the remaining ten PGCs were TAA (for CO1, ATP8, ATP6, ND5, ND6, and ND1) and TAG (for NL4L, Cytb, and ND2). The second and third base composition of codon on 9 PCGs on the L strand in D. celebensis mitogenome showed an anti-G bias (11.0% and 15.0%), respectively.

Transcriptomic changes in western mosquitofish (Gambusia affinis) liver following benzo[a]pyrene exposure

Widely distributed western mosquitofish (Gambusia affinis) has been used as a new model species for hazard assessment of environmental stressors such as polycyclic aromatic hydrocarbons (PAHs). However, most of the PAH studies using G. affinis rely on targeted biomarker-based analysis, and thus may not adequately address the complexity of the toxic mechanisms of the stressors. In the present study, the whole transcriptional sequencing of G. affinis liver after exposure to a PAH model, benzo[a]pyrene (BaP) (100 μg/L), for 20 days was performed by using the HiSeq XTen sequencers. In total, 58,156,233 and 51,825,467 clean nucleotide reads were obtained in the control and BaP-exposed libraries, respectively, with average N50 lengths of 1419 bp. In addition, after G. affinis was exposed for 20 days, 169 genes were upregulated, and 176 genes were downregulated in liver. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were applied to all the genes to determine the genes' biological functions and processes. The results clearly showed that the differentially expressed genes were mainly related to immune pathways and metabolic correlation pathways. Interestingly, almost all the pathways related with the immunity were upregulated, while the metabolism pathways were downregulated. Lastly, quantitative real-time PCR (qRT-PCR) was performed to measure expressional levels of twelve genes confirmed through the DGE analysis. These results demonstrate that BaP damages immunity and enhances the consumption of all available energy storage to activate mechanisms of the detoxification in G. affinis. Up until now, the present study is the first time that a whole transcriptome sequencing analysis in the liver of G. affinis exposed to BaP has been reported.

Cloning and molecular characterization of estrogen-related receptor (ERR) and vitellogenin genes in the brackish water flea Diaphanosoma celebensis exposed to bisphenol A and its structural analogues

Although it has previously been shown that bisphenol (BP) analogues may interfere with the normal hormonal regulation by acting as endocrine disrupting chemicals (EDCs), little information is available on effects of BP analogues in invertebrates, particularly on cladocerans. In the present study, we identified estrogen-related receptors (EER), vitellogenin (VTG), and VTG receptor (VtgR) from the brackish water flea Diaphanosoma celebensis, and examined the effects of BPA and the substitutes, BPF and BPS, in different sublethal concentrations. Gene expression varied with time well matched with brooding, suggesting that DcEER, DcVTG, and DcVtgR play a role in reproduction in D. celebensis. qRT-PCR analysis showed that BPA and its substitutes differently modulated mRNA expressions of DcEER, DcVTG, and DcVtgR, indicating that these compounds adversely affect the normal reproduction-related pathway. This study facilitates better understanding of the molecular mode of action of BP analogues on the reproductive system of D. celebensis.

Acute toxicity of bisphenol A and its structural analogues and transcriptional modulation of the ecdysone-mediated pathway in the brackish water flea Diaphanosoma celebensis

Bisphenol A (BPA) is a representative endocrine disrupting chemical (EDC) that has estrogenic effects in aquatic animals. In recent years, due to the continuing usage of BPA, its analogues have been developed as alternative substances to replace its use. The molting process is a pivotal point in the development and reproduction of crustaceans. However, studies of the effects of EDCs on molting in crustaceans at the molecular level are scarce. In the present study, we examined the acute toxicity of BPA and its analogues bisphenol F (BPF) and S (BPS) to the brackish water flea Diaphanosoma celebensis. We further identified four ecdysteroid pathway - related genes (cyp314a1, EcRA, EcRB, and USP) in D. celebensis, and investigated the transcriptional modulation of these genes during molting and after exposure to BPA and its analogues for 48 h. Sequencing and phylogenetic analyses revealed that these four genes are highly conserved among arthropods and may be involved in development and reproduction in the adult stage. The mRNA expression patterns of cyp314a1, EcRA and USP were matched with the molting cycle, suggesting that these genes play a role in the molting process in the adult stage in cladocerans. Following relative real-time polymerase chain reaction (RT-PCR) analyses, BPA and its analogues were found to modulate the expression of each of these four genes differently, indicating that these compounds can disrupt the normal endocrine system function of D. celebensis. This study improves our understanding of the molecular mode of action of BPA and its analogues in D. celebensis.

Two novel CYP3A isoforms in marine mussel Mytilus coruscus: Identification and response to cadmium and benzo[a]pyrene

CYP3A enzymes play a crucial role in metabolic clearance of a variety of xenobiotics. However, their genetic information and function remain unclear in molluscs. In the present study, two novel CYP3A genes i.e. McCYP3A-1 and McCYP3A-2 were identified and characterized from the thick shell mussel Mytilus coruscus, and their tissue distribution as well as the response to cadmium (Cd) and benzo[a]pyrene (B[α]P) exposure were addressed using real time quantitative RT-PCR (qRT-PCR) and erythromycin N-demethylase (ERND) assay. McCYP3A-1 and McCYP3A-2 possess typically domains of CYP family such as helix-C, helix-I, helix-K, PERF and the heme binding domain as well as the characteristic domains of CYP3s including six SRS motifs. McCYP3A-1 and McCYP3A-2 transcripts were constitutively expressed in all examined tissues with high expression level in digestive glands, hepatopancreas and gonads. Upon B[α]P exposure, McCYP3A-1 and McCYP3A-2 mRNA expression in digestive glands showed a pattern of up-regulation followed by down-regulation, while under Cd exposure, showed a time-dependent induction profile. In addition, ERND activity, generally used as an indicator of CYP3, increased in a time-dependent manner after exposure to Cd and B[α]P. These results collectively indicated that McCYP3A-1 and McCYP3A-2 are CYP3A family member and may play a potential role in metabolic clearance of xenobiotics. Meanwhile, the current results may provide some baseline data to support McCYP3A-1 and McCYP3A-2 as candidate biomarkers for monitoring of PAHs and heavy metal pollution.