Enterovirus Laboratory, Department of Microbiology, National Public Health Institute (KTL), Mannerheimintie 166, FIN-00300 Helsinki, Finland1
Author for correspondence: Tapani Hovi. Fax +358 9 4744 8355. e-mail tapani.hovi{at}ktl.fi
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
HRV serotypes can be grouped in different ways, on the basis of shared structural and biological properties, e.g. receptor specificity (Abraham & Colonno, 1984 ; Uncapher et al., 1991
), sensitivity to antiviral agents (Andries et al., 1990
), immune cross-reactivity (Cooney et al., 1982
) or genetic similarity (Duechler et al., 1987
; Horsnell et al., 1995
). Several groups have published studies on improved methodology for the detection of HRV infection, including molecular amplification methods like PCR and NASBA (Gama et al., 1989
; Hyypiä et al., 1989
; Arruda & Hayden, 1993
; Ireland et al., 1993
; Johnston et al., 1993
; Halonen et al., 1995
; Samuelson et al., 1998
; Andeweg et al., 1999
, Blomqvist et al., 1999
). These studies, however, were designed to recognize HRV as a group, and not the individual serotypes. Molecular methods, e.g. partial genome sequencing, have been used in rapid genetic typing of the closely related group of enteroviruses (Oberste et al., 1999a
, b
). In addition, molecular analysis of wild-type poliovirus strains is routinely used in the Global Poliovirus Eradication Initiative in the follow-up of the programme (Kew et al., 1995
). A rapid molecular identification method might provide a possibility for investigation of potential differences in pathogenicity between the different HRV serotypes and strains.
We have recently shown that HRV strains isolated from a cohort of small children during two successive years show a striking genetic divergence but segregate in the two known genetic clades of HRV (Savolainen et al., 2001 ). Our aim of this study was to explore the possibility of genetic typing of HRV strains. For this purpose we have partially sequenced all of the HRV prototype strains.
![]() |
Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
RNA isolation and RTPCR.
The RNA was isolated from 100 µl cell culture homogenate with the RNeasy Total RNA kit (Qiagen). Upon purification, RNA was eluted from the columns with 30 µl diethyl pyrocarbonate-treated H2O and subsequently stored at -70 °C. RTPCR was performed using primer pair 9565-reverse and 9895-forward [9565 position 10831058 according to HRV1b by Hughes et al. (1988) , 5' GCA TCI GGY ARY TTC CAC CAC CAN CC 3'; 9895 position 534560 in HRV1b, 5' GGG ACC AAC TAC TTT GGG TGT CCG TGT 3'; I=inosine; Y=T, C; R=G, A; N=A, G, C, T], producing a fragment of 549 nt spanning the hypervariable part of the 5'NCR, the entire VP4 gene and the 5' terminus of the VP2 gene. cDNA synthesis and PCR were carried out as described previously (Mulders et al., 2000
). The PCR products were visualized after electrophoresis on an ethidium bromide-stained 2% agarose gel.
Sequence analysis.
PCR products were purified with a PCR purification kit (QIAquick, Qiagen) in the case of a single band in the electrophoresis, or in the case of multiple bands using a gel extraction kit (QIAquick, Qiagen). The purified products were stored at -20 °C or directly used in the cycle sequencing (ABI Prism BigDye Terminator Cycle Sequencing Ready Reaction kit; Applied Biosystems) with the same forward and reverse primers as used in the RTPCR. An automated DNA sequencer was used for sequencing (ABI PRISM, model 377). Sequence data were analysed with Sequencing Analysis (version 3.1, ABI) and Sequence Navigator (version 1, ABI) for pairwise comparisons. Multiple sequence alignments were made using PILEUP, part of the GCG program suite (version 10, Genetics Computer Group, Madison, WI, USA), and CLUSTALX. The analysed fragment was 420 nt in the VP4/VP2 region (207 nt in VP4 and 213 nt in VP2). Distance matrices were estimated using the DNADIST and PROTDIST programs, part of the PHYLIP (Phylogeny Inference) package (version 3.572c; Felsenstein, 1993 ), using the maximum likelihood model of nucleotide substitution. Dendrograms were drawn using the neighbour-joining option in NEIGHBOR (PHYLIP), and were visualized using NJPLOT or TREE VIEW (version 1.5.3). Bootstrap analysis was performed using SEQBOOT (PHYLIP) with 100 or 1000 replicates. For prototype strains HRV31, HRV32, HRV40, HRV66, HRVHanks2102 and HRV79 the analysed sequence was shorter than the standard 420 nt. However, practically identical trees were seen after replacing the terminal deletion with a series of N, or after editing the complete sequences to the length of the shortest one before analysis.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
|
|
Clustering of field strains with the prototype strains
Using aligned sequences in the VP4/VP2 region, most of the clinical isolates could be relatively easily linked to a single prototype strain. In all but one out of the 24 field strains in group A, a prototype strain was closer than 10% (Fig. 2), while four out of the 37 group B field strains remained outside this threshold difference. These strains isolated in September 1996 formed a distinct cluster of their own and did not definitely associate with any of the prototype strains (Fig. 3
). Bootstrap values for demarcation between the recent field isolates on one hand and the closest prototype strain on the other were usually very high, suggesting evolutionary difference. Judging from the observed interserotypic difference between the prototype strains, the 61 field strains could represent as many as 19 different serotypes.
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Our current results are in full agreement with previously published clustering of 21 strains in VP2 (Horsnell et al., 1995 ). The ratio of the prototype strains segregating into the two HRV species was 3:1, while the distribution of group A and B strains among the 61 recent Finnish field isolates was different, with a clear majority in group B. As noted before (Savolainen et al., 2001), the intra-epidemic diversity of HRV strains was remarkable, with several different serotype-like clusters circulating concomitantly in a cohort of about 300 children younger than 2 years of age living in a given suburb.
The clustering of HRV87 with HEV-D was a striking observation but in line with the report of Andeweg et al. (1999) who showed that HRV87 sequences in the 5'NCR cluster separately from those of the remaining HRVs. The identity of HRV87 was confirmed by specific neutralization with a valid monotypic antiserum obtained from the ATCC, and the typical rhinovirus character of acid-lability was also demonstrated by a standard test (Couch, 1992
). HRV87 was definitely acid-labile, typical of rhinoviruses, while EV70 was acid-stable, typical of enteroviruses (data not shown). It has been known for a long time that HRV87 differs from the other HRVs in receptor specificity (Uncapher et al., 1991
). The receptor is sialidase-sensitive, which has also been reported for EV70 (Utagawa et al., 1982
). These observations have prompted further investigations on the relatedness of HRV87 to human enteroviruses (to be reported separately), and may pose taxonomic problems concerning criteria for differentiation between HRV and enteroviruses.
The patterns of antiviral drug sensitivity of HRV strains were previously suggested to follow the genetic grouping (Andries et al., 1990 ), a view partially supported by our current findings. Andries and coworkers divided HRV serotypes, after a very complex multivariate analysis, into antiviral sensitivity clusters A and B. All sequenced genetic group B strains belong to antiviral sensitivity cluster A, which also includes a few definite HRV-A serotypes, HRV8, HRV13, HRV18, HRV32, HRV43, HRV54 and HRV95, close to the intergroup border. According to receptor usage, HRV serotypes have been grouped in the major and minor receptor groups, respectively, and HRV87, which uses yet another molecule for cellular interactions (Uncapher et al., 1991
). Both genetic groups were already previously known to include members of the major receptor group. All the ten minor receptor group members were now partially sequenced and found to belong to genetic group A. At the VP4/VP2 region they did not cluster more closely together (Fig. 2
). However, since VP1 is known to be the counterpart for receptor binding in the case of HRV2 (Hewat et al., 2000
), it would be intriguing to know if alignments made according to VP1 amino acid sequences would change the situation.
With the present data, further genetic subgrouping of HRV serotypes within either of the two established HRV species would appear difficult. Bootstrap values supporting the major branches were very small, especially in genetic group A (Fig. 2). However, differences between distinct serotypes were often, although not always, well supported. Most interserotypic differences in the VP4/VP2 region resemble those between enterovirus serotypes, but for a few pairs of HRV prototypes the difference was of the order of the intergenotypic differences within a given enterovirus serotype, or even less (Mulders et al., 2000
). The strains labelled HRV8 and -95 appear almost identical, and the proposed new serotype represented by strain Hanks (Andries et al., 1990
) clustered very close to the prototype strain of HRV21. Some similarity with previously reported antigenic cross-reactivity (Cooney et al., 1982
) could be found in the clustering of HRV prototype strains. For instance, the reported cross-reactions between HRV2 and -49 and between HRV29, -44 and -62, respectively, are not surprising regarding their relatively close genetic similarity.
Clustering of the clinical isolates with each other on one hand, and with given prototype strains on the other, suggests that genetic typing of HRV strains might also be possible according to the lines previously described for human enteroviruses (Oberste et al., 1999a , b
). For enteroviruses, sequence comparisons in the VP1 region give a more reliable serotype identification than those in the VP4/VP2 region, as the demarcation between interserotypic and intraserotypic variations is greater in VP1 (Oberste et al., 1999a
; Mulders et al., 2000). The current results are promising, but much larger numbers of serotype-confirmed independent field isolates have to be analysed before it is possible to evaluate development of a genomic sequence-based identification system for HRV isolates.
![]() |
Acknowledgments |
---|
![]() |
Footnotes |
---|
b Present address: Laboratoire National de Santé, Department of Immunology, Luxembourg.
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Abzug, M. J., Beam, A. C., Gyorkos, E. A. & Levin, M. J. (1990). Viral pneumonia in the first month of life. Pediatric Infectious Disease Journal 9, 881-885.[Medline]
Andeweg, A. C., Bestebroer, T. M., Huybreghs, M., Kimman, T. G. & de Jong, J. C. (1999). Improved detection of rhinoviruses in clinical samples by using a newly developed nested reverse transcriptionPCR assay. Journal of Clinical Microbiology 37, 524-530.
Andries, K., Dewindt, B., Snoekes, J., Wouters, L., Moereels, H., Lewi, P. J. & Janssen, P. A. J. (1990). Two groups of rhinoviruses revealed by a panel of antiviral compounds present sequence divergence and differential pathogenicity. Journal of Virology 64, 1117-1123.[Medline]
Arola, M., Ziegler, T., Ruuskanen, O., Mertsola, J., Nanto-Salonen, K. & Halonen, P. (1988). Rhinoviruses in acute otitis media. Journal of Pediatrics 113, 693-695.[Medline]
Arruda, E. & Hayden, F. G. (1993). Detection of human rhinovirus RNA in nasal washings by PCR. Molecular and Cellular Probes 7, 373-379.[Medline]
Blomqvist, S., Skyttä, A., Roivainen, M. & Hovi, T. (1999). Rapid detection of human rhinoviruses in nasopharyngeal aspirates by a microwell reverse transcriptionPCR-hybridization assay. Journal of Clinical Microbiology 37, 2813-2816.
Cooney, M. K., Fox, J. P. & Kenny, G. E. (1982). Antigenic groupings of 90 rhinovirus serotypes. Infection and Immunology 37, 642-647.[Medline]
Couch, R. B. (1992). Rhinoviruses. In Laboratory Diagnosis of Viral Infections , pp. 709-729. Edited by E. H. Lennette. New York:Marcel Dekker.
Duechler, M., Skern, T., Sommergruber, W., Neubauer, C., Gruendler, P., Fogy, I., Blaas, D. & Kuechler, E. (1987). Evolutionary relationships within the human rhinovirus genus: comparison of serotypes 89, 2 and 14. Proceedings of the National Academy of Sciences, USA 84, 2605-2609.[Abstract]
Felsenstein, J. (1993). PHYLIP (Phylogeny Inference package) version 3.5c. Distributed by the author, Department of Genetics, University of Washington, Seattle, USA.
Gama, R. E., Horsnell, P. R., Hughes, P. J., North, C., Bruce, C. B., al-Nakib, W. & Stanway, G. (1989). Amplification of rhinovirus specific nucleic acids from clinical samples using the polymerase chain reaction. Journal of Medical Virology 28, 73-77.[Medline]
Halonen, P., Rocha, E., Hierholzer, J., Holloway, B., Hyypiä, T., Hurskainen, P. & Pallansch, M. A. (1995). Detection of enteroviruses and rhinoviruses in clinical specimens by PCR and liquid-phase hybridization. Journal of Clinical Microbiology 33, 648-653.[Abstract]
Hewat, E. A., Neumann, E., Conway, J. F., Moser, R., Ronacher, B., Marlovits, T. C. & Blaas, D. (2000). The cellular receptor to human rhinovirus 2 binds around the 5-fold axis and not in the canyon: a structural view. EMBO Journal 19, 6317-6325.
Horsnell, C., Gama, R. E., Hughes, P. J. & Stanway, G. (1995). Molecular relationships between 21 human rhinovirus serotypes. Journal of General Virology 76, 2549-2555.[Abstract]
Hughes, P. J., North, C., Jellis, C. H., Minor, P. D. & Stanway, G. (1988). The nucleotide sequence of human rhinovirus 1B: molecular relationships within the rhinovirus genus. Journal of General Virology 69, 49-58.[Abstract]
Hyypiä, T., Auvinen, P. & Maaronen, M. (1989). Polymerase chain reaction for human picornaviruses. Journal of General Virology 70, 3261-3268.[Abstract]
Ireland, D. C., Kent, J. & Nicholson, K. G. (1993). Improved detection of rhinoviruses in nasal and throat swabs by seminested RTPCR. Journal of Medical Virology 40, 96-101.[Medline]
Johnston, S. L., Sanderson, G., Pattemore, P. K., Smith, S., Bardin, P. G., Bruce, C. B., Lambden, P. R., Tyrrell, D. A. J. & Holgate, S. T. (1993). Use of polymerase chain reaction for diagnosis of picornavirus infection in subjects with and without respiratory symptoms. Journal of Clinical Microbiology 31, 111-117.[Abstract]
Johnston, S. L., Pattemore, P. K., Sanderson, G., Smith, S., Lampe, F., Josephs, L., Symington, P., OToole, S., Myint, S. H., Tyrrell, D. A. J. & Holgate, S. T. (1995). Community study of role of virus infections in exacerbations of asthma in 911 year old children. British Medical Journal 310, 1225-1229.
Kew, O. M., Mulders, M. N., Lipskaya, G. Y., da Silva, E. E. & Pallansch, M. A. (1995). Molecular epidemiology of enteroviruses. Seminars in Virology 6, 401-414.
Mäkelä, M. J., Puhakka, T., Ruuskanen, O., Leinonen, M., Saikku, P., Kimpimäki, M., Blomqvist, S., Hyypiä, T. & Arstila, P. (1998). Viruses and bacteria in the etiology of the common cold. Journal of Clinical Microbiology 36, 539-542.
Mulders, M. N., Salminen, M., Kalkkinen, N. & Hovi, T. (2000). Molecular epidemiology of coxsackievirus B4 and disclosure of the correct VP1/2Apro cleavage site: evidence for high genomic diversity and long-term endemicity of distinct genotypes. Journal of General Virology 81, 803-812.
Oberste, M. S., Maher, K., Kilpatrick, D. R. & Pallansch, M. A. (1999a). Molecular evolution of the human enteroviruses: correlation of serotype with VP1 sequence and application to picornavirus classification. Journal of Virology 73, 1941-1948.
Oberste, M. S., Maher, K., Kilpatrick, D. R., Flemister, M. R., Brown, B. A. & Pallansch, M. A. (1999b). Typing of human enteroviruses by partial sequencing of VP1. Journal of Clinical Microbiology 37, 1288-1293.
Rueckert, R. R. (1996). Picornaviridae: the viruses and their replication. In Fields Virology , pp. 609-654. Edited by B. N. Fields, D. M. Knipe & P. M. Howley. Philadelphia:LippincottRaven.
Samuelson, A., Westmoreland, D., Eccles, R. & Fox, J. D. (1998). Development and application of a new method for amplification and detection of human rhinovirus RNA. Journal of Virological Methods 71, 197-209.[Medline]
Savolainen, C., Mulders, M. N. & Hovi, T. (2001). Phylogenetic analysis of rhinovirus isolates collected during successive epidemic seasons. Virus Research (in press).
Uncapher, C. R., Dewitt, C. M. & Colonno, R. J. (1991). The major and minor group receptor families contain all but one human rhinovirus serotype. Virology 180, 814-817.[Medline]
Utagawa, E. T., Miyamura, K., Mukoyama, A. & Kono, R. (1982). Neuraminidase-sensitive erythrocyte receptor for enterovirus type 70. Journal of General Virology 63, 141-148.[Abstract]
Vesa, S., Kleemola, M., Blomqvist, S., Takala, A., Kilpi, T. & Hovi, T. (2001). Epidemiology of documented viral respiratory infections and acute otitis media in a cohort of children followed from two to twenty-four months of age. Pediatric Infectious Disease Journal 20, 574-581.[Medline]
Received 21 June 2001;
accepted 25 October 2001.