Mouse models of non-alcoholic fatty liver disease (NAFLD) are required to define therapeutic targets, but detailed time-resolved studies to establish a sequence of events are lacking. Here, we fed male C57Bl/6N mice a Western or standard diet over 48 weeks. Multiscale time-resolved characterization was performed using RNA-seq, histopathology, immunohistochemistry, intravital imaging, and blood chemistry; the results were compared to human disease. Acetaminophen toxicity and ammonia metabolism were additionally analyzed as functional readouts. We identified a sequence of eight key events: formation of lipid droplets; inflammatory foci; lipogranulomas; zonal reorganization; cell death and replacement proliferation; ductular reaction; fibrogenesis; and hepatocellular cancer. Functional changes included resistance to acetaminophen and altered nitrogen metabolism. The transcriptomic landscape was characterized by two large clusters of monotonously increasing or decreasing genes, and a smaller number of ‘rest-and-jump genes’ that initially remained unaltered but became differentially expressed only at week 12 or later. Approximately 30% of the genes altered in human NAFLD are also altered in the present mouse model and an increasing overlap with genes altered in human HCC occurred at weeks 30–48. In conclusion, the observed sequence of events recapitulates many features of human disease and offers a basis for the identification of therapeutic targets.

An in vitro/in silico method that determines the risk of human drug induced liver injury in relation to oral doses and blood concentrations of drugs was recently introduced. This method utilizes information on the maximal blood concentration (C_max) for a specific dose of a test compound, which can be estimated using physiologically-based pharmacokinetic modelling, and a cytotoxicity test in cultured human hepatocytes. In the present study, we analyzed if the addition of an assay that measures the inhibition of bile acid export carriers, like BSEP and/or MRP2, to the existing method improves the differentiation of hepatotoxic and non-hepatotoxic compounds. Therefore, an export assay for 5-chloromethylfluorescein diacetate (CMFDA) was established. We tested 36 compounds in a concentration-dependent manner for which the risk of hepatotoxicity for specific oral doses and the capacity to inhibit hepatocyte export carriers are known. Compared to the CTB cytotoxicity test, substantially lower EC₁₀ values were obtained using the CMFDA assay for several known BSEP and/or MRP2 inhibitors. To quantify if the addition of the CMFDA assay to our test system improves the overall separation of hepatotoxic from non-hepatotoxic compounds, the toxicity separation index (TSI) was calculated. We obtained a better TSI using the lower alert concentration from either the CMFDA or the CTB test (TSI: 0.886) compared to considering the CTB test alone (TSI: 0.775). In conclusion, the data show that integration of the CMFDA assay with an in vitro test battery improves the differentiation of hepatotoxic and non-hepatotoxic compounds in a set of compounds that includes bile acid export carrier inhibitors.

High throughput RNA sequencing experiments are widely conducted and analyzed to identify differentially expressed genes (DEGs). The statistical models calculated for this task are often not clear to practitioners, and analyses may not be optimally tailored to the research hypothesis. Often, interaction effects (IEs) are the mathematical equivalent of the biological research question but are not considered for different reasons. We fill this gap by explaining and presenting the potential benefit of IEs in the search for DEGs using RNA-Seq data of mice that receive different diets for different time periods. Using an IE model leads to a smaller, but likely more biologically informative set of DEGs compared to a common approach that avoids the calculation of IEs.

We extend the scope of application for MCP-Mod (Multiple Comparison Procedure and Modeling) to in vitro gene expression data and assess its characteristics regarding model selection for concentration gene expression curves. Precisely, we apply MCP-Mod on single genes of a high-dimensional gene expression data set, where human embryonic stem cells were exposed to eight concentration levels of the compound valproic acid (VPA). As candidate models we consider the sigmoid $$E_{\max }$$ (four-parameter log-logistic), linear, quadratic, $$E_{\max }$$, exponential, and beta model. Through simulations we investigate the impact of omitting one or more models from the candidate model set to uncover possibly superfluous models and to evaluate the precision and recall rates of selected models. Each model is selected according to Akaike information criterion (AIC) for a considerable number of genes. For less noisy cases the popular sigmoid $$E_{\max }$$ model is frequently selected. For more noisy data, often simpler models like the linear model are selected, but mostly without relevant performance advantage compared to the second best model. Also, the commonly used standard $$E_{\max }$$ model has an unexpected low performance.

The alkaline comet assay is frequently used as in vivo follow-up test within different regulatory environments to characterize the DNA-damaging potential of different test items. The corresponding OECD Test guideline 489 highlights the importance of statistical analyses and historical control data (HCD) but does not provide detailed procedures. Therefore, the working group "Statistics" of the German-speaking Society for Environmental Mutation Research (GUM) collected HCD from five laboratories and >200 comet assay studies and performed several statistical analyses. Key results included that (I) observed large inter-laboratory effects argue against the use of absolute quality thresholds, (II) > 50% zero values on a slide are considered problematic, due to their influence on slide or animal summary statistics, (III) the type of summarizing measure for single-cell data (e.g., median, arithmetic and geometric mean) may lead to extreme differences in resulting animal tail intensities and study outcome in the HCD. These summarizing values increase the reliability of analysis results by better meeting statistical model assumptions, but at the cost of information loss. Furthermore, the relation between negative and positive control groups in the data set was always satisfactorily (or sufficiently) based on ratio, difference and quantile analyses.

Background and Aims: Metabolic dysfunction–associated steatohepatitis (MASH) is associated with c-Jun N-terminal kinase (JNK) activation across various cell types, but its hepatocyte-specific function in steatotic liver disease remains unclear. Our study investigates the role of JNK1/JNK2 during MASH progression and dissects its hepatocyte-specific function. Approach and Results: We showed that UK Biobank patients with a predicted loss-of-function variant of JNK1 presented an increased prevalence of metabolic dysfunction–associated steatotic liver disease and liver damage. Analysis of a pathology cohort of patients with steatotic liver disease revealed increased oxidative stress response and apoptosis. After subjecting mice deficient for Jnk1 and Jnk2 in hepatocytes ( Jnk1/2 ^Δhepa ) to 2 different MASH models, we observed enhanced liver injury, fibrosis, and oxidative stress. RNA sequencing revealed highly upregulated pathways in Jnk1/2 ^Δhepa livers, including inflammatory signals and apoptotic pathways. Additional blocking of Caspase 8 signaling improved high-fat diet–induced liver damage, fibrogenesis, and oxidative stress. Ultimately, a therapeutic approach using lipid nanoparticles containing small interfering RNA targeting Caspase 8 during MASH progression attenuated liver injury and cell death in mice. Conclusions: Our findings define a protective role of JNK1/JNK2 in hepatocytes during the oxidative stress response driving the progression of MASH. This process is mainly mediated by Caspase 8-dependent apoptosis, thereby discovering that Caspase 8 is a downstream target of JNK1/2. Caspase 8-directed therapy in hepatocytes might be a promising treatment for patients with an increased oxidative stress response and MASH.