The effects of combinations of limited ration and diazinon exposure on acetylcholinesterase activity, growth and reproduction in Oryzias latipes, the Japanese medaka.

Environmental contamination can negatively impact fish populations. In addition to acute toxicity leading to death, toxicants can reduce fish growth and lower reproduction. The potential for adverse population level effects of environmental contaminants are estimated to conduct risk assessments from laboratory toxicity tests that most often measure apical endpoints related to growth, survival and reproduction. The relationships between these effect endpoints are being evaluated to predict shifts in fish population demography better after exposure to environmental toxicants. Environmental contaminants can also affect fish populations indirectly by reducing prey biomass. However, estimating the magnitude of the combined effects of prey reduction and direct toxicity is difficult and rarely attempted. Here we describe a toxicity test designed to estimate the effect on Japanese medaka of both reduced food and chronic exposure to diazinon, an acetylcholinesterase inhibiting organophosphate pesticide. Fish were reared with limited food ration and/or diazinon exposure through a full life cycle to assess possible interactions between the two stressors in their effects on growth and reproduction. Diazinon exposure (10 or 20 µg/L), reduced ration (50% and 25% of ad libitum), or combinations of both lowered growth rates and reproductive output of Japanese medaka. In addition, growth and reproduction alone were modeled, and then various relationships between the two stressors (diazinon and ration) and how they relate to growth and reproduction were modeled.

Potential Toxicity of Complex Mixtures in Surface Waters from a Nationwide Survey of United States Streams: Identifying in Vitro Bioactivities and Causative Chemicals.

While chemical analysis of contaminant mixtures remains an essential component of environmental monitoring, bioactivity-based assessments using in vitro systems increasingly are used in the detection of biological effects. Historically, in vitro assessments focused on a few biological pathways, for example, aryl hydrocarbon receptor (AhR) or estrogen receptor (ER) activities. High-throughput screening (HTS) technologies have greatly increased the number of biological targets and processes that can be rapidly assessed. Here we screened extracts of surface waters from a nationwide survey of United States streams for bioactivities associated with 69 different end points using two multiplexed HTS assays. Bioactivity of extracts from 38 streams was evaluated and compared with concentrations of over 700 analytes to identify chemicals contributing to observed effects. Eleven primary biological end points were detected. Pregnane X receptor (PXR) and AhR-mediated activities were the most commonly detected. Measured chemicals did not completely account for AhR and PXR responses. Surface waters with AhR and PXR effects were associated with low intensity, developed land cover. Likewise, elevated bioactivities frequently associated with wastewater discharges included endocrine-related end points ER and glucocorticoid receptor. These results underscore the value of bioassay-based monitoring of environmental mixtures for detecting biological effects that could not be ascertained solely through chemical analyses.

Evidence for Cross Species Extrapolation of Mammalian-Based High-Throughput Screening Assay Results.

High-throughput screening (HTS) and computational technologies have emerged as important tools for chemical hazard identification. The US Environmental Protection Agency (EPA) launched the Toxicity ForeCaster (ToxCast) Program, which has screened thousands of chemicals in hundreds of mammalian-based HTS assays for biological activity. The data are being used to prioritize toxicity testing on those chemicals likely to lead to adverse effects. To use HTS assays in predicting hazard to both humans and wildlife, it is necessary to understand how broadly these data may be extrapolated across species. The US EPA Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS; https://seqapass.epa.gov/seqapass/ ) tool was used to assess conservation of the 484 protein targets represented in the suite of ToxCast assays and other HTS assays. To demonstrate the utility of the SeqAPASS data for guiding extrapolation, case studies were developed which focused on targets of interest to the US Endocrine Disruptor Screening Program and the Organisation for Economic Cooperation and Development. These case studies provide a line of evidence for conservation of endocrine targets across vertebrate species, with few exceptions, and demonstrate the utility of SeqAPASS for defining the taxonomic domain of applicability for HTS results and identifying organisms for suitable follow-up toxicity tests.

Gene transcription ontogeny of hypothalamic-pituitary-thyroid axis development in early-life stage fathead minnow and zebrafish.

The hypothalamic-pituitary-thyroid (HPT) axis is known to play a crucial role in the development of teleost fish. However, knowledge of endogenous transcription profiles of thyroid-related genes in developing teleosts remains fragmented. We selected two model teleost species, the fathead minnow (Pimephales promelas) and the zebrafish (Danio rerio), to compare the gene transcription ontogeny of the HPT axis. Control organisms were sampled at several time points during embryonic and larval development until 33 days post-fertilization. Total RNA was extracted from pooled, whole fish, and thyroid-related mRNA expression was evaluated using quantitative polymerase chain reaction. Gene transcripts examined included: thyrotropin-releasing hormone receptor (trhr), thyroid-stimulating hormone receptor (tshr), sodium-iodide symporter (nis), thyroid peroxidase (tpo), thyroglobulin (tg), transthyretin (ttr), deiodinases 1, 2, 3a, and 3b (dio1, dio2, dio3a and 3b), and thyroid hormone receptors alpha and beta (thrα and ß). A loess regression method was successful in identifying maxima and minima of transcriptional expression during early development of both species. Overall, we observed great similarities between the species, including maternal transfer, at least to some extent, of almost all transcripts (confirmed in unfertilized eggs), increasing expression of most transcripts during hatching and embryo-larval transition, and indications of a fully functional HPT axis in larvae. These data will aid in the development of hypotheses on the role of certain genes and pathways during development. Furthermore, this provides a background reference dataset for designing and interpreting targeted transcriptional expression studies both for fundamental research and for applications such as toxicology.

The effects of continuous diazinon exposure on growth and reproduction in Japanese medaka using a modified Medaka Extended One Generation Reproduction Test (MEOGRT).

The Medaka Extended One Generation Reproduction Test (MEOGRT) is a Tier 2 test within U.S. Environmental Protection Agency's (USEPA) Endocrine Disruptor Screening Program (EDSP), designed to characterize the potential adverse effects to fish of exposure to chemical that can cause disruption of the endocrine system. The MEOGRT focuses primarily on adverse effects to reproduction while collecting information regarding effects on growth, survival, and endocrine-related endpoints. However, the risk assessment process for fish, as mandated by legislation such as the Federal Insecticide, Fungicide, and Rodenticide Act (FIFRA) or the Toxic Substances Control Act (TSCA), could benefit from a more detailed assessment of effects on growth. Typically, fish growth data in support of risk assessment are obtained from full life-cycle tests or early life stage tests using the fathead minnow. As an alternative to these tests, a modified MEOGRT was conducted to assess the effects of diazinon on the various parameters measured in the MEOGRT. Diazinon is an organophosphate insecticide that is detected in the environment, and whose efficacy is a result of inhibition of the acetylcholine esterase enzyme at neuromuscular junctions and synapses of the nervous system. Diazinon (2.9, 5.2, 10.3, 19.8, and 40.2 µg/L) was tested with the MEOGRT protocol, and the lowest observable effect concentrations of 2.9 µg/L for fecundity and 5.2 µg/L for growth were determined. Additional growth measurements were added to the MEOGRT protocol to more robustly define growth rates and to determine the impact size has on reproductive performance. Fish size starting at the first measurement day (i.e. 21 days post-fertilization), and continuing through the duration of the test was reduced with exposure to 5.2 µg/L and higher, and asymptotic size predicted from growth modeling was reduced at 10.3 µg/L and higher. By simply adding non-destructive growth measurements at two additional time points, the MEOGRT provided enough data for the parameterization of growth models, which could be used to characterize the reproductive implications of growth impairment.

An "EAR" on Environmental Surveillance and Monitoring: A Case Study on the Use of Exposure-Activity Ratios (EARs) to Prioritize Sites, Chemicals, and Bioactivities of Concern in Great Lakes Waters.

