General Overview

Shigella sp. are very closely related to E. coli and probably technically belong to that species; however, because of their disease significance, the distinct names are still used.[1] The most common types in developing countries are Shigella flexneri (~60% of cases) and Shigella sonnei (~15% of cases), with S. boydiiand S. dysenteriae contributing ~6% each. However, S. dysenteriae type 1 can spread in epidemics and tends to cause more severe and complicated disease. Nearly all cases occur in developing countries, with 69% of cases in children under 5 years of age.

Shigella tends to cause dysentery (diarrhea with blood or mucus) by invading the mucosa of the colon, leading to ulceration and bleeding. It is usually self-limited after 4-7 days, but it tends to be more serious than watery diarrhea and is not effectively treated by oral rehydration solution. Pandemic waves of shigellosis have been described since the 1960s, and can cause high mortality, especially among refugees. Shigellosis is limited to human hosts. Vaccines against Shigella are currently under development.
http://www.cdc.gov/nczved/divisions/dfbmd/diseases/shigellosis/

Summary of data

Levine et al. (1973)  gave S. dysenteriae strains M 131 and A-1 in milk to healthy young adult male prisoners, and measured illness as the outcome. Four dose levels were used for strain M 131, but only 2 dose levels were used for strain A-1. Data for strain M 131  are fit well by the beta-Poisson model, although the confidence intervals are wide. However, a poor fit is obtained if the data from strain A-1  are pooled with M 131. Powell (2000)  reports a beta-Poisson model fit to the pooled data for strains M 131 and A-1.

DuPont et al. (1969 & 1972) conducted feeding studies of S. flexneri 2a (strain 2457T) in healthy male prisoners. A smaller study  used 5 doses in 31 subjects; a larger study  used 4 doses among 196 subjects. The smaller study only recorded illness as the response; however, the larger study recorded both illness and infection. The smaller study tested a wide dose range, from 104 to 108 CFU; however, the proportion ill did not change greatly over the range of 105 – 108 (68-88% ill). The beta-Poisson model fits better than the exponential model  for each of these 3 datasets. It fits well for the two datasets using illness as a response (DuPont et al. 1969 & 1972), but fits poorly for the infection response.

There have been two feeding studies  in male prisoners of S. sonnei 53G; however, they both used a dose of 500 CFU in all subjects. However, 7/20 subjects became ill in one study, and 19/38 in the other, implying that this dose is close to the ID50.

There have been many additional feeding studies evaluating attenuated strains of Shigella for use in vaccines (e.g., DuPont et al. 1972a , Levine et al. 1973 , Kotloff et al. 1995). However, these are not included here since attenuated strains are deliberately intended to be less infectious or pathogenic than wild strains, and using dose response models based on such experiments might lead to underestimation of risk.

Pooling analysis of all datasets using illness as the response could not disprove the hypothesis that the datasets could be pooled (P > 0.05). However, separate pooled models are also provided for S. dysenteriae and S. flexneri. A previously published pooling analysis  excluded one dose level (1E7 CFU) from experiment 82. Although that dose level contributed disproportionately to the deviance, it was also the dose level with the largest sample size (19 volunteers; the next largest sample size was 8 volunteers); therefore, it was retained in all analyses conducted for this chapter.

Recommended Model

Exponential and betapoisson model.jpg

ID Exposure Route # of Doses Agent Strain Dose Units Host type Μodel LD50/ID50 Optimized parameters Response type Reference
223 oral (in milk) 4.00 2a (strain 2457T) CFU human beta-Poisson 3.11E+03 a = 1.35E-01 N50 = 3.11E+03 illness
DuPont, H. L., Levine, M. M., Hornick, R. B., & Formal, S. B. (1989). Inoculum Size in Shigellosis and Implications for Expected Mode of Transmission. Journal of Infectious Diseases, 159, 6.
82 oral (in milk) 5.00 2a (strain 2457T) CFU human beta-Poisson 3.54E+04 a = 1.43E-01 N50 = 3.54E+04 illness
DuPont, H. L., Hornick, R. B., Snyder, M. J., Libonati, J. P., Formal, S. B., & Gangarosa, E. J. (1972). Immunity in Shigellosis. I. Response of Man to Attenuated Strains of Shigella. Journal of Infectious Diseases, 125, 1.
82, 223 oral (in milk) 9.00 2a (strain 2457T) CFU human beta-Poisson 3.64E+03 a = 1.17E-01 N50 = 3.64E+03 illness
DuPont, H. L., Hornick, R. B., Dawkins, A. T., Snyder, M. J., & Formal, S. B. (1969). The Response of Man to Virulent Shigella flexneri 2a. Journal of Infectious Diseases, 119, 3.
83 oral (in milk) 4.00 2a (strain 2457T) CFU human beta-Poisson 1.48E+03 a = 2.65E-01 N50 = 1.48E+03 positive stool isolation
DuPont, H. L., Hornick, R. B., Snyder, M. J., Libonati, J. P., Formal, S. B., & Gangarosa, E. J. (1972). Immunity in Shigellosis. I. Response of Man to Attenuated Strains of Shigella. Journal of Infectious Diseases, 125, 1.
Exposure Route:
oral (in milk)
# of Doses:
4.00
Agent Strain:
2a (strain 2457T)
Dose Units:
CFU
Host type:
human
Μodel:
beta-Poisson
LD50/ID50:
3.11E+03
Optimized parameters: a = 1.35E-01 N50 = 3.11E+03
Response type:
illness

Model data for Shigella flexneri (2a, strain 2457T) in the human 
Dose Illness Not illness Total
180 9 27 36
5000 28 21 49
1E+04 52 36 88
1E+05 14 10 24

 

