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

The increased availability and use of DNA microarrays has allowed the characterization of gene expression patterns associated with exposure to different toxicants. An important question is whether toxicant induced changes in gene expression in fish are sufficiently diverse to allow for identification of specific modes of action and/or specific contaminants. In theory, each class of toxicant may generate a gene expression profile unique to its mode of toxic action. In this study, isogenic (cloned) rainbow trout Oncorhynchus mykiss were exposed to sublethal levels of a series of model toxicants with varying modes of action, including ethynylestradiol (xeno-estrogen), 2,2,4,4'-tetrabromodiphenyl ether (BDE-47, thyroid active), diquat (oxidant stressor), chromium VI, and benzo[a]pyrene (BaP) for a period of 1-3 weeks. An additional experiment measured trenbolone (anabolic steroid; model androgen) induced gene expression changes in sexually mature female trout. Following exposure, fish were euthanized, livers removed and RNA extracted. Fluorescently labeled cDNA were generated and hybridized against a commercially available Atlantic Salmon/Trout array (GRASP project, University of Victoria) spotted with 16,000 cDNA's. The slides were scanned to measure abundance of a given transcript in each sample relative to controls. Data were analyzed via Genespring (Silicon Genetics) to identify a list of up- and downregulated genes, as well as to determine gene clustering patterns that can be used as "expression signatures". The results indicate each toxicant exposure caused between 64 and 222 genes to be significantly altered in expression. Most genes exhibiting altered expression responded to only one of the toxicants and relatively few were co-expressed in multiple treatments. For example, BaP and Diquat, both of which exert toxicity via oxidative stress, upregulated 28 of the same genes, of over 100 genes altered by either treatment. Other genes associated with steroidogenesis, p450 and estrogen responsive genes appear to be useful for selectively identifying toxicant mode of action in fish, suggesting a link between gene expression profile and mode of toxicity. Our array results showed good agreement with quantitative real time polymerase chain reaction (qRT PCR), which demonstrates that the arrays are an accurate measure of gene expression. The specificity of the gene expression profile in response to a model toxicant, the link between genes with altered expression and mode of toxic action, and the consistency between array and qRT PCR results all suggest that cDNA microarrays have the potential to screen environmental contaminants for biomarkers and mode of toxic action.