Previous Article | Next Article ![]()
Journal of Clinical Microbiology, December 2003, p. 5665-5675, Vol. 41, No. 12
0095-1137/03/$08.00+0 DOI: 10.1128/JCM.41.12.5665-5675.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Animal Health and Wellbeing, Faculty of Veterinary Medicine of Bari, Bari, Italy,1 Bioprocess and Bioanalytical Research, Merck & Co., Inc., West Point, New York,2 Department of Hygiene and Microbiology, University of Palermo, Palermo, Italy,3 Department of Medical Genetics, I.R.C.C.S. "Saverio De Bellis," Castellana Grotte, Italy,4 Institute "Zooprofilattico Sperimentale di Puglia e Basilicata," Sezione di Foggia, Italy5
Received 23 May 2003/ Returned for modification 5 August 2003/ Accepted 10 September 2003
|
|
|---|
|
|
|---|
Serotype designations are based on independent neutralization determinants on the two outer capsid proteins, VP7 (G serotype) and VP4 (P serotype) (19). On the basis of antigenic and genetic characterization, 15 G serotypes have been identified (19, 36). VP4-specific antisera or monoclonal antibodies raised to VP4 have identified 14 P serotypes, while sequence analysis of VP4 or the VP8* trypsin cleavage product of VP4 have recognized 23 P genotypes (30, 36, 42; Liprandi et al., unpublished data). With a few exceptions (10), strains sharing more than 89% amino acid identity are considered to belong to the same P genotype. Following the guidelines of this typing system, the P serotype is indicated with a number immediately after the letter P, while the genotype is represented by a number in square brackets (19, 30, 36).
One nonstructural protein of relevance is NSP4, which has been studied extensively because of its role in viral morphogenesis and its enterotoxic activity (3, 19). The NSP4 gene of group A rotaviruses may be genetically classified into at least five genogroups, KUN (A)-, Wa (B)-, Au-1 (C)-, EW (D)-, and avian (E)-like (12, 29, 31, 47). Within NSP4 genotypes A and B, rotavirus strains isolated from rabbits, horses, cows, and pigs generally cluster according to the species of origin, suggesting a constant pattern of evolution within species (12). In addition, the 11th dsRNA segment of rotaviruses encodes a phosphoprotein with kinase activity, NSP5, and, via an out-of-phase open reading frame, a smaller protein, NSP6; the NSP5/6 gene has been demonstrated to be suitable to trace the origin of rotaviruses (19, 37).
Rotavirus strains belonging to serotypes G6, G8, and G10, in association with VP4 types P6[1], P7[5], and P8[11], are commonly found in cattle, though strains belonging to G1, G2, G3, and G11 have been detected sporadically (18, 36, 57). An unusual P[17]- and G7-like bovine rotavirus has also been isolated from a calf, presumably the result of interspecies transmission from birds to cows (55). In addition, bovine rotaviruses with novel P[21] and G15 types have recently been identified in India (53). Bovine strains usually segregate in the NSP4 B genogroup, KUN-like (12).
In the present study, we report the isolation and molecular characterization of the VP8*, VP7, VP6, NSP4, and NSP5/6 of an unusual rotavirus strain, 10733, from a 7-day-old buffalo calf affected with severe gastroenteritis in southern Italy. Our results show that buffalo rotavirus strain 10733 possesses typical bovine-like VP6 (SGI), VP7 (G6), NSP4 (KUN-like), and NSP5/6 gene products. However, we report for the first time the identification of a rotavirus strain with an unusually high amino acid identity (96%) to the VP8* trypsin cleavage fragment of VP4, predictive of the P type, of rhesus rotavirus (RRV).
|
|
|---|
Electropherotype determination. Viral dsRNA was extracted from infected cells showing 50% cytopathic effect. Cryolysates were extracted with Vertrel XF (Du Pont, Wilmington, Del.) as previously described (44) and clarified twice by centrifugation for 30 min at 4,800 x g. Subsequently, the supernatants were centrifuged for 8 h at 85,000 x g on a 30% sucrose layer in an SW28 rotor (Beckman Coulter, Inc., Fullerton, Calif.). After digestion of the viral pellet with 1% sodium dodecyl sulfate and 1 µg of proteinase K (Sigma-Aldrich) per ml, viral dsRNA was extracted with phenol-chloroform, precipitated with ethanol and 0.3 M sodium acetate, and suspended in TE (Tris-EDTA, pH 8.0). The dsRNA e-type was visualized in a TAE (Tris-EDTA-acetate buffer, pH 8.0) gel with 0.5% agarose for resolving and 1.5% agarose for stacking, after electrophoresis at 15 V for 18 h and ethidium bromide staining.
RNA extraction and PCR amplification of the VP7, VP4, VP6, and NSP4. Viral dsRNA from both the fecal specimen and the 3rd viral passage on MA-104 cells was extracted by adsorption on cellulose CF11 as described previously (59). Viral dsRNA was denatured in dimethyl sulfoxide (Sigma-Aldrich) at 97°C for 5 min. Reverse transcription of dsRNA was carried out with murine leukemia virus reverse transcriptase (Applied Biosystems Italy), while PCR amplification was carried out with AmpliTaq Gold DNA polymerase (Applied Biosystems Italy).
The full-length VP7 gene (1,062 bp) was reverse transcribed and amplified with primer Beg9 (24) and primer End9deg, degenerated on the basis of the 3' end of the gene of 13 VP7 alleles (GGTCACATCDWMCARYTCTAAYYHM). The VP8* subunit of the VP4, the connecting peptide, and the N terminus of the VP5* subunit of VP4 (876 bp) were reverse transcribed and amplified with primer pair Con2-Con3 (20). The nearly full-length gene (725 bp) encoding the NSP4 protein was amplified with primer pair 10Beg16-10End722 (39), while the NSP5/NSP6 full-length gene (667 bp) was amplified with a primer pair described by Krishna Mohan and Atreya (37). The VP6 genogroup, predictive of the VP6 subgroup specificity, was determined by amplification of a 379-bp fragment, spanning from amino acids 241 to 367 of the VP6, with primer pair VP6F-VP6R (32).
Sequence analysis. For direct sequencing of PCR amplicons, three distinct amplicons were analyzed and a consensus sequence was determined. Moreover, the VP7- and VP8*-encoding genes of buffalo strain 10733 were cloned into the vector pCRT7/NT-TOPO (Invitrogen BV, Groningen, The Netherlands), and the sequence was determined with three plasmid clones. Sequences were assembled and analyzed with the Bioedited software package (Department of Microbiology, North Carolina State University) (28).
Phylogenetic and molecular evolutionary analyses were conducted with MEGA version 2.1 (Arizona State University) (38). Phylogenetic trees based on the VP7, VP8*, NSP4, and NSP5/6 gene products were elaborated with both parsimony and distance methods, supplying statistical support with bootstrapping over 100 replicates.
Nucleotide sequence accession numbers. GenBank accession numbers AY281360, AY281359, AY293829, and AY293830 were assigned to VP7, VP8*, NSP4, and NSP5, respectively.
|
|
|---|
Sequence analyses. To molecularly analyze buffalo rotavirus strain 10733, the genes encoding the VP7, VP8* trypsin cleavage product of VP4, VP6 SG-specific fragment, NSP4, and NSP5/6 were determined. The basic structure of the VP7-encoding cognate gene of the Italian buffalo strain 10733 was 1,062 nucleotides long, with two in-phase open reading frames (ORFs) beginning at nucleotides 49 and 136, and a single TAG codon at nucleotide 1027, coding for a VP7 protein of 297 or 326 amino acids, respectively.
The complete deduced amino acid sequence of the VP7 gene from buffalo strain 10733 was determined and compared to those of reference rotavirus strains belonging to all known G serotypes (Table 1). The VP7 amino acid sequence of strain 10733 was 98.7 to 87.4% identical to those of rotavirus strains exhibiting G6 serotype specificity, and the highest amino acid identity (98.7 and 93.5%) was found to P3[9],G6 human rotavirus strains Hun4 (Hungary) and PA151 (Italy), respectively. The VP7 protein of strain 10733 had a potential N-linked glycosylation site located at amino acid 69 (Asn), which tends to be conserved among rotavirus strains with the exception of the G6 bovine rotavirus strains RF and NCDV (10, 13, 48). In addition, buffalo strain 10733 had a second potential glycosylation site at amino acid 238, like most bovine and the human P[9],G6 rotavirus strains (Fig. 1).
|
View this table: [in a new window] |
TABLE 1. Amino acid comparison of VP7 of buffalo 10733, bovine VMRI-29 and NCDV, and human PA169 rotavirus strains with those of well-established G serotypes and with a selection of G6 rotavirusesa
|
![]() View larger version (50K): [in a new window] |
FIG. 1. Deduced amino acid sequence of the VP7 protein of buffalo strain 10733, of a selection of G6 and G3 rotaviruses, and of strain Wa, P1A[8],G1. The VP7 antigenic regions A, B, C, and F are indicated. The glycosylation site NST (amino acids 69 to 71) is indicated by asterisks. The accession numbers of the VP7 sequences are listed in Table 1, footnote a. Abbreviations: Bo, bovine; Bu, buffalo; Go, goat; Hu, human; Po, porcine; Si, simian.
|
When compared to published VP8* trypsin cleavage product sequences, the VP8* of buffalo rotavirus strain 10733 showed the highest nucleotide sequence and amino acid sequence identity, 95.6% and 96.2%, respectively, with those of the P5B[3] simian RRV strain (data not shown and Table 2). Of note, the second best match was to a goat G3 rotavirus, strain GRV, recently identified in Japan (40), with about 83% identity at the nucleotide and 90.8% at the amino acid level. Since it has been established that only rotaviruses that exhibit a VP4 amino acid identity of 89% or greater belong to the same P serotype or subtype (22), strain 10733 can be assigned to the P[3] genotype. The few amino acid differences (less than 3.8%) between the VP8* of strain 10733 and RRV are predictive of a serotype P5B specificity, even if this should be confirmed serologically.
|
View this table: [in a new window] |
TABLE 2. Amino acid comparison of VP8* of buffalo 10733, caprine GRV, and simian RRV rotavirus strains with well-established P genotypesa
|
![]() View larger version (71K): [in a new window] |
FIG. 2. Deduced amino acid sequence of the VP8* trypsin cleavage product of the VP4 protein of Italian buffalo strain 10733 and of a selection of strains representing the remaining P genotypes. The highly conserved cysteine (), prolines ( ), and arginines ( ) are indicated. For optimal alignment, gaps were introduced in the sequences. The first eight residues in the sequence of strain 10733 are inferred from the sequence of primer Con3 (nucleotides 10 to 32).
|
In addition, the 175-amino-acid-long NSP4 protein was deduced from genome segment 10 of buffalo rotavirus strain 10733 and compared with bovine strains of the NSP4 A (KUN-like) genotype and other representative NSP4 sequences of genogroups A, B, C, D, and E (Fig. 3). The NSP4 of strain 10733 showed the highest amino acid identity (97.1%) to the unusual G10 human strain A64 and to the bovine strain UK, P7[5],G6. Common structural features of the NSP4 protein, such as N-linked glycosylation sites at amino acids 8 and 18, the cysteine residues at amino acids 63 and 71, and tyrosine-131 (a key residue for enterotoxic activity) were conserved. Also, the sequence of the enterotoxic peptide was highly conserved (19). There were scattered amino acid variations throughout the protein with respect to the KUN-like strains. In the hypervariable region (amino acids 135 to 141) associated with altered virulence in mice (12), there were at least two amino acid changes for all the other NSP4s, except for human strain E210, which differs in a single residue, Iso-135 to Leu.
![]() View larger version (56K): [in a new window] |
FIG. 3. Deduced amino acid sequence of the NSP4 protein of Italian buffalo strain 10733 and of a selection of strains representing the five distinct NSP4 genotypes (12, 47). The potential glycosylation sites are indicated by asterisks. The highly conserved cysteine () and methionine ( ) residues and the putative VP4 and double-layered particle (DLP)-binding regions are shown. The enterotoxic peptide of simian rotavirus SA-11 strain at amino acids 114 to 135 and tyrosine-131 ( ), considered a key residue for enterotoxigenic activity, are shown (19). Also, the hypervariable region (amino acids 135 to 141) associated with altered virulence in mouse is shown (12). Abbreviations: av, avian; bo, bovine; bu, buffalo; ca, canine; eq, equine; fe, feline; go, goat; hu, human; la, lapine; mu, murine; po, porcine; si, simian.
|
Phylogenetic analyses of VP7, VP8* trypsin cleavage product, NSP4, and NSP5/6 of buffalo rotavirus strain 10733. Among rotavirus strains with serotype G6 specificity, four distinct VP7 branches or groups were observed (Fig. 4A). One group comprised most of the P6[1] and P7[5] serotype G6 bovine rotaviruses from various parts of the world, including the reference strains RF (France), UK (Britain), and NCDV and WC3 (United States). A second group included exclusively P[14],G6 human strains PA169 (Italy), MG6.01, ASG6.02, and MG6 (Australia), and the P[14] caprine strain CAP455 (South Africa), while the third group was formed by only P8[11],G6 bovine strains C8336, MC27, and VMRI-29 (United States), and KN-4 (Japan). The fourth group included buffalo rotavirus strain 10733 and the P3[9],G6 human rotavirus strains Hun4 (Hungary), PA151 (Italy), and Se584 (United States).
![]() ![]() ![]() ![]() View larger version (80K): [in a new window] |
FIG. 4. Phylogenetic trees of the VP7, VP8*, and NSP4 amino acid sequences of Italian buffalo strain 10733. (A) VP7 tree, displaying the relationships among a selection of serotype G6 and G3 animal and human rotavirus strains. The dendrogram is drawn to scale and rooted with the human strain Wa (P1A[8],G1). (B) VP8* tree, displaying the relationships among strains representative of all the VP4 genotypes recognized to date. The dendrogram is drawn to scale and rooted with the avian-like bovine strain Bo/993/83 (P[17],G7). (C) NSP4 tree, displaying the relationships among a selection of animal and human rotavirus strains representative of the five NSP4 genetic groups. The branches of the murine (genogroup D) and avian (genogroup E) strains are not to scale. The prototypes of the NSP4 genotypes A, B, and C (strains KUN, Wa, and Au-1, respectively) are boldfaced and underlined. Buffalo 10733 strain is shown in boldface type. (D) NSP5 tree, displaying the relationships among animal and human rotavirus strains. The tree is unrooted and is inferred from the nucleotide alignment. The P3[9] human rotaviruses (HRVs) are enclosed in a circle. Abbreviations: bo, bovine; bu, buffalo; ca, canine; eq, equine; fe, feline; go, goat; hu, human; la, lapine; mu, murine; po, porcine; si, simian.
|
|
|
|---|
In this study, we isolated a buffalo strain, 10733, displaying a long e-type and SGI specificity. Such a pattern is commonly described for animal rotaviruses (36). Sequence analysis of the VP7, VP8*, NSP4, and NSP5/6 genes revealed G6 and P[3] specificities, as well a bovine-like NSP5 gene and an NSP4 A (KUN-like) genogroup.
High amino acid and nucleotide sequence identities of the VP8* of strain 10733 to that of strain RRV, as well as the close phylogenetic relationship, strongly suggest that strains 10733 and RRV have analogous VP4s and may allow us to predict the specificity of buffalo strain as P5B[3]. While analyzing the sequences in the databases, we found that another strain, goat rotavirus GRV, shared high genetic relatedness with strains 10733 and RRV in the VP4 gene (40). The three strains were clustered together in the VP8* tree, constituting a lineage distinct from the P5A[3] rotavirus strains of canine, feline, and human origin (Fig. 4B). In addition, strain GRV displayed the same NSP4 C (AU-1-like) genogroup and G3 specificity as strain RRV (40), while strain 10733 exhibited an NSP4 A (KUN-like) genogroup, where all the bovine strains are grouped (12).
Analysis of VP7 provided contradictory data on the origin of strain 10733, as its VP7 exhibited the highest similarity (98.7%) to the Hungarian human strain Hun4, P3[9],G6. With a few exceptions, the G6 VP7 allele is one of the most common VP7 antigenic specificities of bovine rotaviruses (19), though it has also been detected in rotaviruses from humans (4, 14, 21). A significant intrahomotypic divergence within G6 rotaviruses from humans and cows has been identified, suggesting the existence of antigenically distinguishable G6 subtypes, as well as a correlation of subtypes with species of origin and P type (6, 14). By phylogenetic analysis of VP7, it was possible to observe that, while there is little or no difference between P[5],G6 and P[1],G6 rotaviruses identified in different parts of the world, the P[14],G6 and P8[11],G6 strains clearly form two distinct groups. It has been speculated previously (10) that the P[14] allele crossed the species barrier via sequential reassortments from rabbits to cows and from cows to humans. The existence of a G6 lineage restricted to human and caprine P[14] rotaviruses strongly supports this hypothesis. Interestingly, there was an additional cluster including the Italian strains 10733, P[3], and PA151, P[9], the American strain Se584, P[9], and the Hungarian strain Hun4, P[9] (4, 21, 26).
The linkage observed between G6 lineage and P type is indicative of multiple reassortment events occurring in the context of a constant pattern of linear evolution. Accordingly, introduction of bovine alleles into the human set of rotavirus gene alleles may have occurred on different occasions, and this would be confirmed by several pieces of evidence: (i) the existence of P3[9],G6 human strains genetically related to bovine strains (21, 26); (ii) the sporadic detection of unusual human strains bearing bovine-like specificities (1, 4, 5, 15); (iii) the high genetic relatedness between bovine and human strains in the NSP4 (12) and P3[9] human rotaviruses in the NSP5 gene (60); and (iv) genetic relatedness of strain 10733 to the unusual G10 human strain A64 (5) in NSP4 and to the unusual human strains Hun4, PA151, and Se584 (P3[9],G6) in the VP7 gene.
Phylogenetic analysis of multiple dsRNA segments has revealed that reassortment between homologous viruses (i.e., viruses of the same species) occurs frequently (33, 41, 58). While the exchange of genes between homologous strains is frequent but difficult to identify, heterologous exchanges are presumably much rarer but lead to the onset of "unusual" strains that may be readily identified because of uncommon G and/or P types. Cosegregation of G6 lineage and P type may be explained as preferential gene selection (49), as development of stable and advantageous conformational alterations in the outer viral capsid of some VP4-VP7 reassortants (51), or as ecological obstacles to free reassortment.
The detection of the rare P[3] RRV-like VP4 allele in ruminants raises questions on the origin and the possible source of introduction of this new VP4 type in the set of alleles of ruminant rotaviruses. The simian rhesus strain RRV was originally isolated from the feces of a 3- to 5-month-old rhesus monkey with diarrhea (56). Strain RRV has been classified both serologically and genetically as P5B[3],G3 (36), with an SGI VP6 and AU1-like NSP4 genotype (9). While the G3 specificity has been reported in a wide range of mammalian species (19), so far, the prototype strain RRV is the only rotavirus possessing the P5B[3] VP4 allele (30). Although it is important as a reference strain in laboratories and as the recipient virus in the human-rhesus vaccine (35, 44, 45), to our knowledge, viruses with RRV-like P type have never been detected, and the P5B[3] VP4 allele is not considered of epidemiological relevance in either humans or animals.
Thus, the finding that ruminant rotaviruses (strains 10733 and GRV) are highly similar in the VP4 allele to strain RRV has only two possible explanations: (i) the RRV-like VP4 type is an allele also present in rotaviruses circulating in ruminants but has not been detected previously because of its rare distribution and/or the typing assays used; and (ii) the RRV-like gene allele has been introduced recently into the genome of ruminant rotaviruses as the result of a reassortment event that occurred with a wild simian rotavirus or during the vaccine trials. For instance, horizontal transmission of a vaccine virus has been demonstrated both in vaccinees and in placebo recipients during vaccine trials in Venezuela (50). As the rhesus rotavirus vaccine has been administered to infants in several countries throughout the world since the mid-1980s (7, 19, 23, 36, 54), environmental contamination with rhesus strain RRV may have occurred several times and in different geographic locations in the last 20 years. The latter hypothesis is extremely intriguing, evocative of the potential biological risks coming from manipulation of live viruses that possess genetic and antigenic plasticity such as RNA segmented viruses. However, divergence within the VP8* sequence of the RRV-like P[3] strains (up to 16 to 17% at the nucleotide level and up to 9% at the amino acid level) seems to be consistent with analogous variation described within strains of the same P type (10, 25, 34, 42) and suggests that a VP4 RRV-like allele is circulating in ruminants, presumably at very low frequency. Furthermore, the chances of environmental contamination would be expected to be far greater in countries where the vaccine trials were made, but the surveillance systems active in those countries have not reported, to date, the existence of either human or animal RRV-like rotaviruses, and therefore, it is highly unlikely that reassortment between strain RRV and other rotaviruses may have occurred in the field.
The detection of a buffalo P5B[3],G6 rotavirus, while providing additional evidence for the genetic and antigenic diversity of group A rotaviruses, is also important from the perspective of understanding the basis of rotavirus host specificity. Although rotaviruses have low host specificity, heterologous infections usually have no clinical signs and are not efficiently propagated (8, 9, 19). Strain 10733 was isolated from a diarrheic buffalo calf, so the presence of this VP4 allele in the background of one bovine-like rotavirus is presumably still permissive for symptomatic infection. Interestingly, the P5B[3],G3 RRV strain is also the only nonlapine rotavirus able to determine productive infection in the rabbit model (8, 9).
Comparative sequence analysis of other genes of strain 10733, as well as experimental infection of buffalo and/or bovine calves, could provide important information on the molecular mechanisms of rotavirus host range restriction as well as help to elucidate the origin of this strain.
We thank Donato Narcisi, Carlo Armenise, Paola Fiorente, and Filomena Cariola for expert technical assistance. We are extremely grateful to Leland Eugene Carmichael for continued encouragement throughout the study.
|
|
|---|
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»