Species Sensitivity Distribution Estimation from Uncertain (QSAR-based) Effects Data

Abstract

In environmental risk assessment, Species Sensitivity Distributions (SSDs) can be applied to estimate a PNEC (Predicted No-Effect Concentration) for a chemical substance, when sufficient data on species toxicities are available. The European Chemicals Agency (ECHA) recommendation is 10 biological species. The question addressed in this paper, is whether QSAR-predicted toxicities can be included in SSD-based PNEC estimates, and whether any modifications need to be made to account for the uncertainty in the QSAR-model estimates. This problem is addressed from a probabilistic modelling point of view. From classical analysis of variation (ANOVA), we review how the error-in-data SSD problem is similar to separation into between-group and within-group variance. ECHA guidance suggests averaging similar endpoint data for a species, which is consistent with group means, as in ANOVA. This exercise reveals that error-in-data reduces the estimation of the between species variation, i.e. the SSD variance, rather than enlarging it. A Bayesian analysis permits the assessment of the uncertainty of the SSD mean and variance parameters for given values of mean species toxicity error. This requires a hierarchical model. Prototyping this model for an artificial five-species data set seems to suggest that the influence of data error is relatively minor. Moreover, when neglecting this data error, a slightly conservative estimate of the SSD results. Hence, we suggest including (model-predicted) data as model point estimates and handling the SSD as usual. The Bayesian simulation of the error-in-data SSD leads to predictive distributions, being an average of posterior spaghetti plot densities or cumulative distributions. We derive new predictive extrapolation constants with several improvements over previous median uncertainty log₁₀HC5 estimates, in that they are easily calculable from spreadsheet Student-t functions and based on a more realistic uniform prior for the SSD standard deviation. Other advantages are that they are single-number extrapolation constants and they are more sensitive to small sample size.

Keywords

Bayesian hierarchical model extrapolation constants predictive uncertainty probabilistic risk assessment QSAR and read-across prediction species sensitivity distribution

References

ECHA (2008). Guidance on information requirements and chemical safety assessment. Chapter R.10: Characterisation of dose [concentration]–response for environment, 65 pp. Helsinki, Finland: European Chemicals Agency.

Aldenberg

, & Slob

(1993). Confidence limits for hazardous concentrations based on logistically distributed NOEC toxicity data. Ecotoxicology & Environmental Safety 25, 48–63.

Aldenberg

, & Jaworska

J.S.

(2000). Uncertainty of the hazardous concentration and fraction affected for normal species sensitivity distributions. Ecotoxicology & Environmental Safety 46, 1–18.

Wagner

, & Løkke

(1991). Estimation of ecotoxicological protection levels from NOEC toxicity data. Water Research 25, 1237–1242.

Van Straalen

N.M.

, & Denneman

C.A.J.

(1989). Ecotoxicological evaluation of soil quality criteria. Ecotoxicology & Environmental Safety 18, 241–251.

Gelman

, Carlin

J.B.

, Stern

H.S.

, & Rubin

D.B.

(2004). Bayesian Data Analysis, 2nd edn, 668 pp. Boca Raton, FL, USA: Chapman & Hall/CRC.

Box

G.E.P.

, & Tiao

G.C.

(1973/1992). Bayesian Inference in Statistical Analysis, 588 pp. New York, NY, USA: Wiley.

Aldenberg

, Jaworska

J.S.

, & Traas

T.P.

(2002). Normal species sensitivity distributions and probabilistic ecological risk assessment. In Species Sensitivity Distributions in Ecotoxicology (ed. Posthuma

, Suter

G.W.

II , & Traas

T.P.

), pp. 49–102. Boca Raton, FL, USA: CRC/Lewis Publishers.

Aldenberg

, & Luttik

(2002). Extrapolation factors for tiny toxicity data sets from Species Sensitivity Distributions with known standard deviation. In Species Sensitivity Distributions in Ecotoxicology (ed. Posthuma

, Suter

G.W.

II , & Traas

T.P.

), pp. 103–118. Boca Raton, FL, USA: CRC/Lewis Publishers.

10.

Rorije

, Aldenberg

, & Peijnenburg

W.J.G.M.

(2013). Read-across estimates of aquatic toxicity for selected fragrances. ATLA 41, 77–90.

11.

Draper

N.R.

, & Smith

(1998). Applied Regression Analysis, 3rd edn, 706 pp. New York, NY, USA: Wiley.

12.

Helsel

D.R.

, & Hirsch

R.M.

(2002). Statistical Methods in Water Resources, 524 pp. Reston, VA, USA: US Geological Survey. Available at: http://pubs.usgs.gov/twri/twri4a3/#pdf (Accessed 02.02.12).

13.

Dalgaard

(2008). Introductory Statistics with R, 2nd edn, 363 pp. New York, NY, USA: Springer.