Hitchhiking and the population genetic structure of avian influenza virus

Mutual information networks reveal evolutionary relationships within the influenza A virus polymerase

The influenza A virus (IAV) RNA polymerase is an essential driver of IAV evolution. Mutations that the polymerase introduces into viral genome segments during replication are the ultimate source of genetic variation, including within the three subunits of the IAV polymerase (polymerase basic protein 2, polymerase basic protein 1, and polymerase acidic protein). Evolutionary analysis of the IAV polymerase is complicated, because changes in mutation rate, replication speed, and drug resistance involve epistatic interactions among its subunits. In order to study the evolution of the human seasonal H3N2 polymerase since the 1968 pandemic, we identified pairwise evolutionary relationships among ∼7000 H3N2 polymerase sequences using mutual information (MI), which measures the information gained about the identity of one residue when a second residue is known. To account for uneven sampling of viral sequences over time, we developed a weighted MI (wMI) metric and demonstrate that wMI outperforms raw MI through simulations using a well-sampled severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) dataset. We then constructed wMI networks of the H3N2 polymerase to extend the inherently pairwise wMI statistic to encompass relationships among larger groups of residues. We included hemagglutinin (HA) in the wMI network to distinguish between functional wMI relationships within the polymerase and those potentially due to hitch-hiking on antigenic changes in HA. The wMI networks reveal coevolutionary relationships among residues with roles in replication and encapsidation. Inclusion of HA highlighted polymerase-only subgraphs containing residues with roles in the enzymatic functions of the polymerase and host adaptability. This work provides insight into the factors that drive and constrain the rapid evolution of influenza viruses.

Mutual information networks reveal evolutionary relationships within the influenza A virus polymerase

The influenza A (IAV) RNA polymerase is an essential driver of IAV evolution. Mutations that the polymerase introduces into viral genome segments during replication are the ultimate source of genetic variation, including within the three subunits of the IAV polymerase (PB2, PB1, and PA). Evolutionary analysis of the IAV polymerase is complicated, because changes in mutation rate, replication speed, and drug resistance involve epistatic interactions among its subunits. In order to study the evolution of the human seasonal H3N2 polymerase since the 1968 pandemic, we identified pairwise evolutionary relationships among ∼7000 H3N2 polymerase sequences using mutual information (MI), which measures the information gained about the identity of one residue when a second residue is known. To account for uneven sampling of viral sequences over time, we developed a weighted MI metric (wMI) and demonstrate that wMI outperforms raw MI through simulations using a well-sampled SARS-CoV-2 dataset. We then constructed wMI networks of the H3N2 polymerase to extend the inherently pairwise wMI statistic to encompass relationships among larger groups of residues. We included HA in the wMI network to distinguish between functional wMI relationships within the polymerase and those potentially due to hitchhiking on antigenic changes in HA. The wMI networks reveal coevolutionary relationships among residues with roles in replication and encapsidation. Inclusion of HA highlighted polymerase-only subgraphs containing residues with roles in the enzymatic functions of the polymerase and host adaptability. This work provides insight into the factors that drive and constrain the rapid evolution of influenza viruses.

Evolutionary features of a prolific subtype of avian influenza A virus in European waterfowl

Avian influenza A virus (AIV) is ubiquitous in waterfowl and is detected annually at high prevalence in waterfowl during the Northern Hemisphere autumn. Some AIV subtypes are globally common in waterfowl, such as H3N8, H4N6, and H6N2, and are detected in the same populations at a high frequency, annually. In order to investigate genetic features associated to the long-term maintenance of common subtypes in migratory ducks, we sequenced 248 H4 viruses isolated across 8 years (2002-9) from mallards (Anas platyrhynchos) sampled in southeast Sweden. Phylogenetic analyses showed that both H4 and N6 sequences fell into three distinct lineages, structured by year of isolation. Specifically, across the 8 years of the study, we observed lineage replacement, whereby a different HA lineage circulated in the population each year. Analysis of deduced amino acid sequences of the HA lineages illustrated key differences in regions of the globular head of hemagglutinin that overlap with established antigenic sites in homologous hemagglutinin H3, suggesting the possibility of antigenic differences among these HA lineages. Beyond HA, lineage replacement was common to all segments, such that novel genome constellations were detected across years. A dominant genome constellation would rapidly amplify in the duck population, followed by unlinking of gene segments as a result of reassortment within 2-3 weeks following introduction. These data help reveal the evolutionary dynamics exhibited by AIV on both annual and decadal scales in an important reservoir host.

Insights on recent approaches in drug discovery strategies and untapped drug targets against drug resistance

BACKGROUND

Despite the various strategies undertaken in the clinical practice, the mortality rate due to antibiotic-resistant microbes has been markedly increasing worldwide. In addition to multidrug-resistant (MDR) microbes, the "ESKAPE" bacteria are also emerging. Of course, the infection caused by ESKAPE cannot be treated even with lethal doses of antibiotics. Now, the drug resistance is also more prevalent in antiviral, anticancer, antimalarial and antifungal chemotherapies.

MAIN BODY

To date, in the literature, the quantum of research reported on the discovery strategies for new antibiotics is remarkable but the milestone is still far away. Considering the need of the updated strategies and drug discovery approaches in the area of drug resistance among researchers, in this communication, we consolidated the insights pertaining to new drug development against drug-resistant microbes. It includes drug discovery void, gene paradox, transposon mutagenesis, vitamin biosynthesis inhibition, use of non-conventional media, host model, target through quorum sensing, genomic-chemical network, synthetic viability to targets, chemical versus biological space, combinational approach, photosensitization, antimicrobial peptides and transcriptome profiling. Furthermore, we optimally briefed about antievolution drugs, nanotheranostics and antimicrobial adjuvants and then followed by twelve selected new feasible drug targets for new drug design against drug resistance. Finally, we have also tabulated the chemical structures of potent molecules against antimicrobial resistance.

CONCLUSION

It is highly recommended to execute the anti-drug resistance research as integrated approach where both molecular and genetic research needs to be as integrative objective of drug discovery. This is time to accelerate new drug discovery research with advanced genetic approaches instead of conventional blind screening.

The Ecology and Evolution of Influenza Viruses

The patterns and processes of influenza virus evolution are of fundamental importance, underpinning such traits as the propensity to emerge in new host species and the ability to rapidly generate antigenic variation. Herein, we review key aspects of the ecology and evolution of influenza viruses. We begin with an exploration of the origins of influenza viruses within the orthomyxoviruses, showing how our perception of the evolutionary history of these viruses has been transformed with metagenomic sequencing. We then outline the diversity of virus subtypes in different species and the processes by which these viruses have emerged in new hosts, with a particular focus on the role played by segment reassortment. We then turn our attention to documenting the spread and phylodynamics of seasonal influenza A and B viruses in human populations, including the drivers of antigenic evolution, and finish with a discussion of virus diversity and evolution at the scale of individual hosts.

H7N6 low pathogenic avian influenza outbreak in commercial turkey farms in Chile caused by a native South American Lineage

In December 2016, low pathogenic avian influenza (LPAI) caused by an H7N6 subtype was confirmed in a grow-out turkey farm located in Valparaiso Region, Chile. Depopulation of exposed animals, zoning, animal movement control and active surveillance were implemented to contain the outbreak. Two weeks later, a second grow-out turkey farm located 70 km north of the first site was also infected by H7N6 LPAI, which subsequently spilled over to one backyard poultry flock. The virus involved in the outbreak shared a close genetic relationship with Chilean aquatic birds' viruses collected in previous years. The A/turkey/Chile/2017(H7N6) LPAI virus belonged to a native South American lineage. Based on the H7 and most of the internal genes' phylogenies, these viruses were also closely related to the ones that caused a highly pathogenic avian influenza outbreak in Chile in 2002. Results from this study help to understand the regional dynamics of influenza outbreaks, highlighting the importance of local native viruses circulating in the natural reservoir hosts.

Phylodynamics and evolutionary epidemiology of African swine fever p72-CVR genes in Eurasia and Africa

African swine fever (ASF) is a complex infectious disease of swine that constitutes devastating impacts on animal health and the world economy. Here, we investigated the evolutionary epidemiology of ASF virus (ASFV) in Eurasia and Africa using the concatenated gene sequences of the viral protein 72 and the central variable region of isolates collected between 1960 and 2015. We used Bayesian phylodynamic models to reconstruct the evolutionary history of the virus, to identify virus population demographics and to quantify dispersal patterns between host species. Results suggest that ASFV exhibited a significantly high evolutionary rate and population growth through time since its divergence in the 18th century from East Africa, with no signs of decline till recent years. This increase corresponds to the growing pig trade activities between continents during the 19th century, and may be attributed to an evolutionary drift that resulted from either continuous circulation or maintenance of the virus within Africa and Eurasia. Furthermore, results implicate wild suids as the ancestral host species (root state posterior probability = 0.87) for ASFV in the early 1700s in Africa. Moreover, results indicate the transmission cycle between wild suids and pigs is an important cycle for ASFV spread and maintenance in pig populations, while ticks are an important natural reservoir that can facilitate ASFV spread and maintenance in wild swine populations. We illustrated the prospects of phylodynamic methods in improving risk-based surveillance, support of effective animal health policies, and epidemic preparedness in countries at high risk of ASFV incursion.

Heterologous Packaging Signals on Segment 4, but Not Segment 6 or Segment 8, Limit Influenza A Virus Reassortment

Influenza A virus (IAV) RNA packaging signals serve to direct the incorporation of IAV gene segments into virus particles, and this process is thought to be mediated by segment-segment interactions. These packaging signals are segment and strain specific, and as such, they have the potential to impact reassortment outcomes between different IAV strains. Our study aimed to quantify the impact of packaging signal mismatch on IAV reassortment using the human seasonal influenza A/Panama/2007/99 (H3N2) and pandemic influenza A/Netherlands/602/2009 (H1N1) viruses. Focusing on the three most divergent segments, we constructed pairs of viruses that encoded identical proteins but differed in the packaging signal regions on a single segment. We then evaluated the frequency with which segments carrying homologous versus heterologous packaging signals were incorporated into reassortant progeny viruses. We found that, when segment 4 (HA) of coinfecting parental viruses was modified, there was a significant preference for the segment containing matched packaging signals relative to the background of the virus. This preference was apparent even when the homologous HA constituted a minority of the HA segment population available in the cell for packaging. Conversely, when segment 6 (NA) or segment 8 (NS) carried modified packaging signals, there was no significant preference for homologous packaging signals. These data suggest that movement of NA and NS segments between the human H3N2 and H1N1 lineages is unlikely to be restricted by packaging signal mismatch, while movement of the HA segment would be more constrained. Our results indicate that the importance of packaging signals in IAV reassortment is segment dependent.IMPORTANCE Influenza A viruses (IAVs) can exchange genes through reassortment. This process contributes to both the highly diverse population of IAVs found in nature and the formation of novel epidemic and pandemic IAV strains. Our study sought to determine the extent to which IAV packaging signal divergence impacts reassortment between seasonal IAVs. Our knowledge in this area is lacking, and insight into the factors that influence IAV reassortment will inform and strengthen ongoing public health efforts to anticipate the emergence of new viruses. We found that the packaging signals on the HA segment, but not the NA or NS segments, restricted IAV reassortment. Thus, the packaging signals of the HA segment could be an important factor in determining the likelihood that two IAV strains of public health interest will undergo reassortment.

Of Ducks and Men: Ecology and Evolution of a Zoonotic Pathogen in a Wild Reservoir Host

A hallmark of disease is that most pathogens are able to infect more than one host species. However, for most pathogens, we still have a limited understanding of how this affects epidemiology, persistence and virulence of infections—including several zoonotic pathogens that reside in wild animal reservoirs and spillover into humans. In this chapter, we review the current knowledge of mallard (Anas platyrhynchos) as host for pathogens. This species is widely distributed, often occupying habitats close to humans and livestock, and is an important game bird species and the ancestor to domestic ducks—thereby being an excellent model species to highlight aspects of the wildlife, domestic animal interface and the relevance for human health. We discuss mallard as host for a range of pathogens but focus more in depth of it as a reservoir host for influenza A virus (IAV). Over the last decades, IAV research has surged, prompted in part to the genesis and spread of highly pathogenic virus variants that have been devastating to domestic poultry and caused a number of human spillover infections. The aim of this chapter is to synthesise and review the intricate interactions of virus, host and environmental factors governing IAV epidemiology and evolution.

Risk of resistant avian influenza A virus in wild waterfowl as a result of environmental release of oseltamivir

Oseltamivir is the best available anti-influenza drug and has therefore been stockpiled worldwide in large quantities as part of influenza pandemic preparedness planning. The active metabolite oseltamivir carboxylate (OC) is stable and is not removed by conventional sewage treatment. Active OC has been detected in river water at concentrations up to 0.86 µg/L. Although the natural reservoir hosts of influenza A virus (IAV) are wild waterfowl that reside in aquatic environments, the ecologic risks associated with environmental OC release and its potential to generate resistant viral variants among wild birds has largely been unknown. However, in recent years a number of in vivo mallard (Anas platyrhynchos) studies have been conducted regarding the potential of avian IAVs to become resistant to OC in natural reservoir birds if these are drug exposed. Development of resistance to OC was observed both in Group 1 (N1) and Group 2 (N2, N9) neuraminidase subtypes, when infected ducks were exposed to OC at concentrations between 0.95 and 12 µg/L in their water. All resistant variants maintained replication and transmission between ducks during drug exposure. In an A(H1N1)/H274Y virus, the OC resistance mutation persisted without selective drug pressure, demonstrating the potential of an IAV with a permissive genetic background to acquire and maintain OC resistance, potentially allowing circulation of the resistant variant among wild birds. The experimental studies have improved the appreciation of the risks associated with the environmental release of OC related to resistance development of avian IAVs among wild birds. Combined with knowledge of efficient methods for improved sewage treatment, the observations warrant implementation of novel efficient wastewater treatment methods, rational use of anti-influenza drugs, and improved surveillance of IAV resistance in wild birds.

Antiviral drug resistance as an adaptive process

Antiviral drug resistance is a matter of great clinical importance that, historically, has been investigated mostly from a virological perspective. Although the proximate mechanisms of resistance can be readily uncovered using these methods, larger evolutionary trends often remain elusive. Recent interest by population geneticists in studies of antiviral resistance has spurred new metrics for evaluating mutation and recombination rates, demographic histories of transmission and compartmentalization, and selective forces incurred during viral adaptation to antiviral drug treatment. We present up-to-date summaries on antiviral resistance for a range of drugs and viral types, and review recent advances for studying their evolutionary histories. We conclude that information imparted by demographic and selective histories, as revealed through population genomic inference, is integral to assessing the evolution of antiviral resistance as it pertains to human health.

Ecosystem Interactions Underlie the Spread of Avian Influenza A Viruses with Pandemic Potential

Despite evidence for avian influenza A virus (AIV) transmission between wild and domestic ecosystems, the roles of bird migration and poultry trade in the spread of viruses remain enigmatic. In this study, we integrate ecosystem interactions into a phylogeographic model to assess the contribution of wild and domestic hosts to AIV distribution and persistence. Analysis of globally sampled AIV datasets shows frequent two-way transmission between wild and domestic ecosystems. In general, viral flow from domestic to wild bird populations was restricted to within a geographic region. In contrast, spillover from wild to domestic populations occurred both within and between regions. Wild birds mediated long-distance dispersal at intercontinental scales whereas viral spread among poultry populations was a major driver of regional spread. Viral spread between poultry flocks frequently originated from persistent lineages circulating in regions of intensive poultry production. Our analysis of long-term surveillance data demonstrates that meaningful insights can be inferred from integrating ecosystem into phylogeographic reconstructions that may be consequential for pandemic preparedness and livestock protection.

Host-Specific and Segment-Specific Evolutionary Dynamics of Avian and Human Influenza A Viruses: A Systematic Review

Understanding the evolutionary dynamics of influenza viruses is essential to control both avian and human influenza. Here, we analyze host-specific and segment-specific Tajima's D trends of influenza A virus through a systematic review using viral sequences registered in the National Center for Biotechnology Information. To avoid bias from viral population subdivision, viral sequences were stratified according to their sampling locations and sampling years. As a result, we obtained a total of 580 datasets each of which consists of nucleotide sequences of influenza A viruses isolated from a single population of hosts at a single sampling site within a single year. By analyzing nucleotide sequences in the datasets, we found that Tajima's D values of viral sequences were different depending on hosts and gene segments. Tajima's D values of viruses isolated from chicken and human samples showed negative, suggesting purifying selection or a rapid population growth of the viruses. The negative Tajima's D values in rapidly growing viral population were also observed in computer simulations. Tajima's D values of PB2, PB1, PA, NP, and M genes of the viruses circulating in wild mallards were close to zero, suggesting that these genes have undergone neutral selection in constant-sized population. On the other hand, Tajima's D values of HA and NA genes of these viruses were positive, indicating HA and NA have undergone balancing selection in wild mallards. Taken together, these results indicated the existence of unknown factors that maintain viral subtypes in wild mallards.

Coupling high-throughput genetics with phylogenetic information reveals an epistatic interaction on the influenza A virus M segment

Background

Epistasis is one of the central themes in viral evolution due to its importance in drug resistance, immune escape, and interspecies transmission. However, there is a lack of experimental approach to systematically probe for epistatic residues.

Results

By utilizing the information from natural occurring sequences and high-throughput genetics, this study established a novel strategy to identify epistatic residues. The rationale is that a substitution that is deleterious in one strain may be prevalent in nature due to the presence of a naturally occurring compensatory substitution. Here, high-throughput genetics was applied to influenza A virus M segment to systematically identify deleterious substitutions. Comparison with natural sequence variation showed that a deleterious substitution M1 Q214H was prevalent in circulating strains. A coevolution analysis was then performed and indicated that M1 residues 121, 207, 209, and 214 naturally coevolved as a group. Subsequently, we experimentally validated that M1 A209T was a compensatory substitution for M1 Q214H.

Conclusions

This work provided a proof-of-concept to identify epistatic residues by coupling high-throughput genetics with phylogenetic information. In particular, we were able to identify an epistatic interaction between M1 substitutions A209T and Q214H. This analytic strategy can potentially be adapted to study any protein of interest, provided that the information on natural sequence variants is available.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2358-7) contains supplementary material, which is available to authorized users.

Genome plasticity of triple-reassortant H1N1 influenza A virus during infection of vaccinated pigs

To gain insight into the evolution of influenza A viruses (IAVs) during infection of vaccinated pigs, we experimentally infected a 3-week-old naive pig with a triple-reassortant H1N1 IAV and placed the seeder pig in direct contact with a group of age-matched vaccinated pigs (n = 10). We indexed the genetic diversity and evolution of the virus at an intra-host level by deep sequencing the entire genome directly from nasal swabs collected at two separate samplings during infection. We obtained 13 IAV metagenomes from 13 samples, which included the virus inoculum and two samples from each of the six pigs that tested positive for IAV during the study. The infection produced a population of heterogeneous alleles (sequence variants) that was dynamic over time. Overall, 794 polymorphisms were identified amongst all samples, which yielded 327 alleles, 214 of which were unique sequences. A total of 43 distinct haemagglutinin proteins were translated, two of which were observed in multiple pigs, whereas the neuraminidase (NA) was conserved and only one dominant NA was found throughout the study. The genetic diversity of IAVs changed dynamically within and between pigs. However, most of the substitutions observed in the internal gene segments were synonymous. Our results demonstrated remarkable IAV diversity, and the complex, rapid and dynamic evolution of IAV during infection of vaccinated pigs that can only be appreciated with repeated sampling of individual animals and deep sequence analysis.

Spatiotemporal dynamics of Puumala hantavirus associated with its rodent host, Myodes glareolus

Many viruses significantly impact human and animal health. Understanding the population dynamics of these viruses and their hosts can provide important insights for epidemiology and virus evolution. Puumala virus (PUUV) is a European hantavirus that may cause regional outbreaks of hemorrhagic fever with renal syndrome in humans. Here, we analyzed the spatiotemporal dynamics of PUUV circulating in local populations of its rodent reservoir host, the bank vole (Myodes glareolus) during eight years. Phylogenetic and population genetic analyses of all three genome segments of PUUV showed strong geographical structuring at a very local scale. There was a high temporal turnover of virus strains in the local bank vole populations, but several virus strains persisted through multiple years. Phylodynamic analyses showed no significant changes in the local effective population sizes of PUUV, although vole numbers and virus prevalence fluctuated widely. Microsatellite data demonstrated also a temporally persisting subdivision between local vole populations, but these groups did not correspond to the subdivision in the virus strains. We conclude that restricted transmission between vole populations and genetic drift play important roles in shaping the genetic structure and temporal dynamics of PUUV in its natural host which has several implications for zoonotic risks of the human population.

Avian influenza virus exhibits distinct evolutionary dynamics in wild birds and poultry

BACKGROUND

Wild birds are the major reservoir hosts for influenza A viruses, occasionally transmitting to other species such as domesticated poultry. Despite an abundance of genomic data from avian influenza virus (AIV), little is known about whether AIV evolves differently in wild birds and poultry, although this is critical to revealing the dynamics and time-scale of viral evolution. In particular, because environmental (water-borne) transmission is more common in wild birds, which may reduce the number of replications per unit time, it is possible that evolutionary rates are systematically lower in wild birds than in poultry.

RESULTS

We estimated rates of nucleotide substitution in two AIV subtypes that are strongly associated with infections in wild birds - H4 and H6 - and compared these to rates in the H5N1 subtype that has circulated in poultry for almost two decades. Our analyses of three internal genes confirm that H4 and H6 viruses are evolving significantly more slowly than H5N1 viruses, suggesting that evolutionary rates of AIV are reduced in wild birds. This result was verified by the analysis of a poultry-associated H6 lineage that exhibited a markedly higher substitution rate than those H6 viruses circulating in wild birds. Interestingly, we also observed a significant difference in evolutionary rate between H4 and H6, despite frequent reassortment rate among them.

