Cow genotyping strategies for genomic selection in a small dairy cattle population.

This study compares how different cow genotyping strategies increase the accuracy of genomic estimated breeding values (EBV) in dairy cattle breeds with low numbers. In these breeds, few sires have progeny records, and genotyping cows can improve the accuracy of genomic EBV. The Guernsey breed is a small dairy cattle breed with approximately 14,000 recorded individuals worldwide. Predictions of phenotypes of milk yield, fat yield, protein yield, and calving interval were made for Guernsey cows from England and Guernsey Island using genomic EBV, with training sets including 197 de-regressed proofs of genotyped bulls, with cows selected from among 1,440 genotyped cows using different genotyping strategies. Accuracies of predictions were tested using 10-fold cross-validation among the cows. Genomic EBV were predicted using 4 different methods: (1) pedigree BLUP, (2) genomic BLUP using only bulls, (3) univariate genomic BLUP using bulls and cows, and (4) bivariate genomic BLUP. Genotyping cows with phenotypes and using their data for the prediction of single nucleotide polymorphism effects increased the correlation between genomic EBV and phenotypes compared with using only bulls by 0.163±0.022 for milk yield, 0.111±0.021 for fat yield, and 0.113±0.018 for protein yield; a decrease of 0.014±0.010 for calving interval from a low base was the only exception. Genetic correlation between phenotypes from bulls and cows were approximately 0.6 for all yield traits and significantly different from 1. Only a very small change occurred in correlation between genomic EBV and phenotypes when using the bivariate model. It was always better to genotype all the cows, but when only half of the cows were genotyped, a divergent selection strategy was better compared with the random or directional selection approach. Divergent selection of 30% of the cows remained superior for the yield traits in 8 of 10 folds.

Genomic evaluation, breed identification, and population structure of Guernsey cattle in North America, Great Britain, and the Isle of Guernsey.

As of December 2015, 2,376 Guernsey bulls and cows had genotypes from collaboration between the United States, Canada, the United Kingdom, and the Isle of Guernsey. Of those, 439 bulls and 504 cows had traditional US evaluations, which provided sufficient data to justify investigation of the possible benefits of genomic evaluation for the Guernsey breed. Evaluation accuracy was assessed using a traditional 4-yr cutoff study. Twenty-two traits were analyzed (5 yield traits, 3 functional traits, and 14 conformation traits). Mean reliability gain over that for parent average was 16.8 percentage points across traits, which compares with 8.2, 18.5, 20.0, and 32.6 percentage points reported for Ayrshires, Brown Swiss, Jerseys, and Holsteins, respectively. Highest Guernsey reliability gains were for rump width (44.5 percentage points) and dairy form (40.5 percentage points); lowest gains were for teat length (1.9 percentage points) and rear legs (side view) (2.3 percentage points). Slight reliability losses (1.5 to 4.5 percentage points) were found for udder cleft, final score, and udder depth as well as a larger loss (13.6 percentage points) for fore udder attachment. Twenty-one single nucleotide polymorphisms were identified for Guernsey breed determination and can be used in routine genotype quality control to confirm breed and identify crossbreds. No haplotypes that affect fertility were identified from the current data set. Principal component analysis showed some divergence of US and Isle of Guernsey subpopulations. However, the overlap of US, Canadian, UK, and Isle of Guernsey subpopulations indicated the presence of gene flow, and the similarities in the subpopulations supports a common genomic evaluation system across the regions.

Increasing the number of single nucleotide polymorphisms used in genomic evaluation of dairy cattle.

