先从 2 套公开转录组资料寻找候选
泊那替尼的心血管毒性涉及多个细胞过程。延世大学原州医学院（Yonsei University Wonju College of Medicine）团队在 International Journal of Molecular Sciences 发表一项公开数据再分析，先整合不同肿瘤细胞的药物反应，再检查相关基因是否也在人源心肌细胞中出现变化。
第一套资料 GSE186341 包含 11 条肿瘤细胞系对 32 种激酶抑制剂的转录反应。研究使用其中 464 个溶剂对照，以及泊那替尼和伊马替尼各 20 个样本。第二套资料 GSE217421 来自 6 个供者背景的人诱导多能干细胞（induced pluripotent stem cells，iPSC）来源心肌细胞，选入 22 个泊那替尼、21 个伊马替尼和 63 个对照样本。
研究先用肿瘤细胞数据提取在不同细胞背景中出现的药物相关候选，再用心肌细胞数据比较这些候选的转录变化。整项研究使用既有数据，没有新增细胞或动物实验。
对照样本的背景差异会影响结果
在肿瘤数据中，对照样本远多于泊那替尼样本，且对照的基因表达并不均一。团队用高变基因进行主成分分析，再将对照分成 5 组，分别估计药物处理与各组对照之间的表达差异。
选用前 20 个主成分时，累计保留约 91.4% 的表达变异。分组数则结合轮廓系数、组内离散程度、重采样稳定性及每组样本量决定。研究讨论了 6 组的部分稳定性指标较好，但最终采用 5 组，兼顾后续比较所需的样本数。
使用不同对照组，泊那替尼显著差异表达基因数从 0 到 6962 个不等；各组的变化幅度相关性也较低。仅取所有组共有的基因，得到 19 个共同下调基因。这个结果促使团队进一步整合各组的效应估计。
整合各组效应后，再检查心肌数据
团队用固定效应荟萃分析汇总各组的表达变化，并按方向一致程度、异质性及统计条件筛选。泊那替尼最终得到 2639 个整合差异表达基因，其中 1645 个没有在单个对照组的差异基因名单中出现。
在肿瘤数据中，这些下调基因主要富集于淋巴细胞与 T 细胞分化、细胞因子生成和细胞黏附等过程。排除伊马替尼也富集的过程后，研究列出 226 个在本次分析中偏向泊那替尼的过程，其中 8 个涉及心脏形态形成、心肌收缩、钙调节或血管结构。
心肌数据单独进行差异表达分析。与前述 2639 个基因取交集后，81 个基因在 2 套数据中均达到差异表达条件，其中 49 个变化方向相同。研究进一步关注 KCNN3、RAPGEF3、FOXO1、HK2、MCL1 和 GCH1 等与离子调节、代谢或线粒体状态有关的候选。
伊马替尼也按同一分析框架处理。在肿瘤数据中得到 1950 个整合差异基因，2 套数据的交集则为 4 个。这 4 个基因没有与泊那替尼的 81 个共同基因重合，显示这次数据分析中的药物反应存在差别。
窄屏可横向滑动查看完整图表。
基因变化提供机制研究线索，功能仍需实验检验
对 81 个共同基因开展通路富集后，没有结果通过多重检验校正。作者因此主要结合变化方向一致的具体基因及既有研究讨论其可能关联，而未据这些富集结果确定新的损伤机制。
肿瘤和心肌数据的细胞类型、药物浓度、暴露时长及测序条件不同。对照分层也可能保留细胞系本身的差别。论文将这些因素作为跨数据集解释时需要考虑的限制。
这项工作的用途是缩小后续实验的候选范围。转录变化还需要结合蛋白水平、离子通道、线粒体功能和心肌收缩检测逐项研究；本次没有给出患者心脏风险的预测性能。