Michael P. Heaton

Development and Characterization of a High Density SNP Genotyping Assay for Cattle

The success of genome-wide association (GWA) studies for the detection of sequence variation affecting complex traits in human has spurred interest in the use of large-scale high-density single nucleotide polymorphism (SNP) genotyping for the identification of quantitative trait loci (QTL) and for marker-assisted selection in model and agricultural species. A cost-effective and efficient approach for the development of a custom genotyping assay interrogating 54,001 SNP loci to support GWA applications in cattle is described. A novel algorithm for achieving a compressed inter-marker interval distribution proved remarkably successful, with median interval of 37 kb and maximum predicted gap of <350 kb. The assay was tested on a panel of 576 animals from 21 cattle breeds and six outgroup species and revealed that from 39,765 to 46,492 SNP are polymorphic within individual breeds (average minor allele frequency (MAF) ranging from 0.24 to 0.27). The assay also identified 79 putative copy number variants in cattle. Utility for GWA was demonstrated by localizing known variation for coat color and the presence/absence of horns to their correct genomic locations. The combination of SNP selection and the novel spacing algorithm allows an efficient approach for the development of high-density genotyping platforms in species having full or even moderate quality draft sequence. Aspects of the approach can be exploited in species which lack an available genome sequence. The BovineSNP50 assay described here is commercially available from Illumina and provides a robust platform for mapping disease genes and QTL in cattle.

A Genome Wide Survey of SNP Variation Reveals the Genetic Structure of Sheep Breeds

The genetic structure of sheep reflects their domestication and subsequent formation into discrete breeds. Understanding genetic structure is essential for achieving genetic improvement through genome-wide association studies, genomic selection and the dissection of quantitative traits. After identifying the first genome-wide set of SNP for sheep, we report on levels of genetic variability both within and between a diverse sample of ovine populations. Then, using cluster analysis and the partitioning of genetic variation, we demonstrate sheep are characterised by weak phylogeographic structure, overlapping genetic similarity and generally low differentiation which is consistent with their short evolutionary history. The degree of population substructure was, however, sufficient to cluster individuals based on geographic origin and known breed history. Specifically, African and Asian populations clustered separately from breeds of European origin sampled from Australia, New Zealand, Europe and North America. Furthermore, we demonstrate the presence of stratification within some, but not all, ovine breeds. The results emphasize that careful documentation of genetic structure will be an essential prerequisite when mapping the genetic basis of complex traits. Furthermore, the identification of a subset of SNP able to assign individuals into broad groupings demonstrates even a small panel of markers may be suitable for applications such as traceability.

Selection and use of SNP markers for animal identification and paternity analysis in U.S. beef cattle

Abstract. DNA marker technology represents a promising means for determining the genetic identity and kinship of an animal. Compared with other types of DNA markers, single nucleotide polymorphisms (SNPs) are attractive because they are abundant, genetically stable, and amenable to high-throughput automated analysis. In cattle, the challenge has been to identify a minimal set of SNPs with sufficient power for use in a variety of popular breeds and crossbred populations. This report describes a set of 32 highly informative SNP markers distributed among 18 autosomes and both sex chromosomes. Informativity of these SNPs in U.S. beef cattle populations was estimated from the distribution of allele and genotype frequencies in two panels: one consisting of 96 purebred sires representing 17 popular breeds, and another with 154 purebred American Angus from six herds in four Midwestern states. Based on frequency data from these panels, the estimated probability that two randomly selected, unrelated individuals will possess identical genotypes for all 32 loci was 2.0 × 10−13 for multi-breed composite populations and 1.9 × 10−10 for purebred Angus populations. The probability that a randomly chosen candidate sire will be excluded from paternity was estimated to be 99.9% and 99.4% for the same respective populations. The DNA immediately surrounding the 32 target SNPs was sequenced in the 96 sires of the multi-breed panel and found to contain an additional 183 polymorphic sites. Knowledge of these additional sites, together with the 32 target SNPs, allows the design of robust, accurate genotype assays on a variety of high-throughput SNP genotyping platforms.

Resolving the evolution of extant and extinct ruminants with high-throughput phylogenomics

The Pecorans (higher ruminants) are believed to have rapidly speciated in the Mid-Eocene, resulting in five distinct extant families: Antilocapridae, Giraffidae, Moschidae, Cervidae, and Bovidae. Due to the rapid radiation, the Pecoran phylogeny has proven difficult to resolve, and 11 of the 15 possible rooted phylogenies describing ancestral relationships among the Antilocapridae, Giraffidae, Cervidae, and Bovidae have each been argued as representations of the true phylogeny. Here we demonstrate that a genome-wide single nucleotide polymorphism (SNP) genotyping platform designed for one species can be used to genotype ancient DNA from an extinct species and DNA from species diverged up to 29 million years ago and that the produced genotypes can be used to resolve the phylogeny for this rapidly radiated infraorder. We used a high-throughput assay with 54,693 SNP loci developed for Bos taurus taurus to rapidly genotype 678 individuals representing 61 Pecoran species. We produced a highly resolved phylogeny for this diverse group based upon 40,843 genome-wide SNP, which is five times as many informative characters as have previously been analyzed. We also establish a method to amplify and screen genomic information from extinct species, and place Bison priscus within the Bovidae. The quality of genotype calls and the placement of samples within a well-supported phylogeny may provide an important test for validating the fidelity and integrity of ancient samples. Finally, we constructed a phylogenomic network to accurately describe the relationships between 48 cattle breeds and facilitate inferences concerning the history of domestication and breed formation.

Worldwide Patterns of Ancestry, Divergence, and Admixture in Domesticated Cattle

The domestication and development of cattle has considerably impacted human societies, but the histories of cattle breeds and populations have been poorly understood especially for African, Asian, and American breeds. Using genotypes from 43,043 autosomal single nucleotide polymorphism markers scored in 1,543 animals, we evaluate the population structure of 134 domesticated bovid breeds. Regardless of the analytical method or sample subset, the three major groups of Asian indicine, Eurasian taurine, and African taurine were consistently observed. Patterns of geographic dispersal resulting from co-migration with humans and exportation are recognizable in phylogenetic networks. All analytical methods reveal patterns of hybridization which occurred after divergence. Using 19 breeds, we map the cline of indicine introgression into Africa. We infer that African taurine possess a large portion of wild African auroch ancestry, causing their divergence from Eurasian taurine. We detect exportation patterns in Asia and identify a cline of Eurasian taurine/indicine hybridization in Asia. We also identify the influence of species other than Bos taurus taurus and B. t. indicus in the formation of Asian breeds. We detect the pronounced influence of Shorthorn cattle in the formation of European breeds. Iberian and Italian cattle possess introgression from African taurine. American Criollo cattle originate from Iberia, and not directly from Africa with African ancestry inherited via Iberian ancestors. Indicine introgression into American cattle occurred in the Americas, and not Europe. We argue that cattle migration, movement and trading followed by admixture have been important forces in shaping modern bovine genomic variation.

Prion gene sequence variation within diverse groups of U.S. sheep, beef cattle, and deer

Prions are proteins that play a central role in transmissible spongiform encephalopathies in a variety of mammals. Among the most notable prion disorders in ungulates are scrapie in sheep, bovine spongiform encephalopathy in cattle, and chronic wasting disease in deer. Single nucleotide polymorphisms in the sheep prion gene (PRNP) have been correlated with susceptibility to natural scrapie in some populations. Similar correlations have not been reported in cattle or deer; however, characterization of PRNP nucleotide diversity in those species is incomplete. This report describes nucleotide sequence variation and frequency estimates for the PRNP locus within diverse groups of U.S. sheep, U.S. beef cattle, and free-ranging deer (Odocoileusvirginianus and O. hemionus from Wyoming). DNA segments corresponding to the complete prion coding sequence and a 596-bp portion of the PRNP promoter region were amplified and sequenced from DNA panels with 90 sheep, 96 cattle, and 94 deer. Each panel was designed to contain the most diverse germplasm available from their respective populations to facilitate polymorphism detection. Sequence comparisons identified a total of 86 polymorphisms. Previously unreported polymorphisms were identified in sheep (9), cattle (13), and deer (32). The number of individuals sampled within each population was sufficient to detect more than 95% of all alleles present at a frequency greater than 0.02. The estimation of PRNP allele and genotype frequencies within these diverse groups of sheep, cattle, and deer provides a framework for designing accurate genotype assays for use in genetic epidemiology, allele management, and disease control.

Single nucleotide polymorphism (SNP) discovery and linkage mapping of bovine cytokine genes

Abstract. Polymorphic markers at bovine gene loci facilitate the integration of cattle genetic maps with those of humans and mice. To this end, 31 single nucleotide polymorphism (SNP) markers were developed for seven bovine chemokine genes. Loci were amplified from bovine genomic DNA by the polymerase chain reaction, and candidate amplicons were sequenced to determine their identity. Amplified loci from 24 founding parents and select progeny from a beef cattle reference population were sequenced and analyzed for SNPs. SNP haplotype alleles were determined by examining segregation patterns and used to establish the locus position on the bovine linkage map. Loci for growth-related proteins (GRO3, GRO1, and GROX) were clustered with the related CXC chemokine genes, interleukin (IL) 8, and epithelial cell inflammatory protein 1, at 84 cM from the centromeric end of the bovine chromosome (BTA) 6 linkage group. Bovine loci for a cluster of IL8 receptors, a stromal cell-derived factor 1, interferon-γ, and tumor necrosis factor-α were mapped at 90, 55, 59, and 34 cM, respectively, from the centromeric ends of the BTA 2, 28, 5, and 23 linkage groups. The positions of these bovine loci were compared with those of orthologous loci on the human map to refine the boundaries of conserved synteny. These seven loci provide examples of SNP development in which the efficiency was largely dependent on the availability of bovine genomic or cDNA sequence. The polymorphic nature of these SNP haplotype markers suggests that they will be useful for mapping complex traits in cattle, such as resistance to infectious disease.

Mobilization of vancomycin resistance by transposon-mediated fusion of a VanA plasmid with an Enterococcus faecium sex pheromone-response plasmid.

A striking feature of recent outbreaks of vancomycin-resistant (VmR) enterococci is the apparent horizontal dissemination of resistance determinants. The plasmids pHKK702 and pHKK703 from Enterococcus faecium clinical isolate R7 have been implicated in the conjugal transfer of VmR. pHKK702 is a 41-kb plasmid that contains an element indistinguishable from the glycopeptide-resistance transposon Tn1546. pHKK703 is an approx. 55-kb putative sex pheromone-response plasmid that is required for conjugative mobilization of pHKK702. During experiments in which strain R7 was used as a donor, a highly conjugative VmR transconjugant was isolated that formed constitutive cellular aggregates. Restriction analyses and DNA hybridizations revealed that the transconjugant harbored a single plasmid of approx. 92 kb and this plasmid (pHKK701) was composed of DNA from both pHKK702 and pHKK703. Results from DNA sequence analyses showed that a 39-kb composite transposon (Tn5506) from pHKK702 had inserted into pHKK703. The left end of Tn5506 contained a single insertion sequence (IS) element, IS1216V2, whereas the right end was composed of a tandem IS structure consisting of the novel 1065-bp IS1252 nested within an IS1216V1 element. Transposition of Tn5506 from pHKK702 to pHKK703 created an 8-bp target sequence duplication at the site of insertion and interrupted an ORF (ORFX) that was 91% identical to that of prgX, a gene proposed to negatively regulate sex pheromone response of the E.faecalis plasmid, pCF10. We propose that the interruption of ORFX by Tn5506 led to the constitutive cellular aggregation phenotype and thereby enhanced the efficiency with which VmR was transferred. Similar IS1216V-mediated transposition events may contribute to the horizontal spread of glycopeptide resistance among enterococci in nature.

Association of Escherichia coli O157:H7 tir polymorphisms with human infection

BackgroundEmerging molecular, animal model and epidemiologic evidence suggests that Shiga-toxigenic Escherichia coli O157:H7 (STEC O157) isolates vary in their capacity to cause human infection and disease. The translocated intimin receptor (tir) and intimin (eae) are virulence factors and bacterial receptor-ligand proteins responsible for tight STEC O157 adherence to intestinal epithelial cells. They represent logical genomic targets to investigate the role of sequence variation in STEC O157 pathogenesis and molecular epidemiology. The purposes of this study were (1) to identify tir and eae polymorphisms in diverse STEC O157 isolates derived from clinically ill humans and healthy cattle (the dominant zoonotic reservoir) and (2) to test any observed tir and eae polymorphisms for association with human (vs bovine) isolate source.ResultsFive polymorphisms were identified in a 1,627-bp segment of tir. Alleles of two tir polymorphisms, tir 255 T>A and repeat region 1-repeat unit 3 (RR1-RU3, presence or absence) had dissimilar distributions among human and bovine isolates. More than 99% of 108 human isolates possessed the tir 255 T>A T allele and lacked RR1-RU3. In contrast, the tir 255 T>A T allele and RR1-RU3 absence were found in 55% and 57%, respectively, of 77 bovine isolates. Both polymorphisms associated strongly with isolate source (p < 0.0001), but not by pulsed field gel electrophoresis type or by stx1 and stx2 status (as determined by PCR). Two eae polymorphisms were identified in a 2,755-bp segment of 44 human and bovine isolates; 42 isolates had identical eae sequences. The eae polymorphisms did not associate with isolate source.ConclusionPolymorphisms in tir but not eae predict the propensity of STEC O157 isolates to cause human clinical disease. The over-representation of the tir 255 T>A T allele in human-derived isolates vs the tir 255 T>A A allele suggests that these isolates have a higher propensity to cause disease. The high frequency of bovine isolates with the A allele suggests a possible bovine ecological niche for this STEC O157 subset.

Estimation of DNA sequence diversity in bovine cytokine genes

Abstract. DNA sequence variation provides the fundamental material for improving livestock through selection. In cattle, single nucleotide polymorphisms and small insertions/deletions (collectively referred to here as SNPs) have been identified in cytokine genes and scored in a reference population to determine linkage map positions. The aim of the present study was twofold: first, to estimate the SNP frequency in a reference population of beef cattle, and second, to determine cytokine haplotypes in a group of sires from commercial populations. Forty-five SNP markers in DNA segments from nine cytokine gene loci were analyzed in 26 reference parents. Comparison of all 52 haploid genomes at each PCR amplicon locus revealed an average of one SNP per 143 bp of sequence, whereas comparison of any two chromosomes identified heterozygous sites, on average, every 443 bp. The combination of these 45 SNP genotypes was sufficient to uniquely identify each of the 26 animals. The average number of haplotype alleles (4.4) per PCR amplicon (688 bp) and the percentage heterozygosity among founding parents (50%) were similar to those for microsatellite markers in the same population. For 49 sires from seven common breeds of beef cattle, SNP genotypes (1225 total) were obtained by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) at three amplicon loci. All three of the amplicon haplotypes were correctly deduced for each sire without the use of parent or progeny genotypes. The latter allows a wide range of genetic studies in commercial populations of cattle where genotypic information from relatives may not be available.