Methods to account for uncertainties in exposure assessment in studies of environmental exposures

A Historical Survey of Key Epidemiological Studies of Ionizing Radiation Exposure

In this article we review the history of key epidemiological studies of populations exposed to ionizing radiation. We highlight historical and recent findings regarding radiation-associated risks for incidence and mortality of cancer and non-cancer outcomes with emphasis on study design and methods of exposure assessment and dose estimation along with brief consideration of sources of bias for a few of the more important studies. We examine the findings from the epidemiological studies of the Japanese atomic bomb survivors, persons exposed to radiation for diagnostic or therapeutic purposes, those exposed to environmental sources including Chornobyl and other reactor accidents, and occupationally exposed cohorts. We also summarize results of pooled studies. These summaries are necessarily brief, but we provide references to more detailed information. We discuss possible future directions of study, to include assessment of susceptible populations, and possible new populations, data sources, study designs and methods of analysis.

Recommendations on statistical approaches to account for dose uncertainties in radiation epidemiologic risk models

PURPOSE

Epidemiological studies of stochastic radiation health effects such as cancer, meant to estimate risks of the adverse effects as a function of radiation dose, depend largely on estimates of the radiation doses received by the exposed group under study. Those estimates are based on dosimetry that always has uncertainty, which often can be quite substantial. Studies that do not incorporate statistical methods to correct for dosimetric uncertainty may produce biased estimates of risk and incorrect confidence bounds on those estimates. This paper reviews commonly used statistical methods to correct radiation risk regressions for dosimetric uncertainty, with emphasis on some newer methods. We begin by describing the types of dose uncertainty that may occur, including those in which an uncertain value is shared by part or all of a cohort, and then demonstrate how these sources of uncertainty arise in radiation dosimetry. We briefly describe the effects of different types of dosimetric uncertainty on risk estimates, followed by a description of each method of adjusting for the uncertainty.

CONCLUSIONS

Each of the method has strengths and weaknesses, and some methods have limited applicability. We describe the types of uncertainty to which each method can be applied and its pros and cons. Finally, we provide summary recommendations and touch briefly on suggestions for further research.

Micro-environmental personal radio-frequency electromagnetic field exposures in Melbourne: A longitudinal trend analysis

BACKGROUND

A knowledge gap exists regarding longitudinal assessment of personal radio-frequency electromagnetic field (RF-EMF) exposures globally. It is unclear how the change in telecommunication technology over the years translates to change in RF-EMF exposure. This study aims to evaluate longitudinal trends of micro-environmental personal RF-EMF exposures in Australia.

METHODS

The study utilised baseline (2015-16) and follow-up (2022) data on personal RF-EMF exposure (88 MHz-6 GHz) measured across 18 micro-environments in Melbourne. Simultaneous quantile regression analysis was conducted to compare exposure data distribution percentiles, particularly median (P50), upper extreme value (P99) and overall exposure trends. RF-EMF exposures were compared across six exposure source types: mobile downlink, mobile uplink, broadcast, 5G-New Radio, Others and Total (of the aforementioned sources). Frequency-specific exposures measured at baseline and follow-up were compared. Total exposure across different groups of micro-environment types were also compared.

RESULTS

For all micro-environmental data, total (median and P99) exposure levels did not significantly change at follow-up. Overall exposure trend of total exposure increased at follow-up. Mobile downlink contributed the highest exposure among all sources showing an increase in median exposure and overall exposure trend. Of seven micro-environment types, five of them showed total exposure levels (median and P99) and overall exposure trend increased at follow-up.

Adopting Mechanistic Molecular Biology Approaches in Exposome Research for Causal Understanding

Through investigating the combined impact of the environmental exposures experienced by an individual throughout their lifetime, exposome research provides opportunities to understand and mitigate negative health outcomes. While current exposome research is driven by epidemiological studies that identify associations between exposures and effects, new frameworks integrating more substantial population-level metadata, including electronic health and administrative records, will shed further light on characterizing environmental exposure risks. Molecular biology offers methods and concepts to study the biological and health impacts of exposomes in experimental and computational systems. Of particular importance is the growing use of omics readouts in epidemiological and clinical studies. This paper calls for the adoption of mechanistic molecular biology approaches in exposome research as an essential step in understanding the genotype and exposure interactions underlying human phenotypes. A series of recommendations are presented to make the necessary and appropriate steps to move from exposure association to causation, with a huge potential to inform precision medicine and population health. This includes establishing hypothesis-driven laboratory testing within the exposome field, supported by appropriate methods to read across from model systems research to human.

