Naturally occurring variation in a cytochrome P450 modifies thiabendazole responses independent of beta-tubulin

Widespread anthelmintic resistance has complicated the management of parasitic nematodes. Resistance to the benzimidazole (BZ) drug class is nearly ubiquitous in many species and is associated with mutations in beta-tubulin genes. However, mutations in beta-tubulin alone do not fully explain all BZ resistance. We performed a genome-wide association study using a genetically diverse panel of Caenorhabditis elegans strains to identify loci that contribute to resistance to the BZ drug thiabendazole (TBZ). We identified a quantitative trait locus (QTL) on chromosome V independent of all beta-tubulin genes and overlapping with two promising candidate genes, the cytochrome P450 gene cyp-35d1 and the nuclear hormone receptor nhr-176, identified by another mapping technique. Both genes were previously demonstrated to play a role in TBZ metabolism. NHR-176 binds TBZ and induces the expression of CYP-35D1, which metabolizes TBZ. We generated single gene deletions of nhr-176 and cyp-35d1 and found that both genes play a role in TBZ response. A predicted high-impact lysine-to-glutamate substitution at position 267 (K267E) in CYP-35D1 was identified in a sensitive parental strain, and reciprocal allele replacement strains in both genetic backgrounds were used to show that the lysine allele conferred increased TBZ resistance. Using competitive fitness assays, we found that neither allele is deleterious, but the lysine allele is selected in the presence of TBZ. Additionally, we found that the lysine allele significantly increased the rate of TBZ metabolism compared to the glutamate allele. Moreover, yeast expression assays showed that the lysine version of CYP-35D1 had twice the enzymatic activity of the glutamate allele. To connect our results to parasitic nematodes, we analyzed four Haemonchus contortus cytochrome P450 orthologs but did not find variation at the 267 position in fenbendazole-resistant populations. Overall, we confirmed that variation in this cytochrome P450 gene is the first locus independent of beta-tubulin to play a role in BZ resistance. Author SummaryBenzimidazoles (BZs) are the most common drug class used to control parasitic nematodes, but because of overuse, resistance is widespread. The known genetic causes of BZ resistance are associated with mutations in beta-tubulin and are the most well understood of any anthelmintic class. However, BZ response varies significantly and differential levels of resistance likely require mutations in genes independent of beta-tubulin. We used the free-living model nematode Caenorhabditis elegans to identify and characterize a novel cytochrome P450 gene, cyp-35d1, associated with natural resistance to the BZ drug thiabendazole (TBZ). We demonstrated that a lysine at position 267 confers TBZ resistance and is selected over multiple generations after TBZ treatment. This allele significantly increased the rate of TBZ metabolism in both C. elegans and yeast. In conclusion, we have characterized the role of variation in a cytochrome P450 that contributes to TBZ resistance, independent of mutations in beta-tubulin.