Deploying digital health data to optimize influenza surveillance at national and local scales

Antigenic drift and subtype interference shape A(H3N2) epidemic dynamics in the United States

Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997-2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.

Antigenic drift and subtype interference shape A(H3N2) epidemic dynamics in the United States

Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997-2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection dynamics, presumably via heterosubtypic cross-immunity.

Molecular characterization of influenza virus circulating in Nepal in the year 2019

Influenza (sometimes referred to as "flu") is a contagious viral infection of the airways in the lungs that affects a significant portion of the world's population. Clinical symptoms of influenza virus infections can range widely, from severe pneumonia to moderate or even asymptomatic sickness. If left untreated, influenza can have more severe effects on the heart, brain, and lungs than on the respiratory tract and can necessitate hospitalization. This study was aimed to investigate and characterize all types of influenza cases prevailing in Nepal and to analyze seasonal occurrence of Influenza in Nepal in the year 2019. A cross sectional, retrospective and descriptive study was carried out at National Influenza Center (NIC), National Public Health Laboratory Kathmandu Nepal for the period of one year (Jan-Dec 2019). A total of 3606 throat swab samples from various age groups and sexes were processed at the NIC. The specimens were primarily stored at 4 °C and processed using ABI 7500 RT PCR system for the identification of Influenza virus types and subtypes. Data accessed for research purpose were retrieved from National Influenza Centre (NIC) on 1st Jan 2020. Of the total 3606 patients suspected of having influenza infection, influenza viruses were isolated from 1213 (33.6%) patients with male predominance. The highest number of infection was caused by Influenza A/Pdm09 strain 739 (60.9%) followed by Influenza B 304 (25.1%) and Influenza A/H3 169 (13.9%) and most remarkable finding of this study was the detection of H5N1 in human which is the first ever case of such infection in human from Nepal. Similar to other tropical nations, influenza viruses were detected year-round in various geographical locations of Nepal. The influenza virus type and subtypes that were in circulation in Nepal were comparable to vaccine candidate viruses, which the currently available influenza vaccine may prevent.

Spatial distribution and determinants of childhood vaccination refusal in the United States

Parental refusal and delay of childhood vaccination has increased in recent years in the United States. This phenomenon challenges maintenance of herd immunity and increases the risk of outbreaks of vaccine-preventable diseases. We examine US county-level vaccine refusal for patients under five years of age collected during the period 2012-2015 from an administrative healthcare dataset. We model these data with a Bayesian zero-inflated negative binomial regression model to capture social and political processes that are associated with vaccine refusal, as well as factors that affect our measurement of vaccine refusal. Our work highlights fine-scale socio-demographic characteristics associated with vaccine refusal nationally, finds that spatial clustering in refusal can be explained by such factors, and has the potential to aid in the development of targeted public health strategies for optimizing vaccine uptake.

Spatial aggregation choice in the era of digital and administrative surveillance data

Traditional disease surveillance is increasingly being complemented by data from non-traditional sources like medical claims, electronic health records, and participatory syndromic data platforms. As non-traditional data are often collected at the individual-level and are convenience samples from a population, choices must be made on the aggregation of these data for epidemiological inference. Our study seeks to understand the influence of spatial aggregation choice on our understanding of disease spread with a case study of influenza-like illness in the United States. Using U.S. medical claims data from 2002 to 2009, we examined the epidemic source location, onset and peak season timing, and epidemic duration of influenza seasons for data aggregated to the county and state scales. We also compared spatial autocorrelation and tested the relative magnitude of spatial aggregation differences between onset and peak measures of disease burden. We found discrepancies in the inferred epidemic source locations and estimated influenza season onsets and peaks when comparing county and state-level data. Spatial autocorrelation was detected across more expansive geographic ranges during the peak season as compared to the early flu season, and there were greater spatial aggregation differences in early season measures as well. Epidemiological inferences are more sensitive to spatial scale early on during U.S. influenza seasons, when there is greater heterogeneity in timing, intensity, and geographic spread of the epidemics. Users of non-traditional disease surveillance should carefully consider how to extract accurate disease signals from finer-scaled data for early use in disease outbreaks.

Addressing the socioeconomic divide in computational modeling for infectious diseases

The COVID-19 pandemic has highlighted how structural social inequities fundamentally shape disease dynamics. Here, the authors provide a set of practical and methodological recommendations to address socioeconomic vulnerabilities in epidemic models.

Network effects in influenza spread: The impact of mobility and socio-economic factors

