Evaluation of 5-day In Vivo Rat Liver and Kidney With High-throughput Transcriptomics for Estimating Benchmark Doses of Apical Outcomes

A 5-day in vivo rat model was evaluated as an approach to estimate chemical exposures that may pose minimal risk by comparing benchmark dose (BMD) values for transcriptional changes in the liver and kidney to BMD values for toxicological endpoints from traditional toxicity studies. Eighteen chemicals, most having been tested by the National Toxicology Program in 2-year bioassays, were evaluated. Some of these chemicals are potent hepatotoxicants (eg, DE71, PFOA, and furan) in rodents, some exhibit toxicity but have minimal hepatic effects (eg, acrylamide and α,β-thujone), and some exhibit little overt toxicity (eg, ginseng and milk thistle extract) based on traditional toxicological evaluations. Male Sprague Dawley rats were exposed once daily for 5 consecutive days by oral gavage to 8-10 dose levels for each chemical. Liver and kidney were collected 24 h after the final exposure and total RNA was assayed using high-throughput transcriptomics (HTT) with the rat S1500+ platform. HTT data were analyzed using BMD Express 2 to determine transcriptional gene set BMD values. BMDS was used to determine BMD values for histopathological effects from chronic or subchronic toxicity studies. For many of the chemicals, the lowest transcriptional BMDs from the 5-day assays were within a factor of 5 of the lowest histopathological BMDs from the toxicity studies. These data suggest that using HTT in a 5-day in vivo rat model provides reasonable estimates of BMD values for traditional apical endpoints. This approach may be useful to prioritize chemicals for further testing while providing actionable data in a timely and cost-effective manner.

Population genetic diversity in zebrafish lines

Toxicological and pharmacological researchers have seized upon the many benefits of zebrafish, including the short generation time, well-characterized development, and early maturation as clear embryos. A major difference from many model organisms is that standard husbandry practices in zebrafish are designed to maintain population diversity. While this diversity is attractive for translational applications in human and ecological health, it raises critical questions on how interindividual genetic variation might contribute to chemical exposure or disease susceptibility differences. Findings from pooled samples of zebrafish support this supposition of diversity yet cannot directly measure allele frequencies for reference versus alternate alleles. Using the Tanguay lab Tropical 5D zebrafish line (T5D), we performed whole genome sequencing on a large group (n = 276) of individual zebrafish embryos. Paired-end reads were collected on an Illumina 3000HT, then aligned to the most recent zebrafish reference genome (GRCz10). These data were used to compare observed population genetic variation across species (humans, mice, zebrafish), then across lines within zebrafish. We found more single nucleotide polymorphisms (SNPs) in T5D than have been reported in SNP databases for any of the WIK, TU, TL, or AB lines. We theorize that some subset of the novel SNPs may be shared with other zebrafish lines but have not been identified in other studies due to the limitations of capturing population diversity in pooled sequencing strategies. We establish T5D as a model that is representative of diversity levels within laboratory zebrafish lines and demonstrate that experimental design and analysis can exert major effects when characterizing genetic diversity in heterogeneous populations.

Systematic determination of the relationship between nanoparticle core diameter and toxicity for a series of structurally analogous gold nanoparticles in zebrafish

