Gene expression patterns in rainbow trout, Oncorhynchus mykiss, exposed to a suite of model toxicants

Biochemical and molecular biomarkers and their association with anthropogenic chemicals in wintering Manx shearwaters (Puffinus puffinus)

Anthropogenic pollution poses a threat to marine conservation by causing chronic toxic effects. Seabirds have contact throughout their lives with pollutants like plastic, metals, polychlorinated biphenyls (PCBs), and organochlorine pesticides such as hexachlorocyclohexanes (HCHs). We assessed 155 Manx shearwaters (Puffinus puffinus) stranded along the Brazilian coast, analyzing associations between organic pollutants, plastic ingestion, biomarkers (transcript levels of aryl hydrocarbon receptor, cytochrome P450-1A-5 [CYP1A5], UDP-glucuronosyl-transferase [UGT1], estrogen receptor alpha-1 [ESR1], and heat shock protein-70 genes) and enzymes activity (ethoxy-resorufin O-deethylase and glutathione S-transferase [GST]). Plastic debris was found in 29 % of the birds. The transcription of UGT1 and CYP1A5 was significantly associated with hexachlorobenzene (HCB) and PCBs levels. ESR1 was associated with HCB and Mirex, and GST was associated with Drins and Mirex. While organic pollutants affected shearwaters more than plastic ingestion, reducing plastic availability remains relevant as xenobiotics are also potentially adsorbed onto plastics.

Contaminants of emerging concern in the Maumee River and their effects on freshwater mussel physiology

Contaminants of emerging concern pose a serious hazard to aquatic wildlife, especially freshwater mussels. The growing number of contaminants in aquatic systems requires scientists and managers to prioritize contaminants that are most likely to elicit a biological response for further monitoring and toxicological testing. The objectives of this study were to identify a sub-category of contaminants most likely to affect Pyganodon grandis and to describe alterations in metabolites and gene expression between various sites. Mussels were deployed in cages for two weeks at four sites along the Maumee River Basin, Ohio, USA. Water samples were analyzed for the presence of 220 contaminants. Hemolymph samples were collected for metabolomics and analyzed using mass spectrometry. Contaminants that significantly covaried with metabolites were identified using partial least-squares (PLS) regression. Tissue samples were collected for transcriptomics, RNA was sequenced using an Illumina HiSeq 2500, and differential expression analysis was performed on assembled transcripts. Of the 220 targeted contaminants, 69 were detected in at least one water sample. Of the 186 metabolites detected in mussel hemolymph, 43 showed significant differences between the four sites. The PLS model identified 44 contaminants that significantly covaried with changes in metabolites. A total of 296 transcripts were differentially expressed between two or more sites, 107 received BLAST hits, and 52 were annotated and assigned to one or more Gene Ontology domains. Our analyses reveal the contaminants that significantly covaried with changes in metabolites and are most likely to negatively impact freshwater mussel health and contribute to ongoing population declines in this group of highly endangered animals. Our integration of "omics" technologies provides a broad and in-depth assessment of the short-term effects of contaminants on organismal physiology. Our findings highlight which contaminants are most likely to be causing these changes and should be prioritized for more extensive toxicological testing.

Effects of heavy metals on fish physiology - A review

The pollution by heavy metals poses a serious threat to the aquatic environment and to the organisms if the concentration of heavy metals in the environment exceeds the safe limits. Due to their non-biodegradable and long persistence nature in the environment, heavy metals cause toxicity in fish by producing oxygen reactive species through oxidizing radical production. In this review, we investigated the effects of heavy metals on fish physiology with special emphasis on hemato-biochemical properties, immunological parameters especially hormones and enzymes, histopathology of different major organs and underlying molecular mechanisms. All those parameters are significantly affected by heavy metal exposure and are found to be important bio-monitoring tools to assess heavy metal toxicity. Hematological and biochemical alterations have been documented including cellular and nuclear abnormalities in different fish species exposed to different concentrations of heavy metals. Major fish organs (gills, liver, kidneys) including intestine, muscles showed different types of pathology specific to organs in acute and chronic exposure to different heavy metals. This study also revealed the expression of different genes involved in oxidative stress and detoxification of heavy metals. In a nutshell, this article shades light on the manipulation of fish physiology by the heavy metals and sought attention in the prevention and maintenance of aquatic environments particularly from heavy metals contaminations.

Therapeutic Intervention with Dietary Chitosan Nanoparticles Alleviates Fish Pathological and Molecular Systemic Inflammatory Responses against Infections

Marine bio-sourced chitosan nanoparticles (CSNP) are antimicrobial and immunomodulatory agents beneficial for fish medicine. Herein, dietary CSNP was investigated for the amelioration of the systemic inflammatory responses of an induced fish model. One hundred and forty-four rainbow trout were assigned to one pathogen-free and non-supplemented group (negative control), and three challenged groups: non-supplemented (positive control), CSNP-preventive, and CSNP-therapeutic. After a feeding experiment extended for 21 days, the organosomatic indices (OSI) and molecular aspects were assessed. After a challenge experiment extended for further 28 days, CSNP-therapeutic intervention was assessed on fish survival and systemic inflammatory responses on pathology, histo-morphology, and molecular aspects. With CSNP administration, OSI nonsignificantly decreased and the relative expression of targeted inflammatory-mediator genes was significantly increased. The CSNP-therapeutic fish showed an RPS of 80% as compared to the positive control group, and CSNP-therapeutic administration retained the highest gene expression augmentation up to 28 days after the challenge. Notably, the splenic reticulin fibers framework of the CSNP-therapeutic group retained the highest integrity among the groups during the infection. After recovery, reticulin fibers density in the CSNP-therapeutic samples was significantly higher than in the negative control group, which indicates high innate immunity. Thus, CSNP showed promising biotherapeutic features enhancing fish resistance against infections.

Gill and Liver Transcript Expression Changes Associated With Gill Damage in Atlantic Salmon (Salmo salar)

Gill damage represents a significant challenge in the teleost fish aquaculture industry globally, due to the gill’s involvement in several vital functions and direct contact with the surrounding environment. To examine the local and systemic effects accompanying gill damage (which is likely to negatively affect gill function) of Atlantic salmon, we performed a field sampling to collect gill and liver tissue after several environmental insults (e.g., harmful algal blooms). Before sampling, gills were visually inspected and gill damage was scored; gill scores were assigned from pristine [gill score 0 (GS0)] to severely damaged gills (GS3). Using a 44K salmonid microarray platform, we aimed to compare the transcriptomes of pristine and moderately damaged (i.e., GS2) gill tissue. Rank Products analysis (5% percentage of false-positives) identified 254 and 34 upregulated and downregulated probes, respectively, in GS2 compared with GS0. Differentially expressed probes represented genes associated with functions including gill remodeling, wound healing, and stress and immune responses. We performed gill and liver qPCR for all four gill damage scores using microarray-identified and other damage-associated biomarker genes. Transcripts related to wound healing (e.g., neb and klhl41b) were significantly upregulated in GS2 compared with GS0 in the gills. Also, transcripts associated with immune and stress-relevant pathways were dysregulated (e.g., downregulation of snaclec 1-like and upregulation of igkv3) in GS2 compared with GS0 gills. The livers of salmon with moderate gill damage (i.e., GS2) showed significant upregulation of transcripts related to wound healing (i.e., chtop), apoptosis (e.g., bnip3l), blood coagulation (e.g., f2 and serpind1b), transcription regulation (i.e., pparg), and stress-responses (e.g., cyp3a27) compared with livers of GS0 fish. We performed principal component analysis (PCA) using transcript levels for gill and liver separately. The gill PCA showed that PC1 significantly separated GS2 from all other gill scores. The genes contributing most to this separation were pgam2, des, neb, tnnt2, and myom1. The liver PCA showed that PC1 significantly separated GS2 from GS0; levels of hsp70, cyp3a27, pparg, chtop, and serpind1b were the highest contributors to this separation. Also, hepatic acute phase biomarkers (e.g., serpind1b and f2) were positively correlated to each other and to gill damage. Gill damage-responsive biomarker genes and associated qPCR assays arising from this study will be valuable in future research aimed at developing therapeutic diets to improve farmed salmon welfare.

Exposure to an androgenic agricultural pollutant does not alter metabolic rate, behaviour, or morphology of tadpoles

Globally, amphibian species are experiencing dramatic population declines, and many face the risk of imminent extinction. Endocrine-disrupting chemicals (EDCs) have been recognised as an underappreciated factor contributing to global amphibian declines. In this regard, the use of hormonal growth promotants in the livestock industry provides a direct pathway for EDCs to enter the environment-including the potent anabolic steroid 17β-trenbolone. Emerging evidence suggests that 17β-trenbolone can impact traits related to metabolism, somatic growth, and behaviour in non-target species. However, far less is known about possible effects of 17β-trenbolone on anuran species, particularly during early life stages. Accordingly, in the present study we investigated the effects of 28-day exposure to 17β-trenbolone (mean measured concentrations: 10 and 66 ng/L) on body size, body condition, metabolic rate, and anxiety-related behaviour of tadpoles (Limnodynastes tasmaniensis). Specifically, we measured rates of O₂ consumption of individual tadpoles as a proxy for metabolic rate and quantified their swimming activity and their time spent in the upper half of the water column as indicators of anxiety-related behaviour. Counter to our predictions based on effects observed in other taxa, we detected no effect of 17β-trenbolone on body size, metabolic rate, or behaviour of tadpoles; although, we did detect a subtle, but statistically significant decrease in body condition at the highest 17β-trenbolone concentration. We hypothesise that 17β-trenbolone may induce taxa-specific effects on metabolic function, growth, and anxiety-related behaviour, with anurans being less sensitive to disruption than fish, and encourage further cross-taxa investigation to test this hypothesis.

