Studying mixture effects on uptake and tissue distribution of PFAS in zebrafish (Danio rerio) using physiologically based kinetic (PBK) modelling

Per- and polyfluoroalkyl substances (PFAS) are ubiquitously distributed in the aquatic environment. They include persistent, mobile, bioaccumulative, and toxic chemicals and it is therefore critical to increase our understanding on their adsorption, distribution, metabolism, excretion (ADME). The current study focused on uptake of seven emerging PFAS in zebrafish (Danio rerio) and their potential maternal transfer. In addition, we aimed at increasing our understanding on mixture effects on ADME by developing a physiologically based kinetic (PBK) model capable of handling co-exposure scenarios of any number of chemicals. All studied chemicals were taken up in the fish to varying degrees, whereas only perfluorononanoate (PFNA) and perfluorooctanoate (PFOA) were quantified in all analysed tissues. Perfluorooctane sulfonamide (FOSA) was measured at concerningly high concentrations in the brain (Cmax over 15 μg/g) but also in the liver and ovaries. All studied PFAS were maternally transferred to the eggs, with FOSA and 6:2 perfluorooctane sulfonate (6,2 FTSA) showing significant (p < 0.02) signs of elimination from the embryos during the first 6 days of development, while perfluorobutane sulfonate (PFBS), PFNA, and perfluorohexane sulfonate (PFHxS) were not eliminated in embryos during this time-frame. The mixture PBK model resulted in >85 % of predictions within a 10-fold error and 60 % of predictions within a 3-fold error. At studied levels of PFAS exposure, competitive binding was not a critical factor for PFAS kinetics. Gill surface pH influenced uptake for some carboxylates but not the sulfonates. The developed PBK model provides an important tool in understanding kinetics under complex mixture scenarios and this use of New Approach Methodologies (NAMs) is critical in future risk assessment of chemicals and early warning systems.

Physiology-informed toxicokinetic model for the zebrafish embryo test developed for bisphenols

Zebrafish embryos (ZFE) is a widely used model organism, employed in various research fields including toxicology to assess e.g., developmental toxicity and endocrine disruption. Variation in effects between chemicals are difficult to compare using nominal dose as toxicokinetic properties may vary. Toxicokinetic (TK) modeling is a means to estimate internal exposure concentration or dose at target and to enable extrapolation between experimental conditions and species, thereby improving hazard assessment of potential pollutants. In this study we advance currently existing TK models for ZFE with physiological ZFE parameters and novel experimental bisphenol data, a class of chemicals with suspected endocrine activity. We developed a five-compartment model consisting of water, plastic, chorion, yolk sack and embryo in which surface area and volume changes as well as the processes of biotransformation and blood circulation influence mass fluxes. For model training and validation, we measured internal concentrations in ZFE exposed individually to BPA, bisphenol AF (BPAF) and Z (BPZ). Bayesian inference was applied for parameter calibration based on the training data set of BPZ. The calibrated TK model predicted internal ZFE concentrations of the majority of external test data within a 5-fold error and half of the data within a 2-fold error for bisphenols A, AF, F, and tetrabromo bisphenol A (TBBPA). We used the developed model to rank the hazard of seven bisphenols based on predicted internal concentrations and measured in vitro estrogenicity. This ranking indicated a higher hazard for BPAF, BPZ, bisphenol B and C (BPB, BPC) than for BPA.

Physiologically Based Toxicokinetic Modeling of Bisphenols in Zebrafish (Danio rerio) Accounting for Variations in Metabolic Rates, Brain Distribution, and Liver Accumulation

