The impact of glucagon to support postabsorptive glucose flux and glycemia in healthy rats and its attenuation in male Zucker diabetic fatty rats.

Hyperglucagonemia is a hallmark of type 2 diabetes (T2DM), yet the role of elevated plasma glucagon (P-GCG) to promote excessive postabsorptive glucose production and contribute to hyperglycemia in patients with this disease remains debatable. We investigated the acute action of P-GCG to safeguard/support postabsorptive endogenous glucose production (EGP) and euglycemia in healthy Zucker control lean (ZCL) rats. Using male Zucker diabetic fatty (ZDF) rats that exhibit the typical metabolic disorders of human T2DM, such as excessive EGP, hyperglycemia, hyperinsulinemia, and hyperglucagonemia, we examined the ability of hyperglucagonemia to promote greater rates of postabsorptive EGP and hyperglycemia. Euglycemic or hyperglycemic basal insulin (INS-BC) and glucagon (GCG-BC) clamps were performed in the absence or during an acute setting of glucagon deficiency (GCG-DF, ∼10% of basal), either alone or in combination with insulin deficiency (INS-DF, ∼10% of basal). Glucose appearance, disappearance, and cycling rates were measured using [2-3H] and [3-3H]-glucose. In ZCL rats, GCG-DF reduced the levels of hepatic cyclic AMP, EGP, and plasma glucose (PG) by 50%, 32%, and 50%, respectively. EGP fell in the presence GCG-DF and INS-BC, but under GCG-DF and INS-DF, EGP and PG increased two- and threefold, respectively. GCG-DF revealed the hyperglucagonemia present in ZDF rats lacked the ability to regulate hepatic intracellular cyclic AMP levels and glucose flux, since EGP and PG levels fell by only 10%. We conclude that the liver in T2DM suffers from resistance to all three major regulatory factors, glucagon, insulin, and glucose, thus leading to a loss of metabolic flexibility.NEW & NOTEWORTHY In postabsorptive state, basal plasma insulin (P-INS) and plasma glucose (PG) act dominantly to increase hepatic glucose cycling and reduce endogenous glucose production (EGP) and PG in healthy rats, which is only counteracted by the acute action of basal plasma glucagon (P-GCG) to support EGP and euglycemia. Hyperglucagonemia, a hallmark of type 2 diabetes (T2DM) present in Zucker diabetic fatty (ZDF) rats, is not the primary mediator of hyperglycemia and high EGP as commonly thought; instead, the liver is resistant to glucagon as well as insulin and glucose.

Probing Liver Injuries Induced by Thioacetamide in Human In Vitro Pooled Hepatocyte Experiments.

Animal studies are typically utilized to understand the complex mechanisms associated with toxicant-induced hepatotoxicity. Among the alternative approaches to animal studies, in vitro pooled human hepatocytes have the potential to capture population variability. Here, we examined the effect of the hepatotoxicant thioacetamide on pooled human hepatocytes, divided into five lots, obtained from forty diverse donors. For 24 h, pooled human hepatocytes were exposed to vehicle, 1.33 mM (low dose), and 12 mM (high dose) thioacetamide, followed by RNA-seq analysis. We assessed gene expression variability using heat maps, correlation plots, and statistical variance. We used KEGG pathways and co-expression modules to identify underlying physiological processes/pathways. The co-expression module analysis showed that the majority of the lots exhibited activation for the bile duct proliferation module. Despite lot-to-lot variability, we identified a set of common differentially expressed genes across the lots with similarities in their response to amino acid, lipid, and carbohydrate metabolism. We also examined efflux transporters and found larger lot-to-lot variability in their expression patterns, indicating a potential for alteration in toxicant bioavailability within the cells, which could in turn affect the gene expression patterns between the lots. Overall, our analysis highlights the challenges in using pooled hepatocytes to understand mechanisms of toxicity.

Assessing Kidney Injury Induced by Mercuric Chloride in Guinea Pigs with In Vivo and In Vitro Experiments.

