Comparing Transcriptomic Points of Departure to Apical Effect Concentrations For Larval Fathead Minnow Exposed to Chemicals with Four Different Modes Of Action.

It is postulated that below a transcriptomic-based point of departure, adverse effects are unlikely to occur, thereby providing a chemical concentration to use in screening level hazard assessment. The present study extends previous work describing a high-throughput fathead minnow assay that can provide full transcriptomic data after exposure to a test chemical. One-day post-hatch fathead minnows were exposed to ten concentrations of three representatives of four chemical modes of action: organophosphates, ecdysone receptor agonists, plant photosystem II inhibitors, and estrogen receptor agonists for 24 h. Concentration response modeling was performed on whole body gene expression data from each exposure, using measured chemical concentrations when available. Transcriptomic points of departure in larval fathead minnow were lower than apical effect concentrations across fish species but not always lower than toxic effect concentrations in other aquatic taxa like crustaceans and insects. The point of departure was highly dependent on measured chemical concentration which were often lower than the nominal concentration. Differentially expressed genes between chemicals within modes of action were compared and often showed statistically significant overlap. In addition, reproducibility between identical exposures using a positive control chemical (CuSO4) and variability associated with the transcriptomic point of departure using in silico sampling were considered. Results extend a transcriptomic-compatible fathead minnow high-throughput assay for possible use in ecological hazard screening.

Characterization of vitellogenin concentration in male fathead minnow mucus compared to plasma, and liver mRNA.

The objective of this study was to characterize vitellogenin (VTG) protein in male fathead minnow (Pimephales promelas) mucus compared with more conventional measures in plasma and mRNA isolated from liver. To assess the intensity and duration of changes in mucus VTG concentrations, male fathead minnows were exposed to 17α-ethinylestradiol (EE2) for 7 days with a subsequent depuration period of 14 days. The experiment was conducted in a flow-through system to maintain a consistent concentration of EE2 at a nominal EC50 concentration of 2.5 ng/L and high concentration of 10 ng/L as a positive control. Mucus, plasma and liver were sampled at regular intervals throughout the study. Relative abundance of vtg mRNA increased after 2 days of exposure and returned to control levels after 4 days of depuration. VTG protein concentration displayed similar induction kinetics in both mucus and plasma, however, it was found to be significantly increased after 2 days of exposure using the mucus-based assays and 7 days with the plasma-based assay. Significantly elevated levels of VTG were detected by both assays throughout the 14-day depuration period. The elimination of the laborious plasma collection step in the mucus-based workflow allowed sampling of smaller organisms where blood volume is limiting. It also resulted in significant gains in workflow efficiency, decreasing sampling time without loss of performance.

The Toxicogenome of Hyalella azteca: A Model for Sediment Ecotoxicology and Evolutionary Toxicology.

Hyalella azteca is a cryptic species complex of epibenthic amphipods of interest to ecotoxicology and evolutionary biology. It is the primary crustacean used in North America for sediment toxicity testing and an emerging model for molecular ecotoxicology. To provide molecular resources for sediment quality assessments and evolutionary studies, we sequenced, assembled, and annotated the genome of the H. azteca U.S. Lab Strain. The genome quality and completeness is comparable with other ecotoxicological model species. Through targeted investigation and use of gene expression data sets of H. azteca exposed to pesticides, metals, and other emerging contaminants, we annotated and characterized the major gene families involved in sequestration, detoxification, oxidative stress, and toxicant response. Our results revealed gene loss related to light sensing, but a large expansion in chemoreceptors, likely underlying sensory shifts necessary in their low light habitats. Gene family expansions were also noted for cytochrome P450 genes, cuticle proteins, ion transporters, and include recent gene duplications in the metal sequestration protein, metallothionein. Mapping of differentially expressed transcripts to the genome significantly increased the ability to functionally annotate toxicant responsive genes. The H. azteca genome will greatly facilitate development of genomic tools for environmental assessments and promote an understanding of how evolution shapes toxicological pathways with implications for environmental and human health.

Fish connectivity mapping: linking chemical stressors by their mechanisms of action-driven transcriptomic profiles.

BACKGROUND: A very large and rapidly growing collection of transcriptomic profiles in public repositories is potentially of great value to developing data-driven bioinformatics applications for toxicology/ecotoxicology. Modeled on human connectivity mapping (Cmap) in biomedical research, this study was undertaken to investigate the utility of an analogous Cmap approach in ecotoxicology. Over 3500 zebrafish (Danio rerio) and fathead minnow (Pimephales promelas) transcriptomic profiles, each associated with one of several dozen chemical treatment conditions, were compiled into three distinct collections of rank-ordered gene lists (ROGLs) by species and microarray platforms. Individual query signatures, each consisting of multiple gene probes differentially expressed in a chemical condition, were used to interrogate the reference ROGLs. RESULTS: Informative connections were established at high success rates within species when, as defined by their mechanisms of action (MOAs), both query signatures and ROGLs were associated with the same or similar chemicals. Thus, a simple query signature functioned effectively as an exposure biomarker without need for a time-consuming process of development and validation. More importantly, a large reference database of ROGLs also enabled a query signature to cross-interrogate other chemical conditions with overlapping MOAs, leading to novel groupings and subgroupings of seemingly unrelated chemicals at a finer resolution. This approach confirmed the identities of several estrogenic chemicals, as well as a polycyclic aromatic hydrocarbon and a neuro-toxin, in the largely uncharacterized water samples near several waste water treatment plants, and thus demonstrates its future potential utility in real world applications. CONCLUSIONS: The power of Cmap should grow as chemical coverages of ROGLs increase, making it a framework easily scalable in the future. The feasibility of toxicity extrapolation across fish species using Cmap needs more study, however, as more gene expression profiles linked to chemical conditions common to multiple fish species are needed.

Effects of Age and Exposure Duration on the Sensitivity of Early Life Stage Fathead Minnow (Pimephales promelas) to Waterborne Propranolol Exposure.

Propranolol is a heavily prescribed, nonspecific beta-adrenoceptor (bAR) antagonist frequently found in wastewater effluents, prompting concern over its potential to adversely affect exposed organisms. In the present study, the transcriptional responses of 4, 5, and 6 days postfertilization (dpf) ±1 h fathead minnow, exposed for 6, 24, or 48 h to 0.66 or 3.3 mg/L (nominal) propranolol were characterized using RNA sequencing. The number of differentially expressed genes (DEGs) was used as an estimate of sensitivity. A trend toward increased sensitivity with age was observed; fish >7 dpf at the end of exposure were particularly sensitive to propranolol. The DEGs largely overlapped among treatment groups, suggesting a highly consistent response that was independent of age. Cluster analysis was performed using normalized count data for unexposed and propranolol-exposed fish. Control fish clustered tightly by age, with fish ≥7 dpf clustering away from younger fish, reflecting developmental differences. When clustering was conducted using exposed fish, in cases where propranolol induced a minimal or no transcriptional response, the results mirrored those of the control fish and did not appreciably cluster by treatment. In treatment groups that displayed a more robust transcriptional response, the effects of propranolol were evident; however, fish <7 dpf clustered away from older fish, despite having similar numbers of DEGs. Increased sensitivity at 7 dpf coincided with developmental milestones with the potential to alter propranolol pharmacokinetics or pharmacodynamics, such as the onset of exogenous feeding and gill functionality as well as increased systemic expression of bAR. These results may have broader implications because toxicity testing often utilizes fish <4 dpf, prior to the onset of these potentially important developmental milestones, which may result in an underestimation of risk for some chemicals. Environ Toxicol Chem 2024;43:807-820. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

Linkage of genomic biomarkers to whole organism end points in a Toxicity Identification Evaluation (TIE).

Aquatic organisms are exposed to many toxic chemicals and interpreting the cause and effect relationships between occurrence and impairment is difficult. Toxicity Identification Evaluation (TIE) provides a systematic approach for identifying responsible toxicants. TIE relies on relatively uninformative and potentially insensitive toxicological end points. Gene expression analysis may provide needed sensitivity and specificity aiding in the identification of primary toxicants. The current work aims to determine the added benefit of integrating gene expression end points into the TIE process. A cDNA library and a custom microarray were constructed for the marine amphipod Ampelisca abdita. Phase 1 TIEs were conducted using 10% and 40% dilutions of acutely toxic sediment. Gene expression was monitored in survivors and controls. An expression-based classifier was developed and evaluated against control organisms, organisms exposed to low or medium toxicity diluted sediment, and chemically selective manipulations of highly toxic sediment. The expression-based classifier correctly identified organisms exposed to toxic sediment even when little mortality was observed, suggesting enhanced sensitivity of the TIE process. The ability of the expression-based end point to correctly identify toxic sediment was lost concomitantly with acute toxicity when organic contaminants were removed. Taken together, this suggests that gene expression enhances the performance of the TIE process.

Pilot testing and optimization of a larval fathead minnow high throughput transcriptomics assay.

Concentrations at which global gene expression profiles in cells or animals exposed to a test substance start to differ significantly from those of controls have been proposed as an alternative point of departure for use in screening level hazard assessment. The present study describes pilot testing of a high throughput compatible transcriptomics assay with larval fathead minnows. One day post hatch fathead minnows were exposed to eleven different concentrations of three metals, three selective serotonin reuptake inhibitors, and four neonicotinoid-like compounds for 24 h and concentration response modeling was applied to whole body gene expression data. Transcriptomics-based points of departure (tPODs) were consistently lower than effect concentrations reported in apical endpoint studies in fish. However, larval fathead minnow-based tPODs were not always lower than concentrations reported to elicit apical toxicity in other aquatic organisms like crustaceans or insects. Random in silico subsampling of data from the pilot assays was used to evaluate various assay design and acceptance considerations such as transcriptome coverage, number of replicate individuals to sequence per treatment, and minimum number of differentially expressed genes to produce a reliable tPOD estimate. Results showed a strong association between the total number of genes for which a concentration response relationship could be derived and the overall variability in the resulting tPOD estimates. We conclude that, for our current assay design and analysis pipeline, tPODs based on fewer than 15 differentially expressed genes are likely to be unreliable for screening and that interindividual variability in gene expression profiles appears to be a more significant driver of tPOD variability than sample size alone. Results represent initial steps toward developing high throughput transcriptomics assays for use in ecological hazard screening.

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

BACKGROUND: Development 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. RESULTS: Multi-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. CONCLUSIONS: The 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.

De Novo Assembly of the Nearly Complete Fathead Minnow Reference Genome Reveals a Repetitive but Compact Genome.

The fathead minnow is a widely used model organism in environmental toxicology. The lack of a high-quality fathead minnow reference genome, however, has severely hampered its uses in toxicogenomics. We present the de novo assembly and annotation of the fathead minnow genome using long PacBio reads, Bionano and Hi-C scaffolding data, and large RNA-sequencing data sets from different tissues and life stages. The new annotated fathead minnow reference genome has a scaffold N50 of 12.0 Mbp and a complete benchmarking universal single-copy orthologs score of 95.1%. The completeness of annotation for the new reference genome is comparable to that of the zebrafish GRCz11 reference genome. The fathead minnow genome, revealed to be highly repetitive and sharing extensive syntenic regions with the zebrafish genome, has a much more compact gene structure than the zebrafish genome. Particularly, comparative genomic analysis with zebrafish, mouse, and human showed that fathead minnow homologous genes are relatively conserved in exon regions but had strikingly shorter intron regions. The new fathead minnow reference genome and annotation data, publicly available from the National Center for Biotechnology Information and the University of California Santa Cruz genome browser, provides an essential resource for aquatic toxicogenomic studies in ecotoxicology and public health. Environ Toxicol Chem 2022;41:448-461. Published 2021. This article is a U.S. Government work and is in the public domain in the USA.

Effects of Metformin and its Metabolite Guanylurea on Fathead Minnow (Pimephales promelas) Reproduction.

Metformin, along with its biotransformation product guanylurea, is commonly observed in municipal wastewaters and subsequent surface waters. Previous studies in fish have identified metformin as a potential endocrine-active compound, but there are inconsistencies with regard to its effects. To further investigate the potential reproductive toxicity of metformin and guanylurea to fish, a series of experiments was performed with adult fathead minnows (Pimephales promelas). First, explants of fathead minnow ovary tissue were exposed to 0.001-100 µM metformin or guanylurea to investigate whether the compounds could directly perturb steroidogenesis. Second, spawning pairs of fathead minnows were exposed to metformin (0.41, 4.1, and 41 µg/L) or guanylurea (1.0, 10, and 100 µg/L) for 23 days to assess impacts on reproduction. Lastly, male fathead minnows were exposed to 41 µg/L metformin, 100 µg/L guanylurea, or a mixture of both compounds, with samples collected over a 96-h time course to investigate potential impacts to the hepatic transcriptome or metabolome. Neither metformin nor guanylurea affected steroid production by ovary tissue exposed ex vivo. In the 23 days of exposure, neither compound significantly impacted transcription of endocrine-related genes in male liver or gonad, circulating steroid concentrations in either sex, or fecundity of spawning pairs. In the 96-h time course, 100 µg guanylurea/L elicited more differentially expressed genes than 41 µg metformin/L and showed the greatest impacts at 96 h. Hepatic transcriptome and metabolome changes were chemical- and time-dependent, with the largest impact on the metabolome observed at 23 days of exposure to 100 µg guanylurea/L. Overall, metformin and guanylurea did not elicit effects consistent with reproductive toxicity in adult fathead minnows at environmentally relevant concentrations. Environ Toxicol Chem 2022;41:2708-2720. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Development of omics biomarkers for estrogen exposure using mRNA, miRNA and piRNAs.

The number of chemicals requiring risk evaluation exceeds our capacity to provide the underlying data using traditional methodology. This has led to an increased focus on the development of novel approach methodologies. This work aimed to expand the panel of gene expression-based biomarkers to include responses to estrogens, to identify training strategies that maximize the range of applicable concentrations, and to evaluate the potential for two classes of small non-coding RNAs (sncRNAs), microRNA (miRNA) and piwi-interacting RNA (piRNA), as biomarkers. To this end larval Pimephales promelas (96 hpf +/- 1h) were exposed to 5 concentrations of 17α- ethinylestradiol (0.12, 1.25, 2.5, 5.0, 10.0 ng/L) for 48 h. For mRNA-based biomarker development, RNA-seq was conducted across all concentrations. For sncRNA biomarkers, small RNA libraries were prepared only for the control and 10.0 ng/L EE2 treatment. In order to develop an mRNA classifier that remained accurate over the range of exposure concentrations, three different training strategies were employed that focused on 10 ng/L, 2.5 ng/L or a combination of both. Classifiers were tested against an independent test set of individuals exposed to the same concentrations used in training and subsequently against concentrations not included in model training. Both random forest (RF) and logistic regression with elastic net regularizations (glmnet) models trained on 10 ng/L EE2 performed poorly when applied to lower concentrations. RF models trained with either the 2.5 ng/L or combination (2.5 + 10 ng/L) treatments were able to accurately discriminate exposed vs. non-exposed across all but the lowest concentrations. glmnet models were unable to accurately classify below 5 ng/L. With the exception of the 10 ng/L treatment, few mRNA differentially expressed genes (DEG) were observed, however, there was marked overlap of DEGs across treatments. Overlapping DEGs have well established linkages to estrogen and several of the 81 DEGs identified in the 10 ng/L treatment have been previously utilized as estrogenic biomarkers (vitellogenin, estrogen receptor-ß). Following multiple test correction, no sncRNAs were found to be differentially expressed, however, both miRNA and piRNA classifiers were able to accurately discriminate control and 10 ng/L exposed organisms with AUCs of 0.83 and 1.0 respectively. We have developed a highly discriminative estrogenic mRNA biomarker that is accurate over a range of concentrations likely to be found in real-world exposures. We have demonstrated that both miRNA and piRNA are responsive to estrogenic exposure, suggesting the need to further investigate their regulatory roles in the estrogenic response.

Global transcriptomic profiling of microcystin-LR or -RR treated hepatocytes (HepaRG).

The canonical mode of action (MOA) of microcystins (MC) is the inhibition of protein phosphatases, but complete characterization of toxicity pathways is lacking. The existence of over 200 MC congeners complicates risk estimates worldwide. This work employed RNA-seq to provide an unbiased and comprehensive characterization of cellular targets and impacted cellular processes of hepatocytes exposed to either MC-LR or MC-RR congeners. The human hepatocyte cell line, HepaRG, was treated with three concentrations of MC-LR or -RR for 2 h. Significant reduction in cell survival was observed in LR1000 and LR100 treatments whereas no acute toxicity was observed in any MR-RR treatment. RNA-seq was performed on all treatments of MC-LR and -RR. Differentially expressed genes and pathways associated with oxidative and endoplasmic reticulum (ER) stress, and the unfolded protein response (UPR) were highly enriched by both congeners as were inflammatory pathways. Genes associated with both apoptotic and inflammatory pathways were enriched in LR1000. We present a model of MC toxicity that immediately causes oxidative stress and leads to ER stress and the activation of the UPR. Differential activation of the three arms of the UPR and the kinetics of JNK activation ultimately determine whether cell survival or apoptosis is favored. Extracellular exosomes were enrichment of by both congeners, suggesting a previously unidentified mechanism for MC-dependent extracellular signaling. The complement system was enriched only in MC-RR treatments, suggesting congener-specific differences in cellular effects. This study provided an unbiased snapshot of the early systemic hepatocyte response to MC-LR and MC-RR congeners and may explain differences in toxicity among MC congeners.

Altered gene expression in the brain and ovaries of zebrafish (Danio rerio) exposed to the aromatase inhibitor fadrozole: microarray analysis and hypothesis generation.

As part of a research effort examining system-wide responses of the hypothalamic-pituitary-gonadal (HPG) axis in fish to endocrine-active chemicals (EACs) with different modes of action, zebrafish (Danio rerio) were exposed to 25 or 100 microg/L of the aromatase inhibitor fadrozole for 24, 48, or 96 h. Global transcriptional response in brain and ovarian tissue of fish exposed to 25 microg/L of fadrozole was compared to that in control fish using a commercially available, 22,000-gene oligonucleotide microarray. Transcripts altered in brain were functionally linked to differentiation, development, DNA replication, and cell cycle. Additionally, multiple genes associated with the one-carbon pool by folate pathway (KEGG 00670) were significantly up-regulated. Transcripts altered in ovary were functionally linked to cell-cell adhesion, extracellular matrix, vasculogenesis, and development. Promoter motif analysis identified GATA-binding factor 2, Ikaros 2, alcohol dehydrogenase gene regulator 1, myoblast-determining factor, and several heat shock factors as being associated with coexpressed gene clusters that were differentially expressed following exposure to fadrozole. Based on the transcriptional changes observed, it was hypothesized that fadrozole elicits neurodegenerative stress in brain tissue and that fish cope with this stress through proliferation of radial glial cells. Additionally, it was hypothesized that changes of gene expression in the ovary of fadrozole-exposed zebrafish reflect disruption of oocyte maturation and ovulation because of impaired vitellogenesis. These hypotheses and others derived from the microarray results provide a foundation for future studies aimed at understanding responses of the HPG axis to EACs and other chemical stressors.

Multigene Biomarkers of Pyrethroid Exposure: Exploratory Experiments.

We describe initial development of microarray-based assays for detecting 4 pyrethroid pesticides (bifenthrin, cypermethrin, esfenvalerate, and permethrin) in water. To facilitate comparison of transcriptional responses with gross apical responses, we estimated concentration-mortality curves for these pyrethroids using flow-through exposures of newly hatched Daphnia magna, Pimephales promelas adults, and 24 h posthatch P. promelas. Median lethal concentration (LC50) estimates were below most reported values, perhaps attributable to the use of flow-through exposures or of measured rather than nominal concentrations. Microarray analysis of whole P. promelas larvae and brains from exposed P. promelas adults showed that assays using either tissue type can detect these pyrethroids at concentrations below LC50 values reported for between 72 and 96% of aquatic species, depending on the pesticide. These estimates are conservative because they correspond to the lowest concentrations tested. This suggests that the simpler and less expensive whole-larval assay provides adequate sensitivity for screening contexts where acute aquatic lethality is observed, but the responsible agent is not known. Gene set analysis (GSA) highlighted several Gene Ontology (GO) terms consistent with known pyrethroid action, but the implications of other GO terms are less clear. Exploration of the sensitivity of results to changes in data processing suggests robustness of the detection assay results, but GSA results were sensitive to methodological variations. Environ Toxicol Chem 2019;38:2436-2446. Published 2019 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.

Characterization of the Fundulus heteroclitus embryo transcriptional response and development of a gene expression-based fingerprint of exposure for the alternative flame retardant, TBPH (bis (2-ethylhexyl)-tetrabromophthalate).

Although alternative Flame Retardant (FR) chemicals are expected to be safer than the legacy FRs they replace, their risks to human health and the environment are often poorly characterized. This study used a small volume, fish embryo system to reveal potential mechanisms of action and diagnostic exposure patterns for TBPH (bis (2-ethylhexyl)-tetrabromophthalate), a component of several widely-used commercial products. Two different concentration of TBPH were applied to sensitive early life stages of an ecologically important test species, Fundulus heteroclitus (Atlantic killifish), with a well-annotated genome. Exposed fish embryos were sampled for transcriptomics or chemical analysis of parent compound and primary metabolite or observed for development and survival through larval stage. Global transcript profiling using RNA-seq was conducted (n = 16 per treatment) to provide a non-targeted and statistically robust approach to characterize TBPH gene expression patterns. Transcriptomic analysis revealed a dose-response in the expression of genes associated with a surprisingly limited number of biological pathways, but included the aryl hydrocarbon receptor signal transduction pathway, which is known to respond to several toxicologically-important chemical classes. A transcriptional fingerprint using Random Forests was developed that was able to perfectly discriminate exposed vs. non-exposed individuals in test sets. These results suggest that TBPH has a relatively low potential for developmental toxicity (at least in fishes), despite concerns related to its structural similarities to endocrine disrupting chemicals and that the early life stage Fundulus system may provide a convenient test system for exposure characterization. More broadly, this study advances the usefulness of a biological testing and analysis system utilizing non-targeted transcriptomics profiling and early developmental endpoints that complements current screening methods to characterize chemicals of ecological and human health concern.

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.

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 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.