Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2007, p. 2116-2122, Vol. 45, No. 7
0095-1137/07/$08.00+0 doi:10.1128/JCM.00027-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

University of Würzburg, Institute of Virology and Immunobiology, Versbacher Str. 7, 97078 Würzburg, Germany,1 University of Würzburg, Children's Hospital, Josef-Schneider-Str. 2, 97080 Würzburg, Germany2
Received 4 January 2007/ Returned for modification 19 March 2007/ Accepted 20 April 2007
|
|
|---|
|
|
|---|
Since its first description in 3.1% of 540 Swedish children with lower respiratory tract disease (1), hBoV has been reported worldwide with frequencies ranging from 1.5 to 18.3% in respiratory samples, mostly from children (2, 3, 5, 9, 10, 14, 18, 24). hBoV has been found in patients with a broad variety of both upper and lower respiratory tract diseases and, in several studies, a large number of coinfections with other respiratory pathogens have been reported (1, 5, 6, 9, 15, 18, 24). Thus, at present it is still unclear whether hBoV is a relevant respiratory pathogen or whether it is only an innocent bystander.
The data available thus far on hBoV phylogeny indicate that there is only limited genetic variation between different hBoV strains. Sequences of PCR products located in the NP-1 region have generally shown nucleotide identity of >98% (1-3, 6, 9, 11, 14, 15, 18). However, sequence heterogeneity in the VP1/VP2 region appears to be greater than in the NP-1 region (1, 11). Therefore, a detailed analysis of hBoV phylogeny in the VP1/VP2 region should provide further knowledge on hBoV DNA variation.
The initial studies on hBoV infections have used conventional single-round PCR for hBoV detection, frequently according to the method originally described by Allander et al. (1). In order to facilitate hBoV DNA analysis and to quantify hBoV DNA in clinical samples, we designed a real-time PCR assay. In the present study, we describe the evaluation of the hBoV real-time PCR and its application for the detection of hBoV DNA in nasopharyngeal aspirates (NPA). In addition, we studied stool and serum samples of hBoV-infected children with the real-time assay in order to investigate which body compartments are affected by hBoV infection in addition to the respiratory tract. Furthermore, we sequenced the VP2 gene of hBoV-positive samples with the aim of analyzing the molecular phylogeny of hBoV infections.
|
|
|---|
Qualitative hBoV PCR. DNA was extracted from 200 µl of the NPA and serum samples using the HighPure viral nucleic acid kit (Roche, Mannheim, Germany) and from 140 µl of stool suspensions using the QIAamp viral RNA minikit (QIAGEN, Hilden, Germany). As shown previously, the viral RNA minikit can be used for the simultaneous extraction of RNA and DNA (16, 23). In exceptional cases when the volume of the archived samples was otherwise insufficient, less than the stated starting material was used for the extraction procedures. Both kits were used according to the instructions of the manufacturers. The elution volume of the extractions was 50 µl in case of the NPA and serum samples and 60 µl for the stool samples.
Amplification of hBoV DNA by conventional qualitative PCR was performed with the NP-1 primers BoV188F and BoV542R (Table 1) (1). Reaction conditions and cycling parameters have been described previously (24). Briefly, the PCR was carried out using HotStarTaq DNA polymerase (QIAGEN), and the cycling conditions were 50 cycles (94°C for 30 s, 53°C for 40 s, and 1 min for 72°C) after a preheating step of 10 min at 95°C. Positive PCRs by agarose gel electrophoresis were sequenced for the confirmation of sequence identity. The lower limit of detection with 95% probability was determined to be 11.7 copies/reaction (95% confidence interval, 8.5 to 21.7 copies/reaction) by probit analysis.
|
View this table: [in a new window] |
TABLE 1. Primer and probes used for the amplification of hBoV
|
The real-time PCR was carried out in a final volume of 20 µl consisting of 5 µl of extracted DNA, primers, and probe at a final concentration of 200 nM each, and 1x Quantitect probe master mix (QIAGEN). Amplification was performed on an ABI7500 real-time PCR system (Applied Biosystems, Darmstadt, Germany). The cycling conditions were 50 cycles with 30 s at 95°C and 60 s at 60°C after a preheating step of 15 min at 95°C.
The plasmid pBoV1 containing the PCR product of the qualitative PCR in the vector pCR2.1-TOPO (Invitrogen, Karlsruhe, Germany) was used as a positive control and for the standard curve. Viral loads were calculated from the threshold cycle (CT) values of the individual samples with respect to the standard curve. To minimize background fluorescence, the baseline was set two cycles below the first logarithmic amplification cycle.
General laboratory procedures to prevent PCR contamination were strictly adhered to. One negative control was extracted and amplified for every five NPA samples. All negative controls were found to be negative for hBoV DNA.
Phylogenetic analysis of the VP2 gene. Positive hBoV samples were amplified with the two overlapping primer pairs BoV3885s-BoV4456a and BoV4287s-BoV4939a (Table 1). The PCR was carried out using HotStarTaq DNA polymerase (QIAGEN) with the following cycling conditions: 95°C for 12 min, followed by 52 cycles of 94°C for 30 s, 57°C for 30 s, and 72°C for 30 s, with a final elongation step at 72°C for 7 min. Positive PCRs were sequenced in both directions using BigDye terminator 3.1 chemistry and the ABI Prism 3100 (Applied Biosystems). The sequences were aligned by CLUSTAL W, and a topology tree was constructed with MEGA version 3.1 (12) by using the Kimura two-parameter estimation to generate the nucleotide distance matrix. To evaluate the phylogenetic relevance, bootstrapping analysis with 1,000 replicates was performed by the neighbor-joining method and the p-distance parameter.
Statistical analysis. Statistical analysis was carried out by using GraphPad Prism version 3.0c for Mac (GraphPad Software, San Diego, CA) and SPSS version 13.01 for windows (SPSS, Chicago, IL).
|
|
|---|
The lower limit of detection was determined by using probit analysis. It was based on three independent runs with a total of 24 replicates at the dilutions 37.5, 18.8,15.0, 9.4, 7.5, 3.8, 3.0, 1.5, 0.8, and 0.4 copies/reaction. A concentration of 14 copies/reaction (95% confidence interval, 10 to 26 copies/reaction) was detectable with 95% probability (data not shown). This corresponds to a lower limit of detection of 700 copies per ml of starting material (serum or NPA).
The reliability of the real-time PCR was assessed by analyzing the intra-assay and interassay variation (Table 2). The intra-assay coefficient of variation of the CT values for three replicates was 2.6% or lower for standard concentrations of 3.8 x 108 to 3.8 x 103. The interassay coefficient of variation of a run control (nominal concentration, 3.8 x 105) that was regularly included in each run was 15.21% (data not shown).
|
View this table: [in a new window] |
TABLE 2. Intra-assay variation of the hBoV real-time PCR assay for different concentrations of the plasmid pBoV1
|
Quantitative real-time PCR of NPA. For clinical evaluation of the hBoV real-time PCR, 11 samples from the period of 1998 to 2001 and 823 samples from 2002 to September 2005 NPA that had previously been tested by the conventional qualitative hBoV DNA (24) were retested with the real-time assay. A total of 84 samples were positive, and 724 samples were negative with both assays (Table 3). Discordant results in initial testing were obtained for 26 samples. All discordant samples were retested twice with both assays. An initially discordant sample was regarded as positive for the qualitative or real-time PCR assay, if at least two of the three results obtained with the assay were positive. After retesting, 16 samples remained with discordant results. In order to resolve the discrepancies, the discordant samples were further tested using the two primer pairs designed for the phylogenetic analysis VP2. In all samples positive only by real-time PCR, hBoV DNA was amplified by at least one additional primer pair, thus confirming the positive real-time PCR result. One of the three samples positive only by conventional PCR was confirmed by the additional primer pairs. Overall, the frequency of hBoV detection in our patient population increased from 10.8 to 12.0% with the real-time PCR assay. However, the proportions of positive results obtained with the two methods were not significantly different (P = 0.44 [chi-square test]).
|
View this table: [in a new window] |
TABLE 3. Comparison of qualitative and real-time PCR for hBoV DNA
|
Because it was found previously that a large number of hBoV DNA-positive children were coinfected with other respiratory pathogens [RSV, adenovirus, influenza virus A/B, parainfluenza virus 1/2/3(24)], viral loads in NPA samples of children with coinfections were compared to viral loads of children without detectable coinfection. Of the children positive for hBoV DNA with the real-time PCR assay, 41% (n = 41) were found to be coinfected with RSV (n = 16), adenoviruses (n = 13), influenza A virus (n = 9), parainfluenza viruses (n = 4), or influenza B virus (n = 1), including two triple infections of RSV, parainfluenza virus 3, and hBoV. The hBoV loads in the NPA samples of these children and in the NPA samples of children without coinfection are presented in Fig. 1. Median values between both groups were not significantly different (median with coinfection = 4.3 x 103 copies/ml; median without coinfection = 4.5 x 103 copies/ml; P = 0.759 [Mann-Whitney test]).
![]() View larger version (17K): [in a new window] |
FIG. 1. Viral load in NPA of children with or without respiratory coinfections. In addition to hBoV DNA, all samples were tested for antigens of RSV, influenza virus A/B, adenoviruses, and parainfluenza viruses 1 to 3 by immunofluorescence assay. Median values are indicated by horizontal bars.
|
![]() View larger version (10K): [in a new window] |
FIG. 2. Bocavirus load in NPA according to the clinical diagnosis. The median values (horizontal bars) were 4.3 x 103 copies/ml for upper respiratory tract disease (URTD), 4.0 x 104 copies/ml for pneumonia, 7.1 x 104 copies/ml for bronchitis, 1.6 x 108 copies/ml for bronchiolitis, 2.1 x 102 copies/ml for febrile seizures, 9.3 x 102 copies/ml for wheezing bronchitis, and 4 x 103 copies/ml for all NPA. The bocavirus load difference between the NPA of children with bronchitis and febrile seizures was statistically significant as determined by Kruskal-Wallis analysis and Dunn's post test (indicated by superscript a; P = 0.0195).
|
![]() View larger version (12K): [in a new window] |
FIG. 3. Comparison of hBoV loads in NPA according to the hBoV status of stool samples from the same patients. Median values as indicated by horizontal bars were significantly different between the two groups (P < 0.001 [Mann Whitney test]).
|
![]() View larger version (14K): [in a new window] |
FIG. 4. Phylogenetic analysis of the VP2 gene. Phylogenetic analysis of hBoV VP2 gene sequences of Germany, Sweden (st1 DQ000495; st2 DQ000496), China (CZ643 DQ457413), and the United States (CRD2 DQ340570). The topology tree was constructed in MEGA3.1. Kimura two-parameter estimation was used to generate the nucleotide distance matrix. The phylogenetic relevance was analyzed by bootstrapping (n = 1,000) and is indicated in nodal confidence values. Only values above 60 are displayed. The scale gives the percent nucleotide substitution with respect to branch length.
|
|
|
|---|
One other study using real-time PCR and probe detection for hBoV amplification has been published thus far (13). Lu et al. established two real-time PCR assays with primers and probes located in the NS1 and in the NP-1 region. Sensitivity, amplification efficiency, and reproducibility were similar to the data obtained with our assay. A total of 1,178 throat swab samples from Thai patients with pneumonia were analyzed with their assay. Most of the samples contained only low copy numbers of hBoV DNA as judged by the CT values. However, a thorough quantitative analysis and a comparison of different disease groups was not performed by Lu et al. (13).
The frequency of hBoV infections in our population remains one of the highest reported thus far and is similar to the 11.3% positive samples in a study from Korea (5). Several other studies have described frequencies ranging from 1.5 to 18.3% (2, 3, 5, 9, 10, 14, 18, 24). Differences in the positivity rate for hBoV DNA may be accounted for by different patient characteristics and seasonal and geographical variation. In addition, the results may be influenced by the sensitivity of the hBoV PCR assays used. Because information about this aspect has been provided only in two studies (13, 15), it is not possible at present to estimate how much the reported detection frequencies for hBoV DNA may be affected by assay sensitivity.
Real-time PCR provides the possibility of quantification of DNA copies per reaction with reference to a standard curve. In contrast to serum or plasma, calculation of a viral load in the original patient sample is only of limited value for NPA because of technical variations in the sample acquisition procedure. For example, the volume of the NPA on arrival in the lab ranged approximately from 0.5 and 10 ml. Although this should be borne in mind, the observed viral loads in patient samples differed by a factor of 108, a range much larger than could be accounted for by variation in NPA volumes obtained.
Surprisingly, the median values of the bocavirus loads in NPA of children with or without known coinfection with other viral pathogens were almost identical. This is in contrast to the situation with coronavirus NL63. In children who were only infected with NL63, the viral load for NL63 was significantly higher than in children who were coinfected with NL63 and other viruses (22). With the limitation that additional coinfections remained undetected in our study because antigen-based methods were used for screening of respiratory viruses other than hBoV, and because several respiratory pathogens such as coronaviruses, rhinoviruses, enteroviruses, and the human metapneumovirus were not tested for, our data do not indicate that hBoV replication is influenced by the presence of other respiratory viruses. Also, the finding of comparatively low viral loads in NPA of more than half of the children without coinfection does not support the assumption that hBoV causes respiratory disease in these children.
Thus, whether hBoV is associated with clinical manifestations or is merely an innocent bystander is still an unresolved question. Based on the hypothesis that high viral loads may be an indication of pathogenic relevance, we compared the hBoV load in NPA according to the final diagnosis. A broad range of viral loads, including values of 108 copies/ml and higher, was observed for children with upper respiratory tract disease, pneumonia, bronchitis, and wheezing bronchitis. In contrast, in children with febrile seizures and without prominent respiratory findings the highest viral load was only about 105 copies/ml, and the median of 2.1 x 103 copies/ml was the lowest of all groups. High viral loads were observed in the two children with bronchiolitis. However, because of the small number these results should be interpreted with much caution until more data from children with bronchiolitis are available. Nevertheless, an association between hBoV and bronchiolitis based on qualitative PCR findings has been indicated in other studies (2, 5, 9, 17).
Because hBoV is distantly related to parvovirus B19, a virus highly viremic for several weeks after primary infection, we analyzed serum samples of hBoV-infected children for the presence of hBoV DNA. Only one of ten samples was weakly positive. However, detection of hBoV in NPA samples is not proof of a primary infection, and therefore, we may have missed a viremic phase with the serum samples examined in the present study. Serological studies to detect hBoV-specific immunoglobulin M and immunoglobulin G seroconversion will be necessary to determine the time point of primary hBoV infection. Once this information is available, elucidation of the question of whether hBoV is generally viremic or not will become possible.
hBoV is most closely related to the other two viruses in the genus Bocavirus, the minute virus of canines (MVC) and the bovine parvovirus (BPV). MVC and BPV are known to cause gastrointestinal infections in dogs and calves, respectively (4, 20). Therefore, we also examined available stool samples of children with hBoV-positive NPA. In contrast to the serum results, the stool samples were frequently positive for hBoV DNA (14 of 31 samples). This finding may be explained by replication of hBoV in gastrointestinal epithelium. The observation that the hBoV load in the NPA samples of the children with hBoV-positive stool samples was significantly higher than the hBoV load in the children with hBoV-negative stool samples, leads to an alternative explanation for the presence of hBoV DNA in stool samples. hBoV may be swallowed during respiratory tract infection and subsequently excreted in the feces without further replication steps in the gastrointestinal tract. Whether the hBoV DNA detected in stool samples is derived from infectious or from degraded virus particles and if hBoV may be a cause of diarrhea, remains to be determined once cell culture and animal model systems for hBoV become available.
Phylogenetic analysis of the hBoV-positive NPA samples in our study population confirmed previous findings that hBoV is highly conserved. Although this finding may be biased by the fact that most studies published thus far have used PCR screening methods based on the NP-1 sequence, which is highly conserved, there is currently no evidence that more than one hBoV genotype exists. In contrast, the sequence variation both on the nucleotide and amino acid level is considerably greater for parvovirus B19. Thus, it may be speculated that hBoV has entered the human population more recently than parvovirus B19. Alternatively, the diverse degrees of variation may be due to different replication strategies.
In conclusion, we have successfully established and used a real-time PCR assay for the detection of hBoV DNA. Application of this assay for quantitative analysis has shown that hBoV DNA can be found in NPA samples at very high titers. In addition to the respiratory tract, hBoV was also found in stool samples. However, the clinical relevance of hBoV infections remains uncertain both for the respiratory and for the gastrointestinal tract.
Published ahead of print on 2 May 2007. ![]()
|
|
|---|
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»