{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,15]],"date-time":"2026-09-15T22:32:58Z","timestamp":1789511578635,"version":"build-2803163510"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1013749","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2025,12,10]],"date-time":"2025-12-10T00:00:00Z","timestamp":1765324800000}}],"reference-count":41,"publisher":"Public Library of Science (PLoS)","issue":"12","license":[{"start":{"date-parts":[[2025,12,3]],"date-time":"2025-12-03T00:00:00Z","timestamp":1764720000000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"ANRS-MIE","award":["ANRS 0160"],"award-info":[{"award-number":["ANRS 0160"]}]},{"DOI":"10.13039\/501100001665","name":"Agence Nationale de la Recherche","doi-asserted-by":"publisher","award":["ANR-24-MIEM-0004 under the France 2030 programme"],"award-info":[{"award-number":["ANR-24-MIEM-0004 under the France 2030 programme"]}],"id":[{"id":"10.13039\/501100001665","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>Repeated waves of emerging variants during the SARS-CoV-2 pandemics have highlighted the urge of collecting longitudinal genomic data and developing statistical methods based on time series analyses for detecting new threatening lineages and estimating their fitness early in time. Most models study the evolution of the prevalence of particular lineages over time and require a prior classification of sequences into lineages which is prone to induce delays and biases. More recently, several authors studied the evolution of the prevalence of mutations over time with alternative clustering approaches, avoiding specific lineage classification. Most existing methods are either non parametric or unsuited to pooled data characterizing, for instance, wastewater samples. The analysis of wastewater samples has recently been pointed out as a valuable complementary approach to clinical sample analysis, however the pooled nature of the data involves specific statistical challenges. In this context, we propose an alternative unsupervised method for clustering mutations according to their frequency trajectory over time and estimating group fitness from time series of pooled mutation prevalence data. Our model is a mixture of observed count data and latent group assignment and we use the expectation-maximization algorithm for model selection and parameter estimation. The application of our method to time series of SARS-CoV-2 sequencing data collected from wastewater treatment plants in France from October 2020 to April 2021 shows its ability to agnostically group mutations in a consistent way with lineages B.1.160, Alpha, B.1.177, Beta, and with selection coefficient estimates per group in coherence with the viral dynamics in France reported by Nextstrain. Moreover, our method detected the Alpha variant as threatening as early as supervised methods (which track specific mutations over time) with the noticeable difference that, since unsupervised, it does not require any prior information on the set of mutations.<\/jats:p>","DOI":"10.1371\/journal.pcbi.1013749","type":"journal-article","created":{"date-parts":[[2025,12,3]],"date-time":"2025-12-03T18:43:26Z","timestamp":1764787406000},"page":"e1013749","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":3,"title":["Unsupervised detection and fitness estimation of emerging SARS-CoV-2 variants: Application to wastewater samples (ANRS0160)"],"prefix":"10.1371","volume":"21","author":[{"ORCID":"https:\/\/orcid.org\/0009-0006-9301-5552","authenticated-orcid":true,"given":"Alexandra","family":"Lefebvre","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Vincent","family":"Mar\u00e9chal","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Arnaud","family":"Gloaguen","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"name":"The Ob\u00e9pine Consortium","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Amaury","family":"Lambert","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yvon","family":"Maday","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"340","published-online":{"date-parts":[[2025,12,3]]},"reference":[{"key":"pcbi.1013749.ref001","doi-asserted-by":"crossref","DOI":"10.1016\/j.lanepe.2021.100202","article-title":"Monitoring SARS-CoV-2 variants alterations in Nice neighborhoods by wastewater nanopore sequencing","volume":"10","author":"G Rios","year":"2021","journal-title":"The Lancet Regional Health\u2013Europe."},{"issue":"50","key":"pcbi.1013749.ref002","doi-asserted-by":"crossref","first-page":"2000776","DOI":"10.2807\/1560-7917.ES.2020.25.50.2000776","article-title":"Evaluation of lockdown effect on SARS-CoV-2 dynamics through viral genome quantification in waste water, Greater Paris, France, 5 March to 23 April 2020","volume":"25","author":"S Wurtzer","year":"2020","journal-title":"Euro Surveill."},{"issue":"10","key":"pcbi.1013749.ref003","doi-asserted-by":"crossref","first-page":"1144","DOI":"10.3390\/v12101144","article-title":"Tracking SARS-CoV-2 in sewage: evidence of changes in virus variant predominance during COVID-19 pandemic","volume":"12","author":"J Martin","year":"2020","journal-title":"Viruses."},{"key":"pcbi.1013749.ref004","doi-asserted-by":"crossref","first-page":"110929","DOI":"10.1016\/j.envres.2021.110929","article-title":"A review of the presence of SARS-CoV-2 RNA in wastewater and airborne particulates and its use for virus spreading surveillance","volume":"196","author":"U Anand","year":"2021","journal-title":"Environ Res."},{"key":"pcbi.1013749.ref005","doi-asserted-by":"crossref","first-page":"150121","DOI":"10.1016\/j.scitotenv.2021.150121","article-title":"SARS-CoV-2 RNA concentrations in wastewater foreshadow dynamics and clinical presentation of new COVID-19 cases","volume":"805","author":"F Wu","year":"2022","journal-title":"Sci Total Environ."},{"issue":"8","key":"pcbi.1013749.ref006","doi-asserted-by":"crossref","first-page":"1151","DOI":"10.1038\/s41564-022-01185-x","article-title":"Early detection and surveillance of SARS-CoV-2 genomic variants in wastewater using COJAC","volume":"7","author":"K Jahn","year":"2022","journal-title":"Nat Microbiol."},{"issue":"1","key":"pcbi.1013749.ref007","doi-asserted-by":"crossref","first-page":"18721","DOI":"10.1038\/s41598-022-22060-z","article-title":"SARS-CoV-2 quasi-species analysis from patients with persistent nasopharyngeal shedding","volume":"12","author":"P Dudouet","year":"2022","journal-title":"Sci Rep."},{"key":"pcbi.1013749.ref008","doi-asserted-by":"crossref","first-page":"138764","DOI":"10.1016\/j.scitotenv.2020.138764","article-title":"First confirmed detection of SARS-CoV-2 in untreated wastewater in Australia: a proof of concept for the wastewater surveillance of COVID-19 in the community","volume":"728","author":"W Ahmed","year":"2020","journal-title":"Sci Total Environ."},{"issue":"7","key":"pcbi.1013749.ref009","doi-asserted-by":"crossref","first-page":"511","DOI":"10.1021\/acs.estlett.0c00357","article-title":"Presence of SARS-coronavirus-2 RNA in sewage and correlation with reported COVID-19 prevalence in the early stage of the epidemic in The Netherlands","volume":"7","author":"G Medema","year":"2020","journal-title":"Environ Sci Technol Lett."},{"key":"pcbi.1013749.ref010","doi-asserted-by":"crossref","first-page":"117104","DOI":"10.1016\/j.watres.2021.117104","article-title":"Rapid screening for SARS-CoV-2 variants of concern in clinical and environmental samples using nested RT-PCR assays targeting key mutations of the spike protein","volume":"197","author":"G La Rosa","year":"2021","journal-title":"Water Res."},{"issue":"10","key":"pcbi.1013749.ref011","doi-asserted-by":"crossref","first-page":"2100806","DOI":"10.2807\/1560-7917.ES.2022.27.10.2100806","article-title":"Inferring transmission fitness advantage of SARS-CoV-2 variants of concern from wastewater samples using digital PCR, Switzerland, December 2020 through March 2021","volume":"27","author":"L Caduff","year":"2022","journal-title":"Euro Surveill."},{"key":"pcbi.1013749.ref012","doi-asserted-by":"crossref","first-page":"106998","DOI":"10.1016\/j.envint.2021.106998","article-title":"A nationwide indicator to smooth and normalize heterogeneous SARS-CoV-2 RNA data in wastewater","volume":"158","author":"N Cluzel","year":"2022","journal-title":"Environ Int."},{"key":"pcbi.1013749.ref013","doi-asserted-by":"crossref","unstructured":"Beerenwinkel N, Dreifuss D, Huisman J, Rusch J, Caduff L, Ganesanandamoorthy P, et al. Estimated transmission dynamics of SARS-CoV-2 variants from wastewater are robust to differential shedding. Research Square. 2023.","DOI":"10.21203\/rs.3.rs-3719342\/v1"},{"issue":"5","key":"pcbi.1013749.ref014","doi-asserted-by":"crossref","first-page":"1405","DOI":"10.3201\/eid2705.204410","article-title":"Monitoring SARS-CoV-2 circulation and diversity through community wastewater sequencing, the Netherlands and Belgium","volume":"27","author":"R Izquierdo-Lara","year":"2021","journal-title":"Emerg Infect Dis."},{"issue":"1","key":"pcbi.1013749.ref015","doi-asserted-by":"crossref","DOI":"10.1128\/mBio.02703-20","article-title":"Genome sequencing of sewage detects regionally prevalent SARS-CoV-2 variants","volume":"12","author":"A Crits-Christoph","year":"2021","journal-title":"mBio."},{"key":"pcbi.1013749.ref016","doi-asserted-by":"crossref","first-page":"117710","DOI":"10.1016\/j.watres.2021.117710","article-title":"High-throughput sequencing of SARS-CoV-2 in wastewater provides insights into circulating variants","volume":"205","author":"RS Fontenele","year":"2021","journal-title":"Water Res."},{"issue":"7","key":"pcbi.1013749.ref017","doi-asserted-by":"crossref","first-page":"1809","DOI":"10.1093\/bioinformatics\/btac047","article-title":"A mixture model for determining SARS-Cov-2 variant composition in pooled samples","volume":"38","author":"R Valieris","year":"2022","journal-title":"Bioinformatics."},{"issue":"8","key":"pcbi.1013749.ref018","doi-asserted-by":"crossref","first-page":"1647","DOI":"10.3390\/v13081647","article-title":"Monitoring SARS-CoV-2 populations in wastewater by amplicon sequencing and using the novel program SAM refiner","volume":"13","author":"DA Gregory","year":"2021","journal-title":"Viruses."},{"issue":"1","key":"pcbi.1013749.ref019","doi-asserted-by":"crossref","first-page":"635","DOI":"10.1038\/s41467-022-28246-3","article-title":"Tracking cryptic SARS-CoV-2 lineages detected in NYC wastewater","volume":"13","author":"DS Smyth","year":"2022","journal-title":"Nat Commun."},{"issue":"6538","key":"pcbi.1013749.ref020","doi-asserted-by":"crossref","DOI":"10.1126\/science.abg3055","article-title":"Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England","volume":"372","author":"NG Davies","year":"2021","journal-title":"Science."},{"issue":"7858","key":"pcbi.1013749.ref021","doi-asserted-by":"crossref","first-page":"266","DOI":"10.1038\/s41586-021-03470-x","article-title":"Assessing transmissibility of SARS-CoV-2 lineage B.1.1.7 in England","volume":"593","author":"E Volz","year":"2021","journal-title":"Nat."},{"issue":"5","key":"pcbi.1013749.ref022","doi-asserted-by":"crossref","first-page":"2100429","DOI":"10.2807\/1560-7917.ES.2022.27.5.2100429","article-title":"Co-circulation of SARS-CoV-2 Alpha, Gamma variants in Italy, February and March 2021","volume":"27","author":"P Stefanelli","year":"2022","journal-title":"Euro Surveill."},{"key":"pcbi.1013749.ref023","doi-asserted-by":"crossref","unstructured":"Althaus CL, Baggio S, Reichmuth ML, Hodcroft EB, Riou J, Neher RA, et al. A tale of two variants: spread of SARS-CoV-2 variants Alpha in Geneva, Switzerland, and Beta in South Africa. MedRxiv. 2021:2021-06.","DOI":"10.1101\/2021.06.10.21258468"},{"key":"pcbi.1013749.ref024","doi-asserted-by":"crossref","unstructured":"Nishiura H, Ito K, Anzai A, Kobayashi T, Piantham C, Rodr\u00edguez-Morales AJ. Relative reproduction number of SARS-CoV-2 Omicron (B. 1.1. 529) compared with Delta variant in South Africa, 2021.","DOI":"10.3390\/jcm11010030"},{"issue":"1","key":"pcbi.1013749.ref025","doi-asserted-by":"crossref","first-page":"7239","DOI":"10.1038\/s41467-021-27369-3","article-title":"Estimating the strength of selection for new SARS-CoV-2 variants","volume":"12","author":"van Dorp CH","year":"2021","journal-title":"Nat Commun."},{"issue":"2","key":"pcbi.1013749.ref026","doi-asserted-by":"crossref","DOI":"10.1093\/ve\/veac089","article-title":"Global estimates of the fitness advantage of SARS-CoV-2 variant Omicron","volume":"8","author":"C van Dorp","year":"2022","journal-title":"Virus Evol."},{"issue":"6599","key":"pcbi.1013749.ref027","doi-asserted-by":"crossref","first-page":"1327","DOI":"10.1126\/science.abm1208","article-title":"Analysis of 6.4 million SARS-CoV-2 genomes identifies mutations associated with fitness","volume":"376","author":"F Obermeyer","year":"2022","journal-title":"Science."},{"issue":"11","key":"pcbi.1013749.ref028","doi-asserted-by":"crossref","first-page":"1403","DOI":"10.1038\/s41564-020-0770-5","article-title":"A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology","volume":"5","author":"A Rambaut","year":"2020","journal-title":"Nat Microbiol."},{"issue":"1","key":"pcbi.1013749.ref029","doi-asserted-by":"crossref","first-page":"21068","DOI":"10.1038\/s41598-021-00496-z","article-title":"Data-driven analysis of amino acid change dynamics timely reveals SARS-CoV-2 variant emergence","volume":"11","author":"A Bernasconi","year":"2021","journal-title":"Sci Rep."},{"key":"pcbi.1013749.ref030","doi-asserted-by":"crossref","first-page":"859241","DOI":"10.3389\/fmicb.2022.859241","article-title":"A new way to trace SARS-CoV-2 variants through weighted network analysis of frequency trajectories of mutations","volume":"13","author":"Q Huang","year":"2022","journal-title":"Front Microbiol."},{"issue":"1","key":"pcbi.1013749.ref031","doi-asserted-by":"crossref","first-page":"9275","DOI":"10.1038\/s41598-022-12442-8","article-title":"Variant-driven early warning via unsupervised machine learning analysis of spike protein mutations for COVID-19","volume":"12","author":"A de Hoffer","year":"2022","journal-title":"Sci Rep."},{"issue":"1","key":"pcbi.1013749.ref032","doi-asserted-by":"crossref","first-page":"3105","DOI":"10.1038\/s41467-023-38809-7","article-title":"Epidemiological inference for emerging viruses using segregating sites","volume":"14","author":"Y Park","year":"2023","journal-title":"Nat Commun."},{"issue":"25","key":"pcbi.1013749.ref033","doi-asserted-by":"crossref","DOI":"10.1073\/pnas.2314262121","article-title":"Estimation of SARS-CoV-2 fitness gains from genomic surveillance data without prior lineage classification","volume":"121","author":"T Donker","year":"2024","journal-title":"Proc Natl Acad Sci U S A."},{"issue":"1","key":"pcbi.1013749.ref034","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1111\/j.2517-6161.1977.tb01600.x","article-title":"Maximum likelihood from incomplete data via the EM algorithm","volume":"39","author":"AP Dempster","year":"1977","journal-title":"Journal of the Royal Statistical Society Series B: Statistical Methodology."},{"issue":"7","key":"pcbi.1013749.ref035","doi-asserted-by":"crossref","first-page":"719","DOI":"10.1109\/34.865189","article-title":"Assessing a mixture model for clustering with the integrated completed likelihood","volume":"22","author":"C Biernacki","year":"2000","journal-title":"IEEE Trans Pattern Anal Machine Intell."},{"key":"pcbi.1013749.ref036","doi-asserted-by":"crossref","first-page":"889811","DOI":"10.3389\/fmicb.2022.889811","article-title":"SARS-CoV-2 whole-genome sequencing using oxford nanopore technology for variant monitoring in wastewaters","volume":"13","author":"L Barb\u00e9","year":"2022","journal-title":"Front Microbiol."},{"issue":"2","key":"pcbi.1013749.ref037","doi-asserted-by":"crossref","first-page":"599","DOI":"10.1534\/genetics.112.140939","article-title":"Estimating allele age and selection coefficient from time-serial data","volume":"192","author":"A-S Malaspinas","year":"2012","journal-title":"Genetics."},{"issue":"2","key":"pcbi.1013749.ref038","doi-asserted-by":"crossref","first-page":"831","DOI":"10.1534\/genetics.115.184598","article-title":"An approximate Markov model for the wright-fisher diffusion and its application to time series data","volume":"203","author":"A Ferrer-Admetlla","year":"2016","journal-title":"Genetics."},{"issue":"12","key":"pcbi.1013749.ref039","doi-asserted-by":"crossref","first-page":"4073","DOI":"10.1534\/g3.119.400778","article-title":"Inference of selection from genetic time series using various parametric approximations to the wright-fisher model","volume":"9","author":"C Paris","year":"2019","journal-title":"G3 (Bethesda)."},{"issue":"1","key":"pcbi.1013749.ref040","doi-asserted-by":"crossref","first-page":"551","DOI":"10.1186\/s12859-022-05100-3","article-title":"VirPool: model-based estimation of SARS-CoV-2 variant proportions in wastewater samples","volume":"23","year":"2022","journal-title":"BMC bioinformatics."},{"issue":"1","key":"pcbi.1013749.ref041","doi-asserted-by":"crossref","first-page":"8","DOI":"10.1186\/s13059-018-1618-7","article-title":"An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar","volume":"20","author":"ND Grubaugh","year":"2019","journal-title":"Genome Biol."}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1013749","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2025,12,10]],"date-time":"2025-12-10T00:00:00Z","timestamp":1765324800000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1013749","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,12,10]],"date-time":"2025-12-10T18:44:46Z","timestamp":1765392286000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1013749"}},"subtitle":[],"editor":[{"given":"Andrey","family":"Rzhetsky","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2025,12,3]]},"references-count":41,"journal-issue":{"issue":"12","published-online":{"date-parts":[[2025,12,3]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1013749","relation":{},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2025,12,3]]}}}