Current environmental monitoring approaches focus primarily on chemical occurrence. However, based on concentration alone, it can be difficult to identify which compounds may be of toxicological concern and should be prioritized for further monitoring, in-depth testing, or management. This can be problematic because toxicological characterization is lacking for many emerging contaminants. New sources of high-throughput screening (HTS) data, such as the ToxCast database, which contains information for over 9000 compounds screened through up to 1100 bioassays, are now available. Integrated analysis of chemical occurrence data with HTS data offers new opportunities to prioritize chemicals, sites, or biological effects for further investigation based on concentrations detected in the environment linked to relative potencies in pathway-based bioassays. As a case study, chemical occurrence data from a 2012 study in the Great Lakes Basin along with the ToxCast effects database were used to calculate exposure-activity ratios (EARs) as a prioritization tool. Technical considerations of data processing and use of the ToxCast database are presented and discussed. EAR prioritization identified multiple sites, biological pathways, and chemicals that warrant further investigation. Prioritized bioactivities from the EAR analysis were linked to discrete adverse outcome pathways to identify potential adverse outcomes and biomarkers for use in subsequent monitoring efforts.

The influence of control group reproduction on the statistical power of the Environmental Protection Agency's Medaka Extended One Generation Reproduction Test (MEOGRT).

Because of various Congressional mandates to protect the environment from endocrine disrupting chemicals (EDCs), the United States Environmental Protection Agency (USEPA) initiated the Endocrine Disruptor Screening Program. In the context of this framework, the Office of Research and Development within the USEPA developed the Medaka Extended One Generation Reproduction Test (MEOGRT) to characterize the endocrine action of a suspected EDC. One important endpoint of the MEOGRT is fecundity of medaka breeding pairs. Power analyses were conducted to determine the number of replicates needed in proposed test designs and to determine the effects that varying reproductive parameters (e.g. mean fecundity, variance, and days with no egg production) would have on the statistical power of the test. The MEOGRT Reproduction Power Analysis Tool (MRPAT) is a software tool developed to expedite these power analyses by both calculating estimates of the needed reproductive parameters (e.g. population mean and variance) and performing the power analysis under user specified scenarios. Example scenarios are detailed that highlight the importance of the reproductive parameters on statistical power. When control fecundity is increased from 21 to 38 eggs per pair per day and the variance decreased from 49 to 20, the gain in power is equivalent to increasing replication by 2.5 times. On the other hand, if 10% of the breeding pairs, including controls, do not spawn, the power to detect a 40% decrease in fecundity drops to 0.54 from nearly 0.98 when all pairs have some level of egg production. Perhaps most importantly, MRPAT was used to inform the decision making process that lead to the final recommendation of the MEOGRT to have 24 control breeding pairs and 12 breeding pairs in each exposure group.

Addition of Phenylmethylsulfonyl Fluoride Increases the Working Lifetime of the Trout Liver S9 Substrate Depletion Assay, Resulting in Improved Detection of Low Intrinsic Clearance Rates.

The activity of a trout liver S9 substrate depletion assay has been shown to decline over time, presumably due to proteolytic degradation of biotransformation enzymes. To address this problem, assay performance was evaluated following the addition of phenylmethylsulfonyl fluoride (PMSF) or a general-purpose protease inhibitor cocktail to liver homogenization buffers and/or S9 reaction mixtures. Addition of PMSF to liver homogenization buffers and/or S9 reaction mixtures had little or no effect on clearance of phenanthrene, a model cytochrome P450 substrate, in short-term (25 or 30 min) depletion experiments but resulted in significant improvements in retention of this initial activity over time. The protease inhibitor cocktail strongly inhibited initial activity when added to homogenization buffers or reaction mixtures. Taking into consideration potential effects on liver carboxylesterases, the treatment approach determined to be optimal was addition of 10 µM PMSF to the S9 reaction mixture. Addition of 10 µM PMSF to the mixture resulted in significantly higher rates of phenanthrene clearance in 2-h incubations relative to those obtained in the absence of PMSF and a 6-fold increase in the working lifetime of the preparation. The results of a statistical power analysis suggest that by increasing the working lifetime of the assay, addition of PMSF to the reaction mixture could result in substantially improved detection of low in vitro clearance rates when compared to current practice. These findings demonstrate the value of adding PMSF to the trout S9 preparation and may have broad implications for use of this assay to support chemical bioaccumulation assessments for fish. Environ Toxicol Chem 2021;40:148-161. © 2020 SETAC. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Differentiating Pathway-Specific From Nonspecific Effects in High-Throughput Toxicity Data: A Foundation for Prioritizing Adverse Outcome Pathway Development.

