Symbiotic bacteria are associated with ivermectin resistance in the parasitic nematode Haemonchus contortus

Abstract

Anthelmintic resistance in gastrointestinal nematodes poses a severe threat to global livestock production and sustainability. While the genomic adaptations of parasites have been extensively characterized, the role of the parasite-associated microbiome in modulating drug susceptibility remains largely unexplored. In this study, we investigated the symbiotic microbiota of ivermectin-resistant and ivermectin-susceptible strains of the blood-feeding parasitic nematode Haemonchus contortus and assessed its potential association with ivermectin resistance. Comparative microbiome analysis revealed Stenotrophomonas maltophilia as a potential resistance biomarker, with a 3.4-fold higher relative abundance in resistant strains. Furthermore, functional prediction analysis revealed that pathways related to genetic information processing and cell motility were significantly enriched in the microbiota of resistant strains. Larval mortality assays demonstrated that gentamicin treatment significantly increased the susceptibility of the larvae to ivermectin. These findings suggest that symbiotic bacteria are correlated with ivermectin resistance in Haemonchus contortus, possibly in combination with genetic factors of the parasite.

Keywords

Haemonchus contortus / Ivermectin resistance / Symbiotic bacteria / Stenotrophomonas maltophilia / Holobiont

References

[1] Aßhauer KP, Wemheuer B, Daniel R, Meinicke P. Tax4Fun: Predicting functional profiles from metagenomic 16S rRNA data. Bioinformatics (Oxford, England), 2015, 31(17): 2882-2884.

