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 and that European-type 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) from 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 Em and E. granulosus s.l. genotypes G8/G10, low intraspecific diversity of Em, and a higher abundance of the European-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.

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.

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.

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.

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

BACKGROUND: 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. CONCLUSIONS: 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.