CONCLUSIONS

AIV experiences markedly different evolutionary dynamics between wild birds and poultry. These results suggest that rate heterogeneity among viral subtypes and ecological groupings should be taken into account when estimating evolutionary rates and divergence times.

Association between Influenza A Virus Infection and Pigs Subpopulations in Endemically Infected Breeding Herds

Influenza A viruses (IAVs) are distributed worldwide in birds, pigs and humans, and cause important endemic disease affecting hosts in all countries. Although pigs play a key role in the ecology of IAVs, the epidemiology of IAVs within swine herds is poorly understood. In this longitudinal study we describe the prevalence of IAVs infection in three subpopulations of pigs in 5 breeding herds in the Midwestern USA. Each herd was sampled monthly for a year and, at each visit, 30 individual nasal swabs were collected from the three subpopulations, namely, a) replacement females, resident on-farm for less than 4 weeks (new gilts), b) replacement females, resident on-farm for more than 4 weeks (gilts), and c) neonatal pigs less than 21 days of age (piglets). Real time reverse transcriptase polymerase chain reaction (RRT-PCR) was used to detect IAVs, and the association between IAVs infection and pig subpopulation was measured using a mixed logistic regression model. Nasal swabs (n = 4,190) were collected from 141 groups of pigs. At least, one IAV-positive nasal swab was found in 19.9% (n = 28) of the sampled groups, and 7.7% (n = 324) of all nasal swabs tested positive. After adjusting by annual quarter and sampling event, the odds of testing IAV positive were 7.9 (95% CI 1.4, 43.9) and 4.4 (95% CI 1.1, 17.1) times higher in groups of new gilts and piglets compared to groups of gilts, respectively. Results indicate that new gilts and piglets had higher odds of testing IAV positive than gilts in swine breeding herds and that season influences IAV infection in pigs. Based on these findings, we recommend that IAV control strategies be aimed at preventing infection before gilts are introduced into the farm, and in pigs prior to weaning.

Spread and persistence of influenza A viruses in waterfowl hosts in the North American Mississippi migratory flyway

While geographic distance often restricts the spread of pathogens via hosts, this barrier may be compromised when host species are mobile. Migratory waterfowl in the order Anseriformes are important reservoir hosts for diverse populations of avian-origin influenza A viruses (AIVs) and are assumed to spread AIVs during their annual continental-scale migrations. However, support for this hypothesis is limited, and it is rarely tested using data from comprehensive surveillance efforts incorporating both the temporal and spatial aspects of host migratory patterns. We conducted intensive AIV surveillance of waterfowl using the North American Mississippi Migratory Flyway (MMF) over three autumn migratory seasons. Viral isolates (n = 297) from multiple host species were sequenced and analyzed for patterns of gene dispersal between northern staging and southern wintering locations. Using a phylogenetic and nucleotide identity framework, we observed a larger amount of gene dispersal within this flyway rather than between the other three longitudinally identified North American flyways. Across seasons, we observed patterns of regional persistence of diversity for each genomic segment, along with limited survival of dispersed AIV gene lineages. Reassortment increased with both time and distance, resulting in transient AIV constellations. This study shows that within the MMF, AIV gene flow favors spread along the migratory corridor within a season, and also that intensive surveillance during bird migration is important for identifying virus dispersal on time scales relevant to pandemic responsiveness. In addition, this study indicates that comprehensive monitoring programs to capture AIV diversity are critical for providing insight into AIV evolution and ecology in a major natural reservoir.

IMPORTANCE

Migratory birds are a reservoir for antigenic and genetic diversity of influenza A viruses (AIVs) and are implicated in the spread of virus diversity that has contributed to previous pandemic events. Evidence for dispersal of avian-origin AIVs by migratory birds is rarely examined on temporal scales relevant to pandemic or panzootic threats. Therefore, characterizing AIV movement by hosts within a migratory season is important for implementing effective surveillance strategies. We conducted surveillance following birds along a major North American migratory route and observed that within a migratory season, AIVs rapidly reassorted and gene lineages were dispersed primarily within the migratory corridor. Patterns of regional persistence were observed across seasons for each gene segment. We show that dispersal of AIV gene lineages by migratory birds occurs quickly along migratory routes and that surveillance for AIVs threatening human and animal health should focus attention on these routes.

Influenza A virus evolution and spatio-temporal dynamics in Eurasian wild birds: a phylogenetic and phylogeographical study of whole-genome sequence data

Low pathogenic avian influenza A viruses (IAVs) have a natural host reservoir in wild waterbirds and the potential to spread to other host species. Here, we investigated the evolutionary, spatial and temporal dynamics of avian IAVs in Eurasian wild birds. We used whole-genome sequences collected as part of an intensive long-term Eurasian wild bird surveillance study, and combined this genetic data with temporal and spatial information to explore the virus evolutionary dynamics. Frequent reassortment and co-circulating lineages were observed for all eight genomic RNA segments over time. There was no apparent species-specific effect on the diversity of the avian IAVs. There was a spatial and temporal relationship between the Eurasian sequences and significant viral migration of avian IAVs from West Eurasia towards Central Eurasia. The observed viral migration patterns differed between segments. Furthermore, we discuss the challenges faced when analysing these surveillance and sequence data, and the caveats to be borne in mind when drawing conclusions from the apparent results of such analyses.

