Influenza Infection Risk and Predominate Exposure Route: Uncertainty Analysis

An effective nonpharmaceutical intervention for influenza interrupts an exposure route that contributes significantly to infection risk. Herein, we use uncertainty analysis (point‐interval method) and Monte Carlo simulation to explore the magnitude of infection risk and predominant route of exposure. We utilized a previously published mathematical model of a susceptible person attending a bed‐ridden infectious person. Infection risk is sensitive to the magnitude of virus emission and contact rates. The contribution of droplet spray exposure to infection risk increases with cough frequency, and decreases with virus concentration in cough particles. We consider two infectivity scenarios: greater infectivity of virus deposited in the upper respiratory tract than virus inhaled in respirable aerosols, based on human studies; and equal infectivity in the two locations, based on studies in guinea pigs. Given that virus have equal probability of infection throughout the respiratory tract, the mean overall infection risk is 9.8 × 10^?2 (95th percentile 0.78). However, when virus in the upper respiratory tract is less infectious than inhaled virus, the overall infection risk is several orders of magnitude lower. In this event, inhalation is a significant exposure route. Contact transmission is important in both infectivity scenarios. The presence of virus in only respirable particles increases the mean overall infection risk by 1–3 orders of magnitude, with inhalation contributing ≥ 99% of the infection risk. The analysis indicates that reduction of uncertainties in the concentration of virus in expiratory particles of different sizes, expiratory event frequency, and infectivity at different sites in the respiratory tract will clarify the predominate exposure routes for influenza.     

The BSE Risk of Processing Meat and Bone Meal in Nonruminant Feed: A Quantitative Assessment for the Netherlands

The total ban on use of meat and bone meal (MBM) in livestock feed has been very successful in reducing bovine spongiform encephalopathy (BSE) spread, but also implies a waste of high-quality proteins resulting in economic and ecological loss. Now that the BSE epidemic is fading out, a partial lifting of the MBM ban might be considered. The objective of this study was to assess the BSE risk for the Netherlands if MBM derived from animals fit for human consumption, i.e., category 3 MBM, would be used in nonruminant feed. A stochastic simulation model was constructed that calculates (1) the probability that infectivity of undetected BSE-infected cows ends up with calves and (2) the quantity of infectivity ( Q_inf ) consumed by calves in case of such an incident. Three pathways were considered via which infectivity can reach cattle: (1) cross-contamination in the feed mill, (2) cross-contamination on the primary farm, and (3) pasture contamination. Model calculations indicate that the overall probability that infectivity ends up with calves is 3.2%. In most such incidents the Q_inf is extremely small (median = 6.5 × 10⁻¹² ID₅₀; mean = 1.8 × 10⁻⁴ ID₅₀), corresponding to an average probability of 1.3 × 10⁻⁴ that an incident results in ≥1 new BSE infections. Cross-contamination in the feed mill is the most risky pathway. Combining model results with Dutch BSE prevalence estimates for the coming years, it can be concluded that the BSE risk of using category 3 MBM derived from Dutch cattle in nonruminant feed is very low.     

Physiologically-Based Pharmacokinetic Model for Trichloroethylene Considering Enterohepatic Recirculation of Major Metabolites

Trichloroacetic acid (TCA) is major metabolite of trichloroethylene (TRI) thought to contribute to its hepatocarcinogenic effects in mice. Recent studies have shown that peak blood concentrations of TCA in rats do not occur until approximately 12 hours following an oral dose of TRI. However, blood concentrations of TRI reach maximum within an hour and are nondetectable after 2 hours.⁽¹⁾ The results of study which examined the enterohepatic recirculation (EHC) of the principle TRI metabolited⁽²⁾ was used to develop physiologically-based pharmacokinetic model for TRI, which includes enterohepatic recirculation of its metabolites. The model quantitatively predicts the uptake, distribution and elimination of TRI, trichloroethanol, trichloroethanol-glucuronide, and TCA and includes production of metabolites through the enterohepatic recirculation pathway. Physiologic parameters used in the model were obtained from the literature.^(3.4) Parameters for TRI metabolism were taken from Fisher et al.⁽⁵⁾ Other kinetic parameters were found in the literature or estimated from experimental data.⁽²⁾ The model was calibrated to data from experiments of an earlier study where TRI was orally administered⁽²⁾ Verification of the model was conducted using data on the enterohepatic recirculation of TCEOH and TCA⁽²⁾ chloral hydrate data (infusion doses) from Merdink,⁽¹⁾ and TRI data from Templin(l) and Larson and Bull.⁽¹⁾     

Development of a Quantitative Microbial Risk Assessment for Human Salmonellosis Through Household Consumption of Fresh Minced Pork Meat in Belgium

A quantitative microbial risk assessment (QMRA) according to the Codex Alimentarius Principles is conducted to evaluate the risk of human salmonellosis through household consumption of fresh minced pork meat in Belgium. The quantitative exposure assessment is carried out by building a modular risk model, called the METZOON-model, which covers the pork production from farm to fork. In the METZOON-model, the food production pathway is split up in six consecutive modules: (1) primary production, (2) transport and lairage, (3) slaughterhouse, (4) postprocessing, (5) distribution and storage, and (6) preparation and consumption. All the modules are developed to resemble as closely as possible the Belgian situation, making use of the available national data. Several statistical refinements and improved modeling techniques are proposed. The model produces highly realistic results. The baseline predicted number of annual salmonellosis cases is 20,513 ( SD 9061.45). The risk is estimated higher for the susceptible population (estimate  4.713 × 10⁻⁵; SD 1.466 × 10⁻⁵  ) compared to the normal population  (estimate 7.704 × 10⁻⁶; SD 5.414 × 10⁻⁶)  and is mainly due to undercooking and to a smaller extent to cross-contamination in the kitchen via cook's hands.     

Cardiac Sensitization Thresholds of Halon Replacement Chemicals Predicted in Humans by Physiologically-Based Pharmacokinetic Modeling

Human exposure to halons and halon replacement chemicals is often regulated on the basis of cardiac sensitization potential. The dose-response data obtained from animal testing are used to determine the no observable adverse effect level (NOAEL) and lowest observable adverse effect level (LOAEL) values. This approach alone does not provide the information necessary to evaluate the cardiac sensitization potential for the chemical of interest under a variety of exposure concentrations and durations. In order to provide a tool for decision-makers and regulators tasked with setting exposure guidelines for halon replacement chemicals, a quantitative approach was established which allowed exposures to be assessed in terms of the chemical concentrations in blood during the exposure. A physiologically-based pharmacokinetic (PBPK) model was used to simulate blood concentrations of Halon 1301 (bromotrifluoromethane, CF₃Br), HFC-125 (pentafluoroethane, CHF₂CF₃), HFC-227ea (heptafluoropropane, CF₃CHFCF₃), HCFC-123 (dichlorotrifluoroethane, CHCl₂CF₃), and CF₃I (trifluoroiodomethane) during inhalation exposures. This work demonstrates a quantitative approach for use in linking chemical inhalation exposures to the levels of chemical in blood achieved during the exposure.     

Cancer Risk Estimation for Mixtures of Coal Tars

