中国血吸虫病防治杂志(中英文) ›› 2026, Vol. 38 ›› Issue (1): 29-37.
王鑫瑶1, 2,彭丹丹2, 3,杨颖1,张键锋2,秦志强1,杨坤2, 3,李石柱1, 4,许静1, 4*
WANG Xinyao1, 2, PENG Dandan2, 3, YANG Ying1, ZHANG Jianfeng2, QIN Zhiqiang1, YANG Kun2, 3, LI Shizhu1, 4, XU Jing1, 4*
摘要: 目的 探究极端温度和干旱胁迫对球形小泡螺生存能力的影响,为无参考基因条件下小泡螺基因组学研究提供理论依据。方法 自桑给巴尔奔巴岛Kiwani Shehia地区现场采集球形小泡螺样本,经实验室饲养、繁殖获取子代螺。选取120只来自同一代的10 周龄小泡螺作为实验螺,随机分配至高温干旱组(HD组)、常温干旱组(D组)、低温干旱组(LD组)和对照组(C组),每组30只螺。将HD组、D组、LD组小泡螺均置于底部为干燥土壤的烧杯中,随后分别置于35、26 ℃和10 ℃的气候培养箱中;C组小泡螺置于装有26 ℃去氯自来水的 500 mL培养皿中。饲养3 d后收集各组活螺,解剖分离软体组织,提取样本总 RNA并构建相应基因库,随后在Illumina HiSeq 4000 测序平台进行高通量测序。采用Trinity 软件进行从头转录组组装并选取最长的转录本作为单基因簇。使用DIAMOND软件对单基因簇进行基因功能注释,用于基因功能标注的数据库包括基因本体论(Gene Ontology,GO)、京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)、NCBI非冗余蛋白序列数据库(NCBI non⁃redundant protein sequences database,NR)、蛋白质家族数据库(Protein Family Database,Pfam)和瑞士⁃Prot蛋白序列数据库(UniProtKB/Swiss⁃Prot,Swiss⁃Prot)。分别使用topGO 软件和clusterProfiler软件对差异表达基因(differentially expressed genes,DEGs)进行GO、KEGG分析。另外选择4个相关基因进行实时荧光定量PCR(real⁃time quantitative PCR,qRT⁃PCR)验证,以确定转录组测序结果的可靠性。结果 饲养3 d后,HD组、LD组、D组和C组分别有7、20、28只和30只小泡螺存活,各组存活率分别为 23.33%(7/30)、66.67%(20/30)、93.33%(28/30)和 100.00%(30/30),差异有统计学意义([χ2] = 52.72,P < 0.001)。通过无参转录组测序装配了176 942个单基因簇,在上述5个数据库中的标注率分别为0.98%、13.49%、26.46%、12.48%和14.39%。基因表达水平分析结果显示,与C组相比,D组共有33个基因上调、 72个基因下调,HD组共有483个基因上调、815个基因下调, LD组共有245个基因上调、172个基因下调。去除各组重叠基因及不匹配基因后,共确定了11个候选基因。GO和KEGG分析结果显示,上述11个DEGs中,热休克蛋白(heat shock protein,HSP)相关DEGs有3个,分别注释为HSP12.2、HSP70、HSP20基因,在各处理组中的表达均显著上调;免疫和神经系统相关DEGs有3个,在各处理组中均显著下调,涉及神经细胞黏附分子L1样蛋白通路、纤维蛋白原结合蛋白通路和白细胞弹性蛋白酶抑制剂样蛋白通路。qRT⁃PCR检测结果显示,4个与温度及干旱胁迫相关的基因表达水平在不同处理组中的表达趋势与转录组测序结果高度一致。结论 极端温度和干旱联合胁迫下球形小泡螺存活率降低,可能与其细胞内稳态调节系统失衡有关。
中图分类号: