Frequency of Benzimidazole Resistance in Haemonchus Contortus Populations Isolated from Buffalo, Goat and Sheep Herds

Anthelmintic resistance is an increasing problem that threatens livestock production worldwide. Understanding of the genetic basis of benzimidazole resistance recently allowed the development of promising molecular diagnostic tools. In this study, isolates of Haemonchus contortus obtained from goats, sheep and buffaloes raised in Brazil were screened for presence of the polymorphism Phe200Tyr in the β-tubulin 1 gene, which confers resistance to benzimidazole. The allelic frequency of the mutation conferring resistance ranged from 7% to 43%, and indicated that resistance to benzimidazole could be found in nematodes isolated from all the ruminant species surveyed. Although significant variation in the frequency of the F200Y mutation was observed between different herds or host species, no significant variation could be found in populations isolated from animals within the same herd. These findings suggest that screening of samples from a few animals has the potential to provide information about the benzimidazole resistance status of the entire herd, which would enable a considerable reduction in the costs of diagnosis for the producer. Molecular diagnosis has practical advantages, since it can guide the choice of anthelmintic drug that will be used, before its application in the herd, thus reducing the economic losses driven by anthelmintic resistance. Resumo A resistência aos anti-helmínticos é um problema crescente que ameaça a produção pecuária em todo o mundo. A compreensão da base genética da resistência ao benzimidazol permitiu, recentemente, o desenvolvimento de métodos diagnósticos moleculares promissores. Neste estudo, isolados de Haemonchus contortus obtidos a partir de rebanhos de caprinos, ovinos e bubalinos criados no Brasil foram avaliados quanto à presença do polimorfismo F200Y no gene da β-tubulina1, o qual confere resistência ao benzimidazol. A frequência alélica da mutação variou de 7% a 43%, indicando que a resistência ao benzimidazol pode ser encontrada em nematoides isolados a partir de todas as espécies de ruminantes pesquisadas. Embora tenha sido observada variação significativa das frequências de mutação F200Y entre


Introduction
Gastrointestinal parasites of the genus Haemonchus constitute a major cause of economic losses in livestock production in tropical and subtropical areas (AMARANTE, 2011).In Brazil, the species Haemonchus contortus prevails in small ruminants, while Haemonchus placei and Haemonchus similis mainly affect cattle (BRASIL et al., 2012).These nematodes, especially the species H. contortus, affect the development of young animals by impairing weight gain and increasing mortality.
Although nematode control was achieved through application of anthelmintics until recent years (AHID et al., 2008;SILVA et al., 2010), the widespread use of these drugs has led to emergence of resistant parasites (KAPLAN, 2004;GETACHEW et al., 2007).Therefore, researchers worldwide have been looking for alternative measures for controlling domestic animal helminthiasis.There is an emerging consensus that efficient control programs are highly dependent on knowledge of the regional particularities coupled with accurate diagnosis of the herd's degree of infection and the environmental contamination by using both clinical and laboratory examinations (EYSKER;PLOEGER, 2000).Furthermore, the efficacy of the anthelmintic chosen for use in the herd should also be proved (WALLER, 1999).
Since anthelmintic resistance has become a problem, a number of tests have been developed for assessing the anthelmintic resistance phenotype or genotype of these parasites.Even though in vivo assays such as the fecal egg count reduction test (FECRT) can easily be applied and do not require an expensive laboratory structure, they have low sensitivity and it usually takes 15 days or more for the results to be presented.The sensitivity of in vitro tests is also relatively poor, since they are considered to allow detection of resistance in isolates only when resistanceconferring alleles comprise at least 25% of the gene pool (VON SAMSON-HIMMELSTJERNA et al., 2009).Furthermore, these methods depend on coprocultures for parasite genus identification (COLES et al., 2006).
On the other hand, molecular tests based on analysis of resistance-associated target gene polymorphisms, through using PCR, are highly sensitive (SILVESTRE; HUMBERT, 2000), do not require coproculture and can complement traditional diagnostic methods.However, one major drawback relating to development of molecular tests is the current poor knowledge of the genetic mutations associated with resistant phenotypes.Nonetheless, in the case of benzimidazole (BZ) drugs, it is now well established that the resistance in trichostrongylid nematodes is mainly associated with replacement of a phenylalanine (Phe, TTC) by a tyrosine (Tyr, TAC) at position 200 (KWA et al., 1994) of the gene encoding β-tubulin isotype 1 (β-tubulin1).Although less frequent, the mutations F167Y (SILVESTRE; CABARET, 2002) and A198G (GHISI et al., 2007) have also been correlated with BZ resistance.
The objective of this study was to analyze the frequency of a BZ-resistance associated mutation in field isolates of H. contortus obtained from different domestic ruminant species in Brazil.In addition, through evaluations on the levels of intra-herd F200Y mutation frequency variation, the results presented here show that it is possible to infer the degree of contamination of the herd from the results obtained from a small group of animals, thereby creating an opportunity for reducing the costs of resistance diagnosis in the field.

