Hostname: page-component-cd9895bd7-gbm5v Total loading time: 0 Render date: 2024-12-24T18:58:25.613Z Has data issue: false hasContentIssue false

Genotyping of benzimidazole resistance using β-tubulin isotype 1 marker in Haemonchus contortus of sheep and goats in Paraná, Southern Brazil

Published online by Cambridge University Press:  10 December 2024

C. Melchior do Prado*
Affiliation:
Laboratory of Veterinary Clinical Parasitology, Federal University of Parana, UFPR. Rua dos Funcionários, 1540, Cabral. CEP: 80035-050. Curitiba, PR, Brazil
J. Ferreira Vasconcelos Rodrigues
Affiliation:
Laboratory of Veterinary Clinical Parasitology, Federal University of Parana, UFPR. Rua dos Funcionários, 1540, Cabral. CEP: 80035-050. Curitiba, PR, Brazil
G.A. Frota
Affiliation:
State University of Acaraú Valley. Av. Padre Francisco S. de Araújo, 850, Alto da Brasília. CEP: 62010-295. Sobral, CE, Brazil
D.L. Vieira
Affiliation:
Laboratory of Veterinary Clinical Parasitology, Federal University of Parana, UFPR. Rua dos Funcionários, 1540, Cabral. CEP: 80035-050. Curitiba, PR, Brazil
J.P. Monteiro
Affiliation:
State University of Acaraú Valley. Av. Padre Francisco S. de Araújo, 850, Alto da Brasília. CEP: 62010-295. Sobral, CE, Brazil Embrapa Caprinos e Ovinos. Estrada Sobral-Groaíras, Km 04, s/n, Zona Rural. CEP: 62010-970. Sobral, CE, Brazil
M. Beltrão Molento*
Affiliation:
Laboratory of Veterinary Clinical Parasitology, Federal University of Parana, UFPR. Rua dos Funcionários, 1540, Cabral. CEP: 80035-050. Curitiba, PR, Brazil
*
Corresponding authors: M.B. Molento and C.M. Prado; Emails: [email protected]; [email protected]
Corresponding authors: M.B. Molento and C.M. Prado; Emails: [email protected]; [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Haemonchus contortus is one of the most pathogenic gastrointestinal parasites that infect small ruminants. The indiscriminate use of anthelmintics (i.e., benzimidazole class, BZ) to control infections has led to the reduction of drug efficacy in H. contortus populations worldwide. Resistance to BZ is associated with high frequencies of single nucleotide polymorphisms at F200Y, F167Y, and E198A positions of the β-tubulin isotype 1 gene. This study aimed to determine the frequency of single nucleotide polymorphisms associated with BZ resistance in H. contortus from 18 farms (545 sheep and 124 goats) in Paraná, Southern Brazil. Health management practices were identified as risk factors from individual farms. Genomic DNA was extracted from 20,000 larvae/farm and used in quantitative polymerase chain reaction assays for the three mutations. We ran a correlation analysis between flock health and quantitative polymerase chain reaction data. H. contortus was the most prevalent parasite in 67% (12/18) of the farms. Resistant allele frequencies were detected for F200Y (var. 46.4 to 72.0%) and F167Y (var. 15.7 to 23.8%). Only (100.0%) susceptible alleles were detected for the E198A. High treatment frequency (15/18), visual weight estimations for anthelmintic dose (15/18), no integration with other farm practices (14/18), treatment of all animals (14/18), and no quarantine period for newly acquired animals (10/18) were considered the most critical risk factors associated with BZ resistance. This is the first systematic prevalence study linking management practices on smallholder farms and the molecular data of BZ resistance of H. contortus in Southern Brazil.

Type
Research Paper
Copyright
© The Author(s), 2024. Published by Cambridge University Press

Key Findings

F200Y mutations were identified in high frequency (72%) in Haemonchus contortus.

F200Y and F167Y mutations were found even on farms with little benzimidazole use.

High treatment frequency was considered the most critical risk factor for drug resistance.

The majority of farmers use visual weight estimations to treat their animals.

Introduction

Small ruminant farming is an alternative for smallholders in low- and middle-income countries. Brazil has an estimated herd of 32 million animals, and the state of Paraná (PR) harbors about 3.0% of the sheep and 1.0% of the goat herd of the country (IBGE, 2021a). Gastrointestinal parasite infections affect the quality of life of animals, whereas Haemonchus contortus is the most prevalent pathogenic helminth (Zajac and Garza, Reference Zajac and Garza2020). Control is primarily carried out using broad-spectrum anthelmintics such as the benzimidazole (BZ) and macrocyclic lactone classes. However, inadequate practices such as the treatment of all animals in the herd, the high frequency of treatment (< 45-day interval) and drug alternation, lack of quarantine of newly arrived animals, and the use of incorrect doses are often practised by farmers and field assistants, accelerating the selection of anthelmintic resistance parasite populations (Cruz et al., Reference Cruz, Rocha, Arruda, Palieraqui, Cordeiro, Santos Junior, Molento and Santos2010; Garcia et al., Reference Garcia, Sprenger, Ortiz and Molento2016; Amaral et al., Reference Amaral, Guerra, Santana, Chicoy-Ramirez, Arenal, Lima, Alves, Molento and Faustino2021).

BZ acts by high-affinity binding to β-tubulin, a microtubule subunit protein, inhibiting its polymerisation (Lacey, Reference Lacey1988). Resistance to BZ is associated with single nucleotide polymorphisms (SNPs) at position 200 (TTC to TAC; F200Y) (Kwa et al., Reference Kwa, Veenstra and Roos1994), 167 (TTC to TAC; F167Y) (Silvestre and Cabaret, Reference Silvestre and Cabaret2002) and 198 (GAA to GCA; E198A) (Ghisi et al., Reference Ghisi, Kaminsky and Maser2007) in the H. contortus β-tubulin isotype 1 gene. In Brazil, a few studies have used molecular diagnosis to determine resistance to BZ. Niciura et al. (Reference Niciura, Veríssimo, Gromboni, Rocha, Mello, Barbosa, Chiebao, Cardoso, Silva, Otsuk, Pereira, Ambrosio, Nardon, Ueno and Molento2012) conducted a survey using nested-PCR and ARMS-PCR to assess the frequency of the F200Y in 33 sheep farms in São Paulo, Brazil. The study found frequencies of resistant homozygous (RR) genotypes to vary from 0.1 to 66.7% and resistance alleles from 9.0 to 74.0%. In the same year, Brasil et al. (Reference Brasil, Nunes, Bastianetto, Drummond, Carvalho, Leite, Molento and Oliveira2012) evaluated the variation of the F200Y, F167Y, and E198A in cattle, buffalo, sheep, and goat herds using nested PCR. Using the AS-PCR technique, Nunes et al. (Reference Nunes, Santos, Bastianetto, Oliveira and Brasil2013) evaluated only the frequency of F200Y in H. contortus in buffaloes, goats, and sheep in Minas Gerais, Brazil. Moreover, Santos et al. (Reference Santos, Monteiro, Ribeiro, Macedo, Camurca-Vasconcelos, Vieira and Bevilaqua2014) reported the first study using quantitative real-time polymerase chain reaction (qPCR) to detect BZ mutations. Chagas et al. (Reference Chagas, Sampaio Junior, Pacheco, Cunha, Cruz, Scofield and Goés-Cavalcante2016) also determined a high frequency (31%) of F200Y-resistant alleles in 305 specimens of H. contortus from 12 sheep farms in Pará, Eastern Amazon, Brazil.

Santos et al. (Reference Santos, Vasconcelos, Frota, Ribeiro, André, Vieria, Teixeira, Bevilaqua and Monteiro2017a) used qPCR to detect all three SNPs after the regular use of oxfendazole, ivermectin (IVM), and their combination. Sheep were also experimentally infected with the H. contortus inbred-susceptible-Edinburgh (ISE) isolate, followed by treatment. The isolates resulting from oxfendazole and IVM selection showed a high frequency of F167Y (85.5 and 63%, respectively). In the same year, as Santos et al. (Reference Santos, Monteiro, Ribeiro, Macedo, Araújo-Filho, Andre, Araújo, Vasconcelos, Freitas, Camurça-Vasconcelos, Vieira and Bevilaqua2017b) expanded the surveillance for BZ resistance, they reported the E198A allele mutation for the first time in Brazil. The authors highlighted that testing all three β-tubulin polymorphisms is vital when determining BZ resistance. Lambert et al. (Reference Lambert, Nishi, Mendonça, Souza, Julião, Gusmão and Almeida2017) found the frequencies of resistant alleles for F167Y to be 50% and F200Y 42% in H. contortus from goat herds in Bahia. Amaral et al. (Reference Amaral, Guerra, Santana, Chicoy-Ramirez, Arenal, Lima, Alves, Molento and Faustino2021) determined the mutation frequency at F200Y in adult worms from seven H. contortus populations in Pernambuco.

In Paraná, studies conducted by Thomaz-Soccol et al. (Reference Thomaz-Soccol, Sotomaior, Souza, Castro, Pessôa Silva and Milezewski1996 and Reference Thomaz-Soccol, Souza, Sotomaior, Castro, Milczewski, Mocelin and Pessôa Silva2004) and Cunha Filho et al. (Reference Cunha Filho, Pereira and Yahamamura1998) evaluated the efficacy of anthelmintics in 6, 42, and 10 sheep flocks, respectively, using the faecal egg count reduction test. In all studies, different levels of drug resistance were detected. Cunha Filho et al. (Reference Cunha Filho, Pereira and Yahamamura1998) reported that in 80% of the farms, animals did not rotate pastures and had an average of eight treatments/year. However, the molecular mechanisms associated with BZ resistance have never been investigated in sheep and goat herds in this state.

qPCR is an advancement of conventional PCR and is one of the most widely used clinical techniques for diagnosing diseases, including H. contortus infection. Many other techniques, such as digital PCR, were developed from this technique. Furthermore, PCR is indispensable to exploring BZ resistance in modern methods, such as next-generation sequencing (Baltrušis et al., Reference Baltrušis, Halvarsson and Höglund2018). Even in the face of more current methodologies such as nemabiome (Avramenko et al., Reference Avramenko, Redman, Lewis, Yazwinski, Wasmuth and Gilleard2015), qPCR allows the detection of SNP with high sensitivity and specificity at 1% in parasite populations (Elmahalawy et al., Reference Elmahalawy, Halvarsson, Skarin and Höglund2018; Costa-Junior et al., Reference Costa-Junior, Chaudhry, Silva, Sousa, Silva, Cutrim-Júnior, Brito and Sargison2021). The technique can also be helpful when associated with herd health practices in preventing parasite infection and delaying parasite selection (Veríssimo et al., Reference Veríssimo, Niciura, Alberti, Rodrigues, Barbosa, Chiebao, Cardoso, Silva, Pereira, Margatho, Costa, Nardon, Ueno, Curci and Molento2012).

Although qPCR has already been used in Brazil, the link between SNP frequencies and parasite control practices has never been thoroughly addressed. The objectives of this work are twofold: 1. to determine the presence of polymorphisms associated with BZ resistance in H. contortus isolated from sheep and goats and 2. to identify management practices that may be related to the development of resistance to BZ.

Material and methods

Sampling and study area

Faecal samples were collected from 545 sheep and 124 goats (Boer goats, and Dorper, Texel, Santa Inês, and Île de France sheep) from 18 farms. The animals were distributed in eight municipalities (Foz do Iguaçu, FI; Matelândia, MA; Medianeira, ME; Ramilândia, RM; Santa Terezinha de Itaipu, ST; São Miguel do Iguaçu, SM; Serranópolis do Iguaçu, SI, and Vera Cruz do Oeste, VC) in Western Paraná state (PR) (Figure 1). Farms were randomly selected from the regional office of the Secretary of Agriculture databank. The inclusion criteria were no anthelmintic treatment for at least 100 days before the sampling.

Figure 1. Location of the municipalities with sheep and goat farms in Western Paraná, Brazil. Obs. Corresponding municipalities: 1. Foz do Iguaçu, 2. Matelândia, 3. Medianeira, 4. Ramilândia, 5. Santa Terezinha de Itaipu, 6. São Miguel do Iguaçu, 7. Serranópolis do Iguaçu, 8. Vera Cruz do Oeste.

Table 1 describes the area of each municipality, the number of inhabitants, and the population of sheep and goats. The region has a subtropical climate with an average annual temperature of 24°C and two distinct seasons, a hot and humid summer and a dry and cold winter, according to Köppen-Geiger (Aparecido et al., Reference Aparecido, Rolim, Richetti, Souza and Johann2016). All municipalities are included in the Atlantic Forest biome (IBGE, 2021b). Except for FI, where the main economic activity is tourism, the other cities have their economy based on the agribusiness of cattle, soybean, and corn plantations (IBGE, 2021a).

Table 1. Municipalities included in the study, total area (km2), population size, main economic activity, and the number of sheep (S) and goats (G) of Western Parana, Brazil

Farm management and health questionnaire

All farmers answered a questionnaire about management and parasite control practices based on Niciura et al. (Reference Niciura, Veríssimo, Gromboni, Rocha, Mello, Barbosa, Chiebao, Cardoso, Silva, Otsuk, Pereira, Ambrosio, Nardon, Ueno and Molento2012). The information included the time the farms have been rearing sheep and goats, total number of animals, breed, whether the pasture is shared with other animals (horses or cattle), and the procedure for handling newly acquired animals (treated and quarantined for a few days, quarantined for a few days, treated and introduced without quarantine, introduced without treatment or quarantine), treatment practices (all herd, after FEC examination, or clinical signs/FAMACHA), frequency of treatments, anthelmintic rotation, history of anthelmintics used, and how the doses of anthelmintics were determined for treatment.

Faecal Sampling and Coproculture

Faeces were collected from at least 50% of the adult animals in each herd. In herds with more than 100 animals, at least 50 animals were sampled. The samples were labeled by the farm’s city name and number and used in the coproculture to obtain third-stage larvae (L3) (Roberts and O’Sullivan, Reference Roberts and O’Sullivan1950). A pooled sample from each property was placed in a 500-ml glass flask, and vermiculite was added and mixed with a spatula to obtain a homogeneous compound. The mixture was humidified, and the flask was covered with a Petri dish and string to allow aeration. The coproculture occurred at room temperature (≥ 22°C) for 10 days. After this period, the flask was filled with warm water (28°C) and inverted over a Petri dish. After 6 h, the L3 was taken with a Pasteur pipette and stored in a 12-ml Falcon tube at 4°C for identification and molecular test.

Morphological identification and quantification of L3

The L3 suspension was centrifuged at 6000×g for 2 min. The supernatant was carefully discarded, and the remaining liquid was homogenized. A 10-μl aliquot of the solution was pipetted onto a microscope slide with 10 μl of Lugol, and a cover slip was added. Genus and species were identified for a minimum of 100 L3. L3 of H. contortus and H. placei have differences of an average length of 730 and 820 μm, tail length of 75 and 100 μm, and the percentage of the filament of the tail sheath extension of 10 and 20 μm, respectively (van Wyk and Mayhew, Reference Van Wyk and Mayhew2013). A pooled sample of each property that presented > 60% of H. contortus was selected for the molecular studies.

DNA Extraction

Genomic DNA was extracted from L3 using a protocol standardized by Santos et al. (Reference Santos, Monteiro, Ribeiro, Macedo, Camurca-Vasconcelos, Vieira and Bevilaqua2014). Larvae (20,000 L3/sample) from each farm were disrupted using 1.0-mm zirconia/silica beads in lysis buffer (0.2% SDS; 50 mM EDTA; 50 mM Tris- HCl, pH 8.0, RNAse at 100 μg/ml and Proteinase K at 0.4 mg/ml) on a BioSpec Mini-BeadBeater-16 (Bartlesville, USA) for 2 min (four 30-s cycles). The lysate was transferred to a new tube and incubated in a dry bath at 56°C for 1 h, followed by adding 5 M potassium acetate and incubating on ice for 2 h. Samples were centrifuged, and the supernatant was transferred to a new tube for DNA precipitation with isopropanol. After another round of centrifugation, the supernatant was discarded, leaving only the DNA pellet. The DNA samples were resuspended in 500 μl TE (1 mM EDTA, 10 mM Tris- HCl, pH, 8.0) and underwent phenol/chloroform purification twice for protein removal. The aqueous phase was collected after centrifugation, and DNA was precipitated with isopropanol, centrifuged, and, after discarding the supernatant, the material was washed with 70% ethanol. The pellet was dried in a dry bath at 37°C. Between 20 to 200 μl of TE (1 mM EDTA, 10 mM Tris-HCl, pH, 8.0) was added and homogenized for DNA resuspension. The samples were stored in a freezer at -20°C until use. The amount of extracted DNA and the presence of protein contaminants were determined by UV spectrophotometry using NanoDrop Thermo Fisher Scientific (Waltham, USA) at 260 and 280 nm, respectively. All centrifugation steps were done at 13,000×g for 15 min. at 4oC.

Quantitative real-time PCR

All tests were performed in triplicates with specific primers for the F200Y, F167Y, and E198A polymorphisms in the β-tubulin isotype 1 gene of H. contortus, according to Santos et al. (Reference Santos, Monteiro, Ribeiro, Macedo, Camurca-Vasconcelos, Vieira and Bevilaqua2014). It is noteworthy that even though the primers were designed for H. contortus, they could probably work for H. placei, as both species interbreed. Quantitative PCR assays were carried out with 12.5 μl 2x Fast Start Universal SYBR Green Master Mix (Roche, West Sussex, UK), 0.3 pMol/μl of each primer (forward and reverse), 25 ng of DNA and distilled water to complete a total volume of 25 μl. Negative controls were performed with distilled water instead of DNA, and genomic DNA from the ISE isolate (Inbred-Susceptible-Edinburgh) was used as a reference for susceptibility. Amplification conditions for F200Y and F167Y were 95°C for 10 min and 35 cycles at 95°C for 15 s and 58°C for 30 s. The amplification of the E198A used only 34 cycles under similar conditions. The dissociation curve was analyzed to detect primer dimers (Santos et al., Reference Santos, Monteiro, Ribeiro, Macedo, Camurca-Vasconcelos, Vieira and Bevilaqua2014).

Analysis of qPCR Data

The cycle in which the detected fluorescence was above the background noise (Cq) for each qPCR reaction was determined by the Eppendorf Realplex 2.2 software (Hamburg, Germany). The allele frequencies were estimated using a previously described formula (Germer et al., Reference Germer, Holland and Higuchi2000), in which the frequency of the resistant allele F = [1/(2ΔCt + 1)]×100, where ΔCt = (Ct resistant allele - Ct susceptible allele). In cases where it was not possible to determine the Cq value for an allele, its frequency was considered zero, and the frequency of its detected counterpart was set to approximately 100%. Parasite populations were considered resistant when the percentage of resistant H. contortus was greater than 10%. Estimations of resistant and susceptible parasites were made assuming that the populations were in the Hardy-Weinberg equilibrium (Barrère et al., Reference Barrère, Keller, Von Samson-Himmelstjerna and Prichard2013).

Descriptive statistics and correlation of management protocols and qPCR

Initially, a descriptive statistical analysis was carried out to visualise the farm data regarding parasite resistance and mutation frequency. The data were plotted for F200Y and F167Y mutations with the frequency of treatments and farm practices. Levene’s test first analysed the association between the frequency of resistant alleles and management to determine whether the relationship between the health variables and the frequency of the SNP had a similar variance (P > 0.05) (https://datatab.net). The Shapiro-Wilk test was performed to determine if the data had a normal distribution (P > 0.05). The one-way analysis of variance ANOVA was performed to verify if there was a difference between the means of the variables (P ≤ 0.05). Statistical analysis was performed using the R software version i388 4.0.2 (www.r-project.org).

Results

Farm Management and Health Questionnaire

The average time the farmers engaged in the activity was 8.5 years, with a herd size averaging 94 animals. Half of the farms raised crossbred sheep, and the other half raised South African Boer goats. Most farms were semi-intensive, and half adopted grazing rotation based on pasture availability. Levamisole and albendazole were used in 55.5 and 38.8% of the farms, respectively. Monepantel and closantel were reported in 5.55 and 33.3%, respectively.

Regarding macrocyclic lactones, doramectin, moxidectin (MOX), and IVM were reported in 33.3, 27.8, and 16.7% of the farms, respectively. Cattle and horses had sporadic interactions with sheep and goats. The FAMACHA method was used by 22.2% of the farmers. Farmers did not use any laboratory method to decide on parasite control (Table 2).

Table 2. Frequency (%) of health measures adopted by 18 sheep and goat farmers in Western Paraná, Brazil

Morphological identification and quantification of L3

H. contortus spp. (66.6%) and Trichostrongylus spp. (45.5%) were the most abundant parasites. Infections exclusively by H. contortus or Trichostrongylus spp. were detected in six (33.3%) and five (27.8%) farms, respectively. Mixed infections of H. contortus and Trichostrongylus spp. were detected in four (22.2%) farms, and H. contortus and Cooperia spp. were found in one (5.55%) farm. One (5.55%) farm had a mixed infection of H. contortus, Trichostrongylus spp., and Teladorsagia sp. In all mixed infections, H. contortus had the highest frequency. H. placei was not found in any samples. No L3 were recovered in the municipality of Vera Cruz do Oeste.

Analysis of qPCR data

The qPCR results showed resistant allelic frequencies of up to 72.0 and 23.8% for F200Y and F167Y, respectively (Table 3). The isolates from ST4 (farm 4 of Santa Terezinha de Itaipu) and SI2 (farm 2 of Serranópolis do Iguaçu) presented the highest frequencies of resistant alleles of F200Y. The SI1 and RM1 showed the highest frequencies of resistant F167Y alleles, with 23.8 and 20.7%, respectively. All studied populations showed frequencies of resistant nematodes above 10% for F200Y, with allelic frequencies ranging from 46.4 to 72.0%. Nine populations showed resistant allelic frequencies above 10% for F167Y, varying from 15.7 to 23.8%. The E198A resulted in the detection of susceptible alleles for all isolates. The ISE isolate showed a frequency of 99.7% of susceptible alleles at F200Y and a 100% frequency of susceptible alleles for the F167Y and E198A alleles (Table 3). The dissociation curve analysis showed one well-defined peak per reaction, indicating that only one PCR product was amplified (data not shown).

Table 3. Average qPCR and cycle quantification (Cq) and allelic frequencies of the F200Y, F167Y, and E198A mutations from Haemonchus contortus of sheep and goat farms in Western Paraná, Brazil

a Municipalities: FI, Foz do Iguaçu; MA, Matelândia; ME, Medianeira; RM, Ramilândia; ST, Santa Terezinha de Itaipu; SM, São Miguel do Iguaçu; SI, Serranópolis do Iguaçu.

b Average Cq (cycle quantification) of the triplicates.

c Susceptible.

d Resistant.

e In cases where it was not possible to determine the Cq value for a given allele, its frequency was zeroed, and the frequency of its detected counterpart was approximately 100%.

Descriptive statistics and correlation of management protocols and qPCR

All examined populations had the F200Y resistant allelic frequencies above 40%, averaging 54.3%. The data revealed that all farms had more than 10% resistant individuals, considered genetically resistant when using the study criteria. Regarding F167Y, farms FA1, ST4, and SI2 showed very low resistant allelic frequencies (< 10%). The only factor these farms had in common was that all animals were mixed breed. This characteristic is noteworthy as the region is hot (> 30-34oC) and dry most of the year. In addition, nine farms showed low frequencies of resistant F167Y, greater than 10%, with a mean of 18.5% but not higher than 25%. It was found that even in farms where BZ was not regularly used, the F200Y was found in high frequencies (Table 4). The frequency of the SNP E198A was not related to any farm management, as only susceptible alleles were found. The correlation between parasite control protocols and the frequency of F200Y and F167Y mutations was not statistically significant (P > 0.05). The comparison of SNP frequencies grouped by general management variables was not statistically significant (P ≥ 0.05).

Table 4. Allelic frequencies (%) of the F200Y, F167Y, and E198A mutations of Haemonchus contortus in sheep and goat farms (isolates) with and without anthelmintic rotation (yes/no), in Paraná, Brazil

ª Municipality/isolate: FI: Foz do Iguaçu, MA: Matelândia, ME: Medianeira, RM: Ramilândia, ST: Santa Terezinha de Itaipu, SM: São Miguel do Iguaçu, SI: Serranópolis do Iguaçu. S: Susceptible, R: Resistant.

* Farms where BZ was not used.

Discussion

The present study reports the molecular detection of BZ resistance in H. contortus of small ruminants in Western Parana, Brazil, for the first time. The qPCR revealed allele frequencies of β-tubulin isotype 1 for the SNPs F200Y and F167Y related to BZ resistance, using the ISE isolate as a control population. The frequency of SNP F200Y resistance in our study was higher than that reported in previous studies in Brazil (Niciura et al., Reference Niciura, Veríssimo, Gromboni, Rocha, Mello, Barbosa, Chiebao, Cardoso, Silva, Otsuk, Pereira, Ambrosio, Nardon, Ueno and Molento2012; Nunes et al., Reference Nunes, Santos, Bastianetto, Oliveira and Brasil2013; Chagas et al., Reference Chagas, Sampaio Junior, Pacheco, Cunha, Cruz, Scofield and Goés-Cavalcante2016; Santos et al., Reference Santos, Monteiro, Ribeiro, Macedo, Araújo-Filho, Andre, Araújo, Vasconcelos, Freitas, Camurça-Vasconcelos, Vieira and Bevilaqua2017b; Lambert et al., Reference Lambert, Nishi, Mendonça, Souza, Julião, Gusmão and Almeida2017; Fávero et al., Reference Fávero, Santos, Araújo, Ramunke, Krucken, Von Samson-Himmelstjerna and Borges2020). This may be related to climate, management practices, and study time. The climate in Paraná is characterized by hot and humid summers and mild winters, which favor the survival of L3 during most of the year. In contrast, Pará, Bahia, and Minas Gerais have dry and hot weather, with periods that are less optimal for the survival of free-living stages. Consequently, animals in Paraná may be treated more often, increasing the selection for resistant alleles. Moreover, some of the previous data were collected up to 10 years before the present field study, reflecting the actual selection process in Parana.

No management practice was significantly associated with high SNP F200Y and SNP F167Y frequencies in the present study. The data reflect the literature, as concrete high-risk factors are limited when correlating health management and parasite control strategies (Niciura et al., Reference Niciura, Veríssimo, Gromboni, Rocha, Mello, Barbosa, Chiebao, Cardoso, Silva, Otsuk, Pereira, Ambrosio, Nardon, Ueno and Molento2012; Odoi et al., Reference Odoi, Gathuma, Gachuiri and Omore2007; Reynecke et al., Reference Reynecke, van Wyk, Gummow, Dorny and Boomker2011). Chagas et al. (Reference Chagas, Sampaio Junior, Pacheco, Cunha, Cruz, Scofield and Goés-Cavalcante2016) found no significant association with health protocols in sheep infected with H. contortus. Similarly, Lambert et al. (Reference Lambert, Nishi, Mendonça, Souza, Julião, Gusmão and Almeida2017) found no association between health practices and the frequency of resistant alleles in H. contortus in goat farms. Reynecke et al. (Reference Reynecke, van Wyk, Gummow, Dorny and Boomker2011) used a multiple regression algorithm to simulate worm burdens in sheep naturally infected with H. contortus during the summer (rainy season) in South Africa. The authors found that the most critical risk factor associated with high worm burden was sex, which is higher in rams than in ewes. A protective factor found in this study was the early detection and deworming of stragglers when using the FAMACHA method. Farmers are expected to adopt sustainable practices to preserve refugia populations, leaving healthy categories (i.e., adult ewes) without anthelmintics (Molento et al., Reference Molento, van Wyk and Coles2004; Vilela et al., Reference Vilela, Bezerra, Bezerra, Dantas, Alcântara, Oliveira, Nóbrega, Calazans, Feitosa, Braga and Molento2021).

Most farmers treat their animals based on visual weight estimation (15/18). The producers described this practice as facilitating management, mainly in farms with large herds (>200 animals). As a result, administering high doses of anthelmintics is essential in increasing parasite selection pressure (Molento, Reference Molento, van Wyk and Coles2004) and may also influence the frequency of resistant mutations in these herds. We have observed a high allele frequency of SNPs and the lack of quarantine in most farms (55%, 10/18). Therefore, the high prevalence of anthelmintic resistance may result from introducing newly acquired animals to the flock without significant health control. Although the recommendation of quarantine of animals has been widely advocated (Leathwick et al., Reference Leathwick, Hosking, Bisset and Mckay2009), this practice needs to be frequently addressed in the field, as reported here. According to Thomaz-Soccol et al. (Reference Thomaz-Soccol, Sotomaior, Souza, Castro, Pessôa Silva and Milezewski1996), from 1989 to 1994, more than 200,000 heads of sheep were imported from Uruguay and Rio Grande do Sul and were distributed to farmers in Parana as part of an agricultural program. This significant animal movement may have resulted in the introduction of resistant parasite populations (Echevarria et al., Reference Echevarria, Armour and Duncan1991) that could have changed the genetic background in the region.

Concerning treatment strategy and the frequency of BZ mutations, it is noteworthy that even in farms where BZ was not used regularly, the SNP F200Y was found in high frequencies. This data may result from animal movement between farms through sale or reproduction spreading BZ-resistant and cross-resistance (BZ and macrocyclic lactone) parasites, as there is evidence that ivermectin interacts directly with H. contortus tubulin (Ashraf et al., Reference Ashraf, Mani, Beech and Prichard2015). Following this, Mottier and Prichard et al. (Reference Mottier and Prichard2008) reported that IVM and MOX alter the frequency of β-tubulin alleles in H. contortus. Santos et al. (2017) also found corresponding results in sheep treated with IVM and increased allele frequency of F200Y and F167Y in H. contortus. Odoi et al. (Reference Odoi, Gathuma, Gachuiri and Omore2007) have determined risk factors in sheep and goat farms in Kenya. The authors reported the farmer’s education level, the animals’ age category, and the deworming during the preceding month as significant predictors of changing FEC. Niciura et al. (Reference Niciura, Veríssimo, Gromboni, Rocha, Mello, Barbosa, Chiebao, Cardoso, Silva, Otsuk, Pereira, Ambrosio, Nardon, Ueno and Molento2012) have also found that inexperienced farmers, sheep breeds, lack of pasture rotation, entire herd treatment, short drug rotation, and visual dose calculation were associated with the high frequency of resistant genotypes.

The high frequencies observed for F200Y (100%) and F167Y (75%) and how they differ from other regions within the country and worldwide are most likely due to the high genetic diversity observed in H. contortus (Parvin et al., Reference Parvin, Dey, Shohana, Anisuzzaman and Alam2024). We are also reporting a high frequency of concomitant (double) mutations in Brazil. The dispersion of such common mutations could have occurred after local fairs increasing animal movement from sales. Although we did not find the SNP E198A, Zhang et al. (Reference Zhang, Gasser, Yang, Yin, Zhao, Bao, Pan, Huang, Wang, Zou, Zhou, Zhao, Fang and Hu2016) evaluating the three mutations in 192 worms from eight provinces in China, found no F167Y. Barrère et al. (Reference Barrère, Alvarez, Suarez, Ceballos, Moreno, Lanusse and Prichard2012), in a study of 24 artificially infected sheep with H. contortus, found no parasites with more than two mutations. Kotze et al. (Reference Kotze, Cowling, Bagnall, Hines, Ruffell, Hunt and Coleman2012) have determined an exciting shift, showing that the SNP F200Y and E198A started at approximately the same level on the Wallangra isolate, but when increasing thiabendazole concentrations, the E198A frequency increased as the F200Y decreased. Ghisi et al. (Reference Ghisi, Kaminsky and Maser2007) studied 25 β-tubulin alleles cloned and sequenced from H. contortus isolates. All but one isolate carried either the SNP F200Y or the SNP E198A, suggesting that parasites with both mutations are uncommon, differing from our data.

Molecular detection methods contribute to the validation of phenotypic studies (in vitro and in vivo). In the present work, the F200Y allele was the most common polymorphism found in H. contortus. However, the presence of the SNP F167Y can also be considered a strong indication of the lack of BZ efficacy. We are concerned about the reported scenario since the use of BZ drugs was confirmed in 55.5% of the farms. These farms may have used BZ with no signs of reversion (genetic dilution) or loss of resistance mutations (Fissiha and Kinde, Reference Fissiha and Kinde2021). The data also represent the intense animal movement amplifying BZ resistance alleles. Further studies are needed using new technologies for resistance diagnostics (i.e., nemabiome metabarcoding). This technology allows in-depth analysis of the composition of gastrointestinal nematode populations and the prevalence of drug-resistance alleles (Avramenko et al., Reference Avramenko, Redman, Lewis, Yazwinski, Wasmuth and Gilleard2015, Reference Avramenko, Redman, Lewis, Bichuette, Palmeira, Yazwinski and Gilleard2017; Workentine et al., Reference Workentine, Chen, Zhu, Gavriliuc, Shaw, Rijke, Redman, Avramenko, Wit, Poissant and Gilleard2020; Costa-Junior et al., Reference Costa-Junior, Chaudhry, Silva, Sousa, Silva, Cutrim-Júnior, Brito and Sargison2021).

Routine molecular diagnostics of an individual or pooled larvae may help determine the population genetic distribution associated with drug failure. Although isolating one condition that would impose parasite selection is extremely difficult, our data indicate that studies with larger sample sizes and more extended periods of observation (i.e., > three years) would be necessary to establish proper risk factors in any region. The diffusion and implementation of more sustainable parasite control strategies would impact the health and welfare of more than 10,000 sheep and 1,800 goats in the area.

Acknowledgements

The authors thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, CAPES, Brazil, for providing an M.Sc. scholarship to C.M. Prado. The authors also thank Embrapa Caprinos e Ovinos for all the technical support and infrastructure.

Author contribution

C.M.P., M.B.M., and J.P.M. conceived and designed the study. C.M.P. performed the data collection. C.M.P., J.F.V., and G.A.F. performed the molecular experiments. C.M.P. performed statistical analyses. C.M.P., M.B.M., and J.P.M. wrote the article.

Financial support

This research received no specific grant from funding agencies.

Competing interest

The authors declare there are no conflicts of interest.

Ethical standard

This study was approved by the Animal Use Ethics Committee of the Agricultural Sciences Campus of the Federal University of Parana (protocol number 014/2020).

References

Amaral, R.L.G., Guerra, N.R., Santana, I.M., Chicoy-Ramirez, Y., Arenal, A., Lima, M.M., Alves, L.C., Molento, M.B., and Faustino, M.A.G. (2021) Molecular diagnostic of Haemonchus contortus resistance to the benzimidazole group in sheep and goats in the state of Pernambuco, Brazil. Archives of Veterinary Science 26(3), 5971.Google Scholar
Aparecido, L.E.O., Rolim, G.S., Richetti, J., Souza, P.S., and Johann, J.A. (2016) Köppen, Thornthwaite and Camargo climate classifications for climatic zoning in the State of Paraná, Brazil. Ciência e Agrotecnologia 40, 405417. Doi: 10.1590/1413-70542016404003916.CrossRefGoogle Scholar
Ashraf, S., Mani, T., Beech, R., and Prichard, R. (2015) Macrocyclic lactones and their relationship to the SNPs related to benzimidazole resistance. Molecular and Biochemical Parasitology 201(2), 128134. doi: 10.1016/j.molbiopara.2015.07.007CrossRefGoogle Scholar
Avramenko, R. W., Redman, E. M., Lewis, R., Yazwinski, T. A., Wasmuth, J.D., and Gilleard, J.S. (2015). Exploring the gastrointestinal “nemabiome”: deep amplicon sequencing to quantify the species composition of parasitic nematode communities. PLoS One 10(12), e0143559. Doi: 10.1371/journal.pone.0143559CrossRefGoogle ScholarPubMed
Avramenko, R.W., Redman, E.M., Lewis, R., Bichuette, M.A., Palmeira, B.M., Yazwinski, T.A., and Gilleard, J.S. (2017). The use of nemabiome metabarcoding to explore gastro-intestinal nematode species diversity and anthelmintic treatment effectiveness in beef calves. International Journal for Parasitology 47(13), 893902. doi: 10.1016/j.ijpara.2017.06.006CrossRefGoogle ScholarPubMed
Baltrušis, P., Halvarsson, P., and Höglund, J. (2018) Exploring benzimidazole resistance in Haemonchus contortus by next-generation sequencing and droplet digital PCR. International Journal for Parasitology: Drugs and Drug Resistance 8, 411419. doi: 10.1016/j.ijpddr.2018.09.003Google ScholarPubMed
Barrère, V., Alvarez, L., Suarez, G., Ceballos, L., Moreno, L., Lanusse, C., and Prichard, R.K. (2012) Relationship between increased albendazole systemic exposure and changes in single nucleotide polymorphisms on the β-tubulin isotype 1 encoding gene in Haemonchus contortus. Veterinary Parasitology 186, 344349. doi: 10.1016/j.vetpar.2011.11.068CrossRefGoogle ScholarPubMed
Barrère, V., Keller, K., Von Samson-Himmelstjerna, G., and Prichard, R.K. (2013) Efficiency of a genetic test to detect benzimidazole resistant Haemonchus contortus nematodes in sheep farms in Quebec, Canada. Parasitology International 62, 464470. doi: 10.1016/j.parint.2013.06.001CrossRefGoogle ScholarPubMed
Brasil, B.S.A.F., Nunes, R.L.A., Bastianetto, E.B., Drummond, M.G., Carvalho, D.C., Leite, R.C., Molento, M.B., and Oliveira, D.A.A. (2012) Genetic diversity patterns of Haemonchus placei and Haemonchus contortus populations isolated from domestic ruminants in Brazil. International Journal for Parasitology: Drugs and Drug Resistance 42, 469479. doi: 10.1016/j.ijpara.2012.03.003CrossRefGoogle ScholarPubMed
Chagas, A.M., Sampaio Junior, F.D., Pacheco, A., Cunha, A.B., Cruz, J.S., Scofield, A., and Goés-Cavalcante, G. (2016) F2OOY polymorphism of the β-tubulin isotype 1 gene in Haemonchus contortus and sheep flock management practices related to anthelmintic resistance in eastern Amazon. Veterinary Parasitology 226, 104108. doi: 10.1016/j.vetpar.2016.06.038CrossRefGoogle ScholarPubMed
Costa-Junior, L.M., Chaudhry, U.N., Silva, C.R., Sousa, D.M., Silva, N.C., Cutrim-Júnior, J.A., Brito, D.R., and Sargison, N.D. (2021). Nemabiome metabarcoding reveals differences between gastrointestinal nematode species infecting co-grazed sheep and goats. Veterinary Parasitology 289, 109339. doi: 10.1016/j.vetpar.2020.109339CrossRefGoogle ScholarPubMed
Cruz, D.G., Rocha, L.O., Arruda, S.S., Palieraqui, J.G.B., Cordeiro, R.C., Santos Junior, E., Molento, M.B., and Santos, C.P. (2010) Anthelmintic efficacy and management practices in sheep farms from the state of Rio de Janeiro, Brazil. Veterinary Parasitology 170, 340343. doi:10.1016/j.vetpar.2010.02.030.CrossRefGoogle ScholarPubMed
Cunha Filho, L.F.C., Pereira, A.B.L., and Yahamamura, M.H. (1998) Resistance anthelmintic in sheep in the region of Londrina, Parana State, Brazil. Semina Ciências Agrarias 19, 3137. doi: 10.5433/1679-0359.1998v19n1p31CrossRefGoogle Scholar
Echevarria, F.A.M., Armour, J., and Duncan, J.L. (1991) Efficacy of some anthelmintics on an ivermectin-resistant strain of Haemonchus contortus in sheep. Veterinary Parasitology 39, 279284. doi: 10.1016/0304-4017(91)90044-VCrossRefGoogle Scholar
Elmahalawy, S.T., Halvarsson, P., Skarin, M., and Höglund, J. (2018) Genetic variants in dyf-7 validated by droplet digital PCR are not drivers for ivermectin resistance in Haemonchus contortus. International Journal for Parasitology: Drugs and Drug Resistance 8, 278286. doi: 10.1016/j.ijpddr.2018.04.005.Google Scholar
Fávero, F.C., Santos, L.B., Araújo, F., Ramunke, S., Krucken, J., Von Samson-Himmelstjerna, G., and Borges, F.A. (2020) Haemonchus sp. in beef cattle in Brazil: species composition and frequency of benzimidazole resistance alleles. Preventive Veterinary Medicine 185, 105162. doi: 10.1016/j.prevetmed.2020.105162CrossRefGoogle ScholarPubMed
Fissiha, W. and Kinde, M. Z. (2021). Anthelmintic resistance and its mechanism: a review. Infection and Drug Resistance, 54035410. doi: 10.2147/IDR.S332378CrossRefGoogle Scholar
Garcia, C., Sprenger, L.K., Ortiz, E.B., and Molento, M.B. First report of multiple anthelmintic resistance in nematodes of sheep in Colombia. Annals of the Brazilian Academy of Sciences, 88(1), 397402. 2016. doi: 10.1590/0001-3765201620140360.CrossRefGoogle ScholarPubMed
Germer, S., Holland, M.J., and Higuchi, R. (2000) High-throughput SNP allele-frequency determination in pooled DNA samples by kinetic PCR. Genome Research 10, 258266. doi: 10.1101/gr.10.2.258CrossRefGoogle ScholarPubMed
Ghisi, M., Kaminsky, R., and Maser, P. (2007) Phenotyping and genotyping of Haemonchus contortus isolates reveal a new putative candidate mutation for benzimidazole resistance in nematodes. Veterinary Parasitology 144, 313320. doi: 10.1016/j.vetpar.2006.10.003.CrossRefGoogle ScholarPubMed
IBGE (2021a) Censo Agropecuário 2021. Retrieved from Instituto Brasileiro de Geografia e Estatistica website: https://sidra.ibge.gov.br/tabela/3939Google Scholar
IBGE (2021b) Cidades. Retrieved from Instituto Brasileiro de Geografia e Estatistica website: https://cidades.ibge.gov.br/ (accessed April 2023).Google Scholar
Kotze, A.C., Cowling, K., Bagnall, N.H., Hines, B.M., Ruffell, A.P., Hunt, P.W., and Coleman, G.T. (2012) Relative level of thiabendazole resistance associated with the E198A and F200Y SNPs in larvae of a multi-drug resistant isolate of Haemonchus contortus. International Journal for Parasitology: Drugs and Drug Resistance 2, 9297. doi: 10.1016/j.ijpddr.2012.02.003.Google ScholarPubMed
Kwa, M.S., Veenstra, J.G., and Roos, M.H. (1994) Benzimidazole resistance in Haemonchus contortus is correlated with a conserved mutation at amino acid 200 in β-tubulin isotype 1. Molecular and Biochemical Parasitology 63(2), 299303. doi: 10.1016/0166-6851(94)90066-3CrossRefGoogle ScholarPubMed
Lacey, E. (1988) The role of the cytoskeletal protein, tubulin, in the mode of action and mechanism of drug resistance to benzimidazoles. International Journal for Parasitology: Drugs and Drug Resistance 18, 885936. doi: 10.1016/0020-7519(88)90175-0.CrossRefGoogle ScholarPubMed
Lambert, S.M., Nishi, S.M., Mendonça, L.R., Souza, B.M.P.S., Julião, F.S., Gusmão, P.S., and Almeida, M.A.O. (2017) Genotypic profile of benzimidazole resistance associated with SNP F167Y and F200Y beta-tubulin gene in Brazilian populations of Haemonchus contortus of goats. Veterinary Parasitology: Regional Studies and Reports 8, 2834. doi: 10.1016/j.vprsr.2017.01.006.Google ScholarPubMed
Leathwick, D.M., Hosking, B.C., Bisset, S.A., and Mckay, C.H. (2009) Managing anthelmintic resistance: is it feasible in New Zealand to delay the emergence of resistance to a new anthelmintic class? New Zealand Veterinary Journal 57, 181192. doi: 10.1080/00480169.2009.36900.CrossRefGoogle ScholarPubMed
Molento, M.B., van Wyk, J.A., and Coles, G.C. Sustainable worm management. The Veterinary Record 155(3), 9596. 2004. PMID: 15311810Google ScholarPubMed
Mottier, M.D.L. and Prichard, R.K. (2008) Genetic analysis of a relationship between macrocyclic lactone and benzimidazole anthelmintic selection on Haemonchus contortus. Pharmacogenetics and Genomics 18, 129140. doi: 10.1097/FPC.0b013e3282f4711d.CrossRefGoogle Scholar
Niciura, S.C.M., Veríssimo, C.J., Gromboni, J.G.G., Rocha, M.I.P., Mello, S.S., Barbosa, C.M.P., Chiebao, D.P., Cardoso, D., Silva, G.S., Otsuk, I.P., Pereira, J.R., Ambrosio, L.A., Nardon, R.F., Ueno, T.E.H., and Molento, M.B. (2012) F200Y polymorphism in the β-tubulin gene in field isolates of Haemonchus contortus and risk factors of sheep flock management practices related to anthelmintic resistance. Veterinary Parasitology 190, 608612. doi: 10.1016/j.vetpar.2012.07.016.CrossRefGoogle ScholarPubMed
Nunes, R.L., Santos, L.L., Bastianetto, E., Oliveira, D.A.A., and Brasil, B.S.A.F. (2013) Frequency of benzimidazole resistance in Haemonchus contortus populations isolated from buffalo, goat and sheep herds. Revista Brasileira de Parasitologia Veterinária 22, 548553. Doi: 10.1590/S1984-29612013000400015.CrossRefGoogle ScholarPubMed
Odoi, A., Gathuma, J.M., Gachuiri, C.K., and Omore, A. (2007) Risk factors of gastrointestinal nematode parasite infections in small ruminants kept in smallholder mixed farms in Kenya. BMC Veterinary Research 3, 111. doi: 10.1186/1746-6148-3-6.CrossRefGoogle ScholarPubMed
Parvin, S., Dey, A.R., Shohana, N.N., Anisuzzaman, Talukder Md.H and Alam, M.Z. (2024) Haemonchus contortus, an obligatory haematophagus worm infection in small ruminants: Population genetics and genetic diversity. Saudi Journal of Biological Sciences 31, 104030. https://doi.org/10.1016/j.sjbs.2024.104030CrossRefGoogle ScholarPubMed
Reynecke, D.P., van Wyk, J.A., Gummow, B., Dorny, P., and Boomker, J. (2011) A stochastic model accommodating the FAMACHA system for estimating worm burdens and associated risk factors in sheep naturally infected with Haemonchus contortus. Veterinary Parasitology 177, 231241. doi: 10.1016/j.vetpar.2011.01.042.CrossRefGoogle ScholarPubMed
Roberts, F.H.S. and O’Sullivan, P.J. (1950) Methods for egg counts and larval cultures for strongyles infesting the gastrointestinal tract of cattle. Australian Journal of Agricultural Research 1, 99102. doi: 10.1071/AR9500099.CrossRefGoogle Scholar
Santos, J.M.L., Monteiro, J.P., Ribeiro, W.L.C., Macedo, I.T.F., Camurca-Vasconcelos, A.L., Vieira, L.S., and Bevilaqua, C.M.L. (2014) Identification and quantification of benzimidazole resistance polymorphisms in Haemonchus contortus isolated in Northeastern Brazil. Veterinary Parasitology 199, 160164. doi: 10.1016/j.vetpar.2013.11.006.CrossRefGoogle ScholarPubMed
Santos, J.M.L., Vasconcelos, J.F., Frota, G.A., Ribeiro, W.L.C., André, W.P.P., Vieria, L.S., Teixeira, M., Bevilaqua, M.L., and Monteiro, J.P. (2017a) Haemonchus contortus β-tubulin isotype 1 gene F200Y and F167Y SNPs are both selected by ivermectin and oxfendazole treatments with differing impacts on anthelmintic resistance. Veterinary Parasitology 248, 9095. doi: 10.1016/j.vetpar.2017.11.003.CrossRefGoogle ScholarPubMed
Santos, J.M.L., Monteiro, J.P., Ribeiro, W.L.C., Macedo, I.T.F., Araújo-Filho, J.V., Andre, W.P.P., Araújo, P.R.M., Vasconcelos, J.F., Freitas, E.P., Camurça-Vasconcelos, A.L.F., Vieira, L.S., and Bevilaqua, C.M.L. (2017b) High levels of benzimidazole resistance and β-tubulin isotype 1 SNP F167Y in Haemonchus contortus populations from Ceará State Brazil. Small Ruminant Research 146, 4852. Doi: 10.1016/j.smallrumres.2016.11.023.CrossRefGoogle Scholar
Silvestre, A. and Cabaret, J. (2002) Mutation in position 167 of isotype 1 β-tubulin gene of Trichostrongylid nematodes: role in benzimidazole resistance?. Molecular and Biochemical Parasitology 120(2), 297300. doi: 10.1016/s0166-6851(01)00455-8CrossRefGoogle ScholarPubMed
Thomaz-Soccol, V., Sotomaior, C., Souza, F.P., Castro, E.A., Pessôa Silva, M.C., and Milezewski, V. (1996) Occurrence of resistance to anthelmintics in sheep in Paraná State, Brazil. Veterinary Record 139, 421422. doi: 10.1136/vr.139.17.421.CrossRefGoogle Scholar
Thomaz-Soccol, V., Souza, F.P., Sotomaior, C., Castro, E.A., Milczewski, V., Mocelin, G., and Pessôa Silva, M.C. (2004) Resistance of gastrointestinal nematodes to anthelmintics in sheep (Ovis aries). Brazilian Archives of Biology and Technology 47, 4147. doi: 10.1590/S1516-89132004000100006.CrossRefGoogle Scholar
Van Wyk, J.A. and Mayhew, E. (2013) Morphological identification of parasitic nematode infective larvae of small ruminants and cattle: a practical lab guide. Onderstepoort Journal of Veterinary Research 80, 114. doi: 10.4102/ojvr.v80i1.539.CrossRefGoogle ScholarPubMed
Veríssimo, C.J., Niciura, S.C.M., Alberti, A.L.L., Rodrigues, C.F.C., Barbosa, C.M.P., Chiebao, D.P., Cardoso, D., Silva, G.S., Pereira, J.R., Margatho, L.F.F., Costa, R.L.D., Nardon, R.F., Ueno, T.E.H., Curci, V.C.L.M., and Molento, M.B. (2012) Multidrug and multispecies resistance in sheep flocks from São Paulo state, Brazil. Veterinary Parasitology 187, 209216. doi: 10.1016/j.vetpar.2012.01.013.CrossRefGoogle ScholarPubMed
Vilela, V.L.R., Bezerra, H.M.F.F., Bezerra, R.A., Dantas, M.O., Alcântara, E.T., Oliveira, L.V.S., Nóbrega, K.S., Calazans, F.B., Feitosa, T.F., Braga, F.R., and Molento, M.B. Sustainable agriculture: the use of FAMACHA method in Santa Ines sheep in the Semi-arid region of Brazil. Semina 42(3), S1, 16471662, 2021. DOI: 10.5433/1679-0359.2021v42n3Supl1p1647Google Scholar
Workentine, M. L., Chen, R., Zhu, S., Gavriliuc, S., Shaw, N., Rijke, J. D., Redman, E. M., Avramenko, R. W., Wit, J., Poissant, J., and Gilleard, J. S. (2020). A database for ITS2 sequences from nematodes. BMC Genetics 21(1), 14. doi: 10.1186/s12863-020-00880-0.CrossRefGoogle ScholarPubMed
Zajac, A.M. and Garza, J. (2020) Biology, epidemiology, and control of gastrointestinal nematodes of small ruminants. Veterinary Clinics of North America: Food Animal Practice 36(1), 7387. doi: 10.1016/j.cvfa.2019.12005.Google ScholarPubMed
Zhang, Z., Gasser, R.B., Yang, X., Yin, F., Zhao, G., Bao, M., Pan, B., Huang, W., Wang, C., Zou, F., Zhou, Y., Zhao, J., Fang, R., and Hu, M. (2016) Two benzimidazole resistance-associated SNPs in the isotype-1 β-tubulin gene predominate in Haemonchus contortus populations from eight regions in China. International Journal for Parasitology. Drugs Drug Resistance 6(3), 199206. doi: 10.1016/j.ijpddr.2016.10.001.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Location of the municipalities with sheep and goat farms in Western Paraná, Brazil. Obs. Corresponding municipalities: 1. Foz do Iguaçu, 2. Matelândia, 3. Medianeira, 4. Ramilândia, 5. Santa Terezinha de Itaipu, 6. São Miguel do Iguaçu, 7. Serranópolis do Iguaçu, 8. Vera Cruz do Oeste.

Figure 1

Table 1. Municipalities included in the study, total area (km2), population size, main economic activity, and the number of sheep (S) and goats (G) of Western Parana, Brazil

Figure 2

Table 2. Frequency (%) of health measures adopted by 18 sheep and goat farmers in Western Paraná, Brazil

Figure 3

Table 3. Average qPCR and cycle quantification (Cq) and allelic frequencies of the F200Y, F167Y, and E198A mutations from Haemonchus contortus of sheep and goat farms in Western Paraná, Brazil

Figure 4

Table 4. Allelic frequencies (%) of the F200Y, F167Y, and E198A mutations of Haemonchus contortus in sheep and goat farms (isolates) with and without anthelmintic rotation (yes/no), in Paraná, Brazil