PLoS ONE
Public Library of Science
image
High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor
DOI 10.1371/journal.pone.0233333 , Volume: 15 , Issue: 5

Table of Contents

Highlights

Notes

Abstract

Goats were initially managed in the Near East approximately 10,000 years ago and spread across Eurasia as economically productive and environmentally resilient herd animals. While the geographic origins of domesticated goats (Capra hircus) in the Near East have been long-established in the zooarchaeological record and, more recently, further revealed in ancient genomes, the precise pathways by which goats spread across Asia during the early Bronze Age (ca. 3000 to 2500 cal BC) and later remain unclear. We analyzed sequences of hypervariable region 1 and cytochrome b gene in the mitochondrial genome (mtDNA) of goats from archaeological sites along two proposed transmission pathways as well as geographically intermediary sites. Unexpectedly high genetic diversity was present in the Inner Asian Mountain Corridor (IAMC), indicated by mtDNA haplotypes representing common A lineages and rarer C and D lineages. High mtDNA diversity was also present in central Kazakhstan, while only mtDNA haplotypes of lineage A were observed from sites in the Northern Eurasian Steppe (NES). These findings suggest that herding communities living in montane ecosystems were drawing from genetically diverse goat populations, likely sourced from communities in the Iranian Plateau, that were sustained by repeated interaction and exchange. Notably, the mitochondrial genetic diversity associated with goats of the IAMC also extended into the semi-arid region of central Kazakhstan, while NES communities had goats reflecting an isolated founder population, possibly sourced via eastern Europe or the Caucasus region.

Hermes, Frachetti, Voyakin, Yerlomaeva, Beisenov, Doumani Dupuy, Papin, Motuzaite Matuzeviciute, Bayarsaikhan, Houle, Tishkin, Nebel, Krause-Kyora, Makarewicz, and Biehl: High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor

Introduction

Domesticated goats (Capra hircus ) are at the subsistence core of pastoralist and farming communities around the world, providing dependable sources of meat, milk, fiber and hides [1]. Goats are remarkably resilient in challenging environments, having a broader foraging spectrum and requiring less water than other livestock taxa [24]. Wild bezoar (Capra aegagrus ) goats were under initial management in various regions of the Near East by the mid-ninth millennium BC [5, 6]. By the eighth millennium BC, goats were closely controlled via husbandry practices that adjusted animal diets, movement, and herd demography favoring female survivorship [711], and also exhibited considerable size diminution and loss of wild phenotypic traits associated with the domestication syndrome, indicating that people had effectively reconfigured the fitness landscape of managed goats and genetically isolated them from sympatric bezoars [12, 13].

Recent paleogenetic research [14] and previous studies [1517] established that goat domestication processes involved multiple maternal origins throughout the Near East and Iranian Plateau, giving rise to early phylogeographic structure of six divergent mtDNA haplogroups (A, B, C, D, F, and G). However, post-Neolithic domesticated goat populations across these regions show a collapse of this initial mtDNA phylogeographic pattern, leading to widespread dominance of A lineages [14], which today characterize 91% of domesticated goats world-wide [1519]. Crucially, as Neolithic technologies and possibly people spread from centers of plant and animal domestication, mtDNA phylogeographic patterns may have been transmitted to adjacent and distant regions across Eurasia, thus, revealing goat transmission routes. Moreover, goat herds of later Bronze and Iron Age peoples might have been characterized by shifts in the relative abundance of mtDNA haplotypes on account of the geographic extent and intensity of social interaction spheres in which goats were bred and exchanged.

Central Asia is a critical region for examining the spread of domesticated plants and animals across the vast expanse of Eurasia, which led to widespread transformations in human subsistence and socio-political structures [20]. In Central Asia, regionally varied pastoralist lifeways developed in continental environments much different from the diverse landscapes of the Near East, where the earliest forms of mobile pastoralism appears to have first emerged in the late seventh-early sixth millennium BC, when specialized herders began residing in semi-permanent settlements located in the more marginal semi-arid eastern steppes of the Levant and possibly practiced residential mobility on a seasonal basis [2124]. There are two hypothesized routes for the transmission of subsistence herding into Central Asia and farther afield. The first is the “Northern Eurasian Steppe” (NES), a distinct ecotone of grasslands at roughly 50–55°N latitude from eastern Europe to southern Siberia (Fig 1) [25]. Archaeologists have long thought that diverse herding communities from the Pontic-Caspian steppe known as the “Yamnaya” engaged in a series of eastward migrations along the NES to the Altai Mountains and then south into Central Asia (ca. 3300–2500 BC), which coincides with pronounced changes in material culture [26], such as a distinct mortuary practice of stone enclosures and ceramic styles resembling that of contemporaneous Yamnaya sites [2729]. Critically, human archaeogenetic research documents close genetic links between these so-called “Afanasievo” and “Yamnaya” groups [3033]. In later periods, the NES arguably functioned as a cultural exchange thoroughfare over which key steppe innovations were transmitted [34], such as bronze technologies [35, 36] and domesticated horses [26, 37, 38]. By the Iron Age, the NES represents as a critical arena for the geographically expansive Iron Age “Scythian” cultural community (ca. 500 cal BC) and subsequent growth of “nomadic empires” [3944].

Map of Eurasia showing sites analyzed in this study in colorful shapes corresponding to analytical regions.
Fig 1
Post-Neolithic distribution of mitochondrial lineages across the Middle East follows Daly et al. [14], and supposed dispersal routes of domesticated goats are shown with dashed arrows. Made with Natural Earth. Free vector and raster map data @ http://www.naturalearthdata.com.Map of Eurasia showing sites analyzed in this study in colorful shapes corresponding to analytical regions.

A second route for the spread of pastoralist subsistence to Central Asia reflects a vector of cultural transmission from the Iranian Plateau via a near continuous series of foothill ecozones to the Altai Mountains, dubbed the “Inner Asian Mountain Corridor” (IAMC) (Fig 1) [20]. Accordingly, pastoralist subsistence (perhaps through population movement or cultural transmission) is thought to have spread into Central Asia during the fourth millennium BC along montane elevational gradients as a strategy to provide graze for livestock by exploiting seasonal rhythms of vegetation growth [45]. Recent archaeological research in the IAMC has documented pastoralist occupations dating to the early third millennium BC that coincided with the initial trans-Eurasian transmissions of domesticated cereals of wheat and barley from the Near East and broomcorn millet from northern China [4648]. Human paleogenetic research suggests that the IAMC in the second millennium BC was a distinct conduit for the northward human gene flow from agro-pastoralist communities from the region of Turan, roughly spanning present-day Turkmenistan and Uzbekistan [31]. In subsequent millennia, pastoralist interactions throughout the IAMC are associated with vibrant cultural exchanges, linking economies and cultures across Eurasia, which scholars commonly refer to as the ancient Silk Roads [45, 4954].

While the NES and IAMC served as major cultural interaction spheres over the past six thousand years, there has been very little systematic evaluation of the dispersals of pastoralist subsistence or animal exchanges within and between these geographic loci. Here, we examine the mitochondrial (mtDNA) diversity of domesticated goats from sites located in the NES and IAMC with Bronze Age and/or Iron Age occupations in order to reconstruct the spread of this key pastoralist animal to Central Asia and subsequent genetic patterns (Fig 1) [55]. Key Bronze Age sites located between the NES and IAMC provide a test for the geographic extent of maternally inherited goat genetic variation potentially characteristic of these macro regions (Fig 1). We analyzed the hypervariable region 1 (HVR1) of the control region and cytochrome b gene (MT-CYB) recovered from ancient goat skeletal remains. HVR1 sequences contain high frequencies of polymorphisms that characterize particular haplogroups, while MT-CYB sequences are highly conserved and generally reflect inter-species variation, which we use here to differentiate between domesticated sheep and goat for ambiguously identified archaeological specimens, in addition to wild Capra species [56], which are likely present in the zooarchaeological assemblages of the sampled sites. Prior to a major study across the Near East [14], a very limited number of ancient mtDNA sequences have been analyzed, which only represent Neolithic Europe [57] and the Near East [5860] and also Iron Age China [61, 62].

Materials and methods

Sample selection, preparation, and DNA extraction

All necessary permits were obtained for the described study, which complied with all relevant regulations. Bones and teeth from 198 probable or possible goats, identified by morphological features [63, 64] were selected from 11 archaeological sites with Bronze and/or Iron Age occupations in Kazakhstan, Kyrgyzstan, Russia, and Mongolia (Fig 1; Table 1, S1 Table). Wear stages of mandibular premolars and molars were scored following Payne [65], in order to differentiate individuals in case of recovering identical haplotypes from multiple specimens at each site. Laboratory work was conducted in the Ancient DNA Laboratory of Kiel University following established protocols [6669]. Negative controls were employed during each set of DNA extractions and PCR preparations.

Table 1
Archaeological descriptions of sites are available in S1 Text.
Analytical regions and sites used in this study.
Region/SiteChronologyCountryElevation (m)DescriptionReference
Northern Eurasian Steppe
BozshakolLate Bronze AgeKazakhstan225Metallurgical complex occupied from 1600 to 1400 cal BC.[76]
Air-TauLate Bronze AgeKazakhstan260Settlement site occupied during the first half of the second millennium BC, based on ceramic typologies.Unpublished
Rublevo-VILate Bronze AgeRussia160Settlement site occupied from 1400 to 1000 cal BC.[7780]
Zamiin-UtugIron AgeMongolia1800Xiongnu cemetery site used from AD 200 to 100.Unpublished
Central Kazakhstan
TaldysaiMiddle-Late Bronze AgeKazakhstan470Metallurgical complex occupied from 1900 to 1200 cal BC.[81]
MyrzhikLate Bronze AgeKazakhstan500Settlement site occupied during the second half of the second millennium BC, based on ceramic typologies.[76, 8184]
AtasuLate Bronze AgeKazakhstan480Settlement site occupied during the second half of the second millennium BC, based on ceramic typologies.[81, 84]
Inner Asian Mountain Corridor
BegashEarly Bronze Age to modernKazakhstan800Settlement site occupied from 2400 cal BC to the early 20th c.[47, 48, 8590]
DaliEarly-Late Bronze AgeKazakhstan1500Settlement site occupied from 2700 to 1100 cal BC.[47, 89]
TasbasLate Bronze AgeKazakhstan1500Settlement site occupied from 2800 to 800 cal BC.[46, 47, 87, 8991]
Uch-KurbuLate Bronze AgeKyrgyzstan1750Cemetery site used from 1700 to 1300 cal BC.[92]

Samples were soaked in consumer bleach solution for 5 min in order to remove external contamination, followed by a soaking for 5 min and two rinses with ddH2 O. Samples were ground to a fine powder with a ball mill and digested in 960 μL 0.45 M EDTA (pH 8.0) and 40 μL proteinase K (25 mg/mL) overnight at 37°C. DNA was extracted with a large-volume extraction protocol using a guanidinium-silica based method [70]. Two independent DNA extractions from sample powder were performed for each specimen yielding an MT-CYB sequence of Capra hircus, and sample amplicons were each sequenced twice, in forward and reverse directions.

mtDNA sequence analysis

A 110-bp fragment in the cytochrome b gene for taxonomic identification among Ovis and Capra spp. was targeted using primers CapFC1 (5’- CTCTGTAACTCACATTTGTC-3’) and CapRB1b (5’- GTTTCATGTTTCTAGAAAGGT-3’ ) [71]. This initial genetic identification of species is crucial, especially since Central Asia is home to wild sheep and goats, including argali (Ovis ammon), urial (Ovis vignei), Siberian ibex (Capra sibirica), and markhor (Capra falconeri), which may be misidentified and confused with domesticated sheep and goat if relying on skeletal morphological traits. Primers CAP-FII (5’- GATCTTCCYCATGCATATAAGCA-3’) and CAP-RIII (5’- GATAAAGTTCATTAAATAGCTAC-3’) targeted a 193-bp fragment in the control region of Capra spp. containing diagnostic polymorphisms for determining haplogroups and haplotypes [60]. In case the 193-bp amplification failed, a shorter 130-bp fragment of the control region was targeted using primers CAP-FII and CAP-RII (5’- CGGGTTGCTGGTTTCAC-3’ ) [71]. Finally, while primer pairs that were used for domesticated goats were also used for Capra sibirica specimens, an additional pair of primers, CSib-F1 (5’- AACATGCGTATCCCGTCCAC-3’) and CSib-R1 (5’- GGCCCTGAAGAAAGAACCAG-3’), were designed for amplifying a 120-bp fragment of the mtDNA control region overlapping with the previously generated 193-bp amplicon, resulting in a final control region sequence of 221-bp. This allowed for better comparison with a limited number of modern Capra sibirica mtDNA sequences.

Amplifications were performed as singleplex PCRs in a final volume of 25 μL (S2 Table). Reactions were performed in a Mastercycler (Eppendorf) with the following conditions: initialization at 94°C for 10 min, 42 cycles of 94°C for 30 s, 50–60°C for 30 s (primer dependent, see S3 Table), 72°C for 30 s, and a final elongation step at 72°C for 10 min. The quality of PCR products was checked on a QIAxcel system, and those with single bands of expected size were subjected to direct sequencing with the BigDye Terminator v3.1 Kit (Applied Biosystems), following the instructions of the manufacturer. Sequencing products were purified applying the DyeEx 2.0 Spin Kit (Qiagen) and analyzed using the ABI Prism 310 and 3130 Genetic Analyzers (Applied Biosystems).

Ancient Capra hircus sequences of the control region and cytochrome b gene were concatenated for each sample, resulting in a total sequence length of 303-bp. Sequences were pair-wise aligned to one another in addition to nine previously published ancient sequences and 10 modern sequences of Capra spp. representing haplogroups A, B, C, D, F, and G (S4 Table). In cases where reference sequences of the control region did not accompany sequences of cytochrome b, we used cytochrome b sequences from sample KG170121. This approach is justified because of very low genetic variation in cytochrome b sequences for each Capra sp., which mainly served to improve phylogenetic clustering of wild Siberian ibex relative to domesticated goats. Sequences were trimmed to exclude the primer binding sites. A phylogenetic tree was constructed using MrBayes v. 3.2.6 [72] in Geneious v. 10.0.7. The simulation was set using the JC69 substitution model, invgamma rate variation, four gamma categories, and a sequence from Ovis aries as an outgroup; MCMC chain length was set to 500,000 with a sampling frequency of 200 and burn-in of 100,000. Additionally, a minimum-spanning haplotype network of sample and reference sequences was generated using PopART v. 1.7 [73]. Nucleotide diversity was calculated for each analytical region from the number of base substitutions per nucleotide position from averaging over all sequence pairs using the Kimura 2-parameter model and gamma distribution shape parameter of 1 [74], with standard error estimates by a bootstrap procedure with 1000 replicates in MEGA X v.10.1.1 [75]. Only sequences represented by 303-bp without gaps were included in the genetic distance analysis.

Results and discussion

MT-CYB sequence recovery and taxonomic identification

Out of 198 ancient skeletal specimens identified as probable goats or ambiguously identified as sheep/goat (caprine), we recovered 50 MT-CYB sequences of Capra hircus, 94 MT-CYB sequences of Ovis aries, and two MT-CYB sequences of Capra sibirica (Table 2). The relatively high overall recovery rate of 73.7% of MT-CYB sequences for Capra and Ovis spp. indicates good preservation of ancient mtDNA in the samples, which likely reflects the generally cool continental climate of sampled sites and high-elevation for sites in the IAMC, although there are appreciable differences in the recovery of MT-CYB sequences between sites (Table 2). The higher frequencies of domesticated sheep versus goat MT-CYB sequences generally echo previously published zooarchaeological data for the sampled sites and other sites across the Eurasian steppe zone, showing a ratio of sheep to goat skeletal remains of approximately 10–15:1 [46, 81, 84, 86, 93, 94].

Table 2
MT-CYB sequences from Begash, Dali, and Tasbas were previously reported [47].
Recovery rates of MT-CYB sequences by site per taxon.
Frequencies
Northern Eurasian SteppeChronologySamples testedMT-CYB Capra hircusMT-CYB Ovis ariesMT-CYB Capra sibiricaMT-CYB FailSuccess Rate
BozshakolLate Bronze Age7510185.7%
Air-TauLate Bronze Age143101028.6%
Rublevo-VILate Bronze Age52300100.0%
Zamiin-UtugIron Age22000100.0%
Central Kazakhstan
TaldysaiMiddle Bronze Age284180678.6%
Late Bronze Age251090676.0%
MyrzhikLate Bronze Age6230183.3%
AtasuLate Bronze Age200020.0%
Inner Asian Mountain Corridor
BegashEarly Bronze Age252170676.0%
Middle Bronze Age2521011252.0%
Late Bronze Age6300350.0%
Iron Age10340370.0%
Mongol Period10100100.0%
Historical Period100010.0%
DaliEarly Bronze Age274220196.3%
Middle Bronze Age53110100.0%
Late Bronze Age10100100.0%
TasbasLate Bronze Age64200100.0%
Uch-KurbuLate Bronze Age21100100.0%
TOTAL198509425273.7%

Two specimens from Begash previously identified as Siberian ibex based on morphological criteria [86] were determined to be domesticated sheep, based on MT-CYB sequences that are specific to Ovis aries [47]. Two additional Siberian ibex specimens did not yield amplicons, likely due to DNA damage resulting from burning that was visible on bone surfaces. One Siberian ibex specimen from Dali dating to the Middle Bronze Age (ca. 1800–1500 cal BC), which was identified based on more reliable cranial morphology, yielded a Capra sibirica MT-CYB sequence, while one specimen previously identified as a domesticated caprine from Begash [86] yielded a Capra sibirica MT-CYB sequence [47].

Misidentifications of Capra sibirica due to shared morphological characteristics with sheep and use of metrical determinations establishing large body-size, a trait typically associated with wild goats, highlights the need to exercise caution when attempting to separate wild and domesticated caprines. Under-identification of Siberian ibex in faunal assemblages limits our understanding of the importance of hunting by ancient pastoralists in Central Asia, which may have been a less common practice than previously suggested [86], and also obscures the significance of large-bodied, male rams in pastoralist livestock herds. In general, morphological criteria used to identify Central Asian sheep and goat draw from those developed for European and Near Eastern populations. Although the skeletal morphology of Siberian ibex is poorly studied, morphological similarity between domesticated sheep and wild ibex in some skeletal elements have been previously identified for both Near Eastern and European populations [95]. While collagen peptide mass fingerprinting (Zooarchaeology by Mass Spectrometry, aka “ZooMS”) may be used to differentiate sheep and goat specimens as a relatively cheap biomolecular method [49, 96, 97], this technique is limited to genus-level resolution and, thus, cannot differentiate between wild and domesticated lineages.

Capra hircus haplogroup diversity

Out of 50 samples that yielded Capra MT-CYB sequences, we recovered 43 control region sequences. One specimen failed to generate the 193-bp amplicon but yielded a 130-bp sequence of HVR1 using primers CAP-FII and CAP-RII. Although the primers used to amplify sequences of the control region are generally Capra-specific, one sample from Taldysai produced a sequence for Ovis aries. Control region sequences were confirmed with duplicates taken from independent DNA extractions of each specimen.

Taken together, the ancient sequences of Capra hircus represent divergent mtDNA lineages, corresponding to haplogroups A, C, and D (Table 3), which are well characterized based on modern datasets consisting of more than 3000 domesticated and wild goats [15, 17, 18]. In our sample of Capra hircus mtDNA sequences, the Late Bronze Age (ca. 1700–1000 BC) is best represented (n = 26), followed by the early-middle Bronze Age (ca. 2700–1700 BC; n = 11) and the Iron Age (ca. 400 BC-AD 200; n = 5).

Table 3
One HVR1 sequence of Siberian ibex was also recovered and is not counted here.
Recovery of Capra hircus mtDNA MT-CYB and HVR1 sequences belonging to major haplogroups by site.
Frequencies
Northern Eurasian SteppeChronologyMT-CYB Capra hircusHaplogroup AHaplogroup CHaplogroup DHVR1 FailSuccess Rate
BozshakolLate Bronze Age55000100.0%
Air-TauLate Bronze Age33000100.0%
Rublevo-VILate Bronze Age22000100.0%
Zamiin-UtugIron Age22000100.0%
Central Kazakhstan
TaldysaiMiddle Bronze Age43100100.0%
Late Bronze Age10611280.0%
MyrzhikLate Bronze Age2100150.0%
AtasuLate Bronze Age00000-
Inner Asian Mountain Corridor
BegashEarly Bronze Age2010150.0%
Middle Bronze Age200020.0%
Late Bronze Age3011166.7%
Iron Age31020100.0%
Mongol Period00000-
Historical Period00000-
DaliEarly Bronze Age4300175.0%
Middle Bronze Age33000100.0%
Late Bronze Age00000-
TasbasLate Bronze Age43010100.0%
Uch-KurbuLate Bronze Age10100100.0%
TOTAL503255884.0%

