Genetic diversity of G9, G3, G8 and G1 rotavirus group A strains circulating among children with acute gastroenteritis in Vietnam from 2016 to 2021

Abstract

Rotavirus group A (RVA) is the most common cause of severe childhood diarrhea worldwide. The introduction of rotavirus vaccination programs has contributed to a reduction in hospitalizations and mortality caused by RVA. From 2016 to 2021, we conducted surveillance to monitor RVA prevalence and genotype distribution in Nam Dinh and Thua Thien Hue (TT Hue) provinces where a pilot Rotavin-M1 vaccine (Vietnam) implementation took place from 2017 to 2020. Out of 6626 stool samples, RVA was detected in 2164 (32.6%) by ELISA. RT-PCR using type-specific primers were used to determine the G and P genotypes of RVA-positive specimens. Whole genome sequences of a subset of 52 specimens randomly selected from 2016 to 2021 were mapped using next-generation sequencing. From 2016 to 2021, the G9, G3 and G8 strains dominated, with detected frequencies of 39%, 23%, and 19%, respectively; of which, the most common genotypes identified were G9P[8], G3P[8] and G8P[8]. G1 strains re-emerged in Nam Dinh and TT Hue (29.5% and 11.9%, respectively) from 2020 to 2021. G3 prevalence decreased from 74% to 20% in TT Hue and from 21% to 13% in Nam Dinh province between 2017 and 2021. The G3 strains consisted of 52% human typical G3 (hG3) and 47% equine-like G3 (eG3). Full genome analysis showed substantial diversity among the circulating G3 strains with different backgrounds relating to equine and feline viruses. G9 prevalence decreased sharply from 2016 to 2021 in both provinces. G8 strains peaked during 2019-2020 in Nam Dinh and TT Hue provinces (68% and 46%, respectively). Most G8 and G9 strains had no genetic differences over the surveillance period with very high nucleotide similarities of 99.2-99.9% and 99.1-99.7%, respectively. The G1 strains were not derived from the RVA vaccine. Changes in the genotype distribution and substantial diversity among circulating strains were detected throughout the surveillance period and differed between the two provinces. Determining vaccine effectiveness against circulating strains over time will be important to ensure that observed changes are due to natural secular variation and not from vaccine pressure.

Keywords: Acute gastroenteritis (AGE); Genotype distribution; Rotavin-M1; Rotavirus group A (RVA); Vietnam; Whole genome sequence.

Fig. 1.

Prevalence of RVA G genotypes in Nam Dinh (a) and TT Hue (b) provinces during 2016–2021.

Fig. 2.

Phylogenetic tree constructed with the nucleotide sequences of the VP7 gene of G1 (A), G3 (B), G8 (C) and G9 (D) strains detection in Vietnam (•) from 2016 to 2021. The number adjacent to the node represents the bootstrap value and values lower than 70% have not been indicated. The scale bar at the bottom of the figure represents genetic distance expressed as the number of substitutions per site. Only some special Vietnam RV strains and the VP7 genes of G3 strains were deposited in GenBank, with accession numbers from OR756342 to OR756448.