《Poultry Science》：Integrating GWAS and transcriptomics to identify candidate genes associated with the crest trait in Yongshun crested duck
编辑推荐：
永顺凤头鸭（YCD）是一种新发现的观赏鸭品种，其特征是头部有一簇球形羽毛。然而，YCD冠羽性状背后的遗传机制尚不清楚。在这项研究中，研究人员整合了全基因组重测序和转录组分析，以确定与冠羽性状相关的候选基因和调控通路。共有175只YCD个体（76只有冠羽，99只
永顺凤头鸭（YCD）是一种新发现的观赏鸭品种，其特征是头部有一簇球形羽毛。然而，YCD冠羽性状背后的遗传机制尚不清楚。在这项研究中，研究人员整合了全基因组重测序和转录组分析，以确定与冠羽性状相关的候选基因和调控通路。共有175只YCD个体（76只有冠羽，99只无冠羽）接受了全基因组关联分析（GWAS）。群体遗传学分析显示有足够的遗传重组且无明显的群体分层，支持该群体适用于关联分析。GWAS在1号染色体上鉴定出显著的关联信号，并突出显示了两个显著的关联峰。关联区域内的候选基因包括MFAP5、APOBEC1、SLC2A3、NANOG和A2ML1，这些基因主要与细胞外基质组织、代谢调节和发育过程相关。为了进一步表征与冠羽形成相关的转录变化，研究人员对5只有冠羽和5只无冠羽雏鸭的头皮组织进行了RNA测序（RNA-seq）。共鉴定出660个差异表达基因（DEGs），包括多个HOX家族基因（HOXA7、HOXB4、HOXB6、HOXA5、HOXB5、HOXC5、HOXC8和HOXC10）和脂质代谢相关基因（PNPLA2、PPARG和CD36）。功能富集分析显示，这些DEGs主要涉及细胞外基质组织、细胞黏附、Wnt、TGF-β、MAPK信号通路和脂质代谢。值得注意的是，GWAS候选基因C3AR1在冠羽相关的头皮组织中也差异表达，为其可能参与冠羽性状提供了基因组和转录组的一致证据。总的来说，研究人员的发现表明，YCD的冠羽性状与1号染色体上的遗传变异以及涉及颅部皮肤发育、细胞外基质重塑、HOX介导的发育调控和脂质代谢的独特转录模式相关。这项研究为鸟类冠羽性状的遗传结构提供了新的见解，并为观赏鸭的保护和分子育种提供了有价值的信息。
**论文解读：永顺凤头鸭冠羽性状的遗传基础研究**
**研究背景与目的**
永顺凤头鸭（YCD）是一种新发现的原产于中国湖南省永顺县的观赏鸭品种，其头部具有一簇球形羽毛，下方存在脂肪组织垫，即冠垫。尽管YCD以冠羽为特征，但在保种群中该性状存在分离，同时存在有冠羽和无冠羽个体。冠羽性状在家鸡、岩鸽、鹦鹉和鸭等多种鸟类中均有报道，其遗传机制在不同物种间存在差异：鸡的冠羽与HOXC8、HOXC10基因的异位表达或调控变异相关，岩鸽的冠羽与EphB2基因的错义突变有关，而中国凤头鸭（RCD）的研究则涉及脂肪生成、颅骨发育及TAS2R40基因的突变。然而，YCD冠羽性状的遗传机制尚不清楚，其是否与其他冠羽禽类共享调控机制或具有品种特异性遗传因子未知。为此，研究人员整合全基因组重测序和转录组分析，旨在鉴定与YCD冠羽性状相关的候选基因和调控通路，为理解鸟类冠羽性状的遗传结构提供新见解。
**研究方法概述**
本研究使用来自永顺凤头鸭保种场的175只180日龄个体（76只有冠羽，99只无冠羽）进行全基因组重测序，样本包含雌雄个体。血液样本提取基因组DNA后，采用酚-氯仿法提取，经Covaris破碎、文库构建后在MGISEQ-2000平台测序。测序数据比对至鸭参考基因组ZJU1.0，经GATK进行SNP calling，并过滤获得高质量SNP。群体结构通过主成分分析（PCA）和遗传关系矩阵（G矩阵）评估，连锁不平衡（LD）衰减分析采用PopLDdecay。全基因组关联分析（GWAS）使用GEMMA软件的线性混合模型，以冠羽性状为二元表型，并将PCA协变量和亲缘关系作为固定和随机效应，显著性阈值通过GEC计算为?log10(P)=7.86。同时，对5只有冠羽和5只无冠羽的1日龄雏鸭头皮组织进行RNA-seq，使用DESeq2鉴定差异表达基因（DEGs），并通过clusterProfiler进行GO和KEGG富集分析。
**研究结果**
*胚胎发育过程中的冠羽形态观察*
研究人员每隔两天从胚胎第4天（E4）到孵化后第1天（D1）观察YCD胚胎，发现头部在早期胚胎阶段出现小突起，并逐渐发育为明显的隆起，表明冠羽性状在胚胎发育期间即开始形成。
*测序质量*
经质量控制后，clean reads的Q20值均≥98.92%，平均测序深度为13.23×。共鉴定出21,234,595个SNP，经软过滤后保留11,303,453个SNP，排除性染色体和不定位scaffold后，最终以9,824,522个常染色体SNP用于关联分析。
*群体结构*
PCA结果显示所有YCD个体聚为一个群体，未根据冠羽表型产生分层；G矩阵显示个体间遗传亲缘关系较弱；LD衰减分析显示r2在0 kb处较高，随遗传距离增加迅速下降并趋于稳定。这些结果表明群体内存在充分的遗传重组，适合进行GWAS。
*GWAS*
以?log10(P)=7.86为阈值，在1号染色体上检测到多个超过阈值的SNP，形成两个显著的关联峰。lead SNP位于Chr1:79,862,378，其原始注释基因LOC119715576在最新NCBI注释中已被撤销，因此该位点所在的基因组区域被视为与冠羽性状显著关联的区域。
*基因注释和功能富集分析*
在显著SNP上下游100 kb范围内注释到多个候选基因，包括MFAP5、APOBEC1、SLC2A3、NANOG和A2ML1（其中A2ML1以多个类似基因模型存在）。功能富集分析发现这些基因显著富集于内肽酶抑制剂活性等分子功能，提示其可能通过细胞外基质组织和发育调控参与冠羽形成。
*冠羽相关头皮组织的转录组分析*
RNA-seq数据的PCA分析显示有冠羽和无冠羽两组头皮组织转录谱明显分离，PC1和PC2共解释45.59%方差。共鉴定出660个DEGs（329个上调，331个下调），其中有冠羽组中多个HOX家族基因（HOXA7、HOXB4、HOXB6、HOXA5、HOXB5、HOXC5、HOXC8、HOXC10）以及脂质代谢相关基因（PNPLA2、PPARG、CD36）显著上调。功能富集分析显示DEGs主要富集于细胞外基质、细胞黏附、Wnt、TGF-β、MAPK信号通路及脂质代谢通路。值得注意的是，GWAS候选基因C3AR1在冠羽头皮组织中表达上调，为基因组关联和转录组变化提供了交叉证据。
**讨论与结论**
讨论部分指出，不同鸟类冠羽性状的遗传基础存在谱系特异性，但可能汇聚于发育信号和颅部组织重塑通路。YCD的GWAS候选基因MFAP5、SLC2A3、APOBEC1、NANOG和A2ML1在功能上涉及细胞外基质、代谢和胚胎发育，但有待实验验证。转录组分析中多个HOX基因的上调提示颅部皮肤区域发育程序改变，而脂质代谢相关基因的表达上调与冠垫的脂肪沉积特征相符。GWAS与RNA-seq结果互补，C3AR1为两者提供了直接关联。结论部分翻译如下：本研究鉴定了一个与YCD冠羽性状相关的1号染色体主要基因组区域，并揭示了冠羽相关颅部皮肤的独特转录变化。基因组和转录组证据整合表明，冠羽形成与发育调控、颅部组织重塑和脂质相关过程的协调变化相关，而非单一孤立的生物学通路。特别是多个HOX家族基因的表达改变突出表明区域颅部皮肤发育与冠羽形成区域之间的潜在联系。这些发现拓宽了对鸭冠羽性状遗传和发育基础的理解，并为未来功能验证提供了候选基因组区域和生物学过程。
