《Frontiers in Genetics》：Epigenetic regulators are preferentially coordinated with protocadherin gene expression across the human brain: a genome-wide co-expression analysis
编辑推荐：
大脑如何构建和维持区分不同神经元的精确布线模式？簇状原钙黏蛋白（clustered protocadherins, PCDHs）——一种赋予每个神经元独特身份标签的细胞表面分子家族——是这一过程的核心。它们的表达受到一个复杂的、位点特异性的表观遗传系统的严格控
大脑如何构建和维持区分不同神经元的精确布线模式？簇状原钙黏蛋白（clustered protocadherins, PCDHs）——一种赋予每个神经元独特身份标签的细胞表面分子家族——是这一过程的核心。它们的表达受到一个复杂的、位点特异性的表观遗传系统的严格控制，该系统涉及DNA甲基化、CTCF（CCCTC结合因子）结合和染色质环化。然而，更广泛的表观遗传调控机器的活性是否与人脑中原钙黏蛋白的表达相协调，此前尚未被系统检验。在此，研究人员发现表观遗传调控因子在多个脑区与原钙黏蛋白优先共表达，表明存在比先前描述的位点特异性机制更广泛的转录协调。利用来自13个脑区的2642个脑样本的GTEx v8 RNA测序（RNA-seq）数据，研究人员进行了全基因组共表达筛选，并观察到在前额叶皮层中，表观遗传调控因子在PCDH协调基因前5%中富集了6.5倍（Fisher精确检验p = 2.8 × 10-10）。该富集在其他脑区独立重复，在多项敏感性分析中持续存在，在校正细胞类型组成后保持不变，并在一个独立脑库队列中重复。排名靠前的表观遗传调控因子汇聚于一组与明确特征性神经发育综合征相关的染色质重塑基因。这些发现通过将原钙黏蛋白的表观遗传协调背景从基因座本身扩展到与神经发育疾病相关的染色质重塑机制共享的更广泛转录程序，重新定义了原钙黏蛋白生物学。这一概念性进展与直接协同调控一致——尽管其本身并未建立直接协同调控关系——并确定了一组可检验的机制假设，用于解释神经发育障碍中破坏的染色质重塑活性如何可能传播到人脑中依赖于PCDH的神经元身份程序。
**论文解读：表观遗传调控因子与大脑原钙黏蛋白基因表达的协调性——一项全基因组共表达分析**
**研究背景与问题**
大脑中数十亿个神经元通过精确的突触连接形成复杂神经网络，而每个神经元需要独特的分子身份标签以实现自我识别和避免自身突触过度连接。簇状原钙黏蛋白（clustered protocadherins, PCDHs）基因家族编码约50种细胞表面分子，通过随机选择单个原钙黏蛋白亚型，为每个神经元生成独特的细胞表面分子条形码，从而驱动神经元自我识别、树突自我回避和突触特异性。该家族位于人类5q31染色体上的三个串联基因簇（PCDHA、PCDHB、PCDHG）中，其表达调控机制已被广泛研究：每个原钙黏蛋白基因拥有独立启动子，启动子选择由DNA甲基化差异调控——活性启动子呈低甲基化，沉默启动子呈高甲基化；CTCF（CCCTC结合因子）结合至未甲基化启动子区域，并与黏连蛋白复合体协同介导选定启动子与下游增强子元件（超级增强子/HS5-1区域）之间的长距离染色质环化，从而激活转录。这一复杂的表观遗传区室使PCDH基因座成为哺乳动物基因组中表观遗传调控最复杂的区域之一。
然而，尽管位点特异性调控机制在细胞系和小鼠模型中已得到充分表征，但个体间表观遗传调控因子（包括染色质重塑因子、组蛋白修饰酶和DNA甲基转移酶）的差异是否与人脑死后组织中PCDH表达水平共变，仍属未知。若存在共表达关联，则表明表观遗传机器与PCDH基因座调控之间的关系可能扩展到人脑群体水平的表达变异。此外，多个表观遗传调控因子（如ARID1A、EP300、CREBBP、CHD8、MECP2等）的失活突变导致智力障碍、自闭症谱系特征等神经发育综合征，这些疾病尽管涉及不同酶活性，却汇聚于重叠的神经发育表型。PCDH网络破坏是否构成这些不同表观遗传扰动的共同下游相关物，此前已被提出但未在群体水平系统检验。因此，该研究旨在通过无偏倚的全基因组共表达分析，检验表观遗传调控因子是否与人脑PCDH表达存在优先协调性。
**研究设计与结论**
研究人员利用GTEx v8 RNA-seq数据（涵盖13个脑区，共2642个脑样本），对16225个脑表达基因进行全基因组共表达筛选，计算每个基因与77个PCDH家族成员的平均Pearson相关系数，并检验表观遗传调控因子（52个脑表达基因）是否在PCDH协调基因前5%中富集。主要发现包括：前额叶皮层（BA9）中表观遗传调控因子富集6.5倍（Fisher精确检验OR=9.40，p=2.8×10
-10）；该富集独立复制于壳核、海马体和伏隔核；通过细胞类型反卷积校正六种主要脑细胞类型比例后，富集程度保持不变；在独立脑库队列（GSE80655）中重复验证；且排名靠前的表观遗传调控因子汇聚于神经发育综合征相关基因（如ARID1A、CHD8、EP300等）。该研究提供了群体水平转录组学证据，表明表观遗传调控因子与PCDH表达存在广泛协调，并提出神经发育疾病中染色质重塑扰动可能通过PCDH依赖性神经元身份程序传播的假设。论文发表在《Frontiers in Genetics》。
**关键技术方法**（不超过250字）
主要关键方法包括：1）全基因组共表达筛选：基于GTEx v8 RNA-seq数据（来自13个脑区，共2642个样本，样本队列来源为GTEx联盟），计算每个脑区内PCDH家族成员间及所有表达基因与PCDH的Pearson相关系数，并以平均相关系数作为PCDH协调评分进行基因排序。2）富集分析：采用Fisher精确检验评估52个表观遗传调控因子在前5%PCDH协调基因中的富集，同时进行无阈值预设的基因集富集分析（GSEA）验证。3）细胞类型反卷积：使用非负最小二乘（NNLS）法基于56个标志基因对六种脑细胞类型（神经元、星形胶质细胞、少突胶质细胞、小胶质细胞、内皮细胞、少突胶质细胞前体细胞）进行比例估计，并通过偏相关校正后进行重复分析。4）独立复制验证：在GSE80655脑库队列（281个样本，三个脑区）中重复共表达分析，并在GSE131706甲基化数据中检验PCDH与DNMT基因的CpG探针甲基化协调性。5）敏感性分析：包括RNA完整性数（RIN）分层、供体人口统计变量校正、阈值选择（1%、5%、10%）和条件特异性（精神疾病）检验。
**研究结果**
**3.1 原钙黏蛋白共表达是人脑区域的一致特征**
在所有13个GTEx脑区中，PCDH家族成员间的成对Pearson相关系数均为强阳性，总平均r=+0.497（区域范围+0.355至+0.595），且所有区域的subject-level bootstrap 95%置信区间排除零（如伏隔核r=+0.595，皮层r=+0.355）。通过比较41个非脑组织和全血（r=+0.233），发现PCDH共表达具有脑特异性（脑-血比2.1:1）。
**3.2 全基因组筛选识别表观遗传调控因子为BA9中PCDH协调基因的顶级成员**
BA9共表达筛选中，16225个脑表达基因的PCDH协调评分范围为-0.486至+0.620，前5%阈值对应r≥0.550。52个脑表达表观遗传调控因子中17个出现在前5%（预期2.6个），富集倍数6.5倍（OR=9.40，p=2.8×10
-10）。排名最高的基因为ARID1A（rank 2，r=0.617）、KDM3B（rank 5，r=0.611）、CHD8（rank 27，r=0.596）、CHD6（rank 33，r=0.594）和EP300（rank 54，r=0.591）。CTCF（PCDH基因座环化组织者）位于第1165位（r=0.540，92.8百分位），低于前5%阈值，但处于协调分布的尾部。富集在1%阈值下更强（13.5倍），GSEA确认无阈值富集（NES=2.27，FDR=0.001）。
**3.3 多区域复制富集**
在壳核（6.5倍）、海马体（3.85倍）和伏隔核（3.5倍）中独立进行全基因组筛选和Fisher精确检验，所有p值均<0.001。表观遗传基因排名在区域间高度一致（Spearman ρ=0.664-0.899，p<10
-7）。仅限28个簇状PCDH的分析显示富集7.3倍（OR=11.17，p=2.5×10
-12），确认结果非由非簇状PCDH驱动。
**3.4 对细胞类型组成和供体协变量的稳健性**
NNLS细胞类型反卷积估计BA9样本中神经元比例均值18.5%，与标志物SYN1（r=+0.749）、SNAP25（r=+0.883）、RBFOX3（r=+0.780）强相关。校正六种细胞类型比例后，富集仍为6.5倍（OR=9.40，p=2.8×10
-10），校正前后基因排名Spearman ρ=0.76-0.83。校正供体年龄、性别和RIN后富集为7.7倍（OR=12.14，p=1.6×10
-13）。RIN分层分析显示PCDH共表达在各四分位数间一致，无单调下降。
**3.5 汇聚于神经发育疾病基因**
排名靠前的表观遗传调控因子汇聚于神经发育综合征相关基因：ARID1A/ARID1B（Coffin-Siris综合征）、CREBBP/EP300（Rubinstein-Taybi综合征）、KMT2D（Kabuki综合征）、NSD1（Sotos综合征）、MECP2（Rett综合征）、KMT2A（Wiedemann-Steiner综合征）、SMARCA2（Nicolaides-Baraitser综合征）、KDM5C（X连锁智力障碍）、CHD8（自闭症谱系障碍）和ASH1L（自闭症谱系障碍）。这提示染色质重塑能力与PCDH依赖的神经元身份之间存在共享转录程序。
**3.6 跨数据集复制及跨条件保持**
GSE80655中三个脑区（前扣带回r=+0.594，伏隔核r=+0.569，DLPFC r=+0.576）的PCDH共表达模式与GTEx高度一致，基因对效应量跨数据集相关系数r=0.615-0.634（p<1×10
-50）。在双相障碍、重度抑郁症和精神分裂症中，PCDH共表达未发现显著变化（Δr范围-0.085至+0.103，所有9个比较FDR不显著），表明其构成性特征而非状态依赖性。GSE131706甲基化分析显示PCDH启动子与DNMT基因CpG探针间存在协调甲基化模式，与转录层协调一致。
**3.7 阴性发现**
SIGMAR1共表达分析未产生可比富集（77百分位）；在较小独立神经退行性疾病脑库队列中复制富集较弱，归因于样本量小及疾病相关变异。
**讨论与结论**
讨论部分指出，该研究发现了表观遗传调控基因与PCDH家族成员在人脑中的统计稳健共表达关联，该关联在四个独立脑区验证、通过无阈值GSEA确认、并在细胞类型反卷积及多项敏感性分析中保持稳健。所有发现为相关性，旨在生成可检验假设而非建立机制关系。细胞类型反卷积排除了粗细胞类型比例作为主要混杂因素，但单细胞验证仍必要。表观遗传调控因子在PCDH协调基因中的富集与已知PCDH基因座调控生物学一致，但共表达不能证明直接协同调控。神经发育疾病基因的汇聚提示这些基因可能通过PCDH网络破坏产生下游效应，但需在疾病相关模型系统中直接检验。PCDH共表达在精神疾病中完全保持，表明其更可能受严重细胞病理影响。RIN分层证明RNA质量不能解释脑库共表达差异，指向生物学原因。研究局限性包括：基于批量RNA-seq不能完全排除细胞类型残留混杂；共表达为相关测量；表观遗传调控因子列表可能不完整；甲基化分析在自闭症脑组织中进行，样本量小；GTEx队列为成人死后组织，发育期模式可能不同。
**结论部分翻译**：个体间原钙黏蛋白基因表达在人脑中的变异优先与表观遗传调控基因表达的变异共表达。该关联在四个脑区复制，在阈值敏感性分析（包括无阈值GSEA）、细胞类型比例校正、独立数据集复制和供体协变量校正中存活。Bootstrap置信区间确认所有13个GTEx脑区均存在正共表达，所有95%置信区间高于零。RIN分层表明仅RNA质量不能解释脑库共表达差异，指向生物学而非技术性解释。神经发育疾病基因在PCDH共表达排名顶部的汇聚，在脑区间一致，支持将表观遗传调控网络与PCDH表达联系起来的假设，值得实验研究。这些发现为未来机制研究提供了转录组学框架，并生成了关于遗传上不同的神经发育疾病突变共享下游后果的可检验假设。
