Dynamics and Control of Diseases in Networks with Community Structure

The dynamics of infectious diseases spread via direct person-to-person transmission (such as influenza, smallpox, HIV/AIDS, etc.) depends on the underlying host contact network. Human contact networks exhibit strong community structure. Understanding how such community structure affects epidemics may provide insights for preventing the spread of disease between communities by changing the structure of the contact network through pharmaceutical or non-pharmaceutical interventions. We use empirical and simulated networks to investigate the spread of disease in networks with community structure. We find that community structure has a major impact on disease dynamics, and we show that in networks with strong community structure, immunization interventions targeted at individuals bridging communities are more effective than those simply targeting highly connected individuals. Because the structure of relevant contact networks is generally not known, and vaccine supply is often limited, there is great need for efficient vaccination algorithms that do not require full knowledge of the network. We developed an algorithm that acts only on locally available network information and is able to quickly identify targets for successful immunization intervention. The algorithm generally outperforms existing algorithms when vaccine supply is limited, particularly in networks with strong community structure. Understanding the spread of infectious diseases and designing optimal control strategies is a major goal of public health. Social networks show marked patterns of community structure, and our results, based on empirical and simulated data, demonstrate that community structure strongly affects disease dynamics. These results have implications for the design of control strategies.     

Adaptive and optimized COVID-19 vaccination strategies across geographical regions and age groups

We evaluate the efficiency of various heuristic strategies for allocating vaccines against COVID-19 and compare them to strategies found using optimal control theory. Our approach is based on a mathematical model which tracks the spread of disease among different age groups and across different geographical regions, and we introduce a method to combine age-specific contact data to geographical movement data. As a case study, we model the epidemic in the population of mainland Finland utilizing mobility data from a major telecom operator. Our approach allows to determine which geographical regions and age groups should be targeted first in order to minimize the number of deaths. In the scenarios that we test, we find that distributing vaccines demographically and in an age-descending order is not optimal for minimizing deaths and the burden of disease. Instead, more lives could be saved by using strategies which emphasize high-incidence regions and distribute vaccines in parallel to multiple age groups. The level of emphasis that high-incidence regions should be given depends on the overall transmission rate in the population. This observation highlights the importance of updating the vaccination strategy when the effective reproduction number changes due to the general contact patterns changing and new virus variants entering.     

Cancer immunotherapy, mathematical modeling and optimal control

Clinical immunologists, among other problems, routinely face a question: what is the best time and dose for a certain therapeutic agent to be administered to the patient in order to decrease/eradicate the pathological condition? In cancer immunotherapies the therapeutic agent is something able to elicit an immune response against cancer. The immune response has its own dynamics that depends on the immunogenicity of the therapeutic agent and on the duration of the immune response. The question then is "how can we decide when and how much of the drug to inject so to have a prolonged and effective immune response to the cancer?". This question can be addressed in mathematical terms in two stages: first one construct a mathematical model describing the cancer-immune interaction and secondly one applies the theory of optimal control to determine when and to which extent to stimulate the immune system by means of an immunotherapeutic agent administered in discrete variable doses within the therapeutic period. The solution of this mathematical problem is described and discussed in this article. We show that the method employed can be applied to find the optimal protocol in a variety of clinical problems where the kinetics of the drug or treatment and its influence on the physiologic/pathologic functions have been described by a system of ordinary differential equations.     

Cost-Effective Control of Plant Disease When Epidemiological Knowledge Is Incomplete: Modelling Bahia Bark Scaling of Citrus

A spatially-explicit, stochastic model is developed for Bahia bark scaling, a threat to citrus production in north-eastern Brazil, and is used to assess epidemiological principles underlying the cost-effectiveness of disease control strategies. The model is fitted via Markov chain Monte Carlo with data augmentation to snapshots of disease spread derived from a previously-reported multi-year experiment. Goodness-of-fit tests strongly supported the fit of the model, even though the detailed etiology of the disease is unknown and was not explicitly included in the model. Key epidemiological parameters including the infection rate, incubation period and scale of dispersal are estimated from the spread data. This allows us to scale-up the experimental results to predict the effect of the level of initial inoculum on disease progression in a typically-sized citrus grove. The efficacies of two cultural control measures are assessed: altering the spacing of host plants, and roguing symptomatic trees. Reducing planting density can slow disease spread significantly if the distance between hosts is sufficiently large. However, low density groves have fewer plants per hectare. The optimum density of productive plants is therefore recovered at an intermediate host spacing. Roguing, even when detection of symptomatic plants is imperfect, can lead to very effective control. However, scouting for disease symptoms incurs a cost. We use the model to balance the cost of scouting against the number of plants lost to disease, and show how to determine a roguing schedule that optimises profit. The trade-offs underlying the two optima we identify—the optimal host spacing and the optimal roguing schedule—are applicable to many pathosystems. Our work demonstrates how a carefully parameterised mathematical model can be used to find these optima. It also illustrates how mathematical models can be used in even this most challenging of situations in which the underlying epidemiology is ill-understood.     

Modeling the Dynamics of a Non-Limited and a Self-Limited Gene Drive System in Structured Aedes aegypti Populations

Recently there have been significant advances in research on genetic strategies to control populations of disease-vectoring insects. Some of these strategies use the gene drive properties of selfish genetic elements to spread physically linked anti-pathogen genes into local vector populations. Because of the potential of these selfish elements to spread through populations, control approaches based on these strategies must be carefully evaluated to ensure a balance between the desirable spread of the refractoriness-conferring genetic cargo and the avoidance of potentially unwanted outcomes such as spread to non-target populations. There is also a need to develop better estimates of the economics of such releases. We present here an evaluation of two such strategies using a biologically realistic mathematical model that simulates the resident Aedes aegypti mosquito population of Iquitos, Peru. One strategy uses the selfish element Medea, a non-limited element that could permanently spread over a large geographic area; the other strategy relies on Killer-Rescue genetic constructs, and has been predicted to have limited spatial and temporal spread. We simulate various operational approaches for deploying these genetic strategies, and quantify the optimal number of released transgenic mosquitoes needed to achieve definitive spread of Medea-linked genes and/or high frequencies of Killer-Rescue-associated elements. We show that for both strategies the most efficient approach for achieving spread of anti-pathogen genes within three years is generally to release adults of both sexes in multiple releases over time. Even though females in these releases should not transmit disease, there could be public concern over such releases, making the less efficient male-only release more practical. This study provides guidelines for operational approaches to population replacement genetic strategies, as well as illustrates the use of detailed spatial models to assist in safe and efficient implementation of such novel genetic strategies.     

Optimal surveillance and eradication of invasive species in heterogeneous landscapes

Cost-effective surveillance strategies are needed for efficient responses to biological invasions and must account for the trade-offs between surveillance effort and management costs. Less surveillance may allow greater population growth and spread prior to detection, thereby increasing the costs of damages and control. In addition, surveillance strategies are usually applied in environments under continual invasion pressure where the number, size and location of established populations are unknown prior to detection. We develop a novel modeling framework that accounts for these features of the decision and invasion environment and determines the long term sampling effort that minimises the total expected costs of new invasions. The optimal solution depends on population establishment and growth rates, sample sensitivity, and sample, eradication, and damage costs. We demonstrate how to optimise surveillance systems under budgetary constraints and find that accounting for spatial heterogeneity in sampling costs and establishment rates can greatly reduce management costs.     

Ontological model of multi-agent Smart-system for predicting drug properties based on modified algorithms of artificial immune systems

Estimating HIV incidence is crucial for monitoring the epidemiology of this infection, planning screening and intervention campaigns, and evaluating the effectiveness of control measures. However, owing to the long and variable period from HIV infection to the development of AIDS and the introduction of highly active antiretroviral therapy, accurate incidence estimation remains a major challenge. Numerous estimation methods have been proposed in epidemiological modeling studies, and here we review commonly-used methods for estimation of HIV incidence. We review the essential data required for estimation along with the advantages and disadvantages, mathematical structures and likelihood derivations of these methods. The methods include the classical back-calculation method, the method based on CD4+ T-cell depletion, the use of HIV case reporting data, the use of cohort study data, the use of serial or cross-sectional prevalence data, and biomarker approach. By outlining the mechanistic features of each method, we provide guidance for planning incidence estimation efforts, which may depend on national or regional factors as well as the availability of epidemiological or laboratory datasets.     

The role of screening and treatment in the transmission dynamics of HIV/AIDS and tuberculosis co-infection: a mathematical study

In this paper, we present a deterministic non-linear mathematical model for the transmission dynamics of HIV and TB co-infection and analyze it in the presence of screening and treatment. The equilibria of the model are computed and stability of these equilibria is discussed. The basic reproduction numbers corresponding to both HIV and TB are found and we show that the disease-free equilibrium is stable only when the basic reproduction numbers for both the diseases are less than one. When both the reproduction numbers are greater than one, the co-infection equilibrium point may exist. The co-infection equilibrium is found to be locally stable whenever it exists. The TB-only and HIV-only equilibria are locally asymptotically stable under some restriction on parameters. We present numerical simulation results to support the analytical findings. We observe that screening with proper counseling of HIV infectives results in a significant reduction of the number of individuals progressing to HIV. Additionally, the screening of TB reduces the infection prevalence of TB disease. The results reported in this paper clearly indicate that proper screening and counseling can check the spread of HIV and TB diseases and effective control strategies can be formulated around ‘screening with proper counseling’.     

Optimising and Communicating Options for the Control of Invasive Plant Disease When There Is Epidemiological Uncertainty

Although local eradication is routinely attempted following introduction of disease into a new region, failure is commonplace. Epidemiological principles governing the design of successful control are not well-understood. We analyse factors underlying the effectiveness of reactive eradication of localised outbreaks of invading plant disease, using citrus canker in Florida as a case study, although our results are largely generic, and apply to other plant pathogens (as we show via our second case study, citrus greening). We demonstrate how to optimise control via removal of hosts surrounding detected infection (i.e. localised culling) using a spatially-explicit, stochastic epidemiological model. We show how to define optimal culling strategies that take account of stochasticity in disease spread, and how the effectiveness of disease control depends on epidemiological parameters determining pathogen infectivity, symptom emergence and spread, the initial level of infection, and the logistics and implementation of detection and control. We also consider how optimal culling strategies are conditioned on the levels of risk acceptance/aversion of decision makers, and show how to extend the analyses to account for potential larger-scale impacts of a small-scale outbreak. Control of local outbreaks by culling can be very effective, particularly when started quickly, but the optimum strategy and its performance are strongly dependent on epidemiological parameters (particularly those controlling dispersal and the extent of any cryptic infection, i.e. infectious hosts prior to symptoms), the logistics of detection and control, and the level of local and global risk that is deemed to be acceptable. A version of the model we developed to illustrate our methodology and results to an audience of stakeholders, including policy makers, regulators and growers, is available online as an interactive, user-friendly interface at http://www.webidemics.com/. This version of our model allows the complex epidemiological principles that underlie our results to be communicated to a non-specialist audience.     

Optimal Treatment Strategies in the Context of ‘Treatment for Prevention’ against HIV-1 in Resource-Poor Settings

