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 formation and integrity, gene expression, and cytokine responses. Gene expression and culture metric analyses revealed that controlling cell passage number reduces variability and maximizes physiological 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-predictive drug-response assays.

Animal models are frequently used as tools for studying gastrointestinal (GI) disease, but they poorly replicate the complexities of the human gut limiting the clinical translation of new therapeutics in development. Human stem cell derived models can better recapitulate human GI physiology, but the inherent dynamic nature of stem cells introduces variability in culture performance. We identified sources of variability in the primary stem-cell derived RepliGut® Planar model to develop robust and reliable assays that can improve preclinical therapeutic development for GI disease. Analysis of barrier formation, gene expression, and cytokine responses demonstrated that controlling cell passage number reduces variability and maximizes physiological relevance of the model. These findings highlight key assay design considerations that can be applied to additional primary-cell derived systems. Availability of reliable and physiologically relevant cell-based models can reduce animal testing, improve research accuracy, and make new treatments more relevant and effective for patients.

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.

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.

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.

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.

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.

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.

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.

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 × 1013/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.

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.