Genetic architectures and selection signatures of body height in Chinese indigenous donkeys revealed by next-generation sequencing.

Donkeys are widely distributed labour animals in the world. During the process of the domestication and artificial selection of domestic donkeys, body sizes show significant differences among different breeds of donkeys. Based on the genome resequencing data of 103 Chinese indigenous donkeys from 11 breeds (Biyang, Dezhou, Guangling, Hetian, Jiami, Kulun, Qingyang, Turfan, Tibetan, Xinjiang, and Yunnan), seven Spanish donkeys from two breeds (Zamorano~Leonés and Andalusian), and three wild donkeys, we investigated the population structures of Chinese domestic donkeys with different body sizes. We used F and XP-EHH analyses to explore the selected regions related to body sizes. The results showed that Chinese indigenous donkeys have a closer relationship with African wild donkeys than with Asian wild donkeys. LCORL/NCAPG, FAM184B, TBX3, and IHH were identified as genes with strong signals in analysis of selection signature (F and XP-EHH) in large and small donkeys. The seven identified variants can be served as candidate loci affecting the body size of Chinese donkeys. Five of seven loci were located in intron 9 of FAM184B and were in a haplotype block, and one of the identified variants (Chr03:112664848) located in the CDS region of the LCORL gene was found to cause stop-loss. These candidate genes and variants shed new light on the molecular basis of donkey body size and will facilitate the breeding activities of donkeys.