An estimated 2.7 million new HIV-1 infections occurred in 2010. `Treatment-for-prevention’ may strongly prevent HIV-1 transmission. The basic idea is that immediate treatment initiation rapidly decreases virus burden, which reduces the number of transmittable viruses and thereby the probability of infection. However, HIV inevitably develops drug resistance, which leads to virus rebound and nullifies the effect of `treatment-for-prevention’ for the time it remains unrecognized. While timely conducted treatment changes may avert periods of viral rebound, necessary treatment options and diagnostics may be lacking in resource-constrained settings. Within this work, we provide a mathematical platform for comparing different treatment paradigms that can be applied to many medical phenomena. We use this platform to optimize two distinct approaches for the treatment of HIV-1: (i) a diagnostic-guided treatment strategy, based on infrequent and patient-specific diagnostic schedules and (ii) a pro-active strategy that allows treatment adaptation prior to diagnostic ascertainment. Both strategies are compared to current clinical protocols (standard of care and the HPTN052 protocol) in terms of patient health, economic means and reduction in HIV-1 onward transmission exemplarily for South Africa. All therapeutic strategies are assessed using a coarse-grained stochastic model of within-host HIV dynamics and pseudo-codes for solving the respective optimal control problems are provided. Our mathematical model suggests that both optimal strategies (i)-(ii) perform better than the current clinical protocols and no treatment in terms of economic means, life prolongation and reduction of HIV-transmission. The optimal diagnostic-guided strategy suggests rare diagnostics and performs similar to the optimal pro-active strategy. Our results suggest that ‘treatment-for-prevention’ may be further improved using either of the two analyzed treatment paradigms.     

The role of pre-emptive culling in the control of foot-and-mouth disease

The 2001 foot-and-mouth disease epidemic was controlled by culling of infectious premises and pre-emptive culling intended to limit the spread of disease. Of the control strategies adopted, routine culling of farms that were contiguous to infected premises caused the most controversy. Here we perform a retrospective analysis of the culling of contiguous premises as performed in 2001 and a simulation study of the effects of this policy on reducing the number of farms affected by disease. Our simulation results support previous studies and show that a national policy of contiguous premises (CPs) culling leads to fewer farms losing livestock. The optimal national policy for controlling the 2001 epidemic is found to be the targeting of all contiguous premises, whereas for localized outbreaks in high animal density regions, more extensive fixed radius ring culling is optimal. Analysis of the 2001 data suggests that the lowest-risk CPs were generally prioritized for culling, however, even in this case, the policy is predicted to be effective. A sensitivity analysis and the development of a spatially heterogeneous policy show that the optimal culling level depends upon the basic reproductive ratio of the infection and the width of the dispersal kernel. These analyses highlight an important and probably quite general result: optimal control is highly dependent upon the distance over which the pathogen can be transmitted, the transmission rate of infection and local demography where the disease is introduced.     

Modeling the joint epidemics of TB and HIV in a South African township

We present a simple mathematical model with six compartments for the interaction between HIV and TB epidemics. Using data from a township near Cape Town, South Africa, where the prevalence of HIV is above 20% and where the TB notification rate is close to 2,000 per 100,000 per year, we estimate some of the model parameters and study how various control measures might change the course of these epidemics. Condom promotion, increased TB detection and TB preventive therapy have a clear positive effect. The impact of antiretroviral therapy on the incidence of HIV is unclear and depends on the extent to which it reduces sexual transmission. However, our analysis suggests that it will greatly reduce the TB notification rate.     

Optimal vaccine allocation for COVID-19 in the Netherlands: A data-driven prioritization

For the control of COVID-19, vaccination programmes provide a long-term solution. The amount of available vaccines is often limited, and thus it is crucial to determine the allocation strategy. While mathematical modelling approaches have been used to find an optimal distribution of vaccines, there is an excessively large number of possible allocation schemes to be simulated. Here, we propose an algorithm to find a near-optimal allocation scheme given an intervention objective such as minimization of new infections, hospitalizations, or deaths, where multiple vaccines are available. The proposed principle for allocating vaccines is to target subgroups with the largest reduction in the outcome of interest. We use an approximation method to reconstruct the age-specific transmission intensity (the next generation matrix), and express the expected impact of vaccinating each subgroup in terms of the observed incidence of infection and force of infection. The proposed approach is firstly evaluated with a simulated epidemic and then applied to the epidemiological data on COVID-19 in the Netherlands. Our results reveal how the optimal allocation depends on the objective of infection control. In the case of COVID-19, if we wish to minimize deaths, the optimal allocation strategy is not efficient for minimizing other outcomes, such as infections. In simulated epidemics, an allocation strategy optimized for an outcome outperforms other strategies such as the allocation from young to old, from old to young, and at random. Our simulations clarify that the current policy in the Netherlands (i.e., allocation from old to young) was concordant with the allocation scheme that minimizes deaths. The proposed method provides an optimal allocation scheme, given routine surveillance data that reflect ongoing transmissions. This approach to allocation is useful for providing plausible simulation scenarios for complex models, which give a more robust basis to determine intervention strategies.     

Comparison of Contact Patterns Relevant for Transmission of Respiratory Pathogens in Thailand and the Netherlands Using Respondent-Driven Sampling

Understanding infection dynamics of respiratory diseases requires the identification and quantification of behavioural, social and environmental factors that permit the transmission of these infections between humans. Little empirical information is available about contact patterns within real-world social networks, let alone on differences in these contact networks between populations that differ considerably on a socio-cultural level. Here we compared contact network data that were collected in the Netherlands and Thailand using a similar online respondent-driven method. By asking participants to recruit contact persons we studied network links relevant for the transmission of respiratory infections. We studied correlations between recruiter and recruited contacts to investigate mixing patterns in the observed social network components. In both countries, mixing patterns were assortative by demographic variables and random by total numbers of contacts. However, in Thailand participants reported overall more contacts which resulted in higher effective contact rates. Our findings provide new insights on numbers of contacts and mixing patterns in two different populations. These data could be used to improve parameterisation of mathematical models used to design control strategies. Although the spread of infections through populations depends on more factors, found similarities suggest that spread may be similar in the Netherlands and Thailand.     

Designing equitable antiretroviral allocation strategies in resource-constrained countries

Background

Recently, a global commitment has been made to expand access to antiretrovirals (ARVs) in the developing world. However, in many resource-constrained countries the number of individuals infected with HIV in need of treatment will far exceed the supply of ARVs, and only a limited number of health-care facilities (HCFs) will be available for ARV distribution. Deciding how to allocate the limited supply of ARVs among HCFs will be extremely difficult. Resource allocation decisions can be made on the basis of many epidemiological, ethical, or preferential treatment priority criteria.

Methods and Findings

Here we use operations research techniques, and we show how to determine the optimal strategy for allocating ARVs among HCFs in order to satisfy the equitable criterion that each individual infected with HIV has an equal chance of receiving ARVs. We present a novel spatial mathematical model that includes heterogeneity in treatment accessibility. We show how to use our theoretical framework, in conjunction with an equity objective function, to determine an optimal equitable allocation strategy (OEAS) for ARVs in resource-constrained regions. Our equity objective function enables us to apply the egalitarian principle of equity with respect to access to health care. We use data from the detailed ARV rollout plan designed by the government of South Africa to determine an OEAS for the province of KwaZulu–Natal. We determine the OEAS for KwaZulu–Natal, and we then compare this OEAS with two other ARV allocation strategies: (i) allocating ARVs only to Durban (the largest urban city in KwaZulu–Natal province) and (ii) allocating ARVs equally to all available HCFs. In addition, we compare the OEAS to the current allocation plan of the South African government (which is based upon allocating ARVs to 17 HCFs). We show that our OEAS significantly improves equity in treatment accessibility in comparison with these three ARV allocation strategies. We also quantify how the size of the catchment region surrounding each HCF, and the number of HCFs utilized for ARV distribution, alters the OEAS and the probability of achieving equity in treatment accessibility. We calculate that in order to achieve the greatest degree of treatment equity for individuals with HIV in KwaZulu–Natal, the ARVs should be allocated to 54 HCFs and each HCF should serve a catchment region of 40 to 60 km.

Conclusion

