Unusual rotavirus genotypes among children with acute diarrhea in Saudi Arabia
BMC Infectious Diseases volume 15, Article number: 192 (2015)
Human rotavirus A (human RV-A) is the most common cause of viral gastroenteritis in infants. The objective of the study was to characterize the G and P genotypes among clinical rotavirus isolates from children with acute diarrhea admitted to a tertiary care hospital in Riyadh, Saudi Arabia.
From 2011 to 2012, 541 pediatric patients with acute diarrhea were tested for rotavirus infection. RNA extractions from the fecal specimens were done by commercial kit. RT-PCR and sequencing techniques were used to detect the prevalent genotypes. Phylogenetic analysis by Maximum Likelihood method was used to study the clustering of the circulating genotypes.
The data showed that 171/541 (31.6%) faecal samples were positive for human RVA and majority were children aged below 2 years. From the G and P [types] detected it was seen that (a) 171 minus 43 ie. 128 rotavirus positives were G typed successfully (b) 171 minus 20 ie. 151 rotavirus positives were P typed successfully; (c) overall G [P] nature was determined for 113 rotavirus positives out of 171. VP4 genotyping showed that majority of the positives 146/151 (96.7%) were P ; 4/151 (2.6%) were P ; 1/151 (0.66%) was P . The dominant strains included G1P  70/113 (61.9%); G9P  19/113 (16.8%); G12P  7/113 (6.2%) and G3P  5/113 (4.4%) while the uncommon strains detected from Saudi Arabia during the study were G1P  1/113 (0.88%) and G12P  1/113 (0.88%). Phylogenetic tree, based on VP4/VP7 sequence analysis, revealed that G1P  was distinctly related to homologous strains included in human RV-A vaccine strains. Nevertheless, the uncommon genotypes G1P  and G12P  were clustered with isolates from other countries such as Bangladesh, China, Japan, Thailand and Philippines.
More studies will be required to further focus on newly emerging genotypes in our region together with the seasonality of rotavirus infection in the region, especially after January 2013 when the rotavirus vaccination has become part of routine childhood immunizations.
Human rotavirus A (human RV-A) is the most common pathogen causing viral gastroenteritis (VGE) . Human RV-A infections are responsible for 500,000 associated deaths annually among the children worldwide. In addition, human RV-A infection causes 25-35% of all gastroenteritis hospitalized cases in children . In previous work, we showed that host susceptibility to viral infections is due to genetic and environmental factors [3,4]. Infants and children are more vulnerable to human RV-A infection. More than 95% of all children under the age of 5 years were found to be infected with RV at least once; half of those were 2 years old or younger [5,6]. The incidence of human RV-A disease in industrial and developing countries is more or less the same [1,7]; however, high mortality rates mainly occur in the developing countries, particularly those with low revenues and poor health care systems .
Rotaviruses (RV) belong to Reoviridae family and are double-stranded (ds) non-enveloped RNA viruses. Their genome consists of 11 ds-RNA segments and they are classified based on 2 of these genome segments called VP4 and VP7, which encode for outer capsid proteins . Group A is the most pathogenic among common human VGE pathogens. RV-As are classified into at least 27 G and 37 P types that are known to be genetically different [10,11]. Most common types, such as G1P , G2P , G3P , and G4P  are causing more than 90% of human RV-A cases in North America and Europe . Despite these global trends observed, there are heterogeneous regional distributions found in terms of human RV-A types.
Currently, the data on the prevalence of human RV-A cases in Saudi Arabia are lacking, while the seasonal increases are observed usually in winter and spring time. Understanding the locally-circulating genotypes will eventually assist in the development of more effective vaccine candidates and prevent the spread of this infectious disease. The main objectives of the study were (1) to determine the prevalence of rotavirus genotypes with a focus on VP4 and VP7 genotypic characterization; and (2) to perform the human RV-A molecular epidemiology study in Saudi Arabia aiming at the possible emergence of new genotypes in the region.
This study was approved by the IRB Committee of King Abdullah International Medical research center (KAIMRC), National Guard Health affaires, Riyadh, KSA (reference number # RC08-112). All the samples were taken as part of standard care.
A total of 541 stool samples were collected from ongoing prospective study assessing the rotavirus epidemiology and genotypes. These isolates were taken from pediatric patients with acute diarrhea in the city of Riyadh between 2011 and 2012, which represents the main urban center in Saudi Arabia. Each stool sample was immediately frozen undiluted at −20°C until use.
Epidemiological data and data analysis
Epidemiological data including age, sex, and symptoms (diarrhoea and/or vomiting and their number of episodes), feeding type, date of onset, and date of sample collection were collected using Epic. Info vr.4 and molecular genotyping data were entered into Excel worksheets. Descriptive analysis, frequencies and percentages were calculated using SPSS vr. 20 statistical software.
ELISA screening for rotavirus
All stool samples were screened using ImmunoCard Stat Rotavirus kit® which detects the presence of rotavirus antigen in stool based on enzyme-linked immunosorbent assay (ELISA) method. Briefly, the stool specimen was diluted 1:15 in sample diluents and thoroughly mixed. The suspension was introduced (50 μl) to the sample port of the device and rotavirus detection was carried out following the manufacturer’s instructions. In addition to the internal (test) control, we included two external controls from the previously tested (known) samples.
All positive ELISA samples were confirmed by RT-PCR. In brief, RNA was extracted from approximately 1 mg of stool sample using MagNA Pure Compact RNA Isolation Kit (Roche Diagnostics GmbH, Mannheim, Germany) and following the manufacturer’s instructions. Extracted RNA was immediately stored at - 80°C until further use.
Reverse transcription (RT) and PCR amplification reactions
RT reactions were carried out in a final volume of 60 μl using Transcriptor First Strand cDNA Synthesis Kit (Roche Diagnostics GmbH, Mannheim, Germany) on a Tetrad 2 Peltier thermal cycler from BIO-RAD. A master mix was prepared containing 5× RT buffer, 40 U/μl RNAse inhibitor, 20U/μl Transcriptor reverse transcriptase, 10 mM each Deoxynucleotide Mix, 600 μM Random Hexamer Primer, RotaA.ext. Fwd primer 5′-TTT AAA ACG AAG TCT TCR ACA TGG AKG TYC TGT A-3′ and RotaA.ext. Rev primer 5′-TAA TTG GTG ATC TAC CAA TTC CTC CAG TTT G-3′ . RNA sample (10 μl) was heated at 99°C for 10 min, placed immediately on ice, and 50 μl of master mix was then added to each tube. The temperature cycling was 50°C for 50 min, 94°C for 7 min and cDNA was amplified using RotaA.ext. forward and reverse primers . The temperature cycling was 95°C for 10 min, 40 cycles of (94°C for 1 min, gradient from 45°C to 65°C for 3 min, 72°C for 1 min), and final extension at 72°C for 10 min.
Genotyping by PCR amplification and sequencing
ELISA positive rotavirus samples (n = 171) were further confirmed by RT-PCR using primers set used for rotavirus genotyping VP4 [P] and VP7 [G] genotypes. Genotyping was performed using the following oligonucleotide primer pairs (Eurofins MWG Operon) targeting the VP4 and VP7 gene regions: VP4-F 5′-TAT GCT CCA GTN AAT TGG-3′, VP4-R 5′-ATT GCA TTT CTT TCC ATA ATG-3′, VP7-F 5′-ATG TAT GGT ATT GAA TAT ACC AC-3′ and VP7-R 5′-AAC TTG CCA CCA TTT TTT CC-3′[13,14]. Thermal cycling for the VP4 region was 95°C for 10 min, 35 cycles of (94°C for 1 min, 50°C for 1 min, 72°C for 1 min), then 72°C for 10 min. Thermal cycling for the VP7 region was 95°C for 10 min, 40 cycles of (94°C for 1 min, gradient from 45°C to 65°C for 3 min, 72°C for 1 min), then 72°C for 10 min. Amplified products were analyzed on 1.5% agarose ethidium bromide-stained gels for genotyping. All PCR products of the appropriate size were purified using 3 M sodium acetate and absolute ethanol. Purified products were confirmed by sequinning using BigDye Terminator v3.1 Cycle Sequencing Kit (Applied biosystems™ Austin, TX, USA). Thermal cycling was 96°C for 1 min, 40 cycles of (96°C for 10 min, 50°C for 5 sec, 60°C for 4 min). Amplified products for VP4 (663 bp) and VP7 (880 bp) were purified a second time using XTerminator™ and SAM™ solutions (Applied Biosystems™ Foster City, CA, USA). Purified products were sequenced using 3730xl DNA Analyzer (Applied Biosystems™, Hitachi, Tokyo, Japan).
Assignment of the genotypes and phylogenetic analysis
All isolates sequences were analyzed in comparison to VP4 and VP7 of the international rotavirus sequences available at NCBI sequence database using BALSTn-2 and RotaC 2.0  automated genotyping tool for Group A rotaviruses. Molecular phylogenetic analysis was carried out using Maximum Likelihood method based on Tamura-Nei model . The tree was drawn to scale with branch lengths measured in the number of substitutions per site and was bootstrapped with 1000 replicates. The analysis involved 635 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 681 positions in the final dataset. Evolutionary analyses were conducted in MEGA6 .
Rotavirus occurrence and demographic data
The total number of cases collected between 2011 and 2012 was 541 and one stool sample was obtained from each child, representing one case. Of all 541 cases, 171 (30.6%) were found rotavirus positive by ELISA. The demographic data of the 541 cases are shown in Table 1.
Seasonal incidence of rotavirus gastroenteritis within the period from 2011 to 2012
The occurrence of acute gastroenteritis and the rotavirus incidence were found to be atypical form the global seasonality. Human RV-A infections started off with steadily increasing numbers of cases in the summer season and peaked at the fall season (around months of October/November). The number of positive cases then tapered off toward winter and after the January, it continued to decrease in the spring months to as low as 24 cases in the month of June (Figure 1).
Human RV-A infections occurred mostly in children below the age of 5 years and the highest incidence rate (69.4%) was observed in infants below the age of one year. More than 96.8% of the patients were under two years old and above 50% aged 6 months or less. The mean of the infection incidence decreased continuously regarding the patients from 3 to 5 years old.
VP4 and VP7 genotyping and molecular epidemiology
The VP4 genotyping showed that most 146/151 (96.7%) of positive samples were of genotype P , followed by genotypes P  4/151 (2.6%), and P  1 (0.66%), however, some samples n = 20/171 (11.7%) were untypeable after 2 or more attempts of amplification. The VP7 genotyping showed the following genotypes: G1 79/128 (61.7%), G2 7/128 (5.5%), G3 5/128 (3.9%), G4 5/128 (3.9%), G9 24/128 (18.8%), G12 8/128 (6.3%) and untypeable G were 43 samples out of 171 (25.1%). Overall, human RV-A genotyping and molecular epidemiology results showed that our regional RV-A combined G and P [types] were: G1P  70/113 (61.9%); G9P  19/113 (16.8%); G12P  7/113 (6.2%); G3P  5/113 (4.4%); G4P  4/113 (3.5); G2P  3/113 (2.6%); G2P  3/113 (2.6%); G1P  1/113 (0.88%); and G12P  1/113 (0.88%) (Table 2). Hence, we identified the infrequent G1P  and G12P  genotypes were observed in 2/113 1.8% of infected Saudi children as a newly emerged genotype (Table 2).
Phylogenetic trees based on the VP4 and VP7 sequence analysis revealed that G1, G2, G3, G4, G9 G12, P  and P  from Riyadh, Saudi Arabia were clustered within their related homologues strains from different countries such as: USA, Japan, China, India, Pakistan and Bangladesh with the following GenBank accession numbers (n = 80) AB118024; AB247943; AB585928; AB741657; AB905463; AF501578; AF528204; AY740736; AY787650; AY787652; D86271; DQ005112; DQ146646; DQ146652; DQ146654; DQ146667; DQ146679; DQ146684; DQ146697; DQ146699; DQ396443; DQ490539; DQ490546; DQ490551; DQ492669; DQ870505; DQ904524; EF554083; EF583046; EF672578; EF672584; EF672591; EF694178; EF694179; EF694184; FJ598026; GQ117006; GQ996883; GU985260; HF952915; HM035518; HM130969; JF490913; JF490979; JN088466; JN232048; JN258874; JN711097; JN711098; JN714986; JN827253; JN849118; JN849146; JQ069522; JQ069530; JQ087450; JQ837882; JX076836; JX195079; JX273720; JX841121; JX841124; KC443350; KC443768; KC580259; KC689824; KF414569; KF414603; KF614667; KF812586; KJ432717; KJ432745; KJ752344; KJ752410; KJ753186; KJ753373; KJ753484; KM008662; KP013536; U08431 (Figures 2 and 3).”
During the specified study period, we were able to identify certain newly emerging rotavirus genotypes introduced to Saudi Arabia, such as the uncommon G1P  and G12P  genotypes. Notably, there are at least 37 P types (P -P ) and 27 G types (G1-G27) [10,11]. The common human RV-A genotypes that are circulating worldwide include: G1P , G2P , G3P , G4P  and G9P . Globally, more than 90% of human RV-A cases have been attributed to G1P , G2P , G3P , and G4P  rotavirus genotypes. In Saudi Arabia, overall, rotavirus infection has a high prevalence (average: 33.7%) . Globally and particularly in the northern hemisphere countries, human RV-A infections occur in clusters in the winter months. In the Middle East and the gulf region several countries reported rotavirus seasonal peak between November and April including Oman, Iran, Tunisia, Morocco and Turkey . While others, Egypt and United Arab Emirates, recorded seasonal peak during July as well . In our country, however, the season of human RV-A infections starts with steadily increasing numbers of cases in July and the incidence of the disease spikes up during chilly winter months. Rotavirus infection cases seem to have increased steadily in the region and the highest numbers were observed in 2012 as compared to the previous years (Balkhy et al. Manuscript in preparation). This consistent increase in the numbers of human RV-A infections can possibly be attributed to the absence of national RV vaccination program in this country. Children less than 2 years of age were the most vulnerable to human RV-A infections. Previously, Kheyami et al. showed that 89% of rotavirus strains were of G1P8 genotype [7,21-23]. Here, we found a consistent increase in the incidence of human RV-A infections during 2012 (30.6%) and the most prevalent genotypes were G1P  (61.5%) and G9P  (16.8%). We speculate that the emergence of the uncommon genotypes (≈5%), such as G1P  and G12P , could be due to either the expression of the natural variation of distinct human RV-A genotypes over time in rotavirus-infected Saudi children or they were circulating without detection. The phylogenetic analysis showed that G1P  and G12P , isolates from Riyadh, Saudi Arabia clustered with isolates from other countries such as Bangladesh, China, Japan, Thailand and Philippines. Nevertheless, it is possible that these genotypes were merging with the immigration of the foreign workers coming from Southeast Asia as the phylogenetic analysis has revealed. The divergent incidence of human RV-A infection and the fluctuating seasonal distribution of its genotypes underscore the need for a national RV vaccination program. As though the available vaccines do not cover all uncommon genotypes, however, their efficiency against various circulating human RV-A genotypes has been demonstrated in several studies [24-26]. Of note, rotavirus vaccines were recommended worldwide by the WHO in 2008; however, these have not been a part of routine Saudi vaccination program until 2013.
Overall, the data show that 80% of our regional human RV-A genotypes included G1P  and G9P  whereas G1P  is the most prevalent (62%) rotavirus genotype in the region. In our region, the uncommon G1P  and G12P  genotypes were identified as the newly emerging RV strains. Further studies will be required to continuously monitor the seasonal prevalence and emergence of new RV genotypes in this country. In this regard, more advanced surveillance and in-depth analysis would be required, especially after the implementation of rotavirus vaccination in Saudi children. Close monitoring and regular reporting of the newly emerging genotypes would lie at the heart of a successful human RV-A infections/gastroenteritis control.
Khoury H, Ogilvie I, El Khoury AC, Duan Y, Goetghebeur MM. Burden of rotavirus gastroenteritis in the Middle Eastern and North African pediatric population. BMC Infect Dis. 2011;11:9.
Tayeb HT, Balkhy HH, Aljuhani SM, Elbanyan E, Alalola S, Alshaalan M. Increased prevalence of rotavirus among children associated gastroenteritis in Riyadh Saudi Arabia. Virol J. 2011;8:548.
Aly M, Wiltshire S, Chahrour G, Osti JC, Vidal SM. Complex genetic control of host susceptibility to coxsackievirus B3-induced myocarditis. Nat Genes Immun. 2007;8:193–204.
Vidovic S, Aly M, Flemming C, Springthorpe S, Sattar SA. First evidence of genotypes Ad3a16 and Ad3a18 in North America, obtained by genetic analysis of infectious human adenovirus from wastewaters of two urban communities in Canada. Appl Environ Microbiol. 2011;77:4256–9.
Ghazi HO, Khan MA, Telmesani AM, Idress B, Mahomed MF. Rotavirus infection in infants and young children in Makkah. Saudi Arabia J Pak Med Assoc. 2005;55:231–4.
Sanchez-Fauquier A, Gonzalez-Galan V, Arroyo S, Cabornero A, Ruiz-Burruecos A, Wilhelmi-De Cal I. Monitoring of children with acute gastroenteritis in Madrid, Spain during 2010–2011: Rotavirus genotype distribution after the vaccines introduction. Enferm Infecc Microbiol Clin. 2013;32(5):280–4.
Kheyami AM, Nakagomi OT, Dove W, Hart CA, Cunliffe NA. Molecular epidemiology of rotavirus diarrhea among children in Saudi Arabia: first detection of G9 and G12 strains. J Clin Microbiol. 2008;46:1185–91.
Telmesani AM. Oral rehydration salts, zinc supplement and rota virus vaccine in the management of childhood acute diarrhea. J Family Community Med. 2010;17:79–82.
Hoshino Y, Jones RW, Ross J, Honma S, Santos N, Gentsch JR, et al. Rotavirus serotype G9 strains belonging to VP7 gene phylogenetic sequence lineage 1 may be more suitable for serotype G9 vaccine candidates than those belonging to lineage 2 or 3. J Virol. 2004;78:7795–802.
Trojnar E, Sachsenroder J, Twardziok S, Reetz J, Otto PH, Johne R. Identification of an avian group a rotavirus containing a novel VP4 gene with a close relationship to those of mammalian rotaviruses. J Gen Virol. 2013;94:136–42.
Matthijnssens J, Ciarlet M, McDonald SM, Attoui H, Banyai K, Brister JR, et al. Uniformity of rotavirus strain nomenclature proposed by the rotavirus classification working group (RCWG). Arch Virol. 2011;156:1397–413.
Chen SC, Tan LB, Huang LM, Chen KT. Rotavirus infection and the current status of rotavirus vaccines. J Formos Med Assoc. 2012;111:183–93.
Gouvea V, Glass RI, Woods P, Taniguchi K, Clark HF, Forrester B, et al. Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens. J Clin Microbiol. 1990;28:276–82.
Ushijima H, Shinozaki T, Fang ZY, Glass RI. Group B rotavirus antibody in Japanese children. J Diarrhoeal Dis Res. 1992;10:41.
Maes P, Matthijnssens J, Rahman M, Van Ranst M. RotaC: a web-based tool for the complete genome classification of group A rotaviruses. BMC Microbiol. 2009;9:238.
Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10:512–26.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Al-Mendalawi MD. Rotavirus gastroenteritis and strain diversity in Saudi Arabia. Current status and future prospects. Saudi Med J. 2011;32:429–30.
Benhafid M, Elomari N, Elqazoui M, Meryem AI, Rguig A, Filali-Maltouf A, et al. Diversity of rotavirus strains circulating in children under 5 years of age admitted to hospital for acute gastroenteritis in Morocco, June 2006 to May 2009. J Med Virol. 2013;85:354–62.
Howidi M, Balhaj G, Yaseen H, Gopala K, Van Doorn LJ, DeAntonio R. Burden and genotyping of rotavirus disease in the United Arab Emirates: A multicenter hospital-based surveillance. Hum Vaccin Immunother. 2014;10(8):2284–9.
Kheyami AM. Rotavirus gastroenteritis and strain diversity in Saudi Arabia. Current status and future prospects. Saudi Med J. 2010;31:276–9.
Kheyami AM, Areeshi MY, Dove W, Nakagomi O, Cunliffe NA, Anthony Hart C. Characterization of rotavirus strains detected among children and adults with acute gastroenteritis in Gizan. Saudi Arabia Saudi Med J. 2008;29:90–3.
Kheyami AM, Cunliffe NA, Hart CA. Rotavirus infection in Saudi Arabia. Ann Saudi Med. 2006;26:184–91.
da Silva Soares L, de Fatima Dos Santos Guerra S, do Socorro Lima de Oliveira A, da Silva Dos Santos F, de Fatima Costa de Menezes EM, Mascarenhas JD, et al. Diversity of rotavirus strains circulating in Northern Brazil after introduction of a rotavirus vaccine: High prevalence of G3P  genotype. J Med Virol. 2013;86(6):1065–72.
Lamabadusuriya SP. Improving child survival through immunisation. Ceylon Med J. 1993;38:170–1.
Rose TL, Miagostovich MP, Leite JP. Rotavirus a genotype G1P : a novel method to distinguish wild-type strains from the Rotarix vaccine strain. Mem Inst Oswaldo Cruz. 2010;105:1068–72.
We thank King Abdullah International Medical Research Center KIMRC for funding #RC08/112, also This work has been partially supported by King Abdulaziz Medical City for Science and Technology, (KACST).
This work was approved by KAIMRC IRB committee.
The authors declare that they have no competing interests.
MA: design of the study analyzing data, performing sequence, and statistical analysis, drafting and reviewing manuscript, AK: collecting samples and data performing molecular and sequencing typing, SJ: collecting samples and data performing ELISA, drafting manuscript, HH: design of the study, data analysis, collecting data, drafting and reviewing manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Aly, M., Al Khairy, A., Al Johani, S. et al. Unusual rotavirus genotypes among children with acute diarrhea in Saudi Arabia. BMC Infect Dis 15, 192 (2015). https://doi.org/10.1186/s12879-015-0923-y