Rotavirus group A genotype circulation patterns across Kenya before and after nationwide vaccine introduction, 2010-2018

Rotavirus genotype dynamics in Pakistan: G9 and G12 emerging as dominant strains in vaccinated children (2019)

Rotavirus A (RVA) is a leading cause of severe gastroenteritis in children worldwide, and vaccination has become a pivotal strategy to reduce the associated morbidity and mortality. This study presents a molecular characterization of RVA genotypes circulating among vaccinated children in Pakistan during the year 2019. A total of 510 stool samples were collected from children of up to five years of age presenting with acute gastroenteritis symptoms in Rawalpindi, Islamabad regions of Pakistan. The RVA antigen was detected using ELISA on these samples. RVA G/P genotyping was performed on ELISA positive samples using Multiplex semi-nested reverse transcriptase PCR. RVA was found in 130 fecal samples, with an overall prevalence of 25.4 %. G9P[8] (20 %) is the most prevalent genotype, followed by G12P[6] (17 %), G3P[8] (14 %), G1P[8] (12 %), G2P[4] (10 %), G12P[8] (7 %), G9P[6] (7 %), G3P[6] (6 %), G3P[4] (4 %) and G1P[6] (3 %) respectively. There is a statistically significant difference (p < 0.05) found in the group age (in months) of RVA gastroenteritis cases as detected by RT-PCR. The highest number of positive cases was found in the age range from 0 to 6 months, followed by 7-12 months, 13-24 months, and 25-60 months, respectively. Dehydration is statistically significantly associated (p˂ 0.05) in RVA gastroenteritis cases compared to those who tested negative. This study emphasizes the significance of maintaining a continuous surveillance system and conducting genomic analysis of RVA genotypes in children upto the age of 5 years. This is essential for tracking the circulation of RVA genotypes. The results from this research enhance our comprehension of how RVA genotypes are changing over time in Pakistan, underscoring the ongoing necessity for improving vaccine coverage and effectiveness. This, in turn, can help reduce the impact of RVA-related illnesses in children.

Long-term impact of rotavirus vaccination on all-cause and rotavirus-specific gastroenteritis and strain distribution in Central Kenya: An 11-year interrupted time-series analysis

BACKGROUND

Kenya introduced a monovalent rotavirus vaccine administered orally at 6 and 10 weeks of age into her National Immunization Program in July 2014. The study evaluated the long-term impact of the vaccine on hospitalization for all-cause and rotavirus-specific acute gastroenteritis (AGE) and strain epidemiology in Kenya.

METHODS

Data on all-cause and rotavirus-specific AGE and strain distribution were derived from an eleven-year hospital-based surveillance of AGE among children aged <5 years at Kiambu County Teaching and Referral Hospital (KCTRH) in Central Kenya between 2009 and 2020. Fecal samples were screened for group A rotavirus using ELISA and genotyped using multiplex semi-nested RT-PCR. Trends in all-cause and rotavirus-related AGE and strain distribution were compared between the pre-vaccine (July 2009-June 2014), early post-vaccine (July 2014-June 2016) and late post-vaccine (February 2019-October 2020) periods.

RESULTS

Rotavirus-specific AGE was detected at 27.5% (429/1546, 95% CI: 25.5-30.1%) in the pre-vaccine period; 13.8% (91/658, 95% CI: 11.3-16.6%) in the early post-vaccine period (July 2014-June 2016); and 12.0% (229/1916, 95% CI: 10.6-13.5%) in the late post-vaccine period (February 2019-October 2020). This amounted to a decline of 49.8% (95% CI: 34.6%-63.7%) in rotavirus-specific AGE in the early post-vaccine period and 53.4% (95% CI: 41.5-70.3%) in the late post-vaccine period when compared to the pre-vaccine period. All-cause AGE hospitalizations declined by 40.2% (95% CI: 30.8%-50.2%) and 75.3% (95% CI: 65.9-83.1%) in the early post-vaccine and late post-vaccine periods, respectively, when compared to the pre-vaccine period. G3P [8] was the predominant strain in the late post-vaccine period, replacing G1P[8] which had predominated in the pre-vaccine and early post-vaccine periods. Additionally, we detected considerable proportions of uncommon strains G3P[6] (4.8%) and G12P[6] (3.5%) in the post-vaccine era.

