发现队列WGS关联分析
发现与WB-LM和A-LM相关的新遗传变异位点
对10,726个WB-LM和10,672个A-LM样本进行约45,000,000个SNP的荟萃分析，使用线性混合模型，先对年龄、年龄平方、性别、体重、身高、总脂肪及研究特异性协变量回归得到残差，再进行逆正态变换和全基因组评分检验。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
多族裔全基因组测序分析去脂体重
Multi-ancestry whole genome sequencing analysis of lean body mass
实验完整度
高包括大规模全基因组测序荟萃分析、罕见变异分析、生物信息学功能注释，以及果蝇和斑马鱼模型的功能验证，形成多层次证据。
主要模型
TOPMed和LOS队列WGS样本（10,729例） 英国生物银行（UKBB）复制样本 果蝇（Drosophila）肌肉特异性CG4984/EMP2敲低模型 斑马鱼（zebrafish）ssuh2基因敲除（F0）模型
重点核对
去脂体重测量方法：双能X射线吸收法（DXA） 全基因组测序覆盖度：TOPMed Freeze 8约30×，LOS约21×，UKBB约32.5× GWAS显著性阈值：p < 5×10^-8 罕见变异分析MAF ≤ 1%，使用SKAT检验和Bonferroni校正 功能验证方法：果蝇CG4984肌肉特异性敲低（Mef2-Gal4）和斑马鱼ssuh2多靶点CRISPR/Cas9敲除
摘要
背景：去脂体重是身体组成的重要生理组成部分。尽管去脂体重具有高遗传度，但迄今为止评估去脂体重遗传决定因素的研究主要局限于全基因组关联研究和常见变异。利用基于全基因组测序的研究，我们旨在多族裔的以人群为基础的队列中发现与去脂体重相关的新遗传变异。结果：我们描述了迄今为止最大规模的全基因组测序去脂体重荟萃分析，涵盖了来自六个TOPMed队列和路易斯安那骨质疏松研究队列的10,729个全基因组测序样本，并通过双能X射线吸收法测量。我们确定了七个先前GWAS未报道的与去脂体重显著相关的全基因组位点。我们在英国生物银行样本中部分复制了这些关联。在罕见变异分析中，我们发现了一个新的蛋白编码基因DMAC1，与女性全身去脂体重及四肢去脂体重相关，以及一个与男性四肢去脂体重相关的长链非编码RNA基因。这两个基因在骨骼肌组织中均表现出显著高表达。我们在动物模型中研究了两个新的去脂体重相关基因EMP2和SSUH2的功能作用。果蝇中EMP2缺陷导致活动能力显著下降，但肌肉组织或体脂形态未见改变，而斑马鱼中SSUH2基因突变刺激肌纤维生长。结论：我们的综合分析涵盖了大规模全基因组测序荟萃分析和功能研究，揭示了与去脂体重性状相关的新基因组位点和基因，为影响肌肉代谢和肌肉质量调控的途径提供了新见解。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 去脂体重（LM）的遗传决定因素是什么，尤其是通过全基因组测序能否发现新的常见和罕见遗传变异？
- 核心机制
- 文献未完全阐明核心机制，但提示EMP2和SSUH2可能参与肌肉代谢和肌肉质量调控，并通过影响肌肉功能或肌纤维生长发挥作用。
- 主要证据
- 大规模WGS荟萃分析（10,729样本）发现七个新位点；罕见变异分析发现DMAC1和ENSG00000273183；果蝇模型中EMP2同源物CG4984缺失导致活动能力下降；斑马鱼中SSUH2突变导致肌纤维面积增大。
- 研究意义
- 为理解肌肉代谢和肌肉质量调控的途径提供新见解，有助于识别肌肉减少症的药物靶点。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
条件分析和复制
验证发现位点的独立性和新颖性，并评估在独立样本中的可重复性
对p<1×10^-6的变异进行逐步条件分析，识别独立信号；同时基于先前GWAS信号进行条件分析，判断新颖性；在UKBB的WGS或SNP-Chip数据中对显著变异进行复制，并进行荟萃分析。
罕见变异关联分析
探索低频和罕见变异对去脂体重的累积效应
基于基因（错义/功能缺失变异，编码/非编码调控变异）、肌肉相关lncRNA、差异甲基化区域和转录起始位点等五种策略分组罕见变异（MAF≤1%），使用SKAT检验和Bonferroni校正。
生物信息学功能注释
识别潜在因果变异及其调控功能，以及相关基因和通路
使用LD>0.7筛选候选因果变异，通过ATAC-seq、DNase-seq、ChromHMM和Hi-C数据评估细胞类型特异性调控功能；使用Enrichr和GO/通路数据库进行富集分析；使用HaploReg和TRANSFAC预测转录因子结合位点。
果蝇模型功能验证
验证EMP2同源基因CG4984在肌肉功能中的作用
使用Mef2-Gal4驱动UAS-CG4984-RNAi实现肌肉特异性敲低，检测幼虫运动能力、蛹化高度和成虫飞行能力；同时观察肌肉形态和脂肪体脂质含量。
斑马鱼模型功能验证
验证SSUH2基因敲除对骨骼肌形态和功能的影响
利用CRISPR/Cas9多靶点技术生成ssuh2 F0敲除斑马鱼，通过H&E染色测量肌纤维面积和周长，通过油红O染色检查脂质含量，并评估游泳行为。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 去脂体重测量 | Hologic 4500W | Hologic | -- |
| Hologic QDR-2000 | Hologic | -- | |
| Lunar/Prodigy | GE-Lunar | -- | |
| Hologic QDR-4500 Discovery DXA扫描仪 | Hologic Inc. | -- | |
| Lunar-GE iDXA | GE-Lunar | -- | |
| 全基因组测序 | BGISEQ-500 | Beijing Genomics Institute (BGI Group) | -- |
| Illumina NovaSeq | Illumina | -- | |
| 基因分型 | Affymetrix UK Biobank Axiom芯片 | Affymetrix | -- |
| Affymetrix UK BiLEVE Axiom芯片 | Affymetrix | -- | |
| 果蝇实验 | 鬼笔环肽染色 | -- | -- |
| 尼罗红染色 | -- | -- | |
| 斑马鱼实验 | 油红O染色 | -- | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 队列样本 | 样本来源（TOPMed六个队列和LOS），样本量、性别比例、年龄等基本特征 阅读提示：Methods: Discovery cohorts: TOPMed and LOS; 表1 |
| 去脂体重测量 | DXA设备型号（Hologic 4500W, Hologic QDR-2000, Lunar/Prodigy, Hologic QDR-4500, Lunar-GE iDXA）及测量协议 阅读提示：Methods: Lean mass phenotype measurements |
| 测序与质控 | 测序平台、平均深度、比对和变异检出流程（e.g., GRCh38, GATK） 阅读提示：Methods: Sequencing data and quality control |
| 单变异关联分析 | 表型调整模型（性别、年龄、年龄平方、体重、身高、总脂肪等），MAF/MAC阈值，显著性阈值 阅读提示：Methods: Single-variant association analysis |
| 罕见变异关联分析 | 分组策略（基因、lncRNA等），MAF≤1%，SKAT检验，Bonferroni校正 阅读提示：Methods: Rare variants aggregated association analysis |
| 果蝇功能验证 | 果蝇品系（Mef2-Gal4, UAS-CG4984-RNAi, Lsp2-Gal4），敲低目标，行为学实验条件（幼虫运动、蛹化高度、负向地性），染色方法 阅读提示：Methods: Drosophila model; 图1及相关图注 |
| 斑马鱼功能验证 | 斑马鱼品系，CRISPR/Cas9多靶点基因敲除方法，F0代，肌纤维面积测量条件，油红O染色，游泳行为分析 阅读提示：Methods: Zebrafish model; 图2和图3 |