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.

Death on the permafrost: Revisiting the 1918-20 influenza pandemic in Alaska using death certificates.

The 1918-20 influenza pandemic devastated Alaska's Indigenous populations. We report on quantitative analyses of pandemic deaths due to pneumonia and influenza (P&I) using information from Alaska death certificates dating between 1915 and 1921 (n=7,147). Goals include a reassessment of pandemic death numbers, analysis of P&I deaths beyond 1919, estimates of excess mortality patterns overall and by age using intercensal population estimates based on Alaska's demographic history, and comparisons between Alaska Native (AN) and non-AN residents. Results indicate that ANs experienced 83% of all P&I deaths and 87% of all-cause excess deaths during the pandemic. AN mortality was 8.1 times higher than non-AN mortality. Analyses also uncovered previously unknown mortality peaks in 1920. Both subpopulations showed characteristically high mortality of young adults, possibly due to imprinting with the 1889-90 pandemic virus, but their age-specific mortality patterns were different: non-AN mortality declined after age 25-29 and stayed relatively low for the elderly, while AN mortality increased after age 25-29, peaked at age 40-44, and remained high up to age 64. This suggests a relative lack of exposure to H1-type viruses pre-1889 among AN persons. In contrast, non-AN persons, often temporary residents, may have gained immunity before moving to Alaska.

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.

SpatialWavePredict: a tutorial-based primer and toolbox for forecasting growth trajectories using the ensemble spatial wave sub-epidemic modeling framework.

BACKGROUND: Dynamical mathematical models defined by a system of differential equations are typically not easily accessible to non-experts. However, forecasts based on these types of models can help gain insights into the mechanisms driving the process and may outcompete simpler phenomenological growth models. Here we introduce a friendly toolbox, SpatialWavePredict, to characterize and forecast the spatial wave sub-epidemic model, which captures diverse wave dynamics by aggregating multiple asynchronous growth processes and has outperformed simpler phenomenological growth models in short-term forecasts of various infectious diseases outbreaks including SARS, Ebola, and the early waves of the COVID-19 pandemic in the US. RESULTS: This tutorial-based primer introduces and illustrates a user-friendly MATLAB toolbox for fitting and forecasting time-series trajectories using an ensemble spatial wave sub-epidemic model based on ordinary differential equations. Scientists, policymakers, and students can use the toolbox to conduct real-time short-term forecasts. The five-parameter epidemic wave model in the toolbox aggregates linked overlapping sub-epidemics and captures a rich spectrum of epidemic wave dynamics, including oscillatory wave behavior and plateaus. An ensemble strategy aims to improve forecasting performance by combining the resulting top-ranked models. The toolbox provides a tutorial for forecasting time-series trajectories, 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, 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 an ensemble spatial wave sub-epidemic wave model. As an epidemic situation or contagion occurs, the tools presented in this tutorial can facilitate policymakers to guide the implementation of containment strategies and assess the impact of control interventions. We demonstrate the functionality of the toolbox with examples, including a tutorial video, and is illustrated using daily data on the COVID-19 pandemic in the USA.

Assessment of fever screening at airports in detecting domestic passengers infected with SARS-CoV-2, 2020-2022, Okinawa prefecture, Japan.