CONCLUSION

Rotavirus vaccination has resulted in a significant decline in all-cause and rotavirus-specific AGE, and thus, provides strong evidence for public health policy makers in Kenya to support the sustained use of the rotavirus vaccine in routine immunization. However, the shift in strain dominance and age distribution of rotavirus AGE in the post-vaccine era underscores the need for continued surveillance to assess any possible vaccine-induced selective pressure that could diminish the vaccine effectiveness over time.

Phylodynamic characteristics of reassortant DS-1-like G3P[8]-strains of rotavirus type A isolated in Nizhny Novgorod (Russia)

Since 2013, there has been an increase in reports of the spread of a double intergroup reassortant strain of rotavirus type A (RVA) with the genotype G3P[8] and other genes belonging to the second genogroup I2-R2-C2-M2-A2-N2-T2-E2-H2. In our study, we provide a molecular genetic characterization of rotaviruses with genotype G3P[8]-I2 isolated in Nizhny Novgorod. In our study, we used RT-PCR, Sanger sequencing, RNA-PAGE methods. Phylogenetic and phylodynamic analysis were performed using the Bayesian approach. According to our study, there was a significant increase in the proportion of G3P[8] from 15% during the period of 2020-2021 to 53% during the period of 2021-2022 in Nizhny Novgorod, Russia. Phylogenetic analysis based on the VP4 gene revealed that DS-1-like RVAs isolated in Nizhny Novgorod belong to different clusters of the P[8]-3.1 lineage, with a level of variation ranging from 1.1% to 1.3%. Based on the VP6 gene, the equine-like RVAs identified by us carry genetic variants belonging to three distinct clusters of the lineage I2-V, with a variation level ranging from 2.0% to 4.5%. These data indicate the genotypic diversity of circulating DS-1-like G3 RVAs. Phylogenetic analysis of the VP7 gene allowed us to assign the isolates identified in our study to the G3-1 lineage. We estimated that the circulation of the most recent common ancestor of the spreading strains dates back to 2002. Additionally, we determined the typical level of mutations in the VP7 gene, which amounted to 2.14*10-3 substitutions/per site/per year.

Classification of group A rotavirus VP7 and VP4 genotypes using random forest

Introduction: Group A rotaviruses are major pathogens in causing severe diarrhea in young children and neonates of many different species of animals worldwide and group A rotavirus sequence data are becoming increasingly available over time. Different methods exist that allow for rotavirus genotyping, but machine learning methods have yet to be explored. Usage of machine learning algorithms such as random forest alongside alignment-based methodology may allow for both efficient and accurate classification of circulating rotavirus genotypes through the dual classification system. Methods: Random forest models were trained on positional features obtained from pairwise and multiple sequence alignment and cross-validated using methods of repeated 10-fold cross-validation thrice and leave one- out cross validation. Models were then validated on unseen data from the testing datasets to observe real-world performance. Results: All models were found to perform strongly in classification of VP7 and VP4 genotypes with high overall accuracy and kappa values during model training (0.975-0.992, 0.970-0.989) and during model testing (0.972-0.996, 0.969-0.996), respectively. Models trained on multiple sequence alignment generally had slightly higher overall accuracy and kappa values than models trained on pairwise sequence alignment method. In contrast, pairwise sequence alignment models were found to be generally faster than multiple sequence alignment models in computational speed when models do not need to be retrained. Models that used repeated 10-fold cross-validation thrice were also found to be much faster in model computational speed than models that used leave-one-out cross validation, with no noticeable difference in overall accuracy and kappa values between the cross-validation methods. Discussion: Overall, random forest models showed strong performance in the classification of both group A rotavirus VP7 and VP4 genotypes. Application of these models as classifiers will allow for rapid and accurate classification of the increasing amounts of rotavirus sequence data that are becoming available.

Genomic epidemiology of the rotavirus G2P[4] strains in coastal Kenya pre- and post-rotavirus vaccine introduction, 2012-8

The introduction of rotavirus vaccines into the national immunization programme in many countries has led to a decline in childhood diarrhoea disease burden. Coincidentally, the incidence of some rotavirus group A (RVA) genotypes has increased, which may result from non-vaccine-type replacement. Here, we investigate the evolutionary genomics of rotavirus G2P[4] which has shown an increase in countries that introduced the monovalent Rotarix® vaccine. We examined sixty-three RVA G2P[4] strains sampled from children (aged below 13 years) admitted to Kilifi County Hospital, coastal Kenya, pre- (2012 to June 2014) and post-(July 2014 to 2018) rotavirus vaccine introduction. All the sixty-three genome sequences showed a typical DS-1-like genome constellation (G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2). Pre-vaccine G2 sequences predominantly classified as sub-lineage IVa-3 and co-circulated with low numbers of sub-lineage IVa-1 strains, whereas post-vaccine G2 sequences mainly classified into sub-lineage IVa-3. In addition, in the pre-vaccine period, P[4] sub-lineage IVa strains co-circulated with low numbers of P[4] lineage II strains, but P[4] sub-lineage IVa strains predominated in the post-vaccine period. On the global phylogeny, the Kenyan pre- and post-vaccine G2P[4] strains clustered separately, suggesting that different virus populations circulated in the two periods. However, the strains from both periods exhibited conserved amino acid changes in the known antigenic epitopes, suggesting that replacement of the predominant G2P[4] cluster was unlikely a result of immune escape. Our findings demonstrate that the pre- and post-vaccine G2P[4] strains circulating in Kilifi, coastal Kenya, differed genetically but likely were antigenically similar. This information informs the discussion on the consequences of rotavirus vaccination on rotavirus diversity.

M. Eljaro Z, Abu-Awwad F, Bozya Y, Hjaija D, Bhat N, Leader T, Ramlawi A, Marzouqa H. Introducing ROTAVAC® to the occupied Palestinian Territories: Impact on diarrhea incidence, rotavirus prevalence and genotype composition

BACKGROUND

Rotavirus infection remains an important cause of morbidity and mortality in children. The introduction of vaccination programs in more than 100 countries has contributed to a decrease in hospitalizations and mortality. This study investigates the epidemiological impact of the rotavirus vaccine ROTAVAC® in the Palestinian Territories, the first country to switch from ROTARIX® to this new vaccine.

METHODS

Clinical surveillance data was collected fromchildren younger than 5attendingoutpatient clinics throughout Gaza withdiarrhea between 2015 and 2020. The incidence of all-cause diarrhea was assessed using an interrupted time-series approach. Rotavirus prevalence was determined at the Caritas Baby Hospital in the West Bank usingELISA on stool specimen of children younger than 5with diarrhea. Genotyping was performed on 325 randomly selected rotavirus-positive samples from January 2015 through December 2020 using multiplex PCR analysis.

RESULTS

Average monthly diarrhea casesdropped by 16.7% annually fromintroduction of rotavirus vaccination in May 2016 to the beginning of the SARS-CoV-2 epidemic in March 2020 for a total of 53%. Case count declines were maintained afterthe switchto ROTAVAC® in October 2018. Rotavirus positivity in stool samples declined by 67.1% over the same period without change followingthe switch to ROTAVAC®. The distribution of predominant genotypes in rotavirus-positive stool samples changed from a pre-vaccination G1P [8] to G9P[8] and G12P[8] during the ROTARIX® period and G2P[4] after the introduction of ROTAVAC®.

CONCLUSION

ROTAVAC® has shown epidemiological impact on par with ROTARIX® after its introduction to the national immunization schedule in the Palestinian Territories. A molecular genotype shift from a pre-vaccination predominance of G1P[8] to a current predominance of G2P[4] requires more long-term surveillance.

Circulating rotavirus P[8]-lineage IV, unlike P[8]-lineage III, significantly related to nonsecretors status in Iranian children

Rotavirus (RV) P[8] strains are responsible for the most of the RV infections globally and are significantly associated with the secretor and Lewis positive status. Among the distinct P[8] lineages, different ligand affinities have been detected which can be linked to differences in secretor status associated histo-blood group antigens (HBGAs). Herein, we report the lineages of P[8] strains and their associated secretor and Lewis antigen phenotypes in Iranian children. The phylogenetic tree and sequence analyses showed that the most common detected RV P[8] strain belonged to P[8]-lineage III (92%) and were significantly associated with secretor and Lewis positive status. In contrast, 8% of P[8] strains clustered into the P[8]-lineage IV and were significantly associated with nonsecretor status, implying that lineage IV tends to infect nonsecretor individuals. Furthermore, protein modeling and amino acid analyses of the VP8* glycan binding site of Iranian P[8]-lineage IV strains indicated two residual substitutions (T184V and N216V/I) compared to the P[8]-lineage III strains that might have affected the glycan affinity among P[8]-lineages IV strains. The corresponding residual changes might permit their continued transmission in nonsecretor children in competition with other P[8]-lineages. Although nonsecretors show natural resistant to P[8] strains, but such residual changes might overcome this natural resistance which in turn might indirectly contribute to the decline in the vaccine efficacy in populations where HBGA polymorphism allows their circulation at high frequency.

Differences in epidemiology of enteropathogens in children pre- and post-rotavirus vaccine introduction in Kilifi, coastal Kenya

Background

Kenya introduced Rotarix^® (GlaxoSmithKline Biologicals, Rixensart, Belgium) vaccination into its national immunization programme beginning July 2014. The impact of this vaccination program on the local epidemiology of various known enteropathogens is not fully understood.

Methods

We used a custom TaqMan Array Card (TAC) to screen for 28 different enteropathogens in 718 stools from children aged less than 13 years admitted to Kilifi County Hospital, coastal Kenya, following presentation with diarrhea in 2013 (before vaccine introduction) and in 2016–2018 (after vaccine introduction). Pathogen positivity rate differences between pre- and post-Rotarix^® vaccination introduction were examined using both univariate and multivariable logistic regression models.

Results

In 665 specimens (92.6%), one or more enteropathogen was detected, while in 323 specimens (48.6%) three or more enteropathogens were detected. The top six detected enteropathogens were: enteroaggregative Escherichia coli (EAggEC; 42.1%), enteropathogenic Escherichia coli (EPEC; 30.2%), enterovirus (26.9%), rotavirus group A (RVA; 24.8%), parechovirus (16.6%) and norovirus GI/GII (14.4%). Post-rotavirus vaccine introduction, there was a significant increase in the proportion of samples testing positive for EAggEC (35.7% vs. 45.3%, p = 0.014), cytomegalovirus (4.2% vs. 9.9%, p = 0.008), Vibrio cholerae (0.0% vs. 2.3%, p = 0.019), Strongyloides species (0.8% vs. 3.6%, p = 0.048) and Dientamoeba fragilis (2.1% vs. 7.8%, p = 0.004). Although not reaching statistical significance, the positivity rate of adenovirus 40/41 (5.8% vs. 7.3%, p = 0.444), norovirus GI/GII (11.2% vs. 15.9%, p = 0.089), Shigella species (8.7% vs. 13.0%, p = 0.092) and Cryptosporidium spp. (11.6% vs. 14.7%, p = 0.261) appeared to increase post-vaccine introduction. Conversely, the positivity rate of sapovirus decreased significantly post-vaccine introduction (7.8% vs. 4.0%, p = 0.030) while that of RVA appeared not to change (27.4% vs. 23.5%, p = 0.253). More enteropathogen coinfections were detected per child post-vaccine introduction compared to before (mean: 2.7 vs. 2.3; p = 0.0025).

Conclusions

In this rural Coastal Kenya setting, childhood enteropathogen infection burden was high both pre- and post-rotavirus vaccination introduction. Children who had diarrheal admissions post-vaccination showed an increase in coinfections and changes in specific enteropathogen positivity rates. This study highlights the utility of multipathogen detection platforms such as TAC in understanding etiology of childhood acute gastroenteritis in resource-limited regions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13099-022-00506-z.

Molecular epidemiology, genetic diversity, and vaccine availability of viral acute gastroenteritis in the middle East and North Africa (MENA) region

Acute gastroenteritis is the cause of considerable mortality and morbidity worldwide, particularly among children under five years in underdeveloped countries. Most acute gastroenteritis (AGE) cases are attributed to viral etiologies, including rotavirus, norovirus, adenovirus, astrovirus, and sapovirus. This paper aimed to determine the prevalence rate of different viral etiologies of AGE in the Middle East and North Africa (MENA) region. Moreover, this paper explored rotavirus phylogenetic relatedness, compared VP7 and VP4 antigenic regions of rotavirus with vaccine strains, and explored the availability of vaccines in the MENA region. The literature search identified 160 studies from 18 countries from 1980 to 2019. The overall prevalence of rotavirus, norovirus, adenovirus, astrovirus, and sapovirus were 29.8 %, 13.9 %, 6.3 %, 3.5 %, and 3.2 % of tested samples, respectively. The most common rotavirus genotype combinations in the MENA region were G1P[8], G9P[9], and G2P[4], whereas GII.4 was the predominant norovirus genotype all of which were reported in almost all the studies with genotyping data. The comparison of VP7 and VP4 between circulating rotavirus in the MENA region and vaccine strains has revealed discrete divergent regions, including the neutralizing epitopes. Rotavirus vaccine was introduced to most of the countries of the MENA region; however, only a few studies have assessed the effectiveness of vaccine introduction. This paper provides a comprehensive update on the prevalence of the different viral agents of AGE in the MENA region.

Genetic characterization of P[8] rotavirus strains circulated in Iran between 2009 and 2017

Group A rotavirus (RVA) is the most common cause of acute gastroenteritis (AGE) worldwide, which is responsible for causing an estimated 120,000 deaths in children under 5 years of age, which mostly occur in the lower income countries of Asia and Africa. The G1P[8] is a common genotype of RVA that has spread throughout the world, including Iran and this genotype is present in two commonly used RVA vaccines, Rotarix^TM and RotaTeq^TM . In this study, we investigated the genetic diversity, viral evolutionary, and differences between antigenic epitopes of Iran's P[8] strains and two licensed vaccines. The phylogenetic and evolutionary analysis was carried out, using MEGA vs 6.0 and BEAST respectively. Antigenic epitopes of VP8* were compared to determine the differences between strains from Iran and Rotarix^TM and RotaTeq^TM . The P[8]-lineages III and IV was found as the predominant P genotype that circulating in Iran. The TMRCA of P[8]-lineages III and IV was estimated at 1987 and 2009 respectively. The P[8]-lineage III strains showed 12 amino acid changes compared to Rotarix^TM and 10 amino acid changes compared to RotaTeq^TM . The P[8]-lineage IV strains showed 10 amino acid variations for both Rotarix^TM and RotaTeq^TM strains. The results revealed that the P[8] strains circulating in Iran differs from Rotarix^TM and RotaTeq^TM strains. To monitor the long-term effects of vaccines on the emergence of P[8] strains with different lineages, routine and successful monitoring of these strains will be crucial. This article is protected by copyright. All rights reserved.

Effect of rotavirus genetic diversity on vaccine impact

Group A rotaviruses (RVAs) are the leading cause of gastroenteritis, causing 0.2 million deaths and several million hospitalisations globally each year. Four rotavirus vaccines (Rotarix^TM , RotaTeq^TM , Rotavac^® and ROTASIIL^® ) have been pre-qualified by the World Health Organization (WHO), but the two newly pre-qualified vaccines (Rotavac^® and ROTASIIL^® ) are currently only in use in Palestine and India, respectively. In 2009, WHO strongly proposed that rotavirus vaccines be included in the routine vaccination schedule of all countries around the world. By the end of 2019, a total of 108 countries had administered rotavirus vaccines, and 10 countries have currently been approved by Gavi for the introduction of rotavirus vaccine in the near future. With 39% of global coverage, rotavirus vaccines have had a substantial effect on diarrhoeal morbidity and mortality in different geographical areas, although efficacy appears to be higher in high income settings. Due to the segmented RNA genome, the pattern of RVA genotypes in the human population is evolving through interspecies transmission and/or reassortment events for which the vaccine might be less effective in the future. However, despite the relative increase in some particular genotypes after rotavirus vaccine use, the overall efficacy of rotavirus mass vaccination worldwide has not been affected. Some of the challenges to improve the effect of current rotavirus vaccines can be solved in the future by new rotavirus vaccines and by vaccines currently in progress.

Tate J, D. Parashar U, Kang G. Rotavirus Strain Distribution before and after Introducing Rotavirus Vaccine in India

In April 2016, an indigenous monovalent rotavirus vaccine (Rotavac) was introduced to the National Immunization Program in India. Hospital-based surveillance for acute gastroenteritis was conducted in five sentinel sites from 2012 to 2020 to monitor the vaccine impact on various genotypes and the reduction in rotavirus positivity at each site. Stool samples collected from children under 5 years of age hospitalized with diarrhea were tested for group A rotavirus using a commercial enzyme immunoassay, and rotavirus strains were characterized by RT-PCR. The proportion of diarrhea hospitalizations attributable to rotavirus at the five sites declined from a range of 56–29.4% in pre-vaccine years to 34–12% in post-vaccine years. G1P[8] was the predominant strain in the pre-vaccination period, and G3P[8] was the most common in the post-vaccination period. Circulating patterns varied throughout the study period, and increased proportions of mixed genotypes were detected in the post-vaccination phase. Continuous long-term surveillance is essential to understand the diversity and immuno-epidemiological effects of rotavirus vaccination.

Multiple Introductions and Predominance of Rotavirus Group A Genotype G3P[8] in Kilifi, Coastal Kenya, 4 Years after Nationwide Vaccine Introduction

Globally, rotavirus group A (RVA) remains a major cause of severe childhood diarrhea, despite the use of vaccines in more than 100 countries. RVA sequencing for local outbreaks facilitates investigation into strain composition, origins, spread, and vaccine failure. In 2018, we collected 248 stool samples from children aged less than 13 years admitted with diarrheal illness to Kilifi County Hospital, coastal Kenya. Antigen screening detected RVA in 55 samples (22.2%). Of these, VP7 (G) and VP4 (P) segments were successfully sequenced in 48 (87.3%) and phylogenetic analysis based on the VP7 sequences identified seven genetic clusters with six different GP combinations: G3P[8], G1P[8], G2P[4], G2P[8], G9P[8] and G12P[8]. The G3P[8] strains predominated the season (n = 37, 67.2%) and comprised three distinct G3 genetic clusters that fell within Lineage I and IX (the latter also known as equine-like G3 Lineage). Both the two G3 lineages have been recently detected in several countries. Our study is the first to document African children infected with G3 Lineage IX. These data highlight the global nature of RVA transmission and the importance of increasing global rotavirus vaccine coverage.

Epidemiological Trends of Five Common Diarrhea-Associated Enteric Viruses Pre- and Post-Rotavirus Vaccine Introduction in Coastal Kenya

Using real-time RT-PCR, we screened stool samples from children aged <5 years presenting with diarrhea and admitted to Kilifi County Hospital, coastal Kenya, pre- (2003 and 2013) and post-rotavirus vaccine introduction (2016 and 2019) for five viruses, namely rotavirus group A (RVA), norovirus GII, adenovirus, astrovirus and sapovirus. Of the 984 samples analyzed, at least one virus was detected in 401 (40.8%) patients. Post rotavirus vaccine introduction, the prevalence of RVA decreased (23.3% vs. 13.8%, p < 0.001) while that of norovirus GII increased (6.6% vs. 10.9%, p = 0.023). The prevalence of adenovirus, astrovirus and sapovirus remained statistically unchanged between the two periods: 9.9% vs. 14.2%, 2.4% vs. 3.2 %, 4.6% vs. 2.6%, (p = 0.053, 0.585 and 0.133), respectively. The median age of diarrhea cases was higher post vaccine introduction (12.5 months, interquartile range (IQR): 7.9–21 vs. 11.2 months pre-introduction, IQR: 6.8–16.5, p < 0.001). In this setting, RVA and adenovirus cases peaked in the dry months while norovirus GII and sapovirus peaked in the rainy season. Astrovirus did not display clear seasonality. In conclusion, following rotavirus vaccine introduction, we found a significant reduction in the prevalence of RVA in coastal Kenya but an increase in norovirus GII prevalence in hospitalized children.