Diversity of group A rotavirus genes detected in the Triângulo Mineiro region, Minas Gerais, Brazil

ABSTRACT

Group A rotaviruses are the main causative agent of infantile gastroenteritis. The segmented nature of the viral genome allows reassortment of genome segments, which can generate genetic variants. In this study, we characterized the diversity of the VP7, VP4 (VP8*), VP6, NSP4, and NSP5 genes of the rotaviruses that circulated from 2005 to 2011 in the Triângulo Mineiro (TM) region of Brazil. Samples with genotypes G2 (sublineages IVa-1 and IVa-3), G1 (sublineage I-A), G9 (lineage III), G12 (lineages II and III), G8 (lineage II), G3 (lineage III), P[4] (sublineages IVa and IVb), P[8] (sublineages P[8]-3.6, P[8]-3.3, and P[8]-3.1), I2 (lineage VII), E2 (lineages VI, XII, and X), and H2 (lineage III) were identified. The associations found in the samples were G1, G9, or G12 with P[8]-I1-E1-H1; G2 or G8 with P[4]-I2-E2-H2; G12 with I3-E3-H6; and G3 with P[4]-I2-E3-H3 (previously unreported combination). Reassortment events in G2P[4] strains and an apparent pattern of temporal segregation within the lineages were observed. Five TM samples contained genes that exhibited high nucleotide and amino acid identities with strains of animal origin. The present study includes a period of pre- and post-introduction of rotavirus vaccination in all Brazilian territories, thereby serving as a basis for monitoring changes in the genetic constitution of rotaviruses. The results also contribute to the understanding of the diversity and evolution of rotaviruses in a global context.

Keywords:
Rotavirus; Infantile diarrhea; Genotype; Genetic variation

Introduction

The group A rotavirus (RVA), a member of the genus Rotavirus in the family Reoviridae, is the main causative agent of severe diarrheal disease in children. In 2008, the virus was associated with the deaths of 453,000 children less than 5 years of age worldwide, with most cases reported in Africa and Asia.¹

The viral particle is non-enveloped and composed of 11 segments of double-stranded RNA surrounded by a triple-layered protein capsid. The viral genome encode six structural (VP1-VP4, VP6 and VP7) and six non-structural (NSP1-NSP6) proteins. The outer capsid proteins, VP7 and VP4 (cleaved by trypsin into VP5* and VP8*), are involved in the attachment and entry during viral replication in host cells, and each can independently induce the production of neutralizing antibodies. The serotypes defined by VP7 and VP4, as well as their encoding genes, are used to classify rotaviruses into G and P genotypes, respectively, based on a widely used binomial system.² Previous studies have demonstrated that RVAs of genotypes G1, G2, G3, G4, G9, and G12, combined with genotypes P[4] or P[8], are most prevalent among the strains that cause human infection.³

Events such as point mutations, rearrangements, interspecies transmissions, and reassortments occurring during co-infections with different strains can generate genetic variants with vast inter- and intragenotypic diversity. This high variability and observed fluctuations in the prevalence of dominant genotypes throughout epidemic seasons represent a major challenge for immunization efforts using vaccines currently available for the prevention of severe infections caused by RVAs in many countries. To better characterize the diversity and evolutionary aspects of rotaviruses, a new classification system based on the 11 segments of the RVA genome has recently been proposed, which emphasizes the importance of analyzing other genes in addition to those that encode for the VP7 and VP4 proteins.⁴ In this latter classification, the VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5/6-encoding gene segments are classified as Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx (x indicates the number of the genotype), respectively. Genes 6, 10, and 11 encode the VP6, NSP4, and NSP5 proteins, respectively. The VP6 protein has been studied as a possible target for vaccine development, and studies in animals have demonstrated protective immune responses induced by VP6 protein immunization.⁵ NSP4 is involved in viral morphogenesis and has enterotoxin-like activity.⁶ Last, NSP5 is essential for viroplasm formation, genome replication, and virion assembly, and its gene is most frequently associated with rearrangement events.²

Since 2005, surveillance of circulating RVA strains has been conducted in the Triângulo Mineiro (TM) region of Minas Gerais State, Brazil. The studies showed that different G-types associated with P[8] co-circulated in 2005 and also that G2P[4] had emerged by mid-2006, continuing as the prevalent genotype until 2010.^7-9 In this study, we characterized the diversity of the VP7, VP4 (VP8* subunit), VP6, NSP4, and NSP5 genes of representative RVA strains that circulated during six seasons (2005-2011) in the TM region.

Materials and methods

Sample collection

Twenty-seven representative RVA strains from a collection of 117 rotavirus-positive fecal samples from Uberaba (2005-2009 and 2011) and Uberlândia (2008-2011) were selected based on their electropherotypes and VP7 and VP4 genotypes, previously characterized in three-laboratory-based surveillance studies (Table 1).^7-9 These previous studies analyzed fecal specimens from children who presented acute gastroenteritis, which were provided by laboratories Jorge Furtado Medicina Diagnóstica (Uberaba) or IPAC Medicina Diagnóstica (Uberlândia). Information on the vaccination history of the patients was not available. The 27 representative samples were collected from outpatients and inpatients from one hospital in Uberaba and two hospitals in Uberlândia. This study was approved by the Human Research Ethics Committee of Universidade Federal do Triângulo Mineiro (protocol CEP/UFTM/672/06).

Table 1
Selected RVA G/P types collected from children with gastroenteritis in Uberaba and Uberlândia, Triângulo Mineiro, Brazil, from 2005 to 2011.

