细胞培养与分化模型建立
建立并验证代表神经发生过程中两个关键阶段（神经祖细胞和神经元）的细胞模型，用于e/sQTL分析。
从14-21孕周胎儿脑组织分离原代人神经祖细胞（phNPC），培养并扩增，分化后通过AAV2-hSyn1-eGFP病毒标记神经元，并通过FACS分选。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
脑特征相关变异在神经发生过程中影响细胞类型特异性的基因调控
Brain-trait-associated variants impact cell-type-specific gene regulation during neurogenesis
实验完整度
高包含原代人神经祖细胞和神经元两种细胞类型，进行了eQTL、sQTL、ASE分析，并结合染色质可及性数据、共定位和TWAS进行机制验证，证据层级丰富。
主要模型
原代人神经祖细胞（phNPC） 分化后的神经元后代
重点核对
细胞来源：14-21孕周胎儿脑组织（推断12-19孕周后受孕） 分化条件：无生长因子分化5周，然后以MOI 20000的AAV2-hSyn1-eGFP转导，再分化3周 神经元分选：FACS分选（BD FACS Aria II或Sony SH800S），分化后56天（8周） RNA-seq建库：KAPA Stranded RNA-seq with Riboerase (HMR)试剂盒，起始RNA 50 ng eQTL/sQTL分析模型：线性混合模型（EMMAX），控制10个MDS、亲缘关系矩阵、及全局表达/剪接主成分
摘要
通过对基因表达和剪接数量性状位点（e/sQTL）分析来解释神经精神疾病和脑相关特征的非编码风险位点的功能，通常是在尸检后的成人bulk组织中进行。然而，遗传风险位点富集在新皮层分化过程中活跃的调控元件中，风险变异的调控效应可能被bulk组织中的异质性所掩盖。在这里，我们在代表两个主要发育阶段的培养细胞——原代人神经祖细胞（n = 85）及其分选后的神经元后代（n = 74）中绘制了e/sQTL和等位基因特异性表达，识别出许多在bulk发育中皮质壁或成人皮质中均未检测到的位点。通过共定位和基于转录组全关联分析的遗传插补，我们揭示了在新皮层分化过程中活跃的、支持脑相关特征风险的细胞类型特异性调控机制。具体而言，我们鉴定了一个祖细胞特异性的CENPW eQTL，与皮层表面积和受教育程度的常见变异关联共定位。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 常见遗传变异如何通过影响神经发生过程中的基因表达和剪接，进而影响脑相关性状和神经精神疾病风险？
- 核心机制
- 在神经祖细胞中，CENPW基因的eQTL变异rs9388486的C等位基因破坏了转录因子（CREM、ATF1、ATF2、ATF4）的结合基序，降低了启动子区域的染色质可及性，导致CENPW表达下降，从而可能通过影响神经祖细胞增殖或凋亡，增加皮层表面积和认知功能。
- 主要证据
- 利用85个供体的神经祖细胞和74个供体的神经元进行eQTL/sQTL分析，发现大量细胞类型特异性e/sQTL；通过共定位分析，CENPW eQTL与皮层表面积和受教育程度GWAS信号共定位；通过染色质可及性数据（caQTL）和转录因子基序分析，支持上述调控机制。
- 研究意义
- 该研究提供了细胞类型和时间特异性e/sQTL资源，有助于理解非编码风险变异在脑发育中的作用，为神经精神疾病的分子机制研究提供了基础，并可能指导未来治疗靶点的开发。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
转录组测序及数据预处理
获取祖细胞和神经元的全转录组表达和剪接数据，并控制技术因素和批次效应。
提取RNA，构建链特异性文库，进行RNA-seq测序；进行比对、定量、质量控制，并识别和校正技术协变量。
eQTL和sQTL定位
分别鉴定神经祖细胞和神经元中与基因表达和剪接相关的遗传变异。
使用线性混合模型（EMMAX）进行cis-eQTL和cis-sQTL分析，控制群体分层和亲缘关系；采用条件分析和分层多重检验校正。
细胞类型特异性分析
评估e/sQTL在祖细胞和神经元之间的共享程度和特异性。
通过LD-based overlap、m-values、π1统计量比较两个细胞类型之间的e/sQTL共享情况。
e/sQTL与GWAS共定位分析
识别与脑相关性状GWAS信号共定位的细胞类型特异性e/sQTL。
对多种GWAS（如精神分裂症、智力、教育程度、皮层表面积等）进行LD阈值共定位分析，采用条件分析验证共定位。
功能注释和调控机制解析
探索e/sQTL所在的功能基因组区域和潜在调控机制。
使用GARFIELD软件进行eQTL/sQTL在染色质状态和RNA结合蛋白结合位点的富集分析；使用motifbreakR预测转录因子结合基序的破坏。
TWAS分析
整合e/sQTL和GWAS，鉴定与脑相关性状相关的基因表达和剪接。
使用FUSION软件，基于基因表达和剪接权重，对多种脑相关性状进行转录组全关联分析（TWAS）。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
验证样本或模型特征
验证分子表达变化
验证剪接变化
验证功能注释
数据分析与统计
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 细胞培养 | AAV2-hSyn1-eGFP病毒 | -- | -- |
| 免疫荧光 | SOX2抗体 | Millipore | AB5603 |
| Ki67抗体 | Invitrogen | 14-5698-82 | |
| HOPX抗体 | Sigma-Aldrich | HPA030180 | |
| TUJ1抗体 | Biolegend | 801202 | |
| GFP抗体 | Millipore | AB16901 | |
| Alexa Fluor 568 | Invitrogen | A11036 | |
| Alexa Fluor 647 | Invitrogen | A21247 | |
| Alexa Fluor 488 | Invitrogen | A11001 | |
| RNA-seq | miRNeasy Minelute试剂盒 | QIAGEN | -- |
| Qubit 2.0荧光计 | Thermo Fisher Scientific | -- | |
| Agilent Tapestation | Agilent | -- | |
| KAPA Stranded RNA-seq with Riboerase (HMR)试剂盒 | Kapa Biosystems | -- | |
| FACS分选 | BD FACS Aria II | BD Biosciences | -- |
| Sony SH800S | Sony | -- | |
| 基因分型 | Illumina HumanOmni2.5 | Illumina | -- |
| Illumina HumanOmni2.5Exome | Illumina | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 细胞培养与分化 | 胎儿脑组织来源（孕周）、培养条件（生长因子、分化时长）、AAV病毒转导MOI、FACS分选时间点 阅读提示：详见Material and methods中的Cell culture部分 |
| RNA-seq | RNA提取试剂盒、建库试剂盒、起始RNA量、测序平台和读长 阅读提示：详见Material and methods中的RNA-seq library preparation和RNA-sequencing data processing部分 |
| eQTL/sQTL分析 | 分析模型（EMMAX）、协变量（MDS、PCs、FACS sorter）、cis窗口大小、多重检验校正方法 阅读提示：详见Material and methods中的Cell-type-specific local eQTL mapping和Splicing QTL mapping部分 |
| 共定位分析 | GWAS汇总统计来源、LD计算人群、共定位阈值（r2>0.8）、条件分析步骤 阅读提示：详见Material and methods中的LD-thresholded colocalization with brain disorders and traits GWAS部分 |
| TWAS分析 | 权重计算模型（FUSION）、LD估计人群、GWAS汇总统计准备 阅读提示：详见Material and methods中的TWAS analysis部分 |