Projecting social contact matrices to different demographic structures

Abstract The modeling of large-scale communicable epidemics has greatly benefited in the last years from the increasing availability of highly detailed data. Particularly, in order to achieve quantitative descriptions of the evolution of epidemics, contact networks and mixing patterns are key. These heterogeneous patterns depend on several factors such as location, socioeconomic conditions, time, and age. This last factor has been shown to encapsulate a large fraction of the observed inter-individual variation in contact patterns, an observation validated by different measurements of age-dependent contact matrices. Recently, several works have studied how to project those matrices to areas where empiric data is not available. However, the dependence of contact matrices on demographic structures and their time evolution has been largely neglected. In this work, we tackle the problem of how to transform an empirical contact matrix that has been obtained for a given demographic structure into a different contact matrix that is compatible with a different demography. The methodology discussed here allows extrapolating a contact structure measured in a particular area to any other whose demographic structure is known, as well as to obtain the time evolution of contact matrices as a function of the demographic dynamics of the populations they refer to. To quantify the effect of considering time-dynamics of contact patterns on disease modeling, we implemented a Susceptible-Exposed-Infected-Recovered (SEIR) model on 16 different countries and evaluated the impact of neglecting the temporal evolution of mixing patterns. Our results show that simulated disease incidence rates, both at the aggregated and age-specific levels, are significantly dependent on contact structures variation driven by demographic evolution. The present work opens the path to eliminate technical biases from model-based impact evaluations of future epidemic threats and warns against the use of contact matrices to model diseases without correcting for demographic evolution or geographic variations. Author summary Large scale epidemic outbreaks represent an ever increasing threat to humankind. In order to anticipate eventual pandemics, mathematical modeling should not only have the capacity to model in real time an ongoing disease, but also to predict the evolution of potential outbreaks in different locations and times. To this end, computational frameworks need to incorporate, among other ingredients, realistic contact patterns into the models. This not only implies anticipating the demographic structure of the populations under study, but also understanding how demographic evolution reshapes social mixing patterns along time. Here we present a mathematical framework to solve this problem and test our modeling approach on 16 different empirical contact matrices. We also evaluate the impact of an eventual future outbreak by simulating a SEIR scenario in the countries analyzed. Our results show that using outdated or imported contact matrices that do not take into account demographic structure or its evolution can lead to largely misleading conclusions.