Dose Response, Dosimetric, and Metabolic Evaluations of Replacement PFAS Perfluoro-(2,5,8-trimethyl-3,6,9-trioxadodecanoic) Acid (HFPO-TeA)

Few studies are available on the environmental and toxicological effects of perfluoroalkyl ether carboxylic acids (PFECAs), such as GenX, which are replacing legacy PFAS in manufacturing processes. To collect initial data on the toxicity and toxicokinetics of a longer-chain PFECA, male and female Sprague Dawley rats were exposed to perfluoro-(2,5,8-trimethyl-3,6,9-trioxadodecanoic) acid (HFPO-TeA) by oral gavage for five days over multiple dose levels (0.3-335.2 mg/kg/day). Clinically, we observed mortality at doses >17 mg/kg/day and body weight changes at doses ≤17 mg/kg/day. For the 17 mg/kg/day dose level, T3 and T4 thyroid hormone concentrations were significantly decreased (p < 0.05) from controls and HFPO-TeA plasma concentrations were significantly different between sexes. Non-targeted analysis of plasma and in vitro hepatocyte assay extractions revealed the presence of another GenX oligomer, perfluoro-(2,5-dimethyl-3,6-dioxanonanoic) acid (HFPO-TA). In vitro to in vivo extrapolation (IVIVE) parameterized with in vitro toxicokinetic data predicted steady-state blood concentrations that were within seven-fold of those observed in the in vivo study, demonstrating reasonable predictivity. The evidence of thyroid hormone dysregulation, sex-based differences in clinical results and dosimetry, and IVIVE predictions presented here suggest that the replacement PFECA HFPO-TeA induces a complex and toxic exposure response in rodents.

Organocatalyst treatment improves variant calling and mutant detection in archival clinical samples

Formalin fixation of biological specimens damages nucleic acids and limits their use in genomic analyses. Previously, we showed that RNA isolation with an organocatalyst (2-amino-5-methylphenyl phosphonic acid, used to speed up reversal of formalin-induced adducts) and extended heated incubation (ORGΔ) improved RNA-sequencing data from formalin-fixed paraffin-embedded (FFPE) tissue samples. The primary goal of this study was to evaluate whether ORGΔ treatment improves DNA-sequencing data from clinical FFPE samples. We isolated RNA and DNA ± ORGΔ from paired FFPE and frozen human renal and ovarian carcinoma specimens collected as part of the National Cancer Institute Biospecimen Pre-analytical Variables program. Tumor types were microscopically confirmed from adjacent tissue sections. Following extraction, DNA was fragmented and sequenced and differences were compared between frozen and FFPE sample pairs. Treatment with ORGΔ improved concurrent SNP calls in FFPE DNA compared to non-ORGΔ FFPE samples and enhanced confidence in SNP calls for all FFPE DNA samples, beyond that of matched frozen samples. In general, the concordant SNPs identified in paired frozen and FFPE DNA samples agreed for both genotype and homozygosity vs. heterozygosity of calls regardless of ORGΔ treatment. The increased confidence in ORGΔ FFPE DNA variant calls relative to the matched frozen DNA suggests a novel application of this method. With further optimization, this method may improve quality of DNA-sequencing data in FFPE as well as frozen tissue samples.

Developmental toxicity of Nafion byproduct 2 (NBP2) in the Sprague-Dawley rat with comparisons to hexafluoropropylene oxide-dimer acid (HFPO-DA or GenX) and perfluorooctane sulfonate (PFOS)

Nafion byproduct 2 (NBP2) is a polyfluoroalkyl ether sulfonic acid that was recently detected in surface water, drinking water, and human serum samples from monitoring studies in North Carolina, USA. We orally exposed pregnant Sprague-Dawley rats to NBP2 from gestation day (GD) 14-18 (0.1-30 mg/kg/d), GD17-21, and GD8 to postnatal day (PND) 2 (0.3-30 mg/kg/d) to characterize maternal, fetal, and postnatal effects. GD14-18 exposures were also conducted with perfluorooctane sulfonate (PFOS) for comparison to NBP2, as well as data previously published for hexafluoropropylene oxide-dimer acid (HFPO-DA or GenX). NBP2 produced stillbirth (30 mg/kg), reduced pup survival shortly after birth (10 mg/kg), and reduced pup body weight (10 mg/kg). Histopathological evaluation identified reduced glycogen stores in newborn pup livers and hepatocyte hypertrophy in maternal livers at ≥ 10 mg/kg. Exposure to NBP2 from GD14-18 reduced maternal serum total T₃ and cholesterol concentrations (30 mg/kg). Maternal, fetal, and neonatal liver gene expression was investigated using RT-qPCR pathway arrays, while maternal and fetal livers were also analyzed using TempO-Seq transcriptomic profiling. Overall, there was limited alteration of genes in maternal or F1 livers from NBP2 exposure with significant changes mostly occurring in the top dose group (30 mg/kg) associated with lipid and carbohydrate metabolism. Metabolomic profiling indicated elevated maternal bile acids for NBP2, but not HFPO-DA or PFOS, while all three reduced 3-indolepropionic acid. Maternal and fetal serum and liver NBP2 concentrations were similar to PFOS, but ∼10-30-fold greater than HFPO-DA concentrations at a given maternal oral dose. NBP2 is a developmental toxicant in the rat, producing neonatal mortality, reduced pup body weight, reduced pup liver glycogen, reduced maternal thyroid hormones, and altered maternal and offspring lipid and carbohydrate metabolism similar to other studied PFAS, with oral toxicity for pup loss that is slightly less potent than PFOS but more potent than HFPO-DA.

Genomic comparisons between hepatocarcinogenic and non-hepatocarcinogenic organophosphate insecticides in the mouse liver

Short-term biomarkers of toxicity have an increasingly important role in the screening and prioritization of new chemicals. In this study, we examined early indicators of liver toxicity for three reference organophosphate (OP) chemicals, which are among the most widely used insecticides in the world. The OP methidathion was previously shown to increase the incidence of liver toxicity, including hepatocellular tumors, in male mice. To provide insights into the adverse outcome pathway (AOP) that underlies these tumors, effects of methidathion in the male mouse liver were examined after 7 and 28 day exposures and compared to those of two other OPs that either do not increase (fenthion) or possibly suppress liver cancer (parathion) in mice. None of the chemicals caused increases in liver weight/body weight or histopathological changes in the liver. Parathion decreased liver cell proliferation after 7 and 28 days while the other chemicals had no effects. There was no evidence for hepatotoxicity in any of the treatment groups. Full-genome microarray analysis of the livers from the 7 and 28 day treatments demonstrated that methidathion and fenthion regulated a large number of overlapping genes, while parathion regulated a unique set of genes. Examination of cytochrome P450 enzyme activities and use of predictive gene expression biomarkers found no consistent evidence for activation of AhR, CAR, PXR, or PPARα. Parathion suppressed the male-specific gene expression pattern through STAT5b, similar to genetic and dietary conditions that decrease liver tumor incidence in mice. Overall, these findings indicate that methidathion causes liver cancer by a mechanism that does not involve common mechanisms of liver cancer induction.

Direct formalin fixation induces widespread transcriptomic effects in archival tissue samples

Sequencing technologies now provide unprecedented access to genomic information in archival formalin-fixed paraffin-embedded (FFPE) tissue samples. However, little is known about artifacts induced during formalin fixation, which could bias results. Here we evaluated global changes in RNA-sequencing profiles between matched frozen and FFPE samples. RNA-sequencing was performed on liver samples collected from mice treated with a reference chemical (phenobarbital) or vehicle control for 7 days. Each sample was divided into four parts: (1) fresh-frozen, (2) direct-fixed in formalin for 18 h, (3) frozen then formalin-fixed, and (4) frozen then ethanol-fixed and paraffin-embedded (n = 6/group/condition). Direct fixation resulted in 2,946 differentially expressed genes (DEGs) vs. fresh-frozen, 98% of which were down-regulated. Freezing prior to formalin fixation had ≥ 95% fewer DEGs vs. direct fixation, indicating that most formalin-derived transcriptional effects in the liver occurred during fixation. This finding was supported by retrospective studies of paired frozen and FFPE samples, which identified consistent enrichment in oxidative stress, mitochondrial dysfunction, and transcription initiation pathways with direct fixation. Notably, direct formalin fixation in the parent study did not significantly impact response profiles resulting from chemical exposure. These results advance our understanding of FFPE samples as a resource for genomic research.

Enhanced Quality Metrics for Assessing RNA Derived From Archival Formalin-Fixed Paraffin-Embedded Tissue Samples

