Molecular and genetic characteristics of group A rotaviruses detected in Moscow in 2015–2020

Cover Page


Cite item

Full Text

Abstract

The aim of the study was to analyze genetic characteristics of strains belonging to group A rotaviruses (RVA) circulating in Moscow in 2015–2020, including rare strains non-typeable by polymerase chain reaction (PCR).

Materials and methods. A total of 289 stool samples were tested; the samples were collected from children aged 1 month to 17 years, hospitalized with acute gastroenteritis. Immunochromatography and real-time reverse transcription-polymerase chain reaction (real-time RT-PCR) assays were used for detection of rotaviruses in the samples. The rotavirus genome sequencing was performed using the Sanger technique and nanopore sequencing.

Results and discussion. RVA RNA was detected in 131 clinical samples, and the G/[P] genotype was identified in 125 samples. The general profile showed prevalence of RVA strains with the G9P[8]I1 genotype (37%) followed by G3P[8]I2, G4P[8]I1, G2P[4]I2, G1P[8]I1, and G3P[8]I1 variants (18, 15, 11, 5, and 2%, respectively). Seven (5%) isolates were identified as GxP[8]. In 2015–2020, the region reported a decline in G4P[8]I1 genotype prevalence (from 39% to 9%) and an increase in the proportion of the G9P[8]I1 genotype (from 6% to 37%) as compared to 2009–2014. In 2018–2020, a large number of cases with the previously unknown DS-1-like reassortant strain with the G3P[8]I2 genotype were reported; the above strain has become widely common worldwide in the recent years. Nanopore sequencing was performed to analyze the genome of the G3P[8]I2 strain and the rare G4P[6]I1 strain. It was found that the G4P[6]I1 strain was phylogenetically related to porcine rotaviruses.

Conclusion. In the recent years, the genetic diversity of RVA circulating in the Moscow Region has changed significantly. The obtained results prove the importance of continuous monitoring of rotavirus infection and selective sequencing of RVA genes to fine-tune data of the type-specific real-time RT-PCR. The ever-changing genetic composition of the circulating RVA strains calls for regular optimization of RVA genotyping systems based on real-time RT-PCR. 

Full Text

Introduction

Group A rotaviruses (RVA) are a major cause of acute gastroenteritis hospitalizations of children aged under 5 years in countries characterized by a low rotavirus vaccination coverage. In 2016, rotavirus infection (RVI) was responsible for 258 million episodes of diarrhea and 128,500 deaths among children younger than 5 years worldwide [1].

In 2020, vaccination against RVA was included in the National Immunization Programs in 107 countries1. In Russia, immunization against RVI is included in the National Vaccination Schedule as required for epidemic reasons. Currently, four rotavirus vaccines have been recommended by WHO and approved for application: pentavalent vaccines — RotaTeq (Merck & Co., Inc.), ROTASIIL (Serum institute of India PVT. Ltd.); and monovalent vaccines — Rotarix (GlaxoSmithKline), ROTAVAC (Bharat Biotech). Mass vaccination against RVA decreases incidence and acute gastroenteritis hospitalization rates in all age groups, especially among infants and people over 65 years [2].

RVAs belong to species Rotavirus A, genus Rotavirus, family Reoviridae. The rotavirus genome consists of 11 segments of double-stranded RNA, which encode 12 proteins [3]. The present-day system of rotavirus classification offers genotyping for all 11 genes (Gх–P[ х]–Iх–Rх–Cх–Mх–Aх–Nх–Tх–Eх–Hх). Most of the circulating human RVAs belong to 3 evolutionary lines different in genome constellations: Wa-like strains (G1–P[ 8]–I1–R1–C1–M1–A1–N1–T1–E1–H1), DS1-like (G2–P[ 4]–I2–R2–C2–M2–A2–N2–T2–E2–H2), and AU-1-like strains (G3–P[ 3]–I3–R3–C3–M3–A3– N3–T3–E3–H3) [3]. The previously used binary system of typing addressed only VP7 (G) and VP4 (P) proteins. G/[P]-genotyping of RVA strains is performed by an reverse transcription followed by multiplex polymerase chain reaction (RT-PCR) [4] or sequencing of the respective genes.

At present, G1P[ 8], G2P[ 4], G3P[ 8], G4P[ 8], G12P[ 8], and G9P[ 8] genotypes, homotypic or partially heterotypic towards the licensed vaccines (RotaTeq and Rotarix), prevail among circulating RVA strains worldwide [5]. There have been described cases of interspecies transmission of RVA from animals to humans as well as cases when humans were infected with strains resulting from reassortment of animal and human rotavirus strains [6–8].

The evolutionary changes of RVA strains over time [9], territorial differences in distribution of circulating strains, their increased diversity after implementation of immunization in some regions, changes in the genotypic structure both in countries with and without scheduled vaccination require continuous epidemiological monitoring of RVI [10].

The aim of the study is the analysis of genetic characteristics of the RVA strains circulating in Moscow in 2015–2020, including rare strains, which are non-typeable by PCR tests.

Materials and methods

Clinical samples

During 5 years (2015–2020), a total of 289 stool samples were collected from children aged 1 month to 17 years, displaying symptoms of acute gastroenteritis and hospitalized to the Vladimirsky Moscow Regional Research Clinical Institute. Children vaccinated against RVI were not included in the study. The informed consent was received from the parents or legal guardians of all tested children.

A total of 45 samples were identified as positive for the RVA antigen using immunochromatography (the RIDA Quick Rotavirus reagent kit, R-Biopharm AG). Stool samples were collected from patients with diarrhea (not later than 72 hours after the onset of symptoms) into sterile containers; before they were shipped to the laboratory, samples were stored at –20ºC, then were shipped to the Mechnikov Research Institute for Vaccines and Sera for genetic analysis of rotaviruses.

RNA extraction

To isolate RNA, we used 10% fecal extracts diluted in sterile saline and cleared through centrifugation (5,000 rpm, 5 min). To isolate RNA from the extracts, we used an AmpliSens® MAGNO-sorb kit (InterLabService) in accordance with the manual. The RNA samples were stored at 80ºС.

Detection of rotavirus RNA

To detect RVA in clinical samples by real-time RT-PCR, we used single-tube reaction with primers and probes described earlier [11]. Genotyping RVA using two-stage multiplex real-time PCR for VP7 proteins (G1, G2, G3, G4, G9), VP4 (P[ 4], P[ 6], P[ 8]), and VP6 (I1, I2) was performed in accordance with the procedure requirements [12]. The analysis of RVA genotyping based on multiplex RT-PCR as well as the subsequent analysis of PCR products using agarose gel electrophoresis were used as reference methods and performed in accordance with the WHO recommendations [4].

Amplification and sequencing of VP7 and VP4 genes

For sequencing of VP7 and VP4 genes of G3P[ 8] I2 strains, we performed amplification using the previously described VP7F/VP7R [13] and VP4F/VP4R [14] primers generating amplicons 881 and 663 base pairs long, respectively. Aliquots of the extracted RNA (10 μl) were mixed with 3 pmol VP7F or VP4F primers, incubated at 95ºC for 1 min and cooled down for 2–3 min to room temperature. The RT test was performed in the 25 μl reaction mixture containing an RT primer, 25 units of MMLV reverse transcriptase (Syntol, Russia). The RT stage included incubation for cDNA synthesis at 45ºC for 30 min and inactivation of the MMLV reverse transcriptase at 95ºC for 5 min. The temperature and time parameters for the real-time PCR were as follows: 95ºC — 120 sec; 95ºC — 60 sec, 52ºC — 40 sec, 72ºC — 40 sec (45 cycles); 72ºC — 40 sec. PCR-amplicons of each gene were purified with a Сleanup Standard kit (Eurogen); the sequencing of both DNA strands was performed using a NANOPHORE ®05 genetic analyzer (Syntol) and a reagent kit from the Syntol Company.

Nanopore sequencing of the rotavirus genome and bioinformatic analysis

The full-genome cDNA of rotavirus gene segments was produced using a mix of unRAf1, unRAf2, unRAf3, unRAr1, unRAr2, and unRAr3 primers (Table 1) for the RT test. For concurrent amplification of all segments of the RVA gene, we used universal primer Up [15] (Table 1). The amplification was performed in the following reaction mixture: TaKaRa LA Taq polymerase (2.5 units), LA PCR Buffer ll (Mg2+plus), 25 mM MgCl2 (final concentration Mg2+ 5 mM) (TaKaRa), dNTP 25 mM (Syntol), 40 pmol Up primer. The amplification parameters: 95ºC — 2 min, then 40 cycles at 95ºC, each for 30 sec, 65ºC — 30 sec, 68ºC — 3.5 min. PCR amplicons were purified using phenol/chloroform extraction. The concentration and purity of amplicons were measured by spectrophotometry (A260/230, 260/280).

 

Table 1. Primers for amplification of full length rotavirus gene segments

Primer

Primer sequence 5’-3’

unRAf1

GCCGGAGCTCTGCAGAATTCGGCTWTWAAA

unRAf2

GCCGGAGCTCTGCAGAATTCGGCTTTTTTT

unRAf3

GCCGGAGCTCTGCAGAATTCGGCTTTTAAT

unRAr1

GCCGGAGCTCTGCAGAATTCGGTCAYATC

unRAr2

GCCGGAGCTCTGCAGAATTCGGTCACAWA

unRAr3

GCCGGAGCTCTGCAGAATTCAGCCACATG

Up

GCCGGAGCTCTGCAGAATTC

The DNA library for nanopore sequencing (NPS) was constructed using a Rapid Sequencing reagent kit with a portable MinION sequencer and standard flowcell R9 (Oxford Nanopore Technologies). We developed a pipeline for accurate classification of NPS data. The Python programming language-based pipeline identified the received reads, actuating the BLAST tool-based analysis for the database of reference sequences of rotavirus genome segments2.

Then, the reads were mapped to the reference sequence using the Minimap2 program [16]; the consensus sequence was created using the scriptsup>3 written in the Python language. In the consensus sequence, we selected a nucleotide with the highest frequency at the alignment position. In rare cases, when two nucleotides were found at the same position in the equal quantity, the second nucleotide was ignored. The RVA genotype was identified by the nucleotide sequence of gene segments using the Rotavirus A Genotype Determination online service4 based on the software from Dan Katzel [17].