Acute kidney injury, which is associated with high levels of morbidity and mortality, affects a significant number of individuals, and can be triggered by multiple factors, such as medications, exposure to toxic chemicals or other substances, disease, and trauma. Because the kidney is a critical organ, understanding and identifying early cellular or gene-level changes can provide a foundation for designing medical interventions. In our earlier work, we identified gene modules anchored to histopathology phenotypes associated with toxicant-induced liver and kidney injuries. Here, using in vivo and in vitro experiments, we assessed and validated these kidney injury-associated modules by analyzing gene expression data from the kidneys of male Hartley guinea pigs exposed to mercuric chloride. Using plasma creatinine levels and cell-viability assays as measures of the extent of renal dysfunction under in vivo and in vitro conditions, we performed an initial range-finding study to identify the appropriate doses and exposure times associated with mild and severe kidney injuries. We then monitored changes in kidney gene expression at the selected doses and time points post-toxicant exposure to characterize the mechanisms of kidney injury. Our injury module-based analysis revealed a dose-dependent activation of several phenotypic cellular processes associated with dilatation, necrosis, and fibrogenesis that were common across the experimental platforms and indicative of processes that initiate kidney damage. Furthermore, a comparison of activated injury modules between guinea pigs and rats indicated a strong correlation between the modules, highlighting their potential for cross-species translational studies.

Genomics and metabolomics of early-stage thioacetamide-induced liver injury: An interspecies study between guinea pig and rat.

To study the complex processes involved in liver injuries, researchers rely on animal investigations, using chemically or surgically induced liver injuries, to extrapolate findings and infer human health risks. However, this presents obvious challenges in performing a detailed comparison and validation between the highly controlled animal models and development of liver injuries in humans. Furthermore, it is not clear whether there are species-dependent and -independent molecular initiating events or processes that cause liver injury before they eventually lead to end-stage liver disease. Here, we present a side-by-side study of rats and guinea pigs using thioacetamide to examine the similarities between early molecular initiating events during an acute-phase liver injury. We exposed Sprague Dawley rats and Hartley guinea pigs to a single dose of 25 or 100 mg/kg thioacetamide and collected blood plasma for metabolomic analysis and liver tissue for RNA-sequencing. The subsequent toxicogenomic analysis identified consistent liver injury trends in both genomic and metabolomic data within 24 and 33 h after thioacetamide exposure in rats and guinea pigs, respectively. In particular, we found species similarities in the key injury phenotypes of inflammation and fibrogenesis in our gene module analysis for liver injury phenotypes. We identified expression of several common genes (e.g., SPP1, TNSF18, SERPINE1, CLDN4, TIMP1, CD44, and LGALS3), activation of injury-specific KEGG pathways, and alteration of plasma metabolites involved in amino acid and bile acid metabolism as some of the key molecular processes that changed early upon thioacetamide exposure and could play a major role in the initiation of acute liver injury.

Glucagon-like peptide-1 receptor agonist inhibits aeroallergen-induced activation of ILC2 and neutrophilic airway inflammation in obese mice.

BACKGROUND: Obesity is a risk factor for the development of asthma. However, pharmacologic therapeutic strategies that specifically target obese asthmatics have not been identified. We hypothesize that glucagon-like peptide-1 receptor agonist (GLP-1RA) treatment inhibits aeroallergen-induced early innate airway inflammation in a mouse model of asthma in the setting of obesity. METHODS: SWR (lean) and TALLYHO (obese) mice were challenged intranasally with Alternaria alternata extract (Alt-Ext) or PBS for 4 consecutive days concurrent with GLP-1RA or vehicle treatment. RESULTS: TALLYHO mice had greater Alt-Ext-induced airway neutrophilia and lung protein expression of IL-5, IL-13, CCL11, CXCL1, and CXCL5, in addition to ICAM-1 expression on lung epithelial cells compared with SWR mice, and all endpoints were reduced by GLP-1RA treatment. Alt-Ext significantly increased BALF IL-33 in both TALLYHO and SWR mice compared to PBS challenge, but there was no difference in the BALF IL-33 levels between these two strains. However, TALLYHO, but not SWR, mice had significantly higher airway TSLP in BALF following Alt-Ext challenge compared to PBS, and BALF TSLP was significantly greater in TALLYHO mice compared to SWR mice following airway Alt-Ext challenge. GLP-1RA treatment significantly decreased the Alt-Ext-induced TSLP and IL-33 release in TALLYHO mice. While TSLP or ST2 inhibition with a neutralizing antibody decreased airway eosinophils, they did not reduce airway neutrophils in TALLYHO mice. CONCLUSIONS: These results suggest that GLP-1RA treatment may be a novel pharmacologic therapeutic strategy for obese persons with asthma by inhibiting aeroallergen-induced neutrophilia, a feature not seen with either TSLP or ST2 inhibition.

