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

Avian influenza viruses in New Zealand wild birds, with an emphasis on subtypes H5 and H7: Their distinctive epidemiology and genomic properties

The rapid spread of highly pathogenic avian influenza (HPAI) A (H5N1) viruses in Southeast Asia in 2004 prompted the New Zealand Ministry for Primary Industries to expand its avian influenza surveillance in wild birds. A total of 18,693 birds were sampled between 2004 and 2020, including migratory shorebirds (in 2004-2009), other coastal species (in 2009-2010), and resident waterfowl (in 2004-2020). No avian influenza viruses (AIVs) were isolated from cloacal or oropharyngeal samples from migratory shorebirds or resident coastal species. Two samples from red knots (Calidris canutus) tested positive by influenza A RT-qPCR, but virus could not be isolated and no further characterization could be undertaken. In contrast, 6179 samples from 15,740 mallards (Anas platyrhynchos) tested positive by influenza A RT-qPCR. Of these, 344 were positive for H5 and 51 for H7. All H5 and H7 viruses detected were of low pathogenicity confirmed by a lack of multiple basic amino acids at the hemagglutinin (HA) cleavage site. Twenty H5 viruses (six different neuraminidase [NA] subtypes) and 10 H7 viruses (two different NA subtypes) were propagated and characterized genetically. From H5- or H7-negative samples that tested positive by influenza A RT-qPCR, 326 AIVs were isolated, representing 41 HA/NA combinations. The most frequently isolated subtypes were H4N6, H3N8, H3N2, and H10N3. Multivariable logistic regression analysis of the relations between the location and year of sampling, and presence of AIV in individual waterfowl showed that the AIV risk at a given location varied from year to year. The H5 and H7 isolates both formed monophyletic HA groups. The H5 viruses were most closely related to North American lineages, whereas the H7 viruses formed a sister cluster relationship with wild bird viruses of the Eurasian and Australian lineages. Bayesian analysis indicates that the H5 and H7 viruses have circulated in resident mallards in New Zealand for some time. Correspondingly, we found limited evidence of influenza viruses in the major migratory bird populations visiting New Zealand. Findings suggest a low probability of introduction of HPAI viruses via long-distance bird migration and a unique epidemiology of AIV in New Zealand.

The H4 subtype of avian influenza virus: a review of its historical evolution, global distribution, adaptive mutations and receptor binding properties

The H4 subtype of avian influenza virus (AIV) exhibits a wide host range and is commonly found in migratory waterfowl. Recent studies have revealed that the H4N6 AIV can infect guinea pigs via aerosol transmission without prior adaptation. Additionally, the Q226L/G228S substitutions in the receptor-binding site have led to structural changes in globular head of H4 AIV, resulting in a configuration similar to that of pandemic H2N2 and H3N2 human influenza viruses. This article provides an updated review of the historical evolution, global distribution, adaptive mutations, receptor-binding preferences, and host range of H4 AIV. The insights presented herein will help in assessing the potential risk of future H4 AIV epidemics.

South American H4N2 influenza A virus improved replication in chicken trachea after low number of passages

Introduction of influenza A viruses (FLUAV) into poultry from waterfowl is frequent, producing economic burden and increasing the probability of human infections. We have previously described the presence of FLUAV in wild birds in Argentina with unique evolutionary trajectories belonging to a South American lineage different from the North American and Eurasian lineages. Adaptability of this South American lineage FLUAV to poultry species is still poorly understood. In the present report, we evaluated the capacity of an H4N2 FLUAV from the South American lineage to adapt to chickens after low number of passages. We found that five mutations were acquired after five passages in 3-days-old chickens. These mutations produced a virus with better infectivity in ex vivo trachea explants but overall lower infection in lung explants. Infection of 3-week-old chickens persisted for a longer period and was detected in more tissues than the parental virus, suggesting adaptation of the H4N2 influenza A virus to chicken.

Evolutionary Dynamics of Avian Influenza Viruses Isolated from Wild Birds in Moscow

Forty-five strains of AIVs were isolated from wild aquatic birds during their autumn migration through Moscow (Russia). The aim of this work is to study the dynamics of AIV genomes in their natural habitat. Viruses were isolated from fecal sample in embryonated chicken eggs; their complete genomes were sequenced, and a phylogenetic analysis was performed. The gene segments of the same lineage persisted over the years in the absence of persistence of complete viral genomes. The genes for internal proteins of the same lineage were often maintained by the viruses over few years; however, they were typically associated with the genes of novel HA and NA subtypes. Although frequent reassortment events were observed for any pair of internal genes, there was no reassortment between HA and NA segments. The differences in the persistence of phylogenetic lineages of surface and internal proteins and the different evolutionary strategy for these two types of genes of AIVs in primary hosts are discussed.

Phylodynamics of Highly Pathogenic Avian Influenza A(H5N1) Virus Circulating in Indonesian Poultry

After its first detection in 1996, the highly pathogenic avian influenza A(H5Nx) virus has spread extensively worldwide. HPAIv A(H5N1) was first detected in Indonesia in 2003 and has been endemic in poultry in this country ever since. However, Indonesia has limited information related to the phylodynamics of HPAIv A(H5N1) in poultry. The present study aimed to increase the understanding of the evolution and temporal dynamics of HPAIv H5N1 in Indonesian poultry between 2003 and 2016. To this end, HPAIv A(H5N1) hemagglutinin sequences of viruses collected from 2003 to 2016 were analyzed using Bayesian evolutionary analysis sampling trees. Results indicated that the common ancestor of Indonesian poultry HPAIv H5N1 arose approximately five years after the common ancestor worldwide of HPAI A(H5Nx). In addition, this study indicated that only two introductions of HPAIv A(H5N1) occurred, after which these viruses continued to evolve due to extensive spread among poultry. Furthermore, this study revealed the divergence of H5N1 clade 2.3.2.1c from H5N1 clade 2.3.2.1b. Both clades 2.3.2.1c and 2.3.2.1b share a common ancestor, clade 1, suggesting that clade 2.3.2.1 originated and diverged from China and other Asian countries. Since there was limited sequence and surveillance data for the HPAIv A(H5N1) from wild birds in Indonesia, the exact role of wild birds in the spread of HPAIv in Indonesia is currently unknown. The evolutionary dynamics of the Indonesian HPAIv A(H5N1) highlight the importance of continuing and improved genomic surveillance and adequate control measures in the different regions of both the poultry and wild birds. Spatial genomic surveillance is useful to take adequate control measures. Therefore, it will help to prevent the future evolution of HPAI A(H5N1) and pandemic threats.

Low Pathogenicity H7N3 Avian Influenza Viruses Have Higher Within-Host Genetic Diversity Than a Closely Related High Pathogenicity H7N3 Virus in Infected Turkeys and Chickens

Within-host viral diversity offers a view into the early stages of viral evolution occurring after a virus infects a host. In recent years, advances in deep sequencing have allowed for routine identification of low-frequency variants, which are important sources of viral genetic diversity and can potentially emerge as a major virus population under certain conditions. We examined within-host viral diversity in turkeys and chickens experimentally infected with closely related H7N3 avian influenza viruses (AIVs), specifically one high pathogenicity AIV (HPAIV) and two low pathogenicity AIV (LPAIVs) with different neuraminidase protein stalk lengths. Consistent with the high mutation rates of AIVs, an abundance of intra-host single nucleotide variants (iSNVs) at low frequencies of 2–10% was observed in all samples collected. Furthermore, a small number of common iSNVs were observed between turkeys and chickens, and between directly inoculated and contact-exposed birds. Notably, the LPAIVs have significantly higher iSNV diversities and frequencies of nonsynonymous changes than the HPAIV in both turkeys and chickens. These findings highlight the dynamics of AIV populations within hosts and the potential impact of genetic changes, including mutations in the hemagglutinin gene that confers the high pathogenicity pathotype, on AIV virus populations and evolution.

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.

Animal Disease Surveillance in the 21st Century: Applications and Robustness of Phylodynamic Methods in Recent U.S. Human-Like H3 Swine Influenza Outbreaks

Emerging and endemic animal viral diseases continue to impose substantial impacts on animal and human health. Most current and past molecular surveillance studies of animal diseases investigated spatio-temporal and evolutionary dynamics of the viruses in a disjointed analytical framework, ignoring many uncertainties and made joint conclusions from both analytical approaches. Phylodynamic methods offer a uniquely integrated platform capable of inferring complex epidemiological and evolutionary processes from the phylogeny of viruses in populations using a single Bayesian statistical framework. In this study, we reviewed and outlined basic concepts and aspects of phylodynamic methods and attempted to summarize essential components of the methodology in one analytical pipeline to facilitate the proper use of the methods by animal health researchers. Also, we challenged the robustness of the posterior evolutionary parameters, inferred by the commonly used phylodynamic models, using hemagglutinin (HA) and polymerase basic 2 (PB2) segments of the currently circulating human-like H3 swine influenza (SI) viruses isolated in the United States and multiple priors. Subsequently, we compared similarities and differences between the posterior parameters inferred from sequence data using multiple phylodynamic models. Our suggested phylodynamic approach attempts to reduce the impact of its inherent limitations to offer less biased and biologically plausible inferences about the pathogen evolutionary characteristics to properly guide intervention activities. We also pinpointed requirements and challenges for integrating phylodynamic methods in routine animal disease surveillance activities.

A Turkey-origin H9N2 Avian Influenza Virus Shows Low Pathogenicity but Different Within-Host Diversity in Experimentally Infected Turkeys, Quail and Ducks

Avian influenza virus (AIV) is a highly diverse and widespread poultry pathogen. Its evolution and adaptation may be affected by multiple host and ecological factors, which are still poorly understood. In the present study, a turkey-origin H9N2 AIV was used as a model to investigate the within-host diversity of the virus in turkeys, quail and ducks in conjunction with the clinical course, shedding and seroconversion. Ten birds were inoculated oculonasally with a dose of 10⁶ EID₅₀ of the virus and monitored for 14 days. Virus shedding, transmission and seroconversion were evaluated, and swabs collected at selected time-points were characterized in deep sequencing to assess virus diversity. In general, the virus showed low pathogenicity for the examined bird species, but differences in shedding patterns, seroconversion and clinical outcome were noted. The highest heterogeneity of the virus population as measured by the number of single nucleotide polymorphisms and Shannon entropy was found in oropharyngeal swabs from quail, followed by turkeys and ducks. This suggests a strong bottleneck was imposed on the virus during replication in ducks, which can be explained by its poor adaptation and stronger selection pressure in waterfowl. The high within-host virus diversity in quail with high level of respiratory shedding and asymptomatic course of infection may contribute to our understanding of the role of quail as an intermediate host for adaptation of AIV to other species of poultry. In contrast, low virus complexity was observed in cloacal swabs, mainly from turkeys, showing that the within-host diversity may vary between different replication sites. Consequences of these observations on the virus evolution and adaptation require further investigation.

Global patterns of avian influenza A (H7): virus evolution and zoonotic threats

Avian influenza viruses (AIVs) continue to impose a negative impact on animal and human health worldwide. In particular, the emergence of highly pathogenic AIV H5 and, more recently, the emergence of low pathogenic AIV H7N9 have led to enormous socioeconomical losses in the poultry industry and resulted in fatal human infections. While H5N1 remains infamous, the number of zoonotic infections with H7N9 has far surpassed those attributed to H5. Despite the clear public health concerns posed by AIV H7, it is unclear why specifically this virus subtype became endemic in poultry and emerged in humans. In this review, we bring together data on global patterns of H7 circulation, evolution and emergence in humans. Specifically, we discuss data from the wild bird reservoir, expansion and epidemiology in poultry, significant increase in their zoonotic potential since 2013 and genesis of highly pathogenic H7. In addition, we analysed available sequence data from an evolutionary perspective, demonstrating patterns of introductions into distinct geographic regions and reassortment dynamics. The integration of all aspects is crucial in the optimisation of surveillance efforts in wild birds, poultry and humans, and we emphasise the need for a One Health approach in controlling emerging viruses such as AIV H7.

Genetic relationship between poultry and wild bird viruses during the highly pathogenic avian influenza H5N6 epidemic in the Netherlands, 2017-2018

In the Netherlands, three commercial poultry farms and two hobby holdings were infected with highly pathogenic avian influenza (HPAI) H5N6 virus in the winter of 2017-2018. This H5N6 virus is a reassortant of HPAI H5N8 clade 2.3.4.4 group B viruses detected in Eurasia in 2016. H5N6 viruses were also detected in several dead wild birds during the winter. However, wild bird mortality was limited compared to the caused by the H5N8 group B virus in 2016-2017. H5N6 virus was not detected in wild birds after March, but in late summer infected wild birds were found again. In this study, the complete genome sequences of poultry and wild bird viruses were determined to study their genetic relationship. Genetic analysis showed that the outbreaks in poultry were not the result of farm-to-farm transmissions, but rather resulted from separate introductions from wild birds. Wild birds infected with viruses related to the first outbreak in poultry were found at short distances from the farm, within a short time frame. However, no wild bird viruses related to outbreaks 2 and 3 were detected. The H5N6 virus isolated in summer shares a common ancestor with the virus detected in outbreak 1. This suggests long-term circulation of H5N6 virus in the local wild bird population. In addition, the pathogenicity of H5N6 virus in ducks was determined, and compared to that of H5N8 viruses detected in 2014 and 2016. A similar high pathogenicity was measured for H5N6 and H5N8 group B viruses, suggesting that biological or ecological factors in the wild bird population may have affected the mortality rates during the H5N6 epidemic. These observations suggest different infection dynamics for the H5N6 and H5N8 group B viruses in the wild bird population.

Zoonotic Potential of Influenza A Viruses: A Comprehensive Overview

Influenza A viruses (IAVs) possess a great zoonotic potential as they are able to infect different avian and mammalian animal hosts, from which they can be transmitted to humans. This is based on the ability of IAV to gradually change their genome by mutation or even reassemble their genome segments during co-infection of the host cell with different IAV strains, resulting in a high genetic diversity. Variants of circulating or newly emerging IAVs continue to trigger global health threats annually for both humans and animals. Here, we provide an introduction on IAVs, highlighting the mechanisms of viral evolution, the host spectrum, and the animal/human interface. Pathogenicity determinants of IAVs in mammals, with special emphasis on newly emerging IAVs with pandemic potential, are discussed. Finally, an overview is provided on various approaches for the prevention of human IAV infections.

Local and relaxed clocks: the best of both worlds

Time-resolved phylogenetic methods use information about the time of sample collection to estimate the rate of evolution. Originally, the models used to estimate evolutionary rates were quite simple, assuming that all lineages evolve at the same rate, an assumption commonly known as the molecular clock. Richer and more complex models have since been introduced to capture the phenomenon of substitution rate variation among lineages. Two well known model extensions are the local clock, wherein all lineages in a clade share a common substitution rate, and the uncorrelated relaxed clock, wherein the substitution rate on each lineage is independent from other lineages while being constrained to fit some parametric distribution. We introduce a further model extension, called the flexible local clock (FLC), which provides a flexible framework to combine relaxed clock models with local clock models. We evaluate the flexible local clock on simulated and real datasets and show that it provides substantially improved fit to an influenza dataset. An implementation of the model is available for download from https://www.github.com/4ment/flc.

Geographical and Historical Patterns in the Emergences of Novel Highly Pathogenic Avian Influenza (HPAI) H5 and H7 Viruses in Poultry

Over the years, the emergence of novel H5 and H7 highly pathogenic avian influenza viruses (HPAI) has been taking place through two main mechanisms: first, the conversion of a low pathogenic into a highly pathogenic virus, and second, the reassortment between different genetic segments of low and highly pathogenic viruses already in circulation. We investigated and summarized the literature on emerging HPAI H5 and H7 viruses with the aim of building a spatio-temporal database of all these recorded conversions and reassortments events. We subsequently mapped the spatio-temporal distribution of known emergence events, as well as the species and production systems that they were associated with, the aim being to establish their main characteristics. From 1959 onwards, we identified a total of 39 independent H7 and H5 LPAI to HPAI conversion events. All but two of these events were reported in commercial poultry production systems, and a majority of these events took place in high-income countries. In contrast, a total of 127 reassortments have been reported from 1983 to 2015, which predominantly took place in countries with poultry production systems transitioning from backyard to intensive production systems. Those systems are characterized by several co-circulating viruses, multiple host species, regular contact points in live bird markets, limited biosecurity within value chains, and frequent vaccination campaigns that impose selection pressures for emergence of novel reassortants. We conclude that novel HPAI emergences by these two mechanisms occur in different ecological niches, with different viral, environmental and host associated factors, which has implications in early detection and management and mitigation of the risk of emergence of novel HPAI viruses.

Changing Role of Wild Birds in the Epidemiology of Avian Influenza A Viruses

Waterbirds are the main reservoir for low pathogenic avian influenza A viruses (LPAIV), from which occasional spillover to poultry occurs. When circulating among poultry, LPAIV may become highly pathogenic avian influenza A viruses (HPAIV). In recent years, the epidemiology of HPAIV viruses has changed drastically. HPAIV H5N1 are currently endemic among poultry in a number of countries. In addition, global spread of HPAIV H5Nx viruses has resulted in major outbreaks among wild birds and poultry worldwide. Using data collected during these outbreaks, the role of migratory birds as a vector became increasingly clear. Here we provide an overview of current data about various aspects of the changing role of wild birds in the epidemiology of avian influenza A viruses.

Genetic variation and co-evolutionary relationship of RNA polymerase complex segments in influenza A viruses

The RNA polymerase complex (RNApc) in influenza A viruses (IVs) is composed of the PB2, PB1 and PA subunits, which are encoded by the three longest genome segments (Seg1-3) and are responsible for the replication of vRNAs and transcription of viral mRNAs. However, the co-evolutionary relationships of the three segments from the known 126 subtypes IVs are unclear. In this study, we performed a detailed analysis based on a total number of 121,191 nucleotide sequences. Three segment sequences were aligned before the repeated, incomplete and mixed sequences were removed for homologous and phylogenetic analyses. Subsequently, the estimated substitution rates and TMRCAs (Times for Most Recent Common Ancestor) were calculated by 175 representative IVs. Tracing the cladistic distribution of three segments from these IVs, co-evolutionary patterns and trajectories could be inferred. The further correlation analysis of six internal protein coding segments reflect the RNApc segments have the closer correlation than others during continuous reassortments. This global approach facilitates the establishment of a fast antiviral strategy and monitoring of viral variation.

The impact of migratory flyways on the spread of avian influenza virus in North America

Background

Wild birds are the major reservoir hosts for influenza A viruses (AIVs) and have been implicated in the emergence of pandemic events in livestock and human populations. Understanding how AIVs spread within and across continents is therefore critical to the development of successful strategies to manage and reduce the impact of influenza outbreaks. In North America many bird species undergo seasonal migratory movements along a North-South axis, thereby providing opportunities for viruses to spread over long distances. However, the role played by such avian flyways in shaping the genetic structure of AIV populations remains uncertain.

Results

To assess the relative contribution of bird migration along flyways to the genetic structure of AIV we performed a large-scale phylogeographic study of viruses sampled in the USA and Canada, involving the analysis of 3805 to 4505 sequences from 36 to 38 geographic localities depending on the gene segment data set. To assist in this we developed a maximum likelihood-based genetic algorithm to explore a wide range of complex spatial models, depicting a more complete picture of the migration network than determined previously.

Conclusions

Based on phylogenies estimated from nucleotide sequence data sets, our results show that AIV migration rates are significantly higher within than between flyways, indicating that the migratory patterns of birds play a key role in viral dispersal. These findings provide valuable insights into the evolution, maintenance and transmission of AIVs, in turn allowing the development of improved programs for surveillance and risk assessment.

Large-Scale Phylogenomic Analysis Reveals the Complex Evolutionary History of Rabies Virus in Multiple Carnivore Hosts

The natural evolution of rabies virus (RABV) provides a potent example of multiple host shifts and an important opportunity to determine the mechanisms that underpin viral emergence. Using 321 genome sequences spanning an unprecedented diversity of RABV, we compared evolutionary rates and selection pressures in viruses sampled from multiple primary host shifts that occurred on various continents. Two major phylogenetic groups, bat-related RABV and dog-related RABV, experiencing markedly different evolutionary dynamics were identified. While no correlation between time and genetic divergence was found in bat-related RABV, the evolution of dog-related RABV followed a generally clock-like structure, although with a relatively low evolutionary rate. Subsequent molecular clock dating indicated that dog-related RABV likely underwent a rapid global spread following the intensification of intercontinental trade starting in the 15^th century. Strikingly, although dog RABV has jumped to various wildlife species from the order Carnivora, we found no clear evidence that these host-jumping events involved adaptive evolution, with RABV instead characterized by strong purifying selection, suggesting that ecological processes also play an important role in shaping patterns of emergence. However, specific amino acid changes were associated with the parallel emergence of RABV in ferret-badgers in Asia, and some host shifts were associated with increases in evolutionary rate, particularly in the ferret-badger and mongoose, implying that changes in host species can have important impacts on evolutionary dynamics.

Zoonoses account for most recently emerged infectious diseases of humans, although little is known about the evolutionary mechanisms involved in cross-species virus transmission. Understanding the evolutionary patterns and processes that underpin such cross-species transmission is of importance for predicting the spread of zoonotic infections, and hence to their ultimate control. We present a large-scale and detailed reconstruction of the evolutionary history of rabies virus (RABV) in domestic and wildlife animal species. RABV is of particular interest as it is capable of infecting many mammals but, paradoxically, is only maintained in distinct epidemiological cycles associated with animal species from the orders Carnivora and Chiroptera. We show that bat-related RABV and dog-related RABV have experienced very different evolutionary dynamics, and that host jumps are sometimes characterized by significant increases in evolutionary rate. Among Carnivora, the association between RABV and particular host species most likely arose from a combination of the historical human-mediated spread of the virus and jumps into new primary host species. In addition, we show that changes in host species are associated with multiple evolutionary pathways including the occurrence of host-specific parallel evolution. Overall, our data indicate that the establishment of dog-related RABV in new carnivore hosts may only require subtle adaptive evolution.