1
|
d'Onofrio A, Iannelli M, Manfredi P, Marinoschi G. Epidemic control by social distancing and vaccination: Optimal strategies and remarks on the COVID-19 Italian response policy. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2024; 21:6493-6520. [PMID: 39176405 DOI: 10.3934/mbe.2024283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/24/2024]
Abstract
After the many failures in the control of the COVID-19 pandemic, identifying robust principles of epidemic control will be key in future preparedness. In this work, we propose an optimal control model of an age-of-infection transmission model under a two-phase control regime where social distancing is the only available control tool in the first phase, while the second phase also benefits from the arrival of vaccines. We analyzed the problem by an ad-hoc numerical algorithm under a strong hypothesis implying a high degree of prioritization to the protection of health from the epidemic attack, which we termed the "low attack rate" hypothesis. The outputs of the model were also compared with the data from the Italian COVID-19 experience to provide a crude assessment of the goodness of the enacted interventions prior to the onset of the Omicron variant.
Collapse
Affiliation(s)
- Alberto d'Onofrio
- Dipartimento di Matematica, Informatica e Geoscienze, Università di Trieste, Via Alfonso Valerio 12, Edificio H2bis, 34127 Trieste, Italy
| | - Mimmo Iannelli
- Department of Mathematics, University of Trento, Via Sommarive 14, 38123 Trento, Italy
| | - Piero Manfredi
- Dipartimento di Economia e Management, University of Pisa, Via Ridolfi 10, 56124 Pisa, Italy
| | - Gabriela Marinoschi
- Gheorghe Mihoc-Caius Iacob Institute of Mathematical Statistics and Applied Mathematics, Romanian Academy, Bucharest, Romania
| |
Collapse
|
2
|
Bolton KJ, McCaw JM, Dafilis MP, McVernon J, Heffernan JM. Seasonality as a driver of pH1N12009 influenza vaccination campaign impact. Epidemics 2023; 45:100730. [PMID: 38056164 DOI: 10.1016/j.epidem.2023.100730] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2023] [Revised: 07/18/2023] [Accepted: 11/16/2023] [Indexed: 12/08/2023] Open
Abstract
Although the most recent respiratory virus pandemic was triggered by a Coronavirus, sustained and elevated prevalence of highly pathogenic avian influenza viruses able to infect mammalian hosts highlight the continued threat of pandemics of influenza A virus (IAV) to global health. Retrospective analysis of pandemic outcomes, including comparative investigation of intervention efficacy in different regions, provide important contributions to the evidence base for future pandemic planning. The swine-origin IAV pandemic of 2009 exhibited regional variation in onset, infection dynamics and annual infection attack rates (IARs). For example, the UK experienced three severe peaks of infection over two influenza seasons, whilst Australia experienced a single severe wave. We adopt a seasonally forced 2-subtype model for the transmission of pH1N12009 and seasonal H3N2 to examine the role vaccination campaigns may play in explaining differences in pandemic trajectories in temperate regions. Our model differentiates between the nature of vaccine- and infection-acquired immunity. In particular, we assume that immunity triggered by infection elicits heterologous cross-protection against viral shedding in addition to long-lasting neutralising antibody, whereas vaccination induces imperfect reduction in susceptibility. We employ an Approximate Bayesian Computation (ABC) framework to calibrate the model using data for pH1N12009 seroprevalence, relative subtype dominance, and annual IARs for Australia and the UK. Heterologous cross-protection substantially suppressed the pandemic IAR over the posterior, with the strength of protection against onward transmission inversely correlated with the initial reproduction number. We show that IAV pandemic timing relative to the usual seasonal influenza cycle influenced the size of the initial waves of pH1N12009 in temperate regions and the impact of vaccination campaigns.
Collapse
Affiliation(s)
- Kirsty J Bolton
- School of Mathematical Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK.
| | - James M McCaw
- School of Mathematics and Statistics, The University of Melbourne, Parkville, Australia; Centre for Epidemiology and Biostatistics, Melbourne School of Population and Global Health, The University of Melbourne, Parkville, Australia
| | - Mathew P Dafilis
- Centre for Epidemiology and Biostatistics, Melbourne School of Population and Global Health, The University of Melbourne, Parkville, Australia
| | - Jodie McVernon
- Peter Doherty Institute for Infection and Immunity, The Royal Melbourne Hospital and The University of Melbourne, Parkville, Australia
| | - Jane M Heffernan
- Centre for Disease Modelling, Mathematics & Statistics, York University, Canada
| |
Collapse
|
3
|
Rysava K, Tildesley MJ. Identification of dynamical changes of rabies transmission under quarantine: Community-based measures towards rabies elimination. PLoS Comput Biol 2023; 19:e1011187. [PMID: 38100528 PMCID: PMC10756519 DOI: 10.1371/journal.pcbi.1011187] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Revised: 12/29/2023] [Accepted: 11/13/2023] [Indexed: 12/17/2023] Open
Abstract
Quarantine has been long used as a public health response to emerging infectious diseases, particularly at the onset of an epidemic when the infected proportion of a population remains identifiable and logistically tractable. In theory, the same logic should apply to low-incidence infections; however, the application and impact of quarantine in low prevalence settings appears less common and lacks a formal analysis. Here, we present a quantitative framework using a series of progressively more biologically realistic models of canine rabies in domestic dogs and from dogs to humans, a suitable example system to characterize dynamical changes under varying levels of dog quarantine. We explicitly incorporate health-seeking behaviour data to inform the modelling of contact-tracing and exclusion of rabies suspect and probable dogs that can be identified through bite-histories of patients presenting at anti-rabies clinics. We find that a temporary quarantine of rabies suspect and probable dogs provides a powerful tool to curtail rabies transmission, especially in settings where optimal vaccination coverage is yet to be achieved, providing a critical stopgap to reduce the number of human and animal deaths due to rabid bites. We conclude that whilst comprehensive measures including sensitive surveillance and large-scale vaccination of dogs will be required to achieve disease elimination and sustained freedom given the persistent risk of rabies re-introductions, quarantine offers a low-cost community driven solution to intersectoral health burden.
Collapse
Affiliation(s)
- Kristyna Rysava
- The Zeeman Institute for Systems Biology & Infectious Disease Epidemiology Research, School of Life Sciences and Mathematics Institute, University of Warwick, Coventry, United Kingdom
| | - Michael J. Tildesley
- The Zeeman Institute for Systems Biology & Infectious Disease Epidemiology Research, School of Life Sciences and Mathematics Institute, University of Warwick, Coventry, United Kingdom
| |
Collapse
|
4
|
Abstract
Recent Covid-19 pandemic has demonstrated the need of efficient epidemic outbreak management. We study the optimal control problem of minimizing the fraction of infected population by applying vaccination and treatment control strategies, while at the same time minimizing the cost of implementing them. We model the epidemic using the degree based Susceptible–Infected–Recovered (SIR) compartmental model. We study the impact of varying network topologies on the optimal epidemic management strategies and present results for the Erdős–Rényi, scale free, and real world networks. For efficient computational modeling we form groups of groups of degree classes, and apply separate vaccination and treatment control signals to each group. This allows us to identify the degree classes that play a significant role in mitigating the epidemic for a given network topology. We compare the optimal control strategy with non optimal strategies (constant control and no control) and study the effect of various model parameters on the system. We identify which strategy (vaccination/treatment) plays a significant role in controlling the epidemic on different network topologies. We also study the effect of the cost of vaccination and treatment controls on the resource allocation. We find that the optimal strategy achieves significant improvements over the non optimal heuristics for all networks studied in this paper. Our results may be of interest to governments and healthcare authorities for devising effective vaccination and treatment campaigns during an epidemic outbreak.
Collapse
|
5
|
Rigatos G, Abbaszadeh M, Cuccurullo G. A nonlinear optimal control method against the spreading of epidemics. INT J BIOMATH 2022. [DOI: 10.1142/s1793524522500267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
To define a vaccination policy and antiviral treatment against the spreading of viral infections a nonlinear optimal (H-infinity) control approach is proposed. Actually, because of the scarcity of the resources for treating infectious diseases in terms of vaccines, antiviral drugs and other medical facilities, there is need to implement optimal control against the epidemics deployment. In this approach, the state-space model of the epidemics dynamics undergoes first approximate linearization around a temporary operating point which is recomputed at each time-step of the control method. The linearization is based on Taylor series expansion and on the computation of the associated Jacobian matrices. Next, an optimal (H-infinity) feedback controller is developed for the approximately linearized model of the epidemics. To compute the controller’s feedback gains an algebraic Riccati equation is solved at each iteration of the control algorithm. Furthermore, the global asymptotic stability properties of the control scheme are proven through Lyapunov stability analysis. This paper’s results confirm that optimal control of the infectious disease dynamics allows for eliminating its spreading while also keeping moderate the consumption of the related medication, that is vaccines and antiviral drugs.
Collapse
Affiliation(s)
- G. Rigatos
- Unit of Industrial Automation, Industrial Systems Institute, 26504, Rion Patras Greece, Greece
| | - M. Abbaszadeh
- Department of ECSE, Rensselaer Polytechnic Institute 12065, NY, USA
| | - G. Cuccurullo
- Department of Industrial Engineering, University of Salerno, Fisciano, 84084, Italy
| |
Collapse
|
6
|
Chang L, Gao S, Wang Z. Optimal control of pattern formations for an SIR reaction-diffusion epidemic model. J Theor Biol 2022; 536:111003. [PMID: 35026213 DOI: 10.1016/j.jtbi.2022.111003] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 12/29/2021] [Accepted: 12/31/2021] [Indexed: 11/19/2022]
Abstract
Patterns arising from the reaction-diffusion epidemic model provide insightful aspects into the transmission of infectious diseases. For a classic SIR reaction-diffusion epidemic model, we review its Turing pattern formations with different transmission rates. A quantitative indicator, "normal serious prevalent area (NSPA)", is introduced to characterize the relationship between patterns and the extent of the epidemic. The extent of epidemic is positively correlated to NSPA. To effectively reduce NSPA of patterns under the large transmission rates, taken removed (recovery or isolation) rate as a control parameter, we consider the mathematical formulation and numerical solution of an optimal control problem for the SIR reaction-diffusion model. Numerical experiments demonstrate the effectiveness of our method in terms of control effect, control precision and control cost.
Collapse
Affiliation(s)
- Lili Chang
- Complex Systems Research Center, Shanxi University, Taiyuan 030006, Shanxi, China.
| | - Shupeng Gao
- School of Mechanical Engineering Northwestern Polytechnical University, Xi'an 710072, Shaanxi, China; School of Artificial Intelligence, Optics, and Electronics (iOPEN), Northwestern Polytechnical University, Xi'an 710072, Shaanxi, China
| | - Zhen Wang
- School of Mechanical Engineering Northwestern Polytechnical University, Xi'an 710072, Shaanxi, China; School of Artificial Intelligence, Optics, and Electronics (iOPEN), Northwestern Polytechnical University, Xi'an 710072, Shaanxi, China.
| |
Collapse
|
7
|
Dynamics of Trachoma Epidemic in Human Contact Network with Seasonally Varying Infectious Medium. PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES INDIA SECTION A-PHYSICAL SCIENCES 2021. [DOI: 10.1007/s40010-021-00754-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
8
|
Saddique A, Adnan S, Bokhari H, Azam A, Rana MS, Khan MM, Hanif M, Sharif S. Prevalence and Associated Risk Factor of COVID-19 and Impacts of Meteorological and Social Variables on Its Propagation in Punjab, Pakistan. EARTH SYSTEMS AND ENVIRONMENT 2021; 5:785-798. [PMID: 34723081 PMCID: PMC8260326 DOI: 10.1007/s41748-021-00218-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Accepted: 04/10/2021] [Indexed: 06/13/2023]
Abstract
The current study identifies the spatial distribution of COVID-19 cases and its association with meteorological and social variables in Punjab (densely populated province of Pakistan). To identify the COVID-19 propagation, the weekly growth, recovery, and deaths rate have also been calculated. The geographic information system (GIS) has used to determine COVID-19 impacts on gender (male/female), age groups, and causalities over an affected population (km-2) for the period of 11th March to 12th August, 2020 in each district of province. Our results show that 43 peak days (where daily positive cases were above 900) have been observed in Punjab during 27th May to 8th July, 2020. The high population density districts, i.e., Lahore and Islamabad, have been affected (five persons per square kilometers) due to COVID-19, whereas the maximum death tolls (> 50 persons per millions) have also been observed in these urban districts. The meteorological variables (temperature, humidity, heat index, and ultraviolet index) show negative significant relationship to basic reproduction number (R0), whereas daily COVID-19 cases are positively correlated to aerosols concentration at 95% confidence level. The government intervention (stringency index) shows a positive impact to reduce the COVID-19 cases over the province. Keeping in view the COVID-19 behavior and climatology of the region, it has been identified that the COVID-19 cases may likely to increase during the dry period (high concentration of aerosols) i.e., October-December, 2020 and post-spring season (April to June), 2021 in urban areas of Pakistan. This study provides an overview on districts vulnerability that would help the policy makers, health agencies to plan their activities to reduce the COVID-19 impacts.
Collapse
Affiliation(s)
- Arbab Saddique
- COMSATS University Islamabad/Kohsar University, Islamabad/Murree, Pakistan
| | - Shahzada Adnan
- Pakistan Meteorological Department, Sector H-8/2, Islamabad, Pakistan
| | - Habib Bokhari
- COMSATS University Islamabad/Kohsar University, Islamabad/Murree, Pakistan
| | - Asima Azam
- Shaheed Benazir Bhutto Women University, Peshawar, Pakistan
| | | | | | - Muhammad Hanif
- Pakistan Meteorological Department, Sector H-8/2, Islamabad, Pakistan
| | - Shawana Sharif
- Shaheed Benazir Bhutto Hospital, Rawalpindi Medical University, Rawalpindi, Pakistan
| |
Collapse
|
9
|
Adnan S, Hanif M, Khan AH, Latif M, Ullah K, Bashir F, Kamil S, Haider S. Impact of Heat Index and Ultraviolet Index on COVID-19 in Major Cities of Pakistan. J Occup Environ Med 2021; 63:98-103. [PMID: 33021515 PMCID: PMC7864608 DOI: 10.1097/jom.0000000000002039] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
INTRODUCTION The world population is under the grip of global pandemic of COVID-19. The present study analyzed relationship between meteorological parameters and COVID-19 in three major cities of Pakistan, that is, Karachi, Lahore, and Peshawar. METHODS The impacts of heat index (HI) and ultraviolet index (UVI) over daily COVID-19 cases have examined to identify its transmission and propagation. The significance of basic reproductive number (R0), growth rate (Gr) and doubling time (Td) of COVID-19 with HI and UVI was determined. RESULTS Both indices show a significant positive correlation (at 5% significance level) to R0, Td, and Gr of COVID-19 patients. Our results showed that the minimum threshold temperature of 33 °C for HI (with a positive variation of 3 °C to 5 °C) put a significant impact on new cases. CONCLUSION HI and UVI impacted significantly to decline COVID-19 cases over the region.
Collapse
Affiliation(s)
- Shahzada Adnan
- Pakistan Meteorological Department (Dr Adnan, Dr Hanif, Dr Khan, Dr Bashir, Dr Kamil, Dr Haider); Department of Meteorology, COMSATS University Islamabad, Chak Shazad (Dr Latif, Dr Ullah), Islamabad, Pakistan
| | | | | | | | | | | | | | | |
Collapse
|
10
|
Iqbal MM, Abid I, Hussain S, Shahzad N, Waqas MS, Iqbal MJ. The effects of regional climatic condition on the spread of COVID-19 at global scale. THE SCIENCE OF THE TOTAL ENVIRONMENT 2020; 739:140101. [PMID: 32531684 PMCID: PMC7280824 DOI: 10.1016/j.scitotenv.2020.140101] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/21/2020] [Revised: 06/07/2020] [Accepted: 06/08/2020] [Indexed: 05/17/2023]
Abstract
The pandemic outbreak of the novel coronavirus epidemic disease (COVID-19) is spreading like a diffusion-reaction in the world and almost 208 countries and territories are being affected around the globe. It became a sever health and socio-economic problem, while the world has no vaccine to combat this virus. This research aims to analyze the connection between the fast spread of COVID-19 and regional climate parameters over a global scale. In this research, we collected the data of COVID-19 cases from the time of 1st reported case to the 5th June 2020 in different affected countries and regional climatic parameters data from January 2020 to 5th June 2020. It was found that most of the countries located in the relatively lower temperature region show a rapid increase in the COVID-19 cases than the countries locating in the warmer climatic regions despite their better socio-economic conditions. A correlation between metrological parameters and COVID-19 cases was observed. Average daylight hours are correlated to total the COVID-19 cases with a coefficient of determination of 0.42, while average high-temperature shows a correlation of 0.59 and 0.42 with total COVID-19 cases and death cases respectively. The finding of the study will help international health organizations and local administrations to combat and well manage the spread of COVID-19.
Collapse
Affiliation(s)
| | - Irfan Abid
- National University of Sciences and Technology (NUST), Islamabad, Pakistan
| | - Saddam Hussain
- Department of Irrigation and Drainage, University of Agriculture, Faisalabad, Pakistan
| | - Naeem Shahzad
- National University of Sciences and Technology (NUST), Islamabad, Pakistan
| | - Muhammad Sohail Waqas
- Soil Conservation Group, Agriculture Department (Field Wing), Government of the Punjab, Pakistan
| | | |
Collapse
|
11
|
Sharp JA, Browning AP, Mapder T, Baker CM, Burrage K, Simpson MJ. Designing combination therapies using multiple optimal controls. J Theor Biol 2020; 497:110277. [PMID: 32294472 DOI: 10.1016/j.jtbi.2020.110277] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2019] [Revised: 02/21/2020] [Accepted: 04/06/2020] [Indexed: 01/31/2023]
Abstract
Strategic management of populations of interacting biological species routinely requires interventions combining multiple treatments or therapies. This is important in key research areas such as ecology, epidemiology, wound healing and oncology. Despite the well developed theory and techniques for determining single optimal controls, there is limited practical guidance supporting implementation of combination therapies. In this work we use optimal control theory to calculate optimal strategies for applying combination therapies to a model of acute myeloid leukaemia. We present a versatile framework to systematically explore the trade-offs that arise in designing combination therapy protocols using optimal control. We consider various combinations of continuous and bang-bang (discrete) controls, and we investigate how the control dynamics interact and respond to changes in the weighting and form of the pay-off characterising optimality. We demonstrate that the optimal controls respond non-linearly to treatment strength and control parameters, due to the interactions between species. We discuss challenges in appropriately characterising optimality in a multiple control setting and provide practical guidance for applying multiple optimal controls. Code used in this work to implement multiple optimal controls is available on GitHub.
Collapse
Affiliation(s)
- Jesse A Sharp
- School of Mathematical Sciences, Queensland University of Technology (QUT), Australia; ARC Centre of Excellence for Mathematical and Statistical Frontiers, QUT, Australia.
| | - Alexander P Browning
- School of Mathematical Sciences, Queensland University of Technology (QUT), Australia; ARC Centre of Excellence for Mathematical and Statistical Frontiers, QUT, Australia
| | - Tarunendu Mapder
- School of Mathematical Sciences, Queensland University of Technology (QUT), Australia; ARC Centre of Excellence for Mathematical and Statistical Frontiers, QUT, Australia
| | - Christopher M Baker
- School of Mathematical Sciences, Queensland University of Technology (QUT), Australia; ARC Centre of Excellence for Mathematical and Statistical Frontiers, QUT, Australia; School of Mathematics and Statistics, The University of Melbourne, Australia
| | - Kevin Burrage
- School of Mathematical Sciences, Queensland University of Technology (QUT), Australia; ARC Centre of Excellence for Mathematical and Statistical Frontiers, QUT, Australia; Department of Computer Science, University of Oxford, UK (Visiting Professor)
| | - Matthew J Simpson
- School of Mathematical Sciences, Queensland University of Technology (QUT), Australia
| |
Collapse
|
12
|
Optimal control of acute myeloid leukaemia. J Theor Biol 2019; 470:30-42. [PMID: 30853393 DOI: 10.1016/j.jtbi.2019.03.006] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2018] [Revised: 03/06/2019] [Accepted: 03/07/2019] [Indexed: 12/14/2022]
Abstract
Acute myeloid leukaemia (AML) is a blood cancer affecting haematopoietic stem cells. AML is routinely treated with chemotherapy, and so it is of great interest to develop optimal chemotherapy treatment strategies. In this work, we incorporate an immune response into a stem cell model of AML, since we find that previous models lacking an immune response are inappropriate for deriving optimal control strategies. Using optimal control theory, we produce continuous controls and bang-bang controls, corresponding to a range of objectives and parameter choices. Through example calculations, we provide a practical approach to applying optimal control using Pontryagin's Maximum Principle. In particular, we describe and explore factors that have a profound influence on numerical convergence. We find that the convergence behaviour is sensitive to the method of control updating, the nature of the control, and to the relative weighting of terms in the objective function. All codes we use to implement optimal control are made available.
Collapse
|