Concordance between Thioacetamide-Induced Liver Injury in Rat and Human In Vitro Gene Expression Data.

The immense resources required and the ethical concerns for animal-based toxicological studies have driven the development of in vitro and in silico approaches. Recently, we validated our approach in which the expression of a set of genes is uniquely associated with an organ-injury phenotype (injury module), by using thioacetamide, a known liver toxicant. Here, we sought to explore whether RNA-seq data obtained from human cells (in vitro) treated with thioacetamide-S-oxide (a toxic intermediate metabolite) would correlate across species with the injury responses found in rat cells (in vitro) after exposure to this metabolite as well as in rats exposed to thioacetamide (in vivo). We treated two human cell types with thioacetamide-S-oxide (primary hepatocytes with 0 (vehicle), 0.125 (low dose), or 0.25 (high dose) mM, and renal tubular epithelial cells with 0 (vehicle), 0.25 (low dose), or 1.00 (high dose) mM) and collected RNA-seq data 9 or 24 h after treatment. We found that the liver-injury modules significantly altered in human hepatocytes 24 h after high-dose treatment involved cellular infiltration and bile duct proliferation, which are linked to fibrosis. For high-dose treatments, our modular approach predicted the rat in vivo and in vitro results from human in vitro RNA-seq data with Pearson correlation coefficients of 0.60 and 0.63, respectively, which was not observed for individual genes or KEGG pathways.

Toxicant-Induced Metabolic Alterations in Lipid and Amino Acid Pathways Are Predictive of Acute Liver Toxicity in Rats.

Liver disease and disorders associated with aberrant hepatocyte metabolism can be initiated via drug and environmental toxicant exposures. In this study, we tested the hypothesis that gene and metabolic profiling can reveal commonalities in liver response to different toxicants and provide the capability to identify early signatures of acute liver toxicity. We used Sprague Dawley rats and three classical hepatotoxicants: acetaminophen (2 g/kg), bromobenzene (0.4 g/kg), and carbon tetrachloride (0.3 g/kg), to identify early perturbations in liver metabolism after a single acute exposure dose. We measured changes in liver genes and plasma metabolites at two time points (5 and 10 h) and used genome-scale metabolic models to identify commonalities in liver responses across the three toxicants. We found strong correlations for gene and metabolic profiles between the toxicants, indicative of similarities in the liver response to toxicity. We identified several injury-specific pathways in lipid and amino acid metabolism that changed similarly across the three toxicants. Our findings suggest that several plasma metabolites in lipid and amino acid metabolism are strongly associated with the progression of liver toxicity, and as such, could be targeted and clinically assessed for their potential as early predictors of acute liver toxicity.

Mechanistic identification of biofluid metabolite changes as markers of acetaminophen-induced liver toxicity in rats.

Acetaminophen (APAP) is the most commonly used analgesic and antipyretic drug in the world. Yet, it poses a major risk of liver injury when taken in excess of the therapeutic dose. Current clinical markers do not detect the early onset of liver injury associated with excess APAP-information that is vital to reverse injury progression through available therapeutic interventions. Hence, several studies have used transcriptomics, proteomics, and metabolomics technologies, both independently and in combination, in an attempt to discover potential early markers of liver injury. However, the casual relationship between these observations and their relation to the APAP mechanism of liver toxicity are not clearly understood. Here, we used Sprague-Dawley rats orally gavaged with a single dose of 2 g/kg of APAP to collect tissue samples from the liver and kidney for transcriptomic analysis and plasma and urine samples for metabolomic analysis. We developed and used a multi-tissue, metabolism-based modeling approach to integrate these data, characterize the effect of excess APAP levels on liver metabolism, and identify a panel of plasma and urine metabolites that are associated with APAP-induced liver toxicity. Our analyses, which indicated that pathways involved in nucleotide-, lipid-, and amino acid-related metabolism in the liver were most strongly affected within 10 h following APAP treatment, identified a list of potential metabolites in these pathways that could serve as plausible markers of APAP-induced liver injury. Our approach identifies toxicant-induced changes in endogenous metabolism, is applicable to other toxicants based on transcriptomic data, and provides a mechanistic framework for interpreting metabolite alterations.

Comparison of the physiological relevance of systemic vs. portal insulin delivery to evaluate whole body glucose flux during an insulin clamp.

