Previous Article | Next Article ![]()
Journal of Virology, August 2008, p. 7336-7345, Vol. 82, No. 15
0022-538X/08/$08.00+0 doi:10.1128/JVI.02371-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Laboratory for Infectious Diseases and Perinatal Screening, National Institute for Public Health and the Environment (RIVM), P.O. Box 1, 3720 BA Bilthoven, The Netherlands,1 Department of Microbiology, Osaka City Institute of Public Health and Environmental Sciences, 8-34 Tojo-cho, Tennoji-ku, Osaka 543-0026, Japan,2 Centre for Molecular and Biomolecular Informatics, Radboud University Medical Center, P.O. Box 9101, 6500HB Nijmegen, The Netherlands,3 Laboratory of Microbiology, Osaka Prefecture University, 1-1 Gakuen-cho, Sakai, Osaka 599-8531, Japan4
Received 2 November 2007/ Accepted 5 May 2008
|
|
|---|
|
|
|---|
Based on the genetic analysis of the RdRp and capsid regions, human NoVs can be divided into three genogroups (Gs), GI, GII, and GIV (2, 14, 39), which further segregate into distinct lineages called genotypes (2, 20, 36, 37). Recently, Kageyama et al. (20) proposed that at least 31 genotypes could be distinguished within GI and GII. The GII.4 genotype, which is represented by the Lordsdale/93/UK strain, has been the most commonly identified genotype worldwide in recent years. Genetic characterization of strains belonging to this genotype have shown a sequence of variants that have arisen over time, suggesting that rapid genetic evolution of GII.4 NoVs may in part explain their successful spread and impact on people of all ages (5, 9, 27, 29, 30, 33, 40).
Viruses belonging to other NoV genotypes are found less consistently, causing sporadic outbreaks or temporary epidemics in a limited geographic region or time period (5, 17, 23, 26). As a result, far less is known about the population structure of these variants (28). The genetic analysis of other genotype NoVs excluding GII.4 will improve our understanding of genetic evolution and its relevance for the epidemiology of NoVs.
During the spring of 2004, an epidemic of GII.2 NoV (which is represented by the Melksham/94/UK strain [Melksham]) occurred in Osaka City, Japan. Our previous study of this regional epidemic described the molecular epidemiology of these GII.2 strains (17). Here, we describe the genetic characterization of GII.2 strains from those outbreaks in comparison with viruses detected over a 12-year period in the GenBank, The Netherlands, and Japan.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. GII.2 NoV strains used in this study
|
![]() View larger version (13K): [in a new window] |
FIG. 1. Phylogenetic analysis of partial RdRp gene sequences of GII.2 NoV strains by the neighbor-joining method (A) and the Bayesian method (B). The genotypes at the RdRp regions that are not assigned numbers are represented as GII-NA. The scale indicates the number of substitutions per site. Reference strains of NoV used in this analysis are given in italics. (A) The bootstrap values are indicated on each branch. (B) Tree topology was evaluated on the basis of 1,500,000 generations. The posterior probabilities are indicated on each branch ( 0.95 of the posterior probability means that the branch has high credibility). The GenBank accession numbers of the additional reference strains used in this analysis are as follows: Hokkaido/133/03/JP, AB212306; Lordsdale/93/UK, X86557; MOH/99/HU, AF397256; Norwalk/68/US, M87661; Pont de Roide 673/04/FR, AY682549; Saitama U1, AB039775; Toronto/93/CA, U02030. For strain abbreviations, see Table 1.
|
Amplification and sequencing for the complete capsid gene of GII.2 strains. Viral RNA was extracted from stool suspensions by using a QIAamp viral RNA Mini kit (Qiagen, Hilden, Germany). RT-PCR was carried out with the reaction mixtures and enzymes as previously described (8). RT was performed at 42°C for 2 to 3 h with reverse primer, N235Rex (17), and enzyme was inactivated at 95°C for 5 min. PCR was performed using several pairs of PCR primers (Table 2) with a GeneAmp PCR system 9700 (Applied Biosystems, Foster City, CA) under the following conditions: denaturation at 95°C for 1 min; 40 cycles of 95°C for 15 s, 50°C for 30 s, and 72°C for 1 min; and a final cycle of incubation at 72°C for 5 min. When a PCR failed to produce strong products, we performed nested PCR. The amplified fragments were sequenced directly with a Big Dye terminator cycle sequencing kit and ABI 3700 sequencer (Applied Biosystems, Foster City, CA). The nucleotide sequences were determined in both orientations using the PCR primers. DNA sequences were edited using SeqManII (DNAStar Inc., Konstanz, Germany).
|
View this table: [in a new window] |
TABLE 2. Primers used to amplify the capsid gene of GII.2 NoV strains
|
Nucleotide sequence accession numbers. The nucleotide sequences determined in this study have been deposited in the DDBJ/EMBL/GenBank with the accession numbers AB281081 to AB281090 (Table 1).
|
|
|---|
83.6% nucleotide and
93.5% amino acid identities among these GII.2 strains. The phylogenetic tree based on the 1,629 nucleotide sequences of the complete capsid gene confirmed that all strains were characterized as GII.2 genotype (Fig. 2).
![]() View larger version (14K): [in a new window] |
FIG. 2. Phylogenetic analysis of complete capsid gene sequences of GII.2 NoV strains by the neighbor-joining method (A) and the Bayesian method (B). Melksham-like strains, which have a matching (GII.2) RdRp sequence, are shown in boldface. Asterisks indicate the GII.2 strains that have different (*) or unknown (**) genotypes in the RdRp region. The scale indicates the number of substitutions per site. (A) The bootstrap values are indicated on each branch. (B) Tree topology was evaluated on the basis of 500,000 generations. The posterior probabilities are indicated on each branch. The GenBank accession numbers of the additional reference strains used in this analysis are as follows: Hawaii/71/US, U07611; Hillingdon/90/UK, AJ277607; Lordsdale/93/UK, X86557; Norwalk/68/US, M87661; Toronto/93/CA, U02030. For strain abbreviations, see Table 1.
|
99.2% nucleotide and
99.1% amino acid identities). To find out if these viruses changed genetically during circulation in the community for a short period of time (about 3 months), the complete nucleotide and amino acid sequences of the capsid gene from these 14 GII.2 strains were compared. A total of 45 nucleotide substitutions were observed (Table 3). The majority were third-base position changes (77.8%) and synonymous substitutions (82.2%). These nucleotide changes resulted in amino acid changes in eight codons, five of which were located in the P2 subdomain. In the alignment of the spring epidemic strains over the 3-month period, of the 45 nucleotide substitutions, only one nucleotide position was fixed (nt 594 in the S domain) at the end of this epidemic. This did not result in an amino acid change. All amino acid changes were sporadic. Eight of the 14 strains had an identical amino acid sequences (Table 1). Of these, strain 04038 was used for further genetic analysis. |
View this table: [in a new window] |
TABLE 3. The numbers of nucleotide and amino acid substitutions in the complete capsid gene among 14 GII.2 strains detected in the spring epidemic of 2004 in Osaka City, Japan
|
|
View this table: [in a new window] |
TABLE 4. The numbers of nucleotide and amino acid substitutions in the complete capsid gene among 29 GII.2 strains collected in the GenBank, The Netherlands, and Japan over a 30-year period
|
![]() View larger version (9K): [in a new window] |
FIG. 3. Amino acid alignment of the P2 subdomains of GII.2 NoV strains showing mutations along the aligned sequences. The upper sequence alignment group (A) includes the Melksham-like strains and the lower group (B)includes the other GII.2 strains, which were recombinant genomes. In each group, sequences are arranged chronologically from top (oldest) to bottom (most recent). The detection years of the strains are indicated in parentheses. The asterisk indicates informative sites among Melksham-like strains. The arrow denotes cumulatively changing amino acid positions (342, 344, and 345) among Melksham-like strains. The numbers above the sequences indicate the sequence position relative to the position in the capsid protein of strain Melksham. Underlined sequences indicate the KGE motif that corresponds to the RGD-like motif of other NoVs and was determined with the amino acid sequence alignment of other NoVs according to Tan et al. (34). For strain abbreviations, see Table 1.
|
91.5% nucleotide and
97.4% amino acid identities among these Melksham-like strains. A total of 301 nucleotide changes were observed (Table 5), the majority of which were third-base changes (87%) and synonymous (88.7%). These nucleotide changes resulted in amino acid changes in 32 codons, half of which were located in the P2 subdomain. Twelve of 32 amino acid positions were informative (9 in the P2 subdomain) (Fig. 3). In contrast with the previous finding, several mutations were fixed: of the 12 informative sites, 2 amino acid substitutions (amino acid positions 245 and 440) in the P1 subdomain and 3 amino acid substitutions (amino acid positions 342, 344, and 345) in P2 subdomain were cumulative (Table 6), segregating the strains into five genetic groups (1994 to 1997, 1999 to 2000, 2001 to 2003, 2004, and 2005) by the neighbor-joining method (Fig. 2A) and Bayesian method (Fig. 2B).The strains detected in the spring epidemic in Osaka City had a unique sequence, with S or P residues at amino acid position 364 (Fig. 3 and Table 6). The other six informative sites were not fixed. |
View this table: [in a new window] |
TABLE 5. The numbers of nucleotide and amino acid substitutions in the complete capsid gene among 20 Melksham-like strains over a 12-year period
|
|
View this table: [in a new window] |
TABLE 6. Amino acid substitutions at six positions in the P domain among Melksham-like strains
|
![]() View larger version (49K): [in a new window] |
FIG. 4. Location of six fixed amino acid residues (positions 245, 342, 344, 345, 364, and 440, shown in red) on the monomer of the capsid protein. This 3D structure for the monomer P domain of the GII.2 NoV capsid protein was made by WHAT IF. The P1 and P2 subdomains are shown in blue and gray, respectively. The S domain is not shown. The KGE motif in the P2 domain is shown in green. (A) Predicted location of amino acid residues 342, 344, 345, and 364 in the P2 subdomain. (B) A view identical to panel A rotated 180° horizontally showing the location of amino acid residues 245 and 440 in the P1 subdomain.
|
![]() View larger version (68K): [in a new window] |
FIG. 5. The 3D structure of a dimer of the GII.2 NoV capsid protein, as predicted by WHAT IF, showing the location of the six informative amino acid residues (red). Monomer A is shown in yellow, while the monomer B is shown in blue (P1) and gray (P2). The S domains are not shown. The KGE motif (green) in the P2 subdomain on monomer B is also shown.
|
|
|
|---|
Molecular characterization of the GII.2 capsids over the full study period showed an interesting difference between recombinant genomes and nonrecombinant (Melksham-like) genomes. The Melksham-like strains clearly evolved over time, with accumulation of mutations particularly in the P2 subdomain. In contrast, GII.2 capsids from recombinant genomes did not fit this pattern and had a seemingly erratic pattern of mutations. Melksham-like strains are occasionally observed in molecular surveillance data from The Netherlands and Osaka City, and strains from Japan and The Netherlands could be consistently grouped in the same alignment, suggesting continuous widespread circulation and an ability to cause disease. The recombinant genomes with GII.2 capsids, however, occur sporadically and did not show evidence for evolution, suggesting that their circulation is limited (1, 4, 17). The finding that strains from such widely separated geographic regions were similar shows that evolution of GII.2 NoVs is a global phenomenon, similar to what has been described for GII.4 NoVs (5, 27, 33).
The molecular characterization of Melksham-like strains showed accumulation of mutations over a 12-year period, segregating them into five genetic groups. Three of these were supported in phylogenetic analysis with high bootstrap values. Two others (Fig. 2, 1999 to 2000 and 2001 to 2003) were not distinguishable in phylogenetic analysis, possibly because of the small number of isolates in the comparison. The present collection, however, contains all GII.2 strains detected in the two countries in the past 12 years of surveillance. Amino acid position 345 was the most variable in all GII.2 strains. This is suggestive of immune pressure, particularly because the amino acid is predicted to be located on a surface-exposed part of the capsid. The fixed amino acid changes were topologically clustered, judging from the 3D structure prediction, indicating that combined they may form an epitope. Immunity to NoVs has not been studied very systematically, but recent work on GII.4 NoVs has suggested that these viruses evolved under immune pressure (9, 30, 33). The lack of a tissue culture model makes it difficult to corroborate these findings by using cross-neutralization tests with hyperimmune serum directed against specific variants. Evolution of strains in an immunocompromised patient with chronic shedding of NoVs has been demonstrated by Nilsson et al. (28), who suggested that the cumulative amino acid substitutions appearing in the P2 subdomain were immune response driven. To clarify the relations between these amino acid mutations and receptor or antibody binding further immunological studies are needed.
Recently, new variants of the predominant NoV, the GII.4 genotype, appeared in 2002 (27), 2004 to 2005 (5), and 2006 (33) with global distribution. These new GII.4 variants were characterized into different clusters by the phylogenetic analysis based on the capsid gene (5, 9, 30). Okada et al. (30) reported that GII.4 strains had cumulative amino acid changes in the P domain and 5% to 18% amino acid diversity in the P2 subdomain between GII.4 subtypes for 5 years study. Bull et al. (5) also described more than 5% amino acid diversity in the complete capsid gene between past GII.4 strains and variants. Siebenga et al. (33) found 9% fixed amino acid mutations across the capsid gene over a 12-year period, but only five of these consistently delineated subsequent epidemic strains. In contrast, here we found only 2.6% amino acid diversity across the complete capsid gene of Melksham-like strains in a 12-year period, showing a clearly lower rate of change. The difference between these genotypes is intriguing and suggests clear differences in epidemiology that are not easy to explain. This study has been done using a limited number of strains in a rare genotype from a limited number of locations. It is necessary to continue molecular surveillance of NoV infections to further the evolutional analysis of each genotype.
In conclusion, the present study showed that NoVs belonging to a rare genotype evolved by the accumulation of mutations in the surface-exposed parts of the capsid, suggesting immune response-driven evolution. The pattern of change similar to what has been observed for GII.4 is a global one, suggesting that these viruses circulate much more than is apparent from their presence in surveillance data. The data also showed for the first time that recombinant genomes of NoVs behave differently, even when they share the same capsid genes as nonrecombinant genomes. Understanding how NoVs evolve is necessary for finding more effective ways to control this disease and particularly its impact in healthcare settings.
This work was supported by the European Commission, DG Research Quality of Life Program, 6th Framework (EVENT, SP22-CT-2004-502571), and by a grant for Research on Food Safety of the Ministry of Health, Labour and Welfare, Japan.
Published ahead of print on 14 May 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»