Background
Haemonchus contortus (H. contortus) is the most common parasitic nematode in ruminants and is prevalent worldwide. H. contortus resistance to albendazole (ABZ) hinders the efficacy of anthelmintic drugs, but little is known about the molecular mechanisms that regulate this of drug resistance. Recent research has demonstrated that long noncoding RNAs (lncRNAs) can exert significant influence as pivotal regulators of the emergence of drug resistance.
Results
In this study, transcriptome sequencing was conducted on both albendazole-sensitive (ABZ-sensitive) and albendazole-resistant (ABZ-resistant) H. contortus strains, with three biological replicates for each group. The analysis of lncRNA in the transcriptomic data revealed that there were 276 differentially expressed lncRNA (DElncRNA) between strains with ABZ-sensitive and ABZ-resistant according to the criteria of |log2Foldchange|≥ 1 and FDR < 0.05. Notably, MSTRG.12969.2 and MSTRG.9827.1 exhibited the most significant upregulation and downregulation, respectively, in the resistant strains. The potential roles of the DElncRNAs included catalytic activity, stimulus response, regulation of drug metabolism, and modulation of the immune response. Moreover, we investigated the interactions between DElncRNAs and other RNAs, specifically MSTRG.12741.1, MSTRG.11848.1, MSTRG.5895.1, and MSTRG.14070.1, involved in regulating drug stimulation through cis/trans/antisense/lncRNA‒miRNA–mRNA interaction networks. This regulation leads to a decrease (or increase) in the expression of relevant genes, consequently enhancing the resistance of H. contortus to albendazole. Furthermore, through comprehensive analysis of competitive endogenous RNAs (ceRNAs) involved in drug resistance-related pathways, such as the mTOR signalling pathway and ABC transporter signalling pathway, the relevance of the MSTRG.2499.1-novel-m0062-3p-HCON_00099610 interaction was identified to mainly involve the regulation of catalytic activity, metabolism, ubiquitination and transcriptional regulation of gene promoters. Additionally, quantitative real-time polymerase chain reaction (qRT–PCR) validation indicated that the transcription profiles of six DElncRNAs and six DEmRNAs were consistent with those obtained by RNA-seq.
Conclusions
The results of the present study allowed us to better understand the changes in the lncRNA expression profile of ABZ-resistant H. contortus. In total, these results suggest that the lncRNAs MSTRG.963.1, MSTRG.12741.1, MSTRG.11848.1 and MSTRG.2499.1 play important roles in the development of ABZ resistance and can serve as promising biomarkers for further study.