This paper introduces new methods of modeling and analyzing social networks that emerge in the context of disease spread. Four methods of constructing informative networks are presented, two of which use. static data and two use temporal data, namely individual citizen mobility observations taken over an extensive period of time. We show how the built networks can be analyzed, and how the numerical results can be interpreted, using network permutation-based surprise analysis. In doing so, we explain the relationship of surprise analysis with conventional network hypothesis testing and Quadratic Assignment Procedure regression. Surprise analysis is more comprehensive, and can be without limitation performed with any form(s) of network subgraphs, including those with multiple nodal attributes, weighted links, and temporal features. To illustrate our methodological work in application, we put them to use for interpreting networks constructed from the data collected over one year in an observational study in Buffalo and Erie counties in New York state during the 2016-2017 influenza season. Even with the limitations in the data size, our methods are able to reveal the global (city- and season-wide) patterns in the spread of influenza, taking into account population mobility and socio-economic factors.

Toward the use of neural networks for influenza prediction at multiple spatial resolutions

Mitigating the effects of disease outbreaks with timely and effective interventions requires accurate real-time surveillance and forecasting of disease activity, but traditional health care-based surveillance systems are limited by inherent reporting delays. Machine learning methods have the potential to fill this temporal "data gap," but work to date in this area has focused on relatively simple methods and coarse geographic resolutions (state level and above). We evaluate the predictive performance of a gated recurrent unit neural network approach in comparison with baseline machine learning methods for estimating influenza activity in the United States at the state and city levels and experiment with the inclusion of real-time Internet search data. We find that the neural network approach improves upon baseline models for long time horizons of prediction but is not improved by real-time internet search data. We conduct a thorough analysis of feature importances in all considered models for interpretability purposes.

Health inequities in influenza transmission and surveillance

The lower an individual's socioeconomic position, the higher their risk of poor health in low-, middle-, and high-income settings alike. As health inequities grow, it is imperative that we develop an empirically-driven mechanistic understanding of the determinants of health disparities, and capture disease burden in at-risk populations to prevent exacerbation of disparities. Past work has been limited in data or scope and has thus fallen short of generalizable insights. Here, we integrate empirical data from observational studies and large-scale healthcare data with models to characterize the dynamics and spatial heterogeneity of health disparities in an infectious disease case study: influenza. We find that variation in social and healthcare-based determinants exacerbates influenza epidemics, and that low socioeconomic status (SES) individuals disproportionately bear the burden of infection. We also identify geographical hotspots of influenza burden in low SES populations, much of which is overlooked in traditional influenza surveillance, and find that these differences are most predicted by variation in susceptibility and access to sickness absenteeism. Our results highlight that the effect of overlapping factors is synergistic and that reducing this intersectionality can significantly reduce inequities. Additionally, health disparities are expressed geographically, and targeting public health efforts spatially may be an efficient use of resources to abate inequities. The association between health and socioeconomic prosperity has a long history in the epidemiological literature; addressing health inequities in respiratory-transmitted infectious disease burden is an important step towards social justice in public health, and ignoring them promises to pose a serious threat.

Optimizing respiratory virus surveillance networks using uncertainty propagation

Infectious disease prevention, control and forecasting rely on sentinel observations; however, many locations lack the capacity for routine surveillance. Here we show that, by using data from multiple sites collectively, accurate estimation and forecasting of respiratory diseases for locations without surveillance is feasible. We develop a framework to optimize surveillance sites that suppresses uncertainty propagation in a networked disease transmission model. Using influenza outbreaks from 35 US states, the optimized system generates better near-term predictions than alternate systems designed using population and human mobility. We also find that monitoring regional population centers serves as a reasonable proxy for the optimized network and could direct surveillance for diseases with limited records. The  proxy method is validated using model simulations for 3,108 US counties and historical data for two other respiratory pathogens - human metapneumovirus and seasonal coronavirus - from 35 US states and can be used to guide systemic allocation of surveillance efforts.

Infectious or Recovered? Optimizing the Infectious Disease Detection Process for Epidemic Control and Prevention Based on Social Media

Detecting the period of a disease is of great importance to building information management capacity in disease control and prevention. This paper aims to optimize the disease surveillance process by further identifying the infectious or recovered period of flu cases through social media. Specifically, this paper explores the potential of using public sentiment to detect flu periods at word level. At text level, we constructed a deep learning method to classify the flu period and improve the classification result with sentiment polarity. Three important findings are revealed. Firstly, bloggers in different periods express significantly different sentiments. Blogger sentiments in the recovered period are more positive than in the infectious period when measured by the interclass distance. Secondly, the optimized disease detection process can substantially improve the classification accuracy of flu periods from 0.876 to 0.926. Thirdly, our experimental results confirm that sentiment classification plays a crucial role in accuracy improvement. Precise identification of disease periods enhances the channels for the disease surveillance processes. Therefore, a disease outbreak can be predicted credibly when a larger population is monitored. The research method proposed in our work also provides decision making reference for proactive and effective epidemic control and prevention in real time.

Differences in Regional Patterns of Influenza Activity Across Surveillance Systems in the United States: Comparative Evaluation

BACKGROUND

The Centers for Disease Control and Prevention (CDC) tracks influenza-like illness (ILI) using information on patient visits to health care providers through the Outpatient Influenza-like Illness Surveillance Network (ILINet). As participation in this system is voluntary, the composition, coverage, and consistency of health care reports vary from state to state, leading to different measures of ILI activity between regions. The degree to which these measures reflect actual differences in influenza activity or systematic differences in the methods used to collect and aggregate the data is unclear.

OBJECTIVE

The objective of our study was to qualitatively and quantitatively compare national and region-specific ILI activity in the United States across 4 surveillance data sources-CDC ILINet, Flu Near You (FNY), athenahealth, and HealthTweets.org-to determine whether these data sources, commonly used as input in influenza modeling efforts, show geographical patterns that are similar to those observed in CDC ILINet's data. We also compared the yearly percentage of FNY participants who sought health care for ILI symptoms across geographical areas.

METHODS

We compared the national and regional 2018-2019 ILI activity baselines, calculated using noninfluenza weeks from previous years, for each surveillance data source. We also compared measures of ILI activity across geographical areas during 3 influenza seasons, 2015-2016, 2016-2017, and 2017-2018. Geographical differences in weekly ILI activity within each data source were also assessed using relative mean differences and time series heatmaps. National and regional age-adjusted health care-seeking percentages were calculated for each influenza season by dividing the number of FNY participants who sought medical care for ILI symptoms by the total number of ILI reports within an influenza season. Pearson correlations were used to assess the association between the health care-seeking percentages and baselines for each surveillance data source.

RESULTS

We observed consistent differences in ILI activity across geographical areas for CDC ILINet and athenahealth data. ILI activity for FNY displayed little variation across geographical areas, whereas differences in ILI activity for HealthTweets.org were associated with the total number of tweets within a geographical area. The percentage of FNY participants who sought health care for ILI symptoms differed slightly across geographical areas, and these percentages were positively correlated with CDC ILINet and athenahealth baselines.

CONCLUSIONS

Our findings suggest that differences in ILI activity across geographical areas as reported by a given surveillance system may not accurately reflect true differences in the prevalence of ILI. Instead, these differences may reflect systematic collection and aggregation biases that are particular to each system and consistent across influenza seasons. These findings are potentially relevant in the real-time analysis of the influenza season and in the definition of unbiased forecast models.

Positive selection pressure on E2 protein of classical swine fever virus drives variations in virulence, pathogenesis and antigenicity: Implication for epidemiological surveillance in endemic areas

Classical swine fever (CSF), caused by CSF virus (CSFV), is considered one of the most important infectious diseases with devasting consequences for the pig industry. Recent reports describe the emergence of new CSFV strains resulting from the action of positive selection pressure, due mainly to the bottleneck effect generated by ineffective vaccination. Even though a decrease in the genetic diversity of the positively selected CSFV strains has been observed by several research groups, there is little information about the effect of this selective force on the virulence degree, antigenicity and pathogenicity of this type of strains. Hence, the aim of the current study was to determine the effect of the positive selection pressure on these three parameters of CSFV strains, emerged as result of the bottleneck effects induced by improper vaccination in a CSF-endemic area. Moreover, the effect of the positively selected strains on the epidemiological surveillance system was assessed. By the combination of in vitro, in vivo and immunoinformatic approaches, we revealed that the action of the positive selection pressure induces a decrease in virulence and alteration in pathogenicity and antigenicity. However, we also noted that the evolutionary process of CSFV, especially in segregated microenvironments, could contribute to the gain-fitness event, restoring the highly virulent pattern of the circulating strains. Besides, we denoted that the presence of low virulent strains selected by bottleneck effect after inefficient vaccination can lead to a relevant challenge for the epidemiological surveillance of CSF, contributing to under-reports of the disease, favouring the perpetuation of the virus in the field. In this study, B-cell and CTL epitopes on the E2 3D-structure model were also identified. Thus, the current study provides novel and significant insights into variation in virulence, pathogenesis and antigenicity experienced by CSFV strains after the positive selection pressure effect.

Incidence estimation from sentinel surveillance data; a simulation study and application to data from the Belgian laboratory sentinel surveillance

Background

Inverse probability weighting (IPW) methods can be used to estimate the total number of cases from the sample collected through sentinel surveillance. Central to these methods are the inverse weights which can be derived in several ways and, in this case, represent the probability that laboratory (lab) sentinel surveillance detects a lab-confirmed case.

Methods

We compare different weights in a simulation study. Weights are obtained from the proportion of participating labs over all labs. We adjust these weights for attractiveness and density of labs over population. The market share of sentinel labs, as estimated by the econometric Huff-model, is also considered. Additionally, we investigate the effect of not recognizing sentinel labs as sentinel labs when they report no cases. We estimate the bias associated with the different weights as the difference between the simulated number of cases and the estimate of this total from the sentinel sample.

As motivating data examples, we apply an extended Huff-model to four pathogens under laboratory sentinel surveillance in Belgium between 2010 and 2015 and discuss the model fit. We estimate the total number of lab-confirmed cases associated with Rotavirus, influenza virus, Y. enterocolitica and Campylobacter spp.. The extended Huff-model takes the lab-concept, the number of reimbursements and the number of departments, lab-density, regional borders, distance and competition between labs in account.

Results

Estimates obtained with the Huff-model were most accurate in the more complex simulation scenarios as compared to other weights. In the data examples, several significant coefficients are identified, but the fit of the Huff-model to the Belgian sentinel surveillance data leaves much variability in market shares unexplained.

Conclusion

The Huff-model allows for estimation of the spatial and population coverage of sentinel surveillance and through IPW-methods also for the estimation of the total number of cases. The Huff-model‘s gravity function allows us to differentiate inside an area while estimating from the full dataset. Our data examples show that additional data on the participation to surveillance and practices of labs is necessary for a more accurate estimation.

Patterns of seasonal influenza activity in U.S. core-based statistical areas, described using prescriptions of oseltamivir in Medicare claims data

Using Medicare claims data on prescriptions of oseltamivir dispensed to people 65 years old and older, we present a descriptive analysis of patterns of influenza activity in the United States for 579 core-based statistical areas (CBSAs) from the 2010-2011 through the 2015-2016 influenza seasons. During this time, 1,010,819 beneficiaries received a prescription of oseltamivir, ranging from 45,888 in 2011-2012 to 380,745 in 2014-2015. For each season, the peak weekly number of prescriptions correlated with the total number of prescriptions (Pearson's r ≥ 0.88). The variance in peak timing decreased with increasing severity (p < 0.0001). Among these 579 CBSAs, neither peak timing, nor relative timing, nor severity of influenza seasons showed evidence of spatial autocorrelation (0.02 ≤ Moran's I ≤ 0.23). After aggregating data to the state level, agreement between the seasonal severity at the CBSA level and the state level was fair (median Cohen's weighted κ = 0.32, interquartile range = 0.26-0.39). Based on seasonal severity, relative timing, and geographic place, we used hierarchical agglomerative clustering to join CBSAs into influenza zones for each season. Seasonal maps of influenza zones showed no obvious patterns that might assist in predicting influenza zones for future seasons. Because of the large number of prescriptions, these data may be especially useful for characterizing influenza activity and geographic distribution during low severity seasons, when other data sources measuring influenza activity are likely to be sparse.

Urbanization and humidity shape the intensity of influenza epidemics in U.S. cities

Influenza epidemics vary in intensity from year to year, driven by climatic conditions and by viral antigenic evolution. However, important spatial variation remains unexplained. Here we show predictable differences in influenza incidence among cities, driven by population size and structure. Weekly incidence data from 603 cities in the United States reveal that epidemics in smaller cities are focused on shorter periods of the influenza season, whereas in larger cities, incidence is more diffuse. Base transmission potential estimated from city-level incidence data is positively correlated with population size and with spatiotemporal organization in population density, indicating a milder response to climate forcing in metropolises. This suggests that urban centers incubate critical chains of transmission outside of peak climatic conditions, altering the spatiotemporal geometry of herd immunity.