Previous Article | Next Article ![]()
Journal of Clinical Microbiology, November 2004, p. 5214-5224, Vol. 42, No. 11
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.11.5214-5224.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Christopher Thomas, and
Linda J. Saif*
Food Animal Health Research Program, Ohio Agricultural Research and Development Center, The Ohio State University, Wooster, Ohio
Received 20 February 2004/ Returned for modification 1 July 2004/ Accepted 7 July 2004
|
|
|---|
|
|
|---|
The enteric caliciviruses include the noroviruses and sapoviruses, which cause acute gastroenteritis in humans and animals (3, 7, 9, 11, 15, 38). Bovine enteric caliciviruses (BECVs) are classified as a third genogroup (GIII) of noroviruses, distinct from GI and GII human noroviruses (1, 5, 23, 29). They include Bo/NVL/Newbury-2/76/UK and Bo/NLV/Jena/80/DE (1). Recently, an unclassified BECV, Nebraska (NB) strain, which is most closely related to Sapovirus and Lagovirus, was recognized as a new NB-like genus or genogroup (34). The capsid genes of GIII noroviruses are composed of 1,560 to 1,569 bp, whereas the NB strain capsid gene is 1,650 bp, which is 81 to 90 bp larger than those of GIII noroviruses (23, 29, 35).
Determination of the antigenic relationship between genogroups and genotypes of enteric caliciviruses is difficult because most viruses, except for Cowden porcine sapovirus, do not grow in cell culture (8, 30). Based on genetic divergence in the RNA-dependent RNA polymerase (RdRp) and capsid genes, Norovirus and Sapovirus are each classified into at least three distinct genogroups, with multiple genotypes, which is indicative of the genetic diversity found within each genogroup (1, 10, 13, 36). However, the antigenic relationships among the genogroups of enteric caliciviruses have not been extensively analyzed. Calicivirus capsid proteins expressed in a baculovirus expression system self-assemble into virus-like particles that are antigenically and morphologically similar to native virions. The virus-like particles of human noroviruses have been used as antigens in enzyme immunoassays for serotyping of human noroviruses (14, 18).
Characterized bovine noroviruses (BoNVs) comprise two genotypes (1 and 2) in Norovirus GIII, of which Bo/Jena/80/DE and Bo/Newbury-2/76/UK are the prototype viruses of GIII/1 and GIII/2, respectively. With the exception of Bo/Jena/80/DE, most other recently reported BoNVs, such as Bo/Aberystwyth24/00/UK, Bo/CH126/98/NET, Bo/CH131/98/NET, Bo/CV95-OH/00/US, Bo/CV186-OH/00/US, Bo/Dumfries/94/UK, and Bo/Penrith55/00/UK, are classified as Norovirus GIII/2 (like Bo/Newbury-2/76/UK) based on the determined sequences of their RdRp and/or capsid genes (1, 29, 35). To further characterize the genetic diversity of genogroups and genotypes of BECVs circulating among calves and to determine their genetic relationships to human noroviruses and sapoviruses, we sequenced the complete capsid genes of 21 BECV strains detected in veal calves in Ohio.
|
|
|---|
RNA extraction and RT-PCR. Viral RNA was extracted from fecal samples diluted 1:10 with phosphate-buffered saline (pH 7.4) by using the TRIzol reagent (Invitrogen Corporation, Carlsbad, Calif.) according to the manufacturer's instructions. RNA pellets were suspended in 50 µl of nuclease-free water and used to amplify the capsid gene by RT-PCR. The remainder of the RNA was stored at 70°C for future use.
Primers used to amplify the complete capsid gene were designed based on previously reported sequence data (5, 23, 34). The primer sets used to amplify the complete capsid gene of BECV were NAcap-F/R and Jenacap-F1/F2/R1/R2/R3 for GIII noroviruses and NBcap-F3/R for NB-like BECV (Table 1). Primers NAcap-F/R and Jenacap-F1/F2/R1/R2/R3 were designed based on published sequences of known GIII/2 BoNVs, including Bo/Newbury-2/76/UK (accession number AF097917) and Bo/Jena/80/DE (AJ011099), respectively. Primer set NBcap-F3/R was designed from the gene sequence of the NB strain (NC_004064). To amplify the partial RdRp and capsid genes of Bo/CV521-OH/02/US, primer RC35/1932 was designed based on its previously reported RdRp gene sequence (AY151257) and the capsid gene sequence determined in this study (AY549161). Internal primers were designed based on the conserved region of the capsid gene sequences of BECV strains characterized in this study and were used to sequence the internal regions of the capsid genes of noroviruses and NB-like BECVs.
|
View this table: [in a new window] |
TABLE 1. Newly designed oligonucleotide primer sequences used for amplification of the capsid genes of BoNVs and NB-like BECVs
|
Cloning and sequencing of the capsid gene. The amplified RT-PCR products were purified with a gel extraction kit (Qiagen Inc., Valencia, Calif.) and cloned into PCR2.1 vectors by using the TOPO TA cloning kit (Invitrogen Corporation) according to the manufacturer's instructions. Plasmid DNA from transformed cells was extracted by using an alkaline lysis method (Qiagen Inc.), and the presence of BECV capsid gene inserts was confirmed by restriction enzyme mapping. Cloned capsid genes were sequenced by using an ABI Prism 377 DNA sequencer (Applied Biosystems, Foster, Calif.). Two to four clones for each BECV strain were sequenced to ascertain the consensus sequence.
Phylogenetic analysis. Multiple alignments of nucleotide and amino acid sequences were performed with Clustal W as implemented in the Lasergene software package (DNASTAR Inc., Madison, Wis.), and nucleotide and amino acid identities were calculated by using the same package. Phylogenetic analyses of Clustal W nucleotide and amino acid alignments generated by the European Bioinformatics Institute (http://www.ebi.ac.uk) were performed by the UPGMA and neighbor-joining methods, and 1,000 bootstrap replicates generated by using the Molecular Evolutionary Genetics Analysis version 2.1 (MEGA2) software (21). Genetic distances among human, porcine, and bovine noroviruses and among Norovirus genogroups and genotypes were calculated by Kimura's two-parameter method (20).
Gene conversion analysis. Genetic recombination was analyzed by using the Recombinant Identification Program (RIP) (33) and the sliding-window genetic diversity plot (SimPlot), version 3.2 beta, which was kindly provided by S. Ray, Johns Hopkins University School of Medicine, Baltimore, Md. (25). Bo/Newbury-2/76/UK and Bo/Jena/80/DE were used for background sequences. The analyzed ORF1 C-terminal region (partial RdRp gene) and ORF2 of Bo/CV521-OH/02/US, whose size was 1,898 bp, corresponded to nucleotides 1651 to 2232 of Bo/Newbury-2/76/UK and to nucleotides 4601 to 6489 of Bo/Jena/80/DE. The aligned sequences were analyzed by RIP with a window size of 100 bp and a threshold for statistical significance of 90% and by SimPlot with a window 200 bp wide and a step size of 20 bp.
Nucleotide sequence accession numbers. The sequences of the complete capsid genes of BECV strains have been assigned GenBank accession numbers. The accession numbers for the BoNVs are as follows: Bo/CV490-OH/02/US, AY549153; Bo/CV499-OH/02/US, AY549154; Bo/CV500-OH/02/US, AY549155; Bo/CV506-OH/02/US, AY549156; Bo/CV509-OH/02/US, AY549157; Bo/CV510-OH/02/US, AY549158; Bo/CV511-OH/02/US, AY549159; Bo/CV511 M-OH/02/US, AY549160; Bo/CV521-OH/02/US, AY549161; Bo/CV533-OH/02/US,AY549162; Bo/CV540-OH/02/US, AY549163; Bo/CV550-OH/02/US, AY549164; Bo/CV551-OH/02/US, AY549165; Bo/CV551 M-OH/02/US, AY549166; and Bo/CV557-OH/02/US, AY549167. The accession numbers for NB-like BECVs are as follows: Bo/CV504-OH/02/US, AY549168; Bo/CV519-OH/02/US, AY549169; Bo/CV526-OH/02/US, AY549170; Bo/CV531-OH/02/US, AY549171; Bo/CV548-OH/02/US, AY549172; and Bo/CV562-OH/02/US, AY549173.
|
|
|---|
Sequence identity and phylogenetic analysis of the capsid genes of BoNV and NB-like BECV. Multiple alignments and phylogenetic analysis of the capsid gene sequences of the BECVs sequenced in this study and of genetically characterized human, bovine, and porcine noroviruses as well as human, porcine, and mink sapoviruses were performed. Of the 21 BECV strains sequenced, 15, including Bo/CV511-OH/02/US, Bo/CV511 M-OH/02/US, Bo/CV551-OH/02/US and Bo/CV551 M-OH/02/US, belonged to Norovirus GIII/2, which includes most known BoNVs, except for Bo/Jena/80/DE (GIII/1) (Fig. 1 and Table 2). Six strains were classified as NB-like BECV (Fig. 2). The capsid genes of all 15 characterized BoNVs were comprised of 1,569 bp and were predicted to encode 522 aa. Each of the six NB-like BECV capsid genes characterized was 1,686 bp in size, which corresponds to the size determined for the BECV NB strain (34).
![]() View larger version (25K): [in a new window] |
FIG. 1. Phylogenetic tree, constructed by the neighbor-joining method, based on the amino acid sequence of the entire capsid gene of BoNV. The bootstrap value (percent) is given at each node. Viruses for which capsid gene sequences were used for phylogenetic analysis (GenBank accession numbers) included (i) GI human noroviruses Hu/NLV/Chiba407/87/JP (AB042808), Hu/NLV/Desert Shield395/90/SA (U04469), Hu/NLV/Hesse3/97/DE (AF093797), Hu/NLV/Musgrove/89/UK (AJ277614), Hu/NLV/Norwalk/68/US (M87661), Hu/NLV/Sindlesham/95/UK (AJ277615), Hu/NLV/Southampton/91/UK (L07418), and Hu/NLV/Winchester/94/UK (AJ277609); (ii) GII human noroviruses Hu/NLV/Bristol/93/UK (X76716), Hu/NLV/Camberwell1101922/94/AUS (AF145896), Hu/NLV/Hawaii/71/US (U07611), Hu/NLV/Hillingdon/90 (AJ277607), Hu/NLV/Leeds/90/UK (AJ277608), Hu/NLV/Lordsdale/93/UK (X86557), Hu/NLV/Melksham/89/UK (X81879), Hu/NLV/Mexico/89/MX (U22498), Hu/NLV/Seacroft/90/UK (AJ277620), Hu/NLV/Toronto24/91/CAN (U02030), Hu/NLV/Snow Mountain/76/US (U70059), and Hu/NLV/VA97207/97/US (AY038599); (iii) GII porcine noroviruses Sw/NLV/Sw43/97/JP (AB074892) and Sw/NLV/Sw918/97/JP (AB074893); and (iv) GIII bovine noroviruses: Bo/Aberystwyth24/00/UK (AY126475), Bo/CH126/98/NET (AF320625), Bo/CH131/98/NET (AF320113), Bo/CV95-OH/00/US (AF542083), Bo/CV186-OH/00/US (AF542084), Bo/Dumfries/94/UK (AY126474), Bo/Jena/80/DE (AJ011099), Bo/Newbury-2/76/UK (AF097917), and Bo/Penrith55/00/UK (AY126476).
|
|
View this table: [in a new window] |
TABLE 2. Capsid gene nucleotide and amino acid identities between BoNV and representative human noroviruses
|
![]() View larger version (24K): [in a new window] |
FIG. 2. Phylogenetic tree, constructed by the neighbor-joining method, based on the amino acid sequence of the entire capsid gene of NB-like BECV. The bootstrap value (percent) is given at each node. Viruses for which capsid gene sequences were used for phylogenetic analysis (GenBank accession numbers) included (i) GI human sapoviruses Hu/SLV/Houston/86/US (U95643), Hu/SLV/Houston/27/90/US (U95644), Hu/SLV/Manchester/93/UK (X86560), Lyon30388/98 (AJ251991), Hu/SLV/Parkville/94/US (U73124), Hu/SLV/Plymouth/92/UK (X86559), Hu/SLV/Sapporo/82/JP (U65427), and Hu/SLV/Stockholm/318/97/SE (AF194182); (ii) GII human sapoviruses Hu/SLV/Bristol/98/UK (AJ249939), Hu/SLV/cruise ship/00/US (AY289804), Hu/London/29845/92/UK (U95645), and Hu/SLV/Lyon/598/97/F (AJ271056); (iii) porcine sapovirus Sw/Cowden//UA (AF182760); (iv) mink sapovirus Canada 151A (AY144337); and (v) unclassified NB-like BECV Bo/NB/80/US (NC_004064).
|
The genetic distances based on the amino acid sequences among the 15 BoNVs were 0.008 to 0.082 (0.032 ± 0.019). The 15 BoNVs had genetic distances of 0.015 to 0.073 (0.035 ± 0.012) from GIII/2 bovine noroviruses and of 0.344 to 0.367 (0.355 ± 0.006) from the GIII/1 bovine norovirus, Bo/Jena/80/DE, compared with 0.659 ± 0.014 and 0.761 ± 0.023 from GI and GII noroviruses, respectively. The genetic distances between Norovirus GI and GII, GI and GIII, and GII and GIII were 0.733 ± 0.025, 0.655 ± 0.013, and 0.755 ± 0.022, respectively. Based on sequence identities and genetic distances, BoNVs were more closely related to GI noroviruses than to GII noroviruses.
The motifs IDPWI, LAGNA, and LYTP, located in the S domain, and QNGR, located in P1A subdomain, were conserved in all noroviruses, including the 15 BoNVs. The motif LAGNA, located at the three- and fivefold axes of symmetry of Hu/NLV/Norwalk/68/US, has been suggested to be of structural importance (29). The stop codon of the 15 BoNVs was TGA, like that of other BoNVs and human GI and GII noroviruses.
Phylogenetic analysis of NB-like BECV. The amino acid sequences of six NB-like BECVs shared 94 to 99.5% identity (87 to 99.8% nucleotide identity) with each other and 93 to 99% identity (88 to 92% nucleotide identity) with the NB strain. The amino acid identity of the six NB-like BECVs to BoNVs, human noroviruses, and sapoviruses, including porcine and mink sapoviruses, was less than 23%. The genetic distances of NB-like BECVs based on the amino acid sequence were 0.002 to 0.058 (0.027 ± 0.024) from each other and 0.08 to 0.013 (0.01 ± 0.002) from the NB strain. However, the genetic distances of NB-like BECVs were higher from human noroviruses and human, mink, and porcine sapoviruses (1.343 to 1.765). Bo/CV519-OH/02/US showed the lowest amino acid (93%) and nucleotide (87%) identities and also the highest genetic distance (0.56 to 0.06) within the six NB-like BECVs.
Comparative phylogenetic analysis of the whole capsid protein and the S and P domains of BoNV. The nucleotide and amino acid identities of the capsid genes of the 15 BoNVs, other characterized BoNVs, and Hu/NLV/Norwalk/68/US are summarized in Table 3. Analyses were conducted on the S (176 aa), P1A (53 aa), P2 (127 aa), and P1B (117 aa) domains, corresponding to those described for Hu/NLV/Norwalk/68/US (32). The S domain was the most conserved region between noroviruses. The 15 BoNVs had the lowest amino acid identity between both genogroups and genotypes in the P2 domain, which includes the hypervariable region of the capsid gene. The amino acid identity of the S domain of 15 BoNVs was highest with Bo/Newbury-2/76/UK. The amino acid identity for the S domain with Bo/Jena/80/DE (Table 3) was higher at 85 to 90% (76 to 78% nucleotide identity), compared to 67 to 69% (68 to 70% nucleotide) identity for the whole capsid gene (Table 3).
|
View this table: [in a new window] |
TABLE 3. Comparison of amino acid and nucleotide sequence identities of the S and P domains of the capsid region in bovine and human noroviruses
|
![]() ![]() View larger version (35K): [in a new window] |
FIG. 3. Comparative phylogenetic trees, constructed by the unweighted-pair-group-method-with-arithmetic method, based on amino acid sequence of bovine noroviruses and representative human noroviruses. (A) S domain; (B) P2 domain. Bootstrap values of greater than 70% are given at each node.
|
From the results of RIP analysis, Bo/CV521-OH/02/US showed higher similarity with Bo/Jena/80/DE in the RdRp region and, in contrast to RdRp, with Bo/Newbury-2/76/UK in the capsid region (Fig. 4). The results of SimPlot were equivalent to those of RIP; graphics from the RIP analysis are presented in this report. Collectively, these results suggested that Bo/CV521-OH/02/US was a recombinant generated by genetic recombination between Jena-like and Newbury-2-like BECVs.
![]() View larger version (34K): [in a new window] |
FIG. 4. Gene conversion analysis of Bo/CV521-OH/02/US compared with Bo/Jena/80/DE (black) and Bo/Newbury-2/76/UK (gray), using RIP (http://hivweb.lanl.gov/RIP/RIPsubmit.html). Program options were a window size of 100 and a threshold for statistical significance of 90%. The RdRp and capsid regions are indicated.
|
|
|
|---|
Based on comparison of the amino acid and nucleotide sequences of the complete capsid gene, the 21 BECV strains detected in fecal samples from two veal calf farms in Ohio were classified as Norovirus GIII/2 (Newbury-2-like, 15 strains) and NB-like BECV (6 strains). These viruses were phylogenetically distinct from human and porcine noroviruses and from human, porcine, and mink sapoviruses. Analyses of the complete capsid genes of the 15 BoNVs showed they shared lower, in comparison to GIII bovine noroviruses, amino acid and nucleotide identities with GI and GII human noroviruses, being more closely related to GI human noroviruses, which is consistent with a previous report (23). This finding also suggests that most bovine and human noroviruses may be genetically distinct. However, as a caveat to this finding, there are few sequence analyses of BoNVs and counterpart human noroviruses from developing countries, where closer contact between cattle and humans is expected. Bovine rotaviruses were thought to be distinctive from human rotaviruses until strains from cattle and children in India were compared and found to share common serotypes (16).
We analyzed the capsid genes of BECVs detected from two veal calf farms in 2002. The calves introduced into each farm were obtained from regional sale barns auctioning calves from four or five states, which suggests that these BECV strains have different geographic origins. Thus, the capsid gene diversity of BECVs detected on these farms is most likely a reflection of this diverse geographic sourcing. It was unclear whether calves were infected with genetically different BECVs before or after they were introduced to the veal farms. Two genetically different BoNV strains, such as (i) Bo/CV511-OH/02/US and Bo/CV511 M-OH/02/US and (ii) Bo/CV551-OH/02/US and Bo/CV551 M-OH/02/US were detected from the same sample, suggesting that calves were coinfected with the same genotype but with genetically different BECVs. Similarly, the RT-PCR results also showed some calves to be coinfected with different genogroups, GIII BoNV and the NB-like BECV. These coinfections with genetically different BECVs in a calf may result in natural genetic recombination between BoNVs, as suggested by our findings (Bo/CV521-OH/02/US), and/or possibly with NB-like BECVs. Genetic recombination may be a common mechanism for the generation of a diversity of genotypes or genogroups, although there is no evidence of recombination between the latter. Several natural recombinants between distinct human norovirus genotypes have been reported, as discussed below (17, 24, 37).
Phylogenetic trees generated from the complete capsid gene sequences of seven of the 15 BoNVs were compared to ones based on the RdRp gene sequences (Smiley and Saif, unpublished data). Except for Bo/CV521-OH/02/US, six BoNVs were classified as Norovirus GIII/2 irrespective of their RdRp or capsid gene sequences. Bo/CV521-OH/02/US was classified as Norovirus GIII/2 based on the whole capsid gene, but Smiley et al. (35) reported Bo/CV521-OH/02/US as Norovirus GIII/I, which was a Jena-like BoNV based on the RdRp gene sequence. This divergent result may be from coinfection by genetically different BECVs in a calf or from genetic recombination of BECVs. To exclude the possibility of amplification of different BECV due to contamination of genetically different BECVs in the same fecal sample, a contiguous portion of the RdRp and capsid genes of Bo/CV521-OH/02/US was amplified with primer set RC35/RC1932. Following analyses of the gene sequence and gene conversion, Bo/CV521-OH/02/US was found to likely be a recombinant, with Norovirus GIII/1-like origin for the RdRp gene and Norovirus GIII/2-like origin for the capsid gene. Natural genetic recombination between various genotypes or clusters of GI and GII human noroviruses has been detected for the WUG1 strain (19), the Arg320 strain (17), Hu/Snow Mountain/76/US (24), and Hu/Wortley/90/UK (37). However, to our knowledge, naturally occurring recombinants for BoNVs have not been previously reported or documented.
The putative region of recombination in human norovirus recombinants has been suggested to be the ORF1-ORF2 junction region (17, 19). For Bo/CV521-OH/02/US, we also suspect that recombination occurred in the ORF1-ORF2 junction (Fig. 4). This finding suggests that such recombination makes it difficult to classify noroviruses based on the C-terminal region of the RdRp or on complete and partial capsid gene sequences, but our subsequent analysis (Table 3) validated this approach. The capsid genes of noroviruses are composed of two domains, S and P (32). The S domain, which is located in the N terminus of the capsid gene, is the most conserved region. The P2 subdomain (aa 279 to 405 of the capsid protein) includes the hypervariable region and is the least conserved region. These patterns of sequence conservation were also consistent for the 15 BoNVs. The S domains of 15 BoNVs shared the highest amino acid identity with Hu/NLV/Norwalk/68/US, in contrast to the P1A, P2, and P1B subdomains (Table 3). The P1B subdomain of GIII/2 BoNVs was more conserved within a genotype than the P1A subdomain. Comparison of phylogenetic trees determined for the complete capsid gene or each domain of BoNVs and human noroviruses showed no significant differences in their placement, since only bootstrap values of greater than 70% were considered significant (Fig. 3). These results suggest that both the C-terminal region of the RdRp gene (ORF1) and the S domain of the capsid gene (ORF2) junction are reliable regions for genotyping noroviruses and discovering a recombinant.
The conserved motifs GXXPGIGKT in the 2C nucleoside triphosphatase, GDGC in the 3C protease, and GLPSG and YGDD in the RdRp regions of caliciviruses were known (4, 26). We identified the one previously described motif, LAGNA (S domain) (33), and the three new motifs, IDPWI (S domain), LYTP (S domain), and QNGR (P1A domain), in the capsid genes of noroviruses, including the 15 GIII BoNVs sequenced in this study. The motifs LAGNA and LYTP of all noroviruses were included in the globular domain as predicted with SMART/Pfam domain and defined by Linding et al. in GlobPlot version 2.1 (http://globplot.embl.de/) (22). However, the motifs IDPWI and QNGR of GII noroviruses were a potential globular domain, but these motifs of GIII BoNVs were not included in the globular domain. The IDPWI and QNGR motifs of some GI noroviruses, such as Hu/NLV/Norwalk/68/US and Hu/NLV/Desert Shield395/90/SA, were located in a globular domain. The role of each motif needs to be clarified for the noroviruses.
In summary, the results of this study demonstrate genetic diversity in the capsid genes of BECVs circulating in Ohio veal calves, provide new data for coinfections with distinct BECV genotypes or genogroups, and describe what is to our knowledge the first natural BoNV recombinant, analogous to the previously reported human norovirus recombinants.
This work was supported by grant R01 AI49716 from the National Institute of Allergy and Infectious Diseases, National Institutes of Health. M. G. Han was partially supported by the Postdoctoral Fellowship Program of the Korea Science and Engineering Foundation.
Present address: University of Kentucky, Division of Laboratory Animal Resources, Chandler Medical Center, Lexington, KY 40536. ![]()
|
|
|---|
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»