Search

In this paper, we suggest a mathematical model of COVID-19 with multiple variants of the virus under optimal control. Mathematical modeling has been used to gain deeper insights into the transmission of COVID-19, and various prevention and control strategies have been implemented to mitigate its spread. Our model is a SEIR-based model for multi-strains of COVID-19 with 7 compartments. We also consider the circulatory structure to account for the termination of immunity for COVID-19. The model is established in terms of the positivity and boundedness of the solution and the existence of equilibrium points, and the local stability of the solution. As a result of fitting data of COVID-19 in Ghana to the model, the basic reproduction number of the original virus and Delta variant was estimated to be 1.9396, and the basic reproduction number of the Omicron variant was estimated to be 3.4905, which is 1.8 times larger than that. We observe that even small differences in the incubation and recovery periods of two strains with the same initial transmission rate resulted in large differences in the number of infected individuals. In the case of COVID-19, infections caused by the Omicron variant occur 1.5 to 10 times more than those caused by the original virus. In terms of the optimal control strategy, we formulate three control strategies focusing on social distancing, vaccination, and testing-treatment. We have developed an optimal control model for the three strategies outlined above for the multi-strain model using the Pontryagin's Maximum Principle. Through numerical simulations, we analyze three optimal control strategies for each strain and also consider combinations of the two control strategies. As a result of the simulation, all control strategies are effective in reducing disease spread, in particular, vaccination strategies are more effective than the other two control strategies. In addition the combination of the two strategies also reduces the number of infected individuals by 1/10 compared to implementing one strategy, even when mild levels are implemented. Finally, we show that if the testing-treatment strategy is not properly implemented, the number of asymptomatic and unidentified infections may surge. These results could help guide the level of government intervention and prevention strategy formulation., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Kim et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The time-varying effective reproduction number Rt is a widely used indicator of transmission dynamics during infectious disease outbreaks. Timely estimates of Rt can be obtained from reported cases counted by their date of symptom onset, which is generally closer to the time of infection than the date of report. Case counts by date of symptom onset are typically obtained from line list data, however these data can have missing information and are subject to right truncation. Previous methods have addressed these problems independently by first imputing missing onset dates, then adjusting truncated case counts, and finally estimating the effective reproduction number. This stepwise approach makes it difficult to propagate uncertainty and can introduce subtle biases during real-time estimation due to the continued impact of assumptions made in previous steps. In this work, we integrate imputation, truncation adjustment, and Rt estimation into a single generative Bayesian model, allowing direct joint inference of case counts and Rt from line list data with missing symptom onset dates. We then use this framework to compare the performance of nowcasting approaches with different stepwise and generative components on synthetic line list data for multiple outbreak scenarios and across different epidemic phases. We find that under reporting delays realistic for hospitalization data (50% of reports delayed by more than a week), intermediate smoothing, as is common practice in stepwise approaches, can bias nowcasts of case counts and Rt, which is avoided in a joint generative approach due to shared regularization of all model components. On incomplete line list data, a fully generative approach enables the quantification of uncertainty due to missing onset dates without the need for an initial multiple imputation step. In a real-world comparison using hospitalization line list data from the COVID-19 pandemic in Switzerland, we observe the same qualitative differences between approaches. The generative modeling components developed in this work have been integrated and further extended in the R package epinowcast, providing a flexible and interpretable tool for real-time surveillance. Author summary: During an infectious disease outbreak, public health authorities require timely indicators of transmission dynamics, such as the effective reproduction number Rt. Since reporting data are delayed and often incomplete, statistical methods must be employed to obtain real-time estimates of case numbers and Rt. Existing methods involve separate steps for imputing missing data, adjusting for reporting delays, and estimating Rt. This stepwise approach impedes uncertainty quantification and can lead to inconsistent smoothing assumptions across steps. In this paper, we propose an alternative approach based on generative Bayesian modeling which integrates all steps into a single nowcasting model that can be directly fit to observed data. Using synthetic and real-world line list data, we demonstrate that the generative approach better captures uncertainty and avoids bias from inconsistent assumptions. The model components of our approach have been integrated into the R package epinowcast for easy use in practice. [ABSTRACT FROM AUTHOR]

The geographical range of schistosomiasis is affected by the ecology of schistosome parasites and their obligate host snails, including their response to temperature. Previous models predicted schistosomiasis' thermal optimum at 21.7°C, which is not compatible with the temperature in sub-Saharan Africa (SSA) regions where schistosomiasis is hyperendemic. We performed an extensive literature search for empirical data on the effect of temperature on physiological and epidemiological parameters regulating the free-living stages of S. mansoni and S. haematobium and their obligate host snails, i.e., Biomphalaria spp. and Bulinus spp., respectively. We derived nonlinear thermal responses fitted on these data to parameterize a mechanistic, process-based model of schistosomiasis. We then re-cast the basic reproduction number and the prevalence of schistosome infection as functions of temperature. We found that the thermal optima for transmission of S. mansoni and S. haematobium range between 23.1–27.3°C and 23.6–27.9°C (95% CI) respectively. We also found that the thermal optimum shifts toward higher temperatures as the human water contact rate increases with temperature. Our findings align with an extensive dataset of schistosomiasis prevalence in SSA. The refined nonlinear thermal-response model developed here suggests a more suitable current climate and a greater risk of increased transmission with future warming for more than half of the schistosomiasis suitable regions with mean annual temperature below the thermal optimum. Author summary: In this research, we explored the complex interplay between temperature and the transmission risk of schistosomiasis, a parasitic disease currently affecting over two hundred million people, predominantly in SSA. We developed a novel mathematical model accounting for the multiple positive and negative ways temperature affects the free-living stages of the parasite and its obligate, non-human host, i.e., specific species of freshwater snails. Our models show that schistosomiasis transmission risk peaks at temperatures 1–6°C higher than previously estimated. This indicates that the impact of climate change on schistosomiasis transmission might be more extensive than previously thought, affecting a wide geographic range where mean annual temperatures are currently below the optimal temperature. Our model projections are consistent with the observed temperatures in locations of SSA where schistosomiasis is endemic and data on infection prevalence in the human population are available. These findings suggest that the current climate is conducive to schistosomiasis transmission, and future warming could escalate the risk further, emphasizing the need for targeted interventions in these regions. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the spread of COVID-19 and the impact of government measures at the early stage of the pandemic (before the introduction of the vaccines) in the Netherlands. We build a multiple linear regression model to predict the effective reproduction rate using key factors and measures and integrate it with a system dynamics model to predict the spread and the impact of measures against COVID-19. Data from February to November 2020 is used to train the model and data until December 2020 is used to validate the model. We use data about the key factors, e.g., disease specific such as basic reproduction rate and incubation period, weather related factors such as temperature, and controllable factors such as testing capacity. We consider particularly the following measures taken by the government: wearing facemasks, event allowance, school closure, catering services closure, and self-quarantine. Studying the strategy of the Dutch government, we control these measures by following four main policies: doing nothing, mitigation, curbing, elimination. We develop a systems dynamic model to simulate the effect of policies. Based on our numerical experiments, we develop the following main insights: It is more effective to implement strict, sharp measures earlier but for a shorter duration than to introduce measures gradually for a longer duration. This way, we can prevent a quick rise in the number of infected cases but also to reduce the number of days under measures. Combining the measures with a high testing capacity and with effective self-quarantine can significantly reduce the spread of COVID-19. [ABSTRACT FROM AUTHOR]

