{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,12,31]],"date-time":"2025-12-31T11:43:44Z","timestamp":1767181424709,"version":"build-2238731810"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1012391","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2025,7,18]],"date-time":"2025-07-18T00:00:00Z","timestamp":1752796800000}}],"reference-count":57,"publisher":"Public Library of Science (PLoS)","issue":"4","license":[{"start":{"date-parts":[[2025,4,24]],"date-time":"2025-04-24T00:00:00Z","timestamp":1745452800000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>\n                    We have developed a robust method for efficiently tracking multiple co-occurring mutations in a sequence database. Evolution often hinges on the interaction of several mutations to produce significant phenotypic changes that lead to the proliferation of a variant. However, identifying numerous simultaneous mutations across a vast database of sequences poses a significant computational challenge. Our approach leverages a matrix factorization technique to automatically and efficiently pinpoint subsets of positions where co-mutations occur, appearing in a substantial number of sequences within the database. We validated our method using SARS-CoV-2 receptor-binding domains, comprising approximately seven hundred thousand sequences of the Spike protein, demonstrating superior performance compared to a reasonably exhaustive brute-force method. Furthermore, we explore the biological significance of the identified\n                    <jats:italic>co-mutational positions<\/jats:italic>\n                    (CMPs) and their potential impact on the virus\u2019s evolution and functionality, identifying key mutations in Delta and Omicron variants. This analysis underscores the significant role of identified CMPs in understanding the evolutionary trajectory. By tracking the \u201cbirth\" and \u201cdeath\" of CMPs, we can elucidate the persistence and impact of specific groups of mutations across different viral strains, providing valuable insights into the virus\u2019 adaptability and thus, possibly aiding vaccine design strategies.\n                  <\/jats:p>","DOI":"10.1371\/journal.pcbi.1012391","type":"journal-article","created":{"date-parts":[[2025,4,24]],"date-time":"2025-04-24T16:02:37Z","timestamp":1745510557000},"page":"e1012391","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":0,"title":["Efficient discovery of frequently co-occurring mutations in a sequence database with matrix factorization"],"prefix":"10.1371","volume":"21","author":[{"ORCID":"https:\/\/orcid.org\/0009-0005-9738-2985","authenticated-orcid":true,"given":"Michael Robert","family":"Kolar","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Valerie","family":"Kobzarenko","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Debasis","family":"Mitra","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"340","published-online":{"date-parts":[[2025,4,24]]},"reference":[{"issue":"6","key":"pcbi.1012391.ref001","doi-asserted-by":"crossref","first-page":"1339","DOI":"10.3390\/biomedicines10061339","article-title":"SARS-CoV-2: a master of immune evasion","volume":"10","author":"A Rubio-Casillas","year":"2022","journal-title":"Biomedicines"},{"issue":"9","key":"pcbi.1012391.ref002","doi-asserted-by":"crossref","first-page":"1141","DOI":"10.1038\/s41401-020-0485-4","article-title":"Structural and functional properties of SARS-CoV-2 spike protein: potential antivirus drug development for COVID-19","volume":"41","author":"Y Huang","year":"2020","journal-title":"Acta Pharmacol Sin"},{"issue":"22","key":"pcbi.1012391.ref003","doi-asserted-by":"crossref","first-page":"12114","DOI":"10.3390\/ijms222212114","article-title":"SARS-CoV-2 variants, RBD mutations, binding affinity, and antibody escape","volume":"22","author":"L Yang","year":"2021","journal-title":"Int J Mol Sci"},{"issue":"15","key":"pcbi.1012391.ref004","doi-asserted-by":"crossref","first-page":"8606","DOI":"10.3390\/ijms23158606","article-title":"COVID-19 and the immune response: a multi-phasic approach to the treatment of COVID-19","volume":"23","author":"T Sapir","year":"2022","journal-title":"Int J Mol Sci"},{"issue":"1","key":"pcbi.1012391.ref005","doi-asserted-by":"crossref","first-page":"164","DOI":"10.1186\/s13073-021-00985-w","article-title":"Comprehensive mapping of binding hot spots of SARS-CoV-2 RBD-specific neutralizing antibodies for tracking immune escape variants","volume":"13","author":"C Yi","year":"2021","journal-title":"Genome Med"},{"issue":"11","key":"pcbi.1012391.ref006","doi-asserted-by":"crossref","first-page":"1535","DOI":"10.1016\/S1473-3099(22)00642-9","article-title":"Further humoral immunity evasion of emerging SARS-CoV-2 BA.4 and BA.5 subvariants","volume":"22","author":"F Jian","year":"2022","journal-title":"Lancet Infect Dis"},{"issue":"4","key":"pcbi.1012391.ref007","doi-asserted-by":"crossref","first-page":"2027","DOI":"10.1534\/genetics.103.023226","article-title":"Estimating the frequency of events that cause multiple-nucleotide changes","volume":"167","author":"S Whelan","year":"2004","journal-title":"Genetics"},{"key":"pcbi.1012391.ref008","doi-asserted-by":"crossref","first-page":"815389","DOI":"10.3389\/fmed.2022.815389","article-title":"SARS-CoV-2 mutations and their impact on diagnostics, therapeutics and vaccines","volume":"9","author":"S Thakur","year":"2022","journal-title":"Front Med (Lausanne)"},{"issue":"1","key":"pcbi.1012391.ref009","doi-asserted-by":"crossref","first-page":"8540","DOI":"10.1038\/s41598-022-12479-9","article-title":"Assessment of mutations on RBD in the Spike protein of SARS-CoV-2 Alpha, Delta and Omicron variants","volume":"12","author":"CHS da Costa","year":"2022","journal-title":"Sci Rep"},{"key":"pcbi.1012391.ref010","doi-asserted-by":"crossref","first-page":"5125","DOI":"10.1016\/j.csbj.2023.10.030","article-title":"Collaborative effects of 2019-nCoV-Spike mutants on viral infectivity","volume":"21","author":"S Fang","year":"2023","journal-title":"Comput Struct Biotechnol J"},{"issue":"5","key":"pcbi.1012391.ref011","doi-asserted-by":"crossref","first-page":"647","DOI":"10.1038\/s41423-022-00853-6","article-title":"A SARS-CoV-2 antibody retains potent neutralization against Omicron by targeting conserved RBM residues","volume":"19","author":"C Yi","year":"2022","journal-title":"Cell Mol Immunol"},{"key":"pcbi.1012391.ref012","doi-asserted-by":"crossref","first-page":"4238","DOI":"10.1016\/j.csbj.2022.07.051","article-title":"Analysis of co-occurring and mutually exclusive amino acid changes and detection of convergent and divergent evolution events in SARS-CoV-2","volume":"20","author":"R Al Khalaf","year":"2022","journal-title":"Comput Struct Biotechnol J"},{"issue":"7","key":"pcbi.1012391.ref013","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1000029","article-title":"Nonnegative matrix factorization: an analytical and interpretive tool in computational biology","volume":"4","author":"K Devarajan","year":"2008","journal-title":"PLoS Comput Biol"},{"issue":"13","key":"pcbi.1012391.ref014","doi-asserted-by":"crossref","DOI":"10.1093\/bioinformatics\/btr206","article-title":"A novel computational framework for simultaneous integration of multiple types of genomic data to identify microRNA-gene regulatory modules","volume":"27","author":"S Zhang","year":"2011","journal-title":"Bioinformatics"},{"issue":"11","key":"pcbi.1012391.ref015","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1011491","article-title":"Identifying promoter sequence architectures via a chunking-based algorithm using non-negative matrix factorisation","volume":"19","author":"S Nikumbh","year":"2023","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1012391.ref016","volume-title":"National Library of Medicine (US)","author":"National Center for Biotechnology Information (NCBI)","year":"2022"},{"issue":"2","key":"pcbi.1012391.ref017","doi-asserted-by":"crossref","DOI":"10.1093\/ve\/veab064","article-title":"Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool","volume":"7","author":"\u00c1 O\u2019Toole","year":"2021","journal-title":"Virus Evol"},{"issue":"1","key":"pcbi.1012391.ref018","doi-asserted-by":"crossref","first-page":"121","DOI":"10.1186\/s12864-022-08358-2","article-title":"Pango lineage designation and assignment using SARS-CoV-2 spike gene nucleotide sequences","volume":"23","author":"\u00c1 O\u2019Toole","year":"2022","journal-title":"BMC Genomics"},{"issue":"17","key":"pcbi.1012391.ref019","doi-asserted-by":"crossref","first-page":"4380","DOI":"10.1016\/j.cell.2021.06.008","article-title":"Identification of novel bat coronaviruses sheds light on the evolutionary origins of SARS-CoV-2 and related viruses","volume":"184","author":"H Zhou","year":"2021","journal-title":"Cell"},{"issue":"1","key":"pcbi.1012391.ref020","doi-asserted-by":"crossref","first-page":"95","DOI":"10.1186\/s12859-022-04637-7","article-title":"3GOLD: optimized Levenshtein distance for clustering third-generation sequencing data","volume":"23","author":"R Logan","year":"2022","journal-title":"BMC Bioinformatics"},{"issue":"46","key":"pcbi.1012391.ref021","doi-asserted-by":"crossref","first-page":"28572","DOI":"10.1073\/pnas.2014567117","article-title":"Evolution of the genetic code; Evidence from serine codon use disparity in Escherichia coli","volume":"117","author":"M Inouye","year":"2020","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"9","key":"pcbi.1012391.ref022","doi-asserted-by":"crossref","first-page":"1445","DOI":"10.1101\/gr.170696.113","article-title":"Error-prone polymerase activity causes multinucleotide mutations in humans","volume":"24","author":"K Harris","year":"2014","journal-title":"Genome Res"},{"issue":"11","key":"pcbi.1012391.ref023","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pgen.1006315","article-title":"Multi-nucleotide de novo Mutations in Humans","volume":"12","author":"S Besenbacher","year":"2016","journal-title":"PLoS Genet"},{"issue":"3","key":"pcbi.1012391.ref024","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pone.0248337","article-title":"Extra base hits: Widespread empirical support for instantaneous multiple-nucleotide changes","volume":"16","author":"AG Lucaci","year":"2021","journal-title":"PLoS One"},{"issue":"1544","key":"pcbi.1012391.ref025","doi-asserted-by":"crossref","first-page":"1203","DOI":"10.1098\/rstb.2009.0305","article-title":"Forces that influence the evolution of codon bias","volume":"365","author":"PM Sharp","year":"2010","journal-title":"Philos Trans R Soc Lond B Biol Sci"},{"issue":"6755","key":"pcbi.1012391.ref026","doi-asserted-by":"crossref","first-page":"788","DOI":"10.1038\/44565","article-title":"Learning the parts of objects by non-negative matrix factorization","volume":"401","author":"DD Lee","year":"1999","journal-title":"Nature"},{"issue":"10","key":"pcbi.1012391.ref027","doi-asserted-by":"crossref","first-page":"2756","DOI":"10.1162\/neco.2007.19.10.2756","article-title":"Projected gradient methods for nonnegative matrix factorization","volume":"19","author":"C-J Lin","year":"2007","journal-title":"Neural Comput"},{"key":"pcbi.1012391.ref028","article-title":"Algorithms for non-negative matrix factorization","volume":"13","author":"D Lee","journal-title":"Adv Neural Inf Process Syst"},{"key":"pcbi.1012391.ref029","first-page":"708","article-title":"Fast local algorithms for large scale nonnegative matrix and tensor factorizations.","author":"A Cichocki","year":"2009"},{"issue":"12","key":"pcbi.1012391.ref030","doi-asserted-by":"crossref","first-page":"4164","DOI":"10.1073\/pnas.0308531101","article-title":"Metagenes and molecular pattern discovery using matrix factorization","volume":"101","author":"J-P Brunet","year":"2004","journal-title":"Proc Natl Acad Sci U S A"},{"key":"pcbi.1012391.ref031","doi-asserted-by":"crossref","first-page":"309","DOI":"10.1016\/j.neunet.2021.03.020","article-title":"A statistical framework for non-negative matrix factorization based on generalized dual divergence","volume":"140","author":"K Devarajan","year":"2021","journal-title":"Neural Netw"},{"issue":"1","key":"pcbi.1012391.ref032","doi-asserted-by":"crossref","first-page":"137","DOI":"10.1007\/s10994-014-5470-z","article-title":"A unified statistical approach to non-negative matrix factorization and probabilistic latent semantic indexing","volume":"99","author":"K Devarajan","year":"2015","journal-title":"Mach Learn"},{"key":"pcbi.1012391.ref033","doi-asserted-by":"crossref","first-page":"50","DOI":"10.1145\/312624.312649","article-title":"Probabilistic latent semantic indexing.","volume-title":"Proceedings of the 22nd annual international ACM SIGIR conference on Research and development in information retrieval","author":"T Hofmann","year":"1999"},{"key":"pcbi.1012391.ref034","doi-asserted-by":"crossref","first-page":"41","DOI":"10.1186\/1471-2105-7-41","article-title":"Discovering semantic features in the literature: a foundation for building functional associations","volume":"7","author":"M Chagoyen","year":"2006","journal-title":"BMC Bioinformatics"},{"issue":"4","key":"pcbi.1012391.ref035","doi-asserted-by":"crossref","first-page":"447","DOI":"10.4025\/actasciagron.v37i4.19746","article-title":"Hierarchical genetic clusters for phenotypic analysis","volume":"37","author":"da Matta LB","year":"2015","journal-title":"Acta Sci Agron"},{"key":"pcbi.1012391.ref036","doi-asserted-by":"crossref","first-page":"63","DOI":"10.1016\/j.neucom.2014.08.095","article-title":"Hierarchical feature extraction by multi-layer non-negative matrix factorization network for classification task","volume":"165","author":"HA Song","year":"2015","journal-title":"Neurocomputing"},{"key":"pcbi.1012391.ref037","doi-asserted-by":"crossref","first-page":"305","DOI":"10.1016\/j.neucom.2021.08.152","article-title":"A survey of deep nonnegative matrix factorization","volume":"491","author":"W-S Chen","year":"2022","journal-title":"Neurocomputing"},{"key":"pcbi.1012391.ref038","article-title":"Deep NMF for speech separation. In: 2015 IEEE International Conference on Acoustics","author":"J Le Roux","year":"2015","journal-title":"Speech and Signal Processing (ICASSP)"},{"issue":"1","key":"pcbi.1012391.ref039","doi-asserted-by":"crossref","DOI":"10.1007\/s43670-023-00077-3","article-title":"Neural nonnegative matrix factorization for hierarchical multilayer topic modeling","volume":"22","author":"J Haddock","year":"2023","journal-title":"Sampl Theory Signal Process Data Anal"},{"issue":"9","key":"pcbi.1012391.ref040","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1009373","article-title":"Genetic dissection of complex traits using hierarchical biological knowledge","volume":"17","author":"H Tanaka","year":"2021","journal-title":"PLoS Comput Biol"},{"issue":"7","key":"pcbi.1012391.ref041","doi-asserted-by":"crossref","first-page":"1104","DOI":"10.1038\/jcbfm.2015.69","article-title":"Clustering-initiated factor analysis application for tissue classification in dynamic brain positron emission tomography","volume":"35","author":"R Boutchko","year":"2015","journal-title":"Journal of Cerebral Blood Flow & Metabolism"},{"issue":"23","key":"pcbi.1012391.ref042","doi-asserted-by":"crossref","first-page":"2684","DOI":"10.1093\/bioinformatics\/btn526","article-title":"Position-dependent motif characterization using non-negative matrix factorization","volume":"24","author":"LN Hutchins","year":"2008","journal-title":"Bioinformatics"},{"key":"pcbi.1012391.ref043","doi-asserted-by":"crossref","first-page":"367","DOI":"10.1186\/1471-2105-11-367","article-title":"A flexible R package for nonnegative matrix factorization","volume":"11","author":"R Gaujoux","year":"2010","journal-title":"BMC Bioinformatics"},{"key":"pcbi.1012391.ref044","doi-asserted-by":"crossref","first-page":"952916","DOI":"10.3389\/fpubh.2022.952916","article-title":"Surface charge changes in spike RBD mutations of SARS-CoV-2 and its variant strains alter the virus evasiveness via HSPGs: A review and mechanistic hypothesis","volume":"10","author":"Z Zhang","year":"2022","journal-title":"Front Public Health"},{"issue":"9","key":"pcbi.1012391.ref045","doi-asserted-by":"crossref","first-page":"831","DOI":"10.1038\/s41422-022-00700-3","article-title":"Hetero-bivalent nanobodies provide broad-spectrum protection against SARS-CoV-2 variants of concern including Omicron","volume":"32","author":"H Ma","year":"2022","journal-title":"Cell Res"},{"issue":"1","key":"pcbi.1012391.ref046","doi-asserted-by":"crossref","first-page":"19089","DOI":"10.1038\/s41598-022-23342-2","article-title":"Rapidly identifying new coronavirus mutations of potential concern in the Omicron variant using an unsupervised learning strategy","volume":"12","author":"LP Zhao","year":"2022","journal-title":"Sci Rep"},{"issue":"4","key":"pcbi.1012391.ref047","doi-asserted-by":"crossref","first-page":"225","DOI":"10.2217\/fvl-2022-0152","article-title":"The impact of mutation sets in receptor-binding domain of SARS-CoV-2 variants on stability of the RBD\u2013ACE2 complex","volume":"18","author":"M Peka","year":"2023","journal-title":"Future Virol"},{"issue":"6","key":"pcbi.1012391.ref048","doi-asserted-by":"crossref","first-page":"403","DOI":"10.1093\/procel\/pwae007","article-title":"Mutations in the SARS-CoV-2 spike receptor binding domain and their delicate balance between ACE2 affinity and antibody evasion","volume":"15","author":"S Xue","year":"2024","journal-title":"Protein Cell"},{"issue":"7898","key":"pcbi.1012391.ref049","doi-asserted-by":"crossref","first-page":"676","DOI":"10.1038\/s41586-021-04388-0","article-title":"Striking antibody evasion manifested by the Omicron variant of SARS-CoV-2","volume":"602","author":"L Liu","year":"2022","journal-title":"Nature"},{"issue":"12","key":"pcbi.1012391.ref050","doi-asserted-by":"crossref","first-page":"1437","DOI":"10.3390\/pathogens11121437","article-title":"SARS-CoV-2 antibody effectiveness is influenced by non-epitope mutation\/binding-induced denaturation of the epitope 3D architecture","volume":"11","author":"MM Malisheni","year":"2022","journal-title":"Pathogens"},{"issue":"9","key":"pcbi.1012391.ref051","doi-asserted-by":"crossref","DOI":"10.3390\/v14092009","article-title":"The rise and fall of SARS-CoV-2 variants and ongoing diversification of omicron","volume":"14","author":"T Wiegand","year":"2009","journal-title":"Viruses"},{"issue":"29","key":"pcbi.1012391.ref052","doi-asserted-by":"crossref","DOI":"10.1073\/pnas.2205784119","article-title":"A broad and potent neutralization epitope in SARS-related coronaviruses","volume":"119","author":"M Yuan","year":"2022","journal-title":"Proc Natl Acad Sci U S A"},{"issue":"1","key":"pcbi.1012391.ref053","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1010822","article-title":"Antibody accessibility determines location of spike surface mutations in SARS-CoV-2 variants","volume":"19","author":"S von B\u00fclow","year":"2023","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1012391.ref054","doi-asserted-by":"crossref","first-page":"39","DOI":"10.3389\/fmicb.2019.00039","article-title":"Influenza Virus Neuraminidase Structure and Functions","volume":"10","author":"JL McAuley","year":"2019","journal-title":"Front Microbiol"},{"issue":"1","key":"pcbi.1012391.ref055","doi-asserted-by":"crossref","first-page":"31","DOI":"10.1186\/s12985-023-01982-8","article-title":"Analysis of 3.5 million SARS-CoV-2 sequences reveals unique mutational trends with consistent nucleotide and codon frequencies","volume":"20","author":"SE Fumagalli","year":"2023","journal-title":"Virol J"},{"issue":"1","key":"pcbi.1012391.ref056","doi-asserted-by":"crossref","first-page":"12331","DOI":"10.1038\/s41598-020-69342-y","article-title":"Human SARS-CoV-2 has evolved to reduce CG dinucleotide in its open reading frames","volume":"10","author":"Y Wang","year":"2020","journal-title":"Sci Rep"},{"issue":"5","key":"pcbi.1012391.ref057","doi-asserted-by":"crossref","first-page":"1092","DOI":"10.1093\/molbev\/msy015","article-title":"Codon usage bias in animals: disentangling the effects of natural selection, effective population size, and gc-biased gene conversion","volume":"35","author":"N Galtier","year":"2018","journal-title":"Mol Biol Evol"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1012391","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2025,7,18]],"date-time":"2025-07-18T00:00:00Z","timestamp":1752796800000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1012391","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,7,18]],"date-time":"2025-07-18T17:55:04Z","timestamp":1752861304000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1012391"}},"subtitle":[],"editor":[{"given":"Jordan","family":"Douglas","sequence":"first","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2025,4,24]]},"references-count":57,"journal-issue":{"issue":"4","published-online":{"date-parts":[[2025,4,24]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1012391","relation":{},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2025,4,24]]}}}