The impact of FASTQ and alignment read order on structural variant calling from long-read sequencing data.

Background: Structural variant (SV) calling from DNA sequencing data has been challenging due to several factors, including the ambiguity of short-read alignments, multiple complex SVs in the same genomic region, and the lack of "truth" datasets for benchmarking. Additionally, caller choice, parameter settings, and alignment method are known to affect SV calling. However, the impact of FASTQ read order on SV calling has not been explored for long-read data. Results: Here, we used PacBio DNA sequencing data from 15 Caenorhabditis elegans strains and four Arabidopsis thaliana ecotypes to evaluate the sensitivity of different SV callers on FASTQ read order. Comparisons of variant call format files generated from the original and permutated FASTQ files demonstrated that the order of input data affected the SVs predicted by each caller. In particular, pbsv was highly sensitive to the order of the input data, especially at the highest depths where over 70% of the SV calls generated from pairs of differently ordered FASTQ files were in disagreement. These demonstrate that read order sensitivity is a complex, multifactorial process, as the differences observed both within and between species varied considerably according to the specific combination of aligner, SV caller, and sequencing depth. In addition to the SV callers being sensitive to the input data order, the SAMtools alignment sorting algorithm was identified as a source of variability following read order randomization. Conclusion: The results of this study highlight the sensitivity of SV calling on the order of reads encoded in FASTQ files, which has not been recognized in long-read approaches. These findings have implications for the replication of SV studies and the development of consistent SV calling protocols. Our study suggests that researchers should pay attention to the input order sensitivity of read alignment sorting methods when analyzing long-read sequencing data for SV calling, as mitigating a source of variability could facilitate future replication work. These results also raise important questions surrounding the relationship between SV caller read order sensitivity and tool performance. Therefore, tool developers should also consider input order sensitivity as a potential source of variability during the development and benchmarking of new and improved methods for SV calling.

IL-33-binding HpARI family homologues with divergent effects in suppressing or enhancing type 2 immune responses.

HpARI is an immunomodulatory protein secreted by the intestinal nematode Heligmosomoides polygyrus bakeri, which binds and blocks IL-33. Here, we find that the H. polygyrus bakeri genome contains three HpARI family members and that these have different effects on IL-33-dependent responses in vitro and in vivo, with HpARI1+2 suppressing and HpARI3 amplifying these responses. All HpARIs have sub-nanomolar affinity for mouse IL-33; however, HpARI3 does not block IL-33-ST2 interactions. Instead, HpARI3 stabilizes IL-33, increasing the half-life of the cytokine and amplifying responses to it in vivo. Together, these data show that H. polygyrus bakeri secretes a family of HpARI proteins with both overlapping and distinct functions, comprising a complex immunomodulatory arsenal of host-targeted proteins.

Identification of potential molecular mimicry in pathogen-host interactions.

Pathogens have evolved sophisticated strategies to manipulate host signaling pathways, including the phenomenon of molecular mimicry, where pathogen-derived biomolecules imitate host biomolecules. In this study, we resurrected, updated, and optimized a sequence-based bioinformatics pipeline to identify potential molecular mimicry candidates between humans and 32 pathogenic species whose proteomes' 3D structure predictions were available at the start of this study. We observed considerable variation in the number of mimicry candidates across pathogenic species, with pathogenic bacteria exhibiting fewer candidates compared to fungi and protozoans. Further analysis revealed that the candidate mimicry regions were enriched in solvent-accessible regions, highlighting their potential functional relevance. We identified a total of 1,878 mimicked regions in 1,439 human proteins, and clustering analysis indicated diverse target proteins across pathogen species. The human proteins containing mimicked regions revealed significant associations between these proteins and various biological processes, with an emphasis on host extracellular matrix organization and cytoskeletal processes. However, immune-related proteins were underrepresented as targets of mimicry. Our findings provide insights into the broad range of host-pathogen interactions mediated by molecular mimicry and highlight potential targets for further investigation. This comprehensive analysis contributes to our understanding of the complex mechanisms employed by pathogens to subvert host defenses and we provide a resource to assist researchers in the development of novel therapeutic strategies.

Transcriptional patterns of sexual dimorphism and in host developmental programs in the model parasitic nematode Heligmosomoides bakeri.

