key: cord-0828879-3o5pupiv authors: da Silva Medeiros, Thais Neris; Lorenzetti, Elis; Alfieri, Alice Fernandes; Alfieri, Amauri Alcindo title: Phylogenetic analysis of a G6P[5] bovine rotavirus strain isolated in a neonatal diarrhea outbreak in a beef cattle herd vaccinated with G6P[1] and G10P[11] genotypes date: 2014-11-07 journal: Arch Virol DOI: 10.1007/s00705-014-2271-4 sha: f0e4fa339fb2f5bbb79b0c78356e384d44000978 doc_id: 828879 cord_uid: 3o5pupiv The aim of this study was to perform the molecular characterization of the eleven genes of a G6P[5] bovine group A rotavirus (RVA) strain detected in a diarrhea outbreak from a vaccinated beef cattle herd. The outbreak affected 80 % of calves between 15–30 days old. RVA was identified by RT-PCR in 12 (70.6 %) out of 17 diarrheic fecal samples evaluated. The rotavirus wild-type strain had the genotype constellation G6(IV)-P[5](IX)-I2c-R2-C2-M2-A3-N2-T6-E2e-H3a. This study confirms the importance of homotypic immunity against the bovine RVA P[5] genotype in neonatal diarrhea in cattle herds that are regularly vaccinated against rotaviruses. Bovine group A rotavirus (RVA) is one of the main etiological agents of neonatal diarrhea in calves worldwide. Morbidity and mortality rates due to bovine RVA infections are high, which in turn cause important direct and indirect economic losses to beef and dairy production [1] . Rotavirus belongs to the family Reoviridae, and it is surrounded by a triple-layered protein capsid. The genome is formed by 11 double-stranded RNA segments, which encode six structural proteins (VP1-VP4, VP6, and VP7) and six non-structural proteins (NSP1-NSP5/6) [2] . Rotaviruses (RVs) are classified into eight distinct groups/ species (A-H) [3, 4] . The VP7 and VP4 proteins of RVA are located in the outer layer of the capsid and induce neutralizing antibodies. The antigenic variation of the genome segments that code for these proteins determines the binary RVA genotype classification system [2] . Currently, 27 G (VP7) and 37 P genotypes (VP4) of RVA have been described in mammals and avian species [5, 6] . In 2008, the Rotavirus Classification Working Group (RCWG) defined the notations Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx (x-Arabic numbers starting from 1) to VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5/ 6 genes of RV strains [5] . Until now, the following genotypes have been described: VP6, I1-I17; VP1, R1-R9; VP2, C1-C9; VP3, M1-M8; NSP1, A1-A18; NSP2, N1-N10; NSP3, T1-T12; NSP4, E1-E15; NSP5/6 , H1-H11 [6] [7] [8] [9] [10] . The most common combinations of G and P genotypes found in bovine RVA strains isolated from diarrhea episodes in calves are G6P [1] (Nebraska calf diarrhea virus [NCDV]-Lincoln), G6P [5] (UK), G8P [1] (A5), and G10P [11] (B223) [11] [12] [13] [14] . These genotypes have been described in Brazil with frequencies of 9.7 % (3/31) to 33.3 % (12/36), 8.6 % (3/35) to 40 % (20/50), 16.7 % (6/ 36), and 12.9 % (4/31) to 16 % (8/50), respectively [15] [16] [17] . G6P [5] has been reported to be a common genotype in bovine RVA strains and has been found in cattle herds from Asia [18] , Europe [14] , Oceania [13] , and America [19, 20] . The RVA G6 genotype is divided into five lineages (I-V). The bovine RVA strains can be found in the lineages G6-II, G6-III, G6-IV, and G6-V. The P [5] [21] . Diarrhea in calves can be prevented by appropriate management practices and programs for vaccination of pregnant cows [16] . However, vaccine failures may occur for many reasons, such as reassortment between strains with different genes [22] , antigenic or genetic diversity [23] , interaction between different or less common genotypes [16, 24, 25] , insufficient heterologous immunity [26] , and inappropriate vaccination management [27] . The aim of this study was to perform a molecular characterization of the eleven genes of a G6P [5] bovine RVA strain detected in a neonatal diarrhea outbreak in a beef cattle herd that had been regularly vaccinated with G6P [1] and G10P [11] strains. The beef cattle farm was located in Mato Grosso do Sul state, Central-West Brazil. The herd was managed on extensive pasture and had good health and nutritional practices. The farm had adopted a breeding and, consequently, birthing season of 90 days. In total, 600 cows were divided in lots ranging from 100 to 120 animals. During the neonatal diarrhea outbreak, there were 170 crossbred calves (Nelore 9 Angus), 15 to 30 days old. Every cow was regularly vaccinated for neonatal diarrhea control with an inactivated commercial vaccine containing the bovine RVA genotypes G6P [1] (NCDV-Lincoln strain) and G10P [11] (B223 strain), bovine coronavirus (Hanson strain), Clostridium perfringens type C (toxoid/beta toxin, NL-1003 strain), and Escherichia coli pili adherence factor (K99 strain). The vaccination schedule was according to the manufacturer's instructions. Diarrheic fecal samples (n = 17) were collected from calves aged less than 30 days old. The fecal samples were stored at -20°C until analysis. The study was submitted to the Ethics Committee on Animal Experiments of the Universidade Estadual de Londrina and approved (46/09). Nucleic acid was extracted using a combination of phenol/chloroform/isoamyl alcohol (25:24:1) and silica/ guanidinium isothiocyanate [28] extraction methods. All diarrheic fecal samples was subjected to RT-PCR assay with previously reported consensus primers to amplify VP7 (1,062 bp) and VP4 (VP8*) (876 bp) genes to determine the G and P genotype, respectively [29] [30] [31] . A diarrheic fecal sample (BRA1532) was selected according to the quality score for analysis of the eleven genes using previously described primers [32] [33] [34] [35] [36] . The RT-PCR products were purified using an Illustra GFX PCR DNA and Gel Band Purification Kit and sequenced using an ABI3500 Genetic Analyzer sequencer with a BigDye Ò Terminator v3.1 Cycle Sequencing Kit. Nucleotide quality analysis and contig assembly of the RVA gene sequences were performed with Phred and CAP3 software, respectively. The identity matrix was constructed using BioEdit software version 7.0.8.0. Phylogenetic trees based on the nucleotides were obtained using the neighbor-joining method in MEGA v6 software. The bootstrapping probabilities were calculated using 1,000 replicates. The sequences described in this study are available in GenBank with the accession numbers JQ943546 to JQ943578. Eighty percent (136/170) of the calves that were 15 to 30 days old had diarrhea. All diarrheic calves were treated with broad-spectrum antibiotics by the parenteral route but were unresponsive to antibiotic therapy, and 2.94 % (4/ 136) of the calves died. Before analysis for bovine RVA, the diarrheic fecal samples were evaluated by the modified Ziehl-Neelsen technique [37] and SN-PCR assay [38] for Cryptosporidium spp. and bovine coronavirus detection, respectively, and the results were negative (data not shown). RT-PCR assays using consensus primers of VP7 and VP4 (VP8*) genes showed that 12 (70.6 %) out of 17 diarrheic fecal samples analyzed were RVA positive. The nucleotide sequence analysis of the VP7 and VP4 (VP8*) genes of 12 Brazilian wild-type bovine RVA strains showed the highest nucleotide sequence identity (94.9 %) to G6 genotype lineage IV (OH-4 equine strain) (Fig. 1a) , and 92 % with the bovine strain VMRI, which belongs to the P[5] genotype lineage VIII (Fig. 1b) . The cutoff value to be considered of the same lineage of the P [5] genotype is 96 % nucleotide sequence identity [21] ; therefore, we propose a new lineage, named P[5]-IX. In the phylogenetic tree, the BRA1532 strain formed a new branch separated from the other lineages of the P [5] genotype. For the VP6, VP1-VP3, and NSP1-NSP5/6 genes, the nucleotide sequence analysis of the RT-PCR products of the bovine RVA BRA1532 field strain revealed the presence of the following genotypes: G6(IV)-P [5] (IX)-I2c-R2-C2-M2-A3-N2-T6-E2e-H3a. The bovine WC3 strain used in the human vaccine and the BRA1532 strain has the same constellation of genotypes. However, the Brazilian field strain showed high nucleotide sequence identity to the WC3 strain only in the VP1 and VP2 genes ( Table 1) . The 70.6 % (12/17) of diarrheic fecal samples that were positive for bovine RVA in this study suggested that rotavirus was the etiological agent involved in this neonatal diarrhea outbreak. Additionally, all samples were negative for common enteropathogens such as bovine coronavirus and Cryptosporidium spp. (data not shown), and the diarrheic calves were unresponsive to two or three doses of broad-spectrum antibiotic therapy [39] . The nucleotide sequence analysis of the VP1-VP3, VP4, VP6-VP7, and NSP1-NSP5/6 genes of the BRA1532 strain showed the same genome constellation described in the bovine WC3 (G6P [5] ) strain; however, the Brazilian bovine RVA strain belonged to the P [5] genotype lineage IX, and the WC3 strain belonged to the P [5] genotype lineage IV. The BRA1532 strain displayed the T6 (NSP3 gene) genotype, which was distinct from the gene found in the bovine UK strain (G6P [5] ), which belongs to the T7 genotype. Additionally, the VP4 (VP8*) and NSP5/6 genes of the BRA1532 strain belonged to different lineages (P[5]-IX and H3a) when compared with the UK strain (P[5]-V and H3e) [6, 40] . According to the NCDV-Lincoln strain, the field Brazilian strain shared the same genotypes and lineages for all genes analyzed, which are present in the commercial vaccine used in this beef cattle herd, with the exception of the VP4 (VP8*) gene. There is significant homotypic divergence within genotype G6 of RVA strains identified in human and bovine hosts, suggesting the existence of antigenic differences between G6 strains [41, 42] . Homologous immunity occurs within the same G6 genotype, and even a small genotypic difference is enough to result in insufficient protection [26] . Some studies have shown that vaccination with the prototype G6P [1] of bovine RVA results in poor heterologous protection against RVA G6 strains containing different P genotypes from the vaccine [27, 43] . Moreover, commercial vaccines containing the genotype G6P [1] may not be as effective, as P [1] may not be the most common genotype depending on the geographic region studied [44] . The vaccine pressure may generate competition, selection or variation between strains [45] , which could be responsible for the introduction or emergence of novel adapted strains or differences in genotypes [23, 45] . Different genotypes have caused diarrhea in calves born to dams vaccinated against rotavirus due to poor heterologous protection against a reassortment between B641 (G6P [5] ) and B223 (G10P [11] ) bovine RVA strains [24] . However, it has been reported that immunity against the different G6 strains is due, at least in part, to heterologous protection against another lineage, even when the P genotype is different [41] . A polyvalent vaccine may not offer protection against all of the currently circulating bovine RVA genotypes, and multiple G and P genotypes can circulate simultaneously in a given population. Their impact on morbidity and mortality rates, high treatment costs, and reduced growth rates would not be minimized [45] . Monitoring the G and P genotypes of bovine RVA strains that are circulating or involved in outbreaks of neonatal calf diarrhea in a vaccinated cattle herd should be performed more often for a comprehensive analysis of the genetic diversity of RVA. It should also be emphasized that a new classification system would also help in understanding the inefficient protection. In conclusion, the G6-IV and P[5]-IX genotypes were identified in the bovine RVA field strain as the cause of a diarrhea outbreak in a vaccinated beef cattle herd. These genotypes were distinct from the G and P genotypes of bovine RVA strains included in the commercial vaccine used in the herd. These results highlight the importance of homotypic immunity and monitoring of the bovine RVA genotypes circulating in cattle populations. Some infectious causes of diarrhea in young farm animals Family: Reoviridae VP6-sequencebased cutoff values as a criterion for rotavirus species demarcation Recommendations for the classification of group A rotaviruses using all 11 genomic RNA segments Uniformity of rotavirus strain nomenclature proposed by the rotavirus classification working group (RCWG) Identification of an avian group A rotavirus containing a novel VP4 gene of close relationship to those of mammalian rotaviruses Full genomic analysis of rabbit rotavirus G3P[14] strain N5 in China: identification of a novel VP6 genotype Novel NSP4 genotype in a Camel G10P[15] rotavirus strain Novel NSP1 genotype characterized in an African camel G8P[11] rotavirus strain Relative frequencies of G and P types among rotaviruses from Indian diarrheic cow and buffalo calves Changes in the prevalence of rotavirus G and P types in diarrheic calves from the Kagoshima prefecture in Japan Detection and analysis of bovine rotavirus strains circulating in Australian calves during Diversity and zoonotic potential of rotaviruses in swine and cattle across Europe G and P genotypes of group A rotavirus strains circulating in calves in Brazil G and P genotypes of group A rotavirus from diarrhoeic calves born to cows vaccinated against the NCDV (P[1], G6) rotavirus strain Molecular characterization of G and P-types bovine rotavirus strains from Goiás, Brazil: high frequency of mixed P-type infections The distribution of G and P types within isolates of bovine rotavirus in Japan Molecular characterization of bovine rotavirus circulating in beef and dairy herds in Argentina during a 10-year period Rotavirus in cattle: risk factors, prevalence and antigenic characterization from dairy calves samples in São Paulo State, Brazil Phylogenetic analysis of typical bovine rotavirus genotypes G6, G10, P[5] and P[11] circulating in Argentinean beef and dairy herds Zoonotic aspects of rotaviruses Detection and genetic characterization of group A rotavirus strains circulating among children with acute gastroenteritis in Japan Serological and genotypic characterization of group A rotavirus reassortants from diarrheic calves born to dams vaccinated against rotavirus Genotypification of bovine group A rotavirus in México Antigenic relationships among some bovine rotaviruses: serum neutralization and cross-protection in gnotobiotic calves Distribution of G (VP7) and P (VP4) genotypes of group A bovine rotaviruses from Turkish calves with diarrhea Frequency of group A rotavirus in diarrhoeic calves in Brazilian cattle herds Polymerase Chain Reaction amplification and typing of rotavirus nucleic acid from stool specimens Identification of group A rotavirus gene 4 types by polymerase chain reaction Identification of a novel VP4 genotype carried by a serotype G5 porcine rotavirus strain Characterization of VP1, VP2 and VP3 gene segments of a human rotavirus closely related to porcine strains Molecular characterization of VP6 genes of human rotavirus isolates: correlation of genogroups with subgroups and evidence of independent segregation Full genomebased classification of rotaviruses reveals a common origin between human Wa-like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains NSP4 gene analysis of rotaviruses recovered from infected children with and without diarrhea Nucleotide sequence analysis of rotavirus gene 11 from two tissue culture-adapted ATCC strains, RRV and Wa Staining of Cryptosporidium by a modified Ziehl-Nielsen technique Improved detection of bovine coronavirus N gene in faeces of calves infected naturally by a semi-nested PCR assay and an internal control Risk factors associated with diarrhoea in newborn calves Molecular characterization of genotype G6 human rotavirus strains detected in Italy from Comparative sequence analysis of the VP7 genes of G6, G8 and G10 bovine group A rotarviruses and further characterization of G6 subtypes Characterization and phylogenetic analysis of the VP7 proteins of serotype G6 and G8 human rotaviruses An unusual group-A rotavirus associated with an epidemic of diarrhea among threemonth-old calves Characterization of mixed infections with different strains of bovine rotavirus in an outbreak of diarrhea in dairy herds in Goiás, Brazil Changing profile of the bovine rotavirus G6 population in the south of Ireland from