基于生物信息学的乳腺癌细胞焦亡相关基因多组学分析及相关中药筛选预测

细胞焦亡（pyroptosis）是一种由炎症小体引发的裂解性程序性细胞死亡形式[1]，该过程由促炎性半胱氨酸蛋白酶（Caspase）-1、4、5、11的底物gasdermin D介导并伴随炎症和免疫反应[2]。研究发现细胞焦亡可以影响肿瘤的增殖、侵袭和转移[3]。因此，有报道认为诱导肿瘤细胞焦亡是一种潜在的癌症治疗策略[4]。

乳腺癌（breast cancer）是目前最常见的恶性肿瘤之一，是女性癌症死亡的最常见原因[5]，其终生患病概率为12.3% [6] ，对女性群体有极大的危害性。受益于诊断及影像技术发展与早期筛查，乳腺癌的死亡率已经大幅降低[7]，但由于存在一定复发及转移风险，乳腺癌患者的预后仍然不佳[8-9]。目前乳腺癌的治疗主要包括手术、内分泌治疗、放疗、化疗、靶向治疗等[10-12]，但这些治疗往往存在较严重的不良反应[13-15]。因此，鉴于细胞焦亡对癌症存在调控作用，以此研究方向为指导或可开发更加安全、有效的新药物。

中药对于乳腺癌治疗的有效性被广泛报道[16]，而挖掘与细胞焦亡相关的中药则更有利于指导相关研究与临床应用。本研究通过生物信息学方法从多组学角度全面分析细胞焦亡相关基因的表达及其与乳腺癌生存、预后、免疫、肿瘤微环境等方面的相关性，并在此基础上筛选出可能对细胞焦亡有调控作用的中药。

1资料与方法

1.1 数据收集及处理

在TCGA数据库（https://portal.gdc.cancer.gov/）中分别获取乳腺癌相关转录组、临床、突变及拷贝数数据。对其数据进行解压、ID转换并将其格式由FPKM（fragments per kilobase per million）转换为TPM（trans per million）。在GEO数据库（https://www.ncbi.nlm.nih.gov/geo/）中检索“breast cancer survival”并限定数据类型（expression profiling by array）及物种（homo sapiens）进行筛选，获取基因表达及临床数据，然后对其数据进行基因符号注释。

1.2 乳腺癌细胞焦亡基因突变、拷贝数变异、差异、生存及预后分析

通过检索细胞焦亡相关文献获得细胞焦亡相关基因，将样品突变数据中的细胞焦亡相关基因利用R语言“maftools”包及Perl语言进行突变类型、频率、数目、碱基改变及负荷的统计分析，并将结果以瀑布图表示；然后检索UCSC Xena（https:// xena.ucsc.edu/）数据库并下载乳腺癌拷贝数相关数据，使用Perl语言提取样品中细胞焦亡相关基因的拷贝数，再使用R语言将样品细胞焦亡相关基因拷贝数结果绘制成拷贝数频率图及圈图；将TCGA及GEO数据库中的样品表达量数据使用R语言中的“limma”和“sva”包进行合并，并提取细胞焦亡基因的表达量，再将2个数据库所得临床数据中的生存时间、生存状态数据进行合并，使用R语言的“survival”和“surminer”包进行生存和预后分析并可视化。

1.3 乳腺癌细胞焦亡分型及分型基因生存分析、基因集变异分析（gene set variation analysis，GSVA）、免疫细胞差异分析（single sample gene set enrichment analysis，ssGSEA）和主成分分析（principal component analysis，PCA）和差异性分析

根据细胞焦亡的表达量使用R语言的“ConsensuClusterPlus”对乳腺癌细胞焦亡基因进行聚类分型，其中聚类算法为K-means，距离类型为欧氏距离；然后使用R语言的“RcolorBrewer”“GSVA”“PRGcluster”等将分型后的两组细胞焦亡基因进行生存分析、细胞焦亡基因表达分析、GSVA、ssGSEA、PCA并可视化，最后将细胞焦亡分型的两组样品基因表达量进行差异分析，并获得分型之间的差异基因。

1.4 乳腺癌细胞焦亡分型差异基因基因本体（gene ontology，GO）富集分析和京都基因与基因组百科全书（Kyoto encyclopedia of genes and genomes，KEGG）通路富集分析

将“1.3”项所得的细胞焦亡分型之间的差异基因使用R语言的“clusterProfiler”和 “enrichplot”包分别进行GO的生物过程（biological process，BP）、分子功能（molecular function，MF）及细胞成分（cellular component，CC）富集分析和KEGG通路富集分析，按照P＜0.005进行筛选，并绘制GO及KEGG分析的气泡图。

1.5 乳腺癌细胞焦亡分型预后相关差异基因再分型及其生存与表达分析

根据“1.3”项所得分型差异基因的生存分析结果探寻预后相关的差异基因，再通过R语言聚类分析将预后相关差异基因进行聚类分型，选取准确性最高的聚类结果作为分型结果。根据分型结果比较不同分型之间的生存情况，并绘制生存曲线与差异基因热图。然后将“1.3”项所得基因分型的结果与“1.2”项所得细胞焦亡基因表达量数据使用R语言进行差异分析并绘制箱线图。

1.6 预后模型的构建、分析及其可视化

①将包含单因素显著基因表达数据和临床生存数据的集合进行随机分组，分为训练（train）组及验证（test）组；然后构建基因表达量及生存时间、状态的套索算法（least absolute shrinkage and selection operator，Lasso）回归模型并进行交叉验证，选择其中误差最小的点对应的基因作为Lasso回归模型显著基因，结合显著基因与其表达量构建COX回归模型。②利用训练组得出模型计算公式，即基因表达量乘以系数再求和，进一步根据公式计算训练组风险得分，根据得分中位数区分高、低风险组；验证组同理，即通过公式计算风险得分，并根据训练组的中位数区分验证组的高、低风险组；进行生存差异分析并得出差异的P值、受试者工作特征曲线（receiver operating characteristic curve，ROC）、风险曲线及列线图。③将细胞焦亡分型、基因分型结果通过R语言分别进行预后分析，比较不同分型之间的差异性及风险得分并可视化；再分析高、低风险组之间细胞焦亡相关基因的差异性，并绘制箱线图。

1.7 免疫细胞浸润与免疫细胞含量分析

将“1.2”项中TCGA和GEO数据库所得数据的并集利用R语言CIBERSORT命令获得总量为1的免疫细胞相对含量；将免疫细胞相对含量的数据结果根据P＜0.05进行过滤，并与“1.6”项风险得分的结果取交集，最后对所有免疫细胞进行循环计算，得到风险得分、基因与免疫细胞的相关性结果。

1.8 肿瘤微环境分析、肿瘤突变负荷分析及干细胞分析

①将“1.2”项中TCGA和GEO数据库所得数据的并集使用R语言“estimate”包对肿瘤微环境进行基质细胞、免疫细胞、肿瘤纯度分析以及综合评分；利用评分结果与“1.6”项所得高、低风险组进行肿瘤微环境差异分析，并采用小提琴图表示。②将所有样品的风险数据与染色体突变位置数据利用Perl语言分析得出其按照高、低风险组区分的肿瘤突变数据，并使用R语言绘制瀑布图；利用R语言的“reshape2”包结合“1.7”项中样品风险数据及“1.5”项中分型数据获取肿瘤突变负荷数据，并以箱线图表示。③将泛癌干细胞打分数据结合“1.7”项中样品风险数据，利用R语言分析获得相关性结果并绘制散点图。

1.9 统计分析

使用R v.4.1.2软件进行统计分析。2组独立样品数据比较采用t检验。3组及以上的数据，则采用单因素方差分析（ANOVA）及Kruskal-Wallis秩和检验，以P＜0.05表示存在统计学意义。

1.10“细胞焦亡相关靶点-有效成分-中草药”网络构建与药物筛选

将细胞焦亡相关靶基因使用TCMSP数据库（https://old.tcmsp-e.com/tcmsp.php）查找靶点相关的有效成分并根据口服生物利用度（oral bioavailability，OB）≥30%、类药性（drug likeness，DL）≥0.18进行筛选，然后通过有效成分筛选相关中药，使用CytoscapeV 3.8.0构建“靶点-成分-中药”网络并通过拓扑学分析获取核心中药。

2 结果

2.1数据收集与处理结果

2.2 乳腺癌细胞焦亡基因突变分析、拷贝数变异分析、差异分析、生存分析及预后分析结果

经过检索文献获得细胞焦亡基因52个，具体如表1所示。通过分析乳腺癌中细胞焦亡相关基因的拷贝数变异和基因突变的发生率，可知在986个突变数据样品中有394个样品发生了变异，突变比例占39.96%，其中TP53突变率最高，达34%，其他为CASP8（2%）、NLRC4（1%）、NLRP3（1%）、NLRP1（1%）、NLRP2（1%）、PLCG1（1%）、NLRP7（1%）、SCAF11（1%）、GSDMC（1%）、NOD1（1%），如图1-A所示；图1-B则显示大多数细胞焦亡相关基因存在拷贝数的变化，而拷贝增加较缺失频率略大；由图1-C可进一步看出拷贝数变化的细胞焦亡相关基因所在染色体的位置；由图1-D可以看出71%的细胞焦亡基因在正常样品和肿瘤样品间表达存在显著性差异，其中BAK1、BAX、CSCP3、CHMP2A、CHMP2B、CHMP4B、CHMP4C、CHMP6、CYCS、CSDMD、IL18、IRF1、IRF2、AIM2、CASP6、CASP8、CSDMC、NLRP6、NLRP7、NOD2、PYCARD、GZMA为肿瘤组表达上调，CASP1、CASP4、CHMP3、ELANE、GSDMD、IL1B、TP63、GPX4、GSDMB、IL6、NLRP1、NLRP3、NOD1、PJVK、SCAF11为肿瘤组表达下调。将TCGA与GEO数据合并后进行生存分析，结果见表2，其中风险率（hazard ratio，HR）＞1的基因为高风险基因，P＜0.05则说明该基因为预后相关，根据图2-A可知，生存曲线中的预后相关基因有CASP1、CASP4、CASP8、CASP9、CHMP2A、CHMP2B、CHMP4B、CHMP4C、CHMP6、CHMP7、CYCS、GPX4、GSDMC、GSDMD、GZMA、GZMB、IL1A、IL18、IRF1、IRF2、NLRC4、NLRP1、NLRP3、NLRP6、NLRP7、PYCRD1、SCAF11、TIRAP、TP53、TP63，其大多为高表达则预后较好；然后进一步构建预后网络，如图2-B所示，其中紫色节点表示高风险基因，绿色节点表示低风险基因、节点大小表示P值，红色连线表示正向共表达，蓝色连线表示负向共表达，由连线密度可知其中存在密切的联系。

2.3 乳腺癌细胞焦亡分型及分型基因生存分析、GSVA、ssGSEA、PCA和差异分析结果

通过对细胞焦亡相关基因表达量数据进行聚类分析，可将样品分为A和B 2个亚型（图3-A）；对2个亚型进行生存分析，如图3-B所示，可知其存在统计学差异（P＜0.05），且A亚型预后较好；对2个亚型基因表达进行分析，如图3-C所示，可得A亚型多数呈现高表达，B亚型则相反；通过GSVA分析，则能看出A亚型主要对免疫相关通路呈现高表达，B亚型主要表现为低表达（图3-D）；从ssGSEA分析结果（图3-E）可知，22个在两亚型间存在统计学差异的免疫细胞，其中活化B细胞、活化的CD4细胞、活化的CD8细胞、活化的树突状细胞、CD56bright自然杀伤细胞、嗜酸性粒细胞、γδT细胞、未成熟B细胞、未成熟的树突状细胞、骨髓源性抑制细胞、巨噬细胞、肥大细胞、单核细胞、自然杀伤T细胞、浆细胞样树突状细胞、调节性T细胞、滤泡辅助T细胞、1型T辅助细胞、17型T辅助细胞、2型T辅助细胞均在A亚型中含量高；唯有CD56dim自然杀伤细胞在B亚型中含量高；PCA分析（图3-F）显示，根据细胞焦亡相关基因的表达量可将A与B亚型明显区分；最后通过设置|log 2FC| ＞0.585、P＜0.05，筛选出细胞焦亡分型的显著差异基因有1223个。

2.4 乳腺癌细胞焦亡分型差异基因GO富集分析和KEGG通路富集分析

通过对1223个细胞焦亡分型的显著差异基因进行GO分析，结果表明，BP主要涉及T细胞活化、细胞间黏附的调节、白细胞-细胞黏附、细胞因子产生的正调控、单核细胞分化、细胞活化的正调控、T细胞黏附的调节等免疫过程；MF主要涉及受体配体活性、信号受体激活剂活性、细胞因子受体结合、细胞因子活性、免疫受体等；CC主要涉及质膜筏的外侧、内吞囊泡、膜筏、膜微区等（图4-A）。KEGG通路富集分析显示，主要包括细胞因子-细胞因子受体相互作用、趋化因子信号通路、细胞黏附分子、爱泼斯坦-巴尔病毒感染、人类T细胞白血病病毒1感染、病毒蛋白与细胞因子和细胞因子受体的相互作用、造血细胞谱系、结核、吞噬体、破骨细胞分化、类风湿关节炎、核因子-κB（nuclear factor-κB，NF-κB）信号通路、Th17细胞分化、TNF信号通路等炎症、免疫相关通路（图4-B）。

2.5 乳腺癌细胞焦亡分型预后相关差异基因再分型及其生存分析与表达分析结果

通过对1223个细胞焦亡分型预后显著差异基因进行再分型，可得到3组分型基因，分别为I、II、III，如图5-A所示；对这3组分型基因进行生存分析，3组分型差异存在显著性差异（P＜0.05），且其中I组预后最好，II组预后最差（图5-B）；对3组分型基因进行表达分析，可见I组主要表现为高表达，II组主要呈现低表达（图5-C），结合生存分析及基因表达分析可知，相关基因呈现高表达则预后较好，反之则预后较差，与预期一致；接下来对3组细胞焦亡基因分型进行差异分析，可知3组之间超过80%细胞焦亡基因存在显著性差异，其中超过60%细胞焦亡基因为I组高表达，超过70%细胞焦亡基因为II组低表达（图5-D）。

2.6预后模型构建、预后分析及其可视化结果

2.6.1预后模型Lasso回归结果及交叉验证结果如图6-A、B所示，COX回归模型得出2组及总体的风险值；预后模型构建过程如图6-C所示，即首先根据细胞焦亡相关基因的表达量将样品分为2个亚型（A、B亚型），然后通过差异分析找出该2亚型中存在的差异基因，通过分析差异基因的表达量对基因进行分型（分型I、分型II和分型III）；对差异基因进行生存分析以找出预后相关的差异基因并构建预后相关模型；再根据预后模型将患者分为高、低风险组，最终根据预后模型预测患者的生存期。

2.6.2 细胞焦亡分型风险得分的差异分析结果表明A、B两亚型间差异有统计学意义，其中分型B的风险评分更高，如图6-D所示；基因分型风险得分的差异分析结果表明分型I、II及分型II、III之间差异有统计学意义，分型I、III之间则无统计学意义，同时可以看出分型II的风险评分最高，如图6-E所示。

2.6.3 高、低风险组细胞焦亡基因表达量的差异分析结果如图6-F所示，高、低风险组中细胞焦亡基因绝大多数（41/52）存在显著差异，且这些基因中除GSDMC、CYCS、BAK13个基因在高风险组为高表达外，其余均表现为低表达；此外由生存曲线（图6-G）可知，不论训练组、验证组还是总体组（训练组及验证组的总和）其高、低风险组的生存时间都存在显著差异，说明构建的预后模型可区分高、低风险组的患者；由ROC曲线（图6-H）可知，训练组线下面积最大而验证组线下面积最小，因此通过预后模型预测生存期训练组准确性最大而验证组最小，而不论训练组、验证组还是总体组线下面积都大于0.65，因此预后模型预测生存期准确性均较高；列线图可提供预后模型对患者不同临床特征生存期的预测，由图6-I可知女性、年龄62岁、高风险评分的患者预后差，其1年、3年、5年生存率分别为96.7%、83.4%、70.9%，图6-J则显示了列线图生存率预测的准确性；最后根据患者高、低风险排序绘制风险曲线（图6-K），不论训练组、验证组还是总体组3者都是随风险值增高而死亡人数增多，与预测相符，而风险热图则提示所列基因皆为低风险基因。

2.7 免疫细胞浸润与免疫细胞含量分析结果

通过CIBERSORT计算求得每种免疫细胞在每个样品中的相对含量作为免疫细胞浸润结果；利用该结果进一步得出风险得分与免疫细胞相关性，其中存在显著负相关性的有幼稚B细胞、静息树突状细胞、巨噬细胞M1、单核细胞、γδT细胞、CD8 T细胞、静息肥大细胞、静息CD4记忆T细胞、浆细胞；存在显著正相关性的有活化树突状细胞、M0巨噬细胞、M2巨噬细胞、活化肥大细胞、静息自然杀伤细胞，见图7-A。参与模型构建的基因与免疫细胞相关性如图7-B所示，其中大部分基因与免疫细胞存在相关性，且呈负相关性较正相关性者多。

2.8 肿瘤微环境分析、肿瘤突变负荷分析及干细胞分析结果

2.8.1肿瘤微环境分析首先获得肿瘤微环境TCGA、GEO两数据库临床样品的评分结果，包含基质细胞评分、免疫细胞评分、总评分以及肿瘤纯度4个方面；然后将肿瘤微环境评分结果按照高、低风险分组进行差异分析，由图8-A可知基质细胞评分、免疫细胞评分及总评分在高、低风险组间均具有统计学差异，且肿瘤微环境评分均为高风险组更低；使用瀑布图展示高、低风险组的突变数据，如图8-B、C所示，2组的突变率较接近（84.12%、85.89%），而高风险组突变率略高于低风险组。

2.8.2肿瘤突变负荷相关性分析结果显示高、低风险组间肿瘤突变负荷存在统计学差异，其中高风险组患者的肿瘤突变负荷较高（图8-D）；对3组基因分型样品的风险评分及肿瘤突变负荷进行相关性分析，由图8-E可知二者相关性有统计学意义，且呈正相关。

2.8.3干细胞相关性分析结果（图8-F）显示干细胞指数与风险评分之间存在相关性（P＜0.05），且二者呈正相关。

2.9 “细胞焦亡靶点-成分-中药”网络构建与中药筛选结果

通过检索TCMSP数据库，找到18个细胞焦亡靶点数据，将与该18个细胞焦亡靶点相关联的34种有效成分（表3）及380味中药进行网络构建（图9-A）并进行分析，分别通过节点度（degree）值筛选核心成分和核心中药，核心成分（度值＞2.5倍中位数）主要包括槲皮素（quercetin）、β-谷甾醇（β-sitosterol）、木犀草素（luteolin）、山柰酚（kaempferol）、β-胡萝卜素（β-carotene）、金合欢素（acacetin）、黄芩素（baicalein）、汉黄芩素（wogonin）、芦荟大黄素（aloe-emodin）、表没食子儿茶精没食子酸酯[(−)-epigallocatechin-3-gallate]、川陈皮素（nobiletin）、人参皂苷Rh 2 （ginsenoside Rh 2 ）、柚皮素（naringenin）；核心中药（度值＞2.5倍中位数）主要包括木蝴蝶、红花、余甘子、白果、半枝莲、牛膝、菊花、香薷、马齿苋、连翘、连钱草、败酱草、枇杷叶、关黄柏、皂角刺、马鞭草、金银花17味。对380味中药进行性味归经统计分析，其中四气以温、寒、微寒、平为主；五味以苦、辛、甘为主；归经主要归肝、肺经，见图9-B、C、D。

3 讨论

细胞焦亡作为一种受调节的促炎形式的细胞死亡，其在形态学、机制学和病理生理学上不同于包括细胞凋亡和坏死在内的其他形式的细胞死亡，其特点是质膜快速破裂，随后释放细胞内容物和促炎介质，包括IL-1β、IL-18和警报素高迁移率族蛋白1（high mobility group protein B1，HMGB-1）等[17]。细胞焦亡已被确认与心血管疾病、传染病、卵巢癌、肺癌等疾病有关，但是其与乳腺癌的关系却仍缺少相关研究证实[18-20]，且相关的生物信息学研究也缺乏涉及中药方面的分析，因此本研究对3者的关系进行了一定程度探究。

本研究首先根据乳腺癌细胞焦亡基因突变分析可知TP53突变频率远高于其他基因，而研究证实TP53的突变可能是除乳腺癌1型易感蛋白（breast cancer type 1 susceptibility protein，BRCA1）、乳腺癌2型易感蛋白（Breast cancer type 2 susceptibility protein，BRCA2）外的另一个导致乳腺癌的关键因素，其与三阴性乳腺癌有关[21]；肿瘤组和正常组的对比则可发现细胞焦亡相关基因的表达大多数在二者中存在显著差异，而肿瘤组表达上调者略多于表达下调者，研究认为[4]细胞焦亡与癌症的关系复杂，不同组织和遗传背景的细胞焦亡对癌症的影响不同而呈现促癌和抑癌2种作用。乳腺癌细胞焦亡基因生存分析显示，生存曲线差异显著且高风险的基因中GSDMC高表达与乳腺癌患者较差的存活率相关[22]，并可介导癌细胞中的非典型细胞焦亡途径，导致肿瘤坏死；CASP3在肿瘤组中低表达，有动物实验证明其表达上调可促进乳腺癌细胞凋亡[23]；PLCG1作为磷脂酶C的主要亚型，是细胞信号传导的重要介质，其在各种癌症中经常富集和突变，并参与肿瘤发生的过程，包括增殖、迁移和侵袭[24]；CYCS有文献报道称其可作为乳腺癌细胞凋亡的生物标志物在化疗后表达增加[25]；分析研究认为BAK1与乳腺癌的骨转移有关[26]；PICAKA所参与的通路则有促成雌激素诱导的乳腺癌细胞增殖和内分泌抵抗的作用[27]；IL-6在包括乳腺癌在内的多种肿瘤组织中过度表达，且在乳腺癌生长、转移及治疗抵抗方面起作用[28]。总体来说细胞焦亡基因在乳腺癌中突变频率不高，但其突变大部分与预后相关。

通过将乳腺癌细胞焦亡基因进行聚类及分型，发现A亚型患者生存结局优于B亚型，且A亚型中细胞焦亡基因多数呈现高表达；GSVA、ssGSEA分析显示，炎症相关通路、免疫细胞均在A亚型中富集，A亚型特征为细胞焦亡基因高表达、免疫激活和免疫细胞浸润、有生存优势，B亚型与之相反。对细胞焦亡基因分型进行预后差异分析，得出差异基因1223个，并将之再次聚类获得3组基因分型，同样得到I组细胞焦亡基因表达最高、生存预期最好，II组细胞焦亡基因表达最低、生存预期最差，2次分型结果一致。通过构建预后模型并进行风险评分，得出高风险组预后差，低风险组细胞焦亡基因高表达，且有较高的准确性。不论免疫浸润还是肿瘤微环境分析都得出同一个结论，即高风险组预后较差、细胞焦亡基因低表达、生存期较短、免疫细胞浸润少；低风险组则与之相反，这也反映了细胞焦亡与乳腺癌及其肿瘤微环境、免疫反应之间的关系。研究认为，肿瘤微环境与乳腺癌的产生与进展密切相关[29-30]，而其则表现为大量炎性细胞浸润，如B细胞浸润与晚期乳腺癌相关、CD4 +、CD8 +T 细胞的浸润与乳腺癌预后相关[30-31]。

