key: cord-0682751-qcpua9ck
authors: Park, Su-Jin; Oh, Eun-Hee; Park, Sang-Ik; Kim, Ha-Hyun; Jeong, Young-Ju; Lim, Guem-Ki; Hyun, Bang-Hun; Cho, Kyoung-Oh
title: Molecular epidemiology of bovine toroviruses circulating in South Korea
date: 2008-01-25
journal: Vet Microbiol
DOI: 10.1016/j.vetmic.2007.07.012
sha: 2d3b6840c1dcd17a5d76835c28a9298416accffe
doc_id: 682751
cord_uid: qcpua9ck

The prevalence of the bovine torovirus (BToV) and its genetic characterization have been reported in North America, Europe and Japan. Therefore, this study examined the prevalence and genetic diversity of the BToV in a total of 645 diarrheic fecal samples from 629 Korean native beef calf herds using RT-PCR and nested PCR with the primer pairs specific to a part of the BToV membrane (M) gene. Overall, 19 (2.9%) out of 645 diarrheic samples from 19 herds (6.9%) tested positive for BToVs by either RT-PCR or nested PCR. A comparison of the nucleotide (nt) and amino acid (aa) sequences of a part of the BToV M gene (409 bp) among the BToVs showed the Korean BToVs to have comparatively higher sequence homology to the Japanese and Dutch BToVs than to the American and Italian BToVs. Generally, the Korean BToV strains clustered with the Japanese and Dutch BToV strains. However, the American and Italian BToV strains clustered on a separate major branch, suggesting that these are more distantly related to other known BToV strains. These results suggest that the BToV infections are sporadic in diarrheic calves in South Korea, and the Korean BToV strains are more closely related to the Japanese and Dutch BToVs than to the American and Italian BToVs.

