Influenza emergence in the face of evolutionary constraints

Different influenza subtypes can evolve at very different rates, but the causes are not well understood. In this paper, we explore whether differences in transmissibility between subtypes can play a role if there are fitness constraints on antigenic evolution. We investigate the problem using a mathematical model that separates the interaction of strains through cross-immunity from the process of emergence for new antigenic variants. Evolutionary constraints are also included with antigenic mutation incurring a fitness cost. We show that the transmissibility of a strain can become disproportionately important in dictating the rate of antigenic drift: strains that spread only slightly more easily can have a much higher rate of emergence. Further, we see that the effect continues when vaccination is considered; a small increase in the rate of transmission can make it much harder to control the frequency at which new strains emerge. Our results not only highlight the importance of considering both transmission and fitness constraints when modelling influenza evolution, but may also help in understanding the differences between the emergence of H1N1 and H3N2 subtypes

A new method for detecting signal regions in ordered sequences of real numbers, and application to viral genomic data.

We present a fast, robust and parsimonious approach to detecting signals in an ordered sequence of numbers. Our motivation is in seeking a suitable method to take a sequence of scores corresponding to properties of positions in virus genomes, and find outlying regions of low scores. Suitable statistical methods without using complex models or making many assumptions are surprisingly lacking. We resolve this by developing a method that detects regions of low score within sequences of real numbers. The method makes no assumptions a priori about the length of such a region; it gives the explicit location of the region and scores it statistically. It does not use detailed mechanistic models so the method is fast and will be useful in a wide range of applications. We present our approach in detail, and test it on simulated sequences. We show that it is robust to a wide range of signal morphologies, and that it is able to capture multiple signals in the same sequence. Finally we apply it to viral genomic data to identify regions of evolutionary conservation within influenza and rotavirus

Epidemic interventions: insights from classic results.

Analytical expressions and approximations from simple models have performed a pivotal role in our understanding of infectious disease epidemiology. During the current COVID-19 pandemic, while there has been proliferation of increasingly complex models, still the most basic models have provided the core framework for our thinking and interpreting policy decisions. Here, classic results are presented that give insights into both the role of transmission-reducing interventions (such as social distancing) in controlling an emerging epidemic, and also what would happen if insufficient control is applied. Though these are simple results from the most basic of epidemic models, they give valuable benchmarks for comparison with the outputs of more complex modelling approaches. This article is part of the theme issue 'Modelling that shaped the early COVID-19 pandemic response in the UK'

Geographic transmission hubs of the 2009 influenza pandemic in the United States

A key issue in infectious disease epidemiology is to identify and predict geographic sites of epidemic establishment that contribute to onward spread, especially in the context of invasion waves of emerging pathogens. Conventional wisdom suggests that these sites are likely to be in densely-populated, well-connected areas. For pandemic influenza, however, epidemiological data have not been available at a fine enough geographic resolution to test this assumption. Here, we make use of fine-scale influenza-like illness incidence data derived from electronic medical claims records gathered from 834 3-digit ZIP (postal) codes across the US to identify the key geographic establishment sites, or “hubs”, of the autumn wave of the 2009 A/H1N1pdm influenza pandemic in the United States. A mechanistic spatial transmission model is fit to epidemic onset times inferred from the data. Hubs are identified by tracing the most probable transmission routes back to a likely first establishment site. Four hubs are identified: two in the southeastern US, one in the central valley of California, and one in the midwestern US. According to the model, 75% of the 834 observed ZIP-level outbreaks in the US were seeded by these four hubs or their epidemiological descendants. Counter-intuitively, the pandemic hubs do not coincide with large and well-connected cities, indicating that factors beyond population density and travel volume are necessary to explain the establishment sites of the major autumn wave of the pandemic. Geographic regions are identified where infection can be statistically traced back to a hub, providing a testable prediction of the outbreak's phylogeography. Our method therefore provides an important way forward to reconcile spatial diffusion patterns inferred from epidemiological surveillance data and pathogen sequence data. Keywords: Pandemic influenza, Transmission hubs, Metapopulation, Gravity model, Phylogeograph

Vaccine escape in a heterogeneous population: insights for SARS-CoV-2 from a simple model.

As a countermeasure to the SARS-CoV-2 pandemic, there has been swift development and clinical trial assessment of candidate vaccines, with subsequent deployment as part of mass vaccination campaigns. However, the SARS-CoV-2 virus has demonstrated the ability to mutate and develop variants, which can modify epidemiological properties and potentially also the effectiveness of vaccines. The widespread deployment of highly effective vaccines may rapidly exert selection pressure on the SARS-CoV-2 virus directed towards mutations that escape the vaccine-induced immune response. This is particularly concerning while infection is widespread. By developing and analysing a mathematical model of two population groupings with differing vulnerability and contact rates, we explore the impact of the deployment of vaccines among the population on the reproduction ratio, cases, disease abundance and vaccine escape pressure. The results from this model illustrate two insights: (i) vaccination aimed at reducing prevalence could be more effective at reducing disease than directly vaccinating the vulnerable; (ii) the highest risk for vaccine escape can occur at intermediate levels of vaccination. This work demonstrates a key principle: the careful targeting of vaccines towards particular population groups could reduce disease as much as possible while limiting the risk of vaccine escape

Five challenges in modelling interacting strain dynamics

Population epidemiological models where hosts can be infected sequentially by different strains have the potential to help us understand many important diseases. Researchers have in recent years started to develop and use such models, but the extra layer of complexity from multiple strains brings with it many technical challenges. It is therefore hard to build models which have realistic assumptions yet are tractable. Here we outline some of the main challenges in this area. First we begin with the fundamental question of how to translate from complex small-scale dynamics within a host to useful population models. Next we consider the nature of so-called "strain space". We describe two key types of host heterogeneities, and explain how models could help generate a better understanding of their effects. Finally, for diseases with many strains, we consider the challenge of modelling how immunity accumulates over multiple exposures

Spatial Transmission of 2009 Pandemic Influenza in the US

The 2009 H1N1 influenza pandemic provides a unique opportunity for detailed examination of the spatial dynamics of an emerging pathogen. In the US, the pandemic was characterized by substantial geographical heterogeneity: the 2009 spring wave was limited mainly to northeastern cities while the larger fall wave affected the whole country. Here we use finely resolved spatial and temporal influenza disease data based on electronic medical claims to explore the spread of the fall pandemic wave across 271 US cities and associated suburban areas. We document a clear spatial pattern in the timing of onset of the fall wave, starting in southeastern cities and spreading outwards over a period of three months. We use mechanistic models to tease apart the external factors associated with the timing of the fall wave arrival: differential seeding events linked to demographic factors, school opening dates, absolute humidity, prior immunity from the spring wave, spatial diffusion, and their interactions. Although the onset of the fall wave was correlated with school openings as previously reported, models including spatial spread alone resulted in better fit. The best model had a combination of the two. Absolute humidity or prior exposure during the spring wave did not improve the fit and population size only played a weak role. In conclusion, the protracted spread of pandemic influenza in fall 2009 in the US was dominated by short-distance spatial spread partially catalysed by school openings rather than long-distance transmission events. This is in contrast to the rapid hierarchical transmission patterns previously described for seasonal influenza. The findings underline the critical role that school-age children play in facilitating the geographic spread of pandemic influenza and highlight the need for further information on the movement and mixing patterns of this age group

Effectiveness of isolation, testing, contact tracing, and physical distancing on reducing transmission of SARS-CoV-2 in different settings: a mathematical modelling study.

BACKGROUND: The isolation of symptomatic cases and tracing of contacts has been used as an early COVID-19 containment measure in many countries, with additional physical distancing measures also introduced as outbreaks have grown. To maintain control of infection while also reducing disruption to populations, there is a need to understand what combination of measures-including novel digital tracing approaches and less intensive physical distancing-might be required to reduce transmission. We aimed to estimate the reduction in transmission under different control measures across settings and how many contacts would be quarantined per day in different strategies for a given level of symptomatic case incidence. METHODS: For this mathematical modelling study, we used a model of individual-level transmission stratified by setting (household, work, school, or other) based on BBC Pandemic data from 40 162 UK participants. We simulated the effect of a range of different testing, isolation, tracing, and physical distancing scenarios. Under optimistic but plausible assumptions, we estimated reduction in the effective reproduction number and the number of contacts that would be newly quarantined each day under different strategies. RESULTS: We estimated that combined isolation and tracing strategies would reduce transmission more than mass testing or self-isolation alone: mean transmission reduction of 2% for mass random testing of 5% of the population each week, 29% for self-isolation alone of symptomatic cases within the household, 35% for self-isolation alone outside the household, 37% for self-isolation plus household quarantine, 64% for self-isolation and household quarantine with the addition of manual contact tracing of all contacts, 57% with the addition of manual tracing of acquaintances only, and 47% with the addition of app-based tracing only. If limits were placed on gatherings outside of home, school, or work, then manual contact tracing of acquaintances alone could have an effect on transmission reduction similar to that of detailed contact tracing. In a scenario where 1000 new symptomatic cases that met the definition to trigger contact tracing occurred per day, we estimated that, in most contact tracing strategies, 15 000-41 000 contacts would be newly quarantined each day. INTERPRETATION: Consistent with previous modelling studies and country-specific COVID-19 responses to date, our analysis estimated that a high proportion of cases would need to self-isolate and a high proportion of their contacts to be successfully traced to ensure an effective reproduction number lower than 1 in the absence of other measures. If combined with moderate physical distancing measures, self-isolation and contact tracing would be more likely to achieve control of severe acute respiratory syndrome coronavirus 2 transmission. FUNDING: Wellcome Trust, UK Engineering and Physical Sciences Research Council, European Commission, Royal Society, Medical Research Council

