《Journal of Hematology & Oncology》：Mitochondrial genetic variation as a novel predictor of graft-versus-host disease after stem-cell transplantation in myelodysplastic neoplasms
编辑推荐：
移植物抗宿主病（Graft-versus-host disease，GvHD）仍然是异基因造血细胞移植后发病和非复发死亡的主要原因。基于先前关于线粒体DNA（mitochondrial DNA，mtDNA）变异在骨髓增生异常肿瘤（myelodysplastic
移植物抗宿主病（Graft-versus-host disease，GvHD）仍然是异基因造血细胞移植后发病和非复发死亡的主要原因。基于先前关于线粒体DNA（mitochondrial DNA，mtDNA）变异在骨髓增生异常肿瘤（myelodysplastic neoplasms，MDS）中作用的证据，研究人员调查了受者和供者mtDNA变异对急性和慢性GvHD的预后影响。研究人员对来自国际血液和骨髓移植研究中心（Center for International Blood and Marrow Transplant Research，CIBMTR）的494对MDS供者-受者配对（发现队列）的移植前样本进行了全基因组测序（whole-genome sequencing，WGS）。使用病因特异性Cox比例风险模型评估mtDNA变异与II-IV级急性GvHD（aGvHD24）、III-IV级急性GvHD（aGvHD34）和慢性GvHD（cGvHD）之间的关联。使用曲线下面积（area under the curve，AUC）量化模型判别能力。研究结果在295例MDS的独立队列中通过靶向mtDNA深度测序进行评估。在发现队列中鉴定出的受者mtDNA基因关联中，与aGvHD34相关的MT-ND2（p = 1.1 × 10?3）和与cGvHD相关的MT-CYB（p = 3.6 × 10?3）被独立验证，而在发现队列中名义上与aGvHD34相关的MT-CO2在受者验证队列中同样显著（p = 3.2 × 10?3）。将受者mtDNA基因纳入CIBMTR aGvHD风险指数后，aGvHD34的判别能力提高了11%（p = 0.009）。总之，mtDNA变异可能为GvHD提供独立的预后信息，并在与标准临床参数整合时增强当前的风险分层能力。
一、研究背景与目的
移植物抗宿主病（Graft-versus-host disease，GvHD）是异基因造血细胞移植（allogeneic hematopoietic cell transplantation，allo-HCT）后发病和非复发死亡的主要原因，而allo-HCT是骨髓增生异常肿瘤（myelodysplastic neoplasms，MDS）唯一的治愈性疗法。除人类白细胞抗原（HLA）外，非HLA遗传因素也可能影响GvHD风险。既往研究已表明线粒体DNA（mitochondrial DNA，mtDNA）变异可预测allo-HCT结局，但受者和供者mtDNA变异对急性和慢性GvHD的具体影响尚未得到系统评估。本研究旨在评估受者和供者mtDNA变异及供者-受者mtDNA相互作用与II-IV级急性GvHD（aGvHD24）、III-IV级急性GvHD（aGvHD34）和慢性GvHD（cGvHD）的关联，并评估mtDNA特征能否在CIBMTR急性GvHD风险指数（Graft-versus-host Disease Risk Index，GRI）基础上进一步改善预测能力。该论文发表在《Journal of Hematology》。
二、研究方法概述
研究人员利用来自国际血液和骨髓移植研究中心（CIBMTR）的494对MDS供者-受者配对作为发现队列，对移植前样本进行全基因组测序（whole-genome sequencing，WGS），鉴定出2,666个受者变异和1,825个供者变异。常见变异（等位基因频率≥1%）进行单独检验；跨16个mtDNA基因/区域的基因水平关联使用SKAT-O方法在病因特异性Cox模型中进行评估，以死亡（未发生GvHD）作为竞争事件，并调整受者年龄、供者年龄、供者类型和GvHD预防方案。研究结果在295例MDS的独立验证队列中通过靶向mtDNA深度测序进行评估。模型判别能力使用曲线下面积（area under the curve，AUC）量化。
三、研究结果
受者mtDNA基因与GvHD的关联
受者常见变异均未通过Bonferroni校正（p < 1.43 × 10?4）。在基因水平上，经Bonferroni校正后，受者SKAT-O分析显示：与aGvHD24相关的基因包括D-Loop、MT-ND1、MT-ND2、MT-ND3和MT-tRNA；与aGvHD34相关的基因包括MT-ND1、MT-ND2、MT-ND4、MT-CYB和MT-ATP6；与cGvHD相关的基因包括D-Loop、MT-ND1、MT-CYB、MT-CO1、MT-CO2和MT-tRNA。
供者mtDNA基因与GvHD的关联
供者MT-ND4、MT-ND5、MT-CO1和MT-CO2与aGvHD24相关；但无一通过aGvHD34或cGvHD的校正。
供者-受者mtDNA相互作用
名义上显著的供者-受者相互作用包括：MT-rRNA（p = 0.038）和MT-tRNA（p = 0.013）与aGvHD24，以及MT-ND5（p = 0.039）和供者MT-ATP8-受者MT-ND2（p = 0.047）与cGvHD。受者MT-ND5变异（HR = 0.41，p = 0.02）和供者MT-ATP8-受者MT-ND2联合变异（HR = 0.46，p = 0.04）与较低的cGvHD风险相关。
线粒体单倍群与GvHD的关联
受者和供者分别观察到16个和18个单倍群。与单倍群H相比，受者单倍群I（HR = 2.67，p = 0.016）和N（HR = 4.11，p = 0.005）与cGvHD风险增加显著相关。aGvHD24或aGvHD34未观察到相关单倍群。
mtDNA基因对GRI预测能力的增量价值
GRI单独模型的AUC为aGvHD24 0.614和aGvHD34 0.597。加入受者mtDNA基因后判别能力适度改善（aGvHD24 AUC = 0.658，p = 6.05 × 10?3；aGvHD34 AUC = 0.689，p = 2.17 × 10?3），供者基因改善更大（aGvHD24 AUC = 0.696，p = 6.61 × 10?5；aGvHD34 AUC = 0.708，p = 5.93 × 10?4）。两者联合改善最大，aGvHD24 AUC达到0.721（p = 2.03 × 10?6），aGvHD34 AUC达到0.749（p = 3.32 × 10?6）。
独立验证队列的验证结果
在295例MDS受者的独立队列中，SKAT-O验证了受者MT-ND2与aGvHD34（p = 1.1 × 10?3）及MT-CYB与cGvHD（p = 3.6 × 10?3）的关联通过Bonferroni校正。MT-CO2在发现队列中名义上与aGvHD34相关（p = 0.004），在验证队列中也显著（p = 3.2 × 10?3）。aGvHD24的关联未得到重复验证。将受者mtDNA基因加入GRI后，aGvHD34的预测显著改善（AUC = 0.722，p = 0.009），但aGvHD24未显著改善（AUC = 0.669，p = 0.083）。
四、讨论与结论
线粒体功能障碍可能通过改变氧化磷酸化、活性氧生成、免疫细胞代谢和炎症信号传导影响GvHD。预处理相关的组织损伤也可能释放mtDNA，作为损伤相关分子模式（damage-associated molecular pattern，DAMP）放大先天免疫反应。MT-ND2和MT-CYB分别编码呼吸链复合物I和III的组分，这些基因的变异可能影响电子传递和氧化还原稳态，但需要功能学研究来确定因果关系。
研究结论：这些发现确定了与GvHD相关的供者和受者mtDNA特征，并独立验证了选定的受者关联。局限性包括缺乏匹配的非造血组织，无法明确区分胚系和体细胞mtDNA变异，以及无法评估与近期开发的免疫遗传学模型的互补性。在mtDNA标志物可用于供者选择或临床风险分层之前，需要进行功能学研究以及在更大、多祖先队列中的验证。
