多层面转录组学整合揭示抽动秽语综合征中以γ-氨基丁酸（GABA）为核心的分子状态及外周—纹状体网络关联

《Clinical and Translational Discovery》：Integrated multi-layer transcriptomics reveals a gamma-aminobutyric acid-centred molecular state in Tourette syndrome

目的：γ-氨基丁酸（GABA）能功能障碍是抽动秽语综合征（Tourette syndrome，TS）的核心，但连接外周信号、上游调控、药理反应和纹状体细胞背景的证据仍较碎片化。研究人员探讨TS中GABA定义的外周分子状态是否可与候选调控信号、药理反应、纹状体细胞类型定位及局部通讯网络相联系。方法：研究人员整合外周血信使RNA（mRNA，GSE30470）、血清微小RNA（miRNA）谱（GSE146509）、慢性地西泮脑转录组（GSE76700）和纹状体单细胞RNA测序（GSE116470）。TS样本按GABA_index分层，GABA_index由GABA相关基因表达推导。分析包括差异表达、miRNA靶基因富集、疾病-药物逆转、细胞类型定位和配体-受体通讯。结果：GABA_index分层将TS分为GABA_Low和GABA_High状态，揭示1290个差异表达基因。唯一高置信度血清miRNA为hsa-miR-483-3p；其验证靶基因在神经营养因子、Ras、Rap1、轴突导向和cAMP通路富集，提示上游皮质-纹状体-丘脑-皮质（CSTC）调控。慢性地西泮在皮质逆转TS转录特征，负向富集TS_Up并正向富集TS_Down基因。单细胞分析将GABA/CSTC程序定位至中型棘状神经元和中间神经元；通讯分析揭示以细胞外基质/黏附为中心的纹状体微环境，并有成纤维细胞和血管软脑膜细胞贡献。结论：TS与以外周GABA为核心的分子状态相关，该状态与上游miRNA调控及纹状体神经元微环境相连。

抽动秽语综合征（Tourette syndrome，TS）是一种儿童期起病的神经发育障碍，以多个运动抽动和至少一个持续一年以上的发声抽动为特征，并常伴注意缺陷多动障碍、强迫症、焦虑等行为问题。其表型异质性明显，现有分子分层标记不足。当前病理生理模型将TS置于功能失调的皮质-纹状体-丘脑-皮质（CSTC）环路中，纹状体是异常抑制控制的核心部位；多巴胺、谷氨酸、5-羟色胺、组胺和γ-氨基丁酸（GABA）等均被涉及，其中GABA因纹状体以GABA能投射神经元和中间神经元为组织核心而居中心地位。既往人神经生理、磁共振、尸检转录组、基底节类器官及小鼠中间神经元耗竭研究支持GABA/中间神经元相关异常，但外周信号、上游调控、药理反应与纹状体细胞背景之间仍缺乏整合。外周血转录组可反映疾病相关变异并与抽动严重度相关，但外周信号不能直接等同于脑转录变化；循环miRNA稳定且可提供转录后调控信息，然既往TS血清/细胞外囊泡miRNA信号有限且不一致。药物治疗仍以对症为主，慢性地西泮暴露可作为增强GABA_A（γ-氨基丁酸A型受体）信号的药理扰动模型，但疾病-药物匹配在TS中较少用于评价GABA相关纠正信号。单细胞RNA测序可提供GABA/CSTC程序所在细胞类型及局部互作背景。因此，研究人员开展多层整合分析，以确定外周GABA定义的分子状态能否连接调控、药理、细胞及微环境特征；其意义在于构建跨外周状态、候选调控、药理敏感性、纹状体细胞定位和局部通讯架构的GABA中心网络，支持该网络在TS中的生物学相关性。

关键技术方法概括（不超过250字）：研究人员基于四个公共数据集开展多层分析：GSE30470为人TS外周血mRNA；GSE146509为人血清miRNA队列，含正常对照、TS、Arnold-Chiari畸形（AC）及AC+TS；GSE76700为慢性地西泮处理与载体对照成年雄性C57BL/6小鼠皮质、海马和杏仁核mRNA；GSE116470为成年雄性C57BL/6小鼠纹状体单细胞RNA测序。研究人员构建GABA_index并按中位数将TS分为GABA_High与GABA_Low，采用limma差异表达、miRTarBase验证靶基因富集、fgsea疾病-药物逆转、Seurat单细胞聚类与模块评分，以及配体-受体通讯分析。

研究结果：3.1 GABA_index分层界定TS的两种外周转录组状态。通过中位数GABA_index将21例TS分为11例GABA_Low和10例GABA_High；差异表达分析识别1290个名义p<0.05基因，其中345个在GABA_High较高、737个在GABA_Low较高，208个未达预设倍数变化阈值。通路分析显示神经活性配体-受体互作和神经活性配体信号向GABA_High富集，核糖体和帕金森病相关通路向GABA_Low富集；GABA_Low上调基因富集线粒体组织、线粒体DNA代谢、线粒体翻译、RNA/tRNA修饰、维生素代谢和电子传递，GABA_High上调基因富集血管生成/血管发育、磷脂酰肌醇3-激酶/蛋白激酶B（PI3K/AKT）信号、轴突发育、神经元迁移和神经元细胞黏附。该分层建立了后续分析锚点。3.2 纹状体细胞类型定位。单细胞分析解析20个转录簇，包括D1/D2中型棘状神经元（MSN）、多种中间神经元、星形胶质、少突胶质、血管、室管膜、成纤维细胞和血管软脑膜细胞（VLMC）；神经元再聚类区分D1-MSN、D2-MSN、D2-MSN样、MSN_Tshz1相关群、活动相关神经元状态及PV、SST和胆碱能中间神经元。GABA相关基因主要集中于神经元，MSN和中间神经元信号最强；模块评分显示GABA代谢、GABA_A受体、GABA突触和CSTC相关特征在D1_MSN、D2_MSN及MSN_Tshz1相关群优先富集，内皮、成纤维-VLMC和多数胶质较低。结论为GABA/CSTC程序主要定位于D1/D2 MSN和中间神经元。3.3 代表转录本紧密追踪GABA_index定义状态。通过差异表达筛选并结合分组比较和Spearman相关，LOC285628、PRADC1在GABA_High较低且与GABA_index负相关；LOC101927112、MEG8、TMPRSS12在GABA_High较高且与GABA_index正相关。上述五个转录本提供外周GABA_index状态的基因级特征。3.4 miR-483-3p靶基因汇聚于上游CSTC调控通路。血清miRNA分析中hsa-miR-483-3p为唯一符合跨对比标准者；其172个实验验证靶基因与Stage 1基因集直接重叠有限，C10orf10与GABA_High较高表达基因重叠，BTBD3、TSPAN3、SORD、CYC1、CKB与GABA_Low较高表达基因重叠，Fisher检验提示富集有限。靶基因功能富集显示形态发生、树突棘发育、突触可塑性、突触组织、神经元投射组织和黑质发育等过程，KEGG富集Ras、Rap1、神经营养因子、轴突导向和cAMP信号。结论为hsa-miR-483-3p构成连接外周GABA状态与CSTC相关神经发育和突触程序的候选调控环节。3.5 慢性地西泮优先在皮质逆转TS相关转录特征。通过比较地西泮签名与TS基因集，基因水平逆转比在皮质为0.518、海马为0.519，均高于0.5（单侧精确二项检验p=5.5×10^?6和1.2×10^?6），杏仁核为0.478且p=1.0，无方向逆转证据。皮质显示最清晰相反模式，TS_Up负富集、TS_Down正富集（标准化富集分数NES=-1.459和1.395）；海马同向但较弱且无显著富集，杏仁核无一致逆转。皮质基因集富集分析（GSEA）确认TS_Up负富集（错误发现率FDR=0.0063）、TS_Down正富集（FDR=0.001）。结论为慢性地西泮提供GABA中心框架的药理支持，皮质最清晰。3.6 通路解析的细胞间通讯突出以细胞外基质/黏附为中心的纹状体微环境。通讯分析显示PTN、NRXN、NCAM、LAMININ和MK信号在神经元与非神经元群体间广泛分布，PTN和LAMININ交互广泛且强度较高，NRXN和NCAM支持密集的神经元中心通讯；EPHA、EPHB和OPIOID信号较弱且更具选择性，COLLAGEN信号集中于成纤维-VLMC及相关非神经元结构轴。结论为MSN/中间神经元GABA/CSTC程序位于以细胞外基质（ECM）/黏附为中心的网络中，成纤维-VLMC和血管相关隔室参与局部细胞背景。3.7 TS中GABA中心分子状态的整合概览。跨层整合识别由外周GABA_index及代表转录本锚定的GABA中心网络，GABA/CSTC程序主要定位于纹状体MSN和中间神经元；调控层为miR-483-3p及神经营养因子、Ras、Rap1、轴突导向和cAMP相关靶网络；药理层显示皮质地西泮签名最清晰相反；通讯层将神经元程序置于涉及PTN、NRXN、NCAM、LAMININ和COLLAGEN信号及成纤维-VLMC、血管相关隔室贡献的ECM/黏附中心纹状体网络。结论为这些跨数据集关系定义了一个涵盖外周状态、候选调控、药理敏感性、纹状体细胞定位和局部通讯架构的GABA中心分子网络。

