《Poultry Science》：Integrated multi-tissue transcriptomic and network analyses reveal systemic crosstalk governing the onset of lay in Wuhua Yellow chickens
编辑推荐：
家禽开产启动是一个需要神经内分泌—代谢紧密协同的系统性转变过程。尽管地方品种保留了宝贵的遗传多样性和抗逆性，但其多组织调控架构仍鲜有阐明。本研究整合了五华黄鸡（Wuhua Yellow chickens）产前（16周龄）和产蛋高峰期（32周龄）四个关键组织（下
家禽开产启动是一个需要神经内分泌—代谢紧密协同的系统性转变过程。尽管地方品种保留了宝贵的遗传多样性和抗逆性，但其多组织调控架构仍鲜有阐明。本研究整合了五华黄鸡（Wuhua Yellow chickens）产前（16周龄）和产蛋高峰期（32周龄）四个关键组织（下丘脑、垂体、卵巢和肝脏）的转录组数据，结合加权基因共表达网络分析（WGCNA）、配体—受体推断和实时荧光定量PCR（qPCR）验证。结果表明，开产启动由协调的下丘脑—垂体—卵巢（HPO）—肝脏表达程序驱动，而非孤立的器官响应。系统性催乳素（PRL）诱导与肝脏PRLR上调（log?FC=4.23）相偶联，同时伴随肝脏卵黄前体基因的大规模诱导：VTG1（log?FC=18.28）、VTG2（log?FC=6.60）和APOB（log?FC=3.31）。垂体脂蛋白受体（VLDLR和LRP2）的相互下调提示可能存在与外周代谢状态一致的反馈相关通路。WGCNA鉴定了与轴系功能分工相关的组织特异性枢纽基因，包括下丘脑（ADGRL3、NRXN1）、垂体（EYA1）、卵巢（MAD2L1、PLK1）和肝脏（HGFAC、PAH）模块。肝脏与中枢/性腺模块之间显著的负相关（r=?0.48至?0.68，P<0.05）凸显了将资源重新分配至卵黄合成的系统性权衡。靶向qPCR证实了核心信号和卵黄前体基因（PRL、VTG1、VTG2、APOB和SPP1）的强烈阶段依赖性变化，以及肝脏PAH在产蛋高峰期的显著下降，而其他代表性枢纽基因（NRXN1、EYA1和PLK1）与RNA-seq趋势一致。总体而言，该多组织框架为理解地方鸡生殖激活的阶段依赖性转录协调提供了基于转录组的分析框架，并提名了表达水平候选基因，供进一步遗传和功能验证，以期在保留本地遗传多样性的同时改善产蛋性能。
论文解读
一、研究背景与科学问题
产蛋性能是家禽的核心经济性状，受神经内分泌、卵巢和代谢调控的复杂协同作用影响。商业蛋鸡经数十年高强度选育，年产蛋量通常超过300枚，而多数中国地方品种年产蛋量往往不足150枚，但保留了宝贵的遗传多样性、抗逆性和优良肉品质。理解未选择群体中开产启动如何在多组织水平上系统组织，对于揭示核心生殖生物学机制和指导兼用型地方鸡的可持续遗传改良均至关重要。
经典的下丘脑—垂体—卵巢（HPO）轴通过促性腺激素释放激素（GnRH）、促性腺激素和卵巢类固醇反馈启动生殖激活。然而在鸟类中，成功的蛋形成还依赖于肝脏——肝脏在雌激素及相关内分泌调控下合成并分泌卵黄前体卵黄蛋白原（VTG）和含载脂蛋白B（APOB）的脂蛋白。催乳素（PRL）及其受体（PRLR）进一步在产蛋转变过程中调节生殖状态和肝脏代谢反应。因此，调控产前到产蛋高峰期转变的功能单元应被概念化为扩展的HPO—肝脏轴，而非孤立的HPO轴。
尽管近年转录组和多组学研究已开始绘制该轴系的组成，但仍存在重要空白：多数研究聚焦单一组织、在单一成年阶段比较高产与低产品系，或主要关注商业品种而非地方遗传资源；因此，真正开产启动阶段的依赖性组织间串扰仍未完全阐明。此外，中枢与外周器官间的配体—受体通讯通常从循环激素测定或候选基因推断，而非从全基因组多组织表达网络解码。将共表达拓扑与细胞间信号图谱及正交qPCR验证相结合的研究在地方家禽中极为有限。
五华黄鸡（WYC）是源自广东的地方兼用型品种，兼具上述挑战与机遇。基因组调查已描述其群体结构、遗传变异和选择信号；同群体GWAS已鉴定产蛋性能候选位点；卵巢m?A和转录组分析已关联表观转录调控与性成熟。然而，产前到产蛋高峰期转变的四组织（下丘脑、垂体、卵巢、肝脏）同步视角，以及连接中枢神经内分泌信号与外周合成输出的组织间通讯路径，仍然缺失。
二、研究内容与主要结论
研究人员对五华黄鸡产前（16周龄）和产蛋高峰期（32周龄）两个阶段的下丘脑、垂体、卵巢和肝脏四组织进行了RNA测序（RNA-seq），整合差异表达分析、全局WGCNA、多阶段配体—受体通讯推断和靶向qPCR验证，旨在：（i）定义开产启动的系统性转录重塑；（ii）鉴定编码HPO—肝脏轴功能分工的组织特异性共表达枢纽；（iii）利用直系同源映射的互作数据库，通过表达推断绘制连接中枢内分泌输出与肝脏卵黄前体合成及卵巢接收的组织间配体—受体轴。
研究得出以下核心结论：开产启动由协调的HPO—肝脏表达程序驱动。系统性PRL诱导与肝脏PRLR上调（log?FC=4.23）相偶联，伴随肝脏卵黄前体基因大规模诱导：VTG1（log?FC=18.28）、VTG2（log?FC=6.60）和APOB（log?FC=3.31）。垂体脂蛋白受体（VLDLR和LRP2）的相互下调提示可能反映外周代谢状态的反馈通路。WGCNA鉴定了组织特异性枢纽基因，包括下丘脑（ADGRL3、NRXN1）、垂体（EYA1）、卵巢（MAD2L1、PLK1）和肝脏（HGFAC、PAH）模块。肝脏与中枢/性腺模块间显著的负相关（r=?0.48至?0.68，P<0.05）凸显了将资源重新分配至卵黄合成的系统性权衡。qPCR验证了核心信号和卵黄前体基因（PRL、VTG1、VTG2、APOB和SPP1）的强烈阶段依赖性变化及肝脏PAH在产蛋高峰期的显著下降。
三、主要关键技术方法
研究人员采用以下关键技术方法：对五华黄鸡产前（16周龄）和产蛋高峰期（32周龄）两阶段各6只母鸡（每阶段n=6）进行采样，其中每组3只用于RNA-seq（每组织—阶段n=3，共24个文库），全部6只用于qPCR验证。RNA-seq基于Illumina NovaSeq X Plus平台进行双端测序（150-bp读长），比对至鸡参考基因组GRCg6a，以TPM定量基因表达。差异表达分析使用DESeq2，采用分层阈值策略（外周组织FDR<0.05且|log?FC|>1；神经内分泌组织未校正P<0.05）。功能富集分析采用KOBAS（GO/KEGG）和GSEApy（MSigDB）。WGCNA在24个样本上以软阈值β=8构建共表达网络，模块划分采用动态树切割，枢纽基因筛选标准为MM>0.9且GS>0.9。配体—受体对从CellChatDB获取并映射至鸡直系同源基因，构建多阶段Sankey通讯网络。qPCR以GAPDH为内参，采用2?ΔΔCt法计算相对表达量，以Student's t检验评估显著性。
四、研究结果
1. 开产启动的表型和系统转录组概览
产前五华黄鸡卵巢呈未发育形态，含小而未分化的卵泡；产蛋高峰期卵巢呈现完全发育的等级卵泡（F1—F4）和小黄卵泡（SYF）。主成分分析（PCA）显示四个组织在产前和产蛋高峰期组间均呈清晰分离，表明整个神经内分泌—生殖—代谢轴存在深刻的阶段特异性转录重编程。分层阈值策略共鉴定下丘脑903个、垂体3419个、卵巢414个、肝脏506个差异表达基因（DEG）。VTG2是少数在四个组织均检测为差异表达的基因之一，且在肝脏中上调最强（log?FC=6.60）。
2. 生殖轴中枢的功能特化
下丘脑和垂体共享484个DEG，反映协调的神经内分泌重塑。关键中枢网络枢纽（NAV3、KIDINS220、NPY）呈同步表达变化。PRL在下丘脑显著上调（log?FC=7.77）。垂体呈现最广泛的转录组变化（3419个DEG），核心生殖通路通过EGR1等基因激活，ESR1被抑制，与调整的反馈控制一致。结构调控因子如MYL2（log?FC=?6.13）大幅下调，表明垂体转变为适合产蛋高峰的功能成熟状态。
3. 终端器官的功能协同与分工
卵巢中DEG协调卵泡发育和排卵，包括卵黄沉积爆发（VTG2 log?FC=5.93、APOB log?FC=4.89）、类固醇合成重构（CYP17A1和STAR下调）以及基质金属蛋白酶13（MMP13，log?FC=5.20）和多种胶原蛋白的上调，表明广泛的细胞外基质（ECM）重塑。肝脏转录组呈现特化代谢状态：VTG1大幅上调（log?FC=18.28），与VTG2和APOB构成蛋生产的主要组分；脂质代谢重构以维持底物流（上调ABHD5、下调LPL log?FC=?6.35）；激素敏感性增强（PRLR和PAQR7上调），分解代谢基因（IDO2）受抑制。
4. 多组织协同的功能通路与GSEA基础
下丘脑将昼夜节律光周期门控（CLOCK、NPAS2）和能量稳态（NPY、ADIPOR2）与神经内分泌输出整合，以神经活性配体—受体互作通路为核心（NES=1.95，FDR<0.01），进一步激活钙、MAPK和GnRH级联。垂体呈现明确的生理优先级：最大化生殖激素合成并确保细胞存活——凋亡显著抑制（NES=?1.57），核糖体（NES=?2.40）和氧化磷酸化（NES=?1.84）显著下调。卵巢呈现从早期细胞命运决定到快速组织扩张和营养接收的发育转换：ECM重塑（ECM-受体互作NES=2.29）、血管生成和类固醇生成激活，Notch信号（NES=?1.83）显著抑制。肝脏在雌激素刺激下优先脂质处理和分泌（蛋白输出NES=1.77），选择性抑制基础翻译机器（核糖体NES=?2.46）、免疫监视（补体级联NES=?2.27）和稳态营养分解代谢。
5. 全局WGCNA共表达网络构建与模块关系拓扑分析
以软阈值β=8构建全局共表达网络，满足无标度拓扑（R2>0.85），共获得18个模块。模块特征基因（ME）相关热图显示特定模块对呈现显著协同（正）或拮抗（负）共表达动态。阶段相关关联主要集中于较小的高度特化子模块，表明产蛋高峰期转变由组织间功能协同而非单一模块主导性转变决定。
6. 组织特异性枢纽基因与功能调控网络
各组织特征模块内MM与GS强正相关（r≥0.76，P=1×10?2??），表明模块内枢纽基因高度代表组织特异性生理特征。下丘脑Blue模块富集肾上腺素能信号、钙信号和细胞黏附分子通路，核心枢纽为ADGRL3、NRXN1和SYT11。垂体Darkmagenta模块富集内吞和糖鞘脂生物合成通路，核心节点为EYA1、EYA2、DACH1和PAX6。卵巢Bisque4模块富集细胞周期、真核生物核糖体生物发生、卵母细胞减数分裂和孕酮介导的卵母细胞成熟通路，核心调控因子为PLK1、MAD2L1、INCENP和EYA2。肝脏Brown模块富集代谢通路和PPAR信号通路，核心驱动基因为F2、PROC、PLG、HGFAC、PAH和NR1H4（FXR）。
7. 模块间串扰与多阶段配体—受体通讯网络整合
肝脏模块与下丘脑（r=?0.68，P<0.001）、垂体（r=?0.63，P<0.01）和卵巢（r=?0.48，P<0.05）模块呈显著负相关，捕捉了将能量重定向至肝脏卵黄前体合成所需的系统性代谢权衡。多阶段Sankey网络映射了跨组织通讯流：系统性PRL（源自下丘脑、垂体、卵巢）靶向肝脏PRLR，肝脏来源的VTG1、VTG2和APOB反馈至垂体受体（VLDLR和LRP2），局部回路（NPY—NPY5R、VIP—VIPR1、GAL—GALR1、SPP1—CD36）精细调控生殖激活。
8. 关键候选基因表达谱验证
qPCR验证显示与RNA-seq TPM谱高度一致：PRL在垂体产蛋高峰期显著升高，肝脏卵黄前体基因（APOB、VTG1、VTG2）多倍上调，卵巢SPP1显著上调（P<0.05）。肝脏PAH在产蛋高峰期显著下降（P<0.01），NRXN1、EYA1和PLK1的表达轨迹与RNA-seq趋势一致（P>0.05）。
五、讨论与结论
研究讨论指出，开产启动是偶联神经内分泌指令与外周代谢能力的系统性生理转变。核心发现是PRL在HPO—肝脏轴多个节点协调上调，同时伴随肝脏PRLR的显著诱导，支持催乳素信号在开产启动中除经典雌激素作用外积极参与合成程序的模型。卵巢呈现互补的"汇"特征，垂体脂蛋白受体下调可能反映外周代谢状态的反馈感知。模块特征基因负相关在拓扑层面独立印证了系统性代谢权衡。组织特异性共表达枢纽编码了轴系功能分工：下丘脑枢纽参与光周期、能量和神经内分泌输入感知；垂体枢纽与激素包装和释放相关；卵巢枢纽锚定细胞周期和卵母细胞成熟通路；肝脏枢纽锚定代谢通路和PPAR信号，确认肝脏作为卵黄前体输出主要位点的角色。PAH的显著下降支持肝脏模块作为多通路合成工厂的广泛观点。
研究局限性包括：RNA-seq每组织—阶段n=3对细微信号检测有限；卵巢采样阶段差异（产前全卵巢vs高峰期小黄卵泡）可能贡献卵巢特异性DEG；配体—受体关联为表达推断且未经循环激素或受体扰动验证；垂体VTG1/VTG2诱导的生物学相关性需蛋白和细胞类型水平确认；仅比较两个阶段，组织间信号的时间顺序未解决；未采样骨组织，无法评估对髓质骨重塑的远端效应；转录组枢纽与遗传信号互补，需eQTL或多组学整合进行因果推断。
研究结论：五华黄鸡的开产启动呈现为同步的神经内分泌—肝脏—卵巢对话。通过推断实施该对话的细胞间通讯路径和组织内枢纽架构，本研究推进了对地方鸡生殖激活的系统层面理解，并提名了优先表达水平候选基因，供进一步验证以可持续改善产蛋性能。
