Embrapa Pecuária Sudeste, Rodovia Washington Luiz, km 234, Fazenda Canchim, São Carlos, SP, CEP 13560-970, Brazil.
NGS Soluções Genômicas, Rua Ajudante Albano, 847, Piracicaba, SP, CEP 13416-030, Brazil.
Parasit Vectors. 2019 Aug 14;12(1):403. doi: 10.1186/s13071-019-3663-9.
Haemonchus contortus, a gastrointestinal nematode parasite of sheep, is mainly controlled by anthelmintics; the occurrence of anthelmintic resistance leads to treatment failures and increases economic burden. Because molecular mechanisms involved in drug resistance can be elucidated by genomic studies, an extreme quantitative trait locus (X-QTL) mapping approach was used to identify co-segregation of the resistance phenotype with genetic markers to detect the genome-wide variants associated with monepantel resistance in H. contortus.
A cross between H. contortus isolates using parental susceptible (Par-S) males and monepantel resistant (Par-R) females resulted in SR progeny, while reciprocal cross resulted in RS progeny. Pools (n = 30,000) of infective larvae (L3) recovered from Par-R, and from SR and RS populations in the F3 generation, collected both before (unselected group) and 7 days after (selected group) selection with monepantel treatment in sheep hosts, were subjected to genome sequencing (Pool-Seq). Pairwise comparisons of allele frequencies between unselected and selected groups were performed for each population by Fisher's exact test (FET) and for both populations combined by a Cochran-Mantel-Haenszel (CMH) test.
Mapping rates varied from 80.29 to 81.77% at a 90.4X mean coverage of aligned reads. After correction for multiple testing, significant (P < 0.05) changes in allele frequencies were detected by FET for 6 and 57 single nucleotide polymorphisms (SNPs) in the SR and RS populations, respectively, and by the CMH test for 124 SNPs in both populations. The significant variants located on chromosome 2 generated a selection signal in a genomic region harboring the mptl-1, deg-3 and des-2 genes, previously reported as candidates for monepantel resistance. In addition, three new variants were identified in the mptl-1 gene.
This study expands knowledge on genome-wide molecular events underlying H. contortus resistance to monepantel. The identification of a genome region harboring major genes previously associated with monepantel resistance supports the results of the employed X-QTL approach. In addition, a deletion in exon 11 of the mptl-1 gene should be further investigated as the putative causal mutation leading to monepantel resistance.
绵羊胃肠道线虫寄生虫捻转血矛线虫主要通过驱虫药控制;驱虫药耐药的发生导致治疗失败并增加经济负担。由于基因组研究可以阐明耐药相关的分子机制,因此采用极端数量性状位点(X-QTL)作图方法来鉴定耐药表型与遗传标记的共分离,以检测与捻转血矛线虫对莫能菌素耐药相关的全基因组变异。
使用对莫能菌素敏感(Par-S)雄性和莫能菌素耐药(Par-R)雌性亲本进行捻转血矛线虫杂交,产生 SR 后代,而反向杂交则产生 RS 后代。从 Par-R 以及 F3 代的 SR 和 RS 群体中回收的感染性幼虫(L3)池(n=30000),分别在未选择(未选择组)和绵羊宿主中用莫能菌素处理后 7 天(选择组)进行基因组测序(Pool-Seq)。通过 Fisher 精确检验(FET)对每个群体的未选择和选择组之间的等位基因频率进行成对比较,通过 Cochran-Mantel-Haenszel(CMH)检验对两个群体进行联合比较。
在对齐读取的 90.4X 平均覆盖率下,映射率从 80.29%到 81.77%不等。经过多次测试校正后,FET 在 SR 和 RS 群体中分别检测到 6 和 57 个单核苷酸多态性(SNP)的等位基因频率发生显著变化(P<0.05),CMH 检验在两个群体中检测到 124 个 SNP 的显著变化。位于染色体 2 上的显著变异在一个包含 mptl-1、deg-3 和 des-2 基因的基因组区域产生了选择信号,这些基因先前被报道为莫能菌素耐药的候选基因。此外,mptl-1 基因中还鉴定到三个新的变异。
本研究扩展了捻转血矛线虫对莫能菌素耐药的全基因组分子事件的认识。鉴定到一个包含先前与莫能菌素耐药相关的主要基因的基因组区域,支持了所采用的 X-QTL 方法的结果。此外,mptl-1 基因外显子 11 中的缺失应进一步研究,作为导致莫能菌素耐药的潜在因果突变。