Burden and characteristics of COVID-19 in the United States during 2020.

The COVID-19 pandemic disrupted health systems and economies throughout the world during 2020 and was particularly devastating for the United States, which experienced the highest numbers of reported cases and deaths during 20201-3. Many of the epidemiological features responsible for observed rates of morbidity and mortality have been reported4-8; however, the overall burden and characteristics of COVID-19 in the United States have not been comprehensively quantified. Here we use a data-driven model-inference approach to simulate the pandemic at county-scale in the United States during 2020 and estimate critical, time-varying epidemiological properties underpinning the dynamics of the virus. The pandemic in the United States during 2020 was characterized by national ascertainment rates that increased from 11.3% (95% credible interval (CI): 8.3-15.9%) in March to 24.5% (18.6-32.3%) during December. Population susceptibility at the end of the year was 69.0% (63.6-75.4%), indicating that about one third of the US population had been infected. Community infectious rates, the percentage of people harbouring a contagious infection, increased above 0.8% (0.6-1.0%) before the end of the year, and were as high as 2.4% in some major metropolitan areas. By contrast, the infection fatality rate fell to 0.3% by year's end.

Multimodeling approach to evaluating the efficacy of layering pharmaceutical and nonpharmaceutical interventions for influenza pandemics.

When an influenza pandemic emerges, temporary school closures and antiviral treatment may slow virus spread, reduce the overall disease burden, and provide time for vaccine development, distribution, and administration while keeping a larger portion of the general population infection free. The impact of such measures will depend on the transmissibility and severity of the virus and the timing and extent of their implementation. To provide robust assessments of layered pandemic intervention strategies, the Centers for Disease Control and Prevention (CDC) funded a network of academic groups to build a framework for the development and comparison of multiple pandemic influenza models. Research teams from Columbia University, Imperial College London/Princeton University, Northeastern University, the University of Texas at Austin/Yale University, and the University of Virginia independently modeled three prescribed sets of pandemic influenza scenarios developed collaboratively by the CDC and network members. Results provided by the groups were aggregated into a mean-based ensemble. The ensemble and most component models agreed on the ranking of the most and least effective intervention strategies by impact but not on the magnitude of those impacts. In the scenarios evaluated, vaccination alone, due to the time needed for development, approval, and deployment, would not be expected to substantially reduce the numbers of illnesses, hospitalizations, and deaths that would occur. Only strategies that included early implementation of school closure were found to substantially mitigate early spread and allow time for vaccines to be developed and administered, especially under a highly transmissible pandemic scenario.

Geographic Variation, Economic Activity, and Labor Market Characteristics in Trajectories of Suicide in the United States, 2008-2020.

Suicide rates in the United States have increased over the past 15 years, with substantial geographic variation in these increases; yet there have been few attempts to cluster counties by the magnitude of suicide rate changes according to intercept and slope or to identify the economic precursors of increases. We used vital statistics data and growth mixture models to identify clusters of counties by their magnitude of suicide growth from 2008 to 2020 and examined associations with county economic and labor indices. Our models identified 5 clusters, each differentiated by intercept and slope magnitude, with the highest-rate cluster (4% of counties) being observed mainly in sparsely populated areas in the West and Alaska, starting the time series at 25.4 suicides per 100,000 population, and exhibiting the steepest increase in slope (0.69/100,000/year). There was no cluster for which the suicide rate was stable or declining. Counties in the highest-rate cluster were more likely to have agricultural and service economies and less likely to have urban professional economies. Given the increased burden of suicide, with no clusters of counties improving over time, additional policy and prevention efforts are needed, particularly targeted at rural areas in the West.

Hindcasts and forecasts of suicide mortality in US: A modeling study.

Deaths by suicide, as well as suicidal ideations, plans and attempts, have been increasing in the US for the past two decades. Deployment of effective interventions would require timely, geographically well-resolved estimates of suicide activity. In this study, we evaluated the feasibility of a two-step process for predicting suicide mortality: a) generation of hindcasts, mortality estimates for past months for which observational data would not have been available if forecasts were generated in real-time; and b) generation of forecasts with observational data augmented with hindcasts. Calls to crisis hotline services and online queries to the Google search engine for suicide-related terms were used as proxy data sources to generate hindcasts. The primary hindcast model (auto) is an Autoregressive Integrated Moving average model (ARIMA), trained on suicide mortality rates alone. Three regression models augment hindcast estimates from auto with call rates (calls), GHT search rates (ght) and both datasets together (calls_ght). The 4 forecast models used are ARIMA models trained with corresponding hindcast estimates. All models were evaluated against a baseline random walk with drift model. Rolling monthly 6-month ahead forecasts for all 50 states between 2012 and 2020 were generated. Quantile score (QS) was used to assess the quality of the forecast distributions. Median QS for auto was better than baseline (0.114 vs. 0.21. Median QS of augmented models were lower than auto, but not significantly different from each other (Wilcoxon signed-rank test, p > .05). Forecasts from augmented models were also better calibrated. Together, these results provide evidence that proxy data can address delays in release of suicide mortality data and improve forecast quality. An operational forecast system of state-level suicide risk may be feasible with sustained engagement between modelers and public health departments to appraise data sources and methods as well as to continuously evaluate forecast accuracy.

Longitudinal Association of COVID-19 Hospitalization and Death with Online Search for Loss of Smell or Taste.

Surveillance of COVID-19 is challenging but critical for mitigating disease, particularly if predictive of future disease burden. We report a robust multiyear lead-lag association between internet search activity for loss of smell or taste and COVID-19-associated hospitalization and deaths. These search data could help predict COVID-19 surges.

Correction: Recommended reporting items for epidemic forecasting and prediction research: The EPIFORGE 2020 guidelines.

[This corrects the article DOI: 10.1371/journal.pmed.1003793.].

Long-Term Fine Particulate Matter Concentrations and Prevalence of Severe Acute Respiratory Syndrome Coronavirus 2: Differential Relationships by Socioeconomic Status Among Pregnant Individuals in New York City.

We aimed to determine whether long-term ambient concentrations of fine particulate matter (particulate matter with an aerodynamic diameter less than or equal to 2.5 µm (PM2.5)) were associated with increased risk of testing positive for coronavirus disease 2019 (COVID-19) among pregnant individuals who were universally screened at delivery and whether socioeconomic status (SES) modified this relationship. We used obstetrical data collected from New-York Presbyterian Hospital/Columbia University Irving Medical Center in New York, New York, between March and December 2020, including data on Medicaid use (a proxy for low SES) and COVID-19 test results. We linked estimated 2018-2019 PM2.5 concentrations (300-m resolution) with census-tract-level population density, household size, income, and mobility (as measured by mobile-device use) on the basis of residential address. Analyses included 3,318 individuals; 5% tested positive for COVID-19 at delivery, 8% tested positive during pregnancy, and 48% used Medicaid. Average long-term PM2.5 concentrations were 7.4 (standard deviation, 0.8) µg/m3. In adjusted multilevel logistic regression models, we saw no association between PM2.5 and ever testing positive for COVID-19; however, odds were elevated among those using Medicaid (per 1-µg/m3 increase, odds ratio = 1.6, 95% confidence interval: 1.0, 2.5). Further, while only 22% of those testing positive showed symptoms, 69% of symptomatic individuals used Medicaid. SES, including unmeasured occupational exposures or increased susceptibility to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) due to concurrent social and environmental exposures, may explain the increased odds of testing positive for COVID-19 being confined to vulnerable pregnant individuals using Medicaid.

Age, period, and cohort effects on suicide death in the United States from 1999 to 2018: moderation by sex, race, and firearm involvement.

The role of sex, race, and suicide method on recent increases in suicide mortality in the United States remains unclear. Estimating the age, period, and cohort effects underlying suicide mortality trends can provide important insights for the causal hypothesis generating process. We generated updated age-period-cohort effect estimates of recent suicide mortality rates in the US, examining the putative roles of sex, race, and method for suicide, using data from all death certificates in the US between 1999 and 2018. After designating deaths as attributable to suicide according to ICD-10 underlying cause of death codes X60-X84, Y87.0, and U03, we (i) used hexagonal grids to describe rates of suicide by age, period, and cohort visually and (ii) modeled sex-, race-, and suicide method-specific age, period, and cohort effects. We found that, while suicide mortality increased in the US between 1999 and 2018 across age, sex, race, and suicide method, there was substantial heterogeneity in age and cohort effects by method, sex, and race, with a first peak of suicide risk in youth, a second peak in older ages-specific to male firearm suicide, and increased rates among younger cohorts of non-White individuals. Our findings should prompt discussion regarding age-specific clinical firearm safety interventions, drivers of minoritized populations' adverse early-life experiences, and racial differences in access to and quality of mental healthcare.

A collaborative multiyear, multimodel assessment of seasonal influenza forecasting in the United States.

Influenza infects an estimated 9-35 million individuals each year in the United States and is a contributing cause for between 12,000 and 56,000 deaths annually. Seasonal outbreaks of influenza are common in temperate regions of the world, with highest incidence typically occurring in colder and drier months of the year. Real-time forecasts of influenza transmission can inform public health response to outbreaks. We present the results of a multiinstitution collaborative effort to standardize the collection and evaluation of forecasting models for influenza in the United States for the 2010/2011 through 2016/2017 influenza seasons. For these seven seasons, we assembled weekly real-time forecasts of seven targets of public health interest from 22 different models. We compared forecast accuracy of each model relative to a historical baseline seasonal average. Across all regions of the United States, over half of the models showed consistently better performance than the historical baseline when forecasting incidence of influenza-like illness 1 wk, 2 wk, and 3 wk ahead of available data and when forecasting the timing and magnitude of the seasonal peak. In some regions, delays in data reporting were strongly and negatively associated with forecast accuracy. More timely reporting and an improved overall accessibility to novel and traditional data sources are needed to improve forecasting accuracy and its integration with real-time public health decision making.

Investigating associations between COVID-19 mortality and population-level health and socioeconomic indicators in the United States: A modeling study.

BACKGROUND: With the availability of multiple Coronavirus Disease 2019 (COVID-19) vaccines and the predicted shortages in supply for the near future, it is necessary to allocate vaccines in a manner that minimizes severe outcomes, particularly deaths. To date, vaccination strategies in the United States have focused on individual characteristics such as age and occupation. Here, we assess the utility of population-level health and socioeconomic indicators as additional criteria for geographical allocation of vaccines. METHODS AND FINDINGS: County-level estimates of 14 indicators associated with COVID-19 mortality were extracted from public data sources. Effect estimates of the individual indicators were calculated with univariate models. Presence of spatial autocorrelation was established using Moran's I statistic. Spatial simultaneous autoregressive (SAR) models that account for spatial autocorrelation in response and predictors were used to assess (i) the proportion of variance in county-level COVID-19 mortality that can explained by identified health/socioeconomic indicators (R2); and (ii) effect estimates of each predictor. Adjusting for case rates, the selected indicators individually explain 24%-29% of the variability in mortality. Prevalence of chronic kidney disease and proportion of population residing in nursing homes have the highest R2. Mortality is estimated to increase by 43 per thousand residents (95% CI: 37-49; p < 0.001) with a 1% increase in the prevalence of chronic kidney disease and by 39 deaths per thousand (95% CI: 34-44; p < 0.001) with 1% increase in population living in nursing homes. SAR models using multiple health/socioeconomic indicators explain 43% of the variability in COVID-19 mortality in US counties, adjusting for case rates. R2 was found to be not sensitive to the choice of SAR model form. Study limitations include the use of mortality rates that are not age standardized, a spatial adjacency matrix that does not capture human flows among counties, and insufficient accounting for interaction among predictors. CONCLUSIONS: Significant spatial autocorrelation exists in COVID-19 mortality in the US, and population health/socioeconomic indicators account for a considerable variability in county-level mortality. In the context of vaccine rollout in the US and globally, national and subnational estimates of burden of disease could inform optimal geographical allocation of vaccines.

Recommended reporting items for epidemic forecasting and prediction research: The EPIFORGE 2020 guidelines.

BACKGROUND: The importance of infectious disease epidemic forecasting and prediction research is underscored by decades of communicable disease outbreaks, including COVID-19. Unlike other fields of medical research, such as clinical trials and systematic reviews, no reporting guidelines exist for reporting epidemic forecasting and prediction research despite their utility. We therefore developed the EPIFORGE checklist, a guideline for standardized reporting of epidemic forecasting research. METHODS AND FINDINGS: We developed this checklist using a best-practice process for development of reporting guidelines, involving a Delphi process and broad consultation with an international panel of infectious disease modelers and model end users. The objectives of these guidelines are to improve the consistency, reproducibility, comparability, and quality of epidemic forecasting reporting. The guidelines are not designed to advise scientists on how to perform epidemic forecasting and prediction research, but rather to serve as a standard for reporting critical methodological details of such studies. CONCLUSIONS: These guidelines have been submitted to the EQUATOR network, in addition to hosting by other dedicated webpages to facilitate feedback and journal endorsement.

Suicide and the agent-host-environment triad: leveraging surveillance sources to inform prevention.

