Characterizing the Risk of Infection from Mycobacterium tuberculosis in Commercial Passenger Aircraft Using Quantitative Microbial Risk Assessment

Quantitative microbial risk assessment was used to predict the likelihood and spatial organization of Mycobacterium tuberculosis ( Mtb ) transmission in a commercial aircraft. Passenger exposure was predicted via a multizone Markov model in four scenarios: seated or moving infectious passengers and with or without filtration of recirculated cabin air. The traditional exponential ( k  = 1) and a new exponential ( k  = 0.0218) dose-response function were used to compute infection risk. Emission variability was included by Monte Carlo simulation. Infection risks were higher nearer and aft of the source; steady state airborne concentration levels were not attained. Expected incidence was low to moderate, with the central 95% ranging from 10⁻⁶ to 10⁻¹ per 169 passengers in the four scenarios. Emission rates used were low compared to measurements from active TB patients in wards, thus a "superspreader" emitting 44 quanta/h could produce 6.2 cases or more under these scenarios. Use of respiratory protection by the infectious source and/or susceptible passengers reduced infection incidence up to one order of magnitude.     

Dose Response for Infection by Escherichia coli O157:H7 from Outbreak Data

In 1996, an outbreak of E. coli O157:H7-associated illness occurred in an elementary school in Japan. This outbreak has been studied in unusual detail, making this an important case for quantitative risk assessment. The availability of stored samples of the contaminated food allowed reliable estimation of exposure to the pathogens. Collection of fecal samples allowed assessment of the numbers infected, including asymptomatic cases. Comparison to other published dose-response studies for E. coli O157:H7 show that the strain that caused the outbreak studied here must have been considerably more infectious. We use this well-documented incident as an example to demonstrate how such information on the response to a single dose can be used for dose-response assessment. In particular, we demonstrate how the high infectivity limits the uncertainty in the low-dose region.     

A Probabilistic Transmission Dynamic Model to Assess Indoor Airborne Infection Risks

The purpose of this article is to quantify the public health risk associated with inhalation of indoor airborne infection based on a probabilistic transmission dynamic modeling approach. We used the Wells-Riley mathematical model to estimate (1) the CO2 exposure concentrations in indoor environments where cases of inhalation airborne infection occurred based on reported epidemiological data and epidemic curves for influenza and severe acute respiratory syndrome (SARS), (2) the basic reproductive number, R0 (i.e., expected number of secondary cases on the introduction of a single infected individual in a completely susceptible population) and its variability in a shared indoor airspace, and (3) the risk for infection in various scenarios of exposure in a susceptible population for a range of R0. We also employ a standard susceptible-infectious-recovered (SIR) structure to relate Wells-Riley model derived R0 to a transmission parameter to implicate the relationships between indoor carbon dioxide concentration and contact rate. We estimate that a single case of SARS will infect 2.6 secondary cases on average in a population from nosocomial transmission, whereas less than 1 secondary infection was generated per case among school children. We also obtained an estimate of the basic reproductive number for influenza in a commercial airliner: the median value is 10.4. We suggest that improving the building air cleaning rate to lower the critical rebreathed fraction of indoor air can decrease transmission rate. Here, we show that virulence of the organism factors, infectious quantum generation rates (quanta/s by an infected person), and host factors determine the risk for inhalation of indoor airborne infection.     

A Multicompartment SIS Stochastic Model with Zonal Ventilation for the Spread of Nosocomial Infections: Detection,Outbreak Management,and Infection Control

In this work, we study the environmental and operational factors that influence airborne transmission of nosocomial infections. We link a deterministic zonal ventilation model for the airborne distribution of infectious material in a hospital ward, with a Markovian multicompartment SIS model for the infection of individuals within this ward, in order to conduct a parametric study on ventilation rates and their effect on the epidemic dynamics. Our stochastic model includes arrival and discharge of patients, as well as the detection of the outbreak by screening events or due to symptoms being shown by infective patients. For each ventilation setting, we measure the infectious potential of a nosocomial outbreak in the hospital ward by means of a summary statistic: the number of infections occurred within the hospital ward until end or declaration of the outbreak. We analytically compute the distribution of this summary statistic, and carry out local and global sensitivity analysis in order to identify the particular characteristics of each ventilation regime with the largest impact on the epidemic spread. Our results show that ward ventilation can have a significant impact on the infection spread, especially under slow detection scenarios or in overoccupied wards, and that decreasing the infection risk for the whole hospital ward might increase the risk in specific areas of the health‐care facility. Moreover, the location of the initial infective individual and the protocol in place for outbreak declaration both form an interplay with ventilation of the ward.     

