Search

Background: Genetic testing for pedigree accuracy is critical for managing genetic diversity in North American (NA) yak ( Bos grunniens), a population expanded mostly from imported zoological park specimens. DNA testing also enhances species conservation by identifying recent B. taurus F1 hybrid ancestors (within three generations). Biallelic single nucleotide polymorphisms (SNPs) can accomplish either task, but increases the marker count and costs necessary to achieve both. Our aim was to identify novel, multifunctional, triallelic yak SNPs (tySNPs), with each having two alleles for yak parentage testing, and a third allele for identifying recent cattle introgression. Methods: Genome sequences were aligned to the cattle UMD3.1 assembly and SNPs were screened for 1) heterozygosity in a NA and a Chinese yak, 2) a third allele at high frequency in cattle, and 3) flanking sequences conserved in both species. Subsequently, tySNPs were filtered for unique alignment to the haplotype-resolved F1 yak assembly. Allele frequencies were estimated in a subset of 87 tySNPs by genotyping 170 NA yak. Results: We identified 610 autosomal tySNPs, distributed in 441 clusters with 5 Mb average genome spacing. The average NA yak minor allele frequency was high (0.296), while average introgressed cattle alleles were low (0.004). In simulations with tySNPs, 28 were sufficient for globally-unique animal identification (P I=5.81x10 -12), 87 were able to exclude 19 random bulls from parentage at the 99% level without using the dam’s genotype (P E=5.3x10 -4), and 87 were able to detect F1 hybridization events after three generations of yak backcrosses (1/16th B. taurus germplasm). Conclusions: Identifying animals, determining parentage and detecting recent hybridization events was efficient with as few as 87 tySNPs. A similar triallelic approach could be used with other bottlenecked Bos species that hybridize with cattle, such as NA plains bison ( B. bison).

Background: Genetic testing for pedigree accuracy is critical for managing genetic diversity in North American (NA) yak ( Bos grunniens), a population expanded mostly from imported zoological park specimens. DNA testing also enhances species conservation by identifying recent B. taurus F1 hybrid ancestors (within three generations). Biallelic single nucleotide polymorphisms (SNPs) can accomplish either task, but increases the marker count and costs necessary to achieve both. Our aim was to identify novel, multifunctional, triallelic yak SNPs (tySNPs), with each having two alleles for yak parentage testing, and a third allele for identifying recent cattle introgression. Methods: Genome sequences were aligned to the cattle UMD3.1 assembly and SNPs were screened for 1) heterozygosity in a NA and a Chinese yak, 2) a third allele at high frequency in cattle, and 3) flanking sequences conserved in both species. Subsequently, tySNPs were filtered for unique alignment to the haplotype-resolved F1 yak assembly. Allele frequencies were estimated in a subset of 87 tySNPs by genotyping 170 NA yak. Results: We identified 610 autosomal tySNPs, distributed in 441 clusters with 5 Mb average genome spacing. The average NA yak minor allele frequency was high (0.296), while average introgressed cattle alleles were low (0.004). In simulations with tySNPs, 28 were sufficient for globally-unique animal identification (P I=5.81x10 -12), 87 were able to exclude 19 random bulls from parentage at the 99% level without using the dam’s genotype (P E=5.3x10 -4), and 87 were able to detect F1 hybridization events after three generations of yak backcrosses (1/16th B. taurus germplasm). Conclusions: Identifying animals, determining parentage and detecting recent hybridization events was efficient with as few as 87 tySNPs. A similar triallelic approach could be used with other bottlenecked Bos species that hybridize with cattle, such as NA plains bison ( B. bison).

Recently, deletions have been identified and published as causal for Familial Adenomatous Polyposis in the 1B promoter region of the APC gene. Those deletions were measured using multiplex ligation-dependent probe amplification. Here, we present and characterize an ~11kb deletion identified by whole genome shotgun sequencing. The deletion occurred in a patient diagnosed with Familial Adenomatous Polyposis, and was located on chr5, between bases 112,034,824 and 112,045,845, fully encompassing the 1B promoter region of the APC gene. Results are presented here that include the sequence evidence supporting the presence of the deletion as well as base level characterization of the deletion site. These results demonstrate the capacity of whole genome sequencing for the detection of large structural variants in single individuals.

Genomics research in mammals has produced reference genome sequences that are essential for identifying variation associated with disease. High quality reference genome sequences are now available for humans, model species, and economically important agricultural animals. Comparisons between these species have provided unique insights into mammalian gene function. However, the number of species with reference genomes is small compared to those needed for studying molecular evolutionary relationships in the tree of life. For example, among the even-toed ungulates there are approximately 300 species whose phylogenetic relationships have been calculated in the 10k trees project. Only six of these have reference genomes: cattle, swine, sheep, goat, water buffalo, and bison. Although reference sequences will eventually be developed for additional hoof stock, the resources in terms of time, money, infrastructure and expertise required to develop a quality reference genome may be unattainable for most species for at least another decade. In this work we mapped 35 Gb of next generation sequence data of a Katahdin sheep to its own species’ reference genome ( Ovis aries Oar3.1) and to that of a species that diverged 15 to 30 million years ago ( Bos taurus UMD3.1). In total, 56% of reads covered 76% of UMD3.1 to an average depth of 6.8 reads per site, 83 million variants were identified, of which 78 million were homozygous and likely represent interspecies nucleotide differences. Excluding repeat regions and sex chromosomes, nearly 3.7 million heterozygous sites were identified in this animal vs. bovine UMD3.1, representing polymorphisms occurring in sheep. Of these, 41% could be readily mapped to orthologous positions in ovine Oar3.1 with 80% corroborated as heterozygous. These variant sites, identified via interspecies mapping could be used for comparative genomics, disease association studies, and ultimately to understand mammalian gene function.

Genomics research in mammals has produced reference genome sequences that are essential for identifying variation associated with disease. High quality reference genome sequences are now available for humans, model species, and economically important agricultural animals. Comparisons between these species have provided unique insights into mammalian gene function. However, the number of species with reference genomes is small compared to those needed for studying molecular evolutionary relationships in the tree of life. For example, among the even-toed ungulates there are approximately 300 species whose phylogenetic relationships have been calculated in the 10k trees project. Only six of these have reference genomes: cattle, swine, sheep, goat, water buffalo, and bison. Although reference sequences will eventually be developed for additional hoof stock, the resources in terms of time, money, infrastructure and expertise required to develop a quality reference genome may be unattainable for most species for at least another decade. In this work we mapped 35 Gb of next generation sequence data of a Katahdin sheep to its own species’ reference genome ( Ovis aries Oar3.1) and to that of a species that diverged 15 to 30 million years ago ( Bos taurus UMD3.1). In total, 56% of reads covered 76% of UMD3.1 to an average depth of 6.8 reads per site, 83 million variants were identified, of which 78 million were homozygous and likely represent interspecies nucleotide differences. Excluding genome repeat regions and sex chromosomes, approximately 3.7 million heterozygous sites were identified in this animal vs. bovine UMD3.1, representing polymorphisms occurring in sheep. Of these, 41% could be readily mapped to orthologous positions in ovine Oar3.1 with 80% corroborated as heterozygous. These variant sites, identified via interspecies mapping could be used for comparative genomics, disease association studies, and ultimately to understand mammalian gene function.