Bayesian inference for the onset time and epidemiological characteristics of emerging infectious diseases.

Background: Emerging infectious diseases pose a significant threat to global public health. Timely detection and response are crucial in mitigating the spread of such epidemics. Inferring the onset time and epidemiological characteristics is vital for accelerating early interventions, but accurately predicting these parameters in the early stages remains challenging. Methods: We introduce a Bayesian inference method to fit epidemic models to time series data based on state-space modeling, employing a stochastic Susceptible-Exposed-Infectious-Removed (SEIR) model for transmission dynamics analysis. Our approach uses the particle Markov chain Monte Carlo (PMCMC) method to estimate key epidemiological parameters, including the onset time, the transmission rate, and the recovery rate. The PMCMC algorithm integrates the advantageous aspects of both MCMC and particle filtering methodologies to yield a computationally feasible and effective means of approximating the likelihood function, especially when it is computationally intractable. Results: To validate the proposed method, we conduct case studies on COVID-19 outbreaks in Wuhan, Shanghai and Nanjing, China, respectively. Using early-stage case reports, the PMCMC algorithm accurately predicted the onset time, key epidemiological parameters, and the basic reproduction number. These findings are consistent with empirical studies and the literature. Conclusion: This study presents a robust Bayesian inference method for the timely investigation of emerging infectious diseases. By accurately estimating the onset time and essential epidemiological parameters, our approach is versatile and efficient, extending its utility beyond COVID-19.

Conventional and frugal methods of estimating COVID-19-related excess deaths and undercount factors.

Across the world, the officially reported number of COVID-19 deaths is likely an undercount. Establishing true mortality is key to improving data transparency and strengthening public health systems to tackle future disease outbreaks. In this study, we estimated excess deaths during the COVID-19 pandemic in the Pune region of India. Excess deaths are defined as the number of additional deaths relative to those expected from pre-COVID-19-pandemic trends. We integrated data from: (a) epidemiological modeling using pre-pandemic all-cause mortality data, (b) discrepancies between media-reported death compensation claims and official reported mortality, and (c) the "wisdom of crowds" public surveying. Our results point to an estimated 14,770 excess deaths [95% CI 9820-22,790] in Pune from March 2020 to December 2021, of which 9093 were officially counted as COVID-19 deaths. We further calculated the undercount factor-the ratio of excess deaths to officially reported COVID-19 deaths. Our results point to an estimated undercount factor of 1.6 [95% CI 1.1-2.5]. Besides providing similar conclusions about excess deaths estimates across different methods, our study demonstrates the utility of frugal methods such as the analysis of death compensation claims and the wisdom of crowds in estimating excess mortality.

Future malaria environmental suitability in Africa is sensitive to hydrology.

Changes in climate shift the geographic locations that are suitable for malaria transmission because of the thermal constraints on vector Anopheles mosquitos and Plasmodium spp. malaria parasites and the lack of availability of surface water for vector breeding. Previous Africa-wide assessments have tended to solely represent surface water using precipitation, ignoring many important hydrological processes. Here, we applied a validated and weighted ensemble of global hydrological and climate models to estimate present and future areas of hydroclimatic suitability for malaria transmission. With explicit surface water representation, we predict a net decrease in areas suitable for malaria transmission from 2025 onward, greater sensitivity to future greenhouse gas emissions, and different, more complex, malaria transmission patterns. Areas of malaria transmission that are projected to change are smaller than those estimated by precipitation-based estimates but are associated with greater changes in transmission season lengths.

Modelling the impact of precaution on disease dynamics and its evolution.

In this paper, we introduce the notion of practically susceptible population, which is a fraction of the biologically susceptible population. Assuming that the fraction depends on the severity of the epidemic and the public's level of precaution (as a response of the public to the epidemic), we propose a general framework model with the response level evolving with the epidemic. We firstly verify the well-posedness and confirm the disease's eventual vanishing for the framework model under the assumption that the basic reproduction number R 0 < 1 . For R 0 > 1 , we study how the behavioural response evolves with epidemics and how such an evolution impacts the disease dynamics. More specifically, when the precaution level is taken to be the instantaneous best response function in literature, we show that the endemic dynamic is convergence to the endemic equilibrium; while when the precaution level is the delayed best response, the endemic dynamic can be either convergence to the endemic equilibrium, or convergence to a positive periodic solution. Our derivation offers a justification/explanation for the best response used in some literature. By replacing "adopting the best response" with "adapting toward the best response", we also explore the adaptive long-term dynamics.

K-Track-Covid: interactive web-based dashboard for analyzing geographical and temporal spread of COVID-19 in South Korea.

The COVID-19 pandemic has necessitated the development of robust tools for tracking and modeling the spread of the virus. We present 'K-Track-Covid,' an interactive web-based dashboard developed using the R Shiny framework, to offer users an intuitive dashboard for analyzing the geographical and temporal spread of COVID-19 in South Korea. Our dashboard employs dynamic user interface elements, employs validated epidemiological models, and integrates regional data to offer tailored visual displays. The dashboard allows users to customize their data views by selecting specific time frames, geographic regions, and demographic groups. This customization enables the generation of charts and statistical summaries pertinent to both daily fluctuations and cumulative counts of COVID-19 cases, as well as mortality statistics. Additionally, the dashboard offers a simulation model based on mathematical models, enabling users to make predictions under various parameter settings. The dashboard is designed to assist researchers, policymakers, and the public in understanding the spread and impact of COVID-19, thereby facilitating informed decision-making. All data and resources related to this study are publicly available to ensure transparency and facilitate further research.

Importance of social inequalities to contact patterns, vaccine uptake, and epidemic dynamics.

Individuals' socio-demographic and economic characteristics crucially shape the spread of an epidemic by largely determining the exposure level to the virus and the severity of the disease for those who got infected. While the complex interplay between individual characteristics and epidemic dynamics is widely recognised, traditional mathematical models often overlook these factors. In this study, we examine two important aspects of human behaviour relevant to epidemics: contact patterns and vaccination uptake. Using data collected during the COVID-19 pandemic in Hungary, we first identify the dimensions along which individuals exhibit the greatest variation in their contact patterns and vaccination uptake. We find that generally higher socio-economic groups of the population have a higher number of contacts and a higher vaccination uptake with respect to disadvantaged groups. Subsequently, we propose a data-driven epidemiological model that incorporates these behavioural differences. Finally, we apply our model to analyse the fourth wave of COVID-19 in Hungary, providing valuable insights into real-world scenarios. By bridging the gap between individual characteristics and epidemic spread, our research contributes to a more comprehensive understanding of disease dynamics and informs effective public health strategies.

A Coupled Statistical and Deterministic Model for Forecasting Climate-Driven Dengue Incidence in Selangor, Malaysia.

The mosquito-borne dengue virus remains a major public health concern in Malaysia. Despite various control efforts and measures introduced by the Malaysian Government to combat dengue, the increasing trend of dengue cases persists and shows no sign of decreasing. Currently, early detection and vector control are the main methods employed to curb dengue outbreaks. In this study, a coupled model consisting of the statistical ARIMAX model and the deterministic SI-SIR model was developed and validated using the weekly reported dengue data from year 2014 to 2019 for Selangor, Malaysia. Previous studies have shown that climate variables, especially temperature, humidity, and precipitation, were able to influence dengue incidence and transmission dynamics through their effect on the vector. In this coupled model, climate is linked to dengue disease through mosquito biting rate, allowing real-time forecast of dengue cases using climate variables, namely temperature, rainfall and humidity. For the period chosen for model validation, the coupled model can forecast 1-2 weeks in advance with an average error of less than 6%, three weeks in advance with an average error of 7.06% and four weeks in advance with an average error of 8.01%. Further model simulation analysis suggests that the coupled model generally provides better forecast than the stand-alone ARIMAX model, especially at the onset of the outbreak. Moreover, the coupled model is more robust in the sense that it can be further adapted for investigating the effectiveness of various dengue mitigation measures subject to the changing climate.

A risk-induced dispersal strategy of the infected population for a disease-free state in the SIS epidemic model.

This article proposes a dispersal strategy for infected individuals in a spatial susceptible-infected-susceptible (SIS) epidemic model. The presence of spatial heterogeneity and the movement of individuals play crucial roles in determining the persistence and eradication of infectious diseases. To capture these dynamics, we introduce a moving strategy called risk-induced dispersal (RID) for infected individuals in a continuous-time patch model of the SIS epidemic. First, we establish a continuous-time n-patch model and verify that the RID strategy is an effective approach for attaining a disease-free state. This is substantiated through simulations conducted on 7-patch models and analytical results derived from 2-patch models. Second, we extend our analysis by adapting the patch model into a diffusive epidemic model. This extension allows us to explore further the impact of the RID movement strategy on disease transmission and control. We validate our results through simulations, which provide the effects of the RID dispersal strategy.

Real-Time Use of a Dynamic Model To Measure the Impact of Public Health Interventions on Measles Outbreak Size and Duration - Chicago, Illinois, 2024.

Measles is a highly infectious, vaccine-preventable disease that can cause severe illness, hospitalization, and death. A measles outbreak associated with a migrant shelter in Chicago occurred during February-April 2024, in which a total of 57 confirmed cases were identified, including 52 among shelter residents, three among staff members, and two among community members with a known link to the shelter. CDC simulated a measles outbreak among shelter residents using a dynamic disease model, updated in real time as additional cases were identified, to produce outbreak forecasts and assess the impact of public health interventions. As of April 8, the model forecasted a median final outbreak size of 58 cases (IQR = 56-60 cases); model fit and prediction range improved as more case data became available. Counterfactual analysis of different intervention scenarios demonstrated the importance of early deployment of public health interventions in Chicago, with a 69% chance of an outbreak of 100 or more cases had there been no mass vaccination or active case-finding compared with only a 1% chance when those interventions were deployed. This analysis highlights the value of using real-time, dynamic models to aid public health response, set expectations about outbreak size and duration, and quantify the impact of interventions. The model shows that prompt mass vaccination and active case-finding likely substantially reduced the chance of a large (100 or more cases) outbreak in Chicago.

Static graph approximations of dynamic contact networks for epidemic forecasting.

Epidemic modeling is essential in understanding the spread of infectious diseases like COVID-19 and devising effective intervention strategies to control them. Recently, network-based disease models have integrated traditional compartment-based modeling with real-world contact graphs and shown promising results. However, in an ongoing epidemic, future contact network patterns are not observed yet. To address this, we use aggregated static networks to approximate future contacts for disease modeling. The standard method in the literature concatenates all edges from a dynamic graph into one collapsed graph, called the full static graph. However, the full static graph often leads to severe overestimation of key epidemic characteristics. Therefore, we propose two novel static network approximation methods, DegMST and EdgeMST, designed to preserve the sparsity of real world contact network while remaining connected. DegMST and EdgeMST use the frequency of temporal edges and the node degrees respectively to preserve sparsity. Our analysis show that our models more closely resemble the network characteristics of the dynamic graph compared to the full static ones. Moreover, our analysis on seven real-world contact networks suggests EdgeMST yield more accurate estimations of disease dynamics for epidemic forecasting when compared to the standard full static method.

Markovian Approach for Exploring Competitive Diseases with Heterogeneity-Evidence from COVID-19 and Influenza in China.

Due to the complex interactions between multiple infectious diseases, the spreading of diseases in human bodies can vary when people are exposed to multiple sources of infection at the same time. Typically, there is heterogeneity in individuals' responses to diseases, and the transmission routes of different diseases also vary. Therefore, this paper proposes an SIS disease spreading model with individual heterogeneity and transmission route heterogeneity under the simultaneous action of two competitive infectious diseases. We derive the theoretical epidemic spreading threshold using quenched mean-field theory and perform numerical analysis under the Markovian method. Numerical results confirm the reliability of the theoretical threshold and show the inhibitory effect of the proportion of fully competitive individuals on epidemic spreading. The results also show that the diversity of disease transmission routes promotes disease spreading, and this effect gradually weakens when the epidemic spreading rate is high enough. Finally, we find a negative correlation between the theoretical spreading threshold and the average degree of the network. We demonstrate the practical application of the model by comparing simulation outputs to temporal trends of two competitive infectious diseases, COVID-19 and seasonal influenza in China.

Threshold dynamics of a reaction-advection-diffusion schistosomiasis epidemic model with seasonality and spatial heterogeneity.

Most water-borne disease models ignore the advection of water flows in order to simplify the mathematical analysis and numerical computation. However, advection can play an important role in determining the disease transmission dynamics. In this paper, we investigate the long-term dynamics of a periodic reaction-advection-diffusion schistosomiasis model and explore the joint impact of advection, seasonality and spatial heterogeneity on the transmission of the disease. We derive the basic reproduction number R 0 and show that the disease-free periodic solution is globally attractive when R 0 < 1 whereas there is a positive endemic periodic solution and the system is uniformly persistent in a special case when R 0 > 1 . Moreover, we find that R 0 is a decreasing function of the advection coefficients which offers insights into why schistosomiasis is more serious in regions with slow water flows.

Transmission dynamics of a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods.

In this paper, we propose a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods. Firstly, we establish the well-posedness of the model. Secondly, we define the basic reproduction number ℜ 0 for this model and show that ℜ 0 is a threshold parameter: if ℜ 0 < 1 , then the disease-free periodic solution is globally attractive; if ℜ 0 > 1 , the system is uniformly persistent. Thirdly, we study the global attractivity of the positive steady state when the spatial environment is homogeneous and the advection of mosquitoes is ignored. As an example, we use the model to investigate the dengue fever transmission case in Guangdong Province, China, and explore the impact of model parameters on ℜ 0 . Our findings indicate that ignoring seasonality may underestimate ℜ 0 . Additionally, the spatial heterogeneity of transmission may increase the risk of disease transmission, while the increase of seasonal developmental durations, intrinsic incubation periods and advection rates can all reduce the risk of disease transmission.

Impact of immune evasion, waning and boosting on dynamics of population mixing between a vaccinated majority and unvaccinated minority.

BACKGROUND: We previously demonstrated that when vaccines prevent infection, the dynamics of mixing between vaccinated and unvaccinated sub-populations is such that use of imperfect vaccines markedly decreases risk for vaccinated people, and for the population overall. Risks to vaccinated people accrue disproportionately from contact with unvaccinated people. In the context of the emergence of Omicron SARS-CoV-2 and evolving understanding of SARS-CoV-2 epidemiology, we updated our analysis to evaluate whether our earlier conclusions remained valid. METHODS: We modified a previously published Susceptible-Infectious-Recovered (SIR) compartmental model of SARS-CoV-2 with two connected sub-populations: vaccinated and unvaccinated, with non-random mixing between groups. Our expanded model incorporates diminished vaccine efficacy for preventing infection with the emergence of Omicron SARS-CoV-2 variants, waning immunity, the impact of prior immune experience on infectivity, "hybrid" effects of infection in previously vaccinated individuals, and booster vaccination. We evaluated the dynamics of an epidemic within each subgroup and in the overall population over a 10-year time horizon. RESULTS: Even with vaccine efficacy as low as 20%, and in the presence of waning immunity, the incidence of COVID-19 in the vaccinated subpopulation was lower than that among the unvaccinated population across the full 10-year time horizon. The cumulative risk of infection was 3-4 fold higher among unvaccinated people than among vaccinated people, and unvaccinated people contributed to infection risk among vaccinated individuals at twice the rate that would have been expected based on the frequency of contacts. These findings were robust across a range of assumptions around the rate of waning immunity, the impact of "hybrid immunity", frequency of boosting, and the impact of prior infection on infectivity in unvaccinated people. INTERPRETATION: Although the emergence of the Omicron variants of SARS-CoV-2 has diminished the protective effects of vaccination against infection with SARS-CoV-2, updating our earlier model to incorporate loss of immunity, diminished vaccine efficacy and a longer time horizon, does not qualitatively change our earlier conclusions. Vaccination against SARS-CoV-2 continues to diminish the risk of infection among vaccinated people and in the population as a whole. By contrast, the risk of infection among vaccinated people accrues disproportionately from contact with unvaccinated people.

A New Auto-Regressive Multi-Variable Modified Auto-Encoder for Multivariate Time-Series Prediction: A Case Study with Application to COVID-19 Pandemics.

The SARS-CoV-2 global pandemic prompted governments, institutions, and researchers to investigate its impact, developing strategies based on general indicators to make the most precise predictions possible. Approaches based on epidemiological models were used but the outcomes demonstrated forecasting with uncertainty due to insufficient or missing data. Besides the lack of data, machine-learning models including random forest, support vector regression, LSTM, Auto-encoders, and traditional time-series models such as Prophet and ARIMA were employed in the task, achieving remarkable results with limited effectiveness. Some of these methodologies have precision constraints in dealing with multi-variable inputs, which are important for problems like pandemics that require short and long-term forecasting. Given the under-supply in this scenario, we propose a novel approach for time-series prediction based on stacking auto-encoder structures using three variations of the same model for the training step and weight adjustment to evaluate its forecasting performance. We conducted comparison experiments with previously published data on COVID-19 cases, deaths, temperature, humidity, and air quality index (AQI) in São Paulo City, Brazil. Additionally, we used the percentage of COVID-19 cases from the top ten affected countries worldwide until May 4th, 2020. The results show 80.7% and 10.3% decrease in RMSE to entire and test data over the distribution of 50 trial-trained models, respectively, compared to the first experiment comparison. Also, model type#3 achieved 4th better overall ranking performance, overcoming the NBEATS, Prophet, and Glounts time-series models in the second experiment comparison. This model shows promising forecast capacity and versatility across different input dataset lengths, making it a prominent forecasting model for time-series tasks.

Informing policy via dynamic models: Cholera in Haiti.

Public health decisions must be made about when and how to implement interventions to control an infectious disease epidemic. These decisions should be informed by data on the epidemic as well as current understanding about the transmission dynamics. Such decisions can be posed as statistical questions about scientifically motivated dynamic models. Thus, we encounter the methodological task of building credible, data-informed decisions based on stochastic, partially observed, nonlinear dynamic models. This necessitates addressing the tradeoff between biological fidelity and model simplicity, and the reality of misspecification for models at all levels of complexity. We assess current methodological approaches to these issues via a case study of the 2010-2019 cholera epidemic in Haiti. We consider three dynamic models developed by expert teams to advise on vaccination policies. We evaluate previous methods used for fitting these models, and we demonstrate modified data analysis strategies leading to improved statistical fit. Specifically, we present approaches for diagnosing model misspecification and the consequent development of improved models. Additionally, we demonstrate the utility of recent advances in likelihood maximization for high-dimensional nonlinear dynamic models, enabling likelihood-based inference for spatiotemporal incidence data using this class of models. Our workflow is reproducible and extendable, facilitating future investigations of this disease system.

Long-term effects of non-pharmaceutical interventions on total disease burden in parsimonious epidemiological models.

The recent global COVID-19 pandemic resulted in governments enacting non-pharmaceutical interventions (NPIs) targeted at reducing transmission of SARS-CoV-2. But the NPIs also affected the transmission of viruses causing non-target seasonal respiratory diseases, including influenza and respiratory syncytial virus (RSV). In many countries, the NPIs were found to reduce cases of such seasonal respiratory diseases, but there is also evidence that subsequent relaxation of NPIs led to outbreaks of these diseases that were larger than pre-pandemic ones, due to the accumulation of susceptible individuals prior to relaxation. Therefore, the net long-term effects of NPIs on the total disease burden of non-target diseases remain unclear. Knowledge of this is important for infectious disease management and maintenance of public health. In this study, we shed light on this issue for the simplified scenario of a set of NPIs that prevent or reduce transmission of a seasonal respiratory disease for about a year and are then removed, using mathematical analyses and numerical simulations of a suite of four epidemiological models with varying complexity and generality. The model parameters were estimated using empirical data pertaining to seasonal respiratory diseases and covered a wide range. Our results showed that NPIs reduced the total disease burden of a non-target seasonal respiratory disease in the long-term. Expressed as a percentage of population size, the reduction was greater for larger values of the basic reproduction number and the immunity loss rate, reflecting larger outbreaks and hence more infections averted by imposition of NPIs. Our study provides a foundation for exploring the effects of NPIs on total disease burden in more-complex scenarios.

Advanced forecasting of COVID-19 epidemic: Leveraging ensemble models, advanced optimization, and decomposition techniques.

In the global effort to address the outbreak of the new coronavirus pneumonia (COVID-19) pandemic, accurate forecasting of epidemic patterns has become crucial for implementing successful interventions aimed at preventing and controlling the spread of the disease. The correct prediction of the course of COVID-19 outbreaks is a complex and challenging task, mainly because of the significant volatility in the data series related to COVID-19. Previous studies have been limited by the exclusive use of individual forecasting techniques in epidemic modeling, disregarding the integration of diverse prediction procedures. The lack of attention to detail in this situation can yield worse-than-ideal results. Consequently, this study introduces a novel ensemble framework that integrates three machine learning methods (kernel ridge regression (KRidge), Deep random vector functional link (dRVFL), and ridge regression) within a linear relationship (L-KRidge-dRVFL-Ridge). The optimization of this framework is accomplished through a distinctive approach, specifically adaptive differential evolution and particle swarm optimization (A-DEPSO). Moreover, an effective decomposition method, known as time-varying filter empirical mode decomposition (TVF-EMD), is employed to decompose the input variables. A feature selection technique, specifically using the light gradient boosting machine (LGBM), is also implemented to extract the most influential input variables. The daily datasets of COVID-19 collected from two countries, namely Italy and Poland, were used as the experimental examples. Additionally, all models are implemented to forecast COVID-19 at two-time horizons: 10- and 14-day ahead (t+10 and t+14). According to the results, the proposed model can yield higher correlation coefficient (R) for both case studies: Italy (t+10 = 0.965, t+14 = 0.961) and Poland (t+10 = 0.952, t+14 = 0.940) than the other models. The experimental results demonstrate that the model suggested in this paper has outstanding results in various kinds of complex epidemic prediction situations. The proposed ensemble model demonstrates exceptional accuracy and resilience, outperforming all similar models in terms of efficacy.

An immuno-epidemiological model with waning immunity after infection or vaccination.

In epidemics, waning immunity is common after infection or vaccination of individuals. Immunity levels are highly heterogeneous and dynamic. This work presents an immuno-epidemiological model that captures the fundamental dynamic features of immunity acquisition and wane after infection or vaccination and analyzes mathematically its dynamical properties. The model consists of a system of first order partial differential equations, involving nonlinear integral terms and different transfer velocities. Structurally, the equation may be interpreted as a Fokker-Planck equation for a piecewise deterministic process. However, unlike the usual models, our equation involves nonlocal effects, representing the infectivity of the whole environment. This, together with the presence of different transfer velocities, makes the proved existence of a solution novel and nontrivial. In addition, the asymptotic behavior of the model is analyzed based on the obtained qualitative properties of the solution. An optimal control problem with objective function including the total number of deaths and costs of vaccination is explored. Numerical results describe the dynamic relationship between contact rates and optimal solutions. The approach can contribute to the understanding of the dynamics of immune responses at population level and may guide public health policies.

Predicting Drug-Drug Interactions Involving Rifampicin Using a Semi-mechanistic Hepatic Compartmental Model.

AIMS: To develop a semi-mechanistic hepatic compartmental model to predict the effects of rifampicin, a known inducer of CYP3A4 enzyme, on the metabolism of five drugs, in the hope of informing dose adjustments to avoid potential drug-drug interactions. METHODS: A search was conducted for DDI studies on the interactions between rifampicin and CYP substrates that met specific criteria, including the availability of plasma concentration-time profiles, physical and absorption parameters, pharmacokinetic parameters, and the use of healthy subjects at therapeutic doses. The semi-mechanistic model utilized in this study was improved from its predecessors, incorporating additional parameters such as population data (specifically for Chinese and Caucasians), virtual individuals, gender distribution, age range, dosing time points, and coefficients of variation. RESULTS: Optimal parameters were identified for our semi-mechanistic model by validating it with clinical data, resulting in a maximum difference of approximately 2-fold between simulated and observed values. PK data of healthy subjects were used for most CYP3A4 substrates, except for gilteritinib, which showed no significant difference between patients and healthy subjects. Dose adjustment of gilteritinib co-administered with rifampicin required a 3-fold increase of the initial dose, while other substrates were further tuned to achieve the desired drug exposure. CONCLUSIONS: The pharmacokinetic parameters AUCR and CmaxR of drugs metabolized by CYP3A4, when influenced by Rifampicin, were predicted by the semi-mechanistic model to be approximately twice the empirically observed values, which suggests that the semi-mechanistic model was able to reasonably simulate the effect. The doses of four drugs adjusted via simulation to reduce rifampicin interaction.