At the close of World War I, the H1N1 influenza A virus swept the world[1]. During the 1917–1918 pandemic, approximately 40% of the human population was infected, with a death rate above 2%. It is estimated that this virus killed more people than died in the world war that was just then ending. Mortality rates have dramatically declined since then[2], but the H1N1 flu has persisted. As a zoonotic pathogen, the influenza virus is able to infect multiple species. It is generally thought that aquatic waterfowl are a primary natural influenza reservoir[3], where there are usually no clinical symptoms[4], and where low level transmission probably perpetuates the viral pool[5]. All 14 influenza subtypes are maintained in waterfowl[5].

H1N1 has had an interesting history. Derivatives of the original virus circulated in humans and swine until 1957, when the human strain went extinct. In 1977, a version identical to those circulating in NE Europe in the early 1950s reappeared in Anshan, China and subsequently spread across the world[5–7]. In 2009, a swine H1N1 jumped to the human population, causing a widespread pandemic. This has increased concern that H1N1 might mutate into a more virulent form. However, since the pandemic of 1917, this has not happened. In fact, H1N1-related human mortality has declined very dramatically and very systematically[2]. Apart from the 1917 pandemic, H1N1 has failed to cause any severe global pandemic, and human H1N1 essentially went extinct from 1957–1977. Since its re-introduction, it has remained a relatively minor cause of influenza mortality[2]. This applies also to the 2009 outbreak, which caused relatively few deaths in those areas with good reporting systems in place[8].

It is therefore reasonable to ask if the striking reduction in H1N1 mortality might be due, in part, to natural attenuation resulting from deleterious mutation accumulation. Herd immunity is undoubtedly an important factor in reduced H1N1 mortality since 1918, but this may not be sufficient to explain the continuous decline in H1N1-related mortality over multiple human generations or the eventual extinction of the viral strain. Likewise, improved medical treatments, such as antibiotic treatment for flu-related pneumonia, were certainly a significant factor reducing H1N1 mortality, but these do not appear to fully explain the nature of the pattern of mortality decline seen for H1N1. For example, the exponential decline in mortality began before the invention of antibiotic treatment.

The literature suggests RNA viruses should be inherently subject to mutational degeneration[9–13]. This includes the bacteriophage MS2[14], the tobacco etch virus[15], HIV[16–19], dengue virus type-2[20], Ebola[21, 22], and SARS[23, 24]. Some have suggested that intentionally increasing the rate of mutation accumulation (“lethal mutagenesis”) may be a way to control viral epidemics by hastening strain extinction[25–30]. There is some long-term historical evidence that supports the concept of natural viral attenuation through mutation accumulation[2], and theoretical studies using numerical simulation strongly support the concept of natural and accelerated genetic attenuation of RNA viruses[13].

Previous genetic studies examining the history of the influenza virus have performed extensive phylogenetic analyses of influenza genomes[8, 31–35]. They have shown considerable nucleotide diversity among circulating strains, given clear evidence for adaptive selection of antigenic variants[36–42], and have shown that most of the major innovations within the flu genome have occurred via reassortment[5], by which one flu strain has recombined with another strain and obtained a segment of RNA from the second strain.

Influenza phylogenies are odd, in some respects, consisting of one main trunk with many short and short-lived side branches[36, 39, 42, 43]. There are several factors that influence strain diversity, including the mutation rate, selection, drift, and cross-immunity-mediated competition between strains[39]. Feguson, et al.[39] felt that strain-transcending immunity was essential in restricting viral diversity. Ito, et al.[42] thought that the many side branches died out because they could not compete with the main-trunk lineage viruses. Others have noted the temporal extinction of circulating strains upon the introduction of new strains or serotypes[5]. Whatever the reason, a hierarchy of strain robustness is evident in the history of influenza viruses.

While phylogenetic studies can build robust family trees, they do so by focusing only on a limited number of “informative” genomic locations[44]. Even though the influenza genome is broken into eight separate RNAs, unless an individual is infected with two strains simultaneously (providing an opportunity for reassortment), all eight sections are inherited as a set in a form of linkage[36]. Thus, neutral and slightly deleterious mutations are carried along with those mutations under positive or purifying selection. This gives us enough information to make many phylogenetic inferences, and we have a wealth of data telling us the history of the various viral lines, but these phylogenetic techniques ignore genetic change within a larger portion of the genome in order to focus on the phylogenetically-informative sites.

