药学学报  2018, Vol. 53 Issue (4): 574-584   PDF    
基于统计建模的壮骨关节丸诱发特异质肝损伤相关易感性细胞因子分析
唐进法1,2, 王晓艳1,2, 杨伟3, 李伟霞1,2, 李宇辉2, 柏兆方4, 王伽伯4, 肖小河4     
1. 河南中医药大学第一附属医院, 河南 郑州 450000;
2. 河南中医药大学, 河南 郑州 450008;
3. 中国中医科学院中医临床基础医学研究所, 北京 100700;
4. 解放军302医院全军中医药研究所, 北京 100039
摘要: 利用数据建模的方式筛选壮骨关节丸(Zhuangguguanjie wan,ZGW)特异质肝损伤相关易感细胞因子。以脂多糖(lipopolysaccharide,LPS)+ZGW大鼠肝脏组织中27种细胞因子水平作为数据源,以肝功能指标谷丙转氨酶(alanine aminotransferase,ALT)活力值高低作为肝损伤结局的评价指标,采用Logistic回归法、分类树及随机森林法、LASSO Logistics回归法、规则发现算法与LASSO相结合的变量筛选方法,筛查ZGW特异质肝损伤相关易感细胞因子,并在THP1细胞上进行再评价。结果表明,初步筛选分析得到ZGW特异质肝损伤密切相关的易感细胞因子组合为白介素-1β(interleukin-1β,IL-1β)、白介素-18(interleukin-18,IL-18)和表皮细胞生长因子(epidermal growth factor,EGF),ZGW醇提物(ethanolic extract of ZGW,EtZ)联合IL-1β或IL-18协同增强THP1细胞肿瘤坏死因子(tumor necrosis factor-α,TNF-α)分泌,EtZ联合IL-1β协同增强THP1细胞白介素-6(interleukin-6,IL-6)的分泌,但EtZ联合EGF协同抑制THP1细胞IL-6的分泌。上述结果表明,能表征ZGW特异质肝损伤相关易感性细胞因子为IL-1β和IL-18,为临床ZGW特异质肝损伤患者的筛查提供依据,为ZGW临床安全用药及风险防控提供了新的实验证据。
关键词: 壮骨关节丸     特异质肝损伤     细胞因子     Logistic回归     随机森林法     LASSO回归法    
Cytokine analysis of Zhuangguguanjie wan-induced idiosyncratic liver injury based on mathematical modeling
TANG Jin-fa1,2, WANG Xiao-yan1,2, YANG Wei3, LI Wei-xia1,2, LI Yu-hui2, BAI Zhao-fang4, WANG Jia-bo4, XIAO Xiao-he4     
1. The First Affiliated Hospital of Henan University of Traditional Chinese Medicine, Zhengzhou 450000, China;
2. Henan University of Chinese Medicine, Zhengzhou 450008, China;
3. Institute of Basic Research in Clinical Medicine, China Academy of Chinese Medical Sciences, Beijing 100700, China;
4. China Military Institute of Chinese Medicine, 302 Military Hospital, Beijing 100039, China
Abstract: In this study, we used a mathematic-based modeling system to screen the cytokines that are sensitive to Zhuangguguanjie wan (ZGW)-induced idiosyncratic liver injury. The values of 27 cytokines were used as the data source in rat liver of lipopolysaccharide (LPS) + ZGW group. The alanine aminotransferase (ALT) activity value of liver function indexes was used as the outcome evaluation index of liver injury. Cytokines of ZGW-induced idiosyncratic liver injury were screened using Logistic regression, random forest method, LASSO Logistics regression and method of combining rule discovery algorithm with LASSO, and cytokines filtered out were revalued in THP1 macrophage. Susceptible cytokine combinations:interleukin-1β (IL-1β), epidermal growth factor (EGF) and interleukin-18 (IL-18) closely related to ZGW-induced idiosyncratic liver injury were obtained after preliminary screening analysis. The result of revalued in THP1 showed that the ethanolic extract of ZGW (EtZ) combined with IL-1β or IL-18 synergistically enhanced tumor necrosis factor-α (TNF-α) secretion in THP1 macrophage, and EtZ combined with IL-1β significantly enhanced interleukin-6 (IL-6) secretion in THP1 macrophage, but EtZ combined with EGF markedly inhibited IL-6 secretion in THP1 macrophage. The results suggest that the sensitive cytokines that can be characterized in the ZGW-induced idiosyncratic liver injury are IL-1β and IL-18, which provides a basis for screening the ZGW-induced idiosyncratic liver injury patients, and a new experimental evidence for clinical safety medication and risk prevention of ZGW.
Key words: Zhuangguguanjie wan     idiosyncratic hepatotoxicity     cytokine     Logistic regression     Random forest method     LASSO regression    

