{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,23]],"date-time":"2026-10-23T20:04:31Z","timestamp":1792785871975,"version":"4.5.0"},"reference-count":42,"publisher":"Oxford University Press (OUP)","license":[{"start":{"date-parts":[[2026,9,21]],"date-time":"2026-09-21T00:00:00Z","timestamp":1789948800000},"content-version":"am","delay-in-days":2,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"Dutch Research Council (NWO) Domain Applied and Engineering Sciences","award":["15494"],"award-info":[{"award-number":["15494"]}],"id":[{"id":"https:\/\/ror.org\/04jsz6e67","id-type":"ROR","asserted-by":"crossref"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>Biomarkers can provide a high-throughput and accurate assessment of the impact of complex chemical mixtures in the environment on organisms, but their identification through gene expression analysis is hindered by noise, synergistic interactions, and nonlinear expression patterns. We generated finely resolved transcriptomic data from the ecotoxicological model species Folsomia candida exposed to two binary pesticide mixtures: One combining two neonicotinoid insecticides (imidacloprid and clothianidin) and the other a neonicotinoid (imidacloprid) with an azole fungicide (cyproconazole). Using these datasets, we developed a Gaussian process (GP) framework to identify robust gene expression biomarkers, accounting for nonlinear and synergistic interaction effects across experiments. Joint analysis of two binary mixtures increased the overlap of differentially expressed genes compared to separate analyses, improving robustness. In simulations, GP models outperformed linear models, accurately fitting complex, nonlinear concentration\u2013response relationships. Four biomarkers, three for neonicotinoids (ARRD, SMCT, and nAchR) and one for azole fungicides (CYP), identified through this framework, were empirically validated and confirmed to be specifically responsive to their target pesticide, even under co-exposure. These findings highlight the effectiveness of GP models for mixture exposure transcriptomics and their broader applicability to other omics data and research fields.<\/jats:p>","DOI":"10.1093\/etojnl\/vgag252","type":"journal-article","created":{"date-parts":[[2026,9,21]],"date-time":"2026-09-21T11:41:22Z","timestamp":1789990882000},"source":"Crossref","is-referenced-by-count":0,"title":["A Gaussian process approach facilitates the identification of robust biomarkers for exposure to complex pesticide mixtures"],"prefix":"10.1093","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-6083-3675","authenticated-orcid":false,"given":"Ruben","family":"Bakker","sequence":"first","affiliation":[{"name":"Faculty of Science, Amsterdam Institute for Life and Environment, Vrije Universiteit Amsterdam , Amsterdam,","place":["The Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-4487-5947","authenticated-orcid":false,"given":"Yuliya","family":"Shapovalova","sequence":"additional","affiliation":[{"name":"Department of Data Science, Institute for Computing and Information Sciences, Radboud Universiteit Nijmegen , Nijmegen,","place":["The Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-4450-701X","authenticated-orcid":false,"given":"Tjeerd M H","family":"Dijkstra","sequence":"additional","affiliation":[{"name":"Institute for Translational Bioinformatics, University Hospital T\u00fcbingen , T\u00fcbingen,","place":["Germany"]},{"name":"Max Planck Institute for Biology , T\u00fcbingen,","place":["Germany"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-3398-5235","authenticated-orcid":false,"given":"Tom","family":"Heskes","sequence":"additional","affiliation":[{"name":"Department of Data Science, Institute for Computing and Information Sciences, Radboud Universiteit Nijmegen , Nijmegen,","place":["The Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5651-0208","authenticated-orcid":false,"given":"Cornelis A M","family":"van Gestel","sequence":"additional","affiliation":[{"name":"Faculty of Science, Amsterdam Institute for Life and Environment, Vrije Universiteit Amsterdam , Amsterdam,","place":["The Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6090-232X","authenticated-orcid":false,"given":"Katja M","family":"Hoedjes","sequence":"additional","affiliation":[{"name":"Faculty of Science, Amsterdam Institute for Life and Environment, Vrije Universiteit Amsterdam , Amsterdam,","place":["The Netherlands"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,9,19]]},"reference":[{"key":"2026102309430672300_vgag252-B1","doi-asserted-by":"publisher","first-page":"2508","DOI":"10.1021\/es2038036","article-title":"Mixture toxicity revisited from a toxicogenomic perspective","volume":"46","author":"Altenburger","year":"2012","journal-title":"Environmental Science & Technology"},{"key":"2026102309430672300_vgag252-B2","doi-asserted-by":"publisher","first-page":"80897","DOI":"10.1007\/s11356-022-21362-z","article-title":"Biomarker development for neonicotinoid exposure in soil under interaction with the synergist piperonyl butoxide in Folsomia candida","volume":"29","author":"Bakker","year":"2022","journal-title":"Environmental Science and Pollution Research International"},{"key":"2026102309430672300_vgag252-B3","doi-asserted-by":"publisher","first-page":"161740","DOI":"10.1016\/j.scitotenv.2023.161740","article-title":"Combining time-resolved transcriptomics and proteomics data for adverse outcome pathway refinement in ecotoxicology","volume":"869","author":"Bakker","year":"2023","journal-title":"Science of the Total Environment"},{"key":"2026102309430672300_vgag252-B4","doi-asserted-by":"publisher","DOI":"10.1007\/978-3-031-78987-8","volume-title":"Statistical methods for environmental mixtures: A primer in environmental epidemiology","author":"Bellavia","year":"2025"},{"key":"2026102309430672300_vgag252-B5","doi-asserted-by":"publisher","first-page":"67","DOI":"10.1186\/s12940-018-0413-y","article-title":"Statistical software for analyzing the health effects of multiple concurrent exposures via Bayesian kernel machine regression","volume":"17","author":"Bobb","year":"2018","journal-title":"Environmental Health"},{"key":"2026102309430672300_vgag252-B6","doi-asserted-by":"publisher","first-page":"252","DOI":"10.1093\/toxsci\/kfx097","article-title":"The role of omics in the application of adverse outcome pathways for chemical risk assessment","volume":"158","author":"Brockmeier","year":"2017","journal-title":"Toxicological Sciences"},{"key":"2026102309430672300_vgag252-B7460580","doi-asserted-by":"publisher","first-page":"67","DOI":"10.1584\/jpestics.D17-01930363948","article-title":"The invertebrate pharmacology of insecticides acting at nicotinic acetylcholine receptors","volume":"42","author":"Crossthwaite","year":"2017","journal-title":"Journal of Pesticide Science"},{"key":"2026102309430672300_vgag252-B7","doi-asserted-by":"publisher","first-page":"D745","DOI":"10.1093\/nar\/gky1113","article-title":"Ensembl 2019","volume":"47","author":"Cunningham","year":"2019","journal-title":"Nucleic Acids Research"},{"key":"2026102309430672300_vgag252-B8","doi-asserted-by":"publisher","first-page":"11","DOI":"10.1186\/1471-2199-12-11","article-title":"High throughput nano-liter RT-qPCR to classify soil contamination using a soil arthropod","volume":"12","author":"De Boer","year":"2011","journal-title":"BMC Molecular Biology"},{"key":"2026102309430672300_vgag252-B9","doi-asserted-by":"publisher","first-page":"54","DOI":"10.1186\/1471-2199-10-54","article-title":"Reference genes for QRT-PCR tested under various stress conditions in Folsomia candida and Orchesella cincta (Insecta, Collembola)","volume":"10","author":"De Boer","year":"2009","journal-title":"BMC Molecular Biology"},{"key":"2026102309430672300_vgag252-B10","doi-asserted-by":"publisher","first-page":"548","DOI":"10.1002\/etc.4634","article-title":"Toxicity in neonicotinoids to Folsomia candida and Eisenia andrei","volume":"39","author":"De Lima E Silva","year":"2020","journal-title":"Environmental Toxicology and Chemistry"},{"key":"2026102309430672300_vgag252-B11","doi-asserted-by":"publisher","first-page":"e52707","DOI":"10.7554\/eLife.52707","article-title":"The transcriptomic response of cells to a drug combination is more than the sum of the responses to the monotherapies","volume":"9","author":"Diaz","year":"2020","journal-title":"eLife"},{"key":"2026102309430672300_vgag252-B12","doi-asserted-by":"publisher","first-page":"1184","DOI":"10.1038\/nprot.2009.97","article-title":"Mapping identifiers for the integration of genomic datasets with the R\/Bioconductor package biomaRt","volume":"4","author":"Durinck","year":"2009","journal-title":"Nature Protocols"},{"key":"2026102309430672300_vgag252-B13","doi-asserted-by":"publisher","first-page":"388","DOI":"10.1126\/science.aay6636","article-title":"Tracking complex mixtures of chemicals in our changing environment","volume":"367","author":"Escher","year":"2020","journal-title":"Science"},{"key":"2026102309430672300_vgag252-B14","doi-asserted-by":"publisher","first-page":"3047","DOI":"10.1093\/bioinformatics\/btw354","article-title":"MultiQC: Summarize analysis results for multiple tools and samples in a single report","volume":"32","author":"Ewels","year":"2016","journal-title":"Bioinformatics"},{"key":"2026102309430672300_vgag252-B15","doi-asserted-by":"publisher","first-page":"493","DOI":"10.1186\/s12864-017-3852-x","article-title":"Coping with living in the soil: The genome of the parthenogenetic springtail Folsomia candida","volume":"18","author":"Faddeeva-Vakhrusheva","year":"2017","journal-title":"BMC Genomics"},{"key":"2026102309430672300_vgag252-B16","doi-asserted-by":"publisher","first-page":"201","DOI":"10.1146\/annurev.ento.50.071803.130331","article-title":"Folsomia candida (COLLEMBOLA): A \u201cStandard\u201d Soil Arthropod","volume":"50","author":"Fountain","year":"2005","journal-title":"Annual Review of Entomology"},{"key":"2026102309430672300_vgag252-B17","doi-asserted-by":"publisher","first-page":"180","DOI":"10.1186\/1471-2105-12-180","article-title":"A simple approach to ranking differentially expressed gene expression time courses through Gaussian process regression","volume":"12","author":"Kalaitzis","year":"2011","journal-title":"BMC Bioinformatics"},{"key":"2026102309430672300_vgag252-B18","doi-asserted-by":"publisher","first-page":"14461","DOI":"10.1021\/acs.est.8b04752","article-title":"DRomics: A Turnkey tool to support the use of the dose\u2013response framework for omics data in ecological risk assessment","volume":"52","author":"Larras","year":"2018","journal-title":"Environmental Science & Technology"},{"key":"2026102309430672300_vgag252-B19","doi-asserted-by":"publisher","first-page":"1340","DOI":"10.3390\/biology10121340","article-title":"Biomarker studies in stress biology: From the gene to population, from the organism to the application","volume":"10","author":"Lemos","year":"2021","journal-title":"Biology"},{"key":"2026102309430672300_vgag252-B20","doi-asserted-by":"publisher","first-page":"1","DOI":"10.2174\/1875318301909010001","article-title":"Pollution biomarkers in environmental and human biomonitoring","volume":"9","author":"Lionetto","year":"2019","journal-title":"Open Biomarkers Journal"},{"key":"2026102309430672300_vgag252-B21","doi-asserted-by":"publisher","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biology"},{"key":"2026102309430672300_vgag252-B22","doi-asserted-by":"publisher","first-page":"143263","DOI":"10.1016\/j.scitotenv.2020.143263","article-title":"Imidacloprid dominates the combined toxicities of neonicotinoid mixtures to stream mayfly nymphs","volume":"761","author":"Macaulay","year":"2021","journal-title":"Science of the Total Environment"},{"key":"2026102309430672300_vgag252-B23","doi-asserted-by":"publisher","first-page":"10","DOI":"10.14806\/ej.17.1.200","article-title":"Cutadapt removes adapter sequences from high-throughput sequencing reads","volume":"17","author":"Martin","year":"2011","journal-title":"EMBnet.journal"},{"key":"2026102309430672300_vgag252-B24","doi-asserted-by":"publisher","first-page":"605","DOI":"10.1038\/embor.2010.80","article-title":"Arrestin domain-containing protein 3 recruits the NEDD4 E3 ligase to mediate ubiquitination of the \u03b22-adrenergic receptor","volume":"11","author":"Nabhan","year":"2010","journal-title":"EMBO Reports"},{"key":"2026102309430672300_vgag252-B25","doi-asserted-by":"publisher","first-page":"1155","DOI":"10.1158\/1535-7163.MCT-15-0843","article-title":"An unbiased oncology compound screen to identify novel combination strategies","volume":"15","author":"O\u2019Neil","year":"2016","journal-title":"Molecular Cancer Therapeutics"},{"key":"2026102309430672300_vgag252-B26","doi-asserted-by":"publisher","first-page":"417","DOI":"10.1038\/nmeth.4197","article-title":"Salmon provides fast and bias-aware quantification of transcript expression","volume":"14","author":"Patro","year":"2017","journal-title":"Nature Methods"},{"key":"2026102309430672300_vgag252-B27","doi-asserted-by":"publisher","first-page":"107167","DOI":"10.1016\/j.agee.2020.107167","article-title":"Residues of currently used pesticides in soils and earthworms: A silent threat?","volume":"305","author":"Pelosi","year":"2021","journal-title":"Agriculture, Ecosystems & Environment"},{"key":"2026102309430672300_vgag252-B28","doi-asserted-by":"publisher","first-page":"68","DOI":"10.1007\/s11356-014-3471-x","article-title":"Effects of neonicotinoids and fipronil on non-target invertebrates","volume":"22","author":"Pisa","year":"2015","journal-title":"Environmental Science and Pollution Research International"},{"key":"2026102309430672300_vgag252-B29","volume-title":"Gaussian processes for machine learning (3. print)","author":"Rasmussen","year":"2008"},{"key":"2026102309430672300_vgag252-B30","volume-title":"R: A language and environment for statistical computing (Version 4.0.0) [Computer software]","author":"R Core Team","year":"2020"},{"key":"2026102309430672300_vgag252-B31","doi-asserted-by":"publisher","first-page":"808335","DOI":"10.3389\/finsc.2021.808335","article-title":"Interaction of insecticides and fungicides in bees","volume":"1","author":"Schuhmann","year":"2021","journal-title":"Frontiers in Insect Science"},{"key":"2026102309430672300_vgag252-B32","doi-asserted-by":"publisher","first-page":"47006","DOI":"10.1289\/EHP7773","article-title":"Transcriptome-wide prediction and measurement of combined effects induced by chemical mixture exposure in zebrafish embryos","volume":"129","author":"Sch\u00fcttler","year":"2021","journal-title":"Environmental Health Perspectives"},{"key":"2026102309430672300_vgag252-B33","doi-asserted-by":"publisher","first-page":"14","DOI":"10.1186\/s12859-021-04508-7","article-title":"Non-parametric synergy modeling of chemical compounds with Gaussian processes","volume":"23","author":"Shapovalova","year":"2022","journal-title":"BMC Bioinformatics"},{"key":"2026102309430672300_vgag252-B34","doi-asserted-by":"publisher","first-page":"108280","DOI":"10.1016\/j.envint.2023.108280","article-title":"Pesticide residues with hazard classifications relevant to non-target species including humans are omnipresent in the environment and farmer residences","volume":"181","author":"Silva","year":"2023","journal-title":"Environment International"},{"key":"2026102309430672300_vgag252-B35","doi-asserted-by":"publisher","first-page":"3002","DOI":"10.1038\/s41598-019-38907-x","article-title":"Comparison of null models for combination drug therapy reveals hand model as biochemically most plausible","volume":"9","author":"Sinzger","year":"2019","journal-title":"Scientific Reports"},{"key":"2026102309430672300_vgag252-B36","doi-asserted-by":"publisher","first-page":"1521","DOI":"10.12688\/f1000research.7563.1","article-title":"Differential analyses for RNA-seq: Transcript-level estimates improve gene-level inferences","volume":"4","author":"Soneson","year":"2015","journal-title":"F1000Research"},{"key":"2026102309430672300_vgag252-B37","doi-asserted-by":"publisher","first-page":"279","DOI":"10.1186\/s12859-022-04775-y","article-title":"Batch effect detection and correction in RNA-seq data using machine-learning-based automated assessment of quality","volume":"23","author":"Sprang","year":"2022","journal-title":"BMC Bioinformatics"},{"key":"2026102309430672300_vgag252-B38","doi-asserted-by":"publisher","first-page":"2124","DOI":"10.1002\/etc.5402","article-title":"Soil ecotoxicology needs robust biomarkers: A meta-analysis approach to test the robustness of gene expression-based biomarkers for measuring chemical exposure effects in soil invertebrates","volume":"41","author":"Swart","year":"2022","journal-title":"Environmental Toxicology and Chemistry"},{"key":"2026102309430672300_vgag252-B39","doi-asserted-by":"publisher","first-page":"452","DOI":"10.1002\/etc.5402","article-title":"GNU parallel: The command-line power tool","volume":"36","author":"Tange","year":"2011","journal-title":"Usenix\u2009&lt;\u2009ag"},{"key":"2026102309430672300_vgag252-B40","doi-asserted-by":"publisher","first-page":"4607","DOI":"10.1038\/s41467-021-24789-z","article-title":"MuSyC is a consensus framework that unifies multi-drug synergy metrics for combinatorial drug discovery","volume":"12","author":"Wooten","year":"2021","journal-title":"Nature Communications"},{"key":"2026102309430672300_vgag252-B41","doi-asserted-by":"publisher","first-page":"R14","DOI":"10.1186\/gb-2010-11-2-r14","article-title":"Gene ontology analysis for RNA-seq: Accounting for selection bias","volume":"11","author":"Young","year":"2010","journal-title":"Genome Biology"}],"container-title":["Environmental Toxicology and Chemistry"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/etc\/advance-article-pdf\/doi\/10.1093\/etojnl\/vgag252\/71254139\/vgag252.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/etc\/advance-article-pdf\/doi\/10.1093\/etojnl\/vgag252\/71254139\/vgag252.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,10,23]],"date-time":"2026-10-23T13:43:19Z","timestamp":1792762999000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/etc\/advance-article\/doi\/10.1093\/etojnl\/vgag252\/8820724"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,9,19]]},"references-count":42,"URL":"https:\/\/doi.org\/10.1093\/etojnl\/vgag252","relation":{},"ISSN":["0730-7268","1552-8618"],"issn-type":[{"value":"0730-7268","type":"print"},{"value":"1552-8618","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,9,19]]},"article-number":"vgag252"}}