The method of nanopore sequencing of the rotavirus genome is described in detail in the article by Faizuloev et al. [18].

GenBank accession numbers (NCBI)

The nucleotide sequences corresponding to 10 segments of the genome of Moscow-40/2020 (G3P[ 8]I2) and Moscow-1P/2015 (G4P[ 6]I1) isolates have been deposited to GenBank under numbers MW558493– MW558502 and MT876633–MT876642, respectively. The partial sequences of three other strains with the G3P[ 8]I2 genotype have been deposited under numbers MT648671, MT648671, MT648671 for VP7 and MT814324, MT814326 for VP4.

Phylogenetic analysis

Phylogenetic trees were built using the MEGA-X program [19] and the maximum likelihood method (the Kimura two-parameter model [20]). The bootstrap test included 1,000 replications. The RVA strains recommended by the Rotavirus Classification Working Group [3] and the strains having, based on the BLAST data, at least 99.5% similarity to the studied strains were used as reference strains.

Results

A total of 289 stool samples were tested for presence of the RNA or RVA antigen; all of them were collected from children hospitalized with symptoms of acute gastroenteritis. RNA was detected in 131 (45%) samples and was further used for G/[P]-genotyping by real-time RT-PCR and/or sequencing. The RVA G/ [P]-genotype was detected in 125 samples (95.4%), and in 7 (5.3%) samples. RVA was defined only by the P gene, while no genotype was identified for 6 (4.6%) samples.

The genetic composition of the studied RVA strains is presented in Fig. 1. In 2015–2020, RVA strains with the G9P[ 8]I1 genotype prevailed in the overall profile (37%), the second place was taken by G3P[ 8]I2 (a new DS-1-like strain) 18%, which was followed by G4P[ 8]I1 (15%), G2P[ 4]I2 (11%), G1P[ 8]I1 (5%), and G3P[ 8]I1 (2%). Seven (5%) strains defined only by the P gene belonged to the P gene variant [8]. We have also found a single case of co-infection with two genotypes (G9P[ 8]I1 and G2P[ 4]I2) and a rare strain with the G4P[ 6]I1 genotype.

 

Fig. 1. Distribution of RVA G/[P]-genotypes isolated from 131 clinical samples in 2015–2020. GxP[ 8]Ix — partially typed PVA; Mix — mixed genotypes; Nt — non-typed samples.
 

Fig. 2 shows the year-to-year distribution of the RVA genotypes identified in the Moscow Region. In the last years, the proportion of the G9P[ 8]I1 genotype has increased significantly, ranging from 36% to 41% during 2015-2020. At the same time, the prevalence of G4P[ 8] I1 genotypes decreased from 38% to 9%, while the prevalence of the G2P[ 4]I2 genotype increased to 14%.

 


Fig. 2. Prevalence of G/[P]-genotypes detected in Moscow during the 10-year monitoring period (2009–2020) [12]. Nt — non-typed samples.

 

Note that the partial sequencing of VP7 (G) and VP4 (P) genes was required for G/[P]-genotyping of some strains with the GxP[ 8]I2 genotype, as the realtime RT-PCR-based laboratory assay was not able to identify their G/[P]-genotype [12]. The sequencing and phylogenetic analysis of RVA G and P genes (the GenBank numbers: MT648671, MT648671, MT648671 for the G gene and MT814324, MT814326 for the P gene) in the sample with detected GxP[ 8]I2 showed that they belong to the G3P[ 8]I2 genotype (Fig. 3). It was found that the absence of a signal during real-time RT-PCR genotyping was caused by the mismatch of nucleotides between the VP7 gene with the G3P[ 8]I2 genotype and the respective probe. The non-typeable samples with the GxP[ 8]I2 genotype were tested with mono-specific PCR with primers to G3 and electrophoretic detection, which identified amplicons with the expected mobility in all samples (the data are not provided), thus confirming that they belonged to the G3P[ 8]I2 genotype.

 

Fig. 3. Phylogenetic trees based on the sequenced genes of the VP4 and VP7 proteins of the new G3P[ 8]I2 strain. ▲ — segments sequenced using the Sanger technique; ♦ — the gene sequenced using the NPS technique; Δ — the equine RVA VP7 gene. The respective GenBank accession number, name of strains and G/[P]-genotype were used for designation of strains.
 

The BLAST search and the phylogenetic analysis demonstrated the similarity of the sequenced VP7 and VP4 genes to the same RVA genes of the new G3DS-1-like constellation detected worldwide. The similar G3-DS-1-like strains were detected in Australia (2013) [21], Spain (2014–2015) [22], Hungary (2016) [23], Brazil (2016) [24], Indonesia (2016) [25], Russia (2019) [28], and other countries (Fig. 3).

In addition, the nanopore sequencing of the Moscow40/2020 isolate with the G3P[ 8]I2 genotype was performed (Table 2). It helped identify gene variants of 10 genome segments (the GenBank numbers MW558493–MW558502) and demonstrated a high degree (91.0–99.8%) of similarity between the consensus sequence and the reference strain RVA/Human-wt/ THA/SKT-281/2013/G3P[ 8] (the GenBank numbers LC086714–LC086724). The rare Moscow-1P/2015 strain with the G4P[ 6]I1 genotype was detected in the clinical sample collected from an 8-year-old patient in 2015. Real-time PCR was able to identify only the P[ 6] gene variant. Sanger sequencing and NPS helped identify the genotype of this strain by 10 genes: G4-P6-I1R1-C1-M1-A1-N1-T1-Ex-H1 (the GenBank numbers MT876633–MT876642, MG271938), accounting for 81.7% of the genome (Table 2).

 

Table 2. Results of nanopore sequencing of the genome of RVA strains with G3P[8]I2 and G4P[6]I1 genotypes

Segment

Viral protein

G4P[6]I1

G3P[8]I2

number of reads

segment coverage, %

genotype

number of reads

segment coverage, %

genotype

1

VP1

1,372

100

R1

12,792

100

R2

2

VP2

3,795

100

C1

11,193

100

C2

3

VP3

1,995

100

M1

-

-

Mx

4

VP4

633

20

P[6]

2,785

100

P[8]

5

NSP1

546

100

A1

53,841

100

A2

6

VP6

15,187

100

11

16,479

100

I2

7

NSP3

14,566

100

T1

80,554

100

T2

8

NSP2

175

36

N1

40,022

100

N2

9

VP7

1,465

100

G4

25,059

100

G3

10

NSP4

-

-

Ex

10,222

100

E2

11

NSP5/6

77

100

H1

11,351

100

H2

The phylogenetic analysis of genes VP7, VP6, and VP4 strain Moscow-1P/2015 (Fig. 4) demonstrates a high degree of similarity of the analyzed sample to genes of porcine RVA (the GenBank numbers: VP4 — KX363402, MK227950, KX363435, MK227948; VP6 — MK227391, MK227402, KX363414, MG066585, KJ126830; VP7 — JX498957, JX498956, MK227392, MN133419, MN133444) or RVA strains (strain RVA/Human-wt/CHN/R1954/2013/G4P[ 6] having the GenBank numbers: KF726066–KF726076, KF726056) and isolated from human feces, though having the confirmed origin from porcine RVA [7]. The BLAST analysis of the other 7 genes also demonstrates a high degree of similarity (92–98%) of the nucleotide sequence to porcine RVA strains. Thus, the phylogenetic analysis indicates that the G4P[ 6] strain is of porcine RVA origin.

 


Fig. 4. Phylogenetic trees based on sequenced genes of RVA VP4, VP7, and VP6 proteins, strain G4P[ 6]I1 (marked by ♦), reference strains (Wa, AU-1, DS-1) representatives of three evolutionary lines RVA (marked by ○) and RVA genes of porcine origin, phylogenetically most closely related to genes of the Moscow-1P/2015 strain based on the BLAST analysis. The respective GenBank number, name of the strain and G/[P]-genotype were used for designation of strains.

 

Discussion

The obtained results are indicative of significant changes, which took place in the "genetic landscape" of RVAs circulating in the Moscow Region. Based on the data from the previous studies [12][26], before 2015, most of the genotyped RVA belonged to genotype G4P[ 8]. Our data indicate a gradual decrease in the prevalence of this genotype, from 38% to 9% in 2017–2018. At the same time, the proportion of the G9P[ 8] genotype increased to 36–41% in the Moscow Region. These data are consistent with the data obtained during the independent studies that were performed in Moscow [27][28], Nizhny Novgorod [29][30], and Orenburg [31].

Our study was focused on clinical samples collected from children hospitalized with acute gastroenteritis. The actual "genetic landscape" and the distribution of RVA genotypes circulating in Moscow can differ from the obtained data, since we did not include patients with mild or moderate gastroenteritis, not requiring hospitalization, in our study. Previously, in the study by E.R. Meskina, it was pointed out that severe rotavirus gastroenteritis could be associated with specific RVA genotypes [32].

The strain with genotype G3P[ 8]I2, which was discovered by us in 2018, is of special interest. It accounted for 18% in the RVA group, coming second only to the G9P[ 8]I1 genotype (Fig. 1). Based on the sequencing and phylogenetic analysis (Fig. 3), we can assume that G3P[ 8]I2 is related to the reassortant strain, which was first detected in Thailand in 2013 [25] and has become common in Europe, Asia, and Australia in the last years [7][22][33][34]. This strain has a DS-1-like constellation of genes (G3–P[ 8]–I2–R2–C2–M2–A2–N2–T2–E2–H2), except for the gene of the VP7 protein, which in the cases typical of this constellation is represented by the G2 genovariant, while the Wa-like constellation is more common for the G3 genovariant [35].

Some researchers [21][22][25] assume that the G3P[ 8] I2 strain owes its emergence to reassortment of the human DS-1-like strain and the equine Erv105 strain (the GenBank number: DQ981479.1), as the gene of the VP7 protein of this strain demonstrates the highest similarity. During this study, we, like some other research groups [22], ran into a problem, finding it impossible to genotype a new virus strain using type-specific PCR. Both our own assay and the primers recommended by WHO [4] failed to detect the gene of the VP7 protein of the G3P[ 8]I2 strain. At the same time, the Sanger sequencing and the BLAST analysis were more efficient and made it possible to identify the genovariant as G3. Point mutations in genes VP4 and VP7 in many cases prevent type-specific PCR primers from identifying their variant [17]. Hence, the importance of sequencing methods in RVI epidemiological monitoring has increased.

Another atypical strain — Moscow-1P/2015 with the G4P[ 6]I1 genotype, which was detected only in one case, may have also resulted from the reassortment of the human and animal RVA or may have the animal origin. In this particular case, it is not clear whether the infection with this strain has a human-to-human transmission route, or the patient was infected with the virus from animals. Based on the literature data, the proportion of RVA with the G4P[ 6]I1 genotype is not high in the total diversity of RVA strains detected in humans. It may be caused by some factor limiting the G4P[ 6]I1 spread in the human population, for example, a species barrier, if we assume that this strain is of animal origin. On the other hand, it should be remembered that most of the related studies tend to focus on children hospitalized with rotavirus enteritis, while mild cases are generally left out of studies. Therefore, the available data do not give any reliable information about the actual prevalence of any of the RVA strains.

Porcine RVA strains are phylogenetically related to strains of human rotaviruses [34]; therefore, it is difficult to find out whether the studied strain resulted from a direct transmission or reassortment [37–39]. Cases of reassortment or direct interspecies transmission of RVA with the G4P[ 6]I1 genotype have been reported and described both in Asia [6][25][40] and in Europe [36].

Conclusion

In 2015–2020, the genetic profile of RVA circulating in the Moscow Region changed significantly: The detection frequency of the G4P[ 8]I1 genotype, which was the most prevalent in the previous years, decreased; at the same time, the number of hospitalizations with RVI caused by the G9P[ 8]I1 strain increased. In addition, the above period is characterized by emergence of RVA strains, presumably of animal origin, both in few numbers (G4P[ 6]I1) and in significant numbers (reassortant G3P[ 8]I2), thus suggesting an important role of interspecies transmission in the evolution of RVA pathogenic to humans. Our study demonstrates the importance of continuous monitoring of RVI. Epidemiological monitoring of RVI provides an effective tool for timely detection of new animal and human reassortant RVAs, which may escape post-vaccinal immunity. The RVA stains that we and other researchers have identified and that are not typeable by real-time RT-PCR (genotypes G3P[ 8]I2 and G4P[ 6]I1) [6][22] prove the significance of selective sequencing of RVA genes and the need to optimize sequences of type-specific primers for real-time RT-PCR.

1. URL: https://preventrotavirus.org/vaccine-introduction/global-introduction-status

2. URL: https://github.com/lioj/bioinformatics/blob/master/py/classificationStat.py

3. URL: https://github.com/lioj/bioinformatics/blob/master/py/bam2consensus.py

4. URL: https://www.viprbrc.org/brc/rvaGenotyper.spg?method=ShowCleanInputPage&decorator=reo

×

About the authors

O. A. Petrusha

I. Mechnikov Research Institute of Vaccines and Sera

Author for correspondence.
Email: petrusha.olga@gmail.com
ORCID iD: 0000-0002-5022-7962

Olga A. Petrusha — junior researcher, Laboratory of molecular virology, 

Moscow

Russian Federation

E. R. Korchevaya

I. Mechnikov Research Institute of Vaccines and Sera

Email: fake@neicon.ru
ORCID iD: 0000-0002-6417-3301

Ekaterina R. Korchevaya — junior researcher, Laboratory of molecular virology,

Moscow

Russian Federation

R. R. Mintaev

I. Mechnikov Research Institute of Vaccines and Sera

Email: fake@neicon.ru
ORCID iD: 0000-0001-5398-3627

Ramil R. Mintaev — junior researcher, Laboratory of genetics of RNA viruses,

Moscow

Russian Federation

A. A. Nikonova

I. Mechnikov Research Institute of Vaccines and Sera

Email: fake@neicon.ru
ORCID iD: 0000-0002-5742-6550

Alexandra A. Nikonova — Cand. Sci. (Biol.), Head, Laboratory of molecular biotechnology, 

Moscow

Russian Federation

I. Yu. Isakov

I. Mechnikov Research Institute of Vaccines and Sera

Email: fake@neicon.ru
ORCID iD: 0000-0001-9610-0935

Igor Yu. Isakov — junior researcher, Laboratory of molecular biotechnology, 

Moscow

Russian Federation

E. R. Meskina

I. Mechnikov Research Institute of Vaccines and Sera;
M. Vladimirsky Moscow Regional Research Clinical Institute

Email: fake@neicon.ru
ORCID iD: 0000-0002-1960-6868

Elena R. Meskina — D. Sci. (Med.), Professor, Department of children's infections, Therapy department, 

Moscow

Russian Federation

A. Yu. Ushakova

M. Vladimirsky Moscow Regional Research Clinical Institute

Email: fake@neicon.ru
ORCID iD: 0000-0001-8438-7609

Anna Yu. Ushakova — Cand. Sci. (Med.), assistant, Department of children's infections, 

Moscow

Russian Federation

M. K. Khadisova

M. Vladimirsky Moscow Regional Research Clinical Institute

Email: fake@neicon.ru
ORCID iD: 0000-0001-8293-6643

Marima K. Khadisova — Cand. Sci. (Med.), Department of children's infections, Therapy department, 

Moscow

Russian Federation

V. V. Zverev

I.M. Sechenov First Moscow State Medical University (Sechenov University)

Email: fake@neicon.ru
ORCID iD: 0000-0002-0017-1892

Vitaly V. Zverev — D. Sci. (Med.), Professor, Academician of the Russian Academy of Sciences, Head, Department of microbiology, virology and immunology, 

Moscow

Russian Federation

E. B. Faizuloev

I. Mechnikov Research Institute of Vaccines and Sera

Email: fake@neicon.ru
ORCID iD: 0000-0001-7385-5083

Evgeny B. Faizuloev — Cand. Sci. (Biol.), Head, Laboratory of molecular virology,

Moscow

Russian Federation

References

  1. Troeger C., Khalil I.A., Rao P.C., Cao S., Blacker B.F., Ahmed T., et al. Rotavirus vaccination and the global burden of rotavirus diarrhea among children younger than 5 years. JAMA Pediatr. 2018; 172(10): 958–65. https://doi.org/10.1001/jamapediatrics.2018.1960
  2. Hoog M.L.A., Vesikari T., Giaquinto C., Huppertz H.I., Martinon-Torres F., Bruijning-Verhagen P. Report of the 5th European expert meeting on rotavirus vaccination (EEROVAC). Hum. Vaccin. Immunother. 2018; 14(4): 1027–34. https://doi.org/10.1080/21645515.2017.1412019
  3. Matthijnssens J., Ciarlet M., Rahman M., Attoui H., Banyai K., Estes M.K., et al. Recommendations for the classification of group A rotaviruses using all 11 genomic RNA segments. Arch. Virol. 2008; 153(8): 1621–9. https://doi.org/10.1007/s00705-008-0155-1
  4. WHO. Manual of rotavirus detection and characterization methods. Geneva; 2009. Available at: https://www.who.int/vaccines-documents
  5. Doro R., Laszlo B., Martella V., Leshem E., Gentsch J., Parashar U., et al. Review of global rotavirus strain prevalence data from six years post vaccine licensure surveillance: is there evidence of strain selection from vaccine pressure? Infect. Genet. Evol. 2014; 28: 446–61. https://doi.org/10.1016/j.meegid.2014.08.017
  6. Zhou X., Wang Y.H., Ghosh S., Tang W.F., Pang B.B., Liu M.Q., et al. Genomic characterization of G3P[6], G4P[6] and G4P[8] human rotaviruses from Wuhan, China: Evidence for interspecies transmission and reassortment events. Infect. Genet. Evol. 2015; 33: 55–71. https://doi.org/10.1016/j.meegid.2015.04.010
  7. Salamunova S., Jackova A., Csank T., Mandelik R., Novotny J., Beckova Z., et al. Genetic variability of pig and human rotavirus group A isolates from Slovakia. Arch. Virol. 2020; 165(2): 463–70. https://doi.org/10.1007/s00705-019-04504-6
  8. Doro R., Farkas S.L., Martella V., Banyai K. Zoonotic transmission of rotavirus: surveillance and control. Expert. Rev. Anti. Infect. Ther. 2015; 13(11): 1337–50. https://doi.org/10.1586/14787210.2015.1089171
  9. Velasquez D.E., Jiang B. Evolution of P[8], P[4], and P[6] VP8* genes of human rotaviruses globally reported during 1974 and 2017: possible implications for rotavirus vaccines in development. Hum. Vaccin. Immunother. 2019; 15(12): 3003–8. https://doi.org/10.1080/21645515.2019.1619400
  10. Hungerford D., Vivancos R., Read J.M., Iturriza-Gomicronmara M., French N., Cunliffe N.A. Rotavirus vaccine impact and socioeconomic deprivation: an interrupted time-series analysis of gastrointestinal disease outcomes across primary and secondary care in the UK. BMC Med. 2018; 16(1): 10. https://doi.org/10.1186/s12916-017-0989-z
  11. Freeman M.M., Kerin T., Hull J., McCaustland K., Gentsch J. Enhancement of detection and quantification of rotavirus in stool using a modified real-time RT-PCR assay. J. Med. Virol. 2008; 80(8): 1489–96. https://doi.org/10.1002/jmv.21228
  12. Kiseleva V., Faizuloev E., Meskina E., Marova A., Oksanich A., Samartseva T., et al. Molecular-genetic characterization of human rotavirus a strains circulating in Moscow, Russia (2009- 2014). Virol. Sin. 2018; 33(4): 304–13. https://doi.org/10.1007/s12250-018-0043-0
  13. Iturriza-Gomara M., Isherwood B., Desselberger U., Gray J. Reassortment in vivo: driving force for diversity of human rotavirus strains isolated in the United Kingdom between 1995 and 1999. J. Virol. 2001; 75(8): 3696–705. https://doi.org/10.1128/JVI.75.8.3696-3705.2001
  14. Simmonds M.K., Armah G., Asmah R., Banerjee I., Damanka S., Esona M., et al. New oligonucleotide primers for P-typing of rotavirus strains: Strategies for typing previously untypeable strains. J. Clin. Virol. 2008; 42(4): 368–73. https://doi.org/10.1016/j.jcv.2008.02.011
  15. Froussard P. rPCR: a powerful tool for random amplification of whole RNA sequences. PCR Methods Appl. 1993; 2(3): 185–90. https://doi.org/10.1101/gr.2.3.185
  16. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34(18): 3094–100. https://doi.org/10.1093/bioinformatics/bty191
  17. Maes P., Matthijnssens J., Rahman M., Van Ranst M. RotaC: a web-based tool for the complete genome classification of group A rotaviruses. BMC Microbiol. 2009; 9: 238. https://doi.org/10.1186/1471-2180-9-238
  18. Faizuloev E., Mintaev R., Petrusha O., Marova A., Smirnova D., Ammour Y., et al. New approach to genetic characterization of group A rotaviruses by the nanopore sequencing method. J. Virol. Methods. 2021; 292: 114114. https://doi.org/10.1016/j.jviromet.2021.114114
  19. Kumar S., Stecher G., Li M., Knyaz C., Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018; 35(6): 1547–9. https://doi.org/10.1093/molbev/msy096
  20. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980; 16(2): 111–20. https://doi.org/10.1007/BF01731581
  21. Cowley D., Donato C.M., Roczo-Farkas S., Kirkwood C.D. Emergence of a novel equine-like G3P[8] inter-genogroup reassortant rotavirus strain associated with gastroenteritis in Australian children. J. Gen. Virol. 2016; 97(2): 403–10. https://doi.org/10.1099/jgv.0.000352
  22. Arana A., Montes M., Jere K.C., Alkorta M., Iturriza-Gomara M., Cilla G. Emergence and spread of G3P[8] rotaviruses possessing an equine-like VP7 and a DS-1-like genetic backbone in the Basque Country (North of Spain), 2015. Infect. Genet. Evol. 2016; 44: 137–44. https://doi.org/10.1016/j.meegid.2016.06.048
  23. Doro R., Marton S., Bartokne A.H., Lengyel G., Agocs Z., Jakab F., et al. Equine-like G3 rotavirus in Hungary, 2015 — Is it a novel intergenogroup reassortant pandemic strain? Acta. Microbiol. Immunol. Hung. 2016; 63(2): 243–55. https://doi.org/10.1556/030.63.2016.2.8
  24. Guerra S.F.S., Soares L.S., Lobo P.S., Penha Junior E.T., Sousa Junior E.C., Bezerra D.A.M., et al. Detection of a novel equinelike G3 rotavirus associated with acute gastroenteritis in Brazil. J. Gen. Virol. 2016; 97(12): 3131–8. https://doi.org/10.1099/jgv.0.000626
  25. Komoto S., Tacharoenmuang R., Guntapong R., Ide T., Haga K., Katayama K., et al. Emergence and characterization of unusual DS-1-Like G1P[8] rotavirus strains in children with diarrhea in Thailand. PLoS One. 2015; 10(11): e0141739. https://doi.org/10.1371/journal.pone.0141739
  26. Kaira A.N., Fayzuloyev E.B., Lavrov V.F., Svitich O.A., Solomay T.V., Nikonova A.A., et al. Epidemiological trends of morbidity and issues of rotavirus vaccination at the present stage. Sanitary Doctor. 2020; 6: 16–25. https://doi.org/0.33920/med-08-2006-02
  27. Ivashechkin A.A., Yuzhakov A.G., Grebennikova T.V., Yuzhakova K.A., Kulikova N.Y., Kisteneva L.B., et al. Genetic diversity of group A rotaviruses in Moscow in 2018-2019. Arch. Virol. 2020; 165(3): 691–702. https://doi.org/10.1007/s00705-020-04534-5
  28. Yuzhakov G., Yuzhakova K., Kulikova N., Kisteneva L., Cherepushkin S., Smetanina S., et al. Prevalence and genetic diversity of group a rotavirus genotypes in Moscow (2019–2020). Pathogens. 2021; 10: 674. https://doi.org/10.3390/pathogens10060674
  29. Sashina T.A., Morozova O.V., Epifanova N.V., Novikova N.A. Predominance of new G9P[8] rotaviruses closely related to Turkish strains in Nizhny Novgorod (Russia). Arch. Virol. 2017; 162(8): 2387–92. https://doi.org/10.1007/s00705-017-3364-7
  30. Новикова Н.А., Сашина Т.A., Солнцев Л.А., Епифанова Н.В., Кашников A.Ю., Погодина Л.В. и др. Проявление эпидемического процесса ротавирусного процесса в Нижнем Новгороде в предвакцинальный период. Журнал микробиологии, эпидемиологии и иммунобиологии. 2017; 94(5): 46–52. https://doi.org/10.36233/0372-9311-2017-5-46-52
  31. Денисюк Н.Б. Генетическая характеристика ротавирусов группы А, циркулирующих в Оренбургском регионе в сезон 2016–2017 гг. Детские инфекции. 2017; 16(4): 42–5. https://doi.org/10.22627/2072-8107-2017-16-4-42-45
  32. Мескина Е.Р., Ушакова А.Ю., Файзулоев Е.Б., Бахтояров Г.Н., Киселева В.В. Сравнительная характеристика гастроэнтерита, вызванного ротавирусами генотипов G4P[8] и G9P[8], у детей, госпитализированных в стационар г. Москвы (эпидсезон 2012–2013 гг.). Инфекционные болезни. 2017; 15(1): 23–8. https://doi.org/10.20953/1729-9225-2017-1-23-28
  33. Katz E.M., Esona M.D., Betrapally N.S., De La Cruz De Leon L.A., Neira Y.R., Rey G.J., et al. Whole-gene analysis of inter-genogroup reassortant rotaviruses from the Dominican Republic: Emergence of equine-like G3 strains and evidence of their reassortment with locally-circulating strains. Virology. 2019; 534: 114–31. https://doi.org/10.1016/j.virol.2019.06.007
  34. Fujii Y., Oda M., Somura Y., Shinkai T. Molecular characteristics of novel mono-reassortant G9P[8] rotavirus a strains possessing the NSP4 gene of the E2 genotype detected in Tokyo, Japan. Jpn J. Infect. Dis. 2020; 73(1): 26–35. https://doi.org/10.7883/yoken.JJID.2019.211
  35. McDonald S.M., Matthijnssens J., McAllen J.K., Hine E., Overton L., Wang S., et al. Evolutionary dynamics of human rotaviruses: balancing reassortment with preferred genome constellations. PLoS Pathog. 2009; 5(10): e1000634. https://doi.org/10.1371/journal.ppat.1000634
  36. Papp H., Borzak R., Farkas S., Kisfali P., Lengyel G., Molnar P., et al. Zoonotic transmission of reassortant porcine G4P[6] rotaviruses in Hungarian pediatric patients identified sporadically over a 15 year period. Infect. Genet. Evol. 2013; 19: 71–80. https://doi.org/10.1016/j.meegid.2013.06.013
  37. Esona M.D., Geyer A., Banyai K., Page N., Aminu M., Armah G.E., et al. Novel human rotavirus genotype G5P[7] from child with diarrhea, Cameroon. Emerg. Infect. Dis. 2009; 15(1): 83–6. https://doi.org/10.3201/eid1501.080899
  38. Wang Y.H., Kobayashi N., Nagashima S., Zhou X., Ghosh S., Peng J.S., et al. Full genomic analysis of a porcine-bovine reassortant G4P[6] rotavirus strain R479 isolated from an infant in China. J. Med. Virol. 2010; 82(6): 1094–102. https://doi.org/10.1002/jmv.21760
  39. Zeller M., Patton J.T., Heylen E., De Coster S., Ciarlet M., Van Ranst M., et al. Genetic analyses reveal differences in the VP7 and VP4 antigenic epitopes between human rotaviruses circulating in Belgium and rotaviruses in Rotarix and RotaTeq. J. Clin. Microbiol. 2012; 50(3): 966–76. https://doi.org/10.1128/JCM.05590-11
  40. Imagawa T., Saito M., Yamamoto D., Saito-Obata M., Masago Y., Ablola A.C., et al. Genetic diversity of species A rotaviruses detected in clinical and environmental samples, including porcine-like rotaviruses from hospitalized children in the Philippines. Infect. Genet. Evol. 2020; 85: 104465. https://doi.org/10.1016/j.meegid.2020.104465