BACKGROUND: While airport screening measures for COVID-19 infected passengers at international airports worldwide have been greatly relaxed, observational studies evaluating fever screening alone at airports remain scarce. The purpose of this study is to retrospectively assess the effectiveness of fever screening at airports in preventing the influx of COVID-19 infected persons. METHODS: We conducted a retrospective epidemiological analysis of fever screening implemented at 9 airports in Okinawa Prefecture from May 2020 to March 2022. The number of passengers covered during the same period was 9,003,616 arriving at 9 airports in Okinawa Prefecture and 5,712,983 departing passengers at Naha Airport. The capture rate was defined as the proportion of reported COVID-19 cases who would have passed through airport screening to the number of suspected cases through fever screening at the airport, and this calculation used passengers arriving at Naha Airport and surveillance data collected by Okinawa Prefecture between May 2020 and March 2021. RESULTS: From May 2020 to March 2021, 4.09 million people were reported to pass through airports in Okinawa. During the same period, at least 122 people with COVID-19 infection arrived at the airports in Okinawa, but only a 10 suspected cases were detected; therefore, the capture rate is estimated to be up to 8.2% (95% CI: 4.00-14.56%). Our result of a fever screening rate is 0.0002% (95%CI: 0.0003-0.0006%) (10 suspected cases /2,971,198 arriving passengers). The refusal rate of passengers detected by thermography who did not respond to temperature measurements was 0.70% (95% CI: 0.19-1.78%) (4 passengers/572 passengers). CONCLUSIONS: This study revealed that airport screening based on thermography alone missed over 90% of COVID-19 infected cases, indicating that thermography screening may be ineffective as a border control measure. The fact that only 10 febrile cases were detected after screening approximately 3 million passengers suggests the need to introduce measures targeting asymptomatic infections, especially with long incubation periods. Therefore, other countermeasures, e.g. preboarding RT-PCR testing, are highly recommended during an epidemic satisfying World Health Organization (WHO) Public Health Emergency of International Concern (PHEIC) criteria with pathogen characteristics similar or exceeding SARS-CoV-2, especially when traveling to rural cities with limited medical resources.

Dynamic prioritization of COVID-19 vaccines when social distancing is limited for essential workers.

COVID-19 vaccines have been authorized in multiple countries, and more are under rapid development. Careful design of a vaccine prioritization strategy across sociodemographic groups is a crucial public policy challenge given that 1) vaccine supply will be constrained for the first several months of the vaccination campaign, 2) there are stark differences in transmission and severity of impacts from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) across groups, and 3) SARS-CoV-2 differs markedly from previous pandemic viruses. We assess the optimal allocation of a limited vaccine supply in the United States across groups differentiated by age and essential worker status, which constrains opportunities for social distancing. We model transmission dynamics using a compartmental model parameterized to capture current understanding of the epidemiological characteristics of COVID-19, including key sources of group heterogeneity (susceptibility, severity, and contact rates). We investigate three alternative policy objectives (minimizing infections, years of life lost, or deaths) and model a dynamic strategy that evolves with the population epidemiological status. We find that this temporal flexibility contributes substantially to public health goals. Older essential workers are typically targeted first. However, depending on the objective, younger essential workers are prioritized to control spread or seniors to directly control mortality. When the objective is minimizing deaths, relative to an untargeted approach, prioritization averts deaths on a range between 20,000 (when nonpharmaceutical interventions are strong) and 300,000 (when these interventions are weak). We illustrate how optimal prioritization is sensitive to several factors, most notably, vaccine effectiveness and supply, rate of transmission, and the magnitude of initial infections.

The Nicaraguan Pediatric Influenza Cohort Study, 2011-2019: Influenza Incidence, Seasonality, and Transmission.

BACKGROUND: Children account for a large portion of global influenza burden and transmission, and a better understanding of influenza in children is needed to improve prevention and control strategies. METHODS: To examine the incidence and transmission of influenza we conducted a prospective community-based study of children aged 0-14 years in Managua, Nicaragua, between 2011 and 2019. Participants were provided with medical care through study physicians and symptomatic influenza was confirmed by reverse-transcription polymerase chain reaction (RT-PCR). Wavelet analyses were used to examine seasonality. Generalized growth models (GGMs) were used to estimate effective reproduction numbers. RESULTS: From 2011 to 2019, 3016 children participated, with an average of ∼1800 participants per year and median follow-up time of 5 years per child, and 48.3% of the cohort in 2019 had been enrolled their entire lives. The overall incidence rates per 100 person-years were 14.5 symptomatic influenza cases (95% confidence interval [CI]: 13.9-15.1) and 1.0 influenza-associated acute lower respiratory infection (ALRI) case (95% CI: .8-1.1). Symptomatic influenza incidence peaked at age 9-11 months. Infants born during peak influenza circulation had lower incidence in the first year of their lives. The mean effective reproduction number was 1.2 (range 1.02-1.49), and we observed significant annual patterns for influenza and influenza A, and a 2.5-year period for influenza B. CONCLUSIONS: This study provides important information for understanding influenza epidemiology and informing influenza vaccine policy. These results will aid in informing strategies to reduce the burden of influenza.

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.

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.

An ensemble n-sub-epidemic modeling framework for short-term forecasting epidemic trajectories: Application to the COVID-19 pandemic in the USA.

We analyze an ensemble of n-sub-epidemic modeling for forecasting the trajectory of epidemics and pandemics. These ensemble modeling approaches, and models that integrate sub-epidemics to capture complex temporal dynamics, have demonstrated powerful forecasting capability. This modeling framework can characterize complex epidemic patterns, including plateaus, epidemic resurgences, and epidemic waves characterized by multiple peaks of different sizes. We systematically assess their calibration and short-term forecasting performance in short-term forecasts for the COVID-19 pandemic in the USA from late April 2020 to late February 2022. We compare their performance with two commonly used statistical ARIMA models. The best fit sub-epidemic model and three ensemble models constructed using the top-ranking sub-epidemic models consistently outperformed the ARIMA models in terms of the weighted interval score (WIS) and the coverage of the 95% prediction interval across the 10-, 20-, and 30-day short-term forecasts. In our 30-day forecasts, the average WIS ranged from 377.6 to 421.3 for the sub-epidemic models, whereas it ranged from 439.29 to 767.05 for the ARIMA models. Across 98 short-term forecasts, the ensemble model incorporating the top four ranking sub-epidemic models (Ensemble(4)) outperformed the (log) ARIMA model 66.3% of the time, and the ARIMA model, 69.4% of the time in 30-day ahead forecasts in terms of the WIS. Ensemble(4) consistently yielded the best performance in terms of the metrics that account for the uncertainty of the predictions. This framework can be readily applied to investigate the spread of epidemics and pandemics beyond COVID-19, as well as other dynamic growth processes found in nature and society that would benefit from short-term predictions.

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.

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.

Comparative epidemiology of five waves of COVID-19 in Mexico, March 2020-August 2022.

BACKGROUND: The Mexican Institute of Social Security (IMSS) is the largest health care provider in Mexico, covering about 48% of the Mexican population. In this report, we describe the epidemiological patterns related to confirmed cases, hospitalizations, intubations, and in-hospital mortality due to COVID-19 and associated factors, during five epidemic waves recorded in the IMSS surveillance system. METHODS: We analyzed COVID-19 laboratory-confirmed cases from the Online Epidemiological Surveillance System (SINOLAVE) from March 29th, 2020, to August 27th, 2022. We constructed weekly epidemic curves describing temporal patterns of confirmed cases and hospitalizations by age, gender, and wave. We also estimated hospitalization, intubation, and hospital case fatality rates. The mean days of in-hospital stay and hospital admission delay were calculated across five pandemic waves. Logistic regression models were employed to assess the association between demographic factors, comorbidities, wave, and vaccination and the risk of severe disease and in-hospital death. RESULTS: A total of 3,396,375 laboratory-confirmed COVID-19 cases were recorded across the five waves. The introduction of rapid antigen testing at the end of 2020 increased detection and modified epidemiological estimates. Overall, 11% (95% CI 10.9, 11.1) of confirmed cases were hospitalized, 20.6% (95% CI 20.5, 20.7) of the hospitalized cases were intubated, and the hospital case fatality rate was 45.1% (95% CI 44.9, 45.3). The mean in-hospital stay was 9.11 days, and patients were admitted on average 5.07 days after symptoms onset. The most recent waves dominated by the Omicron variant had the highest incidence. Hospitalization, intubation, and mean hospitalization days decreased during subsequent waves. The in-hospital case fatality rate fluctuated across waves, reaching its highest value during the second wave in winter 2020. A notable decrease in hospitalization was observed primarily among individuals ≥ 60 years. The risk of severe disease and death was positively associated with comorbidities, age, and male gender; and declined with later waves and vaccination status. CONCLUSION: During the five pandemic waves, we observed an increase in the number of cases and a reduction in severity metrics. During the first three waves, the high in-hospital fatality rate was associated with hospitalization practices for critical patients with comorbidities.

