Disease persistence in epidemiological models: The interplay between vaccination and migration

Abstract

We consider the interplay of vaccination and migration rates on disease persistence in epidemiological systems. We show that short-term and long-term migration can inhibit disease persistence. As a result, we show how migration changes how vaccination rates should be chosen to maintain herd immunity. In a system of coupled SIR models, we analyze how disease eradication depends explicitly on vaccine distribution and migration connectivity. The analysis suggests potentially novel vaccination policies that underscore the importance of optimal placement of finite resources.

Introduction

Countries are increasingly connected by travel and economics. Due to economic disparities and political turmoil, extreme heterogeneity exists in childhood vaccination coverage across the two sides of multiple national boundaries. It has been suggested that the immunization coverage of neighboring countries or those countries well connected by travel can or should be used when crafting national level immunization policy. In the case of hepatitis B, Gay and Edmunds [1] argue that it would be four times more cost effective for the United Kingdom to sponsor a vaccination program in Bangladesh than to introduce its own universal program. When indigenous wild poliovirus was eradicated in all but four endemic countries in 2005: India, Nigeria, Pakistan and Afghanistan, it was exported from northern Nigeria and northern India and subsequently caused >50 outbreaks and paralyzed >1500 children in previously polio-free countries across Asia and Africa [2]. And in 2007, the WHO estimated that there were 197,000 measles deaths, despite the 82% worldwide vaccination coverage. In countries where measles has been largely eliminated, cases imported from other countries remain an important source of infection [3]. It is clear that a country needs to be concerned with the vaccination rate of a neighboring country as well as its own.

On another scale, vaccination policies must also take into consideration the subpopulation dynamics within a country. Wilson et al. [4] models linked urban and rural epidemics of HIV and discusses how to optimize a limited treatment supply to minimize new infections. Cummings et al. [5] uses data to identify a distinct pattern in the periodicity of measles outbreaks in Cameroon before the widespread vaccination efforts of the Measles Initiative. The southern part of Cameroon experienced a significant measles epidemic approximately every three years. In contrast, the three northern provinces contend with annual measles epidemics. In 2000 and 2001, these cyclic outbreaks coincided, exacerbating the situation and causing a much more severe epidemic [5].

Noting that a small contribution of infections from one population to another could drive a new type of epidemic that would not normally occur, we study how migration between populations could change dynamics and respective herd immunity levels in metapopulation models. We analyze a model of a disease imported between subpopulations of a region by short-term and long-term migration with limited vaccination coverage. Our initial study is based on the analysis of a system of canonical SIR compartmental models. The system allows the rigorous proof of the qualitative affects migration has on herd immunity. The model can be enhanced to include more compartments or seasonal forcing, but most of these systems will require numerical exploration of trends in spatial synchrony and bifurcation analysis, which will be explored in future papers. In this article, we revisit the fundamental ways migration is modeled in metapopulation models and how it fundamentally affects herd immunity.

Migration is often treated as a phenomenological input to maintain incidence in a population that might experience local fade-out [6]. Long-term migration has been analyzed by Liebovitch and Schwartz [7], with a thorough derivation of the linear flux term coupling the patches. This approach also agrees with the classes of models proposed by Sattenspiel and Dietz [8] and Lloyd and Jansen [9]. Keeling and Rohani [10] investigated the spatial coupling of dynamics exhibited in models using multiple formulations of migration including mass-action coupling and linear flux terms. However, they did not explore the impact of coupling in the presence of vaccination needed to maintain disease free states. Additional analysis of mixed long-term and short-term migration in transport-related disease spread can be found in [11], [12], [13]. These papers derive the global asymptotic stability of the disease free state for a new disease. Because there is no vaccination, the papers conclude that it is essential to strengthen restrictions of passenger travel as soon as the infectious diseases appear.

Our paper considers how migration directly affects the vaccination levels needed for herd immunity against a known disease and how that would impact optimum usage of limited vaccination supplies. We investigate the dynamics of models that include mass-action coupling, an assumption that assumes mixing occurs at fast time scales, and linear migration, which is more consistent with mixing occurring at long time scales. The organization of this paper is as follows: We introduce a coupled compartmental model in Section2 and perform stability analysis of the disease free state as a function of the migration and vaccination rates. We also consider normal forms of the bifurcations created by the short-term and long-term migration dynamics. Section 3 describes how vaccination rates should be adjusted with respect to short-term and long-term migration levels to preserve herd immunity. Section 4 has a summary of our observations and conclusions.