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

Cadmium accumulation in tropical island paddy soils: From environment and health risk assessment to model prediction

Cultivated soil quality is crucial because it directly affects food safety and human health, and rice is of primary concern because of its centrality to global food networks. However, a detailed understanding of cadmium (Cd) geochemical cycling in paddy soils is complicated by the multiple influencing factors present in many rice-growing areas that overlap with industrial centers. This study analyzed the pollution characteristics and health risks of Cd in paddy soils across Hainan Island and identified key influencing factors based on multi-source environmental data and prediction models. Approximately 27.07% of the soil samples exceeded the risk control standard screening value for Cd in China, posing an uncontaminated to moderate contamination risk. Cd concentration and exposure duration contributed the most to non-carcinogenic and carcinogenic risks to children, teens, and adults through ingestion. Among the nine prediction models tested, Extreme Gradient Boosting (XGBoost) exhibited the best performance for Cd prediction with soil properties having the highest importance, followed by climatic variables and topographic attributes. In summary, XGBoost reliably predicted the soil Cd concentrations on tropical islands. Further research should incorporate additional soil properties and environmental variables for more accurate predictions and to comprehensively identify their driving factors and corresponding contribution rates.

Estimation of toluene exposure in air from BMA (S-benzylmercapturic acid) urinary measures using a reverse dosimetry approach based on physiologically pharmacokinetic modeling

This study aimed to use a reverse dosimetry PBPK modeling approach to estimate toluene atmospheric exposure from urinary measurements of S-benzylmercapturic acid (BMA) in a small group of individuals and to evaluate the uncertainty associated to urinary spot-sampling compared to 24-h collected urine samples. Each exposure assessment technique was developed namely to estimate toluene air exposure from BMA measurements in 24-h urine samples (24-h-BMA) and from distributions of daily urinary BMA spot measurements (DUBSM). Model physiological parameters were described based upon age, weight, size and sex. Monte Carlo simulations with the PBPK model allowed converting DUBSM distribution (and 24-h-BMA) into toluene air levels. For the approach relying on DUBSM distribution, the ratio between the 95% probability of predicted toluene concentration and its 50% probability in each individual varied between 1.2 and 1.4, while that based on 24-h-BMA varied between 1.0 and 1.1. This suggests more variability in estimated exposure from spot measurements. Thus, estimating toluene exposure based on DUBSM distribution generated about 20% more uncertainty. Toluene levels estimated (0.0078-0.0138 ppm) are well below Health Canada's maximum chronic air guidelines. PBPK modeling and reverse dosimetry may be combined to interpret urinary metabolites data of VOCs and assess related uncertainties.

Development and evaluation of a high throughput inhalation model for organic chemicals

Currently it is difficult to prospectively estimate human toxicokinetics (particularly for novel chemicals) in a high-throughput manner. The R software package httk has been developed, in part, to address this deficiency, and the aim of this investigation was to develop a generalized inhalation model for httk. The structure of the inhalation model was developed from two previously published physiologically based models from Jongeneelen and Berge (Ann Occup Hyg 55:841-864, 2011) and Clewell et al. (Toxicol Sci 63:160-172, 2001), while calculated physicochemical data was obtained from EPA's CompTox Chemicals Dashboard. In total, 142 exposure scenarios across 41 volatile organic chemicals were modeled and compared to published data. The slope of the regression line of best fit between log-transformed simulated and observed blood and exhaled breath concentrations was 0.46 with an r2 = 0.45 and a root mean square error (RMSE) of direct comparison between the log-transformed simulated and observed values of 1.11. Approximately 5.1% (n = 108) of the data points analyzed were >2 orders of magnitude different than expected. The volatile organic chemicals examined in this investigation represent small, generally lipophilic molecules. Ultimately this paper details a generalized inhalation component that integrates with the httk physiologically based toxicokinetic model to provide high-throughput estimates of inhalation chemical exposures.

Evaluation of human biomonitoring data in a health risk based context: An updated analysis of population level data from the Canadian Health Measures Survey

In order to characterize exposure of the Canadian population to environmental chemicals, a human biomonitoring component has been included in the Canadian Health Measures Survey (CHMS). This nationally-representative survey, launched in 2007 by the Government of Canada, has measured over 250 chemicals in approximately 30,000 Canadians during the last decade. The capacity to interpret these data at the population level in a health risk context is gradually improving with the development of biomonitoring screening values, such as biomonitoring equivalents (BE) and human biomonitoring (HBM) values. This study evaluates recent population level biomonitoring data from the CHMS in a health risk context using biomonitoring screening values. Nationally representative biomonitoring data for fluoride, selenium, molybdenum, arsenic, silver, thallium, cyfluthrin, 2,4-dichlorophenoxyacetic acid (2,4-D), 3-phenoxybenzoic acid (3-PBA), chlorpyrifos, deltamethrin, bisphenol A, triclosan, acrylamide, cadmium, perfluorooctane sulfonate (PFOS), perfluorooctanoic acid (PFOA), bromoform, chloroform, benzene, toluene, xylene, ethylbenzene, styrene and tetrachloroethylene were screened as part as this study. For non-cancer endpoints, hazard quotients (HQs) were calculated as the ratio of population level concentrations of a specific chemical at the geometric mean and 95th percentile to the corresponding biomonitoring screening value. Cancer risks were calculated at the 5th, 25th, 50th, 75th and 95th percentiles of the population concentration using BEs based on a risk specific dose. Most of the chemicals analyzed had HQs below 1 suggesting that levels of exposure to these chemicals are not a concern at the population level. However, HQs exceeded 1 in smokers for cadmium, acrylamide and benzene, as well as in the general population for inorganic arsenic, PFOS and PFOA, 3-PBA and fluoride. Furthermore, cancer risks for inorganic arsenic, acrylamide, and benzene at most population percentiles of exposure were elevated (>10^-5). Specifically, for inorganic arsenic in the general population, the HQ was 3.13 at the 95th percentile concentration and the cancer risk was 3.4 × 10^-4 at the 50th percentile of population concentrations. These results suggest that the levels of exposure in the Canadian population to some of the environmental chemicals assessed might be of concern. The results of this screening exercise support the findings of previous risk assessments and ongoing efforts to reduce risks from exposure to chemicals evaluated as part of this study. Although paucity of biomonitoring screening values for several environmental contaminants may be a limitation to this approach, our assessment contributes to the prioritization of a number of chemicals measured as part of CHMS for follow-up activities such as more detailed characterization of exposure sources.

Evaluation of potential health effects associated with occupational and environmental exposure to styrene - an update

The potential chronic health risks of occupational and environmental exposure to styrene were evaluated to update health hazard and exposure information developed since the Harvard Center for Risk Analysis risk assessment for styrene was performed in 2002. The updated hazard assessment of styrene's health effects indicates human cancers and ototoxicity remain potential concerns. However, mechanistic research on mouse lung tumors demonstrates these tumors are mouse-specific and of low relevance to human cancer risk. The updated toxicity database supports toxicity reference levels of 20 ppm (equates to 400 mg urinary metabolites mandelic acid + phenylglyoxylic acid/g creatinine) for worker inhalation exposure and 3.7 ppm and 2.5 mg/kg bw/day, respectively, for general population inhalation and oral exposure. No cancer risk value estimates are proposed given the established lack of relevance of mouse lung tumors and inconsistent epidemiology evidence. The updated exposure assessment supports inhalation and ingestion routes as important. The updated risk assessment found estimated risks within acceptable ranges for all age groups of the general population and workers with occupational exposures in non-fiber-reinforced polymer composites industries and fiber-reinforced polymer composites (FRP) workers using closed-mold operations or open-mold operations with respiratory protection. Only FRP workers using open-mold operations not using respiratory protection have risk exceedances for styrene and should be considered for risk management measures. In addition, given the reported interaction of styrene exposure with noise, noise reduction to sustain levels below 85 dB(A) needs be in place.

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.

Human health risk assessment and PAHs in a stretch of river Ganges near Kanpur

The presence of 13 polycyclic aromatic hydrocarbons (PAHs) was measured in a small stretch of river Ganges in Kanpur, using high-performance liquid chromatography equipped with photodiode array detector (HPLC/PDA). Non-carcinogenic risk on human health was calculated in the form of hazardous index. Carcinogenic risk was calculated as chronic daily intake and incremental lifetime cancer risk (ILCR) with the help of monitored PAHs in river water. Due to non-availability of reference dose (RfD) values for all the 13 PAHs, the combined hazardous index for three PAHs was calculated. These values in winter, pre-monsoon, and post-monsoon period varied from 6.37 × 10^-3 (Bithoor) to 1.12 × 10^-2 (Jajmau), 2.89 × 10^-3 (Bithoor) to 8.52 × 10^-3 (Annandeshwar Temple), and 1.80 × 10^-3 (Massacre Ghat) to 6.67 × 10^-3 (Jajmau), respectively. In this study, the carcinogenic risk due to PAHs was calculated in the form of ILCR. ILCR due to benzo(a)pyrene (BaP) varied from 1.39 × 10^-6 (Massacre Ghat) to 1.45 × 10^-5 (Jajmau). ILCR was also calculated with the help of BaP_eq for five age groups of people (adults, teenagers, children, toddlers, and infants). The outcome of the study indicates that there is a need to control pollution of the river water to maintain its quality. Continuous discharge of PAHs into the river poses both human health risk and ecological risk.

Breath volatile organic compounds for the gut-fatty liver axis: Promise, peril, and path forward

The worldwide interest in the gut microbiome and its impact on the upstream liver highlight a critical upside to breath research: it can uniquely measure otherwise unmeasurable biology. Bacteria make gases [volatile organic compounds (VOCs)] that are directly relevant to pathophysiology of the fatty liver and associated conditions, including obesity. Measurement of these VOCs and their metabolites in the exhaled breath, therefore, present an opportunity to safely and easily evaluate, on both a personal and a population level, some of our most pressing public health threats. This is an opportunity that must be pursued. To date, however, breath analysis remains a slowly evolving field which only occasionally impacts clinical research or patient care. One major obstacle to progress is that breath analysis is inherently and emphatically mutli-disciplinary: it connects engineering, chemistry, breath mechanics, biology and medicine. Unbalanced or incomplete teams may produce inconsistent and often unsatisfactory results. A second impediment is the lack of a well-known stepwise structure for the development of non-invasive diagnostics. As a result, the breath research landscape is replete with orphaned single-center pilot studies. Often, important hypotheses and key observations have not been pursued to maturation. This paper reviews the rationale and requirements for breath VOC research applied to the gut-fatty liver axis and offers some suggestions for future development.

Evaluation of biomonitoring data from the CDC National Exposure Report in a risk assessment context: perspectives across chemicals

BACKGROUND

Biomonitoring data reported in the National Report on Human Exposure to Environmental Chemicals [NER; Centers for Disease Control and Prevention (2012)] provide information on the presence and concentrations of > 400 chemicals in human blood and urine. Biomonitoring Equivalents (BEs) and other risk assessment-based values now allow interpretation of these biomonitoring data in a public health risk context.

OBJECTIVES

We compared the measured biomarker concentrations in the NER with BEs and similar risk assessment values to provide an across-chemical risk assessment perspective on the measured levels for approximately 130 analytes in the NER.

METHODS

We identified available risk assessment-based biomarker screening values, including BEs and Human Biomonitoring-I (HBM-I) values from the German Human Biomonitoring Commission. Geometric mean and 95th percentile population biomarker concentrations from the NER were compared to the available screening values to generate chemical-specific hazard quotients (HQs) or cancer risk estimates.

CONCLUSIONS

Most analytes in the NER show HQ values of < 1; however, some (including acrylamide, dioxin-like chemicals, benzene, xylene, several metals, di-2(ethylhexyl)phthalate, and some legacy organochlorine pesticides) approach or exceed HQ values of 1 or cancer risks of > 1 × 10-4 at the geometric mean or 95th percentile, suggesting exposure levels may exceed published human health benchmarks. This analysis provides for the first time a means for examining population biomonitoring data for multiple environmental chemicals in the context of the risk assessments for those chemicals. The results of these comparisons can be used to focus more detailed chemical-specific examination of the data and inform priorities for chemical risk management and research.

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.

Application of physiologically based pharmacokinetic models in chemical risk assessment

Post-exposure risk assessment of chemical and environmental stressors is a public health challenge. Linking exposure to health outcomes is a 4-step process: exposure assessment, hazard identification, dose response assessment, and risk characterization. This process is increasingly adopting "in silico" tools such as physiologically based pharmacokinetic (PBPK) models to fine-tune exposure assessments and determine internal doses in target organs/tissues. Many excellent PBPK models have been developed. But most, because of their scientific sophistication, have found limited field application-health assessors rarely use them. Over the years, government agencies, stakeholders/partners, and the scientific community have attempted to use these models or their underlying principles in combination with other practical procedures. During the past two decades, through cooperative agreements and contracts at several research and higher education institutions, ATSDR funded translational research has encouraged the use of various types of models. Such collaborative efforts have led to the development and use of transparent and user-friendly models. The "human PBPK model toolkit" is one such project. While not necessarily state of the art, this toolkit is sufficiently accurate for screening purposes. Highlighted in this paper are some selected examples of environmental and occupational exposure assessments of chemicals and their mixtures.

Development of screening tools for the interpretation of chemical biomonitoring data

Evaluation of a larger number of chemicals in commerce from the perspective of potential human health risk has become a focus of attention in North America and Europe. Screening-level chemical risk assessment evaluations consider both exposure and hazard. Exposures are increasingly being evaluated through biomonitoring studies in humans. Interpreting human biomonitoring results requires comparison to toxicity guidance values. However, conventional chemical-specific risk assessments result in identification of toxicity-based exposure guidance values such as tolerable daily intakes (TDIs) as applied doses that cannot directly be used to evaluate exposure information provided by biomonitoring data in a health risk context. This paper describes a variety of approaches for development of screening-level exposure guidance values with translation from an external dose to a biomarker concentration framework for interpreting biomonitoring data in a risk context. Applications of tools and concepts including biomonitoring equivalents (BEs), the threshold of toxicologic concern (TTC), and generic toxicokinetic and physiologically based toxicokinetic models are described. These approaches employ varying levels of existing chemical-specific data, chemical class-specific assessments, and generic modeling tools in response to varying levels of available data in order to allow assessment and prioritization of chemical exposures for refined assessment in a risk management context.

Evaluation of NHANES biomonitoring data for volatile organic chemicals in blood: application of chemical-specific screening criteria

Available biomonitoring data for volatile organic chemicals (VOCs) in blood from the National Health and Nutrition Examination Survey (NHANES) (2003-2004) (CDC, 2009) were compared with recently derived screening biomonitoring equivalent (BE) values (Aylward et al., 2010). A BE is defined as the estimated concentration or range of concentrations of a chemical or its metabolites in a biological medium (blood, urine, or other medium) that is consistent with an existing health-based exposure guidance value. Blood concentrations of VOCs from the NHANES data set were compared with predicted screening BE values based upon a hazard quotient (HQ) for individual chemicals, and a hazard index (HI) approach for combined exposures. HI values for detected chemicals were generally at or below a value of 1, suggesting that the potential for deleterious effects is low. However, smoking was an important determinant of HI and HQ values. Detected levels of benzene in non-smokers were within the range of BE values corresponding to a 1 × 10(-6)-1 × 10(-4) range for upper-bound cancer risk; in smokers, levels of benzene were at the upper end of or exceeded this range. For VOCs that were not detected in the NHANES sampling, analytical detection limits were generally sufficiently sensitive to detect concentrations consistent with existing non-cancer and cancer risk-based exposure guidance values. Interpretations of measured blood concentrations of VOCs must be made with caution due to the substantial within-individual, within-day fluctuations in levels expected due to the rapid elimination of VOCs.

Health risk assessment of polycyclic aromatic hydrocarbons in the source water and drinking water of China: Quantitative analysis based on published monitoring data

A carcinogenic risk assessment of polycyclic aromatic hydrocarbons (PAHs) in source water and drinking water of China was conducted using probabilistic techniques from a national perspective. The published monitoring data of PAHs were gathered and converted into BaP equivalent (BaP(eq)) concentrations. Based on the transformed data, comprehensive risk assessment was performed by considering different age groups and exposure pathways. Monte Carlo simulation and sensitivity analysis were applied to quantify uncertainties of risk estimation. The risk analysis indicated that, the risk values for children and teens were lower than the accepted value (1.00E-05), indicating no significant carcinogenic risk. The probability of risk values above 1.00E-05 was 5.8% and 6.7% for adults and lifetime groups, respectively. Overall, carcinogenic risks of PAHs in source water and drinking water of China were mostly accepted. However, specific regions, such as Yellow river of Lanzhou reach and Qiantang river should be paid more attention. Notwithstanding the uncertainties inherent in the risk assessment, this study is the first attempt to provide information on carcinogenic risk of PAHs in source water and drinking water of China, and might be useful for potential strategies of carcinogenic risk management and reduction.