Ensemble bootstrap methodology for forecasting dynamic growth processes using differential equations: application to epidemic outbreaks.

BACKGROUND: Ensemble modeling aims to boost the forecasting performance by systematically integrating the predictive accuracy across individual models. Here we introduce a simple-yet-powerful ensemble methodology for forecasting the trajectory of dynamic growth processes that are defined by a system of non-linear differential equations with applications to infectious disease spread. METHODS: We propose and assess the performance of two ensemble modeling schemes with different parametric bootstrapping procedures for trajectory forecasting and uncertainty quantification. Specifically, we conduct sequential probabilistic forecasts to evaluate their forecasting performance using simple dynamical growth models with good track records including the Richards model, the generalized-logistic growth model, and the Gompertz model. We first test and verify the functionality of the method using simulated data from phenomenological models and a mechanistic transmission model. Next, the performance of the method is demonstrated using a diversity of epidemic datasets including scenario outbreak data of the Ebola Forecasting Challenge and real-world epidemic data outbreaks of including influenza, plague, Zika, and COVID-19. RESULTS: We found that the ensemble method that randomly selects a model from the set of individual models for each time point of the trajectory of the epidemic frequently outcompeted the individual models as well as an alternative ensemble method based on the weighted combination of the individual models and yields broader and more realistic uncertainty bounds for the trajectory envelope, achieving not only better coverage rate of the 95% prediction interval but also improved mean interval scores across a diversity of epidemic datasets. CONCLUSION: Our new methodology for ensemble forecasting outcompete component models and an alternative ensemble model that differ in how the variance is evaluated for the generation of the prediction intervals of the forecasts.

Characterizing all-cause excess mortality patterns during COVID-19 pandemic in Mexico.

BACKGROUND: Low testing rates and delays in reporting hinder the estimation of the mortality burden associated with the COVID-19 pandemic. During a public health emergency, estimating all cause excess deaths above an expected level of death can provide a more reliable picture of the mortality burden. Here, we aim to estimate the absolute and relative mortality impact of COVID-19 pandemic in Mexico. METHODS: We obtained weekly mortality time series due to all causes for Mexico, and by gender, and geographic region from 2015 to 2020. We also compiled surveillance data on COVID-19 cases and deaths to assess the timing and intensity of the pandemic and assembled weekly series of the proportion of tweets about 'death' from Mexico to assess the correlation between people's media interaction about 'death' and the rise in pandemic deaths. We estimated all-cause excess mortality rates and mortality rate ratio increase over baseline by fitting Serfling regression models and forecasted the total excess deaths for Mexico for the first 4 weeks of 2021 using the generalized logistic growth model. RESULTS: We estimated the all-cause excess mortality rate associated with the COVID-19 pandemic in Mexico in 2020 at 26.10 per 10,000 population, which corresponds to 333,538 excess deaths. Males had about 2-fold higher excess mortality rate (33.99) compared to females (18.53). Mexico City reported the highest excess death rate (63.54) and RR (2.09) compared to rest of the country (excess rate = 23.25, RR = 1.62). While COVID-19 deaths accounted for only 38.64% of total excess deaths in Mexico, our forecast estimate that Mexico has accumulated a total of ~ 61,610 [95% PI: 60,003, 63,216] excess deaths in the first 4 weeks of 2021. Proportion of tweets was significantly correlated with the excess mortality (ρ = 0.508 [95% CI: 0.245, 0.701], p-value = 0.0004). CONCLUSION: The COVID-19 pandemic has heavily affected Mexico. The lab-confirmed COVID-19 deaths accounted for only 38.64% of total all cause excess deaths (333,538) in Mexico in 2020. This reflects either the effect of low testing rates in Mexico, or the surge in number of deaths due to other causes during the pandemic. A model-based forecast indicates that an average of 61,610 excess deaths have occurred in January 2021.

