Relationship of plasma sex steroid concentrations in female fathead minnows to reproductive success and population status

Projection of Interspecific Competition (PIC) Matrices: A Conceptual Framework for Inclusion in Population Risk Assessments

Accounting for intraspecific and interspecific competition when assessing the effects of chemical and nonchemical stressors is an important uncertainty in ecological risk assessments. We developed novel projection of interspecific competition (PIC) matrices that allow for analysis of population dynamics of two or more species exposed to a given stressor(s) that compete for shared resources within a landscape. We demonstrate the application of PIC matrices to investigate the population dynamics of two hypothetical fish species that compete with one another and have differences in net reproductive rate and intrinsic rate of population increase. Population status predictions were made under scenarios that included exposure to a chemical stressor that reduced fecundity for one or both species. The results of our simulations demonstrated that measures obtained from the life table and Leslie matrix of an organism, including net reproductive rate and intrinsic rate of increase, can result in erroneous conclusions of population status and viability in the absence of a consideration of resource limitation and interspecific competition. This modeling approach can be used in conjunction with field monitoring efforts and/or laboratory testing to link effects due to stressors to possible outcomes within an ecosystem. In addition, PIC matrices could be combined with adverse outcome pathways to allow for ecosystem projection based on taxonomic conservation of molecular targets of chemicals to predict the likelihood of relative cross-species susceptibility. Overall, the present study shows how PIC matrices can integrate effects across the life cycles of multiple species, provide a linkage between endpoints observed in individual and population-level responses, and project outcomes at the community level for multiple generations for multiple species that compete for limited resources. Environ Toxicol Chem 2024;43:1406-1422. Published 2024. This article is a U.S. Government work and is in the public domain in the USA.

Polystyrene microplastics enhance the microcystin-LR-induced gonadal damage and reproductive endocrine disruption in zebrafish

The coexistence of eutrophication and plastic pollution in the aquatic environment is becoming a realistic water pollution problem worldwide. To investigate the microcystin-LR (MC-LR) bioavailability and the underlying reproductive interferences in the presence of polystyrene microplastic (PSMPs), zebrafish (Danio rerio) were exposed to individual MC-LR (0, 1, 5, and 25 μg/L) and combined MC-LR + PSMPs (100 μg/L) for 60 d. Our results showed that the existence of PSMPs increased the accumulation of MC-LR in zebrafish gonads compared to the MC-LR-only group. In the MC-LR-only exposure group, seminiferous epithelium deterioration and widened intercellular spaces were observed in the testis, and basal membrane disintegration and zona pellucida invagination were noticed in the ovary. Moreover, the existence of PSMPs exacerbated these injuries. The results of sex hormone levels showed that PSMPs enhanced MC-LR-induced reproductive toxicity, which is tightly related to the abnormal increase of 17β-estradiol (E2) and testosterone (T) levels. The changes of gnrh2, gnrh3, cyp19a1b, cyp11a, and lhr mRNA levels in the HPG axis further proved that MC-LR combined with PSMPs aggravated reproductive dysfunction. Our results revealed that PSMPs could increase the MC-LR bioaccumulation by serving as a carrier and exaggerate the MC-LR-induced gonadal damage and reproductive endocrine disruption in zebrafish.

A Multidimensional Matrix Model for Predicting the Effects of Male-Biased Sex Ratios on Fish Populations

Laboratory experiments have established that exposure to certain endocrine-active substances prior to and/or during the period of sexual differentiation can lead to skewed sex ratios in fish. However, the potential long-term population impact of biased sex ratio depends on multiple factors including the life history of the species and whether the ratio is male or female-biased. In the present study, we describe a novel multidimensional, density-dependent matrix model that analyzes age class-structure of both males and females over time, allowing for the quantitative evaluation of the effects of biased sex ratio on population status. This approach can be used in conjunction with field monitoring efforts and/or laboratory testing to link effects on sex ratio due to chemical and/or nonchemical stressors to adverse outcomes in whole organisms and populations. For demonstration purposes, we applied the model to evaluate population trajectories for fathead minnow (Pimephales promelas) exposed to prochloraz, an aromatase inhibitor, during sexual differentiation. The model also was used to explore the population impact in a more realistic exposure scenario in which both adult and early life stages of fish are exposed concurrently to prochloraz, which, in addition to altering sex ratio during development, can decrease vitellogenin and egg production in adult females. For each exposure scenario, the model was used to analyze total population size, numbers of females and of males, and sex specific recruitment of the F1 generation. The present study illustrates the utility of multidimensional matrix population models for ecological risk assessment in terms of integrating effects across a population of an organism even when chemical effects on individuals are manifested via different pathways depending on life stage. Environ Toxicol Chem 2022;41:1066-1077. Published 2022. This article is a U.S. Government work and is in the public domain in the USA.

Demonstration of an aggregated biomarker response approach to assess the impact of point and diffuse contaminant sources in feral fish in a small river case study

The assessment of the exposure of aquatic wildlife to complex environmental mixtures of chemicals originating from both point and diffuse sources and evaluating the potential impact thereof constitutes a significant step towards mitigating toxic pressure and the improvement of ecological status. In the current proof-of-concept study, we demonstrate the potential of a novel Aggregated Biomarker Response (ABR) approach involving a comprehensive set of biomarkers to identify complex exposure and impacts on wild brown trout (Salmo trutta fario). Our scenario used a small lowland river in Germany (Holtemme river in the Elbe river catchment) impacted by two wastewater treatment plants (WWTP) and diffuse agricultural runoff as a case study. The trout were collected along a pollution gradient (characterised in a parallel study) in the river. Compared to fish from the reference site upstream of the first WWTP, the trout collected downstream of the WWTPs showed a significant increase in micronucleus formation, phase I and II enzyme activities, and oxidative stress parameters in agreement with increasing exposure to various chemicals. By integrating single biomarker responses into an aggregated biomarker response, the two WWTPs' contribution to the observed toxicity could be clearly differentiated. The ABR results were supported by chemical analyses and whole transcriptome data, which revealed alterations of steroid biosynthesis and associated pathways, including an anti-androgenic effect, as some of the key drivers of the observed toxicity. Overall, this combined approach of in situ biomarker responses complemented with molecular pathway analysis allowed for a comprehensive ecotoxicological assessment of fish along the river. This study provides evidence for specific hazard potentials caused by mixtures of agricultural and WWTP derived chemicals at sublethal concentrations. Using aggregated biomarker responses combined with chemical analyses enabled an evidence-based ranking of sites with different degrees of pollution according to toxic stress and observed effects.

Simultaneous Measurement of Sex Steroid Hormones in Largemouth Bass (Micropterus salmoides) Plasma by Application of Liquid Chromatography-Tandem Mass Spectrometry

Sex steroid hormones are potential biomarkers of reproductive function in teleost fish, but their measurement continues to rely on antibody-based assays. The objective of this study was to optimize a robust and simultaneous liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for measurement of eight steroid hormones (cortisol, 11-ketotestosterone, estradiol, 17α-ethynyl estradiol, estrone, estriol, progesterone and testosterone) in fish plasma. The extraction was followed by liquid-liquid extraction with tert-butyl methyl ether and time scheduled multi-reaction monitoring (sMRM) was used for quantitation of steroids. Validation of method performance using charcoal-stripped human plasma showed extraction recoveries for eight steroids ranged from 85.5 to 108.2% with matrix effects > 80%. The limits of quantitation were 0.01 pg/μL for testosterone, 0.05 pg/μL for cortisol and progesterone, 0.1 pg/μL for 11-ketotestosterone, estradiol and estrone, 0.125 pg/μL for estriol and 0.25 pg/μL for 17α-ethynyl estradiol. The proposed method was applied to plasma samples of largemouth bass (Micropterus salmoides) collected from contaminated (Lake Apopka) and reference sites (Wildcate Lake) in Florida. Concentrations of testosterone, cortisol, estradiol and estrone were significantly different in female fish, but plasma concentration of cortisol was only statistically different in male fish between two sites (P < 0.05). This study demonstrates the application of a robust LC-MS/MS analysis for a range of sex steroid hormones representative of endocrine function in a top predator, largemouth bass.

Assessing effects of aromatase inhibition on fishes with group-synchronous oocyte development using western mosquitofish (Gambusia affinis) as a model

Exposure to certain anthropogenic chemicals can inhibit the activity to cytochrome P450 aromatase (CYP19) in fishes leading to decreased plasma 17β-estradiol (E2), plasma vitellogenin (VTG), and egg production. Reproductive dysfunction resulting from exposure to aromatase inhibitors has been extensively investigated in several laboratory model species of fish. These model species have ovaries that undergo asynchronous oocyte development, but many fishes have ovaries with group-synchronous oocyte development. Fishes with group-synchronous oocyte development have dynamic reproductive cycles which typically occur annually and are often triggered by complex environmental cues. This has resulted in a lack of test data and uncertainty regarding sensitivities to and adverse effects of aromatase inhibition. The present study used the western mosquitofish (Gambusia affinis) as a laboratory model to investigate adverse effects of chemical aromatase inhibition on group-synchronous oocyte development. Adult female western mosquitofish were exposed to either 0, 2, or 30 μg/L of the model nonsteroidal aromatase inhibiting chemical, fadrozole, for a complete reproductive cycle. Fish were sampled at four time-points representing pre-vitellogenic resting, early vitellogenesis, late vitellogenesis/early ovarian recrudescence, and late ovarian recrudescence. Temporal changes in numerous reproductive parameters were measured, including gonadosomatic index (GSI), plasma sex steroids, and expression of selected genes in the brain, liver, and gonad that are important for reproduction. In contrast to fish from the control treatment, fish exposed to 2 and 30 μg/L of fadrozole had persistent elevated expression of cyp19 in the ovary, depressed expression of vtg in the liver, and a low GSI. These responses suggest that completion of a group-synchronous reproductive cycle was unsuccessful during the assay in fish from either fadrozole treatment. These adverse effects data show that exposure to aromatase inhibitors has the potential to cause reproductive dysfunction in a wide range of fishes with both asynchronous and group-synchronous reproductive strategies.

