key: cord-0984411-9sq69psa
authors: Sun, Huan; Lan, Daoliang; Lu, Lifang; Chen, Molin; Wang, Changsong; Hua, Xiuguo
title: Molecular characterization and phylogenetic analysis of the genome of porcine torovirus
date: 2013-10-12
journal: Arch Virol
DOI: 10.1007/s00705-013-1861-x
sha: 74dc448e76b17fd5deff3ea7ac923d5260a6f486
doc_id: 984411
cord_uid: 9sq69psa

In this study, we amplified and sequenced the first genome of porcine torovirus (PToV SH1 strain). The genome was found to be 28,301 bp in length, sharing 79 % identity with Breda virus. It mainly consists of replicase (20,906 bp) and structural genes: spike (4,722 bp), membrane (702 bp), hemagglutinin-esterase (1,284 bp), and nucleocapsid (492 bp). Sequence alignments and structure prediction suggest genetic differences among toroviruses, mainly in NSP1 (papain-like cysteine proteinase domain). Rooted phylogenetic trees were constructed based on the 3C-like proteinase and RNA-dependent RNA polymerase genes. PToV, Berne virus and Breda virus were clustered together, forming a separate branch from white bream virus that was distant from that of the coronaviruses.

Toroviruses, potential intestinal pathogens that can infect horses, calves, humans, pigs and other animals, belong to the subfamily Torovirinae within the family Coronaviridae, order Nidovirales (ICTV, http://ictvonline.org/ virusTaxonomy.asp). The first identified torovirus was isolated from a horse with diarrhea (equine torovirus, EToV, also called Berne virus) in 1972 [23, 27, 29] . A few years later, torovirus was also found in American calves with diarrhea (bovine torovirus, BToV, also called Breda virus) [29] and in humans with gastroenteritis (human torovirus, HToV) [2, 14, 15] . The presence of torovirus in swine fecal specimens was also verified, but there has not been a report of obvious clinical symptoms or economic losses [16] [17] [18] [19] [20] 22] .

Toroviruses are an enveloped viruses with a positivesense single-stranded RNA genome of about 25-30 kb in length, sharing a similar genomic organization with coronavirus [5, 7, 28] . The genome contains two open reading frames connected by a frameshift site to express a replicase polyprotein ORF1ab (PP1ab), and four open reading frames encoding structural proteins: spike (S), membrane (M), hemagglutinin-esterase (HE) and nucleocapsid (N) [11] [12] [13] 30] .

Epidemiological surveys have shown that porcine torovirus (PToV) has been distributed worldwide and is prevalent in Africa, Europe and Asia [16] [17] [18] [19] [20] 22] . However, research on this virus has been limited, which is partly because of insufficient information as well as a lack of an available in vitro system for conducting further studies. Limited genomic information has also restricted our understanding of the virus.

In this study, we completed the full-length sequencing of the PToV genome based on the genome of the Breda1 virus. All information regarding PToV genomic RNA may be further analyzed to explain the different pathogenicities of toroviruses in horses, calves, humans and pigs.

A PToV-positive stool specimen was collected from a commercial piggery in Shanghai. The sample was diluted 1:1 (v/v) with PBS (0.01 M, pH 7.2-7.4), homogenized by vortexing, and clarified by centrifuging at 5,0009g for 5 min at 4°C. The sample was confirmed to contain only one strain of PToV (PToV SH1 strain) by reverse transcription (RT) PCR and sequencing.

RNA was extracted from 300 ll of the fecal specimen supernatant with TRIzol Reagent (Invitrogen, Carlsbad, USA), following the manufacturer's instructions. The isolated RNA was dissolved in 25 ll RNase-free dH 2 O. The RT reaction was performed using a PrimeScriptÒ RT Reagent Kit (Takara, Dalian, China). Briefly, purified viral RNA 16 ll was mixed with 4 ll of PrimeScriptÒ RT Master Mix to a total volume of 20 ll. The reaction mixture was heated to 42°C for 5 min, chilled on ice and then incubated for 1.5 h at 37°C followed by 5 s at 85°C to inactivate the reverse transcriptase enzyme.