Diversity of mtDNA haplotypes was lowest for sites located in the NES, indicated by all recovered HVR1 sequences belonging to haplogroup A (Table 3) and an average genetic distance of 0.015±0.005. Low goat mtDNA diversity is in the NES is driven by the high occurrence of one HVR1 haplotype that was shared among three individual goats from Bozshakol and two from Air-Tau. This haplotype is the most common one in the dataset and was also present in five goats from sites located in central Kazakhstan (Fig 2). We are confident that the identical haplotypes from each site represent different individuals based on dental specimens being of the same body side or, in cases of specimens representing opposing body sides, individuality is confirmed by age distinction according to tooth wear stages (S1 Table). Interestingly, the remaining sample set from Taldysai is characterized by diverse HVR1 haplotypes, belonging to haplogroups A (n = 4), C (n = 2), and D (n = 1) (Table 3; Fig 2), which together exhibit an average genetic distance of 0.042±0.008.

A minimum spanning haplotype network showing the genetic diversity of domesticated and wild goats analyzed in this study and that represented by reference sequences across the major Capra hircus mtDNA haplogroups.
Fig 2
A minimum spanning haplotype network showing the genetic diversity of domesticated and wild goats analyzed in this study and that represented by reference sequences across the major Capra hircus mtDNA haplogroups.

Mitochondrial diversity was substantially higher for sites located in the IAMC than the NES, shown by sequences belonging to haplogroups A (n = 10), C (n = 3), and D (n = 4) (Table 3; Fig 2), which were unique, except for two C haplotypes from Begash and Uch-Kurbu (Fig 2). Notably, the IAMC is characterized by a high average genetic distance of 0.053±0.010, compared to the NES and central Kazakhstan.

Sequences of domesticated goats recovered from Dali, representing the earliest pastoralist occupations in the study [47], belong to haplogroup A (n = 4), while middle Bronze Age to Iron Age layers of IAMC sites contained an unexpectedly high 3:2 ratio of C and D haplotypes to A haplotypes (Table 3). This haplogroup distribution could indicate an influx of goats with diverse maternal origins into the IAMC after the middle Bronze Age, which strongly contrasts against the predominance of A lineages in post-Neolithic goats across the Near East. Although the sample size for the Iron Age is small, two out of three sequences are of lineage D, suggesting mtDNA diversity may have increased up to the end of the first century BC. This corresponds with the identification of one D lineage goat from Iron Age site of Bancheng in the Inner Mongolia Province of China [61] (Fig 1). Moreover, a recent study documenting mtDNA diversity in modern goats in Kazakhstan found that D lineages only occurred in low numbers in the southeastern region of the country compared to other regions that were dominated by A lineages [98]. Taken together, the diversity in mtDNA lineages from these studies and that in our sample reinforce the IAMC as an important locus connecting ancient Central Asian herders to those in east Asia via western China.

Sequences representing haplogroups B, G, or F were not recovered in our study. We did not expect to identify F haplotypes, since goats belonging to the F lineage are the rarest in modern populations, having only been isolated in three Sicilian domesticated goats [99], in addition to wild bezoars from a wide geographic extent spanning the Caucasus and South Asia [15, 16]. Managed goats belonging to haplogroup F (n = 5) were only present in the southern Levant during the Pre-Pottery Neolithic B, while ancient wild bezoars from southern Turkey and Armenia (n = 2) were also identified as belonging to lineage F [14]. We expected to recover B haplotypes from IAMC sites, since they are the second most abundant in modern distributions (ca. 6.5%) [15, 17], and one B haplotype was found among Bronze Age goats in southern Uzbekistan [14] and two B haplotypes were found among Iron Age goats from northern China [61]. Goats of lineage B must have passed through Central Asia as pastoralists extended the social networks through which animals and crops were initially spread across Asia. The absence of G haplotypes in our sample was difficult to predict, since goats belonging to this lineage were recovered from Neolithic sites in Iran and Turkmenistan [14], which might have been transmitted into Central Asia at similar rates as goats belonging to lineage D. Today, goats of the G lineage are largely restricted to northern and central Africa and have a very low global frequency of 0.9% [15]. The absence of B and G haplotypes in our sample demands further investigation of ancient goat genetic diversity in Iran, the IAMC, and China in order to resolve processes leading to the persistence and extinction of goat lineages in particular times and places.

Network and phylogenetic analysis

Ancient goat sequences recovered in this study are scattered throughout the A clade, which is likely an outcome of A haplotypes being the most dominant in domesticated goats after large population expansions of A groups since domestication [14, 16, 17, 19]. Haplotypes of the C lineage from Begash, Uch-Kurbu, and Taldysai cluster together and are distinct from modern reference haplotypes (Fig 2). Unexpectedly, two out of the four D haplotypes from goats in the IAMC during the Bronze and Iron Ages are newly discovered D haplotypes (Fig 2; Beg_LBA_3 and Beg_IA_1), which cluster phylogenetically as a distinct sub-clade within this haplogroup (Fig 3). This finding suggests the possibility of genetic continuity of rare domesticated goat haplotypes at Begash over a long period of time from 1600 to 20 cal BC, which spans prominent cultural turnovers during the Bronze to Iron Ages. Interestingly, no modern descendant or closely related sequences of this D lineage are known to exist. Haplotypes of lineage D are rare in contemporary world-wide distributions (0.6%) [15, 17], but are more common in bezoar populations on the Iranian plateau (ca. 13%) [15]. Haplotypes of the D lineage have also been recovered from Neolithic and Bronze Age sites in Israel, Iran, Turkmenistan, and Uzbekistan [14], which cluster with the other two Bronze and Iron Age haplotypes from the IAMC on the main D clade (Fig 2). Given our recovery of a relatively high frequency of C and D haplotypes in the IAMC, ancient frequencies of these lineages were likely far higher in ancient Central Asia than in modern distributions. The presence of these lineages further indicates high mtDNA diversity of goats in the mountain regions of Central Asia compared to the NES that is suggestive of different intensities and geographic extents of social interaction between and among these regions.

Phylogenetic tree using MrBayes 3 of ancient mtDNA sequences (HVR1 and MT-CYB regions) identified in this study (colorful symbols reflect site groups show in Fig 1), in addition to previously published ancient and contemporary reference sequences from Eurasia; Ovis aries was used as the outgroup.
Fig 3
Posterior probabilities ≥ 60% are shown on each branch. Haplogroup designations of Capra hircus are listed to the right of the sample ID numbers.Phylogenetic tree using MrBayes 3 of ancient mtDNA sequences (HVR1 and MT-CYB regions) identified in this study (colorful symbols reflect site groups show in Fig 1), in addition to previously published ancient and contemporary reference sequences from Eurasia; Ovis aries was used as the outgroup.

Due to random genetic drift ultimately leading to monomorphism, a single lineage may become dominant in a diverse population after many generations in a population that has little inward gene flow of assorted lineages [100, 101]. While it is impossible to estimate the precise population size of any particular livestock species based on zooarchaeological data, goat skeletal remains are about 10–15 times less frequent than those of sheep and about five times as infrequent as cattle remains from sites in the IAMC during the Bronze and Iron Ages [46, 86, 93], indicating that goat population sizes were relatively small compared to other domesticated species. Without a large population to buffer the effects of genetic drift, the persistence of C and D lineages in goats is likely due to high-rates of gene flow along the IAMC, likely connecting goat herds in the IAMC with diverse goat populations in southern Central Asia or further southwest on the Iranian Plateau at consistent rates. Frequent exchange events over relatively short distances could have linked goat herds of pastoralist communities as a large, global effective population, thus maintaining the prevalence of less common mtDNA lineages compared to the NES. A high degree of community interaction, likely driven by both seasonal mobility of pastoralists [20, 45, 102] and sprawling economic networks for the mining and exchange of copper and tin to produce bronzes during the Bronze Age and Iron Age [36, 103105], could have been a main factor for promoting inter-regional interactions that went hand in hand with gene flows among livestock herds at inter-regional scales.

In contrast, the NES was characterized by goats exhibiting A haplotypes with low genetic diversity, as indicated by identical sequences among NES sites and a higher connectivity within the haplotype network (Fig 2). This pattern could have been the result of domesticated goats spreading eastward across the NES from a monomorphic founding population of goats, which later during the Bronze and Iron Ages was not characterized by inward gene flows of goats exhibiting C or D haplotypes that were present in the IAMC. While the NES has long been conceived as an important ecozone of pastoralist mobility and interaction, our observation of low goat mtDNA variation in this region suggests that, at least in part, there were limited communication and exchange with southern regions in Central Asia, especially the IAMC, where rare goat haplotypes and higher mtDNA diversity persisted in the Bronze and Iron Ages.

Interestingly, Taldysai located in central Kazakhstan was characterized by an intermediate level of goat mtDNA diversity relative to the NES and IAMC (Fig 2). This finding provides an important reconfiguration in how traditional cultural interaction spheres of the Eurasian steppe zone are defined, suggesting that central Kazakhstan served as an arena of interaction between both north and south regions but did not function as a corridor for goat gene flow connecting the north and south together. This finding corresponds to recent research on regional variance in domesticated sheep astragali morphology using geometric morphometrics, showing that localized populations likely were genetically isolated from one another [106]. However, the presence of C and D goat haplotypes at Taldysai suggests repeated exchange events with southern Central Asia and the IAMC, but these lineages likely did not pass along to NES sites. Human paleogenetic studies have described a substantial amount of human gene flow from pastoralist communities in the western NES to southern Central Asia during the second millennium BC [107, 108], which would have most certainly been accompanied by goats predominantly characterized by A haplotypes. Similarly, the recent discovery by Narasimhan et al. [31] of limited human gene flows from Turan and the IAMC into central Kazakhstan during this period would explain the presence of some goat lineages at Taldysai that were recovered in high frequencies at IAMC sites.

The one ancient Capra sibirica sequence recovered from Dali phylogenetically clustered as a sister group to previously published modern Capra sibirica sequences, which is supported by a 99% Mr. Bayes posterior probability (Fig 3). However, the recovered MT-CYB and HVR1 sequences from this Capra sibirica specimen did not cluster with the Capra sibirica reference sequences in the haplotype network (Fig 2), which is likely due to highly divergent mitochondrial variation between ancient and modern populations; the former probably had far greater genetic variation than the latter. To date, our Capra sibirica sequences represent the only ancient genetic sequences recovered from Capra sibirica and demonstrates a newly discovered degree of diversity in this species, while so far being poorly characterized genetically and biogeographically [109113].

Conclusions

Despite small sample sizes, the results of this study suggest higher levels of gene flow among domesticated goats along the IAMC, which is likely a reflection of increased community interaction in the foothill zones, compared to the open steppe landscape. The northern “Northern Eurasian Steppe” is typically viewed as a metaphorical highway of cultural interaction between eastern Europe and Siberia, since the initial spread of pastoralism to the Altai Mountains via purported migrations of Afanasievo communities [26] and through the formation of expansive cultural horizons of the Iron Age [3944]. Based on the low mtDNA diversity present in Bronze Age and later Iron Age domesticated goats in the NES, there may be an unexpectedly low degree of admixture between northern steppe pastoralist communities and those in adjacent regions. The goats analyzed from sites in the NES appear to be consistent with a founder population that was likely present in eastern Europe after an earlier arrival of domesticated goats characterized by A haplotypes from the Near East via Eastern Europe or the Caucasus.

In the IAMC, we observed a possible shift in the frequencies of mtDNA lineages from the early-middle Bronze Age (2700–1600 cal BC) to the late Bronze Age and Iron Age (1600–20 cal BC), which may be a result of reconfigured socio-economic networks that facilitated transfers of goats belonging to C and D lineages to the region from diverse genetic pools from agro-pastoralist communities located in the Iranian plateau. In the absence of individual animals moving long distances along the IAMC, genetic isolation between discrete goat populations may have been greatly diminished due to directed breeding patterns of goats by pastoralists that effectively created a macro-population of goats, thus maintaining diverse maternal lineages. These findings provide additional support that the IAMC was a dynamic arena of cultural activity spanning distant regions of Central and Inner Asia, which ultimately linked these regions to contemporaneous cultural arenas in east Asia [20]. From Taldysai, the recovery of C and D haplotypes and several A haplotypes that were identical to those from NES sites suggests the region of central Kazakhstan was discretely connected to interaction spheres of communities in the NES and IAMC but may have marked a boundary between them.

More research is needed to resolve the precise nature of livestock translocations, inter-regional genetic networks, and selective pressures on livestock. This study offers a much-needed expansion of datasets for characterizing ancient goat mitochondrial diversity in Eurasia, and larger sample sizes spanning more sites along proposed transmission routes will improve the quality of inferences concerning inter-regional human interactions that can be drawn from maternal lineages. The addition of other domesticated taxa, such as sheep and cattle, to future paleogenetic analysis will help establish differences in genetic sourcing from various geographical regions. By also examining genetic variability in wild taxa, ancient biogeographic patterns may be identified to reconstruct a more complete picture of human engagements with the landscape that may reveal critical patterns in the relationship between hunting and herding strategies. The alignment of new evidence for these subsistence pursuits with recent findings of early and intensive cultivation of domesticated crops by Central Asian communities during the Bronze Age will further clarify the role that economic diversification played in shaping regional and trans-regional social interaction spheres.

Acknowledgements

We thank the Institute of Clinical Molecular Biology in Kiel for providing Sanger sequencing as supported by the DFG Clusters of Excellence “Precision Medicine in Chronic Inflammation” and “ROOTS.” We are grateful to Tanja Naujoks for technical support. We also thank Lisa Böhme and Yekaterina Dubyagina for valuable assistance during the course of this research. We acknowledge financial support by DFG within the funding programme Open Access Publizieren.

References

1 

    RW Redding. Theoretical Determinants of a Herder’s Decisions: Modeling Variation in the Sheep/Goat Ratio. In: J Clutton-Brock, C Grigson, editors. Animals and Archaeology: Early Herders and their Flocks.Oxford: BAR International Series; 1984 pp. , pp.223–241.

2 

3 

4 

    M Kaliber, N Koluman, N Silanikove. . Physiological and behavioral basis for the successful adaptation of goats to severe water restriction under hot environmental conditions.Animal. 2016;10: , pp.82–88. , doi: 10.1017/S1751731115001652

5 

    H Hongo, RH Meadow, B Öksüz, G Ilgezdi. . Animal exploitation at Çayönü Tepsei, southeastern Anatolia.Türkiye Bilim Akad Arkeol Derg. 2004; , pp.108–119. , doi: 10.22520/tubaar.2004.0009

6 

    C Makarewicz, N Tuross. . Finding Fodder and Tracking Transhumance: Isotopic Detection of Goat Domestication Processes in the Near East.Curr Anthropol.2012;53: , pp.495–505. , doi: 10.1086/665829

7 

    CA Makarewicz. . Sequential δ13C and δ18O analyses of early Holocene bovid tooth enamel: Resolving vertical transhumance in Neolithic domesticated sheep and goats.Palaeogeogr Palaeoclimatol Palaeoecol.2017;485: , pp.16–29. , doi: 10.1016/j.palaeo.2017.01.028

8 

    MC Stiner, H Buitenhuis, G Duru, SL Kuhn, SM Mentzer, ND Munro, et al. A forager–herder trade-off, from broad-spectrum hunting to sheep management at Aşıklı Höyük, Turkey. Proc Natl Acad Sci. 2014;111: , pp.8404–8409. , doi: 10.1073/pnas.1322723111

9 

    J-D Vigne, C Peters, D Helmer. The First Steps of Animal Domestication: New Archaeological Approaches. Oxford: Oxbow Books; 2005.

10 

    MA Zeder. . The Initial Domestication of Goats (Capra hircus) in the Zagros Mountains 10,000 Years Ago.Science. 2000;287: , pp.2254–2257. , doi: 10.1126/science.287.5461.2254

11 

    BS Arbuckle, L Atici. . Initial diversity in sheep and goat management in Neolithic south-western Asia.Levant. 2013;45: , pp.219–235. , doi: 10.1179/0075891413Z.00000000026

12 

13 

    D Zohary, E Tchernov, LK Horwitz. . The role of unconscious selection in the domestication of sheep and goats. J Zool. 1998;245: , pp.129–135. , doi: 10.1111/j.1469-7998.1998.tb00082.x

14 

    KG Daly, PM Delser, VE Mullin, A Scheu, V Mattiangeli, MD Teasdale, et al. Ancient goat genomes reveal mosaic domestication in the Fertile Crescent. Science. 2018;361: , pp.85–88. , doi: 10.1126/science.aas9411

15 

    L Colli, H Lancioni, I Cardinali, A Olivieri, MR Capodiferro, M Pellecchia, et al. Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability. BMC Genomics. 2015;16: , pp.1115, doi: 10.1186/s12864-015-2342-2

16 

    S Naderi, H-R Rezaei, F Pompanon, MGB Blum, R Negrini, H-R Naghash, et al. The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals. Proc Natl Acad Sci. 2008;105: , pp.17659–17664. , doi: 10.1073/pnas.0804782105

17 

    S Naderi, H-R Rezaei, P Taberlet, S Zundel, S-A Rafat, H-R Naghash, et al. Large-Scale Mitochondrial DNA Analysis of the Domestic Goat Reveals Six Haplogroups with High Diversity.PLoS ONE. 2007;2, doi: 10.1371/journal.pone.0001012

18 

    G Luikart, L Gielly, L Excoffier, J-D Vigne, J Bouvet, P Taberlet. . Multiple maternal origins and weak phylogeographic structure in domestic goats. Proc Natl Acad Sci. 2001;98: , pp.5927–5932. , doi: 10.1073/pnas.091591198

19 

    L Colli, M Milanesi, A Talenti, F Bertolini, M Chen, A Crisà, et al. Genome-wide SNP profiling of worldwide goat populations reveals strong partitioning of diversity and highlights post-domestication migration routes. Genet Sel Evol. 2018;50: , pp.58, doi: 10.1186/s12711-018-0422-x

20 

    MD Frachetti. . Multiregional Emergence of Mobile Pastoralism and Nonuniform Institutional Complexity across Eurasia.Curr Anthropol.2012;53: , pp.2–38. , doi: 10.1086/663692

21 

    GO Rollefson. . The Greening of the Badlands: Pastoral Nomads and the “Conclusion” of Neolithization in the Southern Levant.Paléorient. 2011;37: , pp.101–109. , doi: 10.3406/paleo.2011.5441

22 

    S. Fujii. Chronology of the Jafr Prehistory and Protohistory: a Key to the Process of Pastoral Nomadization in the Southern Levant. Syr Archéologie Art Hist.2013; , pp.49–125. , doi: 10.4000/syria.1723

23 

    L. Martin. Mammal remains from the eastern Jordanian Neolithic, and the nature of caprine herding in the steppe. Paléorient. 1999;25: , pp.87–104. , doi: 10.3406/paleo.1999.4689

24 

    H Miller, D Baird, J Pearson, AL Lamb, M Grove, L Martin, et al. The origins of nomadic pastoralism in the eastern Jordanian steppe: a combined stable isotope and chipped stone assessment.Levant. 2019;0: , pp.1–24. , doi: 10.1080/00758914.2019.1651560

25 

    EN Chernykh. . Formation of the Eurasian “steppe belt” of stockbreeding cultures: viewed through the prism of archaeometallurgy and radiocarbon dating.Archaeol Ethnol Anthropol Eurasia. 2008;35: , pp.36–53. , doi: 10.1016/j.aeae.2008.11.003

26 

    D. AnthonyThe horse, the wheel, and language: how bronze-age riders from the Eurasian steppes shaped the modern world. Princeton: Princeton University Press; 2007.

27 

    VA Posrednikov. . Pit-grave culture migration to the east and the Afanasievo-proto-Tocharian problem.Donezky Arkheologiya Sbornik.1992 pp. , pp.9–20.

28 

    Gryaznov MP, Vadetskaya EB. Afanasievskaya kul’tura. Istoriya Sibiri. Leningrad; 1968. pp. 159–165.

