In vivo alkaline comet assay: Statistical considerations on historical negative and positive control data.

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.

Inhibition of the renal apical sodium dependent bile acid transporter prevents cholemic nephropathy in mice with obstructive cholestasis.

BACKGROUND & AIMS: Cholemic nephropathy (CN) is a severe complication of cholestatic liver diseases for which there is no specific treatment. We revisited its pathophysiology with the aim of identifying novel therapeutic strategies. METHODS: Cholestasis was induced by bile duct ligation (BDL) in mice. Bile flux in kidneys and livers was visualized by intravital imaging, supported by MALDI mass spectrometry imaging and liquid chromatography-tandem mass spectrometry. The effect of AS0369, a systemically bioavailable apical sodium-dependent bile acid transporter (ASBT) inhibitor, was evaluated by intravital imaging, RNA-sequencing, histological, blood, and urine analyses. Translational relevance was assessed in kidney biopsies from patients with CN, mice with a humanized bile acid (BA) spectrum, and via analysis of serum BAs and KIM-1 (kidney injury molecule 1) in patients with liver disease and hyperbilirubinemia. RESULTS: Proximal tubular epithelial cells (TECs) reabsorbed and enriched BAs, leading to oxidative stress and death of proximal TECs, casts in distal tubules and collecting ducts, peritubular capillary leakiness, and glomerular cysts. Renal ASBT inhibition by AS0369 blocked BA uptake into TECs and prevented kidney injury up to 6 weeks after BDL. Similar results were obtained in mice with humanized BA composition. In patients with advanced liver disease, serum BAs were the main determinant of KIM-1 levels. ASBT expression in TECs was preserved in biopsies from patients with CN, further highlighting the translational potential of targeting ASBT to treat CN. CONCLUSIONS: BA enrichment in proximal TECs followed by oxidative stress and cell death is a key early event in CN. Inhibiting renal ASBT and consequently BA enrichment in TECs prevents CN and systemically decreases BA concentrations. IMPACT AND IMPLICATIONS: Cholemic nephropathy (CN) is a severe complication of cholestasis and an unmet clinical need. We demonstrate that CN is triggered by the renal accumulation of bile acids (BAs) that are considerably increased in the systemic blood. Specifically, the proximal tubular epithelial cells of the kidney take up BAs via the apical sodium-dependent bile acid transporter (ASBT). We developed a therapeutic compound that blocks ASBT in the kidneys, prevents BA overload in tubular epithelial cells, and almost completely abolished all disease hallmarks in a CN mouse model. Renal ASBT inhibition represents a potential therapeutic strategy for patients with CN.

Benefit of using interaction effects for the analysis of high-dimensional time-response or dose-response data for two-group comparisons.

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.

Guidance for statistical design and analysis of toxicological dose-response experiments, based on a comprehensive literature review.

The analysis of dose-response, concentration-response, and time-response relationships is a central component of toxicological research. A major decision with respect to the statistical analysis is whether to consider only the actually measured concentrations or to assume an underlying (parametric) model that allows extrapolation. Recent research suggests the application of modelling approaches for various types of toxicological assays. However, there is a discrepancy between the state of the art in statistical methodological research and published analyses in the toxicological literature. The extent of this gap is quantified in this work using an extensive literature review that considered all dose-response analyses published in three major toxicological journals in 2021. The aspects of the review include biological considerations (type of assay and of exposure), statistical design considerations (number of measured conditions, design, and sample sizes), and statistical analysis considerations (display, analysis goal, statistical testing or modelling method, and alert concentration). Based on the results of this review and the critical assessment of three selected issues in the context of statistical research, concrete guidance for planning, execution, and analysis of dose-response studies from a statistical viewpoint is proposed.

Long-term hypercaloric diet exacerbates metabolic liver disease in PNPLA3 I148M animals.

BACKGROUND & AIMS: Nonalcoholic fatty liver disease (NAFLD) is a major health burden associated with the metabolic syndrome leading to liver fibrosis, cirrhosis and ultimately liver cancer. In humans, the PNPLA3 I148M polymorphism of the phospholipase patatin-like phospholipid domain containing protein 3 (PNPLA3) has a well-documented impact on metabolic liver disease. In this study, we used a mouse model mimicking the human PNPLA3 I148M polymorphism in a long-term high fat diet (HFD) experiment to better define its role for NAFLD progression. METHODS: Male mice bearing wild-type Pnpla3 (Pnpla3WT ), or the human polymorphism PNPLA3 I148M (Pnpla3148M/M ) were subjected to HFD feeding for 24 and 52 weeks. Further analysis concerning basic phenotype, inflammation, proliferation and cell death, fibrosis and microbiota were performed in each time point. RESULTS: After 52 weeks HFD Pnpla3148M/M animals had more liver fibrosis, enhanced numbers of inflammatory cells as well as increased Kupffer cell activity. Increased hepatocyte cell turnover and ductular proliferation were evident in HFD Pnpla3148M/M livers. Microbiome diversity was decreased after HFD feeding, changes were influenced by HFD feeding (36%) and the PNPLA3 I148M genotype (12%). Pnpla3148M/M mice had more faecal bile acids. RNA-sequencing of liver tissue defined an HFD-associated signature, and a Pnpla3148M/M specific pattern, which suggests Kupffer cell and monocytes-derived macrophages as significant drivers of liver disease progression in Pnpla3148M/M animals. CONCLUSION: With long-term HFD feeding, mice with the PNPLA3 I148M genotype show exacerbated NAFLD. This finding is linked to PNPLA3 I148M-specific changes in microbiota composition and liver gene expression showing a stronger inflammatory response leading to enhanced liver fibrosis progression.

