Transcriptomic Points of Departure Calculated from Rainbow Trout Gill, Liver, and Gut Cell Lines Exposed to Methylmercury and Fluoxetine

TPD-seq: A high throughput RNA-seq method to derive transcriptomic points of departure from cell lines

There is growing scientific and regulatory interest in transcriptomic points of departure (tPOD) values from high-throughput in vitro experiments. To further help democratize tPOD research, here we outline 'TPD-seq' which links microplate-based exposure methods involving cell lines for human (Caco-2, Hep G2) and environmental (rainbow trout RTgill-W1) health, with a commercially available RNA-seq kit, with a cloud-based bioinformatics tool (ExpressAnalyst.ca). We applied the TPD-seq workflow to derive tPODs for solvents (dimethyl sulfoxide, DMSO; methanol) and positive controls (3,4-dichloroaniline, DCA; hydrogen peroxide, H2O2) commonly used in toxicity testing. The majority of reads mapped to protein coding genes (∼9 k for fish cells; ∼6 k for human cells), and about 50 % of differentially expressed genes were curve-fitted from which 90 % yielded gene benchmark doses. The most robust transcriptomic responses were caused by DMSO exposure, and tPOD values were 31-155 mM across the cell lines. OECD test guideline 249 (RTgill-W1 cells) recommends the use of DCA and here we calculated a tPOD of ∼5 to 76 μM. Finally, exposure of the two human cell lines to H2O2 resulted in tPOD values that ranged from 0.7 to 1.1 mM in Caco-2 cells and 5-30 μM in Hep G2 cells. The methods outlined here are designed to be performed in laboratories with basic molecular and cell culture facilities, and the performance and scalability of the TPD-seq workflow can be determined with additional case studies.

Cytotoxic and Molecular Effects of Soil Extracts from the Agbogbloshie Electronic-Waste Site on Fish and Human Cell Lines

Effect-based methods (EBM) are of growing interest in environmental monitoring programs. Few EBM have incorporated transcriptomics even though these provide a wealth of biological information and can be modeled to yield transcriptomic points of departure (tPODs). The study objectives were to: A) characterize cytotoxic effects of soil extracts on the rainbow trout RTgill-W1 and the human Caco-2 cell lines; B) measure gene expression changes and calculate tPODs; and C) compare in vitro responses to available measures of plastic-related compounds and metals. Extracts were prepared from 35 soil samples collected at the Agbogbloshie E-waste site (Accra, Ghana). Cells were exposed to six soil concentrations (0.3 to 9.4 mg dry weight of extract (eQsed)/ml). Many samples caused cytotoxicity with RTgill cells being more sensitive than Caco-2 cells. Eleven samples were analyzed for transcriptomics in both cell lines, with responses measured in all samples (52 to 5925 differentially expressed genes) even in the absence of cytotoxicity. In RTgill cells there was concordance between cytotoxic measures in tPOD values (spearman = 0.85). Though trends between in vitro measures and contaminant data were observed, more work is needed in this area before definitive conclusions are drawn. Nonetheless, this study helps support ongoing efforts in establishing alternative testing strategies (e.g., alternative to animal methods; toxicogenomics) for the assessment of complex environmental samples.

Transcriptomic Point of Departure (tPOD) of androstenedione in zebrafish embryos as a potential surrogate for chronic endpoints

The transcriptomic Point of Departure (tPOD) is increasingly used in ecotoxicology to derive quantitative endpoints from RNA sequencing studies. Utilizing transcriptomic data in zebrafish embryos as a New Approach Methodology (NAM) is beneficial due to its acknowledgment as an alternative to animal testing under EU Directive 2010/63/EU. Transcriptomic profiles are available in zebrafish for various modes of action (MoA). The limited literature available suggest that tPOD values from Fish Embryo Toxicity (FET) tests align with, but are generally lower than, No Observed Effect Concentrations (NOEC) from long-term chronic fish toxicity tests. In studies with the androgenic hormone androstenedione in a Fish Sexual Development Test (FSDT), a significant shift in the sex ratio towards males was noted at all test concentrations, making it impossible to determine a NOEC (NOEC <4.34 μg/L). To avoid additional animal testing in a repetition of the FSDT and adhere to the 3Rs principle (replacement, reduction, and refinement), a modified zebrafish FET (zFET) was conducted aiming to determine a regulatory acceptable effect threshold. This involved lower concentration ranges (20 to 6105 ng/L), overlapping with the masculinization-observed concentrations in the FSDT. The tPOD analysis in zFET showed consistent results with previous FSDT findings, observing strong expression changes in androgen-dependent genes at higher concentrations but not at lower ones, demonstrating a concentration-response relationship. The tPOD values for androstenedione were determined as 24 ng/L (10th percentile), 60 ng/L (20th gene), and 69 ng/L (1st peak). The 10th percentile tPOD value in zFET was 200 times lower than the lowest concentration in the FSDT. Comparing the tPOD values to literature suggests their potential to inform on the NOEC range in FSDT tests.

Far from Their Origins: A Transcriptomic Investigation on How 2,4-Di-tert-butyl-6-(5-chloro-2H-benzotriazol-2-yl) Phenol Affects Rainbow Trout Alevins

Benzotriazole ultraviolet stabilizers (BUVSs) are a group of widely used chemicals added to a variety of consumer (e.g., plastics) and industrial (e.g., metal coating) goods. Although detected globally as an environmentally persistent pollutant, BUVSs have received relatively little toxicological attention and only recently have been acknowledged to affect development and the endocrine system in vivo. In our previous study, altered behavior, indicative of potential neurotoxicity, was observed among rainbow trout alevins (day 14 posthatching) that were microinjected as embryos with a single environmentally relevant dose of 2,4-di-tert-butyl-6-(5-chloro-2H-benzotriazol-2-yl) phenol (UV-327). In the present follow-up study, we performed whole-transcriptome profiling (RNA sequencing) of newly hatched alevins from the same batch. The primary aim was to identify biomarkers related to behavior and neurology. Dose-specifically, 1 to 176 differentially expressed genes (DEGs) were identified. In the group presenting altered behavior (273.4 ng g-1), 176 DEGs were identified, yet only a fraction was related to neurological functions, including water, calcium, and potassium homeostasis; acetylcholine transmission and signaling; as well insulin and energy metabolism. The second objective was to estimate the transcriptomic point of departure (tPOD) and assess if point estimate(s) are protective of altered behavior. A tPOD was established at 35 to 94 ng UV-327 g-1 egg, making this tPOD protective of behavioral alterations. Holistically, these transcriptomic alterations provide a foundation for future research on how BUVSs can influence rainbow trout alevin development, while providing support to the hypothesis that UV-327 can influence neurogenesis and subsequent behavioral endpoints. The exact structural and functional changes caused by embryonic exposure to UV-327 remain enigmatic and will require extensive investigation before being deciphered and understood toxicologically. Environ Toxicol Chem 2024;43:2026-2038. © 2024 The Author(s). Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Comparing Transcriptomic Points of Departure to Apical Effect Concentrations For Larval Fathead Minnow Exposed to Chemicals with Four Different Modes Of Action

It is postulated that below a transcriptomic-based point of departure, adverse effects are unlikely to occur, thereby providing a chemical concentration to use in screening level hazard assessment. The present study extends previous work describing a high-throughput fathead minnow assay that can provide full transcriptomic data after exposure to a test chemical. One-day post-hatch fathead minnows were exposed to ten concentrations of three representatives of four chemical modes of action: organophosphates, ecdysone receptor agonists, plant photosystem II inhibitors, and estrogen receptor agonists for 24 h. Concentration response modeling was performed on whole body gene expression data from each exposure, using measured chemical concentrations when available. Transcriptomic points of departure in larval fathead minnow were lower than apical effect concentrations across fish species but not always lower than toxic effect concentrations in other aquatic taxa like crustaceans and insects. The point of departure was highly dependent on measured chemical concentration which were often lower than the nominal concentration. Differentially expressed genes between chemicals within modes of action were compared and often showed statistically significant overlap. In addition, reproducibility between identical exposures using a positive control chemical (CuSO4) and variability associated with the transcriptomic point of departure using in silico sampling were considered. Results extend a transcriptomic-compatible fathead minnow high-throughput assay for possible use in ecological hazard screening.

Framework for multi-stressor physiological response evaluation in amphibian risk assessment and conservation

Controlled laboratory experiments are often performed on amphibians to establish causality between stressor presence and an adverse outcome. However, in the field, identification of lab-generated biomarkers from single stressors and the interactions of multiple impacts are difficult to discern in an ecological context. The ubiquity of some pesticides and anthropogenic contaminants results in potentially cryptic sublethal effects or synergistic effects among multiple stressors. Although biochemical pathways regulating physiological responses to toxic stressors are often well-conserved among vertebrates, different exposure regimes and life stage vulnerabilities can yield variable ecological risk among species. Here we examine stress-related biomarkers, highlight endpoints commonly linked to apical effects, and discuss differences in ontogeny and ecology that could limit interpretation of biomarkers across species. Further we identify promising field-based physiological measures indicative of potential impacts to health and development of amphibians that could be useful to anuran conservation. We outline the physiological responses to common stressors in the context of altered functional pathways, presenting useful stage-specific endpoints for anuran species, and discussing multi-stressor vulnerability in the larger framework of amphibian life history and ecology. This overview identifies points of physiological, ecological, and demographic vulnerability to provide context in evaluating the multiple stressors impacting amphibian populations worldwide for strategic conservation planning.

Exploring the effects of experimental parameters and data modeling approaches on in vitro transcriptomic point-of-departure estimates

Multiple new approach methods (NAMs) are being developed to rapidly screen large numbers of chemicals to aid in hazard evaluation and risk assessments. High-throughput transcriptomics (HTTr) in human cell lines has been proposed as a first-tier screening approach for determining the types of bioactivity a chemical can cause (activation of specific targets vs. generalized cell stress) and for calculating transcriptional points of departure (tPODs) based on changes in gene expression. In the present study, we examine a range of computational methods to calculate tPODs from HTTr data, using six data sets in which MCF7 cells cultured in two different media formulations were treated with a panel of 44 chemicals for 3 different exposure durations (6, 12, 24 hr). The tPOD calculation methods use data at the level of individual genes and gene set signatures, and compare data processed using the ToxCast Pipeline 2 (tcplfit2), BMDExpress and PLIER (Pathway Level Information ExtractoR). Methods were evaluated by comparing to in vitro PODs from a validated set of high-throughput screening (HTS) assays for a set of estrogenic compounds. Key findings include: (1) for a given chemical and set of experimental conditions, tPODs calculated by different methods can vary by several orders of magnitude; (2) tPODs are at least as sensitive to computational methods as to experimental conditions; (3) in comparison to an external reference set of PODs, some methods give generally higher values, principally PLIER and BMDExpress; and (4) the tPODs from HTTr in this one cell type are mostly higher than the overall PODs from a broad battery of targeted in vitro ToxCast assays, reflecting the need to test chemicals in multiple cell types and readout technologies for in vitro hazard screening.

Cytotoxicity of 19 Pesticides in Rainbow Trout Gill, Liver, and Intestinal Cell Lines