29 

    AP Pogozheva. Afanas’evskaya Kul’tura. In: AP Pogozheva, MP Rykun, NF Stepanova, SS Tup, editors. Epokha Eneolita i Bronzy Gornogo Altaya: Chast’ 1.Barnaul: Azbuka; 2006 pp. , pp.18–48.

30 

    ME Allentoft, M Sikora, K-G Sjögren, S Rasmussen, M Rasmussen, J Stenderup, et al. Population genomics of Bronze Age Eurasia. Nature. 2015;522: , pp.167–172. , doi: 10.1038/nature14507

31 

    VM Narasimhan, N Patterson, P Moorjani, N Rohland, R Bernardos, S Mallick, et al. The formation of human populations in South and Central Asia. Science. 2019;365: , pp.eaat7487, doi: 10.1126/science.aat7487

32 

    W Haak, I Lazaridis, N Patterson, N Rohland, S Mallick, B Llamas, et al. Massive migration from the steppe was a source for Indo-European languages in Europe. Nature. 2015;522: , pp.207, doi: 10.1038/nature14317

33 

    I Mathieson, I Lazaridis, N Rohland, S Mallick, N Patterson, SA Roodenberg, et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature. 2015;528: , pp.499–503. , doi: 10.1038/nature16152

34 

    D. AnthonyThe opening of the Eurasian steppe at 2000 In: V Mair, editor. The Bronze Age and early Iron Age peoples of Eastern Central Asia. Washington (D.C.): Institute for the Study of Man; 1998 pp. , pp.94–113.

35 

    EN Chernykh. Formation of the Eurasian Steppe Belt Cultures Viewed through the Lens of Archaeometallurgy and Radiocarbon Dating. In: BK Hanks, K Linduff, editors. Social complexity in prehistoric Eurasia: monuments, metals, and mobility.Cambridge: Cambridge University Press; 2009 pp. , pp.115–145.

36 

    EN Chernykh. Ancient metallurgy in the USSR: the early metal age.New York: Cambridge University Press; 1992.

37 

    C Gaunitz, A Fages, K Hanghøj, A Albrechtsen, N Khan, M Schubert, et al. Ancient genomes revisit the ancestry of domestic and Przewalski’s horses. Science. 2018;360: , pp.111–114. , doi: 10.1126/science.aao3297

38 

    A Fages, K Hanghøj, N Khan, C Gaunitz, A Seguin-Orlando, M Leonardi, et al. Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series. Cell. 2019;177: , pp.1419–1435.e31. , doi: 10.1016/j.cell.2019.03.049

39 

    M Unterländer, F Palstra, I Lazaridis, A Pilipenko, Z Hofmanová, M Groß, et al. Ancestry and demography and descendants of Iron Age nomads of the Eurasian Steppe.Nat Commun.2017;8: , pp.14615, doi: 10.1038/ncomms14615

40 

    A Juras, M Krzewińska, AG Nikitin, E Ehler, M Chyleński, S Łukasik, et al. Diverse origin of mitochondrial lineages in Iron Age Black Sea Scythians.Sci Rep.2017;7: , pp.43950, doi: 10.1038/srep43950

41 

    C. BeckwithEmpires of the Silk Road a History of Central Eurasia from the Bronze Age to the Present. Princeton: Princeton University Press; 2009.

42 

    JD Rogers. . Inner Asian States and Empires: Theories and Synthesis.J Archaeol Res. 2012;20: , pp.205–256. , doi: 10.1007/s10814-011-9053-2

43 

    K. MacFarland. Cycles of historical inquiry: evaluating the historical paradigm of nomadism in Iron Age north central Eurasia. World Archaeol. 2019;0: , pp.1–16. , doi: 10.1080/00438243.2019.1604259

44 

    W. Honeychurch. Alternative Complexities: The Archaeology of Pastoral Nomadic States. J Archaeol Res.2014;22: , pp.277–326. , doi: 10.1007/s10814-014-9073-9

45 

    MD Frachetti, CE Smith, CM Traub, T Williams. . Nomadic ecology shaped the highland geography of Asia’s Silk Roads. Nature. 2017;543: , pp.193–198. , doi: 10.1038/nature21696

46 

    PN Doumani, MD Frachetti, R Beardmore, TM Schmaus, RN Spengler, AN Mar’yashev. . Burial ritual, agriculture, and craft production among Bronze Age pastoralists at Tasbas (Kazakhstan).Archaeol Res Asia.2015;, pp.1–2: 17–32. , doi: 10.1016/j.ara.2015.01.001

47 

    TR Hermes, MD Frachetti, PN Doumani Dupuy, A Mar’yashev, A Nebel, CA Makarewicz. . Early integration of pastoralism and millet cultivation in Bronze Age Eurasia. Proc R Soc B. 2019;286, doi: 10.1098/rspb.2019.1273

48 

    MD Frachetti, RN Spengler, GJ Fritz, AN Mar’yashev. . Earliest direct evidence for broomcorn millet and wheat in the central Eurasian steppe region.Antiquity. 2010;84: , pp.993–1010. , doi: 10.1017/S0003598X0006703X

49 

    W Taylor, S Shnaider, A Abdykanova, A Fages, F Welker, F Irmer, et al. Early pastoral economies along the Ancient Silk Road: Biomolecular evidence from the Alay Valley, Kyrgyzstan.PLOS ONE. 2018;13: , pp.e0205646, doi: 10.1371/journal.pone.0205646

50 

    RN Spengler, F Maksudov, E Bullion, A Merkle, T Hermes, M Frachetti. . Arboreal crops on the medieval Silk Road: Archaeobotanical studies at Tashbulak.PLOS ONE. 2018;13: , pp.e0201409, doi: 10.1371/journal.pone.0201409

51 

    S. BarisitzCentral Asia and the Silk Road: Economic Rise and Decline over Several Millennia. Cham: Springer; 2017.

52 

    TR Hermes, MD Frachetti, EA Bullion, F Maksudov, S Mustafokulov, CA Makarewicz. . Urban and nomadic isotopic niches reveal dietary connectivities along Central Asia’s Silk Roads.Sci Rep.2018;8: , pp.5177, doi: 10.1038/s41598-018-22995-2

53 

    GM Matuzeviciute, K Tabaldiev, T Hermes, E Ananyevskaya, M Grikpedis, E Luneau, et al. High-Altitude Agro-Pastoralism in the Kyrgyz Tien Shan: New Excavations of the Chap Farmstead (1065–825 cal b.c.).J Field Archaeol.2020;45: , pp.29–45. , doi: 10.1080/00934690.2019.1672128

54 

    F Maksudov, E Bullion, ER Henry, T Hermes, AM Merkle, MD Frachetti. Nomadic Urbanism at Tashbulak: A New Highland Town of the Qarakhanids. In: C Baumer, M Novák, editors. Urban Cultures of Central Asia from the Bronze Age to the Karakhanids: Learnings and conclusions from new archaeological investigations and discoveries. Wiesbaden: Harrassowitz; 2019 pp. , pp.283–306.

55 

    P Gerbault, A Powell, MG Thomas. . Evaluating demographic models for goat domestication using mtDNA sequences.Anthropozoologica. 2012;47: , pp.64–76. , doi: 10.5252/az2012n2a5

56 

    I Pfeiffer, J Burger, B Brenig. . Diagnostic polymorphisms in the mitochondrial cytochrome b gene allow discrimination between cattle, sheep, goat, roe buck and deer by PCR-RFLP. BMC Genet. 2004;5: , pp.30, doi: 10.1186/1471-2156-5-30

57 

    H Fernández, S Hughes, J-D Vigne, D Helmer, G Hodgins, C Miquel, et al. Divergent mtDNA lineages of goats in an Early Neolithic site, far from the initial domestication areas. Proc Natl Acad Sci. 2006;103: , pp.15375–15379. , doi: 10.1073/pnas.0602753103

58 

    S Kadowaki, K Ohnishi, S Arai, F Guliyev, Y Nishiaki. . Mitochondrial DNA Analysis of Ancient Domestic Goats in the Southern Caucasus: A Preliminary Result from Neolithic Settlements at Göytepe and Hacı Elamxanlı Tepe.Int J Osteoarchaeol. 2017;27: , pp.245–260. , doi: 10.1002/oa.2534

59 

    FH Mazdarani, MT Akbari, RMN Fard, M Hessari, KC Pour, . others. Molecular identification of Capra hircus in East Chia Sabz, an Iranian pre-pottery Neolithic site, Central Zagros, based on mtDNA. J Anim Plant Sci. 2014;24: , pp.945–950.

60 

    I Akis, V Onar, N Toker, O Belli, G Pazvant, K Oztabak. . Ancient DNA Analysis of Anatolian Goat Remains Excavated from a Urartian Castle in Eastern Turkey.Int J Osteoarchaeol. 2016;26: , pp.246–254. , doi: 10.1002/oa.2415

61 

    L Han, H-X Yu, D-W Cai, H-L Shi, H Zhu, H Zhou. . Mitochondrial DNA analysis provides new insights into the origin of the Chinese domestic goat. Small Rumin Res. 2010;90: , pp.41–46. , doi: 10.1016/j.smallrumres.2009.12.011

62 

    O Schröder, M Wagner, S Wutke, Y Zhang, Y Ma, D Xu, et al. Ancient DNA identification of domestic animals used for leather objects in Central Asia during the Bronze Age. The Holocene. 2016;26: , pp.1722–1729. , doi: 10.1177/0959683616641741

63 

    S. Payne. Morphological distinctions between the mandibular teeth of young sheep, Ovis, and goats, Capra. J Archaeol Sci. 1985;12: , pp.139–147. , doi: 10.1016/0305-4403(85)90058-5

64 

    MA Zeder, HA Lapham. . Assessing the reliability of criteria used to identify postcranial bones in sheep, Ovis, and goats, Capra.J Archaeol Sci. 2010;37: , pp.2887–2905. , doi: 10.1016/j.jas.2010.06.032

65 

    S. Payne. Kill-off Patterns in Sheep and Goats: The Mandibles from Aşvan Kale. Anatol Stud.1973;23: , pp.281–303. , doi: 10.2307/3642547

66 

67 

68 

    S Champlot, C Berthelot, M Pruvost, EA Bennett, T Grange, E-M Geigl. . An Efficient Multistrategy DNA Decontamination Procedure of PCR Reagents for Hypersensitive PCR Applications.PLOS ONE. 2010;5: , pp.e13042, doi: 10.1371/journal.pone.0013042

69 

    B Krause-Kyora, M Nutsua, L Boehme, F Pierini, DD Pedersen, S-C Kornell, et al. Ancient DNA study reveals HLA susceptibility locus for leprosy in medieval Europeans. Nat Commun. 2018;9: , pp.1569, doi: 10.1038/s41467-018-03857-x

70 

    N Rohland, M Hofreiter. . Ancient DNA extraction from bones and teeth. Nat Protoc. 2007;2: , pp.1756–1762. , doi: 10.1038/nprot.2007.247

71 

    H Fernández, S Hughes, J-D Vigne, D Helmer, G Hodgins, C Miquel, et al. Divergent mtDNA lineages of goats in an Early Neolithic site, far from the initial domestication areas. Proc Natl Acad Sci. 2006;103: , pp.15375–15379. , doi: 10.1073/pnas.0602753103

72 

    JP Huelsenbeck, F Ronquist. . MRBAYES: Bayesian inference of phylogenetic trees. Bioinforma Oxf Engl. 2001;17: , pp.754–755. , doi: 10.1093/bioinformatics/17.8.754

73 

    JW Leigh, D Bryant. . popart: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6: , pp.1110–1116. , doi: 10.1111/2041-210X.12410

74 

    M. Kimura. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16: , pp.111–120. , doi: 10.1007/BF01731581

75 

    S Kumar, G Stecher, M Li, C Knyaz, K Tamura. . MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35: , pp.1547–1549. , doi: 10.1093/molbev/msy096

76 

    J-S Park, D Voyakin, A Beisenov. . The implication of diachronic changes reflected in LBA bronze assemblages of Central Kazakhstan.Archaeol Anthropol Sci.2020;12: , pp.7, doi: 10.1007/s12520-019-00989-z

77 

    DV Papin. . Khronologiya pamyatnikov epokhi pozdnej bronzy stepnogo i lesostepnogo Altaya.Vestn Kemer Gos Univ. 2015;6: , pp.135–138.

78 

    DV Papin, VG Loman, NF Stepanova, AC Fedoruk. . Rezultaty tekhniko-tekhnologicheskogo analiza keramicheskogo kompleksa poseleniya epokhi pozdnej bronzy Rublevo-VI.Teor Prakt Arkheologiskikh Issled. 2015;2: , pp.115–143.

79 

    DV Papin, AC Fedoruk. Khozyaistvenno-kul’turnyi tsentr epokhi pozdnej bronzy Rublevo-VI.Arkheologicheskiye issledovaniya stepnoj Evrazii. Karaganda: Tengri Ltd.;2013 pp. , pp.129–149.

80 

    YuF Kiryushin, PA Kosintsev, DV Papin, AC Fedoruk. Voprosy khozyajstvennoj deyatel’nosti naseleniya stepnogo Ob’-Irtysh’ya v epokhu pozdnej bronzy. In: YuF Kiryushin, editor. Khozyajstvenno-kul’turnye traditsii Altaya v epokhu bronzy.Barnaul: Izd-vo Alt. un-ta; 2010 pp. , pp.112–127.

81 

    OA Artyukhova, Zh Kurmankulov, A Yermolayeva, A Erzhanova. Kompleks pamyatnikov v urochische Taldysai. Almaty: Margulan Institute of Archaeology; 2013.

82 

    AZ Beisenov. . Pamyarniki verkhov’ev reki Atasu v tsentral’nom Kazakhstane.Vestn Tomsk Gos Univ. 2015;34: , pp.111–122. , doi: 10.17223/19988613/35/17

83 

    MK Kadyrbaev, Zh Kurmankulov. Kul’tura drevnikh skotovodov i metallurgov Sary-Arki. Alma-Ata: Gylym; 1992.

84 

    AK Outram, A Kasparov, NA Stear, V Varfolomeev, E Usmanova, RP Evershed. . Patterns of pastoralism in later Bronze Age Kazakhstan: new evidence from faunal and lipid residue analyses.J Archaeol Sci. 2012;39: , pp.2424–2435. , doi: 10.1016/j.jas.2012.02.009

85 

    MD Frachetti, AN Mar’yashev. . Long-Term Occupation and Seasonal Settlement of Eastern Eurasian Pastoralists at Begash, Kazakhstan.J Field Archaeol.2007;32: , pp.221–242. , doi: 10.1179/009346907791071520

86 

    MD Frachetti, N Benecke. . From sheep to (some) horses: 4500 years of herd structure at the pastoralist settlement of Begash (south-eastern Kazakhstan).Antiquity. 2009;83: , pp.1023–1037. , doi: 10.1017/S0003598X00099324

87 

    TM Schmaus, PN Doumani Dupuy, MD Frachetti. . Variability in seasonal mobility patterns in Bronze and Iron Age Kazakhstan through cementum analysis.Quat Int.2019 [cited 6 May 2020]. , doi: 10.1016/j.quaint.2019.04.018

88 

    R Spengler, M Frachetti, P Doumani, L Rouse, B Cerasetti, E Bullion, et al. Early agriculture and crop transmission among Bronze Age mobile pastoralists of Central Eurasia. Proc R Soc B. 2014;281: , pp.20133382, doi: 10.1098/rspb.2013.3382

89 

    P. DoumaniCeramic Technology of Bronze Age Nomadic Pastoralists of Semirech’ye, Kazakhstan. Doctoral Dissertation, Washington University in St. Louis2014.

90 

    PN Doumani, M Frachetti. . Bronze Age textile evidence in ceramic impressions: weaving and pottery technology among mobile pastoralists of central Eurasia.Antiquity. 2012;86: , pp.368–382. , doi: 10.1017/S0003598X00062827

91 

    RN Spengler, MD Frachetti, PN Doumani. . Late Bronze Age agriculture at Tasbas in the Dzhungar Mountains of eastern Kazakhstan.Quat Int.2014;348: , pp.147–157. , doi: 10.1016/j.quaint.2014.03.039

92 

    G Motuzaite Matuzeviciute, A Abdykanova, S Kume, Y Nishiaki, K Tabaldiev. . The effect of geographical margins on cereal grain size variation: Case study for highlands of Kyrgyzstan.J Archaeol Sci Rep. 2018;20: , pp.400–410. , doi: 10.1016/j.jasrep.2018.04.037

93 

    C Chang, N Benecke, FP Grigoriev, AM Rosen, PA Tourtellotte. . Iron Age society and chronology in South-east Kazakhstan.Antiquity. 2003;77: , pp.298–312. , doi: 10.1017/S0003598X00092280

94 

    A. Haruda. Regional pastoral practice in central and southeastern Kazakhstan in the Final Bronze Age (1300–900BCE).Archaeol Res Asia. 2018;15: , pp.146–156. , doi: 10.1016/j.ara.2017.09.004

95 

    A Von den Driesch, U Wodtke. The Fauna of ‘Ain Ghazal, a Major PPN and Early PN Settlement in Central Jordan. In: HGK Gebel, GO Rollefson, editors. The Prehistory of Jordan, II.Berlin: Ex Oriente; 1997 pp. , pp.511–556.

96 

    M Buckley, S Whitcher Kansa, S Howard, S Campbell, J Thomas-Oates, M Collins. . Distinguishing between archaeological sheep and goat bones using a single collagen peptide.J Archaeol Sci. 2010;37: , pp.13–20. , doi: 10.1016/j.jas.2009.08.020

97 

    WTT Taylor, J Clark, J Bayarsaikhan, T Tuvshinjargal, JT Jobe, W Fitzhugh, et al. Early Pastoral Economies and Herding Transitions in Eastern Eurasia.Sci Rep.2020;10: , pp.1–15. , doi: 10.1038/s41598-019-56847-4

98 

    R Tabata, F Kawaguchi, S Sasazaki, Y Yamamoto, M Bakhtin, P Kazymbet, et al. The Eurasian Steppe is an important goat propagation route: A phylogeographic analysis using mitochondrial DNA and Y-chromosome sequences of Kazakhstani goats. Anim Sci J. 2019;90: , pp.317–322. , doi: 10.1111/asj.13144

99 

    MT Sardina, M Ballester, J Marmi, R Finocchiaro, JBCHMV Kaam, B Portolano, et al. Phylogenetic analysis of Sicilian goats reveals a new mtDNA lineage. Anim Genet. 2006;37: , pp.376–378. , doi: 10.1111/j.1365-2052.2006.01451.x

100 

    JC Avise, JE Neigel, J Arnold. . Demographic influences on mitochondrial DNA lineage survivorship in animal populations. J Mol Evol. 1984;20: , pp.99–105. , doi: 10.1007/BF02257369

101 

    M Nei, F Tajima. . Genetic Drift and Estimation of Effective Population Size. Genetics. 1981;98: , pp.625–640.

102 

    MD Frachetti. . Migration Concepts in Central Eurasian Archaeology.Annu Rev Anthropol.2011;40: , pp.195–212. , doi: 10.2307/41287728

103 

    K Linduff, editor. Metallurgy in ancient eastern Eurasia from the Urals to the Yellow River. Lewiston NY: Edwin Mellen Press; 2004.

104 

    K Linduff, J Mei. . Metallurgy in Ancient Eastern Asia: Retrospect and Prospects.J World Prehistory. 2009;22: , pp.265–281. , doi: 10.1007/s10963-009-9023-5

105 

    EN Chernykh, LB Orlovskaya. . Radiouglerodnaya khronologiya drevneyamnoy obshcnosti i istori kurhannykh kul’tur.Ross Arkheologiya.2004; , pp.84–99.

106 

    AF Haruda, V Varfolomeev, A Goriachev, A Yermolayeva, AK Outram. . A new zooarchaeological application for geometric morphometric methods: Distinguishing Ovis aries morphotypes to address connectivity and mobility of prehistoric Central Asian pastoralists.J Archaeol Sci. 2019;107: , pp.50–57. , doi: 10.1016/j.jas.2019.05.002

107 

    B Damgaard P de, R Martiniano, J Kamm, JV Moreno-Mayar, G Kroonen, M Peyrot, et al. The first horse herders and the impact of early Bronze Age steppe expansions into Asia. Science. 2018; , pp.eaar7711, doi: 10.1126/science.aar7711

108 

    B Damgaard P de, N Marchi, S Rasmussen, M Peyrot, G Renaud, T Korneliussen, et al. 137 ancient human genomes from across the Eurasian steppes. Nature. 2018;557: , pp.369–374. , doi: 10.1038/s41586-018-0094-2

