等位基因特异性CRISPR扰动印迹Dlk1–Dio3结构域揭示BMP–NOTCH–VEGF信号通路在胚胎器官发生中的调控机制

《CELL DEATH AND DIFFERENTIATION》：Allele-specific CRISPR perturbation of the imprinted Dlk1–Dio3 domain reveals regulation of BMP–NOTCH–VEGF signaling in embryonic organogenesis

　　

印迹Dlk1–Dio3结构域对哺乳动物发育至关重要，但其协调多器官胚胎发生的作用尚未完全阐明。研究人员利用CRISPR/Cas9介导的转录终止盒插入技术，构建了涵盖纯合子（HOMO）、母源遗传（MK）、父源遗传（PK）及野生型等位基因（WT）的四种小鼠模型，以阐明亲本来源特异性效应。跨尺度的整合转录组学分析揭示了印迹基因的等位基因特异性失调以及该位点表观遗传稳定性的改变。至关重要的是，对胚胎期14.5天（E14.5）肝脏、心脏和胎盘的单细胞RNA测序（scRNA-seq）显示，Dlk1–Dio3纯合失活或母源靶向的胚胎因BMP–NOTCH–VEGF信号网络崩溃而发育停滞。这表现为器官特异性病理改变：（i）蜕膜-滋养细胞旁分泌通讯受损及BMP/WNT信号通路异常，损害母胎血屏障，从而破坏胎盘血管生成并限制母胎营养交换；（ii）胎肝造血祖细胞表现髓系谱系定向障碍，伴随NOTCH和WNT信号异常，导致肝细胞分化失败；（iii）心脏祖细胞未获得充足的BMP–NOTCH–VEGF信号，产生未成熟心肌细胞及血管缺陷。综合而言，这些发现将Dlk1–Dio3印迹结构域确定为器官发生过程中细胞类型特异性信号通路的关键整合枢纽，并为基因组印迹失调如何引发系统性胚胎发育失败提供了机制性见解。

基因组印迹（genomic imprinting)以亲本特异性单等位基因表达为特征，是胎盘哺乳动物发育的基础性表观遗传机制。Dlk1–Dio3印迹结构域位于小鼠12号染色体及人类14号染色体，跨越约1 Mb区域，包含三个父源表达蛋白编码基因（Dlk1、Rtl1、Dio3）及若干母源表达非编码RNA，如Gtl2/Meg3、Rian、Mirg和一个庞大的微小RNA（microRNA）簇。既往研究表明，Dlk1印迹扰动影响出生后生长及代谢过程，而Rtl1敲除导致胚胎致死或异常。尽管该位点的基因组架构已被明确界定，其在协调多器官发育中的系统性作用仍不清楚。该结构域的双组分印迹调控区域（imprinting control region, ICR）——包括初级基因间差异甲基化区域（IG-DMR）和次级Gtl2-DMR——严格调控其等位基因剂量，而这一剂量的精确维持对不同发育微环境至关重要。在母胎界面，母源microRNA（如miR-127）与父源Rtl1的化学计量平衡驱动迷路区形态发生及毛细血管网络完整性；在胎盘之外，该位点对肝脏及心脏器官发生行使多效性调控。然而，胎盘、肝脏和心脏并非独立发育，而是形成相互关联的发育轴；Dlk1–Dio3结构域作为层级遗传枢纽整合这些多器官信号网络，其复杂跨器官串扰的解码对理解印迹缺陷驱动的广泛发育病理至关重要。

本研究旨在解决的核心问题在于：传统组织整体分析掩盖细胞异质性并遮蔽稀有亚群，而单细胞RNA测序（scRNA-seq）技术虽能精确描绘细胞类型特异性转录程序，但将其与遗传工程动物模型整合以阐明特定遗传扰动如何破坏正常发育细胞机制的研究仍待深入。研究人员此前已建立并表征了Gtl2转录终止小鼠模型，鉴定出纯合子、母源及父源遗传敲入小鼠在胎盘、肝脏和心脏中的不同表型缺陷，但Dlk1–Dio3失活与多器官失败之间的精确细胞机制及细胞间信号破坏尚不明确。

为探究上述问题，研究人员采用CRISPR/Cas9介导的印迹特异性敲入等位基因技术，结合多器官scRNA-seq，构建了HOMO、MK、PK及WT小鼠模型，通过在Dlk1–Dio3区域插入多聚腺苷酸化转录终止子实现等位基因特异性转录终止而不修饰相邻基因组调控序列。表型分析显示，HOMO和MK突变体在E14.5表现胚胎致死，而PK和WT胚胎正常发育，凸显该位点母源表达转录本对胚胎存活的关键功能。

在主要技术方法方面，研究人员采用Easi-CRISPR策略在Gtl2基因第一外显子插入含三个多聚腺苷酸信号（3×polyA）的转录终止盒，获得四种基因型小鼠；通过体视显微镜观察胚胎形态学表型；对E14.5胎盘及胎儿进行整体转录组RNA测序（RNA-seq）；对四种基因型肝脏、心脏和胎盘各12个样本进行scRNA-seq（基于Singleron Matrix平台及GEXSCOPE文库构建，10x Genomics兼容数据）；运用Seurat进行数据预处理与聚类，Harmony消除批次效应，monocle3进行拟时序轨迹分析，CellChat进行细胞间通讯分析，pySCENIC进行转录因子活性推断，并采用基于蛋白质-蛋白质相互作用（PPI）网络的新型基因网络传播算法识别细胞类型特异性发育基因。

研究结果部分如下：

**印迹特异性Dlk1–Dio3转录终止模型的构建**：polyA盒插入导致母源表达转录本（Gtl2、Rian、Mirg）沉默，同时父源表达基因（Dlk1、Rtl1）在MK和HOMO样本中显著上调，证实印迹丢失。HOMO和MK胚胎于E14.5出现生长停滞伴显著胚胎内出血，而PK及WT同窝对照表型正常。批量转录组分析显示MK与WT间差异表达基因（DEGs）数量最多（n=3,194），且致死基因型中DEGs在Dlk1-Dio3位点（chr12qF1）附近呈簇状失调。

**单细胞图谱揭示胎盘滋养细胞及蜕膜分化缺陷**：E14.5胎盘scRNA-seq鉴定出滋养细胞、海绵滋养细胞、蜕膜细胞、基质细胞、内皮细胞、髓系细胞、红细胞及免疫细胞等群体。HOMO/MK胎盘滋养细胞比例增加而蜕膜细胞比例改变。网络传播算法鉴定出Aebp1（蜕膜）和Tnnt1（滋养细胞）为最高权重标志基因。拟时序分析显示Dec_Aebp1_high+成熟状态亚群在HOMO/MK中提前分化停滞，而Trop_Tnnt1_high+亚群在早期拟时相异常积累后随发育时间显著减少。

**致死等位基因相关胎盘中信号通路失调**：SCENIC推断的转录因子活性显示Dec_Aebp1_low-亚群具有高Hand1、Pparg和Irf8活性，其中Hand1调控BMP/WNT信号轴。细胞间通讯分析揭示HOMO/MK胚胎NOTCH相互作用减弱而VEGF、IGF和FGF级联异常激活。KEGG通路分析显示下调基因富集于"核糖体"和"造血细胞谱系"通路，GSEA证实血管生成及核糖体生物发生通路显著富集。

**肝脏和心脏祖细胞发育延迟**：肝脏scRNA-seq鉴定出肝细胞、红细胞、共同髓系祖细胞/造血干细胞（CMP/HSC）、巨核细胞、中性粒细胞、成纤维细胞和库普弗细胞。HOMO/MK肝脏CMP/HSC比例及绝对数量降低（4.35% vs. 9.17%），肝细胞和巨核细胞略减。网络传播算法鉴定Ifitm1（CMP/HSC）、Serpina1b（肝细胞）和Itgb3（巨核细胞）为关键标志。CMP/HSC_Ifitm1_high+亚群在MK中比例升高（73.6% vs. 43.3%），拟时序分析显示该亚群及Hepa_Serpina1b_low+细胞发育停滞于早期阶段。转录因子分析显示Nfyb、Ubp1、Hnf4g、Etv5和Myc活性异常。细胞间通讯显示Cstg-F2r、Cstg-Pard3和Thbs1-Cd47信号破坏。

心脏分析显示心肌细胞（Card）比例增加而平滑肌细胞（SMC）锐减（1.67% vs. 7.35%），心外膜细胞减少。Myl3、Sox17、Myrf和Adamts2分别作为Card、心脏内皮细胞（Cendo）、心外膜细胞（Epicard）和SMC的标志基因。Card_Myl3_high+亚群比例异常增加，SMC_Adamts2_high+亚群降至27.1-32.4%（对照>52%）。Myl3表达随发育时间未能增加，Myrf表达缺失阻碍心外膜-间充质转化（EMT）。心肌细胞中转录因子Gata4、Tbx5、Ets1、Sox11和Hey2活性降低。BMP和VEGF通路相互作用显著减弱，而PK/WT中Cendo_Sox17_high+与SMC_Adamts2_high+群体间NOTCH和ANGPT信号对心室小梁化至关重要。组织学显示HOMO心室发育不良、心肌紊乱，MK心外膜层破坏；胎心率检测证实HOMO/MK严重心动过缓（HOMO 23.7±4.1 bpm；MK 25.3±3.8 bpm vs. WT 39.1±4.9 bpm）。

**多器官信号网络的系统性失败**：整合分析显示胎盘Dec_Aebp1_high+和Trop_Tnnt1_high+分化停滞与肝脏CMP/HSC_Ifitm1_high+祖细胞扩增及心脏心肌细胞成熟阻滞同步发生。差异细胞间相互作用分析揭示BMP/NOTCH/VEGF跨器官信号网络崩溃：正常条件下胎盘滋养细胞分泌VEGF和IGF支持胎儿器官发生，而HOMO/MK中这些输出减弱；BMP信号从滋养细胞至海绵滋养细胞减弱，ANGPT和HSPG信号至肝脏祖细胞和心脏内皮细胞也减弱。肝脏未能向胎盘提供VEGF信号和ANGPTL信号，NOTCH/ANGPTL轴至心脏谱系沉默；心脏未能向胎盘发送CXCL/VEGF/ANGPTL反馈信号。这种级联失败导致胎盘BMP/VEGF/ANGPTL输出减少，胎儿肝脏和心脏失去生长信号，促进未成熟CMP/HSC积累并导致肝脏和心脏信号枢纽崩溃，加剧器官广泛缺氧应激，最终引发多器官衰竭。

在讨论部分，研究人员指出其创新性方法框架将等位基因特异性基因组编辑与多器官单细胞转录组学整合，优于可能导致大规模基因组缺失混肴效应的传统敲除模型。多模态分析表明Dlk1–Dio3印迹破坏的后果超越胎盘缺陷，反映协调器官发生的系统性失败。该结构域作为胚胎中期发育的主调控因子，通过BMP–NOTCH–VEGF回路维持胎盘-肝脏-心脏信号环路；印迹剂量平衡的破坏触发该信号网络进行性崩溃，导致多器官发育失败。这一发现对理解人类印迹疾病及精准医学具有转化意义，为等位基因特异性治疗恢复疾病中的正常印迹提供了新方向。

研究结论总结如下：研究人员的数据确立了印迹Dlk1–Dio3结构域为胚胎中期发育的主调控因子，其通过BMP–NOTCH–VEGF回路发挥血管和器官形成的基础作用。无法维持该印迹剂量导致系统性器官衰竭。该研究不仅阐明了印迹的基础发育生物学，也为表观遗传疾病和基因治疗中的转化进展指明了令人振奋的路径。