RVA detection. From July 2020 to June 2022, in total 4,125 samples of faeces of children hospitalised in an infectious hospital in Nizhny Novgorod with acute intestinal infection were examined.
Rotavirus RNA was detected in 30.4% of cases (n=1253), of which the G[P] type was determined in 81.6% (n=1022). RT-PCR typing revealed the presence of the following combinations of G/P types in the studied period: G1P[8], G2P[4], G4P[8], G3P[8], G9P[8], G3P[6], G3P[9], G2P[8], G8P[8], G12P[6], G6P[9], G9P[4], G3P[4]; as well as partially typed strains with genotypes: G2P[x], G3P[x], G8P[x], G9P[x], GxP[4], GxP[6], GxP[8], GxP[9]. Figure 1 shows that the proportion of G3P[8] RVA increased from 15% in the period 2020-2021 to 53% in the season 2021-2022.
The G3P[8] rotaviruses identified in this work, in addition to typing by RT-PCR, were analysed by RNA-PAGE. Combination of RT-PCR and RNA-PAGE made it possible to establish the presence of G3P[8] RVA with short, DS-1-like profile of segment migration in polyacrylamide gel (PAGE), which clearly indicated intergroup reassortment (Fig. 2).
Phylogenetic analysis of the VP7 gene. Phylogenetic analysis of the VP7 gene of G3P[8] RVAs was performed based on 42 sequences obtained in this study and sequences available in the GenBank database (n=544). On the formed MCC (Maximum Clade Credibility) tree, two large clusters I and II are identified (Fig. 3). A complete phylogenetic tree with the names of sequences and their numbers in the GenBank database is provided in additional materials (Suppl Fig.1).
In cluster I, 3 lines can be reliably distinguished. The I-1 line includes G3P[8] RVA isolated in different regions of the world in 1998-2017. The I-2 and I-3 lines include RVAs with genotypes G3P[6] and G3P[9] isolated in 1997-2016 and 1996-2014, respectively. In addition, cluster I includes a number of G3 rotaviruses isolated from pigs in Europe, Africa and Asia in 2009-2017 in combination with genotypes P[6], P[7], P[13], P[22], P[23], P[32] of the VP4 gene (in Figure 2 marked as Porcine G3).
Four lines are allocated in cluster II. The II-1 line is mainly represented by G3P[8] DS-1-like rotaviruses which have spread widely around the world since 2013. It should be noted that the same line includes RVA with the G3 genotype, which was isolated in India from horses in 2004, as well as G3P[45] rotaviruses isolated in China from rats in 2013. Apparently, the formation of DS-1-like G3P[8] RVA occurred in Asia with the participation of local animal populations. Lines II-2 to II-4 are represented by animal rotaviruses (bat, bovine, rodents (Mouse, Rat), horse, rabbit, simian etc.), which were isolated in different years (Fig. 3).
Philodynamic analysis of the VP7 gene. In this paper, a discrete phylogeographic analysis of the RVA strains belonging to line II-1, shown in Figure 3, was carried out. The estimated mutation rate was 2.14*10-3 substitutions per site per year [95% HPD: 1.38x10-4; 3.18x10-3]. Data of spatio-temporal reconstruction suggest that the spread of the G3 line originated around 2002 [95% HPD: 2000-2007]. Based on the analysis of all relevant sequences available in GenBank, the closest ancestor of the G3 lineage considered and the RVA strain isolated from a horse in India in 2004 was circulating in 1980 [95% HPD: 1974-1989]. The obtained phylogenetic tree with sufficient probability suggests that the origin of DS-1-like G3P[8] rotavirus originated from animals (including horses, rats, bats, cattle) and the source of its distribution is East Asia (India, China) (Fig. 4). The formed variant with the genotype G3P[8] quickly spread around the world. Table 1 shows the values of BF>3 obtained using the SPREAD 3D software. If the value of this factor exceeds 3, this confirms a strong epidemiological connection between the regions being compared. However, since the spread of the DS-1-like variant of G3P[8] occurred in a short time and widely, it is not possible to establish a reliable source of drift for each region (Fig. 3 Plot map).
Phylogenetic analysis of the VP4 gene. According to the classification of da Silva et al., DS-1-like G3P[8] rotaviruses isolated in Nizhny Novgorod by the VP4 gene belong to the line P[8]-3.1 [22]. The nucleotide sequences obtained in this study were clustered together with strains isolated in other countries of Europe, East and Southeast Asia, and Russia in combination with the G3 genotype. Some isolates (a smaller group) were clustered with rotaviruses of G8, G9 and G3 genotypes isolated in Russia, East and Southeast Asia and Europe (Fig. 5).
Phylogenetic analysis of the VP6 gene. During the study, we obtained 25 nucleotide sequences of the VP6 gene of G3 DS-1-like rotaviruses. All of them belonged to genotype I2. In the resulting tree, we identified three lines: I2-I to I2-III (Fig. 6). Interestingly, 23 of the 25 sequenced sequences formed an I2-1 cluster together with G2P[4] RVA strains, which were isolated in Russia and other countries (Fig. 4). G3P[8] RVA isolated in Novosibirsk in 2018 belongs to the same cluster.
The I2-II cluster was formed by RVAs with genotypes G1P[8], G2P[4], G3P[8], G8P[8] and G9P[8], it also included one sequence of the VP6 rotavirus gene with genotype G3P[8] isolated in Nizhny Novgorod in 2021. Cluster I2-III included a single sequence of RVA with the genotype G3P[4], which was isolated by us in 2022. It is noteworthy that the third cluster includes the majority of DS-1-like G3P[8] RVAs described in the literature, including strains isolated in Moscow in 2020. In addition, the third group also included intergroup reassortants with genotype G1P[8] isolated in Nizhny Novgorod in 2014 [10].