To understand the underlying pathology of metabolic diseases, such as diabetes, an accurate determination of whole body glucose flux needs to be made by a method that maintains key physiological features. One such feature is a positive differential in insulin concentration between the portal venous and systemic arterial circulation (P/S-IG). P/S-IG during the determination of the relative contribution of liver and extra-liver tissues/organs to whole body glucose flux during an insulin clamp with either systemic (SID) or portal (PID) insulin delivery was examined with insulin infusion rates of 1, 2, and 5 mU·kg(-1)·min(-1) under either euglycemic or hyperglycemic conditions in 6-h-fasted conscious normal rats. A P/S-IG was initially determined with endogenous insulin secretion to exist with a value of 2.07. During an insulin clamp, while inhibiting endogenous insulin secretion by somatostatin, P/S-IG remained at 2.2 with PID, whereas, P/S-IG disappeared completely with SID, which exhibited higher arterial and lower portal insulin levels compared with PID. Consequently, glucose disappearance rates and muscle glycogen synthetic rates were higher, but suppression of endogenous glucose production and liver glycogen synthetic rates were lower with SID compared with PID. When the insulin clamp was performed with SID at 2 and 5 mU·kg(-1)·min(-1) without managing endogenous insulin secretion under euglycemic but not hyperglycemic conditions, endogenous insulin secretion was completely suppressed with SID, and the P/S-IG disappeared. Thus, compared with PID, an insulin clamp with SID underestimates the contribution of liver in response to insulin to whole body glucose flux.

Glucotoxicity targets hepatic glucokinase in Zucker diabetic fatty rats, a model of type 2 diabetes associated with obesity.

A loss of glucose effectiveness to suppress hepatic glucose production as well as increase hepatic glucose uptake and storage as glycogen is associated with a defective increase in glucose phosphorylation catalyzed by glucokinase (GK) in Zucker diabetic fatty (ZDF) rats. We extended these observations by investigating the role of persistent hyperglycemia (glucotoxicity) in the development of impaired hepatic GK activity in ZDF rats. We measured expression and localization of GK and GK regulatory protein (GKRP), translocation of GK, and hepatic glucose flux in response to a gastric mixed meal load (MMT) and hyperglycemic hyperinsulinemic clamp after 1 or 6 wk of treatment with the sodium-glucose transporter 2 inhibitor (canaglifrozin) that was used to correct the persistent hyperglycemia of ZDF rats. Defective augmentation of glucose phosphorylation in response to a rise in plasma glucose in ZDF rats was associated with the coresidency of GKRP with GK in the cytoplasm in the midstage of diabetes, which was followed by a decrease in GK protein levels due to impaired posttranscriptional processing in the late stage of diabetes. Correcting hyperglycemia from the middle diabetic stage normalized the rate of glucose phosphorylation by maintaining GK protein levels, restoring normal nuclear residency of GK and GKRP under basal conditions and normalizing translocation of GK from the nucleus to the cytoplasm, with GKRP remaining in the nucleus in response to a rise in plasma glucose. This improved the liver's metabolic ability to respond to hyperglycemic hyperinsulinemia. Glucotoxicity is responsible for loss of glucose effectiveness and is associated with altered GK regulation in the ZDF rat.

Central nervous system neuropeptide Y signaling via the Y1 receptor partially dissociates feeding behavior from lipoprotein metabolism in lean rats.

Elevated plasma triglyceride (TG) levels contribute to an atherogenic dyslipidemia that is associated with obesity, diabetes, and metabolic syndrome. Numerous models of obesity are characterized by increased central nervous system (CNS) neuropeptide Y (NPY) tone that contributes to excess food intake and obesity. Previously, we demonstrated that intracerebroventricular (icv) administration of NPY in lean fasted rats also elevates hepatic production of very low-density lipoprotein (VLDL)-TG. Thus, we hypothesize that elevated CNS NPY action contributes to not only the pathogenesis of obesity but also dyslipidemia. Here, we sought to determine whether the effects of NPY on feeding and/or obesity are dissociable from effects on hepatic VLDL-TG secretion. Pair-fed, icv NPY-treated, chow-fed Long-Evans rats develop hypertriglyceridemia in the absence of increased food intake and body fat accumulation compared with vehicle-treated controls. We then modulated CNS NPY signaling by icv injection of selective NPY receptor agonists and found that Y1, Y2, Y4, and Y5 receptor agonists all induced hyperphagia in lean, ad libitum chow-fed Long-Evans rats, with the Y2 receptor agonist having the most pronounced effect. Next, we found that at equipotent doses for food intake NPY Y1 receptor agonist had the most robust effect on VLDL-TG secretion, a Y2 receptor agonist had a modest effect, and no effect was observed for Y4 and Y5 receptor agonists. These findings, using selective agonists, suggest the possibility that the effect of CNS NPY signaling on hepatic VLDL-TG secretion may be relatively dissociable from effects on feeding behavior via the Y1 receptor.

The pancreatic islet ß-cell-enriched transcription factor Pdx-1 regulates Slc30a8 gene transcription through an intronic enhancer.

The SLC30A8 gene encodes the zinc transporter ZnT-8, which provides zinc for insulin-hexamer formation. Genome-wide association studies have shown that a polymorphic variant in SLC30A8 is associated with altered susceptibility to Type 2 diabetes and we recently reported that glucose-stimulated insulin secretion is decreased in islets isolated from Slc30a8-knockout mice. The present study examines the molecular basis for the islet-specific expression of Slc30a8. VISTA analyses identified two conserved regions in Slc30a8 introns 2 and 3, designated enhancers A and B respectively. Transfection experiments demonstrated that enhancer B confers elevated fusion gene expression in both ßTC-3 cells and αTC-6 cells. In contrast, enhancer A confers elevated fusion gene expression selectively in ßTC-3 and not αTC-6 cells. These data suggest that enhancer A is an islet ß-cell-specific enhancer and that the mechanisms controlling Slc30a8 expression in α- and ß-cells are overlapping, but distinct. Gel retardation and ChIP (chromatin immunoprecipitation) assays revealed that the islet-enriched transcription factor Pdx-1 binds enhancer A in vitro and in situ respectively. Mutation of two Pdx-1-binding sites in enhancer A markedly reduces fusion gene expression suggesting that this factor contributes to Slc30a8 expression in ß-cells, a conclusion consistent with developmental studies showing that restriction of Pdx-1 to pancreatic islet ß-cells correlates with the induction of Slc30a8 gene expression and ZnT-8 protein expression in vivo.

Impact of a glycogen phosphorylase inhibitor and metformin on basal and glucagon-stimulated hepatic glucose flux in conscious dogs.

The effects of a glycogen phosphorylase inhibitor (GPI) and metformin (MT) on hepatic glucose fluxes (µmol · kg(-1) · min(-1)) in the presence of basal and 4-fold basal levels of plasma glucagon were investigated in 18-h fasted conscious dogs. Compared with the vehicle treatment, GPI infusion suppressed net hepatic glucose output (NHGO) completely (-3.8 ± 1.3 versus 9.9 ± 2.8) despite increased glucose 6-phosphate (G-6-P) neogenesis from gluconeogenic precursors (8.1 ± 1.1 versus 5.5 ± 1.1). MT infusion did not alter those parameters. In response to a 4-fold rise in plasma glucagon levels, in the vehicle group, plasma glucose levels were increased 2-fold, and NHGO was increased (43.9 ± 5.7 at 10 min and 22.7 ± 3.4 at steady state) without altering G-6-P neogenesis (3.7 ± 1.5 and 5.5 ± 0.5, respectively). In the GPI group, there was no increase in NHGO due to decreased glucose-6-phosphatase flux associated with reduced G-6-P concentration. A lower G-6-P concentration was the result of increased net glycogenesis without altering G-6-P neogenesis. In the MT group, the increment in NHGO (22.2 ± 4.4 at 10 min and 12.1 ± 3.6 at steady state) was approximately half of that of the vehicle group. The lesser NHGO was associated with reduced glucose-6-phosphatase flux but a rise in G-6-P concentration and only a small incorporation of plasma glucose into glycogen. In conclusion, the inhibition of glycogen phosphorylase a activity decreases basal and glucagon-induced NHGO via redirecting glucose 6-phosphate flux from glucose toward glycogen, and MT decreases glucagon-induced NHGO by inhibiting glucose-6-phosphatase flux and thereby reducing glycogen breakdown.

Deletion of the mouse Slc30a8 gene encoding zinc transporter-8 results in impaired insulin secretion.

