Jukuri, open repository of the Natural Resources Institute Finland (Luke) All material supplied via Jukuri is protected by copyright and other intellectual property rights. Duplication or sale, in electronic or print form, of any part of the repository collections is prohibited. Making electronic or print copies of the material is permitted only for your own personal use or for educational purposes. For other purposes, this article may be used in accordance with the publisher’s terms. There may be differences between this version and the publisher’s version. You are advised to cite the publisher’s version. This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): VarGoats Consortium, Isaäc J. Nijman, Benjamin D. Rosen, Philippe Bardou, Thomas Faraut, Tristan Cumer, Kevin G. Daly, Zhuqing Zheng, Yudong Cai, Hojjat Asadollahpour, Bengi Çınar Kul, Wei-Yi Zhang, Guangxin E., A. Ayin, Hayley Baird, Meirat Bakhtin, Valentin A. Bâlteanu, Diana Barfield, Beate Berger …Johannes A. Lenstra Title: Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions Year: 2022 Version: Publisher’s version Copyright: The author(s) 2022 Rights: CC BY-NC-ND 4.0 Rights url: https://creativecommons.org/licenses/by-nc-nd/4.0/ Please cite the original version: VarGoats Consortium, Nijman, I. J., Rosen, B. D., Bardou, P., Faraut, T., Cumer, T., Daly, K. G., Zheng, Z., Cai, Y., Asadollahpour, H., Kul, B. Ç., Zhang, W-Y, E., G., Ayin, A., Baird, H., Bakhtin, M., Bâlteanu, V. A., Barfield, D., Berger, B. … Lenstra, J. A. (2022). Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions. Molecular Ecology, 31, 4364–4380. https://doi.org/10.1111/mec.16579 4364  |  Molecular Ecology. 2022;31:4364–4380.wileyonlinelibrary.com/journal/mec Received: 7 March 2022  | Revised: 12 June 2022  | Accepted: 20 June 2022 DOI: 10.1111/mec.16579 O R I G I N A L A R T I C L E Geographical contrasts of Y- chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions VarGoats Consortiuma | Isaäc J. Nijman1,2 | Benjamin D. Rosen3 | Philippe Bardou4 | Thomas Faraut4 | Tristan Cumer5 | Kevin G. Daly6 | Zhuqing Zheng7 | Yudong Cai7 | Hojjat Asadollahpour7 | Bengi Çınar Kul8 | Wei- Yi Zhang9 | Guangxin E.9  | A. Ayin10 | Hayley Baird11 | Meirat Bakhtin12 | Valentin A. Bâlteanu13 | Diana Barfield14 | Beate Berger15 | Thor Blichfeldt16 | Geert Boink17 | Sri R. A. Bugiwati18 | Zexi Cai19 | Sean Carolan20 | Emily Clark21 | Vlatka Cubric- Curik22 | Muhammad I. A. Dagong18 | Tashi Dorji23 | Louise Drew24 | Jiazhong Guo25 | Jón Hallsson26 | Simon Horvat27 | Juha Kantanen28 | Fuki Kawaguchi10 | Polat Kazymbet12 | Negar Khayatzadeh15 | Namshin Kim29 | Manoj Kumar Shah30 | Yuying Liao31 | Amparo Martínez32 | Joseph S. Masangkay33 | Maho Masaoka10 | Raffaele Mazza34 | John McEwan11 | Marco Milanesi35 | Faruque Md. Omar36 | Yuto Nomura10 | Nadjet- Amina Ouchene- Khelifi37 | Filipe Pereira38  | Goutam Sahana19 | Mazdak Salavati21 | Shinji Sasazaki10 | Anne Da Silva39  | Mojca Simčič27 | Johann Sölkner15 | Alison Sutherland40 | Johannes Tigchelaar17 | Hongpin Zhang25 | Econogene Consortiumb | Paolo Ajmone- Marsan41,42 | Daniel G. Bradley6 | Licia Colli41,43  | Cord Drögemüller44 | Yu Jiang7 | Chuzhao Lei7 | Hideyuki Mannen10 | François Pompanon5 | Gwenola Tosser- Klopp4  | Johannes A. Lenstra1 1Utrecht Univ., Utrecht, Netherlands 2Univ. Medical Center Utrecht, Utrecht Univ., Utrecht, The Netherlands 3Agriculture Research Service, Beltsville, Maryland, USA 4GenPhySE, Univ. Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France 5Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, LECA, Grenoble, France 6Trinity College Dublin, Dublin, Ireland 7College of Animal Science & Technology, Northwest A&F Univ., Yangling, China 8Faculty of Veterinary Medicine, Ankara Univ., Ankara, Turkey 9Southwest Univ., Chongqing, China 10Kobe Univ., Kobe, Japan 11AgResearch, Invermay Agricultural Centre, Mosgiel, New Zealand 12Astana Medical Univ., Astana, Kazakhstan This is an open access article under the terms of the Creative Commons Attribution- NonCommercial- NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non- commercial and no modifications or adaptations are made. © 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd. ahttp://www.goatgenome.org/vargoats.html bSee the list at the end of this paper.     |  4365VarGoats CoNsortIUM et al. 13Inst. of Life Sciences, Univ. Agricultural Sciences and Veterinary Medicine of Cluj- Napoca, Cluj- Napoca, Romania 14Pleystole, Qld, Australia 15Univ. Natural Resources and Life Sciences Vienna (BOKU), Vienna, Austria 16Norwegian Association of Sheep and Goat Breeders, Aas, Norway 17Stichting Zeldzame Huisdierrassen, Wageningen, The Netherlands 18Hasanuddin Univ., Makassar, Indonesia 19Aarhus Univ., Tjele, Denmark 20Old Irish Goat Society, Mulranny, Ireland 21Univ. Edinburgh, Edinburgh, UK 22Univ. Zagreb, Zagreb, Croatia 23International Centre for Integrated Mountain Development, Kathmandu, Nepal 24Finch Hatton, Qld, Australia 25Sichuan Agricultural Univ., Chengdu, China 26Agricultural Univ. Iceland, Reykjavik, Iceland 27Univ. Ljubljana, Biotechnical Faculty, Ljubljana, Slovenia 28Natural Resources Institute Finland (Luke), Jokioinen, Finland 29Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, South Korea 30National Animal Nutrition Research Centre, Lalitpur, Nepal 31Guangxi Key Laboratory of Livestock Genetic Improvement, Guangxi, China 32Universidad de Córdoba, Córdoba, Spain 33Univ. of the Philippines, Los Banos, Philippines 34Laboratorio Genetica e Servizi, Agrotis srl, Cremona, Italy 35Univ. degli studi della Tuscia, Viterbo, Italy 36Bangladesh Agricultural Univ., Mymensingh, Bangladesh 37Science Veterinary Inst., University of Blida, Blida, Algeria 38IDENTIFICA Genetic Testing Maia & Centre for Functional Ecology, Porto, Portugal 39E2LIM, Univ. Limoges, Limoges, France 40New Zealand Arapawa Goat Association, Masterton, New Zealand 41Univ. Cattolica del S. Cuore di Piacenza and BioDNA Biodiversity and Ancient DNA Res. Centre, Piacenza, Italy 42UCSC PRONUTRIGEN Nutrigenomics Res. Centre, Piacenza, Italy 43UCSC BioDNA Biodiversity and Ancient DNA Res. Centre, Piacenza, Italy 44Univ. Bern, Bern, Switzerland Correspondence Johannes A. Lenstra, Faculty of Veterinary Medicine, Utrecht University, IRAS, PO Box 80177, 3508 TD Utrecht, The Netherlands. Email: j.a.lenstra@uu.nl Funding information Croatian Science Foundation, Grant/ Award Number: ANAGRAMS- IP- 2018- 01- 8708; European Union, Grant/Award Number: QLK5- CT2001- 02461; France Génomique, Grant/Award Number: ANR- 10- INBS- 09- 08; the BBSRC Global Challenges Research Fund Data and Resources, Grant/Award Number: BBS/ OS/GC/000012F Handling Editor: David Coltman Abstract By their paternal transmission, Y- chromosomal haplotypes are sensitive markers of population history and male- mediated introgression. Previous studies identified bial- lelic single- nucleotide variants in the SRY, ZFY and DDX3Y genes, which in domestic goats identified four major Y- chromosomal haplotypes, Y1A, Y1B, Y2A and Y2B, with a marked geographical partitioning. Here, we extracted goat Y- chromosomal variants from whole- genome sequences of 386 domestic goats (75 breeds) and seven wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respec- tively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present- day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate an- cient paternal population bottlenecks and expansions during migrations into northern Europe, eastern and southern Asia, and Africa south of the Sahara. In addition, sharing of haplogroups indicates male- mediated introgressions, most notably an early gene 4366  |    VarGoats CoNsortIUM et al. 1  |  INTRODUC TION Goat (Capra hircus), sheep (Ovis aries), cattle (Bos taurus) and pigs (Sus scrofa) are four major livestock species, which after their domestica- tion in southwest Asia about 10,000 years ago (Larson & Fuller, 2014; Stiner et al., 2022) spread to all inhabited continents. Because of their relatively small size, sheep and goats were the earliest domes- ticates, but have become less important than cattle and pigs as sup- pliers of food. However, sheep and goat are suitable for extensive management by smallholders or hobby breeders with goats being favoured in conditions of poverty (Peacock, 2005). Although the high- quality goat cashmere wool and mohair fibres do not attain the volume of sheep wool, the demand for goat milk and cheese has in- creased considerably since the 1960s (Dubeuf et al., 2004; Miller & Lu, 2019). In the last four decades, this has doubled the global num- ber of goats to around one billion (Utaaker et al., 2021), approaching the numbers for sheep and cattle (Hegde, 2019). As for other livestock, genetic isolation, adaptation and selection have created numerous local goat populations, whereas a restricted number of high- performing breeds play a major role in agricultural production. The genetic diversity of goat breeds has been stud- ied extensively (Ajmone- Marsan et al., 2014; Amills et al., 2017; Deniskova et al., 2021; Zheng et al., 2020). This demonstrated for autosomal DNA a geographical partitioning of the diversity (Colli et al., 2018), which is in sharp contrast to similar studies of sheep (Belabdi et al., 2019; Ciani et al., 2020; Kijas et al., 2012), but has been confirmed based on similarities of ancient and modern DNA samples from the same regions (Cai et al., 2020; Daly et al., 2018, 2021). However, the goat Y- chromosome as a marker for paternal lineages has not yet been studied at a worldwide scale. Because of an absence of recombination, the male- specific part of the mammalian Y- chromosome is by far the longest haplotype that is stably transmitted across generations (Hughes et al., 2015). In many species, males have a relatively small male effective pop- ulation size, which makes Y- chromosomal variants highly informa- tive markers for paternal origin that generally show a much stronger phylogeographical differentiation than mitochondrial or autosomal variants. This is now widely exploited in population- genetic stud- ies of humans (Batini & Jobling, 2017; Jobling & Tyler- Smith, 2017; Kivisild, 2017), cattle (Edwards et al., 2011; Ganguly et al., 2020; Xia et al., 2019), horse (Felkel et al., 2019a; Wallner et al., 2017; Wutke et al., 2018), water buffalo (Zhang et al., 2016), sheep (Deng et al., 2020; Meadows & Kijas, 2009; Zhang et al., 2014), camel (Felkel et al., 2019b), llamas and alpacas (Marín et al., 2017), pigs (Choi et al., 2020; Guirao- Rico et al., 2018) and dogs (Natanaelsson et al., 2006; Oetjens et al., 2018). A preliminary analysis of the Y- chromosomal diversity in European and Turkish goats defined the three haplotypes, Y1A, Y1B and Y2, which had a strong geographical differentiation (Lenstra & Econogene Consortium, 2005). The same haplotypes were found in goats from Portugal and North Africa (Pereira et al., 2009), Turkey (Çinar Kul et al., 2015), eastern and southern Asia (Tabata et al., 2018, 2019; Waki et al., 2015), and Switzerland and Spain (Vidal et al., 2017), with the additional haplotypes Y2B in east Asia, Y2C in Turkish Hair and Kilis goats, and Y1B2 as well as Y1C mainly in Switzerland (Table S1). However, these haplotypes are based on a low number of single nucleotide polymorphisms (SNPs) in or near DDX3Y, SRY and ZFY and genotyping in a limited number of domestic goat breeds. Thus, it is not clear if the haplotypes represent major haplogroups or local variants or if other major haplogroups exist. Nor does it inform us on the Y- chromosomal variants that existed in earlier domestic goats or in their wild ancestor, the bezoar (Capra ae- gagrus; Amills et al., 2017). Whole- genome sequencing (WGS), how- ever, has confirmed the differentiation of the Y1 and Y2 haplogroups (Xiao et al., 2021; Zheng et al., 2020). In this study, we used WGS data for a large panel of goat breeds (Denoyelle et al., 2021) to systematically characterize the SNP- level variation in the single- copy male- specific part of the caprine Y- chromosome. In addition, we determined the Y- chromosomal haplogroups in goats originating from several European, Asian or African countries, in ancient goat DNA samples and in the wild be- zoar (Alberto et al., 2018; Cai et al., 2020; Daly et al., 2018, 2021; Zhang et al., 2014; Zheng et al., 2020). We sought to answer the following questions: (i) How are the WGS- based haplogroups related to the previously reported haplotypes? (ii) How are the domestic pa- ternal lineages related to those of bezoars from Iran and Anatolia, respectively? (iii) How strong is the phylogeographical structure of the caprine male lineages? (iv) What does the pattern of diversity tell us about Neolithic and later migrations? (v) Can we also infer other gene flows between or within continents? Answering these flow from Asian goats into Madagascar and the crossbreeding that in the 19th century resulted in the popular Boer and Anglo- Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y- chromosomal variants for reconstructing the history of domestic spe- cies with a wide geographical range. K E Y W O R D S domestication, goat, haplogroup, introgression, migration, phylogeography, Y- chromosome     |  4367VarGoats CoNsortIUM et al. questions will contribute to our understanding of the genetic back- ground of the domestic goat, which is relevant for breed manage- ment and conservation. 2  |  MATERIAL S AND METHODS 2.1  |  WGS data, filtering and phylogenetic analysis We selected as source of the SNPs four scaffolds that together cover 1,567,760 bp of the male- specific part of the caprine Y- chromosome. These are unplaced in the ARS1 assembly but for a large part closely match a recent Y- chromosomal contig of the Saanen_v1 assembly (Table S2; Li et al., 2021) and contain the single- copy Y- chromosomal genes SRY, DDX3Y and ZFY and the SNPs that define the major hap- lotypes Y1A, Y1B, Y2A and Y2B (Çinar Kul et al., 2015; Lenstra & Econogene Consortium, 2005; Waki et al., 2015; see Table S1). The genes USP9Y, UTY, DDX3Y and ZFY are proximate near one of the ends of the male- specific Y- chromosomal region and well separated from SRY (Li et al., 2021). The selected contigs have a low overall level of apparent heterozygosity, indicating a high frequency of hemizygous markers (Table S2). In a preliminary study (https://www.biorx iv.org/conte nt/biorx iv/early/ 2020/02/17/2020.02.17.952051.full.pdf), we used WGS data from the Sequence Read Archive (SRA) for 70 mainly Asian and Moroccan male goats (Alberto et al., 2018; Zheng et al., 2020; Table S3). We extracted the genotypes of 5356 SNPs as described (Zheng et al., 2020), which after filtering yielded 2350 SNPs without female- of male- heterozygous scores, <5% missing scores/SNP and a minor allele frequency (MAF) >0.02. For a more comprehensive global coverage, we used processed and mapped WGS data for 1372 goats generated by the VarGoats project (Denoyelle et al., 2021, www.goatg enome.org/vargo ats. html). These include published data available via the SRA, 61 of which have also been used in the preliminary study. A VCF file of 65,556 variants in the sequences covered by the four selected Y- chromosomal contigs (Table S2) was generated by using the gatk HaplotypeCaller in the GVCF mode and the goat ARS1 whole ge- nome sequence. After import into plink (version 1.90 beta) variants and samples were filtered by the following consecutive steps: 1. After removal of indels, 54,032 Y- chromosomal SNPs were retained. 2. From the 948 female goats, 670 were selected with scores of <1% for the 54,032 Y- chromosomal SNPs in order to minimize scores due to contamination with male DNA. 3. In total, 17,228 SNPs were scored in at least one of the 670 fe- males and were removed. 4. From the 36,804 remaining male- specific SNPs, 7263 SNPs had ≥1 heterozygous score and 506 were monomorphic in 424 male goats, keeping 29,035 hemizygous SNPs. 5. From the 380 domestic male samples, 354 with a call rate of >95% were kept. From the 34 wild goats, two Italian ibexes and four Iranian bezoars had call rates of only 90.7%– 94.4%, but this did not appear to affect their phylogenetic positions or the cor- responding bootstrapping values (see below), so these were re- tained in the data set. 6. From the 29,035 SNPs, 12,540 had a call rate of <99% in the 388 male goats (354 domestic and 34 wild) and were discarded. 7. From the remaining 16,495 SNPs, 552 SNPs had an MAF in the 354 male domestic goats of >1% and from the other 15,943, 9977 SNPs had a at least one score in wild goats, totalling 10,529 SNPs representing male- specific Y- chromosomal variation in domestic and/or wild goats, which differentiate 27 wild and 80 domestic haplotypes. 8. Finally, two domestic goats with unknown breed origins were dis- carded, resulting in a final panel of 386 (352 domestic and 34 wild) goats. Allele- sharing distances between individuals were calcu- lated using plink or mega7 (Tamura et al., 2011), and visualized in neighbour- joining (NJ) trees by using the program splitstree4 (Huson & Bryant, 2006). For calculating bootstrapping values in an NJ tree of the 80 domestic and 20 bezoar haplotypes and as an outgroup one markhor, we selected 2867 SNPs with MAF >5% and used the program mega7. The topology of this tree was essentially identical to a tree of genotypes of the 10,529 SNP panel for the same samples. For construction of median- joining networks (Bandelt et al., 1999), we selected 286 male goats and 27 bezoars, omitting other wild goats as well as transboundary breeds outside their re- gion of origin and balancing the breed representation by analysing ≤18 individuals per breed. With the number of polymorphic SNPs thus reduced to 1734, the program popart (Leigh & Bryant, 2015; http://popart.otago.ac.nz) generated a network of 91 haplotypes (Table S4). In addition to the SNPs previously by dideoxy sequencing (Çinar Kul et al., 2015; Lenstra & Econogene Consortium, 2005; Vidal et al., 2017; Waki et al., 2015), we identified other diagnostic SNPs in the VarGoats data set on the basis of an FST genetic distance (plink version 1.9 – fst) of 1.0 between males with a given haplogroup and all other males. Several Y- chromosomal diagnostic SNPs have been incorporated in the Goat_IGGC_65K_v2 bead array (Table S1). 2.2  |  Haplogroups distribution in 186 goat breeds We differentiated the major Y- chromosomal haplogroups Y1AA, Y1AB, Y1B, Y2A and Y2B for 2212 domestic goats from 186 breeds by combining data from several sources as detailed per breed in Table S5: 1. From the goat panel collected for the Econogene project (Lenstra & Econogene Consortium, 2005), DNA samples of 353 male goats from 38 European or southwestern Asian breeds were analysed by PCR (polymerase chain reaction) amplification and dideoxy- sequencing of DDX3Y, SRY and ZFY segments as 4368  |    VarGoats CoNsortIUM et al. described previously for bovine samples (Edwards et al., 2011; Nijman et al., 2008) using goat- specific primers (Table S6). 2. We used published data from five Portuguese breeds or from Moroccan goats (Pereira et al., 2009), 12 Asian breeds or national populations (Waki et al., 2015, Tabata et al., 2018, 2019, com- bined with unpublished data), eight Turkish breeds (Çinar Kul et al., 2015), and 26 Spanish and Swiss breeds (Vidal et al., 2017). 3. Haplogroup assignments for 354 domestic male goats in the VarGoats data set were derived from their positions in the phy- logenetic trees (Figure 1b; Figure S2), which were fully consistent with the alleles of the diagnostic SNPs (Table S1). This implies that they represent the variation that corresponds to the basal branch of the respective haplogroups in the phylogenetic trees or net- works (Figures 1b and 2). Haplogroups for VarGoats male goat with a call rate of <95% were inferred from the alleles of diagnos- tic SNPs (Table S1). 4. In total, 368 DNA samples from several sources, including the AdaptMap panel (Colli et al., 2018) were genotyped by the KASP assay (Kompetitive allele specific PCR assay) for SNPs NW_017189563.1 g.T280306>A and NW_017189885.1 g. (a) F I G U R E 1  Neighbour- joining tree of allele- sharing distances calculated on the basis of hemizygous male- specific SNPs extracted from the VarGoats WGSs. (a) Fourteen wild goats and five domestic goats; 10,529 SNPs. The wild, bezoar and five domestic haplogroups are indicated. (b) in total, 101 different haplotypes in 352 male domestic goats and 27 male bezoars and (not shown) markhor as an outgroup; 2867 SNPs (MAF >5%). 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     |  4369VarGoats CoNsortIUM et al. A11686>G, carried out at the Van Haeringen Laboratory (Wageningen, Netherlands).These SNPs differentiate Y1 vs. Y2 and Y1A vs. Y1B, respectively. 5. For 31 breeds, genotypes of diagnostic SNPs (Table S1) were obtained by megablast searching of SRA entries with queries of 40– 50 bp overlapping the SNPs. SRA data for pools of individuals were only used if this allowed an unambiguous identification of the haplogroup composition. For samples collected by Vidal et al. (2017) and Çinar Kul et al. (2015) and for samples analysed by KASP before genomic data became available, Y1AA and Y1AB have been both scored as the Y1A haplotype. Several of these samples did not contain Y1A or be- longed to breeds for which additional data are available (Table S5). However, for 17 breeds we only have the Y1A (Y1AA + Y1AB) fre- quency. Likewise, Vidal et al. (2017) and the KASP assays did not dif- ferentiate Y2A and Y2B. Since the VarGoats and Econogene panels with comprehensive coverage of Europe and Africa did not contain a single Y2B- carrying goat, we assigned Y2 scores in other European and African goats to Y2A. 2.3  |  Haplotypes of ancient DNA samples Daly et al. (2021) assigned ancient DNA (aDNA) samples from south- west Asia and southeast Europe to the domestic Y- chromosomal haplogroups on the basis of their positions in a phylogenetic tree. Cai et al. (2020) and Zheng et al. (2020) described eight Chinese aDNA samples and one medieval sample from the northern Caucasian region. Sample YJL2G (coverage 13.4×) clustered with the Y1AB haplogroup, but the coverage of the others ranged from 0.013 to 0.118×. This re- sulted in insufficient overlap with the 10,529 male- specific hemizy- gous SNPs in domestic goats (see above). Therefore, we relaxed our filtering and allowed SNPs scoring in up to 1% of the females. We ex- cluded SNPs with male heterozygote scores only if these occurred in a panel of individuals with call rates >95% representing the 80 different domestic haplotypes (see above) but removing two individuals with a high heterozygosity. This resulted in 5593 SNPs, 1018 of which were also scored in the low- coverage aDNA samples. The combined phylo- genetic signals (Table S7) showed for sample GTM6G a high proportion of inconsistent scores, presumably due to contamination, but allowed plausible haplogroup assignment for samples KA1G, SMG1, SMG7, SMG11 and YJL2G and a tentative assignment for BG3 (Table S7). 3  |  RESULTS 3.1  |  Phylogeny of Y- chromosomal haplogroups A phylogenetic tree of the wild and domestic goats (Figure 1a) shows an intermingling of bezoars and the domestic goat. From the other wild goat species, the markhor is the closest relative of the bezoar and the domestic goat. We found 107 different haplotypes in our panel of 352 domestic goats and 27 bezoars with markhor as an outgroup. A phylogenetic tree (Figure 1b) shows haplogroups that correspond to the haplo- types Y1B, Y2A (Lenstra & Econogene Consortium, 2005) and Y2B (Waki et al., 2015) whereas the Y1A haplotypes are split into hap- logroups Y1AA and Y1AB. This is confirmed in a data set of mainly Asian and Moroccan goats (Figure S1). All domestic haplotypes dif- fer from the Iranian or Anatolian bezoar haplotypes, which also dif- fer from each other. The bezoar haplotypes are associated with the domestic Y1AA or Y2B clusters (CaY1AA and CaY2B, respectively), are linked to the Y1 or Y2 roots (CaY1 and CaY2) or are outside the domestic cluster (Y0). Figure S2 shows subtrees containing all 352 domestic goats. This figure also indicates goats with the previously proposed local hap- lotypes (Vidal et al., 2017; this study, Table S1): Y1AB2 (this study), Y1B2 and Y1C. Y1B2 is represented by Swiss, French and Dutch goats. Our goat panel does not contain goats with the Y2C diagnos- tic allele (Çinar Kul et al., 2015). The phylogenetic relationships of domestic and bezoar haplo- types are confirmed by a median- joining network (MJN, Figure 2). Both the NJ tree and the MJN allow a few interesting observa- tions to be made: (i) a close relationship of Y1B sequences from Switzerland and Korean Native goats, suggesting recent crossbreed- ing; (ii) likewise, a close relationship of Y1AA sequences from Central as well as South Asia, from South African Boer goats and from other Y1AA goats in southern and eastern Africa; and (iii) a clear diver- gence of the Malagasy Y2A haplotypes from the African continental haplotypes. 3.2  |  Geographical distribution of haplogroups Figures 2 and 3 show the haplogroup distributions in 80 breeds from Europe, Asia and Africa. Interestingly, there is clear difference be- tween the ancient and modern DNA samples. All five haplogroups were detected in four ancient goat samples from southeastern Europe and 17 from southwestern Asia (Daly et al., 2018, 2021; Figure 3a; Table S5). In contrast, we observe a strong geographical differentiation of domestic goats: 1. Haplogroup Y1B is predominant in central and northern Europe, but outside Europe and North Africa it has only been found in one Ugandan Karamonja goat, in the Korean native breed and in exported Saanen populations. 2. In northern and central Europe Y2A is only found in a single French des Fossés sample and together with Y1AA in the cross- bred Anglo- Nubian. It is the predominant haplogroup in Spain, Anatolia and Africa south of the Sahara, but it is not found in China or Southeast Asia except in about 20% of the goats in the Philippines and Sulawesi. Remarkably, most Y2A haplotypes on Madagascar are more closely related to Asian than to continental- African Y2A haplotypes. 4370  |    VarGoats CoNsortIUM et al. F I G U R E 1   (Continued) (b)     |  4371VarGoats CoNsortIUM et al. 3. Y2B is absent in Europe, continental Africa and west Asia, but is a major haplotype in east and southeast Asia. It is also observed in one Malagasy Diana sheep. 4. Y1AA in Europe is only represented by three haplotypes in the local Ciocara breed and the Montecristo island population (Somenzi et al., 2022). The Italian haplotypes are outside the clus- ter of closely related south Asian and southeast African Y1AA haplotypes (Figure 2; Figure S2). 5. The available data suggest a contrast of Y1AB dominating in northern China and Y1AA with Y2B in the south. 4  |  DISCUSSION One of the benefits of the current availability of WGS data sets is the access to an abundance of sequence variants, which allow a compar- ison of individuals or populations for several purposes. This is espe- cially useful for the analysis of Y- chromosomal diversity, which was previously restricted by the availability of Y- chromosomal markers. The male- specific part of the Y- chromosome constitutes the long- est haplotype in the mammalian genome and may serve as a marker for mammalian paternal lineages. Here, we combined the data set of the VarGoats project with published data and genotyped diagnostic SNPs in male goat samples from several sources. The Y- chromosomal phylogeny of wild and domestic goats is in agreement with the Y- chromosomal tree on the basis of AMELY and ZFY gene fragments (Pidancier et al., 2006) and with a phylogeny of WGS sequences (Grossen et al., 2020; Cai et al., unpublished). Mitochondrial DNA (mtDNA) trees confirm the close relationship of markhor (Capra falconeri) with bezoar and domestic goat, but do not show the separation of these species and the other wild goats. In addition, mtDNA sequences of some, but not all East Caucasian turs (Capra cylindricornis) cluster with the mtDNA sequences of markhor, bezoar and domestic goat, illustrating a separate history of maternal and paternal lineages in cross- fertile species (Chen et al., 2018; Marín et al., 2017; O'Connell et al., 2014; Zhang et al., 2016, 2020). On the basis of WGS data, Zheng et al. (2020) and Xiao et al. (2021) reproduced the divergence of the domestic Y1 and Y2 haplogroups previously found on the basis of SNPs within or near Y- chromosomal genes (Lenstra & Econogene Consortium, 2005). Here we report a further differentiation of haplogroups, resulting in a phylogeny sup- ported by a largely independent preliminary data set (Figure S1) and two phylogenetic algorithms. We found that the major haplogroups F I G U R E 2  Median- joining network of 286 male domestic goats and 27 bezoars, omitting transboundary breeds outside their region of origin and balancing the breed representation by reducing the sample size of ALP and SAA to 18 individuals. This panel has 1734 SNPs with MAF >0.01, and 91 different haplotypes (Table S4) 4372  |    VarGoats CoNsortIUM et al. correspond to haplotypes defined by SNPs (Çinar Kul et al., 2015; Lenstra & Econogene Consortium, 2005; Pereira et al., 2009; Tabata et al., 2018, 2019; Vidal et al., 2017; Waki et al., 2015), but Y1A hap- lotypes belong to either haplogroup Y1AA or Y1AB. The phylogeny also indicates that these haplogroups diverged after the split of the markhor and the cluster of the wild bezoar and domestic goats. The domestic goats and the two bezoar populations from Anatolia and Iran do not share haplotypes, whereas the bezoar haplotypes are attached to deep nodes in the tree of mainly domes- tic haplotypes. This suggests an absence of male gene flow between the bezoar populations and between the bezoar and domestic goats from the same region. Thus, domestic goats, which possibly were de- rived from bezoar populations not sampled in this study, maintained their paternal lineages during migration from the Fertile Crescent via Anatolia to Europe, despite indications of management of wild goats in central Anatolia (Stiner et al., 2022). Geographical plots of domestic haplogroup frequencies show a considerable spatial differentiation, which resonates with the strong phylogeography displayed by autosomal SNPs (Colli et al., 2018), but is in clear contrast to the weak phylogenetic struc- ture displayed by the major mtDNA haplogroups (Colli et al., 2018; Luikart et al., 2001; Naderi et al., 2007, 2008; Zhao, Zhong, et al., 2014; Zhao, Zhao, et al., 2014). Remarkably, Y- chromosomal haplotypes from all five haplogroups have been found in aDNA samples from Southwest Asia and Southeast Europe. The locations of ancient Y1AA, Y1B and Y2B samples are well outside the range of the corresponding domestic haplogroups (Figure 3). This indi- cates that during the Neolithic and later worldwide migrations a F I G U R E 3   (Continued) BZAl BZZa BZAz Semnan Geor Azer Bulak Acem BZTr Blago n Direkli bezoar Chalow Gans Daher Abdul Hosein BZFa Monjukli BZMaY0 Y1AA, CaY1AA Y1AB, CaY1 Y1B Y2A, CaY2 Y2B, CaY2B 1 sample or 1 pool of samples/breed bezoar Preneolithic bezoar Neolithic Chalcolithic Iron Age/Medieval North- Caucasian DUK CAP ARG BIO BAL BEY HAS LIQ GAR GIR BLA MATGMO ORO FLR MLG PIZ TAS SAR VAL PYY RASVRT MUZ DKL TWZ FIN DPG DUL ALP APP COR PYR ROV GRG SKO BUK BRV SER SGBTOG VAG CRP SAA ABA GUR HAI HNM KLS NDZ TER ICL MUG GDR CHQ SWE CHE MULARR BIL NCG NDG BLO DRZ CRO SSG MOR ARB NDK ANB BCD SRP FSS PVC LOR PTV SAV CCG IST GRS Y1AA or Y1AB Y1AA Y1AB Y1B Y2A Y2B ANK BAG CAG (a) (b)     |  4373VarGoats CoNsortIUM et al. series of bottlenecks and expansions in the domestic male lineage created a strong geographical differentiation of the haplogroup distribution (Figure 4): 1. The dominance of haplogroup Y1B in central and northern Europe may very well reflect population bottlenecks during the Neolithic introduction of agriculture via the Danube route (Cymbron et al., 2005; Rivollat et al., 2015; Tresset & Vigne, 2007) 2. Y2A and Y1AA are almost the only haplogroups in Africa south of the Sahara. The two African Y1AA haplotypes are related to those of Asia, indicating that only Y2A expanded during the first introduction of domesticated goats in central and southern Africa 3. Y2B has been found in two Neolithic Iranian samples whereas re- lated CaYB2 haplotypes are present in Iranian and Anatolian be- zoars. However, as a result of population bottlenecks during the global spread of domestic goats, Y2B now occurs only in Asia east of the Indus River and in one goat from Madagascar (see below). 4. Y1AA was found in Neolithic samples in southeast Europe, but now has a low frequency in Europe. In Asia it expanded together with Y2B and later came to South Africa when Asian goats were used to breed the Boer goat (see below). 5. Y1AB is the most frequent haplogroup in north China. The dis- tributions of Y1AB and Y1AA/Y2B in East Asia correspond to ranges of the north Chinese cashmere goats and the small Southeast Asian “katjang” type, respectively (Porter et al., 2016). (c) (d) F I G U R E 3   (Continued) 4374  |    VarGoats CoNsortIUM et al. This obviously reflects the large difference in climate between northern and southern China, which determined a similar distri- bution of taurine and indicine cattle. These two types of cattle are considered to have entered China via a northern and southern mi- gration route, respectively (Chen et al., 2018; Zhang et al., 2020), supporting the separate eastern expansions of the Y1AB and Y1AA/Y2B goats, respectively. Exceptions to these geographical patterns follow from close re- lationships between haplotypes from different continents, which are probably explained by later major introgressions. Interestingly, in the phylogenetic trees and networks (Figure 2; Figure S2) the Y2A haplotypes on Madagascar are closely related to Asian haplotypes and one Diana goat from northern Madagascar even has an east- ern Asian Y2B haplotype. However, autosomal DNA shows that the Malagasy goats are more closely related to the southern and eastern African continental goats (Colli et al., 2018; Denoyelle et al., 2021). This parallels a recent finding that Malagasy cattle combine Indian and admixed African zebu ancestry (Magnier et al., 2022). The F I G U R E 3  Haplogroup distributions for 186 breeds (2212 goats), 26 aDNA samples and 31 bezoars. (a) Ancient DNA samples and bezoars; (b) European breeds; (c) Asian breeds; (d) African breeds. Breed or bezoar population codes: ABA, Abaza; ABR, Abergelle; ALG, Algarve; ALP, alpine; ANB, Anglo- Nubian; AND, Androy; ANK, angora; APP, Appenzell; ARB, Arabia; ARG, Argentata dell'Etna; ARR, arran; ARW, Arapawa; BAG, Bagot; BAL, Balearic; BBL, black Bengal; BCD, British composite dairy; BEE, Beetal black Pakistan; BEY, Bermeya; Bhutan; BIO, Bionda dell'Amadello; BLA, Blanca Andaluzza; BLB, Bilbery; BLO, Blobe; BOE, Boer; BRI, Barbari (Bari); BRV, Bravia; BUK, polish fawn Coloured; BUL, Balaka- Ulongwe; BZAl, Algorz bezoar; BAAz, Azerbaijan bezoar; BZFa, Fars bezoar; BZMa; Markazi bezoar; BZTr, Anatolian bezoar; BZZa, Zagros bezoar; CAG, Capra Grigia; CAM, Cambodjan; CAP, Capore; CBR, Chengdu Brown; CCG, Ciociara; CDM, Chaidamu Qinhai; CHE, Cheviot; CHM, Chongmin (Shanghai); CHQ, Charnequeira; CHV, Cheviot; COR, Corsican; CPC, Changtangi Pashmira cashmere; CRO, Croatian spotted; CRP, Carpathian; Da'er, Sichuan; DAM, Damani; DDZ, Dedza; DIA, Diana; DJA, Djallonke; DKL, Danish landrace; DPG, Dutch pied goat; DRZ, Dreznica; DSA, Dutch Saanen; DUA, Du'an; DUK, Dukati; DUL, Dutch Nordic goat; ESF, Esfahan; FIN, Finnish; FLR, Florida; FSS, de Fossés; FUE, Fuenteventura; GLL, Galla; GAR, Garganica; GDR, Guadarrama; GIR, Girgentata; GMO, Grigia Molisana; GRG, Greek; GRS, Grisons striped; GUA, Guadeloupe creole; GUE, Guéra; GUM; GUR, Gürcü; GUM, Gumez; HAI, hair (Anatolian black); HAS, Hasi; HMN, Haimen (Jiangsu); HNM, Honamli; ICL, Icelandic; IMC, Inner Mongolia cashmere; IND, Indian; IRA, Iran; IRQ, Iraq; IST, Istrian; JAT, Jatall; JBL, Jintang black; JSA, Japanese Saanen; KCC, Kacchan; KEF, Keffa; KHA, Khazakhstani; KIG, Kigezi; KLS, Kilis; KMO, Karamonja; KMR, Kamorr; KNG, Korean native; KSA, Korean Saanen; LAN, Landim; LAO, Laos; LBA, Lori- Bakhtiari; lCL, Icelandic; LGL, Longlin; LIL, Lilongwe; LIQ, Liqenasi; LOR, Lorraine; MAC, Machen black; MAK, Makatia; MAT, Mati; MAU, Maure; MEN, Menabe; MGL, Mongolian; MLG, Malagueña; MLI, Mali; MLT, Maltese; MLW, Malawi; MLY, Malya; MNT, Montecristo; MOR, Moroccan; MSH, Mashona; MTB, Matebele; MUB, Mubende; MUG, Murciano Granadina; MUL, Mulranny; MUZ, Muzhake; MYA, Myanmar; MZA, M'Zabite; NAC, Nacki; NCG, Norwegian coastal; NDG, Norwegian dairy; NDK, Naine de Kabylie; NDZ, Norduz; NEP, Nepalese; NGD, Nganda; NJI, Nanjiang yellow; ORO, Orobica; PCG, Peacock; PAT, Pateri; PAW, pare- white; PER, Percy Island; PEU, Peulh; PHI, Pilippine; PIZ, Pinzgauer; PRW, pare- white; PTV, Poitevin; PVC, Provençale; PYR, Pyrenean; PYY, Payoya; QIN, Qin; QHI, Qinhai; RAN, rangeland; RAO, Henan Raoshan white; RAS, Rasquera; ROV, rove; RSK, Nigerian Maradi (red Sokoto); SAA, Saanen; SAR, Sarda; SAV, Savoye; SCA, Shaanbei cashmere; SDN, Sudan; SEA, small east African; SEB, Sebei; SER, Serrana; SGB, St Gallen booted; SHI, Shjiba; SHL, Shahel; SHN, Nigerian Sahel; SKO, Skopelos; SNJ, Sanja; SOF, Sofia; SOU, southwest; SRP, Serpentina; SSG, Steirische Schecken; SSU, South Sulawesi; SWE, Swedish landrace (Jämtland); TAS, Tauernschecken; TED, Teddi; TER, Teramo; THA, Thari; TIB, Tibetan; TNF, Tinerfena; TOG, Toggenburg; TWZ, Thuringian Forest; VAG, Valais Blackneck; VAG, Valais Blackneck; VAL, Valdostana; VIE, Vietnamese; VRT, Verata; VZC, Verzasca; WAD, west African dwarf; XJI, Xinjiang; XUH, Xuhuai; YBL, Youzhou black skin; YIM, Shandong Yimeng white; YIW, Yichang white; ZAR, Zaraiba F I G U R E 4  Bottlenecks, migrations and introgressions of domestic goat. This figure does not indicate export of transboundary breeds Y1B Y1AA Y1B 20-21th century After 1990 Expansion after bottleneck introgressions Y1B Y1AA Y1AA     |  4375VarGoats CoNsortIUM et al. Malagasy language has an Austronesian origin, which testifies to the colonization of Madagascar by immigrants from southeastern Asia about 500 CE. Thus, it is likely that these immigrants brought Austronesian goats, cattle and possibly also other livestock from their region of origin to Madagascar. Other introgressions are more recent. The exceptional Y1AA and Y2A in the English Anglo- Nubian is explained by the documented im- port during the 19th century of Indian and African goats to England. These served on the ships as sources of milk and meat, but surviving males were crossed with English goats, which resulted in the emer- gence of a popular transboundary breed. The worldwide popular Boer goat also is of mixed origin (Porter et al., 2016; Vidal et al., 2017) and carries exclusively Y1AA haplo- types. This breed is considered to be a crossbred of local African and Indian goats, possibly mediated by incrossing of Anglo- Nubian males (Porter et al., 2016). The crossbred origin is consistent with the results of Colli et al. (2018): a separate phylogenetic position of the Boer relative the other African and Asian goats and a K = 3 pattern of model- based clustering showing African and west Asian ancestry. The Indian ancestry is entirely in agreement with a close clustering of the Boer and Pakistani Y1AA haplotypes (Figures 1b and 2). Subsequently, the Boer became itself a source of introgression. The same Y1AA haplotypes are closely related to Y1AA haplotypes in local breeds in Uganda, Malawi, Mozambique and Zimbabwe. In these countries crossbreeding with Boer goats from Africa is pop- ular because of its excellent meat production (Banda et al., 1993; Garrine, 2007; Lu, 2011; Onzima et al., 2018). Therefore, it is most likely that the Y1AA haplotype in eastern and southern African goats originates from the Boer goat. There were three out- of- range findings of Y1B, in the Ugandan Karamonja, in the Korean native goat and in the indigenous goats kept on Chongmin Island in Shanghai. Because of the popular- ity of Swiss dairy goats in both Uganda (NAADS, 2005) and Korea (Kim et al., 2019), crossbreeding again is the most likely explana- tion. Although European admixture in the Chongmin goats (Gao et al., 2020) has not been reported, the exotic occurrence of Y chro- mosomal variants appears to be a direct and sensitive indicator of admixture events. These need to be complemented with quantitative admixture tests, such as model- based clustering, the f3 and f4 test or, ideally, identification of introgressed segments across the genome. The latter approach may also lead to clues regarding the phe- notypic consequences of introgression via the identification of the admixed genes (Chen et al., 2018; Lv et al., 2014; Wang et al., 2015; Zheng et al., 2020). A more direct link with Y- chromosomal varia- tion would be provided if this can be linked to male phenotypic traits, but even in human genetics this has scarcely be investigated (Matsunaga et al., 2021; Yang et al., 2018; Zhang et al., 2021). Breeds in which different Y- chromosomal haplogroups occur may allow us to study an association of haplogroups with typically male traits such as male fertility and dominance behaviour. It would be interesting to see if Y- chromosomal variants can be related to climate or other environmental features, because this would imply that the geo- graphical differentiation of the Y- chromosomal variation is driven by regional adaptation. Most introgressions described in this study contribute to the expansion of popular breeds at the expense of the original local breeds. On the one hand, depending on the extent of gene flow this may decrease the diversity of the genetic resources; on the other hand, it does not necessarily disrupt the environmental adaptation, arguably one of the most important components of the phenotypic repertoire. If properly managed, admixture of productive breeds may also contribute to the sustainable conservation of local popu- lations and illustrates that genetic diversity has never been a static phenomenon. We conclude that the Y- chromosomal variation of goats reveals bottlenecks, expansions and introgressions, illustrating the power of Y- chromosomal markers for inferring the genetic origin of mamma- lian populations. AUTHOR CONTRIBUTIONS IJN, PAM and JAL designed the study; GT- K and LC coordinated the VarGoats project; IJN and BÇK carried out the ABI sequenc- ing; BDR, PB, TF, ZZ, YJ, YC, ZC, MMi and GS analysed the WGS data; TC, FPo and CD supplied the bezoar genotypes; KGD, DGB, YC and YJ supplied the aDNA data; HM, FK, SS, MMa, YN, AA, JSM, IAD, SRAB, FMD, TD, MKS, MB and PK provided genotypes for south Asian goats; BDR supplied most of the African samples for KASP genotyping; VAB, DB, BB, TB, SC, VC- C, LD, JG, JH, JK, NKh, NK, AM, RM, JMc, NAOC, FPe, AdS, MS, JS, AS, JT and HZ collected material and/or data for other breeds; JAL and IJN per- formed the downstream analysis; JAL wrote the first draft; and KGD, AM, FPe, BÇK, JMc, MM, MS, PAM, LC, CD and GTK con- tributed to the text. ACKNOWLEDG EMENTS The VarGoats project sequencing effort has been funded mainly by France Génomique “Call for high impact projects” (ANR- 10- INBS- 09- 08) and by the BBSRC Global Challenges Research Fund Data and Resources Grant BBS/OS/GC/000012F awarded to the Roslin Institute for Reference genome and population sequencing of African goats. Detailed funding sources of the VarGoats project are described in Denoyelle et al. (2021). This study was supported by the Croatian Science Foundation (Project ANAGRAMS- IP- 2018- 01- 8708) “Application of NGS in assessment of genomic variability in ruminants” and by the European Union (projects ECONOGENE QLK5– CT2001– 02461). We are grateful to Dr E. Cuppen (Utrecht Medical Centre) for access to dideoxy sequenc- ing facilities. CONFLIC T OF INTERE S T The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. 4376  |    VarGoats CoNsortIUM et al. VARGOATS CONSORTIUM Australia: James Kijas, CSIRO. Denmark: Bernt Guldbrandtsen, Aarhus University; Finland: Juha Kantanen, Luke. France: Gwenola Tosser- Klopp, Philippe Bardou, Laure Denoyelle, Thomas Faraut, Julien Sarry, Estel le Talouarn, GenPhySE, Castanet- Tolosan; Adriana Alberti, Céline Orvain, Stefan Engelen, Université Paris- Saclay, Evry; Dylan Duby, Museum National d'Histoire Naturelle; Pierre Martin, Capgenes; Coralie Danchin, Delphine Duclos, Institut de l'Elevage; Daniel Allain, Remy Arquet, Nathalie Mandonnet, Michel Naves, Isabelle Palhiere, Rachel Rupp, INRAE, and CABRICOOP breeders; Francois Pompanon, LECA. Iran: Hamid R. Rezaei, Gorgan University of Agricultural Sciences and Natural Resources. Ireland: Sean Carolan and Maeve Foran, Old Irish Goats Society. Italy: Alessandra Stella, IBBA- CNR; Paolo Ajmone- Marsan, Licia Colli, Marcello Del Corvo, Universita Cattolica del Sacro Cuore; Alessandra Crisa, Council for Agricultural Research and Economics (CREA); Donata Marletta, University of Catania, Italian Goat Consortium; Paola Crepaldi, University of Milan, Italian Goat Consortium; Michele Ottino, Parco Nazionale del Gran Paradiso; Ettore Randi, ISPRA Istituto Superiore per la Protezione e la Ricerca Ambientale. Kenya: Denis Fidalis Mujibi, International Livestock Research Institute (ILRI), Nairobi; Malawi: Timothy Gondwe, Department of Animal Science, Lilongwe University of Agriculture and Natural Resources. Morocco: Badr Benjelloun, INRA Maroc. Mozambique: Maria Da Gloria Taela, Animal Production Institute, Ministry of Agriculture, Maputo. Netherlands: Johannes A, Lenstra, Utrecht University. Nigeria: Oyekan Nash, National Biotechnology Development Agency. Pakistan: Muhammad Moaeen- ud- Din, PMAS- Arid Agriculture University. South Africa: Carina Visser, Faculty of Natural and Agricultural Sciences. Spain: Felix Goyache, Isabel Alvarez, Área de Genética y Reproducción Animal del Serida, Villaviciosa; Marcel Amills and Armand Sànchez, Centre for Research in Agricultural Genomics (CRAG); Juan Capote, Instituto Canario de Investigaciones Agrarias (ICIA); Jordi Jordana, Universitat Autònoma de Barcelona; Agueda Pons, Serveis de Millora Agrària i Pesquera (SEMILLA), Illes Balears; Amparo Martínez and Antonio Molina, University of Córdoba. Switzerland: Cord Drögemüller, University of Bern. Tanzania: Hassan Ally Mruttu, Ministry of Livestock and Fisheries Development. Uganda: Clet Wandui Masiga, Tropical Institute of Development Innovations (TRIDI); UK: Emily Clark, Mazdak Salavati, University of Edinburgh, Edinburgh. USA: Benjamin Rosen, Curtis P. Van Tassell, USDA/ARS; Jim Reecy, Iowa State University, Ames; Gordon Luikart, Montana Conservation Genomics lab (MCGL), Division of Biological Sciences, University of Montana, Missoula. Zimbabwe: Joseph Sikosana, Department of Research and Specialist Services, Division of Livestock Research. ECONOG ENE CONSORTIUM Albania: Hoda Anila, Dobi Petrit, Fac. Agriculture, Tirana. Austria: Baumung Roswitha, Univ. Natural Resources Applied Life Sciences, Vienna. Belgium: Baret Philippe, Fadlaoui Aziz, Univ. Cath. Louvain, Louvain- la- Neuve. Cyprus: Papachristoforou Christos, Agricult. Research Instit., Nicosia. Egypt: El- Barody M.A.A., Minia Univ. France: Taberlet Pierre, England Phillip, Luikart Gordon, Beja- Pereira Albano, Zundel Stéphanie, Univ. Joseph Fourier, Grenoble; Trommetter Michel, Inst. Recherche Agronomique, Grenoble. Germany: Erhardt Georg, Brandt Horst, Ibeagha- Awemu Eveline, Lühken, Gesine, Daniela Krugmann, Prinzenberg Eva- Maria, Lipsky Shirin, Gutscher Katja, Peter Christina, Justus- Liebig Univ. Giessen; Roosen Jutta, Bertaglia Marco, Univ. Kiel. Greece: Georgoudis Andreas, Al Tarrayrah, Jamil, Kliambas Georgios, Kutita Olga, Karetsou Katerina, Aristotle Univ. Thessaloniki, Thessaloniki; Ligda Christina, National Agric. Research Foundation, Thessaloniki. Hungary: Istvan Anton, Fesus Lazlo, Research Inst. Animal Breeding Nutrition, Herceghalom. Italy: Ajmone- Marsan Paolo, Canali Gabriele, Milanesi Elisabetta, Pellecchia Marco, Università Cattolica S. Cuore, Piacenza; Carta Antonello, Sechi Tiziana, Ist. Zootecnico Caseario Sardegna, Sassari; Cicogna Mario, Fornarelli Francesca, Giovenzana Stefano, Marilli Marta, Univ. Studi di Milano; Marletta Donata, Bordonaro S., D'Urso Giuseppe, Univ. Studi Catania; Pilla Fabio, D'Andrea Mariasilvia, Univ. Molise, Campobasso; Valentini Alessio, Cappuccio Irene, Pariset Lorraine, Univ. Tuscia, Viterbo. Jordan: Abo- Shehada Mahamoud, Jordan Univ. Science Technology. Netherlands: Lenstra Johannes A., Nijman, Isaäc J., Van Cann, Lisette M., Utrecht Univ. Poland: Niznikowski Roman, Dominik Popielarczyk, Strzelec Ewa, Warsaw Agricultural Univ. Romania: Vlaic Augustin, Univ. Cluj- Napoca. Spain: Dunner Susana, Canon Javier, Cortes Oscar, Garcia David, Univ. Complutense Madrid. Switzerland: Caloz Régis, EPFL, Lausanne; Obexer- Ruff Gabriela, Marie- Louise Glowatzki, Universität Bern. Turkey: Ertugrul Okan, Ankara Univ.; Togan Inci, Koban Evren, Middle East Technical Univ., Ankara. UK: Bruford Mike, Perez Trinidad, Juma Gabriela, Cardiff Univ.; Hewitt Godfrey, Dalamitra Stella, Wiskin Louise, Taylor Martin, Univ. East Anglia, Norwich; Jones Sam, The Sheep Trust; Scarpa Riccardo, Univ. York. DATA AVAIL ABILIT Y S TATEMENT Dideoxy sequences of DBY, ZFY and SRY segments representing dif- ferent haplogroups are accessible via the NCBI codes MF741774.1 to MF74182.1 and MG545047.1 to MG545050.1. The preliminary data set or 5356 genotypes in 70 mainly Asian and Moroccan male goats in plink format can be accessed via https://osf.io/ngx7u/ ?view_ only=64b4f 8bf78 794cd db5f4 7a3e3 e3d4534. For the VarGoats data set, see www.goatg enome.org/vargo ats.html and Denoyelle et al. (2021); see PRJEB37507 for the fastq files; vcf files are avail- able upon reasonable request. ORCID Guangxin E. https://orcid.org/0000-0001-7116-5177 Filipe Pereira https://orcid.org/0000-0001-8950-1036 Anne Da Silva https://orcid.org/0000-0001-7854-6642 Licia Colli https://orcid.org/0000-0002-7221-2905 Gwenola Tosser- Klopp https://orcid.org/0000-0003-0550-4673 Johannes A. Lenstra https://orcid.org/0000-0003-1745-7739     |  4377VarGoats CoNsortIUM et al. R E FE R E N C E S Ajmone- Marsan, P., Colli, L., Han, J. L., Achilli, A., Lancioni, H., Joost, S., Crepaldi, P., Pilla, F., Stella, A., Taberlet, P., Boettcher, P., Negrini, R., & Lenstra, J. A. (2014). The characterization of goat genetic diversity: Towards a genomic approach. Small Ruminant Research, 121(1), 58– 72. https://doi.org/10.1016/j.small rumres.2014.06.010 Alberto, F. J., Boyer, F., Orozco- Terwengel, P., Streeter, I., Servin, B., De Villemereuil, P., Benjelloun, B., Librado, P., Biscarini, F., Colli, L., Barbato, M., Zamani, W., Alberti, A., Engelen, S., Stella, A., Joost, S., Ajmone- Marsan, P., Negrini, R., Orlando, L., … Pompanon, F. (2018). Convergent genomic signatures of domestication in sheep and goats. Nature Communications, 9(1), 813– 821. https://doi. org/10.1038/s41467-018-03206-y Amills, M., Capote, J., & Tosser- Klopp, G. (2017). Goat domestication and breeding: A jigsaw of historical, biological and molecular data with missing pieces. Animal Genetics, 48(6), 631– 644. https://doi. org/10.1111/age.12598 Banda, J. W., Ayoude, J. A., Karua, S. K., & Kamwanja, L. A. (1993). The local Malawi goat. World Animal Review, 74(75), 49– 75. Bandelt, H. J., Forster, P., & Röhl, A. (1999). Median- joining networks for inferring intraspecific phylogenies. Molecular Biology and Evolution, 16(1), 37– 48. https://doi.org/10.1093/oxfor djour nals.molbev. a026036 Batini, C., & Jobling, M. A. (2017). Detecting past male- mediated expan- sions using the Y chromosome. Human Genetics, 136(5), 547– 557. https://doi.org/10.1007/s0043 9- 017- 1781- z Belabdi, I., Ouhrouch, A., Lafri, M., Gaouar, S. B. S., Ciani, E., Benali, A. R., Benali, A. R., Ould Ouelhadj, H., Haddioui, A., Pompanon, F., Blanquet, V., Taurisson- Mouret, D., Harkat, S., Lenstra, J. A., Benjelloun, B., & Da Silva, A. (2019). Genetic homogenization of in- digenous sheep breeds in Northwest Africa. Scientific Reports, 9(1), 7920. https://doi.org/10.1038/s41598-019-44137-y Cai, Y., Fu, W., Cai, D., Heller, R., Zheng, Z., Wen, J., Li, H., Wang, X., Alshawi, A., Sun, Z., Zhu, S., Wang, J., Yang, M., Hu, S., Li, Y., Yang, Z., Gong, M., Hou, Y., Lan, T., … Russo, C. (2020). Ancient genomes reveal the evolutionary history and origin of cashmere- producing goats in China. Molecular Biology and Evolution, 37(7), 2099– 2109. https://doi.org/10.1093/molbe v/msaa103 Chen, N., Cai, Y., Chen, Q., Li, R., Wang, K., Huang, Y., Hu, S., Huang, S., Zhang, H., Zheng, Z., Song, W., Ma, Z., Ma, Y., Dang, R., Zhang, Z., Xu, L., Jia, Y., Liu, S., Yue, X., … Lei, C. (2018). Whole- genome resequencing reveals world- wide ancestry and adap- tive introgression events of domesticated cattle in East Asia. Nature Communications, 9(1), 1– 13. https://doi.org/10.1038/ s41467-018-04737-0 Choi, S. K., Kim, K. S., Ranyuk, M., Babaev, E., Voloshina, I., Bayarlkhagva, D., Chong, J. R., Ishiguro, N., Yu, L., Min, M. S., Lee, H., & Markov, N. (2020). Asia- wide phylogeography of wild boar (Sus scrofa) based on mitochondrial DNA and Y- chromosome: Revising the migra- tion routes of wild boar in Asia. PLoS ONE, 15, 1– 17. https://doi. org/10.1371/journ al.pone.0238049 Ciani, E., Mastrangelo, S., Da Silva, A., Marroni, F., Ferenčaković, M., Ajmone- Marsan, P., Baird, H., Barbato, M., Colli, L., Delvento, C., Dovenski, T., Gorjanc, G., Hall, S. J. G., Hoda, A., Li, M.- H., Markovic, B., McEwan, J., Moradi, M. H., Ruiz- Larrañaga, O., … Lenstra, J. A. (2020). On the origin of European sheep as revealed by the diver- sity of the Balkan breeds and by optimizing population- genetic analysis tools. Genetics Selection Evolution, 52(1), 1– 14. https://doi. org/10.1186/s12711-020-00545-7 Çinar Kul, B., Bilgen, N., Lenstra, J. A., Korkmaz Agaoglu, O., Akyuz, B., & Ertugrul, O. (2015). Y- chromosomal variation of local goat breeds of Turkey close to the domestication Centre. Journal of Animal Breeding and Genetics, 132(6), 449– 453. https://doi.org/10.1111/ jbg.12154 Colli, L., Milanesi, M., Talenti, A., Bertolini, F., Chen, M., Crisà, A., Daly, K., Del Corvo, M., Guldbrandtsen, B., Lenstra, J. A., Rosen, B. D., Vajana, E., Catillo, G., Joost, S., Nicolazzi, E. L., Rochat, E., Rothschild, M. F., Servin, B., Sonstegard, T. S., … Stella, A. (2018). Genome- wide SNP profiling of worldwide goat populations reveals strong partitioning of diversity and highlights post- domestication migration routes. Genetics Selection Evolution, 50(1), 1– 20. https:// doi.org/10.1186/s12711-018-0422-x Cymbron, T., Freeman, A. R., Malheiro, M. I., Vigne, J. D., & Bradley, D. G. (2005). Microsatellite diversity suggests different histories for Mediterranean and northern European cattle populations. Proceedings of the Royal Society B: Biological Sciences, 272(1574), 1837– 1843. https://doi.org/10.1098/rspb.2005.3138 Daly, K. G., Delser, P. M., Mullin, V. E., Scheu, A., Mattiangeli, V., Teasdale, M. D., Hare, A. J., Burger, J., Verdugo, M. P., Collins, M. J., Kehati, R., Erek, C. M., Bar- Oz, G., Pompanon, F., Cumer, T., Çakırlar, C., Mohaseb, A. F., Decruyenaere, D., Davoudi, H., … Bradley, D. G. (2018). Ancient goat genomes reveal mosaic domestication in the Fertile Crescent. Science, 361(6397), 85– 88. https://doi. org/10.1126/scien ce.aas9411 Daly, K. G., Mattiangeli, V., Hare, A. J., Davoudi, H., Doost, S. B., Fathi, H., Amiri, S., Khazaeli, R., Decruyenaere, D., Nokandeh, J., Richter, T., Darabi, H., Mortensen, P., Pantos, A., Yeomans, L., Bangsgaard, P., Mashkour, M., Zeder, M. A., & Bradley, D. G. (2021). Herded and hunted goat genomes from the dawn of domestication in the Zagros Mountains. Proceedings of the National Academy of Sciences of the United States of America, 118(25), e2100901118. https://doi. org/10.1073/pnas.21009 01118 Deng, J., Xie, X. L., Wang, D. F., Zhao, C., Lv, F. H., Li, X., Yang, J., Yu, J.- L., Shen, M., Gao, L., Yang, J.- Q., Liu, M.- J., Li, W.- R., Wang, Y.- T., Wang, F., Li, J.- Q., Hehua, E., Liu, Y.- G., Shen, Z.- Q., … Li, M. H. (2020). Paternal origins and migratory episodes of domestic sheep. Current Biology, 30(20), 4085– 4095.e6. https://doi.org/10.1016/j. cub.2020.07.077 Deniskova, T. E., Dotsev, A. V., Selionova, M. I., Reyer, H., Sölkner, J., Fornara, M. S., Aybazov, A. M., Wimmers, K., Brem, G., & Zinovieva, N. A. (2021). SNP- based genotyping provides insight into the west Asian origin of Russian local goats. Frontiers in Genetics, 12, 1– 16. https://doi.org/10.3389/fgene.2021.708740 Denoyelle, L., Talouarn, E., Bardou, P., Colli, L., Alberti, A., Danchin, C., Del Corvo, M., Engelen, S., Orvain, C., Palhière, I., Rupp, R., Sarry, J., Salavati, M., Amills, M., Clark, E., Crepaldi, P., Faraut, T., Masiga, C. W., Pompanon, F., … Nash, O. (2021). VarGoats project: a dataset of 1159 whole- genome sequences to dissect Capra hircus global diversity. Genetics Selection Evolution, 53(1), 1– 14. https://doi. org/10.1186/s12711-021-00659-6 Dubeuf, J.- P., Morand- Fehr, P., & Rubino, R. (2004). Situation, changes and future of goat industry around the world. Small Ruminant Research, 51(2), 165– 173. https://doi.org/10.1016/j.small rumres.2003.08.007 Edwards, C. J., Ginja, C., Kantanen, J., Pérez- Pardal, L., Tresset, A., Stock, F., European Cattle Genetic Diversity Consortium, Gama, L. T., Penedo, M. C., Bradley, D. G., Lenstra, J. A., & Nijman, I. J. (2011). Dual origins of dairy cattle farming - evidence from a comprehensive survey of european Y- chromosomal varia- tion. PLoS ONE, 6(1), e15922. https://doi.org/10.1371/journ al.pone.0015922 Felkel, S., Vogl, C., Rigler, D., Dobretsberger, V., Chowdhary, B. P., Distl, O., Fries, R., Jagannathan, V., Janečka, J. E., Leeb, T., Lindgren, G., McCue, M., Metzger, J., Neuditschko, M., Rattei, T., Raudsepp, T., Rieder, S., Rubin, C. J., Schaefer, R., … Wallner, B. (2019a). The horse Y chromo- some as an informative marker for tracing sire lines. Scientific Reports, 9(1), 1– 12. https://doi.org/10.1038/s41598-019-42640-w Felkel, S., Wallner, B., Chuluunbat, B., Yadamsuren, A., Faye, B., Brem, G., Walzer, C. & Burger, P. A. (2019b). A first Y- chromosomal haplotype 4378  |    VarGoats CoNsortIUM et al. network to investigate Male- driven population dynamics in do- mestic and wild bactrian camels. Frontiers in Genetics, 10, 1– 7. doi: https://doi.org/10.3389/fgene.2019.00423 Ganguly, I., Jeevan, C., Singh, S., Dixit, S. P., Sodhi, M., Ranjan, A., Kumar, S., & Sharma, A. (2020). Y- chromosome genetic diversity of Bos indicus cattle in close proximity to the centre of domesti- cation. Scientific Reports, 10(1), 9992. https://doi.org/10.1038/ s41598-020-66133-3 Gao, J., Lyu, Y., Zhang, D., Reddi, K. K., Sun, F., Yi, J., Liu, C., Li, H., Yao, H., Dai, J., & Xu, F. (2020). Genomic characteristics and selection signatures in indigenous chongming white goat (Capra hircus). Frontiers in Genetics, 11(August), 1– 14. https://doi.org/10.3389/ fgene.2020.00901 Garrine, C. M. L. P. (2007). Genetic characterization of indigenous goat pop- ulations of Mozambique. University of Pretoria. Grossen, C., Guillaume, F., Keller, L. F., & Croll, D. (2020). Purging of highly deleterious mutations through severe bottlenecks in alpine ibex. Nature Communications, 11(1), 1001- 1012. https://doi.org/10.1038/ s4146 7- 020- 14803 - 1 Guirao- Rico, S., Ramirez, O., Ojeda, A., Amills, M., & Ramos- Onsins, S. E. (2018). Porcine Y- chromosome variation is consistent with the occurrence of paternal gene flow from non- Asian to Asian popu- lations. Heredity, 120(1), 63– 76. https://doi.org/10.1038/s4143 7- 017- 0002- 9 Hegde, N. G. (2019). Livestock development for sustainable livelihood of small farmers. Asian Journal of Research in Animal and Veterinary Sciences, 3(2), 1– 17. Hughes, J. F., Skaletsky, H., Koutseva, N., Pyntikova, T., & Page, D. C. (2015). Sex chromosome- to- autosome transposition events counter Y- chromosome gene loss in mammals. Genome Biology, 16(1), 104. https://doi.org/10.1186/s1305 9- 015- 0667- 4 Huson, D. H., & Bryant, D. (2006). Application of phylogenetic networks in evolutionary studies. Molecular Biology and Evolution, 23(2), 254– 267. https://doi.org/10.1093/molbe v/msj030 Jobling, M. A., & Tyler- Smith, C. (2017). Human Y- chromosome varia- tion in the genome- sequencing era. Nature Reviews Genetics, 18(8), 485– 497. https://doi.org/10.1038/nrg.2017.36 Kijas, J. W., Lenstra, J. A., Hayes, B., Boitard, S., Neto, L. R., Cristobal, M. S., Servin, B., McCulloch, R., Whan, V., Gietzen, K., Paiva, S., Barendse, W., Ciani, E., Raadsma, H., McEwan, J., Dalrymple, B., & The International Sheep Genomics Consortium. (2012). Genome- wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biology, 10(2), e1001258. https://doi.org/10.1371/journ al.pbio. 1001258 Kim, J. Y., Jeong, S., Kim, K. H., Lim, W. J., Lee, H. Y., & Kim, N. (2019). Discovery of genomic characteristics and selection signatures in Korean indigenous goats through comparison of 10 goat breeds. Frontiers in Genetics, 10(Jul), 1– 17. https://doi.org/10.3389/ fgene.2019.00699 Kivisild, T. (2017). The study of human Y chromosome variation through ancient DNA. Human Genetics, 136(5), 529– 546. https://doi. org/10.1007/s0043 9- 017- 1773- z Larson, G., & Fuller, D. Q. (2014). The evolution of animal domestication. Annual Review of Ecology, Evolution, and Systematics, 45, 115– 136. https://doi.org/10.1146/annur ev- ecols ys- 11051 2- 135813 Leigh, J. W., & Bryant, D. (2015). POPART: Full- feature software for hap- lotype network construction. Methods in Ecology and Evolution, 6(9), 1110– 1116. https://doi.org/10.1111/2041- 210X.12410 Lenstra, J. A., & Econogene Consortium. (2005). Evolutionary and de- mographic history of sheep and goats suggested by nuclear, mtDNA and Y- chromosomal markers. The Role of Biotechnology for the Characterization of Crop, Forestry, Animal and Fishery Genetic Resources. Turin. Li, R., Yang, P., Dai, X., Asadollahpour Nanaei, H., Fang, W., Yang, Z., Cai, Y., Zheng, Z., Wang, X., & Jiang, Y. (2021). A near complete genome for goat genetic and genomic research. Genetics Selection Evolution, 53(1), 1– 17. https://doi.org/10.1186/s12711-021-00668-5 Lu, C. D. (2011). Boer goat production: Progress and perspective. Proceedings of International Conference on Boer Goats, 1– 11. Anshun, China. Luikart, G., Gielly, L., Excoffier, L., Vigne, J., Bouvet, J., & Taberlet, P. (2001). Multiple maternal origins and weak phylogeographic struc- ture in domestic goats. Proceedings of the National Academy of Sciences of the United States of America, 98(10), 5927– 5932. https:// doi.org/10.1073/pnas.09159 1198 Lv, F. H., Agha, S., Kantanen, J., Colli, L., Stucki, S., Kijas, J. W., Joost, S., Li, M.- H., & Marsan, P. A. (2014). Adaptations to climate- mediated selective pressures in sheep. Molecular Biology and Evolution, 31(12), 3324– 3343. https://doi.org/10.1093/molbe v/msu264 Magnier, J., Druet, T., Naves, M., Ouvrard, M., Raoul, S., Janelle, J., Moazami- Goudarzi, K., Lesnoff, M., Tillard, E., Gautier, M., & Flori, L. (2022). The genetic history of Mayotte and Madagascar cat- tle breeds mirrors the complex pattern of human exchanges in Western Indian Ocean. G3 (Bethesda), 12(4), jkac029. https://doi. org/10.1093/g3jou rnal/jkac029 Marín, J. C., Romero, K., Rivera, R., Johnson, W. E., & González, B. A. (2017). Y- chromosome and mtDNA variation confirms independent domestications and directional hybridization in south American camelids. Animal Genetics, 48(5), 591– 595. https://doi.org/10.1111/ age.12570 Matsunaga, M., Ohtsubo, Y., Masuda, T., Noguchi, Y., Yamasue, H., & Ishii, K. (2021). A genetic variation in the Y chromosome among modern Japanese males related to several physiological and psychological characteristics. Frontiers in Behavioral Neuroscience, 15(December), 1– 9. https://doi.org/10.3389/fnbeh.2021.774879 Meadows, J. R. S., & Kijas, J. W. (2009). Re- sequencing regions of the ovine Y chromosome in domestic and wild sheep reveals novel paternal haplotypes. Animal Genetics, 40(1), 119– 123. https://doi. org/10.1111/j.1365- 2052.2008.01799.x Miller, B. A., & Lu, C. D. (2019). Special issue — Current status of global dairy goat production: An overview. Asian- Australasian Journal of Animal Sciences, 32(8), 1219– 1232. https://doi.org/10.5713/ ajas.19.0253 NAADS. (2005). Goat Production Manual. (January). Naderi, S., Rezaei, H. R., Pompanon, F., Blum, M. G. B., Negrini, R., Naghash, H. R., Balkiz, O., Mashkour, M., Gaggiotti, O. E., Ajmone- Marsan, P., Kenc, A., Vigne, J. D., & Taberlet, P. (2008). The goat do- mestication process inferred from large- scale mitochondrial DNA analysis of wild and domestic individuals. Proceedings of the National Academy of Sciences of the United States of America, 105(46), 17659– 17664. https://doi.org/10.1073/pnas.08047 82105 Naderi, S., Rezaei, H. R., Taberlet, P., Zundel, S., Rafat, S. A., Naghash, H. R., el- Barody, M.A., Ertugrul, O., & Zundel, S. (2007). Large- scale mitochondrial DNA analysis of the domestic goat reveals six hap- logroups with high diversity. PLoS ONE, 2(10), e1012. https://doi. org/10.1371/journ al.pone.0001012 Natanaelsson, C., Oskarsson, M. C. R., Angleby, H., Lundeberg, J., Kirkness, E., & Savolainen, P. (2006). Dog Y chromosomal DNA sequence: Identification, sequencing and SNP discovery. BMC Genetics, 7, 45- 50. https://doi.org/10.1186/1471- 2156- 7- 45 Nijman, I. J., Van Boxtel, D. C. J., Van Cann, L. M., Marnoch, Y., Cuppen, E., & Lenstra, J. A. (2008). Phylogeny of Y chromo- somes from bovine species. Cladistics, 24(5), 723– 726. https://doi. org/10.1111/j.1096- 0031.2008.00201.x O’Connell, J., Gurdasani, D., Delaneau, O., Pirastu, N., Ulivi, S., Cocca, M., Traglia, M., Huang, J., Huffman, J. E., Rudan, I., McQuillan, R., Fraser, R. M., Campbell, H., Polasek, O., Asiki, G., Ekoru, K., Hayward, C., Wright, A. F., Vitart, V., … Marchini, J. (2014). A general approach for haplotype phasing across the full spectrum of relatedness. PLoS Genetics, 10(4), e1004234. https://doi.org/10.1371/journ al.pgen.1004234     |  4379VarGoats CoNsortIUM et al. Oetjens, M. T., Martin, A., Veeramah, K. R., & Kidd, J. M. (2018). Analysis of the canid Y- chromosome phylogeny using short- read sequencing data reveals the presence of distinct haplogroups among Neolithic European dogs. BMC Genomics, 19(1), 1– 9. https://doi.org/10.1186/ s1286 4- 018- 4749- z Onzima, R. B., Upadhyay, M. R., Mukiibi, R., Kanis, E., Groenen, M. A. M., & Crooijmans, R. P. M. A. (2018). Genome- wide population struc- ture and admixture analysis reveals weak differentiation among Ugandan goat breeds. Animal Genetics, 49(1), 59– 70. https://doi. org/10.1111/age.12631 Peacock, C. (2005). Goats – a pathway out of poverty. Small Ruminant Research, 60(1– 2 Spec. Iss), 179– 186. https://doi.org/10.1016/j. small rumres.2005.06.011 Pereira, F., Queirós, S., Gusmão, L., Nijman, I. J., Cuppen, E., Lenstra, J. A., Econogene Consortium, Davis, S. J., Nejmeddine, F., & Amorim, A. (2009). Tracing the history of goat pastoralism: New clues from mitochondrial and y chromosome DNA in North Africa. Molecular Biology and Evolution, 26(12), 2765– 2773. https://doi.org/10.1093/ molbe v/msp200 Pidancier, N., Jordan, S., Luikart, G., & Taberlet, P. (2006). Evolutionary history of the genus Capra (Mammalia, Artiodactyla): Discordance between mitochondrial DNA and Y- chromosome phylogenies. Molecular Phylogenetics and Evolution, 40(3), 739– 749. https://doi. org/10.1016/j.ympev.2006.04.002 Porter, V., Alderson, L., Hall, S. J. G., & Sponenberg, D. P. (2016). Mason's world encyclopedia of livestock breeds and breeding. CABI. Rivollat, M., Mendisco, F., Pemonge, M. H., Safi, A., Saint- Marc, D., Brémond, A., Couture- Veschambre, C., Rottier, S., & Deguilloux, M. F. (2015). When the waves of European neolithization met: First paleogenetic evidence from early farmers in the Southern Paris Basin. PLoS ONE, 10(4), 1– 16. https://doi.org/10.1371/journ al.pone.0125521 Somenzi, E., Senczuk, G., Ciampolini, R., Cortellari, M., Vajana, E., Tosser- Klopp, G., Pilla, F., Ajmone- Marsan, P., Crepaldi, P., & Colli, L. (2022). The SNP- based profiling of montecristo feral goat populations re- veals a history of isolation, bottlenecks, and the effects of man- agement. Genes, 13(2), 213– 226. https://doi.org/10.3390/genes 13020213 Stiner, M. C., Munro, N. D., Buitenhuis, H., Duru, G., & Özbaşaran, M. (2022). An endemic pathway to sheep and goat domestica- tion at Aşıklı Höyük (Central Anatolia, Turkey). Proceedings of the National Academy of Sciences, 119(4), e2110930119. https://doi. org/10.1073/pnas.21109 30119 Tabata, R., Kawaguchi, F., Sasazaki, S., Yamamoto, Y., Bakhtin, M., Kazymbet, P., Meldevekob, A., Suleimenov, M. Z., Nishibori, M., & Mannen, H. (2018). The Eurasian Steppe is an important goat propagation route: A phylogeographic analysis using mitochondrial DNA and Y- chromosome sequences of Kazakhstani goats. Animal Science Journal, 90(3), 317– 322. https://doi.org/10.1111/asj.13144 Tabata, R., Kawaguchi, F., Sasazaki, S., Yamamoto, Y., Rakotondraparany, F., Ratsoavina, F. M., Yonezawa, T., & Mannen, H. (2019). Phylogeographic analysis of Madagascan goats using mtDNA con- trol region and SRY gene sequences. Zoological Science, 36, 294– 298. https://doi.org/10.2108/zs180184 Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., & Kumar, S. (2011). MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsi- mony methods. Molecular Biology and Evolution, 28(10), 2731– 2739. https://doi.org/10.1093/molbe v/msr121 Tresset, A., & Vigne, J.- D. (2007). Substitution of species, techniques and symbols at the Mesolithic- Neolithic transition in Western Europe. In Going over: The Mesolithic- Neolithic transition in North- West Europe. British Academy. https://doi.org/10.5871/bacad/ 97801 97264 140.003.0010 Utaaker, K. S., Chaudhary, S., Kifleyohannes, T., & Robertson, L. J. (2021). Global goat! Is the expanding goat population an important reservoir of cryptosporidium? Frontiers in Veterinary Science, 8, 648500. https://doi.org/10.3389/fvets.2021.648500 Vidal, O., Drögemüller, C., Obexer- Ruff, G., Reber, I., Jordana, J., Martínez, A., Bâlteanu, V. A., Delgado, J. V., Eghbalsaied, S., Landi, V., Goyache, F., Traoré, A., Pazzola, M., Vacca, G. M., Badaoui, B., Pilla, F., D'Andrea, M., Álvarez, I., Capote, J., … Amills, M. (2017). Differential distribution of Y- chromosome haplotypes in Swiss and Southern European goat breeds. Scientific Reports, 7(1), 16161. https://doi.org/10.1038/s41598-017-15593-1 Waki, A., Sasazaki, S., Kobayashi, E., & Mannen, H. (2015). Paternal phylogeography and genetic diversity of East Asian goats. Animal Genetics, 46(3), 337– 339. https://doi.org/10.1111/age.12293 Wallner, B., Palmieri, N., Vogl, C., Rigler, D., Bozlak, E., Druml, T., Jagannathan, V., Leeb, T., Fries, R., Tetens, J., Thaller, G., Metzger, J., Distl, O., Lindgren, G., Rubin, C. J., Andersson, L., Schaefer, R., McCue, M., Neuditschko, M., … Brem, G. (2017). Y chromo- some uncovers the recent oriental origin of modern stallions. Current Biology, 27(13), 2029– 2035.e5. https://doi.org/10.1016/j. cub.2017.05.086 Wang, H., Zhang, L., Cao, J., Wu, M., Ma, X., Liu, Z., Liu, R., Zhao, F., Wei, C., & Du, L. (2015). Genome- wide specific selection in three domes- tic sheep breeds. PLoS ONE, 10(6), 1– 21. https://doi.org/10.1371/ journ al.pone.0128688 Wutke, S., Sandoval- Castellanos, E., Benecke, N., Döhle, H. J., Friederich, S., Gonzalez, J., Hofreiter, M., Lõugas, L., Magnell, O., Malaspinas, A. S., Morales- Muñiz, A., Orlando, L., Reissmann, M., Trinks, A., & Ludwig, A. (2018). Decline of genetic diversity in ancient domestic stallions in Europe. Science Advances, 4(4), eaap9691. https://doi. org/10.1126/sciadv.aap9691 Xia, X., Yao, Y., Li, C., Zhang, F., Qu, K., Chen, H., Huang, B., & Lei, C. (2019). Genetic diversity of Chinese cattle revealed by Y- SNP and Y- STR markers. Animal Genetics, 50(1), 64– 69. https://doi. org/10.1111/age.12742 Xiao, C., Li, J., Xie, T., Chen, J., Zhang, S., Elaksher, S. H., Jiang, F., Jiang, Y., Zhang, L., Zhang, W., Xiang, Y., Wu, Z., Zhao, S., & Du, X. (2021). The assembly of caprine Y chromosome sequence reveals a unique pa- ternal phylogenetic pattern and improves our understanding of the origin of domestic goat. Ecology and Evolution, 11(12), 7779– 7795. https://doi.org/10.1002/ece3.7611 Yang, C., Ba, H., Zhang, W., Zhang, S., Zhao, H., Yu, H., Gao, Z., & Wang, B. (2018). The association of 22 Y chromosome short tandem re- peat loci with initiative– aggressive behavior. Gene, 654(November 2017), 10– 13. https://doi.org/10.1016/j.gene.2018.02.039 Zhang, K., Lenstra, J. A., Zhang, S., Liu, W., & Liu, J. (2020). Evolution and domestication of the Bovini species. Animal Genetics, 51(5), 637– 657. https://doi.org/10.1111/age.12974 Zhang, M., Peng, W. F., Yang, G. L., Lv, F. H., Liu, M. J., Li, W. R., Liu, Y. G., Li, J. Q., Wang, F., Shen, Z. Q., Zhao, S. G., Hehua, E., Marzanov, N., Murawski, M., Kantanen, J., & Li, M. H. (2014). Y chromosome haplotype diversity of domestic sheep (Ovis aries) in northern Eurasia. Animal Genetics, 45(6), 903– 907. https://doi.org/10.1111/ age.12214 Zhang, X., Liao, Z., Tang, S., Yuan, Z., Li, F., & Yue, X. (2021). A muta- tion modulating DDX3Y gene expression cosegregates with the major Y- chromosomal haplogroups and with testis size in hu sheep. Animal Genetics, 53, 193– 202. https://doi.org/10.1111/age.13166 Zhang, Y., Lu, Y., Yindee, M., Li, K. Y., Kuo, H. Y., Ju- Ten, J., Ye, S., Faruque, M. O., Li, Q., Wang, Y., Cuong, V. C., Pham, L. D., Bouahom, B., Yang, B., Liang, X., Cai, Z., Vankan, D., Manatchaiworakul, W., Kowlim, N., … Barker, J. S. F. (2016). Strong and stable geographic differentiation of swamp buffalo maternal and paternal lineages indicates domestication in the China/Indochina border region. Molecular Ecology, 25(7), 1530– 1550. https://doi.org/10.1111/ mec.13518 Zhao, W., Zhong, T., Wang, L. J., Li, L., & Zhang, H. P. (2014). Extensive female- mediated gene flow and low phylogeography among 4380  |    VarGoats CoNsortIUM et al. seventeen goat breeds in Southwest China. Biochemical Genetics, 52(7– 8), 355– 364. https://doi.org/10.1007/s1052 8- 014- 9652- y Zhao, Y., Zhao, R., Zhao, Z., Xu, H., Zhao, E., & Zhang, J. (2014). Genetic diversity and molecular phylogeography of Chinese domes- tic goats by large- scale mitochondrial DNA analysis. Molecular Biology Reports, 41(6), 3695– 3704. https://doi.org/10.1007/s1103 3- 014- 3234- 2 Zheng, Z., Wang, X., Li, M., Li, Y., Yang, Z., Wang, X., Wang, X., Pan, X., Gong, M., Zhang, Y., Guo, Y., Wang, Y., Liu, J., Cai, Y., Chen, Q., Okpeku, M., Colli, L., Cai, D., … Jiang, Y. (2020). The origin of domestication genes in goats. Science Advances, 6(21). eaaz5216 https://doi.org/10.1126/sciadv.aaz5216 SUPPORTING INFORMATION Additional supporting information can be found online in the Supporting Information section at the end of this article. How to cite this article: VarGoats Consortium, Nijman, I. J., Rosen, B. D., Bardou, P., Faraut, T., Cumer, T., Daly, K. G., Zheng, Z., Cai, Y., Asadollahpour, H., Kul, B. Ç., Zhang, W-Y, E., G., Ayin, A., Baird, H., Bakhtin, M., Bâlteanu, V. A., Barfield, D., Berger, B. … Lenstra, J. A. (2022). Geographical contrasts of Y- chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions. Molecular Ecology, 31, 4364–4380. https://doi.org/10.1111/mec.16579