细胞培养与分化
获得人类神经祖细胞及其分化后的神经元,用于后续ATAC-seq分析。
从14-21孕周胎儿脑组织分离原代人类神经祖细胞,培养并分化8周,经FACS分选获得EGFP+神经元。
Cell-type-specific effects of genetic variation on chromatin accessibility during human neuronal differentiation
包括体外细胞模型(原代神经祖细胞和分化神经元)、ATAC-seq组学分析、caQTL/eQTL关联分析、ASCA分析、基序破坏预测、精细定位、遗传力富集分析以及实验验证(荧光素酶报告基因实验)等独立证据层级。
原代人神经祖细胞(phNPCs) 分化后的神经元(经FACS分选的EGFP+神经元)
供体数量:92个;用于caQTL分析的祖细胞样本N=76,神经元样本N=61 分化条件:8周分化,使用AAV2-hSyn1-EGFP标记神经元 caQTL分析模型:线性混合模型(EMMAX),包含亲缘关系矩阵和10个基因型MDS成分作为固定效应 ATAC-seq峰值识别:CSAW以10bp bins,最小平均读数5,合并最小距离100bp 基序破坏分析:motifbreakR,相对熵变化判断破坏,阈值1e-4,仅保留effect="strong"
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
获得人类神经祖细胞及其分化后的神经元,用于后续ATAC-seq分析。
从14-21孕周胎儿脑组织分离原代人类神经祖细胞,培养并分化8周,经FACS分选获得EGFP+神经元。
测量祖细胞和神经元的染色质可及性,生成全基因组图谱。
对分离的细胞核进行ATAC-seq文库制备,然后进行测序。
从测序数据中识别开放染色质区域(峰)。
对ATAC-seq数据进行质量控制和读取比对,使用WASP去除映射偏差,过滤重复和线粒体读段,使用CSAW进行峰值识别。
鉴定祖细胞和神经元之间差异可及的调控元件。
使用limma对批量校正后的VST计数进行线性回归,鉴定差异可及峰。
识别影响染色质可及性的遗传变异(SNP)。
对每个峰,在±100kb范围内检测SNP与可及性水平的关联,使用EMMAX混合线性模型。
在同一位点检测等位基因特异的染色质可及性,确认caQTL的可靠性。
在杂合SNP处比较参考等位基因和替代等位基因的读数,使用DESeq2计算等位基因差异。
预测caSNP对转录因子结合基序的影响。
使用motifbreakR分析caSNP对已知TF基序的相对熵改变,鉴定可能破坏结合的SNP。
验证特定caSNP对增强子活性的影响。
将包含SNP的增强子片段克隆到荧光素酶报告载体(pGL4.23),在人类NPC中测定不同等位基因的荧光素酶活性。
确定caQTL是否与精神神经疾病和其他脑相关性状的GWAS位点共定位。
使用LD回归计算遗传力富集,并进行共定位分析。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 细胞培养 | 纤连蛋白 | Sigma-Aldrich | SLBN9832V |
| 多聚鸟氨酸 | Sigma-Aldrich | P3655-500MG | |
| 层粘连蛋白 | Life Technologies | 23017015 | |
| Neurobasal A培养基 | Life Technologies | 10888-022 | |
| Primocin抗生素 | Invivogen | ant-pm-2 | |
| BIT 9500补充物 | Stemcell Technologies | 09500 | |
| Glutamax添加剂 | Life Technologies | 35050061 | |
| 肝素 | Sigma-Aldrich | H3393-10KU | |
| 表皮生长因子 | Life Technologies | PHG0313 | |
| 成纤维细胞生长因子 | Life Technologies | PHG0023 | |
| 白血病抑制因子 | Life Technologies | PHC9481 | |
| 血小板衍生生长因子 | Life Technologies | PHG1034 | |
| B27补充物 | Life Technologies | 17504-044 | |
| 神经营养因子-3 | Life Technologies | PHC7036 | |
| 脑源性神经营养因子 | Life Technologies | PHC7074 | |
| 病毒标记 | AAV2-hsyn1-eGFP病毒 | UNC Vector Core | -- |
| 细胞消化 | 木瓜蛋白酶 | Worthington | -- |
| 脱氧核糖核酸酶 | Worthington | -- | |
| 细胞分选 | BD FACS Aria II流式细胞分选仪 | BD Biosciences | -- |
| Sony SH800S流式细胞分选仪 | Sony | -- | |
| 免疫细胞化学 | Pax6抗体 | BioLegend | Catalog#:901301 |
| Sox2抗体 | Cell Signaling Technology | Catalog#:3579 | |
| GFP抗体 | Millipore | Catalog#: AB16901 | |
| HOPX抗体 | Sigma-Aldrich | Catalog#:HPA030180 | |
| Nkx2.1抗体 | Millipore | Catalog#: MAB5460 | |
| Tbr2抗体 | eBioscience | Catalog#: 14-4877-82 | |
| Gad67抗体 | EMD Millipore | Catalog#: MAB5406 | |
| TUBB3抗体 | BioLegend | Catalog#: 801202 | |
| Satb2抗体 | Santa-Cruz | Catalog#: sc-81376 | |
| Ctip2抗体 | Abcam | Catalog#: ab18465 | |
| 山羊抗兔AF488二抗 | Thermo Fisher | Catalog#: A-11034 | |
| 山羊抗兔AF568二抗 | Thermo Fisher | Catalog#: A-11011 | |
| 山羊抗鸡AF488二抗 | Thermo Fisher | Catalog#: A-11039 | |
| 山羊抗大鼠AF647二抗 | Thermo Fisher | Catalog#: A-21247 | |
| 山羊抗小鼠AF488二抗 | Thermo Fisher | Catalog#: A-11001 | |
| 山羊抗小鼠AF568二抗 | Thermo Fisher | Catalog#: A-11031 | |
| ATAC-seq文库制备 | KAPA纯化磁珠 | Roche | -- |
| 测序 | Illumina HiSeq2500测序仪 | Illumina | -- |
| MiSeq测序仪 | Illumina | -- | |
| 荧光素酶报告实验 | pGL4.23载体 | Promega | -- |
| Q5定点突变试剂盒 | NEB | -- | |
| Lipofectamine STEM转染试剂 | Thermo Fisher | -- | |
| 双荧光素酶检测系统 | Promega | -- | |
| CLARIOstar Plus酶标仪 | BMG Labtech | -- |
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 细胞培养与分化 | 供体数量(N=92),培养代数(2.5±1.8 s.d.),分化时间(8周),病毒滴度(MOI 20000),接种密度(42,000 cells/cm2),培养基配方(如生长因子浓度) 阅读提示:Methods: Tissue acquisition and cell culture of phNPCs; Extended Data Figure 1a |
| ATAC-seq文库制备 | 细胞核输入数量(50,000),测序深度(平均25.5M),测序平台(HiSeq2500/MiSeq) 阅读提示:Methods: Library Preparation for human neural progenitors and neurons |
| ATAC-seq数据预处理 | 比对软件(bwa mem v0.7.17),WASP用于去除映射偏差,参考基因组(GRCh38/hg38),峰值识别软件(CSAW v1.16.1及参数) 阅读提示:Methods: ATAC-seq data pre-processing |
| caQTL分析 | 混合线性模型(EMMAX),协变量(基因型MDS、PCs),峰值窗口(±100kb),显著阈值(eigenMT-FDR < 0.05) 阅读提示:Methods: caQTL mapping |
| ASCA分析 | 每个SNP的读数阈值(总读数≥10,平均读数≥15),杂合供体数量(≥5),DESeq2模型 阅读提示:Methods: Allele specific chromatin accessibility |
| 基序破坏预测 | motifbreakR参数(阈值1e-4,effect="strong"),数据库(JASPAR2016),相对熵计算方法 阅读提示:Methods: Determining the impact of caSNPs on motifs |
| 功能验证(荧光素酶报告实验) | 载体(pGL4.23),转染量(120 ng/孔报告质粒,30 ng/孔海肾对照),细胞系(人NPC),培养时间(72小时),检测系统(Dual-Glo),供体数量(8个),重复数(至少3个重复/供体) 阅读提示:Methods: Luciferase reporter assay |