2018–2020年人轮状病毒锦州地方株VP4及VP7基因特征

【背景】人A组轮状病毒(Rotavirus Group A,RVA)是婴幼儿胃肠炎的主要病原体及发展中国家婴幼儿死亡的重要原因,目前无特效药物治疗,疫苗预防是唯一可行的预防感染方法。外衣壳蛋白VP7和VP4是疫苗设计的主要靶点,针对该基因加强RVA地方株分子流行病学监测十分必要。【目的】对锦州地方流行RVA株VP7和VP4基因进行型别鉴定和序列特征分析。【方法】收集锦州地区2018-2020年RVA感染腹泻患儿的粪便标本,提取病毒RNA,通过RT-PCR扩增VP7、VP4基因片段并测序,得到7株RVA VP7和VP4序列。使用在线基因分型工具Rota C V2.0对测序结果进行分型分析。应用BLAST、DNAStar、MEGA X、Bio Edit等生物软件与临床流行株及疫苗株进行系统发育分析及氨基酸序列比对分析。【结果】分型结果表明7株锦州地方株均为G9P[8]型,系统发育分析证实其VP7和VP4基因分别属于G9-Ⅵ和P[8]-3谱系,核苷酸序列相似性分别为99.32%-100%与99.41%-100%。JZ株VP7与疫苗株Rotavac和Rotasiil相比,在抗原表位区7-1a、7-1b、7-2中分别存在4个和3个氨基酸替换。JZ株VP4与疫苗株Rotarix和Rota Teq VP4氨基酸序列相比,发现7个和4个氨基酸替换,位于抗原表位区8-1和8-3。【结论】2018-2020年在辽宁锦州地区检测到7株G9P[8]型RVA株,VP7和VP4序列相似性高于99%,G9P[8]型可能是辽宁省锦州地区2018-2020年婴幼儿轮状病毒腹泻的主要流行基因型之一。与同基因型疫苗株比较,位于JZ株VP7和VP4抗原表位区的氨基酸位点差异对于野毒株免疫逃逸机制的研究具有意义。     

Heterogeneity and temporal dynamics of evolution of G1 human rotaviruses in a settled population

A rotavirus sample collection from 19 consecutive years was used to investigate the heterogeneity and the dynamics of evolution of G1 rotavirus strains in a geographically defined population. Phylogenetic analysis of the VP7 gene sequences of G1P[8] human rotavirus strains showed the circulation of a heterogeneous population comprising three lineages and seven sublineages. Increases in the circulation of G1 rotaviruses were apparently associated with the introduction of novel G1 strains that exhibited multiple amino acid changes in antigenic regions involved in rotavirus neutralization compared to the strains circulating in the previous years. The emergence and/or introduction of G1 antigenic variants might be responsible for the continuous circulation of G1 rotaviruses in the local population, with the various lineages and sublineages appearing, disappearing, or cocirculating in an alternate fashion under the influence of immune-pressure mechanisms. Sequence analysis of VP4-encoding genes of the G1 strains revealed that the older strains were associated with a unique VP4 lineage, while a novel VP4 lineage emerged after 1995. The introduction of human rotavirus vaccines might alter the forces and balances that drive rotavirus evolution and determine the spread of novel strains that are antigenically different from those included in the vaccine formulations. The continuous emergence of VP7-VP4 gene combinations in human rotavirus strains should be taken into consideration when devising vaccination strategies.     

G and P genotypes of rotavirus circulating among children with diarrhea in the Colombian northern coast.

A study on the prevalence of rotavirus G and P genotypes was carried out based on 253 stool specimens obtained from children living in the Colombia northern coast region who were less than 3-years-old and who suffered from acute diarrhea. A previous study had detected the presence of rotavirus A in 90 (36.5%) of the 246 samples tested by enzyme immunoassay (EIA), and these strains were investigated in the present study. Of these, 50 strains yielded an RNA electropherotype, most of which (80.0%) had long profiles and 20.0% of which had short profiles. Genotyping of 84 positive samples indicated that 67.9% of the strains could be typed. G1 (57.9%), was the most predominant VP7 genotype, followed by G3 (21.1%), G9 (15.8%) and G2 (5.3%). Among the VP4 genotypes, P[4] (49.1%) was the most prevalent, followed by P[6] 36.4% and P[8] (14.5%). Neither G4 nor G8 nor P[9] types were detected. The most common G-P combinations were G3 P[4] (8.8%) and G9 P[6] (7.0%), followed by G1 P[4] and G1 P[8] (5.3% each). All G1 P[8] strains showed long RNA profiles, whereas G3 P[4] and G9 P[6] displayed both long and short patterns. Mixed infections involved 21.0% of strains. There was a marked diversity among strains collected, and novel strains, including G9, as well as other atypical combinations of G and P genotypes, such as G9 P[6] and G3 P[4], were found.     

Ocurrence of rotavirus and picobirnavirus in wild and exotic avian from amazon forest

The present study reports the occurrence of rotavirus A (RVA), rotavirus D (RVD), rotavirus F (RVF), rotavirus G (RVG), and picobirnavirus (PBV) in fecal specimens of wild (n = 22), and exotic birds (n = 1) from different cities of Pará state. These animals were hospitalized at Veterinary Hospital of the Federal University of Pará, Brazil, in a period from January 2018 to June 2019. The animals exhibited different clinical signs, such as diarrhea, malnutrition, dehydration, and fractures. The results showed 39.1% (9/23) of positivity for RVA by RT-qPCR. Among these, one sample (1/9) for the NSP3 gene of T2 genotype was characterized. About 88.9% (8/9) for the VP7 gene belonging to G1, G3 equine like and G6 genotypes, and 55.5% (5/9) for the VP4 gene of P[2] genotype were obtained. In the current study, approximately 4.5% of the samples (1/23) revealed coinfection for the RVA, RVD and RVF groups. Furthermore, picobirnavirus (PBV) was detected in one of the 23 samples tested, and was classified in the Genogroup I. The findings represent the first report of RVA, RVD, RVF, RVG, and PBV genotypes in wild birds in Brazil, and due to wide distribution it can implies potential impacts of RVs, and PBVs on avian health, and other animals contributing to construction of new knowledge, and care perspectives.     

Distribution of human rotaviruses, especially G9 strains, in Japan from 1996 to 2000

A 4-year (1996-2000) survey of rotavirus infection involving 2,218 diarrheal fecal specimens of children collected from five regions of Japan was conducted. A total of 642 (28.9%) specimens were found to be rotavirus positive. A changed prevalence pattern of rotavirus G serotype was found with an increase of G9 and G2 and a decrease of G1, although G1 remained the prevailing serotype. Serotype G9 was unexpectedly determined to be the prevailing serotype in Sapporo (62.5%) and Tokyo (52.9%) in 1998-1999, and in Saga (78.4%) in 1999-2000. G9 strains isolated from 1998-1999 belonged to the P[8]-NSP4-Wa-group with long RNA pattern, while, G9 strains isolated from 1999-2000 belonged to three groups, the P[8]-NSP4-Wa-group with long RNA pattern, the P[4]-NSP4-KUN-group with short RNA pattern and a mixed-type group (P[4]/P[8]-NSP4-KUN/Wa-group with long RNA pattern). Both sequence and immunological analysis of VP7 revealed that the G9 strains from 1999-2000 were much more closely related to the G9 strains isolated worldwide in the 1990s, including G9 strains found in Thailand in 1997. However, the G9 strains from 1998-1999 were distinct from these and more closely related to the G9 prototype strains F45, AU32 and WI61 discovered in Japan and the US in the 1980s. Thus the G9 strains isolated in 1998-1999 had progenitors common to the G9 prototype strains, while the strains isolated in 1999-2000 did not directly evolve from them but were related to global G9 strains that have emerged in recent years. These data supported our previous report that G9 rotavirus might exist as two or more subtypes with diverse RNA patterns, P-genotype and NSP4 genogroup combinations (Y.M. Zhou et al., J. Med. Virol. 65: 619-628, 2001) and suggested that G9 rotavirus prevalent in Japan during two successive years belonged to different subtypes. The nucleotide sequences presented in this paper were submitted to DDBJ, EMBL and GenBank nucleotide sequence databases. The accession numbers are: 00-Ad2863VP7 (AB091746), 00-OS2986VP7 (AB091747), 00-SG2509VP7 (AB091748), 00-SG2518VP7 (AB091749), 00-SG2541 (AB091750), 00-SG2864 (AB091751), 00-SP2737VP7 (AB091752), 99-SP1542VP7 (AB091753), 99-SP1904VP7 (AB091754), 99-TK2082VP7 (AB091755) and 99-TK2091VP7 (AB091756).     

Frequent reassortments may explain the genetic heterogeneity of rotaviruses: analysis of Finnish rotavirus strains

The predominant rotavirus electropherotypes (e-types) during 17 epidemic seasons (1980 through 1997) in Finland were established, and representative virus isolates were studied by nucleotide sequencing and phylogenetic analysis. The virus isolates were either P[8]G1 or P[8]G4 types. The G1 and G4 strains formed one G1 lineage (VP7-G1-1) and one G4 lineage, respectively. Otherwise, they belonged to two P[8] lineages (VP4-P[8]-1 and -2) unrelated to their G types. Phylogenetic analysis of partial sequences of all 11 RNA segments obtained from the strains also revealed genetic diversity among gene segments other than those defining P and G types. With the exception of segments 1, 3, and 10, the sequences of the other segments could be assigned to 2 to 4 different genetic clusters. The results of this study suggest that, in addition to the RNA segments encoding VP4 and VP7, the other RNA segments may segregate independently as well. In total, the 9 predominant e-types represented 7 different RNA segment combinations when the phylogenetic clusters of their 11 genes were determined. The extensive genetic diversity and number of e-types among rotaviruses are best explained by frequent genetic reassortment.     

Reassortment of Human and Animal Rotavirus Gene Segments in Emerging DS-1-Like G1P[8] Rotavirus Strains

The emergence and rapid spread of novel DS-1-like G1P[8] human rotaviruses in Japan were recently reported. More recently, such intergenogroup reassortant strains were identified in Thailand, implying the ongoing spread of unusual rotavirus strains in Asia. During rotavirus surveillance in Thailand, three DS-1-like intergenogroup reassortant strains having G3P[8] (RVA/Human-wt/THA/SKT-281/2013/G3P[8] and RVA/Human-wt/THA/SKT-289/2013/G3P[8]) and G2P[8] (RVA/Human-wt/THA/LS-04/2013/G2P[8]) genotypes were identified in fecal samples from hospitalized children with acute gastroenteritis. In this study, we sequenced and characterized the complete genomes of strains SKT-281, SKT-289, and LS-04. On whole genomic analysis, all three strains exhibited unique genotype constellations including both genogroup 1 and 2 genes: G3-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2 for strains SKT-281 and SKT-289, and G2-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2 for strain LS-04. Except for the G genotype, the unique genotype constellation of the three strains (P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2) is commonly shared with DS-1-like G1P[8] strains. On phylogenetic analysis, nine of the 11 genes of strains SKT-281 and SKT-289 (VP4, VP6, VP1-3, NSP1-3, and NSP5) appeared to have originated from DS-1-like G1P[8] strains, while the remaining VP7 and NSP4 genes appeared to be of equine and bovine origin, respectively. Thus, strains SKT-281 and SKT-289 appeared to be reassortant strains as to DS-1-like G1P[8], animal-derived human, and/or animal rotaviruses. On the other hand, seven of the 11 genes of strain LS-04 (VP7, VP6, VP1, VP3, and NSP3-5) appeared to have originated from locally circulating DS-1-like G2P[4] human rotaviruses, while three genes (VP4, VP2, and NSP1) were assumed to be derived from DS-1-like G1P[8] strains. Notably, the remaining NSP2 gene of strain LS-04 appeared to be of bovine origin. Thus, strain LS-04 was assumed to be a multiple reassortment strain as to DS-1-like G1P[8], locally circulating DS-1-like G2P[4], bovine-like human, and/or bovine rotaviruses. Overall, the great genomic diversity among the DS-1-like G1P[8] strains seemed to have been generated through reassortment involving human and animal strains. To our knowledge, this is the first report on whole genome-based characterization of DS-1-like intergenogroup reassortant strains having G3P[8] and G2P[8] genotypes that have emerged in Thailand. Our observations will provide important insights into the evolutionary dynamics of emerging DS-1-like G1P[8] strains and related reassortant ones.     

我国1998～1999年流行的婴幼儿腹泻轮状病毒的分型研究

轮状病毒是世界范围内引起婴幼儿腹泻的主要病原。根据病毒外壳蛋白VP4和VP7抗原性的不同可区分为不同型：P（VP4,protease sensitive）型和G（VP7,glycoprotein）型。1998－1999年在中国8个城市（长春、秦皇岛、北京、杭州、福州、广州、成都、昆明）采集了急性腹泻患儿的1093份非细菌性腹泻粪便标本,先进行A组轮状病毒（HRV）的筛选,其中阳性标本433份（39.6％）,电泳型长型占优势（96％）。对HRV标本,再利用血清型特异的MAbELISA和／或RT－PCR进行G分型。结果表明,在1998－1999年,在上述8城市非细菌性腹泻行季节,以HRV G1型为主要流行株,占阳性的83.4％,其次为G3（12.0%）、G4（3.5%)和G2(3.2%）。此外,有3份（0.7％）HRV标本未能分型,12（2.8％）份标本为混合感染,还结合1982－1996年全国12个地区1382份HRV标本的分型资料,分析了我国HVR G血清型的流行规律。实验中又抽样选取了124份GHRV标本,用RT－PCR进行P分型,其中P[8]型76份（61.3％）,P[4]型14份（11.3%）,P[6]型12份(9.7%),P[9]型8份(6.4%)。另外15份（12.1%)未能分出P型,有待进一步检定,实验中HRV分离株除了觉见的P[8]G1(51.4%)、P[4]G2(4.6%)毒株外,还检测到P[8]G3(11.0%)、P[8]G4(6.4%) 和其它较少见的病毒型。以上结果为我国轮状病毒疫苗的应用和开发提供了较系统、清晰的流行病学背景资料。     

Full genome-based classification of rotaviruses reveals a common origin between human Wa-Like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains

Group A rotavirus classification is currently based on the molecular properties of the two outer layer proteins, VP7 and VP4, and the middle layer protein, VP6. As reassortment of all the 11 rotavirus gene segments plays a key role in generating rotavirus diversity in nature, a classification system that is based on all the rotavirus gene segments is desirable for determining which genes influence rotavirus host range restriction, replication, and virulence, as well as for studying rotavirus epidemiology and evolution. Toward establishing such a classification system, gene sequences encoding VP1 to VP3, VP6, and NSP1 to NSP5 were determined for human and animal rotavirus strains belonging to different G and P genotypes in addition to those available in databases, and they were used to define phylogenetic relationships among all rotavirus genes. Based on these phylogenetic analyses, appropriate identity cutoff values were determined for each gene. For the VP4 gene, a nucleotide identity cutoff value of 80% completely correlated with the 27 established P genotypes. For the VP7 gene, a nucleotide identity cutoff value of 80% largely coincided with the established G genotypes but identified four additional distinct genotypes comprised of murine or avian rotavirus strains. Phylogenetic analyses of the VP1 to VP3, VP6, and NSP1 to NSP5 genes showed the existence of 4, 5, 6, 11, 14, 5, 7, 11, and 6 genotypes, respectively, based on nucleotide identity cutoff values of 83%, 84%, 81%, 85%, 79%, 85%, 85%, 85%, and 91%, respectively. In accordance with these data, a revised nomenclature of rotavirus strains is proposed. The novel classification system allows the identification of (i) distinct genotypes, which probably followed separate evolutionary paths; (ii) interspecies transmissions and a plethora of reassortment events; and (iii) certain gene constellations that revealed (a) a common origin between human Wa-like rotavirus strains and porcine rotavirus strains and (b) a common origin between human DS-1-like rotavirus strains and bovine rotaviruses. These close evolutionary links between human and animal rotaviruses emphasize the need for close simultaneous monitoring of rotaviruses in animals and humans.     

Distribution of conserved and specific epitopes on the VP8 subunit of rotavirus VP4.

Three cDNA clones comprising the VP8 subunit of the VP4 of human rotavirus strain KU (VP7 serotype G1; VP4 serotype P1A) G1 were constructed. The corresponding encoded peptides were designated according to their locations in the VP8 subunit as A (amino acids 1 to 102), B (amino acids 84 to 180), and C (amino acids 150 to 246 plus amino acids 247 to 251 from VP5). In addition, cDNA clones encoding peptide B of the VP8 subunit of the VP4 gene from human rotavirus strains DS-1 (G2; P1B) and 1076 (G2; P2) were also constructed. These DNA fragments were inserted into plasmid pGEMEX-1 and expressed in Escherichia coli. Western immunoblot analysis using antisera to rotavirus strains KU (P1A), Wa (P1A), DS-1 (P1B), 1076 (P2), and M37 (P2) demonstrated that peptides A and C cross-reacted with heterotypic human rotavirus VP4 antisera, suggesting that these two peptides represent conserved epitopes in the VP8 subunit. In contrast, peptide B appears to be involved in the VP4 serotype and subtype specificities, because it reacted only with the corresponding serotype- and subtype-specific antiserum. Antiserum raised against peptide A, B, or C of strain KU contained a lower level of neutralizing activity than did that induced by the entire VP8 subunit. In addition, the serotype-specific neutralizing activity of anti-KU VP8 serum was ablated after adsorption with the KU VP8 protein but not with a mixture of peptides A, B, and C of strain KU, suggesting that most of the serotype-specific epitopes in the VP8 subunit are conformational and are dependent on the entire amino acid sequence of VP8.     

