Page 49 - 《水产学报》2025年第6期

金鲟，等 水产学报, 2025, 49(6): 069604

群体 D G Z 群体 D G Z

population population
0.2 0.2

0.1 0.1
PC 2 (1.159 26%) −0.1 0 PC 3 (1.063 85%) −0.1 0

−0.2 −0.1 0 0.1 0.2 −0.2 −0.1 0 0.1 0.2
PC 1 (1.167 89%) PC 1 (1.167 89%)
(a) (b)
图 2 200 尾黄姑鱼关于耐低温时长的群体结构图
(a) PC1 和 PC2 的 PCA 分析，(b) PC1 和 PC3 的 PCA 分析。
Fig. 2 Population structure diagram of 200 N. albiflora based on low-temperature tolerance duration
(a) PCA analysis of PC1 and PC2; (b) PCA analysis of PC1 and PC3.
P-value 的一致性较高 (图 3)，表明 GWAS使用 扫描，寻找候选基因并对其进行功能注释，结
的统计模型能够恰当地对表型数据进行拟合。 果共扫描到 13 个潜在相关的基因，通过功能注
在耐低温全基因组关联分析的曼哈顿图中 释发现主要与免疫、 防御修复、 压力应激和线
筛查出 9 个显著 SNP 位点，其中最显著的位点 粒体能量代谢等过程相关 (图 4，表 4)。

在 18 号染色体上 (P=4.61×10 )。另外，在 5、6、
−8
−6
−7
15、22、23 号染色体上 (P=1.59×10 ~ 1.98×10 ) 3 讨论
均筛查到 1 个显著位点，而 14 号染色体上筛查
研究黄姑鱼应对低温胁迫的分子机制，可
到 3 个位点 (P=9.51×10 ~1.73×10 )(图 3，表 3)。
−6
−8

为黄姑鱼养殖过程中应对极端低温带来的挑战
2.3 候选位点的注释及功能预测 提供理论参考。在本研究中，通过对 200 尾养
对获得的显著位点上下游 50 kb 区域进行 殖黄姑鱼进行 5 °C 的低温胁迫，基于相关性分

8
8 6 6
负 P 值常数对数 −lg(P) 4 2 负观测 P 值常数对数 observed −lg(P) 4

0 2
1 2 3 4 5 6 7 8 9 10 1112 131415161718192021222324 0 1 2 3 4 5 6 7
染色体 负期望 P 值常数对数
chromosone expected −lg(P)
(a) (b)
图 3 黄姑鱼耐低温时长的全基因组关联分析结果
(a) 曼哈顿图，黑色的实线和虚线分别表示−lg(1/n)=6.48 和−lg(10/n)=5.48 时的阈值线，蓝色和橘色的点代表 SNP 位点；(b) Q-Q 图。
Fig. 3 Results of genome-wide association analysis of low temperature tolerance time of N. albiflora
(a) Manhattan plot, where the solid and dashed black lines represent threshold lines at −lg(1/n) = 6.48 and −lg(10/n) = 5.48, respectively, and the blue and
orange dots represent SNP sites; (b) Q-Q plot.

中国水产学会主办 sponsored by China Society of Fisheries https://www.china-fishery.cn
5