Epistatically interacting substitutions are enriched during adaptive protein evolution

Most experimental studies of epistasis in evolution have focused on adaptive changes—but adaptation accounts for only a portion of total evolutionary change. Are the patterns of epistasis during adaptation representative of evolution more broadly? We address this question by examining a pair of protein homologs, of which only one is subject to a well-defined pressure for adaptive change. Specifically, we compare the nucleoproteins from human and swine influenza. Human influenza is under continual selection to evade recognition by acquired immune memory, while swine influenza experiences less such selection due to the fact that pigs are less likely to be infected with influenza repeatedly in a lifetime. Mutations in some types of immune epitopes are therefore much more strongly adaptive to human than swine influenza—here we focus on epitopes targeted by human cytotoxic T lymphocytes. The nucleoproteins of human and swine influenza possess nearly identical numbers of such epitopes. However, mutations in these epitopes are fixed significantly more frequently in human than in swine influenza, presumably because these epitope mutations are adaptive only to human influenza. Experimentally, we find that epistatically constrained mutations are fixed only in the adaptively evolving human influenza lineage, where they occur at sites that are enriched in epitopes. Overall, our results demonstrate that epistatically interacting substitutions are enriched during adaptation, suggesting that the prevalence of epistasis is dependent on the underlying evolutionary forces at play.

Mutations can fix during evolution for two reasons: they can be beneficial and fix for adaptive reasons, or they can be neutral or deleterious and fix solely by chance. Most studies focus on adaptation, where the evolving population is increasing in fitness due to a new selection pressure. Such studies have found an important evolutionary role for epistasis, the phenomenon where the effect of one mutation depends on another mutation. But adaptation only accounts for a fraction of overall evolutionary change. Here we investigate whether epistasis is as common during non-adaptive as adaptive evolution. We do this by comparing the same protein from human and swine influenza. Human influenza is constantly adapting to escape from the immunity that people acquire from previous influenza infections. But swine influenza is under less pressure to escape from acquired immunity since pigs have shorter lifetimes and are less likely to be infected with influenza multiple times. We find that epistasis is less common during the evolution of the swine influenza protein than its human influenza counterpart. Overall, our results suggest that mutations that interact via epistasis are more likely to fix during adaptive evolution.

A synchronized global sweep of the internal genes of modern avian influenza virus

Zoonotic infectious diseases such as influenza continue to pose a grave threat to human health. However, the factors that mediate the emergence of RNA viruses such as influenza A virus (IAV) are still incompletely understood. Phylogenetic inference is crucial to reconstructing the origins and tracing the flow of IAV within and between hosts. Here we show that explicitly allowing IAV host lineages to have independent rates of molecular evolution is necessary for reliable phylogenetic inference of IAV and that methods that do not do so, including 'relaxed' molecular clock models, can be positively misleading. A phylogenomic analysis using a host-specific local clock model recovers extremely consistent evolutionary histories across all genomic segments and demonstrates that the equine H7N7 lineage is a sister clade to strains from birds--as well as those from humans, swine and the equine H3N8 lineage--sharing an ancestor with them in the mid to late 1800s. Moreover, major western and eastern hemisphere avian influenza lineages inferred for each gene coalesce in the late 1800s. On the basis of these phylogenies and the synchrony of these key nodes, we infer that the internal genes of avian influenza virus (AIV) underwent a global selective sweep beginning in the late 1800s, a process that continued throughout the twentieth century and up to the present. The resulting western hemispheric AIV lineage subsequently contributed most of the genomic segments to the 1918 pandemic virus and, independently, the 1963 equine H3N8 panzootic lineage. This approach provides a clear resolution of evolutionary patterns and processes in IAV, including the flow of viral genes and genomes within and between host lineages.

A 4-year study of avian influenza virus prevalence and subtype diversity in ducks of Newfoundland, Canada

The island of Newfoundland, Canada, is at the eastern edge of North America and has migratory bird connections with the continental mainland as well as across the North Atlantic Ocean. Here, we report a 4-year avian influenza virus (AIV) epidemiological study in ducks in the St. John's region of Newfoundland. The overall prevalence of AIV detection in ducks during this study was 7.2%, with American Black Ducks contributing the vast majority of the collected samples and the AIV positives. The juvenile ducks showed a significantly higher AIV detection rate (10.6%) compared with adults (3.4%). Seasonally, AIV prevalence rates were higher in the autumn (8.4%), but positives were still detected in the winter (4.6%). Preliminary serology tests showed a high incidence of previous AIV infection (20/38, 52.6%). A total of 43 viruses were characterized for their HA-NA or HA subtypes, which revealed a large diversity of AIV subtypes and little recurrence of subtypes from year to year. Investigation of the movement patterns of ducks in this region showed that it is a largely non-migratory duck population, which may contribute to the observed pattern of high AIV subtype turnover. Phylogenetic analysis of 4 H1N1 and one H5N4 AIVs showed these viruses were highly similar to other low pathogenic AIV sequences from waterfowl in North America and assigned all gene segments into American-avian clades. Notably, the H1N1 viruses, which were identified in consecutive years, possessed homologous genomes. Such detection of homologous AIV genomes across years is rare, but indicates the role of the environmental reservoir in viral perpetuation.

Protein clustering and RNA phylogenetic reconstruction of the influenza A [corrected] virus NS1 protein allow an update in classification and identification of motif conservation

The non-structural protein 1 (NS1) of influenza A virus (IAV), coded by its third most diverse gene, interacts with multiple molecules within infected cells. NS1 is involved in host immune response regulation and is a potential contributor to the virus host range. Early phylogenetic analyses using 50 sequences led to the classification of NS1 gene variants into groups (alleles) A and B. We reanalyzed NS1 diversity using 14,716 complete NS IAV sequences, downloaded from public databases, without host bias. Removal of sequence redundancy and further structured clustering at 96.8% amino acid similarity produced 415 clusters that enhanced our capability to detect distinct subgroups and lineages, which were assigned a numerical nomenclature. Maximum likelihood phylogenetic reconstruction using RNA sequences indicated the previously identified deep branching separating group A from group B, with five distinct subgroups within A as well as two and five lineages within the A4 and A5 subgroups, respectively. Our classification model proposes that sequence patterns in thirteen amino acid positions are sufficient to fit >99.9% of all currently available NS1 sequences into the A subgroups/lineages or the B group. This classification reduces host and virus bias through the prioritization of NS1 RNA phylogenetics over host or virus phenetics. We found significant sequence conservation within the subgroups and lineages with characteristic patterns of functional motifs, such as the differential binding of CPSF30 and crk/crkL or the availability of a C-terminal PDZ-binding motif. To understand selection pressures and evolution acting on NS1, it is necessary to organize the available data. This updated classification may help to clarify and organize the study of NS1 interactions and pathogenic differences and allow the drawing of further functional inferences on sequences in each group, subgroup and lineage rather than on a strain-by-strain basis.

A case for the ancient origin of coronaviruses

Coronaviruses are found in a diverse array of bat and bird species, which are believed to act as natural hosts. Molecular clock dating analyses of coronaviruses suggest that the most recent common ancestor of these viruses existed around 10,000 years ago. This relatively young age is in sharp contrast to the ancient evolutionary history of their putative natural hosts, which began diversifying tens of millions of years ago. Here, we attempted to resolve this discrepancy by applying more realistic evolutionary models that have previously revealed the ancient evolutionary history of other RNA viruses. By explicitly modeling variation in the strength of natural selection over time and thereby improving the modeling of substitution saturation, we found that the time to the most recent ancestor common for all coronaviruses is likely far greater (millions of years) than the previously inferred range.

Reassortant influenza A viruses in wild duck populations: effects on viral shedding and persistence in water

Wild ducks of the genus Anas represent the natural hosts for a large genetic diversity of influenza A viruses. In these hosts, co-infections with different virus genotypes are frequent and result in high rates of genetic reassortment. Recent genomic data have provided information regarding the pattern and frequency of these reassortant viruses in duck populations; however, potential consequences on viral shedding and maintenance in the environment have not been investigated. On the basis of full-genome sequencing, we identified five virus genotypes, in a wild duck population in northwestern Minnesota (USA), that naturally arose from genetic reassortments. We investigated the effects of influenza A virus genotype on the viral shedding pattern in Mallards (Anas platyrhynchos) and the duration of infectivity in water, under different temperature regimens. Overall, we found that variation in the viral genome composition of these isolates had limited effects on duration, extent and pattern of viral shedding, as well as on the reduction of infectivity in water over time. These results support that, in wild ducks, functionally equivalent gene segments could be maintained in virus populations with no fitness costs when genetic reassortments occur.

A distinct lineage of influenza A virus from bats

Influenza A virus reservoirs in animals have provided novel genetic elements leading to the emergence of global pandemics in humans. Most influenza A viruses circulate in waterfowl, but those that infect mammalian hosts are thought to pose the greatest risk for zoonotic spread to humans and the generation of pandemic or panzootic viruses. We have identified an influenza A virus from little yellow-shouldered bats captured at two locations in Guatemala. It is significantly divergent from known influenza A viruses. The HA of the bat virus was estimated to have diverged at roughly the same time as the known subtypes of HA and was designated as H17. The neuraminidase (NA) gene is highly divergent from all known influenza NAs, and the internal genes from the bat virus diverged from those of known influenza A viruses before the estimated divergence of the known influenza A internal gene lineages. Attempts to propagate this virus in cell cultures and chicken embryos were unsuccessful, suggesting distinct requirements compared with known influenza viruses. Despite its divergence from known influenza A viruses, the bat virus is compatible for genetic exchange with human influenza viruses in human cells, suggesting the potential capability for reassortment and contributions to new pandemic or panzootic influenza A viruses.

Genetic structure of Pacific Flyway avian influenza viruses is shaped by geographic location, host species, and sampling period

The eight gene segments of avian influenza virus (AIV) reassort frequently and rapidly to generate novel genotypes and subtypes that are transmissible to a broad range of hosts. There is evidence that AIV can have a restricted host range and can segregate in space and time. Host-virus relationships at the species, geographic, and spatial scales have not been fully defined for AIV populations of the Pacific Flyway, particularly among the diverse waterfowl that occupy the Flyway in Alaska and California. Using the sequence analysis program Bayesian Tip-association Significance testing (BaTS) created for analysis of phylogeny-trait associations, we determined whether the genetic structure of Pacific Flyway AIVs sampled between 2006 and 2008 was influenced by the host species, geographic location of virus collection, and time of sampling. In posterior sets of trees, genetically similar viruses clustered by host species for thick-billed murres and glaucous gulls (order Charadriiformes), and for northern shovelers, northern pintails, and mallards (order Anseriformes). AIVs from Alaska and California were strongly spatially structured, clustering separately by region across all segments. The timing of sampling influenced the genetic structure of California AIV gene segments, possibly reflecting waves of host species movement into wintering areas. The strength of phylogeny-trait association varied by virus segment and by trait of interest, which we hypothesize is related to the frequent genetic reassortment and interspecies transmission in waterfowl.

Purifying selection can obscure the ancient age of viral lineages

Statistical methods for molecular dating of viral origins have been used extensively to infer the time of most common recent ancestor for many rapidly evolving pathogens. However, there are a number of cases, in which epidemiological, historical, or genomic evidence suggests much older viral origins than those obtained via molecular dating. We demonstrate how pervasive purifying selection can mask the ancient origins of recently sampled pathogens, in part due to the inability of nucleotide-based substitution models to properly account for complex patterns of spatial and temporal variability in selective pressures. We use codon-based substitution models to infer the length of branches in viral phylogenies; these models produce estimates that are often considerably longer than those obtained with traditional nucleotide-based substitution models. Correcting the apparent underestimation of branch lengths suggests substantially older origins for measles, Ebola, and avian influenza viruses. This work helps to reconcile some of the inconsistencies between molecular dating and other types of evidence concerning the age of viral lineages.

Evolutionary dynamics of influenza A nucleoprotein (NP) lineages revealed by large-scale sequence analyses

Influenza A viral nucleoprotein (NP) plays a critical role in virus replication and host adaptation, however, the underlying molecular evolutionary dynamics of NP lineages are less well-understood. In this study, large-scale analyses of 5094 NP nucleotide sequences revealed eight distinct evolutionary lineages, including three host-specific lineages (human, classical swine and equine), two cross-host lineages (Eurasian avian-like swine and swine-origin human pandemic H1N1 2009) and three geographically isolated avian lineages (Eurasian, North American and Oceanian). The average nucleotide substitution rate of the NP lineages was estimated to be 2.4 × 10(-3) substitutions per site per year, with the highest value observed in pandemic H1N1 2009 (3.4 × 10(-3)) and the lowest in equine (0.9 × 10(-3)). The estimated time of most recent common ancestor (TMRCA) for each lineage demonstrated that the earliest human lineage was derived around 1906, and the latest pandemic H1N1 2009 lineage dated back to December 17, 2008. A marked time gap was found between the times when the viruses emerged and were first sampled, suggesting the crucial role for long-term surveillance of newly emerging viruses. The selection analyses showed that human lineage had six positive selection sites, whereas pandemic H1N1 2009, classical swine, Eurasian avian and Eurasian swine had only one or two sites. Protein structure analyses revealed several positive selection sites located in epitope regions or host adaptation regions, indicating strong adaptation to host immune system pressures in influenza viruses. Along with previous studies, this study provides new insights into the evolutionary dynamics of influenza A NP lineages. Further lineage analyses of other gene segments will allow better understanding of influenza A virus evolution and assist in the improvement of global influenza surveillance.

Extensive geographic mosaicism in avian influenza viruses from gulls in the northern hemisphere

Due to limited interaction of migratory birds between Eurasia and America, two independent avian influenza virus (AIV) gene pools have evolved. There is evidence of low frequency reassortment between these regions, which has major implications in global AIV dynamics. Indeed, all currently circulating lineages of the PB1 and PA segments in North America are of Eurasian origin. Large-scale analyses of intercontinental reassortment have shown that viruses isolated from Charadriiformes (gulls, terns, and shorebirds) are the major contributor of these outsider events. To clarify the role of gulls in AIV dynamics, specifically in movement of genes between geographic regions, we have sequenced six gull AIV isolated in Alaska and analyzed these along with 142 other available gull virus sequences. Basic investigations of host species and the locations and times of isolation reveal biases in the available sequence information. Despite these biases, our analyses reveal a high frequency of geographic reassortment in gull viruses isolated in America. This intercontinental gene mixing is not found in the viruses isolated from gulls in Eurasia. This study demonstrates that gulls are important as vectors for geographically reassorted viruses, particularly in America, and that more surveillance effort should be placed on this group of birds.

The nonadaptive nature of the H1N1 2009 Swine Flu pandemic contrasts with the adaptive facilitation of transmission to a new host

BACKGROUND

The emergence of the 2009 H1N1 Influenza pandemic followed a multiple reassortment event from viruses originally circulating in swines and humans, but the adaptive nature of this emergence is poorly understood.

RESULTS

Here we base our analysis on 1180 complete genomes of H1N1 viruses sampled in North America between 2000 and 2010 in swine and human hosts. We show that while transmission to a human host might require an adaptive phase in the HA and NA antigens, the emergence of the 2009 pandemic was essentially nonadaptive. A more detailed analysis of the NA protein shows that the 2009 pandemic sequence is characterized by novel epitopes and by a particular substitution in loop 150, which is responsible for a nonadaptive structural change tightly associated with the emergence of the pandemic.

CONCLUSIONS

Because this substitution was not present in the 1918 H1N1 pandemic virus, we posit that the emergence of pandemics is due to epistatic interactions between sites distributed over different segments. Altogether, our results are consistent with population dynamics models that highlight the epistatic and nonadaptive rise of novel epitopes in viral populations, followed by their demise when the resulting virus is too virulent.

Influenza virus evolution, host adaptation, and pandemic formation

Newly emerging or "re-emerging" viral diseases continue to pose significant global public health threats. Prototypic are influenza viruses that are major causes of human respiratory infections and mortality. Influenza viruses can cause zoonotic infections and adapt to humans, leading to sustained transmission and emergence of novel viruses. Mechanisms by which viruses evolve in one host, cause zoonotic infection, and adapt to a new host species remain unelucidated. Here, we review the evolution of influenza A viruses in their reservoir hosts and discuss genetic changes associated with introduction of novel viruses into humans, leading to pandemics and the establishment of seasonal viruses.