Five new recessive defects were discovered in Holsteins, Jerseys, and Brown Swiss by examining haplotypes that had a high population frequency but were never homozygous. The method required genotypes only from apparently normal individuals and not from affected embryos. Genotypes from the BovineSNP50 BeadChip (Illumina, San Diego, CA) were examined for 58,453 Holsteins, 5,288 Jerseys, and 1,991 Brown Swiss with genotypes in the North American database. Haplotypes with a length of ≤ 75 markers were obtained. Eleven candidate haplotypes were identified, with the earliest carrier born before 1980; 7 to 90 homozygous haplotypes were expected, but none were observed in the genomic data. Expected numbers were calculated using either the actual mating pattern or assuming random mating. Probability of observing no homozygotes ranged from 0.0002 for 7 to 10⁻⁴⁵ for 90 expected homozygotes. Phenotypic effects were confirmed for 5 of the 11 candidate haplotypes using 14,911,387 Holstein, 830,391 Jersey, and 68,443 Brown Swiss records for conception rate. Estimated effect for interaction of carrier service sire with carrier maternal grandsire ranged from -3.0 to -3.7 percentage points, which was slightly smaller than the -3.9 to -4.6 percentage points expected for lethal recessives but slightly larger than estimated effects for previously known lethal alleles of -2.5 percentage points for brachyspina and -2.9 percentage points for complex vertebral malformation. Conception rate was coded as a success only if the gestation went to term or the cow was confirmed to be pregnant. Estimated effect of carrier interaction for stillbirth rate based on 10,876,597 Holstein and 25,456 Jersey records was small. Thus, lethal effects may include conception, gestation, and stillbirth losses. Carrier frequency has been >20% for many years for the confirmed defect in Jerseys and is currently 16% for the defect in Brown Swiss. The 3 defects discovered in Holsteins have carrier frequencies of 2.7 to 6.4% in the current population. For previously known defects, map locations and lack of homozygotes were consistent with the literature and lethal recessive inheritance, but numbers of expected homozygotes for some were small because of low frequency. Very large genotypic and phenotypic data sets allow efficient detection of smaller and less frequent effects. Haplotype tests can help breeders avoid carrier matings for such defects and reduce future frequencies.

Meiotic recombination is an essential biological process that generates genetic diversity and ensures proper segregation of chromosomes during meiosis. From a large USDA dairy cattle pedigree with over half a million genotyped animals, we extracted 186,927 three-generation families, identified over 8.5 million maternal and paternal recombination events, and constructed sex-specific recombination maps for 59,309 autosomal SNPs. The recombination map spans for 25.5 Morgans in males and 23.2 Morgans in females, for a total studied region of 2,516 Mb (986 kb/cM in males and 1,085 kb/cM in females). The male map is 10% longer than the female map and the sex difference is most pronounced in the subtelomeric regions. We identified 1,792 male and 1,885 female putative recombination hotspots, with 720 hotspots shared between sexes. These hotspots encompass 3% of the genome but account for 25% of the genome-wide recombination events in both sexes. During the past forty years, males showed a decreasing trend in recombination rate that coincided with the artificial selection for milk production. Sex-specific GWAS analyses identified PRDM9 and CPLX1 to have significant effects on genome-wide recombination rate in both sexes. Two novel loci, NEK9 and REC114, were associated with recombination rate in both sexes, whereas three loci, MSH4, SMC3 and CEP55, affected recombination rate in females only. Among the multiple PRDM9 paralogues on the bovine genome, our GWAS of recombination hotspot usage together with linkage analysis identified the PRDM9 paralogue on chromosome 1 to be associated in the U.S. Holstein data. Given the largest sample size ever reported for such studies, our results reveal new insights into the understanding of cattle and mammalian recombination.

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.