Sex- and Developmental Stage-Related Differences in the Hepatic Transcriptome of Japanese Quail (Coturnix japonica) Exposed to 17β-Trenbolone

Endocrine-disrupting chemicals can cause transcriptomic changes that may disrupt biological processes associated with reproductive function including metabolism, transport, and cell growth. We investigated effects from in ovo and dietary exposure to 17β-trenbolone (at 0, 1, and 10 ppm) on the Japanese quail (Coturnix japonica) hepatic transcriptome. Our objectives were to identify differentially expressed hepatic genes, assess perturbations of biological pathways, and examine sex- and developmental stage-related differences. The number of significantly differentially expressed genes was higher in embryos than in adults. Male embryos exhibited greater differential gene expression than female embryos, whereas in adults, males and females exhibited similar numbers of differentially expressed genes (>2-fold). Vitellogenin and apovitellenin-1 were up-regulated in male adults exposed to 10 ppm 17β-trenbolone, and these birds also exhibited indications of immunomodulation. Functional grouping of differentially expressed genes identified processes including metabolism and transport of biomolecules, enzyme activity, and extracellular matrix interactions. Pathway enrichment analyses identified as perturbed peroxisome proliferator-activated receptor pathway, cardiac muscle contraction, gluconeogenesis, growth factor signaling, focal adhesion, and bile acid biosynthesis. One of the primary uses of 17β-trenbolone is that of a growth promoter, and these results identify effects on mechanistic pathways related to steroidogenesis, cell proliferation, differentiation, growth, and metabolism of lipids and proteins. Environ Toxicol Chem 2021;40:2559-2570. © 2021 SETAC. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Deciphering influences of testosterone and dihydrotestosterone on lipid metabolism genes using brown trout primary hepatocytes

Despite of physiological and toxicological relevance, the potential of androgens to influence fish lipid metabolism remains poorly explored. Here, brown trout primary hepatocytes were exposed to six concentrations (1 nM to 100 μM) of dihydrotestosterone (DHT) and testosterone (T), to assess changes in the mRNA levels of genes covering diverse lipid metabolic pathways. Acsl1, essential for fatty acid activation, was up-regulated by T and DHT, whereas the lipogenic enzymes FAS and ACC were up-regulated by the highest (100 μM) concentration of T and DHT, respectively. ApoA1, the major component of high-density lipoprotein (HDL), was down-regulated by both androgens. PPARγ, linked to adipogenesis and peroxisomal β-oxidation, was down-regulated by T and DHT, while Acox1-3I, rate-limiting in peroxisomal β-oxidation, was down-regulated by T. Fabp1, StAR and LPL were not altered. Our findings suggest that androgens may impact on lipid transport, adipogenesis and fatty acid β-oxidation and promote lipogenesis in fish liver.

Exposure Impacts of Diquat dibromide herbicide formulation on amphibian larval development

Many anthropogenic chemicals in general, and specifically aquatic herbicide formulations have the potential to modulate the thyroid pathways of the endocrine system of aquatic organisms, because they are normally applied directly into the aquatic system, to manage aquatic weeds. These thyroidal effects have been widely linked with disruption in developmental and reproductive processes. In fact, the exposure impacts of many of these substances on metamorphic organisms could produce a precocious metamorphosis. Using Xenopus Metamorphosis Assay (XEMA) protocol, this study assessed the thyroidal effects of environmentally relevant concentrations of Diquat dibromide at 0.05, 0.11, and 0.14 mg/L on Xenopus laevis metamorphosis. The formulation significantly reduced both the fore and hind limb lengths, and disrupted the developmental stage at concentrations of 0.11 and 0.14 mg/L, with a median at NF-stage 57, while median of NF-stage 60 was recorded in the control. Histopathologically, although there was no significant difference in thyroid gland area, the thyroid colloidal area was significantly reduced at 0.14 mg/L, while the mean height of the thyroid follicle increased at 0.05 mg/L The result indicates an extra-thyroidal pathway, due to the dissociation between stage developmental effects and thyroid histopathology. The role of stress pathway occasioned by oxidative mode of action, involving lipid peroxidation and cell damage observed in this study need further investigation, in order to further characterize the physiological and ecological effects on wildlife.

17α-ethinylestradiol and 4-tert-octylphenol concurrently disrupt the immune response of common carp

The aquatic environment is massively polluted with endocrine-disrupting compounds (EDCs) including synthetic estrogens (e.g. 17α-ethinylestradiol, EE2) and alkylphenols (e.g. 4-tert-octylphenol, 4t-OP). A major mechanism of action for estrogenic EDCs is their interaction with estrogen receptors and consequently their modulation of the action of enzymes involved in steroid conversion e.g. aromatase CYP19. We now studied the effects of EE2 and 4t-OP on the anti-bacterial immune response of common carp. We investigated effects on the number/composition of inflammatory leukocytes and on the gene expression of mediators that regulate inflammation and EDC binding. In vitro we found that high concentrations of both EE2 and 4t-OP down-regulated IFN-γ2 and IFN-γ-dependent immune responses in LPS-stimulated monocytes/macrophages. Similarly, during bacterial infection in fish, in vivo treated with EE2 and 4t-OP, decreased gene expression of il-12p35 and of ifn-γ2 was found in the focus of inflammation. Moreover, during A. salmonicida-induced infection in EE2-treated carp, but not in fish fed with 4t-OP-treated food, we found an enhanced inflammatory reaction manifested by high number of inflammatory peritoneal leukocytes, including phagocytes and higher expression of pro-inflammatory mediators (inos, il-1β, cxcl8_l2). Furthermore, in the liver, EE2 down-regulated the expression of acute phase proteins: CRPs and C3. Importantly, both in vitro and in vivo, EDCs altered the expression of estrogen receptors: nuclear (erα and erβ) and membrane (gpr30). EDCs also induced up-regulation of the cyp19b gene. Our findings reveal that contamination of the aquatic milieu with estrogenic EDCs, may considerably violate the subtle and particular allostatic interactions between the immune response and endogenous estrogens and this may have negative consequences for fish health.

Long-term exposure to low 17α-ethinylestradiol (EE2) concentrations disrupts both the reproductive and the immune system of juvenile rainbow trout, Oncorhynchus mykiss

Estrogenic endocrine disrupting compounds (EEDCs), such as ethinylestradiol (EE2), are well studied for their impact on the reproductive system of fish. EEDCs may also impact the immune system and, as a consequence, the disease susceptibility of fish. It is currently not yet known whether the low concentrations of EEDCs that are able to disrupt the reproductive system of trout are effective in disrupting the immune system and the fish host resistance towards pathogens, too, or whether such immunodisruptive effects would occur only at higher EEDC concentrations. Therefore, in the present study we compare the effect thresholds of low 17α-ethinylestradiol concentrations (1.5 and 5.5 EE2 ng/L) on the reproductive system, the immune system, the energy expenditures and the resistance of juvenile rainbow trout (Oncorhynchus mykiss) against the parasite Tetracapsuloides bryosalmonae - the etiological agent of proliferative kidney disease (PKD) of salmonids. The parasite infection was conducted without injection and under low pathogen exposure concentrations. The disease development was followed over 130 days post infection - in the presence or absence of EE2 exposure. The results show that the long-term EE2 exposure affected, at both concentrations, reproductive parameters like the mRNA levels of hepatic vitellogenin and estrogen receptors. At the same concentrations, EE2 exposure modulated the immune parameters: mRNA levels of several immune genes were altered and the parasite intensity as well as the disease severity (histopathology) were significantly reduced in EE2-exposed fish compared to infected control fish. The combination of EE2 exposure and parasite infection was energetically costly, as indicated by the decreased values of the swim tunnel respirometry. Although further substantiation is needed, our findings suggest that EE2 exerts endocrine disruptive and immunomodulating activities at comparable effect thresholds, since reproductive and immune parameters were affected by the same, low EE2 concentrations.

Morphometric and proteomic responses of early-life stage rainbow trout (Oncorhynchus mykiss) to the aquatic herbicide diquat dibromide

The objective of this study was to examine the acute toxicity and sub-lethal effects of the commercial formulation of diquat dibromide, Reward® Landscape and Aquatic Herbicide, on multiple life stages of rainbow trout. The continuous exposure 96 h LC₅₀ derived for juvenile feeding fry aged 85 d post-hatch was 9.8 mg/L. Rainbow trout eyed embryos and juvenile feeding fry were also exposed to concentrations of Reward® ranging from 0.12 to 10 mg/L during two 24 h pulse exposures separated by 14 d of rearing in fresh water to mimic the manufacturers instructions for direct applications to water bodies. Decreased survival and body morphometrics were evident at 9.3 mg/L during the embryo/alevin exposures, but not in feeding juveniles, indicating a higher sensitivity of the early life stage fish. Quantitative proteomics and subnetwork enrichment analyses were conducted in the livers for both life stages to evaluate protein profiles after exposure to 0.37 mg/L diquat via Reward® exposure. Unique protein profiles were revealed for pre-feeding swim-up fry and for feeding juvenile fish, reflecting differences between the two life stages in sub-cellular responses after diquat dibromide exposure. Hepatic proteome effects were more dramatic in the pre-feeding swim-up fry with 315 proteins differentially expressed between the control and exposed fish while in the later life stage feeding fry, only 84 proteins were different after Reward® exposure. Exposure to Reward® significantly increased RNA/mRNA processes, induced activation of Atk/mTOR and caspase activity, and altered energy homeostasis. Proteomic alterations are associated with reduced growth observed in embryo/alevin at higher exposure concentrations, offering insight into key events underlying growth impairment within the adverse outcome pathway framework. This study is the first to report the sub-cellular and whole organism level effects of diquat dibromide in a commercial formulation and demonstrates that concentrations based on aquatic application rates alter the hepatic proteome.

The acute toxic effects of hexavalent chromium on the liver of marine medaka (Oryzias melastigma)

Chromium is toxic to marine animals and can cause damage to many of their organs, including the liver. To test the toxicity of chromium on marine organisms, we exposed the liver of the marine medaka (Oryzias melastigma) with hexavalent chromium [Cr(VI)]. Our results show that Cr enrichment in the liver demonstrates a positive correlation to the exposure concentration. With the increase of Cr(VI) concentration, pathological changes including nuclear migration, cell vacuolization, blurred intercellular gap, nuclear condensation, become noticeable. To further study changes in gene expression in the liver after Cr(VI) exposure, we used RNA-seq to compare expression profiles before and after Cr(VI) exposure. After acute Cr(VI) exposure (2.61 mg/l) for 96 h, 5862 transcripts significantly changed. It is the first time that the PPAR pathway was found to respond sensitively to Cr(VI) exposure in fish. Finally, combined with other published study, we found that there may be some difference between Cr(VI) toxicity in seawater fish and freshwater fish, due to degree of oxidative stress, distribution patterns and detailed Cr(VI) toxicological mechanisms. Not only does our study explore the mechanisms of Cr(VI) toxicity on the livers of marine medaka, it also points out different Cr(VI) toxicity levels and potential mechanisms between seawater fish and freshwater fish.

A critical review of the environmental occurrence and potential effects in aquatic vertebrates of the potent androgen receptor agonist 17β-trenbolone

Trenbolone acetate is widely used in some parts of the world for its desirable anabolic effects on livestock. Several metabolites of the acetate, including 17β-trenbolone, have been detected at low nanograms per liter concentrations in surface waters associated with animal feedlots. The 17β-trenbolone isomer can affect androgen receptor signaling pathways in various vertebrate species at comparatively low concentrations/doses. The present article provides a comprehensive review and synthesis of the existing literature concerning exposure to and biological effects of 17β-trenbolone, with an emphasis on potential risks to aquatic animals. In vitro studies indicate that, although 17β-trenbolone can activate several nuclear hormone receptors, its highest affinity is for the androgen receptor in all vertebrate taxa examined, including fish. Exposure of fish to nanograms per liter water concentrations of 17β-trenbolone can cause changes in endocrine function in the short term, and adverse apical effects in longer exposures during development and reproduction. Impacts on endocrine function typically are indicative of inappropriate androgen receptor signaling, such as changes in sex steroid metabolism, impacts on gonadal stage, and masculinization of females. Exposure of fish to 17β-trenbolone during sexual differentiation in early development can greatly skew sex ratios, whereas adult exposures can adversely impact fertility and fecundity. To fully assess ecosystem-level risks, additional research is warranted to address uncertainties as to the degree/breadth of environmental exposures and potential population-level effects of 17β-trenbolone in sensitive species. Environ Toxicol Chem 2018;37:2064-2078. Published 2018 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Transcriptome analysis reveals novel insights into the response of low-dose benzo(a)pyrene exposure in male tilapia

Despite a wide number of toxicological studies that describe benzo[a]pyrene (BaP) effects, the metabolic mechanisms that underlie these effects in fish are largely unknown. Of great concern is the presence of BaP in aquatic systems, especially those in close proximity to human activity leading to consumption of potentially contaminated foods. BaP is a known carcinogen and it has been reported to have adverse effects on the survival, development and reproduction of fish. The purpose of this study was to investigate if a low dose of BaP can alter genes and key metabolic pathways in the liver and testis in male adult tilapia, and whether these could be associated with biological endpoints disruption. We used both high-throughput RNA-Sequencing to assess whole genome gene expression following repeated intraperitoneal injections of 3 mg/kg of BaP (every 6 days for 26 days) and morphometric endpoints as indicators of general health. Condition factor (K) along with hepatosomatic and gonadosomatic indices (morphometric parameters) were significantly lower in BaP-treated fish than in controls. BaP exposure induced important changes in the gene expression pattern in liver and testis as revealed by both Pathway and Gene Ontology (GO) analyses. Alterations that were shared by both tissues included arachidonic acid metabolism, androgen receptor to prostate-specific antigen signaling, and insulin-associated effects on lipogenesis. The most salient liver-specific effects included: biological processes involved in detoxification, IL6-associated insulin resistance, mTOR hyperactivation, mitotic cytokinesis, spindle pole and microtubule binding. BaP effects that were confined to the testis included: immune system functions, inflammatory response, estrogen and androgen metabolic pathways. Taken together, gene expression and morphometric end point data indicate that the reproductive success of adult male tilapia could be compromised as a result of BaP exposure. These results constitute new insights on the mechanism of action of low dose BaP in a non-model organism (tilapia).

Gene Expression Profiling in Fish Toxicology: A Review

In this review, we present an overview of transcriptomic responses to chemical exposures in a variety of fish species. We have discussed the use of several molecular approaches such as northern blotting, differential display reverse transcription-polymerase chain reaction (DDRT-PCR), suppression subtractive hybridization (SSH), real time quantitative PCR (RT-qPCR), microarrays, and next-generation sequencing (NGS) for measuring gene expression. These techniques have been mainly used to measure the toxic effects of single compounds or simple mixtures in laboratory conditions. In addition, only few studies have been conducted to examine the biological significance of differentially expressed gene sets following chemical exposure. Therefore, future studies should focus more under field conditions using a multidisciplinary approach (genomics, proteomics and metabolomics) to understand the synergetic effects of multiple environmental stressors and to determine the functional significance of differentially expressed genes. Nevertheless, recent developments in NGS technologies and decreasing costs of sequencing holds the promise to uncover the complexity of anthropogenic impacts and biological effects in wild fish populations.

17β-Trenbolone exposure programs metabolic dysfunction in larval medaka

Here, we used physiological and transcriptomic analyses to evaluate the effects of 17β-trenbolone (TB) on metabolism during the early life stage of medaka (Oryzias latipes). In the physiological experiments, sex reversal rates increased continuously in proportion to TB concentrations (2-100 ng/L), and were 100% (all males) in the 200 ng/L treatment group. TB caused a significant increase in the gonadosomatic index of females at concentrations of 60 and 100 ng/L. These females exhibited swollen abdomens and decreased egg production and fertility. Significant increases were observed in the body mass index of these females. TB caused decreased fertility in males at concentrations >20 ng/L, but no other effects were observed. In the transcriptomic (microarray) experiments, larvae were exposed to TB for up to 7 d. Analyses using the KEGG Orthology Database revealed that predominant categories of significantly upregulated genes included "lipid metabolism" and "metabolism of terpenoids and polyketides." Thirteen genes (including those for hydroxymethylglutaryl-CoA synthase, cytoplasmic synthase, and lanosterol synthase) related to cholesterol biosynthesis via the mevalonate pathway were highlighted in these categories. Reverse transcriptase-polymerase chain reaction analyses were consistent with the microarray results, in terms of the direction and magnitude of change to gene expression. Among the downregulated genes, angiopoietin-like 4 and mitochondrial uncoupling protein 1, which are inversely correlated with obesity, were detected in the TB treatments. In conclusion, the results suggest that the exposure of females to TB during the early life stage may cause metabolic dysfunctions, including obesity and disrupted cholesterol synthesis. © 2015 Wiley Periodicals, Inc. Environ Toxicol 31: 1539-1551, 2016.

Prioritization of contaminated watercourses using an integrated biomarker approach in caged carp

