Anomalies in the influenza virus genome database: new biology or laboratory errors?

Origin of the H1N1 (Russian influenza) pandemic of 1977-A risk assessment using the modified Grunow-Finke tool (mGFT)

In 1977, the Soviet Union (Union of Soviet Socialist Republics [USSR]) notified the World Health Organization (WHO) about an outbreak of H1N1 influenza, which later spread to many countries. The H1N1 strain of 1977 reappeared after being absent from the world for over 20 years. This pandemic simultaneously spread to several cities in the USSR and China. Many theories have been postulated to account for the emergence of this pandemic, including natural and unnatural origins. The purpose of this study was to use the modified Grunow-Finke risk assessment tool (modified Grunow-Finke tool [mGFT]) to investigate the origin of the 1977 H1N1 pandemic. Data was collected from WHO archives and published documents. The assessment of the pandemic's origin involved the utilization of a modified version of the original Grunow-Finke risk assessment tool (GFT). Using the mGFT, the final score was 37 out of 60 points (probability: 62%), indicating a high likelihood that the Russian influenza pandemic of 1977 was of unnatural origin. Several variables supported this finding, including the sudden re-emergence of a previously extinct strain, a genetic signature of laboratory modification for vaccine development, and unusual epidemiology. Inter-rater reliability was moderate to high. By applying the mGFT to the 1977 Russian influenza pandemic, we established a high probability that this pandemic was of unnatural origin. Although this is not definitive, it is consistent with the possibility that it originated from an incompletely attenuated live influenza vaccine. The mGFT is a useful risk analysis tool to evaluate the origin of epidemics.

A universal RT-qPCR assay for "One Health" detection of influenza A viruses

The mutual dependence of human and animal health is central to the One Health initiative as an integrated strategy for infectious disease control and management. A crucial element of the One Health includes preparation and response to influenza A virus (IAV) threats at the human-animal interface. The IAVs are characterized by extensive genetic variability, they circulate among different hosts and can establish host-specific lineages. The four main hosts are: avian, swine, human and equine, with occasional transmission to other mammalian species. The host diversity is mirrored in the range of the RT-qPCR assays for IAV detection. Different assays are recommended by the responsible health authorities for generic IAV detection in birds, swine or humans. In order to unify IAV monitoring in different hosts and apply the One Health approach, we developed a single RT-qPCR assay for universal detection of all IAVs of all subtypes, species origin and global distribution. The assay design was centred on a highly conserved region of the IAV matrix protein (MP)-segment identified by a comprehensive analysis of 99,353 sequences. The reaction parameters were effectively optimised with efficiency of 93–97% and LOD_95% of approximately ten IAV templates per reaction. The assay showed high repeatability, reproducibility and robustness. The extensive in silico evaluation demonstrated high inclusivity, i.e. perfect sequence match in the primers and probe binding regions, established as 94.6% for swine, 98.2% for avian and 100% for human H3N2, pandemic H1N1, as well as other IAV strains, resulting in an overall predicted detection rate of 99% on the analysed dataset. The theoretical predictions were confirmed and extensively validated by collaboration between six veterinary or human diagnostic laboratories on a total of 1970 specimens, of which 1455 were clinical and included a diverse panel of IAV strains.

In silico re-assessment of a diagnostic RT-qPCR assay for universal detection of Influenza A viruses

The ongoing evolution of microbial pathogens represents a significant issue in diagnostic PCR/qPCR. Many assays are burdened with false negativity due to mispriming and/or probe-binding failures. Therefore, PCR/qPCR assays used in the laboratory should be periodically re-assessed in silico on public sequences to evaluate the ability to detect actually circulating strains and to infer potentially escaping variants. In the work presented we re-assessed a RT-qPCR assay for the universal detection of influenza A (IA) viruses currently recommended by the European Union Reference Laboratory for Avian Influenza. To this end, the primers and probe sequences were challenged against more than 99,000 M-segment sequences in five data pools. To streamline this process, we developed a simple algorithm called the SequenceTracer designed for alignment stratification, compression, and personal sequence subset selection and also demonstrated its utility. The re-assessment confirmed the high inclusivity of the assay for the detection of avian, swine and human pandemic H1N1 IA viruses. On the other hand, the analysis identified human H3N2 strains with a critical probe-interfering mutation circulating since 2010, albeit with a significantly fluctuating proportion. Minor variations located in the forward and reverse primers identified in the avian and swine data were also considered.

Positive Selection in CD8+ T-Cell Epitopes of Influenza Virus Nucleoprotein Revealed by a Comparative Analysis of Human and Swine Viral Lineages

Numerous experimental studies have demonstrated that CD8(+) T cells contribute to immunity against influenza by limiting viral replication. It is therefore surprising that rigorous statistical tests have failed to find evidence of positive selection in the epitopes targeted by CD8(+) T cells. Here we use a novel computational approach to test for selection in CD8(+) T-cell epitopes. We define all epitopes in the nucleoprotein (NP) and matrix protein (M1) with experimentally identified human CD8(+) T-cell responses and then compare the evolution of these epitopes in parallel lineages of human and swine influenza viruses that have been diverging since roughly 1918. We find a significant enrichment of substitutions that alter human CD8(+) T-cell epitopes in NP of human versus swine influenza virus, consistent with the idea that these epitopes are under positive selection. Furthermore, we show that epitope-altering substitutions in human influenza virus NP are enriched on the trunk versus the branches of the phylogenetic tree, indicating that viruses that acquire these mutations have a selective advantage. However, even in human influenza virus NP, sites in T-cell epitopes evolve more slowly than do nonepitope sites, presumably because these epitopes are under stronger inherent functional constraint. Overall, our work demonstrates that there is clear selection from CD8(+) T cells in human influenza virus NP and illustrates how comparative analyses of viral lineages from different hosts can identify positive selection that is otherwise obscured by strong functional constraint.

IMPORTANCE

There is a strong interest in correlates of anti-influenza immunity that are protective against diverse virus strains. CD8(+) T cells provide such broad immunity, since they target conserved viral proteins. An important question is whether T-cell immunity is sufficiently strong to drive influenza virus evolution. Although many studies have shown that T cells limit viral replication in animal models and are associated with decreased symptoms in humans, no studies have proven with statistical significance that influenza virus evolves under positive selection to escape T cells. Here we use comparisons of human and swine influenza viruses to rigorously demonstrate that human influenza virus evolves under pressure to fix mutations in the nucleoprotein that promote escape from T cells. We further show that viruses with these mutations have a selective advantage since they are preferentially located on the "trunk" of the phylogenetic tree. Overall, our results show that CD8(+) T cells targeting nucleoprotein play an important role in shaping influenza virus evolution.

The Reemergent 1977 H1N1 Strain and the Gain-of-Function Debate

The 1977-1978 influenza epidemic was probably not a natural event, as the genetic sequence of the virus was nearly identical to the sequences of decades-old strains. While there are several hypotheses that could explain its origin, the possibility that the 1977 epidemic resulted from a laboratory accident has recently gained popularity in discussions about the biosafety risks of gain-of-function (GOF) influenza virus research, as an argument for why this research should not be performed. There is now a moratorium in the United States on funding GOF research while the benefits and risks, including the potential for accident, are analyzed. Given the importance of this historical epidemic to ongoing policy debates, we revisit the evidence that the 1977 epidemic was not natural and examine three potential origins: a laboratory accident, a live-vaccine trial escape, or deliberate release as a biological weapon. Based on available evidence, the 1977 strain was indeed too closely matched to decades-old strains to likely be a natural occurrence. While the origin of the outbreak cannot be conclusively determined without additional evidence, there are very plausible alternatives to the laboratory accident hypothesis, diminishing the relevance of the 1977 experience to the modern GOF debate.

H7N9 influenza A virus in turkeys in Minnesota

Introductions of H7 influenza A virus (IAV) from wild birds into poultry have been documented worldwide, resulting in varying degrees of morbidity and mortality. H7 IAV infection in domestic poultry has served as a source of human infection and disease. We report the detection of H7N9 subtype IAVs in Minnesota (MN) turkey farms during 2009 and 2011. The full genome was sequenced from eight isolates as well as the haemagglutinin (HA) and neuraminidase (NA) gene segments of H7 and N9 virus subtypes for 108 isolates from North American wild birds between 1986 and 2012. Through maximum-likelihood and coalescent phylogenetic analyses, we identified the recent H7 and N9 IAV ancestors of the turkey-origin H7N9 IAVs, estimated the time and geographical origin of the ancestral viruses, and determined the relatedness between the 2009 and 2011 turkey-origin H7N9 IAVs. Analyses supported that the 2009 and 2011 viruses were distantly related genetically, suggesting that the two outbreaks arose from independent introduction events from wild birds. Our findings further supported that the 2011 MN turkey-origin H7N9 virus was closely related to H7N9 IAVs isolated in poultry in Nebraska during the same year. Although the precise origin of the wild-bird donor of the turkey-origin H7N9 IAVs could not be determined, our findings suggested that, for both the NA and HA gene segments, the MN turkey-origin H7N9 viruses were related to viruses circulating in wild birds between 2006 and 2011 in the Mississippi Flyway.

Local-scale diversity and between-year "frozen evolution" of avian influenza A viruses in nature

Influenza A virus (IAV) in wild bird reservoir hosts is characterized by the perpetuation in a plethora of subtype and genotype constellations. Multiyear monitoring studies carried out during the last two decades worldwide have provided a large body of knowledge regarding the ecology of IAV in wild birds. Nevertheless, other issues of avian IAV evolution have not been fully elucidated, such as the complexity and dynamics of genetic interactions between the co-circulating IAV genomes taking place at a local-scale level or the phenomenon of frozen evolution. We investigated the IAV diversity in a mallard population residing in a single pond in the Czech Republic. Despite the relative small number of samples collected, remarkable heterogeneity was revealed with four different IAV subtype combinations, H6N2, H6N9, H11N2, and H11N9, and six genomic constellations in co-circulation. Moreover, the H6, H11, and N2 segments belonged to two distinguishable sub-lineages. A reconstruction of the pattern of genetic reassortment revealed direct parent-progeny relationships between the H6N2, H11N9 and H6N9 viruses. Interestingly the IAV, with the H6N9 subtype, was re-detected a year later in a genetically unchanged form in the close proximity of the original sampling locality. The almost absolute nucleotide sequence identity of all the respective genomic segments between the two H6N9 viruses indicates frozen evolution as a result of prolonged conservation in the environment. The persistence of the H6N9 IAV in various abiotic and biotic environmental components was also discussed.

An experimentally determined evolutionary model dramatically improves phylogenetic fit

All modern approaches to molecular phylogenetics require a quantitative model for how genes evolve. Unfortunately, existing evolutionary models do not realistically represent the site-heterogeneous selection that governs actual sequence change. Attempts to remedy this problem have involved augmenting these models with a burgeoning number of free parameters. Here, I demonstrate an alternative: Experimental determination of a parameter-free evolutionary model via mutagenesis, functional selection, and deep sequencing. Using this strategy, I create an evolutionary model for influenza nucleoprotein that describes the gene phylogeny far better than existing models with dozens or even hundreds of free parameters. Emerging high-throughput experimental strategies such as the one employed here provide fundamentally new information that has the potential to transform the sensitivity of phylogenetic and genetic analyses.

Epistatically interacting substitutions are enriched during adaptive protein evolution

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

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

Systematic phylogenetic analysis of influenza A virus reveals many novel mosaic genome segments

Recombination plays an important role in shaping the genetic diversity of a number of DNA and RNA viruses. Although some recent studies have reported bioinformatic evidence of mosaic sequences in a variety of influenza A viruses, it remains controversial as to whether these represent bona fide natural recombination events or laboratory artifacts. Importantly, mosaic genome structures can create significant topological incongruence during phylogenetic analyses, which can mislead additional phylogeny-based molecular evolutionary analyses such as molecular clock dating, the detection of selection pressures and phylogeographic inference. As a result, there is a strong need for systematic screenings for mosaic structures within the influenza virus genome database. We used a combination of sequence-based and phylogeny-based methods to identify 388 mosaic influenza genomic segments, of which 332 are previously unreported and are significantly supported by phylogenetic methods. It is impossible, however, to ascertain whether these represent natural recombinants. To facilitate the future identification of recombinants, reference sets of non-recombinant sequences were selected for use in an automatic screening protocol for detecting mosaic sequences. Tests using real and simulated mosaic sequences indicate that our screening protocol is both sensitive (average >90%) and accurate (average >77%) enough to identify a range of different mosaic patterns. The relatively high prevalence of mosaic influenza virus sequences implies that efficient systematic screens, such as that proposed here, should be performed routinely to detect natural recombinant strains, potential laboratory artifacts, and sequencing contaminants either prior to sequences being deposited in GenBank or before they are used for phylogenetic analyses.

Isolation and identification of a novel rabies virus lineage in China with natural recombinant nucleoprotein gene

Rabies virus (RABV) causes severe neurological disease and death. As an important mechanism for generating genetic diversity in viruses, homologous recombination can lead to the emergence of novel virus strains with increased virulence and changed host tropism. However, it is still unclear whether recombination plays a role in the evolution of RABV. In this study, we isolated and sequenced four circulating RABV strains in China. Phylogenetic analyses identified a novel lineage of hybrid origin that comprises two different strains, J and CQ92. Analyses revealed that the virus 3′ untranslated region (UTR) and part of the N gene (approximate 500 nt in length) were likely derived from Chinese lineage I while the other part of the genomic sequence was homologous to Chinese lineage II. Our findings reveal that homologous recombination can occur naturally in the field and shape the genetic structure of RABV populations.

Evolution of complementary nucleotides in 5' and 3' untranslated regions of influenza A virus genomic segments

The genome of influenza A virus comprises 8 segments (segments 1-8) of single-stranded RNA (virion RNA: vRNA) with negative-polarity. All vRNAs share 13 and 12 terminal nucleotides in the 5' and 3' untranslated regions (UTRs), respectively, which are partially complementary and constitute panhandle and corkscrew structures. Here, it is shown, from the analysis of genomic sequences for 506 strains of influenza A virus, that the number of contiguous complementary nucleotides in the 5' and 3' UTRs varies from 4 to 7 among segments. Complementary nucleotides were segment specific and highly conserved in all segments except for segment 6, where in the phylogenetic analysis co-evolution was observed to have occurred between and within subtypes of neuraminidase (NA). Mutations in the terminal sequences sometimes appeared to have caused convergence between subtypes, involving changes in multiple nucleotide positions. These observations suggest that intra-segmental (homologous) recombinations may have taken place for transferring terminal sequences in segment 6.

Persistence of avian influenza viruses in various artificially frozen environmental water types

Background. This study investigates the viable persistence of avian influenza viruses (AIVs) in various types of artificially frozen environmental water and evaluates the feasibility of similar occurrence taking place in nature, and allowing for prolonged abiotic virus survival, with subsequent biotic viral recirculation. Methods. Fresh, brackish, and salty water, taken in Japan from aquatic biotopes regularly visited by migratory waterfowl, were seeded with AIVs. We monthly monitored the viability of the seeded viruses in the frozen state at -20°C and -30°C, for 12 months. We also monitored virus viability following repeatedly induced freezing and thawing. Results. The viruses exhibited considerable viable persistence all along that period of time, as well as during freezing-thawing cycles. Appreciable, yet noncrucial variances were observed in relation to some of the parameters examined. Conclusions. As typical waterborne pathogens of numerous northerly aquatic birds, AIVs are innately adapted to both the body temperature of their hosts (40°C to 42°C) and, presumably, to subzero temperatures of frozen lakes (down to -54°C in parts of Siberia) occupied and virus-seeded by subclinically infected birds, prior to freezing. Marked cryostability of AIVs appears to be evident. Preservation in environmental ice has significant ecophylogenetic and epidemiological implications, potentially, and could account for various unexplained phenomena.

Viral replication, persistence in water and genetic characterization of two influenza A viruses isolated from surface lake water

Water-borne transmission has been suggested as an important transmission mechanism for Influenza A (IA) viruses in wild duck populations; however, relatively few studies have attempted to detect IA viruses from aquatic habitats. Water-isolated viruses have rarely been genetically characterized and evaluation for persistence in water and infectivity in natural hosts has never been documented. In this study, we focused on two IA viruses (H3N8 and H4N6 subtypes) isolated from surface lake water in Minnesota, USA. We investigated the relative prevalence of the two virus subtypes in wild duck populations at the sampling site and their genetic relatedness to IA viruses isolated in wild waterbirds in North America. Viral persistence under different laboratory conditions (temperature and pH) and replication in experimentally infected Mallards (Anas platyrhynchos) were also characterized. Both viruses were the most prevalent subtype one year following their isolation in lake water. The viruses persisted in water for an extended time period at constant temperature (several weeks) but infectivity rapidly reduced under multiple freeze-thaw cycles. Furthermore, the two isolates efficiently replicated in Mallards. The complete genome characterization supported that these isolates originated from genetic reassortments with other IA viruses circulating in wild duck populations during the year of sampling. Based on phylogenetic analyses, we couldn't identify genetically similar viruses in duck populations in the years following their isolation from lake water. Our study supports the role for water-borne transmission for IA viruses but also highlights that additional field and experimental studies are required to support inter-annual persistence in aquatic habitats.

Host shifts and molecular evolution of H7 avian influenza virus hemagglutinin

Evolutionary consequences of host shifts represent a challenge to identify the mechanisms involved in the emergence of influenza A (IA) viruses. In this study we focused on the evolutionary history of H7 IA virus in wild and domestic birds, with a particular emphasis on host shifts consequences on the molecular evolution of the hemagglutinin (HA) gene. Based on a dataset of 414 HA nucleotide sequences, we performed an extensive phylogeographic analysis in order to identify the overall genetic structure of H7 IA viruses. We then identified host shift events and investigated viral population dynamics in wild and domestic birds, independently. Finally, we estimated changes in nucleotide substitution rates and tested for positive selection in the HA gene. A strong association between the geographic origin and the genetic structure was observed, with four main clades including viruses isolated in North America, South America, Australia and Eurasia-Africa. We identified ten potential events of virus introduction from wild to domestic birds, but little evidence for spillover of viruses from poultry to wild waterbirds. Several sites involved in host specificity (addition of a glycosylation site in the receptor binding domain) and virulence (insertion of amino acids in the cleavage site) were found to be positively selected in HA nucleotide sequences, in genetically unrelated lineages, suggesting parallel evolution for the HA gene of IA viruses in domestic birds. These results highlight that evolutionary consequences of bird host shifts would need to be further studied to understand the ecological and molecular mechanisms involved in the emergence of domestic bird-adapted viruses.

No observed effect of homologous recombination on influenza C virus evolution

The occurrence of homologous recombination in influenza viruses has been under some debate recently. To determine the extent of homologous recombination in influenza C virus, recombination analyses of all available gene sequences of influenza C virus were carried out. No recombination signal was found. With the previous evidence in influenza A and B viruses, it seems that homologous recombination has minimal or no effect on influenza virus evolution.

The re-emergence of H1N1 influenza virus in 1977: a cautionary tale for estimating divergence times using biologically unrealistic sampling dates

In 1977, H1N1 influenza A virus reappeared after a 20-year absence. Genetic analysis indicated that this strain was missing decades of nucleotide sequence evolution, suggesting an accidental release of a frozen laboratory strain into the general population. Recently, this strain and its descendants were included in an analysis attempting to date the origin of pandemic influenza virus without accounting for the missing decades of evolution. Here, we investigated the effect of using viral isolates with biologically unrealistic sampling dates on estimates of divergence dates. Not accounting for missing sequence evolution produced biased results and increased the variance of date estimates of the most recent common ancestor of the re-emergent lineages and across the entire phylogeny. Reanalysis of the H1N1 sequences excluding isolates with unrealistic sampling dates indicates that the 1977 re-emergent lineage was circulating for approximately one year before detection, making it difficult to determine the geographic source of reintroduction. We suggest that a new method is needed to account for viral isolates with unrealistic sampling dates.

A phylogenetic approach to detecting reassortments in viruses with segmented genomes

When multiple strains of viruses with segmented genomes co-infect a single cell, strains with novel genomic constellations may be created. This mutational process, called reassortment, has caused pandemics of influenza A virus in 1957 and 1968. Here a phylogenetic approach to detecting reassortments, which can be used even when the phylogenetic tree constructed for all strains analyzed is unreliable, is presented. A quartet of strains is examined at a time, where a phylogenetic tree is constructed for each genomic segment and the topology is compared among segments only when all quartet trees are supported with a statistical significance. The occurrence of reassortment and the segments involved in the reassortment event are inferred according to the pattern of topological difference among segments. The reassortment point for a pattern is inferred by superimposing the exterior branches of relevant quartet trees on the all-strains trees. In the analysis of H1N1 and H3N2 human influenza A viruses, a topological difference was observed for all pairs of genomic segments, suggesting that there is no pair of segments that has always co-segregated in reassortment during the evolutionary history of these viruses. When the reassortment point was inferred for the pattern of topological difference that was supported with the largest number of quartets for each virus, the results appeared to be mostly correct, suggesting that the method was largely reliable.

Guidelines for identifying homologous recombination events in influenza A virus

The rapid evolution of influenza viruses occurs both clonally and non-clonally through a variety of genetic mechanisms and selection pressures. The non-clonal evolution of influenza viruses comprises relatively frequent reassortment among gene segments and a more rarely reported process of non-homologous RNA recombination. Homologous RNA recombination within segments has been proposed as a third such mechanism, but to date the evidence for the existence of this process among influenza viruses has been both weak and controversial. As homologous recombination has not yet been demonstrated in the laboratory, supporting evidence, if it exists, may come primarily from patterns of phylogenetic incongruence observed in gene sequence data. Here, we review the necessary criteria related to laboratory procedures and sample handling, bioinformatic analysis, and the known ecology and evolution of influenza viruses that need to be met in order to confirm that a homologous recombination event occurred in the history of a set of sequences. To determine if these criteria have an effect on recombination analysis, we gathered 8307 publicly available full-length sequences of influenza A segments and divided them into those that were sequenced via the National Institutes of Health Influenza Genome Sequencing Project (IGSP) and those that were not. As sample handling and sequencing are executed to a very high standard in the IGSP, these sequences should be less likely to be exposed to contamination by other samples or by laboratory strains, and thus should not exhibit laboratory-generated signals of homologous recombination. Our analysis shows that the IGSP data set contains only two phylogenetically-supported single recombinant sequences and no recombinant clades. In marked contrast, the non-IGSP data show a very large amount of potential recombination. We conclude that the presence of false positive signals in the non-IGSP data is more likely than false negatives in the IGSP data, and that given the evidence to date, homologous recombination seems to play little or no role in the evolution of influenza A viruses.

Evidence for a novel gene associated with human influenza A viruses

BACKGROUND

Influenza A virus genomes are comprised of 8 negative strand single-stranded RNA segments and are thought to encode 11 proteins, which are all translated from mRNAs complementary to the genomic strands. Although human, swine and avian influenza A viruses are very similar, cross-species infections are usually limited. However, antigenic differences are considerable and when viruses become established in a different host or if novel viruses are created by re-assortment devastating pandemics may arise.

RESULTS

Examination of influenza A virus genomes from the early 20th Century revealed the association of a 167 codon ORF encoded by the genomic strand of segment 8 with human isolates. Close to the timing of the 1948 pseudopandemic, a mutation occurred that resulted in the extension of this ORF to 216 codons. Since 1948, this ORF has been almost totally maintained in human influenza A viruses suggesting a selectable biological function. The discovery of cytotoxic T cells responding to an epitope encoded by this ORF suggests that it is translated into protein. Evidence of several other non-traditionally translated polypeptides in influenza A virus support the translation of this genomic strand ORF. The gene product is predicted to have a signal sequence and two transmembrane domains.

CONCLUSION

We hypothesize that the genomic strand of segment 8 of encodes a novel influenza A virus protein. The persistence and conservation of this genomic strand ORF for almost a century in human influenza A viruses provides strong evidence that it is translated into a polypeptide that enhances viral fitness in the human host. This has important consequences for the interpretation of experiments that utilize mutations in the NS1 and NEP genes of segment 8 and also for the consideration of events that may alter the spread and/or pathogenesis of swine and avian influenza A viruses in the human population.

Identifying changes in selective constraints: host shifts in influenza

The natural reservoir of Influenza A is waterfowl. Normally, waterfowl viruses are not adapted to infect and spread in the human population. Sometimes, through reassortment or through whole host shift events, genetic material from waterfowl viruses is introduced into the human population causing worldwide pandemics. Identifying which mutations allow viruses from avian origin to spread successfully in the human population is of great importance in predicting and controlling influenza pandemics. Here we describe a novel approach to identify such mutations. We use a sitewise non-homogeneous phylogenetic model that explicitly takes into account differences in the equilibrium frequencies of amino acids in different hosts and locations. We identify 172 amino acid sites with strong support and 518 sites with moderate support of different selection constraints in human and avian viruses. The sites that we identify provide an invaluable resource to experimental virologists studying adaptation of avian flu viruses to the human host. Identification of the sequence changes necessary for host shifts would help us predict the pandemic potential of various strains. The method is of broad applicability to investigating changes in selective constraints when the timing of the changes is known.

Influenza A's natural reservoir is waterfowl. Sometimes avian virus genomic segments are able to shift to a human host, either in toto or by combining with those that underwent a previous host shift event. Such host shift events can cause worldwide pandemics in their immunologically naive hosts. In order for these host shifts to establish a stable lineage, the virus has to adapt to the new host. Identifying the changes that have occurred in the past can provide important clues about how this process happens, and how surveillance for new influenza threats should be targeted. Unfortunately, it is difficult to determine whether an amino acid has changed due to adaptation to the new host or whether the change occurred through random drift. Here we describe a novel phylogenetic approach to identifying locations where the nature of the selective pressure exerted on the location has changed corresponding to the host shift event. We identify a set of locations on a number of the genomic segments. The approach we describe is of wide applicability when the timing of the change of selective constraints is known in advance.

Reassortment patterns in Swine influenza viruses

Three human influenza pandemics occurred in the twentieth century, in 1918, 1957, and 1968. Influenza pandemic strains are the results of emerging viruses from non-human reservoirs to which humans have little or no immunity. At least two of these pandemic strains, in 1957 and in 1968, were the results of reassortments between human and avian viruses. Also, many cases of swine influenza viruses have reportedly infected humans, in particular, the recent H1N1 influenza virus of swine origin, isolated in Mexico and the United States. Pigs are documented to allow productive replication of human, avian, and swine influenza viruses. Thus it has been conjectured that pigs are the “mixing vessel” that create the avian-human reassortant strains, causing the human pandemics. Hence, studying the process and patterns of viral reassortment, especially in pigs, is a key to better understanding of human influenza pandemics. In the last few years, databases containing sequences of influenza A viruses, including swine viruses, collected since 1918 from diverse geographical locations, have been developed and made publicly available. In this paper, we study an ensemble of swine influenza viruses to analyze the reassortment phenomena through several statistical techniques. The reassortment patterns in swine viruses prove to be similar to the previous results found in human viruses, both in vitro and in vivo, that the surface glycoprotein coding segments reassort most often. Moreover, we find that one of the polymerase segments (PB1), reassorted in the strains responsible for the last two human pandemics, also reassorts frequently.

Origin and phylodynamics of rabbit hemorrhagic disease virus

