Genomic constellation of human Rotavirus A strains identified in Northern Brazil: a 6-year follow-up (2010-2016)

Surveillance of Rotavirus A (RVA) throughout the national territory is important to establish a more complete epidemiological-molecular scenario of this virus circulation in Brazil. The aim of the present study was to investigate the genetic diversity of RVA strains circulating in Tocantins State (Northern Brazil) during six years of post-vaccination follow-up (2010-2016). A total of 248 stool samples were screened by next generation sequencing and 107 (43.1%) nearly full length RVA genome sequences were obtained; one sample was co-infected with two RVA strains (G2/G8P[4]). Six G and P genotypes combinations were detected: G12P[8] strains (78.6%), as well as the G3P[8] (9.3%) and G1P[8] (0.9%) were associated with a Wa-like genogroup backbone. All G2P[4] (5.6%) and G8P[4] (2.8%) strains, including the mixed G2/G8P[4] infection (0.9%) showed the DS-1-like genetic background. The two G12P[4] strains (1.9%) were associated with distinct genetic backbones: Wa-like and DS-1-like. The phylogenetic analysis revealed the circulation of lineages G1-I, G2-IV, G3-III, G8-I and G12-III, and P[4]-V and P[8]-III of the VP7 and VP4 genes, respectively. Conserved clustering pattern and low genetic diversity were observed regarding VP1-VP3 and VP6, as well as NSP1-5 segments. We identified the same RVA circulation pattern reported in other Brazilian regions in the period of 2010-2016, suggesting that rural and low-income areas may not have a different RVA genotypic distribution compared to other parts of the country. The unique presentation of whole-genome data of RVA strains detected in the Tocantins State provides a baseline for monitoring variations in the genetic composition of RVA in this area.