key: cord-0941489-0vam15pt
authors: Li, Hao; Zhang, Bin; Yue, Hua; Tang, Cheng
title: First detection and genomic characteristics of bovine torovirus in dairy calves in China
date: 2020-05-09
journal: Arch Virol
DOI: 10.1007/s00705-020-04657-9
sha: 75beab4b2dac72a75c05fffcd113df9b77edbf18
doc_id: 941489
cord_uid: 0vam15pt

Bovine torovirus (BToV) is a diarrhea-causing pathogen. In this study, 92 diarrheic fecal samples from five farms in four provinces in China were collected and tested for BToV using a RT-PCR assay, and 21.73% samples were found to be BToV positive. Moreover, two complete BToV genome sequences (MN073058 and MN073059) were obtained from the clinical samples, which were 28,297 and 28,301 nucleotides in length, respectively. Sequence analysis showed that the two isolates shared 10 identical amino acid mutations in the S protein compared to the complete S sequences of BToV available in the GenBank database. In addition, seven consecutive amino acid mutations were found from aa 1,486 to 1,492 in the S protein of isolate MN073058. Moreover, the two isolates shared one identical amino acid mutation in the receptor binding sites of the HE protein. To the best of our knowledge, this is the first report on the epidemic and genomic characterization of BToV in China, which is helpful for further understanding the genetic evolution of BToV. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00705-020-04657-9) contains supplementary material, which is available to authorized users.

Toroviruses are members of the family Coronaviridae, order Nidovirales, and include bovine torovirus (BToV) [1, 2] , Berne virus (EToV) [3] , porcine torovirus (PToV) [4] , and human torovirus (HToV) [5] . BToV mainly causes diarrhea in calves and adult cattle. The virus can not only be detected in feces but also in the respiratory tract, indicating that the virus has dual tissue tropism [6] [7] [8] [9] . BToV has been detected in 16 countries with a wide geographical distribution [2, 7, [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] . In addition, a sequence reported to be from goat torovirus is present in the GenBank database, and in 2012, researchers detected the presence of PToV in pig herds by RT-PCR in China [20] .

Presently, three BToV genomic sequence (AY427798, LC088094, LC088095) can be obtained from the GenBank database. The BToV genome encodes four structural proteins: spike glycoprotein (S), membrane glycoprotein (M), hemagglutinin esterase (HE), and nucleocapsid glycoprotein (N) [21] . The S protein is involved in viral infectivity and induces the production of neutralizing antibodies [22] . The M protein plays a role in BToV assembly and nucleocapsid recognition [21] . The HE has a putative F-G-D-S motif and has acetylesterase activity specific for N-acetyl-9-Oacetylneuraminic acid [1, 21] . HE contains three functional domains: a lectin domain (R), an esterase domain (E) and a membrane-proximal domain (MP) [23] . The N protein is the only viral RNA-binding polypeptide found in infected cells. It protects the genome and ensures its timely replication and reliable transmission, as well as playing a role in virus transcription and translation [24, 25] .

Diarrhea is a common disease in dairy cows in China, which leads to serious economic losses. Bovine coronavirus (BCoV), bovine group A rotavirus (BRVA), bovine viral diarrhea virus (BVDV), bovine norovirus (BNoV), and nebovirus (NeV) had been identified as common diarrheacausing viruses in bovines in China [26] [27] [28] [29] [30] . However, there is currently no information regarding BToV in China. The goal of this study was to detect and the genome of characterize of BToV in calves in China.

From September to December 2018, 92 diarrheic fecal samples were collected from calves (aged < 3 months) from five farms in four provinces of China. These included 23 samples from Sichuan Province (one farm), 39 from Liaoning Province (two farms), 20 from Shanxi Province (one farm), and 10 from Henan Province (one farm). All samples were stored at -80 °C in sterile 50-mL centrifuge tubes until further testing.

Total RNA from 300 μL of the fecal suspension was extracted using a QIAamp Viral RNA Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer's instructions. Complementary DNA (cDNA) was synthesized using a Primescript™ reverse transcription kit (TaKaRa, Dalian, China) and stored at − 20 °C.

BToVs were detected using a specific RT-PCR assay established in our laboratory. The specificity and reproducibility of the RT-PCR assay has been validated, and the detection limit is 1.908 × 10 −1 pg/µL. Briefly, a pair of primers 

To investigate coinfection with BRVA, BCoV, BVDV, NeV and BNoV, all BToV-positive samples were subjected to previously described specific RT-PCR assays for these viruses. [27] [28] [29] 31] The detection of BRVA, BCoV, NeV were following our previous reports, BVDV was following the previous report [], BNoV was following the previous report [].

To investigate the genomic characteristics of BToV isolates from China, 29 primer pairs were designed to amplify their genomes (Table S1 ). The 3' and 5' ends of the viral genome were sequenced by rapid amplification of cDNA ends using a Smart RACE cDNA Amplification Kit (TaKaRa). The PCR products were cloned into the pMD19-T simple vector (TaKaRa) for further sequencing (Sangon Biotech).

Sequences were assembled using SeqMan software (version 7.0; DNASTAR, Inc., WI, USA). Putative ORFs in the linear genomes were identified using the ORF Finder tool (http://www.ncbi.nlm.nih.gov/gorf/gorf.html). Nucleotide and deduced amino acid sequences were compared using the MegAlign program of Lasergene software, version 7.1 (DNASTAR, Madison, WI, USA). A multiple sequence alignment was performed, and a neighbor-joining phylogenetic tree was built using the interior branch test method with the aid of MEGA7 software. Possible recombination events were identified using SimPlot software (version 3.5.1) and the Recombination Detection Program (RDP 4.95) by the RDP, Chimaera, BootScan, 3Seq, GENECONV, Max-Chi, and SiScan methods.

Among the 92 diarrheic samples, 21.73% (20/92) tested BToV-positive by RT-PCR. The detection rate was 69.56% 

A portion of the M gene of 20 positive samples was sequenced, and all partial sequences of the M gene (Gen-Bank accession no. MN073058-MN073079) were submitted to the GenBank database. The 20 nucleotide sequences were 98.8% to 100% identical to one another and 96.3% to 99.3% identical to those of previously reported BToV strains. A phylogenetic tree was constructed using the 406bp M fragment sequence, including the 20 BToV sequences from calves as well as 10 other BToV sequences from the GenBank database. The 20 BToV sequences from this study were divided into two different groups (Fig. S1 ). One of these groups included BToV isolates from three different provinces, which formed a unique large branch, but the one remaining strain clustered with the Turkish strains HT1 and HT2.

To obtain more-precise information about the evolutionary relationships of the BToV isolates from calves, six representative BToV isolates were chosen from four positive farms, and two BToV isolates were chosen from each province. However, we were only able to successfully assemble two complete genomic sequences (SC-1 Sichuan / 2018 and SC-2 Sichuan / 2018), both of which were from the same farm in Sichuan province. These sequences were submitted to the GenBank database with accession numbers MN073058 and MN073059. The length of the linear genome was 28,297 and 28,301 nt, respectively, and the G + C content was 36.36%. The complete genome of the strains shared 42.9-99.2% nt sequence identity with the 10 genome sequences of mammalian ToVs in the GenBank database. Analysis of the complete genome sequences revealed that the two Chinese BToVs shared 99.1% nucleotide sequence identity with each other, 96.89% to 97.07% identity with the two genome sequences from Japan, and shared 81.9% nt identities with the prototype strain Bread 1 (Fig. 1 ).

The complete S genes of strains SC-1 and SC-2 were both 4,755 bp in length and encoded a protein of 1,586 amino acids. The complete S genes shared 99.7% nt sequence identity and 99.3% aa sequence identity with each other and shared 92.2-97.6% nt sequence identity (91.3-98.4% aa sequence identity) with the 12 BToV S sequences available in the GenBank database. A phylogenetic analysis based on the complete amino acid sequence of the S protein showed that the BToVs could be separated into four groups (Fig. 2) , designated tentatively as group 1 to group 4. The two BToVs clustered with the BToV strains B145 and HT2 belonging to group 3 and were independently clustered in a small branch. Further analysis showed that, compared to the 12 available BToV S sequences, strains SC-1 and SC-2 had 10 amino acid mutations, of which two were located in the signal peptide, six were located in the putative S1 region, and two were located in the putative S2 region. In addition, SC-1 had a Fig. 1 Phylogenetic tree based on the complete genomic nucleotide sequence of all mammalian toroviruses. Sequence alignments were performed using ClustalW in MEGA 7.0 software. The tree was constructed by the maximum-likelihood method with bootstrap values calculated for 1000 replicates. The scale bar indicates amino acid substitutions per site. Bovine torovirus was found in calves and may be related to members of the genus Coronavirus. The bovine torovirus strains BToV/SC-1/China and BToV /SC-2/China investigated in this study are indicated by black triangles Fig. 2 Phylogenetic tree based on the deduced 1586-aa sequence of the complete S gene. Sequence alignments were performed using ClustalW in MEGA 7.0 software. The tree was constructed by the maximum-likelihood method with bootstrap values calculated for 1,000 replicates. The bovine torovirus strains investigated in this study are indicated by by black triangles continuous 7-amino-acid mutation at the end of the putative S2 region, and SC-2 had a unique amino acid mutation in the putative S2 region (Fig. 3 ).

The complete HE gene of strains SC-1 and SC-2 are 1,260 bp in length and encode a protein of 419 amino acids. F-G-D-S, the putative esterase active site in all HE proteins, was located at aa positions 34-37. This complete HE gene shares 77.1-98.7% nt sequence identity (76.2-98.1% aa sequence identity) with all 16 available BToV HE gene sequences. Phylogenetic trees constructed using full-length HE gene sequences showed that the two HE sequences clustered into the same subgroup and belonged to genotype II (Fig. 4) .

Furthermore, there was a unique amino acid variation (L202F) in strains SC-1 and SC-2 compared to the 16 available BToV HE sequences, and it is interesting that this mutation is located in the neutralizing epitope. In addition, the SC-1 strain had a unique amino acid variation (I289Y), and the SC-2 strain had three unique amino acid variations (K137S, S238Q, W290G) compared to the other available BToV HE sequences (Fig. 5 ).

BToV SC-1and SC-2 were predicted to be recombination strains using the Recombination Detection Program (RDP 4.0) software, and standard similarity plot analysis was performed using Simplot 3.5.1 (Fig. 6) . Notably, both strains SC-1 and SC-2 had the same predicted breakpoints in the crossover region, at nucleotide positions 20,810 and 28,047 in the 3' end of ORF1b and the 5' end of the N coding region, respectively.

In this study, 21.73% (20/92) of the diarrheic stool samples collected from dairy cows were found to be BToV positive. These samples were from four farms in three provinces, and the two more distant farms were separated by more than 1500 kilometers. The results indicate that BToV has circulated in Chinese cows with wide geographical (Table 1) , which is similar to the previous report from the USA [32] . This may lead to an increase in clinical severity and increased difficulty in diagnosing and controlling calf diarrhea.

In this research, we determined the obtained two complete genome sequences of two BToV isolates from the same farm in Sichuan province, increasing the number of BToV genome sequences in the GenBank database to five, thus contributing to a better understanding of the genome structure and genetic evolution of BToV. Phylogenetic analysis indicated that these two BToV isolates had a close genetic relationship to strains from Japan. Strains SC-1 and SC-2 were predicted to be recombinants derived from BToV strains Breda1 and PToV, in agreement with a previous report [8] . Moreover, the two Chinese strains shared identical unique amino acid changes in the S and HE genes when compared to the other strains with sequences available in the GenBank database, indicating the unique evolution in Chinese BToV strains. Interestingly, although these two isolates came from the same farm, their S and HE genes had their own unique amino acid variations. These results will be helpful for understanding the genetic evolution of BToV.

The S protein is involved in the induction of neutralizing antibody production during the infection process and plays an important role in the pathogenesis of BToV, which is similar to that of coronaviruses [22, 24, 33] . Another member of the family Coronaviridae, BCoV, has a genomic structure similar to that of BToV, and mutations in the S protein are associated with viral antigenicity, pathogenicity, tissue tropism, and host range change [34, 35] . It has been reported that cleavage most likely occurs between amino acids 1,002 and 1,003 of the mature S protein in BToV, resulting in the formation of the S1 and S2 proteins [24] . The two Chinese strains shared two identical amino acid mutations in the signal peptide that might affect its function. It is worth noting that these two strains have six identical amino acid mutations in the putative S1 region. Since the S1 protein is responsible for receptor recognition and neutralizing antibody induction, mutations at amino acids 60 and 96 (T60I, L96P) might affect virus replication [36] . Therefore, the effects of these unique amino acid variations on the function of the proteins of strains SC-1 and SC-2 S is worth further study.

Based on the HE gene sequences, BTOV can be divided into three genotypes [18] . In this study, a phylogenetic tree constructed using all 18 complete HE gene sequences available in the GenBank database showed that the two HE Fig. 6 Similarity plot analysis of the complete genome sequences of PToV SHI (blue line), BToV Breda (red line), and PToV NPL (pink line), with BToV SC-1 and SC-2 as query sequences, using a sliding window of 500 nt and a moving step size of 20 nt sequences belonged to a genotype that is prevalent throughout the world [10, 12, 25, 37] . The torovirus HE protein contains three domains: MP, E and R, among which the R domain is involved in receptor recognition and plays an important role in the process of BToV infection [38] . The two Chinese strains shared an identical unique amino acid change in the E region. In addition, strain SC-2 had another unique amino acid change in the E region. Therefore, the ability of these mutations to affect viral receptor binding requires further investigation.

In conclusion, this study is the first to confirm the existence and prevalence of BTOV in Chinese calves with diarrhea, contributing to the diagnosis and control of calf diarrhea in China. Moreover, two complete BToV genome sequences were obtained from the clinical samples, and these two BToV isolates had unique amino acid changes in the S and HE proteins. These findings will enhance our understanding of the epidemic and the genetic evolution of BToV.

Nidovirales: a new order comprising Coronaviridae and Arteriviridae

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

Antibodies to Berne virus in horses and other animals

Identification and characterization of a porcine torovirus

Du Pasquier P (1984) An enveloped virus in stools of children and adults with gastroenteritis that resembles the Breda virus of calves

The complete sequence of the bovine torovirus genome

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

Detection and characterization of bovine torovirus from the respiratory tract in Japanese cattle

Studies with an unclassified virus isolated from diarrheic calves

Detection of bovine torovirus in fecal specimens from calves with diarrhea in Turkey

Detection of bovine torovirus in neonatal calf diarrhoea in Lower Austria and Styria (Austria)

Detection of bovine torovirus in fecal specimens of calves with diarrhea in Japan

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

Detection and molecular characterisation of bovine corona and toroviruses from Croatian cattle

First detection and molecular diversity of Brazilian bovine torovirus (BToV) strains from young and adult cattle

Molecular epidemiology of bovine toroviruses circulating in South Korea

Infectious agents associated with diarrhoea of calves in the canton of Tilarán

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

Breda virus-like particles in calves in South Africa

Molecular characterization and phylogenetic analysis of the genome of porcine torovirus

Bovine torovirus (Breda virus) revisited

Structure, function, and evolution of coronavirus spike proteins

Structural basis for ligand and substrate recognition by torovirus hemagglutinin esterases

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

An interaction between the nucleocapsid protein and a component of the replicase-transcriptase complex is crucial for the infectivity of coronavirus genomic RNA

Critical role of cellular cholesterol in bovine rotavirus infection

Molecular investigation of bovine viral diarrhea virus infection in yaks (Bos gruniens) from Qinghai, China

Detection and molecular characteristics of neboviruses in dairy cows in China

Prevalence of a novel bovine coronavirus strain with a recombinant hemagglutinin/esterase gene in dairy calves in China

Prevalence and complete genome of bovine norovirus with novel VP1 genotype in calves in China

Reverse transcription-PCR assays for detection of bovine enteric caliciviruses (BEC) and analysis of the genetic relationships among BEC and human caliciviruses

Case-control study of microbiological etiology associated with calf diarrhea

A single amino acid change within antigenic domain II of the spike protein of bovine coronavirus confers resistance to virus neutralization

Coronavirus spike proteins in viral entry and pathogenesis

Crystal structure of bovine coronavirus spike protein lectin domain

Genetic and antigenic characterization of newly isolated bovine toroviruses from Japanese cattle

Characterization of epidemic diarrhea outbreaks associated with bovine torovirus in adult cows

Hemagglutinin-esterase, a novel structural protein of torovirus

Acknowledgements This study was supported by funding from the