key: cord-1031327-4vb2boo6
authors: Nagai, Makoto; Omatsu, Tsutomu; Aoki, Hiroshi; Kaku, Yoshihiro; Belsham, Graham J.; Haga, Kei; Naoi, Yuki; Sano, Kaori; Umetsu, Moeko; Shiokawa, Mai; Tsuchiaka, Shinobu; Furuya, Tetsuya; Okazaki, Sachiko; Katayama, Yukie; Oba, Mami; Shirai, Junsuke; Katayama, Kazuhiko; Mizutani, Tetsuya
title: Identification and complete genome analysis of a novel bovine picornavirus in Japan
date: 2015-12-02
journal: Virus Res
DOI: 10.1016/j.virusres.2015.08.001
sha: 1963efabf6787b33386d7589664bd108ad7f01ba
doc_id: 1031327
cord_uid: 4vb2boo6

We identified novel viruses in feces from cattle with diarrhea collected in 2009 in Hokkaido Prefecture, Japan, by using a metagenomics approach and determined the (near) complete sequences of the virus. Sequence analyses revealed that they had a standard picornavirus genome organization, i.e. 5′ untranslated region (UTR) - L- P1 (VP4- VP3- VP2- VP1) - P2 (2A- 2B- 2C) - P3 (3A- 3B- 3C-3D) - 3′UTR- poly(A). They are closely related to other unclassified Chinese picornaviruses; bat picornaviruses group 1–3, feline picornavirus, and canine picornavirus, sharing 45.4–51.4% (P1), 38.0–44.9% (P2), and 49.6–53.3% (P3) amino acid identities, respectively. The phylogenetic analyses and detailed genome characterization showed that they, together with the unclassified Chinese picornaviruses, grouped as a cluster for the P1, 2C, 3CD and VP1 coding regions. These viruses had conserved features (e.g. predicted protein cleavage sites, presence of a leader protein, 2A, 2C, 3C, and 3D functional domains), suggesting they have a common ancestor. Reverse-transcription-PCR assays, using specific primers designed from the 5′UTR sequence of these viruses, showed that 23.0% (20/87) of fecal samples from cattle with diarrhea were positive, indicating the prevalence of these picornavirus in the Japanese cattle population in Hokkaido Prefecture. However, further studies are needed to investigate the pathogenic potential and etiological role of these viruses in cattle.

