Genomic introgression mapping of field-derived multiple-anthelmintic resistance in Teladorsagia circumcincta

Next-generation sequencing technologies for helminth diagnostics and surveillance in ruminants: shifting diagnostic barriers

Helminth infections in grazing ruminants are a major issue for livestock farming globally, but are unavoidable in outdoor grazing systems and must be effectively managed to avoid deleterious effects to animal health, and productivity. Next-generation sequencing (NGS) technologies are transforming our understanding of the genetic basis of anthelmintic resistance (AR) and epidemiological studies of ruminant gastrointestinal parasites. They also have the potential to not only help develop and validate molecular diagnostic tests but to be directly used in routine diagnostics integrating species-specific identification and AR into a single test. Here, we review how these developments have opened the pathway for the development of multi-AR and multispecies identification in a single test, with widespread implications for sustainable livestock farming for the future.

Functional validation of novel levamisole resistance marker S168T in Haemonchus contortus

Recently, a S168T variant in the acetylcholine receptor subunit ACR-8 was associated with levamisole resistance in the parasitic helminth Haemonchus contortus. Here, we used the Xenopus laevis oocyte expression system and two-electrode voltage-clamp electrophysiology to measure the functional impact of this S168T variant on the H. contortus levamisole-sensitive acetylcholine receptor, L-AChR-1.1. Expression of the ACR-8 S168T variant significantly reduced the current amplitude elicited by levamisole compared to acetylcholine, with levamisole changing from a full to partial agonist on the recombinant L-AChR. Functional validation of the S168T mutation on modulating levamisole activity at the receptor level highlights its critical importance as both a mechanism and a marker of levamisole resistance.

The identification of small molecule inhibitors with anthelmintic activities that target conserved proteins among ruminant gastrointestinal nematodes

Gastrointestinal nematode (GIN) infections are a major concern for the ruminant industry worldwide and result in significant production losses. Naturally occurring polyparasitism and increasing drug resistance that potentiate disease outcomes are observed among the most prevalent GINs of veterinary importance. Within the five major taxonomic clades, clade Va represents a group of GINs that predominantly affect the abomasum or small intestine of ruminants. However, the development of effective broad-spectrum anthelmintics against ruminant clade Va GINs has been challenged by a lack of comprehensive druggable genome resources. Here, we first assembled draft genomes for three clade Va species (Cooperia oncophora, Trichostrongylus colubriformis, and Ostertagia ostertagi) and compared them with closely related ruminant GINs. Genome-wide phylogenetic reconstruction showed a relationship among ruminant GINs structured by taxonomic classification. Orthogroup (OG) inference and functional enrichment analyses identified 220 clade Va-specific and Va-conserved OGs, enriched for functions related to cell cycle and cellular senescence. Further transcriptomic analysis identified 61 taxonomically and functionally conserved clade Va OGs that may function as drug targets for new broad-spectrum anthelmintics. Chemogenomic screening identified 11 compounds targeting homologs of these OGs, thus having potential anthelmintic activity. In in vitro phenotypic assays, three kinase inhibitors (digitoxigenin, K-252a, and staurosporine) exhibited broad-spectrum anthelmintic activities against clade Va GINs by obstructing the motility of exsheathed L3 (xL3) or molting of xL3 to L4. These results demonstrate valuable applications of the new ruminant GIN genomes in gaining better insights into their life cycles and offer a contemporary approach to discovering the next generation of anthelmintics.IMPORTANCEGastrointestinal nematode (GIN) infections in ruminants are caused by parasites that inhibit normal function in the digestive tract of cattle, sheep, and goats, thereby causing morbidity and mortality. Coinfection and increasing drug resistance to current therapeutic agents will continue to worsen disease outcomes and impose significant production losses on domestic livestock producers worldwide. In combination with ongoing therapeutic efforts, advancing the discovery of new drugs with novel modes of action is critical for better controlling GIN infections. The significance of this study is in assembling and characterizing new GIN genomes of Cooperia oncophora, Ostertagia ostertagi, and Trichostrongylus colubriformis for facilitating a multi-omics approach to identify novel, biologically conserved drug targets for five major GINs of veterinary importance. With this information, we were then able to demonstrate the potential of commercially available compounds as new anthelmintics.

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.

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.

Mapping resistance-associated anthelmintic interactions in the model nematode Caenorhabditis elegans

Parasitic nematodes infect billions of people and are mainly controlled by anthelmintic mass drug administration (MDA). While there are growing efforts to better understand mechanisms of anthelmintic resistance in human and animal populations, it is unclear how resistance mechanisms that alter susceptibility to one drug affect the interactions and efficacy of drugs used in combination. Mutations that alter drug permeability across primary nematode barriers have been identified as potential resistance mechanisms using the model nematode Caenorhabditis elegans. We leveraged high-throughput assays in this model system to measure altered anthelmintic susceptibility in response to genetic perturbations of potential cuticular, amphidial, and alimentary routes of drug entry. Mutations in genes associated with these tissue barriers differentially altered susceptibility to the major anthelmintic classes (macrocyclic lactones, benzimidazoles, and nicotinic acetylcholine receptor agonists) as measured by animal development. We investigated two-way anthelmintic interactions across C. elegans genetic backgrounds that confer resistance or hypersensitivity to one or more drugs. We observe that genetic perturbations that alter susceptibility to a single drug can shift the drug interaction landscape and lead to the appearance of novel synergistic and antagonistic interactions. This work establishes a framework for investigating combinatorial therapies in model nematodes that can potentially be translated to amenable parasite species.

Importance of ABC Transporters in the Survival of Parasitic Nematodes and the Prospect for the Development of Novel Control Strategies

ABC transporters, a family of ATP-dependent transmembrane proteins, are responsible for the active transport of a wide range of molecules across cell membranes, including drugs, toxins, and nutrients. Nematodes possess a great diversity of ABC transporters; however, only P-glycoproteins have been well-characterized compared to other classes. The ABC transport proteins have been implicated in developing resistance to various classes of anthelmintic drugs in parasitic nematodes; their role in plant and human parasitic nematodes still needs further investigation. Therefore, ABC transport proteins offer a potential opportunity to develop nematode control strategies. Multidrug resistance inhibitors are becoming more attractive for controlling nematodes due to their potential to increase drug efficacy in two ways: (i) by limiting drug efflux from nematodes, thereby increasing the amount of drug that reaches its target site, and (ii) by reducing drug excretion by host animals, thereby enhancing drug bioavailability. This article reviews the role of ABC transporters in the survival of parasitic nematodes, including the genes involved, their regulation and physiological roles, as well as recent developments in their characterization. It also discusses the association of ABC transporters with anthelmintic resistance and the possibility of targeting them with next-generation inhibitors or nutraceuticals (e.g., polyphenols) to control parasitic infections.

Mapping resistance-associated anthelmintic interactions in the model nematode Caenorhabditis elegans

Parasitic nematodes infect billions of people and are mainly controlled by anthelmintic mass drug administration (MDA). While there are growing efforts to better understand mechanisms of anthelmintic resistance in human and animal populations, it is unclear how resistance mechanisms that alter susceptibility to one drug affect the interactions and efficacy of drugs used in combination. Mutations that alter drug permeability across primary nematode barriers have been identified as potential resistance mechanisms using the model nematode Caenorhabditis elegans. We leveraged high-throughput assays in this model system to measure altered anthelmintic susceptibility in response to genetic perturbations of potential cuticular, amphidial, and alimentary routes of drug entry. Mutations in genes associated with these tissue barriers differentially altered susceptibility to the major anthelmintic classes (macrocyclic lactones, benzimidazoles, and nicotinic acetylcholine receptor agonists) as measured by animal development. We investigated two-way anthelmintic interactions across C. elegans genetic backgrounds that confer resistance or hypersensitivity to one or more drugs. We observe that genetic perturbations that alter susceptibility to a single drug can shift the drug interaction landscape and lead to the appearance of novel synergistic and antagonistic interactions. This work establishes a framework for investigating combinatorial therapies in model nematodes that can potentially be translated to amenable parasite species.

Genomic analysis of two phlebotomine sand fly vectors of Leishmania from the New and Old World

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

Intraspecific diversity in the mechanisms underlying abamectin resistance in a cosmopolitan pest

Pesticide resistance relies on a myriad of mechanisms, ranging from single mutations to a complex and polygenic architecture, and it involves mechanisms such as target-site insensitivity, metabolic detoxification, or a combination of these, with either additive or synergistic effects. Several resistance mechanisms against abamectin, a macrocyclic lactone widely used in crop protection, have been reported in the cosmopolitan pest Tetranychus urticae. However, it has been shown that a single mechanism cannot account for the high levels of abamectin resistance found across different mite populations. Here, we used experimental evolution combined with bulked segregant analyses to map quantitative trait loci (QTL) associated with abamectin resistance in two genetically unrelated populations of T. urticae. In these two independent QTL mapping experiments, three and four QTLs were identified, of which three were shared between experiments. Shared QTLs contained genes encoding subunits of the glutamate-gated chloride channel (GluCl) and harboured previously reported mutations, including G314D in GluCl1 and G326E in GluCl3, but also novel resistance candidate loci, including DNA helicases and chemosensory receptors. Surprisingly, the fourth QTL, present only in only one of the experiments and thus unique for one resistant parental line, revealed a non-functional variant of GluCl2, suggesting gene knock-out as resistance mechanism. Our study uncovers the complex basis of abamectin resistance, and it highlights the intraspecific diversity of genetic mechanisms underlying resistance in a cosmopolitan pest.

Variation in anthelmintic responses are driven by genetic differences among diverse C. elegans wild strains

Treatment of parasitic nematode infections in humans and livestock relies on a limited arsenal of anthelmintic drugs that have historically reduced parasite burdens. However, anthelmintic resistance (AR) is increasing, and little is known about the molecular and genetic causes of resistance for most drugs. The free-living roundworm Caenorhabditis elegans has proven to be a tractable model to understand AR, where studies have led to the identification of molecular targets of all major anthelmintic drug classes. Here, we used genetically diverse C. elegans strains to perform dose-response analyses across 26 anthelmintic drugs that represent the three major anthelmintic drug classes (benzimidazoles, macrocyclic lactones, and nicotinic acetylcholine receptor agonists) in addition to seven other anthelmintic classes. First, we found that C. elegans strains displayed similar anthelmintic responses within drug classes and significant variation across drug classes. Next, we compared the effective concentration estimates to induce a 10% maximal response (EC10) and slope estimates of each dose-response curve of each strain to the laboratory reference strain, which enabled the identification of anthelmintics with population-wide differences to understand how genetics contribute to AR. Because genetically diverse strains displayed differential susceptibilities within and across anthelmintics, we show that C. elegans is a useful model for screening potential nematicides before applications to helminths. Third, we quantified the levels of anthelmintic response variation caused by genetic differences among individuals (heritability) to each drug and observed a significant correlation between exposure closest to the EC10 and the exposure that exhibited the most heritable responses. These results suggest drugs to prioritize in genome-wide association studies, which will enable the identification of AR genes.

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

Chromosome-length genome assembly of Teladorsagia circumcincta - a globally important helminth parasite in livestock

BACKGROUND

Gastrointestinal (GIT) helminthiasis is a global problem that affects livestock health, especially in small ruminants. One of the major helminth parasites of sheep and goats, Teladorsagia circumcincta, infects the abomasum and causes production losses, reductions in weight gain, diarrhoea and, in some cases, death in young animals. Control strategies have relied heavily on the use of anthelmintic medication but, unfortunately, T. circumcincta has developed resistance, as have many helminths. Vaccination offers a sustainable and practical solution, but there is no commercially available vaccine to prevent Teladorsagiosis. The discovery of new strategies for controlling T. circumcincta, such as novel vaccine targets and drug candidates, would be greatly accelerated by the availability of better quality, chromosome-length, genome assembly because it would allow the identification of key genetic determinants of the pathophysiology of infection and host-parasite interaction. The available draft genome assembly of T. circumcincta (GCA_002352805.1) is highly fragmented and thus impedes large-scale investigations of population and functional genomics.

RESULTS

We have constructed a high-quality reference genome, with chromosome-length scaffolds, by purging alternative haplotypes from the existing draft genome assembly and scaffolding the result using chromosome conformation, capture-based, in situ Hi-C technique. The improved (Hi-C) assembly resulted in six chromosome-length scaffolds with length ranging from 66.6 Mbp to 49.6 Mbp, 35% fewer sequences and reduction in size. Substantial improvements were also achieved in both the values for N50 (57.1 Mbp) and L50 (5 Mbp). A higher and comparable level of genome and proteome completeness was achieved for Hi-C assembly on BUSCO parameters. The Hi-C assembly had a greater synteny and number of orthologs with a closely related nematode, Haemonchus contortus.

CONCLUSION

This improved genomic resource is suitable as a foundation for the identification of potential targets for vaccine and drug development.

Transcriptomics of ivermectin response in Caenorhabditis elegans: Integrating abamectin quantitative trait loci and comparison to the Ivermectin-exposed DA1316 strain

Parasitic nematodes pose a significant threat to human and animal health, as well as cause economic losses in the agricultural sector. The use of anthelmintic drugs, such as Ivermectin (IVM), to control these parasites has led to widespread drug resistance. Identifying genetic markers of resistance in parasitic nematodes can be challenging, but the free-living nematode Caenorhabditis elegans provides a suitable model. In this study, we aimed to analyze the transcriptomes of adult C. elegans worms of the N2 strain exposed to the anthelmintic drug Ivermectin (IVM), and compare them to those of the resistant strain DA1316 and the recently identified Abamectin Quantitative Trait Loci (QTL) on chromosome V. We exposed pools of 300 adult N2 worms to IVM (10-7 and 10-8 M) for 4 hours at 20°C, extracted total RNA and sequenced it on the Illumina NovaSeq6000 platform. Differentially expressed genes (DEGs) were determined using an in-house pipeline. The DEGs were compared to genes from a previous microarray study on IVM-resistant C. elegans and Abamectin-QTL. Our results revealed 615 DEGs (183 up-regulated and 432 down-regulated genes) from diverse gene families in the N2 C. elegans strain. Of these DEGs, 31 overlapped with genes from IVM-exposed adult worms of the DA1316 strain. We identified 19 genes, including the folate transporter (folt-2) and the transmembrane transporter (T22F3.11), which exhibited an opposite expression in N2 and the DA1316 strain and were deemed potential candidates. Additionally, we compiled a list of potential candidates for further research including T-type calcium channel (cca-1), potassium chloride cotransporter (kcc-2), as well as other genes such as glutamate-gated channel (glc-1) that mapped to the Abamectin-QTL.

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.

The equine ascarids: resuscitating historic model organisms for modern purposes

The equine ascarids, Parascaris spp., are important nematode parasites of juvenile horses and were historically model organisms in the field of cell biology, leading to many important discoveries, and are used for the study of chromatin diminution. In veterinary parasitology, Parascaris spp. are important not only because they can cause clinical disease in young horses but also because they are the only ascarid parasites to have developed widespread anthelmintic resistance. Despite this, much of the general biology and mechanisms of anthelmintic resistance are poorly understood. This review condenses known basic biological information and knowledge on the mechanisms of anthelmintic resistance in Parascaris spp., highlighting the importance of foundational research programs. Although two variants of this parasite were recognized based on the number of chromosomes in the 1870s and suggested to be two species in 1890, one of these, P. univalens, appears to have been largely forgotten in the veterinary scientific literature over the past 100 years. We describe how this omission has had a century-long effect on nomenclature and data analysis in the field, highlighting the importance of proper specimen identification in public repositories. A summary of important basic biology, including life cycle, in vitro maintenance, and immunology, is given, and areas of future research for the improvement of knowledge and development of new systems are given. Finally, the limited knowledge regarding anthelmintic resistance in Parascaris spp. is summarized, along with caution regarding assumptions that resistance mechanisms can be applied across clades.