A multidimensional density dependent matrix population model for assessing risk of stressors to fish populations

Modeling exposure and recovery of fish and wildlife populations after stressor mitigation serves as a basis for evaluating remediation success. Herein, we develop a novel multidimensional density dependent matrix population model that analyzes both size-structure and age class-structure simultaneously. This modeling approach emphasizes application in conjunction with field monitoring efforts (e.g., through effects-based monitoring programs) and/or laboratory analysis to link effects due to stressors to outcomes in populations. We applied the model to investigate Atlantic killifish (Fundulus heteroclitus) exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin with effects on fertility and survival rates. The Atlantic killifish is an important and well-studied model organism for understanding the effects of pollutants and other stressors in estuarine and marine ecosystems. For each exposure concentration, the corresponding plots of total population size, population size structure, and age structure over time were generated. The present study serves as an example of how a multidimensional matrix population model can integrate effects across the life cycle, provide a linkage between endpoints observed in the individual and ecological risk to the population as a whole, and project outcomes for multiple generations.

Differential Sensitivity to In Vitro Inhibition of Cytochrome P450 Aromatase (CYP19) Activity Among 18 Freshwater Fishes

There is significant concern regarding potential impairment of fish reproduction associated with endocrine disrupting chemicals. Aromatase (CYP19) is a steroidogenic enzyme involved in the conversion of androgens to estrogens. Inhibition of aromatase by chemicals can result in reduced concentrations of estrogens leading to adverse reproductive effects. These effects have been extensively investigated in a small number of laboratory model fishes, such as fathead minnow (Pimephales promelas), Japanese medaka (Oryzias latipes), and zebrafish (Danio rerio). But, differences in sensitivity among species are largely unknown. Therefore, this study took a first step toward understanding potential differences in sensitivity to aromatase inhibitors among fishes. Specifically, a standard in vitro aromatase inhibition assay using subcellular fractions of whole tissue homogenates was used to evaluate the potential sensitivity of 18 phylogenetically diverse species of freshwater fish to the nonsteroidal aromatase inhibitor fadrozole. Sensitivity to fadrozole ranged by more than 52-fold among these species. Five species were further investigated for sensitivity to up to 4 additional nonsteroidal aromatase inhibitors, letrozole, imazalil, prochloraz, and propiconazole. Potencies of each of these chemicals relative to fadrozole ranged by up to 2 orders of magnitude among the 5 species. Fathead minnow, Japanese medaka, and zebrafish were among the least sensitive to all the investigated chemicals; therefore, ecological risks of aromatase inhibitors derived from these species might not be adequately protective of more sensitive native fishes. This information could guide more objective ecological risk assessments of native fishes to chemicals that inhibit aromatase.

Stochastically modeling multiscale stationary biological processes

Large scale biological responses are inherently uncertain, in part as a consequence of noisy systems that do not respond deterministically to perturbations and measurement errors inherent to technological limitations. As a result, they are computationally difficult to model and current approaches are notoriously slow and computationally intensive (multiscale stochastic models), fail to capture the effects of noise across a system (chemical kinetic models), or fail to provide sufficient biological fidelity because of broad simplifying assumptions (stochastic differential equations). We use a new approach to modeling multiscale stationary biological processes that embraces the noise found in experimental data to provide estimates of the parameter uncertainties and the potential mis-specification of models. Our approach models the mean stationary response at each biological level given a particular expected response relationship, capturing variation around this mean using conditional Monte Carlo sampling that is statistically consistent with training data. A conditional probability distribution associated with a biological response can be reconstructed using this method for a subset of input values, which overcomes the parameter identification problem. Our approach could be applied in addition to dynamical modeling methods (see above) to predict uncertain biological responses over experimental time scales. To illustrate this point, we apply the approach to a test case in which we model the variation associated with measurements at multiple scales of organization across a reproduction-related Adverse Outcome Pathway described for teleosts.

Quantitative Response-Response Relationships Linking Aromatase Inhibition to Decreased Fecundity are Conserved Across Three Fishes with Asynchronous Oocyte Development

Quantitative adverse outcome pathways (qAOPs) describe quantitative response-response relationships that can predict the probability or severity of an adverse outcome for a given magnitude of chemical interaction with a molecular initiating event. However, the taxonomic domain of applicability for these predictions is largely untested. The present study began defining this applicability for a previously described qAOP for aromatase inhibition leading to decreased fecundity developed using data from fathead minnow (Pimephales promelas). This qAOP includes quantitative response-response relationships describing plasma 17β-estradiol (E2) as a function of plasma fadrozole, plasma vitellogenin (VTG) as a function of plasma E2, and fecundity as a function of plasma VTG. These quantitative response-response relationships simulated plasma E2, plasma VTG, and fecundity measured in female zebrafish (Danio rerio) exposed to fadrozole for 21 days but not these responses measured in female Japanese medaka (Oryzias latipes). However, Japanese medaka had different basal levels of plasma E2, plasma VTG, and fecundity. Normalizing basal levels of each measurement to equal those of female fathead minnow enabled the relationships to accurately simulate plasma E2, plasma VTG, and fecundity measured in female Japanese medaka. This suggests that these quantitative response-response relationships are conserved across these three fishes when considering relative change rather than absolute measurements. The present study represents an early step toward defining the appropriate taxonomic domain of applicability and extending the regulatory applications of this qAOP.

Expression of sf1 and dax-1 are regulated by thyroid hormones and androgens during Silurana tropicalis early development

Thyroid hormones (THs) and androgens have been shown to be extensively involved in sexual development; however, relatively little is known with regard to TH-related and androgenic actions in sex determination. We first established expression profiles of three sex-determining genes (sf1, dax-1, and sox9) during the embryonic development of Western clawed frogs (Silurana tropicalis). Transcripts of sf1 and sox9 were detected in embryos before the period in which embryonic transcription commences indicating maternal transfer, whereas dax-1 transcripts were not detected until later in development. To examine whether TH status affects sex-determining gene expression in embryonic S. tropicalis, embryos were exposed to co-treatments of iopanoic acid (IOP), thyroxine (T4), or triiodothyronine (T3) for 96 h. Expression profiles of TH receptors and deiodinases reflect inhibition of peripheral deiodinase activity by IOP and recovery by T3. Relevantly, elevated TH levels significantly increased the expression of sf1 and dax-1 in embryonic S. tropicalis. Further supporting TH-mediated regulation, examination of the presence and frequency of transcription factor binding sites in the putative promoter regions of sex-determining genes in S. tropicalis and rodent and fish models using in silico analysis also identified TH motifs in the putative promoter regions of sf1 and dax-1. Together these findings advocate that TH actions as early as the period of embryogenesis may affect gonadal fate in frogs. Mechanisms of TH and androgenic crosstalk in relation to the regulation of steroid-related gene expression were also investigated.

Assessing recovery of in vitro steroid production in male rainbow darter (Etheostoma caeruleum) in response to municipal wastewater treatment plant infrastructure changes

The present study examined in vitro 11-ketotestosterone and testosterone production by the testes of rainbow darter (Etheostoma caeruleum) collected from selected reference sites and downstream of 2 municipal wastewater treatment plants (MWWTPs; Waterloo and Kitchener) on the central Grand River (Ontario, Canada), over a 6-yr period (2011-2016). The main objective was to investigate if infrastructure upgrades at the Kitchener MWWTP in 2012 resulted in a recovery of this response in the post-upgrade period (2013-2016). Two supporting studies showed that the fall season is appropriate for measuring in vitro sex steroid production because it provides stable detection of steroid patterns, and that the sample handling practiced in the present study did not introduce a bias. Infrastructure upgrades of the Kitchener MWWTP resulted in significant reductions in ammonia and estrogenicity. After the upgrades, 11-ketotestosterone production by MWWTP-exposed fish increased in 2013 and it continued to recover throughout the study period of 2014 through 2016, returning to levels measured in reference fish. Testosterone production was less sensitive and it lacked consistency. The Waterloo MWWTP underwent some minor upgrades but the level of ammonia and estrogenicity remained variable over time. The production of 11-ketotestosterone and testosterone in rainbow darter below the Waterloo MWWTP was variable and without a clear recovery pattern over the course of the present study. The results of the present study demonstrated that measuring production of sex steroids (especially 11-ketotestosterone) over multiple years can be relevant for assessing responses in fish to environmental changes such as those resulting from major infrastructure upgrades. Environ Toxicol Chem 2018;37:501-514. © 2017 SETAC.