The rainbow trout gill cell line (RTgill-W1), via test guideline 249 of the Organisation for Economic Co-operation and Development, has been established as a promising New Approach Methodology, although to advance confidence in the method more case studies are needed that: 1) expand our understanding of applicability domains (chemicals with diverse properties); 2) increase methodological throughput (96-well format); and 3) demonstrate biological relevance (in vitro to in vivo comparisons; gill vs. other cells). Accordingly, the objective of our study was to characterize the cytotoxicity of 19 pesticides against RTgill-W1 cells, and also liver (RTL-W1) and gut epithelial (RTgutGC) cell lines, and then to compare the in vitro and in vivo data. Of the 19 pesticides tested, 11, 9, and 8 were cytotoxic to the RTgill-W1, RTL-W1, and RTgutGC cells, respectively. Six pesticides (carbaryl, chlorothalonil, chlorpyrifos, dimethenamid-P, metolachlor, and S-metolachlor) were cytotoxic to all three cell lines. Aminomethylphosphonic acid, chlorantraniliprole, dicamba, diquat, imazethapyr, and permethrin exhibited cell-line-specific toxicity. No cytotoxic responses were observed for three herbicides (atrazine, glyphosate, and metribuzin) and four insecticides (clothianidin, diazinon, imidacloprid, and thiamethoxam). When cytotoxicity was measured, there was a strong correlation (rs  = 0.9, p < 0.0001) between in vitro median effect concentration (EC50) values (based on predicted concentrations using the In Vitro Mass Balance Model Equilibrium Partitioning (IV-MBM EQP) Ver. 2.1) derived from RTgill-W1 and RTL-W1 cells with in vivo median lethal concentration (LC50) values from 96-h acute toxicity studies with trout. In all 28 cases, the in vitro EC50 was within 18-fold of the in vivo LC50. These data help increase our understanding of the ecotoxicological domains of applicability for in vitro studies using cultured rainbow trout cells, while also demonstrating that these assays performed well in a 96-well format and have promise to yield data of biological relevance. Environ Toxicol Chem 2024;00:1-13. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Differences in the transcriptome response in the gills of sea lamprey acutely exposed to 3-trifluoromethyl-4-nitrophenol (TFM), niclosamide or a TFM:niclosamide mixture

Sea lamprey (Petromyzon marinus) control in the Laurentian Great Lakes of North America makes use of two pesticides: 3-trifluoromethyl-4-nitrophenol (TFM) and niclosamide, which are often co-applied. Sea lamprey appear to be vulnerable to these agents resulting from a lack of detoxification responses with evidence suggesting that lampricide mixtures produce a synergistic effect. However, there is a lack of information pertaining to the physiological responses of sea lamprey to niclosamide and TFM:niclosamide mixtures. Here, we characterized the transcriptomic responses of the sea lamprey to TFM, niclosamide, and a TFM:niclosamide (1.5 %) mixture in the gill. Along with a control, larval sea lamprey were exposed to each treatment for 6 h, after which gill tissues were extracted for measuring whole-transcriptome responses using RNA sequencing. Differential gene expression patterns were summarized, which included identifying the broad roles of genes and common expression patterns among the treatments. While niclosamide treatment resulted in no differentially expressed genes, TFM- and mixture-treated fish had several differentially expressed genes that were associated with the cell cycle, DNA damage, metabolism, immune function, and detoxification. However, there was no common differential expression among treatments. For the first time, we characterized the transcriptomic response of sea lamprey to niclosamide and a TFM:niclosamide mixture and identified that these agents impact mRNA transcript abundance of genes associated with the cell cycle and cellular death, and immune function, which are likely mediated through mitochondrial dysregulation. These results may help to inform the production of more targeted and effective lampricides in sea lamprey control efforts.

Short-Term Transcriptomic Points of Departure Are Consistent with Chronic Points of Departure for Three Organophosphate Pesticides across Mouse and Fathead Minnow