Formalin-fixed paraffin-embedded (FFPE) tissues provide an important resource for toxicogenomic research. However, variability in the integrity or quality of RNA obtained from archival FFPE specimens can lead to unreliable data and wasted resources, and standard protocols for measuring RNA integrity do not adequately assess the suitability of FFPE RNA. The main goal of this study was to identify improved methods for evaluating FFPE RNA quality for whole-genome sequencing. We examined RNA quality metrics conducted prior to RNA-sequencing in paired frozen and FFPE samples with varying levels of quality based on age in block and time in formalin. RNA quality was measured by the RNA integrity number (RIN), a modified RIN called the paraffin-embedded RNA metric, the percentage of RNA fragments >100-300 nucleotides in size (DV100-300), and 2 quantitative PCR-based methods. This information was correlated to sequencing read quality, mapping, and gene detection. Among fragmentation-based methods, DV and PCR-based metrics were more informative than RIN or paraffin-embedded RNA metric in determining sequencing success. Across low- and high-quality FFPE samples, a minimum of 80% of RNA fragments >100 nucleotides (DV100 > 80) provided the best indication of gene diversity and read counts upon sequencing. The PCR-based methods further showed quantitative reductions in amplifiable RNA of target genes related to sample age and time in formalin that inform input quantity of FFPE RNA for sequencing. These results should aid in screening and prioritizing archival FFPE samples for retrospective analyses of gene expression.

Microbiota alter metabolism and mediate neurodevelopmental toxicity of 17β-estradiol

Estrogenic chemicals are widespread environmental contaminants associated with diverse health and ecological effects. During early vertebrate development, estrogen receptor signaling is critical for many different physiologic responses, including nervous system function. Recently, host-associated microbiota have been shown to influence neurodevelopment. Here, we hypothesized that microbiota may biotransform exogenous 17-βestradiol (E2) and modify E2 effects on swimming behavior. Colonized zebrafish were continuously exposed to non-teratogenic E2 concentrations from 1 to 10 days post-fertilization (dpf). Changes in microbial composition and predicted metagenomic function were evaluated. Locomotor activity was assessed in colonized and axenic (microbe-free) zebrafish exposed to E2 using a standard light/dark behavioral assay. Zebrafish tissue was collected for chemistry analyses. While E2 exposure did not alter microbial composition or putative function, colonized E2-exposed larvae showed reduced locomotor activity in the light, in contrast to axenic E2-exposed larvae, which exhibited normal behavior. Measured E2 concentrations were significantly higher in axenic relative to colonized zebrafish. Integrated peak area for putative sulfonated and glucuronidated E2 metabolites showed a similar trend. These data demonstrate that E2 locomotor effects in the light phase are dependent on the presence of microbiota and suggest that microbiota influence chemical E2 toxicokinetics. More broadly, this work supports the concept that microbial colonization status may influence chemical toxicity.

Demodifying RNA for Transcriptomic Analyses of Archival Formalin-Fixed Paraffin-Embedded Samples

Archival formalin-fixed paraffin-embedded (FFPE) tissue samples offer a vast but largely untapped resource for genomic research. The primary technical issues limiting use of FFPE samples are RNA yield and quality. In this study, we evaluated methods to demodify RNA highly fragmented and crosslinked by formalin fixation. Primary endpoints were RNA recovery, RNA-sequencing quality metrics, and transcriptional responses to a reference chemical (phenobarbital, PB). Frozen mouse liver samples from control and PB groups (n = 6/group) were divided and preserved for 3 months as follows: frozen (FR); 70% ethanol (OH); 10% buffered formalin for 18 h followed by ethanol (18F); or 10% buffered formalin (3F). Samples from OH, 18F, and 3F groups were processed to FFPE blocks and sectioned for RNA isolation. Additional sections from 3F received the following demodification protocols to mitigate RNA damage: short heated incubation with Tris-Acetate-EDTA buffer; overnight heated incubation with an organocatalyst using 2 different isolation kits; or overnight heated incubation without organocatalyst. Ribo-depleted, stranded, total RNA libraries were built and sequenced using the Illumina HiSeq 2500 platform. Overnight incubation (± organocatalyst) increased RNA yield >3-fold and RNA integrity numbers and fragment analysis values by > 1.5- and >3.0-fold, respectively, versus 3F. Postsequencing metrics also showed reduced bias in gene coverage and deletion rates for overnight incubation groups. All demodification groups had increased overlap for differentially expressed genes (77%-84%) and enriched pathways (91%-97%) with FR, with the highest overlap in the organocatalyst groups. These results demonstrate simple changes in RNA isolation methods that can enhance genomic analyses of FFPE samples.

Metabolic Disruption Early in Life is Associated With Latent Carcinogenic Activity of Dichloroacetic Acid in Mice

Early-life environmental factors can influence later-life susceptibility to cancer. Recent evidence suggests that metabolic pathways may mediate this type of latency effect. Previously, we reported that short-term exposure to dichloroacetic acid (DCA) increased liver cancer in mice 84 weeks after exposure was stopped. Here, we evaluated time course dynamics for key events related to this effect. This study followed a stop-exposure design in which 28-day-old male B6C3F1 mice were given the following treatments in drinking water for up to 93 weeks: deionized water (dH2O, control); 3.5 g/l DCA continuously; or 3.5 g/l DCA for 4-52 weeks followed by dH2O. Effects were evaluated at eight interim time points. A short-term biomarker study was used to evaluate DCA effects at 6, 15, and 30 days. Liver tumor incidence was higher in all DCA treatment groups, including carcinomas in 82% of mice previously treated with DCA for only 4 weeks. Direct effects of DCA in the short-term study included decreased liver cell proliferation and marked mRNA changes related to mitochondrial dysfunction and altered cell metabolism. However, all observed short-term effects of DCA were ultimately reversible, and prior DCA treatment did not affect liver cell proliferation, apoptosis, necrosis, or DNA sequence variants with age. Key intermediate events resulting from transient DCA exposure do not fit classical cytotoxic, mitogenic, or genotoxic modes of action for carcinogenesis, suggesting a distinct mechanism associated with early-life metabolic disruption.

Effects of the antimicrobial contaminant triclocarban, and co-exposure with the androgen 17β-trenbolone, on reproductive function and ovarian transcriptome of the fathead minnow (Pimephales promelas)

Triclocarban (TCC) is an antimicrobial agent routinely detected in surface waters that has been hypothesized to interact with the vertebrate endocrine system. The present study examined the effects of TCC alone and in combination with the model endocrine disruptor 17β-trenbolone (TRB) on fish reproductive function. Adult Pimephales promelas were continuously exposed to either 1 µg TCC/L or 5 µg TCC/L, to 0.5 µg TRB/L, or to a mixture (MIX) of 5 µg TCC/L and 0.5 µg TRB/L for 22 d, and a variety of reproductive and endocrine-related endpoints were examined. Cumulative fecundity was significantly reduced in fathead minnows exposed to TRB, MIX, or 5 µg TCC/L. Exposure to 1 µg TCC/L had no effect on reproduction. In general, both TRB and MIX treatments caused similar physiological effects, evoking significant reductions in female plasma vitellogenin, estradiol, and testosterone, and significant increases in male plasma estradiol. Based on analysis of the ovarian transcriptome, there were potential pathway impacts that were common to both TRB- and TCC-containing treatment groups. In most cases, however, those pathways were more plausibly linked to differences in reproductive status than to androgen-specific functions. Overall, TCC was reproductively toxic to fish at concentrations at or near those that have been measured in surface water. There was little evidence that TCC elicits reproductive toxicity through a specific mode of endocrine or reproductive action, nor that it could augment the androgenic effects of TRB. Nonetheless, the relatively small margin of safety between some measured environmental concentrations and effect concentrations suggests that concern is warranted. Environ Toxicol Chem 2017;36:231-242. Published 2016 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.

Editor's Highlight: Dose-Response Analysis of RNA-Seq Profiles in Archival Formalin-Fixed Paraffin-Embedded Samples

Use of archival resources has been limited to date by inconsistent methods for genomic profiling of degraded RNA from formalin-fixed paraffin-embedded (FFPE) samples. RNA-sequencing offers a promising way to address this problem. Here, we evaluated transcriptomic dose responses using RNA-sequencing in paired FFPE and frozen (FROZ) samples from 2 archival studies in mice, one <2 years old and the other >20 years old. Experimental treatments included 3 different doses of di(2-ethylhexyl)phthalate or dichloroacetic acid for the recently archived and older studies, respectively. Total RNA was ribo-depleted and sequenced using the Illumina HiSeq platform. In the recently archived study, FFPE samples had 35% lower total counts compared to FROZ samples but high concordance in fold-change values of differentially expressed genes (DEGs) (r²  =0.99), highly enriched pathways (90% overlap with FROZ), and benchmark dose estimates for preselected target genes (<5% difference vs FROZ). In contrast, older FFPE samples had markedly lower total counts (3% of FROZ) and poor concordance in global DEGs and pathways. However, counts from FFPE and FROZ samples still positively correlated (r² = 0.84 across all transcripts) and showed comparable dose responses for more highly expressed target genes. These findings highlight potential applications and issues in using RNA-sequencing data from FFPE samples. Recently archived FFPE samples were highly similar to FROZ samples in sequencing quality metrics, DEG profiles, and dose-response parameters, while further methods development is needed for older lower-quality FFPE samples. This work should help advance the use of archival resources in chemical safety and translational science.

Phenotype anchoring in zebrafish reveals a potential role for matrix metalloproteinases (MMPs) in tamoxifen's effects on skin epithelium

