High throughput data-based, toxicity pathway-oriented development of a quantitative adverse outcome pathway network linking AHR activation to lung damages

The quantitative adverse outcome pathway (qAOP) is proposed to inform dose-responses at multiple biological levels for the purpose of toxicity prediction. So far, qAOP models concerning human health are scarce. Previously, we proposed 5 key molecular pathways that led aryl hydrogen receptor (AHR) activation to lung damages. The present study assembled an AOP network based on the gene expression signatures of these toxicity pathways, and validated the network using publicly available high throughput data combined with machine learning models. In addition, the AOP network was quantitatively evaluated with omics approaches and bioassays, using 16HBE-CYP1A1 cells exposed to benzo(a)pyrene (BaP), a prototypical AHR activator. Benchmark dose (BMD) analysis of transcriptomics revealed that AHR gene held the lowest BMD value, whereas AHR pathway held the lowest point of departure (PoD) compared to the other 4 pathways. Targeted bioassays were further performed to quantitatively understand the cellular responses, including ROS generation, DNA damage, interleukin-6 production, and extracellular matrix increase marked by collagen expression. Eventually, response-response relationships were plotted using nonlinear model fitting. The present study developed a highly reliable AOP model concerning human health, and validated as well as quantitatively evaluated it, and such a method is likely to be adoptable for risk assessment. [Display omitted] •An AOP network model was constructed concerning AHR-initiated lung damages.•Gene expression signatures of toxicity pathways represented the molecular events.•AOP was validated using high throughput data combined with machine learning models.•AOP quantitative evaluation was proposed using RNA-seq technique and bioassays.