Benchmark Concentration Research Articles

Peritoneal dialysis-related complication: Can diagnostic accuracy be achieved at low iodinated contrast media dose using dual-layer spectral detector CT?

Computed tomography peritoneography (CTp) is pivotal for evaluating peritoneal dialysis (PD)-related complications, yet it comes with drawbacks, specifically exposure to iodinated contrast media (ICM). This study aimed to explore the feasibility of reducing ICM dosage utilizing spectral detector CT (SDCT). 35 rabbits were strategically divided into 7 groups (A - G) according to the ICM concentration ratio in the injection protocol, with respective doses of 10, 15, 20, 25, 30, 40, and 50 mL/2L. The CTp injection protocol involved a 300-mL mixture of non-ionic ICM omnipaque (350 mgI/mL) and peritoneal dialysate (1.5% lactate, 2 L), followed by scans using dual-layer SDCT. Virtual monoenergetic images (VMIs) at 4 distinct energy levels (40 - 70 keV, in 10-keV steps), iodine maps (IMs), and effective atomic number (Zeff) maps were subsequently reconstructed. Both quantitative and qualitative image assessments were conducted, and the parameters from these analyses were compared across images from groups A-G and traditional 50 mL/2L 120-kVp images. In post-determination of the optimal concentration and reconstructions, we illustrated their applications in patients with suspected PD-related complications. The quantitative image quality (IQ) of 15 mL/2L VMIs at 40 keV surpassed that of the 50 mL/2L 120-kVp images (p < 0.05). Furthermore, the diagnostic performance utilizing 15 mL/2L VMIs40 keV, when combined with IMs and Zeff maps, was found to be optimal. The employment of SDCT in CTp allows for a substantial reduction in the ICM dose by 70%, compared to the benchmark concentration of 50 mL/2L, without compromising diagnostic precision.

Histone Methylation-Mediated Reproductive Toxicity to Consumer Product Chemicals in Caenorhabditis elegans: An Epigenetic Adverse Outcome Pathway (AOP).

The significance of histone methylation in epigenetic inheritance underscores its relevance to disease and the chronic effects of environmental chemicals. However, limited evidence of the causal relationships between chemically induced epigenetic changes and organismal-level effects hinders the application of epigenetic markers in ecotoxicological assessments. This study explored the contribution of repressive histone marks to reproductive toxicity induced by chemicals in consumer products in Caenorhabditis elegans, applying the adverse outcome pathway (AOP) framework. Triclosan (TCS) and tetrabromobisphenol A (TBBPA) exposures caused reproductive toxicity and altered histone methyltransferase (HMT) and histone demethylase (HDM) activities, increasing the level of trimethylation of H3K9 and H3K27. Notably, treatment with an H3K27-specific HMT inhibitor alleviated reproductive defects and the transcriptional response of genes related to vitellogenin, xenobiotic metabolism, and oxidative stress. Comparison of points of departure (PODs) based on calculated benchmark concentrations (BMCs) revealed the sensitivity of histone-modifying enzyme activities to these chemicals. Our findings suggest that the 'disturbance of HMT and HDM' can serve as the molecular initiating event (MIE) leading to reproductive toxicity in the epigenetic AOP for TCS and TBBPA. The study extended the biological applicability of these enzymes by identifying model species with analogous protein sequences and functions. This combined approach enhances the essentiality, empirical support, and taxonomic domain of applicability (tDOA), which are crucial considerations for ecotoxicological AOPs. Given the widespread use and environmental distribution of chemicals in consumer products, this study proposes histone-modifying enzyme activity as an effective screening tool for reproductive toxicants and emphasizes the integration of epigenetic mechanisms into a prospective ERA.

Mode of action exploration for prostate epithelial cell injury caused by bisphenol A

Bisphenol A (BPA) is a typical food chemical contaminant with various detrimental effects, especially on reproductive system. Male prostate damage is also one of its major adverse health effects, of which mode of action (MOA) remains unclear. This study aims to explore the MOA for prostate toxicity of BPA using human normal prostate epithelial cell RWPE-1 for 28-day human-relevant-level exposure. A physiological based pharmacokinetic model was used to determine the concentration of BPA based on the actual oral exposure in China. The possible key events were identified by high-throughput transcriptome sequencing and validated by qPCR, Western blot and cell cycle assay, and the benchmark concentration analysis were conducted. The enriched KEGG pathways include the endocytosis, cell cycle, cellular senescence, MAPK and TNF signaling pathways. With increasing BPA concentrations, the increased mRNA and/or protein expressions of MAPKAPK2, c-JUN and c-fos in the MAPK signaling pathway, the increased mRNA expressions of CCND1 and CDKN1A, the decreased mRNA expression of CDC25C, the increased proportion of G0/G1 phase and S phase, as well as the decreased proportion of G2/M phase, were observed. The lowest value of benchmark concentration lower confidence limit (BMCL) was retrieved from G2/M phase ratio, with 110.580 and 175.862 nM for BMCL5 and BMCL10, respectively, much higher than the male gonad maximum concentration of 0.019 nM of BPA at the current exposure level of adult Chinese males. In conclusion, the MOA of BPA induced male prostatic toxicity at human-relevant levels may include: key event (KE)1-MAPK signaling pathway activation, KE2-disorder of cell cycle regulatory gene expression (increased expression of CCND1 and CDKN1A, decreased expression of CDC25C), and KE3-disturbance of cell cycle (increased proportion of G0/G1 and S phases, decreased proportion of G2/M phases). However, more studies are needed to validate and complete the MOA.

Comparison on the mechanism and potency of hepatotoxicity among hemp extract and its four major constituent cannabinoids

Cannabidiol (CBD) has been reported to induce hepatotoxicity in clinical trials and research studies; however, little is known about the safety of other nonintoxicating cannabinoids. New approach methodologies (NAMs) based on bioinformatic analysis of high-throughput transcriptomic data are gaining increasing importance in risk assessment and regulatory decision-making of data-poor chemicals. In the current study, we conducted a concentration response transcriptomic analysis of hemp extract and its four major constituent cannabinoids [CBD, cannabichromene (CBC), cannabigerol (CBG), and cannabinol (CBN)] in hepatocytes derived from human induced pluripotent stem cells (iPSCs). Each compound impacted a distinctive combination of biological functions and pathways. However, all the cannabinoids impaired liver metabolism and caused oxidative stress in the cells. Benchmark concentration (BMC) analysis showed potencies in transcriptional activity of the cannabinoids were in the order of CBN > CBD > CBC > CBG, consistent with the order of their cytotoxicity IC50 values. Patterns of transcriptomic changes induced by hemp extract and its median overall BMC were very similar to CBD but differed significantly from other cannabinoids, suggesting that potential adverse effects of hemp extract were largely due to its major constituent CBD. Lastly, transcriptomic point-of-departure (tPoD) values were determined for each of the compounds, with the value for CBD (0.106 µM) being concordant with a previously reported one derived from apical endpoints of clinical and animal studies. Taken together, the current study demonstrates the potential utility of transcriptomic BMC analysis as a NAM for hazard assessment of data-poor chemicals, improves our understanding of the possible health effects of hemp extract and its constituent cannabinoids, and provides important tPoD data that could contribute to inform human safety assessment of these cannabinoid compounds.

Using Zebrafish to Screen Developmental Toxicity of Per- and Polyfluoroalkyl Substances (PFAS).

Per- and polyfluoroalkyl substances (PFAS) are found in many consumer and industrial products. While some PFAS, notably perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), are developmentally toxic in mammals, the vast majority of PFAS have not been evaluated for developmental toxicity potential. A concentration-response study of 182 unique PFAS chemicals using the zebrafish medium-throughput, developmental vertebrate toxicity assay was conducted to investigate chemical structural identifiers for toxicity. Embryos were exposed to each PFAS compound (≤100 μM) beginning on the day of fertilization. At 6 days post-fertilization (dpf), two independent observers graded developmental landmarks for each larva (e.g., mortality, hatching, swim bladder inflation, edema, abnormal spine/tail, or craniofacial structure). Thirty percent of the PFAS were developmentally toxic, but there was no enrichment of any OECD structural category. PFOS was developmentally toxic (benchmark concentration [BMC] = 7.48 μM); however, other chemicals were more potent: perfluorooctanesulfonamide (PFOSA), N-methylperfluorooctane sulfonamide (N-MeFOSA), ((perfluorooctyl)ethyl)phosphonic acid, perfluoro-3,6,9-trioxatridecanoic acid, and perfluorohexane sulfonamide. The developmental toxicity profile for these more potent PFAS is largely unexplored in mammals and other species. Based on these zebrafish developmental toxicity results, additional screening may be warranted to understand the toxicity profile of these chemicals in other species.