The U.S. Environmental Protection Agency's ToxCast program has screened thousands of chemicals for biological activity, primarily using high-throughput in vitro bioassays. Adverse outcome pathways (AOPs) offer a means to link pathway-specific biological activities with potential apical effects relevant to risk assessors. Thus, efforts are underway to develop AOPs relevant to pathway-specific perturbations detected in ToxCast assays. Previous work identified a "cytotoxic burst" (CTB) phenomenon wherein large numbers of the ToxCast assays begin to respond at or near test chemical concentrations that elicit cytotoxicity, and a statistical approach to defining the bounds of the CTB was developed. To focus AOP development on the molecular targets corresponding to ToxCast assays indicating pathway-specific effects, we conducted a meta-analysis to identify which assays most frequently respond at concentrations below the CTB. A preliminary list of potentially important, target-specific assays was determined by ranking assays by the fraction of chemical hits below the CTB compared with the number of chemicals tested. Additional priority assays were identified using a diagnostic-odds-ratio approach which gives greater ranking to assays with high specificity but low responsivity. Combined, the two prioritization methods identified several novel targets (e.g., peripheral benzodiazepine and progesterone receptors) to prioritize for AOP development, and affirmed the importance of a number of existing AOPs aligned with ToxCast targets (e.g., thyroperoxidase, estrogen receptor, aromatase). The prioritization approaches did not appear to be influenced by inter-assay differences in chemical bioavailability. Furthermore, the outcomes were robust based on a variety of different parameters used to define the CTB.

Reliability of In Vitro Methods Used to Measure Intrinsic Clearance of Hydrophobic Organic Chemicals by Rainbow Trout: Results of an International Ring Trial.

In vitro assays are widely employed to obtain intrinsic clearance estimates used in toxicokinetic modeling efforts. However, the reliability of these methods is seldom reported. Here we describe the results of an international ring trial designed to evaluate two in vitro assays used to measure intrinsic clearance in rainbow trout. An important application of these assays is to predict the effect of biotransformation on chemical bioaccumulation. Six laboratories performed substrate depletion experiments with cyclohexyl salicylate, fenthion, 4-n-nonylphenol, deltamethrin, methoxychlor, and pyrene using cryopreserved hepatocytes and liver S9 fractions from trout. Variability within and among laboratories was characterized as the percent coefficient of variation (CV) in measured in vitro intrinsic clearance rates (CLIN VITRO, INT; ml/h/mg protein or 106 cells) for each chemical and test system. Mean intralaboratory CVs for each test chemical averaged 18.9% for hepatocytes and 14.1% for S9 fractions, whereas interlaboratory CVs (all chemicals and all tests) averaged 30.1% for hepatocytes and 22.4% for S9 fractions. When CLIN VITRO, INT values were extrapolated to in vivo intrinsic clearance estimates (CLIN VIVO, INT; l/d/kg fish), both assays yielded similar levels of activity (<4-fold difference for all chemicals). Hepatic clearance rates (CLH; l/d/kg fish) calculated using data from both assays exhibited even better agreement. These findings show that both assays are highly reliable and suggest that either may be used to inform chemical bioaccumulation assessments for fish. This study highlights several issues related to the demonstration of assay reliability and may provide a template for evaluating other in vitro biotransformation assays.

Summary of the development the US Environmental Protection Agency's Medaka Extended One Generation Reproduction Test (MEOGRT) using data from 9 multigenerational medaka tests.