Bisphenol A (BPA) is an industrial chemical, which has raised human health and environmental concerns due to its endocrine-disrupting properties. BPA analogues are less well-studied despite their wide use in consumer products. These analogues have been detected in water and aquatic organisms around the world, with some analogues showing toxic effects in various species including fish. Here, we present novel organ-specific time-course distribution data of bisphenol Z (BPZ) in female zebrafish (Danio rerio), including concentrations in the ovaries, liver, and brain, a rarely sampled organ with high toxicological relevance. Furthermore, fish-specific in vitro biotransformation rates were determined for 11 selected bisphenols. A physiologically based toxicokinetic (PBTK) model was adapted for four of these bisphenols, which was able to predict levels in the gonads, liver, and brain as well as the whole body within a 2-5-fold error with respect to experimental data, covering several important target organs of toxicity. In particular, predicted liver concentrations improved compared to currently available PBTK models. Predicted data indicate that studied bisphenols mainly distribute to the carcass and gonads and less to the brain. Our model provides a tool to increase our understanding on the distribution and kinetics of a group of emerging pollutants.

Differential expression of enzymes in thymidylate biosynthesis in zebrafish at different developmental stages: implications for dtymk mutation-caused neurodegenerative disorders

Background

Deoxythymidine triphosphate (dTTP) is an essential building block of DNA, and defects in enzymes involved in dTTP synthesis cause neurodegenerative disorders. For instance, mutations in DTYMK, the gene coding for thymidylate kinase (TMPK), cause severe microcephaly in human. However, the mechanism behind this is not well-understood. Here we used the zebrafish model and studied (i) TMPK, an enzyme required for both the de novo and the salvage pathways of dTTP synthesis, and (ii) thymidine kinases (TK) of the salvage pathway in order to understand their role in neuropathology.

Results

Our findings reveal that maternal-stored dNTPs are only sufficient for 6 cell division cycles, and the levels of dNTPs are inversely correlated to cell cycle length during early embryogenesis. TMPK and TK activities are prominent in the cytosol of embryos, larvae and adult fish and brain contains the highest TMPK activity. During early development, TMPK activity increased gradually from 6 hpf and a profound increase was observed at 72 hpf, and TMPK activity reached its maximal level at 96 hpf, and remained at high level until 144 hpf. The expression of dtymk encoded Dtymk protein correlated to its mRNA expression and neuronal development but not to the TMPK activity detected. However, despite the high TMPK activity detected at later stages of development, the Dtymk protein was undetectable. Furthermore, the TMPK enzyme detected at later stages showed similar biochemical properties as the Dtymk enzyme but was not recognized by the Dtymk specific antibody.

Conclusions

Our results suggest that active dNTP synthesis in early embryogenesis is vital and that Dtymk is essential for neurodevelopment, which is supported by a recent study of dtymk knockout zebrafish with neurological disorder and lethal outcomes. Furthermore, there is a novel TMPK-like enzyme expressed at later stages of development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12868-022-00704-0.

Gene co-expression network analysis reveals mechanisms underlying ozone-induced carbamazepine toxicity in zebrafish (Danio rerio) embryos

Sewage effluent ozonation can reduce concentrations of chemical pollutants including pharmaceutical residues. However, the formation of potentially toxic ozonation byproducts (OBPs) is a matter of concern. This study sought to elucidate toxicity mechanisms of ozonated carbamazepine (CBZ), an anti-epileptic drug frequently detected in sewage effluents and surface water, in zebrafish embryos (Danio rerio). Embryos were exposed to ozonated and non-ozonated CBZ from 3 h post-fertilization (hpf) until 144 hpf. Embryotoxicity endpoints (proportion of dead and malformed embryos) were assessed at 24, 48, and 144 hpf. Heart rate was recorded at 48 hpf. Exposure to ozonated CBZ gave rise to cardiovascular-related malformations and reduced heart rate. Moreover, embryo-larvae exposed to ozonated CBZ displayed a lack of swim bladder inflation. Hence, the expression patterns of CBZ target genes involved in cardiovascular and embryonal development were investigated through a stepwise gene co-expression analysis approach. Two co-expression networks and their upstream transcription regulators were identified, offering mechanistic explanations for the observed toxicity phenotypes. The study presents a novel application of gene co-expression analysis elucidating potential toxicity mechanisms of an ozonated pharmaceutical with environmental relevance. The resulting data was used to establish a putative adverse outcome pathway (AOP).

Occurrence and removal of chemicals of emerging concern in wastewater treatment plants and their impact on receiving water systems

