Fossil and genomic evidence constrains the timing of bison arrival in North America.

The arrival of bison in North America marks one of the most successful large-mammal dispersals from Asia within the last million years, yet the timing and nature of this event remain poorly determined. Here, we used a combined paleontological and paleogenomic approach to provide a robust timeline for the entry and subsequent evolution of bison within North America. We characterized two fossil-rich localities in Canada's Yukon and identified the oldest well-constrained bison fossil in North America, a 130,000-y-old steppe bison, Bison cf. priscus We extracted and sequenced mitochondrial genomes from both this bison and from the remains of a recently discovered, ∼120,000-y-old giant long-horned bison, Bison latifrons, from Snowmass, Colorado. We analyzed these and 44 other bison mitogenomes with ages that span the Late Pleistocene, and identified two waves of bison dispersal into North America from Asia, the earliest of which occurred ∼195-135 thousand y ago and preceded the morphological diversification of North American bison, and the second of which occurred during the Late Pleistocene, ∼45-21 thousand y ago. This chronological arc establishes that bison first entered North America during the sea level lowstand accompanying marine isotope stage 6, rejecting earlier records of bison in North America. After their invasion, bison rapidly colonized North America during the last interglaciation, spreading from Alaska through continental North America; they have been continuously resident since then.

Early cave art and ancient DNA record the origin of European bison.

The two living species of bison (European and American) are among the few terrestrial megafauna to have survived the late Pleistocene extinctions. Despite the extensive bovid fossil record in Eurasia, the evolutionary history of the European bison (or wisent, Bison bonasus) before the Holocene (<11.7 thousand years ago (kya)) remains a mystery. We use complete ancient mitochondrial genomes and genome-wide nuclear DNA surveys to reveal that the wisent is the product of hybridization between the extinct steppe bison (Bison priscus) and ancestors of modern cattle (aurochs, Bos primigenius) before 120 kya, and contains up to 10% aurochs genomic ancestry. Although undetected within the fossil record, ancestors of the wisent have alternated ecological dominance with steppe bison in association with major environmental shifts since at least 55 kya. Early cave artists recorded distinct morphological forms consistent with these replacement events, around the Last Glacial Maximum (LGM, ∼21-18 kya).

Past climate changes, population dynamics and the origin of Bison in Europe.

BACKGROUND: Climatic and environmental fluctuations as well as anthropogenic pressure have led to the extinction of much of Europe's megafauna. The European bison or wisent (Bison bonasus), one of the last wild European large mammals, narrowly escaped extinction at the onset of the 20th century owing to hunting and habitat fragmentation. Little is known, however, about its origin, evolutionary history and population dynamics during the Pleistocene. RESULTS: Through ancient DNA analysis we show that the emblematic European bison has experienced several waves of population expansion, contraction, and extinction during the last 50,000 years in Europe, culminating in a major reduction of genetic diversity during the Holocene. Fifty-seven complete and partial ancient mitogenomes from throughout Europe, the Caucasus, and Siberia reveal that three populations of wisent (Bison bonasus) and steppe bison (B. priscus) alternately occupied Western Europe, correlating with climate-induced environmental changes. The Late Pleistocene European steppe bison originated from northern Eurasia, whereas the modern wisent population emerged from a refuge in the southern Caucasus after the last glacial maximum. A population overlap during a transition period is reflected in ca. 36,000-year-old paintings in the French Chauvet cave. Bayesian analyses of these complete ancient mitogenomes yielded new dates of the various branching events during the evolution of Bison and its radiation with Bos, which lead us to propose that the genetic affiliation between the wisent and cattle mitogenomes result from incomplete lineage sorting rather than post-speciation gene flow. CONCLUSION: The paleogenetic analysis of bison remains from the last 50,000 years reveals the influence of climate changes on the dynamics of the various bison populations in Europe, only one of which survived into the Holocene, where it experienced severe reductions in its genetic diversity. The time depth and geographical scope of this study enables us to propose temperate Western Europe as a suitable biotope for the wisent compatible with its reintroduction.

Genetic Analysis of the Henry Mountains Bison Herd.