There is abundant evidence for multiple substitutions in specific places that are, in turn, sites of active selection[42], but what about the rest of the genome? Fitch, et al.[36] noted that the main strategy of influenza viruses was to outrun the immune system of the host by maintaining a high mutation rate. They wondered “why such a virus does not accumulate so many deleterious mutations as to die of its own ineptitude.” They further wondered if the only reason influenza viruses continued to circulate in mammalian species was because of a sort of “gene therapy” through reassortment with stronger, less mutated avian viruses.

The present study is an attempt to elucidate more fully the nature of mutational change in this historically important human pathogen.

Accession numbers for all available complete flu genomes were obtained from the Influenza Research Database[45] as of June 1, 2012. This list was then compared to those available at FluGenome.org[46]. To create a list of the 2009–2010 H1N1 outbreak versions, we used the collated genome list from Kedwaii, et al.[8]. Using these two accession lists, sequence data were obtained from GenBank[47]. We removed non-human/swine, non-H1N1, and incomplete genomes from this list, obtaining 3,755 human and 351 swine genomes. Because there are few indels in these flu genomes, a simple first-pass alignment algorithm, written in Perl with the BioPerl toolkit[48], was used. Essentially, we visually picked out a 6-10-nucleotide candidate region close to the beginning of each of the eight strands of the flu genome that appeared mostly conserved. The algorithm scanned each set of sequences for this segment and added an appropriate number of spaces to the beginning of each for a rough alignment. Using BioEdit[49], we simply looked for sequences that were misaligned and inserted gaps where necessary to complete the alignment, paying close attention to codon usage. In difficult and ambiguous regions, we chose to minimize the number of changes necessary to produce each variant, similar to the methodology of Carter[44].

Once an alignment for each of the eight sections of the flu genome was completed, we used another Perl program to scan through the alignment and compare each sequence to a reference strain. Point mutations were counted individually while indels, irrespective of size, were counted as a single change[44]. There are several different ways to calculate sequence differences (discussed in[50]), but this simple measure is sufficient for the purposes of this paper. Due to variations in technique among the various studies that produced these sequences, the ends of the alignment tended to be ragged. For this reason, we compared only the section of complete alignment between each sequence and the reference genome. The few internal reassortments[32] that have occurred in H1N1 were easily seen in the alignments, but these also made up a small fraction of all sites. We did our best to create a consistent alignment through these regions.

We began by analyzing mutation accumulation during the human H1N1 outbreak of 2009–2010, using strain California/04/2009 as a reference. This was the first sequenced genome of the outbreak included in Kedwaii, et al.[8]. We then calculated the extent of mutational divergence of all subsequent genomes over the duration of the outbreak.

We did a similar, but much larger analysis of mutation accumulation in all human H1N1 strains from 1918 to 2012, using the 1918 “Brevig Mission” reference genome (A/Brevig Mission/1/1918) as a reference. This sequence was missing the last 480 nucleotides of segment 4. We amended this by appending the ending sequence of another 1918 H1N1 virus, A/South Carolina/1/18(H1N1). These two strains had only one nucleotide difference in the ORF of segment 4 (C/T at position 238). Thus, they are essentially identical and we do not expect to have introduced extra mutations in this way.

Using the amended 1918 Brevig Mission virus as a reference and including all human and porcine viruses in the database, we calculated SNPs, indels, transitions, transversions, non-synonymous amino acid changes, dN/dS ratios, predicted protein lengths (for all 11 proteins), the normalized codon scores (NCS) and relative synonymous codon usage (RSCU)[51] score for each predicted protein of each genome. To calculate the normalized codon score, we downloaded the frequency of each of the 64 codons in the human, duck and domestic pig genomes from the Codon Usage Database[52]. We then created a predicted amino acid chain from each of the respective open reading frames. We calculated NCS for each protein by simply summing the frequency of each codon, as it is used in each of the three species, and then dividing by the total number of codons in the proteins produced. RSCU is simply the number of times any particular codon appears, divided by the number of times synonymous codons appear in that genome. Anhlan, et al.[51] were the first to apply RCSU to influenza genomic data, but they only studied one of the eleven protein products of the H1N1 genome and did not analyze changes over time. Here, we extend their techniques to create a time series in which one can better study shifts in codon usage due to mutation accumulation.

During the 2009–2010 H1N1 outbreak, mutations within H1N1 accumulated at a relatively constant rate (Figure1). The strains sequenced at the end of the pandemic had approximately 80 more mutations than the 2009 reference genotype. The line of best fit reveals a slope of 42.1 mutations accumulated per year (3.2 × 10^-3 mutations per site per year). Given that the data were collected from viral strains circulating worldwide, the correlation between mutation count and the date of collection was surprisingly good (r² = 0.77).

The linear accumulation of mutations during the 2009 outbreak clearly also extends to the longer-term accumulation of mutations during the entire history of H1N1 (Figure2). During the last century there has been a remarkably constant increase in mutation count within the H1N1 virus population, except for a striking discontinuity between 1957 and 1976. This discontinuity reflects the extinction of the human H1N1 strain in the mid 1950s, followed by the re-introduction of the strain in 1976, presumably from a researcher’s freezer. Even though most reports (e.g.,[53]) indicate the re-introduction year was 1977, our data suggests the year of re-introduction was 1976. This is based upon the isolate A/New Jersey/1976, which very clearly falls in line with the rest of the re-introduced lineage. Notice that after re-introduction in 1976, the human H1N1 mutation count and the rate of accumulation resumed exactly where it left off just before the extinction occurred in the 1950s.

Our data clearly show that some non-frozen H1N1 genotypes occasionally appeared in the human population after H1N1 dropped from public concern after 1957 (Figure2). Nine H1N1 strains that do not belong to the “frozen” lineage arose in the human population between 1976 and the 2009 H1N1 outbreak (A/Wisconsin/301/1976, A/New Jersey/8/1976, A/California/10/1978, A/California/45/1978, A/Albany/8/1979, A/Memphis/1/1979, A/USSR/46/1979, A/Ohio/3559/1988, and A/Maryland/12/1991). These fall on the same mutation-count trajectory as the H1N1 circulating before the 1957 extinction. However, these nine strains appear to represent repeated transmission events from pigs to humans that failed to cause any pandemic (the porcine lineage had no extinction event, and hence no pause in mutation accumulation). Nearest neighbor calculations (data not shown) indicate these strains are not a continuation of the human lineage. They cluster tightly with the 2009–2010 outbreak porcine viruses, all of which are more closely related to the 1918 isolate than they are to the human lineage viruses. Thus, we included the 2009–2010 viruses and the nine isolated, non-outbreak viruses in the “porcine” category. In this same figure, we can see that genotypes from the 2009 outbreak and after (Figure2, circle) fall directly on the trajectory of the non-frozen lineage. This clearly shows that the 2009 genotype was not derived from the 1976 re-introduction virus. Kedwaii, et al.[8] affirmed earlier conclusions that the 2009 genotype was due to a reassortment between two swine viruses, an H1N2 and an H1N1, from different continents, but this did not affect the mutation accumulation curve. Regardless, general attenuation due to genome-wide mutation accumulation might best explain the very low mortality[54–56] associated with the 2009 pandemic. The earliest 2009 outbreak strain had already accumulated 1,889 mutations compared to the 1918 strain.

The most obvious place to test how reassortment might affect the mutation accumulation curve is with the 2009–2010 “swine flu” outbreak samples. Despite the fact that the viral line has been attributed to reassortment with swine viruses, with partial contribution from earlier reassortments between swine and avian strains[8], it does not display a jump in mutation count. Reassortment can produce novel antigenic variants, but it does not reverse the majority of mutations, for they have accumulated in the non-reassorted areas of the genome.

After we identified the direct 1918 human lineage viruses, we were able to adjust the dates of the later samples to account for a period of dormancy (Figure3). When we shifted the “frozen” samples back 21 years (i.e., sequences from 1976 were assigned a recalibrated date of 1955), we see the data lines up almost perfectly (r² = 0.989; mutation rate = 14.4). The exact amount of time-frame adjustment was not critical; adjustments ranging from 20 to 30 years all gave excellent alignments, with correlation coefficients above 98%. The observed accumulation rate of roughly 14 mutations per year (1.1 × 10^-3 mutations per site per year) is more or less consistent with previous estimates (e.g.,[39, 57]). The accumulation of non-synonymous mutations in human H1N1 genomes is also shown in Figure3. The non-synonymous mutations are not as abundant, as expected, but they are still accumulating at a rapid and constant rate. There were an average of 2.4 amino acid changes per year (0.6 × 10^-3 mutations per codon per year).

Our data show that the lineage re-introduced in 1976 was the only significant human H1N1 strain from 1976 until 2009. The mutation count trajectory of the H1N1 lineage from the 1976 re-introduction stops abruptly in 2009 and does not reappear in any of the 353 H1N1 genomes published from mid-2009 through mid-2012. Human H1N1, the remnant of the 1918 lineage, appears to have gone extinct a second time.

The data in Figure3 can be used to estimate the time of arrival into the human population of the first H1N1 virus. Extrapolating total mutation count linearly backwards to zero yields a year of first introduction into the human population of approximately 1893. Interestingly, this would be in time to explain the 1889–1890 “Russian” flu outbreak, which has been considered by some to possibly be the first H1N1 outbreak (no sequence or serotype data are available for that event). While our data strongly suggest that the ancestral genotype for all human and porcine H1N1 influenza strains (and by extension, all human influenza strains), entered the human population as early as the 1890s, others have estimated the arrival date to be one or more decades later[5]. Our data strongly confirm that all human H1N1 viruses derive from a single, very recent common ancestor, just as others have concluded prior to this study[1]. Given an estimated arrival date of 1893, we infer that, by 1918, H1N1 would have accumulated roughly 375 (25 years × 15 mutations/year) mutations since the initial genotype invaded the human population. This depends, of course, on the exact arrival date.

Within the entire human H1N1 lineage, only 56% of all nucleotide sites were invariant, indicating intense mutation pressure. Before its apparent extinction in 2009, the human H1N1 lineage had accumulated approximately 1400 mutations (mostly point mutations), including 320 non-synonymous mutations, compared to the 1918 genotype. Across all 3,755 H1N1 genomes isolated from humans (i.e., including the outbreak and non-outbreak porcine versions), we found that only 41.8% of all nucleotide sites were invariant. Multiple variant forms (i.e., 3-fold and 4-fold degeneracy, meaning 3 and sometimes 4 different nucleotides observed at a given site) occurred at 21% of the variable loci, further suggesting that the mutations were primarily non-adaptive. When we included the porcine H1N1 genotypes, only 36.1% of sites remained invariant, with a similar distribution of nucleotide degeneracy. Because the re-introduced H1N1 lineage had accumulated fewer mutations than the porcine H1N1 strains, the disappearance of the human lineage and the subsequent dominance of the 2009 porcine strain in the human population resulted in even more divergence from the 1918 stain: over 1900 nucleotide differences (approximately 15% divergence) and 325 amino acids differences (approximately 7% divergence in protein sequence).

The approximately 90-nucleotide PB1-F2 protein is coded on the second genomic segment and overlaps completely the PB1 coding region[58, 59] Table1. PB1-F2 has been implicated in increasing host death rate. Yet, not all influenza viruses code for this protein, so it is not intrinsically necessary for the cycle of infection. In the case of H1N1, the majority (98.2%) of genomes in this study are predicted to fail to produce a full-length version of the protein, including most of the human versions after 1948 and none of the porcine versions prior to 1987. For this reason, only 10 of the 11 influenza proteins were included in RSCU and NCS calculations. All genomes were predicted to produce fully-formed PB2, PB1, PA, HA, NP, NA, M1, M2, and NS1 proteins. Only 14 genomes were expected to fail to produce a functioning NS2 protein due to mutation in the start codon, but all of these had another potential in-frame ATG only a few codons upstream and/or downstream. For this reason, we assumed these viruses would also produce an NS2 protein of proper length.

Mutation accumulation in swine H1N1 since 1918 is shown in Figure4, with the human regression line (with “frozen” linage adjusted for 21 years of dormancy) for comparison. Regression analysis (ANCOVA) indicated a significant difference in the slope (p = 0.0001) of the lines of best fit for the viruses obtained from the two species. Although it is clear that the swine versions have more variability towards the end of the sampling period, they also have a higher mutation rate in general, even after adjusting the human data for the 21 year period of dormancy.

We note that many mutations appeared and disappeared several times over the life of H1N1. The data are too sparse over most of the range to know if this is due to repeat mutations at common loci or incomplete genetic sampling of the circulating H1N1 viruses. However, over the 2009–2010 H1N1 outbreak, Kedwaii, et al.[8] demonstrated near-complete lineage loss among the circulating viruses during a single season. Bedford, et al.[33] had similar results in their study of influenza A (H3N2). It was noted earlier that the phylogeny appears as a main trunk with many short-lived side branches. If most mutations are lost from the population quickly, variations that appear several times over the course of H1N1 are most likely due to mutations occurring multiple times in the same place. We do not have enough sequence data to get a full picture of strain diversity in most years, but this will be interesting to pursue as more data become available.

Even though each genomic segment has a slightly different mutation/substitution rate, strong, genome-wide patterns are evident, including a declining ratio of non-synonymous to synonymous codon changes (dN/dS), a declining ratio of transition to transversion mutations (Ti/Tv), and even changes in codon bias. Additional file1 includes these data on a segment-by-segment basis.

Figure5 illustrates a continuous decline in the Ti/Tv ratio during the entire history of the human H1N1 lineage. This suggests there may be increasing selection against transversions over time, perhaps reflecting some degree of truncation selection as the genome becomes more attenuated and as extinction approaches. A transversion is 2.7 times more likely to cause a non-synonymous change than a transition, and is also more likely to change RNA 3-dimensional configuration. Since we do not see a change in the rate of accumulation of amino acid substitutions over time (Figure3), we suspect that the heightened selection against transversions may be more strongly associated with RNA architectural constraints.

Over the history of the human H1N1 lineage, the accumulation of non-synonymous mutations was highly linear, after correcting the “frozen” samples for 21 years of dormancy (Figure3). While it is true that this simple distance metric applied to amino acid changes would be expected to underestimate evolutionary relationships at large distances due to an increasing likelihood of multiple changes at individual loci, applying the standard Poisson or Gamma corrections[60] made no improvement to the fit of the regression to two significant figures. We therefore conclude that the relationship is indeed linear, at least over the sampled timescale.

Analysis of the normalized codon scores (NCS) clearly shows that codon preference degenerated continuously throughout the entire history of the human H1N1 line (Figure6). All viruses in our database are more duck-like in their codon usage and most (99.8%) are more human-like than pig-like, but codon preference moved toward randomized codon use and away from the codon preference of the human, swine and duck hosts. Thus, the virus is not adapting to use the codon preference of any of the host genomic environments, even when residing in a single species for 70 years (the duration of the human H1N1 lineage, after correcting for 21 years of dormancy) with no evidence for the introduction of new genetic material through reassortment. These data are further evidence that the genome-wide changes are not adaptive, and are in fact degenerative. Non-synonymous changes had little effect on the NCS scores. The 345 non-synonymous amino acid changes between 1918 and 2009 (A/Pensacola/INS235/2009) had a negligible effect (to three significant figures) on the final result for this sequence.

Relative synonymous codon usage (RCSU) patterns support the attenuation hypothesis as well (Figure7, Additional file2). Most codons changed in frequency in a highly consistent manner over the life of the 1918 human H1N1 strain. Those with the greatest rate of positive change over time were CAA (glutamine), ATC (isoleucine), and GTA (valine), and two stop codons, TAA and TAG. Each of these was balanced by a strong decrease in a synonymous codon: TGA (stop), ATT (isoleucine), CAG (glutamine), and GTG (valine).These changes are due to simple transition mutations, possibly associated with some mutational bias. The most stable codons all encoded either serine (TCA, TCC, TCT), or threonine (ACT, ACC, ACA, ACG). More than half of all codons (33 out of 62) were changing in frequency in a direction away from the average human codon usage, including three of the four in Figure7. All but one codon (ACC) had a regression line significantly different from a slope of zero at the p<0.05 level, and most (56 out of 62) scored p<0.0001. Thus, both RCSU and NCS reveal a clear and systematic erosion of codon bias over time.

It has generally been assumed that any non-neutral mutations within the influenza genome have arisen as selective adaptations and generally help drive influenza toward a stronger and more dangerous pathogen (in terms of either pathenogenicity or transmissibility)[41]. This was probably the basis for the extreme caution exhibited during the 2009–2010 H1N1 outbreak. There is a general perception that, given enough time, H1N1 might mutate into a stronger pathogen, and hence might create another catastrophic pandemic, as it did in 1918.

Selection is evident in H1N1 and other influenza genomes. In certain sites of the HA1 genomic segment of H3, non-synonymous substitutions occur at more than twice the rate as in other sites[39]. This is seen as a major signature of adaptive change, but deleterious mutations in the areas not under selection are carried along with the ones under positive selection[36]. In this light, some perceive H1N1 to be a growing threat, with a new outbreak being just a matter of time. Despite this common perception, a more lethal version of H1N1 has not arisen via mutation within the human population during the last 90+ years. This is significant. The two major human influenza pandemics since 1918 did not arise due to mutations within H1N1, but arose via horizontal transmission of new genetic material from bird influenza strains, creating recombinant viruses. They were also less lethal than the 1918 version. It is true that the population had a degree of residual immunity and was not as immunologically naïve as it was in 1917–18, but selection has still not been able to generate a devastating pandemic from the remnants of that which swept the world at the close of WWI.

In this paper, we examine an alternative point of view regarding mutation accumulation within H1N1. We suggest that, while specific adaptive mutations commonly occur within the H1N1 virus, many more deleterious mutations are accumulating than beneficial mutations, even when there is strong selection. Consequently, H1N1 appears to have been in very gradual error catastrophe throughout its history.

Our results strongly confirm the widely recognized fact that all past and present human and swine H1N1 influenza strains derive from the 1918 strain. By extension, this applies to other human influenza strains, including H3N2 and H2N2[1].

Human influenza has such a high mutation rate that, even within a single host individual, the virus quickly becomes a genetically diverse “swarm”[61, 62]. Yet, globally human H1N1 influenza is monophyletic and all current variation is of recent origin. This is only possible if almost all human influenza lineages rapidly go extinct. Moreover, we present strong evidence that the H1N1 genome has been systematically degenerating since 1918. This is evidenced by continuous, systematic, and rapid changes in the H1N1 genome throughout its history. For example, there was an especially rapid and monotonic accumulation of mutations during a single pandemic (Figure1). Similarly, there was a continuous and rapid accumulation of mutations over the entire history of the virus (Figures2 and3), including a similar steady increase in non-synonymous amino acid substitutions (Figure3). While mutations accumulated in the human H1N1s, there was a parallel accumulation of mutations in the porcine H1N1 lineage (Figure4). Fitch, et al.[36] also showed a linear mutation accumulation curve. The “gnarled trunk” of Ito, et al.[42] should be buried in our data as well, but it will be swamped by the pervasive, genome-wide accumulation of mutations not under active selection.

Within the human lineage, there was a systematic decline in the Ti/Tv ratio (Figure5). In addition, there was a very consistent loss of codon specificity over time (Figures6 and7). We show compelling new evidence supporting the extinction of human H1N1 in the 1950s, its subsequent re-introduction in 1976, and an apparent second extinction event of the human H1N1 lineage in 2009. Strain extinction has often occurred when new strains appeared, including H1N1 replacing the circulating H3-like strains in 1917, H2N2 replacing H1N1 in 1957, and H3N2 replacing H2N2 in 1968[5]. To our knowledge, we are the first to document the replacement of the re-introduced human H1N1 with reassorted swine H1N1 in 2009. All this is consistent with the genetic attenuation hypothesis, and we feel this is the most fundamental explanation for the very smooth, systematic, and exponential decline in H1N1 mortality rates since 1918[2].

Sequence analysis of historical and modern genomes of influenza H1N1 reveal a great deal about the history of the virus. The most recent common ancestor existed only about 120 years ago, and there has been universal extinction of all earlier human influenza strains. The rate of mutational divergence from that original genotype has been very constant both in human and porcine H1N1 strains (roughly 14–16 mutations fixed per year). Modern H1N1 strains have diverged from the original genotype by roughly 1,900 fixations (more than 15% of the genome has mutated). Roughly 7% of the amino acids have mutated. Mutation accumulation is also associated with the historical exponential decline in H1N1 human mortalities, which may suggest significant genetic attenuation, as might arise due to very slow error catastrophe.

Sequence analysis confirms that the human H1N1 strain went extinct in 1957, but was reintroduced in 1976, apparently from a specimen frozen in the early 1950s. The resulting “frozen” lineage appeared less mutated compared to all other contemporary H1N1 strains (i.e., porcine strains and a few rare H1N1 strains appearing in humans but not derived from the frozen strain). Consistent with the genetic attenuation model, the frozen human H1N1 lineage disappeared in 2009, and may now be extinct.

It appears that the H1N1 strains currently in circulation are significantly attenuated, and cannot reasonably be expected to back-mutate into a non-attenuated strain. The greatest influenza threat, therefore, is the introduction of a non-attenuated strain from some natural reservoir. This suggests that a better understanding of the origin of such non-attenuated strains should be a priority[5]. Our findings suggest that new strategies that accelerate natural genetic attenuation of RNA viruses may prove useful for managing future pandemics and, perhaps in the long run, precluding the genesis of new influenza strains.