Materials and Methods
Thirty feces samples were obtained directly from the rectum of animals in each of the three herds of different domestic ruminant species analyzed (total: 90 samples).Buffalo samples (Bubalus bubalis) were collected at the farm MG2 (19° 41' S; 46° 10' W), from male Murrah and hybrid buffaloes of ages ranging from 12 to 15 months.Goat samples (Capra hircus) were collected at the farm MG1 (19° 37' S; 44° 26' W), from male Crioula goats of approximately 10 months of age.Sheep samples (Ovis aries) were collected at the farm SP1 (20° 49' S; 49° 22' W), from male and female Santa Inês sheep of approximately 10 months of age.The fecal samples were collected once from each herd between November 2010 and April 2011 and the eggs per gram of feces (EPG) count was determined using the modified McMaster technique (GORDON; WHITLOCK, 1939; UENO; GONÇALVES, 1998).This study was approved by the Ethics Committee of the Federal University of Minas Gerais, Brazil, and followed the legal requirements set by the Genetic Heritage Management Council (CGEN) of the Brazilian Ministry of the Environment.Written consent was provided by the farmers whose farms were used to collect the samples.
Extensive and traditional methods of farming were practiced on all farms.Trading of live animals among farms was conducted at all these farms, but no records of direct trades among these farms were kept.Parasite control through using BZ-class drugs (albendazole) has been routinely performed at SP1 over the last 5 years, among animals that were scored and selected by means of the FAMACHA  method.On the other hand, at MG1 and MG2, BZ drugs have not been used since the establishment of the herds, more than 5 years ago.Treatments using ivermectin have also been performed on all these farms for the last 5 years, in accordance with the animal's age and EPG count and any appearance of clinical signs of infection.
Allele Specific PCR (AS-PCR): Genomic DNA was extracted from single Haemonchus sp.eggs from all the infected animals (5 goats, 6 sheep and 2 buffaloes), using an adaptation of the protocol described for larvae in Silvestre and Humbert (2000).Genotyping of the polymorphism TTC/TAC at position 200 of isotype 1 of the β-tubulin gene was performed on at least 10 specimens from each host, by means of a Haemonchus sp.specific reaction of AS-PCR, as described by Niciura et al. (2012).Plasmids containing cloned β-tubulin1 alleles bearing either the TTC200 (S allele) or the TAC200 (R allele) SNPs previously obtained by Brasil et al. (2012) were included as controls.Chi-square tests were conducted to compare allelic frequencies, using the SAEG 9.1 software (2007).

Results and Discussion
The total prevalence of helminth infection ranged from 6.6% in the buffalo herd (MG2) to 16.6% in the goat herd (MG1) and 20% in the sheep herd (SP1) (Table 1).The average EPG counts were also highest in the sheep herd (7,750), while the buffalo herd (200) and goat herd (1,000) presented lower counts.This was not unexpected, given the different susceptibilities of these ruminant species to nematode infections (KAPLAN, 2004).Sequencing of the ITS-2 region is an efficient tool for identifying nematode species (STEVENSON et al., 1995).Thus, sequences of the ITS-2 region were obtained from the DNA extracted from 30 representative specimens: 10 eggs were selected from a pool from each of the three herd isolates.All the samples analyzed were identified as H. contortus (100% similarity to GenBank sequences JQ342248 or JQ342249).These results reveal that H. contortus is probably the major nematode species infecting the ruminants tested.This finding is compatible with the helminthic fauna described by Borges et al. (2001), Veríssimo et al. (2012) and Brasil et al. (2012).Since only a fraction of the eggs sampled were identified by means of sequencing (n = 30), it is also not surprising that H. contortus appears to be a major species, given its high prolificacy and consequently its greater relative abundance of eggs, in comparison with other helminths.
AS-PCR-based genotyping using Haemonchus sp.specific primers (NICIURA et al., 2012) was performed and presented positive amplification for 182 out of 203 eggs tested, which was in accordance with the high prevalence of H. contortus found in the samples analyzed.The results presented in Table 2 show that the frequency of alleles bearing the TAC200 polymorphism in the β-tubulin isotype 1 gene (R alleles) was lower in the herds of goats (9%) and buffalo (7%), than in the sheep herd (43%).Figure 1 shows a representative result from AS-PCR, including the appropriate controls used for validation of the method.High levels of resistance to BZ have recently been reported in several herds in Brazil (CRUZ et al., 2010;NICIURA et al., 2012;VERÍSSIMO et al., 2012).Furthermore, the higher frequency of ATH treatment in the SP1 herd (4 or more treatments/year) might also explain the higher frequency of the R alleles obtained (Table 2).It is also important to highlight that resistant alleles could be found in nematodes isolated from all the three domestic ruminant species analyzed, even in the herds that were not under BZ selective pressure, i.e.MG1 and MG2.This might be a consequence of the high gene flow and low genetic structure frequently found in trichostrongylid populations, which is caused by frequent animal movements between farms, as well as by sharing of pasture areas and facilities by different ruminant species (BLOUIN et al., 1995;SILVESTRE et al., 2009;BRASIL et al., 2012).
Recent studies have focused on the distribution of the genetic variability and population structure of trichostrongylids, based on neutral markers (BLOUIN et al., 1995;BRAISHER et al., 2004;CERUTTI et al., 2010;ARCHIE;EZENWA, 2011) and genes under ATH selective pressure (SILVESTRE et al., 2009;BRASIL et al., 2012).However, data on the distribution and variation of anthelmintic resistance-conferring polymorphisms  2 demonstrate that there is no significant variation in the frequency of the F200Y SNP among animals within the same herd, irrespective of the host species analyzed.Thus, it should be possible to estimate the frequency of resistance to benzimidazole for the entire herd by sampling only a few individuals.Furthermore, use of molecular techniques for diagnosing resistance presents advantages over the traditional methods, since time-consuming steps such as coproculture are not necessary and it is possible to identify BZ resistance even at low levels (Table 2), directly from DNA extracted from nematode eggs.
Therefore, it is advisable that farmers should monitor the frequencies of any anthelmintic resistance-associated mutations present in the nematode populations that are prevalent in their herds.Von Samson-Himmelstjerna et al. (2009) reported that the critical level indicating phenotypic resistance to BZ drugs (e.g.FECRT ≤ 95% or egg hatching assays with EC50 ≥ 0.1 μg of thiabendazole/ml) corresponded to an allelic frequency ≤ 27.4 ± 5.3% for the susceptible allele (not carrying the F200Y mutation), in the H. contortus population.This would allow farmers to make rational decisions about changing the class of the drug used, as would be recommended in the case of the SP1 farm.Furthermore, rises in resistance-conferring alleles to unmanageable levels could be avoided, since molecular tests  can detect resistance before it becomes phenotypically apparent from FECRT on the nematode population.This will be especially useful in cases where anthelmintic resistance is a recessive trait, such as BZ resistance.

Conclusions
Screening for nematodes isolated from a few animals has the potential to provide information about anthelmintic resistance that is applicable to the entire herd.It is possible to considerably reduce the time and cost of benzimidazole resistance testing in the field by using the molecular methods presented in this study.

Table 1 .
Helminthic infection levels according to EPG counts.

Table 2 .
Allelic and genotypic frequencies of the F200Y polymorphism of the β-tubulin 1 gene in the MG1, MG2 and SP1 herds.Chi-square tests were used to evaluate the statistical differences among the herds' allelic frequencies.f(R): frequency of the BZ resistance-conferring allele (F200Y); f(S): frequency of the BZ sensitivity-conferring allele; RR: resistant homozygous genotype; RS: heterozygous genotype (phenotypically susceptible); SS: susceptible homozygous genotype.