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

A Reference Genome Assembly of American Bison, Bison bison bison

Bison are an icon of the American West and an ecologically, commercially, and culturally important species. Despite numbering in the hundreds of thousands today, conservation concerns remain for the species, including the impact on genetic diversity of a severe bottleneck around the turn of the 20th century and genetic introgression from domestic cattle. Genetic diversity and admixture are best evaluated at genome-wide scale, for which a high-quality reference is necessary. Here, we use trio binning of long reads from a bison-Simmental cattle (Bos taurus taurus) male F1 hybrid to sequence and assemble the genome of the American plains bison (Bison bison bison). The male haplotype genome is chromosome-scale, with a total length of 2.65 Gb across 775 scaffolds (839 contigs) and a scaffold N50 of 87.8 Mb. Our bison genome is ~13× more contiguous overall and ~3400× more contiguous at the contig level than the current bison reference genome. The bison genome sequence presented here (ARS-UCSC_bison1.0) will enable new research into the evolutionary history of this iconic megafauna species and provide a new tool for the management of bison populations in federal and commercial herds.

Were domestic camelids present on the prehispanic South American agricultural frontier? An ancient DNA study

The southern boundary of prehispanic farming in South America occurs in central Mendoza Province, Argentina at approximately 34 degrees south latitude. Archaeological evidence of farming includes the recovery of macrobotanical remains of cultigens and isotopic chemistry of human bone. Since the 1990s, archaeologists have also hypothesized that the llama (Lama glama), a domesticated South American camelid, was also herded near the southern boundary of prehispanic farming. The remains of a wild congeneric camelid, the guanaco (Lama guanicoe), however, are common in archaeological sites throughout Mendoza Province. It is difficult to distinguish bones of the domestic llama from wild guanaco in terms of osteological morphology, and therefore, claims that llama were in geographic areas where guanaco were also present based on osteometric analysis alone remain equivocal. A recent study, for example, claimed that twenty-five percent of the camelid remains from the high elevation Andes site of Laguna del Diamante S4 were identified based on osteometric evidence as domestic llama, but guanaco are also a likely candidate since the two species overlap in size. We test the hypothesis that domesticated camelids occurred in prehispanic, southern Mendoza through analysis of ancient DNA. We generated whole mitochondrial genome datasets from 41 samples from southern Mendoza late Holocene archaeological sites, located between 450 and 3400 meters above sea level (masl). All camelid samples from those sites were identified as guanaco; thus, we have no evidence to support the hypothesis that the domestic llama occurred in prehispanic southern Mendoza.

The complete mitochondrial genome sequence analysis of Tibetan argali (Ovis ammon hodgsoni): implications of Tibetan argali and Gansu argali as the same subspecies

The genus Ovis (Bovidae, Artiodactyla) includes six species, i.e. Ovis ammon, Ovis aries, Ovis canadensis, Ovis dalli, Ovis nivicola and Ovis vignei. Based on morphology, geographical location, habitat, etc., the species O. ammon is divided into nine subspecies. The near threatened Tibetan argali is distributed across the Tibetan Plateau and its peripheral mountains, and believed to be one of the O. ammon subspecies (O. a. hodgsoni). However, considering its morphological features and distributions, a question has been proposed by some researchers about the subspecies status of Tibetan argali. In this study, we employed complete mitochondrial DNA (mtDNA) to explore the phylogenetic relationship and population genetic structure of Tibetan argali. The results revealed that the nucleotide composition, gene arrangement and codon usage pattern of the mitochondrial genome of Tibetan argali are similar to those of other caprines. Phylogenetic analyses showed that Tibetan argali was clustered with O. ammon. Interestingly, five Tibetan argali individuals and one of the three Gansu argali (O. a. dalailamae) individuals were clustered in the same branch, which is a sister group to other two Gansu argali individuals. Together with morphological characteristics, our results suggested that Tibetan argali and Gansu argali may belong to the same subspecies (O. a. hodgsoni) of O. ammon, rather than two different subspecies.

Paternally inherited genetic markers reveal new insights into genetic structuring within Macropus fuliginosus and hybridisation with sympatric Macropus giganteus

