A non-invasive approach to study lifetime exposure and bioaccumulation of PCBs in protected marine mammals: PBPK modeling in harbor porpoises

Recent Progress on Physiologically Based Pharmacokinetic (PBPK) Model: A Review Based on Bibliometrics

Physiologically based pharmacokinetic/toxicokinetic (PBPK/PBTK) models are designed to elucidate the mechanism of chemical compound action in organisms based on the physiological, biochemical, anatomical, and thermodynamic properties of organisms. After nearly a century of research and practice, good results have been achieved in the fields of medicine, environmental science, and ecology. However, there is currently a lack of a more systematic review of progress in the main research directions of PBPK models, especially a more comprehensive understanding of the application in aquatic environmental research. In this review, a total of 3974 articles related to PBPK models from 1996 to 24 March 2024 were collected. Then, the main research areas of the PBPK model were categorized based on the keyword co-occurrence maps and cluster maps obtained by CiteSpace. The results showed that research related to medicine is the main application area of PBPK. Four major research directions included in the medical field were "drug assessment", "cross-species prediction", "drug-drug interactions", and "pediatrics and pregnancy drug development", in which "drug assessment" accounted for 55% of the total publication volume. In addition, bibliometric analyses indicated a rapid growth trend in the application in the field of environmental research, especially in predicting the residual levels in organisms and revealing the relationship between internal and external exposure. Despite facing the limitation of insufficient species-specific parameters, the PBPK model is still an effective tool for improving the understanding of chemical-biological effectiveness and will provide a theoretical basis for accurately assessing potential risks to ecosystems and human health. The combination with the quantitative structure-activity relationship model, Bayesian method, and machine learning technology are potential solutions to the previous research gaps.

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.

Impact of chemical pollution on threatened marine mammals: A systematic review

Marine mammals, due to their long life span, key position in the food web, and large lipid deposits, often face significant health risks from accumulating contaminants. This systematic review examines published literature on pollutant-induced adverse health effects in the International Union for Conservation of Nature (IUCN) red-listed marine mammal species. Thereby, identifying gaps in literature across different extinction risk categories, spatial distribution and climatic zones of studied habitats, commonly used methodologies, researched pollutants, and mechanisms from cellular to population levels. Our findings reveal a lower availability of exposure-effect data for higher extinction risk species (critically endangered 16%, endangered 15%, vulnerable 66%), highlighting the need for more research. For many threatened species in the Southern Hemisphere pollutant-effect relationships are not established. Non-destructively sampled tissues, like blood or skin, are commonly measured for exposure assessment. The most studied pollutants are POPs (31%), metals (30%), and pesticides (17%). Research on mixture toxicity is scarce while pollution-effect studies primarily focus on molecular and cellular levels. Bridging the gap between molecular data and higher-level effects is crucial, with computational approaches offering a high potential through in vitro to in vivo extrapolation using (toxico-)kinetic modelling. This could aid in population-level risk assessment for threatened marine mammals.

Additive Models Reveal Sources of Metals and Organic Pollutants in Norwegian Marine Sediments

We characterized spatial patterns of surface sediment concentrations of seven polychlorinated biphenyls (PCBs), seven polycyclic aromatic hydrocarbons (PAHs), three chlorinated pesticides, and five metals in Norwegian waters and Skagerrak. In total, we analyzed 5036 concentrations of 22 chemical substances that were measured between 1986 and 2014 at 333 sampling sites by means of generalized additive models (GAMs). We found that GAMs with organic carbon content of the sediment and latitude and longitude as co-variates explained as ca. 75% of the variability of the contaminant sediment concentrations. For metals, a predominantly hotspot-driven spatial pattern was found, i.e., we identified historical pollution hotspots (e.g., Sørfjord in western Norway) for mercury, zinc, cadmium, and lead. Highest concentrations of PAHs and PCBs were found close to densely populated and industrialized regions, i.e., in the North Sea and in the Kattegat and Skagerrak. The spatial pattern of the PCBs suggests the secondary and diffuse atmospheric nature of their sources. Atmospheric inputs are the main sources of pollution for most organic chemicals considered, but north of the Arctic circle, we found that concentrations of PAHs increased from south to north most likely related to a combination of coal-eroding bedrock and the biological pump. The knowledge acquired in the present research is essential for developing effective remediation strategies that are consistent with international conventions on pollution control.

Carbon, nitrogen and sulphur isotopic fractionation in captive juvenile hooded seal (Cystophora cristata): Application for diet analysis

RATIONALE

Intrinsic biogeochemical markers, such as stable isotope ratios of carbon, nitrogen and sulphur, are increasingly used to trace the trophic ecology of marine top predators. However, insufficient knowledge of fractionation processes in tissues continues to hamper the use of these markers.

METHODS

We performed a controlled feeding experiment with eight juvenile hooded seals (Cystophora cristata) that were held on a herring-based diet (Clupea harengus) for two years. Stable isotope ratios were measured via isotope ratio mass spectrometry in three of their tissues and related to values of these markers in their diet.

RESULTS

