Skip to main content
  • Research note
  • Open access
  • Published:

Can P-glycoprotein and β-tubulin polymorphisms be used as genetic markers of resistance in Dirofilaria immitis from Rio de Janeiro, Brazil?

Abstract

Objective

Dirofilaria immitis, the causative agent of canine heartworm infection, is worldwide the most important filarid to affect domestic dogs. Prevention of this infection is done by macrocyclic lactones, but some reports on the lack of efficacy have been published. Although the actual cause of resistance is unknown, single nucleotide polymorphisms (SNPs) on a P-glycoprotein ABC transporter and β-tubulin genes have been pointed out as candidates for genetic markers of resistance. We conducted a survey to verify the presence of these suggested genetic markers in microfilariae from 30 naturally infected dogs under macrocyclic lactones treatment living in an endemic area in the state of Rio de Janeiro.

Results

The analysis of these specific SNPs demonstrated no sign of polymorphism on the P-glycoprotein loci, while 72 and 48% of the samples were polymorphic to the first and second SNPs on β-tubulin loci, respectively. This work demonstrates that the P-glycoprotein position 11 and 618 were not polymorphic and, therefore, not suitable as a genetic marker of resistance in Rio de Janeiro whereas both β-tubulin loci were polimorphic. This work points out the difficulty of finding a universal genetic marker for resistance.

Introduction

The filarial nematode, Dirofilaria immitis, the causative agent of canine heartworm infection, is worldwide the most important filarid to affect domestic dogs. Since the 80’s, prevention of this infection is done by macrocyclic lactones (ML), but some reports on the lack of efficacy have been published in the last 12 years. Resistance to macrocyclic lactones (ML) is a major problem in small ruminants, horses and cattle, leading to a discussion on the way these drugs should be used in livestock farms [1]. In dogs, the first report on ML loss of efficacy in D. immitis prevention has been published in 2005 [2]. Resistance in canine worms should be more difficult to develop than in ruminant helminths [3] because stray and feral dogs, as well as wild canids, carry a population of parasites that is not exposed to drug selection and therefore represent a wide refugia. The refugia bring drug-susceptibility alleles that dilute out the resistant ones in the population, reducing the selection pressure for drug resistance [4]. Nevertheless, a study using microfilariae obtained from dogs from US and Japan has shown the genetically variability of D. immitis, suggesting that the development of resistance could be possible under a high selective pressure [5]. Studies conducted using microfilariae obtained from three dogs originally from Louisiana and Arkansas, USA, that were under unsuccessful ML heartworm preventive treatment demonstrated an excess of homozygosity, which might be a sign of selection [6]. In addition, these microfilariae presented high frequency of a GG–GG genotype on two single nucleotide polymorphisms (SNPs) at positions 11 and 618 of a P-glycoprotein (Pgp) ABC transporter gene (i.e. an homozygous G was present on both positions). These findings suggest that the Pgp genotype might be useful as a genetic marker to assay for low responding D. immitis in the field [7]. A total of 9 ABC transporter genes have been identified in D. immitis and both SNPs are located in the Dim-pgp-11 [8]. Interestingly, in Parascaris equorum, resistance to ML was correlated with the presence of 3 individual SNPs and with an increase in the level of Peq-pgp-11 mRNA [9].

β-tubulin (Tub) has also been pointed out as a candidate gene for carrying genetic markers for resistance since it presents two SNPs that are in linkage disequilibrium in resistant microfilariae [5].

In the present study microfilaremic dogs living in the state of Rio de Janeiro, Brazil, received ML monthly for 4 months and a month after the forth dose microfilariae were collected. We aimed to identify the presence of the genotype for the two SNPs in the Pgp and the two others in Tub genes in worms from the studied area.

Main text

The study was conducted in two coastal areas of the state of Rio de Janeiro in sandbanks of the municipalities of Arraial do Cabo and Araruama, where the dogs are exposed to a high challenge, 80.8 and 58.4% respectively [10]. Both areas are located at sea level, apart from one another for 10 km. Each area is located in opposite sides of a salt water lagoon. Yearound climate conditions in the area favor mosquitoes population and therefore heartworm transmission may occur at any time of the year. Microfilaremic and antigenemic dogs were included in the study after positive SNAP 4DX and modified Knott’s [11] tests, once their microfilaraemia was over 300 microfilariae/mL of blood [12]. Dogs’ history of ML treatment was not taken into consideration as there was no trustful record of it. This study was approved by Fundação Oswaldo Cruz (FIOCRUZ) Institutional Animal Care and Use Committee (CEUA—LW-33/11). All dogs received four consecutive monthly doses of avermectin [ivermectin–0.06–0.12 mg/kg (Endogard, Virbac) or selamectin–6 mg/kg (Revolution, Zoetis)] following manufacturer’s guidelines. All doses were administered by the study team. Prior to treatment and a month after the last dose, microfilaremia was blind counted using 2 fixed giemsa-stained blood slides of each dog, each one counted by the same two researchers [12]. Each dog’s microfilaraemia was calculated by the mean of the 4 counts. Statistical tests were performed in GraphPad Prism software, using Wilcoxon test. A month after the last dose, microfilariae were isolated from 1 mL of fresh blood samples by filtration through polycarbonate membranes [13]. Microfilariae were fixed in isopropanol and kept under refrigeration. DNA from pools of 50–100 microfilariae from each dog was extracted using phenol/chloroform method [14] and submitted to PCR amplification with Platinum Taq DNA Polymerase (Invitrogen). The Pgp gene fragment (GenBank HM596853) and the first Tub fragment (position 28 to 719) were amplified using primers described elsewhere [5], while the second Tub fragment (position 2364–2961) was amplified using the primers 5′ TGCTGAGCTCACTCAACAGG 3′ and 5′ GAATGTTGCGCTCATCTTCA 3′. After confirmation by electrophoresis on a 1.5% agarose gel, the amplified fragments were sequenced using a 3730XL DNA analyzer system at the Sequencing Plataform at FIOCRUZ (Plataforma de Sequenciamento de DNA, PDTIS/FIOCRUZ). After sequencing, each individual chromatogram was analyzed in SeqScape 2.1 (Applied Biosystems) and sequences were aligned in BioEdit Sequence Alignment Editor 7.1.11 [15]. Sequences with multiple peaks along the whole fragment at least in two sequencings of the same pool were discarded.

Counting of microfilariae before the first dose and 30 days after the last ML treatment dose showed that, in spite of the four doses, all 30 dogs remained microfilaremic. Microfilaraemia in Araruama canine samples declined (12/12; 100%, P = 0.0005), while in Arraial do Cabo canine samples microfilariae counting remained similar (7/18 declined; 39%, P = 0.2121). Ten dogs from Arraial do Cabo had an increase in microfilaremia, in spite of ML treatment (Fig. 1). Samples of microfilariae from all these dogs, 30 days after the end of treatment, were sequenced. Positions 11 and 618 of the Pgp fragment could be analyzed in 28 and 26 pool samples respectively, while position 561 and 2755 of the Tub gene could be analyzed in 25 and 27 pool samples respectively.

Fig. 1
figure 1

Microfilaria counting—Number of microfilariae/ml calculated from counting blood smears prior to treatment (0 days) and 30 days after treatment end (120 days). Treatment was performed with four consecutive monthly doses of avermectins (ivermectin-0.06–0.12 mg/kg or selamectin—6 mg/kg) following manufacturer’s guidelines (***P < 0.001, Wilcoxon test)

The analysis of Tub chromatograms (Additional file 1) showed the presence of heterozygosity on 72% of samples at position 561 (18/25) and 48% of samples at position 2755 (13/27), showing that these positions are polymorphic in D. immitis circulating among dogs in Rio de Janeiro. On the other hand, none of the Pgp chromatograms showed any sign of polymorphism as a single ‘A’ peak on position 11 and a single ‘G’ peak on position 618 were detected (Table 1). In spite of the SNPs found in the Tub gene, no relation was observed between increase in microfilaremia and presence of SNP (Fisher’s exact test, P = 0.1581) and SNPs could be found in dogs from both areas.

Table 1 Quantitative SNP results

According to our results, it seems that the described Pgp SNPs do not occur in worms of the studied area and therefore are not suited as genetic markers of worm resistance to ML in this area. The Tub SNPs, although present, did not seem to be related to resistance. This information raises the issue that their use might be restricted to a specific geographic region and could not be accepted globally. However, since the limit of detection for alternative alleles with Sanger sequencing is 15–20% [3] it must be considered that an undetectable low frequency of SNPs could be present in these microfilariae as we used pools.

Nevertheless the comparison of the adult worms from US and Japan dogs showed important differences among the polymorphic loci [5]. This fact alone highlights the difficulty in finding a global genetic marker. Likewise in other helminths, such as Fasciola hepatica, it has been demonstrated that the SNP T687G of the Pgp gene is a good genetic marker for drug resistance in some regions of the globe but not in others [16].

On the other hand, the GG–GG genotype described by Bourguinat et al. is not expected to be the actual cause of the resistance, but rather a surrogate marker for resistance [7] and therefore it may be absent in some resistant individuals. Indeed, the sequencing of the Pgp locus in the D. immitis MP3 strain, which is considered to have at least a proportion of resistant individuals, did not reveal the GG–GG genotype neither [7]. Nevertheless, the absence of polymorphism in these Pgp loci by itself argues against the use of these SNPs as genetic markers in D. immitis of the studied area.

It must also be noted that 10 dogs from the municipality of Arraial do Cabo increased the number of microfilariae despite the use of avermectin, whereas none from Araruama presented the same increase. However, the presence of SNPs in the studied loci do not account for that difference. Since 8 years before, in that exact area of Arraial do Cabo [17] microfilaremic dogs received 3 cycles of doxycycline at 6-months intervals, it may be inferred that the antibiotic could have impacted the heartworm Wolbachia population in the area. At the former study the first treatment cycle caused a sharp drop in microfilaraemia. After 6 months microfilaraemia rebounded and the other 2 cycles did not reduce microfilaraemia as much as the first did. Since Wolbachia is known to guarantee worm vital functions [18], to colonize the worm’s hipodermis, to resist to doxycycline although in low numbers [19] and to guarantee the development of infective larvae to the adult stage [20], it may be hypothesized that doxycycline surviving Wolbachia interfere with heartworm ML uptake and, therefore, turn the worms less susceptible to the effects of ML.

Limitations

Analysis was focused in pre-described SNPs, as stated in the main text. We have not looked for other SNPs because our interest was to search for the same polymorphisms described in Lousiana, US. This objective was accomplished but we cannot exclude the presence of SNPs in other positions that would still put Pgp as a candidate gene for bearing a genetic or even functional marker of resistance. At the same time, sequencing was carried out using pools of microfilaria. As Sanger sequencing have a limit of detection of 15–20%, there is a change that we could not visualize an existing SNP if it was present in less than 20% of microfilariaes from a pool. Working in pools was a methodological choice. We preferred to enlarge our data in the number of dogs examined, instead of analysing individual microfilaria in only two or three dogs.

Abbreviations

ML:

macrocyclic lactone

SNP:

single nucleotide polymorphism

PGP:

P-glycoprotein

TUB:

β-tubulin

PCR:

polymerase chain reaction

References

  1. Wolstenholme AJ, Kaplan RM. Resistance to macrocyclic lactones. Curr Pharm Biotechnol. 2012;13(6):873–87.

    Article  CAS  PubMed  Google Scholar 

  2. Hampshire VA. Evaluation of efficacy of heartworm preventive products at the FDA. Vet Parasitol. 2005;133(2–3):191–5. https://0-doi-org.brum.beds.ac.uk/10.1016/j.vetpar.2005.04.004.

    Article  PubMed  Google Scholar 

  3. Tsiatis AC, Norris-Kirby A, Rich RG, Hafez MJ, Gocke CD, Eshleman JR, Murphy KM. Comparison of Sanger sequencing, pyrosequencing, and melting curve analysis for the detection of KRAS mutations: diagnostic and clinical implications. J Mol Diagn. 2010;12(4):425–32. https://0-doi-org.brum.beds.ac.uk/10.2353/jmoldx.2010.090188.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Wolstenholme AJ, Evans CC, Jimenez PD, Moorhead AR. The emergence of macrocyclic lactone resistance in the canine heartworm, Dirofilaria immitis. Parasitology. 2015;142(10):1249–59. https://0-doi-org.brum.beds.ac.uk/10.1017/s003118201500061x.

    Article  CAS  PubMed  Google Scholar 

  5. Bourguinat C, Keller K, Prichard RK, Geary TG. Genetic polymorphism in Dirofilaria immitis. Vet Parasitol. 2011;176(4):368–73. https://0-doi-org.brum.beds.ac.uk/10.1016/j.vetpar.2011.01.023.

    Article  CAS  PubMed  Google Scholar 

  6. Bourguinat C, Keller K, Bhan A, Peregrine A, Geary T, Prichard R. Macrocyclic lactone resistance in Dirofilaria immitis. Vet Parasitol. 2011;181(2–4):388–92. https://0-doi-org.brum.beds.ac.uk/10.1016/j.vetpar.2011.04.012.

    Article  CAS  PubMed  Google Scholar 

  7. Bourguinat C, Keller K, Blagburn B, Schenker R, Geary TG, Prichard RK. Correlation between loss of efficacy of macrocyclic lactone heartworm anthelmintics and P-glycoprotein genotype. Vet Parasitol. 2011;176(4):374–81. https://0-doi-org.brum.beds.ac.uk/10.1016/j.vetpar.2011.01.024.

    Article  CAS  PubMed  Google Scholar 

  8. Bourguinat C, Che H, Mani T, Keller K, Prichard RK. ABC-B transporter genes in Dirofilaria immitis. Int J Parasitol Drugs Drug Resist. 2016;6(2):116–24. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ijpddr.2016.04.001.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Janssen IJ, Krucken J, Demeler J, Basiaga M, Kornas S, von Samson-Himmelstjerna G. Genetic variants and increased expression of Parascaris equorum P-glycoprotein-11 in populations with decreased ivermectin susceptibility. PLoS ONE. 2013;8(4):e61635. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0061635.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Willi LM, Mendes-de-Almeida F, de Souza CD, Laeta T, Paiva JP, de Miranda MG, Knackfuss FB, Labarthe N. Serological evidence of canine exposure to arthropod-borne pathogens in different landscapes in Rio de Janeiro, Brazil. Vet Parasitol Reg Stud Rep. 2016. https://0-doi-org.brum.beds.ac.uk/10.1016/j.vprsr.2016.11.03.

    Google Scholar 

  11. Newton WL, Wright WH. The occurrence of a dog filariid other than Dirofilaria immitis in the United States. J Parasitol. 1956;42(3):246–58.

    Article  CAS  PubMed  Google Scholar 

  12. Schlottauer JC, Stromberg BE, Paul AJ, Todd KS, McCall JW, Dzimianski MT, Blagburn BL, Hendrix CM. Safety and acceptability of ivermectin in dogs with naturally acquired patent infection of Dirofilaria immitis. In: Otto GH, editor. Heartworm symposium 86, 1986, Washington, NW. Proceedings of the heartworm symposium’86. Charleston: American Heartworm Society; 1986. p. 29–35.

  13. Bourguinat C, Lee AC, Lizundia R, Blagburn BL, Liotta JL, Kraus MS, Keller K, Epe C, Letourneau L, Kleinman CL, Paterson T, Gomez EC, Montoya-Alonso JA, Smith H, Bhan A, Peregrine AS, Carmichael J, Drake J, Schenker R, Kaminsky R, Bowman DD, Geary TG, Prichard RK. Macrocyclic lactone resistance in Dirofilaria immitis: failure of heartworm preventives and investigation of genetic markers for resistance. Vet Parasitol. 2015;210(3–4):167–78. https://0-doi-org.brum.beds.ac.uk/10.1016/j.vetpar.2015.04.002.

    Article  CAS  PubMed  Google Scholar 

  14. Sambrook J, Russell DW. Preparation and analysis of eukariotic genomic DNA. In: Sambrook J, Russell DW, editors. Molecular cloning - A laboratory manual. 3rd ed. New York: Cold Spring Harbor Press; 2001, p. 69.

    Google Scholar 

  15. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.

    CAS  Google Scholar 

  16. Elliott TP, Spithill TW. The T687G SNP in a P-glycoprotein gene of Fasciola hepatica is not associated with resistance to triclabendazole in two resistant Australian populations. Mol Biochem Parasitol. 2014;198(1):45–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.molbiopara.2014.11.006.

    Article  CAS  PubMed  Google Scholar 

  17. Bendas A, Paiva J, Rossi MI, Knackfuss F, Silvano D, Mendes-de-Almeida F, Guerrero J, Labarthe N. The use of doxycycline in microfilaremic Dirofilaria immitis naturally infected dogs. Internal Journal Applied Research Veterinary Medicine. 2008;6(1):55–9.

    CAS  Google Scholar 

  18. Bandi C, Trees AJ, Brattig NW. Wolbachia in filarial nematodes: evolutionary aspects and implications for the pathogenesis and treatment of filarial diseases. Vet Parasitol. 2001;98(1–3):215–38.

    Article  CAS  PubMed  Google Scholar 

  19. Genchi C, Simoncini L, McCall JW, Venco L, Sacchi L, Bazzocchi C, Casiraghi M, Kramer L, Passeri B, Bandi C. The Wolbachia endosymbionts of Dirofilaria immitis: implications for pathology, immunology and control. In: Seward RL, Knight DH, editors. Recent advances in heartworm disease. Batavia: American Heartworm Society; 2001. p. 27–33.

    Google Scholar 

  20. McCall JW, Genchi C, Kramer L, Guerrero J, Dzimianski MT, Supakorndej P, Mansour AM, McCall SD, Supakorndej N, Grandi G, Carson B. Heartworm and Wolbachia: therapeutic implications. Vet Parasitol. 2008;158(3):204–14. https://0-doi-org.brum.beds.ac.uk/10.1016/j.vetpar.2008.09.008.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

LMVW participated in the study design, in the acquisition of data, laboratory analysis, interpretation of results and helped to draft the manuscript. NVL conceived the study and participated in the coordination of the study, the acquisition of data, interpretation of results and helped to draft the manuscript. LND participated in the interpretation of results and helped to draft the manuscript. JPP participated in the acquisition of data and performed the statistical analysis. MGNM participated in the acquisition of data and laboratory analysis. FMA participated in the study coordination, acquisition of data, laboratory analysis and helped to draft the manuscript. TZV participated in the study design and coordination, laboratory analysis and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors acknowledge Carolina Haje, Daniel Marques, Daniel Paiva, Daniel Ribeiro, Marcela Machado, Monique Paiva and Thiago Gomes for assistance in the acquisition of data. The authors are grateful to the owners and to the dogs of the study.

Competing interests

LMVW, LND, MGNM and TZV have no competing interests. FMA is a consultant for Bayer Animal Health and Idexx Laboratories in Brazil. JPP is a consultant for Zoetis and Bayer Animal Health in Brazil. NVL is a consultant for Bayer Animal Health, Idexx Laboratories, and Zoetis in Brazil.

Availability of data and materials

The datasets generated during the current study are available from the corresponding author on reasonable request.

Consent to publish

Not applicable.

Ethics approval and consent to participate

This study was approved by Fundação Oswaldo Cruz (FIOCRUZ) Institutional Animal Care and Use Committee (CEUA FIOCRUZ) under the number—LW-33/11. All procedures in dog’s were first authorized by the owner’s written consent.

Funding

The work was financed by Oswaldo Cruz Institute and personal fundings.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Tânia Zaverucha do Valle.

Additional file

Additional file 1.

Sequencing data. Chromatograms obtained by sequencing the samples of microfilariae pools of Dirofilaria immitis. Each figure shows all chromatograms for a given position of either Pgp or Tub gene.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Willi, L.M.V., Labarthe, N.V., d’Escoffier, L.N. et al. Can P-glycoprotein and β-tubulin polymorphisms be used as genetic markers of resistance in Dirofilaria immitis from Rio de Janeiro, Brazil?. BMC Res Notes 11, 152 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s13104-018-3259-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13104-018-3259-z

Keywords