Early detection of emerging viral variants through analysis of community structure of coordinated substitution networks.

The emergence of viral variants with altered phenotypes is a public health challenge underscoring the need for advanced evolutionary forecasting methods. Given extensive epistatic interactions within viral genomes and known viral evolutionary history, efficient genomic surveillance necessitates early detection of emerging viral haplotypes rather than commonly targeted single mutations. Haplotype inference, however, is a significantly more challenging problem precluding the use of traditional approaches. Here, using SARS-CoV-2 evolutionary dynamics as a case study, we show that emerging haplotypes with altered transmissibility can be linked to dense communities in coordinated substitution networks, which become discernible significantly earlier than the haplotypes become prevalent. From these insights, we develop a computational framework for inference of viral variants and validate it by successful early detection of known SARS-CoV-2 strains. Our methodology offers greater scalability than phylogenetic lineage tracing and can be applied to any rapidly evolving pathogen with adequate genomic surveillance data.

SARS-CoV-2 Transmission in Alberta, British Columbia, and Ontario, Canada, January 2020-January 2022.

We estimated COVID-19 transmission potential and case burden by variant type in Alberta, British Columbia, and Ontario, Canada, during January 23, 2020-January 27, 2022; we also estimated the effectiveness of public health interventions to reduce transmission. We estimated time-varying reproduction number (Rt) over 7-day sliding windows and nonoverlapping time-windows determined by timing of policy changes. We calculated incidence rate ratios (IRRs) for each variant and compared rates to determine differences in burden among provinces. Rt corresponding with emergence of the Delta variant increased in all 3 provinces; British Columbia had the largest increase, 43.85% (95% credible interval [CrI] 40.71%-46.84%). Across the study period, IRR was highest for Omicron (8.74 [95% CrI 8.71-8.77]) and burden highest in Alberta (IRR 1.80 [95% CrI 1.79-1.81]). Initiating public health interventions was associated with lower Rt and relaxing restrictions and emergence of new variants associated with increases in Rt.

Parameter estimation and forecasting with quantified uncertainty for ordinary differential equation models using QuantDiffForecast: A MATLAB toolbox and tutorial.

Mathematical models based on systems of ordinary differential equations (ODEs) are frequently applied in various scientific fields to assess hypotheses, estimate key model parameters, and generate predictions about the system's state. To support their application, we present a comprehensive, easy-to-use, and flexible MATLAB toolbox, QuantDiffForecast, and associated tutorial to estimate parameters and generate short-term forecasts with quantified uncertainty from dynamical models based on systems of ODEs. We provide software ( https://github.com/gchowell/paramEstimation_forecasting_ODEmodels/) and detailed guidance on estimating parameters and forecasting time-series trajectories that are characterized using ODEs with quantified uncertainty through a parametric bootstrapping approach. It includes functions that allow the user to infer model parameters and assess forecasting performance for different ODE models specified by the user, using different estimation methods and error structures in the data. The tutorial is intended for a diverse audience, including students training in dynamic systems, and will be broadly applicable to estimate parameters and generate forecasts from models based on ODEs. The functions included in the toolbox are illustrated using epidemic models with varying levels of complexity applied to data from the 1918 influenza pandemic in San Francisco. A tutorial video that demonstrates the functionality of the toolbox is included.

Modeling the health impact of increasing vaccine coverage and nonpharmaceutical interventions against coronavirus disease 2019 in Ghana.

Seroprevalence studies assessing community exposure to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Ghana concluded that population-level immunity remained low as of February 2021. Thus, it is important to demonstrate how increasing vaccine coverage reduces the economic and public health impacts associated with SARS-CoV-2 transmission. To that end, this study used a Susceptible-Exposed-Presymptomatic-Symptomatic-Asymptomatic-Recovered-Dead-Vaccinated compartmental model to simulate coronavirus disease 2019 (COVID-19) transmission and the role of public health interventions in Ghana. The impact of increasing vaccination rates and decline in transmission rates due to nonpharmaceutical interventions (NPIs) on cumulative infections and deaths averted was explored under different scenarios. Latin hypercube sampling-partial rank correlation coefficient (LHS-PRCC) was used to investigate the uncertainty and sensitivity of the outcomes to the parameters. Simulation results suggest that increasing the vaccination rate to achieve 50% coverage was associated with almost 60,000 deaths and 25 million infections averted. In comparison, a 50% decrease in the transmission coefficient was associated with the prevention of about 150,000 deaths and 50 million infections. The LHS-PRCC results indicated that in the context of vaccination rate, cumulative infections and deaths averted were most sensitive to vaccination rate, waning immunity rates from vaccination, and waning immunity from natural infection. This study's findings illustrate the impact of increasing vaccination coverage and/or reducing the transmission rate by NPI adherence in the prevention of COVID-19 infections and deaths in Ghana.

SubEpiPredict: A tutorial-based primer and toolbox for fitting and forecasting growth trajectories using the ensemble n-sub-epidemic modeling framework.

An ensemble n-sub-epidemic modeling framework that integrates sub-epidemics to capture complex temporal dynamics has demonstrated powerful forecasting capability in previous works. This modeling framework can characterize complex epidemic patterns, including plateaus, epidemic resurgences, and epidemic waves characterized by multiple peaks of different sizes. In this tutorial paper, we introduce and illustrate SubEpiPredict, a user-friendly MATLAB toolbox for fitting and forecasting time series data using an ensemble n-sub-epidemic modeling framework. The toolbox can be used for model fitting, forecasting, and evaluation of model performance of the calibration and forecasting periods using metrics such as the weighted interval score (WIS). We also provide a detailed description of these methods including the concept of the n-sub-epidemic model, constructing ensemble forecasts from the top-ranking models, etc. For the illustration of the toolbox, we utilize publicly available daily COVID-19 death data at the national level for the United States. The MATLAB toolbox introduced in this paper can be very useful for a wider group of audiences, including policymakers, and can be easily utilized by those without extensive coding and modeling backgrounds.

GrowthPredict: A toolbox and tutorial-based primer for fitting and forecasting growth trajectories using phenomenological growth models.

Simple dynamic modeling tools can help generate real-time short-term forecasts with quantified uncertainty of the trajectory of diverse growth processes unfolding in nature and society, including disease outbreaks. An easy-to-use and flexible toolbox for this purpose is lacking. This tutorial-based primer introduces and illustrates GrowthPredict, a user-friendly MATLAB toolbox for fitting and forecasting time-series trajectories using phenomenological dynamic growth models based on ordinary differential equations. This toolbox is accessible to a broad audience, including students training in mathematical biology, applied statistics, and infectious disease modeling, as well as researchers and policymakers who need to conduct short-term forecasts in real-time. The models included in the toolbox capture exponential and sub-exponential growth patterns that typically follow a rising pattern followed by a decline phase, a common feature of contagion processes. Models include the 1-parameter exponential growth model and the 2-parameter generalized-growth model, which have proven useful in characterizing and forecasting the ascending phase of epidemic outbreaks. It also includes the 2-parameter Gompertz model, the 3-parameter generalized logistic-growth model, and the 3-parameter Richards model, which have demonstrated competitive performance in forecasting single peak outbreaks. We provide detailed guidance on forecasting time-series trajectories and available software ( https://github.com/gchowell/forecasting_growthmodels ), including the full uncertainty distribution derived through parametric bootstrapping, which is needed to construct prediction intervals and evaluate their accuracy. Functions are available to assess forecasting performance across different models, estimation methods, error structures in the data, and forecasting horizons. The toolbox also includes functions to quantify forecasting performance using metrics that evaluate point and distributional forecasts, including the weighted interval score. This tutorial and toolbox can be broadly applied to characterizing and forecasting time-series data using simple phenomenological growth models. As a contagion process takes off, the tools presented in this tutorial can help create forecasts to guide policy regarding implementing control strategies and assess the impact of interventions. The toolbox functionality is demonstrated through various examples, including a tutorial video, and the examples use publicly available data on the monkeypox (mpox) epidemic in the USA.

Impact of public health policy and mobility change on transmission potential of severe acute respiratory syndrome coronavirus 2 in Rhode Island, March 2020 - November 2021.

To study the SARS-CoV-2 transmission potential in Rhode Island (RI) and its association with policy changes and mobility changes, the time-varying reproduction number, Rt, was estimated. The daily incident case counts (16 March 2020, through 30 November 2021) were bootstrapped within a 15-day sliding window and multiplied by Poisson-distributed multipliers (λ = 4, sensitivity analysis: 11) to generate 1000 estimated infection counts, to which EpiEstim was applied to generate Rt time series. The median Rt percentage change when policies changed was estimated. The time lag correlations were assessed between the 7-day moving average of the relative changes in Google mobility data in the first 90 days, and Rt and estimated infection count, respectively. There were three major pandemic waves in RI in 2020-2021: spring 2020, winter 2020-2021 and fall-winter 2021. The median Rt fluctuated within the range of 0.5-2 from April 2020 to November 2021. Mask mandate (18 April 2020) was associated with a decrease in Rt (-25.99%, 95% CrI: -37.42%, -14.30%). Termination of mask mandates on 6 July 2021 was associated with an increase in Rt (36.74%, 95% CrI: 27.20%, 49.13%). Positive correlations were found between changes in grocery and pharmacy, Rt retail and recreation, transit, and workplace visits, for both Rt and estimated infection count, respectively. Negative correlations were found between changes in residential area visits for both Rt and estimated infection count, respectively. Public health policies enacted in RI were associated with changes in the pandemic trajectory. This ecological study provides further evidence of how non-pharmaceutical interventions and vaccination slowed COVID-19 transmission in RI.

Epidemiological impact of travel enhancement on the inter-prefectural importation dynamics of COVID-19 in Japan, 2020.

Mobility restrictions were widely practiced to reduce contact with others and prevent the spatial spread of COVID-19 infection. Using inter-prefectural mobility and epidemiological data, a statistical model was devised to predict the number of imported cases in each Japanese prefecture. The number of imported cases crossing prefectural borders in 2020 was predicted using inter-prefectural mobility rates based on mobile phone data and prevalence estimates in the origin prefectures. The simplistic model was quantified using surveillance data of cases with an inter-prefectural travel history. Subsequently, simulations were carried out to understand how imported cases vary with the mobility rate and prevalence at the origin. Overall, the predicted number of imported cases qualitatively captured the observed number of imported cases over time. Although Hokkaido and Okinawa are the northernmost and the southernmost prefectures, respectively, they were sensitive to differing prevalence rate in Tokyo and Osaka and the mobility rate. Additionally, other prefectures were sensitive to mobility change, assuming that an increment in the mobility rate was seen in all prefectures. Our findings indicate the need to account for the weight of an inter-prefectural mobility network when implementing countermeasures to restrict human movement. If the mobility rates were maintained lower than the observed rates, then the number of imported cases could have been maintained at substantially lower levels than the observed, thus potentially preventing the unnecessary spatial spread of COVID-19 in late 2020.

Structural identifiability analysis of epidemic models based on differential equations: a tutorial-based primer.

The successful application of epidemic models hinges on our ability to estimate model parameters from limited observations reliably. An often-overlooked step before estimating model parameters consists of ensuring that the model parameters are structurally identifiable from the observed states of the system. In this tutorial-based primer, intended for a diverse audience, including students training in dynamic systems, we review and provide detailed guidance for conducting structural identifiability analysis of differential equation epidemic models based on a differential algebra approach using differential algebra for identifiability of systems (DAISY) and Mathematica (Wolfram Research). This approach aims to uncover any existing parameter correlations that preclude their estimation from the observed variables. We demonstrate this approach through examples, including tutorial videos of compartmental epidemic models previously employed to study transmission dynamics and control. We show that the lack of structural identifiability may be remedied by incorporating additional observations from different model states, assuming that the system's initial conditions are known, using prior information to fix some parameters involved in parameter correlations, or modifying the model based on existing parameter correlations. We also underscore how the results of structural identifiability analysis can help enrich compartmental diagrams of differential-equation models by indicating the observed state variables and the results of the structural identifiability analysis.

Modeling the Impact of Agricultural Mitigation Measures on the Spread of Sharka Disease in Sweet Cherry Orchards.

Sharka is a disease affecting stone fruit trees. It is caused by the Plum pox virus (PPV), with Myzus persicae being one of the most efficient aphid species in transmitting it within and among Prunus orchards. Other agricultural management strategies are also responsible for the spread of disease among trees, such as grafting and pruning. We present a mathematical model of impulsive differential equations to represent the dynamics of Sharka disease in the tree and vector population. We consider three transmission routes: grafting, pruning, and through aphid vectors. Grafting, pruning, and vector control occur as pulses at specific instants. Within the model, human risk perception towards disease influences these agricultural management strategies. Model results show that grafting with infected biological material has a significant impact on the spread of the disease. In addition, detecting infectious symptomatic and asymptomatic trees in the short term is critical to reduce disease spread. Furthermore, vector control to prevent aphid movement between trees is crucial for disease mitigation, as well as implementing awareness campaigns for Sharka disease in agricultural communities that provide a long-term impact on responsible pruning, grafting, and vector control.

Retrospective evaluation of short-term forecast performance of ensemble sub-epidemic frameworks and other time-series models: The 2022-2023 mpox outbreak across multiple geographical scales, July 14th, 2022, through February 26th, 2023.

In May 2022, public health officials noted an unprecedented surge in mpox cases in non-endemic countries worldwide. As the epidemic accelerated, multi-model forecasts of the epidemic's trajectory were critical in guiding the implementation of public health interventions and determining policy. As the case levels have significantly decreased as of early September 2022, evaluating model performance is essential to advance the growing field of epidemic forecasting. Using laboratory-confirmed mpox case data from the Centers for Disease Control and Prevention (CDC) and Our World in Data (OWID) teams through the week of January 26th, 2023, we generated retrospective sequential weekly forecasts (e.g., 1-4-weeks) for Brazil, Canada, France, Germany, Spain, the United Kingdom, the USA, and at the global scale using models that require minimal input data including the auto-regressive integrated moving average (ARIMA), general additive model (GAM), simple linear regression (SLR), Facebook's Prophet model, as well as the sub-epidemic wave (spatial-wave) and n-sub-epidemic modeling frameworks. We assess forecast performance using average mean squared error (MSE), mean absolute error (MAE), weighted interval score (WIS), 95% prediction interval coverage (95% PI coverage), and skill scores. Average Winkler scores were used to calculate skill scores for 95% PI coverage. Overall, the n-sub-epidemic modeling framework outcompeted other models across most locations and forecasting horizons, with the unweighted ensemble model performing best across all forecasting horizons for most locations regarding average MSE, MAE, WIS, and 95% PI coverage. However, many locations had multiple models performing equally well for the average 95% PI coverage. The n-sub-epidemic and spatial-wave frameworks improved considerably in average MSE, MAE, and WIS, and Winkler scores (95% PI coverage) relative to the ARIMA model. Findings lend further support to sub-epidemic frameworks for short-term forecasting epidemics of emerging and re-emerging infectious diseases.

A computational approach to identifiability analysis for a model of the propagation and control of COVID-19 in Chile.

A computational approach is adapted to analyze the parameter identifiability of a compartmental model. The model is intended to describe the progression of the COVID-19 pandemic in Chile during the initial phase in early 2020 when government declared quarantine measures. The computational approach to analyze the structural and practical identifiability is applied in two parts, one for synthetic data and another for some Chilean regional data. The first part defines the identifiable parameter sets when these recover the true parameters used to create the synthetic data. The second part compares the results derived from synthetic data, estimating the identifiable parameter sets from regional Chilean epidemic data. Experiments provide evidence of the loss of identifiability if some initial conditions are estimated, the period of time used to fit is before the peak, and if a significant proportion of the population is involved in quarantine periods.

Evaluation of Bayesian spatiotemporal infectious disease models for prospective surveillance analysis.

BACKGROUND: COVID-19 brought enormous challenges to public health surveillance and underscored the importance of developing and maintaining robust systems for accurate surveillance. As public health data collection efforts expand, there is a critical need for infectious disease modeling researchers to continue to develop prospective surveillance metrics and statistical models to accommodate the modeling of large disease counts and variability. This paper evaluated different likelihoods for the disease count model and various spatiotemporal mean models for prospective surveillance. METHODS: We evaluated Bayesian spatiotemporal models, which are the foundation for model-based infectious disease surveillance metrics. Bayesian spatiotemporal mean models based on the Poisson and the negative binomial likelihoods were evaluated with the different lengths of past data usage. We compared their goodness of fit and short-term prediction performance with both simulated epidemic data and real data from the COVID-19 pandemic. RESULTS: The simulation results show that the negative binomial likelihood-based models show better goodness of fit results than Poisson likelihood-based models as deemed by smaller deviance information criteria (DIC) values. However, Poisson models yield smaller mean square error (MSE) and mean absolute one-step prediction error (MAOSPE) results when we use a shorter length of the past data such as 7 and 3 time periods. Real COVID-19 data analysis of New Jersey and South Carolina shows similar results for the goodness of fit and short-term prediction results. Negative binomial-based mean models showed better performance when we used the past data of 52 time periods. Poisson-based mean models showed comparable goodness of fit performance and smaller MSE and MAOSPE results when we used the past data of 7 and 3 time periods. CONCLUSION: We evaluate these models and provide future infectious disease outbreak modeling guidelines for Bayesian spatiotemporal analysis. Our choice of the likelihood and spatiotemporal mean models was influenced by both historical data length and variability. With a longer length of past data usage and more over-dispersed data, the negative binomial likelihood shows a better model fit than the Poisson likelihood. However, as we use a shorter length of the past data for our surveillance analysis, the difference between the Poisson and the negative binomial models becomes smaller. In this case, the Poisson likelihood shows robust posterior mean estimate and short-term prediction results.

A MATLAB toolbox to fit and forecast growth trajectories using phenomenological growth models: Application to epidemic outbreaks.

Background: Simple dynamic modeling tools can be useful for generating real-time short-term forecasts with quantified uncertainty of the trajectory of diverse growth processes unfolding in nature and society, including disease outbreaks. An easy-to-use and flexible toolbox for this purpose is lacking. Results: In this tutorial-based primer, we introduce and illustrate a user-friendly MATLAB toolbox for fitting and forecasting time-series trajectories using phenomenological dynamic growth models based on ordinary differential equations. This toolbox is accessible to various audiences, including students training in time-series forecasting, dynamic growth modeling, parameter estimation, parameter uncertainty and identifiability, model comparison, performance metrics, and forecast evaluation, as well as researchers and policymakers who need to conduct short-term forecasts in real-time. The models included in the toolbox capture exponential and sub-exponential growth patterns that typically follow a rising pattern followed by a decline phase, a common feature of contagion processes. Models include the 2-parameter generalized-growth model, which has proved useful to characterize and forecast the ascending phase of epidemic outbreaks, and the Gompertz model as well as the 3-parameter generalized logistic-growth model and the Richards model, which have demonstrated competitive performance in forecasting single peak outbreaks.The toolbox provides a tutorial for forecasting time-series trajectories that include the full uncertainty distribution, derived through parametric bootstrapping, which is needed to construct prediction intervals and evaluate their accuracy. Functions are available to assess forecasting performance across different models, estimation methods, error structures in the data, and forecasting horizons. The toolbox also includes functions to quantify forecasting performance using metrics that evaluate point and distributional forecasts, including the weighted interval score. Conclusions: We have developed the first comprehensive toolbox to characterize and forecast time-series data using simple phenomenological growth models. As a contagion process takes off, the tools presented in this tutorial can facilitate policymaking to guide the implementation of control strategies and assess the impact of interventions. The toolbox functionality is demonstrated through various examples, including a tutorial video, and is illustrated using weekly data on the monkeypox epidemic in the USA.

Excess mortality in Ukraine during the course of COVID-19 pandemic in 2020-2021.

In this work, the COVID-19 pandemic burden in Ukraine is investigated retrospectively using the excess mortality measures during 2020-2021. In particular, the epidemic impact on the Ukrainian population is studied via the standardized both all-cause and cause-specific mortality scores before and during the epidemic. The excess mortality counts during the pandemic were predicted based on historic data using parametric and nonparametric modeling and then compared with the actual reported counts to quantify the excess. The corresponding standardized mortality P-score metrics were also compared with the neighboring countries. In summary, there were three "waves" of excess all-cause mortality in Ukraine in December 2020, April 2021 and November 2021 with excess of 32%, 43% and 83% above the expected mortality. Each new "wave" of the all-cause mortality was higher than the previous one and the mortality "peaks" corresponded in time to three "waves" of lab-confirmed COVID-19 mortality. The lab-confirmed COVID-19 mortality constituted 9% to 24% of the all-cause mortality during those three peak months. Overall, the mortality trends in Ukraine over time were similar to neighboring countries where vaccination coverage was similar to that in Ukraine. For cause-specific mortality, the excess observed was due to pneumonia as well as circulatory system disease categories that peaked at the same times as the all-cause and lab-confirmed COVID-19 mortality, which was expected. The pneumonias as well as circulatory system disease categories constituted the majority of all cases during those peak times. The seasonality in mortality due to the infectious and parasitic disease category became less pronounced during the pandemic. While the reported numbers were always relatively low, alcohol-related mortality also declined during the pandemic.

Age-specific mortality and the role of living remotely: The 1918-20 influenza pandemic in Kautokeino and Karasjok, Norway.

The 1918-20 pandemic influenza killed 50-100 million people worldwide, but mortality varied by ethnicity and geography. In Norway, areas dominated by Sámi experienced 3-5 times higher mortality than the country's average. We here use data from burial registers and censuses to calculate all-cause excess mortality by age and wave in two remote Sámi areas of Norway 1918-20. We hypothesise that geographic isolation, less prior exposure to seasonal influenza, and thus less immunity led to higher Indigenous mortality and a different age distribution of mortality (higher mortality for all) than was typical for this pandemic in non-isolated majority populations (higher young adult mortality & sparing of the elderly). Our results show that in the fall of 1918 (Karasjok), winter of 1919 (Kautokeino), and winter of 1920 (Karasjok), young adults had the highest excess mortality, followed by also high excess mortality among the elderly and children. Children did not exhibit excess mortality in the second wave in Karasjok in 1920. It was not the young adults alone who produced the excess mortality in Kautokeino and Karasjok. We conclude that geographic isolation caused higher mortality among the elderly in the first and second waves, and among children in the first wave.

Exploring Social Media Network Connections to Assist During Public Health Emergency Response: A Retrospective Case-Study of Hurricane Matthew and Twitter Users in Georgia, USA.

OBJECTIVE: To assist communities who suffered from hurricane-inflicted damages, emergency responders may monitor social media messages. We present a case-study using the event of Hurricane Matthew to analyze the results of an imputation method for the location of Twitter users who follow school and school districts in Georgia, USA. METHODS: Tweets related to Hurricane Matthew were analyzed by content analysis with latent Dirichlet allocation models and sentiment analysis to identify needs and sentiment changes over time. A hurdle regression model was applied to study the association between retweet frequency and content analysis topics. RESULTS: Users residing in counties affected by Hurricane Matthew posted tweets related to preparedness (n = 171; 16%), awareness (n = 407; 38%), call-for-action or help (n = 206; 19%), and evacuations (n = 93; 9%), with mostly a negative sentiment during the preparedness and response phase. Tweets posted in the hurricane path during the preparedness and response phase were less likely to be retweeted than those outside the path (adjusted odds ratio: 0.95; 95% confidence interval: 0.75, 1.19). CONCLUSIONS: Social media data can be used to detect and evaluate damages of communities affected by natural disasters and identify users' needs in at-risk areas before the event takes place to aid during the preparedness phases.

Real-time forecasting the trajectory of monkeypox outbreaks at the national and global levels, July-October 2022.

BACKGROUND: Beginning May 7, 2022, multiple nations reported an unprecedented surge in monkeypox cases. Unlike past outbreaks, differences in affected populations, transmission mode, and clinical characteristics have been noted. With the existing uncertainties of the outbreak, real-time short-term forecasting can guide and evaluate the effectiveness of public health measures. METHODS: We obtained publicly available data on confirmed weekly cases of monkeypox at the global level and for seven countries (with the highest burden of disease at the time this study was initiated) from the Our World in Data (OWID) GitHub repository and CDC website. We generated short-term forecasts of new cases of monkeypox across the study areas using an ensemble n-sub-epidemic modeling framework based on weekly cases using 10-week calibration periods. We report and assess the weekly forecasts with quantified uncertainty from the top-ranked, second-ranked, and ensemble sub-epidemic models. Overall, we conducted 324 weekly sequential 4-week ahead forecasts across the models from the week of July 28th, 2022, to the week of October 13th, 2022. RESULTS: The last 10 of 12 forecasting periods (starting the week of August 11th, 2022) show either a plateauing or declining trend of monkeypox cases for all models and areas of study. According to our latest 4-week ahead forecast from the top-ranked model, a total of 6232 (95% PI 487.8, 12,468.0) cases could be added globally from the week of 10/20/2022 to the week of 11/10/2022. At the country level, the top-ranked model predicts that the USA will report the highest cumulative number of new cases for the 4-week forecasts (median based on OWID data: 1806 (95% PI 0.0, 5544.5)). The top-ranked and weighted ensemble models outperformed all other models in short-term forecasts. CONCLUSIONS: Our top-ranked model consistently predicted a decreasing trend in monkeypox cases on the global and country-specific scale during the last ten sequential forecasting periods. Our findings reflect the potential impact of increased immunity, and behavioral modification among high-risk populations.

Age-Stratified Model to Assess Health Outcomes of COVID-19 Vaccination Strategies, Ghana.

We assessed the effect of various COVID-19 vaccination strategies on health outcomes in Ghana by using an age-stratified compartmental model. We stratified the population into 3 age groups: <25 years, 25-64 years, and ≥65 years. We explored 5 vaccination optimization scenarios using 2 contact matrices, assuming that 1 million persons could be vaccinated in either 3 or 6 months. We assessed these vaccine optimization strategies for the initial strain, followed by a sensitivity analysis for the Delta variant. We found that vaccinating persons <25 years of age was associated with the lowest cumulative infections for the main matrix, for both the initial strain and the Delta variant. Prioritizing the elderly (≥65 years of age) was associated with the lowest cumulative deaths for both strains in all scenarios. The consensus between the findings of both contact matrices depended on the vaccine rollout period and the objective of the vaccination program.