Toroviruses within the family Coronaviridae are spherical, oval, elongated, or kidney-shaped enveloped viruses that possess a positive-sense single-stranded, polyadenylated RNA genome of approximately 25-30 kb in length (Cornelissen et al., 1997; Horzinek, 1999; Snijder and Horzinek, 1993) . Toroviruses have been detected in normal horses as well as in humans, cattle, pigs, and recently in turkeys with diarrhea (Ali and Reynolds, 2000; Beards et al., 1986; Duckmanton et al., 1997 Duckmanton et al., , 1998a Kroneman et al., 1998; Woode, 1994) . The first report of toroviruses described the identification of an unclassified virus www.elsevier.com/locate/vetmic Veterinary Microbiology 126 (2008) [364] [365] [366] [367] [368] [369] [370] [371] from diarrheic calves (Woode et al., 1982) , which is now known as the Breda virus, or bovine torovirus (BToV) . Fecal shedding of the BToVs have been reported in diarrheic calves in several countries including the United States (Hoet et al., 2002 (Hoet et al., , 2003a Woode et al., 1982 Woode et al., , 1985 , Canada (Duckmanton et al., 1998a) , Costa Rica (Pérez et al., 1998) , Netherlands (Koopmans et al., 1991) , Germany (Liebler et al., 1992) , Hungary (Matiz et al., 2002) , Austria (Haschek et al., 2006) , Japan (Ito et al., 2007) , and South Africa (Vorster and Gerdes, 1993) . However, there are no reports of BToV infections in South Korea. Since the molecular analysis of the BToV has only been carried out in North American, European and Japanese BToV strains (Cornelissen et al., 1997; Draker et al., 2006; Duckmanton et al., 1998b; Gonzálex et al., 2003; Haschek et al., 2006; Hoet et al., 2002 Hoet et al., , 2003a Ito et al., 2007; Smits et al., 2003) , it is unclear if the BToVs circulating in other countries have distinct genetic characteristics. This paper reports the detection of BToV shedding in diarrheic calves using RT-PCR and nested PCR, along with the genetic diversity of the BToV strains based on a partial BToV membrane (M) gene.

According to farmers' request for the treatment of calf diarrhea, local veterinary clinicians visited the cattle farms, and treated diarrheic calves with 2-3 days fasting and the administration of an electrolyte fluid and antibiotics. All the diarrheic calves were aged between 3 and 70 days at the time of sample collection, were fed colostrum from their dams, and housed in the stall. The farm sizes ranged from 5 to 153 animals (average 53). The farms were geographically distributed in Chonnam province, which is the South-Western province in South Korea. In order to diagnose the enteric pathogens, local veterinary clinicians sampled one or two diarrheic fecal specimens from each farm and immediately submitted the sample to the College of Veterinary Medicine, Chonnam National University.

A total of 645 diarrheic fecal specimens were collected from 629 Korean native beef (Hanwoo) calf herds (aged 3-70 days) between January 2004 and December 2005 during the spring (407 samples/ 406 herds), summer (107 samples/98 herds), autumn (73 samples/69 herds) and winter (58 samples/56 herds). Upon arrival of the fecal samples, they were immediately examined for common enteric viral, bacterial and protozoan pathogens including the BToV, bovine coronavirus (BCoV), groups A, B and C bovine rotaviruses (BRV A-C), new bovine enteric calicivirus genus known as becovirus, bovine norovirus (BNoV), bovine viral diarrhea virus (BVDV), Salmonella spp., Clostridium spp., Campylobacter spp., shiga-toxin-producing Escherichia coli, Coccidium spp. and Cryptosporidium spp. (Asakura et al., 1998; Park et al., in press-a,b; Soulsby, 1982; Timoney et al., 1988) . The results of bovine norovirus and bovine coronavirus have been described previously (Park et al., in press-a,b) . Data for other enteric viruses will be reported elsewhere.

The RNA was extracted from a 200 ml starting volume of centrifuged 10% fecal suspensions using the Trizol-LS (Gibco-BRL, Life Tech, Grand Island, USA) procedure. The total RNA recovered was suspended in 50 ml of RNase free water and stored at À80 8C.

The primer sets and RT-PCR and nested PCR conditions used for the detection of the BToV, BCoV, BRV A-C, becovirus, BNoV, and BVDV are described elsewhere (Cho et al., 2001; Park et al., 2006, in pressa,b) . As a negative control, the RNA was extracted from the normal feces of colostrum-deprived calf that had been inoculated with 50 Ml sterile PBS, and used for all RT-PCR assays. In order to confirm the absence of cross-contamination from RT-PCR to nested PCR, RT-PCR was performed with RNA extracted from a normal mock-infected calf fecal sample, and nested PCR was then carried out using 5 ml of this RT-PCR product. The amplification products were analyzed using 1.5 or 2% agarose gel electrophoresis and visualized by irradiating the ethidium bromide stained samples with UV.

The nested PCR products for a portion of the M gene (409 bp) were selected, purified and sequenced directly to verify the reaction specificity as well as to obtain the genomic data for phylogenetic analysis. The nested PCR products were purified using a GenClean II kit (BIO 101, Inc., LaJolla, USA) according to the manufacturer's instructions. DNA sequencing was carried out using an automated DNA sequencer (ABI system 3700, Applied Biosystem Inc., Foster City, USA).

The nucleotide (nt) and deduced amino acid (aa) sequences of the partial BToV M gene were compared with those of the other known toroviruses using the DNA Basic module (DNAsis MAX, Alameda, USA) ( Table 1) . Phylogenetic and bootstrap (1000 replicates) analyses based on the nt and aa alignments were constructed using the neighbor-joining method and the unweighted-pair group method with the average linkages of Molecular Evolutionary Genetics Analysis (MEGA, Version 3.1) with a pairwise distance (Kumar et al., 2004) . A sequence similarity search was carried out for the bovine calicivirus the RdRp protein using the LALIGN Query program of the GENESTREAM network server at the Institut de Génétque Humaine, Montpellier, France (http://www.eng.uiowa.edu/$tscheetz/ sequence-analysis/examples/LALIGN/lalign-guess. html). 

A one-step RT-PCR assay, which targets a 603 bp fragment of the M gene in the BToV, detected three positive fecal samples from three herds. The nested PCR assay, which targeted a 409 bp fragment of the M gene, detected 19 positive fecal samples (2.9%) from 19 herds (6.9%). The calculated percentage of BToV infections in the diarrheic calves in South Korea was comparatively low compared with that reported in other countries, where the reported fecal prevalence of the BToV in calf diarrhea was 5% in Lower Saxony, Germany (Liebler et al., 1992) , 6.5% in Japan (Ito et al., 2007) , 9.7% in the USA (Hoet et al., 2003a) , and 36.4% in southern Ontario, Canada (Duckmanton et al., 1998a) . This suggests that BToV infections are sporadic in diarrheic calves in South Korea.

Of the 19 BToV-positive fecal specimens from the 19 calf herds, 4 fecal samples from 4 herds tested positive to the BToV alone, while the other 15 BToVpositive fecal samples from the 15 herds also tested positive to other enteric pathogens (Table 2 ). This suggests that other bovine enteric pathogens play an important role in the clinical and pathological presentation of this disease. An increase in the severity of the disease after a natural co-infection with other bovine enteric pathogens has been reported (Duckmanton et al., 1998a; Hoet et al., 2003a) .

Seasonally, BToV infections were less prevalent in the diarrheic fecal samples of calves obtained in the spring than in the other seasons (Table 2) ; 4 (1.0%) out of 407 fecal samples tested positive in spring; 1 (1.4%) out of 73 fecal samples tested positive in autumn; 2 (3.4%) out of 58 fecal samples tested positive in winter; and 7 (6.5%) out of 107 fecal samples tested positive in summer. This result is very similar to that reported by Hoet et al. (2003a) , who showed that two peaks of BToV shedding were observed, one during the winter months and the other during the summer season.

Thirteen nested PCR products were selected, purified and sequenced directly, and their nt and aa sequences were compared with those of other known toroviruses (Table 1) . Table 3 shows the sequence relationships of the nt and aa of the partial M gene among the torovirus strains. Interestingly, a comparison of the nt and aa sequences among the BToVs revealed all the Korean BToVs to share 100% nt and aa Table 3 Sequence relationships of the nucleotide and amino acid of the partial membrane gene among torovirus strains identity to each other, except for the K96 strain (99.0% nt and 98.8% aa identity). In order to eliminate cross contamination, the negative control for RT-PCR and nested PCR was used when the field diarrheic fecal samples were analyzed by RT-PCR and nested PCR. However, there was no positive reaction detected in this negative control. Therefore, this identical nt and aa sequence homology of the M gene among the Korean BToVs might be due to the fact that the positive fecal samples used for DNA sequencing originated from the same county, Muan, even though they had been sampled from different farms at different times (Tables 2 and 3 ). The Korean BToVs shared higher nt and aa sequence homology to the Japanese BToVs (98.5-99.5% nt and 98.5-100% aa identity) and the Dutch BToVs (98.3-99.3% nt and 99.3-100% aa identity) than the Italian (96.6% nt and 99.3-100% aa identity) and American (96.6% nt and 98.5-99.3% aa identity) BToVs. In addition, the BToVs including the American, European, Japanese and Korean BToVs had a lower nt and aa identity to the equine torovirus strain, Berne (43.2-44.9% nt and 9.8-10.5% aa identity), than to porcine toroviruses (78.2-81.2% nt and 88.9-91.1% aa identity). From these results, the Korean BToVs were more closely related to the Japanese and Dutch BToVs than to the American and Italian BToVs. In addition, the overall nt and aa sequence homology of the BToV M gene in this study is higher than of the spike (S) protein gene (Ito et al., 2007; Koopmans and Horzinek, 1994) . The relative conservation of the M protein of all BToVs suggests that the structural constraints on this protein are stern, resulting in more limited evolution of this protein than the S protein (Koopmans and Horzinek, Fig. 1 . A phylogenetic tree of the nucleotide sequence of the M gene of the bovine torovirus strains was made using the neighbor-joining method of Molecular Evolutionary Genetics Analysis (Kumar et al., 2004) . The names of the viruses used are listed in Table 1. 1994). Therefore, future studies will be needed to generate more information on the genetic diversity by including other parts of the genome particularly the S gene as well as including the fecal specimens from other regions in South Korea.

Generally, phylogenetic data of nt sequence of the partial BToV M gene revealed that the Korean BToV strains were closely related to the Japanese and Dutch BToV strains (Fig. 1) . Of the Korean BToV strains, all Korean BToV strains clustered together with the exception of the K96 strain, which clustered on a separate branch with the Japanese BToV strains, K-674 and K-676. Other Japanese and Dutch BToVs made separate branches. In addition, the American Breda and the Italian B6 strains clustered on a separate major branch, suggesting these were the most distantly related to the other known BToV strains. Phylogenetic analysis of the aa sequence did not show any significant classification of the BToVs because the homology of the aa sequence among the BToV strains was higher than that of the nt sequence (data not shown). This data supports the above mentioned hypothesis that the Korean BToVs were more closely related to the Japanese and Dutch BToVs than to the American and Italian BToVs. South Korea imports live cows mainly from countries in North American and Oceania, not from Japan and European countries. Hence, the similarity in phylogenetic data on the Korean BToV with the Japanese and Dutch BToVs is quite surprising. Therefore, continuous monitoring and intensive studies will be needed for a more detailed characterization of the BToV with other genes.

In conclusion, BToV infections are sporadic in diarrheic calves in South Korea. Based on the partial sequence of the BToV M gene, the Korean BToVs were more closely related to the Japanese and Dutch BToVs than to the American and Italian BToVs. To our knowledge, this is the first report of the detection of BToV shedding and its genetic diversity in diarrheic calves in South Korea.

Characterization of the stunting syndrome agent: relatedness to known viruses

Detection and genetical characterization of Shiga toxin-producing Escherichia coli from wild deer

Preliminary characterization of torovirus-like particles of humans: comparison with Berne virus of horses and Breda virus of calves

Cross-protection studies between respiratory and calf diarrhea and winter dysentery coronavirus strains in calves and RT-PCR and nested PCR for their detection

Hemagglutinin-esterase, a novel structural protein of torovirus

The complete sequence of the bovine torovirus genome

Detection of bovine torovirus in fecal specimens of calves with diarrhea from Ontario farms

Bovine torovirus: sequencing of the structural genes and expression of the nucleocapsid protein of Breda virus

Characterization of torovirus from human fecal specimens

A comparative sequence analysis to revise the current taxonomy of the family Coronaviridae

Detection of bovine torovirus in neonatal calf diarrhea in lower Austria and Styria (Austria)

Detection of bovine torovirus and other enteric pathogens in feces from diarrhea cases in cattle

Association of enteric shedding of bovine torovirus (Breda virus) and other enteropathogens with diarrhea in veal calves

Enteric and nasal shedding of bovine torovirus (Breda virus) in feedlot cattle

Molecular evolution of corona-and toroviruses

Epidemiological analysis of bovine torovirus in Japan

Toroviruses of animals and humans: a review

Association of diarrhea in cattle with torovirus infections on farms

Identification and characterization of a porcine torovirus

MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment

The significance of bredavirus as a diarrhea agent in calf herds in Lower Saxony

Torovirus detection in faecal specimens of calves and pigs in Hungary: short communication

Molecular epidemiology of bovine noroviruses in South Korea

Detection and molecular characterization of calf diarrhea bovine coronaviruses circulating in South Korea during

Detection and characterization of bovine coronaviruses in fecal specimens of adult cattle with diarrhea during the warmer seasons

Infectious agents associated with diarrhea of calves in the canton of Tilaran Costa Rica. Pre

Phylogenetic and evolutionary relationships among torovirus field variants: evidence for multiple intertypic recombination events

Toroviruses: replication, evolution and comparison with other members of the coronavirus-like superfamily

Primary structure and post-translational processing of the Berne virus peplomer protein

Helminths Arthropods and Protozoa of Domesticated Animals, seventh ed. Lea and Febiger

Hagan and Bruner's Microbiology and Infectious Diseases of Domestic Animals

The toroviruses: bovine (Breda virus) and equine (Berne virus) and the torovirus-like agents of humans and animals

Comparative studies on three isolates of Breda virus of calves

Studies with an unclassified virus isolated from diarrheic calves

This study was supported by grant No. RTI05-01-01 from the Regional Technology Innovation Program of the Ministry of Commerce, Industry and Energy (MOCIE), Republic of Korea. The authors acknowledge a graduate fellowship provided by the Korean Ministry of Education and Human Resources Development through the Brain Korea 21 project.