- Open Access
Virus detection and identification using random multiplex (RT)-PCR with 3'-locked random primers
© Clem et al; licensee BioMed Central Ltd. 2007
- Received: 27 April 2007
- Accepted: 28 June 2007
- Published: 28 June 2007
PCR-based detection and identification of viruses assumes a known, relatively stable genome. Unfortunately, high mutation rates may lead to extensive changes in viral nucleic acid sequences making dedicated PCR primer use problematic. Furthermore, in bioterrorism, viral consensus sequences can be genetically modified as a countermeasure to RT-PCR and DNA chip detection. Accordingly, there is a great need for the development of rapid and universal virus detection and identification technologies.
We report herein that viral genomic DNA or RNA can be separated from host nucleic acids in plasma by filtration and nuclease digestion, and randomly amplified in a single PCR using a mixture of primers designed to be resistant to primer-dimer amplification (5'-VVVVVVVVAA-3', V = A, G or C; 38 or 6561 primers). We have termed this novel PCR method Random Multiplex (RT)-PCR since hundreds of overlapping PCR amplifications occur simultaneously. Using this method, we have successfullydetected and partially sequenced 3 separate viruses in human plasma without using virus-specific reagents (i.e., Adenovirus Type 17, Coxsackievirus A7, and Respiratory Syncytial Virus B). The method is sensitive to ~1000 genome equivalents/ml and may represent the fastest means of detection of unknown viruses.
These studies suggest that the further development of random multiplex (RT)-PCR may lead to a diagnostic assay that can universally detect viruses in donated blood products as well as in patients suffering with idiopathic disease states of possible viral etiology.
- Human Plasma
- Respiratory Syncytial Virus
- Severe Acute Respiratory Syndrome
- Yellow Fever
- Severe Acute Respiratory Syndrome
Relatively benign viruses can be converted into highly virulent viruses via the introduction of genes of interest. For example, Ectromelia virus, a natural pathogen of mice that causes mousepox, recently was recombined with interleukin-4 as part of an effort to develop a live virus immuno-contraceptive vaccine. Surprisingly, the recombined virus caused 60% mortality in 2 strains of mice, whereas the wild type virus caused no death . A credible bioterrorism scenario might entail the release of such a recombined or chimeric virus tailored for maximum infectivity and pathogenicity but not readily detectable using our current "state-of-the-art" diagnostics (i.e., PCR and micro-array chips.) Accordingly, there is a need for methods that can identify unknown viral pathogens and which can reveal extensive genomic information. Such methods would not only prove useful for our defense against bioterrorism, but also would improve our capacities to identify and control outbreaks of naturally occurring pathogenic viruses.
The proven technology to rapidly detect and identify known human pathogens as potential causes of disease or as bioweapons is PCR . A recent example involved a case of fatal yellow fever (YF) in a traveler returning from Amazonas, Brazil in March, 2002 . The otherwise healthy 47-year-old male developed fever, headache, pancytopenia, coagulopathy, renal and liver failure. Pan-cultures were negative, and a peripheral smear yielded no plasmodia. Serological tests (IgG, IgM) performed at the CDC were negative for YF, dengue, St. Louis encephalitis, and other agents, but serum specimens (unfortunately examined only post-mortem) were positive for YF viral RNA as measured by RT-PCR . Subsequent attempts to isolate or culture the virus failed. This example highlights the superiority of PCR over other currently available methods. DNA chips that allow the simultaneous measurement of literally thousands of genes through hybridization are now being developed as the next-generation rapid diagnostic test for all known human pathogens . However, both of these technologies rely on a relatively stable genome, and several human pathogens display a high mutation rate (e.g., HIV 2 × 10-5/base, 9 kB genome ). Moreover, the ability to recombine "non-pathogenic" viruses in vitro introduces the potential not only for de novo pathogenicity but also for enhanced stealth. These considerations suggest that it may be impossible to design DNA micro-arrays which detect nucleic acids from all known and unknown viruses, including less obvious vehicles of bioterrorism such as adenovirus or rhinovirus recombined with a single gene for enhanced virulence. Accordingly, there is a great need for the development of techniques that enable the universal amplification of viral nucleic acids.
There are 2 strategies to amplify genetic sequences with PCR without prior knowledge of the precise sequences. The first strategy relies on the degenerate binding of arbitrarily chosen primers to sample multiple cDNAs or genomic DNA species during PCR (two related methods are known as differential display and RNA Arbitrarily Primed PCR [RAP-PCR] [6–8]). The amplicons range in size from ~50–600 bp and overlap. RAP-PCR was recently used successfully to identify a novel human pneumovirus only after the virus had been cultured . These techniques yield an amplicon "fingerprint" and are generally used to compare two populations of nucleic acids. Accordingly, the pneumovirus was identified by comparing infected cells with non-infected cells . The other strategy is known as Sequence-Independent, Single-Primer Amplification (SISPA) and involves the directional ligation of a linker/adaptor oligonucleotide onto both ends of a target population of either DS DNA or DS cDNA (after reverse transcription from RNA) . The PCR is accomplished with primers specific for the linked adaptor molecule. SISPA, followed by extensive immunoscreening, enabled the cloning of the Norwalk agent genome from stool and Hepatitis G virus from plasma . Unfortunately, these techniques have not found great utility in the direct discovery or identification of novel viruses from human plasma or serum samples (unless accompanied by immunoscreening) because of the relatively high level of host genomic DNA. However, recent studies indicate that DNAse treatment of serum samples can degrade most of the host genomic DNA but not affect viral nucleic acids, which are protected by stable capsids . The combination of DNAse treatment with SISPA was successfully used to identify and sequence in their entirety 2 novel parvoviruses from bovine serum .
While the aforementioned PCR techniques are useful for the discovery of novel viruses, their utility to rapidly detect and identify unknown (or recombined) viruses is quite limited. RAP-PCR, random RT-PCR and SISPA are complicated multi-step procedures requiring infection of cultured cells, days to accomplish, expensive reagents and broad technical expertise. Accordingly, these methods have not found clinical utility for the rapid detection or identification of viral agents.
Here, we present a novel method for rapid virus detection and identification based on random multiplex (RT)-PCR using 3'-locked random primers to avoid primer-dimer amplification. Once detected, virus amplification products can be shot-gun cloned and sequenced for identification. This method may prove useful for the rapid detection and identification of emerging or recombined viruses.
Filtration and Nuclease Treatment Improves Viral DNA Amplification by Removing Host DNA and mRNA
5'-VVVVVVVVAA-3' Primers Enable Random Multiplex Amplification of Plasmid DNA
We next designed a mixture of random primers that would enable the priming and amplification of all DNA or cDNA molecules, but that would be resistant to primer-dimer amplification. We chose 10 base pairs as a compromise length to simultaneously optimize hybridization and possible binding sites. Taq polymerase requires that the three 3' bases of a primer be 100% complementary in order to enable efficient polymerization. We postulated that "locking" the 3' end of the primers with 2 adenosines and not incorporating any thymidines into the upstream 5' portion of the primers would prevent primer-dimer amplification (i.e., 5'-VVVVVVVVAA-3'; where V = A, G or C). This primer design will limit the frequency of binding sites within a particular viral genome, but using the EditSeq DNA Editing application (DNAStar), we observed that primers within this mixture will bind every 200–500 bp in both orientations in the 19 viral genomes that were examined, thus providing ample hybridization and amplification potential.
5'-VVVVVVVVAA-3' Primers Allow For the Random Multiplex Amplification of DNA from Human Plasma Inoculated with Adenovirus Type 17
Sensitivity of the Random Multiplex PCR Using 5'-VVVVVVVVAA-3' Primers
The suspension of Adenovirus Type 17 received from the ATCC had not been quantified and the sensitivity of the random multiplex PCR method thus was uncertain. In order to generate a standard curve, we conducted real-time PCR of a known concentration of PCR-amplified Adenovirus Type 17 genomic template using Adenovirus Type 17-specific primers and found that we could detect as little as 10 copies/ml of template (data not shown). We then subjected the ATCC Adenovirus Type 17 suspension from the ATCC to real-time PCR using the same Adenovirus Type 17-specific primers and found that we could readily detect the genomic DNA (data not shown). Based on the real-time PCR generated from the known concentrations of Adenovirus Type 17 genomic template, we generated a standard curve plotting the concentration of Adenovirus Type 17 genome equivalents versus the cycle number when fluorescence passed the pre-set threshold (Ct; Figure 4A). Based on this standard curve, we determined that the stock solution of Adenovirus Type 17 contained 1 × 107 genome equivalents/ml. We then titrated down the concentration of virus in human plasma and conducted the combination of filtration, DNAse/RNAse treatment and random multiplex PCR using the 5'-VVVVVVVVAA-3' primers. We found that the random multiplex (RT)-PCR method is sensitive to 1000 genome equivalents/ml (Figure 4B).
Application of Random Multiplex PCR to RNA Viruses
In November 2002, an outbreak of atypical pneumonia occurred in Guangdong Province, People's Republic of China [13, 14]. The disease did not respond to empirical antibiotic treatment, and no known bacterial or viral pathogens were identified using serological and (RT)-PCR analyses. The disorder was called Severe Acute Respiratory Syndrome (SARS), and by July 2003 it had caused 8439 reported cases and 812 fatalities [13, 14]. In an effort to identify the causative agent of SARS, the WHO SARS Aetiology Study Group coordinated the distribution of information and infected materials and the analyses of the data amongst eleven separate laboratories . Researchers cultured fetal rhesus kidney cells with infected nasopharyngeal aspirates and then conducted random RT-PCR using the primer, 5'-GCCGGAGCTCTGCAGAATTCNNNNNN-3' for reverse transcription, and 5'-GCCGGAGCTCTGCAGAATTC-3' for PCR. Unique PCR products in the infected cell preparation were cloned and sequenced, and 57% homology was found with bovine Coronavirus. RT-PCR of hundreds of infected clinical specimens, using primers specific for the novel virus, confirmed the identification [13, 14]. While this approach was successful, it is worthwhile to note that it required 5 months and 11 laboratories to identify the virus. By contrast, the procedures described herein are – at least theoretically – capable of identifying an unknown virus without the use of specially designed primers. Furthermore, is possible to at least partially identify an unknown agent within 2–3 days in one laboratory without intermediate culture steps.
In the current study, we have detected and identified 3 separate viruses (i.e., Adenovirus Type 17, Coxsackievirus A7, Respiratory Syncytial Virus B), using the novel approach of random multiplex (RT)-PCR of purified viral nucleic acids. The detection was accomplished in 7 hours using only a micro-centrifuge and PCR machine and did not require any virus-specific reagents. The sensitivity is ~1000 genome equivalents/ml for Adenovirus Type 17 which is not as sensitive as virus-specific PCR (<50 genome equivalents/ml) but should be adequate for the detection of free virus in untreated individuals. Coupled with the potential for direct sequence information, we are confident that this approach may prove useful for our defense against acts of bioterrorism, as well as for the detection and characterization of novel viruses in blood products and in patients displaying hallmarks of infection (e.g., leukocytosis, fever).
Although we believe that the method is adequate in current form to test clinical specimens in a randomized, blinded diagnostic trial, we consider that certain alternative approaches may yield a method that has improved sensitivity and specificity. In particular, we intend to decrease the background amplification that has been variably observed in order to prevent false positives and to increase sensitivity by applying SYBR green fluorescence for DNA detection. Finally, we will apply the optimized approach to blood specimens from virally infected patients in order to detect unsuspected variables and to validate the method within the clinical environment.
Our long-term goal is to provide emergency room and primary care physicians with the ability to universally detect viruses in human plasma. The rationale is, in part, that if an emerging or recombined pathogenic virus were intentionally released into a population, then the ability of physicians to detect the pathogen will be crucial for containment. A second and related application of this technique involves identification of previously undescribed – but naturally occurring – viruses in patients suffering with idiopathic disease states of possible viral etiology. A third application of these procedures may be in the screening of blood bank blood just for the presence of virus. Although the precise identification of the virus requires cloning and sequencing, the amplification pattern will provide an immediate signature indicating the presence of virus. Most clinical laboratories that service emergency rooms and primary care clinics possess the ability to conduct RT-PCR for known viruses such as HIV, and the unmodified random multiplex (RT)-PCR method does not require any further technical expertise. In conclusion, the optimization and testing of the random multiplex (RT)-PCR technology should allow for the development of a universal virus detection assay that will vastly improve our collective defense against viral bioweapons and against emerging, heretofore unidentified natural viruses.
PCR of plasmids using 5'-VVVVVVVVAA-3' and 5'-NNNNNNNNNN-3' primers
Plasmids (10 ng each) pEGFP (Clontech), pSuperNeo (Oligoengine), Pet 30b (Calbiochem), and pBabe (Addgene) were added to a PCR mix containing (per reaction): 1 μl of 40 μM of one of the primers indicated (Oligos Etc), 2.5 μl 10× PCR buffer (Promega), 5 μl (5 mM) MgCl2, 1 μl (5 mM) dNTPs, 14 μl water, and 0.6 μl Taq polymerase (Promega). PCR was performed according to the following protocol on an Eppendorf Mastercycler Gradient: 95°C-2 minutes and then 94°C-30 seconds, 33°C-59 seconds, 34°C-6 seconds (increase in increments of 1°C for 6 seconds each up to 57°C for 6 seconds), 58°C-2 minutes (+ .3°/s), repeat 50 times, 37°C-5 minutes, 65°C-5 minutes (+ .3°/s), and then hold at 4°C. Products were run on a 1% agarose gel containing ethidium bromide and images were captured on a UVP, Bio-Doc It, gel documentation system.
Adenovirus Type 17 was purchased from ATCC (VR-1094). This virus is a dsDNA virus. The source was conjunctival scraping from a 1-year old female, Saudi Arabia, 1955. Coxsackievirus A7 W.P. Parker was purchased from ATCC (VR-166). This virus is a ssRNA virus. The source was stool of a 7-year old with fever, headache, sore throat, stiff neck, and back pain, New York, 1949. Respiratory syncytial virus was purchased from ATCC (VR-1580). This virus is a ssRNA, linear virus. The source was the respiratory secretions from a child with acute respiratory disease.
Blood Plasma Isolation and Viral Separation
10 ml of heparinized whole blood was placed into a 15 ml Falcon tube (BD Biosciences) and centrifuged at 1700 × g for 10 minutes. 250 μl of plasma layer was added to 250 μl room temperature 1× PBS and placed in a SPIN-X .22 micron, 2 ml centrifuge tube (Costar). Plasma was then inoculated with virus and the mixture was centrifuged at 2500 × g for 15 minutes. If any of the mixture remained above the filter, centrifugation was continued at higher speeds until all plasma had passed the filter into the collection tube. 500 μl of PBS was added to the filtered mixture and treated with 7 μl each of un-diluted Riboshredder (Epicentre) and Omnicleave (Epicentre). Treatment of samples with the RNAse and DNAse should cleave any DNA or RNA not protected by viral capsid. All samples were placed in a Thermomixer R (Eppendorf) at 37°C and shaken for 15 seconds every 3 minutes for 2 hours. Plasma was then ready for viral isolation.
Isolation of Virus from Plasma
1 ml of blood plasma ± virus was cooled to room temperature and placed into a 2 ml tube. Viral isolation was carried out according to the manufacturer's protocol (Qiagen). Viral DNA or RNA was eluted off the column in a final volume of 30 μl elution buffer. This product was either used directly for random PCR or used in the reverse transcription reaction.
Reverse Transcription of Viral RNA
Reverse transcription was carried out according to the manufacturer's protocol either with the Omniscript RT (Qiagen) or SuperScript II RT (Invitrogen). The maximum amount of viral RNA was used in place of water for both reactions.
Isolation and cloning of random PCR products
The PCR of viral DNA isolated from whole blood was done according to the method described above utilizing the 5'-VVVVVVVVAA-3' primers. Products were run on a 1% agarose gel containing ethidium bromide. The bands which differed from the water control were excised and the DNA was extracted using the Qiaquick gel extraction kit (Qiagen). DNA was eluted off in a final volume of 20 μl elution buffer. The DNA was cloned into the TOPO cloning vector (Invitrogen) following the manufacturer's protocol and transformation was carried out in TOP 10 chemically competent E. coli. Resulting colonies were picked randomly and mixed with 20 μl PBS. 1 μl of this mixture was used in the PCR reaction using the universal M13 forward and reverse (1 μl of 10 μM) primers and 2 × PCR master mix (Promega). The following parameters were used for colony screening: 95°C-2 minutes, and then 36 cycles of 94°C-30 seconds, 55°C-30 seconds, 72°C-1 minute. PCR products were run out on a 1% agarose gel containing ethidium bromide. Those colonies which yielded products greater than 200 bp were grown in LB broth overnight and the DNA was isolated according to the manufacturer's protocol (Qiagen) eluting the DNA off in a final volume of 25 μl. DNA was sequenced using the universal M13 forward and reverse primers at the University of Louisville's core sequencing facility. Resulting sequences were screened using the BLAST program on the NCBI web site to determine sequence origin.
PCR from Viral DNA or cDNA with Virus-specific Primers
PCR was carried out with either viral DNA or cDNA using a 2× PCR Master Mix (Promega) and 1 μl of viral DNA. 1 μl of 10 μM primers was used for each reaction. The total volume was equal to 25 μl per reaction. The following parameters were used for viral-specific amplification: 94°C-5 minutes, 94°C-30 seconds, 55°C-30 seconds, 72°C-2 minutes, repeat 45 times, 72°C-7 minutes.
Virus-Specific Primers (standard PCR)
For generation of 200 bp fragment:
Coxsackievirus forward primer (923): 5'-TTATCAGAGATGGCAGCACC-3'
Coxsackievirus reverse primer (1105C): 5'-CTTGTCCACCGCTGTAGCCT-3'
For generation of 1.4 Kb fragment:
Coxsackievirus forward primer (923): 5'-TTATCAGAGATGGCAGCACC-3'
Coxsackievirus reverse primer (2161C): 5'-CATTGCCTGCATTCTGTTGG-3'
Virus-Specific Primers (For Real Time PCR)
Coxsackievirus forward primer (937): 5'-AGCACCACTGCAATCACCGA-3'
Coxsackievirus reverse primer (1094C): 5'-CTGTAGCGTCAGTGTCAGGA-3'
Adenovirus forward primer (5770): 5'-TGTAGGTGTAGGCCACGTGA-3'
Adenovirus reverse primer (6084C): 5'-TCGCCAAGCTTCTCTCCAAC-3'
Template Production Primers:
Adenovirus forward primer (5492): 5'-CTCTTACCTCGCGTCTCCAT-3'
Adenovirus reverse primer (6084C): 5'-TCGCCAAGCTTCTCTCCAAC-3'
Coxsackievirus forward primer (923): 5'-TTATCAGAGATGGCAGCACC-3'
Coxsackievirus reverse primer (1105C): 5'-CTTGTCCACCGCTGTAGCCT-3
Generating viral DNA templates for real time PCR
Viral DNA isolated from whole blood plasma was used to generate a template for use in real time PCR. 2 μl of DNA from the virus was used with Platinum PCR Master Mix (Invitrogen) and 2 μl of a 10 μM solution of each of the template primers outlined above. The total volume of the reaction was 50 μl. The same PCR parameters were used as were used for viral-specific amplification. The resulting products were run on a 1% agarose gel containing ethidium bromide. The bands were excised and DNA was eluted according to manufacturer's instructions (Qiagen).
Real Time PCR
All real time PCR samples were run on a Cepheid Smart Cycler. OmniMix Hs Reagent (Cephied, packged by Takara) was re-suspended in 50 μl of water, which was sufficient for two PCR reactions. 1 μl of the stock virus and 1 μl of each of the template primers (10 μM) were added and mixed. Samples were run at 95°C for 30 sec, 54°C for 15 sec, and 72°C for 20 sec for a total of 40 cycles.
Determination of viral molecules in viral template
The molecular weight was determined for each of the dsDNA templates generated by PCR. The OD of the resultant DNA was also recorded. The OD, in grams/μl was divided by the molecular weight of the product and the result was equal to moles/μl. This number was then multiplied by Avogadro's number (6.022 × 1023 molecules/mole) and the resulting number is equal to the number of viral DNA molecules per microliter.
This study was funded by a grant from the Center for Genetics and Molecular Medicine, University of Louisville.
- Jackson RJ, Ramsay AJ, Christensen CD, Beaton S, Hall DF, Ramshaw IA: Expression of mouse interleukin-4 by a recombinant ectromelia virus suppresses cytolytic lymphocyte responses and overcomes genetic resistance to mousepox. J Virol. 2001, 75 (3): 1205-1210. 10.1128/JVI.75.3.1205-1210.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Espy MJ, Cockerill IF, Meyer RF, Bowen MD, Poland GA, Hadfield TL, Smith TF: Detection of smallpox virus DNA by LightCycler PCR. J Clin Microbiol. 2002, 40 (6): 1985-1988. 10.1128/JCM.40.6.1985-1988.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Fatal yellow fever in a traveler returning from Amazonas, Brazil, 2002. Jama. 2002, 287 (19): 2499-2500. 10.1001/jama.287.19.2499.Google Scholar
- Wilson WJ, Strout CL, DeSantis TZ, Stilwell JL, Carrano AV, Andersen GL: Sequence-specific identification of 18 pathogenic microorganisms using microarray technology. Mol Cell Probes. 2002, 16 (2): 119-127. 10.1006/mcpr.2001.0397.PubMedView ArticleGoogle Scholar
- Seo TK, Thorne JL, Hasegawa M, Kishino H: Estimation of effective population size of HIV-1 within a host: a pseudomaximum-likelihood approach. Genetics. 2002, 160 (4): 1283-1293.PubMedPubMed CentralGoogle Scholar
- Liang P, Pardee AB: Differential display of eukaryotic messenger RNA by means of the polymerase chain reaction. Science. 1992, 257 (5072): 967-971. 10.1126/science.1354393.PubMedView ArticleGoogle Scholar
- Nie X, Singh RP: A novel usage of random primers for multiplex RT-PCR detection of virus and viroid in aphids, leaves, and tubers. J Virol Methods. 2001, 91 (1): 37-49. 10.1016/S0166-0934(00)00242-1.PubMedView ArticleGoogle Scholar
- Ralph D, McClelland M, Welsh J: RNA fingerprinting using arbitrarily primed PCR identifies differentially regulated RNAs in mink lung (Mv1Lu) cells growth arrested by transforming growth factor beta 1. Proc Natl Acad Sci U S A. 1993, 90 (22): 10710-10714. 10.1073/pnas.90.22.10710.PubMedPubMed CentralView ArticleGoogle Scholar
- van den Hoogen BG, de Jong JC, Groen J, Kuiken T, de Groot R, Fouchier RA, Osterhaus AD: A newly discovered human pneumovirus isolated from young children with respiratory tract disease. Nat Med. 2001, 7 (6): 719-724. 10.1038/89098.PubMedView ArticleGoogle Scholar
- Reyes GR, Kim JP: Sequence-independent, single-primer amplification (SISPA) of complex DNA populations. Mol Cell Probes. 1991, 5 (6): 473-481. 10.1016/S0890-8508(05)80020-9.PubMedView ArticleGoogle Scholar
- Matsui SM, Kim JP, Greenberg HB, Su W, Sun Q, Johnson PC, DuPont HL, Oshiro LS, Reyes GR: The isolation and characterization of a Norwalk virus-specific cDNA. J Clin Invest. 1991, 87 (4): 1456-1461.PubMedPubMed CentralView ArticleGoogle Scholar
- Allander T, Emerson SU, Engle RE, Purcell RH, Bukh J: A virus discovery method incorporating DNase treatment and its application to the identification of two bovine parvovirus species. Proc Natl Acad Sci U S A. 2001, 98 (20): 11609-11614. 10.1073/pnas.211424698.PubMedPubMed CentralView ArticleGoogle Scholar
- Kuiken T, Fouchier RA, Schutten M, Rimmelzwaan GF, van Amerongen G, van Riel D, Laman JD, de Jong T, van Doornum G, Lim W, Ling AE, Chan PK, Tam JS, Zambon MC, Gopal R, Drosten C, van der Werf S, Escriou N, Manuguerra JC, Stohr K, Peiris JS, Osterhaus AD: Newly discovered coronavirus as the primary cause of severe acute respiratory syndrome. Lancet. 2003, 362 (9380): 263-270. 10.1016/S0140-6736(03)13967-0.PubMedView ArticleGoogle Scholar
- Peiris JS, Lai ST, Poon LL, Guan Y, Yam LY, Lim W, Nicholls J, Yee WK, Yan WW, Cheung MT, Cheng VC, Chan KH, Tsang DN, Yung RW, Ng TK, Yuen KY: Coronavirus as a possible cause of severe acute respiratory syndrome. Lancet. 2003, 361 (9366): 1319-1325. 10.1016/S0140-6736(03)13077-2.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.