大头鲢（Hypophthalmichthys nobilis）全基因组重测序揭示体重与体尺相关头高位点：LG13体重关联获独立群体复制，LG14/LG04仍属发现阶段

《Biology》：Whole-Genome Resequencing Identifies Loci Associated with Body Weight and Body-Size-Related Head Height in Bighead Carp (Hypophthalmichthys nobilis)

　　

生长性状在大头鲢（Hypophthalmichthys nobilis）中具有经济重要性，但其遗传基础仍不清楚。研究人员利用2,891,492个单核苷酸多态性（SNP），在298尾鱼中开展全基因组关联分析（GWAS），采用控制群体结构和基因组亲缘关系的线性混合模型。五个变异达到各性状全基因组阈值。最强的体重关联位于连锁群13（LG13）的13:40990738:C:T（p = 5.09 × 10^?10；效应 = 5.50 ± 0.87 g）。两个LG04变异与体重相关，而两个相邻LG14变异与头高相关（lead p = 8.34 × 10^?9；效应 = 2.59 ± 0.44 mm）。然而，体长校正后头部尺寸关联显著减弱，表明未确立头部特异性遗传效应。区域连锁不平衡（LD）和条件分析将紧凑的LG13和LG14信号与复杂LG04模式区分开，补充敏感性分析支持其稳健性。Sanger测序确认所有118个受检LG13基因型。LG13体重关联在268尾鱼的独立群体中重现，并在体长校正后仍显著。基因型分层逆转录定量聚合酶链反应（RT-qPCR）检测到效应等位基因携带者中mchr2-like gene13.20、syngr1-like gene32.45和cacng2-like gene32.46表达较高，经多重检验校正。这些发现为LG13体重关联提供独立支持，而LG14和LG04仍为发现阶段候选区域，在考虑潜在育种应用前需要独立复制和功能研究。

研究背景与问题：大头鲢（Hypophthalmichthys nobilis）是重要淡水养殖鱼类，体重和头相关性状与生产及产品评价相关，头部在中国烹饪中具有可食价值。既往研究通过高密度连锁图、2b-RAD全基因组关联分析（GWAS）等鉴定生长和头尺寸相关数量性状位点（QTL）及候选基因，提出tp53bp2、fndc5b等候选；但局部连锁不平衡（LD）结构及关联信号的条件依赖性尚未充分表征，绝对头尺寸关联与整体体型关联也需区分。为此，研究人员利用298尾同塘幼鱼的全基因组重测序数据，围绕体重、体长、头长、头高等性状开展关联分析、区域LD与条件分析、体长校正、独立群体验证及候选基因表达分析，以定义并独立评估商业相关尺寸性状的遗传靶点。

关键技术方法：研究人员对298尾同一池塘养殖的大头鲢幼鱼（中国武汉梁子湖实验基地来源）进行全基因组重测序，获得2,891,492个SNP，采用控制主成分（PC1–PC5）和基因组亲缘矩阵的线性混合模型进行GWAS；对LG13、LG14、LG04代表性位点开展局部LD、条件关联、贝叶斯精细定位（SuSiE）及敏感性分析；用Sanger测序复核118尾发现群体样本；在268尾不同亲本、不同池塘的独立群体中验证LG13位点；利用公共RNA-seq和基因型分层RT-qPCR分析候选基因表达，并以体长或体重校正开展头部性状二次分析。

结果部分：
3.1 表型变异与群体结构：298尾个体表型完整，体重均值9.68 ± 5.47 g，变异系数56.47%；头长/体长比均值0.3348 ± 0.0135。绝对尺寸性状强正相关；表型方差解释比例（PVE）估计标准误较大。主成分分析（PCA）PC1–PC5累计解释11.17%，k-means分为C1–C4，C3体重较高，群体固定指数（F_ST）加权均值0.0491。通过表型统计、PCA、聚类和方差分析得出群体结构与表型变异特征。

3.2 全基因组关联分析识别LG13、LG14和LG04主要区域：五个SNP达全基因组显著。LG13 lead 13:40990738:C:T与体重最强相关（p_LRT = 5.09 × 10^?10；效应5.50 ± 0.87 g），并对多个尺寸性状呈提示关联；LG14两个相邻SNP与头高显著相关（lead p_LRT = 8.34 × 10^?9；效应2.59 ± 0.44 mm）；LG04两个SNP与体重显著、一个与头长提示相关；头长/体长比无全基因组显著位点。通过线性混合模型GWAS得出主要关联区域，且五个显著SNP不等同于五个独立信号。

3.3 区域LD和条件分析细化关联结构：LG13 lead位于7,983 bp高LD区间，条件后区域最小p_LRT升至1.08 × 10^?3；LG14高LD区间12,034 bp，条件后最小p_LRT升至1.20 × 10^?3，但LG13条件后LG14头高关联仍保留；LG04三个代表峰LD低，互相条件后正效应保留，联合条件后区域最小p_LRT = 8.85 × 10^?4；SuSiE在LG13获一个纯度合格可信集（2个SNP），LG14获一个10个SNP可信集，LG04可信集纯度低；稀有纯合排除显示LG04-C较敏感。通过LD、条件、精细定位和排除分析得出信号结构与稳健性差异。

3.4 LG13关联的基因型-表型模式与敏感性：LG13 lead的TT、CT、CC体重均值分别为8.97、16.11、17.83 g；C等位基因携带者均属C3；C3内TT个体体重亦高于其他簇TT个体。留一染色体法（LOCO）模型加入簇指示后效应为5.43 ± 0.87 g，C3内效应为5.73 ± 1.17 g；深度/基因型质量掩蔽后效应方向一致但证据减弱；Sanger复核118/118一致。通过基因型分组、LOCO、掩蔽和Sanger测序得出LG13关联较稳健，但遗传背景混杂不能完全排除。

3.5 头尺寸相对体长的二次分析：ln(HL)对ln(BL)斜率0.953，ln(HH)对ln(BL)斜率1.022；体长校正后GWAS无全基因组显著位点，仅见提示信号；加入ln(BL)后五个代表位点头尺寸效应大幅降低，LG13的ln(HL)和ln(HH)校正后p_LRT分别为0.852和0.547，LG14的ln(HH)校正后p_LRT = 0.0572；加体重后亦减弱。通过log–log回归和体长校正混合模型得出头尺寸关联主要反映整体体型。

3.6 区域注释与表达背景优先候选基因：研究人员按距离、蛋白相似性和公共RNA-seq总结十个区域基因；LG13附近gene82.0注释为pIgR-like，gene81.41未解决；LG14附近gene13.20为MCHR2-like，在 hypothalamus–pituitary（HP）中表达较高；LG04-B附近gene32.45为SYNGR1-like、gene32.46为CACNG2-like；High-versus-Low比较未检出显著差异。通过基因组注释和公共表达数据优先候选基因。

3.7 LG13独立关联复制与基因型分层表达分析：独立群体268尾中TT 223、CT 42、CC 3，C等位基因频率8.96%，哈迪-温伯格平衡（HWE）检验p = 0.453；加性模型中每C等位基因体重增加21.52 ± 1.61 g（p = 1.52 × 10^?31，R² = 0.402），ln(BW)效应0.366 ± 0.032；加体长后体重关联仍为+4.76 ± 0.67 g（p = 1.25 × 10^?11）；排除CC后仍为正。RT-qPCR显示LG14 C携带者HP中gene13.20表达高1.95倍（校正p = 0.0362），LG04-B G携带者gene32.45和gene32.46分别高2.57倍和2.14倍（校正p = 0.00482和0.0305），gene82.0肠表达无差异。通过独立群体单标记回归和基因型分层RT-qPCR得出LG13体重关联获复制，候选基因表达关联被优先。

讨论部分总结：4.1 全基因组图谱扩展了大头鲢尺寸相关区域证据；LG13获技术确认和独立群体复制，LG14和LG04仍为发现阶段关联；与Fu et al.、Liu et al.、Zhou et al.、Chen et al.等既往连锁定位和2b-RAD研究互补而非替代；与Luo et al.、Geng et al.、Kong et al.、Ali et al.等跨物种比较提示主效区域与多基因背景并不互斥。4.2 LG13是聚焦的体重关联靶点：C等位基因限于C3，背景混杂存在；LOCO和C3内分析支持关联，Sanger一致，独立群体支持但缺系谱；绝对效应不可泛化，log尺度效应较接近；约8 kb高LD区间为后续单倍型分析起点；lead位于基因间区，重复序列注释和候选基因低表达使机制未明。4.3 LG14头高关联依赖育种表型：绝对头高显著，但体长校正后明显减弱，不能视为头特异位点；MCHR2-like gene13.20为探索候选，RT-qPCR表达差异不能证明调控关系。4.4 LG04和公共转录组：LG04为复杂区域，可信集纯度低，LG04-C对稀有纯合敏感；候选基因涉及代谢和神经功能，RT-qPCR支持gene32.45和gene32.46，但样本小且体重差异不能排除；公共数据未检出显著差异。4.5 育种意义与局限：LG13宜作为候选加入基因组预测或标记辅助决策，而非独立预测因子；LG14和LG04需独立复制；绝对头尺寸与相对头尺寸需区分；发现群体298尾单塘、性别和系谱未知，独立群体仅测LG13且CC仅3尾，RT-qPCR每组5尾，参考基因仅β-actin作初步评估，精细定位依赖模型假设。

研究结论（翻译）：全基因组关联定位在幼年大头鲢中鉴定出LG13和LG04的体重相关区域，以及LG14上与绝对、体型相关头高相关的区域。LG13 lead基因型经118尾发现群体Sanger测序确认，其正体重关联在268尾不同亲本、另一池塘群体中获独立复制，尽管遗传背景和家系相关混杂不能完全排除。相比之下，LG14和LG04关联在标记水平仍未复制，应视为需在独立群体中确认的初步发现。重要的是，LG13与头部尺寸的关联在体长校正后不再获支持（头长p_LRT = 0.852，头高p_LRT = 0.547），而LG14头高关联明显减弱（p_LRT = 0.0572）。这些发现表明观察到的头尺寸信号主要与整体体型相关，未确立头部特异性遗传效应。因此，LG14标记在缺乏独立复制和体尺校正头性状直接评估时，不应被视为适用于头比例选择。区域LD分析、条件检验和贝叶斯精细定位提供了关联模式及其分辨率的补充信息，但未确立因果变异。基因型分层RT-qPCR鉴定出涉及mchr2-like gene13.20、syngr1-like gene32.45和cacng2-like gene32.46的表达关联，而pigr-like gene82.0未显示可检测的肠道表达差异。这些表达结果优先候选基因以供进一步研究，但不构成LG14或LG04关联信号的独立验证，也未证明因果机制。总体而言，LG13代表一个获复制的体重关联标记，值得进一步评估；LG14和LG04代表发现阶段候选区域，在应用于育种前需要独立复制和功能研究。