Excess Mortality in the United States, 2020-21: County-level Estimates for Population Groups and Associations with Social Vulnerability.

To assess the excess mortality burden of Covid-19 in the United States, we estimated sex, age and race stratified all-cause excess deaths in each county of the US during 2020 and 2021. Using spatial Bayesian models trained on all recorded deaths between 2003-2019, we estimated 463,187 (95% uncertainty interval (UI): 426,139 - 497,526) excess deaths during 2020, and 544,105 (95% UI: 492,202 - 592,959) excess deaths during 2021 nationally, with considerable geographical heterogeneity. Excess mortality rate (EMR) nearly doubled for each 10-year increase in age and was consistently higher among men than women. EMR in the Black population was 1.5 times that of the White population nationally and as high as 3.8 times in some states. Among the 25-54 year population excess mortality was highest in the American Indian/Alaskan Native (AI/AN) population among the four racial groups studied, and in a few states was as high as 6 times that of the White population. Strong association of EMR with county-level social vulnerability was estimated, including positive associations with prevalence of disability (standardized effect: 40.6 excess deaths per 100,000), older population (37.6), poverty (23.6), and unemployment (18.5), whereas population density (-50), higher education (-38.6), and income (-35.4) were protective. Together, these estimates provide a more reliable and comprehensive understanding of the mortality burden of the pandemic in the US thus far. They suggest that Covid-19 amplified social and racial disparities. Short-term measures to protect more vulnerable groups in future Covid-19 waves and systemic corrective steps to address long-term societal inequities are necessary.

Age-dependent seasonality in the incidence of stroke: A 21-year population-based study.

INTRODUCTION: Seasonality in the incidence of stroke has been examined in numerous studies, yet data on whether it differs with age are limited. To fill this gap, we utilized a largescale dataset from Israel. PATIENTS AND METHODS: We retrieved data of all hospitalizations for ischemic stroke (IS), transient ischemic attack (TIA) and intra cerebral hemorrhage (ICH) from 2000 to 2020. We maintained separate datasets for IS/TIA and ICH, divided into five age groups: 18-49, 50-59, 60-69, 70-79, and 80+. We modeled the monthly incidence using a generalized additive model. The seasonal effect was defined by the rate ratio (RR) of each month compared to the annual mean. RESULTS: The analysis included 317,586 and 23,789 events of IS/TIA and ICH respectively. We found an interaction between age and seasonality, accounting for a phase shift with age in the seasonal pattern of IS/TIA incidence. For cases under 70 years, the peak was during summertime and the RRs increased with decreasing age, reaching 1.11 (95% CI 1.09-1.13) at the youngest age group. In contrast, among the elderly, a winter peak was observed and the RRs increased with age to 1.07 (95% CI 1.06-1.08) at the oldest age group. For ICH, a winter/autumn peak was identified and the RRs increased with age to 1.20 (95% CI 1.16-1.24). CONCLUSIONS: Our finding of age-dependent seasonal patterns in the occurrence of stroke, suggests closer monitoring of cardiovascular risk factors during wintertime among elderly individuals. The mechanism governing the seasonal phase shift with age in IS/TIA incidence, requires further investigation.

Evaluation of FluSight influenza forecasting in the 2021-22 and 2022-23 seasons with a new target laboratory-confirmed influenza hospitalizations.

Accurate forecasts can enable more effective public health responses during seasonal influenza epidemics. Forecasting teams were asked to provide national and jurisdiction-specific probabilistic predictions of weekly confirmed influenza hospital admissions for one through four weeks ahead for the 2021-22 and 2022-23 influenza seasons. Across both seasons, 26 teams submitted forecasts, with the submitting teams varying between seasons. Forecast skill was evaluated using the Weighted Interval Score (WIS), relative WIS, and coverage. Six out of 23 models outperformed the baseline model across forecast weeks and locations in 2021-22 and 12 out of 18 models in 2022-23. Averaging across all forecast targets, the FluSight ensemble was the 2nd most accurate model measured by WIS in 2021-22 and the 5th most accurate in the 2022-23 season. Forecast skill and 95% coverage for the FluSight ensemble and most component models degraded over longer forecast horizons and during periods of rapid change. Current influenza forecasting efforts help inform situational awareness, but research is needed to address limitations, including decreased performance during periods of changing epidemic dynamics.

Data-driven clustering of infectious disease incidence into age groups.

Understanding the patterns of infectious diseases spread in the population is an important element of mitigation and vaccination programs. A major and common characteristic of most infectious diseases is age-related heterogeneity in the transmission, which potentially can affect the dynamics of an epidemic as manifested by the pattern of disease incidence in different age groups. Currently there are no statistical criteria of how to partition the disease incidence data into clusters. We develop the first data-driven methodology for deciding on the best partition of incidence data into age-groups, in a well defined statistical sense. The method employs a top-down hierarchical partitioning algorithm, with a stopping criteria based on multiple hypotheses significance testing controlling the family wise error rate. The type one error and statistical power of the method are tested using simulations. The method is then applied to Covid-19 incidence data in Israel, in order to extract the significant age-group clusters in each wave of the epidemic.

Population-level implications of the Israeli booster campaign to curtail COVID-19 resurgence.

Israel was one of the first countries to administer mass vaccination against severe acute respiratory syndrome coronavirus 2. Consequently, it was among the first countries to experience substantial breakthrough infections due to the waning of vaccine-induced immunity, which led to a resurgence of the epidemic. In response, Israel launched a booster campaign to mitigate the outbreak and was the first country to do so. Israel's success in curtailing the Delta resurgence while imposing only mild nonpharmaceutical interventions influenced the decision of many countries to initiate a booster campaign. By constructing a detailed mathematical model and calibrating it to the Israeli data, we extend the understanding of the impact of the booster campaign from the individual to the population level. We used the calibrated model to explore counterfactual scenarios in which the booster vaccination campaign is altered by changing the eligibility criteria or the start time of the campaign and to assess the direct and indirect effects in the different scenarios. The results point to the vast benefits of vaccinating younger age groups that are not at a high risk of developing severe disease but play an important role in transmission. We further show that, when the epidemic is exponentially growing, the success of the booster campaign is highly sensitive to the timing of its initiation. Hence, a rapid response is an important factor in reducing disease burden using booster vaccination.

Similarity of Protection Conferred by Previous SARS-CoV-2 Infection and by BNT162b2 Vaccine: A 3-Month Nationwide Experience From Israel.

The worldwide shortage of vaccination against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection while the pandemic still remains uncontrolled has led many countries to the dilemma of whether or not to vaccinate previously infected persons. Understanding the level of protection conferred by previous infection compared with that of vaccination is important for policy-making. We analyzed an updated individual-level database of the entire population of Israel to assess the protection provided by both prior infection and vaccination in preventing subsequent SARS-CoV-2 infection, hospitalization with coronavirus disease 2019 (COVID-19), severe disease, and death due to COVID-19. Outcome data were collected from December 20, 2020, to March 20, 2021. Vaccination was highly protective, with overall estimated effectiveness of 94.5% (95% confidence interval (CI): 94.3, 94.7) for documented infection, 95.8% (95% CI: 95.2, 96.2) for hospitalization, 96.3% (95% CI: 95.7, 96.9) for severe illness, and 96.0% (95% CI: 94.9, 96.9) for death. Similarly, the overall estimated level of protection provided by prior SARS-CoV-2 infection was 94.8% (95% CI: 94.4, 95.1) for documented infection, 94.1% (95% CI: 91.9, 95.7) for hospitalization, and 96.4% (95% CI: 92.5, 98.3) for severe illness. Our results should be considered by policy-makers when deciding whether or not to prioritize vaccination of previously infected adults.

The role of children in the spread of COVID-19: Using household data from Bnei Brak, Israel, to estimate the relative susceptibility and infectivity of children.

