Empirical Characterization of False Discovery Rates of Differentially Expressed Genes and Transcriptomic Benchmark Concentrations in Zebrafish Embryos

Abstract

High-throughput transcriptomics (HTTr) is increasingly applied to zebrafish embryos to survey the toxicological effects of environmental chemicals. Before the adoption of this approach in regulatory testing, it is essential to characterize background noise in order to guide experimental designs. We thus empirically quantified the HTTr false discovery rate (FDR) across different embryo pool sizes, sample sizes, and concentration groups for toxicology studies. We exposed zebrafish embryos to 0.1% dimethyl sulfoxide (DMSO) for 5 days. Pools of 1, 5, 10, and 20 embryos were created (n = 24 samples for each pool size). Samples were sequenced on the TempO-Seq platform and then randomly assigned to mock treatment groups before differentially expressed gene (DEG), pathway, and benchmark concentration (BMC) analyses. Given that all samples were treated with DMSO, any significant DEGs, pathways, or BMCs are false positives. As expected, we found decreasing FDRs for DEG and pathway analyses with increasing pool and sample sizes. Similarly, FDRs for BMC analyses decreased with increasing pool size and concentration groups, with more stringent BMC premodel filtering reducing BMC FDRs. Our study provides foundational data for determining appropriate experiment designs for regulatory toxicity testing with HTTr in zebrafish embryos.