尽管超过100个致病基因具有功能多样性，不同模型间的表型汇聚可能揭示孤独症谱系障碍（Autism Spectrum Disorder, ASD）共同的神经生物学过程。研究人员对11种ASD相关单基因敲除小鼠模型在三个发育阶段、两种性别及两个脑区的251份样本进行单细胞核多组学测序（single-nucleus multi-omic sequencing，同时检测snRNA-seq与snATAC-seq），发现尽管存在遗传异质性，ASD相关突变共同导致放射状胶质细胞（Radial Glial Cell, RG）谱系扰动——表现为短暂的发育延迟而非永久性谱系错误指定，并在出生后阶段恢复正常。分子水平上，最大的转录差异出现在出生后早期神经元中，包括突触及离子通道相关基因下调，符合稳态适应或成熟延迟特征。网络分析显示各发育阶段内模型间存在分子汇聚，提示不同ASD相关突变影响共同、具阶段特异性的发育过程；此汇聚现象在出生后第14天（Postnatal Day 14, P14）减弱，表明ASD相关改变具动态性。跨基因型异质性叠加于阶段特异性效应之上。电生理记录验证了该模式：突变体普遍显示神经元兴奋性及突触特性改变且具模型特异性细微差异。研究还发现性别特异性基因表达改变，雌性小鼠常较雄性表现出更大效应量。综上，本研究全面揭示了ASD模型发育过程中细胞与分子动态变化。
研究背景与意义
孤独症谱系障碍（Autism Spectrum Disorder, ASD）具有高度遗传异质性，已鉴定百余个单基因致病因素，但分散于多基因的突变如何导致相似表型、不同遗传性ASD模型间相似程度如何仍不清楚。既往单细胞研究多聚焦人类样本或单一模型，缺乏跨多模型、多发育阶段、多脑区及双性别的系统比较。该研究通过对11种ASD高危基因单基因敲除小鼠模型在胚胎期E14.5、出生后P4及P14，涵盖前后脑两区域及雌雄两性共251个样本的单细胞核多组学（snRNA-seq＋snATAC-seq）分析，系统描绘细胞状态与分子动态，旨在揭示ASD相关突变是否汇聚于共同发育程序及其时空与性别特征。论文发表于《Nature》。
主要关键技术方法
研究人员构建或获取11种C57BL/6J背景ASD高危基因杂合/纯合敲除小鼠品系（Ash1l+/−、Bckdk−/−、Cul3+/−、HnrnpU+/−、Kdm6b+/−、Kmt5b+/−、Setd5+/−、Trip12+/−、Usp7+/−、Wac+/−及X连锁Ptchd1−/y），取E14.5前脑、P4及P14前脑（含海马）与P4小脑，采用胆固醇修饰寡核苷酸（Cholesterol-Modified Oligonucleotide, CMO）多重标记结合10x Genomics单细胞核多组学（snMultiome-seq）建库测序；经质控、双细胞去除、CMO解复用后整合分析，进行细胞类型注释、拟时序（Pseudotime）分析、伪批量差异表达分析（Pseudobulk DEG with DESeq2）、基因集富集分析（Gene Set Enrichment Analysis, GSEA）、蛋白互作网络共定位（Protein–Protein Interaction network colocalization）、染色质可接近区差异分析（Differentially Accessible Region, DAR）、性别分层分析及全细胞膜片钳电生理记录（Whole-cell Patch-clamp），部分模型辅以EdU掺入增殖标记、RNAscope原位杂交及幼鼠步态分析。
研究结果
Single-nucleus multi-omics profiling（单细胞核多组学图谱绘制）
研究人员优化CMO多重标记snRNA-seq与snATAC-seq同步检测，获得E14.5、P4、P14前脑及P4小脑超20万细胞核，鉴定出主要细胞类型（放射状胶质细胞Radial Glial Cell, RG；中间祖细胞Intermediate Progenitor, IP；皮质板神经元Cortical Plate Neuron, CPN；上层/深层兴奋性神经元；PV+、SST+、VIP+抑制性神经元；少突胶质前体细胞Oligodendrocyte Precursor Cell, OPC；星形胶质细胞Astrocyte；小胶质细胞Microglia等），各样本与基因型均匀分布于聚类，验证数据集可靠性，为后续分析提供发育皮层细胞图谱。
RG progression defects in mutants linked to ASD（ASD相关突变体中RG谱系进程缺陷）
比较突变体与野生型（Wild Type, WT）细胞比例发现E14.5皮质板细胞与OPC减少、RG细胞比例升高，P4未成熟星形胶质细胞与OPC减少，提示RG→神经元/胶质命运转变延迟（gliogenic switch受影响）。拟时序分析显示突变体RG细胞滞留于增殖态RG1亚群（高表达Mki67），分化推迟；EdU标记实验证实E14.5 HnrnpU+/−与Kdm6b+/−胚胎RG增殖池扩大。RG1差异基因有限（仅Rsrp1与Galnt17下调），暗示延迟可能源于更早发育事件或转录后调控，不同突变体具各自RG差异表达谱（如HnrnpU影响最广）。结论：ASD相关突变共享RG谱系发育延迟，但潜在机制汇聚有限。
Convergent early postnatal changes in mutants linked with ASD（ASD相关突变体早期出生后汇聚性改变）
对各细胞类型做突变体合并vs WT伪批量差异分析发现P4上/下层兴奋性神经元差异基因（Differentially Expressed Gene, DEG）最多（715与600个），且多为下调，其他时期与细胞类型较少；单独各突变体分析显示P4各基因型与多细胞类型广泛转录改变（尤以兴奋性神经元显著），E14.5除HnrnpU与Kdm6b外较少，P14 DEG数下降。下采样平衡测序深度后趋势不变。多突变共享约10% P4兴奋性神经元DEG，随机模拟验证非偶然汇聚；频繁DEG富集ASD风险基因及突触蛋白编码基因（34%为ASD关联基因，q＝1.03×10−15），方向一致。结论：ASD突变汇聚于P4神经元突触/离子通道相关基因下调这一阶段性共同响应，叠加各基因型独有改变。
Functional pathway disruptions in mice with mutations linked with ASD（ASD相关突变体功能通路扰动）
GSEA显示E14.5部分突变体富集细胞周期/染色体分离/神经元突起发育条目；P4突出突触与离子通道功能相关基因集扰动，P14仍见突触条目。抑制性神经元受影响弱于兴奋性神经元但定性相似；多数胶质细胞影响较小（Cul3与Ptchd1除外）。各突变体具独特富集条目（如Bckdk富集氨基酸代谢与三羧酸循环相关基因，Trip12富集糖酵解与脂代谢，Kdm6b富集WNT信号）。结论：突触传递与离子通道通路失调是跨模型共性，特定基因型另有 cytoskeleton dynamics、纤毛功能或代谢相关特异性改变。
Network analysis reveals convergent phenotypes（网络分析揭示汇聚表型）
基于蛋白互作网络将DEG聚类为功能模块并评估各突变体对各模块表达改变效应量，发现E14.5富集细胞骨架/细胞周期（HnrnpU、Kdm6b、Cul3显著），P4/P14转向突触过程/离子转运/代谢。同阶段不同突变体DEG网络定位相似性高于随机期望，表明影响共同网络区域；P14相似性降低，尤其兴奋性神经元中。结论：ASD相关突变作用于共同阶段特异性生物过程网络，随发育渐分化为基因型特异性转录程序。
Synaptic and ion channel changes in ASD models（ASD模型中突触与离子通道改变）
重点关注钙/钾/钠通道（如Scn1a下调伴Scn8a下调与Scn9a上调可能为代偿）及突触相关基因（Rims1/2、Stxbp5l下调，Cacng4上调，Camk2a下调等），多属SynGO注释的突触囊泡释放、突触后组织及离子通道类别。电生理选取Bckdk−/−、HnrnpU+/−、Trip12+/−、Usp7+/−于P7层II/III锥体细胞记录：多数突变体固有兴奋性（频率–电流关系）与迷你兴奋性/抑制性突触后电流（miniature Excitatory/Inhibitory Postsynaptic Current, mEPSC/mIPSC）频率降低（匹配突触前释放机器基因下调），离子通道表达偏离度与电生理偏离度Spearman相关；HnrnpU因共存上调补偿基因（如Ap3b2、Tprgl1、Akap12）使miniature电流表型轻微。最小电导模型按转录本倍数调整可复现电生理曲线。结论：离子通道与突触基因改变是ASD突变跨模型分子主题，转录–生理耦合稳健但受转录后缓冲调节。
Chromatin accessibility reveals regulatory dynamics（染色质可接近性揭示调控动态）
snATAC-seq显示E14.5差异可接近区（Differentially Accessible Region, DAR）最多，RG细胞DAR数最高，P4兴奋性神经元仍较多；合并突变体偏向染色质可及性降低，匹配转录下调。HnrnpU与Kdm6b RG细胞DAR涉及神经发生/前体分化相关TF基序（Sp2、Klf4、Zic2等）及WNT/BMP信号通路基因本体（Gene Ontology, GO）条目。snRNA-seq与snATAC-seq log2FC具相关性（RG与CP细胞较强）。结论：早期发育RG细胞染色质景观改变可能是转录延迟上游原因，部分突变直接体现为染色质水平调控异常。
Sex-specific gene expression alterations（性别特异性基因表达改变）
WT雌雄对照差异极小。性别分层分析示P4多个突变体（Ash1l、Bckdk、HnrnpU、Trip12、Usp7、Wac）雌性DEG多于雄性；Trip12主要影响雌性，Cul3两性相近。定义female-biased DEG（fbDEG）与female-specific DEG（fsDEG），fbDEG在雌性中效应量（Cohen's d）更大，fsDEG跨模型富集突触、离子通道与细胞黏附基因（如ApoE在雌性抑制性神经元中Ash1l/Bckdk/Trip12/Usp7突变体特地下调）。Wac突变体fsDEG占比最高，契合女性偏易感性报道。结论：部分ASD基因敲除在雌性小鼠引发更强转录效应，性别是ASD模型分子表型重要修饰因子。
Restricted cerebellar effects in mutants linked with ASD（ASD相关突变体局限于小脑的改变）
P4小脑仅Bckdk突变体广泛DEG（涉及氨基酸代谢与突触，与前脑重叠），Cul3见颗粒细胞/高尔基细胞比例微变；其余突变体无明显共享改变。Bckdk幼鼠（P13–15）表现步长缩短、步速降低、前爪协调减弱等运动缺陷，Kmt5b无此表型。结论：多数ASD突变在P4不广泛重塑小脑转录组，Bckdk例外且伴运动表型，与其代谢通路特异扰动一致。
讨论（Discussion）总结翻译
通过对多ASD模型跨发育阶段检测，本研究强调ASD相关改变的动态本质，勾勒出连贯轨迹：多样ASD相关突变冲击共同阶段依赖性程序（如RG发育延迟、出生后早期离子通道/突触模块下调），同时保留随成熟加剧的基因型特征。网络定位证实汇聚在发育早期最强而P14减弱。功能上转录组–电生理耦联与最小电导建模提示离子通道表达微小偏移足以产生多样表型。结果支持针对共同早期发育改变的阶段与性别感知干预策略，并提出靶向早期共有变化能否重定向后期疾病轨迹的开放问题。
