Dynamics of a multigroup epidemiological model with group-targeted vaccination strategies

Optimal vaccination strategies for a heterogeneous population using multiple objectives: The case of L1- and L2-formulations

The choice of the objective functional in optimization problems coming from biomedical and epidemiological applications plays a key role in optimal control outcomes. In this study, we investigate the role of the objective functional on the structure of the optimal control solution for an epidemic model for sexually transmitted infections that includes a core group with higher sexual activity levels than the rest of the population. An optimal control problem is formulated to find a targeted vaccination program able to control the spread of the infection with minimum vaccine deployment. Both L1- and L2-objectives are considered as an attempt to explore the trade-offs between control dynamics and the functional form characterizing optimality. The results show that the optimal vaccination policies for both the L1- and the L2-formulation share one important qualitative property, that is, immunization of the core group should be prioritized by policymakers to achieve a fast reduction of the epidemic. However, quantitative aspects of this result can be significantly affected depending on the choice of the control weights between formulations. Overall, the results suggest that with appropriate weight constants, the optimal control outcomes are reasonably robust with respect to the L1- or L2-formulation. This is particularly true when the monetary cost of the control policy is substantially lower than the cost associated with the disease burden. Under these conditions, even if the L1-formulation is more realistic from a modeling perspective, the L2-formulation can be used as an approximation and yield qualitatively comparable outcomes.

Implications for infectious disease models of heterogeneous mixing on control thresholds

Mixing among sub-populations, as well as heterogeneity in characteristics affecting their reproduction numbers, must be considered when evaluating public health interventions to prevent or control infectious disease outbreaks. In this overview, we apply a linear algebraic approach to re-derive some well-known results pertaining to preferential within- and proportionate among-group contacts in compartmental models of pathogen transmission. We give results for the meta-population effective reproduction number ([Formula: see text]) assuming different levels of vaccination in the sub-populations. Specifically, we unpack the dependency of [Formula: see text] on the fractions of contacts reserved for individuals within one's own subgroup and, by obtaining implicit expressions for the partial derivatives of [Formula: see text], we show that these increase as this preferential-mixing fraction increases in any sub-population.

Structural identifiability analysis of age-structured PDE epidemic models

Computational and mathematical models rely heavily on estimated parameter values for model development. Identifiability analysis determines how well the parameters of a model can be estimated from experimental data. Identifiability analysis is crucial for interpreting and determining confidence in model parameter values and to provide biologically relevant predictions. Structural identifiability analysis, in which one assumes data to be noiseless and arbitrarily fine-grained, has been extensively studied in the context of ordinary differential equation (ODE) models, but has not yet been widely explored for age-structured partial differential equation (PDE) models. These models present additional difficulties due to increased number of variables and partial derivatives as well as the presence of boundary conditions. In this work, we establish a pipeline for structural identifiability analysis of age-structured PDE models using a differential algebra framework and derive identifiability results for specific age-structured models. We use epidemic models to demonstrate this framework because of their wide-spread use in many different diseases and for the corresponding parallel work previously done for ODEs. In our application of the identifiability analysis pipeline, we focus on a Susceptible-Exposed-Infected model for which we compare identifiability results for a PDE and corresponding ODE system and explore effects of age-dependent parameters on identifiability. We also show how practical identifiability analysis can be applied in this example.

Global dynamical behavior of a multigroup SVIR epidemic model with Markovian switching

In this paper, we are concerned with the global dynamical behavior of a multigroup SVIR epidemic model, which is formulated as a piecewise-deterministic Markov process. We first obtain sufficient criteria for extinction of the diseases. Then we establish sufficient criteria for persistence in the mean of the diseases. Moreover, in the case of persistence, we find a domain which is positive recurrence for the solution of the stochastic system by constructing an appropriate Lyapunov function with regime switching.

Modelling epidemics with fractional-dose vaccination in response to limited vaccine supply