G/P type	City	2005–2006	2007	2008	2009	2010	2011
		Collected/selected	Collected/selected	Collected/selected	Collected/selected	Collected/selected	Collected/selected
G1P[8]	UB	3/1	–	–	–	–	–
G2P[4]	UB	25/4	33/3	5/2	–	–	–
	IP	–	–	19/3	–	16/3	–
G9P[8]	UB	7/4	–	–	–	–	–
G3P[4]a	IP	–	–	–	–	–	1/1
GNTP[8]	UB	2/2	–	–	–	–	–
G3P[NT]	IP	–	–	–	1/1	–	–
GNTP[NT]	UB	2/2	–	–	–	–	1/1
Mixed	UB	–	–	–	1/0	–	–
	IP	–	–	1/0	–	–	–
Total		39/13	33/3	25/5	2/1	16/3	2/2

• ^a
  Typed as G3 by PCR but determined to be G8 by sequence analysis.
• UB, Uberaba; IP, Uberlândia; NT, non-typed.

RNA extraction, reverse transcription polymerase chain reaction (RT-PCR), and nucleotide (nt) sequencing

The fecal specimens were subjected to nucleic acid extraction,¹⁰ and reverse transcription was carried out using random primers (Invitrogen, Carlsbad, CA, USA) and murine leukemia virus reverse transcriptase (Promega, Madison, WI, USA) according to the manufacturers’ instructions. The complementary DNA generated by the RT reaction was then used for the amplification of all genes. The primers and reaction conditions for amplification of each target gene were based on previously described protocols as follows: primers VP7-F and VP7-R for partial (nt positions 49 to 933) VP7 gene¹¹; VP4-F and VP4-R for the VP8* subunit (nt positions 132 to 795) of the VP4 gene¹²; JRG30(+) and JRG31(−) for a full-length gene (nt positions 1 to 738) encoding the NSP4 protein¹³; VP6-F and VP6-R for a partial gene (nt positions 747 to 1126) encoding the VP6 protein¹⁴; and GEN-NSP5-F and GEN-NSP5-R for a full-length gene (nt positions 1 to 667 or to 1043) encoding the NSP5 protein.¹⁵ PCR products were purified using the EasyPrep Gel/PCR Purification Mini Kit (EasyPath, São Paulo, SP, Brazil), separated by agarose gel electrophoresis, and then quantified by comparing the band intensities with those of the Low DNA Mass Ladder (Invitrogen, Carlsbad, CA, USA). The DNA fragments were then sequenced in both directions with the same primers as used for the PCRs using the BigDye^® Terminator Cycle Sequencing Kit v. 3.1 and an ABI 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA).

Genotype assignment and phylogenetic analysis

The VP7, VP8*, VP6, NSP4, and NSP5 gene sequences were curated with the Geneious software (Biomatters Ltd., Auckland, New Zealand). To identify the genotypes and evaluate the genetic diversity, the TM samples were compared with the reference strains representing distinct genotypes, lineages, and sublineages and with most similar field samples from different continents, as well as with vaccine strains. The lineages and sublineages of G2, P[4], I2, E2, and H2 samples were previously suggested by Giammanco et al.,¹⁶ whereas the G1, G9, G12, G8, G3, and P[8] lineages and sublineages were described by Phan et al.,^17,18 Rahman et al.,¹⁹ Ianiro et al.,²⁰ Ndze et al.,²¹ and da Silva et al.,²² respectively. All sequences of the field and reference strains were retrieved from GenBank. Alignments were performed using MUSCLE (Multiple Sequence Comparison by Log-Expectation) and analyzed with the maximum likelihood algorithm.^23,24 The best-fit model of sequence evolution was estimated by jModelTest according to the Bayesian Information Criterion, and phylogenetic trees were constructed for each gene. The maximum-likelihood trees were generated using selected models of nt substitution, GTR+I+G (VP7), HKY+G (VP8* and NSP5), HKY+G+I (VP6), and TrN+G (NSP4), in the Molecular Evolutionary Genetics Analysis 5.1 software (MEGA Team, Japan). The nomenclature of all rotavirus strains in the trees is listed according to Matthijnssens et al.²⁵ The UB, IP and MG abbreviations in the nomenclature of the TM samples represent Uberaba, Uberlândia, and Minas Gerais State, respectively.

Nucleotide sequence accession numbers

The nucleotide sequences were deposited in the GenBank database, and the accession numbers are: KF938698-KF938771 for the G2P[4] samples and KR018449-KR018504 for the other G-types. The accession numbers of the strains retrieved from GenBank are shown at the end of each name in the trees.

Results

A total of 130 nt consensus sequences were obtained for the 27 representative TM samples. In four of the samples, at least one of the five target genes was not amplified with the primers used.

Based on the VP7 gene, six genotypes were characterized: G1, G2, G3, G8, G9, and G12 (Fig. 1a). The TM G2 samples shared 95.3-100% and 96.8-100% nt and amino acid (aa) identity, respectively. They belonged to lineage IV, segregating into sublineages IVa-3 and IVa-1. The TM G1 sample belonged to lineage I, sublineage I-A, and shared 94% nt and 95.3% aa identity with the Rotarix™ vaccine strain of lineage II. The TM G9 samples belonged to lineage III and shared 97.3-100% nt and aa identity. The TM G12 samples shared 96.4-100% and 96.7-100% nt and aa identity, respectively. These samples clustered into lineages II (with strains of G12P[9] combination) and III (with strains of G12P[8] and G12P[6] combinations). The TM G8 sample belonged to lineage II; it exhibited high nt (99.3-99.9%) and aa (99.1-100%) identities with other contemporary human P[4] samples from Italy (SS56), USA (USRo104), Germany (GER1H-09), and Brazil (RN19461-08). The TM G3 sample fell into lineage III, in a separate sub-cluster from the prototype strain (AU-1). They only shared 90.7% nt and 94.9% aa identity, but the TM G3 sample showed higher identity (nt, 99.1%; aa, 98.4%) with a cat sample (BA222) isolated from Italy.

Fig. 1
Maximum-likelihood phylogenetic analysis based on the VP7 (a), VP8* (b), VP6 (c), NSP4 (d), and NSP5 (e) gene nucleotide sequences of Triângulo Mineiro (TM) strains and the sequences from the GenBank database. The numbers at the nodes indicate bootstrap values (1000 replicates), and only values above 70% are shown. The trees are drawn to scale, with the branch lengths measured as the number of substitutions per site. The TM samples are marked with a solid circle, and the reference strains of each lineage and sublineage are marked with an open triangle. The uppercase letters followed by Arabic numerals, Roman numeral, and lowercase letters represent different genotypes, lineages, and sublineages, respectively, except for lineages and sublineages of the P[8] genotype, which are indicated by Arabic numerals, and sublineages of G1 strains, which are indicated by uppercase letters.

Based on the VP8* gene, the TM samples were distributed among P[8]- and P[4]-types (Fig. 1b). The TM P[4] samples shared 96.2-100% and 95.2-100% nt and aa identity, respectively. These samples belonged to lineage IV and split into sublineages IVa and IVb. The TM P[8] samples showed 96.5-100% nt and 95.6-100% aa identity. They belonged to lineage P[8]-3 and were divided between sublineages P[8]-3.1, P[8]-3.3, and P[8]-3.6, regardless of the combinations with the G genotypes.

Based on the VP6 gene, the TM samples belonged to I1, I2, and I3 genotypes (Fig. 1c). The TM I2 strains showed 95.9-100% and 99.1-100% nt and aa identity, respectively, and clustered together with other strains identified within lineage VII worldwide. All I2 strains shared P[4] genotype specificity and formed two main sub-clusters (a and b) with their separation supported by a high bootstrap value (100%). The TM samples that grouped into cluster I3 exhibited G12 specificity and showed high nt (99.4-100%) and aa (100%) identities to the reference sample (T152) of the G12P[9] combination. The TM I1 samples shared an nt identity of 94.7-100% and an aa identity of 98.2-100%. They grouped apart from the prototype (Wa) and vaccine-derived (Rotarix™) strains. The TM I1 samples shared 87.3-100% nt and 95.6-100% aa identity with these reference strains, and the majority shared P[8] genotype specificity.

Based on the NSP4 gene, the TM samples were distributed among the E2, E1, and E3 genotypes (Fig. 1d). The TM E2 samples shared 85.0-100% nt and 94.0-100% aa identity. They were subdivided between the clusters of lineages VI, X, and XII, composed of strains with P[4] specificity. The TM samples of lineage VI formed two phylogenetic sub-clusters (a and b), with their separation supported by a high bootstrap value (99%). TM strains UB-239 and UB-226 displayed high nt (95.7-100%) and aa (97.7-100%) identities with the bovine 1603 strain detected in 2007 in South Africa²⁶ and with a human Brazilian strain, 12389_06RJ (99.0-100% nt and 99.4-100% aa identity), detected in 2006.²⁷ These four strains did not group with any cluster of the previously established lineages.¹⁶ The TM E1 samples shared 91.7-100% and 92.8-100% nt and aa identity, respectively, and the majority of the strains had the P[8] genotype. The TM E3 samples showed 88.7% and 94.0% nt and aa identity, respectively. They grouped with strains that had genotype P[9] combined with G3, G6, or G12.

Based on the NSP5 gene, the TM samples were distributed among the H1, H6, H3, and H2 genotypes (Fig. 1e). The TM H1 samples shared nt and aa identities in the ranges of 97.3-100% and 99.0-100%, respectively. They grouped with field samples of P[8] specificity and clustered apart from the prototypes, Wa and KU, and the vaccine-derived (Rotarix™) strain. The nt identity with these reference strains was 92.3-100%, and the aa identity was 93.9-100%. The TM H6 samples were identical and clustered together with samples exhibiting the G12P[9] combination. They were closely related (nt identity, 94.3%; aa identity, 98.5%) to a simian RVA strain (RRV). The TM H3 sample (IP-218) grouped apart from the prototype (AU-1) and vaccine (RotaTeq™) strains, with which it showed the nt identity of 93.3-95.2% and the aa identity of 92.9-96.0%. At the same time, sample IP-218 showed high nt (98.7%) and aa (97.5%) identities with a cat sample (BA222) of G3P[9] specificity. The TM H2 samples shared the nt identity of 97.2-99.9% and the aa identity of 98.5-100%. They grouped with P[4] strains into two main sub-clusters (a and b) within lineage III.

The classification of the TM samples by genotypes, lineages, and sublineages on the basis of phylogenetic analysis of the five studied genes is summarized in Table 2. When analyzing the TM samples according to their genotypes, four types of combinations (or genetic constellations) could be observed. The first one consisted of samples that exhibited genotypes G1, G9, or G12 associated with P[8]-I1-E1-H1 and had long electropherotypes. The second one was composed of samples that showed genotypes G2 or G8 associated with P[4]-I2-E2-H2 and had short electropherotypes. The third one included samples exhibiting the G12 genotype combined with I3-E3-H6 and having long electropherotypes, and the fourth one was formed by a sample with genotype G3 associated with P[4]-I2-E3-H3 and having a long electropherotype (Table 2).

