Fig. 1. Neighbor-joining tree of allele-sharing distances calculated on the basis of 10529 hemizygous male-specific SNPs extracted from the VarGoats WGSs. A. Fourteen wild goats and five domestic goats with indicated haplogroup. B. 101 different haplotypes in 352 male domestic goats and 27 male bezoars and (not shown) markhor as outgroup. Individual codes consist of (i) Country codes, first two letters: AZ, Azerbaijan; BD, Bangladesh; CH, Switzerland; CN, China, DK, Denmark; ET, Ethiopia; FI, Finland; FR, France; IE, Ireland; IL, Israel; IR, Iran; IS, Iceland; IT, Italy; KE, Kenya; KO, Korea; MG, Madagascar; ML, Mali; MO, Morocco; MZ, Mozambique; NL, Netherlands; NZ, New Zealand; PK, Pakistan; RU, Russia; SD, Sudan, TR, Turkey; TJ, Tadjikistan; TZ, Tanzania; UZ, Uzbekistan; ZA, South Africa; ZW, Zimbabwe); (ii) Species code, 3rd-4th letter: CA (Capra aegagrus ), bezoar; CC (Capra caucasica , Capra cylindricornis ), West- and East-Caucasian tur, respectively; CF (Capra falconeri ), markhor; CH (Capra hircus ), domestic goat; CI (Capra ibex ) Alpine ibex; CN (Capra nubiana ), Nubian ibex; CS (Capra sibirica ), Siberian ibex.(iii) Codes for domestic breeds (three letters), or bezoar populations (four letters): see Figure 3. (iv)Sample number. Numbers near the nodes indicate bootstrapping values (%) after carrying out 100 bootstrap replicates. Y0, CaY1, CaY1AA, CaY2 and CaY2B indicate bezoar haplotypes.