Assessing the bioaccumulation potential of ionizable organic compounds: Current knowledge and research priorities

Accumulation and trophic transfer of per- and polyfluoroalkyl substances (PFAS) in estuarine organisms determined via stable isotopes

Per- and polyfluoroalkyl substances (PFAS) are persistent organic pollutants in estuaries. In this study, 19 PFAS were quantified in surface waters, sediments, marine invertebrates (aquatic worms, Eastern oysters, and blue crab), and forage fish (Atlantic silverside, four-spine stickleback, mummichog, sheepshead minnow, and rainwater killifish) in an aqueous film forming foam (AFFF)-contaminated estuary, Georgica Pond (NY, USA). Carbon and nitrogen stable isotopes (δ13C and δ15N) were used to determine trophic position of organisms and to identify modes of PFAS exposure. The influence of salinity (8 to 26 practical salinity units, PSU) on the relative and absolute abundance of PFAS in all matrices was also investigated. Eleven long- and short-chain perfluoroalkyl acids (PFAAs) were found to have bioaccumulation potential (bioaccumulation factor, BAF; biota-sediment accumulation factor, BSAF) and were positively correlated with relative trophic position. Among these, long-chain PFAAs (perfluorohexanesulfonic acid, PFHxS; perfluorooctane sulfonic acid, PFOS; perfluorooctanoic acid, PFOA; perfluorononanoic acid, PFNA) were the greatest contributors to total body burden and bioaccumulated in all organisms, with PFOS (log BAF = 3.55 ± 0.83) and PFNA (log BAF = 3.17 ± 0.46) having the highest mean values of all compounds. PFOS was present in all biota samples and concentrations significantly increased with food web trophic position (ranging from 0.18 to 777 μg kg-1). Perfluorobutane sulfonic acid (PFBS) was also ubiquitous among all organisms, bioaccumulating in both invertebrate and vertebrate species. Total PFAS concentrations in aquatic worms were significantly higher in lower salinity water while the PFAS profile of Eastern oysters shifted from predominately perfluorocarboxylic acids (66 % of total composition) to perfluorosulfonic acids (62 %) as the ecosystem transitioned from low (9 PSU) to high (25 PSU) salinity. Collectively, this study demonstrates the utility of applying δ13C and δ15N to determine bioaccumulation patterns of both legacy PFAS and short-chain replacement compounds and underscores how shifts in salinity can alter the concentration and speciation of PFAS in estuaries.

Construction of interpretable ensemble learning models for predicting bioaccumulation parameters of organic chemicals in fish

Accurate prediction of bioaccumulation parameters is essential for assessing exposure, hazards, and risks of chemicals. However, the majority of prediction models on bioaccumulation parameters are individual models based on a single algorithm and lack model interpretation, resulting in unsatisfactory prediction accuracy due to inherent constraints of the algorithm and weak interpretability. Ensemble learning (EL) that combine multiple algorithms, coupled with SHapley Additive exPlanation (SHAP) method, may overcome the limitations. Herein, EL models were constructed for three bioaccumulation parameters using datasets covering 2496 chemicals. The EL models demonstrated superior prediction accuracy compared to both individual models developed in this study and those from previous research, achieving a coefficient of determination of up to 0.861 on the validation sets. Applicability domains were characterized using a structure-activity landscape-based (abbreviated as ADSAL) methodology. The optimal EL models, together with the ADSAL, were successfully used to predict bioaccumulation parameters for 4374 chemicals included in the Inventory of Existing Chemical Substances of China. Model interpretation using the SHAP method offered insight into key features influencing bioaccumulation potential, including hydrophobicity, water solubility, polarizability, ionization potential, weight, and volume of molecules. Overall, the study provides data and models to support the sound management and risk assessment of chemicals.

Global Assessment of Emerging Contaminant Removal in Wastewater Treatment Plants: In Silico Hazard Screening and Risk Evaluation

Pharmaceuticals and personal care products (PPCPs) are emerging contaminants (ECs), whose presence in the environment is of increasing concern due to their widespread use and possible detrimental effects on wildlife and humans. These chemicals may present multiple hazardous properties such as environmental persistence, toxicity, high mobility, and the potential for bioaccumulation. In this study, extended bibliographic research was conducted to characterize the removal efficiency (RE) of PPCPs in wastewater treatment plants (WWTPs) considering different technologies. Measured values of RE were collected from the literature or calculated for 251 compounds. The molecular structure of the 245 PPCPs were used as the input to generate predictions of multiple properties using several QSAR tools, such as the OECD Toolbox, OPERA, EPI Suite™, and QSAR-ME Profiler. These predictions were compared to regulatory thresholds to identify hazardous chemicals and to screen persistent, mobile and toxic (PMT) or persistent, bioaccumulative and toxic (PBT) substances. Finally, chemicals were prioritized by combining values of RE and QSAR predictions for multiple properties. A total of 16 out of the 245 molecules were prioritized as the most hazardous compounds to the aquatic environment and, among these, six were associated with potential risk due to their exposure concentrations reported in the literature.

Toxicity of a management bait for grass carp (Ctenopharyngodon idella) incorporated with Antimycin A

No current technology can specifically target grass carp (Ctenopharyngodon idella) for control within aquatic ecosystems. Rotenone and Carbon Dioxide-Carp are currently the only available registered pesticides for grass carp; they are nonselective and typically applied throughout the water, equally exposing target and native species. A more selective control tool or pesticide application could be used by resource managers to support mitigation efforts. Development of delivery systems that exploit carp feeding strategies could increase selectivity of pesticides and minimize effects on native fishes. A pesticide with selective delivery could be less labor intensive and used within an integrative pest management strategy. The present study examined Antimycin A toxicity in juvenile and sub-adult grass carp and rainbow trout (Oncorhynchus mykiss) across two routes of exposure. Water-based toxicity studies were used to calculate the concentration to cause lethality in 50% of treated fish (LC50) at 24-h, while oral gavage toxicity studies were used to calculate the dose to cause lethality in 50% of treated grass carp and rainbow trout (LD50) 24- to 96-h. Although rainbow trout were more sensitive than grass carp to Antimycin A through water-based exposure, oral toxicity was similar between species, even with inherent gastrointestinal morphological differences. Successful delivery of a lethal dose of Antimycin A to grass carp was achieved through an oral route of exposure using the rapeseed bait and shows promise for registration as a control tool and eventual use in pest management plans. Although a lethal dose of Antimycin A could be incorporated into a single bait pellet, more bait was required to achieve desired mortality when fed to fish under laboratory conditions.

Differential tissue distribution of pharmaceuticals in a wild subtropical marine fish

To date, the presence of pharmaceuticals has been extensively documented across a wide range of aquatic systems and biota. Further, substantial progress has been made in transitioning from laboratory assessments of pharmaceutical fate and effects in fish to in situ assessments of exposure and effects; however, certain research areas remain understudied. Among these is investigation of differential accumulation across multiple internal tissues in wild marine fish beyond the species commonly sampled in laboratory and freshwater field settings. This study examined the presence of pharmaceuticals across four tissues (plasma, muscle, brain, and liver) in a wild marine fish, bonefish (Albula vulpes), throughout coastal South Florida, USA. Differential accumulation across tissues was assessed for the number and concentration, identity, and composition of accumulated pharmaceuticals by sampling 25 bonefish and analyzing them for 91 pharmaceuticals. The concentration of pharmaceuticals was highest in plasma > liver > brain > muscle, while the number of pharmaceuticals was highest in liver > brain > plasma > muscle. The identity of detected pharmaceuticals was tissue specific, and there was an inverse relationship between the number of detections for each pharmaceutical and its log Kow. The composition of pharmaceuticals was tissue specific for both pharmaceutical presence/absence and concentration. Across all tissues, the greatest similarity was between brain and liver, which were more similar to plasma than to muscle, and muscle was the most distinct tissue. For tissue compositional variability, muscle was the most diverse in accumulated pharmaceuticals, while plasma, brain, and liver were similarly variable. With the highest concentrations in plasma and highest number in liver, and documented variability in accumulated pharmaceuticals across tissues, our results highlight the importance of tissue selection when surveying exposure in wild fish, suggesting that multi-tissue analysis would allow for a more comprehensive assessment of exposure diversity and risk of adverse effects.

Identifying pathways of pharmaceutical exposure in a mesoconsumer marine fish

Pharmaceutical uptake involves processes that vary across aquatic systems and biota. However, single studies examining multiple environmental compartments, microhabitats, biota, and exposure pathways in mesoconsumer fish are sparse. We investigated the pharmaceutical burden in bonefish (Albula vulpes), pathways of exposure, and estimated exposure to a human daily dose. To evaluate exposure pathways, the number and composition of pharmaceuticals across compartments and the bioconcentration in prey and bonefish were assessed. To evaluate bioaccumulation, we proposed the use of a field-derived bioaccumulation factor (fBAF), due to variability inherent to natural systems. Exposure to a human daily dose was based on bonefish daily energetic requirements and consumption rates using pharmaceutical concentrations in prey. Pharmaceutical number and concentration were highest in prey, followed by bonefish, water and sediment. Fifteen pharmaceuticals were detected in common among bonefish, prey, and water; all of which bioconcentrated in prey and bonefish, and four bioaccumulated in bonefish. The composition of detected pharmaceuticals was compartment specific, and prey were most similar to bonefish. Bonefish were exposed to a maximum of 1.2 % of a human daily dose via prey consumption. Results highlight the need for multicompartment assessments of exposure and consideration of prey along with water as a pathway of exposure.

Integrated Transfer Learning and Multitask Learning Strategies to Construct Graph Neural Network Models for Predicting Bioaccumulation Parameters of Chemicals

Accurate prediction of parameters related to the environmental exposure of chemicals is crucial for the sound management of chemicals. However, the lack of large data sets for training models may result in poor prediction accuracy and robustness. Herein, integrated transfer learning (TL) and multitask learning (MTL) was proposed for constructing a graph neural network (GNN) model (abbreviated as TL-MTL-GNN model) using n-octanol/water partition coefficients as a source domain. The TL-MTL-GNN model was trained to predict three bioaccumulation parameters based on enlarged data sets that cover 2496 compounds with at least one bioaccumulation parameter. Results show that the TL-MTL-GNN model outperformed single-task GNN models with and without the TL, as well as conventional machine learning models trained with molecular descriptors or fingerprints. Applicability domains were characterized by a state-of-the-art structure-activity landscape-based (abbreviated as ADSAL) methodology. The TL-MTL-GNN model coupled with the optimal ADSAL was employed to predict bioaccumulation parameters for around 60,000 chemicals, with more than 13,000 compounds identified as bioaccumulative chemicals. The high predictive accuracy and robustness of the TL-MTL-GNN model demonstrate the feasibility of integrating the TL and MTL strategy in modeling small-sized data sets. The strategy holds significant potential for addressing small data challenges in modeling environmental chemicals.

Comparative in vitro hepatic clearances of commonly used antidepressants, antipsychotics, and anti-inflammatory agents in rainbow trout liver S9 fractions

Residues of human pharmaceuticals are widely detected in surface waters and can be taken up by and bioaccumulate in aquatic organisms, especially fish. One of the key challenges in assessing the bioaccumulation potential of ionizable organic compounds, such as the pharmaceuticals, is the lack of empirical data for biotransformation. In the present study, we assessed the in vitro intrinsic clearances (CLINT) of twelve pharmaceuticals, individually and some additionally as mixtures, in rainbow trout (Oncorhynchus mykiss) liver S9 fractions (RT-S9) adhering to the OECD test guidance 319B. The test substances included four anti-inflammatory agents (diclofenac, ibuprofen, ketoprofen, naproxen), seven antidepressants/antipsychotics (citalopram, haloperidol, levomepromazine, mirtazapine, risperidone, sertraline, venlafaxine) and the O-desmethyl metabolite of venlafaxine. Quantifiable intrinsic clearances were detected for diclofenac, ibuprofen, naproxen, levomepromazine, and sertraline. Apart from diclofenac, the in vitro clearances of the other four pharmaceuticals were shown to be critically dependent on the cytochrome P450 (CYP) metabolism. Therefore, we also determined the half-maximal inhibitory concentrations (IC50) of the same twelve pharmaceuticals toward CYP1A-like (7-ethoxyresorufin-O-deethylation, EROD) and CYP3A-like (benzyloxy-4-trifluoromethylcoumarin-O-debenzyloxylation, BFCOD) activities in RT-S9 using IC50 shift assay. As a result, levomepromazine and sertraline were identified as the most potent inhibitors of both EROD and BFCOD activity (unbound IC50 < 10 µM each), followed by citalopram and haloperidol (10 µM < IC50 < 100 µM). Additionally, mirtazapine was a selective EROD inhibitor (IC50 ∼ 30 µM). The inhibitory impacts of haloperidol and sertraline were indicatively time dependent. Finally, we carried out intrinsic clearance assays with mixtures of diclofenac, ibuprofen, naproxen, levomepromazine, and sertraline to examine the impacts of EROD and BFCOD inhibitions on their in vitro CLINT in RT-S9. Our in vitro data suggests that the intrinsic clearances of ibuprofen, levomepromazine, and sertraline in rainbow trout can be significantly reduced as the result of P450 inhibition by pharmaceutical mixtures, whereas the clearances of diclofenac and naproxen are less impacted.

Evaluation of the EMA log kow trigger for fish BCF testing based on data for several human pharmaceuticals

In the European Medicines Agency (EMA) "Guideline for Environmental Risk Assessment of Medicinal Products for Human Use," a fish bioconcentration factor (BCF) study is triggered in Phase I for pharmaceuticals having log Kow >4.5, to support Persistence, Bioaccumulation and Toxicity (PBT) screening, and in Phase II to assess secondary poisoning and bioaccumulation ('B') potential when log Kow ≥3. The standard sampling schedule outlined in OECD Test Guideline 305 (TG305) may require assessment of approximately 200 fish following exposure to low- and high-test concentrations and a negative control. We report experimental log Kow and BCF values for 64 human pharmaceuticals that were used to evaluate the current BCF testing trigger of log Kow ≥3, and whether a single BCF exposure concentration allows accurate classification of bioaccumulation potential. Our data support raising the BCF testing trigger to log Kow ≥4, and use of a single test concentration. The resulting reduction in the use of fish is consistent with the 3 R s principle and did not adversely affect classification accuracy. An assessment of potential risk of secondary poisoning was also conducted for three drugs classified as either B or vB, and no risks were identified.

Trophic transfer of antibiotics in the benthic-pelagic coupling foodweb in a macrophyte-dominated shallow lake: The importance of pelagic-benthic coupling strength and baseline organism

In lake ecosystems, pelagic-benthic coupling strength (PBCS) is closely related to foodweb structure and pollutant transport. However, the trophic transfer of antibiotics in a benthic-pelagic coupling foodweb (BPCFW) and the manner in which PBCS influences the trophic magnification factor (TMFs) of antibiotics is still not well understood in the whole lake. Herein, the trophic transfer behavior of 12 quinolone antibiotics (QNs) in the BPCFW of Baiyangdian Lake were studied during the period of 2018-2019. It was revealed that 24 dominant species were contained in the BPCFW, and the trophic level was 0.42-2.94. Seven QNs were detected in organisms, the detection frequencies of ofloxacin (OFL), flumequine (FLU), norfloxacin (NOR), and enrofloxacin (ENR) were higher than other QNs. The ∑QN concentration in all species was 11.3-321 ng/g dw. The TMFs for ENR and NOR were trophic magnification, while for FLU/OFL it was trophic dilution. The PBCS showed spatial-temporal variation, with a range of 0.6977-0.7910. The TMFs of ENR, FLU, and OFL were significantly positively correlated with PBCS. Phytoplankton and macrophyte biomasses showed indirect impact on the TMFs of QNs by directly influencing the PBCS. Therefore, the PBCS was the direct influencing factor for the TMFs of chemicals.

Pharmaceutical Residues in Edible Oysters along the Coasts of the East and South China Seas and Associated Health Risks to Humans and Wildlife

The investigation of pharmaceuticals as emerging contaminants in marine biota has been insufficient. In this study, we examined the presence of 51 pharmaceuticals in edible oysters along the coasts of the East and South China Seas. Only nine pharmaceuticals were detected. The mean concentrations of all measured pharmaceuticals in oysters per site ranged from 0.804 to 15.1 ng g-1 of dry weight, with antihistamines being the most common. Brompheniramine and promethazine were identified in biota samples for the first time. Although no significant health risks to humans were identified through consumption of oysters, 100-1000 times higher health risks were observed for wildlife like water birds, seasnails, and starfishes. Specifically, sea snails that primarily feed on oysters were found to be at risk of exposure to ciprofloxacin, brompheniramine, and promethazine. These high risks could be attributed to the monotonous diet habits and relatively limited food sources of these organisms. Furthermore, taking chirality into consideration, chlorpheniramine in the oysters was enriched by the S-enantiomer, with a relative potency 1.1-1.3 times higher when chlorpheniramine was considered as a racemate. Overall, this study highlights the prevalence of antihistamines in seafood and underscores the importance of studying enantioselectivities of pharmaceuticals in health risk assessments.

What Approaches Should be Used to Prioritize Pharmaceuticals and Personal Care Products for Research on Environmental and Human Health Exposure and Effects?

Pharmaceuticals and personal care products (PPCPs) are released from multiple anthropogenic sources and thus have a ubiquitous presence in the environment. The environmental exposure and potential effects of PPCPs on biota and humans has aroused concern within the scientific community and the public. Risk assessments are commonly conducted to evaluate the likelihood of chemicals including PPCPs that pose health threats to organisms inhabiting various environmental compartments and humans. Because thousands of PPCPs are currently used, it is impractical to assess the environmental risk of all of them due to data limitations; in addition, new PPCPs are continually being produced. Prioritization approaches, based either on exposure, hazard, or risk, provide a possible means by which those PPCPs that are likely to pose the greatest risk to the environment are identified, thereby enabling more effective allocation of resources in environmental monitoring programs in specific geographical locations and ecotoxicological investigations. In the present review, the importance and current knowledge concerning PPCP occurrence and risk are discussed and priorities for future research are proposed, in terms of PPCP exposure (e.g., optimization of exposure modeling in freshwater ecosystems and more monitoring of PPCPs in the marine environment) or hazard (e.g., differential risk of PPCPs to lower vs. higher trophic level species and risks to human health). Recommended research questions for the next 10 years are also provided, which can be answered by future studies on prioritization of PPCPs. Environ Toxicol Chem 2024;43:488-501. © 2022 SETAC.

Predicting the Accumulation of Ionizable Pharmaceuticals and Personal Care Products in Aquatic and Terrestrial Organisms

The extent to which chemicals bioaccumulate in aquatic and terrestrial organisms represents a fundamental consideration for chemicals management efforts intended to protect public health and the environment from pollution and waste. Many chemicals, including most pharmaceuticals and personal care products (PPCPs), are ionizable across environmentally relevant pH gradients, which can affect their fate in aquatic and terrestrial systems. Existing mathematical models describe the accumulation of neutral organic chemicals and weak acids and bases in both fish and plants. Further model development is hampered, however, by a lack of mechanistic insights for PPCPs that are predominantly or permanently ionized. Targeted experiments across environmentally realistic conditions are needed to address the following questions: (1) What are the partitioning and sorption behaviors of strongly ionizing chemicals among species? (2) How does membrane permeability of ions influence bioaccumulation of PPCPs? (3) To what extent are salts and associated complexes with PPCPs influencing bioaccumulation? (4) How do biotransformation and other elimination processes vary within and among species? (5) Are bioaccumulation modeling efforts currently focused on chemicals and species with key data gaps and risk profiles? Answering these questions promises to address key sources of uncertainty for bioaccumulation modeling of ionizable PPCPs and related contaminants. Environ Toxicol Chem 2024;43:502-512. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Chemical Space Covered by Applicability Domains of Quantitative Structure-Property Relationships and Semiempirical Relationships in Chemical Assessments

A significant number of chemicals registered in national and regional chemical inventories require assessments of their potential "hazard" concerns posed to humans and ecological receptors. This warrants knowledge of their partitioning and reactivity properties, which are often predicted by quantitative structure-property relationships (QSPRs) and other semiempirical relationships. It is imperative to evaluate the applicability domain (AD) of these tools to ensure their suitability for assessment purpose. Here, we investigate the extent to which the ADs of commonly used QSPRs and semiempirical relationships cover seven partitioning and reactivity properties of a chemical "space" comprising 81,000+ organic chemicals registered in regulatory and academic chemical inventories. Our findings show that around or more than half of the chemicals studied are covered by at least one of the commonly used QSPRs. The investigated QSPRs demonstrate adequate AD coverage for organochlorides and organobromines but limited AD coverage for chemicals containing fluorine and phosphorus. These QSPRs exhibit limited AD coverage for atmospheric reactivity, biodegradation, and octanol-air partitioning, particularly for ionizable organic chemicals compared to nonionizable ones, challenging assessments of environmental persistence, bioaccumulation capability, and long-range transport potential. We also find that a predictive tool's AD coverage of chemicals depends on how the AD is defined, for example, by the distance of a predicted chemical from the centroid of the training chemicals or by the presence or absence of structural features.

Biomagnification and health risks of perflfluoroalkyl acids (PFAAs) in seafood from the Yangtze river estuary of China

Bioaccumulation and human health risk assessment of Perfluoroalkyl acids (PFAAs) is important for pollutant hazard assessment. In this study, 26 aquatic organisms were collected from the Yangtze River estuary, the PFAAs concentrations in organisms were detected by liquid chromatography-mass spectrometry, and the trophic levels of organisms were constructed using nitrogen isotope analysis. The results showed that Perfluorobutane sulfonate (PFBS) was predominant in organisms with the mean concentration of 6.43 ± 8.21 ng/g ww. The biomagnification of organisms along the food chain was widespread, and the biomagnification factor (BMF) of perfluorooctane sulfonic (PFOS) was the most prominent. Trophic magnifcation factors (TMFs) of PFAAs were estimated in the marine food web, and TMFs >1 were observed in Perfluorodecanoic acid (PFDA), Perfluoroundecanoic acid (PFUnDA), Perfluorododecanoic acid (PFDoDA), and PFOS, indicating the biomagnifcation effects of these 4 individual PFAAs in organisms at Yangtze River estuary. The estimated daily intake (EDI) of PFBS was highest in adolescents aged 6-18 years, with EDIs of 18.9 ng/kg·bw/day for males and 14.0 ng/kg·bw/day for females. The hazard ratio (HR) of PFAAs reported in different age and gender groups were lower than 1.

Occurrence, Distribution, and Trophic Transfer of Pharmaceuticals and Personal Care Products in the Bohai Sea

The ubiquitous presence of pharmaceuticals and personal care products (PPCPs) in environments has aroused global concerns; however, minimal information is available regarding their multimedia distribution, bioaccumulation, and trophic transfer in marine environments. Herein, we analyzed 77 representative PPCPs in samples of surface and bottom seawater, surface sediments, and benthic biota from the Bohai Sea. PPCPs were pervasively detected in seawater, sediments, and benthic biota, with antioxidants being the most abundant PPCPs. PPCP concentrations positively correlated between the surface and bottom water with a decreasing trend from the coast to the central oceans. Higher PPCP concentrations in sediment were found in the Yellow River estuary, and the variations in the physicochemical properties of PPCPs and sediment produced a different distribution pattern of PPCPs in sediment from seawater. The log Dow, but not log Kow, showed a linear and positive relationship with bioaccumulation and trophic magnification factors and a parabolic relationship with biota-sediment accumulation factors. The trophodynamics of miconazole and acetophenone are reported for the first time. This study provides novel insights into the multimedia distribution and biomagnification potential of PPCPs and suggests that log Dow is a better indicator of their bioaccumulation and trophic magnification.

Assessment of textile industry effluent (untreated and microbially treated) induced genotoxic, haematological, biochemical, histopathological and ultrastructural alterations in fresh water fish Channa punctata

The unregulated expulsion of untreated textile water into water bodies is a major hazard to aquatic ecosystems. The present investigation was contrived to estimate the impact of textile dye bath effluent (untreated and microbially treated) on fish Channa punctata. Untreated effluent-exposed fish showed extremely altered behaviour (air gulping, erratic and speedy movements, increased opercular activity) and morphology (deposition of dyes on skin and scales, high pigmentation, mucus exudation). Significantly increased micronuclei (1.61-, 1.28-, 1.38-fold) and aberrant cell frequency (1.37-, 1.45-, 1.28-fold) was observed in untreated group as compared to treated group after 15, 30, and 45 days of exposure. Tail length, % tail intensity, tail moment and olive tail moment were also enhanced in all the exposed tissues. However, maximum damage was noticed in gill tissues showing 1.19-, 1.37-, 1.34- and 1.50-fold increased TL, %TI, TM and OTM in untreated group as compared to treated group after 45 days of exposure. On comparing untreated and treated groups, increased blood parameters and significantly reduced white blood cell count (WBC) were noticed in treated group. Significantly enhanced alterations in biochemical parameters were also analysed in untreated group. Reduced alterations in enzymological levels of fishes exposed to treated effluent indicate lesser toxic nature of the degraded metabolites of dye. Histological analysis in fishes exposed to untreated effluent showed several deformities in liver (necrosis, congestion, fusion of cells and melanomacrophage infiltration) and gill tissues (necrosis, bending of lamellae and severe aneurysm). Scanning electron microscopy (SEM) analysis further reaffirmed the pathologies observed in histological analysis. Fewer structural alterations were noticed in treated effluent fishes. The results concluded that untreated effluent inflicted toxicity potential on morphology as well as physiological defects in fish, and the severity increased with increasing duration of exposure, whereas reduction in toxicity in microbially treated groups can be analysed for aquacultural purposes owing to their lesser toxic nature.

Research Priorities for the Environmental Risk Assessment of Per- and Polyfluorinated Substances

Per- and polyfluorinated substances (PFAS) are a group of thousands of ubiquitously applied persistent industrial chemicals. The field of PFAS environmental research is developing rapidly, but suffers from substantial biases toward specific compounds, environmental compartments, and organisms. The aim of our study was therefore to highlight current developments and to identify knowledge gaps and subsequent research needs that would contribute to a comprehensive environmental risk assessment for PFAS. To this end, we consulted the open literature and databases and found that knowledge of the environmental fate of PFAS is based on the analysis of <1% of the compounds categorized as PFAS. Moreover, soils and suspended particulate matter remain largely understudied. The bioavailability, bioaccumulation, and food web transfer studies of PFAS also focus on a very limited number of compounds and are biased toward aquatic biota, predominantly fish, and less frequently aquatic invertebrates and macrophytes. The available ecotoxicity data revealed that only a few PFAS have been well studied for their environmental hazards, and that PFAS ecotoxicity data are also strongly biased toward aquatic organisms. Ecotoxicity studies in the terrestrial environment are needed, as well as chronic, multigenerational, and community ecotoxicity research, in light of the persistency and bioaccumulation of PFAS. Finally, we identified an urgent need to unravel the relationships among sorption, bioaccumulation, and ecotoxicity on the one hand and molecular descriptors of PFAS chemical structures and physicochemical properties on the other, to allow predictions of exposure, bioaccumulation, and toxicity. Environ Toxicol Chem 2023;42:2302-2316. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Insights into the interaction of human serum albumin with ionic liquids - Thermodynamic, spectroscopic and molecular modelling studies

Human serum albumin (HSA) effectively binds different types of low-molecular-weight compounds and thus enables their distribution in living organisms. Recently, it has been reported that the protein-ligand interactions play a crucial role in bioaccumulation processes and provide an important sorption phase, especially for ionogenic compounds. Therefore, the binding interactions of such compounds with proteins are the subject of an ongoing interest in environmental and life sciences. In this paper, the influence of some counter-ions, namely [B(CN)4]- and [C(CN)3]- on the affinity of the [IM1-12]+ towards HSA has been investigated and discussed based on experimental methods (isothermal titration calorimetry and steady-state fluorescence spectroscopy) and molecular dynamics-based computational approaches. Furthermore, the thermal stability of the resulting HSA/ligand complexes was assessed using DSC and CD spectroscopy. As an outcome of the work, it has been ascertained that the protein is able to bind simultaneously the ligands under study but in different regions of HSA. Thus, the presence in the system of [IM1-12]+ does not disturb the binding of [C(CN)3]- and [B(CN)4]-. The presented results provide important information on the presence of globular proteins and some ionogenic compounds in the distribution and bioaccumulation of ILs in the environment and living organisms.

Phytoplankton Biological Pump Controls the Spatiotemporal Bioaccumulation and Trophic Transfer of Antibiotics in a Large Subtropical River

The spatiotemporal bioaccumulation, trophic transfer of antibiotics, and regulation of the phytoplankton biological pump were quantitatively evaluated in the Pearl River, South China. The occurrence of antibiotics in organisms indicated a significant spatiotemporal trend associated with the life cycle of phytoplankton. Higher temporal bioaccumulation factors (BAFs) were found in phytoplankton at the bloom site, while lower BAFs of antibiotics in organisms could not be explained by phytoplankton biomass dilution but were attributed to the low bioavailability of antibiotics, which was highly associated with distribution coefficients (R2 = 0.480-0.595, p < 0.05). Such lower BAFs of antibiotics in phytoplankton at higher biomass sites hampered the entry of antibiotics into food webs, and trophic dilutions were subsequently observed for antibiotics except for ciprofloxacin (CFX) and sulfamerazine (SMZ) at sites with blooms in all seasons. Distribution of CFX, norfloxacin (NFX), and sulfapyridine (SPD) showed further significant positive relationships with the plasma protein fraction (R2 = 0.275-0.216, p < 0.05). Both mean BAFs and trophic magnification factors (TMFs) were significantly negatively correlated with phytoplankton biomass (R2 = 0.661-0.741, p < 0.05). This study highlights the importance of the biological pump in the regulation of spatiotemporal variations in bioaccumulation and trophic transfer of antibiotics in anthropogenic-impacted eutrophic rivers in subtropical regions.

A weight of evidence approach for bioaccumulation assessment

Bioaccumulation assessments conducted by regulatory agencies worldwide use a variety of methods, types of data, metrics, and categorization criteria. Lines of evidence (LoE) for bioaccumulation assessment can include bioaccumulation metrics such as in vivo bioconcentration factor (BCF) and biomagnification factor (BMF) data measured from standardized laboratory experiments, and field (monitoring) data such as BMFs, bioaccumulation factors (BAFs), and trophic magnification factors (TMFs). In silico predictions from mass-balance models and quantitative structure-activity relationships (QSARs) and a combination of in vitro biotransformation rates and in vitro-in vivo extrapolation (IVIVE) models can also be used. The myriad bioaccumulation metrics and categorization criteria and underlying uncertainty in measured or modeled data can make decision-making challenging. A weight of evidence (WoE) approach is recommended to address uncertainty. The Bioaccumulation Assessment Tool (BAT) guides a user through the process of collecting and generating various LoE required for assessing the bioaccumulation of neutral and ionizable organic chemicals in aquatic (water-respiring) and air-breathing organisms. The BAT includes data evaluation templates (DETs) to critically evaluate the reliability of the LoE used in the assessment. The DETs were developed from standardized testing guidance. The approach used in the BAT is consistent with OECD and SETAC WoE principles and facilitates the implementation of chemical policy objectives in chemical assessment and management. The recommended methods are also iterative and tiered, providing pragmatic methods to reduce unnecessary animal testing. General concepts of the BAT are presented and case study applications of the tool for hexachlorobenzene (HCB) and β-hexachlorocyclohexane (β-HCH) are demonstrated. The BAT provides a consistent and transparent WoE framework to address uncertainty in bioaccumulation assessment and is envisaged to evolve with scientific and regulatory developments. Integr Environ Assess Manag 2023;19:1235-1253. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Prioritization of risks posed by synthetic chemicals manufactured in China toward humans and the environment via persistence, bioaccumulation, mobility and toxicity properties

Over a third of the global chemical production and sales occurred in China, which make effective assessment and management for chemicals produced by China's chemical industry essential not just for China but for the world. Here, we systematical assessed the persistence (P), bioaccumulation (B), mobility (M) and toxicity (T) potency properties for the chemicals listed in Inventory of Existing Chemical Substances of China (IECSC) via experimental data retrieved from large scale databases and in silico data generated with well-established models. Potential PBT, PMT and PB&MT substances were identified. High risk potentials were highlighted for groups of synthetic intermediates, raw materials, as well as a series of biocides. The potential PBT and PMT synthetic intermediates and/or raw materials unique to the IECSC were dominated with organofluorines, for example, the intermediates used as electronic light-emitting materials. Meanwhile, the biocides unique to the IECSC were mainly organochlorines. Some conventional classes of insecticides, such as organochlorines and pyrethroids, were classified as being of high concern. We further identified a group of PB&MT substances that were considered to be both "bioaccumulative" and "mobile". Their properties and common substructures for several major clusters were characterized. The present results prioritized groups of substances with high potentials to cause adverse effects to the environment and humans, many of which have not yet been fully recognized.

A dynamic inventory database for assessing age-, gender-, and route-specific chronic internal exposure to chemicals in support of human exposome research

In this study, we proposed a dynamic inventory database to evaluate chronic internal exposure to chemicals at a population level, which enables users to perform modeling exercises specific to a particular chemical, route of exposure, age group, and gender. The database was built based on the steady-state solution of physiologically based kinetic (PBK) models. The biotransfer factors [BTF, the steady-state ratio between the chemical concentration in human tissues and the average daily dose (ADD) of the chemical] of 931 organic chemicals in major organs and tissues were simulated for a total of 14 population age groups for males and females. The results indicated that infants and children had the highest simulated BTFs of chemicals, and middle-aged adults had the lowest simulated BTFs. The route-specific analysis of the simulated BTFs indicated that the biotransformation half-life and octanol-water partition coefficient of chemicals had a profound impact on the BTFs. Organ- and chemical-specific results indicated that the biotransfer potential of chemicals in human bodies was primarily determined by bio-thermodynamic variables (e.g., lipid contents). In conclusion, the proposed inventory database can be conveniently used to access chronic internal exposure doses of chemicals by multiplying the route-specific ADD values for different population groups. In future studies, we recommend incorporating human biotransformation data, partition coefficients of ionizable chemicals, age-specific vulnerable indicators (e.g., the degree of maturation of immune systems), physiological variations within the same age group (e.g., intensity of daily physical activities), growth rates (i.e., the dilution effect on chemical biotransfer), and all possible target organs of carcinogenicity (e.g., bladder) into the proposed dynamic inventory database to help promote human exposome research.

A State-of-the-Science Review of Interactions of Per- and Polyfluoroalkyl Substances (PFAS) with Renal Transporters in Health and Disease: Implications for Population Variability in PFAS Toxicokinetics

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are ubiquitous in the environment and have been shown to cause various adverse health impacts. In animals, sex- and species-specific differences in PFAS elimination half-lives have been linked to the activity of kidney transporters. However, PFAS molecular interactions with kidney transporters are still not fully understood. Moreover, the impact of kidney disease on PFAS elimination remains unclear.

OBJECTIVES

This state-of-the-science review integrated current knowledge to assess how changes in kidney function and transporter expression from health to disease could affect PFAS toxicokinetics and identified priority research gaps that should be addressed to advance knowledge.

METHODS

We searched for studies that measured PFAS uptake by kidney transporters, quantified transporter-level changes associated with kidney disease status, and developed PFAS pharmacokinetic models. We then used two databases to identify untested kidney transporters that have the potential for PFAS transport based on their endogenous substrates. Finally, we used an existing pharmacokinetic model for perfluorooctanoic acid (PFOA) in male rats to explore the influence of transporter expression levels, glomerular filtration rate (GFR), and serum albumin on serum half-lives.

RESULTS

The literature search identified nine human and eight rat kidney transporters that were previously investigated for their ability to transport PFAS, as well as seven human and three rat transporters that were confirmed to transport specific PFAS. We proposed a candidate list of seven untested kidney transporters with the potential for PFAS transport. Model results indicated PFOA toxicokinetics were more influenced by changes in GFR than in transporter expression.

DISCUSSION

Studies on additional transporters, particularly efflux transporters, and on more PFAS, especially current-use PFAS, are needed to better cover the role of transporters across the PFAS class. Remaining research gaps in transporter expression changes in specific kidney disease states could limit the effectiveness of risk assessment and prevent identification of vulnerable populations. https://doi.org/10.1289/EHP11885.

Generalizing routes of plant exposure to pesticides by plant uptake models to assess pesticide application efficiency

Pesticide application techniques are critical not only for integrated pest management (IPM) but also for food and environmental safety. Assessing pesticide application efficiency on plants can help optimize IPM and reduce pesticide environmental impacts. With hundreds of pesticides registered for use in agriculture, this study proposed a modeling approach based on plant uptake models for generalizing routes of plant chemical exposures that can correspond to different types of pesticide application methods and evaluating their respective efficiency on plants. Three representative pesticide application methods (i.e., drip irrigation, foliar spray, and broadcast application) were selected for modeling simulations. The simulation results for three representative pesticides (i.e., halofenozide, pymetrozine, and paraquat) revealed that the soil-based transpiration exposure route facilitated the bioaccumulation of moderately lipophilic compounds in leaves and fruits. While the plant surface-based exposure route (i.e., leaf cuticle penetration) made it easier for highly lipophilic compounds to enter plants, moderately lipophilic pesticides (i.e., log K_OW ∼ 2) were more soluble in phloem sap, which enhanced their subsequent transport within plant tissues. In general, moderately lipophilic pesticides had the highest simulated residue concentrations in plant tissues for the three specific application methods, indicating they had the highest application efficiency due to their enhanced uptake routes (via transpiration and surface penetration) and increased solubility in xylem and phloem saps. Compared to foliar spray and broadcast application, drip irrigation produced higher residue concentrations for a wide variety of pesticides, exhibiting the highest application efficiency for many pesticides, especially for moderately lipophilic compounds. Future research should incorporate plant growth stages, crop safety, pesticide formulations, and multiple application events into the modeling approach for understanding pesticide application efficiency evaluation.

Partitioning into phosphatidylcholine-cholesterol membranes: liposome measurements, coarse-grained simulations, and implications for bioaccumulation

Membrane-water partitioning is an important physical property for the assessment of bioaccumulation and environmental impact. Here, we advance simulation methodology for predicting the partitioning of small molecules into lipid membranes and compare the computational predictions to experimental measurements in liposomes. As a step towards high-throughput screening, we present an automated mapping and parametrization procedure to produce coarse-grained models compatible with the Martini 3 force field. The methodology is general and can also be used for other applications where coarse-grained simulations are appropriate. This article addresses the effect on membrane-water partitioning of adding cholesterol to POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Nine contrasting neutral, zwitterionic and charged solutes are tested. Agreement between experiment and simulation is generally good, with the most challenging cases being permanently charged solutes. For all solutes, partitioning is found to be insensitive to membrane cholesterol concentration up to 25% mole fraction. Hence, for assessment of bioaccumulation into a range of membranes (such as those found in fish), partitioning data measured in pure lipid membranes are still informative.

Investigating the bioaccumulation potential of anionic organic compounds using a permanent rainbow trout liver cell line

Permanent rainbow trout (Oncorhynchus mykiss) cell lines represent potential in vitro alternatives to experiments with fish. We here developed a method to assess the bioaccumulation potential of anionic organic compounds in fish, using the rainbow trout liver-derived RTL-W1 cell line. Based on the availability of high quality in vivo bioconcentration (BCF) and biomagnification (BMF) data and the substances' charge state at physiological pH, four anionic compounds were selected: pentachlorophenol (PCP), diclofenac (DCF), tecloftalam (TT) and benzotriazol-tert-butyl-hydroxyl-phenyl propanoic acid (BHPP). The fish cell line acute toxicity assay (OECD TG249) was used to derive effective concentrations 50 % and non-toxic exposure concentrations to determine exposure concentrations for bioaccumulation experiments. Bioaccumulation experiments were performed over 48 h with a total of six time points, at which cell, medium and plastic fractions were sampled and measured using high resolution tandem mass spectrometry after online solid phase extraction. Observed cell internal concentrations were over-predicted by KOW-derived predictions while pH-dependent octanol-water partitioning (DOW) and membrane lipid-water partitioning (DMLW) gave better predictions of cell internal concentrations. Measured medium and cell internal concentrations at steady state were used to calculate RTL-W1-based BCF, which were compared to DOW- or DMLW-based model approaches and in vivo data. With the exception of PCP, the cell-derived BCF best compared to DOW-based model predictions, which were higher than predictions based on DMLW. All methods predicted the in vivo BCF for diclofenac well. For PCP, the cell-derived BCF was lowest although all BCF predictions underestimated the in vivo BCF by ≥ 1 order of magnitude. The RTL-W1 cells, and all other prediction methods, largely overestimated in vivo BMF, which were available for PCP, TT and BHPP. We conclude that the RTL-W1 cell line can supplement BCF predictions for anionic compounds. For BMF estimations, however, in vitro-in vivo extrapolations need adaptation or a multiple cell line approach.

Influence of habitats and physicochemical factors on trophic transfer processes of antibiotics in a freshwater ecosystem: Application of stable isotopes and human health risks

Habitats of species and physicochemical factors are of great importance in determining the trophic transfer of contaminants in freshwater ecosystems. There is little information on how multiple physicochemical factors and habitats influence the trophic transfer of antibiotics in freshwater food webs. This study investigated the concentrations of 7 sulfonamides (SAs), 4 fluoroquinolones (FQs), 4 tetracyclines (TCs) and 2 macrolides (MLs) in the Lake Dianshan food web. Stable isotope analysis (SIA), and mathematical models were used to assign trophic levels and distinguish between the benthic food web (BFW) and pelagic food web (PFW). Values of stable nitrogen isotope (δ¹⁵N‰) and stable carbon isotope (δ¹³C‰) ranged from 10.2 ± 0.11 to 19.72 ± 0.05 and -33.67 ± 0.18 to -20.79 ± 0.50, respectively. Total concentrations of antibiotics ranged from 36.63 ± 12.73 ng/g dry weight (dw) to 105.85 ± 12.95 ng/g dw for all species. The relative abundance of antibiotics was in the following order: ∑FQs (36.49 %) > ∑SAs (26.70 %), >∑MLs (12.63 %) for all biotas. Trophic magnification factor (TMFs) values for individual antibiotics ranged from 0.10 to 1.20 and 0.31 to 1.82 for PFW and BFW, respectively. Three classes of antibiotics ∑FQs (p < 0.05), ∑TCs (p < 0.05), and ∑MLs (p < 0.05) showed significant trophic dilution in PFW, opposite to non-significant trophic dilution in BFW. The influence of various physicochemical factors was not strong over trophic transfer (e.g., octanol-water partition coefficient-LogK_ow (r = -0.05 in PFW, r = -0.14 in BFW) and distribution coefficient-LogD (r = 0.06 in PWF, r = -0.28 in BFW)) except for aqueous solubility (LogS). Results indicated a significantly higher trophic dilution of antibiotics in the PFW than in the BFW. Among the studied six physicochemical factors, only LogS significantly influences (p < 0.05) the trophic transfer of antibiotics in the freshwater food web. Health risk assessments indicated that currently, there were no serious risks present for urban and rural populations.

Global occurrence and aquatic hazards of antipsychotics in sewage influents, effluent discharges and surface waters

Despite increasing reports of pharmaceuticals in surface waters, aquatic hazard information remains limited for many contaminants, particularly for sublethal, chronic responses plausibly linked to molecular initiation events that are largely conserved across vertebrates. Here, we critically examined available refereed information on the occurrence of 67 antipsychotics in wastewater effluent and surface waters. Because the majority of sewage remains untreated around the world, we also examined occurrence in sewage influents. When sufficient information was available, we developed probabilistic environmental exposure distributions (EEDs) for each compound in each matrix by geographic region. We then performed probabilistic environmental hazard assessments (PEHAs) using therapeutic hazard values (THVs) of each compound, due to limited sublethal aquatic toxicology information for this class of pharmaceuticals. From these PEHAs, we determined predicted exceedances of the respective THVs for each chemical among matrices and regions, noting that THV values of antipsychotic contaminants are typically lower than other classes of human pharmaceuticals. Diverse exceedances were observed, and these aquatic hazards varied by compound, matrix and geographic region. In wastewater effluent discharges and surface waters, sulpiride was the most detected antipsychotic; however, percent exceedances of the THV were minimal (0.6%) for this medication. In contrast, we observed elevated aquatic hazards for chlorpromazine (30.5%), aripiprazole (37.5%), and perphenazine (68.7%) in effluent discharges, and for chlorprothixene (35.4%) and flupentixol (98.8%) in surface waters. Elevated aquatic hazards for relatively understudied antipsychotics were identified, which highlight important data gaps for future environmental chemistry and toxicology research.

Tissue-specific accumulation, depuration, and effects of perfluorooctanoic acid on fish: Influences of aqueous pH and sex

Perfluorooctanoic acid (PFOA) is widely distributed in nature, particularly in aquatic environments. Its bioaccumulation and toxicity in aquatic organisms can be affected by both the chemical status of PFOA in water and the physiology of the organism. However, research on the patterns of these effects is scarce. In this study, we investigated the influence of aqueous pH (pH 6, acidic; pH 7.5, neutral; pH 9, basic) and fish sex on PFOA uptake, clearance, and biochemical effects in crucian carp (C. auratus) using flow-through exposure. In the 17-d kinetic experiment, PFOA bioaccumulation adhered to a uniform first-order model in which PFOA uptake rates from water to blood and liver in acidic conditions were faster than those in other conditions, indicating possible acidic pH influence on PFOA uptake. PFOA clearance rates in these compartments of males were slower than in females, which was attributed to the notably stronger expression of Oat2 (organic anion transporter 2, responsible for reabsorption) in the kidneys of males. Similar responses were observed for peroxisome proliferation-related biomarkers at different pH levels and in different sexes. These biochemical responses were driven by the internal concentrations of PFOA. The inhibition acetylcholinesterase activity in the fish brain was closely linked to changes in P-glycoprotein content, demonstrating a protective relationship. Collectively, an aqueous pH lower than 7.5 might affect the uptake of PFOA by fish. The clearance discrepancies between the sexes highlight the importance of anion carriers for ionizable organic compounds in aquatic organisms.

Perfluorooctanoic acid (PFOA) exposure in relation to the kidneys: A review of current available literature

Perfluorooctanoic acid is an artificial and non-degradable chemical. It is widely used due to its stable nature. It can enter the human body through food, drinking water, inhalation of household dust and contact with products containing perfluorooctanoic acid. It accumulates in the human body, causing potential harmful effects on human health. Based on the biodegradability and bioaccumulation of perfluorooctanoic acid in the human body, there are increasing concerns about the adverse effects of perfluorooctanoic acid exposure on kidneys. Research shows that kidney is the main accumulation organ of Perfluorooctanoic acid, and Perfluorooctanoic acid can cause nephrotoxicity and produce adverse effects on kidney function, but the exact mechanism is still unknown. In this review, we summarize the relationship between Perfluorooctanoic acid exposure and kidney health, evaluate risks more clearly, and provide a theoretical basis for subsequent research.

Review of per- and polyfluoroalkyl substances (PFAS) bioaccumulation in earthworms

Per- and polyfluoroalkyl substances (PFAS) are widely used across the globe in commercial products such textiles, firefighting foams, and surface coatings. Some PFAS, such as perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA), are known to be bioaccumulative. Numerous terrestrial ecosystems including sites near PFAS manufacturing facilities, facilities using PFAS in their manufacturing processes, firefighting training areas, landfills, and agricultural fields treated with some pesticide formulations, have been contaminated with PFAS. Earthworms reside at the base of the terrestrial food chain and to perform risk assessments at terrestrial sites contaminated with PFAS, information on the bioaccumulation of PFAS is needed. To understand the bioaccumulation of PFAS by earthworms, a literature search was performed, and biota-soil accumulation factors (BSAFs), measured in laboratory tests and at field sites contaminated with PFAS, were assembled and evaluated in this review. Based on this review, we conclude that there is enough data available for carboxylic and sulfonic acid PFAS classes to derive useful BSAFs for terrestrial risk assessments. Laboratory tests with PFOS and PFOA will be close to or at steady-state conditions with standardized testing protocols, and for the longer chain carboxylic and sulfonic acids, it is unlikely they will reach steady-state with the completion of the uptake exposure. For PFAS classes beyond the carboxylic and sulfonic acids, data are limited and performing terrestrial risk assessments with these PFAS will be difficult. Lastly, additional measurements are needed for non-acid PFAS classes as well as from field settings for all PFAS classes. Across all studies, PFOS and PFOA had average (standard deviation, count) BSAFs (kg-OC/kg-ww) of 0.167 (0.311, 60) and 0.0413 (0.175, 47), respectively.

Factors Affecting the Binding of Diltiazem to Rainbow Trout Plasma: Implications for the Risk Assessment of Pharmaceuticals in Aquatic Systems

The accumulation of organic toxicants in fish plasma, and how they partition between the bound and unbound fraction once absorbed, are important metrics in models that seek to predict the risk of such contaminants in aquatic settings. Rapid equilibrium dialysis of diltiazem, an ionizable weak base and important human pharmaceutical contaminant of freshwaters, was conducted with rainbow trout (Oncorhynchus mykiss) plasma. The effect of fed state, fish sex, fish strain/size, and dialysis buffer pH on the binding of radiolabeled diltiazem (9 ng ml^-1 ) was assessed. In fed fish, 24.6%-29.5% of diltiazem was free, unbound to plasma proteins. Although starvation of fish resulted in a decrease in plasma protein, the bound fraction of diltiazem remained relatively constant. Consequently, the protein-bound concentration of diltiazem increased with length of starvation. In general, rainbow trout strain was a significant factor affecting plasma binding, although the two strains tested also differed markedly in size. Dialysis buffer pH significantly influenced plasma binding, with a higher unbound diltiazem fraction at pH 6.8 than pH 8.0. These data indicate that empirical measures of plasma binding in fish are important for accurate risk assessment and that the physiological status of a fish is likely to impact its sensitivity to toxicants such as diltiazem. Environ Toxicol Chem 2022;41:3125-3133. © 2022 SETAC.

Predicting Hepatic Clearance of Psychotropic Drugs in Isolated Perfused Fish Livers Using a Combination of Two In Vitro Assays

In vitro biotransformation assays with primary trout hepatocytes (RT-HEP) or liver subcellular fractions (RT-S9) have been proposed as valuable tools to help scientists and regulators better understand the toxicokinetics of chemicals. While both assays have been applied successfully to a diversity of neutral organic chemicals, only the RT-S9 assay has been applied to a large number of ionizable organic chemicals. Here, a combination of an in vitro biotransformation assay with RT-HEP with an active transport assay based on the permanent rainbow trout liver cell line RTL-W1 was used to qualitatively predict the potential hepatic clearance of nine psychotropic drugs with various degrees of ionization. Predictions were compared with rates of clearance measured in isolated perfused rainbow trout livers, and the importance of active transport was verified in the presence of the active transport inhibitor cyclosporin A. For the first time, it was demonstrated that a combination of biotransformation and active transport assays is powerful for the prediction of rates of hepatic clearance of ionizable chemicals. Ultimately, it is expected that this approach will allow for use of fewer animals while at the same time improving our confidence in the use of data from in vitro assays in chemical risk assessment.

The Need for Testing Isomer Profiles of Perfluoroalkyl Substances to Evaluate Treatment Processes

Many environmentally relevant poly-/perfluoroalkyl substances (PFASs) including perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) exist in different isomeric (branched and linear) forms in the natural environment. The isomeric distribution of PFASs in the environment and source waters is largely controlled by the source of contamination and varying physicochemical properties imparted by their structural differences. For example, branched isomers of PFOS are relatively more reactive and less sorptive compared to the linear analogue. As a result, the removal of branched and linear PFASs during water treatment can vary, and thus the isomeric distribution in source waters can influence the overall efficiency of the treatment process. In this paper, we highlight the need to consider the isomeric distribution of PFASs in contaminated matrices while designing appropriate remediation strategies. We additionally summarize the known occurrence and variation in the physicochemical properties of PFAS isomers influencing their detection, fate, toxicokinetics, and treatment efficiency.

Effect of solutes structure and pH on the n-octanol/water partition coefficient of ionizable organic compounds

Liquid-liquid partition coefficient is a useful tool to predict biological and environmental fate of organic compounds, for example bioaccumulation or toxicity of lipophilic contaminants. Conversely, the partitioning of ionizable compounds is poorly studied in contrast to that of neutral compounds. Yet, such topic deserves attention, since numerous organic contaminants are ionizable as well as their degradation products. Hence, the contribution of charged species has to be considered in order to model accurately the mass balance or partition of ionizable compounds. In this context, we investigated the liquid-liquid partition of 13 ionizable compounds (oxalic acid, histidine, benzimidazole, etc.), covering various classes of compounds (carboxylic acids, amino-acids, etc.). The n-octanol/water partition coefficient was measured from pH 1 up to 13, in order to fully gather the distribution of both neutral and charged species. Empirical models describing these results are reviewed and partition parameters adjusted for charged species. The study of benzoic acid derivatives (benzoic, salicylic, ortho- and iso-phthalic acids) provides insights on the influence of chemical groups on the partitioning. In the case of tryptophan, the use of acid/base microconstants allowed to estimate the partition of both the zwitterion and its neutral tautomer. Despite a major zwitterionic form (log P^Z(tryptophan) = -1.58 ± 0.30), the minor but neutral tautomer (log P^N(tryptophan) = +0.03 ± 0.30) drives the partition equilibrium. Overall, the provided data may be useful to assess the retention of contaminants, its dependency on pH and salinity variations, and thus understanding their environmental fate. Such data may also be useful as well for molecular simulation involving solvation of organic ions in aqueous and non-aqueous solvents.

Identifying organic chemicals not subject to bioaccumulation in air-breathing organisms using predicted partitioning and biotransformation properties

Because the respiration processes contributing to the elimination of organic chemicals deviate between air- and water-breathing organisms, existing and widely used procedures for identifying chemicals not subject to bioaccumulation in aquatic organisms based on the octanol-water partition ratio K_OW need to be complemented with similar procedures for organisms respiring air. Here, we propose such a procedure that relies on the comparison of a compound's predicted K_OW , octanol-air partition ratio K_OA , and biotransformation half-life HL_B with three threshold values, below which elimination is judged to be sufficiently rapid to prevent bioaccumulation. The method allows for the consideration of the effect of dissociation on the efficiency of urinary and respiratory elimination. Explicit application of different types of the prediction error, such as the 95% prediction interval or the standard error, allows for variable tolerance for false-negative decisions, that is, the potential to judge a chemical as not bioaccumulative even though it is. A test with a set of more than 1000 diverse organic chemicals confirms the applicability of the prediction methods for a wide range of compounds and the procedure's ability to categorize approximately four-fifth of compounds as being of no bioaccumulation concern, suggesting its usefulness to screen large numbers of commercial chemicals to identify those worthy of further scrutiny. The test also demonstrates that a screening based solely on K_OW and K_OA would be far less effective because the fraction of chemicals that can be judged as sufficiently volatile and/or sufficiently water soluble for rapid respiratory and urinary elimination based on the partitioning properties predicted for their neutral form is relatively small. Future improvements of the proposed procedure depend largely on the development of prediction methods for the biotransformation kinetics in air-breathing organisms and for the potential for renal reabsorption. Integr Environ Assess Manag 2022;18:1297-1312. © 2021 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Advances in computational methods along the exposure to toxicological response paradigm

