Location-specific patterns of exposure to recent pre-pandemic strains of influenza A in southern China

Reconstructed influenza A/H3N2 infection histories reveal variation in incidence and antibody dynamics over the life course

Humans experience many influenza infections over their lives, resulting in complex and varied immunological histories. Although experimental and quantitative analyses have improved our understanding of the immunological processes defining an individual's antibody repertoire, how these within-host processes are linked to population-level influenza epidemiology remains unclear. Here, we used a multi-level mathematical model to jointly infer antibody dynamics and individual-level lifetime influenza A/H3N2 infection histories for 1,130 individuals in Guangzhou, China, using 67,683 haemagglutination inhibition (HI) assay measurements against 20 A/H3N2 strains from repeat serum samples collected between 2009 and 2015. These estimated infection histories allowed us to reconstruct historical seasonal influenza patterns and to investigate how influenza incidence varies over time, space and age in this population. We estimated median annual influenza infection rates to be approximately 18% from 1968 to 2015, but with substantial variation between years. 88% of individuals were estimated to have been infected at least once during the study period (2009-2015), and 20% were estimated to have three or more infections in that time. We inferred decreasing infection rates with increasing age, and found that annual attack rates were highly correlated across all locations, regardless of their distance, suggesting that age has a stronger impact than fine-scale spatial effects in determining an individual's antibody profile. Finally, we reconstructed each individual's expected antibody profile over their lifetime and inferred an age-stratified relationship between probability of infection and HI titre. Our analyses show how multi-strain serological panels provide rich information on long term, epidemiological trends, within-host processes and immunity when analyzed using appropriate inference methods, and adds to our understanding of the life course epidemiology of influenza A/H3N2.

Impact of Shifting University Policies During the COVID-19 Pandemic on Self-Reported Employee Social Networks

Objectives

To ascertain if faculty and staff were the link between the two COVID-19 outbreaks in a rural university county, and if the local university's COVID-19 policies affected contact rates of their employees across all its campuses.

Methods

We conducted two anonymous, voluntary online surveys for faculty and staff of a PAC-12 university on their contact patterns both within and outside the university during the COVID-19 pandemic. One was asked when classes were virtual, and another when classes were in-person but masking. Participants were asked about the individuals they encountered, the type and location of the interactions, what COVID-19 precautions were taken - if any, as well as general questions about their location and COVID-19.

Results

We received 271 responses from the first survey and 124 responses from the second. The first survey had a median of 3 contacts/respondent, with the second having 7 contacts/respondent (p<0.001). During the first survey, most contacts were family contacts (Spouse, Children), with the second survey period having Strangers and Students having the most contact (p<0.001). Over 50% of the first survey contacts happened at their home, while the second survey had 40% at work and 35% at home. Both respondents and contacts masked 42% and 46% of the time for the two surveys respectively (p<0.01).

Conclusion

For future pandemics, it would be wise to take employees into account when trying to plan for the safety of university students, employees, and surrounding communities. The main places to be aware of and potentially push infectious disease precautions would be on campus, especially confined spaces like offices or small classrooms, and the home, as these tend to be the largest areas of non-masked close contact.

Long term intrinsic cycling in human life course antibody responses to influenza A(H3N2): an observational and modeling study

Background

Over a life course, human adaptive immunity to antigenically mutable pathogens exhibits competitive and facilitative interactions. We hypothesize that such interactions may lead to cyclic dynamics in immune responses over a lifetime.

Methods

To investigate the cyclic behavior, we analyzed hemagglutination inhibition titers against 21 historical influenza A(H3N2) strains spanning 47 years from a cohort in Guangzhou, China, and applied Fourier spectrum analysis. To investigate possible biological mechanisms, we simulated individual antibody profiles encompassing known feedbacks and interactions due to generally recognized immunological mechanisms.

Results

We demonstrated a long-term periodicity (about 24 years) in individual antibody responses. The reported cycles were robust to analytic and sampling approaches. Simulations suggested that individual-level cross-reaction between antigenically similar strains likely explains the reported cycle. We showed that the reported cycles are predictable at both individual and birth cohort level and that cohorts show a diversity of phases of these cycles. Phase of cycle was associated with the risk of seroconversion to circulating strains, after accounting for age and pre-existing titers of the circulating strains.

Conclusions

