论文部分内容阅读
昆虫作为地球上种类最多的动物类群,在生态系统上占据着十分重要的地位。昆虫在生长发育过程中会遭遇各种不利的非生物胁迫因子,如高温、低温、重金属、饥饿和紫外线辐射等。当其感受到逆境信号后通过信号传导过程调节细胞内抗逆蛋白的表达,调整自身行为或生理状态适应不利环境。本文以中华蜜蜂(Apis ceranacerana,简称中蜂)为实验对象,研究了饥饿胁迫及紫外线A(UltravioletA,UVA)胁迫不同时间对中蜂存活率、体重及血淋巴蛋白含量的影响,并利用RNA-seq技术,分别对中蜂响应饥饿胁迫及UVA胁迫的免疫应答机制进行深入研究。主要结果如下:1.中蜂经不同时间的饥饿胁迫后,存活率显著下降(P<0.05);短时间饥饿胁迫(0-3 h)对中蜂体重无显著影响(P>0.05),随着饥饿胁迫时间的延长(3 h以后),中蜂的体重显著下降(P<0.05)。中蜂经不同时间的UVA胁迫后,存活率显著下降(P<0.05);短时间胁迫(0-6 h)对中蜂体重无显著影响(P>0.05);随着UVA胁迫时间的延长(6h以后),中蜂的体重显著减轻(P<0.05)。2.中蜂经不同时间饥饿及UVA胁迫后,血淋巴蛋白含量均呈现先升高后下降的趋势,且处理组与对照组的血淋巴蛋白含量有显著性差异(P<0.05)。饥饿胁迫9h时,中蜂血淋巴蛋白含量最高,为76.58±0.60 g/L;UVA胁迫12h时,中蜂血淋巴蛋白含量最高,为93.56±2.71 g/L。3.通过对对照组、饥饿胁迫9 h及UVA胁迫12 h的中蜂进行转录组测序,共获得62.20 GB Clean data,9个样品的Clean data均达到6.65 GB,Q30碱基百分比全部不低于90%,GC含量均在40%左右,9个样品Clean reads 比对上的参考基因组的比例分别为 90.99%,91.14%,92.45%,89.54%,91.32%,92.84%,93.33%,91.14%,92.24%。4.采用DEseq2(以|log2FC|≥1且Qvalue<0.05作为筛选标准)对获得的数据进行差异表达分析:中蜂在饥饿胁迫过程中共筛选出543个差异表达基因,经GO 功能注释结果显示在 Biological Process、Molecular Function、Cellular Component这三大类中分别注释出298、317、321个,其中有102个差异表达基因富集在response to stimulus、58个基因富集到signaling、11个差异表达基因富集到immune system process。对差异表达基因进行KEGG pathway分析结果显示:富集到细胞免疫和体液免疫的差异表达基因绝大部分表现为上调趋势,表明中蜂的免疫系统被显著激活。5.中蜂在UVA胁迫过程中共筛选出585个差异表达基因,经GO功能注释结果显示在 Biological Process、Molecular Function、Cellular Component 这三大类中分别注释出320、337、354个,其中response to stimulus富集了 117个差异表达基因,signaling富集了 61个差异表达基因,immune system process富集了 15个差异表达基因,antioxidant activity富集7个差异表达基因。对差异表达基因进行KEGG pathway分析结果显示:在UVA胁迫下,富集到中蜂细胞免疫的部分差异表达基因表现为上调趋势,更多的表现为下调趋势。富集到体液免疫的差异表达基因绝大部分表现为上调趋势,表明中蜂的细胞免疫通路被部分激活,而体液免疫被显著激活。根据GO和KEGG结果显示,饥饿及UVA胁迫下中蜂的免疫应答反应被显著激活,适应环境的变化。6.两组处理组分别选取8个免疫基因,利用qRT-PCR对其进行验证,结果与RNA-seq一致,证明转录组数据可靠。