Table 2
Epidemiological data, genotypes, lineages, and sublineages of Triângulo Mineiro rotavirus A strains, based on the VP7, VP8*, VP6, NSP4, and NSP5 genes.

Name of the strainsa	Patients age	Origin	Month-year of collection	Electropherotype	Genotypes						Lineages and sublineages
					VP7	VP8*	VP6	NSP4	NSP5		VP7	VP8*	VP6	NSP4	NSP5
Uberaba
RVA/Human-wt/BRA/UB-09 MG/2005/G1P[8]	1 year	Inpatient (HC)	Sep-2005	Long	G1	P[8]	I1	E1	H1		I-A	3.1	–	–	–
RVA/Human-wt/BRA/UB-16 MG/2005/G9P[8]	2 years	Inpatient (HC)	Sep-2005	Long	G9	P[8]	I1	E1	H1		III	3.1	–	–	–
RVA/Human-wt/BRA/UB-142 MG/2006/GxP[8]	2 years	Outpatient	Jul-2006	Long	Gx	P[8]	I1	E1	H1		–	3.1	–	–	–
RVA/Human-wt/BRA/UB-144 MG/2006/G12Px	1 year	Outpatient	Aug-2006	Long	G12	Px	I3	E3	H6		II	–	–	–	–
RVA/Human-wt/BRA/UB-146 MG/2006/G12P[8]	11 months	Outpatient	Aug-2006	Long	G12	P[8]	I1	E1	H1		III	3.1	–	–	–
RVA/Human-wt/BRA/UB-150 MG/2006/G12Px	4 years	Outpatient	Aug-2006	Long	G12	Px	I3	Ex	H6		II	–	–	–	–
RVA/Human-wt/BRA/UB-157 MG/2006/G9P[8]	1 year	Outpatient	Aug-2006	Long	G9	P[8]	I1	E1	H1		III	3.1	–	–	–
RVA/Human-wt/BRA/UB-169 MG/2006/G9P[8]	1 year	Inpatient (HC)	Aug-2006	Long	G9	P[8]	I1	E1	H1		III	3.6	–	–	–
RVA/Human-wt/BRA/UB-208 MG/2006/G9P[8]	1 year	Outpatient	Sep-2006	Long	G9	P[8]	I1	E1	H1		III	3.1	–	–	–
RVA/Human-wt/BRA/UB-558 MG/2011/G9P[8]	23 months	Inpatient (HC)	Dec-2011	Long	G9	P[8]	I1	E1	H1		III	3.3	–	–	–
RVA/Human-wt/BRA/UB-154 MG/2006/G2P[4]	3 years	Outpatient	Aug-2006	Short	G2	P[4]	I2	E2	H2		IVa-1	IVa	VIIa	VIa	IIIa
RVA/Human-wt/BRA/UB-170 MG/2006/G2P[4]	8 months	Inpatient (HC)	Aug-2006	Short	G2	P[4]	I2	E2	H2		IVa-1	IVa	VIIa	VIa	IIIa
RVA/Human-wt/BRA/UB-226 MG/2006/G2P[4]	5 months	Outpatient	Sep-2006	Short	G2	P[4]	I2	E2	H2		IVa-1	IVa	VIIa	–	IIIa
RVA/Human-wt/BRA/UB-239 MG/2006/G2P[4]	1 year	Outpatient	Oct-2006	Short	G2	P[4]	I2	E2	H2		IVa-1	IVa	VIIa	–	IIIa
RVA/Human-wt/BRA/UB-299 MG/2007/G2P[4]	1 year	Inpatient (HC)	Jun-2007	Short	G2	P[4]	I2	E2	H2		IVa-1	IVa	VIIa	VIa	IIIa
RVA/Human-wt/BRA/UB-320 MG/2007/G2P[4]	11 months	Inpatient (HC)	Aug-2007	Short	G2	P[4]	I2	E2	Hx		IVa-3	IVa	VIIb	VIb	–
RVA/Human-wt/BRA/UB-354 MG/2007/G2P[4]	1 year	Inpatient (HC)	Oct-2007	Short	G2	P[4]	I2	E2	H2		IVa-3	IVa	VIIb	VIb	IIIa
RVA/Human-wt/BRA/UB-432 MG/2008/G2P[4]	1 year	Inpatient (HC)	Sep-2008	Short	G2	P[4]	I2	E2	H2		IVa-3	IVa	VIIb	VIb	IIIa
RVA/Human-wt/BRA/UB-433 MG/2008/G2P[4]	4 years	Outpatient	Sep-2008	Short	G2	P[4]	I2	E2	H2		IVa-3	IVb	VIIb	X	IIIb
Uberlândia
RVA/Human-wt/BRA/IP-38 MG/2008/G2P[4]	3 years	Inpatient (HSG)	Aug-2008	Short	G2	P[4]	I2	E2	H2		IVa-3	IVb	VIIb	XII	IIIb
RVA/Human-wt/BRA/IP-47 MG/2008/G2P[4]	3 years	Inpatient (HM)	Sep-2008	Short	G2	P[4]	I2	E2	H2		IVa-3	IVb	VIIb	X	IIIb
RVA/Human-wt/BRA/IP-71 MG/2008/G2P[4]	2 years	Inpatient (HSG)	Nov-2008	Short	G2	P[4]	I2	E2	H2		IVa-3	IVa	VIIa	VIb	IIIa
RVA/Human-wt/BRA/IP-346 MG/2010/G2P[4]	1 year	Inpatient (HM)	May-2010	Short	G2	P[4]	I2	E2	H2		IVa-3	IVb	VIIb	VIb	IIIb
RVA/Human-wt/BRA/IP-364 MG/2010/G2P[4]	2 years	Inpatient (HSG)	Jul-2010	Short	G2	P[4]	I2	E2	H2		IVa-1	IVb	VIIb	X	IIIb
RVA/Human-wt/BRA/IP-393 MG/2010/G2P[4]	6 months	Outpatient	Sep-2010	Short	G2	P[4]	I2	E2	H2		IVa-3	IVb	VIIb	X	IIIb
RVA/Human-wt/BRA/IP-218 MG/2009/G3P[4]	4 years	Outpatient	Sep-2009	Long	G3	P[4]	I2	E3	H3		III	IVa	VIIb	–	–
RVA/Human-wt/BRA/IP-447 MG/2011/G8P[4]	5 years	Outpatient	Jan-2011	Short	G8	P[4]	I2	E2	H2		II	IVa	VIIa	X	IIIa

• ^a
  Name of RV strains according to the nomenclature proposed by Matthijnssens et al.²⁵
• UB, Uberaba; IP, Uberlândia; HC, Hospital da Criança; HSG, Hospital Santa Genoveva; HM, Hospital Madrecor; x, was not possible to amplify.

According to the clustering patterns of the TM G2P[4] samples in the dendrograms for different genes, shown in Fig. 1, some samples always grouped together within the same clusters. However, other samples changed their grouping patterns, clustering with different samples depending on the gene, which is suggestive of intra-genotypic reassortment in particular genes. This was observed for strains IP-364 (VP7 gene), UB-354, UB-432, UB-320, and IP-71 (VP8* gene), IP-71 and IP-393 (VP6 gene), IP-346, UB-239, and UB-226 (NSP4 gene), and UB-354, UB-432, and IP-71 (NSP5 gene). Strains UB-154, UB-170, UB-299, IP-47, UB-433, and IP-38 are probably non-reassortant strains, at least with regard to the genes studied.

Discussion

In the TM samples, genotype G2 was found to be always associated with the P[4], I2, E2, and H2 genotypes, whereas the G1 and G9 genotypes were combined with the P[8]-I1-E1-H1 genotypes (Table 2). The genotype combination G2-P[4]-I2-E2-H2 is characteristic of the DS-1-like genotype constellation, while G1-P[8]-I1-E1-H1 and G9-P[8]-I1-E1-H1 are typical combinations of the Wa-like genotype constellation. These two genotype constellations are found most commonly worldwide.²⁸

Despite having a typical DS-1-like constellation for all genes, the TM G2P[4] samples were shown to belong to lineages distinct from that of the DS-1 prototype strain, as already reported for contemporary samples.^16,27,29 Based on the VP7 gene, the G2P[4] samples were detected within sublineage IVa (Fig. 1a). This sublineage emerged in the 1990s and rapidly became the most common sublineage found around the world. The spread of this sublineage was associated with an amino acid substitution from aspartic acid (found in the DS-1 prototype strain) to asparagine at residue 96 (D96N) and with the occurrence of reassortment events in G2 strains.²⁹ These two mechanisms of variability were also observed in this study (aa data not shown).

Sublineage G2P[4]-IVa-1 of VP7 was found from August 2006 to June 2007 in samples from hospitalized patients and outpatients from Uberaba, followed by the detection of sublineage IVa-3 since August 2007 (Table 2). In Uberlândia, all samples belonged to sublineage IVa-3, except IP-364. Sublineages IVa and IVb of the VP8* gene also seemed to have emerged consecutively in Uberaba. This apparent pattern of temporal segregation of the sublineages of the VP7 and VP8* genes has also been described in other regions of Brazil and for Italian samples.^16,27 An apparent temporal segregation within the G2P[4] lineages from Uberaba was also observed for the VP6, NSP4, and NSP5 genes. Therefore, based on the analyses of the five genes, it seems that a new genetic variant exhibiting the combination of sublineages IVa-3 (VP7), IVb (VP8*), VIIb (VP6), X (NSP4), and IIIb (NSP5) began circulating in Uberaba in September 2008 (UB-433) and was also detected in September 2008 (IP-47) and in September 2010 (IP-393) in Uberlândia, corroborating the results obtained by Gómez et al. who characterized circulation of a new lineage in the southeast region of Brazil.²⁷

Samples UB-226 and UB-239 shared high identity in the NSP4 gene with a cow strain (1603). This result is in agreement with the findings of Matthijnssens et al. who revealed a common origin between human DS-1-like and bovine strains.³⁰

The G8 genotype is believed to be of bovine origin, but strains with this genotype were also isolated from several other hosts, including humans, and have been detected sporadically in Brazil.³ The TM G8 sample (IP-447) is closely related to other human RVA strains, such as SS56 from Italy²⁰ and GER1H-09 from Germany,³¹ sharing the same G8-II lineage. There is evidence that strains of this lineage do not represent a result of recent interspecies transmission,^20,31 contrary to what was seen with some G8 strains identified in Brazilian Indian children.³² The TM G8 sample (IP-447) was genotyped as P[4]-I2-E2-H2, and based on the five genes it was associated with sample GER1H-09, which was described as probably representing an intergenotypic reassortant sample between the G8 genotype and human DS-1-like rotaviruses.³¹

Analysis of the VP7 gene showed that the TM G1P[8] sample (UB-09) detected in 2005, i.e., before the Rotarix™ introduction, belonged to lineage I, which, along with lineage II, represents the most prevalent lineages since the 1990s in Brazil. Through phylogenetic analysis of the VP8* gene, it was possible to verify that UB-09 belonged to lineage P[8]-3, which is believed to be the only lineage currently circulating in Brazil.³³ This sample fell into sublineage P[8]-3.1, which has been shown to be most prevalent in Brazil in recent years.²²

The phylogenetic analysis of the VP7 gene from the TM G9 strains revealed that they belonged to lineage III and were closely related to strains detected on at least three continents. These results are in agreement with previous findings that showed the worldwide spread of this lineage.³⁴ When analyzing the VP8* gene, circulation of three sublineages (P[8]-3.6, P[8]-3.3, and P[8]-3.1) was observed in the TM region. Sublineage P[8]-3.6 was only observed in Brazil after the introduction of the Rotarix™ vaccine.²²

Genotype G12 was first detected in the fecal matter of children from Philippines in 1987. It had not been reported for about 10 years and then rapidly spread globally. It is believed that lineage III was the main lineage responsible for the spread of G12 worldwide.³⁴ This lineage has been reported to be found most often in combination with the P[6] or P[8] genotypes.²⁸ One TM G12 sample (UB-146) was found to belong to this lineage and to have the P[8]-I1-E1-H1 genotype combination. This allele combination (Wa-like) is most commonly found for G12, but this genotype can also be found, although at a much lower frequency, associated with P[9]-I3-E3-H3, which are AU-1-like genotypes.²⁸ Two other TM G12 samples (UB-144 and UB-150) presented AU-1-like genotypes. They also exhibited the H6 genotype, which showed a high nt identity with simian RVA; therefore, it seems that the TM G12 samples have different origins, as has also been observed in another study from Brazil.³⁵ TM samples UB-144 and UB-150 were closely related to reference strain T152 based on all analyzed genes, and there is evidence that the NSP5 gene of this sample is of animal origin.¹⁹

According to the phylogenetic trees, the G12-Px-I3-E3/x-H6 strains probably possess P[9] specificity. The number of non-typable RVA strains depends on the number of genotype-specific primers used for RT-PCR genotyping. The genotyping failure is relatively high, given today's technology, and may represent a low viral load or unusual genotypes (e.g., originating from zoonotic transmission) that are not able to be detected with the primers used.³⁶ Furthermore, the PCR primers may have become outdated due to mutations in the primer-binding sites during co-circulation of wild-type RVAs over the years.³⁷

RVAs of the G3 genotype have several hosts and are commonly found associated with the P[8] genotype in humans. The G3 sample of long electropherotype (IP-218) associated with the P[4] genotype was once detected in the region.⁸ Although this combination is unusual, the G3P[4] genotype was responsible for an outbreak of diarrhea in 1998 in Juiz de Fora, Brazil.³⁸ In the IP-218 sample, a combination with the I2-E3-H3 genotype was found. The evidence of reassortment between G3 samples and typical G2P[4] strains (P[4]-I2-E2-H2) has already been provided³⁹; however, to our knowledge, this is the first time when the combination G3-P[4]-I2-E3-H3 is reported.

Based on the VP7 and NSP5 genes, G3 sample IP-218 fell into the cluster of the prototype strain (AU-1) with an RVA strain isolated from cat (BA222). Based on the VP8* and VP6 genes, it was closely related to human rotaviruses that share DS-1-like genotypes. Therefore, the IP-218 sample probably represents a reassortment between AU-1-like and DS-1-like strains, and it would be interesting to characterize the other genome segments of this sample. A close relationship between human G3 RVAs and feline strains has previously been described in Brazil^40,41 and other countries.⁴²

The phylogenetic analyses of the five genes from the G2P[4] strains showed that various reassortment events likely occurred from 2006 to 2010 in the TM region, which is consistent with the findings from other studies conducted in Brazil and Italy.^16,27 Both co-circulation of different genetic variants and a significant migration flow between countries contribute to the occurrence of reassortment events. Interestingly, a recent study in Nashville, TN, USA,⁴³ described the co-circulation of genetic variants that probably only underwent antigenic drifts without the occurrence of reassortments. This disparity in results may be due to a short period (2010-2011) examined in the above study.

In this study, we characterized, based on five genes, the genetic diversity of RVA strains representative of six seasons in the TM region. A great intragenotypic diversity was observed among the G2P[4] samples, with an apparent pattern of temporal segregation within the lineages or sublineages. The occurrence of intergenotypic reassortment events was evident, leading, in particular, to a combination of genotypes not previously reported (G3-P[4]-I2-E3-H3). Five TM samples contained genes that exhibited high nt and aa identities with strains of animal origin. These results highlight the importance of studying diversity of human and animal rotaviruses in order to better understand the mechanisms and the frequency of interspecies transmission. The data of this study cover a period pre- and post-introduction of rotavirus vaccination in all Brazilian territories, thereby serving as a basis for monitoring changes in the genetic constitution of rotaviruses. These changes can influence the viral fitness and the effectiveness of vaccines over time. The present study also contributes to the understanding of the genetic diversity and evolution of rotaviruses in a global context.

• Associate Editor: Maurício Lacerda Nogueira

Acknowledgements

We acknowledge Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Brazil for granting a PhD fellowship to ACB Dulgheroff; Fundação de Ensino e Pesquisa de Uberaba (FUNEPU), Federal University of Triângulo Mineiro, Brazil for the partial financial support and Programa de Desenvolvimento Tecnológico em Insumos para Saúde (PDTIS) Amazonian Oswaldo Cruz Foundation (FIOCRUZ), Brazil for permitting the use of their sequencing platform (Leonidas e Maria Deane Institute).

REFERENCES

• ¹
  Tate JE, Burton AH, Boschi-Pinto C, et al. 2008 estimate of worldwide rotavirus-associated mortality in children younger than 5 years before the introduction of universal rotavirus vaccination programmes: a systematic review and meta-analysis. Lancet Infect Dis 2012;12(2):136-141.
• ²
  Estes MK, Kapikian AZ. Rotaviruses. In: Knipe DM, Howley PM, eds. Fields Virology Vol 2, 5th ed. Philadelphia: Lippincott Williams & Williams; 2007:1917–1974.
• ³
  Dóró R, László B, Martella V, 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-461.
• ⁴
  Matthijnssens J, Ciarlet M, Rahman M, et al. Recommendations for the classification of group A rotaviruses using all 11 genomic RNA segments. Arch Virol 2008;153(8):1621-1629.
• ⁵
  Lappalainen S, Pastor AR, Tamminen K, et al. Immune responses elicited against rotavirus middle layer protein VP6 inhibit viral replication in vitro and in vivo. Hum Vaccin Immunother 2014;10(7):2039-2047.
• ⁶
  Ball JM, Tian P, Zeng CQ, Morris AP, Estes MK. Age-dependent diarrhea induced by a rotaviral nonstructural glycoprotein. Science 1996;272(5258):101-104.
• ⁷
  Domingues AL, Morais AT, Cruz RL, Moreira LP, Gouvêa VS. Rotavirus-associated infantile diarrhea in Uberaba, Minas Gerais, on the wake of the Brazilian vaccination program. J Clin Virol 2008;43(3):298-301.
• ⁸
  Dulgheroff AC, Figueiredo EF, Moreira LP, et al. Distribution of rotavirus genotypes after vaccine introduction in the Triângulo Mineiro region of Brazil: 4-year follow-up study. J Clin Virol 2012;55(1):67-71.
• ⁹
  Dulgheroff AC, Figueiredo EF, Gouvêa VS, Domingues AL. Changes in epidemiology of rotavirus in the Triângulo Mineiro region of Brazil: lack of two consecutive rotavirus seasons. Braz J Med Biol Res 2014;47(12):1091-1095.
• ¹⁰
  Boom R, Sol CJ, Salimans MM, Jansen CL, Wertheim-van Dillen PM, van der Noordaa J. Rapid and simple method for purification of nucleic acids. J Clin Microbiol 1990;28(3):495-503.
• ¹¹
  Gómara MI, Cubitt D, Desselberger U, Gray J. Amino acid substitution within the VP7 protein of G2 rotavirus strains associated with failure to serotype. J Clin Microbiol 2001;39(10):3796-3798.
• ¹²
  Simmonds MK, Armah G, Asmah R, et al. New oligonucleotide primers for P-typing of rotavirus strains: strategies for typing previously untypeable strains. J Clin Virol 2008;42(4):368-373.
• ¹³
  Cunliffe NA, Woods PA, Leite JP, et al. Sequence analysis of NSP4 gene of human rotavirus allows classification into two main genetic groups. J Med Virol 1997;53(1):41-50.
• ¹⁴
  Iturriza Gómara M, Wong C, Blome S, Desselberger U, Gray J. Molecular characterization of VP6 genes of human rotavirus isolates: correlation of genogroups with subgroups and evidence of independent segregation. J Virol 2002;76(13):6596-6601.
• ¹⁵
  Matthijnssens J, Rahman M, Martella V, et al. Full genomic analysis of human rotavirus strain B4106 and lapine rotavirus strain 30/96 provides evidence for interspecies transmission. J Virol 2006;80(8):3801-3810.
• ¹⁶
  Giammanco GM, Bonura F, Zeller M, et al. Evolution of DS-1-like human G2P[4] rotaviruses assessed by complete genome analyses. J Gen Virol 2014;95(Pt 1):91-109.
• ¹⁷
  Phan TG, Khamrin P, Quang TD, et al. Detection and genetic characterization of group A rotavirus strains circulating among children with acute gastroenteritis in Japan. J Virol 2007;81(9):4645-4653.
• ¹⁸
  Phan TG, Okitsu S, Maneekarn N, Ushijima H. Genetic heterogeneity, evolution and recombination in emerging G9 rotaviruses. Infect Genet Evol 2007;7(5):656-663.
• ¹⁹
  Rahman M, Matthijnssens J, Yang X, et al. Evolutionary history and global spread of the emerging g12 human rotaviruses. J Virol 2007;81(5):2382-2390.
• ²⁰
  Ianiro G, Delogu R, Bonomo P, Castiglia P, Ruggeri FM, Fiore L. Molecular characterization of human G8P[4] rotavirus strains in Italy: proposal of a more complete subclassification of the G8 genotype in three major lineages. Infect Genet Evol 2014;21:129-133.
• ²¹
  Ndze VN, Esona MD, Achidi EA, et al. Full genome characterization of human Rotavirus A strains isolated in Cameroon, 2010-2011: diverse combinations of the G and P genes and lack of reassortment of the backbone genes. Infect Genet Evol 2014;28:537-560.
