Received date: August 17, 2017; Accepted date: August 29, 2017; Published date: August 31, 2017
Citation: Intamaso U, Poomipak W, Chutoam P, Chotchuang P, Sunkkham W, et al. (2017) Genotype Distribution and Phylogenetic Analysis of Rotaviruses in Thailand and Emergence of Uncommon Genotypes. Arch Clin Microbiol. Vol. 8 No.4:60. doi: 10.4172/1989-8436.100060
Copyright: © 2017 Intamaso U, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Rotavirus remains a major cause of diarrhea in infants and young children especially in developing countries. Monitoring changes in rotavirus strains is necessary to assess the potential effectiveness of vaccines in specific geographic locations. This project aimed at monitoring the epidemiological status of rotavirus among Thai children. A total of 73 positive fecal samples for rotavirus were collected from hospitalized infants and children with acute gastroenteritis or diarrhea between 2012 and 2014 and analyzed for G- and P-genotypes by semi-nested multiplex PCR or DNA sequencing. The most prevalent genotype was G1P (61.64%) followed by G8P (21.92%), G2P (8.22%), G3P (4.11%), G9P (2.74%), and G2P (1.37%), Interestingly, the unusual genotypes, G8P and G2P were identified in this study. A phylogenetic analysis of VP7 and VP4 genes of strains detected in the present study revealed that G1 belonged to G1-Ic, and G1-IIc ; G2 belonged to G2-IIa; G3 belonged to G3-Ia; G8 belonged to G8-I; and G9 belonged to G9- III lineages, while P and P were identified as P-Vb and P-III lineages, respectively. The distribution of unusual genotypes that possess neither VP7 nor VP4 specificity with the available rotavirus vaccine currently in use may represent a challenge to the outcome and success of vaccination.
Rotavirus; Prevalence; Genotype; Phylogenetic tree; Thailand
Rotavirus group A (RVA) is one of the most common etiologies of severe diarrhea in infants and young children worldwide. It has been estimated that children under 5 years old are inflicted with gastroenteritis that results in 527,000 fatalities each year, with the highest rate in developing countries in Africa and Asia . In Thailand, between 1997 and 2005, on average 42.6% of pediatric patients admitted to hospitals had severe dehydration due to diarrhea caused by rotavirus infection .
Rotavirus is classified as a member of genus Rotavirus in the Reoviridae family. These viruses are approximately 75 nm icosahedral non-enveloped particles, consisting of 11 segmented double stranded RNA genome enclosed in a triple-layered capsid. The RVA gene segments are monocistronic and encode for six structural (VP1-4, 6 and 7) and four non-structural proteins (NSP1- 4), except for the gene segment 11 which has two overlapping ORFs and codes for NSP5 and NSP6. Besides classification of rotavirus in A-G groups based on the antigenicity of the middle layer VP6, the two outer layer proteins VP7 and VP4 also provide the basis for the dual classification system into G and P genotypes, respectively . To date 27 G-genotypes and 37 P-genotypes have been reported in humans and animals . Globally, G1, G2, G3, G4, and G9 in combination with P, P, and P are the most common rotavirus genotypes. Furthermore, some unusual genotypes of rotaviruses have been reported in several countries [5-7]. Thus the dual classification system is very useful to categorize the viruses in G/P combinations in which up to 120 genotypes could be found . However, only a limited number of the combined G/P genotypes are common in humans, such as G1P, G2P, G3P, G4P and G9P with G1P, which are the prevalent types of human rotaviruses detected in most surveys. Other genotypes such as G5, G6, G8, G10, and G12 may also be co-circulating with P, P, P, and P in various geographic regions of the world.
In Thailand, the fluctuation of predominant genotypes has also been detected since 1997 . Furthermore, the uncommon human rotavirus strains such as G3P, G3P, G3P, G3P, and G12P have been detected in some areas of Thailand [2,10]. Two rotavirus vaccines, Rotarix® (GlaxoSmithKlinne Biologicals), a monovalent G1P1A human attenuated vaccine, and RotaTeq® (Merck&Co., Inc.), a pentavalent human bovine reassorted vaccine expressing human G1, G2, G3, G4, and P1A have already been licensed in many countries. The efficacy of rotavirus vaccines vary among different populations. Following vaccination, the rates of mortality and hospitalization tend to reduce especially in developing countries [11,12] and even lower in populations of the developing countries with the highest burden of the disease . These two licensed rotavirus vaccines were introduced in Thailand in 2006 and 2008, respectively, but are not included as a part of the national immunization program, thereby limiting their access for health care . Because of great genetic diversity and interspecies transmission, the genetic reassortments of rotavirus strains are common. Thus, continuous rotavirus surveillance is important to monitor the new types or new G/P combinations to evaluate whether current vaccines cover the most prevalent genotypes circulating in Thailand.
In this study, we aimed to determine the distribution pattern of group A rotavirus genotypes and perform a phylogenetic analysis of rotaviruses circulating in infants and children with acute diarrhea who were hospitalized in Bamrasnaradura Infectious Diseases Institute, Department of Disease Control, Ministry of Public Health, Nonthaburi, Thailand.
Specimen collection and rotavirus detection
Children under age of five with moderate to severe diarrhea who were admitted to Bamrasnaradura Infectious Diseases Institute, Nonthaburi, Thailand, were screened for rotavirus infection by using a rapid qualitative immunochromatography assays (VIKIA Rota-Adeno, bioMerieux SA, France). 73 samples that were positive for rotavirus were collected anonymously from April 2012 to April 2014 and kept at -80°C for further examination. The study protocol was approved by the Institutional Ethical Committee of Bamrasnaradura Infectious Diseases Institute and Burapha University.
Viral RNA extraction
Fecal samples were diluted 10% in phosphate buffered saline (PBS) and clarified by centrifugation at 3500 g for 10 min. The supernatant fluid was collected and 140 μl of each was processed for viral RNA extraction with QIAamp® Viral RNA Mini Kit (QIAGEN, Germany).
Detection of group A rotavirus by RT-PCR
Partial-length VP7 (881 bp) and the VP8* region (876 bp) of VP4 genes were amplified separately by one step RT-PCR (SuperScript III one-step RT-PCR with Platinum Taq, USA). The reaction mixture consisted of 3 μl of extracted RNA, 12.5 μl of 2x reaction mix (Invitrogen, USA), SuperScriptIII RT/platinum Taq mix, 0.4 μM of each primer, 2.0 mM of MgSO4 in the total of 25 μl. A reaction mixture without a template was used as a negative control, and the reference strain, human rotavirus A strain Wa (ATCC 2018) RNA was used as a positive control. The VP7 gene was amplified with a primer pair VP7F and VP7R  to obtain a partial 881 bpgene fragment. The VP4 gene was amplified with a primer pair VP4 forward con3 and VP4 reverse con2 primer  to generate an 876 bp-partial gene fragment.
G and P genotyping by multiplex semi-nested PCR
The genotyping of group A rotaviruses was characterized by the multiplex semi-nested PCR reaction with a pool of internal type-specific forward primers and the reverse primer. One microliter of VP7 or VP4-gene product was used as the template. RNA of human rotavirus A strain Wa (ATCC 2018) and reaction without template were used as the positive and negative control, respectively.
To characterize G genotype, a mixture of forward primers, aBT1, aCT2, mG3, and mG9 in combination with the VP7 reverse primer  were utilized for amplification of the VP7 G1, G2, G3, and G9 gene fragment, respectively, in the PCR reaction. The reaction was carried out in a total of 25 μl consisting of 0.2 mM dNTPs mix, 2.0 mM MgCl2, 1.0 μM each of forward and reverse primer, 2.5 U Taq DNA polymerase (Thermoscientific, USA). Genotyping was based on size comparison of PCR products after direct visualization on 1.5% agarose gel. Expected DNA fragments of 618 (G1), 521 (G2), 682 (G3), and 179 (G9) bp were detected compared with 100 bp DNA Ladder marker (TrackIt, Invitrogen, USA).
For P genotyping, a mixture of 0.5 μM of each 2T-1, 3T-1, IT-1, and 4T-1 forward primers were used in combination with the VP4 reverse con2 primer  for amplification of the VP4 P, P, P, and P gene fragment, respectively in the PCR reaction containing 0.2 mM dNTPs mix, 2.5 mM MgCl2, and 2.5 U Taq DNA polymerase (Thermoscientific, USA). Expected amplified products of 413 bp, 628 bp, 531 bp, and 502 bp DNA fragments represent P, P, P, and P respectively visualized on 1.5% agarose gel.
Both G and P genotyping were performed under the same thermocycler condition as follows: initial denaturation at 95°C, 3 min; denaturation at 95°C, 30 sec; annealing at 55°C, 30 sec, extension at 72°C, 1 min in the total 30 cycles and followed by the final extension at 72°C, 10 min. The viruses with genotypes that could not be determined by multiplex PCR were subjected to identification by nucleotide sequencing.
Nucleotide sequencing and phylogenetic analysis
Partial-length cDNA were extracted from the agarose gel and purified with PureLink™ Quick Extraction kit (Invitrogen, USA). The nucleotide sequence of each DNA fragment of VP7 and VP4 gene was determined in both directions with the specific oligonucleotide sequencing primers; VP7-F primer or VP7-R primer for VP7 gene; VP4-F (con3) or VP4-R (con2) primers for VP4 gene (First Base Laboratories SDN BHD, Malaysia). The nucleotide sequences of the VP7 and VP4 genes were assembled and analyzed using BioEdit. The sequences were subsequently aligned and compared with those of the reference strains available in the NCBI GenBank database using BLAST search (www.ncbi.nlm.nih.gov/blast) and analyzed using CLUSTAL X program, version 1.83, for multiple sequence alignment, and MEGA4.1 for phylogenetic tree construction. Lineage designation was based on similarity to previously defined lineage reference strain. The statistical significance of the constructed phylogenies was analyzed using neighbor-joining method with the 1000 replicates. Branches with bootstrap number below 50% were considered insignificant and excluded.
Nucleotide sequence accession number
Prototype sequences used in the phylogenetic analyses were retrieved from the GenBank database based on the closest nucleotide sequence relatedness with the strains sequenced in the present study. Prototypes G and P sequences of most common genotypes, lineages, and sub-lineages in human were also retrieved from GenBank and included in the phylogenetic analysis.
The nucleotide sequences of group A rotavirus described in the present study have been deposited in the GenBank database with the access numbers as follows: VP7 G1(KU588886- KU588899, KU588901-KU588902, KU588908, KU588912- KU588915, KU588919, KU588927-KU588928, KU588932- KU588936, KU588939,KU588942-KU588956), VP7 G2(KU588878- KU588884), VP7 G3(KU588859- KU588861) VP7 G8(KU588862- KU588877), and VP7G9 (KU991795-KU991796); VP4P (KU754055-KU754060), and VP4P (KU754064-KU754078, KU754080- KU754084, KU754087- KU754093, KU754098- KU754102, KU754106- KU754107, KU754109- KU754112, KU754116- KU754126, KU754128- KU754132, KU754134- KU754135, KU754137-KU754145).
Distribution of group A rotavirus strains
The G- and P-genotype of 73 strains of group A rotaviruses were identified by RT-multiplex nested-PCR and nucleotide sequencing. Five different G-genotypes, i.e., G1, G2, G3, G8, and G9 were detected during this study period (Table 1). Among the rotavirus strains detected, G1 is the most predominant genotype (61.64%; N=45/73), followed by G8 (21.92%; N=16/73), G2 (9.59%; N=7/73), G3 (4.11%; N=3/73), and G9 (2.74%; N=2/73). G4, a common genotype, was not found in any rotavirus positive samples by nucleotide sequence analysis of partial length VP7 gene. When examined for their P-genotypes, only P and P were detected. The majority of rotavirus isolates were P 91.78% (N=67/73) whereas only a small proportion were P 8.22%. (N=6/73). All of G1, G3, G8, and G9 rotavirus strains were found exclusively with P, while G2 strains were found exclusively with P. However, only one G2 isolate coexisted with P. The G- and P-genotype combinations shown in Table 1 indicate that G1P was the predominant genotype (61.64%; 45/73), followed by G8P (21.92%; 16/73), G2P (8.22%; 6/73), G3P (4.11%; 3/73), G9P (2.74%; 2/73), G2P (1.37%; 1/73).
|P genotypes||G Genotypes|
|2 (2.74%)||67 (91.78%)|
Table 1: Distribution and relative frequencies of various combinations of G- and P-Genotypes GroupA rotavirus detected in hospitalized children with diarrhea.
Phylogenetic analysis of G-type rotaviruses
In order to gain insight into the degree of variability in the rotavirus strains currently circulating in Thailand, nucleotide sequences of rotaviruses obtained from 73 fecal samples were compared with the reference strains reported worldwide. Based on the phylogenetic analysis created for VP7 gene (Figure 1). G1 strains were classified into seven lineages (I-VII). Most G1 strains belonged to lineage II, sub-lineage IIc (Figure 1a). Those strains showed more than 93% similarity between themselves and their genes are closely related (94-99%) to strains previously identified in Japan, Australia, Cuba, the USA, Italy, India, Russia, and Thailand. G1 strains which were clustered into lineage I fell into sub-lineage Ic were segregated into 3 distinct branches, of which the major branch consisted of strains that shared 99% nucleotide similarity with other Thai strains previously reported (GQ996867CU163-KK/07 and GQ 996868 CU246-KK/08). The other three strains (KU588889Ash B-1014, KU588897Ash B-1034, and KU588893Ash B-1019) shared 96-97% identity with previously identified strains (DQ512969 VN-368 Vietnam, DQ207389 6193- 03 Ireland and AY098670 ISO-4 India),while another strain (KU588887 AshB-1011) in a separate branch shared 97-99% similarity with previously identified strains (AF183848 CH9 Taiwan, U26364 Ban-48 Bangladesh, GU985249 PJ13/07 China, GU979204 GER169-09 Germany, AB585918 8735/05 Hong Kong, and HM467952 LB2719 USA).
Figure 1: Phylogenetic analysis of VP7 nucleotide sequences of representative G genotype rotaviruses detected in the study. (a)G1 (b)G2 (c)G8. The phylogenetic tree was constructed based on the neighbor-joining method with the 1000 replicates using the MEGA 5.1 program. Percent bootstrap support is indicated by number at each branch node. The lineage and sub-lineage are indicated on the right side. The representative strains from the present study are indicated with a filled dark circle. Sequences of reference strain were obtained from the GenBank and EMBL databases and named with accession numbers, countries or animals of origin. The bar at the bottom of the trees indicates variation scale.
G2 could be classified into four lineages (I-IV). The VP7 nucleotide sequences of seven strains (KU588883 Ash-B1022, KU588884 Ash-B1027, KU588879 Ash-B1002, KU588881 Ash-B1020, KU588882 Ash-B1021, KU588878 Ash-B1001, and KU588880 Ash-B1012) of G2 belonged to lineage II, sub-lineage IIa and showed 100% similarity to each other (Figure 1b). All strains showed that their genetic backgrounds were closely related to the G2 strains previously identified in the USA, Australian, Uruguay, India, Kenya, Taiwan Japan, Vietnam, China, and Thailand from 2001 to 2003  in the range of 96-99% similarity.
G3 could be classified into 2 lineages (I-II). All three strains of G3 (KU588859 Ash-B2022, KU588860 Ash-B2024, and KU588861 Ash-B2037) belonged to lineage I, sub-lineage Ia, and show 100% identity among themselves. These strains displayed more closely related genetic background to the G3 strains previously identified worldwide, such as in Thailand, China, Hong Kong, Vietnam, South Africa, Argentina, Italy, Germany, and the USA.
Interestingly, very rare G8 strains were also detected in this period of study at a relatively high frequency (21.92%). All 16 strains belonged to lineage I and showed VP7 nucleotide sequence with over 99 % similarity (Figure 1c). Their genetic backgrounds were more closely related to G8 bovine strains previously identified in India (96-98%).
In this period of study only two G9 strains (KU991795 and KU991796) were detected and were classified into lineage III with 99% nucleotide similarity among themselves and 91- 93% similarity with G9 strains previously identified in Belgium, Canada, Japan, China, India, and Thailand (AB436829 CHB058-02 Thailand and AB436825 MS057-02).
Phylogenetic analysis of P-types rotaviruses
Based on VP 4 gene, the total six strains of P were detected in the present study (Figure 2a). All of these strains shared 100% sequence identity with each other, and were more closely related to strains detected in Myanmar with 99% similarity. Based on the analysis phylogenetic tree, all P strains belonged to lineage V, sub-lineage Vb.
Figure 2: Phylogenetic analysis of VP4 nucleotide sequences of P genotype rotaviruses detected in the study: (a) P (b) P. The tree was constructed based on the neighbor-joining method with the 1000 replicates using the MEGA 5.1 program. Percent bootstrap support is indicated by a number at each branch node. The lineage and sub-lineage are indicated on the right side. The representative strains from the present study are indicated with a filled dark circle. Reference strain sequences were obtained from the GenBank and EMBL databases and named with accession numbers and countries or animals of origin. The bar at the bottom of the trees indicates variation scale.
A total of 67 P strains identified in the present study belonged to lineage III and clustered in distinct branches (Figure 2b). These strains shared 96-100% nucleotide sequence identity to P strains previously identified in Thailand.
Rotavirus infection remains an important cause of mortality among young children in the developing countries . Since natural changes in the prevalence of rotavirus strains are common, continuous monitoring of the circulating genotypes in the community is essential to ensure the delivery of efficacious interventions.
When compared to previous studies  from 1993 to 2011, the genotypes of rotavirus prevalence in Thailand have changed over time. The G1 genotype, a common strain worldwide , was still found with high frequency in the present study. In Thailand, G1 expressed a high incidence rate during 1993 to 1999, which subsequently decreased. From 2003 to 2004, the G1 genotype reached its peak and remains to be the predominant strain in Thailand. Sequence analysis of VP7 gene revealed that most of the G1 genotype detected belonged to the lineage II, sub-lineage IIc, and the lineage I, sub-lineage Ic which were closely related to the G1 genotype reported in other regions of Thailand (from 2004 to 2006 and 2007 to 2008) and the strains discovered worldwide. These findings indicated the genetic diversity of G1 and its wide distribution. In addition, G1 strains were observed exclusively in combination with P that shared identical nucleotide sequence with P strains previously identified in most parts of Thailand [22,23].
G2 has been found to be common worldwide. In Thailand the prevalence of G2 rotavirus has also changed over time. During 1993 to 2000 the prevalent rate of G2 in many regions of Thailand was much lower than those of G1[22, 23]. However, from 2002 to 2004, G2 became the dominant genotype in Bangkok  and co-dominant with G9 in other areas of the country . The prevalence rate of G2 showed a fluctuating pattern from1999 to 2004 , then continuously decreased. During 2008 to 2009 the prevalence rate of G2 became co-predominant with G1 again  and then decreased to a very low frequency or was undetectable during 2009 to 2011 . In the present study, G2P was the third most predominant genotype after G1P and G8P. Based on the phylogenetic analysis, all strains belonged to lineage G2, sub-lineage IIa and seemed to be a single strain due to a 100% similarity to each other. A genetic distance analysis also showed that these G2 strains were closely related to G2 strains which were prevalent worldwide. This finding implies that G2 rotavirus has been prevailing and spreading throughout the world . In addition, G2 have been detected exclusively in combination with P genotype. Interestingly, only one case of G2 in combination with P was detected in this study, which was also reported in a previous study in 2006 to 2007 . It suggested that G2P found in this surveillance might re-emerge from the original strains in the area. It is, therefore, important to monitor the trend and evolution of G2 in the next epidemic seasons.
G3 in Thailand have been reported at a far fewer frequency compared to G1 and G2  or even undetectable in many areas of the country . In the present study, the G3 genotype was detected at a low frequency which agreed with the long-term survey in Thailand. However, this trend appeared to be opposite to some areas of the country that showed a prevalence of the G3 genotype in KhonKaen and Bangkok . The phylogenetic analysis revealed that the detected G3 strain belonged to lineage I, sub lineage Ia and is likely to be a single strain which showed a closer genetic relationship with G3 strains previously identified worldwide as well as in Bangkok, Thailand (DQ67493464vp7n) . It suggested that G3P found in this surveillance might have originated from the strains circulating in Bangkok, Thailand.
G9 was first reported in Bangkok at a very low frequency during 1994 to 1995 then disappeared afterwards . During 2000 to 2002, it re-emerged and became the dominant genotype . Subsequently, its prevalent pattern tended to fluctuate or was even undetectable until 2011 [20,25]. In the present study, G9 was found to be of lesser frequency than G3. However, in this study G9 strains from two positive samples were found in exclusive combination with P. The strain was classified into lineage III whose members are spreading worldwide. Because of the close genetic relationship, the G9 rotavirus detected in the present study is likely derived from the G9 previously identified in Thailand . In addition, the nucleotide sequence of their VP7 gene was also closely related to the strains detected in Belgium, Canada, Japan, China, and India, indicating that they might have descended from the same ancestor. Unlike other globally common G1, G2, G3, or G4 strains, which are detected almost exclusively with P or P, the G9 strains have been detected in various P types .
A major finding in this study was that the globally uncommon genotype G8 in association with P was detected at a relatively high frequency. It should be noted that the “Gouvea” primers showed cross-reactivity between types G3 and G8  and DNA sequencing was necessary to confirm the VP7G8 genotyping results. The G8 rotavirus genotype is typical of strains infecting bovine and other animal host species, normally in combination with P, P, and P types. In humans, the first G8 rotavirus was detected in an Indonesian child in 1980 , afterwards its identification in children with diarrhea has been reported sporadically worldwide especially in Africa , also in Sao Tome and Principe, Portugal, , and Italian regions . In Thailand, G8 was first detected at a low percentage from 2003 to 2004 . All G8 strains detected in the present study were closely related to bovine G8 strains previously identified in India (96- 98%) and sheep G8 strains previously identified in Spain (86%). According to the new RV-A classification system, the intragenotype nucleotide identity which utilizes a cutoff value of 80% for the VP7 gene , these two animal G8 strains effectively belong to the same genotype. It is possible that bovines were the main natural reservoir and the source of human rotavirus G8 strains identified in this study. Its appearance in human is likely to be the result of zoonotic transmission in the past. Interspecies transmission between humans and animals or animals and animals might have taken place in nature . These interspecies and zoonotic transmissions and close proximity to humans with domestic animals could result in the reassortment of new rotavirus strains between human P strains bearing G1, G3,G4, or G9 VP7 specificity and cocirculating bovine G8 strains bearing whether the P or P allele [30,34]. In this study the G8P seems almost identical thus it probably excludes putative “epidemic outbreaks”.
However, analysis of the VP7 and VP4 genes may not be sufficient to provide a conclusive date on the possible origin of a rotavirus strain. The complete genome sequence information is necessary to obtain a better understating of the possibility of human-animal reassortments, cross-species transmission. It may also provide evidence of animal rotaviruses acting as reservoirs for one or several genes of human rotavirus strains. Since animal rotavirus or human-animal reassortants circulate in the population and in the environment, the monitoring of new emerging virus strains may help anticipate possible reduced efficacy of current vaccines to protect children against novel unconventional genotypes. Although both Rotarix™ and RotaTeq™ have shown high efficacy against the five major genotypes circulating worldwide and broad cross-protection against serotypes not included in the vaccine composition, efficacy data against genotypes of non-human origin remains limited .
G1P genotype, the most common worldwide, was the most prevalent genotype in this study. The uncommon G8P and G2P genotypes were also identified. This present study has demonstrated the local fluctuation of genotype prevalence that frequently occurs over time, and emphasized the importance of continuous surveillance of rotavirus genotypes circulating in human and other animal species in order to obtain a better understanding of the epidemiology of rotavirus A infection within the region. More extensive epidemiological studies are also required to assess the potential effectiveness of current vaccines in our geographic location.
This work was supported by the office of the National Research Council of Thailand. We would like to thank Prof. Prasert Sobhon and Prof. Somsak Pantuwatana for advice, encouragement, and editing the manuscript.
All Published work is licensed under a Creative Commons Attribution 4.0 International License
Copyright © 2019 All rights reserved. iMedPub LTD Last revised : February 21, 2019