|Year : 2019 | Volume
| Issue : 3 | Page : 129-134
Sequence analysis of Turkish field strains of bovine torovirus shows unique amino acid changes in the partial M gene
Hakan Aydin, Mehmet Ozkan Timurkan, Gulizar Acar Kirmizi
Department of Virology, Faculty of Veterinary, University of Ataturk, Erzurum-25240, Turkey
|Date of Submission||24-Dec-2018|
|Date of Decision||20-Jan-2019|
|Date of Acceptance||13-Mar-2019|
|Date of Web Publication||25-Mar-2019|
Atatürk University, Faculty of Veterinary Medicine, Department of Virology, TR-25240 Erzurum
Source of Support: None, Conflict of Interest: None
Objective: To investigate the presence, prevalence and phylogenetic classification of bovine torovirus (BToV). Methods: Stool samples from 72 calves, which were negative for primary gastroenteritis agents, were examined with the nested PCR method by using BToV M gene-specific primers. Results: BToV was detected in 12 (16.7%) out of 72 samples. Phylogenetic analysis was performed using nucleotide and amino acid sequences. In the phylogenetic tree, European, American, Far East and Turkish strains were found to be divided into different branches. Interestingly, it was observed that Turkish strains were divided into two subgroups. Considering the amino acid sequences of these strains having differences at nucleotide level, the change at the 3rd amino acid of the partial M gene in Turkish strains has made Turkish strains different from all other strains. Similarly, the differences were observed in the 18th, 20th, 63rd and 93rd amino acids of the partial M gene only in Turkish field strains. Conclusions: This study revealed that Turkish strains of BToV constitute a separate phylogenetic group and can be divided into two subgroups. In addition, BToV was found to be a common pathogen causing diarrhea in calves in Turkey, and it is a necessity to consider BToV in cases of diarrhea with unknown cause.
Keywords: Calf, Diarrhea, Torovirus, Turkey
|How to cite this article:|
Aydin H, Timurkan MO, Kirmizi GA. Sequence analysis of Turkish field strains of bovine torovirus shows unique amino acid changes in the partial M gene. Asian Pac J Trop Biomed 2019;9:129-34
|How to cite this URL:|
Aydin H, Timurkan MO, Kirmizi GA. Sequence analysis of Turkish field strains of bovine torovirus shows unique amino acid changes in the partial M gene. Asian Pac J Trop Biomed [serial online] 2019 [cited 2019 Jul 15];9:129-34. Available from: http://www.apjtb.org/text.asp?2019/9/3/129/254607
| 1. Introduction|| |
Bovine torovirus (BToV) is an enteropathogen that causes mild to moderate diarrhea, affecting cattle, horses, swines and humans. Toroviruses, a member of the Coronaviridae family in the Nidovirales order, have positive-sense, enveloped, and single-stranded RNA genome,,. Genome of the torovirus is 25-30 kb long, which has a pleomorfic appearance and fringe-like peplomers on its surface, and encodes five different proteins consisting of RNA polymerase, spike (S), membrane (M), hemagglutinin-esterase (HE) and nucleocapsid (N) proteins. Although the torovirus is very similar to the coronavirus in electron microscope visually, there is no antigenic relationship between them,.
Toroviruses have first been diagnosed in a foal with diarrhea in Bern, Switzerland in 1972, and named as Berne virus. Seven years later, in the city of Breda which is in the state of Iowa in the United States, particles similar to Berne virus were detected in the electron microscope in the samples from calves with severe diarrhea, and named as Bredavirus. These viruses, which are prototypes of toroviruses, are now called as toroviruses and are expressed with the species name of the infected animal in taxonomy.
Toroviruses cause severe diarrhea in horses and cattle, whereas they have an asymptomatic characteristic in swines. In humans, it has been reported to cause nosocomial viral gastroenteritis, especially in the pediatric age group.
Among the causes of viral diarrhea, rotavirus and coronavirus are the first agents that come to mind. However, advances in technology and the development of various diagnostic methods have led to the identification of adenovirus, enterovirus, norovirus, and astrovirus as other diarrhea agents,,,. Although BToV has been detected many years ago for the first time, its prevalence and biology in the world have not been fully elucidated yet. Presence of torovirus in calves with diarrhea has been reported in various countries, such as USA,, Austria, Brazil, Germany,, Hungary, Japan, and Korea. However, there is only one study that investigates torovirus in calves in Turkey. Horse breeding is widespread in eastern region of Turkey, due to javelin sport in particular, and horses and cattle share the same barn or environment in this region. For this reason, it is necessary to identify the species of BToV, circulating in Turkey, and reveal these toroviruses with different genetic characteristics.
This pilot study aims to investigate the presence of BToV in eastern Turkey and to reveal its genetic diversity.
| 2. Materials and methods|| |
Stool samples were obtained between 2016 and 2017 from small family farms in eight different centers, namely Erzurum center and Yakutiye, Palandöken, Aziziye, Aşkale, Tortum, Pasinler, and Hasankale [Figure 1]. According to the rapid test kit results, stool samples that have none of the common enteric pathogens (rotavirus, coronavirus, Escherichia coli, Clostridium and Cryptosporidium) were included in the study. Stool samples from a total of 72 calves with diarrhea were used in the study.
|Figure 1: Locations of samples collected from the Province of Erzurum, and surroundings in Turkey.|
White points represent the centers where the samples are collected.
Click here to view
2.2. RNA extraction
Stool specimens that were negative for the initial test (rapid test) were diluted 1/10 and vortexed. Then, they were centrifuged at 3 000 g for 20 min. After centrifugation, 200 μL supernatant was transferred to a new tube and stored at -80 °C until the day of study. The GF-1 Viral Nucleic Acid Extraction Kit (Vivantis, Malaysia) was used for the viral nucleic acid isolation in the stool samples according to the procedure recommended by the manufacturer.
2.3. cDNA synthesis and RT-PCR
Extracted viral nucleic acids were first converted into cDNA. For this purpose, RevertAid™ First Strand cDNA Synthesis Kit (Fermentas Life Sciences, USA) was used. The procedure provided by the manufacturer company was used. The cDNA of the samples was used for detection of BToV by PCR. For the detection of toroviruses, nested PCR was performed using membrane (M) gene-specific primers as defined by the Park et al [Table 1]. After gel electrophoresis was performed, the PCR products of 409 bp size under UV light were considered as positive. Positive samples with strong gel image were subjected to sequence analysis using the ABI Prism BigDye Terminator version 3.1 sequencing kit and Applied Biosystems 310 DNA analyzer (Applied Biosystems Inc., CA, USA) bidirectionally for the BToV M gene using PCR primers.
2.4. Sequencing and phylogenetic analyses
The raw data obtained after the sequencing was confirmed by GenBank (NCBI National Center for Biotechnology Information/ BLAST). Sequences of reference strains and studied samples were aligned with BioEdit (version 7.0.5) program. Phylogenetic and bootstrap (1000 replicates) analysis, based on the M gene nucleotide and amino acid sequence, was performed using the neighbor-joining method under the Molecular Evolutionary Genetics Analysis (MEGA, Version 6.0) program. [Table 2] shows the accession numbers from GenBank for the partial M gene sequences obtained in this study.
| 3. Results|| |
After nested-PCR performed by targeting the partial M gene sequence of 409 bp size of BToV, 12 (16.7%) out of 72 stool samples were found to be positive. Of the 12 positive samples, 9 samples with strong gel images were subjected to sequence analysis and phylogenetic analysis in order to determine their genetic diversity and similarity to other genogroups. Four different lineages were formed after phylogenetic analysis of the positive samples with the reference sequences obtained from GenBank. While all the strains from Turkey formed a group, strains from Italy and the Netherlands formed the European strains group. Similarly, Japanese and Korean strains and American Breda strains were distinctly clustered under different groups. And, the BToV strains identified in Turkey were further divided into two different subgroups (Sub-group 1, Sub-group 2) [Figure 2]. When we look at the two subgroups of Turkish strains, it is seen that the strains collected from the eastern and western regions of Erzurum are distinct from each other, and formed these two subgroups.
|Figure 2: Neighbor-joining phylogenetic trees of 406-nt BToV membrane (M) gene sequences from 9 strains studied and 16 sequences form GenBank.|
Accession numbers and countries of the sequences obtained from GenBank are shown on the tree. BToV strains from this study are shown in bold. The bar graph shows the distance.
Click here to view
When we look at [Figure 3], in which the similarity percentages of nucleotide sequences of Turkish strains and the reference strains are evaluated, we see that our study strains are similar among themselves by 96.4-100.0%. BToV strains of Turkey, which is considered a bridge between Europe and Asia, were found to be similar to the European strains by 97.4-99.0%, as expected. Nucleotide similarity between the studied strains and the Far East (Japan & Korean) strains was 95.1-98.7%. And, the relationship between the strains obtained in our study and the American origin Breda virus (Accession number: Ay427798) was most distant by 94.3-96.1% [Figure 3].
Considering the amino acid sequences of the M genes of BToV, it is noteworthy that Turkish strains clearly have different amino acid sequences compared to Europe, America and Far East strains [Figure 4]. When the amino acid sequences obtained from the Turkish strains of BToV are examined, it is observed that methionine (M), the 3rd amino acid in partial M gene, is replaced by asparagine (N) amino acid. In addition, the amino acid sequences obtained in our study were found to differentiate from other Turkish strains through asparagine (N) to isoleucine (I) transformation in 18th amino acid, tryptophan (W) to glycine (G) in the 20th amino acid, and valine (V) to isoleucine (I) in 63rd and 93rd amino acids. When we look at the American strain “Breda”, we see that the serine amino acid (S) has been exchanged with arginine amino acid (R) change at the 70th amino acid of the partial M gene, constituting a unique sequence.
|Figure 4: Amino acid sequences of the BToV M gene.|
Mutation regions are shown in the box.
Click here to view
| 4. Discussion|| |
There is no sufficient information about the prevalence, biology and genetics of BToV associated with the calf diarrhea. Although there are studies reporting that BToV is causing coronavirus-like winter dysentery in the bovine population, it has not yet been placed in the routine tests,. In our study, calf stool samples which were negative by routine diagnostic tests (cause not determined) were examined for BToV and their contribution in diarrhea etiology was evaluated for our region. To our knowledge, there is only one study on BToV in Turkey, and it has reported that BToV infection is sporadic with 4.7% prevalence in Turkey. In our study which constitutes the second report from Turkey, the BToV prevalence was determined to be 16.7%. This rate is more than three times higher from the previous study, which suggests that the disease is widespread, rather than sporadic, in our regions of study. Furthermore, when we consider that the stool samples used in our study were negative for primary pathogens associated with diarrhea, the prevalence of BToV infection in diarrhea aetiology becomes more prominent. When we look at various studies on the prevalence of BToV in calf diarrhea in the world, we find prevalence rates of 36.4% in the USA, 5.2% in Austria, 43.2% in Crotia, 6.5% in Japan, 6.4% in the Netherlands and 2.9% in South Korea. When compared with the findings from our study, we see that the prevalence of BToV is relatively higher than European countries and Far East (Japan & Korean).
Amino acid and nucleotide sequences of BToV positive amplicons were analysed between themselves and with previously reported sequences. Phylogenetic analysis of the M gene sequences of BToV field strains showed that America, Europe, Turkey and Far Eastern strains are grouped together, and divided into phylogenetic branches. However, it was observed that Turkish strains can be divided into two subgroups. This is thought to be a result of mutations frequently observed in RNA viruses. When we examined the nucleotide sequences of BToV M gene fragments obtained in our study, it was found that they are identical among themselves. Comparison of our study strains with BToV strains previously reported by Gulacti et al. in Turkey showed a sequence similarity rate of 96.1-99.7%. Turkish strains, which form the most distant sequence relationship with American Breda strains by 94.3%, have closer nucleotide sequence similarities with European strains as expected. It is noteworthy that European strains have clustered together in the phylogenetic tree in various studies in Europe,. This shows that BToV strains circulated in Europe are in genetic stability, but there were differences in nucleotide levels observed in Turkish strains.
When we examine the amino acid sequence analysis of the reference strains and our study strains phylogenetically, it is observed that American Breda strain behaves as an out-group probably due to the arginine and serine amino acid exchange in the 70th amino acid. Supporting the hypothesis we have emphasized, Turkey strains are separated from the America, Europe and the Far East strains, especially due to the exchange of the methionine amino acid at the 3rd amino acid with the asparagine. In addition to creating a different cluster, Turkish strains that stand out in our study are divided into sub-groups within themselves. This is because, changes were found in 18th (N→I), 20th (W→G), 63rd (V→I), and 93rd (V→I) amino acids according to all other isolates. These differences identified at the M gene level of BToV cannot be neglected, and it is necessary to clarify this with multi-genomic studies that should be conducted in the future. This has the potential to affect the pathogenicity and virulence of the virus or the diagnosis of the disease.
Mutations are common in RNA viruses. These mutations, which are also common in coronaviruses, can eventually result in virulence-altered strains or different phylogenetic groups and even new species. The virus host spectrum can change and even gain zoonotic properties as a result of mutations,. As an example to mutations from the Coronaviridae family, the feline coronavirus type 2 virus was emerged as a result of recombination between feline coronavirus type 1 and canine coronavirus type 2. The amino acid changes occurred in Turkish BToV field strains in our study are examples of small-scale mutations.
Phylogenetic analysis, nucleotide homology and amino acid sequence analysis all support relationship among different strains. It is noteworthy that clusters are based on the continents or regions. Europe, Far East, Turkey and the America strains formed separate clusters. In addition, strains from western districts of Erzurum (Aziziye and Aşkale) and Samsun strains can be called as sub-group 2, and strains from other districts of Erzurum (Palandöken, Yakutiye, Tortum and Pasinler) and strains from provinces of Adana and Van can be named as sub-group 1. However, this investigation should not only be conducted on a single level of gene region, and further studies should be conducted to support this finding.
As a result, the presence and prevalence of torovirus in north-eastern Anatolia region was investigated for the first time in the study, and positive strains were examined at amino acid and nucleotide level. The only study on BToV in Turkey emphasizes that BToV infection is a sporadic and rare infection in Turkey, however, our rates show that it is far from being sporadic. Future studies should clarify the effect of BToV strains with different nucleotide/amino acid sequences on the clinical and pathological manifestation of the disease. Currently, there is no vaccine for protection against this infection. It was concluded that a vaccine to be developed at the M gene level in the future might be protective both in Europe and in Turkey. In addition, no diarrheal infections were found in horses during the study. Therefore, transition of toroviruses between cattle and horses in our region could not be clarified. It would be beneficial to clarify the cattle-horse and torovirus relationship in future studies.
Although BToV is known to cause an infection in the respiratory and gastrointestinal system in humans and animals, we do not yet have enough knowledge of the transition between species. Epidemiological and genetic studies are a necessity in order to establish a firm basis for the infections. The data obtained in this study provides not only information about the prevalence of BToV, but also important data in the genetic sense.
Conflict of interest statement
We declare no conflict of interest.
| References|| |
Ávila-Pérez G, Rejas MT, Rodríguez D. Ultrastructural characterization of membranous torovirus replication factories. Cellular Microbiol
Atli K, Yildirim Y, Hasircioglu S, Bulut O, Avcı O. The viral agents that cause diarrhea in calves. Ayrinti Derg
Liu D. Toroviruses. In: Liu D (ed). Molecular detection of animal viral pathogens
. USA: CRC Press; 2016, p. 347-352.
Aita T, Kuwabara M, Murayama K, Sasagawa Y, Yabe S, Higuchi R, et al. Characterization of epidemic diarrhea outbreaks associated with bovine torovirus in adult cows. Arch Virol
Haschek B, Klein D, Benetka V, Herrera C, Sommerfeld-Stur I, Vilcek Š, et al. Detection of bovine torovirus in neonatal calf diarrhoea in Lower Austria and Styria (Austria). Zoonoses Public Health
Koopmans M, Horzinek MC. Toroviruses of animals and humans: A review. Adv Virus Res
Gubbay J, Al-Rezqi A, Hawkes M, Williams L, Richardson S, Matlow A. The role of torovirus in nosocomial viral gastroenteritis at a large tertiary pediatric centre. Can J Infect Dis Med Microbiol
Timurkan MO, Aydin H, Aktas O. Frequency and molecular characterization of human norovirus in Erzurum, Turkey. Turk J Med Sci
Alkan F, Karayel I, Catella C, Bodnar L, Lanave G, Banyai K, et al. Identification of a bovine enteric calicivirus, kirklareli virus, distantly related to neboviruses, in calves with enteritis in Turkey. J Clin Microbiol
Gulacti I, Isidan H, Sozdutmaz I. Detection of bovine torovirus in fecal specimens from calves with diarrhea in Turkey. Arch Virol
Gulacti I, Sozdutmaz I, Isıdan H. Molecular characterization of the bovine noroviruses from diarrheic calves in Turkey. Turk J Vet Anim Sci
Duckmanton L, Carman S, Nagy E, Petric M. Detection of bovine torovirus in fecal specimens of calves with diarrhea from Ontario farms. J Clin Microbiol
Hoet AE, Cho KO, Chang KO, Loerch SC, Wittum TE, Saif LJ. Enteric and nasal shedding of bovine torovirus (Breda virus) in feedlot cattle. Am J Vet Res
Nogueira JS, Asano KM, de Souza SP, Brandão PE, Richtzenhain LJ. First detection and molecular diversity of Brazilian bovine torovirus (BToV) strains from young and adult cattle. Res Vet Sci
Liebler E, Klüver S, Pohlenz J, Koopmans M. The significance of bredavirus as a diarrhea agent in calf herds in Lower Saxony. Dtsch Tierarztl Wochenschr
Bidokhti MR, Tråvén M, Ohlson A, Baule C, Hakhverdyan M, Belák S, et al. Tracing the transmission of bovine coronavirus infections in cattle herds based on S gene diversity. Vet J
Matiz K, Kecskemeti S, Kiss I, Adám Z, Tanyi J, Nagy B. Torovirus detection in faecal specimens of calves and pigs in Hungary. Acta Vet Hung
Ito T, Okada N, Okawa M, Fukuyama S, Shimizu M. Detection and characterization of bovine torovirus from the respiratory tract in Japanese cattle. Vet Microbiol
Kuwabara M, Wada K, Maeda Y, Miyazaki A, Tsunemitsu H. First isolation of cytopathogenic bovine torovirus in cell culture from a calf with diarrhea. Clin Vaccine Immunol
Park SJ, Oh EH, Park SI, Kim HH, Jeong YJ, Lim GK, et al. Molecular epidemiology of bovine toroviruses circulating in South Korea. Vet Microbiol
Park SI, Jeong C, Kim HH, Park SH, Park SJ, Hyun BH, et al. Molecular epidemiology of bovine noroviruses in South Korea. Vet Microbiol
Hall TA. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol
Lojkic I, Kresic N, Simic I, Bedekovic T. Detection and molecular characterisation of bovine corona and toroviruses from Croatian cattle. BMC Vet Res
Ito T, Okada N, Fukuyama S. Epidemiological analysis of bovine torovirus in Japan. Virus Res
Dhama K, Pawaiya RVS, Chakraborty S, Tiwari R, Verma AK. Toroviruses affecting animals and humans: A review. Asian J Anim Vet Adv
Aydin H, Timurkan MO. A pilot study on feline astrovirus and feline panleukopenia virus in shelter cats in Erzurum, Turkey. Rev Med Vet
Aydin H, Aktas O. Rotavirus genotypes in children with gastroenteritis in Erzurum: first detection of G12P 
and G12P 
genotypes in Turkey. Prz Gastroenterol
Smits S, Lavazza A, Matiz K, Horzinek M, Koopmans M, De Groot R. Phylogenetic and evolutionary relationships among torovirus field variants: Evidence for multiple intertypic recombination events. J Virol
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2]