Applying adverse outcome pathways and species sensitivity-weighted distribution to predicted-no-effect concentration derivation and quantitative ecological risk assessment for bisphenol A and 4-nonylphenol in aquatic environments: A case study on Tianjin City, China

Adverse outcome pathways (AOPs) are a novel concept that effectively considers the toxic modes of action and guides the ecological risk assessment of chemicals. To better use toxicity data including biochemical or molecular responses and mechanistic data, we further developed a species sensitivity-weighted distribution (SSWD) method for bisphenol A and 4-nonylphenol. Their aquatic predicted-no-effect concentrations (PNECs) were derived using the log-normal statistical extrapolation method. We calculated aquatic PNECs of bisphenol A and 4-nonylphenol with values of 4.01 and 0.721 µg/L, respectively. The ecological risk of each chemical in different aquatic environments near Tianjin, China, a coastal municipality along the Bohai Sea, was characterized by hazard quotient and probabilistic risk quotient assessment techniques. Hazard quotients of 7.02 and 5.99 at 2 municipal sewage sites using all of the endpoints were observed for 4-nonylphenol, which indicated high ecological risks posed by 4-nonylphenol to aquatic organisms, especially endocrine-disrupting effects. Moreover, a high ecological risk of 4-nonylphenol was indicated based on the probabilistic risk quotient method. The present results show that combining the SSWD method and the AOP concept could better protect aquatic organisms from adverse effects such as endocrine disruption and could decrease uncertainty in ecological risk assessment. Environ Toxicol Chem 2018;37:551-562. © 2017 SETAC.

Quantitative Adverse Outcome Pathways and Their Application to Predictive Toxicology

A quantitative adverse outcome pathway (qAOP) consists of one or more biologically based, computational models describing key event relationships linking a molecular initiating event (MIE) to an adverse outcome. A qAOP provides quantitative, dose-response, and time-course predictions that can support regulatory decision-making. Herein we describe several facets of qAOPs, including (a) motivation for development, (b) technical considerations, (c) evaluation of confidence, and (d) potential applications. The qAOP used as an illustrative example for these points describes the linkage between inhibition of cytochrome P450 19A aromatase (the MIE) and population-level decreases in the fathead minnow (FHM; Pimephales promelas). The qAOP consists of three linked computational models for the following: (a) the hypothalamic-pitutitary-gonadal axis in female FHMs, where aromatase inhibition decreases the conversion of testosterone to 17β-estradiol (E2), thereby reducing E2-dependent vitellogenin (VTG; egg yolk protein precursor) synthesis, (b) VTG-dependent egg development and spawning (fecundity), and (c) fecundity-dependent population trajectory. While development of the example qAOP was based on experiments with FHMs exposed to the aromatase inhibitor fadrozole, we also show how a toxic equivalence (TEQ) calculation allows use of the qAOP to predict effects of another, untested aromatase inhibitor, iprodione. While qAOP development can be resource-intensive, the quantitative predictions obtained, and TEQ-based application to multiple chemicals, may be sufficient to justify the cost for some applications in regulatory decision-making.

Molecular initiating events of the intersex phenotype: Low-dose exposure to 17α-ethinylestradiol rapidly regulates molecular networks associated with gonad differentiation in the adult fathead minnow testis

Intersex, or the presence of oocytes in the testes, has been documented in fish following exposure to wastewater effluent and estrogenic compounds. However, the molecular networks underlying the intersex condition are not completely known. To address this, we exposed male fathead minnows to a low, environmentally-relevant concentration of 17alpha-ethinylestradiol (EE2) (15ng/L) and measured the transcriptome response in the testis after 96h to identify early molecular initiating events that may proceed the intersex condition. The short-term exposure to EE2 did not affect gonadosomatic index and proportion of gametes within the testes. However, the production of 11-ketotestosterone and testosterone from the testis in vitro was decreased relative to controls. Expression profiling using a 8×60K fathead minnow microarray identified 10 transcripts that were differentially expressed in the testes, the most dramatic change being that of coagulation factor XIII A chain (20-fold increase). Transcripts that included guanine nucleotide binding protein (Beta Polypeptide 2), peroxisome proliferator-activated receptor delta, and WNK lysine deficient protein kinase 1a, were down-regulated by EE2. Subnetwork enrichment analysis revealed that EE2 suppressed transcriptional networks associated with steroid metabolism, hormone biosynthesis, and sperm mobility. Most interesting was that gene networks associated with doublesex and mab-3 related transcription factor 1 (dmrt1) were suppressed in the adult testis, despite the fact that dmrt1 itself was not different in expression from control males. Transcriptional networks involving forkhead box L2 (foxl2) (transcript involved in ovarian follicle development) were increased in expression in the testis. Noteworthy was that a gene network associated to granulosa cell development was increased over 100%, suggesting that this transcriptome network may be important for monitoring estrogenic exposures. Other cell processes rapidly downregulated by EE2 at the transcript level included glucose homeostasis, response to heavy metal, amino acid catabolism, and the cyclooxygenase pathway. Conversely, lymphocyte chemotaxis, intermediate filament polymerization, glucocorticoid metabolism, carbohydrate utilization, and anterior/posterior axis specification were increased. These data provide new insight into the transcriptional responses that are perturbed prior to gonadal remodeling and intersex following exposure to estrogens. These data demonstrate that low concentrations of EE2 (1) rapidly suppresses male hormone production, (2) down-regulate molecular networks related to male sex differentiation, and (3) induce transcriptional networks related to granulosa cell development in the adult testis. These responses are hypothesized to be key molecular initiating events that occur prior to the development of the intersex phenotype following estrogenic exposures.

Multimedia screening of contaminants of emerging concern (CECS) in coastal urban watersheds in southern California (USA)

To examine the occurrence and fate of contaminants of emerging concern (CECs) and inform future monitoring of CECs in coastal urban waterways, water, sediment, and fish tissue samples were collected and analyzed for a broad suite of pharmaceuticals and personal care products (PPCPs), commercial and/or household chemicals, current use pesticides, and hormones in an effluent-dominated river and multiple embayments in southern California (USA). In the Santa Clara River, which receives treated wastewater from several facilities, aqueous phase CECs were detectable at stations nearest discharges from municipal wastewater treatment plants but were attenuated downstream. Sucralose and the chlorinated phosphate flame retardants tris(1-chloro-2-propyl) phosphate (TCPP), tris(1,3-dichloro-2-propyl) phosphate (TDCPP), and tris(2-chloroethyl) phosphate (TCEP) were most abundant in water, with maximum concentrations of 35 μg/L, 3.3 μg/L, 1.4 μg/L, and 0.81 μg/L, respectively. Triclocarban, an antimicrobial agent in use for decades, was more prevalent in water than triclosan or nonylphenol. Maximum concentrations of bifenthrin, permethrin, polybrominated diphenyl ethers (PBDEs), and degradates of fipronil exceeded CEC-specific monitoring trigger levels recently established for freshwater and estuarine sediments by factors of 10 to 1000, respectively. Maximum fish tissue concentrations of PBDEs varied widely (370 ng/g and 7.0 ng/g for the Santa Clara River and coastal embayments, respectively), with most species exhibiting concentrations at the lower end of this range. These results suggest that continued monitoring of pyrethroids, PBDEs, and degradates of fipronil in sediment is warranted in these systems. In contrast, aqueous pharmaceutical concentrations in the Santa Clara River were not close to exceeding current monitoring trigger levels, suggesting a lower priority for targeted monitoring in this medium. Environ Toxicol Chem 2016;35:1986-1994. © 2016 SETAC.

Vitellogenin, sex steroid levels and gonadal biomarkers in wild Solea solea and Solea senegalensis from NW Mediterranean fishing grounds

Specimens of Solea solea and Solea senegalenesis at different developmental stages were obtained from seven fishing grounds along the NW Mediterranean. Gonad development in males was classified into five stages, from early spermatogenesis to recovery, while four stages were considered in females, from growth to maturation. Vitellogenin (VTG) and sex steroid levels including an estrogen (estradiol, E2), two androgens (testosterone, T and 11-ketotestosterone, 11KT) and a progestin (17,20β-dihydroxy pregn-4-en-3-one, 17,20β-P or maturation inducing steroid, MIS) were analysed in plasma. Their levels were more clearly related to the developmental stage of the gonads than to the sampling site characteristics. In addition, enzyme activities in gonads, such as acetylcholinesterase (AChE) and carboxylesterase (CbE) were gender-dependent and higher in males than in females. Gonadal glutathione S-transferase (GST) activity was enhanced in the most anthropogenic impacted sites. VTG was absent in males and very low or undetectable in immature females, while mature females exhibited high VTG levels, clearly related to the gonado-somatic index. Sex steroid levels (ng/ml) varied in males and females regardless of the species. E2 levels in females ranged from 0.22 to 6.98 while in males ranged from 0.11 to 0.27. T varied from 0.12 to 0.93 in females and from 0.56 to 1.36 in males, while 11KT in females fluctuated from 0.03 to 0.57 and from 0.26 to 6.42 in males. Similarly, MIS in females ranged from 0.75 to 3.71 and from 1.12 to 5.61 in males. The lack of endocrine disturbances was confirmed by histological examination of the gonads. This study informs on basal sex hormone levels and enzyme activities during gonadal maturation of wild Solea spp. that can be useful in the identification and further remediation of possible pollution events.