Goodness of fit and model selection
Model Deviance Δ Degrees 
of freedom
χ20.95,1 
p-value
χ20.95,m-k 
p-value
Exponential 154 153 3 3.84 
0
7.81 
0
Beta Poisson 1.69 2 5.99 
0.429
Beta-Poisson fits better than exponential; cannot reject good fit for beta-Poisson.

 

Optimized parameters for the beta-Poisson model, from 10000 bootstrap iterations
Parameter MLE estimate Percentiles
0.5% 2.5% 5% 95% 97.5% 99.5%
α 1.35E-01 9.86E-04 1.22E-02 6.85E-02 2.06E-01 2.22E-01 2.58E-01
N50 3.11E+03 5.09E+01 4.31E+02 9.64E+02 7.49E+03 9.14E+03 1.89E+04

 

Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters
Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters
beta Poisson model plot, with confidence bounds around optimized model
beta Poisson model plot, with confidence bounds around optimized model
Exposure Route:
oral (in milk)
# of Doses:
5.00
Agent Strain:
2a (strain 2457T)
Dose Units:
CFU
Host type:
human
Μodel:
beta-Poisson
LD50/ID50:
3.54E+04
Optimized parameters: a = 1.43E-01 N50 = 3.54E+04
Response type:
illness

Model data for Shigella flexneri (2a, strain 2457T) in the human 
Dose Illness Not illness Total
1E+04 1 3 4
1E+05 3 1 4
1E+06 7 1 8
1E+07 13 6 19
1E+08 7 1 8

 

Goodness of fit and model selection
Model Deviance Δ Degrees 
of freedom
χ20.95,1 
p-value
χ20.95,m-k 
p-value
Exponential 73.6 70.2 4 3.84 
0
9.49 
3.89e-15
Beta Poisson 3.44 3 7.81 
0.329
Beta-Poisson fits better than exponential; cannot reject good fit for beta-Poisson.

 

Optimized parameters for the beta-Poisson model, from 10000 bootstrap iterations
Parameter MLE estimate Percentiles
0.5% 2.5% 5% 95% 97.5% 99.5%
α 1.43E-01 1.32E-02 2.12E-02 2.63E-02 2.90E-01 3.32E-01 4.27E-01
N50 3.54E+04 1.29E-11 8.16E-07 3.11E-04 3.34E+05 4.65E+05 9.68E+05
Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters.
Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters.

 

beta Poisson model plot, with confidence bounds around optimized model
beta Poisson model plot, with confidence bounds around optimized model
Exposure Route:
oral (in milk)
# of Doses:
9.00
Agent Strain:
2a (strain 2457T)
Dose Units:
CFU
Host type:
human
Μodel:
beta-Poisson
LD50/ID50:
3.64E+03
Optimized parameters: a = 1.17E-01 N50 = 3.64E+03
Response type:
illness

Model data for Shigella flexneri (2a, strain 2457T) in the human 
Dose Illness Not illness Total
180 9 27 36
5000 28 21 49
1E+04 1 3 4
1E+04 52 36 88
1E+05 3 1 4
1E+05 14 10 24
1E+06 7 1 8
1E+07 13 6 19
1E+08 7 1 8

 

Goodness of fit and model selection
Model Deviance Δ Degrees 
of freedom
χ20.95,1 
p-value
χ20.95,m-k 
p-value
Exponential 677 671 8 3.84 
0
15.5 
0
Beta Poisson 6.04 7 14.1 
0.535
Beta-Poisson fits better than exponential; cannot reject good fit for beta-Poisson.

 

Optimized parameters for the beta-Poisson model, from 10000 bootstrap iterations
Parameter MLE estimate Percentiles
0.5% 2.5% 5% 95% 97.5% 99.5%
α 1.17E-01 4.86E-02 6.43E-02 7.25E-02 1.68E-01 1.78E-01 2.02E-01
N50 3.64E+03 3.32E+02 8.92E+02 1.22E+03 9.16E+03 1.12E+04 1.73E+04

 

Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters.
Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters.
beta Poisson model plot, with confidence bounds around optimized model
beta Poisson model plot, with confidence bounds around optimized model
Highest quality
Exposure Route:
oral (in milk)
# of Doses:
4.00
Agent Strain:
2a (strain 2457T)
Dose Units:
CFU
Host type:
human
Μodel:
beta-Poisson
LD50/ID50:
1.48E+03
Optimized parameters: a = 2.65E-01 N50 = 1.48E+03
Response type:
positive stool isolation

Model data for Shigella flexneri (2a, strain 2457T) in the human 
Dose Positive stool isolation No positive stool isolation Total
180 6 30 36
5000 33 16 49
1E+04 66 21 87
1E+05 15 9 24

 

Goodness of fit and model selection
Model Deviance Δ Degrees 
of freedom
χ20.95,1 
p-value
χ20.95,m-k 
p-value
Exponential 169 160 3 3.84 
0
7.81 
0
Beta Poisson 8.73 2 5.99 
0.0127
Neither the exponential nor beta-Poisson fits well; beta-Poisson is less bad.

 

Optimized parameters for the beta-Poisson model, from 10000 bootstrap iterations
Parameter MLE estimate Percentiles
0.5% 2.5% 5% 95% 97.5% 99.5%
α 2.65E-01 1.47E-01 1.74E-01 1.87E-01 3.82E-01 4.09E-01 4.70E-01
N50 1.48E+03 5.60E+02 7.33E+02 8.26E+02 2.55E+03 2.78E+03 3.38E+03
Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters.
Parameter scatter plot for beta Poisson model ellipses signify the 0.9, 0.95 and 0.99 confidence of the parameters.
beta Poisson model plot, with confidence bounds around optimized model
beta Poisson model plot, with confidence bounds around optimized model

References