J. W. Keele

Analysis of copy number variations among diverse cattle breeds

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.

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.

A second-generation linkage map of the sheep genome.

A genetic map of Ovis aries (haploid n = 27) was developed with 519 markers (504 microsatellites) spanning ∼3063 cM in 26 autosomal linkage groups and 127 cM (female specific) of the X Chromosome (Chr). Genotypic data were merged from the IMF flock (Crawford et al., Genetics 140, 703, 1995) and the USDA mapping flock. Seventy-three percent (370/504) of the microsatellite markers on the map are common to the USDA-ARS MARC cattle linkage map, with 27 of the common markers derived from sheep. The number of common markers per homologous linkage group ranges from 5 to 22 and spans a total of 2866 cM (sex average) in sheep and 2817 cM in cattle. Marker order within a linkage group was consistent between the two species with limited exceptions. The reported translocation between the telomeric end of bovine Chr 9 (BTA 9) and BTA 14 to form ovine Chr 9 is represented by a 15-cM region containing 5 common markers. The significant genomic conservation of marker order will allow use of linkage maps in both species to facilitate the search for quantitative trait loci (QTLs) in cattle and sheep.

Genome-wide association study of growth in crossbred beef cattle

Chromosomal regions harboring variation affecting cattle birth weight and BW gain to 1 yr of age were identified by marker association using the highly parallel BovineSNP50 BeadChip (50K) assay composed of 54,001 individual SNP. Genotypes were obtained from progeny (F(1); 590 steers) and 2-, 3-, and 4-breed cross grandprogeny (F(1)(2) = F(1) x F(1); 1,306 steers and 707 females) of 150 AI sires representing 7 breeds (22 sires per breed; Angus, Charolais, Gelbvieh, Hereford, Limousin, Red Angus, and Simmental). Genotypes and birth, weaning, and yearling BW records were used in whole-genome association analyses to estimate effects of individual SNP on growth. Traits analyzed included growth component traits: birth weight (BWT), 205-d adjusted birth to weaning BW gain (WG), 160-d adjusted postweaning BW gain (PWG); cumulative traits: 205-d adjusted weaning weight (WW = BWT + WG) and 365-d adjusted yearling weight (YW = BWT + WG + PWG); and indexes of relative differences between postnatal growth and birth weight. Modeled fixed effects included additive effects of calf and dam SNP genotype, year-sex-management contemporary groups, and covariates for calf and dam breed composition and heterosis. Direct and maternal additive polygenic effects and maternal permanent environment effects were random. Missing genotypes, including 50K genotypes of most dams, were approximated with a single-locus BLUP procedure from pedigree relationships and known 50K genotypes. Various association criteria were applied: stringent tests to account for multiple testing but with limited power to detect associations with small effects, and relaxed nominal P that may detect SNP associated with small effects but include excessive false positive associations. Genomic locations of the 231 SNP meeting stringent criteria generally coincided with described previously QTL affecting growth traits. The 12,425 SNP satisfying relaxed tests were located throughout the genome. Most SNP associated with BWT and postnatal growth affected components in the same direction, although detection of SNP associated with one component independent of others presents a possible opportunity for SNP-assisted selection to increase postnatal growth relative to BWT.

Sequence analysis of a rainbow trout cDNA library and creation of a gene index

Expressed sequence tag (EST) projects have produced extremely valuable resources for identifying genes affecting phenotypes of interest. A large-scale EST sequencing project for rainbow trout was initiated to identify and functionally annotate as many unique transcripts as possible. Over 45,000 5′ ESTs were obtained by sequencing clones from a single normalized library constructed using mRNA from six tissues. The production of this sequence data and creation of a rainbow trout Gene Index eliminating redundancy and providing annotation for these sequences will facilitate research in this species.

A small-insert bovine genomic library highly enriched for microsatellite repeat sequences

A bovine genomic phagemid library was constructed with randomly sheared DNA. Enrichment of this single-stranded DNA library with CA or GT primers resulted in 45% positive clones. The 14% of positive clones with (CA · GT)>12, and not containing flanking repetitive elements, were sequenced, and the efficiency of marker production was compared with random M13 bacteriophage libraries. Primer sequences and genotyping information are presented for 390 informative bovine microsatellite markers. The genomic frequency for 11 tri- and tetranucleotide repeats was estimated by hybridization to a lambda genomic library. Only GCT, GGT, and GGAT were estimated to have a frequency of >100 per genome. Enrichment of the phagemid library for these repeats failed to provide a viable source of microsatellite markers in the bovine. Comparison of map interval lengths between 100 markers from the enriched library prepared from randomly sheared DNA and M13 bacteriophage libraries prepared from Mbo1 restriction digests suggested no bias in skeletal genomic coverage based on source of small insert DNA. In conclusion, enrichment of the bovine phagemid library provides a sufficient source of microsatellites so that small repeat lengths and flanking repetitive sequences common in the bovine can be eliminated, resulting in a high percentage of informative markers.

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.

Porcine gene discovery by normalized cDNA-library sequencing and EST cluster assembly.

Genetic and environmental factors affect the efficiency of pork production by influencing gene expression during porcine reproduction, tissue development, and growth. The identification and functional analysis of gene products important to these processes would be greatly enhanced by the development of a database of expressed porcine gene sequence. Two normalized porcine cDNA libraries (MARC 1PIG and MARC 2PIG), derived respectively from embryonic and reproductive tissues, were constructed, sequenced, and analyzed. A total of 66,245 clones from these two libraries were 5?-end sequenced and deposited in GenBank. Cluster analysis revealed that within-library redundancy is low, and comparison of all porcine ESTs with the human database suggests that the sequences from these two libraries represent portions of a significant number of independent pig genes. A Porcine Gene Index (PGI), comprising 15,616 tentative consensus sequences and 31,466 singletons, includes all sequences in public repositories and has been developed to facilitate further comparative map development and characterization of porcine genes (http://www.tigr.org/tdb/ssgi/). The clones and sequences from these libraries provide a catalog of expressed porcine genes and a resource for development of high-density hybridization arrays for transcriptional profiling of porcine tissues. In addition, comparison of porcine ESTs with sequences from other species serves as a valuable resource for comparative map development. Both arrayed cDNA libraries are available for unrestricted public use.

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.

A physical map of the bovine genome

BackgroundCattle are important agriculturally and relevant as a model organism. Previously described genetic and radiation hybrid (RH) maps of the bovine genome have been used to identify genomic regions and genes affecting specific traits. Application of these maps to identify influential genetic polymorphisms will be enhanced by integration with each other and with bacterial artificial chromosome (BAC) libraries. The BAC libraries and clone maps are essential for the hybrid clone-by-clone/whole-genome shotgun sequencing approach taken by the bovine genome sequencing project.ResultsA bovine BAC map was constructed with HindIII restriction digest fragments of 290,797 BAC clones from animals of three different breeds. Comparative mapping of 422,522 BAC end sequences assisted with BAC map ordering and assembly. Genotypes and pedigree from two genetic maps and marker scores from three whole-genome RH panels were consolidated on a 17,254-marker composite map. Sequence similarity allowed integrating the BAC and composite maps with the bovine draft assembly (Btau3.1), establishing a comprehensive resource describing the bovine genome. Agreement between the marker and BAC maps and the draft assembly is high, although discrepancies exist. The composite and BAC maps are more similar than either is to the draft assembly.ConclusionFurther refinement of the maps and greater integration into the genome assembly process may contribute to a high quality assembly. The maps provide resources to associate phenotypic variation with underlying genomic variation, and are crucial resources for understanding the biology underpinning this important ruminant species so closely associated with humans.