Estimate of Novel Influenza A/H1N1 cases in Mexico at the early stage of the pandemic with a spatially structured epidemic model

Statistical inference using GLEaM model with spatial heterogeneity and correlation between regions

A better understanding of various patterns in the coronavirus disease 2019 (COVID-19) spread in different parts of the world is crucial to its prevention and control. Motivated by the previously developed Global Epidemic and Mobility (GLEaM) model, this paper proposes a new stochastic dynamic model to depict the evolution of COVID-19. The model allows spatial and temporal heterogeneity of transmission parameters and involves transportation between regions. Based on the proposed model, this paper also designs a two-step procedure for parameter inference, which utilizes the correlation between regions through a prior distribution that imposes graph Laplacian regularization on transmission parameters. Experiments on simulated data and real-world data in China and Europe indicate that the proposed model achieves higher accuracy in predicting the newly confirmed cases than baseline models.

Community movement and COVID-19: a global study using Google's Community Mobility Reports

Google's 'Community Mobility Reports' (CMR) detail changes in activity and mobility occurring in response to COVID-19. They thus offer the unique opportunity to examine the relationship between mobility and disease incidence. The objective was to examine whether an association between COVID-19-confirmed case numbers and levels of mobility was apparent, and if so then to examine whether such data enhance disease modelling and prediction. CMR data for countries worldwide were cross-correlated with corresponding COVID-19-confirmed case numbers. Models were fitted to explain case numbers of each country's epidemic. Models using numerical date, contemporaneous and distributed lag CMR data were contrasted using Bayesian Information Criteria. Noticeable were negative correlations between CMR data and case incidence for prominent industrialised countries of Western Europe and the North Americas. Continent-wide examination found a negative correlation for all continents with the exception of South America. When modelling, CMR-expanded models proved superior to the model without CMR. The predictions made with the distributed lag model significantly outperformed all other models. The observed relationship between CMR data and case incidence, and its ability to enhance model quality and prediction suggests data related to community mobility could prove of use in future COVID-19 modelling.

Rich do not rise early: spatio-temporal patterns in the mobility networks of different socio-economic classes

We analyse the urban mobility in the cities of Medellín and Manizales (Colombia). Each city is represented by six mobility networks, each one encoding the origin-destination trips performed by a subset of the population corresponding to a particular socio-economic status. The nodes of each network are the different urban locations whereas links account for the existence of a trip between two different areas of the city. We study the main structural properties of these mobility networks by focusing on their spatio-temporal patterns. Our goal is to relate these patterns with the partition into six socio-economic compartments of these two societies. Our results show that spatial and temporal patterns vary across these socio-economic groups. In particular, the two datasets show that as wealth increases the early-morning activity is delayed, the midday peak becomes smoother and the spatial distribution of trips becomes more localized.

Smart monitoring and controlling of Pandemic Influenza A (H1N1) using Social Network Analysis and cloud computing

H1N1 is an infectious virus which, when spread affects a large volume of the population. It is an airborne disease that spreads easily and has a high death rate. Development of healthcare support systems using cloud computing is emerging as an effective solution with the benefits of better quality of service, reduced costs and flexibility. In this paper, an effective cloud computing architecture is proposed which predicts H1N1 infected patients and provides preventions to control infection rate. It consists of four processing components along with secure cloud storage medical database. The random decision tree is used to initially assess the infection in any patient depending on his/her symptoms. Social Network Analysis (SNA) is used to present the state of the outbreak. The proposed architecture is tested on synthetic data generated for two million users. The system provided 94% accuracy for the classification and around 81% of the resource utilization on Amazon EC2 cloud. The key point of the paper is the use of SNA graphs to calculate role of an infected user in spreading the outbreak known as Outbreak Role Index (ORI). It will help government agencies and healthcare departments to present, analyze and prevent outbreak effectively.

Estimating the incidence reporting rates of new influenza pandemics at an early stage using travel data from the source country