There are several aspects of biology in which the contribution of males and females is unequal. In these instances the examination of Y chromosome markers may be used to elucidate male-specific attributes. Here, male dispersal patterns and genetic structuring were examined using four Y-microsatellite loci in 186 male western grey kangaroos, Macropus fuliginosus, from throughout the species’ trans-continental distribution. In addition, 52 male grey kangaroos were examined to investigate hybridisation between M. fuliginosus and the eastern grey kangaroo, Macropus giganteus, in their region of sympatry in eastern Australia. Detected Y chromosome diversity was low, resulting from low effective male population size due to skewed sex ratios and a polygynous mating system. As expected, male dispersal was high across the range. However, the Lake Torrens–Flinders Ranges region appears to have significantly restricted male movement between eastern and central/western Australia. There was little evidence to suggest that other barriers (Nullarbor Plain and Swan River Valley) previously identified by nuclear and mitochondrial DNA marker studies restrict male movement. Hence, the admixture events previously identified may be associated with high male dispersal. Within the region of sympatry between M. fuliginosus and M. giganteus in eastern Australia, four M. giganteus individuals were found to possess M. fuliginosus Y-haplotypes. These results confirm the occurrence of hybridisation between male M. fuliginosus and female M. giganteus. Additionally, the introgression of M. fuliginosus Y-haplotypes into M. giganteus populations indicates that at least some male hybrids are fertile, despite evidence to the contrary from captive studies. This study has provided insights into the male contribution to population history, structure and hybridisation in M. fuliginosus, which were not predicted by comparisons between biparentally and maternally inherited markers. This highlights the importance of direct examination of the Y chromosome to provide novel insights into male-mediated processes, especially where the contribution of the sexes may differ.

Molecular detection of hybridization between sympatric kangaroo species in south-eastern Australia

Introgressive hybridization has traditionally been regarded as rare in many vertebrate groups, including mammals. Despite a propensity to hybridize in captivity, introgression has rarely been reported between wild sympatric macropodid marsupials. Here we investigate sympatric populations of western (Macropus fuliginosus) and eastern (Macropus giganteus) grey kangaroos through 12 autosomal microsatellite loci and 626 bp of the hypervariable mitochondrial DNA (mtDNA) control region. M. fuliginosus and M. giganteus within the region of sympatry corresponded, both genetically and morphologically, to their respective species elsewhere in their distributions. Of the 223 grey kangaroos examined, 7.6% displayed evidence of introgression, although no F1 hybrids were detected. In contrast to captive studies, there was no evidence for unidirectional hybridization in sympatric grey kangaroos. However, a higher portion of M. giganteus backcrosses existed within the sample compared with M. fuliginosus. Hybridization in grey kangaroos is reflective of occasional breakdowns in species boundaries, occurring throughout the region and potentially associated with variable conditions and dramatic reductions in densities. Such rare hybridization events allow populations to incorporate novel diversity while still retaining species integrity.

Maternal and paternal genealogy of Eurasian taurine cattle (Bos taurus)

Maternally inherited mitochondrial DNA (mtDNA) has been used extensively to determine origin and diversity of taurine cattle (Bos taurus) but global surveys of paternally inherited Y-chromosome diversity are lacking. Here, we provide mtDNA information on previously uncharacterised Eurasian breeds and present the most comprehensive Y-chromosomal microsatellite data on domestic cattle to date. The mitochondrial haplogroup T3 was the most frequent, whereas T4 was detected only in the Yakutian cattle from Siberia. The mtDNA data indicates that the Ukrainian and Central Asian regions are zones where hybrids between taurine and zebu (B. indicus) cattle have existed. This zebu influence appears to have subsequently spread into southern and southeastern European breeds. The most common Y-chromosomal microsatellite haplotype, termed here as H11, showed an elevated frequency in the Eurasian sample set compared with that detected in Near Eastern and Anatolian breeds. The taurine Y-chromosomal microsatellite haplotypes were found to be structured in a network according to the Y-haplogroups Y1 and Y2. These data do not support the recent hypothesis on the origin of Y1 from the local European hybridization of cattle with male aurochsen. Compared with mtDNA, the intensive culling of breeding males and male-mediated crossbreeding of locally raised native breeds has accelerated loss of Y-chromosomal variation in domestic cattle, and affected the contribution of genetic drift to diversity. In conclusion, to maintain diversity, breeds showing rare Y-haplotypes should be prioritised in the conservation of cattle genetic resources.

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.

An examination of positive selection and changing effective population size in Angus and Holstein cattle populations (Bos taurus) using a high density SNP genotyping platform and the contribution of ancient polymorphism to genomic diversity in Domestic cattle

Background

Identifying recent positive selection signatures in domesticated animals could provide information on genome response to strong directional selection from domestication and artificial selection. With the completion of the cattle genome, private companies are now providing large numbers of polymorphic markers for probing variation in domestic cattle (Bos taurus). We analysed over 7,500 polymorphic single nucleotide polymorphisms (SNP) in beef (Angus) and dairy (Holstein) cattle and outgroup species Bison, Yak and Banteng in an indirect test of inbreeding and positive selection in Domestic cattle.

