- Open Access
Metagenomic analysis of the turkey gut RNA virus community
Virology Journalvolume 7, Article number: 313 (2010)
Viral enteric disease is an ongoing economic burden to poultry producers worldwide, and despite considerable research, no single virus has emerged as a likely causative agent and target for prevention and control efforts. Historically, electron microscopy has been used to identify suspect viruses, with many small, round viruses eluding classification based solely on morphology. National and regional surveys using molecular diagnostics have revealed that suspect viruses continuously circulate in United States poultry, with many viruses appearing concomitantly and in healthy birds. High-throughput nucleic acid pyrosequencing is a powerful diagnostic technology capable of determining the full genomic repertoire present in a complex environmental sample. We utilized the Roche/454 Life Sciences GS-FLX platform to compile an RNA virus metagenome from turkey flocks experiencing enteric disease. This approach yielded numerous sequences homologous to viruses in the BLAST nr protein database, many of which have not been described in turkeys. Our analysis of this turkey gut RNA metagenome focuses in particular on the turkey-origin members of the Picornavirales, the Caliciviridae, and the turkey Picobirnaviruses.
Enteric disease syndromes such as Poult Enteritis Complex (PEC) in young turkeys and Runting-Stunting Syndrome (RSS) in chickens are a continual economic burden for poultry producers. The only reliable method to reproduce the clinical signs of these syndromes in experimental birds is oral inoculation with crude preparations of intestinal contents from naturally infected birds. Further, the full spectrum of the field signs observed associated with these syndromes is difficult to reproduce experimentally with isolated viruses [1, 2]. Numerous viruses are known to be circulating in turkey flocks in the United States, with recent research efforts targeting RNA viruses such as the turkey astroviruses, novel turkey-origin reovirus, and avian rotavirus, and DNA viruses such as the recently described turkey parvovirus [3–6]. However, there remains a possibility that an unidentified virus or combination of viruses may play a role in poultry enteric disease. Despite the isolation and characterization of many of these suspect viruses, the etiology of the poultry enteric disease syndromes remains elusive, and many enteric viruses can be detected in otherwise healthy turkey and chicken flocks [3, 4]. Regional and national enteric virus surveys have revealed the ongoing presence of avian reoviruses, rotaviruses and astroviruses in turkey and chicken flocks, with combinations of viruses often present in the poultry gut [3, 4]. A non-biased, comprehensive approach to virus discovery that would not require viral cultivation would reveal a great deal about the complex viral community in the turkey gut. Further, a community-based understanding of the viruses in the poultry gut will be an invaluable asset in ongoing studies of the enteric disease syndromes and would be welcome knowledge to poultry producers who rely upon efficient conversion of feed in the gut to produce an economically important commodity. A recent study utilizing a sequence-independent molecular screen of virus particle associated nucleic acid (PAN) in chicken enteric samples identified a novel chicken parvovirus (ChPV). This parvovirus is a member of the Parvovirinae sub-family within the Parvoviridae, and a PCR-based diagnostic test has been developed that targets the ChPV non-structural (NS) gene [6, 7].The success of this PAN procedure suggests that similar approaches can be used for virus discovery in the poultry gut . Ultra high-throughput nucleic acid pyrosequencing has emerged as a powerful diagnostic technology that can be applied to determine the full genomic repertoire present in a complex environmental sample . Viral metagenomics can be specifically utilized to analyze viral sequences in just about any sample type, and is a powerful tool for virus discovery [9–13]. Further, viral metagenomics can be specifically applied to the problem of determining etiology in diseases and disease syndromes with no known cause [14–16]. In order to characterize the un-described viruses present in the turkey gut, we utilized the Roche/454 Life Sciences GS-FLX pyrosequencing platform to compile an RNA virus metagenome from turkeys experiencing enteric disease. The present analysis focused on RNA viruses in the turkey gut due to the large number of RNA viruses that have been identified to date as possibly contributing to enteric disease and poultry production problems. This approach yielded numerous sequences homologous to viruses in the National Center for Biotechnology Information (NCBI) BLAST non-redundant (nr) protein database, many of which have not been described in turkeys. These results validate this metagenomic approach to identifying known and novel RNA viruses in the poultry gut. The sequence data generated via this approach will prove useful in the molecular characterization of the viral constituency of the poultry gut, and will inform the selection of molecular diagnostic tests for enteric viruses. This will facilitate the development of updated molecular diagnostic tests, and a more thorough knowledge of the viral constituency in the poultry gut will lead to a better understanding of the role viruses play in enteric disease and in the performance of poultry flocks in general.
Results and Discussion
The initial pyrosequencing runs produced in excess of 139,000,000 bases of high quality nucleotide sequence with an average read length of 362. The sequence data was used to assemble 6526 contigs ranging in size from 97 to 2578 bp, with the majority of contigs falling in the range of approximately 250 to 450 bp. 4563 contigs produced no hits in the nr protein database using the blastx search parameters and the MEGAN default settings. 724 contigs had similarity to sequences from cellular organisms, including bacteria, fungi and avian species. 788 contigs had similarity to RNA viral sequences, including sequences from the dsRNA viruses (Reoviridae and Picobirnaviruses), and the ssRNA viruses (Caliciviridae, Leviviridae, Picornavirales, and Astroviridae) (Figure 1). The number of cellular sequences in the present dataset are likely due to the use of intestinal homogenates, which included intestinal tissue in order to ensure the discovery of cell-associated viruses in the submitted samples. The tblastx search output produced a MEGAN taxon tree similar to the one presented in Figure 1 and revealed that many of the unassigned contigs were similar to avian sequences.
The majority of the assigned viral contigs (620) showed similarity to database sequences from the Picornavirales order and other picorna-like viruses, viruses that, as a group, contain a positive sense single-stranded RNA genome and a virion approximately 30 nm in diameter . Recently, a retrospective study of electron micrographs of enteric viruses from California turkeys experiencing enteric disease revealed a large number of "small round viruses" ranging in size from 15 to 30 nm, like most members of the Picornavirales . These small round viruses are present in turkeys across a range of ages, but they have only been identified morphologically, making specific identification difficult. Avian enterovirus-like viruses have been described for years in domestic poultry; again this designation has historically been made based primarily upon morphological characterization, and little is known about their pathogenicity or their transmission characteristics. It is unclear what role these picornaviruses and picorna-like viruses may play in turkey enteric disease or in turkey performance in general, but the presence of picornaviruses in other agricultural species has been closely associated with enteric disease, namely in pigs and cattle [20–23]. Members of the Picornavirales also infect avian species, with the etiologic agents of duck hepatitis (DH) in ducklings and avian encephalomyelitis (AE) in several poultry species both being picornaviruses [24, 25]. The present metagenomic analysis has identified RNA sequences with homology to seven of the nine recognized picornavirus genera , with the largest proportion of the sequences bearing homology to the Kobuvirus genus.
The picobirnaviruses (PBVs) are a relatively recently described group of viruses that contain dsRNA, bi-segmented genomes and have non-enveloped capsids generally around 35 nm in diameter . Since their initial description, the PBVs have been detected in enteric samples from several mammalian hosts, including humans [28–30]. A PBV has been described in chickens based upon morphological characterization and electropherotyping, along with a similar virus with an apparent tri-segmented genome . The chicken PBV was not specifically associated with enteric disease. PBVs have been associated with gastroenteritis in humans [32, 33], but a recent metagenomic analysis of the viral community in feces from healthy human volunteers revealed a relatively large number of PBV sequences . Interestingly, PBVs have been detected in humans together with rotaviruses and astroviruses [32, 34]. Phylogenetic analysis of a portion of the RNA-dependent RNA polymerase gene (RdRp) reveals that the putative turkey-origin PBV identified in the present metagenomic analysis is unique among the available PBV sequences, which include PBVs from humans, pigs, dogs, rats, snakes and municipal raw sewage (Figure 2). At the nucleotide level, the portion of the turkey-origin PBV used for the alignment and subsequent phylogenetic analysis shared from 49.5 to 70% sequence identity with the PBV sequences selected from the database; the highest identity was with the PBV detected in raw sewage in Washington state.
The family Caliciviridae includes four genera: Norovirus, Sapovirus, Lagovirus, and Vesivirus . The present metagenomic analysis revealed nucleic acid with homology to database sequences from the Sapovirus and Lagovirus genera. The Sapovirus genus includes viruses that cause enteritis in swine and humans, and the Lagoviruses infect rabbits and hares (lagomorphs) . In a phylogenetic analysis using a ~300 amino acid portion of ORF1 polyprotein homologous to the conserved P-loop NTPase superfamily, the putative turkey-origin Calicivirus grouped with the porcine enteric Sapoviruses (Figure 3). In general, the Sapoviruses are a very genetically heterogeneous genus , with numerous genogroups recognized, and many porcine Sapoviruses are related to human strains ; it will be interesting to determine the place the novel turkey caliciviruses hold in this genus, and it is notable that the portion of the ORF polyprotein analyzed in this study shared only 23.2 to 35% amino acid sequence identity with the Calicivirus isolates included in the alignment and subsequent phylogenetic analysis.
The appearance of avian reovirus and avian astrovirus sequences in the present analysis is not surprising, since these viruses are widespread in turkeys and chickens [3–5], and a good deal of recent research has helped to characterize the roles these viruses play in the poultry enteric syndromes [1, 2, 37–39]. The turkey astrovirus revealed using this metagenomic approach was most similar to the previously described type 2 turkey astroviruses (TAstV-2), which are very common in the turkey intestine . The avian reovirus genes revealed using this approach--namely the lambdaA core protein gene, the non-structural protein muC, and the structural protein muB--have not been described in turkeys previously, but they have homology to previously described avian reovirus genome segments [40, 41].
This first look at the turkey gut viral community using metagenomics has revealed a great deal of novel RNA virus sequence, and this analysis is a step toward identifying some of these undescribed, small enteric viruses. The dataset includes samples from areas that historically have had problems with enteric disease and includes a range of flock ages. The samples were collected regardless of the present enteric disease status of the flocks to ensure that the viral flora had not been perturbed by advanced enteric disease signs and due to the numerous observations that enteric viruses are often present in healthy birds [3, 4]. The sequence data generated via this approach will prove useful in the molecular characterization of the viral constituency of the poultry gut, and will inform the selection of molecular diagnostic tests for enteric viruses. Further, this study sets the stage for subsequent comparative metagenomic analyses to determine the viruses commonly found in flocks with enteric syndromes versus the viral constituency of healthy flocks, for regional comparisons of circulating enteric viruses, and for comparing specific management and nutritional techniques and their effect on the gut microbiome.
Materials and methods
Gut samples receipt and preparation
With the cooperation of industry stakeholders, complete intestinal tracts (from duodenum/pancreas to cloaca, including cecal tonsils) from five turkey farms in North Carolina, U.S.A. with histories of enteric disease problems were received at the Southeast Poultry Research Laboratory in October 2008. Five complete intestinal tracts from each farm were collected. The intestines were processed promptly via blending into a ~20% homogenate in sterile phosphate-buffered saline (PBS) and were pooled into a single sample. This pooled intestinal homogenate represented turkeys ranging in age from 7 days to 34 days. After 5000 rpm (2400 × G) and 7500 rpm (5500 × G) centrifugation steps (15 min, 4°C) to clarify the homogenate (SLA 1500 SuperLite rotor, Sorvall), a stepwise filtration process involving 0.8 μm, 0.45 μm, and 0.2 μm cutoff bottle filters (Nalgene) was used to remove large particles and bacteria. Virus-sized particles were pelleted by ultracentrifugation (5 hr., 113,000 × G, 4°C).
Viral RNA isolation and cDNA synthesis
The virus particle pellet was resuspended in Tris-HCl buffer (pH 7.5) and treated with RNAse A (Invitrogen) to remove unencapsidated (non-viral) RNA. Total RNA was extracted from the pellet using TRIZOL-LS reagent (Invitrogen) and RNA was further purified from the TRIZOL-LS aqueous phase using the MagMax Viral RNA isolation kit (Ambion). This RNA sample was treated with Turbo DNAse (Ambion) to minimize any remaining DNA. cDNA was generated with random hexamers using the Invitrogen SuperScript Choice System and was ligated to the included EcoRI/NotI double-stranded oligonucleotide adapters and sephacryl column purified per the manufacturer's instructions. PCR using an EcoRI/NotI adapter primer (5'-CGG CCG CGT CGA C-3') was used to amplify the cDNA. An initial PCR reaction of 50 μl was split into 5 equal reactions prior to thermal cycling (94°C for 30s, 55°C for 30s, 72°C for 2 min, times 35 cycles, followed by an extension at 72°C for 7 min). The reactions were pooled and purified (Qiagen MinElute PCR cleanup kit).
High throughput sequencing and analysis
The amplified cDNA was utilized in high-throughput nucleic acid sequencing using Genome Sequencer FLX Titanium pyrosequencing technology and reagents (Roche). Contigs were assembled using the gsAssembler software (454 Life Sciences) using stringent parameters (50 bp overlap with 95% identity). Using the assembled contigs as query sequences, the blast non-redundant (nr) protein database (GenBank) was searched using the blastx and tblastx programs. The blastx and tblastx output was compiled and contigs were assigned to taxa with MEGAN using the default LCA algorithm parameters . Nucleotide and amino acid sequences were aligned using ClustalW and phylogenetic trees were prepared using MEGA4; the relationship of the sequences was inferred using the Neighbor-Joining method, which was determined to be a computationally efficient method to deal with datasets of this size and produced single trees to illustrate an initial placement of these novel viruses among available sequences . In order to confirm directly the presence of picobirnavirus and calicivirus sequences in gut samples, the metagenomic contigs were utilized to design RT-PCR primers to amplify portions of the picobirnavirus RNA-dependent RNA polymerase (RdRp) and calicivirus ORF1 polyprotein. These primers were subsequently used to amplify viral sequences in the original RNA prep used to create the metagenome and in archived turkey intestinal samples from North Carolina, U.S.A (Figure 4).
Pantin-Jackwood MJ, Spackman E, Day JM: Pathogenesis of type 2 turkey astroviruses with variant capsid genes in 2-day-old specific pathogen free poults. Avian Pathol 2008, 37: 193-201. 10.1080/03079450801932200
Spackman E, Pantin-Jackwood M, Day J, Sellers H: The pathogenesis of turkey origin reoviruses in turkeys and chickens. Avian Pathology 2005, 34: 291-296. 10.1080/03079450500178501
Pantin-Jackwood M, Day JM, Jackwood MW, Spackman E: Enteric viruses detected by molecular methods in commercial chicken and turkey flocks in the United States between 2005 and 2006. Avian Dis 2008, 52: 235-244. 10.1637/8174-111507-Reg.1
Pantin-Jackwood MJ, Spackman E, Day JM, Rives D: Periodic monitoring of commercial turkeys for enteric viruses indicates continuous presence of astrovirus and rotavirus on the farms. Avian Dis 2007, 51: 674-680. 10.1637/0005-2086(2007)51[674:PMOCTF]2.0.CO;2
Pantin-Jackwood MJ, Spackman E, Woolcock PR: Molecular characterization and typing of chicken and turkey astroviruses circulating in the United States: implications for diagnostics. Avian Dis 2006, 50: 397-404. 10.1637/7512-020606R.1
Zsak L, Strother KO, Day JM: Development of a polymerase chain reaction procedure for detection of chicken and turkey parvoviruses. Avian Dis 2009, 53: 83-88. 10.1637/8464-090308-Reg.1
Zsak L, Strother KO, Kisary J: Partial genome sequence analysis of parvoviruses associated with enteric disease in poultry. Avian Pathol 2008, 37: 435-441. 10.1080/03079450802210648
Margulies M, Egholm M, Altman W, Attiya S, Bader J, Bemben L, Berka J, Braverman M, Chen YJ, Chen Z, et al.: Genome sequencing in microfabricated high-density picolitre reactors. Nature 2005, 437: 376-380.
Edwards RA, Rohwer F: Viral Metagenomics. Nat Rev Microbiol 2005, 3: 504-510. 10.1038/nrmicro1163
Breitbart M, Hewson I, Felts B, Mahaffy J, Nulton J, Salamon P, Rohwer F: Metagenomic Analyses of an Uncultured Viral Community from Human Feces. J Bacteriol 2003, 185: 6220-6223. 10.1128/JB.185.20.6220-6223.2003
Zhang T, Breitbart M, Lee W, Run JQ, Wei C, Soh S, Hibberd M, Liu E, Rohwer F, Ruan Y: RNA Viral Community in Human Feces: Prevalence of Plant Pathogenic Viruses. PLoS Biology 2006, 4: e3. 10.1371/journal.pbio.0040003
Valles SM, Hashimoto Y: Isolation and characterization of Solenopsis invicta virus 3, a new positive-strand RNA virus infecting the red imported fire ant, Solenopis invicta. Virology 2009, 388: 354-361. 10.1016/j.virol.2009.03.028
Culley A, Lang A, Suttle C: Metagenomic Analysis of Coastal RNA Virus Communities. Science 2006, 312: 1795-1798. 10.1126/science.1127404
Victoria JG, Kapoor A, Dupuis K, Schnurr DP, Delwart EL: Rapid identification of known and new RNA viruses from animal tissues. PLoS Pathog 2008, 4: e1000163. 10.1371/journal.ppat.1000163
Cox-Foster DLCS, Holmes EC, Palacios G, Evans JD, Moran NA, Quan PL, Briese T, Hornig M, Geiser DM, Martinson V, vanEngelsdorp D, Kalkstein AL, Drysdale A, Hui J, Zhai J, Cui L, Hutchison SK, Simons JF, Egholm M, Pettis JS, Lipkin WI: A metagenomic survey of microbes in honey bee colony collapse disorder. Science 2007, 318: 283-287. 10.1126/science.1146498
Ng TFF, Manire C, Borrowman K, Langer T, Ehrhart L, Breitbart M: Discovery of a novel single-stranded DNA virus from a sea turtle fibropapilloma by using viral metagenomics. J Virol 2009, 83: 2500-2509. 10.1128/JVI.01946-08
Le Gall O, Christian P, Fauquet CM, King AMQ, Knowles NJ, Nakashima N, Stanway G, Gorbalenya AE: Picornavirales, a proposed order of positive-sense single-stranded RNA viruses with a pseudo-T = 3 virion architecture. Archives of Virology 2008, 153: 715-727. 10.1007/s00705-008-0041-x
Woolcock PA, Shivaprasad HL: Electron Microscopic Identification of Viruses Associated with Poult Enteritis in Turkeys Grown in California 1993-2003. Avian Diseases 2008, 52: 209-213. 10.1637/8106-090607-Reg.1
Guy JS, McNulty MS, Hayhow CS: Avian Enterovirus-like viruses. In Diseases of Poultry. 12th edition. Edited by: Saif YM. Blackwell Publishing; 2008.
Reuter G, Egyed L: Bovine kobuvirus in Europe. Emerging and Infectious Diseases 2009, 15: 822-823. 10.3201/eid1505.081427
Khamrin P, Maneekarn N, Peerakome S, Okitsu S, Mizuguchi M, Ushijima H: Bovine kobuviruses from cattle with diarrhea. Emerging and Infectious Diseases 2008, 14: 985-986. 10.3201/eid1406.070784
Reuter G, Boldizsar A, Kiss I, Pankovics P: Candidate new species of Kobuvirus in porcine hosts. Emerging and Infectious Diseases 2008, 14: 1968-1970. 10.3201/eid1412.080797
Yu J, Jin M, Zhang Q, Li H, Li D, Xu Z, Li J, Cui S, Yang S, Liu N, Duan Z: Candidate porcine Kobuvirus, China. Emerging and Infectious Diseases 2009, 15: 823-825. 10.3201/eid1505.081518
Woolcock PR: Duck Hepatitis. In Diseases of Poultry. 12th edition. Edited by: Saif YM. Blackwell Publishing; 2008.
Calnek BW: Avian Encephalomyelitis. In Diseases of Poultry. 12th edition. Edited by: Saif YM. Blackwell Publishing; 2008.
Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA: Virus Taxonomy: Eighth report of the International Committee on Taxonomy of Viruses. Elsevier Academic Press; 2005.
Pereira HG, Flewett TH, Candeias JAN, Barth OM: A virus with a bisegmented double-stranded RNA genome in rat (Oryzomys nigripes) intestines. Journal of General Virologh 1988, 69: 2749-2754. 10.1099/0022-1317-69-11-2749
Fregolente MC, Gatti MS: Nomenclature proposal for picobirnavirus. Archives of Virology 2009, 154: 1953-1954. 10.1007/s00705-009-0537-z
Chandra R: Picobirnavirus, a novel group of undescribed viruses of mammals and birds: a minireview. Acta Virologica 1997, 41: 59-62.
Fregolente MC, de Castro-Dias E, Martins SS, Spilki FR, Allegretti SM, Gatti MS: Molecular characterization of picobirnaviruses from new hosts. Virus Research 2009, 143: 134-136. 10.1016/j.virusres.2009.03.006
Leite JPG, Monteiro SP, Fialho AM, Pereira HG: A novel avian virus with trisegmented double-stranded RNA and further observations on previously described similar viruses with bisegmented genome. Virus Research 1990, 16: 119-126. 10.1016/0168-1702(90)90016-5
Bhattacharya R, Sahoo GC, Nayak MK, Rajendran K, Dutta P, Mitra U, Bhattacharya MK, Naik TN, Bhattacharya SK, Krishnan T: Detection of Genogroup I and II human picobirnaviruses showing small genomic RNA profile causing acute watery diarrhoea among children in Kolkata, India. Infection, Genetics and Evolution 2007, 7: 229-238. 10.1016/j.meegid.2006.09.005
Grohmann GS, Glass RI, Pereira HG, Monroe SS, Hightower AW, Weber R, Bryan RT: Enteric Viruses and Diarrhea in HIV-Infected Patients. New England Journal of Medicine 1993, 329: 14-20. 10.1056/NEJM199307013290103
Bhattacharya R, Sahoo GC, Nayak MK, Ghosh S, Dutta P, Bhattacharya MK, Mitra U, Gangopadhyay D, Dutta S, Niyogi SK, et al.: Molecular epidemiology of human astrovirus infections in Kolkata, India. Infection, Genetics and Evolution 2006, 6: 425-435. 10.1016/j.meegid.2006.02.001
L'Homme Y, Brassard J, Ouardani M, Gagne MJ: Characterization of novel porcine sapoviruses. Archives of Virology 2010.
Reuter G, Zimsek-Mijovski J, Poljsak-Prijatelj M, Di Bartolo I, Ruggeri FM, Kantala T, Maunula L, Kiss I, Kecskemeti S, Halaihel N, et al.: Incidence, diversity, and molecular epidemiology of sapoviruses in swine across Europe. Journal of Clinical Microbiology 2010, 48: 363-368. 10.1128/JCM.01279-09
Day JM, Spackman E, Pantin-Jackwood M: Turkey origin reovirus induced immune dysfunction in specific pathogen free and commercial turkey poults. Avian Dis 2008, 52: 387-391. 10.1637/8190-120607-Reg
Pantin-Jackwood MJ, Spackman E, Day JM: Pathology and virus tissue distribution of Turkey origin reoviruses in experimentally infected Turkey poults. Vet Pathol 2007, 44: 185-195. 10.1354/vp.44-2-185
Spackman E, Day JM, Pantin-Jackwood MJ: Astrovirus, reovirus, and rotavirus concomitant infection causes decreased weight gain in broad-breasted white poults. Avian Dis 2010, 54: 16-21. 10.1637/8986-070909-Reg.1
Noad L, Shou J, Coombs KM, Duncan R: Sequences of avian reovirus M1, M2 and M3 genes and predicted structure/function of the encoded mu proteins. Virus Res 2006, 116: 45-57. 10.1016/j.virusres.2005.08.014
Xu W, Coombs KM: Avian reovirus L2 genome segment sequences and predicted structure/function of the encoded RNA-dependent RNA polymerase protein. Virol J 2008, 5: 153. 10.1186/1743-422X-5-153
Das A, Spackman E, Pantin-Jackwood MJ, Suarez DL: Removal of real-time reverse transcription polymerase chain reaction (RT-PCR) inhibitors associated with cloacal swab samples and tissues for improved diagnosis of Avian influenza virus by RT-PCR. J Vet Diagn Invest 2009, 21: 771-778.
Huson DH, Auch AF, Qi J, Schuster SC: MEGAN analysis of metagenomic data. Genome Research 2007, 17: 377-386. 10.1101/gr.5969107
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analyisis (MEGA) software version 4.0. Molecular Biology and Evolution 2007, 24: 1596-1599. 10.1093/molbev/msm092
The authors thank Laura Ferguson and the USDA ARS/SAA sequencing facility for excellent technical support, and industry stakeholders for their continued interest and support. This research was supported by USDA/ARS CRIS project 6612-32000-054-00 and by a grant from the United States Poultry and Egg Association.
The authors declare that they have no competing interests.
JMD conceived and coordinated the study, performed the virus isolation and cDNA production, and wrote the paper. LLB performed the bioinformatic analyses. MVD prepared the samples and performed the 454 pyrosequencing. BES coordinated the bioinformatic analyses and pyrosequencing. LZ provided technical input and intellectually contributed to the study design. All authors read and approved the final manuscript.