Urbanization prolongs hantavirus epidemics in cities.

Urbanization and rural-urban migration are two factors driving global patterns of disease and mortality. There is significant concern about their potential impact on disease burden and the effectiveness of current control approaches. Few attempts have been made to increase our understanding of the relationship between urbanization and disease dynamics, although it is generally believed that urban living has contributed to reductions in communicable disease burden in industrialized countries. To investigate this relationship, we carried out spatiotemporal analyses using a 48-year-long dataset of hemorrhagic fever with renal syndrome incidence (HFRS; mainly caused by two serotypes of hantavirus in China: Hantaan virus and Seoul virus) and population movements in an important endemic area of south China during the period 1963-2010. Our findings indicate that epidemics coincide with urbanization, geographic expansion, and migrant movement over time. We found a biphasic inverted U-shaped relationship between HFRS incidence and urbanization, with various endemic turning points associated with economic growth rates in cities. Our results revealed the interrelatedness of urbanization, migration, and hantavirus epidemiology, potentially explaining why urbanizing cities with high economic growth exhibit extended epidemics. Our results also highlight contrasting effects of urbanization on zoonotic disease outbreaks during periods of economic development in China.

Estimating Risk for Death from Coronavirus Disease, China, January-February 2020.

Since December 2019, when the first case of coronavirus disease (COVID-19) was identified in the city of Wuhan in the Hubei Province of China, the epidemic has generated tens of thousands of cases throughout China. As of February 28, 2020, the cumulative number of reported deaths in China was 2,858. We estimated the time-delay adjusted risk for death from COVID-19 in Wuhan, as well as for China excluding Wuhan, to assess the severity of the epidemic in the country. Our estimates of the risk for death in Wuhan reached values as high as 12% in the epicenter of the epidemic and ≈1% in other, more mildly affected areas. The elevated death risk estimates are probably associated with a breakdown of the healthcare system, indicating that enhanced public health interventions, including social distancing and movement restrictions, should be implemented to bring the COVID-19 epidemic under control.

Temporary Fertility Decline after Large Rubella Outbreak, Japan.

Japan experienced 2 large rubella epidemics in 2004 and 2012-2014. Because of suboptimal immunization levels, the country has been experiencing a third major outbreak during 2018-2020. We conducted time series analyses to evaluate the effect of the 2012-2014 nationwide rubella epidemic on prefecture-level natality in Japan. We identified a statistically significant decline in fertility rates associated with rubella epidemic activity and increased Google searches for the term "rubella." We noted that the timing of fertility declines in 2014 occurred 9-13 months after peak rubella incidence months in 2013 in 4 prefectures with the highest rubella incidence. Public health interventions should focus on enhancing vaccination campaigns against rubella, not only to protect pregnant women from infection but also to mitigate declines in population size and birth rates.

Severe Acute Respiratory Syndrome Coronavirus 2 Transmission Potential, Iran, 2020.

To determine the transmission potential of severe acute respiratory syndrome coronavirus 2 in Iran in 2020, we estimated the reproduction number as 4.4 (95% CI 3.9-4.9) by using a generalized growth model and 3.5 (95% CI 1.3-8.1) by using epidemic doubling time. The reproduction number decreased to 1.55 after social distancing interventions were implemented.

Doubling Time of the COVID-19 Epidemic by Province, China.

In China, the doubling time of the coronavirus disease epidemic by province increased during January 20-February 9, 2020. Doubling time estimates ranged from 1.4 (95% CI 1.2-2.0) days for Hunan Province to 3.1 (95% CI 2.1-4.8) days for Xinjiang Province. The estimate for Hubei Province was 2.5 (95% CI 2.4-2.6) days.