Hepatotoxicity of silver nanoparticles: Benchmark concentration modeling of an in vitro transcriptomics study in human iPSC-derived hepatocytes

Despite two decades of research on silver nanoparticle (AgNP) toxicity, a safe threshold for exposure has not yet been established, albeit being critically needed for risk assessment and regulatory decision-making. Traditionally, a point-of-departure (PoD) value is derived from dose response of apical endpoints in animal studies using either the no-observed-adverse-effect level (NOAEL) approach, or benchmark dose (BMD) modeling. To develop new approach methodologies (NAMs) to inform human risk assessment of AgNPs, we conducted a concentration response modeling of the transcriptomic changes in hepatocytes derived from human induced pluripotent stem cells (iPSCs) after being exposed to a wide range concentration (0.01–25 μg/ml) of AgNPs for 24 h. A plausible transcriptomic PoD of 0.21 μg/ml was derived for a pathway related to the mode-of-action (MOA) of AgNPs, and a more conservative PoD of 0.10 μg/ml for a gene ontology (GO) term not apparently associated with the MOA of AgNPs. A reference dose (RfD) could be calculated from either of the PoDs as a safe threshold for AgNP exposure. The current study illustrates the usefulness of in vitro transcriptomic concentration response study using human cells as a NAM for toxicity study of chemicals that lack adequate toxicity data to inform human risk assessment.

Modeling HepaRG metabolome responses to pyrrolizidine alkaloid exposure for insight into points of departure and modes of action

The new challenges in toxicology demand novel and innovative in vitro approaches for deriving points of departure (PODs) and determining the mode of action (MOA) of chemicals. Therefore, the aim of this original study was to couple in vitro studies with untargeted metabolomics to model the concentration-response of extra- and intracellular metabolome data on human HepaRG cells treated for 48 h with three pyrrolizidine alkaloids (PAs): heliotrine, retrorsine and lasiocarpine. Modeling revealed that the three PAs induced various monotonic and, importantly, biphasic curves of metabolite content. Based on unannotated metabolites, the endometabolome was more sensitive than the exometabolome in terms of metabolomic effects, and benchmark concentrations (BMCs) confirmed that lasiocarpine was the most hepatotoxic PA. Regarding its MOA, impairment of lipid metabolism was highlighted at a very low BMC (first quartile, 0.003 µM). Moreover, results confirmed that lasiocarpine targets bile acids, as well as amino acid and steroid metabolisms. Analysis of the endometabolome, based on coupling concentration-response and PODs, gave encouraging results for ranking toxins according to their hepatotoxic effects. Therefore, this novel approach is a promising tool for next-generation risk assessment, readily applicable to a broad range of compounds and toxic endpoints.

Water Quality and the First-Flush Effect in Roof-Based Rainwater Harvesting, Part I: Water Quality and Soil Accumulation

Rainfall runoff may be captured and stored for later use, but the quality of this water can be detrimental in some uses without the use of appropriately designed first-flush diverters. The rainfall runoff water quality was measured on nineteen new small-scale and two aged commercial roofs located near high traffic highways. Roof coverings included asphalt shingles, sheet metal, clay tiles, and tar and gravel. Runoff samples were evaluated for polycyclic aromatic hydrocarbons (PAHs), phosphorus flame retardants (PFRs), and pyrethroid insecticides. Eighteen small-scale roofs were subjected to a range of simulated rainfall events, while natural runoff was sampled on the commercial roofs and one small-scale roof. Runoff was analyzed for pH, conductivity, turbidity, total suspended solids, boron, iron, copper, zinc, manganese, sodium adsorption ratio, nitrate-nitrogen, seventeen PAHs, tris(2-chloroethyl) phosphate, tris(1,3-dichloro-2-propyl)phosphate, bifenthrin, cypermethrin, and lambda-cyhalothrin. Samples from four natural storm events were also analyzed for total coliforms and Escherichia coli. In addition, soils below seventeen existing gutter downspouts were sampled to determine long-term pollutant accumulation. Atmospheric deposition was the main contributor of pollutants in the roof runoff. A majority of samples fell within the U.S. EPA guidelines for non-potable urban and agricultural water reuse. Trace levels of PAHs, PFRs, and insecticides were detected, but all detections were three orders of magnitude below the USGS health-based screening level benchmark concentrations. Results indicate that diverting the first flush, based on turbidity, total suspended solids, or conductivity, can improve the overall water quality and reduce the concentrations of PAHs in harvested rainwater. Downspout soil sampling showed potential for the long-term accumulation of PAHs at concentrations exceeding the minimum human-health risk-based screening levels at these high runoff-loading locations.

Probabilistic Reference and 10% Effect Concentrations for Characterizing Inhalation Non-cancer and Developmental/Reproductive Effects for 2,160 Substances.

Chemicals assessment and management frameworks rely on regulatory toxicity values, which are based on points of departure (POD) identified following rigorous dose-response assessments. Yet, regulatory PODs and toxicity values for inhalation exposure (i.e., reference concentrations [RfCs]) are available for only ∼200 chemicals. To address this gap, we applied a workflow to determine surrogate inhalation route PODs and corresponding toxicity values, where regulatory assessments are lacking. We curated and selected inhalation in vivo data from the U.S. EPA's ToxValDB and adjusted reported effect values to chronic human equivalent benchmark concentrations (BMCh) following the WHO/IPCS framework. Using ToxValDB chemicals with existing PODs associated with regulatory toxicity values, we found that the 25th %-ile of a chemical's BMCh distribution () could serve as a suitable surrogate for regulatory PODs (Q2 ≥ 0.76, RSE ≤ 0.82 log10 units). We applied this approach to derive for 2,095 substances with general non-cancer toxicity effects and 638 substances with reproductive/developmental toxicity effects, yielding a total coverage of 2,160 substances. From these , we derived probabilistic RfCs and human population effect concentrations. With this work, we have expanded the number of chemicals with toxicity values available, thereby enabling a much broader coverage for inhalation risk and impact assessment.

Development of an Automated Morphometric Approach to Assess Vascular Outcomes following Exposure to Environmental Chemicals in Zebrafish.

Disruptions in vascular formation attributable to chemical insults is a pivotal risk factor or potential etiology of developmental defects and various disease settings. Among the thousands of chemicals threatening human health, the highly concerning groups prevalent in the environment and detected in biological monitoring in the general population ought to be prioritized because of their high exposure risks. However, the impacts of a large number of environmental chemicals on vasculature are far from understood. The angioarchitecture complexity and technical limitations make it challenging to analyze the entire vasculature efficiently and identify subtle changes through a high-throughput in vivo assay. We aimed to develop an automated morphometric approach for the vascular profile and assess the vascular morphology of health-concerning environmental chemicals. High-resolution images of the entire vasculature in Tg(fli1a:eGFP) zebrafish were collected using a high-content imaging platform. We established a deep learning-based quantitative framework, ECA-ResXUnet, combined with MATLAB to segment the vascular networks and extract features. Vessel scores based on the rates of morphological changes were calculated to rank vascular toxicity. Potential biomarkers were identified by vessel-endothelium-gene-disease integrative analysis. Whole-trunk blood vessels and the cerebral vasculature in larvae exposed to 150 representative chemicals were automatically segmented as comparable to human-level accuracy, with sensitivity and specificity of 95.56% and 95.81%, respectively. Chemical treatments led to heterogeneous vascular patterns manifested by 31 architecture indexes, and the common cardinal vein (CCV) was the most affected vessel. The antipsychotic medicine haloperidol, flame retardant 2,2-bis(chloromethyl)trimethylenebis[bis(2-chloroethyl) phosphate], and tert-butylphenyl diphenyl phosphate ranked as the top three in vessel scores. Pesticides accounted for the largest group, with a vessel score of , characterized by a remarkable inhibition of subintestinal venous plexus and delayed development of CCV. Multiple-concentration evaluation of nine per- and polyfluoroalkyl substances (PFAS) indicated a low-concentration effect on vascular impairment and a positive association between carbon chain length and benchmark concentration. Target vessel-directed single-cell RNA sequencing of cells from larvae treated with , perfluorohexanesulfonic acid, or benzylbutyl phthalate, along with vessel-endothelium-gene-disease integrative analysis, uncovered potential associations with vascular disorders and identified biomarker candidates. This study provides a novel paradigm for phenotype-driven screenings of vascular-disrupting chemicals by converging morphological and transcriptomic profiles at a high-resolution level, serving as a powerful tool for large-scale toxicity tests. Our approach and the high-quality morphometric data facilitate the precise evaluation of vascular effects caused by environmental chemicals. https://doi.org/10.1289/EHP13214.