HIV/AIDS and COVID-19 co-infection is a common global health and socio-economic problem. In this paper, a mathematical model for the transmission dynamics of HIV/AIDS and COVID-19 co-infection that incorporates protection and treatment for the infected (and infectious) groups is formulated and analyzed. Firstly, we proved the non-negativity and boundedness of the co-infection model solutions, analyzed the single infection models steady states, calculated the basic reproduction numbers using next generation matrix approach and then investigated the existence and local stabilities of equilibriums using Routh-Hurwiz stability criteria. Then using the Center Manifold criteria to investigate the proposed model exhibited the phenomenon of backward bifurcation whenever its effective reproduction number is less than unity. Secondly, we incorporate time dependent optimal control strategies, using Pontryagin's Maximum Principle to derive necessary conditions for the optimal control of the disease. Finally, we carried out numerical simulations for both the deterministic model and the model incorporating optimal controls and we found the results that the model solutions are converging to the model endemic equilibrium point whenever the model effective reproduction number is greater than unity, and also from numerical simulations of the optimal control problem applying the combinations of all the possible protection and treatment strategies together is the most effective strategy to drastically minimizing the transmission of the HIV/AIDS and COVID-19 co-infection in the community under consideration of the study. [ABSTRACT FROM AUTHOR]

In a viral epidemic, the emergence of a novel strain with increased transmissibility (larger value of basic reproduction number R0) sparks the fear that the increase in transmissibility is likely to lead to an increase in disease severity. It is required to investigate if a new, more contagious strain will be necessarily dominant in the population and resulting in more disease severity. In this paper, the impact of the asymptomatic transmission and the emergence time of a more transmissible variant of a multi-strain viral disease on the disease prevalence, disease severity, and the dominant variant in an epidemic was investigated by a proposed 2-strain epidemic model. The simulation results showed that considering only R0, is insufficient to predict the outcome of a new, more contagious strain in the population. A more transmissible strain with a high fraction of asymptomatic cases can substantially reduce the mortality rate. If the emergence time of the new strain is closer to the start of the epidemic, the new, more contagious variant has more chance to win the viral competition and be the dominant strain; otherwise, despite being more contagious, it cannot dominate previous strains. In conclusion, three factors of R0, the fraction of asymptomatic transmission, and the emergence time of the new strain are required to correctly determine the prevalence, disease severity, and the winner of the viral competition. [ABSTRACT FROM AUTHOR]

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) took the world by surprise. Following the first outbreak of COVID-19 in December 2019, several models have been developed to study and understand its transmission dynamics. Although the spread of COVID-19 is being slowed down by vaccination and other interventions, there is still a need to have a clear understanding of the evolution of the pandemic across countries, states and communities. To this end, there is a need to have a clearer picture of the initial spread of the disease in different regions. In this project, we used a simple SEIR model and a Bayesian inference framework to estimate the basic reproduction number of COVID-19 across Africa. Our estimates vary between 1.98 (Sudan) and 9.66 (Mauritius), with a median of 3.67 (90% CrI: 3.31-4.12). The estimates provided in this paper will help to inform COVID-19 modeling in the respective countries/regions., Competing Interests: The authors have declared that no competing interests exist.

To control the HIV/AIDS epidemics in Guangxi Zhuang Autonomous Region in China, Guangxi government launched the 5-year Guangxi AIDS Conquering Project (GACP, Phase I: 2010-2014, Phase II: 2015-2020). In the project, three measures are implemented, such as great improvements of the coverage of HIV/AIDS education, promotion of HIV voluntary counseling and testing, and enhancement of antiretroviral treatment. In this paper, we explore the effects of the three measures of GACP by construction of a Susceptible-Infected-Diagnosed-Treated population compartments model and via evaluation of the basic reproduction number derived from the model. A computational framework is developed for estimating the model parameters based on the HIV surveillance data, with application of the Markov-Chain Monte-Carlo method and Nonlinear Least Squares method. By estimating the new infections and evaluating the basic reproduction number, we find that the implementation of the three measures of GACP has a significant effect on controlling the rise of HIV/AIDS cases and the epidemic trend. Compared with HIV voluntary counseling and testing, strengthening HIV/AIDS education and expanding the coverage of antiretroviral treatment show a greater impact on HIV/AIDS epidemic control, which provides a reference project for other provinces with a similar epidemic situation in Guangxi Zhuang Autonomous Region. At the same time, our research fills the current research gap for the evaluation of large-scale AIDS prevention and control projects in developing areas. [ABSTRACT FROM AUTHOR]

Information transmission significantly impacts social stability and technological advancement. This paper compares the phenomenon of "Super transmission" and "Asymptomatic infection" in COVID-19 transmission to information transmission. The former is similar to authoritative information transmission individuals, whereas the latter is similar to individuals with low acceptance in information transmission. It then constructs an S2EIR model with transmitter authority and individual acceptance levels. Then, it analyzes the asymptotic stability of information-free and information-existence equilibrium on a local and global scale, as well as the model's basic reproduction number, R0. Distinguished with traditional studies, the population density function and Hamiltonian function are constructed by taking proportion of "Super transmitter" and proportion of hesitant group turning into transmitters as optimization control variables. Based on the Pontryagin maximum principle, an optimal control strategy is designed to effectively facilitate information transmission. The numerical simulation corroborates the theoretical analysis results and the system's sensitivity to control parameter changes. The research results indicate that the authoritative "Super transmitter" has a beneficial effect on information transmission. In contrast, the "Asymptomatic infected individual" with poor individual acceptance level negatively affects information transmission. [ABSTRACT FROM AUTHOR]

Estimating the basic reproduction number of a pandemic and the changes that appear on this value over time provide a good understanding of the contagious nature of the virus and efficiency of the controlling strategies. In this paper, we focus on studying the basic reproduction number (R0) for two important variants of COVID-19 pandemic in Iran: Alpha and Delta variants. We use four different methods, three statistical models and one mathematical model, to compute R0: Exponential Growth Rate (EGR), Maximum Likelihood (ML), Sequential Bayesian (SB), and time-dependent SIR model. Alpha variant of COVID-19 was active in Iran from March 10, 2021 until June 10, 2021. Our computations indicate that total R0 of this variant according to EGR, ML, SB, and SIR model is respectively 0.9999 (95% CI: 0.9994-1), 1.046 (95% CI: 1.044-1.049), 1.06 (95% CI: 1.03-1.08), and 2.79 (95% CI: 2.77-2.81) in the whole active time interval. Moreover, during the time interval from April 3, 2021 to April 9, 2021 in which this variant was in its exponential growth in Iran, R0 of Alpha variant in Iran according to SB, EGR, ML, and SIR model is respectively 2.26 (95% CI: 2.04-2.49), 2.64 (95% CI: 2.58-2.7), 11.38 (95% CI: 11.28-11.48), and 12.13 (95% CI: 12.12-12.14). Delta variant was active in Iran during the time interval from June 22, 2021 until September 22, 2021. Our computations show that during the time interval from July 3, 2021 to July 8, 2021 in which this variant was in its exponential growth in Iran, R0 of Delta variant in Iran according to SB, EGR, ML, and SIR model is respectively 3 (95% CI: 2.34-3.66), 3.1 (95% CI: 3.02-3.17), 12 (95% CI: 11.89-12.12), and 23.3 (95% CI: 23.19-23.41). Further, total R0 of Delta variant in Iran in the whole active time interval according to EGR, ML, SB, and SIR model is respectively 1.042 (95% CI: 1.04-1.043), 1.053 (95% CI: 1.051-1.055), 0.79 (95% CI: 0.63-0.95), and 5.65 (95% CI: 5.6-5.7). As the results show Delta variant was more severe than Alpha variant in Iran. Chasing the changes in R0 during each variant shows that the controlling strategies applied were effective in controlling the virus spread. [ABSTRACT FROM AUTHOR]