Because of the ever increasing complexity of environmental contamination profiles, there are limitations to the use of analytical pollutant measurements for monitoring and prioritization of watercourses. The potential of biomarkers has been debated for many years, especially in laboratory settings, but there is a need for studies evaluating these approaches in the field. We evaluated the usefulness of a selection of biomarkers, mostly indicators of general physiological status and common stress responses such as oxidative stress, to discriminate among environmental pollution profiles, with the aim of prioritizing contaminated watercourses for targeted remediation efforts. To this end, juvenile common carp (Cyprinus carpio Lin.) were exposed in cages in the field to Flemish watercourses with varying pollution profiles. After six weeks of exposure, the bioaccumulation of key pollutants was measured, and a set of organismal, biochemical and transcriptional endpoints was determined in several tissue types. After data integration a discrete set of 14 parameters was identified, that could successfully distinguish all watercourses from each other. We show that an integrated biomarker approach, mainly targeting common stress responses, can offer the resolving power to discriminate among environmentally relevant exposure scenarios, and a means to prioritize watercourses for targeted remediation.

Contaminant-induced oxidative stress in fish: a mechanistic approach

The presence of reactive oxygen species (ROS) in living organisms was described more than 60 years ago and virtually immediately it was suggested that ROS were involved in various pathological processes and aging. The state when ROS generation exceeds elimination leading to an increased steady-state ROS level has been called "oxidative stress." Although ROS association with many pathological states in animals is well established, the question of ROS responsibility for the development of these states is still open. Fish represent the largest group of vertebrates and they inhabit a broad range of ecosystems where they are subjected to many different aquatic contaminants. In many cases, the deleterious effects of contaminants have been connected to induction of oxidative stress. Therefore, deciphering of molecular mechanisms leading to such contaminant effects and organisms' response may let prevent or minimize deleterious impacts of oxidative stress. This review describes general aspects of ROS homeostasis, in particular highlighting its basic aspects, modification of cellular constituents, operation of defense systems and ROS-based signaling with an emphasis on fish systems. A brief introduction to oxidative stress theory is accompanied by the description of a recently developed classification system for oxidative stress based on its intensity and time course. Specific information on contaminant-induced oxidative stress in fish is covered in sections devoted to such pollutants as metal ions (particularly iron, copper, chromium, mercury, arsenic, nickel, etc.), pesticides (insecticides, herbicides, and fungicides) and oil with accompanying pollutants. In the last section, certain problems and perspectives in studies of oxidative stress in fish are described.

Defining the role of omics in assessing ecosystem health: Perspectives from the Canadian environmental monitoring program

Scientific reviews and studies continue to describe omics technologies as the next generation of tools for environmental monitoring, while cautioning that there are limitations and obstacles to overcome. However, omics has not yet transitioned into national environmental monitoring programs designed to assess ecosystem health. Using the example of the Canadian Environmental Effects Monitoring (EEM) program, the authors describe the steps that would be required for omics technologies to be included in such an established program. These steps include baseline collection of omics endpoints across different species and sites to generate a range of what is biologically normal within a particular ecosystem. Natural individual variability in the omes is not adequately characterized and is often not measured in the field, but is a key component to an environmental monitoring program, to determine the critical effect size or action threshold for management. Omics endpoints must develop a level of standardization, consistency, and rigor that will allow interpretation of the relevance of changes across broader scales. To date, population-level consequences of routinely measured endpoints such as reduced gonad size or intersex in fish is not entirely clear, and the significance of genome-wide molecular, proteome, or metabolic changes on organism or population health is further removed from the levels of ecological change traditionally managed. The present review is not intended to dismiss the idea that omics will play a future role in large-scale environmental monitoring studies, but rather outlines the necessary actions for its inclusion in regulatory monitoring programs focused on assessing ecosystem health.

Retene causes multifunctional transcriptomic changes in the heart of rainbow trout (Oncorhynchus mykiss) embryos

Fish are particularly sensitive to aryl hydrocarbon receptor (AhR)-mediated developmental toxicity. The molecular mechanisms behind these adverse effects have remained largely unresolved in salmonids, and for AhR-agonistic polycyclic aromatic hydrocarbons (PAHs). This study explored the cardiac transcriptome of rainbow trout (Oncorhynchus mykiss) eleuteroembryos exposed to retene, an AhR-agonistic PAH. The embryos were exposed to retene (nominal concentration 32 μg/L) and control, their hearts were collected before, at and after the onset of the visible signs of developmental toxicity, and transcriptomic changes were studied by microarray analysis. Retene up- or down-regulated 122 genes. The largest Gene Ontology groups were signal transduction, transcription, apoptosis, cell growth, cytoskeleton, cell adhesion/mobility, cardiovascular development, xenobiotic metabolism, protein metabolism, lipid metabolism and transport, and amino acid metabolism. Together these findings suggest that retene affects multiple signaling cascades in the heart of rainbow trout embryos, and potentially disturbs processes related to cardiovascular development and function.

Challenges for using quantitative PCR test batteries as a TIE-type approach to identify metal exposure in benthic invertebrates

The epibenthic amphipod Melita plumulosa shows unique gene expression profiles when exposed to different contaminants. We hypothesized that specific changes in transcript abundance could be used in a battery of quantitative polymerase chain reaction (qPCR) assays as a toxicity identification evaluation (TIE)-like approach to identify the most relevant stressor in field-contaminated sediments. To test this hypothesis, seven candidate transcriptomic markers were selected, and their specificity following metal exposure was confirmed. The performance of these markers across different levels of added metals was verified. The ability of these transcripts to act as markers was tested by exposing amphipods to metal-contaminated field-collected sediments and measuring changes in transcript abundance via qPCR. For two of the three sediments tested, at least some of the transcriptomic patterns matched our predictions, suggesting that they would be effective in helping to identify metal exposure in field sediments. However, following exposure to the third sediment, transcriptomic patterns were unlike our predictions. These results suggest that the seven transcripts may be insufficient to discern individual contaminants from complex mixtures and that microarray or RNA-Seq global gene expression profiles may be more effective for TIE. Changes in transcriptomics based on laboratory exposures to single compounds should be carefully validated before the results are used to analyze mixtures.

Dose-dependent effects of metals on gene expression in the sydney rock oyster, Saccostrea glomerata

In the current study, we tested the effects of common environmental contaminants (the metals zinc and lead) on gene expression in Sydney rock oysters (Saccrostrea glomerata). Oysters were exposed to a range of metal concentrations under controlled laboratory conditions. The expression of 14 putative stress response genes was then measured using quantitative, real-time (q) PCR. The expression of all 14 genes was significantly affected (p < 0.05 vs. nonexposed controls) by at least one of the metals, and by at least one dose of metal. For 5 of the 14 target genes (actin, calmodulin, superoxide dismutase, topoisomerase I, and tubulin) the alteration of expression relative to controls was highest at intermediate (rather than high) doses of metals. Such responses may reflect adaptive (acclimation) reactions in gene expression at low to intermediate doses of contaminants, followed by a decline in expression resulting from exposure at higher doses. The data are discussed in terms of the intracellular pathways affected by metal contamination, and the relevance of such gene expression data to environmental biomonitoring.

Adaptation of a Bioinformatics Microarray Analysis Workflow for a Toxicogenomic Study in Rainbow Trout

Sex steroids play a key role in triggering sex differentiation in fish, the use of exogenous hormone treatment leading to partial or complete sex reversal. This phenomenon has attracted attention since the discovery that even low environmental doses of exogenous steroids can adversely affect gonad morphology (ovotestis development) and induce reproductive failure. Modern genomic-based technologies have enhanced opportunities to find out mechanisms of actions (MOA) and identify biomarkers related to the toxic action of a compound. However, high throughput data interpretation relies on statistical analysis, species genomic resources, and bioinformatics tools. The goals of this study are to improve the knowledge of feminisation in fish, by the analysis of molecular responses in the gonads of rainbow trout fry after chronic exposure to several doses (0.01, 0.1, 1 and 10 μg/L) of ethynylestradiol (EE2) and to offer target genes as potential biomarkers of ovotestis development. We successfully adapted a bioinformatics microarray analysis workflow elaborated on human data to a toxicogenomic study using rainbow trout, a fish species lacking accurate functional annotation and genomic resources. The workflow allowed to obtain lists of genes supposed to be enriched in true positive differentially expressed genes (DEGs), which were subjected to over-representation analysis methods (ORA). Several pathways and ontologies, mostly related to cell division and metabolism, sexual reproduction and steroid production, were found significantly enriched in our analyses. Moreover, two sets of potential ovotestis biomarkers were selected using several criteria. The first group displayed specific potential biomarkers belonging to pathways/ontologies highlighted in the experiment. Among them, the early ovarian differentiation gene foxl2a was overexpressed. The second group, which was highly sensitive but not specific, included the DEGs presenting the highest fold change and lowest p-value of the statistical workflow output. The methodology can be generalized to other (non-model) species and various types of microarray platforms.

The trenbolone acetate affects the immune system in rainbow trout, Oncorhynchus mykiss

In aquatic systems, the presence of endocrine-disrupting chemicals (EDC) can disrupt the reproductive function but also the immune system of wildlife. Some studies have investigated the effects of androgens on the fish immune parameters but the mechanisms by which the xenoandrogens alter the immunity are not well characterized. In order to test the effects of trenbolone acetate (TbA) on fish immune system, we exposed rainbow trout male juveniles during three weeks to TbA levels at 0.1 and 1μg/L. The present results suggest that TbA impacts, in a tissue-dependent manner, the rainbow trout immunity by affecting primarily the humoral immunity. Indeed, TbA inhibited lysozyme activity in plasma and liver and enhanced the alternative complement pathway activity (ACH50) in kidney. In plasma, the modulation of the complement system was time-dependent. The mRNA expression of genes encoding some cytokines such as renal TGF-β1, TNF-α in skin and hepatic IL-1β was also altered in fish exposed to TbA. Regarding the cellular immunity, no effect was observed on the leucocyte population. However, the expression of genes involved in the development and maturation of lymphoid cells (RAG-1 and RAG-2) was decreased in TbA-treated fish. Among those effects, we suggest that the modulation of RAG-1 and mucus apolipoprotein-A1 gene expression as well as plasma and hepatic lysozyme activities are mediated through the action of the androgen receptor. All combined, we conclude that trenbolone affects the rainbow trout immunity.

Influence of short-term exposure to low levels of 17α-ethynylestradiol on expression of genes involved in immunity and on immune parameters in rainbow trout, Oncorhynchus mykiss

Fish are exposed to endocrine-disrupting chemicals (EDC), which are well known to disturb not only the reproductive system but also the immune system in vertebrates. However, the mechanisms by which these compounds are able to modify fish immunity are not well understood. In order to test the EE2 effects on immunity in selected organs, we exposed rainbow trout male juveniles for 3 weeks to EE2 concentrations ranging from 0.01 to 1 μg/L. The results of this study suggest that EE2 affects the immunity of rainbow trout in a tissue dependent manner. This molecule affects both cellular and humoral immune systems. Indeed, blood leukocyte populations, as well as hepatic and plasma lysozyme, plasma MPO and renal complement activities, are modulated by EE2. Moreover, EE2 alters the gene expression of some mucus compounds, hepatic expression of complement sub-unit and lysozyme, or genes involved in the hepatic phagocytosis and transport of immunoglobulin across the liver.

The role of biomarkers in the assessment of aquatic ecosystem health

Ensuring the health of aquatic ecosystems and identifying species at risk from the detrimental effects of environmental contaminants can be facilitated by integrating analytical chemical analysis with carefully selected biological endpoints measured in tissues of species of concern. These biological endpoints include molecular, biochemical, and physiological markers (i.e., biomarkers) that when integrated, can clarify issues of contaminant bioavailability, bioaccumulation, and ecological effects while enabling a better understanding of the effects of nonchemical stressors. In the case of contaminant stressors, an understanding of chemical modes of toxicity can be incorporated with diagnostic markers of aquatic animal physiology to help understand the health status of aquatic organisms in the field. Furthermore, new approaches in functional genomics and bioinformatics can help discriminate individual chemicals, or groups of chemicals among complex mixtures that may contribute to adverse biological effects. Although the use of biomarkers is not a new paradigm, such approaches have been underused in the context of ecological risk assessment and natural resource damage assessment. From a regulatory standpoint, these approaches can help better assess the complex effects from coastal development activities to assessing ecosystem integrity pre- and post development or site remediation.

Correlation of gene expression and contaminant concentrations in wild largescale suckers: a field-based study

Toxic compounds such as organochlorine pesticides (OCs), polychlorinated biphenyls (PCBs), and polybrominated diphenyl ether flame retardants (PBDEs) have been detected in fish, birds, and aquatic mammals that live in the Columbia River or use food resources from within the river. We developed a custom microarray for largescale suckers (Catostomus macrocheilus) and used it to investigate the molecular effects of contaminant exposure on wild fish in the Columbia River. Using Significance Analysis of Microarrays (SAM) we identified 72 probes representing 69 unique genes with expression patterns that correlated with hepatic tissue levels of OCs, PCBs, or PBDEs. These genes were involved in many biological processes previously shown to respond to contaminant exposure, including drug and lipid metabolism, apoptosis, cellular transport, oxidative stress, and cellular chaperone function. The relation between gene expression and contaminant concentration suggests that these genes may respond to environmental contaminant exposure and are promising candidates for further field and laboratory studies to develop biomarkers for monitoring exposure of wild fish to contaminant mixtures found in the Columbia River Basin. The array developed in this study could also be a useful tool for studies involving endangered sucker species and other sucker species used in contaminant research.

Assessing mechanisms of toxicant response in the amphipod Melita plumulosa through transcriptomic profiling

This study describes the function of transcripts with altered abundance in the epibenthic amphipod, Melita plumulosa, following whole-sediment exposure to a series of common environmental contaminants. M. plumulosa were exposed for 48 h to sediments spiked and equilibrated with the following contaminants at concentrations predicted to cause sublethal effects to reproduction: porewater ammonia 30 mg L(-1); bifenthrin at 100 μg kg(-1); fipronil at 50 μg kg(-1); 0.6% diesel; 0.3% crude oil; 250 mg Cu kg(-1); 400 mg Ni kg(-1); and 400 mg Zn kg(-1). RNA was extracted and hybridized against a custom Agilent microarray developed for this species. Although the microarray represented a partial transcriptome and not all features on the array could be annotated, unique transcriptomic profiles were generated for each of the contaminant exposures. Hierarchical clustering grouped the expression profiles together by contaminant class, with copper and zinc, the petroleum products and nickel, and the pesticides each forming a distinct cluster. Many of the transcriptional changes observed were consistent with patterns previously described in other crustaceans. The changes in the transcriptome demonstrated that contaminant exposure caused changes in digestive function, growth and moulting, and the cytoskeleton following metal exposure, whereas exposure to petroleum products caused changes in carbohydrate metabolism, xenobiotic metabolism and hormone cycling. Functional analysis of these gene expression profiles can provide a better understanding of modes of toxic action and permits the prediction of mixture effects within contaminated ecosystems.

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.

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.

Using transcriptomic profiles in the diatom Phaeodactylum tricornutum to identify and prioritize stressors

The transcriptomic profile of the marine diatom, Phaeodactylum tricornutum, exposed to several ecologically relevant stressors, was used to develop toxicity identification evaluation (TIE)-like gene expression assays. Algal growth inhibition was measured by flow cytometry to determine exposure concentrations that elicited a sublethal toxic response. P. tricornutum was exposed to concentrations of copper (2 μg L⁻¹), cadmium (5 μg L⁻¹), silver (20 μg L⁻¹), simazine (75 μg L⁻¹), the water accommodated fraction (WAF) of weathered crude oil (5 mg L⁻¹), 50 μg L⁻¹ ammonia, a decreased salinity treatment (15‰), and a mixture exposure of ammonia, decreased salinity and cadmium (10 μg L⁻¹). Analysis of the gene expression via microarray indicated that unique transcriptomic signals were generated for each of the individual treatments. Transcriptomic profiles of ammonia and the mixture treatment overlapped substantially. Photosynthesis related transcripts were altered in the simazine (herbicide) treatment. A transcript involved in degrading hydrocarbons, dioxygenase, had increased abundance after crude oil exposure. Overall, transcriptomic responses in the different treatments were associated with stress responses, membrane transport, transcription and translation and could be linked to contaminant mode of action. The transcriptomic profiles were used to design real-time (quantitative) polymerase chain reaction (qPCR) assays that would link changes in transcript abundance to a particular stressor in a TIE-based approach. At least one transcript for each contaminant tested (copper, cadmium, silver, salinity and ammonia) responded exclusively to that contaminant. With further development of additional transcriptomic markers for each contaminant, this new approach has potential to enhance traditional toxicology bioassays by providing additional lines of evidence to identify biologically relevant stressors within a contaminated ecosystem based on changes in the transcriptomic profile.

Induction and recovery of estrogenic effects after short-term 17β-estradiol exposure in juvenile rainbow trout (Oncorhynchus mykiss)

Estrogenic compounds found in the aquatic environment include natural and synthetic estrogen hormones as well as other less potent estrogenic xenobiotics. In this study, a comprehensive approach was used to examine effects on fish endocrine system endpoints during a short-term xenoestrogen exposure as well as after post-exposure recovery. Rainbow trout (Oncorhynchus mykiss) were exposed to an aqueous 17β-estradiol (E2) concentration of 0.473 μg l(-1) for 2 and 7 days (d) followed by a 14-d recovery period. At d2 and d7, plasma E2 concentrations in treated fish were 458- and 205-fold higher than in control fish and 23- and 16-fold higher than the exposure water concentration. E2 treatment resulted in significant increases in hepatosomatic index (HSI), plasma vitellogenin (VTG) protein concentrations, and liver VTG and estrogen receptor alpha mRNA levels. All of these parameters, with the exception of plasma VTG protein, returned to baseline values during the recovery period. Plasma cortisol concentrations were unaffected by treatment. This research shows varied time frames of the estrogen-responsive molecular-, biochemical-, and tissue-level alterations, as well as their persistence, in juvenile rainbow trout treated with aqueous E2. These results have implications for feral rainbow trout exposed to xenoestrogens and indicate the importance of evaluating a comprehensive suite of endpoints in assessing the impact of this type of environmental contaminant.

The effects of 17-α-ethinylestradiol (EE2) on molecular signaling cascades in mummichog (Fundulus heteroclitus)

Exposures to ≤10 ng/L of 17-α-ethinylestradiol (EE2) will reduce or shut down egg production in freshwater fish models, while mummichog (Fundulus heteroclitus), an estuarine species, are able to produce eggs at EE2 concentrations >3000 ng/L. The objective of this study was to gain mechanistic insight into how mummichog are able to produce eggs during exposures to high EE2. Mummichog were exposed to 0, 50 or 250 ng/L of EE2 for 14 d. There were no changes in gonadosomatic index, liversomatic index, gonad development, or plasma estradiol levels after exposure to EE2. However, testosterone significantly decreased with EE2 exposures (50, 250 ng/L). Microarray analysis in the liver revealed that cell processes associated with lipids were affected by EE2 at the transcriptome level. Based on the transcriptomics data, we hypothesize that mummichog are able to maintain lipid transport and uptake into the ovary and this may be associated with apolipoproteins, facilitating normal oocyte development. Novel gene regulatory networks for protein modification targets were also constructed to learn more about the potential roles of estrogens in the teleost liver. Although post-translational modifications (PTMs) are important regulatory mechanisms, the roles of PTMs in protein regulation in fish and the susceptibility of PTMs to aquatic pollutants are largely unexplored and may offer novel insight into mechanisms of endocrine disruption.

Transcriptomic analyses of sexual dimorphism of the zebrafish liver and the effect of sex hormones

The liver is one of the most sex-dimorphic organs in both oviparous and viviparous animals. In order to understand the molecular basis of the difference between male and female livers, high-throughput RNA-SAGE (serial analysis of gene expression) sequencing was performed for zebrafish livers of both sexes and their transcriptomes were compared. Both sexes had abundantly expressed genes involved in translation, coagulation and lipid metabolism, consistent with the general function of the liver. For sex-biased transcripts, from in addition to the high enrichment of vitellogenin transcripts in spawning female livers, which constituted nearly 80% of total mRNA, it is apparent that the female-biased genes were mostly involved in ribosome/translation, estrogen pathway, lipid transport, etc, while the male-biased genes were enriched for oxidation reduction, carbohydrate metabolism, coagulation, protein transport and localization, etc. Sexual dimorphism on xenobiotic metabolism and anti-oxidation was also noted and it is likely that retinol x receptor (RXR) and liver x receptor (LXR) play central roles in regulating the sexual differences of lipid and cholesterol metabolisms. Consistent with high ribosomal/translational activities in the female liver, female-biased genes were significantly regulated by two important transcription factors, Myc and Mycn. In contrast, Male livers showed activation of transcription factors Ppargc1b, Hnf4a, and Stat4, which regulate lipid and glucose metabolisms and various cellular activities. The transcriptomic responses to sex hormones, 17β-estradiol (E2) or 11-keto testosterone (KT11), were also investigated in both male and female livers and we found that female livers were relatively insensitive to sex hormone disturbance, while the male livers were readily affected. E2 feminized male liver by up-regulating female-biased transcripts and down-regulating male-biased transcripts. The information obtained in this study provides comprehensive insights into the sexual dimorphism of zebrafish liver transcriptome and will facilitate further development of the zebrafish as a human liver disease model.

Impact of the redox-cycling herbicide diquat on transcript expression and antioxidant enzymatic activities of the freshwater snail Lymnaea stagnalis

The presence of pesticides in the environment results in potential unwanted effects on non-target species. Freshwater organisms inhabiting water bodies adjacent to agricultural areas, such as ditches, ponds and marshes, are good models to test such effects as various pesticides may reach these habitats through several ways, including aerial drift, run-off, and drainage. Diquat is a non-selective herbicide used for crop protection or for weed control in such water bodies. In this study, we investigated the effects of diquat on a widely spread aquatic invertebrate, the holarctic freshwater snail Lymnaea stagnalis. Due to the known redox-cycling properties of diquat, we studied transcript expression and enzymatic activities relative to oxidative and general stress in the haemolymph and gonado-digestive complex (GDC). As diquat is not persistent, snails were exposed for short times (5, 24, and 48 h) to ecologically relevant concentrations (22.2, 44.4, and 222.2 μg l(-1)) of diquat dibromide. RT-qPCR was used to quantify the transcription of genes encoding catalase (cat), a cytosolic superoxide dismutase (Cu/Zn-sod), a selenium-dependent glutathione peroxidase (gpx), a glutathione reductase (gred), the retinoid X receptor (rxr), two heat shock proteins (hsp40 and hsp70), cortactin (cor) and the two ribosomal genes r18S and r28s. Enzymatic activities of SOD, Gpx, Gred and glutathione S-transferase (GST) were investigated in the GDC using spectrophoto/fluorometric methods. Opposite trends were obtained in the haemolymph depending on the herbicide concentration. At the lowest concentration, effects were mainly observed after 24 h of exposure, with over-transcription of cor, hsp40, rxr, and sod, whereas higher concentrations down-regulated the expression of most of the studied transcripts, especially after 48 h of exposure. In the GDC, earlier responses were observed and the fold-change magnitude was generally much higher: transcription of all target genes increased significantly (or non-significantly for cat) after 5 h of exposure, and went back to control levels afterwards, suggesting the onset of an early response to oxidative stress associated to the unbalance of reactive oxygen species (ROS) in hepatocytes. Although increases obtained for Gred and SOD activities were globally consistent with their respective transcript expressions, up-regulation of transcription was not always correlated with increase of enzymatic activity, indicating that diquat might affect steps downstream of transcription. However, constitutive levels of enzymatic activities were at least maintained. In conclusion, diquat was shown to affect expression of the whole set of studied transcripts, reflecting their suitability as markers of early response to oxidative stress in L. stagnalis.

Alteration of immune function endpoints and differential expression of estrogen receptor isoforms in leukocytes from 17β-estradiol exposed rainbow trout (Oncorhynchus mykiss)

While the endocrine system is known to modulate immune function in vertebrates, the role of 17β-estradiol (E2) in cellular immune function of teleosts is poorly understood. The cellular and molecular responses of juvenile rainbow trout (Oncorhynchus mykiss) to E2 treatment were evaluated by exposing fish to 0.47±0.02μg/L E2 (mean±SEM) for either 2 or 7d, with a subsequent 14d recovery period. After 2 and 7d of exposure to E2, hematocrit was significantly lower than in control fish. Lipopolysaccharide-induced lymphocyte proliferation was elevated on day 2 and concanavalin A-induced lymphocyte proliferation was reduced following 7d of E2 exposure. Four estrogen receptor (ER) transcripts were identified in purified trout head kidney leukocytes (HKL) and peripheral blood leukocytes (PBL). While the mRNA abundance of ERβ1 and ERβ2 was unaffected by treatment, ERα1 was up-regulated in HKL and PBL following 7d of E2 exposure. ERα2 was up-regulated in HKL after 7d of E2 exposure, but down-regulated in PBL after 2 and 7d of treatment. All parameters that were altered during the E2 exposure period returned to baseline levels following the recovery period. This study reports the presence of the full repertoire of ERs in purified HKL for the first time, and demonstrates that ERα transcript abundance in leukocytes can be regulated by waterborne E2 exposure. It also demonstrated that physiologically-relevant concentrations of E2 can modulate several immune functions in salmonids, which may have widespread implications for xenoestrogen-associated immunotoxicity in feral fish populations inhabiting contaminated aquatic environments.

Systematic evaluation of medium-throughput mRNA abundance platforms

Profiling of mRNA abundances with high-throughput platforms such as microarrays and RNA-seq has become an important tool in both basic and biomedical research. However, these platforms remain prone to systematic errors and have challenges in clinical and industrial applications. As a result, it is standard practice to validate a subset of key results using alternate technologies. Similarly, clinical and industrial applications typically involve transitions from a high-throughput discovery platform to medium-throughput validation ones. These medium-throughput validation platforms have high technical reproducibility and reduced sample input needs, and low sensitivity to sample quality (e.g., for processing FFPE specimens). Unfortunately, while medium-throughput platforms have proliferated, there are no comprehensive comparisons of them. Here we fill that gap by comparing two key medium-throughput platforms--NanoString's nCounter Analysis System and ABI's OpenArray System--to gold-standard quantitative real-time RT-PCR. We quantified 38 genes and positive and negative controls in 165 samples. Signal:noise ratios, correlations, dynamic range, and detection accuracy were compared across platforms. All three measurement technologies showed good concordance, but with divergent price/time/sensitivity trade-offs. This study provides the first detailed comparison of medium-throughput RNA quantification platforms and provides a template and a standard data set for the evaluation of additional technologies.

Toxicity of atrazine and nonylphenol in juvenile rainbow trout (Oncorhynchus mykiss): effects on general health, disease susceptibility and gene expression

Atrazine (ATZ) and nonylphenol (NP) are commonly identified contaminants in aquatic habitats; however, few studies have considered the impact of these endocrine disrupters on immune function and resistance to disease. This study examined the immunotoxicological effects of ATZ and NP at multiple levels of biological organization. Juvenile rainbow trout (Oncorhynchus mykiss) were exposed to a solvent control (0.00625%, v/v anhydrous ethanol), or sub-lethal concentrations of ATZ (59 μg/L and 555 μg/L) or NP (2.3 μg/L or 18 μg/L) for 4d. At the end of exposure, fish were assessed for a number of physiological endpoints, including a host resistance challenge, and liver gene expression was assessed using a salmonid microarray (cGRASP, 32K version 1). While the low ATZ and low NP treatments had no measurable effects on the physiological endpoints measured, fish exposed to the high ATZ concentration (555 μg/L) exhibited significantly elevated plasma cortisol, a decrease in SSI, and decreased lymphocytes and increased monocytes in peripheral blood, with suppression of early immune system processes apparent at the molecular level. In contrast, fish exposed to the high NP concentration (18 μg/L) showed physiological (e.g. significantly elevated LSI) and gene expression changes (e.g. induction of vitellogenin) consistent with estrogenic effects, as well as decreased lymphocytes in the peripheral blood and more limited alterations in immune system related pathways in the liver transcriptome. Fish exposed to high ATZ or NP concentrations incurred higher mortality than control fish following a disease challenge with Listonella anguillarum, while fish exposed to the lower concentrations were unaffected. Microarray analysis of the liver transcriptome revealed a total of 211 unique, annotated differentially regulated genes (DRGs) following high ATZ exposure and 299 DRGs following high NP exposure. Functional (enrichment) analysis revealed effects on immune system function, metabolism, oxygen homeostasis, cell cycle, DNA damage, and other processes affected by ATZ or NP exposure. Overall, this study provides evidence at multiple levels of biological organization that both ATZ and NP are immunotoxic at sub-lethal concentrations and highlights the potential risk posed by these chemicals to wild fish populations.

Comparison of toxicity and transcriptomic profiles in a diatom exposed to oil, dispersants, dispersed oil

Dispersants are commonly used to mitigate the impact of oil spills, however, the ecological cost associated with their use is uncertain. The toxicity of weathered oil, dispersed weathered oil, and the hydrocarbon-based dispersant Slickgone NS(®), to the diatom Phaeodactylum tricornutum has been examined using standardized toxicity tests. The assumption that most toxicity occurs via narcosis was tested by measuring membrane damage in diatoms after exposure to one of the petroleum products. The mode of toxic action was determined using microarray-based gene expression profiling in diatoms after exposure to one of the petroleum products. The diatoms were found to be much more sensitive to dispersants than to the water accommodated fraction (WAF), and more sensitive to the chemically enhanced WAF (CEWAF) than to either the WAF itself or the dispersants. Exposure to dispersants and CEWAF caused membrane damage, while exposure to WAF did not. The gene expression profiles resulting from exposure to all three petroleum mixtures were highly similar, suggesting a similar mode of action for these compounds. The observed toxicity bore no relationship to PAH concentrations in the water column or to total petroleum hydrocarbon (TPH), suggesting that an undescribed component of the oil was causing toxicity. Taken together, these results suggest that the use of dispersants to clean up oil spills will dramatically increase the oil toxicity to diatoms, and may have implications for ecological processes such as the timing of blooms necessary for recruitment.

Exploring androgen-regulated pathways in teleost fish using transcriptomics and proteomics

In the environment, there are aquatic pollutants that disrupt androgen signaling in fish. Laboratory and field-based experiments have utilized omics technologies to characterize the molecular mechanisms underlying androgen-receptor agonism/antagonism. Transcriptomics and proteomics studies with 17β-trenbolone, a growth-promoting pharmaceutical found in water systems surrounding cattle feed lots, and androgens such as 17α-methyltestosterone and 17α-methyldihydrotestosterone, have been conducted in ovary and liver of fish that include the fathead minnow (FHM) (Pimephales promelas), common carp (Cyprinus carpio), Qurt medaka (Oryzias latipes), and zebrafish (Danio rerio). In this mini-review, we survey recent omics studies in fish and reveal that, despite the diversity of species and tissues examined, there are common cellular responses that are observed with waterborne androgenic treatments. Recurring themes in gene ontology include apoptosis, transport and oxidation of lipids, synthesis and transport of hormones, immune response, protein metabolism, and cell proliferation. However, we also discuss other mechanisms other than androgen receptor (AR) activation, such as responses to toxicant stress, estrogen receptor agonism, aromatization of androgens into estrogens, and inhibitory feedback mechanisms by high levels of androgens that may also explain molecular responses in fish. To further explore androgen-responsive protein networks, a sub-network enrichment analysis was performed on protein data collected from the livers of female FHMs exposed to 17β-trenbolone. We construct a putative AR-regulated protein/cell process network in the liver that includes B-lymphocyte differentiation, xenobiotic clearance, low-density lipoprotein oxidation, proliferation of smooth muscle cells, and permeability of blood vessels. We demonstrate that construction of protein networks can offer insight into cell processes that are potentially regulated by androgens.

BDE 49 and developmental toxicity in zebrafish

The polybrominated diphenyl ethers (PBDEs) are a group of brominated flame retardants. Human health concerns of these agents have largely centered upon their potential to elicit reproductive and developmental effects. Of the various congeners, BDE 49 (2,2',4,5'-tetrabromodiphenyl ether) has been poorly studied, despite the fact that it is often detected in the tissues of fish and wildlife species. Furthermore, we have previously shown that BDE 49 is a metabolic debromination product of BDE 99 hepatic metabolism in salmon, carp and trout, underscoring the need for a better understanding of biological effects. In the current study, we investigated the developmental toxicity of BDE 49 using the zebrafish (Danio rerio) embryo larval model. Embryo and larval zebrafish were exposed to BDE 49 at either 5 hours post fertilization (hpf) or 24 hpf and monitored for developmental and neurotoxicity. Exposure to BDE 49 at concentrations of 4iμ-32 μM caused a dose-dependent loss in survivorship at 6 days post fertilization (dpf). Morphological impairments were observed prior to the onset of mortality, the most striking of which included severe dorsal curvatures of the tail. The incidence of dorsal tail curvatures was dose and time dependent. Exposure to BDE 49 caused cardiac toxicity as evidenced by a significant reduction in zebrafish heart rates at 6 dpf but not earlier, suggesting that cardiac toxicity was non-specific and associated with physiological stress. Neurobehavioral injury from BDE 49 was evidenced by an impairment of touch-escape responses observed at 5 dpf. Our results indicate that BDE 49 is a developmental toxicant in larval zebrafish that can cause morphological abnormalities and adversely affect neurobehavior. The observed toxicities from BDE 49 were similar in scope to those previously reported for the more common tetrabrominated congener, BDE 47, and also for other lower brominated PBDEs, suggest that these compounds may share similarities in risk to aquatic species.

Transcriptional evidence for low contribution of oil droplets to acute toxicity from dispersed oil in first feeding Atlantic cod (Gadus morhua) larvae

We evaluated the potential contribution of oil droplets to the toxicity of dispersed oil to first feeding fish larvae. Atlantic cod larvae were exposed to five concentrations of either artificially weathered (200°C residue) dispersed oil (D1-D5) containing oil droplets [medium size 11-13 μm based on volume] and water-soluble fraction [WSF] or the filtered dispersions containing only the corresponding equilibrium WSFs only (W1-W5). The larvae were exposed for 4 days and harvested for transcriptional analysis at 13 days post hatching. The most significant differently expressed genes were observed in cod larvae exposed to the highest concentration of the dispersed oil (containing 10.41 ± 0.46 μg ∑PAH/L), with CYP1A showing the strongest response. Functional analysis further showed that the top scored network as analyzed with Ingenuity Pathway Analysis was "Drug Metabolism, Endocrine System Development and Function, Lipid Metabolism". Oil exposure also increased the expression of genes involved in bone resorption and decreased the expression of genes related to bone formation. In conclusion, oil exposure affects drug metabolism, endocrine regulation, cell differentiation and proliferation, apoptosis, fatty acid biosynthesis and tissue development in Atlantic cod larvae. The altered gene transcription was dominated by the WSF and the corresponding oil droplet fraction only had a moderate contribution to the observed changes.

Progress and promises in toxicogenomics in aquatic toxicology: is technical innovation driving scientific innovation?

In the last decade, new technologies have been invented to analyze large amounts of information such as gene transcripts (transcriptomics), proteins (proteomics) and small cellular molecules (metabolomics). Many studies have been performed in the last few years applying these technologies to aquatic toxicology, mainly in fish. In this article, we summarize the current state of knowledge and question whether the application of modern technology for descriptive purposes truly represents scientific advancement in aquatic toxicology. We critically discuss the advantages and disadvantages of these technologies and emphasize the importance of these critical aspects. To date, these techniques have been used mainly as a proof of principle, demonstrating effects of model compounds. The potential to use these techniques to better analyze the mode-of-action of a toxicant or the effects of a compound within organisms has rarely been met. This is partly due to a lack of baseline data and the fact that the expression of mRNA and protein profiles is rarely linked to physiology or toxicologically meaningful outcomes. It seems premature to analyze mixtures or environmental samples until more is known about the expression profiles of individual toxicants. Gene transcription, protein, or metabolic data give only a partial view of these effects. Thus, we emphasize that data obtained by these technologies must be linked to physiological changes to fully understand their significance. The use of these techniques in aquatic toxicology is still in its infancy, data cannot yet be applied to environmental risk assessment or regulation until more emphasis is placed on interpreting the data within their physiological and toxicological contexts.

Hepatic in vitro toxicity assessment of PBDE congeners BDE47, BDE153 and BDE154 in Atlantic salmon (Salmo salar L.)

The brominated flame retardant congeners BDE47, BDE153 and BDE154 are among the congeners accumulating to the highest degree in fish. In order to gain knowledge about the toxicological effects of PBDEs in fish, microarray-based transcriptomic and 2D-DIGE/MALDI-TOF/TOF proteomic approaches were used to screen for effects in primary Atlantic salmon hepatocytes exposed to these congeners alone or in combination (PBDE-MIX). A small set of stress related transcripts and proteins were differentially expressed in the PBDE exposed hepatocytes. The PBDE-MIX, and BDE153 to a lesser degree, seems to have induced metabolic disturbances by affecting several pathways related to glucose homeostasis. Further, effects on cell cycle control and proliferation signal pathways in PBDE-MIX-exposed hepatocytes clearly suggest that the PBDE exposure affected cell proliferation processes. CYP1A was 7.41- and 7.37-fold up-regulated in hepatocytes exposed to BDE47 and PBDE-MIX, respectively, and was the only biotransformation pathway affected by the PBDE exposure. The factorial design and PLS regression analyses of the effect of the PBDE-MIX indicated that BDE47 contributed the most to the observed CYP1A response, suggesting that this congener should be incorporated in the toxic equivalent (TEQ) concept in future risk assessment of dioxin-like chemicals. Additionally, a significant up-regulation of the ER-responsive genes VTG and ZP3 was observed in cells exposed to BDE47 and PBDE-MIX. Further analyses suggested that BDE47 and BDE154 have an estrogenic effect in male fish. The data also suggested an antagonistic interaction between BDE153 and BDE154. In conclusion, this study shows that PBDEs can affect several biological systems in Atlantic salmon cells, and demonstrates the need for more studies on the simultaneous exposure to chemical mixtures to identify combined effects of chemicals.

Transcriptional responses in neonate and adult Daphnia magna in relation to relative susceptibility to genotoxicants

Little information is available on the responses of lower animals to genotoxic chemicals or on their sensitivity for detecting genotoxic chemicals, especially at different life-stages, despite the established use of the water flea Daphnia magna in ecotoxicity testing. Comet assay methodology was developed and applied to daphnid cells but only limited, non-statistically significant responses to the genotoxicants sodium dichromate (0.2-1 μM), chrysoidine (0.1-2 μM), and mixtures of benzo-a-pyrene (BaP) and sodium dichromate were found (from 0.01 μM BaP & 0.1 μM sodium dichromate to 0.25 μM BaP & 0.75 μM sodium dichromate). Transcriptomic analyses using Agilent D. magna oligonucleotide microarrays were undertaken to assess the effect of a mixture of sodium dichromate and BaP (designed to produce both adducted and oxidised DNA) on gene transcription. Neonates (<24h) and adults (day 7) were exposed for 6h and 24h at two combination concentration levels (0.02 μM BaP & 0.15 μM sodium dichromate and 0.1 μM BaP & 0.75 μM sodium dichromate). The greatest differences in transcriptional profile occurred between adults and neonates. Subsets of the transcriptional profiles distinguished genotoxicant-exposed animals from controls, both for neonates and adults. Higher transcript levels of DNA repair genes were found in adults and adults also displayed significant induction of DNA repair gene transcripts in response to exposure whereas neonates did not. Transcriptional changes in response to genotoxicant exposure proved more sensitive than measurement of DNA strand breaks by the Comet assay and the extensive differences in transcription between adults and neonates emphasized the importance of life stage in toxicant testing with Daphnia.

Transcriptional response of hepatic largemouth bass (Micropterus salmoides) mRNA upon exposure to environmental contaminants

Microarrays enable gene transcript expression changes in near-whole genomes to be assessed in response to environmental stimuli. We utilized oligonucleotide microarrays and subsequent gene set enrichment analysis (GSEA) to assess patterns of gene expression changes in male largemouth bass (Micropterus salmoides) hepatic tissues after a 96 h exposure to common environmental contaminants. Fish were exposed to atrazine, cadmium chloride, PCB 126, phenanthrene and toxaphene via intraperitoneal injection with target body burdens of 3.0, 0.00067, 2.5, 50 and 100 µg g(-1), respectively. This was conducted in an effort to identify potential biomarkers of exposure. The expressions of 4, 126, 118, 137 and 58 mRNA transcripts were significantly (P ≤ 0.001, fold change ≥2×) affected by exposure to atrazine, cadmium chloride, PCB 126, phenanthrene and toxaphene exposures, respectively. GSEA revealed that none, four, five, five and three biological function gene ontology categories were significantly influenced by exposure to these chemicals, respectively. We observed that cadmium chloride elicited ethanol metabolism responses, and along with PCB 126 and phenanthrene affected transcripts associated with protein biosynthesis. PCB 126, phenanthrene and toxaphene also influenced one-carbon compound metabolism while PCB 126 and phenanthrene affected mRNA transcription and mRNA export from the nucleus and may have induced an antiestrogenic response. Atrazine was found to alter the expression of few hepatic transcripts. This work has highlighted several biological processes of interest that may be helpful in the development of gene transcript biomarkers of chemical exposure in fish.

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.

Use of gene expression, biochemical and metabolite profiles to enhance exposure and effects assessment of the model androgen 17β-trenbolone in fish

The impact of exposure by water to a model androgen, 17β-trenbolone (TRB), was assessed in fathead minnows using an integrated molecular approach. This included classical measures of endocrine exposure such as impacts on testosterone (T), 17β-estradiol (E2), and vitellogenin (VTG) concentrations in plasma, as well as determination of effects on the hepatic metabolome using proton nuclear magnetic resonance spectroscopy. In addition, the rates of production of T and E2 in ovary explants were measured, as were changes in a number of ovarian gene transcripts hypothesized to be relevant to androgen exposure. A temporally intensive 16-d test design was used to assess responses both during and after the TRB exposure (i.e., depuration/recovery). This strategy revealed time-dependent responses in females (little impact was seen in the males), in which changes in T and E2 production in the ovary, as well as levels in plasma, declined rapidly (within 1 d), followed shortly by a return to control levels. Gene expression measurements revealed dynamic control of transcript levels in the ovary and suggested potential mechanisms for compensation during the exposure phase of the test. Proton nuclear magnetic resonance spectroscopy revealed a number of hepatic metabolite changes that exhibited strong time and dose dependence. Furthermore, TRB appeared to induce the hepatic metabolome of females to become more like that of males at both high test concentrations of TRB (472 ng/L) and more environmentally relevant levels (33 ng/L).

Relationships between the transcriptome and physiological indicators of reproduction in female rainbow trout over an annual cycle

Normal transcriptomic patterns along the brain-pituitary-gonad-liver (BPGL) axis should be better characterized if endocrine-disrupting compound-induced changes in gene expression are to be understood. Female rainbow trout were studied over a complete year-long reproductive cycle. Tissue samples from pituitary, ovary, and liver were collected for microarray analysis using the 16K Genomic Research on Atlantic Salmon Project (GRASP) microarray and for quantitative polymerase chain reaction measures of estrogen receptor (ER) isoform messenger RNA (mRNA) levels. Plasma was collected to determine levels of circulating estradiol-17β (E2), follicle-stimulating hormone (FSH), and luteinizing hormone (LH). As an a priori hypothesis, changes in gene expression were correlated to either circulating levels of E2, FSH, and LH, or ER mRNAs quantified by quantitative polymerase chain reaction. In the liver, most transcriptomic patterns correlated to levels of either E2, LH, or ERs. Fewer ovarian transcripts could be correlated to levels of E2, ERα, or FSH. No significant associations were obvious in the pituitary. As a post hoc hypothesis, changes in transcript abundance were compared with microarray features with known roles in gonadal maturation. Many altered transcripts in the ovary correlated to transcript levels of estradiol 17-beta-dehydrogenase 8 or 17 B HSD12, or to glycoprotein alpha chain 1 or 2. In the pituitary, genes involved with the growth axis (e.g., growth hormone, insulin-related growth factor binding protein) correlated with the most transcripts. These results suggest that transcriptional networks along the BPGL axis may be regulated by factors other than circulating steroid hormones.