Traditional Chinese Medicine Transcriptome Analysis of Characteristic Distribution of Scorpion Genes in Venomous Glands
-
摘要:
目的 研究东亚钳蝎Buthus martensii Karsch转录组信息特征。 方法 以山东省临沂地区特色东亚钳蝎为研究对象,利用二代高通量测序平台Illumina HiSeq 4000 150PE分别对东亚钳蝎的前腹部、后腹部、毒液腺进行转录组测序。利用测定的转录组数据将毒液腺基因分别与前腹部和后腹部基因集进行比较,然后取两两比较之后的交集进行分析,筛选差异基因集进行GO和KO分析。 结果 GO功能分类共有60个大项,含有23个生物学途径,13个细胞组分和24个分子功能。富集最显著的前20个KEGG通路中与毒素密切相关的为钙离子信号通路、磷脂酶信号通路。 结论 通过无参转录组测序,构建了东亚钳蝎转录组序列数据库,为今后东亚钳蝎功能基因的挖掘提供了序列基础并为代谢产物的生物合成提供了数据支撑,为传统中药在临床应用乃至精准医学方面提供了新的思路。 Abstract:Objective To obtain the transcriptome dataset of Buthus martensii. Methods Taking the characteristic Buthus martensii in Linyi area of Shandong Province as the research object, the second-generation high-throughput sequencing platform Illumina HiSeq 4000 150PE was used to sequence the transcriptome of the anterior abdomen, posterior abdomen and venom gland of Buthus martensii. The venom gland genes were compared with the preabdomen and the postabdomen gene sets, and then the intersection of the two comparison was analyzed. The screening differential gene sets were analyzed by GO and KO. Results There were 60 major categories of GO functional classification, containing 23 biological pathways, 13 cellular components, and 24 molecular functions. Among the first 20 KEGG pathways that were most significantly enriched, the calcium ion signaling pathway and phospholipase signaling pathway were closely related to toxins. Conclusion In this study, we have constructed the transcriptome sequence database of B. martensii by nonparticipating transcriptome sequencing, which provids the sequence basis for the future research on the gene mining function of Buthus martensii and the research basis for the biosynthesis of metabolites. -
Key words:
- Buthus martensii /
- Venom gland /
- Preabdomen /
- Postabdomen /
- Transcriptome sequencing
-
越来越多老年人需要行无痛胃肠镜检查,无痛胃肠镜检查是通过内窥镜检查消化道腔体内有无病变的一种检查,为许多消化道疾病的首选检查。在2019年全球老龄化人口占世界总人口的9.1%,据估计将在2030年达到11.7%,2050年将达到15.9%,在中国2019年老龄化人口占总人口的12.0%,将在分别在2030年与2050年达到20.7%与26.1%[1]。很多研究表明年龄是发生术后认知功能障碍(postoperative cognitive dysfunction,POCD)的危险因素,年龄≥65岁老年病人术后精神障碍发生率是年轻病人的2~10倍[2]。无痛胃肠镜检查多使用非气管插管的全凭静脉麻醉,药物多选择丙泊酚复合咪达唑仑镇静[3]。丙泊酚是一种酚类衍生物,其通过抑制γ-氨基丁酸(γ-aminobutyric acid,GABA)受体减少GABA 的释放频率,从而具有镇静和催眠作用,多数学者认为其无镇痛作用,所以需要复合镇痛药物使用[4]。咪达唑仑(Midazolam)为麻醉常用苯二氮卓类镇静药物,具有抗焦虑、镇静、催眠、抗惊厥及肌肉松弛作用[5]。右美托咪定(Dexmedetomidine)是一种高选择性α2肾上腺素能受体激动剂,通过作用于中枢和外周神经系统的α2受体产生抗焦虑、降低应激反应、稳定血流动力学、镇痛、抑制唾液腺分泌等作用,其优势在于镇静的同时无呼吸抑制的风险[6]。在心血管麻醉和围术期应用时,与咪达唑仑相比,右美托咪定无呼吸抑制的风险,并减少术后不良反应,体现出一定优势[7]。但在无痛胃肠镜检查时与丙泊酚配合使用,右美托咪定麻醉效果及最佳剂量尚无定论。本研究旨在观察不同剂量丙泊酚复合右美托咪定与丙泊酚复合咪达唑仑对老年人无痛胃肠镜的麻醉效果及不良反应的影响,探索右美托咪定最佳剂量。
1. 资料与方法
1.1 一般资料
本研究经昆明医科大学第一附属医院伦理委员会批准,患者本人与家属同意并签署知情同意书。前瞻性招募2018年9月至2019年4月行无痛胃肠镜联合检查的患者340例,年龄60~75岁,BMI值18.9~27.3 kg/m2,ASA Ⅰ~Ⅱ级。排除标准[2-3]:重要脏器失代偿,对计划用药过敏,或无法完成用蒙特利尔认知评估量表(MOCA)评估的患者。先用随机数字表法将前250例患者分成5组,每组50例;美托咪定0.25 μg/kg组(D1组)、右美托咪定0.5 μg/kg组(D2组)、右美托咪定0.75 μg/kg组(D3组)、咪达唑仑0.03 mg/kg组(M组),对照组(C组)。选出最优右美托咪定剂量后,用随机数字表法将后90例患者分成三组,每组30例,右美托咪定组(D组),咪达唑仑组(M2),对照组(C2)。
1.2 研究方法
1.2.1 麻醉方法
各组患者检测前常规胃肠道准备与禁饮、禁食。术前检查无明显特殊。签署无痛胃肠镜检查和麻醉知情同意书,准备监护仪,氧气,鼻导管吸氧套件,备用麻醉机,面罩,气管插管及急救药品。入室常规监测血压、心电,血氧饱和度,开放上肢外周静脉,予以鼻导管吸氧(3 L/min)。右美托咪定组(D1、D2、D3组),分别予以0.25 μg/kg、0.5 μg/kg、0.75 μg/kg右美托咪定静脉泵注(15 min泵注完毕)。咪达唑仑组(M组)0.03 mg/kg静脉推注。对照组(C组)予以静脉推注等量安慰剂(0.9%氯化钠溶液)。之后各组均静脉推注芬太尼(0.1 μg/kg)与丙泊酚(1.5~3 mg/kg),直至患者睫毛反射消失后行胃肠镜联检查。开始先行胃镜检查,结束后立即行肠镜检查。若术中发生三级体动时,追加丙泊酚0.2~0.5 mg/kg。若术中HR < 50次/min为术中心动过缓,予以阿托品处理;若MAP < 65 mmHg或MAP下降幅度大于30%为低血压,予以麻黄碱处理;若SpO2 < 90%为呼吸抑制,予以简易呼吸器辅助呼吸,必要时暂停检查,严重时予气管插管处理。全程监护患者直至Stward苏醒评分恢复至6分时,各项生命体征平稳后即可在家属陪同下离开。
1.2.2 观察指标
记录患者记录体动、呃逆、心动过缓、术中低血压、呼吸抑制等不良反应。体动分三级,一级为无体动;二级为患者稍体动不影响检查;三级为患者体动影响检查。最佳剂量右美托咪定组(D组),咪达唑仑组(M2),对照组(C2),用蒙特利尔认知评估量表(MOCA)评估术前(T0)、苏醒后5 min (T1)、30 min (T2)、1 h (T3)、2 h (T4)、6 h(T5)认知功能情况。MOCA评分 < 26分者,判定为POCD。用10分法(0分最不满意,10分最满意)评估患者满意度与消化内科医师满意度。D、M2和C2组记录吸烟史、饮酒史、糖尿病、高血压、高胆固醇史,是否有认知功能疾病家族史。
1.3 统计学处理
采用SPSS统计软件分析实验结果。计量资料服从正态分布用(
${{\bar x}}\pm s$ )描述,组间比较用方差分析,不同时间点的比较用重复测量的方差分析,两两比较用LSD及SNK检验,P < 0.05为差异有统计学意义。不服从正态分布用M(P25,P75)描述,组间比较用秩和检验;计数资料用n(%)描述,组间比较用χ2检验,两两比较采用调整的检验水准(α = 0.05)。2. 结果
本研究招募340例患者,D1、D2、D3、C、M组各50例,D、M2、C2组各30例。各组间年龄、性别组成、BMI、检查时间,差异无统计学意义(P > 0.05),见表1。
表 1 八组患者一般情况比较[M(P25,P75),($ \bar x \pm s$ )]Table 1. Comparison of general conditions of eight groups of patients[M(P25,P75),($ \bar x \pm s$ )]组别 年龄(岁) 男/女(n) BMI(kg/m2) 操作时间(min) 术前MOCA C组 67(62.75,71) 14/36 22.12 ± 2.86 25.16 ± 2.91 − M组 67.5(64,72.25) 32/18 23.24 ± 2.91 24.36 ± 3.35 − D1组 68.5(65,71.25) 28/22 22.70 ± 3.00 25.28 ± 2.89 − D2组 67(62.75,71.25) 30/20 23.07 ± 2.80 24.70 ± 3.30 − D3组 68(65,71.25) 26/24 23.43 ± 2.67 24.60 ± 3.08 − C2组 66.53 ± 4.68 17/13 22.95 ± 2.21 28.60 ± 4.44 28.47 ± 1.14 D组 67.80 ± 4.88 18/12 22.71 ± 2.14 27.63 ± 3.71 28.43 ± 1.07 M2组 67.97 ± 4.75 16/14 22.30 ± 1.85 28.10 ± 3.51 28.63 ± 1.13 2.1 丙泊酚用量、苏醒时间和满意度比较
丙泊酚用量在C、M、D1、D2、D3 5组中差异有统计学意义(P < 0.05),其中C组用量最多(P < 0.05);苏醒时间在C、M、D1、D2、D3 5组中,差异有统计学意义(P < 0.05),其中M组最长,D3组次之(P < 0.05)。消化内科医师满意度在C、M、D1、D2、D3 5组中,差异有统计学意义(P < 0.05),其中D2组满意度最高(P < 0.05)。患者满意度各组,差异无统计学意义(P > 0.05),见表2。
表 2 丙泊酚用量、苏醒时间、患者满意度和消化内科医师满意度比较[M(P25,P75)]Table 2. Comparison of propofol dosage,recovery time,patient satisfaction and gastroenterologist satisfaction[M(P25,P75)]组别 丙泊酚用量 苏醒时间 患者满意度 消化内科医师满意度 C组 205(195,215) 4(3,5)▲# 8(7,9) 8(7,9) △ M组 180(170,191.25) * 6(4,7) 8(7,9) 8(7,9) △ D1组 190(185,195) * 4(3,4) ▲# 8(7,9) 8(7,9) △ D2组 180(170,190) * 3(2,5) ▲# 8(7,9) 10(9,10) D3组 165(160,175) * 5.5(3.75,7) 7(6,9) 8(7,10) △ 与C组比较,*P < 0.05;与D3组比较,▲P < 0.05;与D2组比较, △P < 0.05;与M组比较,#P < 0.05。 2.2 不良反应比较
呼吸抑制在C、M、D1、D2、D3 5组中差异有统计学意义(P < 0.05),其中C组最多,M组次之(P < 0.05);心动过缓在C、M、D1、D2、D3 5组中,差异有统计学意义(P < 0.05),其中D3组最多,D2组次之(P < 0.05);体动在C、M、D1、D2、D3 5组中,差异有统计学意义(P < 0.05),其中C组最多,D3组最少(P < 0.05);各组间呃逆、术中低血压,差异无统计学意义(P > 0.05),见表3。
表 3 不良反应频数表[n(%)]Table 3. Frequency of adverse reactions [n(%)]组别 体动 呃逆 心动过缓 术中低血压 呼吸抑制 C组 38(76.0) ▲ 15(30.0) 2(4.0) ▲△ 15(30.0) 22(44.0) M组 27(54.0) ▲* 14(28.0) 3(6.0) ▲△ 9(18.0) 16(32.0) *△ D1组 29(58.0) ▲ 15(30.0) 5(10.0) ▲△ 8(16.0) 12(24.0) * D2组 26(52.0) ▲* 12(24.0) 12(24.0) ▲ 11(22.0) 7(14.0) * D3组 16(32.0) 7(14.0) 22(44.0) 18(36.0) 11(22.0) * 与C组比较,* P < 0.05;与D2组比较,△P < 0.05;与D3组比较,▲P < 0.05。 2.3 MOCA评分比较
由此得出,丙泊酚复合0.5 ug/kg右美托咪定为最优丙泊酚剂量。
MOCA评分在各组T0、T1、T2、T3、T4、T5时间点中,差异有统计学意义(P < 0.05),其中T1、T2、T3 MOCA评分低于T0(P < 0.05),且随着苏醒时间的增加,患者MOCA评分提升,在T4时评分基本恢复到术前水平,T5高于T4这可能与评分表重复测量导致患者学习有关。组间比较T1时刻(苏醒5 min),D组与C组患者评分高于M组 (P < 0.05),见表4。
表 4 C2、M2和D组各时刻MOCA评分表(${{\bar x}} \pm s$ )Table 4. MOCA score table at each moment in groups C2,M2 and D(${{\bar x}} \pm s$ )组别 T0 T1 T2 T3 T4 T5 C2组 28.47 ± 1.10 25.67 ± 2.11#▲ 27.17 ± 1.80# 27.43 ± 1.70# 28.40 ± 1.33 28.57 ± 1.28△ M2组 28.63 ± 1.13 24.60 ± 2.13# 26.30 ± 2.17# 27.20 ± 1.69# 28.10 ± 1.47 28.77 ± 1.28△ D组 28.43 ± 1.07 25.80 ± 1.92#▲ 26.87 ± 1.61# 27.37 ± 1.47# 28.27 ± 1.36 28.57 ± 1.22△ 与TO组比较,# P < 0.05;与T4组比较,△P < 0.05;与M2组比较,▲P < 0.05。 2.4 POCD发生率比较及危险因素分析
M2组POCD发生率高于D组(P < 0.05),见表5,为了探老年人POCD危险因素,对老年人是否发生POCD行二元Logistic回归分析,自变量包括:性别、BMI、吸烟史、饮酒史、糖尿病、高血压、高胆固醇史,是否有家族史、低血压、呼吸抑制。得出低血压及呼吸抑制是POCD发生的危险因素,术中发生低血压后发生POCD的概率增加6.930倍,且术中发生呼吸抑制后发生POCD的概率增加3.463倍,见表6。
表 5 C2、M2和D组因素与POCD阳性频数表[n(%)]Table 5. C2,M2 and D group factors and POCD positive frequency table [n(%)]组别 吸烟史 饮酒史 糖尿病史 高血压史 高胆固醇血症 家族史 术中低血压 呼吸抑制 POCD C2组 13(43.33) 12(40.00) 13(43.33) 13(43.33) 14(46.67) 3(10.00) 9(30.00) 13(43.33) 12(40.00) D组 15(50.00) 10(33.33) 10(3333) 15(50.00) 12(40.00) 2(6.67) 7(23.33) 5(16.67) 10(33.33) ▲ M2组 13(43.33) 10(33.33) 12(40.00) 14(46.67) 16(53.33) 2(6.67) 7(23.33) 8(26.67) 20(66.67) 与M2组比较,▲P < 0.05。 表 6 POCD危险因素logistic回归分析结果Table 6. Logistic regression analysis results of POCD risk factors应变量 变量 B P OR(95%CI) POCD 低血压 1.936 0.001 6.930(2.290,20.978) POCD 呼吸抑制 1.242 0.018 3.463(1.240,9.672) 3. 讨论
目前无痛胃肠镜已是临床上的常用检查手段,人口老龄化及生活水平的提高使得无痛胃肠镜检查的需求更为突出,这就给临床麻醉工作提出来更高的要求和更多值得探究的问题。
右美托咪定是一种α2-肾上腺素受体激动剂,可以抑制交感神经及减少术中术后并发症的发生率,有辅助麻醉的作用[7]。右美托咪定与咪达唑仑都具有协同镇静的效果,虽然各自作用位点不同,但都能减少丙泊酚用量,达到镇静效果[8-9]。本研究发现,在无痛胃肠镜检查中,丙泊酚用量在复合右美托咪定和咪达唑仑时减少。右美托咪定的镇静、镇痛效果存在剂量依赖性。本研究发现随着右美托咪定剂量的增加,苏醒时间延长,这可能与右美托咪定的剂量依赖性相关。我们发现心动过缓的发生率随着右美托咪定的剂量增加而增加,这与右美托咪定抑制交感神经活性有关,且存在一定的剂量依赖性,与其他研究相符[10]。右美托咪较咪达唑仑对呼吸功能的抑制作用小,有助于患者在麻醉期间保留自主呼吸,Inatomi等人的研究也得出相同的结果[11]。C组听体动反应最多可能是因丙泊酚较短的清除半衰期,术中体动发生较其他组多。这使丙泊酚剂量加大,丙泊酚的不良反应增加。
POCD是指患者在接受麻醉手术后,出现认知功能下降,尤其是记忆力和执行功能的下降,常表现为患者的焦虑、精神异常、人格行为改变[12]。老年人中枢神经发生退行性改变,在接受静脉麻醉后,PCOD的发生率增加[13]。本研究发现,术后早期MOCA评分降低,发生了认知功能障碍,在苏醒2 h后恢复到术前。说明丙泊酚麻醉或复合咪达唑仑或右美托咪定,术后2 h后的认知功能可较好恢复,三者都是安全的给药方式,这与目前一些研究结论是一致的[14-15]。最近的研究表明,丙泊酚会导致大鼠海马的自噬,从而导致认知功能障碍的发生[16]。丙泊酚与咪达唑仑会引起大鼠海马中晚期糖基化终产物的表达,从而使其发生认知功能障碍[17]。Yamamnaka等[18]的研究表明右美托咪定能减轻老年大鼠海马的炎症反应,从而减少认知功能障碍的发生率。陈年平等[19]学者的研究表明右美托咪定能抑制大鼠海马中白介素1β(IL-1β),肿瘤坏死因子-α(TNF-α)和核因子κB(NF-κB)表达水平,改善了大鼠海马的炎症反应,从而改善了大鼠的术后认知功能。本研究发现Dex较咪达唑仑能减少早期POCD的发生。Xu-YH等[20]的研究也揭示了,咪达唑仑的轻度镇静作用会导致高龄患者发生选择性认知障碍和长期认知障碍,与本研究相符。
本研究的发现右美托咪定0.5 μg/kg复合丙泊酚是老年人无痛胃肠镜检查中更好的给药方式,可稳定术中生命体征,提供良好的镇静效果,减少呼吸抑制的发生率,减少术后POCD的发生。
-
表 1 测序数据的统计结果
Table 1. Summary of transcriptome sequencing date and transcriptome assembly
样本 原始数据 有效数据 碱基数目(G) 错误率(%) Q20(%) Q30(%) 前腹部 46977080 46593018 6945967490 0.02 98.26 94.29 后腹部 45682534 45222188 6735593318 0.02 98.27 94.3 毒液腺 51173286 50641874 7463605213 0.02 98.26 94.36 表 2 差异表达基因的 GO 功能分类
Table 2. GO functional categories of Buthus martensii Unigenes
GO序列号 数目 功能描述 GO分类 GO:0044260 8 cellular macromolecule metabolic process BP GO:0009987 26 cellular process BP GO:0090304 1 nucleic acid metabolic process BP GO:0044237 22 cellular metabolic process BP GO:0043170 15 macromolecule metabolic process BP GO:0016055 7 Wnt signaling pathway BP GO:0034641 14 cellular nitrogen compound metabolic process BP GO:0006177 4 GMP biosynthetic process BP GO:0046037 4 GMP metabolic process BP GO:0006807 19 nitrogen compound metabolic process BP GO:0008150 71 biological_process BP GO:0044763 19 single-organism cellular process BP GO:0007275 10 multicellular organism development BP GO:1901360 13 organic cyclic compound metabolic process BP GO:0006725 13 cellular aromatic compound metabolic process BP GO:0009059 1 macromolecule biosynthetic process BP GO:0046483 13 heterocycle metabolic process BP GO:0006139 13 nucleobase-containing compound metabolic process BP GO:1901070 4 guanosine-containing compound biosynthetic process BP GO:0009225 4 nucleotide-sugar metabolic process BP GO:0034645 1 cellular macromolecule biosynthetic process BP GO:1901135 17 carbohydrate derivative metabolic process BP GO:0005975 15 carbohydrate metabolic process BP GO:0005615 25 extracellular space CC GO:0044421 29 extracellular region part CC GO:0044424 17 intracellular part CC GO:0032991 5 macromolecular complex CC GO:0044464 31 cell part CC GO:0044446 4 intracellular organelle part CC GO:0044422 4 organelle part CC GO:0044444 6 cytoplasmic part CC GO:0043226 8 organelle CC GO:0043229 8 intracellular organelle CC GO:0005578 6 proteinaceous extracellular matrix CC GO:0043234 5 protein complex CC GO:0031012 6 extracellular matrix CC GO:0030414 23 peptidase inhibitor activity MF GO:0061134 23 peptidase regulator activity MF GO:0030234 23 enzyme regulator activity MF GO:0004866 22 endopeptidase inhibitor activity MF GO:0061135 22 endopeptidase regulator activity MF GO:0004857 23 enzyme inhibitor activity MF GO:0003676 11 nucleic acid binding MF GO:0008238 12 exopeptidase activity MF GO:0098772 23 molecular function regulator MF GO:0097159 35 organic cyclic compound binding MF GO:1901363 35 heterocyclic compound binding MF GO:0004177 9 aminopeptidase activity MF GO:0004553 10 hydrolase activity,hydrolyzing O-glycosyl compounds MF GO:0016798 10 hydrolase activity,acting on glycosyl bonds MF GO:0008237 14 metallopeptidase activity MF GO:0003938 3 IMP dehydrogenase activity MF GO:0042813 3 Wnt-activated receptor activity MF GO:0009032 2 thymidine phosphorylase activity MF GO:0016154 2 pyrimidine-nucleoside phosphorylase activity MF GO:0005102 8 receptor binding MF GO:0004697 3 protein kinase C activity MF GO:0046914 25 transition metal ion binding MF GO:0008241 2 peptidyl-dipeptidase activity MF GO:0003824 102 catalytic activity MF -
[1] 李钟杰,吴文澜,王勇. 东亚钳蝎多肽资源的药理活性研究进展[J]. 中草药,2018,49(21):5179-5190. doi: 10.7501/j.issn.0253-2670.2018.21.033 [2] 毛歆,刘辉,李文鑫,等. 马氏正钳蝎的形态和结构特征[J]. 动物学杂志,2017,52(3):468-477. [3] Chippaux J P,Goyffon M. Epidemiology of scorpionism:A global appraisal[J]. Acta Tropica,2008,107(2):71-79. doi: 10.1016/j.actatropica.2008.05.021 [4] 杨彬,高文远,张艳军. 基于转录组学-蛋白质组学-多肽组学整合关联分析策略的动物药蛋白多肽类成分研究思路及方法[J]. 中草药,2019,50(5):1033-1038. doi: 10.7501/j.issn.0253-2670.2019.05.001 [5] 王巍杰,王涛. 蓖麻基因组毒素家族的生物信息学分析[J]. 基因组学与应用生物学,2017,36(8):3083-3090. [6] 张博,苏杭,任献青,等. 基于高通量转录组测序的菟丝子黄酮改善雷公藤多苷片致大鼠生殖损伤的机制研究[J]. 中国中药杂志,2019,44(16):3478-3485. [7] 吴昕怡,严媛,刘小莉. 基于高通量测序的青叶胆转录组研究[J]. 中国现代应用药学,2018,35(3):363-369. [8] 吴超,彭娟,向林,等. 基于高通量测序的铁皮石斛叶片转录组分析[J]. 分子植物育种,2016,14(12):3334-3346. [9] Grabherr M G,Haas B J,Yassour M,et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome[J]. Nature Biotechnology,2011,29(7):644-652. doi: 10.1038/nbt.1883 [10] 贾岩,张福生,肖淑贤,等. 款冬花不同发育阶段的代谢组学和比较转录组学分析[J]. 中国生物化学与分子生物学报,2017,33(6):615-623. [11] 冯陈晨,王文龙,红梅,等. 骆驼斯氏副柔线虫雌虫和雄虫比较转录组学分析[J]. 中国农业大学学报,2017,22(7):61-70. doi: 10.11841/j.issn.1007-4333.2017.07.008 [12] 江灵敏,王碧霞,孟蕾. 小叶女贞叶片转录组研究[J]. 中国现代应用药学,2018,35(11):1670-1674. [13] Ali Mortazavi,Brian A Williams,Kenneth McCue,et al. Mapping and quantifying mammalian transcriptomes by RNA-Seq[J]. Nature Methods,2008,5(7):621-628. doi: 10.1038/nmeth.1226 [14] Ye J,Fang L,Zheng H,et al. WEGO:A web tool for plotting GO annotations[J]. Nucleic Acids Res,2006,34(2):293-297. [15] Anders S,Huber W. Differential expression analysis for sequence count data[J]. Genome Biol,2010,11(10):R106. doi: 10.1186/gb-2010-11-10-r106 [16] Conesa A,Gotz S,et al. Blast2GO:A universal tool for annotation,visualization and analysis in functional genomics research[J]. Bioinformatics,2005,21(18):3674-3676. doi: 10.1093/bioinformatics/bti610 [17] 罗辉,叶华,肖世俊,等. 转录组学技术在水产动物研究中的运用[J]. 水产学报,2015,39(4):598-607. [18] 马赛箭,安超,薛文娇,等. 基于高通量测序的出芽短梗霉转录组学研究[J]. 基因组学与应用生物学,2016,35(4):931-941. [19] 蒋可人,马峥,郑航,等. 转录组与蛋白质组整合分析在生物学研究中的应用[J]. 生物技术通报,2018,34(12):50-55. [20] 展波,高媛媛,林维平,等. 中华大蟾蜍皮肤Cathelicidin家族新型抗菌肽的鉴定及其抗菌活性[J]. 中国中药杂志,2016,41(4):630-635. [21] RA Dueñas-Cuellar,Santana C,ACM Magalhães,et al. Scorpion toxins and ion channels:Potential applications in cancer therapy[J]. Toxins,2020,12(5):326-327. doi: 10.3390/toxins12050326 [22] 陈惠,辛丽丽,龚婕宁. 基于全转录组测序技术的转录组学在中医药领域的应用前景分析[J]. 环球中医药,2013,6(10):759-763. doi: 10.3969/j.issn.1674-1749.2013.10.014 [23] 张晓萌,李健春,王琼,等. 转录组测序技术在中医药领域的应用[J]. 中国现代中药,2016,18(8):1084-1087. [24] 肖炜卓,张龙,王晓冬,等. 慢性牙周炎患者龈下菌斑宏转录组学分析[J]. 中国医药导报,2019,16(21):21-23. -