We propose a method to estimate pairwise strain interactions from population-level frequencies across different endemic settings. We apply the framework of replicator dynamics, derived from a multi-strain SIS model with co-colonization, to extract from 5 datasets the fundamental backbone of strain interactions. In our replicator, each pairwise invasion fitness explicitly arises from local environmental context and trait variations between strains. We adopt the simplest formulation for multi-strain coexistence, where context is encoded in basic reproduction number R 0 and mean global susceptibility to co-colonization k, and trait variations α ij capture pairwise deviations from k. We integrate Streptococcus pneumoniae serotype frequencies and serotype identities collected from 5 environments: epidemiological surveys in Denmark, Nepal, Iran, Brazil and Mozambique, and mechanistically link their distributions. Our results have twofold implications. First, we offer a new proof-of-concept in the inference of multi-species interactions based on cross-sectional data. We also discuss 2 key aspects of the method: the site ordering for sequential fitting, and stability constraints on the dynamics. Secondly, we effectively estimate at high-resolution more than 70% of the 92 × 92 pneumococcus serotype interaction matrix in co-colonization, allowing for further projections and hypotheses testing. We show that, in these bacteria, both within- and between- serotype interaction coefficients' distribution emerge to be unimodal, their difference in mean broadly reflecting stability assumptions on serotype coexistence. This framework enables further model calibration to global data: cross-sectional across sites, or longitudinal in one site over time, - and should allow a more robust and integrated investigation of intervention effects in such biodiverse ecosystems.

Center for Computational and Stochastic Mathematics Instituto Superior Tecnico Lisbon Portugal

Department of Mathematics and Statistics Masaryk University Brno Czech Republic

Institut Denis Poisson University of Tours Tours France

Zobrazit více v PubMed

Adebanjo T, Lessa FC, Mucavele H, Moiane B, Chauque AP, Massora S, Carvalho MDG, Whitney CG, Sigauque B (2018) Pneumococcal carriage and serotype distribution among children with and without pneumonia in Mozambique, 2014–2016. PloS one 13(6):e0199363 PubMed PMC

Aguas R, Dorigatti I, Coudeville L, Luxemburger C, Ferguson N (2019) Cross-serotype interactions and disease outcome prediction of dengue infections in Vietnam. Sci Rep 9(1):9395 PubMed PMC

Akman D, Akman O, Schaefer E (2018) Parameter estimation in ordinary differential equations modeling via particle swarm optimization. J Appl Math 2018(1):9160793

Alizon S (2013) Co-infection and super-infection models in evolutionary epidemiology. Interface Focus 3(6):20130031 PubMed PMC

Alizon S, de Roode JC, Michalakis Y (2013) Multiple infections and the evolution of virulence. Ecol Lett 16(4):556–567 PubMed

Alizon S, Murall CL, Saulnier E, Sofonea MT (2019) Detecting within-host interactions from genotype combination prevalence data. Epidemics 29:100349 PubMed PMC

Allesina S, Levine JM (2011) A competitive network theory of species diversity. Proc Natl Acad Sci 108(14):5638–5642 PubMed PMC

Allesina S, Tang S (2012) Stability criteria for complex ecosystems. Nature 483(7388):205–208 PubMed

Alshawaqfeh M, Serpedin E, Younes AB (2017) Inferring microbial interaction networks from metagenomic data using SgLV-EKF algorithm. BMC Genomics 18:1–12 PubMed PMC

Aubin J-P, Sigmund K (1988) Permanence and viability. J Comput Appl Math 22(2–3):203–209

Auranen K, Mehtälä J, Tanskanen A, Kaltoft MS (2010) Between-strain competition in acquisition and clearance of pneumococcal carriage? epidemiologic evidence from a longitudinal study of day-care children. Am J Epidemiol 171(2):169–176 PubMed PMC

Bedford T, Rambaut A, Pascual M (2012) Canalization of the evolutionary trajectory of the human influenza virus. BMC Biol 10:1–12 PubMed PMC

Bentley SD, Aanensen DM, Mavroidi A, Saunders D, Rabbinowitsch E, Collins M, Donohoe K, Harris D, Murphy L, Quail MA et al (2006) Genetic analysis of the capsular biosynthetic locus from all 90 pneumococcal serotypes. PLoS Genet 2(3):e31 PubMed PMC

Chaturvedi AK, Katki HA, Hildesheim A, Rodríguez AC, Quint W, Schiffman M, Van Doorn L-J, Porras C, Wacholder S, Gonzalez P et al (2011) Human papillomavirus infection with multiple types: pattern of coinfection and risk of cervical disease. J Infect Dis 203(7):910–920 PubMed PMC

Chesson P (1994) Multispecies competition in variable environments. Theor Popul Biol 45(3):227–276

Cobey S, Baskerville EB, Colijn C, Hanage W, Fraser C, Lipsitch M (2017) Host population structure and treatment frequency maintain balancing selection on drug resistance. J R Soc Interface 14(133):20170295 PubMed PMC

Cody ML, Diamond JM (1975) Ecology and evolution of communities. Harvard University Press, Cambridge

Coyte KZ, Schluter J, Foster KR (2015) The ecology of the microbiome: networks, competition, and stability. Science 350(6261):663–666 PubMed

Dekaj E, Gjini E (2024) Pneumococcus and the stress-gradient hypothesis: a trade-off links R0 and susceptibility to co-colonization across countries. Theor Popul Biol 156:77–92 PubMed

Eberhart R, Kennedy J (1995) Particle swarm optimization. In Proceedings of the IEEE international conference on neural networks, volume 4, pages 1942–1948. Citeseer

Garcia Quesada M, Peterson ME, Bennett JC, Hayford K, Zeger SL, Yang Y, Hetrich MK, Feikin DR, Von Gottberg A, van der Linden M et al (2021) Serotype distribution of remaining pneumococcal meningitis in the mature PCV10/13 period: findings from the PSERENADE project. Microorganisms 9(4):738 PubMed PMC

Geno KA, Gilbert GL, Song JY, Skovsted IC, Klugman KP, Jones C, Konradsen HB, Nahm MH (2015) Pneumococcal capsules and their types: past, present, and future. Clin Microbiol Rev 28(3):871–899 PubMed PMC

Gjini E (2017) Geographic variation in pneumococcal vaccine efficacy estimated from dynamic modeling of epidemiological data post-PCV7. Sci Rep 7(1):1–16 PubMed PMC

Gjini E, Madec S (2021) The ratio of single to co-colonization is key to complexity in interacting systems with multiple strains. Ecol Evol 11(13):8456–8474 PubMed PMC

Gjini E, Madec S (2023) Towards a mathematical understanding of invasion resistance in multispecies communities. R Soc Open Sci 10(11):231034 PubMed PMC

Gjini E, Valente C, Sa-Leao R, Gomes MGM (2016) How direct competition shapes coexistence and vaccine effects in multi-strain pathogen systems. J Theor Biol 388:50–60 PubMed

Gog J, Swinton J (2002) A status-based approach to multiple strain dynamics. J Math Biol 44(2):169–184 PubMed

Grilli J (2020) Macroecological laws describe variation and diversity in microbial communities. Nat Commun 11(1):1–11 PubMed PMC

Harboe ZB, Slotved H-C, Konradsen HB, Kaltoft MS (2012) A pneumococcal carriage study in Danish pre-school children before the introduction of pneumococcal conjugate vaccination. Open Microbiol J 6:40 PubMed PMC

Hofbauer J, Sigmund K (1998) Evolutionary games and population dynamics. Cambridge University Press, Cambridge

Hu J, Amor DR, Barbier M, Bunin G, Gore J (2022) Emergent phases of ecological diversity and dynamics mapped in microcosms. Science 378(6615):85–89 PubMed

Hubbell SP (2011) The unified neutral theory of biodiversity and biogeography (MPB-32). Princeton University Press, New Jersey

Hui C, Richardson DM, Landi P, Minoarivelo HO, Roy HE, Latombe G, Jing X, CaraDonna PJ, Gravel D, Beckage B et al (2021) Trait positions for elevated invasiveness in adaptive ecological networks. Biol Invasions 23:1965–1985

Imöhl M, Reinert RR, van der Linden M (2010) Regional differences in serotype distribution, pneumococcal vaccine coverage, and antimicrobial resistance of invasive pneumococcal disease among German federal states. Int J Med Microbiol 300(4):237–247 PubMed

Jacquet C, Moritz C, Morissette L, Legagneux P, Massol F, Archambault P, Gravel D (2016) No complexity-stability relationship in empirical ecosystems. Nat Commun 7(1):12573 PubMed PMC

Kandasamy R, Gurung M, Thapa A, Ndimah S, Adhikari N, Murdoch DR, Kelly DF, Waldron DE, Gould KA, Thorson S et al (2015) Multi-serotype pneumococcal nasopharyngeal carriage prevalence in vaccine naive Nepalese children, assessed using molecular serotyping. PLoS ONE 10(2):e0114286 PubMed PMC

Koelle K, Cobey S, Grenfell B, Pascual M (2006) Epochal evolution shapes the phylodynamics of interpandemic influenza a (H3N2) in humans. Science 314(5807):1898–1903 PubMed

Le TMT, Madec S (2025) Spatiotemporal evolution of coinfection dynamics: a reaction-diffusion model. J Dyn Diff Equat 37(1):317–362

Le TMT, Madec S, Gjini E (2022) Disentangling how multiple traits drive 2 strain frequencies in SIS dynamics with coinfection. J Theor Biol 538:111041 PubMed

Le TMT, Gjini E, Madec S (2023) Quasi-neutral dynamics in a coinfection system with N strains and asymmetries along multiple traits. J Math Biol 87(3):48 PubMed

Lewnard JA, Hanage WP (2019) Making sense of differences in pneumococcal serotype replacement. Lancet Infect Dis 19(6):e213–e220 PubMed

Lin J, Andreasen V, Levin SA (1999) Dynamics of influenza a drift: the linear three-strain model. Math Biosci 162(1–2):33–51 PubMed

Lipsitch M, Colijn C, Cohen T, Hanage WP, Fraser C (2009) No coexistence for free: neutral null models for multistrain pathogens. Epidemics 1(1):2–13 PubMed PMC

Lipsitch M, Abdullahi O, D’Amour A, Xie W, Weinberger DM, Tchetgen ET, Scott JAG (2012) Estimating rates of carriage acquisition and clearance and competitive ability for pneumococcal serotypes in Kenya with a Markov transition model. Epidemiology (Cambridge, Mass.), 23(4):510 PubMed PMC

Lo C, Marculescu R (2017) Inferring microbial interactions from metagenomic time-series using prior biological knowledge. In Proceedings of the 8th ACM international conference on bioinformatics, computational biology, and health informatics, pages 168–177

MacArthur R (1960) On the relative abundance of species. Am Nat 94(874):25–36

Madec S, Gjini E (2020) Predicting n-strain coexistence from co-colonization interactions: epidemiology meets ecology and the replicator equation. Bull Math Biol 82:142 PubMed PMC

Man I, Auranen K, Wallinga J, Bogaards JA (2019) Capturing multiple-type interactions into practical predictors of type replacement following human papillomavirus vaccination. Philos Trans R Soc B 374(1773):20180298 PubMed PMC

Man I, Benincà E, Kretzschmar ME, Bogaards JA (2023) Reconstructing multi-strain pathogen interactions from cross-sectional survey data via statistical network inference. J R Soc Interface 20(205):20220912 PubMed PMC

Marino S, Baxter NT, Huffnagle GB, Petrosino JF, Schloss PD (2014) Mathematical modeling of primary succession of murine intestinal microbiota. Proc Natl Acad Sci 111(1):439–444 PubMed PMC

May RM (1972) Will a large complex system be stable? Nature 238(5364):413–414 PubMed

McCann KS (2000) The diversity-stability debate. Nature 405(6783):228–233 PubMed

Mejlhede N, Pedersen BV, Frisch M, Fomsgaard A (2010) Multiple human papilloma virus types in cervical infections: competition or synergy? APMIS 118(5):346–352 PubMed

Moore MR (2009) Rethinking replacement and resistance PubMed

Motomura I (1932) On the statistical treatment of communities. Zoological Magazine, Tokyo 44:379–383

Raß A, Schmitt M, Wanka R (2015) Explanation of stagnation at points that are not local optima in particle swarm optimization by potential analysis. In Proceedings of the companion publication of the 2015 annual conference on genetic and evolutionary computation, GECCO companion ’15, page 1463-1464. Association for Computing Machinery

Rodrigues H, Pinto T, Barros R, Teixeira L, Neves F (2017) Pneumococcal nasopharyngeal carriage among children in Brazil prior to the introduction of the 10-valent conjugate vaccine: a culture-and pcr-based survey. Epidemiol Infection 145(8):1720–1726 PubMed PMC

Shrestha S, Foxman B, Weinberger DM, Steiner C, Viboud C, Rohani P (2013) Identifying the interaction between influenza and pneumococcal pneumonia using incidence data. Sci Translational Med 5(191):19ra184-19ra184 PubMed PMC

Stein RR, Bucci V, Toussaint NC, Buffie CG, Rätsch G, Pamer EG, Sander C, Xavier JB (2013) Ecological modeling from time-series inference: insight into dynamics and stability of intestinal microbiota. PLoS Comput Biol 9(12):e1003388 PubMed PMC

Tabatabaei SR, Fallah F, Shiva F, Shamshiri A, Hajia M, Navidinia M, Karimi A, Rahbar M (2014) Multiplex pcr assay for detection of pneumococcal serotypes in nasopharyngeal samples of healthy children; Tehran, 2009-2010. Annual Research & Review in Biology, pages 3780–3790

Tokita K (2023) May’s dream: the interactions in biological networks were random matrices after all. In: The eighth international symposium on biocomplexity

Tokita K (2004) Species abundance patterns in complex evolutionary dynamics. Phys Rev Lett 93(17):178102 PubMed

Tokita K (2006) Statistical mechanics of relative species abundance. Eco Inform 1(3):315–324

van Baalen M, Sabelis MW (1995) The dynamics of multiple infection and the evolution of virulence. Am Nat 146(6):881–910

Weinberger DM, Trzciński K, Lu Y-J, Bogaert D, Brandes A, Galagan J, Anderson PW, Malley R, Lipsitch M (2009) Pneumococcal capsular polysaccharide structure predicts serotype prevalence. PLoS Pathog 5(6):e1000476 PubMed PMC

Wikramaratna PS, Kucharski A, Gupta S, Andreasen V, McLean AR, Gog JR (2015) Five challenges in modelling interacting strain dynamics. Epidemics 10:31–34 PubMed

Xu Y, Zhou X, Zheng W, Cui B, Xie C, Liu Y, Qin X, Liu J (2024) Serotype distribution, antibiotic resistance, multilocus sequence typing, and virulence factors of invasive and non-invasive streptococcus pneumoniae in Northeast China from 2000 to 2021. Med Microbiol Immunol 213(1):12 PubMed

Yoshino Y, Galla T, Tokita K (2008) Rank abundance relations in evolutionary dynamics of random replicators. Phys Rev E-Stat Nonlinear Soft Matter Phys 78(3):031924 PubMed