壮骨关节丸(Zhuangguguanjie wan, ZGW)由熟地、补骨脂、淫羊藿、狗脊、独活、续断、鸡血藤、桑寄生、木香、骨碎补、乳香和没药共12味中药组成。用于肝肾不足、血淤气滞、脉络闭阻所致的骨性关节炎和腰肌劳损。上市以来, 虽然临床疗效较好, 但国家药品不良反应监测中心分别于2001年和2008年报道了ZGW的肝损伤现象。药物性肝损伤分为固有型肝损伤和特异质肝损伤, 而特异质肝损伤又可分为免疫特异质和代谢特异质[1, 2]。课题组前期基于中药特异质肝损伤“免疫应激三因致毒”机制的观点[3, 4], 利用脂多糖(lipopolysaccharide, LPS)介导的药物特异质肝损伤模型阐明了ZGW致肝损伤与何首乌诱导的特异质肝损伤一致[5, 6], 为免疫特异质肝损伤[7]。近年来, 诸多研究发现许多细胞因子介导或参与免疫性肝损伤过程, 例如, 在刀豆蛋白A (concanavalin A, Con A)诱导的急性免疫性肝损伤中, 则主要是由Th1类细胞、巨噬细胞和Th2类细胞分泌的炎性细胞因子所造成的[8], 肝损伤小鼠血清中白介素-2 (interleukin-2, IL-2)、白介素-4 (interleukin-4, IL-4)、白介素-6 (interleukin-6, IL-6)、白介素-10 (interleukin-10, IL-10)、干扰素-γ (interferon-γ, IFNγ)和肿瘤坏死因子(tumor necrosis factor, TNF-α)含量较正常对照组明显升高[9, 10]。但特异质肝损伤相关细胞因子鲜有人报道, 尤其是ZGW特异质肝损伤相关细胞因子。本文旨在通过统计建模的方法筛选ZGW特异质肝损伤相关易感细胞因子, 并在THP1巨噬细胞上进行再评价, 以期为ZGW的风险评估及合理用药提供参考。

材料与方法

动物   SD 大鼠, SPF 级, 体重180~190 g, 雄性, 购自军事医学科学院实验动物中心, 合格证号为SCXK-(军) 2012-0004。

细胞株   THP1细胞(人白血病单核细胞), 来自解放军302医院全军中医药研究所。

药物与试剂   ZGW (批号: 1512007S, 华润三九医药股份有限公司); LPS (批号: 046M4045V)、豆蔻酸佛波乙酯(phorbol 12-myristate 13-acetate, PMA, 批号: MKCC8966) (Sigma公司); 水合氯醛(批号: 20150923, 宇龙海藻有限公司); Protease Inhibitor Cocktail I (20-201)、BCA Protein Assay Kit (71285-3CN)、MILLIPLEX MAP Lysis buffer for Multiplexing (43-040)、MILLIPLEX Catalog ID (RECYTMAG-65K-27.Rat) [默克化工技术(上海)有限公司]; 表皮细胞生长因子(epidermal growth factor, EGF, 批号: 1116AFC05)、白介素-1β (interleukin-1β, IL-1β, 批号: 0103B95R1) (Peprotech公司); 白介素-18 (interleukin-18, IL-18, 批号: EHK0317012, R & D Systems公司); IL-6、TNF-α酶联免疫试剂盒(批号: 11/2017, 上海酶联生物科技有限公司)。

仪器  FW-500高速万能粉碎机(北京科伟永兴仪器有限公司); ME204电子天平(梅特勒−托利多仪器有限公司); 迈瑞BS-300全自动生化分析仪(上海富众生物科学有限公司); Neofuge1600R低温台式高速离心机(上海力申科学仪器有限公司); Luminex 200蛋白芯片多功能液相芯片分析仪[美国默克化工技术(上海)有限公司]; Milliplex磁珠检测系统、Milliplex多功能液相芯片分析平台、Multiskan FC酶标仪(美国Thermo Scientific公司)。

药物制备  ①动物实验所需ZGW制备方法:将ZGW用粉碎机粉碎, 过筛(100目)备用, 实验时用0.5%羧甲基纤维素钠(carboxymethylcellulose sodium, CMC)配制为混悬液, 备用。②细胞实验所需ZGW醇提物(EtZ)制备方法:取ZGW 200 g, 粉碎, 用5倍量75%乙醇浸泡0.5 h, 回流提取2次, 合并提取液, 回收乙醇, 浸膏进行冷冻干燥后, 称重, 得到ZGW醇提物粉末(EtZ) 56.3 g, 即3.55 g·g−1 ZGW成方(每克提取物相当于3.55克ZGW原药), 实验时用细胞培养基溶解即可。

动物模型制备及给药  实验前将SD大鼠禁食12 h, 实验分为2组, 正常对照组和ZGW (3.8 g·kg−1) + LPS (2.8 mg·kg−1)组。对照组灌胃0.5% CMC, ZGW + LPS组灌胃ZGW (0.5% CMC配制的ZGW混悬液), 2 h后对照组尾静脉注射生理盐水, ZGW + LPS组尾静脉注射LPS, 10 h后采用10%水合氯醛(3 mL·kg−1)麻醉大鼠, 腹主动脉取血, 低温离心(4 ℃、3 500 r·min−1) 10 min, 取上清液。取血后, 立即剪取大鼠肝脏组织, 放入−80 ℃冰箱冻存, 备用。

肝功能指标检测  利用全自动生化分析仪按照使用说明书测定血清谷丙转氨酶(ALT)活性。

细胞因子检测  将肝组织从−80 ℃冰箱中取出, 冰上操作, 剪取一定重量肝组织, 加入含蛋白抑制剂的裂解液(9 mL:1 g), 利用高速粉碎机粉碎为匀浆液, 4 ℃, 4 500 ×g, 离心20 min, 取上清液。按照BCA试剂盒说明书进行蛋白定量并进行蛋白调平, 采用Milliplex蛋白芯片Kit对大鼠肝组织裂解液中的嗜酸性粒细胞趋化因子(eotaxin)、趋化因子-10 (inter feron-inducible protein-10, IP-10)、单核细胞趋化蛋白1 (monocyte chemotactic protein 1, MCP-1)、调节活化正常T细胞表达和分泌的趋化因子(regulated upon activation normal T cell expressed and secreted factor, RANTES)、嗜中性趋化因子(GRO/KC)、CX3CL1趋化因子(fractalkine)、脂多糖诱导的cxc趋化因子(LIX)、巨噬细胞炎性蛋白1α (macrophage inflame matory protein 1α, MIP-lα)、MIP-2、肿瘤坏死因子-α (tumor necrosis factor-α, TNF-α)、粒细胞集落刺激因子(granulocyte colony stimulating factor, G-CSF)、粒细胞−巨噬细胞集落刺激因子(granulocyte-macro phage colony stimulating factor, GM-CSF)、表皮细胞生长因子(epidermal growth factor, EGF)、血管内皮生长因子(vascular endothelial growth factor, VEGF)、白介素-1α ((interleukin-1α, IL-1α)、IL-1β、IL-2、IL-4、IL-5、IL-6、IL-10、IL-12 (p70)、IL-13、IL-17A、IL-18、γ-干扰素(interferon-γ, IFNγ)和瘦素(leptin)等27种细胞因子进行检测与筛选。

易感细胞因子再评价  THP1细胞培养于含10%胎牛血清、1%青霉素和链霉素的1640培养液, 置于37 ℃、5% CO2培养箱中。取对数生长期的THP1细胞制成悬液, 以每孔1×104个细胞接种于96孔培养板, 实验前将THP1细胞先用100 ng·mL−1 PMA处理72 h, 诱导其为巨噬细胞。EtZ浓度取无毒剂量, IL-1β[11]、IL-18[12]和EGF[13]参考文献浓度设置3个剂量, 每孔加入200 µL不同质量浓度的药物溶液, 分别为EtZ (2.5 µg·mL−1) + IL-18 (1/10/100 ng·mL−1), EtZ (25 µg·mL−1) + IL-18 (1/10/100 ng·mL−1), EtZ (250 µg·mL−1) + IL-18 (1/10/100 ng·mL−1); EtZ (2.5 µg·mL−1) + EGF (1/10/200 ng·mL−1), EtZ (25 µg·mL−1) + EGF (1/10/200 ng·mL−1), EtZ (250 µg·mL−1) + EGF (1/10/200 ng·mL−1); EtZ (2.5 µg·mL−1) + IL-1β (0.1/1/2 ng·mL−1), EtZ (25 µg·mL−1) + IL-1β (0.1/1/2 ng·mL−1), EtZ (250 µg·mL−1) + IL-1β (0.1/1/2 ng·mL−1); 对照组EtZ (2.5/25/250 µg·mL−1), IL-18 (1/10/100 ng·mL−1), EGF (1/10/200 ng·mL−1), IL-1β (0.1/1/2 ng·mL−1)。每组设3个复孔, 培养24 h后收集细胞上清液, 采用ELISA方法检测细胞上清液中TNF-α、IL-6的含量。

统计学分析

建模数据分析方法运用R3.2.4软件, 数据中主要包括计量资料, 根据它们的分布特征选用恰当的检验方法。每个变量的具体分布特征见如下分析。组间差异比较采用3种方式及检验方法: ①正态分布连续变量用均值(标准差)描述, 两组比较用t检验, 多组比较用方差分析(ANOVA), 成组多重比较采用Tukey HSD检验法(Tukey honestly significant differences testing); ②非正态分布连续变量或等级(有序)变量用中位数(下四分位数; 上四分位数)描述, 组间比较用非参数Kruskall-Wallis秩和检验; 成组多重比较采用Benjamini-Hochberg检验法; ③如果无法确定连续变量是否正态分布, 则利用Shapiro-Wilk检验是否服从正态分布, 取显著水平为0.01, 然后再根据①和②方法进行统计学检验; ④模型主要采用多元Logistic回归、分类树及随机森林、LASSO Logistics回归法、规则发现算法与LASSO相结合的变量筛选方法[14]

再评价数据分析方法采用SPSS 17.0统计软件进行分析, 剂量资料多组间比较采用单因素方差分析(one-way ANOVA), 两组间比较用t检验, P < 0.05表示具有显著性差异。

结果 1 ZGW+LPS组各原始指标在肝损伤结局之间的差异

肝损伤动物的各指标以ALT高于2倍正常值为肝损伤, 取值为1, 低于2倍正常值为未发生肝损伤, 取值为0;以是否肝损伤为结局变量, 对与肝损伤可能相关的肝细胞因子进行分析。在ZGW+LPS组中, 对各原始指标在肝损伤结局之间的差异, 采用“统计学分析”中的第①种方法进行统计学检验, 以均值(标准差)的形式表示, 结果见表 1

Table 1 The comparison 1 of original indexes between liver injury groups.  x± s.   0 represent samples in the Zhuangguguanjie wan (ZGW)+lipopolysaccharide (LPS) group with alanine aminotransferase (ALT) activity less than 2 times when compared with that of the normal group. 1 represent samples in the ZGW+LPS group with ALT activity more than 2 times when compared with that of the normal group. IP-10: Interferon-inducible protein-10; MCP-1: Monocyte chemotactic protein 1; RANTES: Regulated upon activation normal T cell expressed and secreted factor; MIP-lα: Macrophage inflammatory protein-lα; MIP-2: Macrophage inflammatory protein-2; TNF-α: Tumor necrosis factor-α; G-CSF: Granulocyte colony stimulating factor; GM-CSF: Granulocyte-macrophage colony stimulating factor; EGF: Epidermal growth factor; VEGF: Vascular endothelial growth factor; IL-1α: Interleukin-1α; IFNγ: Interferon-γ

用“统计学分析”中的第③种方法进行统计学检验, 以均值(标准差)或中位数(下四分位数, 上四分位数)的形式表示, 结果见表 2

Table 2 The comparison 2 of original indexes between liver injury groups

在ZGW+LPS组中, 对各指标标度转换后在肝损伤结局之间的差异, 采用“统计学分析”中的第①种方法进行统计检验, 以均值(标准差)的形式表示, 结果见表 3

Table 3 The comparison 1 of converted data between liver injury groups.   x±s

采用“统计学分析”中的第③种方法进行统计学检验, 以均值(标准差)或中位数(下四分位数, 上四分位数)的形式表示, 结果见表 4

Table 4 The comparison 2 of converted data between liver injury groups
2 结局间各指标差异的可视化分析

根据以上相关系数及结局间比较的检验结果, 利用散点图、叠加密度图和箱线图(图 1~3), 对可能与肝损伤结局有较显著相关性的原始指标及标度转化后进行可视化分析。

Figure 1 Scatter plot of significant correlation cytokines from raw data without (A) or with (B) conversion. The red dots represent samples in the ZGW+LPS group with ALT activity less than 2 times when compared with that of the normal group. The blue square represent samples in the ZGW+LPS group with ALT activity more than 2 times when compared with that of the normal group

Figure 2 Overlap density map of significant correlation cytokines from raw data without (A) or with (B) conversion. Red solid line represent samples in the ZGW+LPS group with ALT activity less than 2 times when compared with that of the normal group. Blue dotted line represent samples in the ZGW+LPS group with ALT activity more than 2 times when compared with that of the normal group

Figure 3 Box chart of significant correlation cytokines from raw data without (A) or with (B) conversion
3 细胞因子筛选 3.1 Logistic回归法

假定一组协变量X, 对p(X) = P(Y = 1|X)进行建模, 利用链接函数Logit, 即$\log \frac{{p(X)}}{{1 -p(X)}} = \frac{{P(Y = 1|X)}}{{1 -P(Y = 1|X)}}$, 以最大似然估计来拟合此模型, 通过简单变换, $\frac{{p(X)}}{{1-p(X)}} = {e^{{\beta _0} + {\beta _1}X}}$称其为优势(odds), 且进一步由$\log (\frac{{p(X)}}{{1-p(X)}}) = {\beta _0} + {\beta _1}X$, 其中$p(X) = P(Y = 1|X) = \frac{{{e^{{\beta _0} + {\beta _1}X}}}}{{1 + {e^{{\beta _0} + {\beta _1}X}}}}$则获得Logistic回归方程: $\log \frac{{p(X)}}{{1-p(X)}} = {\beta _0} + {\beta _1}X$

以是否肝损伤为结局, 27种肝脏细胞因子为预测变量或影响因素。利用标度转换后的各指标数据进行Logistic回归模型分析。设定显著水平为0.2, 则单因素分析结果显示, 系数不为0的Wald检验P值小于0.2的指标如下: IL-1α、MIP-1α、IL-4、IL-1β、IL-2、IL-12p70、IL-18、MCP-1、IP-10和VEGF。然后再放入多因素Logistic回归模型中, 采用逐步回归法, 根据AIC评分准则, 选择最优模型筛选易感细胞因子, 结果见表 5

Table 5 Optimal model by Logistic regression method
3.2 分类树及随机森林法

采用分类树模型对此数据进行分析, 从各肝细胞因子中推理出合理树结构的表示形式。随机森林是利用bootstrap方法, 在分类树模型的基础上, 同时不断随机抽取样本, 以构建新的树, 构造一定数量的分类树后, 将得到的多棵树组合, 对模型进行综合评价的方法, 它能够大大提高预测精度。本数据采用随机森林的方法对变量进行建模。以是否肝损伤为结局, 27种肝脏细胞因子为预测变量或影响因素。将数据随机分为训练集和测试集, 比例为8:2, 训练集中的数据用来建立模型, 测试集的数据用来检验预测效果, 对训练集重复抽样10-折交叉验证, 选出最优模型和重要变量。随机森林选取算法特征重要性的排序结果见图 4

Figure 4 The importance of 27 cytokines for liver injury
3.3 LASSO Logistics的变量选择法

对27种肝脏细胞因子与肝损伤的关系进行分析, 主要应用带LASSO惩罚的Logistics回归模型。此处也对原始指标进行标度转换, 采用重复抽样10-折交叉验证技术筛选并构建模型。筛选出具有主效应的指标变量及其系数估计, 结果见表 6

Table 6 The variable and coefficient estimation of the main effective indexes
3.4 规则发现算法与LASSO相结合变量筛选

利用随机森林的规则发现优势与LASSO的变量筛选优势相结合的模型对此数据进行规则发现和变量筛选[14], 结果见表 7。由表中可以清晰看到, 规则1中的系数为正值最大值(0.863 1), 贡献度最高。

Table 7 Regular discovery and variable screening table of the data analyzed by combined model between random forest and LASSO methods
4 易感细胞因子再评价

与空白对照组相比, 单独给予EtZ和IL-1β对TNF-α含量无明显影响; 与EtZ (2.5 µg·mL−1)相比, EtZ (2.5 µg·mL−1) + IL-1β (0.1/1 ng·mL−1)可显著增加TNF-α含量; 与IL-1β (0.1 ng·mL−1)相比, EtZ (2.5 µg·mL−1) + IL-1β (0.1 ng·mL−1)可显著增加TNF-α含量; 与IL-1β (1 ng·mL−1)相比, EtZ (2.5/25 µg·mL−1) + IL-1β (1 ng·mL−1)均可显著增加TNF-α含量。结果见图 5A

Figure 5 Influence of co-treatment with EtZ and IL-1β, EGF, IL-18 on TNF-α and IL-6 in the THP1 macrophages. n = 3, x±s. **P < 0.01 vs control group; P < 0.05, △△P < 0.01 vs EtZ (2.5/25/250 µg·mL−1 EtZ); #P < 0.05, ##P < 0.01 vs cytokines (IL-1β, EGF, IL-18). A, C and E represent the contents of TNF-α after different treatment with EtZ and IL-1β, EGF, IL-18; B, D and F represent the contents of IL-6 after different treatment with EtZ and IL-1β, EGF, IL-18

与空白对照组相比, 单独给予EtZ和IL-1β对IL-6含量无明显影响; 与EtZ (2.5 µg·mL−1)相比, EtZ (2.5 µg·mL−1) + IL-1β (0.1/1/2 ng·mL−1)均可显著增加IL-6含量; 与EtZ (25 µg·mL−1)相比, EtZ (25 µg·mL−1) + IL-1β (1 ng·mL−1)可显著增加IL-6含量; 与EtZ (250 µg·mL−1)相比, EtZ (250 µg·mL−1) + IL-1β (0.1 ng·mL−1)可显著增加IL-6含量; 与IL-1β (0.1 ng·mL−1)相比, EtZ (2.5/250 µg·mL−1) + IL-1β (0.1 ng·mL−1)可显著增加IL-6含量; 与IL-1β (1 ng·mL−1)相比, EtZ (2.5/25/250 µg·mL−1) + IL-1β (1 ng·mL−1)均可显著增加IL-6含量; 与IL-1β (2 ng·mL−1)相比, EtZ (2.5 µg·mL−1) + IL-1β (2 ng·mL−1)可显著增加IL-6含量。结果见图 5B

与空白对照组相比, 单独给予EtZ和EGF对TNF-α含量无明显影响; 与EtZ相比, EtZ+EGF对TNF-α含量无明显影响; 与EGF (1 ng·mL−1)相比, EtZ (2.5/25 µg·mL−1) + EGF (1 ng·mL−1)可显著增加TNF-α含量。结果见图 5C

与空白对照组相比, 单独给予EtZ和EGF对IL-6含量无明显影响; 与EtZ相比, EtZ+EGF对IL-6含量无明显影响; 与EGF (1 ng·mL−1)相比, EtZ (2.5 µg·mL−1) + EGF (1 ng·mL−1)可显著降低IL-6含量; 与EGF (200 ng·mL−1)相比, EtZ (25/250 µg·mL−1) + EGF (200 ng·mL−1)可显著降低IL-6含量。结果见图 5D

与空白对照组相比, 单独给予IL-18对TNF-α含量无明显影响, 单独给予EtZ (25 µg·mL−1)可显著降低TNF-α含量; 与EtZ (25 µg·mL−1)相比, EtZ (25 µg·mL−1) + IL-18 (1/10/100 ng·mL−1)可显著增加TNF-α含量; 与IL-18 (1 ng·mL−1)相比, EtZ (2.5 µg·mL−1) + IL-18 (1 ng·mL−1)可显著增加TNF-α含量。结果见图 5E

与空白对照组比, 单独给予EtZ和IL-18对IL-6含量无明显影响; 与EtZ/IL-18相比, EtZ+IL-18对IL-6含量无明显影响。结果见图 5F

讨论

药物肝损伤包括固有型和特异质肝损伤两类[15, 16], 而药物特异质肝损伤又分免疫特异质和代谢特异质两种类型, 前者与机体的免疫应激状态有关[17, 18]。诸多研究表明免疫细胞分泌的细胞因子在免疫性肝损伤过程中发挥重要作用。IL-1β由单核巨噬细胞分泌的促炎细胞因子, 是IL-1的主要分泌形式, 与肝纤维化密切相关[19]。IL-2由CD4+T (Th1)细胞分泌, 是免疫调节网络中的核心物质, 其在重型肝炎患者血清中的水平显著下降[20]; IL-6即是免疫调节因子, 又是炎症介质, 参与炎症反应过程, 其在慢性肝病和重型肝炎患者血清中处于高活性状态[21]; IL-10是由Th2细胞分泌, 是重要的抗炎细胞因子和免疫抑制因子, 其水平在慢性肝炎患者的肝纤维化程度中随之降低; IL-18由活化的巨噬细胞、B细胞、Th1细胞和NK细胞等多种细胞产生的一种细胞因子[22], 可刺激和促进外周单核细胞产生IFNγ、IL-2、G-CSF和GM-CSF等细胞因子, 促进T细胞和NK细胞增殖, 并诱导Th1细胞分化成熟, 促进炎症的发展; G-CSF可抑制IL-1β、IL-12、IFNγ、IL-18及TNF-α的产生和分泌, 在体内具有较强的抗炎、抑制免疫效应[23], 健康人体内G-CSF可动员树突状淋巴系细胞(DC2), 进而促进T细胞产生IL-4和IL-10, 并促进血清肝细胞生长因子水平升高, 在自身免疫性肝病中有重要的调节作用。

本研究利用R3.2.4中软件包CompareGroups、AppliedPredictiveModeling[24]、caret、glinternet和RuleFit等对所得数据进行分析与模型构建。利用Logistics回归构建筛选模型, 其模型简单且解释性较好[25], 但只能考虑较少的变量, 且假定结局变量与预测变量之间是线性关系。进而采用分类树构建预测模型, 分类树算法对于数据假设很少, 适用面广泛, 且模型的结果非常适合用于解释变量之间的关系, 是一种常见且有效的非线性关系模型, 计算量相对较小且易解释。但此模型易过度拟合且只有随机单棵树, 缺乏伸缩性, 难于处理大数据集, 即使改进算法也降低了分类的准确性。而随机森林是一个包含多个分类树的模型[26], 其预测结果是由单棵树输出类别的众数而定, 它是利用bootstrap方法, 在分类树模型的基础上, 不断随机抽取样本, 以构建新的树, 构造一定数量的分类树后, 将得到的多棵树组合, 对模型进行综合评价的方法, 能够大大提高筛选精度。而且, 随机森林适合做多分类问题, 训练和预测速度快, 能有效估计缺失值, 有效处理大数据集, 在不做特征选择的情况下处理高维数据, 识别指标之间的相互影响及重要性程度, 还不易过度拟合。但此模型已经被证明在某些噪音较大的分类或回归问题上会过度拟合, 对于有多个不同取值的指标数据, 取值划分较多的指标会对随机森林产生更大的影响。因此, 根据27种指标的不同范围的取值, 利用随机森林先作为指标范围组合的发现算法, 对这些指标范围进行彻底搜索候选组合, 对其形成的有可能是多余且可能与结局变量不相关的组合(即所谓的过拟合问题), 再利用LASSO (稀疏)回归等变量选择方法选择最小风险组合集, 最后使用筛选的重要指标组合作为变量。

以是否肝损伤为结局, 27种细胞因子为变量或影响因素。直接利用标度转换后的各指标数据进行Logistic回归模型分析。设定显著水平为0.2, 先进行单因素分析, 纳入系数Wald检验P值小于0.2的指标放入多因素Logistic回归模型中, 采用逐步回归法, 根据AIC评分准则, 选择最优模型的系数估计结果如表 5。利用10-折交叉验证法, 筛选出最优随机森林模型, 根据各指标对模型Gini系数整体改善的贡献, 每个指标对肝损伤的重要程度如图 4所示, 可看出IL-18、IL-1α、VEGF、MIP-2、IL-1β、IL-10和MIP-1α对是否肝损伤具有很强的相关性。最后, 采用随机森林与LASSO相结合的模型对此数据进行变量筛选, 筛选出的易感细胞因子规则见表 7, 根据系数coeff = 0.863 1 > 0和importance = 74.86, 可认为此项指标取值组合是ZGW特异质肝损伤的易感细胞因子组合: IL-1β > 834.3、EGF > 0.970 0和IL-18 > 604.3。

TNF-α和IL-6合为机体的正性致炎因子, 使细胞发生趋化运动, 激活中性粒细胞释放弹性蛋白, 生成活性氧代谢物及引起组织浸润反应。已有研究表明, TNF-α和IL-6在一定程度上可以反映肝损伤的程度[27]。本研究在THP1细胞水平上针对上述方法筛选得到的IL-1β、EGF和IL-18细胞因子进行再评价, 联合3种细胞因子和ZGW的醇提物(EtZ)给予THP1细胞, 检测其细胞上清液中TNF-α和IL-6的含量。结果发现, EtZ联合IL-1β或IL-18协同增强THP1细胞TNF-α分泌, EtZ联合IL-1β协同增强THP1细胞IL-6的分泌, 但EtZ联合EGF协同抑制THP1细胞IL-6的分泌, 说明运用上述筛选与ZGW特异质肝损伤相关细胞因子的方法可靠, 筛选出的相关易感细胞因子(IL-1β和IL-18)可信, 可作为ZGW特异质肝损伤风险评估及合理用药的参考指标。

参考文献
[1] Amacher DE. The primary role of hepatic metabolism in idiosyncratic drug-induced liver injury[J]. Expert Opin Drug Metab Toxicol, 2012, 8: 335–347. DOI:10.1517/17425255.2012.658041
[2] Zhang X, Liu F, Chen X, et al. Involvement of the immune system in idiosyncratic drug reactions[J]. Drug Metab Pharmacokinet, 2011, 26: 47–59. DOI:10.2133/dmpk.DMPK-10-RV-085
[3] Wang JB, Cui HR, Bai ZF, et al. Precision medicine-oriented safety assessment strategy for traditional Chinese medicines:disease-syndrome-based toxicology[J]. Acta Pharm Sin (药学学报), 2016, 51: 1681–1688.
[4] Wang JB, LI CY, Zhu Y, et al. Integrated evidence chain-based identification of Chinese herbal medicine-induced hepa-totoxicity and rational usage:exemplification by Polygonum multiflorum (He shou wu)[J]. Chin Sci Bull (科学通报), 2016, 61: 971–980.
[5] He LZ, Meng YK, Zhang ZF, et al. Protective effect of Liuweiwuling tablets on idiosyncratic liver injury induced by Polygonum multiflorum in rats based on immunomodulation[J]. Chin Tradit Herbal Drugs (中草药), 2017, 48: 136–142.
[6] Li CY, Li XF, Tu C, et al. The idiosyncratic hepatotoxicity of Polygonum multiflorum based on endotoxin model[J]. Acta Pharm Sin (药学学报), 2015, 50: 28–33.
[7] Tang JF, Wang XY, Wen Q, et al. Idiosyncratic hepatotoxicity evaluation of Zhuangguguanjie wan mediated by immune stress[J]. Acta Pharm Sin (药学学报), 2017, 52: 1033–1040.
[8] Hei CC, Ma NF, Wang YR. The effect of Th1/Th2 cytokines on the mice with autoimmune acute liver injury[J]. J Histochem Cytochem (中国组织化学与细胞化学杂志), 2011, 20: 634–641.
[9] Li J, Wang XF. Liver-targeted expression of chemokine CXCL10 inhibits Con A-induced liver injury[J]. Chin J Immun (中国免疫学杂志), 2017, 33: 828–832.
[10] Xu JL, Zhang YC, Sheng GG, et al. Effect of Haizhu Yigan Jiawei formula on cell factor in immunity hepatic injury mice[J]. Chin J Integrat Tradit Western Med Liver Dis (中西医结合肝病杂志), 2017, 27: 168–169.
[11] Sahni A, Guo M, Sahni SK, et al. Interleukin-1β but not IL-1α binds to fibrinogen and fibrin and has enhanced activity in the bound form[J]. Blood, 2004, 104: 409–414. DOI:10.1182/blood-2004-01-0126
[12] Valente AJ, Yoshida T, Izadpanah R, et al. Interleukin-18 enhances IL-18R/Nox1 binding, and mediates TRAF3IP2-dependent smooth muscle cell migration. Inhibition by simvastatin[J]. Cell Signal, 2013, 25: 1447–1456. DOI:10.1016/j.cellsig.2013.03.007
[13] Xu Z, Zhang Y, Jiang J, et al. Epidermal growth factor induces HCCR expression via PI3K/Akt/mTOR signaling in PANC-1 pancreatic cancer cells[J]. BMC Cancer, 2010, 10: 161. DOI:10.1186/1471-2407-10-161
[14] Friedman JH, Popescu BE. Predictive learning via rule ensembles[J]. Ann Appl Stat, 2008, 2: 916–954. DOI:10.1214/07-AOAS148
[15] Zhang SY, Zhu XG, Zhang GP, et al. Study on acute and chronic toxicities of Polygonimu ltiflori Radix praeparata extracts[J]. J Toxicol (毒理学杂志), 2013, 27: 261–264.
[16] Navarro VJ, Senior JR. Drug-related hepatotoxicity[J]. New Engl J Med, 2006, 354: 731–739. DOI:10.1056/NEJMra052270
[17] Tailor A, Faulkner L, Naisbitt DJ, et al. The chemical, genetic and immunological basis of idiosyncratic drug-induced liver injury[J]. Hum Exp Toxicol, 2015, 34: 1310–1317. DOI:10.1177/0960327115606529
[18] Woodhead JL, Watkins PB, Howell BA, et al. The role of quantitative systems pharmacology modeling in the prediction and explanation of idiosyncratic drug-induced liver injury[J]. Drug Metab Pharmacokinet, 2017, 32: 40–45. DOI:10.1016/j.dmpk.2016.11.008
[19] Kamari Y, Shaish A, Vax E, et al. Lack of Interleukin-1α or interleukin-1β inhibits transformation of steatosis to steato-hepatitis and liver fibrosis in hypercholesterolemic mice[J]. J Hepatol, 2011, 55: 1086–1094. DOI:10.1016/j.jhep.2011.01.048
[20] Wang M, Shi Q, Zhang R, et al. Herbal compound "Jiedu Huayu" reduces liver injury in rats via regulation of IL-2, TLR4, and PCNA expression levels[J]. Evid Based Complement Alternat Med, 2017, 2017: 9819350.
[21] Zhao D, Wu T, Yi D, et al. Dietary supplementation with Lactobacillus casei alleviates lipopolysaccharide-induced liver injury in a porcine model[J]. Int J Mol Sci, 2017, 18: 2535. DOI:10.3390/ijms18122535
[22] Denton AE, Doherty PC, Turner SJ, et al. IL-18, but not IL-12, is required for optimal cytokine production by influenza virus-specific CD8+ T cells[J]. Eur J Immunol, 2007, 37: 368–375. DOI:10.1002/(ISSN)1521-4141
[23] Rutella S, Zavala F, Danese S, et al. Granulocyte colony-stimulating factor:a novel mediator of T cell tolerance[J]. J Immunol, 2005, 175: 7085–7091. DOI:10.4049/jimmunol.175.11.7085
[24] Kuhn M, Johnson K. Applied Predictive Modeling[M]. New York: Springer, 2013: 27-35.
[25] Hyndman RJ, Athanasopoulos G. Forecasting:Principles and Practice[M]. London: Bowker-Saur, 2014: 1-10.
[26] Kuhn M. Building predictive models in R using the caret package[J]. J Stat Softw, 2008, 28: 1–26.
[27] Yeo W, Mo FKF, Chan SL, et al. Hepatitis B viral load predicts survival of HCC patients undergoing systemic chemotherapy[J]. Hepatology, 2007, 45: 1382–1389. DOI:10.1002/(ISSN)1527-3350