Conservation genomics: disequilibrium mapping of domestic cattle chromosomal segments in North American bison populations

Introgressive hybridization is one of the major threats to species conservation, and is often induced by human influence on the natural habitat of wildlife species. The ability to accurately identify introgression is critical to understanding its importance in evolution and effective conservation management of species. Hybridization between North American bison (Bison bison) and domestic cattle (Bos taurus) as a result of human activities has been recorded for over 100 years, and domestic cattle mitochondrial DNA was previously detected in bison populations. In this study, linked microsatellite markers were used to identify domestic cattle chromosomal segments in 14 genomic regions from 14 bison populations. Cattle nuclear introgression was identified in five populations, with an average frequency per population ranging from 0.56% to 1.80%. This study represents the first use of linked molecular markers to examine introgression between mammalian species and the first demonstration of domestic cattle nuclear introgression in bison. To date, six public bison populations have been identified with no evidence of mitochondrial or nuclear domestic cattle introgression, providing information critical to the future management of bison genetic resources. The ability to identify even low levels of introgression resulting from historic hybridization events suggests that the use of linked molecular markers to identify introgression is a significant development in the study of introgressive hybridization across a broad range of taxa.     

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.     

Differential introgression of uniparentally inherited markers in bison populations with hybrid ancestries

Historical hybridization between Bison bison (bison) and Bos taurus (cattle) has been well documented and resulted in cattle mitochondrial DNA (mtDNA) introgression, previously identified in six different bison populations. In order to examine Y chromosome introgression, a microsatellite marker (BYM-1) with non-overlapping allele size distributions in bison and cattle was isolated from a bacterial artificial chromosome (BAC) clone, and was physically assigned to the Y chromosome by fluorescence in situ hybridization. BYM-1 genotypes for a sample of 143 male bison from 10 populations, including all six populations where cattle mtDNA haplotypes were previously identified, indicated that cattle Y chromosome introgression had not occurred in these bison populations. The differential permeability of uniparentally inherited markers to introgression is consistent with observations of sterility among first generation hybrid males and a sexual asymmetry in the direction of hybridization favouring matings between male bison and female cattle.     

Paternally inherited markers in bovine hybrid populations

The genetic integrity of crossfertile bovine- or cattle-like species may be endangered by species hybridization. Previously, amplified fragment length polymorphism, satellite fragment length polymorphism and microsatellite assays have been used to analyze the species composition of nuclear DNA in taurine cattle, zebu, banteng and bison populations, while mitochondrial DNA reveals the origin of the maternal lineages. Here, we describe species-specific markers of the paternally transmitted Y-chromosome for the direct detection of male-mediated introgression. Convenient PCR-restriction fragment length polymorphism and competitive PCR assays are shown to differentiate the Y-chromosomes of taurine cattle, American bison and European bison, and to detect the banteng origin of Indonesian Madura and Bali cattle bulls.     

Bison PRNP genotyping and potential association with Brucella spp. seroprevalence

The implication that host cellular prion protein (PrP(C)) may function as a cell surface receptor and/or portal protein for Brucella abortus in mice prompted an evaluation of nucleotide and amino acid variation within exon 3 of the prion protein gene (PRNP) for six US bison populations. A non-synonymous single nucleotide polymorphism (T50C), resulting in the predicted amino acid replacement M17T (Met --> Thr), was identified in each population. To date, no variation (T50; Met) has been detected at the corresponding exon 3 nucleotide and/or amino acid position for domestic cattle. Notably, 80% (20 of 25) of the Yellowstone National Park bison possessing the C/C genotype were Brucella spp. seropositive, representing a significant (P = 0.021) association between seropositivity and the C/C genotypic class. Moreover, significant differences in the distribution of PRNP exon 3 alleles and genotypes were detected between Yellowstone National Park bison and three bison populations that were either founded from seronegative stock or previously subjected to test-and-slaughter management to eradicate brucellosis. Unlike domestic cattle, no indel polymorphisms were detected within the corresponding regions of the putative bison PRNP promoter, intron 1, octapeptide repeat region or 3'-untranslated region for any population examined. This study provides the first evidence of a potential association between nucleotide variation within PRNP exon 3 and the presence of Brucella spp. antibodies in bison, implicating PrP(C) in the natural resistance of bison to brucellosis infection.     

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.     

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.     

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.     

Phylogenies using mtDNA and SRY provide evidence for male-mediated introgression in Asian domestic cattle

Using nucleotide sequences of the mitochondrial DNA (mtDNA) cytochrome b and SRY genes, we examined the genetic status of two major groups of domestic cattle, the humpless taurine (Bos taurus) and humped zebu (B. indicus), using 10 cattle populations in Asia. Several sequence polymorphisms specific for each major group were found, although the frequency of these polymorphisms varied in each population. Six major mtDNA-SRY composite types were observed. The Mishima, Mongolian, Korean, Chinese Yellow and Sri Lanka cattle populations had a full match between the mtDNA and SRY sequences, specifically the taurine/taurine type or zebu/zebu type. A non-match type (zebu/taurine type) was found at a high frequency in the Bangladesh (83.4%) and Nepal populations (83.3%). Our results suggest that these non-match type populations developed from genetic hybridization of different strains. Also, the domestication history of modern Asian domestic cattle could be explained by male-mediated introgression. Additionally, our results suggest the occurrence of introgression of mtDNA from other Bibos or Poephagus species into native cattle populations. The existence of other mtDNA-SRY composite types, such as the Bali-zebu and yak-zebu types in Indonesia (85.7%) and Nepal (16.7%), respectively, suggests that genetic introgression also occurred from other genera into domestic cattle during the process of domestication.     

Bidirectional mitochondrial introgression between Korean cobitid fish mediated by hybridogenetic hybrids

Genomic introgression through interspecific hybridization has been observed in some species of the freshwater fish family Cobitidae. Within this family, a Cobitis hankugensis–Iksookimia longicorpa diploid–triploid hybrid species complex on the Korean peninsula is unique in displaying hybridogenesis, a unisexual reproduction mode that allows hybrids to mediate the transfer of mitochondrial DNA (but not nuclear DNA) between the two parent species. However, populations of the parental species in the wild have never been examined for the potential effect of introgression on their genomes. To address the genetic consequences of unisexual hybridization on the parental species, we examined genetic structure of the two parental species, C. hankugensis and I. longicorpa, in three independent natural habitats where they coexist with their hybrid complex using DNA sequence data of one mitochondrial gene and three nuclear genes. We found that mitochondrial introgression between the two species was extensive in all the examined localities, while there was no evidence of nuclear introgression across the species boundary. This result indicates that the hybridogenetic individuals mediate mitochondrial introgression from one species to the other, producing mito‐nuclear mosaic genomes such as C. hankugensis nuclear genomes associated with I. longicorpa mitochondrial DNA and the reverse. The direction and degree of introgression varied among the three localities, but the underlying mechanisms for this observation proved elusive. Introgression might depend on which species serves as the predominant sperm or ovum donor or the environmental conditions of the localities. The present study suggests that introgressive hybridization between pure C. hankugensis and I. longicorpa species is highly likely where the two species co‐occur with hybridogenetic individuals, but the consequence of introgression could be variable due to the history and environmental characteristics of particular populations across the parental species’ ranges.     

Cattle ancestry in bison: explanations for higher mtDNA than autosomal ancestry

Understanding and documenting the process of hybridization and introgression between related species is a major focus of recent evolutionary research using molecular techniques. Many North American bison herds have cattle ancestry introduced by crossbreeding over a century ago. Molecular estimates of this ancestry have shown much higher levels for cattle mtDNA than for autosomal cattle genes. A large part of this difference appears to be the result of partial reproductive isolation between the two species where only bison bull × domestic cow crosses are successful, and all the surviving progeny are females. In addition, selection against autosomal cattle genes in bison may have contributed to differential levels of cattle ancestry. The impact of selection against cattle mtDNA and gene flow of bison mtDNA are examined to explain particular combinations of mtDNA and autosomal cattle ancestry. A bottleneck, after the level of cattle ancestry in bison was reduced to a low level, is consistent with the high variance over autosomal loci observed for cattle ancestry, and differential selection among cattle loci in bison does not need to be invoked. Further examination of the cattle genome in bison may shed light on whether these markers, or their associated regions, are indeed neutral.     

The influence of nuclear content on developmental competence of gaur x cattle hybrid in vitro fertilized and somatic cell nuclear transfer embryos

In nondomestic and endangered species, the use of domestic animal oocytes as recipients for exotic donor nuclei causes the normal pattern of cytoplasmic inheritance to be disrupted, resulting in the production of nuclear-cytoplasmic hybrids. Evidence suggests that conflict between nuclear and cytoplasmic control elements leads to a disruption of normal cellular processes, including metabolic function and cell division. This study investigated the effects of nuclear-cytoplasmic interactions on the developmental potential of interspecies embryos produced by in vitro fertilization and somatic cell nuclear transfer: cattle x cattle, gaur x cattle, hybrid x cattle. Cattle control and hybrid embryos were examined for development to the blastocyst stage and blastocyst quality, as determined by cell number and allocation, apoptosis incidence, and expression patterns of mitochondria-related genes. These analyses demonstrated that a 100% gaur nucleus within a domestic cattle cytoplasmic environment was not properly capable of directing embryo development in the later preimplantation stages. Poor blastocyst development accompanied by developmental delay, decreased cell numbers, and aberrant apoptotic and related gene expression profiles, all signs of disrupted cellular processes associated with mitochondrial function, were observed. Developmental potential was improved when at least a portion of the nuclear genome corresponded to the inherited cytoplasm, indicating that recognition of cytoplasmic components by the nucleus is crucial for proper cellular function and embryo development. A better understanding of the influence of the cytoplasmic environment on embryonic processes is necessary before interspecies somatic cell nuclear transfer can be considered a viable alternative for endangered species conservation.     

Ruminal ciliated protozoa in bison.

Ruminal contents from 79 slaughtered bison and 2 ruminally cannulated bison were collected to obtain information on total numbers and species distribution of ciliated protozoa. The bison originated from numerous herds throughout the Great Plains and were grouped into three dietary categories: (i) only forage; (ii) forage with moderate levels of supplementation; and (iii) feedlot concentrate-silage diet. Total ciliate counts were highest in bison receiving grain supplementation (210.1 x 10(4)/g) and lowest in bison consuming only forage (27.1 x 10(4)/g). All protozoan species found in bison have been reported in domestic livestock, although Ophryoscolex sp., a relatively common protozoan in cattle, was detected at low concentrations in only eight bison. The uncommon holotrich Microcetus lappus was present in five bison in concentrations reaching 8.4% of the total ciliate population. Charonina ventriculi, another infrequently observed species, was present in 18 bison, with the highest concentrations in forage-fed animals. Thirty bison possessed a type B protozoan population, characterized by Epidinium sp., Eudiplodinium maggii, and Eudiplodinium bovis. Thirty-eight bison possessed a mixed A-B population, characterized by Polyplastron sp. coexisting with low numbers of Eudiplodinium maggii or Epidinium sp. or both. Thirteen bison possessed populations lacking any remnant type B ciliate species. At least 29 of the bison possessing Polyplastron sp. were known to have been in contact with cattle, whereas all bison isolated from cattle had type B populations. The reduction of type B populations in bison becomes increasingly likely as bison production expands into areas inhabited by domestic livestock.     

Ruminal ciliated protozoa in bison

Ruminal contents from 79 slaughtered bison and 2 ruminally cannulated bison were collected to obtain information on total numbers and species distribution of ciliated protozoa. The bison originated from numerous herds throughout the Great Plains and were grouped into three dietary categories: (i) only forage; (ii) forage with moderate levels of supplementation; and (iii) feedlot concentrate-silage diet. Total ciliate counts were highest in bison receiving grain supplementation (210.1 x 10(4)/g) and lowest in bison consuming only forage (27.1 x 10(4)/g). All protozoan species found in bison have been reported in domestic livestock, although Ophryoscolex sp., a relatively common protozoan in cattle, was detected at low concentrations in only eight bison. The uncommon holotrich Microcetus lappus was present in five bison in concentrations reaching 8.4% of the total ciliate population. Charonina ventriculi, another infrequently observed species, was present in 18 bison, with the highest concentrations in forage-fed animals. Thirty bison possessed a type B protozoan population, characterized by Epidinium sp., Eudiplodinium maggii, and Eudiplodinium bovis. Thirty-eight bison possessed a mixed A-B population, characterized by Polyplastron sp. coexisting with low numbers of Eudiplodinium maggii or Epidinium sp. or both. Thirteen bison possessed populations lacking any remnant type B ciliate species. At least 29 of the bison possessing Polyplastron sp. were known to have been in contact with cattle, whereas all bison isolated from cattle had type B populations. The reduction of type B populations in bison becomes increasingly likely as bison production expands into areas inhabited by domestic livestock.     

Comparing the impact of a grazing regime with European bison versus one with free-ranging cattle on coastal dune vegetation in the Netherlands

