Search

5-Fluorouracil (5-FU) is a widely used chemotherapeutical that induces acute toxicity in the small and large intestine of patients. Symptoms can be severe and lead to the interruption of cancer treatments. However, there is limited understanding of the molecular mechanisms underlying 5-FU-induced intestinal toxicity. In this study, well-established 3D organoid models of human colon and small intestine (SI) were used to characterize 5-FU transcriptomic and metabolomic responses. Clinically relevant 5-FU concentrations for in vitro testing in organoids were established using physiologically based pharmacokinetic simulation of dosing regimens recommended for cancer patients, resulting in exposures to 10, 100 and 1000 µM. After treatment, different measurements were performed: cell viability and apoptosis; image analysis of cell morphological changes; RNA sequencing; and metabolome analysis of supernatant from organoids cultures. Based on analysis of the differentially expressed genes, the most prominent molecular pathways affected by 5-FU included cell cycle, p53 signalling, mitochondrial ATP synthesis and apoptosis. Short time-series expression miner demonstrated tissue-specific mechanisms affected by 5-FU, namely biosynthesis and transport of small molecules, and mRNA translation for colon; cell signalling mediated by Rho GTPases and fork-head box transcription factors for SI. Metabolomic analysis showed that in addition to the effects on TCA cycle and oxidative stress in both organoids, tissue-specific metabolic alterations were also induced by 5-FU. Multi-omics integration identified transcription factor E2F1, a regulator of cell cycle and apoptosis, as the best key node across all samples. These results provide new insights into 5-FU toxicity mechanisms and underline the relevance of human organoid models in the safety assessment in drug development., (© 2021. The Author(s).)

Unravelling gene regulatory networks (GRNs) influenced by chemicals is a major challenge in systems toxicology. Because toxicant-induced GRNs evolve over time and dose, the analysis of global gene expression data measured at multiple time points and doses will provide insight in the adverse effects of compounds. Therefore, there is a need for mathematical methods for GRN identification from time-over-dose-dependent data. One of the current approaches for GRN inference is Time Series Network Identification (TSNI). TSNI is based on ordinary differential equations (ODE), describing the time evolution of the expression of each gene, which is assumed to be dependent on the expression of other genes and an external perturbation (i.e. chemical exposure). Here, we present Dose-Time Network Identification (DTNI), a method extending TSNI by including ODE describing how the expression of each gene evolves with dose, which is supposed to depend on the expression of other genes and the exposure time. We also adapted TSNI in order to enable inclusion of time-over-dose-dependent data from multiple compounds. Here, we show that DTNI outperforms TSNI in inferring a toxicant-induced GRN. Moreover, we show that DTNI is a suitable method to infer a GRN dose- and time-dependently induced by a group of compounds influencing a common biological process. Applying DTNI on experimental data from TG-GATEs, we demonstrate that DTNI provides in-depth information on the mode of action of compounds, in particular key events and potential molecular initiating events. Furthermore, DTNI also discloses several unknown interactions which have to be verified experimentally.