Population-based dose-response analysis of liver transcriptional response to trichloroethylene in mouse

Intra- and Inter-Species Variability in Urinary N7-(1-Hydroxy-3-buten-2-yl)guanine Adducts Following Inhalation Exposure to 1,3-Butadiene

1,3-Butadiene is a known carcinogen primarily targeting lymphoid tissues, lung, and liver. Cytochrome P450 activates butadiene to epoxides which form covalent DNA adducts that are thought to be a key mechanistic event in cancer. Previous studies suggested that inter-species, -tissue, and -individual susceptibility to adverse health effects of butadiene exposure may be due to differences in metabolism and other mechanisms. In this study, we aimed to examine the extent of inter-individual and inter-species variability in the urinary N7-(1-hydroxy-3-buten-2-yl)guanine (EB-GII) DNA adduct, a well-known biomarker of exposure to butadiene. For a population variability study in mice, we used the collaborative cross model. Female and male mice from five strains were exposed to filtered air or butadiene (590 ppm, 6 h/day, 5 days/week for 2 weeks) by inhalation. Urine samples were collected, and the metabolic activation of butadiene by DNA-reactive species was quantified as urinary EB-GII adducts. We quantified the degree of EB-GII variation across mouse strains and sexes; then, we compared this variation with the data from rats (exposed to 62.5 or 200 ppm butadiene) and humans (0.004-2.2 ppm butadiene). We show that sex and strain are significant contributors to the variability in urinary EB-GII levels in mice. In addition, we find that the degree of variability in urinary EB-GII in collaborative cross mice, when expressed as an uncertainty factor for the inter-individual variability (UF_H), is relatively modest (≤threefold) possibly due to metabolic saturation. By contrast, the variability in urinary EB-GII (adjusted for exposure) observed in humans, while larger than the default value of 10-fold, is largely consistent with UF_H estimates for other chemicals based on human data for non-cancer endpoints. Overall, these data demonstrate that urinary EB-GII levels, particularly from human studies, may be useful for quantitative characterization of human variability in cancer risks to butadiene.

Sex-dependent effects of preconception exposure to arsenite on gene transcription in parental germ cells and on transcriptomic profiles and diabetic phenotype of offspring

Chronic exposure to inorganic arsenic (iAs) has been linked to diabetic phenotypes in both humans and mice. However, diabetogenic effects of iAs exposure during specific developmental windows have never been systematically studied. We have previously shown that in mice, combined preconception and in utero exposures to iAs resulted in impaired glucose homeostasis in male offspring. The goal of the present study was to determine if preconception exposure alone can contribute to this outcome. We have examined metabolic phenotypes in male and female offspring from dams and sires that were exposed to iAs in drinking water (0 or 200 μg As/L) for 10 weeks prior to mating. The effects of iAs exposure on gene expression profiles in parental germ cells, and pancreatic islets and livers from offspring were assessed using RNA sequencing. We found that iAs exposure significantly altered transcript levels of genes, including diabetes-related genes, in the sperm of sires. Notably, some of the same gene transcripts and the associated pathways were also altered in the liver of the offspring. The exposure had a more subtle effect on gene expression in maternal oocytes and in pancreatic islets of the offspring. In female offspring, the preconception exposure was associated with increased adiposity, but lower blood glucose after fasting and after glucose challenge. HOMA-IR, the indicator of insulin resistance, was also lower. In contrast, the preconception exposure had no effects on blood glucose measures in male offspring. However, males from parents exposed to iAs had higher plasma insulin after glucose challenge and higher insulinogenic index than control offspring, indicating a greater requirement for insulin to maintain glucose homeostasis. Our results suggest that preconception exposure may contribute to the development of diabetic phenotype in male offspring, possibly mediated through germ cell-associated inheritance. Future research can investigate role of epigenetics in this phenomenon. The paradoxical outcomes in female offspring, suggesting a protective effect of the preconception exposure, warrant further investigation.

Using Collaborative Cross Mouse Population to Fill Data Gaps in Risk Assessment: A Case Study of Population-Based Analysis of Toxicokinetics and Kidney Toxicodynamics of Tetrachloroethylene

BACKGROUND

Interindividual variability in susceptibility remains poorly characterized for environmental chemicals such as tetrachloroethylene (PERC). Development of population-based experimental models provide a potential approach to fill this critical need in human health risk assessment.

OBJECTIVES

In this study, we aimed to better characterize the contribution of glutathione (GSH) conjugation to kidney toxicity of PERC and the degree of associated interindividual toxicokinetic (TK) and toxicodynamic (TD) variability by using the Collaborative Cross (CC) mouse population.

METHODS

Male mice from 45 strains were intragastrically dosed with PERC ([Formula: see text]) or vehicle (5% Alkamuls EL-620 in saline), and time-course samples were collected for up to 24 h. Population variability in TK of S-(1,2,2-trichlorovinyl)GSH (TCVG), S-(1,2,2-trichlorovinyl)-L-cysteine (TCVC), and N-acetyl-S-(1,2,2-trichlorovinyl)-L-cysteine (NAcTCVC) was quantified in serum, liver, and kidney, and analyzed using a toxicokinetic model. Effects of PERC on kidney weight, fatty acid metabolism-associated genes [ Acot1 (Acyl-CoA thioesterase 1), Fabp1 (fatty acid-binding protein 1), and Ehhadh (enoyl-coenzyme A, hydratase/3-hydroxyacyl coenzyme A dehydrogenase)], and a marker of proximal tubular injury [KIM-1 (kidney injury molecule-1)/Hepatitis A virus cellular receptor 1 ( Havcr1)] were evaluated. Finally, quantitative data on interstrain variability in both formation of GSH conjugation metabolites of PERC and its kidney effects was used to calculate adjustment factors for the interindividual variability in both TK and TD.

RESULTS

Mice treated with PERC had significantly lower kidney weight, higher kidney-to-body weight (BW) ratio, and higher expression of fatty acid metabolism-associated genes ( Acot1, Fabp1, and Ehhadh) and a marker of proximal tubular injury (KIM-1/ Havcr1). Liver levels of TCVG were significantly correlated with KIM-1/ Havcr1 in kidney, consistent with kidney injury being associated with GSH conjugation. We found that the default uncertainty factor for human variability may be marginally adequate to protect 95%, but not more, of the population for kidney toxicity mediated by PERC.

DISCUSSION

Overall, this study demonstrates the utility of the CC mouse population in characterizing metabolism-toxicity interactions and quantifying interindividual variability. Further refinement of the characterization of interindividual variability can be accomplished by incorporating these data into in silico population models both for TK (such as a physiologically based pharmacokinetic model), as well as for toxicodynamic responses. https://doi.org/10.1289/EHP5105.

Population-Based Analysis of DNA Damage and Epigenetic Effects of 1,3-Butadiene in the Mouse

Metabolism of 1,3-butadiene, a known human and rodent carcinogen, results in formation of reactive epoxides, a key event in its carcinogenicity. Although mice exposed to 1,3-butadiene present DNA adducts in all tested tissues, carcinogenicity is limited to liver, lung, and lymphoid tissues. Previous studies demonstrated that strain- and tissue-specific epigenetic effects in response to 1,3-butadiene exposure may influence susceptibly to DNA damage and serve as a potential mechanism of tissue-specific carcinogenicity. This study aimed to investigate interindividual variability in the effects of 1,3-butadiene using a population-based mouse model. Male mice from 20 Collaborative Cross strains were exposed to 0 or 635 ppm 1,3-butadiene by inhalation (6 h/day, 5 days/week) for 2 weeks. We evaluated DNA damage and epigenetic effects in target (lung and liver) and nontarget (kidney) tissues of 1,3-butadiene-induced carcinogenesis. DNA damage was assessed by measuring N-7-(2,3,4-trihydroxybut-1-yl)-guanine (THB-Gua) adducts. To investigate global histone modification alterations, we evaluated the trimethylation and acetylation of histones H3 and H4 across tissues. Changes in global cytosine DNA methylation were evaluated from the levels of methylation of LINE-1 and SINE B1 retrotransposons. We quantified the degree of variation across strains, deriving a chemical-specific human variability factor to address population variability in carcinogenic risk, which is largely ignored in current cancer risk assessment practice. Quantitative trait locus mapping identified four candidate genes related to chromatin remodeling whose variation was associated with interstrain susceptibility. Overall, this study uses 1,3-butadiene to demonstrate how the Collaborative Cross mouse population can be used to identify the mechanisms for and quantify the degree of interindividual variability in tissue-specific effects that are relevant to chemically induced carcinogenesis.

Bradykinin contributes to immune liver injury via B2R receptor-mediated pathways in trichloroethylene sensitized mice: A role in Kupffer cell activation

We have previously shown trichloroethylene (TCE) induced occupational medicamentosa-like dermatitis due to TCE (OMLDT) with immune liver injury, and kallikrein-kinin system (KKS) activation as a probably mechanism underlying the immune damage. Bradykinin (BK) is an important active component of KKS system function, but the specific role of BK in the immune liver injury has never been examined. The present study aimed to explore the important role of BK and mechanisms of action in immune liver injury induced by TCE. TCE sensitization significantly increased the expression of BK receptor (B2R) in the liver. Compared to blank and vehicle control group, TCE sensitization positive mice developed exacerbated liver injury evidenced by elevated AST, ALT levels and hepatocyte damage. TCE sensitization also stimulated MAPK and STAT3 activation in liver tissue. B2R antagonist HOE140 ameliorated these changes. Kupffer cells (KCs) of the liver were also activated following TCE sensitization; both CD68⁺ KCs and CD16/CD32⁺ M1 type KCs were increased in TCE positive group. Further experiments isolated the KCs from the liver in each group and showed that TCE sensitization resulted activation of MAPK signal pathway which in turn caused release of the pro-inflammatory cytokines, IL-1β, IL-6, TNF-α, in KCs; the antagonist HOE140 again decreased these changes in KCs. These results uncover a novel role of BK and B2R cross-talk in KCs activation in TCE sensitized mice, mediated by pro-inflammatory cytokine release via MAPK and STAT3 activation, contributing to the immune liver injury.

Multicellular Models Bridging Intracellular Signaling and Gene Transcription to Population Dynamics

Cell signaling and gene transcription occur at faster time scales compared to cellular death, division, and evolution. Bridging these multiscale events in a model is computationally challenging. We introduce a framework for the systematic development of multiscale cell population models. Using message passing interface (MPI) parallelism, the framework creates a population model from a single-cell biochemical network model. It launches parallel simulations on a single-cell model and treats each stand-alone parallel process as a cell object. MPI mediates cell-to-cell and cell-to-environment communications in a server-client fashion. In the framework, model-specific higher level rules link the intracellular molecular events to cellular functions, such as death, division, or phenotype change. Cell death is implemented by terminating a parallel process, while cell division is carried out by creating a new process (daughter cell) from an existing one (mother cell). We first demonstrate these capabilities by creating two simple example models. In one model, we consider a relatively simple scenario where cells can evolve independently. In the other model, we consider interdependency among the cells, where cellular communication determines their collective behavior and evolution under a temporally evolving growth condition. We then demonstrate the framework’s capability by simulating a full-scale model of bacterial quorum sensing, where the dynamics of a population of bacterial cells is dictated by the intercellular communications in a time-evolving growth environment.

Comparative analysis of metabolism of trichloroethylene and tetrachloroethylene among mouse tissues and strains

Trichloroethylene (TCE) and tetrachloroethylene (PCE) are structurally similar chemicals that are metabolized through oxidation and glutathione conjugation pathways. Both chemicals have been shown to elicit liver and kidney toxicity in rodents and humans; however, TCE has been studied much more extensively in terms of both metabolism and toxicity. Despite their qualitative similarities, quantitative comparison of tissue- and strain-specific metabolism of TCE and PCE has not been performed. To fill this gap, we conducted a comparative toxicokinetic study where equimolar single oral doses of TCE (800 mg/kg) or PCE (1000 mg/kg) were administered to male mice of C57BL/6J, B6C3F1/J, and NZW/LacJ strains. Samples of liver, kidney, serum, brain, and lung were obtained for up to 36 h after dosing. For each tissue, concentrations of parent compounds, as well as their oxidative and glutathione conjugation metabolites were measured and concentration-time profiles constructed. A multi-compartment toxicokinetic model was developed to quantitatively compare TCE and PCE metabolism. As expected, the flux through oxidation metabolism pathway predominated over that through conjugation across all mouse strains examined, it is 1,200-3,800 fold higher for TCE and 26-34 fold higher for PCE. However, the flux through glutathione conjugation, albeit a minor metabolic pathway, was 21-fold higher for PCE as compared to TCE. The degree of inter-strain variability was greatest for oxidative metabolites in TCE-treated and for glutathione conjugation metabolites in PCE-treated mice. This study provides critical data for quantitative comparisons of TCE and PCE metabolism, and may explain the differences in organ-specific toxicity between these structurally similar chemicals.

gQTL: A Web Application for QTL Analysis Using the Collaborative Cross Mouse Genetic Reference Population

Multi-parental recombinant inbred populations, such as the Collaborative Cross (CC) mouse genetic reference population, are increasingly being used for analysis of quantitative trait loci (QTL). However specialized analytic software for these complex populations is typically built in R that works only on command-line, which limits the utility of these powerful resources for many users. To overcome analytic limitations, we developed gQTL, a web accessible, simple graphical user interface application based on the DOQTL platform in R to perform QTL mapping using data from CC mice.