Search

The pathogenic bacteria Neisseria meningitidis, which causes invasive meningococcal disease (IMD), predominantly colonizes humans asymptomatically; however, invasive disease occurs in a small proportion of the population. Here, we explore the seasonality of IMD and develop and validate a suite of models for simulating and forecasting disease outcomes in the United States. We combine the models into multi-model ensembles (MME) based on the past performance of the individual models, as well as a naive equally weighted aggregation, and compare the retrospective forecast performance over a six-month forecast horizon. Deployment of the complete vaccination regimen, introduced in 2011, coincided with a change in the periodicity of IMD, suggesting altered transmission dynamics. We found that a model forced with the period obtained by local power wavelet decomposition best fit and forecast observations. In addition, the MME performed the best across the entire study period. Finally, our study included US-level data until 2022, allowing study of a possible IMD rebound after relaxation of non-pharmaceutical interventions imposed in response to the COVID-19 pandemic; however, no evidence of a rebound was found. Our findings demonstrate the ability of process-based models to retrospectively forecast IMD and provide a first analysis of the seasonality of IMD before and after the complete vaccination regimen. Author summary: This paper presents a time series analysis of Invasive Meningococcal Disease caused by the pathogenic bacteria Neisseria meningitidis as well as the development and validation of a suite of models for simulating and forecasting disease outcomes in the United States. We present the models and their performance and construct a multi-model ensemble (MME) forecast that combines individual models based on past forecast performance, as well as an equally weighted aggregation. Deployment of a complete vaccination regimen, introduced for adolescents in 2011, coincided with a change in the periodicity of IMD, suggesting altered transmission dynamics. We found that a model forced with the period obtained by local power wavelet decomposition best fits and forecasts observations. In addition, the MME forecasts performed the best across the entire study period. Finally, our study included US-level data until 2022, allowing study of a possible IMD rebound after the relaxation of non-pharmaceutical interventions imposed in response to the COVID-19 pandemic; however, no evidence of a rebound was found. Our findings demonstrate the ability of process-based models to retrospectively forecast IMD and provide a first analysis of the seasonality of IMD before and after the complete vaccination regimen. [ABSTRACT FROM AUTHOR]

Influenza affects millions of people every year. It causes a considerable amount of medical visits and hospitalisations as well as hundreds of thousands of deaths. Forecasting influenza prevalence with good accuracy can significantly help public health agencies to timely react to seasonal or novel strain epidemics. Although significant progress has been made, influenza forecasting remains a challenging modelling task. In this paper, we propose a methodological framework that improves over the state-of-the-art forecasting accuracy of influenza-like illness (ILI) rates in the United States. We achieve this by using Web search activity time series in conjunction with historical ILI rates as observations for training neural network (NN) architectures. The proposed models incorporate Bayesian layers to produce associated uncertainty intervals to their forecast estimates, positioning themselves as legitimate complementary solutions to more conventional approaches. The best performing NN, referred to as the iterative recurrent neural network (IRNN) architecture, reduces mean absolute error by 10.3% and improves skill by 17.1% on average in nowcasting and forecasting tasks across 4 consecutive flu seasons. Author summary: Modelling the prevalence of an infectious disease enables public health organisations to prepare for and minimise its impact. However, traditional disease indicators are often quite restrictive as they provide information with a significant delay. Recent research efforts have provided evidence of the value of alternative information sources such as Web search activity trends. Our work incorporates this information into influenza forecasting models to achieve state-of-the-art accuracy. In addition, the proposed Bayesian neural network architectures also provide associated uncertainty estimates for the forecasts, positioning our methodology as a practical complementary tool for disease surveillance and policy making. [ABSTRACT FROM AUTHOR]

Measles is one the best-documented and most-mechanistically-studied non-linear infectious disease dynamical systems. However, systematic investigation into the comparative performance of traditional mechanistic models and machine learning approaches in forecasting the transmission dynamics of this pathogen are still rare. Here, we compare one of the most widely used semi-mechanistic models for measles (TSIR) with a commonly used machine learning approach (LASSO), comparing performance and limits in predicting short to long term outbreak trajectories and seasonality for both regular and less regular measles outbreaks in England and Wales (E&W) and the United States. First, our results indicate that the proposed LASSO model can efficiently use data from multiple major cities and achieve similar short-to-medium term forecasting performance to semi-mechanistic models for E&W epidemics. Second, interestingly, the LASSO model also captures annual to biennial bifurcation of measles epidemics in E&W caused by susceptible response to the late 1940s baby boom. LASSO may also outperform TSIR for predicting less-regular dynamics such as those observed in major cities in US between 1932–45. Although both approaches capture short-term forecasts, accuracy suffers for both methods as we attempt longer-term predictions in highly irregular, post-vaccination outbreaks in E&W. Finally, we illustrate that the LASSO model can both qualitatively and quantitatively reconstruct mechanistic assumptions, notably susceptible dynamics, in the TSIR model. Our results characterize the limits of predictability of infectious disease dynamics for strongly immunizing pathogens with both mechanistic and machine learning models, and identify connections between these two approaches. Author summary: Machine learning techniques in infectious disease modeling have grown in popularity in recent years. However, systematic investigation into the comparative performance of these approaches with traditional mechanistic models are still rare. In this paper, we compare one of the most widely used semi-mechanistic models for measles (TSIR) with a commonly used machine learning approach (LASSO), comparing performance and limits in predicting short to long term outbreaks of measles, one of the best-documented and most-mechanistically-studied non-linear infectious disease dynamical systems. Our results show that in general the LASSO outperform TSIR for predicting less-regular dynamics, and it can achieve similar performance in other scenarios when compared to the TSIR. The LASSO also has the advantages of not requiring explicit demographic data in model training. Finally, we identify connections between these two approaches and show that the LASSO model can both qualitatively and quantitatively reconstruct mechanistic assumptions, notably susceptible dynamics, in the TSIR model. [ABSTRACT FROM AUTHOR]

Reported COVID-19 cases and deaths provide a delayed and incomplete picture of SARS-CoV-2 infections in the United States (US). Accurate estimates of both the timing and magnitude of infections are needed to characterize viral transmission dynamics and better understand COVID-19 disease burden. We estimated time trends in SARS-CoV-2 transmission and other COVID-19 outcomes for every county in the US, from the first reported COVID-19 case in January 13, 2020 through January 1, 2021. To do so we employed a Bayesian modeling approach that explicitly accounts for reporting delays and variation in case ascertainment, and generates daily estimates of incident SARS-CoV-2 infections on the basis of reported COVID-19 cases and deaths. The model is freely available as the covidestim R package. Nationally, we estimated there had been 49 million symptomatic COVID-19 cases and 404,214 COVID-19 deaths by the end of 2020, and that 28% of the US population had been infected. There was county-level variability in the timing and magnitude of incidence, with local epidemiological trends differing substantially from state or regional averages, leading to large differences in the estimated proportion of the population infected by the end of 2020. Our estimates of true COVID-19 related deaths are consistent with independent estimates of excess mortality, and our estimated trends in cumulative incidence of SARS-CoV-2 infection are consistent with trends in seroprevalence estimates from available antibody testing studies. Reconstructing the underlying incidence of SARS-CoV-2 infections across US counties allows for a more granular understanding of disease trends and the potential impact of epidemiological drivers. Author summary: Because many COVID-19 infections go undetected, reported numbers of COVID-19 cases and deaths underestimate the true size of the epidemic. To address this problem, we built a model to estimate the number of new SARS-CoV-2 infections over time in each U.S. state and county. In this paper, we present time trends of infections and other disease outcomes from the first reported case in the U.S. until January 1, 2021, for each state and county. The time series of infection estimates suggest that the US epidemic is best described a series of related epidemics, varying in their timing and intensity. We estimate that over a quarter of the US population was infected with SARS-CoV-2 in 2020 and by the end of 2020 0.12% of the population had died from COVID-19. State-level results were consistent with external measures of disease burden, including estimates of SARS-CoV-2 seroprevalence and excess mortality. Our findings help better understand the epidemic in the pre-vaccine era and demonstrate the feasibility of estimating SARS-CoV-2 infections at local levels using routinely reported case and death data. [ABSTRACT FROM AUTHOR]

The ability to produce timely and accurate flu forecasts in the United States can significantly impact public health. Augmenting forecasts with internet data has shown promise for improving forecast accuracy and timeliness in controlled settings, but results in practice are less convincing, as models augmented with internet data have not consistently outperformed models without internet data. In this paper, we perform a controlled experiment, taking into account data backfill, to improve clarity on the benefits and limitations of augmenting an already good flu forecasting model with internet-based nowcasts. Our results show that a good flu forecasting model can benefit from the augmentation of internet-based nowcasts in practice for all considered public health-relevant forecasting targets. The degree of forecast improvement due to nowcasting, however, is uneven across forecasting targets, with short-term forecasting targets seeing the largest improvements and seasonal targets such as the peak timing and intensity seeing relatively marginal improvements. The uneven forecasting improvements across targets hold even when “perfect” nowcasts are used. These findings suggest that further improvements to flu forecasting, particularly seasonal targets, will need to derive from other, non-nowcasting approaches. [ABSTRACT FROM AUTHOR]

Prophylactic interventions such as vaccine allocation are some of the most effective public health policy planning tools. The supply of vaccines, however, is limited and an important challenge is to optimally allocate the vaccines to minimize epidemic impact. This resource allocation question (which we refer to as VID) has multiple dimensions: when, where, to whom, etc. Most of the existing literature in this topic deals with the latter (to whom), proposing policies that prioritize individuals by age and disease risk. However, since seasonal influenza spread has a typical spatial trend, and due to the temporal constraints enforced by the availability schedule, the when and where problems become equally, if not more, relevant. In this paper, we study the VID problem in the context of seasonal influenza spread in the United States. We develop a national scale metapopulation model for influenza that integrates both short and long distance human mobility, along with realistic data on vaccine uptake. We also design GA, a greedy algorithm for allocating the vaccine supply at the state level under temporal constraints and show that such a strategy improves over the current baseline of pro-rata allocation, and the improvement is more pronounced for higher vaccine efficacy and moderate flu season intensity. Further, the resulting strategy resembles a ring vaccination applied spatiallyacross the US. [ABSTRACT FROM AUTHOR]

Ancestral state reconstruction in Bayesian phylogeography of virus pandemics have been improved by utilizing a Bayesian stochastic search variable selection (BSSVS) framework. Recently, this framework has been extended to model the transition rate matrix between discrete states as a generalized linear model (GLM) of genetic, geographic, demographic, and environmental predictors of interest to the virus and incorporating BSSVS to estimate the posterior inclusion probabilities of each predictor. Although the latter appears to enhance the biological validity of ancestral state reconstruction, there has yet to be a comparison of phylogenies created by the two methods. In this paper, we compare these two methods, while also using a primitive method without BSSVS, and highlight the differences in phylogenies created by each. We test six coalescent priors and six random sequence samples of H3N2 influenza during the 2014–15 flu season in the U.S. We show that the GLMs yield significantly greater root state posterior probabilities than the two alternative methods under five of the six priors, and significantly greater Kullback-Leibler divergence values than the two alternative methods under all priors. Furthermore, the GLMs strongly implicate temperature and precipitation as driving forces of this flu season and nearly unanimously identified a single root state, which exhibits the most tropical climate during a typical flu season in the U.S. The GLM, however, appears to be highly susceptible to sampling bias compared with the other methods, which casts doubt on whether its reconstructions should be favored over those created by alternate methods. We report that a BSSVS approach with a Poisson prior demonstrates less bias toward sample size under certain conditions than the GLMs or primitive models, and believe that the connection between reconstruction method and sampling bias warrants further investigation. [ABSTRACT FROM AUTHOR]

Background: Treatment of surgical pain is a common reason for opioid prescriptions. Being able to predict which patients are at risk for opioid abuse, dependence, and overdose (opioid-related adverse outcomes [OR-AE]) could help physicians make safer prescription decisions. We aimed to develop a machine-learning algorithm to predict the risk of OR-AE following surgery using Medicaid data with external validation across states. Methods: Five machine learning models were developed and validated across seven US states (90–10 data split). The model output was the risk of OR-AE 6-months following surgery. The models were evaluated using standard metrics and area under the receiver operating characteristic curve (AUC) was used for model comparison. We assessed calibration for the top performing model and generated bootstrap estimations for standard deviations. Decision curves were generated for the top-performing model and logistic regression. Results: We evaluated 96,974 surgical patients aged 15 and 64. During the 6-month period following surgery, 10,464 (10.8%) patients had an OR-AE. Outcome rates were significantly higher for patients with depression (17.5%), diabetes (13.1%) or obesity (11.1%). The random forest model achieved the best predictive performance (AUC: 0.877; F1-score: 0.57; recall: 0.69; precision:0.48). An opioid disorder diagnosis prior to surgery was the most important feature for the model, which was well calibrated and had good discrimination. Conclusions: A machine learning models to predict risk of OR-AE following surgery performed well in external validation. This work could be used to assist pain management following surgery for Medicaid beneficiaries and supports a precision medicine approach to opioid prescribing. Author summary: The opioid epidemic continues in the United States and prescription opioids contribute to opioid-related overdose deaths. The study aimed to develop a machine-learning algorithm to predict the risk of opioid-related adverse outcomes (OR-AE) following surgery using Medicaid data with external validation across seven US states. This retrospective study of 96,974 patients presents seven machine learning models to predict opioid-related adverse outcomes (OAO) following surgery. The best predictive performance was achieved by the random forest model, with an AUC of 0.877. Patients with depression, diabetes, or obesity had higher rates of OR-AE. The model was well calibrated and had good discrimination, with an opioid disorder diagnosis prior to surgery being the most important feature. This work provides evidence that an automated risk prediction can potentially guide postoperative opiate therapy and monitoring. [ABSTRACT FROM AUTHOR]

Tumor heterogeneity is a large obstacle for cancer study and treatment. Different cancer patients may involve different combinations of gene mutations or the distinct regulatory pathways for inducing the progression of tumor. Investigating the pathways of gene mutations which can cause the formation of tumor can provide a basis for the personalized treatment of cancer. Studies suggested that KRAS, APC and TP53 are the most significant driver genes for colorectal cancer. However, it is still an open issue regarding the detailed mutation order of these genes in the development of colorectal cancer. For this purpose, we analyze the mathematical model considering all orders of mutations in oncogene, KRAS and tumor suppressor genes, APC and TP53, and fit it on data describing the incidence rates of colorectal cancer at different age from the Surveillance Epidemiology and End Results registry in the United States for the year 1973–2013. The specific orders that can induce the development of colorectal cancer are identified by the model fitting. The fitting results indicate that the mutation order with KRAS → APC → TP53, APC → TP53 → KRAS and APC → KRAS → TP53 explain the age–specific risk of colorectal cancer with very well. Furthermore, eleven pathways of gene mutations can be accepted for the mutation order of genes with KRAS → APC → TP53, APC → TP53 → KRAS and APC → KRAS → TP53, and the alternation of APC acts as the initiating or promoting event in the colorectal cancer. The estimated mutation rates of cells in the different pathways demonstrate that genetic instability must exist in colorectal cancer with alterations of genes, KRAS, APC and TP53. Author summary: Cumulative mutations in driver genes are the essential cause of cancer disease. For the colorectal cancer, KRAS, APC and TP53 are the common driver genes, and approximately 15% patients with colorectal cancer carry all mutations of the three genes. Exploring the pathway of mutations in these gene is extremely useful for the diagnosis and treatment of cancer. However, the mutation orders of these genes may vary in different patients due to the heterogeneity of tumor. Hence, we discuss all possible mutation orders in the genes, KRAS, APC and TP53 by using the model with five hits and find out the mutation pathways of genes that can effectively fit the incidence rate of colorectal cancer at different age in this article. In addition, we give the estimated values of mutation rates in each pathway that can explain the procession of colorectal cancer. The results obtained can offer guidelines to the treatment strategy of colorectal cancer. [ABSTRACT FROM AUTHOR]

Population contact patterns fundamentally determine the spread of directly transmitted airborne pathogens such as SARS-CoV-2 and influenza. Reliable quantitative estimates of contact patterns are therefore critical to modeling and reducing the spread of directly transmitted infectious diseases and to assessing the effectiveness of interventions intended to limit risky contacts. While many countries have used surveys and contact diaries to collect national-level contact data, local-level estimates of age-specific contact patterns remain rare. Yet, these local-level data are critical since disease dynamics and public health policy typically vary by geography. To overcome this challenge, we introduce a flexible model that can estimate age-specific contact patterns at the subnational level by combining national-level interpersonal contact data with other locality-specific data sources using multilevel regression with poststratification (MRP). We estimate daily contact matrices for all 50 US states and Washington DC from April 2020 to May 2021 using national contact data from the US. Our results reveal important state-level heterogeneities in levels and trends of contacts across the US over the course of the COVID-19 pandemic, with implications for the spread of respiratory diseases. Author summary: Reliable quantitative estimates of contact patterns are in high demand because they are a critical input for modeling and reducing the spread of directly transmitted infectious diseases such as COVID-19 and influenza. But, while national-level contact data are available from surveys collected in many countries—including the United States—local-level estimates of age-specific contact patterns are almost never available. Yet these local-level contact patterns are crucial, since disease dynamics and public health policy vary by geography. Here, we introduce a new method for estimating levels and trends in local-level rates of disease-relevant contacts by combining information from national-level contact data and other locality-specific data sources. We illustrate our method by estimating novel state-level estimates of contact patterns in the US over the course of the COVID-19 pandemic. These estimates are an important resource for researchers and policymakers who wish to understand the state-level dynamics of the COVID-19 pandemic in the United States. The new methods we introduce are applicable to any directly transmitted infectious disease and can be applied in many different settings all over the world. [ABSTRACT FROM AUTHOR]

Schizophrenia is a debilitating psychiatric disorder, leading to both physical and social morbidity. Worldwide 1% of the population is struggling with the disease, with 100,000 new cases annually only in the United States. Despite its importance, the goal of finding effective treatments for schizophrenia remains a challenging task, and previous work conducted expensive large-scale phenotypic screens. This work investigates the benefits of Machine Learning for graphs to optimize drug phenotypic screens and predict compounds that mitigate abnormal brain reduction induced by excessive glial phagocytic activity in schizophrenia subjects. Given a compound and its concentration as input, we propose a method that predicts a score associated with three possible compound effects, i.e., reduce, increase, or not influence phagocytosis. We leverage a high-throughput screening to prove experimentally that our method achieves good generalization capabilities. The screening involves 2218 compounds at five different concentrations. Then, we analyze the usability of our approach in a practical setting, i.e., prioritizing the selection of compounds in the SWEETLEAD library. We provide a list of 64 compounds from the library that have the most potential clinical utility for glial phagocytosis mitigation. Lastly, we propose a novel approach to computationally validate their utility as possible therapies for schizophrenia. Author summary: Phagocytosis is a fundamental biological process to protect biological organisms from exogenous infectious particles as well as to preserve equilibrium and efficiency of the host by removing its unwanted cells. A dysregulation of the phagocytic activity can lead to severe consequences for the host. In this study, we focus on a recent theory that relates an excessive phagocytic activity in brain cells, and a consequent abnormal reduction in brain volume, to the development of schizophrenia. Our working hypothesis is that pharmaceutical compounds that can reduce excessive of phagocytic activity might prove effective as a schizophrenia treatment. Rather than attempting to develop ex-novo such a chemical compound, we rely on a more cost-effective and efficient approach that seeks candidate therapies in a set of approved chemical compounds. To achieve this, we train a machine learning model capable of predicting, with good accuracy, the ability of a molecular compound to increase or decrease phagocytosis in the target brain cells. Our approach leverages learning models capable of directly processing the molecular graph of the compound, leading to the identification of 64 candidate drugs of potential clinical utility. [ABSTRACT FROM AUTHOR]

We present artificial neural networks as a feasible replacement for a mechanistic model of mosquito abundance. We develop a feed-forward neural network, a long short-term memory recurrent neural network, and a gated recurrent unit network. We evaluate the networks in their ability to replicate the spatiotemporal features of mosquito populations predicted by the mechanistic model, and discuss how augmenting the training data with time series that emphasize specific dynamical behaviors affects model performance. We conclude with an outlook on how such equation-free models may facilitate vector control or the estimation of disease risk at arbitrary spatial scales. Author summary: Aedes aegypti mosquitoes affect millions of people each year through infectious diseases such as chikungunya, dengue, and Zika. Because local vector levels need to be sufficiently high for associated outbreaks to occur, the ability to estimate mosquito abundance is a central component of assessing disease risk. The mosquito landscape model (MoLS) is a mechanistic model that estimates Aedes aegypti abundance from local weather time series, and is able to reproduce trends observed in surveillance data. However, scaling this up to a large number of locations is resource intensive, requiring a high-performance computing system. In this article, we develop artificial neural network models that are significantly faster than MoLS and can produce abundance estimates directly from local weather data. This approach reduces the computational time associated with estimating local mosquito levels, thereby allowing for a corresponding increase in the spatiotemporal resolution of these predictions. We compare network design choices, including architecture and training data, in their ability to accurately reproduce MoLS estimates and analyze model performance in locations across the contiguous United States. [ABSTRACT FROM AUTHOR]

Broad scale remote sensing promises to build forest inventories at unprecedented scales. A crucial step in this process is to associate sensor data into individual crowns. While dozens of crown detection algorithms have been proposed, their performance is typically not compared based on standard data or evaluation metrics. There is a need for a benchmark dataset to minimize differences in reported results as well as support evaluation of algorithms across a broad range of forest types. Combining RGB, LiDAR and hyperspectral sensor data from the USA National Ecological Observatory Network's Airborne Observation Platform with multiple types of evaluation data, we created a benchmark dataset to assess crown detection and delineation methods for canopy trees covering dominant forest types in the United States. This benchmark dataset includes an R package to standardize evaluation metrics and simplify comparisons between methods. The benchmark dataset contains over 6,000 image-annotated crowns, 400 field-annotated crowns, and 3,000 canopy stem points from a wide range of forest types. In addition, we include over 10,000 training crowns for optional use. We discuss the different evaluation data sources and assess the accuracy of the image-annotated crowns by comparing annotations among multiple annotators as well as overlapping field-annotated crowns. We provide an example submission and score for an open-source algorithm that can serve as a baseline for future methods. Author summary: Combining RGB, LiDAR and hyperspectral sensor data from the USA National Ecological Observatory Network's Airborne Observation Platform with multiple types of evaluation data, we created a benchmark dataset to assess crown detection and delineation methods for canopy trees covering dominant forest types in the United States. This benchmark dataset includes an R package to standardize evaluation metrics and simplify comparisons between methods. The benchmark dataset contains over 6,000 image-annotated crowns, 400 field-annotated crowns, and 3,000 canopy stem points from a wide range of forest types. [ABSTRACT FROM AUTHOR]

Effectively designing and evaluating public health responses to the ongoing COVID-19 pandemic requires accurate estimation of the prevalence of COVID-19 across the United States (US). Equipment shortages and varying testing capabilities have however hindered the usefulness of the official reported positive COVID-19 case counts. We introduce four complementary approaches to estimate the cumulative incidence of symptomatic COVID-19 in each state in the US as well as Puerto Rico and the District of Columbia, using a combination of excess influenza-like illness reports, COVID-19 test statistics, COVID-19 mortality reports, and a spatially structured epidemic model. Instead of relying on the estimate from a single data source or method that may be biased, we provide multiple estimates, each relying on different assumptions and data sources. Across our four approaches emerges the consistent conclusion that on April 4, 2020, the estimated case count was 5 to 50 times higher than the official positive test counts across the different states. Nationally, our estimates of COVID-19 symptomatic cases as of April 4 have a likely range of 2.3 to 4.8 million, with possibly as many as 7.6 million cases, up to 25 times greater than the cumulative confirmed cases of about 311,000. Extending our methods to May 16, 2020, we estimate that cumulative symptomatic incidence ranges from 4.9 to 10.1 million, as opposed to 1.5 million positive test counts. The proposed combination of approaches may prove useful in assessing the burden of COVID-19 during resurgences in the US and other countries with comparable surveillance systems. Author summary: Accurate estimates of the weekly incidence of COVID-19 in the United States is essential for planning and researching effective public health responses. Because of systematic testing shortages across the United States, official positive COVID-19 test counts are an unreliable indicator of true incidence. In this study, we present four alternative approaches for estimating cumulative incidence, which leverage different data sources and assumptions. Nationally, our estimates of COVID-19 symptomatic cases as of April 4 have a likely range of 2.3 to 4.8 million, with possibly as many as 7.6 million cases, up to 25 times greater than the cumulative confirmed cases of about 311,000. We emphasize that comparing multiple models rather than relying on a single method gives more reliable estimates of COVID-19 incidence. Our approaches could be useful for tracking the resurgence of COVID-19 in the United States as well as in other countries. [ABSTRACT FROM AUTHOR]

