Incorporating population-level genetic variability within laboratory models in toxicology: From the individual to the population

Transgenerational effects of triazole fungicides on gene expression and egg compounds in non-exposed offspring: A case study using Red-Legged Partridges (Alectoris rufa)

Triazole fungicides are widely used to treat cereal seeds before sowing. Granivorous birds like the Red-legged Partridge (Alectoris rufa) have high exposure risk because they ingest treated seeds that remain on the field surface. As triazole fungicides can act as endocrine disruptors, affecting sterol synthesis and reproduction in birds several months after exposure, we hypothesized that these effects could also impact subsequent generations of exposed birds. To test this hypothesis, we exposed adult partridges (F0) to seeds treated at commercial doses with four different formulations containing triazoles as active ingredients (flutriafol, prothioconazole, tebuconazole, and a mixture of the latter two), simulating field exposure during late autumn sowing. During the subsequent reproductive season, two to four months after exposure, we examined compound allocation of steroid hormones, cholesterol, vitamins, and carotenoids in eggs laid by exposed birds (F1), as well as the expression of genes encoding enzymes involved in sterol biosynthesis in one-day-old chicks of this F1. One year later, F1 animals were paired again to investigate the expression of the same genes in the F2 chicks. We found changes in the expression of some genes for all treatments and both generations. Additionally, we observed an increase in estrone levels in eggs from partridges treated with flutriafol compared to controls, a decrease in tocopherol levels in partridges exposed to the mixture of tebuconazole and prothioconazole, and an increase in retinol levels in partridges exposed to prothioconazole. Despite sample size limitations, this study provides novel insights into the mechanisms of action of the previously observed effects of triazole fungicide-treated seeds on avian reproduction with evidence that the effects can persist beyond the exposure windows, affecting unexposed offspring of partridges fed with treated seeds. The results highlight the importance of considering long-term chronic effects when assessing pesticide risks to wild birds.

Utilizing a Population-Genetic Framework to Test for Gene-Environment Interactions between Zebrafish Behavior and Chemical Exposure

Individuals within genetically diverse populations display broad susceptibility differences upon chemical exposures. Understanding the role of gene-environment interactions (GxE) in differential susceptibility to an expanding exposome is key to protecting public health. However, a chemical's potential to elicit GxE is often not considered during risk assessment. Previously, we've leveraged high-throughput zebrafish (Danio rerio) morphology screening data to reveal patterns of potential GxE effects. Here, using a population genetics framework, we apportioned variation in larval behavior and gene expression in three different PFHxA environments via mixed-effect modeling to assess significance of GxE term. We estimated the intraclass correlation (ICC) between full siblings from different families using one-way random-effects model. We found a significant GxE effect upon PFHxA exposure in larval behavior, and the ICC of behavioral responses in the PFHxA exposed population at the lower concentration was 43.7%, while that of the control population was 14.6%. Considering global gene expression data, a total of 3746 genes showed statistically significant GxE. By showing evidence that heritable genetics are directly affecting gene expression and behavioral susceptibility of individuals to PFHxA exposure, we demonstrate how standing genetic variation in a heterogeneous population such as ours can be leveraged to test for potential GxE.

Genetic background and sex control the outcome of high-fat diet feeding in mice

The sharp increase in obesity prevalence worldwide is mainly attributable to changes in physical activity and eating behavior but the metabolic and clinical impacts of these obesogenic conditions vary between sexes and genetic backgrounds. This warrants personalized treatments of obesity and its complications, which require a thorough understanding of the diversity of metabolic responses to high-fat diet intake. By analyzing nine genetically diverse mouse strains, we show that much like humans, mice exhibit a huge variety of physiological and biochemical responses to high-fat diet. The strains exhibit various degrees of alterations in their phenotypic makeup. At the transcriptome level, we observe dysregulations of immunity, translation machinery, and mitochondrial genes. At the biochemical level, the enzymatic activity of mitochondrial complexes is affected. The diversity across mouse strains, diets, and sexes parallels that found in humans and supports the use of diverse mouse populations in future mechanistic or preclinical studies on metabolic dysfunctions.

•

Strain- and sex-specific profile of metabolic dysfunction in mice

