Introduction {#s0005}
============

Rotaviruses (RVs) are a major etiological agent of acute viral gastroenteritis in young animals and children worldwide. RVs belong to the *Reoviridae* family, possess a genome composed of 11 segments of double-stranded RNA (dsRNA) and are currently classified into eight groups (A--H) based on antigenic properties and sequence based classification of the inner viral capsid protein 6 (VP6) ([@bib4], [@bib15], [@bib34], [@bib49]). RVA is considered the most important RV group because of its high prevalence and pathogenicity in both humans and a variety of domestic and wild animals ([@bib15], [@bib30], [@bib33]). While human RVA and RVC have been described around the world, current reports indicate that human RVB strains have been described only in China ([@bib9], [@bib11], [@bib16], [@bib20]), India ([@bib22], [@bib25]), and Bangladesh ([@bib1], [@bib39], [@bib40], [@bib42]). While porcine RVB were first identified in the 1980s ([@bib6], [@bib47]) four other rotavirus groups, RVA, RVC, RVE, and RVH have also been described in pigs ([@bib8], [@bib21], [@bib53]). In addition to pigs, RVB strains have been also detected in cattle ([@bib3], [@bib7], [@bib17], [@bib52]), lambs ([@bib44]), and rats ([@bib14]). No genetic information exists for porcine RVE strains, only one paper has reported swine RVH ([@bib53]), and a limited number of RVC strains have been published ([@bib10], [@bib26], [@bib28], [@bib50]).

RVA have been well characterized and hundreds of complete genomes have been sequenced. In 2008, a sequenced-based classification system was proposed for RVA strains, and for each of the 11 gene segments, a nucleotide percent identity cut-off value was determined to discriminate between genotypes based on the phylogenetic analyses and pairwise sequence identity profiles ([@bib31]). Subsequently, a Rotavirus Classification Working Group (RCWG) was formed to set the RVA classification guidelines and maintain the proposed classification system ([@bib32]). To avoid misclassification of RVA strains not belonging to any of the established RVA genotypes, new sequences are submitted to the RCWG. The correct genotype is assigned, or a novel genotype is created. Currently, only RVA classification has been developed and maintained by the RCWG. Several other viruses such as astroviruses ([@bib37]), sapoviruses ([@bib45]), noroviruses ([@bib56]), hantavirus ([@bib27]), and papillomaviruses ([@bib12]) have been successfully classified using similar methodology.

Serological and molecular characterization of RVB strains is limited due to the difficulty of adapting RVB strains to cell culture ([@bib5], [@bib43]). Complete genome sequencing has been performed on several human RVB strains from Southeast Asia ([@bib1], [@bib23], [@bib54], [@bib55]). A limited number of RVB gene segments from a rat and bovine have also been sequenced ([@bib3], [@bib7], [@bib14], [@bib17], [@bib22], [@bib38], [@bib39], [@bib52]). Due to the increased ease of sequencing, Kuga and colleagues (2009) sequenced the VP7 of 38 porcine RVB strains and constructed phylogenetic trees and pairwise identity frequency graphs for G genotype classification purposes ([@bib24]). Based on their analyses, they proposed 5 genotypes which were further divided into 12 clusters, using 67% and 76% nucleotide cut-off values (66% and 79%, respectively, on the amino acid level) ([@bib24]).

In order to better understand the role of RV in porcine gastroenteritis cases, stool samples and small intestines from pigs of different ages were tested for RVA, RVB, and RVC as well as common pathogens associated with porcine diarrhea. The VP7 open reading frame (ORF) for 68 swine RVB samples from the United States was determined, revealing a substantial genetic diversity of porcine RVB strains. The genotype classification cut-off values as proposed by Kuga and colleagues ([@bib24]) were modified after including the novel sequence data generated in this study. Based on phylogenetic and molecular evolutionary analyses, a modified classification system is proposed currently describing 20 RVB VP7 genotypes based on the nucleotide cut-off value of 80%.

Results {#s0010}
=======

Small intestines or fecal samples from pigs with current or previous evidence of gastrointestinal disease (with or without diarrhea) are routinely examined by light microscopy for RV-like lesions and tested for RVA, RVC, and transmissible gastroenteritis coronavirus (TGEV) by RT-PCR at the University of Minnesota Veterinary Diagnostic Laboratory. Initially, one sample revealed RV-like lesions including villus atrophy, fusion, and necrosis by light microscopic but was negative for RVA and RVC by RT-PCR. Electron microscopy confirmed the presences of RV-like particles. The sample was found positive using an RT-PCR specially designed to detect the VP7 of RVB, and sequencing of the VP7 revealed a genetic identity of approximately 70% to bovine RVB strains, since porcine VP7 sequences were unavailable in GenBank at that time (August 2009).

Subsequently, 173 samples from pigs of all ages and residing in 119 different sites across the United States were retested for RVB by RT-PCR. Of the 173 samples, 156 samples (90.2%) were positive for RVA, RVB, and/or RVC while 17 samples (9.8%) were negative for RVs ([Fig. 2](#f0010){ref-type="fig"}). RVA and RVC were found in 110 samples (63.6%) and 100 samples (57.8%), respectively. RVB was found in 81 samples (46.8%). The majority (70 of 81) of RVB positive intestines were simultaneously positive for either RVA or RVC. Among all of the pathogens tested, only RVB was identified in the intestinal samples from nine pigs. These nine pig samples had typical RV lesions, such as those in [Fig. 1](#f0005){ref-type="fig"}, and were confirmed to contain RV particles by electron microscopy, and were negative for RVA and RVC by RT-PCR. Three of the 173 samples (1.7%) were positive for TGEV. Of these three samples, one sample was also found to be positive for RVA.Fig. 1Sections of jejunum from 2--3 day old pigs with typical lesions of acute RV enteritis by 200× magnification. Villi are characterized by edema, early epithelial necrosis, and sloughing to overt epithelial necrosis with denuded villus tips. (a) RVA infection (b) RVB infection (c) RVC infection (d) Normal villi.

In an attempt to estimate prevalence of RV detection by age of pig, the samples were divided by age groups as follows: \<21, 21--55, and \>55 days ( [Fig. 2](#f0010){ref-type="fig"}). The highest percentage of RVB positive samples was observed in the \>55 days age group (72.7%), followed by the 21--55 days age group (64.2%), and the \<21 days age group (18.6%). The highest percentage of RVB-only positives were seen in the \<21 days age group (12.9%), followed by the 21--55 days age group (2.5%). RVB only positive samples were absent in the \>55 days age group. The highest percentage of RVB strains found in a mixed infection with either RVA or RVC was in the \>55 days age group with RVA (27.3%). The highest percentage of a RVB mixed infection with RVA and RVC was present in the \>55 days age group (45.5%). Porcine enteric pathogens *S. enterica* spp., *Escherichia coli*, *Clostridium difficile* and *C. perfringens*, coccidia, *Lawsonia intracellularis* and weakly hemolytic *Brachyspira* species were detected in different age groups; however, the correlations between these pathogens and RV infections were outside the scope of this study.Fig. 2(a) Histogram of distribution of RVA, RVB and RVC RT-PCR positive porcine samples by 3 age groups: \<21 days, 21--55 days, and \>55 days. (b) Numeric values of positive samples by age group.

We attempted to determine the complete RVB VP7 ORF from all 81 positive samples; however, 13 of the 81 samples contained more than one RVB strain, which resulted in the inability to generate an accurate full ORF sequence of the VP7 gene. Therefore, the VP7 of only 68 RVB strains identified in 51 different sites was successfully determined. [Table 1](#t0005){ref-type="table"} lists each of the 68 RVB strain names, clinical signs of diarrhea or lack thereof, RVA, RVC, and TGEV results, the age of the host, G genotype, state of origin, site and GenBank accession number. Of the 68 sequenced positive RVB samples, 36 had clinical signs of diarrhea, whereas in six samples, the presence of diarrhea was unknown. Single RVB G genotypes were detected in 47 sites while two G genotypes were found to co-circulate in four sites (G8/G16, G12/G16, G12/G17, and G12/G18).Table 1Distribution of RVB sequences by diarrhea status, PCR results (RVA, RVC, and TGEV), age, G genotype, state, site, and accession number (*n*=68).**NCBI nameDiarrhea statusRVA resultRVC resultTGE resultAgeG genotypeStateSiteAccession number**RVB/Pig-wt/USA/IL09-1/2009/G18P\[X\]+−−−1--10 dayG18Illinois1[JQ043748](ncbi-n:JQ043748)RVB/Pig-wt/USA/IL09-2/2009/G18P\[X\]+−−−1--10 dayG18Illinois1[JQ043749](ncbi-n:JQ043749)RVB/Pig-wt/USA/IL09-3/2009/G18P\[X\]+−−−1--10 dayG18Illinois1[JQ043750](ncbi-n:JQ043750)RVB/Pig-wt/USA/IL09-4/2009/G18P\[X\]+−−−1--10 dayG18Illinois1[JQ043751](ncbi-n:JQ043751)RVB/Pig-wt/USA/IL09-5/2009/G18P\[X\]+−−−1--10 dayG18Illinois1[JQ043752](ncbi-n:JQ043752)RVB/Pig-wt/USA/MO09-32/2009/G16P\[X\]+−−−4 dayG16Missouri23[JQ043779](ncbi-n:JQ043779)RVB/Pig-wt/USA/MN09-42/2009/G20P\[X\]+−−−3 dayG20Minnesota31[JQ043789](ncbi-n:JQ043789)RVB/Pig-wt/USA/MN09-59/2009/G14P\[X\]+−−−3--10 dayG14Minnesota43[JQ043806](ncbi-n:JQ043806)RVB/Pig-wt/USA/IL09-25/2009/G16P\[X\]+−+−4 weeksG16Illinois16[JQ043772](ncbi-n:JQ043772)RVB/Pig-wt/USA/MN09-68/2009/G12P\[X\]+−+−4 weeksG12Minnesota21[JQ043777](ncbi-n:JQ043777)RVB/Pig-wt/USA/MN09-58/2009/G14P\[X\]+−+−3--10 dayG14Minnesota43[JQ043805](ncbi-n:JQ043805)RVB/Pig-wt/USA/MN09-6/2009/G8P\[X\]++−−6 weeksG8Minnesota2[JQ043753](ncbi-n:JQ043753)RVB/Pig-wt/USA/NC09-14/2009/G16P\[X\]++−−5 weeksG16North Carolina8[JQ043761](ncbi-n:JQ043761)RVB/Pig-wt/USA/MO09-17/2009/G16P\[X\]++−−3 weeksG16Missouri10[JQ043764](ncbi-n:JQ043764)RVB/Pig-wt/USA/MN09-35/2009/G6P\[X\]++−−4 weeksG6Minnesota26[JQ043782](ncbi-n:JQ043782)RVB/Pig-wt/USA/NC09-41/2009/G16P\[X\]++−−4 weeksG16North Carolina30[JQ043788](ncbi-n:JQ043788)RVB/Pig-wt/USA/MN09-43/2009/G8P\[X\]++−−15 weeksG8Minnesota32[JQ043790](ncbi-n:JQ043790)RVB/Pig-wt/USA/MO09-56/2009/G8P\[X\]++−−5.5 weeksG8Missouri42[JQ043803](ncbi-n:JQ043803)RVB/Pig-wt/USA/OH09-60/2009/G11P\[X\]++−−4 weeksG11Ohio44[JQ043807](ncbi-n:JQ043807)RVB/Pig-wt/USA/MN09-64/2009/G8P\[X\]++−−18 dayG8Minnesota48[JQ043811](ncbi-n:JQ043811)RVB/Pig-wt/USA/IA09-67/2009/G16P\[X\]++−−4 weeksG16Iowa50[JQ043814](ncbi-n:JQ043814)RVB/Pig-wt/USA/MN09-68/2009/G12P\[X\]++−−5 weeksG12Minnesota51[JQ043815](ncbi-n:JQ043815)RVB/Pig-wt/USA/MN09-20/2009/G11P\[X\]+++−24 dayG11Minnesota12[JQ043767](ncbi-n:JQ043767)RVB/Pig-wt/USA/MO09-21/2009/G12P\[X\]+++−4--5 weeksG12Missouri13[JQ043768](ncbi-n:JQ043768)RVB/Pig-wt/USA/MO09-22/2009/G16P\[X\]+++−4--5 weeksG16Missouri13[JQ043769](ncbi-n:JQ043769)RVB/Pig-wt/USA/IL09-23/2009/G20P\[X\]+++−4 weeksG20Illinois14[JQ043770](ncbi-n:JQ043770)RVB/Pig-wt/USA/MN09-24/2009/G12P\[X\]+++−4 weeksG12Minnesota15[JQ043771](ncbi-n:JQ043771)RVB/Pig-wt/USA/NE09-26/2009/G8P\[X\]+++−14--16 weeksG8Nebraska17[JQ043773](ncbi-n:JQ043773)RVB/Pig-wt/USA/MN09-27/2009/G12P\[X\]+++−6 weeksG12Minnesota18[JQ043774](ncbi-n:JQ043774)RVB/Pig-wt/USA/NC09-28/2009/G16P\[X\]+++−4 weeksG16North Carolina19[JQ043775](ncbi-n:JQ043775)RVB/Pig-wt/USA/TX09-46/2009/G12P\[X\]+++−10 weeksG12Texas35[JQ043793](ncbi-n:JQ043793)RVB/Pig-wt/USA/TX09-47/2009/G18P\[X\]+++−10 weeksG18Texas35[JQ043794](ncbi-n:JQ043794)RVB/Pig-wt/USA/MO09-57/2009/G16P\[X\]+++−5.5 weeksG16Missouri42[JQ043804](ncbi-n:JQ043804)RVB/Pig-wt/USA/NC09-61/2009/G16P\[X\]+++−4 weeksG16North Carolina45[JQ043808](ncbi-n:JQ043808)RVB/Pig-wt/USA/MN09-62/2009/G12P\[X\]+++−3 weeksG12Minnesota46[JQ043809](ncbi-n:JQ043809)RVB/Pig-wt/USA/MN09-65/2009/G8P\[X\]+++−18 dayG8Minnesota48[JQ043812](ncbi-n:JQ043812)RVB/Pig-wt/USA/MN09-7/2009/G8P\[X\]−+−−20 weeksG8Minnesota3[JQ043754](ncbi-n:JQ043754)RVB/Pig-wt/USA/CO09-12/2009/G8P\[X\]−+−−14 weeksG8Colorado6[JQ043759](ncbi-n:JQ043759)RVB/Pig-wt/USA/MO09-34/2009/G10P\[X\]−+−−3 weeksG10Missouri25[JQ043781](ncbi-n:JQ043781)RVB/Pig-wt/USA/OK09-40/2009/G16P\[X\]−+−−4 weeksG16Oklahoma29[JQ043787](ncbi-n:JQ043787)RVB/Pig-wt/USA/MN09-55/2009/G8P\[X\]--+----5 weeksG8Minnesota41[JQ043802](ncbi-n:JQ043802)RVB/Pig-wt/USA/MN09-63/2009/G18P\[X\]−+−−5 weeksG18Minnesota47[JQ043810](ncbi-n:JQ043810)RVB/Pig-wt/USA/MN09-8/2009/G8P\[X\]−++−20 weeksG8Minnesota3[JQ043755](ncbi-n:JQ043755)RVB/Pig-wt/USA/AR09-9/2009/G14P\[X\]−++−3 weeksG14Arkansas4[JQ043756](ncbi-n:JQ043756)RVB/Pig-wt/USA/PA09-10/2009/G12P\[X\]--++--6 weeksG12Pennsylvania5[JQ043757](ncbi-n:JQ043757)RVB/Pig-wt/USA/CO09-11/2009/G8P\[X\]--++--14 weeksG8Colorado6[JQ043758](ncbi-n:JQ043758)RVB/Pig-wt/USA/MO09-13/2009/G20P\[X\]--++--3 weeksG20Missouri7[JQ043760](ncbi-n:JQ043760)RVB/Pig-wt/USA/NC09-15/2009/G16P\[X\]--++--4 weeksG16North Carolina8[JQ043762](ncbi-n:JQ043762)RVB/Pig-wt/USA/IL09-16/2009/G20P\[X\]−++−4 weeksG20Illinois9[JQ043763](ncbi-n:JQ043763)RVB/Pig-wt/USA/OK09-18/2009/G16P\[X\]−++−5 weeksG16Oklahoma11[JQ043765](ncbi-n:JQ043765)RVB/Pig-wt/USA/OK09-19/2009/G16P\[X\]−++−5 weeksG16Oklahoma11[JQ043766](ncbi-n:JQ043766)RVB/Pig-wt/USA/NC09-29/2009/G12P\[X\]−++−35 dayG12North Carolina20[JQ043776](ncbi-n:JQ043776)RVB/Pig-wt/USA/OK09-38/2009/G16P\[X\]−++−4 weeksG16Oklahoma29[JQ043785](ncbi-n:JQ043785)RVB/Pig-wt/USA/OK09-39/2009/G16P\[X\]−++−4 weeksG16Oklahoma29[JQ043786](ncbi-n:JQ043786)RVB/Pig-wt/USA/KS09-44/2009/G16P\[X\]−++−28 dayG16Kansas33[JQ043791](ncbi-n:JQ043791)RVB/Pig-wt/USA/KS09-48/2009/G16P\[X\]−++−3 weeksG16Kansas36[JQ043795](ncbi-n:JQ043795)RVB/Pig-wt/USA/KS09-49/2009/G16P\[X\]−++−3 weeksG16Kansas36[JQ043796](ncbi-n:JQ043796)RVB/Pig-wt/USA/OK09-50/2009/G12P\[X\]−++−3 weeksG12Oklahoma37[JQ043797](ncbi-n:JQ043797)RVB/Pig-wt/USA/OK09-51/2009/G17P\[X\]−++−3 weeksG17Oklahoma37[JQ043798](ncbi-n:JQ043798)RVB/Pig-wt/USA/OK09-52/2009/G8P\[X\]−++−5 weeksG8Oklahoma38[JQ043799](ncbi-n:JQ043799)RVB/Pig-wt/USA/OK09-53/2009/G8P\[X\]−++−8 weeksG8Oklahoma39[JQ043800](ncbi-n:JQ043800)RVB/Pig-wt/USA/MN09-54/2009/G12P\[X\]−++−5 weeksG12Minnesota40[JQ043801](ncbi-n:JQ043801)RVB/Pig-wt/USA/OK09-33/2009/G16P\[X\]unknown−+−10 dayG16Oklahoma24[JQ043780](ncbi-n:JQ043780)RVB/Pig-wt/USA/MN09-37/2009/G20P\[X\]unknown+−−3 weeksG20Minnesota28[JQ043784](ncbi-n:JQ043784)RVB/Pig-wt/USA/SD09-45/2009/G14P\[X\]unknown+−−14 dayG14South Dakota34[JQ043792](ncbi-n:JQ043792)RVB/Pig-wt/USA/MN09-66/2009/G8P\[X\]unknown+−−18 dayG8Missouri49[JQ043813](ncbi-n:JQ043813)RVB/Pig-wt/USA/MO09-31/2009/G18P\[X\]unknown++−21 dayG18Missouri22[JQ043778](ncbi-n:JQ043778)RVB/Pig-wt/USA/MO09-36/2009/G8P\[X\]unknown++−3 weeksG8Missouri27[JQ043783](ncbi-n:JQ043783)

Using the 57 VP7 ORF nucleotide sequences of RVB available from GenBank, together with the 68 new VP7 sequences of RVB strains from this study, a phylogenetic tree ( [Fig. 3](#f0015){ref-type="fig"}) was constructed and a nucleotide sequence based pairwise identity frequency graph was generated ( [Fig. 4](#f0020){ref-type="fig"}). In the pairwise identity frequency graph, the typical pattern of clearly resolved peaks is apparent. The use of a 68% cut-off value resulted in 5 distinct clades in the phylogenetic tree ([Fig. 3](#f0015){ref-type="fig"}, left dashed line), with a limited overlap between the inter- (red) and intra-cluster (orange and green) identities ([Fig. 4](#f0020){ref-type="fig"}). Four of these clades contained RVB strains from only one host species (murine, human and two clusters of porcine RVB), while clade 3 contained both bovine and porcine RVB strains ([Fig. 3](#f0015){ref-type="fig"}). A second suitable cut-off value of 80% resulted in the identification of 20 G genotypes ([Fig. 3](#f0015){ref-type="fig"}, the right dashed line) in the phylogenetic tree, with a very limited overlap between inter- (red and orange) and intra-genotype (green) identities ([Fig. 4](#f0020){ref-type="fig"}). These 20 different genotypes contained RVB isolated from a single host species (G1: murine, G2: human, G3 and G5: bovine, and G4, G6--G20: porcine) ([Fig. 3](#f0015){ref-type="fig"}). [Table 2](#t0010){ref-type="table"} show the nucleotide and amino acid identity ranges among RVB strains belonging to the same G genotype. As could be expected the nucleotide identity ranges were above the determined 80% cut-off value for the majority of the genotypes. Only the range among RVB strains belonging to the G12 genotype extended slightly below 80%.Fig. 3Phylogenetic tree of RVB VP7 nucleotide sequences. Numbers left of the node characterize the bootstrap (1000 replicates) for each cluster. Bootstrap values less than 50% are not shown. The left dashed line represent the 68% cut-off value, resulting in 5 clades (or genotypes in the old classification of Kuga and colleagues), and the right dashed line represents the 80% nucleotide cut-off value, resulting in the 20 new G genotypes. The new and old clasification system are shown on the right side of the figures in brackets.Fig. 4Pairwise identity frequency graph using 125 available RVB VP7 ORF nucleotide sequences. Two suitable nucleotide cut-off values (68% and 80%) are depicted by vertical solid lines. (For interpretation of the references to color in this figure, the reader is referred to the web version of this article.)Table 2RVB VP7 nucleotide and amino acid pairwise frequencies comparison within genotypes.**G genotypeNucleotideAmino acid**G1nanaG291.0--99.995.5--100G392.5--95.392.0--99.0G4nanaG599.4--99.999.5--100G680.585.6G782.1--87.087.6--94.0G889.4--10093.0--100G9nanaG10nanaG1187.3--96.693.0--97.0G1278.3--99.985.1--100G13nanaG1481.2--10092.0--100G15nanaG1683.6--10087.1--100G17nanaG1882.3--10089.6--100G19nanaG2080.0--98.190.0--98.0

Alignment of the deduced amino acid VP7 sequences from the 20 genotypes is shown in [Fig. 5](#f0025){ref-type="fig"}. Six conserved cysteine residues (29, 84, 108, 137, 148, and 189) were present in the selected representative RVB strains, suggesting that these may be important to the folding, stability, and function of the VP7 protein. An additional cysteine residue was present at residue 13 of genotypes G3--G5 and G11--G20. RVB VP7 sequences contained 1--4 potential N-linked glycosylation sites: a conserved site at residues 47--49, and additional sites including residues 64--66, 107--109, and 242--245 for genotypes G11--G20. Putative signal cleavage sites were located at residues 15--16 (Ala and Gln, respectively) and were conserved except for G1 (IDIR) and G8. The VP7 amino acid sequences for the 125 RVB strains (grouped by genotype) used in the phylogenetic analysis are presented in [Fig. 6](#f0030){ref-type="fig"}. The data show that the genetic diversity (both conserved and highly variable amino acid positions) across the VP7 amino acid sequence was evenly distributed with an apparent slightly higher level of conservation in the region at the C-terminal.Fig. 5Multiple amino acid sequence alignment of selected RVB strains including the new and previously described VP7 genotypes. Matching amino acids are indicated with dots; gaps are shown by dashes; potential N-linked glycosylation sites are underlined; cysteine residues are boxed; insertion sites are indicated by addition symbols; ambiguous amino acids are represented by question mark, and putative signal cleavage sites are indicated by carets. Conserved residues within the 125 strains are shaded, while asterisks indicate highly divergent positions.Fig. 6Visual representation of the alignment for the 125 RVB amino acid sequence divided into the 20 proposed genotypes. The designated G genotype is listed on the left side. The consensus amino acid sequence is listed at the top of the chart where each color represents a specific amino acid. The histogram represents diversity within each amino acid position where peaks (dark green) represent conserved residues while valleys (red) represent divergent amino acid residues.

Discussion {#s0015}
==========

Kuga and colleagues ([@bib24]) and Collins and colleagues ([@bib10]) reported higher prevalence of RVA and RVC than RVB in pigs. However, our results also indicate that nearly half (46.8%; 81/173) of the tested samples contained RVB, which is higher than reported previously. Interestingly, the data reported in this study suggest that co-infection of RVB with RVA, and/or RVC can appear in pigs without concurrent clinical diarrhea. While porcine RVA or RVC infections can appear as both symptomatic and asymptomatic, it is important to remember that asymptomatic animals can still have remarkable gastrointestinal lesions of varying severity. Diarrhea is only one clinical sign associated with a RV infection. Less obvious clinical signs are associated with malabsorption of nutrients by villi with lesions, and may be manifested by loss of body weight, insufficient weight gain, lethargy, and/or mild dehydration. Our findings suggest that RVB can be a primary enteric pathogen in pigs. To that end, we found nine samples (eight of which yielded accurate gene sequences) positive by the RVB VP7 RT-PCR that were not only associated with diarrhea and had the presence of intestinal lesions typical of a RV infection, but also had RV particles identified by electron microscopy. These same samples were negative for other enteric viral or bacterial pathogens tested (*E. coli*, *C. difficile* and *C. perfringens*, coccidian, *L. intracellularis* and *Brachyspira species*, details of testing not shown), supporting the possibility that RVB can be an important cause of porcine enteritis. In addition, since all the samples were positive by the RVB VP7 RT-PCR designed specifically to detect only RVB, it is highly unlikely that the intestinal lesions and RV particles by electron microscopy were associated with other RV groups. However, we cannot absolutely rule out the possible presence of other RV groups such as RVE or RVH. These data suggest that along with RVA and RVC, screening for RVB is important and should be actively incorporated into testing schemes for etiological agents of gastrointestinal disease in pigs.

In piglets \<21 days of age, 65.7% of the positive samples appeared as single infections either with RVA (15.7%), RVB (12.9%), or RVC (37.1%). A few studies have investigated the prevalence of RVA and RVC in young piglets and have indicated that RVA infection in piglets is more common ([@bib28]). However, in the \<21 days old age group, we found a higher percentage of RVC (50.0%) than RVA (27.1%). Although this observation could be a sampling artifact, the difference may be due to differences in control measures between RVA and RVC. A RVA vaccine (ProSystems RCE, Merck, Whitehouse Station, NJ) is commercially available for swine in the United States, but there is no RVC vaccine available. While the vaccination rates are unknown, the use of vaccination may impact the detection and prevalence of RVA in pigs. Because of the lack of routine detection and commercial RVC vaccine, other control measures such as feedback have been used to prevent RVC infections. Feedback material contains feces or intestinal content from sows and piglets, which is fed to sows generally 4--8 weeks before farrowing to stimulate maternal antibodies to virus therein, providing passive immunity to piglets against RV ([@bib19]).

Multiple infections with genetically variant RVA strains have been described before, but the high detection rate (24.3%) of RVA/RVB/RVC co-infections in pigs described in our study is novel. A Japanese study suggested that comingling weaned piglets with different immunological histories and fecal shedding profiles exposes the piglets to other RVA strains, possibly leading to the multiple RVA infections within a group of nursery pigs ([@bib35]). We identified four sites from which samples positive for RVB contained more than one RVB genotype. In addition, the RT-PCR assay identified 13 RVB positive samples co-infected with more than one RVB genotypes, suggesting that like RVA, multiple RVB genotypes can co-circulate within one site. In addition to multiple RVB infections, our data shows that mixed RVA/RVB/RVC infections occur rather frequently as well. In samples from the 21--55 days old age group, there was a higher percentage of mixed RVA, RVB, and/or RVC infections than single RV infections. RVA was most prevalent; however, the prevalence of mixed RVA/RVB/RVC infections was higher than expected. In the \>55 days old age group, all samples tested were positive for RVA. The continuous infection with RVA even in more mature pigs may be due to the limited cross protection between RVA genotypes ([@bib19]). Nearly half (45.5%) of the samples from the \>55 days old age group were positive for all three RV groups tested, indicating that pigs continue to be infected throughout their early life, as suggested by [@bib35]. Interestingly, the presence of RVB only, RVC only, or the combination of these two were not found in the \>55 days old age group, and this leads to several intriguing hypotheses, all of which require further investigation. The presence of multiple RV group infections in older pigs might suggest that the ability of RVB or RVC to cause disease in older animals as the only infecting pathogen is more limited than the ability of RVA strains. Furthermore, this observation might suggest that pigs \>55 days of age do not get gastrointestinal disease upon infection with a single RV infection, and that a co-infection of multiple RV pathogens might be needed to cause clinical disease in these older animals. A primary confounding factor in interpreting any of these hypotheses is that there is no current method to measure existing RVB or RVC immunity. Until immunity measures are established, the significance, pathogenesis, and ecology of multiple RV infections in older pigs will remain unclear.

An 89% amino acid cut-off value was found to give a good correlation between G serotypes and G genotypes of RVA strain ([@bib31]). The standardized classification system for RVA developed by the RCWG, also demonstrated a good correlation between the 89% amino acid cut-off value and the 80% nucleotide cut-off value which is now routinely used to classify VP7 gene segments of RVA strains into 35 G genotypes ([@bib32]). In 2007, Martella and colleagues proposed a classification system for the RVC VP7 gene segment based on phylogenetic data, and an 89% amino acid cut-off value was used to distinguish at least 6 G genotypes ([@bib29]). Kuga and colleagues used phylogenetic analyses and pairwise identity frequency graphs of 50 VP7 sequences from RVB strains available at that time for classification purposes ([@bib24]). Their pairwise identity frequency graphs revealed the presence of two suitable nucleotide cut-off values: 67% and 76% (66% and 79% on the amino acid level). Our new analyses with a total of 125 available RVB VP7 sequences, including the 68 porcine RVB VP7 sequences determined in this study, revealed a similar picture as was observed by Kuga and colleagues, with slightly higher nucleotide cut-off values of 68% and 80%. Kuga and colleagues used their calculated 67% and 76% nucleotide cut-off values to define five G genotypes and 12 clusters inside these genotypes, respectively. Based on our most recent analyses, we propose to use the higher cut-off value (80%) to define 20 G genotypes for RVB VP7 sequences instead of the lower cut-off value (68%). The rationale for this choice is that recent data revealed that the genetic diversity observed among RVA and RVB strains is in the same order of magnitude ([@bib34], [@bib49]), opposed to what was believed in the past, which suggested that the genetic diversity of RVB strains was significantly larger than that of RVA strains ([@bib14], [@bib24], [@bib51]). Our calculations showed that an 80% nucleotide cut-off value for VP7 of RVB was suitable, which is the same cut-off used to distinguish G genotype for RVA strains. In addition, a correlation between RVB G genotypes based on a 68% nucleotide-cut-off value and G serotypes is unlikely. Although serotyping analyses cannot be performed currently due to the lack of methods to easily adapt RVB strains to cell culture, it is a reasonable assumption that a roughly similar number of amino acid changes (and corresponding nucleotide changes) have to be present into VP7 of RVA and RVB strains before two strains can be considered to belong to the distinct serotype ([@bib18]). Since the amino acid residues, which form the epitopes responsible for the serotypes specificity are currently unknown for VP7 of RVB, it is difficult to speculate on the correlation between RVB VP7 serotypes and genotypes, and future serological assays will have to be performed to study these relationships. However, for the majority of the RVB genotypes, the amino acid identity ranges among RVB strains belonging to the same G genotype is also above the 89% cut-off value ([Table 2](#t0010){ref-type="table"}) as previously determined for RVA strains to discriminate serotypes ([@bib36]). For a few RVB genotypes (G6, G7, G12 and G16) the diversity extends below the 89% amino acid cut-off value, which was also observed for a few RVA G genotypes such as G1--G4 and G6 ([@bib31]).

To update and maintain this classification system for RVB VP7 sequences, as has been done for RVA strains; we propose that assignment of VP7 sequences of RVB strains to potentially new genotypes would be done by the RCWG.

In conclusion, the genetic diversity observed in VP7 sequences of porcine RVB strains suggests that RVB strains have been circulating in the porcine population for a prolonged time (at least since the 1980s) and may be more prevalent than originally thought. Failure to detect the true prevalence of RVB strains in the swine population may have been overshadowed by the high prevalence, pathogenicity, and zoonotic properties of RVA. Most likely, the lack of proper diagnostic assays has concealed the prevalence of RVB in the swine population. The newly developed RT-PCR can be used to further increase our diagnostic capabilities and future studies on the molecular evolution and epidemiology of RVB strains. Moreover, the increased number of VP7 sequences of RVB strains has permitted its classification into 20 G genotypes based on an 80% nucleotide cut-off value by using the approach and guidelines used by the RCWG to classify RVA strains, leading to further comprehensive standardization of RV genotyping.

Materials and Methods {#s0020}
=====================

Origin of samples {#s0025}
-----------------

Animal intestinal samples, with a history or current outbreak of diarrhea, are routinely submitted to the University of Minnesota Veterinary Diagnostic Laboratory for determination of etiological agents of gastrointestinal disease. Clinical signs of gastrointestinal disease are most commonly weight-loss and diarrhea. Between the months of August and September 2009, 173 porcine samples (119 different sites) from 14 US states were submitted for diagnostic investigation. Porcine samples from all age groups are routinely screened for *S. enterica spp.*, TGEV, RVA, RVB, and RVC. Screening for other pathogens in addition to the former is dependent on the age of the pig: \<10 days, *E. coli*, *C. difficile* and *C. perfringens*; 11--55 days, *E. coli*, coccidia; \>55 days, *L. intracellularis* and *Brachyspira* species.

Electron microscopy {#s0030}
-------------------

Ammonium molybdate (1%) prepared in water (pH=5.5) was used to perform negative stain electron microscopy to visualize RV particles in samples.

Histopathology {#s0035}
--------------

Small intestine was fixed in 10% neutral buffered formalin, routinely processed, paraffin embedded, and stained using Harris\'s hematoxylin and eosin (HE).

Extraction of genomic material, reverse transcriptase-polymerase chain reaction (RT-PCR) amplification, and sequence analysis {#s0040}
-----------------------------------------------------------------------------------------------------------------------------

Approximately, 3 g of each sample was homogenized with 2 mL of gamma-irradiated Hyclone donor equine serum (Thermo Fisher Scientific, Waltham, MA) and centrifuged at 4200 rpm for 1 h. Subsequently, the viral RNA was extracted from the homogenate supernatant using an Ambion MagMax extraction kit (Applied Biosystems, Foster City, CA), according to the manufacturer\'s instructions.

Molecular detection by RT-PCR of RVA, RVC and TGEV were carried out according to controlled Standard Operating Procedures at the Veterinary Diagnostic Laboratory at the University of Minnesota (unpublished). The TGEV RT-PCR amplifies approximately a 120-nucleotide region of the spike protein while both the RVA and RVC RT-PCR amplify a region of approximately 100 nucleotides at the 3′ end of VP6.

For the molecular detection of RVB strains, a modified RT-PCR protocol was used as described by Matthijnssens et al., 2008 ([@bib31]). Oligonucleotide primers: Forward RVB VP7 (5′--GGA AAT AAT CAG AGA TGG CGT-3′, nucleotides 1--21 and Reverse RVB VP7 (5′--TCG CCT AGT CYT CTT TAT GC-3′, nucleotides 759--778) were designed based on sequence alignments of one human and five bovine RVB VP7 sequences listed in GenBank ([DQ869567](ncbi-n:DQ869567), [AB016818](ncbi-n:AB016818), [U84141](ncbi-n:U84141), [U84472](ncbi-n:U84472), [U84473](ncbi-n:U84473), and [AF531910](ncbi-n:AF531910) respectively). Briefly, 5 μL of extracted RNA was denatured at 97 °C for 3 min. RT-PCR was carried out using the Qiagen OneStep RT-PCR Kit (Qiagen/Westburg) in a GeneAmp PCR System 9700 thermal cycler (Applied Biosystems). Thermal cycling conditions were as follow: initial reverse transcription, 45 °C for 30 min; PCR activation, 95 °C for 15 min., followed by 40 cycles of 45 s at 94 °C, 45 s at 48 °C, and 2 min at 68 °C, and a final extension of 7 min at 72 °C.

PCR amplicons, approximately 800 nucleotide in length, were purified using the QIAquick PCR purification kit (Qiagen/Westburg) and sequencing was performed using a fully automated ABI 3730xl DNA Analyzer (Perkin-Elmer) with ABI BigDye Terminator version 3.1 chemistry (Perkin-Elmer). Primer sequence walking was performed to provide double coverage of each sequence. Sequences were analyzed using the Seqman 8.0 program of the Lasergene software (DNASTAR, Madison, WI). Sequences were aligned using Clustal W ([@bib48]) and phylogenetic analysis was performed using MEGA 5 ([@bib46]). Kimura 2-parameter correction at the nucleotide level and the Poisson correction parameter at the amino acid level were utilized to calculate the genetic distances ([@bib32]). The Neighbor-joining method was used to create the phylogenetic trees ([@bib41]). Amino acid profile and identity figures (sliding window of 1) were generated using Geneious Pro ([@bib13])

To obtain suitable cut-off values for evolution based classification of the RVB VP7 genome segment, the percentage identities between the complete ORFs of the available RVB VP7 genome segments in GenBank ([Supplemental Data 1](#s0045)), as well as the sequences determined in this study ([Table 1](#t0005){ref-type="table"}), were calculated using the pairwise distances program MEGA 5 ([@bib46]). Pairwise identity frequency graphs were constructed by plotting all the calculated pairwise identities in a graph with the percentage identities in the abscissa (*X*-axis) and the frequency of each of the calculated pairwise identities in the ordinate (*Y*-axis) ([@bib2]).

Appendix A. Supplementary materials {#s0050}
===================================

Supplementary Material

We thank Mary Thurn and Dr. Yin Jiang at the University of Minnesota Veterinary Diagnostic Laboratory (UMVDL) for their technical assistance. Research was supported by the UMVDL J.M was supported by an FWO ('Fonds voor Wetenschappelijk Onderzoek') post-doctoral fellowship.

Supplementary data associated with this article can be found in the online version at [10.1016/j.virol.2012.07.006](http://dx.doi.org/10.1016/j.virol.2012.07.006).