A Computational Model of the Rainbow Trout Hypothalamus-Pituitary-Ovary-Liver Axis

Reproduction in fishes and other vertebrates represents the timely coordination of many endocrine factors that culminate in the production of mature, viable gametes. In recent years there has been rapid growth in understanding fish reproductive biology, which has been motivated in part by recognition of the potential effects that climate change, habitat destruction and contaminant exposure can have on natural and cultured fish populations. New approaches to understanding the impacts of these stressors are being developed that require a systems biology approach with more biologically accurate and detailed mathematical models. We have developed a multi-scale mathematical model of the female rainbow trout hypothalamus-pituitary-ovary-liver axis to use as a tool to help understand the functioning of the system and for extrapolation of laboratory findings of stressor impacts on specific components of the axis. The model describes the essential endocrine components of the female rainbow trout reproductive axis. The model also describes the stage specific growth of maturing oocytes within the ovary and permits the presence of sub-populations of oocytes at different stages of development. Model formulation and parametrization was largely based on previously published in vivo and in vitro data in rainbow trout and new data on the synthesis of gonadotropins in the pituitary. Model predictions were validated against several previously published data sets for annual changes in gonadotropins and estradiol in rainbow trout. Estimates of select model parameters can be obtained from in vitro assays using either quantitative (direct estimation of rate constants) or qualitative (relative change from control values) approaches. This is an important aspect of mathematical models as in vitro, cell-based assays are expected to provide the bulk of experimental data for future risk assessments and will require quantitative physiological models to extrapolate across biological scales.

Reproduction in fishes and other vertebrates represents the timely coordination of many endocrine factors that culminate in the production of mature, viable gametes. Improving the ability to estimate reproductive performance in fish is important, due to the growth of the aquaculture industry and the need to maintain adequate broodstock and concerns over the effects of anthropogenic stressors on feral fish populations. We present here a quantitative, mathematical model of the female rainbow trout reproductive cycle. We show how the model is able to accurately describe experimentally measured data associated with pituitary, ovarian and liver reproductive performance. We also use the model to describe similar data sets collected in rainbow trout by other researchers. An important value of quantitative biological models is the ability to simulate various physiological conditions, real or hypothetical. We demonstrate this by predicting the effects of exposure to an endocrine disruptor on oocyte growth. The need to limit cost and animal usage will encourage future experimental studies to use in vitro methods. The model presented here can assist with the extrapolation of in vitro effects to the whole fish.

Reproduction impairment and endocrine disruption in female zebrafish after long-term exposure to MC-LR: A life cycle assessment

Microcystin-LR (MC-LR) has been found to cause reproductive and developmental impairments as well as to disrupt sex hormone homeostasis of fish during acute and sub-chronic toxic experiments. However, fish in natural environments are continuously exposed to MC-LR throughout their entire life cycle as opposed to short-term exposure. Here, we tested the hypothesis that the mechanism by which MC-LR harms female fish reproduction and development within natural water bodies is through interference of the reproductive endocrine system. In the present study, zebrafish hatchlings (5 d post-fertilization) were exposed to 0, 0.3, 3 and 30 μg/L MC-LR for 90 d until reaching sexual maturity. Female zebrafish were selected, and the changes in growth and developmental indicators, ovarian ultrastructure as well as the levels of gonadal steroid hormones and vitellogenin (VTG) were examined along with the transcription of related genes in the hypothalamic-pituitary-gonadal-liver axis (HPGL-axis). The results showed for the first time, a life cycle exposure to MC-LR caused growth inhibition, decreased ovary weight and ovarian ultra-pathological lesions. Decreased ovarian testosterone levels indicated that MC-LR disrupted sex steroid hormone balance. Significantly up-regulated transcription of brain FSHβ and LHβ along with ovarian ERα, FSHR and LHR suggested positive feedback regulation in the HPGL-axis was induced as a compensatory mechanism for MC-LR damage. It was also noted that ovarian VTG content and hepatic ERα and VTG1 expression were all down-regulated, which might be responsible for reduced vitellus storage noted in our histological observations. Our findings indicate that a life cycle exposure to MC-LR impairs the development and reproduction of female zebrafish by disrupting the transcription of related HPGL-axis genes, suggesting that MC-LR has potential adverse effects on fish reproduction and thus population dynamics in MCs-contaminated aquatic environment.

Developmental abnormalities and differential expression of genes induced in oil and dispersant exposed Menidia beryllina embryos

Exposure of fish embryos to relatively low concentrations of oil has been implicated in sub-lethal toxicity. The objective of this study was to determine the effects of the exposure of Menidia beryllina embryos at 30-48h post-fertilization to the water accommodated fractions of oil (WAF, 200ppm, v/v), dispersants (20ppm, v/v, Corexit 9500 or 9527), and mixtures of oil and each of the dispersants to produce chemically enhanced water accommodated fractions (CEWAFs) over a 72-hour period. The polyaromatic hydrocarbon (PAH) and benzene, toluene, ethylene and xylene (BTEX) constituents of the 5X concentrated exposure solutions (control, WAF, dispersants and CEWAFs) were determined and those of the 1× exposures were derived using a dilution factor. PAH, BTEX and low molecular weight PAH constituents greater than 1ppb were observed in WAF and the dispersants, but at much higher levels in CEWAFs. The WAF and CEWAFs post-weathering were diluted at 1:5 (200ml WAF/CEWAF: 800ml 25ppt saltwater) for embryo exposures. Mortality, heartbeat, embryo normalcy, abnormality types and severities were recorded. The qPCR assay was used to quantify abundances of transcripts of target genes for sexual differentiation and sex determination (StAR, dmrt-1, amh, cyp19b, vtg and chg-L,), growth regulation (ghr) and stress response (cyp1a and Hsp90); and gapdh served as the housekeeping gene. Temperature was 21±1.5°C throughout the experimental period, while mortality was low and not significantly different (p=0.68) among treatments. Heartbeat was significantly different (0.0034) with the lowest heartbeats recorded in Corexit 9500 (67.5beats/min) and 9527 (67.1beats/min) exposed embryos compared with controls (82.7beats/min). Significantly more treated embryos were in a state of deterioration, with significantly more embryos presenting arrested tissue differentiation compared with controls (p=0.021). Exposure to WAF, dispersants and CEWAF induced aberrant expression of all the genes, with star, dmrt-1, ghr and hsp90 being significantly down-regulated in CEWAF and cyp19b in Corexit 9527. The cyp1a and cyp19b were significantly up-regulated in CEWAFs and WAF, respectively. The molecular endpoints were most sensitive, especially the expression of star, cyp19b, cyp1a, hsp90 and could therefore be used as early indicators of long term effects of Corexit 9500 and 9527 usage in oil spill management on M. beryllina, a valid sentinel for oil pollution events.

The effects of model androgen 5α-dihydrotestosterone on mummichog (Fundulus heteroclitus) reproduction under different salinities

Endocrine disrupting substances (EDSs) have the potential to disturb sensitive hormone pathways, particularly those involved in development and reproduction. Both fresh and estuarine water bodies receive inputs of EDSs from a variety of sources, including sewage effluent, industrial effluent and agricultural runoff. Based on current literature, freshwater species appear to respond to lower levels of EDSs than estuarine or marine species. Therefore, effects elicited by EDSs in freshwater teleosts may not be an accurate representation of how EDSs affect teleosts in estuarine and marine environments. To address this potential difference, a short-term reproductive bioassay was conducted under conditions of low and high salinity using mummichog (Fundulus heteroclitus), a euryhaline species that is native to the east coast of North America. The goals of this study were to determine the response of mummichog when exposed to an androgenic EDS and whether salinity affected the response. A model androgen, 5α-dihydrotestosterone (DHT), was selected for this experiment. Impacts on reproduction were evaluated at multiple biological levels, including physiological (sex steroid levels), organismal (gonad size and gonad morphology), and functional (egg production) endpoints. Under conditions of high salinity, egg production was significantly reduced at all exposure concentrations. Under conditions of low salinity, there were no significant differences based on DHT treatment; however, egg production in all treatment groups including the control were significantly reduced relative to the high salinity control group. Other reproductive endpoints, such as sex steroid production, showed stronger correlation to fecundity in females than males. This study demonstrates that mummichog fecundity is sensitive to androgenic endocrine disruption while also underscoring the importance of how changes in salinity, an environmental variable, can impact reproduction.

Linking mechanistic toxicology to population models in forecasting recovery from chemical stress: A case study from Jackfish Bay, Ontario, Canada

