Abstract
The generation of omics data sets has become an important approach in modern pharmacological and toxicological research as it can provide mechanistic and quantitative information on a large scale. Analyses of these data frequently revealed a non-linear dose-response relationship underscoring the importance of the modeling process to infer biological exposure limits. A number of tools have been developed for dose-response modeling and various thresholds have been defined as a quantitative representation of the effect of a substance, such as effective concentrations or benchmark doses (BMD). Here we present DoseRider an easy-to-use web application and a companion R package for linear and non-linear dose-response modeling and assessment of BMD at the level of biological pathways or signatures using generalized mixed effect models. This approach allows to analyze custom or provided multi-omics data such as RNA sequencing or metabolomics data and its annotation of a collection of pathways and gene sets from various species. Moreover, we introduce the concept of the trend change doses (TCDs) as a numerical descriptor of effects derived from complex dose-response curves. The usability of DoseRider was demonstrated by analyses of RNA sequencing data of bisphenol AF (BPAF) treatment of a human breast cancer cell line (MCF-7) at 8 different concentrations using gene sets for chemical and genetic perturbations (MSigDB). The BMD for BPAF and a set of genes upregulated by estrogen in breast cancer was 0.2 µM (95 %-CI 0.1-0.5 µM) and the lowest TCD (TCD1) was 0.003 µM (95 %-CI 0.0006-0.01 µM). The comprehensive presentation of the results underlines the suitability of the system for pharmacogenomics, toxicogenomics, and applications beyond.