Quantifying fungal growth underpins our ability to effectively treat severe fungal infections. Current methods quantify fungal growth rates from time-course morphology-specific data, such as hyphal length data. However, automated large-scale collection of such data lies beyond the scope of most clinical microbiology laboratories. In this paper, we propose a mathematical model of fungal growth to estimate morphology-specific growth rates from easy-to-collect, but indirect, optical density (OD600) data of Aspergillus fumigatus growth (filamentous fungus). Our method accounts for OD600 being an indirect measure by explicitly including the relationship between the indirect OD600 measurements and the calibrating true fungal growth in the model. Therefore, the method does not require de novo generation of calibration data. Our model outperformed reference models at fitting to and predicting OD600 growth curves and overcame observed discrepancies between morphology-specific rates inferred from OD600 versus directly measured data in reference models that did not include calibration.

MOTIVATION: Approximate Bayesian computation (ABC) is an important framework within which to infer the structure and parameters of a systems biology model. It is especially suitable for biological systems with stochastic and nonlinear dynamics, for which the likelihood functions are intractable. However, the associated computational cost often limits ABC to models that are relatively quick to simulate in practice. RESULTS: We here present a Julia package, GpABC, that implements parameter inference and model selection for deterministic or stochastic models using (i) standard rejection ABC or sequential Monte Carlo ABC or (ii) ABC with Gaussian process emulation. The latter significantly reduces the computational cost. AVAILABILITY AND IMPLEMENTATION: https://github.com/tanhevg/GpABC.jl.