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.
The Sardinian dhole (Cynotherium sardous)1 was an iconic and unique canid species that was endemic to Sardinia and Corsica until it became extinct at the end of the Late Pleistocene.2-5 Given its peculiar dental morphology, small body size, and high level of endemism, several extant canids have been proposed as possible relatives of the Sardinian dhole, including the Asian dhole and African hunting dog ancestor.3,6-9 Morphometric analyses3,6,8-12 have failed to clarify the evolutionary relationship with other canids.We sequenced the genome of a ca-21,100-year-old Sardinian dhole in order to understand its genomic history and clarify its phylogenetic position. We found that it represents a separate taxon from all other living canids from Eurasia, Africa, and North America, and that the Sardinian dhole lineage diverged from the Asian dhole ca 885 ka. We additionally detected historical gene flow between the Sardinian and Asian dhole lineages, which ended approximately 500-300 ka, when the land bridge between Sardinia and mainland Italy was already broken, severing their population connectivity. Our sample showed low genome-wide diversity compared to other extant canids-probably a result of the long-term isolation-that could have contributed to the subsequent extinction of the Sardinian dhole.
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 has suffered a remarkable decline in the last 150 years due to human activities. Its subspecies, the northern (NWR) and the southern white rhinoceros (SWR), however, underwent opposite fates: the NWR vanished quickly, while the SWR recovered after the severe decline. Such demographic events are predicted to have an erosive effect at the genomic level, linked to the extirpation of diversity, and increased genetic drift and inbreeding. However, there is little empirical data available to directly reconstruct the subtleties of such processes in light of distinct demographic histories. Therefore, we generated a whole-genome, temporal data set consisting of 52 resequenced white rhinoceros genomes, representing both subspecies at two time windows: before and during/after the bottleneck. Our data reveal previously unknown population structure within both subspecies, as well as quantifiable genomic erosion. Genome-wide heterozygosity decreased significantly by 10% in the NWR and 36% in the SWR, and inbreeding coefficients rose significantly by 11% and 39%, respectively. Despite the remarkable loss of genomic diversity and recent inbreeding it suffered, the only surviving subspecies, the SWR, does not show a significant accumulation of genetic load compared to its historical counterpart. Our data provide empirical support for predictions about the genomic consequences of shrinking populations, and our findings have the potential to inform the conservation efforts of the remaining white rhinoceroses.