Integrated environmental health assessment: Proposed approaches to exposure during chemical incidents

An integrated environmental health exposure assessment (IEHA) refers to the integration of human biomonitoring data (HBM) and environmental measurements and aims to optimize the exposure assessment process. Due to lack of data, this approach remains an issue during chemical incidents. This study aims to explore integrated exposure approaches for assessing human health risks during chemical incidents. Based on the Preferred Reporting Items of Systematic reviews and Meta-Analyses statement, a literature analysis was performed. A level of confidence ranging from 1 to 4 was established to define the quality and strength of data used to undertake an IEHA approach. Twenty-seven articles (n = 18) and texts (n = 9) from Europe (41%) and the United States (37%) were analyzed. Among the 18 scientific articles, 61% (n = 11) presented a quantitative approach and 17% (n = 3) presented a qualitative approach. Quantitative approaches must be based on accurate data, coupled with predictive models. Of all the scientific papers, 40% (n = 7) responded to a confidence level greater than or equal to 2. Uncertainties detected through the integrated exposure approaches were related to input data, analytical methods, and HBM reference value interpretations. During chemical incidents, direct measurements were the most relevant data. Few scientific studies have developed an integrated approach during emergency situations. However, when this was used, they presented a high level of confidence by defining levels of exposure that support decision-making processes. Despite the multiple approaches, there was a lack of guidelines allowing an integrated risk assessment to be performed during an emergency chemical exposure. Integr Environ Assess Manag 2024;20:481-497. © 2023 SETAC.

Optimal-design domain-adaptation for exposure prediction in two-stage epidemiological studies

BACKGROUND

In the first stage of a two-stage study, the researcher uses a statistical model to impute the unobserved exposures. In the second stage, imputed exposures serve as covariates in epidemiological models. Imputation error in the first stage operate as measurement errors in the second stage, and thus bias exposure effect estimates.

OBJECTIVE

This study aims to improve the estimation of exposure effects by sharing information between the first and second stages.

METHODS

At the heart of our estimator is the observation that not all second-stage observations are equally important to impute. We thus borrow ideas from the optimal-experimental-design theory, to identify individuals of higher importance. We then improve the imputation of these individuals using ideas from the machine-learning literature of domain adaptation.

RESULTS

Our simulations confirm that the exposure effect estimates are more accurate than the current best practice. An empirical demonstration yields smaller estimates of PM effect on hyperglycemia risk, with tighter confidence bands.

SIGNIFICANCE

Sharing information between environmental scientist and epidemiologist improves health effect estimates. Our estimator is a principled approach for harnessing this information exchange, and may be applied to any two stage study.

A generalisation of the method of regression calibration

There is direct evidence of risks at moderate and high levels of radiation dose for highly radiogenic cancers such as leukaemia and thyroid cancer. For many cancer sites, however, it is necessary to assess risks via extrapolation from groups exposed at moderate and high levels of dose, about which there are substantial uncertainties. Crucial to the resolution of this area of uncertainty is the modelling of the dose-response relationship and the importance of both systematic and random dosimetric errors for analyses in the various exposed groups. It is well recognised that measurement error can alter substantially the shape of this relationship and hence the derived population risk estimates. Particular attention has been devoted to the issue of shared errors, common in many datasets, and particularly important in occupational settings. We propose a modification of the regression calibration method which is particularly suited to studies in which there is a substantial amount of shared error, and in which there may also be curvature in the true dose response. This method can be used in settings where there is a mixture of Berkson and classical error. In fits to synthetic datasets in which there is substantial upward curvature in the true dose response, and varying (and sometimes substantial) amounts of classical and Berkson error, we show that the coverage probabilities of all methods for the linear coefficient [Formula: see text] are near the desired level, irrespective of the magnitudes of assumed Berkson and classical error, whether shared or unshared. However, the coverage probabilities for the quadratic coefficient [Formula: see text] are generally too low for the unadjusted and regression calibration methods, particularly for larger magnitudes of the Berkson error, whether this is shared or unshared. In contrast Monte Carlo maximum likelihood yields coverage probabilities for [Formula: see text] that are uniformly too high. The extended regression calibration method yields coverage probabilities that are too low when shared and unshared Berkson errors are both large, although otherwise it performs well, and coverage is generally better than these other three methods. A notable feature is that for all methods apart from extended regression calibration the estimates of the quadratic coefficient [Formula: see text] are substantially upwardly biased.

A generalisation of the method of regression calibration

There is direct evidence of risks at moderate and high levels of radiation dose for highly radiogenic cancers such as leukaemia and thyroid cancer. For many cancer sites, however, it is necessary to assess risks via extrapolation from groups exposed at moderate and high levels of dose, about which there are substantial uncertainties. Crucial to the resolution of this area of uncertainty is the modelling of the dose-response relationship and the importance of both systematic and random dosimetric errors for analyses in the various exposed groups. It is well recognised that measurement error can alter substantially the shape of this relationship and hence the derived population risk estimates. Particular attention has been devoted to the issue of shared errors, common in many datasets, and particularly important in occupational settings. We propose a modification of the regression calibration method which is particularly suited to studies in which there is a substantial amount of shared error, and in which there may also be curvature in the true dose response. This method can be used in settings where there is a mixture of Berkson and classical error. In fits to synthetic datasets in which there is substantial upward curvature in the true dose response, and varying (and sometimes substantial) amounts of classical and Berkson error, we show that the coverage probabilities of all methods for the linear coefficient \(\alpha\) are near the desired level, irrespective of the magnitudes of assumed Berkson and classical error, whether shared or unshared. However, the coverage probabilities for the quadratic coefficient \(\beta\) are generally too low for the unadjusted and regression calibration methods, particularly for larger magnitudes of the Berkson error, whether this is shared or unshared. In contrast Monte Carlo maximum likelihood yields coverage probabilities for \(\beta\) that are uniformly too high. The extended regression calibration method yields coverage probabilities that are too low when shared and unshared Berkson errors are both large, although otherwise it performs well, and coverage is generally better than these other three methods. A notable feature is that for all methods apart from extended regression calibration the estimates of the quadratic coefficient \(\beta\) are substantially upwardly biased.

Dose estimates and their uncertainties for use in epidemiological studies of radiation-exposed populations in the Russian Southern Urals

Many residents of the Russian Southern Urals were exposed to radioactive environmental pollution created by the operations of the Mayak Production Association in the mid- 20th century. There were two major releases: the discharge of about 1x1017 Bq of liquid waste into the Techa River between 1949 and 1959; and the atmospheric release of 7.4 * 1016 Bq as a result an explosion in the radioactive waste-storage facility in 1957. The releases into the Techa River resulted in the exposure of more than 30,000 people who lived in riverside villages between 1950 and 1961. The 1957 accident contaminated a larger area with the highest exposure levels in an area that is called the East Urals Radioactive Trace (EURT). Current epidemiologic studies of the exposed populations are based on dose estimates obtained using a Monte-Carlo dosimetry system (TRDS-2016MC) that provides multiple realizations of the annual doses for each cohort member. These dose realizations provide a central estimate of the individual dose and information on the uncertainty of these dose estimates. In addition, the correlation of individual annual doses over realizations provides important information on shared uncertainties that can be used to assess the impact of shared dose uncertainties on risk estimate uncertainty.This paper considers dose uncertainties in the TRDS-2016MC. Individual doses from external and internal radiation sources were reconstructed for 48,036 people based on environmental contamination patterns, residential histories, individual 90Sr body-burden measurements and dietary intakes. Dietary intake of 90Sr resulted in doses accumulated in active bone marrow (or simply, marrow) that were an order of magnitude greater than those in soft tissues. About 84% of the marrow dose and 50% of the stomach dose was associated with internal exposures. The lognormal distribution is well-fitted to the individual dose realizations, which, therefore, could be expressed and easily operated in terms of geometric mean (GM) and geometric standard deviation (GSD). Cohort average GM for marrow and stomach cumulative doses are 0.21 and 0.03 Gy, respectively. Cohort average dose uncertainties in terms of GSD are as follows: for marrow it is 2.93 (90%CI: 2.02-4.34); for stomach and the other non-calcified tissues it is 2.32 (90% CI: 1.78-2.9).