A Bayesian Approach to Analyse Genetic Variation within RNA Viral Populations

The development of modern and affordable sequencing technologies has allowed the study of viral populations to an unprecedented depth. This is of particular interest for the study of within-host RNA viral populations, where variation due to error-prone polymerases can lead to immune escape, antiviral resistance and adaptation to new host species. Methods to sequence RNA virus genomes include reverse transcription (RT) and polymerase chain reaction (PCR). RT-PCR is a molecular biology technique widely used to amplify DNA from an RNA template. The method itself relies on the in vitro synthesis of copy DNA from RNA followed by multiple cycles of DNA amplification. However, this method introduces artefactual errors that can act as confounding factors when the sequence data are analysed. Although there are a growing number of published studies exploring the intra- and inter-host evolutionary dynamics of RNA viruses, the complexity of the methods used to generate sequences makes it difficult to produce probabilistic statements about the likely sources of observed sequence variants. This complexity is further compounded as both the depth of sequencing and the length of the genome segment of interest increase. Here we develop a Bayesian method to characterise and differentiate between likely structures for the background viral population. This approach can then be used to identify nucleotide sites that show evidence of change in the within-host viral population structure, either over time or relative to a reference sequence (e.g. an inoculum or another source of infection), or both, without having to build complex evolutionary models. Identification of these sites can help to inform the design of more focussed experiments using molecular biology tools, such as site-directed mutagenesis, to assess the function of specific amino acids. We illustrate the method by applying to datasets from experimental transmission of equine influenza, and a pre-clinical vaccine trial for HIV-1

SARS-CoV-2 infection in UK university students: lessons from September-December 2020 and modelling insights for future student return.

Funder: Isaac Newton Institute for Mathematical Sciences; Id: http://dx.doi.org/10.13039/501100005347Funder: Wellcome Trust; Id: http://dx.doi.org/10.13039/100004440Funder: Medical Research Council; Id: http://dx.doi.org/10.13039/501100000265Funder: UKRIFunder: University of Nottingham; Id: http://dx.doi.org/10.13039/501100000837In this paper, we present work on SARS-CoV-2 transmission in UK higher education settings using multiple approaches to assess the extent of university outbreaks, how much those outbreaks may have led to spillover in the community, and the expected effects of control measures. Firstly, we found that the distribution of outbreaks in universities in late 2020 was consistent with the expected importation of infection from arriving students. Considering outbreaks at one university, larger halls of residence posed higher risks for transmission. The dynamics of transmission from university outbreaks to wider communities is complex, and while sometimes spillover does occur, occasionally even large outbreaks do not give any detectable signal of spillover to the local population. Secondly, we explored proposed control measures for reopening and keeping open universities. We found the proposal of staggering the return of students to university residence is of limited value in terms of reducing transmission. We show that student adherence to testing and self-isolation is likely to be much more important for reducing transmission during term time. Finally, we explored strategies for testing students in the context of a more transmissible variant and found that frequent testing would be necessary to prevent a major outbreak.This work was supported by EPSRC grant no EP/R014604/1. The authors would also like to thank the Virtual Forum for Knowledge Exchange in Mathematical Sciences (V-KEMS) for the support during the workshop Unlocking higher education Spaces – What Might Mathematics Tell Us? where work on this paper was undertaken. K.J.B. acknowledges support from a University of Nottingham Anne McLaren Fellowship. E.L.F. acknowledges support via K.J.B.’s fellowship and the Nottingham BBSRC Doctoral Training Partnership. M.L.T. was supported by the UK Engineering and Physical Sciences Research Council (grant no. EP/N509620/1). E.B.-P., E.J.N., L.D., J.R.G. and M.J.T. were supported by UKRI through the JUNIPER modelling consortium (grant no. MR/V038613/1). E.M.H., L.D. and M.J.T. were supported by the Medical Research Council through the COVID-19 Rapid Response Rolling Call (grant no. MR/V009761/1). H.B.S. is funded by the Wellcome Trust and the Royal Society (grant no. 202562/Z/16/Z). J.E. is partially funded by the UK Engineering and Physical Sciences Research Council (grant no. EP/T004878/1)