Evolutionary dynamics of the N1 neuraminidases of the main lineages of influenza A viruses

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.

Molecular footprints of selective pressure in the neuraminidase gene of currently circulating human influenza subtypes and lineages

Influenza neuraminidase (NA) is under selective pressure (SP) of both host immune system and drug use. Here, we assembled large datasets of NA sequences of worldwide circulating viruses to estimate the global and site-specific SP acting on all current subtypes/lineages of human influenza NA. An overall negative SP of similar magnitude and a prevalence of negatively selected sites were observed for all subtypes/lineages. Positively selected sites varied according to the subtype/lineage, including N1-NA sites 247 and 275, N2-NA sites 148 and 151, and B/Victoria-NA site 395 associated with drug-resistance or reduced susceptibility. These results evidenced a potential role of positive selection in the low-level spread of A(H1N1)pdm09-H275Y drug-resistant viruses, and alerted for a potential higher risk of spread of a synergistic A(H1N1)pdm09 drug-resistant variant (H275Y/S247N). The positive selection detected at N2-NA sites 148 and 151 was probably an artefact from cell-culture. Overall mapping revealed six potential new druggable regions.

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.

Time dependency of foamy virus evolutionary rate estimates

Background

It appears that substitution rate estimates co-vary very strongly with their timescale of measurement; the shorter the timescale, the higher the estimated value. Foamy viruses have a long history of co-speciation with their hosts, and one of the lowest estimated rates of evolution among viruses. However, when their rate of evolution is estimated over short timescales, it is more reminiscent of the rapid rates seen in other RNA viruses. This discrepancy between their short-term and long-term rates could be explained by the time-dependency of substitution rate estimates. Several empirical models have been proposed and used to correct for the time-dependent rate phenomenon (TDRP), such as a vertically-translated exponential rate decay model and a power-law rate decay model. Nevertheless, at present, it is still unclear which model best describes the rate dynamics. Here, we use foamy viruses as a case study to empirically describe the phenomenon and to determine how to correct rate estimates for its effects. Four empirical models were investigated: (i) a vertically-translated exponential rate decay model, (ii) a simple exponential rate decay model, (iii) a vertically-translated power-law rate decay model, and (iv) a simple power-law rate decay model.

Results

Our results suggest that the TDRP is likely responsible for the large discrepancy observed in foamy virus short-term and long-term rate estimates, and the simple power-law rate decay model is the best model for inferring evolutionary timescales. Furthermore, we demonstrated that, within the Bayesian phylogenetic framework, currently available molecular clocks can severely bias evolutionary date estimates, indicating that they are inadequate for correcting for the TDRP. Our analyses also suggest that different viral lineages may have different TDRP dynamics, and this may bias date estimates if it is unaccounted for.

Conclusions

As evolutionary rate estimates are dependent on their measurement timescales, their values must be used and interpreted under the context of the timescale of rate estimation. Extrapolating rate estimates across large timescales for evolutionary inferences can severely bias the outcomes. Given that the TDRP is widespread in nature but has been noted only recently the estimated timescales of many viruses may need to be reconsidered and re-estimated. Our models could be used as a guideline to further improve current phylogenetic inference tools.

Electronic supplementary material

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

Intrasubtype reassortments cause adaptive amino acid replacements in H3N2 influenza genes

Reassortments and point mutations are two major contributors to diversity of Influenza A virus; however, the link between these two processes is unclear. It has been suggested that reassortments provoke a temporary increase in the rate of amino acid changes as the viral proteins adapt to new genetic environment, but this phenomenon has not been studied systematically. Here, we use a phylogenetic approach to infer the reassortment events between the 8 segments of influenza A H3N2 virus since its emergence in humans in 1968. We then study the amino acid replacements that occurred in genes encoded in each segment subsequent to reassortments. In five out of eight genes (NA, M1, HA, PB1 and NS1), the reassortment events led to a transient increase in the rate of amino acid replacements on the descendant phylogenetic branches. In NA and HA, the replacements following reassortments were enriched with parallel and/or reversing replacements; in contrast, the replacements at sites responsible for differences between antigenic clusters (in HA) and at sites under positive selection (in NA) were underrepresented among them. Post-reassortment adaptive walks contribute to adaptive evolution in Influenza A: in NA, an average reassortment event causes at least 2.1 amino acid replacements in a reassorted gene, with, on average, 0.43 amino acid replacements per evolving post-reassortment lineage; and at least ~9% of all amino acid replacements are provoked by reassortments.