New approach methods (NAMs) can reduce the need for chronic animal studies. Here, we apply benchmark dose (concentration) (BMD(C))-response modeling to transcriptomic changes in the liver of mice and in fathead minnow larvae after short-term exposures (7 days and 1 day, respectively) to several dose/concentrations of three organophosphate pesticides (OPPs): fenthion, methidathion, and parathion. The mouse liver transcriptional points of departure (TPODs) for fenthion, methidathion, and parathion were 0.009, 0.093, and 0.046 mg/Kg-bw/day, while the fathead minnow larva TPODs were 0.007, 0.115, and 0.046 mg/L, respectively. The TPODs were consistent across both species and reflected the relative potencies from traditional chronic toxicity studies with fenthion identified as the most potent. Moreover, the mouse liver TPODs were more sensitive than or within a 10-fold difference from the chronic apical points of departure (APODs) for mammals, while the fathead minnow larva TPODs were within an 18-fold difference from the chronic APODs for fish species. Short-term exposure to OPPs significantly impacted acetylcholinesterase mRNA abundance (FDR p-value <0.05, |fold change| ≥2) and canonical pathways (IPA, p-value <0.05) associated with organism death and neurological/immune dysfunctions, indicating the conservation of key events related to OPP toxicity. Together, these results build confidence in using short-term, molecular-based assays for the characterization of chemical toxicity and risk, thereby reducing reliance on chronic animal studies.

Considerations for application of benchmark dose modeling in radiation research: workshop highlights

BACKGROUND

Exposure to different forms of ionizing radiation occurs in diverse occupational, medical, and environmental settings. Improving the accuracy of the estimated health risks associated with exposure is therefore, essential for protecting the public, particularly as it relates to chronic low dose exposures. A key aspect to understanding health risks is precise and accurate modeling of the dose-response relationship. Toward this vision, benchmark dose (BMD) modeling may be a suitable approach for consideration in the radiation field. BMD modeling is already extensively used for chemical hazard assessments and is considered statistically preferable to identifying low and no observed adverse effects levels. BMD modeling involves fitting mathematical models to dose-response data for a relevant biological endpoint and identifying a point of departure (the BMD, or its lower bound). Recent examples in chemical toxicology show that when applied to molecular endpoints (e.g. genotoxic and transcriptional endpoints), BMDs correlate to points of departure for more apical endpoints such as phenotypic changes (e.g. adverse effects) of interest to regulatory decisions. This use of BMD modeling may be valuable to explore in the radiation field, specifically in combination with adverse outcome pathways, and may facilitate better interpretation of relevant in vivo and in vitro dose-response data. To advance this application, a workshop was organized on June 3rd, 2022, in Ottawa, Ontario that brought together BMD experts in chemical toxicology and the radiation scientific community of researchers, regulators, and policy-makers. The workshop's objective was to introduce radiation scientists to BMD modeling and its practical application using case examples from the chemical toxicity field and demonstrate the BMDExpress software using a radiation dataset. Discussions focused on the BMD approach, the importance of experimental design, regulatory applications, its use in supporting the development of adverse outcome pathways, and specific radiation-relevant examples.

CONCLUSIONS

Although further deliberations are needed to advance the use of BMD modeling in the radiation field, these initial discussions and partnerships highlight some key steps to guide future undertakings related to new experimental work.

Characterizing Variability and Uncertainty Associated with Transcriptomic Dose-Response Modeling

Transcriptomics dose-response analysis (TDRA) has emerged as a promising approach for integrating toxicogenomics data into a risk assessment context; however, variability and uncertainty associated with experimental design are not well understood. Here, we evaluated n = 55 RNA-seq profiles derived from Japanese quail liver tissue following exposure to chlorpyrifos (0, 0.04, 0.1, 0.2, 0.4, 1, 2, 4, 10, 20, and 40 μg/g; n = 5 replicates per group) via egg injection. The full dataset was subsampled 637 times to generate smaller datasets with different dose ranges and spacing (designs A-E) and number of replicates (n = 2-5). TDRA of the 637 datasets revealed substantial variability in the gene and pathway benchmark doses, but relative stability in overall transcriptomic point-of-departure (tPOD) values when tPODs were calculated with the "pathway" and "mode" methods. Further, we found that tPOD values were more dependent on the dose range and spacing than on the number of replicates, suggesting that optimal experimental designs should use fewer replicates (n = 2 or 3) and more dose groups to reduce uncertainty in the results. Finally, tPOD values ranged by over ten times for all surveyed experimental designs and tPOD types, suggesting that tPODs should be interpreted as order-of-magnitude estimates.