Our findings reveal the existence of long-term periodicities in individual antibody responses to A(H3N2). We hypothesize that these cycles are driven by preexisting antibody responses blunting responses to antigenically similar pathogens (by preventing infection and/or robust antibody responses upon infection), leading to reductions in antigen-specific responses over time until individual's increasing risk leads to an infection with an antigenically distant enough virus to generate a robust immune response. These findings could help disentangle cohort effects from individual-level exposure histories, improve our understanding of observed heterogeneous antibody responses to immunizations, and inform targeted vaccine strategy.

Funding

This study was supported by grants from the NIH R56AG048075 (DATC, JL), NIH R01AI114703 (DATC, BY), the Wellcome Trust 200861/Z/16/Z (SR), and 200187/Z/15/Z (SR). This work was also supported by research grants from Guangdong Government HZQB-KCZYZ-2021014 and 2019B121205009 (YG and HZ). DATC, JMR and SR acknowledge support from the National Institutes of Health Fogarty Institute (R01TW0008246). JMR acknowledges support from the Medical Research Council (MR/S004793/1) and the Engineering and Physical Sciences Research Council (EP/N014499/1). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

The interplay between disease spreading and awareness diffusion in multiplex networks with activity-driven structure

The interplay between disease and awareness has been extensively studied in static networks. However, most networks in reality will evolve over time. Based on this, we propose a novel epidemiological model in multiplex networks. In this model, the disease spreading layer is a time-varying network generated by the activity-driven model, while the awareness diffusion layer is a static network, and the heterogeneity of individual infection and recovery ability is considered. First, we extend the microscopic Markov chain approach to analytically obtain the epidemic threshold of the model. Then, we simulate the spread of disease and find that stronger heterogeneity in the individual activities of a physical layer can promote disease spreading, while stronger heterogeneity of the virtual layer network will hinder the spread of disease. Interestingly, we find that when the individual infection ability follows Gaussian distribution, the heterogeneity of infection ability has little effect on the spread of disease, but it will significantly affect the epidemic threshold when the individual infection ability follows power-law distribution. Finally, we find the emergence of a metacritical point where the diffusion of awareness is able to control the onset of the epidemics. Our research could cast some light on exploring the dynamics of epidemic spreading in time-varying multiplex networks.

Influenza and other respiratory viral infections associated with absence from school among schoolchildren in Pittsburgh, Pennsylvania, USA: a cohort study

Background

Information on the etiology and age-specific burden of respiratory viral infections among school-aged children remains limited. Though school aged children are often recognized as driving the transmission of influenza as well as other respiratory viruses, little detailed information is available on the distribution of respiratory infections among children of different ages within this group. Factors other than age including gender and time spent in school may also be important in determining risk of infection but have been little studied in this age group.

Methods

We conducted a cohort study to determine the etiology of influenza like illness (ILI) among 2519 K–12 students during the 2012–13 influenza season. We obtained nasal swabs from students with ILI-related absences. Generalized linear mixed-effect regressions determined associations of outcomes, including ILI and laboratory-confirmed respiratory virus infection, with school grade and other covariates.

Results

Overall, 459 swabs were obtained from 552 ILI–related absences. Respiratory viruses were found in 292 (63.6%) samples. Influenza was found in 189 (41.2%) samples. With influenza B found in 134 (70.9%). Rates of influenza B were significantly higher in grades 1 (10.1, 95% CI 6.8–14.4%), 2 (9.7, 6.6–13.6%), 3 (9.3, 6.3–13.2%), and 4 (9.9, 6.8–13.8%) than in kindergarteners (3.2, 1.5–6.0%). After accounting for grade, sex and self-reported vaccination status, influenza B infection risk was lower among kindergarteners in half-day programs compared to kindergarteners in full-day programs (OR = 0.19; 95% CI 0.08–0.45).

Conclusions

ILI and influenza infection is concentrated in younger schoolchildren. Reduced infection by respiratory viruses is associated with a truncated school day for kindergarteners but this finding requires further investigation in other grades and populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-05922-1.

Using serological measures to estimate influenza incidence in the presence of secular trends in exposure and immuno-modulation of antibody response

BACKGROUND

Influenza infection is often measured by a fourfold antibody titer increase over an influenza season (ie seroconversion). However, this approach may fail when influenza seasons are less distinct as it does not account for transient effects from recent infections. Here, we present a method to determine seroconversion for non-paired sera, adjusting for changes in individuals' antibody titers to influenza due to the transient impact of recent exposures, varied sampling times, and laboratory processes.

METHODS

We applied our method using data for five H3N2 strains collected from 942 individuals, aged 2-90 years, during the first two study visits of the Fluscape cohort study (2009-2012) in Guangzhou, China.

RESULTS

After adjustment, apparent seroconversion rates for non-circulating strains decreased while we observed a 20% increase in seroconversion rates to recently circulating strains. When examining seroconversion to the most recently circulating strain (A/Brisbane/20/2007) in our study, participants aged under 18, and over 64 had the highest seroconversion rates compared to other age groups.

CONCLUSIONS

Our results highlight the need for improved methods when using antibody titers as an endpoint in settings where there is no clear influenza "off" season. Methods, like those presented here, that use titers from circulating and non-circulating strains may be key.

Life course exposures continually shape antibody profiles and risk of seroconversion to influenza

Complex exposure histories and immune mediated interactions between influenza strains contribute to the life course of human immunity to influenza. Antibody profiles can be generated by characterizing immune responses to multiple antigenically variant strains, but how these profiles vary across individuals and determine future responses is unclear. We used hemagglutination inhibition titers from 21 H3N2 strains to construct 777 paired antibody profiles from people aged 2 to 86, and developed novel metrics to capture features of these profiles. Total antibody titer per potential influenza exposure increases in early life, then decreases in middle age. Increased titers to one or more strains were seen in 97.8% of participants during a roughly four-year interval, suggesting widespread influenza exposure. While titer changes were seen to all strains, recently circulating strains exhibited the greatest titer rise. Higher pre-existing, homologous titers at baseline reduced the risk of seroconversion to recent strains. After adjusting for homologous titer, we also found an increased frequency of seroconversion against recent strains among those with higher immunity to older previously exposed strains. Including immunity to previously exposures also improved the deviance explained by the models. Our results suggest that a comprehensive quantitative description of immunity encompassing past exposures could lead to improved correlates of risk of influenza infection.

Timescales of influenza A/H3N2 antibody dynamics

Human immunity influences the evolution and impact of influenza strains. Because individuals are infected with multiple influenza strains during their lifetime, and each virus can generate a cross-reactive antibody response, it is challenging to quantify the processes that shape observed immune responses or to reliably detect recent infection from serological samples. Using a Bayesian model of antibody dynamics at multiple timescales, we explain complex cross-reactive antibody landscapes by inferring participants' histories of infection with serological data from cross-sectional and longitudinal studies of influenza A/H3N2 in southern China and Vietnam. We find that individual-level influenza antibody profiles can be explained by a short-lived, broadly cross-reactive response that decays within a year to leave a smaller long-term response acting against a narrower range of strains. We also demonstrate that accounting for dynamic immune responses alongside infection history can provide a more accurate alternative to traditional definitions of seroconversion for the estimation of infection attack rates. Our work provides a general model for quantifying aspects of influenza immunity acting at multiple timescales based on contemporary serological data and suggests a two-armed immune response to influenza infection consistent with competitive dynamics between B cell populations. This approach to analysing multiple timescales for antigenic responses could also be applied to other multistrain pathogens such as dengue and related flaviviruses.

Evidence for history-dependence of influenza pandemic emergence

Influenza A viruses have caused a number of global pandemics, with considerable mortality in humans. Here, we analyse the time periods between influenza pandemics since 1700 under different assumptions to determine whether the emergence of new pandemic strains is a memoryless or history-dependent process. Bayesian model selection between exponential and gamma distributions for these time periods gives support to the hypothesis of history-dependence under eight out of nine sets of modelling assumptions. Using the fitted parameters to make predictions shows a high level of variability in the modelled number of pandemics from 2010–2110. The approach we take here relies on limited data, so is uncertain, but it provides cheap, safe and direct evidence relating to pandemic emergence, a field where indirect measurements are often made at great risk and cost.

Synchronicity of influenza activity within Phoenix, AZ during the 2015-2016 seasonal epidemic

Background

Variability in the timing of influenza epidemics has been observed across global and regional scales, but this variability has not been studied extensively at finer spatial scales. As such, the aim of this study was to test whether influenza cases were synchronized across sites and/or age-groups within a major city.

Methods

We used influenza cases identified by rapid influenza tests from a network of clinics across Phoenix, AZ during the 2015–2016 influenza A season. We used a combination of KS tests and a bootstrapping approach to evaluate whether the temporal distribution of cases varied by site and/or age group.

Results

Our analysis indicates that the timing of influenza cases during the 2015–2016 seasonal influenza epidemic were generally synchronized across sites and age groups. That said, we did observe some statistically significant differences in the timing of cases across some sites, and by site and age group. We found no evidence that influenza activity consistently begins or peaks earlier in children than in adults.

Conclusions

To our knowledge, this is the first study to investigate differences in the intra-urban timing of influenza using influenza-specific case data. We were able to show evidence that influenza cases are not entirely synchronized across an urban area, but the differences we observed were relatively minor. It is important to understand the geographic scale at which influenza is synchronized in order to gain a better understanding of local transmission dynamics, and to determine the appropriate geographic scale that influenza surveillance data should be aggregated for prediction and warning systems.

A comparison of hemagglutination inhibition and neutralization assays for characterizing immunity to seasonal influenza A

BACKGROUND

Serum antibody to influenza can be used to identify past exposure and measure current immune status. The two most common methods for measuring this are the hemagglutination inhibition assay (HI) and the viral neutralization assay (NT), which have not been systematically compared for a large number of influenza viruses.

METHODS

A total of 151 study participants from near Guangzhou, China, were enrolled in 2009 and provided serum. HI and NT assays were performed for 12 historic and recently circulating strains of seasonal influenza A. We compared titers using Spearman correlation and fit models to predict NT using HI results.

RESULTS

