Standardization and application of a modified RFLP-PCR methodology for analysis of polymorphisms linked to treatment resistance in Ancylostoma braziliense
Parasites & Vectors volume 11, Article number: 540 (2018)
Single nucleotide polymorphisms (SNPs) in codons 167, 198 and 200 of the beta-tubulin isotype 1 gene are associated with benzimidazoles resistance in many helminths. Codon 167 mutation has never been described in hookworms; however, polymorphisms in codons 198 and 200 have been described for Ancylostoma caninum and Necator americanus. These mutations have never been investigated in Ancylostoma braziliense; therefore, it is not known if they are present in this species and whether they are correlated with treatment resistance. The RFLP-PCR technique has been used to analyze these polymorphisms in some nematodes, but depending on the species, these alterations do not create or eliminate any restriction enzyme cleavage site, making it impossible to use this technique. Here, we describe the standardization and application of a modified RFLP-PCR technique for detecting polymorphisms in individual A. braziliense worms recovered from naturally infected dogs in two Brazilian states.
The molecular techniques used were sensitive, specific, and easy to apply. To our knowledge, we report for the first time the presence of a polymorphism at codon 198 of the beta-tubulin gene of A. braziliense (1/81; 95% CI: 0–3.69%).
It is not known whether the presence of the mutation in codon 198 of the beta-tubulin gene of A. braziliense has importance for this parasite. However, based on studies of other helminths, it is possible that this polymorphism is directly related to the resistance to benzimidazoles. This may be a major concern, since this nematode has considerable relevance as a parasite of canids and felids and as one of the agents of cutaneous larva migrans in humans. Standardized methodologies will be useful for screening for polymorphisms in the beta-tubulin gene of canine hookworms in a broader population. The method could also be adapted for the analysis of other SNPs in other nematode species.
Ancylostoma braziliense is an intestinal nematode of dogs, cats and some wild animals, and can cause severe damage, such as anemia, to its hosts . This parasite has considerable zoonotic importance as an agent of erythematous and pruriginous inflammation in the dermis, known as cutaneous larva migrans .
The periodic treatment of helminths with benzimidazoles can select for treatment-resistant parasites . Single nucleotide polymorphisms (SNPs) at codons 167 (TTC/Phe→TAC/Tyr), 198 (GAG/Glu→GCG/Ala) and 200 (TTC/Phe→TAC/Tyr) in the beta-tubulin isotype 1 gene have been linked to benzimidazole resistance in some helminths [3,4,5]. Codon 167 mutation has never been described for hookworms . Conversely, SNPs have been described in codon 200 for Ancylostoma caninum  and in codons 198  and 200  for Necator americanus. These mutations have never been investigated in A. braziliense; therefore, it is not known if they are present in this species.
Some molecular methodologies have been described in the literature for genotyping SNPs in hookworms, such as ARMS-PCR [7, 10], Tetraprimer ARMS-PCR , real-time PCR (qPCR)  and Smart Amplification Process 2 (SmartAmp2) . Restriction fragment length polymorphism-polymerase chain reaction (RFLP-PCR) has been used successfully to detect SNPs linked to resistance, and it is a sensitive, specific, low cost and easy to apply technique. Moreover, the process of restriction enzyme digestion yields a quick result by gel analysis. According to Ota et al. , this technique is simple, sensitive, and reliable and requires minimal investment in instrumentation; however, it has some limitations, such as in cases where the DNA sequences are not recognized by commercial restriction enzymes or harbor many recognition sites for a single enzyme. Here, we standardize a simple and economical approach based on RFLP-PCR and site-directed mutagenesis followed by RFLP-PCR for genotyping SNPs in two Brazilian populations of A. braziliense. We did not have information on the dogs’ treatment history or on the possibility of drug resistance. To our knowledge, we describe for the first time the presence of a mutation at codon 198 of the beta-tubulin isotype 1 gene for this species.
Parasite collection and DNA extraction
Worms were recovered from dogs from two Brazilian state capitals, Teresina (Piauí) and Belo Horizonte (Minas Gerais). The dogs were routinely subjected to euthanasia following the approved procedures of the Municipality Health Centers of the respective Brazilian cities by the Zoonoses Control Center (CCZ). Adult worms were collected from 9 dogs from Piauí and 10 dogs from Minas Gerais during necropsy. Species identification was performed by examining the buccal capsule under a stereomicroscope . DNA was extracted employing a technique described by Waldschmidt et al.  from 81 individual worms, including 39 worms collected from Piauí (18 males and 21 females) and 42 worms collected from Minas Gerais (15 males and 27 females).
There was no A. braziliense beta-tubulin isotype 1 gene sequence available in the databases. Therefore, we tested primers designed for the beta-tubulin gene of A. caninum  and sequenced the A. braziliense regions flanking the codons of interest. The beta-tubulin isotype 1 gene sequence of A. braziliense was deposited in GenBank, under accession number MG866017. The primers used in this study were designed using the program Oligo Explorer 1.4™ (Gene link, USA). Additional file 1: Table S1 shows all the primers designed in this study.
Construction of controls
To standardize molecular analyses, we first synthesized controls for the presence and absence of mutations for all codons, according to Furtado & Rabelo . To construct a wild-type control allele for codon 167 (N167Ab), an initial PCR amplification was performed with the primers Fa167Ab and Ra167Ab (415 bp) using genomic DNA from A. braziliense. Since codons 198 and 200 are close to each other in the genome, it was possible to construct a single unmutated control for both codons using the primers Fa198/200Ab and Ra198/200Ab (424 bp). Next, nested PCR was performed using the primers Fc167Ab and Rc167Ab (306 bp, to codon 167) and the primers Fc198/200Ab and Rc198/200Ab (308 bp, to codons 198 and 200), and the fragments obtained were sequenced to confirm the absence of mutations in both amplified fragments. The fragments were subsequently cloned using the pGEM™-T Easy Vector System (Promega, Madison, USA), transformed into XL1-blue™ cells (Phoneutria, Belo Horizonte, Brazil) and recovered via miniprep (Wizard™ Plus Miniprep DNA Purification System, Promega, USA).
To construct mutated positive controls for codons 167 (M167Ab), 198 (M198Ab) and 200 (M200Ab), site-directed mutagenesis was performed using the Megaprimer-PCR technique. Cloned N167Ab was used as a template for PCR amplification of codon 167 using the primers Fc167Ab and Rm167Ab (176 bp), while cloned N198/200Ab was used as a template for PCR amplification of codon 198 and 200 using the primer combinations Fc198/200Ab and Rm198Ab (132 bp) and Fc198/200Ab and Rm200Ab (138 bp), respectively. The Rm167Ab, Rm198Ab and Rm200Ab primers were designed to include a mismatch at position 8 of the 5′ end of the primer (see Additional file 1: Table S1) to mimic the mutated sequence. The reaction products were subjected to electrophoresis on 1.0% agarose gels (w/v) (Midsci, Valley Park, USA) with 0.5× TAE buffer, and the gels were stained with GelRed™ (Biotium, Fremont, USA). The fragments were gel excised and purified (Illustra GFX™ PCR DNA and Gel Band Purification Kit, GE Healthcare, Chicago, USA), and their concentrations were determined. Approximately 20 ng of the first reaction product was used as a forward mega-primer in the second reaction in combination with the reverse Rc167Ab (306 bp, to codon 167) and Rc198/200Ab (308 bp to codon 198 and 200) primers. The products of the three codons were purified, sequenced, cloned and recovered in the same way as described for the negative controls. All PCR amplifications were performed using GoTaq Green Master Mix (Promega), with a final concentration of 0.2 μM for each primer. All amplification steps were performed in a Mastercycler (Eppendorf, Hamburg, Germany) thermocycler according to the following program: 94 °C for 5 min; 30 cycles of 94 °C for 1 min, 57 °C for 1 min, and 72 °C for 1 min; and 72 °C for 8 min. In all amplification runs, a “blank” sample was included in which the DNA was replaced with water.
Analysis of codons
For the analysis of SNPs in codons 167, 198 and 200 of the beta-tubulin gene of A. braziliense, the RFLP-PCR technique was performed. Sequences were analyzed using the NEBcutter V2.0 tool (http://www.labtools.us/nebcutter-v2-0/). The mutation at codon 167 created a restriction site for the enzyme RsaI. Therefore, this enzyme was used to evaluate this codon by RFLP-PCR. However, no restriction sites were found that could differentiate mutated from unmutated sequences for codons 198 and 200. Thus, to analyze codons 198 and 200, site-directed mutagenesis was required to create a site for the restriction enzymes DdeI and RsaI, respectively, to differentiate the genotypes by RFLP-PCR.
For codon 167, initial PCR using primers Fa167Ab and Ra167Ab (415 bp) was performed, followed by nested PCR with the primers Fc167Ab and Rc167Ab (306 bp). For codons 198 and 200, initial PCR with the primers Fa198/200Ab and Ra198/200Ab (415 bp) was performed. To introduce the restriction sites, primers for the nested reaction were designed with an alteration at position 3 from the 3′ end for codon 198, which in association with the remaining sequence created a site for the restriction enzyme DdeI in the unmutated allele. For codon 200, a mismatch was introduced in the primer Fsite200Ab at position 2 from the 3′ end that in combination with the remaining sequence, created an additional restriction site for the enzyme RsaI in the mutated allele. Additional file 2: Figure S1 and Additional file 3: Figure S2 show the methodology used for codons 198 and 200, respectively. Therefore, using the product of this first reaction, nested-PCR was performed for codon 198 with the primers Fsite198Ab and Rc198/200Ab (214 bp), and another nested PCR was performed for codon 200 with the primers Fsite200Ab and Rc198/200Ab (208 bp). The three PCR products were subjected to 1.0% (w/v) agarose gel electrophoresis (Midsci) using 0.5× TAE buffer and GelRed® staining (Biotium) to confirm the amplicon sizes. Restriction enzyme digestion of all of the codons was performed with 5 U (codons 167 and 200) or 2 U (codon 198) of each enzyme in individual reactions according to the manufacturer's recommendations. After digestion, the products were subjected to 15.0% (w/v) polyacrylamide gel electrophoresis in 0.5× TBE buffer, followed by GelRed® staining (Biotium). Additional file 4: Table S2 shows the expected fragment sizes after the digestion of each genotype. The confidence interval (CI) for the results was obtained using Graph Prism version 7.0 (GraphPad Software, USA).
Mutations in codons 167 and 200 in worms from dogs in both states were not found. For codon 198, only one of the specimens from Minas Gerais presented a mutation at codon 198, corresponding to 1.2% (1/81) of the total samples analyzed. This single mutation was present in a heterozygous worm. Figure 1 shows gels demonstrating the mutation screening at codons 167 (1a), 198 (1b) and 200 (1c).
SNPs in the beta-tubulin isotype 1 gene at codons 167 (TTC/Phe→TAC/Tyr), 198 (GAG/Glu→GCG/Ala) and 200 (TTC/Phe→TAC/Tyr) have been linked to benzimidazole resistance in some helminths [3,4,5]. Earlier detection of these SNPs increases the possibility of intervening to delay the establishment of a drug-resistant worm population. Therefore, there is a need for tests that are sensitive, specific and easy to apply. Here, we describe the standardization and application of a modified RFLP-PCR technique for detecting polymorphisms in individual A. braziliense worms from two Brazilian states.
None of the samples evaluated presented a mutation at codon 167. The polymorphism in this codon was first described for Haemonchus contortus by Prichard , and although this mutation was investigated in dog hookworms  and humans , it was never described for this group of parasites. In cyathostomes, the mutation at codon 167 has been suggested as the main mechanism of resistance to benzimidazoles . Due to the phylogenetic proximity , it was conceivable that this SNP would occur in hookworms; however, this hypothesis was not confirmed by the present study, but the SNP may be found in a larger study or in different geographical localities.
This report is the first description of mutation in codon 198 of the beta-tubulin gene of A. braziliense (1/81, 95% CI: 0–3.69). This mutation has been previously investigated in populations of A. caninum from Brazil  and the USA  but has never been reported for this species. According to Silvestre & Cabaret , it is possible that this single heterozygous mutation is sufficient to render the worm resistant to benzimidazole treatment, which is concerning. The frequency of this mutation found in the present work can be considered low compared to that found in studies of parasites of livestock . In general, research involving cattle parasites tends to show a high frequency of SNPs in the beta-tubulin isotype 1 gene . We suggest that this fact may be attributed to the control strategies used for different hosts; for example, cattle herds are periodically treated massively with anthelmintics for economic reasons, while this is not usually the case with dogs.
No mutation was observed at codon 200 in any sample analyzed. Rashwan et al.  and Schwenkenbecher et al.  also did not detect mutation in this codon in N. americanus, while Furtado et al.  and Diawara et al.  detected a mutation in this codon in A. caninum and N. americanus, respectively. There is no reason to believe that a polymorphism at codon 200 is not present in A. braziliense. The history of dogs used in this work is unknown, and there are no reports on the effectiveness of A. braziliense treatment in Brazilian populations of dogs or cats. It was not possible to know if these animals had already been treated with an anthelmintic. In fact, this work proposed to carry out a molecular scan for the natural presence of these SNPs in populations of worms isolated from dogs that had not necessarily received previous treatment to identify the potential risk of the large-scale use of this anthelmintic. We suggest that in vitro phenotypic tests should be performed to evaluate the sensitivity of hookworm populations in Brazil, which would help to more fully elucidate the resistance and sensitivity of this species to different anthelmintics.
A. braziliense, similar to A. caninum, is a zoonotic parasite, and although it is not able to develop in humans, it may be responsible for a cutaneous infection called cutaneous larva migrans . Considering that this class of drugs has the same mechanism of action against different developmental stages of the helminth, polymorphisms in the beta-tubulin gene may confer treatment resistance to larvae present in cutaneous infections. Accordingly, there are reports of failure with standard albendazole treatment for this type of infection, highlighting the possibility of treatment resistance . Ideally, to confirm the relationship between failures of benzimidazole treatment for larva migrans and the presence of these SNPs, it would be necessary to perform molecular analyses of resistant larvae directly collected from human skin lesions.
The molecular methodologies standardized here were sensitive and specific for the detection of mutations related to benzimidazole resistance. Since SNPs at codons 198 and 200 of A. braziliense do not create or eliminate a cleavage site for any restriction enzyme, it was necessary for the analysis of both codons to introduce a change in one of the PCR primer by site-directed mutagenesis, creating a cleavage site in association with the codon sequence of interest. In this methodology, larger primers for introducing the change in the DNA will better distinguish the genotypes in the gel. However, in the analyses for A. braziliense, the intentionally altered primers (Fsite198Ab and Fsite200Ab) were 30 bp, only slightly larger than the usual primers, although they were able to differentiate the genotypes of interest.
Our group performed analyses of these SNPs in A. caninum populations collected from different Brazilian states. These analyses were performed for codons 167 and 200 by ARMS-PCR [7, 10], whereas codon 198 was analyzed by Tetraprimer ARMS-PCR . Although standardized techniques produced reliable results, with high sensitivity and specificity, the standardization was very laborious. The main drawback of these methodologies is that the primers need to strictly pair with the target DNA sequence, without allowing for adjustments in primer pair compatibility. Once the primers pair is not compatible, it becomes impossible to design new primers, rendering these techniques infeasible . However, this type of problem is not observed for the use of RFLP-PCR since, depending on the objective, amplification can be performed using primers that are designed to any sequence that flanks the target region. Furthermore, this method can be easily applied in DNA genotyping of multiple organisms.
Studies based on qPCR and sequencing are described in the literature for the analysis of SNPs in the beta-tubulin isotype 1 gene of helminths [12, 24]. However, these methods, although sensitive, require specific equipment, while RFLP-PCR techniques only require the use of a conventional thermocycler, providing a direct result without the need for graph analysis or chromatograms. Schwenkenbecher et al.  detected a mutation in codons 167 and 200 at a very low frequency in N. americanus from children who received treatment periodically and associated these data with a possible experimental error in qPCR. Rashwan et al.  developed a genotyping assay for SNPs at codons 167, 198 and 200 in N. americanus beta-tubulin isotype 1 gene using the SmartAmp2 method. According to these authors, the advantages of this technique are that sophisticated equipment is not needed and only a short time is required for the reaction. On the other hand, the combination of several initiators (five to six) in the same reaction may hinder the proper functioning of these reactions.
We conclude that a mutation in codon 198 in the beta-tubulin gene of A. braziliense exists in the population evaluated, although at a low frequency. The presence of the mutation in this population associated with recurrent treatment can select for resistant strains of parasites, suggesting the need for constant monitoring. This may represent a major problem since this nematode has considerable relevance as a parasite of canids and felids and as an agent of cutaneous larva migrans in humans.
Amplification-refractory mutation system
Polymerase chain reaction
Restriction fragment length polymorphism
Rep BH. Pathogenicity of Ancylostoma braziliense. IV. Blood loss caused by the worms in the prepatent period. Trop Geogr Med. 1966;18:329–52.
Suber WJ. Creeping eruption due to Ancylostoma braziliense. A case report. J Am Podiatry Assoc. 1962;52:128–9.
Ghisi M, Kaminsky R, Mäser P. Phenotyping and genotyping of Haemonchus contortus isolates reveals a new putative candidate mutation for benzimidazole resistance in nematodes. Vet Parasitol. 2007;144:313–20.
Kwa MS, Veenstra JG, Roos MH. Benzimidazole resistance in Haemonchus contortus is correlated with a conserved mutation at amino acid 200 in beta-tubulin isotype 1. Mol Biochem Parasitol. 1994;63:299–303.
Prichard RK. Genetic variability following selection of Haemonchus contortus with anthelmintics. Trends Parasitol. 2001;17:445–53.
Furtado LF, Bello AC, Rabelo ÉM. Benzimidazole resistance in helminths: From problem to diagnosis. Acta Trop. 2016;162:95–102.
Furtado LF, Bello AC, Santos HA, Carvalho MR, Rabelo ÉM. First identification of the F200Y SNP in the β-tubulin gene linked to benzimidazole resistance in Ancylostoma caninum. Vet Parasitol. 2014;206:313–6.
Rashwan N, Bourguinat C, Keller K, Gunawardena NK, de Silva N, Prichard R. Isothermal diagnostic assays for monitoring single nucleotide polymorphisms in Necator americanus associated with benzimidazole drug resistance. PLoS Negl Trop Dis. 2016;10:e0005113.
Diawara A, Halpenny CM, Churcher TS, Mwandawiro C, Kihara J, Kaplan RM, et al. Association between response to albendazole treatment and β-tubulin genotype frequencies in soil-transmitted helminths. PLoS Negl Trop Dis. 2013;7:e2247.
Furtado LF, Rabelo ÉM. Molecular analysis of the F167Y SNP in the β-tubulin gene by screening genotypes of two Ancylostoma caninum populations. Vet Parasitol. 2015;210:114–7.
Furtado LF, Alves WP, Moreira TB, Costa Junior LM, Miranda RR, Rabelo ÉM. Standardization and application of the tetraprimer ARMS-PCR technique for screening of the E198A SNP in the β-tubulin gene of hookworm populations in Brazil. Vet Parasitol. 2016;210:114–7.
Schwenkenbecher JM, Kaplan RM. Real-time PCR assays for monitoring benzimidazole resistance-associated mutations in Ancylostoma caninum. Exp Parasitol. 2009;122:6–10.
Ota M, Fukushima H, Kulski JK, Inoko H. Single nucleotide polymorphism detection by polymerase chain reaction-restriction fragment length polymorphism. Nat Protoc. 2007;2:2857–64.
Silva LM, Miranda RR, Santos HA, Rabelo EM. Differential diagnosis of dog hookworms based on PCR-RFLP from the ITS region of their rDNA. Vet Parasitol. 2006;140:373–7.
Waldschmidt AM, Salomão TM, Barros EG, Campos LA. Extraction of genomic DNA from Melipona quadrifasciata (Hymenoptera: Apidae, Meliponinae). Braz J Genet. 1997;20:421–3.
Furtado LF, Rabelo ÉM. Development of a new amplification-refractory mutation system for detection of a single nucleotide polymorphism linked to drug resistance in Ancylostoma caninum. Genet Mol Res. 2015;14:5103–11.
Ishii JB, Arenal A, Felix A, Yoshitani U, Beech R, Molento MB. Diagnosis of resistance alleles in codon 167 of the beta-tubulin (Cya-tbb-1) gene from third-stage larvae of horse cyathostomins. Res Vet Sci. 2017;115:92–5.
Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, Vierstraete A, et al. A molecular evolutionary framework for the phylum Nematoda. Nature. 1998;392(6671):71–5.
Silvestre A, Cabaret J. Mutation in position 167 of isotype 1 beta-tubulin gene of trichostrongylid nematodes: role in benzimidazole resistance? Mol Biochem Parasitol. 2002;120:297–300.
Dixit AK, Das G, Dixit P, Singh AP, Kumbhakar NK, Sankar M, et al. An assessment of benzimidazole resistance against caprine nematodes in Central India. Trop Anim Health Prod. 2017;49:1471–8.
Chaudhry U, Miller M, Yazwinski T, Kaplan R, Gilleard J. The presence of benzimidazole resistance mutations in Haemonchus placei from US cattle. Vet Parasitol. 2014;204:411–5.
Schwenkenbecher JM, Albonico M, Bickle Q, Kaplan RM. Characterization of beta-tubulin genes in hookworms and investigation of resistance-associated mutations using real-time PCR. Mol Biochem Parasitol. 2007;156:167–74.
Quashie NB, Tsegah E. An unusual recurrence of pruritic creeping eruption after treatment of cutaneous larva migrans in an adult Ghanaian male: a case report with a brief review of literature. Pan Afr Med J. 2015;21:285.
Lambert SM, Nishi SM, Mendonça LR, Souza BM, Julião FS, Gusmão PS, et al. Genotypic profile of benzimidazole resistance associated with SNP F167Y and F200Y beta-tubulin gene in Brazilian populations of Haemonchus contortus of goats. Vet Parasitol. 2017;8:28–34.
This work received financial support from Fundação de Amparo à Pesquisa de Minas Gerais - FAPEMIG (Process number APQ-02417-16) and Pró-Reitoria de Pesquisa from Universidade Federal de Minas Gerais (PRPq/UFMG).
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. The datasets for the present study are available from the corresponding author upon request.
Ethics approval and consent to participate
All animal procedures were approved by the Animal Care Ethics Committee at the Universidade Federal de Minas Gerais (Protocol number 328/12) and were performed under the guidelines of the Brazilian Council of Animal Experimentation (CONCEA), strictly following the Brazilian law for “Procedures for the Scientific Use of Animals” (11.794/2008).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primers used for SNP analysis in this work with their respective annealing temperatures and base substitutions (when applicable). The bases that have been replaced are in bold. (PDF 76 kb)
Figure S1. Representation of the methodology used to analyze codon 198 of the beta-tubulin isotype 1 gene from Ancylostoma braziliense. The Fsite198Ab primer was designed to add a mutation in the amplicon, regardless of whether there was a mutation at codon 198. The absence of the mutation at codon 198 in combination with the primer-introduced change creates a site for DdeI. If the allele is mutated, even the altered primer cannot create a DdeI site. Codon 198 is underlined, with the base of interest in bold. (PDF 615 kb)
Figure S2. Representation of the methodology used to analyze codon 200 of the beta-tubulin isotype 1 gene from Ancylostoma braziliense. The Fsite200Ab primer was designed to add a mutation in the amplicon, regardless of whether there was mutation at codon 200. The presence of the mutation at codon 200 in combination with the primer-introduced change creates a site for RsaI. If the allele does not harbor a mutation, even the altered primer cannot create a RsaI site. Codon 200 is underlined, with the base of interest in bold. (PDF 88 kb)
Table S2. Amplicon sizes after RFLP-PCR for the analysis of codons 167, 198 and 200 of the beta-tubulin isotype 1 gene of Ancylostoma braziliense. (PDF 30 kb)
About this article
Cite this article
Furtado, L.F.V., Magalhães, J.G.S. & Rabelo, É.M.L. Standardization and application of a modified RFLP-PCR methodology for analysis of polymorphisms linked to treatment resistance in Ancylostoma braziliense. Parasites Vectors 11, 540 (2018). https://doi.org/10.1186/s13071-018-3125-9