•

Liver mitochondrial complex activity in vivo associates with metabolic traits

•

Open data source for evaluating different mouse strains for metabolic disease

•

Interactive data exploration through an online application

Systemic review of genetic and epigenetic factors underlying differential toxicity to environmental lead (Pb) exposure

Lead (Pb) poisoning is a major public health concern in environmental justice communities of the USA and in many developing countries. There is no identified safety threshold for lead in blood, as low-level Pb exposures can lead to severe toxicity in highly susceptible individuals and late onset of diseases from early-life exposure. However, identifying "susceptibility genes" or "early exposure biomarkers" remains challenging in human populations. There is a considerable variation in susceptibility to harmful effects from Pb exposure in the general population, likely due to the complex interplay of genetic and/or epigenetic factors. This systematic review summarizes current state of knowledge on the role of genetic and epigenetic factors in determining individual susceptibility in response to environmental Pb exposure in humans and rodents. Although a number of common genetic and epigenetic factors have been identified, the reviewed studies, which link these factors to various adverse health outcomes following Pb exposure, have provided somewhat inconsistent evidence of main health effects. Acknowledging the compelling need for new approaches could guide us to better characterize individual responses, predict potential adverse outcomes, and identify accurate and usable biomarkers for Pb exposure to improve mitigation therapies to reduce future adverse health outcomes of Pb exposure.

Human Chemical Exposure from Background Emissions in the United States and the Implication for Quantifying Risks from Marginal Emission Increase

The linear dose-response relationship has long been assumed in assessments of health risk from an incremental chemical emission relative to background emissions. In this study, we systematically examine the relevancy of such an assumption with real-world data. We used the reported emission data, as background emissions, from the 2017 U.S. National Emission Inventory for 95 organic chemicals to estimate the central tendencies of exposures of the general U.S. population. Previously published nonlinear dose-response relationships for chemicals were used to estimate health risk from exposure. We also explored and identified four intervals of exposure in which the nonlinear dose-response relationship may be linearly approximated with fixed slopes. Predicted rates of exposure to these 95 chemicals are all within the lowest of the four intervals and associated with low health risk. The health risk may be overestimated if a slope on the dose-response relationship extrapolated from toxicological assays based on high response rates is used for a marginal increase in emission not substantially higher than background emissions. To improve the confidence of human health risk estimates for chemicals, future efforts should focus on deriving a more accurate dose-response relationship at lower response rates and interface it with exposure assessments.

CRISPR approach in environmental chemical screening focusing on population variability

A significant barrier to include population variability in risk assessment is our incomplete understanding of inter-individual variability and the differential susceptibility to environmental exposures induced adverse outcomes. By combining genome editing tools with the population diversity model, this article intended to highlight a potential strategy to identify and characterize the inter-individual variability factors, the determinant gene anchoring to a particular phenotype. The goal could be achieved by integrating the perturbed CRISPR-based unbiased functional genomics screening, genome-wide or a focused subset of genes, in a population-based in vitro model system (such as the lymphoblastoid cell lines, LCL, available from HapMap and 1000 Genomes project). Then data can be translated to genetic variability and individual (or subpopulation) susceptibility by incorporating ethnicity and corresponding genome-wide association studies (GWAS) with functional genomics screening results. This approach can provide complementary data for next-generation risk assessment, in particular, for environmental stressors. The current paper outlined the previous work conducted with a population-based in vitro model system, perturbed CRISPR-based functional toxicogenomic screening of environmental chemicals, and finally, the potential strategies to combine these two platforms with their opportunities and challenges to achieve a mechanistic understanding of population variability.

Host genetics and gut microbiota cooperatively contribute to azoxymethane-induced acute toxicity in Collaborative Cross mice

Azoxymethane (AOM) is a widely used carcinogen to study chemical-induced colorectal carcinogenesis and is an agent for studying fulminant hepatic failure. The inter-strain susceptibility to acute toxicity by AOM has been reported, but its association with host genetics or gut microbiota remains largely unexplored. Here a cohort of genetically diverse Collaborative Cross (CC) mice was used to assess the contribution of host genetics and the gut microbiome to AOM-induced acute toxicity. We observed variation in AOM-induced acute liver failure across CC strains. Quantitative trait loci (QTL) analysis revealed three chromosome regions significantly associated with AOM toxicity. Genes located within these QTL, including peroxisome proliferator-activated receptor alpha (Ppara), were enriched for enzyme activator and nucleoside-triphosphatase regulator activity. We further demonstrated that the protein level of PPARα in liver tissues from sensitive strains was remarkably lower compared to levels in resistant strains, consistent with protective role of PPAR family in liver injury. We discovered that the abundance levels of gut microbial families Anaeroplasmataceae, Ruminococcaceae, Lactobacillaceae, Akkermansiaceae and Clostridiaceae were significantly higher in the sensitive strains compared to the resistant strains. Using a random forest classifier method, we determined that the relative abundance levels of these microbial families predicted AOM toxicity with the area under the receiver-operating curve (AUC) of 0.75. Combining the three genetic loci and five microbial families increased the predictive accuracy of AOM toxicity (AUC of 0.99). Moreover, we found that Ruminococcaceae and Lactobacillaceae acted as mediators between host genetics and AOM toxicity. In conclusion, this study shows that host genetics and specific microbiome members play a critical role in AOM-induced acute toxicity, which provides a framework for analysis of the health effects from environmental toxicants.

Peanut butter as an alternative dose delivery method to prevent strain-dependent orogastric gavage-induced stress in mouse teratogenicity studies

INTRODUCTION

Animal-based studies are essential for assessing toxicity to environmental pollutants, especially when the potential targets are specific developmental time points, teratogenic, or multi-organ systems that cannot be modeled in vitro. Orogastric gavage is a widely used technique for exposure because of its increased accuracy of dose administration over free feeding. However, repeated use of this method has been reported to cause physiological stress on the exposed animals that could interfere with interpretation of results. Previous studies have shown that genetic background also contributes to the level of stress and can affect individual response.

METHODS

To evaluate the impact of stress on repeated orogastric gavage, we exposed C67BL/6J and 129S1/SvImJ inbred mouse strains to 2,3,7,8 tetrachlorodibenzo-p-dioxin (TCDD), a potent xenobiotic that has been extensively studied in vivo. Pregnant females were dosed for ten days after mating using orogastric gavage with olive oil as vehicle or through diet using peanut butter as vehicle. Serum corticosterone levels, body weight, and reproduction endpoints were measured to evaluate levels of stress induced by the dosing technique.

RESULTS

The levels of stress caused by orogastric gavage was strongly dependent on strain background and on the phenotypic endpoint. Orogastric gavage-induced stress was more detrimental in 129S1/SvlmJ pregnant female mice than in C57BL/6J.

CONCLUSION

These results show that administration of xenobiotics via controlled diet can improve the reproducibility and rigor of exposure studies requiring orogastric delivery.

Engineered tissues and strategies to overcome challenges in drug development

Current preclinical studies in drug development utilize high-throughput in vitro screens to identify drug leads, followed by both in vitro and in vivo models to predict lead candidates' pharmacokinetic and pharmacodynamic properties. The goal of these studies is to reduce the number of lead drug candidates down to the most likely to succeed in later human clinical trials. However, only 1 in 10 drug candidates that emerge from preclinical studies will succeed and become an approved therapeutic. Lack of efficacy or undetected toxicity represents roughly 75% of the causes for these failures, despite these parameters being the primary exclusion criteria in preclinical studies. Recently, advances in both biology and engineering have created new tools for constructing new preclinical models. These models can complement those used in current preclinical studies by helping to create more realistic representations of human tissues in vitro and in vivo. In this review, we describe current preclinical models to identify their value and limitations and then discuss select areas of research where improvements in preclinical models are particularly needed to advance drug development. Following this, we discuss design considerations for constructing preclinical models and then highlight recent advances in these efforts. Taken together, we aim to review the advances as of 2020 surrounding the prospect of biological and engineering tools for adding enhanced biological relevance to preclinical studies to aid in the challenges of failed drug candidates and the burden this poses on the drug development enterprise and thus healthcare.

Hepatocyte-like cells derived from human induced pluripotent stem cells using small molecules: implications of a transcriptomic study

