Previous Article | Next Article ![]()
Journal of Clinical Microbiology, March 2007, p. 990-997, Vol. 45, No. 3
0095-1137/07/$08.00+0 doi:10.1128/JCM.01992-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Beatrice Karlsson,2,
Johan Nordgren,2
Margarita Paniagua,1
Alcides González,3
Juan Jose Amador,3
Felix Espinoza,1 and
Lennart Svensson2*
Department of Microbiology, National Autonomous University of Nicaragua, León, Nicaragua,1 Division of Molecular Virology, Medical Faculty, Linköping University, Linköping, Sweden,2 Ministry of Health, Managua, Nicaragua3
Received 26 September 2006/ Returned for modification 29 November 2006/ Accepted 4 January 2007
|
|
|---|
64,000 individuals and causing
56 deaths, predominantly in children under 5 years of age. Through a nationwide laboratory-based study, stool samples were collected and investigated for rotavirus. Of 108 stool samples examined, 72 (67%) were positive for rotavirus. While 69% (50/72) of the positive samples were found in children less than 2 years of age, 50% (6/12) of the adult samples were positive. A mutated G4P[8] strain was the most commonly recognized strain (85%), followed by mixed G strains (8%) and G9P[8] (7%) strains. Phylogenetic analysis of the VP7 gene revealed that the G4 strains belonged to the emerging lineage Ic and was distantly related to the ST3 and VA70 G4 strains. Secondary structure predictions of the VP7 G4 protein revealed an insert of an asparagine residue in position 76, which, combined with additional mutations, surprisingly modified two downstream ß-sheets at amino acid positions 80 to 85 and 115 to 119. The 2005 G4P[8] strain compared to a G4P[8] strain from 2002 had a substitution of an asparagine residue for threonine (Asn
Thr) at position 96 within antigenic region A, thus eliminating a potential glycosylation site. The mutated G4 virus was introduced in Nicaragua after 2002 and probably emerged from Brazil, Argentina, or Uruguay. |
|
|---|
It is estimated that all children will be infected at least once by the age of 5 years (5), which leads to protective immunity later in life. In adults rotavirus infections are usually asymptomatic, but there are reports of disease in elderly and immunocompromised patients (1, 4, 26, 28, 33).
Reports from different countries suggest an epidemiological shift of HRV strains over the years (20, 30). It seems that one specific variant or lineage within the serotypes might be responsible for outbreaks in well-defined geographic regions; for instance, the lineage Ic of the strain G4P[8] increased in prevalence in some cities of Argentina in 1998 and 1 year later in Paraguay (7, 9). The lineage Ic of the G4P[8] strain was also found to circulate in Italy in 1999 and 2000 (3), and during the rotavirus season in 1997 the emergence of a G4P[6] strain with a different RNA profile was observed in South Africa (34). Furthermore, in Nicaragua the uncommon strain G4P[6] was observed to circulate at very low frequency in 2001 and 2002 (12). The emergence of these novel strains and new lineages is due to different evolutionary mechanisms including gene rearrangements, accumulation of point mutations, and reassortment of genome segments (18).
During February and March 2005, one of the largest recorded outbreaks of severe acute gastroenteritis occurred in Nicaragua, affecting
64,000 individuals and causing at least 56 deaths, with children <5 years of age being most affected. We have carried out a molecular epidemiology study to investigate properties of the rotavirus strains associated with this unique and large epidemic outbreak. Our main observation was that the outbreak was associated with a mutated G4P[8] virus not previously detected in Nicaragua and most likely introduced from South America. The knowledge from this study may have implications for rotavirus prevention including vaccine introduction in Nicaragua.
|
|
|---|
The stool specimens were stored at 4°C until transported to the Department of Microbiology of the National Autonomous University of Nicaragua-Leon. A suspension of 10% (vol/vol) phosphate-buffered saline (pH 7.0) was prepared for HRV antigen detection and two aliquots of stool-phosphate-buffered saline suspensions were stored at 20°C for further analysis.
The laboratories of the health care facilities from the health system participated in the survey. Samples were geographically representative from all over the country and collected as follows: Pacific region, 68% (including Leon, 40% [43/108]; Managua, 21% [23/108]; Chinandega, Granada, Masaya, and Jinotepe, 7% [8/108]); northern and central regions, 29% (including Matagalpa-Jinotega, 11% [12/108]; Esteli, 10% [11/108]; Madriz and Nueva Segovia, 5% [5/108]; and Chontales, 3% [3/108]); and the Caribbean coast, 3% (4/108). While age was not reported in 6% (7/108) of the cases, 40% (43/108) of the samples were identified as coming from children <12 months old, 82% (89/108) from children <5 years old, and 11% (12/108) from young people and adults (Table 1).
|
View this table: [in a new window] |
TABLE 1. Distribution of G and P types of HRV strains isolated during an epidemic outbreak of acute gastroenteritis in Nicaragua 2005
|
RNA extraction. Viral double-stranded RNA (dsRNA) was extracted from stool suspensions following the manufacturer's instructions using a QIAamp Viral RNA Mini Kit (QIAGEN, Hilden, Germany). A total of 60 µl of viral dsRNA was collected and stored at 20°C until RNA gel electrophoresis or reverse transcription was carried out.
Polyacrylamide gel electrophoresis. RNA gel electrophoresis and negative silver staining were performed as previously described (35).
G and P multiplex genotyping. G and P genotyping was performed as described previously (21). This method is a modification of the original methods described by Gouvea et al. for G typing and Gentsch et al. for P typing (15, 16) and includes random primers [pd(N)6] to produce the cDNA and modified G and P primers.
RT. Reverse transcription (RT) was essentially carried out as described previously (19). Briefly, 28 µl of dsRNA was mixed with 50 pmol of random hexadeoxynucleotides [pd(N)6], denatured at 97°C for 5 min, and quickly chilled on ice for 2 min, followed by addition of one RT-PCR bead (Amersham Biosciences, United Kingdom) and RNase-free water to a final volume of 50 µl. The RT reaction was carried out for 30 min at 42°C to produce the cDNA used for genotyping and sequence analysis.
PCR. A PCR mix was prepared using 5 µl of 10x Native plus PFU buffer (Stratagene, La Jolla, CA), 1 µl of 10 mM deoxynucleoside triphosphate (dNTP) mix (Applied Biosystems, Warrington, United Kingdom), 4 pmol of each consensus primer (VP7-F and VP7-R) (17), 2.5 U of Native DNA polymerase (Stratagene, La Jolla, CA), and RNase-free water to a final volume of 50 µl. The PCR was performed at 94°C for 2 min, followed by 35 cycles of 94°C for 1 min, 50°C for 1 min, and 72°C for 1 min, with a final extension of 72°C for 7 min.
The amplicons were analyzed by gel electrophoresis using 2% agarose gel and ethidium bromide staining, and a consensus region of 881 bp of the VP7 encoding gene was obtained. The PCR products were purified with spin column purification (QIAprep Spin Miniprep Kit; QIAGEN, Hilden, Germany) and the amount of DNA was determined by a NanoDrop ND-1000 UV-visible light spectrophotometer (Saveen Werner AB, Malmö, Sweden).
Cloning. Purified PCR products were cloned into the pPCR-Script Amp SK(+) cloning vector, followed by transformation into XL10-Gold Kan ultracompetent cells according to the manufacturer's instructions (Stratagene, La Jolla, CA). The transformed bacteria were examined for recombinant plasmids by blue and white screening and PCR. The PCR mixture contained 2.5 µl of 10x PCR buffer (Invitrogen, Carlsbad, CA), 1 µl of 50 mM MgCl2 (Invitrogen, Carlsbad, CA), 5 pmol of each standard sequencing M13 forward and reverse primer, 2 µl of 2.5 mM dNTP (Invitrogen, Carlsbad, CA), 1 U of Taq DNA polymerase (Invitrogen, Carlsbad, CA), and RNase-free water to a final volume of 25 µl. The PCR was performed at 94° for 5 min, followed by 35 cycles of 94° for 30 s, 50°C for 30 s, and 72°C for 1 min, with a final extension of 72°C for 10 min.
Nucleotide sequencing. One independent colony from each cloned sample containing the desired insert was collected for sequencing. However, before sequencing, the PCR products were cleaned from excess dNTP and primers by using ExoSap-IT (GE Healthcare, Chalfont St. Giles, United Kingdom), in a reaction mixture containing 2.5 µl of PCR product and 1 µl of ExoSap-IT enzyme. This reaction was performed in 37°C for 15 min, followed by inactivation of the enzyme at 80°C for 15 min. Nucleotide sequences were obtained according to the manufacturer's instructions using a DYEnamic Dye Terminator Kit (GE Healthcare, Chalfont St. Giles, United Kingdom), M13 forward and reverse primers, and a Mega BACE 500 automated sequencer (GE Healthcare, Chalfont St. Giles, United Kingdom). For accuracy, clones with inserts were sequenced four times in both the forward and reverse directions. Complete sequences were obtained by assembling overlapping contigs with DNASTAR (Madison, WI).
Sequence analysis. Multiple sequence alignment of VP7 G4 proteins and nucleotide sequencing were performed, using the ClustalW algorithm, version 1.8, with default parameters on the European Bioinformatics Institute server. A phylogenetic analysis of the VP7 nucleotide alignment of the G4 and G9 strains was performed using the MEGA 3.1 software package. A tree was constructed using the neighbor-joining and Kimura two-parameter methods. The statistical significance of the relationships obtained was estimated by bootstrap resampling analysis (1,000 replications).
Secondary structure prediction. Secondary structure prediction of the VP7 G4 protein was performed using the PSIPRED protein structure prediction server (University College London) (8, 29).
|
|
|---|
High prevalence of rotavirus in children and adults. Rotavirus antigen screening revealed that 72/108 (67%) samples examined from children and adults suffering from diarrhea were positive. Furthermore 69% (50/72) of the rotavirus-positive samples were obtained from children less than 2 years of age, and 14% (10/72) were from children between 2 and 5 years of age. Eight percent of the positive samples (6/72) were derived from young people and adults (12 to 72 years of age; average, 29 years of age), but more important, 50% (6/12) of all the young and adult patients were rotavirus positive (Table 1). Four of the six were from the same city, and three of those from the same neighborhood.
Nationwide rotavirus outbreak predominantly associated with G4P[8]. To obtain information about the G and P type distribution, all 72 HRV-positive isolates were processed for G and P genotyping. G4 strains were found in 85% (61/72) of the samples, G9 in 7% (5/72), and mixed G infections in 8% (6/72) of the samples, including G4-G1 (two samples), G4-G9 (two samples), G3-G9 (one sample) and G4-G3-G9 (one sample). To eliminate any possibility of cross-priming occurring in these mixed infections, all the mixed infections were reanalyzed by single-locus PCR with genotype-specific primers to confirm the presence of the constituent G types. Only one sample was not consistent with the multiplex RT-PCR and the single-locus PCR analysis. That sample was the G4-G8 mixed infection in the first analysis, but no result was observed in the single-locus PCR analysis with the G8 primer. However, a clear G4 genotype was observed with the G4 primer (data not shown). Thus, five mixed infections occurred and all in children less than 5 years of age.
The P typing revealed that all strains belonged to the P[8] genotype (data not shown). The most common G and P combinations were G4P[8], representing 85% (61/72), and G9P[8], with 7% (5/72). All mixed infections were P[8], and a G9P[8] isolate came from a child less than 1 year of age living in Esteli. The G4P[8] strains were detected in all age groups and all over the country, including the Pacific, central, and Caribbean regions (Table 2).
|
View this table: [in a new window] |
TABLE 2. Description of isolates used for VP7 gene nucleotide sequencing
|
Nucleotide sequence analysis. Seven G4 strains, including a sample isolated in 2002, representing different age groups and geographic locations, were selected for VP7 sequence analysis. In addition a G9P[8] strain from Esteli was sequenced (Table 2).
Bioinformatic analysis revealed that the G9P[8] strain (Nic03/05 [the year of isolation is indicated after the slash]) circulating in Esteli during the 2005 outbreak shared >99% homology at the nucleotide level with strains belonging to common G9 strains and particularly with strain IN826 (USAIN826/00) circulating in the United States in the years 2000 to 2001 (25) (Fig. 1). The G4P[6] virus (NicBH63/02) that circulated with very low (6%) frequency in 2002 in Nicaragua (12, 13) shared 94% homology at the nucleotide level with the G4P[8] strains circulating during the outbreak. Nucleotide sequence alignment of the VP7 gene from the G4P[8] strains (Table 2 and Fig. 2) indicated that they all belonged to the same cluster and shared >99% homology at the nucleotide level with strains isolated in Argentina, Uruguay, and Brazil between 1996 and 2004.
![]() View larger version (18K): [in a new window] |
FIG. 1. Phylogenetic analysis of VP7 nucleotide sequences of G4 and G9 strains. The tree was constructed using the Kimura two-parameter and neighbor-joining methods using MEGA 3.1 software. Bootstrap values are shown at the branch nodes (values of <50% are not shown). The Nicaraguan strains are marked with a bold triangle and the lineages (I) and genotypes (G) are indicated at the right. The number of substitutions per site is indicated by the scale bar. Abbreviations for locations: Aus, Australia; Bel, Belgium; Bra, Brazil; Fin, Finland; Mvd, Montevideo (Uruguay); Nic, Nicaragua; Jap, Japan; Py, Paraguay; SA, South Africa; Swe, Sweden; Tha, Thailand.
|
![]() View larger version (62K): [in a new window] |
FIG. 2. Unique asparagine insert in G4P[8] virus from 2005 in Nicaragua. Amino acid alignment of the rotavirus VP7 gene using ClustalW, version 1.8, with default parameters. The alignment reveals an asparagine insert at position 76 in the 2005 Nicaraguan isolates, which also is present in isolates from Argentina (1996 to 1998), Uruguay (1999), and Brazil (2000 to 2004). *, conserved residue; :, position with conserved substitutions; ., position with semiconserved substitutions.
|
|
View this table: [in a new window] |
TABLE 3. Amino acid evolution in rotavirus G4 VP7 in Nicaragua between 2002 and 2005
|
![]() View larger version (29K): [in a new window] |
FIG. 3. Amino acid insertions and substitutions affect downstream ß-sheets in VP7. Schematic representations of secondary structure predictions, using the PSIPRED protein structure prediction server (8). Both panels depict aa 41 to 160. The residue number is shifted 1 aa in the G4 VP7 sequence (Nic118/05) from position 76 compared to the 2002 strain (NicBH63/02). (Boxes a) Beginning at amino acid 80, the 6-amino-acid ß-sheet in the 2002 isolate has mutated into a 7-amino-acid ß-sheet in the 2005 isolate. (Boxes b) Located in the vicinity of amino acid 120, the 5-amino-acid ß-sheet in the 2002 isolate has evolved into a 6-amino-acid ß-sheet in the 2005 isolate. Conf, confirmed; Pred, predicted.
|
Ile and Asn
Thr) and B (Arg
Lys and Ala
Thr) (Tables 3 and 4). Most interesting was the observation that the amino acid sequences in all three antigenic regions of G4P[8] strains isolated in Argentina, Uruguay, and Brazil between 1998 and 2004 were identical to the G4P[8] virus isolated in Nicaragua in 2005 but not in 2002. |
View this table: [in a new window] |
TABLE 4. The VP7 amino acid sequences in the antigenic regions of the G4P[8] outbreak virus are identical to those of the G4 viruses isolated in Argentina (1998), Uruguay (1999), and Brazil (2000 to 2004)
|
![]() View larger version (107K): [in a new window] |
FIG. 4. Amino acid alignment using ClustalW, version 1.8, with default parameters on the European Bioinformatics Institute server of the rotavirus VP7 gene showing the number of amino acid substitutions in the 2005 Nicaraguan isolates compared to the isolate from 2002. The substitutions occurring at positions A to I are all conserved or semiconserved. The three important antigenic sites (VR-5, aa 86 to 101; VR-7, aa 142 to 152; and VR-8, aa 211 to 223) are underlined. The location of the asparagine insert at position 76 is indicated with a circle. *, conserved residue; :, position with conserved substitutions; ., position with semiconserved substitutions.
|
The G4P[8] strain of 2002 (NicBH63/02) was identical to strain ST3 in antigenic region C, while having substitutions in regions A and B. Compared to the VA75 strain, the G4P[8] 2002 strain as well as the G4P[8] strain isolated in 2005 had substitutions in all the compared regions. However, the VA75 strain and the G4P[8] strain isolated in 2002 have two identical substitutions compared to the 2005 G4P[8] strain, where there is one substitution of an asparagine residue for a threonine (Asn
Thr) at position 96 within antigenic region A, thus eliminating a potential glycosylation site (Table 4).
|
|
|---|
Sequence analysis of the VP7 gene from G4P[8] strains circulating in Argentina in 1998 and 1 year later in Paraguay and Uruguay demonstrated that an emerging lineage of this serotype was circulating in those countries (6, 7, 9). Interestingly, the G4P[8] virus identified in this outbreak shares >99% homology at the nucleotide level with strains isolated in Argentina, Uruguay, and Brazil from 1998 to 2004 and belongs to the G4 lineage Ic (6, 7, 36). This may indicate that introduction of the mutated and highly virulent G4P[8] variant in Nicaragua probably originated from the South American region. In support of this suggestion is the finding that the G4P[8] outbreak strain is less related to strains Hochi and Odelia isolated in Japan in 1980 and 1984, respectively, and to the South African G4 strains isolated in 1986 and the Nicaraguan G4 strains isolated in 2002 (13, 24, 31). The outbreak of G4P[8] virus in 2005 may thus have been a result of infection of a population without preexisting immunity. The fact that 50% of the examined adults were symptomatically infected with G4P[8] may indicate that the G4P[8] virus has been absent for some time in Nicaragua or that the virus more recently has mutated. The Nicaraguan G4 VP7 strains all had a unique asparagine residue insertion at position 76, an insertion also found in strains isolated in Argentina, Uruguay, and Brazil during 1998 to 2004. This finding also supports the suggestion of an emergence of G4P[8] virus from South America to a highly susceptible population in Nicaragua. Another, albeit less likely, explanation is that the G4P[8] outbreak strain may have originated from reassortment events between G4P[6] and G1P[8] both circulating at low (5%) and high (45%) frequencies, respectively, during 2001 and 2002 (13).
The asparagine insertion at position 76 in the VP7 gene is close to the glycosylation motif Asn-X-Thr at residues 69 to 71, a site conserved in most human G4 strains, and it cannot be ruled out that the insertion influences glycosylation and, hence, alters antigenic properties. Residue 76 is located in a hydrophilic region, and an insertion of an asparagine in this area may cause an increase in hydrophilicity in the region (6). Most surprising was the observation that the asparagine insertion together with other mutations resulted in two minor structural modifications by altering two downstream ß-sheets in amino acid positions 80 to 85 and 115 to 119 (Fig. 3). The observation that the mutations affected downstream ß-sheets is interesting and warrants further investigation.
The amino acid substitutions that occurred between 2002 and 2005 in the Nicaraguan strains influenced the amino acid composition of important antigenic sites (10, 11). As illustrated in Table 3, the substitutions are located at position 94 (Val
Ile) in antigenic region A and positions 144 (Arg
Lys) and 146 (Ala
Thr) in antigenic region B. Furthermore, we found that G4P[8] strains isolated in Brazil, Uruguay, and Argentina between 1998 and 2002 had amino acid sequences in antigenic regions A to C that are identical to those of the G4P[8] strain isolated in this outbreak but different from the Nicaraguan strain isolated in 2002. These findings suggest that the reemergence of G4 strains in Nicaragua could be due to immune evasion as a consequence of altered antigenicity conferred not only by the asparagine residue insertion in VP7 at position 76 but also by the substitution of the asparagine (Asn
Thr) residue at position 96 in antigenic region A, which eliminated an N-linked glycosylation site and thus possibly altered antigenicity (11). An older population, presumably with preexisting immunity but with a high incidence of HRV infection due to immune evasion, has also been suggested by Iturriza-Gomara et al. in a report describing a G2 strain causing gastroenteritis between 1995 and 1998 in the United Kingdom (20, 35).
It has previously been observed that the antigenic regions in the VP7 protein, even though distant in the linear molecule, interact closely together in the folded form of the VP7 molecule and thus influence, for instance, antibody binding and immune responses (11). Indeed, a single amino acid substitution in VP7 and predominantly at antigenic regions A, B, and C has been observed to alter antigenicity (22). Thus, it cannot be ruled out that the identified mutations result in a mutated virus displaying a different molecular makeup, conferring an increased virulence enabling this particular virus strain to escape neutralization by G4 antisera. Indeed, in support for this hypothesis, it has been suggested by Iturriza-Gómara et al. that the epidemic reemergence of G2 rotavirus strains in Taiwan was due to immune evasion as a consequence of altered antigenicity conferred by an amino acid substitution at position 96 in antigenic region A of the VP7 gene (17, 37). Interestingly, an identical substitution seen in the Nicaraguan strain (Asn
Thr) at position 96 was also observed in G4P[6] symptomatic strains isolated from neonates in South Africa (31). This suggests, but does not prove, that substitution of an asparagine to a threonine at position 96 may lead to increased virulence as seen in G4P[8] strains of this outbreak. Another observation in favor of increased virulence through immune evasion is the fact that 50% of the adults with diarrhea were rotavirus positive. However, it cannot be ruled out that the lack of immunity might have been associated with the P type.
Apart from El Salvador, the G9P[8] virus had not previously been detected in Central America (2, 12). In this study it represented 7% of the HRV-positive cases in children attending the health care centers of Managua and Esteli. An interesting observation was that no apparent difference in RNA profiles from samples collected in different geographic settings or belonging to different G types (G4P[8] or G9P[8]) was found. There is a possibility that the G9P[8] strain examined might have originated from the United States because of its high nucleotide homology (>99%) with strain IN826 (USAIN826/00) isolated in Indianapolis (Indiana) in 2000 to 2001. In support of this explanation is also the fact that both strains share a long RNA profile (25). Considering that serotype G9 is emerging and may present a progressive increase in occurrence, as happened in Australia from 1997 until 2001 (23), it is reasonable to suggest that this G type strain will continue to increase in Nicaragua and elsewhere in the Central American region.
Published ahead of print on 17 January 2007. ![]()
F.B. and B.K. contributed equally to this work. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»