Viruses in the family Picornaviridae are small, icosahedral and nonenveloped, with a polyadenylated, single-stranded, positivesense RNA genome and these viruses are currently divided into 29 officially recognized genera: Aphthovirus, Aquamavirus, Avihepatovirus, Avisivirus, Cardiovirus, Cosavirus, Dicipivirus, Enterovirus, Erbovirus, Gallivirus, Hepatovirus, Hunnivirus, Kobuvirus, Kunsagivirus, Megrivirus, Mischivirus, Mosavirus, Oscivirus, Parechovirus, Pasivirus, Passerivirus, Rosavirus, Sakobuvirus, Salivirus, Sapelovirus, Senecavirus, Sicinivirus, Teschovirus and Tremovirus (Knowles et al., 2012; Adams et al., 2013) plus several candidate genera. Recently, a number of novel picornaviruses were discovered from different hosts Pankovics et al., 2012; Reuter et al., 2012; Sauvage et al., 2012; Boros et al., 2013; Liao et al., 2014) and the list of picornaviruses is rapidly expanding (http:// www.picornaviridae.com). In the Asia-Pacific region, several new picornaviruses, including three bat picornaviruses, one feline picornavirus, and one canine picornavirus were discovered in Hong Kong (Lau et al., 2011 Woo et al., 2012) . However, their clinical relevance remains to be elucidated. Picornaviruses can cause various symptoms including diarrhea and they are widely distributed in animals (Whitton et al., 2005) . In the course of surveillance for unexplained bovine diarrhea, using a metagenomics approach, we identified some novel picornaviruses. The present study describes the identification of these new picornaviruses (termed bovine Japanese picornaviruses) from fecal samples of calves with diarrhea in Hokkaido Prefecture, Japan and the genomic characterization of these viruses.

Forty fecal samples obtained from cattle with diarrhea in Hokkaido Prefecture in 2009 were analyzed using deep sequencing. The most common bovine viral pathogens inducing diarrhea, including group A, B, and C rotaviruses, bovine torovirus, and bovine coronavirus had been ruled out as the cause of the diarrhea by using reverse transcription (RT)-PCR diagnostic assays (20%; 10 of 50 samples). Samples were collected directly from the rectums to prevent contamination from the environment or between samples. Fecal samples were diluted 1:9 (w/v) in sterile phosphate-buffered saline with 1% antibiotics, and centrifuged at 10,000 × g for 10 min followed by filtration through a 0.22 m filter. The filtrates were stored at −80 • C until used for genomic analysis and virus isolation.

Viral RNAs were extracted from the filtrates using TRIzol ® LS Reagent (Life technologies, Carlsbad, CA, USA) and treated with DNase I (Takara Bio, Shiga, Japan). cDNA libraries for deep sequencing were constructed from RNA using the NEBNext Ultra Directional RNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, USA) according to the manufacturer's guidelines. Whole genome sequencing was determined by deep sequencing and the rapid amplification cDNA end method (RACE) (Roche Diagnostics GmbH, Mannheim, Germany) and 5 -Full RACE Core Set (Takara Bio). After assessing the library quality and quantity on a Bioanalyzer ® and Qubit ® 2.0 Fluorimeter (Invitrogen, Carlsbad, CA, USA), deep sequencing was conducted by using a MiSeq benchtop sequencer (Illumina, San Diego, CA, USA) using 51 nucleotide (nt) single reads. The runs were performed three times. Sequence data analysis was performed using the MiSeq Reporter v1 (Illumina) to generate FASTQ formatted sequence data. Trimmed reads were assembled into contigs using de novo assembly in CLC Genomics Workbench 6.0 (CLC bio, Aarhus, Denmark).

Viral sequence searches for assembled contigs were performed with the BLAST program in the National Center for Biotechnology Information non-redundant nucleotide database. Nucleotide sequences were aligned using ClustalW (Thompson et al., 1997) , and phylogenetic analysis was performed by the maximum likelihood method using MEGA5.22 (Tamura et al., 2011) . The trees were statistically supported by bootstrapping with 1000 replicates (Felsenstein, 1985) . Sequence similarity calculations were performed by using CLC Genomics Workbench 6.0 (CLC bio). Hypothetical polyprotein cleavage sites in the picornavirus were predicted using the NetPicoRNA program (Blom et al., 1996) . Secondary structure elements within the 5 UTR were predicted using Mfold (Zuker, 2003) .

Step RT-PCR Kit Ver.2 (TaKaRa Bio) using RNA extracted from fecal samples. The primers F (5 -CTT TTT CCC CCT CTT GYA AC -3 ) and R (5 -TTA GCC GCA TTC AGG GKC CTG GAG -3 ) were designed using the nucleotide sequence from the 5 UTR of the bovine Japanese picornaviruses. The production and amplification of the cDNA fragment was achieved under the following conditions: 50 • C for 30 min and 94 • C for 2 min, followed by 40 cycles of 94 • C for 30 s, 55 • C for 30 s, and 72 • C for 1 min, and a final extension at 72 • C for 5 min. The RT-PCR products were electrophoresed on a 2% agarose gel and purified.

The nucleotide sequences of the bovine Japanese picornavirus (Bo-11-39/2009/JPN; Bo-11-39, Bo-12-3/2009/JPN; Bo-12-3, Bo-12-7/2009/JPN; Bo-12-7, Bo-12-11/2009/JPN; Bo-12-11, Bo-12-38/2009/JPN; Bo-12-38) were deposited in the DNA Data Bank of Japan, the DDBJ/EMBL/GenBank, under accession numbers LC006971 and LC036579 to LC036582.

