Evidence from two independent backcross experiments supports genetic linkage of microsatellite Hcms8a20, but not other candidate loci, to a major ivermectin resistance locus in Haemonchus contortus

Haemonchus contortus is the leading parasitic nematode species used to study anthelmintic drug resistance. A variety of candidate loci have been implicated as being associated with ivermectin resistance in this parasite but definitive evidence of their importance is still lacking. We have previously performed two independent serial backcross experiments to introgress ivermectin resistance loci from two H. contortus ivermectin-resistant strains – MHco4(WRS) and MHco10(CAVR) – into the genetic background of the ivermectin-susceptible genome reference strain MHco3(ISE). We have interrogated a number of candidate ivermectin resistance loci in the resulting backcross populations and assessed the evidence for their genetic linkage to an ivermectin resistance locus. These include the microsatellite marker Hcms8a20 and six candidate genes Hco-glc-5, Hco-avr-14, Hco-lgc-37 (previously designated Hco-hg-1), Hco-pgp-9 (previously designated Hco-pgp-1), Hco-pgp-2 and Hco-dyf-7. We have sampled the haplotype diversity of amplicon markers within, or adjacent to, each of these loci in the parental strains and fourth generation backcross populations to assess the evidence for haplotype introgression from the resistant parental strain into the genomic background of the susceptible parental strain in each backcross. The microsatellite Hcms8a20 locus showed strong evidence of such introgression in both independent backcrosses, suggesting it is linked to an important ivermectin resistance mutation in both the MHco4(WRS) and MHco10(CAVR) strains. In contrast, Hco-glc-5, Hco-avr-14, Hco-pgp-9 and Hco-dyf-7 showed no evidence of introgression in either backcross. Hco-lgc-37 and Hco-pgp-2 showed only weak evidence of introgression in the MHco3/4 backcross but not in the MHco3/10 backcross. Overall, these results suggest that microsatellite marker Hcms8a20, but not the other candidate genes tested, is linked to a major ivermectin resistance locus in the MHco4(WRS) and MHco10(CAVR) strains. This work also emphasises the need for genome-wide approaches to identify mutations responsible for the ivermectin resistance in this parasite

The cytochrome P450 family in the parasitic nematode <i>Haemonchus contortus</i>

&lt;i&gt;Haemonchus contortus&lt;/i&gt;, a highly pathogenic and economically important parasitic nematode of sheep, is particularly adept at developing resistance to the anthelmintic drugs used in its treatment and control. The basis of anthelmintic resistance is poorly understood for many commonly used drugs with most research being focused on mechanisms involving drug targets or drug efflux. Altered or increased drug metabolism is a possible mechanism that has yet to receive much attention despite the clear role of xenobiotic metabolism in pesticide resistance in insects. The cytochrome P450s (CYPs) are a large family of drug-metabolising enzymes present in almost all living organisms, but for many years thought to be absent from parasitic nematodes. In this paper, we describe the CYP sequences encoded in the &lt;i&gt;H. Contortus&lt;/i&gt; genome and compare their expression in different parasite life-stages, sexes and tissues. We developed a novel real-time PCR approach based on partially assembled CYP sequences “tags” and confirmed findings in the subsequent draft genome with RNA-seq. Constitutive expression was highest in larval stages for the majority of CYPs, although higher expression was detected in the adult male or female for a small subset of genes. Many CYPs were expressed in the worm intestine. A number of &lt;i&gt;H. Contortus&lt;/i&gt; genes share high identity with &lt;i&gt;Caenorhabditis elegans&lt;/i&gt; CYPs and the similarity in their expression profiles supports their classification as putative orthologues. Notably, &lt;i&gt;H. Contortus&lt;/i&gt; appears to lack the dramatic CYP subfamily expansions seen in &lt;i&gt;C. elegans&lt;/i&gt; and other species, which are typical of CYPs with exogenous roles. However, a small group of &lt;i&gt;H. Contortus&lt;/i&gt; genes cluster with the &lt;i&gt;C. elegans&lt;/i&gt; CYP34 and CYP35 subfamilies and may represent candidate xenobiotic metabolising genes in the parasite

Deep amplicon sequencing highlights low intra-host genetic variability of Echinococcus multilocularis and high prevalence of the European-type haplotypes in coyotes and red foxes in Alberta, Canada

Echinococcus multilocularis (Em) is a zoonotic parasite considered a global emergent pathogen. Recent findings indicate that the parasite is expanding its range in North America (NA) and that European-type (EU) haplotypes are circulating in western Canada. However, genetic analyses are usually conducted only on a few parasites out of thousands of individuals within each definitive host, likely underestimating the prevalence of less common haplotypes. Moreover, mixed infections with several mtDNA haplotypes in the same host have been reported, but their relative abundance within the host was never estimated. We aimed to 1) estimate the frequency of co-infections of different Em haplotypes in coyotes (Canis latrans) and red foxes (Vulpes vulpes) of western Canada and their relative abundance within the definitive hosts, 2) detect less prevalent haplotypes by sampling a larger proportion of the parasite subpopulation per host, and 3) investigate differences in the distribution of Em haplotypes in these main definitive hosts; foxes and coyotes. We extracted DNA from ~10% of the worm subpopulation per host (20 foxes and 47 coyotes) and used deep amplicon sequencing (NGS technology) on four loci, targeting the most polymorphic regions from the mitochondrial genes cox1 (814 bp), nad1 (344 bp), and cob (387 bp). We detected the presence of mixed infections with multiple Em haplotypes and with different Echinococcus species including E. granulosus s.l. genotypes G8/G10, low intraspecific diversity of Em, and a higher abundance of the EU-type haplotypes in both hosts. Our results suggest a population expansion of the European over the North American strain in Alberta and a limited distribution of some European-type haplotypes. Our findings indicate that deep amplicon sequencing represents a valuable tool to characterize Em in multiple hosts, to assess the current distribution and possible origins of the European strain in North America. The potential use of next-generation sequencing technologies is particularly important to understand the patterns of geographic expansion of this parasite

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

The confounding effects of high genetic diversity on the determination and interpretation of differential gene expression analysis in the parasitic nematode Haemonchus contortus

Differential expression analysis between parasitic nematode strains is commonly used to implicate candidate genes in anthelmintic resistance or other biological functions. We have tested the hypothesis that the high genetic diversity of an organism such as Haemonchus contortus could complicate such analyses. First, we investigated the extent to which sequence polymorphism affects the reliability of differential expression analysis between the genetically divergent H. contortus strains MHco3(ISE), MHco4(WRS) and MHco10(CAVR). Using triplicates of 20 adult female worms from each population isolated under parallel experimental conditions, we found that high rates of sequence polymorphism in RNAseq reads were associated with lower efficiency read mapping to gene models under default TopHat2 parameters, leading to biased estimates of inter-strain differential expression. We then showed it is possible to largely compensate for this bias by optimising the read mapping single nucleotide polymorphism (SNP) allowance and filtering out genes with particularly high single nucleotide polymorphism rates. Once the sequence polymorphism biases were removed, we then assessed the genuine transcriptional diversity between the strains, finding ≥824 differentially expressed genes across all three pairwise strain comparisons. This high level of inter-strain transcriptional diversity not only suggests substantive inter-strain phenotypic variation but also highlights the difficulty in reliably associating differential expression of specific genes with phenotypic differences. To provide a practical example, we analysed two gene families of potential relevance to ivermectin drug resistance; the ABC transporters and the ligand-gated ion channels (LGICs). Over half of genes identified as differentially expressed using default TopHat2 parameters were shown to be an artifact of sequence polymorphism differences. This work illustrates the need to account for sequence polymorphism in differential expression analysis. It also demonstrates that a large number of genuine transcriptional differences can occur between H. contortus strains and these must be considered before associating the differential expression of specific genes with phenotypic differences between strains

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.Wellcome, BBSR

Resurrection and redescription of Varestrongylus alces (Nematoda; Protostrongylidae), a lungworm of the Eurasian moose (Alces alces), with report on associated pathology

Varestrongylus alces, a lungworm in Eurasian moose from Europe has been considered a junior synonym of Varestrongylus capreoli, in European roe deer, due to a poorly detailed morphological description and the absence of a type-series. Methods Specimens used in the redescription were collected from lesions in the lungs of Eurasian moose, from Vestby, Norway. Specimens were described based on comparative morphology and integrated approaches. Molecular identification was based on PCR, cloning and sequencing of the ITS-2 region of the nuclear ribosomal DNA. Phylogenetic analysis compared V. alces ITS-2 sequences to these of other Varestrongylus species and other protostrongylids. Results Varestrongylus alces is resurrected for protostrongylid nematodes of Eurasian moose from Europe. Varestrongylus alces causes firm nodular lesions that are clearly differentiated from the adjacent lung tissue. Histologically, lesions are restricted to the parenchyma with adult, egg and larval parasites surrounded by multinucleated giant cells, macrophages, eosinophilic granulocytes, lymphocytes. The species is valid and distinct from others referred to Varestrongylus, and should be separated from V. capreoli. Morphologically, V. alces can be distinguished from other species by characters in the males that include a distally bifurcated gubernaculum, arched denticulate crura, spicules that are equal in length and relatively short, and a dorsal ray that is elongate and bifurcated. Females have a well-developed provagina, and are very similar to those of V. capreoli. Morphometrics of first-stage larvae largely overlap with those of other Varestrongylus. Sequences of the ITS-2 region strongly support mutual independence of V. alces, V. cf. capreoli, and the yet undescribed species of Varestrongylus from North American ungulates. These three taxa form a well-supported crown-clade as the putative sister of V. alpenae. The association of V. alces and Alces or its ancestors is discussed in light of host and parasite phylogeny and host historical biogeography. Varestrongylus alces is a valid species, and should be considered distinct from V. capreoli. Phylogenetic relationships among Varestrongylus spp. from Eurasia and North America are complex and consistent with faunal assembly involving recurrent events of geographic expansion, host switching and subsequent speciation. Cervidae, Cryptic species, Historical biogeography, ITS-2, Metastrongyloidea, Parasite biodiversity, Varestrongylinae, Varestrongylus capreoli, Verminous pneumoniapublishedVersio

Investigating the genetic basis of ivermectin resistance in Haemonchus contortus

Parasitic nematodes have a major impact on human and animal health and their control is threatened by the emergence of resistance to the anthelmintic drugs on which control depends. Although resistance to ivermectin is widespread in many livestock parasites, the genetic mechanisms remain elusive. In this thesis, a number of approaches are undertaken to investigate the genetics of ivermectin resistance in Haemonchus contortus, an important parasitic nematode research model. A transcriptomic comparison of the susceptible MHco3(ISE) and two ivermectin resistant strains, MHco4(WRS) and MHco10(CAVR), was undertaken to investigate if differences in gene expression revealed insights into ivermectin resistance mechanisms. Although this approach did not reveal clear ivermectin resistance gene candidates, the results have some important implications for RNAseq analysis in genetically diverse organisms. First, extremely high levels of genetic diversity in H. contortus had a major impact on RNAseq analysis and needed to be accounted for reliable identification of differentially expressed genes. Second, there were a remarkably large number of genes differentially expressed between the three H. contortus strains. In a novel genetic crossing approach, two independent serially backcrossed H. contortus populations were utilized in which ivermectin resistance loci had been introgressed into the MHco3(ISE) genome. Markers for six leading candidate ivermectin resistance genes from the literature were examined and none were found to show evidence of introgression in the backcross populations. A seventh marker, Hcms8a20 showed clear evidence of introgression and was later confirmed by collaborators using a whole genome sequencing approach to show its location in a 11.2 Mb introgressed region. A set of 25 in vivo or in vitro phenotyped ivermectin resistant H. contortus field populations were then used to investigate evidence of selection on the Hcms8a20 marker and four other loci distributed across the 11.2 Mb introgressed region. A deep amplicon sequencing approach revealed that, of the loci tested, only Hcms8a20 showed consistent evidence of purifying selection in the ivermectin resistant field populations. Overall, the results provide strong evidence that a major ivermectin resistance locus is located within a few megabases of the Hcms8a20 locus in many different ivermectin resistant H. contortus field populations

Microsattilite genotypes for all individuals at all 16 loci

This file contains the microsattilite data that was attained for this study. It is in genepop format. This file contains all loci that were genotyped including those that were removed from the analysis. The 16 numbers at the top of the file represent loci names. Each entry under a "POP" header represents an individual of that population, where A = Alaska, B = Brannen Lake, E = East Coast, M = Mayliewan (our Alberta population), and N = Norway. Each of the 16 sets of 4 digits represent diploid genotypes at each of the 16 loci for that individual. For example, 0414 would represent a heterozygous call for allele 4 and allele 14 for that individual for whatever locus it is under. The first column of 4 digits represents locus 51 and the last represents locus 222 (the rest are ordered accordingly). 00s represent missing data