Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2004, p. 2048-2053, Vol. 42, No. 5
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.5.2048-2053.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Pediatrics, Sapporo Medical University School of Medicine, Sapporo, Hokkaido 060-8543, Japan
Received 23 November 2003/ Returned for modification 31 December 2003/ Accepted 3 February 2004
|
|
|---|
|
|
|---|
RSV was initially found to have two distinct antigenic groups, designated A and B, by their different reactivity with monoclonal antibodies against the viral antigens (1, 20). Epidemiology studies of RSV demonstrated that both antigenic groups circulate concurrently or alternately in a community during epidemics (18, 32). Among viral surface antigens, the attachment G glycoprotein (G protein) has the greatest antigenic diversity between these two groups and among strains within each group. The G-protein variability is concentrated in its ectodomain, which contains two hypervariable regions separated by a conserved 11-amino-acid motif (4, 7, 19, 29). Recent molecular epidemiological studies suggest that many distinct RSV genotypes circulate worldwide and that similar genetic variants are clustered by time rather than by geographic location (6, 8, 17, 26).
Genetic variability of RSV has most often been studied by using restriction fragment length polymorphism (RFLP) analysis and DNA sequencing (4, 30). RFLP analysis is relatively easy to perform, although genotypes obtained by this method provide limited information. DNA sequencing, in contrast, shows the most precise genetic variability, although this method is not available to all laboratories because of costs and equipment.
Heteroduplex mobility assay (HMA) was originally elaborated to detect genetic diversity of the human immunodeficiency virus type 1 in 1994 (10, 11). Since then, it has been widely applied to other viruses for subtyping a viral species (2, 14) and for molecular epidemiology studies (13) and to detect viral genome quasispecies in a single patient (12). The advantages of HMA for molecular epidemiology viral studies are that it is technically simple, it costs less, and it is less time-consuming than the other methods mentioned above. In our study, HMA provided more details of RSV genetic diversity than did RFLP and was compatible with the results of DNA sequencing.
|
|
|---|
12 months old: 51 (47.1%), 34 (31.7%), 19 (17.3%), and 4 (3.9%) patients were diagnosed as having bronchiolitis, tracheobronchitis, pneumonia, and upper respiratory tract infection, respectively. RNA extraction and cDNA synthesis. RNA extraction and cDNA synthesis were described previously (26). Briefly, the frozen RSV isolates were inoculated on HEp-2 cells in 24-well plates and total RNA was extracted by adding 0.8 ml of RNAzol B (TEL-TEST, Inc., Friendswood, Tex.) to each well when the extensive cytopathic effect was observed. For cDNA synthesis, 100 ng of the total RNA was reverse transcribed with Moloney murine leukemia virus reverse transcriptase (Invitrogen) and either random hexamer (Takara) or oligo(dT)12-18 primer (Invitrogen) at 37°C for 1 h.
PCR. PCR was performed as described previously with minor modifications (26). Primers used in this PCR amplified the sequence between nucleotides (nt) 1 and 584 of the G gene of the Long strain, which covers the first hypervariable region of the G gene. The sequence of each primer was as follows: primer 1 (5'-GGGGCAAATGCAAACATGTCC-3'; nt 1 to 21) and primer 2 (5'-GGTATTCTTTTGCAGATAGC-3'; nt 565 to 584) (6, 35). PCR was carried out in a total volume of 100 µl containing a 1 µM concentration of each primer, a 250 µM deoxynucleoside triphosphate concentration (Takara), 10 mM Tris-HCl (pH 9.0), 500 mM KCl, 1% Triton X-100, 1.5 mM MgCl2, 2 U of Taq DNA polymerase (Promega, Madison, Wis.), and 2 to 4 µl of cDNA template. The mixture was heated at 94°C for 2 min, followed by 30 cycles of denaturing at 94°C for 45 s, annealing at 54 to 57°C for 45 s, and extension at 74°C for 45 s, with a final extension at 74°C for 7 min by using GeneAmp PCR System 9700 (Applied Biosystems, Foster City, Calif.).
RFLP analysis. The procedure of RFLP was described previously (26). Briefly, 8 µl of each PCR product was digested with AluI, TaqI, and NdeII at 37°C for 90 min. The digested PCR products were analyzed on 2% agarose gels stained with ethidium bromide.
HMA. To form heteroduplexes, the Sap/645/80 strain of RSV which was the earliest isolate from the study community was used throughout as a reference. An equal volume (4.5 µl) of reference PCR product was mixed with a sample PCR product and 1 µl of 10x annealing buffer (1 M NaCl, 100 mM Tris-HCl [pH 7.8], 20 mM EDTA) (11). The fragment mixture was denatured at 94°C for 2 min and annealed by rapid cooling at 4°C in the thermal cycler. After it cooled, the mixture was loaded onto a 5% polyacrylamide gel. Electrophoresis was performed in a vertical gel apparatus (ATTO, Tokyo, Japan) at 150 V for 3 h. DNA fragments were visualized by ethidium bromide staining.
Nucleotide sequencing and computer analysis. The PCR products were sequenced directly by a BigDye Terminator v3.0 cycle sequence kit and an ABI Prism 310 Genetic Analyzer (Applied Biosystems). Phylogenetic analysis was done by the neighbor-joining method with CLUSTAL W in DDBJ, and a phylogenetic tree was made by TreeView (version 1.6.6).
|
|
|---|
![]() View larger version (53K): [in a new window] |
FIG. 1. Restriction fragment analysis of amplified PCR products derived from cDNAs of the G protein gene between nt 1 and 584. The PCR products were digested with AluI, TaqI, and NdeII and separated by 2% agarose gel electrophoresis. Five, three, and four restriction patterns were obtained with AluI (A to E), TaqI (A to C), and NdeII (A to D), respectively. The molecular weight marker (100-bp DNA ladder) is shown on the left.
|
|
View this table: [in a new window] |
TABLE 1. Distribution of RFLP genotypes of RSV subgroup A isolates during 19 epidemics
|
|
View this table: [in a new window] |
TABLE 2. Distribution of HMA genotypes of RSV subgroup A isolates during 19 epidemics
|
|
View this table: [in a new window] |
TABLE 3. Comparison of 9 RFLP and 31 HMA genotypes of RSV subgroup A isolates during 19 epidemics
|
Genotype distribution determined by RFLP analysis and HMA. One to four different RFLP genotypes emerged in each epidemic season (Table 1). Six genotypes appeared during three consecutive seasons, and then the genotype changed to another. Genotypes AAA, ABA, and BAA reemerged several years after the first dominance disappeared. Genotype AAA, which is the earliest genotype in the present study, was observed intermittently over a 20-year-period. Likewise, genotypes ABA and BAA also circulated intermittently for 13 and 11 years, respectively. The other six genotypes were detected in only one to four consecutive seasons, and these genotypes never reemerged.
One to six different HMA types occurred in each epidemic season, as was seen in the RFLP genotype distribution (Table 2). However, one HMA genotype lasted only one to three epidemic seasons and was then replaced by another genotype. The maximum duration of circulation of one genotype was 5 years, and individual genotypes usually did not reemerge after the initial circulation period in this observation, although genotype d' surfaced again exceptionally after a 2-year break.
Correlation between HMA and nucleotide sequencing. We looked for any correlation between the mobility of heteroduplexes and the degree of nucleotide differences in the corresponding sequence. We selected 21 isolates belonging to different HMA types during the 15 epidemic seasons from 1980 to 1995. Their nucleotide sequences (nt 1 to 584) were compared to that of the reference strain (Sap/645/80) which belongs to genotype a. These strains were aligned in a polyacrylamide gel in order of faster (left) to slower (right) migration (Fig. 2). The nucleotide changes of these isolates compared to those of the reference strain were all within the range of 43 to 51 nt. Therefore, in the present study, a correlation between the mobility of heteroduplexes and the number of point mutations could not be fully evaluated. However, it appeared that strains with a similar amount of G-gene nucleotide differences from the reference strain were often distinct HMA types. For example, the strains belonging to distinct HMA genotypes k, b, o, n, and f had 45- nt changes compared to a reference strain.
![]() View larger version (33K): [in a new window] |
FIG. 2. Heteroduplex mobility assay with amplicons of cDNAs of the G gene (nt 1 to 584) of 21 isolates during 15 epidemic seasons from 1980 to 1995. Heteroduplexes with 21 distinct mobility patterns were aligned from faster (left) and slower (right) migrations, and alphabetical assignment was provided as shown in Table 2. The number of nucleotide changes compared to those of the reference strain are indicated below. The molecular weight marker (100-bp DNA ladder) is shown on both sides.
|
RFLP, HMA, and phylogenetic analysis. We selected one isolate from each of 31 HMA genotypes and performed phylogenetic analysis using whole sequences from nt 1 to 584. HMA and RFLP genotypes of each strains were also indicated. The phylogenetic tree showed seven clusters (clusters 1 to 7), and each cluster consisted of two to nine isolates that have close isolation times (Fig. 3). Genetic heterogeneity within RSV field strains detected by HMA analysis was confirmed phylogenetically. That is, all 31 strains with individual HMA genotypes were discriminated from each other by nucleotide sequence analysis. Strains within the same branches usually exhibited similar but distinct HMA types (types f, l, and g in C-1; types m and i in C-2; types p, s, and n in C-3; and types x and z in C-4), although there were some exceptions, i.e., types k and j in C-1 and types q and v in C-4. On the other hand, strains belonging to distant clusters sometimes presented similar, although distinct HMA patterns, e.g., types e, h, and p in C-1, C-7, and C-3, respectively (Fig. 2 and 3).
![]() View larger version (25K): [in a new window] |
FIG. 3. Phylogenetic tree of 31 isolates across all HMA genotypes based on sequences of the first variable region of the G gene of RSV. Names of HMA genotypes (a to e') and RFLP genotypes (three capital letters) are indicated next to the strain names. The phylogenetic clusters (C-1 to C-7) are also indicated.
|
|
|
|---|
Trincado et al. (31) used HMA for genotyping human cytomegalovirus isolates and compared this method to RFLP, single-stranded conformation polymorphism, and DNA sequencing analyses. These authors concluded that HMA was the most accurate and concise of the four methods used to genotype human cytomegalovirus isolates.
Genetic variability of the G gene of RSV has been investigated during the past 10 years in both developed and developing countries by using RFLP and nucleotide sequencing (5, 8, 9, 22-25, 30, 34). These investigations revealed that multiple genotypes cocirculated in each epidemic and that the dominant genotype changes year by year. The present study demonstrates that the results of HMA for the G-gene variability correlate well with those of RFLP analysis and, moreover, that HMA can detect more detailed genotype manifestations. We and others (8) found that some RFLP genotypes emerged for very long periods. In a study of the molecular epidemiology of RSV in Korea, one genotype (GP-A2) of subgroup A was detected over eight consecutive epidemics, and one subgroup B genotype (GP-B1) was also found for six epidemics (8). In our study, the genotype AAA that was epidemic in the 1980 and 1981 season emerged repeatedly over 20 years. However, the genotype AAA strains isolated after these two seasons were never assigned to the same HMA genotype as that of the former AAA strains. The same observation was found in RFLP genotype ABA, which had appeared over 13 years. The genotype ABA corresponded to seven HMA genotypes and two of them, genotypes d and f, appeared for only 3 and 4 years, respectively. All strains with different HMA genotypes were confirmed as discrete strains by nucleotide sequence analysis. Our findings show that HMA is precise and can detect more subtle genetic changes in RSV isolates than can RFLP. In this regard, HMA is more discriminatory than RFLP. Our observations also indicate that in RSV epidemics new genotypes are continuously emerging, and these new genotypes prevail for only a limited period, at most 4 to 5 years.
We confirmed that RSV field strains with similar sequence divergence from the reference strain could often display distinct HMA patterns. However, it also appeared that the locations of nucleotide changes in each strain compared to the reference strain varied, although the total number of mutations was the same or very similar between the strains. Thus, our results suggest that the mobility of heteroduplexes is not only based on the number of point mutations but also correlates with the location of the mismatch, i.e., structural differences in heteroduplexes reduce the electrophoretic mobility. Therefore, HMA may be more sensitive with respect to structural deformations in double-stranded DNA than nucleotide sequence. The nucleotide sequences of isolates in the same HMA genotype had at least 98% homology. This was comparable to the previous observation that a sequence divergence of more than 1 to 2% was necessary for detectable heteroduplex retardation (3, 10).
The HMA had some correlation with phylogenetic analysis of the G gene by nucleotide sequencing. In the present study, the phylogenetic tree showed seven clusters, and each cluster consisted of two to nine HMA genotypes. The strains within each branch represent similar but distinct HMA types. However, the degree of genetic divergence evident in the phylogenetic analysis as distinct clusters could not be represented precisely by HMA analysis. Therefore, the HMA analysis in the present study could distinguish genetic divergence among strains with a higher sensitivity, but it might not be possible to express these quantitatively.
In a study of the differentiation of influenza A viruses from different animal species, such as avian and swine viruses, there was an excellent correlation between HMA results and nucleotide divergence, as determined by direct sequencing and phylogenetic analysis (14). In our study, the RSV field strains with various sequence divergence from reference DNA could not be evaluated; however, we confirmed that HMA could differentiate between strains with a similar sequence divergence from reference DNA.
In conclusion, it appears that HMA is useful for elucidating the genetic variability of RSV field strains. This method can demonstrate a clear distribution of many RSV genotypes and provides more precise molecular epidemiology than RFLP analysis. In addition, the HMA genotypes correlated to some extent with the clustering of strains in a phylogenetic tree.
HMA is technically easy and inexpensive compared to other conventional methods for the assessment of genetic diversity of viral strains. In particular, in developing countries, the molecular epidemiology of RSV epidemics is becoming more important, but equipment for DNA sequencing is not always available in many laboratories. For this reason, HMA is an alternative and important method for identifying genetic variants and investigating the molecular epidemiology of RSV infection. Further prospective studies are needed to support and confirm these observations.
|
|
|---|
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»