Email updates

Keep up to date with the latest news and content from Virology Journal and BioMed Central.

Open Access Research

Molecular epidemiology of Porcine torovirus (PToV) in Sichuan Province, China: 2011–2013

Lu Zhou1, Haoche Wei13, Yuancheng Zhou1, Zhiwen Xu12*, Ling Zhu12 and Jim Horne3

Author Affiliations

1 Animal Biotechnology Center, College of Veterinary Medicine, Ya’an, China

2 Key Laboratory of Animal Disease and Human Health, College of Veterinary Medicine, Sichuan Agricultural University, Ya’an, China

3 Faculty of Bio and Biochemistry, University of Bath, Bath spa, England

For all author emails, please log on.

Virology Journal 2014, 11:106  doi:10.1186/1743-422X-11-106


The electronic version of this article is the complete one and can be found online at: http://www.virologyj.com/content/11/1/106


Received:5 January 2014
Accepted:30 May 2014
Published:5 June 2014

© 2014 Zhou et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.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.

Abstract

Background

Porcine torovirus (PToV) is a member of the genus Torovirus which is responsible for gastrointestinal disease in both human beings and animals with particular prevalence in youth. Torovirus infections are generally asymptomatic, however, their presence may worsen disease consequences in concurrent infections with other enteric pathogens.

Methods

A total of 872 diarrheic fecal samples from pigs of different ages were collected from 12 districts of Sichuan Province in the southwest of China. RT-PCR was done with PToV S gene specific primers to detect the presence of PToV positive samples. M gene specific primers were used with the PToV positive samples and the genes were sequenced. A phylogenetic tree was constructed based on the M gene nucleotide sequences from the 19 selected novel Sichuan strains and 21 PToV and BToV M gene sequences from GenBank.

Results

A total of 331 (37.96%, 331/872) samples were found to be positive for PToV and the highest prevalence was observed in piglets aged from 1 to 3 weeks old. Through phylogenetic inference the 40 PToV M gene containing sequences were placed into two genotypes (I & II). The 19 novel Sichuan strains of genotype I showed strong correlations to two Korean gene sequences (GU-07-56-11 and GU-07-56-22). Amino-acid sequence analysis of the 40 PToV M gene strains revealed that the M gene protein was highly conserved.

Conclusions

This study uncovered the presence of PToV in Sichuan Province, and demonstrated the need for continuous surveillance PToV of epidemiology.

Keywords:
Porcine torovirus; Membrane gene; Phylogenetic analysis

Introduction

Porcine torovirus (PToV) is a member of the genus Torovirus, family Coronaviridae, and order Nidovirales[1]. Toroviruses have been detected in fecal samples from cows, horses, humans, pigs and turkeys suffering from diarrheal illness [2-7]. The first identified torovirus was isolated from a horse with diarrhea (equine torovirus, EToV, also called Berne virus) in 1972 [8].

To date, research on toroviruses has been limited because of difficulties associated with their propagation in vitro. In 1983, EToV was found to be able to adapt and grow in equine cells. As a consequence, much of what is known of torovirus at the molecular level has been based on EToV [9].

PToV particles were initially observed by electron microscopy from pig fecal samples. They were described as spherical, oval, elongated, or kidney-shaped enveloped viruses possessing a positive-sense, single-stranded, polyadenylated RNA genome of approximately 25–30 kb in length [10-12]. In 2013, a study reported that the genome was found to be 28,301 bp long [13], sharing 79% identity with bovine torovirus (BToV). It mainly consists of a replicase (20,906 bp) and structural genes: spike, S (4,722 bp); membrane, M (702 bp); hemagglutinin-esterase, HE (1,284 bp); and nucleocapsid, N (492 bp) [13-15]. Sequence data for PToV strains are based on the S, M, HE, and N genes; M is more conserved compared with other genes, sharing 98% amino acid sequence identity between the described PToVs M genes [13], which have only been reported in Spain and South Korea.

PToV infections have been detected in many countries, including the Netherlands, Canada, South Africa, the United States, Hungary, the UK and South Korea. Serum samples, obtained from farms in Spain and screened via the ELISA method indicated very high seroprevalence against PToV in some studies [16,17], while other studies showed low seroprevalence [18]. In this article, we investigate the infection and characteristics of PToV on the basis of nucleotide and amino acid sequencing and phylogenetic analysis of M genes, for PToV isolated from porcine fecal samples in Sichuan Province, China.