[2] Atif M, Smith JJ, Estrada-Mondragon A, Xiao X, Salim AA, Capon RJ, Lynch JW, Keramidas A. GluClR-mediated inhibitory postsynaptic currents reveal targets for ivermectin and potential mechanisms of ivermectin resistance. PLoS Pathogens, 2019, 15(1. ArticleID: e1007570

[3] Benjamini Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society: Series B, 1995, 57(1): 289-300.

[4] Bhat AH, Tak H, Ganai BA, Malik IM, Bhat TA. Bacteria associated with ovine gut parasites Trichuris ovis and Haemonchus contortus. Journal of Helminthology, 2023, 97. ArticleID: e75

[5] Bhat AH, Tak H, Ganai BA, Malik IM, Bambou JC. Beyond parasitism: Exploring the microbial profile of Haemonchus contortus and its predilection site (abomasum) in Kashmir Merino sheep. Veterinary Parasitology, 2024, 330. ArticleID: 110243

[6] Bouchet C, Deng Q, Umair S. Bacteria associated with the parasitic nematode Haemonchus contortus and its control using antibiotics. Parasitologia, 2022, 2(2): 63-70.

[7] Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecological Monographs, 1957, 27(4): 325-349.

[8] Brown MB, Forsythe AB. Robust tests for the equality of variances. Journal of the American Statistical Association, 1974, 69(4): 364-367.

[9] Bubrig LT, Sutton JM, Fierst JL. Caenorhabditis elegans dauers vary recovery in response to bacteria from natural habitat. Ecology and Evolution, 2020, 10(18): 9886-9895.

[10] Bukin YS, Galachyants YP, Morozov IV, Bukin SV, Zakharenko AS, Zemskaya TI. The effect of 16S rRNA region choice on bacterial community metabarcoding results. Scientific Data, 2019, 6. ArticleID: 190007

[11] Cabreiro F, Au C, Leung KY, Vergara-Irigaray N, Cocheme HM, Noori T, Weinkove D, Schuster E, Greene ND, Gems D. Metformin retards aging in Caenorhabditis elegans by altering microbial folate and methionine metabolism. Cell, 2013, 153(1): 228-239.

[12] Callahan BJ, Wong J, Heiner C, Oh S, Theriot CM, Gulati AS, McGill SK, Dougherty MK. High-throughput amplicon sequencing of the full-length 16S rRNA gene with single-nucleotide resolution. Nucleic Acids Research, 2019, 47(18. ArticleID: e103

[13] Campbell WC. History of avermectin and ivermectin, with notes on the history of other macrocyclic lactone antiparasitic agents. Current Pharmaceutical Biotechnology, 2012, 13(6): 853-865.

[14] Clarke KR. Nonparametric multivariate analyses of changes in community structure. Australian Journal of Ecology, 1993, 18(1): 117-143.

[15] Dirksen P, Marsh SA, Braker I, Heitland N, Wagner S, Nakad R, Mader S, Petersen C, Kowallik V, Rosenstiel P, Félix MA, Schulenburg H. The native microbiome of the nematode Caenorhabditis elegans: Gateway to a new host-microbiome model. BMC Biology, 2016, 14: 38.

[16] Doyle SR, Laing R, Bartley D, Morrison A, Holroyd N, Maitland K, Antonopoulos A, Chaudhry U, Flis I, Howell S, McIntyre J, Gilleard JS, Tait A, Mable B, Kaplan R, Sargison N, Britton C, Berriman M, Devaney E, Cotton JA. Genomic landscape of drug response reveals mediators of anthelmintic resistance. Cell Reports, 2022, 41(3. ArticleID: 111522

[17] Du AF, Li XJ, Wang SH. Cloning and sequence analysis of the 24,000 excretory-secretory antigen gene of Haemonchus contortus ZJ strain.. Chinese Journal of Veterinary Science, 2004, 24(6): 553-555. In Chinese.

[18] Dunnett CW. A multiple comparison procedure for comparing several treatments with a control. Journal of the American Statistical Association, 1955, 50: 1096-1121.

[19] Edgar RC. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nature Methods, 2013, 10(10): 996-998.

[20] El-Ashram S, Suo X. Exploring the microbial community (microflora) associated with ovine Haemonchus contortus (macroflora) field strains. Scientific Reports, 2017, 7(1. ArticleID: 70

[21] Encalada-Mena LA, Torres-Acosta JF, Sandoval-Castro CA, Reyes-Guerrero DE, Mancilla-Montelongo MG, López-Arellano R, Olmedo-Juárez A, López-Arellano ME. Comparison of P-glycoprotein gene expression of two Haemonchus contortus isolates from Yucatan, Mexico, with resistant or susceptible phenotype to ivermectin in relation to a susceptible reference strain. Veterinary Parasitology: Regional Studies and Reports, 2024, 52. ArticleID: 101047

[22] Fercoq F, Cormerais C, Remion E, Gal J, Plisson J, Fall A, Alonso J, Lhermitte-Vallarino N, Hübner MP, Kohl L, Landmann F, Martin C. Host environment shapes filarial parasite fitness and Wolbachia endosymbionts dynamics. PLOS Pathogens, 2025, 21(7. ArticleID: e1013301

[23] Fu N, Luo X, Li J, Zhou X, Du J, Zhai B, Li A, Yang X, Liu Y. Role of autophagy in resistance to ivermectin in Haemonchus contortus. Veterinary Parasitology, 2024, 330. ArticleID: 110360

[24] Gao T, Zhang Y, Sun W, Li Q, Huang X, Zhi D, Zi H, Ji R, Long Y, Gong C, Yang Y. The symbiont Wolbachia increases resistance to bifenthrin in Ectropis grisescens by regulating the host detoxification function. Ecotoxicology and Environmental Safety, 2025, 289. ArticleID: 117666

[25] Gilleard JS. Haemonchus contortus as a paradigm and model to study anthelmintic drug resistance. Parasitology, 2013, 140(12): 1506-1522.

[26] Han R, Wang Y, Deng Y, Zhang Y, Zhang L, Niu Q. Stenotrophomonas strain CPCC 101271, an intestinal lifespan-prolonging bacterium for Caenorhabditiselegans that assists in host resistance to "Bacillus nematocida” colonization. Archives of Microbiology, 2021, 203(8): 4951-4960.

[27] Han S, Akhtar MR, Xia X. Functions and regulations of insect gut bacteria. Pest Management Science, 2024, 80(10): 4828-4840.

[28] Jakobs N, Yilmaz E, Krücken J. Transgenic expression of Haemonchus contortus cytochrome P450 Hco-cyp-13A11 decreases susceptibility to particular but not all macrocyclic lactones in the model organism Caenorhabditis elegans. International Journal of Molecular Sciences, 2022, 23(16. ArticleID: 9155

[29] Jakobs N, Andreotti S, Ramünke S, von Samson-Himmelstjerna G, Krücken J. Differences in constitutive gene expression of cytochrome P450 enzymes and ATP-binding cassette transporter gene expression between a susceptible and a highly macrocyclic lactone-resistant Haemonchus contortus isolate in the absence of drug-inducible expression. Parasites & Vectors, 2024, 17(1. ArticleID: 505

[30] Johnson JS, Spakowicz DJ, Hong BY, Petersen LM, Demkowicz P, Chen L, Leopold SR, Hanson BM, Agresta HO, Gerstein M, Sodergren E, Weinstock GM. Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis. Nature Communications, 2019, 10(1. ArticleID: 5029

[31] Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research, 2000, 28: 27-30.

[32] Khan S, Nisar A, Yuan J, Luo X, Dou X, Liu F, Zhao X, Li J, Ahmad H, Mehmood SA, Feng X. A whole genome resequencing based GWA analysis reveals candidate genes associated with ivermectin resistance in Haemonchus contortus. Genes, 2020, 11(4. ArticleID: 367

[33] Kikuchi Y, Hayatsu M, Hosokawa T, Nagayama A, Tago K, Fukatsu T. Symbiont-mediated insecticide resistance. Proceedings of the National Academy of Sciences of the United States of America, 2012, 109(22): 8618-8622.

[34] Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glöckner FO. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Research, 2013, 41(1. ArticleID: e1

[35] Kotze AC, Prichard RK. Anthelmintic resistance in Haemonchus contortus: History, mechanisms and diagnosis. Advances in Parasitology, 2016, 93: 397-428.

[36] Kruskal WH, Wallis WA. Use of ranks in one-criterion variance analysis. Journal of the American Statistical Association, 1952, 47: 583-621.

[37] Laing R, Maitland K, Lecová L, Skuce PJ, Tait A, Devaney E. Analysis of putative resistance gene loci in UK field populations of Haemonchus contortus after 6 years of macrocyclic lactone use. International Journal for Parasitology, 2016, 46(10): 621-630.

[38] Laing R, Gillan V, Devaney E. Ivermectin – Old drug, new tricks?. Trends in Parasitology, 2017, 33(6): 463-472.

[39] Liu Y, Tang L, Zheng W, Tuersong W, Li J, Zhang H. Transcriptomics and proteomics of Haemonchus contortus in response to ivermectin treatment. Animals, 2023, 13(5): 919.

[40] Liu L, Zhang Z, Liu F, Liu H, Ye L, Liu F, Gupta N, Wang C, Hu M. In vitro culture of the parasitic stage larvae of hematophagous parasitic nematode Haemonchus contortus. International Journal for Parasitology, 2025, 55(5): 263-271.

[41] Louca S, Polz MF, Mazel F, Albright MBN, Huber JA, O’Connor MI, Ackermann M, Hahn AS, Srivastava DS, Crowe SA, Doebeli M, Parfrey LW. Function and functional redundancy in microbial systems. Nature Ecology & Evolution, 2018, 2(6): 936-943.

[42] Luo, X. 2021. Investigation of anthelmintic resistance and screening of resistant candidate genes in Haemonchus contortus. Doctoral dissertation, Inner Mongolia Agricultural University, Hohhot, China. https://doi.org/10.27229/d.cnki.gnmnu.2021.000018.

[43] Mafuna T, Soma P, Tsotetsi-Khambule AM, Hefer CA, Muchadeyi FC, Thekisoe OMM, Pierneef RE. Bacterial profiling of Haemonchus contortus gut microbiome infecting Dohne Merino sheep in South Africa. Scientific Reports, 2021, 11(1. ArticleID: 5905

[44] Mann HB, Whitney DR. On a test of whether one of two random variables is stochastically larger than the other. The Annals of Mathematical Statistics, 1947, 18(1): 50-60.

[45] Martin RJ, Robertson AP, Choudhary S. Ivermectin: An anthelmintic, an insecticide, and much more. Trends in Parasitology, 2021, 37(1): 48-64.

[46] Muchiut SM, Fernández AS, Steffan PE, Riva E, Fiel CA. Anthelmintic resistance: Management of parasite refugia for Haemonchus contortus through the replacement of resistant with susceptible populations. Veterinary Parasitology, 2018, 254: 43-48.

[47] Omura S, Crump A. The life and times of ivermectin - A success story. Nature Reviews Microbiology, 2004, 2(12): 984-989.

[48] Otranto D, Wall R. Veterinary parasitology, 2024, 5th EditionWiley-Blackwell.

[49] Papadopoulos E, Gallidis E, Ptochos S. Anthelmintic resistance in sheep in Europe: A selected review. Veterinary Parasitology, 2012, 189(1): 85-88.

[50] Reyes-Guerrero DE, Higuera-Piedrahita RI, Maza-Lopez J, Mendoza-de-Gives P, Camas-Pereyra R, López-Arellano ME. Analysis of P-gp genes relative expression associated with ivermectin resistance in Haemonchus contortus larval stages from in vitro cultures (L3 and xL3) and from gerbils (Meriones unguiculatus) (L4) as models of study. Journal of Helminthology, 2024, 98. ArticleID: e19

[51] Rodrigues JFV, Santos JMLD, Frota GA, Vieira LDS, Teixeira M, Monteiro MS, Monteiro JP. Expression of transporter genes in anthelmintic resistant isolates of Haemonchus contortus. Genetics and Molecular Biology, 2024, 47(3. ArticleID: e20230350

[52] Rodrigues JFV, Monteiro JP, Almeida TM, Molento MB. A systematic review of the molecular mechanisms related to anthelmintic resistance in Haemonchus contortus: A contemporary narrative. Veterinary Parasitology, 2025, 334. ArticleID: 110394

[53] Salgado JA, Santos CDP. Overview of anthelmintic resistance of gastrointestinal nematodes of small ruminants in Brazil. Revista Brasileira De Parasitologia Veterinária, 2016, 25(1): 3-17.

[54] Scheirer CJ, Ray WS, Hare N. The analysis of ranked data derived from completely randomized factorial designs. Biometrics, 1976, 32: 429-434.

[55] Shapiro SS, Wilk MB. An analysis of variance test for normality: Complete samples. Biometrika, 1965, 52(3-4): 591-611.

[56] Sinnathamby G, Henderson G, Umair S, Janssen P, Bland R, Simpson H. The bacterial community associated with the sheep gastrointestinal nematode parasite Haemonchus contortus. PLoS ONE, 2018, 13(2. ArticleID: e0192164

[57] Student. 1908. The probable error of a mean. Biometrika, 6, 1–25. https://doi.org/10.1093/biomet/6.1.1.

[58] Tipthara P, Jittamala P, Hanpithakpong W, Tarning J, Pukrittayakamee S, Hoglund RM. Identification of the metabolites of ivermectin in humans. Pharmacology Research & Perspectives, 2021, 9(1. ArticleID: e00712

[59] Tuersong W, Zhou C, Wu S, Qin P, Wang C, Di W, Liu L, Liu H, Hu M. Comparative analysis on transcriptomics of ivermectin resistant and susceptible strains of Haemonchus contortus. Parasites & Vectors, 2022, 15(1): 159.

[60] Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and Environmental Microbiology, 2007, 73(16): 5261-5267.

[61] Wang YS, Zheng XC, Hu QW, Zheng YG. Degradation of abamectin by newly isolated Stenotrophomonas maltophilia ZJB-14120 and characterization of its abamectin-tolerance mechanism. Research in Microbiology, 2015, 166(5): 408-418.

[62] Wasimuddin W, Schlaeppi K, Ronchi F, Leib SL, Erb M, Ramette A. Evaluation of primer pairs for microbiome profiling from soils to humans within the One Health framework. Molecular Ecology Resources, 2020, 20(6): 1558-1571.

[63] Wen H, Zhang Y, Yang L, Wang W, Liu C. Study of the drug resistance function of ivermectin-resistance-related miRNAs in Haemonchus contortus. Acta Parasitologica, 2025, 70(1. ArticleID: 7

[64] Yang X, Kui L, Tang M, Li D, Wei K, Chen W, Miao J, Dong Y. High-throughput transcriptome profiling in drug discovery. Frontiers in Genetics, 2020, 11: 19.

[65] Zeng T, Fu Q, Luo F, Dai J, Fu R, Qi Y, Deng X, Lu Y, Xu Y. Lactic acid bacteria modulate the CncC pathway to enhance resistance to β-cypermethrin in the oriental fruit fly.. The ISME Journal, 2024, 18(1. ArticleID: wrae058

Funding

National Key Research and Development Program of China(2023YFD1801201)

National Natural Science Foundation of China(32172881)

RIGHTS & PERMISSIONS

The Author(s)