Previous Article | Next Article ![]()
Journal of Clinical Microbiology, July 2008, p. 2167-2174, Vol. 46, No. 7
0095-1137/08/$08.00+0 doi:10.1128/JCM.00081-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention, WHO Collaborating Center for Malaria, Schistosomiasis, and Filariasis, 207 Rui Jin Er Road, Shanghai 200025, China,1 National Engineering Center for Biochip at Shanghai, 151 Li Bing Road, Pudong, Shanghai 201203, China2
Received 31 December 2007/ Returned for modification 14 March 2008/ Accepted 17 April 2008
|
|
|---|
|
|
|---|
Antimalarial drug resistance has been associated with single-nucleotide polymorphisms (SNPs), in particular P. falciparum genes. The P. falciparum CQ resistance transporter (Pfcrt) gene K76T mutation has been linked to CQ and possibly AQ resistance (6, 17, 20); the P. falciparum multidrug resistance 1 (Pfmdr1) gene N86Y mutation has been associated with the response of the parasite to amino quinolines such as MQ (9, 25, 27) and LUM (26); mutations N51I, C59R, S108N/T, and I164L in the P. falciparum dihydrofolate reductase (Pfdhfr) gene have been associated with resistance to pyrimethamine (13, 14, 19, 21, 22); mutations S436A, A437G, K540E, A581G, and A613 S/T in the P. falciparum dihydropteroate synthetase (Pfdhps) gene have been shown to confer resistance to sulfadoxine (13, 14, 21, 22); and P. falciparum ATPase-6 (PfATPase6) S769N mutation has been associated with increased 50% inhibitory concentrations of artemether (11, 16, 18). These SNPs are believed to represent molecular epidemiology surveillance tools of antimalarial drug resistance, which can complement the more conventional and logistically complex in vitro or in vivo test (7, 24, 33, 34).
Existing molecular methods for analysis of these SNPs include allele-specific PCR, PCR-restriction fragment length polymorphism (RFLP) analysis, multiplex PCR (mPCR)-RFLP, DNA sequencing, dot blot hybridization techniques, the molecular beacon method, real-time PCR, PCR-enzyme-linked immunosorbent assay, and pyrosequencing (1, 2, 6, 8, 10, 28, 31, 37, 38). Each of these techniques offers its own advantages and limitations. Among the techniques, PCR-RFLP protocols constitute the most common methodological approach for the analyses of these SNPs, but this approach is laborious and time-consuming. A rapid and high-throughput genotyping method would be ideal for large-scale population-based studies.
The oligonucleotide microarray technology, also known as DNA chip, has provided an incredible technical development in rapid, large-scale, and automatic analysis of SNPs (5, 36). Briefly, the principle of oligonucleotide microarray is based on reverse Southern hybridization on a glass or silicon substrate. The fluorescently labeled products hybridized to the immobilized probes on the substrate can be easily detected by using a fluorescence-detecting scanner, and the readable fluorescent signals obtained are combined to determine the genotype of genes (12). In recent years, a microarray technique has been successfully applied to the detection of resistance genes or resistance-associated SNPs in various pathogens (3, 4, 23, 29).
Here, we introduce a one-step mPCR in combination with an oligonucleotide microarray method for rapid, high-throughput detection of the antimalarial resistance-associated SNPs in Pfcrt, Pfmdr1, Pfdhps, Pfdhfr, and PfATPase6.
|
|
|---|
DNA extraction. Genomic DNAs of the P. falciparum laboratory strains 3D7, HB3, and Dd2 were provided by the Malaria Research and Reference Reagent Resource Center (http://www.mr4.org/). Extraction of DNA from bloodspots on filter paper was carried out by the Chelex-100 (Bio-Rad Laboratories, Hercules, CA) method described by Wooden et al. (32) with some modifications described by Pearce et al. (22). The quality of the DNA samples was tested by determining the optical density at 260 and 280 nm.
SNPs studied. The following 21 SNPs were analyzed: Pfcrt 391T/A, 392G/C, 399G/T, 400A/G, 402T/A, and 404A/C (corresponding to codons C72S, M74I, N75E, and K76T); Pfmdr1 256A/T and 257A/T (N86Y/F); Pfdhps 1482 T/G, 1483C/T/G, 1486 C/G, 1794 A/G, 1918 C/G, and 2013G/T/A (S436A/F/C, A437G, K540E, A581G, and A613S/T); Pfdhfr 148T/C, 152A/T, 153T/C, 175T/C, 323G/A/C, and 490A/T (C50R, N51I, C59R, S108N/T, and I164L); and PfATPase6 2306G/A (S769N). These SNPs are located at 11 nucleotide positions since some positions include two or more SNPs.
mPCR. To reduce the time and cost of DNA amplification, an mPCR was developed to amplify the DNA fragments containing all of the SNPs studied in one reaction. Five pairs of primers were designed by using Primer 5.0 software and checked for specificity in a BLAST search available through the National Center for Biotechnology Information (NCBI; http://www.ncbi.nlm.nih.gov/BLAST/). The primers used in the mPCR are shown in Table 1. mPCR was performed with the PCR amplification kit of Genscript (Piscataway, NJ) according to the manufacturer's instructions. Briefly, the mPCR was performed by using the PTC 200 cycler (MJ Research, Inc., Waltham, MA) in a total volume of 50 µl of reaction mix containing 2.5 µl of template DNA, 25 µl of 2x Genscript master mix, and 1.0 µl of primer mixture (consisting of 3.2 µl of 12.5 µM P1F, 3.2 µl of 12.5 µM P1R, 2.0 µl of 12.5 µM P2F, 2.0 µl of 12.5 µM P2R, 2.0 µl of 12.5 µM P3F, 2.0 µl of 12.5 µM P3R, 2.0 µl of 12.5 µM P4F, 2.0 µl of 12.5 µM P4R, 2.0 µl of 12.5 µM P5F, and 2.0 µl of 12.5 µM P5R). The mPCR cycling conditions were 94°C for 15 min for 1 cycle; followed by 30 cycles of denaturation at 94°C for 40 s, annealing at 50°C for 2 min, and extension (ramping from 50 to 72°C by 0.1°C/s) for about 3 min 40 s; followed in turn by 72°C for 3 min for 1 cycle and holding at 4°C. The resulting amplicon covers all of the 21 observed SNPs. A total of 3 µl of PCR product was identified by electrophoresis on a 20-g/liter agarose gel to confirm the successful amplification.
|
View this table: [in a new window] |
TABLE 1. Primers used in the mPCR
|
Oligonucleotide probe design and microarray manufacture. For each of the polymorphic sites associated with antimalarial drug resistance, a group of probes was designed, differentiated by the interrogated base(s) positioned in the central region of the probes. To reduce the number of potential ambiguous hybridization results, some oligonucleotide probes were designed to represent combinations of SNPs rather than one SNP. In all, 11 groups of probes were designed. The specificity of the oligonucleotide probes was verified through a BLAST search. Based on initial hybridization experiments, some probes were optimized by modifying their length and sequence to achieve more uniform fluorescence intensities across the microarray. Probes were commercially synthesized and modified (Table 2). Each probe contained a 5'-amino group for immobilization chemistry and a 15-mer poly(dT) spacer, followed by the nucleotide hybridization sequence.
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide probes used in the microarray
|
![]() View larger version (67K): [in a new window] |
FIG. 1. Schematic diagram of the probe positions on the microarray. A1 to K5, specific probes for detection of the SNPs; S, spotting control probe; N, negative control probe; P, positive control probe; B, blank control probe. Each probe was spotted three times.
|
Microarray scanning and analysis. Slides were scanned by using a scanner GenePix 4000B (Axon Instruments), and data acquisition and processing were performed by using GenePix Pro 3.0 software. The median of local background-corrected feature intensities (henceforth referred to as net intensities) and their signal-to-noise ratios (net intensities/standard deviation of the local background) were used for further analyses. Each signal was ranked from highest to lowest in each probe group according to their net intensities. The highest ranked signal in each group was regarded as a positive signal, and a signal with an intensity more than one-third of the strongest signal in that group was also considered a positive signal, whereas signals with intensities less than one-third of the strongest signal were considered negative. When hybridized to single-infection samples, only one probe in each group produced a positive signal as the perfect match. In multiple infections, two or more probes with the highest signal intensities in a group might be selected as positive.
Microarray analyses of synthesized oligonucleotides. Because genotype variations of the five genes from the samples studied here were limited, we also used a set of 37 synthesized 50-bp oligonucleotides (Table 3) harboring all known genotypes of the loci studied to evaluate the microarray. Twenty combinations (Table 4) of the synthesized oligonucleotides were detected by the microarray to test whether the method can correctly identify all possible genotypes in single or mixed infection. The labeling and hybridization conditions were as described above.
|
View this table: [in a new window] |
TABLE 3. Thirty-seven synthesized oligonucleotides
|
|
View this table: [in a new window] |
TABLE 4. Twenty combinations of the synthesized oligonucleotides
|
|
|
|---|
![]() View larger version (30K): [in a new window] |
FIG. 2. Agarose gel electrophoresis of mPCR products. M, 100-bp DNA ladder. Lane 1, P. falciparum 3D7; lane 2, P. falciparum Dd2; lane 3, P. falciparum HB3; lanes 4 to 5, P. falciparum samples collected from field; lane 6, P. vivax; lane 7, P. berghei; lane 8, P. cynomolgi; lane 9, blank control (H2O).
|
![]() View larger version (12K): [in a new window] |
FIG. 3. Signal-to-noise ratios and mismatch discrimination. The data are based on three independent hybridizations of P. falciparum 3D7. (A) Signal-to-noise ratios (means and standard errors) of perfectly matching probes; (B) relative intensities of mismatched probes (means and standard errors).
|
![]() View larger version (28K): [in a new window] |
FIG. 4. Scanning images of P. falciparum 3D7 and Dd2. For 3D7, the positive probes were A1, B1, C1, D1, E1, F1, G4, H1, I1, J1, K1, P, and S (genotypes: Pfdhfr, 50/51CN, 59C, 108S, and 164I; Pfmdr1, 86N; PfATPase6, 769S; Pfdhps, 436/437 SG, 540K, 581A, and 613A; Pfcrt, 72-76CVMNK). For Dd2, the positive probes were A2, B2, C2, D1, E3, F1, G6, H1, I1, J2, K2, P, and S (genotypes: Pfdhfr, 50/51CI, 59R, 108N, and 164I; Pfmdr1, 86F; PfATPase6, 769S; Pfdhps, 436/437 FG, 540K, 581A, and 613S; Pfcrt, 72-76CVIET).
|
The specificity of microarray assay was tested by using mPCR products from P. berghei, P. cynomolgi, and P. vivax. For these three parasites no positive signal was obtained with any specific probe except the spotting control probe and the positive hybridization control probe.
To determine at what level minor SNPs in mixed infections can be detected, we mixed the genomic DNA of the 3D7 and Dd2 strains. These two strains have different genotypes at Pfdhfr codons 51, 59, and 108, Pfmdr1 codon 86, Pfdhps codons 436 and 613, and Pfcrt codons 74, 75, and 76. Although the DNA concentration of one parasite strain was fixed at 3 ng/µl, the DNA concentration of the other strain was serially 10-fold diluted from 30 ng/µl to 0.3 pg/µl. When the DNA concentration of 3D7 was fixed, Dd2 could be satisfactorily analyzed at 30 pg/µl. When the DNA concentration of Dd2 was fixed, 3D7 could also be satisfactorily analyzed at 30 pg/µl.
There was 98.2% (926/943) agreement among the genotypes obtained from three identical subarrays in the same slide for 87 samples. Thirty samples were investigated in two separate tests by the microarray method and interpreted blindly, and the reproducibility between the repeated typing was 96.1% (317/330).
Microarray analyses of the synthesized oligonucleotides. All SNPs in the 20 combinations of the synthesized oligonucleotides were identified without ambiguity by the microarray except for two genotyped positions of Pfdhps 540 that could not give usable signals due to spotting failure. The net intensity ratios of matching probes to the mismatch ones were more than five times for 83.9% of the total and between 2.7 and 5.0 times for the others. Figure 5 shows a scanning image of the combination of 18 of the synthesized oligonucleotides (simulation of mixed infection of P. falciparum 3D7 and P. falciparum Dd2).
![]() View larger version (52K): [in a new window] |
FIG. 5. Scanning image of simulated mixed infection of P. falciparum 3D7 and P. falciparum Dd2 (combination of 18 synthesized oligonucleotides). The positive probes were A1, B1, C1, D1, E1, F1, G4, H1, I1, J1, K1, A2, B2, C2, E3, G6, J2, K2, P, and S.
|
|
|
|---|
A high degree of concordance was observed when we compared the new microarray assay with the existing nested PCR-sequencing method, probably the most credible method currently used for detection of these SNPs. This indicates that the microarray technique for drug resistance-associated genes genotyping is reliable and accurate.
Traditional nested PCR method requires 10 PCRs for application of the five genes, which may take five workdays to accomplish. Here, the use of mPCR to amplify the five genes in one reaction leads to a significant decrease in time, cost, and number of manipulations. It takes a total of about 8 h to perform the entire detection procedure, from mPCR reaction (4.0 h), DNA labeling (2.0 h), hybridization (1.5 h), and microarray washing (0.25 h) to microarray scanning and analysis (0.25 h). There are 30 subarrays in a slide, allowing the genotyping of 10 different samples thrice by each assay, so approximately 100 samples could be analyzed for the 21 resistance-associated SNPs in one workday by two technicians. Therefore, this method is especially suitable for large-scale surveillance of the molecular markers. Moreover, the method requires less genomic DNA, which is propitious for analyzing samples with low parasite density.
The high throughput of the microarray method substantially reduces the cost. In estimating the cost, we included expenses associated with DNA isolation, multiplex-PCR, PCR product labeling, microarray production, microarray hybridization, and microarray washing. The equipment and labor costs were not included. We calculated a price of 2.9 Chinese yuan ($0.41 [United States]) per SNP for the microarray method; this cost was lower than that for pyrosequencing ($2.28 per SNP), conventional sequencing ($3.66 per SNP), and RFLP ($6.58 per SNP) (37).
Filter paper samples used in the present study offered considerable advantages for field collection, transportation, and storage over frozen liquid samples, but eight samples collected from Hainan province in 2001 failed to be analyzed by microarray. We found that these samples were not stored in a dry environment and became moldy before tests were conducted. Thus, it is probable that the genomic DNA of the parasite in these samples was destroyed or degraded. Therefore, it is important to store samples properly to ensure high-quality DNA. In addition, since false microarray determinations were mainly caused by nonspecific signals and poor probe spotting, the microarray silylated slides should be carefully kept from dust contamination and more efforts should be made to further improve the spotting quality.
The testing results of the 20 combinations of the synthesized oligonucleotides show that the microarray is able to identify all of the possible genotypes at the loci studied. Furthermore, it indicates that the technique can also differentiate the SNPs of mixed infections with P. falciparum strains of different genotypes, which allows the application of microarrays in various regions for different genetic background of P. falciparum.
As a limitation, the microarray could only be used to analyze the known SNPs in the resistance-related genes. Further studies are needed to develop a microarray assay for detection of the Pfmdr1 copy number, which has been associated with resistance against chloroquine, mefloquine, and probably artemisinins. Since new SNPs associated with antimalarial drug resistance continue to be discovered, the present microarray should be regularly updated by adding new probes.
It should be noted that microarray technology requires a microarray spotter and a microarray scanner, which may not be available in ordinary laboratories. We hereby suggest that a microarray laboratory be established in one or several countries where P. falciparum is endemic so that every site within the area can send filter paper samples to that laboratory for testing. Moreover, in the next phase, we will try to adopt chemiluminescence in the microarray technology so that the microarray assay can be performed in more laboratories.
In conclusion, mPCR-based microarray provides a promising tool for resistance-associated SNP detection in P. falciparum. This approach not only saves time and cost but also permits high-throughput detection of the SNPs. However, we consider that molecular techniques, including microarray, will complement rather than replace conventional in vivo or in vitro testing. Conventional resistance testing will remain indispensable for the detection of resistance because the molecular results do not always correlate with the clinical response associated with P. falciparum. Therefore, microarray-based analysis will be most useful in combination with conventional methods. We expect that the method developed here will considerably facilitate the molecular surveillance of antimalarial resistance in countries and regions where malaria is endemic.
We are grateful to Yi Zhong Duan and Ying Xue Lin for providing P. falciparum samples. We thank Yao Yu Feng and Yi Chang Ni for critically reading the manuscript. We also thank the Malaria Research and Reference Reagent Resource Center for providing the genomic DNA for P. falciparum 3D7, HB3, and Dd2.
Published ahead of print on 30 April 2008. ![]()
Supplemental material for this article may be found at http://jcm.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»