热带爪蟾后肢芽中皮质酮与三碘甲状腺原氨酸信号互作的转录组学及骨骼发育调控研究

《Developmental Biology》：Crosstalk between corticosterone and triiodothyronine signaling in Xenopus tropicalis hindlimb buds.

本研究针对甲状腺激素(T3)与糖皮质激素(皮质酮，CORT)在脊椎动物胚胎后发育中的互作机制尚不明确的问题，通过体外培养热带爪蟾后肢芽并联合RNA-Seq技术，系统解析了T3与CORT协同调控骨骼发育关键基因的表达网络。研究发现CORT通过增强T3信号通路中ZBTB16、TXN等靶点的表达，加速后肢软骨向骨组织分化，并揭示该互作可能以免疫应答抑制和基因组稳定性降低为代价。成果为理解环境应激激素调控发育可塑性的分子基础提供了重要依据。

在脊椎动物发育过程中，甲状腺激素（TH）和糖皮质激素（GC）是调控胚胎后过渡阶段的两大关键内分泌信号。甲状腺激素中最活跃的形式三碘甲状腺原氨酸（T₃）主导器官成熟与代谢转换，而皮质酮（CORT）作为两栖类主要的GC，不仅介导应激反应，更将环境信号转化为发育调控指令。尽管二者在哺乳动物围产期和两栖类变态过程中均呈现协同升高的趋势，且功能交错（如共同影响骨骼、免疫和神经系统成熟），但它们的分子互作机制至今仍缺乏系统性解析。这一空白限制了人们对环境压力如何通过激素“重编程”发育轨迹的理解，也阻碍了相关发育源性疾病的机制探索。

为揭示T₃与CORT在组织自主性层面的互作模式，研究团队以热带爪蟾（Xenopus tropicalis）为模型，聚焦其后肢芽这一形态发生活跃且易观察的器官。通过离体培养NF-53至NF-54阶段后肢芽，分别施加T₃（10 nM）、CORT（100 nM）及两者联合处理（T₃+CORT），24小时后采集样本进行RNA-Seq分析。实验同步设立体内处理组（全身暴露）以排除中枢反馈干扰，并通过RT-qPCR验证测序结果可靠性。在表型层面，研究人员对激素脉冲处理后的蝌蚪进行阿利新蓝-阿尔辛红双染色，量化骨骼成熟度，并在变态完成后测量幼蛙后肢长度。

关键实验技术概述

本研究核心方法包括：①热带爪蟾后肢芽离体培养与激素处理；②RNA-Seq转录组测序（Illumina HiSeq 1500平台，50 bp单端测序）及生物信息学分析（DESeq差异表达鉴定、Fuzzy c-mean聚类、KEGG通路富集）；③RT-qPCR验证靶基因表达；④全骨架染色（Alcian Blue/Alizarin Red）评估软骨-骨转化；⑤网络拓扑分析（Cytoscape构建蛋白互作网络，识别枢纽基因与子网络）。

研究结果

1. 后肢芽转录组特征与激素响应模式

RNA-Seq分析显示，T₃单独处理可显著调控3449个基因，而CORT单独作用仅影响29个基因，表明早期后肢芽对CORT的响应高度依赖T₃背景。联合处理时，4199个基因发生差异表达，其中274个基因（214个“增强型”和60个“减弱型”）表现出T₃-CORT互作特征（图3）。增强型基因中74%为T₃诱导上调基因的进一步强化（如ZBTB16、TXN），减弱型则主要表现为T₃抑制效应的逆转。

2. 互作基因富集于骨骼发育与免疫通路

GO与KEGG分析表明，互作基因显著富集于骨骼发育（如“骨化”“软骨代谢”）、免疫炎症反应和细胞周期调控等相关术语（图4）。手动注释进一步识别出16个与骨骼病理（如骨质疏松、关节炎）密切相关的基因，其中12个经RT-qPCR验证确认表达趋势一致（图4C）。

3. 网络分析揭示枢纽基因与功能子网络

通过整合96条KEGG通路构建蛋白互作网络，发现9个枢纽基因（如PIK3R1、CDKN1A、ZBTB16）在互作中起核心作用（表1）。关键子网络包括：

• •
  免疫调控子网络（TLR2-PIK3R1-CD40-TRAF3-TANK）：联合处理下调该通路，提示免疫应答抑制；
• •
  骨骼发育子网络（NCF1-TRAF4-ZBTB16-TXNIP-TXN）：ZBTB16（促骨分化转录因子）表达增强，TXN（抗氧化蛋白）上调且其抑制因子TXNIP下调，共同促进成骨分化（图5，表2）。

4. CORT协同T₃加速骨骼成熟但导致远期发育代价

表型分析显示，T₃+CORT联合处理48小时后，后肢芽出现局部阿尔辛红强染色（图6C），标志骨化提前1.5个发育阶段（相当于NF-57至NF-58过渡）。然而，在NF-58/59阶段短期暴露CORT的个体，变态完成后后肢长度显著缩短（图6G），表明CORT虽促进短期骨分化，但可能以牺牲远期生长为代价。

结论与意义

本研究首次在器官水平系统描绘了T₃与CORT信号在后肢发育中的转录组互作图谱，证实CORT通过增强ZBTB16等成骨关键因子表达，协同T₃推动骨骼提前成熟。然而，这种加速伴随免疫通路抑制与DNA复制相关基因（如MCM2/MCM4）下调，暗示发育可塑性存在“代价”——可能损害基因组稳定性与免疫适应能力。成果为理解环境应激（如水体干涸）如何通过激素互作驱动发育表型提供了分子依据，对探讨哺乳动物围产期GC暴露引发的骨骼、免疫远期效应（如DOADH假说）具有重要启示。论文发表于《Developmental Biology》，为发育内分泌学领域建立了可推广的研究范式。