GeneSeek (Neogen Corp., Lexington, KY) designed a new version of the GeneSeek Genomic Profiler HD BeadChip for Dairy Cattle, which originally had >77,000 single nucleotide polymorphisms (SNP). A set of >140,000 SNP was selected that included all SNP on the existing GeneSeek chip, all SNP used in US national genomic evaluations, SNP that were possible functional mutations, and other informative SNP. Because SNP with a lower minor allele frequency might track causative variants better, 30,000 more SNP were selected from the Illumina BovineHD Genotyping BeadChip (Illumina Inc., San Diego, CA) by choosing SNP to maximize differences in minor allele frequency between a SNP being considered for the new chip and the 2 SNP that flanked it. Single-gene tests were included if their location was known and bioinformatics indicated relevance for dairy cattle. To determine which SNP from the new chip should be included in genomic evaluations, genotypes available from chips already in use were used to impute and evaluate the SNP set. Effects for 134,511 usable SNP were estimated for all breed-trait combinations; SNP with the largest absolute values for effects were selected (5,000 for Holsteins, 1,000 for Jerseys, and 500 each for Brown Swiss and Ayrshires for each trait). To increase overlap with the 60,671 SNP currently used for genomic evaluation, 12,094 more SNP with the largest effects were added. After removing SNP with many parent-progeny conflicts, 84,937 SNP remained. Three cutoff studies were conducted with 3 SNP sets to determine reliability gain over that for parent average when evaluations based on August 2011 data were used to predict December 2014 performance. Across all traits, mean Holstein reliability gains were 32.5, 33.4, and 32.0 percentage points for 60,671, 84,937, and 134,511 SNP, respectively. After genotypes from the new chip became available, the proposed set was reduced from 84,937 to 77,321 SNP to remove SNP that were not included during manufacture, reduce computing time, and improve imputation performance. The set of 77,321 SNP was evaluated using August 2011 data to predict April 2015 performance. Reliability gain over 60,671 SNP was 1.4 percentage points across traits for Holsteins. Improvement over 84,937 SNP was partially the result of 4mo of additional data and genotypes from the new chip. Revision of the SNP set used for genomic evaluation is expected to be an ongoing process to increase evaluation accuracy.

Using genomics to enhance selection of novel traits in North American dairy cattle.

The objectives of this paper were to briefly review progress in the genetic evaluation of novel traits in Canada and the United States, assess methods to predict selection accuracy based on cow reference populations, and illustrate how the use of indicator traits could increase genomic selection accuracy. Traits reviewed are grouped into the following categories: udder health, hoof health, other health traits, feed efficiency and methane emissions, and other novel traits. The status of activities expected to lead to national genetic evaluations is indicated for each group of traits. For traits that are more difficult to measure or expensive to collect, such as individual feed intake or immune response, the development of a cow reference population is the most effective approach. Several deterministic methods can be used to predict the reliability of genomic evaluations based on cow reference population size, trait heritability, and other population parameters. To provide an empirical validation of those methods, predicted accuracies were compared with observed accuracies for several cow reference populations and traits. Reference populations of 2,000 to 20,000 cows were created through random sampling of genotyped Holstein cows in Canada and the United States. The effects of single nucleotide polymorphisms (SNP) were estimated from those cow records, after excluding the dams of validation bulls. Bulls that were first progeny tested in 2013 and 2014 were then used to carry out a validation and estimate the observed accuracy of genomic selection based on those SNP effects. Over the various cow population sizes and traits considered in the study, even the best prediction methods were found, on average, to either under-evaluate observed accuracy by 0.20 or over-evaluate it by 0.22, depending on the approach used to estimate the number of independently segregating chromosome segments. In some instances, differences between observed and predicted accuracies were as large as 0.47. Indicator traits can be very useful for the selection of novel traits. To illustrate this, protein yield, body weight, and mid-infrared data were used as indicator traits for feed efficiency. Using those traits in conjunction with 5,000 cow records for dry matter intake increased the reliability of genomic predictions for young animals from 0.20 to 0.50.

Short communication: Improving accuracy of Jersey genomic evaluations in the United States and Denmark by sharing reference population bulls.

The effect on prediction accuracy for Jersey genomic evaluations of Danish and US bulls from using a larger reference population was assessed. Each country contributed genotypes from 1,157 Jersey bulls to the reference population of the other. Data were separated into reference (US only, Danish only, and combined US-Danish) and validation (US only and Danish only) populations. Depending on trait (milk, fat, and protein yields and component percentages; productive life; somatic cell score; daughter pregnancy rate; 14 conformation traits; and net merit), the US reference population included 2,720 to 4,772 bulls and cows with traditional evaluations as of August 2009; the Danish reference population included 635 to 996 bulls. The US validation population included 442 to 712 bulls that gained a traditional evaluation between August 2009 and December 2013; the Danish validation population included 105 to 196 bulls with multitrait across-country evaluations on the US scale by December 2013. Genomic predicted transmitting abilities (GPTA) were calculated on the US scale using a selection index that combined direct genomic predictions with either traditional predicted transmitting ability for the reference population or traditional parent averages (PA) for the validation population and a traditional evaluation based only on genotyped animals. Reliability for GPTA was estimated from the reference population and August 2009 traditional PA and PA reliability. For prediction of December 2013 deregressed daughter deviations on the US scale, mean August 2009 GPTA reliability for Danish validation bulls was 0.10 higher when based on the combined US-Danish reference population than when the reference population included only Danish bulls; for US validation bulls, mean reliability increased by 0.02 when Danish bulls were added to the US reference population. Exchanging genotype data to increase the size of the reference population is an efficient approach to increasing the accuracy of genomic prediction when the reference population is small.

Short communication: Analysis of genomic predictor population for Holstein dairy cattle in the United States--Effects of sex and age.

Increased computing time for the ever-growing predictor population and linkage decay between the ancestral population and current animals have become concerns for genomic evaluation systems. The effects on reliability of US genomic evaluations from including cows and bulls in the Holstein predictor population and also from excluding older bulls from the predictor population were examined. Holstein data collected for December 2013 US genomic evaluations were used in cutoff studies to determine reliability gains, regression coefficients, and bias for 5 yield, 3 fitness, 2 fertility, and 18 conformation traits. Three predictor populations were examined based on animal sex: 30,852 cows with traditional evaluations as of August 2012, 21,883 bulls with traditional evaluations as of August 2012, and a combined group of all bulls and cows. Three subsets of the bull predictor population were examined to determine effect of age: bulls born before 1996 excluded (25% of bulls excluded), bulls born before 2001 excluded (50%), and bulls born before 2005 excluded (75%). The validation set for all predictor populations was either bulls or cows first receiving a traditional evaluation between August 2012 and December 2013. Across all traits, the addition of cows to the bull predictor population increased reliability gains by 0.4 percentage points for validation bulls and 4.4 points for validation cows. Across all traits, excluding bulls born before 1996 from the bull-only predictor population decreased gains in genomic reliability by 1.8 percentage points. For 19 of 28 traits, excluding bulls born before 2005 from the predictor population resulted in lower bias in genomic evaluations of validation bulls. Although the contribution of cows and older bulls to improved accuracy of US genomic evaluations is small, a plateau of achievable gain has not yet been reached.

Technical note: Rapid calculation of genomic evaluations for new animals.

A method was developed to calculate preliminary genomic evaluations daily or weekly before the release of official monthly evaluations by processing only newly genotyped animals using estimates of single nucleotide polymorphism effects from the previous official evaluation. To minimize computing time, reliabilities and genomic inbreeding are not calculated, and fixed weights are used to combine genomic and traditional information. Correlations of preliminary and September official monthly evaluations for animals with genotypes that became usable after the extraction of genotypes for August 2014 evaluations were >0.99 for most Holstein traits. Correlations were lower for breeds with smaller population size. Earlier access to genomic evaluations benefits producers by enabling earlier culling decisions and genotyping laboratories by making workloads more uniform across the month.

Decreased efficacy of twice-weekly intravaginal dehydroepiandrosterone on vulvovaginal atrophy.

OBJECTIVE: While daily intravaginal administration of 0.50% (6.5 mg) dehydroepiandrosterone (DHEA, prasterone) for 12 weeks has shown clinically and statistically significant effects on moderate to severe (MS) dyspareunia as the most bothersome symptom (MBS), the present study analyzes the effect of a reduced dosing regimen on MBS vaginal dryness. METHOD: Daily intravaginal 0.50% prasterone for 2 weeks followed by twice weekly for 10 weeks versus placebo. RESULTS: Maximal beneficial changes in vaginal parabasal and superficial cells and pH were observed at 2 weeks as observed for intravaginal 10 µg estradiol (E2). This was followed by a decrease or lack of efficacy improvement after switching to twice-weekly dosing. The decrease in percentage of parabasal cells, increase in percentage of superficial cells and decrease in vaginal pH were all highly significant (p < 0.0001 to 0.0002 over placebo) at 12 weeks. In parallel, the statistical significance over placebo (p value) on MBS vaginal dryness at 6 weeks was 0.09 followed by an increase to 0.198 at 12 weeks. For MBS dyspareunia, the p value of 0.008 at 6 weeks was followed by a p value of 0.077 at 12 weeks, thus illustrating a decrease of efficacy at the lower dosing regimen. The improvements of vaginal secretions, color, epithelial integrity and epithelial surface thickness were observed at a p value < 0.01 or 0.05 over placebo at 2 weeks, with a similar or loss of statistical difference compared to placebo at later time intervals. No significant adverse event was observed. Vaginal discharge related to the melting of Witepsol was reported in 1.8% of subjects. CONCLUSION: The present data show that daily dosing with 0.50% DHEA for 2 weeks followed by twice-weekly dosing is a suboptimal treatment of the symptoms/signs of vulvovaginal atrophy resulting from a substantial loss of the efficacy achieved at daily dosing.

Genomic evaluation, breed identification, and discovery of a haplotype affecting fertility for Ayrshire dairy cattle.

Genomic evaluations of dairy cattle in the United States have been available for Brown Swiss, Holsteins, and Jerseys since 2009. As of January 2013, 1,023 Ayrshires had genotypes in the North American database. Evaluation accuracy was assessed using genomic evaluations based on 646 bulls with 2008 traditional evaluations to predict daughter performance of up to 180 bulls in 2012. Mean gain in reliability over parent average for all traits was 8.2 percentage points. The highest gains were for protein yield (16.9 percentage points), milk yield (16.6 percentage points), and stature (16.2 percentage points). Twelve single nucleotide polymorphisms were useful for Ayrshire breed determination. Fewer breed-determining SNP were available for Ayrshires than for Holsteins, Jerseys, and Brown Swiss because of the similarity of Ayrshires and Holsteins. A haplotype that affects fertility was identified on chromosome 17 and traces back in the genotyped population to the bull Selwood Betty's Commander (born in 1953). The haplotype carrier frequency for genotyped Ayrshires was 26.1%. Sire conception rate was decreased by 4.3 ± 2.5 percentage points for carriers of the haplotype as determined by 618 matings of carrier sire by carrier maternal grandsire. Genomic evaluations for Ayrshires were officially implemented in the United States in April 2013.

A single nucleotide polymorphism in the dopamine receptor D2 gene may be informative for resistance to fescue toxicosis in angus-based cattle.

Fescue toxicosis (FT) reduces beef animal growth and fertility. Animals afflicted with FT typically have decreased circulating prolactin concentrations and thicker summer hair coats. Preliminary experiments examined the informativeness of a novel Dopamine Receptor 2 (DRD2) G/A SNP for resistance to FT. Steers grazed tall fescue containing a toxic (E+) or non-toxic (NTE) strain of endophyte. Decreased serum prolactin concentrations were observed in GG steers in May compared to AA steers when grazing E+ pastures (P < 0.02). In a second study, GG steers had decreased prolactin concentrations (P = 0.004) and increased hair coat scores (P = 0.01) relative to AA steers when grazing E+ pastures. Allele and genotypic frequencies were different (P = 0.016 and 0.026, respectively) between spring-calving and fall-calving herds grazing E+ pastures, such that the A allele and the AA genotype were more prevalent in spring-calving herds, suggesting active selection for the A allele. Regardless of calving season, AA heifers tended toward fewer days to first calf (733.6 ± 4.4 d) than did GG heifers (756.6 ± 9.2 days; P = 0.055). These results suggest that the DRD2 SNP may have use in selecting animals resistant to FT.

Short communication: relationship of call rate and accuracy of single nucleotide polymorphism genotypes in dairy cattle.

Call rates on both a single nucleotide polymorphism (SNP) basis and an animal basis are used as measures of data quality and as screening tools for genomic studies and evaluations of dairy cattle. To investigate the relationship of SNP call rate and genotype accuracy for individual SNP, the correlation between percentages of missing genotypes and parent-progeny conflicts for each SNP was calculated for 103,313 Holsteins. Correlations ranged from 0.14 to 0.38 for the BovineSNP50 and BovineLD (Illumina Inc., San Diego, CA) and GeneSeek Genomic Profiler (Neogen Corp., Lincoln, NE) chips, with lower correlations for newer chips. For US genomic evaluations, genotypes are excluded for animals with a call rate of <90% across autosomal SNP or <80% across X-specific SNP. Mean call rate for 220,175 Holstein, Jersey, and Brown Swiss genotypes was 99.6%. Animal genotypes with a call rate of ≤99% were examined from the US Department of Agriculture genotype database to determine how genotype call rate is related to accuracy of calls on an animal basis. Animal call rate was determined from SNP used in genomic evaluation and is the number of called autosomal and X-specific SNP genotypes divided by the number of SNP from that type of chip. To investigate the relationship of animal call rate and parentage validation, conflicts between a genotyped animal and its sire or dam were determined through a duo test (opposite homozygous SNP genotypes between sire and progeny; 1,374 animal genotypes) and a trio test (also including conflicts with dam and heterozygous SNP genotype for the animal when both parents are the same homozygote; 482 animal genotypes). When animal call rate was ≤ 80%, parentage validation was no longer reliable with the duo test. With the trio test, parentage validation was no longer reliable when animal call rate was ≤ 90%. To investigate how animal call rate was related to genotyping accuracy for animals with multiple genotypes, concordance between genotypes for 1,216 animals that had a genotype with a call rate of ≤ 99% (low call rate) as well as a genotype with a call rate of >99% (high call rate) were calculated by dividing the number of identical SNP genotype calls by the number of SNP that were called for both genotypes. Mean concordance between low- and high-call genotypes was >99% for a low call rate of >90% but decreased to 97% for a call rate of 86 to 90% and to 58% for a call rate of <60%. Edits on call rate reduce the use of incorrect SNP genotypes to calculate genomic evaluations.

Confirmation and discovery of maternal grandsires and great-grandsires in dairy cattle.

Selection, mating, and improvement of dairy animals have required accurate pedigrees. Genomic tools allow paternal ancestors to be easily confirmed or discovered because most sires are genotyped for many markers, but maternal ancestors are more difficult to discover because most female ancestors are not genotyped. Three methods to discover maternal grandsires (MGS) were developed and compared. Conflicts were counted one single nucleotide polymorphism (SNP) at a time between genotypes of the animal and potential MGS (duo method) or also using the sire's genotype (trio method). Alternatively, haplotypes of a potential MGS were matched to the animal's maternal haplotype, obtained by using linkage across loci (HAP method). The duo and trio methods can be performed as soon as a genotype is received because no imputation is required. The HAP method improved accuracy because genotypes with 2,683 (3 K) SNP were imputed to the 45,187 (50K) SNP used for genomic evaluation. The HAP method was tested using modified pedigrees with 5% of true MGS replaced by a random genotyped bull from the same birth year and 5% of MGS set to missing for 4,134 Holsteins, 552 Jerseys, and 142 Brown Swiss that had confirmed, genotyped sires. Those same animals were used to test the duo and trio methods, except that some animals had multiple genotypes and imputed dams were excluded. Accuracy measured how often the correct MGS was selected from among 12,152 genotyped Holstein, 2,265 Jersey, and 1,605 Brown Swiss potential MGS. Accuracies were 61, 60, and 65%, respectively, with the duo method; 95, 91, and 94% with the trio method; and 97, 95, and 97% with the HAP method. Accuracy of the duo method was poor (only 52% for animals genotyped with 3 K and 65% with 50K) because additional information from the paternal genotype is not used. Accuracy of the trio method was 97% with 50K but only 78% with 3K because the missing SNP were not imputed. Accuracy of the HAP method was 94% with 3 K genotypes, 98% with 50K, and 92% with nongenotyped, imputed dams. When the HAP method was extended to great-grandsires, the accuracy of maternal great-grandsire discovery was 92% for 652 Holsteins, 95% for 33 Jerseys, and 85% for 20 Brown Swiss. Accuracy was even higher using simulated genotypes. Because most dairy bulls over several generations have been genotyped, percentages of haplotypes shared with candidate males can accurately confirm, correct, or discover the sires, MGS, and even more distant ancestors of most animals.