Recovery of fish and wildlife populations after stressor mitigation serves as a basis for evaluating remediation success. Unfortunately, effectively monitoring population status on a routine basis can be difficult and costly. In the present study, the authors describe a framework that can be applied in conjunction with field monitoring efforts (e.g., through effects-based monitoring programs) to link chemically induced alterations in molecular and biochemical endpoints to adverse outcomes in whole organisms and populations. The approach employs a simple density-dependent logistic matrix model linked to adverse outcome pathways (AOPs) for reproductive effects in fish. Application of this framework requires a life table for the organism of interest, a measure of carrying capacity for the population of interest, and estimation of the effect of stressors on vital rates of organisms within the study population. The authors demonstrate the framework using linked AOPs and population models parameterized with long-term monitoring data for white sucker (Catostomus commersoni) collected from a study site at Jackfish Bay, Lake Superior, Canada. Individual responses of fish exposed to pulp mill effluent were used to demonstrate the framework's capability to project alterations in population status, both in terms of ongoing impact and subsequent recovery after stressor mitigation associated with process changes at the mill. The general approach demonstrated at the Jackfish Bay site can be applied to characterize population statuses of other species at a variety of impacted sites and can account for effects of multiple stressors (both chemical and nonchemical) and dynamics within complex landscapes (i.e., meta-populations including emigration and immigration processes).

Effects of 17 α-methyltestosterone on transcriptome, gonadal histology and sex steroid hormones in rare minnow Gobiocypris rarus

The 17α-methyltestosterone (MT), a synthetic androgen, is known for its interference effects on the endocrine system. Aiming to investigate the transcriptome profiling of gonads induced by MT and to understand the molecular mechanism by which MT causes adverse effects in fish, transcriptome profiling of gonads, gonadal histology and the sex steroid hormones in response to MT were analyzed in Gobiocypris rarus. Eight libraries, 4 from the ovary and 4 from the testis, were constructed and sequenced and then a total number of clean reads per sample ranging from 7.03 to 9.99 million were obtained. In females, a total of 191 transcripts were differentially regulated by MT, consisting of 102 up-regulated transcripts and 89 down-regulated transcripts. In males, 268 differentially expressed genes with 108 up-regulated and 160 down-regulated were detected upon MT exposure. Testosterone serves as the major sex steroid hormone content in G. rarus of both sexes. The concentrations of 17β-estradiol, testosterone and 11-ketotestosterone were significantly increased in females and decreased in males after MT exposure. Interestingly, MT caused a decreased number of vitellogenic oocytes in the ovary and spermatozoa in the testis. After MT exposure, four differentially expressed genes (ndufa4, slc1a3a, caskin-2 and rpt3) were found in G. rarus of both sexes. Overall, we suggest that MT seemed to affect genes involved in pathways related to physiological processes in the gonads of G. rarus. These processes include the electron transfer of Complex IV, endothelial cell activation, axon growth and guidance, and proteasome assembly and glutamate transport metabolic.

Linking feeding inhibition with reproductive impairment in Gammarus confirms the ecological relevance of feeding assays in environmental monitoring

The in situ feeding bioassay in Gammarus fossarum is recognized as a reliable tool for monitoring the toxicity of freshwater contamination. However, whether recorded feeding inhibitions can potentially provoke population-level adverse outcomes remains an open question. In the present study, the authors present an experimental study in G. fossarum, which contributes to the quantitative description of the links between feeding inhibitions and impacts on female reproductive performance. The authors studied the impacts of food deprivation on reproductive endpoints (i.e., fecundity, fertility, molt cycle) during 2 successive molting cycles. Among the main results, the authors found that food deprivation triggered a slowdown of the molting process and a reduction in fertility but no alteration to embryonic development. These reproductive impairments appeared for feeding inhibition values usually recorded in monitoring programs of environmental pollution. Using a population model translating Gammarus life-history, the authors predicted that the observed reproductive alterations predict a strong degradation of population dynamics. The present study underlines the importance of feeding inhibition in population-level risk assessment and discusses how establishing upscaling schemes based on quantitative mechanistic links between impacts at different levels of biological organization can be applied in environmental monitoring to propose an ecotoxicological assessment of water quality, which would be sensitive, specific, and ecologically relevant.

Which coastal and marine environmental contaminants are truly emerging?

To better understand the past and present impact of contaminants of emerging concern (CECs) in coastal and marine ecosystems, archived samples were analyzed for a broad suite of analytes, including pharmaceuticals and personal care products (PPCPs), flame retardants (including PBDEs), perfluorinated compounds (PFCs), and current-use pesticides. Surface sediment, mussels (Mytilus spp.) and sediment core samples collected from the California (USA) coast were obtained from environmental specimen banks. Selected CECs were detected in recent surface sediments, with nonylphenol (4-NP), its mono- and di-ethoxylates (NP1EO and NP2EO), triclocarban, and pyrethroid insecticides in the greatest abundance. Alkylphenols, triclocarban, and triclosan were present in sediment core segments from the 1970s, as well as in Mytilus tissue collected during the 1990s. Increasing concentrations of some CECs (e.g., miconazole, triclosan) were observed in the surface layers (ca. 2007) of a sediment core, in contrast to peak concentrations of 4-NP and triclocarban corresponding to input during the 1970s, and an apparent peak input for PBDEs during the 1990s. These results suggest that chemicals sometimes referred to as "emerging" (e.g., alkylphenols, triclocarban) have been present in the aquatic environment for several decades and are decreasing in concentration, whereas others (e.g., miconazole, triclosan) are increasing.

Low concentrations of 17β-trenbolone induce female-to-male reversal and mortality in the frog Pelophylax nigromaculatus

Trenbolone, as a growth promoter in animal agriculture, has become an environmental androgen in surface water. Here, we aimed to reveal the effects of 17β-trenbolone on survival, growth, and gonadal differentiation in the frog Pelophylax nigromaculatus, which is widespread in East Asia and undergoing population decline. P. nigromaculatus tadpoles were exposed to 17β-trenbolone (0.1, 1, 10 μg/L) from Gosner stage 24/25 to complete metamorphosis. We found that 17β-trenbolone resulted in significantly high mortality in a concentration-dependent manner, with a decrease in body weight in the high concentration group compared with the solvent control. Based on gross gonadal morphology, no females were observed, instead of about 15% ambiguous sexes and 85% males, in all 17β-trenbolone treatment groups. Like normal testes, the gonads with sex-ambiguous morphology exhibited testicular histology, showing that the sex-ambiguous gonads were incomplete ovary-to-testis reversals (IOTTRs) with certain ovarian morphological features. In the IOTTRs, the transcriptional levels of ovary-biased genes decreased drastically relative to normal ovaries, and even declined to the levels in normal testes. These observations confirmed that all test concentrations of 17β-trenbolone resulted in 100% sex reversal, although some sex-reversed testes retained some ovarian characteristics at the morphological level. To our knowledge, this is the first report strongly demonstrating that trenbolone can cause female-to-male reversal in amphibians. Given that the lowest concentration tested is environmentally relevant, our study highlights the risks of trenbolone and other environmental androgens for P. nigromaculatus and other amphibians, in particular the species with high sensitivity of gonadal differentiation to androgenic chemicals.

Benzo[a]pyrene effects on reproductive endpoints in Fundulus heteroclitus

Benzo[a]pyrene (BaP) is a polycyclic aromatic hydrocarbon (PAH) that has been implicated in modulating aromatase enzyme function with the potential to interrupt normal reproductive function. The aim of this study was to use a fish model, Fundulus heteroclitus, to assess whether BaP exposure adversely impacts reproduction. Adult fish were exposed to waterborne BaP nominal concentrations of (0, 1, or 10 μg/l) for 28 days. Males and females were combined for the second half of the exposure (days 14-28) in order to quantitate egg production and fertilization success. Egg fertilization and subsequent hatching success of F1 embryos was significantly decreased by the high dose of BaP. In males, both gonad weight and plasma testosterone concentrations were significantly reduced compared to controls by 10 μg/l BaP. Histopathological examination of testes including spermatogonia, spermatocyte and spermatid cyst areas, percentage of cysts per phase, and area of spermatozoa per seminiferous tubule were not significantly affected. Other biomarkers, including male liver weight, liver vitellogenin (vtg) mRNA expression and sperm concentrations, were also not affected. In females, estradiol concentrations were significantly reduced after BaP exposure, but egg production, gonad weight, liver weight, vtg expression and oocyte maturation were not altered. Steroid concentrations in Fundulus larvae from exposed parents at 1 and 3 weeks posthatch were not significantly changed. BaP exposure at these environmentally relevant concentrations caused negative alterations particularly in male fish to both biochemical and phenotypic biomarkers associated with reproduction and multigenerational embryo survival.

An inter-laboratory study on the variability in measured concentrations of 17β-estradiol, testosterone, and 11-ketotestosterone in white sucker: implications and recommendations

Endocrine-disrupting chemicals are exogenous substances that can impact the reproduction of fish, potentially by altering circulating concentrations of 17β-estradiol (E2), testosterone (T), and 11-ketotestosterone (11-KT). Common methods to measure steroids in plasma samples include radioimmunoassays (RIAs) and enzyme-linked immunosorbant assays (ELISAs). The present study examines variability in E2, T, and 11-KT across 8 laboratories measuring reference and pulp mill effluent-exposed white sucker (Catostomus commersoni) plasma. We examine the contribution of assay type (RIA vs ELISA), standardized hormone extraction, location of values on the standard curve (upper and lower limits), and other variables on the ability to distinguish hormone levels between reference and exposed fish and the impact of these variables on quantitation of hormones in different laboratories. Of the 8 participating laboratories, 7 of 8 and 7 of 7 identified differences between sites for female E2 and female T, respectively, and 7 of 7 and 4 of 5 identified no differences between male T and male 11-KT. Notably, however, the ng/mL concentration of steroids measured across laboratories varied by factors of 10-, 6-, 14-, and 10-fold, respectively. Within laboratory intra-assay variability was generally acceptable and below 15%. Factors contributing to interlaboratory variability included calculation errors, assay type, and methodology. Based on the interlaboratory variability detected, we provide guidelines and recommendations to improve the accuracy and precision of steroid measurements in fish ecotoxicology studies.