Evolutionary history and phylodynamics of influenza A and B neuraminidase (NA) genes inferred from large-scale sequence analyses

Background

Influenza neuraminidase (NA) is an important surface glycoprotein and plays a vital role in viral replication and drug development. The NA is found in influenza A and B viruses, with nine subtypes classified in influenza A. The complete knowledge of influenza NA evolutionary history and phylodynamics, although critical for the prevention and control of influenza epidemics and pandemics, remains lacking.

Methodology/Principal findings

Evolutionary and phylogenetic analyses of influenza NA sequences using Maximum Likelihood and Bayesian MCMC methods demonstrated that the divergence of influenza viruses into types A and B occurred earlier than the divergence of influenza A NA subtypes. Twenty-three lineages were identified within influenza A, two lineages were classified within influenza B, and most lineages were specific to host, subtype or geographical location. Interestingly, evolutionary rates vary not only among lineages but also among branches within lineages. The estimated tMRCAs of influenza lineages suggest that the viruses of different lineages emerge several months or even years before their initial detection. The d_N/d_S ratios ranged from 0.062 to 0.313 for influenza A lineages, and 0.257 to 0.259 for influenza B lineages. Structural analyses revealed that all positively selected sites are at the surface of the NA protein, with a number of sites found to be important for host antibody and drug binding.

Conclusions/Significance

The divergence into influenza type A and B from a putative ancestral NA was followed by the divergence of type A into nine NA subtypes, of which 23 lineages subsequently diverged. This study provides a better understanding of influenza NA lineages and their evolutionary dynamics, which may facilitate early detection of newly emerging influenza viruses and thus improve influenza surveillance.

One misdated sequence of rabbit hemorrhagic disease virus prevents accurate estimation of its nucleotide substitution rate

BACKGROUND

The literature is ripe with phylogenetic estimates of nucleotide substitution rates, especially of measurably evolving species such as RNA viruses. However, it is not known how robust these rate estimates are to inaccuracies in the data, particularly in sampling dates that are used for molecular clock calibration. Here we report on the rate of evolution of the emerging pathogen Rabbit hemorrhagic disease virus (RHDV), which has significantly different rates of evolution for the same outer capsid (VP60) gene published in the literature. In an attempt to reconcile the conflicting data and further elucidate details of RHDV 's evolutionary history, we undertook fresh Bayesian analyses and employed jackknife control methods to produce robust substitution rate and time to most recent common ancestor (TMRCA) estimates for RHDV based on the VP60 and RNA-dependent RNA polymerase genes.

RESULTS

Through these control methods, we were able to identify a single misdated taxon, a passaged lab strain used for vaccine production, which was responsible for depressing the RHDV capsid gene's rate of evolution by 65%. Without this isolate, the polymerase and the capsid protein genes had nearly identical rates of evolution: 1.90x10-3 nucleotide substitutions/site/year, ns/s/y, (95% highest probability density (HPD) 1.25x10-3-2.55x10-3) and 1.91x10-3 ns/s/y (95% HPD 1.50x10-3-2.34x10-3), respectively.

CONCLUSIONS

After excluding the misdated taxon, both genes support a significantly higher substitution rate as well as a relatively recent emergence of RHDV, and obviate the need for previously hypothesized decades of unobserved diversification of the virus. The control methods show that using even one misdated taxon in a large dataset can significantly skew estimates of evolutionary parameters and suggest that it is better practice to use smaller datasets composed of taxa with unequivocal isolation dates. These jackknife controls would be useful for future tip-calibrated rate analyses that include taxa with ambiguous dates of isolation.

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.