During the surveillance of influenza pandemics, underreported data are a public health challenge that complicates the understanding of pandemic threats and can undermine mitigation efforts. We propose a method to estimate incidence reporting rates at early stages of new influenza pandemics using 2009 pandemic H1N1 as an example. Routine surveillance data and statistics of travellers arriving from Mexico were used. Our method incorporates changes in reporting rates such as linearly increasing trends due to the enhanced surveillance. From our results, the reporting rate was estimated at 0·46% during early stages of the pandemic in Mexico. We estimated cumulative incidence in the Mexican population to be 0·7% compared to 0·003% reported by officials in Mexico at the end of April. This method could be useful in estimation of actual cases during new influenza pandemics for policy makers to better determine appropriate control measures.

Predictive validation of an influenza spread model

BACKGROUND

Modeling plays a critical role in mitigating impacts of seasonal influenza epidemics. Complex simulation models are currently at the forefront of evaluating optimal mitigation strategies at multiple scales and levels of organization. Given their evaluative role, these models remain limited in their ability to predict and forecast future epidemics leading some researchers and public-health practitioners to question their usefulness. The objective of this study is to evaluate the predictive ability of an existing complex simulation model of influenza spread.

METHODS AND FINDINGS

We used extensive data on past epidemics to demonstrate the process of predictive validation. This involved generalizing an individual-based model for influenza spread and fitting it to laboratory-confirmed influenza infection data from a single observed epidemic (1998-1999). Next, we used the fitted model and modified two of its parameters based on data on real-world perturbations (vaccination coverage by age group and strain type). Simulating epidemics under these changes allowed us to estimate the deviation/error between the expected epidemic curve under perturbation and observed epidemics taking place from 1999 to 2006. Our model was able to forecast absolute intensity and epidemic peak week several weeks earlier with reasonable reliability and depended on the method of forecasting-static or dynamic.

CONCLUSIONS

Good predictive ability of influenza epidemics is critical for implementing mitigation strategies in an effective and timely manner. Through the process of predictive validation applied to a current complex simulation model of influenza spread, we provided users of the model (e.g. public-health officials and policy-makers) with quantitative metrics and practical recommendations on mitigating impacts of seasonal influenza epidemics. This methodology may be applied to other models of communicable infectious diseases to test and potentially improve their predictive ability.

Real-time numerical forecast of global epidemic spreading: case study of 2009 A/H1N1pdm

BACKGROUND

Mathematical and computational models for infectious diseases are increasingly used to support public-health decisions; however, their reliability is currently under debate. Real-time forecasts of epidemic spread using data-driven models have been hindered by the technical challenges posed by parameter estimation and validation. Data gathered for the 2009 H1N1 influenza crisis represent an unprecedented opportunity to validate real-time model predictions and define the main success criteria for different approaches.

METHODS

We used the Global Epidemic and Mobility Model to generate stochastic simulations of epidemic spread worldwide, yielding (among other measures) the incidence and seeding events at a daily resolution for 3,362 subpopulations in 220 countries. Using a Monte Carlo Maximum Likelihood analysis, the model provided an estimate of the seasonal transmission potential during the early phase of the H1N1 pandemic and generated ensemble forecasts for the activity peaks in the northern hemisphere in the fall/winter wave. These results were validated against the real-life surveillance data collected in 48 countries, and their robustness assessed by focusing on 1) the peak timing of the pandemic; 2) the level of spatial resolution allowed by the model; and 3) the clinical attack rate and the effectiveness of the vaccine. In addition, we studied the effect of data incompleteness on the prediction reliability.

RESULTS

Real-time predictions of the peak timing are found to be in good agreement with the empirical data, showing strong robustness to data that may not be accessible in real time (such as pre-exposure immunity and adherence to vaccination campaigns), but that affect the predictions for the attack rates. The timing and spatial unfolding of the pandemic are critically sensitive to the level of mobility data integrated into the model.

CONCLUSIONS

Our results show that large-scale models can be used to provide valuable real-time forecasts of influenza spreading, but they require high-performance computing. The quality of the forecast depends on the level of data integration, thus stressing the need for high-quality data in population-based models, and of progressive updates of validated available empirical knowledge to inform these models.

Seroprevalence of antibodies to influenza A/H1N1/2009 among transmission risk groups after the second wave in Mexico, by a virus-free ELISA method

OBJECTIVE

No serological studies have been performed in Mexico to assess the seroprevalence of influenza A/H1N1/2009 in groups of people according to the potential risk of transmission. The aim of this study was to determine the seroprevalence of antibodies against influenza A/H1N1/2009 in subjects in Mexico grouped by risk of transmission.

METHODS

Two thousand two hundred and twenty-two subjects were categorized into one of five occupation groups according to the potential risk of transmission: (1) students, (2) teachers, (3) healthcare workers, (4) institutional home residents aged >60 years, and (5) general population. Seroprevalence by potential transmission group and by age grouped into decades was determined by a virus-free ELISA method based on the recombinant receptor-binding domain of the hemagglutinin of influenza A/H1N1/2009 virus as antigen (85% sensitivity; 95% specificity). The Wilson score, Chi-square test, and logistic regression models were used for the statistical analyses.

RESULTS

Seroprevalence for students was 47.3%, for teachers was 33.9%, for older adults was 36.5%, and for the general population was 33.0%, however it was only 24.6% for healthcare workers (p=0.011). Of the students, 56.6% of those at middle school, 56.4% of those at high school, 52.7% of those at elementary school, and 31.1% of college students showed positive antibodies (p<0.001). Seroprevalence was 44.6% for college teachers, 31.6% for middle school teachers, and 29.8% for elementary school teachers, but was only 20.3% for high school teachers (p=0.002).

CONCLUSIONS

The student group was the group most affected by influenza A/H1N1/2009, while the healthcare worker group showed the lowest prevalence. Students represent a key target for preventive measures.

Early outbreak of 2009 influenza A (H1N1) in Mexico prior to identification of pH1N1 virus

BACKGROUND

In the aftermath of the global spread of 2009 influenza A (pH1N1) virus, still very little is known of the early stages of the outbreak in Mexico during the early months of the year, before the virus was identified.

METHODOLOGY/MAIN FINDINGS

We fit a simple mathematical model, the Richards model, to the number of excess laboratory-confirmed influenza cases in Mexico and Mexico City during the first 15 weeks in 2009 over the average influenza case number of the previous five baseline years of 2004-2008 during the same period to ascertain the turning point (or the peak incidence) of a wave of early influenza infections, and to estimate the transmissibility of the virus during these early months in terms of its basic reproduction number. The results indicate that there may have been an early epidemic in Mexico City as well as in all of Mexico during February/March. Based on excess influenza cases, the estimated basic reproduction number R₀ for the early outbreak was 1.59 (0.55 to 2.62) for Mexico City during weeks 5-9, and 1.25 (0.76, 1.74) for all of Mexico during weeks 5-14.

CONCLUSIONS

We established the existence of an early epidemic in Mexico City and in all of Mexico during February/March utilizing the routine influenza surveillance data, although the location of seeding is unknown. Moreover, estimates of R₀ as well as the time of peak incidence (the turning point) for Mexico City and all of Mexico indicate that the early epidemic in Mexico City in February/March had been more transmissible (larger R₀) and peaked earlier than the rest of the country. Our conclusion lends support to the possibility that the virus could have already spread to other continents prior to the identification of the virus and the reporting of lab-confirmed pH1N1 cases in North America in April.

Comparability of different methods for estimating influenza infection rates over a single epidemic wave

Estimation of influenza infection rates is important for determination of the extent of epidemic spread and for calculation of severity indicators. The authors compared estimated infection rates from paired and cross-sectional serologic surveys, rates of influenza like illness (ILI) obtained from sentinel general practitioners (GPs), and ILI samples that tested positive for influenza using data from similar periods collected during the 2009 H1N1 epidemic in Singapore. The authors performed sensitivity analyses to assess the robustness of estimates to input parameter uncertainties, and they determined sample sizes required for differing levels of precision. Estimates from paired seroconversion were 17% (95% Bayesian credible interval (BCI): 14, 20), higher than those from cross-sectional serology (12%, 95% BCI: 9, 17). Adjusted ILI estimates were 15% (95% BCI: 10, 25), and estimates computed from ILI and laboratory data were 12% (95% BCI: 8, 18). Serologic estimates were least sensitive to the risk of input parameter misspecification. ILI-based estimates were more sensitive to parameter misspecification, though this was lessened by incorporation of laboratory data. Obtaining a 5-percentage-point spread for the 95% confidence interval in infection rates would require more than 1,000 participants per serologic study, a sentinel network of 90 GPs, or 50 GPs when combined with laboratory samples. The various types of estimates will provide comparable findings if accurate input parameters can be obtained.

Estimates of the true number of cases of pandemic (H1N1) 2009, Beijing, China

During 2009, a total of 10,844 laboratory-confirmed cases of pandemic (H1N1) 2009 were reported in Beijing, People’s Republic of China. However, because most cases were not confirmed through laboratory testing, the true number is unknown. Using a multiplier model, we estimated that ≈1.46–2.30 million pandemic (H1N1) 2009 infections occurred.

Modelling and analysis of influenza A (H1N1) on networks

Background

In April 2009, a new strain of H1N1 influenza virus, referred to as pandemic influenza A (H1N1) was first detected in humans in the United States, followed by an outbreak in the state of Veracruz, Mexico. Soon afterwards, this new virus kept spreading worldwide resulting in a global outbreak. In China, the second Circular of the Ministry of Health pointed out that as of December 31, 2009, the country’s 31 provinces had reported 120,000 confirmed cases of H1N1.

Methods

We formulate an epidemic model of influenza A based on networks. We calculate the basic reproduction number and study the effects of various immunization schemes. The final size relation is derived for the network epidemic model. The model parameters are estimated via least-squares fitting of the model solution to the observed data in China.

Results

For the network model, we prove that the disease-free equilibrium is globally asymptotically stable when the basic reproduction is less than one. The final size will depend on the vaccination starting time, T, the number of infective cases at time T and immunization schemes to follow. Our theoretical results are confirmed by numerical simulations. Using the parameter estimates based on the observation data of the cumulative number of hospital notifications, we estimate the basic reproduction number R₀ to be 1.6809 in China.

Conclusions

Network modelling supplies a useful tool for studying the transmission of H1N1 in China, capturing the main features of the spread of H1N1. While a uniform, mass-immunization strategy helps control the prevalence, a targeted immunization strategy focusing on specific groups with given connectivity may better control the endemic.

Networks and Models with Heterogeneous Population Structure in Epidemiology

Heterogeneous population structure can have a profound effect on infectious disease dynamics, and is particularly important when investigating “tactical” disease control questions. At times, the nature of the network involved in the transmission of the pathogen (bacteria, virus, macro-parasite, etc.) appears to be clear; however, the nature of the network involved is dependent on the scale (e.g. within-host, between-host, or between-population), the nature of the contact, which ranges from the highly specific (e.g. sexual acts or needle sharing at the person-to-person level) to almost completely non-specific (e.g. aerosol transmission, often over long distances as can occur with the highly infectious livestock pathogen foot-and-mouth disease virus—FMDv—at the farm-to-farm level, e.g. Schley et al. in J. R. Soc. Interface 6:455–462, 2008), and the timescale of interest (e.g. at the scale of the individual, the typical infectious period of the host). Theoretical approaches to examining the implications of particular network structures on disease transmission have provided critical insight; however, a greater challenge is the integration of network approaches with data on real population structures. In this chapter, some concepts in disease modelling will be introduced, the relevance of selected network phenomena discussed, and then results from real data and their relationship to network analyses summarised. These include examinations of the patterns of air traffic and its relation to the spread of SARS in 2003 (Colizza et al. in BMC Med., 2007; Hufnagel et al. in Proc. Natl. Acad. Sci. USA 101:15124–15129, 2004), the use of the extensively documented Great Britain livestock movements network (Green et al. in J. Theor. Biol. 239:289–297, 2008; Robinson et al. in J. R. Soc. Interface 4:669–674, 2007; Vernon and Keeling in Proc. R. Soc. Lond. B, Biol. Sci. 276:469–476, 2009) and the growing interest in combining contact structure data with phylogenetics to identify real contact patterns as they directly relate to diseases of interest (Cottam et al. in PLoS Pathogens 4:1000050, 2007; Hughes et al. in PLoS Pathogens 5:1000590, 2009).