《Organogenesis》：Functional characterization of a non-coding RUNX1T1 variant in a kindred with esophageal atresia
编辑推荐：
食管闭锁（EA）是一种先天性前肠畸形，具有显著的遗传异质性。尽管编码区变异可解释部分综合征型病例，但非综合征型EA的分子基础在很大程度上仍不清楚。研究人员调查了一个非综合征型EA高发的以色列贝都因家系。对10例受累个体中的3例进行了全外显子组/全基因组测序。经
食管闭锁（EA）是一种先天性前肠畸形，具有显著的遗传异质性。尽管编码区变异可解释部分综合征型病例，但非综合征型EA的分子基础在很大程度上仍不清楚。研究人员调查了一个非综合征型EA高发的以色列贝都因家系。对10例受累个体中的3例进行了全外显子组/全基因组测序。经变异过滤并利用人群数据库和计算机模拟工具进行优先级排序后，对候选非编码变异在所有5例可获得的受累个体中进行了Sanger测序。通过小基因剪接试验、定量PCR、双荧光素酶报告基因试验和果蝇过表达模型评估功能效应。对3例受累个体的测序数据变异分析鉴定出一个共享的杂合深内含子RUNX1T1变异（chr8:91994544G > A，hg38），随后在5例参与研究的受累个体中得到验证。尽管该变异被纳入一个罕见的RUNX1T1转录本，但并未影响剪接。相反，它增加了经典RUNX1T1亚型的表达并增强了报告基因活性，与调控性功能获得效应一致。RUNX1T1直系同源基因（nvy）在果蝇中的全身过表达导致摄食器官的形态学改变，包括硬度降低和上唇感觉器官伸长。在一个表现为不完全外显的显性非综合征型EA家系中，研究人员鉴定出一个新型非编码RUNX1T1变异，并具有调控性功能获得效应的功能证据。这些发现提示RUNX1T1失调是EA发病机制的贡献因素，并强调了非编码调控变异在先天性前肠畸形中的重要性。
研究背景方面，食管闭锁（EA，OMIM 189960）是一种罕见的先天性畸形，发病率约为每3000例活产中1例，以食管连续性中断为特征，常伴有气管食管瘘（TEF）。EA被认为源于早期前肠发育紊乱，即呼吸管与消化管分离形成腹侧气管和背侧食管的過程受到干扰。根据Gross分型，EA分为A型（无TEF，8%）、B型（伴近端TEF，3%）、C型（伴远端TEF，85%）、D型（伴近端和远端TEF，>1%），E型（4%）为非综合征型TEF不伴EA。EA出生后面临无法喂养、食管狭窄和运动障碍、吸入性肺炎等危及生命的风险，治疗以手术为主。EA和TEF属于遗传异质性疾病，仅约11%–12%的病例确定了遗传病因，涵盖染色体非整倍体、拷贝数变异和罕见单基因综合征。尽管少数家系可见孤立性孟德尔遗传，但以散发病例为主、同胞复发风险低、单卵双胎一致性增高及显著的表型变异性提示大多数受累个体具有强遗传基础但未必为单基因遗传。已明确遗传学病因的病例具有一个共同特征即基因剂量敏感性，杂合功能缺失变异或拷贝数改变破坏精密调控的发育程序，如转录和染色质调控因子SOX2、CHD7、MYCN和ZIC3，以及涉及SHOX、ABLIM1和LPP的复发性拷贝数变异。结合广泛的解剖表现谱和较低的总体诊断率，这些观察支持多基因和多因子疾病架构。在此模型中，多个遗传变异各自具有独立效应，以外显率可变的方式与环境、随机或发育修饰因素相互作用决定疾病表现。鉴于该大家系内EA的家族性发生，研究人员开展了基因组和功能研究以探究促成该表型的潜在遗传因素。研究纳入一个大型隔离以色列贝都因家系中10例诊断为非综合征型C型EA的个体。临床遗传学评估未发现提示可识别综合征型EA的额外共享表型特征，尤其无受累个体符合VACTERL联合征或CHARGE综合征的临床标准，亦无一致的面部畸形或发育迟缓/智力障碍。但部分受累个体观察到可变的心脏和脑部异常。研究结论为：在一个表现为不完全外显的显性非综合征型EA家系中，研究人员鉴定出一个新型非编码RUNX1T1变异，并具有调控性功能获得效应的功能证据；这些发现提示RUNX1T1失调是EA发病机制的贡献因素，并强调了非编码调控变异在先天性前肠畸形中的重要性。研究意义在于将RUNX1T1确定为可能对EA发生起贡献作用的候选遗传因子，并鉴定出一个先前未注释的内含子顺式调控元件，其活性被所研究变异增强，呈现调控性功能获得机制，为前肠形态发生调控因子在EA发病机制中的作用增添了新的遗传学证据。
关键技术方法方面，研究人员对3例受累个体行全外显子组测序（WES）和全基因组测序（WGS），经VARista平台过滤并以CADD评分优先级排序，采用Sanger测序在家系成员中验证变异；通过小基因剪接试验（minigene splicing assay）联合丝氨酸/精氨酸富集剪接因子（SRSF1、SRSF2、SRSF5、SRSF6）过表达及定量PCR评估变异对选择性剪接的影响；采用双荧光素酶增强子报告基因试验评估该区域的顺式调控活性；并构建果蝇（Drosophila）UAS-nvy（RUNX1T1直系同源基因）全身过表达模型，利用扫描电子显微镜（SEM）进行形态学分析。样本来源为以色列贝都因家系，仅5例受累个体参与分子研究。
研究结果方面，临床特征分析显示，该扩展家系多个分支共10例EA个体中5例参与临床和分子研究，均为C型EA，部分个体伴有可变先天畸形，未发现额外共享综合征表型；家系中受累个体分布最符合常染色体显性遗传伴不完全外显，估计外显率约40%。遗传分析显示，3例受累个体WES/WGS分析未发现EA/TEF已知相关基因的变异，亦无共享的编码候选变异；扩展至非编码变异分析后，保留25个共享杂合变异，其中8个CADD评分超过15，RUNX1T1深内含子变异（hg38: chr8:91994544G > A; rs112217168）排名最高，该变异在5例受累个体中经Sanger测序验证；该变异位于绝大多数RUNX1T1亚型的深内含子区，但被纳入独特"esoteric"亚型ENST00000520978.5的第2外显子。g.8:91994544G > A变异上调经典RUNX1T1亚型转录，小基因剪接试验联合剪接因子过表达的定量PCR显示esoteric亚型表达未受变异影响，而经典RUNX1T1亚型表达在突变组显著增加约2倍，提示该变异不影响剪接而是对经典RUNX1T1转录产生调控性功能获得效应。g.8:91994544G > A变异所在顺式调控元件具有功能获得效应，双荧光素酶报告基因试验显示含变异序列区域较空载显著提高转录输出，引入变异后转录活性进一步增强。果蝇nvy全身过表达导致摄食器官形态改变，actin5C-GAL4驱动的nvy过表达果蝇虽发育至形态成熟成虫阶段但未能羽化而死于蛹壳内，扫描电子显微镜观察显示摄食器官（喙、吻突和唇瓣）呈现松弛、硬度降低及表面质地改变，定量形态学分析显示上唇感觉器官（LSO）显著伸长。
讨论部分总结：贝都因人群高度近亲结婚导致变异等位基因频率增高，所研究的部分外显常染色体显性修饰变异即体现此类人群特异性遗传动态，按ACMG/AMP框架归类为意义不明变异（VUS），反映了当前框架解读非编码调控变异的局限。研究人员报告了非编码RUNX1T1变异与扩展近亲家系中EA的关联，体外功能研究证明该深内含子变异位于调控元件内并上调RUNX1T1表达，结合动物模型中RUNX1T1直系同源基因操作所观察到的表型，提示RUNX1T1为可能贡献EA病因的候选基因；但现有数据未确立该变异与EA的因果关系，也未提供RUNX1T1与C型EA特异性关系的证据。家系中2例受累个体观察到脑畸形，其与EA表型或该变异的关系未知；鉴于人群和家系内高度近亲结婚，部分成员的心血管和脑异常可能代表共存的遗传疾病而非EA或该杂合变异的相关特征。功能分析表明该变异位于先前未识别的顺式调控元件内，不影响剪接或编码序列；虽明确显示变异上调RUNX1T1表达，但不能排除其通过局部调控效应影响邻近基因（PIP4P2、OTUD6B、LRRC69、SLC26A7和TRIQK）表达的可能，唯这些邻近基因的表达模式或已知功能与EA表型无明显相关性。Runx1t1缺陷小鼠的中肠闭锁为变异优先级排序提供了生物学支持，但不应解释为中肠闭锁与EA共享共同发育机制的证据，因前肠和中肠经历不同的区域特化和形态发生，中肠闭锁形态不能重现C型EA。果蝇中nvy过表达后LSO伸长不应解释为人类EA的解剖或发育对应物，LSO为果蝇摄食器官的感觉组件，与脊椎动物食管非同源。nvy全身过表达致死原因尚不确定，摄食器官异常不能假定解释致死表型。RUNX1T1编码转录辅阻遏物，主要在急性髓系白血病（AML）中研究，涉及t(8;21)产生致病性RUNX1–RUNX1T1融合蛋白；受累个体及家系分支中未报告AML或其他血液系统恶性肿瘤，EA与AML关联证据亦不明确。RUNX1T1天然形式与RBPJ和SHARP组装成复合物，抑制Notch信号通路靶基因转录，并可能作为表观遗传修饰因子。迄今与EA相关的基因多为综合征型，均表现为不完全外显，与本队列非综合征型RUNX1T1相关EA一致；增强子胚系变异常与可变外显率/表达相关，与本研究发现相符。基于已知综合征型EA遗传修饰因子，可尝试勾勒RUNX1T1功能获得致EA的分子通路：Sox2在小鼠中被证明为Notch靶基因，由Notch胞内结构域直接激活并被RBPJ复合物抑制，而RUNX1T1与RBPJ和SHARP组装抑制Notch靶基因，故RUNX1T1功能获得可能下调SOX2表达驱动EA发生；Nkx2-1在气管食管形态发生中起关键作用，Notch胞内结构域过表达可减少斑马鱼nkx2-1表达空间域，故RUNX1T1功能获得可能影响nkx2-1表达。CHD7与RBPJ直接相互作用，RUNX1T1亦直接影响MYCN活性。上述数据提示RUNX1T1活性改变与EA的可能分子通路，但确切机制尚待阐明，其在中胚层和血管发育中的作用（调控VEGFA、BMP4和TGF-β2表达）可能提供潜在联系。研究利用多种先进基因组、生物信息学和功能方法阐明所研究变异的分子后果，体现整合性变异解读的潜力，强调临床医生与基础研究者紧密合作的价值； cohort规模较小和外显不全限制了基因-疾病关联强度，但数据证明该变异对基因表达具有明确上调效应，鉴定出一个先前未注释的增强子或顺式调控元件，该效应可能源于阻遏物结合位点破坏、激活物结合基序引入或激活物位点亲和力增强。该调控元件的时空表达模式及变异改变其活性的方式尚待确定，斑马鱼体内转基因增强子报告试验未显示野生型或突变型构建体可检测活性。现有功能证据提示该元件主要在多种癌细胞系及若干胚胎组织中活跃，在至少三种患者来源白血病细胞系（包括携带t(8;21)并表达RUNX1–RUNX1T1融合蛋白的Kasumi-1）中观察到调控活性，值得进一步研究其在肿瘤发生尤其是AML中的潜在作用。研究局限性包括：因医疗和遗传服务可及性有限及参与研究性基因检测意愿有限，10例EA个体中仅5例参与分子研究；尽管近亲结婚程度高，受累个体分布不符合简单常染色体隐性遗传模式，提示显性伴不完全外显；仅少数未受累亲属可能可及采样，无法可靠估计变异外显率。
