高SNP密度基因与功能通路分析揭示肉鸡饲料转化率的关键调控因子

《Journal of Genetic Engineering and Biotechnology》：High-SNP-density genes and functional pathway analysis reveal key regulators of feed conversion ratio in broiler chickens

　　

饲料转化率（FCR）是肉鸡生产中的关键经济性状，影响生产力与可持续性。研究人员从动物数量性状位点数据库（Animal QTLdb）中检索到2272条与FCR相关的初始高峰单核苷酸多态性（SNP）记录，经过滤与去重后保留513个唯一非冗余SNP，利用Ensembl变异效应预测器（VEP）进行功能注释，映射至543个基因并分析其基因组分布与功能效应。多数SNP为内含子变异（54.0%），其次为下游（17.7%）、上游（10.2%）与同义变异（6.5%）。错义变异（1.4%）出现于ST3GAL6、NOX4与TGFBR3等生物学相关基因中，2.5%的SNP位于调控区。生物型分类显示93.7%的变异发生于蛋白质编码区，少部分位于长链非编码RNA（lncRNA）、增强子与微RNA（miRNA）区域。SNP聚类揭示潜在基因组热点，ST3GAL6与DGKZ呈现异常高密度，提示其在FCR相关调控通路中的作用。通过g:Profiler的功能富集分析识别出与转录因子结合基序（如PU.1、NF-E4、MAFA）的显著关联，以及催化活性（GO:0003824）与磷酸酶相关复合体的富集。上述发现突出了ST3GAL6、DGKZ等潜在遗传靶标及包含催化活性的调控网络，可能显著影响FCR。

研究背景与立项依据：饲料转化率（Feed Conversion Ratio, FCR）定义为特定时期内采食量体重增量的比值，数值越低代表效率越高。该性状在肉鸡养殖中占据核心经济地位，饲料成本通常占生产总支出70%以上，且FCR遗传力中等偏高（0.36–0.45），具备通过遗传选择改良的空间。然而FCR属典型多基因复杂性状，受大量微效基因共同作用，传统选择手段面临精度瓶颈。尽管全基因组关联研究（Genome-Wide Association Study, GWAS）已定位众多与FCR显著关联的SNP，但统计关联如何转化为生物学功能仍待解析。为此，研究人员整合Animal QTLdb高峰SNP、VEP功能注释与g:Profiler富集分析，系统挖掘肉鸡FCR的基因组架构与调控网络，相关成果发表于《Journal of Genetic Engineering and Biotechnology》。

主要技术方法概述：研究人员从Animal QTLdb获取Marshall肉鸡FCR相关2272条高峰SNP记录，保留在至少两项独立研究中以p<1×10^?5复现的513个唯一非冗余SNP；使用基于Gallus gallus GRCg7b组装的VEP预测分子后果与生物型；将注释所得543个基因提交g:Profiler，联合Gene Ontology（GO）、KEGG、Reactome与TRANSFAC进行功能与转录因子基序富集，以Benjamini–Hochberg FDR<0.05为显著阈值，低置信度修饰型变异不进入富集。

研究结果

1. 3.
   Results：513个SNP映射至543基因。分子后果分布上，内含子变异占54.0%，下游17.7%，上游10.2%，同义6.5%，错义1.4%（落于ST3GAL6、NOX4、BLM、RPL22、MVP、TXLNB、WNT11B、MBOAT4、ANKRD1、TGFBR3），3′UTR占6.9%，5′UTR 0.57%，调控区变异2.5%，剪接多嘧啶区变异1.2%（落于SLC38A2、HMGA2、HAL等），48个SNP具多重后果。生物型上93.7%在蛋白编码区，3.25% lncRNA，1.36%增强子，0.66%开放染色质，0.44% miRNA，0.13%启动子，余为snoRNA等。基因内SNP重叠分析显示512基因含1–10个SNP，23基因含11–20个；CDK2AP1、COL6A2、KIAA1522、MPDZ、STAU2、VPS4B含21–30个；DGKZ含40个、ST3GAL6含50个，呈极端高密度热点。富集分析中PU.1基序关联459基因，NF-E4关联348，E2F-4关联303，E2F关联248，SOX关联54，MAFA高置信基序TF:M01709_1关联293基因；分子功能层面催化活性（GO:0003824，FDR=0.0005）涉及147基因最显著，丝氨酸/苏氨酸磷酸酶复合体、磷酸酶复合体各含7基因（ITPR1、PPP1CB、PPP1CC、PPP1R12A、PPP2CA、PTPA、WDR82），PTW/PP1与PP1型复合体各含4基因。
2. 4.
   Discussion：内含子主导（约54%）印证非编码变异经剪接或调控元件互作影响表达；ST3GAL6、NOX4、WNT11B错义变异提示少量编码区改变亦参与；SLC38A2、HMGA2等多嘧啶区变异或扰剪接。93.7% SNP落蛋白编码区与既往肉鸡GWAS一致。lncRNA、增强子、miRNA、snoRNA等非零散存在，说明转录与转录后多层调控并存。ST3GAL6（50 SNP）与DGKZ（40 SNP）高密度符合育种选择热点特征，可作标记辅助或基因组选择候选。转录因子PU.1、NF-E4、E2F-4、MAFA、SOX跨数百基因出现，指向代谢与胰岛发育相关协同转录控制；催化活性最强富集呼应空肠十二指肠转录组与胸肌蛋白组中能量代谢、氧化还原与氧化磷酸化差异；磷酸酶复合体（PP1/PP2A家族）富集显示可逆去磷酸化开关联结糖元、糖酵解与线粒体功能。
3. 5.
   Conclusions：研究人员确认513 SNP定位于543基因，蛋白编码区为主，ST3GAL6与DGKZ为SNP密度热点；富集凸显催化活性与磷酸酶复合体所中枢的代谢/信号过程。上述基因与通路可纳入基因组预测模型或靶向SNP面板以提升肉鸡FCR选择准确度，未来需经转录组、蛋白组与基因组编辑验证因果。

（作者：Heydar Ghiasi、Navid Ghavi Hossein-Zadeh；伦理声明None；基金None）