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

Department of Pediatrics, Fukushima Medical University School of Medicine, Fukushima 960-1295,1 Research and Development Department, Mitsubishi Kagaku Bio-Clinical Laboratories, Inc., Tokyo 174-8555,2 Fukushima Institute for Public Health, Fukushima 960-8560,Japan3
Received 5 April 2006/ Returned for modification 5 July 2006/ Accepted 21 October 2006
|
|
|---|
|
|
|---|
HEV71-related illness is more severe, with a significantly greater frequency of serious complications and fatality, than illness caused by CV-A16 (2). In 1997, deaths associated with epidemics of HEV71-associated HFMD in Sarawak, Malaysia, followed by outbreaks with high mortality in Taiwan in 1998 and 2000, raised considerable public concern about the virulence of this virus. Since then, several groups have attempted to describe the molecular epidemiology of HEV71 in the Asia-Pacific region and have reported the relationships of HEV71 epidemics with the genetic diversity of HEV71 strains (1, 5, 10). The results indicate that HEV71 strains causing HFMD outbreaks were genetically changed.
On the other hand, the molecular epidemiology of CV-A16 associated with HFMD epidemics has not been fully described (9). In the present study, we evaluated the relationship between the chronologic CV-A16 epidemics in a restricted area, i.e., Fukushima Prefecture, Japan, and the genetic diversity of the CV-A16 strains. We also examined the geographic genetic relationship between the CV-A16 strains isolated in Fukushima and those isolated in other areas of Japan and China, using phylogenetic analyses constructed using the neighbor-joining method on the basis of the VP4 and VP1 sequences.
|
|
|---|
PCR and sequence determination of VP4 gene. Randomly selected isolates (63 of 241 strains) from 1983 to 1999 and all isolates (69 strains) from 2000 to 2003 were used for further genetic analysis. A total of 132 strains were isolated from patients with HFMD. The methods of molecular diagnosis of enteroviruses by nested reverse transcription-PCR (RT-PCR) and phylogeny-based classification using the VP4 sequences are described elsewhere (7). Briefly, viral RNA was directly extracted from 100 µl of the stock virus samples using the Smitest R kit (Genome Science Laboratories) according to the manufacturer's instructions. The RNA was dissolved with 10 µl of RNase-free distilled water containing 40 U of RNase inhibitor (RNasin; Promega) and 50 pmol of a reverse primer, OL68-1 [nucleotides (nt) 1178 to 1197, 5'-GGTAA(C/T)TTCCACCACCA(A/G/C/T)CC-3']. The positions of the primers for RT-PCR were numbered according to the complete nucleotide sequence of the attenuated poliovirus Sabin 1 strain (11). The RNA was subjected to heat denaturation for 15 s at 100°C, followed by snap-cooling in an ice-water bath. Reaction mixture (10 µl; 200 U of Moloney murine leukemia virus reverse transcriptase [Life Technology], 2.5 mM deoxynucleoside triphosphates, and 40 U of RNasin [Promega]) was added to each RNA sample. cDNA synthesis was performed for 1 h at 37°C. In total, 5 µl of the cDNA reaction mixture was added to 45 µl of 1x Taq buffer containing 12.5 pmol of a forward primer, MD91 (nt 444 to 468, 5'-CCTCCGGCCCCTGAATGCGGCTAAT-3') and 2.5 U of Taq DNA polymerase (Roche Diagnostic Systems). Seminested PCR was performed using 5 µl of the PCR product with a pair of primers, EVP4 (nt 541 to 560, 5'-CTACTTTGGGTGTCCGTGTT-3') and OL68-1. After initial denaturation at 94°C for 5 min, 40 cycles of amplification were performed by using the GeneAmp PCR System 9600 (PE-Applied Biosystems). Each cycle consisted of denaturation at 95°C for 30 s, primer annealing at 55°C for 30 s, and an extension reaction at 72°C for 1 min. The 40 amplification cycles were followed by a final extension at 72°C for 7 min. The first PCR primer pair allowed for amplification of 753 bp, and the seminested PCR primer pair allowed for amplification of 656 bp. The final PCR products encoded the 3' one-third of the 5' nontranslated region, the entire VP4, and the 5' one-third of VP2 of each CV-A16 strain. The PCR products, including the entire VP4 sequences, were separated in 1% agarose gels and purified with a QIAquick gel extraction kit (QIAGEN). The nucleotide sequence was determined by using a 373A DNA autosequencer (PE-Applied Biosystems) with fluorescent dideoxy chain terminations (PE-Applied Biosystems) and EVP4 and OL68-1 primers.
Phylogenetic analysis based on the VP4 gene. The entire VP4 nucleotide sequences of the 129 strains isolated in Fukushima were determined and used for phylogeny-based analysis, along with those of 64 prototype enterovirus strains. We estimated the evolutionary distances using the Kimura two-parameter method (8) and constructed unrooted phylogenetic trees with the neighbor-joining method (14). Bootstrap analysis was performed by resampling the data sets 1,000 times. Bootstrap values greater than 70% were considered to be statistically significant for the grouping. The VP4 sequences of representative 54 CV-A16 strains isolated in Fukushima were also compared to those of 20 strains isolated in other areas of Japan, 29 strains in China, and 1 strain in the United Kingdom taken from international databases (GenBank) by using phylogenetic analysis.
PCR and phylogenetic analysis based on the VP1 gene. To confirm the relevance of genogrouping based on the VP4 region, we used the phylogeny-based classification methods using the VP1 sequences that were reported by Oberste et al. (12). Twenty-five CV-A16 strains isolated in Fukushima were randomly selected, and their VP1 regions were amplified using two sets of primer pairs, i.e., 055/VP3-011/2A and 055/VP3-009/2A. The VP1 sequences were determined and compared by using phylogenetic analysis with those of 28 strains isolated in China and 1 Taiwan strain taken from international databases.
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 1. Numbers
of CV-A16 isolates in Fukushima Prefecture and in Japan between 1983
and 2003. Numbers of CV-A16 isolates in Fukushima (closed circles and
solid line) and in Japan (open circles and broken line) are expressed,
as reported to the Infectious Disease Surveillance Center in Japan by
prefectural and municipal public health institutes through the Japanese
infectious agents surveillance program. A, B, and C refer to genogroups
designated based on phylogenetic analysis of the VP4 sequences of
CV-A16 strains. Solid bars and broken bars indicate HFMD
outbreaks and epidemic periods, respectively, due to corresponding
genogroups.
|
![]() View larger version (35K): [in a new window] |
FIG. 2. Phylogram
depicting the phylogenetic relationships on the basis of the VP4
sequence among 129 CV-A16 strains isolated in Fukushima from 1983 to
2003. Bootstrap analysis was performed by resampling the data sets
1,000 times. Bootstrap values greater than 70% were considered to be
statistically significant for the grouping and were denoted in the
figure. Isolated place, strain name, and isolated year were indicated.
The GenBank accession number of one isolate to represent isolates with
identical sequences is also indicated in parenthesis. CA-A16/G-10/51 is
the prototype CV-A16 strain. The VP4 nucleotide sequence of prototype
HEV71/BrCr/71 was used as an outgroup in the
analysis.
|
![]() View larger version (45K): [in a new window] |
FIG. 3. Phylogram
depicting the phylogenetic relationships on the basis of the VP4
sequence among representative 54 CV-A16 strains isolated in Fukushima
and 50 CV-A16 strains isolated in Japan, China, and United Kingdom
taken from international databases (GenBank). Bootstrap analysis was
performed by resampling the data sets 1,000 times. Bootstrap values
greater than 70% were considered to be statistically significant for
the grouping and were denoted in the figure. Isolated place, strain
name, and isolated year were indicated. GenBank accession numbers of
CV-A16 strains isolated in Fukushima and the strains taken from
international databases are indicated in parentheses. CA-A16/G-10/51 is
the prototype CV-A16 strain. The VP4 nucleotide sequence of prototype
HEV71/BrCr/71 was used as an outgroup in the
analysis.
|
![]() ![]() ![]() View larger version (176K): [in a new window] |
FIG.4. Alignment
of the entire VP4 sequences. Common genetic diversities were observed
in each genogroup. Nucleotide differences from the prototype
CV-A16/G-10/51 strain observed in genogroup A, genogroup A and B, or
genogroup A and C are indicated by shading and those in genogroup B,
genogroups B and C, or genogroup C are indicated in
boxes.
|
|
|
|---|
Phylogeny-based classification using the VP4 sequence is useful for the identification of human enteroviruses (3, 4). The method takes advantage of the detection of the divergence in VP4 sequences both between and within serotypes, and thus is also of use for global epidemiologic studies of enteroviruses (1, 5, 7). We investigated the genetic diversity of CV-A16 strains associated with HFMD epidemics in the Fukushima Prefecture, Japan, from 1983 to 2003 and compared their genetic relation to those isolated in other areas of Japan and in China using the same method. CV-A16 strains isolated in Fukushima, from 1983 to 2003, made at least three distinct clusters on the phylogenetic tree. The three clusters were designated genogroups A, B, and C and were associated with the HFMD epidemics from 1984 to 1994 (including the 1985 and 1991 outbreaks), 1987 to 1998 (including the 1988 and 1998 outbreaks), and 1995 to 2003 (including the 1995 and 2002 outbreaks), respectively. The predominant genogroup was replaced with a new genogroup. CV-A16 strains within a genogroup gradually became genetically divergent in several epidemics. These results demonstrated that the introduction of a new genogroup in addition to the genetic divergence within a genogroup resulted in repeated HFMD outbreaks in Fukushima, Japan.
Each genogroup formed the same cluster with CV-A16 strains isolated during essentially the same time period in other areas of Japan and in China. The clustering seemed to be more closely related to the date of isolation than to the geographic location. Genetic diversities appeared randomly throughout the VP4 gene and were common to each genogroup. These indicated that a new CV-A16 genogroup might derive from other regions of the world, be predominant for several years with genetic divergence, and then disappear.
The protein encoded by the VP1 gene is the most exposed and immunodominant of the capsid proteins. Serotypic identification and classification rely on antigenic methods, and VP1 gene sequence data are likely to give the most useful information in molecular epidemiologic studies (13). Shimizu et al. analyzed the phylogeny of HEV71 isolates, which is the other major etiologic agent of HFMD, based on the nucleotide sequence alignment of both the VP1 and VP4 regions (15). Phylogenetic trees based on VP4 sequences revealed that the HEV71 strains isolated from the Western Pacific Region formed two major genogroups, B and C, and were similar to those based on the VP1 sequences. Cardosa et al. reported that the VP1 and VP4 gene sequences both provide similar phylogenetic information, but the higher bootstrap values in the VP1 dendrograms provide greater confidence, particularly when elucidating new genotypes. Thus, the use of the shorter VP4 gene might be helpful for HEV71 surveillance, but the VP1 gene should be used for confirming data obtained with VP4-based analysis (1). Therefore, we confirmed the data obtained from the analysis based on VP4 genes of CV-A16 isolates by using VP1-based analysis. Phylogenetic analysis based on the VP1 gene demonstrated the existence of three genetically distinct groups, and the grouping was completely identical to the result obtained with VP4-based analysis. Bootstrap values for the grouping of A, B, and C based on the VP1 gene were higher than 80%, whereas those for the grouping of A and C based on the VP4 gene were higher than 70%. Although bootstrap values are slightly higher in the VP1 dendrograms than in the VP4 dendrograms, the VP4 gene as well as the VP1 gene seems to be appropriate for the epidemiologic study of CV-A16. In conclusion, phylogenetic analysis based on the VP1 gene confirmed that CV-A16 strains isolated in Fukushima formed three genogroups, and each genogroup formed the same cluster with CV-A16 strains isolated during essentially the same time period in China.
In summary, we describe the evident genetic diversity and changes in the VP4 protein region of CV-A16 strains that were isolated in a restricted region through more than 20 successive epidemics. Our results indicate that CV-A16 strains causing HFMD had genetically changed twice during the period and that these CV-A16 strains might have been transmitted overseas. We conclude that the repeated outbreaks of CV-A16-related HFMD might be caused, in part, by the worldwide transmission of the genetically changed CV-A16 strains, as well as a large cohort of nonimmune individuals. To test this hypothesis, a worldwide surveillance system for HFMD and genetic analysis of isolated CV-A16 strains is necessary.
Published ahead of print on 8 November 2006. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»