Our OEAS would substantially improve equality in treatment accessibility in comparison with other allocation strategies. Furthermore, our OEAS is extremely different from the currently planned strategy. We suggest that our novel methodology be used to design optimal ARV allocation strategies for resource-constrained countries.     

Mathematical models of HIV pathogenesis and treatment

We review mathematical models of HIV dynamics, disease progression, and therapy. We start by introducing a basic model of virus infection and demonstrate how it was used to study HIV dynamics and to measure crucial parameters that lead to a new understanding of the disease process. We discuss the diversity threshold model as an example of the general principle that virus evolution can drive disease progression and the destruction of the immune system. Finally, we show how mathematical models can be used to understand correlates of long-term immunological control of HIV, and to design therapy regimes that convert a progressing patient into a state of long-term non-progression.     

Modeling the impact of random screening and contact tracing in reducing the spread of HIV

Mathematical models can help predict the effectiveness of control measures on the spread of HIV and other sexually transmitted diseases (STDs) by reducing the uncertainty in assessing the impact of intervention strategies such as random screening and contact tracing. Even though contact tracing is one of the most effective methods used for controlling treatable STDs, it is still a controversial strategy for controlling HIV because of cost and confidentiality issues. To help estimate the effectiveness of these control measures, we formulate two models with random screening and contact tracing based on the differential infectivity (DI) model and the staged-progression (SP) model. We derive formulas for the reproductive numbers and the endemic equilibria and compare the impact that random screening and contact tracing have in slowing the epidemic in the two models. In the DI model the infected population is divided into groups according to their infectiousness, and HIV is largely spread by a small, highly infectious, group of superspreaders. In this model contact tracing is an effective approach to identifying the superspreaders and has a large effect in slowing the epidemic. In the SP model every infected individual goes through a series of infection stages and the virus is primarily spread by individuals in an initial highly infectious stage or in the late stages of the disease. In this model random screening is more effective than for the DI model, and contact tracing is less effective. Thus the effectiveness of the intervention strategy strongly depends on the underlying etiology of the disease transmission.     

An intuitive formulation for the reproductive number for the spread of diseases in heterogeneous populations

The thresholds for mathematical epidemiology models specify the critical conditions for an epidemic to grow or die out. The reproductive number can provide significant insight into the transmission dynamics of a disease and can guide strategies to control its spread. We define the mean number of contacts, the mean duration of infection, and the mean transmission probability appropriately for certain epidemiological models, and construct a simplified formulation of the reproductive number as the product of these quantities. When the spread of the epidemic depends strongly upon the heterogeneity of the populations, the epidemiological models must account for this heterogeneity, and the expressions for the reproductive number become correspondingly more complex. We formulate several models with different heterogeneous structures and demonstrate how to define the mean quantities for an explicit expression for the reproductive number. In complex heterogeneous models, it seems necessary to define the reproductive number for each structured subgroup or cohort and then use the average of these reproductive numbers weighted by their heterogeneity to estimate the reproductive number for the total population.     

Approaches to antiretroviral therapy in China

INTRODUCTION Advances in antiretroviral therapy (ART) have resulted in dramatic declines in hospitalizations and death rates in the U.S. and Europe, with concomitant steady increases in life expectancy. In China, the government has recog- nized the imminent threat that HIV/AIDS poses to its popu- lation and responded with a national antiretroviral (ARV) treatment program providing ARV drugs free to those most in need [1]. However, given the potentially catastrophic impact of HIV …     

Mathematical model for optimal control in wastewater discharges: the global performance

In this work we show how mathematical models and optimal control techniques can help us to solve some problems of environmental engineering, more precisely, water pollution problems arising from wastewater discharges into coastal areas or rivers. We deal with a complete two-dimensional mathematical model for the evolution of pollutant concentration in a shallow water domain. By integrating this model we obtain a zero-dimensional one and we use it to study the global performance of the system in a realistic situation. Finally, by using the two-dimensional model, we recall two optimal control problems related to the wastewater disposal problem.