《Animal Genetics》：Untangling the Copy Number Variation at the Basis of Belted Phenotypes in Cattle Using Long-Read Sequencing
编辑推荐：
牛中的腰带表型，即在荷兰腰带牛（Dutch Belted）和腰带加洛韦牛（Belted Galloway）等品种中观察到的横跨中部的独特白色带，是由TWIST2（一种参与黑素细胞发育的基因）上游的拷贝数变异（CNV）引起的。虽然该CNV是腰带形成所必需的，但腰
牛中的腰带表型，即在荷兰腰带牛（Dutch Belted）和腰带加洛韦牛（Belted Galloway）等品种中观察到的横跨中部的独特白色带，是由TWIST2（一种参与黑素细胞发育的基因）上游的拷贝数变异（CNV）引起的。虽然该CNV是腰带形成所必需的，但腰带宽度和完整度的变异遗传基础仍不清楚。为了细化对这些表型的理解，研究人员生成了一个荷兰腰带牛的从头基因组组装（de novo genome assembly），并分析了来自正常宽度、小宽度和不完整腰带动物的短读长和长读长测序数据。标准变异检测工具未能识别该CNV，突显了检测大插入片段的挑战性。相反，覆盖度分析和qPCR确认该CNV是一个位于TWIST2上游约16.5 kb处的6 kb串联重复阵列（tandem array）。长读长证据显示，正常宽度、小宽度和不完整腰带个体的单倍型完全相同（四个拷贝），表明CNV拷贝数并不影响腰带宽度或完整度。此外，在一个小宽度腰带个体中，未在CNV或TWIST2内检测到特有变异，提示影响腰带宽度的修饰因子很可能位于基因组的其他位置。全基因组扫描显示，CNV区域核苷酸多样性极低，这与荷兰腰带牛品种特有的强烈人工选择驱动的选择性清除（selective sweep）一致。综上所述，这些结果表明，虽然该CNV是腰带形成所必需的，但额外的修饰位点（modifier loci）可能影响腰带宽度和完整度。
论文解读文章
**研究背景与问题**
牛拥有多种毛色和斑纹，其中腰带表型（belted phenotype）尤为突出，表现为环绕动物中部的白色带，其位置、宽度和完整度在不同个体间存在差异。该表型在牛、小鼠、猪和山羊等物种中均有出现，但遗传机制不同。在牛中，早期研究将腰带表型定位到3号染色体端粒区域，并最终发现TWIST2（一种调节黑素细胞发育的转录因子）上游约16.5 kb处的一个6 kb拷贝数变异（CNV）是致病变异（Awasthi Mishra et al., 2017）。然而，过去的研究主要依赖短读长测序、SNP芯片强度和qPCR，这些方法难以准确解析复杂或重复序列丰富的结构变异（structural variant, SV），导致CNV的真实结构、单倍型组成及可能的致病变异尚未完全阐明。此外，腰带宽度和完整度的遗传基础仍不清楚，CNV拷贝数是否影响这些表型、是否存在额外的修饰位点（modifier loci）以及该区域是否受到选择压力，均是悬而未决的问题。由于荷兰腰带牛（Dutch Belted）品种有效群体规模小，且育种者对该品种的腰带表型有严格选择（要求腰带完整且宽度≥25 cm），导致近交严重，种群被列为“濒危”。因此，解析腰带表型背后的遗传机制，对于理解色素遗传学和支持育种实践至关重要。
**研究内容与结论**
为填补上述知识空白，研究人员生成了一个荷兰腰带牛的从头基因组组装（de novo genome assembly），并结合短读长（Illumina）和长读长（Oxford Nanopore Technologies, ONT）测序，对正常宽度腰带、小宽度腰带和不完整腰带个体进行了分析。主要结论包括：（1）标准结构变异检测工具（如smoove、delly、CNVnator、sniffles2）均未能识别该CNV，表明其检测难度大；（2）覆盖度分析和qPCR确认CNV位于ARS-UCD2.0参考基因组BTA3:118,028,480-118,034,892区域，为6 kb串联重复阵列；（3）长读长单倍型分析显示，所有可解析的单倍型均包含4个CNV拷贝，且不同表型（正常宽度、小宽度、不完整腰带）个体间单倍型完全一致，表明CNV拷贝数本身并不决定腰带宽度或完整度；（4）在小宽度腰带个体中，未在CNV或TWIST2内发现任何特有变异（SNP或SV），提示影响腰带宽度的修饰因子位于基因组其他区域；（5）核苷酸多样性分析显示，荷兰腰带牛中CNV区域的多样性极低（0.04百分位），远低于其他牛品种（如荷兰弗里生红牛为29.4百分位，格罗宁根白头牛为2.85百分位），且该低多样性区域延伸至CNV上游约20 Mb，符合由强烈人工选择驱动的选择性清除（selective sweep）特征。该研究发表在《Animal Genetics》。
**关键技术方法**
研究人员主要采用以下关键技术方法：（1）从头基因组组装（de novo genome assembly）：使用Flye v2.9.3对ONT长读长进行组装，并通过NextPolish v1.4.1结合ONT和Illumina短读长进行抛光，再用RagTag v2.1.0支架到ARS-UCD2.0参考基因组，得到荷兰腰带牛的高质量基因组组装（N50=102.8 Mb）。（2）CNV检测与定量：通过短读长覆盖度分析（Mosdepth）、qPCR（Comparative Ct法）和长读长BLAST比对，对CNV拷贝数进行多重验证。（3）变异发现：使用Clair3 v1.0.5和Sniffles2 v2.2分别对短读长和长读长数据进行SNP和SV检测，并筛选仅存在于小宽度腰带个体中而正常宽度个体中缺失的变异。（4）核苷酸多样性分析：使用VCFtools v0.1.17以25 kb窗口计算来自荷兰腰带牛、荷兰弗里生红牛和格罗宁根白头牛的短读长变异数据，得到π值。
**研究结果**
**3.1 在荷兰腰带牛中发现CNV的挑战**
通过多种结构变异检测工具（smoove、delly、CNVnator、sniffles2）对短读长和长读长数据进行调用，均未能检测到CNV。转而采用读段覆盖度分析，发现TWIST2上游约16.5 kb处存在一个约6 kb的区域，其覆盖度在短读长中增加约4倍，在长读长中增加约2倍，据此将CNV位置锚定到ARS-UCD2.0的BTA3:118,028,480-118,034,892。
**3.2 组装荷兰腰带牛基因组**
从头组装获得了高质量的荷兰腰带牛基因组，BUSCO完整单拷贝评分94.0%，N50达到102.8 Mb，高于参考基因组ARS-UCD2.0的26.4 Mb。重复序列分析显示CNV区域重复比例高达63.15%，高于基因组整体水平（43.18%）。由于拷贝间高度相似（尤其是第2和第3拷贝共享断裂点），组装仅重建出4个拷贝中的3个。
**3.3 拷贝数变异对腰带表型的影响**
通过覆盖度法、qPCR和长读长单倍型分析，评估了不同表型个体的CNV拷贝数。正常宽度腰带个体（n=12）拷贝数范围为5-11（均值8.0±1.3）；小宽度腰带个体（n=4）qPCR显示均值8.0±1.7（7或10拷贝）；不完整腰带个体（n=3）qPCR显示均值5.0±1.4（4或6拷贝）。但长读长分析显示，所有可解析单倍型（正常宽度、小宽度和不完整腰带个体）均含有4个CNV拷贝，未发现少于或多于4拷贝的单倍型。因此，CNV拷贝数本身不能解释腰带宽度或完整度的变异。
**3.4 宽度相关变异**
通过比较小宽度腰带个体与正常宽度腰带的从头组装，识别出201,669个SNPs和79,938个indels仅存在于小宽度腰带个体中，但这些变异均不位于CNV或TWIST2及其周围区域，也未发现基因组范围的大片段变异（如CNV、插入或倒位）。因此，宽度相关变异可能位于CNV区域之外，或涉及细微调控机制。
**3.5 CNV核苷酸多样性**
对11个荷兰腰带牛个体的短读长数据以25 kb窗口计算核苷酸多样性（π），结果显示CNV区域π值极低，位于全基因组第0.04百分位。对比荷兰弗里生红牛（29.4百分位）和格罗宁根白头牛（2.85百分位），表明该低多样性是荷兰腰带牛品种特有的。此外，CNV上游约20 Mb的染色体区域同样呈现低多样性，符合由强烈人工选择导致的选择性清除特征。
**总结与讨论**
研究人员确认TWIST2上游的CNV是腰带形成所必需的，但其拷贝数并不决定腰带宽度或完整度。长读长测序揭示所有腰带个体均携带4拷贝单倍型，未发现单倍型水平差异。在小宽度腰带个体中未发现CNV或TWIST2内的特有变异，提示修饰腰带宽度的因子可能位于基因组其他位点。CNV区域极度低的核苷酸多样性，且该现象为荷兰腰带牛品种所特有，表明该区域经历了一次由人工选择驱动的选择性清除，可能是一个祖先单倍型被选择性地扩增到高频率。未来工作需借助单倍型分辨的长读长测序和更大样本量来鉴定修饰位点，以平衡表型选择与遗传多样性维持。
**结论翻译**
本研究细化了对荷兰腰带牛腰带表型遗传架构的理解。研究人员确认TWIST2上游的CNV是腰带形成所必需的，但总CNV拷贝数本身并不能解释腰带宽度或完整度的变异。相反，研究人员推测基因组其他位置的额外修饰位点（modifier loci）可能影响腰带宽度和完整度。重要的是，研究人员未发现CNV或TWIST2内存在宽度决定变异的证据，提示其他基因组元件必须为此育种相关性状负责。CNV区域在荷兰腰带牛中表现出极低的核苷酸多样性，与封闭群体中强烈人工选择驱动的选择信号一致。因此，该模式表明腰带群体中发生了强烈的选择性清除（selective sweep），其中单一CNV单倍型被选择性地传播。未来工作应聚焦于单倍型分辨的长读长测序以全面表征CNV架构，并利用更大样本量鉴定腰带宽度的潜在修饰因子。这些见解不仅将推动对牛色素遗传学的理解，还将支持育种者开发遗传工具，以平衡腰带表型选择与遗传多样性维持。
