Nowcasting and forecasting the 2022 U.S. mpox outbreak: Support for public health decision making and lessons learned.

In June of 2022, the U.S. Centers for Disease Control and Prevention (CDC) Mpox Response wanted timely answers to important epidemiological questions which can now be answered more effectively through infectious disease modeling. Infectious disease models have shown to be valuable tools for decision making during outbreaks; however, model complexity often makes communicating the results and limitations of models to decision makers difficult. We performed nowcasting and forecasting for the 2022 mpox outbreak in the United States using the R package EpiNow2. We generated nowcasts/forecasts at the national level, by Census region, and for jurisdictions reporting the greatest number of mpox cases. Modeling results were shared for situational awareness within the CDC Mpox Response and publicly on the CDC website. We retrospectively evaluated forecast predictions at four key phases (early, exponential growth, peak, and decline) during the outbreak using three metrics, the weighted interval score, mean absolute error, and prediction interval coverage. We compared the performance of EpiNow2 with a naïve Bayesian generalized linear model (GLM). The EpiNow2 model had less probabilistic error than the GLM during every outbreak phase except for the early phase. We share our experiences with an existing tool for nowcasting/forecasting and highlight areas of improvement for the development of future tools. We also reflect on lessons learned regarding data quality issues and adapting modeling results for different audiences.

Updating Reproduction Number Estimates for Mpox in the Democratic Republic of Congo Using Surveillance Data.

Incidence of human monkeypox (mpox) has been increasing in West and Central Africa, including in the Democratic Republic of Congo (DRC), where monkeypox virus (MPXV) is endemic. Most estimates of the pathogen's transmissibility in the DRC are based on data from the 1980s. Amid the global 2022 mpox outbreak, new estimates are needed to characterize the virus' epidemic potential and inform outbreak control strategies. We used the R package vimes to identify clusters of laboratory-confirmed mpox cases in Tshuapa Province, DRC. Cases with both temporal and spatial data were assigned to clusters based on the disease's serial interval and spatial kernel. We used the size of the clusters to infer the effective reproduction number, Rt, and the rate of zoonotic spillover of MPXV into the human population. Out of 1,463 confirmed mpox cases reported in Tshuapa Province between 2013 and 2017, 878 had both date of symptom onset and a location with geographic coordinates. Results include an estimated Rt of 0.82 (95% CI: 0.79-0.85) and a rate of 132 (95% CI: 122-143) spillovers per year assuming a reporting rate of 25%. This estimate of Rt is larger than most previous estimates. One potential explanation for this result is that Rt could have increased in the DRC over time owing to declining population-level immunity conferred by smallpox vaccination, which was discontinued around 1982. Rt could be overestimated if our assumption of one spillover event per cluster does not hold. Our results are consistent with increased transmissibility of MPXV in Tshuapa Province.

Marburg virus disease outbreaks, mathematical models, and disease parameters: a systematic review.

The 2023 Marburg virus disease outbreaks in Equatorial Guinea and Tanzania highlighted the importance of better understanding this lethal pathogen. We did a systematic review (PROSPERO CRD42023393345) of peer-reviewed articles reporting historical outbreaks, modelling studies, and epidemiological parameters focused on Marburg virus disease. We searched PubMed and Web of Science from database inception to March 31, 2023. Two reviewers evaluated all titles and abstracts with consensus-based decision making. To ensure agreement, 13 (31%) of 42 studies were double-extracted and a custom-designed quality assessment questionnaire was used for risk of bias assessment. We present detailed information on 478 reported cases and 385 deaths from Marburg virus disease. Analysis of historical outbreaks and seroprevalence estimates suggests the possibility of undetected Marburg virus disease outbreaks, asymptomatic transmission, or cross-reactivity with other pathogens, or a combination of these. Only one study presented a mathematical model of Marburg virus transmission. We estimate an unadjusted, pooled total random effect case fatality ratio of 61·9% (95% CI 38·8-80·6; I2=93%). We identify epidemiological parameters relating to transmission and natural history, for which there are few estimates. This systematic review and the accompanying database provide a comprehensive overview of Marburg virus disease epidemiology and identify key knowledge gaps, contributing crucial information for mathematical models to support future Marburg virus disease epidemic responses.

Possible Exposures Among Mpox Patients Without Reported Male-to-Male Sexual Contact - Six U.S. Jurisdictions, November 1-December 14, 2022.

The extent to which the 2022 mpox outbreak has affected persons without a recent history of male-to-male sexual contact (MMSC) is not well understood. During November 1-December 14, 2022, CDC partnered with six jurisdictional health departments to characterize possible exposures among mpox patients aged ≥18 years who did not report MMSC during the 3 weeks preceding symptom onset. Among 52 patients included in the analysis, 14 (27%) had a known exposure to a person with mpox, including sexual activity and other close intimate contact (eight) and household contact (six). Among 38 (73%) patients with no known exposure to a person with mpox, self-reported activities before illness onset included sexual activity and other close intimate contact (17; 45%), close face-to-face contact (14; 37%), attending large social gatherings (11; 29%), and being in occupational settings involving close skin-to-skin contact (10; 26%). These findings suggest that sexual activity remains an important route of mpox exposure among patients who do not report MMSC.

Risk of Rabies and Implications for Postexposure Prophylaxis Administration in the US.

Importance: In the US, rabies postexposure prophylaxis (PEP) is often administered without a comprehensive and regionally appropriate rabies risk assessment. For low-risk exposures, this can result in patients incurring out-of-pocket expenses or experiencing adverse effects of PEP unnecessarily. Objective: To use a model to estimate (1) the probability that an animal would test positive for rabies virus (RABV) given that a person was exposed, and (2) the probability that a person would die from rabies given that they were exposed to a suspect rabid animal and did not receive PEP, and to propose a risk threshold for recommending PEP according to model estimates and a survey. Design, Setting, and Participants: In this decision analytical modeling study, positivity rates were calculated using more than 900 000 animal samples tested for RABV between 2011 and 2020. Other parameters were estimated from a subset of the surveillance data and the literature. Probabilities were estimated using Bayes' rule. A survey was administered among a convenience sample of state public health officials in all US states (excluding Hawaii) plus Washington, DC and Puerto Rico to determine a risk threshold for PEP recommendation. Respondents were asked whether they would recommend PEP given 24 standardized exposure scenarios while accounting for local rabies epidemiology. Main Outcomes and Measures: A quantitative and regionally appropriate approach for helping health care practitioners and public health professionals determine whether to recommend and/or administer rabies PEP. Results: A total of 1728 unique observations were obtained from the model for the probability that an animal would test positive for RABV given that a person was exposed, and 41 472 for ) the probability that a person would die from rabies given that they were exposed to a suspect rabid animal and did not receive PEP. The median probability that an animal would test positive for RABV given that a person was exposed ranged from 3 × 10-7 to 0.97, while the probability that a person would die from rabies given that they were exposed to a suspect rabid animal and did not receive PEP ranged from 1 × 10-10 to 0.55. Fifty public health officials out of a target sample size of 102 responded to the survey. Using logistic regression, a risk threshold was estimated for PEP recommendation of 0.0004; PEP may not be recommended for exposures with probabilities below this threshold. Conclusions and Relevance: In this modeling study of rabies in the US, the risk of death|exposure was quantified and a risk threshold was estimated. These results could be used to inform the decision-making process as to the appropriateness of recommending rabies PEP.

Serial Interval and Incubation Period Estimates of Monkeypox Virus Infection in 12 Jurisdictions, United States, May-August 2022.

Using data from 12 US health departments, we estimated mean serial interval for monkeypox virus infection to be 8.5 (95% credible interval 7.3-9.9) days for symptom onset, based on 57 case pairs. Mean estimated incubation period was 5.6 (95% credible interval 4.3-7.8) days for symptom onset, based on 35 case pairs.

Correction: The epidemiology of Mayaro virus in the Americas: A systematic review and key parameter estimates for outbreak modelling.

[This corrects the article DOI: 10.1371/journal.pntd.0009418.].

Estimating Zika virus attack rates and risk of Zika virus-associated neurological complications in Colombian capital cities with a Bayesian model.

Zika virus (ZIKV) is a mosquito-borne pathogen that caused a major epidemic in the Americas in 2015-2017. Although the majority of ZIKV infections are asymptomatic, the virus has been associated with congenital birth defects and neurological complications (NC) in adults. We combined multiple data sources to improve estimates of ZIKV infection attack rates (IARs), reporting rates of Zika virus disease (ZVD) and the risk of ZIKV-associated NC for 28 capital cities in Colombia. ZVD surveillance data were combined with post-epidemic seroprevalence data and a dataset on ZIKV-associated NC in a Bayesian hierarchical model. We found substantial heterogeneity in ZIKV IARs across cities. The overall estimated ZIKV IAR across the 28 cities was 0.38 (95% CrI: 0.17-0.92). The estimated ZVD reporting rate was 0.013 (95% CrI: 0.004-0.024), and 0.51 (95% CrI: 0.17-0.92) cases of ZIKV-associated NC were estimated to be reported per 10 000 ZIKV infections. When we assumed the same ZIKV IAR across sex or age group, we found important spatial heterogeneities in ZVD reporting rates and the risk of being reported as a ZVD case with NC. Our results highlight how additional data sources can be used to overcome biases in surveillance data and estimate key epidemiological parameters.

