Bovine group A rotavirus (BoRV-A) is an important pathogen that causes neonatal calf diarrhea ([@bib30], [@bib13], [@bib9], [@bib11], [@bib1]). Viruses classified as group A rotavirus (RV-A), including BoRV-A, possess two independent neutralization antigens, namely, VP4 and VP7, on the outer capsid, that specify the P type (for protease-sensitive protein) and G type (for glycoprotein), respectively ([@bib10], [@bib32]). To date, six P types (P6\[1\], P7\[5\], P8\[11\], P11\[14\], P\[17\], and P\[21\]) and 8 G types (G1, G3, G5, G6, G7, G8, G10 and G15) have been reported among BoRV-A ([@bib37], [@bib3], [@bib4], [@bib29], [@bib14], [@bib15], [@bib44], [@bib1]). Types G2 and G11 in BoRV-A have been detected using PCR-derived cDNA probes ([@bib27]). Reports on the distribution patterns of these types reveal that P7\[5\], P8\[11\], G6, and G10 are the most common ([@bib13], [@bib11], [@bib24], [@bib40], [@bib1]). G8 ranks third among the most common G type in BoRV-A, and several G8 strains reported to date have been detected and/or isolated sporadically, or their prevalence has been reported ([@bib46], [@bib11], [@bib14], [@bib40], [@bib12]).

Bovine group B and C rotaviruses cause diarrhea in adult cattle ([@bib42], [@bib5], [@bib48], [@bib49], [@bib25], [@bib38]); in the case of BoRV-A, however, there were only few reports associating it with diarrhea in adult cattle ([@bib46]). There were several reports of adult diarrhea caused by RV-A in humans ([@bib26], [@bib31], [@bib8], [@bib23]). We discovered the prevalence of adult cattle diarrhea in Tochigi Prefecture, Japan, in 2004, and detected BoRV-As from these diarrheic specimens. The P and G types in these diarrheic specimens could not be identified using RT-PCR. The purpose of this study was to genetically characterize the detected BoRV-A by the partial sequence analysis of the genome segments encoding the proteins VP4 and VP7.

The prevalence of adult cattle diarrhea was noted in a dairy farm in Tochigi Prefecture, Japan, in October 2004. This dairy farm raised approximately 90 dairy cattle, and approximately 8 adult cattle suffered from watery diarrhea for nearly 9 days. In contrast, no calves raised on this farm suffered from diarrhea over this period. Unfortunately, we did not record the incidence of diarrhea in adult cattle and calves raised on this farm after the observation period. The diarrheic specimens were obtained from four adult cattle aged 3, 4, 5, and 8 years, and RV-A was detected from these specimens by using the commercial kit Rotascreen (Denka Seiken, Tokyo, Japan) according to the manufacturer\'s instructions. BoRV-A was attempted to isolate from the diarrheic specimens in MA-104 cells in the presence of trypsin as described previously ([@bib37]). Indirect immunofluorescence assay using anti-BoRV-A hyperimmune guinea pig serum, which has been described previously, was used to detect BoRV-As in MA-104 cells inoculated with the diarrheic specimens, culture supernatants, and infected cell culture lysates ([@bib33]). Viral genomic double-stranded RNA (dsRNA) was extracted from the diarrheic specimens by using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer\'s instructions. The extracted dsRNA was then analyzed by polyacrylamide gel electrophoresis (PAGE) and silver staining, as described previously ([@bib13]). The dsRNA was also used for first-stranded cDNA synthesis by MuLV-reverse transcriptase (Applied Biosystems, Foster City, CA, USA), and the cDNA was amplified by PCR using the ReadyMix Taq PCR Reaction Mix (Sigma--Aldrich, Saint Louis, MO, USA) as described previously ([@bib13]). In order to perform P typing, Bov4Com5, Bov4Com3, P1, P5, P10, con2, con3, pNCDV, pUK, and pB223 primers were synthesized as described by [@bib28] and [@bib19]. To perform G typing, sBeg9, End9(UK), DT6, HT8, and ET10 primers were synthesized as described by [@bib17], [@bib18], [@bib20]. Subsequently, the genome segments encoding the VP7s of the detected BoRV-As were amplified by PCR, using the primers VP7-MS (sense, 5′-CRGARYTAGATATGTCAGAA-3′, 491-510) and VP7-MA (antisense, 5′-AAYGTTATGTCCATYGGATT-3′, 563--544). Among the four diarrheic specimens obtained, the PCR-amplified products of the specimen termed Tak2 were selected and sequenced. The products were ligated into the pMOS*Blue* vector (Amersham Biosciences, Piscataway, NJ, USA), and both strands were sequenced using the BigDye Terminator Cycle Sequencing Kit (Applied Biosystems) on an ABI PRISM 3100 DNA analyzer (Applied Biosystems). The nucleotide (nt) and deduced amino acid (aa) sequences of the genome segments encoding VP4 and VP7 of Tak2 were compared with the corresponding RV-A sequences, which have been published previously. Phylogenetic trees were constructed by the neighbor-joining method ([@bib45]) using Clustal W ([@bib47]), and the bootstrap probabilities of each node were calculated using 1000 replicates. The final tree was obtained using the TreeView program, version 1.6.6 ([@bib41]). The other causative agents of adult cattle diarrhea, such as bovine group B and C rotaviruses, bovine coronavirus, bovine viral diarrhea virus, *Clostridium* spp., *Salmonella* spp. and *Coccidium* spp. were examined as described previously ([@bib6], [@bib50], [@bib48], [@bib49], [@bib2], [@bib7], [@bib38]).

RV-A was detected from all the collected diarrheic specimens by the commercial kit and indirect immunofluorescence assay, while the other causative agents of adult cattle diarrhea were not detected from these specimens. However, marked fluorescence and cytopathic effect were not observed in the MA-104 cells inoculated with the culture supernatants and infected cell culture lysates of the second and third passages (data not shown); thus, BoRV-A could not be isolated from these specimens. In addition, the dsRNAs of all diarrheic specimens could not be detected successfully using PAGE (data not shown).

In P typing, no amplified signal was observed in the second PCR using Bov4Com5, P1, P5, and P10 primers or in the first and second PCR using con2, con3, pNCDV, pUK, and pB223 primers, although weak signal was observed in the first PCR using Bov4Com5 and Bov4Com3 primers (data not shown). For the genome segment encoding VP4, an 817-bp fragment in the partial sequence corresponding to its VP5\* region was determined. The sequence reported here was deposited in the DDBJ/EMBL/GenBank databases (accession no. [AB259664](AB259664)). When the partial VP4 sequence of Tak2 was compared with those of the currently identified 24 P genotypes, it was found that these sequences share low identities ranging between 61.6% and 67.6% at the nt level and 58.0% and 71.4% at the aa level with the established 24 P genotypes ([Table 1](#tbl1){ref-type="table"} ). On the other hand, for the P\[22\] and P\[23\] types, these sequences could not be compared with those of Tak2 because only the sequences corresponding to their VP8\* regions have been published ([@bib35], [@bib39]). Phylogenetic analysis revealed that Tak2 formed a branch separate from the 24 established P genotypes ([Fig. 1](#fig1){ref-type="fig"} ).Table 1Comparison of the nt and aa sequence identities of the genome segment encoding VP5\* of Tak2 with strains representing the established twenty-four P typesStrainSpecies of originP genotypeP serotypeTak2Accession no.nt identity (%)aa identity (%)NCDV-LincolnBovine1665.767.8[AB119636](AB119636)SA11Simian25B66.769.2[D16346](D16346)RRVSimian35B65.770.0[M18736](M18736)L26Human41B65.367.8[M58292](M58292)UKBovine5765.268.9[M22306](M22306)ST3Human62A63.366.7[L33895](L33895)OSUPorcine7966.468.9[X13190](X13190)WaHuman81A65.765.9[L34161](L34161)AU-1Human9364.467.0[D10970](D10970)69MHuman10466.071.1[M60600](M60600)B223Bovine11861.658.0[D13394](D13394)H2Equine12467.671.4[D13397](D13397)MDR-13Porcine13Not tested65.265.9[L07886](L07886)Sun9Bovine141166.767.0[AB158430](AB158430)Lp14Ovine15Not tested64.867.8[L11599](L11599)EbMurine161064.068.5[L18992](L18992)993-83Bovine17Not tested63.263.4[D16352](D16352)L338Equine18Not tested65.268.5[D13399](D13399)4FPorcine191266.568.1[L10359](L10359)EHPMurine201363.166.7[U08424](U08424)Hg18Bovine21Not tested66.468.9[AF237665](AF237665)TUCHSimian24Not tested65.967.0[AY596189](AY596189)Dhaka6Human25Not tested64.465.8[AY773004](AY773004)134/04-15Porcine26Not tested65.064.8[DQ061053](DQ061053)Fig. 1The phylogenetic tree constructed by the neighbor-joining method, using the nucleotide sequences of the genome segments encoding VP5\* of the reference strains and Tak2, which is indicated in boldface. The accession numbers of the reference strains are listed in [Table 1](#tbl1){ref-type="table"}. The percentage bootstrap values are shown at the branch nodes. The scale bar represents 10% nucleotide difference. P genotypes are indicated after the strain names.

In G typing, no amplified signal was observed in the first and second PCR using sBeg9, End9(UK), DT6, HT8, and ET10 primers (data not shown). Therefore, two primers were synthesized according to the nt sequences of the genome segments encoding the VP7 of previously published BoRV-As, and amplification of the genome segment encoding the VP7 of Tak2 was attempted by PCR. Among the first PCR performed using sBeg9 and End9(UK) and the second PCR using sBeg9 and VP7-MA as well as VP7-MS and End9(UK), approximately 600 bp of the product was detected only in the second PCR performed using sBeg9 and VP7-MA (data not shown). Unfortunately, no amplified signal was observed in the second PCR using VP7-MS and End9(UK) (data not shown). For the genome segment encoding VP7, 522 bp of the partial sequence corresponding to approximately half of the amino terminal portion was determined. The sequence reported here was deposited in the DDBJ/EMBL/GenBank databases (accession no. [AB259665](AB259665)). The partial VP7 sequence of Tak2 was compared with those of strains representing the 15 previously established G types ([Table 2](#tbl2){ref-type="table"} ). The percentages of nt and aa identities with all the established 15 G types were low, ranging from 57.8% to 73.5% and 61.2% to 70.9%, respectively. Phylogenetic analysis revealed that Tak2 formed a branch separate from all the 15 established G types ([Fig. 2](#fig2){ref-type="fig"} ).Table 2Comparison of the nt and aa sequence identities of the genome segment encoding the half of the amino terminal side of VP7 of Tak2 with strains representing the established fifteen G typesStrainSpecies of originG sero/genotypeTak2Accession no.nt identity (%)aa identity (%)WaHuman171.170.3[K02033](K02033)S2Human269.764.8[M11164](M11164)HO-5Equine373.070.9[AB046464](AB046464)HOCHIHuman467.863.0[AB012078](AB012078)OSUPorcine571.667.9[X04613](X04613)NCDV-LincolnBovine672.468.5[M12394](M12394)PO-13Pigeon757.867.9[D82979](D82979)Sun9Bovine872.465.5[AB158431](AB158431)116EHuman969.366.7[L14072](L14072)KK3Bovine1073.566.7[D01056](D01056)YMPorcine1169.767.3[M23194](M23194)L26Human1267.864.2[M58290](M58290)L338Equine1371.767.3[D13549](D13549)CH3Equine1469.065.5[D25229](D25229)Hg18Bovine1571.761.2[AF237666](AF237666)Fig. 2The phylogenetic tree constructed by the neighbor-joining method, using the nucleotide sequences of the genome segments encoding half of the amino terminal portion of VP7 of reference strains and Tak2, which is indicated in boldface. The accession numbers of the reference strains are listed in [Table 2](#tbl2){ref-type="table"}. The percentage bootstrap values are shown at the branch nodes. The scale bar represents 10% nucleotide difference. The G types are indicated after the strain names.

RV-As are classified within the same P genotype if the relevant aa sequences share \>89% identity ([@bib16]). Considering this, our results suggested that BoRV-A Tak2 might belong to a novel P genotype, although only partial sequence of the genome segment encoding VP4 of Tak2 was determined, and the sequence could not be compared with those of the P\[22\] and P\[23\] types. However, phylogenetic analysis of the aa sequences of the VP8\* of 26 P genotypes in RV-As revealed that the P\[22\] and P\[23\] types were closely related to the P\[13\] and P\[24\] types, respectively ([@bib36]). In contrast, Tak2 was observed to share a common origin with the P\[11\] and P\[17\] types, which appeared to have a distinct genetic origin as compared to the VP4s belonging to the other P genotypes in previous study ([@bib44]). These results suggested that Tak2 could be a novel P genotype. Despite numerous efforts, the amplification of the PCR product, including the complete genome segment encoding VP4, of all four diarrheic specimens contained Tak2 and the determination of those sequences have been unsuccessful to date. The present study suggests that Tak2 is a BoRV-A with a novel P genotype specificity, which has not been reported among BoRV-A thus far.

RV-As with different G types share \<85% aa sequence conservation in the VP7, while strains sharing \>90% aa identity are considered to belong to the same G type ([@bib21], [@bib22], [@bib32]). Based on sequence and phylogenetic analysis of the partial sequence of the VP7 gene, strain Tak2 could not be assigned to any of the previously established G types. Accordingly, although the full-length VP7 gene remains to be determined, the data suggest that strain Tak2 may be a BoRV-A with a novel G type specificity, never reported among animal RV-A so far.

The feces of the animals were screened against a broad panel of viruses, bacteria and parasites, but the presence of other undetected pathogens, such as caliciviruses, parvoviruses, *Escherichia coli*, *Yersinia* spp., *Campylobacter* spp., or *Giardia* spp. was not investigated. Accordingly, it is not possible to associate the outbreak of diarrhea in adult cattle with the novel RV-A strain Tak2. However, BoRV-A in adult cattle with diarrhea have been already detected ([@bib46]). The factors responsible for symptomatic infections by BoRV-A in adult cattle are unknown. Immunity raised against BoRV-A carrying unusual P and G type specificities, such as Tak2, is estimated to be low because P7\[5\], P8\[11\], G6, and G10 are the most common P and G types in BoRV-A, according to previous studies ([@bib13], [@bib11], [@bib24], [@bib40], [@bib1]). Accordingly, the lack of homotypic immunity (specific for Tak2-like RV-As) may be one of the factors that have helped the virus to infect and cause disease in adult animals. The isolation of BoRV-A from the specimens of adult cattle diarrhea in this study needs to be repeated, and serological surveillance on these isolates in cattle is necessary in the future. However, the fact that no calf raised in this farm was affected with diarrhea is inconsistent with the above hypothesis. These questions warrant further studies to evaluate the pathogenic attitude of the strain Tak2.

Since 2000, RV-As carrying new P and G type specificities have been detected and/or isolated in quick succession, and almost all of those RV-As were of animal origin ([@bib44], [@bib34], [@bib35], [@bib36], [@bib39], [@bib43]). The discovery of BoRV-A Tak2, which may have a unique combination of novel VP7 and VP4 types, is an example of the enormous diversity among the circulating RV-A strains. Complete genome sequencing of all 11 segments may help to understand the origin of the unusual BoRV-A strain Tak2. The evolution of RV-As is the result of interspecies transmissions, reassortment of genome segments between diverse strains, genome rearrangements and accumulation of point mutations. The recent discovery of novel P and G types has important implications in future Human RV-A vaccine strategies. Since domestic animals, including cattle, may serve as an important reservoir of novel RV-A strains, epidemiological RV-A surveillance needs to be performed regularly.
