| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, February 2007, p. 521-527, Vol. 51, No. 2
0066-4804/07/$08.00+0 doi:10.1128/AAC.00902-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Department of Epidemiology, University of North Carolina, Chapel Hill, North Carolina
Received 21 July 2006/ Returned for modification 25 September 2006/ Accepted 7 November 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Using these molecular tools, a higher COI has been associated with symptomatic malaria manifestation in pregnant women (3). In a recent study, the risk of treatment failure in children was positively correlated with the COI (13).
The current World Health Organization protocol recommends that all clinical efficacy studies with more than 14 days of follow-up use PCR-based methods to determine whether recurrent infection is due to drug failure (recrudescence) or a newly acquired infection (reinfection) (25). The inherent limitations of conventional nested PCR genotyping, such as sensitivity only to size differences, and different interpretations of the PCR comparisons can both have a dramatic impact on the treatment failure rate (18). Therefore, more definitive methods are needed.
We have previously used the heteroduplex tracking assay (HTA) in an area of low malaria transmission to characterize the COI in Southeast Asia (15, 26). In the HTA, radiolabeled P. falciparum merozoite surface protein 1 (PfMSP1) probe is annealed to amplicons from a patient whose PfMSP1 block-2 region has undergone a single round of PCR amplification. Double-stranded complexes are formed between the probe and host amplicon(s); these complexes, termed heteroduplexes, migrate on a nondenaturing gel at various speeds based on the complementarity of the patient-probe complexes. Clustered base-pair mismatches, insertions, and deletions all change the migration speed of the heteroduplexes and, after electrophoresis on a nondenaturing gel, the host PfMSP1 sequence variants can be quantified. Because they are simple to perform, sensitive to size and clustered sequence polymorphisms, and allow quantification of each sequence variant, HTAs are commonly used by virologists to characterize polyclonal human immunodeficiency virus (HIV) infections (5, 10). We report here the first use of the HTA to characterize P. falciparum sequence diversity in an area of high malaria transmission and discuss a new approach to determine the rate of treatment failure.
(These results were presented in part at the 54th Annual American Society for Tropical Medicine and Hygiene Conference, Washington, DC, 2005 [abstr. 1007].)
| MATERIALS AND METHODS |
|---|
|
|
|---|
Malaria was diagnosed by microscopic inspection of Giemsa-stained thick blood films. The total number of parasites detected per 200 leukocytes was multiplied by 30 to estimate the parasite count per microliter of blood (assuming a leukocyte count of 6,000 leukocytes/ml of blood). A thick film was considered negative if no parasites were detected after examining 100 microscopic fields. Parasite density was dichotomized at the median value (930 parasites/µl) into categories of high and low parasite density.
Women were offered counseling and testing for HIV and nevirapine according to the HIVNET-012 protocol if they were HIV positive. HIV testing was performed by using the Determine HIV-1/2 rapid test (Abbott Laboratories) and Unigold Test (Trinity Biotech, Dublin, Ireland). Discordant results were resolved with the Hemastrip Rapid Test (Saliva Diagnostic Systems).
The hemoglobin concentration was measured by using a HemoCue machine (HemoCue, Inc., Angelholm, Sweden). Women were considered anemic if their hemoglobin concentration was <11 g/dl.
DNA isolation. Venous blood was spotted onto filter paper and stored under desiccation. DNA was extracted from the filter paper by using the QIAMP DNA minikit (QIAGEN, Valencia, CA) according to the manufacturer's instructions.
PCR. Block-2 of PfMSP1 (22) was amplified with a Peltier thermal cycler (MJ Research, Waltham, MA) according to the protocol of Ngrenngarmlert et al. (15) using the primers C1F (5'-GAAGATGCAGTATTGACAGG-3') and C3R (5'-TGATTGGTTAAATCAAAGAG-3') (8); the C1F primer starts at nucleotide 106 in the 3D7 MSP1 gene (9). PCR products were run on a 1% agarose gel, stained with ethidium bromide, visualized by UV light, and sized against a 100-bp molecular weight marker (Promega, Madison, WI). Successful amplification resulted in a product of approximately 400 bp. All PCRs were done in duplicate.
Plasmids and probes.
The COI was determined using a single, heterologously expressed PfMSP1 block-2 probe (MHP [described in reference 15]). Ten micrograms of the MHP-pT7Blue construct was digested with BamHI at 37°C for 1 h and then end labeled in the same buffer by filling in the BamHI overhang for15 min at 25°C (1x BamHI buffer [New England Biolabs], 10 mM dithiothreitol, 50 µM dGTP, 50 µCi of [
-35S]dATP [1,250 Ci/mmol; ICN], 10 U of Klenow fragment). The reaction was stopped by the addition of EDTA (1 mM final concentration) and heat inactivation at 75°C for 15 min. The radiolabeled probe was released from the vector by digestion with PstI at 37°C for 1 h, separated on an agarose gel, excised, and gel purified with the QIAquick gel extraction kit (QIAGEN).
Heteroduplex tracking assay. A total of 8 µl of PCR-amplified block-2 of PfMSP1 from patient samples was mixed with 1.2 µl of 10x annealing buffer (1 M NaCl, 100 mM Tris-HCl [pH 7.5], 20 mM EDTA), 0.1 µM C1F primer, 0.1 µM C3R primer, 1 µl of 35S-labeled probe, and 2 µl of 6x loading dye (15% Ficoll, 0.25% bromophenol blue, 0.25% xylene xylanol [Sigma, St. Louis, MO]) in a total volume of 12 µl. The mixture was denatured at 95°C for 2 min and allowed to reanneal at 25°C for 5 min. Heteroduplexes were electrophoresed on a 6% polyacrylamide gel (acrylamide-bisacrylamide, 37.5:1) in 1x Tris-borate-EDTA buffer at 17 mA per gel for 5.5 h. The gels were dried onto filter paper (Whatman) and exposed to BioMax MR X-ray film (Kodak, Rochester, NY). For quantitative data, the dried gels were exposed to a phosphorimager screen for 7 days, and the band intensities were determined by quantifying the area under the curve according to the manufacturer's instructions (Molecular Dynamics). HTA bands were considered unique PfMSP1 sequence variants if they were not in the probe alone lane and they were present in both PCR replicates. The relative migration distance (Rf) was calculated for each variant by dividing the distance migrated by the variant by the distance migrated by the probe homoduplex. To test the precision of the Rf measurement, a single variant was run on 16 separate gels and the Rf was calculated. The cloned variant had a coefficient of variation of 14% (the average Rf was 0.16, and the standard deviation was 0.022).
Colony PCR. The PfMSP1 amplicon from participant 1/019 was cloned into the pCR2.1TOPO by using the Invitrogen TOPO TA cloning kit and transformed into Escherichia coli (Invitrogen Corp., Carlsbad, CA). Then, 52 colonies were selected and PCR amplified using the primers and thermal cycler settings described above. PCR products from representative variants were cleaned by using the QIAquick PCR purification kit (QIAGEN) and sequenced using the C1F primer at the University of North Carolina-Chapel Hill Genome Analysis Facility. DNA sequences were hand cleaned, aligned with CLUSTAL X (http://www2.ebi.ac.uk/clustalw), and displayed with GeneDoc (http://www.psc.edu/biomed/genedoc/).
Statistics and modeling. The association between the COI and clinical features was tested by using the Wilcoxon rank-sum test. The relationship between the COI and treatment failure was monotonic (data not shown), and therefore COI was coded as a continuous variable in the multivariable regression. Primigravidity and parasite burden (dichotomized at the median) were used as covariates in the multivariable regression. All statistical analysis was done by using STATA v.8.2, (STATACorp, College Station, TX).
Ethical approval and consent. This study was approved by the College of Medicine Research Committee (University of Malawi, Blantyre, Malawi) and the Institutional Review Boards of the University of North Carolina (Chapel Hill, NC). Separate consent forms in the local language were administered for blood and/or data collection and for HIV testing.
GenBank accession numbers. Six MSP1 block-2 sequences from participant 1/019 are deposited in GenBank under accession numbers DQ855130 to DQ855135. The MHP probe is deposited under accession number DQ858166.
| RESULTS |
|---|
|
|
|---|
PfMSP1 sequence variation. HTAs were run, in duplicate, on the 125 patients, and the numbers of variants were determined. A representative gel is shown in Fig. 1; sequence variants migrated both above and below the single-stranded probe. The participants had a mean COI of 3.82 (median = 3, range = 1 to 11). The distance migrated by each variant was measured, and the variants were divided into 20 bins of 0.05 Rf units. The prevalence of individual variants ranged from 0.03 to 15.75%, and three variants were present in >10% of the women (Fig. 2).
|
|
|
|
Classification of recurrences. One important potential use of HTA would be to distinguish new infections (reinfections) from treatment failures (recrudescences). Many reports use PCR to distinguish reinfection from treatment failure but do not specify the criteria used to distinguish each outcome. Often, if the initial and recurrent infections share any variant, then the recurrent infection is classified as a treatment failure (12, 17). However, in an area of high malaria transmission, such as Malawi, if specific PfMSP1 variants are frequent in the parasite population, then a new infection containing one of these variants might be misclassified as a failure. In the present study, three variants were common, with prevalences between 11 and 16% (Fig. 2). Since the average infection comprised four variants, any recurrent infection had a 37 to 50% chance of containing each of these variants coincidentally. To eliminate the possibility of misclassifying patients with these common variants as failures, we consider a second episode indeterminate if the initial and recurrent episodes only shared one band with a prevalence >10%.
Comparison of first and second malaria episodes. Table 2 summarizes the HTA patterns of 35 women with second malaria episodes. Five patients failed treatment within 14 days and were considered early treatment failures. Thirty patients had a second infection 14 to 56 days posttreatment. Six of these shared multiple and/or uncommon PfMSP1 variants with the first episode with no new variants and were considered late treatment failures. Four shared no variants with the first episode and were considered reinfections. The remaining 20 women had recurrent infections that both shared variants with the initial infection and contained new ones. For 10 of these participants, either the shared variant was uncommon (<10% prevalence) or two common variants were shared, and the second episode was classified as a late treatment failure. For the remaining 10 patients, the single shared variant was common (>10% prevalence), and the second episode was considered indeterminate.
|
Estimation of true failure rate. Thus, there were 10 second episodes and 4 third episodes that were classified as indeterminate. Fortunately, since we know the prevalence of the common bands and the number of variants in the recurrent episode, we can calculate a probability of reinfection by each variant. For the recurrent patient with x variants and sharing a single variant of prevalence y, the binomial probability (16) that this variant is found by chance in a recurrent infection is calculated as 1 (1 y)x [and therefore the probability that it represents a failure is (1 y)x]. The probabilities of reinfection for the 14 indeterminate episodes are shown in Table 3. The mean probabilities of failure are 0.59 for the 10 indeterminate second episodes (or a predicted 5.93 patients) and 0.48 for the 4 indeterminate third episodes (or a predicted 1.91 patients). Thus, the estimated failure rates are (21 + 5.93)/125 or 21.5% for the first treatment and (4 + 1.91)/35 or 16.9% for the second treatment.
|
| DISCUSSION |
|---|
|
|
|---|
The HTA method has several advantages. First, it is more sensitive to minority variants than nested PCR and is less prone to cross-contamination (15). Second, as we demonstrate here, it is quantitative and capable of measuring variants representing as little as 1% of the parasite population in an individual host. Third, because it generates highly reproducible Rf values, one can determine the prevalence of individual variants in a given population. Fourth, the combination of HTA and colony PCR allows for the rapid cloning and sequencing of variants.
The major disadvantage of the HTA is that it requires the use of radioactive probes. However, once HTA methods are firmly established, related nonradioactive methods, such as the heteroduplex migration assay (HMA) can be developed (4, 24); whether the HMA can detect low-abundance variants as well as the HTA remains unknown. In addition, the HTA is also potentially limited in its resolution, in that multiple variants, with similar sequences, could run at the same position on the gel; these variants could only be detected by cloning and sequencing. A third disadvantage of the HTA is the possibility that observed variants are actually recombination artifacts. Because the probability of recombination increases with the size of the amplicon, the degree of recombination of a 400-bp PfMSP1 block-2 amplicon would most likely be less than the recombination observed by Tanabe et al. (23). However, this drawback holds for all PCR-based methods.
HTAs, combined with a binomial probabilistic approach, represent a new way to estimate drug efficacy. Alternative categorical methods for "PCR correction" suffer from several problems. First, since non-HTA methods may not be able to detect minority variants, sharing of such variants by enrollment and recurrent samples might be missed and failures misclassified as reinfections. Second, since specific variants may predominate in certain areas, reinfections, by these common variants, might be misclassified as failures. In contrast, HTAs can detect minority variants. Thus, inferring dichotomous outcomes (failure versus reinfection) for PCR correction methods can be misleading. However, by calculating the prevalence of each variant in a population, we can determine the probability that each individual patient originally scored as indeterminate failed and, by summation, the estimated overall failure rate.
There are several limitations to this method. First, if a treated patient were reinfected but bore gametocytes from his initial infection, then this patient might be misclassified as a failure (25). This is a problem facing all genotyping correction methods. A second potential limitation of the binomial method, which considers the prevalence of the recurrent variants, is the assumption that the variants are transmitted independently (i.e., a single mosquito bite only contains a single parasite variant). Although this assumption has not been directly assessed in Malawi, in 2003 in the Republic of Guinea Bissau Arez et al. (1) found single PfMSP1/MSP2/GLURP alleles in ca. 50% of the P. falciparum-infected mosquitoes. A third potential limitation is the choice of bin size. Although it is possible that different variants may be counted within the same bin and smaller bin sizes might be preferable, in our study we observed only one individual out of 125 for whom this may have happened. A fourth limitation is that this method does not lend itself to survival analysis, a method commonly used in drug efficacy studies.
In summary, by HTAs the COI was found to be elevated in Malawian women with HIV and in those who failed therapy. Assigning a probability of failure to specific recurrent episodes mitigates misclassification for clinical trials in Africa and other high transmission areas. The method, developed here, will now be used to compare the efficacy of SP, SP plus artesunate, and SP plus azithromycin within this study population.
| ACKNOWLEDGMENTS |
|---|
This research was supported by CDC/ASPH grant S1935-21/21.
| FOOTNOTES |
|---|
Published ahead of print on 20 November 2006. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Clin. Vaccine Immunol. | Clin. Microbiol. Rev. |
|---|---|
| J. Clin. Microbiol. | ALL ASM JOURNALS |