- Open Access
Molecular epidemiology of Porcine torovirus (PToV) in Sichuan Province, China: 2011–2013
Virology Journalvolume 11, Article number: 106 (2014)
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.
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.
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.
This study uncovered the presence of PToV in Sichuan Province, and demonstrated the need for continuous surveillance PToV of epidemiology.
Porcine torovirus (PToV) is a member of the genus Torovirus, family Coronaviridae, and order Nidovirales. 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 .
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 .
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 , 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 , 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 . 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.
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/). Multiple alignments were achieved using the Clustal_W method of the MegAlign 5.01 program (DNASTAR Inc., Madison, WI, USA) . The phylogenetic tree was constructed using the Neighbor-Joining method on Molecular Evolutionary Genetics Analysis (MEGA) software version 5.0  and carried out using the Kimura 2-parameter model . 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% . 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 . 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).
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.
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.
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.
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 . 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 . In contrast, the lowest prevalence and lowest anti-PToV Ig G levels were found among 3-week-old piglets . Both ELISA reactivity values and rates of seroprevalence rose markedly in piglets from 3 to 11 weeks . 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 . 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 . 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 , however, there is a potential threat that a much more serious outcome will appear in concurrent infections with other enteric pathogens.
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.
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.
Cavanagh D: Nidovirales: a new order comprising Coronaviridae and Arteriviridae. Arch Virol 1997, 142: 629-633.
Ali A, Reynolds D: Characterization of the stunting syndrome agent: relatedness to known viruses. Avian diseases 2000, 44: 45-50. 10.2307/1592506
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. 10.1007/s00705-010-0595-2
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.
Kroneman A, Cornelissen L, Horzinek M, De Groot R, Egberink H: Identification and characterization of a porcine torovirus. J Virol 1998, 72: 3507-3511.
Duckmanton L, Luan B, Devenish J, Tellier R, Petric M: Characterization of torovirus from human fecal specimens. Virology 1997, 239: 158-168. 10.1006/viro.1997.8879
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. 10.1002/jmv.1890200109
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.
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.
Horzinek MC: Molecular evolution of corona- and toroviruses. Adv Exp Med Biol 1999, 473: 61-72. 10.1007/978-1-4615-4143-1_5
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.
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.
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. 10.1007/s00705-013-1861-x
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. 10.1007/s11262-013-0926-y
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. 10.1006/smvy.1997.0104
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. 10.1016/j.vetmic.2010.05.023
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. 10.1186/1756-0500-5-675
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. 10.1016/j.jviromet.2010.06.001
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. 10.1177/104063870101300611
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. 10.1046/j.1439-0450.2002.00510.x
Higgins DG, Sharp PM: CLUSTAL: a package for performing multiple sequence alignment on a microcomputer. Gene 1988, 73: 237-244. 10.1016/0378-1119(88)90330-7
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 1987, 4: 406-425.
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. 10.1007/BF01731581
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. 10.1093/molbev/msr121
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.
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. 10.1007/s00705-010-0774-1
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.
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.
Maestre AM, Garzon A, Rodriguez D: Equine torovirus (BEV) induces caspase-mediated apoptosis in infected cells. PLoS One 2011, 6: e20972. 10.1371/journal.pone.0020972
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).
The authors declare that they have no competing interests.
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.
Lu Zhou, Haoche Wei contributed equally to this work.