目录
研究背景：为什么关注蝙蝠 HoxD 增强子体内功能验证
在发育生物学和演化发育学研究中，一个核心问题是：基因组中的非编码调控序列如何影响器官形态演化？蝙蝠前肢演化为翼，是研究肢体发育、Hox基因簇调控和增强子功能变化的典型模型。相比只观察序列差异，研究人员更需要回答的是：这些候选非编码序列是否能在胚胎体内真实驱动组织特异性表达。
本篇365英国上市集团生物 文献解读聚焦PLOS Genetics论文“Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus”。该研究将比较基因组学与转基因小鼠增强子检测结合起来，在E12.5小鼠胚胎中评估候选蝙蝠加速区域（bat accelerated regions，BARs）的Hsp68-LacZ报告活性，从而把“加速演化信号”转化为可观察的发育表达证据。
文献发表于《PLOS Genetics》
文献概览：PLOS Genetics 研究与365英国上市集团相关模型信息
| 文献题目 | Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus |
|---|---|
| 发表信息 | Booker BM 等，PLOS Genetics，2016，12(3): e1005738；DOI: 10.1371/journal.pgen.1005738 |
| 365英国上市集团相关模型/服务 | 原文方法学说明，所有转基因小鼠均由 Cyagen Biosciences 按标准流程构建。 |
| 核心用途 | 用于Hsp68-LacZ小鼠增强子检测，在E12.5胚胎中验证候选蝙蝠加速区域（BAR）的体内组织特异性增强子活性。 |
该研究围绕蝙蝠前肢如何演化为翼这一经典问题，结合多物种基因组比较、发育期小鼠转基因增强子检测和胚胎原位杂交，寻找可能参与蝙蝠前肢特异性发育的非编码调控元件。
研究的核心发现是：研究人员鉴定出166个蝙蝠加速区域。这些区域在多数脊椎动物中高度保守，却在蝙蝠祖先分支中出现快速序列变化。进一步的体内增强子检测显示，5个来自小棕蝠（Myotis lucifugus）的BAR均可在E12.5小鼠胚胎肢体中驱动LacZ表达。其中，位于HoxD基因簇端粒侧的BAR116最为突出，蝙蝠序列能够驱动明显的前肢表达，而小鼠同源序列及嵌合序列未能稳定复现这一模式。
这篇论文的亮点不只是发现新的候选调控序列，更重要的是把基因组层面的加速演化信号与胚胎体内的组织表达模式联系起来。通过转基因小鼠胚胎中的LacZ报告系统，研究者将候选非编码序列的功能差异转化为可观察、可比较的发育证据。
研究材料与技术方法：从BAR候选序列到Hsp68-LacZ报告系统
该研究的材料和方法围绕“候选调控区域发现—体内增强子检测—发育表达验证”展开。比较基因组学用于筛选在脊椎动物中保守、但在蝙蝠祖先分支中加速演化的非编码区域；Hsp68-LacZ转基因报告系统则用于判断候选序列是否具备体内增强子活性。
| 研究环节 | 关键内容 | 验证目的 |
|---|---|---|
| 比较基因组学数据 | 纳入多物种基因组比对数据，筛选脊椎动物保守、蝙蝠分支加速演化的非编码区域 | 发现候选蝙蝠加速区域（BAR） |
| 候选序列 | 重点选择BAR2、BAR4、BAR61、BAR97、BAR116，并检测部分小鼠同源序列 | 比较蝙蝠序列与小鼠同源序列的功能差异 |
| 转基因报告系统 | 候选增强子序列被克隆至Hsp68-LacZ报告载体 | 检测候选序列是否能在体内驱动组织特异性表达 |
| 动物模型 | 单细胞小鼠胚胎注射后获得E12.5转基因小鼠胚胎 | 在真实胚胎发育环境中观察LacZ表达模式 |
| 表达验证 | Miniopterus natalensis蝙蝠胚胎、野生型Parkes小鼠胚胎及Hoxd10-Hoxd13原位杂交探针 | 评估BAR116是否可能参与HoxD基因簇发育调控 |
在技术方法上，研究团队先通过PhastCons识别脊椎动物保守序列，再利用PhyloP评估蝙蝠祖先分支中的加速演化信号。同时，研究整合发育期小鼠肢体H3K27ac和p300 ChIP-seq峰，以富集潜在肢体发育增强子。候选BAR序列或小鼠同源序列被克隆至Hsp68-LacZ载体，注射至单细胞小鼠胚胎后，在E12.5阶段进行LacZ染色和组织表达观察。
对于关注非编码调控元件、胚胎期组织特异性表达或LacZ报告系统的研究团队，这类动物模型定制服务能够帮助将序列预测结果转化为体内表达证据。若研究目标涉及报告基因插入或基因位点标记，也可结合基因敲入小鼠等策略进行后续验证设计。
技术路线：从比较基因组学到E12.5胚胎LacZ检测
本文解读整理的技术路线示意图
研究整体思路是先从基因组序列差异出发，再通过体内转基因小鼠增强子检测回答这些序列差异是否具有真实的发育调控功能。计算分析负责锁定候选调控区域，E12.5转基因小鼠胚胎中的LacZ表达则用于判断这些序列是否能在真实发育环境中驱动组织特异性表达。
研究结果一：比较基因组学筛选出166个蝙蝠加速区域
研究人员首先整合发育期小鼠肢体H3K27ac与p300 ChIP-seq数据，并与脊椎动物保守序列重叠，从而筛选出具有潜在肢体增强子特征的候选区域。随后，研究团队基于四种蝙蝠基因组进行PhyloP分析，最终获得166个蝙蝠加速区域。
这些BAR既保留了“在多数脊椎动物中保守”的特征，又在蝙蝠谱系中出现显著快速演化，因此被认为可能与蝙蝠翼发育过程中的基因表达变化有关。
识别蝙蝠加速区域的计算流程（图源：Booker et al., PLOS Genetics, 2016）
在候选区域中，BAR2邻近Twist2，BAR4邻近Spry1，BAR61与Shh的ZRS调控区重叠，BAR97邻近Spg20，BAR116位于HoxD基因簇端粒侧。这些基因或区域均与肢体发育、骨骼形成或肢体形态异常相关，因此具备进一步进行体内功能验证的生物学合理性。
候选BAR在小鼠基因组中的位置及邻近发育相关基因（图源：Booker et al., PLOS Genetics, 2016）
研究结果二：5个蝙蝠BAR在小鼠胚胎中具有肢体增强子活性
为验证候选序列是否真正具有体内增强子功能，研究人员将5个M. lucifugus BAR序列克隆至Hsp68-LacZ报告载体，并注射入单细胞小鼠胚胎。原文方法学说明，所有转基因小鼠均由Cyagen Biosciences按标准流程生成。胚胎在E12.5阶段收集并进行LacZ染色，以观察候选序列能否驱动特定组织中的报告基因表达。
结果显示，5个蝙蝠BAR均在发育期小鼠肢体中表现出增强子活性：BAR2在肢体远端自足区出现信号，BAR4主要在近端肢体区域表达，BAR61在前后肢远端后半部分表达，BAR97呈较弱而弥散的后部自足区信号，BAR116则在前肢近端和远端区域表现出强烈表达，并在后肢近端区域表现出较弱活性。
M. lucifugus BAR在发育期小鼠肢体中的增强子活性（图源：Booker et al., PLOS Genetics, 2016）
研究结果三：蝙蝠与小鼠同源序列呈现不同增强子表达模式
为了判断这些调控元件的功能差异是否与蝙蝠谱系的序列变化有关，研究人员进一步比较了部分BAR的蝙蝠序列和小鼠同源序列。结果显示，BAR61/ZRS在蝙蝠和小鼠序列中均保持类似的后部自足区表达，提示该Shh相关调控元件功能较为保守；而BAR4、BAR97、BAR116则表现出明显差异，尤其是小鼠BAR116序列未能驱动与蝙蝠BAR116相似的肢体增强子活性。
这一结果提示，蝙蝠谱系中发生的非编码序列变化并非只是“序列差异”，而可能已经转化为真实的体内表达调控差异。对于非编码突变、增强子功能以及物种特异性形态形成研究而言，小鼠体内增强子检测的价值正在于此：它可以将候选调控序列与发育过程中的组织表达模式直接连接起来。
蝙蝠BAR与小鼠同源序列在前肢和后肢中的表达差异（图源：Booker et al., PLOS Genetics, 2016）
研究结果四：BAR116核心序列是驱动肢体表达的关键
在所有候选区域中，BAR116最受关注。它位于HoxD基因簇附近，而HoxD基因簇是调控脊椎动物肢体发育和骨骼模式形成的关键区域。为进一步确认BAR116的功能来源，研究人员构建了蝙蝠-小鼠BAR116嵌合序列：将蝙蝠BAR116核心PhastCons序列替换为小鼠同源序列，同时保留蝙蝠侧翼序列。
嵌合构建体在E12.5转基因小鼠胚胎中的肢体表达不稳定，不能像完整蝙蝠BAR116序列那样稳定驱动肢体增强子活性。这说明决定肢体表达的关键并不是侧翼区域，而是M. lucifugus BAR116核心序列本身。换言之，BAR116内的蝙蝠特异性序列变化可能是其获得或维持前肢增强子功能的关键。
蝙蝠-小鼠BAR116嵌合序列导致肢体组织特异性减弱（图源：Booker et al., PLOS Genetics, 2016）
研究结果五：HoxD表达支持BAR116参与蝙蝠前肢发育调控
为进一步判断BAR116与HoxD基因簇之间的功能关联，研究人员检测了蝙蝠和小鼠胚胎中Hoxd10、Hoxd11、Hoxd12、Hoxd13的表达模式。结果显示，在与E12.5小鼠相当的蝙蝠发育阶段，Hoxd10和Hoxd11在前肢远端区域出现明显表达，而在后肢远端区域表达降低。这一“前肢强、后肢弱”的表达趋势与蝙蝠BAR116在转基因小鼠胚胎中的增强子活性方向相似。
虽然作者也指出，Hoxd10和Hoxd11的表达并不能完全复现BAR116的所有空间表达细节，但两者在前肢/后肢差异上的一致性支持了一个重要推断：BAR116可能是调控蝙蝠HoxD基因簇、并参与蝙蝠前肢特殊形态形成的候选顺式调控元件。
蝙蝠与小鼠胚胎中Hoxd10-Hoxd13表达模式比较（图源：Booker et al., PLOS Genetics, 2016）
研究价值：非编码调控元件功能验证的启发
综上所述，本研究以蝙蝠翼发育为切入点，建立了从多物种比较基因组学筛选到转基因小鼠体内功能验证的完整研究路径。研究人员发现166个蝙蝠加速区域，并证明其中5个M. lucifugus BAR均可在发育期小鼠肢体中发挥增强子功能；其中BAR116位于HoxD基因簇附近，蝙蝠序列具有强烈前肢增强子活性，而小鼠同源序列及蝙蝠-小鼠嵌合序列无法稳定驱动相同表达，提示BAR116本身可能是调控蝙蝠前肢发育的重要候选序列。
从论文整体来看，这项研究展示了一条清晰的非编码调控元件研究路径：先通过比较基因组学发现候选区域，再借助转基因小鼠胚胎模型进行体内增强子活性验证，最后结合HoxD基因簇表达模式提出可能的发育调控解释。相比单纯的序列预测，体内表达结果使BAR116等候选区域的功能意义更加直观。
这类研究思路对发育生物学、演化发育学、肢体形态形成、Hox基因簇调控以及疾病相关非编码变异研究都有参考价值。尤其是在非编码区域功能解释越来越受到关注的背景下，能够在胚胎体内观察候选序列的组织特异性表达，是连接“基因组信号”与“生物学功能”的关键一步。
引用来源与文献原文链接
1. Booker BM, Friedrich T, Mason MK, VanderMeer JE, Zhao J, Eckalbar WL, Logan M, Illing N, Pollard KS, Ahituv N. Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus. PLOS Genetics. 2016;12(3):e1005738. DOI: 10.1371/journal.pgen.1005738. PMID: 27019019; PMCID: PMC4809552.
2. 文献原文链接：NCBI PMC 原文页面。
3. 文中模型来源信息：论文“Mouse transgenic enhancer assays”方法部分写明，PCR产物和合成序列被克隆至Hsp68-LacZ载体并完成序列验证，所有转基因小鼠均由Cyagen Biosciences按标准流程生成，并在E12.5阶段收集后进行LacZ染色。