Inferring pandemic growth rates from sequence data

Using sequence data to infer population dynamics is playing an increasing role in the analysis of outbreaks. The most common methods in use, based on coalescent inference, have been widely used but not extensively tested against simulated epidemics. Here, we use simulated data to test the ability of both parametric and non-parametric methods for inference of effective population size (coded in the popular BEAST package) to reconstruct epidemic dynamics. We consider a range of simulations centred on scenarios considered plausible for pandemic influenza, but our conclusions are generic for any exponentially growing epidemic. We highlight systematic biases in non-parametric effective population size estimation. The most prominent such bias leads to the false inference of slowing of epidemic spread in the recent past even when the real epidemic is growing exponentially. We suggest some sampling strategies that could reduce (but not eliminate) some of the biases. Parametric methods can correct for these biases if the infected population size is large. We also explore how some poor sampling strategies (e.g. that over-represent epidemiologically linked clusters of cases) could dramatically exacerbate bias in an uncontrolled manner. Finally, we present a simple diagnostic indicator, based on coalescent density and which can easily be applied to reconstructed phylogenies, that identifies time-periods for which effective population size estimates are less likely to be biased. We illustrate this with an application to the 2009 H1N1 pandemic

Assessing the role of live poultry trade in community-structured transmission of avian influenza in China

The live poultry trade is thought to play an important role in the spread and maintenance of highly pathogenic avian influenza A viruses (HP AIVs) in Asia. Despite an abundance of small-scale observational studies, the role of the poultry trade in disseminating AIV over large geographic areas is still unclear, especially for developing countries with complex poultry production systems. Here we combine virus genomes and reconstructed poultry transportation data to measure and compare the spatial spread in China of three key subtypes of AIV: H5N1, H7N9, and H5N6. Although it is difficult to disentangle the contribution of confounding factors, such as bird migration and spatial distance, we find evidence that the dissemination of these subtypes among domestic poultry is geographically continuous and likely associated with the intensity of the live poultry trade in China. Using two independent data sources and network analysis methods, we report a regional-scale community structure in China that might explain the spread of AIV subtypes in the country. The identification of this structure has the potential to inform more targeted strategies for the prevention and control of AIV in China

Antibody responses to avian influenza viruses in wild birds broaden with age

For viruses such as avian influenza, immunity within a host population can drive the emergence of new strains by selecting for viruses with novel antigens that avoid immune recognition. The accumulation of acquired immunity with age is hypothesized to affect how influenza viruses emerge and spread in species of different lifespans. Despite its importance for understanding the behaviour of avian influenza viruses, little is known about age-related accumulation of immunity in the virus's primary reservoir, wild birds. To address this, we studied the age structure of immune responses to avian influenza virus in a wild swan population (Cygnus olor), before and after the population experienced an outbreak of highly pathogenic H5N1 avian influenza in 2008. We performed haemagglutination inhibition assays on sampled sera for five avian influenza strains and show that breadth of response accumulates with age. The observed age-related distribution of antibody responses to avian influenza strains may explain the age-dependent mortality observed during the highly pathogenic H5N1 outbreak. Age structures and species lifespan are probably important determinants of viral epidemiology and virulence in birds

Short-Range Correlations and the Nuclear EMC Effect in Deuterium and Helium-3

The EMC effect in deuterium and helium-3 is studied using a convolution formalism that allows isolating the impact of high-momentum nucleons in short-ranged correlated (SRC) pairs. We assume that the modification of the structure function of bound nucleons is given by a universal (i.e. nucleus independent) function of their virtuality, and find that the effect of such modifications is dominated by nucleons in SRC pairs. This SRC-dominance of nucleon modifications is observed despite the fact that the bulk of the nuclear inelastic scattering cross-section comes from interacting with low-momentum nucleons. These findings are found to be robust to model details including nucleon modification function parametrization, free nucleon structure function and treatment of nucleon motion effects. While existing data cannot discriminate between such model details, we present predictions for measured, but not yet published, tritium EMC effect and tagged nucleon structure functions in deuterium that are sensitive to the neutron structure functions and bound nucleon modification functions.Comment: 9 pages, 9 figures, 1 table, and online supplementary material

Unifying the spatial epidemiology and molecular evolution of emerging epidemics

We introduce a conceptual bridge between the previously unlinked fields of phylogenetics and mathematical spatial ecology, which enables the spatial parameters of an emerging epidemic to be directly estimated from sampled pathogen genome sequences. By using phylogenetic history to correct for spatial autocorrelation, we illustrate how a fundamental spatial variable, the diffusion coefficient, can be estimated using robust nonparametric statistics, and how heterogeneity in dispersal can be readily quantified. We apply this framework to the spread of the West Nile virus across North America, an important recent instance of spatial invasion by an emerging infectious disease. We demonstrate that the dispersal of West Nile virus is greater and far more variable than previously measured, such that its dissemination was critically determined by rare, long-range movements that are unlikely to be discerned during field observations. Our results indicate that, by ignoring this heterogeneity, previous models of the epidemic have substantially overestimated its basic reproductive number. More generally, our approach demonstrates that easily obtainable genetic data can be used to measure the spatial dynamics of natural populations that are otherwise difficult or costly to quantify

Molecular evolution, diversity and adaptation of H7N9 viruses in China during 2013-2017

The substantial increase in prevalence and emergence of antigenically divergent or highly pathogenic influenza A(H7N9) viruses during 2016–17 raises concerns about the epizootic potential of these viruses. We investigated the evolution and adaptation of H7N9 viruses by analyzing available data and newly generated virus sequences isolated in Guangdong Province, China, during 2015–2017. Phylogenetic analyses showed that circulating H7N9 viruses belong to distinct lineages with differing spatial distributions. Hemagglutination inhibition assays performed on serum samples from patients infected with these viruses identified 3 antigenic clusters for 16 strains of different virus lineages. We used ancestral sequence reconstruction to identify parallel amino acid changes on multiple separate lineages. We inferred that mutations in hemagglutinin occur primarily at sites involved in receptor recognition or antigenicity. Our results indicate that highly pathogenic strains likely emerged from viruses circulating in eastern Guangdong Province during March 2016 and are associated with a high rate of adaptive molecular evolution

A standardized framework for accurate, high-throughput genotyping of recombinant and non-recombinant viral sequences

Human immunodeficiency virus type-1 (HIV-1), hepatitis B and C and other rapidly evolving viruses are characterized by extremely high levels of genetic diversity. To facilitate diagnosis and the development of prevention and treatment strategies that efficiently target the diversity of these viruses, and other pathogens such as human T-lymphotropic virus type-1 (HTLV-1), human herpes virus type-8 (HHV8) and human papillomavirus (HPV), we developed a rapid high-throughput-genotyping system. The method involves the alignment of a query sequence with a carefully selected set of pre-defined reference strains, followed by phylogenetic analysis of multiple overlapping segments of the alignment using a sliding window. Each segment of the query sequence is assigned the genotype and sub-genotype of the reference strain with the highest bootstrap (>70%) and bootscanning (>90%) scores. Results from all windows are combined and displayed graphically using color-coded genotypes. The new Virus-Genotyping Tools provide accurate classification of recombinant and non-recombinant viruses and are currently being assessed for their diagnostic utility. They have incorporated into several HIV drug resistance algorithms including the Stanford (http://hivdb.stanford.edu) and two European databases (http://www.umcutrecht.nl/subsite/spread-programme/ and http://www.hivrdb.org.uk/) and have been successfully used to genotype a large number of sequences in these and other databases. The tools are a PHP/JAVA web application and are freely accessible on a number of servers including

Comparative micro-epidemiology of pathogenic avian influenza virus outbreaks in a wild bird population

Understanding the epidemiological dynamics of highly pathogenic avian influenza virus (HPAIV) in wild birds is crucial for guiding effective surveillance and control measures. The spread of H5 HPAIV has been well characterized over large geographical and temporal scales. However, information about the detailed dynamics and demographics of individual outbreaks in wild birds is rare and important epidemiological parameters remain unknown. We present data from a wild population of long-lived birds (mute swans; Cygnus olor) that has experienced three outbreaks of related H5 HPAIVs in the past decade, specifically, H5N1 (2007), H5N8 (2016) and H5N6 (2017). Detailed demographic data were available and intense sampling was conducted before and after the outbreaks; hence the population is unusually suitable for exploring the natural epidemiology, evolution and ecology of HPAIV in wild birds. We show that key epidemiological features remain remarkably consistent across multiple outbreaks, including the timing of virus incursion and outbreak duration, and the presence of a strong age-structure in morbidity that likely arises from an equivalent age-structure in immunological responses. The predictability of these features across a series of outbreaks in a complex natural population is striking and contributes to our understanding of HPAIV in wild birds

First Complete Genome Sequences of Zika Virus Isolated from Febrile Patient Sera in Ecuador

Here, we present the complete genome sequences of two Zika virus (ZIKV) strains, EcEs062_16 and EcEs089_16, isolated from the sera of febrile patients in Esmeraldas City, in the northern coastal province of Esmeraldas, Ecuador, in April 2016. These are the first complete ZIKV genomes to be reported from Ecuador