1. How consistent are we? Interlaboratory comparison study in fathead minnows using the model estrogen 17α-ethinylestradiol to develop recommendations for environmental transcriptomics.
- Author
-
Feswick A, Isaacs M, Biales A, Flick RW, Bencic DC, Wang RL, Vulpe C, Brown-Augustine M, Loguinov A, Falciani F, Antczak P, Herbert J, Brown L, Denslow ND, Kroll KJ, Lavelle C, Dang V, Escalon L, Garcia-Reyero N, Martyniuk CJ, and Munkittrick KR
- Subjects
- Animals, Cyprinidae metabolism, Enzyme-Linked Immunosorbent Assay, Liver drug effects, Male, Models, Chemical, Oligonucleotide Array Sequence Analysis, RNA isolation & purification, RNA metabolism, Vitellogenins blood, Cyprinidae genetics, Endocrine Disruptors toxicity, Ethinyl Estradiol toxicity, Laboratories standards, Liver metabolism, Transcriptome drug effects
- Abstract
Fundamental questions remain about the application of omics in environmental risk assessments, such as the consistency of data across laboratories. The objective of the present study was to determine the congruence of transcript data across 6 independent laboratories. Male fathead minnows were exposed to a measured concentration of 15.8 ng/L 17α-ethinylestradiol (EE2) for 96 h. Livers were divided equally and sent to the participating laboratories for transcriptomic analysis using the same fathead minnow microarray. Each laboratory was free to apply bioinformatics pipelines of its choice. There were 12 491 transcripts that were identified by one or more of the laboratories as responsive to EE2. Of these, 587 transcripts (4.7%) were detected by all laboratories. Mean overlap for differentially expressed genes among laboratories was approximately 50%, which improved to approximately 59.0% using a standardized analysis pipeline. The dynamic range of fold change estimates was variable between laboratories, but ranking transcripts by their relative fold difference resulted in a positive relationship for comparisons between any 2 laboratories (mean R
2 > 0.9, p < 0.001). Ten estrogen-responsive genes encompassing a fold change range from dramatic (>20-fold; e.g., vitellogenin) to subtle (∼2-fold; i.e., block of proliferation 1) were identified as differentially expressed, suggesting that laboratories can consistently identify transcripts that are known a priori to be perturbed by a chemical stressor. Thus, attention should turn toward identifying core transcriptional networks using focused arrays for specific chemicals. In addition, agreed-on bioinformatics pipelines and the ranking of genes based on fold change (as opposed to p value) should be considered in environmental risk assessment. These recommendations are expected to improve comparisons across laboratories and advance the use of omics in regulations. Environ Toxicol Chem 2017;36:2593-2601. © 2017 SETAC., (© 2017 SETAC.)- Published
- 2017
- Full Text
- View/download PDF