The zebrafish is a powerful alternative model used to link phenotypes with molecular effects to discover drug mode of action. Using a zebrafish embryo-larval toxicity bioassay, we evaluated the effects of tamoxifen--a widely used anti-estrogen chemotherapeutic. Zebrafish exposed to ≥ 10 μM tamoxifen exhibited a unique necrotic caudal fin phenotype that was rapidly induced regardless of developmental life-stage when treatment was applied. To define tamoxifen's bioactivity resulting in this phenotype, targeted gene expression was used to evaluate 100 transcripts involved in tissue remodeling, calcium signaling, cell cycle and cell death, growth factors, angiogenesis and hypoxia. The most robustly misregulated transcripts in the tail were matrix metalloproteinases mmp9 and mmp13a, induced 127 and 1145 fold, respectively. Expression of c-fos, c-jun, and ap1s1 were also moderately elevated (3-7 fold), consistent with AP-1 activity--a transcription factor that regulates MMP expression. Immunohistochemistry confirmed high levels of induction for MMP13a in affected caudal fin skin epithelial tissue. The necrotic caudal fin phenotype was significantly attenuated or prevented by three functionally unique MMP inhibitors: EDTA (metal chelator), GM 6001 (broad MMP inhibitor), and SR 11302 (AP-1 transcription factor inhibitor), suggesting MMP-dependence. SR 11302 also inhibited induction of mmp9, mmp13a, and a putative MMP target, igfbp1a. Overall, our studies suggest that tamoxifen's effect is the result of perturbation of the MMP system in the skin leading to ectopic expression, cytotoxicity, and the necrotic caudal fin phenotype. These studies help advance our understanding of tamoxifen's non-classical mode of action and implicate a possible role for MMPs in tissues such as skin.

Toward safer multi-walled carbon nanotube design: Establishing a statistical model that relates surface charge and embryonic zebrafish mortality

Given the increased utility and lack of consensus regarding carbon nanotube (CNT) environmental and human health hazards, there is a growing demand for guidelines that inform safer CNT design. In this study, the zebrafish (Danio rerio) model is utilized as a stable, sensitive biological system to evaluate the bioactivity of systematically modified and comprehensively characterized multi-walled carbon nanotubes (MWNTs). MWNTs were treated with strong acid to introduce oxygen functional groups, which were then systematically thermally reduced and removed using an inert temperature treatment. While 25 phenotypic endpoints were evaluated at 24 and 120 hours post-fertilization (hpf), high mortality at 24 hpf prevented further resolution of the mode of toxicity leading to mortality. Advanced multivariate statistical methods are employed to establish a model that identifies those MWNT physicochemical properties that best estimate the probability of observing an adverse outcome. The physicochemical properties considered in this study include surface charge, percent surface oxygen, dispersed aggregate size and morphology and electrochemical activity. Of the five physicochemical properties, surface charge, quantified as the point of zero charge (PZC), was determined as the best predictor of mortality at 24 hpf. From a design perspective, the identification of this property-hazard relationship establishes a foundation for the development of design guidelines for MWNTs with reduced hazard.

Comparative Metal Oxide Nanoparticle Toxicity Using Embryonic Zebrafish

Engineered metal oxide nanoparticles (MO NPs) are finding increasing utility in the medical field as anticancer agents. Before validation of in vivo anticancer efficacy can occur, a better understanding of whole-animal toxicity is required. We compared the toxicity of seven widely used semiconductor MO NPs made from zinc oxide (ZnO), titanium dioxide, cerium dioxide and tin dioxide prepared in pure water and in synthetic seawater using a five-day embryonic zebrafish assay. We hypothesized that the toxicity of these engineered MO NPs would depend on physicochemical properties. Significant agglomeration of MO NPs in aqueous solutions is common making it challenging to associate NP characteristics such as size and charge with toxicity. However, data from our agglomerated MO NPs suggests that the elemental composition and dissolution potential are major drivers of toxicity. Only ZnO caused significant adverse effects of all MO particles tested, and only when prepared in pure water (point estimate median lethal concentration = 3.5-9.1 mg/L). This toxicity was life stage dependent. The 24 h toxicity increased greatly (~22.7 fold) when zebrafish exposures started at the larval life stage compared to the 24 hour toxicity following embryonic exposure. Investigation into whether dissolution could account for ZnO toxicity revealed high levels of zinc ion (40-89% of total sample) were generated. Exposure to zinc ion equivalents revealed dissolved Zn²⁺ may be a major contributor to ZnO toxicity.

Developing predictive approaches to characterize adaptive responses of the reproductive endocrine axis to aromatase inhibition: I. Data generation in a small fish model

Adaptive or compensatory responses to chemical exposure can significantly influence in vivo concentration-duration-response relationships. This study provided data to support development of a computational dynamic model of the hypothalamic-pituitary-gonadal axis of a model vertebrate and its response to aromatase inhibitors as a class of endocrine active chemicals. Fathead minnows (Pimephales promelas) were either exposed to the aromatase inhibitor fadrozole (0.5 or 30 μg/l) continuously for 1, 8, 12, 16, 20, 24, or 28 days or exposed for 8 days and then held in control water (no fadrozole) for an additional 4, 8, 12, 16, or 20 days. The time course of effects on ovarian steroid production, circulating 17β-estradiol (E2) and vitellogenin (VTG) concentrations, and expression of steroidogenesis-related genes in the ovary was measured. Exposure to 30 μg fadrozole/l significantly reduced plasma E2 and VTG concentrations after just 1 day and those effects persisted throughout 28 days of exposure. In contrast, ex vivo E2 production was similar to that of controls on day 8-28 of exposure, whereas transcripts coding for aromatase and follicle-stimulating hormone receptor were elevated, suggesting a compensatory response. Following cessation of fadrozole exposure, ex vivo E2 and plasma E2 concentrations exceeded and then recovered to control levels, but plasma VTG concentrations did not, even after 20 days of depuration. Collectively these data provide several new insights into the nature and time course of adaptive responses to an aromatase inhibitor that support development of a computational model (see companion article).

A time-course analysis of effects of the steroidogenesis inhibitor ketoconazole on components of the hypothalamic-pituitary-gonadal axis of fathead minnows

The objective of this study was to evaluate temporal effects of the model steroidogenesis inhibitor ketoconazole (KTC) on aspects of reproductive endocrine function controlled by the hypothalamic-pituitary-gonadal (HPG) axis in the fathead minnow (Pimephales promelas). Ketoconazole inhibits the activity of two cytochrome P450s (CYPs) key to sex steroid production in vertebrates, CYP11a (cholesterol side chain cleavage) and CYP17 (c17α-hydroxylase/17, 20-lyase). Sexually mature fish were exposed to water-borne KTC (30 or 300 μg/L) in a flow-through system for up to 8d, following which animals were allowed to recover in clean water. Fish were sampled after 1, 4 and 8d of exposure, and after 1, 8 and 16d of recovery. A shorter-term time-course experiment also was conducted in which females were sampled on seven occasions during a 12h KTC exposure. Ketoconazole consistently depressed ex vivo gonadal synthesis of testosterone (T) in both sexes, and 17β-estradiol (E2) in females during both exposure and recovery phases of the time-course studies. Effects on ex vivo steroidogenesis in females occurred within as little as 1h of exposure. Plasma concentrations of T in males and E2 in females also were depressed by exposure to KTC, but these decreases did not persist to the same degree as observed for the ex vivo effects. In females, after decreases within 12h, plasma E2 concentrations were similar to (or greater than) controls at 24h of exposure, while in males, plasma T returned to levels comparable to controls within 1d of cessation of KTC exposure. The discrepancy between the ex vivo and in vivo data at later stages in the test is consistent with some type of compensatory response to KTC in fish. However, we were unable to ascertain the mechanistic basis for such a response. For example, although a number of genes related to steroid synthesis (e.g., cyp11a, cyp17) were up-regulated in the gonads of both males and females during the exposure and early recovery phases of the experiment, this did not seem to account for the resurgence in plasma steroid concentrations in KTC-exposed fish. Further studies focused on metabolism and clearance of steroids might lend insights as to the effects of KTC on plasma steroid concentrations. Overall, our results demonstrate the complex, temporally dynamic nature of the vertebrate HPG system in response to chemical stressors.

Effects of a glucocorticoid receptor agonist, dexamethasone, on fathead minnow reproduction, growth, and development

Synthetic glucocorticoids are pharmaceutical compounds prescribed in human and veterinary medicine as anti-inflammatory agents and have the potential to contaminate natural watersheds via inputs from wastewater treatment facilities and confined animal-feeding operations. Despite this, few studies have examined the effects of this class of chemicals on aquatic vertebrates. To generate data to assess potential risk to the aquatic environment, we used fathead minnow 21-d reproduction and 29-d embryo-larvae assays to determine reproductive toxicity and early-life-stage effects of dexamethasone. Exposure to 500 µg dexamethasone/L in the 21-d test caused reductions in fathead minnow fecundity and female plasma estradiol concentrations and increased the occurrence of abnormally hatched fry. Female fish exposed to 500 µg dexamethasone/L also displayed a significant increase in plasma vitellogenin protein levels, possibly because of decreased spawning. A decrease in vitellogenin messenger ribonucleic acid (mRNA) expression in liver tissue from females exposed to the high dexamethasone concentration lends support to this hypothesis. Histological results indicate that a 29-d embryo-larval exposure to 500 µg dexamethasone/L caused a significant increase in deformed gill opercula. Fry exposed to 500 µg dexamethasone/L for 29 d also exhibited a significant reduction in weight and length compared with control fry. Taken together, these results indicate that nonlethal concentrations of a model glucocorticoid receptor agonist can impair fish reproduction, growth, and development.

Temporal evaluation of effects of a model 3β-hydroxysteroid dehydrogenase inhibitor on endocrine function in the fathead minnow

Inhibition of enzymes involved in the synthesis of sex steroids can substantially impact developmental and reproductive processes controlled by the hypothalmic-pituitary-gonadal (HPG) axis. A key steroidogenic enzyme that has received little attention from a toxicological perspective is 3β-hydroxysteroid dehydrogenase (3β-HSD). In these studies, we exposed reproductively-active fathead minnows (Pimephales promelas) to the model 3β-HSD inhibitor trilostane at two test concentrations (300 and 1,500 µg/L) over a 16-d period that included both 8-d exposure and 8-d recovery phases. Plasma concentrations of 17β-estradiol (E2) in females were depressed within hours of exposure to the drug and remained decreased at the highest trilostane concentration throughout the 8-d exposure. Reductions in E2 were accompanied by decreases in plasma concentrations of the estrogen-responsive protein vitellogenin (VTG). During the recovery phase of the test, plasma E2 and VTG concentrations returned to levels comparable to those of controls, in the case of E2 within 1 d. Up-regulation of ovarian expression of gene products for follicle-stimulating hormone receptor (fshr) and aromatase (cyp19a1a) suggested active compensation in trilostane-exposed animals. Effects of trilostane on HPG-related endpoints in exposed males were less pronounced, although, as in females, up-regulation of gonadal fshr was seen. Data from these time-course studies provide insights as to direct impacts, compensatory responses, and recovery from effects associated with perturbation of a comparatively poorly characterized enzyme/pathway critical to sex steroid synthesis. This information is important to the design and interpretation of approaches for assessing the occurrence and effects of HPG-active chemicals in both the laboratory and the field.

A method for the determination of genetic sex in the fathead minnow, Pimephales promelas, to support testing of endocrine-active chemicals

Certain endocrine-active toxicants have been reported to completely sex reverse both male and female individuals in amphibian, avian, fish, invertebrate, and reptile species, resulting in a phenotype indistinguishable from unaffected individuals. Detection of low-level sex reversal often requires large numbers of organisms to achieve the necessary statistical power, especially in those species with predominantly genetic sex determination and cryptic/homomorphic sex chromosomes. Here we describe a method for determining the genetic sex in the commonly used ecotoxicological model, the fathead minnow (Pimephales promelas). Analysis of amplified fragment length polymorphisms (AFLP) in a spawn of minnows resulted in detection of 10 sex-linked AFLPs, which were isolated and sequenced. No recombination events were observed with any sex-linked AFLP in the animals examined (n=112). A polymerase chain reaction (PCR) method was then developed that determined the presence of one of these sex-linked polymorphisms for utilization in routine toxicological testing. Analyses of additional spawns from our in-house culture indicate that fathead minnows utilize a XY sex determination strategy and confirm that these markers can be used to genotype sex; however, this method is currently limited to use in laboratory studies in which breeders possess a defined genetic makeup. The genotyping method described herein can be incorporated into endocrine toxicity assays that examine the effects of chemicals on gonad differentiation.

Screening complex effluents for estrogenic activity with the T47D-KBluc cell bioassay: assay optimization and comparison with in vivo responses in fish

Wastewater treatment plant (WWTP) effluents can contain estrogenic chemicals, which potentially disrupt fish reproduction and development. The current study focused on the use of an estrogen-responsive in vitro cell bioassay (T47D-KBluc), to quantify total estrogenicity of WWTP effluents. We tested a novel sample preparation method for the T47D-KBluc assay, using powdered media prepared with direct effluent. Results of the T47D-KBluc assay were compared with the induction of estrogen receptor-regulated gene transcription in male fathead minnows (Pimephales promelas) exposed to the same effluents. Effluent samples for the paired studies were collected over the course of three months. According to the T47D-KBluc assay, the effluent estrogenicity ranged from 1.13 to 2.00 ng 17β-estradiol (E2) equivalents/L. Corresponding in vivo studies exposing male fathead minnows to 0, 10, 50, and 100% effluent dilutions demonstrated that exposure to 100% effluent significantly increased hepatic vitellogenin (VTG) and estrogen receptor α subunit transcripts relative to controls. The induction was also significant in males exposed to 250 ng E2/L or 100 ng E2/L. The in vitro and in vivo results support the conclusion that the effluent contains significant estrogenic activity, but there was a discrepancy between in vitro- and in vivo-based E2 equivalent estimates. Our results suggest that the direct effluent preparation method for the T47D-KBluc assay is a reasonable approach to estimate the estrogenicity of wastewater effluent.

Use of chemical mixtures to differentiate mechanisms of endocrine action in a small fish model

Various assays with adult fish have been developed to identify potential endocrine-disrupting chemicals (EDCs) which may cause toxicity via alterations in the hypothalamic-pituitary-gonadal (HPG) axis. These assays can be sensitive and highly diagnostic for key mechanisms such as agonism of the estrogen and androgen receptors (ERs, ARs) and inhibition of steroid synthesis. However, most of the tests do not unambiguously identify AR antagonists. The purpose of this work was to explore the utility of a mixture test design with the fathead minnow (Pimephales promelas) for detecting different classes of EDCs including AR antagonists. Adults of both sexes were exposed via the water to EDCs with diverse mechanisms of action in the absence or presence of 17beta-trenbolone (TB), a potent AR agonist which masculinizes female fathead minnows. Similar to previous studies with the model AR antagonists flutamide and vinclozolin, exposure of females to the AR antagonist cyproterone acetate in the presence of TB decreased expression of an easily-observed masculinization response, nuptial tubercle formation. Mixture studies with TB and the model ER agonists, 17alpha-ethinylestradiol and bisphenol A, also showed inhibition of tubercle formation in the females, but unlike the AR antagonists, the estrogens markedly induced synthesis of vitellogenin (VTG: egg yolk protein), particularly in males. The ER agonists also offset TB-induced depressions in plasma VTG concentrations in female fish. Additional mixture experiments were conducted with TB and triclocarban, an anti-microbial reported to enhance AR-mediated responses, or ammonia, a "negative control" with no known direct effects on HPG function. Neither chemical affected VTG status in males or females in the absence or presence of TB; however, both slightly enhanced TB-induced tubercle formation in females. Based on studies described herein and elsewhere with the fathead minnow, a TB co-exposure assay appears to be an effective approach for clearly identifying AR antagonists as well as potential EDCs with other relevant mechanisms of action.

Influence of ovarian stage on transcript profiles in fathead minnow (Pimephales promelas) ovary tissue

Interpretation of toxicogenomic experiments conducted with ovary tissue from asynchronous-spawning small fish species is complicated by background variation in the relative abundance and proportion of follicles at different stages within the ovary tissue sample. This study employed both real-time quantitative polymerase chain reaction and a 15,000 gene oligonucleotide microarray to examine variation in the fathead minnow (Pimephales promelas) ovarian transcriptional profile as a function of quantitative and qualitative differences in ovarian histology. The objectives were to provide data that could potentially aid interpretation of future toxicogenomics experiments, identify putative stage-related transcriptional markers, and generate insights into basic biological regulation of asynchronous oocyte development. Multiple lines of evidence from the present study indicate that variation in the transcriptional profile is primarily dependent on the relative abundance of previtellogenic versus vitellogenic follicles in the ovary. Due to the relatively small proportions of mature ovulated follicles or atretic follicles in the overall follicle population, few potential transcriptional markers of maturation, ovulation, or atresia could be identified. However, among the 460 differentially expressed genes identified in the present study, several targets, including HtrA serine peptidase 3 (htra3), tissue inhibitor of metalloproteinase 3 (timp3), aquaporin 8 (aqp8), transgelin 2 like (tagln2), Nedd4 family interacting protein 2 (ndfip2), chemokine ligand 12a (cxcl12a), midkine-related growth factor (mdka), and jagged 1b (jag 1b) exhibited responses and functional properties that support them as candidate molecular markers of significant shift in gross ovarian stage. Genes associated with a diversity of functions including cellular development, morphogenesis, coated vesicle transport, sexual reproduction, and neuron development, among others, were statistically enriched within the list of 460 genes differentially expressed among different ovarian classes. Overall, results of this study provide insights into background variation in ovary transcript profiles that should aid and enhance the interpretation of toxicogenomic data generated in experiments conducted with small, asynchronous-spawning fish species.