《Frontiers in Bioengineering and Biotechnology》：HNRNPH2-R114W remodels RNA splicing and neurodevelopmental gene programs in isogenic cerebral organoids
编辑推荐：
自闭症谱系障碍（autism spectrum disorder, ASD）与皮质发生过程中的RNA加工异常有关，但特定变异相关的机制仍难以阐明。本研究利用携带HNRNPH2-R114W变异的患者来源诱导多能干细胞（induced pluripotent st
自闭症谱系障碍（autism spectrum disorder, ASD）与皮质发生过程中的RNA加工异常有关，但特定变异相关的机制仍难以阐明。本研究利用携带HNRNPH2-R114W变异的患者来源诱导多能干细胞（induced pluripotent stem cell, iPSC）系及其CRISPR校正后的同基因对应体所衍生的同基因脑类器官（isogenic cerebral organoids），作为人源发育模型，用以探究变异相关的RNA失调。在单一发育时间点对每种基因型的混合类器官RNA进行成像定量和探索性批量RNA测序（bulk RNA sequencing），研究人员发现神经元和代谢程序的协同重塑，以及涉及6,310个差异使用连接（junction）的连接中心剪接失调。计算结构对接将HNRNPH2的RNA识别基序（RNA recognition motifs, RRMs）1–3定位至剪接体背景下，并将残基R114置于预测的RNA接触界面附近，表明R114W替换可能改变RNA结合。将差异表达分析与事件水平和异构体水平分析整合后，研究人员识别出富含轴突导向和突触信号通路的调控枢纽。基于单细胞参考的去卷积分析提示存在适度的组成变化，但这些变化不太可能解释所观察到的广泛转录和剪接改变。综上，这些发现表明HNRNPH2-R114W参与广泛的RNA重连，并在人源皮质类器官中定义了可用于机制和治疗研究的候选异构体水平终点。这些分子特征为研究HNRNPH2相关RNA失调如何在发育时间进程中演变以及在独立生物学背景下重复出现开辟了新途径，为该罕见神经发育障碍的进一步机制研究奠定了基础。
研究背景与科学问题
自闭症谱系障碍（autism spectrum disorder, ASD）及相关神经发育障碍（neurodevelopmental disorders, NDDs）由多种遗传和环境因素共同驱动，影响早期皮层发育。近年人类转录组研究表明，RNA加工尤其是可变剪接（alternative splicing）是ASD和NDD风险的关键贡献者，因为动态剪接程序调控着皮质发生、神经元亚型特化、突触成熟和环路功能。然而，ASD中的RNA失调无法仅通过基因水平的表达变化充分解释，需要同时考虑转录本、外显子和剪接连接（splice junction）层面的结构变化。
在大脑发育过程中，剪接调控作为一个协调的发育网络运行，负责产生神经元特异性异构体、调节神经元兴奋性并在特定发育时期塑造突触信号。RNA结合蛋白（RNA-binding proteins, RBPs）如RBFOX、NOVA和异质性核糖核蛋白（heterogeneous nuclear ribonucleoproteins, HNRNPs）在胎儿皮质发生中调控精确时序的剪接转换。这些网络的破坏主要影响突触和投射相关基因程序，而这些程序在ASD和NDD中富集。具体而言，HNRNP H家族蛋白调控剪接位点选择并阻止G-rich基序处的异常连接，提示其破坏会导致连接聚焦的神经元异构体程序重组，而非单纯的外显子获得或丢失。
HNRNPH2的罕见致病性变异直接影响RNA加工机器，导致一种X连锁神经发育综合征，表现为发育迟缓、癫痫和高比例的自闭症特征。该基因编码HNRNP H家族蛋白，包含准RNA识别基序（quasi-RNA recognition motifs, qRRMs）和一个PY型核定位信号，促进与G-rich RNA的结合并通过karyopherin-β2通路进入细胞核。疾病相关的错义变异位于PY-NLS和qRRM结构域，预计会损害核运输或改变RNA结合特异性，从而导致HNRNP H/F靶向基序处剪接位点选择和连接使用的改变。尽管临床和机制证据强烈提示HNRNPH2与NDD和ASD相关，但在剪接程序和RBP网络活跃的组织相关环境中，从HNRNPH2变异到人皮层RNA表型的因果通路仍未完全阐明。
为此，研究人员利用同基因HNRNPH2-R114W和CRISPR校正的脑类器官，在遗传匹配背景下测量基因表达、连接和异构体水平的剪接以及推断的细胞类型组成，整合这些数据层以展示该变异如何改变人皮层转录组，描绘连接主导的剪接病（spliceopathy）及其与差异表达的关系，并识别可能为RBP失调相关ASD治疗策略提供信息的异构体水平终点。
研究方法概述
研究人员从携带HNRNPH2-R114W变异的患者外周血单个核细胞（peripheral blood mononuclear cells, PBMCs）中建立iPSC系，并通过CRISPR-Cas9在同基因背景下校正该变异获得对照系。采用改良Lancaster方案结合Stem Cell Technologies脑类器官试剂盒生成脑类器官，于分化后第20天收集每个基因型约12个类器官的混合RNA样本，每个基因型进行两次技术重复测序。RNA-seq数据分析包括：使用DESeq2进行差异表达分析；使用SGSeq进行连接聚焦的差异剪接分析；使用featureCounts结合voom/diffSplice进行外显子聚焦的差异使用分析；基于已知HNRNPH2结构模板（PDB: 1WG5, 1WEZ, 6DG1, 8SGH等）进行计算结构建模；使用MuSiC方法结合人神经类器官单细胞图谱（human neural organoid single-cell atlas, HNOCA）进行细胞类型去卷积分析。
研究结果
2.1 iPSC系的验证及HNRNPH2-R114W罕见变异背景下类器官发育的评估
免疫荧光证实两种基因型的iPSC均维持多能性标志物POU5F1（OCT4）和NANOG的表达。两种基因型均能生成具有中央腔和放射状排列祖细胞区的神经上皮结构，在形态学上无显著差异，PAX6阳性祖细胞和TUBB3阳性神经元分布正常。
2.2 HNRNPH2定向免疫荧光揭示基因型依赖性染色模式
在Plus类器官中，HNRNPH2信号主要定位于神经上皮区域，而Minus类器官呈现较弱且弥散的染色，提示突变导致分布或细胞定位改变。由于该检测仅依赖免疫荧光而无正交蛋白定量，这些发现被解释为HNRNPH2染色模式的变化而非蛋白丰度降低的直接证据。
2.3 HNRNPH2-R114W类器官中的全局转录重塑
主成分分析显示Minus与Plus样本沿PC1（解释99.5%方差）显著分离。差异表达分析鉴定出大量基因表达改变，其中Minus中下调基因富集于神经元分化、轴突发育和突触传递等神经发育过程，上调基因则富集于WNT/NOTCH信号和色素通路（包括黑色素生成、酪氨酸代谢和视黄醇生物合成）。
2.4 HNRNPH2-R114W类器官中广泛的可变剪接重连
计算结构模型将HNRNPH2的RRM1-3定位至剪接体B复合物附近，并将残基R114置于预测的RNA接触界面，R114W替换预计改变这一界面。连接聚焦分析鉴定出超过6,000个显著差异使用的连接（约占所有检测连接的6%-7%），而外显子聚焦分析仅发现17个显著差异外显子。这种连接主导的模式与HNRNP H家族蛋白结合G-rich RNA和调控剪接位点的已知功能一致。基因水平分析识别出九个在外显子使用上有显著变化的基因，包括CDH6、EFS、EVL、HNRNPC、UTRN、CRMP1、RUFY3、SPAG9和PTPRS，这些基因编码参与神经元突起、轴突发育、胞外基质组织和细胞骨架重塑的蛋白。
2.5 整合的RNA层面指标识别优先候选基因与功能趋同
研究人员开发了整合差异表达、事件水平和异构体水平证据的基因聚焦优先化框架。耦合强度热图揭示了不同基因在各层面的贡献程度差异，如CNTN2、PMEL、XIST、TUBA1A、CRABP1和CKB表现出事件主导、异构体主导或混合耦合模式。信号几何分析区分了表达变化大但RNA加工信号中等与表达变化小但RNA层面扰动强的基因。优先候选基因包括TRPM3、WNT2B、NR2F1、CLDN1、FLRT2、ALDH1A3、HMGCS1、DST和RTN4。功能富集分析显示这些基因主要涉及轴突导向、神经元突起形态发生、翻译和肽生物合成过程以及胞外基质和纤毛组织等结构程序。
2.6 HNRNPH2-R114W类器官中的细胞类型组成变化
MuSiC去卷积分析表明，Minus类器官中神经上皮和胶质丰富群体（成胶质细胞、少突胶质前体细胞）富集，而某些神经元群体轻度减少，平均神经元比例比Plus类器官低约5%。标志物基因热图支持这些推断的组成变化。然而，组成差异的相对较小幅度表明，主要的基因表达和剪接变化不能仅由细胞组成改变来解释。
讨论与结论总结
本研究证明，HNRNPH2-R114W罕见变异脑类器官表型的特征是早期神经元分化得以保持、广泛而有组织的表达变化、连接主导的剪接重连以及关键基因枢纽的汇聚，细胞组成变化作为修饰因子而非分化的主要驱动因素。这些发现支持如下模型：HNRNPH2变异通过改变蛋白定位、剪接和核糖核蛋白组装来破坏转录后调控，产生毒性或显性负效应而非简单的功能丧失。PY-NLS和qRRM突变可能损害核输入和RNA结合，从而限制HNRNPH1的代偿并诱导神经剪接网络中的毒性功能获得效应。该研究结论指出，类器官模型为这一观点提供了实验支持，可作为评估等位基因选择性反义寡核苷酸策略、HNRNPH2调控或旁系同源补偿方法的平台。这种整合差异表达、多层剪接分析和基于单细胞的批量RNA-seq去卷积的方法学，很可能适用于其他RBP驱动的神经发育障碍。
