AUTHOR=Li Guangzhen , Luo Jing , Wang Fuwen , Xu Donghui , Ahmed Zulfiqar , Chen Shengmei , Li Ruizhe , Ma Zhijie 

TITLE=Whole-genome resequencing reveals genetic diversity, differentiation, and selection signatures of yak breeds/populations in Qinghai, China

JOURNAL=Frontiers in Genetics

VOLUME=13

YEAR=2023

URL=https://www.frontiersin.org/journals/genetics/articles/10.3389/fgene.2022.1034094

DOI=10.3389/fgene.2022.1034094

ISSN=1664-8021

ABSTRACT=<p>The Qinghai Province of China is located in the northeast region of the Qinghai–Tibetan Plateau (QTP) and carries abundant yak genetic resources. Previous investigations of archaeological records, mitochondrial DNA, and Y chromosomal markers have suggested that Qinghai was the major center of yak domestication. In the present study, we examined the genomic diversity, differentiation, and selection signatures of 113 Qinghai yak, including 42 newly sequenced Qinghai yak and 71 publicly available individuals, from nine yak breeds/populations (wild, Datong, Huanhu, Xueduo, Yushu, Qilian, Geermu, Tongde, and Huzhu white) using high-depth whole-genome resequencing data. We observed that most of Qinghai yak breeds/populations have abundant genomic diversity based on four genomic parameters (nucleotide diversity, inbreeding coefficients, linkage disequilibrium decay, and runs of homozygosity). Population genetic structure analysis showed that Qinghai yak have two lineages with two ancestral origins and that nine yak breeds/populations are clustered into three distinct groups of wild yak, Geermu yak, and seven other domestic yak breeds/populations. <italic>F</italic><sub>ST</sub> values showed moderate genetic differentiation between wild yak, Geermu yak, and the other Qinghai yak breeds/populations. Positive selection signals were detected in candidate genes associated with disease resistance (<italic>CDK2AP2</italic>, <italic>PLEC</italic>, and <italic>CYB5B</italic>), heat stress (<italic>NFAT5</italic>, <italic>HSF1</italic>, and <italic>SLC25A48</italic>), pigmentation (<italic>MCAM</italic>, <italic>RNF26</italic>, and <italic>BOP1</italic>), vision (<italic>C1QTNF5</italic>, <italic>MFRP</italic>, and <italic>TAX1BP3</italic>), milk quality (<italic>OPLAH</italic> and <italic>GRINA</italic>), neurodevelopment (<italic>SUSD4</italic>, <italic>INSYN1</italic>, and <italic>PPP1CA</italic>), and meat quality (<italic>ZRANB1</italic>), using the integrated PI, composite likelihood ratio (CLR), and <italic>F</italic><sub>ST</sub> methods. These findings offer new insights into the genetic mechanisms underlying target traits in yak and provide important information for understanding the genomic characteristics of yak breeds/populations in Qinghai.</p>