In the course of surveillance for unexplained bovine diarrhea using a metagenomics approach, we found 18 large sequence contigs that had a length of approximately 6000-8000 nucleotides (nt) from 14 samples. Ten, two and one sequence contigs showed high sequence homology to published bovine kobuvirus, bovine astrovirus and enterovirus F genomes respectively using nucleotide BLAST searches. Five of 18 sequence contigs, that had a length of approximately 7500 nt, were derived from five calves with diarrhea kept on five dairy farms and exhibited similarity to nucleotide sequences of other unclassified picornavirus genomes. Three of the five samples had more than two contigs (two samples: unclassified picornavirus and bovine kobuvirus, one sample; unclassified picornavirus, bovine kobuvirus and bovine astrovirus). Since the nucleotide sequence homology of the five contigs to other picornavirus genomes was rather low (<50%), we performed further analyses of these five contigs which were termed Bo-11-39, Bo-12-3, Bo-12-7, Bo-12-11, and Bo-12-38.

The complete genome length of Bo-11-39 was 7570 nt with G + C content of 0.41, excluding the poly(A) tail. We determined nearly complete genome sequences (including the complete open reading frames) of Bo-12-3, Bo-12-7, Bo-12-11 and Bo-12-38 but could not completely determine the sequence due to insufficient amount of samples. The genome organization of each is similar to that of other picornaviruses. They include a single large open reading frame (ORF) containing 6915 nt (Bo-11-38), or 7023 nt (Bo-12-3, 7, 11 and 38), which encodes a polyprotein precursor of 2305 amino acids (aa) (Bo-11-39), or 2341 aa (Bo-12-3, 7, 11 and 38), flanked by the 5'UTR and the 3'UTR with a poly(A) tail. Cleavage sites within the polyprotein were predicted by NetPicoRNA analysis and alignments with other picornaviruses. The junction sites appeared to be mainly Gln(Q)/Gly(G), except for predicted cleavage sites of the L/VP4 (Arg(R)/Gly(G), Bo-11-38 and Thr(T)/Gly(G), Bo-12-3, 7, 11 and 38), VP4/VP2 (Lys(K)/Ser(S)), and VP1/2A (Asp(D)/ Gly(G), Bo-12-3, 7, 11 and 38), respectively (Fig. 1) . The coding region consists of a leader protein (L), a structural protein region P1 (833 aa), and the non-structural protein regions P2 (645 aa in Bo-11-39 and 654 aa in Bo-12-3, 7, 11 and 38), and P3 (768 aa in Bo-11-39 and 793 aa in Bo-12-3, 7, 11 and 38). Pair-wise alignment of the amino acid sequences of Bo-11-39 and Bo-12-3, 7, 11 and 38 with those of other picornaviruses revealed that they shared higher sequence identities with bovine enterovirus 1 1) . These viruses are clearly very different from all previously described picornaviruses.

The complete 5'UTR of Bo-11-39 and near complete 5 UTR of Bo-12-3, 7, 11 and 38 had 38.5-41.7% and 38.5-41.9% identities to those of bat picornavirus group 3 and canine picornavirus, which have a type I internal ribosomal entry site (IRES). In contrast, the only 20.1-24.2% nt identities were found with the bat picornavirus group 1 and 2, and feline picornavirus sequences, which have type IV IRES elements (Table 2) . Within the picornavirus family it is now known that there are four classes of IRES element which vary in size and properties (Belsham, 2009) . Secondary structure prediction of the 5 -UTR of Bo-11-39 using the Mfold program revealed a stable stem-loop structure at the extreme 5 end of the genome, many picornaviruses have such a feature ( Fig. 2(a) ). Bat picornavirus group 3 has a section of the 5 -UTR with high nucleotide identity to the corresponding region of poliovirus (Lau et al., 2011) . Similarly, the 5 -UTR of Bo-11-39 and Bo-12-3, 7, 11 and 38 also has a region (position from nt 353 to nt 527 in Bo-11-39 and nt 317 to nt 490 in Bo-12-3, 7, 11 and 38) with high similarity to that of human enteroviruses (76-81%) including poliovirus (75-81%). A potential structures for domain V of the Bo-11-39 IRES, which by analogy to the poliovirus IRES, should bind to eIF4G (de Breyne et al., 2009 ) is shown in Fig. 2(b) . Both the base of the domain and the apex were extremely similar to other type I IRES elements and the domain is followed by a pyrimidine rich tract and then a short stem loop (domain VI), which contains an AUG and this is probably the site at which ribosomes start scanning although protein synthesis starts at the next AUG (at nt 554) ( Fig. 2(c) ). We predict that the translation initiation sites of Bo-11-39 and Bo-12-3, 7, 11 and 38 are at position nt 554 and nt 513, respectively, which is contained in the sequence GCCA 554 UGU in Bo-11-39 ( Fig. 2(c) ) and AACA 513 UGC in Bo-12-3, 7, 11 and 38. The initiation codon is not in an optimal Kozak context (Kozak, 1986 ) RNNAUGG with one nucleotide substitution (from G to U after AUG (RNNAUGU) in Bo-11-39 and G to C after AUG (RNNAUGC) in Bo-12-3, 7, 11 and 38). The 3 -UTR of bovine Japanese picornavirus was not closely related to any nucleotide sequences in the DDBJ/EMBL/GenBank database. There is less than 22.8% sequence identity within the 3 UTR between them and the unclassified Chinese picornaviruses.