• ²²
  da Silva MF, Gómez MM, Rose TL, et al. VP8*P[8] lineages of group A rotaviruses circulating over 20 years in Brazil: proposal of six different sub-lineages for P[8]-3 clade. Infect Genet Evol 2013;16:200-205.
• ²³
  Felsenstein J. Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol 1981;17(6):368-376.
• ²⁴
  Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004;32(5):1792-1797.
• ²⁵
  Matthijnssens J, Ciarlet M, McDonald SM, et al. Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG). Arch Virol 2011;156(8):1397-1413.
• ²⁶
  Jere KC, Mlera L, O’Neill HG, Peenze I, van Dijk AA. Whole genome sequence analyses of three African bovine rotaviruses reveal that they emerged through multiple reassortment events between rotaviruses from different mammalian species. Vet Microbiol 2012;159(1-2):245-250.
• ²⁷
  Gómez MM, de Mendonça MC, Volotão Ede M, et al. Rotavirus A genotype P[4]G2: genetic diversity and reassortment events among strains circulating in Brazil between 2005 and 2009. J Med Virol 2011;83(6):1093-1106.
• ²⁸
  Matthijnssens J, Van Ranst M. Genotype constellation and evolution of group A rotaviruses infecting humans. Curr Opin Virol 2012;2(4):426-433.
• ²⁹
  Doan YH, Nakagomi T, Cunliffe NA, Pandey BD, Sherchand JB, Nakagomi O. The occurrence of amino acid substitutions D96N and S242N in VP7 of emergent G2P[4] rotaviruses in Nepal in 2004-2005: a global and evolutionary perspective. Arch Virol 2011;156(11):1969-1978.
• ³⁰
  Matthijnssens J, Ciarlet M, Heiman E, et al. Full genome-based classification of rotaviruses reveals a common origin between human Wa-Like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains. J Virol 2008;82(7):3204-3219.
• ³¹
  Pietsch C, Petersen L, Patzer L, Liebert UG. Molecular characteristics of German G8P[4] rotavirus strain GER1H-09 suggest that a genotyping and subclassification update is required for G8. J Clin Microbiol 2009;47(11):3569-3576.
• ³²
  Luchs A, Cilli A, Morillo SG, Ribeiro CD, Carmona Rde C, Timenetsky Mdo C. Rotavirus genotypes and the indigenous children of Brazilian midwest in the vaccine era, 2008-2012: footprints of animal genome. J Med Virol 2015;87(11):1881-1889.
• ³³
  Silva MF, Rose TL, Gómez MM, et al. G1P[8] species A rotavirus over 27 years - pre- and post-vaccination eras - in Brazil: full genomic constellation analysis and no evidence for selection pressure by Rotarix™® vaccine. Infect Genet Evol 2015;30:206-218.
• ³⁴
  Matthijnssens J, Heylen E, Zeller M, Rahman M, Lemey P, Van Ranst M. Phylodynamic analyses of rotavirus genotypes G9 and G12 underscore their potential for swift global spread. Mol Biol Evol 2010;27(10):2431-2436.
• ³⁵
  Gómez MM, Resque HR, Volotão Ede M, et al. Distinct evolutionary origins of G12P[8] and G12P[9] group A rotavirus strains circulating in Brazil. Infect Genet Evol 2014;28:385-388.
• ³⁶
  Iturriza-Gómara M, Dallman T, Bányai K, et al. Rotavirus genotypes co-circulating in Europe between 2006 and 2009 as determined by EuroRotaNet, a pan-European collaborative strain surveillance network. Epidemiol Infect 2011;139(6):895-909.
• ³⁷
  Iturriza-Gómara M, Green J, Brown DW, Desselberger U, Gray JJ. Diversity within the VP4 gene of rotavirus P[8] strains: implications for reverse transcription-PCR genotyping. J Clin Microbiol 2000;38(2):898-901.
• ³⁸
  Rosa E, Silva ML, Pires De Carvalho I, Gouvea V. 1998-1999 rotavirus seasons in Juiz de Fora, Minas Gerais, Brazil: detection of na unusual G3P[4] epidemic strain. J Clin Microbiol 2002;40(8):2837-2842.
• ³⁹
  Tran TN, Nakagomi T, Nakagomi O. Evidence for genetic reassortment between human rotaviruses by full genome sequencing of G3P[4] and G2P[4] strains co-circulating in India. Trop Med Health 2013;41(1):13-20.
• ⁴⁰
  Maestri RP, Kaiano JH, Neri DL, et al. Phylogenetic analysis of probable non-human genes of group A rotaviruses isolated from children with acute gastroenteritis in Belém, Brazil. J Med Virol 2012;84(12):1993-2002.
• ⁴¹
  Luchs A, Timenetsky Mdo C. G8P[6] rotaviruses isolated from Amerindian children in Mato Grosso do Sul, Brazil, during 2009: close relationship of the G and P genes with those of bovine and bat strains. J Gen Virol 2014;95(Pt 3):627-641.
• ⁴²
  Matthijnssens J, De Grazia S, Piessens J, et al. Multiple reassortment and interspecies transmission events contribute to the diversity of feline, canine and feline/canine-like human group A rotavirus strains. Infect Genet Evol 2011;11(6):1396-1406.
• ⁴³
  Dennis AF, McDonald SM, Payne DC, et al. Molecular epidemiology of contemporary G2P[4] human rotaviruses cocirculating in a single U.S. community: footprints of a globally transitioning genotype. J Virol 2014;88(7):3789-3801.

Publication Dates

• Publication in this collection
  Jul-Sep 2016

History

• Received
  23 June 2015
• Accepted
  8 Dec 2015