Development of physiologically based toxicokinetic models for improving the human indoor exposure assessment to water contaminants: trichloroethylene and trihalomethanes

Comparison of Trihalomethane exposure assessment metrics in epidemiologic analyses of reproductive and developmental outcomes

BACKGROUND

Researchers have developed exposure assessment metrics for disinfection by-products (DBPs) utilizing drinking water monitoring data and accounting for spatial and temporal variability, water consumption, and showering and bathing time with an expectation of decreasing exposure misclassification compared to the use of measured concentrations at public water supply (PWS) monitoring locations alone.

OBJECTIVE

We used exposure data collected for a previous study of DBPs to evaluate how different sources of information impact trihalomethane (THM) exposure estimates.

METHODS

We compared gestational exposure estimates to THMs based on water utility monitoring data alone, statistical imputation of daily concentrations to incorporate temporal variability, and personal water consumption and use (bathing and showering). We used Spearman correlation coefficients and ranked kappa statistics to compare exposure classifications.

RESULTS

Exposure estimates based on measured or imputed daily THM concentrations, self-reported consumption, or bathing and showering differed substantially from estimates based solely on concentrations from PWS quarterly monitoring reports. Ranked exposure classifications, high to low quartiles or deciles, were generally consistent across each exposure metric (i.e., a subject with "high" exposure based on measured or imputed THM concentrations generally remained in the "high" category across exposure metrics.) The measured concentrations and imputed daily (i.e., spline regression) concentrations were highly correlated (r = 0.98). The weighted kappa statistics comparing exposure estimates using different exposure metrics ranged from 0.27 to 0.89, with the highest values for the ingestion + bathing/showering metrics compared to metrics for bathing/showering only (0.76 and 0.89). Bathing and showering contributed the most to "total" THM exposure estimates.

IMPACT STATEMENT

We compare exposure metrics capturing temporal variability and multiple estimates of personal THM exposure with THM concentrations from PWS monitoring data. Our results show exposure estimates based on imputed daily concentrations accounting for temporal variability were very similar to the measured THM concentrations. We observed low agreement between imputed daily concentrations and ingestion-based estimates. Considering additional routes of exposure (e.g., inhalation and dermal) slightly increased agreement with the measured PWS exposure estimate in this population. Overall, the comparison of exposure assessment metrics allows researchers to understand the added value of additional data collection for future epidemiologic analyses of DBPs.

Inter-individual exposure variability interpretation through reflection of biological age algorithm in physiologically based toxicokinetic model: Application to human risk assessment of di-isobutyl-phthalate

Age-related changes and interindividual variability in the degree of exposure to hazardous substances in the environment are pertinent factors to be considered in human risk assessment. Existing risk assessments remain in a one-size-fits-all approach, often without due consideration of inter-individual toxicokinetic variability factors, such as age. The purpose of this study was to advance from the existing risk assessment of hazardous substances based on toxicokinetics to a precise human risk assessment by additionally considering the effects of physiologic and metabolic fluctuations and interindividual variability in age. Qualitative age-associated physiologic and metabolic changes in humans, obtained through a meta-analysis, were quantitatively modeled to produce the final biological age algorithm (BAA). The developed BAAs (for males) were extended and applied to the reported testicular reproductive toxicity-focused di-isobutyl-phthalate (DiBP)-mono-isobutyl-phthalate (MiBP) physiologically based toxicokinetic (PBTK) model in males. The advanced PBTK model combined with the BAA was applied to the human risk assessment based on MiBP biomonitoring data. As a result, the specialized DiBP external exposure values for each age could be estimated. Additionally, by applying the Monte Carlo simulation, the distribution of internal exposure diversity among individuals according to the same external exposure dose could be estimated. The contributions of physiologic and metabolic factors to the age-dependent toxicokinetic changes were approximately 93.41-99.99 and 0.01-6.59%, respectively. In addition, the relative contribution of metabolic factors was major in infants and continued to decrease as age increased (up to about age 30 years). This study provides a step-by-step platform that can be widely applied to overcome the limitations of existing toxicokinetic models that still require interindividual pharmacokinetic variability explanations. This will be important for the rationalization and explanation of inter-individual variability in the pharmacokinetics of many substances.

Bromoform exposure is associated with non-melanoma skin cancer: evidence from NHANES 2011-2020

Background

Non-melanoma skin cancer (NMSC) is a prevalent skin malignancy. It has been indicated in many studies that trihalomethanes (THMs) exposure has a strong association with tumors but has not been associated with NMSC. Our investigation aims to explore the association between THMs exposure and NMSC.

Methods

Cross-sectional data from the 2011 to 2020 National Health and Nutrition Examination Survey (NHANES) was collected. Poisson regression and subgroup analyses were performed to evaluate the association between individual THMs components and NMSC. Fitted smoothing curves and generalized additive models were also used.

Results

This study involved 5,715 individuals, 98 (1.7%) of whom self-reported NMSC. After adjusting for covariates, Poisson regression showed that higher blood TBM levels were associated with an increased likelihood of NMSC (OR = 1.03; 95% CI: 1.01-1.05, p = 0.002). However, the correlation between the blood levels of TCM, DBCM, and BDCM and the likelihood of NMSC was not statistically significant (all p > 0.05). Subgroup analysis and interaction tests showed no significant differences between blood TBM concentration and the likelihood of NMSC, indicating that age, gender, and race were significantly independent of this positive association (all p < 0.05).

Conclusions

Our results implied that among adults older than 65 years old in the U.S., elevated blood TBM concentrations were positively associated with NMSC. More prospective investigations are required to validate this relationship with the early prevention of NMSC.

Rapid determination and health risk assessment of neonicotinoids in source water and tap water of the tropical Hainan Island, China

Neonicotinoids (NEOs) pesticides are widely used around the world, especially in the tropics with greater frequency and intensity. However, little is known about NEOs residue in drinking water of tropics. In this study, a highly efficient method using ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) was established for determining eight NEOs in source water and tap water of Hainan Island, China. The method adopted a high-throughput direct aqueous injection without sample concentration steps, with a rapid analyzing period of 5.0 min, method detection limits (MDLs) in the range of 0.84-1.82 ng/L and the average recoveries ranged from 83% to 116%. NEOs were detected in all source water samples and at an upper level as compared with other parts of China. The most frequently detected NEO was imidacloprid with a detection frequency of 94%, followed by clothianidin (88%) and thiamethoxam (78%), with maximum concentrations of 86.4, 164, and 188 ng/L, respectively. Moreover, seasonal and spatial variations had remarkable impacts on NEO contamination in source water. Drinking water treatment processes removed approximately 20% of NEOs from surface water. However, 90% of tap water samples contained at least one NEO, With 3 samples' concentration of single NEO exceeding the acceptable value recommended by the European Union (100 ng/L). Therefore, the risk of human exposure through drinking water was evaluated for 4 age group and 2 genders. Young children aged 9 months to 3 years old were found to have the highest risk, with the median exposure up to 4 times greater than teenagers and adults. Next, water intake is likely only a small part of the daily intake of these individuals, thus the potential health problems caused by NEOs present in the tap water of Hainan should not be ignored.

Urinary trihalomethane concentrations and liver function indicators: a cross-sectional study in China

While it is known that exposure to disinfection by-products (DBPs), including trihalomethanes (THMs), impairs liver function, few epidemiological studies have explored this association. Here, we determined the concentrations of four urinary trihalomethanes (chloroform [TCM], and three Br-THMs, bromodichloromethane [BDCM], dibromochloromethane [DBCM], and bromoform [TBM]), and nine serum liver function indicators in 182 adults ≥ 18 years of age, examined at a medical examination center in Wuxi, China, in 2020 and 2021. Generalized linear model analysis revealed positive associations between urinary DBCM and alanine aminotransferase (ALT), aspartate aminotransferase (AST), total bilirubin (TBIL), total protein (TP), and albumin (ALB). Urinary Br-THMs and total THMs (TTHMs) were positively associated with ALT, AST, TBIL, indirect bilirubin (IBIL), TP, and ALB (all P < 0.05). Urinary THMs were not associated with alkaline phosphatase (ALP) or glutamine transaminase (GGT) (all P > 0.05). Generalized additive model-based penalized regression splines were used to confirm these associations. In conclusion, THM exposure was associated with altered serum biomarkers of liver function.

Cardiopulmonary parameters and organ blood flows for workers expressed in terms of VO2 for use in physiologically based toxicokinetic modeling

Minute ventilation rates (VE), alveolar ventilation rates (VA), cardiac outputs (Q), liver blood flow (LBF) and kidneys blood flows (KBF) for physiologically based toxicokinetic modeling and occupational health risk assessment in active workers have apparently not been determined. Minute energy expenditure rates (E) and oxygen consumption rates (VO₂) in workers during exertions and their aggregate daytime activities are obtained by using open-circuit wearable devices for indirect calorimetry measurements and the doubly labeled water method respectively. Hundreds of E (in kcal/min) and VO₂ (in L of O₂/min) were previously reported for workers. The oxygen uptake factors of 0.2059 ± 0.0019 and 0.2057 ± 0.0018 L of O₂/kcal during postprandial and fasting phases respectively enabled conversion of E into VO₂. Equations determined in this study based upon more than 25 000 published measurements enable the calculation of 15 parameters in the same worker only by using the VO₂ reflecting workload. These parameters, notably VE, VA, VE/VO₂ VA/Q, Q, LBF and KBF were found to be interrelated. Altering one of these changes the order of magnitude of the others. Q, LBF and KBF decrease when supine adults at rest switch to an upright position. This effect of gravity diminished when VO₂ increased. The fall in LBF and KBF during exertion might enhance muscle blood flow as reported previously. Taken together these equations and data may improve the accuracy of physiologically based toxicokinetic modeling as well as occupational health assessment studies in active workers exposed to xenobiotics.List of main abbreviations: AVOD: arterioveinous oxygen content difference.BMI: body mass index (in kg/m²).BSA: body surface area (in m²).BTPS: body temperature and saturated with water vapor.Bw: body weight (in kg).E: minute energy expenditure rate (in kcal/min).F_GE: organ blood flow factor for the gravitational effect on blood circulation.H: oxygen uptake factor, volume of oxygen (at STPD) consumed to produce 1 kcal of energy expended.KBF: kidneys blood flow (in ml/min).LBF: liver blood flow (in ml/min).P_BF: liver or kidneys blood flows expressed in terms of percentages (in %) of Q_{sup C} values: namely P_BF = (LBF or KBF/Q_{sup C}) x 100.Q: cardiac output (in L/min or ml/min).Q_{sup C:} cardiac output for the cohort of males or females in supination (in ml/min).STPD: standard temperature and pressure, dry air.sup: values measured when adults are in the supine position.up: values measured when adults are in the upright position.VDphys: physiological dead space at BTPS (in L).VT: tidal volume at BTPS (in L).VA: alveolar ventilation rate at BTPS (in L/min).VA/Q: ventilation-perfusion ratio (unitless).VE: minute ventilation rate at BTPS (in L/min).VO₂: oxygen consumption rate (i.e. the oxygen uptake) at STPD (in L/min).VQ: ventilatory equivalent for VO₂ (VE at BTPS /VO₂ at STPD).

Contamination of drinking water by neonicotinoid insecticides in China: Human exposure potential through drinking water consumption and percutaneous penetration

Neonicotinoids (NEOs) are the most widely used pesticides and have posed a serious threat to human health. However, data on human exposure to NEOs are extremely scarce. To bridge this gap, human exposure potential of NEOs through drinking water consumption and percutaneous penetration was evaluated with the influences of 17 age groups, 4 seasons, 6 regions, and 2 genders. The results showed that drinking water in the present study had an upper middle level of NEO contamination. Anthropogenic activity and weather condition played important roles in the regional distribution of NEOs in tap water. For both children and adults, NEOs intake from drinking water exposure (NDE) and percutaneous exposure (NPE) in the south regions of China are significantly higher than those in the north regions, while the order of NDE and NPE by season is summer > spring = autumn > winter. Furthermore, human age and gender also have remarkable impacts on NDE and NPE. The age groups of children subjected to the highest NDE and NPE were 9 months - 2 years old and 9-12 years old, respectively. This study provides insights into the role of seasonal and regional influence, age and gender in the risk of drinking water and percutaneous exposure to NEOs.

High-throughput PBTK models for in vitro to in vivo extrapolation

INTRODUCTION

Toxicity data are unavailable for many thousands of chemicals in commerce and the environment. Therefore, risk assessors need to rapidly screen these chemicals for potential risk to public health. High-throughput screening (HTS) for in vitro bioactivity, when used with high-throughput toxicokinetic (HTTK) data and models, allows characterization of these thousands of chemicals.

AREAS COVERED

This review covers generic physiologically based toxicokinetic (PBTK) models and high-throughput PBTK modeling for in vitro-in vivo extrapolation (IVIVE) of HTS data. We focus on 'httk', a public, open-source set of computational modeling tools and in vitro toxicokinetic (TK) data.

EXPERT OPINION

HTTK benefits chemical risk assessors with its ability to support rapid chemical screening/prioritization, perform IVIVE, and provide provisional TK modeling for large numbers of chemicals using only limited chemical-specific data. Although generic TK model design can increase prediction uncertainty, these models provide offsetting benefits by increasing model implementation accuracy. Also, public distribution of the models and data enhances reproducibility. For the httk package, the modular and open-source design can enable the tool to be used and continuously improved by a broad user community in support of the critical need for high-throughput chemical prioritization and rapid dose estimation to facilitate rapid hazard assessments.

Exposure to disinfection byproducts and risk of type 2 diabetes: a nested case-control study in the HUNT and Lifelines cohorts

INTRODUCTION

Environmental chemicals acting as metabolic disruptors have been implicated with diabetogenesis, but evidence is weak among short-lived chemicals, such as disinfection byproducts (trihalomethanes, THM composed of chloroform, TCM and brominated trihalomethanes, BrTHM).

OBJECTIVES

We assessed whether THM were associated with type 2 diabetes (T2D) and we explored alterations in metabolic profiles due to THM exposures or T2D status.

METHODS

A prospective 1:1 matched case-control study (n = 430) and a cross-sectional 1:1 matched case-control study (n = 362) nested within the HUNT cohort (Norway) and the Lifelines cohort (Netherlands), respectively, were set up. Urinary biomarkers of THM exposure and mass spectrometry-based serum metabolomics were measured. Associations between THM, clinical markers, metabolites and disease status were evaluated using logistic regressions with Least Absolute Shrinkage and Selection Operator procedure.

RESULTS

Low median THM exposures (ng/g, IQR) were measured in both cohorts (cases and controls of HUNT and Lifelines, respectively, 193 (76, 470), 208 (77, 502) and 292 (162, 595), 342 (180, 602). Neither BrTHM (OR = 0.87; 95% CI: 0.67, 1.11 | OR = 1.09; 95% CI: 0.73, 1.61), nor TCM (OR = 1.03; 95% CI: 0.88, 1.2 | OR = 1.03; 95% CI: 0.79, 1.35) were associated with incident or prevalent T2D, respectively. Metabolomics showed 48 metabolites associated with incident T2D after adjusting for sex, age and BMI, whereas a total of 244 metabolites were associated with prevalent T2D. A total of 34 metabolites were associated with the progression of T2D. In data driven logistic regression, novel biomarkers, such as cinnamoylglycine or 1-methylurate, being protective of T2D were identified. The incident T2D risk prediction model (HUNT) predicted well incident Lifelines cases (AUC = 0.845; 95% CI: 0.72, 0.97).

CONCLUSION

Such exposome-based approaches in cohort-nested studies are warranted to better understand the environmental origins of diabetogenesis.

The Impact of Scaling Factor Variability on Risk-Relevant Pharmacokinetic Outcomes in Children: A Case Study Using Bromodichloromethane (BDCM)

Biotransformation rates extrapolated from in vitro data are used increasingly in human physiologically based pharmacokinetic (PBPK) models. This practice requires use of scaling factors, including microsomal content (mg of microsomal protein/g liver, MPPGL), enzyme specific content, and liver mass as a fraction of body weight (FVL). Previous analyses indicated that scaling factor variability impacts pharmacokinetic (PK) outcomes used in adult population dose-response studies. This analysis was extended to pediatric populations because large inter-individual differences in enzyme ontogeny likely would further contribute to scaling factor variability. An adult bromodichloromethane (BDCM) model (Kenyon, E. M., Eklund, C., Leavens, T. L., and Pegram, R. A. (2016a). Development and application of a human PBPK model for bromodichloromethane (BDCM) to investigate impacts of multi-route exposure. J. Appl. Toxicol. 36, 1095-1111) was re-parameterized for neonates, infants, and toddlers. Monte Carlo analysis was used to assess the impact of pediatric scaling factor variation on model-derived PK outcomes compared with adult findings. BDCM dose metrics were estimated following a single 0.05-liter drink of water or a 20-min bath, under typical (5 µg/l) and plausible higher (20 µg/l) BDCM concentrations. MPPGL, CYP2E1, and FVL values reflected the distribution of reported pediatric population values. The impact of scaling factor variability on PK outcome variation was different for each exposure scenario, but similar for each BDCM water concentration. The higher CYP2E1 expression variability during early childhood was reflected in greater variability in predicted PK outcomes in younger age groups, particularly for the oral exposure route. Sensitivity analysis confirmed the most influential parameter for this variability was CYP2E1, particularly in neonates. These findings demonstrate the importance of age-dependent scaling factor variation used for in vitro to in vivo extrapolation of biotransformation rates.

Health risk assessment of trihalomethanes mixtures from daily water-related activities via multi-pathway exposure based on PBPK model

In this study, the concentrations of trihalomethanes (THMs) in tap water and direct drinking water were analyzed, and based on the human behavior patterns and building parameters, the concentrations of THMs in indoor air were simulated with the water-air concentration conversion model. In addition, concentrations of THMs in human tissues were predicted based on physiologically based pharmacokinetic (PBPK) model, and the health risk of THMs for participants were estimated. Furthermore, the carcinogenic risk of mixtures according to the method proposed by USEPA and PBPK model based method were calculated and compared. The concentrations of chloroform, bromodichloromethane, dibromochloromethane and bromoform in tap water were 11.28-16.21, 4.83-6.28, 0.81-1.32 and 0.08-0.21 μg/L, and those in direct drinking water were 3.29-6.88, 0.35-0.47, 0.03-0.08 and 0.04-0.08 μg/L, respectively. The results of water-air concentration conversion model demonstrated that pollutants in air had a strong correlation with water-related activities. Multi-pathway PBPK model showed that THMs concentrations in liver, kidney and richly perfused tissue were higher than those in other tissues. The results of risk assessment showed that the mean risk levels of mixtures were 1.69 × 10^-4 and 1.72 × 10^-4 calculated by the USEPA recommended method and PBPK based method, which seriously exceeded the acceptable level. TCM and BDCM were the major risk factors, and inhalation was the main exposure route of THMs.

Formation of trihalomethanes in swimming pool waters using sodium dichloroisocyanurate as an alternative disinfectant

Disinfection is essential to preventing infection caused by microbial pathogens in swimming pool water. The most commonly used disinfection methods are chlorine disinfectant, including sodium hypochlorite (NaOCl) and sodium dichloroisocyanurate (SDIC, C₃O₃N₃Cl₂Na) with characteristics of available chlorine formation and stability. In this study, we estimate the formation of trihalomethanes (THM₄) in indoor swimming pools filled with seawater that adopt these disinfection methods, and we investigated the factors influencing the formation of THM₄. Formation of THM₄ by free chlorine (FC) and SDIC respectively is 327.8 μg L^-1 and 307.6 μg L^-1; Br-THMs is 226.7 μg L^-1 for FC, 198 μg L^-1 for SDIC. SDIC has less THM₄ formation than FC in the same molar Br^- dosage (10 μM) and total chlorine dosage (20 μM), and bromoform is the main Br-THMs species. The occurrence of THM₄ is inhibited at high dosages of natural organic matter and Br^-. The total Br-THMs increase from 75 μg L^-1 to 189.7 μg L^-1 and from 64.6 μg L^-1 to 190 μg L^-1 by FC and SDIC at pH 5.0-9.0, both of which are highly dependent on pH. In real water, similar results were found in Br^--containing water (1 mg-Br^- L^-1).

Reverse dosimetry modeling of toluene exposure concentrations based on biomonitoring levels from the Canadian health measures survey

Biomonitoring might provide useful estimates of population exposure to environmental chemicals. However, data uncertainties stemming from interindividual variability are common in large population biomonitoring surveys. Physiologically based pharmacokinetic (PBPK) models might be used to account for age- and gender-related variability in internal dose. The objective of this study was to reconstruct air concentrations consistent with blood toluene measures reported in the third Canadian Health Measures Survey using reverse dosimetry PBPK modeling techniques. Population distributions of model's physiological parameters were described based upon age, weight, and size for four subpopulations (12-19, 20-39, 40-59, and 60-79 years old). Monte Carlo simulations applied to PBPK modeling allowed converting the distributions of venous blood measures of toluene obtained from CHMS into related air levels. Based upon blood levels observed at the 50^th, 90^th and 95^th percentiles, corresponding air toluene concentrations were estimated for teenagers aged 12-19 years as being, respectively, 0.009, 0.04 and 0.06 ppm. Similarly, values were computed for adults aged 20-39 years (0.007, 0.036, and 0.06 ppm), 40-59 years (0.007, 0.036 and 0.06 ppm) and 60-79 years (0.006, 0.022 and 0.04 ppm). These estimations are well below Health Canada's maximum recommended chronic air guidelines for toluene. In conclusion, PBPK modeling and reverse dosimetry may be combined to help interpret biomonitoring data for chemical exposure in large population surveys and estimate the associated toxicological health risk.

A Physiologically-Based Pharmacokinetic Modeling Approach Using Biomonitoring Data in Order to Assess the Contribution of Drinking Water for the Achievement of an Optimal Fluoride Dose for Dental Health in Children

Due to an optimal fluoride concentration in drinking water advised for caries prevention purposes, the population is now exposed to multiple sources of fluoride. The availability of population biomonitoring data currently allow us to evaluate the magnitude of this exposure. The objective of this work was, therefore, to use such data in order to estimate whether community water fluoridation still represents a significant contribution toward achieving a suggested daily optimal fluoride (external) intake of 0.05 mg/kg/day. Therefore, a physiologically-based pharmacokinetic model for fluoride published in the literature was used and adapted in Excel for a typical 4-year-old and 8-year-old child. Biomonitoring data from the Canadian Health Measures Survey among people living in provinces with very different drinking water fluoridation coverage (Quebec, 2.5%; Ontario, 70% of the population) were analyzed using this adapted model. Absorbed doses for the 4-year-old and 8-year-old children were, respectively, 0.03 mg/kg/day and 0.02 mg/kg/day in Quebec and of 0.06 mg/kg/day and 0.05 mg/kg/day in Ontario. These results show that community water fluoridation contributes to increased fluoride intake among children, which leads to reaching, and in some cases even exceeding, the suggested optimal absorbed dose of 0.04 mg/kg/day, which corresponds to the suggested optimal fluoride intake mentioned above. In conclusion, this study constitutes an incentive to further explore the multiple sources of fluoride intake and suggests that a new balance between them including drinking water should be examined in accordance with the age-related physiological differences that influence fluoride metabolism.

Time of the day dictates the variability of biomarkers of exposure to disinfection byproducts

Non-persistent environmental chemicals (NOPEC) are xenobiotics with short half-lives of elimination (<7h). Similar to chronopharmacokinetics, NOPEC metabolism may follow diurnal patterns of cytochrome P450 activity. The role of circadian liver clock in shaping NOPEC metabolism and their concomitant measurements of biomarkers of exposure and effect remains poorly understood in real-life human settings. Metabolic activation (toxication) by CYP2E1 converts trihalomethanes (THM) to harmful metabolites. We investigated the diurnal variation of urinary THM exposures and their metabolism patterns as catalyzed by CYP2E1 redox activity, using the surrogate marker of 4-hydroxynonenal (4HNE). We implemented three time-series trials with adult volunteers conducting specific household cleaning activities at predefined times of a day. Circadia variation of 4HNE was assessed with a cosinor model and its mesor levels increased with THM exposure. The time of exposure within the day dictated the magnitude of urinary THM levels and their toxication effect; in all three trials and relative to urinary THM levels before the activity, lower and higher median THM were measured right after the activity in morning and afternoon/night, respectively. This is consistent with higher reported CYP2E1 redox activity in light/active phase. Population health studies should incorporate time-stamped biomarker data to improve the understanding of chronic disease processes.

A framework and case studies for evaluation of enzyme ontogeny in children's health risk evaluation

Knowledge of the ontogeny of Phase I and Phase II metabolizing enzymes may be used to inform children's vulnerability based upon likely differences in internal dose from xenobiotic exposure. This might provide a qualitative assessment of toxicokinetic (TK) variability and uncertainty pertinent to early lifestages and help scope a more quantitative physiologically based toxicokinetic (PBTK) assessment. Although much is known regarding the ontogeny of metabolizing systems, this is not commonly utilized in scoping and problem formulation stage of human health risk evaluation. A framework is proposed for introducing this information into problem formulation which combines data on enzyme ontogeny and chemical-specific TK to explore potential child/adult differences in internal dose and whether such metabolic differences may be important factors in risk evaluation. The framework is illustrated with five case study chemicals, including some which are data rich and provide proof of concept, while others are data poor. Case studies for toluene and chlorpyrifos indicate potentially important child/adult TK differences while scoping for acetaminophen suggests enzyme ontogeny is unlikely to increase early-life risks. Scoping for trichloroethylene and aromatic amines indicates numerous ways that enzyme ontogeny may affect internal dose which necessitates further evaluation. PBTK modeling is a critical and feasible next step to further evaluate child-adult differences in internal dose for a number of these chemicals.

Prenatal exposure to drinking-water chlorination by-products, cytochrome P450 gene polymorphisms and small-for-gestational-age neonates

Genetic susceptibility may modulate chlorination by-products (CBPs) effects on fetal growth, especially genes coding for the cytochrome P450 involved in the metabolism of CBPs and steroidogenesis. In a case-control study of 1432 mother-child pairs, we assessed the association between maternal and child single nucleotide polymorphisms (SNPs) within CYP1A2, CYP2A6, CYP2D6 and CYP17A1 genes and small-for-gestational-age neonates (SGA<10th percentile) as well as interaction between these SNPs and maternal exposure to trihalomethanes or haloacetic acids (HAAs) during the third trimester of pregnancy. Interactions were found between mother and neonate carrying CYP17A1 rs4919687A and rs743572G alleles and maternal exposure to total trihalomethanes or five regulated HAAs species. However, these interactions became non statistically significant after correction for multiple testing. There is some evidence, albeit weak, of a potential effect modification of the association between CBPs and SGA by SNPs in CYP17A1 gene. Further studies are needed to validate these observations.

A decision tree approach to screen drinking water contaminants for multiroute exposure potential in developing guideline values

The consideration of inhalation and dermal routes of exposures in developing guideline values for drinking water contaminants is important. However, there is no guidance for determining the eligibility of a drinking water contaminant for its multiroute exposure potential. The objective of the present study was to develop a 4-step framework to screen chemicals for their dermal and inhalation exposure potential in the process of developing guideline values. The proposed framework emphasizes the importance of considering basic physicochemical properties prior to detailed assessment of dermal and inhalation routes of exposure to drinking water contaminants in setting guideline values.

Development and application of a human PBPK model for bromodichloromethane to investigate the impacts of multi-route exposure

As a result of its presence in water as a volatile disinfection byproduct, bromodichloromethane (BDCM), which is mutagenic, poses a potential health risk from exposure via oral, dermal and inhalation routes. We developed a refined human physiologically based pharmacokinetic (PBPK) model for BDCM (including new chemical-specific human parameters) to evaluate the impact of BDCM exposure during showering and bathing on important measures of internal dose compared with oral exposure. The refined model adequately predicted data from the published literature for oral, dermal and bathing/showering exposures. A liter equivalency approach (L-eq) was used to estimate BDCM concentration in a liter of water consumed by the oral route that would be required to produce the same internal dose of BDCM resulting from a 20-min bath or a 10-min shower in water containing 10 µg l(-1) BDCM. The oral liter equivalent concentrations for the bathing scenario were 605, 803 and 5 µg l(-1) BDCM for maximum venous blood concentration (Cmax), the area under the curve (AUCv) and the amount metabolized in the liver per hour (MBDCM), respectively. For a 10-min showering exposure, the oral L-eq concentrations were 282, 312 and 2.1 µg l(-1) for Cmax, AUC and MBDCM, respectively. These results demonstrate large contributions of dermal and inhalation exposure routes to the internal dose of parent chemical reaching the systemic circulation, which could be transformed to mutagenic metabolites in extrahepatic target tissues. Thus, consideration of the contribution of multiple routes of exposure when evaluating risks from water-borne BDCM is needed, and this refined human model will facilitate improved assessment of internal doses from real-world exposures. Published 2015. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Assessing human variability in kinetics for exposures to multiple environmental chemicals: a physiologically based pharmacokinetic modeling case study with dichloromethane, benzene, toluene, ethylbenzene, and m-xylene

The objective of this study was to compare the magnitude of interindividual variability in internal dose for inhalation exposure to single versus multiple chemicals. Physiologically based pharmacokinetic models for adults (AD), neonates (NEO), toddlers (TODD), and pregnant women (PW) were used to simulate inhalation exposure to "low" (RfC-like) or "high" (AEGL-like) air concentrations of benzene (Bz) or dichloromethane (DCM), along with various levels of toluene alone or toluene with ethylbenzene and xylene. Monte Carlo simulations were performed and distributions of relevant internal dose metrics of either Bz or DCM were computed. Area under the blood concentration of parent compound versus time curve (AUC)-based variability in AD, TODD, and PW rose for Bz when concomitant "low" exposure to mixtures of increasing complexities occurred (coefficient of variation (CV) = 16-24%, vs. 12-15% for Bz alone), but remained unchanged considering DCM. Conversely, AUC-based CV in NEO fell (15 to 5% for Bz; 12 to 6% for DCM). Comparable trends were observed considering production of metabolites (AMET), except for NEO's CYP2E1-mediated metabolites of Bz, where an increased CV was observed (20 to 71%). For "high" exposure scenarios, Cmax-based variability of Bz and DCM remained unchanged in AD and PW, but decreased in NEO (CV= 11-16% to 2-6%) and TODD (CV= 12-13% to 7-9%). Conversely, AMET-based variability for both substrates rose in every subpopulation. This study analyzed for the first time the impact of multiple exposures on interindividual variability in toxicokinetics. Evidence indicates that this impact depends upon chemical concentrations and biochemical properties, as well as the subpopulation and internal dose metrics considered.

Use of urinary trichloroacetic acid as an exposure biomarker of disinfection by-products in cancer studies

Urinary trichloroacetic acid (TCAA) has been proposed as a valid exposure biomarker for ingested disinfection by-products (DBP) for reproductive studies. However, it has never been used in epidemiologic studies on cancer. We investigate the performance of urinary TCAA as a biomarker of DBP exposure in the framework of an epidemiologic study on cancer. We conducted home visits to collect tap water, first morning void urine, and a 48h fluid intake diary among 120 controls from a case-control study of colorectal cancer in Barcelona, Spain. We measured urine TCAA and creatinine, and 9 haloacetic acids and 4 trihalomethanes (THM) in tap water. Lifetime THM exposure was estimated based on residential history since age 18 plus routine monitoring data. Robust linear regressions were used to estimate mean change in urinary TCAA adjusted by covariates. Among the studied group, mean age was 74 years (range 63-85) and 41 (34%) were females. Mean total tap water consumption was 2.2l/48h (standard error, 0.1l/48h). Geometric mean urine TCAA excretion rate was 17.3pmol/min [95%CI: 14.0-21.3], which increased 2% for a 10% increase in TCAA ingestion and decreased with total tap water consumption (-17%/l), water intake outside home (-32%), plasmatic volume (-64%/l), in smokers (-79%), and in users of non-steroidal anti-inflammatory drugs (-50%). Urinary TCAA levels were not associated with lifetime THM exposure. In conclusion, our findings support that urine TCAA is not a valid biomarker in case-control studies of adult cancer given that advanced age, comorbidites and medication use are prevalent and are determinants of urine TCAA levels, apart from ingested TCAA levels. In addition, low TCAA concentrations in drinking water limit the validity of urine TCAA as an exposure biomarker.

Spatial and seasonal variability of urinary trihalomethanes concentrations in urban settings

A complex network of sources and routes of exposure to disinfection by-products (DBP), such as trihalomethanes (THM) has been driving the wide variability of daily THM intake estimates in environmental epidemiological studies. We hypothesized that the spatiotemporal variability of THM exposures could be differentially expressed with their urinary levels among residents whose households are geographically clustered in district-metered areas (DMA) receiving the same tap water. Each DMA holds unique drinking-water pipe network characteristics, such as pipe length, number of pipe leaking incidences, number of water meters by district, average minimum night flow and average daily demand. The present study assessed the spatial and seasonal variability in urinary THM levels among residents (n=310) of geocoded households belonging to two urban DMA of Nicosia, Cyprus, with contrasting water network properties. First morning urine voids were collected once in summer and then in winter. Results showed that the mean sum of the four urinary THM analytes (TTHM) was significantly higher during summer for residents of both areas. Linear mixed effects models adjusted for age, season and gender, illustrated spatially-resolved differences in creatinine-adjusted urinary chloroform and TTHM levels between the two studied areas, corroborated by differences observed in their pipe network characteristics. Additional research is warranted to shed light on the contribution of spatially-resolved and geographically-clustered environmental exposures coupled with internal biomarker of exposure measurements towards better understanding of health disparities within urban centers.

Characterization of the human kinetic adjustment factor for the health risk assessment of environmental contaminants

A default uncertainty factor of 3.16 (√10) is applied to account for interindividual variability in toxicokinetics when performing non-cancer risk assessments. Using relevant human data for specific chemicals, as WHO/IPCS suggests, it is possible to evaluate, and replace when appropriate, this default factor by quantifying chemical-specific adjustment factors for interindividual variability in toxicokinetics (also referred to as the human kinetic adjustment factor, HKAF). The HKAF has been determined based on the distributions of pharmacokinetic parameters (e.g., half-life, area under the curve, maximum blood concentration) in relevant populations. This article focuses on the current state of knowledge of the use of physiologically based algorithms and models in characterizing the HKAF for environmental contaminants. The recent modeling efforts on the computation of HKAF as a function of the characteristics of the population, chemical and its mode of action (dose metrics), as well as exposure scenario of relevance to the assessment are reviewed here. The results of these studies, taken together, suggest the HKAF varies as a function of the sensitive subpopulation and dose metrics of interest, exposure conditions considered (route, duration, and intensity), metabolic pathways involved and theoretical model underlying its computation. The HKAF seldom exceeded the default value of 3.16, except in very young children (i.e., <≈ 3 months) and when the parent compound is the toxic moiety. Overall, from a public health perspective, the current state of knowledge generally suggest that the default uncertainty factor is sufficient to account for human variability in non-cancer risk assessments of environmental contaminants.

Accounting for the impact of short-term variations in the levels of trihalomethane in drinking water on exposure assessment for epidemiological purposes. Part II: biological aspects

The variability of trihalomethane (THM) levels in drinking water raises the question of whether or not short-term variations (within-day) should be accounted for when assessing exposure to contaminants suspected of being carcinogenic and reprotoxic agents. The purpose of this study was to determine the magnitude of the impact on predicted biological levels of THMs (internal doses) exerted by within-day variations of THMs in drinking water. A database extracted from a campaign in the Québec City distribution system served to produce 81, 79 and 64 concentration profiles for the three most abundant THMs, namely chloroform (TCM), dichlorobromomethane (DCBM) and chlorodibromomethane (CDBM), respectively. Using a physiologically based toxicokinetic modeling approach, we simulated exposures (1.5 l water per day and a 10-min shower) based on each of these profiles and predicted, for 2000 individuals (Monte-Carlo simulations), maximum blood concentrations (Cmax), areas under the time versus blood concentrations curve (24 h-AUCcv) and total absorbed doses (ADs). Three different hypotheses were tested: [A] assuming a constant THM concentration in water (e.g., mean value of a day); [B] accounting for within-day variations in THM levels; and [C] a worst-case scenario assuming within-day variations and showering while THM levels were maximal. For each exposure profile, exposure indicator and individual, we calculated the ratios of values obtained according to each hypothesis (e.g., CmaxB/CmaxA and CmaxC/CmaxA) and the values corresponding to the 5th and 95th percentiles of these ratios. The closer these percentiles are to the value of 1, the smaller the error associated with assuming constant THM concentrations rather than their actual variability. Results showed that the minimal gap between these percentiles was TCM-AD(B)/TCM-AD(A) (5th=0.91; 95th=1.09), whereas the maximal gap was CDBM-Cmax(C)/CDBM-Cmax(A) (5th=0.50; 95th=3.40). Overall, TCM and ADs were the less affected (TCM<DCBM<CDBM and AD<AUCcv<Cmax) when accounting for within-day variations in water levels.

Assessing the impact of child/adult differences in hepatic first-pass effect on the human kinetic adjustment factor for ingested toxicants

The objective of this study was to evaluate the impact of interindividual differences in hepatic first-pass effect (FPE) on the magnitude of the human kinetic adjustment factor (HKAF) for ingested toxicants. This factor aims at replacing a default value of 3.2 used in non-cancer risk assessment. Coupled with Monte Carlo simulations, steady-state equations that account for FPE were used to obtain distributions of arterial blood concentrations (CAss) and rates of metabolism in adults, neonates, infants and toddlers continuously exposed to an oral dose of 1 μg/kg/d of theoretical CYP2E1 and CYP1A2 substrates. For such substrates exhibiting a range of blood:air partition coefficients (Pb: 1-10,000) and hepatic extraction ratios in an average adult (E(ad): 0.01-0.99), HKAFs were computed as the ratio of the 95th percentile of dose metrics for each subpopulation over the 50th percentile value in adults. The reduced hepatic enzyme content in neonates as compared to adults resulted in correspondingly diminished FPE. Consequently, HKAFs greater than 3.2 could be observed, based on CAss only, in the following cases: for some CYP2E1 substrates with E(ad) ≤ 0.3, in neonates (max.: 6.3); and for some CYP1A2 substrates with E(ad) ≤ 0.1 and 0.7, in, respectively, neonates and infants (max.: 28.3). Overall, this study pointed out the importance of accounting for child/adult differences in FPE when determining the HKAF for oral exposure.

Occurrence and spatial and temporal variations of disinfection by-products in the water and air of two indoor swimming pools

In order to improve disinfection by-product (DBP) exposure assessment, this study was designed to document both water and air levels of these chemical contaminants in two indoor swimming pools and to analyze their within-day and day-to-day variations in both of them. Intensive sampling was carried out during two one-week campaigns to measure trihalomethanes (THMs) and chloramines (CAMs) in water and air, and haloacetic acids (HAAs) in water several times daily. Water samples were systematically collected at three locations in each pool and air samples were collected at various heights around the pool and in other rooms (e.g., changing room) in the buildings. In addition, the ability of various models to predict air concentrations from water was tested using this database. No clear trends, but actual variations of contamination levels, appeared for both water and air according to the sampling locations and times. Likewise, the available models resulted in realistic but imprecise estimates of air contamination levels from water. This study supports the recommendation that suitable minimal air and water sampling should be carried out in swimming pools to assess exposure to DBPs.

Maternal exposure to drinking-water chlorination by-products and small-for-gestational-age neonates

BACKGROUND

There is concern about possible effects of disinfection by-products on reproductive outcomes. The purpose of this study was to evaluate the association between maternal exposure to chlorination by-products and the risk of delivering a small for-gestational-age (SGA) neonate.

METHODS

We conducted a population-based case-control study in the Québec City (Canada) area. Term newborn cases with birth weights <10th percentile (n = 571) were compared with 1925 term controls with birth weights ≥10th percentile. Concentrations of trihalomethanes and haloacetic acids in the water-distribution systems of participants were monitored during the study period, and a phone interview on maternal habits was completed within 3 months after childbirth. We estimated chlorination by-products ingestion during the last trimester of pregnancy and trihalomethanes doses resulting from inhalation and dermal exposure. We evaluated associations between chlorination by-products in utero exposure and SGA by means of unconditional logistic regression with control of potential confounders.

RESULTS

When total trihalomethanes and the 5 regulated haloacetic acids concentrations were divided into quartiles, no clear dose-response relationship was found with SGA. However, increased risk was observed when haloacetic concentrations were above the fourth quartile and when either trihalomethanes or haloacetic acids concentrations were above current water standards (adjusted OR= 1.5 [95% confidence interval = 1.1-1.9] and 1.4 [1.1-1.9], respectively). Inhalation and dermal absorption of trihalomethanes did not contribute to this risk, but a monotonic dose-response was found with haloacetic acids ingestion.

CONCLUSION

Oral exposure to high levels of chlorination by-products in drinking water could be a risk factor for term SGA.

Quantitative Property-Property Relationship for Screening-Level Prediction of Intrinsic Clearance of Volatile Organic Chemicals in Rats and Its Integration within PBPK Models to Predict Inhalation Pharmacokinetics in Humans

The objectives of this study were (i) to develop a screening-level Quantitative property-property relationship (QPPR) for intrinsic clearance (CL_int) obtained from in vivo animal studies and (ii) to incorporate it with human physiology in a PBPK model for predicting the inhalation pharmacokinetics of VOCs. CL_int, calculated as the ratio of the in vivo V_max (μmol/h/kg bw rat) to the K_m (μM), was obtained for 26 VOCs from the literature. The QPPR model resulting from stepwise linear regression analysis passed the validation step (R² = 0.8; leave-one-out cross-validation Q² = 0.75) for CL_int normalized to the phospholipid (PL) affinity of the VOCs. The QPPR facilitated the calculation of CL_int (L PL/h/kg bw rat) from the input data on log P_ow, log blood: water PC and ionization potential. The predictions of the QPPR as lower and upper bounds of the 95% mean confidence intervals (LMCI and UMCI, resp.) were then integrated within a human PBPK model. The ratio of the maximum (using LMCI for CL_int) to minimum (using UMCI for CL_int) AUC predicted by the QPPR-PBPK model was 1.36 ± 0.4 and ranged from 1.06 (1,1-dichloroethylene) to 2.8 (isoprene). Overall, the integrated QPPR-PBPK modeling method developed in this study is a pragmatic way of characterizing the impact of the lack of knowledge of CL_int in predicting human pharmacokinetics of VOCs, as well as the impact of prediction uncertainty of CL_int on human pharmacokinetics of VOCs.

Comparison of trihalomethanes in tap water and blood: a case study in the United States

BACKGROUND

Epidemiological studies have used various measures to characterize trihalomethane (THM) exposures, but the relationship of these indicators to exposure biomarkers remains unclear.

OBJECTIVES

We examined temporal and spatial variability in baseline blood THM concentrations and assessed the relationship between these concentrations and several exposure indicators (tap water concentration, water-use activities, multiroute exposure metrics).

METHODS

We measured water-use activity and THM concentrations in blood and residential tap water from 150 postpartum women from three U.S. locations.

RESULTS

Blood ΣTHM [sum of chloroform (TCM), bromodichloromethane (BDCM), dibromo-chloromethane (DBCM), and bromoform (TBM)] concentrations varied by site and season. As expected based on variable tap water concentrations and toxicokinetic properties, the proportion of brominated species (BDCM, DBCM, and TBM) in blood varied by site (site 1, 24%; site 2, 29%; site 3, 57%) but varied less markedly than in tap water (site 1, 35%; site 2, 75%; site 3, 68%). The blood-water ΣTHM Spearman rank correlation coefficient was 0.36, with correlations higher for individual brominated species (BDCM, 0.62; DBCM, 0.53; TBM, 0.54) than for TCM (0.37). Noningestion water activities contributed more to the total exposure metric than did ingestion, but tap water THM concentrations were more predictive of blood THM levels than were metrics that incorporated water use.

CONCLUSIONS

Spatial and temporal variability in THM concentrations was greater in water than in blood. We found consistent blood-water correlations across season and site for BDCM and DBCM, and multivariate regression results suggest that water THM concentrations may be an adequate surro-gate for baseline blood levels.

Adverse birth outcomes and maternal exposure to trichloroethylene and tetrachloroethylene through soil vapor intrusion in New York State

BACKGROUND

Industrial spills of volatile organic compounds (VOCs) in Endicott, New York (USA), have led to contamination of groundwater, soil, and soil gas. Previous studies have reported an increase in adverse birth outcomes among women exposed to VOCs in drinking water.

OBJECTIVE

We investigated the prevalence of adverse birth outcomes among mothers exposed to trichloroethylene (TCE) and tetrachloroethylene [or perchloroethylene (PCE)] in indoor air contaminated through soil vapor intrusion.

METHODS

We examined low birth weight (LBW), preterm birth, fetal growth restriction, and birth defects among births to women in Endicott who were exposed to VOCs, compared with births statewide. We used Poisson regression to analyze births and malformations to estimate the association between maternal exposure to VOCs adjusting for sex, mother's age, race, education, parity, and prenatal care. Two exposure areas were identified based on environmental sampling data: one area was primarily contaminated with TCE, and the other with PCE.

RESULTS

In the TCE-contaminated area, adjusted rate ratios (RRs) were significantly elevated for LBW [RR = 1.36; 95% confidence interval (CI): 1.07, 1.73; n = 76], small for gestational age (RR = 1.23; 95% CI: 1.03, 1.48; n = 117), term LBW (RR = 1.68; 95% CI: 1.20, 2.34; n = 37), cardiac defects (RR = 2.15; 95% CI: 1.27, 3.62; n = 15), and conotruncal defects (RR = 4.91; 95% CI: 1.58, 15.24; n = 3). In the PCE-contaminated area, RRs for cardiac defects (five births) were elevated but not significantly. Residual socioeconomic confounding may have contributed to elevations of LBW outcomes.

CONCLUSIONS

Maternal residence in both areas was associated with cardiac defects. Residence in the TCE area, but not the PCE area, was associated with LBW and fetal growth restriction.

An assessment of the interindividual variability of internal dosimetry during multi-route exposure to drinking water contaminants

The objective of this study was to evaluate inter-individual variability in absorbed and internal doses after multi-route exposure to drinking water contaminants (DWC) in addition to the corresponding variability in equivalent volumes of ingested water, expressed as liter-equivalents (LEQ). A multi-route PBPK model described previously was used for computing the internal dose metrics in adults, neonates, children, the elderly and pregnant women following a multi-route exposure scenario to chloroform and to tri- and tetra-chloroethylene (TCE and PERC). This scenario included water ingestion as well as inhalation and dermal contact during a 30-min bathroom exposure. Monte Carlo simulations were performed and distributions of internal dose metrics were obtained. The ratio of each of the dose metrics for inhalation, dermal and multi-route exposures to the corresponding dose metrics for the ingestion of drinking water alone allowed computation of LEQ values. Mean BW-adjusted LEQ values based on absorbed doses were greater in neonates regardless of the contaminant considered (0.129–0.134 L/kg BW), but higher absolute LEQ values were obtained in average adults (3.6–4.1 L), elderly (3.7–4.2 L) and PW (4.1–5.6 L). LEQ values based on the parent compound’s AUC were much greater than based on the absorbed dose, while the opposite was true based on metabolite-based dose metrics for chloroform and TCE, but not PERC. The consideration of the 95th percentile values of BW-adjusted LEQ did not significantly change the results suggesting a generally low intra-subpopulation variability during multi-route exposure. Overall, this study pointed out the dependency of the LEQ on the dose metrics, with consideration of both the subpopulation and DWC.

What's in the pool? A comprehensive identification of disinfection by-products and assessment of mutagenicity of chlorinated and brominated swimming pool water

BACKGROUND

Swimming pool disinfectants and disinfection by-products (DBPs) have been linked to human health effects, including asthma and bladder cancer, but no studies have provided a comprehensive identification of DBPs in the water and related that to mutagenicity.

OBJECTIVES

We performed a comprehensive identification of DBPs and disinfectant species in waters from public swimming pools in Barcelona, Catalonia, Spain, that disinfect with either chlorine or bromine and we determined the mutagenicity of the waters to compare with the analytical results.

METHODS

We used gas chromatography/mass spectrometry (GC/MS) to measure trihalomethanes in water, GC with electron capture detection for air, low- and high-resolution GC/MS to comprehensively identify DBPs, photometry to measure disinfectant species (free chlorine, monochloroamine, dichloramine, and trichloramine) in the waters, and an ion chromatography method to measure trichloramine in air. We assessed mutagenicity with the Salmonella mutagenicity assay.

RESULTS

We identified > 100 DBPs, including many nitrogen-containing DBPs that were likely formed from nitrogen-containing precursors from human inputs, such as urine, sweat, and skin cells. Many DBPs were new and have not been reported previously in either swimming pool or drinking waters. Bromoform levels were greater in brominated than in chlorinated pool waters, but we also identified many brominated DBPs in the chlorinated waters. The pool waters were mutagenic at levels similar to that of drinking water (approximately 1,200 revertants/L-equivalents in strain TA100-S9 mix).

CONCLUSIONS

This study identified many new DBPs not identified previously in swimming pool or drinking water and found that swimming pool waters are as mutagenic as typical drinking waters.

Evaluation of the impact of the exposure route on the human kinetic adjustment factor

The objective of this study was to assess the impact of the exposure route on the human kinetic adjustment factor (HKAF), for which a default value of 3.16 is used in non-cancer risk assessment. A multi-route PBPK model was modified from the literature and used for computing the internal dose metrics in adults, neonates, children, elderly and pregnant women following three route-specific scenarios to chloroform, bromoform, tri- or per-chloroethylene (TCE or PERC). These include 24-h inhalation exposure, body-weight adjusted oral exposure and 30 min dermal exposure to contaminated drinking water. Distributions for body weight (BW), height (BH) and hepatic cytochrome P450 2E1 (CYP2E1) content were obtained from the literature, whereas model parameters (flows, volumes) were calculated from BW and BH. Monte Carlo simulations were performed and the HKAF was calculated as the ratio of the 95th percentile value of internal dose metrics in subpopulation to the 50th percentile value in adults. On the basis of the area under the parent compound's arterial blood concentration vs time curve (AUC(pc)), highest HKAFs were obtained in neonates for every scenario considered, and were the highest for bromoform (range: 3.6-7.4). Exceedance of the default value based on AUC(PC) was also observed for an oral exposure to chloroform in neonates (4.9). In all other cases, HKAFs remained below the default value. Overall, this study has pointed out the dependency of the HKAF on the exposure route, dose metrics and subpopulation considered, as well as characteristics of the chemicals investigated.

A Bayesian population PBPK model for multiroute chloroform exposure

A Bayesian hierarchical model was developed to estimate the parameters in a physiologically based pharmacokinetic (PBPK) model for chloroform using prior information and biomarker data from different exposure pathways. In particular, the model provides a quantitative description of the changes in physiological parameters associated with hot-water bath and showering scenarios. Through Bayesian inference, uncertainties in the PBPK parameters were reduced from the prior distributions. Prediction of biomarker data with the calibrated PBPK model was improved by the calibration. The posterior results indicate that blood flow rates varied under two different exposure scenarios, with a two-fold increase of the skin's blood flow rate predicted in the hot-bath scenario. This result highlights the importance of considering scenario-specific parameters in PBPK modeling. To demonstrate the application of a probability approach in toxicological assessment, results from the posterior distributions from this calibrated model were used to predict target tissue dose based on the rate of chloroform metabolized in liver. This study demonstrates the use of the Bayesian approach to optimize PBPK model parameters for typical household exposure scenarios.

Chemical-specific screening criteria for interpretation of biomonitoring data for volatile organic compounds (VOCs)--application of steady-state PBPK model solutions

The National Health and Nutrition Examination Survey (NHANES) generates population-representative biomonitoring data for many chemicals including volatile organic compounds (VOCs) in blood. However, no health or risk-based screening values are available to evaluate these data from a health safety perspective or to use in prioritizing among chemicals for possible risk management actions. We gathered existing risk assessment-based chronic exposure reference values such as reference doses (RfDs), reference concentrations (RfCs), tolerable daily intakes (TDIs), cancer slope factors, etc. and key pharmacokinetic model parameters for 47 VOCs. Using steady-state solutions to a generic physiologically-based pharmacokinetic (PBPK) model structure, we estimated chemical-specific steady-state venous blood concentrations across chemicals associated with unit oral and inhalation exposure rates and with chronic exposure at the identified exposure reference values. The geometric means of the slopes relating modeled steady-state blood concentrations to steady-state exposure to a unit oral dose or unit inhalation concentration among 38 compounds with available pharmacokinetic parameters were 12.0 microg/L per mg/kg-d (geometric standard deviation [GSD] of 3.2) and 3.2 microg/L per mg/m(3) (GSD=1.7), respectively. Chemical-specific blood concentration screening values based on non-cancer reference values for both oral and inhalation exposure range from 0.0005 to 100 microg/L; blood concentrations associated with cancer risk-specific doses at the 1E-05 risk level ranged from 5E-06 to 6E-02 microg/L. The distribution of modeled steady-state blood concentrations associated with unit exposure levels across VOCs may provide a basis for estimating blood concentration screening values for VOCs that lack chemical-specific pharmacokinetic data. The screening blood concentrations presented here provide a tool for risk assessment-based evaluation of population biomonitoring data for VOCs and are most appropriately applied to central tendency estimates for such datasets.

Public health interpretation of trihalomethane blood levels in the United States: NHANES 1999-2004

Trihalomethanes (THMs) can form as byproducts during drinking water disinfection, which is crucial for limiting human exposure to disease-causing pathogens. The US Environmental Protection Agency (USEPA), recognizing both the importance of water disinfection for public health protection and potential risks associated with THM exposure, developed disinfection byproduct rules with the parallel goals of ensuring safe drinking water and limiting the levels of THMs in public water systems. The National Health and Nutrition Examination Survey (NHANES) THM blood data can be used as a means for assessing US population exposures to THMs; biomonitoring equivalents (BEs) can provide human health risk-based context to those data. In this paper, we examine the blood THM levels in the 1999-2004 NHANES data to (i) determine weighted population percentiles of blood THMs, (ii) explore whether gender and/or age are associated with blood THM levels, (iii) determine whether temporal trends can be discerned over the 6-year timeframe, and (iv) draw comparisons between population THM blood levels and BEs. A statistically significant decrease in blood chloroform levels was observed across the 1999-2004 time period. Age-related differences in blood chloroform levels were not consistent and no gender-related differences in blood chloroform levels were observed. The concentrations of all four THMs in the blood of US residents from the 2003 to 2004 NHANES dataset are below BEs consistent with the current US EPA reference doses. For bromodichloromethane and dibromochloromethane, the measured median blood concentrations in the United States are within the BEs for the 10(-6) and 10(-4) cancer risk range, whereas measured values for bromoform generally fall below the 10(-6) cancer risk range. These assessments indicate that general population blood concentrations of THMs are in a range considered to be a low to medium priority for risk assessment follow-up, according to the guidelines for interpretation of biomonitoring data using BEs.

Probabilistic human health risk assessment for quarterly exposure to high chloroform concentrations in drinking-water distribution network of the Province of Quebec, Canada

Because quarterly concentrations of total trihalomethanes (THM) exceeding the 80 μg/L guideline are often tolerated by the public health authorities of the Province of Quebec (Canada), this study examined whether quarterly episodes of high concentrations of THM may pose a risk to the health of its population. Using Monte Carlo simulations, a probabilistic risk assessment was performed for infants (0-<6 mo), toddlers (6 mo-<5 yr) and adults (≥20 yr). Multiroute exposure including ingestion of drinking water as well as inhalation and dermal exposure while showering or bathing was considered. The resulting absorbed doses were compared to short-term reference values for chloroform, used as surrogate for THM, by calculating risk quotients (RQ). On the basis of THM concentrations values in Quebec's drinking water distribution systems during the months of July to October and exceeding the guideline value (>80 μg/L), the 95th percentile value of RQ were 0.65, 0.46, and 0.24 for infants, toddlers, and adults, respectively. Back-calculation allowed determining that a chloroform concentration of 330 μg/L would result in RQ ≤ 1 for 99% of infants, the subgroup considered the most susceptible among the general population. Overall, this study showed that episodes of high THM concentration encountered in Quebec drinking-water distribution network need not be considered as an immediate health concern for the general population. However, these results should not be interpreted as an authorization to exceed the 80 μg/L standard but rather as a risk management tool for public health authorities.

Estimation of chlorination by-products presence in drinking water in epidemiological studies on adverse reproductive outcomes: a review

Chlorination of drinking water is essential to prevent waterborne disease. However, chlorine reacts with organic matter present in surface waters to form various by-products. In the last decade, several epidemiological studies have been conducted to determine the connection between exposure to these chlorination by-products (CBPs) and human health defects, such as adverse reproductive outcomes. However, the methodology used to assess exposure of pregnant women in these studies had serious limitations, particularly in relation to determining CBP presence in the subject's tap water. The purpose of this paper is to critically review of methods used to evaluate the CBP presence in a subject's tap water for exposure assessment purposes in epidemiological studies focused on adverse reproductive outcomes and CBPs in drinking water. Interest is directed more precisely at space-time features related to CBPs for an optimal estimation of their presence in a subject's tap water.

The effect of workload on biological monitoring of occupational exposure to toluene and n-Hexane: contribution of physiologically based toxicokinetic modeling

A physiologically based toxicokinetic model was used to examine the impact of work load on the relationship between the airborne concentrations and exposure indicator levels of two industrial solvents, toluene and n-Hexane. The authors simulated occupational exposure (8 hr/day, 5 days/week) at different concentrations, notably 20 ppm and 50 ppm, which are the current threshold limit values recommended by ACGIH for toluene and n-hexane, respectively. Different levels of physical activity, namely, rest, 25 W, and 50 W (for 12 hr followed by 12 hr at rest) were simulated to assess the impact of work load on the recommended biological exposure indices: toluene in blood prior to the last shift of the workweek, urinary o-cresol (a metabolite of toluene) at the end of the shift, and free (nonhydrolyzed) 2,5-hexanedione (a metabolite of n-hexane) at the end of the shift at the end of the workweek. In addition, urinary excretion of unchanged toluene was simulated. The predicted biological concentrations were compared with the results of both experimental studies among human volunteers and field studies among workers. The highest predicted increase with physical exercise was noted for toluene in blood (39 microg/L at 50 W vs. 14 microg/L at rest for 20 ppm, i.e., a 2.8-fold increase). The end-of-shift urinary concentrations of o-cresol and toluene were two times higher at 50 W than at rest (for 20 ppm, 0.65 vs. 0.33 mg/L for o-cresol and 43 vs. 21 microg/L for toluene). Urinary 2,5-hexanedione predicted for 50 ppm was 1.07 mg/L at 50 W and 0.92 mg/L at rest (+16%). The simulations that best describe the concentrations among workers exposed to toluene are those corresponding to 25 W or less. In conclusion, toxicokinetic modeling confirms the significant impact of work load on toluene exposure indicators, whereas only a very slight effect is noted on n-hexane kinetics. These results highlight the necessity of taking work load into account in risk assessment relative to toluene exposure.

Assessment of trichloroethylene (TCE) exposure in murine strains genetically-prone and non-prone to develop autoimmune disease

There is increasing laboratory and epidemiologic evidence relating exposure to trichloroethylene (TCE) with autoimmune disease including scleroderma and lupus. New Zealand Black/New Zealand White (NZBWF1) and B6C3F1 mice were exposed to TCE (0, 1, 400 or 14,000 ppb) via drinking water for 27 or 30 weeks, respectively. NZBWF1 mice spontaneously develop autoimmune disease while B6C3F1 mice, a standard strain used in immunotoxicology testing, are not genetically prone to develop autoimmune disease. During the TCE exposure period, serum levels of total IgG, and autoantibodies (anti-ssDNA, -dsDNA, and -glomerular antigen [GA]) were monitored. At the termination of the study, renal pathology, natural killer (NK) cell activity, total IgG levels, autoantibody production, T-cell activation, and lymphocytic proliferative responses were evaluated. TCE did not alter NK cell activity, or T- and B-cell proliferation in either strain. Numbers of activated T-cells (CD4+/CD44+) were increased in the B6C3F1 mice but not in the NZBWF1 mice. Renal pathology, as indicated by renal score, was significantly increased in the B6C3F1, but not in the NZBWF1 mice. Serum levels of autoantibodies to dsDNA and ssDNA were increased at more time points in B6C3F1, as compared to the NZBWF1 mice. Anti-GA autoantibodies were increased by TCE treatment in early stages of the study in NZBWF1 mice, but by 23 weeks of age, control levels were comparable to those of TCE-exposed animals. Serum levels anti-GA autoantibodies in B6C3F1 were not affected by TCE exposure. Overall, these data suggest that TCE did not contribute to the progression of autoimmune disease in autoimmune-prone mice during the period of 11-36 weeks of age, but rather lead to increased expression of markers associated with autoimmune disease in a non-genetically prone mouse strain.

A physiologically based pharmacokinetic model for the assessment of infant exposure to persistent organic pollutants in epidemiologic studies

BACKGROUND

It has been suggested that pre- and postnatal exposure to persistent organic pollutants (POPs) can promote several adverse effects in children, such as altered neurodevelopment. Epidemiologic studies to date have relied on the analysis of biological samples drawn pre- or post-natally for exposure assessment, an approach that might not capture some key events in the toxicokinetics of POPs.

OBJECTIVES

We aimed to build a generic physiologically based pharmacokinetic (PBPK) modeling framework for neutral POPs to assess infant toxicokinetic profiles and to validate the model using data on POP levels measured in mothers and infants from a Northern Québec Inuit population.

METHODS

The PBPK model developed herein was based upon a previously published model to which an infant submodel was added. Using the model and maternal blood levels at the time of delivery, exposure to 1,1-dichloro-2,2-bis(p-chlorophenyl)ethylene (p,p'-DDE), 1,1,1-trichloro-2,2-bis(p-chlorophenyl)ethane (p,p'-DDT), hexachlorobenzene (HCB), beta-hexachlorocyclohexane (beta-HCH), 2,2',3,4,4',5'-hexachlorobiphenyl (PCB-138), 2,2',4,4',5,5'-hexachlorobiphenyl (PCB-153), and 2,2',3,4,4',5,5'-heptachlorobiphenyl (PCB-180) in mothers was estimated to subsequently simulate infant blood, breast milk, and cord blood POP concentration. Simulations were then compared with corresponding measured levels through Spearman correlation analyses.

RESULTS

Predictions were highly correlated with measured concentrations for PCB-153, PCB-180, PCB-138, HCB, and p,p'-DDE (r = 0.83-0.96). Weaker correlations were observed for p,p'-DDT and beta-HCH for which levels were near the limits of detection.

CONCLUSION

This is the first study to validate a PBPK model of POPs in infants on an individual basis. This approach will reduce sampling efforts and enable the use of individualized POP toxicokinetic profiles in the epidemiologic studies of POP adverse effects on child development.