|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
The mode of action of benzimidazoles was found to be their interference with tubulin polymerization into microtubules by binding to ß-tubulin.57 Resistance against benzimidazoles in a number of nematode species has been shown to be caused by a phenylalanine (TTT/TTC) to tyrosine (TAT/TAC) substitution at either position 167 or 200 of nematode ß-tubulin.816 The position 200 tyrosine mutation appears to be more common in parasitic nematodes and was found to be recessive.11,16 Resistance to benzimidazoles was also found to be associated with 200Tyr or 167Tyr substitutions in normally benzimidazole-susceptible fungi,1720 and the 200Tyr mutation in ß-tubulin was also found in benzimidazole-resistant strains of protozoa.21 However, the position 200 mutation was not found in Necator americanus hookworms from Pembla Island in which benzimidazole resistance was suspected.22 Eukaryotes that normally code for tyrosine at codon 200 appear to be uniformly tolerant for benzimidazole. Phenylalanine at codons 167 and 200 appears to be critical for the high-affinity binding of benzimidazole anthelmintics to nematode tubulin and for sensitivity to these anthelmintics.12
In a study carried out in India,23 pooled W. bancrofti microfilariae from 14 patients from five geographic regions were sequenced to detect the codon 200 resistance mutation. The consensus sequences of the 14 pooled samples failed to show the codon 200 resistance mutation. However, sequencing the pooled samples would not have shown the mutation unless it was abundant, since sequencing pooled DNA samples is not a sensitive tool for the detection of rare single nucleotide polymorphisms (SNPs). The sequencing result of a pooled sample will show the most abundant sequence.
The objectives of this study were to assess whether mutations associated with benzimidazole resistance in other nematodes are present in W. bancrofti, to examine their frequency in untreated populations from Ghana and Burkina Faso, and to determine their frequency in microfilaria obtained from patients treated with albendazole and ivermectin.
| MATERIALS AND METHODS |
|
|
|---|
Amplification of Brugia malayi and W. bancrofti ß-tubulin. Primers were designed based on the B. pahangi ß-tubulin sequence (GenBank accession no. M36380) to amplify the full length ß-tubulin cDNA from B. malayi. Total RNA was extracted from bulk adult worms (obtained from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD) using TRIzol reagent (Invitrogen, Carlsbad, CA) and methods described by the manufacturer. Total RNA was reverse transcribed with the adaptor primer (5'-3') GGC CAC GCG TCG ACT AGT AC(T)17 and 200 units of murine moloney leukemia virus reverse transciptase (Invitrogen) according to the manufacturers instructions. cDNA was used to amplify a fragment of ß-tubulin using a reaction containing 2.5 µL of 10 x polymerase chain reaction (PCR) buffer (200 mM Tris-HCl, pH 8.4, 500 mM KCl), 1 µL (50 µM) of MgCl2, 0.5 µL each of primers sense (5'-3') GGT ACC ATG GAT TCT ATT CG (AY705382 position 211-230) and antisense (5'-3') GAT CGG CGT TCA ACT GTC CA (position 729-748) at a concentration of 12.5 µM, 1 unit of Taq polymerase, 1 µL (10 µM) of deoxynucleotide triphosphates (dNTPs), and sterile water to give a final volume of 25 µL. The amplification conditions were 94°C for 2 minutes, followed by 30 cycles at 94°C for 55 seconds, 55°C for 55 seconds, and 72°C for 55 seconds, and a final extension at 72°C for 10 minutes.
The full-length sequence was obtained using a 3' rapid amplification of cDNA ends (RACE) procedure. Briefly, 5 µL of cDNA was amplified in a nested reaction containing 2.5 µL of 10 x PCR buffer (200 mM Tris-HCl, pH 8.4, 500 mM KCl), 1.5 µL (50 µM), of MgCl2, 0.5 µL each of the sense primers (5'-3') GGC AAA TAT GTG CCA CGA GC (position 169-188) and (5'-3') GGT ACC ATG GAT TCT ATT CG (position 211-230) at a concentration of 20 µM, 1 unit of Taq polymerase, 1 µL (10 µM) of dNTPs, and sterile water to give a final volume of 25 µL. The amplification conditions for both steps were 94°C for 2 minutes, followed by 40 cycles at 94°C for 30 seconds, 56°C for 30 seconds, and 72°C for 75 seconds, and a final extension at 72°C for 10 minutes. The 5' end of the gene was amplified using the same conditions as the 3' RACE, but with different primers: sense primer, which corresponded to the nematode splice leader sequence SL1 (5'-3') GGT TTA ATT ACC AAG TTT GAG, and antisense primers (5'-3') AGG GCT CGA TAA GCA GCA GC (position 832-851) and (5'-3') GAT CGG CGT TCA ACT GTC CA (position 729-748). The resulting fragments were cloned into a PCR 2.1 vector using a TA cloning kit (Invitrogen) as per the manufacturers instructions and sequenced with a Beckman (Palo Alto, CA) CEQ DNA sequencer.
Primers to amplify genomic DNA from W. bancrofti were designed based on the full-length B. malayi ß- tubulin cDNA to amplify by PCR a region surrounding the mutations to be diagnosed from W. bancrofti microfilaria in a reaction containing 2.5 µL of 10 x PCR buffer, 1 µL (50 µM) of MgCl2, 0.5 µL each of the primers sense (5'-3') GGC AAA TAT GTG CCA CGA GC (AY705382 position 169-188) and (5'-3') AGG GCT CGA TAA GCA GCA GC (AY705382 position 832-851) at a concentration of 12.5 µM, 1 unit of Taq polymerase, 1 µL (10 µM) of dNTPs, and sterile water to give a final volume of 25 µL. The amplification conditions were 94°C for 2 minutes, followed by 30 cycles at 94°C for 55 seconds, 55°C for 55 seconds, and 72°C for 55 seconds, and a final extension at 72°C for 10 minutes. All reactions were carried out in a PTC-200 Thermal Cycler (MJ Research, Inc., Waltham, MA). Resulting fragments were cloned and sequenced as described earlier in this report.
Diagnosis of resistance-associated mutations in W. bancrofti. Blood samples (10 mL) containing microfilaria were transferred to centrifuge tubes and centrifuged at 1,100 rpm (300 x g) for 10 minutes at 4°C. The pellet was resuspended in 25 mL of red blood cell lysis solution, incubated on ice for 10 minutes, and centrifuged again as above. The final pellet was resuspended in 1.5 mL of 70% isopropyl alcohol. To separate individual worms, a small amount of the sample was transferred to a small Petri dish. Single microfilariae were washed in 70% ethanol, separated into individual PCR-microcentrifuge tubes using a micropipette under an inverted microscope, and identified visually as W. bancrofti. Extraction of DNA was carried out directly in the PCR tube using the Qiagen (Valencia, CA) DNeasy tissue kit at one-fourth the recommended volume and following the manufacturers instructions.
To detect mutations using a pyrosequencer, a smaller fragment surrounding the SNPs was amplified using microfilarial DNA or, when the DNA content was small, the product of the PCR amplification performed earlier was used as a template. The reaction contained 5 µL of 10 x PCR buffer, 2 µL (50 µM) of MgCl2, 0.5 µL each of primers sense (5'-3') (AY705383 position 308-328) GGG AAC ATT GCT GAT CTC GAA and antisense (5'-3') GGA AGC AGA TGT CAT ACA AAG CC (AY705383 position 672-694) at a concentration of 20 µM, 1 unit of Taq polymerase, 2 µL (10 µM) of dNTPs, and sterile water to give a final volume of 50 µL. The second primer had biotin attached to its 3' end to generate a single-stranded template using streptavidin beads. The amplification conditions were 94°C for 2 minutes, followed by 40 cycles at 94°C for 15 seconds, 55°C for 30 seconds, and 72°C for 30 seconds, and a final extension at 72°C for 10 minutes. The sequencing primers used for SNP analysis in a Pyrosequencer® (Pyrosequencing, Inc., Westborough, MA) were (5'-3') CGG ATC GAA TTA TGA GCT CT (position 346-365) for position 167 and (5'-3') GAA AAC ACT GAC GAA ACT T (position 637-655) for position 200.
Real-time diagnosis of mutations at position 200 and 167 was done in a single reaction in a Roche (Mississauga, Ontario, Canada) Light Cycler®. The reaction contained 1.6 µL (25 µM) of MgCl2, 0.4 µL each of the primers sense (5'-3') GGT ACC ATG GAT TCT ATT CG (AY705383 position 1-20) and antisense (5'-3') GAT CGG CGT TCA ACT GTC CA (position 907-926) at a concentration of 12.5 µM, 0.2 µL each of the probes 5'-CCG TGA GGA GTA TCC GGA TCG AAT TAT G (fluorescein)-3' (position 332-359) and 5'-(LC Red 640)CTC TTT TTC GGT TGT GCC GTC G (phosphate)-3' (position 362-382) for position 167 and 5'-ACT TTC TGC ATT GAT AAC GAG GC (fluorescine)-3' (position 652-674) and 5'-(LC Red 705) TGT ATG ACA TCT GCT TCC GAA CGT (phosphate)-3' (position 677-700) for position 200 at a concentration of 20 µM. The reaction conditions were 20 sec at 95°C, followed by 40 cycles at 95°C for 10 seconds, 55°C for 10 seconds, and 72°C for 40 seconds, a melting step from 50°C to 80°C, and a final cooling step to 40°C. All temperature changes were with a ramp speed of 20°C/second. Melting temperatures for the position 200 probes are 67°C for the wild-type gene and 63°C for the mutant. Melting temperatures for the position 167 probes are 63 °C for the wild-type gene and 58.5°C for the mutant.
Both diagnostic tests were optimized with a plasmid containing the wild-type gene and plasmids containing the mutation(s), generated by site-directed mutagenesis, using the QuickChange kit (Stratagene, La Jolla, CA) and the following mutagenesis primers: sense (5'-3') ATG AGC TCT TAT TCG GTT GTG CCG TCG (AY705383 position 357-383) and antisense (5'-3') ACA ACC GAA TAA GAG CTC ATA ATT CG (position 351-377) for the mutation at codon 167, and sense (5'-3') CGA AAC TTA CTG CAT TGA TAA CG (position 648-670) and antisense (5'-3) CGT TAT CAA TGC AGT AAG TTT CG (position 648-670) for the mutation at codon 200, following the manufacturers instructions.
Experimental design and statistical analyses. Microfilariae from patients infected with W. bancrofti were collected in Burkina Faso and Ghana. Treatment programs for lymphatic filariasis have been instituted in communities in Ghana and Burkina Faso and records and patient questionnaires were available to establish the treatment histories in the individuals sampled. We analyzed microfilaria from patients before treatment, after one round of treatment, and after two rounds of treatment with 400 mg of albendazole and 200 µg/kg of ivermectin. In the samples obtained from Burkina Faso, 515 microfilariae per patient were genotyped. Microfilariae from Ghana were obtained as pooled samples from several patients. Two pooled samples were obtained that had been collected at different times in 2001. These two samples were treated separately because they were sampled from different populations. A chi-square test was used to detect significant differences in genotype frequencies and in allele frequencies between the three treatment groups within each country. It was determined whether genotype frequencies of the whole population and within each patient were in Hardy-Weinberg equilibrium using a chi-square analysis to compare observed and expected values.
| RESULTS |
|
|
|---|
We obtained several hundred microfilaria from Ghana and Burkina Faso. Ghanaian samples consisted of microfilaria from patients with lymphatic filariasis who had not been treated with albendazole and ivermectin as part of the Global Alliance to Eliminate Lymphatic Filariasis. Blood samples containing W. bancrofti collected in Burkina Faso were from either untreated patients or patients treated once or twice with albendazole and ivermectin. Samples were taken at least seven days after the last anthelmintic treatment. Microfilaremia was low in the treated patients, although the patients were still microfilaria positive despite recent treatment. The first two batches of microfilaria from Ghana were analyzed using the FRET assay. We were able to identify several microfilariae that had the position 200 tyrosine genotype. The frequency of the mutant allele in the first group of 130 microfilariae was 2.7%. The frequency of the mutant allele in the second group of 300 microfilariae was 0.33% (Table 1
). Because of a lack of information regarding the size of the human population from which the microfilariae were pooled, we cannot make a reliable estimate of the frequency of these resistance-associated mutations in the entire parasite population. However, these findings indicate that the allele, which potentially causes albendazole resistance, is present in the population prior to treatment with albendazole. This suggests that with the advent of the mass treatment program in Ghana, selection for this allele in the population could occur.
|
|
2 = 139.8). The genotype frequencies for all three groups from Burkina Faso are shown in Figure 1
|
| DISCUSSION |
|
|
|---|
We have shown that there was a significantly higher allele frequency of the 200 tyrosine genotype in worms collected from patients who had been treated with albendazole in combination with ivermectin. This indicates that this mutation may be selected for with drug treatment. The total allele frequency was 31.6% higher in the treated versus the non-treated worm population. In worms collected from patients who had been exposed to two rounds of albendazole/ ivermectin treatment, the allele frequency was an additional 25.9% higher. Such selection is a strong indication of the development of drug resistance. Because pretreatment microfilarial counts were not taken, we cannot make conclusions about the success of treatment. However, it was only those microfilaria that were still present in the treated patients after treatment that were genotyped. In the case of the patients treated twice, microfilaria were harvested seven days after their last treatment and could therefore represent a population of microfilarial survivors of treatment or progeny of adult worms that were able to maintain microfilarial production despite treatment.
In some parasitic nematodes, a second ß-tubulin gene appears to be present,26,27 which may contribute to benzimidazole sensitivity and resistance.24 We did not detect a second ß-tubulin isotype in W. bancrofti. However, if present, it could contribute to maintaining susceptibility to albendazole.
It must also be noted that this was not a prospective study because parasites from different treatment groups were collected from different patients residing in separate villages. Therefore, we cannot rule out natural genetic variation, which may occur between populations, as being responsible for the differences in polymorphism between the treated and untreated populations.
We tested whether alleles of this gene were in Hardy-Weinberg equilibrium. When considering all worms in a treatment group (including untreated patients) as part of one population, the frequencies were not in balance, but showed an excess of homozygotes. When considering all worms sampled from an individual patient as being one population, the sum of the resulting allele frequencies from all patients was in Hardy-Weinberg equilibrium. The resistant-associated allele appeared to be aggregated within certain patients, rather than evenly distributed throughout the population. This result may be due to the fact that random mating does not occur between worms in different hosts. An excess of homozygotes may strongly accelerate selection for resistance, since it is believed to be a completely recessive trait,10 thus a larger proportion of the 200tyr as homozygotes could lead to a greater proportion of treatment failures.
Combination treatment with ivermectin and albendazole could be expected to reduce the rate of selection for resistance to either of these anthelmintics provided that the resistance mechanisms involve different genes. However, it has recently been found that ivermectin selects on ß-tubulin in the filarial nematode Onchocerca volvulus (Eng KLJ and Prichard RK, unpublished data). Although ivermectin selection is not associated with the phe200tyr SNP examined here, selection by the two anthelmintics can involve the same gene. Thus, there could be linkage between the region implicated in ivermectin selection and possible albendazole selection for the phe200tyr SNP, which could modify the rate of selection for albendazole resistance. The implications of combination treatment on genetic selection need to be examined in future studies. Based on these findings, it is imperative to continue monitoring for the presence of this mutation, in other treatment zones, to detect early evidence of possible resistance selection and to correlate this with responses to treatment.
Received September 2, 2004. Accepted for publication January 11, 2005.
Acknowledgments: We thank Drs. Brian Bagnall and Mark Bradley (GlaxoSmithKline) for their interest and encouragement.
Financial support: This study was supported by the Formation de Chercheurs et lAide à la Recherche (Quebec, Canada), the Natural Sciences and Engineering Research Council of Canada (Ottawa, Ontario, Canada), the Centre for Host Parasite Interactions (Sainte-Anne-de-Bellevue, Quebec, Canada), and GlaxoSmithKline (Brentford, Middlesex, United Kingdom).
* Address correspondence to Roger K. Prichard, Institute of Parasitology, McGill University, 21111 Lakeshore Road, Sainte-Anne-Bellevue, Quebec H9X 3V9, Canada. E-mail: roger.prichard{at}mcgill.ca ![]()
Authors addresses: Anne E. Schwab, Institute of Parasitology, McGill University, 21111 Lakeshore Road, Sainte-Anne-de-Bellevue, Quebec H9X 3V9, Canada, Telephone: 514-398-8767, Fax: 514-398-7857, E-mail: aschwa4{at}po-box.mcgill.ca. Daniel A. Boakye, Noguchi Memorial Institute for Medical Research, PO Box LG581, Legon, Ghana, Telephone: 233-21-501-178, Fax: 233-21-502-182, E-mail: dboakye{at}noguchi.mimcom. Dominique Kyelem, Ministry of Health, 01 BP 2935, Ouagadougou 01, Burkina Faso, Telephone: 226-30-87-90, Fax: 226-326-335, E-mail: dominiquekyelem{at}yahoo.fr. Roger K. Prichard, Institute of Parasitology, McGill University, 21111 Lakeshore Road, Sainte-Anne-de-Bellevue, Quebec H9X 3V9, Canada, Telephone: 514-398-7961, Fax: 514-398-7857, E-mail: roger.prichard{at}mcgill.ca.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
E. Ghedin, S. Wang, D. Spiro, E. Caler, Q. Zhao, J. Crabtree, J. E. Allen, A. L. Delcher, D. B. Guiliano, D. Miranda-Saavedra, et al. Draft Genome of the Filarial Nematode Parasite Brugia malayi Science, September 21, 2007; 317(5845): 1756 - 1760. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. MCCARTHY IS ANTHELMINTIC RESISTANCE A THREAT TO THE PROGRAM TO ELIMINATE LYMPHATIC FILARIASIS? Am J Trop Med Hyg, August 1, 2005; 73(2): 232 - 233. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |