路易体痴呆、阿尔茨海默病与帕金森病痴呆的死后前额叶皮层转录组比较及LBD遗传证据整合分析

《International Journal of Molecular Sciences》：Analysis of Melatonin Secretion Patterns in Buffalo Milk and Its Association with AANAT Polymorphisms

　　

路易体痴呆（dementia with Lewy bodies，DLB）与阿尔茨海默病（Alzheimer’s disease，AD）及帕金森病痴呆（Parkinson’s disease dementia，PDD）共享大量临床和分子特征，使疾病特异性分子改变的识别变得复杂。本研究旨在表征这些疾病间的转录组差异，并将表达发现与已发表的路易体痴呆（Lewy body dementia，LBD）遗传证据整合。研究人员采用经年龄和性别校正的经验贝叶斯（empirical-Bayes）模型，配合错误发现率（false discovery rate，FDR）校正，分析了GSE150696中42例死后前额叶皮层样本（包括对照、AD、DLB和PDD）的基因水平表达谱；随后进行了敏感性分析、事后探索性候选基因筛选，以及与外部LBD多组学基因优先级排序数据集的整合。在DLB与AD的比较中，没有基因达到FDR<0.05；相比之下，SLC19A2、VAMP1和COX5B在DLB中表达低于PDD，并在主要年龄和性别校正模型中获得FDR支持，但其统计学显著性部分依赖于模型设定。一项事后跨疾病分析确定了60个探索性候选基因，其中SYT11是与已发表的LBD优先级基因集重叠的唯一基因；总的60个基因集并未显著富集LBD相关基因。这些发现未展示广泛稳健的疾病间差异，与大量分子重叠一致，同时也可归因于有限的统计功效。有限数量的转录组候选基因（包括探索性的SYT11重叠）需独立验证。

解读文章：

路易体痴呆（DLB）是一种以痴呆、认知波动、视幻觉和自发帕金森综合征为主要表现的进行性神经退行性疾病，其病理标志为α-突触核蛋白在路易小体中的异常聚集。DLB与帕金森病痴呆（PDD）共同构成路易体痴呆谱系，二者区分主要依赖痴呆相对于运动症状出现的时间关系；同时，DLB与阿尔茨海默病（AD）在临床和病理上存在大量重叠，包括混合病理和趋同的临床表现。这种复杂的临床病理重叠阻碍了疾病特异性分子分类、准确诊断和生物标志物开发。遗传学研究显示LBD具有异质性遗传架构，涉及突触核蛋白、溶酶体、脂质、免疫和神经退行性相关通路。已有整合性研究如Zhu等的大型LBD全基因组关联分析（GWAS）和多组学整合，优先排序了85个LBD相关基因，为转录组发现提供了外部框架。然而，直接比较DLB、AD和PDD三种疾病的转录组差异，并在严格多重检验控制下进行遗传证据整合的研究仍然有限。为此，研究人员对GSE150696数据集进行了基因水平重新分析，并将结果与已发表的LBD遗传证据和单核转录组背景相结合，以评估这些疾病之间的分子差异。

研究人员使用GSE150696中42例死后前额叶皮层（Brodmann 9区）样本，包括9例神经学正常对照、9例AD、12例DLB和12例PDD。通过基因符号水平的中位数概括，获得23,934个基因的表达谱。采用年龄和性别校正的经验贝叶斯线性模型，对五个诊断对比（DLB与AD、DLB与PDD、DLB与对照、AD与对照、PDD与对照）进行差异表达分析，并使用Benjamini–Hochberg错误发现率（BH-FDR）控制多重检验。随后进行敏感性分析、事后探索性候选基因筛选、与Zhu等85个LBD优先基因的整合，以及利用GSE178146单核转录组研究进行细胞类型背景解读。该框架严格区分了FDR支持的确认性发现与事后探索性候选，并进行了统计功效和前瞻性复制样本量评估。

主要技术方法包括：从GEO获取GSE150696标准化表达矩阵，经注释过滤保留399,440个外显子探针选择区域（PSR）特征，以基因符号为单位的样本内中位数概括得到基因水平表达谱；使用经验贝叶斯方差调节模型评估差异表达，并以BH-FDR<0.05为显著性阈值；通过诊断-only模型、排除低相关样本、增加死后间隔（PMI）协变量、非调节普通线性模型、表达过滤和单转录簇限制等方法进行敏感性分析；事后定义60个探索性跨疾病候选基因，并与外部LBD优先基因集进行Fisher精确检验；同时使用已发表的GSE178146单核转录组结果作为细胞类型背景参考，但未将其作为独立验证数据集。样本队列来源为GSE150696（42例）和GSE178146（28例）。

结果部分按以下小标题总结：

2.1 分析工作流程和发现队列：研究按预设流程从数据选择、基因水平预处理、质量控制、主差异表达分析、敏感性分析、事后探索候选、外部遗传整合、细胞类型背景到整合解释逐步推进。发现队列包含42例样本，各组年龄和性别分布已在文中列表描述。

2.2 全局表达结构和诊断组间差异表达：主成分分析显示四组样本在全局表达空间上广泛重叠，第一主成分解释26.86%方差，中位样本间Pearson相关系数范围为0.9533至0.9757。在主要年龄和性别校正模型中，DLB与AD对比无任何基因达到FDR<0.05；DLB与PDD对比中，SLC19A2（log₂FC=-0.670，FDR=0.02064）、VAMP1（log₂FC=-1.299，FDR=0.02064）和COX5B（log₂FC=-0.869，FDR=0.02064）在DLB中表达显著较低。DLB与对照对比中，1806个基因达到FDR显著，其中1499个下调；AD与对照对比中612个基因显著；PDD与对照对比则无显著基因。