Predictive models for the impact of nanomaterials on biological systems remain elusive. Although there is agreement that physicochemical properties (particle diameter, shape, surface chemistry, and core material) influence toxicity, there are limited and often contradictory, data relating structure to toxicity, even for core diameter. Given the importance of size in determining nanoscale properties, we aimed to address this data gap by examining the biological effects of a defined series of gold nanoparticles (AuNPs) on zebrafish embryos. Five AuNPs samples with narrowly spaced core diameters (0.8-5.8 nm) were synthesized and functionalized with positively charged N,N,N-trimethylammonium ethanethiol (TMAT) ligands. We assessed the bioactivity of these NPs in a high-throughput developmental zebrafish assay at eight concentrations (0.5-50 µg/mL) and observed core diameter-dependent bioactivity. The smaller diameter AuNPs were the most toxic when expressing exposures based on an equal mass. However, when expressing exposures based on total surface area, toxicity was independent of the core diameter. When holding the number of nanoparticles per volume constant (at 6.71 × 10¹³/mL) in the exposure medium across AuNPs diameters, only the 5.8 nm AuNPs exhibited toxic effects. Under these exposure conditions, the uptake of AuNPs in zebrafish was only weakly associated with core diameter, suggesting that differential uptake of TMAT-AuNPs was not responsible for toxicity associated with the 5.8 nm core diameter. Our results indicate that larger NPs may be the most toxic on a per particle basis and highlight the importance of using particle number and surface area, in addition to mass, when evaluating the size-dependent bioactivity of NPs.

Elucidating Gene-by-Environment Interactions Associated with Differential Susceptibility to Chemical Exposure

BACKGROUND

Modern societies are exposed to vast numbers of potentially hazardous chemicals. Despite demonstrated linkages between chemical exposure and severe health effects, there are limited, often conflicting, data on how adverse health effects of exposure differ across individuals.

OBJECTIVES

We tested the hypothesis that population variability in response to certain chemicals could elucidate a role for gene-environment interactions (GxE) in differential susceptibility.

METHODS

High-throughput screening (HTS) data on thousands of chemicals in genetically heterogeneous zebrafish were leveraged to identify a candidate chemical (Abamectin) with response patterns indicative of population susceptibility differences. We tested the prediction by generating genome-wide sequence data for 276 individual zebrafish displaying susceptible (Affected) vs. resistant (Unaffected) phenotypes following identical chemical exposure.

RESULTS

We found GxE associated with differential susceptibility in the sox7 promoter region and then confirmed gene expression differences between phenotypic response classes.

CONCLUSIONS

The results for Abamectin in zebrafish demonstrate that GxE associated with naturally occurring, population genetic variation play a significant role in mediating individual response to chemical exposure. https://doi.org/10.1289/EHP2662.

Transcriptional pathways linked to fetal and maternal hepatic dysfunction caused by gestational exposure to perfluorooctanoic acid (PFOA) or hexafluoropropylene oxide-dimer acid (HFPO-DA or GenX) in CD-1 mice

Per- and polyfluoroalkyl substances (PFAS) comprise a diverse class of chemicals used in industrial processes, consumer products, and fire-fighting foams which have become environmental pollutants of concern due to their persistence, ubiquity, and associations with adverse human health outcomes, including in pregnant persons and their offspring. Multiple PFAS are associated with adverse liver outcomes in adult humans and toxicological models, but effects on the developing liver are not fully described. Here we performed transcriptomic analyses in the mouse to investigate the molecular mechanisms of hepatic toxicity in the dam and its fetus after exposure to two different PFAS, perfluorooctanoic acid (PFOA) and its replacement, hexafluoropropylene oxide-dimer acid (HFPO-DA, known as GenX). Pregnant CD-1 mice were exposed via oral gavage from embryonic day (E) 1.5-17.5 to PFOA (0, 1, or 5 mg/kg-d) or GenX (0, 2, or 10 mg/kg-d). Maternal and fetal liver RNA was isolated (N = 5 per dose/group) and the transcriptome analyzed by Affymetrix Array. Differentially expressed genes (DEG) and differentially enriched pathways (DEP) were obtained. DEG patterns were similar in maternal liver for 5 mg/kg PFOA, 2 mg/kg GenX, and 10 mg/kg GenX (R2: 0.46-0.66). DEG patterns were similar across all 4 dose groups in fetal liver (R2: 0.59-0.81). There were more DEGs in fetal liver compared to maternal liver at the low doses for both PFOA (fetal = 69, maternal = 8) and GenX (fetal = 154, maternal = 93). Upregulated DEPs identified across all groups included Fatty Acid Metabolism, Peroxisome, Oxidative Phosphorylation, Adipogenesis, and Bile Acid Metabolism. Transcriptome-phenotype correlation analyses demonstrated > 1000 maternal liver DEGs were significantly correlated with maternal relative liver weight (R2 >0.92). These findings show shared biological pathways of liver toxicity for PFOA and GenX in maternal and fetal livers in CD-1 mice. The limited overlap in specific DEGs between the dam and fetus suggests the developing liver responds differently than the adult liver to these chemical stressors. This work helps define mechanisms of hepatic toxicity of two structurally unique PFAS and may help predict latent consequences of developmental exposure.

