CONSERVATION GENETIC ANALYSIS OF THE TEXAS STATE BISON HERD

When Theory Meets Practice: Balancing Genetic Diversity and Behaviour When Choosing Founders for a Recently Reintroduced Bison (Bison bison) Herd in Banff National Park, Canada

The tragic decline of plains bison (Bison bison bison) in North America is well known. By the turn of the 20th century, plains bison, once numbering in the tens of millions, were extirpated from the wild in Canada. Banff National Park, AB, Canada, is within the original range of bison in North America, and still contains suitable habitat for the species. In 2017, we used measures of genetic diversity, as well as age and reproductive status, to select 16 plains bison for reintroduction to an unoccupied corner of their former historic range in Banff National Park. However, five of the preferred founders from this genetically diverse group needed to be substituted for less ideal candidates when aggressive behavior was observed, which could have compromised animal welfare during translocation. The level of genetic diversity in the founding population was therefore lower than if all proposed founders had been used, although the difference was not substantial. As expected, adults in this small, reintroduced herd had lower observed heterozygosity and number of alleles than the larger source herd, but expected heterozygosity was higher, likely a result of preferentially choosing individuals with rare alleles. Another substantial contribution to genetic diversity was in unborn fetuses: sired by males that were not translocated, the calves born a few months later provided a boost to the genetic diversity of this small founder herd. Where opportunities allow, the inclusion of pregnant females can significantly increase the genetic diversity of small founder populations. When compared to other founded herds of bison, the techniques described here have ensured that the newly established Banff NP herd is one of the most diverse bison herds in Canada, despite its relatively small number of founders.

Opportunities and Limitations for Reproductive Science in Species Conservation

Reproductive science in the context of conservation biology is often understood solely in terms of breeding threatened species. Although technologies developed primarily for agriculture or biomedicine have a potentially important role in species conservation, their effectiveness is limited if we regard the main objective of animal conservation as helping to support populations rather than to breed a small number of individuals. The global threats facing wild species include the consequences of climate change, population growth, urbanization, atmospheric and water pollution, and the release of chemicals into the environment, to cite but a few. Reproductive sciences provide important and often unexpected windows into many of these consequences, and our aim here is both to demonstrate the breadth of reproductive science and the importance of basic knowledge and to suggest where some of the insights might be useful in mitigating the problems. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Development of SNP-Based Genomic Tools for the Canadian Bison Industry: Parentage Verification and Subspecies Composition

Genomic technologies have been increasingly applied in livestock production due to their utility in production management and animal genetic improvement. The current project aimed to develop genomic resources for the Canadian bison industry, specifically a parentage verification tool and a subspecies composition tool. Both products stand to help with building and maintaining purebred and crossbred bison populations, and in turn bison conservation and production. The development of this genomic toolkit proceeded in two stages. In the single-nucleotide polymorphism (SNP) discovery and selection stage, raw sequence information from 41 bison samples was analyzed, and approximately 52.5 million candidate biallelic SNPs were discovered from 21 samples with high sequence quality. A set of 19,954 SNPs (2,928 for parentage verification and 17,026 for subspecies composition) were then selected for inclusion on an Axiom myDesign custom array. In the refinement and validation stage, 480 bison were genotyped using the custom SNP panel, and the resulting genotypes were analyzed to further filter SNPs and assess tool performance. In various tests using real and simulated genotypes, the two genomic tools showed excellent performance for their respective tasks. Final SNP sets consisting of 191 SNPs for parentage and 17,018 SNPs for subspecies composition are described. As the first SNP-based genomic toolkit designed for the Canadian bison industry, our results may provide a new opportunity in improving the competitiveness and profitability of the industry in a sustainable manner.

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.

Genetic variation and mating success in managed american plains bison

The American plains bison (Bison bison) was pushed to the brink of extinction in the late 1800s but has since rebounded. Less than 5% of animals currently exist in conservation herds that are critical for maintaining genetic variability. Here, we use 25 microsatellite loci to assess genetic diversity and patterns of mating success over a 3-year period in a managed conservation herd at Konza Prairie Biological Station, Kansas (total number of individuals genotyped = 587). Heterozygosity was comparable to and allelic diversity higher than that in 11 other wild and managed herds for which similar estimates are available. Parentage analyses revealed that males within the oldest age classes (5-7 years) sired >90% of calves over the study period, consistent with a polygynous breeding system. Asymmetries in siring success also were observed within age classes, with the same males enjoying high siring success over multiple seasons. Empirical results of paternity will facilitate future modeling and empirical efforts to determine how demographic factors, population size, and variation in siring success interact to determine the retention (or loss) of genetic diversity in natural and managed herds, thus allowing informed recommendations for management practices and conservation efforts of this symbolic North American species.

Polymorphism and haplotype structure in River Buffalo (Bubalus bubalis) toll-like receptor 5 (TLR5) coding sequence

Most of the 160 million river buffalo in the world are in Asia where they are used extensively, both as a food source and for draught power. Only recently have investigations begun exploring the buffalo genome for variation that might influence health and productivity of these economically important animals. This paper describes the sequence variability of the toll-like receptor 5 (TLR5) gene, which recognizes bacterial flagellin and is a key player in the immune system. TLR5 is comprised of a single exon that is 2577 bp and codes 858 amino acids. We examined single-nucleotide polymorphisms (SNPs) located within the coding region. Overall, 17 SNPs were discovered, seven of which are non-synonymous. Our study population yielded four different haplotypes. We examined predicted protein domain structure and found that river buffalo, swamp buffalo, and African Forest buffalo shared the same protein domain structure and are more similar to each other than they are to cattle and American bison, which are similar to each other. PolyPhen 2 analysis revealed one amino acid substitution in the river buffalo population with potential functional significance.

