key: cord-0001242-4cniomio authors: Phan, Tung Gia; Vo, Nguyen Phung; Boros, Ákos; Pankovics, Péter; Reuter, Gábor; Li, Olive T. W.; Wang, Chunling; Deng, Xutao; Poon, Leo L. M.; Delwart, Eric title: The Viruses of Wild Pigeon Droppings date: 2013-09-04 journal: PLoS One DOI: 10.1371/journal.pone.0072787 sha: 6ad08282a281d311d1a9e87b8085a6f0c968ad50 doc_id: 1242 cord_uid: 4cniomio Birds are frequent sources of emerging human infectious diseases. Viral particles were enriched from the feces of 51 wild urban pigeons (Columba livia) from Hong Kong and Hungary, their nucleic acids randomly amplified and then sequenced. We identified sequences from known and novel species from the viral families Circoviridae, Parvoviridae, Picornaviridae, Reoviridae, Adenovirus, Astroviridae, and Caliciviridae (listed in decreasing number of reads), as well as plant and insect viruses likely originating from consumed food. The near full genome of a new species of a proposed parvovirus genus provisionally called Aviparvovirus contained an unusually long middle ORF showing weak similarity to an ORF of unknown function from a fowl adenovirus. Picornaviruses found in both Asia and Europe that are distantly related to the turkey megrivirus and contained a highly divergent 2A1 region were named mesiviruses. All eleven segments of a novel rotavirus subgroup related to a chicken rotavirus in group G were sequenced and phylogenetically analyzed. This study provides an initial assessment of the enteric virome in the droppings of pigeons, a feral urban species with frequent human contact. Many infectious diseases in humans are caused by pathogens originating from a wide variety of animals. More than 60% of emerging diseases are estimated to originate from wildlife [1] [2] [3] [4] . Public awareness of zoonoses has recently increased because of their public health and economic impacts. Birds are recognized as frequent reservoirs for viruses that are of concern to humans; notably influenza A which is capable of infecting other mammals thereby facilitating genome segment reassortments and changes in tropism and transmission efficiency [5] [6] [7] [8] . Sporadic human infections of the virulent H5N1 resulting from direct contact with infected poultry or wild birds have been reported in 15 countries, mainly in Asia [9] [10] [11] [12] [13] , and H7N9 has recently emerged as a virus of concern. The prevalence of avian influenza viruses was 12% of oropharyngeal and 20% of cloacal swab specimens collected from urban pigeons in Slovakia [14] . H5N1 was found in a dead feral pigeon in Hong Kong [15] but is generally apathogenic in this host species and the overall risk of H5N1 transmission from pigeons to humans or chickens appears low [16, 17] . West Nile virus (WNV) and Saint Louis encephalitis (SLE) virus, two arboviruses in the Flavivirus genus transmitted by mosquitoes bites, are disseminated by wild birds [18] [19] [20] [21] [22] [23] . WNV-specific antibody and viremia was found in 25.7% and 11% of rock pigeons, respectively in the United States [24] . WNV was also isolated in pools of brains, kidneys, heart and spleen of feral pigeons and mapgies [25] . Pigeons developed low levels of WNV viremia; insufficient to infect mosquitoes [24, 26] . Avian paramyxoviruses, including Newcastle disease virus, are common domestic and wild bird pathogens [27] [28] [29] [30] . Paramyxovirus type-1 can be found in pigeons worldwide [31] [32] [33] [34] but the clinical signs vary depending on the immunity of the host and virulence of the specific isolates [35] . While human infection with Newcastle disease virus is rare, at least two outbreaks of conjunctivitis due to Newcastle disease virus have been reported in poultry workers [28, 36, 37] . Chicken anemia virus (CAV), until recently the only member of the gyrovirus genus, is highly contagious and causes severe anemia, hemorrhage and depletion of lymphoid tissue in chickens [38] [39] [40] . Related gyroviruses were recently characterized in human feces, blood and on healthy human skin [41] [42] [43] [44] indicating possible human tropism. Gyrovirus DNA was also detected in three blood samples of solid organ transplant patients and in one HIV-infected person [44] as well as in 0.85% of healthy French blood donations [45] . Pigeons are therefore natural reservoirs for pathogens that have caused emerging and re-emerging diseases in humans. In order to better understand the viruses shed by pigeons to which humans are frequently exposed, we genetically characterized the viral community in droppings from wild pigeons in Hong Kong and Hungary following an unbiased amplification method and deep sequencing. Biological Samples and viral Metagenomics 51 fecal specimens were collected from feral pigeons (Columba livia) in Hong Kong (N = 50) in August 2011 and in Pécs, Hungary (N = 1) in April 2011, then stored at 280uC. Droppings were collected without any physical contact with wild urban pigeons. Procedure therefore did not require animal committee ethical review. The field studies did not involve endangered or protected species. Fresh and well-separated droppings were sampled using sterile swabs which were stored in two ml of viral transport media as described [46] . The fecal sample from Pécs, Hungary was collected from a nesting 5-8 days old, clinically healthy pigeon into Hanks' buffered saline solution (Gibco BRL). Prior to the shipment on dry ice, samples were heat treated at 74uC for 50 minutes to inactivate viruses. The suspensions were vigorously vortexed and clarified by 15,000 6g centrifugation for 10 minutes. 200 ml of the supernatant was filtered through a 0.45-mm filter (Millipore) to remove bacterium-sized particles. The filtrate was then treated with a mixture of DNases (Turbo DNase from Ambion, Baseline-ZERO from Epicentre, and Benzonase from Novagen) and RNase (Fermentas) to digest unprotected nucleic acids [47] . Viral nucleic acids protected from digestion within viral capsids and other small particles were then extracted using a QIAamp spin-column technique according to the manufacturer's instructions (Qiagen). For the single fecal specimen from Hungary, non-specific RNA and DNA amplification was then performed by random RT-PCR using a primer with randomized 39 ends [47] . Amplicons were then pooled, 454 libraries generated according to the manufacturer's protocol (GS FLX Titanium General Library Preparation Kit, Roche) and pyrosequenced using the 454 Titanium FLX+ sequencer. For the 50 fecal specimens from Hong Kong, pools of nucleic acids from five fecal specimens were generated resulting in ten separate pools. A library was then constructed using ScriptSeq TM v2 RNA-Seq Library Preparation Kit (Epicentre) and then sequenced using the Miseq Illumina platform with 250 bases paired ends with a distinct molecular tag for each pool. Because fecal specimens were analyzed in pools of five, the viruses identified are reported per specimen pool rather than for individual specimens. The 454 and Illumina singleton reads and assembled contigs greater than 100-bp were compared to the GenBank protein databases using BLASTx. An E value of 10 25 was used as the cutoff value for significant hits. The sequence data from the MiSeq run is in the short read archive at GenBank accession number SRX263026. Sequences showing significant but divergent BLASTx hits to parvoviruses, picornaviruses or rotaviruses in the same sample pool were linked together using RT-PCR or PCR and primers based on the initial short sequence reads or contigs. 59 and 39 rapid amplification of cDNA end (59 and 39 RACE) was used to acquire the 59 and 39 extremities of viral genome. Phylogenetic analyses were performed using novel virus sequences, their best BLASTx hits, and representative members of related viral species or genera. All alignments and phylogenetic analyses were based on the translated amino acid sequences. Sequence alignment was performed using CLUSTAL X (version 2.0.3) with the default settings [48] . Aligned sequences were trimmed to match the genomic regions of the viral sequences obtained in the study. A phylogenetic tree with 100 bootstrap resamples of the alignment data sets was generated using the neighbor-joining method based on the Jones-Taylor-Thornton matrix-based model in MEGA version 5 [49] . Bootstrap values (based on 100 replicates) for each node are given if .70%. The generated phylogenetic trees were visualized using the program MEGA version 5. Resulting trees were examined for consistency with published phylogenetic trees. Sequence identity was measured using BioEdit [50] . The sequence distance comparison was calculated by using SSE [51] . Conserved amino acid analyses were performed over the alignments produced by CLUSTAL X (version 2.0.3). The hypothetical cleavage map of the picornavirus polyprotein was derived from alignments with other closest picornaviruses and NetPicoRNA prediction. Putative ORFs in the genome were predicted by NCBI ORF finder. The secondary structure of 59 and 39 UTR was predicted using the Mfold program. We received 36.5 million 250-bp paired-end reads generated on the Illumina Mi-Seq platform. Short reads were debarcoded using vendor software from Illumina. A total of 6,818 reads from 454 pyrosequencing were trimmed of their primer sequence and adjacent eight nucleotides corresponding to the randomized part of the primer. In-house analysis pipelines were developed to process both datasets. Adaptors were trimmed by in-house code. The cleaned reads were de-novo assembled using SOAPdenovo2 (38) . The assembled contigs, along with singlets, were aligned to an in-house viral protein database using NCBI BLASTx. The significant hits to viral sequences were then aligned to an inhouse non-virus-non-redundant (NVNR) universal protein database using BLASTx. Hits with more significant adjusted E-value to NVNR than to virus were removed. The distributions of the viral hits among the ten pools of five fecal samples from Hong Kong are shown in Table 1 based on BLASTx E scores ,10 25 . A range of best hits (BLASTx E scores ,10 25 ) to different viral families infecting plants (Germinividae, Nanoviridae, Caulimoviridae, Hypoviridae, Umbravirus, Mycovirus, Totiviridae, Tombusviridae) and insects (Iflavirus, Tetnovirus, Chronic bee paralysis virus, Laem Singh virus) were detected ( Table 1 ). The mammalian viral hits belonged to the Circoviridae, Parvovirinae subfamily, Picornavirus, Rotavirus, Adenovirus, Astrovirus, and Calicivirus. The Circoviridae-related sequences showed best hits to three different clades, Circovirus (1824 reads), Cyclovirus (20 reads) and Gyrovirus (2 reads). 77% of the circovirus hits showed best matches to pigeon circovirus of which 108 sequences had .90% aaidentity to pigeon circovirus. 74% of the Parvovirinae sequences showed best hits to AAV (Adeno-associated virus). A majority of the picornavirus sequences (85%; 759/889) showed best hits to Turkey hepatitis virus [52] . The adenovirus sequences showed best hits to three difference genera Mastadenovirus (3 reads), Aviadenovirus (77 reads) and Siadenovirus (296 reads). The pigeon fecal sample from Hungary was analyzed using 6,818 unique reads from 454 pyrosequencing and showed viral hits to picornaviruses, astroviruses, and tobamoviruses. We further characterized a subset of divergent mammalian viruses by acquiring their nearly complete or full genomes and phylogenetically comparing them to their closest viral relatives. Table 1 . Distribution of sequence reads to different viral species/families in 10 sample pools from Hong Kong and the single sample from Hungary. [57, 58] . Parvoviruses are widespread pathogens that can cause a range of diseases in a variety of mammals as well as birds and reptiles [57] . One hundred and five out of a total of 39,193 reads from one specimen pool (Table 1 pair 2) were related to parvovirus proteins (BLASTx E scores: 7610 26 to 5610 241 ). Gaps between short sequence reads were filled by PCR and then directly Sangersequenced. The nearly complete genome (5492-bp) of pigeon parvovirus, including a partial 59 UTR (252-bp), complete NS (631-aa), complete VP (696-aa) and a partial 39 UTR (54-bp) was acquired (GenBank KC876004). As expected the genome contained two large open reading frames (ORFs) with the left and right ORFs encoding non-structural (NS) and viral capsid protein (VP), respectively ( Figure 1A ). Distinct from its closest genetic relatives (chicken and turkey parvoviruses), pigeon parvovirus A contained an unusually long middle ORF (477-aa) ( Figure 1A) . The carboxy termini of this ORF (135-aa; position 343-477) showed best although weak similarity (BLASTx to NR E score: 8.9610 22 ) to a segment of a 221-aa ORF43 predicted encoded protein from the large dsDNA fowl adenovirus genome [59] , sharing 28% amino acid identity ( Figure 1A) . A smaller middle ORF (155-aa) was also detected, which did not show Table S4 . doi:10.1371/journal.pone.0072787.g002 similarity to any other sequence in GenBank. The start codon of NS was located in a strong Kozak sequence context, CAG-CATGGC and contained the ATP or GTP binding Walker loop motif 401 GPANTGKT 408 [GXXXXGK(T/S)] [60] . Additionally, two conserved replication initiator motifs 118 RCHVHIMLI 126 and 145 ITKYVTEALT 154 were also found [61] . Similar to chicken parvovirus, pigeon parvovirus's VP protein did not possess the phospholipase A 2 (PLA 2 ) motif with its highly conserved calcium-binding site (YLGPF) as well as the phospholipase catalytic residues (HD and D). The N-terminus of its VP protein contained glycine rich sequence (GGGGSVGSGGGGGVG) also present in other parvoviral VP proteins. Pair-wise amino acid sequence analysis showed that the NS and VP proteins in pigeon parvovirus shared the highest aa-identities of 41% and 34% to those of chicken parvovirus [62] , and less than 20% to the proteins of other parvovirus genera (Table S1 ). Figure 1B shows that pigeon parvovirus NS and VP proteins shared a monophyletic root with chicken and turkey parvoviruses, indicating a common origin. Based on the NS and VP phylogenetic analysis and genetic distance calculations, parvoviruses from pigeons, chickens and turkeys may be included as members of as a novel genus with a proposed name of Aviparvovirus (for avian parvovirus) in the subfamily Parvovirinae. Picornaviruses are small, non-enveloped, single-stranded RNA viruses whose prototype is poliovirus. The family Picornaviridae currently consists of seventeen genera (http://www. picornaviridae.com) that infect a wide range of hosts including mammalian and bird species. Many other recently characterized picornavirus genomes have increased the number of potential Picornaviridae genera [63] [64] [65] [66] . The genome of a typical picornavirus ranges from 7 Kb to 9 Kb and, with one exception, contains a single long ORF coding for a polyprotein [67] . Picornaviruses have been reported in fecal and respiratory specimens from numerous vertebrates including humans, bats, rodents, pigs, fish and reptiles [65, [68] [69] [70] . In this study we characterize two complete genomes of a novel picornavirus from pigeon fecal samples collected from Hong Kong (strain HK-21, KC876003) and Hungary (strain GALII-5/2011/ HUN, KC811837) with a 9,072/9,192 nucleotide (nt) long genome [excluding the poly(A) tail] (Figure 2A ). The P1, P2 and P3 regions of GALII-5/2011/HUN showed 74%, 92% and 97% amino acid identity to strain HK-21 indicating a presence of two genotypes of the same picornavirus species. Polyprotein alignments revealed that strain GALII-5/2011/HUN had an alternative upstream start codon MREY instead of the putative start codon of MATF seen in strain HK-21 ( Figure S2) . A stop codon was found between MREY and MATF in strain HK-21 but none was found in strain GALII-5/2011/HUN. The start codon in the upstream MREY in GALII-5/2011/HUN was also located in a weaker Kozak consensus sequence [GGGGATGCG] which plays a major role in the initiation of the translation process than MATF [GGAGATGGC] ( Figure S2 ). We therefore selected the second methionine codon (MATF) as the start codon for the analysis. The single 2,707/2,711-aa-long (of HK-21/GALII-5/2011/HUN) polyprotein coding region was flanked by the 620/727 nt-long 59 UTR and the 320/332-nt-long 39 UTR. By BLASTx, its polyprotein showed the closest identity to the turkey hepatitis virus belonging to Megrivirus genus [52] . Therefore, we provisionally named these piconaviruses as Mesivirus-1 (strain HK-21) and Mesivirus-2 (strain GALII-5/2011/HUN) for Megrivirus sisterclade virus. The mesivirus 59 UTR contained a highly conserved, 20 nucleotides long motif and apical ''8''-like structure 407 GGAGGTGGTGCTGAAATATTGCAAGCCACT 437 (with unknown function) that were also seen in avian origin picornaviruses, including turkey hepatitis virus (genus Megrivirus), duck hepatitis A virus-1 (genus Avihepatovirus), quail picornavirus (unassigned genus) and pigeon picornavirus B (unassigned genus) [52, 71, 72] . All these picornaviruses had type IV-like IRES [http://www.picornaviridae.com], however, we could not draw the potential secondary structure of mesivirus IRES. The hypothetical cleavage map of the polyproteins of mesiviruses were derived from alignments with other known picornaviruses and NetPicoRNA prediction [73] . Mesiviruses did not have a putative L protein preceding the capsid region. The P1 regions were 807/ 809-aa (Mesivirus-1/Mesivirus-2) in length containing rhv-like superfamily domains (Figure 2A ) revealed by Conserved Domain Database (CDD) search [74] . The conserved motif GXXXT/S for myristylation was also absent. Similar to the megrivirus, there was no potential VP2/VP4 cleavage site in mesiviruses. The P1 polyproteins of mesivirus-1 and -2 could be theoretically cleaved at VP0/VP3 (Q 389-390 QY/T), VP3/VP1 (Q 557/558 QG) and VP1/ 2A1 (Q 807/809 QD/E). BLASTx showed that the P1 region had the highest aa-identity of 43%-41% to the megrivirus (strain 0091.1, HQ189775) ( [77] . Phylogenetic analysis of P1-P3 also showed that mesiviruses were most closely related to megrivirus ( Figure 2B and Figure S1B ) and in pair-wise comparison mesivirus-1/mesivirus-2 had identities of 43%/41%, 34%/30% and 47%/47% to turkey megrivirus (strain 0091.1) at the amino acid level for the P1, P2 and P3 regions respectively (Table S2 ). According to the ICTV, the member of a picornavirus genus should share .40%, .40% and Table S5 . doi:10.1371/journal.pone.0072787.g003 .50% aa-identities in P1, P2 and P3 regions respectively (http:// www.picornaviridae.com). Although the identities of the P2 and P3 regions qualified mesiviruses as belonging to a novel genus, its P1 region shared more than the 40% aa-identity threshold to megrivirus. Mesiviruses might therefore be considered as a new picornavirus species in the genus Megrivirus or given their distinct 2A1 and high P2 and P3 divergence to megrivirus as members of a separate genus. Rotaviruses consist of at least eight groups or species (A through H) with multiple P and G genotypes which together comprise a genus in the family Reoviridae. Rotavirus has a non-enveloped, triple-layered icosahedral capsid containing eleven segments of double-stranded RNA, encoding for six structural and five nonstructural proteins [78] [79] [80] . We successfully acquired 17,827 nucleotides of this highly divergent rotavirus (HK18), encoding the complete proteins of all 11 genome segments (GenBank KC876005-KC876015). The six structural proteins (VP1-VP7) showed best hits by BLASTx to group G rotavirus, first described in chicken feces from Northern Ireland [81] . While VP1-VP3 shared 84%-94% aa-identities with group G rotavirus, other structural proteins showed lower aa-identities, ranging from 35% to 67% ( Table 2) . Classification of rotavirus is based mainly on the inner capsid protein VP6 and according to ICTV definition the members of the same rotavirus species should share .60% nt-identity in VP6 region [82] . The sequence (1176-nt) of the VP6-encoding genome segment shared 68% nt-identity to group G rotavirus ( Figure 3A ). Direct observation of VP6 alignment between HK18 and its closest relatives revealed numerous mismatches in the groupspecific antigenic regions ( Figure S3 ), indicating that HK18 represents a candidate prototype for a novel subgroup within group G. Two structural proteins VP4 (the protease-cleaved polyprotein) and VP7 (the glycoprotein) on the outermost surface of rotavirus are targets for neutralizing antibodies and have also been used to classify rotaviruses into their P and G genotypes, respectively [83] . The nucleotide sequence alignments of HK18 VP4 and VP7 showed low identity to other rotaviruses ( Figure 4A ). In addition, HK18 and group G rotavirus shared only 35% for VP4 and 51% for VP7 at the amino acid level, suggesting that HK18 should be considered as novel P and G genotypes, tentatively proposed as G2P [2] within group G rotavirus. The VP1 protein of HK18 contained three conserved RNA-dependent RNA polymerase motifs 558 AEKIILYTDVSQWDAS 573 , 638 LKI RYLGVASGEKTTKIGNSFANVALI 664 and 683 MRVDGDD NVVT 693 [84] . Similar to group G rotavirus, the VP3 protein of HK18 did not possess the NTP-binding sequence motif present in group C rotavirus [85] . However, an alignment of VP3 sequences from different rotavirus groups showed that the VP3 protein of HK18 shared a high degree of conservation in ALYXLSN [ALYSLSN] motif, which has unknown function [86] . In order to phylogenetically classify pigeon rotavirus HK18, its six structural proteins were aligned to other group rotavirus representatives. While the VP1-VP3, VP6 and VP7 of HK18 were clustered with those of group G rotavirus ( Figure 3B and Figure S4 ), its VP4 was located on a branch diverging closer to the root ( Figure 4B ). The pair-wise amino acid sequence analysis demonstrated that nonstructural proteins NSP2-NSP5 shared top identities to group G rotavirus rather than other groups, ranging from 50% to 92% at the amino acid level (Table 2 ). However, NSP1 had a very low aaidentity, less than 31% to group G rotavirus. Like group G rotavirus, the NSP1 gene of HK18 had two ORFs encoding two minor and major peptides labeled as NSP1-1 (104-aa) and NSP1-2 (310-aa), respectively. While NSP1-1 had the best identity of 30% to group G rotavirus, NSP1-2 shared the closet match to group B rotavirus with 23% identity versus 19% to group G. The NSP1 zinc-binding domain could not be identified in HK18. The NSP2 was found to have a conserved sequence HGXGHXRXV and histidine triad (His-X-His-X-His-XX) located at the RNA binding domain [87] . Phylogenetic analysis of NSP2-NSP5 showed that the group G rotavirus was the closest relative to HK18 ( Figure S4 ) in all fragments except for NSP1 where it appeared basal to both groups G and B. A similar observation was also made for avian group A rotaviruses where their NSP1 could not be grouped into taxonomic species [88] . Metagenomics has been used to analyze viral nucleic acids in feces collected from humans and a growing range of animals including primates, horse, bats, rodents, pigs, dogs, and turkeys [42, 58, 65, 69, [89] [90] [91] [92] . We investigated viral sequences in 51 wild pigeon droppings followed by viral sequence similarity searches revealing sequences closely or distantly related to viral genomes in GenBank. The eukaryotic viruses detected included insect and plant viruses likely reflecting the pigeon diet as well as viruses already known to infect pigeons such as pigeon circovirus and adenoviruses. Also detected and characterized were novel genomes in viral families known to infect vertebrates. Circovirus-like inclusion bodies were initially identified in the bursa of pigeons [93, 94] and the molecular epidemiology of pigeon circovirus worldwide has confirmed its importance as a pathogen associated with a wide range of illnesses in pigeon populations, including weight loss, respiratory distress and diarrhea [94] [95] [96] ; possibly as a result of induced immunodeficiency aggravating the pathogenicity of co-infections such as those of adenoviruses [97] . Both pigeon circovirus and diverse adenovirus sequences were identified here. Parvoviruses cause a variety of mild to severe symptoms in birds [98] [99] [100] [101] . The Parvoviridae subfamily has recently undergone a large expansion in the number of know genera and species [58, 62, 102, 103] . Here a novel parvovirus was characterized in pigeon feces that was distantly related to chicken and turkey parvoviruses found in feces of farm animals with signs of enteric diseases [62, 98] . The pigeon parvovirus genome, the first from that host species, possessed an unusually long middle ORF that encoded a protein showing similarity to that encoded by an ORF of unknown function in the fowl adenovirus genome. This observation might reflect past lateral gene exchange between avian DNA viruses replicating in the nucleus. The pigeon, chicken and turkey parvoviruses cluster together phylogenetically, but have not yet been assigned to a Parvovirinae genus. We therefore propose a new genus provisionally called Aviparvovirus containing at least two species, pigeon parvovirus and the closely related turkey/ chicken parvoviruses. Two related picornaviruses (81% nucleotide similarity) from pigeons in Hong Kong and Hungary were also characterized and provisionally named Mesivirus-1 and -2 reflecting a wide geographic distribution of this viral clade. Based on the phylogenetic analysis of the most conserved region (P3) many of the bird picornaviruses clustered together (Gallivirus, Orthoturdovirus, Paraturdivirus, Megrivirus, and Mesivirus) in a supported clade that includes only the mammalian infecting Kobuvirus and Salivirus genera. It is therefore possible to speculate on the past existence of a strictly avian clade a member of which adapted to mammals later resulting in the Salivirus/Kobuvirus clade. Rotaviruses were first described as a causative agent of gastroenteritis in humans in 1973 by using electron microscopy to examine biopsies of duodenal mucosa from children with acute non-bacterial gastroenteritis [104] . The first pigeon rotavirus was isolated from feces in 1983 and belongs to group A [105] . We describe here the complete coding regions of all eleven segments of a group G rotavirus, the first from pigeons, distinct from chicken rotavirus [85] . The addition of a second pigeon rotavirus genome to the viral database also shows that the diversity of known avian rotaviruses is likely to continue to expand with multiple rotavirus species infecting the same host species. Animal species living in or around human habitation, including some bats and rodents species and pigeons, are recognized as the reservoirs of multiple zoonotic pathogens. Because of increasing contact between such peridomestic animals and humans a better understanding of these animals' virome can inform future studies of their cross-species transmissions potential. Monitoring viral exchanges between pigeons and highly exposed humans using nucleic acids or serological assays for avian viruses will be facilitated by improved knowledge of their virome. Table S1 Pairwise amino acid sequence identities (%) between NS and VP regions of the novel pigeon parvovirus, turkey parvovirus and representatives of parvovirus genera. GenBank numbers of these viruses are available in Table S3 . (PDF) Table S2 Coding potential/putative proteins of the genome of Mesivirus-1 and comparison of amino acid sequence identity (%) of its P1-3 and other closelyrelated picornaviruses in the Picornaviridae family. (PDF) Table S3 Representative members in the subfamily Parvovirinae for the phylogenetic trees in Figure 1 and their GenBank numbers. *Parvoviruses were used for pairwise calculation in Table S1 . (PDF) Table S4 Representative members in the family Picornaviridae for the phylogenetic tree in Figure 2 and their GenBank numbers. (PDF) Reservoirs and vectors of emerging viruses The role of wildlife in emerging and re-emerging zoonoses Public health threat of new, reemerging, and neglected zoonoses in the industrialized world Wildlife, exotic pets, and emerging zoonoses Recent zoonoses caused by influenza A viruses Is avian influenza virus A(H5N1) a real threat to human health? Influenza a virus entry: implications in virulence and future therapeutics Perspectives on the global threat: the challenge of avian influenza viruses for the world's veterinary community A global perspective on avian influenza Control strategies for highly pathogenic avian influenza: a global perspective The next influenza pandemic: lessons from Hong Kong Epidemiology of St. Louis encephalitis virus History and evolution of HPAI viruses in southeast Asia Determination of hemagglutinin and neuraminidase subtypes of avian influenza A viruses in urban pigeons by a new nested RT-PCR Investigation of outbreaks of highly pathogenic H5N1 avian influenza in waterfowl and wild birds in Hong Kong in late Limited susceptibility of pigeons experimentally inoculated with H5N1 highly pathogenic avian influenza viruses Influenza A virus monitoring in urban and free-ranging pigeon populations in Germany West Nile virus infection and immunity The virology, epidemiology, and clinical impact of West Nile virus: a decade of advancements in research since its introduction into the Western Hemisphere A global perspective on the epidemiology of West Nile virus West nile virus disease and other arboviral diseases -United States West Nile virus: North American experience Unifying the spatial epidemiology and molecular evolution of emerging epidemics West Nile virus viremia in wild rock pigeons Reemergence of West Nile virus in Italy Experimental infection of North American birds with the New York 1999 strain of West Nile virus A meta-analysis of observational epidemiological studies of Newcastle disease in African agro-systems The long view: a selective review of 40 years of Newcastle disease research Newcastle disease in the European Union Newcastle disease: a review of field recognition and current methods of laboratory detection Characteristics of pigeon paramyxovirus serotype-1 isolates (PPMV-1) from the Russian Federation from Phylogenetic analysis and comparison of eight strains of pigeon paramyxovirus type 1 (PPMV-1) isolated in China between Biological and molecular characterization of a pigeon paramyxovirus type-1 isolate found in Argentina Pigeon paramyxovirus-1 (P-group) as the cause of severe outbreaks in fancy Columba livia in Saudi Arabia A comparative infection study of pigeon and avian paramyxovirus type 1 viruses in pigeons: evaluation of clinical signs, virus shedding and seroconversion An outbreak of conjunctivitis due to Newcastle disease virus (NDV) occurring in poultry workers Outbreaks of conjunctivitis due to the Newcastle disease virus among workers in chicken-broiler factories Chicken anemia virus Immunopathogenesis of chicken anemia virus infection Economically important non-oncogenic immunosuppressive viral diseases of chicken-current status Characterization of a novel gyrovirus in human stool and chicken meat A third gyrovirus species in human faeces Discovery of a genome of a distant relative of chicken anemia virus reveals a new member of the genus Gyrovirus Human gyrovirus DNA in human blood Human gyrovirus in healthy blood donors, france Avian coronavirus in wild aquatic birds Metagenomic analyses of viruses in stool samples from children with acute flaccid paralysis The neighbor-joining method: a new method for reconstructing phylogenetic trees MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT SSE: a nucleotide and amino acid sequence analysis platform Novel picornavirus in Turkey poults with hepatitis The expanding range of parvoviruses which infect humans The emergence of parvoviruses of carnivores Animal bocaviruses: a brief review Densoviruses for control and genetic manipulation of mosquitoes A newly identified bocavirus species in human stool Acute diarrhea in West African children: diverse enteric viruses and a novel parvovirus genus Coding potential and transcript analysis of fowl adenovirus 4: insight into upstream ORFs as common sequence features in adenoviral transcripts Distantly related sequences in the alpha-and beta-subunits of ATP synthase, myosin, kinases and other ATP-requiring enzymes and a common nucleotide binding fold Conserved sequence motifs in the initiator proteins for rolling circle DNA replication encoded by diverse replicons from eubacteria, eucaryotes and archaebacteria Determination and analysis of the full-length chicken parvovirus genome Genetic characterization of a novel picornavirus in turkeys (Meleagris gallopavo) distinct from turkey galliviruses and megriviruses and distantly related to the members of the genus Avihepatovirus Two closely related novel picornaviruses in cattle and sheep in Hungary from 2008 to 2009, proposed as members of a new genus in the family Picornaviridae The fecal viral flora of wild rodents Identification and complete genome characterization of a novel picornavirus in turkey (Meleagris gallopavo) Natural occurrence and characterization of two internal ribosome entry site elements in a novel virus, canine picodicistrovirus, in the picornavirus-like superfamily Bat guano virome: predominance of dietary viruses from insects and plants plus novel mammalian viruses The fecal virome of pigs on a high-density farm Southeast Asian foot-and-mouth disease viruses in Eastern Asia Complete genome sequence of a novel duck hepatitis A virus discovered in southern China Novel picornavirus in domesticated common quail (Coturnix coturnix) in Hungary Cleavage site analysis in picornaviral polyproteins: discovering cellular targets by neural networks CDD: conserved domains and protein three-dimensional structure The 2A proteins of three diverse picornaviruses are related to each other and to the H-rev107 family of proteins involved in the control of cell proliferation Cysteine proteases of positive strand RNA viruses and chymotrypsin-like serine proteases. A distinct protein superfamily with a common structural fold Comparative analysis of six genome sequences of three novel picornaviruses, turdiviruses 1, 2 and 3, in dead wild birds, and proposal of two novel genera, Orthoturdivirus and Paraturdivirus, in the family Picornaviridae Rotavirus proteins: structure and assembly Rotavirus genome replication and morphogenesis: role of the viroplasm Rotavirus infection: an update on management and prevention Epidemiology of rotavirus infection in broiler chickens: recognition of four serogroups Sequence analysis of the VP6-encoding genome segment of avian group F and G rotaviruses Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG) Conservation of a potential metal binding motif despite extensive sequence diversity in the rotavirus nonstructural protein NS53 Analysis of rotavirus species diversity and evolution including the newly determined full-length genome sequences of rotavirus F and G Complete nucleotide sequence of a group A avian rotavirus genome and a comparison with its counterparts of mammalian rotaviruses Nucleotide and amino acid sequence analysis of the rotavirus nonstructural RNA-binding protein NS35 The first complete genome sequence of a chicken group A rotavirus indicates independent evolution of mammalian and avian strains High variety of known and new RNA and DNA viruses of diverse origins in untreated sewage Metagenomic analysis of the turkey gut RNA virus community Phylogenetic and gene-centric metagenomics of the canine intestinal microbiome reveals similarities with humans and mice Analysis of the virus population present in equine faeces indicates the presence of hundreds of uncharacterized virus genomes A retrospective study of circovirus infection in pigeons: nine cases (1986-1993) Quantification of pigeon circovirus in serum, blood, semen and different tissues of naturally infected pigeons using a real-time polymerase chain reaction Epidemiological investigation of selected pigeon viral infections in Poland Prevalence of pigeon circovirus infections in feral pigeons in Ljubljana Viral infections in pigeons Partial genome sequence analysis of parvoviruses associated with enteric disease in poultry Prevalence of parvoviruses in commercial turkey flocks Viral infections in goose flocks in Poland High prevalence of turkey parvovirus in turkey flocks from Hungary experiencing enteric disease syndromes Simian immunodeficiency virus infects follicular helper CD4 T cells in lymphoid tissues during pathogenic infection of pigtail macaques Two novel parvoviruses in frugivorous New and Old World bats Letter: Evidence for viral gastroenteritis Isolation and characterization of rotavirus from feral pigeon in mammalian cell cultures