《Frontiers in Cell and Developmental Biology》：Genotype-specific molecular and functional trajectories in NGLY1 deficiency: a comparative multi-omics study of patient-derived iPSC neurons
编辑推荐：
背景： NGLY1缺乏症是一种超罕见的先天性去糖基化疾病，以显著的神经发育障碍和明显的等位基因异质性为特征。此前的神经元模型大多在单一遗传背景下表征NGLY1缺失，尚不清楚不同患者基因型是汇聚于共同的分子状态还是遵循不同的路径。由于不同基因型的临床表现相互重叠
背景： NGLY1缺乏症是一种超罕见的先天性去糖基化疾病，以显著的神经发育障碍和明显的等位基因异质性为特征。此前的神经元模型大多在单一遗传背景下表征NGLY1缺失，尚不清楚不同患者基因型是汇聚于共同的分子状态还是遵循不同的路径。由于不同基因型的临床表现相互重叠，解答这一问题需要对多种基因型进行系统内直接比较，而非仅对单一细胞系进行表征。
方法： 研究人员从三名携带不同双等位NGLY1基因型（R401X/R401X、Q208X/G310G、L318P/R390P）的患者及两名健康对照中生成iPSC来源的谷氨酸能神经元，并在第28天使用匹配的转录组学与蛋白质组学进行头对头分析，同时开展RNA-蛋白质一致性分析、差异外显子使用分析、通路与表型网络整合，以及早期网络成熟过程中的纵向微电极阵列（MEA）记录。
结果： 在两个分子层面，样本均按基因型强烈分离，而非形成统一疾病特征。差异表达分析显示基因型间存在显著异质性，仅有一个较小的共享分子核心。RNA-蛋白质一致性为中等水平且具有基因型依赖性，差异外显子使用在不一致基因中呈突变依赖性富集，提示剪接水平调控参与转录-蛋白质差异的形成。整合网络将基因型特异性一致变化与突触信号传导、钙处理、代谢及细胞外基质重塑相关联，并通过不同分子路径连接到重叠的临床表型。MEA记录显示网络成熟延迟及基因型依赖性的网络兴奋性变化，其功能聚类与分子分离相呼应。
结论： 不同的NGLY1基因型产生不同的分子与功能轨迹，通过非相同的路径汇聚于网络成熟受损。这些发现将NGLY1缺乏症重新定义为一系列基因型特异性分子状态，而非单一表型，对机制亚组划分及临床前模型和验证研究的设计具有直接意义。
研究解读：NGLY1缺乏症的基因型特异性分子与功能轨迹——一项患者来源iPSC神经元的多组学比较研究
一、研究背景与科学问题
NGLY1缺乏症是一种超罕见的常染色体隐性遗传性去糖基化疾病，由NGLY1基因双等位功能缺失变异引起。NGLY1编码胞质N-聚糖酶1，参与内质网相关降解（ERAD）途径，负责将错误折叠糖蛋白从内质网逆转运至胞质后去除其N-连接聚糖。临床上该病表现为多系统神经发育综合征，包括发育迟缓、运动异常、无泪症、癫痫、周围神经病变及程度不一的肝脏受累。在分子层面，NGLY1缺失破坏蛋白质稳态，并损害NFE2L1的激活——该转录因子是蛋白酶体在蛋白毒性应激后恢复所必需的。
该疾病表现出明显的等位基因和临床异质性，基因型-表型关系在患者队列中尚未完全阐明。既往人类神经元模型大多在单一遗传背景下表征NGLY1缺失的后果，例如R401X突变患者来源神经干细胞与CRISPR校正的同源对照比较显示内质网应激、氧化应激和自噬通量升高；运动神经元和神经肌肉接头模型则显示电生理改变和突触功能受损。然而，这些基于单一背景的研究留下了一个关键问题：鉴于NGLY1缺乏症显著的等位基因异质性和患者间重叠的临床表现，不同基因型是汇聚于共享的神经元分子状态，还是通过不同路径走向共同的临床谱系？回答这一问题需要对多种基因型进行系统内直接比较。
二、研究设计与主要发现
研究人员首次在单一神经元系统内对三种不同NGLY1基因型进行了头对头比较。研究使用三名患者来源的iPSC系——分别携带纯合无义突变R401X/R401X（592C，预测为无功能蛋白）、复合杂合突变Q208X/G310G（519D，预测为部分功能缺失）和复合错义突变L318P/R390P（594D，预测为蛋白不稳定或活性降低）——以及两名健康对照（WC31为主要对照，BD3为补充对照），在体外第28天进行匹配的转录组学和蛋白质组学分析，并结合RNA-蛋白质一致性分析、差异外显子使用分析、通路-表型网络整合，以及第11至29天的纵向微电极阵列（MEA）记录。
核心发现是：不同NGLY1基因型产生不同的分子和功能轨迹，仅共享有限的分子核心，却汇聚于网络成熟受损这一共同功能终点。转录组和蛋白质组的PCA分析均显示三个疾病系与对照清晰分离，且各疾病系占据不同的PCA空间区域。差异表达呈梯度响应：592C受影响最严重（RNA层面9,040个差异表达基因，蛋白质层面1,081个差异蛋白），594D居中（5,573个DEG，236个差异蛋白），519D最轻（2,509个DEG，113个差异蛋白）。这些变化在两种数据集中均以细胞系特异性为主——592C独有5,112个DEG和927个差异蛋白，而三条系共享的仅404个DEG和7个蛋白。
三、关键技术方法
研究采用的主要技术方法包括：1）患者来源iPSC向谷氨酸能神经元的分化培养体系，样本来源于三名携带不同双等位NGLY1基因型的患者和两名健康对照；2）RNA测序（RNA-seq）与串联质量标签（TMT）标记的液相色谱-串联质谱（LC-MS/MS）蛋白质组学匹配分析；3）edgeR和limma差异表达分析及RNA-蛋白质一致性整合分析（含Deming回归、Bootstrap重采样和模拟分析）；4）DEXSeq差异外显子使用分析；5）clusterProfiler通路富集分析及突变-通路-症状网络整合；6）Axion Maestro系统多电极阵列（MEA）纵向电生理记录。
四、研究结果
神经元表型背景： 免疫荧光分析显示所有细胞系均能形成MAP2阳性的树突过程和相互连接的神经元网络，但NGLY1疾病系表现出更明显的细胞聚集和MAP2阳性过程分布不均，提示突变依赖性的神经元形态和网络组织改变。
多组学表达谱区分疾病细胞系与对照： PCA和层次聚类在转录组和蛋白质组层面均将三个疾病系与对照清晰分离，且各疾病系占据不同空间区域。差异表达呈梯度响应，592C受影响最重，变化以细胞系特异性为主。
转录失调在蛋白质层面显著衰减： 在4,128个共同检测的基因产物中，大多数属于"两者均无变化"类别。RNA-蛋白质Spearman相关性从519D的0.30到592C的0.57不等，Deming回归斜率均低于1（0.21-0.44），表明蛋白质倍数变化相对于转录本持续压缩。方向一致性在双检测显著的特征中超过97%。模拟分析表明蛋白质层面发现不足不能仅由统计功效解释，提示存在真实的转录后衰减。
差异外显子使用标记不一致基因： 强不一致基因在519D中显著富集差异外显子使用（OR 2.05），但在592C和594D中不富集。NGLY1本身在519D和592C中显示显著基因水平DEU，XBP1仅在592C中显示局灶性显著变化。DEU负担与整体转录组响应一致，592C携带最多显著基因（7,332个）和外显子bin（49,752个）。
通路分析揭示趋同的细胞外基质重塑与细胞系特异性程序： 细胞外基质、胶原和粘附程序在所有三条系和所有分析层面均显著富集，是唯一在519D中达到一致显著性的程序。592C独特富集蛋白酶体、ERAD和泛素质量控制及线粒体能量代谢，同时强烈上调翻译和核糖体程序并广泛下调突触和神经元程序。594D以RNA驱动的突触传递、神经递质受体和GPCR信号下调为主。519D响应最受限，集中于细胞外基质组织。
整合一致标记与通路连接临床表型： 网络整合显示跨三条系一致改变的标记包括CXCL14、PTN、S100A6、GPNMB、SORCS2、SFRP2和MTTP等分泌和基质相关因子。这些标记通过基质相关通路连接到NGLY1缺乏症的多系统症状模块，包括生长障碍、神经发育损害、肌张力低下、运动障碍、癫痫等。
功能分析揭示网络活动改变： MEA记录显示对照培养表现出进行性的网络成熟，而NGLY1缺陷神经元显示降低的平均放电率、活性电极增加缓慢、爆发起始延迟和网络爆发参数改变。594D显示网络爆发数量增加但每次网络爆发尖峰数减少，提示混合的兴奋性表型。DIV29的49个MEA参数的PCA显示每个患者系形成不重叠的聚类，与分子数据的分离模式相呼应。
五、讨论与结论
讨论部分指出，该研究将NGLY1缺乏症重新定义为一系列基因型特异性分子状态而非单一表型。共享分子核心有限，而基因型特异性特征显著扩展。RNA-蛋白质不一致部分归因于转录后调控和蛋白质稳定性差异，部分归因于技术覆盖限制，但模拟分析证实存在真实的转录后衰减。差异外显子使用提供了不一致的机制解释——基因水平RNA-seq量化聚合了异构体读数，向无功能或NMD靶向异构体的转变可在总转录信号不成比例变化的情况下减少蛋白质输出。
研究结论指出，通过头对头比较三种不同NGLY1基因型的患者来源谷氨酸能神经元，研究人员证明NGLY1缺乏症不产生统一的神经元表型，而是产生仅共享有限共同核心却汇聚于网络成熟受损的不同轨迹。这些发现为NGLY1缺乏症的机制亚组划分和患者来源模型及验证靶点的合理选择建立了框架。由于该比较设计中基因型与遗传背景混杂，同源基因校正对照是确定哪些优先特征直接归因于NGLY1变异的必要下一步。
