《International Journal of Molecular Sciences》：Integrative Genetic and Transcriptomic Analysis Reveals Candidate Molecular Differences Among Dementia with Lewy Bodies, Alzheimer’s Disease, and Parkinson’s Disease Dementia
编辑推荐：
路易体痴呆（DLB）与阿尔茨海默病（AD）和帕金森病痴呆（PDD）在临床和分子特征上存在大量重叠，这增加了识别疾病特异性分子改变的难度。本研究旨在表征这些疾病之间的转录组差异，并将表达发现与已发表的路易体痴呆（LBD）遗传证据进行整合。研究人员分析了GSE15
路易体痴呆（DLB）与阿尔茨海默病（AD）和帕金森病痴呆（PDD）在临床和分子特征上存在大量重叠，这增加了识别疾病特异性分子改变的难度。本研究旨在表征这些疾病之间的转录组差异，并将表达发现与已发表的路易体痴呆（LBD）遗传证据进行整合。研究人员分析了GSE150696中42例死后前额叶皮层样本的基因水平表达谱，包括对照组、AD、DLB和PDD，使用经年龄和性别校正的经验贝叶斯模型（empirical-Bayes models）并采用错误发现率（FDR）校正，随后进行敏感性分析、事后探索性候选基因筛选，以及与外部LBD多组学基因优先级排序数据集的整合。在DLB与AD的比较中，没有基因达到FDR < 0.05。相比之下，SLC19A2、VAMP1和COX5B在DLB中的表达低于PDD，并在主要年龄和性别校正模型中获得FDR支持，尽管其统计学显著性部分依赖于模型设定。一项事后跨疾病分析确定了60个探索性候选基因，其中SYT11是唯一与已发表LBD优先级基因集重叠的基因。总体60个基因集并未显著富集LBD相关基因。这些发现未显示广泛稳健的疾病间差异，与大量分子重叠一致，同时仍与有限的统计功效相容。少数转录组候选基因，包括探索性的SYT11重叠，值得独立验证。
论文解读
一、研究背景与问题
路易体痴呆（DLB）是一种进行性神经退行性疾病，以痴呆伴随认知波动、反复发作的视幻觉、快速眼动睡眠行为障碍和自发性帕金森综合征为主要临床特征。其神经病理学标志是错误折叠的α-突触核蛋白在路易小体和路易神经突中的积累，并常伴有阿尔茨海默型病理改变。这种临床病理复杂性对准确诊断、分子分层和疾病特异性生物标志物的开发构成了重大障碍。
帕金森病痴呆（PDD）与DLB在临床和神经病理学特征上高度重叠，两者通常被视为路易体痴呆谱系的一部分。其临床区分主要依赖于痴呆相对于帕金森综合征出现的时间顺序：DLB中痴呆先于或接近运动症状出现，而PDD中痴呆在已确立的帕金森病之后发展。这种操作性区分并不必然意味着两者是截然不同的分子实体。DLB的鉴别诊断还因与AD的重叠而进一步复杂化，AD主要以淀粉样蛋白-β和tau病理为特征。混合性AD和路易体病理常见，临床表型可能在疾病各阶段趋同。这些重叠使疾病特异性分子分类变得困难，也促使研究者开展直接的疾病间比较。
遗传学研究进一步支持LBD的异质性架构，涉及与突触核蛋白、溶酶体、脂质、免疫和更广泛的神经退行性生物学相关的共享及疾病特异性位点和通路。Zhu等人整合了大型LBD全基因组关联研究（GWAS）荟萃分析与功能基因组数据集，优先排序了85个候选LBD相关基因，为转录组学发现的优先排序提供了外部框架。然而，直接比较DLB、PDD和AD需要严格的多次检验控制和谨慎的生物学解释。
GSE150696数据集具有独特价值，因为它在一个数据集中同时包含了来自对照组和AD、DLB、PDD个体的死后前额叶皮层样本。既往研究虽已对该数据集进行了差异表达、通路、蛋白质-蛋白质相互作用和共表达网络分析，但采用了不同的分析流程，其结果应被视为互补而非直接复制或矛盾。
二、研究内容与结论
研究人员对GSE150696数据集进行了基因符号水平的重新分析，采用年龄和性别校正的经验贝叶斯模型、全转录组BH-FDR校正以及预设的敏感性分析。事后跨疾病候选基因与Zhu等人优先排序的LBD基因进行了比较，并利用GSE178146的细胞类型分辨证据进行了解读。该框架明确结合了严格的多次检验控制、模型稳健性评估、前瞻性功效和样本量评估、外部遗传优先排序以及细胞类型背景，同时保持了确证性证据与探索性证据之间的层级关系。
研究发现，在DLB与AD的直接比较中，没有基因达到FDR < 0.05，表明在该数据集内缺乏全转录组显著的单基因差异。在DLB与PDD的比较中，SLC19A2、VAMP1和COX5B三个基因在主要年龄和性别校正模型中获得了FDR支持，且均在DLB中表达较低。DLB与对照组的比较则显示出更广泛的转录组改变，共1806个基因达到FDR < 0.05。事后跨疾病分析确定了60个探索性候选基因，其中SYT11是唯一与已发表LBD优先级基因集重叠的基因，但总体60个基因集并未显著富集LBD相关基因。
这些发现表明，DLB与PDD之间的直接转录组分离有限，DLB与AD之间未检测到全转录组显著的单基因差异，这与大量分子重叠一致，同时也与有限统计功效相容。研究未检测到广泛稳健的疾病间差异，不能将阴性结果解释为生物学等同性。
三、关键技术方法
研究人员使用了以下主要技术方法：首先，从GEO数据库获取GSE150696数据集，包含42例死后前额叶皮层样本（9例对照组、9例AD、12例DLB、12例PDD），对沉积的标准化表达矩阵进行注释过滤和基因符号水平中位数汇总，生成23,934个基因的表达谱。其次，采用基因水平线性模型结合经验贝叶斯方差调节进行差异表达分析，以诊断组、年龄和性别为解释变量，评估五个成对对比，并使用Benjamini-Hochberg程序进行全转录组FDR校正。第三，进行多项敏感性分析，包括仅诊断模型、排除低相关性样本、完全病例分析、PMI额外校正、非调节普通线性模型、表达过滤和单转录簇限制。第四，将转录组发现与Zhu等人报告的85个LBD优先排序基因进行Fisher精确检验重叠分析。第五，利用GSE178146已发表的单核RNA测序研究提供细胞类型分辨的生物学背景。统计分析和功效计算使用Python 3.13.5和statsmodels 0.14.6完成。
四、研究结果
分析流程与发现队列。 研究基于GSE150696前额叶皮层数据集，包含42例死后样本。经注释过滤和基因符号水平汇总后，获得23,934个独特基因的表达谱，构成后续差异表达和整合分析的发现表达空间。
全局表达结构与跨诊断组差异表达。 主成分分析显示四个诊断组之间存在大量重叠，第一主成分解释26.86%的方差，第二主成分解释7.71%，未见清晰分离的疾病特异性聚类。在DLB与AD比较中，无基因达到FDR < 0.05，最小校正p值为0.5549。在DLB与PDD比较中，SLC19A2（log2FC = ?0.670，FDR = 0.02064）、VAMP1（log2FC = ?1.299，FDR = 0.02064）和COX5B（log2FC = ?0.869，FDR = 0.02064）三个基因获得FDR支持，均在DLB中表达较低。DLB与对照组比较显示最广泛的差异表达模式，共1806个基因达到FDR < 0.05，其中307个表达较高、1499个表达较低。AD与对照组比较识别出612个FDR支持基因，而PDD与对照组比较无基因达到FDR < 0.05。
DLB与PDD候选基因的跨疾病评估。 三个FDR支持基因在DLB与对照组比较中也均显示较低表达。SLC19A2在DLB与AD比较中显示名义显著的低表达（p = 0.01506），但未通过全转录组校正。VAMP1和COX5B在DLB与AD比较中均未显示显著差异。
DLB与PDD关联的稳健性。 三个候选基因在仅诊断模型和排除低相关性样本后仍保持FDR支持。在41例具有完整PMI信息的样本中，年龄和性别校正模型中三个基因仍显著（FDR = 0.04356）。但额外调整PMI后，校正p值分别增加至0.12095和0.18140，效应方向保持不变。在非调节普通线性模型分析中，三个基因均未低于FDR阈值0.05，但仍是DLB与PDD最强的信号。
探索性跨疾病候选基因集。 由于没有基因同时在DLB与AD和DLB与PDD比较中达到FDR < 0.05且方向一致，研究人员定义了事后探索性分析。要求基因在两个疾病间比较中均达到名义p < 0.05、方向一致，且在DLB与对照组比较中FDR < 0.05，识别出79个基因。进一步要求|log2FC| ≥ 0.25后，最终确定60个探索性候选基因，其中7个表达较高、53个表达较低。SLC19A2同时属于此集合和三个主要FDR支持基因之一。
与已发表LBD遗传证据的整合。 85个LBD优先排序基因中有69个可映射到发现表达空间。三个主要DLB与PDD候选基因均未与外部85个LBD优先排序基因集重叠。SYT11是60个探索性转录组候选基因与LBD遗传参考集之间唯一共享的基因，但60基因集与映射LBD基因集之间的重叠在Fisher精确检验中不显著（p = 0.1592）。SYT11在DLB相关比较中显示一致的较低表达，其中DLB与对照组比较达到全转录组显著（log2FC = ?1.371，FDR = 0.00301）。
来自GSE178146的细胞类型分辨转录组证据。 已发表的单核分析显示DLB和PDD具有大量整体分子相似性，但直接的DLB与PDD比较揭示了细胞类型分辨的差异，在报告的细胞类型特异性比较中，DLB相对于PDD有303个基因下调和87个基因上调。转录改变在兴奋性神经元群体和少突胶质细胞中尤为显著。
转录组与遗传证据的整合总结。 研究确定了两个互补的证据类别：SLC19A2、VAMP1和COX5B是主要年龄和性别校正DLB与PDD模型中FDR支持的三个基因，在DLB中表达较低；SYT11是60基因探索性转录组集合与已发表LBD优先排序基因集之间唯一共享的基因。
统计功效与前瞻性复制样本量估计。 在α = 0.05时，当前DLB与PDD比较（12对12）对Cohen's d ≥ 1.20的标准化效应提供80%功效。基于发现效应并采用50%收缩的保守规划假设，SLC19A2、VAMP1和COX5B的靶向独立复制研究在α = 0.0167下，80%功效需要每组约19至20个样本，90%功效需要每组约24至25个样本。
五、讨论与结论
讨论部分指出，DLB与PDD之间有限的转录组分离与其密切的临床、病理和遗传关系一致。DLB与PDD的区分主要基于痴呆与帕金森综合征的时间关系，而非截然不同的神经病理学标准。Zhu等人的LBD GWAS荟萃分析报告了LBD与AD和帕金森病均存在正向遗传相关，支持部分共享的病因架构。SLC19A2编码硫胺素转运蛋白1，与细胞能量代谢和神经代谢疾病相关；VAMP1是囊泡相关SNARE蛋白，参与突触囊泡融合和神经递质释放；COX5B编码细胞色素c氧化酶复合物IV的核编码Vb亚基，与线粒体功能障碍相关。SYT11在ATP13A2相关通路中调节溶酶体功能和自噬体降解，其表达水平与神经元功能的关系可能具有背景依赖性。
研究结论为：本研究识别了DLB、AD和PDD之间有限的候选分子差异集，而非稳健的全转录组疾病判别性特征。SLC19A2、VAMP1和COX5B在主要年龄和性别校正DLB与PDD模型中获得FDR支持，但其统计学显著性依赖于模型设定。SYT11是60基因转录组候选集与已发表LBD优先排序基因集之间唯一的探索性重叠，该重叠未获得显著富集支持，需要独立验证。总体而言，研究未检测到广泛稳健的疾病间差异，这一发现与大量分子重叠一致，但也与有限统计功效相容。需要更大规模、独立、经神经病理学表征和细胞类型分辨的队列来验证这些观察结果。保守的前瞻性计算表明，三个主要DLB与PDD候选基因的靶向复制需要每组约20名参与者（80%功效）或25名（90%功效）。未来整合单细胞多组学、蛋白质组基因组学和可解释机器学习方法的研究可能有助于解析细胞类型特异性疾病机制，并改善路易体和阿尔茨海默相关痴呆的分子分层。
