Whole‐genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats

Whole‐genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.Deutscher Akademischer Austauschdienst http://dx.doi.org/10.13039/501100001655Javna Agencija za Raziskovalno Dejavnost RS http://dx.doi.org/10.13039/501100004329Ministry of Science and Higher Education of the Russian Federation http://dx.doi.org/10.13039/501100012190Peer Reviewe

Range shifts or extinction? Ancient DNA and distribution modelling reveal past and future responses to climate warming in cold-adapted birds.

Global warming is predicted to cause substantial habitat rearrangements, with the most severe effects expected to occur in high-latitude biomes. However, one major uncertainty is whether species will be able to shift their ranges to keep pace with climate-driven environmental changes. Many recent studies on mammals have shown that past range contractions have been associated with local extinctions rather than survival by habitat tracking. Here, we have used an interdisciplinary approach that combines ancient DNA techniques, coalescent simulations and species distribution modelling, to investigate how two common cold-adapted bird species, willow and rock ptarmigan (Lagopus lagopus and Lagopus muta), respond to long-term climate warming. Contrary to previous findings in mammals, we demonstrate a genetic continuity in Europe over the last 20 millennia. Results from back-casted species distribution models suggest that this continuity may have been facilitated by uninterrupted habitat availability and potentially also the greater dispersal ability of birds. However, our predictions show that in the near future, some isolated regions will have little suitable habitat left, implying a future decrease in local populations at a scale unprecedented since the last glacial maximum

Nonreceding hare lines: genetic continuity since the Late Pleistocene in European mountain hares (Lepus timidus)

Throughout time, climate changes have caused substantial rearrangements of habitats which have alternately promoted and disfavoured different types of taxa. At first glance, the mountain hare (Lepus timidus) shows the typical hallmarks of a cold-adapted species that has retreated to refugia since the onset of the current Holocene interglacial. In contrary to expectations, however, the species has a high contemporary genetic diversity with no clear differentiation between geographically isolated populations. In order to clarify the phylogeographic history of European mountain hares, we here analysed ancient DNA from the glacial populations that inhabited the previous midlatitude European tundra region. Our results reveal that the Ice Age hares had similar levels of genetic variation and lack of geographic structure as observed today, and the ancient samples were intermingled with modern individuals throughout the reconstructed evolutionary tree. This suggest a temporal genetic continuity in Europe, where the mountain hares were able to keep pace with the rapid changes at the last glacial/interglacial transition, and successfully track their shifting habitat to northern and alpine regions. Further, the temporal demographic analyses showed that the species’ population size in Europe appear to have been tightly linked with palaeoclimatic fluctuations, with increases and declines occurring during periods of global cooling and warming, respectively. Taken together, our results suggest that neither habitat shifts nor demographic fluctuations have had any substantial impact on the genetic diversity of European mountain hares. This remarkable resilience, which contrasts to a majority of previously investigated cold-adapted species, is likely due to its generalist nature which makes it less vulnerable to environmental changes

Spotted phenotypes in horses lost attractiveness in the Middle Ages

Horses have been valued for their diversity of coat colour since prehistoric times; this is especially the case since their domestication in the Caspian steppe in ~3,500 BC. Although we can assume that human preferences were not constant, we have only anecdotal information about how domestic horses were influenced by humans. Our results from genotype analyses show a significant increase in spotted coats in early domestic horses (Copper Age to Iron Age). In contrast, medieval horses carried significantly fewer alleles for these phenotypes, whereas solid phenotypes (i.e., chestnut) became dominant. This shift may have been supported because of (i) pleiotropic disadvantages, (ii) a reduced need to separate domestic horses from their wild counterparts, (iii) a lower religious prestige, or (iv) novel developments in weaponry. These scenarios may have acted alone or in combination. However, the dominance of chestnut is a remarkable feature of the medieval horse population.Peer Reviewe

Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries)

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Aşıklı Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Aşıklı Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Aşıklı Höyük’s occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations

Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries)

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Aşıklı Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Aşıklı Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Aşıklı Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations

Demographic History of Indigenous Populations in Mesoamerica Based on mtDNA Sequence Data

The genetic characterization of Native American groups provides insights into their history and demographic events. We sequenced the mitochondrial D-loop region (control region) of 520 samples from eight Mexican indigenous groups. In addition to an analysis of the genetic diversity, structure and genetic relationship between 28 Native American populations, we applied Bayesian skyline methodology for a deeper insight into the history of Mesoamerica. AMOVA tests applying cultural, linguistic and geographic criteria were performed. MDS plots showed a central cluster of Oaxaca and Maya populations, whereas those from the North and West were located on the periphery. Demographic reconstruction indicates higher values of the effective number of breeding females (Nef) in Central Mesoamerica during the Preclassic period, whereas this pattern moves toward the Classic period for groups in the North and West. Conversely, Nef minimum values are distributed either in the Lithic period (i.e. founder effects) or in recent periods (i.e. population declines). The Mesomerican regions showed differences in population fluctuation as indicated by the maximum Inter-Generational Rate (IGRmax): i) Center-South from the lithic period until the Preclassic; ii) West from the beginning of the Preclassic period until early Classic; iii) North characterized by a wide range of temporal variation from the Lithic to the Preclassic. Our findings are consistent with the genetic variations observed between central, South and Southeast Mesoamerica and the North-West region that are related to differences in genetic drift, structure, and temporal survival strategies (agriculture versus hunter-gathering, respectively). Interestingly, although the European contact had a major negative demographic impact, we detect a previous decline in Mesoamerica that had begun a few hundred years before

CoalDist: An application to infer a population's demographic history from coalescent times of genetic data

<p><span>CoalMCMC is a software application for inferring a natural population demographic history out of data extracted from a coalescent genealogy in turn inferred by other means such as an intra-specific phylogenetic reconstruction. The analysis requires a list of coalescent times, samples ages, and a population design. The ultimate purpose of CoalMCMC is, however, not the demographic history <em>per se</em>, for which more sophisticated methods exist, such as Bayesian Skyline Plots but to assist in the interpretation of results, specifically the effects that genetic structure and temporal sampling patterns have in the demographic inference. The rationale to do this is that, for instance, genetic structure or a temporal agregation of samples may produce spurious signals of a population growth or a demogaphic bottleneck, respectively. If such artifacts are suspected, the inference made with CoalMCMC can tell if this was the case. For instance, if CoalMCMC inference, which takes into account the genetic structure or temporal sampling agregation, infers a flat demography (no growth or no bottleneck), then such events inferred by other methods can be assumed spurious. </span></p> <p><span>The reason for not relying only on a Skyline Plot for demographic inference is that Skyline Plots implement smoothing algorithms that mask the presence of sudden changes, and such a complex inferential procedure doesn't allow to know if the coalescent times contain information suggesting e.g. a bottleneck or not.  </span></p> <p><span>The input of CoalMCMC partially follows the input design of the simulations module of the software BaySICS (10.5281/zenodo.10210929, </span><span><a href="https://doi.org/10.1371/journal.pone.0098011"><span>https://doi.org/10.1371/journal.pone.0098011</span></a></span><span>). The design of BaySICS is intended for defining coalescent simulations with flexible scenarios that include arbitrary definitions of historical genetic structuring (including split and merging of populations), migration, and change of effective population sizes (N<sub>e</sub>). Although the primary use of CoalMCMC is to infer historical demography, additional elements can be inferred as well. The inferred data is the set of parameter values, established as priors in the input file (and initially sampled from their prior probability distributions). </span></p&gt