Relating suborganismal processes to ecotoxicological and population level endpoints using a bioenergetic model

Ecological effects of environmental stressors are commonly evaluated using organismal or suborganismal data, such as standardized toxicity tests that characterize responses of individuals (e.g., mortality and reproduction) and a rapidly growing body of “omics” data. A key challenge for environmental risk assessment is relating such information to population dynamics. One approach uses dynamic energy budget (DEB) models that relate growth and reproduction of individuals to underlying flows of energy and elemental matter. We hypothesize that suborganismal information identifies DEB parameters that are most likely impacted by a particular stressor and that the DEB model can then project suborganismal effects on life history and population endpoints. We formulate and parameterize a model of growth and reproduction for the water flea Daphnia magna. Our model resembles previous generic bioenergetic models, but has explicit representation of discrete molts, an important feature of Daphnia life history. We test its ability to predict six endpoints commonly used in chronic toxicity studies in specified food environments. With just one adjustable parameter, the model successfully predicts growth and reproduction of individuals from a wide array of experiments performed in multiple laboratories using different clones of D. magna raised on different food sources. Fecundity is the most sensitive endpoint, and there is broad correlation between the sensitivities of fecundity and long‐run growth rate, as is desirable for the default metric used in chronic toxicity tests. Under some assumptions, we can combine our DEB model with the Euler‐Lotka equation to estimate long‐run population growth rates at different food levels. A review of Daphnia gene‐expression experiments on the effects of contaminant exposure reveals several connections to model parameters, in particular a general trend of increased transcript expression of genes involved in energy assimilation and utilization at concentrations affecting growth and reproduction. The sensitivity of fecundity to many model parameters was consistent with frequent generalized observations of decreased expression of genes involved in reproductive physiology, but interpretation of these observations requires further mechanistic modeling. We thus propose an approach based on generic DEB models incorporating few essential species‐specific features for rapid extrapolation of ecotoxicogenomic assays for Daphnia‐based population risk assessment.