捻转血矛线虫耐药虫株与敏感虫株circRNAs差异表达分析

摘要

为分析捻转血矛线虫伊维菌素敏感虫株与耐药虫株之间差异表达的circRNAs在调控耐药性中的作用，本研究以捻转血矛线虫敏感虫株和耐药虫株为研究对象，采用Illumina Hiseq 2 000平台进行高通量测序及cDNA文库的构建，利用相关生物信息学软件对测得的circRNAs进行差异性分析，同时对来源基因进行GO和KEGG pathway功能富集。运用TargetScan和Miranda软件预测关键circRNAs靶向的miRNAs，探究circRNAs-miRNAs的作用。随后又随机选取5个显著差异的circRNAs，进行RT-qPCR检测并与转录组测序结果进行一致性评定。结果显示：1）敏感虫株与耐药虫株共筛选出677个差异表达的circRNAs，其中47个显著差异，23个显著上调，24个显著下调。2）所有差异circRNAs的来源基因GO富集分析显示，共富集到生物过程336个、细胞组分59个和分子功能110个，主要与代谢过程、细胞膜结构和催化活性等相关；KEGG富集分析显示，富集到257条通路，主要与代谢（脂肪酸生物合成、甘油磷脂代谢和丙酸代谢）、胰高血糖素信号通路和AMPK信号通路等相关。3)circRNAs-miRNAs调控网络分析显示，627个差异circRNAs靶向结合187个miRNAs，组成了10 602组调控关系。综上，本研究成功筛选了捻转血矛线虫伊维菌素敏感虫株和耐药虫株间差异表达的circRNAs，对所有差异circRNAs的来源基因进行功能注释和富集分析，并构建了部分显著差异circRNAs-miRNAs调控网络。本研究结果为捻转血矛线虫耐药性产生的分子机制研究提供了新思路。

关键词

捻转血矛线虫 / 伊维菌素 / 耐药性 / 转录组 / circRNAs / 差异分析