The excessive and uncontrolled use of anthelmintics, e.g. ivermectin (IVM) for the treatment of livestock parasites has led to widespread resistance in gastrointestinal nematodes, such as Haemonchus contortus. There is an urgent need for better management of drug-use in nematode control and development of novel anthelmintics. Discovery and identification of anthelmintic resistance-associate molecules/markers can provide a basis for rational anthelmintics-use and development of novel drugs. Recent studies have shown that ivermectin resistance in H. contortus is likely to be multi-genic in nature except for several genes coding for IVM target and efflux pump. However, no other IVM resistance-associated genes were characterized by conventional methods or strategies. In the present study we adopted a new strategy, i.e. using genome-wide single nucleotide polymorphism (SNP) analysis based on 2b-RAD sequencing, for discovering SNPs markers across the genomes in both IVM susceptible and resistant isolates of H. contortus and identifying potential IVM resistance-associated genes. We discovered 2962 and 2667 SNPs within both susceptible and resistant strains of H. contortus, respectively. A relative lower and similar genetic variations were observed within both resistant and susceptible strains (average π values were equal to 0.1883 and 0.1953, respectively); whereas a high genetic variation was found across both strains (average π value was equal to 0.3899). A significant differentiation across 2b-RAD tags nucleotide sites was also observed between the two strains (average FST value was equal to 0.3076); the larger differences in average FST were observed at SNPs loci between coding and noncoding (including intronic) regions. Comparison between resistant and susceptible strains revealed that 208 SNPs loci exhibited significantly elevated FST values, 24 SNPs of those loci were located in the CDS regions of the nine genes and were likely to have signature of IVM directional selection. Seven of the nine candidate genes were predicted to code for some functional proteins such as potential IVM target and/or efflux pump proteins, component proteins of receptor complex in membrane on neuromuscular cells, and transcriptional regulation proteins. Those genes might be involved in resistance to IVM. Our data suggest that candidate genes putatively associated with resistance to IVM in H. contortus may be identified by genome-wide SNP analysis using 2b-RAD sequencing.

A major problem faced by the agricultural industry is the resistance of Haemonchus contortus to anthelmintic drugs. For a better understanding of the response of H. contortus to IVM and for the screening of drug-resistance-related genes, we used RNA sequencing and isobaric tags for relative and absolute quantification (iTRAQ) technology to detect the transcriptomic and proteomic changes in H. contortus after ivermectin treatment. An integrated analysis of the two omics showed that the differentially expressed genes and proteins were significantly enriched in the pathways of amino acid degradation, the metabolism of xenobiotics by cytochrome P450, the biosynthesis of amino acids, and the tricarboxylic acid cycle. We found that the upregulated UDP-glycosyltransferases (UGT), glutathione S-transferase (GST), cytochrome P450 (CYP), and p-glycoprotein (Pgp) genes play important roles in drug resistance in H. contortus. Our work will help in the understanding of the transcriptome and proteome changes in H. contortus after IVM and will facilitate the discovery of genes related to drug resistance. This information can be further applied to increase the understanding of the response of IVM in relation to H. contortus.

Nematodirus causes major economic losses in the development of the livestock industry, as they are common gastrointestinal parasites of cattle, sheep, and other ruminants. The present study investigated Nematodirus infections in sheep in the Hulunbuir and Xilingol Grasslands, two well-known natural pastures in Inner Mongolia, China. The results indicated that the average infection rate was more than 50%. Furthermore, a naturally occurring strain of Nematodirus was obtained using numerous isolation methods in the survey region. Conventional morphology and molecular biology were used to identify the strain. The larval migration inhibition test was used to determine the ivermectin level in the isolated strain. The results indicated that the larvae belonged to the species Nematodirus oiratianus. The strain was found to be ivermectin-resistant. Thus, these results recorded the detailed development processes of Nematodirus eggs, providing novel insights regarding the prevention and control of nematodes by using conventional anthelmintic regimens and by offering valuable species resources for studying anthelmintic resistance in nematodes.

Haemonchus contortus poses a severe hazard to the healthy development of the sheep industry and threatens the welfare of sheep. Ivermectin is the primary anthelmintic used for the prevention and treatment of H. contortus parasitism. However, the widespread and uncontrolled application of ivermectin has resulted in the development and spread of resistant strains of H. contortus. P-glycoprotein (P-gp) plays important roles in the pharmacology and toxicology of ivermectin, and changes in P-gp expression levels can be used to analyze the resistance of H. contortus to ivermectin. This study aimed to analyze the effects of ivermectin on P-gp expression in H. contortus L3 larvae isolated from China and to evaluate whether changes in P-gp expression levels can be used to analyze resistant H. contortus strains. In the absence of drug treatment, the ivermectin-resistant strains isolated in China showed increased expression of P-gp11 (p < 0.01) compared with sensitive strains from elsewhere, whereas the expressions of P-gp2 and P-gp9.1 were downregulated (p < 0.01). When the same strain was compared before and after drug treatment, obvious differences in expression were observed between the different strains. Ivermectin-induced P-gp expression was found to be very complex among the L3 larvae of different strains. In addition, it was confirmed that using P-gp to determine ivermectin resistance in H. contortus strains from different geographic environments can yield different results.