尼安德特人信息标记（NIM）的鉴定与关联分析
探究尼安德特人来源的遗传变异是否与PBC风险相关。
从1000 Genomes Project中提取235,592个NIMs，在CANUK队列中进行关联检验，并校正群体分层。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
与原发性胆汁性胆管炎风险相关的遗传风险变异遗传自尼安德特人
A Genetic Risk Variant Associated With the Risk of Primary Biliary Cholangitis Is Inherited From Neanderthals
实验完整度
中基于四个队列的基因型关联分析、RHE-mc遗传力估计、古代DNA时间动态分析及RegulomeDB功能注释，但无体内外功能验证实验。
主要模型
CANUK队列（4615例病例，9233例对照） OLD_IT队列（444例病例，901例对照） NEW_IT队列（255例病例，579例对照） MAYO队列（891例病例，621例对照）
重点核对
rs12531711在CANUK队列中p=2.04×10−26，OR=1.50 rs12531711在OLD_IT、NEW_IT和MAYO队列中复制，p值均<0.05，OR范围1.44-1.65 NIMs遗传力估计：h2=0.0019，SE=0.0075，富集度0.11±0.42，无显著贡献 古代DNA分析：青铜时代（5000-3000 BP）等位基因频率欧洲从2.0%升至7.7%，亚洲从3.4%升至17.4% RegulomeDB评分1f，提示强调控活性
摘要
背景：原发性胆汁性胆管炎（PBC）是一种多基因结构的自身免疫性胆管病，病因不明。有证据将尼安德特人来源的遗传变异与自身免疫性疾病联系起来；然而，它们对PBC易感性的贡献仍未被探索。目的：我们研究古老等位基因是否影响PBC风险。设计：分析了来自四个PBC队列的基因型数据：CANUK（4615例病例，9233例对照）、OLD IT（444例病例，901例对照）、NEW IT（255例病例，579例对照）和MAYO（891例病例，621例对照）。对235 592个尼安德特人信息标记（NIM）进行了关联检验，并校正了群体分层。使用RHE-mc估计遗传力，并使用来自欧洲和亚洲的古代DNA数据评估了时间单倍型动态。结果：TNPO3中一个尼安德特人来源的变异（rs12531711），此前由Cordell等人报道，在CANUK中显示出强关联（p = 2.04 × 10−26），并在所有验证队列中重复。该变异作为eQTL发挥作用，在血液中上调IRF5并下调TNPO3。它在非洲人群中罕见，但在欧洲人和亚洲人中常见。对携带rs12531711的单倍型的时间分析显示，其频率在青铜时代（5000–3000 BP）大幅增加，在欧洲从2.0%上升到7.7%，在亚洲从3.4%上升到17.4%，随后趋于稳定。尽管存在这种强单变异关联，但NIM解释的遗传力与非NIM变异没有显著差异。结论：一个尼安德特人来源的变异导致PBC风险，其频率在青铜时代增加，可能是由于病原体驱动的选择。古老基因渗入总体上对PBC遗传力的影响有限。需要功能研究来阐明rs12531711的作用。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 尼安德特人来源的等位基因是否影响原发性胆汁性胆管炎（PBC）的遗传风险？
- 核心机制
- rs12531711作为eQTL，上调IRF5表达并下调TNPO3表达，从而可能影响免疫调节和自身免疫风险。
- 主要证据
- 在四个欧洲血统队列中，rs12531711与PBC显著关联（CANUK p=2.04×10−26，验证队列p<0.05）；该变异位于TNPO3内含子，与IRF5和TNPO3表达相关；古代DNA显示该单倍型频率在青铜时代增加。
- 研究意义
- 揭示了PBC的进化起源，表明尼安德特人来源的变异可能通过免疫失调影响疾病风险，尽管对整体遗传力贡献有限，但为理解PBC遗传结构提供了新视角。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
rs12531711在多个PBC队列中的验证
验证rs12531711与PBC的关联是否可重复。
在OLD_IT、NEW_IT和MAYO三个独立队列中检验rs12531711与PBC的关联。
遗传力分析
评估尼安德特人变异对PBC遗传力的贡献。
使用RHE-mc方法在CANUK队列中估计NIMs的遗传力，并比较其与非NIM变异的差异。
功能注释
探究rs12531711的调控功能。
使用RegulomeDB v2对rs12531711进行功能注释，评估其调控潜力。
时间动态分析
追踪尼安德特人单倍型在欧亚人群中的频率变化。
利用古代DNA数据（AADR）分析rs10488631等位基因的频率随时间的变化，rs10488631与rs12531711完全连锁不平衡。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 基因型质量控制 | GWASTools | -- | -- |
| 遗传力分析 | RHE-mc | -- | -- |
| 连锁不平衡分析 | LDlink | -- | -- |
| 功能注释 | RegulomeDB | -- | -- |
| eQTL分析 | GTEx Portal | -- | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 关联分析 | 队列组成、样本量、遗传主成分数量、NIM主成分数量、显著性阈值 阅读提示：Methods 2.1-2.4，Table 1 |
| 遗传力分析 | 队列选择、MAF阈值、MHC区域排除、随机向量数、jackknife块数、校正主成分 阅读提示：Methods 2.5 |
| 时间动态分析 | 古代DNA样本数量、地理边界、时间分段、等位基因频率计算方法 阅读提示：Methods 2.6，Figure 3 |
| 功能注释 | RegulomeDB版本、评分标准 阅读提示：Results 3.2，Figure S2 |