IF31.7｜R-loop区变异分析揭示RNU2-2与RNU5B-1为神经发育障碍新基因

从R-loop区域富集的非编码变异中，锁定两个剪接体snRNA基因，为大量未确诊NDD提供遗传解释。

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

研究背景

R-loop是转录过程中形成的DNA-RNA杂交结构，其单链DNA暴露可能增加突变风险，但R-loop与人类孟德尔遗传病的关联尚未被系统探索。已有研究显示R-loop区域突变率升高，且某些非编码RNA基因（如RNU4-2）的变异可导致神经发育障碍（NDD）。然而，R-loop区域变异在罕见病队列中的整体分布模式，以及是否还有其他非编码RNA基因参与NDD发病，仍不清楚。

科学问题

卡在哪里： 此前研究未将R-loop区域变异与疾病队列进行系统性关联分析，也未充分利用R-loop特征来发现新的非编码致病基因。特别是剪接体snRNA基因家族中，除RNU4-2外，其他成员是否也携带致病性变异并导致NDD，缺乏证据。

本文要回答： 本文旨在分析R-loop区域变异在罕见病队列中的富集模式，并鉴定新的NDD相关非编码RNA基因。

技术路线

作者主要用了：全基因组测序、de novo变异分析、R-loop区域注释、小RNA-seq、表型关联分析。

作者首先利用100KGP罕见病队列的de novo变异数据，与已发表的R-loop共识区域取交集，比较R-loop区与非R-loop区的突变率。随后按基因生物型统计de novo R-loop区域变异（RRV）的富集情况，发现snRNA等非编码RNA基因显著富集。进一步聚焦于无已知疾病关联且存在多个新发RRV的基因，锁定RNU2-2和RNU5B-1。通过扩展多个国际罕见病队列验证，并利用小RNA-seq确认其表达，最终结合表型分析确定两种新的NDD综合征。

图 1　Fig. 1。图1展示了R-loop区域突变率升高及RRV在snRNA基因中的富集。b图小提琴图显示在启动子、外显子、内含子和随机区域中，R-loop重叠区域的de novo变异率均显著高于非重叠区域。c图点图显示核酶、snoRNA和snRNA基因在疾病队列中显著富集de novo RRV。e和f图分别展示RNU2-2和RNU5B-1的变异位置及gnomADv4中的变异耗竭区域，表明变异集中在高度受限的功能区。

核心亮点

亮点 1｜R-loop区域de novo变异显著富集

在100KGP罕见病队列中，5.4%的de novo变异位于R-loop区域，且在启动子、外显子、内含子及随机区域中，R-loop区的de novo变异率均显著高于非R-loop区（P<2.22×10^-16）。在独立冰岛对照队列中，除启动子外，其余区域也得到验证。进一步分析发现，实验验证的R-loop形成序列的de novo变异率显著高于仅生物信息学预测的序列，而两组序列的平均可突变性无显著差异，提示R-loop结构本身促进突变。

图 2　Fig. 2。图2展示了RNU2-2和RNU5B-1变异携带者的表型富集及临床特征。a和b图显示RNU2-2变异个体显著富集全面发育迟缓、肌张力低下和小头畸形等HPO术语；c和d图显示RNU5B-1变异个体显著富集肌张力低下、大头畸形和生长迟缓等。面部照片和MRI图像显示RNU2-2患者有小头畸形和脑体积减少，RNU5B-1患者有大头畸形和胼胝体发育不全等。

亮点 2｜snRNA基因在疾病队列中显著富集RRV

按基因生物型比较疾病队列与对照队列的de novo RRV富集情况，发现核酶、snoRNA和snRNA基因仅在疾病队列中显著富集（>2 log2 fold change）。其中snRNA基因包含18个基因的86个变异，包括已知的隐性致病基因和显性ReNU综合征相关RNU4-2（36个变异）。这一结果验证了该方法的有效性，并提示其他snRNA基因可能携带致病变异。

图 3　Fig. 3。图3展示了RNU2-2和RNU5B-1的表达及变异对RNA结构的影响。a图气球图显示在ENCODE发育中大脑各区域，RNU2-2表达高于RNU2-1，RNU5B-1也高表达。b图箱线图显示在人类脉络膜、神经感觉视网膜和RPE中，RNU2和RNU5旁系同源基因的表达。c和d图分别展示RNU2-2变异位于U2-U6相互作用区域和分支位点识别序列，RNU5B-1变异位于U5 loop I保守区，提示可能破坏剪接体功能。

亮点 3｜RNU2-2与RNU5B-1为新NDD基因

在100KGP中，RNU2-2的5个de novo变异（n.3C>A, n.4G>A, n.7_8insA, n.35A>C, n.35A>G）位于1-60bp的受限区域，见于9例患者；RNU5B-1的3个de novo变异（n.37G>C, n.42_43insA, n.44A>G）位于30-50bp的受限区域，见于4例患者。通过扩展搜索，共收集到27例RNU2-2变异和9例RNU5B-1变异个体。所有21例100KGP先证者此前均未确诊。

亮点 4｜RNU2-2相关障碍表型特征

RNU2-2变异个体显著富集严重全面发育迟缓（OR=33.6, P<10^-8）、全身性肌张力低下（OR=15.8, P=1.2×10^-8）和小头畸形（OR=9.0, P=1.1×10^-5）。详细临床信息显示，该综合征特征为全面发育迟缓、显著言语障碍、癫痫，部分表现类似Pitt-Hopkins或Rett综合征。

亮点 5｜RNU5B-1相关障碍表型特征

RNU5B-1变异个体显著富集全身性肌张力低下（OR=37.3, P=7.4×10^-5）、大头畸形（OR=30.6, P=1.8×10^-4）、生长迟缓（OR=24.7, P=4.5×10^-4）、全面发育迟缓（OR=22.1, P=5.6×10^-3）和眼部异常（OR=18.7, P=1.4×10^-3）。详细临床信息显示，该综合征特征为全面发育迟缓、相对大头畸形、癫痫发作和生长迟缓。

亮点 6｜RNU2-2与RNU5B-1在发育中大脑高表达

利用ENCODE人类发育中大脑的小RNA-seq数据，采用严格的多重比对处理策略，发现RNU2-2的表达高于其多拷贝旁系同源基因RNU2-1，且RNU5B-1也高表达。在人类视网膜、脉络膜和RPE组织中同样观察到高表达。值得注意的是，RNU2-2此前被注释为假基因，但本研究证实其编码的U2 snRNA可被整合进剪接体，且高表达，因此并非假基因。

生信可带走

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

• 方法栈： 全基因组测序、de novo变异分析、R-loop区域注释、小RNA-seq、表型关联分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事罕见病基因组学、非编码变异解读及神经发育障碍研究的科研人员。

带走一句

将R-loop区域信息整合进变异分析，可有效富集非编码致病基因，尤其适用于剪接体snRNA等非经典基因的发现。