Technical note: Characteristics and use of the Illumina BovineLD and GeneSeek Genomic Profiler low-density bead chips for genomic evaluation.

The GeneSeek Genomic Profiler (GGP) BeadChip (GeneSeek, Lincoln, NE), which became available commercially in February 2012, is based on the Illumina BovineLD Genotyping BeadChip (Illumina Inc., San Diego, CA), with 1,745 additional single nucleotide polymorphisms (SNP) for genomic evaluation and SNP for proprietary single-gene tests. The BovineLD chip with 6,909 SNP, which replaced the Illumina GoldenGate Bovine3K Genotyping BeadChip, has been available since October 2011. The GGP's additional SNP for genomic evaluation were selected to improve imputation by filling SNP gaps on chromosomes and including more Bovine3K SNP than were on the BovineLD chip and to impute microsatellite alleles to facilitate parentage validation. The SNP for single-gene tests were included to minimize the number of separate tests required for those genes, particularly for bulls. The September 2012 US national genomic evaluation included genotypes from BovineLD and GGP chips for 82,510 animals. For those data, BovineLD and GGP performance was similar. The call rate for SNP on these chips that were used in genomic evaluation was 99.6%. The 9 Y-chromosome SNP in common on the BovineLD and GGP chips were highly effective in sex validation (call rate of 99% for males and 0.01% for females). For both chips, the rate of parent-progeny conflicts on a SNP basis (≤ 0.004%) was similar to that for SNP on the Illumina BovineSNP50 Genotyping BeadChip. Imputation accuracy for 45,187 BovineSNP50 SNP averaged 99.4% for H olsteins. Imputation accuracy was slightly higher for the GGP chip compared with the BovineLD chip because of its additional SNP. Reliability for genomic evaluations using BovineLD and GGP genotypes was 3 percentage points higher than that for Bovine3K genotypes.

Technical note: adjustment of all cow evaluations for yield traits to be comparable with bull evaluations.

