Transcriptomic analyses implicate neuronal plasticity and chloride homeostasis in ivermectin resistance and response to treatment in a parasitic nematode

Efficacy of fenbendazole against gastrointestinal nematodes in naturally infected goats in Maputo Province, Mozambique using in vivo, in vitro and molecular assessment

Anthelmintic resistance occurs worldwide in strongyles of ruminants but data from low-income countries are sparse and rarely apply most up to date methods, while effects of management practices in these countries are poorly documented. In Mozambique, benzimidazole resistance has been previously reported; the present study followed this up in detail, applying in vivo faecal egg count (FEC) reduction test (FECRT), in vitro egg hatch test (EHT) and molecular deep amplicon sequencing approaches targeting the internal transcribed spacer 2 (ITS-2, nemabiome) and the isotype 1 β-tubulin gene to determine the resistance status on farms and the strongyle species involved. Adult Landim goats (433) from six semi-intensive and ten extensive farms (22-30 animals/farm) from Maputo Province were visited April 2021 to February 2022. Fenbendazole (5 mg/kg bw, Panacur®) was administered orally and FEC determined using Mini-FLOTAC. The eggCounts package was used to calculate FECRs with 90% confidence intervals from paired day 0 and 14 data. In vivo and in vitro tests detected AR on 5/16 (31%) farms. This included 1/10 extensive and 4/6 semi-intensive farms. The odds of finding resistant strongyles on a semi-intensive commercial farm was 40-fold higher than on an extensive farm (p = 0.016, logistic regression). A strong, negative correlation was observed between FECRT and EHT EC50 values (Pearson's R = -0.83, P = 0.001; Cohen's κ coefficient 1.0). Nemabiome data showed that Haemonchus contortus, Trichostrongylus colubriformis and unclassified Oesophagostomum closely related to Oesophagostomum columbianum were most abundant before treatment and in particular H. contortus frequencies increased after treatment. Benzimidazole resistance associated polymorphisms were detected in H. contortus and T. colubriformis. Moreover, there were hints that resistance alleles were present in Trichostrongylus axei and Teladorsagia circumcincta. Farmers should regularly test the efficacy of anthelmintics used and consider more sustainable worm control approaches to reduce selection for resistance.

Chromosomal genome assembly resolves drug resistance loci in the parasitic nematode Teladorsagia circumcincta

The parasitic nematode Teladorsagia circumcincta is one of the most important pathogens of sheep and goats in temperate climates worldwide and can rapidly evolve resistance to drugs used to control it. To understand the genetics of drug resistance, we have generated a highly contiguous genome assembly for the UK T. circumcincta isolate, MTci2. Assembly using PacBio long-reads and Hi-C long-molecule scaffolding together with manual curation resulted in a 573 Mb assembly (N50 = 84 Mb, total scaffolds = 1,286) with five autosomal and one sex-linked chromosomal-scale scaffolds consistent with its karyotype. The genome resource was further improved via annotation of 22,948 genes, with manual curation of over 3,200 of these, resulting in a robust and near complete resource (96.3% complete protein BUSCOs) to support basic and applied research on this important veterinary pathogen. Genome-wide analyses of drug resistance, combining evidence from three distinct experiments, identified selection around known candidate genes for benzimidazole, levamisole and ivermectin resistance, as well as novel regions associated with ivermectin and moxidectin resistance. These insights into contemporary and historic genetic selection further emphasise the importance of contiguous genome assemblies in interpreting genome-wide genetic variation associated with drug resistance and identifying key loci to prioritise in developing diagnostic markers of anthelmintic resistance to support parasite control.

A systematic review of the molecular mechanisms related to anthelmintic resistance in Haemonchus contortus: A contemporary narrative

Haemonchus contortus is a gastrointestinal parasite that affects ruminants (cattle, sheep, etc.), having a significant welfare impact worldwide. The rise of anthelmintic resistance poses a growing challenge to adequate control, compromising the success of treatments. This study presents a systematic review of the molecular mechanisms involved in the resistance of H. contortus to anthelmintic drugs. Following an extensive literature search (9075 total articles/excluding duplications), 61 articles were examined. From these, benzimidazoles (BZD) and macrocyclic lactone (ML) were the most reported drug classes (17 and 29, respectively). The mutations in the β-tubulin gene were the primary mechanism of BZD resistance. Important comparisons from early reports of resistance mechanisms to ML (published before 2020) mainly based on parasite-population gene expression (e.g., ligand-gated chloride channels, LGCC, and P-glycoproteins, PGP) with more recent genomic and transcriptomic data (e.g., transcription factor, cky-1 gene) are highlighted. Additionally, resistance mechanisms to levamisole (LEV) and monepantel are discussed, showing evidence of polymorphisms in genes related to the nicotinic acetylcholine receptors (nAChR). Considering the available data, it is possible to divide the reports into two technological periods, evidencing that the availability of a chromosome-level genome for H. contortus in association with experiments based on controlled genetic crosses and transcriptome-wide data allowed for the visualization of genes and polymorphisms that were previously indistinguishable from unrelated genetic variation (i.e., genetic noise). Therefore, the study of anthelmintic resistance mechanisms is facing new possibilities, reflecting the large data banks and the speed at which this information is being processed. We suggest that new publications on drug resistance should adopt the approaches and refer to this new era of scientific discoveries. Consistent data interpretation, including artificial intelligence (AI) support, will help us to suggest novel biological mechanisms involved in drug resistance and predict its evolution, allowing a more comprehensible approach toward sustainable parasite control strategies.

Differences in constitutive gene expression of cytochrome P450 enzymes and ATP-binding cassette transporter gene expression between a susceptible and a highly macrocyclic lactone-resistant Haemonchus contortus isolate in the absence of drug-inducible expression

BACKGROUND

Anthelmintic resistance in ruminants is a widespread problem that has a severe impact on productivity and animal welfare. The helminth Haemonchus contortus is generally considered the most important parasite in small ruminants due to its high pathogenicity and the widespread occurrence of anthelmintic resistance in it. Although the molecular mechanisms associated with resistance against the anthelmintics benzimidazoles (BZs) and levamisole are relatively well understood, the resistance mechanisms against the widely used anthelmintic macrocyclic lactones (MLs) ivermectin (IVM) and moxidectin (MOX) remain poorly understood. Detoxifying enzymes and xenobiotic transporters have been frequently proposed to play a role in ML resistance in multiple organisms, including nematodes.

METHODS

The reference genome of H. contortus was screened for cytochrome P450 genes (cyp genes) by using the Basic Local Alignment Search Tool, and maximum-likelihood phylogenetic analysis was used to assign the sequences to gene families. Fourth-stage larvae of the susceptible (McMaster) and the ML-resistant (Berlin-selected) H. contortus isolates were generated in vitro and compared regarding basal expression levels of cyp genes and ATP-binding cassette (ABC) transporters by using RNA sequencing. The resistant isolate was further incubated with 100 nM IVM or MOX for 3, 6 and 12 h, and the effects of incubation time and drugs were evaluated.

RESULTS

Twenty-five cyp genes were identified in the H. contortus genome and assigned to 13 different families. The ML-resistant isolate showed significantly higher and lower constitutive expression of 13 and four cyp genes, respectively. Out of the 50 ABC transporter genes, only six showed significantly higher expression in the ML-resistant isolate, while 12 showed lower expression. The fold changes were in general low (range 0.44-5.16). Only pgp-13 showed significant downregulation in response to IVM (0.77 fold change at 6 h, 0.96 fold change at 12 h) and MOX (0.84 fold change at 12 h). In contrast, mrp-5 was significantly, albeit minimally, upregulated in the presence of IVM, but not MOX, after 12 h (1.02 fold change).

CONCLUSIONS

Despite little observable ML-inducible gene expression in the isolate examined here, some of the changes in the baseline expression levels might well contribute to ML resistance in the context of additional changes in a multigenic resistance model. However, neither cyp genes nor the ABC transporters appear to be the main drivers that can explain the high levels of resistance observed in the resistant isolate examined here.

P-glycoproteins in anthelmintic safety, efficacy, and resistance

P-glycoprotein (PGP) is a pivotal transmembrane transporter governing the cellular flux of diverse substances shielding mammals from toxics. It can thwart the effectiveness of medicines such as ivermectin (IVM) and other macrocyclic lactone (ML) anthelmintics, undermining therapeutic efforts. We analyze the role of PGPs in limiting the toxicity of these drugs in hosts, and their potential contribution to anthelmintic resistance in nematodes. Targeting nematode PGPs to increase drug sensitivity to MLs seems interesting, but is hampered by the lack of selective inhibitors. The nuclear hormone receptor (NHR)-8 should be seriously considered as a target because it upregulates multiple PGPs involved in anthelmintic resistance and it is specific to nematodes. This would advance our understanding of host-pathogen dynamics and foster innovative therapeutic strategies.

Genome-wide RNA interference of the nhr gene family in barber's pole worm identified members crucial for larval viability in vitro

Nuclear hormone receptors (NHRs) are emerging target candidates against nematode infection and resistance. However, there is a lack of comprehensive information on NHR-coding genes in parasitic nematodes. In this study, we curated the nhr gene family for 60 major parasitic nematodes from humans and animals. Compared with the free-living model organism Caenorhabditis elegans, a remarkable contraction of the nhr family was revealed in parasitic species, with genetic diversification and conservation unveiled among nematode Clades I (10-13), III (16-42), IV (33-35) and V (25-64). Using an in vitro biosystem, we demonstrated that 40 nhr genes in a blood-feeding nematode Haemonchus contortus (clade V; barber's pole worm) were responsive to host serum and one nhr gene (i.e., nhr-64) was consistently stimulated by anthelmintics (i.e., ivermectin, thiabendazole and levamisole); Using a high-throughput RNA interference platform, we knocked down 43 nhr genes of H. contortus and identified at least two genes that are required for the viability (i.e., nhr-105) and development (i.e., nhr-17) of the infective larvae of this parasitic nematode in vitro. Harnessing this preliminary functional atlas of nhr genes for H. contortus will prime the biological studies of this gene family in nematode genetics, infection, and anthelmintic metabolism within host animals, as well as the promising discovery of novel intervention targets.

Comparison of P-glycoprotein gene expression of two Haemonchus contortus isolates from Yucatan, Mexico, with resistant or susceptible phenotype to ivermectin in relation to a susceptible reference strain

The variability in the expression of different P-glycoprotein (P-gp) genes in parasitic nematodes of ruminants such as Haemonchus contortus (Hco-pgp) may be caused by different factors including nematode biology, geographical region and anthelmintic pressure. This study analysed the relative expression level of 10 P-gp genes in two H. contortus (Hco-pgp) field isolates from Yucatan, Mexico: 1) PARAISO (IVM-resistant) and 2) FMVZ-UADY (IVM-susceptible). These isolates were compared with a susceptible reference isolate from Puebla, Mexico, namely "CENID-SAI". In all cases H. contortus adult males were used. The Hco-pgp genes (1, 2, 3, 4, 9, 10, 11, 12, 14 and 16) were analysed for each isolate using the RT-qPCR technique. The Hco-pgp expressions were pairwise compared using the 2-ΔΔCt method and a t-test. The PARAISO isolate showed upregulation compared to the CENID-SAI isolate for Hco-pgp 1, 3, 9, 10 and 16 (P < 0.05), and the PARAISO isolate showed upregulation vs. FMVZ-UADY isolate for Hco-pgp 2 and 9 (P < 0.05), displaying 6.58- and 5.93-fold differences (P < 0.05), respectively. In contrast, similar Hco-pgp gene expression levels were recorded for FMVZ-UADY and CENID-SAI isolates except for Hco-pgp1 (P <0.1), which presented a significant upregulation (6.08-fold). The relative expression of Hco-pgp allowed confirming the IVM-resistant status of the PARAISO isolate and the IVM-susceptible status of the FMVZ-UADY isolate when compared to the CENID-SAI reference isolate. Therefore, understanding the association between the Hco-pgp genes expression of H. contortus and its IVM resistance status could help identifying the genes that could be used as molecular markers in the diagnosis of IVM resistance. However, it is important to consider the geographic origin of the nematode isolate and the deworming history at the farm of origin.

The response of Anisakis simplex (s. s.) to anthelmintics - Specific changes in xenobiotic metabolic processes

Anisakiasis is a parasitic disease transmitted through the consumption of raw or undercooked fish and cephalopods that are infected with larvae of Anisakis simplex (sensu stricto) or Anisakis pegreffii. The purpose of this study was to investigate how A. simplex (s. s.) responds to the influence of anthelmintics such as ivermectin (IVM) and pyrantel (PYR). In vitro experiments were conducted using larvae at two developmental stages of A. simplex (s. s.) (L3 and L4) obtained from Baltic herring (Clupea harengus membras). Larvae were cultured with different concentrations of IVM or PYR (1.56, 3.125, and 6.25 μg/mL) for various durations (3, 6, 9, and 12 h) under anaerobic conditions (37 °C, 5% CO2). The gene expression of actin, ABC transporter, antioxidant enzymes, γ-aminobutyric acid receptors, and nicotinic acetylcholine receptors, as well as the oxidative status were analyzed. The results showed that A. simplex (s. s.) L3 stage had lower mobility when cultured with PYR compared to IVM. The analysis of relative gene expression revealed significant differences in the mRNA level of ABC transporters after treatment with IVM and PYR, compared to the control group. Similar patterns were observed in the gene expression of antioxidant enzymes in response to both drugs. Furthermore, the total antioxidant capacity (TAC) and glutathione S-transferase (GST) activity were higher in the treatment groups than in the control group. These findings suggest a relationship between the expression of the studied genes, including those related to oxidative metabolism, and the effectiveness of the tested drugs.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans

Albendazole (a benzimidazole) and ivermectin (a macrocyclic lactone) are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channels (GluCls), but it is unknown whether GluCl genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of-function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required the loss of two GluCl genes (avr-14 and avr-15). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole conditions and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance.

Generation of IgG antibodies against Strongyloides stercoralis in mice via immunization with recombinant antigens A133 and Ss-IR

Strongyloidiasis, caused by the nematode Strongyloides stercoralis, remains a threat to global public health, and a vaccine would be useful to control the disease, especially in developing countries. This study aimed to evaluate the efficacy of recombinant proteins, A133 and Ss-IR, as potential vaccine candidates against strongyloidiasis by investigating the humoral and cellular immune responses in immunized mice. Respective antigens were adjuvanted with Complete Freund's Adjuvant (prime) and Incomplete Freund's Adjuvant (boost) and administered intraperitoneally (prime) and subcutaneously (boost) to female BALB/c mice. For antigen-only doses, only antigens were injected without adjuvants. Altogether, 1 prime dose, 4 booster doses, and 2 antigen-only doses were administered successively. ELISAs were conducted to assess the antibody responses, along with flow cytometry and cytokine ELISA to elucidate the cellular immune responses. Results showed that A133 and Ss-IR induced the production of IgG1 and IgG2a, with A133 generating more robust IgG2a responses than Ss-IR. Flow cytometry findings indicated that effector CD8+T-cells and memory B-cells activity were upregulated significantly for A133 only, whereas cytokine ELISA demonstrated that a Th1/Th2/Th17 mixed cell responses were triggered upon vaccination with either antigen. This preliminary study illustrated the good potential of recombinant A133 and Ss-IR as vaccine candidates against S. stercoralis. It provided information on the probable immune mechanism involved in host defence and the elicitation of protection against S. stercoralis.

Altered larval activation response associated with multidrug resistance in the canine hookworm Ancylostoma caninum

Parasitic gastrointestinal nematodes pose significant health risks to humans, livestock, and companion animals, and their control relies heavily on the use of anthelmintic drugs. Overuse of these drugs has led to the emergence of resistant nematode populations. Herein, a naturally occurring isolate (referred to as BCR) of the dog hookworm, Ancylostoma caninum, that is resistant to 3 major classes of anthelmintics is characterized. Various drug assays were used to determine the resistance of BCR to thiabendazole, ivermectin, moxidectin and pyrantel pamoate. When compared to a drug-susceptible isolate of A. caninum, BCR was shown to be significantly resistant to all 4 of the drugs tested. Multiple single nucleotide polymorphisms have been shown to impart benzimidazole resistance, including the F167Y mutation in the β-tubulin isotype 1 gene, which was confirmed to be present in BCR through molecular analysis. The frequency of the resistant allele in BCR was 76.3% following its first passage in the lab, which represented an increase from approximately 50% in the founding hookworm population. A second, recently described mutation in codon 134 (Q134H) was also detected at lower frequency in the BCR population. Additionally, BCR exhibits an altered larval activation phenotype compared to the susceptible isolate, suggesting differences in the signalling pathways involved in the activation process which may be associated with resistance. Further characterization of this isolate will provide insights into the mechanisms of resistance to macrocyclic lactones and tetrahydropyrimidine anthelmintics.

Getting around the roundworms: Identifying knowledge gaps and research priorities for the ascarids

The ascarids are a large group of parasitic nematodes that infect a wide range of animal species. In humans, they cause neglected diseases of poverty; many animal parasites also cause zoonotic infections in people. Control measures include hygiene and anthelmintic treatments, but they are not always appropriate or effective and this creates a continuing need to search for better ways to reduce the human, welfare and economic costs of these infections. To this end, Le Studium Institute of Advanced Studies organized a two-day conference to identify major gaps in our understanding of ascarid parasites with a view to setting research priorities that would allow for improved control. The participants identified several key areas for future focus, comprising of advances in genomic analysis and the use of model organisms, especially Caenorhabditis elegans, a more thorough appreciation of the complexity of host-parasite (and parasite-parasite) communications, a search for novel anthelmintic drugs and the development of effective vaccines. The participants agreed to try and maintain informal links in the future that could form the basis for collaborative projects, and to co-operate to organize future meetings and workshops to promote ascarid research.

Gene co-expression network analysis reveal core responsive genes in Parascaris univalens tissues following ivermectin exposure

Anthelmintic resistance in equine parasite Parascaris univalens, compromises ivermectin (IVM) effectiveness and necessitates an in-depth understanding of its resistance mechanisms. Most research, primarily focused on holistic gene expression analyses, may overlook vital tissue-specific responses and often limit the scope of novel genes. This study leveraged gene co-expression network analysis to elucidate tissue-specific transcriptional responses and to identify core genes implicated in the IVM response in P. univalens. Adult worms (n = 28) were exposed to 10-11 M and 10-9 M IVM in vitro for 24 hours. RNA-sequencing examined transcriptional changes in the anterior end and intestine. Differential expression analysis revealed pronounced tissue differences, with the intestine exhibiting substantially more IVM-induced transcriptional activity. Gene co-expression network analysis identified seven modules significantly associated with the response to IVM. Within these, 219 core genes were detected, largely expressed in the intestinal tissue and spanning diverse biological processes with unspecific patterns. After 10-11 M IVM, intestinal tissue core genes showed transcriptional suppression, cell cycle inhibition, and ribosomal alterations. Interestingly, genes PgR028_g047 (sorb-1), PgB01_g200 (gmap-1) and PgR046_g017 (col-37 & col-102) switched from downregulation at 10-11 M to upregulation at 10-9 M IVM. The 10-9 M concentration induced expression of cuticle and membrane integrity core genes in the intestinal tissue. No clear core gene patterns were visible in the anterior end after 10-11 M IVM. However, after 10-9 M IVM, the anterior end mostly displayed downregulation, indicating disrupted transcriptional regulation. One interesting finding was the non-modular calcium-signaling gene, PgR047_g066 (gegf-1), which uniquely connected 71 genes across four modules. These genes were enriched for transmembrane signaling activity, suggesting that PgR047_g066 (gegf-1) could have a key signaling role. By unveiling tissue-specific expression patterns and highlighting biological processes through unbiased core gene detection, this study reveals intricate IVM responses in P. univalens. These findings suggest alternative drug uptake of IVM and can guide functional validations to further IVM resistance mechanism understanding.

Genome structure and population genomics of the canine heartworm Dirofilaria immitis

The heartworm, Dirofilaria immitis, is a filarial parasitic nematode responsible for significant morbidity and mortality in wild and domesticated canids. Resistance to macrocyclic lactone drug prevention represents a significant threat to parasite control and has prompted investigations to understand the genetic determinants of resistance. This study aimed to improve the genomic resources of D. immitis to enable a more precise understanding of how genetic variation is distributed within and between parasite populations worldwide, which will inform the likelihood and rate by which parasites, and in turn, resistant alleles, might spread. We have guided the scaffolding of a recently published genome assembly for D. immitis (ICBAS_JMDir_1.0) using the chromosomal-scale reference genomes of Brugia malayi and Onchocerca volvulus, resulting in an 89.5 Mb assembly composed of four autosomal- and one sex-linked chromosomal-scale scaffolds representing 99.7% of the genome. Publicly available and new whole-genome sequencing data from 32 D. immitis samples from Australia, Italy and the USA were assessed using principal component analysis, nucleotide diversity (Pi) and absolute genetic divergence (Dxy) to characterise the global genetic structure and measure within- and between-population diversity. These population genetic analyses revealed broad-scale genetic structure among globally diverse samples and differences in genetic diversity between populations; however, fine-scale subpopulation analysis was limited and biased by differences between sample types. Finally, we mapped single nucleotide polymorphisms previously associated with macrocyclic lactone resistance in the new genome assembly, revealing the physical linkage of high-priority variants on chromosome 3, and determined their frequency in the studied populations. This new chromosomal assembly for D. immitis now allows for a more precise investigation of selection on genome-wide genetic variation and will enhance our understanding of parasite transmission and the spread of genetic variants responsible for resistance to treatment.

Understanding anthelmintic resistance in livestock using "omics" approaches

Widespread and improper use of various anthelmintics, genetic, and epidemiological factors has resulted in anthelmintic-resistant (AR) helminth populations in livestock. This is currently quite common globally in different livestock animals including sheep, goats, and cattle to gastrointestinal nematode (GIN) infections. Therefore, the mechanisms underlying AR in parasitic worm species have been the subject of ample research to tackle this challenge. Current and emerging technologies in the disciplines of genomics, transcriptomics, metabolomics, and proteomics in livestock species have advanced the understanding of the intricate molecular AR mechanisms in many major parasites. The technologies have improved the identification of possible biomarkers of resistant parasites, the ability to find actual causative genes, regulatory networks, and pathways of parasites governing the AR development including the dynamics of helminth infection and host-parasite infections. In this review, various "omics"-driven technologies including genome scan, candidate gene, quantitative trait loci, transcriptomic, proteomic, and metabolomic approaches have been described to understand AR of parasites of veterinary importance. Also, challenges and future prospects of these "omics" approaches are also discussed.

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of parasitic nematodes.

Recent Advances in the Control of Endoparasites in Ruminants from a Sustainable Perspective

Consumer awareness of animal welfare and environmental health has led to a plateau level of global consumption putting serious pressure on the livestock industry [...].

The turkey ascarid, Ascaridia dissimilis, as a model genetic system

Parasitic nematodes cause significant effects on humans each year, with the most prevalent being Ascaris lumbricoides. Benzimidazoles (BZ) are the most widely used anthelmintic drug in humans, and although the biology of resistance to this drug class is understood in some species, resistance is poorly characterized in ascarids. Models such as Caenorhabditis elegans were essential in developing our current understanding of BZ resistance, but more closely related model nematodes are needed to understand resistance in ascarids. Here, we propose a new ascarid model species that infects turkeys, Ascaridia dissimilis, to develop a better understanding of BZ resistance.

Pharmacokinetic and pharmacodynamic considerations for treating sarcoptic mange with cross-relevance to Australian wildlife

Sarcoptes scabiei is the microscopic burrowing mite responsible for sarcoptic mange, which is reported in approximately 150 mammalian species. In Australia, sarcoptic mange affects a number of native and introduced wildlife species, is particularly severe in bare-nosed wombats (Vombatus ursinus) and an emerging issue in koala and quenda. There are a variety of acaricides available for the treatment of sarcoptic mange which are generally effective in eliminating mites from humans and animals in captivity. In wild populations, effective treatment is challenging, and concerns exist regarding safety, efficacy and the potential emergence of acaricide resistance. There are risks where acaricides are used intensively or inadequately, which could adversely affect treatment success rates as well as animal welfare. While reviews on epidemiology, treatment strategies, and pathogenesis of sarcoptic mange in wildlife are available, there is currently no review evaluating the use of specific acaricides in the context of their pharmacokinetic and pharmacodynamic properties, and subsequent likelihood of emerging drug resistance, particularly for Australian wildlife. This review critically evaluates acaricides that have been utilised to treat sarcoptic mange in wildlife, including dosage forms and routes, pharmacokinetics, mode of action and efficacy. We also highlight the reports of resistance of S. scabiei to acaricides, including clinical and in vitro observations.

Transcriptome profile of halofuginone resistant and sensitive strains of Eimeria tenella

The antiparasitic drug halofuginone is important for controlling apicomplexan parasites. However, the occurrence of halofuginone resistance is a major obstacle for it to the treatment of apicomplexan parasites. Current studies have identified the molecular marker and drug resistance mechanisms of halofuginone in Plasmodium falciparum. In this study, we tried to use transcriptomic data to explore resistance mechanisms of halofuginone in apicomplexan parasites of the genus Eimeria (Apicomplexa: Eimeriidae). After halofuginone treatment of E. tenella parasites, transcriptome analysis was performed using samples derived from both resistant and sensitive strains. In the sensitive group, DEGs associated with enzymes were significantly downregulated, whereas the DNA damaging process was upregulated after halofuginone treatment, revealing the mechanism of halofuginone-induced parasite death. In addition, 1,325 differentially expressed genes (DEGs) were detected between halofuginone resistant and sensitive strains, and the DEGs related to translation were significantly downregulated after halofuginone induction. Overall, our results provide a gene expression profile for further studies on the mechanism of halofuginone resistance in E. tenella.

Assembly and Analysis of Haemonchus contortus Transcriptome as a Tool for the Knowledge of Ivermectin Resistance Mechanisms

Haemonchus contortus (Hc) is an important parasitic nematode of small ruminants. In this study we assembled the transcriptome of Hc as a model to contribute to the knowledge about the profile of the differential gene expression between two Mexican Hc strains under different anthelmintic resistance statuses, one susceptible and the other resistant to ivermectin (IVMs and IVMr, respectively), in order to improve and/or to have new strategies of control and diagnosis. The transcript sequence reads were assembled and annotated. Overall, ~127 Mbp were assembled and distributed into 77,422 transcript sequences, and 4394 transcripts of the de novo transcriptome were matched base on at least one of the following criteria: (1) Phylum Nemathelminthes and Platyhelminthes, important for animal health care, and (2) ≥55% of sequence identity with other organisms. The gene ontology (GO) enrichment analysis (GOEA) was performed to study the level of gene regulation to IVMr and IVMs strains using Log Fold Change (LFC) filtering values ≥ 1 and ≥ 2. The upregulated-displayed genes obtained via GOEA were: 1993 (for LFC ≥ 1) and 1241 (for LFC ≥ 2) in IVMr and 1929 (for LFC ≥ 1) and 835 (for LFC ≥ 2) in IVMs. The enriched GO terms upregulated per category identified the intracellular structure, intracellular membrane-bounded organelle and integral component of the cell membrane as some principal cellular components. Meanwhile, efflux transmembrane transporter activity, ABC-type xenobiotic transporter activity and ATPase-coupled transmembrane transporter activity were associated with molecular function. Responses to nematicide activity, pharyngeal pumping and positive regulation of synaptic assembly were classified as biological processes that might be involved in events related to the anthelmintic resistance (AR) and nematode biology. The filtering analysis of both LFC values showed similar genes related to AR. This study deepens our knowledge about the mechanisms behind the processes of H. contortus in order to help in tool production and to facilitate the reduction of AR and promote the development of other control strategies, such as anthelmintic drug targets and vaccines.