The control strategies of emergency infectious diseases are constrained by limited medical resources. The fractional dose vaccination strategy as one of feasible strategies was proposed in response to global shortages of vaccine stockpiles. Although a variety of epidemic models have been developed under the circumstances of limited resources in treatment, few models particularly investigated vaccination strategies in resource-limited settings. In this paper, we develop a two-group SIR model with incorporation of proportionate mixing patterns and n-fold fractional dose vaccination related parameters to evaluate the efficiency of fractional dose vaccination on disease control at the population level. The existence and uniqueness of the final size of the two-group SIR epidemic model, the formulation of the basic reproduction number and the relationship between them are established. Moreover, numerical simulations are performed based on this two-group vector-free model to investigate the effectiveness of n-fold fractional dose vaccination by using the emergency outbreaks of yellow fever in Angola in 2016. By employing linear and nonlinear dose-response relationships, we compare the resulting fluctuations of four characteristics of the epidemics, which are the outbreak size, the peak time of the outbreak, the basic reproduction number and the infection attack rate (IAR). For both types of dose-response relationships, dose-fractionation takes positive effects in lowering the outbreak size, delay the peak time of the outbreak, reducing the basic reproduction number and the IAR of yellow fever only when the vaccine efficacy is high enough. Moreover, five-fold fractional dose vaccination strategy may not be the optimal vaccination strategy as proposed by the World Health Organization if the dose-response relationship is nonlinear.

A general theory for target reproduction numbers with applications to ecology and epidemiology

A general framework for threshold parameters in population dynamics is developed using the concept of target reproduction numbers. This framework identifies reproduction numbers and other threshold parameters in the literature in terms of their roles in population control. The framework is applied to the analysis of single and multiple control strategies in ecology and epidemiology, and this provides new biological insights.

Global dynamics of a mathematical model for the possible re-emergence of polio

Motivated by studies warning about a possible re-emergence of poliomyelitis in Europe, we analyse a compartmental model for the transmission of polio describing the possible effect of unvaccinated people arriving to a region with low vaccination coverage. We calculate the basic reproduction number, and determine the global dynamics of the system: we show that, depending on the parameters, one of the two equilibria is globally asymptotically stable. The main tools applied are Lyapunov functions and persistence theory. We illustrate the analytic results by numerical examples, which also suggest that in order to avoid the risk of polio re-emergence, vaccinating the immigrant population might result insufficient, and also the vaccination coverage of countries with low rates should be increased.

Global threshold dynamics of SIQS epidemic model in time fluctuating environment

The paper characterizes the global threshold dynamics of an epidemic model of SIQS type in environments with fluctuations, where the quarantine class is explicitly involved. Criteria are established for the permanence and extinction of the infective in environments with time oscillations. In particular, we further consider an environment which varies periodically in time. The global threshold dynamic scenarios i.e. the existence and global asymptotic stability of the disease-free periodic solution, the existence of the endemic periodic solution and the permanence of the infective are completely characterized by the basic reproduction number defined by the spectral radius of an associated linear integral operator.

Assessing the Potential Impact of Hormonal-Based Contraceptives on HIV Transmission Dynamics Among Heterosexuals

