Improving Inference from Reported Concentrations in Environmental Surveillance by Modeling the Statistical Features of Digital PCR
7 affiliations
- 1Department of Biosystems Science and Engineering
- 2ETH Zurich
- 3SIB Swiss Institute of Bioinformatics
- 4Eawag
- 5Swiss Federal Institute of Aquatic Science and Technology
- 6Swiss Tropical and Public Health Institute
- 7University of Basel
Abstract
Digital polymerase chain reaction (dPCR) is a powerful technique for quantifying gene targets in environmental samples, with applications ranging from species monitoring to wastewater-based epidemiology. Although accurate statistical models for analyzing dPCR measurements exist, these require exact assay parameters and the number of positive and total PCR partitions. In practice, however, many environmental studies and monitoring programs analyze only concentration estimates, assuming normally or log-normally distributed measurements. Such assumptions ignore key statistical features of PCR assays, including concentration-dependent measurement noise and nondetects, leading to biased environmental estimates. In this work, we present a Bayesian model with a dPCR-specific likelihood that can be fitted directly to reported concentrations, while incorporating uncertainty in assay parameters through interpretable priors. Using real-world case studies of free-eDNA decay in seawater and pathogen transmission from wastewater, we show that our approach yields similar estimates as a fully informed model with partition counts, while avoiding biases from normal or log-normal approximations. This enables accurate inference from dPCR measurements even when partition count data and assay parameters are unavailable. The method is implemented in the R packages "dPCRfit" for regression analyses and "EpiSewer" for wastewater surveillance.
Related
- Efficacy of statistical and design interventions for improving reproducibility in biomedical research
- Serum NT-proBNP Concentrations in Pediatric Patients with a History of Idiopathic Cardiac Arrhythmias: A Pilot Study
- Sex differences dominate urinary extracellular DNA concentrations and variability in humans and rodents
- Accurate quantification of protein concentrations and state populations by 2D 1H-N Q-BEST NMR
- Extreme Slowdown of dsDNA and RNA:DNA Hybrids in Solid-state Nanopores for LiCl concentrations above 4 mol/L
BibTeX
@article{lison2026improving,
title = {{Improving Inference from Reported Concentrations in Environmental Surveillance by Modeling the Statistical Features of Digital PCR}},
author = {Adrian Lison and Timothy Ryan Julian and Tanja Stadler},
journal = {ACS ES\&T Water},
year = {2026},
month = jun,
doi = {10.1021/acsestwater.5c01051},
url = {https://doi.org/10.1021/acsestwater.5c01051},
note = {CC BY},
}
RIS
TY - JOUR
TI - Improving Inference from Reported Concentrations in Environmental Surveillance by Modeling the Statistical Features of Digital PCR
AU - Adrian Lison
AU - Timothy Ryan Julian
AU - Tanja Stadler
PY - 2026
DA - 2026/06/10
JO - ACS ES&T Water
DO - 10.1021/acsestwater.5c01051
UR - https://doi.org/10.1021/acsestwater.5c01051
AB - Digital polymerase chain reaction (dPCR) is a powerful technique for quantifying gene targets in environmental samples, with applications ranging from species monitoring to wastewater-based epidemiology. Although accurate statistical models for analyzing dPCR measurements exist, these require exact assay parameters and the number of positive and total PCR partitions. In practice, however, many environmental studies and monitoring programs analyze only concentration estimates, assuming normally or log-normally distributed measurements. Such assumptions ignore key statistical features of PCR assays, including concentration-dependent measurement noise and nondetects, leading to biased environmental estimates. In this work, we present a Bayesian model with a dPCR-specific likelihood that can be fitted directly to reported concentrations, while incorporating uncertainty in assay parameters through interpretable priors. Using real-world case studies of free-eDNA decay in seawater and pathogen transmission from wastewater, we show that our approach yields similar estimates as a fully informed model with partition counts, while avoiding biases from normal or log-normal approximations. This enables accurate inference from dPCR measurements even when partition count data and assay parameters are unavailable. The method is implemented in the R packages "dPCRfit" for regression analyses and "EpiSewer" for wastewater surveillance.
SN - 2690-0637
ER -