Modeling the brain-pituitary-gonad axis in salmon

PBK-TD modelling of the gonadotropic axis: Case study with two azole fungicides in female zebrafish

Endocrine disruptors (EDs) can disrupt the gonadotropic axis, which consists of the Hypothalamus-Pituitary-Gonads (HPG), notably by altering aromatase (cyp19a), a key enzyme regulating the endocrine system and reproductive function in fish. The effects of EDs can be predicted by integrating both toxicokinetic (TK) and toxicodynamic (TD) processes in order to relate adverse outcomes to external exposures. In this study, we developed a physiologically based kinetic-toxicodynamic model to simulate the disruption of the HPG axis (PBK-TD, hereafter named PBK-HPG) in female zebrafish exposed to either of two aromatase inhibitors, imazalil or prochloraz. The model was calibrated using Bayesian methods and supported by novel experimental data, including measurements of vitellogenin, 17β-estradiol, and 11-ketotestosterone levels, along with in vivo monitoring of the cyp19a1a gene in transgenic cyp19a1a-GFP ebrafish. Seamless integration of a PBK model within a TD model of the HPG-axis, provided the link between external exposure and internal levels of imazalil and prochloraz in key organs, allowing for mechanistic predictions of their inhibitory effects on gonadal aromatase. Our PBK-HPG model accurately predicted both baseline homeostasis and the effects of aromatase inhibition, with all endocrine endpoints including reproductive disruption, i.e., decreased egg production, falling within a twofold range of both experimental and literature data. Therefore, our PBK-HPG model could further support the development of a mechanistic qAOP with TK considerations. The model offers significant potential for improving environmental risk assessments of EDs and possibly other stressors across species.

In Silico Computational Transcriptomics Reveals Novel Endocrine Disruptors in Largemouth Bass ( Micropterus salmoides)

In recent years, decreases in fish populations have been attributed, in part, to the effect of environmental chemicals on ovarian development. To understand the underlying molecular events we developed a dynamic model of ovary development linking gene transcription to key physiological end points, such as gonadosomatic index (GSI), plasma levels of estradiol (E2) and vitellogenin (VTG), in largemouth bass ( Micropterus salmoides). We were able to identify specific clusters of genes, which are affected at different stages of ovarian development. A subnetwork was identified that closely linked gene expression and physiological end points and by interrogating the Comparative Toxicogenomic Database (CTD), quercetin and tretinoin (ATRA) were identified as two potential candidates that may perturb this system. Predictions were validated by investigation of reproductive associated transcripts using qPCR in ovary and in the liver of both male and female largemouth bass treated after a single injection of quercetin and tretinoin (10 and 100 μg/kg). Both compounds were found to significantly alter the expression of some of these genes. Our findings support the use of omics and online repositories for identification of novel, yet untested, compounds. This is the first study of a dynamic model that links gene expression patterns across stages of ovarian development.

Transcriptomic analysis reveals transgenerational effect of hypoxia on the neural control of testicular functions

There are over 400 hypoxic zones in the ocean worldwide. Both laboratory and field studies have shown that hypoxia causes endocrine disruption and reproductive impairments in vertebrates. More importantly, our recent study discovered that parental (F0) hypoxia exposure resulted in the transgenerational impairment of sperm quality in the F2 generation through the epigenetic regulation of germ cells. In the present study, we aim to test the hypothesis that the brain, as the major regulator of the brain-pituitary-gonad (BPG) axis, is also involved in the observed transgenerational effect. Using comparative transcriptomic analysis on brain tissues of marine medaka Oryzias melastigma, 45 common differentially expressed genes caused by parental hypoxia exposure were found in the hypoxic group of the F0 and F2 generations, and the transgenerational groups of the F2 generation. The bioinformatic analysis on this deregulated gene cluster further highlighted the possible involvement of the brain in the transgenerational effect of hypoxia on testicular structure, including abnormal morphologies of the epididymis and the seminal vesicle, and degeneration of the seminiferous tubule. This finding is concordant to the result of hematoxylin and eosin staining, which showed the reduction of testicular lobular diameter in the F0 and F2 generations. Our study demonstrated for the first time the involvement of the brain in the transgenerational effect of hypoxia.