Suicide in the US has increased in the last decade, across virtually every age and demographic group. Parallel increases have occurred in non-fatal self-harm as well. Research on suicide across the world has consistently demonstrated that suicide shares many properties with a communicable disease, including person-to-person transmission and point-source outbreaks. This essay illustrates the communicable nature of suicide through analogy to basic infectious disease principles, including evidence for transmission and vulnerability through the agent-host-environment triad. We describe how mathematical modeling, a suite of epidemiological methods, which the COVID-19 pandemic has brought into renewed focus, can and should be applied to suicide in order to understand the dynamics of transmission and to forecast emerging risk areas. We describe how new and innovative sources of data, including social media and search engine data, can be used to augment traditional suicide surveillance, as well as the opportunities and challenges for modeling suicide as a communicable disease process in an effort to guide clinical and public health suicide prevention efforts.

Forecasting the spatial transmission of influenza in the United States.

Recurrent outbreaks of seasonal and pandemic influenza create a need for forecasts of the geographic spread of this pathogen. Although it is well established that the spatial progression of infection is largely attributable to human mobility, difficulty obtaining real-time information on human movement has limited its incorporation into existing infectious disease forecasting techniques. In this study, we develop and validate an ensemble forecast system for predicting the spatiotemporal spread of influenza that uses readily accessible human mobility data and a metapopulation model. In retrospective state-level forecasts for 35 US states, the system accurately predicts local influenza outbreak onset,-i.e., spatial spread, defined as the week that local incidence increases above a baseline threshold-up to 6 wk in advance of this event. In addition, the metapopulation prediction system forecasts influenza outbreak onset, peak timing, and peak intensity more accurately than isolated location-specific forecasts. The proposed framework could be applied to emergent respiratory viruses and, with appropriate modifications, other infectious diseases.

Reappraising the utility of Google Flu Trends.

Estimation of influenza-like illness (ILI) using search trends activity was intended to supplement traditional surveillance systems, and was a motivation behind the development of Google Flu Trends (GFT). However, several studies have previously reported large errors in GFT estimates of ILI in the US. Following recent release of time-stamped surveillance data, which better reflects real-time operational scenarios, we reanalyzed GFT errors. Using three data sources-GFT: an archive of weekly ILI estimates from Google Flu Trends; ILIf: fully-observed ILI rates from ILINet; and, ILIp: ILI rates available in real-time based on partial reporting-five influenza seasons were analyzed and mean square errors (MSE) of GFT and ILIp as estimates of ILIf were computed. To correct GFT errors, a random forest regression model was built with ILI and GFT rates from the previous three weeks as predictors. An overall reduction in error of 44% was observed and the errors of the corrected GFT are lower than those of ILIp. An 80% reduction in error during 2012/13, when GFT had large errors, shows that extreme failures of GFT could have been avoided. Using autoregressive integrated moving average (ARIMA) models, one- to four-week ahead forecasts were generated with two separate data streams: ILIp alone, and with both ILIp and corrected GFT. At all forecast targets and seasons, and for all but two regions, inclusion of GFT lowered MSE. Results from two alternative error measures, mean absolute error and mean absolute proportional error, were largely consistent with results from MSE. Taken together these findings provide an error profile of GFT in the US, establish strong evidence for the adoption of search trends based 'nowcasts' in influenza forecast systems, and encourage reevaluation of the utility of this data source in diverse domains.

Accuracy of real-time multi-model ensemble forecasts for seasonal influenza in the U.S.

Seasonal influenza results in substantial annual morbidity and mortality in the United States and worldwide. Accurate forecasts of key features of influenza epidemics, such as the timing and severity of the peak incidence in a given season, can inform public health response to outbreaks. As part of ongoing efforts to incorporate data and advanced analytical methods into public health decision-making, the United States Centers for Disease Control and Prevention (CDC) has organized seasonal influenza forecasting challenges since the 2013/2014 season. In the 2017/2018 season, 22 teams participated. A subset of four teams created a research consortium called the FluSight Network in early 2017. During the 2017/2018 season they worked together to produce a collaborative multi-model ensemble that combined 21 separate component models into a single model using a machine learning technique called stacking. This approach creates a weighted average of predictive densities where the weight for each component is determined by maximizing overall ensemble accuracy over past seasons. In the 2017/2018 influenza season, one of the largest seasonal outbreaks in the last 15 years, this multi-model ensemble performed better on average than all individual component models and placed second overall in the CDC challenge. It also outperformed the baseline multi-model ensemble created by the CDC that took a simple average of all models submitted to the forecasting challenge. This project shows that collaborative efforts between research teams to develop ensemble forecasting approaches can bring measurable improvements in forecast accuracy and important reductions in the variability of performance from year to year. Efforts such as this, that emphasize real-time testing and evaluation of forecasting models and facilitate the close collaboration between public health officials and modeling researchers, are essential to improving our understanding of how best to use forecasts to improve public health response to seasonal and emerging epidemic threats.

Spatiotemporal clustering of suicides in the US from 1999 to 2016: a spatial epidemiological approach.

PURPOSE: This study aims to describe and characterize the spatial and temporal clustering patterns of suicide in the ten states with the greatest suicide burden in the United States from 1999 to 2016. METHODS: All suicide deaths from January 1, 1999 to December 31, 2016 in the United States were identified using data from the Wide-ranging Online Data for Epidemiologic Research (WONDER) dataset. The ten states with the highest age-adjusted suicide rates were Montana, Alaska, Wyoming, New Mexico, Nevada, Utah, Idaho, Colorado, Arizona, and Oklahoma. A spatiotemporal scan statistic using a discrete Poisson model was employed to retrospectively detect spatiotemporal suicide clusters. RESULTS: From 1999 to 2016, a total of 649,843 suicides were recorded in the United States. Nineteen statistically significant spatiotemporal suicide mortality clusters were identified in the states with the greatest suicide rates, and 13.53% of the suicide cases within these states clustered spatiotemporally. The risk ratio of the clusters ranged from 1.45 to 3.64 (p < 0.001). All states had at least one cluster, with three clusters spanning multiple states, and four clusters were found in Arizona. While there was no clear secular trend in the average size of suicide clusters, the number of clusters increased from 1999 to 2016. CONCLUSIONS: Hot spots for suicidal behavior in the United States warrant public health intervention and continued surveillance. As suicide rates in the US continue to increase annually, public health efforts could be maximized by focusing on regions with substantial clustering.

Individual versus superensemble forecasts of seasonal influenza outbreaks in the United States.

Recent research has produced a number of methods for forecasting seasonal influenza outbreaks. However, differences among the predicted outcomes of competing forecast methods can limit their use in decision-making. Here, we present a method for reconciling these differences using Bayesian model averaging. We generated retrospective forecasts of peak timing, peak incidence, and total incidence for seasonal influenza outbreaks in 48 states and 95 cities using 21 distinct forecast methods, and combined these individual forecasts to create weighted-average superensemble forecasts. We compared the relative performance of these individual and superensemble forecast methods by geographic location, timing of forecast, and influenza season. We find that, overall, the superensemble forecasts are more accurate than any individual forecast method and less prone to producing a poor forecast. Furthermore, we find that these advantages increase when the superensemble weights are stratified according to the characteristics of the forecast or geographic location. These findings indicate that different competing influenza prediction systems can be combined into a single more accurate forecast product for operational delivery in real time.

The use of ambient humidity conditions to improve influenza forecast.

Laboratory and epidemiological evidence indicate that ambient humidity modulates the survival and transmission of influenza. Here we explore whether the inclusion of humidity forcing in mathematical models describing influenza transmission improves the accuracy of forecasts generated with those models. We generate retrospective forecasts for 95 cities over 10 seasons in the United States and assess both forecast accuracy and error. Overall, we find that humidity forcing improves forecast performance (at 1-4 lead weeks, 3.8% more peak week and 4.4% more peak intensity forecasts are accurate than with no forcing) and that forecasts generated using daily climatological humidity forcing generally outperform forecasts that utilize daily observed humidity forcing (4.4% and 2.6% respectively). These findings hold for predictions of outbreak peak intensity, peak timing, and incidence over 2- and 4-week horizons. The results indicate that use of climatological humidity forcing is warranted for current operational influenza forecast.

Type- and Subtype-Specific Influenza Forecast.

Prediction of the growth and decline of infectious disease incidence has advanced considerably in recent years. As these forecasts improve, their public health utility should increase, particularly as interventions are developed that make explicit use of forecast information. It is the task of the research community to increase the content and improve the accuracy of these infectious disease predictions. Presently, operational real-time forecasts of total influenza incidence are produced at the municipal and state level in the United States. These forecasts are generated using ensemble simulations depicting local influenza transmission dynamics, which have been optimized prior to forecast with observations of influenza incidence and data assimilation methods. Here, we explore whether forecasts targeted to predict influenza by type and subtype during 2003-2015 in the United States were more or less accurate than forecasts targeted to predict total influenza incidence. We found that forecasts separated by type/subtype generally produced more accurate predictions and, when summed, produced more accurate predictions of total influenza incidence. These findings indicate that monitoring influenza by type and subtype not only provides more detailed observational content but supports more accurate forecasting. More accurate forecasting can help officials better respond to and plan for current and future influenza activity.