Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2006, p. 1484-1489, Vol. 44, No. 4
0095-1137/06/$08.00+0 doi:10.1128/JCM.44.4.1484-1489.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Central Virology Laboratory, Public Health Services, Ministry of Health, Chaim Sheba Medical Center, Tel-Hashomer,1 Pediatric Pulmonary Unit, The Edmond and Lily Safra Children's Hospital, Chaim Sheba Medical Center, Tel Hashomer,2 Faculty of life Sciences, Bar Ilan University, Ramat-Gan, Israel3
Received 19 May 2005/ Returned for modification 1 September 2005/ Accepted 1 January 2006
|
|
|---|
|
|
|---|
Sequence analysis of hMPV isolates and infected patients revealed that there are two main genetic genotypes (1 and 2) (30, 41). Subsequently these two genotypes were found to represent four subgroups, 1A, 2A, 1B, and 2B (42). The presence of these two genotypes have been confirmed in many parts of the world (3, 4, 6, 43). More recently, Schildgen et al. have reported a new variant of hMPV isolated from a 6.5-year-old child in Germany (33).
hMPV infections have worldwide distribution. The virus has been detected in patient samples from countries all over the world. hMPV has been reported from Canada (6), Australia (29), Denmark (44), Tunisia (14), France (15), Italy (28), Hong Kong (31), Japan (10), Brazil (8), the United States of America (11), Argentina (16), South Africa (20), Thailand (39), and Israel (46). The frequency of hMPV detection in patients with respiratory tract infections ranged from 2 to 25% (28, 36).
In Southern Israel, Wolf et al. reported the seroprevalence of hMPV among children less than 2 years of age to be 52% (46). To gain a better understanding of the epidemiology of hMPV in Israel, we evaluated the prevalence of hMPV in archived and prospective clinical samples from 388 hospitalized patients with defined respiratory clinical syndromes. Two patient populations were studied: one was negative to most common respiratory viruses for which laboratory diagnosis was sought (influenza A and B, RSV A and B, parainfluenza viruses 1, 2, 3, and adenovirus) and the other population was positive for either RSV or influenza viruses A or B. In addition, we performed phylogenetic analysis on hMPV-positive patient samples to understand which hMPV genotypes have circulated in Israel. Finally, we performed hMPV seroprevalence analysis on the sera of 204 individuals with a wide age distribution who had no clinical respiratory symptoms when the blood samples were obtained.
|
|
|---|
Determination of hMPV seroprevalence by indirect immunofluorescent assay (IFA) was performed on 204 serum samples obtained during 2002 from hospitalized patients with no respiratory clinical symptoms. The distribution of patient samples in each of the groups was as follows: <1 year, 25; 1 to 4 years, 23; 5 to 14 years, 31; 15 to 19 years, 20; 20 to 49 years, 60; 50 to 59 years, 20; >60 years, 25. Patient serum was stored at 70°C.
Cell line and virus strain. Vero cell monolayers in T-25 flasks were maintained in Dulbecco's modified Eagle's medium supplemented with 10% fetal bovine serum, 200 µg/ml streptomycin, 100 U/ml penicillin, and 121.5 U/ml nystatin in a 37°C incubator. hMPV reference strain virus 200026583 (CAN 97-83) passage 3, which was used as a positive control, was graciously supplied by Dean D. Erdman from the Centers for Disease Control and Prevention, Atlanta, Ga. Vero cells (ATCC CCL8) were used to grow the hMPV reference strain. hMPV reference virus was vortexed for 30 s to homogenize samples, and 200 µl was used to inoculate 70% confluent Vero cells in duplicate. Inoculated cells were maintained in serum-free Dulbecco's modified Eagle's medium supplemented with 200 µg/ml streptomycin, 100 U/ml penicillin, 121.5 U/ml nystatin, and 2 µg/ml trypsin type IX (T-0.134; Sigma) and incubated at 34°C. Cell cultures were observed daily for 21 days for the appearance of cytopathic effects. Medium was replaced weekly.
IFA.
IFA was performed as previously described (24, 41) with some modifications. Briefly, Vero cells infected with the hMPV reference strain were mixed with uninfected Vero cells at a ratio of 70% and 30%, respectively. The cells were then spotted on acetone-cleaned glass slides and left to dry at room temperature. The cells on each slide were then fixed with ice-cold acetone for 10 min and stored at 70°C until used. Patient's sera were serially diluted in phosphate-buffered saline and incubated with either hMPV-infected or uninfected (negative control) fixed Vero cells for 45 min at 37°C. The cells were then washed three times with phosphate-buffered saline and incubated with fluorescein isothiocyanate-conjugated AffiniPure goat anti-human immunoglobulin G, Fc
fragment specific (Jackson ImmunoResearch Laboratories, Inc.) for 35 min at 37°C in a humidified box. Cells were then washed and dried, and mounting fluid was added. A Zeiss inverted fluorescent microscope was used to visualize all the samples. Patient samples with a titer greater than 1:32 was considered positive for hMPV antibodies as described by van den Hoogen et al. (41).
RNA extraction. Viral genomic RNA was extracted from patients' sample supernatants using either the QIAamp (QIAGEN GmbH, Hilden, Germany) RNA extraction kit during the 2002 to 2003 respiratory season or the High Pure viral RNA extraction kit (Roche Diagnostics GmbH, Mannheim, Germany) during the 2003 to 2004 respiratory season. Both kits were validated for viral RNA extraction from patient samples. In our hands, the efficiency for viral RNA extraction was similar for both kits used (data not shown). RNA was extracted by following the manufacturers' suggested protocols. Briefly, clinical samples were homogenized by vortexing for 30 s, and 140 µl was used to extract viral RNA by QIAGEN (QIAamp kit) or 200 µl was used to extract viral RNA by Roche (High Pure viral RNA). The RNA was eluted from the columns after several washes and stored at 70°C. A 5-µl aliquot was used for the one-step reverse-transcription (RT)-PCRs.
RT-PCR analysis and sequencing. Two different sets of primers were used in this study. Primer sequences encoding part of the L (polymerase) gene, which were generously provided by B. G. van den Hoogen from Erasmus University Medical Center, Rotterdam, The Netherlands, were used for the diagnosis of patient clinical samples, for determining the presence of hMPV genome, and for phylogenetic analysis. The sequence of the forward (L6) and reverse (L7) primers used were 5'-CAT GCC CAC TAT AAA AGG TCA G-3' (nucleotide positions 11336 to 11357) and 5'-CAC CCC AGT CTT TCT TGA AA-3' (nucleotide positions 11506 to 11486), respectively. A positive control was prepared from the plasmid pCR 2.1, which contains a 171-bp-long insert derived from the hMPV L gene. The plasmid was kindly donated by B. G. van den Hoogen. The insert "RAP PCR fragment 10" was amplified using primers L6 and L7 after digestion of plasmid preps with HindIII and agarose gel purification. The purified fragment was diluted serially to serve as a marker for the sensitivity of the RT-PCR.
Primers sequences MPV01.2 (forward) and MPV02.2 (reverse) derived from the N gene, which were designed by Mackay et al. (25), were used for the phylogenetic analysis of hMPV-positive patient samples and to confirm hMPV patient samples positive by the L gene. The sequence of the forward primer for the N gene was 5'-AAC CGT GTA CTA AGT GAT GCA CTC-3' (MPV01.2), and the sequence of the reverse primer was 5'-CAT TGT TTG ACC GGC CCC ATA A-3' (MPV02.2). The resulting fragment length is 213 bp.
RT-PCR was performed using QIAGEN OneStep RT-PCR kit (QIAGEN GmbH, D-40724 Hilden, Germany). Briefly, 5 µl extracted RNA was added to a master mix composed of enzyme mix (mixture of heterodimeric recombinant reverse transcriptases Omniscript, Sensiscript, and HotStart Taq DNA polymerase), 400 µM concentrations of each deoxynucleoside triphosphates, 20 U RNase inhibitor (CPG, Lincoln Park, N.J.) and the hMPV primers sets at 20 pM final concentration (each). Amplification used the optimized profile provided by B. G. van den Hoogen in the thermal cycler (PTC-100; MJ Research, Watertown, MA): 42°C for 45 min, 95°C for 7 min, followed by 35 amplification cycles (denaturation at 95°C for 1 min, annealing at 45°C for 2 min, and synthesis at 72°C for 3 min) (personal communication). Amplification was completed with a prolonged synthesis at 72°C for 10 min. Amplicons were visualized by ethidium bromide staining following electrophoresis on a 2% agarose gel.
Real-time RT-PCR analysis. TaqMan chemistry was utilized to detect the hMPV genome in patient samples that were positive for either RSV or influenza viruses A and B. hMPV assays conditions were previously described by Maertzdorf et al. (27). The assay was at least as sensitive as the regular RT-PCR with gel detection.
Phylogenetic analysis. RT-PCR products were gel purified or purified directly from the RT-PCR mix and sequenced using ABI PRISM Dye Deoxy Terminator cycle sequencing kit (Applied Biosystems, Foster City, CA). Reaction mixtures were analyzed on Applied Biosystems model 373 DNA automatic sequencing systems. The Sequencher (Gencodes Corporation, Ann Arbor, MI) and University of Wisconsin GCG gene analysis programs (GCG analysis programs, Aclerrys, San Diego, Calif.) were used to compare nucleotide sequences. Phylogenetic trees were prepared by either (i) nearest neighbor analysis using Clustal X with 1,000 bootstraps (40) or (ii) maximum-likelihood phylogenetic analysis using the Dnaml application of PHYLIP (http://evolution.genetics.washington.edu/phylip/getme.html; J. Felsenstein, Department of Genetics, University of Washington) with data bootstrapped 100 times and sequence order randomized thrice for each bootstrap. Trees were visualized using TreeView or NJ plot.
Nucleotide sequence accession number. Partial N and L gene sequences for individual hMPV isolates have been submitted to the EMBL database as a noncontiguous sequence. Sequences have been assigned accession numbers AM056029 to AM056048. The alignment of these sequences compared to isolate 00-1 (AF371337) for reference has been deposited in the EMBL alignment database as ALIGN_000919.
|
|
|---|
Of the 274 negative samples, 42 (15.3%) were positive by RT-PCR for the hMPV L gene. Confirmation of the hMPV-positive samples was done with a different set of primers derived from the N gene. No discrepancy in hMPV-positive results was observed between the two sets of primers. Analyses of each of the two respiratory seasons revealed that 25 (19%) of the 131 samples analyzed during the 2002 to 2003 season and 17 (12%) of the 143 samples analyzed during 2003 to 2004 season were positive for hMPV. Interestingly, none of the patient samples positive for RSV or influenza viruses (n = 114) were positive for hMPV when tested by a real-time RT-PCR assay. This assay was at least as sensitive as the regular RT-PCR assay, and random repeat testing for RSV and influenza viruses for the positive samples was positive with cycle threshold values similar to those initially detected on the original sample, indicating that the RNA was not degraded. Since none of the 114 samples positive for either RSV or influenza were hMPV positive, we estimate that the rate of double infection was below 0.9%. Thus, in our study population, the percentage of positive hMPV RNA from the total specimens submitted to CVL (n = 613) for respiratory viral testing was estimated to be about 6.9%.
hMPV circulated in central Israel between November and April in each of the two seasons studied. Figure 1 shows the monthly distribution of positive hMPV samples. Increased hMPV incidence was observed during February/March for the 2002 to 2003 season, which corresponds to the RSV season in Israel. Main RSV activity was between the months of January and March 2003. On the other hand, hMPV activity was mainly during the months March and April for the 2003 to 2004 season, while the main RSV activity during this season was recorded earlier, between January and March 2004.
![]() View larger version (20K): [in a new window] |
FIG. 1. Histogram of monthly distribution of RT-PCR-positive clinical cases for hMPV RNA among hospitalized patients during the winter seasons 2002 to 2003 (n = 25) and 2003 to 2004 (n = 17).
|
hMPV phylogenetic analysis. Phylogenetic analysis was performed for the N (164 bp) and or L gene (113 bp) from 20 of 42 (48%) hMPV RNA-positive samples isolated between January and December 2003. Sequences from all 20 isolates were highly similar (98 to 100%) to hMPV sequences in the EMBL/GenBank/DDBJ databases. The alignment of sequences from all 20 isolates compared to isolate 00-1 (AF371337) for reference has been deposited in EMBL alignment database as ALIGN_000919. Two main genotypes (groups 1 and 2), with two subgroups each, were reported in several areas of the world (3, 5, 21). Three of the four subgroups also circulated in Israel during 2003, as illustrated by the phylogenetic tree for N gene sequences from 12 randomly selected Israeli isolates and representative subgroup sequences obtained from the EMBL/GenBank/DDBJ databases (Fig. 2). The distribution was 25% subgroup 1A, 65% subgroup 1B, and 10% subgroup 2B. No representatives of subgroup 2A were observed among the 20 samples analyzed. Bootstrapped maximum likelihood phylogenetic trees were prepared for the 11 isolates in which both the L and N gene amplicons were sequenced (Fig. 3). The topology of the phylogenetic tree for the L gene was very similar to the tree of the N gene.
![]() View larger version (40K): [in a new window] |
FIG. 2. Phylogenetic comparison of hMPV N gene nucleotides detected during the respiratory seasons 2002 to 2003 and 2003 to 2004 with sequences from EMBL/GenBank database. The Clustal X nearest neighbor-joining method (number of bootstraps, 1,000) was used to compare 164-nucleotide sequences encoding the hMPV N gene with seven sequences from the EMBL/GenBank database. Subclades 1A, 1B, 2A, and 2B are indicated at the right hand side of the figure.
|
![]() View larger version (30K): [in a new window] |
FIG. 3. Phylogenic analysis of hMPV-positive specimens detected in Israel. The nucleotide sequences of hMPV L gene (113 nucleotides) and the N gene (164 nucleotides) were determined. The maximum likelihood phylogenetic analysis using the Dnaml application of PHYLIP (100 bootstraps) was used to compare different Israeli sequences encoding hMPV L and N genes. Subgroup 1A, 1B, and 2B are indicated.
|
![]() View larger version (10K): [in a new window] |
FIG. 4. Age distribution of patients in central Israel (n = 204) seropositive for hMPV as determined by indirect immunofluorescence assay on samples obtained in 2002. The number of individuals in each age group is indicated on the top of each bar.
|
|
|
|---|
The discovery of hMPV in 2001 by van den Hoogen et al., has helped to identify a human pathogen that is causing some of the undiagnosed respiratory tract infections (41). In our study, by using hMPV RNA detection assays, we have found hMPV to be associated with about 15% of the respiratory illness in hospitalized children negative for other common respiratory pathogens over the two-season study period. At Israel CVL, the incidence of hMPV RNA detection was 6.9% when determined from the total patient samples analyzed in the laboratory. Thus, during the two respiratory seasons studied in central Israel, hMPV is likely to have been the causative agent of at least 6.9% of the upper respiratory tract infections which led to hospitalization, although we cannot rule out other viral or bacterial pathogens as the actual causative agents of the illness. This rate would not change dramatically, since no double infections were detected in 114 (one-third) of the RSV- or influenza-positive patient samples available for analysis from central Israel. Dual infection did not appear to play a major role in the patient populations analyzed during 2002 to 2003 and 2003 to 2004 winter seasons in central Israel. Our data are consistent with the study reported by Lazar from Connecticut, who did not detect any dual infections of hMPV and RSV in their hospitalized patient population with severe or mild RSV infection (23). On the contrary, Semple et al. (34) reported a 10% rate of hMPV and RSV double infection in children hospitalized in the general wards and as high as 72% in infected infants admitted to the pediatric intensive care unit. The dramatic difference in the incidence of dual infection reported by Semple et al. and us could not be due to technical issues, since we used similar molecular techniques. Indeed, our evaluation of dual infections was with a highly sensitive real-time molecular assay that has a sensitivity and specificity at least similar to those of regular RT-PCR. Other groups which examined the rate of hMPV coinfection with other viruses among hospitalized children with acute RTI reported rates of 15 to 40% of the hMPV infections (1, 22). We have not detected hMPV coinfection with other viruses. However, we have not tested our samples for enteroviruses and rhinoviruses, which are common respiratory pathogens which only recently were recognized as contributing to acute expiratory wheezing in infants and children (22).
The data on the incidence of hMPV infections from other countries are inconsistent due to different calculation methods and year-to-year variation within a given region (35). The hMPV positivity rate ranged from 1.5% in Australia, 7.1% in Canada, 20% in the United States, and 25% in Italy (5, 28, 29, 45).
The clinical picture of hMPV-infected children was similar to what was reported in the world (31, 35, 37, 38). The Israeli children infected with hMPV presented with bronchiolitis symptoms similar to those caused by RSV. The symptoms of the hMPV-infected patients ranged from mild self-limiting respiratory illness, with cough, rhinorrhea, and wheezing as the major symptoms, to much more severe illness requiring mechanical ventilation.
Interestingly, there was a month difference between the hMPV seasons during the study period. During the season 2002 to 2003, the RSV and the hMPV season peaked at the same time between January and March 2002. However, during the 2003 to 2004, hMPV season peaked between March and April 2003, while RSV major activity was between the months January and March. Several countries reported different hMPV seasonality. In Canada, Boivin et al. reported hMPV activity in April (6), while a Dutch group reported hMPV activity during December (30).
In our study, patient's specimens were analyzed for the presence of hMPV mainly by RT-PCR to detect the hMPV L gene using conditions which were optimized by van den Hoogen et al. (41). Patient samples positive with the L gene were confirmed by a second PCR using a set of primers for the N gene as described by Mackay et al. (25). Sequence analysis of 20 patients' N and L hMPV genes revealed that the two main hMPV groups 1 and 2 circulated in Israel. Many studies reported the detection of hMPV group 1 as the only or the predominant serotype circulating (2, 3, 6, 28, 31). In our study population, hMPV group 1 also had the highest circulation rate (92% of the sequenced samples). Of the four subgroups reported (3, 4, 21), three were identified in central Israel (1A, 1B, and 2B). hMPV subgroup 1B had the highest circulation rate (65%), followed by 1A (25%) and 2B (10%), of the sequenced samples. No circulation of subgroup 2A was detected during the study period. However, subgroup 2A could be present in central Israel but at a frequency too low to have been detected in the small number of samples analyzed.
The similar topological clustering within the phylogenetic trees of N, P, M, and F genes of individual Canadian and Japanese hMPV strains (2, 21, 23) indicates parallel evolution of these genes. In general, this means that when the subgroup of an isolate is defined for one gene, the other genes will also belong to the same subgroup. The similar topological distribution of N and L genes in the bootstrapped phylogenetic trees of our isolates presented in Fig. 3, especially after collapsing branches with bootstrap values of less than 90%, indicates that a subgroup assignment for an hMPV isolate based on the N gene can be extrapolated to the L gene, as has been done for the P, M, and F genes (2, 21, 23). However, it should be noted that there appears to be less variability in the L gene subgroup 1B than in the corresponding N gene 1B subgroup (Fig. 3).
To determine the seroprevalence of hMPV in central Israel, the serum of hospitalized Israeli patients admitted for symptoms other than respiratory tract-related illness was checked for hMPV antibodies (Fig. 4). The seroprevalence of hMPV in children less than 1 year of age was about 28%. This is similar to the 30% that was reported by Wolf et al. from children less than 13 months of age in the southern part of Israel (46). Moreover, a similar observation was reported by van den Hoogen et al., where the seropositivity of children less than 1 year of age from The Netherlands was about 25% (41). hMPV seroprevalence in Japanese children less than 5 years of age showed a pattern similar to that of the Israeli children (9). As the children's ages increase, hMPV seroprevalence increases reaching 100% within the age group of 5 to 14 years olds. The increase in the seroprevalence is similar to what was observed in The Netherlands (41). Unlike The Netherlands, where the percent seropositivity was at 100% in all age groups older than 5 years, in central Israel, the percent seropositivity after 14 years was less than 100%, ranging from a low of 72% in the age group of 20 to 49 to a high of 80% in individuals older than 60 years of age. The decrease in seroprevalence which occurs between 14 to 60 years might be due to a decline in immune responsiveness with age or due to reduced exposure to hMPV. This is similar to the pattern for RSV and may help explain why some adults are prone to reinfection, but the small sample size limits the interpretation of this finding.
In conclusion, this study is the first report that describes the viral genetic lineages and the prevalence of hMPV infections in central Israel. The discovery of hMPV and its incorporation to routine diagnostic procedures has allowed us to identify the cause of a significant number of respiratory tract infections. More upper RTI will be identified when sensitive diagnostic techniques are developed and routinely applied from testing common respiratory pathogens such as rhinoviruses, and coronavirus. However, there still remain a significant number of infections that have still not been associated with any known pathogen.
|
|
|---|
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»