讨论部分总结：该研究中心发现是以GABA为中心的分子网络连接外周GABA_index分层与候选调控、药理敏感性、纹状体细胞定位和局部通讯架构。既往血GABA、乙酰胆碱和儿茶酚胺相关基因表达与抽动严重度相关，本研究分层进一步显示GABA聚焦签名可将TS分为具有协调通路改变的不同状态；该解释与TS儿童和成人抑制性张力或区域GABA异常、尸检和类器官中间神经元病理及纹状体转录紊乱证据一致。miR-483-3p作为唯一符合预设标准的血清miRNA被视为候选上游调控因子；循环miRNA为上游转录后调控因子，跨组织比较不预期直接对应，现有TS外周生物标志物证据异质，但其靶基因富集于神经营养因子、Ras、Rap1、轴突导向和cAMP通路，与突触发育和活动依赖可塑性相关，miR-483家族亦与神经调节相关。慢性地西泮在皮质最清晰逆转TS相关签名，与TS皮质抑制受损及皮质GABA异常报道一致；单细胞分析提供外周数据不能提供的解剖背景，GABA/CSTC签名集中于MSN和中间神经元，符合纹状体细胞组织及单细胞图谱，并得到PV阳性/胆碱能中间神经元减少、纹状体转录组异常、基底节类器官发育错配和中间神经元缺陷、实验性中间神经元耗竭产生抽动样或重复行为、局灶阻断纹状体GABA_A受体诱导抽动样运动及正向调节GABA_A受体信号减少抽动相关行为的支持。通讯分析将神经元程序置于ECM/黏附中心微环境和局部互作网络，成纤维-VLMC、层粘连蛋白、NCAM、neurexin、胶原和pleiotrophin等信号突出，CNS血管周围成纤维细胞支持血管稳定、基质组织和细胞间信号，黏附与基质分子影响突触稳定和可塑性；因此TS中GABA功能障碍的解释扩展至环路-微环境互作、结构支持和局部信号架构。局限性包括跨不同队列、组织和物种整合，血清miRNA信号有限，地西泮作为机制性GABA_A探针且临床治疗疗效不在分析范围，miR-483-3p和代表GABA_index相关转录本的基因特异性因果作用需未来在TS相关神经元或动物模型中进行扰动研究；功能支持来自跨数据集一致性、实验验证miR-483-3p靶基因、地西泮相关转录逆转、纹状体细胞类型定位、配体-受体网络映射及已发表纹状体GABA能信号和中间神经元功能扰动研究。

研究结论翻译：本多层面分析识别TS中以GABA为中心的分子网络，该网络锚定于外周GABA_index定义的状态。调控层识别miR-483-3p及神经营养因子、Ras、Rap1、轴突导向和cAMP信号作为与CSTC相关神经发育和突触程序的候选联系。慢性地西泮在皮质产生最清晰的相反转录模式，为GABA_A受体介导的疾病相关状态调节提供药理支持。单细胞定位将GABA/CSTC相关程序主要定位至纹状体中型棘状神经元和中间神经元，细胞间通讯分析将这些群体与PTN、NRXN、NCAM、LAMININ和COLLAGEN相关信号相连。这些分子、药理、细胞和通讯层的跨数据集一致性支持所提出GABA中心网络在TS中的生物学相关性。