Exposure to air pollution and gains in body weight and waist circumference among middle-aged and older adults

INTRODUCTION

Emerging research suggested a nexus between air pollution exposure and risks of overweight and obesity, while existing longitudinal evidence was extensively sparse, particularly in densely populated regions. This study aimed to quantify concentration-response associations of changes in weight and waist circumference (WC) related to air pollution in Chinese adults.

METHODS

We conceived a nationally representative longitudinal study from 2011 to 2015, by collecting 34,854 observations from 13,757 middle-aged and older adults in 28 provincial regions of China. Participants' height, weight and WC were measured by interviewers using standardized devices. Concentrations of major air pollutants including fine particulate matter (PM_2.5), nitrogen dioxide (NO₂) and ozone (O₃) predicted by well-validated spatiotemporal models were assigned to participants according to their residential cities. Possible exposure biases were checked through 1000 random simulated exposure at individual level, using a Monte Carlo simulation approach. Linear mixed-effects models were applied to estimate the relationships of air pollution with weight and WC changes, and restricted cubic spline functions were adopted to smooth concentration-response (C-R) curves.

RESULTS

Each 10-μg/m³ rise in PM_2.5, NO₂ and O₃ was associated with an increase of 0.825 (95% confidence interval: 0.740, 0.910), 0.921 (0.811, 1.032) and 1.379 (1.141, 1.616) kg in weight, respectively, corresponding to WC gains of 0.688 (0.592, 0.784), 1.189 (1.040, 1.337) and 0.740 (0.478, 1.002) cm. Non-significant violation for linear C-R relationships was observed with exception of NO₂-weight and PM_2.5/NO₂-WC associations. Sex-stratified analyses revealed elevated vulnerability in women to gain of weight in exposure to PM_2.5 and NO₂. Sensitive analyses largely supported our primary findings via assessing exposure estimates from 1000 random simulations, and performing reanalysis based on non-imputed covariates and non-obese participants, as well as alternative indicators (i.e., body mass index and waist-to-height ratio).

CONCLUSIONS

We found positively robust associations of later-life exposure to air pollutants with gains in weight and WC based on a national sample of Chinese adult men and women. Our findings suggested that mitigation of air pollution may be an efficient intervention to relieve obesity burden.

Human health risk assessment for contaminated sites: A retrospective review

Soil contamination is a serious global hazard as contaminants can migrate to the human body through the soil, water, air, and food, threatening human health. Human Health Risk Assessment (HHRA) is a commonly used method for estimating the magnitude and probability of adverse health effects in humans that may be exposed to contaminants in contaminated environmental media in the present or future. Such estimations have improved for decades with various risk assessment frameworks and well-established models. However, the existing literature does not provide a comprehensive overview of the methods and models of HHRA that are needed to grasp the current status of HHRA and future research directions. Thus, this paper aims to systematically review the HHRA approaches and models, particularly those related to contaminated sites from peer-reviewed literature and guidelines. The approaches and models focus on methods used in hazard identification, toxicity databases in dose-response assessment, approaches and fate and transport models in exposure assessment, risk characterization, and uncertainty characterization. The features and applicability of the most commonly used HHRA tools are also described. The future research trend for HHRA for contaminated sites is also forecasted. The transition from animal experiments to new methods in risk identification, the integration and update and sharing of existing toxicity databases, the integration of human biomonitoring into the risk assessment process, and the integration of migration and transformation models and risk assessment are the way forward for risk assessment in the future. This review provides readers with an overall understanding of HHRA and a grasp of its developmental direction.

Uncertainty and sensitivity analyses of human health risk from bioaccessible arsenic exposure via rice ingestion in Bangkok, Thailand

BACKGROUND

Rice can be a source of arsenic (As) exposure, causing health impacts after ingestion.

OBJECTIVE

This study analyzed health risks due to As exposure through rice consumption, focusing on both bioaccessible (bAs) and total (tAs) As levels.

METHODS

Monte Carlo simulations were applied to determine health risk uncertainties and to analyze factors influencing health risks.

RESULTS

