高质量组装揭示连城鸭独特羽色表型的PMEL基因

A high-quality assembly revealing the PMEL gene for the unique plumage phenotype in Liancheng ducks.

作者信息Zhen Wang, Zhanbao Guo, Hongfei Liu, Tong Liu, Dapeng Liu, Simeng Yu, Hehe Tang, He Zhang, Qiming Mou, Bo Zhang, Junting Cao, Martine Schroyen, Shuisheng Hou, Zhengkui Zhou
PMID39804725
发布时间2025-01-06
DOI10.1093/gigascience/giae114

实验完整度

高

包含基因组组装、GWAS、转录组、qPCR、免疫荧光、荧光素酶报告基因及Hi-C等多层级实验验证。

主要模型

连城鸭(白色羽毛) 北京鸭(白色羽毛) F2分离群体(连城鸭×北京鸭) 鸭胚成纤维细胞(DEF)和A375黑色素瘤细胞

重点核对

F2群体表型分离比符合孟德尔独立分配定律(3:6:3:4) PMEL基因在黑白羽鸭中的差异表达 Chr33: 5,303,994A>G和5,303,997A>G两个SNP与表型共分离 荧光素酶报告基因验证增强子活性 Hi-C显示候选SNP与PMEL位于同一TAD

摘要

背景 羽色是鸭的一个显著特征,连城鸭以其白色羽毛和黑色喙及蹼足为特征,是此类研究的绝佳对象。然而,学术界对鸭羽色遗传机制的理解仍然有限。为此,本研究利用HiFi reads从头组装了连城鸭基因组(GCA_039998735.1),并由连城鸭和北京鸭生成了F2分离群体。目的是鉴定连城鸭白色羽毛的遗传机制。 结果 本研究从头组装了1.29 Gb的连城鸭基因组,其中contig N50为12.17 Mb,scaffold N50为83.98 Mb。除了MITF基因的上位效应外,全基因组关联分析定位到一个包含PMEL基因的0.8 Mb基因组区域。该基因编码一种色素细胞特异性蛋白,对于黑素体(负责色素沉着的细胞器)内纤维片层的形成至关重要。此外,连锁不平衡分析揭示了2个候选单核苷酸多态性(Chr33: 5,303,994A>G; 5,303,997A>G)可能改变PMEL转录,从而可能影响连城鸭的羽色。 结论 我们的研究组装了连城鸭的高质量基因组,并提供了令人信服的证据表明该品种的白色羽毛特征归因于PMEL基因。总的来说,这些发现为未来研究和育种计划提供了重要的见解和方向,以理解和操纵禽类羽色。 关键词 鸭,基因组组装,羽色,PMEL,黑色素

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
连城鸭白色羽毛的遗传机制是什么?
核心机制
PMEL基因上游两个连锁SNP(Chr33: 5,303,994A>G和5,303,997A>G)通过影响增强子活性,远程调控PMEL基因表达,导致黑色素沉积减少,从而产生白色羽毛。
主要证据
基于高质量连城鸭基因组,通过GWAS、IBD分析、转录组、qPCR、免疫荧光和荧光素酶报告基因实验证实PMEL为Rr基因,且两个SNP增强子活性差异显著。
研究意义
研究为鸭羽色遗传机制提供新见解,为家鸭育种和禽类羽色研究提供分子标记和方向。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

基因组组装与注释

获得连城鸭高质量参考基因组,用于后续变异检测和关联分析。

利用PacBio HiFi测序和Hi-C技术组装连城鸭基因组,并进行基因注释。

2

群体构建与表型鉴定

构建分离群体以研究羽色遗传模式。

用30只北京鸭和120只连城鸭杂交,获得F1和F2群体,记录羽色表型。

3

全基因组重测序与关联分析

定位与羽色相关的基因组区域和候选基因。

对366只鸭进行重测序,进行GWAS和IBD分析,定位到Chr33上的候选区域。

4

候选基因表达分析

筛选候选基因中与羽色相关的差异表达基因。

通过RNA-seq和qPCR分析候选区域基因在黑白羽鸭中的表达,以及免疫荧光检测PMEL蛋白分布。

5

因果突变筛选

鉴定导致羽色差异的具体变异。

基于IBD和FST分析筛选候选SNP,并通过多品种基因型-表型共分离验证。

6

功能验证

验证候选SNP对PMEL基因表达的调控作用。

构建pGL3荧光素酶报告载体,在DEF和A375细胞中检测启动子和增强子活性。

7

染色质构象分析

探讨候选SNP与PMEL基因的空间位置关系。

通过Hi-C分析检测PMEL及其上游区域是否位于同一TAD。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
Trizol试剂Vazyme--
HiScript III 一步法逆转录预混液(qPCR专用)Vazyme--
Tap Pro通用SYBR qPCR预混液Vazyme--
EDTAServicebio--
杜氏改良Eagle培养基(DMEM)Pricella--
胎牛血清(FBS)Pricella--
Lipofectamine 8000转染试剂Beyotime--
Veritas微孔板发光检测仪Promega--
PMEL抗体ABclonal--
牛血清白蛋白(BSA)Solarbio--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
基因组组装
HiFi测序覆盖度88×,Hi-C数据量584.72 Gb,组装软件hifiasm版本v0.19.3
阅读提示:Materials and Methods - Genome assembly and gene annotation; Results - Analyses
群体构建
亲本数量:30只北京鸭和120只连城鸭;F2群体数量1281只;表型分类标准
阅读提示:Materials and Methods - Ducks and sampling; Results - Inheritance of F2 population traits
全基因组重测序
样本数366,平均深度5×,测序平台Illumina HiSeq X-Ten,read长度150 bp
阅读提示:Materials and Methods - Whole genome resequencing
GWAS分析
样本数188,混合线性模型,固定效应包含性别和PCA前3个主成分,Bonferroni阈值0.01/N(N=8,887,194)
阅读提示:Materials and Methods - Genome-wide association analysis; Results - Genome-wide association analysis
转录组分析
样本数37,包括不同发育阶段皮肤和不同羽色羽球,测序平台Illumina,150 bp双端
阅读提示:Materials and Methods - Transcriptome sequencing and analysis; Supplementary Table S2
qPCR
引物序列,内参基因β-Actin,反应体系10 μL,重复3次
阅读提示:Materials and Methods - Quantitative PCR analysis; Supplementary Table S3
免疫荧光
组织固定于4%多聚甲醛,切片厚度5 µm,抗原修复EDTA 100°C 20分钟,一抗孵育4°C 12小时
阅读提示:Materials and Methods - Immunofluorescence experiment
荧光素酶报告基因
细胞密度0.5×10^5每孔,转染DNA量237.5 ng,内参pRL-TK 12.5 ng,重复3次
阅读提示:Materials and Methods - Luciferase reporter assay
Hi-C分析
分辨率为5 kb,TAD识别使用HiCPlotter默认算法
阅读提示:Materials and Methods - Hi-C sequencing and analysis; Results - Hi-C data显示TAD