Next Article in Journal
Heavy Metals Content in Selected Medicinal Plants Produced and Consumed in Serbia and Their Daily Intake in Herbal Infusions
Previous Article in Journal
Successful Use of Continuous Veno-Venous Haemodialysis in a Case of Potential Lethal Caffeine Intoxication
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxics
Volume 11
Issue 2
10.3390/toxics11020197
toxics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis of Diethylstilbestrol in Daphnia Magna: Energy Metabolism and Growth Inhibition

by
Qi Li
*,
Qian Zhao
,
Jiahua Guo
,
Xi Li
and
Jinxi Song
Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, College of Urban and Environmental Sciences, Northwest University, Xi’an 710127, China
*
Author to whom correspondence should be addressed.
Toxics 2023, 11(2), 197; https://doi.org/10.3390/toxics11020197
Submission received: 7 January 2023 / Revised: 1 February 2023 / Accepted: 9 February 2023 / Published: 20 February 2023
(This article belongs to the Section Emerging Contaminants)
Browse Figures
Versions Notes

Abstract

:
With the widespread use of diethylstilbestrol (DES), it has become a common contaminant in the aquatic environment. It is toxic to a wide range of aquatic organisms, disrupting the water flea growth and further interfering with several ecosystem services. Nevertheless, the molecular mechanism of DES in water fleas is still unexplicit. In this study, the 21-day chronic test showed that a negative effect of growth and reproduction can be observed with DES exposure. Subsequently applied transcriptomic analysis illustrated the molecular mechanism in mode freshwater invertebrate Daphnia magna (D. magna) exposed to 2, 200, and 1000 μg·L−1 of DES for 9 days. Meanwhile, exposure to DES at 200 and 1000 μg·L−1 significantly restrains the growth (body length) and reproduction (first spawning time) of D. magna. Identified differentially expressed genes (DEGs) are majorly enriched relative to energy metabolism, lipid metabolism, the digestive system, transport and catabolism pathways which were remarkably changed. These repressed and up-regulated pathways, in relation to energy synthesis and metabolism, may be the reasons for the reduced body length and delayed first spawning time. Taken together, this study revealed that DES is a threat to D. magna in the aquatic environment and clarifies the molecular mechanism of the toxicity.

Graphical Abstract

1. Introduction

Freshwater crustacean Daphnia magna (D. magna) is the primary consumer of algae, bacteria, and protozoans, and is the primary forage for fish that delivers important ecological services in an aquatic ecosystem [1,2]. Because D. magna is susceptible to various toxic substances, the adverse effects on D. magna population and health probably disturb the multiple ecosystem services it provides [3]. In view of this, D. magna is frequently applied to hazard assessment and the classification of emerging chemicals, owing to its widespread distribution, short life cycle, sensitivity to chemical substances, parthenogenesis, and genetic homology [4,5,6]. With the application of D. magna in various chemical treatments, excellent ecotoxicological data have therefore been collected. For instance, most ecotoxicity data include the median effective concentration (EC50) and physiological indicators based on growth and reproduction endpoints. Yet, the molecular mechanism of chemicals in D. magna remains unknown, and concerns have been raised about the underlying mechanisms of adverse effects. Recent advances in ecotoxicology have shown that acute toxicity tests of chemicals have been transformed into chronic tests at environmentally relevant concentrations in order to elucidate the causes of these changes at the level of gene expression patterns [7]. Accompanied with the utilization of the quantification of total genome gene expression analysis, known as transcriptome analysis, it is possible to identify biological pathways disrupted by toxic substances [8,9]. This technique has been used to study the transient transcriptome alteration in D. magna after acute exposure to 17β-estradio to evaluate the dysregulated pathways [10]. Nevertheless, investigations related to the mechanism of action of diethylstilbestrol (DES) at the genetic level are still lacking.
DES is a typical synthetic estrogen, which is prescribed for preventing miscarriages, regulating estrogen secretion and animal growth [11,12]. Although DES have been prohibited for utilization in aquaculture and animal husbandry in many countries [13], it remains a threat to human health due to its illegal use in livestock production which leads to growing annual emissions [14,15]. The massive use has resulted in the content of DES in aquatic products at a μg/kg level; this accumulation might be a health threat to people and organisms [14]. Moreover, concentrations of DES as high as 24.9–102 ng·L−1 and 7.2–16.9 μg·L−1 were detected in some Chinese rivers and fisheries, respectively; this phenomenon showed that DES cannot be completely removed by a treatment plant but exists in natural water bodies for a long time [16,17]. Hence, DES existing in an aquatic environment may disrupt the health of organism. In D. magna, 17β-estradio interferes with RNA transport and signal transduction by competitively binding receptors, and affects the pathways related to steroid biosynthesis, lysosome, the intestinal immune system, and metabolism of active compounds [10]. Since DES has a stronger estrogenic activity than 17β-estradio, it can be inferred that DES treatment may cause similar genetic changes [18,19]. Laboratory studies demonstrated that exposure to DES caused reproductive effects in daphnids, and altered steroid metabolism capacities [20]. Successive DES treatments of F0 and F1 significantly inhibited the reproductive ability of offspring and molting frequency [21], but there was no research studying the chemical potential toxicity mechanism. At present, the application of transcriptome analysis has been a common method for identifying the toxic mechanism of the contaminant at a genome expression level [22,23]. While the chronic toxicological effects of DES on D. magna have been partially elucidated, the molecular mechanism of its inhibition of growth and reproduction based on the transcriptomic method remains poor.
In this study, the adverse effects of DES on the growth and reproduction of D. magna were evaluated on physiological and molecular levels by analyzing the alterations in physiological parameters and gene expression patterns. The chronic treatments were exposed to the series concentrations (20, 200, 1000 μg·L−1) of DES for 21 days to obtain physiological changes (body length, first spawning time, number of neonates, and molting frequency) of D. magna. Moreover, the molecular mechanisms of toxic effects were revealed by the changes in the transcriptome. Our hypothesis was that the exposure to range concentrations of DES may reduce body length, delay the reproduction time, and eventually result in growth and reproduction toxicity. In addition to the physiological changes, signaling pathways involved in energy metabolism, lipid metabolism, the digestive system, transport, and catabolism would be remarkably hindered. The objectives of this study were (1) to determine whether DES affects the life history parameters of D. magna at exposure concentrations; (2) to confirm the hypothesis by preliminarily clarifying the molecular mechanism of DES underlying the growth and reproduction inhibition.

2. Materials and Methods

2.1. D. magna Culture Maintenance

D. magna with homologous genes were provided by the Guangdong Provincial Laboratory Animal Research Institute. The D. magna were cultured in 2 L beakers containing artificial M4 medium with pH 7.8 ± 0.02, and the medium was renewed three times a week. All testing organisms were cultured in an incubator with the following conditions: a constant temperature of 21 ± 1 ℃, photoperiod of 16 h lightness:8 h darkness, and a light intensity of 1000–1500 lux. The D. magna were fed with green alga chlorella at a level of 2 × 106 cells/mL (0.16 mg total organic carbon) daily per organism.

2.2. Chemicals

The DES (CAS No. 6998-97-1; purity ≥ 99%) and dimethylsulfoxide (DMSO; CAS No. 67-68-5; purity ≥ 98%) were purchased from Shanghai Yuanye Bio-Technology Co., Ltd. (Shanghai, China). A previous study showed that less than 0.1% (v/v) DMSO was not toxic to D. magna [24]. The DES stock solutions were made by dissolving in DMSO and were sonicated in a water bath for 30 min. All experiments have solvent control groups. Other chemicals used for preparing M4 medium were at least reagent grade.

2.3. Acute and Chornic Toxicity Test

Referring to the OECD 202 guideline for daphnid, acute immobilization test [25], acute toxicity tests was conducted to study the reference toxicant (potassium dichromate) and DES. Using the range-finding method to assess the EC50, the DES experimental concentration was determined to be 1, 1.5, 2, 2.4, 2.8, and 3 mg·L−1 (four replicates per concentration), and 10 mg·L−1 stock solution was prepared with M4 medium applied for an acute test. During the experiment, test neonates (<24 h) were not fed and per replicate consisted of five neonates. After exposure for 48 h, the immobilized neonates were counted to draw a concentration–response curve to then calculate the 48 h EC50.
Following the OECD 211 guideline [26], we performed the chronic toxicity test. Neonates were exposed to a solvent control group and three treatment groups with concentration of 20, 200, 1000 μg·L−1 for 21 days. Each neonate was cultured in a 50 mL glass beaker containing 20 mL test solution, and per group consisted of 20 replicates. The test solutions were renewed; meanwhile, the test organisms were fed every two days. During the 21 days of exposure, the time of first spawning, molting frequency, and number of neonates were recorded every day per organism. The newborn offspring were removed from the solutions after being counted daily. The body length, from the top of the head to the base of the caudal spine, of each test organism was measured by a stereomicroscope at 7, 14, and 21 days. In the experimental groups, all physiological alterations of organisms were recorded.
With reference to the above chronic test, a short-term experiment exposed for 9 days was conducted to analyze the alterations of the transcriptome. The experiment method was the same as the chronic test; the only difference was that each group consisted of 36 replicates to guarantee the biomass for sequencing. After test organisms completed their first spawning at day 9, daphnia in each treatment were randomly divided into three groups that each contained 12 individuals placed into 1.5 mL cryogenic tubes, snap frozen in liquid nitrogen. All samples were stored at −80 °C before total RNA isolation.

2.4. RNA Sequencing

Total RNA in daphnia samples were extracted and the concentration and purity of total RNA were measured by Nanodrop 2000 (Thermo Scientific, Wilmington, DE, USA). 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA) was used to measure the concentration and purity of total RNA, and the RIN value of each sample was above 7.5. Total RNA samples were submitted to next-generation sequencing analysis (NGS) services with the application of Illumina HiSeq X Ten (Nanjing Personal Gene Technology Co., Ltd., Nanjing, China). The mapped reads of each sample varied from 52,595,459-57,663,891, with the coverage being from 94.78% to 95.77%.

2.5. RNA-Seq Data Analyses

