Previous Article | Next Article ![]()
Journal of Clinical Microbiology, January 2008, p. 97-102, Vol. 46, No. 1
0095-1137/08/$08.00+0 doi:10.1128/JCM.01117-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Centre for Translational and Applied Genomics, PHSA Labs,1 Children's and Women's Health Centre of BC, PHSA, Vancouver, British Columbia, Canada2
Received 4 June 2007/ Returned for modification 19 July 2007/ Accepted 24 October 2007
|
|
|---|
|
|
|---|
While approximately 35% of all VIRAP specimens have a pathogen identified by DFA, the extended turn-around time and limited detection spectra of those
65% of specimens referred to viral cell culture have led us to examine alternative second-line diagnostic methods. Our prior successful pilot study of the multiplex molecular assays sold by Genaco Biomedical Products (now Qiagen Inc.) (8) led us to a larger prospective study of the method, with parallel specimen processing by VIRAP and ResPlex I/II panels.
(Some preliminary aspects of these data were presented previously, at the Third International Symposium of the Chinese American Association for Clinical Microbiology, Shenzhen, China, October 2006; the First Asian Pacific Symposium on Advanced Molecular Technologies, Hong Kong SAR, China, October 2006; and the CACMID/AMMI Canada Annual Conference, Halifax, Nova Scotia, Canada, 14 to 18 March 2007.)
|
|
|---|
Stored samples were thawed once and immediately processed for nucleic acid extraction. Two-hundred-microliter aliquots of NPW were mixed with an equal volume of phosphate-buffered saline and extracted by a Qiagen BioRobot M48 machine (Qiagen Inc., Hilden, Germany), using the manufacturer's MagAttract M48 virus mini kit. In accordance with the instructions for processing of phosphate-buffered saline-containing samples, Qiagen protease was resuspended in Qiagen AVE buffer instead of protease resuspension buffer. Purified nucleic acids were eluted in a 50-µl volume and frozen at –80°C prior to further study.
Extracts were thawed a single time and analyzed by ResPlex I and ResPlex II panels (Genaco Biomedical Products, Huntsville, AL) in accordance with the supplied protocols. Extracts were analyzed simultaneously by both panels, with 10-µl aliquots of extract serving as templates in 50-µl PCR (ResPlex I) and reverse transcription-PCR (RT-PCR) (ResPlex II) mixtures, using the supplied primer mixes. Qiagen HotStarTaq master mix was used for the PCRs, and Qiagen OneStep RT-PCR master mix was used for the RT-PCRs, without the inclusion of RNasin. Both PCR and RT-PCR mixtures were thermocycled together under the provided RT-PCR conditions (identical to the PCR conditions except for an initial 42°C 30-min RT step, during which the PCRs were effectively idled prior to the hot start). Following thermocycling, amplicons were detected by mixing 5-µl portions of amplification products with either ResPlex I beads (PCR products) or ResPlex II beads (RT-PCR products) in hybridization buffer at 52°C for 10 min. Streptavidin-phycoerythrin conjugate was added, and mixtures were incubated at 52°C for a further 5 min prior to the addition of stop buffer and analysis on a Bio-Rad BioPlex instrument (Bio-Rad Laboratories, Hercules, CA) running Bio-Rad BPM 4.1 software. The BioPlex instrument was calibrated to the high RP1 value, with assays run at 50 beads per class minimum, at a plate temperature set to 52°C. Raw mean fluorescence intensity data from each run were exported to Excel for storage and data analysis. Normally, samples were batch processed from (RT)-PCR through data collection in sets of 46 specimens, with one negative control (distilled H2O amplification template) and one positive control (supplied with the ResPlex reagents) for each run of the ResPlex I and ResPlex II panels.
The assay cutoff value for each of the 21 targets was set at a mean fluorescence intensity of 1,250, as validated by our prior study (8). Both ResPlex I and ResPlex II panels contain a human X chromosome marker as an internal control for the sample collection, nucleic acid extraction, amplification, and detection/classification steps; the lower bound limit for this signal as a control for each specimen was set at the average plus 10 standard deviations of the negative control X chromosome signals for 36 duplicates each for ResPlex I and ResPlex II. A specimen was thus required to generate a minimum X chromosome signal of 415.8 (ResPlex I) or 942.6 (ResPlex II) for data from the respective panel to be accepted for further analysis.
|
|
|---|
For six of the seven pathogen types detected by VIRAP DFA, concordance between the assay methods was good, ranging from 100% agreement (PIV-2) to 84% agreement (RSV-B), with the VIRAP DFA method being taken as the reference; for most of the shared targets, concordance was above 90%. The seventh shared pathogen (adenovirus) gave very low concordance values (<10%) and is discussed further below. As expected for an amplification-based method, the ResPlex panels detected more positive samples for all target species than those detected by VIRAP; these likely originated from low-pathogen-titer specimens below the DFA detection limit. Based on the observed concordance, our earlier validation study including some of the bacterial as well as viral targets (8), preliminary data on ResPlex assay specificity from other groups (25), and specificity evaluation of related Templex-based assays (12), we believe that false-positive results arising from the ResPlex I/II assays are, at most, rare events and, if present at all, are unlikely to skew the large data set analyzed here. A single likely exception to this was observed in the case of PIV-3-positive samples; while there was good concordance between VIRAP and ResPlex results (21/22 [96%] PIV-3-positive samples by VIRAP were also PIV-3 positive by ResPlex II), ResPlex II also reported PIV-1 positivity for 20/21 samples and PIV-4 positivity for 19/21 of the same specimens (not supported by our VIRAP data), indicating a cross-reactivity between the PIV-1, -3, and -4 targets in the ResPlex assay. To ensure that this did not introduce errors into our data set, all 92 samples with any PIV-1, -3, or -4 signal above the cutoff were filtered so as to score positivity only for the highest signal among these three targets. In the majority of cases, one of the three target values was very substantially elevated over the other two, suggesting that this was a valid filter.
A total of 1,177 specimens out of 1,742 scored positive for one or more pathogens by ResPlex assay (68%); the remaining 565 specimens scored negative. A histogram of the number of pathogens detected per sample is shown in Fig. 1. The observed prevalence for all pathogens detected by the ResPlex I/II assay is shown in Fig. 2, in comparison with the prevalence by VIRAP for the seven shared targets. A number of pathogens detected by ResPlex assay but not VIRAP were found to be frequent in our sample population (particularly Streptococcus pneumoniae, Haemophilus influenzae, human metapneumovirus [hMPV], the coxsackie virus/echovirus [CVEV] group, and rhinovirus). None of these is particularly surprising, and our ResPlex observed prevalence values appear to be generally in line with published incidences for all of the target species (1, 6, 7, 13, 17, 19, 20, 26, 29), although some differences among the values reported in these references suggest that regional variation is to be expected.
![]() View larger version (12K): [in a new window] |
FIG. 1. Pathogen prevalence per sample. The histogram shows the number of pathogens detected per sample by ResPlex I/II assay of 1,742 total specimens.
|
![]() View larger version (15K): [in a new window] |
FIG. 2. Pathogen prevalence by ResPlex I/II assay and VIRAP DFA (as obtained with the SimulFluor assay according to the manufacturer's protocols). Light bars, ResPlex I/II prevalence values for all samples; dark bars, VIRAP DFA prevalence values for the same sample set. Note that VIRAP examines only adenovirus, RSV, influenza A virus, influenza B virus, and PIV-1 to PIV-3 and that it groups RSV-A and RSV-B together; the indicated VIRAP RSV-A prevalence should thus be compared against the ResPlex RSV-A and RSV-B prevalence values combined. Abbreviations: M. pneumo, Mycoplasma pneumoniae; C. pneumo, Chlamydia pneumoniae; L. pneumo, Legionella pneumophila; S. pneumo, Streptococcus pneumoniae; N. meningo, Neisseria meningitidis; H. flu 1, Haemophilus influenzae (all types); H. flu 2, Haemophilus influenzae (strains a, b, c, and d); H. flu 3, Haemophilus influenzae (strains e and f); AdV, adenovirus; RSV, respiratory syncytial virus; Flu, influenza virus; PIV, parainfluenza virus; hMPV, human metapneumovirus; CVEV, coxsackie virus/echovirus family; Rhino, rhinovirus.
|
|
View this table: [in a new window] |
TABLE 1. Cross-correlation of pathogen prevalence rates by ResPlex assaya
|
In order to select only those values for which meaningful interpretations can be made, all pathogen pairs in Table 1 were evaluated through 2 x 2 Pearson's chi-square analysis with Bonferroni-Holm P value correction for multiple comparisons (22); in cases where Pearson's chi-square analysis was constrained by low expect values, Fisher's exact test was employed. Eight pathogen pairs resulting in significant P values (<0.05) by this method are indicated in Table 1 and given in Table 2, along with the test method employed, the uncorrected P value observed, and the Bonferroni-Holm-corrected P value for the interaction. The uncorrected P value results were also analyzed with the Benjamini-Hochberg adjustment (5), which indicated a further five significant interactions, as indicated in Tables 1 and 2. Following this analysis, in order to examine our data for weaker but still potentially suggestive interactions, we took an approach of calculating the 90% confidence intervals for each Table 1 cell's prevalence value and for the overall prevalence frequency of each target pathogen; if these values were nonoverlapping, the result was considered suggestive. While this method lacks some statistical rigor due to the fact that the confidence intervals were not calculated for totally independent populations, it is simple and analogous to how one might consider these data if they were presented as a bar graph. Cell values (high and low) flagged as suggestive by this approach are also highlighted in Table 1. As expected, all interactions determined as significant by the more rigorous first approach were also detected by this method, and a small number of additional pathogen interaction pairs were highlighted as being of potential interest. We comment further on the relevance of this observation in Discussion.
|
View this table: [in a new window] |
TABLE 2. Pathogen pairs with significant codetection ratesa
|
|
|
|---|
Firstly, codetection of pathogens was the norm, not the exception, among our acute-phase respiratory samples. Our data averaged 1.02 pathogens per sample, including 565 specimens reporting no pathogen; it would not be unreasonable to assume that these nonclassified specimens had, as a minimum, one unidentified respiratory pathogen responsible for the symptomatic presentation. If this conservative estimate of single-pathogen infection is true for these 565 samples, our actual average number of pathogens per specimen value rises to 1.34. The fact that such unidentified pathogens are present in our sample set and would by inclusion increase the average number of pathogens per sample is hardly speculative; the ResPlex I/II assay panels do not include several important pathogens, including bocavirus and members of the coronavirus family (such as strains OC43, 229E, and NL63), which are reported to occur at frequencies of approximately 12% and 5 to 10%, respectively (13, 14, 15, 26), and would thus certainly be expected to be present within our samples. Our results in this regard differ markedly from those of Weigl and coworkers (27), who observed coinfections in a remarkably small percentage of specimens in their study. This may be attributable, to a large extent, to two factors. The first is methodological differences, as Weigl et al. employed a classical high-multiplicity RT-PCR technique; aside from the difficulty in optimizing all of the individual reactions within such a multiplex reaction to work under the same thermocycling conditions, these assays can suffer from a squelching effect whereby one successful target amplification in a sample can suppress amplification of other targets present and thus have a tendency to underrepresent codetection frequencies. In contrast, the Templex ("QiaPlex") method used here is effectively a simplex PCR which avoids this squelching effect (12) and can thus more accurately detect the simultaneous presence of multiple assay targets. The second factor may be that Weigl et al. were primarily examining viral targets, for which we also observed low codetection frequencies, as commented on below.
Our data provide statistical evidence demonstrating that coinfections are not random; there are clear correlations for the occurrence of certain pathogens, as shown in Tables 1 and 2. The major infectivity-enhancing interactions we observed included an enhancement between S. pneumoniae and H. influenzae, a very large (and to our knowledge previously unreported) enhancement between Neisseria meningitidis and influenza B virus (on the order of eightfold, either way), and an enhancement between the CVEV group and the rhinovirus group. While the mechanisms behind these enhancements may vary from case to case, a number of plausible candidates have been discussed widely (for reviews, see references 4 and 10).
Beyond the 13 strongly supported interactions shown in Table 2, we saw a number of weaker but suggestive interactions indicated by our second analysis method, as indicated in Table 1. While our data set is not large enough to give strong statistical support to any of these, we believe the example of S. pneumoniae in an RSV-B-infected background is potentially meaningful in this regard and warrants additional discussion. Our data show that in a background of RSV-B positivity, S. pneumoniae prevalence is 26%, which is greater than its overall prevalence of 20% in our samples. This interaction approached but did not cross the threshold for significance (P < 0.05) by our first analysis method and escaped, by a tiny margin (1.3%), being flagged as suggestive by our second confidence interval approach. The reverse situation—RSV-B positivity in an S. pneumoniae background—does not meet significance criteria for either approach. A review of the literature not only reveals long-standing epidemiological evidence for coinfection by these two agents (2, 10, 21) but also provides convincingly demonstrated mechanisms centered around the ability of S. pneumoniae adherence to the RSV-B glycoprotein G (gG) expressed either on the surfaces of infected cells or on free virions to enhance S. pneumoniae adherence and infection (11). Either route enables previously present RSV-B to increase the infectivity of S. pneumoniae. Our elevated prevalence value in this case thus detects a documented phenomenon, despite not being quite large enough to be picked by either of our statistical analysis methods. We believe that this example highlights that even a relatively large data set such as the one presented here may not be large enough to detect some real interactions and thus that beyond the well-supported interactions shown in Table 2, the weaker interactions indicated as suggestive by our second screening method may be worthy of further directed study.
In all cases of observed enhancement, a major concern is whether these may be spurious results arising from cross-detection of the involved targets; however, as discussed previously, we have not seen evidence of cross-reaction of these target sets. In the case of the observed CVEV/rhinovirus targets, the large number of serotypes represented in each group means that our prior tests for cross-reactivity have not been exhaustive, and it is conceivable that some serovars of each group could cross-react with the target for the other. We can, however, state that in many samples where a strong positive signal was seen for one or the other of these two targets, the signal was restricted to that target only; thus, our data indicate that any such cross-reactivity cannot be a general occurrence, if it happens at all.
Our data include cases of clearly supported or suggestive pathogen cosuppression, where infection with one pathogen reduces the risk for infection by the other. Interestingly, most of these interactions (Table 1) are reciprocal and occur between single-stranded RNA viruses. One possible explanation for this could be activation of nonspecific antiviral functions (including protein kinase R, the 2-5A system, Mx proteins, and apoptotic pathways) (3, 18, 28) by the first infecting agent; inhibition of initiation of infection by a second RNA virus in the face of multiple already activated antiviral responses would be a plausible outcome.
All of these data beg the question of whether codetection of multiple respiratory pathogens in a sample has clinical relevance. In particular, while both H. influenzae and S. pneumoniae can be important pathogens, they are also well known to occur in apathogenic carrier states; more detailed clinical analysis of outcomes when they occur simultaneously with other respiratory pathogens is clearly desirable for our understanding of how to interpret their detection in these cases. An initial line of evidence comes from a review of particularly sick pediatric patients in our intensive care unit who had samples included in this study. Seven such patients with clinical diagnoses of RSV bronchiolitis were selected based purely on severity of presentation; when the ResPlex data for these cases were examined, we found that five of seven patients had multiple pathogens present (one patient had no pathogen detected; one patient was positive for RSV alone; one patient was positive for S. pneumoniae plus influenza A virus; two patients were positive for H. influenzae, S. pneumoniae, and RSV; and two patients were positive for H. influenzae, RSV, and CVEV). While preliminary, these data suggest that flow-scheme diagnostics, which would have ceased with diagnosis of a single viral pathogen in these examples, might have led to exclusion of antibiotic therapy when it was potentially warranted.
We hope that the data presented herein will be a starting point for further examination of this topic. While our total specimen set was of appreciable size, some of the individual positivity subgroups were too small to yield meaningful data; there are several cell values in Table 1 which appear to be highly suggestive yet do not reach our definition of statistical relevance with the current sample size. Until a much larger data set examining respiratory specimens for multiple pathogens is available, a more thorough interpretation of the results presented here is therefore difficult. Based on our results, we concur with the recent opinion of Freymuth et al. (9) and feel that DFA backed with a multiplex molecular method is the current best practice in respiratory diagnostics. As such, we hope that the needed data will increasingly be available in the future as we and other clinical diagnostic laboratories move toward this approach.
We thank Jacob Rozmus for his assistance in selecting the subset of pediatric intensive care unit cases based solely on severity of clinical presentation, as described in the Discussion.
Published ahead of print on 31 October 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»