Development of a Zebrafish S1500+ Sentinel Gene Set for High-Throughput Transcriptomics

Sentinel gene sets have been developed with the purpose of maximizing the information from targeted transcriptomic platforms. We recently described the development of an S1500+ sentinel gene set, which was built for the human transcriptome, utilizing a data- and knowledge-driven hybrid approach to select a small subset of genes that optimally capture transcriptional diversity, correlation with other genes based on large-scale expression profiling, and known pathway annotation within the human genome. While this detailed bioinformatics approach for gene selection can in principle be applied to other species, the reliability of the resulting gene set depends on availability of a large body of transcriptomics data. For the model organism zebrafish, we aimed to create a similar sentinel gene set (Zf S1500+ gene set); however, there is insufficient standardized expression data in the public domain to train the gene correlation model. Therefore, our strategy was to use human-zebrafish ortholog mapping of the human S1500+ genes and nominations from experts in the zebrafish scientific community. In this study, we present the bioinformatics curation and refinement process to produce the final Zf S1500+ gene set, explore whole transcriptome extrapolation using this gene set, and assess pathway-level inference. This gene set will add value to targeted high-throughput transcriptomics in zebrafish for toxicogenomic screening and other research domains.

Utility of Extrapolating Human S1500+ Genes to the Whole Transcriptome: Tunicamycin Case Study

The TempO-Seq S1500+ platform(s), now available for human, mouse, rat, and zebrafish, measures a discrete number of genes that are representative of biological and pathway co-regulation across the entire genome in a given species. While measurement of these genes alone provides a direct assessment of gene expression activity, extrapolating expression values to the whole transcriptome (~26 000 genes in humans) can estimate measurements of non-measured genes of interest and increases the power of pathway analysis algorithms by using a larger background gene expression space. Here, we use data from primary hepatocytes of 54 donors that were treated with the endoplasmic reticulum (ER) stress inducer tunicamycin and then measured on the human S1500+ platform containing ~3000 representative genes. Measurements for the S1500+ genes were then used to extrapolate expression values for the remaining human transcriptome. As a case study of the improved downstream analysis achieved by extrapolation, the “measured only” and “whole transcriptome” (measured + extrapolated) gene sets were compared. Extrapolation increased the number of significant genes by 49%, bringing to the forefront many that are known to be associated with tunicamycin exposure. The extrapolation procedure also correctly identified established tunicamycin-related functional pathways reflected by coordinated changes in interrelated genes while maintaining the sample variability observed from the “measured only” genes. Extrapolation improved the gene- and pathway-level biological interpretations for a variety of downstream applications, including differential expression analysis, gene set enrichment pathway analysis, DAVID keyword analysis, Ingenuity Pathway Analysis, and NextBio correlated compound analysis. The extrapolated data highlight the role of metabolism/metabolic pathways, the ER, immune response, and the unfolded protein response, each of which are key activities associated with tunicamycin exposure that were unrepresented or underrepresented in one or more of the analyses of the original “measured only” dataset. Furthermore, the inclusion of the extrapolated genes raised “tunicamycin” from third to first upstream regulator in Ingenuity Pathway Analysis and from sixth to second most correlated compound in NextBio analysis. Therefore, our case study suggests an approach to extend and enhance data from the S1500+ platform for improved insight into biological mechanisms and functional outcomes of diseases, drugs, and other perturbations.

