1
|
Pant B, Gumel AB. Mathematical assessment of the roles of age heterogeneity and vaccination on the dynamics and control of SARS-CoV-2. Infect Dis Model 2024; 9:828-874. [PMID: 38725431 PMCID: PMC11079469 DOI: 10.1016/j.idm.2024.04.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 04/10/2024] [Accepted: 04/11/2024] [Indexed: 05/12/2024] Open
Abstract
The COVID-19 pandemic, caused by SARS-CoV-2, disproportionately affected certain segments of society, particularly the elderly population (which suffered the brunt of the burden of the pandemic in terms of severity of the disease, hospitalization, and death). This study presents a generalized multigroup model, with m heterogeneous sub-populations, to assess the population-level impact of age heterogeneity and vaccination on the transmission dynamics and control of the SARS-CoV-2 pandemic in the United States. Rigorous analysis of the model for the homogeneous case (i.e., the model with m = 1) reveal that its disease-free equilibrium is globally-asymptotically stable for two special cases (with perfect vaccine efficacy or negligible disease-induced mortality) whenever the associated reproduction number is less than one. The model has a unique and globally-asymptotically stable endemic equilibrium, for special a case, when the associated reproduction threshold exceeds one. The homogeneous model was fitted using the observed cumulative mortality data for the United States during three distinct waves (Waves A (October 17, 2020 to April 5, 2021), B (July 9, 2021 to November 7, 2021) and C (January 1, 2022 to May 7, 2022)) chosen to align with time periods when the Alpha, Delta and Omicron were, respectively, the predominant variants in the United States. The calibrated model was used to derive a theoretical expression for achieving vaccine-derived herd immunity (needed to eliminate the disease in the United States). It was shown that, using the one-group homogeneous model, vaccine-derived herd immunity is not attainable during Wave C of the pandemic in the United States, regardless of the coverage level of the fully-vaccinated individuals. Global sensitivity analysis was carried out to determine the parameters of the model that have the most influence on the disease dynamics and burden. These analyses reveal that control and mitigation strategies that may be very effective during one wave may not be so very effective during the other wave or waves. However, strategies that target asymptomatic and pre-symptomatic infectious individuals are shown to be consistently effective across all waves. To study the impact of the disproportionate effect of COVID-19 on the elderly population, we considered the heterogeneous model for the case where the total population is subdivided into the sub-populations of individuals under 65 years of age and those that are 65 and older. The resulting two-group heterogeneous model, which was also fitted using the cumulative mortality data for wave C, was also rigorously analysed. Unlike for the case of the one-group model, it was shown, for the two-group model, that vaccine-derived herd immunity can indeed be achieved during Wave C of the pandemic if at least 61% of the populace is fully vaccinated. Thus, this study shows that adding age heterogeneity into a SARS-CoV-2 vaccination model with homogeneous mixing significantly reduces the level of vaccination coverage needed to achieve vaccine-derived herd immunity (specifically, for the heterogeneous model, herd-immunity can be attained during Wave C if a moderate proportion of susceptible individuals are fully vaccinated). The consequence of this result is that vaccination models for SARS-CoV-2 that do not explicitly account for age heterogeneity may be overestimating the level of vaccine-derived herd immunity threshold needed to eliminate the SARS-CoV-2 pandemic.
Collapse
Affiliation(s)
- Binod Pant
- Department of Mathematics, University of Maryland, College Park, MD, 20742, USA
| | - Abba B. Gumel
- Department of Mathematics, University of Maryland, College Park, MD, 20742, USA
- Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria, 0002, South Africa
| |
Collapse
|
2
|
Aldila D, Fardian BL, Chukwu CW, Hifzhudin Noor Aziz M, Kamalia PZ. Improving tuberculosis control: assessing the value of medical masks and case detection-a multi-country study with cost-effectiveness analysis. ROYAL SOCIETY OPEN SCIENCE 2024; 11:231715. [PMID: 39100187 PMCID: PMC11296204 DOI: 10.1098/rsos.231715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 02/08/2024] [Accepted: 03/18/2024] [Indexed: 08/06/2024]
Abstract
Tuberculosis (TB) remains a significant global health concern, necessitating effective control strategies. This article presents a mathematical model to evaluate the comparative effectiveness of medical mask usage and case detection in TB control. The model is constructed as a system of ordinary differential equations and incorporates crucial aspects of TB dynamics, including slow-fast progression, medical mask use, case detection, treatment interventions and differentiation between symptomatic and asymptomatic cases. A key objective of TB control is to ensure that the reproduction number,R c , remains below unity to achieve TB elimination or persistence ifR c exceeds 1. Our mathematical analysis reveals the presence of a transcritical bifurcation when theR c = 1 signifies a critical juncture in TB control strategies. These results confirm that the effectiveness of case detection in diminishing the endemic population of symptomatic individuals within a TB-endemic equilibrium depends on exceeding a critical threshold value. Furthermore, our model is calibrated using TB yearly case incidence data per 100 000 population from Indonesia, India, Lesotho and Angola. We employed the bootstrap resampling residual approach to assess the uncertainty inherent in our parameter estimates which provides a comprehensive distribution of the parameter values. Despite a declining trend in new incidence, these four countries exhibit a reproduction number greater than 1, indicating persistent TB cases in the presence of ongoing TB control programmes. We employ the partial rank correlation coefficient in conjunction with the Latin hypercube sampling method to conduct a global sensitivity analysis of theR c parameter for each fitted parameter in every country. We find that the medical mask use is more sensitive to reduceR c compared with the case detection implementation. To further gain insight into the necessary control strategy, we formulated an optimal control and studied the cost-effectiveness analysis of our model to investigate the impact of case detection and medical mask use as control measures in TB spread. Cost-effectiveness analysis demonstrates that combining these interventions emerges as the most cost-effective strategy for TB control. Our findings highlight the critical importance of medical masks and their efficacy coupled with case detection in shaping TB control dynamics, elucidating the primary parameter of concern for managing the control reproduction number. We envisage our findings to have implications and be vital for TB control if implemented by policymakers and healthcare practitioners involved in TB control efforts.
Collapse
Affiliation(s)
- Dipo Aldila
- Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Indonesia, Depok16424, Indonesia
| | - Basyar Lauzha Fardian
- Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Indonesia, Depok16424, Indonesia
| | | | | | - Putri Zahra Kamalia
- Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Indonesia, Depok16424, Indonesia
| |
Collapse
|
3
|
Lamba S, Das T, Srivastava PK. Impact of infectious density-induced additional screening and treatment saturation on COVID-19: Modeling and cost-effective optimal control. Infect Dis Model 2024; 9:569-600. [PMID: 38558959 PMCID: PMC10978547 DOI: 10.1016/j.idm.2024.03.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 02/18/2024] [Accepted: 03/11/2024] [Indexed: 04/04/2024] Open
Abstract
This study introduces a novel SI2HR model, where "I2" denotes two infectious classes representing asymptomatic and symptomatic infections, aiming to investigate and analyze the cost-effective optimal control measures for managing COVID-19. The model incorporates a novel concept of infectious density-induced additional screening (IDIAS) and accounts for treatment saturation. Furthermore, the model considers the possibility of reinfection and the loss of immunity in individuals who have previously recovered. To validate and calibrate the proposed model, real data from November-December 2022 in Hong Kong are utilized. The estimated parameters obtained from this calibration process are valuable for prediction purposes and facilitate further numerical simulations. An analysis of the model reveals that delays in screening, treatment, and quarantine contribute to an increase in the basic reproduction number R0, indicating a tendency towards endemicity. In particular, from the elasticity of R0, we deduce that normalized sensitivity indices of baseline screening rate (θ), quarantine rates (γ, αs), and treatment rate (α) are negative, which shows that delaying any of these may cause huge surge in R0, ultimately increases the disease burden. Further, by the contour plots, we note the two-parameter behavior of the infectives (both symptomatic and asymptomatic). Expanding upon the model analysis, an optimal control problem (OCP) is formulated, incorporating three control measures: precautionary interventions, boosted IDIAS, and boosted treatment. The Pontryagin's maximum principle and the forward-backward sweep method are employed to solve the OCP. The numerical simulations highlight that enhanced screening and treatment, coupled with preventive interventions, can effectively contribute to sustainable disease control. However, the cost-effectiveness analysis (CEA) conducted in this study suggests that boosting IDIAS alone is the most economically efficient and cost-effective approach compared to other strategies. The CEA results provide valuable insights into identifying specific strategies based on their cost-efficacy ranking, which can be implemented to maximize impact while minimizing costs. Overall, this research offers significant insights for policymakers and healthcare professionals, providing a framework to optimize control efforts for COVID-19 or similar epidemics in the future.
Collapse
Affiliation(s)
- Sonu Lamba
- Department of Mathematics, Indian Institute of Technology Patna Bihta – 801106, Patna, Bihar, India
| | - Tanuja Das
- Department of Mathematics, Indian Institute of Technology Patna Bihta – 801106, Patna, Bihar, India
- Department of Mathematics and Statistics, University of New Brunswick Fredericton, NB, E3B 5A3, Canada
| | - Prashant K. Srivastava
- Department of Mathematics, Indian Institute of Technology Patna Bihta – 801106, Patna, Bihar, India
| |
Collapse
|
4
|
Tollett Q, Safdar S, Gumel AB. Dynamics of a two-group model for assessing the impacts of pre-exposure prophylaxis, testing and risk behaviour change on the spread and control of HIV/AIDS in an MSM population. Infect Dis Model 2024; 9:103-127. [PMID: 38187461 PMCID: PMC10770619 DOI: 10.1016/j.idm.2023.11.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 11/09/2023] [Accepted: 11/10/2023] [Indexed: 01/09/2024] Open
Abstract
Although much progress has been made in reducing the public health burden of the human immunodeficiency virus (HIV), which causes acquired immunodeficiency syndrome (AIDS), since its emergence in the 1980s (largely due to the large-scale use and availability of potent antiviral therapy, improved diagnostic and intervention and mitigation measures), HIV remains an important public health challenge globally, including in the United States. This study is based on the use of mathematical modeling approaches to assess the population-level impact of pre-exposure prophylaxis (PrEP), voluntary testing (to detect undetected HIV-infected individuals), and changes in human behavior (with respect to risk structure), on the spread and control of HIV/AIDS in an MSM (men-who-have sex-with-men) population. Specifically, a novel two-group mathematical model, which stratifies the total MSM population based on risk (low or high) of acquisition of HIV infection, is formulated. The model undergoes a PrEP-induced backward bifurcation when the control reproduction number of the model is less than one if the efficacy of PrEP to prevent a high-risk susceptible MSM individual from acquiring HIV infection is not perfect (the consequence of which is that, while necessary, having the reproduction number of the model less than one is no longer sufficient for the elimination of the disease in the MSM population). For the case where the efficacy of PrEP is perfect, this study shows that the disease-free equilibrium of the two-group model is globally-asymptotically stable when the associated control reproduction number of the model is less than one. Global sensitivity analysis was carried out to identify the main parameters of the model that have the highest influence on the value of the control reproduction number of the model (thereby, having the highest influence on the disease burden in the MSM population). Numerical simulations of the model, using a plausible range of parameter values, show that if half of the MSM population considered adhere strictly to the specified PrEP regimen (while other interventions are maintained at their baseline values), a reduction of about 22% of the new yearly HIV cases recorded at the peak of the disease could be averted (compared to the worst-case scenario where PrEP-based intervention is not implemented in the MSM population). The yearly reduction at the peak increases to about 50% if the PrEP coverage in the MSM population increases to 80%. This study showed, based on the parameter values used in the simulations, that the prospects of elimination of HIV/AIDS in the MSM community are promising if high-risk susceptible individuals are no more than 15% more likely to acquire HIV infection, in comparison to their low-risk counterparts. Furthermore, these prospects are significantly improved if undetected HIV-infected individuals are detected within an optimal period of time.
Collapse
Affiliation(s)
- Queen Tollett
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ, 85287, USA
| | - Salman Safdar
- Department of Mathematics, University of Karachi, University Road, 75270, Pakistan
| | - Abba B. Gumel
- Department of Mathematics, University of Maryland, College Park, MD, 20742, USA
- Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria, 0002, South Africa
| |
Collapse
|
5
|
Aldila D, Awdinda N, Fatmawati, Herdicho FF, Ndii MZ, Chukwu CW. Optimal control of pneumonia transmission model with seasonal factor: Learning from Jakarta incidence data. Heliyon 2023; 9:e18096. [PMID: 37519720 PMCID: PMC10375561 DOI: 10.1016/j.heliyon.2023.e18096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Revised: 07/05/2023] [Accepted: 07/06/2023] [Indexed: 08/01/2023] Open
Abstract
Pneumonia is a dangerous disease that can lead to death without proper treatment. It is caused by a bacterial infection that leads to the inflammation of the air sacs in human lungs and potentially results in a lung abscess if not properly untreated. Here in this article we introduced a novel mathematical model to investigate the potential impact of Pneumonia treatments on disease transmission dynamics. The model is then validated against data from Jakarta City, Indonesia. In the model, the infection stage in infected individuals is categorized into three stages: the Exposed, Congestion and Hepatization, and the Resolution stage. Mathematical analysis shows that the disease-free equilibrium is always locally asymptotically stable when the basic reproduction number is less than one and unstable when larger than one. The endemic equilibrium only exists when the basic reproduction number is larger than one. Our proposed model always exhibits a forward bifurcation when the basic reproduction number is equal to one, which indicates local stability of the endemic equilibrium when the basic reproduction number is larger than one but close to one. A global sensitivity analysis shows that the infection parameter is the most influential parameter in determining the size of the total infected individual in the endemic equilibrium point. Furthermore, we also found that the hospitalization and the acceleration of the treatment duration can be used to control the level of endemic size. An optimal control problem was constructed from the earlier model and analyzed using the Pontryagin Maximum Principle. We find that the implementation of treatment in the earlier stage of infected individuals is needed to avoid a more significant outbreak of Pneumonia in a long-term intervention.
Collapse
Affiliation(s)
- Dipo Aldila
- Department of Mathematics, Universitas Indonesia, Depok 16424, Indonesia
| | - Nadya Awdinda
- Department of Mathematics, Universitas Indonesia, Depok 16424, Indonesia
| | - Fatmawati
- Department of Mathematics, Faculty of Science and Technology, Universitas Airlangga, Surabaya 60115, Indonesia
| | - Faishal F. Herdicho
- Department of Mathematics, Faculty of Science and Technology, Universitas Airlangga, Surabaya 60115, Indonesia
| | - Meksianis Z. Ndii
- Department of Mathematics, University of Nusa Cendana, Kupang-NTT 85361, Indonesia
| | - Chidozie W. Chukwu
- Department of Mathematics, Wake Forest University, Winston-Salem, NC 27109, USA
| |
Collapse
|