《Molecular Psychiatry》：Multi-omic data integration improves the resolution of the molecular etiology of autism in a mouse model
编辑推荐：
自闭症谱系障碍（ASD）是一种具有复杂分子病因学的多因素神经发育障碍。由于遗传原因占ASD病例的不到50%，因此需要新的方法来克服表征基因和分子通路（将风险等位基因与表型联系起来）的挑战。研究人员使用RNA测序、基于蛋白质的三维染色质构象（Hi-ChIP）和全
自闭症谱系障碍（ASD）是一种具有复杂分子病因学的多因素神经发育障碍。由于遗传原因占ASD病例的不到50%，因此需要新的方法来克服表征基因和分子通路（将风险等位基因与表型联系起来）的挑战。研究人员使用RNA测序、基于蛋白质的三维染色质构象（Hi-ChIP）和全基因组DNA甲基化测序方法，研究来自ASD小鼠模型（Cntnap2敲除[Cntnap2 KO]）的海马组织，以确定多组学数据整合是否能提高对导致复杂ASD表型的关键分子通路的分辨率。每个单独的组学数据集都识别出众多基因和通路的紊乱，提供了广泛的ASD相关见解，例如：1）699个下调基因，富集了神经元本体论术语；2）Cntnap2 KO小鼠中独特的染色质相互作用；3）大多数差异甲基化基因（1220/1659）具有神经元功能。多组学数据整合降低了异质性，并将数据精炼至43个与ASD相关的基因（例如Zbtb18、Cttnbp2、Gabbr2）。对这43个基因的通路分析识别出一个基因本体论术语：神经元突触可塑性的调节。这些发现与人类ASD死后脑组织的大规模基因表达研究一致，表明多组学数据整合可用于实现对异质性ASD分子病因学的更高分辨率。
**论文解读：多组学数据整合揭示自闭症小鼠模型中分子病因学的精细分辨率**
**研究背景与问题**
自闭症谱系障碍（ASD）是一种复杂的异质性神经发育疾病，其遗传病因涉及约1000个风险基因，但已知遗传因素仅能解释不到50%的病例。大量死后脑组织转录组研究提示突触功能与免疫通路异常，而表观遗传修饰（如DNA甲基化、组蛋白乙酰化）及非编码区调控元件（如增强子）的长程染色质相互作用在神经发育中起关键作用。然而，单一组学数据（如转录组、甲基化组或染色质构象）各自揭示了广泛的分子紊乱，却难以从异质性中精确识别核心致病通路。因此，研究人员提出假设：整合多组学数据（转录组、染色质构象和DNA甲基化组）能提高对ASD复杂分子病因学关键通路的分辨率。本研究利用已建立的Cntnap2敲除（Cntnap2 KO）小鼠模型——该模型表现出ASD样表型（重复行为增加、社交障碍）及神经病理特征（皮质神经元迁移缺陷、纹状体GABA能中间神经元减少）——对海马组织进行多组学分析，以验证上述假设。论文发表在《Molecular Psychiatry》。
**主要关键技术方法**
1. **RNA测序（RNA-seq）**：对3月龄雌性Cntnap2 KO与野生型（WT）小鼠（每组N=3）的海马组织进行poly-A富集文库构建及Illumina NovaSeq双端150bp测序，使用DESeq2进行差异表达分析（FDR<0.05）。
2. **蛋白质中心染色质构象捕获（Hi-ChIP）**：另取一批小鼠（每组10只，海马组织混合后约350mg），分别用抗RNA聚合酶II（Pol II）和抗组蛋白H3赖氨酸27乙酰化（H3K27ac）抗体进行免疫沉淀，然后进行邻近连接和测序；使用FitHiChIP识别显著染色质相互作用（FDR<0.05），并筛选两种抗体均富集的“可重复染色质相互作用”。
3. **全基因组DNA甲基化测序**：使用同一批小鼠的DNA，通过dmrseq识别差异甲基化区域（DMR，>5%甲基化差异，置换P<0.05）。
4. **多组学整合**：将差异表达基因、基因特异性染色质相互作用（特别是启动子-增强子相互作用，PEI）及DMR相关基因进行重叠分析，最终获得43个在三个组学层面均受扰动的基因。
**研究结果**
**1. Cntnap2 KO小鼠中ASD相关基因下调**
RNA-seq鉴定出1440个差异表达基因（741个上调、699个下调）。下调基因显著富集于神经元本体论术语（如神经元发育、突触传递正调控），且其中ASD相关基因（SFARI数据库，综合征及基因评分1级）占比8.3%（58/699），包括Grin2a、Nrxn3、Reln等。与人类ASD死后脑组织转录组比较，发现237个（前额叶皮质）和250个（海马）重叠基因，86个为三向重叠，提示小鼠与人类ASD转录变化存在趋同。
**2. Cntnap2 KO海马中染色质互作景观改变**
Hi-ChIP分析在Cntnap2 KO小鼠中鉴定出8985个可重复染色质相互作用（WT为7439个）。注释到基因后，发现177个高置信度ASD相关基因在Cntnap2 KO中具有更高的互作强度（累积FDR和归一化接触计数均显著高于WT，P<0.001），表明ASD相关基因的染色质构象广泛失调。
**3. Cntnap2 KO特异性启动子-增强子相互作用（PEI）相关基因**
筛选位于基因启动子区域（转录起始位点TSS周围3个5kb bin）的PEI，发现Cntnap2 KO特异性PEI相关基因（4217个）富集于神经系统发育术语（如神经发生、轴突发生），而WT特异性PEI相关基因（2750个）富集于细胞代谢过程。转录因子（TF）结合基序分析显示，Cntnap2 KO特异性PEI的启动子区域富集了13个与神经元过程及ASD相关的TF基序（如NRF1、E2F3、Sp1、Egr1）。
**4. Cntnap2 KO特异性PEI与已知ASD基因下调相关**
整合转录组与PEI数据，发现600个差异表达且PEI基因型特异性的基因（317个下调）。多数高置信度ASD相关基因（如Grin2a、Nlgn2、Gabbr2）在Cntnap2 KO中表现出PEI强度增加伴随表达水平降低，与下调基因的神经元本体论富集一致。
**5. DNA甲基化改变与差异表达相关**
同一批小鼠的甲基化组与转录组比较发现174个差异表达基因包含DMR，其中78%（130/166）表现为高甲基化伴随低表达，通路分析富集于突触组织、神经元分化等神经过程。
**6. 多组学数据整合优化分子分辨率**
将DMR注释到PEI后获得507个PEI相关基因（含DMR），再整合转录组得到43个差异表达基因（如App、Prkce、Nell2）。通路分析仅识别出一个基因本体论术语：神经元突触可塑性调节。进一步分类发现，34%的基因（14个）表现为下调、PEI减少和高甲基化，42%的基因（17个）表现为下调、PEI增多和高甲基化，提示染色质相互作用与DNA甲基化可能存在协同或竞争关系。
**讨论与结论**
研究表明，单一组学数据揭示了广泛的分子紊乱，但多组学整合将异质性降低至43个基因，并精炼至一个核心通路（突触可塑性调节），与人类ASD死后脑组织研究一致。讨论部分指出，尽管存在局限性（单一模型、脑区、性别及细胞类型异质性），但该研究证明了多组学整合在解析复杂ASD分子病因学中的价值。结论翻译：总之，研究人员在单基因敲除ASD小鼠模型海马组织中的多组学方法提高了复杂ASD分子病因学的分辨率，并与人类ASD死后脑组织研究一致。这些发现为理解与ASD异质性遗传架构和分子病因学相关的趋同生物学过程提供了见解。
