The recent identification of isolates of D. immitis with confirmed resistance to the macrocyclic lactone preventatives presents an opportunity for comparative genomic studies using these isolates, and examining the genetic diversity within and between them. We studied the genomes of Wolbachia endosymbionts of five isolates of D. immitis maintained at the University of Georgia.Missouri and Georgia-2 are maintained as drug susceptible isolates, and JYD-27, Yazoo-2013 and Metairie-2014 are resistant to the macrocyclic lactone preventatives. We used whole genome amplification followed by Illumina-based sequencing from 8-12 individual microfilariae from each of the five isolates, obtaining a depth of coverage of approximately 40-75 fold for each. The Illumina sequences were used to create new genome assemblies for all the Wolbachia isolates studied. Comparisons of the Wolbachia sequences revealed more than 3,000 sequence variations in each isolate. We identified 67 loci specific in resistant isolates but not in susceptible isolates, including 18 genes affected. Phylogenetic analysis suggested that the endosymbionts of the drugsusceptible isolates are more closely related to each other than to those from any of the resistant parasites. This level of variation in the Wolbachia endosymbionts of D. immitis isolates suggests a potential for selection for resistance against drugs targeting them.

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 X-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 SNPs previously associated with macrocyclic lactone resistance in the new genome assembly, revealing 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.

No abstract