BACKGROUND: Heligmosomoides bakeri (often mistaken for Heligmosomoides polygyrus) is a promising model for parasitic nematodes with the key advantage of being amenable to study and manipulation within a controlled laboratory environment. While draft genome sequences are available for this worm, which allow for comparative genomic analyses between nematodes, there is a notable lack of information on its gene expression. METHODS: We generated biologically replicated RNA-seq datasets from samples taken throughout the parasitic life of H. bakeri. RNA from tissue-dwelling and lumen-dwelling worms, collected under a dissection microscope, was sequenced on an Illumina platform. RESULTS: We find extensive transcriptional sexual dimorphism throughout the fourth larval and adult stages of this parasite and identify alternative splicing, glycosylation, and ubiquitination as particularly important processes for establishing and/or maintaining sex-specific gene expression in this species. We find sex-linked differences in transcription related to aging and oxidative and osmotic stress responses. We observe a starvation-like signature among transcripts whose expression is consistently upregulated in males, which may reflect a higher energy expenditure by male worms. We detect evidence of increased importance for anaerobic respiration among the adult worms, which coincides with the parasite's migration into the physiologically hypoxic environment of the intestinal lumen. Furthermore, we hypothesize that oxygen concentration may be an important driver of the worms encysting in the intestinal mucosa as larvae, which not only fully exposes the worms to their host's immune system but also shapes many of the interactions between the host and parasite. We find stage- and sex-specific variation in the expression of immunomodulatory genes and in anthelmintic targets. CONCLUSIONS: We examine how different the male and female worms are at the molecular level and describe major developmental events that occur in the worm, which extend our understanding of the interactions between this parasite and its host. In addition to generating new hypotheses for follow-up experiments into the worm's behavior, physiology, and metabolism, our datasets enable future more in-depth comparisons between nematodes to better define the utility of H. bakeri as a model for parasitic nematodes in general.

Trickle infection with Heligmosomoides polygyrus results in decreased worm burdens but increased intestinal inflammation and scarring.

Introduction: Intestinal roundworms cause chronic debilitating disease in animals, including humans. Traditional experimental models of these types of infection use a large single-dose infection. However, in natural settings, hosts are exposed to parasites on a regular basis and when mice are exposed to frequent, smaller doses of Heligmosomoides polygyrus, the parasites are cleared more quickly. Whether this more effective host response has any negative consequences for the host is not known. Results: Using a trickle model of infection, we found that worm clearance was associated with known resistance-related host responses: increased granuloma and tuft cell numbers, increased levels of granuloma IgG and decreased intestinal transit time, as well as higher serum IgE levels. However, we found that the improved worm clearance was also associated with an inflammatory phenotype in and around the granuloma, increased smooth muscle hypertrophy/hyperplasia, and elevated levels of Adamts gene expression. Discussion: To our knowledge, we are the first to identify the involvement of this protein family of matrix metalloproteinases (MMPs) in host responses to helminth infections. Our results highlight the delicate balance between parasite clearance and host tissue damage, which both contribute to host pathology. When continually exposed to parasitic worms, improved clearance comes at a cost.

Different structural variant prediction tools yield considerably different results in Caenorhabditis elegans.

The accurate characterization of structural variation is crucial for our understanding of how large chromosomal alterations affect phenotypic differences and contribute to genome evolution. Whole-genome sequencing is a popular approach for identifying structural variants, but the accuracy of popular tools remains unclear due to the limitations of existing benchmarks. Moreover, the performance of these tools for predicting variants in non-human genomes is less certain, as most tools were developed and benchmarked using data from the human genome. To evaluate the use of long-read data for the validation of short-read structural variant calls, the agreement between predictions from a short-read ensemble learning method and long-read tools were compared using real and simulated data from Caenorhabditis elegans. The results obtained from simulated data indicate that the best performing tool is contingent on the type and size of the variant, as well as the sequencing depth of coverage. These results also highlight the need for reference datasets generated from real data that can be used as 'ground truth' in benchmarks.

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.

Genotyping Canadian Cyclospora cayetanensis Isolates to Supplement Cyclosporiasis Outbreak Investigations.

Cyclospora cayetanensis is an emerging foodborne parasite that causes cyclosporiasis, an enteric disease of humans. Domestically acquired outbreaks have been reported in Canada every spring or summer since 2013. To date, investigations into the potential sources of infection have relied solely on epidemiological data. To supplement the epidemiological data with genetic information, we genotyped 169 Canadian cyclosporiasis cases from stool specimens collected from 2010 to 2021 using an existing eight-marker targeted amplicon deep (TADS) scheme specific to C. cayetanensis as previously described by the US Centers for Disease Control and Prevention (CDC). This is the first study to genotype Canadian Cyclospora cayetanensis isolates, and it focuses on evaluating the genotyping performance and genetic clustering. Genotyping information was successfully collected with at least part of one of the markers in the TADS assay for 97.9% of specimens, and 81.1% of cyclosporiasis cases met the minimum requirements to genetically cluster into 20 groups. The performance of the scheme suggests that examining cyclosporiasis cases genetically will be a valuable tool for supplementing epidemiological outbreak investigations and to minimize further infections. Further research is required to expand the number of discriminatory markers to improve genetic clustering.

Draft Hybrid Genome Assembly of a Canadian Cyclospora cayetanensis Isolate.

The apicomplexan parasite Cyclospora cayetanensis causes foodborne gastrointestinal disease in humans. Here, we report the first hybrid assembly for C. cayetanensis, which uses both Illumina MiSeq and Oxford Nanopore Technologies MinION platforms to generate genomic sequence data. The final genome assembly consists of 44,586,677 bases represented in 313 contigs.

Characterization of the ß-tubulin gene family in Ascaris lumbricoides and Ascaris suum and its implication for the molecular detection of benzimidazole resistance.

BACKGROUND: The treatment coverage of control programs providing benzimidazole (BZ) drugs to eliminate the morbidity caused by soil-transmitted helminths (STHs) is unprecedently high. This high drug pressure may result in the development of BZ resistance in STHs and so there is an urgent need for surveillance systems detecting molecular markers associated with BZ resistance. A critical prerequisite to develop such systems is an understanding of the gene family encoding ß-tubulin proteins, the principal targets of BZ drugs. METHODOLOGY AND PRINCIPAL FINDINGS: First, the ß-tubulin gene families of Ascaris lumbricoides and Ascaris suum were characterized through the analysis of published genomes. Second, RNA-seq and RT-PCR analyses on cDNA were applied to determine the transcription profiles of the different gene family members. The results revealed that Ascaris species have at least seven different ß-tubulin genes of which two are highly expressed during the entire lifecycle. Third, deep amplicon sequencing was performed on these two genes in more than 200 adult A. lumbricoides (Ethiopia and Tanzania) and A. suum (Belgium) worms, to investigate the intra- and inter-species genetic diversity and the presence of single nucleotide polymorphisms (SNPs) that are associated with BZ resistance in other helminth species; F167Y (TTC>TAC or TTT>TAT), E198A (GAA>GCA or GAG>GCG), E198L (GAA>TTA) and F200Y (TTC>TAC or TTT>TAT). These particular SNPs were absent in the two investigated genes in all three Ascaris populations. SIGNIFICANCE: This study demonstrated the presence of at least seven ß-tubulin genes in Ascaris worms. A new nomenclature was proposed and prioritization of genes for future BZ resistance research was discussed. This is the first comprehensive description of the ß-tubulin gene family in Ascaris and provides a framework to investigate the prevalence and potential role of ß-tubulin sequence polymorphisms in BZ resistance in a more systematic manner than previously possible.

Small Intestinal Levels of the Branched Short-Chain Fatty Acid Isovalerate Are Elevated during Infection with Heligmosomoides polygyrus and Can Promote Helminth Fecundity.

Heligmosomoides polygyrus is a helminth which naturally infects mice and is widely used as a laboratory model of chronic small intestinal helminth infection. While it is known that infection with H. polygyrus alters the composition of the host's bacterial microbiota, the functional implications of this alteration are unclear. We investigated the impact of H. polygyrus infection on short-chain fatty acid (SCFA) levels in the mouse intestine and sera. We found that helminth infection resulted in significantly upregulated levels of the branched SCFA isovaleric acid, exclusively in the proximal small intestine, which is the site of H. polygyrus colonization. We next set out to test the hypothesis that elevating local levels of isovaleric acid was a strategy used by H. polygyrus to promote its own fitness within the mammalian host. To test this, we supplemented the drinking water of mice with isovalerate during H. polygyrus infection and examined whether this affected helminth fecundity or chronicity. We did not find that isovaleric acid supplementation affected helminth chronicity; however, we found that it did promote helminth fecundity, as measured by helminth egg output in the feces of mice. Through antibiotic treatment of helminth-infected mice, we found that the bacterial microbiota was required in order to support elevated levels of isovaleric acid in the proximal small intestine during helminth infection. Overall, our data reveal that during H. polygyrus infection there is a microbiota-dependent localized increase in the production of isovaleric acid in the proximal small intestine and that this supports helminth fecundity in the murine host.

Benchmarking hybrid assemblies of Giardia and prediction of widespread intra-isolate structural variation.

BACKGROUND: Currently available short read genome assemblies of the tetraploid protozoan parasite Giardia intestinalis are highly fragmented, highlighting the need for improved genome assemblies at a reasonable cost. Long nanopore reads are well suited to resolve repetitive genomic regions resulting in better quality assemblies of eukaryotic genomes. Subsequent addition of highly accurate short reads to long-read assemblies further improves assembly quality. Using this hybrid approach, we assembled genomes for three Giardia isolates, two with published assemblies and one novel, to evaluate the improvement in genome quality gained from long reads. We then used the long reads to predict structural variants to examine this previously unexplored source of genetic variation in Giardia. METHODS: With MinION reads for each isolate, we assembled genomes using several assemblers specializing in long reads. Assembly metrics, gene finding, and whole genome alignments to the reference genomes enabled direct comparison to evaluate the performance of the nanopore reads. Further improvements from adding Illumina reads to the long-read assemblies were evaluated using gene finding. Structural variants were predicted from alignments of the long reads to the best hybrid genome for each isolate and enrichment of key genes was analyzed using random genome sampling and calculation of percentiles to find thresholds of significance. RESULTS: Our hybrid assembly method generated reference quality genomes for each isolate. Consistent with previous findings based on SNPs, examination of heterozygosity using the structural variants found that Giardia BGS was considerably more heterozygous than the other isolates that are from Assemblage A. Further, each isolate was shown to contain structural variant regions enriched for variant-specific surface proteins, a key class of virulence factor in Giardia. CONCLUSIONS: The ability to generate reference quality genomes from a single MinION run and a multiplexed MiSeq run enables future large-scale comparative genomic studies within the genus Giardia. Further, prediction of structural variants from long reads allows for more in-depth analyses of major sources of genetic variation within and between Giardia isolates that could have effects on both pathogenicity and host range.

Horizontal transfer of a retrotransposon between parasitic nematodes and the common shrew.

BACKGROUND: As the genomes of more metazoan species are sequenced, reports of horizontal transposon transfers (HTT) have increased. Our understanding of the mechanisms of such events is at an early stage. The close physical relationship between a parasite and its host could facilitate horizontal transfer. To date, two studies have identified horizontal transfer of RTEs, a class of retrotransposable elements, involving parasites: ticks might act as vector for BovB between ruminants and squamates, and AviRTE was transferred between birds and parasitic nematodes. RESULTS: We searched for RTEs shared between nematode and mammalian genomes. Given their physical proximity, it was necessary to detect and remove sequence contamination from the genome datasets, which would otherwise distort the signal of horizontal transfer. We developed an approach that is based on reads instead of genomic sequences to reliably detect contamination. From comparison of 43 RTEs across 197 genomes, we identified a single putative case of horizontal transfer: we detected RTE1_Sar from Sorex araneus, the common shrew, in parasitic nematodes. From the taxonomic distribution and evolutionary analysis, we show that RTE1_Sar was horizontally transferred. CONCLUSION: We identified a new horizontal RTE transfer in host-parasite interactions, which suggests that it is not uncommon. Further, we present and provide the workflow a read-based method to distinguish between contamination and horizontal transfer.

MIPhy: identify and quantify rapidly evolving members of large gene families.

After transitioning to a new environment, species often exhibit rapid phenotypic innovation. One of the fastest mechanisms for this is duplication followed by specialization of existing genes. When this happens to a member of a gene family, it tends to leave a detectable phylogenetic signature of lineage-specific expansions and contractions. These can be identified by analyzing the gene family across several species and identifying patterns of gene duplication and loss that do not correlate with the known relationships between those species. This signature, termed phylogenetic instability, has been previously linked to adaptations that change the way an organism samples and responds to its environment; conversely, low phylogenetic instability has been previously linked to proteins with endogenous functions. With the increase in genome-level data, there is a need to identify and quantify phylogenetic instability. Here, we present Minimizing Instability in Phylogenetics (MIPhy), a tool that solves this problem by quantifying the incongruence of a gene's evolutionary history. The motivation behind MIPhy was to produce a tool to aid in interpreting phylogenetic trees. It can predict which members of a gene family are under adaptive evolution, working only from a gene tree and the relationship between the species under consideration. While it does not conduct any estimation of positive selection-which is the typical indication of adaptive evolution-the results tend to agree. We demonstrate the usefulness of MIPhy by accurately predicting which members of the mammalian cytochrome P450 gene superfamily metabolize xenobiotics and which metabolize endogenous compounds. Our predictions correlate very well with known substrate specificities of the human enzymes. We also analyze the Caenorhabditis collagen gene family and use MIPhy to predict genes that produce an observable phenotype when knocked down in C. elegans, and show that our predictions correlate well with existing knowledge. The software can be downloaded and installed from https://github.com/dave-the-scientist/miphy and is also available as an online web tool at http://www.miphy.wasmuthlab.org.

Improved enzyme annotation with EC-specific cutoffs using DETECT v2.

Summary: We present DETECT v2-an enzyme annotation tool which considers the effect of sequence diversity when assigning enzymatic function [as an Enzyme Commission (EC) number] to a protein sequence. In addition to capturing more enzyme classes than the previous version, we now provide EC-specific cutoffs that greatly increase precision and recall of assignments and show its performance in the context of pathways. Availability and implementation: https://github.com/ParkinsonLab/DETECT-v2. Supplementary information: Supplementary data are available at Bioinformatics online.

Expanding the view on the evolution of the nematode dauer signalling pathways: refinement through gene gain and pathway co-option.

BACKGROUND: Signalling pathways underlie development, behaviour and pathology. To understand patterns in the evolution of signalling pathways, we undertook a comprehensive investigation of the pathways that control the switch between growth and developmentally quiescent dauer in 24 species of nematodes spanning the phylum. RESULTS: Our analysis of 47 genes across these species indicates that the pathways and their interactions are not conserved throughout the Nematoda. For example, the TGF-ß pathway was co-opted into dauer control relatively late in a lineage that led to the model species Caenorhabditis elegans. We show molecular adaptations described in C. elegans that are restricted to its genus or even just to the species. Similarly, our analyses both identify species where particular genes have been lost and situations where apparently incorrect orthologues have been identified. CONCLUSIONS: Our analysis also highlights the difficulties of working with genome sequences from non-model species as reliance on the published gene models would have significantly restricted our understanding of how signalling pathways evolve. Our approach therefore offers a robust standard operating procedure for genomic comparisons.

Exploring the Gastrointestinal "Nemabiome": Deep Amplicon Sequencing to Quantify the Species Composition of Parasitic Nematode Communities.

Parasitic helminth infections have a considerable impact on global human health as well as animal welfare and production. Although co-infection with multiple parasite species within a host is common, there is a dearth of tools with which to study the composition of these complex parasite communities. Helminth species vary in their pathogenicity, epidemiology and drug sensitivity and the interactions that occur between co-infecting species and their hosts are poorly understood. We describe the first application of deep amplicon sequencing to study parasitic nematode communities as well as introduce the concept of the gastro-intestinal "nemabiome". The approach is analogous to 16S rDNA deep sequencing used to explore microbial communities, but utilizes the nematode ITS-2 rDNA locus instead. Gastro-intestinal parasites of cattle were used to develop the concept, as this host has many well-defined gastro-intestinal nematode species that commonly occur as complex co-infections. Further, the availability of pure mono-parasite populations from experimentally infected cattle allowed us to prepare mock parasite communities to determine, and correct for, species representation biases in the sequence data. We demonstrate that, once these biases have been corrected, accurate relative quantitation of gastro-intestinal parasitic nematode communities in cattle fecal samples can be achieved. We have validated the accuracy of the method applied to field-samples by comparing the results of detailed morphological examination of L3 larvae populations with those of the sequencing assay. The results illustrate the insights that can be gained into the species composition of parasite communities, using grazing cattle in the mid-west USA as an example. However, both the technical approach and the concept of the 'nemabiome' have a wide range of potential applications in human and veterinary medicine. These include investigations of host-parasite and parasite-parasite interactions during co-infection, parasite epidemiology, parasite ecology and the response of parasite populations to both drug treatments and control programs.

Does the arthropod microbiota impact the establishment of vector-borne diseases in mammalian hosts?

The impact of the microbiota on the immune status of its host is a source of intense research and publicity. In comparison, the effect of arthropod microbiota on vector-borne infectious diseases has received little attention. A better understanding of the vector microbiota in relation to mammalian host immune responses is vital, as it can lead to strategies that affect transmission and improve vaccine design in a field of research where few vaccines exist and effective treatment is rare. Recent demonstrations of how microbiota decrease pathogen development in arthropods, and thus alter vector permissiveness to vector-borne diseases (VBDs), have led to renewed interest. However, hypotheses on the interactions between the arthropod-derived microbiota and the mammalian hosts have yet to be addressed. Advances in DNA sequencing technology, increased yield and falling costs, mean that these studies are now feasible for many microbiologists and entomologists. Here, we distill current knowledge and put forward key questions and experimental designs to shed light on this burgeoning research topic.

Kablammo: an interactive, web-based BLAST results visualizer.

MOTIVATION: Kablammo is a web-based application that produces interactive, vector-based visualizations of sequence alignments generated by BLAST. These visualizations can illustrate many features, including shared protein domains, chromosome structural modifications and genome misassembly. AVAILABILITY AND IMPLEMENTATION: Kablammo can be used at http://kablammo.wasmuthlab.org. For a local installation, the source code and instructions are available under the MIT license at http://github.com/jwintersinger/kablammo. CONTACT: jeff@wintersinger.org.

Morphological variability and molecular identification of Uncinaria spp. (Nematoda: Ancylostomatidae) from grizzly and black bears: new species or phenotypic plasticity?

The hookworms Uncinaria rauschi Olsen, 1968 and Uncinaria yukonensis ( Wolfgang, 1956 ) were formally described from grizzly ( Ursus arctos horribilis) and black bears ( Ursus americanus ) of North America. We analyzed the intestinal tracts of 4 grizzly and 9 black bears from Alberta and British Columbia, Canada and isolated Uncinaria specimens with anatomical traits never previously documented. We applied morphological and molecular techniques to investigate the taxonomy and phylogeny of these Uncinaria parasites. The morphological analysis supported polymorphism at the vulvar region for females of both U. rauschi and U. yukonensis. The hypothesis of morphological plasticity for U. rauschi and U. yukonensis was confirmed by genetic analysis of the internal transcribed spacers (ITS-1 and ITS-2) of the nuclear ribosomal DNA. Two distinct genotypes were identified, differing at 5 fixed sites for ITS-1 (432 base pairs [bp]) and 7 for ITS-2 (274 bp). Morphometric data for U. rauschi revealed host-related size differences: adult U. rauschi were significantly larger in black bears than in grizzly bears. Interpretation of these results, considering the historical biogeography of North American bears, suggests a relatively recent host-switching event of U. rauschi from black bears to grizzly bears which likely occurred after the end of the Wisconsin glaciation. Phylogenetic maximum parsimony (MP) and maximum likelihood (ML) analyses of the concatenated ITS-1 and ITS-2 datasets strongly supported monophyly of U. rauschi and U. yukonensis and their close relationship with Uncinaria stenocephala (Railliet, 1884), the latter a parasite primarily of canids and felids. Relationships among species within this group, although resolved by ML, were unsupported by MP and bootstrap resampling. The clade of U. rauschi, U. yukonensis, and U. stenocephala was recovered as sister to the clade represented by Uncinaria spp. from otariid pinnipeds. These results support the absence of strict host-parasite co-phylogeny for Uncinaria spp. and their carnivore hosts. Phylogenetic relationships among Uncinaria spp. provided a framework to develop the hypothesis of similar transmission patterns for the closely related U. rauschi, U. yukonensis, and U. stenocephala.