Gene co-expression network analysis reveal core responsive genes in Parascaris univalens tissues following ivermectin exposure.
Abstract: Anthelmintic resistance in equine parasite Parascaris univalens, compromises ivermectin (IVM) effectiveness and necessitates an in-depth understanding of its resistance mechanisms. Most research, primarily focused on holistic gene expression analyses, may overlook vital tissue-specific responses and often limit the scope of novel genes. This study leveraged gene co-expression network analysis to elucidate tissue-specific transcriptional responses and to identify core genes implicated in the IVM response in P. univalens. Adult worms (n = 28) were exposed to 10-11 M and 10-9 M IVM in vitro for 24 hours. RNA-sequencing examined transcriptional changes in the anterior end and intestine. Differential expression analysis revealed pronounced tissue differences, with the intestine exhibiting substantially more IVM-induced transcriptional activity. Gene co-expression network analysis identified seven modules significantly associated with the response to IVM. Within these, 219 core genes were detected, largely expressed in the intestinal tissue and spanning diverse biological processes with unspecific patterns. After 10-11 M IVM, intestinal tissue core genes showed transcriptional suppression, cell cycle inhibition, and ribosomal alterations. Interestingly, genes PgR028_g047 (sorb-1), PgB01_g200 (gmap-1) and PgR046_g017 (col-37 & col-102) switched from downregulation at 10-11 M to upregulation at 10-9 M IVM. The 10-9 M concentration induced expression of cuticle and membrane integrity core genes in the intestinal tissue. No clear core gene patterns were visible in the anterior end after 10-11 M IVM. However, after 10-9 M IVM, the anterior end mostly displayed downregulation, indicating disrupted transcriptional regulation. One interesting finding was the non-modular calcium-signaling gene, PgR047_g066 (gegf-1), which uniquely connected 71 genes across four modules. These genes were enriched for transmembrane signaling activity, suggesting that PgR047_g066 (gegf-1) could have a key signaling role. By unveiling tissue-specific expression patterns and highlighting biological processes through unbiased core gene detection, this study reveals intricate IVM responses in P. univalens. These findings suggest alternative drug uptake of IVM and can guide functional validations to further IVM resistance mechanism understanding.
Copyright: © 2024 Dube et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Publication Date: 2024-02-15 PubMed ID: 38359071PubMed Central: PMC10868809DOI: 10.1371/journal.pone.0298039Google Scholar: Lookup
The Equine Research Bank provides access to a large database of publicly available scientific literature. Inclusion in the Research Bank does not imply endorsement of study methods or findings by Mad Barn.
- Journal Article
Summary
This research summary has been generated with artificial intelligence and may contain errors and omissions. Refer to the original study to confirm details provided. Submit correction.
The study investigates the resistance mechanisms of the equine parasite, Parascaris univalens, towards the common anthelmintic drug ivermectin (IVM) by exploring tissue-specific transcriptional responses and identifying critical genes involved.
Research Methodology
- P. univalens (adult worms, n = 28) were exposed to two various concentrations of IVM (10-11 M & 10-9 M) in a controlled lab setting over a 24-hour period.
- Post-exposure, the researchers utilized RNA-sequencing to observe transcriptional changes in two specific worm tissues: the anterior end and the intestine.
- Upon analyzing differential expression, they detected pronounced tissue differences, noting greater IVM-induced transcriptional activity in the intestinal tissue compared to the anterior end.
Findings
- Through gene co-expression network analysis, the researchers identified seven modules that exhibited significant associations with IVM response.
- Within these modules, 219 core genes were found, most of which exhibited expression in the intestinal tissue and were involved in diverse biological processes.
- When exposed to the lower concentration of IVM (10-11M), these intestinal core genes showed a decrease in transcription activity, cell cycle inhibition, and changes in the ribosomal structure.
- However, certain genes (PgR028_g047, PgB01_g200, and PgR046_g017) displayed behavior contrary to the general pattern, switching from downregulation to upregulation when the IVM concentration was increased to 10-9 M.
- Under higher IVM concentration (10-9M), core genes associated with cuticle and membrane integrity exhibited enhanced expression in the intestinal tissue. No consistent patterns were observed in the anterior end under 10-11M IVM; however, it mainly showed downregulation after exposure to 10-9M IVM.
- Another unique discovery was the calcium-signaling gene, PgR047_g066 (gegf-1), which was found to connect with 71 other genes across four different modules. These connected genes were generally enriched for transmembrane signaling activity, suggesting a critical signaling role for the PgR047_g066 (gegf-1) gene.
Conclusion
- The study highlights tissue-specific expression patterns in response to IVM, thus shedding light on the complexities of IVM responses in P. univalens.
- The findings indicate potential alternate methods of drug uptake of IVM by the parasite, promising to guide further research into functional validation and understanding resistance mechanisms.
- P. univalens (adult worms, n = 28) were exposed to two various concentrations of IVM (10-11 M & 10-9 M) in a controlled lab setting over a 24-hour period.
- Post-exposure, the researchers utilized RNA-sequencing to observe transcriptional changes in two specific worm tissues: the anterior end and the intestine.
- Upon analyzing differential expression, they detected pronounced tissue differences, noting greater IVM-induced transcriptional activity in the intestinal tissue compared to the anterior end.
Findings
- Through gene co-expression network analysis, the researchers identified seven modules that exhibited significant associations with IVM response.
- Within these modules, 219 core genes were found, most of which exhibited expression in the intestinal tissue and were involved in diverse biological processes.
- When exposed to the lower concentration of IVM (10-11M), these intestinal core genes showed a decrease in transcription activity, cell cycle inhibition, and changes in the ribosomal structure.
- However, certain genes (PgR028_g047, PgB01_g200, and PgR046_g017) displayed behavior contrary to the general pattern, switching from downregulation to upregulation when the IVM concentration was increased to 10-9 M.
- Under higher IVM concentration (10-9M), core genes associated with cuticle and membrane integrity exhibited enhanced expression in the intestinal tissue. No consistent patterns were observed in the anterior end under 10-11M IVM; however, it mainly showed downregulation after exposure to 10-9M IVM.
- Another unique discovery was the calcium-signaling gene, PgR047_g066 (gegf-1), which was found to connect with 71 other genes across four different modules. These connected genes were generally enriched for transmembrane signaling activity, suggesting a critical signaling role for the PgR047_g066 (gegf-1) gene.
Conclusion
- The study highlights tissue-specific expression patterns in response to IVM, thus shedding light on the complexities of IVM responses in P. univalens.
- The findings indicate potential alternate methods of drug uptake of IVM by the parasite, promising to guide further research into functional validation and understanding resistance mechanisms.
Cite This Article
APA
Dube F, Delhomme N, Martin F, Hinas A, Åbrink M, Svärd S, Tydén E.
(2024).
Gene co-expression network analysis reveal core responsive genes in Parascaris univalens tissues following ivermectin exposure.
PLoS One, 19(2), e0298039.
https://doi.org/10.1371/journal.pone.0298039 Publication
Researcher Affiliations
- Department of Animal Biosciences, Swedish University of Agricultural Sciences, Uppsala, Sweden.
- Umeå Plant Science Centre (UPSC), Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, Umeå, Sweden.
- Department of Animal Biosciences, Swedish University of Agricultural Sciences, Uppsala, Sweden.
- Department of Cell and Molecular Biology, Uppsala University, Uppsala Sweden.
- Department of Animal Biosciences, Swedish University of Agricultural Sciences, Uppsala, Sweden.
- Department of Cell and Molecular Biology, Uppsala University, Uppsala Sweden.
- Department of Animal Biosciences, Swedish University of Agricultural Sciences, Uppsala, Sweden.
MeSH Terms
- Horses / genetics
- Animals
- Ivermectin / pharmacology
- Anthelmintics / pharmacology
- Gene Expression Regulation
- Gene Expression Profiling
- Ascaridoidea / genetics
- Drug Resistance / genetics
Conflict of Interest Statement
The authors have declared that no competing interests exist.
References
This article includes 72 references
- Gilleard JS, Kotze AC, Leathwick D, Nisbet AJ, McNeilly TN, Besier B. A journey through 50 years of research relevant to the control of gastrointestinal nematodes in ruminant livestock and thoughts on future directions.. Int J Parasitol 2021;51(13–14):1133–51.
- Veesenmeyer AF. Important Nematodes in Children.. Pediatr Clin North Am 2022;69(1):129–39.
- Feng XP, Hayashi J, Beech RN, Prichard RK. Study of the nematode putative GABA type-A receptor subunits: evidence for modulation by ivermectin.. J Neurochem 2002;83(4):870–8.
- Wolstenholme AJ, Neveu C. The avermectin/milbemycin receptors of parasitic nematodes.. Pestic Biochem Physiol 2022;181:105010.
- Wolstenholme AJ, Rogers AT. Glutamate-gated chloride channels and the mode of action of the avermectin/milbemycin anthelmintics.. Parasitology 2005;131 Suppl:S85–95.
- Reinemeyer CR. Diagnosis and control of anthelmintic-resistant Parascaris equorum.. Parasit Vectors 2009;2 Suppl 2(Suppl 2):S8.
- Cribb NC, Cote NM, Bouré LP, Peregrine AS. Acute small intestinal obstruction associated with Parascaris equorum infection in young horses: 25 cases (1985–2004).. N Z Vet J 2006;54(6):338–43.
- Clayton HM, Duncan JL. Clinical signs associated with parascaris equorum infection in worm-free pony foals and yearlings.. Veterinary Parasitology 1978;4(1):69–78.
- Boersema JH, Eysker M, Nas JW. Apparent resistance of Parascaris equorum to macrocylic lactones.. Vet Rec 2002;150(9):279–81.
- Peregrine AS, Molento MB, Kaplan RM, Nielsen MK. Anthelmintic resistance in important parasites of horses: does it really matter?. Vet Parasitol 2014;201(1–2):1–8.
- Whittaker JH, Carlson SA, Jones DE, Brewer MT. Molecular mechanisms for anthelmintic resistance in strongyle nematode parasites of veterinary importance.. J Vet Pharmacol Ther 2017;40(2):105–15.
- Doyle SR, Laing R, Bartley D, Morrison A, Holroyd N, Maitland K. Genomic landscape of drug response reveals mediators of anthelmintic resistance.. Cell Rep 2022;41(3):111522.
- Laing R, Doyle SR, McIntyre J, Maitland K, Morrison A, Bartley DJ. Transcriptomic analyses implicate neuronal plasticity and chloride homeostasis in ivermectin resistance and response to treatment in a parasitic nematode.. PLoS Pathog 2022;18(6):e1010545.
- Ménez C, Alberich M, Courtot E, Guegnard F, Blanchard A, Aguilaniu H. The transcription factor NHR-8: A new target to increase ivermectin efficacy in nematodes.. PLoS Pathog 2019;15(2):e1007598.
- Dent JA, Smith MM, Vassilatis DK, Avery L. The genetics of ivermectin resistance in Caenorhabditis elegans.. Proceedings of the National Academy of Sciences 2000;97(6):2674–9.
- Njue AI, Hayashi J, Kinne L, Feng XP, Prichard RK. Mutations in the extracellular domains of glutamate-gated chloride channel alpha3 and beta subunits from ivermectin-resistant Cooperia oncophora affect agonist sensitivity.. J Neurochem 2004;89(5):1137–47.
- McCavera S, Rogers AT, Yates DM, Woods DJ, Wolstenholme AJ. An ivermectin-sensitive glutamate-gated chloride channel from the parasitic nematode Haemonchus contortus.. Mol Pharmacol 2009;75(6):1347–55.
- Williamson SM, Storey B, Howell S, Harper KM, Kaplan RM, Wolstenholme AJ. Candidate anthelmintic resistance-associated gene expression and sequence polymorphisms in a triple-resistant field isolate of Haemonchus contortus.. Mol Biochem Parasitol 2011;180(2):99–105.
- Ardelli BF, Prichard RK. Inhibition of P-glycoprotein enhances sensitivity of Caenorhabditis elegans to ivermectin.. Vet Parasitol 2013;191(3–4):264–75.
- Dicker AJ, Nisbet AJ, Skuce PJ. Gene expression changes in a P-glycoprotein (Tci-pgp-9) putatively associated with ivermectin resistance in Teladorsagia circumcincta.. Int J Parasitol 2011;41(9):935–42.
- Janssen IJ, Krücken J, Demeler J, Basiaga M, Kornaś S, von Samson-Himmelstjerna G. Genetic variants and increased expression of Parascaris equorum P-glycoprotein-11 in populations with decreased ivermectin susceptibility.. PLoS One 2013;8(4):e61635.
- Martin F, Dube F, Karlsson Lindsjö O, Eydal M, Höglund J, Bergström TF. Transcriptional responses in Parascaris univalens after in vitro exposure to ivermectin, pyrantel citrate and thiabendazole.. Parasit Vectors 2020;13(1):342.
- Clayton HM, Duncan JL. The migration and development of Parascaris equorum in the horse.. International Journal for Parasitology 1979;9(4):285–92.
- . Comparative genomics of the major parasitic worms.. Nat Genet 2019;51(1):163–74.
- Thompson DP, Ho NFH, Sims SM, Geary TG. Mechanistic approaches to quantitate anthelmintic absorption by gastrointestinal nematodes.. Parasitology Today 1993;9(1):31–5.
- Yin Y, Martin J, Abubucker S, Scott AL, McCarter JP, Wilson RK. Intestinal transcriptomes of nematodes: comparison of the parasites Ascaris suum and Haemonchus contortus with the free-living Caenorhabditis elegans.. PLoS Negl Trop Dis 2008;2(8):e269.
- Jasmer DP, Rosa BA, Tyagi R, Mitreva M. Omics Driven Understanding of the Intestines of Parasitic Nematodes.. Front Genet 2019;10:652.
- Schafer W. Nematode nervous systems.. Current Biology 2016;26(20):R955–R9.
- Wang J. Genome Analysis of Programmed DNA Elimination in Parasitic Nematodes.. Methods Mol Biol 2021;2369:251–61.
- Ballesteros C, Tritten L, O’Neill M, Burkman E, Zaky WI, Xia J. The Effects of Ivermectin on Brugia malayi Females In Vitro: A Transcriptomic Approach.. PLoS Negl Trop Dis 2016;10(8):e0004929.
- Liu Y, Wang X, Luo X, Wang R, Zhai B, Wang P. Transcriptomics and Proteomics of Haemonchus contortus in Response to Ivermectin Treatment.. Animals (Basel) 2023;13(5).
- Dube F, Hinas A, Roy S, Martin F, Åbrink M, Svärd S. Ivermectin-induced gene expression changes in adult Parascaris univalens and Caenorhabditis elegans: a comparative approach to study anthelminthic metabolism and resistance in vitro.. Parasit Vectors 2022;15(1):158.
- Dube F, Hinas A, Delhomme N, Åbrink M, Svärd S, Tydén E. Transcriptomics of ivermectin response in Caenorhabditis elegans: Integrating abamectin quantitative trait loci and comparison to the Ivermectin-exposed DA1316 strain.. PLoS One 2023;18(5):e0285262.
- Karim MR, Beyan O, Zappa A, Costa IG, Rebholz-Schuhmann D, Cochez M. Deep learning-based clustering approaches for bioinformatics.. Briefings in bioinformatics 2021;22(1):393–415.
- Cao Y, Geddes TA, Yang JYH, Yang P. Ensemble deep learning in bioinformatics.. Nature Machine Intelligence 2020;2(9):500–8.
- Sohail A, Arif F. Supervised and unsupervised algorithms for bioinformatics and data science.. Progress in biophysics and molecular biology 2020;151:14–22.
- Schiffthaler B, van Zalen E, Serrano AR, Street NR, Delhomme N. Seiðr: Efficient calculation of robust ensemble gene networks.. Heliyon 2023;9(6):e16811.
- Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis.. BMC Bioinformatics 2008;9(1):559.
- van Dam S, Võsa U, van der Graaf A, Franke L, de Magalhães JP. Gene co-expression analysis for functional classification and gene-disease predictions.. Brief Bioinform 2018;19(4):575–92.
- Martin F, Eydal M, Höglund J, Tydén E. Constitutive and differential expression of transport protein genes in Parascaris univalens larvae and adult tissues after in vitro exposure to anthelmintic drugs.. Vet Parasitol 2021;298:109535.
- Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data.. Bioinformatics 2012;28(24):3211–7.
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data.. Bioinformatics 2014;30(15):2114–20.
- de Sena Brandine G, Smith AD. Falco: high-speed FastQC emulation for quality control of sequencing data.. F1000Res 2019;8:1874.
- Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression.. Nat Methods 2017;14(4):417–9.
- Howe KL, Bolt BJ, Shafie M, Kersey P, Berriman M. WormBase ParaSite − a comprehensive resource for helminth genomics.. Molecular and Biochemical Parasitology 2017;215:2–10.
- Dube F. ruqse/Parascaris-IVM-GeneNetwork: Parascaris-IVM-GeneNetwork.. Zenodo 2024.
- R Core Team A, Team RC. R: A language and environment for statistical computing.. R Foundation for Statistical Computing, Vienna, Austria 2012. 2022.
- Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences.. F1000Res 2015;4:1521.
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.. Genome Biol 2014;15(12):550.
- Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W. limma powers differential expression analyses for RNA-sequencing and microarray studies.. Nucleic Acids Res 2015;43(7):e47.
- Schurch NJ, Schofield P, Gierliński M, Cole C, Sherstnev A, Singh V. How many biological replicates are needed in an RNA-seq experiment and which differential expression tool should you use?. Rna 2016;22(6):839–51.
- Kolberg L, Raudvere U, Kuzmin I, Vilo J, Peterson H. gprofiler2—an R package for gene list functional enrichment analysis and namespace conversion toolset g:Profiler.. F1000Res 2020;9.
- Zhong R, Allen JD, Xiao G, Xie Y. Ensemble-based network aggregation improves the accuracy of gene network reconstruction.. PLoS One 2014;9(11):e106319.
- Pages H, Carlson M, Falcon S, Li N, Maintainer MBP. Package ‘AnnotationDbi’.. Bioconductor Packag Maint 2017.
- Traag VA, Waltman L, Van Eck NJ. From Louvain to Leiden: guaranteeing well-connected communities.. Scientific reports 2019;9(1):5233.
- Csardi G, Nepusz T. The igraph software package for complex network research.. InterJournal, complex systems 2006;1695(5):1–9.
- Scare JA, Dini P, Norris JK, Steuer AE, Scoggin K, Gravatte HS. Ascarids exposed: a method for in vitro drug exposure and gene expression analysis of anthelmintic naïve Parascaris spp.. Parasitology 2020;147(6):659–66.
- Gerhard AP, Krücken J, Heitlinger E, Janssen IJI, Basiaga M, Kornaś S. The P-glycoprotein repertoire of the equine parasitic nematode Parascaris univalens.. Sci Rep 2020;10(1):13586.
- Page A, Johnstone I. The cuticle, WormBook, ed. The C. elegans Research Community WormBook.. 2007.
- Njume FN, Razzauti A, Soler M, Perschin V, Fazeli G, Bourez A. A lipid transfer protein ensures nematode cuticular impermeability.. iScience 2022;25(11):105357.
- Zhang ZY. Protein-tyrosine phosphatases: biological function, structural characteristics, and mechanism of catalysis.. Crit Rev Biochem Mol Biol 1998;33(1):1–52.
- Dicker AJ, Nath M, Yaga R, Nisbet AJ, Lainson FA, Gilleard JS. Teladorsagia circumcincta: the transcriptomic response of a multi-drug-resistant isolate to ivermectin exposure in vitro.. Exp Parasitol 2011;127(2):351–6.
- Ruiz-May E, Álvarez-Sánchez ME, Aguilar-Tipacamú G, Elizalde-Contreras JM, Bojórquez-Velázquez E, Zamora-Briseño JA. Comparative proteome analysis of the midgut of Rhipicephalus microplus (Acari: Ixodidae) strains with contrasting resistance to ivermectin reveals the activation of proteins involved in the detoxification metabolism.. J Proteomics 2022;263:104618.
- Stasiuk SJ, MacNevin G, Workentine ML, Gray D, Redman E, Bartley D. Similarities and differences in the biotransformation and transcriptomic responses of Caenorhabditis elegans and Haemonchus contortus to five different benzimidazole drugs.. International Journal for Parasitology: Drugs and Drug Resistance 2019;11:13–29.
- Khan S, Nisar A, Yuan J, Luo X, Dou X, Liu F. A Whole Genome Re-Sequencing Based GWA Analysis Reveals Candidate Genes Associated with Ivermectin Resistance in Haemonchus contortus.. Genes (Basel) 2020;11(4).
- Hartman JH, Widmayer SJ, Bergemann CM, King DE, Morton KS, Romersi RF. Xenobiotic metabolism and transport in Caenorhabditis elegans.. J Toxicol Environ Health B Crit Rev 2021;24(2):51–94.
- Loveless T, Qadota H, Benian GM, Hardin J. Caenorhabditis elegans SORB-1 localizes to integrin adhesion sites and is required for organization of sarcomeres and mitochondria in myocytes.. Mol Biol Cell 2017;28(25):3621–33.
- Hammarlund M, Hobert O, Miller DM, 3rd, Sestan N. The CeNGEN Project: The Complete Gene Expression Map of an Entire Nervous System.. Neuron 2018;99(3):430–3.
- Airs PM, Vaccaro K, Gallo KJ, Dinguirard N, Heimark ZW, Wheeler NJ. Spatial transcriptomics reveals antiparasitic targets associated with essential behaviors in the human parasite Brugia malayi.. PLoS Pathog 2022;18(4):e1010399.
- Rosa BA, Jasmer DP, Mitreva M. Genome-Wide Tissue-Specific Gene Expression, Co-expression and Regulation of Co-expressed Genes in Adult Nematode Ascaris suum.. PLOS Neglected Tropical Diseases 2014;8(2):e2678.
- Alvarez J, Alvarez-Illera P, García-Casas P, Fonteriz RI, Montero M. The Role of Ca2+ Signaling in Aging and Neurodegeneration: Insights from Caenorhabditis elegans Models.. Cells 2020;9(1):204.
- Choi YJ, Bisset SA, Doyle SR, Hallsworth-Pepin K, Martin J, Grant WN. Genomic introgression mapping of field-derived multiple-anthelmintic resistance in Teladorsagia circumcincta.. PLoS Genet 2017;13(6):e1006857.
Citations
This article has been cited 0 times.Use Nutrition Calculator
Check if your horse's diet meets their nutrition requirements with our easy-to-use tool Check your horse's diet with our easy-to-use tool
Talk to a Nutritionist
Discuss your horse's feeding plan with our experts over a free phone consultation Discuss your horse's diet over a phone consultation
Submit Diet Evaluation
Get a customized feeding plan for your horse formulated by our equine nutritionists Get a custom feeding plan formulated by our nutritionists