Background <p>Coat color was an economic trait in sheep, and identifying the causal variations affecting its diverse phenotypes had long been of interest. We performed whole-genome sequencing (WGS) to investigate the genetic composition and conducted a genome-wide association study (GWAS) for coat color using our lab data for 115 sheep (80 newly sequenced) from four Liangshan populations. Additionally, we retrieved public data for 104 sheep (from six other breeds) and 44 wild sheep.</p> Results <p><?tk 4?>Our study revealed that the Liangshan sheep population exhibited moderate genetic diversity and a closer genetic relationship with Yunnan sheep. Using WGS data, we performed a GWAS to identify genetic loci and causal variants influencing the black coat pigmentation in 115 Liangshan sheep (i.e., 85 white coat LSS vs. 30 black coat LSB). GWAS identified a 1.526&#xa0;Mb region (including the <i>MC1R</i> gene) on CHI14, spanning from 12,982,495&#xa0;bp to 14,508,891&#xa0;bp, as significantly associated with black coat (− log<sub>10</sub><i>P</i> &gt; 8.67). <i>MC1R</i> was highly expressed in the LSB group (mean = 15.394), significantly higher than in the LSS group (mean = 2.371; <i>P</i> = 0.032). A contingency table analysis revealed a highly significant difference in genotype frequencies at the <i>MC1R</i> mutation site c.361G &gt; A between the LSB and LSS populations (Fisher’s exact test, <i>P</i> = 8.268 × 10<sup>–25</sup>). Molecular docking results further demonstrated that the c.361G &gt; A mutation in <i>MC1R</i> of the LSB population disrupted a hydrogen bond between MC1R residue Asn121 and ASIP residue Phe119, reducing their molecular interaction. Haplotype analysis showed that 96.67% of LSB individuals carried the H2 haplotype, with strong linkage disequilibrium (<i>r</i><sup>2</sup> = 1) maintained across the region.</p> Conclusions <p>This study clarified the population structure of Liangshan sheep and confirmed the <i>MC1R</i> c.361G &gt; A mutation as a key locus influencing the black coat of LSB, potentially acting via the H2 haplotype.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Whole-genome resequencing revealed genetic diversity and the haplotype containing MC1R associated with black coat color in Liangshan sheep

  • Xueliang Sun,
  • Jiazhong Guo,
  • Jiahan Li,
  • Jing Luo,
  • Yiming Gong,
  • Guo Wu,
  • Mancheng Zhang,
  • Feng Xu,
  • Yuan Xu,
  • Chengqi Wei,
  • Li Li,
  • Ran Li,
  • Quzhe Emu,
  • Hongping Zhang

摘要

Background

Coat color was an economic trait in sheep, and identifying the causal variations affecting its diverse phenotypes had long been of interest. We performed whole-genome sequencing (WGS) to investigate the genetic composition and conducted a genome-wide association study (GWAS) for coat color using our lab data for 115 sheep (80 newly sequenced) from four Liangshan populations. Additionally, we retrieved public data for 104 sheep (from six other breeds) and 44 wild sheep.

Results

Our study revealed that the Liangshan sheep population exhibited moderate genetic diversity and a closer genetic relationship with Yunnan sheep. Using WGS data, we performed a GWAS to identify genetic loci and causal variants influencing the black coat pigmentation in 115 Liangshan sheep (i.e., 85 white coat LSS vs. 30 black coat LSB). GWAS identified a 1.526 Mb region (including the MC1R gene) on CHI14, spanning from 12,982,495 bp to 14,508,891 bp, as significantly associated with black coat (− log10P > 8.67). MC1R was highly expressed in the LSB group (mean = 15.394), significantly higher than in the LSS group (mean = 2.371; P = 0.032). A contingency table analysis revealed a highly significant difference in genotype frequencies at the MC1R mutation site c.361G > A between the LSB and LSS populations (Fisher’s exact test, P = 8.268 × 10–25). Molecular docking results further demonstrated that the c.361G > A mutation in MC1R of the LSB population disrupted a hydrogen bond between MC1R residue Asn121 and ASIP residue Phe119, reducing their molecular interaction. Haplotype analysis showed that 96.67% of LSB individuals carried the H2 haplotype, with strong linkage disequilibrium (r2 = 1) maintained across the region.

Conclusions

This study clarified the population structure of Liangshan sheep and confirmed the MC1R c.361G > A mutation as a key locus influencing the black coat of LSB, potentially acting via the H2 haplotype.