Effect of Whole-Body Vibration Training on Selected Intrinsic Risk Factors in Women Aged 60+ at Fall Risk: A Randomized Controlled Trial.

The aim of the study was to determine whether Whole Body Vibration Training (WBVT) affects intrinsic risk factors for falls in women aged 60+ at fall risk. DESIGN: Randomized controlled clinical trial. Blinding was applied to the persons in charge of evaluating the intervention's clinical results and statistical analysis. METHODS: Forty-two women over 60 years old were randomly assigned to an experimental group (EG-12-week WBVT; n = 22) and a control group (CG-no additional physical activities; n = 20). Fear of falling was measured by the FES-I questionnaire, gait and dynamic balance using the Time-Up and Go test (TUG), aerobic endurance with the 6-Minute Walk Test (6MWT), and the functional strength of the lower body muscles with the 30-s Chair Stand Test (30SCST) at baseline and post-intervention. Additionally assayed were participants' blood concentrations of interleukin-6 (IL-6). RESULTS: The 12-week WBVT improves gait and balance (TUG, p = 0.009), exercise tolerance (6MWT, p = 0.001), and functional strength (30SCST; p = 0.027) but does not reduce the intensity of fear of falling (FES-I, p = 0.655) and the IL-6 serum concentration (p = 0.377). CONCLUSIONS: WBVT affects selected fall risk factors in women aged 60+ at fall risk.

In vitro/in silico prediction of drug induced steatosis in relation to oral doses and blood concentrations by the Nile Red assay.

The accumulation of lipid droplets in hepatocytes is a key feature of drug-induced liver injury (DILI) and can be induced by a subset of hepatotoxic compounds. In the present study, we optimized and evaluated an in vitro technique based on the fluorescent dye Nile Red, further named Nile Red assay to quantify lipid droplets induced by the exposure to chemicals. The Nile Red assay and a cytotoxicity test (CTB assay) were then performed on cells exposed concentration-dependently to 60 different compounds. Of these, 31 were known to induce hepatotoxicity in humans, and 13 were reported to also cause steatosis. In order to compare in vivo relevant blood concentrations, pharmacokinetic models were established for all compounds to simulate the maximal blood concentrations (Cmax) at therapeutic doses. The results showed that several hepatotoxic compounds induced an increase in lipid droplets at sub-cytotoxic concentrations. To compare how well (1) the cytotoxicity test alone, (2) the Nile Red assay alone, and (3) the combination of the cytotoxicity test and the Nile Red assay (based on the lower EC10 of both assays) allow the differentiation between hepatotoxic and non-hepatotoxic compounds, a previously established performance metric, the Toxicity Separation Index (TSI) was calculated. In addition, the Toxicity Estimation Index (TEI) was calculated to determine how well blood concentrations that cause an increased DILI risk can be estimated for hepatotoxic compounds. Our findings indicate that the combination of both assays improved the TSI and TEI compared to each assay alone. In conclusion, the study demonstrates that inclusion of the Nile Red assay into in vitro test batteries may improve the prediction of DILI compounds.

Influence of bile acids on the cytotoxicity of chemicals in cultivated human hepatocytes.

Bile acids (BA) are known to influence the susceptibility of hepatocytes to chemicals. We investigated the cytotoxicity of 18 compounds with known hepatotoxicity status and pharmacokinetics in cultivated primary human hepatocytes with and without the addition of a BA mix to the cell culture medium. This BA mix consisted of physiological ratios of the most abundant human BA at a cholestatic sum concentration of 0.5 mM, which corresponds to 50% of the EC10 (cytotoxicity) of the mix. The BA mix decreased the EC10 of 7 compounds by a factor greater than 1.5, but also increased the EC10 of 5 compounds. The compounds with increased susceptibility include the known hepatotoxicants and BSEP/MRP2 inhibitors rifampicin, ketoconazole, atorvastatin, and cyclosporin A. However, the cytotoxicity of some non-hepatotoxic compounds was also enhanced, among them glucose, which is not known to be an inhibitor of canalicular bile acid export. A recently established technique to quantify how well hepatotoxic and non-hepatotoxic compounds are separated by an in vitro test indicated that the addition of the BA mix did not improve separation. In conclusion, the addition of BA to cultivated hepatocytes leads to a complex situation with increased and decreased susceptibilities depending on the specific compound.

Model selection characteristics when using MCP-Mod for dose-response gene expression data.

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 Emax$E_{\max }$ (four-parameter log-logistic), linear, quadratic, Emax$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 Emax$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 Emax$E_{\max }$ model has an unexpected low performance.

The hepatocyte export carrier inhibition assay improves the separation of hepatotoxic from non-hepatotoxic compounds.

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 (Cmax) 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 EC10 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.

Spatio-Temporal Multiscale Analysis of Western Diet-Fed Mice Reveals a Translationally Relevant Sequence of Events during NAFLD Progression.

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.