HIV susceptibility linked to hormonal contraception (HC) has been studied before, but with mixed results. Reports from some of the recent findings have prompted the World Health Organisation to encourage women who use HC to concurrently use condoms in order to prevent HIV infection in the light of possible increased HIV risk of infection associated with hormone-based contraceptives. A two-sex HIV model classifying women into three risk groups consisting of individuals who use condoms, natural methods, and hormone-based contraceptives is formulated and analysed to assess the possible effects of various birth control strategies on the transmission dynamics of the disease. Our model results showed that women who use HC could be key drivers of the epidemic and that their increased infectivity may be critical in driving the epidemic. Women who use hormone-based contraceptives potentially act as a core group from which men get infected and in turn transmit the disease to other population groups. We fitted the model to HIV prevalence data for Zimbabwe reported by UNAIDS and Zimbabwe Ministry of Health and Child Care and used the model fit to project HIV prevalence. Predictions using HIV data for Zimbabwe suggest that a hypothesised increase in susceptibility and infectivity of two-, three-, and fourfold would result in a 25, 50, and 100% increase in baseline HIV prevalence projection, respectively, thus suggesting possible increased disease burden even in countries reporting plausible HIV prevalence declines. Although a possible causal relationship between HIV susceptibility and HC use remains subject of continuing scientific probe, its inclusion as part of birth control strategy has been shown in this study, to possibly increase HIV transmission. If proven, HC use may potentially explain the inordinate spread of HIV within the sub-Saharan Africa region and therefore compel for urgent assessment with a view to reorienting birth control methods in use in settings with generalised epidemics.

An elaboration of theory about preventing outbreaks in homogeneous populations to include heterogeneity or preferential mixing

The goal of many vaccination programs is to attain the population immunity above which pathogens introduced by infectious people (e.g., travelers from endemic areas) will not cause outbreaks. Using a simple meta-population model, we demonstrate that, if sub-populations either differ in characteristics affecting their basic reproduction numbers or if their members mix preferentially, weighted average sub-population immunities cannot be compared with the proportionally-mixing homogeneous population-immunity threshold, as public health practitioners are wont to do. Then we review the effect of heterogeneity in average per capita contact rates on the basic meta-population reproduction number. To the extent that population density affects contacts, for example, rates might differ in urban and rural sub-populations. Other differences among sub-populations in characteristics affecting their basic reproduction numbers would contribute similarly. In agreement with more recent results, we show that heterogeneous preferential mixing among sub-populations increases the basic meta-population reproduction number more than homogeneous preferential mixing does. Next we refine earlier results on the effects of heterogeneity in sub-population immunities and preferential mixing on the effective meta-population reproduction number. Finally, we propose the vector of partial derivatives of this reproduction number with respect to the sub-population immunities as a fundamentally new tool for targeting vaccination efforts.

The efficiency of targeted intervention in limiting the spread of HIV and Hepatitis C Virus among injecting drug users

BACKGROUND

Interventions aimed at minimizing the spread of blood borne infections among Injecting Drug Users (IDU) are impeded by limitations in resources. To enhance their efficiency, it may be beneficial to target specific behavioural subpopulations, distinguished by syringe sharing tendencies.

METHODS

We used mathematical modelling to explore the effects of two types of intervention: removal of individuals from the injecting population and risk decrease at group-level (e.g. distribution of syringes). We computed the direct effects of intervention on the probability of obtaining and spreading infection as a function of baseline risk behaviour. Population level effects of (targeted) intervention were explored using a differential equations model, which incorporated two levels of risk.

RESULTS

Within most scenarios of risk distribution considered, HIV could be substantially reduced or eliminated by targeting high risk IDU only. Conversely, higher incidence reductions for HCV were reached in many scenarios when targeting low risk IDU. The potential for preventing infections by removal of uninfected IDU increases with baseline risk, but so does the probability that an IDU is already infected before being reached by intervention. Decreasing risk is likely to only delay rather than prevent infection for IDU borrowing many syringes, especially for a very infectious disease such as HCV.

CONCLUSIONS

The efficiency of intervention on injecting drug users may be much enhanced by targeting specific risk subgroups. However, the optimal targeting policy depends strongly on the infection under consideration.

Extending the type reproduction number to infectious disease control targeting contacts between types

A new quantity called the target reproduction number is defined to measure control strategies for infectious diseases with multiple host types such as waterborne, vector-borne and zoonotic diseases. The target reproduction number includes as a special case and extends the type reproduction number to allow disease control targeting contacts between types. Relationships among the basic, type and target reproduction numbers are established. Examples of infectious disease models from the literature are given to illustrate the use of the target reproduction number.