《Biologia Futura》：Comparative analysis of conserved non-coding elements identifies gene regulatory networks rewired during the water-to-land transition in vertebrates
编辑推荐：
脊椎动物征服陆地是进化史上的一个关键时刻。适应新的栖息地 necessitates 脊椎动物解剖学和生理学的众多变化，在四足动物的发育基因调控网络(GRNs)上留下了持久的印记。过去十年间高质量基因组资源的增加，使得研究水生向陆生转变的基因组遗产成为可能。虽然
脊椎动物征服陆地是进化史上的一个关键时刻。适应新的栖息地 necessitates 脊椎动物解剖学和生理学的众多变化，在四足动物的发育基因调控网络(GRNs)上留下了持久的印记。过去十年间高质量基因组资源的增加，使得研究水生向陆生转变的基因组遗产成为可能。虽然人们将大量注意力集中在基因组中高度保守的非编码元件(CNEs)上，这些元件在进化上分化的支系间共享高度相似性，但最近的证据表明，或许应该给予“缺失的”CNEs以相当的关注——这些保守序列片段存在于现存的有颌类干群和辐鳍鱼类中，但在四足动物适应陆地生活的过程中变得无法检测，无论是通过真实的序列丢失还是分歧超出了可比对性。这些序列可以帮助研究人员揭示与水生生活方式相关的某些发育约束的松弛，这些约束的松弛使得在发育景观中达到新的适应性峰值成为可能。在本论文中，研究人员以斑马鱼(Danio rerio)基因组为参考，搜索此类CNEs并将其与泛有颌类CNEs进行比较表征。研究结果表明，与色素沉着和肌肉结构形成相关的发育网络重布线留下了最大的基因组印记。研究人员还发现，典型Wnt和Hedgehog信号通路的组分在鱼类中保留的CNEs中富集。
**论文解读：脊椎动物水生向陆生转变过程中基因调控网络重布线的保守非编码元件比较研究**
**一、研究背景与科学问题**
脊椎动物从水生环境征服陆地是进化史上具有里程碑意义的事件。这一生态位的转变 necessitated 脊椎动物解剖学和生理学上的大量革新，包括角质化皮肤的出现以防止水分流失、渗透调节方式的转变（鳃失去渗透调节功能而肾脏作用增强）、造血部位从肾髓向骨髓的迁移，以及头部结构的广泛重塑。最为重要的变革体现在中轴骨骼的转化和鳍到四肢的转变上：鱼类特有的分节W形肌节被四足动物的多节段肌肉束所取代，快缩糖酵解纤维为主向慢氧化纤维与浅层快缩纤维的布局转变，鳍条被能够承托体重的内软骨性自肢所替代。
这些深刻的解剖学变革得以实现，依赖于基因组层面的转型变化——基因内容的重排与发育基因调控网络(GRNs)的重布线。尽管过去的许多研究集中于基因组中高度保守的非编码元件(CNEs)在进化中的保守性，但近期证据表明，那些在现存有颌类干群和辐鳍鱼类中保守、却在四足动物中“缺失”的序列片段同样值得关注。这些“缺失”的CNEs可能揭示了与水生生活方式相关的发育约束的松弛，正是这种松弛使得四足动物在发育景观中达到新的适应性峰值成为可能。然而，系统性地鉴定并表征这类鱼类特异性CNEs的研究仍然缺乏。基于此，研究人员开展了本项研究，旨在通过比较基因组学方法识别和分析这些在四足动物中丢失的古老CNEs。
**二、主要技术方法与数据来源**
研究人员构建了基于snakemake的定制化分析流程，对7个物种的基因组进行了全基因组比对分析。样本队列包括：四足动物——人(Homo sapiens)、鸡(Gallus gallus)、西爪蟾(Xenopus tropicalis)；辐鳍鱼类——斑马鱼(Danio rerio)、青鳉(Oryzias latipes)、斑点雀鳝(Lepisosteus oculatus)；以及外群有颌类——小斑猫鲨(Scyliorhinus canicula)。所有基因组数据均从NCBI获取。核心方法包括：使用progressive cactus比对工具进行全基因组多序列比对；基于4D简并位点构建中性替换模型；运用phastCons软件分别在有颌类、辐鳍鱼类和肉鳍鱼类三个系统发育树上鉴定保守元件；整合ATAC-seq和RNA-seq数据（来源于Bogdanovi?等人24小时斑马鱼胚胎的公开数据）评估开放染色质状态和基因活性；利用rGREAT进行基因本体(GO)富集分析；基于SwissRegulon数据库的转录因子结合基序进行motif富集分析。
**三、主要研究结果**
**1. CNE的鉴定与基本特征**
通过比较基因组分析，研究人员共鉴定出3,227个有颌类保守(Gnath-cons)CNEs和626个辐鳍鱼类特异性(Acti-spec)CNEs，并成功定位到斑马鱼基因组上。染色体分布分析显示，4号染色体右臂几乎完全缺乏CNEs——这与该区域包含快速进化的性别决定位点和最大的异染色质区块的解释一致；14号染色体上的Acti-spec CNE富集峰与原钙粘蛋白簇重合；12号和24号染色体上的Gnath-cons CNE热点则位于已知的超保守非编码元件中心（ebf3a、esrrga和zfhx4基因附近）。
**2. 序列与基因组特征比较**
长度分析显示Acti-spec CNEs的中位长度显著长于Gnath-cons CNEs（81 bp对57 bp）。尽管Acti-spec CNEs仅局限于较短的进化窗口，但其保守性分数（中位phastCons值18）与深度保守的Gnath-cons集合（中位值20）相当，表明这些元件确实受到真实的纯化选择作用。两组CNEs在基因组的分布模式上与已知的ATAC-seq峰和转录活跃基因相关联，证实了这些元件是真正的顺式调控序列而非比对伪迹。
**3. CNE热点基因分析**
大多数基因仅携带单个CNE（76%的Acti-spec基因和59%的Gnath-cons基因），CNE热点代表真正的离群值。Acti-spec数据集中排名首位的是原钙粘蛋白基因pcdh2g1（13个CNEs），其他热点包括多种锌指和叉头转录因子(zfhx4、esrrga、zfhx3b、znf423、foxp1b、sall1a、znf536)以及细胞骨架和粘附基因(macf1a、pcdh2aa1、rbfox1、nrxn3a)。转录因子在两组热点基因中均占相当比例（Acti-spec的top 15中占47%，Gnath-cons的top 15中占33%）。
**4. 基因本体富集分析**
比较GO分析揭示了若干重要的差异模式。发育色素沉着和色素细胞分化在两组CNEs中均富集但显著偏向Acti-spec（发育色素沉着：FE=3.7对1.6；色素细胞分化：FE=4.0对1.9），色素沉着本身仅在Acti-spec数据集中显著。在信号通路层面，典型Wnt信号通路及其调控几乎仅在Acti-spec数据集中富集（FE=3.3），Hedgehog/Smoothened通路组分同样仅在该数据集中显著富集；而激活素受体信号通路和视黄酸信号通路仅在Gnath-cons数据集中显著。肌肉器官发育和肌肉细胞分化在两组中均富集但Acti-spec略高。鱼鳍和附肢形态发生在两组中具有相近的富集度，而更广泛的骨骼系统发育和骨矿化在Gnath-cons中显著性更高。
**5. 转录因子结合基序富集分析**
motif扫描结果显示，仅有一个基序——眼区转录因子对Six6a/Six6b在Acti-spec CNEs中显著富集，而Gnath-cons CNEs中有16个基序显著富集，包括眼区因子Rx1/Rx2、Pdx1、Mecp2、Neurod1、多个Hox因子和Myod1等。这一不对称性提示鱼类特异性CNE层虽然在功能上具有一致性，但并非由一套相对应的统一转录因子所调控。
**四、讨论总结与研究结论**
本研究的核心发现是：驱动基础发育过程的基因相关CNEs在有颌类中深度保守，但相同的发育程序通常还携带来自第二层鱼类特异性CNEs的调控输入。这一鱼类特异性调控层分布于大多数所检视的发育程序中。特定过程和信号通路表现出一致的偏向性：Wnt和Hedgehog信号偏向鱼类特异性CNE层，激活素和视黄酸信号偏向深度保守层，FGF受体信号在两层中均有代表性。鱼类特异CNE层中Wnt信号的富集与斑马鱼鳍再生的芽基形成能力相关，而这一再生能力在四足动物肢体中已丧失；色素沉着方面的偏向与鱼类拥有多种神经嵴来源色素细胞类型而四足动物仅保留黑色素细胞谱系的事实一致。内耳方面，耳囊形成和基板诱导在Acti-spec中富集度更高，而耳石发育和半规管发育在Gnath-cons中更显著。
研究结论指出，与色素沉着和肌肉结构形成相关的发育网络重布线在水生向陆生转变过程中留下了最大的基因组印记。需要强调的是，“缺失”这一操作定义并不区分真实的序列删除和分歧超出检测阈值。斑马鱼作为比对参考的选用引入了以鱼类为中心的系统偏差，且当前流程无法单独区分上述两种场景。尽管如此，本研究的预测旨在作为假设生成的资源：Actinopterygii特异性与Gnathostome保守CNEs之间的候选调控差异，为在实验层面测试脊椎动物GRN重布线提供了具体位点。其中Six6a/Six6b位点具有直接的转化医学相关性——人类SIX6是原发性开角型青光眼的已验证遗传风险位点，解析其周围顺式调控元件的保守与谱系特异性，有助于更好地利用斑马鱼模型研究人类青光眼和小眼畸形。
