Intrasubtype reassortments cause adaptive amino acid replacements in H3N2 influenza genes

A novel reassorted swine H3N2 influenza virus demonstrates an undetected human-to-swine spillover in Latin America and highlights zoonotic risks

Influenza A virus (FLUAV) affects a wide range of hosts, including humans and animals, posing a threat to public health. In swine, H3N2 subtype is associated with human-to-swine spillovers of seasonal viruses. In Latin America, the molecular and antigenic characteristics of swine FLUAV H3N2, as well as its phylogenetic origin, are poorly understood. Therefore, the objective of this study was to characterize the first swine H3N2 detected in Colombia. The origin and lineage of the virus were estimated through phylogenetic and molecular clock analyses. Antigenic characterization was achieved by comparing the amino acid constitution of the HA with previously reported swine FLUAVs and seasonal vaccine strains using a sequence-based method. In addition to HA and NA, internal genes were also characterized. The results showed that the Colombian H3N2 corresponded to a novel phylogenetic and antigenic swine FLUAV variant that emerged due to an independent reverse zoonotic event, likely occurring in Colombia in the early 2000s. The immunodominant epitope in the virus was predominantly present in antigenic epitope A, which showed the highest amino acid variation. Some mutations that alter the N-Glycosylation of antigenic sites at the HA were detected. Internally, the virus exhibited pandemic configuration. This study provides the first evidence of a novel FLUAV in Colombia and describes its origin, variability, and persistence in geographically restricted populations, highlighting the need for strengthen molecular surveillance of the virus in animal populations.

Pandemic risk stemming from the bovine H5N1 outbreak: an account of the knowns and unknowns

H5N1 subtype influenza A viruses represent a long-standing pandemic concern. Owing to their global occurrence in poultry, humans are routinely exposed to these viruses, and hundreds of human cases have been documented worldwide since 2003. The relevant viral lineages are not static, however, and have recently undergone a massive expansion of host range and geographic distribution. Within this expansion, the introduction of H5N1 viruses into dairy cattle in the United States has spawned a novel animal-human interface. In response, public health agencies have sought to evaluate the risk of an H5N1 pandemic stemming from the bovine outbreak. These assessments draw on evidence from the field and the laboratory to score a series of recognized risk factors. As such, their utility hinges on fundamental understanding of the processes that drive pandemic emergence and the availability of relevant data. Advancing this understanding and gathering data prior to and during an outbreak are primary missions of the NIAID Centers of Excellence for Influenza Research and Response (CEIRR) Network. To further these goals and highlight the need for an invigorated response across US agencies, here, we review gaps in understanding of the dairy cattle outbreak and identify constraints on efforts to close these gaps.

Multiple transatlantic incursions of highly pathogenic avian influenza clade 2.3.4.4b A(H5N5) virus into North America and spillover to mammals

Highly pathogenic avian influenza (HPAI) viruses have spread at an unprecedented scale, leading to mass mortalities in birds and mammals. In 2023, a transatlantic incursion of HPAI A(H5N5) viruses into North America was detected, followed shortly thereafter by a mammalian detection. As these A(H5N5) viruses were similar to contemporary viruses described in Eurasia, the transatlantic spread of A(H5N5) viruses was most likely facilitated by pelagic seabirds. Some of the Canadian A(H5N5) viruses from birds and mammals possessed the PB2-E627K substitution known to facilitate adaptation to mammals. Ferrets inoculated with A(H5N5) viruses showed rapid, severe disease onset, with some evidence of direct contact transmission. However, these viruses have maintained receptor binding traits of avian influenza viruses and were susceptible to oseltamivir and zanamivir. Understanding the factors influencing the virulence and transmission of A(H5N5) in migratory birds and mammals is critical to minimize impacts on wildlife and public health.

Molecular characterization of the whole genome of H9N2 avian influenza virus isolated from Egyptian poultry farms

H9N2 avian influenza viruses (AIVs) affect both poultry and humans on a global level, and they are especially prevalent in Egypt. In this study, we sequenced the entire genome of AIV H9N2 isolated from chickens in Egypt in 2021, using next-generation sequencing (NGS) technology. Phylogenetic analysis of the resulting sequences showed that the studied strain was generally monophyletic and grouped within the G1 sublineage of the Eurasian lineage. Four segments (polymerase basic 2 [PB2], polymerase basic 1 [PB1], polymerase acidic [PA], and non-structural [NS]) were related to Egyptian genotype II, while the nucleoprotein (NP), neuraminidase (NA), matrix (M), and haemagglutinin (HA) segments were related to Egyptian genotype I. Molecular analysis revealed that HA protein contained amino acid residues (191H and 234L) that suggested a predilection for attaching to human-like receptors. The antigenic sites of HA had two nonsynonymous mutations: V194I at antigenic site A and M40K at antigenic site B. Furthermore, the R403W and S372A mutations, which have been observed in H3N2 and H2N2 strains that caused human pandemics, were found in the NA protein of the detected strain. The internal proteins contained virulence markers: 504V in the PB2 protein, 622G, 436Y, 207K, and 677T in the PB1 protein, 127V, 550L, and 672L in PA protein, and 64F and 69P in the M protein. These results show that the detected strain had undergone intrasubtype reassortment. Furthermore, it contains changes in the viral proteins that make it more likely to be virulent, raising a question about the tendency of AIV H9N2 to become highly pathogenic in the future for both poultry and humans.

Characterizing a century of genetic diversity and contemporary antigenic diversity of N1 neuraminidase in influenza A virus from North American swine

Influenza A viruses (IAVs) of the H1N1 classical swine lineage became endemic in North American swine following the 1918 pandemic. Additional human-to-swine transmission events after 1918, and a spillover of H1 viruses from wild birds in Europe, potentiated a rapid increase in genomic diversity via reassortment between introductions and the endemic classical swine lineage. To determine mechanisms affecting reassortment and evolution, we conducted a phylogenetic analysis of N1 and paired HA swine IAV genes in North America between 1930 and 2020. We described fourteen N1 clades within the N1 Eurasian avian lineage (including the N1 pandemic clade), the N1 classical swine lineage, and the N1 human seasonal lineage. Seven N1 genetic clades had evidence for contemporary circulation. To assess antigenic drift associated with N1 genetic diversity, we generated a panel of representative swine N1 antisera and quantified the antigenic distance between wild-type viruses using enzyme-linked lectin assays and antigenic cartography. Within the N1 genes, the antigenic similarity was variable and reflected shared evolutionary history. Sustained circulation and evolution of N1 genes in swine had resulted in a significant antigenic distance between the N1 pandemic clade and the classical swine lineage. Between 2010 and 2020, N1 clades and N1-HA pairings fluctuated in detection frequency across North America, with hotspots of diversity generally appearing and disappearing within 2 years. We also identified frequent N1-HA reassortment events (n = 36), which were rarely sustained (n = 6) and sometimes also concomitant with the emergence of new N1 genetic clades (n = 3). These data form a baseline from which we can identify N1 clades that expand in range or genetic diversity that may impact viral phenotypes or vaccine immunity and subsequently the health of North American swine.

TreeKnit: Inferring ancestral reassortment graphs of influenza viruses

When two influenza viruses co-infect the same cell, they can exchange genome segments in a process known as reassortment. Reassortment is an important source of genetic diversity and is known to have been involved in the emergence of most pandemic influenza strains. However, because of the difficulty in identifying reassortment events from viral sequence data, little is known about their role in the evolution of the seasonal influenza viruses. Here we introduce TreeKnit, a method that infers ancestral reassortment graphs (ARG) from two segment trees. It is based on topological differences between trees, and proceeds in a greedy fashion by finding regions that are compatible in the two trees. Using simulated genealogies with reassortments, we show that TreeKnit performs well in a wide range of settings and that it is as accurate as a more principled bayesian method, while being orders of magnitude faster. Finally, we show that it is possible to use the inferred ARG to better resolve segment trees and to construct more informative visualizations of reassortments.

Influenza viruses evolve quickly and escape immune defenses which requires frequent update of vaccines. Understanding this evolution is key to an effective public health response. The genome of influenza viruses is made up of 8 pieces called segments, each coding for different viral proteins. Within each segment, evolution is an asexual process in which genetic diversity is generated by mutations. But influenza also diversifies through reassortment which can occur when two different viruses infect the same cell: offsprings can then contain a combination of segments from both viruses. Reassortment is akin to sexual reproduction and can generate viruses that combine segments from diverged viral lineages. Reassortment is a crucial component of viral evolution, but it is challenging to reconstruct where reassortments happened and which segments share history. Here, we develop a method called TreeKnit to detect reassortment events. TreeKnit is based on genealogical trees of single segments that can be reconstructed using standard bioinformatics tools. Inconsistencies between these trees are then used as signs of reassortment. We show that TreeKnit is as accurate as other recent methods, but runs much faster. Our method will facilitate the study of reassortment and its consequences for influenza evolution.

Genome-Wide Reassortment Analysis of Influenza A H7N9 Viruses Circulating in China during 2013-2019

Reassortment with the H9N2 virus gave rise to the zoonotic H7N9 avian influenza virus (AIV), which caused more than five outbreak waves in humans, with high mortality. The frequent exchange of genomic segments between H7N9 and H9N2 has been well-documented. However, the reassortment patterns have not been described and are not yet fully understood. Here, we used phylogenetic analyses to investigate the patterns of intersubtype and intrasubtype/intralineage reassortment across the eight viral segments. The H7N9 virus and its progeny frequently exchanged internal genes with the H9N2 virus but rarely with the other AIV subtypes. Before beginning the intrasubtype/intralineage reassortment analyses, five Yangtze River Delta (YRD A-E) and two Pearl River Delta (PRD A-B) clusters were divided according to the HA gene phylogeny. The seven reset segment genes were also nomenclatured consistently. As revealed by the tanglegram results, high intralineage reassortment rates were determined in waves 2–3 and 5. Additionally, the clusters of PB2 c05 and M c02 were the most dominant in wave 5, which could have contributed to the onset of the largest H7N9 outbreak in 2016–2017. Meanwhile, a portion of the YRD-C cluster (HP H7N9) inherited their PB2, PA, and M segments from the co-circulating YRD-E (LP H7N9) cluster during wave 5. Untanglegram results revealed that the reassortment rate between HA and NA was lower than HA with any of the other six segments. A multidimensional scaling plot revealed a robust genetic linkage between the PB2 and PA genes, indicating that they may share a co-evolutionary history. Furthermore, we observed relatively more robust positive selection pressure on HA, NA, M2, and NS1 proteins. Our findings demonstrate that frequent reassortment, particular reassorted patterns, and adaptive mutations shaped the H7N9 viral genetic diversity and evolution. Increased surveillance is required immediately to better understand the current state of the HP H7N9 AIV.

RF-Net 2: fast inference of virus reassortment and hybridization networks

MOTIVATION

A phylogenetic network is a powerful model to represent entangled evolutionary histories with both divergent (speciation) and convergent (e.g. hybridization, reassortment, recombination) evolution. The standard approach to inference of hybridization networks is to (i) reconstruct rooted gene trees and (ii) leverage gene tree discordance for network inference. Recently, we introduced a method called RF-Net for accurate inference of virus reassortment and hybridization networks from input gene trees in the presence of errors commonly found in phylogenetic trees. While RF-Net demonstrated the ability to accurately infer networks with up to four reticulations from erroneous input gene trees, its application was limited by the number of reticulations it could handle in a reasonable amount of time. This limitation is particularly restrictive in the inference of the evolutionary history of segmented RNA viruses such as influenza A virus (IAV), where reassortment is one of the major mechanisms shaping the evolution of these pathogens.

RESULTS

Here, we expand the functionality of RF-Net that makes it significantly more applicable in practice. Crucially, we introduce a fast extension to RF-Net, called Fast-RF-Net, that can handle large numbers of reticulations without sacrificing accuracy. In addition, we develop automatic stopping criteria to select the appropriate number of reticulations heuristically and implement a feature for RF-Net to output error-corrected input gene trees. We then conduct a comprehensive study of the original method and its novel extensions and confirm their efficacy in practice using extensive simulation and empirical IAV evolutionary analyses.

AVAILABILITY AND IMPLEMENTATION

RF-Net 2 is available at https://github.com/flu-crew/rf-net-2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Whole-genome analysis to determine the rate and patterns of intra-subtype reassortment among influenza type-A viruses in Africa

Influenza type-A viruses (IAVs) present a global burden of human respiratory infections and mortality. Genome reassortment is an important mechanism through which epidemiologically novel influenza viruses emerge and a core step in the safe reassortment-incompetent live-attenuated influenza vaccine development. Currently, there are no data on the rate, spatial and temporal distribution, and role of reassortment in the evolution and diversification of IAVs circulating in Africa. We aimed to detect intra-subtype reassortment among Africa pandemic H1N1pdm09 (2009-10), seasonal H1N1pdm09 (2011-20), and seasonal H3N2 viruses and characterize the genomic architecture and temporal and spatial distribution patterns of the resulting reassortants. Our study was nested within the Uganda National Influenza Surveillance Programme. Next-generation sequencing was used to generate whole genomes (WGs) from 234 H1N1pdm09 (n = 116) and H3N2 (n = 118) viruses sampled between 2010 and 2018 from seven districts in Uganda. We combined our newly generated WGs with 658 H1N1pdm09 and 1131 H3N2 WGs sampled between 1994 and 2020 across Africa and identified reassortants using an automated Graph Incompatibility Based Reassortment Finder software. Viral reassortment rates were estimated using a coalescent reassortant constant population model. Phylogenetic analysis was used to assess the effect of reassortment on viral genetic evolution. We observed a high frequency of intra-subtype reassortment events, 12 · 4 per cent (94/758) and 20 · 9 per cent (256/1,224), and reassortants, 13 · 3 per cent (101/758) and 38 · 6 per cent (472/1,224), among Africa H1N1pdm09 and H3N2 viruses, respectively. H1N1pdm09 reassorted at higher rates (0.1237-0.4255) than H3N2 viruses (0 · 00912-0.0355 events/lineage/year), a case unique to Uganda. Viral reassortants were sampled in 2009 through 2020, except in 2012. 78 · 2 per cent (79/101) of H1N1pdm09 reassortants acquired new non-structural, while 57 · 8 per cent (273/472) of the H3N2 reassortants had new hemagglutinin (H3) genes. Africa H3N2 viruses underwent more reassortment events involving larger reassortant sets than H1N1pdm09 viruses. Viruses with a specific reassortment architecture circulated for up to five consecutive years in specific countries and regions. The Eastern (Uganda and Kenya) and Western Africa harboured 84 · 2 per cent (85/101) and 55 · 9 per cent (264/472) of the continent's H1N1pdm09 and H3N2 reassortants, respectively. The frequent reassortment involving multi-genes observed among Africa IAVs showed the intracontinental viral evolution and diversification possibly sustained by viral importation from outside Africa and/or local viral genomic mixing and transmission. Novel reassortant viruses emerged every year, and some persisted in different countries and regions, thereby presenting a risk of influenza outbreaks in Africa. Our findings highlight Africa as part of the global influenza ecology and the advantage of implementing routine whole-over partial genome sequencing and analyses to monitor circulating and detect emerging viruses. Furthermore, this study provides evidence and heightens our knowledge on IAV evolution, which is integral in directing vaccine strain selection and the update of master donor viruses used in recombinant vaccine development.

Spatial and Temporal Coevolution of N2 Neuraminidase and H1 and H3 Hemagglutinin Genes of Influenza A Virus in United States Swine

The neuraminidase (NA) and hemagglutinin (HA) are essential surface glycoproteins of influenza A virus (IAV). In this study, the evolution of subtype N2 NA paired with H1 and H3 subtype HA in swine was evaluated to understand if genetic diversity of HA and NA were linked. Using time-scaled Bayesian phylodynamic analyses, the relationships of paired swine N2 with H1 or H3 from 2009 to 2018 were evaluated. These data demonstrated increased relative genetic diversity within the major N2 clades circulating in swine in the United States (N2.1998 between 2014-2017 and N2.2002 between 2010-2016). Preferential pairing was observed among specific NA and HA genetic clades. Gene reassortment between cocirculating influenza A strains resulted in novel pairings that persisted. The changes of genetic diversity in the NA gene were quantified using Bayesian phylodynamic analyses and increases in diversity were observed subsequent to novel NA-HA reassortment events. The rate of evolution among NA-N2 clades and HA-H1 and HA-H3 clades were similar. Bayesian phylodynamic analyses demonstrated strong spatial patterns in N2 genetic diversity, but frequent interstate movement of rare N2 clades provided opportunity for reassortment and emergence of new N2-HA pairings. The frequent regional movement of pigs and their influenza viruses is an explanation for the documented patterns of reassortment and subsequent changes in gene diversity. The reassortment and evolution of NA and linked HA evolution may result in antigenic drift of both major surface glycoproteins, reducing vaccine efficacy, with subsequent impact on animal health.

Whole-genome-based phylogenomic analysis of the Belgian 2016-2017 influenza A(H3N2) outbreak season allows improved surveillance

Seasonal influenza epidemics are associated with high mortality and morbidity in the human population. Influenza surveillance is critical for providing information to national influenza programmes and for making vaccine composition predictions. Vaccination prevents viral infections, but rapid influenza evolution results in emerging mutants that differ antigenically from vaccine strains. Current influenza surveillance relies on Sanger sequencing of the haemagglutinin (HA) gene. Its classification according to World Health Organization (WHO) and European Centre for Disease Prevention and Control (ECDC) guidelines is based on combining certain genotypic amino acid mutations and phylogenetic analysis. Next-generation sequencing technologies enable a shift to whole-genome sequencing (WGS) for influenza surveillance, but this requires laboratory workflow adaptations and advanced bioinformatics workflows. In this study, 253 influenza A(H3N2) positive clinical specimens from the 2016–2017 Belgian season underwent WGS using the Illumina MiSeq system. HA-based classification according to WHO/ECDC guidelines did not allow classification of all samples. A new approach, considering the whole genome, was investigated based on using powerful phylogenomic tools including beast and Nextstrain, which substantially improved phylogenetic classification. Moreover, Bayesian inference via beast facilitated reassortment detection by both manual inspection and computational methods, detecting intra-subtype reassortants at an estimated rate of 15 %. Real-time analysis (i.e. as an outbreak is ongoing) via Nextstrain allowed positioning of the Belgian isolates into the globally circulating context. Finally, integration of patient data with phylogenetic groups and reassortment status allowed detection of several associations that would have been missed when solely considering HA, such as hospitalized patients being more likely to be infected with A(H3N2) reassortants, and the possibility to link several phylogenetic groups to disease severity indicators could be relevant for epidemiological monitoring. Our study demonstrates that WGS offers multiple advantages for influenza monitoring in (inter)national influenza surveillance, and proposes an improved methodology. This allows leveraging all information contained in influenza genomes, and allows for more accurate genetic characterization and reassortment detection.

Diversity and Reassortment Rate of Influenza A Viruses in Wild Ducks and Gulls

Influenza A viruses (IAVs) evolve via point mutations and reassortment of viral gene segments. The patterns of reassortment in different host species differ considerably. We investigated the genetic diversity of IAVs in wild ducks and compared it with the viral diversity in gulls. The complete genomes of 38 IAVs of H1N1, H1N2, H3N1, H3N2, H3N6, H3N8, H4N6, H5N3, H6N2, H11N6, and H11N9 subtypes isolated from wild mallard ducks and gulls resting in a city pond in Moscow, Russia were sequenced. The analysis of phylogenetic trees showed that stable viral genotypes do not persist from year to year in ducks owing to frequent gene reassortment. For comparison, similar analyses were carried out using sequences of IAVs isolated in the same period from ducks and gulls in The Netherlands. Our results revealed a significant difference in diversity and rates of reassortment of IAVs in ducks and gulls.

A brief outline of respiratory viral disease outbreaks: 1889-till date on the public health perspectives

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is currently causing the respiratory illness termed as the coronavirus disease 2019 or the COVID-19 pandemic. Indeed, the significant increase in deaths in the current days due to influenza around the world started in 1889 is a continued public health threat because of its intermittent style of pandemic outbreaks. An array of research on the influenza viruses has been conducted especially pointing on (1) the development of the anti-viral drugs and the design of probable vaccines on trial basis, (2) the biochemical and genetic aspects underlying the viral pathogenicity, (3) the viral epidemiology, and on (4) the protective immunity against the influenza viruses. Current review briefly discussed the epidemic/ pandemic history of influenza and correlated with the current epidemiology, the possible preventive measures that may be taken by the public health professionals as well as to increase the protective awareness among the general people. The viral reassortments during the initiation of pandemics have also been focused based on the previous literatures.

Evolution and rapid spread of a reassortant A(H3N2) virus that predominated the 2017-2018 influenza season

The 2017-2018 North American influenza season caused more hospitalizations and deaths than any year since the 2009 H1N1 pandemic. The majority of recorded influenza infections were caused by A(H3N2) viruses, with most of the virus's North American diversity falling into the A2 clade. Within A2, we observe a subclade which we call A2/re that rose to comprise almost 70 per cent of A(H3N2) viruses circulating in North America by early 2018. Unlike most fast-growing clades, however, A2/re contains no amino acid substitutions in the hemagglutinin (HA) segment. Moreover, hemagglutination inhibition assays did not suggest substantial antigenic differences between A2/re viruses and viruses sampled during the 2016-2017 season. Rather, we observe that the A2/re clade was the result of a reassortment event that occurred in late 2016 or early 2017 and involved the combination of the HA and PB1 segments of an A2 virus with neuraminidase (NA) and other segments a virus from the clade A1b. The success of this clade shows the need for antigenic analysis that targets NA in addition to HA. Our results illustrate the potential for non-HA drivers of viral success and necessitate the need for more thorough tracking of full viral genomes to better understand the dynamics of influenza epidemics.

Allele-specific nonstationarity in evolution of influenza A virus surface proteins

Influenza A virus (IAV) is a major public health problem and a pandemic threat. Its evolution is largely driven by diversifying positive selection so that relative fitness of different amino acid variants changes with time due to changes in herd immunity or genomic context, and novel amino acid variants attain fitness advantage. Here, we hypothesize that diversifying selection also has another manifestation: the fitness associated with a particular amino acid variant should decline with time since its origin, as the herd immunity adapts to it. By tracing the evolution of antigenic sites at IAV surface proteins, we show that an amino acid variant becomes progressively more likely to become replaced by another variant with time since its origin-a phenomenon we call "senescence." Senescence is particularly pronounced at experimentally validated antigenic sites, implying that it is largely driven by host immunity. By contrast, at internal sites, existing variants become more favorable with time, probably due to arising contingent mutations at other epistatically interacting sites. Our findings reveal a previously undescribed facet of adaptive evolution and suggest approaches for prediction of evolutionary dynamics of pathogens.

Mutation and Epistasis in Influenza Virus Evolution

Influenza remains a persistent public health challenge, because the rapid evolution of influenza viruses has led to marginal vaccine efficacy, antiviral resistance, and the annual emergence of novel strains. This evolvability is driven, in part, by the virus's capacity to generate diversity through mutation and reassortment. Because many new traits require multiple mutations and mutations are frequently combined by reassortment, epistatic interactions between mutations play an important role in influenza virus evolution. While mutation and epistasis are fundamental to the adaptability of influenza viruses, they also constrain the evolutionary process in important ways. Here, we review recent work on mutational effects and epistasis in influenza viruses.

Predictive Modeling of Influenza Shows the Promise of Applied Evolutionary Biology

Seasonal influenza is controlled through vaccination campaigns. Evolution of influenza virus antigens means that vaccines must be updated to match novel strains, and vaccine effectiveness depends on the ability of scientists to predict nearly a year in advance which influenza variants will dominate in upcoming seasons. In this review, we highlight a promising new surveillance tool: predictive models. Based on data-sharing and close collaboration between the World Health Organization and academic scientists, these models use surveillance data to make quantitative predictions regarding influenza evolution. Predictive models demonstrate the potential of applied evolutionary biology to improve public health and disease control. We review the state of influenza predictive modeling and discuss next steps and recommendations to ensure that these models deliver upon their considerable biomedical promise.

Datamonkey 2.0: A Modern Web Application for Characterizing Selective and Other Evolutionary Processes

Inference of how evolutionary forces have shaped extant genetic diversity is a cornerstone of modern comparative sequence analysis. Advances in sequence generation and increased statistical sophistication of relevant methods now allow researchers to extract ever more evolutionary signal from the data, albeit at an increased computational cost. Here, we announce the release of Datamonkey 2.0, a completely re-engineered version of the Datamonkey web-server for analyzing evolutionary signatures in sequence data. For this endeavor, we leveraged recent developments in open-source libraries that facilitate interactive, robust, and scalable web application development. Datamonkey 2.0 provides a carefully curated collection of methods for interrogating coding-sequence alignments for imprints of natural selection, packaged as a responsive (i.e. can be viewed on tablet and mobile devices), fully interactive, and API-enabled web application. To complement Datamonkey 2.0, we additionally release HyPhy Vision, an accompanying JavaScript application for visualizing analysis results. HyPhy Vision can also be used separately from Datamonkey 2.0 to visualize locally executed HyPhy analyses. Together, Datamonkey 2.0 and HyPhy Vision showcase how scientific software development can benefit from general-purpose open-source frameworks. Datamonkey 2.0 is freely and publicly available at http://www.datamonkey.org, and the underlying codebase is available from https://github.com/veg/datamonkey-js.

Integrating patient and whole-genome sequencing data to provide insights into the epidemiology of seasonal influenza A(H3N2) viruses

Genetic surveillance of seasonal influenza is largely focused on sequencing of the haemagglutinin gene. Consequently, our understanding of the contribution of the remaining seven gene segments to the evolution and epidemiological dynamics of seasonal influenza is relatively limited. The increased availability of next-generation sequencing technologies allows rapid and economic whole-genome sequencing (WGS) of influenza virus. Here, 150 influenza A(H3N2) positive clinical specimens with linked epidemiological data, from the 2014/15 season in Scotland, were sequenced directly using both Sanger sequencing of the HA1 region and WGS using the Illumina MiSeq platform. Sequences generated by the two methods were highly correlated, and WGS provided on average >90 % whole genome coverage. As reported in other European countries during 2014/15, all strains belonged to genetic group 3C, with subgroup 3C.2a predominating. Multiple inter-subgroup reassortants were identified, including three 3C.3 viruses descended from a single reassortment event, which had persisted in the population. Cases of severe acute respiratory illness were significantly clustered on phylogenies of multiple gene segments indicating potential genetic factors warranting further investigation. Severe cases were also more likely to be associated with reassortant viruses and to occur later in the season. These results suggest that WGS provides an opportunity to develop our understanding of the relationship between the influenza genome and disease severity and the epidemiological consequences of within-subtype reassortment. Therefore, increased levels of WGS, linked to clinical and epidemiological data, could improve influenza surveillance.

Implications of segment mismatch for influenza A virus evolution

Influenza A virus (IAV) is an RNA virus with a segmented genome. These viral properties allow for the rapid evolution of IAV under selective pressure, due to mutation occurring from error-prone replication and the exchange of gene segments within a co-infected cell, termed reassortment. Both mutation and reassortment give rise to genetic diversity, but constraints shape their impact on viral evolution: just as most mutations are deleterious, most reassortment events result in genetic incompatibilities. The phenomenon of segment mismatch encompasses both RNA- and protein-based incompatibilities between co-infecting viruses and results in the production of progeny viruses with fitness defects. Segment mismatch is an important determining factor of the outcomes of mixed IAV infections and has been addressed in multiple risk assessment studies undertaken to date. However, due to the complexity of genetic interactions among the eight viral gene segments, our understanding of segment mismatch and its underlying mechanisms remain incomplete. Here, we summarize current knowledge regarding segment mismatch and discuss the implications of this phenomenon for IAV reassortment and diversity.

Selection on non-antigenic gene segments of seasonal influenza A virus and its impact on adaptive evolution

Most studies on seasonal influenza A/H3N2 virus adaptation have focused on the main antigenic gene, hemagglutinin. However, there is increasing evidence that the genome-wide genetic background of novel antigenic variants can influence these variants’ emergence probabilities and impact their patterns of dominance in the population. This suggests that non-antigenic genes may be important in shaping the viral evolutionary dynamics. To better understand the role of selection on non-antigenic genes in the adaptive evolution of seasonal influenza viruses, we have developed a simple population genetic model that considers a virus with one antigenic and one non-antigenic gene segment. By simulating this model under different regimes of selection and reassortment, we find that the empirical patterns of lineage turnover for the antigenic and non-antigenic gene segments are best captured when there is both limited viral coinfection and selection operating on both gene segments. In contrast, under a scenario of only neutral evolution in the non-antigenic gene segment, we see persistence of multiple lineages for long periods of time in that segment, which is not compatible with observed molecular evolutionary patterns. Further, we find that reassortment, occurring in coinfected individuals, can increase the speed of viral adaptive evolution by primarily reducing selective interference and genetic linkage effects. Together, these findings suggest that, for influenza, with six internal or non-antigenic gene segments, the evolutionary dynamics of novel antigenic variants are likely to be influenced by the genome-wide genetic background as a result of linked selection among both beneficial and deleterious mutations.

Fitness cost of reassortment in human influenza

Reassortment, which is the exchange of genome sequence between viruses co-infecting a host cell, plays an important role in the evolution of segmented viruses. In the human influenza virus, reassortment happens most frequently between co-existing variants within the same lineage. This process breaks genetic linkage and fitness correlations between viral genome segments, but the resulting net effect on viral fitness has remained unclear. In this paper, we determine rate and average selective effect of reassortment processes in the human influenza lineage A/H3N2. For the surface proteins hemagglutinin and neuraminidase, reassortant variants with a mean distance of at least 3 nucleotides to their parent strains get established at a rate of about 10⁻² in units of the neutral point mutation rate. Our inference is based on a new method to map reassortment events from joint genealogies of multiple genome segments, which is tested by extensive simulations. We show that intra-lineage reassortment processes are, on average, under substantial negative selection that increases in strength with increasing sequence distance between the parent strains. The deleterious effects of reassortment manifest themselves in two ways: there are fewer reassortment events than expected from a null model of neutral reassortment, and reassortant strains have fewer descendants than their non-reassortant counterparts. Our results suggest that influenza evolves under ubiquitous epistasis across proteins, which produces fitness barriers against reassortment even between co-circulating strains within one lineage.

The genome of the human influenza virus consists of 8 disjoint RNA polymer segments. These segments can undergo reassortment: when two viruses co-infect a host cell, they can produce viral offspring with a new combination of segments. In this paper, we show that reassortment within a given influenza lineage induces a fitness cost that increases in strength with increasing genetic distance of the parent viruses. Our finding suggests that evolution continuously produces viral proteins whose fitness depends on each other; reassortment reduces fitness by breaking up successful combinations of proteins. Thus, selection across proteins constrains viral evolution within a given lineage, and it may be an important factor in defining a viral species.

Adaptive evolution during the establishment of European avian-like H1N1 influenza A virus in swine

An H1N1 subtype influenza A virus with all eight gene segments derived from wild birds (including mallards), ducks and chickens, caused severe disease outbreaks in swine populations in Europe beginning in 1979 and successfully adapted to form the European avian‐like swine (EA‐swine) influenza lineage. Genes of the EA‐swine lineage that are clearly segregated from its closest avian relatives continue to circulate in swine populations globally and represent a unique opportunity to study the adaptive process of an avian‐to‐mammalian cross‐species transmission. Here, we used a relaxed molecular clock model to test whether the EA‐swine virus originated through the introduction of a single avian ancestor as an entire genome, followed by an analysis of host‐specific selection pressures among different gene segments. Our data indicated independent introduction of gene segments via transmission of avian viruses into swine followed by reassortment events that occurred at least 1–4 years prior to the EA‐swine outbreak. All EA‐swine gene segments exhibit greater selection pressure than avian viruses, reflecting both adaptive pressures and relaxed selective constraints that are associated with host switching. Notably, we identified key amino acid mutations in the viral surface proteins (H1 and N1) that play a role in adaptation to new hosts. Following the establishment of EA‐swine lineage, we observed an increased frequency of intrasubtype reassortment of segments compared to the earlier strains that has been associated with adaptive amino acid replacements, disease severity and vaccine escape. Taken together, our study provides key insights into the adaptive changes in viral genomes following the transmission of avian influenza viruses to swine and the early establishment of the EA‐swine lineage.

Seasonal H3N2 and 2009 Pandemic H1N1 Influenza A Viruses Reassort Efficiently but Produce Attenuated Progeny

Reassortment of gene segments between coinfecting influenza A viruses (IAVs) facilitates viral diversification and has a significant epidemiological impact on seasonal and pandemic influenza. Since 1977, human IAVs of H1N1 and H3N2 subtypes have cocirculated with relatively few documented cases of reassortment. We evaluated the potential for viruses of the 2009 pandemic H1N1 (pH1N1) and seasonal H3N2 lineages to reassort under experimental conditions. Results of heterologous coinfections with pH1N1 and H3N2 viruses were compared to those obtained following coinfection with homologous, genetically tagged, pH1N1 viruses as a control. High genotype diversity was observed among progeny of both coinfections; however, diversity was more limited following heterologous coinfection. Pairwise analysis of genotype patterns revealed that homologous reassortment was random while heterologous reassortment was characterized by specific biases. pH1N1/H3N2 reassortant genotypes produced under single-cycle coinfection conditions showed a strong preference for homologous PB2-PA combinations and general preferences for the H3N2 NA, pH1N1 M, and H3N2 PB2 except when paired with the pH1N1 PA or NP. Multicycle coinfection results corroborated these findings and revealed an additional preference for the H3N2 HA. Segment compatibility was further investigated by measuring chimeric polymerase activity and growth of selected reassortants in human tracheobronchial epithelial cells. In guinea pigs inoculated with a mixture of viruses, parental H3N2 viruses dominated but reassortants also infected and transmitted to cage mates. Taken together, our results indicate that strong intrinsic barriers to reassortment between seasonal H3N2 and pH1N1 viruses are few but that the reassortants formed are attenuated relative to parental strains.IMPORTANCE The genome of IAV is relatively simple, comprising eight RNA segments, each of which typically encodes one or two proteins. Each viral protein carries out multiple functions in coordination with other viral components and the machinery of the cell. When two IAVs coinfect a cell, they can exchange genes through reassortment. The resultant progeny viruses often suffer fitness defects due to suboptimal interactions among divergent viral components. The genetic diversity generated through reassortment can facilitate the emergence of novel outbreak strains. Thus, it is important to understand the efficiency of reassortment and the factors that limit its potential. The research described here offers new tools for studying reassortment between two strains of interest and applies those tools to viruses of the 2009 pandemic H1N1 and seasonal H3N2 lineages, which currently cocirculate in humans and therefore have the potential to give rise to novel epidemic strains.

M Gene Reassortment in H9N2 Influenza Virus Promotes Early Infection and Replication: Contribution to Rising Virus Prevalence in Chickens in China

Segment reassortment and base mutagenesis of influenza A viruses are the primary routes to the rapid evolution of high-fitness virus genotypes. We recently described a predominant G57 genotype of avian H9N2 viruses that caused countrywide outbreaks in chickens in China during 2010 to 2013, which led to the zoonotic emergence of H7N9 viruses. One of the key features of the G57 genotype is the replacement of the earlier A/chicken/Beijing/1/1994 (BJ/94)-like M gene with the A/quail/Hong Kong/G1/1997 (G1)-like M gene of quail origin. We report here the functional significance of the G1-like M gene in H9N2 viruses in conferring increased infection severity and infectivity in primary chicken embryonic fibroblasts and chickens. H9N2 virus housing the G1-like M gene, in place of the BJ/94-like M gene, showed an early surge in viral mRNA and viral RNA (vRNA) transcription that was associated with enhanced viral protein production and with an early elevated release of progeny virus comprising largely spherical rather than filamentous virions. Importantly, H9N2 virus with the G1-like M gene conferred extrapulmonary virus spread in chickens. Five highly represented signature amino acid residues (37A, 95K, 224N, and 242N in the M1 protein and 21G in the M2 protein) encoded by the prevalent G1-like M gene were demonstrated to be prime contributors to enhanced infectivity. Therefore, the genetic evolution of the M gene in H9N2 virus increases reproductive virus fitness, indicating its contribution to the rising virus prevalence in chickens in China.IMPORTANCE We recently described the circulation of a dominant genotype (genotype G57) of H9N2 viruses in countrywide outbreaks in chickens in China, which was responsible, through reassortment, for the emergence of H7N9 viruses that cause severe human infections. A key feature of the genotype G57 H9N2 virus is the presence of the quail-origin G1-like M gene, which had replaced the earlier BJ/94-like M gene. We found that H9N2 virus with the G1-like M gene, but not the BJ/94-like M gene, showed an early surge in progeny virus production and more severe pathology and extrapulmonary virus spread in chickens. Five highly represented amino acid residues in the M1 and M2 proteins derived from the G1-like M gene were shown to mediate enhanced virus infectivity. These observations enhance what we currently know about the roles of reassortment and mutations in virus fitness and have implications for assessing the potential of variant influenza viruses that can cause a rising prevalence in chickens.

Frequency of influenza H3N2 intra-subtype reassortment: attributes and implications of reassortant spread

BACKGROUND

Increasing evidence suggests that influenza reassortment not only contributes to the emergence of new human pandemics but also plays an important role in seasonal influenza epidemics, disease severity, evolution, and vaccine efficacy. We studied this process within 2091 H3N2 full genomes utilizing a combination of the latest reassortment detection tools and more conventional phylogenetic analyses.

RESULTS

We found that the amount of H3N2 intra-subtype reassortment depended on the number of sampled genomes, occurred with a steady frequency of 3.35%, and was not affected by the geographical origins, evolutionary patterns, or previous reassortment history of the virus. We identified both single reassortant genomes and reassortant clades, each clade representing one reassortment event followed by successful spread of the reassorted variant in the human population. It was this spread that was mainly responsible for the observed high presence of H3N2 intra-subtype reassortant genomes. The successfully spread variants were generally sampled within one year of their formation, highlighting the risk of their rapid spread but also presenting an opportunity for their rapid detection. Simultaneous spread of several different reassortant lineages was observed, and despite their limited average lifetime, second and third generation reassortment was detected, as well as reassortment between viruses belonging to different vaccine-associated clades, likely displaying differing antigenic properties. Some of the spreading reassortants remained confined to certain geographical regions, while others, sharing common properties in amino acid positions of the HA, NA, and PB2 segments, were found throughout the world.

CONCLUSIONS

Detailed surveillance of seasonal influenza reassortment patterns and variant properties may provide unique information needed for prediction of spread and construction of future influenza vaccines.

Changing preferences: deformation of single position amino acid fitness landscapes and evolution of proteins

The fitness landscape-the function that relates genotypes to fitness-and its role in directing evolution are a central object of evolutionary biology. However, its huge dimensionality precludes understanding of even the basic aspects of its shape. One way to approach it is to ask a simpler question: what are the properties of a function that assigns fitness to each possible variant at just one particular site-a single position fitness landscape-and how does it change in the course of evolution? Analyses of genomic data from multiple species and multiple individuals within a species have proved beyond reasonable doubt that fitness functions of positions throughout the genome do themselves change with time, thus shaping protein evolution. Here, I will briefly review the literature that addresses these dynamics, focusing on recent genome-scale analyses of fitness functions of amino acid sites, i.e. vectors of fitnesses of 20 individual amino acid variants at a given position of a protein. The set of amino acids that confer high fitness at a particular position changes with time, and the rate of this change is comparable with the rate at which a position evolves, implying that this process plays a major role in evolutionary dynamics. However, the causes of these changes remain largely unclear.

Genome-scale phylodynamics and evolution analysis of global H7N7 influenza viruses

Previous studies lacked of comprehensive analysis about the evolutionary history and phylogeography of global H7N7 viruses. In this study, it is essential to undertake a genome-scale analysis to investigate the evolutionary processes in a global perspective. There was local phylogenetic divergence among eight trees based on individual segments of 132 strains. We detected four reassortments between four distinct groups of viruses divided by HA gene, suggesting intrasubtype reassortment could accelerate the emergence of highly pathogenic virus. The molecular clock estimated that H7N7 virus evolved at a slower evolutionary rate ranged from 1.03E-03 to 2.81E-03subs/site/year. And we also showed that all gene segments of the virus were under strong purifying selection. A total of 11 positively selected sites were detected by at least two out of three methods. We reconstructed the population dynamics of global H7N7 viruses spanning over a century, revealing that temporal trends of the effective population size were consistent with the major epidemics previously reported. Our study adopt a Bayesian phylogeographic approach to investigate the geographic spread of H7N7 viruses, which combined with temporal and spatial information of all sequences. We have confirmed several migration events between different geographic locations supported by higher values of Bayes factor. The diffusion patterns of H7N7 viruses reveal that the virus is more likely to evolve to expand their host ranges even cross the species.

Reassortment in segmented RNA viruses: mechanisms and outcomes

Segmented RNA viruses are widespread in nature and include important human, animal and plant pathogens, such as influenza viruses and rotaviruses. Although the origin of RNA virus genome segmentation remains elusive, a major consequence of this genome structure is the capacity for reassortment to occur during co-infection, whereby segments are exchanged among different viral strains. Therefore, reassortment can create viral progeny that contain genes that are derived from more than one parent, potentially conferring important fitness advantages or disadvantages to the progeny virus. However, for segmented RNA viruses that package their multiple genome segments into a single virion particle, reassortment also requires genetic compatibility between parental strains, which occurs in the form of conserved packaging signals, and the maintenance of RNA and protein interactions. In this Review, we discuss recent studies that examined the mechanisms and outcomes of reassortment for three well-studied viral families - Cystoviridae, Orthomyxoviridae and Reoviridae - and discuss how these findings provide new perspectives on the replication and evolution of segmented RNA viruses.

Coordinated Evolution of Influenza A Surface Proteins

The surface proteins hemagglutinin (HA) and neuraminidase (NA) of human influenza A virus evolve under selection pressures to escape adaptive immune responses and antiviral drug treatments. In addition to these external selection pressures, some mutations in HA are known to affect the adaptive landscape of NA, and vice versa, because these two proteins are physiologically interlinked. However, the extent to which evolution of one protein affects the evolution of the other one is unknown. Here we develop a novel phylogenetic method for detecting the signatures of such genetic interactions between mutations in different genes – that is, inter-gene epistasis. Using this method, we show that influenza surface proteins evolve in a coordinated way, with mutations in HA affecting subsequent spread of mutations in NA and vice versa, at many sites. Of particular interest is our finding that the oseltamivir-resistance mutations in NA in subtype H1N1 were likely facilitated by prior mutations in HA. Our results illustrate that the adaptive landscape of a viral protein is remarkably sensitive to its genomic context and, more generally, that the evolution of any single protein must be understood within the context of the entire evolving genome.

The fitness of an organism depends on the coordinated function of many genes. Thus, how a mutation in one gene affects fitness often depends on what mutations are present in other genes. This dependence is called “genetic interaction” or “epistasis”. The prevalence and type of such interactions are not well understood. Epistasis can be inferred from time-series sequencing data when a mutation in one gene is observed to facilitate the spread of a mutation in another gene. However, the situation is much more complicated when new combinations of genes are formed by processes such as recombination or reassortment. In such cases, deducing the time and order of genetic changes is difficult. Here, we devise a method to infer pairs of mutations in different genes which closely follow one another in the presence of reassortment. We apply it to evolution of two surface proteins of influenza A virus, hemagglutinin and neuraminidase, which are important targets for the human immune system and drugs. We show that mutations in one of these proteins are often facilitated by prior mutations, or compensated by subsequent mutations, in the other protein. In particular, drug-resistance mutations in neuraminidase were likely made possible by prior mutation in hemagglutinin. Knowledge of such interactions is necessary to fully understand and predict evolution.

Positive Selection Drives Preferred Segment Combinations during Influenza Virus Reassortment

Influenza A virus (IAV) has a segmented genome that allows for the exchange of genome segments between different strains. This reassortment accelerates evolution by breaking linkage, helping IAV cross species barriers to potentially create highly virulent strains. Challenges associated with monitoring the process of reassortment in molecular detail have limited our understanding of its evolutionary implications. We applied a novel deep sequencing approach with quantitative analysis to assess the in vitro temporal evolution of genomic reassortment in IAV. The combination of H1N1 and H3N2 strains reproducibly generated a new H1N2 strain with the hemagglutinin and nucleoprotein segments originating from H1N1 and the remaining six segments from H3N2. By deep sequencing the entire viral genome, we monitored the evolution of reassortment, quantifying the relative abundance of all IAV genome segments from the two parent strains over time and measuring the selection coefficients of the reassorting segments. Additionally, we observed several mutations coemerging with reassortment that were not found during passaging of pure parental IAV strains. Our results demonstrate how reassortment of the segmented genome can accelerate viral evolution in IAV, potentially enabled by the emergence of a small number of individual mutations.

Haemagglutinin mutations and glycosylation changes shaped the 2012/13 influenza A(H3N2) epidemic, Houston, Texas

While the early start and higher intensity of the 2012/13 influenza A virus (IAV) epidemic was not unprecedented, it was the first IAV epidemic season since the 2009 H1N1 influenza pandemic where the H3N2 subtype predominated. We directly sequenced the genomes of 154 H3N2 clinical specimens collected throughout the epidemic to better understand the evolution of H3N2 strains and to inform the H3N2 vaccine selection process. Phylogenetic analyses indicated that multiple co-circulating clades and continual antigenic drift in the haemagglutinin (HA) of clades 5, 3A, and 3C, with the evolution of a new 3C subgroup (3C-2012/13), were the driving causes of the epidemic. Drift variants contained HA substitutions and alterations in the potential N-linked glycosylation sites of HA. Antigenic analysis demonstrated that viruses in the emerging subclade 3C.3 and subgroup 3C-2012/13 were not well inhibited by antisera generated against the 3C.1 vaccine strains used for the 2012/13 (A/Victoria/361/2011) or 2013/14 (A/Texas/50/2012) seasons. Our data support updating the H3N2 vaccine strain to a clade 3C.2 or 3C.3-like strain or a subclade that has drifted further. They also underscore the challenges in vaccine strain selection, particularly regarding HA and neuraminidase substitutions derived during laboratory passage that may alter antigenic testing accuracy.

Influenza virus non-structural protein NS1: interferon antagonism and beyond

Most viruses express one or several proteins that counter the antiviral defences of the host cell. This is the task of non-structural protein NS1 in influenza viruses. Absent in the viral particle, but highly expressed in the infected cell, NS1 dramatically inhibits cellular gene expression and prevents the activation of key players in the IFN system. In addition, NS1 selectively enhances the translation of viral mRNAs and may regulate the synthesis of viral RNAs. Our knowledge of the virus and of NS1 has increased dramatically during the last 15 years. The atomic structure of NS1 has been determined, many cellular partners have been identified and its multiple activities have been studied in depth. This review presents our current knowledge, and attempts to establish relationships between the RNA sequence, the structure of the protein, its ligands, its activities and the pathogenicity of the virus. A better understanding of NS1 could help in elaborating novel antiviral strategies, based on either live vaccines with altered NS1 or on small-compound inhibitors of NS1.

Genetic dynamic analysis of the influenza A H5N1 NS1 gene in China

The direct precursors of the A/Goose/Guangdong/1/1996 (GS/GD) virus lineage and its reassortants have been established geographically and ecologically. To investigate the variation and evolutionary dynamics of H5N1 viruses, whole-genome viral sequences (n = 164) were retrieved from the NCBI Influenza Virus Resource. Here, we present phylogenetic evidence for intrasubtype reassortments among H5N1 viruses isolated from China during 1996–2012. On the basis of phylogenetic analysis, we identified four major groups and further classified the reassortant viruses into three subgroups. Putative mosaic structures were mostly found in the viral ribonucleoprotein (vRNP) complexes and 91.0% (10/11) mosaics were obtained from terrestrial birds. Sequence variability and selection pressure analyses revealed that both surface glycoproteins (HA and NA) and nonstructural protein 1 (NS1) have higher dN/dS ratio and variability than other internal proteins. Furthermore, we detected 47 positively selected sites in genomic segments with the exception of PB2 and M1 genes. Hemagglutinin (HA) and neuraminidase (NA) are considered highly variable due to host immune pressure, however, it is not known what drives NS1 variability. Therefore, we performed a thorough analysis of the genetic variation and selective pressure of NS1 protein (462 available NS1 sequences). We found that most of positively selected sites and variable amino acids were located in the C-terminal effector domain (ED) of NS1. In addition, we focused on the NS1–RNA and NS1–protein interactions that were involved in viral replication mechanisms and host immune response. Transcriptomic analysis of H5N1-infected monkey lungs showed that certain PI3K-related genes were up-regulated.

Influenza A virus reassortment

Reassortment is the process by which influenza viruses swap gene segments. This genetic exchange is possible due to the segmented nature of the viral genome and occurs when two differing influenza viruses co-infect a cell. The viral diversity generated through reassortment is vast and plays an important role in the evolution of influenza viruses. Herein we review recent insights into the contribution of reassortment to the natural history and epidemiology of influenza A viruses, gained through population scale phylogenic analyses. We describe methods currently used to study reassortment in the laboratory, and we summarize recent progress made using these experimental approaches to further our understanding of influenza virus reassortment and the contexts in which it occurs.