Previous Article | Next Article ![]()
Journal of Clinical Microbiology, February 2009, p. 421-429, Vol. 47, No. 2
0095-1137/09/$08.00+0 doi:10.1128/JCM.01852-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Laboratory of Infectious Diseases and Biological Agents, Faculty of Pharmacy, TU-5000 Monastir, Tunisia,1 National Reference Center for Enteric Viruses, Laboratory of Virology, University Hospital of Dijon, F-21070 Dijon, France,2 Laboratory of Immunology, University Hospital Fattouma Bourguiba, TU-5000 Monastir, Tunisia,3 Pediatric Department, University Hospital Fattouma Bourguiba, TU-5000 Monastir, Tunisia4
Received 25 September 2008/ Returned for modification 15 November 2008/ Accepted 12 December 2008
|
|
|---|
|
|
|---|
In many countries, numerous works have studied the incidence and the molecular epidemiology of NoVs causing gastroenteritis. Indeed, each year, NoVs may cause more than 1.1 million hospitalizations and 200,000 deaths in children <5 years of age in developing countries (28). However, in Tunisia, as in many other African countries, while diarrhea is a major cause of illness among children, little is known about the molecular epidemiology of NoVs in gastroenteritis. A first study conducted by our laboratory between January 2003 and June 2005 showed NoVs as the second most common agent responsible for sporadic cases of viral gastroenteritis in Tunisian children but equal to rotaviruses in terms of necessity of hospitalization and severity of the clinical symptoms (32). This present study completes the previous survey by focusing on the molecular epidemiology of NoV strains circulating in Tunisia. Indeed, we screened 788 stool specimens collected from Tunisian children under 12 years of age who had acute diarrhea between January 2003 and April 2007. Using reverse transcription-PCR (RT-PCR) with GI- and GII-specific primer sets and sequencing, we were able to phylogenetically analyze the strains. The genetic characterization of the NoVs detected in our study showed a diversity of strains circulating in Tunisia but with a high predominance of the GGII.4 genotype. Surprisingly, our results show the presence of the GGII.4 variant Hunter from the beginning of our study in January 2003, whereas its first isolation was reported in the literature to be in February 2004, which is 1 year later.
|
|
|---|
Specimens. Fecal specimens were collected within 24 h of admission. Four hundred eight (51.8%) specimens were collected from hospitalized children (inpatients) throughout the study, and 380 (48.2%) were collected from outpatient children presenting in dispensaries from January 2003 to May 2004. They were all negative for bacteria and parasites. The specimens were frozen, sent to the laboratory, and then stored at –40°C until being tested for NoVs.
RNA extraction. Viral RNA was extracted from 10% stool suspensions in phosphate-buffered saline (pH 7.5) with a QIAamp viral RNA kit (Qiagen, Hilden, Germany), according to the manufacturer's instructions, and was stored at –40°C.
RT-PCR and sequencing. NoVs were detected by RT-PCR using several sets of primers in separate reactions. The primer set JV12/JV13 (40) was used to amplify a fragment of the RNA polymerase gene. The primers sets G1SKF/G1SKR and G2SKF/G2SKR (17) were used to detect a fragment of the capsid gene of NoV genogroups I and II, respectively. RT-PCR was performed with the One Step RT-PCR kit (Qiagen) according to the manufacturer's instructions and to the cycles of amplification given by the authors of each primer set.
All the amplified cDNA samples were systematically purified from the gel by using the QIAquick gel extraction kit (Qiagen). Sequencing of the PCR products was performed with the same primers as those for amplification by using the ABI Prism BigDye Terminator cycle sequencing ready reaction kit (Applera Corporation, Foster City, CA) and an ABI 3100 automated sequencer (PE Biosystems). The entire sequences of ORF2 were amplified for seven samples (five from 2003, one from 2004, and one from 2007) using primer pair FW1/RT5 and the Titan one-tube RT-PCR kit (Roche Applied Biosystems) (A. H. Kamel, M. A. Ali, P. Pothier, and G. Belliot, submitted for publication).
The ORF1/ORF2 junction of the viral genome was sequenced for all positive samples (n = 14) whose genotype in the capsid gene differs from the genotype in the polymerase gene. The PCR products obtained with the primer set JV12/G2SKR (1,112 bp) were cloned into the pGEM-T Easy vector system (Promega Corporation, Madison, WI), according to the manufacturer's instructions, and then sequenced using the same primers.
Characterization and phylogenetic analysis. The nucleotide sequences of the amplicons were compared to corresponding sequences of NoV strains available in the GenBank database using Fasta program version 3, available from the European Bioinformatics Institute (http://www.ebi.ac.uk). For the RNA polymerase, analyses were also performed with reference strains available in the database of the European Food-Borne Viruses Network (https://hypocrates.rivm.nl/bnwww/Divine-Event/index.html). Sequence alignments were carried out by using Clustal W software. RNA polymerase nucleotide sequence identities within a genotype range from 87 to 99% and 91 to 100% for GGI and GGII NoVs, respectively. Capsid nucleotide sequence identities range from 88 to 99% within a genotype for both genogroups (41).
For the phylogenetic analyses, sequence alignment and clustering were performed by the UPGMA (unweighted-pair group method using average linkages) using Bionumerics software (Applied Maths).
Nucleotide sequence accession numbers. The sequences obtained in this study have been submitted to GenBank under the following accession numbers (isolate numbers and dates are given in parentheses): EU650205 (Monastir 18/2003), EU650206 (Monastir 294/2003), EU650207 (Monastir 3/2003), EU650208 (Monastir 273/2003), EU650209 (Monastir 430/2003), EU650210 (Monastir 17285/2004), EU650211 (Monastir 127/2003), EU650212 (Monastir 474/2003), EU650213 (Monastir 1795/2004), EU650214 (Monastir 578/2003), EU650215 (Monastir 17291/2006), EU650216 (Monastir 519/2003), EU650217 (Monastir 862/2003), EU650218 (Monastir 39/2004), EU650219 (Monastir 249/2003), EU650220 (Monastir 310/2003), EU650221 (Monastir 375/2003), EU650223 (Monastir 8655/2007), EU650224 (Monastir 389/2003), EU650225 (Monastir 17736/2005), EU650226 (Monastir 13562/2004), and EU650227 (Monastir 4456/2005). For the entire ORF2 sequences, accession numbers are as follows: EU916955 (Monastir 3968/2004), EU916956 (Monastir 8655/2007), EU916957 (Monastir 127/2003), EU916958 (Monastir 113/2003), EU916959 (Monastir 529/2003), EU916960 (Monastir 493/2003), and EU916961 (Monastir 715/2003).
|
|
|---|
The distribution of the sampling during the years of our study is indicated in Table 1: 76.6% (n = 604) of the samples were collected during 2003 and 2004, and 83.6% of NoV-positive samples were detected in these years. However, the ratio between the number of NoV-positive samples and the number of samples tested decreased over time.
|
View this table: [in a new window] |
TABLE 1. Annual distribution of NoV-positive samples detected in Monastir, Tunisia, between January 2003 and April 2007
|
![]() View larger version (11K): [in a new window] |
FIG. 1. Monthly distribution of NoV genotypes and numbers of samples tested in Tunisian children between January 2003 and April 2007. J through D, January through December, respectively.
|
Based on the nucleotide sequences, representative phylogenetic trees of the polymerase and capsid regions were constructed by the UPGMA method using Bionumerics software (Fig. 2 and 3). Phylogenetic analysis of both regions showed that the GGII.4 strains identified in this study from 2003 to 2007 were all closely related to the GGII.4 2004 variant Hunter (GenBank accession numbers DQ078801 and DQ078794), first described in Australia in February 2004 (9). The sequence analysis based on 285 bp from the polymerase region, which is the most conserved region of the genome and the most widely analyzed in the literature, showed that our strains differ from Hunter by 1 bp. The analysis of 302-bp sequences of the capsid gene showed 100% identity with Hunter.
![]() View larger version (20K): [in a new window] |
FIG. 2. Phylogenetic analysis based on the partial nucleotide sequences (285 bp) of the RNA-dependent RNA polymerase coding gene of the NoV strains associated with pediatric gastroenteritis in Monastir, Tunisia, between January 2003 and April 2007. The tree was constructed using UPGMA clustering. GenBank accession numbers for reference strains in this figure are as follows: Hunter virus, DQ078794; Southampton virus, L07418; Chiba virus, AB042808; Amsterdam virus, AB025907; GGIIb virus, AJ487794; Arg 320 virus, AF190817; Farmington Hills virus, AY502023. For more clarity, the sequences that present 100% nucleotide identity are not all represented on the tree. Only some of them, representative of each genotype, are indicated. In the same way, only some representative sequences were submitted to GenBank. Their GenBank accession numbers are indicated in Materials and Methods, and they are marked with asterisks in the tree. The percent bootstrap values in which the major groupings were observed among 100 replicates are indicated in each branch.
|
![]() View larger version (20K): [in a new window] |
FIG. 3. Phylogenetic analysis based on the partial nucleotide sequences (302 bp) of the capsid coding gene of NoV isolates detected in Tunisia. Reference strains of NoV were selected from the GenBank database under the following accession numbers: Farmington Hills virus, AY502023; Hunter virus, DQ078794; Amsterdam virus, AF195848; Hawaii virus, U07611; Melksham virus, X81879; Arg 320 virus, AF190817; M7/1999/US virus, AY130761. The GenBank accession numbers of strains of this study are indicated in the text, and they are marked with asterisks in the tree. The numbers in the branches indicate the bootstrap values.
|
![]() View larger version (34K): [in a new window] |
FIG. 4. Phylogenetic tree constructed from complete amino acid sequences of ORF2 of the seven Tunisian isolates and the GGII.4 representative strains from GenBank. For the samples in this study, accession numbers are as follows: EU916955 (Monastir-3968/2004), EU916956 (Monastir-8655/2007), EU916957 (Monastir-127/2003), EU916958 (Monastir-113/2003), EU916959 (Monastir-529/2003), EU916960 (Monastir-493/2003), and EU916961 (Monastir-715/2003). The numbers on each branch indicate the bootstrap values for the clusters supported by that branch.
|
The 10 GGI.2 strains showed 94.7% nucleotide identity in the polymerase region with the reference strain Southampton/1991/UK (GenBank accession number L07418), and the strain GGI.4 presented 90.8% nucleotide identity in the polymerase gene with Chiba/1987/JP (accession number AB042808). The capsid gene of these strains failed to amplify by our method. Among the four GGII.8 strains, two were amplified in the polymerase gene and displayed 99 to 100% nucleotide identity with strains referenced in the database of the European Food-Borne Viruses Network. The two other GGII.8 strains were amplified in the capsid gene and showed 96% nucleotide identity with Amsterdam/1998/NL (accession number AF195848). The four strains of GGII.14 displayed 95.7% nucleotide identity in the capsid gene with the strain designated M7/1999/US (accession number AY130761). The 24 GGIIb/GGII.3 strains found in our study were all sequenced in the polymerase and capsid genes. They displayed 97.5% nucleotide identity in the polymerase gene with the GGIIb strain described by Buesa et al. (6) (accession number AJ487794). Their capsid region showed 97.4% nucleotide identity with Mexico/1995 (accession number U22498). The 11 GGIIb/GGII.2 strains were divided into two groups that differed by five nucleotides according to the alignment of 267-bp sequences of the capsid gene; these two groups of six and five strains displayed 97.4% and 98.5% nucleotide identity with Melksham/1994/UK (accession number X81879), respectively. Their polymerase gene showed 97.9% nucleotide identity with the GGIIb strain (accession number AJ487794). All the GGII.1 strains shared 95.3% nucleotide identity in the capsid gene to the Hawaii strain (accession number U07611), except for one that differed by six nucleotides from the others but presented also 95.3% nucleotide identity with the Hawaii strain. The polymerase gene of these strains could not be amplified by our method.
|
|
|---|
Our results have shown the presence of eight distinct genotypes of NoVs. This is consistent with the genetic diversity found in other studies (1, 3, 6, 14, 24, 31) and indicates that gastroenteritis can be due to genetically different NoVs in the same area. However, the prevalence of the different genotypes changes over time, with a period of cocirculation of numerous strains, followed by periods of circulation of almost one unique genotype. Indeed, we observed a cocirculation of seven distinct strains until May 2004 (including GGIIb/GGII.3 and GGII.4), and then strains classified as GGIIb/GGII.3 were almost the only ones detected in the second half of 2004. Again in 2005 the pattern of circulation changed, with the disappearance of these strains and their almost-exclusive replacement by the GGII.4 genotype. These changes in the predominance of a particular genotype have already been observed in several countries; the GGII.3 strain, predominant in Japan in the winter season 2001-2002 (26), and the GGII.1 strain, predominant in the Italian region of Parma and in Hungary in 2000 (24, 31), were all replaced the next season by the emerging strain GGII.4. In our study the end of cocirculation of diverse strains in May 2004 coincides with the end of collecting stools from children presenting in dispensaries in addition to hospitalized children. However, no correlation was found between the two events (data not shown), challenging the hypothesis that some genotypes could be responsible for more severe symptoms and therefore be more often detected in the hospitalized population than in the community.
The natural emergent recombinant GGIIb strains were detected for the first time in France in August 2000 (2) and then spread throughout France and Europe in the following seasons, with outbreaks and sporadic cases reported in Spain, Sweden, and Hungary (6, 20, 31). They were also reported in other parts of the world, such as in Japan in the 2003-2004 season (29). Our study showed that the recombinant GGIIb strains were the second most frequently detected NoVs in Tunisian children, as previously found in Hungary between 2001 and 2003 (31). However, it is interesting that these strains were observed in our study only in 2003 and 2004 and have not diffused since. The same pattern has been observed in France, with a high prevalence until 2004 (5) and then only episodic gastroenteritis cases due to the GGIIb strains occurring (personal data).
Except between June and December 2004 when GGIIb/GGII.3 was the predominant circulating genotype, GGII.4 was the most prevalent cluster isolated during the whole study. This genotype is considered the predominant genotype worldwide, responsible for the majority of outbreaks of gastroenteritis since 1995 (9, 19, 22, 25, 35), and is also shown in recent molecular epidemiological studies as the most prevalent genotype in sporadic community cases in recent years (1, 3, 10, 45). The emergence and global spread of this virus coincided with a worldwide increase in the number of gastroenteritis outbreaks, suggesting that this strain could be more virulent or better adapted for survival in the environment or may be more contagious than other NoVs (43). It is interesting that this overall increase in the number of outbreaks is punctuated with epidemic seasons due to GGII.4, followed by postepidemic years with non-GGII.4 strains being more commonly found in outbreaks (34). This pattern of evolution has been termed epochal evolution. Several studies focused on the GGII.4 genotype and have shown that several epidemic variants can be distinguished within the genotype that succeed one another over time: one variant emerges, replaces the previously circulating variant to become the predominant strain for a few years, and is ultimately replaced with another new variant (7, 13, 21, 33, 34). This phenomenon has been explained by an accumulation of point mutations, essentially in the protruding P2 region of the capsid protein which conduces to genetic drift and an immune evasion. Five lineages have been identified since 1995 that are responsible for the four epidemic peaks observed worldwide in 1995-1996, 2002, 2004, and 2006. The 95/96US cluster observed in the United Kingdom, the United States, Brazil, Canada, Australia, The Netherlands, China, and Germany (25, 43) has been displaced by the variant Farmington Hills, isolated in 2002 in Europe (23) and in the United States (44). Then a new variant termed Hunter (9) became the most common circulating strain worldwide in 2004, displacing the Farmington Hills variant (18). Lately two new variants, 2006a and 2006b, have been described as the global predominant strains (7, 13, 15, 19, 34, 35, 36, 37, 39). Our results show the presence of the 2004 variant Hunter from the very start of our study, in January 2003, as confirmed by the sequencing of the complete ORF2 of two samples collected in that month. This is contrary to the literature that dates the first isolation of this variant to February 2004 in New South Wales (9), 1 year later. The presence of a specific strain in sporadic infantile cases before its global spread has been already pointed out by Medici et al. (24), who proposed that children with sporadic gastroenteritis may act as a reservoir for emerging epidemic NoV strains. In fact our results do not correspond to the general schema of succession of epidemic variants over time or to the worldwide succession of epidemic seasons and global increase of the number of gastroenteritis cases explained by this succession of epidemic variants. Indeed, the systematic sequencing of the polymerase and capsid fragments of all our GGII.4 samples, as well as the sequencing of the complete ORF2 of seven samples, shows that the same variant is found throughout our study, even in the three samples collected in 2006 and 2007, after the global spread of the 2006 variants.
Even if we observed peaks of infection in winter and in summer, we globally observed a decrease of the number of NoV-positive samples, essentially since January 2005. This decrease can be explained in part by the decrease of the number of samplings after May 2004. Nevertheless, a real decrease exists of the ratio between the number of samples tested and the number of samples positive for NoV over time. The reason for this phenomenon is unclear, but the small number of positive samples does not allow us draw conclusions about the circulating strains, especially GGII.4 variants, in 2006 and 2007 in the child population of Monastir and needs more investigation. In addition, it would be interesting to analyze samples from 2002 to see what variant was circulating at that time.
In conclusion, our study shows the predominance of the GGII.4 genotype in sporadic community cases of gastroenteritis in Tunisian children, interestingly with the appearance of the Hunter variant as early as 2003, 1 year before its worldwide emergence. Even if this delay of 1 year between the two events is somewhat surprising, especially considering the importance of the contacts between North Africa and Europe particularly, this reinforces the idea that sporadic gastroenteritis cases may be a reservoir for emerging epidemic NoV strains.
Published ahead of print on 24 December 2008. ![]()
|
|
|---|
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»