Qualitative and quantitative concentration-response modelling of gene co-expression networks to unlock hepatotoxic mechanisms for next generation chemical safety assessment.

ABSTRACT

Next generation risk assessment of chemicals revolves around the use of mechanistic information without animal experimentation. In this regard, toxicogenomics has proven to be a useful tool to elucidate the underlying mechanisms of adverse effects of xenobiotics. In the present study, two widely used human in vitro hepatocyte culture systems, namely primary human hepatocytes (PHH) and human hepatoma HepaRG cells, were exposed to liver toxicants known to induce liver cholestasis, steatosis or necrosis. Benchmark concentration-response modelling was applied to transcriptomics gene co-expression networks (modules) to derive benchmark concentrations (BMCs) and to gain mechanistic insight into the hepatotoxic effects. BMCs derived by concentration-response modelling of gene co-expression modules recapitulated concentration-response modelling of individual genes. Although PHH and HepaRG cells showed overlap in deregulated genes and modules by the liver toxicants, PHH demonstrated a higher responsiveness, based on the lower BMCs of co-regulated gene modules. Such BMCs can be used as transcriptomics point of departure (tPOD) for assessing module-associated cellular (stress) pathways/processes. This approach identified clear tPODs of around maximum systemic concentration (Cmax) levels for the tested drugs, while for cosmetics ingredients the BMCs were 10-100-fold higher than the estimated plasma concentrations. This approach could serve next generation risk assessment practice to identify early responsive modules at low BMCs, that could be linked to key events in liver adverse outcome pathways. In turn, this can assist in delineating potential hazards of new test chemicals using in vitro systems and used in a risk assessment when BMCs are paired with chemical exposure assessment.

Risk assessment of chemicals has traditionally been focused on animal experiments. In contrast, next generation risk assessment uses biological information obtained from experiments in cell culture models without animals to identify potential hazards. Since the liver is the main target organ of toxicity, many liver cell (hepatocyte) models have been developed and applied for hazard assessment. In this study, two widely used human hepatocyte cell models, PHH and HepaRG, were exposed to liver toxic chemicals. Biological changes in gene expression were measured in a concentration range to identify the concentration at which a biological response was perturbed using concentration response modelling. Genes belonging to the same biological process were joined based on co-expression to derive an average concentration of this process. This animal-free approach could be applied for risk assessment when biological response concentrations were related to the expected human exposure to identify potential hazard of the test chemicals.