MR 4.全基因组关联与孟德尔随机化联合分析送代码-CSDN博客

本文链接：https://blog.csdn.net/weixin_41368414/article/details/146168279

1.1 文章基本信息

标题
: Genome-wide association study unravels mechanisms of brain glymphatic activity
作者
: Huang, SY., Ge, YJ., Ren, P. et al.
期刊
: nature communications
影响因子
: 14+
发表日期
: 2025年1月
DOI
: 10.1126/sciadv.adr4606

1.2 文章摘要

根据沿血管周围间隙(ALPS)指数的弥散张量成像分析显示，脑淋巴活性与发育性神经精神和神经退行性疾病有关，但其遗传结构尚不清楚。在这里，我们从英国生物银行的31,021个个体的发现样本中鉴定出17个独特的全基因组显著位点和161个与ALPS指数相关的候选基因。七个基因座在两个独立的数据集中被复制。位于2p23.3位点的遗传信号在年轻和老年队列中产生了显著一致的效应。遗传相关和多基因重叠分析表明，ALPS指数、脑室容积和脑脊液TAU水平之间存在共同的潜在遗传机制，其中GMNC(3q28)和C16orf95(16q24.2)是共同的遗传基础。我们的研究结果增强了对ALPS指数遗传学的理解，并为进一步研究整个生命周期中淋巴清除活性的神经生物学机制及其与神经精神表型的关系提供了见解。

1.3 研究背景与科学问题

1.3.1 背景

类淋巴系统（Glymphatic System）
：近年发现的中枢神经系统废物清除系统，通过脑脊液-间质液循环清除代谢废物（如β-淀粉样蛋白），与神经退行性疾病（如阿尔茨海默病）密切相关。
研究空白
：类淋巴活动的遗传调控机制尚不明确，缺乏大规模遗传学研究。

1.3.2 科学问题

哪些遗传位点与类淋巴系统活动相关？
这些位点如何通过分子机制调控类淋巴功能？
类淋巴系统与神经退行性疾病的遗传关联性如何？

1.4 研究方法与技术路线

数据来源
：

样本量
：N=19542（如UK Biobank或定制队列）。
表型定义
：通过MRI技术（如动态对比增强成像/DCE-MRI）量化类淋巴活动。

GWAS分析

：

基因分型与质控
：标准GWAS质控流程（MAF>1%, HWE检验等）。
关联分析
：线性回归模型（调整年龄、性别、遗传背景等协变量）。
显著性阈值
：全基因组水平（P<5×10⁻⁸）。

功能注释与机制解析

：

功能基因组学
：使用GTEx、ENCODE等数据库进行eQTL/染色质互作分析。
通路富集
：GO、KEGG、Reactome通路分析。
跨疾病关联
：通过LDSC或MAGMA评估与神经退行性疾病的遗传相关性。

实验验证

：

小鼠模型
：靶向基因敲除/过表达，评估类淋巴活动变化（如荧光示踪剂清除率）。
细胞实验
：星形胶质细胞中候选基因的分子功能验证。

1.5 数据筛选

1.5.1 数据库整理

数据名称	用途描述	链接
GWAS汇总统计数据（ALPS-indexes）	类淋巴活动（ALPS-index）的全基因组关联分析（GWAS）结果数据	figshare 或 GWAS Catalog (GCP001081)
UK Biobank (UKB) 个体数据	类淋巴活动GWAS研究的原始基因型与表型数据来源（需申请访问）	UKB申请号19542
ABCD数据集	青少年脑认知发育研究的神经影像与遗传数据	ABCD
IMAGEN数据集	青少年神经行为与遗传关联研究的纵向数据	IMAGEN
HCP数据集	人类连接组计划的脑影像与遗传数据	Human Connectome Project
ADNI数据集	阿尔茨海默病神经影像学研究的临床与影像数据	ADNI
GWAS Catalog	存储全球GWAS研究结果的公共数据库	GWAS Catalog
GWAS Atlas	表型-基因关联的跨表型分析平台	GWAS Atlas
DSigDB数据库	药物-基因相互作用及靶点预测数据库	DSigDB
Human Protein Atlas	人类蛋白质组织表达谱数据库（用于验证基因的蛋白表达水平）	Human Protein Atlas

1.5.2 样品和表型描述

1.6 使用的软件与R包整理

名称	用途领域	类型	主要功能描述
PLINK	基因型数据处理	独立软件	全基因组关联分析（GWAS）基础工具，支持基因型质控、关联分析和基础统计计算。
GCTA	遗传力估计	独立软件	基于GREML方法估计复杂性状的遗传力，构建遗传关系矩阵（GRM）。
SAIGE	GWAS分析	R包	混合模型关联分析工具，适用于大规模队列数据，支持二元/连续表型和样本相关性校正。
ggplot2	数据可视化	R包	基于图形语法系统绘制高质量图表（曼哈顿图、QQ图等），支持高度自定义可视化。
locuszoomr	区域关联可视化	R包	生成LocusZoom风格区域关联图，整合SNP关联信号、基因注释和连锁不平衡（LD）信息。
metafor	Meta分析	R包	提供多种效应量模型和异质性检验，支持GWAS结果的Meta分析及森林图绘制。
coloc	共定位分析	R包	评估不同表型（如GWAS与eQTL）是否共享因果变异，支持贝叶斯共定位概率计算。
MAGMA	基因集富集分析	独立软件	将SNP关联信号映射到基因或通路水平，评估基因集与表型的遗传相关性（常通过R脚本调用）。
FUMA	功能注释	在线工具	自动化GWAS功能注释流程，整合eQTL、染色质互作和基因集富集分析。
TwoSampleMR	孟德尔随机化	R包	利用公开GWAS数据开展两样本孟德尔随机化分析，评估暴露与结局的因果关系。
SNPRelate	基因型数据处理	R包	高效处理大规模SNP数据，支持主成分分析（PCA）、亲缘关系计算和基因型质控。
GENESIS	遗传关联分析	R包	基于混合模型的GWAS框架，适用于复杂家系数据或存在群体分层的队列（如UK Biobank）。
qqman	可视化	R包	快速生成曼哈顿图和QQ图，支持全基因组显著性阈值标记和多染色体数据展示。
Bedtools	基因组区间操作	独立软件	处理BED/VCF文件，支持区间合并、交集和覆盖度分析，常用于功能基因组学数据整合。
ANNOVAR	变异注释	独立软件	对遗传变异进行功能注释（如基因区域、保守性、临床表型关联等），支持多数据库集成。

1.7 主要发现

目前的研究确定了17个独特的遗传风险位点和161个与ALPS指数相关的基因，ALPSs指数是一种有希望的淋巴清除活性标志物，研究了整个生命周期中年龄依赖和年龄独立的遗传效应，并证明了遗传与神经精神表型的关联。我们的发现提高了对ALPS指数遗传结构的理解，并为进一步研究整个生命周期中淋巴清除活性的生物学基础及其与大脑健康的相关性提供了线索。

1.8 研究成果

1.8.1 ALPS指数计算示意图

分析过程从初始的弥散加权成像图像到沿血管周围空间指数的计算分析。

1.8.2 与ALPS指数相关的遗传位点

1.8.3 基于SNP与平均ALPS指数和功能注释的关联

A.全基因组关联研究的曼哈顿图显示了UKB发现队列的平均ALPS指数。

B.基因组风险位点的大小和与平均ALPS指数相关的变异数量的概述。

C.饼状图显示LD中SNP的功能后果分布。

1.8.4 遗传对ALPS指数的影响及整个生命周期年增率

A.固定效应荟萃分析基因座2p23.3 rs4665946对平均ALPS指数的影响。

B-C.采用线性和二次元模型对SNP效应按年龄进行Meta回归。

D.森林图显示rs6012259和rs12146713在位点20q13.12和12q23.3的遗传效应。

1.8.5 鉴定的ALPS指数位点的遗传相关性

A.平均ALPS指数与脑容量之间的遗传相关性。

B.火山图显示了平均ALPS指数与白质微结构之间的遗传相关估计。

C.气泡图显示了三个ALPS指数与假定的危险因素、认知措施和其他脑衰老的MRI标记物之间的遗传相关性。

1.8.6 ALPS指数与神经精神和神经疾病及生物标志物的多基因重叠

A.通过condFDR分析，增强了对神经精神和神经疾病和生物标志物的三个ALPS指数的遗传位点的发现。

B.精神分裂症、脑脊液pTau、HIP-PVS、WMH条件下平均ALPS指数的条件Q-Q图。

1.8.7 通过GWGAS和基因定位策略鉴定的ALPS指数相关基因

A.GWGAS在平均ALPS指数上的曼哈顿图。

B.Venn显示了与三个ALPS指数中的任何一个相关的基因数量。

C.Upset图显示了三个ALPS指数中独特和共有基因的数量。

1.9 研究亮点

创新性
：首个大规模GWAS揭示类淋巴系统遗传基础。
多组学整合
：从遗传位点到分子机制、动物模型的系统性解析。
临床转化潜力
：为AD等疾病的治疗靶点开发提供新方向。

1.10 局限性与未来方向

1.10.1 局限性

表型测量限制
：MRI技术可能无法完全捕捉类淋巴活动的动态变化。
人群偏倚
：样本主要来源于欧洲血统，需跨族群验证。

1.10.2 未来方向

开发更高时空分辨率的类淋巴活动成像技术。
探索基因-环境交互作用（如睡眠、衰老的影响）。
靶向候选基因的干预研究（如小分子药物或基因治疗）。

1.11 临床意义

早期诊断
：类淋巴活动相关遗传标志物或助力AD风险预测。
治疗靶点
：GMNC和C16orf95调控可能成为改善类淋巴功能的新策略。
个性化干预
：基于遗传风险的睡眠调节或药物干预方案。

1.12 总结

本研究通过GWAS结合多组学分析，揭示了类淋巴系统活动的遗传调控网络，并建立了其与神经退行性疾病的分子联系，为理解脑废物清除机制及开发新型疗法提供了重要依据。

2 Reference

Huang, SY., Ge, YJ., Ren, P. et al. Genome-wide association study unravels mechanisms of brain glymphatic activity. Nat Commun 16, 626 (2025)关注桓峰基因发消息'ALPS'下载链接。