白耳黄鸡三叉冠性状分析

张文武，邱峰芳，许继国

（1.东华理工大学化学生物与材料科学学院，江西 南昌 330013；2.南昌师范学院生物技术研究所/江西地方鸡种遗传改良重点实验室，江西 南昌 330029）

【研究意义】了解物种内和物种间表型多样性的遗传基础是现代生物学研究的重大挑战之一。家鸡经过数千年的选择性繁殖聚集了丰富的表型变异，是研究和理解表型多样性遗传基础的重要研究材料。【前人研究进展】鸡冠作为易见的装饰性状，无论在为人类提供肉蛋的商业化品种还是作为观赏用品种上，均受到强烈的人工选择［1］。作为鸡的第二性征，鸡冠与性成熟之间的关系已经得到很好的研究。科学家还发现鸡冠与交配决策、社会等级、体温调节、生殖能力和骨量有重要关系［2-7］。研究表明，多因子调控鸡冠不同形态。通过转录组测序技术已鉴定到与鸡冠发育相关的候选基因有BMP2基因和CHADL基因［8-9］。三叉冠是单冠后端两侧分裂出肉质隆起，形态像叉状（图1）。Pearl 等［10］最早提出三叉冠性状，认为三叉冠是纯单冠和纯豆冠的中间态。有学者通过白来航鸡的杂交试验发现两个常染色体上有两对等位基因控制单冠和分叉，但试验群体不大，其结论仍有待进一步验证［11］。有学者通过杂交试验发现三叉冠基因型为Sss1s1、SSs1s1、ssS1S1、ssS1s1和sss1s1，非三叉基因型为SSS1S1、SsS1S1、SSS1s1和SsS1s1最符合交配试验后代的分布比例，但该假设未能很好解析左叉和右叉等冠型的遗传规律，推测三叉冠是多基因遗传控制的性状［12］。目前对其他冠型方面的研究成果很多，有研究者发现细齿冠宁都黄公鸡的AR与VIPR1 蛋白免疫组织化学阳性细胞比率显著低于其他冠型个体，倒冠个体的AR 与VIPR1 蛋白的免疫组织化学阳性细胞比率均显著高于其他冠型个体，并且在鸡冠窦状血管网内皮细胞及生发层细胞中具有较高免疫组织化学阳性反应［13］。研究人员发现玫瑰冠纯合子的公鸡个体生育能力较低［14］。携带豆冠等位基因的个体其肉垂也比较小且龙骨皮肤变厚［15］。由于进化保守区的重复导致SOX5基因的异位表达，而异位表达的SOX5基因又改变了SSH基因的表达，从而导致鸡冠异常发育［16］。Somes［17］认为该基因座位还存在一种控制蝴蝶冠的等位基因，其显隐性位于V 形冠等位基因（DV）和双冠等位基因（DC）之间。基因组重复与鸡胚鸡冠和两种双鸡冠表型的异位异胚层表达有关［18］。

图1 白耳黄鸡三叉冠性状Fig.1 Trigeminal comb traits of white-ear yellow chickens

【本研究切入点】随着近几年分子生物学技术手段的进步，玫瑰冠、豆冠和双冠等主要冠型的致因突变已经鉴定且性状形成机制也已被初步阐明，但目前尚无有关三叉冠性状的分子遗传学研究。【拟解决的关键问题】生产中发现，白耳黄鸡群体有约5%的个体为三叉冠。本研究以白耳黄鸡为研究材料，采用杂交实验、受精率实验和全基因组关联分析（Genome-wide Association Analysis,GWAS）对白耳黄鸡三叉冠性状进行分析，以期为三叉冠的进一步研究提供参考。

1 材料与方法

1.1 试验材料

供试白耳黄鸡来自江西省上饶市广丰区白耳黄鸡原种场。供试白耳黄鸡的饲养管理流程按原种场管理；饲料营养方面，育雏阶段（1～30 日龄）饲喂正大慢速品种小鸡料，育成阶段（30～170 日龄）饲喂正大育成料。

1.2 试验方法

1.2.1 正反交试验和受精率试验 以60 只120 日龄的白耳黄鸡进行正反交试验，15 只三叉冠公鸡与15 只单冠母鸡交配，以及15 只单冠公鸡与15只三叉冠母鸡交配。以120 只120 日龄的白耳黄鸡进行受精率试验，30 只单冠公鸡和30 只三叉冠公鸡分别与30 只单冠母鸡交配。

1.2.2 GWAS 分析 对82 只150 日龄的白耳黄公鸡（23 只三叉冠、23 只单冠）于翅静脉采血各5 mL，用EDTA 抗凝处理，提取基因组DNA 用于GWAS 分析。基因分型采用Illumina 平台的鸡京芯一号60 k 芯片。数据质控使用PLINK 1.90 软件。基因型质控过程将删除达到以下标准的个体或位点：分型缺失数＞10%的位点、最小等位基因频率＜ 0.01 的位点和哈代-温伯格检验P值＜1×10-4的位点。

在GWAS 分析中，试验结果的假阳性部分来自群体分层［19］。主成分分析（PCA）方法能用于遗传学的聚类分析，根据遗传差异可将个体分为不同亚组。本研究采用PCA 方法进行群体分层分析，使用GCTA 软件对有效SNP 进行PCA，获得基因组数据特征值和特征向量，并用R 程序作PCA 图。GWAS 分析采用混合线性模型GEMMA。采用Bonferroni 多重检验方法确定关联显性阈值，阈值为1/有效SNP 位点数量。

在NCBI 网站上下载对应物种的参考基因组信息（https://www.ncbi.nlm.nih.gov/genome/111?genome_assembly_id=1543395），使用ANNOVAR 软件将显著SNP 注释到其对应的基因上，并结合文献注释候选基因的功能。