Two-year chronic bioassays were conducted by using B6C3F1 female mice fed several concentrations of two different mixtures of coal tars from manufactured gas waste sites or benzo(a)pyrene (BaP). The purpose of the study was to obtain estimates of cancer potency of coal tar mixtures, by using conventional regulatory methods, for use in manufactured gas waste site remediation. A secondary purpose was to investigate the validity of using the concentration of a single potent carcinogen, in this case benzo(a)pyrene, to estimate the relative risk for a coal tar mixture. The study has shown that BaP dominates the cancer risk when its concentration is greater than 6,300 ppm in the coal tar mixture. In this case the most sensitive tissue site is the forestomach. Using low-dose linear extrapolation, the lifetime cancer risk for humans is estimated to be: Risk < 1.03 × 10⁻⁴ (ppm coal tar in total diet) + 240 × 10⁻⁴ (ppm BaP in total diet), based on forestomach tumors. If the BaP concentration in the coal tar mixture is less than 6,300 ppm, the more likely case, then lung tumors provide the largest estimated upper limit of risk, Risk < 2.55 × 10⁻⁴ (ppm coal tar in total diet), with no contribution of BaP to lung tumors. The upper limit of the cancer potency (slope factor) for lifetime oral exposure to benzo(a)pyrene is 1.2 × 10⁻³ per μg per kg body weight per day from this Good Laboratory Practice (GLP) study compared with the current value of 7.3 × 10⁻³ per μg per kg body weight per day listed in the U.S. EPA Integrated Risk Information System.     

On the Effect of Probability Distributions of Input Variables in Public Health Risk Assessment

A central part of probabilistic public health risk assessment is the selection of probability distributions for the uncertain input variables. In this paper, we apply the first-order reliability method (FORM)^(1–3) as a probabilistic tool to assess the effect of probability distributions of the input random variables on the probability that risk exceeds a threshold level (termed the probability of failure) and on the relevant probabilistic sensitivities. The analysis was applied to a case study given by Thompson et al. ⁽⁴⁾ on cancer risk caused by the ingestion of benzene contaminated soil. Normal, lognormal, and uniform distributions were used in the analysis. The results show that the selection of a probability distribution function for the uncertain variables in this case study had a moderate impact on the probability that values would fall above a given threshold risk when the threshold risk is at the 50th percentile of the original distribution given by Thompson et al. ⁽⁴⁾ The impact was much greater when the threshold risk level was at the 95th percentile. The impact on uncertainty sensitivity, however, showed a reversed trend, where the impact was more appreciable for the 50th percentile of the original distribution of risk given by Thompson et al. ⁴ than for the 95th percentile. Nevertheless, the choice of distribution shape did not alter the order of probabilistic sensitivity of the basic uncertain variables.     

A Trichloroethylene Risk Assessment Using a Monte Carlo Analysis of Parameter Uncertainty in Conjunction with Physiologically-Based Pharmacokinetic Modeling

A Monte Carlo simulation is incorporated into a risk assessment for trichloroethylene (TCE) using physiologically-based pharmacokinetic (PBPK) modeling coupled with the linearized multistage model to derive human carcinogenic risk extrapolations. The Monte Carlo technique incorporates physiological parameter variability to produce a statistically derived range of risk estimates which quantifies specific uncertainties associated with PBPK risk assessment approaches. Both inhalation and ingestion exposure routes are addressed. Simulated exposure scenarios were consistent with those used by the Environmental Protection Agency (EPA) in their TCE risk assessment. Mean values of physiological parameters were gathered from the literature for both mice (carcinogenic bioassay subjects) and for humans. Realistic physiological value distributions were assumed using existing data on variability. Mouse cancer bioassay data were correlated to total TCE metabolized and area-under-the-curve (blood concentration) trichloroacetic acid (TCA) as determined by a mouse PBPK model. These internal dose metrics were used in a linearized multistage model analysis to determine dose metric values corresponding to 10^-6 lifetime excess cancer risk. Using a human PBPK model, these metabolized doses were then extrapolated to equivalent human exposures (inhalation and ingestion). The Monte Carlo iterations with varying mouse and human physiological parameters produced a range of human exposure concentrations producing a 10^-6 risk.     

On the Relationship Between Carcinogenicity and Acute Toxicity