Human health risk assessment is a function of chemical toxicity, bioavailability to reach target biological tissues, and potential environmental exposure. These factors are complicated by many physiological, biochemical, physical and lifestyle factors. Furthermore, chemical health risk assessment is challenging in view of the large, and continually increasing, number of chemicals found in the environment. These challenges highlight the need to prioritize resources for the efficient and timely assessment of those environmental chemicals that pose greatest health risks. Computational methods, either predictive or investigative, are designed to assist in this prioritization in view of the lack of cost prohibitive in vivo experimental data. Computational methods provide specific and focused toxicity information using in vitro high throughput screening (HTS) assays. Information from the HTS assays can be converted to in vivo estimates of chemical levels in blood or target tissue, which in turn are converted to in vivo dose estimates that can be compared to exposure levels of the screened chemicals. This manuscript provides a review for the landscape of computational methods developed and used at the U.S. Environmental Protection Agency (EPA) highlighting their potentials and challenges.

A food web bioaccumulation model for the accumulation of per- and polyfluoroalkyl substances (PFAS) in fish: how important is renal elimination?

Per- and polyfluoroalkyl substances (PFAS) are a large class of highly fluorinated anthropogenic chemicals. Some PFAS bioaccumulate in aquatic food webs, thereby posing risks for seafood consumers. Existing models for persistent organic pollutants (POPs) perform poorly for ionizable PFAS. Here we adapt a well-established food web bioaccumulation model for neutral POPs to predict the bioaccumulation behavior of six perfluoroalkyl acids (PFAAs) and two perfluoroalkyl ether acids (HFPO-DA, 9-Cl-PF3ONS) produced as PFAA replacements. The new model includes sorption to blood plasma proteins and phospholipids, empirically parameterized membrane transport, and renal elimination for PFAAs. Improved performance relative to prior models without these updates is shown by comparing simulations to field and lab measurements. PFAS with eight or more perfluorinated carbons (η_pfc ≥ 8, i.e., C8 perfluorosulfonic acid, C10-C11 perfluorocarboxylic acid, 9-Cl-PF3ONS) are often the most abundant in aquatic food webs. The new model reproduces their observed bioaccumulation potential within a factor of two for >80% of fish species, indicating its readiness to support development of fish consumption advisories for these compounds. Results suggest bioaccumulation of η_pfc ≥ 8 PFAS is primarily driven by phospholipid partitioning, and that renal elimination is negligible for these compounds. However, specific protein binding mechanisms are important for reproducing the observed tissue concentrations of many shorter-chain PFAAs, including protein transporter-mediated renal elimination. Additional data on protein-binding and membrane transport mechanisms for PFAS are needed to better understand the biological behavior of shorter-chain PFAAs and their alternatives.

Human Apparent Volume of Distribution Predicts Bioaccumulation of Ionizable Organic Chemicals in Zebrafish Embryos

Chemicals with elevated bioaccumulation profiles present potential hazards to public health and the environment. Ionizable organic compounds (IOCs) increasingly represent a large proportion of commercial chemicals; however, historical approaches for bioaccumulation determinations are mainly developed for neutral chemicals, which were not appropriate for IOCs. Herein, we employed the zebrafish embryo, a common vertebrate model in environmental and biomedical studies, to elucidate toxicokinetics and bioconcentration of eight IOCs with diverse physicochemical properties and pharmacokinetic parameters. At an environmentally relevant pH (7.5), most IOCs exhibited rapid uptake and depuration in zebrafish, suggesting the ionized forms of IOCs are readily bioavailable. Bioconcentration factors (BCF) of these IOCs ranged from 0.0530 to 250 L·kg^-1 wet weight. The human pharmacokinetic proportionality factor, apparent volume of distribution (V_D), better predicted the BCF of selected IOCs than more commonly used hydrophobicity-based parameters (e.g., pH-dependent octanol-water distribution ratio, D_ow). Predictive bioaccumulation models for IOCs were constructed and validated using V_D alone or with D_ow. Significant relationships between fish BCF and human V_D, which is readily available for pharmaceuticals, highlighted the utility of biologically based "read-across" approaches for predicting bioaccumulative potential of IOCs. Our novel findings thus provided an understanding of the partitioning behavior and improved predictive bioconcentration modeling for IOCs.

Bioaccumulation of per- and polyfluoroalkyl substance in fish from an urban river: Occurrence, patterns and investigation of potential ecological drivers

Per- and polyfluoroalkyl substances (PFAS) are ubiquitous in aquatic environments and a recent shift toward emerging PFAS is calling for new data on their occurrence and fate. In particular, understanding the determinants of their bioaccumulation is fundamental for risk assessment purposes. However, very few studies have addressed the combined influence of potential ecological drivers of PFAS bioaccumulation in fish such as age, sex or trophic ecology. Thus, this work aimed to fill these knowledge gaps by performing a field study in the Seine River basin (France). Composite sediment and fish (European chub, Squalius Cephalus) samples were collected from four sites along a longitudinal transect to investigate the occurrence of 36 PFAS. Sediment molecular patterns were dominated by fluorotelomer sulfonamidoalkyl betaines (i.e. 6:2 and 8:2 FTAB, 46% of ∑PFAS on average), highlighting the non-negligible contribution of PFAS of emerging concern. C₉-C₁₄ perfluoroalkyl carboxylic acids, perfluorooctane sulfonic acid (PFOS), perfluorooctane sulfonamide (FOSA) and 10:2 fluorotelomer sulfonate (10:2 FTSA) were detected in all fish samples. Conversely, 8:2 FTAB was detected in a few fish from the furthest downstream station only, suggesting the low bioaccessibility or the biotransformation of FTABs. ∑PFAS in fish was in the range 0.22-3.8 ng g^-1 wet weight (ww) and 11-140 ng g^-1 ww for muscle and liver, respectively. Fish collected upstream of Paris were significantly less contaminated than those collected downstream, pointing to urban and industrial inputs. The influence of trophic ecology and biometry on the interindividual variability of PFAS burden in fish was examined through analyses of covariance (ANCOVAs), with sampling site considered as a categorical variable. While the latter was highly significant, diet was also influential; carbon sources and trophic level (i.e. estimated using C and N stable isotope ratios, respectively) equally explained the variability of PFAS levels in fish.

Metabolome adaptation and oxidative stress response of common carp (Cyprinus carpio) to altered water pollution levels

Treated wastewater ponds (TWPs) serve as recipients and passive tertiary treatment mediators for recycled water. These nutrient-rich habitats are increasingly utilised in aquaculture, nevertheless multiple loads of various contaminants with adverse effects on aquatic fauna, including fish, have been recorded. In the present study, we investigated the effects of fish transfer in response to altered levels of pollution on liver metabolic profiles and tissue-specific oxidative stress biomarkers during short- and long-term exposure. In a field experiment, common carp (Cyprinus carpio) originating in severely polluted TWP were restocked after one year to a reference pond with a background pollutant concentration typical of the regional river. In contrast, fish that originated in the reference pond were restocked to TWP. Fish were sampled 0, 7, 14, 60, and 180 days after restocking and fish liver, kidney, intestine, and gill tissues were subjected to biomarker analysis. Pharmaceutically active compounds (PhACs) and metabolic profiles were determined in fish liver using liquid chromatography high-resolution mass spectrometry (LC-HRMS). Fish transferred from reference to polluted pond increased the antioxidant response and absorbed PhACs into metabolism within seven days. Fish liver metabolic profiles were shifted rapidly, but after 180 days to a lesser extent than profiles in fish already adapted in polluted water. Restocked fish from polluted to reference pond eliminated PhACs during the short phase within 14 days, and the highest antioxidant response accompanied the depuration process. Numerous elevated metabolic compounds persisted in such exposed fish for at least 60 days. The period of two weeks was suggested as sufficient for PhACs depuration, but more than two months after restocking is needed for fish to stabilise their metabolism. This study contributed to determining the safe handling with marketed fish commonly restocked to wastewaters and clarified that water pollution irreversibly altered fish metabolic profile.

pH-Dependent Partitioning of Ionizable Organic Chemicals between the Silicone Polymer Polydimethylsiloxane (PDMS) and Water

The silicone polymer polydimethysiloxane (PDMS) is a popular passive sampler for in situ and ex situ sampling of hydrophobic organic chemicals. Despite its limited sorptive capacity for polar and ionizable organic chemicals (IOC), IOCs have been found in PDMS when extracting sediment and suspended particulate matter. The pH-dependent partitioning of 190 organics and IOCs covering a range of octanol-water partition constants log K _ow from -0.3 to 7.7 was evaluated with a 10-day shaking method using mixtures composed of all chemicals at varying ratios of mass of PDMS to volume of water. This method reproduced the PDMS-water partition constant K _PDMS/w of neutral chemicals from the literature and extended the dataset by 93 neutral chemicals. The existing quantitative structure-activity relationship between the log K _ow and K _PDMS/w could be extended with the measured K _PDMS/w linearly to a log K _ow of -0.3. Fully charged organics were not taken up into PDMS. Thirty-eight monoprotic organic acids and 42 bases showed negligible uptake of the charged species, and the pH dependence of the apparent D _PDMS/w(pH) could be explained by the fraction of neutral species multiplied by the K _PDMS/w of the neutral species of these IOCs. Seventeen multiprotic chemicals with up to three acidity constants pK _a also showed a pH dependence of D _PDMS/w(pH) with the tendency that the neutral and zwitterionic forms showed the highest D _PDMS/w(pH). D _PDMS/w(pH) of charged species of more hydrophobic multiprotic chemicals such as tetrabromobisphenol A and telmisartan was smaller but not negligible. Since these chemicals show high bioactivity, their contribution to mixture effects has to be considered when testing passive sampling extracts with in vitro bioassays. This work has further implications for understanding the role of microplastic as a vector for organic micropollutants.

A Generalized Physiologically Based Kinetic Model for Fish for Environmental Risk Assessment of Pharmaceuticals

An increasing number of pharmaceuticals found in the environment potentially impose adverse effects on organisms such as fish. Physiologically based kinetic (PBK) models are essential risk assessment tools, allowing a mechanistic approach to understanding chemical effects within organisms. However, fish PBK models have been restricted to a few species, limiting the overall applicability given the countless species. Moreover, many pharmaceuticals are ionizable, and fish PBK models accounting for ionization are rare. Here, we developed a generalized PBK model, estimating required parameters as functions of fish and chemical properties. We assessed the model performance for five pharmaceuticals (covering neutral and ionic structures). With biotransformation half-lives (HLs) from EPI Suite, 73 and 41% of the time-course estimations were within a 10-fold and a 3-fold difference from measurements, respectively. The performance improved using experimental biotransformation HLs (87 and 59%, respectively). Estimations for ionizable substances were more accurate than any of the existing species-specific PBK models. The present study is the first to develop a generalized fish PBK model focusing on mechanism-based parameterization and explicitly accounting for ionization. Our generalized model facilitates its application across chemicals and species, improving efficiency for environmental risk assessment and supporting an animal-free toxicity testing paradigm.

New insights on in vitro biotransformation of anticoagulant rodenticides in fish

The assessment of the bioaccumulation potential of chemicals is an essential and mandatory part of their regulatory environmental risk and hazard assessment. So far, in vitro data on fish metabolism is rarely available for biocidal active substances such as anticoagulant rodenticides. In this case study we present in vitro biotransformation rates of eight biocidal and one pharmaceutical anticoagulants in rainbow trout (Oncorhynchus mykiss) liver subcellular S9 fraction (RT-S9) determined following the Organisation for Economic Co-operation and Development test guideline 319B method at two different incubation temperatures (i.e., 12 ± 1 °C and 23 ± 2 °C). Furthermore, we address challenges associated with the usability and interpretation of in vitro data to support the decision making within the regulatory bioaccumulation assessment in bridging the gap between in silico methods and in vivo studies. According to our results, four of the tested substances (i.e., chlorophacinone, coumatetralyl, bromadiolone, and difenacoum) exhibited significant intrinsic clearance (p < .001) in the RT-S9 assay. Overall, the observed metabolism was (very) slow and clearance rates were temperature-dependent. Whether the determined in vitro biotransformation rate had a substantial influence on the predicted bioconcentration factor during extrapolation was subject to the lipophilicity of the test substance. Further improvements of existing concepts are needed to overcome uncertainties in the prediction of bioconcentration factors for chemicals such as anticoagulants.

Sorption and Mobility of Charged Organic Compounds: How to Confront and Overcome Limitations in Their Assessment

Permanently charged and ionizable organic compounds (IOC) are a large and diverse group of compounds belonging to many contaminant classes, including pharmaceuticals, pesticides, industrial chemicals, and natural toxins. Sorption and mobility of IOCs are distinctively different from those of neutral compounds. Due to electrostatic interactions with natural sorbents, existing concepts for describing neutral organic contaminant sorption, and by extension mobility, are inadequate for IOC. Predictive models developed for neutral compounds are based on octanol-water partitioning of compounds (Kow) and organic-carbon content of soil/sediment, which is used to normalize sorption measurements (KOC). We revisit those concepts and their translation to IOC (Dow and DOC) and discuss compound and soil properties determining sorption of IOC under water saturated conditions. Highlighting possible complementary and/or alternative approaches to better assess IOC mobility, we discuss implications on their regulation and risk assessment. The development of better models for IOC mobility needs consistent and reliable sorption measurements at well-defined chemical conditions in natural porewater, better IOC-, as well as sorbent characterization. Such models should be complemented by monitoring data from the natural environment. The state of knowledge presented here may guide urgently needed future investigations in this field for researchers, engineers, and regulators.

Enhancing Scientific Support for the Stockholm Convention's Implementation: An Analysis of Policy Needs for Scientific Evidence

The Stockholm Convention is key to addressing the global threats of persistent organic pollutants (POPs) to humanity and the environment. It has been successful in identifying new POPs, but its national implementation remains challenging, particularly by low- and middle-income Parties. Concerted action is needed to assist Parties in implementing the Convention's obligations. This analysis aims to identify and recommend research and scientific support needed for timely implementation of the Convention. We aim this analysis at scientists and experts from a variety of natural and social sciences and from all sectors (academia, civil society, industry, and government institutions), as well as research funding agencies. Further, we provide practical guidance to scientists and experts to promote the visibility and accessibility of their work for the Convention's implementation, followed by recommendations for sustaining scientific support to the Convention. This study is the first of a series on analyzing policy needs for scientific evidence under global governance on chemicals and waste.

Trout and Human Plasma Protein Binding of Selected Pharmaceuticals Informs the Fish Plasma Model

Concerns are increasing that pharmaceuticals released into the environment pose a risk to nontarget organism such as fish. The fish plasma model is a read-across approach that uses human therapeutic blood plasma concentrations for estimating likely effects in fish. However, the fish plasma model neglects differences in plasma protein binding between fish and humans. Because binding data for fish plasma are scarce, the binding of 12 active pharmaceutical ingredients (APIs; acidic, basic, and neutral) to rainbow trout (Oncorhynchus mykiss) and human plasma was measured using solid-phase microextraction (SPME). The plasma/water distribution ratios (D _plasma/w ) of neutral and basic APIs were similar for trout and human plasma, differing by no more than a factor of 2.7 for a given API. For the acidic APIs, the D _plasma/w values of trout plasma were much lower than for human plasma, by up to a factor of 71 for naproxen. The lower affinity of the acidic APIs to trout plasma compared with human plasma suggests that the bioavailability of these APIs is higher in trout. Read-across approaches like the fish plasma model should account for differences in plasma protein binding to avoid over- or underestimation of effects in fish. For the acidic APIs, the effect ratio of the fish plasma model would increase by a factor of 5 to 60 if the unbound plasma concentrations were used to calculate the effect ratio. Environ Toxicol Chem 2022;41:559-568. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Evidence for Specific Receptor-Mediated Toxicity of Pharmaceuticals in Aquatic Organisms Derived from Acute and Chronic Standard Endpoints

The toxicity of 17 active pharmaceutical ingredients (APIs) was investigated using standardized acute and chronic tests with Daphnia magna and 2 algae species. Chronic toxicity was generally greater for Daphnia than for algae. Compilation of additional data resulted in 100 APIs for which the acute-to-chronic ratio (ACR) was determined for Daphnia. The frequency of high ACRs (~20% with ACRs > 100) indicates that specific receptor-mediated toxicity toward D. magna is rather common among APIs. The 11 APIs with ACRs > 1000 included lipid-modifying agents, immunosuppressants, antibiotics, antineoplastics, antiobesics, antivirals, and antihistamines. There was no consistent association between ACR and chronic toxicity, ionization status, or lipophilicity. High ACRs were not exclusively associated with the presence of orthologs of the pharmacological target in Daphnia. Statins, acetylcholinesterase inhibitors, and antihistamines are discussed in more detail regarding the link between targets and toxic mode of action. For acetylcholinesterase inhibitors, receptor-mediated toxicity was already apparent after acute exposure, whereas the high ACR and chronic toxicity of some antihistamines probably related to interaction with a secondary rather than the primary pharmacological target. Acute or modeled chronic toxicity estimates have often been used for prioritizing pharmaceuticals. This may be seriously misleading because chronic effects are currently not predictable for APIs with specific receptor-mediated toxicity. However, it is exactly these APIs that are the most relevant in terms of environmental risks. Environ Toxicol Chem 2022;41:601-613. © 2021 SETAC.

Distribution, sources, and dietetic-related health risk assessment of perfluoroalkyl acids (PFAAs) in the agricultural environment of an industrial-agricultural interaction region (IAIR), Changshu, East China

The exploration of the distribution and dietetic-related health risks of perfluoroalkyl acids (PFAAs) in industrial-agricultural interaction regions (IAIRs) is of significant importance, due to the transfer of many PFAA-related factories to developing countries with intensive agricultural activities. In the present study, based on the local diet, edible parts of rice, vegetables, fish, and their corresponding soils and irrigation/aquaculture water were investigated in a typical Chinese city (Changshu). The concentrations of total perfluoroalkyl acids (ΣPFAAs) in the edible parts of rice /vegetables and fish tissues ranged from 26.69 to 37.09 ng/g dw, 12.93 to 40.77 ng/g dw, and 13.27 to 29.82 ng/g ww, with perfluorohexanoic acid (PFPeA) and perfluorooctane sulfonic acid (PFOS) as the most dominant compounds. The PFAA concentrations in the corresponding rice soils, vegetable soils, irrigation water, and aquaculture water ranged from 11.99 to 26.33 ng/g dw, 14.06 to 36.19 ng/g dw, 141.36 to 297.00 ng/L, and 179.23 to 235.82 ng/L, respectively. Biota-sediment accumulation factor (BSAF) values for the plant-soil system were far greater than those for bioaccumulation factor (BAF) values for the plant-irrigation water system. PFAAs were more inclined to accumulate in the gills of fish as determined by their highest BAF values. Correlation analysis showed that PFAAs in root vegetables had a stronger correlation with those in soil compared with those in irrigation water. Source analysis showed that emissions from fluoride industries, textiles, and food industries may be the dominant sources of PFAAs in agricultural environments. The estimated dietary intake (EDI) for the selected diet was lower than that for rice/vegetables but was higher than that found in fish. Toddlers (2-5 years) had the highest exposure risk, and rural residents were more exposed to PFAAs than urban residents under the selected diet.