Jim Lazorchak

Gene expression profiling of the androgen receptor antagonists flutamide and vinclozolin in zebrafish (Danio rerio) gonads

The studies presented in this manuscript focus on characterization of transcriptomic responses to anti-androgens in zebrafish (Danio rerio). Research on the effects of anti-androgens in fish has been characterized by a heavy reliance on apical endpoints, and molecular mechanisms of action (MOA) of anti-androgens remain poorly elucidated. In the present study, we examined effects of a short term exposure (24-96h) to the androgen receptor antagonists flutamide (FLU) and vinclozolin (VZ) on gene expression in gonads of sexually mature zebrafish, using commercially available zebrafish oligonucleotide microarrays (4×44K platform). We found that VZ and FLU potentially impact reproductive processes via multiple pathways related to steroidogenesis, spermatogenesis, and fertilization. Observed changes in gene expression often were shared by VZ and FLU, as demonstrated by overlap in differentially-expressed genes and enrichment of several common key pathways including: (1) integrin and actin signaling, (2) nuclear receptor 5A1 signaling, (3) fibroblast growth factor receptor signaling, (4) polyamine synthesis, and (5) androgen synthesis. This information should prove useful to elucidating specific mechanisms of reproductive effects of anti-androgens in fish, as well as developing biomarkers for this important class of endocrine-active chemicals.

Quantification and associated variability of induced vitellogenin gene transcripts in fathead minnow (Pimephales promelas) by quantitative real‐time polymerase chain reaction assay

Ecological risk assessors have a growing need for sensitive and rapid indicators of environmental exposures in aquatic ecosystems resulting from natural and synthetic estrogen-like compounds. Investigators developing subcellular exposure markers in traditional sentinel organisms must be vigilant about inherent variability of analyses, especially regarding regulatory and policy statements. Here, we report a quantitative real-time polymerase chain reaction (QPCR) assay for the detection of vitellogenin transcripts environmentally triggered in fathead minnows (Pimephales promelas). We demonstrate that our QPCR assay exhibits little inter- or intra-assay variability (21.7 and 11.9%, respectively). This method appears to be robust in terms of variability stemming from extrinsic sources, indicating that it may be readily transferable to laboratories having the requisite equipment. Our primary focus in development of this method derived from the observation that transcriptional responses of the vitellogenin gene (vtg) in fathead minnows demonstrated high biological variability between identically treated individuals, even under controlled laboratory conditions (coefficient of variation, > 100%). This variability was not seen in other genes from the same RNA preparations that we examined, suggesting that it is specific to the vitellogenin response. Our data and those of others suggest that variability in vtg expression is common to a number of aquatic vertebrates, which is indicative of genetic causation. Despite a relatively high degree of variability in vtg transcription, this method is sensitive enough to detect exposures of 5.0 ng 17alpha-ethinylestradiol (EE2)/L within 24 h of exposure, and it has the ability to discriminate 10.0 and 5.0 ng EE2/L within 48 h. The vitellogenin QPCR assay is a highly sensitive, comparatively rapid, and inexpensive method for the detection and characterization of exposure to environmental estrogens and estrogen mimics.

DNA Microarray‐based ecotoxicological biomarker discovery in a small fish model species

As potential biomarkers, gene classifiers are gene expression signatures or patterns capable of distinguishing biological samples belonging to different classes or conditions. This is the second of two papers on profiling gene expression in zebrafish (Danio rerio) treated with endocrine-disrupting chemicals of different modes of action, with a focus on comparative analysis of microarray data for gene classifier discovery. Various combinations of gene feature selection/class prediction algorithms were evaluated, with the use of microarray data organized by a chemical stressor or tissue type, for their accuracy in determining the class memberships of independent test samples. Two-way clustering of gene classifiers and treatment conditions offered another alternative to assess the performance of these potential biomarkers. Both gene feature selection methods and class prediction algorithms were shown to be important in identifying successful gene classifiers. The genetic algorithm and support vector machine yielded classifiers with the best prediction accuracy, regardless of sample size, nature of class prediction, and data complexity. A chemical stressor significantly altering the expression of a greater number of genes tended to generate gene classifiers with better performance. All combinations of gene feature selection/class prediction algorithms performed similarly well with data of high signal to noise ratio. Gene classifier discovery and application on the basis of individual sampling and sample data pooling, respectively, were found to enhance class predictions. Gene expression profiles of the top gene classifiers, identified from both microarray and quantitative polymerase chain reaction assays, displayed greater similarity between fadrozole and 17beta-trenbolone than either one to 17alpha-ethinylestradiol. These gene classifiers could serve as potential biomarkers of exposure to specific classes of endocrine disruptors.