Genetic population substructure in bison at Yellowstone National Park

The Yellowstone National Park bison herd is 1 of only 2 populations known to have continually persisted on their current landscape since pre-Columbian times. Over the last century, the census size of this herd has fluctuated from around 100 individuals to over 3000 animals. Previous studies involving radiotelemetry, tooth wear, and parturition timing provide evidence of at least 2 distinct groups of bison within Yellowstone National Park. To better understand the biology of Yellowstone bison, we investigated the potential for limited gene flow across this population using multilocus Bayesian clustering analysis. Two genetically distinct and clearly defined subpopulations were identified based on both genotypic diversity and allelic distributions. Genetic cluster assignments were highly correlated with sampling locations for a subgroup of live capture individuals. Furthermore, a comparison of the cluster assignments to the 2 principle winter cull sites revealed critical differences in migration patterns across years. The 2 Yellowstone subpopulations display levels of differentiation that are only slightly less than that between populations which have been geographically and reproductively isolated for over 40 years. The identification of cryptic population subdivision and genetic differentiation of this magnitude highlights the importance of this biological phenomenon in the management of wildlife species.

Conservation genetics and North American bison (Bison bison)

The many millions of North American bison in the mid-19th century were reduced to near extinction by the middle 1880s. Plains bison, the subspecies found in the United States, were saved from extinction primarily by 5 private ranchers and the survival of a small herd in what is now Yellowstone National Park. This bottleneck resulted in the present-day plains bison population being descended from less than 100 founders. In addition, many conservation herds have cattle ancestry because of hybridization promoted by these ranchers in the late 1800s and early 1900s. Today, although there are around 500,000 plains bison in North America, only 4% (20,000) are in conservation herds. Only 1 conservation herd with no known ancestry from cattle has an effective population size of more than 1000. Here I review and evaluate this situation and provide recommendations for the reduction of cattle ancestry, avoidance of inbreeding depression, and maintenance of genetic variation in the conservation herds of bison.

Patterns of genetic variation in US federal bison herds

Like many wide-ranging mammals, American bison (Bison bison) have experienced significant range contraction over the past two centuries and are maintained in artificially isolated populations. A basic understanding of the distribution of genetic variation among populations is necessary to facilitate long-term germplasm preservation and species conservation. The 11 herds maintained within the US federal system are a critically important source of germplasm for bison conservation, as they include many of the oldest herds in the USA and have served as a primary resource for the establishment of private and public herds worldwide. In this study, we used a panel of 51 nuclear markers to investigate patterns of neutral genetic variation among these herds. Most of these herds have maintained remarkably high levels of variation despite the severe bottleneck suffered in the late 1800s. However, differences were noted in the patterns of variation and levels of differentiation among herds, which were compared with historical records of establishment, supplementation, herd size, and culling practices. Although some lineages have been replicated across multiple herds within the US federal system, other lineages with high levels of genetic variation exist in isolated herds and should be considered targets for the establishment of satellite herds. From this and other studies, it is clear that the genetic variation represented in the US federal system is unevenly distributed among National Park Service and Fish and Wildlife Service herds, and that these resources must be carefully managed to ensure long-term species conservation.

Patterns of male reproductive success in a highly promiscuous whale species: the endangered North Atlantic right whale

Parentage analyses of baleen whales are rare, and although mating systems have been hypothesized for some species, little data on realized male reproductive success are available and the patterns of male reproductive success have remained elusive for most species. Here we combine over 20 years of photo-identification data with high-resolution genetic data for the majority of individual North Atlantic right whales to assess paternity in this endangered species. There was significant skew in male reproductive success compared to what would be expected if mating was random (P < 0.001). The difference was due to an excess of males assigned zero paternities, a deficiency of males assigned one paternity, and an excess of males assigned as fathers for multiple calves. The variance in male reproductive success was high relative to other aquatically mating marine mammals, but was low relative to mammals where the mating system is based on resource- and/or mate-defence polygyny. These results are consistent with previous data suggesting that the right whale mating system represents one of the most intense examples of sperm competition in mammals, but that sperm competition on its own does not allow for the same degree of polygyny as systems where males can control access to resources and/or mates. The age distribution of assigned fathers was significantly biased towards older males (P < 0.05), with males not obtaining their first paternity until approximately 15 years of age, which is almost twice the average age of first fertilization in females (8 years), suggesting that mate competition is preventing younger males from reproducing. The uneven distribution of paternities results in a lower effective population size in this species that already has one of the lowest reported levels of genetic diversity, which may further inhibit reproductive success through mate incompatibility of genetically similar individuals.

Marker genotypes and population admixture and their association with body weight, height and relative body mass in United States federal bison herds

Elucidating genetic influences on bison growth and body composition is of interest, not only because bison are important for historical, cultural, and agricultural reasons, but also because their unusual population history makes them valuable models for finding influential loci in both domestic cattle and humans. We tested for trait loci associated with body weight, height, and bison mass index (BMI) while controlling for estimated ancestry to reduce potential confounding effects due to population admixture in 1316 bison sampled from four U.S. herds. We used 60 microsatellite markers to model each phenotype as a function of herd, sex, age, marker genotypes, and individual ancestry estimates. Statistical significance for genotype and its interaction with ancestry was evaluated using the adaptive false discovery rate. Of the four herds, two appeared to be admixed and two were nonadmixed. Although none of the main effects of the loci were significant, estimated ancestry and its interaction with marker loci were significantly associated with the phenotypes, illustrating the importance of including ancestry in the models and the dependence of genotype-phenotype associations on background ancestry. Individual loci contributed approximately 2.0% of variation in weight, height, and BMI, which confirms the utility and potential importance of adjusting for population stratification.

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.