构建肥胖相关细胞类型的调控图谱
在脂肪细胞和下丘脑神经元分化过程中获取染色质可及性、基因表达和增强子-启动子相互作用图谱，为解读GWAS变异提供注释基础。
将人iPSC分化为下丘脑神经元，收集D12、D16、D27三个时间点；将SGBS前脂肪细胞分化为成熟白色脂肪细胞，收集四个关键时间点；对各时间点进行ATAC-seq、RNA-seq和原位启动子捕获Hi-C（cHi-C）。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
功能基因组学流程在肥胖相关GWAS位点中鉴定多效性和跨组织效应
A functional genomics pipeline identifies pleiotropy and cross-tissue effects within obesity-associated GWAS loci.
实验完整度
高多层级证据：体外细胞模型分化时间序列、MPRA功能筛选、cHi-C、CRISPR缺失与CRISPRi功能验证
主要模型
人iPSC来源的下丘脑神经元分化时间序列 SGBS人前脂肪细胞及成熟白色脂肪细胞分化时间序列 MPRA细胞系（SGBS、3T3-L1、GT1-7、HT22） 人iPSC（NA19101）及增强子缺失克隆
重点核对
MPRA中2396个变异在五种细胞系中的增强子活性及等位基因差异（双侧Mann-Whitney U检验，q<0.05） cHi-C在脂肪细胞和神经元分化多个时间点的增强子-启动子相互作用（CHiCAGO score>5） CRISPR-cas9删除rs2650492或rs9972768所在增强子后，iPSC分化四个时间点RNA-seq检测SBK1等基因表达 HEK293t中CRISPRi靶向rs2650492增强子后RNA-seq检测cis基因表达变化 荧光素酶报告基因验证EMVar等位基因特异性增强子活性
摘要
全基因组关联研究（GWAS）已鉴定出许多疾病相关变异，但这些关联背后的机制仍不清楚。为了理解肥胖相关变异，我们在脂肪细胞和下丘脑神经元中跨细胞分化阶段生成基因调控注释。然后，我们使用大规模并行报告基因检测（MPRA）测试了97个肥胖相关位点中的变异，并鉴定出具有细胞类型特异性或跨组织增强子调节特性的推定因果变异。将这些变异与基因调控信息整合，提示了肥胖GWAS关联背后的基因。我们还研究了16p11.2上一个复杂的基因组区间，其中两个独立位点表现出兆碱基范围的跨位点染色质相互作用。我们证明这两个位点内的变异调控一个共享的基因集。总之，我们的数据支持一个模型：GWAS位点包含的变异可改变跨组织的增强子活性，可能具有时间限制性效应，从而影响多个基因的表达。这种复杂模型对当前理解GWAS的努力具有广泛意义。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 肥胖相关GWAS位点中的非编码变异如何通过调控增强子活性影响靶基因表达，以及这些调控是否具有细胞类型特异性和时间特异性。
- 核心机制
- GWAS位点包含多个功能性变异，可改变增强子活性，通过长距离染色质相互作用调控一个或多个靶基因，效应可具有细胞类型特异性或跨组织多效性，并可能受时间窗口限制。
- 主要证据
- 在脂肪细胞和iPSC来源下丘脑神经元分化时间序列中生成ATAC-seq、RNA-seq和cHi-C图谱；MPRA在97个位点中鉴定出94个EMVar，其中多个位点含多个EMVar；CRISPR删除16p11.2两个独立位点的增强子后SBK1表达在特定分化阶段降低；CRISPRi靶向rs2650492增强子降低NUPR1表达。
- 研究意义
- 提供了在疾病相关细胞类型中功能解读GWAS位点的流程，支持GWAS位点常包含多个功能性变异并通过跨组织、时间限制性增强子调控多个基因的模型，对理解GWAS机制具有广泛意义。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
MPRA筛选肥胖GWAS位点的功能性变异
鉴定97个肥胖GWAS位点中与lead SNP高度连锁的变异是否具有增强子活性及等位基因特异性增强子调节效应。
选择97个独立肥胖GWAS位点的lead SNP及CEU人群中MAF≥5%、r2>0.8的高连锁双等位SNP共2396个；合成175 bp片段，每个等位基因连接18-19个独特10 bp条形码，克隆入pMPRA1载体；转染五种细胞系并进行MPRA测序分析。
整合cHi-C与eQTL将功能性变异分配至靶基因
利用三维基因组相互作用和eQTL证据，将EMVar与候选靶基因连接，并建立基因支持等级。
将EMVar位置与脂肪细胞和神经元分化时间序列cHi-C数据交集，定义EMVar与启动子相互作用（距离≤1 kb）；整合GTEx脂肪和脑eQTL数据；根据证据强度将基因分为I-IV类。
CRISPR-cas9删除16p11.2位点增强子并评估基因表达
验证16p11.2两个独立GWAS位点中的增强子是否调控共享基因，并确定效应的时间窗口。
在人iPSC（NA19101，非风险单倍型纯合）中分别删除含rs2650492的750 bp区域（SBK1 3'UTR内）和含rs9972768的1.3 kb区域（TUFM与SH2B1之间）；每个缺失系获得四个独立克隆，分化为下丘脑谱系，收集四个早期发育时间点（iPSC、腹侧化细胞、神经元前体、下丘脑前体）进行RNA-seq。
CRISPRi靶向rs2650492增强子评估其他细胞类型中的调控
在神经元谱系之外验证rs2650492增强子是否可调控其他基因，探究多效性。
在HEK293t细胞中使用dCas9-KRAB和四条靶向rs2650492区域的guide RNA进行CRISPRi，FACS分选双阳性细胞，提取RNA进行RNA-seq，检测SBK1-ATP2A1位点内cis基因表达变化。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 细胞培养 | DMEM/F12 (1:1)培养基 | Life Technologies | 11330-032 |
| 青霉素-链霉素溶液（10,000 U/ml） | Gibco | 15140122 | |
| 脂肪细胞分化 | 泛酸 | Sigma | B4639 |
| 生物素 | Sigma | B4639 | |
| 人转铁蛋白 | Sigma | T2252 | |
| 3-异丁基-1-甲基黄嘌呤 | -- | -- | |
| 罗格列酮 | -- | -- | |
| iPSC培养 | Matrigel hESC合格基质 | Corning | 354277 |
| mTeSR1培养基 | Stemcell | 85850 | |
| 下丘脑神经元分化 | 脑源性神经营养因子 | R&D Bioscience | 248‐BD |
| B27添加剂 | Thermo Fisher | 17504044 | |
| RNA提取 | RNeasy试剂盒 | Qiagen | 74104 |
| RNA-seq文库构建 | NEB Next Ultra II定向RNA-seq试剂盒 | NEB | -- |
| MPRA转染 | Lipofectamine LTX & Plus试剂 | Invitrogen | 15338100 |
| Polyplus jetPEI DNA转染试剂 | Polyplus | 101-10N | |
| MPRA克隆 | pMPRA1载体 | addgene | 49349 |
| Gibson组装预混液 | NEB | E2611S | |
| T4 DNA连接酶 | NEB | B0202S | |
| MegaX DH10B T1R电感受态细菌 | Invitrogen | C640003 | |
| MPRA mRNA分离 | Invitrogen Dynabeads | ThermoFisher | 61006 |
| Promega RQ1 DNA酶 | Promega | M6101 | |
| MPRA RNA定量 | Promega QuantiFluor RNA系统 | Promega | E3310 |
| MPRA cDNA处理 | RNA酶A | Invitrogen | 12091-021 |
| RNA酶T1 | ThermoFisher | EN0541 | |
| MPRA文库扩增 | Q5热启动高保真2×预混液 | NEB | M0494S |
| MPRA文库纯化 | Agencourt AMPure XP磁珠 | Beckman Coulter | A63882 |
| 荧光素酶报告基因检测 | 双荧光素酶报告基因检测系统 | Promega | E1910 |
| cHi-C交联 | 甲醛 | Sigma | 252549 |
| cHi-C淬灭 | 甘氨酸 | Sigma | G8898 |
| CRISPR-cas9编辑 | crRNA-Atto550和cas9蛋白 | IDT | -- |
| CRISPRi | dCas9-KRAB载体 | -- | -- |
| eGFP基因上游的guide载体 | -- | -- | |
| 流式分选 | BD FACS Aria Fusion 5–18流式细胞仪 | BD | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| MPRA变异筛选与增强子活性检测 | 变异纳入标准（CEU MAF≥5%，r2>0.8）、五种细胞系及转染试剂、每个等位基因条形码数量18-19、转染后48小时收集 阅读提示：Methods中Massively parallel reporter assay小节及Figure 2图注 |
| cHi-C与ATAC-seq时间序列 | 下丘脑神经元分化时间点D12/D16/D27及脂肪细胞四个时间点、每个时间点技术重复数、cHi-C显著性阈值CHiCAGO score>5 阅读提示：Methods中In situ promoter capture HiC及ATAC-seq小节，Figure 1图注 |
| CRISPR-cas9增强子缺失 | iPSC细胞系NA19101、删除区域大小（750 bp和1.3 kb）、每个缺失系独立克隆数（四个）、分化四个时间点及RNA-seq克隆重复数 阅读提示：Methods中CRISPR–cas9 editing of iPSC cells小节及Figure 5图注 |
| CRISPRi靶向rs2650492增强子 | HEK293t细胞、四条guide RNA、FACS分选双阳性细胞、独立实验次数n=4-5、Bonferroni校正阈值p<0.05/14 阅读提示：Methods中HEK293t CRISPRi小节及Figure 5d图注 |
| 荧光素酶验证EMVar等位基因效应 | 测试区域大小（约1 kb或470 bp上下文）、细胞系（SGBS前脂肪细胞和/或HT22）、独立实验次数（n=3-7）、内参共转染比例 阅读提示：Methods中MPRA enhancer validation luciferase assays及ATP2A1-SBK1 EMVar luciferase assays小节，Figure 4d图注 |