A Risk Analysis for Airborne Pathogens with Low Infectious Doses: Application to Respirator Selection Against Coccidioides immitis Spores

Probability models incorporating a deterministic versus stochastic infectious dose are described for estimating infection risk due to airborne pathogens that infect at low doses. Such pathogens can be occupational hazards or candidate agents for bioterrorism. Inputs include parameters for the infectious dose model, distribution parameters for ambient pathogen concentrations, the breathing rate, the duration of an exposure period, the anticipated number of exposure periods, and, if a respirator device is used, distribution parameters for respirator penetration values. Application of the models is illustrated with a hypothetical scenario involving exposure to Coccidioides immitis, a fungus present in soil in areas of the southwestern United States Inhaling C. immitis spores causes a respiratory tract infection and is a recognized occupational hazard in jobs involving soil dust exposure in endemic areas An uncertainty analysis is applied to risk estimation in the context of selecting respiratory protection with a desired degree of efficacy.     

Risk Assessment of Virus in Drinking Water

The reevaluation of drinking water treatment practices in a desire to minimize the formation of disinfection byproducts while assuring minimum levels of public health protection against infectious organisms has caused it to become necessary to consider the problem of estimation of risks posed from exposure to low levels of microorganisms, such as virus or protozoans, found in treated drinking water. This paper outlines a methodology based on risk assessment principles to approach the problem. The methodology is validated by comparison with results obtained in a prospective epidemiological study. It is feasible to produce both point and interval estimates of infection, illness and perhaps mortality by this methodology. Areas of uncertainty which require future data are indicated.     

Time-Dose-Response Models for Microbial Risk Assessment

While microbial risk assessment (MRA) has been used for over 25 years, traditional dose-response analysis has only predicted the overall risk of adverse consequences from exposure to a given dose. An important issue for consequence assessment from bioterrorist and other microbiological exposure is the distribution of cases over time due to the initial exposure. In this study, the classical exponential and beta-Poisson dose-response models were modified to include exponential-power dependency of time post inoculation (TPI) or its simplified form, exponential-reciprocal dependency of TPI, to quantify the time of onset of an effect presumably associated with the kinetics of in vivo bacterial growth. Using the maximum likelihood estimation approach, the resulting time-dose-response models were found capable of providing statistically acceptable fits to all tested pooled animal survival dose-response data. These new models can consequently describe the development of animal infectious response over time and represent observed responses fairly accurately. This is the first study showing that a time-dose-response model can be developed for describing infections initiated by various pathogens. It provides an advanced approach for future MRA frameworks.     

Effect of detection methods on risk estimates of exposure to biosolids-associated human enteric viruses

This study illustrates the effect of virus detection methods on estimates of risks of infection of biosolids-associated viruses for occupational workers and residential population during a hypothetical exposure of biosolids. Five gastroenteritis-associated human enteric viruses--enteroviruses (echovirus-12, enteroviruse types 68-71), adenoviruses, rotaviruses, and noroviruses genotype--I-were considered to represent human enteric viruses for risk estimation purposes. Ingested viral doses were calculated using literature-reported total infectious virus concentrations (based on BGM and A549 cell lines) and genome copies (GCs) in Michigan dewatered and class B biosolids. Cell-line-based infectivity parameters (i.e., ratio of total infectious virus concentration to GCs) were developed for different viruses in biosolids to use GCs for calculating ingested viral dose, addressing the issue of integration of molecular methods with biosolids-based virus risk assessment. Use of virus concentrations from molecular methods (with and without using cell-line-based infectivity parameter) resulted in higher risk estimates than culture methods, indicating the effect of the virus detection method on risk estimates. Further, use of virus concentrations from A549 cell lines resulted in higher risk estimates compared to those from BGM cell lines, suggesting the need for a proper choice of cell lines in determining infectious viral dose. The Monte Carlo uncertainty analyses of estimates for risk of infection due to enteroviruses showed that enteroviruses concentration was the most important parameter influencing risk estimates, indicating the need for reducing associated uncertainty. More work is required to develop cell-line-based infectivity parameters for different virus concentration levels and sample matrix types using a cut-off-based approach.     

Application of a Risk Analysis Tool to Middle East Respiratory Syndrome Coronavirus (MERS-CoV) Outbreak in Saudi Arabia

The Grunow–Finke assessment tool (GFT) is an accepted scoring system for determining likelihood of an outbreak being unnatural in origin. Considering its high specificity but low sensitivity, a modified Grunow–Finke tool (mGFT) has been developed with improved sensitivity. The mGFT has been validated against some past disease outbreaks, but it has not been applied to ongoing outbreaks. This study is aimed to score the outbreak of Middle East respiratory syndrome coronavirus (MERS-CoV) in Saudi Arabia using both the original GFT and mGFT. The publicly available data on human cases of MERS-CoV infections reported in Saudi Arabia (2012–2018) were sourced from the FluTrackers, World Health Organization, Saudi Ministry of Health, and published literature associated with MERS outbreaks investigations. The risk assessment of MERS-CoV in Saudi Arabia was analyzed using the original GFT and mGFT criteria, algorithms, and thresholds. The scoring points for each criterion were determined by three researchers to minimize the subjectivity. The results showed 40 points of total possible 54 points using the original GFT (likelihood: 74%), and 40 points of a total possible 60 points (likelihood: 67%) using the mGFT, both tools indicating a high likelihood that human MERS-CoV in Saudi Arabia is unnatural in origin. The findings simply flag unusual patterns in this outbreak, but do not prove unnatural etiology. Proof of bioattacks can only be obtained by law enforcement and intelligence agencies. This study demonstrated the value and flexibility of the mGFT in assessing and predicting the risk for an ongoing outbreak with simple criteria.     

A Risk Assessment Scheme of Infection Transmission Indoors Incorporating the Impact of Resuspension

A new risk assessment scheme was developed to quantify the impact of resuspension to infection transmission indoors. Airborne and surface pathogenic particle concentration models including the effect of two major resuspension scenarios (airflow‐induced particle resuspension [AIPR] and walking‐induced particle resuspension [WIPR]) were derived based on two‐compartment mass balance models and validated against experimental data found in the literature. The inhalation exposure to pathogenic particles was estimated using the derived airborne concentration model, and subsequently incorporated into a dose‐response model to assess the infection risk. Using the proposed risk assessment scheme, the influences of resuspension towards indoor infection transmission were examined by two hypothetical case studies. In the case of AIPR, the infection risk increased from 0 to 0.54 during 0–0.5 hours and from 0.54 to 0.57 during 0.5–4 hours. In the case of WIPR, the infection risk increased from 0 to 0.87 during 0–0.5 hours and from 0.87 to 1 during 0.5–4 hours. Sensitivity analysis was conducted based on the design‐of‐experiments method and showed that the factors that are related to the inspiratory rate of viable pathogens and pathogen virulence have the most significant effect on the infection probability under the occurrence of AIPR and WIPR. The risk assessment scheme could serve as an effective tool for the risk assessment of infection transmission indoors.     

Combining QMRA and Epidemiology to Estimate Campylobacteriosis Incidence

The disease burden of pathogens as estimated by QMRA (quantitative microbial risk assessment) and EA (epidemiological analysis) often differs considerably. This is an unsatisfactory situation for policymakers and scientists. We explored methods to obtain a unified estimate using campylobacteriosis in the Netherlands as an example, where previous work resulted in estimates of 4.9 million (QMRA) and 90,600 (EA) cases per year. Using the maximum likelihood approach and considering EA the gold standard, the QMRA model could produce the original EA estimate by adjusting mainly the dose‐infection relationship. Considering QMRA the gold standard, the EA model could produce the original QMRA estimate by adjusting mainly the probability that a gastroenteritis case is caused by Campylobacter. A joint analysis of QMRA and EA data and models assuming identical outcomes, using a frequentist or Bayesian approach (using vague priors), resulted in estimates of 102,000 or 123,000 campylobacteriosis cases per year, respectively. These were close to the original EA estimate, and this will be related to the dissimilarity in data availability. The Bayesian approach further showed that attenuating the condition of equal outcomes immediately resulted in very different estimates of the number of campylobacteriosis cases per year and that using more informative priors had little effect on the results. In conclusion, EA was dominant in estimating the burden of campylobacteriosis in the Netherlands. However, it must be noted that only statistical uncertainties were taken into account here. Taking all, usually difficult to quantify, uncertainties into account might lead to a different conclusion.     

Knowledge-Based Risk Assessment Under Uncertainty for Species Invasion

Management of invasive species depends on developing prevention and control strategies through comprehensive risk assessment frameworks that need a thorough analysis of exposure to invasive species. However, accurate exposure analysis of invasive species can be a daunting task because of the inherent uncertainty in invasion processes. Risk assessment of invasive species under uncertainty requires potential integration of expert judgment with empirical information, which often can be incomplete, imprecise, and fragmentary. The representation of knowledge in classical risk models depends on the formulation of a precise probabilistic value or well-defined joint distribution of unknown parameters. However, expert knowledge and judgments are often represented in value-laden terms or preference-ordered criteria. We offer a novel approach to risk assessment by using a dominance-based rough set approach to account for preference order in the domains of attributes in the set of risk classes. The model is illustrated with an example showing how a knowledge-centric risk model can be integrated with the dominance-based principle of rough set to derive minimal covering "if ... , then...," decision rules to reason over a set of possible invasion scenarios. The inconsistency and ambiguity in the data set is modeled using the rough set concept of boundary region adjoining lower and upper approximation of risk classes. Finally, we present an extension of rough set to evidence a theoretic interpretation of risk measures of invasive species in a spatial context. In this approach, the multispecies interactions in an invasion risk are approximated with imprecise probability measures through a combination of spatial neighborhood information of risk estimation in terms of belief and plausibility.     

Insight into “Calculated Risk”: An Application to the Prioritization of Emerging Infectious Diseases for Blood Transfusion Safety

Increasing identification of transmissions of emerging infectious diseases (EIDs) by blood transfusion raised the question which of these EIDs poses the highest risk to blood safety. For a number of the EIDs that are perceived to be a threat to blood safety, evidence on actual disease or transmission characteristics is lacking, which might render measures against such EIDs disputable. On the other hand, the fact that we call them “emerging” implies almost by definition that we are uncertain about at least some of their characteristics. So what is the relative importance of various disease and transmission characteristics, and how are these influenced by the degree of uncertainty associated with their actual values? We identified the likelihood of transmission by blood transfusion, the presence of an asymptomatic phase of infection, prevalence of infection, and the disease impact as the main characteristics of the perceived risk of disease transmission by blood transfusion. A group of experts in the field of infectious diseases and blood transfusion ranked sets of (hypothetical) diseases with varying degrees of uncertainty associated with their disease characteristics, and used probabilistic inversion to obtain probability distributions for the weight of each of these risk characteristics. These distribution weights can be used to rank both existing and newly emerging infectious diseases with (partially) known characteristics. Analyses show that in case there is a lack of data concerning disease characteristics, it is the uncertainty concerning the asymptomatic phase and the disease impact that are the most important drivers of the perceived risk. On the other hand, if disease characteristics are well established, it is the prevalence of infection and the transmissibility of the disease by blood transfusion that will drive the perceived risk. The risk prioritization model derived provides an easy to obtain and rational expert assessment of the relative importance of an (emerging) infectious disease, requiring only a limited amount of information. Such a model might be used to justify a rational and proportional response to an emerging infectious disease, especially in situations where little or no specific information is available.     

Development of Dose‐Response Models to Predict the Relationship for Human Toxoplasma gondii Infection Associated with Meat Consumption

Toxoplasma gondii is a protozoan parasite that is responsible for approximately 24% of deaths attributed to foodborne pathogens in the United States. It is thought that a substantial portion of human T. gondii infections is acquired through the consumption of meats. The dose‐response relationship for human exposures to T. gondii‐infected meat is unknown because no human data are available. The goal of this study was to develop and validate dose‐response models based on animal studies, and to compute scaling factors so that animal‐derived models can predict T. gondii infection in humans. Relevant studies in literature were collected and appropriate studies were selected based on animal species, stage, genotype of T. gondii, and route of infection. Data were pooled and fitted to four sigmoidal‐shaped mathematical models, and model parameters were estimated using maximum likelihood estimation. Data from a mouse study were selected to develop the dose‐response relationship. Exponential and beta‐Poisson models, which predicted similar responses, were selected as reasonable dose‐response models based on their simplicity, biological plausibility, and goodness fit. A confidence interval of the parameter was determined by constructing 10,000 bootstrap samples. Scaling factors were computed by matching the predicted infection cases with the epidemiological data. Mouse‐derived models were validated against data for the dose‐infection relationship in rats. A human dose‐response model was developed as P (d) = 1–exp (–0.0015 × 0.005 × d) or P (d) = 1–(1 + d × 0.003 / 582.414)^?1.479. Both models predict the human response after consuming T. gondii‐infected meats, and provide an enhanced risk characterization in a quantitative microbial risk assessment model for this pathogen.     

Incorporating Time‐Dose‐Response into Legionella Outbreak Models

A novel method was used to incorporate in vivo host–pathogen dynamics into a new robust outbreak model for legionellosis. Dose‐response and time‐dose‐response (TDR) models were generated for Legionella longbeachae exposure to mice via the intratracheal route using a maximum likelihood estimation approach. The best‐fit TDR model was then incorporated into two L. pneumophila outbreak models: an outbreak that occurred at a spa in Japan, and one that occurred in a Melbourne aquarium. The best‐fit TDR from the murine dosing study was the beta‐Poisson with exponential‐reciprocal dependency model, which had a minimized deviance of 32.9. This model was tested against other incubation distributions in the Japan outbreak, and performed consistently well, with reported deviances ranging from 32 to 35. In the case of the Melbourne outbreak, the exponential model with exponential dependency was tested against non‐time‐dependent distributions to explore the performance of the time‐dependent model with the lowest number of parameters. This model reported low minimized deviances around 8 for the Weibull, gamma, and lognormal exposure distribution cases. This work shows that the incorporation of a time factor into outbreak distributions provides models with acceptable fits that can provide insight into the in vivo dynamics of the host‐pathogen system.     

A Simulation Study of Quantitative Risk Assessment for Bivariate Continuous Outcomes

The neurotoxic effects of chemical agents are often investigated in controlled studies on rodents, with binary and continuous multiple endpoints routinely collected. One goal is to conduct quantitative risk assessment to determine safe dose levels. Yu and Catalano (2005) describe a method for quantitative risk assessment for bivariate continuous outcomes by extending a univariate method of percentile regression. The model is likelihood based and allows for separate dose‐response models for each outcome while accounting for the bivariate correlation. The approach to benchmark dose (BMD) estimation is analogous to that for quantal data without having to specify arbitrary cutoff values. In this article, we evaluate the behavior of the BMD relative to background rates, sample size, level of bivariate correlation, dose‐response trend, and distributional assumptions. Using simulations, we explore the effects of these factors on the resulting BMD and BMDL distributions. In addition, we illustrate our method with data from a neurotoxicity study of parathion exposure in rats.     

Quantitative Assessment of the Risk of Lung Cancer Associated with Occupational Exposure to Refractory Ceramic Fibers

We present the results of a quantitative assessment of the lung cancer risk associated with occupational exposure to refractory ceramic fibers (RCF). The primary sources of data for our risk assessment were two long-term oncogenicity studies in male Fischer rats conducted to assess the potential pathogenic effects associated with prolonged inhalation of RCF. An interesting feature of the data was the availability of the temporal profile of fiber burden in the lungs of experimental animals. Because of this information, we were able to conduct both exposure–response and dose–response analyses. Our risk assessment was conducted within the framework of a biologically based model for carcinogenesis, the two-stage clonal expansion model, which allows for the explicit incorporation of the concepts of initiation and promotion in the analyses. We found that a model positing that RCF was an initiator had the highest likelihood. We proposed an approach based on biological considerations for the extrapolation of risk to humans. This approach requires estimation of human lung burdens for specific exposure scenarios, which we did by using an extension of a model due to Yu. Our approach acknowledges that the risk associated with exposure to RCF depends on exposure to other lung carcinogens. We present estimates of risk in two populations: (1) a population of nonsmokers and (2) an occupational cohort of steelworkers not exposed to coke oven emissions, a mixed population that includes both smokers and nonsmokers.     

