A Bayesian approach to analyse genetic variation within RNA viral populations

PIQMEE: Bayesian Phylodynamic Method for Analysis of Large Data Sets with Duplicate Sequences

Next-generation sequencing of pathogen quasispecies within a host yields data sets of tens to hundreds of unique sequences. However, the full data set often contains thousands of sequences, because many of those unique sequences have multiple identical copies. Data sets of this size represent a computational challenge for currently available Bayesian phylogenetic and phylodynamic methods. Through simulations, we explore how large data sets with duplicate sequences affect the speed and accuracy of phylogenetic and phylodynamic analysis within BEAST 2. We show that using unique sequences only leads to biases, and using a random subset of sequences yields imprecise parameter estimates. To overcome these shortcomings, we introduce PIQMEE, a BEAST 2 add-on that produces reliable parameter estimates from full data sets with increased computational efficiency as compared with the currently available methods within BEAST 2. The principle behind PIQMEE is to resolve the tree structure of the unique sequences only, while simultaneously estimating the branching times of the duplicate sequences. Distinguishing between unique and duplicate sequences allows our method to perform well even for very large data sets. Although the classic method converges poorly for data sets of 6,000 sequences when allowed to run for 7 days, our method converges in slightly more than 1 day. In fact, PIQMEE can handle data sets of around 21,000 sequences with 20 unique sequences in 14 days. Finally, we apply the method to a real, within-host HIV sequencing data set with several thousand sequences per patient.

Nanopore sequencing for rapid diagnostics of salmonid RNA viruses

Analysis of pathogen genome variation is essential for informing disease management and control measures in farmed animals. For farmed fish, the standard approach is to use PCR and Sanger sequencing to study partial regions of pathogen genomes, with second and third-generation sequencing tools yet to be widely applied. Here we demonstrate rapid and accurate sequencing of two disease-causing viruses affecting global salmonid aquaculture, salmonid alphavirus (SAV) and infectious salmon anaemia virus (ISAV), using third-generation nanopore sequencing on the MinION platform (Oxford Nanopore Technologies). Our approach complements PCR from infected material with MinION sequencing to recover genomic information that matches near perfectly to Sanger-verified references. We use this method to present the first SAV subtype-6 genome, which branches as the sister to all other SAV lineages in a genome-wide phylogenetic reconstruction. MinION sequencing offers an effective strategy for fast, genome-wide analysis of fish viruses, with major potential applications for diagnostics and robust investigations into the origins and spread of disease outbreaks.

Seven challenges in modeling pathogen dynamics within-host and across scales

The population dynamics of infectious disease is a mature field in terms of theory and to some extent, application. However for microparasites, the theory and application of models of the dynamics within a single infected host is still an open field. Further, connecting across the scales--from cellular to host level, to population level--has potential to vastly improve our understanding of pathogen dynamics and evolution. Here, we highlight seven challenges in the following areas: transmission bottlenecks, heterogeneity within host, dynamic fitness landscapes within hosts, making use of next-generation sequencing data, capturing superinfection and when and how to model more than two scales.

Genetic imprint of vaccination on simian/human immunodeficiency virus type 1 transmitted viral genomes in rhesus macaques

Understanding the genetic, antigenic and structural changes that occur during HIV-1 infection in response to pre-existing immunity will facilitate current efforts to develop an HIV-1 vaccine. Much is known about HIV-1 variation at the population level but little with regard to specific changes occurring in the envelope glycoprotein within a host in response to immune pressure elicited by antibodies. The aim of this study was to track and map specific early genetic changes occurring in the viral envelope gene following vaccination using a highly controlled viral challenge setting in the SHIV macaque model. We generated 449 full-length env sequences from vaccinees, and 63 from the virus inoculum. Analysis revealed a different pattern in the distribution and frequency of mutations in the regions of the envelope gene targeted by the vaccine as well as different patterns of diversification between animals in the naïve control group and vaccinees. Given the high stringency of the model it is remarkable that we were able to identify genetic changes associated with the vaccination. This work provides insight into the characterization of breakthrough viral populations in less than fully efficacious vaccines and illustrates the value of HIV-1 Env SHIV challenge model in macaques to unravel the mechanisms driving HIV-1 envelope genetic diversity in the presence of vaccine induced-responses.

Inferring the inter-host transmission of influenza A virus using patterns of intra-host genetic variation

Influenza A viruses (IAVs) cause acute, highly transmissible infections in a wide range of animal species. Understanding how these viruses are transmitted within and between susceptible host populations is critical to the development of effective control strategies. While viral gene sequences have been used to make inferences about IAV transmission dynamics at the epidemiological scale, their utility in accurately determining patterns of inter-host transmission in the short-term--i.e. who infected whom--has not been strongly established. Herein, we use intra-host sequence data from the viral HA1 (hemagglutinin) gene domain from two transmission studies employing different IAV subtypes in their natural hosts--H3N8 in horses and H1N1 in pigs-to determine how well these data recapitulate the known pattern of inter-host transmission. Although no mutations were fixed over the course of either experimental transmission chain, we show that some minor, transient alleles can provide evidence of host-to-host transmission and, importantly, can be distinguished from those that cannot.

Evolution of an Eurasian avian-like influenza virus in naïve and vaccinated pigs

Influenza viruses are characterized by an ability to cross species boundaries and evade host immunity, sometimes with devastating consequences. The 2009 pandemic of H1N1 influenza A virus highlights the importance of pigs in influenza emergence, particularly as intermediate hosts by which avian viruses adapt to mammals before emerging in humans. Although segment reassortment has commonly been associated with influenza emergence, an expanded host-range is also likely to be associated with the accumulation of specific beneficial point mutations. To better understand the mechanisms that shape the genetic diversity of avian-like viruses in pigs, we studied the evolutionary dynamics of an Eurasian Avian-like swine influenza virus (EA-SIV) in naïve and vaccinated pigs linked by natural transmission. We analyzed multiple clones of the hemagglutinin 1 (HA1) gene derived from consecutive daily viral populations. Strikingly, we observed both transient and fixed changes in the consensus sequence along the transmission chain. Hence, the mutational spectrum of intra-host EA-SIV populations is highly dynamic and allele fixation can occur with extreme rapidity. In addition, mutations that could potentially alter host-range and antigenicity were transmitted between animals and mixed infections were commonplace, even in vaccinated pigs. Finally, we repeatedly detected distinct stop codons in virus samples from co-housed pigs, suggesting that they persisted within hosts and were transmitted among them. This implies that mutations that reduce viral fitness in one host, but which could lead to fitness benefits in a novel host, can circulate at low frequencies.

The latest human influenza pandemic highlights the ability of influenza viruses to jump species barriers and emerge in new hosts, as well as the role of pigs in generating viruses with pandemic potential. The mutational power of influenza virus, caused by intrinsically error-prone viral polymerases, has been directly linked to viral emergence, as adaptive mutations present in the reservoir host are likely to be key to the evolution of sustained transmission in new hosts. Hence, studying how mutations are generated, maintained and transmitted in and among pigs is critical to understanding how novel viruses could emerge. Here we characterized the evolution and mutational spectra of influenza virus populations within naïve and vaccinated pigs linked by natural transmission, by analyzing multiple viral sequences obtained at different times post-infection. We show that the genetic make-up of influenza viruses in pigs is highly dynamic: the frequency of particular mutations, including those that could potentially alter host specificity or result in vaccine escape, fluctuated markedly, including one rapid fixation event. We also show that co-infections are common and multiple viruses – even defective ones – were transmitted between pigs despite being vaccinated. Our results provide empirical evidence of the complex dynamics of influenza viral populations in pigs and provide insight on the evolutionary basis of RNA viral emergence.

Binding of more than one Tva800 molecule is required for ASLV-A entry

BACKGROUND

Understanding the mechanism by which viruses enter their target cell is an essential part of understanding their infectious cycle. Previous studies have focussed on the multiplicity of viral envelope proteins that need to bind to their cognate receptor to initiate entry. Avian sarcoma and leukosis virus Envelope protein (ASLV Env) mediates entry via a receptor, Tva, which can be attached to the cell surface either by a phospholipid anchor (Tva800) or a transmembrane domain (Tva950). In these studies, we have now investigated the number of target receptors necessary for entry of ASLV Env-pseudotyped virions.

RESULTS

Using titration and modelling experiments we provide evidence that binding of more than one receptor, probably two, is needed for entry of virions via Tva800. However, binding of just one Tva950 receptor is sufficient for successful entry.

CONCLUSIONS

The different modes of attachment of Tva800 and Tva950 to the cell membrane have important implications for the utilisation of these proteins as receptors for viral binding and/or uptake.

Molecular evolution analysis of the human immunodeficiency virus type 1 envelope in simian/human immunodeficiency virus-infected macaques: implications for challenge dose selection

Since the demonstration that almost 80% of human immunodeficiency virus type 1 (HIV-1) infections result from the transmission of a single variant from the donor, biological features similar to those of HIV mucosal transmission have been reported for macaques inoculated with simian immunodeficiency virus (SIV). Here we describe the early diversification events and the impact of challenge doses on viral kinetics and on the number of variants transmitted in macaques infected with the chimeric simian/human immunodeficiency virus SHIV(sf162p4). We show that there is a correlation between the dose administered and the number of variants transmitted and that certain inoculum variants are preferentially transmitted. This could provide insight into the viral determinants of transmission and could aid in vaccine development. Challenge through the mucosal route with high doses results in the transmission of multiple variants in all the animals. Such an unrealistic scenario could underestimate potential intervention measures. We thus propose the use of molecular evolution analysis to aid in the determination of challenge doses that better mimic the transmission dynamics seen in natural HIV-1 infection.