2.3 跨疾病评估DLB与PDD候选基因：三个候选基因在DLB相对对照也表现出较低表达，且调整后p值均显著（SLC19A2 p=0.01193，VAMP1 p=0.04159，COX5B p=0.03871）。其中SLC19A2在DLB相对AD中显示名义显著（p=0.01506），但未通过转录组范围校正（FDR=0.7674）；VAMP1和COX5B在DLB-AD对比中无显著差异。

2.4 DLB与PDD关联的稳健性：SLC19A2、VAMP1和COX5B在诊断-only模型和排除低相关样本（GSM4556852）后仍保持FDR支持；在41例具有完整PMI信息的样本中，年龄性别校正模型下仍显著（FDR=0.04356）。但额外校正PMI后，调整p值升高至0.12095（SLC19A2和VAMP1）和0.18140（COX5B），方向不变；在未调节方差的普通线性模型中，三基因均未低于FDR阈值，但仍是DLB与PDD间最强信号（约0.052）。表达过滤和单转录簇限制未改变结论。

2.5 探索性跨疾病候选集：由于没有基因在DLB-AD和DLB-PDD两个对比中同时达到FDR显著且方向一致，研究人员定义了事后探索标准：在DLB-AD和DLB-PDD中均名义p<0.05、方向一致，并在DLB-对照中FDR<0.05，同时|log₂FC|≥0.25。最终获得60个探索性候选基因，其中7个高表达、53个低表达，包括SLC29A4、PTS、THOC7、ADSS、TUBB4B、KCTD8、NXPE2和SLC19A2。该集合被解释为假设生成，而非验证的分子特征。

2.6 与已发表LBD遗传证据的整合：将60个探索基因和3个主要候选基因与Zhu等85个LBD优先基因比较。85个基因中69个可映射到23,934基因背景（62个精确匹配、7个别名映射）。3个主要候选基因均不在85个基因中；SYT11是60个探索基因集中唯一与LBD优先基因重叠的基因。60基因集与映射后LBD基因集的Fisher分析p=0.1592，校正后p=0.3404，无显著富集。SYT11在外部框架中受GWAS、基因定位、基因集富集等多个组件支持，但在GSE150696中仅在DLB与对照比较时达到FDR显著（log₂FC=-1.371，FDR=0.00301），DLB-AD和DLB-PDD对比中仅名义显著。

2.7 GSE178146的细胞类型分辨转录组证据：已发表的单核转录组研究分析了28例前扣带皮层（7对照、7帕金森病、7 PDD、7 DLB），显示DLB和PDD具有广泛的分子相似性，直接比较有303个基因下调和87个上调；改变在兴奋性神经元和少突胶质细胞中尤为突出。该研究未被用作独立验证，因为其表达数据也纳入了Zhu等的多组学框架。

2.8 转录组和遗传证据整合摘要：研究形成两个互补证据类别：一是FDR支持的DLB与PDD差异基因（SLC19A2、VAMP1、COX5B），其效应方向和大小在多项敏感性分析中稳定，但正式显著性部分依赖协变量和方差模型设定；二是探索性SYT11重叠，该基因是唯一同时出现在转录组候选和外部LBD优先基因集中的基因，但60基因集整体未显著富集。

2.9 统计功效和前瞻性复制样本量估计：当前样本量（DLB 12 vs PDD 12）在α=0.05下只能对约Cohen’s d≥1.20的标准化效应提供80%功效，DLB vs AD（12 vs 9）对应d≥1.30。三个主要候选基因的观察标准化效应约为d=2.16-2.21；若按50%收缩进行保守规划，达到80%功效需每组19-20个样本，90%功效需每组24-25个样本。

讨论部分总结：有限的疾病间转录组分离与DLB和PDD密切的临床、病理及遗传关系一致。SLC19A2编码硫胺素转运体1，其下调可能涉及DLB相关皮层代谢改变；VAMP1是突触囊泡融合相关的SNARE蛋白，可能与突触囊泡循环异常有关；COX5B是细胞色素c氧化酶（复合体IV）的核编码亚基，其低表达可能与线粒体功能障碍有关。但bulk组织转录组不能区分细胞内在调控、神经元丢失或构成改变。SYT11与溶酶体功能和多巴胺释放相关，实验模型显示其过表达或缺失均可能产生有害效应，因此应谨慎解释。未发现60基因集显著富集外部LBD优先基因，支持将其保留为假设生成。敏感性分析强调结果对协变量和方差模型设定的依赖性。主要局限包括队列规模小、PMI数据缺失、未评估剪接异构体、探索性集合为事后确定等。

研究结论翻译：本研究确定了有限的候选分子差异，而非稳健的转录组范围疾病判别特征。SLC19A2、VAMP1和COX5B在主要年龄和性别校正模型中为FDR支持，但其统计学显著性依赖于模型设定。SYT11是60基因探索性转录组候选集与已发表LBD优先基因集重叠的唯一基因，该重叠未获得显著富集支持，需要独立验证。总体而言，研究未检测到广泛稳健的疾病间差异，这一结果与大量分子重叠一致，也与有限统计功效相容。更大规模的独立、神经病理学表征及细胞类型分辨的队列将有助于验证上述发现；保守的前瞻性计算表明，靶向复制三个主要DLB与PDD候选基因需要每组约20例（80%功效）或25例（90%功效）。未来整合单细胞多组学、蛋白质基因组学和可解释机器学习方法，可能有助于阐明细胞类型特异性的疾病机制并改善路易体相关痴呆的分子分层。