For the coding regions, the predicted leader (L) protein of Bo-11-39 and Bo-12-3, 7, 11 and 38 are 59 and 61 aa long, respectively, and exhibit very low aa identities (less than 11.6%) to those of other representative picornaviruses, including unclassified Chinese picornaviruses. These L proteins are very short compared to viruses within the Aphthovirus and Erbovirus genera. The putative zing finger motifs which play a role in viral genome translation contained in the Cardiovirus genus were not found in the L proteins of the new bovine viruses (Chen et al., 1995; Dvorak et al., 2001) . The function of these bovine Japanese picornavirus L proteins remains to be identified.

The P1 region of the genomes of Bo-11-39 and Bo-12-3, 7, 11 and 38 encode the capsid proteins VP4, VP2, VP3, and VP1. Myristic acid is covalently linked to a glycine at the amino terminus of VP4 of most picornaviruses (Chow et al., 1987) . A potential myristoylation motif (GxxxT/S) is present in VP4 of Bo-11-39 (G 62 vsnS) and Bo-12-3, 7, 11 and 38 (G 62 ltvS) (Fig. 1) .

The P2 region of the genomes of Bo-11-39 and Bo-12-3, 7, 11 and 38 encode the nonstructural proteins 2A, 2B, and 2C. The conserved GxCG motif, which is considered to form part of the active site of the protease (Gorbalenya et al., 1989a) , in the 2A protein of Bo-11-39 (G 1056 wCG) and Bo-12-3, 7, 11 and 38 (G 1060 wCG) suggests that the 2A protein of these viruses cleaves at its own N-terminus, i.e. the VP1/2A junction. The conserved H-box/NC motif involved in cell proliferation control (Hughes and Stanway, 2000) was not present in the 2A. Similar to all the other picornaviruses, the NTPbinding motif GXXGXGKS (Gorbalenya et al., 1989b) was found in the 2C of Bo-11-39 and Bo-12-3, 7, 11 and 38 as G 1337 gpGsGKS and G 1348 apGtGKS, respectively. The DDLXQ motif is believed to be important for putative helicase activity (Gorbalenya et al., 1989a (Gorbalenya et al., , 1990 ) but the third amino acid of this motif in the 2C of Bo-11-39 and Bo-12-3, 7, 11 and 38 was substituted by Val (D 1385 DVGQ and D 1396 DVGQ). The P3 genome region of Bo-11-39 Bo-12-3, 7, 11 and 38 encodes the non-structural proteins 3A, 3B, 3C pro , and 3D pol . The conserved catalytic triad of H-D/E-C was seen in 3C pro of Bo-11-39 (H 1700 -D 1737 -C 1794 ) and Bo-12-3, 7, 11 and 38 (H 1736 D 1773 -C 1830 ). The conserved GxCG motif was also found in 3C pro of Bo- 11-39 (G 1803 qCG) and Bo-12-3, 7, 11 and 38 (G 1839 qCG) , as with other picornaviruses. Two highly conserved motifs of RNAdependent RNA polymerases (Kamer and Argos, 1984 ), K 2001 DELR, G 2128 GmPSG, Y 2170 GDD, and F 2217 LKR (Bo-11-39) and K 2037 DELR, G 2164 GmPSG, Y 2206 GDD, and F 2253 , were also present in the 3D pol of these viruses.

Phylogenetic trees were constructed from the complete amino acid cording sequences for P1, 2C, 3CD, and VP1 of bovine Japanese picornaviruses together with representative strains of the different genera and some unclassified picornaviruses. 7, 11 and 38 were closely related to each other and formed a distinct cluster in the 3CD and VP1 trees with high bootstrap support (98% and 79%). The bovine Japanese picornaviruses were related to enteroviruses, Sapeloviruses, and unclassified Chinese picornaviruses and formed an independent cluster with unclassified Chinese picornaviruses supported by high bootstrap value (76%-99%) in all four phylogenetic trees (Fig. 3) .

The prevalence of bovine Japanese picornaviruses in Hokkaido Prefecture in 2009 was studied using RT-PCR assays performed using specific primers designed from the 5 UTR sequences. Fecal samples were collected from 87 cattle with diarrhea (calves: 47; heifers: 14; cows: 26) from 72 farms. Using standard diagnostic assays, the most common bovine viral pathogens that induce diarrhea were ruled out as being responsible. However, amplicons of the expected size (326 base pairs) consistent with the presence of the bovine Japanese picornaviruses were obtained from twenty of these samples (from 15 calves (32%), a heifer (7%) and four cows (15%)) from nineteen different farms.

In recent years, the advances in the metagenomic field have permitted the definition of new viral genomic sequences for clinical and epidemiological research (Beerenwinkel et al., 2012; Mokili et al., 2011; Rosario and Breitbart, 2011) . Deep sequencing technologies have facilitated virus discovery due to the potential to detect viruses having an unknown genetic background. In this study, we have discovered a new picornavirus from cattle using metagenomic technology. We tried to isolate viruses using primary bovine testis cells, primary bovine fetal muscles cells, the Madin-Darby Bovine Kidney cells, Vero cells, baby hamster kidney cells, cloned porcine kidney cells, and HRT-18 cells. Neither cytopathic effect nor the expected RT-PCR products from RNA extracted from cell culture supernatants were observed in any of the inoculated cells, indicating a lack of virus replication in these cell lines. However, using RNA-seq based strategy with DNase treatment before construction of cDNA libraries, we could obtain sufficient viral genomic sequences of RNA viruses from fecal samples not only for picornaviruses but also, in earlier studies, for other viruses such as rotavirus (Minami-Fukuda et al., 2013; Masuda et al., 2014; Nagai et al., 2015) .

Genomic characterization showed that these bovine Japanese picornaviruses had a genome organization similar to other picornaviruses, i.e. the presence of single large open reading frame flanked by UTRs. The 5 UTR of bovine Japanese picornavirus begins with a stem loop structure of about 50 nt, which is different from entero-/rhinoviruses that have a cloverleaf structure. However these viruses appear to have a type I IRES with a predicted domain IV (including a loop B), plus the domains V and VI including the eIF4G binding site and an AUG codon respectively, the latter is preceded by a polypyrimidine tract. Comparisons with representative strains of the different defined genera in the family Picornaviridae revealed that amino acid identities in the P1, P2, and P3 regions were up to 39.5, 30.5, and 45.9%, respectively. The Picornaviridae Study Group of the International Committee on Taxonomy of Viruses (http://www.picornastudygroup.com/ definitions/genus definition.htm) has proposed that different picornavirus genera are defined by amino acid identity in the P1, P2, and P3 regions being <40, <40, and <50%, respectively. Based on these criteria, bovine Japanese picornaviruses do not fall into the current official genera. On the other hand, bovine Japanese picornaviruses shared 45.4 to 51.4% (P1), 38.0 to 44.9% (P2), and 49.6 to 55.2% (P3) amino acid identities with certain unclassified Chinese picornaviruses; bat picornaviruses group 1-3 (Lau et al., 2011) , feline picornavirus , and canine picornavirus , respectively. Furthermore, phylogenetic analyses demonstrated that bovine Japanese picornaviruses formed a cluster with unclassified Chinese picornaviruses in P1, 2C, 3CD, and VP1 coding regions, respectively, suggesting that bovine Japanese picornaviruses may be classified together with the currently unclassified Chinese picornaviruses into a novel genus in the family Picornaviridae.

Detailed analyses showed that bovine Japanese picornaviruses possessed conserved characteristic motifs, similar to unclassified Chinese picornaviruses (e.g. protein cleavage sites, presence of a leader protein, 2A, 2C, 3C, and 3D functional domains), with minor variance. Genomic comparison of bovine Japanese picornaviruses with unclassified Chinese picornaviruses revealed that although amino acid identities in L, 2A, 2B, and 3A were low (<40%), amino acid identities in the other regions were higher and indicated that these viruses were moderately related to each other. Phylogenetic analysis of the VP1 coding regions indicated that bovine Japanese picornaviruses and unclassified Chinese picornaviruses formed a cluster with high bootstrap support, even though VP1 is the most diverse capsid protein of picornaviruses (Oberste et al., 1999; Rossmann et al., 1985) . These findings suggest that they have a common ancestor. Bovine Japanese picornaviruses seem most related to bat picornavirus group 3 and canine picornaviruses which also contain a type I IRES. RNA viruses have a remarkable potential to evolve due to high mutation and recombination rates (Cook et al., 2013) . Bovine Japanese picornaviruses may be derived from the unclassified Chinese picornavirus and have evolved independently in cattle in Japan.

Finally, RT-PCR assays identified the presence of RNA from bovine Japanese picornaviruses in 15 calves (32%), one heifer (7%), and four cows (15%) with diarrhea from 19 farms. These finding indicates that these viruses are circulating in Hokkaido Prefecture in Japan. The potential transmission of these viruses to other domestic animals or human should be clarified. Since we did not perform analyses using fecal samples of cattle without diarrhea, the association of bovine Japanese picornaviruses with bovine diarrhea is unclear. To clarify whether they are associated or causal for bovine diarrhea, further studies are needed.

B14 Human Rhinovirus B (K02121 )

CV-B1 Enterovirus B (M16560 ) [Enterovirus] ne picornavirus (JN572117 )

Sape lovirus] Simian sapelovirus (AY064708 )

Quail picornavirus (JN674502 )

Hunn ivirus (JQ941880 ) [Huun ivirus] Porcine tesc hovirus (AJ011380 )

Equ ine rhinitis B virus (X96871 )

Foo t-and-mouth virus (AF308157 )

Avian sapelovirus (AY563023 )

Mosavirus M-7 (KF958461 )

Galli virus A (JQ691613 ) [Galli virus] Pass erivirus A1 (GU182406 )

Osc ivirus A (GU182408 )

Sali virus (GQ179640 )

Di cipivirus] Rosavirus M-7 (JF973686 )

Hepatitis A virus 1 (M14707 )

B14 Human Rhinovirus B (K02121 )

Simian sapelovirus (AY064708 )

Cali fornia sea li on sapelovirus 1 (JN420368 )

Avian sapelovirus (AY563023 ) [Sape lovirus] Pigeo n picornavirus (FR727144 )

Quail picornavirus (JN674502 )

Mosavirus M-7 (KF958461 ) [Mosavirus] Hunn ivirus (JQ941880 ) [Huun ivirus] Porcine tesc hovirus (AJ011380 )

Equ ine rhinitis B virus (X96871 )

Enceph alomyocarditis virus 1 (M81861 ) [Cardiovirus] Seneca vall ey virus 1 (DQ641247 )

Foo t-and-mouth virus (AF308157 )

Galli virus A (JQ691613 ) [Galli virus] Pass erivirus A1 (GU182406 )

Osc ivirus A (GU182408 )

Sali virus (GQ179640 )

Rosavirus M-7 (JF973686 )

Aqu amavirus A (EU142040 )

Canine picornavirus

Avian sapelovirus (AY563023 )

B14 Human Rhinovirus B (K02121 )

Osc ivirus A (GU182408 )

Rosavirus M-7 (JF973686 )

Mosavirus M-7 (KF958461 ) [Mosavirus] Hunn ivirus (JQ941880 )

Tesc hovirus] Equ ine rhinitis B virus (X96871 ) [Erbo virus] Foo t-and-mouth virus (AF308157 )

Hepatitis A virus 1 (M14707 )

Swine/ K23 /2008/HUN Ente rovirus G (HQ702854 )

Unclass ified] Canine picornavirus (JN831356 ) [Unclass ified] picornavirus (FR727144 )

Sape lovirus] Foo t-and-mouth virus (AF308157 )

Avian enceph alomyeli tis virus 1 (AJ225173 ) [Tremovirus] Galli virus A (JQ691613 ) [Galli virus] Pass erivirus A1 (GU182406 )

Equ ine rhinitis B virus (X96871 ) [Erbo virus] Cadicivirus (JN819202 )

Rosavirus M-7 (JF973686 )

open box) with reference picornaviruses obtained from the DDBJ/EMBL/GenBank database. Phylogenetic trees were constructed using the maximum likelihood method in MEGA5.22 with bootstrap values (1000 replicates) above 70 being shown. The scale bar indicates amino acid substitutions per site. References Adams

Challenges and opportunities in estimating viral genetic diversity from next-generation sequencing data

Divergent picornavirus IRES elements

Cleavage site analysis in picornaviral polyproteins: discovering cellular targets by neural networks

Identification and complete genome characterization of a novel picornavirus in turkey (Meleagris gallopavo)

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

The leader peptide of Theiler's murine encephalomyelitis virus is a zinc-binding protein

Myristylation of picornavirus capsid protein VP4 and its structural significance

Novel virus discovery and genome reconstruction from field RNA samples reveals highly divergent viruses in dipteran hosts

Direct functional interaction of initiation factor eIF4G with type 1 internal ribosomal entry sites

Leader protein of encephalomyocarditis virus binds zinc, is phosphorylated during viral infection, and affects the efficiency of genome translation

Confidence limits on phylogenies: an approach using the bootstrap

Cysteine proteases of positive strand RNA viruses and chymotrypsin-like serine proteases. A distinct protein superfamily with a common structural fold

Two related superfamilies of putative helicases involved in replication, recombination, repair and expression of DNA and RNA genomes

A new superfamily of putative NTP-binding domains encoded by genomes of small DNA and RNA viruses

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

Primary structural comparison of RNA-dependent polymerases from plant, animal and bacterial viruses

Picornaviridae

Point mutations define a sequence flanking the AUG initiator codon that modulates translation by eukaryotic ribosomes

Complete genome analysis of three novel picornaviruses from diverse bat species

Identification of a novel feline picornavirus from the domestic cat

Genomic characterization of a novel picornavirus in Pekin ducks

Identification of novel bovine group A rotavirus G15P[14] strain from epizootic diarrhea of adult cows by de novo sequencing using a next-generation sequencer

Detection of bovine group a rotavirus using rapid antigen detection kits, rt-PCR and next-generation DNA sequencing

Metagenomics and future perspectives in virus discovery

Molecular evolution of the human enteroviruses: correlation of serotype with VP1 sequence and application to picornavirus classification

Novel picornavirus in domesticated common quail (Coturnix coturnix) in Hungary

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

Exploring the viral world through metagenomics

Structure of a human common cold virus and functional relationship to other picornaviruses

A member of a new Picornaviridae genus is shed in pig feces

MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods

The CLUSTAL X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools

Host and virus determinants of picornavirus pathogenesis and tropism

Complete genome sequence of a novel picornavirus, canine picornavirus, discovered in dogs

Mfold web server for nucleic acid folding and hybridization prediction

This work was supported by the Grants from the Ministry of Health, Labor and Welfare of Japan.