Cooked white and brown rice contained lower tAs and bAs than FAO/WHO standards of 0.20 and 0.35 mg/kg, respectively. As became less bioaccessible after cooking (14.0% in white rice and 18.5% in brown rice). Non-carcinogenic effects (MOS < 1) were found in 5% of children. Carcinogenic effects (MOE<100), especially lung cancer, were found in 75% of adults, with a probable incidence of 7 in 1,000,000. The lowest and highest annual cancer cases were 18 in 10,000,000 adolescents and 15 in 1,000,000 adults, respectively. The risks were mainly affected by body weight and bAs concentration.

SIGNIFICANCE

The results identified a certain risk level of non-carcinogenic effects in children and adolescents as well as carcinogenic effects in adults. The per capita consumption of rice in Thai adults should be reduced to prevent incidences of lung cancer.

Mathematical characterization of dose uncertainty effects on functions summarizing findings of community noise attitudinal surveys

Previous Monte Carlo simulations have quantified the extent to which dose (sound level) uncertainty in community noise dose-response surveys can bias the shape of inferred dose-response functions. The present work extends the prior findings to create a mathematical model of the biasing effect. The exact effect on any particular data set depends on additional attributes (situational variables) beyond dose uncertainty itself. Several variables and their interaction effects are accounted for in the model. The model produced identical results to the prior Monte Carlo simulations and thereby demonstrated the same slope reduction effect. This model was further exercised to demonstrate the nature and extent of situational variable interaction effects related to the range of doses employed and their distribution across the range. One manifestation was a false asymptotic behavior in the observed dose-response relationship. The mathematical model provides a means to not only predict dose uncertainty effects but also to serve as a foundation for correcting for such effects in regression analyses of transportation noise dose-response relationships.

An integrated exposure and pharmacokinetic modeling framework for assessing population-scale risks of phthalates and their substitutes

To effectively incorporate in vitro-in silico-based methods into the regulation of consumer product safety, a quantitative connection between product phthalate concentrations and in vitro bioactivity data must be established for the general population. We developed, evaluated, and demonstrated a modeling framework that integrates exposure and pharmacokinetic models to convert product phthalate concentrations into population-scale risks for phthalates and their substitutes. A probabilistic exposure model was developed to generate the distribution of multi-route exposures based on product phthalate concentrations, chemical properties, and human activities. Pharmacokinetic models were developed to simulate population toxicokinetics using Bayesian analysis via the Markov chain Monte Carlo method. Both exposure and pharmacokinetic models demonstrated good predictive capability when compared with worldwide studies. The distributions of exposures and pharmacokinetics were integrated to predict the population distributions of internal dosimetry. The predicted distributions showed reasonable agreement with the U.S. biomonitoring surveys of urinary metabolites. The "source-to-outcome" local sensitivity analysis revealed that food contact materials had the greatest impact on body burden for di(2-ethylhexyl) adipate (DEHA), di-2-ethylhexyl phthalate (DEHP), di(isononyl) cyclohexane-1,2-dicarboxylate (DINCH), and di(2-propylheptyl) phthalate (DPHP), whereas the body burden of diethyl phthalate (DEP) was most sensitive to the concentration in personal care products. The upper bounds of predicted plasma concentrations showed no overlap with ToxCast in vitro bioactivity values. Compared with the in vitro-to-in vivo extrapolation (IVIVE) approach, the integrated modeling framework has significant advantages in mapping product phthalate concentrations to multi-route risks, and thus is of great significance for regulatory use with a relatively low input requirement. Further integration with new approach methodologies will facilitate these in vitro-in silico-based risk assessments for a broad range of products containing an equally broad range of chemicals.

Systematic Review Methodologies and Endocrine Disrupting Chemicals: Improving Evaluations of the Plastic Monomer Bisphenol A

BACKGROUND

Endocrine disrupting chemicals (EDCs) are found in plastics, personal care products, household items, and other consumer goods. Risk assessments are intended to characterize a chemical's hazards, identify the doses at which adverse outcomes are observed, quantify exposure levels, and then compare these doses to determine the likelihood of risk in a given population. There are many problems with risk assessments for EDCs, allowing people to be exposed to levels that are later associated with serious health outcomes in epidemiology studies.

OBJECTIVE

In this review, we examine issues that affect the evaluation of EDCs in risk assessments (e.g., use of insensitive rodent strains and absence of disease-oriented outcomes in hazard assessments; inadequate exposure assessments). We then review one well-studied chemical, Bisphenol A (BPA; CAS #80-05-7) an EDC found in plastics, food packaging, and other consumer products. More than one hundred epidemiology studies suggest associations between BPA exposures and adverse health outcomes in environmentally exposed human populations.

FINDINGS

We present support for the use of systematic review methodologies in the evaluation of BPA and other EDCs. Systematic reviews would allow studies to be evaluated for their reliability and risk of bias. They would also allow all data to be used in risk assessments, which is a requirement for some regulatory agencies.

CONCLUSION

Systematic review methodologies can be used to improve evaluations of BPA and other EDCs. Their use could help to restore faith in risk assessments and ensure that all data are utilized in decision-making. Regulatory agencies are urged to conduct transparent, well-documented and proper systematic reviews for BPA and other EDCs.

Fluoroscopy X-Ray Organ-Specific Dosimetry System (FLUXOR) for Estimation of Organ Doses and Their Uncertainties in the Canadian Fluoroscopy Cohort Study

As part of ongoing efforts to assess lifespan disease mortality and incidence in 63,715 patients from the Canadian Fluoroscopy Cohort Study (CFCS) who were treated for tuberculosis between 1930 and 1969, we developed a new FLUoroscopy X-ray ORgan-specific dosimetry system (FLUXOR) to estimate radiation doses to various organs and tissues. Approximately 45% of patients received medical procedures accompanied by fluoroscopy, including artificial pneumothorax (air in pleural cavity to collapse of lungs), pneumoperitoneum (air in peritoneal cavity), aspiration of fluid from pleural cavity and gastrointestinal series. In addition, patients received chest radiographs for purposes of diagnosis and monitoring of disease status. FLUXOR utilizes age-, sex- and body size-dependent dose coefficients for fluoroscopy and radiography exams, estimated using radiation transport simulations in up-to-date computational hybrid anthropomorphic phantoms. The phantoms include an updated heart model, and were adjusted to match the estimated mean height and body mass of tuberculosis patients in Canada during the relevant time period. Patient-specific data (machine settings, exposure duration, patient orientation) used during individual fluoroscopy or radiography exams were not recorded. Doses to patients were based on parameter values inferred from interviews with 91 physicians practicing at the time, historical literature, and estimated number of procedures from patient records. FLUXOR uses probability distributions to represent the uncertainty in the unknown true, average value of each dosimetry parameter. Uncertainties were shared across all patients within specific subgroups of the cohort, defined by age at treatment, sex, type of procedure, time period of exams and region (Nova Scotia or other provinces). Monte Carlo techniques were used to propagate uncertainties, by sampling alternative average values for each parameter. Alternative average doses per exam were estimated for patients in each subgroup, with the total average dose per individual determined by the number of exams received. This process was repeated to produce alternative cohort vectors of average organ doses per patient. This article presents estimates of doses to lungs, female breast, active bone marrow and heart wall. Means and 95% confidence intervals (CI) of average organ doses across all 63,715 patients were 320 (160, 560) mGy to lungs, 250 (120, 450) mGy to female breast, 190 (100, 340) mGy to heart wall and 92 (47, 160) mGy to active bone marrow. Approximately 60% of all patients had average doses to the four studied organs of less than 10 mGy, 10% received between 10 and 100 mGy, 25% between 100 and 1,000 mGy, and 5% above 1,000 mGy. Pneumothorax was the medical procedure that accounted for the largest contribution to cohort average doses. The major contributors to uncertainty in estimated doses per procedure for the four organs of interest are the uncertainties in exposure duration, tube voltage, tube output, and patient orientation relative to the X-ray tube, with the uncertainty in exposure duration being most often the dominant source. Uncertainty in patient orientation was important for doses to female breast, and, to a lesser degree, for doses to heart wall. The uncertainty in number of exams was an important contributor to uncertainty for ∼30% of patients. The estimated organ doses and their uncertainties will be used for analyses of incidence and mortality of cancer and non-cancer diseases. The CFCS cohort is an important addition to existing radio-epidemiological cohorts, given the moderate-to-high doses received fractionated over several years, the type of irradiation (external irradiation only), radiation type (X rays only), a balanced combination of both genders and inclusion of people of all ages.