109 

    N Pidancier, S Jordan, G Luikart, P Taberlet. . Evolutionary history of the genus Capra (Mammalia, Artiodactyla): Discordance between mitochondrial DNA and Y-chromosome phylogenies.Mol Phylogenet Evol. 2006;40: , pp.739–749. , doi: 10.1016/j.ympev.2006.04.002

110 

    B Khan, A Ablimit, G Khan, AW Jasra, H Ali, R Ali, et al. Abundance, distribution and conservation status of Siberian ibex, Marco Polo and Blue sheep in Karakoram-Pamir mountain area.J King Saud Univ—Sci.2016;28: , pp.216–225. , doi: 10.1016/j.jksus.2015.02.007

111 

    RP Reading, S Amgalanbaatar, D Kenny, A DeNicola, E Tuguldur. . Siberian ibex (Capra sibirica) home ranges in Ikh Nart Nature Reserve, Mongolia: preliminary findings.Mong J Biol Sci. 2007;5: , pp.29–37. , doi: 10.22353/mjbs.2007.05.05

112 

    A Hassanin, A Ropiquet, A Couloux, C Cruaud. . Evolution of the Mitochondrial Genome in Mammals Living at High Altitude: New Insights from a Study of the Tribe Caprini (Bovidae, Antilopinae).J Mol Evol. 2009;68: , pp.293–310. , doi: 10.1007/s00239-009-9208-7

113 

    A Ropiquet, A Hassanin. . Hybrid origin of the Pliocene ancestor of wild goats. Mol Phylogenet Evol. 2006;41: , pp.395–404. , doi: 10.1016/j.ympev.2006.05.033

10 Dec 2019

PONE-D-19-29815

High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor

PLOS ONE

Dear Dr. Hermes,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.

==============================

all comments need to be addressed.

==============================

We would appreciate receiving your revised manuscript by Jan 24 2020 11:59PM. When you are ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter.

To enhance the reproducibility of your results, we recommend that if applicable you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. For instructions see: http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols

Please include the following items when submitting your revised manuscript:

    A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). This letter should be uploaded as separate file and labeled 'Response to Reviewers'.
    A marked-up copy of your manuscript that highlights changes made to the original version. This file should be uploaded as separate file and labeled 'Revised Manuscript with Track Changes'.
    An unmarked version of your revised paper without tracked changes. This file should be uploaded as separate file and labeled 'Manuscript'.

Please note while forming your response, if your article is accepted, you may have the opportunity to make the peer review history publicly available. The record will include editor decision letters (with reviews) and your responses to reviewer comments. If eligible, we will contact you to opt in or out.

We look forward to receiving your revised manuscript.

Kind regards,

Peter F. Biehl, PhD

Academic Editor

PLOS ONE

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

1) Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

http://www.journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and http://www.journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2) In your manuscript, please provide additional information regarding the specimens used in your study. Ensure that you have reported specimen numbers and complete repository information, including museum name and geographic location. If permits were required, please ensure that you have provided details for all permits that were obtained, including the full name of the issuing authority. For more information on PLOS ONE's requirements for paleontology and archaeology research, see https://journals.plos.org/plosone/s/submission-guidelines#loc-paleontology-and-archaeology-research.

3) We note that Figure 1 in your submission contains map/satellite images which may be copyrighted. All PLOS content is published under the Creative Commons Attribution License (CC BY 4.0), which means that the manuscript, images, and Supporting Information files will be freely available online, and any third party is permitted to access, download, copy, distribute, and use these materials in any way, even commercially, with proper attribution. For these reasons, we cannot publish previously copyrighted maps or satellite images created using proprietary data, such as Google software (Google Maps, Street View, and Earth). For more information, see our copyright guidelines: http://journals.plos.org/plosone/s/licenses-and-copyright.

We require you to either (1) present written permission from the copyright holder to publish these figures specifically under the CC BY 4.0 license, or (2) remove the figures from your submission:

    • You may seek permission from the original copyright holder of Figure(s) [#] to publish the content specifically under the CC BY 4.0 license. 

We recommend that you contact the original copyright holder with the Content Permission Form (http://journals.plos.org/plosone/s/file?id=7c09/content-permission-form.pdf) and the following text:

“I request permission for the open-access journal PLOS ONE to publish XXX under the Creative Commons Attribution License (CCAL) CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). Please be aware that this license allows unrestricted use and distribution, even commercially, by third parties. Please reply and provide explicit written permission to publish XXX under a CC BY license and complete the attached form.”

Please upload the completed Content Permission Form or other proof of granted permissions as an "Other" file with your submission.

In the figure caption of the copyrighted figure, please include the following text: “Reprinted from [ref] under a CC BY license, with permission from [name of publisher], original copyright [original copyright year].”

    • If you are unable to obtain permission from the original copyright holder to publish these figures under the CC BY 4.0 license or if the copyright holder’s requirements are incompatible with the CC BY 4.0 license, please either i) remove the figure or ii) supply a replacement figure that complies with the CC BY 4.0 license. Please check copyright information on all replacement figures and update the figure caption with source information. If applicable, please specify in the figure caption text when a figure is similar but not identical to the original image and is therefore for illustrative purposes only.

The following resources for replacing copyrighted map figures may be helpful:

USGS National Map Viewer (public domain): http://viewer.nationalmap.gov/viewer/

The Gateway to Astronaut Photography of Earth (public domain): http://eol.jsc.nasa.gov/sseop/clickmap/

Maps at the CIA (public domain): https://www.cia.gov/library/publications/the-world-factbook/index.html and https://www.cia.gov/library/publications/cia-maps-publications/index.html

NASA Earth Observatory (public domain): http://earthobservatory.nasa.gov/

Landsat: http://landsat.visibleearth.nasa.gov/

USGS EROS (Earth Resources Observatory and Science (EROS) Center) (public domain): http://eros.usgs.gov/#

Natural Earth (public domain): http://www.naturalearthdata.com/

4) Thank you for stating the following in the Competing Interests section:

"The authors have declared that no competing interests exist.".

We note that one or more of the authors are employed by a commercial company:Archaelogical Expertise, LLC, Almaty, Kazakhstan.

    • Please provide an amended Funding Statement declaring this commercial affiliation, as well as a statement regarding the Role of Funders in your study. If the funding organization did not play a role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript and only provided financial support in the form of authors' salaries and/or research materials, please review your statements relating to the author contributions, and ensure you have specifically and accurately indicated the role(s) that these authors had in your study. You can update author roles in the Author Contributions section of the online submission form.

Please also include the following statement within your amended Funding Statement.

“The funder provided support in the form of salaries for authors [insert relevant initials], but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the ‘author contributions’ section.”

If your commercial affiliation did play a role in your study, please state and explain this role within your updated Funding Statement.

2. Please also provide an updated Competing Interests Statement declaring this commercial affiliation along with any other relevant declarations relating to employment, consultancy, patents, products in development, or marketed products, etc.  

Within your Competing Interests Statement, please confirm that this commercial affiliation does not alter your adherence to all PLOS ONE policies on sharing data and materials by including the following statement: "This does not alter our adherence to  PLOS ONE policies on sharing data and materials.” (as detailed online in our guide for authors http://journals.plos.org/plosone/s/competing-interests) . If this adherence statement is not accurate and  there are restrictions on sharing of data and/or materials, please state these. Please note that we cannot proceed with consideration of your article until this information has been declared.

Please include both an updated Funding Statement and Competing Interests Statement in your cover letter. We will change the online submission form on your behalf.

Please know it is PLOS ONE policy for corresponding authors to declare, on behalf of all authors, all potential competing interests for the purposes of transparency. PLOS defines a competing interest as anything that interferes with, or could reasonably be perceived as interfering with, the full and objective presentation, peer review, editorial decision-making, or publication of research or non-research articles submitted to one of the journals. Competing interests can be financial or non-financial, professional, or personal. Competing interests can arise in relationship to an organization or another person. Please follow this link to our website for more details on competing interests: http://journals.plos.org/plosone/s/competing-interests

5) Please include captions for your Supporting Information files at the end of your manuscript, and update any in-text citations to match accordingly. Please see our Supporting Information guidelines for more information: http://journals.plos.org/plosone/s/supporting-information.

Additional Editor Comments (if provided):

Your manuscript has now been seen by two referees, whose comments are appended below. You will see from these comments that the referees find your work of high interest. In light of these comments, only minor revisions are needed - see especially reviewer 2 with the suggestion to move the site description to the supplement.

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: This manuscript is extremely well written in a very clear and well-organized style. It will require very little correcting, as I failed to find any typos or grammatical errors. The authors are highly qualified, with extensive years of experience in Asian steppe archaeology, focused on the Bronze and Iron Age. The ancient mtDNA, stable isotope and cementum analysis, as well as standard morphological analysis of plant and animal remains form a solid, holistic picture of goat herding in the Djungar Mountains and surrounding regions of Siberia, northern and central Kazakhstan and northwestern Mongolia. The sites from which samples were derived were selected in a meaningful way to capture the evidence to develop their results on subsistence and caprine breeding and maintenance. The sample sizes were ample, although future expansion of data is certainly welcome. Here I would like to briefly summarize their most significant results.

The found that diversity of mtDNA haplotypes was lower for NES sites and higher for IAMC sites, which has major implications for these two regions. They interpret some of the genetic findings as an indication of connections with western China. The absence of B and G haplotypes is interesting and shows the need to increase sample sizes from Iran and China before the migrations of goat populations through Central Asia can be better understood.

One interesting and potentially very significant finding was that morphological identification of wild species of caprines conflicts with their ancient mtDNA results. This casts doubt on previous publications that report hunting of wild species when, in fact, the remains may just be from large males in the domestic livestock. This would mean that hunting was less significant than originally thought, at least for caprines, and that size variation within domestic herds is something worth considering. It could imply that selection for particular characteristics and perhaps even breed development was being initiated in these livestock populations or, as they indicate, the influx of goats from outside after the middle Bronze Age.

Another important development is the assessment of the low frequencies of goats vs. sheep and cattle. This can shed light on the importance of their uses. Goats tend to be bigger milk producers than sheep, and, though the hair of certain breeds is highly valued, sheep would yield wool. The amount of meat they produce is basically the same, so perhaps these populations depended more on cow milk than goat and considered wool production more important than goat hair (if the goats had cashmere or angora-type coats, which is unknown at present-they may have had a short coat).

