IF31.7｜系统性分析snRNA基因揭示RNU2-2变异在显性和隐性发育性癫痫脑病中高频出现

RNU2-2变异是神经发育障碍的主要贡献者，隐性形式至少是显性的两倍。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

小核RNA（snRNA）是剪接体的核心组分。2024年，RNU4-2（ReNU综合征）、RNU5B-1和RNU2-2的de novo变异被报道与显性神经发育障碍（NDD）相关，揭示了非编码RNA基因在遗传病中的意外贡献。然而，snRNA基因高度冗余且注释不完整，许多功能性基因被误标为假基因，导致变异解读困难。此外，snRNA基因具有高突变性，de novo变异率是其他区域的十倍，这为区分功能基因与假基因提供了线索，但也增加了显性与隐性遗传模式交织的复杂性。

科学问题

卡在哪里： 此前研究主要聚焦于显性de novo变异，对snRNA基因的隐性遗传贡献缺乏系统评估。RNU2-2虽已被关联到显性NDD，但其隐性致病变异谱、临床表型及与显性形式的关系尚不清楚。

本文要回答： 系统分析200个潜在功能性snRNA基因在大型罕见病队列中的变异，明确RNU2-2在显性和隐性NDD中的贡献及其临床和分子特征。

技术路线

作者主要用了：全基因组测序、RNA-seq、DNA甲基化芯片、结构预测、变异富集分析。

作者首先通过整合ENCODE cCRE、超突变性和HGNC注释，从2094个snRNA基因中筛选出200个潜在功能性基因。随后在法国PFMG队列（34,329例）中提取这些基因的de novo和双等位变异，并通过与已解决/未解决病例及NDD/非NDD病例的比较进行富集分析。对富集显著的RNU2-2，进一步扩大国际队列验证，结合临床表型、结构预测、淋巴细胞转录组和血液甲基化分析，解析变异的功能影响。

图 1　Fig. 1。图1展示了从2094个snRNA基因中筛选200个潜在功能性基因的策略，以及它们在人类脑组织中的表达情况。过滤后，主要剪接体snRNA基因拷贝数远多于次要剪接体，U6最多。ENCODE小RNA-seq数据显示，87个snRNA在脑中有表达，包括39个注释为假基因的基因，提示更多NDD可能由snRNA变异驱动。

核心亮点

亮点 1｜RNU2-2双等位变异是高频隐性NDD病因

在PFMG队列中，双等位RNU2-2变异在NDD病例中显著富集（P=8.9×10^-5）。结合国际协作，共发现73个家系91例隐性患者，其中54例复合杂合、17例纯合、2例半合子合并基因缺失。隐性形式至少是显性的两倍，且常涉及一个de novo变异与一个遗传变异反式组合，反映了snRNA基因的高突变性。

图 2　Fig. 2。图2比较了已解决与未解决病例及NDD与非NDD病例中稀有变异的比例。de novo变异在RNU4-2和RNU2-2中显著富集于未解决和NDD病例；双等位变异仅在RNU2-2中显著富集（NDD vs 非NDD，P=8.9×10^-5），表明RNU2-2同时贡献显性和隐性NDD。

亮点 2｜显性与隐性RNU2-2 NDD表型重叠

112例患者的临床分析显示，显性（n.4G>A、n.35A>G）和隐性患者均表现为全面发育迟缓、重度/极重度智力障碍（77%）和高癫痫发生率（85%），起病多在3岁前。隐性患者肌阵挛发作更常见（75% vs 29%），且所有8例死亡均发生在隐性组。家族内表型一致，提示变异组合影响严重程度。

图 3　Fig. 3。图3展示了U2-2 snRNA的二级结构及本研究发现的点变异位置和遗传模式。显性致病变异n.4G>A和n.35A>G位于5'端关键功能域，而隐性变异分布更广，包括BSL、SLIIa、Sm位点和3'茎环，提示不同区域变异的功能后果差异。

亮点 3｜变异位置决定功能影响梯度

结构预测表明，5'端变异（如n.4G>A、n.35A>G）破坏U2/U6 helix II或BSL，可能被整合进剪接体但功能异常；而Sm位点或3'茎环变异可能阻碍U2 snRNP成熟，行为类似无效等位基因。这支持“影响梯度”模型：显性变异单独足以致病，而其他变异需与第二个等位基因组合才能达到致病阈值。

图 4　Fig. 4。图4显示了16个家系中双等位变异的共分离情况，以及显性和隐性患者的典型面部特征。隐性患者面部畸形轻微且非特异，与显性患者相似，支持两者表型重叠。家系分析证实变异与疾病共分离，支持致病性。

亮点 4｜血液转录组和甲基化改变微弱且异质

对19例携带者的淋巴细胞RNA-seq分析未发现一致剪接特征，但线性回归显示显性变异主要导致外显子跳跃增加，隐性变异则双向改变。DNA甲基化分析在24例患者中仅检测到微弱的变异特异性信号，n.4G>A携带者呈低甲基化，n.35A>G呈高甲基化，隐性携带者异质性大，整体效应小于RNU4-2。

图 5　Fig. 5。图5概述了RNU2-2变异在U2 snRNA结构上的分布及其对剪接体组装的预测影响。5'端变异破坏U2/U6 helix II或BSL，可能影响剪接体组装和分支点识别；3'端变异可能阻碍snRNP成熟，导致功能丧失。这为理解显性与隐性遗传的分子基础提供了框架。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 全基因组测序、RNA-seq、DNA甲基化芯片、结构预测、变异富集分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事罕见病基因组分析、非编码变异解读和神经发育障碍研究的生信人员与临床遗传学家。

带走一句

分析snRNA基因变异时，必须同时考虑显性和隐性模式，de novo变异应积极寻找反式第二个等位基因；高突变性要求更严格的致病性评估。