A dose-based modeling approach for accumulation and toxicity of arsenic in tilapia Oreochromis mossambicus

Arsenic bioaccumulation and biotransformation in different tissues of Nile tilapia (Oreochromis niloticus): A comparative study between As(III) and As(V) exposure and evaluation of antagonistic effects of selenium

The speciation of arsenic in fish has been widely investigated, but bioaccumulation and biotransformation of inorganic As in different tissues of Nile tilapia (Oreochromis niloticus) are not fully understood. The present study aimed to investigate the bioaccumulation of As in Nile tilapia, as well as to evaluate the distribution of the main arsenic species (As(III), As(V), MMA, DMA, and AsB) in liver, stomach, gill, and muscle, after controlled exposures to As(III) and As(V) at concentrations of 5.0 and 10.0 mg L-1 during periods of 1 and 7 days. Total As was determined by inductively coupled plasma mass spectroscopy (ICP-MS). For both exposures (As(III) and As(V)), the total As levels after 7-day exposure were highest in the liver and lowest in the muscle. Overall, the Nile tilapia exposed to As(III) showed higher tissue levels of As after the treatments, compared to As(V) exposure. Speciation of arsenic present in the tissues employed liquid chromatography coupled to ICP-MS (LC-ICP-MS), revealing that the biotransformation of As included As(V) reduction to As(III), methylation to monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA), and subsequent conversion to nontoxic arsenobetaine (AsB), which was the predominant arsenic form. Finally, the interactions and antagonistic effects of selenium in the bioaccumulation processes were tested by the combined exposure to As(III), the most toxic species of As, together with tetravalent selenium (Se(IV)). The results indicated a 4-6 times reduction of arsenic toxicity in the tilapia.

Intestinal uptake and low transformation increase the bioaccumulation of inorganic arsenic in freshwater zebrafish

Arsenate [As(V)] is the main form of arsenic (As) present in freshwater taken up by freshwater fish. Data on the main uptake tissue, biotransformation, and bioaccumulation in freshwater fish exposed to As(V) were limited, and the reasons for its bioaccumulation in the muscle tissue of freshwater fish remain undetermined. Accordingly, we simulated bioaccumulation and depuration in zebrafish (Danio rerio) exposed to waterborne As(V) by employing a six-compartment physiologically based pharmacokinetic model and As speciation analysis. Modeling and biotransformation suggested that intestines were the main uptake site for waterborne As(V), instead of the gills. This novel finding was evidenced by the higher As transfer constant from water to intestines (k₀₃ = 1.52 × 10^-4 L d^-1) compared to gills (k₀₂ = 5.28 × 10^-5 L d^-1). The low concentration and percentage of arsenobetaine (AsB) in the intestines suggested a weak ability to synthesize AsB. Our results showed a substantial proportion of inorganic As in intestines and a relatively substantial percentage in muscle tissue. Therefore, high As(V) uptake in the intestines and lack of biotransformation contributed to high bioaccumulation of inorganic As in freshwater fish. Inorganic As posed concerns due to the human health risks associated with consuming As(V)-contaminated fish and should be addressed.

An Arsenic-76 radiotracer to study the routes of assimilation, hemolymph distribution, and tissue inventories in the bioindicator organism Pomacea canaliculata

The aim of this work was to study the absorption, distribution through the hemolymph, and bioaccumulation of arsenic by the freshwater Pomacea canaliculata using a short-lived tracer (⁷⁶As, t_1/2: 1.07 d) with high specific activity. Arsenic travels mainly dissolved in the plasma of the snail's hemolymph. This element is transferred from the hemolymph to the tissues (87%) 4 h after the inoculation of 50 μL of a 0.04 g/L of ⁷⁶As radiotracer solution, being the digestive gland, kidney, and head-foot the main places of arsenical inventories. Snails exhibited a rapid arsenic accumulation response in a wide range of concentrations (from 1 to 1000 μg/L) of the metalloid dissolved in water and in a concentration-dependent manner. Also, snails incorporated As from the digestive system when they received a single safe dose of ~2 μg of ⁷⁶As inoculated in a fish food pellet. The (semi) physiologically based toxicokinetic model developed in this study is based on anatomical and physiological parameters (blood flow, irrigation, tissue volume and other). Together, these findings make P. canaliculata an excellent sentinel organism to evaluate freshwater bodies naturally contaminated with As.

Physiologically Based Pharmacokinetic Model for the Biotransportation of Arsenic in Marine Medaka (Oryzias melastigma)

The toxicity of arsenic (As) targets specific tissues of organisms, while the biotransportation of As among the tissues of fish remains poorly understood. In the present study, radiotracer techniques followed by a physiologically based pharmacokinetic (PBPK) modeling were applied to simulate the biotransportation (absorption, distribution, and elimination) of ⁷³As(V) and biotransformation of As(V) in the marine medaka Oryzias melastigma after waterborne As exposure. Fish were simulated by a six-compartment model by assuming that blood was the intermediate exchange among different compartments (gill, intestine, liver, head, and carcass). Modeling suggested that intestine and gill were the uptake, exchange, as well as elimination sites of waterborne As, while carcass and head were the main storage sites. Intestine played a vital role in the metabolism of As(V) by biotransforming inorganic As into arsenobetaine (AsB), possibly because of the important role of gut microbiota. The correlation between the PBPK model constants and the As speciation (e.g., AsB %, inorganic As %, and methylated As %) indicated that AsB tended to be stored in the tissues rather than being depurated, while inorganic and methylated As were more easily transferred from tissues to the blood and eliminated. Modeling simulation coupling with biotransformation for the first time demonstrated that the fish intestine was the main metabolic site, and synthesis of AsB as mediated by the microbiota in the intestine contributed to the high As bioaccumulation in marine fish.

Applications of dynamic models in predicting the bioaccumulation, transport and toxicity of trace metals in aquatic organisms

This review evaluates the three dynamic models (biokinetic model: BK, physiologically based pharmacokinetic model: PBPK, and toxicokinetic-toxicodynamic model: TKTD) in our understanding of the key questions in metal ecotoxicology in aquatic systems, i.e., bioaccumulation, transport and toxicity. All the models rely on the first-order kinetics principle of metal uptake and elimination. The BK model basically treats organisms as a single compartment, and is both physiologically and geochemically based. With a good understanding of each kinetic parameter, bioaccumulation of metals in any aquatic organisms can be studied holistically and mechanistically. Modeling efforts are not merely restrained from the prediction of metal accumulation in the tissues, but instead provide the direction of the key processes that need to be addressed. PBPK is more physiologically based since it mainly addresses the transportation, transformation and distribution of metals in the organisms. It can be treated conceptually as a multi-compartmental kinetic model, whereas the physiology is driving the development of any good PBPK model which is no generic for aquatic animals and contaminants. There are now increasingly applications of the PBPK modeling specifically in metal studies, which reveal many important processes that are impossible to be teased out by direct experimental measurements without adequate modeling. TKTD models further focus on metal toxicity in addition to metal bioaccumulation. The TK part links exposure and bioaccumulation, while the TD part links bioaccumulation and toxic effects. The separation of TK and TD makes it possible to model processes, e.g., toxicity modification by environmental factors, interaction between different metals, at both the toxicokinetic and toxicodynamic levels. TKTD models provide a framework for making full use of metal toxicity data, and thus provide more information for environmental risk assessments. Overall, the three models reviewed here will continue to provide guiding principles in our further studies of metal bioaccumulation and toxicity in aquatic organisms.

The effect of thermal pre-incubation and exposure on sensitivity of zebrafish (Danio rerio) to copper and cadmium single and binary exposures