Primer 5.0 was used for primer design (synthesized by Sangon, Shanghai, China). Seventeen pairs of primers were designed based on Breda1 virus and the most similar PToV sequences for obtaining the majority of the PToV genome. Full-length genome synthesis was achieved by using six additional pairs of primers from sequenced fragments.

PCR was performed using 2 ll of cDNA and a master mix containing 5 ll 109 LA PCR buffer, 2 ll of each primer (20 lm), 8 ll dNTP mixture (2.5 mM each), 0.5 ll LA Taq (Takara, Dalian, China), and 30.5 ll of ddH 2 O in a total reaction volume of 50 ll. The mixture was denatured at 94°C for 1 min, followed by 35 cycles of 45 s at 94°C, 45 s at 50°C and 2 min at 72°C, and a final extension at 72°C for 10 min. There were minor modifications to the cycling conditions depending on the primers and the length of the products. PCR products were extracted from an agarose gel using an AxyPrep TM DNA Gel Extraction Kit. Purified PCR products were then ligated into a PMD-18T vector (Takara, Dalian, China). For each product, three to five positive colonies were selected and sequenced. Several purified products with much longer nucleotide sequences were sequenced directly using PCR primers.

Polyprotein ORF1ab (pp1ab) was compared with the sequence databases using HMMER [8] and HHpred [3, 10] . MUSCLE [6] was used for producing alignments of nidovirus proteins. Alignments were prepared in publication quality with Jalview [26] . Maximum-likelihood phylogenetic trees were constructed using the WAG model by MUSCLE alignments employing MEGA.v5.05 software [24] . Antigen analysis was performed using Protean in DNAStar [4] . The program KnoTInFrame was used for prediction of ribosomal frameshifts and slippery sequences [25] . The pseudoknot structure was viewed using pknotsRG [21] .

The complete genome of the PToV SH1 strain was composed of 28,301 bp, with a GC content of 35.02 %, while that of Breda virus was 28,475 bp in length and had a GC The genome sequence has been deposited in GenBank and is named porcine torovirus SH1 strain. The accession number is JQ860350.

The Conserved domains of the 3C-like serine protease (3CLpro) and RNA-dependent RNA polymerase (RdRp) were predicted from sequence alignments employing HMMER, Pfam, and MUSCLE and structure prediction using HHpred with PSI-BLAST. PToV 3CLpro was 287 aa in length, the same length with that of Breda virus, and 5 aa shorter than that of Berne virus [17] . The length of RdRp in PToV was 1,019 aa, which was longer than those of Breda virus (1,012 aa) and Berne virus (1,010 aa).

There is a slippery sequence and an RNA pseudoknot at position 14,031, where a ribosome frameshift occurs and forms a fusion protein, ORF1ab. The PToV heptanucleotide slippery sequence, UUUAAAC (nt 14,031-14,037), is identical to those of BToV and EToV [6] .

The PToV HE was 10 aa longer than those of Breda virus and HToV. In antigenicity analysis, these 10 amino acids formed an antigen site and a specific I-E region of an MHC motif at position 6 ( Fig. 2) . In contrast to Breda virus and Berne virus, PToV S contained three discontinuous deletions, with the sizes of 3, 1 and 4 amino acids at position 195, 240 and 679, respectively. There is also a 2-aa insertion at position 743. Whether the amino acid deletions and insertions in PToV are associated with its low virulence needs to be evaluated in future studies [1, 9] .

To infer the direction of torovirus evolution among toroviruses in the family Coronaviridae, phylogenetic trees were constructed using the domains 3CLpro and RdRp, which are conserved in members of the order Nidovirales. Arteriviruses were used as an outgroup to root the trees.

Despite similar topologies in the two trees, there were somewhat incongruent phylogenies (Fig. 3) . For the 3CLpro tree, PToV was clustered with Breda virus, forming a separate branch with Berne virus. For the RdRp tree, the phylogenies were represented more closely within groups. In the torovirus group, PToV formed a branch with Breda virus and Berne virus, separate from white bream virus. However, both trees showed that PToV was more distantly related to Berne virus than to Breda virus.

Although there have been reports on torovirus infection in horses, calves, pigs and humans, evidence of pathogenicity was mainly focused on BToV, EToV and HToV. There is no evidence that PToV causes diarrhea [1, [16] [17] [18] [19] [20] 22] . The main difference in the non-structural genes among toroviruses was found in replicase ORF1a, where bases were lost or replaced frequently. This variation hampered our efforts to amplify the segment. Compared with Breda and Berne virus, PToV SH1 strain pp1ab lacked 28 and 152 amino acids, respectively and had numerous amino acid substitutions, which could possibly result in functional changes in some domains and adaption to the host.

Our primers, specific for non-structural genes, were designed based on Breda1 virus. The genome is highly similar except for the first 26 nucleotides of the 5 0 UTR and the last 22 nucleotides of the 3 0 UTR. The 5 0 UTR and 3 0 UTR were obtained easily and are conserved, giving us confidence that most of the bases in the primers are correct.

The HE protein is located on the surface of the virion. Antibodies against the HE protein can neutralize the virus. The inserted 9-10 amino acids in PToV HE formed an antigen site and a specific I-E region of the MHC motif, which is a characteristic of porcine torovirus only. The spike protein is responsible for the antigenic properties of the virus and contains the binding site for the cell-surface receptor. As the values of the antigen index did not show obvious divergence in comparison with Breda virus and Berne virus, the discontinuous deletions and insertions in the PToV spike protein might be the result of host specificity or low virulence. However, this hypothesis needs to be studied further with more evidence from genomic information and experimental data.

Phylogenetic analysis demonstrated that there is a rather large evolutionary distance between toroviruses and arteriviruses and coronaviruses. Rooted RdRp phylogenic tree analysis revealed that toroviruses were closer to the outgroup node. However, in the 3CLpro tree analysis, the torovirus group was farther away from the outgroup. The topologies of toroviruses were also somewhat different. Nonetheless, the evolutionary distance among Berne virus, PToV and Breda virus was similar. On the other hand, the length of 3Clpro and RdRp proteins varied frequently in toroviruses. The changes were in accordance with the phylogenies.

In conclusion, PToV SH1 genome sequencing not only provides more information on porcine torovirus but also contributes to a better understanding of toroviruses at the molecular level. The genome sequence, combined with the BToV genome and partial EToV genome sequence, is expected to facilitate further research on virus antigenicity, virulence, infection, and replication in animals.

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

An enveloped virus in stools of children and adults with gastroenteritis that resembles the Breda virus of calves

The MPI bioinformatics toolkit for protein sequence analysis

DNASTAR's Lasergene sequence analysis software

The complete sequence of the bovine torovirus genome

MUSCLE: multiple sequence alignment with high accuracy and high throughput

A morphological study of the replication of Breda virus (proposed family Toroviridae) in bovine intestinal cells

HMMER web server: interactive sequence similarity searching

Coronavirus spike proteins in viral entry and pathogenesis

Fast and accurate automatic structure prediction with HHpred

The nucleocapsid of Berne virus

The peplomers of Berne virus

Surface proteins of Breda virus

Association of torovirus with acute and persistent diarrhea in children

Electronmicroscopic evidence of torovirus like particles in children with diarrhoea

Identification and characterization of a porcine torovirus

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

Molecular characterization of a new PToV strain. Evolutionary implications

Longitudinal serological and virological study on porcine torovirus (PToV) in piglets from Spanish farms

Detection of porcine torovirus by real time RT-PCR in piglets from a Spanish farm

Design, implementation and evaluation of a practical pseudoknot folding algorithm based on thermodynamics

Detection and molecular characterization of porcine toroviruses in Korea

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

Genome sequence of porcine torovirus 777

MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods

KnotInFrame: prediction of -1 ribosomal frameshift events

Jalview version 2-a multiple sequence alignment editor and analysis workbench

Purification and partial characterization of a new enveloped RNA virus (Berne virus)

Morphogenesis of Berne virus (proposed family Toroviridae)

Studies with an unclassified virus isolated from diarrheic calves

The haemagglutinating activity of Berne virus