We observed high positive mean correlation between HI and NT assays (Spearman's rank correlation, ρ=.86) across all strains. Correlation was highest within subtypes and within close proximity in time. Overall, an HI=20 corresponded to NT=10, and HI=40 corresponded to NT=20. Linear regression of log(NT) on log(HI) was statistically significant, with age modifying this relationship. Strain-specific area under a curve (AUC) indicated good accuracy (>80%) for predicting NT with HI.

CONCLUSIONS

While we found high overall correspondence of titers between NT and HI assays for seasonal influenza A, no exact equivalence between assays could be determined. This was further complicated by correspondence between titers changing with age. These findings support generalized comparison of results between assays and give further support for use of the hemagglutination inhibition assay over the more resource intensive viral neutralization assay for seasonal influenza A, although attention should be given to the effect of age on these assays.

Concurrent epidemics of influenza A/H3N2 and A/H1N1pdm in Southern China: A serial cross-sectional study

OBJECTIVES

This study aimed to elucidate the antibody response pattern of multiple influenza subtypes through a 4-year serological study of a general population in Shenzhen, Southern China.

METHODS

A serial cross-sectional serological survey was conducted at eight time points between 2009 and 2012. A total number of 5876 subjects were recruited from all age groups. The influenza subtypes tested were A/H1N1, A/H3N2, B/Yamagata, B/Victoria, and A/H1N1pdm. Genetic sequencing and phylogenetic analysis were performed on 127 H3 genes and 28 H1pdm genes.

RESULTS

We found concurrent epidemics of A/H3N2 and A/H1N1pdm with simultaneous peak times at March 2011. A/H3N2 was the dominant subtype. Ten residue substitutions (S61N, T64I, K78E, K160N, N161S, A214S, T228A, A229V, V239I, N294K, and N328S) were found in the H3 gene that might underlie its epidemic. The elderly group showed an antibody response cycle that was weaker in magnitude and slower in peak time than in younger groups.

CONCLUSIONS

The study provides a broad transmission picture and epidemiological characteristics of the major flu subtypes. The findings suggest that it may be necessary to include the A/H1N1pdm strain to the current trivalent or quadrivalent vaccine design. The delayed antibody response cycle in the elderly group indicates the need for better protection of elderly people that might be achieved by an earlier vaccination at a higher dose.

Coupled disease-behavior dynamics on complex networks: A review

It is increasingly recognized that a key component of successful infection control efforts is understanding the complex, two-way interaction between disease dynamics and human behavioral and social dynamics. Human behavior such as contact precautions and social distancing clearly influence disease prevalence, but disease prevalence can in turn alter human behavior, forming a coupled, nonlinear system. Moreover, in many cases, the spatial structure of the population cannot be ignored, such that social and behavioral processes and/or transmission of infection must be represented with complex networks. Research on studying coupled disease-behavior dynamics in complex networks in particular is growing rapidly, and frequently makes use of analysis methods and concepts from statistical physics. Here, we review some of the growing literature in this area. We contrast network-based approaches to homogeneous-mixing approaches, point out how their predictions differ, and describe the rich and often surprising behavior of disease-behavior dynamics on complex networks, and compare them to processes in statistical physics. We discuss how these models can capture the dynamics that characterize many real-world scenarios, thereby suggesting ways that policy makers can better design effective prevention strategies. We also describe the growing sources of digital data that are facilitating research in this area. Finally, we suggest pitfalls which might be faced by researchers in the field, and we suggest several ways in which the field could move forward in the coming years.

Capturing the dynamics of pathogens with many strains

Pathogens that consist of multiple antigenic variants are a serious public health concern. These infections, which include dengue virus, influenza and malaria, generate substantial morbidity and mortality. However, there are considerable theoretical challenges involved in modelling such infections. As well as describing the interaction between strains that occurs as a result cross-immunity and evolution, models must balance biological realism with mathematical and computational tractability. Here we review different modelling approaches, and suggest a number of biological problems that are potential candidates for study with these methods. We provide a comprehensive outline of the benefits and disadvantages of available frameworks, and describe what biological information is preserved and lost under different modelling assumptions. We also consider the emergence of new disease strains, and discuss how models of pathogens with multiple strains could be developed further in future. This includes extending the flexibility and biological realism of current approaches, as well as interface with data.

Statistical identifiability and sample size calculations for serial seroepidemiology

•

We investigate whether disease dynamics can be inferred by repeated serum collections.

•

Measuring antibody waning is critical for inference in serological time series.

•

Collecting 200 samples every 2 months allows for inference of transmission parameters.

•

Low-level seasonality is difficult to detect statistically.

Inference on disease dynamics is typically performed using case reporting time series of symptomatic disease. The inferred dynamics will vary depending on the reporting patterns and surveillance system for the disease in question, and the inference will miss mild or underreported epidemics. To eliminate the variation introduced by differing reporting patterns and to capture asymptomatic or subclinical infection, inferential methods can be applied to serological data sets instead of case reporting data. To reconstruct complete disease dynamics, one would need to collect a serological time series. In the statistical analysis presented here, we consider a particular kind of serological time series with repeated, periodic collections of population-representative serum. We refer to this study design as a serial seroepidemiology (SSE) design, and we base the analysis on our epidemiological knowledge of influenza. We consider a study duration of three to four years, during which a single antigenic type of influenza would be circulating, and we evaluate our ability to reconstruct disease dynamics based on serological data alone. We show that the processes of reinfection, antibody generation, and antibody waning confound each other and are not always statistically identifiable, especially when dynamics resemble a non-oscillating endemic equilibrium behavior. We introduce some constraints to partially resolve this confounding, and we show that transmission rates and basic reproduction numbers can be accurately estimated in SSE study designs. Seasonal forcing is more difficult to identify as serology-based studies only detect oscillations in antibody titers of recovered individuals, and these oscillations are typically weaker than those observed for infected individuals. To accurately estimate the magnitude and timing of seasonal forcing, serum samples should be collected every two months and 200 or more samples should be included in each collection; this sample size estimate is sensitive to the antibody waning rate and the assumed level of seasonal forcing.

Estimating the life course of influenza A(H3N2) antibody responses from cross-sectional data

The immunity of a host population against specific influenza A strains can influence a number of important biological processes, from the emergence of new virus strains to the effectiveness of vaccination programmes. However, the development of an individual's long-lived antibody response to influenza A over the course of a lifetime remains poorly understood. Accurately describing this immunological process requires a fundamental understanding of how the mechanisms of boosting and cross-reactivity respond to repeated infections. Establishing the contribution of such mechanisms to antibody titres remains challenging because the aggregate effect of immune responses over a lifetime are rarely observed directly. To uncover the aggregate effect of multiple influenza infections, we developed a mechanistic model capturing both past infections and subsequent antibody responses. We estimated parameters of the model using cross-sectional antibody titres to nine different strains spanning 40 years of circulation of influenza A(H3N2) in southern China. We found that "antigenic seniority" and quickly decaying cross-reactivity were important components of the immune response, suggesting that the order in which individuals were infected with influenza strains shaped observed neutralisation titres to a particular virus. We also obtained estimates of the frequency and age distribution of influenza infection, which indicate that although infections became less frequent as individuals progressed through childhood and young adulthood, they occurred at similar rates for individuals above age 30 y. By establishing what are likely to be important mechanisms driving epochal trends in population immunity, we also identified key directions for future studies. In particular, our results highlight the need for longitudinal samples that are tested against multiple historical strains. This could lead to a better understanding of how, over the course of a lifetime, fast, transient antibody dynamics combine with the longer-term immune responses considered here.

Epidemic spreading in metapopulation networks with heterogeneous infection rates

In this paper, we study epidemic spreading in metapopulation networks wherein each node represents a subpopulation symbolizing a city or an urban area and links connecting nodes correspond to the human traveling routes among cities. Differently from previous studies, we introduce a heterogeneous infection rate to characterize the effect of nodes' local properties, such as population density, individual health habits, and social conditions, on epidemic infectivity. By means of a mean-field approach and Monte Carlo simulations, we explore how the heterogeneity of the infection rate affects the epidemic dynamics, and find that large fluctuations of the infection rate have a profound impact on the epidemic threshold as well as the temporal behavior of the prevalence above the epidemic threshold. This work can refine our understanding of epidemic spreading in metapopulation networks with the effect of nodes' local properties.

Determinants of influenza transmission in South East Asia: insights from a household cohort study in Vietnam

To guide control policies, it is important that the determinants of influenza transmission are fully characterized. Such assessment is complex because the risk of influenza infection is multifaceted and depends both on immunity acquired naturally or via vaccination and on the individual level of exposure to influenza in the community or in the household. Here, we analyse a large household cohort study conducted in 2007–2010 in Vietnam using innovative statistical methods to ascertain in an integrative framework the relative contribution of variables that influence the transmission of seasonal (H1N1, H3N2, B) and pandemic H1N1pdm09 influenza. Influenza infection was diagnosed by haemagglutination-inhibition (HI) antibody assay of paired serum samples. We used a Bayesian data augmentation Markov chain Monte Carlo strategy based on digraphs to reconstruct unobserved chains of transmission in households and estimate transmission parameters. The probability of transmission from an infected individual to another household member was 8% (95% CI, 6%, 10%) on average, and varied with pre-season titers, age and household size. Within households of size 3, the probability of transmission from an infected member to a child with low pre-season HI antibody titers was 27% (95% CI 21%–35%). High pre-season HI titers were protective against infection, with a reduction in the hazard of infection of 59% (95% CI, 44%–71%) and 87% (95% CI, 70%–96%) for intermediate (1∶20–1∶40) and high (≥1∶80) HI titers, respectively. Even after correcting for pre-season HI titers, adults had half the infection risk of children. Twenty six percent (95% CI: 21%, 30%) of infections may be attributed to household transmission. Our results highlight the importance of integrated analysis by influenza sub-type, age and pre-season HI titers in order to infer influenza transmission risks in and outside of the household.

Influenza causes an estimated three to five million severe illnesses worldwide each year. In order to guide control policies it is important to determine the key risk factors for transmission. This is often done by studying transmission in households but in the past, analysis of such data has suffered from important simplifying assumptions (for example not being able to account for the effect of biological markers of protection like pre-season antibody titers). We have developed new statistical methods that address these limitations and applied them to a large household cohort study conducted in 2007–2010 in Vietnam. By comparing a large range of model variants, we have obtained unique insights into the patterns and determinants of transmission of seasonal (H1N1, H3N2, B) and pandemic H1N1pdm09 influenza in South East Asia. This includes estimating the proportion of cases attributed to household transmission, the average household transmission probability, the protection afforded by pre-season HI titers, and the effect of age on infection risk after correcting for pre-season HI titers.

Spatial epidemiology of networked metapopulation: an overview

An emerging disease is one infectious epidemic caused by a newly transmissible pathogen, which has either appeared for the first time or already existed in human populations, having the capacity to increase rapidly in incidence as well as geographic range. Adapting to human immune system, emerging diseases may trigger large-scale pandemic spreading, such as the transnational spreading of SARS, the global outbreak of A(H1N1), and the recent potential invasion of avian influenza A(H7N9). To study the dynamics mediating the transmission of emerging diseases, spatial epidemiology of networked metapopulation provides a valuable modeling framework, which takes spatially distributed factors into consideration. This review elaborates the latest progresses on the spatial metapopulation dynamics, discusses empirical and theoretical findings that verify the validity of networked metapopulations, and the sketches application in evaluating the effectiveness of disease intervention strategies as well.

Social mixing patterns in rural and urban areas of southern China

A dense population, global connectivity and frequent human-animal interaction give southern China an important role in the spread and emergence of infectious disease. However, patterns of person-to-person contact relevant to the spread of directly transmitted infections such as influenza remain poorly quantified in the region. We conducted a household-based survey of travel and contact patterns among urban and rural populations of Guangdong, China. We measured the character and distance from home of social encounters made by 1821 individuals. Most individuals reported 5-10 h of contact with around 10 individuals each day; however, both distributions have long tails. The distribution of distance from home at which contacts were made is similar: most were within a kilometre of the participant's home, while some occurred further than 500 km away. Compared with younger individuals, older individuals made fewer contacts which tended to be closer to home. There was strong assortativity in age-based contact rates. We found no difference between the total number or duration of contacts between urban and rural participants, but urban participants tended to make contacts closer to home. These results can improve mathematical models of infectious disease emergence, spread and control in southern China and throughout the region.

How human location-specific contact patterns impact spatial transmission between populations?

The structured-population model has been widely used to study the spatial transmission of epidemics in human society. Many seminal works have demonstrated the impact of human mobility on the epidemic threshold, assuming that the contact pattern of individuals is mixing homogeneously. Inspired by the recent evidence of location-related factors in reality, we introduce two categories of location-specific heterogeneous human contact patterns into a phenomenological model based on the commuting and contagion processes, which significantly decrease the epidemic threshold and thus favor the outbreak of diseases. In more detail, we find that a monotonic mode presents for the variance of disease prevalence in dependence on the contact rates under the destination-driven contact scenario; while under the origin-driven scenario, enhancing the contact rate counterintuitively weakens the disease prevalence in some parametric regimes. The inclusion of heterogeneity of human contacts is expected to provide valuable support to public health implications.

Estimating the reproductive number in the presence of spatial heterogeneity of transmission patterns

Background

Estimates of parameters for disease transmission in large-scale infectious disease outbreaks are often obtained to represent large groups of people, providing an average over a potentially very diverse area. For control measures to be more effective, a measure of the heterogeneity of the parameters is desirable.

Methods

We propose a novel extension of a network-based approach to estimating the reproductive number. With this we can incorporate spatial and/or demographic information through a similarity matrix. We apply this to the 2009 Influenza pandemic in South Africa to understand the spatial variability across provinces. We explore the use of five similarity matrices to illustrate their impact on the subsequent epidemic parameter estimates.

Results

When treating South Africa as a single entity with homogeneous transmission characteristics across the country, the basic reproductive number, R₀, (and imputation range) is 1.33 (1.31, 1.36). When fitting a new model for each province with no inter-province connections this estimate varies little (1.23-1.37). Using the proposed method with any of the four similarity measures yields an overall R₀ that varies little across the four new models (1.33 to 1.34). However, when allowed to vary across provinces, the estimated R₀ is greater than one consistently in only two of the nine provinces, the most densely populated provinces of Gauteng and Western Cape.

Conclusions

Our results suggest that the spatial heterogeneity of influenza transmission was compelling in South Africa during the 2009 pandemic. This variability makes a qualitative difference in our understanding of the epidemic. While the cause of this fluctuation might be partially due to reporting differences, there is substantial evidence to warrant further investigation.

The consortium for the standardization of influenza seroepidemiology (CONSISE): a global partnership to standardize influenza seroepidemiology and develop influenza investigation protocols to inform public health policy

CONSISE - The consortium for the Standardization of Influenza Seroepidemiology - is a global partnership to develop influenza investigation protocols and standardize seroepidemiology to inform health policy. This international partnership was formed in 2011 and was created out of a need, identified during the 2009 H1N1 pandemic, for timely seroepidemiological data to better estimate pandemic virus infection severity and attack rates to inform policy decisions. CONSISE has developed into a consortium of two interactive working groups: epidemiology and laboratory, with a steering committee composed of individuals from several organizations. CONSISE has had two international meetings with more planned for 2013. We seek additional members from public health agencies, academic institutions and other interested parties.

The role of social contacts and original antigenic sin in shaping the age pattern of immunity to seasonal influenza

Recent serological studies of seasonal influenza A in humans suggest a striking characteristic profile of immunity against age, which holds across different countries and against different subtypes of influenza. For both H1N1 and H3N2, the proportion of the population seropositive to recently circulated strains peaks in school-age children, reaches a minimum between ages 35–65, then rises again in the older ages. This pattern is little understood. Variable mixing between different age classes can have a profound effect on disease dynamics, and is hence the obvious candidate explanation for the profile, but using a mathematical model of multiple influenza strains, we see that age dependent transmission based on mixing data from social contact surveys cannot on its own explain the observed pattern. Instead, the number of seropositive individuals in a population may be a consequence of ‘original antigenic sin’; if the first infection of a lifetime dominates subsequent immune responses, we demonstrate that it is possible to reproduce the observed relationship between age and seroprevalence. We propose a candidate mechanism for this relationship, by which original antigenic sin, along with antigenic drift and vaccination, results in the age profile of immunity seen in empirical studies.

The way in which a population builds immunity to influenza affects outbreak size and the emergence of new strains. However, although age-specific immunity has been widely discussed for the 2009 influenza pandemic, the age profile of immunity to seasonal influenza remains little understood. In contrast to many infections, the proportion of people immune to recent strains peaks in school-age children then reaches a minimum between ages 35–65, before rising again in older age groups. Our results suggest that rather than variable mixing between different age groups being solely responsible, the pattern may be shaped by an effect known as ‘original antigenic sin’, by which the first infection of a lifetime dictates subsequent immune responses: instead of developing antibodies to every new virus that is encountered, the immune system may reuse the response to a similar virus it has already seen. The framework we describe, which extends theoretical models to allow for comparison with data, also opens the possibility of investigating the mechanisms behind patterns of immunity to other evolving pathogens.

Evidence for antigenic seniority in influenza A (H3N2) antibody responses in southern China

A key observation about the human immune response to repeated exposure to influenza A is that the first strain infecting an individual apparently produces the strongest adaptive immune response. Although antibody titers measure that response, the interpretation of titers to multiple strains – from the same sera – in terms of infection history is clouded by age effects, cross reactivity and immune waning. From July to September 2009, we collected serum samples from 151 residents of Guangdong Province, China, 7 to 81 years of age. Neutralization tests were performed against strains representing six antigenic clusters of H3N2 influenza circulating between 1968 and 2008, and three recent locally circulating strains. Patterns of neutralization titers were compared based on age at time of testing and age at time of the first isolation of each virus. Neutralization titers were highest for H3N2 strains that circulated in an individual's first decade of life (peaking at 7 years). Further, across strains and ages at testing, statistical models strongly supported a pattern of titers declining smoothly with age at the time a strain was first isolated. Those born 10 or more years after a strain emerged generally had undetectable neutralization titers to that strain (<1∶10). Among those over 60 at time of testing, titers tended to increase with age. The observed pattern in H3N2 neutralization titers can be characterized as one of antigenic seniority: repeated exposure and the immune response combine to produce antibody titers that are higher to more ‘senior’ strains encountered earlier in life.

The human immune response to an influenza infection is not the same for every infection. It has often been observed that we tend to have the highest antibody titer (and presumably our strongest immune response) against strains of influenza that we were exposed to early in life. In this study, we obtained blood samples from 151 people between 7 and 81 years of age and tested the samples for the concentration of antibodies to many different (H3N2) strains. We chose strains according to when they first circulated, starting with a strain isolated just after the 1968 pandemic and going all the way through to very recent strains. We found that a participant's age at the time a strain first circulated was very predictive of the strength of their antibody against that strain. Not just for the first strain they were likely to have seen, but also for the second, third and all subsequent strains circulating during their lifetime. This suggests to us that antibody titers to influenza A H3N2 follow a pattern of antigenic seniority, suggesting that we produce progressively fewer specific antibodies to each subsequent infection as we age.