An efficient genome sequencing method for equine influenza [H3N8] virus reveals a new polymorphism in the PA-X protein
© Rash et al.; licensee BioMed Central Ltd. 2014
Received: 4 June 2014
Accepted: 20 August 2014
Published: 2 September 2014
H3N8 equine influenza virus (EIV) has caused disease outbreaks in horses across the world since its first isolation in 1963. However, unlike human, swine and avian influenza, there is relatively little sequence data available for this virus. The majority of published sequences are for the segment encoding haemagglutinin (HA), one of the two surface glycoproteins, making it difficult to study the evolution of the other gene segments and determine the level of reassortment occurring between sub-lineages.
To facilitate the generation of full genome sequences for EIV, we developed a simple, cost-effective and efficient method. M13-tagged primers were used to amplify short, overlapping RT-PCR products, which were then sequenced using Sanger dideoxynucleotide sequencing technology. We also modified a previously published method, developed for human H3N2 and avian H5N1 influenza viruses, which was based on the ligation of viral RNA and subsequent amplification by RT-PCR, to sequence the non-coding termini (NCRs). This necessitated the design of novel primers for an N8 neuraminidase segment.
Two field isolates were sequenced successfully, A/equine/Lincolnshire/1/07 and A/equine/Richmond/1/07, representative of the Florida sublineage clades 1 and 2 respectively. A total of 26 PCR products varying in length from 400–600 nucleotides allowed full coverage of the coding sequences of the eight segments, with sufficient overlap to allow sequence assembly with no primer-derived sequences. Sequences were also determined for the non-coding regions and revealed cytosine at nucleotide 4 in the polymerase segments. Analysis of EIV genomes sequenced using these methods revealed a novel polymorphism in the PA-X protein in some isolates.
These methods can be used to determine the genome sequences of EIV, including the NCRs, from both clade 1 and clade 2 of the Florida sublineage. Full genomes were covered efficiently using fewer PCR products than previously reported methods for influenza A viruses, the techniques used are affordable and the equipment required is available in most research laboratories. The adoption of these methods will hopefully allow for an increase in the number of full genomes available for EIV, leading to improved surveillance and a better understanding of EIV evolution.
KeywordsEquine influenza virus H3N8 Genome sequencing Non-coding regions M13 PA-X
Equine influenza virus (EIV) is an influenza A virus belonging to the Orthomyxoviridae family. These viruses have a negative sense, single-stranded RNA genome consisting of eight viral gene segments . Originally thought to have transmitted from birds, H3N8 EIV was first isolated during a widespread outbreak in the United States in 1963 , and has since spread worldwide causing multiple major outbreaks of disease in horses. During the 1980s the virus diverged into two antigenically distinct lineages , American and Eurasian, and since then the American lineage has further evolved into the Florida sublineage clades 1 and 2, which continue to co-circulate today . These lineages have historically been based on antigenic and genetic data for HA. A phylogenetic study by Murcia et al.  showed that phylogenetic trees produced for each of the viral gene segments also supported division into the American and Eurasian lineages, and all but segment 7 divided into the two clades of the Florida sublineage. However, less than 100 complete viral genomes covering the 46 years from 1963 to 2008 were available at the time of the study.
Other groups have studied the evolution of individual influenza A virus genes. H3N8 EIV PB2  and matrix proteins  were found to belong to the same lineage as North American avian strains, whilst PB1, PA, HA and NP were found to have evolved independently from other influenza A viruses [8–10]. Equine NS was suggested as being restricted to subtype , as the NS segments of the H3N8 viruses were close to one another but not to that of the H7N7 viruses, however very limited numbers of EIV genomes were available at the time these studies were performed, and in some cases only two different EIV strains were used. More recently a study found that the internal genes of the 1963 EIV pandemic virus were of western hemispheric avian influenza origin . This study also showed that the virus shared a most recent common ancestor with avian influenza viruses from South America shortly before its emergence. The time to most recent common ancestor for avian/equine NP was calculated as being 1954, which agreed with the hypothesis that the virus emerged in South America prior to its introduction into the USA in 1963 by horses imported by air from Argentina .
Following an extensive outbreak in 1989 affecting a highly vaccinated population of racehorses in the UK, it became clear that, like human influenza virus, EIV undergoes antigenic drift and therefore vaccine strains need to be kept up to date [13, 14]. A formal process for vaccine strain selection, overseen by the World Organisation for Animal Health (OIE) was put in place. This process relies heavily on surveillance data collected from the field, of which most is focussed solely on the HA gene and the protein it encodes. Therefore the majority of published sequences are for HA, which makes it difficult to study the evolution of either the other gene segments or the virus as a whole. It is also known that reassortment between the different lineages of EIV has occurred [4–8, 15, 16] but the full extent is unknown due to a lack of data available. Avian H3N8 influenza viruses have been shown to frequently exchange internal gene segments, and it has been suggested that the extensive reassortment within the H3 subtype poses a threat to human and animal health . Next generation sequencing technologies have made it easier to sequence whole viral genomes, however these technologies are not readily available to all as considerable investment in equipment and bioinformatics expertise are needed. We aimed to develop a simple and robust method to sequence whole EIV genomes from all H3N8 lineages using Sanger dideoxynucleotide sequencing technology.
Each of the eight influenza virus gene segments contains two non-coding regions (NCRs), one at the 5′ terminus containing 13 conserved nucleotides, and the other at the 3′ terminus, which contains 12 nucleotides . Unlike the 5′ end, the 3′ terminus exhibits variation at the fourth nucleotide. This variation in the fourth nucleotide has been shown to affect the rescue of virus from a reverse genetics system . A second objective, to implement a method previously described for sequencing the NCRs of influenza viruses , was adapted and carried out on an EIV, as well as an N8 subtype neuraminidase, for the first time.
Here we describe the genome sequencing method and highlight the sequence differences found between representatives of the two circulating clades of the Florida sublineage.
Genome sequencing of equine influenza viruses
Primer sequences and annealing temperatures used to sequence the genome of H3N8 EIV
Primer sequence (5′-3′)
Approximate nucleotide coverage (5′-3′)
Annealing temperature used (°C)
GC GTAAAACGACGGCCAGT AGCGAAAGCAGGCAAATATATTCAATATG
GC AACAGCTATGACCATG CTCTTTCTAGCATGTAT
GC GTAAAACGACGGCCAGT CACAACTAACAATAACCAA
GC AACAGCTATGACCATG CCTCAAGAGTTGATG
GC GTAAAACGACGGCCAGT GCAATAATTGTAGCC
GC AACAGCTATGACCATG ATTATTTGAGCAGTATC
GC GTAAAACGACGGCCAGT GAAGCCAATACAGCGGT
GC AACAGCTATGACCATG AGTAGAAACAAGG TCGTTTTTAAACAATTC
GC GTAAAACGACGGCCAGT AGCGAAAGCAGG CAAACCATTTGAATGG
GC AACAGCTATGACCATG CAGCGTCCTTGGTCATTG
GC GTAAAACGACGGCCAGT CTTCCAACGAAAGAGAA
GC AACAGCTATGACCATG GGTTTAATATGGATACACC
GC GTAAAACGACGGCCAGT GCGGCTTCACTGAGTCCTGGC
GC AACAGCTATGACCATG CATTTTAAACAAACTTC
GC GTAAAACGACGGCCAGT CAAAGACTGGTCTACTG
GC AACAGCTATGACCATG AGTAGAAACAAGG CATTTTTTCATGAAGATC
GC GTAAAACGACGGCCAGT AGCAAAAGCAGG TACTGATCCAAAATGG
GC AACAGCTATGACCATG GCCTCTCTCGGACTGAC
GC GTAAAACGACGGCCAGT GCCAGAATCAAGACCAGG
GC AACAGCTATGACCATG CTCACTTGGAATCCAACTTGC
GC GTAAAACGACGGCCAGT GAGAGAAAGTGGATTTTGAGGATTG
GC AACAGCTATGACCATG CTGAAGGAGGCAGCGCC
GC GTAAAACGACGGCCAGT GACCCATGTTTTTGTATG
GC AACAGCTATGACCATG AGTAGAAACAAGG TACTTTTTTGGACAG
GC GTAAAACGACGGCCAGT AGCGAAAGCAGGGGACGATATT
GC AACAGCTATGACCATG GATTTGTTAGCCAATTCAG
GC GAAAACGACGGCCAGT CAGGTGTCACTCAAAAC G
GC AACAGCTATGACCATG GGATTTGCTTTTCTGGTAC
GC GTAAAACGACGGCCAGT GGTTACATATGGAAAATGCC
GC AACAGCTATGACCATG GAGCCACCAGCAATTCT
GC GTAAAACGACGGCCAGT GAAGGAAGAATTCAGGA
GC AACAGCTATGACCATG GAGTAGAAACAAGGGTGTTTTTAAC
GC GTAAAACGACGGCCAGT AGCGAAAGCAGGGTAGATAATC
GC AACAGCTATGACCATG CCGTGGGAGGGTTGAGCC
GC GTAAAACGACGGCCAGT GACACCACATACCAAAC
GC AACAGCTATGACCATG CTCTCAGGTCCTCAAAT
GC GTAAAACGACGGCCAGT CCAGCACACAAGAGCCAG
GC AACAGCTATGACCATG AGTAGAAACAAGGGTATTTTTC
GC GTAAAACGACGGCCAGT AGCAAAAGCAGGAGTTT
GC AACAGCTATGACCATG GCCCTATTTTGACACTC
GC GTAAAACGACGGCCAGT CACACAGGGCTCATTAC
GC AACAGCTATGACCATG CCGAAACCTTTTACACCG
GC GTAAAACGACGGCCAGT CACAGTTGGATATTTGTG
GC AACAGCTATGACCATG AGTAGAAACAAGGAGTT
GC GTAAAACGACGGCCAGT AGCGAAAGCAGGTAGATATTTAAAG
GC AACAGCTATGACCATG CTAGCCTTACTAGCAAC
GC GTAAAACGACGGCCAGT CAGTACCACGGCTAAAG
GC AACAGCTATGACCATG AGTAGAAACAAGGTAGTTTTTTAC
GC GTAAAACGACGGCCAGT AGCGAAAGCAGGGTGACAAAAAC
GC AACAGCTATGACCATG CTGCTCCTTCTTCGGTG
GC GTAAAACGACGGCCAGT CATCATACTTAAAGCAAAC
GC AACAGCTATGACCATG AGTAGAAACAAGGTAGTGTTTTTTAT
Nucleotide and amino acid differences between A/equine/Richmond/1/07 and A/equine/Lincolnshire/1/07
Amino acid changes
Amino acid changes from Richmond/1/07 to Lincolnshire/1/07
A105T, K251R, I398V, K660R
F94L, M119V, V149I, M179I, R329Q, E377D,D618E, K621R
D64E, I86M, M210A, K237E, G240E, P259S, N321S, L348I, S409N, I465V, T476A, I500V, R626K
K-14T, F-11L, I-9a1, F-9b1, N7G, R62K, V78A, D104N, A138S, N159S, E291D
T43A, E50G, I198V
I131M, T257I, A359T, S450N
T9A, S12F, V35A, E40K, G42D, H66Y, P78S, I191V, N235D, S337N, I410V, G416E
I15V, I80V, K95R
I48S, I84V, Y207H, G210R
The two smaller segments, M and NS, contained a total of 10 amino acid changes within the four predicted polypeptides (M1, M2, NS1 and NS2/NEP) that they encode.
Nucleotide changes in three segments resulted in different lengths for their predicted polypeptides. A duplication of six nucleotides in A/equine/Richmond/1/07 resulted in a two amino acid insertion within the putative signal peptide of the precursor HA protein extending its length from 15 to 17 amino acids, as observed in recent FC2 isolates . In contrast, when compared to earlier EIV isolates from 1963–2000 the predicted amino acid sequences for the NS1 protein from both viruses were truncated by 11 amino acids, as seen for other recent isolates . This was caused by a premature stop codon at position 220, resulting in a predicted polypeptide length of 219 amino acids. The open reading frame of NS2/NEP however, was unaffected by this nucleotide substitution. In addition, a novel truncation in the recently discovered PA-X gene was identified in A/equine/Richmond/1/07, caused by an early stop codon at position 20 of the +1 reading frame. The truncation of PA-X by 42 amino acids has not been described before, with the majority of strains having either a full length version (252 amino acids) or are truncated by 19 amino acids. To investigate further and to study the evolution of the truncation, the PA-X region of segment 3 from an additional 29 EIV isolated in the UK between 2005 and 2013, including 9 from 2007, were sequenced using the method described here (see Additional file 3). The 42 amino acid truncation in A/equine/Richmond/1/07 PA-X was identified in three of the isolates, one from another horse in the same outbreak that A/equine/Richmond/1/07 was isolated from, and the other two from a separate outbreak in 2007. The remaining isolates did not share the truncation, and all had a full length PA-X of 252 amino acids (61 amino acids following the +1 frameshift) (see Additional file 3).
Sequencing of the non-coding regions of equine influenza viruses
Primer sequences used for sequencing 3′ and 5′ NCRs of H3N8 EIV
5′-GG G TATTTCATTGCCATCATCC-3′
5′-GAC A GT A TCCATGGTGTATCCTGT-3′
5′-G T ATGGT T GAGGCCATGGTGTC C-3′
5′-A T C CCTGTGATTTCAAATCTTTCTTC-3′
5′-GA A TGCCTGATTAATGATCCCTG-3′
5′-A AT GTTCC A TTTG CT ACTGCATG-3′
5′-GGATTTC A TT C GCCATATCATG-3′
5′-TCAAA T GC C GA A AG T-3′
5′-CCGATCGTGCC T TC C TTTGACA T-3′
5′-A T AAAGCGTCTACGCTGCAGTCC-3′
5′-C C GT A TTATCATTCCATT T AAG-3′
5′-TGATAA T ACG GTTA GAA TCTCT-3′
We and others have previously shown that reassortment has occurred between different EIV [4–8, 15, 16], but a lack of full genome sequences for EIV makes it difficult to ascertain the extent of reassortment between them, and whether reassortment has occurred between EIV and influenza A viruses from other species. We therefore developed a simple method for sequencing viral genomes that included primers with M13 sequence tags to improve the sequencing efficiency. This was based upon the approach recommended by the WHO for sequencing swine influenza virus isolates in 2009 , however our method used only 26 PCR fragments to cover the segment-specific regions of EIV, rather than 46 fragments. Alternative methods have been employed for sequencing influenza A viruses, such as using universal primers to simultaneously amplify all eight genome segments, or segment specific primers to amplify entire segments; however, in our hands such protocols result in poor amplification of the three largest genome segments (data not shown). Other methods based on amplification of small PCR fragments do not include the M13 sequences in the primers, which makes the method described here simple and efficient. A method previously described for sequencing the NCRs of the influenza gene segments was also modified and successfully used for the first time on an equine influenza virus, with novel primers designed for an N8 subtype NA.
The two segments encoding the surface glycoproteins, HA and NA, contained a large number of amino acid differences between the two viruses. This was expected as these two proteins are under constant immune-driven selection pressure to undergo antigenic drift. Interestingly a high number of amino acid differences were found in PA, especially when compared to the other two polymerase subunits PB2 and PB1, and a similar finding was observed by Murcia et al. . The other internal segments contained fewer changes, which is not surprising as they are both smaller and may be under less immune pressure than the surface proteins.
Interestingly, a mutation in the +1 reading frame of PA, causing a premature stop codon in the translated amino acid sequence of PA-X, was observed in A/equine/Richmond/1/07. PA-X is a recently discovered protein containing the N-terminal 191 amino acids of PA and, in the majority of strains, a further 61 amino acids derived from a frameshift to the +1 reading frame of PA . PA-X has been implicated in the modulation of influenza virus pathogenicity and virulence in a mouse model, whereby PA-X deficient viruses caused greater clinical signs and were less able to shut off host cell responses compared to wild-type viruses with full length PA-X . The premature stop codon in A/equine/Richmond/1/07 would lead to a truncation of the protein by 42 amino acids. Truncated forms of PA-X have been described previously, however the majority of these are due to a nonsense mutation at codon 42 in the +1 reading frame . Sequencing of PA, as described here, revealed that several other virus isolates from different outbreaks in 2007 as well as from the same yard as A/equine/Richmond/1/07, had the same truncated form of PA-X, however the truncated form did not persist in the UK.
Sequence analysis of the NCRs from each segment showed that EIV strain A/equine/Richmond 1/07 had cytosine at position 4 of the 3′ vRNA in the three polymerase segments, as found in other influenza viruses, and uracil at this position in the remaining 5 segments. This is the same pattern seen in the majority of other influenza A viruses for which the promoter sequences have been determined, including the prototype avian influenza virus, A/chicken/Rostock/34 (H7N1) .
The methods outlined here can be used to determine the genome sequences of EIV, including the NCRs, from both clade 1 and clade 2 of the Florida sublineage. The techniques described here are affordable, and the equipment required is available in most research laboratories. The sequence assembly process is simple and does not require in depth bioinformatics, unlike next generation sequencing methodology. Due to the small genome size and small sample numbers usually associated with EIV, this method is therefore highly cost effective and straightforward. Amplicon sequencing has also been shown to be less labour intensive and more affordable than plasmid cloning methods . This method also permits the sequencing of individual gene segments with relative ease, as was the case with PA described here to investigate the frequency of the truncated form of PA-X.
We have developed a simple, efficient and affordable method for sequencing whole genomes of EIV that offers an improvement compared with previously published methods. The adoption of these methods should facilitate an increase in the number of full genome sequences available for EIV. This will benefit surveillance programmes for EIV and improve understanding of the evolutionary paths taken by the virus, including the level of reassortment.
EIV A/equine/Richmond/1/07 and A/equine/Lincolnshire/1/07 had previously been isolated and passaged twice in embryonated chicken eggs . RNA was isolated from 140 μl virus stocks containing ~107 EID50/ml using a QIAamp Viral RNA Mini Kit (Qiagen), according to manufacturer’s directions. RNA was eluted in 50 μl elution buffer.
Using a UNI-12 primer, as described by Hoffmann et al. , cDNA was transcribed by denaturing 2 μl RNA in the presence of UNI-12 (1 μM final concentration) and 7 μl water at 70°C for 10 minutes and then cooling on ice. Following this, dNTP mix (final concentration each 0.5 μM), 1 × First Strand buffer, 200U Superscript II reverse transcriptase (Invitrogen) and water to a final volume of 20 μl were added. The reaction mixture was incubated at 42°C for 45 minutes.
PCR amplification of gene segments
Viral gene segments were amplified in 50 μl PCR reactions consisting of 2 μl cDNA (representing 10% of the reverse transcription reaction), dNTP mix (0.2 mM each final concentration) (Qiagen), 1 × Pfu buffer, 2.5U Native Pfu DNA polymerase (Stratagene), water and oligonucleotide pairs (final concentration each of 0.2 μM) as listed in Table 1. The cycling conditions were as follows: initial denaturation at 96°C for 1 minute, followed by 25 cycles of denaturation at 96°C for 15 seconds, primer annealing at 50-60°C (see Table 1) for 10 seconds and elongation at 60°C for 5 minutes. PCR reactions were analysed on a 1% agarose gel containing GelRed nucleic acid stain (Biotium) according to manufacturer’s directions. PCR products were purified using a QIAquick PCR purification kit (Qiagen) according to manufacturer’s directions.
PCR amplification of non-coding regions
The method described by de Wit et al.  was used with modifications, as detailed in Table 2. Novel primers were designed to amplify the N8 subtype NA segment. Briefly, following an initial denaturation at 65°C for 5 minutes in the presence of T4 RNA ligase buffer and 20U RNAsin RNase inhibitor (Promega), 15 μl RNA was ligated using 40U T4 RNA ligase (New England Biolabs) at 37°C for 1 hour. The ligation reaction was stopped by heat inactivation at 65°C for 10 minutes. cDNA was transcribed from 4 μl of the ligated RNA by incubating the RNA in a mixture consisting of 0.5 μg random primer (Promega), dNTP mix (0.5 mM each final concentration) (Qiagen) and 20U RNasin at 65°C for 5 minutes, then cooling to 4°C, before adding 20U RNasin, 5 mM DTT, 200U Superscript II (Invitrogen) and 1 × First Strand buffer in a total reaction volume of 20 μl. The reaction mixture was subsequently incubated at 25°C for 5 minutes, followed by 50°C for 1 hour. 50 μl PCR reactions consisting of 4 μl cDNA (representing 20% of the reverse transcription reaction), dNTP mix (20 mM final concentration), universal primer (3′- or 5′- final concentration 0.2 μM) (Table 2), gene segment specific primer (3′- or 5′- final concentration 0.2 μM) (Table 2), 1 × Pfu buffer and 2.5U Native Pfu DNA polymerase (Stratagene) were made. The cycling conditions were as follows: initial denaturation at 96°C for 6 minutes, followed by 40 cycles of denaturation at 96°C for 30 seconds, primer annealing at 37°C for 1 minute and elongation at 72°C for 2 minutes. PCR reactions were analysed on a 2.5% agarose gel containing GelRed nucleic acid stain (Biotium) according to manufacturer’s directions. PCR products were purified using a QIAquick PCR purification kit (Qiagen) according to manufacturer’s directions. Where multiple bands were present in the gel, bands of the correct size were excised and purified using a QIAquick gel extraction kit (Qiagen) according to manufacturer’s directions.
Sequencing reactions were performed using the BigDye terminator sequencing kit version 3.1 (Applied Biosystems). M13 forward and reverse primers were used for gene segment PCR products, whilst for the non-coding regions the primers used for the PCR stage were reused, both at a final concentration of 80nM. The sequencing reactions were run on a 3130xl genetic analyzer (Applied Biosystems), and the resulting nucleotide sequences were visualised, assembled and edited using SeqMan II version 5.03 (DNAStar, Inc) and BioEdit version 188.8.131.52 (Ibis Pharmaceuticals Inc.).
This work was supported by the Horserace Betting Levy Board, EIP 2011–2012.
- McGeoch D, Fellner P, Newton C: Influenza virus genome consists of eight distinct RNA species. Proc Natl Acad Sci 1976, 73: 3045-3049. 10.1073/pnas.73.9.3045PubMedPubMed CentralView Article
- Waddell GH, Teigland MB, Sigel MM: A new influenza virus associated with equine respiratory disease. J Am Vet Med Assoc 1963, 143: 587-590.PubMed
- Daly JM, Lai AC, Binns MM, Chambers TM, Barrandeguy M, Mumford JA: Antigenic and genetic evolution of equine H3N8 influenza A viruses. J Gen Virol 1996, 77: 661-671. 10.1099/0022-1317-77-4-661PubMedView Article
- Bryant NA, Rash AS, Russell CA, Ross J, Cooke A, Bowman S, MacRae S, Lewis NS, Paillot R, Zanoni R, Meier H, Griffiths LA, Daly JM, Tiwari A, Chambers TM, Newton JR, Elton DM: Antigenic and genetic variations in European and North American equine influenza virus strains (H3N8) isolated from 2006 to 2007. Vet Microbiol 2009, 138: 41-52. 10.1016/j.vetmic.2009.03.004PubMedView Article
- Murcia PR, Wood JL, Holmes EC: Genome-scale evolution and phylodynamics of equine H3N8 influenza A virus. J Virol 2011, 85: 5312-5322. 10.1128/JVI.02619-10PubMedPubMed CentralView Article
- Gorman OT, Donis RO, Kawaoka Y, Webster RG: Evolution of influenza A virus PB2 genes: implications for evolution of the ribonucleoprotein complex and origin of human influenza A virus. J Virol 1990, 64: 4893-4902.PubMedPubMed Central
- Ito T, Gorman OT, Kawaoka Y, Bean WT, Webster RG: Evolutionary analysis of the influenza A virus M gene with comparison of the M1 and M2 proteins. J Virol 1991, 65: 5491-5498.PubMedPubMed Central
- Gorman OT, Bean WJ, Kawaoka Y, Webster RG: Evolution of the nucleoprotein gene of influenza A virus. J Virol 1990, 64: 1487-1497.PubMedPubMed Central
- Kawaoka Y, Krauss S, Webster RG: Avian-to-human transmission of the PB1 gene of influenza A viruses in the 1957 and 1968 pandemics. J Virol 1989, 63: 4603-4608.PubMedPubMed Central
- Okazaki K, Kawaoka Y, Webster RG: Evolutionary pathways of the PA genes of influenza A viruses. Virology 1989, 172: 601-608. 10.1016/0042-6822(89)90202-XPubMedView Article
- Nakajima K, Nobusawa E, Nakajima S: Evolution of the NS genes of the influenza A viruses. I. The genetic relatedness of the NS genes of animal influenza viruses. Virus Genes 1990, 4: 5-13. 10.1007/BF00308561PubMedView Article
- Worobey M, Han G-Z, Rambaut A: A synchronized global sweep of the internal genes of modern avian influenza virus. Nature 2014, 508: 254-257. 10.1038/nature13016PubMedPubMed CentralView Article
- Livesay GJ, O’Neill T, Hannant D, Yadav MP, Mumford JA: The outbreak of equine influenza (H3N8) in the United Kingdom in 1989: diagnostic use of an antigen capture ELISA. Vet Rec 1993, 133: 515-519. 10.1136/vr.133.21.515PubMedView Article
- Binns MM, Daly JM, Chirnside ED, Mumford JA, Wood JM, Richards CM, Daniels RS: Genetic and antigenic analysis of an equine influenza H3 isolate from the 1989 epidemic. Arch Virol 1993, 130: 33-43. 10.1007/BF01318994PubMedView Article
- Adeyefa CAO, Quayle K, McCauley JW: A rapid method for the analysis of influenza virus genes: application to the reassortment of equine influenza virus genes. Virus Res 1994, 32: 391-399. 10.1016/0168-1702(94)90087-6PubMedView Article
- Bean WJ: Correlation of influenza A virus nucleoprotein genes with host species. Virology 1984, 133: 438-442. 10.1016/0042-6822(84)90410-0PubMedView Article
- Dong BB, Xu CL, Dong LB, Cheng HJ, Yang L, Zou SM, Chen M, Bai T, Zhang Y, Gao RB, Li XD, Shi JH, Yuan H, Chen T, Zhu Y, Xiong Y, Yang S, Shu YL: A novel reassortant H3N8 influenza virus isolated from drinking water for duck in a domestic farm in Poyang Lake area. Biomed Environ Sci 2013, 26: 546-551.PubMed
- Skehel JJ, Hay AJ: Nucleotide sequences at the 5′ termini of influenza virus RNAs and their transcripts. Nucleic Acids Res 1978, 5: 1207-1219. 10.1093/nar/5.4.1207PubMedPubMed CentralView Article
- de Wit E, Spronken MI, Bestebroer TM, Rimmelzwaan GF, Osterhaus AD, Fouchier RA: Efficient generation and growth of influenza A/PR/8/34 from eight cDNA fragments. Virus Res 2004, 103: 155-161. 10.1016/j.virusres.2004.02.028PubMedView Article
- de Wit E, Bestebroer TM, Spronken MIJ, Rimmelzwaan GF, Osterhaus ADME, Fouchier RAM: Rapid sequencing of the non-coding regions of influenza A virus. J Virol Methods 2007, 139: 85-89. 10.1016/j.jviromet.2006.09.015PubMedView Article
- WHO: Sequencing primers and protocol. [http://www.who.int/csr/resources/publications/swineflu/sequencing_primers/en/] 
- GISAID EpiFlu database [http://platform.gisaid.org] 
- Lee K, Seong BK: The position 4 nucleotide at the 3′ end of influenza virus neuraminidase vRNA is involved in temporal regulation of transcription and replication of neuraminidase RNAs and affects the repertoire of influenza virus surface antigens. J Gen Virol 1998, 79: 1923-1934.PubMedView Article
- Jagger BW, Wise HM, Kash JC, Walters K-A, Wills NM, Xiao Y-L, Dunfee RL, Schwartzman LM, Ozinsky A, Bell GL, Dalton RM, Lo A, Efstathiou S, Atkins JF, Firth AE, Taubenberger JK, Digard P: An overlapping protein-coding region in influenza A virus segment 3 modulates the host response. Science 2012, 337: 199-204. 10.1126/science.1222213PubMedPubMed CentralView Article
- Shi M, Jagger BW, Wise HM, Digard P, Holmes EC, Taubenberger JK: Evolutionary conservation of the PA-X open reading frame in segment 3 of influenza A virus. J Virol 2012, 86: 12411-12413. 10.1128/JVI.01677-12PubMedPubMed CentralView Article
- Robertson JS: 5′ and 3′ terminal nucleotide sequences of the RNA genome segments of influenza virus. Nucleic Acids Res 1979, 6: 3745-3757. 10.1093/nar/6.12.3745PubMedPubMed CentralView Article
- Lee HK, Tang JW-T, Kong DH-L, Koay ES-C: Simplified large-scale Sanger genome sequencing for influenza A/H3N2 virus. PLoS One 2013, 8: e64785. 10.1371/journal.pone.0064785PubMedPubMed CentralView Article
- Hoffmann E, Stech J, Guan Y, Webster RG, Perez DR: Universal primer set for the full-length amplification of all influenza A viruses. Arch Virol 2001, 146: 2275-2289. 10.1007/s007050170002PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.