Whole genome amplification and sequencing of individual Dirofilaria immitis microfilariae.

Dirofilaria immitis is a filarial parasitic nematode of veterinary significance. With the emergence of drug-resistant isolates in the USA, it is imperative to determine the likelihood of resistance occurring in other regions of the world. One approach is to conduct population genetic studies across an extensive geographical range, and to sequence the genomes of individual worms to understand genome-wide genetic variation associated with resistance. The immature life stages of D. immitis found in the host blood are more accessible and less invasive to sample compared to extracting adult stages from the host heart. To assess the use of immature life stages for population genetic analyses, we have performed whole genome amplification and whole-genome sequencing on nine (n = 9) individual D. immitis microfilaria samples isolated from dog blood. On average, less than 1% of mapped reads aligned to each D. immitis genome (nuclear, mitochondrial, and Wolbachia endosymbiont). For the dog genome, an average of over 99% of mapped reads aligned to the nuclear genome and less than 1% aligned to the mitochondrial genome. The average coverage for all D. immitis genomes and the dog nuclear genome was less than 1, while the dog mitochondrial genome had an average coverage of 2.87. The overwhelming proportion of sequencing reads mapping to the dog host genome can be attributed to residual dog blood cells in the microfilariae samples. These results demonstrate the challenges of conducting genome-wide studies on individual immature parasite life stages, particularly in the presence of extraneous host DNA.

Genome structure and population genomics of the canine heartworm Dirofilaria immitis.

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