Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2006, p. 2541-2546, Vol. 44, No. 7
0095-1137/06/$08.00+0 doi:10.1128/JCM.00054-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Epidemiology,1 Department of Pathobiology, School of Public Health and Community Medicine,4 Department of Medicine,3 Department of Laboratory Medicine, School of Medicine, University of Washington,5 Program in Infectious Diseases, Fred Hutchinson Cancer Research Center, Seattle, Washington2
Received 10 January 2006/ Returned for modification 15 March 2006/ Accepted 1 May 2006
|
|
|---|
|
|
|---|
To develop sequence-based methods for distinguishing between HSV-2 strains, we performed RFLP, HMA, and DNA sequence analysis of selected regions of the HSV-2 genome. We attempted to maximize the degree of sequence variability through our selection of genome regions for study. Noncoding genome regions were chosen for sequencing because we hypothesized a lack of purifying selection in these regions, and glycoprotein-coding regions were chosen based on previous sequence variation data in analogous segments in the HSV-1 genome (4, 13, 16). Our goals were to distinguish between strains from unrelated individuals and to identify epidemiologically linked HSV-2 strains.
(Some of the data in this work were presented at the International Herpesvirus Workshop, 25 to 30 July 2004 [abstr. 4.25].)
|
|
|---|
RFLP. Approximately 750 ng of HSV-2 DNA from each individual was digested to completion with SalI and separately with XhoI (Invitrogen, Carlsbad, CA) by methods similar to those reported previously (12). Digests were then electrophoresed on a 0.6% agarose gel overnight. DNA fragments were detected with SYBR green.
PCR. Seven segments of the HSV-2 genome (Table 1) were amplified by PCR from each isolate. Segments gC, noncoding region 1 (NC1), and NC3 were amplified using the Stratagene (La Jolla, CA) PfuUltra enzyme. Segments from gB, gG, NC2, and NC4 were amplified using the Applied Biosystems (Foster City, CA) Amplitaq enzyme when they failed to amplify efficiently with the PfuUltra enzyme. PCRs included 8% glycerol to increase reaction efficiency, 0.3 µM of each primer (Table 1), and 104 to 108 copies of DNA template based on standard quantitative real-time TaqMan PCR (Applied Biosystems, Foster City, CA) (18). Because of high GC content in the regions amplified, 5% dimethyl sulfoxide was added to the PCRs for all regions except NC1. SYBR green in dimethyl sulfoxide and ROX (Invitrogen, Carlsbad, CA) were added to enable real-time detection of amplicons. Three-step PCR cycling conditions included denaturing at 94°C for 20 s, annealing at 68°C for 30 s for NC1 and NC3 and at 60°C for 30 s for gG, gB, gC, NC2, and NC4, and extension at 72°C for 1 min for 35 cycles. Amplicons were then purified using DNA Clean & Concentrator (Zymo Research, Orange, CA).
|
View this table: [in a new window] |
TABLE 1. Regions of HSV-2 genome amplified by PCR
|
DNA sequencing. All PCR amplicons were routinely sequenced bidirectionally, except for that in region NC3, using BigDye Terminator (Applied Biosystems, Foster City, CA). Second-strand sequencing of NC3 was done only if polymorphisms were suspected on unidirectional sequencing. Sequencing primers were identical to those used for PCR (Table 1). The sequences were analyzed using Sequencher, version 4.1.4 (Gene Codes, Ann Arbor, MI). HSV-2 strain HG52 (GenBank accession no. NC001798) was used as the reference sequence. Visual genotypes were generated with polymorphism data using VG2 software (20). If the sequences from related pairs did not match, a second confirmatory PCR and sequence for the questionable individual and region were obtained using the original DNA product from culture. No unmatched polymorphisms were found to be caused by MDA or PCR artifacts through the study.
To assess the stability of the genomic HSV-2 regions during limited passage in culture, we compared DNA prepared as published elsewhere (18) directly from swabs of genital lesions to DNA prepared from HSV-2 cultured from the same lesions. The cultures were initially isolated on mink lung cells and were passed once on mink lung cells and then once in Vero cells prior to DNA isolation (QIAGEN [Valencia, CA] blood kit). The samples used for this stability analysis were isolated on mink lung cells instead of HDF cells because mink lung cells are currently used in our laboratory for this procedure, whereas HDF cells were used for isolation before our era of PCR specimen analysis. DNA was amplified and bidirectional sequencing performed as described for the other specimens in this report. For this analysis of stability in culture, genital HSV-2 lesion DNA specimens from four additional individuals before and after culture were aligned, and the original chromatograms were examined at any candidate loci of discrepancy.
Nucleotide sequence accession numbers. Sequences from this study have been submitted to GenBank under accession no. AY827201 to AY827401 and DQ236133 to DQ236152. The similarity of these sequences to the HSV-2 genome has been confirmed by comparison to the GenBank database by using megaBLAST (December 2005).
|
|
|---|
RFLP. Nine couples had identical RFLP patterns after digestion with each enzyme (Table 2; Fig. 1). Minor banding differences were observed between partners in five couples with SalI and between partners in three couples with XhoI. Related individuals showed no more than a 2-band difference by use of either enzyme. However, comparisons of unrelated individuals showed that no particular number of band differences could be used as a threshold for determining relatedness with one enzyme. For example, between pairs 2, 4, 7, 8, 13, and 14, the mean number of band differences between unrelated individuals was 3.9, with some unrelated individuals differing by only 1 band (Fig. 1). However, each individual was more closely related to his or her partner than to any unrelated individual.
|
View this table: [in a new window] |
TABLE 2. Number of band differences within pairs, by restriction enzyme
|
![]() View larger version (131K): [in a new window] |
FIG. 1. Results of RFLP for six selected pairs. HSV-2 genomic DNA was digested using the SalI enzyme and electrophoresed on a 0.6% agarose gel overnight with a 1-kb molecular weight marker. RFLP pattern differences between partners can be seen with pairs 13 and 4.
|
|
View larger version (16K): [in a new window] |
FIG. 2. HMA results for HSV-2 region NC1 amplicons by pair and partner. A radiolabeled amplicon from HSV-2 strain 333 was hybridized to each sample and then electrophoresed on a 0.6% polyacrylamide gel overnight. Slight mobility shifts are visible for pair 2, pair 8, partner B of pair 1, and partner A of pair 3.
|
Noncoding regions showed the most variation. Among all 28 samples, one or more strains showed variation at a total of 22.44 nucleotides/1,000, 30.04 nucleotides/1,000, and 33.33 nucleotides/1,000 for NC1, NC3, and NC4, respectively. No variation was found in 190 bp of NC2. In coding regions, one single-nucleotide polymorphism each was found in 152 bp of gG and 228 bp of gC, and no variation was present in 217 bp of gB. Overall, 21 sites of sequence variation were found in 1,482 bp of analyzed sequence (Table 3). Five of those sites consisted of variable-length repeats of one or two nucleotides. Two of these repeat mutation sites were in NC1, two were in NC2, and one was in NC4 (Table 3). The remaining sites were single-nucleotide polymorphisms (Fig. 3). In addition to the variations listed, one individual had the sequence GGAGAGGGAGAGGGA in place of a variable-length G repeat located near position 84794 in the NC4 region.
|
View this table: [in a new window] |
TABLE 3. HSV-2 sequence polymorphisms
|
![]() View larger version (68K): [in a new window] |
FIG. 3. Individual genotypes at each polymorphism. The consensus alleles (blue) and minor alleles (yellow) determined by the population studied are listed for each individual at each locus where a sequence polymorphism was observed by DNA sequence analysis. Missing data are indicated by gray squares or the letter N. The number of C repeats at locus 48802 and the number of GT repeats at locus 84796 are listed for each individual. Asterisks indicate regions where no sequence variation was found. (A) Related pairs with matching sequences; (B) related pairs with differing sequences.
|
Taking into account all sites of variation, it was possible to differentiate between all unrelated infections (Fig. 3). All sequences were also distinguishable from the published HG52 sequence. Eleven of the 14 pairs had sequences that were identical between the related partners at all sites. Pair 8 differed at site 48958, a variable-length C repeat in a noncoding region (Table 3). Pair 12 differed at site 84796, a variable-length GT repeat in a noncoding region. Pair 10 differed at six sites. A clinical chart review revealed multiple concurrent sexual partners for partner B from pair 10. This sexual history, in combination with the sequence data, indicated that partner A is an unlikely source of partner B's genital herpes infection.
To assess the stability of the sequences during tissue culture, we evaluated DNA samples from HSV-2 lesions from four individuals before and after culture. All genomic regions were evaluated for specimens from all but one individual. We found no differences in sequence between DNA amplified directly from lesions and DNA obtained after virus was cultured from the same clinical lesions and passaged twice. During this sequencing of additional strains, occasional 3-nucleotide polymorphisms were uncovered, two in the noncoding region NC1 and one in the noncoding region NC4. These can be accessed from the GenBank submissions accompanying this report. These data are consistent with stability in HSV-2 genomic DNA during a limited number of in vitro passages in these loci and the identification of bona fide polymorphisms in our study.
|
|
|---|
Past HSV sequencing efforts have focused on coding regions, particularly glycoprotein sequences (4, 13, 16). Two publications have reported that large sections of HSV-1 sequence, including complete coding regions, could be used for examining strain differences between populations (4) and individuals (16). However, a small segment (179 bp) of the HSV-2 DNA polymerase gene was found to have no variability between two samples and a standard strain (13). Nevertheless, we found that similarly small segments could be of use if they were generated from intergenic spacers that do not code for a product.
Five of the 21 polymorphisms observed were variable-length repeats of one or two bases. Similar variations in the length of GT or C repeats in HSV sequences due to DNA polymerase error have been reported elsewhere (9, 10). These polymerase errors are possible at several points in the pathogenesis of HSV-2. Mutations may have occurred during viral replication in the recurrent lesion of partner A or the primary lesion of partner B. Alternatively, a variant strain may have colonized the ganglia of partner A in the intervening time between sample collection from partner A and partner B, which was 31 days for pair 8 and 51 days for pair 12.
We find it unlikely that HSV-2 sequence variation was introduced by the methods used in our study. First, even though mutations in DNA sequences can be introduced by the polymerases used in PCR and sequencing, all polymorphisms were confirmed by second-strand sequencing, and nucleotide differences between partners were confirmed by a second PCR and sequencing using the DNA harvested from culture. Second, we found the sequences to be stable when DNA samples taken directly from the lesion swab were compared to samples obtained after 2 passages in cell culture. There are no data to suggest that genetic drift in 2 passages in mink lung cells would be any more or less likely than drift in passages in HDF cells. Sequence drift during culture has been reported in noncoding hypervariable regions of HSV-1 (14); however, we found no such drift in our study.
We expected to observe increased variability between individuals in noncoding regions due to our hypothesized absence of purifying selection in these regions. Our study confirmed that noncoding regions exhibit a greater level of variability than similar-length segments of glycoprotein-coding regions. It is possible that, during HSV-2 replication, DNA polymerase errors in noncoding regions may not decrease the fitness of the virus. This is in contrast to the coding regions, which evolve with functional restraints and therefore are expected to be more conserved. The increased variability of noncoding regions facilitates the detection of variation with a much smaller, and thus more easily obtainable, amplicon. Overall, noncoding regions may be preferable for distinguishing HSV-2 strains at an individual level.
While we believe our major conclusions are firm, we do recognize that we sampled a small portion of the HSV-2 genome. While our source partners did not acquire HSV-2 exclusively in Seattle, we did analyze only a limited number of persons. It is possible that sequencing of HSV strains from a larger geographical area may show greater variability. However, it is noteworthy that HG52 differed from the consensus sequence only at three sites, even though its origin is geographically distant from that of our samples (7). In addition, the amplified segments of gC, gB, and gG were chosen by ease of primer selection and amplification. These factors may have influenced the selection of a more conserved region. Consequently, the variation we observed in those segments may not be indicative of the level of variation present across each entire coding region, and further research is needed to determine to what extent these polymorphisms, as well as others, can be used to define and compare HSV-2 strains across individuals and populations. Specific attention must be given to sequencing efforts that analyze large numbers of individuals to refine the use of HSV-2 sequencing for molecular epidemiology. Furthermore, DNA sequencing of consecutive isolates from an individual can confirm whether an HSV-2 infection in a single individual is limited to a single strain or whether more than one strain has colonized the ganglia. While existing literature suggests that most infections occur with a single strain, infection with more than a single strain has been described using RFLP (5, 19, 21).
In conclusion, we found that DNA sequencing was more powerful than both RFLP and HMA in distinguishing HSV-2 strains. The sequences obtained from noncoding regions provided the most information on differences between unrelated samples. Our results suggest that a multilocus assay based on several DNA sequences has the potential to be an informative tool for identifying specific strains in an HSV-2 transmission event.
We acknowledge Lizette Embuscado for her technical assistance with the RFLP analysis.
|
|
|---|
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»