Abstract
The diversity of livestock coat color results from human positive selection and is an indispensable part of breed registration. As an important biodiversity resource, Asiatic wild ass has many special characteristics, including the most visualized feature, its yellowish-brown coat color, and excellent adaptation. To explore the genetic mechanisms of phenotypic characteristics in Asiatic wild ass and its hybrids, we resequenced the whole genome of one Mongolian Kulan (a subspecies of Asiatic wild ass) and 29 Kulan hybrids (Mongolian Kulan ♂×Xinjiang♀), and the ancestor composition indicated the true lineage of the hybrids. XP-EHH (Cross Population Extended Haplotype Homozygosity), θπ-ratio (Nucleotide Diversity Ratio), CLR (Composite Likelihood Ratio) and θπ (Nucleotide Diversity) methods were used to detect the candidate regions of positive selection in Asiatic wild ass and its hybrids. Several immune genes (DEFA1, DEFA5, DEFA7, GIMAP4, GIMAP1, IGLC1, IGLL5, GZMB and HLA) were observed by the CLR and θπ methods. XP-EHH and θπ-ratio revealed that these genes are potentially responsible for coat color (KITLG) and meat quality traits (PDE1B and MYLK2). Furthermore, the heatmap was able to show the clear difference in the haplotype of the KITLG gene between the Kulan hybrids and Asiatic wild ass group and the Guanzhong black donkey group, which is a powerful demonstration of the key role of KITLG in donkey color. Therefore, our study may provide new insights into the genetic basis of coat color, meat quality traits and immunity of Asiatic wild ass and its hybrids.