Wild American plains bison (Bison bison) populations virtually disappeared in the late 1800s, with some remnant animals retained in what would become Yellowstone National Park and on private ranches. Some of these private bison were intentionally crossbred with cattle for commercial purposes. This forced hybridization resulted in both mitochondrial and nuclear introgression of cattle genes into some of the extant bison genome. As the private populations grew, excess animals, along with their history of cattle genetics, provided founders for newly established public bison populations. Of the US public bison herds, only those in Yellowstone and Wind Cave National Parks (YNP and WCNP) appear to be free of detectable levels of cattle introgression. However, a small free-ranging population (~350 animals) exists on public land, along with domestic cattle, in the Henry Mountains (HM) of southern Utah. This isolated bison herd originated from a founder group translocated from YNP in the 1940s. Using genetic samples from 129 individuals, we examined the genetic status of the HM population and found no evidence of mitochondrial or nuclear introgression of cattle genes. This new information confirms it is highly unlikely for free-living bison to crossbreed with cattle, and this disease-free HM bison herd is valuable for the long-term conservation of the species. This bison herd is a subpopulation of the YNP/WCNP/HM metapopulation, within which it can contribute significantly to national efforts to restore the American plains bison to more of its native range.

Hunting the Extinct Steppe Bison (Bison priscus) Mitochondrial Genome in the Trois-Frères Paleolithic Painted Cave.

Despite the abundance of fossil remains for the extinct steppe bison (Bison priscus), an animal that was painted and engraved in numerous European Paleolithic caves, a complete mitochondrial genome sequence has never been obtained for this species. In the present study we collected bone samples from a sector of the Trois-Frères Paleolithic cave (Ariège, France) that formerly functioned as a pitfall and was sealed before the end of the Pleistocene. Screening the DNA content of the samples collected from the ground surface revealed their contamination by Bos DNA. However, a 19,000-year-old rib collected on a rock apart the pathway delineated for modern visitors was devoid of such contaminants and reproducibly yielded Bison priscus DNA. High-throughput shotgun sequencing combined with conventional PCR analysis of the rib DNA extract enabled to reconstruct a complete mitochondrial genome sequence of 16,318 bp for the extinct steppe bison with a 10.4-fold coverage. Phylogenetic analyses robustly established the position of the Bison priscus mitochondrial genome as basal to the clade delineated by the genomes of the modern American Bison bison. The extinct steppe bison sequence, which exhibits 93 specific polymorphisms as compared to the published Bison bison mitochondrial genomes, provides an additional resource for the study of Bovinae specimens. Moreover this study of ancient DNA delineates a new research pathway for the analysis of the Magdalenian Trois-Frères cave.

Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: Application to the tribe Bovini (Cetartiodactyla, Bovidae).

Mitochondrial sequences are widely used for species identification and for studying phylogenetic relationships among closely related species or populations of the same species. However, many studies of mammals have shown that the maternal history of the mitochondrial genome can be discordant with the true evolutionary history of the taxa. In such cases, the analyses of multiple nuclear genes can be more powerful for deciphering interspecific relationships. Here, we designed primers for amplifying 13 new exon-primed intron-crossing (EPIC) autosomal loci for studying shallow phylogeny and taxonomy of Laurasiatherian mammals. Three criteria were used for the selection of the markers: gene orthology, a PCR product length between 600 and 1200 nucleotides, and different chromosomal locations in the bovine genome. Positive PCRs were obtained from different species representing the orders Carnivora, Cetartiodactyla, Chiroptera, Perissodactyla and Pholidota. The newly developed markers were analyzed in a phylogenetic study of the tribe Bovini (the group containing domestic and wild cattle, bison, yak, African buffalo, Asian buffalo, and saola) based on 17 taxa and 18 nuclear genes, representing a total alignment of 13,095 nucleotides. The phylogenetic results were compared to those obtained from analyses of the complete mitochondrial genome and Y chromosomal genes. Our analyses support a basal divergence of the saola (Pseudoryx) and a sister-group relationship between yak and bison. These results contrast with recent molecular studies but are in better agreement with morphology. The comparison of pairwise nucleotide distances shows that our nuDNA dataset provides a good signal for identifying taxonomic levels, such as species, genera, subtribes, tribes and subfamilies, whereas the mtDNA genome fails because of mtDNA introgression and higher levels of homoplasy. Accordingly, we conclude that the genus Bison should be regarded as a synonym of Bos, with the European bison relegated to a subspecies rank within Bos bison. We compared our molecular dating estimates to the fossil record in order to propose a biogeographic scenario for the evolution of Bovini during the Neogene.