Zebrafish (Danio rerio) is a prominent model organism in a wide range of biological studies including toxicology. However, toxicological studies are often performed at species specific optimum temperature, and knowledge on the effect of different temperature regimes on the toxicity of metal ions is rather limited. To address this knowledge gap, present study investigates the effect of various thermal scenarios (simultaneous and sequential; acute and chronic) on the toxicity of Cu and Cd in zebrafish. For this purpose we assessed mortality and whole body metal burdens as indicators of toxicity and bioavailability, respectively, and whole body electrolyte concentrations and body condition as the indicators of physiological condition. Thermal pre-incubations (for 12 or 96 h or 28 days) and subsequent metal ion exposures (for 10 days) were conducted at 17, 22, 25, 28, 32 and 34 °C. The metal exposures were performed at Cu concentrations of 1.2 μM and Cd concentrations of 0.2 μM, both singly and in binary mixtures. Irrespective of thermal treatments, Cu exposures resulted in greater mortality than Cd exposures at the given concentrations. Moreover, the Cu and Cd mixture indicated a synergistic effect. While acute pre-incubation for 12 or 96 h at elevated temperatures increased mortality in the subsequent metal exposure at the optimum temperature (28 °C), pre-incubation at cold temperatures in this scenario appeared to increase tolerance towards the subsequent metal exposure. Chronic thermal pre-incubation of zebrafish to a range of temperatures for 28 days moderated the effect of temperature fluctuations on subsequent metal toxicity at the optimum temperature. Chronic thermal pre-incubation at a range of temperatures followed by metal exposure at the same temperature showed that environmental temperature variations (higher or lower than optimal temperature) coupled with metal exposure, led to increased mortality, furthermore, the highest whole body metal burdens were measured in this scenario. Nevertheless, neither the whole body burden of metals, nor the metal accumulation rate, were predictors of mortality, i.e. these two values were not higher in dead fish in comparison to those that survived the exposures. Finally, we observed a significant decrease in the whole body Na⁺ level of dead fish in comparison to fish which survived the exposure conditions, suggesting that survival depends on maintaining Na⁺ homeostasis under the applied multi-stress conditions. Overall, our results show that thermal pre-history and ambient temperature play an important role in determining the tolerance of zebrafish towards metal ion stress.

Bioaccumulation and acute toxicity of As(III) and As(V) in Nile tilapia (Oreochromis niloticus)

The distribution of arsenic in the different tissues of tilapia fish is determined by the exposure time and the depuration rate. The mechanisms of toxicity/carcinogenicity depend on the arsenic species involved in the biotransformation processes. After a 7-day exposure period, the accumulation of inorganic arsenic (iAs) in the tilapia tissues studied was in the order: liver > stomach > gills > muscles. In bioaccumulation assays, the values of the organ uptake constant (k_a) ranged from 0.06 to 0.51 mL g^-1 d^-1, while the depuration rate constant (k_d) values were in the range 0.03-1.15 d^-1. Higher iAs bioaccumulation factor (BCF) values were observed for the stomach (3.1 mL g^-1) and the liver (1.6 mL g^-1), reflecting their high capacity to accumulate iAs species. These organs act as long-term storage sites for iAs, following chronic exposure. The LC₅₀ values were determined considering the average iAs concentration and the cumulative fish mortality. For As(III), the LC₅₀ values indicated fish mortality at concentrations above 30 mg L^-1. The fish showed greater tolerance to exposure to As(V), compared to As(III), with fish mortality after the second day of exposure requiring an As(V) concentration 7-fold higher than As(III).

Toxicokinetic models and related tools in environmental risk assessment of chemicals

Environmental risk assessment of chemicals for the protection of ecosystems integrity is a key regulatory and scientific research field which is undergoing constant development in modelling approaches and harmonisation with human risk assessment. This review focuses on state-of-the-art toxicokinetic tools and models that have been applied to terrestrial and aquatic species relevant to environmental risk assessment of chemicals. Both empirical and mechanistic toxicokinetic models are discussed using the results of extensive literature searches together with tools and software for their calibration and an overview of applications in environmental risk assessment. These include simple tools such as one-compartment models, multi-compartment models to physiologically-based toxicokinetic (PBTK) models, mostly available for aquatic species such as fish species and a number of chemical classes including plant protection products, metals, persistent organic pollutants, nanoparticles. Data gaps and further research needs are highlighted.

Homeostatic regulation of copper in a marine fish simulated by a physiologically based pharmacokinetic model

Copper (Cu) is an essential yet potentially toxic metal, thus delicate homeostatic controls are developed in the fish. In this study, a physiologically based pharmacokinetic (PBPK) model was developed to simulate the homeostatic regulation of Cu in a marine fish (Terapon jarbua) under dietary and waterborne exposures. In this model, fish were schematized as a six-compartment model, with the intestine being divided into two sub-compartments (chyme and gut wall). The blood was assumed to be the "carrier" distributing Cu into different compartments. The transfer rates between different compartments were determined in fish during Cu exposure (20 d) and depuration (20 d). The differences in Cu transfer from chyme to gut wall between dietary and waterborne treatments suggested that the intestine regulated the dietary uptake and re-absorption of Cu from the chyme. The extremely low uptake rate constant (0.0013 d^-1) for gills under waterborne exposure indicated that gills strongly restricted Cu uptake from the ambient water. For both treatments, the liver had considerable input rate through the enterohepatic circulation and comparably high exchange rate with the blood, suggesting that the liver can efficiently accumulate newly absorbed Cu. The differences in Cu output from the liver between dietary and waterborne treatments suggested that it can effectively regulate the redistribution of Cu. All of these observations demonstrated that the liver played the central role in Cu homeostasis by serving as the main depository and distributing center. Modeling results also indicated that renal and branchial excretion was of minor importance, whereas biliary excretion combined with defecation played the most important role in whole-body Cu elimination in marine fish. The effective regulation by the "Blood-Liver-Intestine" cycle could be the main reason for the relatively low levels of Cu in fish.

Development and Validation of a Biodynamic Model for Mechanistically Predicting Metal Accumulation in Fish-Parasite Systems

Because of different reported effects of parasitism on the accumulation of metals in fish, it is important to consider parasites while interpreting bioaccumulation data from biomonitoring programmes. Accordingly, the first step is to take parasitism into consideration when simulating metal bioaccumulation in the fish host under laboratory conditions. In the present study, the accumulation of metals in fish-parasite systems was simulated by a one-compartment toxicokinetic model and compared to uninfected conspecifics. As such, metal accumulation in fish was assumed to result from a balance of different uptake and loss processes depending on the infection status. The uptake by parasites was considered an efflux from the fish host, similar to elimination. Physiological rate constants for the uninfected fish were parameterised based on the covalent index and the species weight while the parameterisation for the infected fish was carried out based on the reported effects of parasites on the uptake kinetics of the fish host. The model was then validated for the system of the chub Squalius cephalus and the acanthocephalan Pomphorhynchus tereticollis following 36-day exposure to waterborne Pb. The dissolved concentration of Pb in the exposure tank water fluctuated during the exposure, ranging from 40 to 120 μg/L. Generally, the present study shows that the one-compartment model can be an effective method for simulating the accumulation of metals in fish, taking into account effects of parasitism. In particular, the predicted concentrations of Cu, Fe, Zn, and Pb in the uninfected chub as well as in the infected chub and the acanthocephalans were within one order of magnitude of the measurements. The variation in the absorption efficiency and the elimination rate constant of the uninfected chub resulted in variations of about one order of magnitude in the predicted concentrations of Pb. Inclusion of further assumptions for simulating metal accumulation in the infected chub led to variations of around two orders of magnitude in the predictions. Therefore, further research is required to reduce uncertainty while characterising and parameterising the model for infected fish.

Toxicokinetics/toxicodynamics with damage feedback improves risk assessment for tilapia and freshwater clam exposed to arsenic

It has been proposed that irreversible responses of organisms exposed to contaminants are due to a systems-level feedback. Here we tested this hypothesis by reanalyzing the published data on toxicokinetics and survival probability based on a systems-level threshold damage model (TDM) incorporating with a positive damage feedback to explore the steady-state response and dynamic behavior of damage for tilapia and freshwater clam exposed to waterborne arsenic (As). We found that ultrasensitivity appeared in As-tilapia and freshwater clam systems with Hill coefficient n ≥ 4, indicating that the positive damage feedback mechanism has been triggered. We confirmed that damage can trigger a positive feedback loop that together with As stressor increases irreversibility. This study also showed that TDM with positive feedback gave a much better predictability than that of TDM at As concentrations ranging from 100 to 500 mg l(-1) for freshwater clam, whereas for tilapia, two models had nearly same performance on predictability. We suggested that mortality-time profile derived Hill coefficient could be used as a new risk indicator to assess the survival probability for species exposed to waterborne metals. We anticipated that the proposed toxicokinetics/toxicodynamics with a positive damage feedback may facilitate our understanding and manipulation of complex mechanisms of metal susceptibility among species and improve current risk assessment strategies.

Bioavailability links mode of action can improve the long-term field risk assessment for tilapia exposed to arsenic

