Reassortment Patterns in Swine Influenza Viruses

Abstract

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.

About the Authors:

Hossein Khiabanian

* E-mail: [e-mail]

Affiliation: Department of Biomedical Informatics and Center for Computational Biology and Bioinformatics, Columbia University College of Physicians and Surgeons, New York, New York, United States of America

Vladimir Trifonov

Affiliation: Department of Biomedical Informatics and Center for Computational Biology and Bioinformatics, Columbia University College of Physicians and Surgeons, New York, New York, United States of America

Raul Rabadan

Affiliation: Department of Biomedical Informatics and Center for Computational Biology and Bioinformatics, Columbia University College of Physicians and Surgeons, New York, New York, United States of America

Introduction

Pandemics are epidemics that rapidly spread on a worldwide scale, caused by pathogens against which humans have no immunity that infect a large part of the population and lead to associated serious illnesses. Human influenza pandemics are caused by emerging influenza viruses from non-human reservoirs. From the three influenza pandemics of the twentieth century, the 1918 pandemic was possibly caused by an influenza virus with an avian origin [1], [2] and the other two, in 1957 and 1968, were caused by new strains that were combinations of avian and human viruses through the process of reassortment [3], [4].

There also have been many cases of swine influenza viruses infecting humans [5], [6]. In particular, in March 2009, a new human H1N1 influenza A virus of swine origin was isolated in Mexico and the United States [7], [8]. Preliminary analysis of the genome of this strain indicated that it is a descendant of common reassortant swine influenza A viruses [9]. Moreover, since 2003, a highly pathogenic H5N1 avian virus has been successfully infecting more than 400 humans with a mortality rate of 60% [10]. It is not clear whether any of these viruses will be the cause of the next human influenza pandemic, however, it is vital to understand the mechanisms behind the genomic evolution of influenza virus and its adaption to new hosts, in particular through the process of reassortment.

Influenza A virus can be found in humans and a variety of animals with aquatic birds being considered as its main reservoir. Influenza viruses do not usually transmit between different hosts. However, pigs are documented to be infected with avian and human viruses, in addition to the swine viruses. Furthermore, multiple reassortment events are found to happen under natural conditions [11]. Hence, it has been postulated that swine are the mixing vessel for inter-host influenza viruses [12].

The influenza A virus genome consists of eight single-stranded RNA segments that code for eleven known proteins. The PB2, PB1, and PA segments encode the RNA polymerase, and HA, NP, NA, and M encode hemagglutinin, nucleoprotein, neuraminidase, and the matrix proteins, respectively. Two distinct non-structural proteins are also coded by the NS segment. The subtypes of influenza A viruses are determined based on their antigenic surface glycoproteins, hemagglutinin and neuraminidase. Hemagglutinin binds to α2,3-galactose- and α2,6-galactose-linked sialic acids. The former is more preferential in avian viruses and the later in human viruses. However they are both present on the tracheal epithelium surface in pigs, making them susceptible to both avian and human viruses.

In addition to the genomic drift of influenza A virus that is caused by the high error rate in the process of replication of its genome, and the antigenic pressure on the HA and NA segments, the evolution of the virus is shaped by the reassortment process. When two different strains of influenza virus co-infect the same cell, new virions can be created that contain a mix of segments from both original strains. This phenomenon was responsible for the 1957 pandemic when the human H1N1 strain that had been circulating since 1918 reassorted to become a human H2N2 strain with new PB1, HA, and NA segments of avian origin [3], [4]. Also, in 1968, the reassortment of the PB1 and HA segments created a new human H3N2 strain which is currently co-circulating with the human H1N1 strain that reappeared in 1977 [4], [13], [14].

Swine classical H1N1 strains have been circulating in pigs since the human influenza pandemic in 1918 and were the dominant strains in the United States until 1998, when two new swine H3N2 strains were identified. These new strains were the results of a double reassortment of swine classical H1N1 with the PB1, HA, and NA segments from a human H3N2 strain, and a triple reassortment of swine classical H1N1, with the PB1, HA, and NA segments of a human H3N2 strain and the PB2 and PA segments of avian lineage [15]–[17]. So far, multiple strains of influenza virus (with various subtypes such as H1N2, H3N1, H2N3, H4N6, H5N1, etc.) have been isolated in pigs around the world, including both inter-host reassortments and whole genome adaptations of human and/or avian viruses [11], [18]–[26].

In this paper, we employ the temporally and geographically diverse information deposited in the Influenza Virus Resource of the National Center for Biotechnology Information [27] to study the reassortment phenomena in swine influenza A viruses. By integrating the information from the publicly available sequences, we investigate patterns in the reassortment events. Applying several statistical techniques, we identify the differential variability of the segments in the influenza genome and enumerate the independent reassortment events. These techniques include diversity/entropy measures of each segment and correlations between them. We find that the reassortment patterns in swine viruses are similar to the previously reported results from human viruses that HA and NA reassort more frequently than the other segments [28], [29]. Surprisingly, we also find that one of the polymerase segments, PB1, reassorts quite frequently, reiterating similar experimental results from human viruses reported by Downie [30].

Methods

To compare the diversity within the segments of swine influenza A virus, we use strains deposited in the Influenza Virus Resource of the NCBI that have all eight segments completely sequenced. We include 150 sequences, containing 99 H1N1, 25 H1N2, 23 H3N2, and 3 H3N1 strains (see Appendix S1). For each segment, we align the sequences of their coding regions using the Smith-Waterman algorithm and calculate the normalized Hamming distances only at the third codon positions, to eliminate the effects of evolutionary pressure due to positive selection. For the M and NS segments, we only consider the coding regions of the M1 and NS1 genes, as they are the longest and the most frequently sequenced sections of the M and NS segments. Because homologous recombination is very rare or absent in influenza viruses [31], [32], this restriction does not alter the results of our analysis.

To measure the diversity of a segment i, we calculate Di, Rao's quadratic entropy [33], according to where N is the total number of strains in the dataset and diab is the Hamming distance between strains a and b at the third codon positions of their corresponding segment i. We estimate the confidence intervals for the diversity measurements via 1000 bootstrap re-samplings of the dataset.

To find the possible reassortant strains, we primarily follow the method introduced by Rabadan et al. (2008), which was initially applied to complete sequences of human influenza A strains [29]. Briefly, in this method, the number of nucleotide differences between the segments of any two strains is calculated. Assuming that the segments have proportional substitution rates at the third codon positions, the differences between two segments of two strains should be proportional if the two segments have a common origin. A violation of this rule indicates that the histories of the two segments are different, i.e. there has been a reassortment event. Therefore, when the distances between two segments of different strains are plotted against each other, the points corresponding to possible reassortment events lie off the diagonal (Figure 1).

[Figure omitted. See PDF.]

Figure 1. Pair-wise hamming distances at the third codon positions in PB2 vs. PB1.

The colors demonstrate the logarithm of the cumulative probability for the points, among which the ones with a cumulative probability of less than 10−7 indicate possible reassortment events. Left: The results from 150 strains in the dataset, where there are candidates for reassortment events in both PB2 and PB1. Right: The results when the dataset is limited to the classical H1N1 strains isolated in the 70's, 80's, and 90's, where there are distinctively more candidates for reassortment events in PB1.

https://doi.org/10.1371/journal.pone.0007366.g001

Given two strains a and b and two segments i and j, the probability to obtain hamming distances equal to diab and djab by random chance only is given by the hypergeometric distribution: where Li and Lj are the respective lengths of the segments divided by three. Hence, fixing the total distance between segments i and j of the two strains, the probability of observing a distance in segment i at most diab is the cumulative of the hypergeometric distribution. Maintaining the assumption of similar average substitution rates at the third codon positions in all segments, in this model the lower the cumulative probability, the more likely it is that the two segments do not have a common ancestor. To correct for multiple hypotheses testing, for every two segments of each strain we generate 100 pairs of segments by randomly permuting their third codon positions. We observe that the cumulative probabilities for distances of pairs from the generated data are at least 10−7. Thus, a cumulative probability of at most 10−7 for two given segments of two strains indicates a reassortment event.

Finally, for each of the 150 strains, we generate a list of strains with which they have low probabilities of having common ancestors, hinting to reassortment events. For further investigation of the origin of the segments, we compile a large target database of more than 10,800 strains of influenza A virus that includes all completely sequenced human and avian isolates, in addition to all swine isolates deposited in the Influenza Virus Resource of the NCBI. We use this database to compare the histories of two segments of a given swine strain. First, we align with NCBI BLAST [34] the two segments to the sequences in the target database, which precede in time the strain of interest. Second, we define the history overlap of the two segments as a function of the alignment identity in the following way. For a given alignment identity x, let Ix be the set of target strains with which the first segment has identity at least x. Similarly define Jx for the second segment. Then the history overlap for alignment identity x is the number of strains common to Ix and Jx over the number of strains included in either one of them. In general, low values of the history overlap function indicate distinct histories of the segments and high values correspond to common history. A decrease in the values of the history overlap function could indicate a potential fork in the lineage of one of the segments. Conversely, an increase can be the result of a merge of the lineages of the two segments, i.e. a reassortment event. Those observations allow us to confirm in an alternative and independent manner the reassortment events predicted by the hypergeometric probability analysis. The history overlap analysis is limited by the sequences present in the target database, but when enough data is available and the converging/diverging lineages are sufficiently different, it can provide a good indicator of the corresponding event.

For a demonstration of the analyses described above consider the strain A/swine/Tennesse/23/1976. When compared to A/swine/Iowa/1/1976 the hamming distance in the PB1 segments is 11% and the hamming distance in the NP segments is 3%. The cumulative hypergeometric probability of this event is less than 10−7, which indicates a reassortment event in at least one of those strains at either segment PB1 or NP. The history overlaps for those two segments and the rest of the segments of A/swine/Tennesse/23/1976 are shown in Figure 2. The figure shows that NP and all the other segments except PB1 share a common recent history, whereas the recent history of PB1 is different from the other seven segments. This allows us to assert that the PB1 segment of A/swine/Tennesse/23/1976 is the result of a reassortment. An interesting feature apparent in Figure 2 is that at lower identities M1 shares fewer strains with PB1 and NP. This observation can be attributed to a possible slower evolutionary rate of M1 and a fork in its lineage to a line of human viruses. Similar considerations show that the PB1 segment of A/swine/Iowa/1/1976 is also the result of a reassortment, however the PB1 segments of A/swine/Iowa/1/1976 and A/swine/Tennessee/23/1976 are from different lineages and the target database contains isolates close to the former, but not the latter.

[Figure omitted. See PDF.]

Figure 2. History overlap for segments PB1 and NP of the swine influenza strain A/Swine/Tennessee/23/1976.

NP and all the other segments except PB1 share a common recent history, whereas the recent history of PB1 is different from the other seven segments, indicating a reassortment event at PB1. The small history overlap of M1 with PB1 and NP at lower identities can be attributed to a possible slower evolutionary rate of M1 and a fork in its lineage to a line of human viruses. The fluctuations in the history overlap of NP at 99% identity are due to small number of sample points at that level of identity.

https://doi.org/10.1371/journal.pone.0007366.g002

Results and Discussion

Viruses present an enormous diversity due to their high mutation rate, short replication time, and high number of replicates. There are several ways of measuring the diversity of a viral population: richness, evenness, Rao's entropy [33], Shannon entropy, other Renyi entropies, etc. When applied to actual viral populations, all these measures encounter similar problems: sampling bias (for instance, most of human influenza samples come from a few studies in New York State and New Zealand [27]), exponential growth and bottleneck structures of viral populations, population stratification, etc. Although the exact interpretation of these measures applied to highly structured populations is not clear, they can be used to compare the variation of diversity in different sections of the genome of a particular organism. Since similar histories imply similar diversity measures, strong differences in these measures for different sections of the genome point to their different histories.