This paper uses publicly available data and various statistical models to estimate the basic reproduction number (R0) and other disease parameters for Ghana's early COVID-19 pandemic outbreak. We also test the effectiveness of government imposition of public health measures to reduce the risk of transmission and impact of the pandemic, especially in the early phase. R0 is estimated from the statistical model as 3.21 using a 0.147 estimated growth rate [95% C.I.: 0.137–0.157] and a 15-day time to recovery after COVID-19 infection. This estimate of the initial R0 is consistent with others reported in the literature from other parts of Africa, China and Europe. Our results also indicate that COVID-19 transmission reduced consistently in Ghana after the imposition of public health interventions—such as border restrictions, intra-city movement, quarantine and isolation—during the first phase of the pandemic from March to May 2020. However, the time-dependent reproduction number (Rt) beyond mid-May 2020 does not represent the true situation, given that there was not a consistent testing regime in place. This is also confirmed by our Jack-knife bootstrap estimates which show that the positivity rate over-estimates the true incidence rate from mid-May 2020. Given concerns about virus mutations, delays in vaccination and a possible new wave of the pandemic, there is a need for systematic testing of a representative sample of the population to monitor the reproduction number. There is also an urgent need to increase the availability of testing for the general population to enable early detection, isolation and treatment of infected individuals to reduce progression to severe disease and mortality. [ABSTRACT FROM AUTHOR]

In this paper, a new mathematical model is formulated that describes the interaction between uninfected cells, infected cells, viruses, intracellular viral RNA, Cytotoxic T-lymphocytes (CTLs), and antibodies. Hence, the model contains certain biological relations that are thought to be key factors driving this interaction which allow us to obtain precise logical conclusions. Therefore, it improves our perception, that would otherwise not be possible, to comprehend the pathogenesis, to interpret clinical data, to control treatment, and to suggest new relations. This model can be used to study viral dynamics in patients for a wide range of infectious diseases like HIV, HPV, HBV, HCV, and Covid-19. Though, analysis of a new multiscale HCV model incorporating the immune system response is considered in detail, the analysis and results can be applied for all other viruses. The model utilizes a transformed multiscale model in the form of ordinary differential equations (ODE) and incorporates into it the interaction of the immune system. The role of CTLs and the role of antibody responses are investigated. The positivity of the solutions is proven, the basic reproduction number is obtained, and the equilibrium points are specified. The stability at the equilibrium points is analyzed based on the Lyapunov invariance principle. By using appropriate Lyapunov functions, the uninfected equilibrium point is proven to be globally asymptotically stable when the reproduction number is less than one and unstable otherwise. Global stability of the infected equilibrium points is considered, and it has been found that each equilibrium point has a specific domain of stability. Stability regions could be overlapped and a bistable equilibria could be found, which means the coexistence of two stable equilibrium points. Hence, the solution converges to one of them depending on the initial conditions. [ABSTRACT FROM AUTHOR]

This study presents a reliable mathematical model to explain the spread of typhoid fever, covering stages of susceptibility, infection, carrying, and recovery, specifically in the Sheno town community. A detailed analysis is done to ensure the solutions are positive, stay within certain limits, and are stable for both situations where the disease is absent and where it is consistently present. The Routh-Hurwitz stability criterion has been used and applied for the purpose of stability analysis. Using the next-generation matrix, we determined the intrinsic potential for disease transmission. It showing that typhoid fever is spreading actively in Sheno town, with cases above a critical level. Our findings reveal the instability of the disease-free equilibrium point alongside the stability of the endemic equilibrium point. We identified two pivotal factors for transmission of the disease: the infectious rate, representing the speed of disease transmission, and the recruitment rate, indicating the rate at which new individuals enter the susceptible population. These parameters are indispensable for devising effective control measures. It is imperative to keep these parameters below specific thresholds to maintain a basic reproduction number favorable for disease control. Additionally, the study carefully examines how different factors affect the spread of typhoid fever, giving a detailed understanding of its dynamics. At the end, this study provides valuable insights and specific strategies for managing the disease in the Sheno town community. [ABSTRACT FROM AUTHOR]

Hepatitis B virus (HBV) infection is a global public health issue. We offer a comprehensive analysis of the dynamics of HBV, which can be successfully controlled with vaccine and treatment. Hepatitis B virus (HBV) causes a significantly more severe and protracted disease compared to hepatitis A. While it initially presents as an acute disease, in approximately 5 to 10% of cases, it can develop into a chronic disease that causes permanent damage to the liver. The hepatitis B virus can remain active outside the body for at least seven days. If the virus penetrates an individual's body without immunization, it may still result in infection. Upon exposure to HBV, the symptoms often last for a duration ranging from 10 days to 6 months. In this study, we developed a new model for Hepatitis B Virus (HBV) that includes asymptomatic carriers, vaccination, and treatment classes to gain a comprehensive knowledge of HBV dynamics. The basic reproduction number R0 is calculated to identify future recurrence. The local and global stabilities of the proposed model are evaluated for values of R0 that are both below and above 1. The Lyapunov function is employed to ensure the global stability of the HBV model. Further, the existence and uniqueness of the proposed model are demonstrated. To look at the solution of the proposed model graphically, we used a useful numerical strategy, such as the non-standard finite difference method, to obtain more thorough numerical findings for the parameters that have a significant impact on disease elimination. In addition, the study of treatment class in the population, we may assess the effectiveness of alternative medicines to treat infected populations can be determined. Numerical simulations and graphical representations are employed to illustrate the implications of our theoretical conclusions. [ABSTRACT FROM AUTHOR]

We construct a model to investigate HIV/AIDS dynamics in real cases and study its mathematical analysis. The study examines the qualitative outcomes and confirms the local and global asymptotic stability of both the endemic equilibrium and the disease-free equilibrium. The model's criteria for exhibiting both local and global asymptotically stable behavior are examined. We compute the endemic equilibria and obtain the existence of a unique positive endemic equilibrium. The data is fitted to the model using the idea of nonlinear least-squares fitting. Accurate parameter values are achieved by fitting the data to the model using a 95% confidence interval. The basic reproduction number is computed using parameters that have been fitted or estimated. Sensitivity analysis is performed to discover the influential parameters that impact the reproduction number and the eradication of the disease. The results show that implementing preventive measures can reduce HIV/AIDS cases. [ABSTRACT FROM AUTHOR]

Epidemic or pathogen emergence is the phenomenon by which a poorly transmissible pathogen finds its evolutionary pathway to become a mutant that can cause an epidemic. Many mathematical models of pathogen emergence rely on branching processes. Here, we discuss pathogen emergence using Markov chains, for a more tractable analysis, generalizing previous work by Kendall and Bartlett about disease invasion. We discuss the probability of emergence failure for early epidemics, when the number of infected individuals is small and the number of the susceptible individuals is virtually unlimited. Our formalism addresses both directly transmitted and vector-borne diseases, in the cases where the original pathogen is 1) one step-mutation away from the epidemic strain, and 2) undergoing a long chain of neutral mutations that do not change the epidemiology. We obtain analytic results for the probabilities of emergence failure and two features transcending the transmission mechanism. First, the reproduction number of the original pathogen is determinant for the probability of pathogen emergence, more important than the mutation rate or the transmissibility of the emerged pathogen. Second, the probability of mutation within infected individuals must be sufficiently high for the pathogen undergoing neutral mutations to start an epidemic, the mutation threshold depending again on the basic reproduction number of the original pathogen. Finally, we discuss the parameterization of models of pathogen emergence, using SARS-CoV1 as an example of zoonotic emergence and HIV as an example for the emergence of drug resistance. We also discuss assumptions of our models and implications for epidemiology. [ABSTRACT FROM AUTHOR]

West Nile virus (WNV) is a vector-borne flavivirus that causes an increasing number of human and equine West Nile fever cases in Europe. While the virus has been present in the Mediterranean basin and the Balkans since the 1960s, recent years have witnessed its northward expansion, with the first human cases reported in Germany in 2018 and the Netherlands in 2020. WNV transmission and amplification within mosquitoes are temperature-dependent. This study applies a mathematical modelling approach to assess the conditions under which WNV circulation occurs based on the proportion of mosquito bites on WNV-competent birds (dilution), vector-host ratios, mosquito season length and the observed daily temperature data. We modelled five distinct European regions where previous WNV circulation has been observed within the Netherlands, Germany, Spain, Italy, and Greece. We observed that the number of days in which the basic reproduction number (R0) is above one, increased over the last 40 years in all five regions. In the Netherlands, the number of days in which the R0 is above one, is 70% lower than in Spain. The temperature in Greece, Spain and Italy allowed for circulation under low vector-host ratios, and at a high dilution. On the other hand in the Netherlands and Germany, given the observed daily temperature, the thresholds for circulation requires a lower dilution and higher vector-host ratios. For the Netherlands, a short window of introductions between late May and mid-June would result in detectable outbreaks. Our findings revealed that the temperate maritime climate of the Netherlands allows WNV circulation primarily during warmer summers, and only under high vector-host ratios. This research contributes valuable insights into the dynamic relationship between temperature, vector properties, and WNV transmission, offering guidance for proactive strategies in addressing this emerging health threat in Europe. Author summary: West Nile virus (WNV) is a virus that is transmitted by mosquitoes, and that leads to West Nile fever cases in humans and horses. Although the virus has been in certain parts of Europe since the 1960s, it has recently spread northward, with the first human cases reported in Germany in 2018 and the Netherlands in 2020. Here, we use mathematical modelling to understand the conditions necessary for WNV transmission. We consider the proportion of mosquito bites on birds that can carry the virus, the number of hosts and vectors, the length of mosquito season, and temperature data in specific European regions where WNV has been observed previously (Netherlands, Germany, Spain, Italy, and Greece). We found that the number of days suitable for WNV transmission (determined by the basic reproduction number, R0) has increased over the last 40 years in all these regions. Furthermore, we found differences among these regions. In the Netherlands, for instance, the conditions for virus circulation requires higher vector-host ratios, while in warmer regions like Greece, Spain, and Italy, circulation occurs under lower ratios. Understanding the connection between temperature, mosquito traits, and WNV transmission is crucial for implementing surveillance and preventive measures. [ABSTRACT FROM AUTHOR]

During the COVID-19 pandemic, forecasting COVID-19 trends to support planning and response was a priority for scientists and decision makers alike. In the United States, COVID-19 forecasting was coordinated by a large group of universities, companies, and government entities led by the Centers for Disease Control and Prevention and the US COVID-19 Forecast Hub (https://covid19forecasthub.org). We evaluated approximately 9.7 million forecasts of weekly state-level COVID-19 cases for predictions 1–4 weeks into the future submitted by 24 teams from August 2020 to December 2021. We assessed coverage of central prediction intervals and weighted interval scores (WIS), adjusting for missing forecasts relative to a baseline forecast, and used a Gaussian generalized estimating equation (GEE) model to evaluate differences in skill across epidemic phases that were defined by the effective reproduction number. Overall, we found high variation in skill across individual models, with ensemble-based forecasts outperforming other approaches. Forecast skill relative to the baseline was generally higher for larger jurisdictions (e.g., states compared to counties). Over time, forecasts generally performed worst in periods of rapid changes in reported cases (either in increasing or decreasing epidemic phases) with 95% prediction interval coverage dropping below 50% during the growth phases of the winter 2020, Delta, and Omicron waves. Ideally, case forecasts could serve as a leading indicator of changes in transmission dynamics. However, while most COVID-19 case forecasts outperformed a naïve baseline model, even the most accurate case forecasts were unreliable in key phases. Further research could improve forecasts of leading indicators, like COVID-19 cases, by leveraging additional real-time data, addressing performance across phases, improving the characterization of forecast confidence, and ensuring that forecasts were coherent across spatial scales. In the meantime, it is critical for forecast users to appreciate current limitations and use a broad set of indicators to inform pandemic-related decision making. Author summary: As SARS-CoV-2 began to spread throughout the world in early 2020, modelers played a critical role in predicting how the epidemic could take shape. Short-term forecasts of epidemic outcomes (for example, infections, cases, hospitalizations, or deaths) provided useful information to support pandemic planning, resource allocation, and intervention. Yet, infectious disease forecasting is still a nascent science, and the reliability of different types of forecasts is unclear. We retrospectively evaluated COVID-19 case forecasts, which were often unreliable. For example, forecasts did not anticipate the speed of increase in cases in early winter 2020. This analysis provides insights on specific problems that could be addressed in future research to improve forecasts and their use. Identifying the strengths and weaknesses of forecasts is critical to improving forecasting for current and future public health responses. [ABSTRACT FROM AUTHOR]

Background: Leishmaniasis are a group of diseases caused by more than 20 species of the protozoan that are transmitted through the bite of female sand fly. The disease is endemic to 98 countries of the world. It affects most commonly the poorest of the poor and mainly males. Several research has been conducted to propose disease control strategies. Effective medical care, vector control, environmental hygiene, and personal protection are the mainstays of the current preventative and control methods. The mathematical models for the transmission dynamics of the disease studied so far did not consider the sex-biased burden of the disease into consideration. Methodology: Unlike the previous VL works, this study introduces a new deterministic sex-structured model for understanding the transmission dynamics of visceral leishmaniasis. Basic properties of the model including basic reproduction number (R0), and conditions for the existence of backward bifurcation of the model are explored. Baseline parameter values were estimated after the model was fitted to Ethiopia's VL data. Sensitivity analysis of the model was performed to identify the parameters that significantly impact the disease threshold. Numerical simulations were performed using baseline parameter values, and scenario analysis is performed by changing some of these parameters as appropriate. Conclusion: The analysis of the model shows that there is a possibility for a backward bifurcation for R0<1 , which means bringing R0 to less than unity may not be enough to eradicate the disease. Our numerical result shows that the implementation of disease-preventive strategies, as well as effectively treating the affected ones can significantly reduce the disease prevalence if applied for more proportion of the male population. Furthermore, the implementation of vector management strategies also can considerably reduce the total prevalence of the disease. However, it is demonstrated that putting more effort in treating affected reservoir animals may not have any significant effect on the overall prevalence of the disease as compared to other possible mechanisms. The numerical simulation infers that a maximum of 60% of extra preventative measures targeted to only male population considerably reduces the total prevalence of VL by 80%. It is also possible to decrease the total prevalence of VL by 69.51% when up to 50% additional infected males receive treatment with full efficacy. Moreover, applying a maximum of 15% additional effort to reduce the number of vectors, decreases the total VL prevalence by 57.71%. Therefore, in order to reduce the disease burden of visceral leishmaniasis, public health officials and concerned stakeholders need to give more emphasis to the proportion of male humans in their intervention strategies. [ABSTRACT FROM AUTHOR]

We formulate a mathematical model of social media addiction and depression (SMAD) in this study. Key aspects, such as social media addiction and depression disease-free equilibrium point (SMADDFEP), social media addiction and depression endemic equilibrium point (SMADEEP), and basic reproduction number (R0), have been analyzed qualitatively. The results indicate that if R0 < 1, the SMADDFEP is locally asymptotically stable. The global asymptotic stability of the SMADDFEP has been established using the Castillo-Chavez theorem. On the other hand, if R0 > 1, the unique endemic equilibrium point (SMADEEP) is locally asymptotically stable by Lyapunov theorem, and the model exhibits a forward bifurcation at R0 = 1 according to the Center Manifold theorem. To examine the model's sensitivity, we calculated the normalized forward sensitivity index and conducted a Partial Rank Correlation Coefficient (PRCC) analysis to describe the influence of parameters on the SMAD. The numerical results obtained using the Fourth-order Runge-Kutta (RK-4) scheme show that increasing the number of addicted individuals leads to an increase in the number of depressed individuals. [ABSTRACT FROM AUTHOR]

Assessing public health intervention strategies is crucial for effectively managing dengue. While numerous studies have explored the impact of dengue interventions on its transmission dynamics, limited research has focused on the combined effects of implementing multiple therapeutic interventions for disease control. This study presents an epidemic model for understanding dengue transmission dynamics, incorporating two critical therapeutic measures: vaccination and treatment of infected individuals. The model is characterized by ordinary differential equations involving seven-state variables. The investigation encompasses both disease-free and endemic equilibria of the model. The findings reveal that the disease-free equilibrium (only) is globally stable when the basic reproduction number is below one. Interestingly, when the vaccine's effectiveness is low, treatment emerges as a more successful approach in reducing dengue cases than vaccination. In contrast, a highly effective vaccine alone significantly curtails dengue occurrences. Moreover, the study introduces an optimal control problem, featuring an objective function integrating two control mechanisms: vaccination and treatment. The analysis strongly suggests that implementing two control strategies outweighs the efficacy of a single approach in effectively mitigating the spread of the disease. [ABSTRACT FROM AUTHOR]

Despite use of tecovirimat since the beginning of the 2022 outbreak, few data have been published on its antiviral effect in humans. We here predict tecovirimat efficacy using a unique set of data in nonhuman primates (NHPs) and humans. We analyzed tecovirimat antiviral activity on viral kinetics in NHP to characterize its concentration–effect relationship in vivo. Next, we used a pharmacological model developed in healthy volunteers to project its antiviral efficacy in humans. Finally, a viral dynamic model was applied to characterize mpox kinetics in skin lesions from 54 untreated patients, and we used this modeling framework to predict the impact of tecovirimat on viral clearance in skin lesions. At human-recommended doses, tecovirimat could inhibit viral replication from infected cells by more than 90% after 3 to 5 days of drug administration and achieved over 97% efficacy at drug steady state. With an estimated mpox within-host basic reproduction number, R0, equal to 5.6, tecovirimat could therefore shorten the time to viral clearance if given before viral peak. We predicted that initiating treatment at symptom onset, which on average occurred 2 days before viral peak, could reduce the time to viral clearance by about 6 days. Immediate postexposure prophylaxis could not only reduce time to clearance but also lower peak viral load by more than 1.0 log10 copies/mL and shorten the duration of positive viral culture by about 7 to 10 days. These findings support the early administration of tecovirimat against mpox infection, ideally starting from the infection day as a postexposure prophylaxis. Despite its use since the beginning of the 2022 mpox outbreak, few data have been published on the antiviral effect of tecovirimat in humans. This study predicts the impact of early tecovirimat administration on the time to viral clearance in patients with mpox infection, using an integrative modeling approach that combines pre-clinical and clinical data. [ABSTRACT FROM AUTHOR]

Since 2014, Brazil has experienced an unprecedented epidemic caused by chikungunya virus (CHIKV), with several waves of East-Central-South-African (ECSA) lineage transmission reported across the country. In 2018, Rio de Janeiro state, the third most populous state in Brazil, reported 41% of all chikungunya cases in the country. Here we use evolutionary and epidemiological analysis to estimate the timescale of CHIKV-ECSA-American lineage and its epidemiological patterns in Rio de Janeiro. We show that the CHIKV-ECSA outbreak in Rio de Janeiro derived from two distinct clades introduced from the Northeast region in mid-2015 (clade RJ1, n = 63/67 genomes from Rio de Janeiro) and mid-2017 (clade RJ2, n = 4/67). We detected evidence for positive selection in non-structural proteins linked with viral replication in the RJ1 clade (clade-defining: nsP4-A481D) and the RJ2 clade (nsP1-D531G). Finally, we estimate the CHIKV-ECSA's basic reproduction number (R0) to be between 1.2 to 1.6 and show that its instantaneous reproduction number (Rt) displays a strong seasonal pattern with peaks in transmission coinciding with periods of high Aedes aegypti transmission potential. Our results highlight the need for continued genomic and epidemiological surveillance of CHIKV in Brazil, particularly during periods of high ecological suitability, and show that selective pressures underline the emergence and evolution of the large urban CHIKV-ECSA outbreak in Rio de Janeiro. Author summary: Chikungunya is a mosquito-borne viral disease that has emerged as a significant public health concern in many regions worldwide. The state of Rio de Janeiro, a key economic and tourism hub in Brazil, has experienced multiple outbreaks of chikungunya in recent years, resulting in significant morbidity and economic burden. However, our understanding of the establishment and epidemiology of chikungunya in Rio de Janeiro remains limited. We conducted an analysis of chikungunya epidemiology in Rio de Janeiro, focusing on the first four years of virus circulation in the state. We estimated the magnitude of chikungunya transmission during this period, providing evidence that the trends of incidence and transmissibility can be influenced by climatic fluctuations, which impact vector abundance. Through analysis of novel genetic data, we inferred the evolutionary history of the virus, providing evidence of molecular adaptation in the lineage circulating in Rio de Janeiro. These findings highlight the need for designing effective vector control strategies to prevent and mitigate chikungunya in Rio de Janeiro and similar settings. Further research is needed to continue monitoring viral epidemiological trends and to better understand the complex interactions between climatic factors, vector dynamics, and viral evolution in shaping the transmission of this emerging disease. [ABSTRACT FROM AUTHOR]

The urgency to develop vaccines during the COVID-19 pandemic has resulted in the acceleration of clinical trials. Specifically, a broad spectrum of efficacy levels has been reported for various vaccines based on phase III cohort studies. Our study demonstrates that conducting large cohort phase III clinical trials during the peak of an epidemic leads to a significant underestimation of vaccine efficacy, even in the absence of confounding factors. Furthermore, we find that this underestimation increases with the proportion of infectious individuals in the population during the experiment and the severity of the epidemic, as measured by its basic reproduction number. [ABSTRACT FROM AUTHOR]

African horse sickness is an equine orbivirus transmitted by Culicoides Latreille biting midges. In the last 80 years, it has caused several devastating outbreaks in the equine population in Europe, the Far and Middle East, North Africa, South-East Asia, and sub-Saharan Africa. The disease is endemic in South Africa; however, a unique control area has been set up in the Western Cape where increased surveillance and control measures have been put in place. A deterministic metapopulation model was developed to explore if an outbreak might occur, and how it might develop, if a latently infected horse was to be imported into the control area, by varying the geographical location and months of import. To do this, a previously published ordinary differential equation model was developed with a metapopulation approach and included a vaccinated horse population. Outbreak length, time to peak infection, number of infected horses at the peak, number of horses overall affected (recovered or dead), re-emergence, and Rv (the basic reproduction number in the presence of vaccination) were recorded and displayed using GIS mapping. The model predictions were compared to previous outbreak data to ensure validity. The warmer months (November to March) had longer outbreaks than the colder months (May to September), took more time to reach the peak, and had a greater total outbreak size with more horses infected at the peak. Rv appeared to be a poor predictor of outbreak dynamics for this simulation. A sensitivity analysis indicated that control measures such as vaccination and vector control are potentially effective to manage the spread of an outbreak, and shortening the vaccination window to July to September may reduce the risk of vaccine-associated outbreaks. Author summary: African Horse Sickness is a devastating global equine disease which results in high numbers of deaths. In South Africa it is considered endemic, and this poses implications for trade and exportation of horses out of the country. As a result, South Africa has developed a "control area" where there are increased surveillance and control strategies to keep it free from African horse sickness. Our research explores where and how an outbreak may occur in the "control area" if a horse harbouring the virus were to be imported into the region. We developed a mathematical simulation which considered the differences in horse population numbers, midge population numbers, environmental temperature, and vaccination status of horses in different areas across the "control area". This helped to determine how many horses would be affected at the end of the outbreak as well as at the peak, how long the outbreak would occur, and the reproduction number (Rv), and how this varied across the study area. Our research uncovers insights into the African horse sickness outbreak dynamics of the "control area" and helps provide information to further develop control strategies. [ABSTRACT FROM AUTHOR]

In recent years measles has been one of the most critical public health problem in Bangladesh. Although the Ministry of Health in Bangladesh employs a broad extension of measles control policies, logistical challenges exist, and there is significant doubt regarding the disease burden. Mathematical modelling of measles is considered one of the most effective ways to understand infection transmission and estimate parameters in different countries, such as Bangladesh. In this study, a mathematical modelling framework is presented to explore the dynamics of measles in Bangladesh. We calibrated the model using cumulative measles incidence data from 2000 to 2019. Also, we performed a sensitivity analysis of the model parameters and found that the contact rate had the most significant influence on the basic reproduction number R0. Four hypothetical intervention scenarios were developed and simulated for the period from 2020 to 2035. The results show that the scenario which combines enhanced treatment for exposed and infected population, first and second doses of vaccine is the most effective at rapidly reducing the total number of measles incidence and mortality in Bangladesh. Our findings also suggest that strategies that focus on a single interventions do not dramatically affect the decline in measles incidence cases; instead, those that combine two or more interventions simultaneously are the most effective in decreasing the burden of measles incidence and mortality. In addition, we also evaluated the cost-effectiveness of varying combinations of three basic control strategies including distancing, vaccination and treatment, all within the optimal control framework. Our finding suggested that combines distancing, vaccination and treatment control strategy is the most cost-effective for reducing the burden of measles in Bangladesh. Other strategies can be comprised to measles depending on the availability of funds and policymakers' choices. [ABSTRACT FROM AUTHOR]

Epidemics, such as COVID-19, have caused significant harm to human society worldwide. A better understanding of epidemic transmission dynamics can contribute to more efficient prevention and control measures. Compartmental models, which assume homogeneous mixing of the population, have been widely used in the study of epidemic transmission dynamics, while agent-based models rely on a network definition for individuals. In this study, we developed a real-scale contact-dependent dynamic (CDD) model and combined it with the traditional susceptible-exposed-infectious-recovered (SEIR) compartment model. By considering individual random movement and disease spread, our simulations using the CDD-SEIR model reveal that the distribution of agent types in the community exhibits spatial heterogeneity. The estimated basic reproduction number R0 depends on group mobility, increasing logarithmically in strongly heterogeneous cases and saturating in weakly heterogeneous conditions. Notably, R0 is approximately independent of virus virulence when group mobility is low. We also show that transmission through small amounts of long-term contact is possible due to short-term contact patterns. The dependence of R0 on environment and individual movement patterns implies that reduced contact time and vaccination policies can significantly reduce the virus transmission capacity in situations where the virus is highly transmissible (i.e., R0 is relatively large). This work provides new insights into how individual movement patterns affect virus spreading and how to protect people more efficiently. [ABSTRACT FROM AUTHOR]

Background: In March 2022, the Omicron variant of SARS-CoV-2 spread rapidly in Shanghai, China. The city adopted strict non-pharmacological intervention (NPI) measures, including lockdown (implemented on March 28 in Pudong and April 1 in Puxi) and blanket PCR testing (April 4). This study aims to understand the effect of these measures. Methods: We tabulated daily case counts from official reports and fitted a two-patch stochastic SEIR model to the data for the period of March 19 to April 21. This model considered two regions in Shanghai, namely Pudong and Puxi, as the implementation of control measures in Shanghai was carried out on different dates in these regions. We verified our fitting results using the data from April 22 to June 26. Finally, we applied the point estimate of parameter values to simulate our model while varying the dates of control measure implementation, and studied the effectiveness of the control measures. Results: Our point estimate for the parameter values yields expected case counts that agree well the data for both the periods from March 19 to April 21 and from April 22 to June 26. Lockdown did not significantly reduce the intra-region transmission rates. Only about 21% cases were reported. The underlying basic reproduction number R0 was 1.7, and the control reproduction number with both lockdown and blanket PCR testing was 1.3. If both measures were implemented on March 19, only about 5.9% infections would be prevented. Conclusions: Through our analysis, we found that NPI measures implemented in Shanghai were not sufficient to reduce the reproduction number to below unity. Thus, earlier intervention only has limited effect on reducing cases. The outbreak dies out because of only 27% of the population were active in disease transmission, possibly due to a combination of vaccination and lockdown. [ABSTRACT FROM AUTHOR]

For vector-borne diseases the basic reproduction number R0 , a measure of a disease's epidemic potential, is highly temperature-dependent. Recent work characterizing these temperature dependencies has highlighted how climate change may impact geographic disease spread. We extend this prior work by examining how newly emerging diseases, like Zika, will be impacted by specific future climate change scenarios in four diverse regions of Brazil, a country that has been profoundly impacted by Zika. We estimated a R0(T) , derived from a compartmental transmission model, characterizing Zika (and, for comparison, dengue) transmission potential as a function of temperature-dependent biological parameters specific to Aedes aegypti. We obtained historical temperature data for the five-year period 2015–2019 and projections for 2045–2049 by fitting cubic spline interpolations to data from simulated atmospheric data provided by the CMIP-6 project (specifically, generated by the GFDL-ESM4 model), which provides projections under four Shared Socioeconomic Pathways (SSP). These four SSP scenarios correspond to varying levels of climate change severity. We applied this approach to four Brazilian cities (Manaus, Recife, Rio de Janeiro, and São Paulo) that represent diverse climatic regions. Our model predicts that the R0(T) for Zika peaks at 2.7 around 30°C, while for dengue it peaks at 6.8 around 31°C. We find that the epidemic potential of Zika will increase beyond current levels in Brazil in all of the climate scenarios. For Manaus, we predict that the annual R0 range will increase from 2.1–2.5, to 2.3–2.7, for Recife we project an increase from 0.4–1.9 to 0.6–2.3, for Rio de Janeiro from 0–1.9 to 0–2.3, and for São Paulo from 0–0.3 to 0–0.7. As Zika immunity wanes and temperatures increase, there will be increasing epidemic potential and longer transmission seasons, especially in regions where transmission is currently marginal. Surveillance systems should be implemented and sustained for early detection. Author summary: Rising temperatures through climate change are expected to increase arboviral disease pressure, so understanding the impact of climate change on newly emerging diseases such as Zika is essential to prepare for future outbreaks. However, because disease transmission may be less effective at very high temperatures, it is uncertain whether risk will uniformly increase in different regions. Given the nonlinear relationship between temperature and many important biological vector traits, mathematical modeling is a useful tool for predicting the impact of temperature on arbovirus risk. We used a temperature-dependent infectious disease transmission model to derive a temperature-dependent basic reproduction number. We then used historical temperature data and temperature projections for the years 2045–2049 to forecast Zika risk in four cities in Brazil under various climate change scenarios. We predict an overall increase in arbovirus risk, as well as extended risk seasons in cities that are not currently suitable for year-round spread, such as Rio de Janeiro. We also found little-to-no protective effect of increasing temperatures even in warmer climates like Manaus. Our results indicate that preparation for future Zika outbreaks (and of those of other arboviruses including dengue) should include the implementation of national disease surveillance and early detection systems. [ABSTRACT FROM AUTHOR]

Superspreading events play an important role in the spread of several pathogens, such as SARS-CoV-2. While the basic reproduction number of the original Wuhan SARS-CoV-2 is estimated to be about 3 for Belgium, there is substantial inter-individual variation in the number of secondary cases each infected individual causes—with most infectious individuals generating no or only a few secondary cases, while about 20% of infectious individuals is responsible for 80% of new infections. Multiple factors contribute to the occurrence of superspreading events: heterogeneity in infectiousness, individual variations in susceptibility, differences in contact behavior, and the environment in which transmission takes place. While superspreading has been included in several infectious disease transmission models, research into the effects of different forms of superspreading on the spread of pathogens remains limited. To disentangle the effects of infectiousness-related heterogeneity on the one hand and contact-related heterogeneity on the other, we implemented both forms of superspreading in an individual-based model describing the transmission and spread of SARS-CoV-2 in a synthetic Belgian population. We considered its impact on viral spread as well as on epidemic resurgence after a period of social distancing. We found that the effects of superspreading driven by heterogeneity in infectiousness are different from the effects of superspreading driven by heterogeneity in contact behavior. On the one hand, a higher level of infectiousness-related heterogeneity results in a lower risk of an outbreak persisting following the introduction of one infected individual into the population. Outbreaks that did persist led to fewer total cases and were slower, with a lower peak which occurred at a later point in time, and a lower herd immunity threshold. Finally, the risk of resurgence of an outbreak following a period of lockdown decreased. On the other hand, when contact-related heterogeneity was high, this also led to fewer cases in total during persistent outbreaks, but caused outbreaks to be more explosive in regard to other aspects (such as higher peaks which occurred earlier, and a higher herd immunity threshold). Finally, the risk of resurgence of an outbreak following a period of lockdown increased. We found that these effects were conserved when testing combinations of infectiousness-related and contact-related heterogeneity. Author summary: To investigate the effect of different sources of superspreading on disease dynamics, we implemented superspreading driven by heterogeneity in infectiousness and heterogeneity in contact behavior into an individual-based model for the transmission of SARS-CoV-2 in the Belgian population. We compared the impact of both forms of superspreading in a scenario without interventions as well as in a scenario in which a period of strict social distancing (i.e. a lockdown) is followed by a period of partial release. We found that both forms of superspreading have very different effects. On the one hand, increasing the level of infectiousness-related heterogeneity led to less outbreaks being observed following the introduction of one infected individual in the population. Furthermore, final outbreak sizes decreased, and outbreaks became slower, with lower and later peaks, and a lower herd immunity threshold. Finally, the risk for resurgence of an outbreak following a period of lockdown also decreased. On the other hand, when contact-related heterogeneity was high, this also led to smaller final sizes, but caused outbreaks to be more explosive regarding other aspects (such as higher peaks that occurred earlier). The herd immunity threshold also increased, as did the risk of resurgence of outbreaks. [ABSTRACT FROM AUTHOR]

During the COVID-19 pandemic authorities have been striving to obtain reliable predictions for the spreading dynamics of the disease. We recently developed a multi-"sub-populations" (multi-compartments: susceptible, exposed, pre-symptomatic, infectious, recovered) model, that accounts for the spatial in-homogeneous spreading of the infection and shown, for a variety of examples, how the epidemic curves are highly sensitive to location of epicenters, non-uniform population density, and local restrictions. In the present work we test our model against real-life data from South Carolina during the period May 22 to July 22 (2020). During this period, minimal restrictions have been employed, which allowed us to assume that the local basic reproduction number is constant in time. We account for the non-uniform population density in South Carolina using data from NASA's Socioeconomic Data and Applications Center (SEDAC), and predict the evolution of infection heat-maps during the studied period. Comparing the predicted heat-maps with those observed, we find high qualitative resemblance. Moreover, the Pearson's correlation coefficient is relatively high thus validating our model against real-world data. We conclude that the model accounts for the major effects controlling spatial in-homogeneous spreading of the disease. Inclusion of additional sub-populations (compartments), in the spirit of several recently developed models for COVID-19, can be easily performed within our mathematical framework. [ABSTRACT FROM AUTHOR]

Multi-strain pandemics have emerged as a major concern. We introduce a new model for assessing the connection between multi-strain pandemics and mortality rate, basic reproduction number, and maximum infected individuals. The proposed model provides a general mathematical approach for representing multi-strain pandemics, generalizing for an arbitrary number of strains. We show that the proposed model fits well with epidemiological historical world health data over a long time period. From a theoretical point of view, we show that the increasing number of strains increases logarithmically the maximum number of infected individuals and the mean mortality rate. Moreover, the mean basic reproduction number is statistically identical to the single, most aggressive pandemic strain for multi-strain pandemics. [ABSTRACT FROM AUTHOR]

Background: During the Coronavirus Disease 2019 (COVID-19) pandemic, the United Kingdom government imposed public health policies in England to reduce social contacts in hopes of curbing virus transmission. We conducted a repeated cross-sectional study to measure contact patterns weekly from March 2020 to March 2021 to estimate the impact of these policies, covering 3 national lockdowns interspersed by periods of less restrictive policies.Methods and Findings: The repeated cross-sectional survey data were collected using online surveys of representative samples of the UK population by age and gender. Survey participants were recruited by the online market research company Ipsos MORI through internet-based banner and social media ads and email campaigns. The participant data used for this analysis are restricted to those who reported living in England. We calculated the mean daily contacts reported using a (clustered) bootstrap and fitted a censored negative binomial model to estimate age-stratified contact matrices and estimate proportional changes to the basic reproduction number under controlled conditions using the change in contacts as a scaling factor. To put the findings in perspective, we discuss contact rates recorded throughout the year in terms of previously recorded rates from the POLYMOD study social contact study. The survey recorded 101,350 observations from 19,914 participants who reported 466,710 contacts over 53 weeks. We observed changes in social contact patterns in England over time and by participants' age, personal risk factors, and perception of risk. The mean reported contacts for adults 18 to 59 years old ranged between 2.39 (95% confidence interval [CI] 2.20 to 2.60) contacts and 4.93 (95% CI 4.65 to 5.19) contacts during the study period. The mean contacts for school-age children (5 to 17 years old) ranged from 3.07 (95% CI 2.89 to 3.27) to 15.11 (95% CI 13.87 to 16.41). This demonstrates a sustained decrease in social contacts compared to a mean of 11.08 (95% CI 10.54 to 11.57) contacts per participant in all age groups combined as measured by the POLYMOD social contact study in 2005 to 2006. Contacts measured during periods of lockdowns were lower than in periods of eased social restrictions. The use of face coverings outside the home has remained high since the government mandated use in some settings in July 2020. The main limitations of this analysis are the potential for selection bias, as participants are recruited through internet-based campaigns, and recall bias, in which participants may under- or overreport the number of contacts they have made.Conclusions: In this study, we observed that recorded contacts reduced dramatically compared to prepandemic levels (as measured in the POLYMOD study), with changes in reported contacts correlated with government interventions throughout the pandemic. Despite easing of restrictions in the summer of 2020, the mean number of reported contacts only returned to about half of that observed prepandemic at its highest recorded level. The CoMix survey provides a unique repeated cross-sectional data set for a full year in England, from the first day of the first lockdown, for use in statistical analyses and mathematical modelling of COVID-19 and other diseases. [ABSTRACT FROM AUTHOR]

COVID-19 vaccines have been approved for children of age five and older in many countries. However, there is an ongoing debate as to whether children should be vaccinated and at what priority. In this work, we use mathematical modeling and optimization to study how vaccine allocations to different age groups effect epidemic outcomes. In particular, we consider the effect of extending vaccination campaigns to include the vaccination of children. When vaccine availability is limited, we consider Pareto-optimal allocations with respect to competing measures of the number of infections and mortality and systematically study the trade-offs among them. In the scenarios considered, when some weight is given to the number of infections, we find that it is optimal to allocate vaccines to adolescents in the age group 10-19, even when they are assumed to be less susceptible than adults. We further find that age group 0-9 is included in the optimal allocation for sufficiently high values of the basic reproduction number. Author summary: The management of a vaccination campaign requires constant evaluation, and, as conditions change, redesign of the campaign for the future. One example of a change that leads to a decision point is the extension of vaccine eligibility to broader segments of the population. Indeed, one of the questions public health experts and policymakers confront is whether children of age 5 and older should be vaccinated against COVID-19 and at what priority. In this work, we use mathematical and computational methods to study this question in a systematic, quantitative way, by considering optimal outcomes, under different criteria for optimality, that can be achieved with and without vaccination of children, under conditions of limited vaccine availability. To explore tradeoffs among different goals, such as reducing the number of infections and reducing mortality, we use the the idea of Pareto optimization. Our results show that, under a broad range of conditions, optimal vaccine allocations include vaccination of age group 10–19, while the population-level benefits of vaccination of age group 0–9 depends on the reproductive number and vaccine availability. [ABSTRACT FROM AUTHOR]

Inferring the dynamics of pathogen transmission during an outbreak is an important problem in infectious disease epidemiology. In mathematical epidemiology, estimates are often informed by time series of confirmed cases, while in phylodynamics genetic sequences of the pathogen, sampled through time, are the primary data source. Each type of data provides different, and potentially complementary, insight. Recent studies have recognised that combining data sources can improve estimates of the transmission rate and the number of infected individuals. However, inference methods are typically highly specialised and field-specific and are either computationally prohibitive or require intensive simulation, limiting their real-time utility. We present a novel birth-death phylogenetic model and derive a tractable analytic approximation of its likelihood, the computational complexity of which is linear in the size of the dataset. This approach combines epidemiological and phylodynamic data to produce estimates of key parameters of transmission dynamics and the unobserved prevalence. Using simulated data, we show (a) that the approximation agrees well with existing methods, (b) validate the claim of linear complexity and (c) explore robustness to model misspecification. This approximation facilitates inference on large datasets, which is increasingly important as large genomic sequence datasets become commonplace. Author summary: Mathematical epidemiologists typically study time series of cases, ie the epidemic curve, to understand the spread of pathogens. Genetic epidemiologists study similar problems but do so using observations of the genetic sequence of the pathogen, which also contains information about the transmission process. There have been many attempts to unite these approaches and utilise both data sources. However, striking a suitable balance between model flexibility and fidelity, in a computationally tractable way, has proven challenging. There are several competing methods, but they are often intractable when applied to a large dataset. As sequencing of pathogen genomes becomes common, and an increasing amount of epidemiological data are collected, this situation will only be exacerbated. To bridge the gap between the time series and genomic methods we developed an approximation scheme, called TimTam. TimTam can accurately and efficiently estimate key features of an epidemic, such as the prevalence of infection (how many people are currently infected) and the basic reproduction number (a measure of the transmissibility of the infection.). [ABSTRACT FROM AUTHOR]

Many plant viruses are transmitted by insect vectors. Transmission can be described as persistent or non-persistent depending on rates of acquisition, retention, and inoculation of virus. Much experimental evidence has accumulated indicating vectors can prefer to settle and/or feed on infected versus noninfected host plants. For persistent transmission, vector preference can also be conditional, depending on the vector's own infection status. Since viruses can alter host plant quality as a resource for feeding, infection potentially also affects vector population dynamics. Here we use mathematical modelling to develop a theoretical framework addressing the effects of vector preferences for landing, settling and feeding–as well as potential effects of infection on vector population density–on plant virus epidemics. We explore the consequences of preferences that depend on the host (infected or healthy) and vector (viruliferous or nonviruliferous) phenotypes, and how this is affected by the form of transmission, persistent or non-persistent. We show how different components of vector preference have characteristic effects on both the basic reproduction number and the final incidence of disease. We also show how vector preference can induce bistability, in which the virus is able to persist even when it cannot invade from very low densities. Feedbacks between plant infection status, vector population dynamics and virus transmission potentially lead to very complex dynamics, including sustained oscillations. Our work is supported by an interactive interface https://plantdiseasevectorpreference.herokuapp.com/. Our model reiterates the importance of coupling virus infection to vector behaviour, life history and population dynamics to fully understand plant virus epidemics. Author summary: Plant virus diseases–which cause devastating epidemics in plant populations worldwide–are most often transmitted by insect vectors. Recent experimental evidence indicates how vectors do not choose between plants at random, but instead can be affected by whether plants are infected (or not). Virus infection can cause plants to "smell" different, because they produce different combinations of volatile chemicals, or "taste" different, due to chemical changes in infected tissues. Vector reproduction rates can also be affected when colonising infected versus uninfected plants. Potential effects on epidemic spread through a population of plants are not yet entirely understood. There are also interactions with the mode of virus transmission. Some viruses can be transmitted after only a brief probe by a vector, whereas others are only picked up after an extended feed on an infected plant. Furthermore there are differences in how long vectors remain able to transmit the virus. This ranges from a matter of minutes, right up to the entire lifetime of the insect, depending on the plant-virus-vector combination under consideration. Here we use mathematical modelling to synthesise all this complexity into a coherent theoretical framework. We illustrate our model via an online interface https://plantdiseasevectorpreference.herokuapp.com/. [ABSTRACT FROM AUTHOR]

In the past year, the global epidemic situation is still not optimistic, showing a trend of continuous expansion. With the research and application of vaccines, there is an urgent need to develop some optimal vaccination strategies. How to make a reasonable vaccination strategy to determine the priority of vaccination under the limited vaccine resources to control the epidemic and reduce human casualties? We build a dynamic model with vaccination which is extended the classical SEIR model. By fitting the epidemic data of three countries—China, Brazil, Indonesia, we have evaluated age-specific vaccination strategy for the number of infections and deaths. Furthermore, we have evaluated the impact of age-specific vaccination strategies on the number of the basic reproduction number. At last, we also have evaluated the different age structure of the vaccination priority. It shows that giving priority to vaccination of young people can control the number of infections, while giving priority to vaccination of the elderly can greatly reduce the number of deaths in most cases. Furthermore, we have found that young people should be mainly vaccinated to reduce the number of infections. When the emphasis is on reducing the number of deaths, it is important to focus vaccination on the elderly. Simulations suggest that appropriate age-specific vaccination strategies can effectively control the epidemic, both in terms of the number of infections and deaths. [ABSTRACT FROM AUTHOR]

There is a large burden of norovirus disease in child-care centers in Australia and around the world. Despite the ubiquity of norovirus outbreaks in child-care centers, little is known about the extent of this burden within the child-care center and the surrounding household clusters. Therefore, we performed an in-depth analysis of a gastroenteritis outbreak to examine the patterns of transmissions, household attack rates and the basic reproduction number (R0) for Norovirus in a child-care facility. We used data from parental interviews of suspected cases sent home with gastroenteritis at a child-care center between 24th of August and 18th of September 2020. A total of 52 persons in 19 household clusters were symptomatic in this outbreak investigation. Of all transmissions, 23 (46.9%) occurred in the child-care center, the rest occurring in households. We found a household attack rate of 36.5% (95% CI 27.3, 47.1%). Serial intervals were estimated as mean 2.5 ± SD1.45 days. The R0, using time-dependent methods during the growth phase of the outbreak (days 2 to 8) was 2.4 (95% CI 1.50, 3.50). The count of affected persons of a child-care center norovirus outbreak is approximately double the count of the total symptomatic staff and attending children. In the study setting, each symptomatic child-care attendee likely infected one other child-care attendee or staff and just over one household contact on average. [ABSTRACT FROM AUTHOR]

We present a detailed set-based analysis of the well-known SIR and SEIR epidemic models subjected to hard caps on the proportion of infective individuals, and bounds on the allowable intervention strategies, such as social distancing, quarantining and vaccination. We describe the admissible and maximal robust positively invariant (MRPI) sets of these two models via the theory of barriers. We show how the sets may be used in the management of epidemics, for both perfect and imperfect/uncertain models, detailing how intervention strategies may be specified such that the hard infection cap is never breached, regardless of the basic reproduction number. The results are clarified with detailed examples. [ABSTRACT FROM AUTHOR]

Due to COVID-19, many countries including Japan have implemented a suspension of economic activities for infection control. It has contributed to reduce the transmission of COVID-19 but caused severe economic losses. Today, several promising vaccines have been developed and are already being distributed in some countries. Therefore, we evaluated various vaccine and intensive countermeasure strategies with constraint of economic loss using SEIR model to obtain knowledge of how to balance economy with infection control in Japan. Our main results were that the vaccination strategy that prioritized younger generation was better in terms of deaths when a linear relationship between lockdown intensity and acceptable economic loss was assumed. On the other hand, when a non-linearity relationship was introduced, implying that the strong lockdown with small economic loss was possible, the old first strategies were best in the settings of small basic reproduction number. These results indicated a high potential of remote work when prioritizing vaccination for the old generation. When focusing on only the old first strategies as the Japanese government has decided to do, the strategy vaccinating the young next to the old was superior to the others when a non-linear relationship was assumed due to sufficient reduction of contact with small economic loss. [ABSTRACT FROM AUTHOR]