Background

Hepatocyte-like cells (HLCs) derived from human induced pluripotent stem cells (iPSCs) hold great promise in toxicological applications as well as in regenerative medicine. Previous efforts on hepatocyte differentiation have mostly relied on the use of growth factors (GFs) to recapitulate developmental signals under in vitro conditions. Recently, the use of small molecules (SMs) has emerged as an attractive tool to induce cell fate transition due to its superiority in terms of both quality and cost. However, HLCs derived using SMs have not been well characterized, especially on the transcriptome level.

Methods

HLCs were differentiated from human iPSCs using a protocol that only involves SMs and characterized by transcriptomic analysis using whole genome microarrays.

Results

HLCs derived using the SM protocol (HLC_SM) displayed specific hepatic marker expression and demonstrated key hepatic functions. Transcriptomic analysis of the SM-driven differentiation defined a hepatocyte differentiation track and characterized the expression of some key marker genes in major stages of hepatocyte differentiation. In addition, HLC_SM were scored with CellNet, a bioinformatics tool quantifying how closely engineered cell populations resemble their target cell type, and compared to primary human hepatocytes (PHHs), adult liver tissue, fetal liver tissue, HLCs differentiated using GFs (HLC_GF), and commercially available HLCs. Similar to HLC_GF, HLC_SM displayed a mixed phenotype of fetal and adult hepatocytes and had relatively low expression of metabolic enzymes, transporters, and nuclear receptors compared to PHHs. Finally, the differentially expressed genes in HLC_SM compared to HLC_GF and to PHHs were analyzed to identify pathways and upstream transcription regulators which could potentially be manipulated to improve the differentiation of HLCs.

Conclusions

Overall, the present study demonstrated the usefulness of the SM-based hepatocyte differentiation method, offered new insights into the molecular basis of hepatogenesis and associated gene regulation, and suggested ways for further improvements in hepatocyte differentiation in order to obtain more mature HLCs that could be used in toxicological studies.

Toxic Responses Induced at High Doses May Affect Benchmark Doses

To derive reference points (RPs) for health-based guidance values, the benchmark dose (BMD) approach increasingly replaces the no-observed-adverse-effect level approach. In the BMD approach, the RP corresponds to the benchmark dose lower confidence bounds (BMDLs) of a mathematical dose-response model derived from responses of animals over the entire dose range applied. The use of the entire dose range is seen as an important advantage of the BMD approach. This assumes that responses over the entire dose range are relevant for modeling low-dose responses, the basis for the RP. However, if part of the high-dose response was unnoticed triggered by a mechanism of action (MOA) that does not work at low doses, the high-dose response distorts the modeling of low-dose responses. Hence, we investigated the effect of high-dose specific responses on BMDLs by assuming a low- and a high-dose MOA. The BMDLs resulting from modeling fictitious quantal data were scattered over a broad dose range overlapping with the toxic range. Hence, BMDLs are sensitive to high-dose responses even though they might be irrelevant to low-dose response modeling. When applying the BMD approach, care should be taken that high-dose specific responses do not unduly affect the BMDL that derives from low doses.

Accuracy, discriminative properties and reliability of a human ESC-based in vitro toxicity assay to distinguish teratogens responsible for neural tube defects

The poor correlation of developmental toxicity studies in animals with human outcome data has emphasized the need for complementary assays based on human cells and tissues. As neural tube defects represent an important proportion of congenital malformations, we evaluated here the accuracy of a human embryonic stem cell (hESC)-based assay to predict chemically induced disruption of neural tube formation. As teratogenic compounds, we used cyclopamine (CPA), valproic acid (VPA), ochratoxin A (OTA) and mycophenolic acid (MMF), all suspected or known inducers of human neural tube defects, as well as theophylline and saccharin as negative control compounds. We analyzed their effects on the ability of hES cells to give rise to neural precursors (expressing specific marker Nestin), to form neural tube-like structures (rosettes), and to express specific markers (Sox1, Otx2, Lix1, EvI1, Rspo3) during rosette formation. The results showed that various effects of the selected compounds on early neural development could be specifically revealed in vitro through related alterations of neurogenic differentiation of hESC. Furthermore, it was possible to discriminate toxicants acting at different time points during embryonic development and, therefore, responsible for distinct adverse effects on neural tube formation. By comparing four different hESC lines, we observed a significant (up to fivefold) variability of the line-dependent response to toxicants. We highlight at least two sources of variability: one related to the heterogeneity of hESC lines in culture (stemness/commitment profiles); the second to possible genetically determined differences in individual sensitivity to teratogens.

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.

Initial Assessment of Variability of Responses to Toxicants in Donor-Specific Endothelial Colony Forming Cells

There is increased interest in using high throughput in vitro assays to characterize human population variability in response to toxicants and drugs. Utilizing primary human endothelial colony-forming cells (ECFCs) isolated from blood would be highly useful for this purpose because these cells are involved in neonatal and adult vasculogenesis. We characterized the cytotoxicity of four known toxic chemicals (NaAsO2, CdCl2, tributyltin [TBT], and menadione) and their four relatively nontoxic counterparts (Na2HAsO4, ZnCl2, SnCl2, and phytonadione, respectively) in eight ECFC clones representing four neonatal donors (2 male and 2 female donors, 2 clones per donor). ECFCs were exposed to 9 concentrations of each chemical in duplicate; cell viability was evaluated 48 h later using the fluorescent vital dye fluorescent dye 5-Carboxyfluorescein Diacetate (CFDA), yielding concentration-effect curves from each experiment. Technical (day-to-day) variability of the assay, assessed from three independent experiments, was low: p-values for the differences of results were 0.74 and 0.64 for the comparison of day 2 vs. day 1 and day 3 vs. day 1, respectively. The statistical analysis used to compare the entire concentration-effect curves has revealed significant differences in levels of cytotoxicity induced by the toxic and relatively nontoxic chemical counterparts, demonstrating that donor-specific ECFCs can clearly differentiate between these two groups of chemicals. Partitioning of the total variance in the nested design assessed the contributions of between-clone and between-donor variability for different levels of cytotoxicity. Individual ECFC clones demonstrated highly reproducible responses to the chemicals. The most toxic chemical was TBT, followed by NaAsO2, CdCl2, and Menadione. Nontoxic counterparts exhibited low cytotoxicity at the higher end of concentration ranges tested. Low variability was observed between ECFC clones obtained from the same donor or different donors for CdCl2, NaAsO2, and TBT, but for menadione, the between-donor variability was much greater than the between-clone variability. The low between-clone variability indicates that an ECFC clone may represent an individual donor in cell-based assays, although this finding must be confirmed using a larger number of donors. Such confirmation would demonstrate that an in vitro ECFC-based testing platform can be used to characterize the inter-individual variability of neonatal ECFCs exposed to drugs and/or environmental toxicants.

Characterizing Serpinb2 as a Modulator of TCDD-Induced Suppression of the B Cell

2,3,7,8-Tetrachlordibenzo- p-dioxin (TCDD) is an environmental pollutant that can cause various toxic effects, including chloracne, metabolic syndrome, and immune suppression. Most of the toxicity associated with TCDD is mediated through activation of the aryl hydrocarbon receptor (AHR). Recent research has suggested the presence of a wide-range of interindividual variability in TCDD-mediated suppression of the Immunoglobulin-M (IgM) response across the human population. In an attempt to identify putative modifiers of AHR-mediated immunosuppression beyond the AHR, B cells were isolated from a panel of genetically diverse mouse strain to scan for modulators that drive interstrain differences in TCDD-mediated suppression of the IgM response. Results implicated a region of mouse Chromosome 1 near a gene encoding serine peptidase inhibitor, clade B, member 2 ( Serpinb2) whose human ortholog is plasminogen activator inhibitor 2 (PAI2). Further downstream analyses indicated that Serpinb2 is dysregulated by TCDD and, furthermore, that B cells from Serpinb2 ^-/- mice are significantly more sensitive to TCDD-mediated suppression as compared to littermate controls. This study suggests a protective role of Serpinb2 within TCDD-mediated immunosuppression and, furthermore, a novel function of Serpinb2-related activity in the IgM response.