One of the significant unanswered questions about COVID-19 epidemiology relates to the role of children in transmission. This study uses data on infections within households in order to estimate the susceptibility and infectivity of children compared to those of adults. The data were collected from households in the city of Bnei Brak, Israel, in which all household members were tested for COVID-19 using PCR (637 households, average household size of 5.3). In addition, serological tests were performed on a subset of the individuals in the study. Inspection of the PCR data shows that children are less likely to be tested positive compared to adults (25% of children positive over all households, 44% of adults positive over all households, excluding index cases), and the chance of being positive increases with age. Analysis of joint PCR/serological data shows that there is under-detection of infections in the PCR testing, which is more substantial in children. However, the differences in detection rates are not sufficient to account for the differences in PCR positive rates in the two age groups. To estimate relative transmission parameters, we employ a discrete stochastic model of the spread of infection within a household, allowing for susceptibility and infectivity parameters to differ among children and adults. The model is fitted to the household data using a simulated maximum likelihood approach. To adjust parameter estimates for under-detection of infections in the PCR results, we employ a multiple imputation procedure using estimates of under-detection in children and adults, based on the available serological data. We estimate that the susceptibility of children (under 20 years old) is 43% (95% CI: [31%, 55%]) of the susceptibility of adults. The infectivity of children was estimated to be 63% (95% CI: [37%, 88%]) relative to that of adults.

A two-stage approach for estimating the parameters of an age-group epidemic model from incidence data.

Age-dependent dynamics is an important characteristic of many infectious diseases. Age-group epidemic models describe the infection dynamics in different age-groups by allowing to set distinct parameter values for each. However, such models are highly nonlinear and may have a large number of unknown parameters. Thus, parameter estimation of age-group models, while becoming a fundamental issue for both the scientific study and policy making in infectious diseases, is not a trivial task in practice. In this paper, we examine the estimation of the so-called next-generation matrix using incidence data of a single entire outbreak, and extend the approach to deal with recurring outbreaks. Unlike previous studies, we do not assume any constraints regarding the structure of the matrix. A novel two-stage approach is developed, which allows for efficient parameter estimation from both statistical and computational perspectives. Simulation studies corroborate the ability to estimate accurately the parameters of the model for several realistic scenarios. The model and estimation method are applied to real data of influenza-like-illness in Israel. The parameter estimates of the key relevant epidemiological parameters and the recovered structure of the estimated next-generation matrix are in line with results obtained in previous studies.

Characterizing risk of Ebola transmission based on frequency and type of case-contact exposures.

During the initial months of the 2013-2016 Ebola epidemic, rapid geographical dissemination and intense transmission challenged response efforts across West Africa. Contextual behaviours associated with increased risk of exposure included travel to high-transmission settings, caring for sick and preparing the deceased for traditional funerals. Although such behaviours are widespread in West Africa, high-transmission pockets were observed. Superspreading and clustering are typical phenomena in infectious disease outbreaks, as a relatively small number of transmission chains are often responsible for the majority of events. Determining the characteristics of contacts at greatest risk of developing disease and of cases with greatest transmission potential could therefore help curb propagation of infection. Our analysis of contact tracing data from Montserrado County, Liberia, suggested that the probability of transmission was 4.5 times higher for individuals who were reported as having contact with multiple cases. The probability of individuals developing disease was not significantly associated with age or sex of their source case but was higher when they were in the same household as the infectious case. Surveillance efforts for rapidly identifying symptomatic individuals and effectively messaged campaigns encouraging household members to bring the sick to designated treatment centres without administration of home care could mitigate transmission.This article is part of the themed issue 'The 2013-2016 West African Ebola epidemic: data, decision-making and disease control'.

Modeling the spread of polio in an IPV-vaccinated population: lessons learned from the 2013 silent outbreak in southern Israel.

BACKGROUND: Polio eradication is an extraordinary globally coordinated health program in terms of its magnitude and reach, leading to the elimination of wild poliovirus (WPV) in most parts of the world. In 2013, a silent outbreak of WPV was detected in Israel, a country using an inactivated polio vaccine (IPV) exclusively since 2005. The outbreak was detected using environmental surveillance (ES) of sewage reservoirs. Stool surveys indicated the outbreak to be restricted mainly to children under the age of 10 in the Bedouin population of southern Israel. In order to curtail the outbreak, a nationwide vaccination campaign using oral polio vaccine (OPV) was conducted, targeting all children under 10. METHODS: A transmission model, fitted to the results of the stool surveys, with additional conditions set by the ES measurements, was used to evaluate the prevalence of WPV in Bedouin children and the effectiveness of the vaccination campaign. Employing the parameter estimates of the model fitting, the model was used to investigate the effect of alternative timings, coverages and dosages of the OPV campaign on the outcome of the outbreak. RESULTS: The mean estimate for the mean reproductive number was 1.77 (95 % credible interval, 1.46-2.30). With seasonal variation, the reproductive number maximum range was between zero and six. The mean estimate for the mean infectious periods was 16.8 (8.6-24.9) days. The modeling indicates the OPV campaign was effective in curtailing the outbreak. The mean estimate for the attack rate in Bedouin children under 10 at the end of 2014 was 42 % (22-65 %), whereas without the campaign the mean projected attack rate was 57 % (35-74 %). The campaign also likely shortened the duration of the outbreak by a mean estimate of 309 (2-846) days. A faster initiation of the OPV campaign could have reduced the incidence of WPV even if a lower coverage was reached, at the risk of prolonging the outbreak. CONCLUSIONS: OPV campaigns are essential for interrupting WPV transmission, even in a developed country setting with a high coverage of IPV. In this setting, establishing ES of WPV circulation is particularly crucial for early detection and containment of an outbreak.

Multiannual forecasting of seasonal influenza dynamics reveals climatic and evolutionary drivers.

Human influenza occurs annually in most temperate climatic zones of the world, with epidemics peaking in the cold winter months. Considerable debate surrounds the relative role of epidemic dynamics, viral evolution, and climatic drivers in driving year-to-year variability of outbreaks. The ultimate test of understanding is prediction; however, existing influenza models rarely forecast beyond a single year at best. Here, we use a simple epidemiological model to reveal multiannual predictability based on high-quality influenza surveillance data for Israel; the model fit is corroborated by simple metapopulation comparisons within Israel. Successful forecasts are driven by temperature, humidity, antigenic drift, and immunity loss. Essentially, influenza dynamics are a balance between large perturbations following significant antigenic jumps, interspersed with nonlinear epidemic dynamics tuned by climatic forcing.

Modeling and statistical analysis of the spatio-temporal patterns of seasonal influenza in Israel.

BACKGROUND: Seasonal influenza outbreaks are a serious burden for public health worldwide and cause morbidity to millions of people each year. In the temperate zone influenza is predominantly seasonal, with epidemics occurring every winter, but the severity of the outbreaks vary substantially between years. In this study we used a highly detailed database, which gave us both temporal and spatial information of influenza dynamics in Israel in the years 1998-2009. We use a discrete-time stochastic epidemic SIR model to find estimates and credible confidence intervals of key epidemiological parameters. FINDINGS: Despite the biological complexity of the disease we found that a simple SIR-type model can be fitted successfully to the seasonal influenza data. This was true at both the national levels and at the scale of single cities.The effective reproductive number R(e) varies between the different years both nationally and among Israeli cities. However, we did not find differences in R(e) between different Israeli cities within a year. R(e) was positively correlated to the strength of the spatial synchronization in Israel. For those years in which the disease was more "infectious", then outbreaks in different cities tended to occur with smaller time lags. Our spatial analysis demonstrates that both the timing and the strength of the outbreak within a year are highly synchronized between the Israeli cities. We extend the spatial analysis to demonstrate the existence of high synchrony between Israeli and French influenza outbreaks. CONCLUSIONS: The data analysis combined with mathematical modeling provided a better understanding of the spatio-temporal and synchronization dynamics of influenza in Israel and between Israel and France. Altogether, we show that despite major differences in demography and weather conditions intra-annual influenza epidemics are tightly synchronized in both their timing and magnitude, while they may vary greatly between years. The predominance of a similar main strain of influenza, combined with population mixing serve to enhance local and global influenza synchronization within an influenza season.

Modelling seasonal influenza in Israel.

Mathematical modeling approaches are used to study the epidemic dynamics of seasonal influenza in Israel. The recent availability of highly resolved ten year timeseries of influenza cases provides an opportunity for modeling and estimating important epidemiological parameters in the Israeli population. A simple but well known SIR discrete-time deterministic model was fitted to consecutive epidemics allowing estimation of the initial number of susceptibles in the population S0, as well as the reproductive number R0 each year. The results were corroborated by implementing a stochastic model and using a maximum likelihood approach. The paper discusses the difficulties in estimating these important parameters especially when the reporting rate of influenza cases might only be known with limited accuracy, as is generally the case. In such situations invariant parameters such as the percentage of susceptibles infected, and the effective reproductive rate might be preferred, as they do not depend on reporting rate. Results are given based on the Israeli timeseries.

Onset of a pandemic: characterizing the initial phase of the swine flu (H1N1) epidemic in Israel.

BACKGROUND: The swine influenza H1N1 first identified in Mexico, spread rapidly across the globe and is considered the fastest moving pandemic in history. The early phase of an outbreak, in which data is relatively scarce, presents scientific challenges on key issues such as: scale, severity and immunity which are fundamental for establishing sound and rapid policy schemes. Our analysis of an Israeli dataset aims at understanding the spatio-temporal dynamics of H1N1 in its initial phase. METHODS: We constructed and analyzed a unique dataset from Israel on all confirmed cases (between April 26 to July 7, 2009), representing most swine flu cases in this period. We estimated and characterized fundamental epidemiological features of the pandemic in Israel (e.g. effective reproductive number, age-class distribution, at-risk social groups, infections between sexes, and spatial dynamics). Contact data collected during this stage was used to estimate the generation time distribution of the pandemic. RESULTS: We found a low effective reproductive number (Re=1.06), an age-class distribution of infected individuals (skewed towards ages 18-25), at-risk social groups (soldiers and ultra Orthodox Jews), and significant differences in infections between sexes (skewed towards males). In terms of spatial dynamics, the pandemic spread from the central coastal plain of Israel to other regions, with higher infection rates in more densely populated sub-districts with higher income households. CONCLUSIONS: Analysis of high quality data holds much promise in reducing uncertainty regarding fundamental aspects of the initial phase of an outbreak (e.g. the effective reproductive number Re, age-class distribution, at-risk social groups). The formulation for determining the effective reproductive number Re used here has many advantages for studying the initial phase of the outbreak since it neither assumes exponential growth of infectives and is independent of the reporting rate. The finding of a low Re (close to unity threshold), combined with identification of social groups with high transmission rates would have enabled the containment of swine flu during the summer in Israel. Our unique use of contact data provided new insights into the differential dynamics of influenza in different ages and sexes, and should be promoted in future epidemiological studies. Thus our work highlights the importance of conducting a comprehensive study of the initial stage of a pandemic in real time.

[Mathematical models as tools for studying and developing strategies in the case of a pandemic influenza outbreak].

The current spread of swine flu H1N1 raises serious concerns for public health worldwide. Mathematical modelling has proved to be an essential tool for both developing strategies in preparation for an outbreak and for predicting and evaluating the effectiveness of control policies during an outbreak. Given its growing importance, this article outlines some of the fundamental contributions of mathematical modelling in the study of infectious diseases. The authors review the classical SIR model which has become central to epidemiology, demonstrating basic concepts such as outbreak threshold, the reproductive number Ro and herd immunity. The authors show how the model can be expanded to include different intervention and mitigation strategies, and discuss other biological and social complexities that may be introduced. Finally, the paper illustrates different scenarios for the spread of swine flu in Israel and provides estimates for Reproductive rate (Ro).