Parodi et al. ⁽¹⁾ and Zeise et al. ⁽²⁾ found a surprising statistical correlation (or association) between acute toxicity and carcinogenic potency. In order to shed light on the questions of whether or not it is a causal correlation, and whether or not it is a statistical or tautological artifact, we have compared the correlations for the NCI/NTP data set with those for chemicals not in this set. Carcinogenic potencies were taken from the Gold et al. database. We find a weak correlation with an average value of TD₅₀/LD₅₀= 0.04 for the non-NCI data set, compared with TD₅₀/LD₅₀= 0.15 for the NCI data set. We conclude that it is not easy to distinguish types of carcinogens on the basis of whether or not they are acutely toxic.     

The Possibility for "Natural" Generation of Chlorinated Organic Compounds

There is continuing concern for the exposure of persons to various chlorinated organics via the environment, for example, chlorinated disinfection byproducts in drinking water.⁽¹⁾ Some of these may be carcinogenic,⁽²⁾ although the evidence is far from strong.⁽³⁾ There is an accumulating body of evidence that one of the normal human immunological responses to foreign agents is the generation of hypochlorous acid. This evidence will be summarized. The possibility that this HOCl generated in vivo could result in the formation of organo-chlorine compounds does not appear to have been seriously considered. Based on best available information, the amount of such byproduct formation will be estimated.     

An Evaluation of Alternative Safety Criteria for Nuclear Power Plants

Safety criteria for frequency of nuclear-reactor accidents and for reactor-induced risk to individuals and to society are evaluated on the basis of their comprehensiveness, clarity, recognition of uncertainty, practicability, defensibility, simplicity, and internal consistency. Many criteria were found to be comprehensive and practicable; few completely satisfied the other evaluation standards. A consensus inferred from the most favorably evaluated criteria would allow between 1.0time10^-4 and 1.0time10^-3 core melts per reactor year, between 1.0time10^-6 and 2.0time10^-5 fatalities per reactor year per individual, and a total exposure in the United States of between 1000 and 10,000 person rems per reactor year. This consensus is consistent with the criteria proposed in NUREG0880.     

Benzene Risk Estimation Using Radiation Equivalent Coefficients

We estimated benzene risk using a novel framework of risk assessment that employed the measurement of radiation dose equivalents to benzene metabolites and a PBPK model. The highest risks for 1 μg/m³ and 3.2 mg/m³ life time exposure of benzene estimated with a linear regression were 5.4 × 10⁻⁷ and 1.3 × 10⁻³, respectively. Even though these estimates were based on in vitro chromosome aberration test data, they were about one-sixth to one-fourteenth that from other studies and represent a fairly good estimate by using radiation equivalent coefficient as an "internal standard."     

Risk Assessment for Aflatoxin: An Evaluation Based on the Multistage Model

Lifetime cancer potency of alfatoxin was assessed based on the Yeh et al. study from China in which both aflatoxin exposure and hepatitis B prevalence were measured. This study provides the best available information for estimating the carcinogenic risk posed by aflatoxin to the U.S. population. Cancer potency of aflatoxin was estimated using a biologically motivated risk assessment model. The best estimate of aflatoxin potency was 9 (mg/kg/day)⁻¹ for individuals negative for hepatitis B and 230 (mg/kg/day)⁻¹ for individuals positive for hepatitis B.     

Evacuation Risks: A Tentative Approach for Quantification

This study tries to assess the risk of deaths and injuries from motor vehicle accidents associated with an evacuation of population groups in case of nuclear plant accidents. The risk per person–km is evaluated using: (a) data from previous evacuation: information from Soufriere evacuation (Guadeloupe Island 1976) and Mississauga (1979), added to Hans and Sell's data: no road accident occurred for a sample of 1,500,000 persons; (b) national recording system for motor vehicle accident: the rates of 2.2 10 ^-8 deaths per person–km and 32 10^-8 injuries per person–km is calculated as an average. These last rates in France overestimate the number of casualties. A reasonable hypothesis is to assume that the probability of road accident occurrence follows a Poisson distribution, as these events are independent and unfrequent, as no accident was observed in a sample of 1,500,000 persons the probability is between 0 and an upper value of 0.24 10^-8 deaths per person-km and 3.29 10^-8 injuries per person–km. The average and maximum population involved within different radii around French and U.S. Nuclear power sites are taken as a sample size in order to study the total risk of deaths and injuries in the hypothesis of an evacuation being necessary to protect the populations.     