Diet-tissue isotope enrichment (trophic enrichment factor, TEF) values between dietary herring and seal tissues for carbon (Δ¹³ C) were +0.7 ‰ for red blood cells, +1.9 ‰ for hair and +1.1 ‰ for muscle. The TEFs for nitrogen trophic (Δ¹⁵ N) were +3.3 ‰ for red blood cells, +3.6 ‰ for hair and +4.3 ‰ for muscle. For sulphur, the Δ³⁴ S values were +1.1 ‰ for red blood cells, +1.0 ‰ for hair and +0.9 ‰ for muscle.

CONCLUSIONS

These enrichment values were greater than those previously measured in adult seals. This increase may be related to the higher rate of protein synthesis and catabolism in growing animals. This study is the first report on sulphur isotope enrichment values for a marine mammal species.

Accumulation of polybrominated diphenyl ethers in the brain compared with the levels in other tissues among different vertebrates from an e-waste recycling site

This study aimed to investigate the accumulation of polybrominated diphenyl ethers (PBDEs) in the brain compared with that in other tissues among different vertebrates. We collected mice, chickens, ducks, frogs, and fish from an e-waste recycling region in Taizhou, China, and measured PBDE concentrations in brain, liver and muscle tissues. The levels of PBDE in the tissues of mice, chickens, ducks, frogs and fish ranged 0.45-206, 0.06-18.8, 1.83-112, 2.75-108, and 0.02-32.0 ng/g wet weight, respectively. Preferential distribution in the liver and muscle relative to the brain was observed for PBDEs in mice, chickens, ducks and frogs. However, a high retention in the brain compared to the liver and muscle was observed in fish. Comparison of the brain/liver concentration (B/L) ratios revealed differences in PBDEs accumulation in the brain among these vertebrates. PBDEs accumulation in the brain was greatest in fish, followed by frogs, while the lowest accumulation occurred in the brains of mammals and birds. The findings apparently coincided with the evolution of the blood-brain barrier (BBB) across vertebrates, i.e. the BBB of fish might be less efficient than those of mammals, birds and amphibian. Low brominated congeners (such as BDE-28, BDE-47 and BDE-99) were predominant in the brains of investigated vertebrates, whereas BDE-209 was most abundant in liver and muscle tissues of mice, chickens and ducks. Significant differences in B/L ratios among PBDE congeners were found in both mice and chickens (p < 0.05). Particularly in mice, the B/L ratios of PBDE congeners presented a declining trend with increased bromine number. Our findings suggested that low brominated congeners might have a higher capacity to penetrate the BBB and accumulate in the brain, whereas high brominated congeners such as BDE-209 might have less potency to pass through the barrier. Further experimental studies are needed to confirm our findings.

Additive modelling reveals spatiotemporal PCBs trends in marine sediments

We developed generalised additive mixed models (GAMMs) to infer spatiotemporal trends of environmental PCB concentrations from an extensive dataset (n=1219) of PCB concentrations measured between 1991 and 2010 in sediments of the Belgian Coastal Zone (BCZ) and the Western Scheldt estuary. A GAMM with time, geographical zone, periodicity and the organic carbon - water partition coefficient as covariates explained 49% of the variability in the log transformed PCB sediment concentrations. The time trends unraveled two to threefold PCB concentration decreases in the BCZ during the last 20 years. However, in the Western Scheldt estuary, time trends were spatially heterogeneous and not significantly decreasing. These results demonstrate that international efforts to cut down emissions of PCBs have been effective to reduce concentrations in open water ecosystems like the BCZ but had little effect in the urbanised and industrialised area of the Scheldt estuary.

Lifetime PCB 153 bioaccumulation and pharmacokinetics in pilot whales: Bayesian population PBPK modeling and Markov chain Monte Carlo simulations

Physiologically based pharmacokinetic (PBPK) models for wild animal populations such as marine mammals typically have a high degree of model uncertainty and variability due to the scarcity of information and the embryonic nature of this field. Parameters values used in marine mammals models are usually taken from other mammalian species (e.g. rats or mice) and might not be entirely suitable to properly explain the kinetics of pollutants in marine mammals. Therefore, several parameters for a PBPK model for the bioaccumulation and pharmacokinetics of PCB 153 in long-finned pilot whales were estimated in the present study using the Bayesian approach executed with Markov chain Monte Carlo (MCMC) simulations. This method uses 'prior' information of the parameters, either from the literature or from previous model runs. The advantage is that this method uses such 'prior' parameters to calculate probability distributions to determine 'posterior' values that best explain the field observations. Those field observations or datasets were PCB 153 concentrations in blubber of long-finned pilot whales from Sandy Cape and Stanley, Tasmania, Australia. The model predictions showed an overall decrease in PCB 153 levels in blubber over the lifetime of the pilot whales. All parameters from the Sandy Cape model were updated using the Stanley dataset, except for the concentration of PCB 153 in the milk. The model presented here is a promising and preliminary start to PBPK modeling in long-finned pilot whales that would provide a basis for non-invasive studies in these protected marine mammals.

Modeling PCB-bioaccumulation in the Bottlenose Dolphin ( Tursiops truncatus ): estimating a dietary threshold concentration

An individually based (IB) model to predict PCB concentrations in the bottlenose dolphin population of Charleston, SC, USA, was developed with the aim to gain a better understanding of the bioaccumulation behavior and health risk of dietary PCBs across the population and their prey. PCB concentrations predicted in male and female bottlenose dolphin were in good agreement with observed tissue concentrations corroborating the reliability of the model performance and its utility in gaining a more complete view of risk. The modeled cumulative distribution of ΣPCB concentrations for the population with a breakdown into juvenile, adult male, and female subclasses ranged from 3600 to 144,400 ng/g lipid with 66% to >80% of the population exceeding the established threshold for adverse health effects of 17,000 ng/g lipid. The model estimated that a dietary PCB concentration not exceeding 5.1 ng/g wet wt would be required to reach a condition where 95% of the population would have tissue levels below the health effect threshold. The IB model for PCBs in bottlenose dolphins provides a novel approach to estimating the maximum acceptable dietary concentration for PCBs, a central and important factor to protect these apex predators. The model also enables effective prediction of concentrations in dolphins from fish contaminant surveys which are logistically easier and less costly to collect.

Relationships between in vitro lymphoproliferative responses and levels of contaminants in blood of free-ranging adult harbour seals (Phoca vitulina) from the North Sea

In vitro culture of peripheral blood leucocytes (PBLs) is currently used in toxicological studies of marine mammals. However, blood cells of wild individuals are exposed in vivo to environmental contaminants before being isolated and exposed to contaminants in vitro. The aim of this study was to highlight potential relationships between blood contaminant levels and in vitro peripheral blood lymphocyte proliferation in free-ranging adult harbour seals (Phoca vitulina) from the North Sea. Blood samples of 18 individuals were analyzed for trace elements (Fe, Zn, Se, Cu, Hg, Pb, Cd) and persistent organic contaminants and metabolites (ΣPCBs, ΣHO-PCBs, ΣPBDEs, 2-MeO-BDE68 and 6-MeO-BDE47, ΣDDXs, hexachlorobenzene, oxychlordane, trans-nonachlor, pentachlorophenol and tribromoanisole). The same samples were used to determine the haematology profiles, cell numbers and viability, as well as the in vitro ConA-induced lymphocyte proliferation expressed as a stimulation index (SI). Correlation tests (Bravais-Pearson) and Principal Component Analysis with multiple regression revealed no statistically significant relationship between the lymphocyte SI and the contaminants studied. However, the number of lymphocytes per millilitre of whole blood appeared to be negatively correlated to pentachlorophenol (r=-0.63, p=0.005). In adult harbour seals, the interindividual variations of in vitro lymphocyte proliferation did not appear to be directly linked to pollutant levels present in the blood, and it is likely that other factors such as age, life history, or physiological parameters have an influence. In a general manner, experiments with in vitro immune cell cultures of wild marine mammals should be designed so as to minimize confounding factors in which case they remain a valuable tool to study pollutant effects in vitro.

Application of Bayesian population physiologically based pharmacokinetic (PBPK) modeling and Markov chain Monte Carlo simulations to pesticide kinetics studies in protected marine mammals: DDT, DDE, and DDD in harbor porpoises

Physiologically based pharmacokinetic (PBPK) modeling in marine mammals is a challenge because of the lack of parameter information and the ban on exposure experiments. To minimize uncertainty and variability, parameter estimation methods are required for the development of reliable PBPK models. The present study is the first to develop PBPK models for the lifetime bioaccumulation of p,p'-DDT, p,p'-DDE, and p,p'-DDD in harbor porpoises. In addition, this study is also the first to apply the Bayesian approach executed with Markov chain Monte Carlo simulations using two data sets of harbor porpoises from the Black and North Seas. Parameters from the literature were used as priors for the first "model update" using the Black Sea data set, the resulting posterior parameters were then used as priors for the second "model update" using the North Sea data set. As such, PBPK models with parameters specific for harbor porpoises could be strengthened with more robust probability distributions. As the science and biomonitoring effort progress in this area, more data sets will become available to further strengthen and update the parameters in the PBPK models for harbor porpoises as a species anywhere in the world. Further, such an approach could very well be extended to other protected marine mammals.

Computational toxicology: Physiologically based pharmacokinetic models (PBPK) for lifetime exposure and bioaccumulation of polybrominated diphenyl ethers (PBDEs) in marine mammals

Due to migration of harbour porpoises towards more polluted areas like the North Sea and their sensitivity towards pollution, there is a need for proper conservation measures for this species. As a consequence, knowledge about the pollutant's kinetics is required. The present study is the first to investigate the kinetics of PBDEs in marine mammals using PBPK modeling as a non-destructive tool for describing the chemical's kinetics in a protected animal species. The models were developed and parameterized using data from the literature and Black Sea harbour porpoises through computer optimization. The predictability of these models in time was assessed by reverse dosimetry modeling using data from North Sea porpoises (1990-2008). From these predictions, PBDE 99 levels were found to decrease the fastest, followed by PBDE 153, 47 and 100. Results show that the PBPK models can be applied for harbour porpoises from different regions and also simulate time trends.