Supplementary files

Supplementary Files
Action
1. JATS XML
Download
2. Fig. 1. Distribution of RVA G/[P]-genotypes isolated from 131 clinical samples in 2015–2020. GxP[8]Ix — partially typed PVA; Mix — mixed genotypes; Nt — non-typed samples.

Download (69KB)
Indexing metadata
3. Fig. 2. Prevalence of G/[P]-genotypes detected in Moscow during the 10-year monitoring period (2009–2020) [12]. Nt — non-typed samples.

Download (48KB)
Indexing metadata
4. Fig. 3. Phylogenetic trees based on the sequenced genes of the VP4 and VP7 proteins of the new G3P[8]I2 strain. ▲ — segments sequenced using the Sanger technique; ♦ — the gene sequenced using the NPS technique; ∆ — the equine RVA VP7 gene. The respective GenBank accession number, name of strains and G/[P]-genotype were used for designation of strains.

Download (704KB)
Indexing metadata
5. Fig. 4. Phylogenetic trees based on sequenced genes of RVA VP4, VP7, and VP6 proteins, strain G4P[6]I1 (marked by ♦), reference strains (Wa, AU-1, DS-1) representatives of three evolutionary lines RVA (marked by ○) and RVA genes of porcine origin, phylogenetically most closely related to genes of the Moscow-1P/2015 strain based on the BLAST analysis. The respective GenBank number, name of the strain and G/[P]-genotype were used for designation of strains.

Download (189KB)
Indexing metadata

Copyright (c) 2022 Petrusha O.A., Korchevaya E.R., Mintaev R.R., Nikonova A.A., Isakov I.Y., Meskina E.R., Ushakova A.Y., Khadisova M.K., Zverev V.V., Faizuloev E.B.

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.
СМИ зарегистрировано Федеральной службой по надзору в сфере связи, информационных технологий и массовых коммуникаций (Роскомнадзор).
Регистрационный номер и дата принятия решения о регистрации СМИ: ПИ № ФС77-75442 от 01.04.2019 г.


This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies