Previous Article | Next Article ![]()
Journal of Virology, February 2003, p. 2789-2798, Vol. 77, No. 4
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.4.2789-2798.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
L. El-Attar,1 C. Gallimore,2 D. W. G. Brown,2 J. Green,2 and J. C. Bridger1*
Department of Pathology & Infectious Diseases, Royal Veterinary College, London NW1 OTU,1 Enteric and Respiratory Virus Laboratory, Central Public Health Laboratory, London NW9 5HT, United Kingdom2
Received 16 August 2002/ Accepted 12 November 2002
|
|
|---|
|
|
|---|
The United Kingdom reference BoCV, Newbury agent 2 (Bo/Newbury2/76/UK), was identified in southern England, near Newbury, in 1976 (48), and Dumfries virus (Bo/Dumfries/94/UK) was identified near Dumfries, Scotland, in 1994. Both were available as diarrheic fecal samples from experimentally infected calves (6; D. R. Snodgrass, personal communication) and contained calicivirus-like particles as shown by electron microscopy. Between 1998 and 2000, 476 bovine diarrheic fecal samples from dairy and beef farms were supplied from the Veterinary Laboratory Agency laboratories at Aberystwyth and Carmarthen (western Wales), Shrewsbury (central England), Penrith, Preston, and Thirsk (northern England), Starcross (southwestern England), and Winchester (southern England). Samples were stored at +4°C until tested and were then stored at -70°C. RNA was extracted from fecal samples by using 1,1,2-trichloro-trifluoroethane (Sigma-Aldrich Company Ltd.) and by a method modified from that of Boom et al., with guanidine isothiocyanate and silica (3). To identify positive samples, extracted RNA was reverse transcribed at 37°C for 60 min by using 200 IU of Moloney murine leukemia virus reverse transcriptase (Promega) and 25 pmol of random hexamer primers (Promega). The resulting cDNA was amplified by PCR by using 2 IU of HotStar Taq (Qiagen, Ltd.) in a 50-µl reaction volume containing 70 mM deoxynucleoside trisphosphates and 3 mM MgCl2 with 20 pmol of one of three primer pairs designed to detect human NLVs, namely, NI/E3 (10), GLPSG1/YGDD1, or GLPSG2/YGDD1 (15). A DNA Thermal Cycler (Perkin-Elmer) was programmed for 95°C for 15 min and 40 cycles of 94°C for 15 s, 40°C for 45 s, and 72°C for 1 min, with a final extension at 72°C for 10 min. PCR amplicons were resolved on 2.5% agarose gels. A positive signal for primer pair NI/E3 was a 114-bp amplicon, and positive signals for primer pairs GLPSG1/YGDD1 and GLPSG2/YGDD1 were 155-bp amplicons. Nuclease-free water was used as a negative control for every five fecal samples in both the reverse transcription (RT) and PCR steps and was handled simultaneously. Bo/Dumfries/94/UK was used as a positive virus control for the RT step. cDNA from Bo/Dumfries/94/UK was used as a positive control for the PCR. Positive samples were selected for long RT-PCR to produce amplicons of approximately 2.8 kbp from the 3' poly-A tail. Positive samples were selected by geographical location, by the primer pair used for the initial amplification from the fecal samples, or by the relationship in a heteroduplex mobility assay (30; paper in preparation) to Bo/Newbury2/76/UK and two human NLVs currently circulating in the United Kingdom, Hu/Grimsby/95/UK and Hu/RHB/93/UK. RT was conducted at 42°C with the primer NA-2 poly-T EcoRI (GCGCTTAAGT(15)NMCAATGC) (LifeTechnologies Ltd.) by using 200 IU of Moloney murine leukemia virus (H-) reverse transcriptase (Promega). PCR was conducted by using the high-fidelity Z-Taq polymerase (BioWhittaker) and primer pairs NI/NA-2 poly-T EcoRI or GLPSG1/NA-2 poly-T EcoRI. Thirty-five cycles of 98°C for 5 s, 50°C for 10 s, and 72°C for 1 min were performed, with final extension at 72°C for 5 min, and PCR amplicons resolved on 1% agarose gels. The
2.8-kbp amplicons were gel purified and cloned by using a TOPO-TA XL-PCR cloning kit (Invitrogen Ltd.). Cycle sequencing was performed by using a dye terminator cycle sequencing kit and a CEQ 2000 XL DNA analysis system (Beckman Coulter). Nucleotide sequences for the polymerase, capsid, and open reading frame 3 (ORF3) genes were derived from cloned amplicons. Consensus nucleotide sequences were prepared by using the STADEN sequencing package (43). Multiple alignments of nucleotide and translated amino acid sequences were prepared by using Clustal X version 1.8 (44), and identities were calculated by using GeneDoc (www.psc.edu/biomed/genedoc). Open reading frames (ORFs) were predicted by using the NCBI ORF finder. Phylogenetic analyses were performed with nucleotide and amino acid sequence alignments by parsimony (DNApars and Protpars), Fitch-Margoliash and bootstrap analyses were performed by using PHYLIP (J. Felsenstein, Department of Genetics, University of Washington, Seattle; Phylogeny Inference Package, version 3.5c), and maximum likelihood with quartet puzzling was determined by using TreePuzzle 5.0 (38). Phylogenetic trees were prepared by using TreeView (35). Nucleotide and amino acid sequences of the polymerase, capsid, and ORF3 genes of the UK BoCVs were subjected to a GenBank BLAST search. An amino acid identity plot for the capsid proteins was generated by using SimPlot (29) with a window of 80 amino acids at 20-amino-acid intervals.
Nucleotide sequence accession numbers. The nucleotide sequence accession numbers referred to in the text are human genogroup I NLVs: Hu/NLV/Norwalk/68/US, M87661; Hu/NLV/Chiba/87/JP, AB042808; Hu/NLV/Musgrove/89/UK,AJ277614; Hu/NLV/Seto124/89/JP, AB031013; Hu/NLV/Desert Shield395/90/SA, U04469; Hu/NLV/Thistlehall/90/UK, AJ277621; Hu/NLV/Southampton/91/UK, L07418; Hu/NLV/Birmingham/93/UK, AJ277612; Hu/NLV/Winchester/94/UK, AJ277609; Hu/NLV/Sindlesham/95/UK, AJ277615; Hu/NLV/Valetta/95/Malta, AJ277616; Hu/NLV/Whiterose/96/UK, AJ277610; and Hu/NLV/Hesse3/97/DE, AF093797. The following nucleotide sequence accession numbers are human genogroup II NLVs: Hu/NLV/Hawaii/71/US, U07611; Hu/NLV/Snow Mountain/76/US, U75682; Hu/NLV/Melksham/89/UK, X81879; Hu/NLV/Mexico/89/MX, U22498; Hu/NLV/Hillingdon/90/UK, AJ277607; Hu/NLV/Leeds/90/UK AJ277608; Hu/NLV/Seacroft/90/UK, AJ277620; Hu/NLV/Wortley/90/UK, AJ277618 (given as genogroup I in the GenBank database but genogroup II by J. Green et al. (11); Hu/NLV/Toronto24/91/CAN, U02030; Hu/NLV/Girlington/93/UK, AJ277606; Hu/NLV/Lordsdale/93/UK, X86557; Hu/NLV/Rbh/93/UK, AJ277617 (given as genogroup I in the GenBank database but genogroup II by J. Green et al. (11); Hu/NLV/Camberwell/94/AUS, U46500; Hu/NLV/Bham132/95/UK, AJ277611; Hu/NLV/Grimsby/95/UK, AJ004864; Hu/NLV/Mendoza320/95/ARG, AF190817; Hu/NLV/Parkroyal/95/UK, AJ277613; Hu/NLV/Symgreen/95/UK, AJ277619; Hu/NLV/Alphatron/98-2/1998/NET, AF195847; and Hu/NLV/Amsterdam/98-18/1998/NET, AF195848. The accession number for sapovirus Hu/SLV/Manchester/93/UK is X86560. Accession numbers for BoCVs reported previously are as follows: Bo/Newbury2/76/UK, AF097917 and AF097918; Bo/Jena/80/DE, AJ011099; Bo/NLV/176/98/NET, AF194183; Bo/NLV/CH126/98/NET, AF320625; and Bo/NLV/CH131/98/NET, AF320113. Accession numbers for the BoCVs reported in the present paper are as follows: Bo/Aberystwyth24/00/UK, AY126475; Bo/Aberystwyth58/00/UK, AY126462; Bo/Aberystwyth65/00/UK, AY126463; Bo/Carmarthen42/98/UK, AY126459;Bo/Carmarthen10/99/UK, AY126460; Bo/Carmarthen28/99/UK, AY126461; Bo/Dumfries/94/UK, AY126474; Bo/Newbury2/76/UK, AF097917; Bo/Penrith9/00/UK, AY126469; Bo/Penrith24/00/UK, AY126470; Bo/Penrith33/00/UK, AY126464; Bo/Penrith41/00/UK, AY126471; Bo/Penrith55/00/UK, AY126476; Bo/Shrewsbury24/00/UK, AY126472; Bo/Shrewsbury39/00/UK, AY126465; Bo/Shrewsbury46/00/UK, AY126466; Bo/Shrewsbury54/00/UK, AY126473; Bo/StarCross37/00/UK, AY126467; and Bo/Thirsk10/00/UK, AY126468.
Analysis of the polymerase, capsid and ORF3 gene sequences of the UK reference strain, Bo/Newbury2/76/UK. A 581-nucleotide fragment of the Bo/Newbury2/76/UK polymerase gene (from nucleotides 4482 to 5064 of the previously described Bo/Jena/80/DE sequence [28]) had 89% amino acid (77% nucleotide) identity to Bo/Jena/80/DE isolated in Germany in 1980. The complete Bo/Newbury2/76/UK capsid gene (nucleotides 5051 to 6610 of the Bo/Jena/80/DE sequence) (Fig. 1) had 68% amino acid (67% nucleotide) identity to Bo/Jena/80/DE, which had two deletions of 3 and 6 nucleotides at positions 871 to 873 (amino acid 291) and 1270 to 1275 (amino acids 424 and 425) compared to Bo/Newbury2/76/UK. The complete Bo/Newbury2/76/UK ORF3 gene (excluding the extended amino terminus [see below]) had 63% amino acid (59% nucleotide) identity to Bo/Jena/80/DE and two deletions of 18 and 3 nucleotides at positions 355 to 372 (amino acids 119 to 124) and 562 to 564 (amino acid 188) compared to the Bo/Jena/80/DE ORF3 gene. The overlap of the Bo/Newbury2/76/UK ORF1 and ORF2 genes was identical in size to that of Bo/Jena/80/DE (nucleotide positions 5051 to 5064) (Fig. 2; Table 1), but the overlap of the Bo/Newbury2/76/UK ORF2 and ORF3 genes was considerably longer (nucleotide positions 6392 to 6610). However, visual inspection showed the presence of a second AUG codon in the Bo/Newbury2/76/UK sequence downstream of the predicted ORF3 initiation codon and in an analogous location to the initiation codon found in Bo/Jena/80/DE. The Bo/Newbury2/76/UK capsid and ORF3 genes were longer than that of Bo/Jena/80/DE, but the 3' nontranslated region was shorter.
![]() View larger version (67K): [in a new window] |
FIG. 1. Multiple amino acid sequence alignment of the four UK BoCV capsid proteins, the German Bo/Jena/80/DE, and a UK representative of human NLV genogroups I (Hu/NLV/Southampton/91/UK) and II (Hu/NLV/Melksham/89/UK). The dotted line beneath the sequences represents amino acids 50 to 225 of the S domain, the solid lines represent the P1 domain, and the dashed line represents the P2 domain of the Hu/Norwalk/68/US (36). The conserved 5-amino-acid LAGNA motifs located at the 3- and 5-fold axes of symmetry of Hu/Norwalk/68/US are highlighted in gray.
|
![]() View larger version (36K): [in a new window] |
FIG. 2. Schematic representation and nucleotide sequence of the reference UK BoCV (Bo/Newbury2/76/UK), the German BoCV (Bo/Jena/80/DE), the human NLV genogroup I Hu/NLV/Norwalk/68/US, and the human NLV genogroup II Hu/NLV/Melksham/89/UK viruses in the region of the ORF1-2 and ORF2-3 overlaps, showing start (in bold) and stop (underlined) codons as predicted by the NCBI ORF finder. The boxes on the left depict the 3' ends of the ORF1s, the pairs of black boxes in the center depicts the 5' and 3' ends of the ORF2s, and the boxes on the right depict the complete ORF3s. The second initiation codon of the Bo/Newbury2/76/UK sequence is shown (bold italicized). The 173 bases not shown in the Bo/Newbury2/76/UK ORF2-3 overlap are labeled 173. The scale indicates nucleotide positions in relation to Bo/Jena/80/DE.
|
|
View this table: [in a new window] |
TABLE 1. Comparison of the genomic organizations of the UK BoCVs sequenced in the present study, the German Bo/Jena/80/DE, and human genogroup I and II NLVs, as predicted by the NCBI ORF finderk
|
Partial polymerase sequences of 76 nucleotides from 18 samples from geographically diverse locations (Aberystwyth24/00/UK, Aberystwyth58/00/UK, Aberystwyth65/00/UK, Carmarthen42/98/UK, Carmarthen10/99/UK, Carmarthen28/99/UK, Dumfries/94/UK, Penrith9/00/UK, Penrith24/00/UK, Penrith33/00/UK, Penrith41/00/UK, Penrith55/00/UK, Shrewsbury24/00/UK, Shrewsbury39/00/UK, Shrewsbury46/00/UK, Shrewsbury54/00/UK, Starcross37/00/UK, and Thirsk10/00/UK), including the samples from the 4-month-old calf and the two cows, showed a close relationship to each other, to Bo/Newbury2/76/UK, and to Bo/Jena/80/DE virus (28) (92 to 100% amino acid and 68 to 100% nucleotide identity). Nucleotide diversity showed no correlation with year of isolation, age, or geographical location within the United Kingdom. Longer partial polymerase sequences of 416 nucleotides of 7 samples, namely, Aberystwyth24/00/UK, Dumfries/94/UK, Penrith9/00/UK, Penrith24/00/UK, Penrith55/00/UK, Shrewsbury24/00/UK, Shrewsbury54/00/UK (identified using the NI/E3 primer pair), and Penrith41/00/UK (identified with the GLPSGI/YGDDI primer pair), also showed a close relationship to each other, to Bo/Newbury2/76/UK, and to Bo/Jena/80/DE (88 to 100% amino acid and 74 to 99% nucleotide identity). In addition, the 61 nucleotides (20 translated amino acids) available from the polymerase gene of the Dutch BoCV, Bo/NLV/176/98/NET (45), which overlapped with the sequences obtained in the present study, had 90 to 100% amino acid (70 to 100% nucleotide) identity to the UK BoCVs.
The complete capsid genes of Newbury2/76/UK, Dumfries/94/UK, Aberystwyth24/00/UK, and Penrith55/00/UK had 95 to 97% amino acid (87 to 89% nucleotide) identity to each other and to two unpublished BoCV sequences from The Netherlands (Bo/NLV/CH126/98/NET and Bo/NLV/CH131/98/NET) but were less closely related to the German Bo/Jena/80/DE virus (68% amino acid identity and 66 to 67% nucleotide identity). Alignment of translated amino acid sequences for the conserved shell (S) domain (36) showed total conservation between the UK BoCVs, apart from one amino acid substitution in the Aberystwyth24/00/UK sequence (Fig. 1), but showed lower identity to the German Bo/Jena/80/DE virus (Fig. 3). In the protruding (P) domain (36), single amino acid substitutions and pairs of substitutions were found between the UK BoCVs which had 93 to 96% amino acid (85 to 88% nucleotide) identities in this region, but substantial variation occurred between the UK BoCVs and Bo/Jena/80/DE (Fig. 3). The complete ORF3 genes (starting from the Bo/Jena/80/DE ORF3 initiation codon) of the four UK BoCVs, Newbury2/76/UK, Dumfries/94/UK, Aberystwyth24/00/UK, and Penrith55/00/UK, had 92 to 96% amino acid (87 to 90% nucleotide) identity to each other but had 61 to 63% amino acid (58 to 59% nucleotide) identity to Bo/Jena/80/DE.
![]() View larger version (22K): [in a new window] |
FIG. 3. Amino acid identity plot generated by using SimPlot (29) with a window of 80 amino acids at 20-amino-acid intervals to compare the complete capsid protein of Bo/Newbury2/76/UK with the capsid protein of Bo/Jena/80/DE and the genogroup I human NLV Hu/Southampton/91/UK and the genogroup II NLV Hu/Melksham/89/UK. The dotted line beneath the plot represents the S domain, the solid line represents the P1 domain, and the dashed line represents the P2 domain of the Hu/NLV/Norwalk/68/US, as determined by X-ray crystallography (36).
|
Comparison of BoCVs and the human NLVs. Analysis of the polymerase (76 and 416 nucleotides), complete capsid, and ORF3 gene sequences of the UK BoCVs, plus the German Bo/Jena/80/DE virus, showed that the BoCVs were distinct from representative human genogroup I and II NLVs and the NLVs currently associated with gastroenteritis in the UK human population from 1989 to 2000 (11; personal observations). The highest identities between the BoCVs and the human NLVs were observed with the polymerase gene sequences, but the identities seen between the BoCVs and genogroup I and II NLVs (Table 2) were comparable to those seen between genogroup I and II NLV, which had 62 to 68% amino acid identity in the 416-nucleotide region examined in the present study. The amino acid identities of the complete BoCV capsid proteins to representatives of the seven genetic clusters of human NLV genogroups I and II (41 to 48%) were similar to the identities between the two human NLV genogroups (41 to 46%) (Table 3). The greatest divergence between the BoCVs and the human NLVs was seen with the BoCV ORF3 genes (Table 2). A BLAST search with UK BoCV sequences from all three genes failed to identify human NLVs with identities greater than those between the UK BoCVs and the representative human genogroup I and II NLVs described above. Phylogenetic analyses of the nucleotide and translated amino acid sequences of all three genes by the independent methods of parsimony, distance, and maximum likelihood each showed that the UK BoCVs, the German Bo/Jena/80/DE, and two unpublished capsid sequences from bovine viruses identified in The Netherlands (Bo/NLV/CH126/98/NET and Bo/NLV/CH131/98/NET) formed a distinct genogroup with quartet puzzling support values of 100% for the polymerase and capsid sequences and 76% for the ORF3 sequences (Fig. 4 and 5). Heterogeneity between the BoCVs and the human NLVs was distributed throughout the capsid protein but was concentrated in the P2 domain (Fig. 1 and 3). A 5-amino-acid motif, LAGNA, was conserved in the S domain of the UK BoCVs, the German Bo/Jena/80/DE, 2 Dutch BoCVs (Bo/CH126/1998/NL and Bo/CH131/1998/NL), and 21 genogroup I and 49 genogroup II human NLVs but was absent in the caliciviruses classified in the Vesivirus, Lagovirus and Sapporo-like genera. By using the molecular imaging program RASMOL, the LAGNA motif was located at the 3- and 5-fold axes of symmetry of Hu/Norwalk/68/US, suggesting structural importance (not shown). The UK BoCVs had an overall genomic organization similar to that of the human NLVs but with some consistent differences (Table 1; Fig. 2). The BoCV-predicted ORF1-2 overlaps were at least 3 nucleotides shorter, but the ORF2-3 overlaps were considerably longer than those of the human NLVs because of the additional initiation codon (see above). The termination codon for the UK and German BoCV ORF2 genes was UGA and not the UAA termination codon which has been commonly observed with human NLVs. The UK BoCVs capsid proteins were at least 8 amino acids shorter, their ORF3 genes were longer, and their 3' untranslated regions were shorter than those of any known human NLV.
|
View this table: [in a new window] |
TABLE 2. Percent amino acid (nucleotide) identities of BoCV partial polymerase and ORF3 genes to representative human genogroup I and II NLVs
|
|
View this table: [in a new window] |
TABLE 3. Percent amino acid identities between the complete capsid proteins of BoCVs and the seven genetic clusters of human genogroup I and II NLVs
|
![]() View larger version (32K): [in a new window] |
FIG. 4. Phylogenetic analysis generated by using the TreePuzzle maximum likelihood method for the translated amino acid sequences of (A) the long polymerase sequence (138 amino acids) and (B) 113 amino acids from the 5' end of the ORF3 gene (starting from the Bo/Jena/80/DE ORF3 initiation codon) of the UK BoCVs, the German Bo/NLV/Jena/80/DE, and representatives of the human NLV genogroups I and II, with the Sapporo-like virus Hu/Manchester/93/UK used as the outgroup. The numbers at the nodes (quartet puzzling values) indicate the frequencies of occurrence for 1,000 replicate trees. The asterisk indicates the sample identified with the GLPSG1/YGDD1 primer pair.
|
![]() View larger version (30K): [in a new window] |
FIG. 5. Phylogenetic analysis generated by using the TreePuzzle maximum likelihood method for the translated amino acid sequences of the complete capsid protein gene, the UK BoCVs, the German Bo/Jena/80/DE, the two unpublished Dutch BoCVs Bo/CH126/98/NET and Bo/CH131/98/NET, and representatives of the seven clusters in the human NLV genogroups I and II, with the Sapporo-like virus Hu/Manchester/93/UK used as the outgroup. The numbers at the nodes (quartet puzzling values) indicate the frequencies of occurrence for 1,000 replicate trees. The superscripts represent the genetic clusters of the genogroup I and II viruses, as determined by Ando et al. (1).
|
The choice of primer pairs used in the present study was deliberate in order to identify the BoCVs most closely related to human NLVs. The three sets of primer pairs are currently used successfully for diagnostic RT-PCR on human diarrheic fecal samples at the Central Public Health Laboratory by the authors. We may have detected a subset of BoCVs, but it seems unlikely that the use of additional primer pairs would identify BoCVs with a closer relationship to human NLVs. We have not detected all BoCVs circulating in the United Kingdom cattle population, as we have failed to amplify Newbury agent 1 (4, 7; personal observations). Other primer pairs will be needed to identify the complete range of BoCVs circulating in cattle, but we predict that these yet undefined BoCVs are unlikely to be closely related to the human NLVs.
The close relationship of the BoCVs identified in 1976, 1994, and between 1998 and 2000 in the United Kingdom and their close relationship to the German (28) and Dutch BoCVs (45) identified in 1980 and 1998 contrasted with the diversity identified with human NLVs over a similar time period (reviewed in references 1 and 13). However, there have been instances of closely related human NLVs identified many years apart from different continents (14, 19). We have no evidence to indicate that laboratory contamination accounted for the close relationship between the UK BoCVs. All RTs and PCRs were strictly controlled with negative samples, and no cross-contamination occurred. Although little heterogeneity was identified between the amino acid sequences, it should be noted that considerable heterogeneity was identified between some UK BoCVs at the nucleotide level. Ando et al. (1) suggested that BoCVs might be separated into two clusters represented by viruses similar to Bo/Newbury2/76/UK and Bo/Jena/80/DE. The GLSPG1/YGDD1 primer pair, which was used by Liu et al. (28) to identify Bo/Jena/80/DE, detected several UK BoCVs, but we failed to amplify viruses closely related to Bo/Jena/80/DE. The 68% amino acid identity of the capsid protein of Bo/Jena/80/DE to that of the UK BoCVs was comparable to the relationship of human NLVs within the two human genogroups (Table 3), confirming that the Bo/Jena/80/DE virus should be classified in the bovine genogroup III. BoCVs caused enteric disease in experimental calves (4, 48), caused pathology in the small intestine (18), and were found in association with outbreaks of calf diarrhea in the United Kingdom (37). However, van der Poel et al. (45) found BoCVs commonly in veal calves, highlighting the need for further studies to determine their role in bovine diarrhea outbreaks.
We thank the Biotechnology and Biological Sciences Research Council (BBSRC) for studentship support to S.L.O.
Present address: Merial Asia Pte. Ltd., Singapore 118238, Republic of Singapore. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»