Methods and materials

A total of 872 fecal samples were collected from pigs with diarrhea across different districts of Sichuan Province during the winter from 2011 to 2013. Most of the sampled piglets were 1 to 3 weeks old, with the reminder being 3 to 7 weeks old, 7 to 11 weeks old or over 11 weeks old. Antibiotic treatment was invalid in all sampled pigs. Samples were analyzed based on collection time and geographical location in either the east or west sections of Sichuan Province. A portion of each sample was homogenized for RNA extraction, and the remaining samples were stored at −70°C.

According to the manufacturer’s instructions, total RNA was extracted from fecal homogenate using TRIzol reagent (Life Technologies, Gaithersburg, MD, USA). Viral cDNA was synthesized from RNA using reverse transcriptase (M-MLV, Takara, Kyoto, Japan) according to the manufacturer’s instructions, and stored at −20°C until used in RT-PCR reactions.

The cDNA was screened by PCR using the following method we had established. A 451-bp fragment of the conserved region of the S gene was amplified with the forward primer (5' - ACCCCTGCCTGAGGTTTCYTT - 3'), and reverse primer (5' – AGCACGACGTTGTCTRCGTGT - 3’). Amplification was carried out in PCR buffer containing 200 mM of each dNTP, 10 pmol of each primer, 1.0 U Taq DNA polymerase (Promega, Madison, WI, USA), and 1.5 mM MgCl2, in a total volume of 40 μl. PCR was performed at 94°C for 2 min, followed by 30 cycles of amplification (94°C for 30 s, 57°C for 30 s, and 72°C for 30 s), and a final extension of 72°C for 7 min. The PCR product was resolved using 1% agarose gel electrophoresis, stained with ethidium bromide (Invitrogen, Carlsbad, CA, USA), and visualized under ultraviolet light with the Bio-Rad gel imaging system (Hercules, CA, USA). All specimens were also tested for the presence of porcine epidemic diarrhea virus (PEDV), transmissible gastroenteritis coronavirus (TGEV) and group A rotavirus (RVA) in terms of the methods described in previous studies [19,20]. The PCR and Nested PCR specific primers used are listed in Table 1.

Table 1. Oligonucleotide primers used for the detection of PEDV, TGEV and RVA in fecal samples obtained from pigs with diarrhea

PToV genomic cDNA was obtained from 19 positive samples as described above. The 702 bp fragment of the complete M gene was amplified with primers M1 (5' - ATGTTTGATACAAATTTTTGGCCTT - 3') and M2 (5' – CTACTCAAACTTAACA CTTGACAACTGC - 3'). PCR amplification was carried out as described above, and the PCR products were visualized using 1% agarose gel electrophoresis under ultraviolet light.

The PCR products were gel-purified using a Gel Extraction Kit (Tiangen Biotech, Beijing, China). The purified target fragments were ligated into a linear vector pMD19-T (Takara, Dalian, China), and the recombinant plasmids were transformed into Escherichia coli DH5α competent cells (Invitrogen). The identity of the constructs was confirmed by sequencing (Invitrogen).

Genetic distance was initially measured by pairwise comparisons of nucleotide sequences to reference PToV sequences available in GenBank using the Basic Local Alignment Search Tool (BLAST) (http://blast.ncbi.nlm.nih.gov/ webcite). Multiple alignments were achieved using the Clustal_W method of the MegAlign 5.01 program (DNASTAR Inc., Madison, WI, USA) [21]. The phylogenetic tree was constructed using the Neighbor-Joining method on Molecular Evolutionary Genetics Analysis (MEGA) software version 5.0 [22] and carried out using the Kimura 2-parameter model [23]. The transition/transversion bias (R) and substitution rates were estimated using MEGA 5.0. The substitution pattern and rates were estimated using the Kimura 2-parameter model, with nucleotide frequencies of A = 25.00%, T/U = 25.00%, C = 25.00%, and G = 25.00% [23]. For estimating ML- maximum likelihood values, a user-specified topology was used. The maximum log likelihood for this computation was −1680.662 and −2731.225. Codon positions included were 1st + 2nd + 3rd + Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 702 positions in the final dataset. Evolutionary analyses were conducted in MEGA5 [24]. All analyses were based on PToV M gene sequences.

The novel PToV M nucleotide sequences have been submitted to NCBI GenBank and assigned accession numbers KF727566-KF727584. Other toroviruses full-length sequences for M gene were obtained from GenBank (http://www.ncbi.nlm.nih.gov/nuccore webcite).

Results

The samples from 12 districts in Sichuan Province, China are shown in Figure 1. Of these samples, 331 (37.96%, 331/872) were positive for PToV, as detected by RT-PCR. PToV infection rates were greater than 37.96% for Suining, Deyang, Yibin, Chengdu, Mianyang, and Meishan districts. PToV infection rates of 36.52% (42/115), 5.77% (3/52), 9.52% (2/21), 27.78% (10/36), 4% (1/25), and 10.26% (4/39) were detected for Ya’an, Leshan, Dazhou, Zigong, Luzhou and Ziyang, respectively. For other tested enteric pathogens, the positive rates of PEDV, TGEV, RVA were 62.5% (545/872), 12.5% (109/872) and 66.9% (583/872), respectively, while the PToV Co-infected with RVA, TGEV and PEDV positive rate was 4.1% (36/872). And the details of PToV, RVA, TGEV and PEDV infection and co-infection rates in different districts of Sichuan Province were in Table 2. The PToV positive rates at different growth stages of pigs revealed that 1 to 3 week old piglets have the highest PToV infection rate 42.47% (296/697), with PToV infection rates of 25.53% (24/94), 16.07% (9/56) and 8% (2/25) in 3 to 7 weeks, 7 to 11 weeks and over 11 weeks, respectively. The complete M sequence of 42 of the 331 PToV-positive strains were amplified and sequenced. Of the 42 PToV strains, 19 strains were selected, based on geographic factors and their isolation from different farms in the same district.

thumbnailFigure 1. Geographical locations of samples collected in this study. The green parts indicate the districts of Sichuan Province where the porcine fecal samples were collected.

Table 2. Summary of enteric pathogens present in the fecal samples obtained from pigs with diarrhea

The complete M gene sequences of the 19 strains were compared with each other, and against reference sequences for 21 strains of toroviruses in the GenBank database (Table 3). As shown in Figure 2, phylogenetic and homology studies at the nucleotide (702 bp) level demonstrated that PToV strains belonged to two genotypes (I and II), separated with bootstrap values of 90 and 60, respectively. All 19 novel PToV strains within genotype I were further divided into three genotypes (A, B and C). With only two exceptions (GU-07-56-11 and GU-07-56-22), all other representative strains from NCBI belonged to genotype II. Strains GU-07-56-11 and GU-07-56-22 shared high homology with four novel PToV strains, SC-2012-L1, SC-2012-L2, SC-2012-K and SC-2012-J, and were classified as genotype B. Nucleotide and the deduced amino acid sequence homologies for the 19 Chinese PToV M sequences identified in this study were between 97.3–99.9% and 99.1–99.6% respectively. The Chinese PToV strains were 92.7–99.9% and 97.4–99.6% identical to the PToV strains in other countries while only 80.1%–82.8% and 89.3%–91.0% to those of bovine torovirus.

Table 3. GenBank accession numbers of Chinese PToVs and the reference toroviruses used in molecular and phylogenetic analyses”

thumbnailFigure 2. Phylogenetic analysis based on the complete nucleotide sequence (702 bp) of gene M for novel PToV isolated from Sichuan, China, and reference viruses (GenBank). The reference PToV M nucleotide sequences were obtained from GenBank: Korea strain (accession no. GU181240.1, GU181241.1, GU181242.1, GU181243.1,GU181244.1, GU181245.1 ,GU181246.1, GU181247.1, GU181248.1, GU181249.1, GU181250.1); Netherlands strain(accession no. AJ575374.1, AJ575376.1, AJ575377.1, AJ575369.1, AJ575368.1, AJ575371.1); Canada strain (accession no. AF076621.1, AY427798.1); The 19 strains analyzed in this study are named by SC-.

The 19 novel M nucleotide sequences and their deduced amino acid sequences were aligned. The estimated transition/transversion bias (R) was 1.15 among the 19 PToV strains from Sichuan in China, and 1.84 among global PToV strains. The Sichuan PToV nucleotide sequences presented 54 mutations (Table 4); 26 of which were commonly mutated sites (at nucleotide positions 17, 18, 27, 39, 72, 121, 276, 288, 306, 495, 512, 513, 537, 552, 561, 579, 618, 633, 636, 652, 654, 657, 663, 666, 678, 690) in the gene M nucleotide sequence. The Sichuan strains had three common amino-acid substitutions(Y6 → F, N171 → S, N181 → S), while the rest were silent mutations.

Table 4. Comparison of nucleotide and deduced amino acid sequences of M genes between 19 Chinese PToVs and reference PToV (GenBank accession No. GU181244)

Discussion

Based on the availability of an effective and feasible RT-PCR diagnostic tool for PToV, a survey of PToV was performed across 12 districts in Sichuan Province Southwestern China. This RT-PCR methodology had been previously proved to be sensitive and specific for PToV [25]. RT-PCR displayed that PToV infections were widespread in large-scale breeding bases and on rural farms in Sichuan. The overall positive rate was 37.96%, and six districts, Suining, Deyang, Yibin, Chengdu, Mianyang and Meishan, had higher infection rates than this. These results may be due to the increased urban agglomeration and greater levels of prosperity within these districts leading to greater movement of, and poorer storage and feeding of livestock with in these areas. According to the PToV infection rates at different growth stages, piglets of 1 to 3 weeks of age were more vulnerable to PToV than the others. Nonetheless, the result is likely to be an incomplete reflection the total number of infections in Sichuan province. The reason for this is that the highest prevalence and greatest anti-PToV Ig G titers were observed in adults [17]. In contrast, the lowest prevalence and lowest anti-PToV Ig G levels were found among 3-week-old piglets [17]. Both ELISA reactivity values and rates of seroprevalence rose markedly in piglets from 3 to 11 weeks [17]. These results suggested that RT-PCR cannot reflect the infections completely due to swine developing their own immune response to the virus resulting in an effective cure in some infected adults. Consequently, further immunological methods should be established to develop a more robust serological surveillance of PToV.

According to the co-infection rates, we can deduce PToV-infected piglets were probably always co-infected with PEDV, RVA and TGEV. Meanwhile, its interesting to note that diarrhea samples from Luzhou tested positive for PToV alone, agreeing with a previous study showing the same results when a survey for enteric pathogens in diarrheic pigs was carried out [26]. However, we cannot make any conclusions about the association between PToV and other enteric pathogens, because of the difficulty of growing the virus in cultured cells.

The comparison of the nucleotide and deduced amino acid sequences of the full-length PToV M nucleotide sequence confirmed that the Sichuan PToV strains were highly conserved for the M gene and most mutations were silent. The phylogenetic tree showed that the 19 Sichuan strains and 2 Korean strains belonged to genotype I, while all other strains resided in genotype II. One special case was that 4 PToV strains within genotype I from the study were classified with 2 strains from Korea in genotype B. Therefore, we speculated that the 4 Chinese PToV strains in subroup B may have originated in Korea, but further study on the origin of circulating PToV in China should be assessed in the future.

Torovirus is widely depicted as a pathogen which is responsible for diarrheal illnesses in animals and human beings [27-29]. In recent years, studies conducted on pigs farmed in Spain, reported that the seroprevalence against PToV in piglets over 11 weeks of age was > 99% indicating that PToV was endemic in piglets in Spain [17]. Epidemiological data on PToV has also been reported from in other countries [3,16], suggesting that it may be circulating in swine populations globally, and with high prevalence. Most PToV infections appear to be subclinical [5], however, there is a potential threat that a much more serious outcome will appear in concurrent infections with other enteric pathogens.

Conclusion

This is a molecular epidemiological investigation of PToV in Sichuan Province, China. A PToV infection rate of 37.96% (331/872) supported the infectious situation of PToV in southwest China. In addition, nucleotide and amino acid sequence analyses confirmed that the PToV M protein was highly conserved. Since the relationship between PToV infection and diarrhea remains unclear, the virus is still a potential threat for us. Further studies to reveal the epidemiological status of PToV infection in China are essential, and the research should focus on the epidemiology and pathogenic potential of PoTV. Finally, the rapid and continual changes within the torovirus genome highlight the need for persistent surveillance and regular reviews of phylogeny so as to better understand the evolution of this pathogen.

Consent

Written informed consent was obtained from patients for publication of the data in this manuscript and any accompanying images. A copy of the written consent is available for review by the Editor-in-Chief of this journal.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: LZ, ZX. Performed the experiments: LZ, YZ. Analyzed the data: LZ, YZ, HW. Contributed reagents/ materials/analysis tools: LZ, HW, YZ, LZ. Wrote the paper: LZ. Checked the grammar: JH. All authors read and approved the final manuscript.

Acknowledgments

We thank Jim Horne for checking and correcting the grammar. This study was supported by Sichuan Province Science and Technology Support Project (2012NZ0001) and Fostering Excellent Master Dissertation Foundation of Sichuan Agricultural University (Grant No.04310736), program for New Century Excellent Talents of Ministry of Education of China (NCET 11–1059).

References

  1. Cavanagh D: Nidovirales: a new order comprising Coronaviridae and Arteriviridae.

    Arch Virol 1997, 142:629-633. PubMed Abstract | Publisher Full Text OpenURL

  2. Ali A, Reynolds D: Characterization of the stunting syndrome agent: relatedness to known viruses.

    Avian diseases 2000, 44:45-50. PubMed Abstract | Publisher Full Text OpenURL

  3. Shin D-J, Park S-I, Jeong Y-J, Hosmillo M, Kim H-H, Kim H-J, Kwon H-J, Kang M-I, Park S-J, Cho K-O: Detection and molecular characterization of porcine toroviruses in Korea.

    Arch Virol 2010, 155:417-422. PubMed Abstract | Publisher Full Text OpenURL

  4. Duckmanton L, Carman S, Nagy É, Petric M: Detection of bovine torovirus in fecal specimens of calves with diarrhea from Ontario farms.

    J Clin Microbiol 1998, 36:1266-1270. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Kroneman A, Cornelissen L, Horzinek M, De Groot R, Egberink H: Identification and characterization of a porcine torovirus.

    J Virol 1998, 72:3507-3511. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Duckmanton L, Luan B, Devenish J, Tellier R, Petric M: Characterization of torovirus from human fecal specimens.

    Virology 1997, 239:158-168. PubMed Abstract | Publisher Full Text OpenURL

  7. Beards D, Green B, Flewett T: Preliminary characterisation of Torovirus‒like particles of humans: comparison with Berne virus of horses and Breda virus of calves.

    J Med Virol 1986, 20:67-78. PubMed Abstract | Publisher Full Text OpenURL

  8. Weiss M, Steck F, Horzinek MC: Purification and partial characterization of a new enveloped RNA virus (Berne virus).

    J Gen Virol 1983, 64(Pt 9):1849-1858. PubMed Abstract | Publisher Full Text OpenURL

  9. Brown D, Beards G, Flewett T: Detection of Breda virus antigen and antibody in humans and animals by enzyme immunoassay.

    J Clin Microbiol 1987, 25:637-640. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Horzinek MC: Molecular evolution of corona- and toroviruses.

    Adv Exp Med Biol 1999, 473:61-72. PubMed Abstract | Publisher Full Text OpenURL

  11. Cornelissen LA, Wierda CM, van der Meer FJ, Herrewegh AA, Horzinek MC, Egberink HF, de Groot RJ: Hemagglutinin-esterase, a novel structural protein of torovirus.

    J Virol 1997, 71:5277-5286. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Snijder EJ, Horzinek MC: Toroviruses: replication, evolution and comparison with other members of the coronavirus-like superfamily.

    J Gen Virol 1993, 74(Pt 11):2305-2316. PubMed Abstract | Publisher Full Text OpenURL

  13. Sun H, Lan D, Lu L, Chen M, Wang C, Hua X: Molecular characterization and phylogenetic analysis of the genome of porcine torovirus.

    Arch Virol 2014, 159:773-778. PubMed Abstract | Publisher Full Text OpenURL

  14. Cong Y, Zarlenga DS, Richt JA, Wang X, Wang Y, Suo S, Wang J, Ren Y, Ren X: Evolution and homologous recombination of the hemagglutinin-esterase gene sequences from porcine torovirus.

    Virus Genes 2013, 47:66-74. PubMed Abstract | Publisher Full Text OpenURL

  15. De Vries AA, Horzinek MC, Rottier PJ, de Groot RJ: The genome organization of the Nidovirales: similarities and differences between arteri-, toro-, and coronaviruses.

    Semin Virol 1997, 8:33-47. Publisher Full Text OpenURL

  16. Pignatelli J, Grau-Roma L, Jimenez M, Segales J, Rodriguez D: Longitudinal serological and virological study on porcine torovirus (PToV) in piglets from Spanish farms.

    Vet Microbiol 2010, 146:260-268. PubMed Abstract | Publisher Full Text OpenURL

  17. Alonso-Padilla J, Pignatelli J, Simon-Grife M, Plazuelo S, Casal J, Rodriguez D: Seroprevalence of porcine torovirus (PToV) in Spanish farms.

    BMC Res Notes 2012, 5:675. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Hosmillo MD, Jeong Y-J, Kim H-J, Collantes TM, Alfajaro MM, Park J-G, Kim H-H, Kwon H-J, Park S-J, Kang M-I: Development of universal SYBR Green real-time RT-PCR for the rapid detection and quantitation of bovine and porcine toroviruses.

    J Virol Methods 2010, 168:212-217. PubMed Abstract | Publisher Full Text OpenURL

  19. Kim SY, Song DS, Park BK: Differential detection of transmissible gastroenteritis virus and porcine epidemic diarrhea virus by duplex RT-PCR.

    J Vet Diagn Invest 2001, 13:516-520. PubMed Abstract | Publisher Full Text OpenURL

  20. Elschner M, Prudlo J, Hotzel H, Otto P, Sachse K: Nested reverse transcriptase‒polymerase chain reaction for the detection of group A rotaviruses.

    J Vet Med B Infect Dis Vet Public Health 2002, 49:77-81. PubMed Abstract | Publisher Full Text OpenURL

  21. Higgins DG, Sharp PM: CLUSTAL: a package for performing multiple sequence alignment on a microcomputer.

    Gene 1988, 73:237-244. PubMed Abstract | Publisher Full Text OpenURL

  22. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees.

    Mol Biol Evol 1987, 4:406-425. PubMed Abstract | Publisher Full Text OpenURL

  23. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences.

    J Mol Evol 1980, 16:111-120. PubMed Abstract | Publisher Full Text OpenURL

  24. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.

    Mol Biol Evol 2011, 28:2731-2739. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Lu Z, Song L, Wei-wei X, Yuan-cheng Z, Rui Z, Yu-han C, Jang-jiang W, Ling Z, Zhi-wen X, Wan-zhu G: Development and application of a RT-PCR assay for detection of Porcine Torovirus(in Chinese).

    Chinese Vet Sci 2013, 42:1264-1267. OpenURL

  26. Park S-J, Kim H-K, Moon H-J, Song D-S, Rho S-M, Han J-Y, Nguyen V-G, Park B-K: Molecular detection of porcine kobuviruses in pigs in Korea and their association with diarrhea.

    Arch Virol 2010, 155:1803-1811. PubMed Abstract | Publisher Full Text OpenURL

  27. Oberste MS, Maher K, Kilpatrick DR, Pallansch MA: Molecular evolution of the human enteroviruses: correlation of serotype with VP1 sequence and application to picornavirus classification.

    J Virol 1999, 73:1941-1948. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Oberste MS, Maher K, Kilpatrick DR, Flemister MR, Brown BA, Pallansch MA: Typing of human enteroviruses by partial sequencing of VP1.

    J Clin Microbiol 1999, 37:1288-1293. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Maestre AM, Garzon A, Rodriguez D: Equine torovirus (BEV) induces caspase-mediated apoptosis in infected cells.

    PLoS One 2011, 6:e20972. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL