|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
Allelic discrimination is possible by molecular means and multiple concurrent infections are usually found in blood samples from areas endemic for malaria.24 Highly polymorphic marker genes amplified by polymerase chain reaction (PCR) can be used to track individual clones in longitudinal sample sets. Despite the lower detection limit of PCR relative to microscopy, some clones still seem to disappear but are in fact persisting, and genotyping often detects clones recurring periodically at intervals of 48 hours.5,6 We have previously analyzed infection dynamics of individual clones as determined by PCR-restriction fragment length polymorphism (RFLP) using transition models based on the relative frequencies of different patterns of infection defined by presence/absence at successive surveys.7,8 These models take into account the imperfect detectability of clones sequestered at the time of sampling.
We now compare two molecular techniques for identifying individual clones. The first is RFLP of PCR fragments (PCR-RFLP), which has been used in a number of studies to track parasite clones longitudinally. The marker gene merozoite surface protein 2 (msp2) is highly polymorphic because of intragenic repeats, with > 50 genotypes identified by RFLP in each of its two allelic families 3D7 and Fc27.3,9 The pattern of restriction fragments clearly identifies individual clones, and particular patterns are recognized also within mixed infections. However, in samples with more than five concurrent genotypes, the superimposed patterns are increasingly more difficult to analyze. The second technique is based on sizing PCR fragments by an automated sequencer using the Gene-Scan program. This technique uses fluorescent-labeled PCR primers specific for the msp2 allelic families. FC27- and 3D7-type msp2 PCR fragments are identified by the two fluorescent markers 6-Fam and VIC, respectively. The use of a different dye for each allelic family increases resolution of fragment sizing. GeneScan analysis software uses a size standard added to each sample after PCR to create an internal calibration curve to determine the size of each PCR fragment. In comparison with PCR-RFLP, this technique increases throughput and avoids subjectivity in analyzing the readout.
We report analyses of P. falciparum infections in 100 individuals of all ages from the Kassena-Nankana District (KND) of Ghana sampled at 2-monthly intervals over 1 year. We compare the performance of genotyping by PCR-RFLP versus GeneScan for the analysis of the force of infection (clonal acquisition rate) and duration of individual infections based on the patterns of persistence of individual genotypes.
| MATERIALS AND METHODS |
|
|
|---|
DNA isolation and genotyping. Whole blood (10 µL) was dotted onto ISOCode Stix. Irrespective of microscopy results, six samples from each of 100 individuals were screened for presence of P. falciparum by PCR. Processing of Stix, PCR conditions, and RFLP procedures have been described in detail.3,12 In brief, primary and nested PCR were performed for RFLP analysis using primer pair S2/3 (5'-GAAGGTAATTAAAACATTGTC-3'/ 5'-GAGGGATGTTGCTGCTCCACAG-3') and S1/4 (5'-GAGTATAAGGAGAAGTATG-3'/5'-CTAGAACCATGCATATGTCC-3'), respectively, followed by restriction digest with HinfI and DdeI as described elsewhere.3,12
For GeneScan analysis, 2 µL of primary PCR product was amplified in nested PCR with the fluorescent-labeled family-specific primer M5 (FC27-specific: 5'-6-FAM-GCAT-TGCCAGAACTTGAA-3') or N5 (3D7-specific: 5'-VIC-CTGAAGAGGTACTGGTAGA-3') at a concentration of 100 and 200 nM, respectively. The non-fluorescentlabeled forward primer STail (5'-CTTATAATATGAGTATAAGGAGAA-3') was modified at the 5' end by adding a 7-bp tail to avoid non-templatedirected addition of a single nucleotide to the 3' end of a blunt-end double-stranded DNA ("plus-A-artefact").13 The cycle conditions were 5 min at 94°C followed by 30 cycles of 30 s at 94°C, 1 min at 50°C, and 1 min at 70°C, and a final elongation for 7 min at 70°C.
nPCR product (0.5 µL) was combined with 10 µL ROX-labeled size standard (diluted 1:10 with H2O to minimize pipetting errors). Samples were dried and sent to the Genomics Core Laboratory of the MRC Clinical Science Center in London. Highly deionized formamide was added, and after denaturation, samples were analyzed on an ABI PRISM 3700 genetic analyzer.
Determination of detection limits. Genomic DNA was isolated from P. falciparum in vitro cultures (FC27 and 3D7 strains) as described previously.14 DNA concentration was determined photometrically. Tenfold dilutions of genomic DNA were amplified in PCR reactions containing 1000, 100, 50, 10, 5, 1, and 0.1 genomes of either 3D7 or FC27 DNA or a mixture of both.
Data analysis. An in-housegenerated computer program was used to process the output of the GeneScan analyzer. The main tasks of this program were as follows: 1) Cut-off determination: a sample-specific cut-off was used to separate real signals from noise and to allow for variability between GeneScan runs. The cut-off was established by GeneScan analysis of sequenced msp2 reference alleles (cloned msp2 fragments or single clone infections). The presence of a single template per reaction made it possible to determine background fluorescence levels to 300 arbitrary fluorescence units. A sample-specific cut-off was determined by multiplying the arithmetic mean of peak heights of the size standard signals per sample by a constant (the empirically chosen cut-off of 300 units divided by the mean peak height of size standard peaks of all samples). Because with our primers none of the msp2 sequences available at GenBank would give rise to a PCR fragment of < 216 bp, peaks with a measured size of < 200 bp were not considered. 2) Elimination of bleeding and "plus-A-artifacts": spectral overlap of the fluorescent dyes labeling the family-specific primers caused bleeding in case of strong signals caused by PCR fragments present at very high concentrations (> 5,000 fluorescent units). Although the dyes emit light at different wavelengths, some overlap exists despite using a GeneScan software matrix file to remove spectral overlaps. Taq polymerase has terminal deoxynucleotidyl transferase activity to add an extra nucleotide, usually adenine ("plus-A"), at the 3' end of PCR products. This results in two populations of amplified products with a size difference of one nucleotide. Despite the use of a tailed primer that supports "plus-A" addition, we observed in case of high intensity peaks also a small peak about one nucleotide before the actual peak. Our software eliminated peaks caused by bleeding and "plus-A-artifact." 3) Elimination of PCR artifacts: allele-specific PCR artifacts were detectable when particular alleles were present at a very high concentration. Such artifacts are likely caused by the intragenic repeats of msp2 alleles, which can facilitate aberrant annealing of an incompletely synthesized strand to the repeat region. We determined allele-specific PCR artifacts by analyzing cloned msp2 alleles at high DNA concentrations. Artifacts were omitted from subsequent analyses. 4) Genotype calling: For analysis of longitudinal sets of samples, a persisting genotype must be accurately identified in sequential blood samples. There were slight variations among repeated size determinations of identical fragments depending on the concentration of the amplified fragment. To allow for inaccuracies in size determination, peaks were assigned to size bins with a width of 2.4 bp. Because a coding region is genotyped, fragment sizes must differ by multiples of 3 bp.
Statistical analysis. For comparing the two genotyping methods, SAS statistical software Release 8.2 (SAS Institute, Cary, NC) was used. Infection dynamics were analyzed by calculating the frequency of gains, losses, and persistence of infecting clones. An infection present in survey at time t, but not seen in the subsequent survey t + 1 was considered as "loss" (+ ), whereas "gain" ( +) was noted when an infection was observed in round t but not in the previous round t 1. Where infections were observed in consecutive surveys this was recorded as "persistence."
We analyzed the infection dynamics using methods that allow for imperfect detection using the method of Smith et al.7 In further analysis, we fitted an immigration-death model to the full sequences of six observations.8 This provided estimates of the rate of new infection,
, clearance rate, µ, and detectability.
| RESULTS |
|
|
|---|
|
Before application to field samples, the technique was validated on a panel of sequenced reference alleles (data not shown). The SD of repeated sizing of the same fragment was only about 0.06 and 0.17 nucleotides for FC27 and 3D7 DNA, respectively. For the most frequent allele in the Ghanaian field samples, an SD of 0.14 was calculated.
When comparing all msp2 sequences submitted to Gen-Bank until 2004, we detected two genotype pairs of the same allelic family sharing the same amplicon size (accession numbers U07001/U16842 and AY534506/AF010461). Thus, Gene-Scan cannot discriminate these alleles. This represents the limitation of our method.
Longitudinal genotyping in field samples from Ghana. From 100 individuals enrolled at baseline, 550 blood samples were collected during the 1-year follow-up in 2-monthly intervals; 78.2% were found to be positive for P. falciparum by PCR. For 99 individuals, GeneScan analysis was successfully performed, and 96 individuals were analyzed by PCR-RFLP, accounting for 1,405 and 1,084 observed clonal infections, respectively. The discrepancy is caused by differences in the resolution of these methods. Problems occurred in particular with the PCR-RFLP method when interpreting superimposed RFLP patterns.
A total of 164 different msp2 alleles were distinguished by GeneScan analysis; 116 and 48 belonged to the 3D7 and FC27 allelic family, respectively. Frequencies of 3D7-type alleles were all < 3% (Figure 1a
). Some FC27-alleles occurred at very high allelic frequencies, the most frequent one reaching 14% (Figure 1b
).
|
|
|
|
|
, and death refers to the clearance of a parasitic genotype with rate µ. The model assumes that corresponding to each observed process that is the detection or failure to detect a parasitic genotype is an underlying true process, which is hidden as a result of imperfect detection by the parasitological diagnostic tool (GeneScan or RFLP). For GeneScan, this approach estimates a detectability of 0.35 (95% CI: 0.310.39) that is an average of 35% of the parasite present in the host are detected in a finger-prick blood sample. For RFLP, the detectability was 0.47 (0.420.51). This estimation of the detectability was also done based on a simpler approach,7 which led to comparable results (data not shown).
Based on the genotype acquisition and clearance rates
and µ, respectively, we determined the average duration of an infection. We estimated a value of 19.6 (17.322.0) per year for
and 1.7 (1.42.1) per year for µ, which corresponds to an average duration of infection of 210 (176256) days when samples were analyzed by GeneScan. For RFLP analysis,
and µ were found to be 16.3 (14.818.0) and 2.4 (2.12.7), respectively, which corresponds to an average duration of infection of 152 (133177) days. Our comparison of both genotyping techniques and both statistical approaches showed that the turnover rate for 3D7- type alleles was almost twice as high as for the FC27-type alleles and that infections with FC27-type alleles persisted longer.
| DISCUSSION |
|
|
|---|
In our study, we investigated several parameters: MOI by age, different transition states by allelic family, and seasonal effects on these transitions. In both data sets, mean MOI was low among infants, increased steadily until 59 years, and started to decrease again during adulthood. This finding is consistent with the previously published cross-sectional results from the first survey.11 The distribution of the three transition types (gain, loss, and persistence of clones) was similar in both methods, with acquisition of infections being most frequent within the first and second transition interval, reflecting the increase in transmission rate during the second half of the wet season. The number of transitions, however, was different between the two laboratory methods. By Gene-Scan, more genotypes were gained, whereas by PCR-RFLP, the clearance rate was higher. Both these findings might be caused by higher sensitivity of GeneScan. In particular, 3D7-type alleles are difficult to determine in polyacrylamide gels used for PCR-RFLP but easy to identify by using family-specific fluorescent primers in GeneScan analysis, resulting in a higher number of acquisitions. The higher percentage of losses by PCR-RFLP can also be explained by higher sensitivity of GeneScan, because the transition pattern + by RFLP analysis was found to be + + by GeneScan.
Another parameter calculated from transition rates was the duration of infections. The turnover rate of 3D7 infections was higher in both data sets. This suggests that FC27 alleles are more stable over time and more resistant to elimination under selective pressure by the hosts immune system than 3D7 alleles.
Sensitivity and detectability. The sensitivity of PCR-RFLP versus the GeneScan technique and their performance in a molecular epidemiologic field study can be compared by the total number of parasite clones detected by each technique. Compared with PCR-RFLP, we found an increase by 22.8% for GeneScan-based genotyping. This is also reflected by the parameter "mean MOI," which by GeneScan showed an increase of 1.6 infections at the peak of the age distribution. We found that the increased sensitivity in field samples was caused by higher precision in discriminating 3D7-type genotypes. This can be explained by the scrambled repeat structure of 3D7-type alleles, giving rise to only small differences in fragment sizes that can hardly be discriminated by PCR-RFLP in gel electrophoresis. We conclude that the increase in detected clones by GeneScan typing was accounted for by problems of PCR-RFLP in resolving the complex 3D7-type patterns in multiple infections.
Even optimal sensitivity of P. falciparum detection does not reflect the whole parasite population in an infected individual. Sequestration of late stage parasites causes a 48-hour periodicity in detectability. Such fluctuations were observed when monitoring the daily dynamics of P. falciparum clones.5 Thus, imperfect detection of some of the P. falciparum clones concurrently present in a host is a consequence of the parasites life cycle, and the estimation of persistence is complicated even using molecular methods. Therefore, mathematical models need to be applied to estimate the molecular parameters based on such imperfect data.
Both the approaches we used7,8 take into account the frequencies of transitions. A substantial frequency of the pattern {+ +} indicates imperfect detection. The detectability is our estimate of that proportion of the duration of an infection during which the densities are high enough to be detected, averaged over all the infections known to be present. Paradoxically, we found a lower detectability by GeneScan (35%) than by PCR-RFLP (47%), despite the fact that GeneScan had detected more genotypes. We arrived at the same conclusion when the method of Smith et al.7 was used. This can be explained by different detection limits. Figure 5
shows the density profile of three hypothetical infections. Curve A represents an infection that persists at detectable density and is detected in almost the same proportion by both techniques. Curve B represents an infection that occasionally reaches detectable density, but persists sub-patently for a long time. A proportion of such infections is detected by GeneScan, but not by RFLP. Curve C represents an infection that persists for a short time and is also more likely to be detected by Gene-Scan. It is likely that patterns B and C contribute to a reduced overall estimate of GeneScan detectability compared with RFLP, which never sees infections such as B or C.
|
The potential of GeneScan-based msp2 genotyping has recently been shown for the discrimination of recrudescent from new infections in drug efficacy studies.22 We anticipate further applications of molecular parameters in a range of intervention studies against malaria.23 In malaria vaccine trials, molecular monitoring has been applied to detect selective effects if the vaccine was polymorphic.17,24 Molecular infection dynamics parameters could further contribute to describe the parasitological outcomes of vaccine trials and perhaps elucidate biologic effects of candidate vaccines.
Received July 11, 2005. Accepted for publication December 8, 2005.
Acknowledgments: We thank the community members of the KND, especially the participants and the parents of the children consenting on their behalf, the staff of the Navrongo Health Research Center for their field assistance, especially Lucas Amenga-Etego and Victor Asoala, and Béatrice Glinz-Szára and André Tiaden for genotyping.
Financial support: This work was supported by Swiss National Science Foundation Grant 3300C0-105994).
* Address correspondence to Ingrid Felger, Swiss Tropical Institute, Socinstrasse 57, PO Box, CH-4002 Basel, Switzerland. E-mail: ingrid.felger{at}unibas.ch ![]()
Authors addresses: Nicole Falk, Nicolas Maire, Wilson Sama, Tom Smith, Hans-Peter Beck, and Ingrid Felger, Swiss Tropical Institute, Socinstrasse 57, PO Box, CH-4002 Basel, Switzerland, E-mails: nicole.falk{at}stud.unibas.ch, nicolas.maire{at}unibas.ch, Wilson.Sama{at}unibas.ch, Thomas-A.Smith{at}unibas.ch, hans-peter.beck{at}unibas.ch, and ingrid.felger{at}unibas.ch. Seth Owusu-Agyei, Kintampo Health Research Centre, Ghana Health Service, PO Box 200, Kintampo, Ghana, E-mail: seth.owusu-agyei{at}ghana-khrc.org.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
B. Greenhouse, C. Dokomajilar, A. Hubbard, P. J. Rosenthal, and G. Dorsey Impact of Transmission Intensity on the Accuracy of Genotyping To Distinguish Recrudescence from New Infection in Antimalarial Clinical Trials Antimicrob. Agents Chemother., September 1, 2007; 51(9): 3096 - 3103. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |