This is an open access article under the terms of the Creat ive Commo ns Attri bution-NonCo mmercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Background The contemporary Italian wolf (Canis lupus italicus) represents a case of morphological and genetic uniqueness. Today, Italian wolves are also the only documented population to fall exclusively within the mitochondrial haplogroup 2, which was the most diffused across Eurasian and North American wolves during the Late Pleistocene. However, the dynamics leading to such distinctiveness are still debated. Methods In order to shed light on the ancient genetic variability of this wolf population and on the origin of its current diversity, we collected 19 Late Pleistocene-Holocene samples from northern Italy, which we analyzed at a short portion of the hypervariable region 1 of the mitochondrial DNA, highly informative for wolf and dog phylogenetic analyses. Results Four out of the six detected haplotypes matched the ones found in ancient wolves from northern Europe and Beringia, or in modern European and Chinese wolves, and appeared closely related to the two haplotypes currently found in Italian wolves. The haplotype of two Late Pleistocene samples matched with primitive and contemporary dog sequences from the canine mitochondrial clade A. All these haplotypes belonged to haplogroup 2. The only exception was a Holocene sample dated 3,250 years ago, affiliated to haplogroup 1. Discussion In this study we describe the genetic variability of the most ancient wolf specimens from Italy analyzed so far, providing a preliminary overview of the genetic make-up of the population that inhabited this area from the last glacial maximum to the Middle Age period. Our results endorsed that the genetic diversity carried by the Pleistocene wolves here analyzed showed a strong continuity with other northern Eurasian wolf specimens from the same chronological period. Contrarily, the Holocene samples showed a greater similarity only with modern sequences from Europe and Asia, and the occurrence of an haplogroup 1 haplotype allowed to date back previous finding about its presence in this area. Moreover, the unexpected discovery of a 24,700-year-old sample carrying a haplotype that, from the fragment here obtained, falls within the canine clade A, could represent the oldest evidence in Europe of such dog-rich clade. All these findings suggest complex population dynamics that deserve to be further investigated based on mitochondrial or whole genome sequencing.

Summary The evolution of the genera Bos and Bison , and the nature of gene flow between wild and domestic species, is poorly understood, with genomic data of wild species being limited. We generated two genomes from the likely extinct kouprey ( Bos sauveli ) and analyzed them alongside other Bos and Bison genomes. We found that B. sauveli possessed genomic signatures characteristic of an independent species closely related to Bos javanicus and Bos gaurus . We found evidence for extensive incomplete lineage sorting across the three species, consistent with a polytomic diversification of the major ancestry in the group, potentially followed by secondary gene flow. Finally, we detected significant gene flow from an unsampled Asian Bos -like source into East Asian zebu cattle, demonstrating both that the full genomic diversity and evolutionary history of the Bos complex has yet to be elucidated and that museum specimens and ancient DNA are valuable resources to do so.

During historical times many local grey wolf (Canis lupus) populations underwent a substantial reduction of their sizes or became extinct. Among these, the wolf population once living in Sicily, the biggest island of the Mediterranean Sea, was completely eradicated by human persecution in the early decades of the XX century. In order to understand the genetic identity of the Sicilian wolf, we applied ancient DNA techniques to analyse the mitochondrial DNA of six specimens actually stored in Italian museums. We successfully amplified a diagnostic mtDNA fragment of the control region (CR) in four of the samples. Results showed that two samples shared the same haplotype, that differed by two substitutions from the currently most diffused Italian wolf haplotype (W14) and one substitution from the only other Italian haplotype (W16). The third sample showed a wolf-like haplotype never described before and the fourth a haplotype commonly found in dogs. Furthermore, all the wolf haplotypes detected in this study belonged to the mitochondrial haplogroup that includes haplotypes detected in all the known European Pleistocene wolves and in several modern southern European populations. Unfortunately, this endemic island population, bearing unique mtDNA variability, was definitively lost before it was possible to understand its taxonomic uniqueness and conservational value.