For raw reads in the Fast Q format, the 3′ end _band adapter sequence and the low quality value (QV < 20) were deleted utilizing Cutadapt. The trimmed reads were mapped to the reference genome sequence of D. magna (GCA_003990815.1) using the HISAT2 (http://ccb.jhu.edu./software/hisat2/index.shtml (accessed on 16 November 2020)) software [27,28]. Following this, the read count value for each gene was estimated using the HTSeq statistics. The fragments per kilobase per million fragments (FPKM), serving as the gene expression level of read counts, were normalized to compare the expression levels. In order to verify the consistency of the three replications between groups, Pearson’s correlation analysis was estimated to examine the correlation of gene expression levels between samples. On the basis of gene expression, the similarities between solvent control and treatments were indicated by the utilization of principal component analysis (PCA) in the “DESeq” package of R software. The expression levels altered in DES-treated groups compared to solvent controls named differentially exposed genes (DEGs) were identified using the DESeq package, satisfying the cutoff: |log2 Fold Change (FC)| > 1 and adj p value < 0.05. The visual volcano plots and heat maps were generated, respectively, using the “ggplots2” and “Pheatmap” software packages. The Venn diagrams of different treatment groups were plotted to show the number of DEGs. The biological functions of DEGs were manifested by proceeding Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genome (KEGG) pathways analysis.
With respect to GO enrichment, analyses were conducted with application of the “TopGO” package; the annotations of DEGs contained the molecular function (MF), biological process (BP), and cell component (CC). Significant enrichments of GO terms were deemed when the p value < 0.05. With regards to KEGG pathway analysis, the p < 0.05 was considered as a necessary condition to determine the enriched pathways.

2.6. Quantitative Real-Time PCR

In order to verify the gene expression profile obtained by RNA-seq, the expression levels of genes related to the digestive system and amino acid metabolism pathways were studied, and the following genes were measured for verification, including nieman C2 protein (npc2) and 4-aminobutyrate aminotransferase (abat), by real-time quantitative polymerase chain reaction (qRT-PCR). The expression levels of the above genes were subjected to normalization of a housekeeping gene, β-actin butler (β-actin). The special PCR process can be found in the Supplementary Data.

2.7. Statistical Analyses

In this paper, GraphPad Prism 8 software (San Diego, CA, USA) was applied to determine the statistical significance between treatments. The differences in physiological parameters between the DES treatment groups and solvent control groups were evaluated by one-way analysis of variance (ANOVA), followed by Dunnett’s post hoc test and p < 0.05 was regarded as statistically significant. To determine the correlation between NGS and qRT-PCR, Pearson correlation analysis was adopted, and the significant correlation was limited with p < 0.05.

3. Results

3.1. Effects of Acute and Chronic Toxicity Test of DES

In an acute exposure test, no immobile behavior was observed in solvent control groups, while the swimming behavior of daphnia treated with DES were inhibited after exposure for 48 h. The concentration–response curve for the immobilization of DES was present in Figure 1; the EC50 of 1.89 mg·L−1 with a 95% confidence interval range from 1.81 to 1.91 mg·L−1 was assessed for D. magna. Based on the above EC50, the range of concentrations for the chronic test and transcriptome analysis were determined. In the chronic toxicity test, no immobile organisms were observed in all treatments. Exposure to the low (2 μg/L), medium (200 μg/L), and high (1000 μg/L) concentrations of DES remarkably delayed the first spawning time relative to solvent control groups (Figure 2a). Moreover, the DES treatment at the highest level significantly reduced the body length of D. magna compared with solvent control groups at 7, and 21 days, meaning that the growth of the organisms was inhibited during development (Figure 2b). In the medium to highest level of DES treatments, a tendency to decrease the number of neonates was observed but with no significance (Figure 2c). Conversely, an exposure to low concentrations of DES significantly increased the total number of neonates (Figure 2c). In addition, no significance was observed in all treatments of molting frequency (Figure S1).

3.2. Transcriptome Analysis

3.2.1. Differentially Expressed Genes

A total of 13,137 transcripts were examined in D. magna in all treatments. With regard to gene expression patterns, a correlation coefficient no less than 0.93 was detected among the replications in each treatment, demonstrating a high co-correlation between samples (Figure 3a). In the PCA analysis (Figure S2), gene expression showed a resemblance to low treatments and solvent control treatments, while the medium and high treatment differed from the solvent control treatment, indicating that DES at medium to high levels may result in hazardous effects that are consistent with the results in heatmaps (Figure 3b).
With exposure to DES, at three treatment levels, about 41 (2 up-regulated and 38 down-regulated), 31 (6 up-regulated and 25 down-regulated), and 26 (10 up-regulated and 16 down-regulated) genes were identified as DEGs in daphnia, respectively (Figure 3c and Figure S3–S5). To confirm these results, the selective DEGs were conducted with qRT-PCR. The FPKM value of selective DEGs was in accordance with the mRNA expression levels measured by qRT-PCR, showing that the mRNA expression profiles in D. magna supplied by transcriptomic sequencing were quantified with confidence. For instance, the FPKM values of npc2 and abat were remarkably in correlation with the mRNA expression levels measured by qRT-PCR with p < 0.0001 and R2 = 0.82, and p = 0.0002 and R2 = 0.78, respectively (Figure S6, Table S1).

3.2.2. Analyses of Gene Ontology and Functional Pathway

GO analysis was conducted to study the main biological functions of DEGs in the CC, MF, and BP. The top 20 GO enrichment terms in the low, medium, and high DES treatment groups were summarized (Figures S7–S9). In short, the DEs detected at low concentrations were mainly related to lipid biosynthesis processes (e.g., GO: 0008610 lipid biosynthesis process, 0006631 fatty acid metabolism process). Some GO terms were associated with energy metabolism in the medium treatment groups (e.g., GO: 0006486 protein glycosylation, GO: 0016798 hydrolase activity). DEGs in the high treatment groups were primarily enriched in GO terms related to the biological processes of wax (GO: 0010025 wax biosynthesis, GO: 0010166 wax metabolism) and energy metabolism (GO: 0009058 biosynthesis, GO: 0071704 organic metabolism). Therefore, these biological processes played an important role in the inhibition of DES on the growth of D. magna. The results of GO analysis were proven by enriched pathways connected with cholesterol metabolism, lipid metabolism, cutin, suberine and wax biosynthesis, and lysosome; the details are shown in Table 1. Special pathways in conjunction with other foundational biological processes are shown in Table S2.

4. Discussion

The ecotoxicity of DES to D. magna raised attention to its residues in the environment. However, studies related to transcriptome analysis and the molecular mechanism of DES were still deficient. The present study showed that chronic exposure to DES adversely affected the growth and reproduction of D. magna and revealed the toxic mechanism on a genetic level. Since estrogen drugs were highly toxic to D. magna, they served as references to further clarify the molecular mechanism of DES. In this study, the hypothesis that DES treatments may interfere with the pathways associated with energy metabolism, lipid metabolism, digestive system, transport, and catabolism was confirmed by transcriptome analysis. Consistent with this hypothesis, the down-regulation of the genes involved in the digestive system was closely related to the delay of the first spawning time and the decrease in body length after short-term exposure to DES. Thus, the potential associations between enriched signaling pathways and phenotypic changes were discussed in the following sections.

4.1. Genes Related to Energy Metabolism

Nitrogen is an important component of amino acids, and plays an important role in the structure and metabolism of organisms [29]. In this study, nitrogen metabolism was significantly inhibited by DES at low concentrations, in which the expression level of the carbonic anhydrase (ca) gene was down-regulated. Carbonic anhydrase (ca) is an active metal enzyme ubiquitous in organisms; its active site contains zinc ions bound to hydroxide, which can capture and absorb carbon dioxide and then catalyze the hydration reaction of carbon dioxide to produce bicarbonate [30]. Bicarbonate and ammonia can produce a precursor of the arginine, which is named the ammonia formate acid ester; this complex reaction mainly includes three steps: phosphorylation of bicarbonate, formation of carbamate, and phosphorylation of carbamate [31]. It can be speculated that the inhibition of ca activity can lead to a decrease in the bicarbonate content and inhibition of arginine synthesis. In addition, it has been found that bicarbonate catalyzed by ca plays an important role in maintaining the balance of intracellular pH value [32]. Since arginine deficiency can reduce growth and protein deposition, the arginine supplement can promote the growth of fish [33]. In view of this, the phenomenon of delayed first reproduction time of D. magna at a low concentration may be related to the reduction in arginine content caused by ca down-regulation; the lack of vital proteins greatly weakens the process of energy metabolism.
Amino acids are the precursors of bioactive molecules such as neurotransmitters, second messengers, and cytokines, which control various cellular processes, and its metabolism disorders can lead to various pathological phenomena [34]. In this study, the up-regulated 4-aminobutyrate transaminase (abat) gene is involved in amino acid metabolism pathways under a high DES level. It has been documented that abat is mainly responsible for the conversion of the inhibitory neurotransmitter 4-aminobutyric acid into succinate hemaldehydes, and finally succinate [35]. Succinate is regarded as a danger signal to promote the inflammatory response in the immune system, with increased levels in chronic inflammatory and metabolic diseases [36]. Therefore, the increase in succinate level induced by up-regulated abat in vivo indicated that the D. magna might be involved in immune metabolism, which proved that DES had a toxic effect on the organism at a high concentration. In addition, as reported when the articular cartilage of mice was damaged, the expression of abat in articular chondrocytes was significantly up-regulated. The induction of abat expression not only leads to metabolic changes in chondrocytes, but also increases the content of succinic acid which can result in increased mitochondrial respiration and consumption of a lot of energy [37]. Studies have shown that the overexpression of abat can lead to the catabolism of chondrocytes, but inhibition of abat expression can maintain chondrocyte homeostasis. It can be inferred that the significant decrease in the body length of D. magna at the high concentration treatment may be attributed to the overexpression of abat, which reduces chondrocytes and might cause bone deformation.

4.2. Genes Related to Lipid Metabolism

Sterol is an important component of biofilms, which is ubiquitous in crustaceans and has various biological functions [38]. In the low concentration treatment in this study, the gene of methylsterol monooxygenase (erg25), involved in the steroid biosynthesis pathway, was significantly down-regulated, which was consistent with the phenotypic changes of D. magna with the delayed first reproduction time. In the process of sterol synthesis, the initial lanosterol undergoes three demethylation before the final product of lanosterol is formed. The first demethylation occurs directly on lanosterol, leading to the removal of the C-14 methyl group, and the remaining demethylation is catalyzed by erg25, which eventually transforms dimethyllanosterol into the precursor of ergosterol called yeast sterol [39]. The lack of erg25 can lead to the accumulation of the dimethyllanosterol intermediate and its metabolites in vivo [40], ultimately resulting in toxic effects on organisms. In addition, ergosterol plays an important role in maintaining intracellular sterol homeostasis, and its deficiency leads to reduced plasma membrane mobility and intracellular ATP levels [41]. In short, the toxicity of DES on D. magna at a low concentration may be related to the abnormal expression of the erg25 gene, due to the homeostasis disrupted by massive intermediates and metabolites accumulated in cells.
Wax esters are macromolecules composed of long-chain fatty alcohols esterified into fatty acids, which are low-density neutral lipids used to store energy and contribute to the buoyance of copepods [41]. In the medium and high concentration treatments, the expression of aliphaloyl-coA reductase (far) genes involved in the wax biosynthesis pathway was up-regulated. The production of wax esters consists of two catalytic steps: fatty acids are reduced to fatty alcohols by far at first, and then fatty alcohols are trans-esterified to fatty acids to generate wax esters [42]. With the combined actions of far and wax ester hydrolase, the contents of fatty alcohols and wax esters in organisms maintain the dynamic balance. The generation rate of wax esters is limited only by fatty alcohols, which has been confirmed [43]. Thus, it is speculated that the induction of far expression can promote the production of wax esters, but some studies have shown that wax esters in adult copepods can be transferred to eggs during reproduction to ensure the development of eggs [44]. Therefore, it was inferred that the wax esters in D. magna may be transferred to the eggs during the reproduction process to meet the energy demands of offspring development, resulting in the energy deficiency in parental organisms, which is further in line with the negative changes in body length and number of offspring in the first brood.

4.3. Genes Related to Digestive System

The biliary acid is formed from cholesterol in the liver which plays an important role in the digestion and absorption of dietary lipids [45]. In the medium and high concentration treatments, the gene expressions of Nieman C2 protein (npc2) and cholesterol ester hydrolase (lipa) involved in the cholesterol metabolism pathway were significantly changed. The major function of lipa is to hydrolyze cholesterol esters into cholesterol, and its induction contributes to a large accumulation of free cholesterol in the late lysosome and inhibits the lysosome [46]. Npc2 is a key transporter driven by electrochemical potential, which regulates the transfer of cholesterol from lysosomes to the endoplasmic reticulum at the late stage. The efficiency of transporting cholesterol to the endoplasmic reticulum by vesicle and non-vesicle pathways was reduced by the down-regulation of npc2 [47]. In addition, the distribution and content of cholesterol in cells being in a highly dynamic balance is an important guarantee to maintain the integrity of cell function [48]. Thus, DES exposure may disrupt the metabolic process of cholesterol in cells and even interfere with the cell function. Bile acid is an important hydroxylated steroid synthesized from cholesterol in the liver and used to facilitate the absorption of nutrients, such as lipids and fat-soluble vitamins in the gut [49]. It has been proven that normal cholesterol metabolism is necessary to maintain a balance of lipids and bile acids in liver cells [44]. Hence, it was speculated that the passive effects of growth and reproduction on D. magna might be closely related to the inhibition of npc2. The reduction of cholesterol in the endoplasmic reticulum would significantly affect the secretion of bile acid, leading to a decreased ability of D. magna to digest algae and absorb the fat soluble cellulose.
Protein is eventually digested into amino acids, which are essential for normal growth, development, repair, and energy absorption of organisms [50]. In this study, gene expression of trypsin (prss1_2_3) and collagen type I (col1a) involved in protein digestion and the absorption pathway were both up-regulated in the high concentration treatments. Pancreatic digestive enzymes play a central role in digestion and usually are transported to the small intestine to hydrolyze complex nutrients, such as the endopeptidase of prss1_2_3 which is capable of splitting polypeptides into oligopeptides and amino acids [50]. The up-regulation of the prss1_2_3 gene may have aggravated the energy deficiency due to increased energy consumption. Col1a is an extracellular matrix protein as well an important part of the tumor mesenchymal environment, which plays a vital role in disrupting cell adhesion and activating β1-integrin to promote the epithelial–mesenchymal transformation in pancreatic cancer [51]. The expression of Col1a in liver cancer cells was significantly up-regulated, meaning that its overexpression could have affected the integrity of tissues and could lead to the diffusion and metastasis of tumor cells [52]. In the present study, the decreased body length of D. magna observed with a high level of DES treatment may be connected with abnormal liver cells; diseased liver cells may lead to the intake of green algae still not being digested, which might aggravate the problem of energy deficiency.

4.4. Genes Related to Transport and Catabolism

In the present study, the expression of lysosomal acid lipase (lipa) and Niemannc1 protein (npc) genes involved in the transport and catabolism pathway were significantly changed at medium and high concentrations. Lipa is an important regulatory factor of lysosomal lipolysis; meanwhile, the hydrolysis of cholesterol esters and triglycerides in lysosomes under acidic conditions is closely related to the production of lipid mediators, macrophage M2 activation, and inflammation [53]. As it was proven that the dysfunction of the lysosome is connected with the abnormal expression of lysosomal autophagy genes, the lysosomal somatic autophagy could be effectively blocked by down-regulated lipa [54]. Thus, the overexpression of lipa may be related to the toxic effect of the lysosome in the D. magna response to DES, exhibiting inflammatory reactions and consuming massive energy, ultimately leading to impaired energy transport. Moreover, npc acts as a transport protein on the lysosome membrane, and down-regulation of the gene can cause a gradual accumulation of cholesterol triglycerides, non-metabolic substrates such as proteins, and can result in lysosomal storage disorders; the continuous and vast accumulation of undegraded substrates in the lysosome will eventually lead to cell dysfunction and even death [55]. Therefore, lysosomes may be hydrolyzed and apoptotic in response to DES stress. Sequentially, the normal energy transport was restrained, and we observed the phenomenon that the growth and reproduction of D. magna were inhibited. Hence, DES exposure may inhibit the transport and catabolism of energy, which can account for the adverse physiological effect.

5. Conclusions

Overall, chronic exposure to DES remarkably affects the growth and reproduction of D. magna at low, medium, and high treatment groups; whereas, no mortality was observed in all treatment groups. DES treatment resulted in a difference in the D. magna transcriptome, demonstrating that the down-regulation of pathways in relation to energy metabolism, lipid metabolism, digestive system, transport, and catabolism may limit food digestion, energy absorption, and metabolic processes, ultimately leading to reduced body length and delayed reproduction of D. magna. In summary, DES treatments can observably affect the physiological phenotypes and molecular pathways of D. magna. Applying a transcriptomic approach, the differentially expressed genes and disordered molecular pathways were identified. They are recommended to serve as biomarkers to fully elucidate their toxic mechanism to aquatic organisms and provide a warning of health risk.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxics11020197/s1, Figure S1: Effects of DES exposure to D. magna in molting frequency; Figure S2: Principal component analysis (PCA) of DEGs of D. magna exposed to DES for 9d in transcriptomic analysis; Figure S3: Gene expression profiles of D. magna after 9 d exposure to DES at 2 μg L−1 volcano plot compared with solvent control groups; Figure S4: Gene expression profiles of D. magna after 9 d exposure to DES at 200 μg L−1 volcano plot compared with solvent control groups; Figure S5: Gene expression profiles of D. magna after 9 d exposure to DES at 1000 μg L−1 volcano plot compared with solvent control groups; Figure S6: qRT-PCR and RNA-seq correlation analysis results diagram of selected genes. NPC2: Niemann-Pick C2 protein; abat: 4-aminobutyrate aminotransferase; Figure S7: After a 9-d of exposure to different concentrations of DES, the top 20 enrichment terms for each category were significantly higher in the 2 μg L−1 treatment group than in the solvent control; Figure S8: After a 9-d of exposure to different concentrations of DES, the top 20 enrichment terms for each category were significantly higher in the 200 μg L−1 treatment group than in the solvent control; Figure S9: After a 9-d of exposure to different concentrations of DES, the top 20 enrichment terms for each category were significantly higher in the 200 μg L−1 treatment group than in the solvent control; Table S1: Primers utilised for quantitative polymerase chain reaction (qPCR) validation of gene expression; Table S2: List of fundamental pathways affected by a 9-day exposure to DES.

Author Contributions

Conceptualization, Q.L. and J.G. and J.S.; methodology, J.G. and Q.Z. and X.L.; software, X.L. and Q.Z.; validation, Q.L. and J.S.; formal analysis, J.G.; investigation, Q.Z. and X.L.; resources, Q.L. and J.S.; data curation, X.L.; writing—original draft preparation, Q.Z.; writing—review and editing, Q.L. and J.G. and J.S.; visualization, J.G.; supervision, Q.L.; project administration, Q.L.; funding acquisition, Q.L.. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China [Grant No. 42277414];The Second Tibetan Plateau Scientific Expedition and Research Program [Grant No. 2019QZKK1003]; Key Research and Development Plan of Shaanxi Province [Grant No. 2020SF-400]; Scientific Research Foundation of the Water Conservancy Department of Shaanxi Province [Grant No. 2020slkj-13]. And The APC was funded by the National Natural Science Foundation of China.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

We declare no competing interests.

References

  1. Carpenter, S.R.; Kitchell, J.F.; Hodgson, J.R.; Cochran, P.A.; Elser, J.J.; Elser, M.M.; Lodge, D.M.; Kretchmer, D.; He, X.; von Ende, C.N. Regulation of Lake Primary Productivity by Food Web Structure. Ecology 1987, 68, 1863–1876. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Tessier, A.J.; Leibold, M.A.; Tsao, J. A Fundamental Trade-Off in Resource Exploitation by Daphnia and Consequences to Plankton Communities. Ecology 2000, 81, 826. [Google Scholar] [CrossRef]
  3. Akbar, S.; Gu, L.; Sun, Y.; Zhang, L.; Lyu, K.; Huang, Y.; Yang, Z. Understanding host-microbiome-environment interactions: Insights from Daphnia as a model organism. Sci. Total. Environ. 2022, 808, 152093. [Google Scholar] [CrossRef] [PubMed]
  4. Toyota, K.; McNabb, N.A.; Spyropoulos, D.D.; Iguchi, T.; Kohno, S. Toxic effects of chemical dispersant Corexit 9500 on water flea Daphnia magna. J. Appl. Toxicol. 2017, 37, 201–206. [Google Scholar] [CrossRef]
  5. Bownik, A.; Kowalczyk, M.; Bańczerowski, J. Lambda-cyhalothrin affects swimming activity and physiological responses of Daphnia magna. Chemosphere 2019, 216, 805–811. [Google Scholar] [CrossRef]
  6. Gowler, C.D.; Rogalski, M.A.; Shaw, C.L.; Hunsberger, K.K.; Duffy, M.A. Density, parasitism, and sexual reproduction are strongly correlated in lake Daphnia populations. Ecol. Evol. 2021, 11, 10446–10456. [Google Scholar] [CrossRef]
  7. Brockmeier, E.K.; Hodges, G.; Hutchinson, T.H.; Butler, E.; Hecker, M.; Tollefsen, K.E.; Garcia-Reyero, N.; Kille, P.; Becker, D.; Chipman, K.; et al. The Role of Omics in the Application of Adverse Outcome Pathways for Chemical Risk Assessment. Toxicol. Sci. 2017, 158, 252–262. [Google Scholar] [CrossRef] [Green Version]
  8. Guo, J.; Ito, S.; Nguyen, H.T.; Yamamoto, K.; Tanoue, R.; Kunisue, T.; Iwata, H. Effects of prenatal exposure to triclosan on the liver transcriptome in chicken embryos. Toxicol. Appl. Pharmacol. 2018, 347, 23–32. [Google Scholar] [CrossRef]
  9. Fuertes, I.; Jordão, R.; Piña, B.; Barata, C. Time-dependent transcriptomic responses of Daphnia magna exposed to metabolic disruptors that enhanced storage lipid accumulation. Environ. Pollut. 2019, 249, 99–108. [Google Scholar] [CrossRef]
  10. Zheng, Y.; Yuan, J.; Gu, Z.; Yang, G.; Li, T.; Chen, J. Transcriptome alterations in female Daphnia (Daphnia magna) exposed to 17β-estradiol. Environ. Pollut. 2020, 261, 114208. [Google Scholar] [CrossRef]
  11. Hunt, P.A.; Sathyanarayana, S.; Fowler, P.A.; Trasande, L. Female Reproductive Disorders, Diseases, and Costs of Exposure to Endocrine Disrupting Chemicals in the European Union. J. Clin. Endocrinol. Metab. 2016, 101, 1562–1570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, X.; Ding, H.; Yu, X.; Shi, X.; Sun, A.; Li, D.; Zhao, J. Characterization and application of molecularly imprinted polymer-coated quantum dots for sensitive fluorescent determination of diethylstilbestrol in water samples. Talanta 2019, 197, 98–104. [Google Scholar] [CrossRef] [PubMed]
  13. Reed, C.E.; Fenton, S.E. Exposure to diethylstilbestrol during sensitive life stages: A legacy of heritable health effects. Birth Defects Res. Part C Embryo Today Rev. 2013, 99, 134–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Li, F.; Sun, J.; Jia, H.; Zhang, Y.; Pang, L.; He, L.; Chai, T. A disposable paper-based sample clean-up slides for the sensitive determination of trace diethylstilbestrol residues in aquatic products. Microchem. J. 2019, 151, 104243. [Google Scholar] [CrossRef]
  15. Król, J.; Pobłocki, W.; Bockenheimer, T.; Hliwa, P. Effect of diethylstilbestrol (DES) and 17 β-estradiol (E2) on growth, survival and histological structure of the internal organs in juvenile European catfish Silurus glanis (L.). Aquac. Int. 2014, 22, 53–62. [Google Scholar] [CrossRef] [Green Version]
  16. Chen, H.-C.; Kuo, H.-W.; Ding, W.-H. Determination of estrogenic compounds in wastewater using liquid chromatography–tandem mass spectrometry with electrospray and atmospheric pressure photoionization following desalting extraction. Chemosphere 2009, 74, 508–514. [Google Scholar] [CrossRef]
  17. Zhang, X.; Zhang, D.; Zhang, H.; Luo, Z.; Yan, C. Occurrence, distribution, and seasonal variation of estrogenic compounds and antibiotic residues in Jiulongjiang River, South China. Environ. Sci. Pollut. Res. 2012, 19, 1392–1404. [Google Scholar] [CrossRef]
  18. Nikov, G.N.; Eshete, M.; Rajnarayanan, R.V.; Alworth, W.L. Interactions of synthetic estrogens with human estrogen receptors. J. Endocrinol. 2001, 170, 137–145. [Google Scholar] [CrossRef] [Green Version]
  19. Shang, G.; Xue, J.; Li, M.; Hu, H.-Y.; Lu, Y. Estrogen receptor affinity chromatography: A new method for characterization of novel estrogenic disinfection by-products. Chemosphere 2014, 104, 251–257. [Google Scholar] [CrossRef]
  20. Baldwin, W.S.; Milam, D.L.; Leblanc, G.A. Physiological and biochemical perturbations inDaphnia magnafollowing exposure to the model environmental estrogen diethylstilbestrol. Environ. Toxicol. Chem. 1995, 14, 945–952. [Google Scholar] [CrossRef]
  21. Brennan, S.J.; Brougham, C.A.; Roche, J.J.; Fogarty, A.M. Multi-generational effects of four selected environmental oestrogens on Daphnia magna. Chemosphere 2006, 64, 49–55. [Google Scholar] [CrossRef] [PubMed]
  22. Watanabe, H.; Kobayashi, K.; Kato, Y.; Oda, S.; Abe, R.; Tatarazako, N.; Iguchi, T. Transcriptome profiling in crustaceans as a tool for ecotoxicogenomics. Cell Biol. Toxicol. 2008, 24, 641–647. [Google Scholar] [CrossRef] [PubMed]
  23. David, R.M.; Dakic, V.; Williams, T.D.; Winter, M.J.; Chipman, J.K. Transcriptional responses in neonate and adult Daphnia magna in relation to relative susceptibility to genotoxicants. Aquat. Toxicol. 2011, 104, 192–204. [Google Scholar] [CrossRef] [PubMed]
  24. Palma, P.; Palma, V.; Fernandes, R.; Bohn, A.; Soares, A.; Barbosa, I. Embryo-toxic effects of environmental concentrations of chlorpyrifos on the crustacean Daphnia magna. Ecotoxicol. Environ. Saf. 2009, 72, 1714–1718. [Google Scholar] [CrossRef]
  25. OECD. Test No. 202: Daphnia sp. Acute Immobilisation Test, OECD Guidelines for the Testing of Chemicals, Section 2; OECD Publishing: Paris, France, 2004. [Google Scholar] [CrossRef] [Green Version]
  26. OECD. Test No. 211: Daphnia magna Reproduction Test, OECD Guidelines for the Testing of Chemicals, Section 2; OECD Publishing: Paris, France, 2012. [Google Scholar] [CrossRef] [Green Version]
  27. Federico, A.; Serra, A.; Ha, M.K.; Kohonen, P.; Choi, J.-S.; Liampa, I.; Nymark, P.; Sanabria, N.; Cattelani, L.; Fratello, M.; et al. Transcriptomics in Toxicogenomics, Part II: Preprocessing and Differential Expression Analysis for High Quality Data. Nanomaterials 2020, 10, 903. [Google Scholar] [CrossRef]
  28. Kowal, K.; Tkaczyk, A.; Pierzchała, M.; Bownik, A.; Ślaska, B. Identification of Mitochondrial DNA (NUMTs) in the Nuclear Genome of Daphnia magna. Int. J. Mol. Sci. 2020, 21, 8725. [Google Scholar] [CrossRef]
  29. Supuran, C.T.; Simone, G.D. Carbonic Anhydrases as Biocatalysts: From Theory to Medical and Industrial Applications; Elsevier: Amsterdam, The Netherlands, 2015; ISBN 9780444632586. [Google Scholar] [CrossRef]
  30. Hu, L.; Diez-Fernandez, C.; Ruefenacht, V.; Hismi, B.Ö.; Ünal, Ö.; Soyucen, E.; Çoker, M.; Bayraktar, B.T.; Gunduz, M.; Kiykim, E.; et al. Recurrence of carbamoyl phosphate synthetase 1 (CPS1) deficiency in Turkish patients: Characterization of a founder mutation by use of recombinant CPS1 from insect cells expression. Mol. Genet. Metab. 2014, 113, 267–273. [Google Scholar] [CrossRef]
  31. Shi, J.; Jiang, Y.; Jiang, Z.; Wang, X.; Wang, X.; Zhang, S.; Han, P.; Yang, C. Enzymatic conversion of carbon dioxide. Chem. Soc. Rev. 2015, 44, 5981–6000. [Google Scholar] [CrossRef]
  32. Tu, Y.; Xie, S.; Han, D.; Yang, Y.; Jin, J.; Zhu, X. Dietary arginine requirement for gibel carp (Carassis auratus gibelio var. CAS III) reduces with fish size from 50g to 150g associated with modulation of genes involved in TOR signaling pathway. Aquaculture 2015, 449, 37–47. [Google Scholar] [CrossRef]
  33. Suzuki, A.; Iwata, J. Amino acid metabolism and autophagy in skeletal development and homeostasis. Bone 2021, 146, 115881. [Google Scholar] [CrossRef]
  34. Kim, Y.T.; Song, Y.H.; Churchich, J.E. Recombinant brain 4-aminobutyrate aminotransferases overexpression, purification, and identification of Lys-330 at the active site. Biochim. et Biophys. Acta (BBA)-Protein Struct. Mol. Enzym. 1997, 1337, 248–256. [Google Scholar] [CrossRef]
  35. Fernández-Veledo, S.; Ceperuelo-Mallafré, V.; Vendrell, J. Rethinking succinate: An unexpected hormone-like metabolite in energy homeostasis. Trends Endocrinol. Metab. 2021, 32, 680–692. [Google Scholar] [CrossRef]
  36. Basoli, V.; Li, Z.; Traweger, A.; Sanchez-Antequera, Y.; Plank, C.; Rip, J.; Alini, M.; Grad, S. Effect of nanoparticle based mrna delivery on modulation of inflammation in an osteochondral inflammation model. Osteoarthr. Cartil. 2021, 29, S13. [Google Scholar] [CrossRef]
  37. Ballantine, J.A.; Roberts, J.C.; Morris, R.J. The sterols of crustaceans: Decapods (Sub-order macrura). Comp. Biochem. Physiol. Part B Comp. Biochem. 1980, 67, 75–79. [Google Scholar] [CrossRef]
  38. Bard, M.; A Bruner, D.; A Pierson, C.; Lees, N.D.; Biermann, B.; Frye, L.; Koegel, C.; Barbuch, R. Cloning and characterization of ERG25, the Saccharomyces cerevisiae gene encoding C-4 sterol methyl oxidase. Proc. Natl. Acad. Sci. USA 1996, 93, 186–190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Ward, D.M.; Chen, O.S.; Li, L.; Kaplan, J.; Bhuiyan, S.A.; Natarajan, S.K.; Bard, M.; Cox, J.E. Altered sterol metabolism in budding yeast affects mitochondrial iron–sulfur (Fe-S) cluster synthesis. J. Biol. Chem. 2018, 293, 10782–10795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Suchodolski, J.; Muraszko, J.; Bernat, P.; Krasowska, A. A Crucial Role for Ergosterol in Plasma Membrane Composition, Localisation, and Activity of Cdr1p and H+-ATPase in Candida albicans. Microorganisms 2019, 7, 378. [Google Scholar] [CrossRef] [Green Version]
  41. Hofvander, P.; Doan, T.T.; Hamberg, M. A prokaryotic acyl-CoA reductase performing reduction of fatty acyl-CoA to fatty alcohol. FEBS Lett. 2011, 585, 3538–3543. [Google Scholar] [CrossRef] [Green Version]
  42. Sargent, J.; Gatten, R.; McIntosh, R. Biosynthesis of wax esters in cell-free preparations of euchaeta norvegica. Comp. Biochem. Physiol. Part B Comp. Biochem. 1974, 47, 217–227. [Google Scholar] [CrossRef]
  43. Bauermeister, A.; Sargent, J. Wax esters: Major metabolites in the marine environment. Trends Biochem. Sci. 1979, 4, 209–211. [Google Scholar] [CrossRef]
  44. Norlin, M.; Wikvall, K. Enzymes in the conversion of cholesterol into bile acids. Curr. Mol. Med. 2007, 7, 199–218. [Google Scholar] [CrossRef]
  45. Dubland, J.A.; Francis, G.A. Lysosomal acid lipase: At the crossroads of normal and atherogenic cholesterol metabolism. Front. Cell Dev. Biol. 2015, 3, 3. [Google Scholar] [CrossRef] [Green Version]
  46. Subramanian, K.; Balch, W.E. NPC1/NPC2 function as a tag team duo to mobilize cholesterol. Proc. Natl. Acad. Sci. USA 2008, 105, 15223–15224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Maxfield, F.R.; Menon, A.K. Intracellular sterol transport and distribution. Curr. Opin. Cell Biol. 2006, 18, 379–385. [Google Scholar] [CrossRef] [PubMed]
  48. Xie, A.-J.; Mai, C.-T.; Zhu, Y.-Z.; Liu, X.-C.; Xie, Y. Bile acids as regulatory molecules and potential targets in metabolic diseases. Life Sci. 2021, 287, 120152. [Google Scholar] [CrossRef]
  49. Whitcomb, D.C.; Lowe, M.E. Human Pancreatic Digestive Enzymes. Dig. Dis. Sci. 2007, 52, 1–17. [Google Scholar] [CrossRef]
  50. Koenig, A.; Mueller, C.; Hasel, C.; Adler, G.; Menke, A. Collagen Type I Induces Disruption of E-Cadherin–Mediated Cell-Cell Contacts and Promotes Proliferation of Pancreatic Carcinoma Cells. Cancer Res. 2006, 66, 4662–4671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Liu, G.; Liu, R.; Shan, Y.; Sun, C. Marine bacterial exopolysaccharide EPS11 inhibits migration and invasion of liver cancer cells by directly targeting collagen I. J. Biol. Chem. 2021, 297, 101133. [Google Scholar] [CrossRef]
  52. Emanuel, R.; Sergin, I.; Bhattacharya, S.; Turner, J.N.; Epelman, S.; Settembre, C.; Diwan, A.; Ballabio, A.; Razani, B. Induction of Lysosomal Biogenesis in Atherosclerotic Macrophages Can Rescue Lipid-Induced Lysosomal Dysfunction and Downstream Sequelae. Arter. Thromb. Vasc. Biol. 2014, 34, 1942–1952. [Google Scholar] [CrossRef] [Green Version]
  53. Hu, Z.-Q.; Li, Q.; Hu, Z.-H.; Liu, H.-C.; Rao, C.-L.; Zhang, M.-J.; Xia, Y.-P.; Deng, L.; Mao, X.-H.; Fang, Y. MicroRNA-146a inhibits autophagy to maintain the intracellular survival of Burkholderia pseudomallei by targeting LIPA. Microb. Pathog. 2021, 158, 104969. [Google Scholar] [CrossRef] [PubMed]
  54. Klein, A.; Futerman, A.H. Lysosomal storage disorders: Old diseases, present and future challenges. Pediatr. Endocrinol. Rev. 2013, 11 (Suppl. S1), 59–63. [Google Scholar]
  55. Lei, K.; Lin, C.-Y.; Zhu, Y.; Chen, W.; Pan, H.-Y.; Sun, Z.; Sweetman, A.; Zhang, Q.; He, M.-C. Estrogens in municipal wastewater and receiving waters in the Beijing-Tianjin-Hebei region, China: Occurrence and risk assessment of mixtures. J. Hazard. Mater. 2020, 389, 121891. [Google Scholar] [CrossRef] [PubMed]
Toxics 11 00197 g001 550
Figure 1. The immobilization rate of DES during 48 h test. Data are mean ± SD (n = 8).
Figure 1. The immobilization rate of DES during 48 h test. Data are mean ± SD (n = 8).
Toxics 11 00197 g001
Toxics 11 00197 g002 550
Figure 2. Physiological parameters of D. magna exposed to DES during chronic test for 21 d. (a) The first spawning time of D. magna; (b) the body length of D. magna at 7, 14, 21 d; (c) the total number of neonates of D. magna. “*” suggesting the parameters are statistically different between treatment groups and solvent control groups (p < 0.05 and n = 20).
Figure 2. Physiological parameters of D. magna exposed to DES during chronic test for 21 d. (a) The first spawning time of D. magna; (b) the body length of D. magna at 7, 14, 21 d; (c) the total number of neonates of D. magna. “*” suggesting the parameters are statistically different between treatment groups and solvent control groups (p < 0.05 and n = 20).
Toxics 11 00197 g002
Toxics 11 00197 g003a 550Toxics 11 00197 g003b 550
Figure 3. Transcriptomic profiles of D. magna after 9 d exposure to diethylstilbestrol (DES). (a) Correlation analysis of patterns of gene expression in solvent control groups and DES treatment groups; (b) a heatmap of centered and scaled FPKM value of DEGs in solvent control groups and DES treatment groups; (c) Venn diagram of the number of DEGs in each DES treatment group. SC: solvent control; L: low; M: medium; H: high.
Figure 3. Transcriptomic profiles of D. magna after 9 d exposure to diethylstilbestrol (DES). (a) Correlation analysis of patterns of gene expression in solvent control groups and DES treatment groups; (b) a heatmap of centered and scaled FPKM value of DEGs in solvent control groups and DES treatment groups; (c) Venn diagram of the number of DEGs in each DES treatment group. SC: solvent control; L: low; M: medium; H: high.
Toxics 11 00197 g003aToxics 11 00197 g003b
Table
Table 1. KEGG pathways significantly enriched (p < 0.05) in D. magna exposed to DES.
Table 1. KEGG pathways significantly enriched (p < 0.05) in D. magna exposed to DES.
PathwaysCategoryUp-GenesDown-Genes
SC vs. L
Steroid biosynthesisLipid metabolism-meso1, erg25
Terpenoid backbone biosynthesisMetabolism of terpenoids and polyketides-pdss1
Nitrogen metabolismEnergy metabolism-ca
Cell adhesion molecules (CAMs)Signaling molecules and interaction-cntnap2
SC vs. M
Cholesterol metabolismDigestive system-npc2
LysosomeTransport and catabolismhgsnatnpc
Cutin, suberine, and wax biosynthesisLipid metabolismfar-
Glycosaminoglycan degradationGlycan biosynthesis and metabolismhgsnat-
SC vs. H
Cholesterol metabolismDigestive systemlalnpc1, npc2
LysosomeTransport and catabolismlipanpc1, npc2
Cutin, suberine, and wax biosynthesisLipid metabolismfar-
Alanine, aspartate, and glutamate metabolismAmino acid metabolismabat-
Protein digestion and absorptionDigestive systemprss1_2_3, col1a-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Q.; Zhao, Q.; Guo, J.; Li, X.; Song, J. Transcriptomic Analysis of Diethylstilbestrol in Daphnia Magna: Energy Metabolism and Growth Inhibition. Toxics 2023, 11, 197. https://doi.org/10.3390/toxics11020197

AMA Style

Li Q, Zhao Q, Guo J, Li X, Song J. Transcriptomic Analysis of Diethylstilbestrol in Daphnia Magna: Energy Metabolism and Growth Inhibition. Toxics. 2023; 11(2):197. https://doi.org/10.3390/toxics11020197

Chicago/Turabian Style

Li, Qi, Qian Zhao, Jiahua Guo, Xi Li, and Jinxi Song. 2023. "Transcriptomic Analysis of Diethylstilbestrol in Daphnia Magna: Energy Metabolism and Growth Inhibition" Toxics 11, no. 2: 197. https://doi.org/10.3390/toxics11020197

APA Style

Li, Q., Zhao, Q., Guo, J., Li, X., & Song, J. (2023). Transcriptomic Analysis of Diethylstilbestrol in Daphnia Magna: Energy Metabolism and Growth Inhibition. Toxics, 11(2), 197. https://doi.org/10.3390/toxics11020197

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop