Physiologically based pharmacokinetic (PBPK) model for residues of lipophilic pesticides in poultry

A brief review on models for birds exposed to chemicals

"A Who's Who of pesticides is therefore of concern to us all. If we are going to live so intimately with these chemicals eating and drinking them, taking them into the very marrow of our bones - we had better know something about their nature and their power."-Rachel Carson, Silent Spring. In her day, Rachel Carson was right: plant protection products (PPP), like all the other chemical substances that humans increasingly release into the environment without further precaution, are among our worst enemies today (Bruhl and Zaller, 2019; Naidu et al., 2021; Tang et al., 2021; Topping et al., 2020). All compartments of the biosphere, air, soil and water, are potential reservoirs within which all species that live there are impaired. Birds are particularly concerned: PPP are recognized as a factor in the decline of their abundance and diversity predominantly in agricultural landscapes. Due to the restrictions on vertebrates testing, in silico-based approaches are an ideal choice alternative given input data are available. This is where the problem lies as we will illustrate in this paper. We performed an extensive literature search covering a long period of time, a wide diversity of bird species, a large range of chemical substances, and as many model types as possible to encompass all our future need to improve environmental risk assessment of chemicals for birds. In the end, we show that poultry species exposed to pesticides are the most studied at the individual level with physiologically based toxicokinetic models. To go beyond, with more species, more chemical types, over several levels of biological organization, we show that observed data are crucially missing (Gilbert, 2011). As a consequence, improving existing models or developing new ones could be like climbing Everest if no additional data can be gathered, especially on chemical effects and toxicodynamic aspects.

Interspecies scaling of toxicity reference values in human health versus ecological risk assessments: A critical review

Risk assessments that focus on anthropogenic chemicals in environmental media-whether considering human health or ecological effects-often rely on toxicity data from experimentally studied species to estimate safe exposures for species that lack similar data. Current default extrapolation approaches used in both human health risk assessments and ecological risk assessments (ERAs) account for differences in body weight between the test organisms and the species of interest, but the two default approaches differ in important ways. Human health risk assessments currently employ a default based on body weight raised to the three-quarters power. Ecological risk assessments for wildlife (i.e., mammals and birds) are typically based directly on body weight, as measured in the test organism and receptor species. This review describes differences in the experimental data underlying these default practices and discusses the many factors that affect interspecies variability in chemical exposures. The interplay of these different factors can lead to substantial departures from default expectations. Alternative methodologies for conducting more accurate interspecies extrapolations in ERAs for wildlife are discussed, including tissue-based toxicity reference values, physiologically based toxicokinetic and/or toxicodynamic modeling, chemical read-across, and a system of categorical defaults based on route of exposure and toxic mode of action. Integr Environ Assess Manag 2024;20:749-764. © 2023 SETAC.

Design of a generic model based on physiology for persistent organic pollutants in laying hens: Applications on chlordecone and chlorinated paraffins

A Physiology Based Pharmacokinetic (PBPK) model has been developed to predict the kinetics of Persistent Organic Pollutants (POPs) in laying hens. Different datasets have enabled the calibration of the model for chlordecone (CLD), an organochlorine pesticide used in the French West Indies between 1972 and 1993, as well as for chlorinated paraffins (CPs), widely used for various industrial applications worldwide. For this purpose, the sensitivity analysis showed that intake parameters, laying rate, partition coefficients of yolk, hepatic clearance, percentage of metabolism and age were key parameters. Applied to CLD and CPs, this model shows a good capacity for prediction, with 88 % of the experimental values ranging within 1.5-fold of the predicted value at steady state for CPs and 100 % for CLD. The fine modelling of the physiology and the laying process contributes to precision of the model and gives genericity, enabling the switch from one bird species to another. The model can be implemented with other POPs if the clearance and partition coefficient are known.

Determination of safe levels and toxic levels for feed hazardous materials in broiler chickens: a review

Feed safety is needed to produce and provide safe animal feeds for consumers, animals, and the environment. Although feed safety regulations have been set for each country, there is a lack of clear feed safety regulations for each livestock. Feed safety regulations are mainly focused on heavy metals, mycotoxins, and pesticides. Each country has different safe levels of hazardous materials in diets. Safe levels of hazardous materials in diets are mostly set for mixed diets of general livestock. Although there is a difference in the metabolism of toxic materials among animals, the safe level of feed is not specific for individual animals. Therefore, standardized animal testing methods and toxicity studies for each animal are needed to determine the correct safe and toxic levels of hazardous materials in diets. If this goal is achieved, it will be possible to improve livestock productivity, health, and product safety by establishing appropriate feed safety regulations. It will also provide an opportunity to secure consumer confidence in feed and livestock products. Therefore, it is necessary to establish a scientific feed safety evaluation system suitable for each country's environment. The chance of outbreaks of new hazardous materials is increasing. Thus, to set up appropriate toxic levels or safe levels in feed, various toxicity methods have been used to determine toxic levels of hazardous materials for humans and animals. Appropriate toxic testing methods should be developed and used to accurately set up and identify toxicity and safe levels in food and feed.

Development and Application of a Physiologically Based Pharmacokinetic Model for Diclazuril in Broiler Chickens

Withdrawal periods for diclazuril in broilers have traditionally been determined through regression analysis. However, over the last two decades, the physiologically based pharmacokinetic (PBPK) model has gained prominence as a predictive tool for veterinary drug residues, which offers an alternative method for establishing appropriate withdrawal periods for veterinary drugs. In this current study, a flow-limited PBPK model was developed to predict diclazuril concentrations in broilers following long-duration administration via medicated feed and water. This model consists of nine compartments, including arterial and venous plasma, lung, muscle, skin + fat, kidney, liver, intestine contents, and the rest of the body compartment. Physiological parameters such as tissue weights (V_cxx) and blood flow (Q_cxx) were gathered from published studies, and tissue/plasma partition coefficients (P_xx) were calculated through the area method or parameter optimization. Published diclazuril concentrations were compared to the predicted values, indicating the accuracy and validity of the model. The sensitivity analysis showed that parameters associated with cardiac output, drug absorption, and elimination significantly affected diclazuril concentrations in the muscle. Finally, a Monte Carlo analysis, consisting of 1000 iterations, was conducted to calculate the withdrawal period. Based on the Chinese MRL values, we calculated a withdrawal period of 0 days for both recommended dosing regimens (through mediated water and feed at concentrations of 0.5-1 mg/L and 1 mg/kg, respectively). However, based on the European MRLs, longer periods were determined for the mediated feed dosing route. Our model provides a foundation for scaling other coccidiostats and poultry species.

A critical review of effect modeling for ecological risk assessment of plant protection products

A wide diversity of plant protection products (PPP) is used for crop protection leading to the contamination of soil, water, and air, which can have ecotoxicological impacts on living organisms. It is inconceivable to study the effects of each compound on each species from each compartment, experimental studies being time consuming and cost prohibitive, and animal testing having to be avoided. Therefore, numerous models are developed to assess PPP ecotoxicological effects. Our objective was to provide an overview of the modeling approaches enabling the assessment of PPP effects (including biopesticides) on the biota. Six categories of models were inventoried: (Q)SAR, DR and TKTD, population, multi-species, landscape, and mixture models. They were developed for various species (terrestrial and aquatic vertebrates and invertebrates, primary producers, micro-organisms) belonging to diverse environmental compartments, to address different goals (e.g., species sensitivity or PPP bioaccumulation assessment, ecosystem services protection). Among them, mechanistic models are increasingly recognized by EFSA for PPP regulatory risk assessment but, to date, remain not considered in notified guidance documents. The strengths and limits of the reviewed models are discussed together with improvement avenues (multigenerational effects, multiple biotic and abiotic stressors). This review also underlines a lack of model testing by means of field data and of sensitivity and uncertainty analyses. Accurate and robust modeling of PPP effects and other stressors on living organisms, from their application in the field to their functional consequences on the ecosystems at different scales of time and space, would help going toward a more sustainable management of the environment. Graphical Abstract Combination of the keyword lists composing the first bibliographic query. Columns were joined together with the logical operator AND. All keyword lists are available in Supplementary Information at https://doi.org/10.5281/zenodo.5775038 (Larras et al. 2021).

Fipronil and fipronil sulfone in chicken: From in vitro experiments to in vivo PBK model predictions

In 2017 a large-scale fipronil contamination in eggs occurred in several European countries. Fipronil and its metabolites have the potential to be transferred into the eggs of laying hens, thereby entering the human food chain. Here, first the metabolism of fipronil was measured in vitro using chicken liver S9. The results show that fipronil is mainly metabolised into fipronil sulfone and the clearance obtained in vitro was extrapolated to in vivo liver clearance. In a second step a physiologically based kinetic model was developed with a focus on fipronil and its major sulfone metabolite and the model outcome was compared to available in vivo data in eggs from the literature. The experimentally obtained clearance was used as model input to evaluate whether such an in vitro-based model can be used in an early phase of a contamination incident to predict the time-concentration curves. Overall, all model predictions were within a 10-fold difference and the estimated elimination half-life for fipronil equivalents was 14 days. In vitro experiments are definitely recommended compared to in vivo studies, since they provide a fast first insight into the behaviour of a chemical in an organism.

Acute toxicity and metabolism of pesticides in birds

The median lethal dose of pesticide in acute oral toxicity, used as a conservative index in avian risk assessment, varies by the species with differences of less than one order of magnitude, depending on body size, feeding habit, and metabolic enzyme activity. The profiles of pesticide metabolism in birds with characteristic conjugations are basically common to those in mammals, but less information is available on their relevant enzymes. The higher toxicity of some pesticides in birds than in mammals is due to the lower activity of avian metabolic enzymes. The bioaccumulation in birds is limited for very hydrophobic pesticides resistant to metabolic degradation. Several in silico approaches using the descriptors of a pesticide molecule have recently been employed to estimate the profiles of acute oral toxicity and bioaccumulation.

Regulation of pesticide soil standards for protecting human health based on multiple uses of residential soil

To help environmental agencies manage pesticides in residential soil and reduce the associated risks to human health, we developed a screening-level framework that derives pesticide soil standards (PSSs) while considering the multiple uses of residential soil. Our screening models simulated the risk from exposure to soil pesticides via direct and three major indirect (i.e., tuber crops, animal-sourced food, and groundwater) exposure pathways. Based on these models, we derived PSSs for five types of residential soil. Our results showed that, in general, indirect pathways contributed more than the direct pathway to the overall exposure to soil pesticides. Consequently, in rural environments, where residential soil is also subjected to activities such as agriculture, animal grazing, and groundwater consumption, the derived PSSs were low. In addition, we compared the derived PSSs to the current worldwide standards for 13 commonly used pesticides. We found that the current global PSSs were appropriate only for urban residential soil. In many rural environments where the boundaries between different soil uses may be indistinct, the current PSSs are insufficient to protect humans from exposure to soil pesticides. Based on this analysis and the proposed PSSs, we provide regulatory recommendations for the management of pesticides in various types of residential soils.

In vitro metabolism of lidocaine in subcellular post-mitochondrial fractions and precision cut slices from cattle liver

In vitro models are widely used to study the biotransformation of xenobiotics and to provide input parameters to physiologically based kinetic models required to predict the kinetic behavior in vivo. For farm animals this is not common practice yet. The use of slaughterhouse-derived tissue material may provide opportunities to study biotransformation reactions in farm animals. The goal of the present study was to explore the potential of slaughterhouse-derived bovine liver S9 (S9) and precision cut liver slices (PCLSs) to capture observed biotransformation reactions of lidocaine in cows. The in vitro data obtained with both S9 and PCLSs confirm in vivo findings that 2,6-dimethylaniline (DMA) is an important metabolite of lidocaine in cows, being for both PCLSs and S9 the end-product. In case of S9, also conversion of lidocaine to lidocaine-N-oxide and monoethylglycinexylidine (MEXG) was observed. MEGX is considered as intermediate for DMA formation, given that this metabolite was metabolized to DMA by both PLCSs and S9. In contrast to in vivo, no in vitro conversion of DMA to 4-OH-DMA was observed. Further work is needed to explain this lack of conversion and to further evaluate the use of slaughterhouse-derived tissue materials to predict the biotransformation of xenobiotics in farm animals.

Physiological parameter values for physiologically based pharmacokinetic models in food-producing animals. Part II: Chicken and turkey

Physiologically based pharmacokinetic (PBPK) models are growing in popularity due to human food safety concerns and for estimating drug residue distribution and estimating withdrawal intervals for veterinary products originating from livestock species. This paper focuses on the physiological and anatomical data, including cardiac output, organ weight, and blood flow values, needed for PBPK modeling applications for avian species commonly consumed in the poultry market. Experimental and field studies from 1940 to 2019 for broiler chickens (1-70 days old, 40 g - 3.2 kg), laying hens (4-15 months old, 1.1-2.0 kg), and turkeys (1 day-14 months old, 60 g -12.7 kg) were searched systematically using PubMed, Google Scholar, ProQuest, and ScienceDirect for data collection in 2019 and 2020. Relevant data were extracted from the literature with mean and standard deviation (SD) being calculated and compiled in tables of relative organ weights (% of body weight) and relative blood flows (% of cardiac output). Trends of organ or tissue weight growth during different life stages were calculated when sufficient data were available. These compiled data sets facilitate future PBPK model development and applications, especially in estimating chemical residue concentrations in edible tissues to calculate food safety withdrawal intervals for poultry.

An open source physiologically based kinetic model for the chicken (Gallus gallus domesticus): Calibration and validation for the prediction residues in tissues and eggs

Xenobiotics from anthropogenic and natural origin enter animal feed and human food as regulated compounds, environmental contaminants or as part of components of the diet. After dietary exposure, a chemical is absorbed and distributed systematically to a range of organs and tissues, metabolised, and excreted. Physiologically based kinetic (PBK) models have been developed to estimate internal concentrations from external doses. In this study, a generic multi-compartment PBK model was developed for chicken. The PBK model was implemented for seven compounds (with log K_ow range -1.37-6.2) to quantitatively link external dose and internal dose for risk assessment of chemicals. Global sensitivity analysis was performed for a hydrophilic and a lipophilic compound to identify the most sensitive parameters in the PBK model. Model predictions were compared to measured data according to dataset-specific exposure scenarios. Globally, 71% of the model predictions were within a 3-fold change of the measured data for chicken and only 7% of the PBK predictions were outside a 10-fold change. While most model input parameters still rely on in vivo experiments, in vitro data were also used as model input to predict internal concentration of the coccidiostat monensin. Future developments of generic PBK models in chicken and other species of relevance to animal health risk assessment are discussed.

A PBPK model to study the transfer of α-hexabromocyclododecane (α-HBCDD) to tissues of fast- and slow-growing broilers

A physiologically based pharmacokinetic (PBPK) model was developed to investigate the production-specific factors involved in the transfer of α-hexabromocyclododecane (α-HBCDD) to broiler meat. The model describes growth and lipid deposition in tissues of fast- (FG) and slow- (SG) growing broilers from hatching to slaughter and simulates the exposure through the ingestion of contaminated feed or expanded polystyrene insulation material. Growth parameters were obtained from the literature while parameters relative to uptake, distribution, and elimination of α-HBCDD were adjusted using results of a previous experiment involving broilers exposed through feed throughout the rearing period or allowed to depurate before slaughter. The model was used to compare the two main edible tissues, breast and leg meat, as well as skin, and to investigate the variability within strain. Between strains and within strain, α-HBCDD assimilation efficiency (AE) is higher when the animals are slaughtered young or heavy. However, increasing slaughter age will lower α-HBCDD concentration in tissues, due to dilution. Based on fresh weight, the concentration of α-HBCDD in breast muscles and skin tends to be lower in SG than in FG broilers (-30 to +10%), while it is 10% to 80% higher in leg muscles. Compared to breast muscles, consuming leg muscles would elicit an exposure 9 and 16 times higher in FG and SG broilers, respectively. The consumption of skin together with muscles would multiply the exposure by up to 36 times compared to breast muscle alone. In case of acute exposure, the α-HBCDD concentration in tissues increased sharply, all the more since the animals are lighter in weight, and then decreased rapidly. In FG broilers, dilution through growth contributed for up to 37%, 28% and 97% to the decontamination of breast muscles, leg muscles and skin, respectively, depending on the duration of depuration before slaughter.

Assessing Global Human Exposure to T-2 Toxin via Poultry Meat Consumption Using a Lifetime Physiologically Based Pharmacokinetic Model

Residue depletion of T-2 toxin in chickens after oral gavage at 2.0 mg/kg twice daily for 2 days was determined in this study. A flow-limited physiologically based pharmacokinetic (PBPK) model was developed for lifetime exposure assessment in chickens. The model was calibrated with data from the residue depletion study and then validated with independent data. A local sensitivity analysis was performed, and 16 sensitive parameters were subjected to Monte Carlo analysis. The population PBPK model was applied to estimate daily intake values of T-2 toxin in different countries based on reported consumption factors and the guidance value of 0.25 mg/kg in feed for chickens by the European Food Safety Authority (EFSA). The predicted daily intakes in different countries were all lower than the EFSA's total daily intake, suggesting that the EFSA's guidance value has minimal risk. This model provides a foundation for scaling to other mycotoxins and other food animal species.

Accumulation of α-hexabromocyclododecane (α-HBCDD) in tissues of fast- and slow-growing broilers (Gallus domesticus)

The aim of the current study was to describe the fate of ingested α-hexabromocyclododecane (α-HBCDD) in fast-growing (FG) and slow-growing (SG) broilers, through an exposure to a dietary concentration of 50 ng α-HBCDD g^-1 feed during 42 and 84 days, respectively. Depuration parameters were assessed in SG broilers successively exposed during 42 days and depurated during 42 days. At market age, SG broilers had ingested 42% more feed than FG broilers, while their body weight gain per g of feed ingested was 34% lower. No isomerization of α- to β- or γ-HBCDD forms occurred, while OH-HBCDD was identified as a product of α-HBCDD metabolism. Irrespective of the strain, abdominal fat displayed the highest α-HBCDD concentration on a lipid weight basis, followed leg muscles and then breast muscle, liver and plasma. The accumulation ratios of α-HBCDD were slightly higher in SG (6.7, 2.1, 2.6 and 9.9 in leg muscles, breast muscle, liver and abdominal fat, respectively) than in FG broilers (5.2, 2.2, 1.1 and 8.4, respectively). The elimination half-lives in SG broilers were 20, 12 and 19 d in leg muscles, breast muscle and abdominal fat, respectively, to which dilution through growth contributed for around 50%. The overall assimilation efficiency of α-HBCDD was estimated at 58 and 50% in FG and SG broilers, respectively, while 22 and 17% of α-HBCDD ingested were estimated to be eliminated in excreta as metabolites.

Modeling PCB transfer into hen eggs: influence of physiological characteristics of the animal

Laying hens are likely to be exposed to a wide range of persistent organic pollutants (POPs), including polychlorinated biphenyls (PCBs). To improve the safety of poultry farming systems in terms of POPs, the present research focused on assessing the impact of physiological characteristics of the hen on the transfer of ingested PCBs to eggs. Modeling was used as a research tool to explore the impact of some physiological characteristics on the transfer of PCBs in the laying hen. The mathematical model simulates the dynamics of the size of the lipid compartments in the animal and the frequency of laying, with the PCB concentrations in egg yolk and adipose tissue being model outputs. Simulations were run to assess effects of animal characteristics on the transfer of PCBs to eggs. Laying rate proportionally influenced the PCB level of eggs and adipose tissue at steady state. Body fat diluted absorbed PCBs in the absence of laying and significantly influenced the decontamination rate of tissues during depuration after an exposure period. Application of the present model to actual exposure cases highlights its value in improving the support of risk management in livestock farming.

Estimating tulathromycin withdrawal time in pigs using a physiologically based pharmacokinetics model

A physiologically based pharmacokinetics model was developed to predict tulathromycin concentrations in edible swine tissues. Physiological parameters included volumes of and plasma flows through different tissues which were obtained from the literatures. The tissue/plasma partition coefficient was calculated according to the area method, and the model was validated through a comparison of predicted and observed concentrations. Withdrawal times in different tissues were predicted. The physiologically based pharmacokinetics model presented here provided accurate predictions of the observed concentrations in all tissues. The results showed that the injection site had the longest withdrawal time (21 days), followed by skin together with fat (19 days) and then kidney (10 days), lung (6 days), liver (4 days) and muscle (1 day). A withdrawal time of 21 days was finally predicted for tulathromycin in swine after a single intramuscular injection at 2.5 mg/kg body weight.

A physiologically based pharmacokinetic model for the prediction of the depletion of methyl-3-quinoxaline-2-carboxylic acid, the marker residue of olaquindox, in the edible tissues of pigs

To estimate the consumer exposure to olaquindox (OLA) residues in porcine edible tissues, a physiologically based pharmacokinetic (PBPK) model for methyl-3-quinoxaline-2-carboxylic acid (MQCA), the marker residue of OLA, was developed in pigs based on the assumptions of the flow-limited distribution, hepatic metabolism, and renal excretion. The model included separate compartments corresponding to blood, muscle, liver, kidney, adipose, and an extra compartment representing the remaining carcass. Physiological parameters were determined from literatures. Plasma protein binding, partition coefficients, and renal clearance for MQCA were determined in in vitro and in vivo studies. The metabolic conversion of OLA to MQCA was assumed as a simple, one-step process, and an apparent first-order rate constant (k) was employed to describe this metabolic process. The PBPK model was optimized and validated with plasma and tissue data from literatures and our study. Sensitivity analysis and Monte Carlo simulation were also implemented to estimate the influence of model parameters on the goodness of fit. When compared with the observed data, the PBPK model underestimated the MQCA level in all compartments at the early time points, whereas gave excellent predictions of MQCA concentration in porcine edible tissues at later time points. The correlation coefficients between the predicted and observed values were over 0.88. The consistency between the model predictions and the real residues of OLA in pigs proved the good applicability of our model in food safety risk assessment.