Large vertebrates are extremely sensitive to anthropogenic pressure, and their populations are declining fast. The white rhinoceros (Ceratotherium simum) is a paradigmatic case: this African megaherbivore suffered a remarkable population reduction in the last 150 years due to human activities. The two white rhinoceros subspecies, the northern (NWR) and the southern white rhinoceros (SWR), however, underwent opposite fates: the NWR vanished quickly after the onset of the decline, while the SWR recovered after a severe bottleneck. Such demographic events are predicted to have an erosive effect at the genomic level, in connection with the extirpation of diversity, and increased genetic drift and inbreeding. However there is currently little empirical data available that allows us to directly reconstruct the subtleties of such processes in light of distinct demographic histories. Therefore to assess these effects, we generated a whole-genome, temporal dataset consisting of 52 re-sequenced white rhinoceros genomes, that represents both subspecies at two time windows: before and during/after the bottleneck. Our data not only reveals previously unknown population substructure within both subspecies, but allowed us to quantify the genomic erosion undergone by both, with post-bottleneck white rhinoceroses harbouring significantly fewer heterozygous sites, and showing higher inbreeding coefficients than pre-bottleneck individuals. Moreover, the effective population size suffered a decrease of two and three orders of magnitude in the NWR and SWR respectively, due to the recent bottleneck. Our data therefore provides much needed empirical support for theoretical predictions about the genomic consequences of shrinking populations, information that is relevant for understanding the process of population extinction. Furthermore, our findings have the potential to inform management approaches for the conservation of the remaining white rhinoceroses.

Dog domestication is still largely unresolved due to time-gaps in the sampling of regions. Ancient Italian canids are particularly understudied, currently represented by only a few specimens. In the present study, we sampled 27 canid remains from Northern Italy dated between the Late Pleistocene and Bronze Age to assess their genetic variability, and thus add context to dog domestication dynamics. They were targeted at four DNA fragments of the hypervariable region 1 of mitochondrial DNA. A total of 11 samples had good DNA preservation and were used for phylogenetic analyses. The dog samples were assigned to dog haplogroups A, C and D, and a Late Pleistocene wolf was set into wolf haplogroup 2. We present our data in the landscape of ancient and modern dog genetic variability, with a particular focus on the ancient Italian samples published thus far. Our results suggest there is high genetic variability within ancient Italian canids, where close relationships were evident between both a ~24,700 years old Italian canid, and Iberian and Bulgarian ancient dogs. These findings emphasize that disentangling dog domestication dynamics benefits from the analysis of specimens from Southern European regions.

No abstract

The Sicilian wolf represented the only population of wolves living on a Mediterranean island until the first half of the twentieth century (1930s-1960s). Previous studies hypothesised that they remained isolated from mainland wolves from the end of the Last Glacial Maximum (LGM), until human persecutions led them to extinction. There are only seven known Sicilian wolf specimens from the 19th and 20th century preserved in museums in Italy and recent morphometric analyses assigned them to the new subspecies Canis lupus cristaldii 10. To better understand the origins of the Sicilian wolf, and its relationship to other wolf populations, we sequenced four whole genomes (3.8x-11.6x) and five mitogenomes. We investigated the relationship between Sicilian wolves and other modern breeds to identify potential admixture. Furthermore, considering that the last land-bridge between Sicily and Italy disappeared after the LGM, around 17 kya, we explored the possibility that the Sicilian wolf retained ancestry from ancient wolf and dog lineages. Additionally, we explored whether the long-term isolation might have affected the genomic diversity, inbreeding levels and genetic load of the Sicilian wolf. Our findings show that the Sicilian wolves shared most ancestry with the modern Italian wolf population but are better modelled as admixed with European dog breeds, and shared traces of Eneolithic and Bronze age European dogs. We also find signatures of severe inbreeding and low genomic diversity at population and individual levels due to long-term isolation and drift, suggesting also low effective population size.