Results

Outgroup species: Bison, Yak and Banteng, were genotyped with high levels of success (90%) and used to determine ancestral and derived allele states in domestic cattle. Frequency spectrums of the derived alleles in Angus and Holstein were examined using Fay and Wu's H test. Significant divergences from the predicted frequency spectrums expected under neutrality were identified. This appeared to be the result of combined influences of positive selection, inbreeding and ascertainment bias for moderately frequent SNP. Approximately 10% of all polymorphisms identified as segregating in B. taurus were also segregating in Bison, Yak or Banteng; highlighting a large number of polymorphisms that are ancient in origin.

Conclusion

These results suggest that a large effective population size (N_e) of approximately 90,000 or more existed in B. taurus since they shared a common ancestor with Bison, Yak and Banteng ~1–2 million years ago (MYA). More recently N_e decreased sharply probably associated with domestication. This may partially explain the paradox of high levels of polymorphism in Domestic cattle and the relatively small recent N_e in this species. The period of inbreeding caused Fay and Wu's H statistic to depart from its expectation under neutrality mimicking the effect of selection. However, there was also evidence for selection, because high frequency derived alleles tended to cluster near each other on the genome.

Horse domestication and conservation genetics of Przewalski's horse inferred from sex chromosomal and autosomal sequences

Despite their ability to interbreed and produce fertile offspring, there is continued disagreement about the genetic relationship of the domestic horse (Equus caballus) to its endangered wild relative, Przewalski's horse (Equus przewalskii). Analyses have differed as to whether or not Przewalski's horse is placed phylogenetically as a separate sister group to domestic horses. Because Przewalski's horse and domestic horse are so closely related, genetic data can also be used to infer domestication-specific differences between the two. To investigate the genetic relationship of Przewalski's horse to the domestic horse and to address whether evolution of the domestic horse is driven by males or females, five homologous introns (a total of approximately 3 kb) were sequenced on the X and Y chromosomes in two Przewalski's horses and three breeds of domestic horses: Arabian horse, Mongolian domestic horse, and Dartmoor pony. Five autosomal introns (a total of approximately 6 kb) were sequenced for these horses as well. The sequences of sex chromosomal and autosomal introns were used to determine nucleotide diversity and the forces driving evolution in these species. As a result, X chromosomal and autosomal data do not place Przewalski's horses in a separate clade within phylogenetic trees for horses, suggesting a close relationship between domestic and Przewalski's horses. It was also found that there was a lack of nucleotide diversity on the Y chromosome and higher nucleotide diversity than expected on the X chromosome in domestic horses as compared with the Y chromosome and autosomes. This supports the hypothesis that very few male horses along with numerous female horses founded the various domestic horse breeds. Patterns of nucleotide diversity among different types of chromosomes were distinct for Przewalski's in contrast to domestic horses, supporting unique evolutionary histories of the two species.

Reduced genetic structure of north Ethiopian cattle revealed by Y-chromosome analysis

Ethiopia is considered to be a putative migratory corridor for both Near-East Bos taurine and Arabian and Indian B. indicus cattle into East Africa. African pastoralism, which is associated with adaptation to specific habitats and farming systems, has contributed to the composite constitution of Ethiopian cattle. We analyse, for the first time, five Y-chromosome microsatellite markers from seven north Ethiopian cattle populations, using a European Holstein-Friesian population as a reference, to assess the paternal gene pool and to explore the mechanisms behind the genetic structure. Our results reveal that the indicine alleles predominate in the present populations, with only one animal in the Arado carrying the taurine alleles. The north Ethiopian cattle populations with one exception (Abergelle) are characterized by a general low Y-chromosome haplotype diversity, as well as by a reduced interpopulation variance (Phi(ST)=4.0%), which can be a result of strong male-mediated selective sweeps. Population structure revealed by multidimensional-scaling analysis differentiates two populations (Arado and Abergelle) from the rest. Analysis of molecular variance does not lend support to the traditional classification for the populations, which is mainly based on physical characteristics. A network analysis indicates two closely related founding haplotypes accounting for a large proportion (50.0% in Abergelle and 85.0-94.7% in others) of north Ethiopian cattle Y-chromosomes. Our findings point to a common, but limited, paternal origin of the north Ethiopian cattle populations and strong male-mediated gene flow among them. The findings also provide insight into the historical immigration of cattle into East Africa.

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.

Fecal androgens of bison bulls during the rut