2 结果与分析

2.1 白耳黄鸡三叉冠表型遗传方式

用三叉冠白耳黄鸡个体与单冠白耳黄鸡个体进行正反交试验，结果（表1）显示后代全为单冠。在遗传杂交试验中通常采用正交和反交来判断某性状的基因是位于常染色体上还是性染色体上。根据后代全为单冠初步推测，白耳黄鸡三叉冠为常染色体遗传。

表1 白耳黄鸡正反交试验Table 1 Orthogonal and reverse cross tests of white-ear yellow chickens

2.2 单冠和三叉冠白耳黄鸡的种蛋受精率比较

比较30 只单冠公鸡和30 只三叉冠公鸡与配单冠母鸡的受精能力，结果显示，与配单冠群体的种蛋受精率为94.9%（258/272），与配三叉冠群体的种蛋受精率为92.7%（241/260），差异不显著。

2.3 白耳黄鸡群体有效SNP 数和群体结构分析

GWAS 分析首先要对基因型及表型进行一定标准的质控，以排除异常基因型和异常表型对分析结果产生的影响。经质控白耳黄鸡群体得到55 023 个有效SNP。通过对标记矩阵进行PCA 分析，提取解释方差最大的前几位PCA，绘制对应散点图可反映群体结构。从图2 可以看出，白耳黄鸡群体的点相对集中，表明白耳黄鸡群体内不存在明显的群体分层，可以进行后续分析。

图2 白耳黄鸡群体结构分析Fig.2 Population structure analysis of white-ear yellow chickens

2.4 白耳黄鸡三叉冠显著关联SNP 位点分析

图3 和图4 分别是对白耳黄鸡三叉冠GWAS分析所得结果的QQ 图和曼哈顿图。QQ 图中X轴和Y 轴分别为期望的P值和实际的P值，两者一致时，由随机漂变导致；两者不一致时，存在真实变异，对应的点与表型相关。在曼哈顿图中每个基因组SNP 位点沿染色体序列排列，Y 轴的P值可理解为每个SNP 与表型的关联程度。结果（表2）显示，与白耳黄鸡三叉冠显著关联SNP位点共有10 个，其中有1 个SNP 位于第1 号染色体，2 个SNP 位于第2 号染色体，2 个SNP 位于第3 号染色体，1 个SNP 位于第4 号染色体，1 个SNP 位于第5 号染色体，1 个SNP 位于第12 号染色体，2 个SNP 位于第14 号染色体。这10 个SNP 分别邻近或座落于SPRY2、NDFIP2、ITGA9、EMC2、TMEM17、EHBP1、SCAF8、TIAM2、HTT、STON2、PRKCD、RBBP6、TNRC6A、MAPK8IP3基因上。

图3 白耳黄鸡三叉冠关联分析QQ 图Fig.3 QQ plot of association analysis of white-ear yellow chicken with trigeminal comb

图4 白耳黄鸡三叉冠关联分析曼哈顿图Fig.4 Manhattan diagram of association analysis of white-ear yellow chicken with trigeminal comb

表2 白耳黄鸡三叉冠显著关联SNP 位点Table 2 SNP loci significantly associated with trigeminal comb of white-ear yellow chicken

3 讨论

3.1 白耳黄鸡三叉冠遗传方式

在遗传杂交试验中通常采用正交和反交来判断某性状的基因位于常染色体上还是性染色体上。有学者通过白色来航鸡的杂交试验发现两个常染色体上有两对等位基因控制单冠和分叉冠［11］。有学者通过杂交试验发现三叉冠基因型是常染色体多基因遗传控制的性状［12］。本研究发现白耳黄鸡三叉冠呈现常染色体遗传，这与前人研究结果相一致。

3.2 白耳黄鸡三叉冠与受精率的关系

有研究报道，玫瑰冠与豆冠、单冠的精子活率和受精率差异不显著，但玫瑰冠与单冠在精子密度和受精蛋孵化率上显著高于豆冠［20］。有学者研究玫瑰冠个体的精液质量与不同基因型个体的生育能力和生育持续时间，发现精子活力缺陷是导致低生育能力以及生育时间短的原因［21］。但还未有研究三叉冠与受精率关系的报道，我们发现与配三叉冠公鸡的母鸡群体的种蛋受精率低于与配单冠公鸡个体，但差异不显著，具体原因还需要进一步研究。

3.3 白耳黄鸡三叉冠的GWAS 分析

在鸡生产中，GWAS 被用来鉴定与生产相关性状的候选基因。本研究对白耳黄鸡三叉冠进行GWAS 分析，检测到10 个显著SNP，分别邻近或座落于SPRY2、NDFIP2、ITGA9、EMC2、TMEM17、EHBP1、SCAF8、TIAM2、HTT、STON2、PRKCD、RBBP6、TNRC6A和MAPK8IP3基因上。但根据前人研究结果［22-36］，这些基因并非与白耳黄鸡三叉冠形成直接相关的候选基因。

4 结论

本研究对白耳黄鸡三叉冠性状进行分析，结果显示，白耳黄鸡三叉冠呈现常染色体遗传；GWAS 分析发现10 个与白耳黄鸡三叉冠显著关联的SNP，分别邻近或座落于SPRY2、NDFIP2、ITGA9、EMC2、TMEM17、EHBP1、SCAF8、TIAM2、HTT、STON2、PRKCD、RBBP6、TNRC6A、MAPK8IP3基因上，但这些基因并非与白耳黄鸡三叉冠形成直接相关的候选基因，这些基因是否参与三叉冠的形成需要进一步证实。