A transcriptomics-based biological framework for studying mechanisms of endocrine disruption in small fish species.

This study sought to construct a transcriptomics-based framework of signal transduction pathways, transcriptional regulatory networks, and the hypothalamic-pituitary gonadal (HPG) axis in zebrafish (Danio rerio) to facilitate formulation of specific, testable hypotheses regarding the mechanisms of endocrine disruption in fish. For the analyses involved, we used data from a total of more than 300 microarrays representing 58 conditions, which encompassed 4 tissue types from zebrafish of both genders exposed for 1 of 3 durations to 10 different test chemicals (17alpha-ethynyl estradiol, fadrozole, 17beta-trenbolone, fipronil, prochloraz, flutamide, muscimol, ketoconazole, trilostane, and vinclozolin). Differentially expressed genes were identified by one class t-tests for each condition, and those with false discovery rates of less than 40% and treatment/control ratios > or =1.3-fold were mapped to orthologous human, mouse, and rat pathways by Ingenuity Pathway Analysis to look for overrepresentation of known biological pathways. To complement the analysis of known biological pathways, the genes regulated by approximately 1800 transcription factors were inferred using the ARACNE mutual information-based algorithm. The resulting gene sets for all transcriptional factors, along with a group of compiled HPG-axis genes and approximately 130 publicly available biological pathways, were analyzed for their responses to the 58 treatment conditions by Gene Set Enrichment Analysis (GSEA) and its variant, Extended-GSEA. The biological pathways and transcription factors associated with multiple distinct treatments showed substantial interactions among the HPG-axis, TGF-beta, p53, and several of their cross-talking partners. These candidate networks/pathways have a variety of profound impacts on such cellular functions as stress response, cell cycle, and apoptosis.

A quantitative real‐time polymerase chain reaction method for the analysis of vitellogenin transcripts in model and nonmodel fish species

The measurement of vitellogenin (vtg) gene transcription has been shown to be a reliable indicator of exposure to estrogenic compounds. Unfortunately, the relatively poor molecular characterization of North American fish species has hindered its application to a larger number of ecologically important species. The current research aimed to demonstrate specific amplification of vtg gene transcripts in three model (zebrafish, rainbow trout, and medaka) and six nonmodel (emerald shiner, pearl dace, smallmouth bass, creek chub, white sucker, and golden redhorse) fish species. Quantitative polymerase chain reaction (QPCR) primers for model species were designed from publicly available vtg sequences. Successful amplification of vtg was demonstrated in fish exposed to 17alpha-ethinylestradiol (EE(2)) for all model species. Vitellogenin primers for selected nonmodel species were designed from published sequences of closely related species. Multiple primers were developed targeting different regions of the vtg gene. The successful amplification of vtg was confirmed through size and sequence analysis for all nonmodel species with the exception of the white sucker, in which amplifications failed. Furthermore, QPCR primers and conditions were quantitative over five orders of magnitude in at least one species (pearl dace) exposed to 5 ng/L of EE(2) for 24 h. The selected species are found in a wide array of ecological habitats that span the United States. Inclusion of vtg transcriptional analysis for wild, ecologically relevant fish in monitoring studies may aid in understanding the extent of estrogenic exposure in aquatic ecosystems across the United States.

DNA microarray application in ecotoxicology: experimental design, microarray scanning, and factors affecting transcriptional profiles in a small fish species.

