General overview

Pseudomonas aeruginosa is a gram-negative bacterium capable of causing keratitis, bacteremia, and acute and chronic respiratory infections[1]. P. aeruginosa is more commonly associated with nosocomial, ventilator-associated, and community-acquired pneumonia, however it may also cause lung infections in immunocomponent hosts[2]. The bacteria develop biofilms and has been known to colonize moist, humid environments within engineered water systems.

To better understand the risk of infection associated with the inhalation route of exposure, Dean and Mitchell (2020)[3] fit traditional and nontraditional microbial dose response models to animal data evaluating the course of P. aeruginosa pneumonia. Ojielo et al. (2003)[4] inoculated groups of 10 wild-type, specific pathogen-free B6D2F1/J mice with seven graded doses of P. aeruginosa. The endpoint evaluated was death. Notably, Dean and Mitchell (2020)[3] also identified 19 additional point estimates of the dose associated with a death response from a range of studies. These 19 point estimates were used to benchmark the dose response models fit to the Ojielo et al. (2003)[4] data and support the utility of the mulit-hit model for this pathogen and exposure scenario. A more nuanced discussion of the single-hit versus cooperativity theory for P. aeruginosa can be found in the published article: https://doi.org/10.1016/j.mran.2020.100115

Recommended Model

Neither the exponential or beta-Poisson model provided a good fit to the Ojielo et al. (2003)[4] data. The exponential model faciliates a steeper curve than the beta-Poisson, but the Ojileo et al. (2003)[4] data and the additional dose response point estimates identifed in the literature suggest a more extreme slope may be a feature of P. aeruginosa dose response relationships for the inhalation route of exposure. The multi-hit model is represented by the incomplete gamma function, and has two parameters: the probability that pathogens survives to initiate an infection, k, and the number of organisms needed to initiate an infection, kmin. The multi-hit model did provide a good fit to the Ojielo et al. (2003)[4] data, with median estimates of 4.12 x 10-6 and 11 for k and kmin, respectively. 

Although the multi-hit model provided a good fit to the data, as opposed to the traditionally applied microbial dose response models, it is important to acknowledge that there has historcially been more evidence for the biological plausiblity of the single-hit models and the multi-hit model estimates de minimis risk at low doses. Any application of this dose response model should explore the impact of model selection on risk estimates, and in lieu of additional primary datasets, users are encouraged to conduct uncertainty analyses with the exponential and multi-hit models presented, the distributions explored with the hierarchical analysis, and surrogate models for similar pathogens and the inhalation exposure route. 


References

ID Exposure Route # of Doses Agent Strain Dose Units Host type Μodel LD50/ID50 Optimized parameters Response type Reference
Ojielo2003 intratracheal 7 PA01 CFU C57BL/6 mice exponential
k = 3.22E-7

death Ojielo, Charles I., et al. "Defective Phagocytosis and Clearance of Pseudomonas aeruginosa in the Lung Following Bone Marrow Transplantation." The Journal of Immunology. 171.8 (2003): 4416-4424.
intratracheal 1 TBCF10839 CFU C57BL/6 mice exponential
k = 2.8*10^-8

death Wölbeling, Florian, et al. "Lung function and inflammation during murine Pseudomonas aeruginosa airway infection." Immunobiology. 216 (2011): 901-908.
Experiment ID:
Ojielo2003
# of Doses:
7
Agent Strain:
PA01
Dose Units:
CFU
Host type:
C57BL/6 mice
Μodel:
exponential
Optimized parameters:
k = 3.22E-7


Reference:
Table 1: Goodness of fit
Model Deviance Δ Degrees 
of freedom
χ20.95,1 
p-value
χ20.95,m-k 
p-value
Exponential 16.88 0.001 6 3.84 
0.98
12.6
0.009
Beta Poisson 16.88 5 11.1
0.004
Neither model provides a good fit to the data

 

Table 2: Multi-hit dose response statistics
Model Deviance Δ DF

χ20.95,df

χ20.95,1

Good fit? Parameters LD50
Multi-hit 1.09 15.69 5 11.1 3.84 Yes k = 4.12 × 10−6 kmin=11 2,588,047

 

Fig. 1. The multi-hit dose response model with 95% and 99% confidence bands

 

Fig 2. Histogram of the k parameter estimates after bootstrapping

 

Experiment ID:
# of Doses:
1
Agent Strain:
TBCF10839
Dose Units:
CFU
Host type:
C57BL/6 mice
Μodel:
exponential
Optimized parameters:
k = 2.8*10^-8


Reference: