摘要
为了探究饥饿胁迫对日本医蛭唾液腺基因表达水平的影响,本实验通过Illumina Hiseq 2500高通量测序平台分别对饥饿30 d处理组(D30)和饥饿0 d (D0)对照组的日本医蛭的唾液腺组织进行双端测序,对获得的原始数据进行质量控制及从头组装,获得145 981个unigenes,平均长度为675 bp,N50为1 127 bp。针对CDD、KOG、COG、NR、NT、PFAM、Swissprot、TrEMBL、GO和KEGG等数据库的序列比对分析,145 981个unigenes均能得到注释。通过引入TPM(transcripts per million)来估算基因表达水平,并通过DEGseq进行基因表达差异分析。以饥饿0 d为对照组,显著差异基因筛选条件设置为Q value<0.05且差异倍数|Fold Change|>2,获得2 650个差异基因,其中667个基因显著上调,1 983个基因显著下调。选取日本医蛭唾液腺4个重要功能基因进行荧光定量PCR(qRT-PCR)验证,结果证明转录组测序分析可靠。将所有差异表达基因进行GO功能富集分析,显著富集到175条途径,其中胞质核糖体途径富集程度最为显著,核糖体途径次之,且参与这些途径的差异基因基本均下调。KEGG通路富集分析进一步证明,差异基因在核糖体通路中富集程度最为显著。此外,差异基因中4个基因被预测参与抗凝、抗血栓、抗菌、抗炎和抗肿瘤过程,这可能在各种疾病的治疗中发挥重要作用。参与核糖体通路的基因表达量显著降低,表明日本医蛭唾液腺通过降低蛋白质代谢以应对饥饿环境。该实验结果为深入研究日本医蛭唾液腺饥饿胁迫适应的分泌调控机制以及药用价值基因的发掘提供了重要的参考材料,并为其他医学蛭类研究提供参考依据。
- 单位