A major locus confers triclabendazole resistance in Fasciola hepatica and shows dominant inheritance

Fasciola hepatica infection is responsible for substantial economic losses in livestock worldwide and poses a threat to human health in endemic areas. The mainstay of control in livestock and the only drug licenced for use in humans is triclabendazole (TCBZ). TCBZ resistance has been reported on every continent and threatens effective control of fasciolosis in many parts of the world. To date, understanding the genetic mechanisms underlying TCBZ resistance has been limited to studies of candidate genes, based on assumptions of their role in drug action. Taking an alternative approach, we combined a genetic cross with whole-genome sequencing to localise a ~3.2Mbp locus within the 1.2Gbp F. hepatica genome that confers TCBZ resistance. We validated this locus independently using bulk segregant analysis of F. hepatica populations and showed that it is the target of drug selection in the field. We genotyped individual parasites and tracked segregation and reassortment of SNPs to show that TCBZ resistance exhibits Mendelian inheritance and is conferred by a dominant allele. We defined gene content within this locus to pinpoint genes involved in membrane transport, (e.g. ATP-binding cassette family B, ABCB1), transmembrane signalling and signal transduction (e.g. GTP-Ras-adenylyl cyclase and EGF-like protein), DNA/RNA binding and transcriptional regulation (e.g. SANT/Myb-like DNA-binding domain protein) and drug storage and sequestration (e.g. fatty acid binding protein, FABP) as prime candidates for conferring TCBZ resistance. This study constitutes the first experimental cross and genome-wide approach for any heritable trait in F. hepatica and is key to understanding the evolution of drug resistance in Fasciola spp. to inform deployment of efficacious anthelmintic treatments in the field.

Genomic landscape of drug response reveals mediators of anthelmintic resistance

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field.

Genomic signatures of selection associated with benzimidazole drug treatments in Haemonchus contortus field populations

Genome-wide methods offer a powerful approach to detect signatures of drug selection. However, limited availability of suitable reference genomes and the difficulty of obtaining field populations with well-defined, distinct drug treatment histories mean there is little information on the signatures of selection in parasitic nematodes and on how best to detect them. This study addresses these knowledge gaps by using field populations of Haemonchus contortus with well-defined benzimidazole treatment histories, leveraging a recently completed chromosomal-scale reference genome assembly. We generated a panel of 49,393 genomic markers to genotype 20 individual adult worms from each of four H. contortus populations: two from closed sheep flocks with an approximate 20 year history of frequent benzimidazole treatment, and two populations with a history of little or no treatment. Sampling occurred in the same geographical region to limit genetic differentiation and maximise the detection sensitivity. A clear signature of selection was detected on chromosome I, centred on the isotype-1 β-tubulin gene. Two additional, but weaker, signatures of selection were detected; one near the middle of chromosome I spanning 3.75 Mbp and 259 annotated genes, and one on chromosome II spanning a region of 3.3 Mbp and 206 annotated genes, including the isotype-2 β-tubulin locus. We also assessed how sensitivity was impacted by sequencing depth, worm number, and pooled versus individual worm sequence data. This study provides the first known direct genome-wide evidence for any parasitic nematode, that the isotype-1 β-tubulin gene is quantitatively the single most important benzimidazole resistance locus. It also identified two additional genomic regions that likely contain benzimidazole resistance loci of secondary importance. This study provides an experimental framework to maximise the power of genome-wide approaches to detect signatures of selection driven by anthelmintic drug treatments in field populations of parasitic nematodes.