Whole-Genome Characterization of Rotavirus G9P[6] and G9P[4] Strains That Emerged After Rotavirus Vaccine Introduction in Mozambique

Mozambique introduced the Rotarix vaccine into the National Immunization Program in September 2015. Following vaccine introduction, rotavirus A (RVA) genotypes, G9P[4] and G9P[6], were detected for the first time since rotavirus surveillance programs were implemented in the country. To understand the emergence of these strains, the whole genomes of 47 ELISA RVA positive strains detected between 2015 and 2018 were characterized using an Illumina MiSeq-based sequencing pipeline. Of the 29 G9 strains characterized, 14 exhibited a typical Wa-like genome constellation and 15 a DS-1-like genome constellation. Mostly, the G9P[4] and G9P[6] strains clustered consistently for most of the genome segments, except the G- and P-genotypes. For the G9 genotype, the strains formed three different conserved clades, separated by the P type (P[4], P[6] and P[8]), suggesting different origins for this genotype. Analysis of the VP6-encoding gene revealed that seven G9P[6] strains clustered close to antelope and bovine strains. A rare E6 NSP4 genotype was detected for strain RVA/Human-wt/MOZ/HCN1595/2017/G9P[4] and a genetically distinct lineage IV or OP354-like P[8] was identified for RVA/Human-wt/MOZ/HGJM0644/2015/G9P[8] strain. These results highlight the need for genomic surveillance of RVA strains detected in Mozambique and the importance of following a One Health approach to identify and characterize potential zoonotic strains causing acute gastroenteritis in Mozambican children.