Complete mitochondrial DNA sequence analysis of Bison bison and bison-cattle hybrids: function and phylogeny.

Complete mitochondrial DNA (mtDNA) genomes from 43 bison and bison-cattle hybrids were sequenced and compared with other bovids. Selected animals reflect the historical range and current taxonomic structure of bison. This study identified regions of potential nuclear-mitochondrial incompatibilities in hybrids, provided a complete mtDNA phylogenetic tree for this species, and uncovered evidence of bison population substructure. Seventeen bison haplotypes defined by 66 polymorphic sites were discovered, whereas 728 fixed differences and 86 non-synonymous mutations were identified between bison and bison-cattle hybrid sequences. The potential roles of the mtDNA genome in the function of hybrid animals and bison taxonomy are discussed.

Unraveling bovin phylogeny: accomplishments and challenges.

The phylogenetic systematics of bovin species forms a common basis for studies at multiple scales, from the level of domestication in populations to major cladogenesis. The main big-picture accomplishments of this productive field, including two recent works, one in BMC Genomics, are reviewed with an eye for some of the limitations and challenges impeding progress.

Phylogenetic position of the saola (Pseudoryx nghetinhensis) inferred from cytogenetic analysis of eleven species of Bovidae.

Previous morphological and molecular analyses failed to resolve the phylogenetic position of the critically endangered saola (Pseudoryx nghetinhensis) with respect to its placement in Bovina (cattle, bison, and yak) or Bubalina (Asian and African buffaloes). In the present study, G- and C-banding, Ag-staining and FISH with 28S and telomeric probes was undertaken for 17 bovid species. An analysis of these data allowed us to identify 49 structural rearrangements that included autosomes, gonosomes and 17 different NOR sites. The combined data set was subjected to a cladistic analysis aimed at: (i) providing new insights on phylogenetic relationships of the saola and other species within the subfamily Bovinae, and (ii) testing the suitability of different classes of chromosomal characters for phylogenetic reconstruction of the family Bovidae. The study revealed that nucleolar organizing regions (NORs) are phylogenetically informative. It was shown that at least one, or sometimes two of these characters punctuate divergences that include nodes that are the most basal in the tree, to those that are the most recent. In this context, the shared presence of three NORs in saola and species of Syncerus and Bubalus strongly suggests the saola's placement within the subtribe Bubalina. This contrasts with Robertsonian rearrangements which are informative only at the generic level. These findings suggest that NORs are an important and frequently overlooked source of additional phylogenetic information within the Bovidae that may also have applicability at higher taxonomic levels, possibly even for Pecora.

Rise and fall of the Beringian steppe bison.

The widespread extinctions of large mammals at the end of the Pleistocene epoch have often been attributed to the depredations of humans; here we present genetic evidence that questions this assumption. We used ancient DNA and Bayesian techniques to reconstruct a detailed genetic history of bison throughout the late Pleistocene and Holocene epochs. Our analyses depict a large diverse population living throughout Beringia until around 37,000 years before the present, when the population's genetic diversity began to decline dramatically. The timing of this decline correlates with environmental changes associated with the onset of the last glacial cycle, whereas archaeological evidence does not support the presence of large populations of humans in Eastern Beringia until more than 15,000 years later.

Maternal and paternal lineages in cross-breeding bovine species. Has wisent a hybrid origin?