In vivo assessment of respiratory burst inhibition by xenobiotic exposure using larval zebrafish

Currently, assessment of the potential immunotoxicity of a given agent involves a tiered approach for hazard identification and mechanistic studies, including observational studies, evaluation of immune function, and measurement of susceptibility to infectious and neoplastic diseases. These studies generally use costly low-throughput mammalian models. Zebrafish, however, offer an excellent alternative due to their rapid development, ease of maintenance, and homology to mammalian immune system function and development. Larval zebrafish also are a convenient model to study the innate immune system with no interference from the adaptive immune system. In this study, a respiratory burst assay (RBA) was utilized to measure reactive oxygen species (ROS) production after developmental xenobiotic exposure. Embryos were exposed to non-teratogenic doses of chemicals and at 96 h post-fertilization, the ability to produce ROS was measured. Using the RBA, 12 compounds with varying immune-suppressive properties were screened. Seven compounds neither suppressed nor enhanced the respiratory burst; five reproducibly suppressed global ROS production, but with varying potencies: benzo[a]pyrene, 17β-estradiol, lead acetate, methoxychlor, and phenanthrene. These five compounds have all previously been reported as immunosuppressive in mammalian innate immunity assays. To evaluate whether the suppression of ROS by these compounds was a result of decreased immune cell numbers, flow cytometry with transgenic zebrafish larvae was used to count the numbers of neutrophils and macrophages after chemical exposure. With this assay, benzo[a]pyrene was found to be the only chemical that induced a change in the number of immune cells by increasing macrophage but not neutrophil numbers. Taken together, this work demonstrates the utility of zebrafish larvae as a vertebrate model for identifying compounds that impact innate immune function at non-teratogenic levels and validates measuring ROS production and phagocyte numbers as metrics for monitoring how xenobiotic exposure alters the innate immune system.

Impact of Aligner, Normalization Method, and Sequencing Depth on TempO-seq Accuracy

High-throughput transcriptomics has advanced through the introduction of TempO-seq, a targeted alternative to traditional RNA-seq. TempO-seq platforms use 50 nucleotide probes, each specifically designed to target a known transcript, thus allowing for reduced sequencing depth per sample compared with RNA-seq without compromising the accuracy of results. Thus far, studies using the TempO-seq method have relied on existing tools for processing the resulting short read data. However, these tools were originally designed for other data types. While they have been used for processing of early TempO-seq data, they have not been systematically assessed for accuracy or compared to determine an optimal framework for processing and analyzing TempO-seq data. In this work, we re-analyze several publicly available TempO-seq data sets covering a range of experimental designs and use corresponding RNA-seq data sets as a gold standard to rigorously assess accuracy at multiple levels. We compare 6 aligners and 5 normalization methods across various accuracy and performance metrics. Our results demonstrate the overall robust accuracy of the TempO-seq platform, independent of data processing methods. Complex aligners and advanced normalization methods do not appear to have any general advantage over simpler methods when it comes to analyzing TempO-seq data. The reduced complexity of the sequencing space, and the fact that TempO-seq probes are all equal length, appears to reduce the need for elaborate bioinformatic or statistical methods used to address these factors in RNA-seq data.

Quantifying liver-toxic responses from dose-dependent chemical exposures using a rat genome-scale metabolic model

Because the liver plays a vital role in the clearance of exogenous chemical compounds, it is susceptible to chemical-induced toxicity. Animal-based testing is routinely used to assess the hepatotoxic potential of chemicals. Although large-scale high-throughput sequencing data can indicate the genes affected by chemical exposures, we need system-level approaches to interpret these changes. To this end, we developed an updated rat genome-scale metabolic model to integrate large-scale transcriptomics data and utilized a chemical structure similarity-based ToxProfiler tool to identify chemicals that bind to specific toxicity targets to understand the mechanisms of toxicity. We used high-throughput transcriptomics data from a 5-day in vivo study where rats were exposed to different non-toxic and hepatotoxic chemicals at increasing concentrations and investigated how liver metabolism was differentially altered between the non-toxic and hepatotoxic chemical exposures. Our analysis indicated that the genes identified via toxicity target analysis and those mapped to the metabolic model showed a distinct gene expression pattern, with the majority showing upregulation for hepatotoxicants compared with non-toxic chemicals. Similarly, when we mapped the metabolic genes at the pathway level, we identified several pathways in carbohydrate, amino acid, and lipid metabolism that were significantly upregulated for hepatotoxic chemicals. Furthermore, using our system-level integration of gene expression data with the rat metabolic model, we could differentiate metabolites in these pathways that were systematically elevated or suppressed due to hepatotoxic versus non-toxic chemicals. Thus, using our combined approach, we were able to identify a set of potential gene signatures that clearly differentiated liver toxic responses from non-toxic chemicals, which helped us identify potential metabolic pathways and metabolites that are systematically associated with the toxicant exposure.

Characterization and optimization of variability in a human colonic epithelium culture model

Animal models have historically been poor preclinical predictors of gastrointestinal (GI) directed therapeutic efficacy and drug-induced GI toxicity. Human stem and primary cell-derived culture systems are a major focus of efforts to create biologically relevant models that enhance preclinical predictive value of intestinal efficacy and toxicity. The inherent variability in stem cell-based cultures makes development of useful models a challenge; the stochastic nature of stem cell differentiation interferes with the ability to build and validate reproducible assays that query drug responses and pharmacokinetics. In this study, we aimed to characterize and reduce sources of variability in a complex stem cell-derived intestinal epithelium model, termed RepliGut® Planar, across cells from multiple human donors, cell lots, and passage numbers. Assessment criteria included barrier for­mation and integrity, gene expression, and cytokine responses. Gene expression and culture metric analyses revealed that controlling cell passage number reduces variability and maximizes physi­ological relevance of the model. In a case study where passage number was optimized, distinct cytokine responses were observed among four human donors, indicating that biological variability can be detected in cell cultures originating from diverse human sources. These findings highlight key considerations for designing assays that can be applied to additional primary cell-derived systems, as well as establish utility of the RepliGut® Planar platform for robust development of human-pre­dictive drug-response assays.

Correction: Pannala et al. High-Throughput Transcriptomics Differentiates Toxic versus Non-Toxic Chemical Exposures Using a Rat Liver Model. Int. J. Mol. Sci. 2023, 24, 17425

In the original publication [...].

Characterization and optimization of variability in a human colonic epithelium culture model

Animal models have historically been poor preclinical predictors of gastrointestinal (GI) directed therapeutic efficacy and drug-induced GI toxicity. Human stem and primary cell-derived culture systems are a major focus of efforts to create biologically relevant models that enhance preclinical predictive value of intestinal efficacy and toxicity. The inherent variability in stem-cell-based complex cultures makes development of useful models a challenge; the stochastic nature of stem-cell differentiation interferes with the ability to build and validate robust, reproducible assays that query drug responses and pharmacokinetics. In this study, we aimed to characterize and reduce potential sources of variability in a complex stem cell-derived intestinal epithelium model, termed RepliGut® Planar, across cells from multiple human donors, cell lots, and passage numbers. Assessment criteria included barrier formation and integrity, gene expression, and cytokine responses. Gene expression and culture metric analyses revealed that controlling for stem/progenitor-cell passage number reduces variability and maximizes physiological relevance of the model. After optimizing passage number, donor-specific differences in cytokine responses were observed in a case study, suggesting biologic variability is observable in cell cultures derived from multiple human sources. Our findings highlight key considerations for designing assays that can be applied to additional primary-cell derived systems, as well as establish utility of the RepliGut® Planar platform for robust development of human-predictive drug-response assays.

Unraveling Human Hepatocellular Responses to PFAS and Aqueous Film-Forming Foams (AFFFs) for Molecular Hazard Prioritization and In Vivo Translation

Aqueous film-forming foams (AFFFs) are complex product mixtures that often contain per- and polyfluorinated alkyl substances (PFAS) to enhance fire suppression and protect firefighters. However, PFAS have been associated with a range of adverse health effects (e.g., liver and thyroid disease and cancer), and innovative approach methods to better understand their toxicity potential and identify safer alternatives are needed. In this study, we investigated a set of 30 substances (e.g., AFFF, PFAS, and clinical drugs) using differentiated cultures of human hepatocytes (HepaRG, 2D), high-throughput transcriptomics, deep learning of cell morphology images, and liver enzyme leakage assays with benchmark dose analysis to (1) predict the potency ranges for human liver injury, (2) delineate gene- and pathway-level transcriptomic points-of-departure for molecular hazard characterization and prioritization, (3) characterize human hepatocellular response similarities to inform regulatory read-across efforts, and (4) introduce an innovative approach to translate mechanistic hepatocellular response data to predict the potency ranges for PFAS-induced hepatomegaly in vivo. Collectively, these data fill important mechanistic knowledge gaps with PFAS/AFFF and represent a scalable platform to address the thousands of PFAS in commerce for greener chemistries and next-generation risk assessments.

Exposure to PFAS chemicals induces sex-dependent alterations in key rate-limiting steps of lipid metabolism in liver steatosis

Toxicants with the potential to bioaccumulate in humans and animals have long been a cause for concern, particularly due to their association with multiple diseases and organ injuries. Per- and polyfluoro alkyl substances (PFAS) and polycyclic aromatic hydrocarbons (PAH) are two such classes of chemicals that bioaccumulate and have been associated with steatosis in the liver. Although PFAS and PAH are classified as chemicals of concern, their molecular mechanisms of toxicity remain to be explored in detail. In this study, we aimed to identify potential mechanisms by which an acute exposure to PFAS and PAH chemicals can induce lipid accumulation and whether the responses depend on chemical class, dose, and sex. To this end, we analyzed mechanisms beginning with the binding of the chemical to a molecular initiating event (MIE) and the consequent transcriptomic alterations. We collated potential MIEs using predictions from our previously developed ToxProfiler tool and from published steatosis adverse outcome pathways. Most of the MIEs are transcription factors, and we collected their target genes by mining the TRRUST database. To analyze the effects of PFAS and PAH on the steatosis mechanisms, we performed a computational MIE-target gene analysis on high-throughput transcriptomic measurements of liver tissue from male and female rats exposed to either a PFAS or PAH. The results showed peroxisome proliferator-activated receptor (PPAR)-α targets to be the most dysregulated, with most of the genes being upregulated. Furthermore, PFAS exposure disrupted several lipid metabolism genes, including upregulation of fatty acid oxidation genes (Acadm, Acox1, Cpt2, Cyp4a1-3) and downregulation of lipid transport genes (Apoa1, Apoa5, Pltp). We also identified multiple genes with sex-specific behavior. Notably, the rate-limiting genes of gluconeogenesis (Pck1) and bile acid synthesis (Cyp7a1) were specifically downregulated in male rats compared to female rats, while the rate-limiting gene of lipid synthesis (Scd) showed a PFAS-specific upregulation. The results suggest that the PPAR signaling pathway plays a major role in PFAS-induced lipid accumulation in rats. Together, these results show that PFAS exposure induces a sex-specific multi-factorial mechanism involving rate-limiting genes of gluconeogenesis and bile acid synthesis that could lead to activation of an adverse outcome pathway for steatosis.