Wastewater treatment plants (WWTPs) are considered the main sources of chemicals of emerging concern (CECs) in aquatic environments, and can negatively impact aquatic ecosystems. In this study, WWTP influent, effluent, and sludge, and upstream and downstream waters from the WWTP recipient were investigated at 15 locations for a total of 164 CECs, including pharmaceuticals, personal care products, industrial chemicals, per- and polyfluoroalkyl substances (PFASs), and pesticides. In addition, zebrafish (Danio rerio) embryo toxicity tests (ZFET) were applied to WWTP influent and effluent, and upstream and downstream waters from WWTP recipients. A total of 119 CECs were detected in at least one sample, mean concentrations ranging from 0.11 ng/L (propylparaben) to 64,000 ng/L (caffeine), in wastewater samples and from 0.44 ng/L (ciprofloxacin) to 19,000 ng/L (metformin) in surface water samples. Large variations of CEC concentrations were found between the selected WWTPs, which can be explained by differences in CEC composition in influent water and WWTP treatment process. The sludge-water partitioning coefficient (K_d) of CECs showed a significant linear correlation to octanol/water partition coefficient (K_OW) (p < 0.001), and thus could be used for predicting their fate in the aqueous and solid phase. The ΣCEC concentrations in WWTPs declined by on average 60%, based on comparisons of WWTP influent and effluent concentrations. The high concentrations of CECs in WWTP effluent resulted in, on average, 50% higher concentrations of CECs in water downstream of WWTPs compared with upstream. Some WWTP samples showed toxicity in ZFET compared with the respective control group, but no individual CECs or groups of CECs could explain this toxicity. These results could provide a theoretical basis for optimization of existing treatment systems of different designs, and could significantly contribute to protecting recipient waters.

Behavioural effects and bioconcentration of per- and polyfluoroalkyl substances (PFASs) in zebrafish (Danio rerio) embryos

Per- and polyfluoroalkyl substances (PFASs) are contaminants of emerging concern (CECs) that cause concern regarding their environmental impact and risk to human health. In this study, zebrafish (Danio rerio) embryos were exposed to PFASs for six days, to investigate behaviour toxicity and bioconcentration factor (BCF). Nine individual PFASs (five C₄-C₈ perfluoroalkyl carboxylates (PFCAs) (PFPeA, PFHxA, PFHpA, PFOA, PFNA), three C₄, C₆ and C₈ perfluoroalkyl sulfonates (PFSAs) (PFBS, PFHxS, PFOS) and 6:2 fluorotelomersulfonate (6:2 FTSA)) and a mixture of these were investigated at seven concentrations ranging from environmentally relevant to acutely toxic levels. In exposed embryos, significant differences were found in total swimming distance (PFHpA, PFOA, PFNA, PFHxS, PFOS, 6:2 FTSA, PFAS mixture), burst activity (PFOA, PFNA, PFHxS, PFOS, PFAS mixture) and startle response (PFNA, PFHxS, PFOS, PFAS mixture). Toxicity was only observed at concentrations well above environmental levels. The toxicity of the PFAS mixture generally followed that of the individual substances, but the mixture reduced the potencies of individual PFASs. BCF was determined for all nine PFASs and ranged between 0.9 (PFPeA) and 2700 (PFOS). Long-chain PFASs (C₈) and PFASs with sulfonate as an active group showed the greatest toxic potential, while short-chain PFASs (C₆ and C₇) also caused significant behaviour alterations and accumulated in the embryos. To our knowledge, this is the first study to compare the behaviour toxicity of a PFAS mixture with that of the individual PFASs. Follow-up studies are needed to identify the mechanistic responses to PFAS mixtures.

Carbamazepine Ozonation Byproducts: Toxicity in Zebrafish (Danio rerio) Embryos and Chemical Stability

Carbamazepine (CBZ) is an anticonvulsant medication with highly persistent properties in the aquatic environment, where it has the potential to affect nontarget biota. Because CBZ and many other pharmaceuticals are not readily removed in conventional sewage treatment plants (STP), additional STP effluent treatment technologies are being evaluated and implemented. Whole effluent ozonation is a prospective method to remove pharmaceuticals such as CBZ, yet knowledge on the toxicity of CBZ ozonation byproducts (OBPs) is lacking. This study presents, for the first time, in vivo individual and mixture toxicity of four putative OBPs, that is, carbamazepine 10,11-epoxide, 10,11-Dihydrocarbamazepine, 1-(2-benzaldehyde)-4-hydro-(1H,3H)-quinazoline-2-one (BQM), and 1-(2-benzaldehyde)-(1H,3H)-quinazoline-2,4-dione (BQD) in developing zebrafish (Danio rerio) embryos. BQM and BQD were isolated from the ozonated solution as they were not commercially available. The study confirmed that the OBP mixture caused embryotoxic responses comparable to that of ozonated CBZ. Individual compound embryotoxicity assessment further revealed that BQM and BQD were the drivers of embryotoxicity. OBP chemical stability in ozonated CBZ water solution during 2 week dark storage at 22 °C was also assessed. The OBP concentrations remained over time, except for BQD which decreased by 94%. Meanwhile, ozonated CBZ persistently induced embryotoxicity over 2 week storage, potentially illustrating environmental concern.

Effects of antidepressants with different modes of action on early life stages of fish and amphibians

Drugs are excreted from the human body as both original substances and as metabolites and enter aquatic environment through waste water. The aim of this study was to widen the current knowledge considering the effects of waterborne antidepressants with different modes of action-amitriptyline, venlafaxine, sertraline-on embryos of non-target aquatic biota-fish (represented by Danio rerio) and amphibians (represented by Xenopus tropicalis). The tested concentrations were 0.3; 3; 30; 300 and 3000 μg/L in case of amitriptyline and venlafaxine and 0.1; 1; 10; 100 and 1000 μg/L for sertraline. Test on zebrafish embryos was carried out until 144 h post fertilization, while test on Xenopus embryos was terminated after 48 h. Lethal and sublethal effects as well as swimming alterations were observed at higher tested concentrations that are not present in the environment. In contrast, mRNA expression of genes related to heart, eye, brain and bone development (nkx2.5, otx 2, bmp4 and pax 6) seems to be impacted also at environmentally relevant concentrations. In a wider context, this study reveals several indications on the ability of antidepressants to affect non target animals occupying environments which may be contaminated by such compounds.

Embryotoxicity of ozonated diclofenac, carbamazepine, and oxazepam in zebrafish (Danio rerio)

Pharmaceutical residues are polluting the surface water environments worldwide. Sewage and wastewater treatment, therefore, needs to be improved in order to remove pharmaceutical residues from the effluent. One such treatment improvement is effluent ozonation. Even though ozonation has proven to be very efficient in reducing pharmaceutical parent compound concentrations in wastewater effluents, much remains unclear regarding potentially toxic ozonation by-product (OBP) formation. In this study, we sought to elucidate the aquatic toxicity of ozonated pharmaceuticals in zebrafish (Danio rerio) embryos in a static 144 h post fertilization (hpf) fish embryotoxicity (ZFET) assay. Three pharmaceuticals commonly detected in wastewater effluents, i.e. carbamazepine, diclofenac, and oxazepam, were selected for testing. Toxicity was assessed before and after 1 min ozonation (0.053 mg L^-1 peak O₃ concentration) and 10 min ozonation (0.147 mg L^-1 peak O₃ concentration). Chemical analysis showed that carbamazepine and diclofenac were largely removed by ozone (90 ± 11% and 97 ± 3.8%), whereas oxazepam was removed to a lesser extent (19 ± 5.7%). The ZFET assay revealed diverging toxicities. Diclofenac embryotoxicity decreased with increasing ozonation. Oxazepam did not cause embryotoxicity in the ZFET assay either pre- or post ozonation, but larvae swimming activity was affected at 144 hpf. Carbamazepine embryotoxicity, on the other hand, increased with increasing ozonation. Chemical analysis showed the formation of two OBPs (carbamazepine-10,11-epoxide and 10,11-dihydrocarbamazepine), possibly explaining the increased embryotoxicity. The results of this study highlight the importance of new chemical and toxicological knowledge regarding the formation of OBPs in post-ozonated effluents.

Swimming activity in zebrafish larvae exposed to veterinary antiparasitic pharmaceuticals

Veterinary antiparasitic pharmaceuticals have been detected in surface waters and several of these pharmaceuticals act on the nervous system on the target organisms implying that neurological effects also might be of concern in non-target animals such as fish. We tested if exposure to antiparasitic pharmaceuticals affect swimming activity in 6 days old zebrafish larvae. The results revealed that most pharmaceuticals did not cause any effects in swimming activity. However, larvae exposed to 0.58 mg/L doramectin displayed reduced swimming activity even though they were classified as normal, having no morphological abnormalities.

Effects of ozonated sewage effluent on reproduction and behavioral endpoints in zebrafish (Danio rerio)

Pharmaceutical residues and other micro-contaminants may enter aquatic environments through effluent from sewage treatment plants (STPs) and could cause adverse effects in wild fish. One strategy to alleviate this situation is to improve wastewater treatment by ozonation. To test the effectiveness of full-scale wastewater effluent ozonation at a Swedish municipal STP, the added removal efficiency was measured for 105 pharmaceuticals. In addition, gene expression, reproductive and behavioral endpoints were analyzed in zebrafish (Danio rerio) exposed on-site over 21 days to ozonated or non-ozonated effluents as well as to tap water. Ozone treatment (7 g O₃/m³) removed pharmaceuticals by an average efficiency of 77% in addition to the conventional treatment, leaving 11 screened pharmaceuticals above detection limits. Differences in biological responses of the exposure treatments were recorded in gene expression, reproduction and behavior. Hepatic vitellogenin gene expression was higher in male zebrafish exposed to the ozonated effluent compared to the non-ozonated effluent and tap water treatments. The reproductive success was higher in fish exposed to ozonated effluent compared to non-ozonated effluent and to tap water. The behavioral measurements showed that fish exposed to the ozonated STP effluent were less active in swimming the first minute after placed in a novel vessel. Ozonation is a capable method for removing pharmaceuticals in effluents. However, its implementation should be thoroughly evaluated for any potential biological impact. Future research is needed for uncovering the factors which produced the in vivo responses in fish.

Toxicity bioassays with concentrated cell culture media-a methodology to overcome the chemical loss by conventional preparation of water samples

The use of in vitro bioassays for studies of toxic activity in environmental water samples is a rapidly expanding field of research. Cell-based bioassays can assess the total toxicity exerted by a water sample, regardless whether the toxicity is caused by a known or unknown agent or by a complex mixture of different agents. When using bioassays for environmental water samples, it is often necessary to concentrate the water samples before applying the sample. Commonly, water samples are concentrated 10-50 times. However, there is always a risk of losing compounds in the sample in such sample preparation. We have developed an alternative experimental design by preparing a concentrated cell culture medium which was then diluted in the environmental water sample to compose the final cell culture media for the in vitro assays. Water samples from five Swedish waste water treatment plants were analyzed for oxidative stress response, estrogen receptor (ER), and aryl hydrocarbon receptor (AhR) activity using this experimental design. We were able to detect responses equivalent to 8.8-11.3 ng/L TCCD for AhR activity and 0.4-0.9 ng/L 17β-estradiol for ER activity. We were unable to detect oxidative stress response in any of the studied water samples. In conclusion, we have developed an experimental design allowing us to examine environmental water samples in toxicity in vitro assays at a concentration factor close to 1, without the risk of losing known or unknown compounds during an extraction procedure.

Sexual disruption in zebrafish (Danio rerio) exposed to mixtures of 17α-ethinylestradiol and 17β-trenbolone

Environmental estrogens and androgens can be present simultaneously in aquatic environments and thereby interact to disturb multiple physiological systems in organisms. Studies on interaction effects in fish of androgenic and estrogenic chemicals are limited. Therefore, the aim of the present study was to evaluate feminization and masculinization effects in zebrafish (Danio rerio) exposed to combinations of two synthetic steroid hormones detected in environmental waters: the androgen 17β-trenbolone (Tb) and the oestrogen 17α-ethinylestradiol (EE2). Juvenile zebrafish were exposed between days 20 and 60 post-hatch to different binary mixtures of Tb (1, 10, and 50 ng/L) and EE2 (2 and 5 ng/L). The endpoints studied were whole-body homogenate vitellogenin concentration at 40 days post-hatch, and sex ratio including gonad maturation at 60 days post-hatch. The feminizing potency of 5 ng/L of EE2, alone as well as in combination with Tb, was clear in the present study, with exposures resulting in almost all-female populations and females being sexually immature. Masculinization effects with male-biased sex ratios were observed when fish were exposed to 2 ng/L of EE2 in combination with Tb concentrations. Intersex fish were observed after exposure to mixtures of 2 ng/L EE2 with 50 ng/L Tb. Sexual maturity generally increased among males at increasing concentrations of Tb. The results of the present study show that exposure to environmentally relevant mixtures of an oestrogen and androgen affects the process of gonad differentiation in zebrafish and lead to sexual disruption.

Tissue uptake, distribution and elimination of (14)C-PFOA in zebrafish (Danio rerio)

Perfluorooctanoic acid (PFOA) is a long-chain perfluorinated chemical that has been shown to be non-degradable and persistent in the environment. Laboratory studies on bioconcentration and compound-specific tissue distribution in fish can be valuable for prediction of the persistence and environmental effects of the chemicals. In the present study male and female zebrafish (Danio rerio) were continuously exposed to 10μg/L of radiolabeled perfluorooctanoic acid ((14)C-PFOA) for 40 days, after which the exposed fish were transferred to fresh clean water for another 80 days wash-out period. At defined periodic intervals during the uptake and wash-out, fish were sampled for liquid scintillation counting and whole body autoradiography to profile the bioconcentration and tissue distribution of PFOA. The steady-state concentration of (14)C-PFOA in the zebrafish was reached within 20-30 days of exposure. The concentration-time course of (14)C-PFOA displayed a bi-exponential decline during washout, with a terminal half-life of approximately 13-14 days. At steady-state the bioconcentration of (14)C-PFOA into whole-body fish was approximately 20-30 times greater than that of the exposure concentration, with no differences between females and males. The bioconcentration factors for liver and intestine were approximately 100-fold of the exposure medium, while in brain, ovary and gall bladder the accumulation factors were in the range 15-20. Whole-body autoradiograms confirmed the highest labeling of PFOA in bile and intestines, which implies enterohepatic circulation of PFOA. The (14)C-PFOA was also observed in maturing vitellogenic oocytes, suggesting chemical accumulation via yolk proteins into oocytes with plausible risk for adverse effects on early embryonic development and offspring health. The bioconcentration at several (14)C-PFOA exposure concentrations were also investigated (0.3-30μg/L). This showed that bioconcentration increased linearly with tank exposure in the present in vivo model under steady-state conditions. From this model tissue concentrations of PFOA can be predicted when the external exposure level is known. The present study has generated experimental data on PFOA kinetics in zebrafish that can be valuable for aquatic environmental risk assessment.

Comparison of developmental toxicity of seven perfluoroalkyl acids to zebrafish embryos

The toxicity of individual perfluoroalkyl acids (PFAAs) has been suggested to be determined by the carbon chain length as well as the functional group attached. We tested seven different PFAAs including both sulfonic and carboxylic PFAAs with different chain length to evaluate the developmental toxicity in zebrafish embryos. Generally, the acute toxicity of PFAAs is relatively low to zebrafish embryos. The EC50 values ranged from 1.5 to 2200mg/L. We observed a relationship between higher toxicity with longer carbon chain. In addition, we also observed a higher toxicity for sulfonic PFAAs than for carboxylic PFAAs.