Systematic review on the laboratory methodology for conducting wastewater and environmental surveillance for Salmonella
Lucky Sangal, Vishesh Sood, Karin Haar, Takana Mary Silubonde, Yuka Jinnai, Suman Rijal

TL;DR
This systematic review evaluates lab methods for detecting Salmonella in wastewater and surface waters, highlighting the need for standardized protocols to improve surveillance.
Contribution
The study identifies six methodological pathways and highlights the lack of standardized protocols and quality control in Salmonella wastewater surveillance.
Findings
Grab sampling is the most commonly used method for wastewater and surface water sampling.
Fewer than 14% of studies reported comprehensive quality control measures.
Six distinct methodological pathways for Salmonella detection were identified.
Abstract
Wastewater and environmental surveillance (WES) is a valuable supplementary tool to clinical surveillance for infectious diseases, especially in low- and middle-income countries. This systematic review evaluates laboratory methods for detecting Salmonella spp. in wastewater and contaminated surface waters, focusing on methodological diversity, feasibility, and the need for standardized protocols. The review was performed using protocol registered with PROSPERO (ID: CRD42024573052) following PRISMA 2020 guidelines. The review was funded by the European Commission’s Health Emergency Preparedness and Response Authority (HERA) and the World Health Organization (WHO). Search was conducted in PubMed, EMBASE, and Web of Science (last update: May 31, 2025). Studies describing sampling and laboratory methods for Salmonella detection in wastewater or contaminated surface waters were included.…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6| Paper ID (Reference) | Sample | Processing | Culture | Enumeration | Biotyping | Serotyping | AST | Pheno. other | Bacteriophage | Molecular | Mol. ARG | Mol. other | Sequencing |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Paper_001 ( | G | ✓ | ✓ | ||||||||||
| T | ✓ | ✓ | |||||||||||
| Paper_002 ( | T | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_003 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_004 ( | G | ✓ | ✓ | ||||||||||
| Paper_005 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_006 ( | T | ✓ | ✓ | ✓ | |||||||||
| Paper_007 ( | T | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_008 ( | C | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_009 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_010 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_011 ( | T | ✓ | ✓ | ✓ | |||||||||
| Paper_012 ( | G | ✓ | ✓ | ||||||||||
| Paper_013 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_014 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||
| Paper_015 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_016 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_017 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_018 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_019 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_020 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_021 ( | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_022 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_023 ( | C | ✓ | ✓ | ||||||||||
| Paper_024 ( | C | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_025 ( | C | ✓ | ✓ | ||||||||||
| Paper_026 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_027 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_028 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_029 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_030 ( | G | ✓ | ✓ | ||||||||||
| Paper_031 ( | G | ✓ | ✓ | ||||||||||
| Paper_032 ( | T | ✓ | ✓ | ||||||||||
| Paper_033 ( | G | ✓ | ✓ | ||||||||||
| Paper_034 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_035 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_036 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||
| Paper_037 ( | ✓ | ✓ | ✓ | ||||||||||
| Paper_038 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_039 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_040 ( | G | ✓ | ✓ | ||||||||||
| Paper_041 ( | G | ✓ | ✓ | ||||||||||
| Paper_042 ( | G | ✓ | |||||||||||
| Paper_043 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_044 ( | S | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_045 ( | C | ✓ | ✓ | ||||||||||
| Paper_046 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_047 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_048 ( | T | ✓ | ✓ | ||||||||||
| G | ✓ | ✓ | |||||||||||
| Paper_049 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_050 ( | G | ✓ | ✓ | ||||||||||
| Paper_051 ( | T | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_052 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_053 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_054 ( | C | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_055 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_056 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_057 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_058 ( | C | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_059 ( | G | ✓ | ✓ | ||||||||||
| T | ✓ | ||||||||||||
| Paper_060 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_061 ( | G | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_062 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_063 ( | G | ✓ | ✓ | ||||||||||
| T | ✓ | ✓ | |||||||||||
| Paper_064 ( | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_065 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_066 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| T | ✓ | ✓ | ✓ | ✓ | |||||||||
| Paper_067 ( | G | ✓ | ✓ | ||||||||||
| Paper_068 ( | G | ✓ | ✓ | ✓ | |||||||||
| G | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_069 ( | C | ✓ | ✓ | ||||||||||
| Paper_070 ( | T | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_071 ( | T | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_072 ( | T | ✓ | ✓ | ✓ | ✓ | ||||||||
| Paper_073 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_074 ( | G | ✓ | ✓ | ||||||||||
| Paper_075 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_076 ( | G | ✓ | ✓ | ||||||||||
| Paper_077 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_078 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_079 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | |||||
| Paper_080 ( | G | ✓ | ✓ | ||||||||||
| Paper_081 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | |||||
| Paper_082 ( | G | ✓ | ✓ | ||||||||||
| Paper_083 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | |||||||
| Paper_084 ( | ✓ | ✓ | ✓ | ||||||||||
| Paper_85 ( | G | ✓ | ✓ | ||||||||||
| T | ✓ | ✓ | ✓ | ||||||||||
| Paper_86 ( | G | ✓ | ✓ | ||||||||||
| T | ✓ | ✓ | ✓ | ||||||||||
| Paper_87 ( | T | ✓ | |||||||||||
| Paper_88 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_89 ( | G | ✓ | ✓ | ✓ | |||||||||
| Paper_90 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_91 ( | ✓ | ✓ | |||||||||||
| Paper_92 ( | C | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ||||||
| Paper_093 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | |||||
| Paper_094 ( | G | ✓ | ✓ | ✓ | ✓ | ✓ |
| Public health scope | Domains | Output type | Testing pathways |
|---|---|---|---|
|
Monitoring for the detection of importations Monitoring of community-level baselines Monitoring of changes in risk factors Monitoring of epidemiological changes Monitoring the effect of changes in healthcare practices Optimization of resource and budget allocation Evidence generation for developing and implementing public health policies Evidence generation for evaluation of public health policies Evidence generation for calibration of public health interventions |
Outbreak detection or identification (Niche) Disease Prevalence (Core) AMR Prevalence (Core) and/or Mechanisms (Supporting) Monitoring wastewater usage (Supporting) Environmental Health (Core) One Health (Cross-cutting) Analytical Method Validation (Cross-cutting) |
Cross-sectional (single site, single time) Time-series (single site, multiple times) Longitudinal (multiple sites, multiple times) Spatial (multiple sites) Spatio-temporal (multiple sites, multiple times) Hierarchical (nested structure of sites, e.g., small to large drains) Network (interconnected site, e.g., multiple drains linked to WWTP) |
P1 P2 P3 P4 P5 P6 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsSARS-CoV-2 detection and testing · Fecal contamination and water quality · Salmonella and Campylobacter epidemiology
Introduction
1
Wastewater and environmental surveillance (WES) is well-established for poliovirus surveillance as a core part of the Global Polio Eradication Initiative (1–4), and gained further prominence during the COVID-19 pandemic as a valuable tool for understanding the disease burden in communities (5) The success of WES in COVID-19 has sparked interest in its use to monitor other pathogens of public health concern (6). Prioritization of pathogens for WES remains a fundamental question for optimizing resource utilization and ensuring operational flexibility and adaptability in the event of outbreaks caused by new pathogens (7–9). Proposed prioritization frameworks emphasize several key factors for successful WES adaptation, including the public health significance of the pathogen, the usefulness of WES data for public health actions, and the analytical feasibility of conducting WES (8, 9). Prioritizing a pathogen within this framework helps align WES efforts with broader public health goals.
Salmonella infections, particularly with typhoid and paratyphoid serovars, continue to pose a substantial burden on global public health systems, and remain a significant challenge in low- and middle-income countries (LMICs), including within the World Health Organization’s South-East Asia Region (SEAR) (10, 11). Non-typhoid and non-paratyphoid Salmonella serovars. Infect both humans and various animals, making them a significant concern for zoonotic transmission and for both animal husbandry and the food industry (12, 13). Asymptomatic carriers and subclinical infections play a key role in maintaining the transmission chain of Salmonella infections (14, 15). Therefore, accurately determining the true prevalence of Salmonella-related diseases requires supplementing clinical surveillance with serosurveys or contact tracing during outbreaks (16, 17). However, current methods for additional surveillance have limited sensitivity. For instance, the clinical diagnosis of typhoid and paratyphoid often relies on non-specific Widal tests or blood cultures, both of which have low sensitivity due to suboptimal sampling times post-incubation period (18, 19). Despite the World Health Organization not recommending the Widal test, it remains widely used in clinical practice in our region. In addition to its limited sensitivity, the Widal test is prone to cross-reactivity with other pathogens, further reducing its diagnostic specificity. Additionally, estimating HlyE IgG antibodies using ELISA is the preferred method for serosurveys on typhoid and paratyphoid prevalence; however, HlyE antibodies can exhibit cross-reactivity, as many other bacteria also express HlyE (20, 21). As a result, there is a need to establish supplementary tools to accurately estimate the true prevalence of the infections caused by Salmonella Typhi and Paratyphi.
Since Salmonella is present in wastewater due to shedding in the feces of both symptomatic and asymptomatic individuals, depending on the stage of infection, WES has proven effective in evaluating its community burden in endemic countries, complementing existing clinical surveillance and serosurveys (22–24). Besides assessing community-level disease burdens, WES can also generate data on circulating strains and antimicrobial resistance—provided the bacteria can be cultured—which directly informs vaccination and antimicrobial resistance (AMR) strategies, offering significant public health benefits for Salmonella monitoring (25, 26). Therefore, the public health importance of Salmonella and the utility of WES make it a priority pathogen for WES implementation (27).
The primary challenge of Salmonella WES lies in analytical feasibility, due to the heterogeneous and variable factors outside the laboratory that affect sample collection and quality. Unlike high-income countries with centralized and closed sewage systems, most Salmonella Typhi and Paratyphi-endemic countries, particularly in the SEAR region (24, 28, 29), frequently rely on decentralized, informal, or mixed drainage networks, including open drains, septic tanks, and combined sewer-stormwater systems, which are often poorly maintained and vulnerable to contamination during monsoons (30, 31). Thus, the pre-examination factors such as variability in infrastructure, flow dynamics, and ambient conditions complicate sample collection and pathogen recovery and demand context-specific adaptations to sampling and testing protocols. Laboratory capacity constraints, including cold chain logistics and molecular diagnostics standardization and result interpretations, further limit the applicability of WES protocols (28, 29, 32).
Inherent variability in sampling site characteristics, public health goals, and laboratory capacity underscores the urgent need to develop harmonized methodologies that are scientifically sound, operationally practical, and adaptable to diverse infrastructure contexts, particularly in low- and middle-income countries. This systematic review was conducted to evaluate the scientific and operational feasibility of laboratory methods for detecting Salmonella in wastewater, aiming to guide the development of harmonized, context-specific field and laboratory protocols that can support regional public health goals like integrated disease surveillance and inform the deployment of typhoid conjugate vaccine (TCV) in endemic areas and LMICs. The review assessed the completeness of methodological reporting; such as site selection, sample handling, and quality control to identify critical gaps in methodological reporting that might impede reproducibility and scalability. The identified methodologies were further categorized as pathways to match protocol steps with the wastewater sampling and socio-economic status of reporting countries. Finally, the review also identified the primary public health domains that researchers utilize to develop a framework for aligning surveillance objectives with laboratory capacities and infrastructure realities.
Methods
2
Study design
2.1
This qualitative systematic review was designed and reported in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) 2020 guidelines (33). A detailed protocol outlining the objectives, eligibility criteria, and methodological approach was developed before the initiation of the review and registered with the International Prospective Register of Systematic Reviews (PROSPERO) on August 5, 2024 (registration ID: CRD42024573052). The PICOS framework was adapted to develop a research question and inclusion and exclusion criteria were defined to align with PICOS framework (Supplementary Table 7). Based on the PICOS framework, the objective of this systematic review was to map the diversity in the laboratory workflows and thematic domains used for the detection of Salmonella spp. from wastewater and wastewater-impacted surface waters, and to synthesize these findings into evidence-based heuristic framework for resource planning. Further, the systematic review also aimed at evaluating the completeness of the reporting of laboratory methodologies and work on the recommendations to ensure reporting of reproducible methodology.
Search strategy
2.2
A comprehensive literature search was conducted to identify peer-reviewed studies describing laboratory methodologies for the detection or isolation of Salmonella spp. From wastewater and wastewater contaminated surface waters. The scope of methodologies was kept at species levels as isolation and characterization at species level can be extended to specific sub-species. The initial search was conducted on September 10, 2024, using three major scientific databases: PubMed, EMBASE, and Web of Science. It was updated on May 31, 2025, to include recent publications. The specific combination of databases was chosen because it achieves a 90–95% recall rate in over 80% of reported systematic reviews (34). Google Scholar was not included in the search strategy to improve the recall rate further as Google search uses proprietary algorithms that can result in variation in search results based on time, place, and person, thereby making Google Scholar an unreliable tool when repeatability is required (35, 36). The search strategy was designed to retrieve studies relevant to Salmonella Typhi surveillance in environmental matrices, with a focus on wastewater and surface waters. Briefly, a naïve search was conducted on PubMed after identifying relevant MESH and MAJR terms for WES of Salmonella spp. The query used was - “Salmonella”[Mesh] AND (“Wastewater-Based Epidemiological Monitoring” [Mesh] OR “Environmental Monitoring”[Mesh] OR “Sewage/microbiology”[MAJR] OR “Wastewater/microbiology”[MAJR]). The easyPubMed package in R was used to import the naïve search results, which were analyzed with the litsearchr package in R to identify keywords in an unbiased manner (Supplementary Figure 1; Supplementary Table 3) (37, 38). The identified keywords were used to create PubMed search queries using Boolean operators and wildcard symbols (e.g., *). The search query was optimized to evaluate its sensitivity against a benchmark set of 16 studies on Salmonella spp. (Supplementary Table 4). The search terms were refined until all 16 benchmark studies were captured (Supplementary Figure 2; Supplementary Table 5). Once 100% sensitivity was achieved for the benchmark studies, the final PubMed query was translated into EMBASE and Web of Science formats using the polyglot application (39). The final search queries for each database are provided in (Supplementary Table 6).
Systematic review process
2.3
The final search was conducted in May 2025, and results were imported into the Covidence platform, which automatically removed duplicates. Remaining duplicates were manually reviewed and excluded. Title and abstract screening was then performed, with inclusion and exclusion criteria detailed in Supplementary Table 7. The proportional agreement for the screening step by the authors is provided in Supplementary Table 13. Two authors independently reviewed the studies, resolving conflicts by consensus among all authors. Full texts were retrieved for eligible studies and screened again against the criteria. The Covidence platform was also used to prepare templates for data extraction and completness of reported methodolgy assessment, as described in Supplementary Tables 8, 9. Data extraction was performed independently by two authors, with disagreements resolved by consensus. Supplementary Tables 1, 2 present the 2020 PRISMA checklist and PRISMA Abstract checklist for this systematic review, respectively.
Data analysis and visualization
2.4
A narrative synthesis was conducted due to the variability in study designs, sampling strategies, and laboratory methods. Data visualization involved subgroup analysis based on the methodological areas (e.g., sampling, processing, testing), and patterns were identified across the studies Further, an exploratory analysis was performed to understand commonalities between reported methods to identify major pathways and to identify thematic domains and their interconnectedness in the study dataset of selected manuscripts.
Python (version 3.12.7) was used within the Spyder IDE (version 6.0.7) to perform data analysis and visualization, using a collection of specialized libraries. Pandas handled data import and transformation; NumPy supported numerical calculations; and SciPy was used for estimating Jaccard distances and performing chi-squared tests. The silhouette score was calculated with Scikit-learn, while NetworkX enabled network analysis. GeoPandas managed geospatial polygon data for countries, Matplotlib produced standard plots, Seaborn generated heatmaps, and UpSetPlot was used to create UpSet diagrams.
Geographical mapping of studies was conducted using cultural raster map shapefiles obtained from Natural Earth and analyzed with the GeoPandas library in Python. Since some studies reported using multiple methods or samples, the methods were identified to aid analysis. To find commonalities among the methods, pathway analysis was performed. During this process, methods were clustered using the Jaccard distance approach based on factors such as the economic status of countries, sample types, and the presence or absence of specific protocol steps. For simplicity, the LMIC classification included all countries from the LIC (low-income countries), LMIC, and UMIC (upper middle-income countries) categories. The samples analyzed included grab samples, trap samples, and composite samples. Protocol steps included processing, culture, biotyping (using biochemicals and other biotyping techniques), serotyping (via the Kauffman-White scheme and PCR), antimicrobial susceptibility testing, genotyping (using molecular assays and ARGs), and genomics methods, including targeted sequencing, whole-genome sequencing, and metagenomics. For understanding the commonalities in the reported methodologies, the identified methods were clustered. For clustering, qualitative descriptors were converted into a binary matrix. Each unique protocol step was defined as a binary variable with ‘explicit encoding’ of the presence and absence of the step as ‘1’ and ‘0’, respectively. To empirically identify common methodological pathways, a hierarchical clustering on the binary protocol matrix was performed. Jaccard distance was used to measure the dissimilarity between the studies as it prioritizes the presence of shared features while ignoring the shared absence. Average linkage was used to minimize the variance between the clusters. The optimal number of clusters was determined through the silhouette score and the elbow method. The results were reported as the probability of presence of a step in the pathway. To verify that the pathways identified were not a result of poor reporting, a stability analysis was performed with respect to the completeness of methodology reporting (Supplementary Table 10; Supplementary material 2). First, the reporting scores of pathways excluded from clustering was compared with those of clustered methods using a Welch’s t-test to justify their exclusions from the pathway clustering. The stability of clustering of retained methods was tested using the one-way ANOVA to test variability in reporting scores across pathways. The one-way ANNOVA was also repeated at the study level by aggregating the scores by Paper ID. Finally, the internal cluster stability was quantified using a bootstrap approach (n = 1,000 iterations) to calculate the coefficient of variation (CV) of quality scores. To further explore the thematic domains in the selected papers, two authors identified the domains based on the titles, keywords, and abstracts of the selected studies and any disagreement was resolved based on all authors conensus. Eight thematic domains were recognized: (A) outbreak detection/investigation, (B) disease prevalence, (C) AMR prevalence, (D) mechanisms of AMR, (E) wastewater monitoring, (F) environmental health, (G) One Health, and (H) method validation. Some studies encompassed more than one domain. The co-occurrence of domains was calculated using NumPy. To visualize the core domains and their connections to other domains, the co-occurrence matrix was interpreted as an undirected weighted graph with NetworkX. The network comprises nodes representing individual domains, with node size proportional to the number of studies in each domain. Edges indicate co-occurrences between domains, with edge width reflecting the strength of these co-occurrences. Additionally, a force-directed spring layout was used to position strongly related domains based on the study dataset closer together.
Reporting bias
2.5
The review aimed to eliminate bias during both the searching and review stages. To minimize search bias, multiple databases were utilized, and the search string was refined through unbiased keyword selection and by evaluating the search strategy against a set of benchmark studies. For study selection bias, two authors independently reviewed the abstracts during the screening process to determine eligibility for full-text review, and two authors independently conducted the data extraction.
Results
3
Literature search results
3.1
A total of 2,007 articles were identified across PubMed, Embase, and Web of Science. After removing 686 duplicates, 1,321 records were screened by title and abstract, resulting in 1,143 exclusions. Of the 178 full-text articles assessed, 94 met the eligibility criteria and were included in the review (22–24, 26, 28, 29, 32, 40–126). Studies were excluded for reasons such as incomplete methodology, non-peer-reviewed status, or publication in a language other than English. Included studies were further classified by methodological quality assessment into five categories: excellent (n = 20), robust (n = 22), good (n = 22), fair (n = 22), and low (n = 8). The PRISMA workflow for the systematic review is presented in Figure 1. The extracted data and quality assessment data are provided as Supplementary materials 1, 2, respectively. As some studies used multiple samples and multiple methods, the extracted data were further refined to obtain the methods used for each sample in different studies. This resulted in the identification of 102 methods. Table 1 summarizes the sample types and methodology used by the included studies.
PRISMA flow diagram for database searches, deduplication, screening, retrieval, exclusion, inclusion, and quality-based classification of studies on WES for Salmonella spp.
Geographical distribution of reported studies
3.2
The selected studies represented a wide geographic distribution across 36 countries, five of them in SEAR countries (Figure 2A). Most (n = 89) were single-country studies, while five involved multiple countries. According to the World Bank’s income classification, the studies spanned two LICs, ten LMICs, eight UMICs, and 16 high-income countries (HICs). The United States of America and India contributed the highest number of studies published after 2020 (Figure 2B).
(A) Geographical distribution of countries reporting environmental surveillance of Salmonella spp. (n = 94); (B) Number of studies from different countries over the years 1980 to 2020 in 10-year intervals, and after 2020, along with the socioeconomic status. LIC: Low-income countries, LMIC: low- and middle-income countries, UMIC: upper middle-income countries, and HIC: high-income countries.
Reported laboratory methods
3.3
The included studies employed a range of sampling methods to collect wastewater, each reflecting different operational contexts and surveillance goals (Table 1). Most studies reported using grab sampling (n = 62), which involves collecting a predefined volume of wastewater in a sterile container at a single point in time. This was followed by trap sampling (n = 10), a passive technique in which a receptacle, often a Moore swab, is exposed to flowing wastewater over a set period to capture microorganisms. Composite sampling (n = 9) was also employed, involving autosamplers that collect wastewater at regular intervals over an extended period. Additionally, some studies (n = 7) employed a combination of grab and trap sampling, while five studies did not specify the type of sample collected. One study uniquely used sewage sludge as the sample type.
Sample handling after collection was reported in 61 studies (Table 1). Among these, 56 studies described the transportation of samples under cold chain conditions to preserve microbial integrity. Additionally, 52 studies provided details on the time elapsed between sample collection and laboratory processing. Of these, 49 studies processed samples within 24 h, while three studies reported delays exceeding 24 h.
Sample processing, defined as the steps taken before initiating the microbiological procedure, was also mentioned by various methods (n = 74) (Table 1). The methods reported processing using filtration (n = 32) to trap microorganisms or remove large debris, centrifugation (n = 8) to pellet the microorganisms, dilution or serial dilution (n = 5) to reduce inhibitory components, processing of the Moore swab to extract its contents (n = 3), using specially designed magnetic particles that bind with bacteria (n = 1), or a combination of methods (n = 23) to remove debris and inhibitors and trap microorganisms.
The laboratory testing could be further resolved in different protocol steps, including bacterial culture, bacterial enumeration, phenotypic characterization of observable bacterial traits (biochemical identification, serotyping, antimicrobial susceptibility, and other methods), genotypic characterization based on bacterial nucleic acid (PCR and other methods), and sequencing (Figure 3A). Two studies used a novel bacteriophage-based method that employed detecting Salmonella-specific bacteriophages as an indicator for Salmonella. It was noted that the use of different protocol steps is influenced by the central question or hypothesis that these studies aimed to address.
(A) Count of steps used by each method (n = 102); (B) Upset plot for culture steps used in the methods, showing intersection size, that is, the count of methods using culture steps either alone or in combination.
Bacterial culture was attempted by 70% of the methods (n = 71), as shown in Figure 3A. This step typically involved enrichment in non-selective or selective media, followed by selective culture or isolation of Salmonella spp., using standard culture media. However, not all studies included a culture step as part of their methodology, some relied solely on molecular or alternative detection techniques. Figure 3B illustrates the combination of culture steps used in different methods (n = 67), excluding four methods that reported using standard methods (ISO 6579, ISO 19250, FDA Bioanalytical manual protocol for Salmonella, and APHA standard method) and one method that exclusively used bacteriophage-specific methods. Among those that performed culture, the most common protocol involved enrichment, selective enrichment, and selective culture (n = 32), followed by selective enrichment and culture (n = 13). Additionally, four methods reported bacterial enumeration using serial dilutions and the most probable number (MPN) method to estimate bacterial load.
Figure 3A also illustrates how different phenotypic characterization assays were primarily used to characterize isolated bacteria based on observable traits, such as growth in specific media, serotype, or antimicrobial susceptibility. The phenotypic methods included biotyping with biochemical media (n = 48), serotyping (n = 26), antimicrobial susceptibility testing (n = 47), and other phenotypic techniques mainly involving phage typing (n = 3), Matrix-Assisted Laser Desorption/Ionization Time-of-Flight, MALDI-TOF (n = 2), and both MALDI-TOF and phage typing (n = 1). Biotyping utilized standard biochemical identification techniques, either manual (n = 32), automated (n = 12), or a combination of both (n = 3). One study did not specify the biochemicals used. For serotyping, most studies employed the Kauffman-White serotyping scheme to characterize Salmonella isolates (n = 24). Two studies reported the use of a PCR-based serotyping scheme. Regarding antimicrobial susceptibility, the majority of studies used the disc diffusion assay (n = 36), followed by broth microdilution (n = 4), automated systems (n = 3), a combination of disc diffusion and automated systems (n = 2), and a combination of disc diffusion and broth microdilution (n = 1). One of the studies also utilized resistance transfer testing to understand the mechanism of AMR gene transfer to a susceptible host.
Genotypic characterization of bacterial nucleic acid by molecular assays primarily included variants of PCR (n = 65), as well as other molecular assays, such as pulsed field gel electrophoresis (n = 10), plasmid analysis (n = 2), DNA fingerprinting (n = 1), fluorescence in situ hybridization (n = 1), and sequencing (n = 28). The methods reported included qPCR (n = 22), PCR (n = 18), molecular assays for antimicrobial resistance genes (ARG) (n = 12), and a combination of qPCR and ARG PCR (n = 2). Other than that, one study each reported using multiplex PCR, RT-PCR BioFire FilmArray® panel, crystal digital PCR, high-throughput qPCR, PCR with virulence marker PCR, culture PCR with high-throughput qPCR, 16S RNA PCR with PCR, qPCR with high-throughput qPCR, PCR with qPCR, and qPCR with denaturing gradient gel electrophoresis. Most used genomic method was whole genome sequencing, WGS (n = 11), followed by untargeted or shotgun metagenomics (n = 6), and 16S rDNA targeted sequencing (n = 2). One study each reported targeted sequencing, 16S rRNA sequencing, ARG genes sequencing, untargeted metagenomics, long and short read sequencing using Nanopore and Illumina sequencing, 16S rRNA sequencing with WGS, 16S RNA sequencing with metagenomics, and 16S rDNA sequencing with biomarker sequencing.
Identification of laboratory testing pathways
3.4
To further understand the commonalities in the reported methodologies of WES for Salmonella spp. In the selected papers an exploratory analysis was performed to identify major testing pathways from sample collection to testing. In this approach, 87 different methods from 79 studies were employed (Supplementary Table 10). Fifteen methods from 14 studies were excluded from this analysis because no sample was specified, standard procedures were not discussed in detail, bacteriophage surveillance was not included, or methods used for stored isolates were not specified.
Based on the clustering, a total of six pathways were identified (Figure 4A; Supplementary Figure 3). The details of each method mapped to a pathway are provided in Supplementary Table 10. Eight methods shared the pathway P1, mainly involving a culture step followed by identification using molecular methods, with a slight association with a processing step. Pathway P2 was the most used pathway, with 43 methods. P2 has a strong association with the processing step, culture, biotyping of isolates, and antimicrobial susceptibility testing. The pathway was mildly associated with molecular assays, with low association to other steps. P3 was shared by seven methods and was associated with culture, biotyping, serotyping, and AST. Seven studies shared pathway P4, which was strongly associated with a processing step and the use of molecular assays for characterization. P4 also had a moderate association with culture and genomics methods. Pathway P5 was the second most used pathway, with 15 methods, primarily involving molecular characterization after a processing step. Lastly, pathway P6 was shared by seven methods and included a processing step and testing using a genomics method. Figure 4B illustrates the fractional distribution of methodological pathways across sample types in LMIC and HIC settings. Both LMICs and HICs studies commonly employed pathways P2, P5, and P6 for grab samples. Pathway P2 was used at similar rates in both income groups, while P5 was more prevalent in LMICs and P6 in HICs settings. In LMICs, P5 was also applied to composite and trap samples. For trap sampling, pathway P1 was used equally by both LMIC and HIC studies. Pathway P3 was preferred in HICs, whereas LMICs applied it to both trap and composite samples. Pathway P4 was predominantly associated with composite sampling across studies. Stability analysis of pathways confirmed that the identified pathways represent methodological choices by the authors. The exclusion of unclassified methods showed significantly lower reporting completeness score compared to the clustered pathways (p = 0.02) confirming filtering out these methods from the pathway clustering due to insufficient reported details. For the clustered methods, reporting scores were uniform with no difference in completeness scores detected by ANOVA (p = 0.17). To control for potential pseudoreplication, we performed a sensitivity analysis by aggregating completeness score at study level (n = 94), which yielded consistent result (p = 0.17). A stricter sensitivity analysis excluding removing papers with multiple methods was ruled out as it resulted in dropping of Pathway P1 sample size below the threshold required for reliable variance estimation (n = 3). Bootstrap analysis further demonstrated high internal stability, with CV remaining below 0.08 for all pathways (Supplementary Table 12). These results, although exploratory, indicate that the methodologies utilized for Salmonella WES can be generalized at the level of protocol steps, thereby providing an opportunity to identify critical steps and define quality control protocols for Salmonella WES.
(A) Clustering of methods used to identify the pathways involved in wastewater sample processing and testing, based on the findings from 79 extracted studies that included 87 different methods. For each pathway cluster, the number of methods is indicated; (B) Fraction plot showing the use of pathways for testing by sample type and country category (LMIC group includes LIC, LMIC, and UMIC).
Domain mapping for selected studies
3.5
To further understand the thematic domains of the studies to provide context for planning of Salmonella WES, an exploratory analysis was done to assign thematic domainsto all 94 unique studies included in the systematic review (Figure 5A). Two authors identified the thematic domains examining titles, keywords, and abstracts; where structured abstracts were unavailable, the first page was reviewed. The details of the identified study domains are provided in Supplementary Table 11. Each study could be mapped to one or more domains. A total of eight domains were identified: Domain A - Outbreak detection and investigation (n = 33): studies using WES to detect, investigate, or retrospectively link outbreaks. Domain B - Disease prevalence (n = 79): studies using WES to supplement clinical or sentinel surveillance for estimating disease burden. Domain C - AMR prevalence (n = 73): studies assessing the spread of AMR organisms or genes in the community. Domain D - Mechanisms of AMR (n = 31): studies exploring the physiological, molecular, or genetic mechanisms of transmission of resistance via wastewater. Domain E - Wastewater usage monitoring (n = 19): studies evaluating microbial diversity in reclaimed or contaminated water used for agriculture or irrigation. Domain F - Environmental health (n = 62): studies investigating links between wastewater microbial diversity and anthropogenic or environmental factors. Domain G - One Health (n = 23): studies addressing human-animal-environment interactions, including cross-species transmission and intersectoral AMR evidence. Domain H - Method validation (n = 20): studies focused on validating WES methods for sample collection or testing, including assessment of assay sensitivity and standardization of protocols.
(A) Study domain co-occurrence map (n = 94); (B) Network of co-occurrence map of thematic domains (n = 94).
The co-occurrence matrix of the identified domains was mapped onto a network to understand relatedness between the identified thematic domains from the dataset of 94 included studies (Figure 5B). The network analysis revealed that domains B (Disease Prevalence), C (AMR Prevalence), and F (Environmental Health) are closely grouped. These may be considered as core domains that guide key questions in the field, such as how WES can detect the presence or absence of Salmonella spp., assess AMR, and understand the role of environmental factors in pathogen establishment. Domains D (Mechanisms of AMR) and E (Wastewater Usage Monitoring) may act as bridge or support domains, connecting strongly with the core domains and providing important context, such as understanding resistance mechanisms or identifying sources of contamination through wastewater reuse. Domains G (One Health) and H (Method Validation) may be linked to multiple domains but in smaller numbers, indicating their cross-cutting nature. These domains may contribute to broader perspectives, such as intersectoral collaboration or methodological rigor. Finally, domain A (Outbreak Detection and Investigation) is relatively isolated, with fewer connections to other domains. This suggests that studies in this domain may often require specialized approaches and may not overlap extensively with broader surveillance objectives. The identified domains can be used to provide important context to the Salmonella WES, especially in cases where public health context is not sufficiently provided by Salmonella Typhi.
Completeness of reported methods
3.6
To evaluate the completeness of reporting on the methods, a structured template (Supplementary Table T9) was used to determine whether the studies thoroughly documented the wastewater methodology. The extracted data is provided in the Supplementary material 2. Briefly, the manuscripts were reviewed by two authors and graded as per the rubric outlined in Supplementary Table T9. All the disagreements were resolved by consensus between all authors. The methodological quality assessment criteria also classified the studies based on the study methodology (Figure 1), however, this metric only classifies the studies based on completeness of the reporting of the methodology and not on the quality of the method. Supplementary Figure 4 illustrates how studies (n = 94) reported the methodology across various assessment criteria. Nearly all studies clearly provided information on sampling site details (n = 79), sample processing (n = 86), and testing methods (n = 79). Additionally, the choice of site, based on the hypothesis or central question posed by the study (n = 60), sample collection details (n = 60), details of testing procedures, including reagents (n = 53), and sample transport conditions with transient times (n = 46), were reported inconsistently. The reporting of quality control procedures used for laboratory methods remains the only criterion that is not frequently reported (n = 13). The completeness assessment of reported methodology highlighted that due to lack of widely accepted standardized methodology, most of the included studies documented aspects such as site selection, sample handling, transport, and testing with reasonable consistency, while under-reporting the quality control measures needed to ensure reproducibility of results.
Recommendations on protocols and public health question framing
3.7
In view of the lack of standardized methodology or guidance on Salmonella wastewater surveillance, we have attempted to outline recommendations that could support countries in initiating WES for Salmonella (Figure 6). Figure 6 provides detailed technical guidance on various steps that should be considered while implementing WES for Salmonella. It is also essential for the reproducibility of the selected method that all information for the chosen steps is included, either within the manuscript, Supplementary material, or as an online-published protocol during publication.
Recommendations for details to be included while reporting the WES methodologies for Salmonella.
The knowledge gained in this systematic review is synthesized into a heuristic framework for resource planning to assist all stakeholders, including public health program managers, laboratory directors, donors, and policy makers, in framing the right questions (Table 2). The developed framework provides essential context for defining the public health goals of establishing WES for Salmonella, including its scope and domain. The core domains may help shape the main questions, while the supporting domains add additional objectives may strengthen the impact of the central questions. Cross-cutting domains can potentially provide context when multi-sectoral engagement is necessary. Niche domains may be used in conjunction with other domains to frame questions but will often require specific objectives that may not be relevant to different sectors. The framed question then guides the selection of suitable output data and sites to meet the goals, as well as justifying the testing pathway based on the available resources and infrastructure. The chosen testing pathway can be made reproducible by following the recommendations in Figure 6. Therefore, the proposed framework serves as a foundational evidence-based instrument for resource planning while simultaneously facilitating consensus building for among various stakeholders for initiating Salmonella WES.
Discussion
4
Summary of key findings
4.1
This systematic review identified substantial heterogeneity in laboratory methodologies used for WES of Salmonella spp. across 94 studies spanning 36 countries. A total of 102 distinct methodological approaches were documented, reflecting wide variation in sampling strategies, sample types, and laboratory testing protocols. Grab sampling was the most common method, although trap and composite sampling were also used, often without consistent reporting on sample handling or transport conditions. Culture-based methods were frequently employed, yet many studies relied solely on molecular or genomic techniques, underscoring the lack of standardized testing pathways. Six distinct methodological pathways were identified, each reflecting different combinations of protocol steps and resource contexts.
Domain mapping showed that studies often addressed multiple public health objectives, with disease prevalence, AMR, and environmental health emerging as core domains. Outbreak detection and method validation were underrepresented, suggesting a need for targeted investment in these areas. The interdisciplinary nature of Salmonella WES, encompassing One Health, environmental monitoring, and epidemiology, reflects its evolving role in integrated public health action (127, 128). The identified domains emphasize the importance of Salmonella WES in providing community-level signals for disease prevalence, information on pathogen importation, characterization of AMR emergence and spread, and assessment of environmental transmission risk, which directly align with the WHO-defined potential use cases for routine WES (27).
Context for interpretation of results
4.2
The findings must be interpreted in the context of infrastructural and epidemiological realities in LMICs, particularly in the WHO South-East Asia Region. Rapid urbanization has outpaced the development of sanitation infrastructure, resulting in fragmented wastewater systems that complicate the recovery and surveillance of pathogens (129, 130). Wastewater reuse for agriculture and urban needs has increased exposure to waterborne diseases, including typhoid fever (131–137).
Historically, Salmonella Typhi was among the first pathogens monitored in sewage to identify asymptomatic carriers (138), and this approach has been used to locate transmission hotspots (139). Despite its early promise, WES for typhoid has remained underutilized, with poliovirus being the only pathogen for which environmental surveillance is widely institutionalized (2–4, 140). The infrastructure established for wastewater sample collection and molecular testing for polio ES could be leveraged to initiate Salmonella WES. The COVID-19 pandemic catalyzed renewed interest in wastewater-based surveillance, demonstrating its utility for early detection and public health decision-making (140, 141).
The post-pandemic surge in publications reflects this shift, with studies emerging from both LMICs and HICs. Unconventional sampling sources; such as aircraft, refugee ships, and border entry points, have expanded the scope of surveillance (142–145). However, the dominance of grab sampling, limited use of Moore swabs (146), and inconsistent reporting of sample handling suggest that feasibility often outweighs methodological rigor.
Implications of results in the WES domain
4.3
The lack of standardization has direct implications for reproducibility and comparability. Our quality assessment revealed that while most studies reported basic methodological components, critical details such as quality control procedures, reagent specifications, and validation criteria were frequently omitted (139, 147–150). This gap limits the utility of published protocols for replication or scale-up in other settings.
To address this, we categorized the methods into six distinct pathways based on protocol steps and resource contexts. This classification offers a pragmatic framework for selecting appropriate methodologies aligned with laboratory capacity and surveillance goals. Importantly, the heuristic framework for resource planning developed from this synthesis (Table 2) enables stakeholders to align methodological choices with public health objectives, whether for outbreak detection, disease burden estimation, or AMR monitoring.
The interdisciplinary nature of Salmonella WES calls for integrated policy frameworks. Surveillance programs should be embedded within broader public health strategies that facilitate cross-sectoral collaboration among human, animal, and environmental health agencies (127, 128). This aligns with the Quadripartite One Health Joint Plan of Action and supports the development of multisectoral early warning systems for emerging pathogens (151).
Limitations of the study
4.4
Despite efforts to optimize the search strategy, the review may have missed relevant studies published in non-indexed journals or in languages other than English. The reliance on published literature means that methodological details were often incomplete or inconsistently reported, particularly regarding quality control procedures, reagent specifications, and validation criteria. This limited the ability to fully assess reproducibility and operational feasibility. The review also did not include grey literature, internal reports, or unpublished protocols, which may contain valuable insights into real-world implementation challenges. While the heuristic framework for resource planning, thematic domain identification, and pathway classification are grounded in extracted data, they have not yet been validated through field testing or stakeholder consultation, which we aim to do as a next step. The review focused exclusively on wastewater or contaminated surface waters sources and excluded other environmental matrices such as surface waters, which may also be relevant for Salmonella surveillance in certain contexts such as one health.
Recommendations and future directions
4.5
To advance WES for typhoid control, we recommend the development and adoption of standardized protocols that are both scientifically robust and operationally feasible across varied infrastructure settings. These protocols should include clear specifications for sample collection (e.g., volume, timing, and type), transport conditions, and laboratory testing workflows along with validation criteria for result interpretation. The consistent use of Moore swabs, which have demonstrated superior sensitivity in flowing wastewater environments, should be encouraged in typhoid-endemic regions (146). Furthermore, quality control procedures must be embedded throughout the surveillance process, with explicit documentation of reagents, test conditions, and validation criteria to ensure methodological transparency and reliability (139, 147–150).
Surveillance systems must be tailored to local resource contexts. The six methodological pathways identified in this review provide a flexible framework for laboratories with differing capacities. These pathways can be matched to surveillance objectives using the heuristic framework developed herein, which links public health goals to appropriate testing strategies and site selection models. Such alignment is particularly critical in LMICs, where infrastructure constraints necessitate pragmatic and cost-effective approaches.
Strategically, Salmonella WES should be integrated into national disease control programs and linked to typhoid conjugate vaccine (TCV) deployment. Environmental data can complement clinical surveillance and inform the timing and targeting of vaccination campaigns, especially in settings where blood culture-based diagnostics are limited (27, 152–154).
Future studies should evaluate the sensitivity and cost-effectiveness of composite sampling approaches, particularly in decentralized and low-flow wastewater systems common in LMICs. Comparative assessments of the six identified methodological pathways under field conditions are needed to determine which combinations of protocol steps yield the most reliable results across diverse environmental contexts.
Additionally, research should explore the integration of WES data with clinical and AMR surveillance systems to enhance early warning capabilities and inform vaccine deployment strategies. The operational feasibility of implementing surveillance in unconventional settings; such as border crossings, refugee camps, and transportation hubs, also requires further study, especially in light of emerging global health threats.
Finally, interdisciplinary implementation models that align with One Health frameworks should be piloted and evaluated. These models must address not only technical performance but also governance, stakeholder engagement, and sustainability in resource-limited settings.
Conclusion
5
WES for Salmonella spp., particularly Salmonella Typhi, presents a promising yet underutilized tool for public health action in resource-limited settings. Standardized protocol and its harmonized implementation are thus a key success factor for Salmonella WES, however, it necessitates scientific and operational research to achieve the outcome. Leveraging existing systems for polio ES has potential to expedite the implementation of Salmonella WES. This review offers a foundation for methodological harmonization, strategic integration, and interdisciplinary collaboration. By aligning surveillance design with public health objectives and local capacities, stakeholders can advance robust, scalable systems that support typhoid control and broader One Health goals.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1GPEI. Global polio eradication initiative Geneva, Switzerland: World Health Organization (2025) Available online at: https://polioeradication.org/.
- 2World Health Organization. Guidelines for environmental surveillance of poliovirus circulation. Geneva, Switzerland: World Health Organization, Department of Vaccines and Biologicals (2003).
- 3GPEI. Guidelines on environmental surveillance for detection of poliovirus. Guideline. Geneva: Global Polio Eradication Initiative, World Health Organization, Initiative. GPE (2015).
- 4GPEI. Field guidance for the implementation of environmental surveillance for poliovirus. Geneva: Global Polio Eradication Initiative, World Health Organization (2023).
- 5Levy JI Andersen KG Knight R Karthikeyan S. Wastewater surveillance for public health. Science. (2023) 379:26–7. doi: 10.1126/science.ade 2503, 36603089 PMC 10065025 · doi ↗ · pubmed ↗
- 6Grassly NC Shaw AG Owusu M. Global wastewater surveillance for pathogens with pandemic potential: opportunities and challenges. Lancet Microbe. (2025) 6:100939. doi: 10.1016/j.lanmic.2024.07.002, 39222653 · doi ↗ · pubmed ↗
- 7Tiwari A Radu E Kreuzinger N Ahmed W Pitkänen T. Key considerations for pathogen surveillance in wastewater. Sci Total Environ. (2024) 945:173862. doi: 10.1016/j.scitotenv.2024.173862, 38876348 · doi ↗ · pubmed ↗
- 8Toro L de Valk H Zanetti L Huot C Tarantola A Fournet N . Pathogen prioritisation for wastewater surveillance ahead of the Paris 2024 Olympic and Paralympic games, France. Euro Surveill. (2024) 29:2400231. doi: 10.2807/1560-7917.ES.2024.29.28.2400231, 38994605 PMC 11241851 · doi ↗ · pubmed ↗