Initial interaction of rotavirus strains with N-acetylneuraminic (sialic) acid residues on the cell surface correlates with VP4 genotype, not species of origin

We examined 41 human and animal rotavirus strains representative of all known P genotypes for their dependency on cellular N-acetylneuraminic (sialic) acid (SA) residues for infectivity. Our results showed that all rotaviruses studied, whether of animal or human origin, belonging to P genotypes [1], [2], [3], and [7] depended on SA residues on the cell surface for efficient infectivity but that all human and animal rotavirus strains representative of the remaining known P genotypes were SA independent. The SA residue requirement for efficient infectivity did not change for reassortant rotavirus strains with altered VP4-VP7 combinations. The initial interaction of rotavirus strains with SA residues on the cell surface correlated with VP4 genotype specificity, not with species of origin or VP7 G serotype specificity (P = 0.001; r2 = 1.00, Pearson's correlation coefficient). In addition to being a requirement for infectivity, the presence of SA residues on the cell surface is a requirement for efficient growth in cell culture; recognition of the association of specific P genotypes with the binding of rotavirus to SA residues will facilitate our understanding of the molecular basis of the early events of rotavirus-cell interactions in cell culture models and of pathogenicity in vivo.     

Genetic variation in the VP4 and NSP4 genes of human rotavirus serotype 3 (G3 type) isolated in China and Japan

Sequence analyses of the VP4 and NSP4 genes were performed on twenty human isolates of serotype G3 rotavirus obtained from China and Japan. One isolate from China, CHW17, possessed P[4] genotype VP4 and KUN group NSP4 genes which are associated with G2. One isolate (02/92) from Japan, which was shown to have a wider spacing between RNA segments 10 and 11 by RNA polyacrylamide gel electrophoretic analysis like AU-1, possessed P[9] genotype VP4 and AU-1 group NSP4 genes. The other isolates had P[8] genotype VP4 and Wa group NSP4 genes. While the nucleotide sequence conservation among the G3 VP7 genes was more than 79% (Wen et al, Arch. Virol., 1997, 142: 1481-1489), the conservation of VP4 and NSP4 genes in the same genotypes or groups was more than 85%.     

Genotype distribution and evolutionary analysis of rotavirus associated with acute diarrhea outpatients in Hubei,China, 2013–2016

Group A human rotaviruses (RVAs) annually cause the deaths of 215,000 infants and young children. To understand the epidemiological characteristics and genetic evolution of RVAs, we performed sentinel surveillance on RVA prevalence in a rotavirus-surveillance network in Hubei, China. From 2013 to 2016, a total of 2007 fecal samples from hospital outpatients with acute gastroenteritis were collected from four cities of Hubei Province. Of the 2007 samples, 153 (7.62%) were identified positive for RVA by real-time RT-PCR. RVA infection in Hubei mainly occurred in autumn and winter. The highest detection rate of RVA infection was in 1–2 years old of outpatients (16.97%). No significant difference of RVA positive rate was observed between females and males. We performed a phylogenetic analysis of the G/P genotypes based on the partial VP7/VP4 gene sequences of RVAs. G9P[8] was the most predominant strain in all four years but the prevalence of G2P[4] genotype increased rapidly since 2014. We reconstructed the evolutionary time scale of RVAs in Hubei, and found that the evolutionary rates of the G9, G2, P[8], and P[4] genotypes of RVA were 1.069×10^-3, 1.029×10^-3, 1.283×10^-3 and 1.172×10^-3 nucleotide substitutions/site/year, respectively. Importantly, using a molecular clock model, we showed that most G9, G2, P[8], and P[4] genotype strains dated from the recent ancestor in 2005, 2005, 1993, and 2013, respectively. The finding of the distribution of RVAs in infants and young children in Hubei Province will contribute to the understanding of the epidemiological characteristics and genetic evolution of RVAs in China.     

Isolation of human monoclonal antibodies that neutralize human rotavirus

A human antibody library constructed by utilizing a phage display system was used for the isolation of human antibodies with neutralizing activity specific for human rotavirus. In the library, the Fab form of an antibody fused to truncated cp3 is expressed on the phage surface. Purified virions of strain KU (G1 serotype and P[8] genotype) were used as antigen. Twelve different clones were isolated. Based on their amino acid sequences, they were classified into three groups. Three representative clones-1-2H, 2-3E, and 2-11G-were characterized. Enzyme-linked immunosorbent assay with virus-like particles (VLP-VP2/6 and VLP-VP2/6/7) and recombinant VP4 protein produced from baculovirus recombinants indicated that 1-2H and 2-3E bind to VP4 and that 2-11G binds to VP7. The neutralization epitope recognized by each of the three human antibodies might be human specific, since all of the antigenic mutants resistant to mouse monoclonal neutralizing antibodies previously prepared were neutralized by the human antibodies obtained here. After conversion from the Fab form of an antibody into immunoglobulin G1, the neutralizing activities of these three clones toward various human rotavirus strains were examined. The 1-2H antibody exhibited neutralizing activity toward human rotaviruses with either the P[4] or P[8] genotype. Similarly, the 2-3E antibody showed cross-reactivity against HRVs with the P[6], as well as the P[8] genotype. In contrast, the 2-11G antibody neutralized only human rotaviruses with the G1 serotype. The concentration of antibodies required for 50% neutralization ranged from 0.8 to 20 micro g/ml.     

Full genomic analysis of human rotavirus strain B4106 and lapine rotavirus strain 30/96 provides evidence for interspecies transmission

The Belgian rotavirus strain B4106, isolated from a child with gastroenteritis, was previously found to have VP7 (G3), VP4 (P[14]), and NSP4 (A genotype) genes closely related to those of lapine rotaviruses, suggesting a possible lapine origin or natural reassortment of strain B4106. To investigate the origin of this unusual strain, the gene sequences encoding VP1, VP2, VP3, VP6, NSP1, NSP2, NSP3, and NSP5/6 were also determined. To allow comparison to a lapine strain, the 11 double-stranded RNA segments of a European G3P[14] rabbit rotavirus strain 30/96 were also determined. The complete genome similarity between strains B4106 and 30/96 was 93.4% at the nucleotide level and 96.9% at the amino acid level. All 11 genome segments of strain B4106 were closely related to those of lapine rotaviruses and clustered with the lapine strains in phylogenetic analyses. In addition, sequence analyses of the NSP5 gene of strain B4106 revealed that the altered electrophoretic mobility of NSP5, resulting in a super-short pattern, was due to a gene rearrangement (head-to-tail partial duplication, combined with two short insertions and a deletion). Altogether, these findings confirm that a rotavirus strain with an entirely lapine genome complement was able to infect and cause severe disease in a human child.     

Dominance of Emerging G9 and G12 Genotypes and Polymorphism of VP7 and VP4 of Rotaviruses from Bhutanese Children with Severe Diarrhea Prior to the Introduction of Vaccine

A prospective study was performed to determine the molecular characteristics of rotaviruses circulating among children aged <5 years in Bhutan. Stool samples were collected from February 2010 through January 2011 from children who attended two tertiary care hospitals in the capital Thimphu and the eastern regional headquarters, Mongar. The samples positive for rotavirus was mainly comprised genotype G1, followed by G12 and G9. The VP7 and VP4 genes of all genotypes clustered mainly with those of neighboring countries, thereby indicating that they shared common ancestral strains. The VP7 gene of Bhutanese G1 strains belonged to lineage 1c, which differed from the lineages of vaccine strains. Mutations were also identified in the VP7 gene of G1 strains, which may be responsible for neutralization escape strains. Furthermore, we found that lineage 4 of P[8] genotype differed antigenically from the vaccine strains, and mutations were identified in Bhutanese strains of lineage 3. The distribution of rotavirus genotypes varies among years, therefore further research is required to determine the distribution of rotavirus strain genotypes in Bhutan.