Donkey (Equus asinus) is an important livestock animal in China because of its draft and medicinal value. After a long period of natural and artificial selection, the variety and phenotype of donkeys have become abundant. We clarified the genetic and demographic characteristics of Chinese domestic donkeys and the selection pressures by analyzing 78 whole genomes from 12 breeds. According to population structure, most Chinese domestic donkeys showed a dominant ancestral type. However, the Chinese donkeys still represented a significant geographical distribution trend. In the selective sweep, gene annotation, functional enrichment, and differential expression analyses between large and small donkey groups, we identified selective signals, including NCAPG and LCORL, which are related to rapid growth and large body size. Our findings elucidate the evolutionary history and formation of different donkey breeds and provide theoretical insights into the genetic mechanism underlying breed characteristics and molecular breeding programs of donkey clades.
Al-Mamun, H.A., Kwan, P., Clark, S.A., Ferdosi, M.H., Tellam, R., Gondro, C., 2015. Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight. Genet. Sel. Evol. 47, 66.
|
Alexander, D.H., Novembre, J., Lange, K., 2009. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655-1664.
|
Allen, H.L., Estrada, K., Lettre, G., Berndt, S.I., Weedon, M.N., Rivadeneira, F., Willer, C.J., Jackson, A.U., Vedantam, S., Raychaudhuri, S., 2010. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature 467, 832-838.
|
Ardestani, S.S., Aminafshar, M., Maryam, M.B.Z.B., Banabazi, M.H., Sargolzaei, M., Miar, Y., 2020. A genome-wide signatures of selection study of Welsh ponies and draft horses revealed five genes associated with horse type variation. Gene Rep. 100833.
|
Beja-Pereira, A., England, P.R., Ferrand, N., Jordan, S., Bakhiet, A.O., Abdalla, M.A., Mashkour, M., Jordana, J., Taberlet, P., Luikart, G., African origins of the domestic donkey. Science 304, 1781-1781.
|
Boessneck, J., 1986. Eine equidenbestattung in spätfrühdynastischer zeit. Mitteilungen der Deutschen Orient-Gesellschaft zu Berlin 45-50.
|
Boessneck, J., von den Driesch, A., 1992. Tell el-dab'a: Tiere und historische umwelt im nordost-delta im 2. Jahrtausend v. Chr. Anhand der knochenfunde der ausgrabungen 1975-1986. Vii. Verlag der ₠ Osterreichischen Akademie der Wissenschaften 47-48.
|
Boessneck, J.v., von den Driesch, A., Eissa, A., 1992. Eine eselsbestattung der 1. Dynastie in abusir. Mitt. Dtsch. Archäol. Inst. Abt. Kairo 48, 1-10.
|
Boyko, A.R., Brooks, S.A., Behan-Braman, A., Castelhano, M., Corey, E., Oliveira, K.C., Swinburne, J.E., Todhunter, R.J., Zhang, Z., Ainsworth, D.M., 2014. Genomic analysis establishes correlation between growth and laryngeal neuropathy in thoroughbreds. BMC Genom. 15, 1-9.
|
Carty, C.L., Johnson, N.A., Hutter, C.M., Reiner, A.P., Peters, U., Tang, H., Kooperberg, C., 2012. Genome-wide association study of body height in african americans: the women’s health initiative SNP health association resource (SHARe). Hum. Mol. Genet. 21, 711-720.
|
Chang, H., 2011. Ainmal Genetic Resource in China: Horses, Donkeys and Camels. China Agricultural Press, China.
|
Chen, H., Patterson, N., Reich, D., 2010. Population differentiation as a test for selective sweeps. Genome Res. 20, 393-402.
|
Chen, J., Sun, Y., Manglai, D., Min, L., Pan, Q., 2010. Maternal genetic diversity and population structure of four Chinese donkey breeds. Livest. Sci. 131, 272-280.
|
Danecek, P., Auton, A., Abecasis, G., Albers, C.A., Banks, E., DePristo, M.A., Handsaker, R.E., Lunter, G., Marth, G.T., Sherry, S.T., 2011. The variant call format and VCFtools. Bioinformatics 27, 2156-2158.
|
Eberlein, A., Takasuga, A., Setoguchi, K., Pfuhl, R., Flisikowski, K., Fries, R., Klopp, N., Fürbass, R., Weikard, R., Kühn, C., 2009. Dissection of genetic factors modulating fetal growth in cattle indicates a substantial role of the non-SMC condensin I complex, subunit G (NCAPG) gene. Genetics 183, 951-964.
|
Falconer, D., Mackay, T., 1996. Introduction to quantitative genetics pearson education. Ltd, Essex, UK.
|
Ferenčaković, M., Hamzić, E., Gredler, B., Solberg, T., Klemetsdal, G., Curik, I., Sölkner, J., 2013. Estimates of autozygosity derived from runs of homozygosity: empirical evidence from selected cattle populations. J. Anim. Breed. Genet. 130, 286-293.
|
Fernandes, I., Bastien, Y., Wai, T., Nygard, K., Lin, R., Cormier, O., Lee, H.S., Eng, F., Bertos, N.R., Pelletier, N., 2003. Ligand-dependent nuclear receptor corepressor LCoR functions by histone deacetylase-dependent and-independent mechanisms. Mol. Cell 11, 139-150.
|
Gibson, J., Morton, N.E., Collins, A., 2006. Extended tracts of homozygosity in outbred human populations. Hum. Mol. Genet. 15, 789-795.
|
Gudbjartsson, D.F., Walters, G.B., Thorleifsson, G., Stefansson, H., Halldorsson, B.V., Zusmanovich, P., Sulem, P., Thorlacius, S., Gylfason, A., Steinberg, S., 2008. Many sequence variants affecting diversity of adult human height. Nat. Genet. 40, 609-615.
|
Horikoshi, M., Yaghootkar, H., Mook-Kanamori, D.O., Sovio, U., Taal, H.R., Hennig, B.J., Bradfield, J.P., St Pourcain, B., Evans, D.M., Charoen, P., 2013. New loci associated with birth weight identify genetic links between intrauterine growth and adult height and metabolism. Nat. Genet. 45, 76-82.
|
Huang, J., Zhao, Y., Bai, D., Shiraigol, W., Li, B., Yang, L., Wu, J., Bao, W., Ren, X., Jin, B., 2015. Donkey genome and insight into the imprinting of fast karyotype evolution. Sci. Rep. 5, 14106.
|
Kimura, B., Marshall, F.B., Chen, S., Rosenbom, S., Moehlman, P.D., Tuross, N., Sabin, R.C., Peters, J., Barich, B., Yohannes, H., 2011. Ancient DNA from Nubian and Somali wild ass provides insights into donkey ancestry and domestication. Proc. Biol. Sci. 278, 50-57.
|
Li, H., 2018. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094-3100.
|
Li, H., Durbin, R., 2009. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics 25, 1754-1760.
|
Li, Y., Fan, Y., Shaikh, A.S., Wang, Z., Wang, D., Tan, H., 2020. Dezhou donkey (Equus asinus) milk a potential treatment strategy for type 2 diabetes. J. Ethnopharmacol. 246, 112221.
|
Lindholm-Perry, A.K., Sexten, A.K., Kuehn, L.A., Smith, T.P., King, D.A., Shackelford, S.D., Wheeler, T.L., Ferrell, C.L., Jenkins, T.G., Snelling, W.M., 2011. Association, effects and validation of polymorphisms within the NCAPG-LCORL locus located on BTA6 with feed intake, gain, meat and carcass traits in beef cattle. BMC Genet. 12, 103.
|
Marras, G., Gaspa, G., Sorbolini, S., Dimauro, C., Ajmone-Marsan, P., Valentini, A., Williams, J.L., Macciotta, N.P., 2015. Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy. Anim. Genet. 46, 110-121.
|
Metzger, J., Schrimpf, R., Philipp, U., Distl, O., Expression levels of lcorl are associated with body size in horses. PLoS One e56497.
|
Miyake, Y., Yasui, M., Ikeda, K., KONDO, T., TSUKAMOTO, S., HORI, C., OKEMOTO, N., MASHOU, K., BANDO, R., NAKAMURA, N., 2009. Molecular cloning, a laboratory manual molecular cloning, a laboratory manual, 1989. Bioscience, biotechnology, and biochemistry, 73 1550-1560.
|
N'Diaye, A., Chen, G.K., Palmer, C.D., Ge, B., Tayo, B., Mathias, R.A., Ding, J., Nalls, M.A., Adeyemo, A., Adoue, V., Identification, replication, and fine-mapping of loci associated with adult height in individuals of African ancestry. PLoS Genet. e1002298.
|
Nekrutenko, A., Taylor, J., 2012. Next-generation sequencing data interpretation: enhancing reproducibility and accessibility. Nat. Rev. Genet. 13, 667-672.
|
Nilsson, M.A., Arnason, U., Spencer, P.B., Janke, A., 2004. Marsupial relationships and a timeline for marsupial radiation in south gondwana. Gene 340, 189-196.
|
Okada, Y., Kamatani, Y., Takahashi, A., Matsuda, K., Hosono, N., Ohmiya, H., Daigo, Y., Yamamoto, K., Kubo, M., Nakamura, Y., 2010. A genome-wide association study in 19 633 Japanese subjects identified LHX3-QSOX2 and IGF1 as adult height loci. Hum. Mol. Genet. 19, 2303-2312.
|
Paradis, E., 2010. Pegas: an R package for population genetics with an integrated-modular approach. Bioinformatics 26, 419-420.
|
Patterson, N., Price, A.L., Reich, D., Population structure and eigenanalysis. PLoS Genet. e190.
|
Petrie, W.M.F., 2013. Roman Portraits and Memphis Iv, Tarkhan I and Memphis V, Tarkhan Ii. Cambridge University Press, UK.
|
Pickrell, J., Pritchard, J., Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 8, e1002967.
|
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M.A., Bender, D., Maller, J., Sklar, P., De Bakker, P.I., Daly, M.J., 2007. Plink: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559-575.
|
Rossel, S., Marshall, F., Peters, J., Pilgram, T., Adams, M.D., O’Connor, D., 2008. Domestication of the donkey: timing, processes, and indicators. Proc. Natl. Acad. Sci. U. S. A. 105, 3715-3720.
|
Rubin, C.-J., Megens, H.-J., Barrio, A.M., Maqbool, K., Sayyab, S., Schwochow, D., Wang, C., Carlborg, Ö., Jern, P., Jørgensen, C.B., 2012. Strong signatures of selection in the domestic pig genome. Proc. Natl. Acad. Sci. U. S. A. 109, 19529-19536.
|
Setoguchi, K., Furuta, M., Hirano, T., Nagao, T., Watanabe, T., Sugimoto, Y., Takasuga, A., 2009. Cross-breed comparisons identified a critical 591-kb region for bovine carcass weight QTL (CW-2) on chromosome 6 and the Ile-442-Met substitution in NCAPG as a positional candidate. BMC Genet. 10, 1-12.
|
Setoguchi, K., Watanabe, T., Weikard, R., Albrecht, E., Kühn, C., Kinoshita, A., Sugimoto, Y., Takasuga, A., 2011. The SNP C.1326T>G in the non-SMC condensin I complex, subunit G (NCAPG) gene encoding a p. Ile442Met variant is associated with an increase in body frame size at puberty in cattle. Anim. Genet. 42, 650-655.
|
Signer-Hasler, H., Flury, C., Haase, B., Burger, D., Simianer, H., Leeb, T., Rieder, S., A genome-wide association study reveals loci influencing height and other conformation traits in horses. PLoS One e37282.
|
Silió, L., Rodríguez, M., Fernández, A., Barragán, C., Benítez, R., Óvilo, C., Fernández, A., 2013. Measuring inbreeding and inbreeding depression on pig growth from pedigree or SNP-derived metrics. J. Anim. Breed. Genet. 130, 349-360.
|
Snelling, W., Allan, M., Keele, J., Kuehn, L., Mcdaneld, T., Smith, T., Sonstegard, T., Thallman, R., Bennett, G., 2010. Genome-wide association study of growth in crossbred beef cattle. J. Anim. Sci. 88, 837-848.
|
Soranzo, N., Rivadeneira, F., Chinappen-Horsley, U., Malkina, I., Richards, J.B., Hammond, N., Stolk, L., Nica, A., Inouye, M., Hofman, A., Meta-analysis of genome-wide scans for human adult stature identifies novel loci and associations with measures of skeletal frame size. PLoS Genet. e1000445.
|
Sovio, U., Bennett, A.J., Millwood, I.Y., Molitor, J., O'Reilly, P.F., Timpson, N.J., Kaakinen, M., Laitinen, J., Haukka, J., Pillas, D., Genetic determinants of height growth assessed longitudinally from infancy to adulthood in the northern Finland birth cohort 1966. PLoS Genet. e1000409.
|
Takasuga, A., 2016. Plag1 and ncapg-lcorl in livestock. Anim. Sci. J. 87, 159-167.
|
Tetens, J., Widmann, P., Kühn, C., Thaller, G., 2013. A genome-wide association study indicates LCORL/NCAPG as a candidate locus for withers height in German Warmblood horses. Anim. Genet. 44, 467-471.
|
Uerpmann, H.-P., 1991. Equids in the ancient world. Naturwissenschaften 19, 2.
|
van der Valk, R.J., Kreiner-Møller, E., Kooijman, M.N., Guxens, M., Stergiakouli, E., Sääf, A., Bradfield, J.P., Geller, F., Hayes, M.G., Cousminer, D.L., 2015. A novel common variant in DCST2 is associated with length in early life and height in adulthood. Hum. Mol. Genet. 24, 1155-1168.
|
Vila, E., 2006. Equids in Time and Space. Oxbow Books, Oxford. von den Driesch, A., 1997. Tierreste aus buto im nidelta. Archaeofauna, 6, 23-29.
|
von den Driesch, A., 1997. Tierreste aus buto im nidelta. Archaeofauna, 6, 23-29.
|
Weedon, M.N., Lango, H., Lindgren, C.M., Wallace, C., Evans, D.M., Mangino, M., Freathy, R.M., Perry, J.R., Stevens, S., Hall, A.S., 2008. Genome-wide association analysis identifies 20 loci that influence adult height. Nat. Genet. 40, 575-583.
|
Weikard, R., Altmaier, E., Suhre, K., Weinberger, K.M., Hammon, H.M., Albrecht, E., Setoguchi, K., Takasuga, A., Kühn, C., 2010. Metabolomic profiles indicate distinct physiological pathways affected by two loci with major divergent effect on bos taurus growth and lipid deposition. Physiol. Genom. 42, 79-88.
|
Weir, B.S., Cockerham, C.C., 1984. Estimating F-statistics for the analysis of population structure. Evolution, 38, 1358-1370.
|
Zhang, Q., Su, R., Shan, C., Gao, C., Wu, P., 2018. Non-SMC condensin I complex, subunit G (NCAPG) is a novel mitotic gene required for hepatocellular cancer cell proliferation and migration. Oncol. Res. 26, 269-276.
|
Zhang, Y.-S., Yang, X.-Y., Wang, X.-B., Zhang, C.-M., Qin, F., Zhou, Z.-H., Lan, X.-Y., Chen, H., Lei, C.-Z., 2010. Cytochrome b genetic diversity and maternal origin of Chinese domestic donkey. Biochem. Genet. 48, 636-646.
|