Next-generation ecological risk assessment: Predicting risk from molecular initiation to ecosystem service delivery

Ecological risk assessment is the process of evaluating how likely it is that the environment may be impacted as the result of exposure to one or more chemicals and/or other stressors. It is not playing as large a role in environmental management decisions as it should be. A core challenge is that risk assessments often do not relate directly or transparently to protection goals. There have been exciting developments in in vitro testing and high-throughput systems that measure responses to chemicals at molecular and biochemical levels of organization, but the linkage between such responses and impacts of regulatory significance - whole organisms, populations, communities, and ecosystems - are not easily predictable. This article describes some recent developments that are directed at bridging this gap and providing more predictive models that can make robust links between what we typically measure in risk assessments and what we aim to protect.

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.

Progesterone increases ex vivo testosterone production and decreases the expression of progestin receptors and steroidogenic enzymes in the fathead minnow (Pimephales promelas) ovary

Progesterone (P4) is a metabolic precursor for a number of steroids, including estrogens and androgens. P4 also has diverse roles within the vertebrate ovary that include oocyte growth and development. The objectives of this study were to measure the effects of P4 on testosterone (T) and 17β-estradiol (E2) production in the fathead minnow (FHM) ovary and on the mRNA abundance of transcripts involved in steroidogenesis and steroid receptor signaling. Ovary explants were treated with P4 (10(-6)M) for 6 and 12h. P4 administration significantly increased T production ∼3-fold at both 6 and 12h, whereas E2 production was not affected, consistent with the hypothesis that excess P4 is not converted to terminal estrogens in the mature ovary. Nuclear progesterone receptor mRNA was decreased at 6h and membrane progesterone receptor gamma-2 mRNA was significantly down-regulated at both 6 and 12h; however there was no change in membrane progesterone receptor alpha or beta mRNA levels. Androgen receptor (ar) and estrogen receptor 2a (esr2a) mRNA were significantly reduced at 6h with P4 treatment, but there was no change in esr2b mRNA at either time point. Transcripts for enzymes in the steroid pathway (star, hsd11b2) were significantly lower at 6h compared to controls, whereas cyp17a and cyp19a mRNA abundance did not change with treatments at either time point. These data suggest that P4 incubation can lead to increased T production in the FHM ovary without a concomitant change in E2, and that the membrane bound progestin receptors are differentially regulated by P4 in the teleost ovary. As environmental progestins have received increased attention due to their suspected role as endocrine disruptors, mechanistic data on the role of exogenous P4 treatments in the male and female gonad is warranted.

High-throughput sequencing and pathway analysis reveal alteration of the pituitary transcriptome by 17α-ethynylestradiol (EE2) in female coho salmon, Oncorhynchus kisutch

Considerable research has been done on the effects of endocrine disrupting chemicals (EDCs) on reproduction and gene expression in the brain, liver and gonads of teleost fish, but information on impacts to the pituitary gland are still limited despite its central role in regulating reproduction. The aim of this study was to further our understanding of the potential effects of natural and synthetic estrogens on the brain-pituitary-gonad axis in fish by determining the effects of 17α-ethynylestradiol (EE2) on the pituitary transcriptome. We exposed sub-adult coho salmon (Oncorhynchus kisutch) to 0 or 12 ng EE2/L for up to 6 weeks and effects on the pituitary transcriptome of females were assessed using high-throughput Illumina(®) sequencing, RNA-Seq and pathway analysis. After 1 or 6 weeks, 218 and 670 contiguous sequences (contigs) respectively, were differentially expressed in pituitaries of EE2-exposed fish relative to control. Two of the most highly up- and down-regulated contigs were luteinizing hormone β subunit (241-fold and 395-fold at 1 and 6 weeks, respectively) and follicle-stimulating hormone β subunit (-3.4-fold at 6 weeks). Additional contigs related to gonadotropin synthesis and release were differentially expressed in EE2-exposed fish relative to controls. These included contigs involved in gonadotropin releasing hormone (GNRH) and transforming growth factor-β signaling. There was an over-representation of significantly affected contigs in 33 and 18 canonical pathways at 1 and 6 weeks, respectively, including circadian rhythm signaling, calcium signaling, peroxisome proliferator-activated receptor (PPAR) signaling, PPARα/retinoid x receptor α activation, and netrin signaling. Network analysis identified potential interactions between genes involved in circadian rhythm and GNRH signaling, suggesting possible effects of EE2 on timing of reproductive events.

Defining and modeling known adverse outcome pathways: Domoic acid and neuronal signaling as a case study

An adverse outcome pathway (AOP) is a sequence of key events from a molecular-level initiating event and an ensuing cascade of steps to an adverse outcome with population-level significance. To implement a predictive strategy for ecotoxicology, the multiscale nature of an AOP requires computational models to link salient processes (e.g., in chemical uptake, toxicokinetics, toxicodynamics, and population dynamics). A case study with domoic acid was used to demonstrate strategies and enable generic recommendations for developing computational models in an effort to move toward a toxicity testing paradigm focused on toxicity pathway perturbations applicable to ecological risk assessment. Domoic acid, an algal toxin with adverse effects on both wildlife and humans, is a potent agonist for kainate receptors (ionotropic glutamate receptors whose activation leads to the influx of Na(+) and Ca²(+)). Increased Ca²(+) concentrations result in neuronal excitotoxicity and cell death, primarily in the hippocampus, which produces seizures, impairs learning and memory, and alters behavior in some species. Altered neuronal Ca²(+) is a key process in domoic acid toxicity, which can be evaluated in vitro. Furthermore, results of these assays would be amenable to mechanistic modeling for identifying domoic acid concentrations and Ca²(+) perturbations that are normal, adaptive, or clearly toxic. In vitro assays with outputs amenable to measurement in exposed populations can link in vitro to in vivo conditions, and toxicokinetic information will aid in linking in vitro results to the individual organism. Development of an AOP required an iterative process with three important outcomes: a critically reviewed, stressor-specific AOP; identification of key processes suitable for evaluation with in vitro assays; and strategies for model development.

Predicting chemical impacts on vertebrate endocrine systems

Animals have evolved diverse protective mechanisms for responding to toxic chemicals of both natural and anthropogenic origin. From a governmental regulatory perspective, these protective responses complicate efforts to establish acceptable levels of chemical exposure. To explore this issue, we considered vertebrate endocrine systems as potential targets for environmental contaminants. Using the hypothalamic-pituitary-thyroid (HPT), hypothalamic-pituitary-gonad (HPG), and hypothalamic-pituitary-adrenal (HPA) axes as case examples, we identified features of these systems that allow them to accommodate and recover from chemical insults. In doing so, a distinction was made between effects on adults and those on developing organisms. This distinction was required because endocrine system disruption in early life stages may alter development of organs and organ systems, resulting in permanent changes in phenotypic expression later in life. Risk assessments of chemicals that impact highly regulated systems must consider the dynamics of these systems in relation to complex environmental exposures. A largely unanswered question is whether successful accommodation to a toxic insult exerts a fitness cost on individual animals, resulting in adverse consequences for populations. Mechanistically based mathematical models of endocrine systems provide a means for better understanding accommodation and recovery. In the short term, these models can be used to design experiments and interpret study findings. Over the long term, a set of validated models could be used to extrapolate limited in vitro and in vivo testing data to a broader range of untested chemicals, species, and exposure scenarios. With appropriate modification, Tier 2 assays developed in support of the U.S. Environmental Protection Agency's Endocrine Disruptor Screening Program could be used to assess the potential for accommodation and recovery and inform the development of mechanistically based models.