Effects of glycyrrhizin on apoptosis and inflammatory response of human embryonic lung fibroblasts infected with human cytomegalovirus
-
摘要:目的
探讨甘草酸对感染人巨细胞病毒(HCMV)的人胚肺成纤维细胞MRC-5增殖、凋亡和炎性反应的影响。
方法体外培养MRC-5细胞，并将其分为对照组、病毒组、更昔洛韦组(50 μg/mL)、甘草酸低剂量组(0.25 mg/mL)、甘草酸中剂量组(0.50 mg/mL)、甘草酸高剂量组(1.00 mg/mL)。接种HCMV不同时间后，分别使用CCK8试剂盒检测细胞增殖活性，使用流式细胞仪检测细胞凋亡情况，采用蛋白质印迹法(Western blot)检测细胞中活化含半胱氨酸的天冬氨酸蛋白水解酶3(Cleaved-caspase3)等蛋白表达情况，采用酶联免疫吸附试验(ELISA)检测细胞培养上清液中肿瘤坏死因子-α(TNF-α)和白细胞介素-6(IL-6)水平。
结果CCK8实验结果显示，经HCMV感染48、72 h时，病毒组增殖活性均低于对照组，更昔洛韦组、甘草酸中剂量组、甘草酸高剂量组增殖活性均高于病毒组，差异有统计学意义(P < 0.05)；甘草酸中剂量组与甘草酸高剂量组增殖活性比较，差异无统计学意义(P>0.05)，故选择甘草酸中剂量组(0.50 mg/mL)进行后续实验，并简称甘草酸组。流式细胞仪检测结果显示，经HCMV感染48、72 h时，病毒组细胞凋亡率均高于对照组，更昔洛韦组、甘草酸组则低于病毒组，差异有统计学意义(P < 0.05)。Western blot检测结果显示，经HCMV感染48 h后，病毒组细胞中Cleaved-caspase3、TNF-α、IL-6蛋白相对表达量均高于对照组，更昔洛韦组、甘草酸组则低于病毒组，差异有统计学意义(P < 0.05)。ELISA检测结果显示，经HCMV感染24、48、72 h时，病毒组MRC-5细胞TNF-α、IL-6水平均高于对照组，更昔洛韦组、甘草酸组则低于病毒组，差异有统计学意义(P < 0.05)。
结论甘草酸具有体外抗HCMV效应，可促进感染HCMV的人胚肺成纤维细胞MRC-5增殖，并抑制HCMV感染诱导的细胞凋亡和细胞炎性反应。
Abstract:ObjectiveTo investigate the effects of glycyrrhizin on the proliferation, apoptosis and inflammatory response of human embryonic lung fibroblast MRC-5 infected with human cytomegalovirus (HCMV).
MethodsMRC-5 cells were cultured in vitro and were divided into control group, virus group, ganciclovir group(50 μg/mL), low dose glycyrrhizin group(0.25 mg/mL), medium dose glycyrrhizin group(0.50 mg/mL) and high dose glycyrrhizin group(1.00 mg/mL). After different periods of virus inoculation, cell proliferation activity was detected by CCK8 assay, cell apoptosis was detected by flow cytometry, cleaved-caspase3 caspase was detected by western blot, and tumor necrosis factor α(TNF-α) and interleukin-6 (IL-6) levels in cell culture supernatant were detected by enzyme-linked immunosorbent assay (ELISA).
ResultsAt 24 and 48 h of HCMV infection, the results of CCK8 experiment showed that the proliferation activity of cells in virus group was significantly reduced compared with the control group, and was higher in the ganciclovir group, medium dose glycyrrhizin group and high dose glycyrrhizin group than that in the virus group(P < 0.05). There were no statistically significant difference in proliferative activity between medium dose glycyrrhizic acid group and high dose glycyrrhizic acid group (P>0.05), so the medium dose glycyrrhizic acid group (0.50 mg/mL) was selected for follow-up experiment and was called glycyrrhizic acid group for short. Flow cytometry showed that after 48 and 72 h of HCMV infection, the apoptosis rate of the virus group was higher than that of the control group, but the apoptosis rates in the ganciclovir group and glycyrrhizic acid group were lower than that of the virus group (P < 0.05). Western blot analysis showed that 48 h after HCMV infection, the relative expression levels of cleaved caspase3, TNF-α, and IL-6 in the virus group were higher than those in the control group, but were lower in the ganciclovir group and glycyrrhizic acid group than those in the virus group (P < 0.05). ELISA results showed that the levels of TNF-α and IL-6 in MRC-5 cells in the virus group were higher than those in the control group at 24, 48 and 72 h after HCMV infection, while those in ganciclovir group and glycyrrhizic acid group were lower than those in the virus group (P < 0.05).
ConclusionGlycyrrhizic acid has anti-HCMV effect in vitro, which could promote the proliferation of infected MRC-5 cells, inhibit cell apoptosis and inflammatory response induced by HCMV infection.
-
在人口老龄化和代谢危险因素的双重压力下，中国居民的心血管疾病风险持续增加，心肌梗死患者的病死率亦逐渐上升[1]。心肌梗死是指心肌缺血性坏死，通常是在冠状动脉病变的基础上发生冠状动脉血供减少或中断，心肌细胞因严重且持续的缺血缺氧而发生坏死[2]。心肌细胞的死亡会引发一系列免疫炎症反应，这些免疫炎症反应介导的受损心肌修复在梗死后心室重构、心功能和患者结局等方面起着关键作用[3]。急性心肌梗死的诊断主要依据临床表现、血清生物学标志物和心电图，其中临床常用的生物学标志物包括肌红蛋白、肌酸激酶-MB、心肌肌钙蛋白I和心肌肌钙蛋白T[4-5]。然而，这些标志物水平在心力衰竭、肾衰竭和甲状腺疾病患者中也同样升高，故结合多组学平台探寻更多特异性的生物学标志物非常必要[6-7]。外周血单个核细胞(PBMCs)包括T细胞、B细胞、单核细胞、树突状细胞和自然杀伤细胞，在免疫应答、免疫监测和免疫治疗中发挥着重要作用[8]。本研究基于PBMCs分析心肌梗死和梗死恢复过程中潜在的关键基因和生物学改变，旨在为心肌梗死的诊断、评估、管理和预后分析提供新的见解。
1. 材料和方法
1.1 数据集获取
从GEO数据库中下载关于心肌梗死PBMCs的转录组测序芯片数据集GSE59867[9-10], 芯片注释平台为GPL6244, 芯片信息为Affymetrix Human Gene 1.0 ST Array。该数据集收集了111例ST段抬高型心肌梗死(STEMI)患者入院、出院、心肌梗死后1个月、心肌梗死后6个月时的外周血样本数据。对照组为46例无心肌梗死病史的稳定型冠心病患者。基于芯片注释平台对数据集进行基因名转换，用于后续共表达网络的构建。
1.2 加权基因共表达网络分析(WGCNA)
通过WGCNA R包筛选与心肌梗死有关的易感基因模块[11]。WGCNA是一种常用的生物信息学方法，可以识别具有相似表达模式的基因模块，分析基因模块与样本表型之间的关系，绘制基因模块中的调控网络，识别关键模块和调控基因[12-13]。首先，选取数据集中方差变异程度最高的前3 000个基因构建共表达网络; 其次，利用R函数pickSoftThreshold计算软阈值1~20; 然后，对这些基因进行层次聚类和动态树切割来划分模块，并合并相似的模块(阈值为0.25); 最后，计算模块与疾病之间的Spearman相关系数，进而确定易感基因模块。为了探讨所选模块的潜在生物学功能，通过ClusterProfiler R包对模块中的基因进行基因本体论(GO)分析和京都基因与基因组百科全书(KEGG)通路分析[14]。GO分析包括生物过程(BP)、分子功能(MF)和细胞成分(CC)。通过Ggplot2 R包进行结果展示。
1.3 基因差异分析和转录因子预测
使用Limma R包进行基因差异分析[15]。差异表达基因(DEGs)的筛选标准为矫正后P < 0.05, |log2(倍数变化)|>0.5。基于ClusterProfiler R包对DEGs进行GO分析和KEGG通路分析, P < 0.05表示显著富集。对上述易感基因模块和DEGs进行交集分析，通过ChEA3数据库对交集基因进行转录因子预测。ChEA3是一种转录因子富集分析工具，整合了多个数据库，如ENCODE、GTEx、ARCHS4等数据库[16]。由于ChEA3数据库并未对预测结果可靠度进行严格划分，参考相关研究[16-17]结论，本研究对平均排名分≤30分的共同转录因子进行结果展示。基于STRING数据库构建交集基因的蛋白相互作用网络，并将交互评分>0.4分的数据导入Cytoscape 3.8.2软件进行可视化处理[18-19]。
1.4 关键基因的确定和受试者工作特征(ROC) 曲线分析
LASSO回归分析是常用的筛选变量的压缩估计方法，通过构造惩罚函数获得较为精炼的模型，使得其压缩一些回归系数，同时设定一些回归系数为0，因此保留了子集收缩的优点，是一种处理具有复共线性数据的有偏估计，常被用于协变量筛选[20]。基于glmnet R包对交集基因中存在蛋白相互作用的基因进行LASSO回归分析，进一步筛选关键基因[21]。采用10倍交叉验证法选择惩罚项(λ), 并选择最小二项偏差在1个标准误差内的最简单模型的λ值(lambda.1se)。通过pROC R包对关键基因进行ROC曲线分析，计算曲线下面积(AUC), 评估关键基因对心肌梗死的诊断价值[22]。
1.5 基因集富集分析(GSEA)和外部数据集验证
根据每个关键基因的表达量将数据集划分为高表达组和低表达组，通过GSEA探讨高表达组和低表达组之间潜在的生物学功能变化，将特征基因集作为预先定义的基因集。作为一种常见的生物信息学方法, GSEA可评估一个预先定义的基因集的基因在与表型相关度排序的基因表中的分布趋势，从而判断其对表型的影响[23]。从GEO数据库中选取另一个数据集GSE123342, 将其用于验证关键基因的表达。GSE123342数据集包括急性心肌梗死(65例)、梗死后30 d(64例)、梗死后1年(37例)、稳定型冠心病(22例)样本数据以及4例技术重复样本数据。
2. 结果
2.1 易感基因模块的确定
根据无标度网络原理，选取10作为共表达网络的软阈值，见图 1A。通过层次聚类和动态树切割，将前3 000个基因划分为15个基因模块，见图 1B。剔除聚类失败的灰色模块，棕色模块的总体表达与心肌梗死的相关性最高(r=0.4), 且与出院至心肌梗死后6个月的演变过程呈负相关(r=-0.33), 见图 1C。棕色模块中的共表达基因伴随着心肌梗死的发生呈现高表达趋势，而在心肌梗死恢复过程中则呈现低表达趋势。由此提示，棕色模块可能是易感基因模块，其中包含243个基因，即ABHD5、ACER3、ACP3、ACSL1、ADAM17、ADAM9、ADM、AHR、ALCAM、ALDH2、ALPK1、ANKRD50、ANO5、AQP9、ARHGAP24、ARHGAP29、ARRDC4、ASPH、ATP6V0A1、ATP6V1A、B3GNT5、BACH1、BLVRB、BST1、C3AR1、C5、C9orf72、CALCRL、CAPG、CAPZA2、CARD6、CASP1、CASP5、CCDC88A、CCR2、CD14、CD163、CD1D、CD33、CD36、CD63、CLEC1A、CLEC4A、CLEC4D、CLEC4E、CLEC7A、CLMN、CNTLN、CPED1、CPM、CPNE8、CPVL、CR1、CR1L、CREG1、CRISPLD2、CSTA、CTNNA1、CTSH、CYBRD1、CYP1B1、CYP1B1-AS1、DACH1、DDIAS、DOCK4、DOCK5、DRAM1、DSC2、DSE、DUSP6、EDNRB、ENTPD1、ERLIN1、EVI5、F5、FAM114A1、FAM13A、FAM151B、FAM20A、FAM198B、FAR2、FBN2、FBP1、FCGR2A、FGD4、FGD6、FLT3、FLVCR2、FMO5、FPR1、FPR2、FRRS1、FUCA1、FUCA2、GAPT、GAS2L3、GCA、GIMAP8、GLT1D1、GPAT3、GPR141、HAL、HGF、HMGB2、HNMT、HORMAD1、HP、HPSE、IDH1、IFNGR1、IL13RA1、IL15、IL18、IMPA2、IRAK3、JAK2、KCNJ15、KCNJ2、KIF13A、KLHL8、KYNU、LHFPL2、LILRA5、LIN7A、LIPN、LMNB1、LPCAT2、LRMDA、LRRK2、LTA4H、LY86、LYVE1、MAP3K20、MCEMP1、ME1、METTL7A、METTL7B、MFSD1、MGST1、MGST2、MILR1、MNDA、 MOSPD2、MS4A4A、MSRB1、MTARC1、MYOF、NAAA、NAIP、NETO2、NLN、NLRC4、NPL、OLFML2B、P2RX7、P2RY13、PGD、PIP4P2、PLA2G4A、PLBD1、PLD1、PLIN2、PLSCR1、PPARG、PPT1、PRRG4、PSTPIP2、PYGL、QPCT、RAB39A、RALB、RBP7、RGL1、RNASE2、RNASE6、RNF141、RNF217、RRAGD、S100A12、S100A9、S100Z、SAT2、SEMA3C、SEPTIN10、SERPINB1、SERPINB2、SESTD1、SGMS2、SH3PXD2B、SHTN1、SIPA1L2、SIRPD、SLC15A2、SLC1A3、SLC22A15、SLC22A4、SLC26A8、SLC31A1、SLC36A4、SLC7A7、SLC8A1、SLITRK4、SMAD1、SMPDL3A、SNX10、SOD2、SORT1、SPATA6、ST3GAL6、ST6GALNAC3、STEAP4、SULT1B1、TASL、TBC1D12、TCN2、TDRD9、TFEC、TGFBI、TLR1、TLR2、TLR4、TLR5、TLR6、TLR7、TLR8、TM6SF1、TMEM144、TMEM167A、TMTC2、TNFAIP6、TNFSF13B、TPST1、TSPO、UBE2D1、UGGT2、VNN2、VNN3、WASF1、WDFY3、WLS、ZC3H12C、ZFYVE16、ZNF438。
![]()
图 1 易感基因模块的确定与功能注释A: 共表达网络的无标度拟合指数; B: 基因聚类树状图; C: 模块特征基因与临床特征相关性热图; D、E: GO分析结果、KEGG分析结果的气泡图展示。棕色模块GO分析最显著的结果分别是骨髓白细胞激活(BP)、囊泡(CC)和模式识别受体活性(MF), 见图 1D; KEGG分析结果则主要为感染、免疫炎症等相关生物学过程的改变，包括Toll样受体信号通路等，见图 1E。
2.2 基因差异分析
进一步对数据集进行基因差异分析，根据筛选标准共得到142个DEGs(图 2A), 其中上调的DEGs有77个，下调的DEGs有65个，见表 1。BP分析结果主要包括免疫反应、免疫系统过程和免疫效应过程, CC分析结果包括分泌颗粒、分泌小泡和细胞质小泡部分等, MF分析结果主要包括信号受体活性、分子传感器活性和碳水化合物结合，见图 2B。KEGG分析提示，这些DEGs主要涉及免疫、炎症和感染的相关生物学改变，见图 2C。DEGs与上述棕色模块存在35个交集基因，见图 2D。基于ChEA3数据库对35个交集基因进行共同转录因子预测，其中平均排名分≤30分的共同转录因子分别是CREB5、MTF1、NFE4、SPI1、ZNF467、NFE2、TFEC、MXD1、NR1H3、BORCS8MEF2B、NFIL3和CEBPE, 见表 2。构建蛋白互作网络，其中有18个基因存在蛋白质相互作用，见图 2E。
![]()
图 2 基因差异分析和蛋白互作网络A: 差异基因火山图; B、C: 差异基因的GO分析、KEGG分析结果; D: DEGs与棕色模块的交集基因; E: 18个交集基因的蛋白相互作用网络以及各基因的log2(倍数变化)改变量。表 1 差异基因与交集基因上调DEGs 下调DEGs 与棕色模块交集的DEGs ACSL1、ADM、AQP9、ASGR2、BLVRB、BST1、CCR1、CCR2、CD14、CD163、CES1、CES1P1、CR1、CYP1B1、CYP27A1、DSC2、DYSF、ECRP、EGR1、EGR2、FAM20A、FAM20C、FLVCR2、FMN1、FOLR3、FPR1、FPR2、FAM198B、GPR162、HLX、HP、IFITM3、IL1R2、KCNJ15、LILRA3、LILRA5、LILRB4、LRP1、MAFB、MCEMP1、MERTK、MIR21、MS4A4A、MTARC1、NFE2、NRG1、PADI2、PLB1、PLBD1、PPARG、QPCT、RNASE1、RNASE2、S100A12、S100A9、SASH1、SERPINB10、SERPINB2、SH3PXD2B、SIGLEC16、SIGLEC9、SLED1、SOCS3、ST14、STAB 1、STEAP4、TBC1D2、TCN2、TDRD9、TMEM150B、TMEM176A、TMEM176B、TNFAIP6、TREM1、TRIB1、TSPO、VSIG4 ARAP2、BIRC3、C12orf75、CD160、CD3G、CENPK、CEP78、CLC、CLEC2D、CPA3、DTHD1、EPB41L4A、FAM169A、FCRL6、FGFBP2、GBP4、GNLY、HRH4、IL18RAP、JAKMIP2、KLRA1P、KLRC1、KLRC3、KLRC4-KLRK1、KLRD1、KLRF1、KLRG1、KPNA5、MIR29A、MS4A3、MYBL1、MYO6、ODF2L、PDGFD、PRSS23、RNU5A-1、RORA、SAMD3、SCARNA9、SH2D1A、SLC4A4、SLCO4C1、SNORA16B、SNORA20、SNORA75、SNORD20、SNORD28、SNORD30、SNORD45C、SNORD50A、SNORD54、SNORD58A、SNORD59A、SNORD59B、SNORD61、SNORD81、SNORD82、SYTL2、TARP、TBC1D19、TFDP2、TGFBR3、VPS13A、ZNF675、ZNF737 ACSL1、ADM、AQP9、BLVRB、BST1、CCR2、CD14、CD163、CR1、CYP1B1、DSC2、FAM20A、FLVCR2、FPR1、FPR2、FAM198B、HP、KCNJ15、LILRA5、MCEMP1、MS4A4A、MTARC1、PLBD1、PPARG、QPCT、RNASE2、S100A12、S100A9、SERPINB2、SH3PXD2B、STEAP4、TCN2、TDRD9、TNFAIP6、TSPO 表 2 共同转录因子预测结果(平均排名分≤30分)排名 转录因子 平均得分/分 转录因子数据库计算方法与得分 重叠基因 1 CREB5 7.667 ARCHS4 Coexpression, 5分; Enrichr Queries, 14分; GTEx Coexpression, 4分 CR1、STEAP4、TNFAIP6、ACSL1、AQP9、FPR1、MCEMP1、ADM、FPR2、LILRA5、BST1、QPCT、PLBD1、CYP1B1、CD14、S100A9、CCR2 2 MTF1 8.500 ARCHS4 Coexpression, 6分; GTEx Coexpression, 11分 CR1、SERPINB2、TNFAIP6、ACSL1、AQP9、FPR1、FPR2、LILRA5、BST1、PLBD1、CYP1B1、FLVCR2、CD14、S100A9 3 NFE4 9.000 GTEx Coexpression, 9分 BST1、CR1、AQP9、FPR1、PLBD1、MCEMP1、FPR2、LILRA5、CCR2 4 SPI1 9.667 Literature ChIP-seq, 5分; ARCHS4 Coexpression, 11分; ENCODE ChIP-seq, 8分; Enrichr Queries, 11分; ReMap ChIP-seq, 5分; GTEx Coexpression, 18分 CD163、CR1、STEAP4、SERPINB2、ACSL1、AQP9、FPR1、TDRD9、MCEMP1、FPR2、RNASE2、LILRA5、MS4A4A、BST1、QPCT、PLBD1、CYP1B1、CD14、S100A9、CCR2 5 ZNF467 19.333 ARCHS4 Coexpression, 18分; Enrichr Queries, 30分; GTEx Coexpression, 10分 CR1、AQP9、FPR1、MCEMP1、ADM、FPR2、LILRA5、BST1、PLBD1、CYP1B1、BLVRB、CD14、S100A9、CCR2 6 NFE2 21.000 ARCHS4 Coexpression, 16分; ENCODE ChIP-seq, 78分; Enrichr Queries, 1分; ReMap ChIP-seq, 8分; GTEx Coexpression, 2分 CR1、SERPINB2、AQP9、FPR1、MCEMP1、ADM、FPR2、RNASE2、LILRA5、BST1、TCN2、QPCT、PLBD1、S100A12、BLVRB、TSPO、CYP1B1、PPARG、CD14、S100A9、CCR2 7 TFEC 21.667 ARCHS4 Coexpression, 21分; Enrichr Queries, 4分; GTEx Coexpression, 40分 CD163、CR1、SERPINB2、TNFAIP6、ACSL1、AQP9、FPR1、MS4A4A、BST1、TCN2、FLVCR2、PPARG、CD14、S100A9、CCR2 8 MXD1 23.000 ARCHS4 Coexpression, 4分; Enrichr Queries, 9分; GTEx Coexpression, 56分 CR1、SERPINB2、STEAP4、ACSL1、AQP9、FPR1、ADM、FPR2、LILRA5、BST1、PLBD1、CYP1B1、CD14、S100A9 9 NR1H3 26.600 Literature ChIP-seq, 10分; ARCHS4 Coexpression, 17分; Enrichr Queries, 10分; ReMap ChIP-seq, 36分; GTEx Coexpression, 60分 CD163、STEAP4、TCN2、ACSL1、FLVCR2、CYP1B1、PPARG、ADM、CD14、S100A9、MS4A4A 10 BORCS8MEF2B 29.000 ARCHS4 Coexpression, 29分 AQP9、FPR1、TSPO、FPR2、S100A9、LILRA5 11 NFIL3 30.000 ARCHS4 Coexpression, 1分; Enrichr Queries, 21分; GTEx Coexpression, 68分 CR1、TNFAIP6、ACSL1、AQP9、FPR1、ADM、FPR2、LILRA5、MS4A4A、BST1、CYP1B1、PLBD1、CD14、S100A9 12 CEBPE 30.000 ARCHS4 Coexpression, 31分; Enrichr Queries, 12分; GTEx Coexpression, 47分 BST1、FPR1、PLBD1、ADM、MCEMP1、FPR2、CD14、RNASE2、S100A9、CCR2、LILRA5 计算方法包括ARCHS4 Coexpression、ENCODE ChIP-seq、Enrichr Queries、GTEx Coexpression、Literature ChIP-seq、ReMap ChIP-seq; 平均得分为各计算方法得分之和除以计算方法数量所得之商。 2.3 关键基因的确定
对18个存在蛋白相互作用的基因进一步行LASSO回归分析，这些基因不同惩罚参数值所对应的系数见图 3A。本研究选取最小二项偏差在1个标准误差内的λ值(lambda.1se, 6个基因)，该λ值提供更精简的模型，见图 3B。ROC曲线分析显示, CD163、RNASE2、HP、FAM20A、MCEMP1和FAM198B基因表达水平对心肌梗死的发生均具有良好的诊断价值, AUC分别为0.831、0.798、0.775、0763、0.866和0.829, 见图 3C。与对照组稳定型冠心病患者相比，这些关键基因的表达水平在入院时心肌梗死患者中显著上升，并在出院时、心肌梗死后1个月、心肌梗死后6个月逐步下降，见图 3D。
![]()
图 3 LASSO回归与ROC分析A、B: 基于惩罚项的LASSO回归系数图和交叉验证图; C: 基于ROC分析评估关键基因表达在预测心肌梗死中的作用; D: 关键基因在不同时点的表达折线图(两者比较, ***P < 0.001)。2.4 GSEA结果和验证
GSEA结果提示，这些关键基因涉及的生物学改变主要与糖脂代谢、活性氧、免疫炎症等有关，见图 4。这些基因在外部数据集中大多也存在差异表达。基于外部验证数据集GSE123342, 本研究同样发现，与稳定型冠心病相比, HP、FAM198B、CD163、FAM20A、MCEMP1表达水平在急性心肌梗死发生时显著上升，而RNASE2仅表现出上升趋势，见图 5。
3. 讨论
本研究基于PBMCs转录组学的变化探讨心肌梗死及梗死恢复过程中可能存在的生物学改变和潜在标志物, WGCNA和DEGs分析结果显示，心肌梗死过程伴随着免疫炎症紊乱。基于LASSO分析，本研究从上述2种分析方法得到的共同基因中鉴定出6个关键基因，即 CD163、RNASE2、HP、FAM20A、MCEMP1和FAM198B。ROC曲线分析结果提示，这些关键基因对心肌梗死的发生均具有较高的诊断价值。本研究还发现，在梗死恢复过程中，这些基因的表达呈现明显下降趋势。GSEA分析结果表明，这些关键基因涉及的生物学改变主要与糖脂代谢、活性氧、免疫炎症等有关。
CD163分子是Ⅰ型膜蛋白，蛋白表达限于单核细胞/巨噬细胞系, CD163抗原特异性释放机制可能在炎症调节过程中起重要作用。在动脉粥样硬化过程中, CD163+巨噬细胞能促进血管生成和增加血管通透性，并伴随炎症反应；此外，破裂的冠状动脉斑块中 CD163 的表达增加了心肌梗死和冠心病发生风险[24]。 RNASE2 编码的蛋白质是非分泌型核糖核酸酶，属于胰核糖核酸酶家族[25]。 FAM20A 是激酶编码基因家族的成员，本身不具有激酶活性，通过与家族成员 FAM20C 形成复合物，增强 FAM20C 的激酶活性，从而使分泌通路内的蛋白磷酸化[26]。 MCEMP1 又称 C19ORF59 , 可编码表达于肥大细胞、巨噬细胞等的跨膜蛋白[27]。研究[27]表明, MCEMP1 基因可能是卒中诊断和预后评估的新生物标志物。生物信息学研究[28-29]提示, RNASE2、FAM20A和 MCEMP1 可能是心肌梗死的关键基因，然而这些基因在心肌梗死及梗死恢复过程中的具体作用机制有待进一步研究。HP基因编码一种前蛋白，经处理后产生α链和β链，结合为四聚体产生触珠蛋白，其基因型与急性心肌梗死的发病风险密切相关，并且能决定心肌梗死面积[30-31]。 FAM198B(GASK1B) 是目前未知功能的新基因，可能编码定位在高尔基体上的膜结合糖蛋白，被认为参与癌症的转移与进展[32]。本研究分析结果显示, FAM198B 可能是心肌梗死和反映心肌梗死恢复的关键基因。
本研究尚存在一定局限性: 首先，本研究基于GEO数据库中下载的心肌梗死患者PBMCs的测序数据集进行分析，还需要进一步开展分子生物学实验来验证这些基因在心肌梗死中的表达及其可能的生物学机制; 其次，本研究仅分析了现有的样本数据，未来还需基于更大的样本规模和更详细的样本类型进一步深入研究; 最后，由于缺乏临床资料，本研究未分析这些关键基因的表达与心肌梗死患者年龄、性别、射血分数等临床指标的相关性。
综上所述，本研究基于PBMCs确定了6个可能与心肌梗死密切相关的关键基因，即 CD163、RNASE2、HP、FAM20A、MCEMP1和FAM198B, 为心肌梗死的诊断、评估、管理和预后分析提供了新的思路。
-
![]()
图 2 Western blot检测细胞中Cleaved-caspase3和Total-caspase3蛋白表达情况
A: Western blot检测结果; B: Cleaved-caspase3和Total-caspase3表达情况
(与对照组比较, *P < 0.05; 与病毒组比较, #P < 0.05)。![]()
图 3 Western blot检测各组细胞中TNF-α和IL-6蛋白表达情况
A: Western blot检测结果; B: TNF-α和IL-6表达情况(与对照组比较, *P < 0.05; 与病毒组比较, #P < 0.05)。
表 1 HCMV感染后不同时点各组细胞增殖活性比较(x±s)
组别 n OD450 nm 24 h 48 h 72 h 对照组 9 0.31±0.04 0.53±0.06 0.84±0.09 病毒组 9 0.36±0.04 0.42±0.05* 0.56±0.06* 更昔洛韦组 9 0.35±0.03 0.55±0.03# 0.86±0.05# 甘草酸低剂量组 9 0.33±0.02 0.47±0.04 0.62±0.03 甘草酸中剂量组 9 0.34±0.04 0.51±0.05# 0.79±0.08# 甘草酸高剂量组 9 0.32±0.03 0.59±0.06# 0.91±0.07# 与对照组比较, * P < 0.05; 与病毒组比较, #P < 0.05。 
下载: 导出CSV
表 2 HCMV感染后不同时点各组细胞总凋亡率比较(x±s)
% 组别 n 24 h 48 h 72 h 对照组 9 7.86±0.94 8.09±0.85 8.13±0.76 病毒组 9 7.51±0.83 22.38±2.14* 25.97±2.49* 更昔洛韦组 9 8.04±0.91 12.74±0.98# 11.68±1.05# 甘草酸组 9 7.74±0.79 13.39±1.16# 12.71±1.11# 与对照组比较, * P < 0.05; 与病毒组比较, #P < 0.05。
下载: 导出CSV
表 3 各组细胞Cleaved-caspase3和Total-caspase3蛋白相对表达量比较(x±s)
组别 n Cleaved-caspase3 Total-caspase3 对照组 9 0.99±0.10 1.00±0.12 病毒组 9 2.09±0.22* 0.98±0.07 更昔洛韦组 9 1.23±0.11# 1.01±0.09 甘草酸组 9 1.38±0.14# 0.99±0.08 Cleaved-caspase3: 活化含半胱氨酸的天冬氨酸蛋白水解酶3; Total-caspase3: 含半胱氨酸的天冬氨酸蛋白水解酶3总蛋白。与对照组比较, * P < 0.05; 与病毒组比较, #P < 0.05。
下载: 导出CSV
表 4 HCMV感染后不同时点各组细胞培养液中TNF-α、IL-6水平比较(x±s)
pg/mL 时点 组别 n TNF-α IL-6 24 h 对照组 9 14.37±2.09 76.94±16.82 病毒组 9 22.63±2.74* 154.69±20.67* 更昔洛韦组 9 16.98±1.91# 106.77±15.39# 甘草酸组 9 15.91±1.63# 117.84±12.71# 48 h 对照组 9 15.11±2.26 80.09±19.74 病毒组 9 63.94±7.15* 328.96±47.95* 更昔洛韦组 9 32.74±4.31# 121.39±13.28# 甘草酸组 9 35.68±3.86# 132.74±13.95# 72 h 对照组 9 17.81±3.11 78.82±20.03 病毒组 9 74.08±8.24* 371.52±51.33* 更昔洛韦组 9 35.86±3.92# 136.84±14.82# 甘草酸组 9 40.11±4.24# 141.79±15.08# TNF-α: 肿瘤坏死因子-α; IL-6: 白细胞介素-6。
与对照组比较, * P < 0.05; 与病毒组比较, #P < 0.05。
下载: 导出CSV
表 5 各组细胞TNF-α、IL-6蛋白相对表达量比较(x±s)
组别 n TNF-α蛋白相对表达量 IL-6蛋白相对表达量 对照组 9 1.00±0.09 1.00±0.11 病毒组 9 2.24±0.29* 2.46±0.31* 更昔洛韦组 9 1.29±0.14# 1.25±0.11# 甘草酸组 9 1.42±0.16# 1.51±0.17# 与对照组比较, * P < 0.05; 与病毒组比较, #P < 0.05。
下载: 导出CSV
-
[1] RODRÍGUEZ-SÁNCHEZ I, MUNGER J. Meal for two: human Cytomegalovirus-induced activation of cellular metabolism[J]. Viruses, 2019, 11(3): 273. doi: 10.3390/v11030273
[2] GRIFFITHS P, REEVES M. Pathogenesis of human Cytomegalovirus in the immunocompromised host[J]. Nat Rev Microbiol, 2021, 19(12): 759-773. doi: 10.1038/s41579-021-00582-z
[3] 窦冬冬, 张楠, 潘峰. 先天性人巨细胞病毒感染对肝脏的损伤作用及相关血清标志物的分析[J]. 肝脏, 2019, 24(10): 1161-1163. doi: 10.3969/j.issn.1008-1704.2019.10.022 [4] SELYUTINA O Y, POLYAKOV N E. Glycyrrhizic acid as a multifunctional drug carrier-From physicochemical properties to biomedical applications: a modern insight on the ancient drug[J]. Int J Pharm, 2019, 559: 271-279. doi: 10.1016/j.ijpharm.2019.01.047
[5] 王永洪, 兰婉莹, 陈恬. 甘草酸体外抗人巨细胞病毒AD169的作用[J]. 中成药, 2016, 38(10): 2121-2127. https://www.cnki.com.cn/Article/CJFDTOTAL-ZCYA201610006.htm [6] SHAO Y Z, QI Y, HUANG Y J, et al. Human Cytomegalovirus miR-US4-5p promotes apoptosis via downregulation of p21-activated kinase 2 in cultured cells[J]. Mol Med Rep, 2017, 16(4): 4171-4178. doi: 10.3892/mmr.2017.7108
[7] 马南兰, 茅惠娟, 张亚平, 等. 体外膜肺氧合技术治疗血行播散性肺结核合并巨细胞病毒性肺炎并发急性呼吸窘迫综合征1例报道[J]. 实用临床医药杂志, 2021, 25(4): 29-33. doi: 10.7619/jcmp.20201705 [8] 贾阳杰, 吕亚丽, 刘丽宏. 姜黄素对人巨细胞病毒感染人胚肺成纤维细胞后细胞凋亡的影响[J]. 中国医药导报, 2018, 15(34): 4-9, 22. https://www.cnki.com.cn/Article/CJFDTOTAL-YYCY201834002.htm [9] 刘杨, 黄媛, 廖毅, 等. 黄芩苷体外抗人巨细胞病毒活性及其对感染人胚肺成纤维细胞凋亡的影响[J]. 医药导报, 2017, 36(10): 1124-1129. https://www.cnki.com.cn/Article/CJFDTOTAL-YYDB201710010.htm [10] SUN Z G, ZHAO T T, LU N, et al. Research progress of glycyrrhizic acid on antiviral activity[J]. Mini Rev Med Chem, 2019, 19(10): 826-832. doi: 10.2174/1389557519666190119111125
[11] LI J J, XU D G, WANG L L, et al. Glycyrrhizic acid inhibits SARS-CoV-2 infection by blocking spike protein-mediated cell attachment[J]. Molecules, 2021, 26(20): 6090. doi: 10.3390/molecules26206090
[12] 廖卫波, 申宝德, 沈成英, 等. 绿原酸与甘草酸联合应用体外抗呼吸道合胞病毒作用研究[J]. 解放军药学学报, 2017, 33(1): 12-16. https://www.cnki.com.cn/Article/CJFDTOTAL-JFJN201701004.htm [13] ROGERS C, FERNANDES-ALNEMRI T, MAYES L, et al. Cleavage of DFNA5 by caspase-3 during apoptosis mediates progression to secondary necrotic/pyroptotic cell death[J]. Nat Commun, 2017, 8: 14128.
[14] ZEDTWITZ-LIEBENSTEIN K, DIAB-ELSCHAHAW M, FRASS M. Human Cytomegalovirus infection in nonimmunocompromised patients: a retrospective analysis and review of the literature[J]. Intervirology, 2016, 59(3): 159-162.
[15] LAU B, KERR K, GU Q, et al. Human Cytomegalovirus long non-coding RNA1. 2 suppresses extracellular release of the pro-inflammatory cytokine IL-6 by blocking NF-κB activation[J]. Front Cell Infect Microbiol, 2020, 10: 361.
-
期刊类型引用(2)
1. 于琪力，葛雪，刘丽. 冠状动脉粥样硬化性心脏病患者心力衰竭危险因素分析. 新乡医学院学报. 2024(06): 560-565 . 
百度学术
2. 孙宁，李飞. 心肌梗死患者血清miR-34a-5p和C1q/肿瘤坏死因子相关蛋白-9水平及其与心肌重构的相关性. 检验医学与临床. 2024(19): 2785-2789 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 194
- HTML全文浏览量: 46
- PDF下载量: 12
- 被引次数: 2
苏公网安备 32100302010246号