Avian Influenza Risk Perception and Live Poultry Purchase in Guangzhou, China, 2006

Human H5N1 highly pathogenic avian influenza (HPAI) infection is associated with intimate exposure to live poultry. Perceptions of risk can modify behaviors, influencing actual exposure. However, greater hazard is not necessarily followed by perception of greater risk and more precautionary behavior because self-serving cognitive biases modulate precautionary and hazardous behaviors. We examined risk perception associated with avian influenza. A total of 1,550 face-to-face within-household interviews and 1,760 telephone interviews were derived to study avian influenza risk perception and live poultry use in Guangzhou and Hong Kong, respectively. Chi-square and Mann-Whitney tests assessed bivariate associations and risk distributions, respectively, and fully adjusted multivariate logistic models determined independent risk associations. Relative to Hong Kong, perceived "generalized" risk from buying live poultry (GZ, 58%, 95% confidence interval 55–60% vs. HK, 41%, 39–43%; χ²= 86.95, df  = 1, p < 0.001) and perceived self/family risk from buying ( z  =−2.092, p  = 0.036) were higher in Guangzhou. Higher perceived "generalized" risk was associated with not buying live poultry (OR = 0.65, 0.49–0.85), consistent with the pattern seen in Hong Kong, while perceived higher self/family risk was associated with buying ("likely/very likely/certain" OR = 1.74, 1.18–2.59); no such association was seen in Hong Kong. Multivariate adjustment indicated older age was associated with buying live poultry in Guangzhou (OR = 2.91, 1.36–6.25). Guangzhou respondents perceived greater risk relative to Hong Kong. Buying live poultry was associated with perceptions of less "generalized" risk but more self/family risk. Higher generalized risk was associated with fewer live poultry purchases, suggesting generalized risk may be a useful indicator of precautionary HPAI risk behavior.     

On the Correlation Coefficient Between the TD50 and the MTD

The existence of correlation between the carcinogenic potency and the maximum tolerated dose has been the subject of many investigations in recent years. Several attempts have been made to quantify this correlation in different bioassay experiments. By using some distributional assumptions, Krewski et al .⁽¹⁾ derive an analytic expression for the coefficient of correlation between the carcinogenic potency TD₅₀ and the maximum tolerated dose. Here, we discuss the deviation that may result in using their analytical expression. By taking a more general approach we derive an expression for the correlation coefficient which includes the result of Krewski et al .⁽¹⁾ as a special case, and show that their expression may overestimate the correlation in some instances and yet underestimate the correlation in other instances. The proposed method is illustrated by application to a real dataset.     

A Physiologically Based Model for Gastrointestinal Absorption and Excretion of Chemicals Carried by Lipids

Pharmacokinetic models which incorporate independently measured anatomical characteristics and physiological flows have been widely used to predict the pharmacokinetic behavior of drugs, anesthetics, and other chemicals. Models appearing in the literature have included as many as 18,⁽¹⁾ or as few as 5 tissue compartments.⁽²⁾ With the exception of the multiple-compartment delay trains used by Bischoff⁽³⁾ to model the delays inherent to the appearance of drug metabolites in bile and segments of the intestinal lumen, very little effort has been made to incorporate the available information on gastrointestinal anatomy and physiology into more accurate gastrointestinal absorption/enterohepatic recirculation submodels. Since several authors have shown that the lymphatic system is the most significant route of absorption for highly lipophilic chemicals, we have constructed a model of gastrointestinal absorption that emphasizes chylomicron production and transport as the most significant route of absorption for nonvolatile, lipophilic chemicals. The absorption and distribution of hexachlorobenzene after intravenous vs. oral dosing are used to demonstrate features of this model.     

Dose-response assessment for influenza A virus based on data sets of infection with its live attenuated reassortants

Reported data sets on infection of volunteers challenged with wild-type influenza A virus at graded doses are few. Alternatively, we aimed at developing a dose-response assessment for this virus based on the data sets for its live attenuated reassortants. Eleven data sets for live attenuated reassortants that were fit to beta-Poisson and exponential dose-response models. Dose-response relationships for those reassortants were characterized by pooling analysis of the data sets with respect to virus subtype (H1N1 or H3N2), attenuation method (cold-adapted or avian-human gene reassortment), and human age (adults or children). Furthermore, by comparing the above data sets to a limited number of reported data sets for wild-type virus, we quantified the degree of attenuation of wild-type virus with gene reassortment and estimated its infectivity. As a result, dose-response relationships of all reassortants were best described by a beta-Poisson model. Virus subtype and human age were significant factors determining the dose-response relationship, whereas attenuation method affected only the relationship of H1N1 virus infection to adults. The data sets for H3N2 wild-type virus could be pooled with those for its reassortants on the assumption that the gene reassortment attenuates wild-type virus by at least 63 times and most likely 1,070 times. Considering this most likely degree of attenuation, 10% infectious dose of H3N2 wild-type virus for adults was estimated at 18 TCID50 (95% CI = 8.8-35 TCID50). The infectivity of wild-type H1N1 virus remains unknown as the data set pooling was unsuccessful.     

Morphometric Approaches for Evaluating Pulmonary Toxicity in Mammals: Implications for Risk Assessment

Recent advances in quantitative morphology provide all the tools necessary to obtain structural information in the lung that can be quantified and interpreted in the three-dimensional world of toxicology. Structural hierarchies of conducting airways and parenchyma of the lung provide: (1) numbers of cells per airway, lobe, or lung; (2) surface areas of cells, airways, and alveoli; (3) length of airways and vessels; and (4) volumes of cells, alveoli, airways, vessels, and individual lobes or the entire lung. Unbiased sampling of these subcompartments of the lung requires fractionation of lobes or individual airways. Individual airways of proximal and distal generations are obtained by airway microdissection along one axial pathway and comparisons made between airway generations. Vertical sections of selected airways are used to sample epithelium and interstitium. Using this unbiased approach of quantitative morphology, we have shown that inhalation of low ambient concentrations of ozone ([O₃]0.15 ppm) near or at the United States National Ambient Air Quality Standard (NAAQS) (0.12 ppm O₃) induces significant alterations in bronchiolar epithelium and interstitium in nonhuman primates but not rats. The alterations do not appear to be concentration- or time-dependent, thereby bringing into question the current NAAQS that may be at or above the threshold for distal airway injury in primates. Unbiased morphometric methods are critical in a quantitative evaluation of toxicological injury of mammalian tracheobronchial airways.     

A Method for Computing the Fraction of Attributable Risk Related to Climate Damages

The recent decision of the U.S. Supreme Court on the regulation of CO₂ emissions from new motor vehicles^{( 1 )} shows the need for a robust methodology to evaluate the fraction of attributable risk from such emissions. The methodology must enable decisionmakers to reach practically relevant conclusions on the basis of expert assessments the decisionmakers see as an expression of research in progress, rather than as knowledge consolidated beyond any reasonable doubt.^{( 2,3,4 )} This article presents such a methodology and demonstrates its use for the Alpine heat wave of 2003. In a Bayesian setting, different expert assessments on temperature trends and volatility can be formalized as probability distributions, with initial weights (priors) attached to them. By Bayesian learning, these weights can be adjusted in the light of data. The fraction of heat wave risk attributable to anthropogenic climate change can then be computed from the posterior distribution. We show that very different priors consistently lead to the result that anthropogenic climate change has contributed more than 90% to the probability of the Alpine summer heat wave in 2003. The present method can be extended to a wide range of applications where conclusions must be drawn from divergent assessments under uncertainty.