编辑推荐：
为解决气候变暖对全球家禽生产和粮食安全的威胁，研究人员针对14个地理和气候多样性鸡种，开展了整合基因组选择扫描、多组织表达数据、eQTL、TWAS和跨物种PheWAS分析的研究。他们发现了一个由ATP1A1、PLCB4、RYR2和AKT3组成的调控枢纽，该枢纽通过协调心血管、钙离子和生存信号，形成了一个整合的体温调节轴。这项研究揭示了驱动热适应的集成心脏-神经内分泌钙离子网络，为培育耐热鸡提供了潜在的遗传靶点。
在全球气候变化的背景下，日益升高的温度和频繁的热浪已成为全球家禽生产和粮食安全面临的紧迫问题。热带和亚热带地区尤为脆弱，高温可导致鸡的生长受损、免疫力下降、产蛋和产肉量减少，并显著提高死亡率。鸡作为对热应激高度敏感的物种，因其缺乏汗腺，且现代集约化育种为提高生长和生产性能，反而因其更高的代谢率而加剧了其脆弱性。一个引人注目的现象是，热带地区的本土鸡种能够在40-45°C的高温环境中常规生存，而这在商业品种中会造成严重损失。然而，鸡对高温耐受的遗传机制仍然知之甚少。理解种群耐受高温的能力，对于认识其地理分布边界和对气候相关挑战的易感性至关重要。为了破解这一科学谜题，一项整合了多组学数据和跨物种验证的研究在《Poultry Science》杂志上发表，揭示了驱动鸡热适应的核心分子网络。
研究人员综合运用了多项关键技术。他们收集了来自巴基斯坦、泰国和中国的4个耐热本土鸡种（Aseel, Desi, Pradu Hang Dam, Thai Native）的血液样本，并整合了来自公开数据库的另外10个品种（包括耐热品种和耐寒品种）以及30个红原鸡的全基因组数据，构建了一个包含14个地理气候多样性品种的272个样本的队列。研究采用了商业FPGA硬件加速平台GTX-One进行变异检测。核心分析包括：1） 通过四种互补的群体遗传统计量（FST、π比率、Tajima's D差异、XP-EHH）进行全基因组选择扫描，以识别受选择的区域和基因；2） 整合鸡GTEx数据库的多组织RNA-seq数据，分析候选基因的表达谱和共表达网络；3） 利用鸡GTEx的eQTL（表达数量性状位点）和TWAS（转录组范围关联研究）数据，验证候选基因的组织特异性调控和基因-性状关联；4） 进行跨物种的PheWAS（表型全基因组关联研究），利用人类GWAS（全基因组关联研究）数据评估关键候选基因的普适性功能。
群体结构、遗传分化和连锁不平衡模式
主成分分析显示，鸡种按照气候分类（热-干、热-湿、寒冷和红原鸡）出现了明显的遗传分化。热区品种和寒冷区品种在PC1上被清晰分开。遗传分化分析（FST）显示，热区品种间的遗传分化较低，而热区与寒冷区品种间的FST值更高，表明存在与气候适应相关的遗传分化。连锁不平衡衰减分析和近交系数也揭示了不同品种间的遗传多样性差异，为后续的选择信号分析提供了合理的群体背景。
全基因组选择扫描分析（耐热 vs 红原鸡）
比较耐热家鸡与祖先红原鸡，通过四种选择检测方法，在基因组顶部2%的区域中鉴定出大量候选基因。交叉验证发现了9个被所有四种方法共同检测到的基因。值得注意的是，研究中还识别出了经典的驯化基因，如甲状腺刺激激素受体（TSHR）、脑源性神经营养因子（BDNF）和TBC1域家族成员1（TBC1D1），这反映了从红原鸡中分离出来后，后续适应不同气候的过程共同塑造了家鸡的遗传和生理多样性。
全基因组选择扫描分析（耐热 vs 耐寒）
为了聚焦于温度适应本身，研究人员直接比较了耐热和耐寒鸡种。分析同样鉴定了数百个候选基因。在至少三种方法中共同检测到的基因有69个，而其中三个基因（AGMO、ARHGAP5和 BMPR2）被所有四种方法识别。这个对比帮助区分了普遍的驯化信号和特异性的温度适应信号。
功能富集分析
对从上述选择扫描中鉴定出的候选基因进行功能富集分析，发现在耐热 vs 耐寒、耐热 vs 红原鸡以及热-干 vs 热-湿三组比较中，有五条关键通路被一致且显著地富集。这五条通路是：心肌细胞中的肾上腺素能信号通路、GnRH（促性腺激素释放激素）信号通路、apelin信号通路、钙离子信号通路和血管平滑肌收缩。这些通路共同指向了心血管调控、神经内分泌和钙离子稳态在热适应中的核心作用。
通路汇聚、组织特异性基因表达和网络连通性揭示了一个协调轴
基于选择信号的强度、在五条核心热响应通路中的富集情况以及组织表达模式，研究最终锁定了25个高置信度的候选基因。其中，PLCB4（磷脂酶C beta 4）与所有五条核心信号通路都有联系，而ATP1A1（Na?/K? ATP酶亚基）、RYR2（兰尼碱受体2）和AKT3（蛋白激酶B）构成了一个关键的调控枢纽。蛋白-蛋白互作网络分析显示，这些基因形成了一个紧密的相互作用模块。组织表达谱分析进一步表明，这些基因在心脏、下丘脑、肝脏等关键体温调节器官中协调表达，形成了一个整合的心脏-神经内分泌轴。例如，钙离子调控相关基因（如RYR2、CACNA2D1）在心脏和骨骼肌中高表达，而信号转导基因（如RAMP2、RAMP3）和应激反应基因（如AKT3、MAPK14）则分别在下丘脑/视网膜和肝脏/气管等组织中高表达。
基因相关矩阵和调控变异
对四个核心枢纽基因（ATP1A1、PLCB4、RYR2、AKT3）的等位基因频率分析发现，在耐热和耐寒群体之间存在大量高度分化的单核苷酸多态性，表明了强烈的选择压力。利用鸡GTEx数据库，研究人员在这些基因附近注释了359个组织特异性的顺式eQTLs，这为这些基因的表达水平受到遗传调控提供了证据，暗示选择可能直接作用于调控变异。
TWAS和PheWAS关联揭示核心基因是热适应的关键调控因子
整合TWAS数据发现，ATP1A1表现出最广泛的性状关联谱，在13个组织和14个性状类别中具有145个显著关联，其中最强的是与肝脏中蛋壳调节的关联。这直接将离子稳态与耐热繁殖力联系了起来。跨物种PheWAS分析利用人类GWAS数据进一步验证了这些基因功能的保守性和多效性。例如，RYR2与静息心率显著相关，ATP1A1与氯离子稳态和免疫细胞计数相关，AKT3与身体质量指数相关，PLCB4与血压和身体成分相关。这些结果在人类复杂性状中印证了这些基因在心血管、代谢和离子平衡中的关键作用。
研究结论与讨论
这项研究揭示了一个协调心血管功能、生殖能力和细胞稳态的整合信号网络，该网络驱动了鸡的热适应。独立的选择统计数据与基因表达、组织特异性调控模式相结合，共同将ATP1A1、PLCB4、RYR2和AKT3确定为持续进化压力下的调控枢纽。尤为值得注意的是，研究发现心血管性能与产蛋量之间存在共享的遗传架构。自然选择似乎倾向于通过钙离子稳态、心脏性能和蛋壳钙化等调控通路，在热应激期间维持生殖功能。这为理解热带本土鸡种如何在极端高温下保持繁殖力提供了关键的遗传学解释。这项研究的意义在于，它通过结合群体基因组学、组织分辨率表达、调控变异和跨物种表型验证，证明了生物对长期环境压力的适应是通过重新连接现有的信号通路（如钙信号、肾上腺素能信号和GnRH信号）来实现的，而非创造全新的机制。这种方法不仅深化了我们对复杂性状适应机制的理解，更重要的是，它为应对气候变化的挑战提供了具体、可操作的基因组靶点，可用于未来的功能研究和耐热家禽育种计划。