Although the third codon evolutionary rates in influenza A viruses are thought to be similar in all segments, the analysis of the genomic diversity of the strains in our dataset reveals a very inhomogeneous pattern. Figure 3 (left), shows Rao's quadratic entropy, measured at the third codon positions, for 150 swine influenza A viruses that have all 8 fully sequenced segments deposited in NCBI, along with the 95% bootstrap percentile confidence intervals. This figure indicates a statistically significant difference between NA, HA, and PB1 compared to PB2, PA, NP, M, and NS. Moreover, within a particular subtype, where the variations in HA and NA are fixed, PB1 appears as the most diverse segment. Figure 3 (right) shows the diversity in the swine classical H1N1 strains that were isolated in the 70's, 80's and 90's. This analysis shows that the eight segments do not have a common history, with PB1, HA, and NA presenting a higher level of diversity.

[Figure omitted. See PDF.]

Figure 3. Diversity measurements in swine influenza viruses and the corresponding 95% bootstrap percentile confidence intervals.

Left: Considering the 150 strains in the dataset, NA, HA, and PB1 present a higher diversity than the rest. Right: When the dataset is limited to the classical H1N1 strains isolated in the 70's, 80's, and 90's, which fixes the HA and NA variations in the population, shows a higher diversity in PB1 than the rest of the segments.

https://doi.org/10.1371/journal.pone.0007366.g003

We further investigate the sources of variation in diversity via the pair-wise Pearson correlation of the distances at the third codon positions of the viral segments. Correlations, linear or non-linear, or any other measure of dependence, such as mutual information, encounter the same problems as those of the measures of diversity (sampling bias, bottleneck structures, population stratification, etc.). Nonetheless, they are revealing indicators of the origins of diversity in a population. When all the 150 strains in the dataset are considered, the correlations are lower between the surface glycoprotein coding segments and the other segments. More interestingly, the PB1 segment also has a low correlation with all segments that are not polymerase coding (Figure 4, left). Especially, when the strains from a particular subtype are considered and the variations in HA and NA segments are fixed in the dataset, PB1 presents the least correlation relative to the other segments. This is evident, among the classical swine H1N1 strains isolated in the 70's, 80's and 90's (Figure 4, right).

[Figure omitted. See PDF.]

Figure 4. Pair-wise Pearson correlation of the distances at the third codon positions of the viral segments.

Left: The HA, NA, and PB1 segments have the least correlation in regards to the rest of the segments. Right: When the HA and NA variations are fixed in the population by limiting the dataset to the classical H1N1 strains isolated in the 70's, 80's, and 90's, PB1 presents a distinctively lower correlation relative to the other segments.

https://doi.org/10.1371/journal.pone.0007366.g004

The above observations from the diversity and pair-wise correlation measures hint to distinct evolutionary patterns in the HA, NA, and PB1 segments. To elucidate the role of the process of reassortment in these patterns, we have enumerated the independent reassortment events in swine viruses that we identify through the hypergeometric distribution analysis of Rabadan et al. (2008) [29] and confirm via history overlap analysis, described in the Methods section. Table 1 lists these events, represented by different strains and a simple inspection reveals the frequent role of NA, HA, and PB1 in the reassortment process. Because the reassortment events are frequent in swine viruses and the sampling is not, it is difficult to determine their exact history. However, especially in cases where there are multiple reassortments, we have attempted to identify the fully sequenced strains that are the earliest independent isolations of the reassortant viruses. As indicated in the last column of Table 1, some of the listed strains have been already published independently. In addition, for a more comprehensive list, we have included the reported reassortment events for which there are no completely sequenced isolates available, so that they could not be identified by our method. Finally, we have listed other published reassortant strains, which according to Krasnitz et al. (2008) either are “frozen in time” or show evidence of homologous recombination [31]. In addition to independent confirmation of the known cases of reassortment in swine viruses, our methods have succeeded in more than doubling the number of cases, as shown in Table 1.

[Figure omitted. See PDF.]

Table 1. Summery of the reassortment events in swine influenza A viruses.

https://doi.org/10.1371/journal.pone.0007366.t001

To summarize, our analyses show that not every segment of the swine influenza virus reassorts in equal fashion. In accordance with the previous results from human influenza A viruses, both in vitro [28] and in vivo [29], we find that the surface glycoproteins coding segments (HA and NA) of swine influenza viruses reassort at a higher rate as well. Perhaps, the most intriguing conclusion of our analyses is the characteristic role of one of the polymerase coding segments (PB1) that appears frequently in both inter-host and intra-host reassortment events among swine viruses and we believe that our overall analysis is the first to quantify this role. Interestingly, this is the same pattern observed in the strains responsible for the 1957 and 1968 pandemics, when human viruses also obtained PB1 segments of avian origin.

The mechanisms behind the preferential reassortments are not clear, however several hypotheses can be advanced. There is substantial evidence for biases in the packaging mechanism of the viral RNA into the virion for influenza A viruses, which can impose a selective pressure on segments that can be exchanged between strains [35]. Another constraint on the reassortment events can be associated with compensatory mutations due to interactions between the different proteins.

Supporting Information

[Figure omitted. See PDF.]

Appendix S1.

https://doi.org/10.1371/journal.pone.0007366.s001

(0.04 MB DOC)

Author Contributions

Conceived and designed the experiments: HK VT RR. Performed the experiments: HK VT RR. Analyzed the data: HK VT RR. Contributed reagents/materials/analysis tools: HK VT RR. Wrote the paper: HK VT RR.

Citation: Khiabanian H, Trifonov V, Rabadan R (2009) Reassortment Patterns in Swine Influenza Viruses. PLoS ONE4(10): e7366. https://doi.org/10.1371/journal.pone.0007366

References

1. Taubenberger JK, Reid AH, Lourens RM, Wang R, Jin G, et al. (2005) Characterization of the 1918 influenza virus polymerase genes. Nature 437(7060): 889–93.JK TaubenbergerAH ReidRM LourensR. WangG. Jin2005Characterization of the 1918 influenza virus polymerase genes.Nature437(7060)88993

2. Rabadan R, Levine AJ, Robins H (2006) Comparison of avian and human influenza A viruses reveals a mutational bias on the viral genomes. J Virol 80(23): 11887–91.R. RabadanAJ LevineH. Robins2006Comparison of avian and human influenza A viruses reveals a mutational bias on the viral genomes.J Virol80(23)1188791

3. Lindstrom SE, Cox NJ, Klimov A (2004) Genetic analysis of human H2N2 and early H3N2 influenza viruses, 1957–1972: evidence for genetic divergence and multiple reassortment events. Virology 328(1): 101–19.SE LindstromNJ CoxA. Klimov2004Genetic analysis of human H2N2 and early H3N2 influenza viruses, 1957–1972: evidence for genetic divergence and multiple reassortment events.Virology328(1)10119

4. Scholtissek C, Rohde W, Von Hoyningen V, Rott R (1978) On the origin of the human influenza virus subtypes H2N2 and H3N2. Virology 87(1): 13–20.C. ScholtissekW. RohdeV. Von HoyningenR. Rott1978On the origin of the human influenza virus subtypes H2N2 and H3N2.Virology87(1)1320

5. Myers KP, Olsen CW, Gray GC (2007) Cases of Swine Influenza in Humans: A Review of the Literature. Clin Infect Dis 44(8): 1084–1088.KP MyersCW OlsenGC Gray2007Cases of Swine Influenza in Humans: A Review of the Literature.Clin Infect Dis44(8)10841088

6. Shinde V, Bridges CB, Uyeki TM, Shu B, Balish A, et al. (2009) N Engl J Med 360(25): 2616–25.V. ShindeCB BridgesTM UyekiB. ShuA. Balish2009N Engl J Med360(25)261625

7. Center for Disease Control (2009) Swine Influenza A (H1N1) Infection in Two Children - Southern California, March - April 2009. MMWR 58(15): 400–402.Center for Disease Control2009Swine Influenza A (H1N1) Infection in Two Children - Southern California, March - April 2009.MMWR58(15)400402

8. Center for Disease Control (2009) Outbreak of Swine-Origin Influenza A (H1N1) Virus Infection - Mexico, March–April 2009. MMWR 58(Dispatch): 1–3.Center for Disease Control2009Outbreak of Swine-Origin Influenza A (H1N1) Virus Infection - Mexico, March–April 2009.MMWR58(Dispatch)13

9. Trifonov V, Khiabanian H, Greenbaum B, Rabadan R (2009) The origin of the recent swine influenza A(H1N1) virus infecting humans. Euro Surveill 14(17): pii = 19193.V. TrifonovH. KhiabanianB. GreenbaumR. Rabadan2009The origin of the recent swine influenza A(H1N1) virus infecting humans.Euro Surveill14(17)pii = 19193

10. Komar N, Olsen B (2008) Avian Influenza Virus (H5N1) Mortality Surveillance. Emerg Infect Dis 14(7): 1176–1178.N. KomarB. Olsen2008Avian Influenza Virus (H5N1) Mortality Surveillance.Emerg Infect Dis14(7)11761178

11. Brown IA (2000) The epidemiology and evolution of influenza viruses in pigs. Vet Microbiol 74(1–2): 29–46.IA Brown2000The epidemiology and evolution of influenza viruses in pigs.Vet Microbiol74(1–2)2946

12. Scholtissek C (1990) Pigs as the ‘mixing vessel’ for the creation of new pandemic influenza A viruses. Med Princip Prac 2: 65–71.C. Scholtissek1990Pigs as the ‘mixing vessel’ for the creation of new pandemic influenza A viruses.Med Princip Prac26571

13. Nakajima K, Desselberger U, Palese P (1978) Recent human influenza A (H1N1) viruses are closely related genetically to strains isolated in 1950. Nature 274(5669): 334–339.K. NakajimaU. DesselbergerP. Palese1978Recent human influenza A (H1N1) viruses are closely related genetically to strains isolated in 1950.Nature274(5669)334339

14. Scholtissek C, von Hoyningen V, Rott R (1978) Genetic relatedness between the new 1977 epidemic strains (H1N1) of influenza and human influenza strains isolated between 1947 and 1957 (H1N1). Virology 89(2): 613–617.C. ScholtissekV. von HoyningenR. Rott1978Genetic relatedness between the new 1977 epidemic strains (H1N1) of influenza and human influenza strains isolated between 1947 and 1957 (H1N1).Virology89(2)613617

15. Webby RJ, Swenson SL, Krauss SL, Gerrish PJ, Goyal SM, et al. (2000) Evolution of swine H3N2 influenza viruses in the United States. J Virol 74(18): 8243–51.RJ WebbySL SwensonSL KraussPJ GerrishSM Goyal2000Evolution of swine H3N2 influenza viruses in the United States.J Virol74(18)824351

16. Zhou NN, Senne DA, Landgraf JS, Swenson SL, Erickson G, et al. (1999) Genetic reassortment of avian, swine, and human influenza A viruses in American pigs. J Virol 73(10): 8851–6.NN ZhouDA SenneJS LandgrafSL SwensonG. Erickson1999Genetic reassortment of avian, swine, and human influenza A viruses in American pigs.J Virol73(10)88516

17. Vincent AL, Ma W, Lager KM, Janke BH, Richt JA (2008) Swine influenza viruses: a North American perspective. Adv Virus Res 72: 127–54.AL VincentW. MaKM LagerBH JankeJA Richt2008Swine influenza viruses: a North American perspective.Adv Virus Res7212754

18. Marozin S, Gregory V, Cameron K, Bennett M, Valette M, et al. (2002) Antigenic and genetic diversity among swine influenza A H1N1 and H1N2 viruses in Europe. J Gen Virol 83(Pt 4): 735–745.S. MarozinV. GregoryK. CameronM. BennettM. Valette2002Antigenic and genetic diversity among swine influenza A H1N1 and H1N2 viruses in Europe.J Gen Virol83(Pt 4)735745

19. Qi X, Lu CP (2006) Genetic characterization of novel reassortant H1N2 influenza A viruses isolated from pigs in southeastern. Arch Virol 151(11): 2289–99.X. QiCP Lu2006Genetic characterization of novel reassortant H1N2 influenza A viruses isolated from pigs in southeastern.Arch Virol151(11)228999

20. Zell R, Motzke S, Krumbholz A, Wutzler P, Herwig V, et al. (2008) Novel reassortant of swine influenza H1N2 virus in Germany. J Gen Virol 89(PT 1): 271–276.R. ZellS. MotzkeA. KrumbholzP. WutzlerV. Herwig2008Novel reassortant of swine influenza H1N2 virus in Germany.J Gen Virol89(PT 1)271276

21. Saito T, Suzuki H, Maeda K, Inai K, Takemae N, et al. (2006) Molecular characterization of an H1N2 swine influenza virus isolated in Miyazaki, Japan, in 2006. J Vet Med Sci 70(4): 423–7.T. SaitoH. SuzukiK. MaedaK. InaiN. Takemae2006Molecular characterization of an H1N2 swine influenza virus isolated in Miyazaki, Japan, in 2006.J Vet Med Sci70(4)4237

22. Ludwig S, Stitz L, Planz O, Van H, Fitch WM, et al. (1995) European swine virus as a possible source for the next influenza pandemic? Virology 212(2): 555–61.S. LudwigL. StitzO. PlanzH. VanWM Fitch1995European swine virus as a possible source for the next influenza pandemic?Virology212(2)55561

23. Castrucci MR, Donatelli I, Sidoli L, Barigazzi G, Kawaoka Y, et al. (1993) Genetic reassortment between avian and human influenza A viruses in Italian pigs. Virology 193(1): 503–6.MR CastrucciI. DonatelliL. SidoliG. BarigazziY. Kawaoka1993Genetic reassortment between avian and human influenza A viruses in Italian pigs.Virology193(1)5036

24. Chutinimitkul S, Thippamom N, Damrongwatanapokin S, Payungporn S, Thanawongnuwech R, et al. (2008) Genetic characterization of H1N1, H1N2 and H3N2 swine influenza virus in Thailand. Arch Virol 153(6): 1049–56.S. ChutinimitkulN. ThippamomS. DamrongwatanapokinS. PayungpornR. Thanawongnuwech2008Genetic characterization of H1N1, H1N2 and H3N2 swine influenza virus in Thailand.Arch Virol153(6)104956

25. Yoo GJ, Kim CJ, Choi YK (2008) Seroprevalence and genetic evolutions of swine influenza viruses under vaccination pressure in Korean swine herds. Virus Res 138 (1–2): 43–49.GJ YooCJ KimYK Choi2008Seroprevalence and genetic evolutions of swine influenza viruses under vaccination pressure in Korean swine herds.Virus Res 138(1–2)4349

26. Karasin AI, Brown IH, Carman S, Olsen CW (2000) Isolation and characterization of H4N6 avian influenza viruses from pigs with pneumonia in Canada. J Virol 74(19): 9322–7.AI KarasinIH BrownS. CarmanCW Olsen2000Isolation and characterization of H4N6 avian influenza viruses from pigs with pneumonia in Canada.J Virol74(19)93227

27. Bao Y, Bolotov P, Dernovoy D, Kiryutin B, Zaslavsky L, et al. (2008) The Influenza Virus Resource at the National Center for Biotechnology Information. J Virol 82(2): 596–601.Y. BaoP. BolotovD. DernovoyB. KiryutinL. Zaslavsky2008The Influenza Virus Resource at the National Center for Biotechnology Information.J Virol82(2)596601

28. Lubeck MD, Palese P, Schulman JL (1979) Nonrandom association of parental genes in influenza A virus recombinants. Virology 95(1): 269–74.MD LubeckP. PaleseJL Schulman1979Nonrandom association of parental genes in influenza A virus recombinants.Virology95(1)26974

29. Rabadan R, Levin AJ, Krasnitz M (2008) Non-random reassortment in human influenza A viruses. Influenza and Other Respiratory Viruses 2(1): 9–22.R. RabadanAJ LevinM. Krasnitz2008Non-random reassortment in human influenza A viruses.Influenza and Other Respiratory Viruses2(1)922

30. Downie JC (2004) Reassortment of influenza A virus genes linked to PB1 polymerase gene. International Congress Series 1263: 714–718.JC Downie2004Reassortment of influenza A virus genes linked to PB1 polymerase gene.International Congress Series1263714718

31. Krasnitz M, Levine AJ, Rabadan R (2008) Anomalies in the Influenza Virus Genome Database: New Biology or Laboratory Errors? J Virol 82(17): 8947–8950.M. KrasnitzAJ LevineR. Rabadan2008Anomalies in the Influenza Virus Genome Database: New Biology or Laboratory Errors?J Virol82(17)89478950

32. Boni MF, Zhou Y, Taubenberger JK, Holmes EC (2008) Homologous Recombination Is Very Rare or Absent in Human Influenza A Virus. J Virol 82(10): 4807–4811.MF BoniY. ZhouJK TaubenbergerEC Holmes2008Homologous Recombination Is Very Rare or Absent in Human Influenza A Virus.J Virol82(10)48074811

33. Rao CR (1982) Diversity and dissimilarity coefficients: a unified approach. Theor Popul Biol 21(1): 24–43.CR Rao1982Diversity and dissimilarity coefficients: a unified approach.Theor Popul Biol21(1)2443

34. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402.SF AltschulTL MaddenAA SchafferJ. ZhangZ. Zhang1997Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.Nucleic Acids Res2533893402

35. Marsh GA, Rabadan R, Levine AJ, Palese P (2008) Highly conserved regions of influenza a virus polymerase gene segments are critical for efficient viral RNA packaging. J Virol 82(5): 2295–304.GA MarshR. RabadanAJ LevineP. Palese2008Highly conserved regions of influenza a virus polymerase gene segments are critical for efficient viral RNA packaging.J Virol82(5)2295304

36. Karasin AI, Schutten MM, Cooper LA, Smith CB, Subbarao K, et al. (2000) Genetic characterization of H3N2 influenza viruses isolated from pigs in North America, 1977–1999: evidence for wholly human and reassortant virus genotypes. Virus Res 68(1): 71–85.AI KarasinMM SchuttenLA CooperCB SmithK. Subbarao2000Genetic characterization of H3N2 influenza viruses isolated from pigs in North America, 1977–1999: evidence for wholly human and reassortant virus genotypes.Virus Res68(1)7185

37. Karasin AI, Landgraf J, Swenson S, Erickson G, Goyal S, et al. (2002) Genetic characterization of H1N2 influenza A viruses isolated from pigs throughout the United States. J Clin Microbiol 40(3): 1073–1079.AI KarasinJ. LandgrafS. SwensonG. EricksonS. Goyal2002Genetic characterization of H1N2 influenza A viruses isolated from pigs throughout the United States.J Clin Microbiol40(3)10731079

38. Lekcharoensuk P, Lager KM, Vemulapalli R, Woodruff M, Vincent AL, et al. (2006) Novel swine influenza virus subtype H3N1, United States. Emerg Infect Dis 12(5): 787–94.P. LekcharoensukKM LagerR. VemulapalliM. WoodruffAL Vincent2006Novel swine influenza virus subtype H3N1, United States.Emerg Infect Dis12(5)78794

39. Ma W, Vincent AL, Gramer MR, Brockwell CB, Lager KM, et al. (2007) Identification of H2N3 influenza A viruses from swine in the United States. Proc Natl Acad Sci U S A 104(52): 20949–20954.W. MaAL VincentMR GramerCB BrockwellKM Lager2007Identification of H2N3 influenza A viruses from swine in the United States.Proc Natl Acad Sci U S A104(52)2094920954

40. Karasin AI, Carman S, Olsen CW (2006) Identification of human H1N2 and human-swine reassortant H1N2 and H1N1 influenza A viruses among pigs in Ontario, Canada (2003 to 2005). J Clin Microbiol 44(3): 1123–6.AI KarasinS. CarmanCW Olsen2006Identification of human H1N2 and human-swine reassortant H1N2 and H1N1 influenza A viruses among pigs in Ontario, Canada (2003 to 2005).J Clin Microbiol44(3)11236

41. Olsen CW, Karasin AI, Carman S, Li Y, Bastien N, et al. (2006) Triple reassortant H3N2 influenza A viruses, Canada. Emerg Infect Dis 12 (7): 1132–1135.CW OlsenAI KarasinS. CarmanY. LiN. Bastien2006Triple reassortant H3N2 influenza A viruses, Canada.Emerg Infect Dis 12(7)11321135