Comparison of transcriptomic profiles between HFPO-DA and prototypical PPARα, PPARγ, and cytotoxic agents in mouse, rat, and pooled human hepatocytes.

Like many per- or polyfluorinated alkyl substances (PFAS), toxicity studies with HFPO-DA (ammonium, 2,3,3,3-tetrafluoro-2-(heptafluoropropoxy)-propanoate), a short-chain PFAS used in the manufacture of some types of fluorinated polymers, indicate that the liver is the primary target of toxicity in rodents following oral exposure. Although the current weight of evidence supports the PPARα mode of action (MOA) for liver effects in HFPO-DA-exposed mice, alternate MOAs have also been hypothesized including PPARγ or cytotoxicity. To further evaluate the MOA for HFPO-DA in rodent liver, transcriptomic analyses were conducted on samples from primary mouse, rat, and pooled human hepatocytes treated for 12, 24, or 72 h with various concentrations of HFPO-DA, or agonists of PPARα (GW7647), PPARγ (rosiglitazone), or cytotoxic agents (ie, acetaminophen or d-galactosamine). Concordance analyses of enriched pathways across chemicals within each species demonstrated the greatest concordance between HFPO-DA and PPARα agonist GW7647-treated hepatocytes compared with the other chemicals evaluated. These findings were supported by benchmark concentration modeling and predicted upstream regulator results. In addition, transcriptomic analyses across species demonstrated a greater transcriptomic response in rodent hepatocytes treated with HFPO-DA or agonists of PPARα or PPARγ, indicating rodent hepatocytes are more sensitive to HFPO-DA or PPARα/γ agonist treatment. These results are consistent with previously published transcriptomic analyses and further support that liver effects in HFPO-DA-exposed rodents are mediated through rodent-specific PPARα signaling mechanisms as part of the MOA for PPARα activator-induced rodent hepatocarcinogenesis. Thus, effects observed in mouse liver are not appropriate endpoints for toxicity value development for HFPO-DA in human health risk assessment.

High-Throughput Transcriptomics of Nontumorigenic Breast Cells Exposed to Environmentally Relevant Chemicals.

There is a suite of chemicals, including metals, pesticides, and personal care product compounds, which are commonly detected at high levels in US Center for Disease Control's National Health and Nutrition Examination Survey (NHANES) chemical biomarker screens. Whether these chemicals influence development of breast cancer is not well understood. The objectives were to perform an unbiased concentration-dependent assessment of these chemicals, to quantify differences in cancer-specific genes and pathways, to describe if these differences occur at human population-relevant concentrations, and to specifically test for differences in markers of stemness and cellular plasticity. We treated nontumorigenic mammary epithelial cells, MCF10A, with 21 chemicals at four concentrations (, , , and ) for 48 h. We conducted RNA-sequencing for these 408 samples, adapting the plexWell plate-based RNA-sequencing method to analyze differences in gene expression. We calculated gene and biological pathway-specific benchmark concentrations (BMCs) using BMDExpress3, identifying differentially expressed genes and generating the best fit benchmark concentration models for each chemical across all genes. We identified enriched biological processes and pathways for each chemical and tested whether chemical exposures change predicted cell type distributions. We contextualized benchmark concentrations relative to human population biomarker concentrations in NHANES. We detected chemical concentration-dependent differences in gene expression for thousands of genes. Enrichment and cell type distribution analyses showed benchmark concentration responses correlated with differences in breast cancer-related pathways, including induction of basal-like characteristics for some chemicals, including arsenic, lead, copper, and methyl paraben. Comparison of benchmark data to NHANES chemical biomarker (urine or blood) concentrations indicated an overlap between exposure levels and levels sufficient to cause a gene expression response. These analyses revealed that many of these 21 chemicals resulted in differences in genes and pathways involved in breast cancer in vitro at human exposure-relevant concentrations. https://doi.org/10.1289/EHP12886.

Empirical Characterization of False Discovery Rates of Differentially Expressed Genes and Transcriptomic Benchmark Concentrations in Zebrafish Embryos.

High-throughput transcriptomics (HTTr) is increasingly applied to zebrafish embryos to survey the toxicological effects of environmental chemicals. Before the adoption of this approach in regulatory testing, it is essential to characterize background noise in order to guide experimental designs. We thus empirically quantified the HTTr false discovery rate (FDR) across different embryo pool sizes, sample sizes, and concentration groups for toxicology studies. We exposed zebrafish embryos to 0.1% dimethyl sulfoxide (DMSO) for 5 days. Pools of 1, 5, 10, and 20 embryos were created (n = 24 samples for each pool size). Samples were sequenced on the TempO-Seq platform and then randomly assigned to mock treatment groups before differentially expressed gene (DEG), pathway, and benchmark concentration (BMC) analyses. Given that all samples were treated with DMSO, any significant DEGs, pathways, or BMCs are false positives. As expected, we found decreasing FDRs for DEG and pathway analyses with increasing pool and sample sizes. Similarly, FDRs for BMC analyses decreased with increasing pool size and concentration groups, with more stringent BMC premodel filtering reducing BMC FDRs. Our study provides foundational data for determining appropriate experiment designs for regulatory toxicity testing with HTTr in zebrafish embryos.

A Bayesian benchmark concentration analysis for urinary fluoride and intelligence in adults in Guizhou, China

Environmental fluoride exposure has been linked to numerous cases of fluorosis worldwide. Previous studies have indicated that long-term exposure to fluoride can result in intellectual damage among children. However, a comprehensive health risk assessment of fluorosis-induced intellectual damage is still pending. In this research, we utilized the Bayesian Benchmark Dose Analysis System (BBMD) to investigate the dose-response relationship between urinary fluoride (U-F) concentration and Raven scores in adults from Nayong, Guizhou, China. Our research findings indecate a dose-response relationship between the concentration of U-F and intelligence scores in adults. As the benchmark response (BMR) increased, both the benchmark concentration (BMCs) and the lower bound of the credible interval (BMCLs) increased. Specifically, BMCs for the association between U-F and IQ score were determined to be 0.18 mg/L (BMCL1 = 0.08 mg/L), 0.91 mg/L (BMCL5 = 0.40 mg/L), 1.83 mg/L (BMCL10 = 0.83 mg/L) when using BMRs of 1 %, 5 %, and 10 %. These results indicate that U-F can serve as an effective biomarker for monitoring the loss of IQ in population. We propose three interim targets for public policy in preventing interllectual harm from fluoride exposure.

Mapping Interindividual Variability of Toxicodynamics Using High-Throughput Transcriptomics and Primary Human Hepatocytes from Fifty Donors.

Understanding the variability across the human population with respect to toxicodynamic responses after exposure to chemicals, such as environmental toxicants or drugs, is essential to define safety factors for risk assessment to protect the entire population. Activation of cellular stress response pathways are early adverse outcome pathway (AOP) key events of chemical-induced toxicity and would elucidate the estimation of population variability of toxicodynamic responses. We aimed to map the variability in cellular stress response activation in a large panel of primary human hepatocyte (PHH) donors to aid in the quantification of toxicodynamic interindividual variability to derive safety uncertainty factors. High-throughput transcriptomics of over 8,000 samples in total was performed covering a panel of 50 individual PHH donors upon 8 to 24 h exposure to broad concentration ranges of four different toxicological relevant stimuli: tunicamycin for the unfolded protein response (UPR), diethyl maleate for the oxidative stress response (OSR), cisplatin for the DNA damage response (DDR), and tumor necrosis factor alpha () for signaling. Using a population mixed-effect framework, the distribution of benchmark concentrations (BMCs) and maximum fold change were modeled to evaluate the influence of PHH donor panel size on the correct estimation of interindividual variability for the various stimuli. Transcriptome mapping allowed the investigation of the interindividual variability in concentration-dependent stress response activation, where the average of BMCs had a maximum difference of 864-, 13-, 13-, and 259-fold between different PHHs for UPR, OSR, DDR, and signaling-related genes, respectively. Population modeling revealed that small PHH panel sizes systematically underestimated the variance and gave low probabilities in estimating the correct human population variance. Estimated toxicodynamic variability factors of stress response activation in PHHs based on this dataset ranged between 1.6 and 6.3. Overall, by combining high-throughput transcriptomics and population modeling, improved understanding of interindividual variability in chemical-induced activation of toxicity relevant stress pathways across the human population using a large panel of plated cryopreserved PHHs was established, thereby contributing toward increasing the confidence of in vitro-based prediction of adverse responses, in particular hepatotoxicity. https://doi.org/10.1289/EHP11891.

Unlocking the Power of Transcriptomic Biomarkers in Qualitative and Quantitative Genotoxicity Assessment of Chemicals.

To modernize genotoxicity assessment and reduce reliance on experimental animals, new approach methodologies (NAMs) that provide human-relevant dose-response data are needed. Two transcriptomic biomarkers, GENOMARK and TGx-DDI, have shown a high classification accuracy for genotoxicity. As these biomarkers were extracted from different training sets, we investigated whether combining the two biomarkers in a human-derived metabolically competent cell line (i.e., HepaRG) provides complementary information for the classification of genotoxic hazard identification and potency ranking. First, the applicability of GENOMARK to TempO-Seq, a high-throughput transcriptomic technology, was evaluated. HepaRG cells were exposed for 72 h to increasing concentrations of 10 chemicals (i.e., eight known in vivo genotoxicants and two in vivo nongenotoxicants). Gene expression data were generated using the TempO-Seq technology. We found a prediction performance of 100%, confirming the applicability of GENOMARK to TempO-Seq. Classification using TGx-DDI was then compared to GENOMARK. For the chemicals identified as genotoxic, benchmark concentration modeling was conducted to perform potency ranking. The high concordance observed for both hazard classification and potency ranking by GENOMARK and TGx-DDI highlights the value of integrating these NAMs in a weight of evidence evaluation of genotoxicity.

Qualitative and quantitative concentration-response modelling of gene co-expression networks to unlock hepatotoxic mechanisms for next generation chemical safety assessment.

Next generation risk assessment of chemicals revolves around the use of mechanistic information without animal experimentation. In this regard, toxicogenomics has proven to be a useful tool to elucidate the underlying mechanisms of adverse effects of xenobiotics. In the present study, two widely used human in vitro hepatocyte culture systems, namely primary human hepatocytes (PHH) and human hepatoma HepaRG cells, were exposed to liver toxicants known to induce liver cholestasis, steatosis or necrosis. Benchmark concentration-response modelling was applied to transcriptomics gene co-expression networks (modules) to derive benchmark concentrations (BMCs) and to gain mechanistic insight into the hepatotoxic effects. BMCs derived by concentration-response modelling of gene co-expression modules recapitulated concentration-response modelling of individual genes. Although PHH and HepaRG cells showed overlap in deregulated genes and modules by the liver toxicants, PHH demonstrated a higher responsiveness, based on the lower BMCs of co-regulated gene modules. Such BMCs can be used as transcriptomics point of departure (tPOD) for assessing module-associated cellular (stress) pathways/processes. This approach identified clear tPODs of around maximum systemic concentration (Cmax) levels for the tested drugs, while for cosmetics ingredients the BMCs were 10-100-fold higher than the estimated plasma concentrations. This approach could serve next generation risk assessment practice to identify early responsive modules at low BMCs, that could be linked to key events in liver adverse outcome pathways. In turn, this can assist in delineating potential hazards of new test chemicals using in vitro systems and used in a risk assessment when BMCs are paired with chemical exposure assessment.

Investigating the dependency of in vitro benchmark concentrations on exposure time in transcriptomics experiments

There is increasing interest to employ in vitro transcriptomics experiments in toxicological testing, for example to determine a point-of-departure (PoD) for chemical safety assessment. However current practices to derive PoD tend to utilise a single exposure time despite the importance of exposure time on the manifestation of toxicity caused by a chemical. Therefore it is important to investigate both concentration and exposure time to determine how these factors affect biological responses, and as a consequence, the derivation of PoDs. In this study, metabolically competent HepaRG cells were exposed to five known toxicants over a range of concentrations and time points for subsequent gene expression analysis, using a targeted RNA expression assay (TempO-Seq). A non-parametric factor-modelling approach was used to model the collective response of all significant genes, which exploited the interdependence of differentially expressed gene responses. This in turn allowed the determination of an isobenchmark response (isoBMR) curve for each chemical in a reproducible manner. For 2 of the 5 chemicals tested, the PoD was observed to vary by 0.5–1 log-order within the 48-h timeframe of the experiment. The approach and findings presented here clearly demonstrate the need to take both concentration and exposure time into account when designing in vitro toxicogenomics experiments to determine PoD. Doing so also provides a means to use concentration-time-response modelling as a basis to extrapolate a PoD from shorter to longer exposure durations, and to identify chemicals of concern that can cause cumulative effects over time.

Toxicological Mechanisms and Potencies of Organophosphate Esters in KGN Human Ovarian Granulosa Cells as Revealed by High-throughput Transcriptomics.

Despite the growing number of studies reporting potential risks associated with exposure to organophosphate esters (OPEs), their molecular mechanisms of action remain poorly defined. We used the high-throughput TempO-Seq™ platform to investigate the effects of frequently detected OPEs on the expression of ∼3000 environmentally responsive genes in KGN human ovarian granulosa cells. Cells were exposed for 48 h to one of five OPEs (0.1 to 50 μM): tris(methylphenyl) phosphate (TMPP), isopropylated triphenyl phosphate (IPPP), tert-butylphenyl diphenyl phosphate (BPDP), triphenyl phosphate (TPHP), or tris(2-butoxyethyl) phosphate (TBOEP). The sequencing data indicate that four OPEs induced transcriptional changes, whereas TBOEP had no effect within the concentration range tested. Multiple pathway databases were used to predict alterations in biological processes based on differentially expressed genes. At lower concentrations, inhibition of the cholesterol biosynthetic pathway was the predominant effect of OPEs; this was likely a consequence of intracellular cholesterol accumulation. At higher concentrations, BPDP and TPHP had distinct effects, primarily affecting pathways involved in cell cycle progression and other stress responses. Benchmark concentration (BMC) modelling revealed that BPDP had the lowest transcriptomic point of departure. However, in vitro to in vivo extrapolation modeling indicated that TMPP was bioactive at lower concentrations than the other OPEs. We conclude that these new approach methodologies provide information on the mechanism(s) underlying the effects of data-poor compounds and assist in the derivation of protective points of departure for use in chemical read-across and decision-making.

Dose dependence of prenatal fluoride exposure associations with cognitive performance at school age in three prospective studies.

Fluoride may be a developmental neurotoxicant at elevated exposures. We merged new data from a prospective Odense Child Cohort (OCC) with results from two previous birth cohort studies from Mexico and Canada to characterize the dose-effect relationship in greater detail. The OCC contributed 837 mother-child pairs to the total of >1500. We measured creatinine-adjusted urine-fluoride concentrations in maternal urine samples obtained during late pregnancy. Child IQ was determined at age 7 years using an abbreviated version of the Wechsler Intelligence Scales for Children. Findings from the three cohorts were used to calculate the joint benchmark concentration (BMC) and the lower confidence limit (BMCL) after adjustment for covariables. In the OCC, urine-fluoride concentrations varied between 0.08 and 3.04 mg/l (median 0.52 mg/l) but were not significantly associated with full-scale IQ at age 7 years (β = 0.08; 95% confidence interval -1.14 to 1.30 for a doubling in exposure). No difference was apparent between boys and girls. In the OCC, the BMC was 0.92 mg/l, with a BMCL of 0.30 mg/l. The joint analysis of all three cohorts showed a statistically significant association between urine-fluoride and IQ, with a BMC of 0.45 mg/l (BMCL, 0.28 mg/l), slightly higher than the BMC previously reported for the two North American cohorts alone. As the BMCL reflects an approximate threshold for developmental neurotoxicity, the results suggest that pregnant women and children may need protection against fluoride toxicity.