The tribe Bovini comprises cattle and cattle-like species. Reconstructions of their phylogeny have so far been incomplete and have yielded conflicting conclusions about the relationship of American bison and wisent (European bison). We have compared the sequences of three mitochondrial and two Y-chromosomal DNA segments. Mitochondrial DNA indicates that four distinct maternal lineages diverged after an early split-off of the buffalo species, leading to (1) taurine cattle and zebu, (2) wisent, (3) American bison and yak, and (4) banteng, gaur, and gayal, respectively. At a higher level, lineages (1) and (2) and lineages (3) and (4) are probably associated. In contrast, Y-chromosomal sequences indicate a close association of American and European bison, which is in agreement with their morphological similarity, complete fertility of hybrid offspring, and amplified fragment length polymorphism (AFLP) fingerprints of nuclear DNA. One explanation for the anomalous divergence of the mitochondrial DNA from the two bison species is lineage sorting, which implies that two distinct mitochondrial lineages coexisted in the bison-yak branch until the recent divergence of American bison and wisent. Alternatively, the wisent may have emerged by species hybridization initiated by introgression of bison bulls in another ancestral species. This "transpatric" mode of species formation would be consistent with the recent appearance of the wisent in the fossil record without clearly identifiable ancestors.

Phylogeny of bovine species based on AFLP fingerprinting.

The Bovini species comprise both domestic and wild cattle species. Published phylogenies of this tribe based on mitochondrial DNA contain anomalies, while nuclear sequences show only low variation. We have used amplified fragment length polymorphism (AFLP) fingerprinting in order to detect variation in loci distributed over the nuclear genome. Computer-assisted scoring of electrophoretic fingerprinting patterns yielded 361 markers, which provided sufficient redundancy to suppress stochastic effects of intraspecies polymorphisms and length homoplasies (comigration of non-homologous fragments). Tree reconstructions reveal three clusters: African buffalo with water buffalo, ox with zebu, and bison with wisent. Similarity values suggest a clustering of gaur and banteng, but bifurcating clustering algorithms did not assign consistent positions to these species and yak. We propose that because of shared polymorphisms and reticulations, tree topologies are only partially adequate to represent the phylogeny of the Bovini. Principal-coordinate analysis positions zebu between a gaur/banteng cluster and taurine cattle. This correlates with the region of origin of these species and suggests that genomic distances between the cattle species have been influenced by genetic exchange between neighbouring ancestral populations.

Genetic variation within and relatedness among wood and plains bison populations.

There are two recognized subspecies of bison, wood (Bison bison athabascae) and plains (Bison bison bison) bison. The establishment of most bison populations from a small number of individuals has raised concerns about their genetic variation. To this end, 11 bison populations were surveyed with 11 microsatellite loci in order to calculate genetic variation and genetic distances. Mean number of alleles ranged between 3.18 at Antelope Island State Park (Utah) and 6.55 at Wood Buffalo National Park (Alberta and Northwest Territories). Mean heterozygosity ranged from 0.295 at Antelope Island State Park to 0.669 at Custer State Park (South Dakota). The amount of genetic variability present in the bison populations as measured by mean number of alleles and overall probability of identity was found to correlate with the number of founders for all sampled populations. The G-test for heterogeneity revealed some evidence for the existence of subpopulations at Wood Buffalo National Park, however very small genetic distances between these subpopulations suggest that nuclear material from the plains bison introduced into Wood Buffalo National Park has diffused throughout the park. Genetic distances between the sampled populations were generally larger between than within the two bison subspecies.

[The estimation of the degree of genetic divergence in Bovidae by DNA restrictase analysis]. / Otsenka stepeni geneticheskoi divergentsii nekotorykh polorogikh metodom restriktanogo analiza DNK.

Restrictase DNA analysis was performed in seven species of Bovidae family--European (Belovezh) bison, Belovezh-Caucasian bison, European-American bison, mountain Caucasian bison, American bison, domestic bull, and domestic sheep. DNA was digested with MspI and TaqI restriction endonucleases. The products obtained were end-labeled and electrophoretically separated in poliacrilamid gel. DNA fragments consisting of highly repetitive genomic sequences were detected as a set of bands corresponding to fragments between 464 and 67 bp in length. All forms of Bison genus analysed were characterized by the identical sets of bands. The band patterns of bulls is similar to these of the bison, whereas the patterns of sheep is rather different. These data indicate that American bison and European bison geographically isolated forms of the same polytypic species.