Previous Article | Next Article ![]()
Journal of Virology, August 2004, p. 8349-8358, Vol. 78, No. 15
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.15.8349-8358.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
National VZV Laboratory, Herpesvirus Section, Division of Viral and Rickettsial Diseases, National Center for Infectious Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia,1 National Microbiology Laboratory, Winnipeg, Canada,2 Virology Department, Ruprecht-Karls University, Heidelberg, Germany,3 Department of Medical Virology, Landspitali University Hospital, Reyjavik, Iceland4
Received 16 October 2003/ Accepted 24 March 2004
|
|
|---|
|
|
|---|
During the past 2 decades, several groups attempted to evaluate VZV phylogeny. Early efforts in VZV typing used DNA restriction fragment length polymorphism (RFLP) analysis (13, 37, 38), an approach that confirmed the identity of the VZV strain that causes varicella on primary infection and later reactivates to cause zoster. Relatively consistent restriction enzyme digestion profiles for different VZV strains were observed, providing the first evidence that VZV has a highly conserved genome. Intrastrain variation in restriction enzyme fragment profiles among wild-type VZV isolates was observed (22, 36, 37). However, the most prominent differences were linked to variation in the number and composition of VZV genome repeat elements (23, 39, 40, 41). A number of RFLP strategies to distinguish VZV isolates from different geographic regions and to differentiate between the live-attenuated Oka vaccine strain and wild-type viruses have been proposed (1, 5, 23-28, 39, 41). It has since been determined that neither the repeat regions R1 to R5 nor the combination of ORF38 and ORF54 single-nucleotide polymorphisms (SNPs) reliably discriminate between VZV wild-type strains and the Oka vaccine strain (20, 30).
Nevertheless, two restriction enzyme cleavage sites corresponding to these SNPs have served as powerful reference points for characterizing both epidemiologic and geographic strain variation in VZV. Most VZV isolates in the United States contain a PstI restriction enzyme site in ORF38 that distinguishes them from a subset of Japanese strains, such as the Oka vaccine strain, that lack this restriction site (24). In addition, a C-to-T change at position 95241 leads to the elimination of a BglI restriction enzyme site in ORF54 that is found in all Japanese strains and that is absent in most U.S. VZV isolates (1, 28-31). As such, Japanese strains are either PstI+ BglI+ or PstI BglI+, while most isolates in the United States are PstI+ BglI. The latter genotype clearly predominates in eastern Australia, the United Kingdom, and the United States (21, 29, 36, 38). Strains bearing the Japanese PstI+ BglI+ genotype are also common in countries with a history of European colonization, and this is the predominant genotype in equatorial Africa, India, Bangladesh, China, and western Australia (29-31, 33-34). Furthermore, in eastern London, in which a large number of Indian and Bangladashi immigrants reside, the percentage of BglI+ VZV strains increased from 10% in the 1980s to more than 30% in the 1990s (21). At first, it was believed that the PstI BglI+ (Oka vaccine strain) genotype was no longer circulating, but it has recently been identified among isolates from Japan, the western United States, and Hawaii (30, 31, 40). VZV strains with a PstI BglI profile have never been observed. In summary, RFLP analysis of VZV strains using the above approach identified two major groups. The first group includes strain Dumas and the typical European, North American, and Australian PstI+ BglI isolates, whereas the second group harbors Asian (including Japanese), African, and western Australian PstI BglI+ isolates.
RFLP analysis has also been performed using long-distance PCR products, facilitating the classification of 14 Japanese clinical specimens into nine subgroups (41). A separate study of VZV strains isolated from patients in Thailand used restriction endonuclease analysis with BglI, PstI, EcoRI, SmaI, and BamHI to show that VZV strains circulating in Thailand are distinct from those circulating in Japan. Sixteen different genotypes were evident, including those of strains with the BglI restriction site characteristic of Japanese (42) and African-Indian (30, 34) isolates. The high number of variants identified using this approach undoubtedly resulted from its reliance on genome repeat element copy number and the use of multiple endonucleases. Nevertheless, Japanese and Thai VZV strains displayed both common and distinctive markers. A milestone finding obtained by using antipeptide antibodies to discriminate Japanese and European strains was also reported by Kinchington and Turse (26).
RFLP studies of VZV mutability in serial culture have produced inconsistent observations. While some investigators have documented changes in VZV restriction profiles on serial passage (11, 12, 20, 45, 47), others failed to observe any differences after as many as 33 passages in cultured cells (20, 45). The serial passage of the Oka parental virus in tissue culture resulted in a well-documented accumulation of mostly single point mutations (20). The completion of the genomic sequencing for VZV Oka vaccine and parental strains (20) has opened new horizons not only for the definitive differentiation of VZV vaccine strains from wild-type viruses but also for the genotypic analysis of wild-type VZV strains. The direct comparison of these data with the published genome sequence for the European Dumas strain has yielded important information. For example, the sequences that vary between the Oka vaccine and parental strains are concentrated in a single ORF, ORF62 (20). In addition, sequence polymorphisms were shown to be extremely rare among VZV wild-type strains and to be generally limited to individual point mutations distributed across the entire genome. Last, the differences between the Japanese strain (Oka) and the European strain (Dumas) provided a blueprint of mutations valuable for genotyping. Given these observations, it is not surprising that conventional RFLP methodology was ineffective at distinguishing these minor variations. Almost all of the sequence changes in the VZV Oka vaccine strain genome are SNPs, the detection of which requires the use of a series of specific restriction enzymes.
To date, three research groups, including our own, have advanced the development of a strategy for practical VZV genotyping (14, 33, 34, 44; V. N. Loparev and D. S. Schmid, Abstr. 26th Int. Herpesvirus Workshop, abstr. 3.15, 2001). As a first step toward a refined strain surveillance methodology, we evaluated VZV genome variability to establish a basis for the stable classification of wild-type strains into major viral genotypes. Genotype-based classification schemes are feasible even in the absence of distinct phenotypes, as has been amply demonstrated for other human herpesviruses (1, 4, 7, 9, 32, 43, 46). Nevertheless, careful genotypic analysis could lead to the identification of virulence factors and hence to an improved understanding of VZV biology (35). Another important issue that may be addressed using this approach is the assessment of the direction and predisposition of virus phylogeny in VZV-vaccinated and unvaccinated populations. We report here the development of a novel test for the simple typing of wild-type strains of VZV into three major circulating genotypes, which we have designated E (European), J (Japanese), and M (mosaic), using SNPs identified in a 447-bp region of ORF22. Altogether, 323 globally distributed VZV isolates obtained from recent clinical infection and representing currently circulating strains in Europe, Asia, Africa, Oceania, and the Americas could all be readily assigned to one of three genotypes. Furthermore, regional specificity for each genotype was evident from this analysis. These data shed light on evolutionary relationships among these major genotypes and their geographic clustering. We anticipate that these observations will help to predict future evolutionary patterns and to track the changes in VZV epidemiology that inevitably accompany active vaccination programs.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. VZV strains analyzed for this study
|
All samples initially were characterized by the use of conventional PCR targeting VZV-specific sequences in ORF38, -54, and -62 in a 50-µl reaction mixture, followed by restriction enzyme analysis as previously described (28, 29) or by fluorescence resonance energy transfer (FRET) probe-based annealing and melting curve evaluation (30). A PCR/RFLP method for detecting SNPs in ORF38 and -54 was adapted to FRET-based real-time PCR, using the published primers for amplification (28). Ten-microliter aliquots of the resulting amplicon were mixed with 10 µl of LightCycler probes (CDC) in a mixture of 1x SSC (0.15 M NaCl plus 0.015 M sodium citrate; Difco Laboratories, Detroit, Mich.), 1x bovine serum albumin (New England Biolabs, Beverly, Mass.), and 3 mM MgCl2 (Applied Biosystems, Foster City, Calif.), and the melting point for each probe was detected. The probes 5'-GGA CTT GAA GAT GAA CTT AAT GAA GCC CGT GA-3' (anchor probe; positions 69393 to 69362) and 5'-ACG ATA TAT ACC GCA GTT GTT GCG GTA-3' (sensor probe; positions 69360 to 69334) were used for sequence evaluation of the PstI site region, in which the melting curve peak differed from that of the PstI+ VZV strain Webster by 3°C. The probes 5'-CAT TTT GCA TAC ACT CAA CTA GGC TTG TGA-3' (anchor probe; positions 95201 to 95230) and 5'-AAC CGC CGC TCC TCT GG-3' (sensor probe; positions 95232 to 95248) were used for evaluating the BglI site region, with a 4°C difference between the melting curve peaks of the BglI+ VZV Oka vaccine strain and the BglI VZV strain Webster. The five pairs of FRET probes for genotyping are depicted in Table 2. The anchor probes were 5' labeled with LC640 dye and 3' phosphorylated. The sensor probes were 3' labeled with fluorescein, and the resulting data were used for preliminary genotyping (Fig. 1). All nucleotide positions for the sequence polymorphisms, primers, and genes indicated in this paper correspond to the reference Dumas strain sequence (GenBank accession number gi: 9625875). For all melting curve analyses with FRET probes, the LightCycler (Roche) real-time PCR instrument was used.
|
View this table: [in a new window] |
TABLE 2. Probes used for genotyping
|
![]() View larger version (81K): [in a new window] |
FIG. 1. Analysis of genotypic variation using data from multiple VZV ORFs. The alignment of SNP in ORF1, -22, -31, -37, -54, -67, and -62 was performed for clinical isolates of VZV from around the world. RFLP of Pst+ Bgl, Pst+ Bgl+, and Pst Bgl isolates was also performed. Numbers in the top rows indicate ORF and position of the SNPs in the genome, with the European VZV strain Dumas as a reference for the base pair numbers. The ORF54 SNP is the site that has been used in some protocols to discriminate the Oka vaccine strain from most wild-type strains. Green indicates mutations common to the E genotype reference strain; yellow indicates mutations in common with the J genotype Oka parental strain; red indicates strain-specific individual mutations; white indicates bases common to both the J and E reference strains.
|
For DNA purification, PCR products were electrophoresed on 2% precast agarose gels (BMA, Rockland, Maine); the 447-bp band was cut out, and DNA was extracted with a gel extraction kit (QIAGEN, Valencia, Calif.). To prepare amplicons for the sequencing procedure without gel purification, we used ExoSAP-IT reagent (USB Corporation, Cleveland, Ohio).
DNA sequencing and analysis. Automated DNA sequencing was performed with the ABI PRISM 377 genetic analyzer (Perkin-Elmer, Foster City, Calif.). Sequencing reactions of gel-purified PCR products were performed by using the BigDye Terminator, version 3.0, cycle sequencing ready reaction mixture (Applied Biosystems) according to the manufacturer's instructions. Either the p22R1f or the p22R1r primer (3.2 pmol) and 10 ng of PCR product as a template were used. The GeneAmp 9700 thermal cycler (Applied Biosystems) was set to 25 cycles of 96°C for 10 s, 50°C for 5 s, and 60°C for 4 min. The resulting products were purified on spin columns (Princeton Separations, Adelphia, N.J.), dried down, and resuspended in template suppression buffer (Applied Biosystems). Cycle sequencing reaction products were fractionated on an ABI 377 analyzer (Applied Biosystems). Primary DNA sequence assembly and analysis were performed with Sequencher (Gene Codes Corp., Ann Arbor, Mich.), and the sequence was compared with that of VZV strain Dumas.
Altogether for this study, 326 strains were used to amplify and sequence the 447-bp fragment, for a total 145,722 bp (Table 1), and about 10,000 bp was sequenced in various ORFs for 16 strains (160,000 bp), for a grand total of 305,722 bp worth of DNA sequence data used in these analyses. This is the equivalent of 2.446 complete VZV genomes. Targeted sequences included those of ORF1 (positions 915 to 592), ORF22 (positions 34083 to 42371), ORF31 (positions 57008 to 59611), ORF37 (positions 66074 to 68596), ORF54 (positions 95984 to 93678), ORF67 (positions 114496 to 115560), and ORF62 (positions 109133 to 105204).
The obtained DNA sequences were compared with the Dumas and Oka parental strain sequences (data not shown) in order to locate SNPs. The sample-derived sequences were aligned with the corresponding regions of the reference VZV strains by the PileUp program of the Wisconsin GCG Sequence Analysis Package, version 10 (Genetics Computer Group, Madison, Wis.) or MegAlign (DNAstar package; DNASTAR Incorporated) into pairwise and multiple sequence alignments of DNA with the Clustal V method. MegAlign was also used to create phylogenetic trees. The Clustal method groups sequences into clusters by examining sequence distances between all pairs. Clusters are aligned as pairs and then collectively as sequence groups to produce the overall alignment. After the multiple alignment is completed, a neighbor-joining method is employed to reconstruct phylogeny for the putative alignment. A neighbor-joining tree was calculated to illustrate the distribution of intergenotypic differences.
In a final short version of VZV genotyping we described VZV genotypes as lineages of VZV distinguished by differing in least two out of four reference mutations in the 447-bp region of ORF22.
|
|
|---|
Preliminary analysis indicated that Pst+ BglI viruses represent most (75 to 90%) strains collected in United States and Canada (Table 1). In a representative subset of 52 U.S. strains presented in Table 1, 86% were determined to have the PstI+ BglI profile. This genotype was also the most commonly isolated genotype in Canada (76%), Europe (95%), and the Asian portion of the Russian Federation (100%). All specimens collected from eastern Australia were also PstI+ BglI strains.
In contrast the PstI+ BglI+ genotype represented a minority of U.S. isolates (8%). As illustrated in Table 1, this genotype was common among isolates from India (100%), south China (100%), Bangladesh (100%), Nepal (100%), DRC (100%), and Chad (100%). Seventy percent of the strains collected in western Australia were also of this genotype.
We also confirmed earlier reports (28, 29) that most modern Japanese strains have the PstI+ BglI+ restriction enzyme site profile (100% of 20 isolates examined). Of interest, three wild-type VZV strains with the PstI BglI+ profile of the Japanese Oka vaccine strain genotype were isolated in Hawaii in 2002. All PstI/BglI genotyping was performed in parallel, by both RFLP analysis and melting point analysis with FRET probes, with 100% observed agreement between the two methods. Given that 326 specimens of diverse geographic origins were evaluated, this study confirms the validity of using LightCycler-based FRET analysis to identify SNPs.
Targeted sequence analysis in multiple VZV ORFs for selected representative VZV isolates. From the CDC repository of more than 400 VZV isolates obtained from countries all over the world, a subset of 16 VZV isolates representing all of the PstI/BglI genotypes (identified through RFLP analysis) were selected for extensive targeted sequence analysis. Direct comparison of newly acquired and previously reported VZV DNA sequence data (14, 20) documented very limited variation in genes encoding the envelope glycoproteins gB, gE, gI, gH, and gL and in ORF1, -4, -10, -22, -33, and -62. However, some common patterns became apparent and are displayed in Fig. 1.
Most VZV strains collected in the United States had the greatest overall similarity to the reference European strain Dumas. In a like manner, all Japanese strains exhibited maximum similarity to other Japanese isolates and differed the most from Dumas and most U.S. isolates. This was true for both Pst and Pst+ Japanese strains.
A number of strains isolated from patients in various countries, particularly Africa and Indochina, carry a combination of Japanese (J) genotype-like and European (E) genotype-like mutations alternately across the entire genome. Thus, they feature a mosaic of J genotype-like and E genotype-like SNPs. This pattern was consistently observed for specific mutations analyzed in ORF22, -31, -37, and -62 (Fig. 1) and for 27 additional mutations in different VZV genomic loci (data not shown). ORF1 sequences were identical for Japanese and African-Chinese strains (Fig. 1) and, in addition, invariably shared the J-genotype-specific SNP (Bgl+) in ORF54. In contrast, the ORF62 sequence of those same African-Chinese strains had closest homology to E genotype isolates.
Several strain-specific mutations were also observed (positions 560, 38019, and 38229) and are presented in Fig. 1. Further analysis of VZV strains and additional targeted sequence analysis will be required to clarify how useful these more restricted markers will be for characterizing regional genotypic variation.
VZV genotyping using sequence variation in ORF22. Genomic analysis of mutations located in seven VZV ORFs permitted the categorization of all of our clinical isolates into three distinct genotypes: E, J, and M (mosaic, combined type). Given that the quantity of DNA recovered from clinical isolates is often limited and that assessing such a large number of SNPs is laborious, we sought a simplified approach to performing the preliminary genotyping of VZV strains.
We identified four polymorphic positions in a short 447-bp fragment of ORF22 that could be used to robustly categorize circulating strains of VZV into one of the three genotypes. In a manner similar to the expanded-genotyping approach used in Fig. 1, a simplified protocol using a single ORF22 amplicon unambiguously segregated strains into the E, J, and M major genotypes (Fig. 2). All of the VZV isolates that appear in this study were analyzed by this method. Single representatives of every variation that we observed are listed in Fig. 2. Several variants of M genotype were also identified by this method. As illustrated in Fig. 2, the M genotype strain California123 contained more E type variations (three of four) than J type variations (one of four). By contrast, the Mexico_68, DR (Morocco), and Iceland_006 strains contained more J type variations (three of four) and only one of four E type variations in the analyzed region. Additional strain-specific mutations at position 38019 and position 38229 in Mexico_68 and Iceland_006, respectively, may provide the basis for further differentiation into M subtypes.
![]() View larger version (43K): [in a new window] |
FIG. 2. Analysis of genotypic variation using the 447-bp ORF22 site. Nucleotide diversity found in the 447-bp fragment of ORF22 was used to type all 326 VZV isolates displayed in Table 1. Dumas and Oka parental strains are included as the reference standard for E and J genotypes, respectively. All E and J strains displayed the same SNP patterns. Some variability was apparent among M genotype strains. The examples displayed represent all of the VZV variations that were observed; the figure has been condensed to make it more readable. The color scheme for the boxes is as for Fig. 1.
|
![]() View larger version (23K): [in a new window] |
FIG. 3. Genetic relatedness of global VZV strains. Phylogenetic analysis of ORF31 and -22 of VZV strains was performed by the neighbor-joining method (Clustal). All different structural variants detected in this study are included (16).
|
The E genotype was also most prevalent in Asian Russia, Kazakhstan, eastern Australia, Canada, and the United States. E type viruses were also isolated in Jordan, Chile, and Nicaragua (Fig. 1). On the basis of analysis of numerous isolates obtained from widely separated geographic locales, the E genotype was shown to form a conservative group based on the profile of ORF22 SNPs, and all E strains also displayed the PstI+ BglI marker pattern. One surprising observation was the discovery of geographic dominance of M genotype strains in tropical and subtropical areas such as Africa (DRC, Chad, Morocco) and in Asia (India, Bangladesh, Nepal, and China). For this aspect of the study, we added a number of additional strains more recently obtained from Brazil, Argentina, Cote D'Ivoire, Ethiopia, Zimbabwe, South Africa, Thailand, Vietnam, north China, Great Britain, Tajikistan, Kazakhstan, Uzbekistan, Kyrgyzstan, Turkmenistan, and Azerbaijan. Above 30oN and below 30oS, E strains and J strains represent the vast majority of circulating VZV strains, whereas in the midlatitudes framed by those lines of latitude M strains predominate (Fig. 4). A substantial fraction of western Australian isolates also had the M genotype, exclusively in the portion of Australia that lies above the 30th parallel. All of these M strains contained the BglI site in ORF38.
![]() View larger version (35K): [in a new window] |
FIG. 4. Global distribution of the three major genotypes. Marked areas indicate the major circulating genotype (>70%) for the region. The 326 strains used in previous analysis (Table 1) are all represented; additionally, a number of other more recently acquired strains were analyzed only by ORF22 genotyping: Brazil, 10 isolates (8 M, 2 E); Argentina, 6 E isolates; Cote D'Ivoire, 5 M isolates; Ethiopia, 5 M isolates; Zimbabwe, 3 M isolates; South Africa, 3 E isolates; Thailand, 1 M isolate; Vietnam, 2 M isolates; north China, 3 E isolates; Great Britain, 5 E isolates; Tajikistan, 2 E isolates; Kazakhstan, 2 E isolates; Uzbekistan, 2 E isolates; Kyrgyzstan 2 E isolates; Turkmenistan, 2 E isolates; Azerbaijan, 2 E isolates. (Modified with permission, Corel Corp., copyright [copy] 2004.)
|
All amplicons from the targeted region in ORF22 were tested with the FRET LyghtCycler probes. As a result of analyzing 326 samples (Table 1), we identified 201 E type, 88 M type, and 28 J type strains and confirmed the practicality of VZV genotyping without resorting to DNA sequencing.
|
|
|---|
Genotyping is a useful epidemiologic tool for many human herpesviruses, although it has not generally led to correlations between sequence variation and virus pathology (7, 8, 9, 17, 32, 43, 46). Genotyping protocols intended for general laboratory use must be simple to perform. They should also be parsimonious and able to yield results from the limited amount of DNA typically recovered from clinical specimens, e.g., vesicular swabs or scabs. An ideal format permits reliable VZV genotyping using a single short amplicon or, if necessary, a very limited number of them.
We searched for short sequences containing multiple SNPs that consistently identified the E, J, and M genotypes observed by using broader SNP data (ORF1, -22, -31, and -62). The VZV genome is extremely stable; all three genotypes of VZV differ in DNA sequence by only about 0.2%. The biological properties of VZV help to confer genomic stability. VZV is highly infectious and causes large outbreaks in naive populations, a characteristic that flattens the molecular divergence of VZV strains. VZV establishes lifelong latency, but subclinical reactivation is believed to occur less frequently than with other alphaherpesviruses (14). The low number of replicative cycles in the infected host (14) is thus thought to restrict opportunities to introduce mutations. Reactivated VZV is sometimes transmitted to susceptible individuals, causing varicella and periodically reintroducing older strains into the population.
Fortuitously, we identified a 447-bp sequence proximal to the carboxy terminus-encoding portion of ORF22 that satisfied our need for detecting viral DNA in low-titer clinical specimens and for reliably discriminating the three major genotypes. Both the multiple-ORF SNP analysis and the ORF22 fragment SNP patterns unerringly resolved the same stable VZV genotypes, all three of which are likely to have emerged in the distant past (19). We also noted that the ORF22 region selected for genotyping was stable in highly passaged European (Webster) and Japanese (Oka) strains.
Geographic distribution of the three ORF22 genotypes. Regional VZV strain dominance could reflect the immunogenetic composition of the local human population in which it originally evolved (2, 6, 19). This might explain the occurrence of a uniform genotype in relatively isolated and conservative populations, e.g., Japan and the former Soviet Union. In geographic locations with populations in which large-scale immigration from multiple regions occurred, such as the United States, western Europe, and Australia, the distribution of VZV genotypes reflects that demographic quality. In support of this view, we have isolated J strains in North America exclusively from areas with the highest rates of immigrants from Asia.
We could also differentiate European and North American strains from Japanese, Asian, and African strains based on the ORF54 BglI site alone. Thus, ORF54-based sequence analysis could be useful for preliminary genotyping in the United States, Canada, much of Europe, eastern Australia, and South America, where M and J VZV strains are not as common. The routine differentiation of J and M strains requires DNA sequencing, RFLP analysis, or real-time hybridization methods with fluorescent probes, as described in this report.
In addition to discriminating the major circulating genotypes of VZV, the analysis of the SNPs in the 447-bp fragment also permitted further subtyping of M genotype strains. Further analysis of genomic sequence is likely to reveal additional regions that will enable finer discrimination of regionally distributed VZV strains, a process that should be accelerated by our recent determination of a genomic sequence for an M strain.
Correlation with other polymorphic genomic markers. Genotypic analysis of the N-terminal part of the VZV gB gene (ORF31) revealed a highly conserved amino acid sequence and 99.88% homology at the nucleotide sequence level, in agreement with previously published results (7). This has held true for the 212 isolates we have examined for sequence variation in gB thus far. Comparison of the Dumas (E) and Oka parental (J) strains identified three prominent SNPs in ORF31 (Table 1): one was at codon 73 (position 57224), resulting in an amino acid substitution (T versus P), and one silent mutation each was at codons 99 (position 57301) and 131 (position 57397). A curious finding was that all ORF22-classified E and J strains were allocated to the same genotypic clusters by analysis of the ORF31 5'-terminal region. Two of the M genotypic clusters could be resolved with the ORF22 fragment SNPs, specifically, a cluster with more J type variations (M1; African Morocco_01 and DRA_1 strains) or more E type variations (M2; USOR_12 and USOR_17) identified through ORF31 SNPs. Analysis using ORF62 variations yielded the same results (Fig. 1).
We also determined that the J and E gB alleles were each significantly linked with specific allelic forms of ORF22. All these data provide additional confirmation for the robustness of the J, E, and M genotyping. The extent to which these observed differences in selected fragments of the VZV genome are stabilized by environmental or host-related selective pressures remains to be investigated.
Toward a common strategy for VZV genotyping. The limited number of VZV SNPs are disseminated evenly over the genome but can nonetheless be used to classify strains into genotypes having a specific regional distribution. Evaluation of phylogenetic trees produced high bootstrap values, confirming the integrity of the genotypic analyses used in several reports (14, 33, 34). However, the use of this approach can be misleading, particularly for the identification of subgenotypes within the major circulating types, and must be confirmed by evaluating large numbers of clinically isolated specimens obtained from diverse geographic locales.
Although other investigators have attempted to describe the phylogenetic origin of currently circulating VZV strains, no consensus for a genotyping system has been reached. Here, we propose a simple methodology for elucidating the phylogenetic relationship between the VZV strains responsible for recent outbreaks and their evolutionary and/or transmission history. To prevent unnecessary confusion arising from the differing genotyping strategies and nomenclatures adopted by different laboratories, it is useful to summarize and correlate those findings here.
Two broad, nonintersecting studies focused on multiple SNPs positioned in various loci, comprehensively representing several thousand base pairs of sequence. Recent publications on ORF62 (tegument, transactivator), ORF31 (gB), ORF67 (gI), ORF68 (gE), ORF60 (gL), and ORF37 (gH) indicate that all are well conserved, and no single small region of the genome displays a genetic signature common to all VZV isolates in the United States (20, 33, 34).
Faga et al. and Wagenaar et al. (14, 44) proposed the division of VZV strains into four groups: (i) an Asian cluster (chiefly Japanese), (ii) a U.S. west coast genotype, (iii) a Dumas-like European cluster, and (iv) a central U.S. cluster. This represented the only effort previous to ours to characterize U.S. VZV genotypes but was limited by its examination of a relatively small number of strains and by a heavy reliance on laboratory-cultured strains.
Muir et al. and Quinlivan et al. (33, 34) examined a large number of VZV isolates from the United Kingdom, Asia, and Africa and, using sequence data obtained from ORF1, -21, -38, -50, and -54, identified three major genotypes, designated A, B, and C. B and C isolates were most commonly isolated in the United Kingdom and the United States, and A isolates were most common in Asia and Africa. However, the data in Fig. 1 indicate that all four of the targets in ORF1 selected for the Breuer study (33) coincidentally carry mainly J genotype-like SNPs at positions common to both J and M strains and as such could not clearly differentiate African strains from Asian.
Virtually all of our specimens were clinical isolates or low-passage strains, with only limited resort to established laboratory strains. It was observed (14) that laboratory strains of VZV can establish selected sequence variations in culture. For example, the T-to-C change at position 106262 that we have used to discriminate the vaccine strain from the wild type (30, 31) has been found only in the extensively cultured VZV Ellen and Oka vaccine strains (14). European genotype strains were most commonly isolated in Europe and in countries with a history of European colonization. The Japanese genotype was almost exclusively found in Japan, with the exceptions of Hawaii, the west coast of North America, and eastern Australia. Isolates from Africa and South America, as well as western Australia were of the M genotype.
The SNPs in the 447-bp fragment were able to discriminate several subtypes of M in a geographically distinctive fashion. These results refine the current understanding of VZV sequence variation and should permit both practical strain surveillance and the identification of the probable geographic origins of clinically isolated strains.
Speculation on mechanisms for establishing the M genotypes. VZV genomes characterized by alternating J-like and E-like regions detected in M strains can be explained by at least two mechanisms, neither of which is excluded or favored by the present state of knowledge. One mechanism is strain recombination in dually infected cells, and the other mechanism involves independently arising, nonrandom point mutations.
Mixed infections with two different strains have been reported for other human herpesviruses, notably human cytomegalovirus (7), Epstein-Barr virus (15), and human herpesvirus type 6 (8). Second-episode varicella has been observed in recent years, increasing the possibility of intrastrain recombination in the host (3, 18, 21). However, there is as yet no direct evidence that VZV recombination can occur in coinfected cells. A proposed recombinant genotype from Brazil described by Muir et al. (33) is likely an M genotype strain. Of 10 independently obtained isolates from Brazil tested in our laboratory, 8 strains were M (Fig. 4) and had SNP patterns similar to those described by Muir et al. (33). Since the VZV genome is so stable and observed DNA sequence variation is virtually all individual point mutations, it just as plausible that M strains reflect a small number of favored variations that have occurred independently among wild-type strains in more than one geographic location. We also cannot rule out the possibility that M virus occurred first in tropical regions and mutated following the migration of humans into temperate climates. J and E strains may have arisen independently from M as a result of the isolation of the Japanese islands and the broad explorations of the Europeans.
In conclusion, this study demonstrates that VZV genotypes can be classified into subgenotypes on the basis of limited DNA sequence variation in ORF22. Most important, this study introduces a new approach to the molecular epidemiology of VZV infection that should be especially advantageous for the study of the host-VZV interaction. The amplification of ORF22 sequences directly from infected tissues produces a broader image of VZV heterogeneity than do other methods described thus far. The application of this tool to the analysis of large collections of worldwide isolates should help to elucidate the mechanism or mechanisms through which VZV variability develops. Last, the method should provide a useful tool for the detection of VZV coinfection and could provide new insights into the mechanisms of VZV persistence.
We thank colleagues who provided samples for this study, including Neli N. Maltceva, Richard Garceau, Inger Damon, Michaela Schmidt, Herman Meyer, Mikhail Sousloparov, Karin Galil, Naoki Inoue, Michiko Takayama, Ann Arvin, Leigh Zerboni, John Zaia, and Sharka Nemechkova.
|
|
|---|
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»