The Bovine HapMap Consortium* The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.T he emergence of modern civilization was accompanied by adaptation, assimilation, and interbreeding of captive animals. In cattle (Bos taurus), this resulted in the development of individual breeds differing in, for example, milk yield, meat quality, draft ability, and tolerance or resistance to disease and pests. However, despite mapping and diversity studies (1-5) and the identification of mutations affecting some quantitative phenotypes (6-8), the detailed genetic structure and history of cattle are not known.Cattle occur as two major geographic types, the taurine (humpless-European, African, and Asian) and indicine (humped-South Asian, and East African), which diverged >250 thousand years ago (Kya) (3). We sampled individuals representing 14 taurine (n = 376), three indicine (n = 73) (table S1), and two hybrid breeds (n = 48), as well as two individuals each of Bubalus quarlesi and Bubalus bubalis, which diverged from Bos taurus~1.25 to 2.0 Mya (9, 10). All breeds except Red Angus (n = 12) were represented by at least 24 individuals. We preferred individuals that were unrelated for ≥4 generations; however, each breed had one or two sire, dam, and progeny trios to allow assessment of genotype quality.Single-nucleotide polymorphisms (SNPs) that were polymorphic in many populations were primarily derived by comparing whole-genome sequence reads representing five taurine and one indicine breed to the reference genome assembly obtained from a Hereford cow (10) (table S2). This led to the ascertainment of SNPs with high minor allele frequencies (MAFs) within the discovery breeds (table S5). Thus, as expected, with trio progeny removed, SNPs discovered within the taurine breeds had higher average MAFs

Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high-density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (FROH ), which were compared with inbreeding coefficients estimated from pedigree information (FPED ) and using the genomic relationship matrix (FGRM ). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The FPED were moderately correlated with FROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with FGRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of FROH with FPED and FGRM . ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium-density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.

For about 10 000 years, farmers have been managing cattle, sheep, and goats in a sustainable way, leading to animals that are well adapted to the local conditions. About 200 years ago, the situation started to change dramatically, with the rise of the concept of breed. All animals from the same breed began to be selected for the same phenotypic characteristics, and reproduction among breeds was seriously reduced. This corresponded to a strong fragmentation of the initial populations. A few decades ago, the selection pressures were increased again in order to further improve productivity, without enough emphasis on the preservation of the overall genetic diversity. The efficiency of modern selection methods successfully increased the production, but with a dramatic loss of genetic variability. Many industrial breeds now suffer from inbreeding, with effective population sizes falling below 50. With the development of these industrial breeds came economic pressure on farmers to abandon their traditional breeds, and many of these have recently become extinct as a result. This means that genetic resources in cattle, sheep, and goats are highly endangered, particularly in developed countries. It is therefore important to take measures that promote a sustainable management of these genetic resources; first, by in situ preservation of endangered breeds; second, by using selection programmes to restore the genetic diversity of industrial breeds; and finally, by protecting the wild relatives that might provide useful genetic resources.

Genomic evaluations for 161,341 Holsteins were computed by using 311,725 of 777,962 markers on the Illumina BovineHD Genotyping BeadChip (HD). Initial edits with 1,741 HD genotypes from 5 breeds revealed that 636,967 markers were usable but that half were redundant. Holstein genotypes were from 1,510 animals with HD markers, 82,358 animals with 45,187 (50K) markers, 1,797 animals with 8,031 (8K) markers, 20,177 animals with 6,836 (6K) markers, 52,270 animals with 2,683 (3K) markers, and 3,229 nongenotyped dams (0K) with >90% of haplotypes imputable because they had 4 or more genotyped progeny. The Holstein HD genotypes were from 1,142 US, Canadian, British, and Italian sires, 196 other sires, 138 cows in a US Department of Agriculture research herd (Beltsville, MD), and 34 other females. Percentages of correctly imputed genotypes were tested by applying the programs findhap and FImpute to a simulated chromosome for an earlier population that had only 1,112 animals with HD genotypes and none with 8K genotypes. For each chip, 1% of the genotypes were missing and 0.02% were incorrect initially. After imputation of missing markers with findhap, percentages of genotypes correct were 99.9% from HD, 99.0% from 50K, 94.6% from 6K, 90.5% from 3K, and 93.5% from 0K. With FImpute, 99.96% were correct from HD, 99.3% from 50K, 94.7% from 6K, 91.1% from 3K, and 95.1% from 0K genotypes. Accuracy for the 3K and 6K genotypes further improved by approximately 2 percentage points if imputed first to 50K and then to HD instead of imputing all genotypes directly to HD. Evaluations were tested by using imputed actual genotypes and August 2008 phenotypes to predict deregressed evaluations of US bulls proven after August 2008. For 28 traits tested, the estimated genomic reliability averaged 61.1% when using 311,725 markers vs. 60.7% when using 45,187 markers vs. 29.6% from the traditional parent average. Squared correlations with future data were slightly greater for 16 traits and slightly less for 12 with HD than with 50K evaluations. The observed 0.4 percentage point average increase in reliability was less favorable than the 0.9 expected from simulation but was similar to actual gains from other HD studies. The largest HD and 50K marker effects were often located at very similar positions. The single-breed evaluation tested here and previous single-breed or multibreed evaluations have not produced large gains. Increasing the number of HD genotypes used for imputation above 1,074 did not improve the reliability of Holstein genomic evaluations.

This work reports on use of the recently described amplified fragment length polymorphism (AFLP) technology for DNA fingerprinting in cattle. The AFLP technology produces molecular markers through the high-stringency polymerase chain reaction (PCR)-amplification of restriction fragments that are ligated to synthetic adapters and amplified using primers, complementary to the adapters, which carry selective nucleotides at their 3' ends. While, for plants, the double digestion of genomic DNA with EcoRI and MseI is suggested, in mammals the enzyme combination EcoRI/TaqI produces clearer and more polymorphic AFLP patterns. In a sample of 47 Italian Holstein genotypes, 16 EcoRI/TaqI primer combinations identified 248 polymorphic bands in a species known for its low level of restriction polymorphism. In spite of the low information content carried by each AFLP polymorphism (average polymorphism information content = 0.31), the number of fragments revealed by each primer combination increased significantly the level of genetic information gained in each experiment. AFLP patterns are reproducible in independent experiments and polymorphic fragments segregate in cattle families according to Mendelian rules.