家蚕感染二分浓核病毒（镇江株）的数字基因表达谱分析

• 作者: 尚梦珂 覃光星 郭锡杰 钱荷英 高坤
• 来源期刊: 中国农业科学
• 关键词: 家蚕 数字基因表达谱 二分浓核病毒 qRT-PCR
• 年，卷（期）: 2016,（17）
• 所属期刊栏目: 畜牧?兽医?资源昆虫
• 页码范围: 3452-3463
• 页数: 12页
• 语种: 中文
• DOI: 10.3864/j.issn.0578-1752.2016.17.018

摘要

【目的】筛选家蚕与二分浓核病毒（Bombyx mori bidensovirus Zhenjiang strain，BmBDV-ZJ）感染相关的差异表达基因，鉴定家蚕与病毒感染有关的调控基因，为进一步阐明家蚕抗BmBDV-ZJ的分子机制提供理论依据。【方法】采用Illumina高通量测序技术，构建家蚕品种JS口服感染BmBDV-ZJ的数字基因表达谱，为排除个体间差异，以10头蚕作为一个样本用于DGE检测。样本中基因的差异表达检测通过严格的运算法则进行，对差异检验的P值（P value）作多重假设检验校正，通过控制FDR（false discovery rate）来决定P值的域值。本研究中，差异表达基因定义为FDR≤0.001且差异倍数在2倍及以上（|log2ratio|≥1）的基因。采用基因本体论（GO）分类体系确定所有差异表达基因可能的功能。用 GO 计算 P值和 bonferoni校正。选用校正P值≤0.05作为基因组显著富集的阈值。WEGO软件用来视化、比较和绘制GO注释结果。利用KEGG数据库进行通路富集分析,进一步确定显著富集代谢途径或信号传导途径，Q值≤0.05的通路指定为DGEs 中的显著富集通路。通过 qRT-PCR 方法对部分差异表达基因进行验证。【结果】感染组和对照组分别得到4850663和4875307个原始标签,去除低质量标签后，分别得到4757934和4788406个清洁标签，对应的标签种类数量分别为62436和63680种。两个文库间的清洁标签和清洁标签种类的数量在不同拷贝区间分布类似，感染组和对照组样本的测序量分别为3.5 M和3.7 M，测序深度符合试验的要求，两样本的DGE数据是可信的。将这两个DGE数据库的所有清洁标签与家蚕参考基因库进行比对，在对照组与感染组中，分别有36.39%和45.30%的清洁标签可以比对到基因。另有50.02%和43.34%的清洁标签可以比对到家蚕参考基因组，剩余的未知标签分别占清洁标签总数的13.59%和12.35%。共发现了447个差异表达基因，其中306个上调表达，141个下调表达。分别有218、147和179个差异表达基因涉及GO 3个本体中的分子功能、细胞组分和生物过程。利用KEGG公共数据库进行Pathway显著性富集分析，注释到的基因总数为8473个。447个差异表达基因经鉴定后，其中的330个基因被归类到151个KEGG路径中。差异表达基因显著富集的Pathway （Q值≤0.05）有19个，其中最显著富集的是细胞质中DNA识别通路。挑选了24个差异表达基因进行qRT-PCR验证，其中20个基因的差异表达趋势与DGE的结果一致。其中在DNA识别通路中共检测到9个差异表达基因， BGIBMGA009408-TA、BGIBMGA004913-TA、BGIBMGA011753-TA均为编码RNA聚合酶III的基因，表达量均上调，是对照组的4.3、2.3、3.4倍。【结论】构建了3龄家蚕JS感染BmBDV-ZJ后28 h感染组及对照组幼虫的数字基因表达谱，Pathway显著性富集分析和qRT-PCR验证显示，家蚕感染BmBDV-ZJ后可能通过启动胞质内DNA识别通路来感应入侵病毒的异源DNA成分并迅速启动天然免疫抵御BmBDV-ZJ病毒感染,为研究BmBDV-ZJ侵染家蚕和家蚕抵御BmBDV感染的分子机制打下了基础。