The research presented here is part of a larger study of the molecular mode of action of endocrine-disrupting chemicals targeting the hypothalamic-pituitary-gonadal axis in zebrafish (Danio rerio). It addresses several issues critical to microarray application in aquatic ecotoxicology: experimental design, microarray scanning, gene expression intensity distribution, and the effect of experimental parameters on the zebrafish transcriptome. Expression profiles from various tissues of individual zebrafish exposed to 17alpha-ethinylestradiol (30 ng/L), fadrozole (25 micro.g/L), or 17beta-trenbolone (3.0 microg/L) for 48 or 96 h were examined with the Agilent Oligo Microarray (G2518A). As a flexible and efficient alternative to the designs commonly used in microarray studies, an unbalanced incomplete block design was found to be well suited for this work, as evidenced by high data reproducibility, low microarray-to-microarray variability, and little gene-specific dye bias. Random scanner noise had little effect on data reproducibility. A low-level, slightly variable Cyanine 3 (Cy3) contaminant was revealed by hyperspectral imaging, suggesting fluorescence contamination as a potential contributor to the large variance associated with weakly expressed genes. Expression intensities of zebrafish genes were skewed toward the lower end of their distribution range, and more weakly expressed genes tended to have larger variances. Tissue type, followed in descending order by gender, chemical treatment, and exposure duration, had the greatest effect on the overall gene expression profiles, a finding potentially critical to experimental design optimization. Overall, congruence was excellent between quantitative polymerase chain reaction results and microarray profiles of 13 genes examined across a subset of 20 pairs of ovarian samples. These findings will help to improve applications of microarrays in future ecotoxicological studies.

Evaluation of microsomal and cytosolic biomarkers in a seven-day larval trout sediment toxicity test

Rainbow trout (Oncorhynchus mykiss) sac fry (larvae) were exposed to River Po sediments for 7 days. The sediments were collected in the River Po at two sites located upstream and downstream of the confluence of a polluted tributary, the River Lambro. An additional sediment treatment was also tested, in which trout larvae were kept from direct contact with the downstream sediment by interposing a Teflon net. Benzo[a]pyrene hydroxylase (AHH), ethoxyresorufin-O-deethylase (EROD), aminopyrine-N-demethylase (APDM) and UDP glucuronyl transferase (UDPGT) activities were found to be significantly induced in whole-body microsomal preparations of sac fry exposed to the downstream sediment. No significant modification was evident in any of the tested cytosolic biomarkers, i.e. glutathione reductase (GR) and glutathione peroxidase (GPx), glucose 6-phosphate dehydrogenase (G6PD), 6-phosphogluconate dehydrogenase (6PGD) and the content of nonprotein thiols (SH). With the exception of a slight induction of AHH enzyme activity, no difference could be found between fry exposed to control sediment and those screened from the downstream sediment, suggesting that direct contact with sediment was the major route of exposure to contaminants. This study demonstrates that several enzyme activities, which are known to occur in juvenile and adult rainbow trout, are also detectable at the sac-fry (larval) stage, and some of these activities can be induced by a short-term exposure to a contaminated sediment.

Transcriptional regulatory dynamics of the hypothalamic–pituitary–gonadal axis and its peripheral pathways as impacted by the 3-beta HSD inhibitor trilostane in zebrafish (Danio rerio)

To study mechanisms underlying generalized effects of 3β hydroxysteroid dehydrogenase (HSD3B) inhibition, reproductively mature zebrafish (Danio rerio) were exposed to trilostane at two dosages for 24, 48, or 96 h and their gonadal RNA samples profiled with Agilent zebrafish microarrays. Trilostane had substantial impact on the transcriptional dynamics of zebrafish, as reflected by a number of differentially expressed genes (DEGs) including transcription factors (TFs), altered TF networks, signaling pathways, and Gene Ontology (GO) biological processes. Changes in gene expression between a treatment and its control were mostly moderate, ranging from 1.3 to 2.0 fold. Expression of genes coding for HSD3B and many of its transcriptional regulators remained unchanged, suggesting transcriptional up-regulation is not a primary compensatory mechanism for HSD3B enzyme inhibition. While some trilostane-responsive TFs appear to share cellular functions linked to endocrine disruption, there are also many other DEGs not directly linked to steroidogenesis. Of the 65 significant TF networks, little similarity, and therefore little cross-talk, existed between them and the hypothalamic-pituitary-gonadal (HPG) axis. The most enriched GO biological processes are regulations of transcription, phosphorylation, and protein kinase activity. Most of the impacted TFs and TF networks are involved in cellular proliferation, differentiation, migration, and apoptosis. While these functions are fairly broad, their underlying TF networks may be useful to development of generalized toxicological screening methods. These findings suggest that trilostane-induced effects on fish endocrine functions are not confined to the HPG-axis alone. Its impact on corticosteroid synthesis could also have contributed to some system wide transcriptional changes in zebrafish observed in this study.

Discovery and validation of gene classifiers for endocrine-disrupting chemicals in zebrafish (danio rerio).

BackgroundDevelopment and application of transcriptomics-based gene classifiers for ecotoxicological applications lag far behind those of biomedical sciences. Many such classifiers discovered thus far lack vigorous statistical and experimental validations. A combination of genetic algorithm/support vector machines and genetic algorithm/K nearest neighbors was used in this study to search for classifiers of endocrine-disrupting chemicals (EDCs) in zebrafish. Searches were conducted on both tissue-specific and tissue-combined datasets, either across the entire transcriptome or within individual transcription factor (TF) networks previously linked to EDC effects. Candidate classifiers were evaluated by gene set enrichment analysis (GSEA) on both the original training data and a dedicated validation dataset.ResultsMulti-tissue dataset yielded no classifiers. Among the 19 chemical-tissue conditions evaluated, the transcriptome-wide searches yielded classifiers for six of them, each having approximately 20 to 30 gene features unique to a condition. Searches within individual TF networks produced classifiers for 15 chemical-tissue conditions, each containing 100 or fewer top-ranked gene features pooled from those of multiple TF networks and also unique to each condition. For the training dataset, 10 out of 11 classifiers successfully identified the gene expression profiles (GEPs) of their targeted chemical-tissue conditions by GSEA. For the validation dataset, classifiers for prochloraz-ovary and flutamide-ovary also correctly identified the GEPs of corresponding conditions while no classifier could predict the GEP from prochloraz-brain.ConclusionsThe discrepancies in the performance of these classifiers were attributed in part to varying data complexity among the conditions, as measured to some degree by Fisher’s discriminant ratio statistic. This variation in data complexity could likely be compensated by adjusting sample size for individual chemical-tissue conditions, thus suggesting a need for a preliminary survey of transcriptomic responses before launching a full scale classifier discovery effort. Classifier discovery based on individual TF networks could yield more mechanistically-oriented biomarkers. GSEA proved to be a flexible and effective tool for application of gene classifiers but a similar and more refined algorithm, connectivity mapping, should also be explored. The distribution characteristics of classifiers across tissues, chemicals, and TF networks suggested a differential biological impact among the EDCs on zebrafish transcriptome involving some basic cellular functions.

An integrated assessment of sediment remediation in a midwestern U.S. stream using sediment chemistry, water quality, bioassessment, and fish biomarkers

A comprehensive biological, sediment, and water quality study of the lower Little Scioto River near Marion, Ohio, USA, was undertaken to evaluate the changes or improvements in biotic measurements following the removal of creosote-contaminated sediment. The study area covered 7.5 river miles (RMs), including a remediated section between RMs 6.0 and 6.8. Fish and macroinvertebrate assemblages, fish biomarkers (i.e., polycyclic aromatic hydrocarbon [PAH] metabolite levels in white sucker [Castostomus commersoni] and common carp [Cyprinus carpio] bile and DNA damage), sediment chemistry, and water quality were assessed at five locations relative to the primary source of historical PAH contamination-upstream (RM 9.2), adjacent (RM 6.5), and downstream (RMs 5.7, 4.4, and 2.7). Overall, the biomarker results were consistent with the sediment PAH results, showing a pattern of low levels of PAH bile metabolites and DNA damage at the upstream (reference or background location), as well as the remediated section, high levels at the two immediate downstream sites, and somewhat lower levels at the furthest downstream site. Results show that remediation was effective in reducing sediment contaminant concentrations and exposure of fish to PAHs and in improving fish assemblages (60% increase in index of biotic integrity scores) in remediated river sections. Additional remedial investigation and potentially further remediation is needed to improve the downstream benthic fish community, which is still heavily exposed to PAH contaminants.