Woody plant encroachment has increased across the globe and threatens biodiversity associated with open habitats. In order to prevent or reduce woody encroachment, conservation managers across Europe introduce large mammalian herbivores. While up to recently, managers were mostly using free-ranging domestic cattle and horses for this, there is an increasing interest in the use of European bison for nature management. However, we lack studies that compare the impact of these different grazers on vegetation. We report results from a unique grazing pilot in the National Park Zuid-Kennemerland, a heterogeneous coastal dune landscape in the Netherlands, where European bison, horses, and cattle were introduced to reverse the encroachment of grass and shrub species. We present results of an 8-year study on the development of woody and grassy vegetation on fixed transects in three different grazing areas within the national park; one area with European bison and horses, one area with cattle and horses, and one area where these large grazers were excluded. In all three areas, rabbit, fallow deer, and roe deer were present. Over time, we observed strong reductions in the vitality of several woody species, such as spindle tree, and this decline was similar across all areas. Grass height and cover also declined and the proportion of herbs increased in all three grazing areas in similar ways. However, the type of herbivore use (debarking, foraging on buds, branches) of several woody species differed significantly among areas. For instance, maple tree was only debarked in the E. bison area, while hawthorn branches were eaten significantly more in the cattle than in the bison area. Due to differences in herbivore densities among areas, it was difficult to draw strong conclusions on how the different herbivore species differed in their impact, but, importantly, we found that grazing regimes with bison can lead to as strong effects on vegetation structure and composition as grazing regimes with cattle. This is an important result since certain conditions, such as legal aspects, may motivate managers to introduce a wild large grazer rather than a domesticated one.     

Microsatellite DNA Variation and the Evolution, Domestication and Phylogeography of Taurine and Zebu Cattle (Bos Taurus and Bos Indicus)

Genetic variation at 20 microsatellite loci was surveyed to determine the evolutionary relationships and molecular biogeography of 20 different cattle populations from Africa, Europe and Asia. Phylogenetic reconstruction and multivariate analysis highlighted a marked distinction between humpless (taurine) and humped (zebu) cattle, providing strong support for a separate origin for domesticated zebu cattle. A molecular clock calculation using bison (Bison sp.) as an outgroup gave an estimated divergence time between the two subspecies of 610,000-850,000 years. Substantial differences in the distribution of alleles at 10 of these loci were observed between zebu and taurine cattle. These markers subsequently proved very useful for investigations of gene flow and admixture in African populations. When these data were considered in conjunction with previous mitochondrial and Y chromosomal studies, a distinctive male-mediated pattern of zebu genetic introgression was revealed. The introgression of zebu-specific alleles in African cattle afforded a high resolution perspective on the hybrid nature of African cattle populations and also suggested that certain West African populations of valuable disease-tolerant taurine cattle are under threat of genetic absorption by migrating zebu herds.     

Detection of mitochondrial DNA from domestic cattle in bison on Santa Catalina Island

In 1924, 14 American bison (Bison bison) were introduced to Santa Catalina Island, California and sporadically supplemented thereafter with additional animals. To reduce the herd and its impact on native vegetation, over 2000 animals have been exported during the past four decades. Today, the herd is estimated to contain around 250 individuals. Genetic analysis was performed on 98 animals removed from the island in 2004. Forty-four samples (45%) had domestic cattle mitochondrial DNA (mtDNA), 12 (12%) had previously reported bison haplotypes and 42 (43%) had a new haplotype differing by one base pair from a previously reported bison haplotype. A complement of five restriction enzymes was found to be useful in identifying bison with domestic cattle mtDNA.     

Intense harvesting of eastern wolves facilitated hybridization with coyotes

Despite ethical arguments against lethal control of wildlife populations, culling is routinely used for the management of predators, invasive or pest species, and infectious diseases. Here, we demonstrate that culling of wildlife can have unforeseen impacts that can be detrimental to future conservation efforts. Specifically, we analyzed genetic data from eastern wolves (Canis lycaon) sampled in Algonquin Provincial Park (APP), Ontario, Canada from 1964 to 2007. Research culls in 1964 and 1965 killed the majority of wolves within a study region of APP, accounting for approximately 36% of the park's wolf population at a time when coyotes were colonizing the region. The culls were followed by a significant decrease in an eastern wolf mitochondrial DNA (mtDNA) haplotype (C1) in the Park's wolf population, as well as an increase in coyote mitochondrial and nuclear DNA. The introgression of nuclear DNA from coyotes, however, appears to have been curtailed by legislation that extended wolf protection outside park boundaries in 2001, although eastern wolf mtDNA haplotype C1 continued to decline and is now rare within the park population. We conclude that the wolf culls transformed the genetic composition of this unique eastern wolf population by facilitating coyote introgression. These results demonstrate that intense localized harvest of a seemingly abundant species can lead to unexpected hybridization events that encumber future conservation efforts. Ultimately, researchers need to contemplate not only the ethics of research methods, but also that future implications may be obscured by gaps in our current scientific understanding.     

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.