Traditional evaluations of cows with genotypes have been adjusted since April 2010 to be comparable with evaluations of bulls so that their value for estimation of single nucleotide polymorphism effects in genomic evaluation programs would be improved. However, that adjustment made them not comparable with traditional evaluations of nongenotyped cows. To create an adjustment for all cows with an evaluation based on US data, Mendelian sampling, which is the difference between predicted transmitting ability (PTA) and parent average (PA), was calculated for milk, fat, and protein yields and divided by a deregression factor. Standard deviations for the deregressed Mendelian sampling (DMS) were grouped by reliability with PA contribution removed (REL(no PA)). A multiplicative adjustment to reduce the DMS standard deviation for cows so that it would be the same as for bulls with similar REL(no PA) was represented as a linear function of REL(no PA). Mean cow PA by birth year was subtracted from individual bull and cow PA to create within-year PA deviation groups, and mean DMS was calculated by PA deviation group. Means decreased for bulls and increased for cows with increasing deviation. The differences were fit by linear regression on PA deviation and used to adjust cow DMS. The adjustment reduced PTA of cows with a high PA and increased PTA of cows with a low PA but did not change estimated genetic trend because adjustment was within birth year. The adjustment also reduced variance of cow evaluations within birth year. Traditional evaluations of genotyped cows with a REL(no PA) of ≥55% were further adjusted so that the difference between those evaluations and direct genomic values calculated using only bulls as predictors was similar to that for bulls. The second adjustment was small compared with a 2010 adjustment and, therefore, had little effect on the comparability of evaluations for genotyped and nongenotyped cows. Cows with converted evaluations from other countries were excluded from the predictor population, and their converted evaluations were adjusted so that the difference between their mean PTA and direct genomic value was the same as the corresponding difference for bulls. For cows with converted evaluations, the adjustment amount differed depending on REL(no PA) (<55% or ≥55%). The new adjustment was implemented by USDA in April 2011 and permits a fairer comparison of estimated genetic merit between nongenotyped and genotyped cows.

Use of the Illumina Bovine3K BeadChip in dairy genomic evaluation.

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%.

Technical note: adjustment of traditional cow evaluations to improve accuracy of genomic predictions.

Genomic evaluations are calculated using deregressed predicted transmitting abilities (PTA) from traditional evaluations to estimate effects of single nucleotide polymorphisms. The direct genomic value (sum of an animal's marker effects) should be consistent with traditional PTA, which is the case for bulls. However, traditional PTA of yield traits (milk, fat, and protein) for genotyped cows are higher than their direct genomic values. To ensure that characteristics of cow PTA for yield traits were more similar to those for bull PTA, mean and variance of cow Mendelian sampling (PTA minus parent average) were adjusted to be similar to those of bulls. The same adjustments were used for all genotyped cows in a breed. To determine gains in reliabilities, predictions were made for bulls with August 2010 evaluations that did not have traditional evaluations in August 2006. By adjusting cow PTA and parent averages of genotyped animals, Holstein and Jersey regressions of August 2010 deregressed PTA on genomic evaluations based on August 2006 data became closer to 1 for the adjusted predictor population compared with the unadjusted predictor population. Evaluation bias was decreased for Holsteins when the predictor population was adjusted. Mean gain in reliability over parent average increased 3.5 percentage points across yield traits for Holsteins and 0.9 percentage points for Jerseys when the predictor population was adjusted. The accuracy of genomic evaluations for Holsteins and Jerseys was increased through better use of information from cows.

The genomic evaluation system in the United States: past, present, future.

Implementation of genomic evaluation has caused profound changes in dairy cattle breeding. All young bulls bought by major artificial insemination organizations now are selected based on such evaluation. Evaluation reliability can reach approximately 75% for yield traits, which is adequate for marketing semen of 2-yr-old bulls. Shortened generation interval from using genomic evaluations is the most important factor in increasing the rate of genetic improvement. Genomic evaluations are based on 42,503 single nucleotide polymorphisms (SNP) genotyped with technology that became available in 2007. The first unofficial USDA genomic evaluations were released in 2008 and became official for Holsteins, Jerseys, and Brown Swiss in 2009. Evaluation accuracy has increased steadily from including additional bulls with genotypes and traditional evaluations (predictor animals). Some of that increase occurs automatically as young genotyped bulls receive a progeny test evaluation at 5 yr of age. Cow contribution to evaluation accuracy is increased by decreasing mean and variance of their evaluations so that they are similar to bull evaluations. Integration of US and Canadian genotype databases was critical to achieving acceptable initial accuracy and continues to benefit both countries. Genotype exchange with other countries added predictor bulls for Brown Swiss. In 2010, a low-density chip with 2,900 SNP and a high-density chip with 777,962 SNP were released. The low-density chip has increased greatly the number of animals genotyped and is expected to replace microsatellites in parentage verification. The high-density chip can increase evaluation accuracy by better tracking of loci responsible for genetic differences. To integrate information from chips of various densities, a method to impute missing genotypes was developed based on splitting each genotype into its maternal and paternal haplotypes and tracing their inheritance through the pedigree. The same method is used to impute genotypes of nongenotyped dams based on genotyped progeny and mates. Reliability of resulting evaluations is discounted to reflect errors inherent in the process. Further increases in evaluation accuracy are expected because of added predictor animals and more SNP. The large population of existing genotypes can be used to evaluate new traits; however, phenotypic observations must be obtained for enough animals to allow estimation of SNP effects with sufficient accuracy for application to the general population.

Differences among methods to validate genomic evaluations for dairy cattle.

Two methods of testing predictions from genomic evaluations were investigated. Data used were from the August 2006 and April 2010 official USDA genetic evaluations of dairy cattle. The training data set consisted of both cows and bulls that were proven (had own or daughter information) as of August 2006 and included 8,022, 1,959, and 1,056 Holsteins, Jerseys, and Brown Swiss, respectively. The validation data set consisted of bulls that were unproven as of August 2006 and were proven by April 2010 with 2,653, 411, and 132 Holsteins, Jerseys, and Brown Swiss for the production traits. Method 1 used the training animal's predicted transmitting ability (PTA) from August of 2006. Method 2 used the training animal's April 2010 PTA to estimate single nucleotide polymorphism effects. Both methods were tested using several regressions with the same validation animals. In both cases, the validation animals were tested using the deregressed April 2010 PTA. All traits that had genomic evaluations from the official USDA April 2010 genetic evaluations were tested. Results included bias, differences from expected regressions (calculated using selection intensities), and the coefficients of determination. The genomic information increased the predictive ability for most of the traits in all of the breeds. The 2 methods of testing resulted in some differences that would affect interpretation of results. The coefficient of determination was higher for all traits using method 2. This was the expected result as the data were not independent because evaluations of the validation bulls contributed to their sires' evaluations. The regression coefficients from method 2 were often higher than the regression coefficients from method 1. Many traits had regression coefficients that were higher than 2 standard deviations from the expected regressions when using method 2. This was partially due to the lack of independence of the training and validation data sets. Most traits did have some level of bias in the prediction equations, regardless of breed. The use of method 1 made it possible to evaluate the increased accuracy in proven first-crop bull evaluations by using genomic information. Proven first-crop bulls had an increase in accuracy from the addition of genomic information. It is advised to use method 1 for validation of genomic evaluations.

Effects of weaning and syndyphalin-33 on expression of melanocortinergic appetite-regulating genes in swine.

Syndyphalin-33 (SD-33) increases feed intake in sheep and recently weaned pigs. To assess the effects of SD-33 on hypothalamic gene expression, hypothalami were collected from unweaned pigs (n=19; 21±3 d of age) on day 0. Remaining pigs received an intramuscular injection of 0.5 µmole/kg SD-33 (SD) or saline (VEH) and weaned into individual pens. On days 1, 4, and 7 after weaning, hypothalami were collected from subsets of pigs (n=8 or 9) within each treatment group. Expression of µ-opioid receptor (MOR) was less in SD pigs than in VEH pigs on day 1 and day 4, suggesting down-regulation of the receptor by SD-33. Expression of hypothalamic melanocortin 4 receptor (MC4R) at 1 d after weaning was increased in VEH pigs (but not SD pigs) relative to levels before weaning. Expression of AGRP was not significantly altered by weaning or treatment at 1 d after weaning. At 4 d after weaning, expression of AGRP was greater in SD pigs than in VEH pigs, but at day 7 expression was less in SD pigs than in VEH pigs. A strong positive correlation was noted between expression levels of MOR and MC4R across treatment and time. Treatment with SD-33 appeared to partially abrogate the effects of weaning on expression of two key appetite-regulating genes within 24 h. Effects of SD-33 appear to be mediated at least in part by the µ-opioid receptor and include actions on the melanocortinergic pathway.

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.