Spatial and temporal invasion dynamics of the 2014-2017 Zika and chikungunya epidemics in Colombia.

Zika virus (ZIKV) and chikungunya virus (CHIKV) were recently introduced into the Americas resulting in significant disease burdens. Understanding their spatial and temporal dynamics at the subnational level is key to informing surveillance and preparedness for future epidemics. We analyzed anonymized line list data on approximately 105,000 Zika virus disease and 412,000 chikungunya fever suspected and laboratory-confirmed cases during the 2014-2017 epidemics. We first determined the week of invasion in each city. Out of 1,122, 288 cities met criteria for epidemic invasion by ZIKV and 338 cities by CHIKV. We analyzed risk factors for invasion using linear and logistic regression models. We also estimated that the geographic origin of both epidemics was located in Barranquilla, north Colombia. We assessed the spatial and temporal invasion dynamics of both viruses to analyze transmission between cities using a suite of (i) gravity models, (ii) Stouffer's rank models, and (iii) radiation models with two types of distance metrics, geographic distance and travel time between cities. Invasion risk was best captured by a gravity model when accounting for geographic distance and intermediate levels of density dependence; Stouffer's rank model with geographic distance performed similarly well. Although a few long-distance invasion events occurred at the beginning of the epidemics, an estimated distance power of 1.7 (95% CrI: 1.5-2.0) from the gravity models suggests that spatial spread was primarily driven by short-distance transmission. Similarities between the epidemics were highlighted by jointly fitted models, which were preferred over individual models when the transmission intensity was allowed to vary across arboviruses. However, ZIKV spread considerably faster than CHIKV.

Descriptive analysis of surveillance data for Zika virus disease and Zika virus-associated neurological complications in Colombia, 2015-2017.

Zika virus (ZIKV) is a mosquito-borne pathogen that recently caused a major epidemic in the Americas. Although the majority of ZIKV infections are asymptomatic, the virus has been associated with birth defects in fetuses and newborns of infected mothers as well as neurological complications in adults. We performed a descriptive analysis on approximately 106,000 suspected and laboratory-confirmed cases of Zika virus disease (ZVD) that were reported during the 2015-2017 epidemic in Colombia. We also analyzed a dataset containing patients with neurological complications and recent febrile illness compatible with ZVD. Females had higher cumulative incidence of ZVD than males. Compared to the general population, cases were more likely to be reported in young adults (20 to 39 years of age). We estimated the cumulative incidence of ZVD in pregnant females at 3,120 reported cases per 100,000 population (95% CI: 3,077-3,164), which was considerably higher than the incidence in both males and non-pregnant females. ZVD cases were reported in all 32 departments. Four-hundred and eighteen patients suffered from ZIKV-associated neurological complications, of which 85% were diagnosed with Guillain-Barré syndrome. The median age of ZIKV cases with neurological complications was 12 years older than that of ZVD cases. ZIKV-associated neurological complications increased with age, and the highest incidence was reported among individuals aged 75 and older. Even though neurological complications and deaths due to ZIKV were rare in this epidemic, better risk communication is needed for people living in or traveling to ZIKV-affected areas.

The epidemiology of Mayaro virus in the Americas: A systematic review and key parameter estimates for outbreak modelling.

Mayaro virus (MAYV) is an arbovirus that is endemic to tropical forests in Central and South America, particularly within the Amazon basin. In recent years, concern has increased regarding MAYV's ability to invade urban areas and cause epidemics across the region. We conducted a systematic literature review to characterise the evolutionary history of MAYV, its transmission potential, and exposure patterns to the virus. We analysed data from the literature on MAYV infection to produce estimates of key epidemiological parameters, including the generation time and the basic reproduction number, R0. We also estimated the force-of-infection (FOI) in epidemic and endemic settings. Seventy-six publications met our inclusion criteria. Evidence of MAYV infection in humans, animals, or vectors was reported in 14 Latin American countries. Nine countries reported evidence of acute infection in humans confirmed by viral isolation or reverse transcription-PCR (RT-PCR). We identified at least five MAYV outbreaks. Seroprevalence from population based cross-sectional studies ranged from 21% to 72%. The estimated mean generation time of MAYV was 15.2 days (95% CrI: 11.7-19.8) with a standard deviation of 6.3 days (95% CrI: 4.2-9.5). The per-capita risk of MAYV infection (FOI) ranged between 0.01 and 0.05 per year. The mean R0 estimates ranged between 2.1 and 2.9 in the Amazon basin areas and between 1.1 and 1.3 in the regions outside of the Amazon basin. Although MAYV has been identified in urban vectors, there is not yet evidence of sustained urban transmission. MAYV's enzootic cycle could become established in forested areas within cities similar to yellow fever virus.