Acetylcholinesterase activity in the host-parasite system of the cod Gadus morhua and acanthocephalan Echinorhynchus gadi from the southern Baltic Sea

Acetylcholinesterase (AChE) activity measurement is widely used as a specific biomarker of neurotoxic effects. The aim of this study was to evaluate AChE activity in a host fish (the cod) and its acanthocephalan parasite Echinorhynchus gadi from the southern Baltic. AChE activity in hosts and parasites was inversely related: the highest cod AChE activity corresponded to the lowest E. gadi enzymatic activity and vice versa ("mirror effect"). This is the first report on the simultaneous application of this biomarker in cod and its acanthocephalan parasites. Results obtained for the host-parasite system are complementary and provide comprehensive information about the response of this biomarker. Analysis of the system allows for detection of a greater number of factors influencing AChE activity in the marine environment than separate analysis of the host and parasites. Thus, AChE activity measurement in a host-parasite system may be considered to be a promising tool for biomonitoring.

The occurrence and fate of chemicals of emerging concern in coastal urban rivers receiving discharge of treated municipal wastewater effluent

To inform future monitoring and assessment of chemicals of emerging concern (CECs) in coastal urban watersheds, the occurrence and fate of more than 60 pharmaceuticals and personal care products (PPCPs), commercial/household chemicals, current-use pesticides, and hormones were characterized in 2 effluent-dominated rivers in southern California (USA). Water samples were collected during 2 low-flow events at locations above and below the discharge points of water reclamation plants (WRPs) and analyzed using gas chromatography-mass spectrometry and liquid chromatography-tandem mass spectrometry. Approximately 50% of targeted CECs were detectable at stations downstream from WRPs, compared with <31% and <10% at the reference stations above the WRPs. Concentrations of chlorinated phosphate flame retardants were highest among the CECs tested, with mean total aggregate concentrations of tris(2-chloroethyl) phosphate (TCEP), tris(1-chloro-2-propyl) phosphate (TCPP), and tris(1,3-dichloro-2-propyl) phosphate (TDCPP) of 3400 ng/L and 2400 ng/L for the 2 rivers. Maximum in-stream concentrations of pyrethroids (bifenthrin and permethrin), diclofenac, and galaxolide exceeded risk-based thresholds established for monitoring of CECs in effluent-dominated receiving waters. In contrast, maximum concentrations of PPCPs commonly detected in treated wastewater (e.g., acetaminophen, N,N,diethyl-meta-toluamide [DEET], and gemfibrozil) were less than 10% of established thresholds. Attenuation of target CECs was not observed downstream of WRP discharge until dilution by seawater occurred in the tidal zone, partly because of the short hydraulic residence times in these highly channelized systems (<3 d). In addition to confirming CECs for future in-stream monitoring, these results suggest that conservative mass transport is an important boundary condition for assessment of the input, fate, and effects of CECs in estuaries at the bottom of these watersheds.

Investigating alternatives to the fish early-life stage test: a strategy for discovering and annotating adverse outcome pathways for early fish development

The fish early-life stage (FELS) test (Organisation for Economic Co-operation and Development [OECD] test guideline 210) is the primary test used internationally to estimate chronic fish toxicity in support of ecological risk assessments and chemical management programs. As part of an ongoing effort to develop efficient and cost-effective alternatives to the FELS test, there is a need to identify and describe potential adverse outcome pathways (AOPs) relevant to FELS toxicity. To support this endeavor, the authors outline and illustrate an overall strategy for the discovery and annotation of FELS AOPs. Key events represented by major developmental landmarks were organized into a preliminary conceptual model of fish development. Using swim bladder inflation as an example, a weight-of-evidence-based approach was used to support linkage of key molecular initiating events to adverse phenotypic outcomes and reduced young-of-year survival. Based on an iterative approach, the feasibility of using key events as the foundation for expanding a network of plausible linkages and AOP knowledge was explored and, in the process, important knowledge gaps were identified. Given the scope and scale of the task, prioritization of AOP development was recommended and key research objectives were defined relative to factors such as current animal-use restrictions in the European Union and increased demands for fish toxicity data in chemical management programs globally. The example and strategy described are intended to guide collective efforts to define FELS-related AOPs and develop resource-efficient predictive assays that address the toxicological domain of the OECD 210 test.

An adaptive, comprehensive monitoring strategy for chemicals of emerging concern (CECs) in California's Aquatic Ecosystems

A scientific advisory panel was convened by the State of California to recommend monitoring for chemicals of emerging concern (CECs) in aquatic systems that receive discharge of municipal wastewater treatment plant (WWTP) effluent and stormwater runoff. The panel developed a risk-based screening framework that considered environmental sources and fate of CECs observed in receiving waters across the State. Using existing occurrence and risk threshold data in water, sediment, and biological tissue, the panel applied the framework to identify a priority list of CECs for initial monitoring in three representative receiving water scenarios. The initial screening list of 16 CECs identified by the panel included consumer and commercial chemicals, flame retardants, pesticides, pharmaceuticals and personal care products, and natural hormones. The panel designed an iterative, phased strategy with interpretive guidelines that direct and update management actions commensurate with potential risk identified using the risk-based framework and monitoring data. Because of the ever-changing nature of chemical use, technology, and management practices, the panel offered recommendations to improve CEC monitoring, including development of bioanalytical screening methods whose responses integrate exposure to complex mixtures and that can be linked to higher-order effects; development or refinement of models that predict the input, fate, and effects of future chemicals; and filling of key data gaps on CEC occurrence and toxicity. Finally, the panel stressed the need for adaptive management, allowing for future review of, and if warranted, modifications to the strategy to incorporate the latest science available to the water resources community.

Selection of reliable biomarkers from PCR array analyses using relative distance computational model: methodology and proof-of-concept study

It is increasingly evident about the difficulty to monitor chemical exposure through biomarkers as almost all the biomarkers so far proposed are not specific for any individual chemical. In this proof-of-concept study, adult male zebrafish (Danio rerio) were exposed to 5 or 25 µg/L 17β-estradiol (E2), 100 µg/L lindane, 5 nM 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) or 15 mg/L arsenic for 96 h, and the expression profiles of 59 genes involved in 7 pathways plus 2 well characterized biomarker genes, vtg1 (vitellogenin1) and cyp1a1 (cytochrome P450 1A1), were examined. Relative distance (RD) computational model was developed to screen favorable genes and generate appropriate gene sets for the differentiation of chemicals/concentrations selected. Our results demonstrated that the known biomarker genes were not always good candidates for the differentiation of pair of chemicals/concentrations, and other genes had higher potentials in some cases. Furthermore, the differentiation of 5 chemicals/concentrations examined were attainable using expression data of various gene sets, and the best combination was the set consisting of 50 genes; however, as few as two genes (e.g. vtg1 and hspa5 [heat shock protein 5]) were sufficient to differentiate the five chemical/concentration groups in the present test. These observations suggest that multi-parameter arrays should be more reliable for biomonitoring of chemical exposure than traditional biomarkers, and the RD computational model provides an effective tool for the selection of parameters and generation of parameter sets.

From 'omics to otoliths: responses of an estuarine fish to endocrine disrupting compounds across biological scales

Endocrine disrupting chemicals (EDCs) cause physiological abnormalities and population decline in fishes. However, few studies have linked environmental EDC exposures with responses at multiple tiers of the biological hierarchy, including population-level effects. To this end, we undertook a four-tiered investigation in the impacted San Francisco Bay estuary with the Mississippi silverside (Menidia audens), a small pelagic fish. This approach demonstrated links between different EDC sources and fish responses at different levels of biological organization. First we determined that water from a study site primarily impacted by ranch run-off had only estrogenic activity in vitro, while water sampled from a site receiving a combination of urban, limited ranch run-off, and treated wastewater effluent had both estrogenic and androgenic activity. Secondly, at the molecular level we found that fish had higher mRNA levels for estrogen-responsive genes at the site where only estrogenic activity was detected but relatively lower expression levels where both estrogenic and androgenic EDCs were detected. Thirdly, at the organism level, males at the site exposed to both estrogens and androgens had significantly lower mean gonadal somatic indices, significantly higher incidence of severe testicular necrosis and altered somatic growth relative to the site where only estrogens were detected. Finally, at the population level, the sex ratio was significantly skewed towards males at the site with measured androgenic and estrogenic activity. Our results suggest that mixtures of androgenic and estrogenic EDCs have antagonistic and potentially additive effects depending on the biological scale being assessed, and that mixtures containing androgens and estrogens may produce unexpected effects. In summary, evaluating EDC response at multiple tiers is necessary to determine the source of disruption (lowest scale, i.e. cell line) and what the ecological impact will be (largest scale, i.e. sex ratio).

Assessment of status of white sucker (Catostomus commersoni) populations exposed to bleached kraft pulp mill effluent

Credible ecological risk assessments often need to include analysis of population-level impacts. In the present study, a predictive model was developed to investigate population dynamics for white sucker (Catostomus commersoni) exposed to pulp mill effluent at a well-studied site in Jackfish Bay, Lake Superior, Canada. The model uniquely combines a Leslie population projection matrix and the logistic equation to translate changes in the fecundity and the age structure of a breeding population of white sucker exposed to pulp mill effluent to alterations in population growth rate. Application of this density-dependent population projection model requires construction of a life table for the organism of interest, a measure of carrying capacity, and an estimation of the effect of stressors on vital rates. A white sucker population existing at carrying capacity and subsequently exposed to pulp mill effluent equivalent to a documented exposure experienced during the period 1988 to 1994 in Jackfish Bay would be expected to exhibit a 34% to 51% annual decrease in recruitment during the first 5 yr of exposure and approach a population size of 71% of carrying capacity. The Jackfish Bay study site contains monitoring data for biochemical endpoints in white sucker, including circulating sex steroid concentrations, that could be combined with population modeling to utilize the model demonstrated at the Jackfish Bay study site for investigation of other white sucker populations at sites that are less data-rich.

Inhibition of vitellogenin gene induction by 2,3,7,8-tetrachlorodibenzo-p-dioxin is mediated by aryl hydrocarbon receptor 2 (AHR2) in zebrafish (Danio rerio)

Vitellogenins are hepatically derived yolk-protein precursors required for oogenesis in all oviparous teleosts. Altered gene-regulation of vitellogenesis by environmental contaminants can have profound effects on reproductive success, and ultimately population sustainability. To better understand chemical effects on vitellogenin gene regulation, we tested the hypothesis that activation of the aryl hydrocarbon receptor 2 (AHR2) by dioxin inhibits the estrogen receptor pathway regulation of 3 vitellogenin genes (vtg1-3) in vivo, using zebrafish (Danio rerio) as a model teleost. Using an embryo-larval bioassay, embryos were either treated with 1000 pptr (parts-per-trillion, pg/mL) 17α-ethynylestradiol (EE2) alone from 6h post fertilization (hpf) to 4 days post fertilization (dpf), or pre-treated with dioxin (4-5 hpf) prior to EE2. Pre-treatment with 400 pptr 2,3,7,8-tetrachlorodibenzo-p-dioxin (2,3,7,8-TCDD) or 1,2,3,7,8-pentachlorodibenzo-p-dioxin inhibited the EE2 induction of vtg1, vtg2 and vtg3 by >95% (p≤0.05). In comparison, a splice-blocking AHR2 morpholino used to down-regulate ahr2 expression significantly reduced the inhibition of vtg1, vtg2 and vtg3 by 400 pptr 2,3,7,8-TCDD (20.7-27.4% rescue). These studies demonstrate that 2,3,7,8-TCDD directly inhibits the vitellogenin pathway in vivo through activation of the AHR2. This work provides evidence for AHR2 dependent cross-talk inhibition of vitellogenin genes and offers insight into anti-estrogenic reproductive effects observed in oviparous species exposed to AHR agonist contaminants.

Differences in the induction of cyp1A and related genes in cultured rainbow trout Oncorhynchus mykiss. Additional considerations for the use of EROD activity as a biomarker

Two rainbow trout Oncorhynchus mykiss fish farms were repeatedly sampled in order to observe the variability of ethoxyresorufin-O-deethylase (EROD) activity and of related genes in the liver. Fish coming from fish farm A exhibited EROD levels that could be considered as basal according to the scientific literature, however, EROD activity in fish coming from fish farm B was significantly increased. This was accompanied by augmented aryl hydrocarbon receptor (ahr) and cytochrome P4501A (cyp1A) messenger RNA expression and reduced oestrogen receptor (er) and vitellogenin (vtg) transcription. Only sediment extracts from the entry channel of fish farm B induced EROD activity in O. mykiss cultured cells, however, this induction could not be explained by the levels of polyaromatic hydrocarbons (PAH) and polychlorinated biphenyls (PCB) measured in the sediments. The results of this study point out that O. mykiss cultured in fish farms could be used as sentinels for indication of pollution. In this particular work, however, no conclusive evidence has been found for a relationship between the presence of PAHs and PCBs and the observed EROD induction.

Evaluation of the fish short term reproduction assay for detecting endocrine disrupters

In a fish testing strategy, positive results of the fish short term reproduction assay (FSTR), often trigger a definitive test like the fish sexual development test (FSDT) or the fish full life cycle test (FFLC), entailing ethical and economic problems. This study analysed 137 studies encompassing 35 chemicals with different modes of actions (MOAs). Variability is quantified for MOA endpoints vitellogenin (VTG) and secondary sex characteristics (SSCs) as well as for apical endpoints. Two MOA endpoints could indicate estrogenic, anti-estrogenic, androgenic, anti-androgenic and steroidogenesis activities. Great variability, however, has been observed for chemicals with anti-androgenic and steroidogenesis activities, suggesting that TG229/230 may not be sensitive enough to detect these types of chemicals and may produce false negatives. Changes in apical endpoints like fecundity are not limited to endocrine disrupting chemicals (EDCs). Non-EDCs could induce the similar effects on these apical endpoints. If elucidating MOA is needed, targeted in vitro MOA tests are suggested. Positive in vitro MOA results trigger a definitive test, which could be used for confirmation of the MOA in vivo and for deriving a no observed effect concentration (NOEC). Based on positive MOA results of TG229, a definitive test such as the FSDT or the FFLC is still needed, because the current TG229 has limitation on the derivation of a NOEC. An extended TG229 with more power to detect reproduction effects, as recently proposed in the OECD test guideline program, would improve the possibility to derive a NOEC and increase its usefulness in risk assessment.

Decreased vitellogenin inducibility and 17β-estradiol levels correlated with reduced egg production in killifish (Fundulus heteroclitus) from Newark Bay, NJ

Aquatic species inhabiting polluted estuaries are exposed to complex mixtures of xenobiotics which can alter normal reproduction. We previously reported that female Atlantic killifish (Fundulus heteroclitus) from the highly contaminated Newark Bay, NJ (USA) exhibited an inhibition of oocyte development due to reduced vitellogenin (egg-yolk precursor) levels. Our hypothesis was that the inhibition of oocyte development in Newark Bay killifish is due to (1) deficient levels of circulating 17β-estradiol, and (2) a decreased sensitivity of the vitellogenin pathway to physiological doses of 17β-estradiol. In the first study, adult naïve killifish from Tuckerton, NJ (reference) were caged at Tuckerton and Newark Bay. After 1 month, males caged at Newark Bay exhibited inductions of hepatic vitellogenin and estrogen receptor α, which were transient and returned to basal levels after 2 months (p≤0.05). In the second study, fecundity and 17β-estradiol levels were measured in reproductively active adult females from Tuckerton and Newark Bay. Tuckerton females produced 140 eggs per female and Newark Bay females produced 11 eggs per female. Embryos from Newark Bay had 34% greater mortality and 28% less hatch, relative to Tuckerton. In addition, embryo mass and yolk-volume of Newark Bay embryos compared to Tuckerton embryos was 16% and 25% lower, respectively. Circulating 17β-estradiol levels in Newark Bay females (0.26 ng/mL) were measured to be 8-fold lower than Tuckerton females (2.25 ng/mL). In the third study, adult killifish from both sites were dosed with 17β-estradiol to assess the sensitivity of the vitellogenin pathway. At doses of 0.01, 0.1, 1 and 10 ng/g body weight, induction levels of circulating vitellogenin in Newark Bay males were significantly inhibited by 97, 99, 98 and 44%, respectively, compared to Tuckerton males. At doses of 0.01, 0.1, 1, 10 and 100 ng/g body weight, induction levels of circulating vitellogenin in Newark Bay females were inhibited by 89, 79, 61, 40 and 30%, respectively, compared to Tuckerton females. These differences in inducibility could not be explained by altered hepatic expression of estrogen receptors α, βa or βb. Based on the caged and dose-response studies, contaminants that down-regulate vitellogenin would interfere with its ability to be used as a biomarker for xeno-estrogen exposures. These studies demonstrate that contaminants within Newark Bay exert both estrogenic and anti-estrogenic responses which results in an overtly anti-estrogenic phenotype (reduced egg production due to inhibition of vitellogenesis).

Endpoint sensitivity in fish endocrine disruption assays: regulatory implications

Identifying potential endocrine disrupting chemicals (EDCs) needs screening and testing for mode of action (MOA) and intrinsic toxicological properties. MOA is often indicated by biomarker endpoints, whereas toxicity by apical endpoints. Risk assessment is mainly based on apical but not on biomarker endpoints. The 21-day fish assay (OECD TG229) is considered a screening test. But it includes both biomarker and apical endpoints. This study explores the utility of results of the 21-day fish assay for risk assessment purposes. Endpoint sensitivity was analysed by compiling 142 data sets for 21-day fish assays and 38 data sets for the fish sexual development test (FSDT), encompassing 62 chemicals with different MOAs. Conclusions from this analysis include: (1) vitellogenin (VTG), fecundity and gonad histology are the most sensitive endpoints for fathead minnow, medaka and zebrafish in 21-day fish assays; secondary sex characteristics (SSC) are a less sensitive endpoint and is likely inadequate to detect all known MOAs. (2) Biomarker endpoints like VTG and apical endpoints like fecundity from the 21-day fish assay can be used for risk assessment. (3) Lowest observed effect concentrations (LOECs) of the most chemicals are comparable for the 21-day fish assay and for the FSDT, further supporting that results of 21-day fish assays can be used for risk assessment. However, a significant difference in LOECs was observed for some chemicals, suggesting that chemical specific effects should be taken into account. This paper emphasizes that a weight of evidence approach is important for interpretation of results of the 21-day fish assay.

Systems biology: leading the revolution in ecotoxicology

The rapid development of new technologies such as transcriptomics, proteomics, and metabolomics (Omics) are changing the way ecotoxicology is practiced. The data deluge has begun with genomes of over 65 different aquatic species that are currently being sequenced, and many times that number with at least some level of transcriptome sequencing. Integrating these top-down methodologies is an essential task in the field of systems biology. Systems biology is a biology-based interdisciplinary field that focuses on complex interactions in biological systems, with the intent to model and discover emergent properties of the system. Recent studies demonstrate that Omics technologies provide valuable insight into ecotoxicity, both in laboratory exposures with model organisms and with animals exposed in the field. However, these approaches require a context of the whole animal and population to be relevant. Powerful approaches using reverse engineering to determine interacting networks of genes, proteins, or biochemical reactions are uncovering unique responses to toxicants. Modeling efforts in aquatic animals are evolving to interrelate the interacting networks of a system and the flow of information linking these elements. Just as is happening in medicine, systems biology approaches that allow the integration of many different scales of interaction and information are already driving a revolution in understanding the impacts of pollutants on aquatic systems.

Pharmaceuticals as neuroendocrine disruptors: lessons learned from fish on Prozac

Pharmaceuticals are increasingly detected in a variety of aquatic systems. One of the most prevalent environmental pharmaceuticals in North America and Europe is the antidepressant fluoxetine, a selective serotonin reuptake inhibitor (SSRI) and the active ingredient of Prozac. Usually detected in the range below 1 μg/L, fluoxetine and its active metabolite norfluoxetine are found to bioaccumulate in wild-caught fish, particularly in the brain. This has raised concerns over potential disruptive effects of neuroendocrine function in teleost fish, because of the known role of serotonin (5-HT) in the modulation of diverse physiological processes such as reproduction, food intake and growth, stress and multiple behaviors. This review describes the evolutionary conservation of the 5-HT transporter (the therapeutic target of SSRIs) and reviews the disruptive effects of fluoxetine on several physiological endpoints, including involvement of neuroendocrine mechanisms. Studies on the goldfish, Carassius auratus, whose neuroendocrine regulation of reproduction and food intake are well characterized, are described and represent a reliable model to study neuroendocrine disruption. In addition, fish studies investigating the effects of fluoxetine, not only on reproduction and food intake, but also on stress and behavior, are discussed to complement the emerging picture of neuroendocrine disruption of physiological systems in fish exposed to fluoxetine. Environmental relevance and key lessons learned from the effects of the antidepressant fluoxetine on fish are highlighted and may be helpful in designing targeted approaches for future risk assessments of pharmaceuticals disrupting the neuroendocrine system in general.

Reverse engineering adverse outcome pathways

The toxicological effects of many stressors are mediated through unknown, or incompletely characterized, mechanisms of action. The application of reverse engineering complex interaction networks from high dimensional omics data (gene, protein, metabolic, signaling) can be used to overcome these limitations. This approach was used to characterize adverse outcome pathways (AOPs) for chemicals that disrupt the hypothalamus-pituitary-gonadal endocrine axis in fathead minnows (FHM, Pimephales promelas). Gene expression changes in FHM ovaries in response to seven different chemicals, over different times, doses, and in vivo versus in vitro conditions, were captured in a large data set of 868 arrays. Potential AOPs of the antiandrogen flutamide were examined using two mutual information-based methods to infer gene regulatory networks and potential AOPs. Representative networks from these studies were used to predict network paths from stressor to adverse outcome as candidate AOPs. The relationship of individual chemicals to an adverse outcome can be determined by following perturbations through the network in response to chemical treatment, thus leading to the nodes associated with the adverse outcome. Identification of candidate pathways allows for formation of testable hypotheses about key biological processes, biomarkers, or alternative endpoints that can be used to monitor an AOP. Finally, the unique challenges facing the application of this approach in ecotoxicology were identified and a road map for the utilization of these tools presented.

Challenges and opportunities with the use of biomarkers to predict reproductive impairment in fishes exposed to endocrine disrupting substances

Biomarkers are commonly used as signposts to evaluate the potential of contaminants to disrupt the endocrine system. However, the relationship between responses in these biomarkers and whole organism endpoints that directly affect population status is not clearly understood. In this study, the relationship between egg production (a whole-organism endpoint which has been directly linked to population-level responses) and biomarkers (sex steroids, vitellogenin (VTG) and gonad size) is examined. Data were collected from short-term reproductive tests in which a wide variety of fish species were exposed to a suite of contaminants with known or unknown modes/mechanisms of action (MOA). The potential to use biomarkers as signposts was evaluated by determining the occurrence of false negatives (i.e., an effect in egg production was not accompanied by a biomarker response) and false positives (i.e., an effect in biomarkers was not followed by an effect in egg production). The quantitative relationships between biomarkers and egg production, and the ability to use these quantitative relationships to predict population-level responses based on modeling was also assessed. A suite of female biomarkers resulted in a relatively low occurrence of both false positives and negatives, indicating the potential for their use as signposts for reproductive effects via endocrine disruption. Egg production in short-term adult fish reproductive tests showed significant relationships to 17β-estradiol (E2), changes in female VTG levels, and relative female gonad size (gonadosomatic index; GSI). Weaker significant relationships were found between egg production and both VTG levels and GSI in males. However, use of these quantitative relationships to predict population-level effects are cautioned because of high levels of uncertainty. This study demonstrates that there are qualitative and quantitative relationships among biomarkers, regardless of fish species used or the MOA of contaminants and concludes that a suite of female reproductive biomarkers can be used as effective signposts to screen chemicals and assess waste streams for endocrine disrupting substances with different MOA.

Environmentally relevant exposure to 17alpha-ethinylestradiol affects the telencephalic proteome of male fathead minnows

Estrogens are key mediators of neuronal processes in vertebrates. As such, xenoestrogens present in the environment have the potential to alter normal central nervous system (CNS) function. The objectives of the present study were (1) to identify proteins with altered abundance in the male fathead minnow telencephalon as a result of low-level exposure to 17alpha-ethinylestradiol (EE(2)), and (2) to better understand the underlying mechanisms of 17beta-estradiol (E(2)) feedback in this important neuroendocrine tissue. Male fathead minnows exposed to a measured concentration of 5.4 ng EE(2)/L for 48 h showed decreased plasma E(2) levels of approximately 2-fold. Of 77 proteins that were quantified statistically, 14 proteins were down-regulated after EE(2) exposure, including four histone proteins, ATP synthase, H+ transporting subunits, and metabolic proteins (lactate dehydrogenase B4, malate dehydrogenase 1b). Twelve proteins were significantly induced by EE(2) including microtubule-associated protein tau (Mapt), astrocytic phosphoprotein, ependymin precursor, and calmodulin. Mapt showed an increase in protein abundance but a decrease in mRNA expression after EE(2) exposure(,) suggesting there may be a negative feedback response in the telencephalon to decreased mRNA transcription with increasing Mapt protein abundance. These results demonstrate that a low, environmentally relevant exposure to EE(2) can rapidly alter the abundance of proteins involved in cell differentiation and proliferation, neuron network morphology, and long-term synaptic potentiation. Together, these findings provide a better understanding of the molecular responses underlying E(2) feedback in the brain and demonstrate that quantitative proteomics can be successfully used in ecotoxicology to characterize affected cellular pathways and endocrine physiology.

Effects of monocrotophos on the reproductive axis in the female goldfish (Carassius auratus)

Monocrotophos (MCP) is a highly toxic organophosphorus pesticide. To elucidate the influence of MCP on female fish reproduction, plasma 17beta-estradiol, testosterone and gonadotropin levels and aromatase and gonadotropin beta subunit gene expression levels were examined in female goldfish (Carassius auratus) following a 21-day exposure to 0.01, 0.10 and 1.00 mg L-(1) 40% MCP-based pesticide in a semi-static exposure system. The results indicated that MCP induced increases in plasma 17beta-estradiol levels and the 17beta-estradiol/testosterone ratio via interference with the reproductive axis at multiple potential sites through two mechanisms: (a) MCP exposure enhanced the mRNA expression of gonadal aromatase, the enzyme that converts androgens into estrogens, consequently reducing plasma testosterone levels and increasing plasma concentrations of 17beta-estradiol; and (b) MCP treatment increased follicle-stimulating hormone beta subunit mRNA expression and secretion and decreased luteinizing hormone beta subunit mRNA expression and secretion, leading to the disruption of reproductive endocrine control and androgen and estrogen balance. This study provided convincing evidence for reproductive toxicology of MCP by disrupting of the HPG axis at multiple sites in female goldfish.