Transgenic Expression of Haemonchus contortus Cytochrome P450 Hco-cyp-13A11 Decreases Susceptibility to Particular but Not All Macrocyclic Lactones in the Model Organism Caenorhabditis elegans

The number of reported macrocyclic lactones (ML) resistance cases across all livestock hosts is steadily increasing. Different studies in the parasitic nematode Haemonchus contortus assume the participation of cytochrome P450s (Cyps) enzymes in ML resistance. Still, functional data about their individual contribution to resistance or substrate specificity is missing. Via microinjection, transgenic Caenorhabditis elegans expressing HCON_00141052 (transgene-Hco-cyp-13A11) from extrachromosomal arrays were generated. After 24 h of exposure to different concentrations of ivermectin (IVM), ivermectin aglycone (IVMa), selamectin (SEL), doramectin (DRM), eprinomectin (EPR), and moxidectin (MOX), motility assays were performed to determine the impact of the H. contortus Cyp to the susceptibility of the worms against each ML. While transgene-Hco-cyp-13A11 significantly decreased susceptibility to IVM (four-fold), IVMa (2-fold), and SEL (3-fold), a slight effect for DRM and no effect for MOX, and EPR was observed. This substrate specificity of Hco-cyp-13A11 could not be explained by molecular modeling and docking studies. Hco-Cyp-13A11 molecular models were obtained for alleles from isolates with different resistance statuses. Although 14 amino acid polymorphisms were detected, none was resistance specific. In conclusion, Hco-cyp-13A11 decreased IVM, IVMa, and SEL susceptibility to a different extent, but its potential impact on ML resistance is not driven by polymorphisms.

Allele specific PCR for a major marker of levamisole resistance in Haemonchus contortus

Haemonchus contortus is a haematophagous parasitic nematode that infects small ruminants and causes significant animal health concerns and economic losses within the livestock industry on a global scale. Treatment primarily depends on broad-spectrum anthelmintics, however, resistance is established or rapidly emerging against all major drug classes. Levamisole (LEV) remains an important treatment option for parasite control, as resistance to LEV is less prevalent than to members of other major classes of anthelmintics. LEV is an acetylcholine receptor (AChR) agonist that, when bound, results in paralysis of the worm. Numerous studies implicated the AChR sub-unit, ACR-8, in LEV sensitivity and in particular, the presence of a truncated acr-8 transcript or a deletion in the acr-8 locus in some resistant isolates. Recently, a single non-synonymous SNP in acr-8 conferring a serine-to-threonine substitution (S168T) was identified that was strongly associated with LEV resistance. Here, we investigate the role of genetic variation at the acr-8 locus in a controlled genetic cross between the LEV susceptible MHco3(ISE) and LEV resistant MHco18(UGA2004) isolates of H. contortus. Using single worm PCR assays, we found that the presence of S168T was strongly associated with LEV resistance in the parental isolates and F3 progeny of the genetic cross surviving LEV treatment. We developed and optimised an allele-specific PCR assay for the detection of S168T and validated the assay using laboratory isolates and field samples that were phenotyped for LEV resistance. In the LEV-resistant field population, a high proportion (>75%) of L₃ encoded the S168T variant, whereas the variant was absent in the susceptible isolates studied. These data further support the potential role of acr-8 S168T in LEV resistance, with the allele-specific PCR providing an important step towards establishing a sensitive molecular diagnostic test for LEV resistance.

Genome-wide analysis of the response to ivermectin treatment by a Swedish field population of Haemonchus contortus

Haemonchus contortus is a pathogenic gastrointestinal nematode of small ruminants and, in part due to its capacity to develop resistance to drugs, contributes to significant losses in the animal production sector worldwide. Despite decades of research, comparatively little is known about the specific mechanism(s) driving resistance to drugs such as ivermectin in this species. Here we describe a genome-wide approach to detect evidence of selection by ivermectin treatment in a field population of H. contortus from Sweden, using parasites sampled from the same animals before and seven days after ivermectin exposure followed by whole-genome sequencing. Despite an 89% reduction in parasites recovered after treatment measured by the fecal egg count reduction test, the surviving population was highly genetically similar to the population before treatment, suggesting that resistance has likely evolved over time and that resistance alleles are present on diverse haplotypes. Pairwise gene and SNP frequency comparisons indicated the highest degree of differentiation was found at the terminal end of chromosome 4, whereas the most striking difference in nucleotide diversity was observed in a region on chromosome 5 previously reported to harbor a major quantitative trait locus involved in ivermectin resistance. These data provide novel insight into the genome-wide effect of ivermectin selection in a field population as well as confirm the importance of the previously established quantitative trait locus in the development of resistance to ivermectin.

Anthelmintic resistance in ruminants: challenges and solutions

Anthelmintic resistance (AR) is a growing concern for effective parasite control in farmed ruminants globally. Combatting AR will require intensified and integrated research efforts in the development of innovative diagnostic tests to detect helminth infections and AR, sustainable anthelmintic treatment strategies and the development of complementary control approaches such as vaccination and plant-based control. It will also require a better understanding of socio-economic drivers of anthelmintic treatment decisions, in order to support a behavioural shift and develop targeted communication strategies that promote the uptake of evidence-based sustainable solutions. Here, we review the state-of-the-art in these different fields of research activity related to AR in helminths of livestock ruminants in Europe and beyond. We conclude that in the advent of new challenges and solutions emerging from continuing spread of AR and intensified research efforts, respectively, there is a strong need for transnational multi-actor initiatives. These should involve all key stakeholders to develop indicators of infection and sustainable control, set targets and promote good practices to achieve them.

A journey through 50 years of research relevant to the control of gastrointestinal nematodes in ruminant livestock and thoughts on future directions

This review article provides an historical perspective on some of the major research advances of relevance to ruminant livestock gastrointestinal nematode control over the last 50 years. Over this period, gastrointestinal nematode control has been dominated by the use of broad-spectrum anthelmintic drugs. Whilst this has provided unprecedented levels of successful control for many years, this approach has been gradually breaking down for more than two decades and is increasingly unsustainable which is due, at least in part, to the emergence of anthelmintic drug resistance and a number of other factors discussed in this article. We first cover the remarkable success story of the discovery and development of broad-spectrum anthelmintic drugs, the changing face of anthelmintic drug discovery research and the emergence of anthelmintic resistance. This is followed by a review of some of the major advances in the increasingly important area of non-pharmaceutical gastrointestinal nematode control including immunology and vaccine development, epidemiological modelling and some of the alternative control strategies such as breeding for host resistance, refugia-based methods and biological control. The last 50 years have witnessed remarkable innovation and success in research aiming to improve ruminant livestock gastrointestinal nematode control, particularly given the relatively small size of the research community and limited funding. In spite of this, the growing global demand for livestock products, together with the need to maximise production efficiencies, reduce environmental impacts and safeguard animal welfare - as well as specific challenges such as anthelmintic drug resistance and climate change- mean that gastrointestinal nematode researchers will need to be as innovative in the next 50 years as in the last.

Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration

Control and elimination of the parasitic disease schistosomiasis relies on mass administration of praziquantel. Whilst these programmes reduce infection prevalence and intensity, their impact on parasite transmission and evolution is poorly understood. Here we examine the genomic impact of repeated mass drug administration on Schistosoma mansoni populations with documented reduced praziquantel efficacy. We sequenced whole-genomes of 198 S. mansoni larvae from 34 Ugandan children from regions with contrasting praziquantel exposure. Parasites infecting children from Lake Victoria, a transmission hotspot, form a diverse panmictic population. A single round of treatment did not reduce this diversity with no apparent population contraction caused by long-term praziquantel use. We find evidence of positive selection acting on members of gene families previously implicated in praziquantel action, but detect no high frequency functionally impactful variants. As efforts to eliminate schistosomiasis intensify, our study provides a foundation for genomic surveillance of this major human parasite.

Comparative evaluation of different molecular methods for DNA extraction from individual Teladorsagia circumcincta nematodes

Background

The purpose of this study was to develop a reliable DNA extraction protocol to use on individual Teladorsagia circumcincta nematode specimens to produce high quality DNA for genome sequencing and phylogenetic analysis. Pooled samples have been critical in providing the groundwork for T. circumcincta genome construction, but there is currently no standard method for extracting high-quality DNA from individual nematodes. 11 extraction kits were compared based on DNA quality, yield, and processing time.

Results

11 extraction protocols were compared, and the concentration and purity of the extracted DNA was quantified. Median DNA concentration among all methods measured on NanoDrop 2000™ ranged between 0.45–11.5 ng/μL, and on Qubit™ ranged between undetectable – 0.962 ng/μL. Median A260/280 ranged between 0.505–3.925, and median A260/230 ranged − 0.005 – 1.545. Larval exsheathment to remove the nematode cuticle negatively impacted DNA concentration and purity.

Conclusions

A Schistosoma sp. DNA extraction method was determined as most suitable for individual T. circumcincta nematode specimens due to its resulting DNA concentration, purity, and relatively fast processing time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00695-6.

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses

Parasitic nematodes cause a massive worldwide burden on human health along with a loss of livestock and agriculture productivity. Anthelmintics have been widely successful in treating parasitic nematodes. However, resistance is increasing, and little is known about the molecular and genetic causes of resistance for most of these drugs. The free-living roundworm Caenorhabditis elegans provides a tractable model to identify genes that underlie resistance. Unlike parasitic nematodes, C. elegans is easy to maintain in the laboratory, has a complete and well annotated genome, and has many genetic tools. Using a combination of wild isolates and a panel of recombinant inbred lines constructed from crosses of two genetically and phenotypically divergent strains, we identified three genomic regions on chromosome V that underlie natural differences in response to the macrocyclic lactone (ML) abamectin. One locus was identified previously and encodes an alpha subunit of a glutamate-gated chloride channel (glc-1). Here, we validate and narrow two novel loci using near-isogenic lines. Additionally, we generate a list of prioritized candidate genes identified in C. elegans and in the parasite Haemonchus contortus by comparison of ML resistance loci. These genes could represent previously unidentified resistance genes shared across nematode species and should be evaluated in the future. Our work highlights the advantages of using C. elegans as a model to better understand ML resistance in parasitic nematodes.

Pharyngeal Pumping and Tissue-Specific Transgenic P-Glycoprotein Expression Influence Macrocyclic Lactone Susceptibility in Caenorhabditis elegans

Macrocyclic lactones (MLs) are widely used drugs to treat and prevent parasitic nematode infections. In many nematode species including a major pathogen of foals, Parascaris univalens, resistance against MLs is widespread, but the underlying resistance mechanisms and ML penetration routes into nematodes remain unknown. Here, we examined how the P-glycoprotein efflux pumps, candidate genes for ML resistance, can modulate drug susceptibility and investigated the role of active drug ingestion for ML susceptibility in the model nematode Caenorhabditis elegans. Wildtype or transgenic worms, modified to overexpress P. univalens PGP-9 (Pun-PGP-9) at the intestine or epidermis, were incubated with ivermectin or moxidectin in the presence (bacteria or serotonin) or absence (no specific stimulus) of pharyngeal pumping (PP). Active drug ingestion by PP was identified as an important factor for ivermectin susceptibility, while moxidectin susceptibility was only moderately affected. Intestinal Pun-PGP-9 expression elicited a protective effect against ivermectin and moxidectin only in the presence of PP stimulation. Conversely, epidermal Pun-PGP-9 expression protected against moxidectin regardless of PP and against ivermectin only in the absence of active drug ingestion. Our results demonstrate the role of active drug ingestion by nematodes for susceptibility and provide functional evidence for the contribution of P-glycoproteins to ML resistance in a tissue-specific manner.

Susceptible trichostrongyloid species mask presence of benzimidazole-resistant Haemonchus contortus in cattle

Background

Benzimidazole (BZ) anthelmintics are widely used to control infections with parasitic nematodes, but BZ resistance is an emerging threat among several nematode species infecting humans and animals. In Sudan, BZ-resistant Haemonchus contortus populations were recently reported in goats in South Darfur State. The objective of this study was to collect data regarding the situation of BZ resistance in cattle parasitic nematodes in South Darfur using phenotypic and molecular approaches, besides providing some epidemiological data on nematodes in cattle.

Methods

The faecal egg count reduction test and the egg hatch test (EHT) were used to evaluate benzimidazole efficacy in cattle nematodes in five South Darfur study areas: Beleil, Kass, Nyala, Rehed Al-Birdi and Tulus. Genomic DNA was extracted from pools of third-stage larvae (L3) (n = 40) during trials, before and after treatment, and pools of adult male Haemonchus spp. (n = 18) from abattoirs. The polymorphisms F167Y, E198A and F200Y in isotype 1 β-tubulin genes of H. contortus and H. placei were analysed using Sanger and pyrosequencing.

Results

Prevalence of gastro-intestinal helminths in cattle was 71% (313/443). Reduced albendazole faecal egg count reduction efficacy was detected in three study areas: Nyala (93.7%), Rehed Al-Birdi (89.7%) and Tulus (88.2%). In the EHT, EC₅₀ values of these study areas ranged between 0.032 and 0.037 µg/ml thiabendazole. Genus-specific PCRs detected the genera Haemonchus, Trichostrongylus and Cooperia in L3 samples collected after albendazole treatment. Sanger sequencing followed by pyrosequencing assays did not detect elevated frequencies of known BZ resistance-associated alleles in codon F167Y, E198A and F200Y in isotype 1 β-tubulin gene of H. placei (≤ 11.38%). However, polymorphisms were detected in H. contortus and in samples with mixed infections with H. contortus and H. placei at codon 198, including E198L (16/58), E198V (2/58) and potentially E198Stop (1/58). All pooled L3 samples post-albendazole treatment (n = 13) were identified as H. contortus with an E198L substitution at codon 198.

Conclusions

To the knowledge of the authors, this is the first report of reduced albendazole efficacy in cattle in Sudan and is the first study describing an E198L substitution in phenotypically BZ-resistant nematodes collected from cattle.

Challenges and opportunities for the adoption of molecular diagnostics for anthelmintic resistance

Anthelmintic resistance is a significant threat to livestock production systems worldwide and is emerging as an important issue in companion animal parasite management. It is also an emerging concern for the control of human soil-transmitted helminths and filaria. An important aspect of managing anthelmintic resistance is the ability to utilise diagnostic tests to detect its emergence at an early stage. In host-parasite systems where resistance is already widespread, diagnostics have a potentially important role in determining those drugs that remain the most effective. The development of molecular diagnostics for anthelmintic resistance is one focus of the Consortium for Anthelmintic Resistance and Susceptibility (CARS) group. The present paper reflects discussions of this issue that occurred at the most recent meeting of the group in Wisconsin, USA, in July 2019. We compare molecular resistance diagnostics with in vivo and in vitro phenotypic methods, and highlight the advantages and disadvantages of each. We assess whether our knowledge on the identity of molecular markers for resistance towards the different drug classes is sufficient to provide some expectation that molecular tests for field use may be available in the short-to-medium term. We describe some practical aspects of such tests and how our current capabilities compare to the requirements of an 'ideal' test. Finally, we describe examples of drug class/parasite species interactions that provide the best opportunity for commercial use of molecular tests in the near future. We argue that while such prototype tests may not satisfy the requirements of an 'ideal' test, their potential to provide significant advances over currently-used phenotypic methods warrants their development as field diagnostics.

Ivermectin: An Anthelmintic, an Insecticide, and Much More

Here we tell the story of ivermectin, describing its anthelmintic and insecticidal actions and recent studies that have sought to reposition ivermectin for the treatment of other diseases that are not caused by helminth and insect parasites. The standard theory of its anthelmintic and insecticidal mode of action is that it is a selective positive allosteric modulator of glutamate-gated chloride channels found in nematodes and insects. At higher concentrations, ivermectin also acts as an allosteric modulator of ion channels found in host central nervous systems. In addition, in tissue culture, at concentrations higher than anthelmintic concentrations, ivermectin shows antiviral, antimalarial, antimetabolic, and anticancer effects. Caution is required before extrapolating from these preliminary repositioning experiments to clinical use, particularly for Covid-19 treatment, because of the high concentrations of ivermectin used in tissue-culture experiments.

Characterization of the Complete Mitochondrial Genomes of Two Sibling Species of Parasitic Roundworms, Haemonchus contortus and Teladorsagia circumcincta

Haemonchus contortus and Teladorsagia circumcincta are among the two most pathogenic internal parasitic nematodes infecting small ruminants, such as sheep and goats, and are a global animal health issue. Accurate identification and delineation of Haemonchidae species is essential for development of diagnostic and control strategies with high resolution for Trichostrongyloidea infection in ruminants. Here, we describe in detail and compare the complete mitochondrial (mt) genomes of the New Zealand H. contortus and T. circumcincta field strains to improve our understanding of species- and strain-level evolution in these closely related roundworms. In the present study, we performed extensive comparative bioinformatics analyses on the recently sequenced complete mt genomes of the New Zealand H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP field strains. Amino acid sequences inferred from individual genes of each of the two mt genomes were compared, concatenated and subjected to phylogenetic analysis using Bayesian inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP). The AT-rich mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP are 14,001 bp (A+T content of 77.4%) and 14,081 bp (A+T content of 77.3%) in size, respectively. All 36 of the typical nematode mt genes are transcribed in the forward direction in both species and comprise of 12 protein-encoding genes (PCGs), 2 ribosomal RNA (rrn) genes, and 22 transfer RNA (trn) genes. The secondary structures for the 22 trn genes and two rrn genes differ between H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP, however the gene arrangements of both are consistent with other Trichostrongylidea sequenced to date. Comparative analyses of the complete mitochondrial nucleotide sequences, PCGs, A+T rich and non-coding repeat regions of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP further reinforces the high levels of diversity and gene flow observed among Trichostrongylidea, and supports their potential as ideal markers for strain-level identification from different hosts and geographical regions with high resolution for future studies. The complete mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP presented here provide useful novel markers for further studies of the meta-population connectivity and the genetic mechanisms driving evolution in nematode species.

Teladorsagia circumcincta beta tubulin: the presence of the E198L polymorphism on its own is associated with benzimidazole resistance

BACKGROUND

Benzimidazole resistance is associated with isotype-1 β-tubulin gene F200Y, E198A and F167Y SNPs. In this study, the recently described polymorphism E198L was reported and analysed in Teladorsagia circumcincta.

METHODS

The benzimidazole phenotypic resistance was measured by the faecal egg count reduction test (FECRT) and the egg hatch test (EHT) using a discriminating dose (DD) in 39 sheep flocks. Around 1000 larvae collected before and after treatment were used for DNA extraction. The resistant species identified in all flocks was T. circumcincta. The resistance alleles frequencies were measured for F200Y and E198A. A 371-bp fragment of the isotype-1 β-tubulin gene was analysed, including the three codons of interest, and a new pyrosequencing assay was designed for testing E198L.

RESULTS

The percentage of resistant flocks was 35% by FECRT or 26% by EHT; however, F200Y and E198A SNPs were absent in T. circumcincta. The amplification of a 371-bp fragment confirmed the absence of F167Y and F200Y in 6 resistant flocks. Regarding codon 198, all samples after treatment carried a leucine (CTA). A pyrosequencing assay analysed the allele frequencies for the first two bases at codon 198 independently, G/C and A/T. The correlation between C and T frequencies was almost 1 (r = 0.929, P < 0.0001) and the mean value of both was calculated to measure the leucine frequency; this value ranged between 10.4-80.7% before treatment, and 82.3-92.8% after treatment. High and similar correlations were reported between the genotypic variables (C frequency, T frequency or mean of both frequencies) and phenotypic resistance (r > 0.720, P < 0.0001), although negatively associated with the FECRT and positively with the EHT. According to multivariate linear regression analysis, the T frequency was the most significant variable influencing the phenotypic resistance (FECRT or EHT; P < 0.0001). In the EHT, 67.1% of the phenotypic variability is associated with the T frequency but in the FECRT only 33.4%; therefore, the EHT using a DD seems to detect the genotypic resistance more accurately than the FECRT.

CONCLUSIONS

The E198L polymorphism can confer BZ resistance on its own in T. circumcincta.

The P-glycoprotein repertoire of the equine parasitic nematode Parascaris univalens

P-glycoproteins (Pgp) have been proposed as contributors to the widespread macrocyclic lactone (ML) resistance in several nematode species including a major pathogen of foals, Parascaris univalens. Using new and available RNA-seq data, ten different genomic loci encoding Pgps were identified and characterized by transcriptome-guided RT-PCRs and Sanger sequencing. Phylogenetic analysis revealed an ascarid-specific Pgp lineage, Pgp-18, as well as two paralogues of Pgp-11 and Pgp-16. Comparative gene expression analyses in P. univalens and Caenorhabditis elegans show that the intestine is the major site of expression but individual gene expression patterns were not conserved between the two nematodes. In P. univalens, PunPgp-9, PunPgp-11.1 and PunPgp-16.2 consistently exhibited the highest expression level in two independent transcriptome data sets. Using RNA-Seq, no significant upregulation of any Pgp was detected following in vitro incubation of adult P. univalens with ivermectin suggesting that drug-induced upregulation is not the mechanism of Pgp-mediated ML resistance. Expression and functional analyses of PunPgp-2 and PunPgp-9 in Saccharomyces cerevisiae provide evidence for an interaction with ketoconazole and ivermectin, but not thiabendazole. Overall, this study established reliable reference gene models with significantly improved annotation for the P. univalens Pgp repertoire and provides a foundation for a better understanding of Pgp-mediated anthelmintic resistance.

New codon 198 β-tubulin polymorphisms in highly benzimidazole resistant Haemonchus contortus from goats in three different states in Sudan

Background

Benzimidazole (BZ) resistance in gastrointestinal nematodes is a worldwide problem for livestock production, particularly in small ruminants. Assignment of the emergence of resistance using sensitive and reliable methods is required to adopt the correct strategies for control. In Sudan, BZ resistant Haemonchus contortus populations were recently reported in goats in South Darfur. This study aimed to provide additional data regarding albendazole efficacy and to describe the prevailing molecular BZ resistance mechanisms.

Methods

Faecal egg count reduction and egg hatch tests (EHT) were used to evaluate albendazole efficacy in three different areas of South Darfur using naturally (Rehed Al-Birdi and Tulus) and experimentally infected (Tulus and Um Dafuq) goats. Using samples from Central, East and South Darfur, pyro- and Sanger sequencing were used to detect the polymorphisms F167Y, E198A and F200Y in H. contortus isotype 1 β-tubulin in DNA extracted from pooled third-stage larval (L3) samples (n = 36) on days 0 and 10 during trials, and from pooled adult male H. contortus (treated goats, n = 14; abattoirs, n = 83) including samples from populations previously found to be resistant in South Darfur.

Results

Albendazole efficacies at 5, 7.5 and 10 mg/kg doses were 73.5–90.2% on day 14 in natural and experimental infections while 12.5 mg/kg showed > 96.6% efficacy. EC₅₀ in the EHT were 0.8 and 0.11 µg/ml thiabendazole in natural and experimental infection trials, respectively. PCRs detected Haemonchus, Trichostrongylus and Cooperia in L3 samples from albendazole-treated goats. Haemonchus contortus allele frequencies in codons 167 and 200 using pyrosequencing assays were ≤ 7.4% while codon 198 assays failed. Sanger sequencing revealed five novel polymorphisms at codon 198. Noteworthy, an E198L substitution was present in 82% of the samples (L3 and adults) including all post-treatment samples. Moreover, E198V, E198K and potentially E198I, and E198Stop were identified in a few samples.

Conclusions

To our knowledge, this is the first report of E198L in BZ resistant H. contortus and the second where this is the predominant genotype associated with resistance in any strongyle species. Since this variant cannot be quantified using pyrosequencing, the results highlight important limitations in the general applicability of pyrosequencing to quantify BZ resistance genotypes.

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

De novo Assembly of the Brugia malayi Genome Using Long Reads from a Single MinION Flowcell

Filarial nematode infections cause a substantial global disease burden. Genomic studies of filarial worms can improve our understanding of their biology and epidemiology. However, genomic information from field isolates is limited and available reference genomes are often discontinuous. Single molecule sequencing technologies can reduce the cost of genome sequencing and long reads produced from these devices can improve the contiguity and completeness of genome assemblies. In addition, these new technologies can make generation and analysis of large numbers of field isolates feasible. In this study, we assessed the performance of the Oxford Nanopore Technologies MinION for sequencing and assembling the genome of Brugia malayi, a human parasite widely used in filariasis research. Using data from a single MinION flowcell, a 90.3 Mb nuclear genome was assembled into 202 contigs with an N50 of 2.4 Mb. This assembly covered 96.9% of the well-defined B. malayi reference genome with 99.2% identity. The complete mitochondrial genome was obtained with individual reads and the nearly complete genome of the endosymbiotic bacteria Wolbachia was assembled alongside the nuclear genome. Long-read data from the MinION produced an assembly that approached the quality of a well-established reference genome using comparably fewer resources.

Genotypic characterisation of monepantel resistance in historical and newly derived field strains of Teladorsagia circumcincta

Recent reports of monepantel (MPTL) resistance in UK field isolates of Teladorsagia circumcincta has highlighted the need for a better understanding of the mechanism of MPTL-resistance in order to preserve its anthelmintic efficacy in this economically important species. Nine discrete populations of T. circumcincta were genotypically characterised; three MPTL-susceptible isolates, three experimentally selected MPTL-resistant strains and three field derived populations. Full-length Tci-mptl-1 gene sequences were generated and comparisons between the MPTL-susceptible isolates, MPTL-resistant strains and one field isolate, showed that different putative MPTL-resistance conferring mutations were present in different resistant isolates. Truncated forms of the Tci-mptl-1 gene were also observed. The genetic variability of individual larvae, within and between populations, was examined using microsatellite analyses at 10 'neutral' loci (presumed to be unaffected by MPTL). Results confirmed that there was little background genetic variation between the populations, global FST <0.038. Polymorphisms present in exons 7 and 8 of Tci-mptl-1 enabled genotyping of individual larvae. A reduction in the number of genotypes was observed in all MPTL-resistant strains compared to the MPTL-susceptible strains that they were derived from, suggesting there was purifying selection at Tci-mptl-1 as a result of MPTL-treatment. The potential link between benzimidazole (BZ)-resistance and MPTL-resistance was examined by screening individual larvae for the presence of three SNPs associated with BZ-resistance in the β-tubulin isotype-1 gene. The majority of larvae were BZ-susceptible homozygotes at positions 167 and 198. Increased heterozygosity at position 200 was observed in the MPTL-resistant strains compared to their respective MPTL-susceptible population. There was no decrease in the occurrence of BZ-resistant genotypes in larvae from each population. These differences, in light of the purifying selection at this locus in all MPTL-resistant isolates, suggests that Tci-mptl-1 confers MPTL-resistance in T. circumcincta, as in Haemonchus contortus, but that different mutations in Tci-mptl-1 can confer resistance in different populations.

The complete mitochondrial genome of the New Zealand parasitic roundworm Teladorsagia circumcincta (Trichostrongyloidea: Haemonchidae) field strain NZ_Teci_NP

The complete mitochondrial genome of the New Zealand parasitic nematode Teladorsagia circumcincta field strain NZ_Teci_NP was sequenced and annotated. The 14,083 bp-long mitogenome contains 12 protein-coding genes (atp8 gene missing), two ribosomal RNAs (rRNAs), and 22 transfer RNAs (tRNAs). Phylogenetic analysis showed that T. circumcincta NZ_Teci_NP forms a monophyletic cluster with the remaining Haemonchidae species and further reinforces the high levels of diversity and gene flow observed among Trichostrongylidae.

Extreme-QTL mapping of monepantel resistance in Haemonchus contortus

Background

Haemonchus contortus, a gastrointestinal nematode parasite of sheep, is mainly controlled by anthelmintics; the occurrence of anthelmintic resistance leads to treatment failures and increases economic burden. Because molecular mechanisms involved in drug resistance can be elucidated by genomic studies, an extreme quantitative trait locus (X-QTL) mapping approach was used to identify co-segregation of the resistance phenotype with genetic markers to detect the genome-wide variants associated with monepantel resistance in H. contortus.

Methods

A cross between H. contortus isolates using parental susceptible (Par-S) males and monepantel resistant (Par-R) females resulted in SR progeny, while reciprocal cross resulted in RS progeny. Pools (n = 30,000) of infective larvae (L3) recovered from Par-R, and from SR and RS populations in the F3 generation, collected both before (unselected group) and 7 days after (selected group) selection with monepantel treatment in sheep hosts, were subjected to genome sequencing (Pool-Seq). Pairwise comparisons of allele frequencies between unselected and selected groups were performed for each population by Fisher’s exact test (FET) and for both populations combined by a Cochran-Mantel-Haenszel (CMH) test.

Results

Mapping rates varied from 80.29 to 81.77% at a 90.4X mean coverage of aligned reads. After correction for multiple testing, significant (P < 0.05) changes in allele frequencies were detected by FET for 6 and 57 single nucleotide polymorphisms (SNPs) in the SR and RS populations, respectively, and by the CMH test for 124 SNPs in both populations. The significant variants located on chromosome 2 generated a selection signal in a genomic region harboring the mptl-1, deg-3 and des-2 genes, previously reported as candidates for monepantel resistance. In addition, three new variants were identified in the mptl-1 gene.

Conclusions

This study expands knowledge on genome-wide molecular events underlying H. contortus resistance to monepantel. The identification of a genome region harboring major genes previously associated with monepantel resistance supports the results of the employed X-QTL approach. In addition, a deletion in exon 11 of the mptl-1 gene should be further investigated as the putative causal mutation leading to monepantel resistance.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3663-9) contains supplementary material, which is available to authorized users.

Perspectives on the utility of moxidectin for the control of parasitic nematodes in the face of developing anthelmintic resistance

Macrocyclic lactone (ML) anthelmintics are the most important class of anthelmintics because of our high dependence on them for the control of nematode parasites and some ectoparasites in livestock, companion animals and in humans. However, resistance to MLs is of increasing concern. Resistance is commonplace throughout the world in nematode parasites of small ruminants and is of increasing concern in horses, cattle, dogs and other animals. It is suspected in Onchocerca volvulus in humans. In most animals, resistance first arose to the avermectins, such as ivermectin (IVM), and subsequently to moxidectin (MOX). Usually when parasite populations are ML-resistant, MOX is more effective than avermectins. MOX may have higher intrinsic potency against some parasites, especially filarial nematodes, than the avermectins. However, it clearly has a significantly different pharmacokinetic profile. It is highly distributed to lipid tissues, less likely to be removed by ABC efflux transporters, is poorly metabolized and has a long half-life. This results in effective concentrations persisting for longer in target hosts. It also has a high safety index. Limited data suggest that anthelmintic resistance may be overcome, at least temporarily, if a high concentration can be maintained at the site of the parasites for a prolonged period of time. Because of the properties of MOX, there are reasonable prospects that strains of parasites that are resistant to avermectins at currently recommended doses will be controlled by MOX if it can be administered at sufficiently high doses and in formulations that enhance its persistence in the host. This review examines the properties of MOX that support this contention and compares them with the properties of other MLs. The case for using MOX to better control ML-resistant parasites is summarised and some outstanding research questions are presented.

Minimal modulation of macrocyclic lactone susceptibility in Caenorhabditis elegans following inhibition of cytochrome P450 monooxygenase activity

Anthelmintic and in particular macrocyclic lactone (ML) resistance is a widespread problem in trichostrongyloid parasitic nematodes, yet mechanisms of ML resistance are still poorly understood. In the absence of target-site changes in resistant parasite field populations, increased drug extrusion and xenobiotic metabolism have been implicated in modification of susceptibility to MLs. In addition to P-glycoproteins, cytochrome P450 monooxygenases (CYPs) were considered to be involved in ML resistance. CYPs are highly divergent in nematodes with about 80 genes in the model organism Caenorhabditis elegans. Using larval development assays in the C. elegans model, piperonyl butoxide (PBO) and a temperature-sensitive variant of the emb-8 cytochrome reductase were used for chemical and genetic ablation of CYP activity. Additionally, a loss-of-function variant of cyp-14A5 was characterized to determine whether increased expression of this CYP in an ivermectin (IVM)-tolerant C. elegans line might be related to the phenotype. In a preliminary experiment with PBO, susceptibility to 5 nM IVM was synergistically increased by PBO. However, effects of genetic ablation of CYP activity on the EC₅₀ values were small (1.5-fold decrease) for IVM and not significant for moxidectin (MOX). However, due to the steep concentration-response-curves, there were again strong differences between the wild-type and the CYP deficient genotype at individual IVM but not MOX concentrations. Although these results suggest small but significant effects on the susceptibility level of C. elegans to IVM, the cyp14A5 gene proposed by a previous study as candidate was ruled out since it was neither IVM/MOX inducible nor did a strain with a loss-of-function allele show increased susceptibility to either drug. In conclusion, the effect of the CYP system on IVM susceptibility in C. elegans is at best low while effects on MOX susceptibility were not detected. The previously suggested candidate cyp14A5 could be excluded to be involved in ML metabolism.

Population genomic and evolutionary modelling analyses reveal a single major QTL for ivermectin drug resistance in the pathogenic nematode, Haemonchus contortus

BACKGROUND

Infections with helminths cause an enormous disease burden in billions of animals and plants worldwide. Large scale use of anthelmintics has driven the evolution of resistance in a number of species that infect livestock and companion animals, and there are growing concerns regarding the reduced efficacy in some human-infective helminths. Understanding the mechanisms by which resistance evolves is the focus of increasing interest; robust genetic analysis of helminths is challenging, and although many candidate genes have been proposed, the genetic basis of resistance remains poorly resolved.

RESULTS

Here, we present a genome-wide analysis of two genetic crosses between ivermectin resistant and sensitive isolates of the parasitic nematode Haemonchus contortus, an economically important gastrointestinal parasite of small ruminants and a model for anthelmintic research. Whole genome sequencing of parental populations, and key stages throughout the crosses, identified extensive genomic diversity that differentiates populations, but after backcrossing and selection, a single genomic quantitative trait locus (QTL) localised on chromosome V was revealed to be associated with ivermectin resistance. This QTL was common between the two geographically and genetically divergent resistant populations and did not include any leading candidate genes, suggestive of a previously uncharacterised mechanism and/or driver of resistance. Despite limited resolution due to low recombination in this region, population genetic analyses and novel evolutionary models supported strong selection at this QTL, driven by at least partial dominance of the resistant allele, and that large resistance-associated haplotype blocks were enriched in response to selection.

CONCLUSIONS

We have described the genetic architecture and mode of ivermectin selection, revealing a major genomic locus associated with ivermectin resistance, the most conclusive evidence to date in any parasitic nematode. This study highlights a novel genome-wide approach to the analysis of a genetic cross in non-model organisms with extreme genetic diversity, and the importance of a high-quality reference genome in interpreting the signals of selection so identified.

The population genetics of parasitic nematodes of wild animals

Parasitic nematodes are highly diverse and common, infecting virtually all animal species, and the importance of their roles in natural ecosystems is increasingly becoming apparent. How genes flow within and among populations of these parasites - their population genetics - has profound implications for the epidemiology of host infection and disease, and for the response of parasite populations to selection pressures. The population genetics of nematode parasites of wild animals may have consequences for host conservation, or influence the risk of zoonotic disease. Host movement has long been recognised as an important determinant of parasitic nematode population genetic structure, and recent research has also highlighted the importance of nematode life histories, environmental conditions, and other aspects of host ecology. Commonly, factors influencing parasitic nematode population genetics have been studied in isolation, such that an integrated view of the drivers of population genetic structure of parasitic nematodes is still lacking. Here, we seek to provide a comprehensive, broad, and integrative picture of these factors in parasitic nematodes of wild animals that will be a useful resource for investigators studying non-model parasitic nematodes in natural ecosystems. Increasingly, new methods of analysing the population genetics of nematodes are becoming available, and we consider the opportunities that these afford in resolving hitherto inaccessible questions of the population genetics of these important animals.

Clonal amplification of Fasciola hepatica in Galba truncatula: within and between isolate variation of triclabendazole-susceptible and -resistant clones

Background

Fasciola hepatica is of worldwide significance, impacting on the health, welfare and productivity of livestock and regarded by WHO as a re-emerging zoonosis. Triclabendazole (TCBZ), the drug of choice for controlling acute fasciolosis in livestock, is also the drug used to treat human infections. However TCBZ-resistance is now considered a major threat to the effective control of F. hepatica. It has yet to be demonstrated whether F. hepatica undergoes a genetic clonal expansion in the snail intermediate host, Galba truncatula, and to what extent amplification of genotypes within the snail facilitates accumulation of drug resistant parasites. Little is known about genotypic and phenotypic variation within and between F. hepatica isolates.

Results

Six clonal isolates of F. hepatica (3× triclabendazole-resistant, TCBZ-R and 3× triclabendazole-susceptible, TCBZ-S) were generated. Snails infected with one miracidium started to shed cercariae 42–56 days post-infection and shed repeatedly up to a maximum of 11 times. A maximum of 884 cercariae were shed by one clonally-infected snail (FhLivS1) at a single time point, with > 3000 clonal metacercariae shed over its lifetime. Following experimental infection all 12 sheep were FEC positive at the time of TCBZ treatment. Sheep infected with one of three putative TCBZ-S clones and treated with TCBZ had no parasites in the liver at post-mortem, whilst sheep each infected with putative TCBZ-R isolates had 35–165 adult fluke at post-mortem, despite TCBZ treatment. All six untreated control animals had between 15–127 parasites. A single multi-locus genotype was reported for every fluke from each of the six clonal isolates. Adult F. hepatica showed considerable variation in weight, ranging from 20–280 mg, with variation in weight evident within and amongst clonal isolates.

Conclusions

A genetic clonal expansion occurs within G. truncatula, highlighting the potential for amplification of drug resistant genotypes of F. hepatica. Variation in the weight of parasites within and between clonal isolates and when comparing isolates that are either susceptible or resistant to TCBZ represent inherent variation in liver fluke and cannot be attributed to their resistance or susceptibility traits.

