key: cord-0837749-pos7xvu5 authors: Kim, Soo Joung; Han, Jeong Hee; Kwon, Hyuk Moo title: Partial sequence of the spike glycoprotein gene of transmissible gastroenteritis viruses isolated in Korea date: 2003-07-17 journal: Vet Microbiol DOI: 10.1016/s0378-1135(03)00101-9 sha: 413227edb89564cc432e35e54b4455e72f44e8d1 doc_id: 837749 cord_uid: pos7xvu5 The spike (S) glycoprotein of transmissible gastroenteritis virus (TGEV) is the predominant inducer of neutralizing antibodies and has been implicated in virulence and host cell tropism. In this study, the nucleotide and deduced amino acid sequences of the amino terminal half of the S glycoprotein gene of one Korean field TGEV strain (133) isolated in 1997 and three Korean field TGEV strains (KT2, KT3 and KT4) isolated in 2000 and HKT2 strain, KT2 passaged 104 times in ST cells, were determined. The amino terminal half of the S glycoprotein gene including antigenic sites A, B, C and D, were amplified by reverse transcriptase-polymerase chain reaction (RT-PCR). Amplified PCR products were cloned, sequenced, and compared with published sequences for non-Korean TGEV strains. Korea TGEV field strains had 98.5–99.5% nucleotide sequence and 97.2–99.0% amino acid sequence similarity with each other. They had 96.5–99.0% nucleotide sequence similarity and 94.9–97.6% amino acid sequence similarity compared to non-Korean TGEV strains. Korean TGEV strains had several specific nucleotide and amino acid sequences which were not found in foreign TGEV or PRCV strains. HKT2 strain differed by 0.89% in nucleotide and 2.03% amino acid sequences compared to original KT2 strain although the regions forming four antigenic sites were not changed. By phylogenetic tree analysis, Korean field TGEV strains were branched into different groups from non-Korean TGEV or PRCV strains. Korean TGEV field strains KT2 and 133 were branched in separate groups that were differentiated from the other Korean TGEV strains. The Korean TGEV strains seemed to be evolved from a separate lineage of TGEV strain. Transmissible gastroenteritis virus (TGEV) is the etiological agent of transmissible gastroenteritis (TGE), which is a highly contagious enteric disease in piglets causing vomiting, severe yellowish diarrhea, weight loss, dehydration and high mortality and resulting in severe economical losses to the affected farms (Saif and Wesley, 1999) . TGEV belongs to the family Coronaviridae. It is a pleomorphic enveloped virus that contains a large positive-sense, single-stranded RNA genome (Siddell et al., 1983) . TGEV has three major structural proteins which are spike (S) glycoprotein, membrane (M) glycoprotein and nucleocapsid (N) protein. The S glycoprotein forms the peplomers on the virion envelope, induces neutralizing antibodies, contains receptor binding regions, and contains four major antigenic sites (A, B, C and D) Delmas et al., 1986; Gebauer et al., 1991) . The sites A and D are known to be regions inducing major neutralizing antibodies Gebauer et al., 1991) . The M glycoprotein is mostly embedded in the lipid virion envelope and the N phosphoprotein is associated with the genomic RNA to form the nucleocapsid (Kapke and Brian, 1986; Laude et al., 1987) . Porcine respiratory coronavirus (PRCV) was first reported in Europe (Pensaert et al., 1986; Callebaut et al., 1988) and has been reported in US (Wesley et al., 1990b; Vaughn et al., 1994) . By sequence data, PRCV has a large deletion in the 5 -region of the S glycoprotein gene which may be related to the differences observed in tissue tropism or pathogenicity between PRCV and TGEV strains (Rasschaert et al., 1990; Wesley et al., 1990c; Britton et al., 1991) . The sequence variation in the S glycoprotein gene and antigenic diversity of TGEV isolates have been reported (Hohdatsu et al., 1987; Sanchez et al., 1992; Vaughn and Paul, 1993; Paton and Lowings, 1997) . Based on partial sequence analysis of S glycoprotein gene, TGEV strains isolated in Korea were different from foreign TGEV isolates (Kwon et al., 1998) . The purpose of the present study was to genetically characterize TGEV strains (KT2, KT3, KT4 and 133 strains) isolated in Korea and compare them to a highly passaged TGEV KT2 strain (HKT2) as well as other non-Korean published TGEV and PRCV sequences. The 5 -end of the S glycoprotein genes including antigenic sites A, B, C and D of the Korean TGEV isolates were amplified by RT-PCR, cloned, and sequenced. Comparisons with the sequences of other published TGEVs showed that the nucleotide and amino acid sequences of Korean TGEV isolates differed by 1.0-3.5 and 2.4-5.1%, respectively, from non-Korean TGEV strains in the amino terminal half of the S glycoprotein gene. The Korean TGEV field strain 133 and three other TGEV field strains (KT2, KT3 and KT4) were isolated from small intestines of TGE suspected piglets using swine testicular (ST) cells in 1997 and 2000, respectively. The ST cells were grown in Dulbecco's modified Eagle medium (DMEM, GIBCOBRL, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (FBS, GIBCOBRL) and maintained in maintenance medium (MM) which was DMEM supplemented with 2% FBS at 37 • C in a 5% CO 2 incubator. The Korean TGEV KT2 field strain was passaged serially in ST cells. The KT2 strain (HKT2) passaged 104 times on ST cells was analyzed in this study. Virus-infected ST cells showing cytopathic effects were frozen, thawed and cell debris pelleted by centrifugation. The clarified cell culture supernatant was collected and used for preparation of viral RNA. Viral RNAs were extracted and purified as described previously (Kwon et al., 1993) . Briefly, sodium dodecyl sulfate (final concentration 2% w/v) and proteinase K (final concentration 250 g/ml) were added to clarified cell culture supernatant, incubated for 5 min at 55 • C, extracted with acid phenol (pH 4.0) and chloroform/isoamyl alcohol (49:1), and further purified using the RNaid kit (BIO101, La Jolla, CA, USA). Finally, the RNA was resuspended in diethyl-pyrocarbonate (DEP) treated water and stored at −70 • C until used in the reverse transcriptase reaction. The amino-terminal half of the S gene, including antigenic sites A, B, C and D of TGEV, was amplified with forward 2 (5 -AAGGAAGGGTAAGTTGCTCA-3 , −71 to −59)/reverse 3 (5 -GGTCCATCAGTTACGCCGAA-3 , 1159-1184) primer pair (1255 bp PCR product) and forward 4 (5 -AAGGAAGGGTAAGTTGCTCA-3 , 1064-1183)/reverse 6 (5 -GGTCCATCAGTTCGCCGAA-3 , 2344-2363) primer pair (1300 bp PCR product), which were partially overlapped . The RT reaction to synthesize cDNA contained 3 g random primers (cat. no. 48190-011, GIBCOBRL), 10 mM dNTPs (Promega, Madison, WI, USA) and purified RNA. The mixture was heated at 65 • C for 5 min. The 5× first-strand buffer (250 mM Tris-HCl (pH 8.3 at room temperature), 375 mM KCl, 15 mM MgCl 2 ), 0.2 M DTT (GIBCOBRL) and 20-40 units RNase inhibitor (Promega) were added to the mixture and incubated at 25 • C for 10 min and 42 • C for 2 min. Two hundred units Superscript II reverse transcriptase (GIBCOBRL) was added, heated at 42 • C for 50 min then at 70 • C for 15 min. One unit ribonuclease H (GIBCOBRL) was added, heated at 37 • C for 20 min and 80 • C for 10 min. The PCR reaction mixture to amplify cDNA contained 25 l PCR premix (1.25 unit/25 ul TaKaRa Ex Taq polymerase, 0.4 mM dNTPs, Ex TaqTM buffer including 4 mM Mg 2+ ) (TaKaRa, Shiga, Japan), forward and reverse primers and cDNA. A 50 l total reaction volume was obtained by adding distilled water. The PCR was performed by 35 cycles of denaturation at 94 • C for 1 min, annealing at 51 • C for 90 s, and polymerization at 74 • C for 90 s. The initial denaturation and polymerization steps were at 94 • C for 5 min and 74 • C for 6 min, respectively. The final annealing and polymerization steps were at 51 • C for 2 min and 74 • C for 10 min, respectively. The PCR products amplified with each primer pair were purified using GENECLEAN Turbo kit (BIO101, La Jolla, CA, USA) according to the manufacturer's recommendation. The purified DNA was ligated into the pCR2.1-TOPO (Invitrogen, Calsbad, CA, USA) cloning vector and transformed into competent cells (TOP10) (Invitrogen). Cells carrying recombinant plasmid were selected on LB agar plates containing kanamycin and X-gal. Plasmid DNA for sequencing was prepared by Plasmid Maxi kit (Qiagen, Santa Clarita, CA, USA). The sequences were determined using the ALFexpress automatic DNA sequencer (Pharmacia Biotech, Uppsala, Sweden). For each TGEV strain, two or three independent clones originated from different PCR products were sequenced because of the possibility of errors arising during RT, PCR, or the cloning procedure. Nucleotide and predicted amino acid sequences were analyzed with the MegAlign software (DNAStar Inc., Madison, WI, USA). Phylogenetic tree for partial S glycoprotein was generated by the maximum parsimony method with 100 bootstrap replicates in a heuristic search with the PAUP 4.0b software program (Sinauer Associates Inc., Sunderland, MA, USA). The tree was rooted to an sequence of the TGEV KT2 isolate. The sequence data were deposited in the GenBank database under the accession numbers AF481360 (KT2), AF481361 (KT3), AF481362 (KT4), AF481365 (133) and AF481366 (HKT2). The origin and history of sequences used for comparison in this study were described elsewhere (Britton and Page, 1990; Chen et al., 1995; Sanchez et al., 1992; Wesley, 1990a) . Sequences used for comparison in this study were from the following GenBank database accession numbers: TFI83 (Z35758) (Chen et al., 1995) , PUR46-MAD (M94101), NEB72 (M94099), TOY56 (M94103), HOL87 (M94097) (Sanchez et al., 1992) , BRI70-FS (X53128) (Britton and Page, 1990) , MIL65-AME (S51223) (Wesley, 1990a) . The nucleotide and deduced amino acid sequences of the S glycoprotein gene including four antigenic sites (A, B, C and D) of four TGEV strains (133, KT2, KT3 and KT4) isolated in Korea and KT2 strain (HKT2) passaged 104 times on ST cells were determined and compared with the sequences of published TGEV and PRCV strains (Fig. 1) . Korean TGEV isolates had nucleotide sequence similarity between 98.5% (KT2 and 133) and 99.5% (KT3 and KT4) with each other and they had nucleotide sequence similarity between 96.5% (KT2 and TFI83) and 99.0% (133 and NEB7 and PUR46-MAD) with non-Korean TGEV strains (Table 1 ). Korean TGEV isolates had amino acid sequence similarity between 97.2% (KT2 and 133) and 99.0% (KT3 and KT4) with each other and they had amino acid sequence similarity between 94.9% (KT2 and TFI83) and 97.6% (133 and NEB72) with foreign TGEV strains (Table 1) . Korean TGEV strains had six and five specific nucleotide and amino acid sequence differences, respectively, compared to foreign TGEV or PRCV strains. The seven Korean Fig. 1 . The deduced amino acid sequence comparisons of the amino terminal half region of S glycoprotein gene in Korean TGEV isolates KT2, KT3, KT4 and HKT2 strains with Korea field TGEV 133 strain isolated in 1997 (Kwon et al., 1998) , TFI83 (Chen et al., 1995) , MIL65-AME and BRI70-FS (Britton and Page, 1990; Wesley, 1990a) , PUR46-MAD, NEB72, TOY56 and PRCV isolate HOL87 strains (Sanchez et al., 1992) . The shaded regions indicate four antigenic site amino acids. Symbol (-) indicates positions where the sequences are identical to those of KT2. Deletions within the sequences are shown with asterisks. The five Korean TGEV-specific amino acids changes were underlined. TGEV-specific nucleotide changes were 1212 (T to C), 1234 (T or C to A), 1241 (C to T), 1321 (T to C), 1689 (C to A) and 1798 (G to C), and five amino acid changes were 412 (L to I), 414 (T to I), 441 (F to L), 563 (N to K) and 600 (D to H) (Fig. 1) . Korean TGEV KT2 strain had nine and seven specific nucleotide and amino acid sequence changes, respectively, compared to other Korean TGEV isolates and non-Korean TGEV or PRCV isolates. The nine Korean TGEV KT2-specific nucleotide changes were 926 (T), 1369 (T), Phylogenetic analysis of the S glycoprotein genes of Korean TGEV isolates and non-Korean TGEV or PRCV strain is shown in Fig. 2 . Four distinct branches were observed in the phylogenetic trees. Korean TGEV KT3, KT4 and HKT2 strains were branched into one group. Non-Korean TGEV or PRCV strains (TFI83, HOL87, BRI70-FS, MIL65-AME, TOY56, PUR46-MAD and NEB72) were branched into another group in which NEB72 and PUR46-MAD strains formed a subgroup different from the other foreign strains. Korean TGEV KT2 and 133 strains was separated distinctly from the other two groups. The amino-terminal half of the S glycoprotein gene of recent Korean TGEV isolates were sequenced to determine the genetic diversity among Korean TGEV isolates and differences with foreign TGEV isolates. The three TGEV strains (KT2, KT3 and KT4) isolated in 2000 had high similarity, showing 99.2-99.5% nucleotide sequence similarity and 98.4-99.0% amino acid sequence similarity among themselves. But Korean TGEV field strain 133 isolate had low similarity compared to TGEV strains isolated in 2000. TGEV strain 133 showed 98.5% (KT2) to 98.8% (KT3) nucleotide sequence similarity compared to other Korean TGEV isolates and 97.0% (TFI83) to 99.0% (NEB72 and PUR46-MAD) nucleotide sequence homology compared to non-Korean TGEV strains. In amino acid sequences, TGEV 133 strain 97.2% (KT2) to 97.8% (KT3) similarity compared to other TGEV isolates and 95.4% (TFI83) to 97.6% (NEB72) similarity. These data showed that TGEV strains differentiated from non-Korean TGEV strains in genetic composition have existed in Korea. Korean TGEV strains have several unique characteristics. All Korean TGEV strains including the HKT2 strain had six unique nucleotides between positions 1212 and 1798 and five unique amino acid sequences between positions 412 and 600, which were not found in non-Korean TGEV or PRCV isolates. Those nucleotide or amino acid sequences may be used to differentiate Korean TGEV strains from foreign TGEV strains although more TGEV isolates need to be analyzed. Paton and Lowings (1997) reported the nucleotide sequences (250-648 nucleotides) near the amino terminus of the S glycoprotein gene of several TGEV isolates, which are regions deleted in PRCV. TGEV Erica and Slagharen (Dutch isolates) had 15 nucleotide insertions between positions 496 and 497, and TGEV V63 (Belgium isolate) had 5 nucleotide insertions compared to other TGEV isolates (Paton and Lowings, 1997) . And TGEV 83-3289 (English isolate) had a deletion of six nucleotides between positions 286 and 291, which is a region known to be one of two B antigenic sites Gebauer et al., 1991) . The additions or deletions described above were not observed in Korean TGEV isolates. Sanchez et al. (1992) reported that the protein domain including antigenic subsites Aa and Ab and site D (Gebauer et al., 1991) showed more amino acid changes than other areas of the S glycoprotein. Six Korean TGEV-specific amino acids were found in that domain. A higher number of amino acid changes in Korean and non-Korean TGEV strains were fixed in that antigenic domain because it is involved in the neutralization of virus infectivity. TGEV isolates showed biologic and antigenic diversity (Vaughn and Paul, 1993) . Some TGEV isolates had cytopathic effect differences in swine testis cell culture and virus neutralization differences were observed when monoclonal antibodies were used in the virus neutralization. The virus neutralization test using polyclonal antibodies revealed only one serotype among the isolates as reported by Kemeny (1976) . Although the difference in amino acids forming antigenic sites were not found between Korean TGEV isolates and non-Korean TGEV strains, the virus neutralization test using monoclonal or polyclonal antibodies against TGEV need to be used to test antigenicity of Korean TGEV strains. Differences in the virulence of TGEV strains have been reported (Furuuchi et al., 1979; Cubero et al., 1992; Kim and Chae, 2002) . One Korean TGEV strain was less virulent than two American TGEV strains (Purdue and Miller) (Kim and Chae, 2002) . The relationship between nucleotide and amino acid sequence differences and virulence may be revealed when several TGEV strains are tested under identical conditions. The Korean TGEV KT2 strain was highly passaged in ST cells. The HKT2 strain, which was the KT2 strain passaged 104 times on ST cells, had several changes in nucleotide and amino acid sequences when compared to the parent TGEV KT2 strain. However, nucleotides or amino acids known to form four antigenic sites were not changed Gebauer et al., 1991) . The HKT2 strain had a deletion in genes 3a and 3b, which may be involved in virulence of TGEV (Woods, 1978; Wesley et al., 1990c) , and it was non-pathogenic in piglets (unpublished data). Therefore, the HKT2 strain may be used for a live vaccine to prevent TGE in piglets after testing the pathogenicity and immunogenicity in piglets and sows. In the phylogenetic tree, all Korean TGEV strains formed a different cluster from foreign TGEV or PRCV strains. Although Korean TGEV strains except for the KT2 and 133 strains formed one cluster, several sequence differences were found among Korean TGEV strains isolated in the same year. The Korean TGEV KT2 and 133 strains formed distinct clusters that were different from other Korean TGEV strains or non-Korean TGEV or PRCV strains. Recent European TGEV isolates were quite different from earlier TGEV types because they appear to originate from a different ancestor and result from accumulation of mutations but some isolates showed genetic stability over years (Paton and Lowings, 1997) . It has been suggested that reemergence of similar viruses into the field, laboratory cross-contamination, point mutations or recombination could affect the appearance of various TGEV strains (Banner and Lae, 1991; Paton and Lowings, 1997; Sanchez et al., 1999) . Those mechanisms may contribute to the emergence of TGEV strains with the differences in nucleotides and amino acids in Korea. In conclusion, four Korean TGEV field strains were different from published non-Korean TGEV or PRCV strains in nucleotide and amino acid sequences and were clustered into different groups from non-Korean TGEV or PRCV strains in phylogenetic tree. And four Korean TGEV isolates were clustered into three groups according to isolation years in phylogenetic tree. The highly passaged TGEV HKT2 strain had changes in S glycoprotein gene but did not contain a deletion in the S glycoprotein gene, which was found in PRCV. The pathogenicity and immunogenicity of the HKT2 strain are currently being tested in vivo. Random nature of coronavirus RNA recombination in the absence of selection pressure Sequence of the S gene from a virulent British field isolate of transmissible gastroenteritis virus The cloning and sequencing of the virion protein genes from a British isolate of porcine respiratory coronavirus, comparison with transmissible gastroenteritis virus genes Antigenic differentiation between transmissible gastroenteritis virus of swine and a related porcine respiratory coronavirus Cloning and sequencing of a 8.4-kb region from the 3 -end of a Taiwanese virulent isolate of the coronavirus transmissible gastroenteritis virus Antigenic structure of the E2 glycoprotein from transmissible gastroenteritis coronavirus Pathogenicity and antigen detection of the Nouzilly strain of transmissible gastroenteritis coronavirus, in 1-week-old piglets Antigenic structure of transmissible gastroenteritis virus. II. Domains in the peplomer glycoprotein Multiplication of low and high cell culture passaged strains of transmissible gastroenteritis virus in organs of newborn piglets Residues involved in the antigenic sites of transmissible gastroenteritis coronavirus S glycoprotein Antigenic variation of porcine transmissible gastroenteritis virus detected by monoclonal antibodies Sequence analysis of the porcine transmissible gastroenteritis coronavirus nucleocapsid protein gene Antibody response in pigs inoculated with transmissible gastroenteritis virus and cross reactions among ten isolates Experimental infection of piglets with transmissible gastroenteritis virus: a comparison of three strains (Korean, Purdue and Miller) Differentiation of infectious bronchitis virus serotypes using the polymerase chain reaction and restriction fragment length polymorphism analysis Sequence of the spike gene containing antigenic sites A and D of transmissible gastroenteritis virus isolated in Korea Sequence and N-terminal processing of the transmembrane protein E1 of the coronavirus transmissible gastroenteritis virus Discrimination between transmissible gastroenteritis virus isolates Isolation of a porcine respiratory, non-enteric coronavirus related to transmissible gastroenteritis The predicted primary structure of the peplomer protein E2 of the porcine coronavirus transmissible gastroenteritis virus Porcine respiratory coronavirus differs from transmissible gastroenteritis virus by a few genomic deletions Transmissible gastroenteritis Genetic evolution and tropism of transmissible gastroenteritis coronavirues Targeted recombination demonstrates that the spike gene of transmissible gastroenteritis coronavirus is a determinant of its enteric tropism and virulence Antigenic and biological diversity among transmissible gastroenteritis virus isolates of swine Three new isolates of porcine respiratory coronavirus with various pathogenicities and spike (S) gene deletions Nucleotide sequence of the E2-peplomer protein gene partial nucleotide sequence of the upstream polymerase gene of transmissible gastroenteritis virus (Miller strain) Evidence for a porcine respiratory coronavirus, antigenically similar to transmissible gastroenteritis virus, in the United States Genetic basis for the pathogenesis of transmissible gastroenteritis virus Small plaque variant transmissible gastroenteritis virus a Percentage of nucleotide similarity in upper triangle. b Percentage of amino acid similarity in lower triangle. The authors thank Dr. Chang-Won Lee for his excellent assistance to make phylogenetic tree and Dr. Mark W. Jackwood for editorial comments. This study was supported by Agricultural R&D Promotion Center, Korea.