温度胁迫下灵芝（Ganoderma lucidum）中HSF基因家族的全基因组鉴定与表达分析 胡金玉（Jinyu Hu） 李一红（Yihong Li） 吴少华（Shaohua Wu） 刘立伟（Liwei Liu） 周家伟（Jiawei Zhou） 李伟（Wei Li） 张瑞（Rui Zhang） 梁宗索（Zongsuo Liang） 杨东峰（Dongfeng Yang） 杨宗琪（Zongqi Yang）

《Genes》：The Genome-Wide Identification and Expression Profiling of the HSF Gene Family in Ganoderma lucidum Under Temperature Stress Jinyu Hu, Yihong Li, Shaohua Wu, Liwei Liu, Jiawei Zhou, Wei Li, Rui Zhang, Zongsuo Liang, Dongfeng Yang and Zongqi Yang

本研究通过全外显子测序分析247例早发性和家族性帕金森病患者及1622名对照的467个线粒体代谢相关核基因的罕见变异负担，发现mtDNA维持基因集的罕见变异负担在PD患者中显著增加，支持罕见变异累积影响PD易感性的假说。

　　
帕金森病（PD）的线粒体代谢相关遗传变异研究

帕金森病作为第二常见的神经退行性疾病，其病理机制涉及多维度生物学过程。近年研究聚焦于线粒体功能障碍在PD发生中的作用，本文通过大样本基因组学分析，系统考察了早发及家族性PD患者中与线粒体代谢相关的罕见遗传变异负担。

研究采用斯洛文尼亚、克罗地亚和塞尔维亚的247例PD患者与1622例对照组进行全外显子测序。样本选择严格遵循临床诊断标准，排除继发性PD及非神经退行性疾病患者。通过整合Genomics England的线粒体疾病诊断面板，确定467个与线粒体代谢密切相关的核基因作为研究对象。

在数据质量控制方面，采用双平台测序（Illumina NovaSeq与NextSeq）确保深度≥60×，并通过多维度质控（分子性别匹配、异质性比、覆盖率等）剔除低质量样本。最终获得6036例罕见错义突变和648例功能丧失性突变用于分析。

基因集负担分析显示，线粒体DNA维护相关基因集具有显著统计学差异（p<0.05）。该基因集涵盖mtDNA复制调控、动态平衡维持及核苷酸代谢等关键过程，包括TOP3A、FBXL4等14个已知与线粒体功能相关的基因。其中，mtDNA复制复合体核心基因如PNPT1（poly(A) polymerase I）和AARS2（α-氨甲酰磷酸合成酶II）的变异负担显著高于对照组。

单基因分析发现，在错义突变层面，NDUFB8（复合体I亚基）、SDHA（复合体II亚基）等8个基因在显性遗传模式下呈现趋势性关联；TOP3A（DNA拓扑异构酶III）和MRPL3（复合体I亚基）等6个基因在隐性遗传模式下达到显著性水平。值得注意的是，这些基因多涉及氧化磷酸化关键环节，其中NDUFB8与复合体I功能直接相关，SDHA参与三羧酸循环。

研究特别关注mtDNA维护机制。通过整合现有文献发现，PD患者中普遍存在mtDNA拷贝数减少及基因组稳定性下降。D?lle团队的前期研究揭示，纹状体神经元mtDNA拷贝数在PD患者中呈现年龄相关增长停滞，提示DNA维护系统受损。本研究的基因集分析进一步支持这一观点，TOP3A和FBXL4等基因的变异可能干扰mtDNA复制复合体的组装，导致基因组不稳定累积。

功能分析显示，PD相关变异多集中于线粒体翻译调控（AARS2、PUS1）、呼吸链组装（NDUFB8、ATPAF2）及代谢整合（SLC25A19、DLAT）三大模块。其中：
1. 翻译相关基因（PNPT1、AARS2）变异可能影响mtDNA编码蛋白的合成效率，导致复合体I亚基缺失；
2. 复合体组装基因（NDUFB8、SDHA）突变直接削弱呼吸链功能；
3. 代谢调控基因（SLC25A19、DLAT）变异影响丙酮酸代谢流，加剧能量危机。

研究创新性地将线粒体基因组维护机制纳入分析框架。通过比较PD患者与对照的mtDNA相关基因变异负担，发现TOP3A基因的隐性遗传模式显著相关（OR=1.89，p=0.023）。该基因编码的DNA拓扑异构酶III在mtDNA复制中起关键作用，其功能缺失可能引发复制压力，导致mtDNA缺失/重复突变累积。

值得注意的是，研究排除了已知PD致病基因（如SNCA、LRRK2等），聚焦于非编码区域及表观遗传调控相关的罕见变异。这种设计有效规避了已知致病基因的干扰，更精准地揭示线粒体代谢通路的潜在致病机制。

方法学上采用分层分析策略：首先通过基因集分析识别关键通路，再在特定通路内进行单基因负担检验。这种多层次分析方法有效提高了统计效力，特别是在样本量有限的情况下（n=247），仍能检测到具有临床意义的变异信号。

研究局限性主要在于人群来源单一（巴尔干地区），可能影响结果的泛化性。后续研究建议采用多中心、跨种族的大规模队列（目标样本量≥5000例）进行验证。此外，虽未检测到独立致病基因，但结合多基因风险评分系统，可更全面评估变异累积效应。

临床转化方面，研究提出三个方向：
1. 开发针对mtDNA维护的靶向疗法，如小分子激活剂（如mmitophagy促进剂）；
2. 构建基于多基因负担的生物标志物体系，用于早期诊断及疗效监测；
3. 探索线粒体-核基因组互作机制，建立整合型病理模型。

该研究为理解PD的线粒体致病机制提供了新视角，证实多罕见变异协同作用可能通过干扰mtDNA维护-呼吸链稳态轴，导致神经元能量代谢衰竭。后续研究可结合蛋白质互作组学及表观遗传修饰分析，深入解析变异如何改变线粒体代谢网络，为开发靶向治疗提供理论依据。