The objective of this paper was to develop a mechanistic-based framework to explicitly incorporate the factors controlling the bioavailability, toxicodynamics and mode of action to enhance predictive ability of arsenic (As) toxicity to protect the health of farmed tilapia Oreochromis mossambicus. We linked the biotic ligand model and damage assessment model to develop a toxicokinetic model for elucidating the site-specific temporal changes of As bioavailability and to characterize how the fish regulate the metal toxicity. We built a bioavailability-mode of action-based growth toxicity model by linking a bioenergetic growth model and damage assessment model to predict how the As affects on the tilapia growth in the entire life span in site-specific field ecosystems. Here we show that the proposed model well describes the water-chemistry-dependent toxicokinetics and toxicodynamics variations of As to tilapia. We selected two local tilapia farms with different water chemistries located at southwestern Taiwan coast region to implement the proposed algorithm to predict the risk of As exposure. Results indicate that the growth toxicity of O. mossambicus in Taihsi is more sensitive than that in Peimen. We found that the effect of ion competition on the As bioavailability and their ecotoxicological effects on tilapia are more obvious in Taihsi comparing with that in Peimen. We suggested that the proposed bioavailability- and mode of action-based framework can be used to capture the biological response and regulation of tilapia to As exposures. It is applicable for a site-specific and long-term ecotoxicological risk assessment.

Sodium Gill Potential as a Tool to Monitor Valve Closure Behavior in Freshwater Clam Corbicula fluminea in Response to Copper

Valve closure behavior in freshwater clam Corbicula fluminea is a biologically sensitive endpoint. The purpose of this paper was to derive an electrophysiological response model of C. fluminea to assess copper (Cu)–sodium (Na) interactions in gill membrane, whereby valve closure behavior and Cu toxicity could be monitored. The proposed model was based on the integration of Cu bioavailability, Na and Cu internalizations, and electrochemically-based gill potentials. Based on Na active transport under non-equilibrium conditions, predicted gill potential of −8.2 mV agreed reasonably well with published the measured transepithelial potential of −7 mV in C. fluminea. Our proposed framework captured the general features observed in model applications including: (i) 50% inhibitory Cu²⁺ activities for Na membrane potential (E_Na) and uptake rate (J_Na) were estimated to be 0.072 and 0.043 μM, respectively, with a stoichiometry of 3Cu²⁺: 1E_Na and 1J_Na; (ii) the external Cu²⁺–dependent internal Na concentration could be parsimoniously estimated, and (iii) the site-specific clam gill potentials could be monitored. Here we provided a new approach to monitor waterborne metal toxicity to reduce the nationwide economic losses due to bans on harvesting of contaminated clam and the potential risks to the health of clams.

Arsenic cancer risk posed to human health from tilapia consumption in Taiwan

Ingested inorganic arsenic is strongly associated with a wide spectrum of adverse health outcomes. We propose a bioaccumulation and the Weibull model-based epidemiological framework to accurately estimate the reference arsenic intake guideline for tilapia consumption and tilapia-cultured water arsenic concentration based on bioaccumulations of tilapia and gender/age/cancer-specific epidemiological data from the arseniasis-endemic area in Taiwan. Our results show a positive relationship between arsenic exposure and age/gender- and cancer-specific cumulative incidence ratio using Weibull dose-response model. Based on male bladder cancer with an excess lifetime cancer risk of 10(-4), we estimate the reference tilapia inorganic arsenic guideline value to be 0.084 microg g(-1) dry wt based on the suggested daily consumption rate of 120 gd(-1). Our findings show that consumption of tilapia in a blackfoot disease (BFD)-endemic area poses no significant cancer risk (excess cancer risks ranging from 3.4 x 10(-5) to 9.3 x 10(-5)), implying that people in BFD-endemic areas are not readily associated with higher fatalities for bladder cancer exposed from tilapia consumption. We are confident that our model can be easily adapted for other aquaculture species, and encourage risk managers to use the model to evaluate the potential population-level long-term low-dose cancer risks. We conclude that, by integrating the bioaccumulation concept and epidemiological investigation of humans exposed to arsenic, we can provide a scientific basis for risk analysis to enhance risk management strategies.

Effects of lead, molybdenum, rubidium, arsenic and organochlorines on spermatogenesis in fish: monitoring at Mekong Delta area and in vitro experiment

To estimate the influence of water contaminants on fish reproduction in the Mekong Delta area, we sampled cultivated male catfish (Pangasianodon hypophthalmus), investigated testicular development, and measured persistent organic pollutants (POPs) and trace element levels in muscle and liver, respectively. Various testes sizes were observed although sampling took place during a short period. Histological analysis revealed that all developmental stages of germ cells were observed in catfish with large testis, whereas only necrotic spermatogonia but no other germ cells were observed in catfish with small testis. In small testis, furthermore, vacuolization and hypertrophy of Sertoli cells were observed. Measurement of POPs in muscle and trace elements in liver demonstrated that there were negative correlations between GSI and the concentrations of Pb, Mo, Rb and As. To clarify possible direct effects of Pb, Mo, Rb and As on spermatogenesis in fish, we investigated the effects of these trace elements on spermatogenesis using in vitro testicular organ culture of Japanese eel (Anguilla japonica). Treatment with each of the trace elements alone did not affect spermatogenesis. However, treatment with 10(-7)M of Pb, 10(-5) and 10(-4)M of Mo, 10(-5)-10(-3)M of Rb or 10(-5)M of As inhibited the spermatogenesis induced by 11-ketotestosterone (11KT). Furthermore, treatment with 10(-4)M of As in combination with 11KT caused necrosis of testicular fragments. Taken together, these results are consistent with the hypothesis that Pb, Mo, Rb and As can exert inhibitory effects on spermatogenesis in catfish inhabiting the Mekong Delta area.

Linking valve closure behavior and sodium transport mechanism in freshwater clam Corbicula fluminea in response to copper

The purpose of this study is to develop a mechanistic model to describe a conceptually new "flux-biological response" approach based on biotic ligand model (BLM) and Michaelis-Menten (M-M) kinetics to allow the linkage between valve closure behavior and sodium (Na) transport mechanism in freshwater clam Corbicula fluminea in response to waterborne copper (Cu). We test the proposed model against published data regarding Na uptake kinetics in rainbow trout and Na uptake profile in C. fluminea, confirming that the predictive model is robust. Here, we show that the predicted M-M maximum Cu internalization flux in C. fluminea is 0.369 micromolg(-1)h(-1) with a half-saturation affinity constant of 7.87x10(-3) microM. Dynamics of Na uptake and valve closure daily rhythm driven by external Cu can also be predicted simultaneously. We suggest that this "Na transport-valve closure behavior" approach might provide the basis of a future design of biomonitoring tool.

Toxicokinetics/toxicodynamics of arsenic for farmed juvenile milkfish Chanos chanos and human consumption risk in BFD-endemic area of Taiwan

This paper presents a toxicokinetic/toxicodynamic analysis to appraise arsenic (As) bioaccumulation in farmed juvenile milkfish Chanos chanos at blackfoot disease (BFD)-endemic area in Taiwan, whereas probabilistic incremental lifetime cancer risk (ILCR) and hazard quotient (HQ) models are also employed to assess the range of exposures for the fishers and non-fishers who eat the contaminated fish. We conducted a 7-day exposure experiment to obtain toxicokinetic parameters, whereas a simple critical body burden toxicity model was verified with LC50(t) data obtained from a 7-day acute toxicity bioassay. Acute toxicity bioassay indicates that 96-h LC50 for juvenile milkfish exposed to As is 7.29 (95% CI: 3.10-10.47) mg l(-1). Our risk analysis for milkfish reared in BFD-endemic area indicates a low likelihood that survival is being affected by waterborne As. Human risk analysis demonstrates that 90%-tile probability exposure ILCRs for fishers in BFD-endemic area have orders of magnitude of 10(-3), indicating a high potential carcinogenic risk, whereas there is no significant cancer risk for non-fishers (ILCRs around 10(-5)). All predicted 90%-tiles of HQ are less than 1 for non-fishers, yet larger than 10 for fishers which indicate larger contributions from farmed milkfish consumptions. Sensitivity analysis indicates that to increase the accuracy of the results, efforts should focus on a better definition of probability distributions for milkfish daily consumption rate and As level in milkfish. Here we show that theoretical human health risks for consuming As-contaminated milkfish in the BFD-endemic area are alarming under a conservative condition based on a probabilistic risk assessment model.