Disentangling how multiple traits drive 2 strain frequencies in SIS dynamics with coinfection

A general theory for competitive dynamics among many strains at the epidemiological level is required to understand polymorphisms in virulence, transmissibility, antibiotic resistance and other biological traits of infectious agents. Mathematical coinfection models have addressed specific systems, focusing on the criteria leading to stable coexistence or competitive exclusion, however, due to their complexity and nonlinearity, analytical solutions in coinfection models remain rare. Here we study a 2-strain SIS compartmental model with co-infection/co-colonization, incorporating multiple fitness dimensions under the same framework: variation in transmissibility, duration of carriage, pairwise susceptibilities to coinfection, coinfection duration, and transmission priority effects from mixed coinfection. Taking advantage of a singular perturbation approach, under the assumption of strain similarity, we expose how strain dynamics on a slow timescale are explicitly governed by a replicator equation which encapsulates all traits and their interplay. This allows us to predict explicitly not only the final epidemiological outcome of a given 2-player competition, but moreover, their entire frequency dynamics as a direct function of their relative variation and of strain-transcending global parameters. Based on mutual invasion fitnesses, we analyze and report rigorous results on transition phenomena in the 2-strain system, strongly mediated via coinfection prevalence. We show that coinfection is not always a promoter of coexistence; instead, its effect to favour or prevent polymorphism is non-monotonic and depends on the type and level of phenotypic differentiation between strains. This framework offers a deeper analytical understanding of 2-strain competitive games in coinfection, with theoretical and practical applications in epidemiology, ecology and evolution.