从上述分析可知，细胞焦亡与乳腺癌有显著的相关性，因此挖掘可调控细胞焦亡靶基因的中药则可为乳腺癌治疗药物的研发提供参考。乳腺癌属中医“乳岩”范畴，其发病主要为情志失调、脾胃损伤、冲任失调，痰浊、瘀血乘结乳房而生[32]。从数据挖掘结果来看，所得核心中药木蝴蝶、红花、余甘子、白果、半枝莲、牛膝、菊花、香薷、马齿苋、连翘、连钱草、败酱草、枇杷叶、关黄柏、皂角刺、马鞭草、金银花[33-49]等均具有抗肿瘤、调节免疫及抗炎作用，其所含成分中山柰酚可通过抑制Caspase-1的表达和活性、增加Toll样受体4和NLRP3的蛋白表达以及抑制GSDMD的分解来减少细胞焦亡[50]；动物实验证明槲皮素可阻断氧化应激诱导的软骨细胞焦亡来抑制骨关节炎[51]；木犀草素则被报道有助于减轻脓毒症诱导的肺损伤小鼠模型肺组织中Caspase-11依赖的细胞焦亡[52]；黄芩素可上调微小RNA-192-5p（miR-192-5p）/硫氧还蛋白互作蛋白（thioredoxin-interacting protein，TXNIP）轴而抑制NLRP3/Caspase-1通路，从而调控细胞焦亡和炎症[53]；汉黄芩素可通过抑制GSDMD介导的细胞焦亡，减轻顺铂诱导的小鼠心脏毒性[54]；表没食子儿茶精没食子酸酯可抑制Caspase-1激活和IL-1β分泌[55]；柚皮素可通过调节微小RNA-200b（miR-200b）/锌指蛋白802（zinc finger protein 802，JAZF1）轴抑制乳腺癌的发生[56]。根据生物信息学分析表明，乳腺癌中细胞焦亡过程涉及免疫浸润与炎性因子聚集，而炎性因子聚集多造成局部炎症发热，该过程与中医的“火毒”存在一定的相似性[57]，因此其治疗用药以苦、寒为主；而内火多因“郁”而成，毒则为“聚而不散”所致，故用药应当辛散温通；中医历来强调顾护正气，正气存内，邪不可干，故甘平补益之品当为辅。从归经来看细胞焦亡相关药物主要归肝、肺、脾、胃经，因龙虎回环、脾升胃降皆为周身气机之枢，故细胞焦亡、乳腺癌此类“郁”所致的病理过程或疾患应以此为归经。

4 结论

本研究采用生物信息学方法对TCGA、GEO数据库的乳腺癌数据进行细胞焦亡相关多组学分析，从而得出细胞焦亡与乳腺癌有显著的相关性，其涉及肿瘤微环境免疫细胞浸润相关机制。而对细胞焦亡起调控作用的中药主要为木蝴蝶、红花、余甘子、白果、半枝莲、牛膝、菊花、香薷、马齿苋、连翘、连钱草、败酱草、枇杷叶、关黄柏、皂角刺、马鞭草、金银花，其主要为苦寒，次为辛温，辅以甘平之品，而主要调节肝、肺、脾、胃等脏腑，发挥清热解毒、调和肝脾之功。本研究结果可为中药的临床应用与进一步研究提供指导与参考。

利益冲突所有作者均声明不存在利益冲突

参考文献（略）