M.E. Tooker

Genomic inbreeding and relationships among Holsteins, Jerseys, and Brown Swiss.

Genomic measures of relationship and inbreeding within and across breeds were compared with pedigree measures using genotypes for 43,385 loci of 25,219 Holsteins, 3,068 Jerseys, and 872 Brown Swiss. Adjustment factors allow genomic and pedigree relationships to match more closely within breeds and in multibreed populations and were estimated using means and regressions of genomic on pedigree relationships and allele frequencies in base populations. Correlations of genomic relationships with pedigree inbreeding were higher within each breed when an allele frequency of 0.5, rather than base population frequencies, was used, whereas correlations of average genomic relationships with average pedigree relationships and also reliabilities of genomic evaluations were higher using base population frequencies. Allele frequencies differed in the 3 breeds and were correlated by 0.65 to 0.67 when estimated from genotyped animals compared with 0.72 to 0.74 when estimated from breed base populations. The largest difference in allele frequency was between Holstein and the other breeds on chromosome Bos taurus autosome 4 near a gene affecting appearance of white skin patches (vitiligo) in humans. Each animals breed composition was predicted very accurately with a standard deviation of <3% using regressions on genotypes at all loci or less accurately with a standard deviation of <6% using subsets of loci. Genomic future inbreeding (half an animals mean genomic relationship to current animals of the same breed) was correlated by 0.75 to 0.94 with expected future inbreeding (half the average pedigree relationship). Correlations of both were slightly higher with parent averages than with genomic evaluations for net merit of young Holstein bulls. Thus, rates of increase in genomic and pedigree inbreeding per generation should be slightly reduced with genomic selection, in agreement with previous simulations. Genomic inbreeding and future inbreeding have been provided with individual genomic predictions since 2008. New methods to adjust pedigree and genomic relationship matrices so that they match may provide an improved basis for multibreed genomic evaluation. Positive definite matrices can be obtained by adjusting pedigree relationships for covariances among base animals within breed, whereas adjusting genomic relationships to match pedigree relationships can introduce negative eigenvalues. Pedigree relationship matrices ignore common ancestry shared by base animals within breed and may not approximate genomic relationships well in multibreed populations.

Multibreed genomic evaluations using purebred Holsteins, Jerseys, and Brown Swiss.

Multibreed models are currently used in traditional US Department of Agriculture (USDA) dairy cattle genetic evaluations of yield and health traits, but within-breed models are used in genomic evaluations. Multibreed genomic models were developed and tested using the 19,686 genotyped bulls and cows included in the official August 2009 USDA genomic evaluation. The data were divided into training and validation sets. The training data set comprised bulls that were daughter proven and cows that had records as of November 2004, totaling 5,331 Holstein, 1,361 Jersey, and 506 Brown Swiss. The validation data set had 2,508 Holstein, 413 Jersey, and 185 Brown Swiss bulls that were unproven (no daughter information) in November 2004 and proven by August 2009. A common set of 43,385 single nucleotide polymorphisms (SNP) was used for all breeds. Three methods of multibreed evaluation were investigated. Method 1 estimated SNP effects separately within breed and then applied those breed-specific SNP estimates to the other breeds. Method 2 estimated a common set of SNP effects from combined genotypes and phenotypes of all breeds. Method 3 solved for correlated SNP effects within each breed estimated jointly using a multitrait model where breeds were treated as different traits. Across-breed genomic predicted transmitting ability (GPTA) and within-breed GPTA were compared using regressions to predict the deregressed validation data. Method 1 worked poorly, and coefficients of determination (R(2)) were much lower using training data from a different breed to estimate SNP effects. Correlations between direct genomic values computed using training data from different breeds were less than 30% and sometimes negative. Across-breed GPTA from method 2 had higher R(2) values than parent average alone but typically produced lower R(2) values than the within-breed GPTA. The across-breed R(2) exceeded the within-breed R(2) for a few traits in the Brown Swiss breed, probably because information from the other breeds compensated for the small numbers of Brown Swiss training animals. Correlations between within-breed GPTA and across-breed GPTA ranged from 0.91 to 0.93. The multibreed GPTA from method 3 were significantly better than the current within-breed GPTA, and adjusted R(2) for protein yield (the only trait tested for method 3) were highest of all methods for all breeds. However, method 3 increased the adjusted R(2) by only 0.01 for Holsteins, ≤0.01 for Jerseys, and 0.01 for Brown Swiss compared with within-breed predictions.

Selection and management of DNA markers for use in genomic evaluation

To facilitate routine genomic evaluation, a database was constructed to store genotypes for 50,972 single nucleotide polymorphisms (SNP) from the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA). Multiple samples per animal are allowed. All SNP genotypes for a sample are stored in a single row. An indicator specifies whether the genotype for a sample was selected for use in genomic evaluation. Samples with low call rates or pedigree conflicts are designated as unusable. Among multiple samples that qualify for use in genomic evaluation, the one with the highest call rate is designated as usable. When multiple samples are stored for an animal, a composite is formed during extraction by using SNP genotypes from other samples to replace missing genotypes. To increase the number of SNP available, scanner output for approximately 19,000 samples was reprocessed. Any SNP with a minor allele frequency of > or = 1% for Holsteins, Jerseys, or Brown Swiss was selected, which was the primary reason that the number of SNP used for USDA genomic evaluations increased. Few parent-progeny conflicts (< or = 1%) and a high call rate (> or = 90%) were additional requirements that eliminated 2,378 SNP. Because monomorphic SNP did not degrade convergence during estimation of SNP effects, a single set of 43,385 SNP was adopted for all breeds. The use of a database for genotypes, detection of conflicts as genotypes are stored, online access for problem resolution, and use of a single set of SNP for genomic evaluations have simplified tracking of genotypes and genomic evaluation as a routine and official process.

Use of the Illumina Bovine3K BeadChip in dairy genomic evaluation1

Genomic evaluations using genotypes from the Illumina Bovine3K BeadChip (3K) became available in September 2010 and were made official in December 2010. The majority of 3K-genotyped animals have been Holstein females. Approximately 5% of male 3K genotypes and between 3.7 and 13.9%, depending on registry status, of female genotypes had sire conflicts. The chemistry used for the 3K is different from that of the Illumina BovineSNP50 BeadChip (50K) and causes greater variability in the accuracy of the genotypes. Approximately 2% of genotypes were rejected due to this inaccuracy. A single nucleotide polymorphism (SNP) was determined to be not usable for genomic evaluation based on percentage missing, percentage of parent-progeny conflicts, and Hardy-Weinberg equilibrium discrepancies. Those edits left 2,683 of the 2,900 3K SNP for use in genomic evaluations. The mean minor allele frequencies (MAF) for Holstein, Jersey, and Brown Swiss were 0.32, 0.28, and 0.29, respectively. Eighty-one SNP had both a large number of missing genotypes and a large number of parent-progeny conflicts, suggesting a correlation between call rate and accuracy. To calculate a genomic predicted transmitting ability (GPTA) the genotype of an animal tested on a 3K is imputed to the 45,187 SNP included in the current genomic evaluation based on the 50K. The accuracy of imputation increases as the number of genotyped parents increases from none to 1 to both. The average percentage of imputed genotypes that matched the corresponding actual 50K genotypes was 96.3%. The correlation of a GPTA calculated from a 3K genotype that had been imputed to 50K and GPTA from its actual 50K genotype averaged 0.959 across traits for Holsteins and was slightly higher for Jerseys at 0.963. The average difference in GPTA from the 50K- and 3K-based genotypes across trait was close to 0. The evaluation system has been modified to accommodate the characteristics of the 3K. The low cost of the 3K has greatly increased genotyping of females. Prior to the availability of the 3K (August 2010), female genotyping accounted for 38.7% of the genotyped animals. In the past year, the portion of total genotypes from females across all chip types rose to 59.0%.

Selecting sequence variants to improve genomic predictions for dairy cattle

BackgroundMillions of genetic variants have been identified by population-scale sequencing projects, but subsets of these variants are needed for routine genomic predictions or genotyping arrays. Methods for selecting sequence variants were compared using simulated sequence genotypes and real July 2015 data from the 1000 Bull Genomes Project.MethodsCandidate sequence variants for 444 Holstein animals were combined with high-density (HD) imputed genotypes for 26,970 progeny-tested Holstein bulls. Test 1 included single nucleotide polymorphisms (SNPs) for 481,904 candidate sequence variants. Test 2 also included 249,966 insertions-deletions (InDels). After merging sequence variants with 312,614 HD SNPs and editing steps, Tests 1 and 2 included 762,588 and 1,003,453 variants, respectively. Imputation quality from findhap software was assessed with 404 of the sequenced animals in the reference population and 40 randomly chosen animals for validation. Their sequence genotypes were reduced to the subset of genotypes that were in common with HD genotypes and then imputed back to sequence. Predictions were tested for 33 traits using 2015 data of 3983 US validation bulls with daughters that were first phenotyped after August 2011.ResultsThe average percentage of correctly imputed variants across all chromosomes was 97.2 for Test 1 and 97.0 for Test 2. Total time required to prepare, edit, impute, and estimate the effects of sequence variants for 27,235 bulls was about 1 week using less than 33 threads. Many sequence variants had larger estimated effects than nearby HD SNPs, but prediction reliability improved only by 0.6 percentage points in Test 1 when sequence SNPs were added to HD SNPs and by 0.4 percentage points in Test 2 when sequence SNPs and InDels were included. However, selecting the 16,648 candidate SNPs with the largest estimated effects and adding them to the 60,671 SNPs used in routine evaluations improved reliabilities by 2.7 percentage points.ConclusionsReliabilities for genomic predictions improved when selected sequence variants were added; gains were similar for simulated and real data for the same population, and larger than previous gains obtained by adding HD SNPs. With many genotyped animals, many data sources, and millions of variants, computing strategies must efficiently balance costs of imputation, selection, and prediction to obtain subsets of markers that provide the highest accuracy.

Determination of quantitative trait nucleotides by concordance analysis between quantitative trait loci and marker genotypes of US Holsteins

Experimental designs that exploit family information can provide substantial predictive power in quantitative trait nucleotide discovery projects. Concordance between quantitative trait locus genotype as determined by the a posteriori granddaughter design and marker genotype was determined for 30 trait-by-chromosomal segment effects segregating in the US Holstein population with probabilities of <10-20 to accept the null hypotheses of no segregating gene affecting the trait within the chromosomal segment. Genotypes for 83 grandsires and 17,217 sons were determined by either complete sequence or imputation for 3,148,506 polymorphisms across the entire genome; 471 Holstein bulls had a complete genome sequence, including 64 of the grandsires. Complete concordance was obtained only for stature on chromosome 14 and daughter pregnancy rate on chromosome 18. For each quantitative trait locus, effects of the 30 polymorphisms with highest concordance scores for the analyzed trait were computed by stepwise regression for predicted transmitting abilities of 26,750 bulls with progeny test and imputed genotypes. Effects for stature on chromosome 11, daughter pregnancy rate on chromosome 18, and protein percentage on chromosome 20 met 3 criteria: complete or almost complete concordance, nominal significance of the polymorphism effect after correction for all other polymorphisms, and marker coefficient of determination >40% of total multiple-regression coefficient of determination for the 30 polymorphisms with highest concordance. An intronic variant marker on chromosome 5 at 93,945,738 bp explained 7% of variance for fat percentage and 74% of total multiple-marker regression variance but was concordant for only 24 of 30 families. The missense polymorphism Phe279Tyr in GHR at 31,909,478 bp on chromosome 20 was confirmed as the causative mutation for fat and protein concentration. For effect on fat percentage on chromosome 14, 12 additional missense polymorphisms were found that had almost complete concordance with the suggested causative polymorphism (missense mutation Ala232Glu in DGAT1). The only polymorphism found likely to improve predictive power for genomic evaluation of dairy cattle was on chromosome 15; that polymorphism had a frequency of 0.45 for the allele with economically positive effects on all production traits.