In response to various legislative mandates, the US Environmental Protection Agency (USEPA) formed its Endocrine Disruptor Screening Program (EDSP), which in turn, formed the basis of a tiered testing strategy to determine the potential of pesticides, commercial chemicals, and environmental contaminants to disrupt the endocrine system. The first tier of tests is intended to detect the potential for endocrine disruption mediated through estrogen, androgen, or thyroid pathways, whereas the second tier is intended to further characterize the effects on these pathways and to establish a dose-response relationship for adverse effects. One of these tier 2 tests, the Medaka Extended One Generation Reproduction Test (MEOGRT), was developed by the USEPA for the EDSP and, in collaboration with the Japanese Ministry of the Environment, for the Guidelines for the Testing of Chemicals of the Organisation for Economic Co-operation and Development (OECD). The MEOGRT protocol was iteratively modified based on knowledge gained after the successful completion of 9 tests with variations in test protocols. The present study describes both the final MEOGRT protocol that has been published by the USEPA and the OECD, and the iterations that provided valuable insights into nuances of the protocol. The various tests include exposure to 17ß-estradiol, 4-t-octylphenol, o,p'- dichlorodiphenyltrichloroethane, 4-chloro-3-methylphenol, tamoxifen, 17ß-trenbolone, vinclozolin, and prochloraz. Environ Toxicol Chem 2017;36:3387-3403. Published 2017 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Statistical analysis of histopathological endpoints.

Histopathological assessments of fish from aquatic ecotoxicology studies are being performed with increasing frequency. Aquatic ecotoxicology studies performed for submission to regulatory agencies are usually conducted with multiple subjects (e.g., fish) in each of multiple vessels (replicates) within a water control and within each of several concentrations of a test substance. A number of histopathological endpoints are evaluated in each fish, and a severity score is generally recorded for each endpoint. The severity scores are often recorded using a nonquantitative scale of 0 to 4, with 0 indicating no effect, 1 indicating minimal effect, through 4 for severe effect. Statistical methods often used to analyze these scores suffer from several shortcomings: computing average scores as though scores were quantitative values, considering only the frequency of abnormality while ignoring severity, ignoring any concentration-response trend, and ignoring the possible correlation between responses of individuals within test vessels. A new test, the Rao-Scott Cochran-Armitage by Slices (RSCABS), is proposed that incorporates the replicate vessel experimental design and the biological expectation that the severity of the effect tends to increase with increasing doses or concentrations, while retaining the individual subject scores and taking into account the severity as well as frequency of scores. A power simulation and examples demonstrate the performance of the test. R-based software has been developed to carry out this test and is available free of charge at www.epa.gov/med/Prods_Pubs/rscabs.htm. The SAS-based RSCABS software is available from the first and third authors.

Use of gene expression data to determine effects on gonad phenotype in Japanese medaka after exposure to trenbolone or estradiol.

Various aquatic bioassays using one of several fish species have been developed or are in the process of being developed by organizations like the US Environmental Protection Agency and the Office of Economic Cooperation and Development for testing potential endocrine-disrupting chemicals (EDCs). Often, these involve assessment of the gonad phenotype of individuals as a key endpoint that is inputted into a risk or hazard assessment. Typically, gonad phenotype is determined histologically, which involves specialized and time-consuming techniques. The methods detailed here utilize an entirely different methodology, reverse-transcription quantitative polymerase chain reaction, to determine the relative expression levels of 4 genes after exposure to either 17ß-estradiol or 17ß-trenbolone and, by extension, the effects of EDCs on the phenotypic status of the gonad. The 4 genes quantified, Sox9b, protamine, Fig1α, and ZPC1, are all involved in gonad development and maintenance in Japanese medaka (Oryzias latipes); these data were then inputted into a permutational multivariate analysis of variance to determine whether significant differences exist between treatment groups. This information in conjunction with the sexual genotype, which can be determined in medaka, can be used to determine adverse effects of exposure to EDCs in a similar fashion to the histologically determined gonad phenotype.