Comparison of Six Dose-Response Models for Use with Food-Borne Pathogens

Food-related illness in the United States is estimated to affect over six million people per year and cost the economy several billion dollars. These illnesses and costs could be reduced if minimum infectious doses were established and used as the basis of regulations and monitoring. However, standard methodologies for dose-response assessment are not yet formulated for microbial risk assessment. The objective of this study was to compare dose-response models for food-borne pathogens and determine which models were most appropriate for a range of pathogens. The statistical models proposed in the literature and chosen for comparison purposes were log-normal, log-logistic, exponential, -Poisson and Weibull-Gamma. These were fit to four data sets also taken from published literature, Shigella flexneri, Shigella dysenteriae,Campylobacter jejuni, and Salmonella typhosa, using the method of maximum likelihood. The Weibull-gamma, the only model with three parameters, was also the only model capable of fitting all the data sets examined using the maximum likelihood estimation for comparisons. Infectious doses were also calculated using each model. Within any given data set, the infectious dose estimated to affect one percent of the population ranged from one order of magnitude to as much as nine orders of magnitude, illustrating the differences in extrapolation of the dose response models. More data are needed to compare models and examine extrapolation from high to low doses for food-borne pathogens.     

A Tiered Approach for Risk‐Benefit Assessment of Foods

Risk‐benefit analyses are introduced as a new paradigm for old problems. However, in many cases it is not always necessary to perform a full comprehensive and expensive quantitative risk‐benefit assessment to solve the problem, nor is it always possible, given the lack of required date. The choice to continue from a more qualitative to a full quantitative risk‐benefit assessment can be made using a tiered approach. In this article, this tiered approach for risk‐benefit assessment will be addressed using a decision tree. The tiered approach described uses the same four steps as the risk assessment paradigm: hazard and benefit identification, hazard and benefit characterization, exposure assessment, and risk‐benefit characterization, albeit in a different order. For the purpose of this approach, the exposure assessment has been moved upward and the dose‐response modeling (part of hazard and benefit characterization) is moved to a later stage. The decision tree includes several stop moments, depending on the situation where the gathered information is sufficient to answer the initial risk‐benefit question. The approach has been tested for two food ingredients. The decision tree presented in this article is useful to assist on a case‐by‐case basis a risk‐benefit assessor and policymaker in making informed choices when to stop or continue with a risk‐benefit assessment.     

Quantifying Water Pathogen Risk in an Epidemiological Framework

Traditionally, microbial risk assessors have used point estimates to evaluate the probability that an individual will become infected. We developed a quantitative approach that shifts the risk characterization perspective from point estimate to distributional estimate, and from individual to population. To this end, we first designed and implemented a dynamic model that tracks traditional epidemiological variables such as the number of susceptible, infected, diseased, and immune, and environmental variables such as pathogen density. Second, we used a simulation methodology that explicitly acknowledges the uncertainty and variability associated with the data. Specifically, the approach consists of assigning probability distributions to each parameter, sampling from these distributions for Monte Carlo simulations, and using a binary classification to assess the output of each simulation. A case study is presented that explores the uncertainties in assessing the risk of giardiasis when swimming in a recreational impoundment using reclaimed water. Using literature-based information to assign parameters ranges, our analysis demonstrated that the parameter describing the shedding of pathogens by infected swimmers was the factor that contributed most to the uncertainty in risk. The importance of other parameters was dependent on reducing the a priori range of this shedding parameter. By constraining the shedding parameter to its lower subrange, treatment efficiency was the parameter most important in predicting whether a simulation resulted in prevalences above or below non outbreak levels. Whereas parameters associated with human exposure were important when the shedding parameter was constrained to a higher subrange. This Monte Carlo simulation technique identified conditions in which outbreaks and/or nonoutbreaks are likely and identified the parameters that most contributed to the uncertainty associated with a risk prediction.