Discovery of genomic intervals that underlie nematode responses to benzimidazoles

Parasitic nematodes impose a debilitating health and economic burden across much of the world. Nematode resistance to anthelmintic drugs threatens parasite control efforts in both human and veterinary medicine. Despite this threat, the genetic landscape of potential resistance mechanisms to these critical drugs remains largely unexplored. Here, we exploit natural variation in the model nematodes Caenorhabditis elegans and Caenorhabditis briggsae to discover quantitative trait loci (QTL) that control sensitivity to benzimidazoles widely used in human and animal medicine. High-throughput phenotyping of albendazole, fenbendazole, mebendazole, and thiabendazole responses in panels of recombinant lines led to the discovery of over 15 QTL in C. elegans and four QTL in C. briggsae associated with divergent responses to these anthelmintics. Many of these QTL are conserved across benzimidazole derivatives, but others show drug and dose specificity. We used near-isogenic lines to recapitulate and narrow the C. elegans albendazole QTL of largest effect and identified candidate variants correlated with the resistance phenotype. These QTL do not overlap with known benzimidazole target resistance genes from parasitic nematodes and present specific new leads for the discovery of novel mechanisms of nematode benzimidazole resistance. Analyses of orthologous genes reveal conservation of candidate benzimidazole resistance genes in medically important parasitic nematodes. These data provide a basis for extending these approaches to other anthelmintic drug classes and a pathway towards validating new markers for anthelmintic resistance that can be deployed to improve parasite disease control.

The treatment of roundworm (nematode) infections in both humans and animals relies on a small number of anti-parasitic drugs. Resistance to these drugs has appeared in veterinary parasite populations and is a growing concern in human medicine. A better understanding of the genetic basis for parasite drug resistance can be used to help maintain the effectiveness of anti-parasitic drugs and to slow or to prevent the spread of drug resistance in parasite populations. This goal is hampered by the experimental intractability of nematode parasites. Here, we use non-parasitic model nematodes to systematically explore responses to the critical benzimidazole class of anti-parasitic compounds. Using a quantitative genetics approach, we discovered unique genomic intervals that control drug effects, and we identified differences in the effects of these intervals across compounds and doses. We were able to narrow a major-effect genomic region associated with albendazole resistance and to establish that candidate genes discovered in our genetic mappings are largely conserved in important human and animal parasites. This work provides new leads for understanding parasite drug resistance and contributes a powerful template that can be extended to other anti-parasitic drug classes.

A Genome Resequencing-Based Genetic Map Reveals the Recombination Landscape of an Outbred Parasitic Nematode in the Presence of Polyploidy and Polyandry

The parasitic nematode Haemonchus contortus is an economically and clinically important pathogen of small ruminants, and a model system for understanding the mechanisms and evolution of traits such as anthelmintic resistance. Anthelmintic resistance is widespread and is a major threat to the sustainability of livestock agriculture globally; however, little is known about the genome architecture and parameters such as recombination that will ultimately influence the rate at which resistance may evolve and spread. Here, we performed a genetic cross between two divergent strains of H. contortus, and subsequently used whole-genome resequencing of a female worm and her brood to identify the distribution of genome-wide variation that characterizes these strains. Using a novel bioinformatic approach to identify variants that segregate as expected in a pseudotestcross, we characterized linkage groups and estimated genetic distances between markers to generate a chromosome-scale F1 genetic map. We exploited this map to reveal the recombination landscape, the first for any helminth species, demonstrating extensive variation in recombination rate within and between chromosomes. Analyses of these data also revealed the extent of polyandry, whereby at least eight males were found to have contributed to the genetic variation of the progeny analyzed. Triploid offspring were also identified, which we hypothesize are the result of nondisjunction during female meiosis or polyspermy. These results expand our knowledge of the genetics of parasitic helminths and the unusual life-history of H. contortus, and enhance ongoing efforts to understand the genetic basis of resistance to the drugs used to control these worms and for related species that infect livestock and humans throughout the world. This study also demonstrates the feasibility of using whole-genome resequencing data to directly construct a genetic map in a single generation cross from a noninbred nonmodel organism with a complex lifecycle.

P-glycoprotein-9 and macrocyclic lactone resistance status in selected strains of the ovine gastrointestinal nematode, Teladorsagia circumcincta

The Teladorsagia circumcincta P-glycoprotein-9 (Tci-pgp-9) gene has previously been implicated in multiple-anthelmintic resistance in this parasite. Here we further characterise genetic diversity in Tci-pgp-9 and its possible role in ivermectin (IVM) and multi-drug resistance using two UK field isolates of T. circumcincta, one susceptible to anthelmintics (MTci2) and the other resistant to most available anthelmintics including IVM (MTci5). A comparison of full-length Tci-pgp-9 cDNA transcripts from the MTci2 and MTci5 isolates (∼3.8 kb in both cases) indicated that they shared 95.6% and 99.5% identity at the nucleotide and amino acid levels, respectively. Nine non-synonymous SNPs were found in the MTci5 sequences relative to their MTci2 counterparts. Twelve genomic sequence variants of the first internucleotide binding domain of Tci-pgp-9 were identified and up to 10 of these were present in some individual worms, strongly supporting previous evidence that amplification of this gene has occurred in T. circumcincta. On average, fewer distinct sequence variants of Tci-pgp-9 were present in individual worms of the MTci5 isolate than in those of the MTci2 isolate. A further reduction in the number of sequence variants was observed in individuals derived from an IVM-treated sub-population of MTci5. These findings suggest that Tci-pgp-9 was under purifying selection in the face of IVM treatment in T. circumcincta, with some sequence variants being selected against.

•

Tci-pgp-9 cDNA sequences from the MTci2 and MTci5 isolates shared 95.6% identity.

•

Comparison of Tci-pgp-9 sequences from MTci5 and MTci2 showed 9 non-synonymous SNPs.

•

Multiple Tci-pgp-9-IBDA sequence variants occurred in the majority of individual larvae.

•

Tci-pgp-9 appears to be under purifying selection after anthelmintics are applied.

Increased Expression of a MicroRNA Correlates with Anthelmintic Resistance in Parasitic Nematodes

Resistance to anthelmintic drugs is a major problem in the global fight against parasitic nematodes infecting humans and animals. While previous studies have identified mutations in drug target genes in resistant parasites, changes in the expression levels of both targets and transporters have also been reported. The mechanisms underlying these changes in gene expression are unresolved. Here, we take a novel approach to this problem by investigating the role of small regulatory RNAs in drug resistant strains of the important parasite Haemonchus contortus. microRNAs (miRNAs) are small (22 nt) non-coding RNAs that regulate gene expression by binding predominantly to the 3′ UTR of mRNAs. Changes in miRNA expression have been implicated in drug resistance in a variety of tumor cells. In this study, we focused on two geographically distinct ivermectin resistant strains of H. contortus and two lines generated by multiple rounds of backcrossing between susceptible and resistant parents, with ivermectin selection. All four resistant strains showed significantly increased expression of a single miRNA, hco-miR-9551, compared to the susceptible strain. This same miRNA is also upregulated in a multi-drug-resistant strain of the related nematode Teladorsagia circumcincta. hco-miR-9551 is enriched in female worms, is likely to be located on the X chromosome and is restricted to clade V parasitic nematodes. Genes containing predicted binding sites for hco-miR-9551 were identified computationally and refined based on differential expression in a transcriptomic dataset prepared from the same drug resistant and susceptible strains. This analysis identified three putative target mRNAs, one of which, a CHAC domain containing protein, is located in a region of the H. contortus genome introgressed from the resistant parent. hco-miR-9551 was shown to interact with the 3′ UTR of this gene by dual luciferase assay. This study is the first to suggest a role for miRNAs and the genes they regulate in drug resistant parasitic nematodes. miR-9551 also has potential as a biomarker of resistance in different nematode species.