To determine the origin, phylogenetic relationships, and evolutionary dynamics of rabbit hemorrhagic disease virus (RHDV), we examined 210 partial and complete capsid gene nucleotide sequences. Using a Bayesian Markov chain Monte Carlo approach, we estimated that these sequences evolved at a rate of 3.9 x 10(-4) to 11.9 x 10(-4) nucleotide substitutions per site per year. This rate was consistent across subsets of data, was robust in response to recombination, and casts doubt on the provenance of viral strains isolated from the 1950s to the 1970s, which share strong sequence similarity to modern isolates. Using the same analysis, we inferred that the time to the most recent common ancestor for a joint group of RHDV and rabbit calicivirus sequences was <550 years ago and was <150 years ago for the RHDV isolates that have spread around the world since 1984. Importantly, multiple lineages of RHDV were clearly circulating before the major Chinese outbreak of 1984, a finding indicative of an early evolution of RHDV virulence. Four phylogenetic groups within RHDV were defined and analyzed separately. Each group shared a common ancestor in the mid-1960s or earlier, and each showed an expansion of populations starting before 1984. Notably, the group characterized by the antigenic variant RHDVa harbors the greatest genetic diversity, compatible with an elevated fitness. Overall, we contend that the high virulence of RHDV likely evolved once in the early part of the 20th century, well before the documented emergence of rabbit hemorrhagic disease in 1984.

Panorama phylogenetic diversity and distribution of type A influenza viruses based on their six internal gene sequences

Background

Type A influenza viruses are important pathogens of humans, birds, pigs, horses and some marine mammals. The viruses have evolved into multiple complicated subtypes, lineages and sublineages. Recently, the phylogenetic diversity of type A influenza viruses from a whole view has been described based on the viral external HA and NA gene sequences, but remains unclear in terms of their six internal genes (PB2, PB1, PA, NP, MP and NS).

Methods

In this report, 2798 representative sequences of the six viral internal genes were selected from GenBank using the web servers in NCBI Influenza Virus Resource. Then, the phylogenetic relationships among the representative sequences were calculated using the software tools MEGA 4.1 and RAxML 7.0.4. Lineages and sublineages were classified mainly according to topology of the phylogenetic trees and distribution of the viruses in hosts, regions and time.

Results

The panorama phylogenetic trees of the six internal genes of type A influenza viruses were constructed. Lineages and sublineages within the type based on the six internal genes were classified and designated by a tentative universal numerical nomenclature system. The diversity of influenza viruses circulating in different regions, periods, and hosts based on the panorama trees was analyzed.

Conclusion

This study presents the first whole views to the phylogenetic diversity and distribution of type A influenza viruses based on their six internal genes. It also proposes a tentative universal nomenclature system for the viral lineages and sublineages. These can be a candidate framework to generalize the history and explore the future of the viruses, and will facilitate future scientific communications on the phylogenetic diversity and evolution of the viruses. In addition, it provides a novel phylogenetic view (i.e. the whole view) to recognize the viruses including the origin of the pandemic A(H1N1) influenza viruses.

Reassortment patterns in Swine influenza viruses

Previous human influenza pandemics were the results of emerging viruses from non-human reservoirs, with at least two caused by strains of mixed human and avian origin. Also, many cases of swine influenza viruses have reportedly infected humans, including the recent human H1N1 strain, isolated in Mexico and the United States. Pigs are documented to get infected with human, avian, and swine viruses and allow productive replication, thus it has been conjectured that they are the "mixing vessel" that create reassortant strains, causing the human pandemics. In this paper, we apply several statistical techniques to an ensemble of publicly available swine viruses to study the reassortment phenomena. The reassortment patterns in swine viruses confirm previous results found in human viruses that the glycoprotein coding segments reassort most often. Moreover, one of the polymerase segments (PB1), reassorted in the strains responsible for the last two human pandemics of 1957 and 1968, also reassorts frequently.

Homologous recombination is apparent in infectious bursal disease virus

Infectious bursal disease virus (IBDV) is a non-enveloped double-stranded RNA virus belonging to the Birnaviridae family. It shows substantial variation in the major antigen region of the viral capsid protein VP2, where a hypervariable region plays a key role in the virulence of IBDV and its epitope. This study identifies several putative recombinants from previously published data to suggest that homologous recombination may naturally occur between different IBDV strains. In addition, a novel very virulence sublineage emerges in the VP2 phylogenic tree, comprising three putative recombination strains isolated in Korea and China, KSH, KK1 and SH-h. The major putative parents of the three mosaics are descended from the vaccine lineage while their hypervariable regions from vvIBDV. These findings also suggest that vaccine coverage may have influence on the evolution and genetic diversity of IBDV, resulting in a novel group with vvIBDV phenotype through recombination with wild IBDV.