The influence of sex hormones is a key proximate factor underlying male reproductive behavior in mammals. Effective conservation policies for the remaining purebred plains bison (Bison bison bison) herds require knowledge of the physiology underlying bison reproductive biology. We used fecal steroid analysis to characterize androgen levels in adult bison bulls before, during, and after the rut, and to examine androgen levels of bulls differing in reproductive status, age, and mating success. Fieldwork was carried out at the Fort Niobrara National Wildlife Refuge in north-central Nebraska. All adult bison in the herd were individually known by unique brands. Fecal samples were collected during 2003 from bulls during pre-rut (June), rut (July-August), and post-rut (September), and behavioral observations focused on reproductive status and mating success during the rut. Matched sample data indicated that androgen levels (ng/g feces) of bulls peaked during the rut, doubling from pre-rut to rut and then declining by 75% during post-rut. Dominant bulls that tended (guarded) cows maintained higher androgen levels than bulls that were not tending. There was a positive correlation between bull age (associated with mating success) and androgens, with higher androgen levels in prime-aged bulls compared with younger bulls. Nonetheless, there was no correlation between mating success (measured by number of copulations observed) and androgen level. This suggests that while androgens may provide the proximate motivation to compete for matings, other factors determine the mating success of bison bulls.

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.

Genetic variation in Przewalski’s horses, with special focus on the last wild caught mare, 231 Orlitza III

In our continuing efforts to document genetic diversity in Przewalski's horses and relatedness with domestic horses, we report genetic variation at 22 loci of blood group and protein polymorphisms and 29 loci of DNA (microsatellite) polymorphisms. The loci have been assigned by linkage or synteny mapping to 20 autosomes and the X chromosome of the domestic horse (plus four loci unassigned to a chromosome). With cumulative data from tests of 568 Przewalski's horses using blood, hair or tooth samples, no species-defining markers were identified, however a few markers were present in the wild species but not in domestic horses. Inheritance patterns and linkage relationships reported in domestic horses appeared to be conserved in Przewalski's horses. A derived type for the last wild caught mare 231 Orlitza III provided evidence for markers apparently not found in (or not currently available by descent from) the other species founders that were captured at the end of the nineteenth century. This information has been critical to the development of parentage analyses in the studbook population of Przewalski's horses at Askania Nova, at one time the largest herd of captive animals and the source of stock for reintroduction efforts. Some horses in the study showed genetic incompatibilities with their sire or dam, contradicting published studbook information. In many cases alternative parentage could be assigned from living animals. To assist in identification of correct parentage, DNA marker types for deceased horses were established from archived materials (teeth) or derived from offspring. Genetic markers were present in pedigreed animals whose origin could not be accounted for from founders. Genetic distance analysis of erythrocyte protein, electrophoretic and microsatellite markers in Przewlaski's horses and ten breeds of domestic horse place the Przewalski's horse as an outgroup to domestic horses, introgression events from domestic horses not withstanding.

Identification of bovine and novel interferon-? alleles in the American plains bison (bison bison) by analysis of hybrid cattle x bison blastocysts

The objective of this study was to generate bison x cattle hybrid embryos by in vitro fertilization, to assess their developmental potential, to determine the pattern of secretion of the embryonic signaling molecule interferon-tau (IFN-tau), and to identify novel IFN-tau mRNA polymorphism in the American plains bison. A total of 600 bovine oocytes were inseminated with frozen-thawed bison semen. Of these, 40.7% cleaved and 14.8% proceeded to the blastocyst stage. Individual blastocysts were cultured on a basement membrane (Matrigel) and their ability to attach and form outgrowths was monitored. A total of 36 blastocysts were cultured of which 22 formed outgrowths. During individual culture, medium samples were collected and their IFN-tau concentration was measured. On day 6 after onset of individual culture, attached outgrowths produced significantly more IFN-tau than unattached viable or degenerate blastocysts. At this time, female conceptuses also produced significantly more IFN-tau than their male cohorts. However, by day 12 this difference had disappeared. Total mRNA was extracted from three individual outgrowths and analyzed by RT-PCR. Subsequent sequencing of 28 clones showed several known bovine IFN-tau sequences as well as two novel sequences termed bisIFN-tau1 and 2. To determine the origin of these, DNA was extracted from bison semen and analyzed by PCR. One bovine IFN-tau sequence (bovIFN-tau1d) as well as bisIFN-tau2 and a third novel sequence bisIFN-tau3 were detected. This study demonstrates the feasibility of using hybrid embryos for the analysis of developmentally regulated gene expression in species where embryos may not be available.