捻转血矛线虫丙硫咪唑敏感虫株和耐药虫株miRNA-mRNA关联比较分析

摘要

旨在比较分析捻转血矛线虫敏感虫株和耐药虫株miRNA表达谱,探讨miRNA与捻转血矛线虫丙硫咪唑耐药性相关基因之间的调控机制。运用Illumina Hiseq2000平台进行测序并进行cDNA文库的构建,利用RNAhybrid、Miranda和TargetScan对测序得到的数据进行靶基因预测,并取交集作为miRNA的靶基因预测结果,使用Bowtie、RepeatMasker、MIREAP、Rfam、miRBase和DAVID等生物信息学分析软件和权威数据库系统,筛选出差异的miRNA以及与捻转血矛线虫耐药相关的miRNA并进行靶基因预测,同时对KEGG pathway和GO功能富集到的靶基因和可能参与调控的通路进行预测分析。结果表明,敏感虫株和耐药虫株共筛选显著差异表达的miRNA共294个,其中113个上调,181个下调。对miRNA和其靶向调节的mRNA进行关联分析,共关联到1770个miRNA-mRNA对为负调控关系,其中涉及到274个差异的miRNA和603个差异的mRNA。显著差异表达的miRNA靶基因被注释到了1430条GO terms和327条KEGG pathway,富集到FoxO信号通路(FoxO signaling pathway),mTOR信号通路(mTOR signaling pathway),PI3K-Akt信号通路(PI3K-Akt signaling pathway)等与耐药相关的一些抗性通路。综上,通过对miRNA表达谱分析以及miRNA-mRNA靶向分析,发现捻转血矛线虫敏感虫株和耐药虫株显著差异表达的miRNA和其对应的靶基因,并对部分耐药相关通路中存在的靶基因及显著差异的靶基因所在的通路进行筛选和分析,为后续进一步通过转录组学深层次挖掘捻转血矛线虫产生耐药性的关键基因及相关调控分子提供科学依据。