人类神经元分化过程中遗传变异对染色质可及性的细胞类型特异性效应

Cell-type-specific effects of genetic variation on chromatin accessibility during human neuronal differentiation

作者信息Dan Liang, Angela L Elwell, Nil Aygün, Oleh Krupa, Justin M Wolter, Felix A Kyere, Michael J Lafferty, Kerry E Cheek, Kenan P Courtney, Marianna Yusupova, Melanie E Garrett, Allison Ashley-Koch, Gregory E Crawford, Michael I Love, Luis de la Torre-Ubieta, Daniel H Geschwind, Jason L Stein
PMID34017130
发布时间2021-07
DOI10.1038/s41593-021-00858-w

实验完整度

包括体外细胞模型(原代神经祖细胞和分化神经元)、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"

摘要

精神神经疾病的常见遗传风险富集在皮层神经发生过程中活跃的调控元件中。然而,这些变异如何影响基因调控仍知之甚少。为了模拟常见遗传变异在人类皮层发育过程中对非编码基因组的功能影响,我们对来自92名供体的培养的人神经祖细胞及其分化后的神经元后代进行了ATAC-seq,并分析了染色质可及性数量性状位点。我们在988/1839个神经元/祖细胞调控元件中鉴定出显著的遗传效应,具有高度的细胞类型和时间特异性。约30%的caQTL还与基因表达变化相关。转录激活因子的基序破坏等位基因通常导致染色质可及性降低,而抑制因子的基序破坏等位基因则导致染色质可及性增加。通过整合细胞类型特异性caQTL和大脑相关的全基因组关联数据,我们能够精细定位并识别非编码精神神经疾病风险位点背后的调控机制。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
非编码基因组中的常见遗传变异如何影响人类皮层发育过程中的染色质可及性,并进而影响基因表达和神经精神疾病风险?
核心机制
非编码风险变异通过破坏转录因子结合基序,影响染色质可及性,进而调控基因表达,该机制具有细胞类型和时间特异性。
主要证据
基于ATAC-seq的caQTL映射发现数千个细胞类型特异性caQTL;约30%的caQTL与eQTL存在共享效应,且与eQTL效应方向一致;基序破坏分析发现激活子基序的破坏通常导致可及性降低,抑制子基序的破坏导致可及性增加;通过荧光素酶报告实验验证了多个具体位点(如rs185220、rs11544037、rs9930307)。
研究意义
该研究提供了人类皮层神经发生过程中遗传变异对基因调控影响的资源,有助于理解非编码变异在精神神经疾病风险中的功能,并有助于精细定位因果变异和细胞类型。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

细胞培养与分化

获得人类神经祖细胞及其分化后的神经元,用于后续ATAC-seq分析。

从14-21孕周胎儿脑组织分离原代人类神经祖细胞,培养并分化8周,经FACS分选获得EGFP+神经元。

2

ATAC-seq文库制备与测序

测量祖细胞和神经元的染色质可及性,生成全基因组图谱。

对分离的细胞核进行ATAC-seq文库制备,然后进行测序。

3

ATAC-seq数据处理与峰值识别

从测序数据中识别开放染色质区域(峰)。

对ATAC-seq数据进行质量控制和读取比对,使用WASP去除映射偏差,过滤重复和线粒体读段,使用CSAW进行峰值识别。

4

差异可及性分析

鉴定祖细胞和神经元之间差异可及的调控元件。

使用limma对批量校正后的VST计数进行线性回归,鉴定差异可及峰。

5

caQTL定位

识别影响染色质可及性的遗传变异(SNP)。

对每个峰,在±100kb范围内检测SNP与可及性水平的关联,使用EMMAX混合线性模型。

6

ASCA分析

在同一位点检测等位基因特异的染色质可及性,确认caQTL的可靠性。

在杂合SNP处比较参考等位基因和替代等位基因的读数,使用DESeq2计算等位基因差异。

7

基序破坏预测

预测caSNP对转录因子结合基序的影响。

使用motifbreakR分析caSNP对已知TF基序的相对熵改变,鉴定可能破坏结合的SNP。

8

功能验证(荧光素酶报告实验)

验证特定caSNP对增强子活性的影响。

将包含SNP的增强子片段克隆到荧光素酶报告载体(pGL4.23),在人类NPC中测定不同等位基因的荧光素酶活性。

9

与GWAS的整合分析

确定caQTL是否与精神神经疾病和其他脑相关性状的GWAS位点共定位。

使用LD回归计算遗传力富集,并进行共定位分析。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
纤连蛋白Sigma-AldrichSLBN9832V
多聚鸟氨酸Sigma-AldrichP3655-500MG
层粘连蛋白Life Technologies23017015
Neurobasal A培养基Life Technologies10888-022
Primocin抗生素Invivogenant-pm-2
BIT 9500补充物Stemcell Technologies09500
Glutamax添加剂Life Technologies35050061
肝素Sigma-AldrichH3393-10KU
表皮生长因子Life TechnologiesPHG0313
成纤维细胞生长因子Life TechnologiesPHG0023
白血病抑制因子Life TechnologiesPHC9481
血小板衍生生长因子Life TechnologiesPHG1034
B27补充物Life Technologies17504-044
神经营养因子-3Life TechnologiesPHC7036
脑源性神经营养因子Life TechnologiesPHC7074
AAV2-hsyn1-eGFP病毒UNC Vector Core--
木瓜蛋白酶Worthington--
脱氧核糖核酸酶Worthington--
BD FACS Aria II流式细胞分选仪BD Biosciences--
Sony SH800S流式细胞分选仪Sony--
Pax6抗体BioLegendCatalog#:901301
Sox2抗体Cell Signaling TechnologyCatalog#:3579
GFP抗体MilliporeCatalog#: AB16901
HOPX抗体Sigma-AldrichCatalog#:HPA030180
Nkx2.1抗体MilliporeCatalog#: MAB5460
Tbr2抗体eBioscienceCatalog#: 14-4877-82
Gad67抗体EMD MilliporeCatalog#: MAB5406
TUBB3抗体BioLegendCatalog#: 801202
Satb2抗体Santa-CruzCatalog#: sc-81376
Ctip2抗体AbcamCatalog#: ab18465
山羊抗兔AF488二抗Thermo FisherCatalog#: A-11034
山羊抗兔AF568二抗Thermo FisherCatalog#: A-11011
山羊抗鸡AF488二抗Thermo FisherCatalog#: A-11039
山羊抗大鼠AF647二抗Thermo FisherCatalog#: A-21247
山羊抗小鼠AF488二抗Thermo FisherCatalog#: A-11001
山羊抗小鼠AF568二抗Thermo FisherCatalog#: A-11031
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