Chromosomal genome assembly resolves drug resistance loci in the parasitic nematode Teladorsagia circumcincta

The parasitic nematode Teladorsagia circumcincta is one of the most important pathogens of sheep and goat farming in temperate climates worldwide and can rapidly evolve resistance to drugs used to control it. To understand the genetics of drug resistance, we have generated a highly contiguous genome assembly for the UK T. circumcincta isolate, MTci2. Assembly using PacBio long-reads and Hi-C long-molecule scaffolding together with manual curation resulted in a 573 Mb assembly (N50 = 84 Mb, n = 1,286) with five autosomal and one sex-linked chromosomal-scale scaffolds consistent with its karyotype. The genome resource was further improved via annotation of 22,948 genes, with manual curation of over 3,200 of these, resulting in a robust and near complete resource (96.3% complete protein BUSCOs) to support basic and applied research on this important veterinary pathogen. Genome-wide analyses of drug resistance, combining evidence from three distinct experiments, identified selection around known candidate genes for benzimidazole, levamisole and ivermectin resistance, as well as novel regions associated with ivermectin and moxidectin resistance. These insights into contemporary and historic genetic selection further emphasise the importance of contiguous genome assemblies in interpreting genome-wide genetic variation associated with drug resistance and identified key loci to prioritise in developing diagnostic markers of anthelmintic resistance to support parasite control. AUTHOR SUMMARYUnderstanding the genetics of anthelmintic resistance is a critical part of the sustainable control of parasitic worms. Here, we have generated a chromosome-scale genome assembly for the gastrointestinal nematode, Teladorsagia circumcincta, one of the most important pathogens of sheep in temperate climates worldwide. This genome and its annotation offers a substantial improvement over existing genetic resources, and will enable new insight into the biology for this parasite and new opportunities to characterise therapeutic targets such as drug and vaccine candidates. We used this resource to map genetic variation associated with resistance to multiple anthelmintic drug classes used as the primary means of parasite control, confirming known candidate genes and variants (i.e., beta-tubulin isotypes 1 and 2 associated with benzimidazole resistance, and acr-8 associated with levamisole resistance, and pgp-9 copy number variation, of which overexpression is associated broadly with anthelmintic resistance) and revealing new regions of the genome associated with drug treatment responses (i.e. 34-38 Mb on chromosome 5 associated with ivermectin resistance). Our results highlight the importance of a contiguous genome assembly and the use of complementary experimental approaches to reveal the impact of drug-mediated selection, both historically and directly in response to treatment, on genome-wide genetic variation.