The Slc30a8 gene encodes the islet-specific zinc transporter ZnT-8, which provides zinc for insulin-hexamer formation. Polymorphic variants in amino acid residue 325 of human ZnT-8 are associated with altered susceptibility to Type 2 diabetes and ZnT-8 autoantibody epitope specificity changes in Type 1 diabetes. To assess the physiological importance of ZnT-8, mice carrying a Slc30a8 exon 3 deletion were analysed histologically and phenotyped for energy metabolism and pancreatic hormone secretion. No gross anatomical or behavioural changes or differences in body weight were observed between wild-type and ZnT-8-/- mice, and ZnT-8-/- mouse islets were indistinguishable from wild-type in terms of their numbers, size and cellular composition. However, total zinc content was markedly reduced in ZnT-8-/- mouse islets, as evaluated both by Timm's histochemical staining of pancreatic sections and direct measurements in isolated islets. Blood glucose levels were unchanged in 16-week-old, 6 h fasted animals of either gender; however, plasma insulin concentrations were reduced in both female (approximately 31%) and male (approximately 47%) ZnT-8-/- mice. Intraperitoneal glucose tolerance tests demonstrated no impairment in glucose clearance in male ZnT-8-/- mice, but glucose-stimulated insulin secretion from isolated islets was reduced approximately 33% relative to wild-type littermates. In summary, Slc30a8 gene deletion is accompanied by a modest impairment in insulin secretion without major alterations in glucose metabolism.

A toxicogenomic approach to assess kidney injury induced by mercuric chloride in rats.

Kidney injury caused by disease, trauma, environmental exposures, or drugs may result in decreased renal function, chronic kidney disease, or acute kidney failure. Diagnosis of kidney injury using serum creatinine levels, a common clinical test, only identifies renal dysfunction after the kidneys have undergone severe damage. Other indicators sensitive to kidney injury, such as the level of urine kidney injury molecule-1 (KIM-1), lack the ability to differentiate between injury phenotypes. To address early detection as well as detailed categorization of kidney-injury phenotypes in preclinical animal or cellular studies, we previously identified eight sets (modules) of co-expressed genes uniquely associated with kidney histopathology. Here, we used mercuric chloride (HgCl2)-a model nephrotoxicant-to chemically induce kidney injuries as monitored by KIM-1 levels in Sprague Dawley rats at two doses (0.25 or 0.50 mg/kg) and two exposure lengths (10 or 34 h). We collected whole transcriptome RNA-seq data derived from five animals at each dose and time point to perform a toxicogenomics analysis. Consistent with documented injury phenotypes for HgCl2 toxicity, our kidney-injury-module approach identified the onset of necrosis and dilation as early as 10 h after a dose of 0.50 mg/kg that produced only mild injury as judged by urinary KIM-1 excretion. The results of these animal studies highlight the potential of the kidney-injury-module approach to provide a sensitive and histopathology-specific readout of renal toxicity.

Genome-Scale Model-Based Identification of Metabolite Indicators for Early Detection of Kidney Toxicity.

Identifying early indicators of toxicant-induced organ damage is critical to provide effective treatment. To discover such indicators and the underlying mechanisms of toxicity, we used gentamicin as an exemplar kidney toxicant and performed systematic perturbation studies in Sprague Dawley rats. We obtained high-throughput data 7 and 13 h after administration of a single dose of gentamicin (0.5 g/kg) and identified global changes in genes in the liver and kidneys, metabolites in the plasma and urine, and absolute fluxes in central carbon metabolism. We used these measured changes in genes in the liver and kidney as constraints to a rat multitissue genome-scale metabolic network model to investigate the mechanism of gentamicin-induced kidney toxicity and identify metabolites associated with changes in tissue gene expression. Our experimental analysis revealed that gentamicin-induced metabolic perturbations could be detected as early as 7 h postexposure. Our integrated systems-level analyses suggest that changes in kidney gene expression drive most of the significant metabolite alterations in the urine. The analyses thus allowed us to identify several significantly enriched injury-specific pathways in the kidney underlying gentamicin-induced toxicity, as well as metabolites in these pathways that could serve as potential early indicators of kidney damage.

Mechanism-based identification of plasma metabolites associated with liver toxicity.

Early diagnosis of liver injuries caused by drugs or occupational exposures is necessary to enable effective treatments and prevent liver failure. Whereas histopathology remains the gold standard for assessing hepatotoxicity in animals, plasma aminotransferase levels are the primary measures for monitoring liver dysfunction in humans. In this study, using Sprague Dawley rats, we investigated whether integrated analyses of transcriptomic and metabolomic data with genome-scale metabolic models (GSMs) could identify early indicators of injury and provide new insights into the mechanisms of hepatotoxicity. We obtained concurrent measurements of gene-expression changes in the liver and kidneys, and expression changes along with metabolic profiles in the plasma and urine, from rats 5 or 10 h after exposing them to one of two classical hepatotoxicants, acetaminophen (2 g/kg) or bromobenzene (0.4 g/kg). Global multivariate analyses revealed that gene-expression changes in the liver and metabolic profiles in the plasma and urine of toxicant-treated animals differed from those of controls, even at time points much earlier than changes detected by conventional markers of liver injury. Furthermore, clustering analysis revealed that both the gene-expression changes in the liver and the metabolic profiles in the plasma induced by the two hepatotoxicants were highly correlated, indicating commonalities in the liver toxicity response. Systematic GSM-based analyses yielded metabolites associated with the mechanisms of toxicity and identified several lipid and amino acid metabolism pathways that were activated by both toxicants and those uniquely activated by each. Our findings suggest that several metabolite alterations, which are strongly associated with the mechanisms of toxicity and occur within injury-specific pathways (e.g., of bile acid and fatty acid metabolism), could be targeted and clinically assessed for their potential as early indicators of liver damage.

Assessing Chemical-Induced Liver Injury In Vivo From In Vitro Gene Expression Data in the Rat: The Case of Thioacetamide Toxicity.

Consumers are exposed to thousands of chemicals with potentially adverse health effects. However, these chemicals will never be tested for toxicity because of the immense resources needed for animal-based (in vivo) toxicological studies. Today, there are no viable in vitro alternatives to these types of animal studies. To develop an in vitro approach, we investigated whether we could predict in vivo organ injuries in rats with the use of RNA-seq data acquired from tissues early in the development of toxicant-induced injury, by comparing gene expression data from RNA isolated from these rat tissues with those obtained from in vitro exposure of primary liver and kidney cells. We collected RNA-seq data from the liver and kidney tissues of Sprague-Dawley rats 8 or 24 h after exposing them to vehicle (control), low (25 mg/kg), or high (100 mg/kg) doses of thioacetamide, a known liver toxicant that promotes fibrosis; we used these doses and exposure times to cause only mild toxicant-induced injury. For the in vitro study, we treated two cell types from Sprague-Dawley rats, primary hepatocytes (vehicle; low, 0.025 mM; or high, 0.125 mM dose), and renal tube epithelial cells (vehicle; low, 0.125 mM; or high, 0.500 mM) dose) with the thioacetamide metabolite, thioacetamide-S-oxide, selecting in vitro doses and exposure times to recreate the early-stage toxicant-induced injury model that we achieved in vivo. RNA-seq data were collected 9 or 24 h after application of vehicle or thioacetamide-S-oxide. We found that our modular approach for the analysis of gene expression data derived from in vivo RNA-seq strongly correlated (R2 > 0.6) with the in vitro results at two different dose levels of thioacetamide/thioacetamide-S-oxide after 24 h of exposure. The top-ranked liver injury modules in vitro correctly identified the ensuing development of liver fibrosis.

Network Modeling of Liver Metabolism to Predict Plasma Metabolite Changes During Short-Term Fasting in the Laboratory Rat.

The liver-a central metabolic organ that integrates whole-body metabolism to maintain glucose and fatty-acid regulation, and detoxify ammonia-is susceptible to injuries induced by drugs and toxic substances. Although plasma metabolite profiles are increasingly investigated for their potential to detect liver injury earlier than current clinical markers, their utility may be compromised because such profiles are affected by the nutritional state and the physiological state of the animal, and by contributions from extrahepatic sources. To tease apart the contributions of liver and non-liver sources to alterations in plasma metabolite profiles, here we sought to computationally isolate the plasma metabolite changes originating in the liver during short-term fasting. We used a constraint-based metabolic modeling approach to integrate central carbon fluxes measured in our study, and physiological flux boundary conditions gathered from the literature, into a genome-scale model of rat liver metabolism. We then measured plasma metabolite profiles in rats fasted for 5-7 or 10-13 h to test our model predictions. Our computational model accounted for two-thirds of the observed directions of change (an increase or decrease) in plasma metabolites, indicating their origin in the liver. Specifically, our work suggests that changes in plasma lipid metabolites, which are reliably predicted by our liver metabolism model, are key features of short-term fasting. Our approach provides a mechanistic model for identifying plasma metabolite changes originating in the liver.