With an estimated $10.4 billion in medical costs and 31.4 million outpatient visits each year, influenza poses a serious burden of disease in the United States. To provide insights and advance warning into the spread of influenza, the U.S. Centers for Disease Control and Prevention (CDC) runs a challenge for forecasting weighted influenza-like illness (wILI) at the national and regional level. Many models produce independent forecasts for each geographical unit, ignoring the constraint that the national wILI is a weighted sum of regional wILI, where the weights correspond to the population size of the region. We propose a novel algorithm that transforms a set of independent forecast distributions to obey this constraint, which we refer to as probabilistically coherent. Enforcing probabilistic coherence led to an increase in forecast skill for 79% of the models we tested over multiple flu seasons, highlighting the importance of respecting the forecasting system's geographical hierarchy. Author summary: Seasonal influenza causes a significant public health burden nationwide. Accurate influenza forecasting may help public health officials allocate resources and plan responses to emerging outbreaks. The U.S. Centers for Disease Control and Prevention (CDC) reports influenza data at multiple geographical units, including regionally and nationally, where the national data are by construction a weighted sum of the regional data. In an effort to improve influenza forecast accuracy across all models submitted to the CDC's annual flu forecasting challenge, we examined the effect of imposing this geographical constraint on the set of independent forecasts, made publicly available by the CDC. We developed a novel method to transform forecast densities to obey the geographical constraint that respects the correlation structure between geographical units. This method showed consistent improvement across 79% of models and that held when stratified by targets and test seasons. Our method can be applied to other forecasting systems both within and outside an infectious disease context that have a geographical hierarchy. [ABSTRACT FROM AUTHOR]

Mechanical interactions between flowing and coagulated blood (thrombus) are crucial in dictating the deformation and remodeling of a thrombus after its formation in hemostasis. We propose a fully-Eulerian, three-dimensional, phase-field model of thrombus that is calibrated with existing in vitro experimental data. This phase-field model considers spatial variations in permeability and material properties within a single unified mathematical framework derived from an energy perspective, thereby allowing us to study effects of thrombus microstructure and properties on its deformation and possible release of emboli under different hemodynamic conditions. Moreover, we combine this proposed thrombus model with a particle-based model which simulates the initiation of the thrombus. The volume fraction of a thrombus obtained from the particle simulation is mapped to an input variable in the proposed phase-field thrombus model. The present work is thus the first computational study to integrate the initiation of a thrombus through platelet aggregation with its subsequent viscoelastic responses to various shear flows. This framework can be informed by clinical data and potentially be used to predict the risk of diverse thromboembolic events under physiological and pathological conditions. Author summary: Thromboembolism is associated with detachment of small thrombus pieces from the bulk in the blood vessel. These detached pieces, also known as emboli, travel through the blood flow and may block other vessels downstream, e.g., they may plug the deep veins of the leg, groin or arm, leading to venous thromboembolism (VTE). VTE is a significant cause of morbidity and mortality and it affects more than 900,000 people in the United States and result in approximately 100,000 deaths every year. Mechanical interaction between flowing blood and a thrombus is crucial in determining the deformation of the thrombus and the possibility of releasing emboli. In this study, we develop a phase-field model that is capable of describing the structural properties of a thrombus and its biomechanical properties under different blood flow conditions. Moreover, we combine this thrombus model with a particle-based model which simulates the initiation of the thrombus. This combined framework is the first computational study to simulate the development of a thrombus from platelet aggregation to its subsequent viscoelastic responses to various shear flows. Informed by clinical data, this framework can be used to predict the risk of diverse thromboembolic events under physiological and pathological conditions. [ABSTRACT FROM AUTHOR]

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. Author summary: Seasonal influenza outbreaks cause millions of infections and tens of thousands of deaths in the United States each year. Forecasting the track of an influenza season can help public health officials, business leaders, and the general public decide how to respond to an ongoing or emerging outbreak. Our team assembled over 20 unique forecasting models for seasonal influenza and combined them together into a single "ensemble" model. We made predictions of the 2017/2018 influenza season, each week sending real-time forecasts to the US Centers for Disease Control and Prevention (CDC). In the 2017/2018 influenza season, one of the largest seasonal outbreaks in the last 15 years, our ensemble model performed better on average than all individual forecast models in the ensemble. Based on results from this study, the CDC used forecasts from our ensemble model in public communication and internal reports in the subsequent 2018/2019 influenza season. [ABSTRACT FROM AUTHOR]

Aedes (Stegomyia) aegypti (L.) and Ae. (Stegomyia) albopictus (Skuse) mosquitoes can transmit dengue, chikungunya, yellow fever, and Zika viruses. Limited surveillance has led to uncertainty regarding the geographic ranges of these vectors globally, and particularly in regions at the present-day margins of habitat suitability such as the contiguous United States. Empirical habitat suitability models based on environmental conditions can augment surveillance gaps to describe the estimated potential species ranges, but model accuracy is unclear. We identified previously published regional and global habitat suitability models for Ae. aegypti (n = 6) and Ae. albopictus (n = 8) for which adequate information was available to reproduce the models for the contiguous U.S. Using a training subset of recently updated county-level surveillance records of Ae. aegypti and Ae. albopictus and records of counties conducting surveillance, we constructed accuracy-weighted, probabilistic ensemble models from these base models. To assess accuracy and uncertainty we compared individual and ensemble model predictions of species presence or absence to both training and testing data. The ensemble models were among the most accurate and also provided calibrated probabilities of presence for each species. The quantitative probabilistic framework enabled identification of areas with high uncertainty and model bias across the U.S. where improved models or additional data could be most beneficial. The results may be of immediate utility for counties considering surveillance and control programs for Ae. aegypti and Ae. albopictus. Moreover, the assessment framework can drive future efforts to provide validated quantitative estimates to support these programs at local, national, and international scales. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Bioinformatics has become an indispensable part of life science over the past 2 decades. However, bioinformatics education is not well integrated at the undergraduate level, especially in liberal arts colleges and regional universities in the United States. One significant obstacle pointed out by the Network for Integrating Bioinformatics into Life Sciences Education is the lack of faculty in the bioinformatics area. Most current life science professors did not acquire bioinformatics analysis skills during their own training. Consequently, a great number of undergraduate and graduate students do not get the chance to learn bioinformatics or computational biology skills within a structured curriculum during their education. To address this gap, we developed a module-based, week-long short course to train small college and regional university professors with essential bioinformatics skills. The bioinformatics modules were built to be adapted by the professor-trainees afterward and used in their own classes. All the course materials can be accessed at . [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Seasonal influenza epidemics offer unique opportunities to study the invasion and re-invasion waves of a pathogen in a partially immune population. Detailed patterns of spread remain elusive, however, due to lack of granular disease data. Here we model high-volume city-level medical claims data and human mobility proxies to explore the drivers of influenza spread in the US during 2002–2010. Although the speed and pathways of spread varied across seasons, seven of eight epidemics likely originated in the Southern US. Each epidemic was associated with 1–5 early long-range transmission events, half of which sparked onward transmission. Gravity model estimates indicate a sharp decay in influenza transmission with the distance between infectious and susceptible cities, consistent with spread dominated by work commutes rather than air traffic. Two early-onset seasons associated with antigenic novelty had particularly localized modes of spread, suggesting that novel strains may spread in a more localized fashion than previously anticipated. [ABSTRACT FROM AUTHOR]

Infectious diseases are one of the leading causes of morbidity and mortality around the world; thus, forecasting their impact is crucial for planning an effective response strategy. According to the Centers for Disease Control and Prevention (CDC), seasonal influenza affects 5% to 20% of the U.S. population and causes major economic impacts resulting from hospitalization and absenteeism. Understanding influenza dynamics and forecasting its impact is fundamental for developing prevention and mitigation strategies. We combine modern data assimilation methods with Wikipedia access logs and CDC influenza-like illness (ILI) reports to create a weekly forecast for seasonal influenza. The methods are applied to the 2013-2014 influenza season but are sufficiently general to forecast any disease outbreak, given incidence or case count data. We adjust the initialization and parametrization of a disease model and show that this allows us to determine systematic model bias. In addition, we provide a way to determine where the model diverges from observation and evaluate forecast accuracy. Wikipedia article access logs are shown to be highly correlated with historical ILI records and allow for accurate prediction of ILI data several weeks before it becomes available. The results show that prior to the peak of the flu season, our forecasting method produced 50% and 95% credible intervals for the 2013-2014 ILI observations that contained the actual observations for most weeks in the forecast. However, since our model does not account for re-infection or multiple strains of influenza, the tail of the epidemic is not predicted well after the peak of flu season has passed. [ABSTRACT FROM AUTHOR]

Over the past ten years the incidence of pertussis in the United States (U.S.) has risen steadily, with 2012 seeing the highest case number since 1955. There has also been a shift over the same time period in the age group reporting the largest number of cases (aside from infants), from adolescents to 7–11 year olds. We use epidemiological modelling and a large case incidence dataset to explain the upsurge. We investigate several hypotheses for the upsurge in pertussis cases by fitting a suite of dynamic epidemiological models to incidence data from the National Notifiable Disease Surveillance System (NNDSS) between 1990–2009, as well as incidence data from a variety of sources from 1950–1989. We find that: the best-fitting model is one in which vaccine efficacy and duration of protection of the acellular pertussis (aP) vaccine is lower than that of the whole-cell (wP) vaccine, (efficacy of the first three doses 80% [95% CI: 78%, 82%] versus 90% [95% CI: 87%, 94%]), increasing the rate at which disease is reported to NNDSS is not sufficient to explain the upsurge and 3) 2010–2012 disease incidence is predicted well. In this study, we use all available U.S. surveillance data to: 1) fit a set of mathematical models and determine which best explains these data and 2) determine the epidemiological and vaccine-related parameter values of this model. We find evidence of a difference in efficacy and duration of protection between the two vaccine types, wP and aP (aP efficacy and duration lower than wP). Future refinement of the model presented here will allow for an exploration of alternative vaccination strategies such as different age-spacings, further booster doses, and cocooning. [ABSTRACT FROM AUTHOR]

Infectious disease forecasting is of great interest to the public health community and policymakers, since forecasts can provide insight into disease dynamics in the near future and inform interventions. Due to delays in case reporting, however, forecasting models may often underestimate the current and future disease burden. In this paper, we propose a general framework for addressing reporting delay in disease forecasting efforts with the goal of improving forecasts. We propose strategies for leveraging either historical data on case reporting or external internet-based data to estimate the amount of reporting error. We then describe several approaches for adapting general forecasting pipelines to account for under- or over-reporting of cases. We apply these methods to address reporting delay in data on dengue fever cases in Puerto Rico from 1990 to 2009 and to reports of influenza-like illness (ILI) in the United States between 2010 and 2019. Through a simulation study, we compare method performance and evaluate robustness to assumption violations. Our results show that forecasting accuracy and prediction coverage almost always increase when correction methods are implemented to address reporting delay. Some of these methods required knowledge about the reporting error or high quality external data, which may not always be available. Provided alternatives include excluding recently-reported data and performing sensitivity analysis. This work provides intuition and guidance for handling delay in disease case reporting and may serve as a useful resource to inform practical infectious disease forecasting efforts., Competing Interests: The authors have declared that no competing interests exist.

Predictions of COVID-19 case growth and mortality are critical to the decisions of political leaders, businesses, and individuals grappling with the pandemic. This predictive task is challenging due to the novelty of the virus, limited data, and dynamic political and societal responses. We embed a Bayesian time series model and a random forest algorithm within an epidemiological compartmental model for empirically grounded COVID-19 predictions. The Bayesian case model fits a location-specific curve to the velocity (first derivative) of the log transformed cumulative case count, borrowing strength across geographic locations and incorporating prior information to obtain a posterior distribution for case trajectories. The compartmental model uses this distribution and predicts deaths using a random forest algorithm trained on COVID-19 data and population-level characteristics, yielding daily projections and interval estimates for cases and deaths in U.S. states. We evaluated the model by training it on progressively longer periods of the pandemic and computing its predictive accuracy over 21-day forecasts. The substantial variation in predicted trajectories and associated uncertainty between states is illustrated by comparing three unique locations: New York, Colorado, and West Virginia. The sophistication and accuracy of this COVID-19 model offer reliable predictions and uncertainty estimates for the current trajectory of the pandemic in the U.S. and provide a platform for future predictions as shifting political and societal responses alter its course., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: GLW, JAZ, PS, TB, and report personal fees from Private Health Management during the conduct of the study. CMR reports grants and personal fees from Private Health Management. MAS reports grants from US National Institutes of Health, grants from IQVIA, personal fees from Janssen Research and Development, and personal fees from Private Health Management during the conduct of the study. DX, LZ, JS, and AWR declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Dengue and influenza-like illness (ILI) are two of the leading causes of viral infection in the world and it is estimated that more than half the world's population is at risk for developing these infections. It is therefore important to develop accurate methods for forecasting dengue and ILI incidences. Since data from multiple sources (such as dengue and ILI case counts, electronic health records and frequency of multiple internet search terms from Google Trends) can improve forecasts, standard time series analysis methods are inadequate to estimate all the parameter values from the limited amount of data available if we use multiple sources. In this paper, we use a computationally efficient implementation of the known variable selection method that we call the Autoregressive Likelihood Ratio (ARLR) method. This method combines sparse representation of time series data, electronic health records data (for ILI) and Google Trends data to forecast dengue and ILI incidences. This sparse representation method uses an algorithm that maximizes an appropriate likelihood ratio at every step. Using numerical experiments, we demonstrate that our method recovers the underlying sparse model much more accurately than the lasso method. We apply our method to dengue case count data from five countries/states: Brazil, Mexico, Singapore, Taiwan, and Thailand and to ILI case count data from the United States. Numerical experiments show that our method outperforms existing time series forecasting methods in forecasting the dengue and ILI case counts. In particular, our method gives a 18 percent forecast error reduction over a leading method that also uses data from multiple sources. It also performs better than other methods in predicting the peak value of the case count and the peak time., Competing Interests: The authors have declared that no competing interests exist.