This paper sheds light on the degree of human population migration owing to pastoral mobility and the bronze trade in the Bronze and Iron Ages. For the NES, there seems to have been less mobility and exchange than in the IAMC. This research marks a significant beginning in understanding connections between human and livestock populations, most notably the activities in Central Kazakhstan as a possible intermediate region, but also a border between the north and south. They interpret the low diversity in the NEC as reflecting the eastward migration of Near Eastern goats from Eastern Europe via the Afanasievo culture (and perhaps others) to the Altai region. The IAMC, on the other hand, was more dynamic, with connections potentially with Iran and China. This exciting research is definitely a great start, but, as they point out, their tantalizing results show that much more work needs to be done in the future to fine-tune our understanding of mobility of herds and people in Central Asia and connections between Eastern Europe, the Near East, and the far East. I look forward to their continuing research and new results and hope that they will also be able to address issues such as breed development and product exploitation (meat, coat, and dairy) through genetic and traditional archaeological information.

Reviewer #2: I think this paper presents novel and interesting data and is professionally prepared and presented. I suggest a few minor changes and reorganizations but otherwise supports its acceptance for publication.

In this paper, the authors present new mt DNA sequences for archaeological (Bronze and Iron Age) goats from two regions 'Central Asia' including the North Eurasian Steppe (NAS) and the Inner Asian Mountain Corridor (IAMC). The authors present mt sequences and haplogroup assignments for 50 domestic goat specimens from 10 sites. This itself is a valuable contribution to the understanding of the historical genetics of goat populations as well as the migration and interaction histories of these two regions. Even more interesting is the patterns the authors identify in the mt data. They point out that goats from the NAS region have low mt diversity, all characterized as mt haplogroup A. In contrast, the goats from the IAMC region reflect higher variation with low frequencies of haplogroup A while also exhibiting haplogroups C and D. From these patterns the authors link the IAMC region with goat diversity on the Iranian plateau whereas the A lineages in the NAS perhaps derive from eastern Europe or the Caucasus region and thus reflect divergent (goat) colonization histories.

These are very interesting data and patterns that fit into the very high profile current discussions of the role of mobility, trade and exchange, and population movements in and through these central Asian corridor regions.

My main concern with reading the paper comes from the organization. I feel like the section describing the sites is intrusive and disrupts the narrative. I suggest moving the site descriptions to a supplement rather than including them in the main text of the paper. I also think it should be easier to identify which specimens represent which haplogroup--this should be included in the main specimen table in Supplement table 1.

Specific and generally minor comments/questions/suggestions:

Line 85: what do the authors mean by saying that mobile pastoralism emerged in the 5th millennium BC? What does ‘mobile’ mean?

Line 102: descriptions of the routes for the spread of domestic animals make it sound like a simple process of migration. Does not address the fact that these regions were occupied by local hunter-gatherers and there are active discussion concerning the mechanisms by which domesticates moved into this broad region.

Line 196-330: this section describing the sites is usually moved to a supplemental file in aDNA publications. I would recommend moving this text and adding a Table to the main text with the site names, dates and perhaps a very brief description of the most relevant details.

Line 365: I am not sure that this sentence is relevant here since it is not referring specifically to

Capra sibirica which, based on Pidancier et al 2006, seems to be a more basal member of the Capra genus, although Sardina et al 2006’s analysis shows Capra sibirica closely related to bezoar and domestic goat lineages. Either way, it is not clear that this comment is relevant unless it is specifically directed at Capra sibirica morphology rather than ‘ibexes’ more generally.

Line 420: lineage D is described as ‘rare’. It is rare based on a study of modern goats in Naderi et al 2007 but based on Daly et al 2018, lineage D is not all that rare—it is present in Neolithic and Bronze Age Iran as well as Chalcolithic Israel. Based on ancient samples in Daly et al 2018 and Colli et al 2015, C is much more rare and is found in wild Iranian bezoar as well as Neolithic domestic goats in the Balkans, modern domestic goats in central Europe, and modern domestic goats in Iran. As a result, C is a less good indicator of the migration direction.

Line 431: Distribution of F is greater than stated here. Based on Colli et al 2015, haplogroup F goats have also been found in wild bezoar in Iran. The authors also report 1 modern Iranian domestic goat exhibiting the F haplogroup (Supplement Table 4). Colli et al 2015 also report that F is also present in wild bezoar from the Caucasus to Pakistan which I think is referencing data from Naderi et al 2008 also report that F is present in modern bezoar populations in southern Turkey and eastern Turkey as well.

Line 475: change wording of “more infrequent” to less frequent

Line 477: this statement is especially true since goats have higher reproductive potential than both sheep and cattle.

Line 521: reword “…while for now…”

Why are the mt haplogroups not labeled in the list of specimens in the Supplemental Tables?

Figure 2: the CKZ label is not identified on the map in Figure 1.

Figure 2: why do the 2 Capra sibirica samples not cluster together? Does this suggest hybridizing with domestic goat populations or just divergent wild Capra mt haplogroups?

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: Yes: Sandra L. Olsen

Reviewer #2: No

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files to be viewed.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org. Please note that Supporting Information files do not need this step.


11 Feb 2020

Thank you for the opportunity to revise our manuscript “High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor.” I would also like to thank the two reviewers for carefully reading this work and providing constructive criticism. I have applied all the changes suggested by you and reviewer #2, and below I will briefly summarize these modifications to the manuscript. (Reviewer #1 did not request corrections.)

I made changes to the formatting of the manuscript to conform to the journal’s style requirements, including the naming of attached files for figures and supporting information. I reviewed S1 Table and filled in missing specimen information, such as skeletal element identifications and archival location. Regarding Figure 1, which was flagged for displaying map imagery, I would like to point out that this imagery is in the public domain provided by naturalearthdata.org. I added a note about its origin and public domain status to the caption of Figure 1. For the supporting information materials, I added captions at the end of the manuscript that describe their contents.

In response to reviewer #2, I applied all suggested changes. Most substantively, I moved the descriptions of the archaeological sites to the supplementary materials (S1 Text) and replaced these with a concise table that includes basic information about the sites and in-text citations. The reviewer aptly noticed that the haplogroup designations of domesticated goats were not listed in S1 Table, and I have added these.

Reviewer #2’s suggestion to define mobile pastoralism in the introduction of the manuscript was good, and I have added a clause to the sentence describing the emergence of this lifeway in the Near East.

Reviewer #2’s idea to modify the sentence in the introduction about the spread of pastoralists into Inner Asia in order to reflect this phenomenon as a process rather than a homogenous event was on point. I made a simple change here to describe this as the spread of “pastoralist subsistence”, which could be either people moving or cultural transmissions of subsistence technologies.

Reviewer #2 noticed that we used a citation for the confusing skeletal morphology between ibex and domesticated sheep and goat that targets European and Near Eastern populations. This was intended to cover Siberian ibex. I agree that this was an improper citation, so I corrected the sentence to reflect a likelihood that Siberian ibex may also be confused with domesticated caprines, given what we know about these scenarios for European and Near Eastern populations.

Reviewer #2’s suggestion to move away from using “rare” to describe D lineage haplotypes of domesticated goats was good. I corrected this throughout the text, in addition to noting that C haplotypes are also notably rare in ancient distributions. Along these lines, I modified the text that describes F haplotypes in domesticated goats and ibexes following Reviewer #2’s suggestions.

Reviewer #2 kindly gave very minor suggestions for wording of a couple phrases , which I changed accordingly. I also corrected Figure 1 to display the “Central Kazakhstan” label for this analytical group.

Lastly, Reviewer #2 noticed that our ancient Capra sibirica sequences did not cluster with modern Capra sibirica sequences in the haplotype network shown in Figure 2. I added text to the section describing this species in the manuscript that suggests the disagreement is the result of divergence between modern and ancient populations. Hybridization between Siberian ibex and domesticated goats would not explain this observation, since mitochondrial DNA is inherited as a whole unit.

I would like to sincerely thank you and the reviewers for considering this manuscript for publication. I look forward to hearing from you soon.


5 May 2020

High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor

PONE-D-19-29815R1

Dear Dr. Hermes,

We are pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it complies with all outstanding technical requirements.

Within one week, you will receive an e-mail containing information on the amendments required prior to publication. When all required modifications have been addressed, you will receive a formal acceptance letter and your manuscript will proceed to our production department and be scheduled for publication.

Shortly after the formal acceptance letter is sent, an invoice for payment will follow. To ensure an efficient production and billing process, please log into Editorial Manager at https://www.editorialmanager.com/pone/, click the "Update My Information" link at the top of the page, and update your user information. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to enable them to help maximize its impact. If they will be preparing press materials for this manuscript, you must inform our press team as soon as possible and no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

With kind regards,

Peter F. Biehl, PhD

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

Reviewers' comments:


8 May 2020

PONE-D-19-29815R1

High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor

Dear Dr. Hermes:

I am pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please notify them about your upcoming paper at this point, to enable them to help maximize its impact. If they will be preparing press materials for this manuscript, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

For any other questions or concerns, please email plosone@plos.org.

Thank you for submitting your work to PLOS ONE.

With kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Peter F. Biehl

Academic Editor

PLOS ONE

This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

https://www.researchpad.co/tools/openurl?pubtype=article&doi=10.1371/journal.pone.0233333&title=High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor&author=&keyword=&subject=Research Article,Biology and Life Sciences,Organisms,Eukaryota,Animals,Vertebrates,Amniotes,Mammals,Ruminants,Goats,Biology and Life Sciences,Genetics,Heredity,Genetic Mapping,Haplotypes,Biology and Life Sciences,Organisms,Eukaryota,Animals,Animal Types,Domestic Animals,Biology and Life Sciences,Zoology,Animal Types,Domestic Animals,Biology and life sciences,Genetics,DNA,Forms of DNA,Mitochondrial DNA,Biology and life sciences,Biochemistry,Nucleic acids,DNA,Forms of DNA,Mitochondrial DNA,Biology and Life Sciences,Paleontology,Paleogenetics,Earth Sciences,Paleontology,Paleogenetics,Biology and Life Sciences,Evolutionary Biology,Population Genetics,Haplogroups,Biology and Life Sciences,Genetics,Population Genetics,Haplogroups,Biology and Life Sciences,Population Biology,Population Genetics,Haplogroups,People and Places,Geographical Locations,Asia,Biology and Life Sciences,Agriculture,Livestock,