New mutation within a common haplotype is associated with calf muscle weakness in Holsteins

A recessive haplotype resulting in elevated calf mortality but with apparent incomplete penetrance was previously linked to the end of chromosome 16 (78.7-80.7 Mbp). Genotype analysis of 5.6 million Holsteins indicated that the haplotype was common and traced back to 1952, with a key ancestor born in 1984 (HOUSA1964484, Southwind) identified from chip genotypes as homozygous for the suspect haplotype. Sequence data from Southwind (an affected calf) and the sire of the affected calf were scanned for candidate mutations. A missense mutation with a deleterious projected impact at 79,613,592 bp was homozygous in the affected calf and heterozygous in the calf's sire and Southwind. Sequence data available from the Cooperative Dairy DNA Repository for 299 other Holsteins indicated a 97% concordance with the haplotype and an 89% call rate. The exon amino acid sequence appears to be broadly conserved in the CACNA1S gene, and mutations in humans and mice can cause phenotypes of temporary or permanent paralysis analogous to those in calves with the haplotype causing muscle weakness (HMW). Improved methods for using pedigree to track new mutations within existing haplotypes were developed and applied to the haplotypes for both muscle weakness and Holstein cholesterol deficiency (HCD). For HCD, concordance of the gene test with its haplotype status was greatly improved. For both defects, haplotype status was matched to heifer livability records for 558,000 calves. For HMW, only 46 heifers with livability records were homozygous and traced only to Southwind on both sides. Of those, 52% died before 18 mo at an average age of 1.7 ± 1.6 mo, but that death rate may be underestimated if only healthier calves were genotyped. The death rate was 2.4% for noncarriers. Different reporting methods or dominance effects may be needed to include HMW and other partially lethal effects in selection and mating. Direct tests are needed for new mutations within existing common haplotypes because tracking can be difficult even with accurate pedigrees when the original haplotype has a high frequency.

Inheritance of a mutation causing neuropathy with splayed forelimbs in Jersey cattle

A new undesirable genetic factor, neuropathy with splayed forelimbs (JNS), has been identified recently in the Jersey breed. Calves affected with JNS are unable to stand on splayed forelimbs that exhibit significant extensor rigidity and excessive lateral abduction at birth. Affected calves generally are alert at birth but exhibit neurologic symptoms, including spasticity of head and neck and convulsive behavior. Other symptoms reported include dislocated shoulders, congenital craniofacial anomalies, and degenerative myelopathy. Inheritance of an undesirable genetic factor was determined from a study of 16 affected calves reported by Jersey breeders across the United States. All of their pedigrees traced back on both paternal and maternal sides to a common ancestor born in 1995. Genotypes revealed that JNS is attributable to a specific haplotype on Bos taurus autosome 6. Currently 8.2% of the genotyped US Jersey population are carriers of the haplotype. Sequencing of the region of shared homozygosity revealed missense variant rs1116058914 at base 60,158,901 of the ARS-UCD1.2 reference map as the most concordant with the genetic condition and the most likely cause. The single-base G to A substitution is in the coding region of the last exon of UCHL1, which is conserved across species. Mutations in humans and gene knockouts in mice cause similar recessive symptoms and muscular degeneration. Since December 2020, carrier status has been tracked using the identified haplotype and reported for all 459,784 genotyped Jersey animals. With random mating, about 2,200 affected calves per year with losses of about $250,000 would result from the 1.3 million US Jersey cows in the national population. Selection and mating programs can reduce numbers of JNS-affected births using either the haplotype status or a direct gene test in the future. Breeders should report calf abnormalities to their breed association to help discover new defects such as JNS.

Development of polymorphic markers in the immune gene complex loci of cattle

The addition of cattle health and immunity traits to genomic selection indices holds promise to increase individual animal longevity and productivity, and decrease economic losses from disease. However, highly variable genomic loci that contain multiple immune-related genes were poorly assembled in the first iterations of the cattle reference genome assembly and underrepresented during the development of most commercial genotyping platforms. As a consequence, there is a paucity of genetic markers within these loci that may track haplotypes related to disease susceptibility. By using hierarchical assembly of bacterial artificial chromosome inserts spanning 3 of these immune-related gene regions, we were able to assemble multiple full-length haplotypes of the major histocompatibility complex, the leukocyte receptor complex, and the natural killer cell complex. Using these new assemblies and the recently released ARS-UCD1.2 reference, we aligned whole-genome shotgun reads from 125 sequenced Holstein bulls to discover candidate variants for genetic marker development. We selected 124 SNPs, using heuristic and statistical models to develop a custom genotyping panel. In a proof-of-principle study, we used this custom panel to genotype 1,797 Holstein cows exposed to bovine tuberculosis (bTB) that were the subject of a previous GWAS study using the Illumina BovineHD array. Although we did not identify any significant association of bTB phenotypes with these new genetic markers, 2 markers exhibited substantial effects on bTB phenotypic prediction. The models and parameters trained in this study serve as a guide for future marker discovery surveys particularly in previously unassembled regions of the cattle genome.

Marker selection and genomic prediction of economically important traits using imputed high-density genotypes for 5 breeds of dairy cattle

Marker sets used in US dairy genomic predictions were previously expanded by including high-density (HD) or sequence markers with the largest effects for Holstein breed only. Other non-Holstein breeds lacked enough HD genotyped animals to be used as a reference population at that time, and thus were not included in the genomic prediction. Recently, numbers of non-Holstein breeds genotyped using HD panels reached an acceptable level for imputation and marker selection, allowing HD genomic prediction and HD marker selection for Holstein plus 4 other breeds. Genotypes for 351,461 Holsteins, 347,570 Jerseys, 42,346 Brown Swiss, 9,364 Ayrshires (including Red dairy cattle), and 4,599 Guernseys were imputed to the HD marker list that included 643,059 SNP. The separate HD reference populations included Illumina BovineHD (San Diego, CA) genotypes for 4,012 Holsteins, 407 Jerseys, 181 Brown Swiss, 527 Ayrshires, and 147 Guernseys. The 643,059 variants included the HD SNP and all 79,254 (80K) genetic markers and QTL used in routine national genomic evaluations. Before imputation, approximately 91 to 97% of genotypes were unknown for each breed; after imputation, 1.1% of Holstein, 3.2% of Jersey, 6.7% of Brown Swiss, 4.8% of Ayrshire, and 4.2% of Guernsey alleles remained unknown due to lower density haplotypes that had no matching HD haplotype. The higher remaining missing rates in non-Holstein breeds are mainly due to fewer HD genotyped animals in the imputation reference populations. Allele effects for up to 39 traits were estimated separately within each breed using phenotypic reference populations that included up to 6,157 Jersey males and 110,130 Jersey females. Correlations of HD with 80K genomic predictions for young animals averaged 0.986, 0.989, 0.985, 0.992, and 0.978 for Jersey, Ayrshire, Brown Swiss, Guernsey, and Holstein breeds, respectively. Correlations were highest for yield traits (about 0.991) and lowest for foot angle and rear legs-side view (0.981and 0.982, respectively). Some HD effects were more than twice as large as the largest 80K SNP effect, and HD markers had larger effects than nearby 80K markers for many breed-trait combinations. Previous studies selected and included markers with large effects for Holstein traits; the newly selected HD markers should also improve non-Holstein and crossbred genomic predictions and were added to official US genomic predictions in April 2020.

Major quantitative trait loci influencing milk production and conformation traits in Guernsey dairy cattle detected on Bos taurus autosome 19

The goal of this study was to identify potential quantitative trait loci (QTL) for 27 production, fitness, and conformation traits of Guernsey cattle through genome-wide association (GWA) analyses, with extra emphasis on BTA19, where major QTL were observed for several traits. Animals' de-regressed predicted transmitting abilities (PTA) from the December 2018 traditional US evaluation were used as phenotypes. All of the Guernsey cattle included in the QTL analyses were predictor animals in the reference population, ranging from 1,077 to 1,685 animals for different traits. Single-trait GWA analyses were carried out by a mixed-model approach for all 27 traits using imputed high-density genotypes. A major QTL was detected on BTA19, influencing several milk production traits, conformation traits, and livability of Guernsey cattle, and the most significant SNP lie in the region of 26.2 to 28.3 Mb. The myosin heavy chain 10 (MYH10) gene residing within this region was found to be highly associated with milk production and body conformation traits of dairy cattle. After the initial GWA analyses, which suggested that many significant SNP are in linkage with one another, conditional analyses were used for fine mapping. The top significant SNP on BTA19 were fixed as covariables in the model, one at a time, until no more significant SNP were detected on BTA19. After this fine-mapping approach was applied, only 1 significant SNP was detected on BTA19 for most traits, but multiple, independent significant SNP were found for protein yield, dairy form, and stature. In addition, the haplotype that hosts the major QTL on BTA19 was traced to a US Guernsey born in 1954. The haplotype is common in the breed, indicating a long-term influence of this QTL on the US Guernsey population.

High-density genome-wide association study for residual feed intake in Holstein dairy cattle

Improving feed efficiency (FE) of dairy cattle may boost farm profitability and reduce the environmental footprint of the dairy industry. Residual feed intake (RFI), a candidate FE trait in dairy cattle, can be defined to be genetically uncorrelated with major energy sink traits (e.g., milk production, body weight) by including genomic predicted transmitting ability of such traits in genetic analyses for RFI. We examined the genetic basis of RFI through genome-wide association (GWA) analyses and post-GWA enrichment analyses and identified candidate genes and biological pathways associated with RFI in dairy cattle. Data were collected from 4,823 lactations of 3,947 Holstein cows in 9 research herds in the United States. Of these cows, 3,555 were genotyped and were imputed to a high-density list of 312,614 SNP. We used a single-step GWA method to combine information from genotyped and nongenotyped animals with phenotypes as well as their ancestors' information. The estimated genomic breeding values from a single-step genomic BLUP were back-solved to obtain the individual SNP effects for RFI. The proportion of genetic variance explained by each 5-SNP sliding window was also calculated for RFI. Our GWA analyses suggested that RFI is a highly polygenic trait regulated by many genes with small effects. The closest genes to the top SNP and sliding windows were associated with dry matter intake (DMI), RFI, energy homeostasis and energy balance regulation, digestion and metabolism of carbohydrates and proteins, immune regulation, leptin signaling, mitochondrial ATP activities, rumen development, skeletal muscle development, and spermatogenesis. The region of 40.7 to 41.5 Mb on BTA25 (UMD3.1 reference genome) was the top associated region for RFI. The closest genes to this region, CARD11 and EIF3B, were previously shown to be related to RFI of dairy cattle and FE of broilers, respectively. Another candidate region, 57.7 to 58.2 Mb on BTA18, which is associated with DMI and leptin signaling, was also associated with RFI in this study. Post-GWA enrichment analyses used a sum-based marker-set test based on 4 public annotation databases: Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, Reactome pathways, and medical subject heading (MeSH) terms. Results of these analyses were consistent with those from the top GWA signals. Across the 4 databases, GWA signals for RFI were highly enriched in the biosynthesis and metabolism of amino acids and proteins, digestion and metabolism of carbohydrates, skeletal development, mitochondrial electron transport, immunity, rumen bacteria activities, and sperm motility. Our findings offer novel insight into the genetic basis of RFI and identify candidate regions and biological pathways associated with RFI in dairy cattle.

Short communication: Phenotypic and genetic effects of the polled haplotype on yield, longevity, and fertility in US Brown Swiss, Holstein, and Jersey cattle

Phenotypes from the December 2018 US national genetic evaluations were used to compute effects of the polled haplotype in US Brown Swiss (BS), Holstein (HO), and Jersey (JE) cattle on milk, fat, and protein yields, somatic cell score, single-trait productive life, daughter pregnancy rate, heifer conception rate, and cow conception rate. Lactation records pre-adjusted for nongenetic factors and direct genomic values were used to estimate phenotypic and genetic effects of the polled haplotype, respectively. No phenotypic or direct genomic values effects were different from zero for any trait in any breed. Genomic PTA (gPTA) for the lifetime net merit (NM$) selection index of bulls born since January 1, 2012, that received a marketing code from the National Association of Animal Breeders (Madison, WI), and cows born on or after January 1, 2015, were compared to determine whether there was a systematic benefit to polled or horned genetics. Horned bulls had the highest average gPTA for NM$ in all 3 breeds, but that difference was significant only in HO and JE (HO: 615.4 ± 1.9, JE: 402.3 ± 3.4). Homozygous polled BS cows had significantly higher average gPTA for NM$ than their heterozygous polled or horned contemporaries (PP = 261.4 ± 43.5, Pp = 166.1 ± 13.7, pp = 174.1 ± 1.8), but the sample size was very small (n = 9). In HO and JE, horned cows had higher gPTA for NM$ (HO = 378.3 ± 0.2, JE = 283.3 ± 0.3). Selection for polled cattle should not have a detrimental effect on yield, fertility, or longevity, but these differences show that, in the short term, selection for polled over horned cattle will result in lower rates of genetic gain.

Genetic and nongenetic profiling of milk pregnancy-associated glycoproteins in Holstein cattle

Pregnancy-associated glycoproteins (PAG) are secreted by the trophoblast and are detectable in maternal circulation around the time of attachment of the fetal placenta, as well as in blood and milk throughout gestation. The understanding of the genetic mechanisms controlling PAG levels can confer advantages for livestock breeding programs given the precocity and the ease of obtaining this phenotype from routine pregnancy diagnosis. The aim of this study was to characterize PAG levels by estimating genetic parameters and correlations with other dairy traits, and to identify genomic regions and candidate genes associated with PAG levels in milk. The PAG data consisted of pregnancy diagnoses using commercial assays from 2012 to 2017, and genotype data consisted of 54,123 SNP markers for 2,352 individuals (embryos and dams). The model included contemporary group (herd, year, and season) and embryo age as fixed effects, and random embryonic (direct) and maternal (indirect) genetic effects. Using genomic data, the estimated heritability for direct and maternal genetic effects (± standard deviations) were 0.23 ± 0.05 and 0.11 ± 0.05, respectively. The genetic correlation between these effects was almost zero (0.001 ± 0.06). A preliminary analysis revealed low correlations between milk PAG levels and other dairy traits. The genome-wide association analysis was performed using 2 approaches: single-marker and single-step using all markers. Four genomic regions with direct genetic effects were detected on Bos taurus autosome (BTA) 6, BTA7, BTA19, and BTA29 of the embryonic genome. The BTA29 locus was within the bovine PAG gene cluster, suggesting a cis-regulatory quantitative trait locus on the PAG expression. However, other associations, without an obvious link to PAG expression, could be related to the transportation of PAG and their concentration in milk. Only 1 region from the maternal genome, on BTA4, had a significant indirect effect, where WNT2 is a candidate gene related to placenta vascularization and gestation health. Collectively, our results suggest a moderate genetic control of milk PAG levels from both maternal and fetal genomes, but larger studies are needed to fully evaluate the usefulness of milk PAG in the genetic evaluation of fetal growth and cow fertility.

Explorations in genome-wide association studies and network analyses with dairy cattle fertility traits

The objective of this study was to identify single nucleotide polymorphisms and gene networks associated with 3 fertility traits in dairy cattle-daughter pregnancy rate, heifer conception rate, and cow conception rate-using different approaches. Deregressed predicted transmitting abilities were available for approximately 24,000 Holstein bulls and 36,000 Holstein cows sampled from the National Dairy Database with high-density genotypes. Of those, 1,732 bulls and 375 cows had been genotyped with the Illumina BovineHD Genotyping BeadChip (Illumina Inc., San Diego, CA). The remaining animals were genotyped with various chips of lower density that were imputed to high density. Univariate and trivariate genome-wide association studies (GWAS) with both medium- (60,671 markers) and high-density (312,614 markers) panels were performed for daughter pregnancy rate, heifer conception rate, and cow conception rate using GEMMA (version 0.94; http://www.xzlab.org/software.html). Analyses were conducted using bulls only, cows only, and a sample of both bulls and cows. The partial correlation and information theory algorithm was used to develop gene interaction networks. The most significant markers were further investigated to identify putatively associated genes. Little overlap in associated genes could be found between GWAS using different reference populations of bulls only, cows only, and combined bulls and cows. The partial correlation and information theory algorithm was able to identify several genes that were not identified by ordinary GWAS. The results obtained herein will aid in further dissecting the complex biology underlying fertility traits in dairy cattle, while also providing insight into the nuances of GWAS.

Reducing animal sequencing redundancy by preferentially selecting animals with low-frequency haplotypes

Many studies leverage targeted whole-genome sequencing (WGS) experiments to identify rare and causal variants within populations. As a natural consequence of their experimental design, many of these surveys tend to sequence redundant haplotype segments due to their high frequency in the base population, and the variants discovered within sequencing data are difficult to phase. We propose a new algorithm, called inverse weight selection (IWS), that preferentially selects individuals based on the cumulative presence of rare frequency haplotypes to maximize the efficiency of WGS surveys. To test the efficacy of this method, we used genotype data from 112,113 registered Holstein bulls derived from the US national dairy database. We demonstrate that IWS is at least 6.8% more efficient than previously published methods in selecting the least number of individuals required to sequence all haplotype segments ≥4% frequency in the US Holstein population. We also suggest that future surveys focus on sequencing homozygous haplotype segments as a first pass to achieve a 50% reduction in cost with an added benefit of phasing variant calls efficiently. Together, this new selection algorithm and experimental design suggestion significantly reduce the overall cost of variant discovery through WGS experiments, making surveys for causal variants influencing disease and production even more efficient.

256 GENOMIC EVALUATION OF FERTILITY TRAITS AND DISCOVERY OF HAPLOTYPES THAT AFFECT FERTILITY OF US DAIRY CATTLE

Genomic evaluations of dairy cattle became official in the United States in January 2009 for Holsteins and Jerseys, and later for Brown Swiss, Ayrshires, and Guernseys. Up to 33 yield, fitness, calving, and conformation traits are evaluated, and the fertility traits included daughter pregnancy rate and heifer and cow conception rates. Additional fertility traits, such as age at first calving and days from calving to first insemination, also are being studied. Male fertility (sire conception rate) is evaluated phenotypically rather than through genomics. Over 1 million animals have genotypes in the national database, which reflects collaboration with Canada and Europe. Most of the genotypes are from females and are from genotyping chips with <30 000 single nucleotide polymorphisms (SNP). To combine data across chips, genotypes are imputed to a set of >77 000 SNP. The imputation process involves dividing the chromosome into segments of approximately equal length and determining the paternal or maternal origin of the alleles. Because some segments were never homozygous, they were assumed to contain an abnormality that resulted in early embryonic death. If a decrease in sire conception rate could be associated with a bull that was a carrier of such a chromosomal segment, the haplotype was designated as affecting fertility. Once the region was identified, bioinformatic analysis was used to discover the causative variant for many of those haplotypes. Accuracy of genomic evaluations is determined by size of the reference population and heritability of the trait. The reference population for Holsteins includes >180 000 bulls and cows. Because fertility traits have low heritabilities, genomic information is particularly useful in improving evaluation accuracy. Accuracy of fertility evaluations is expected to increase further by discovering causative variants for various aspects of conception and gestation through investigation of sequence data.

5 A SINGLE NUCLEOTIDE POLYMORPHISM IN COQ9 AFFECTS MITOCHONDRIAL FUNCTION, BODY WEIGHT CHANGE AFTER CALVING, AND FERTILITY IN HOLSTEIN COWS

A single nucleotide polymorphism for COQ9 has been associated with genetic merit for fertility in 2 separate populations of Holstein cattle, with the A allele associated with higher fertility. COQ9 is necessary for the synthesis of coenzyme Q10, a component of the electron transport system of the mitochondria. We evaluated the effect of COQ9 genotype on the electron transport system, body weight changes after calving, and phenotypic measurements of fertility and production in Holstein cows. The single nucleotide polymorphism in COQ9 was genotyped using a Sequenom MassARRAY® (Sequenom Inc., San Diego, CA, USA). In the first study, cows ≥200 days in milk were selected for analysis of mitochondrial oxygen consumption [COQ9 genotype: AA (n = 12), AG (n = 12), and GG (n = 12)]. Peripheral blood mononuclear cells were isolated and respiration assessed using the Oroboros O2k high-resolution respirometer to evaluate routine respiration, R; leak respiration, L; and electron transport system capacity, E. There were additive effects of genotype on respiratory function (P < 0.05): R was 3.4 ± 0.3, 4.7 ± 0.3, and 4.9 ± 0.3 pmol of O2/s per 106 cells, L was 1.9 ± 0.3, 2.7 ± 0.3, and 3.0 ± 0.3 pmol of O2/s per 106 cells, and the uncoupling control ratio (E/R) was 3.4 ± 0.2, 2.5 ± 0.2, and 2.1 ± 0.2 for AA, AG, and GG, respectively. In a second study, body weight was recorded for AA (n = 106), AG (n = 223), and GG (n = 86) cows during the first 20 weeks postpartum for 2 consecutive lactations. In both lactations, body weight postpartum was affected by genotype × time postpartum (P < 0.001), with cows of the AA genotype experiencing less weight loss than AG (second lactation only) and GG cows. Days open, services per conception, and 305-day milk yield (MY) for the first 2 lactations were evaluated in a population of 2273 Holstein cows grouped based on predicted transmitting ability for daughter pregnancy rate: ≤–1 (n = 1220) and ≥1.5 (n = 1053). Continuous data were analysed using the MIXED procedure of SAS, and categorical data were analysed using the GLIMMIX procedure. The model included farm, genotype, and the numerator relationship matrix to account for (co)variances among animals. Additive and dominance effects were estimated. Genotype affected each trait (P < 0.05). Values for AA, AG, and GG for the first lactation were as follows: days open, 123.6 ± 3.5, 134.3 ± 2.8, and 139.4 ± 3.5 days; services per conception, 2.4 ± 0.1, 2.5 ± 0.1, and 2.7 ± 0.1; and MY, 11 278 ± 65, 11 416 ± 51, and 11 478 ± 65 kg. For the second lactation COQ9 affected (P < 0.05) days open (133.2 ± 4.7, 142.9 ± 3.1, and 147.9 ± 3.9 days) and services per conception (2.5 ± 0.1, 2.6 ± 0.8, and 2.7 ± 0.1), but there was no effect (P = 0.63) on MY (11 486 ± 66, 11 502 ± 52, and 11 526 ± 57 kg). Results indicate that the same genotype associated with genetic merit for fertility (AA) is associated with more efficient respiratory function and less body-weight loss postpartum. Moreover, the favourable genotype was associated with higher phenotypic measurements of fertility and lower MY. Results indicate the single nucleotide polymorphism in COQ9 could be a potential marker for fertility and that allelic variants may affect fertility by altering respiratory efficiency. Study was supported by USDA AFRI 2013–68004–20365.

Single nucleotide polymorphisms associated with thermoregulation in lactating dairy cows exposed to heat stress

Dairy cows with increased rectal temperature experience lower milk yield and fertility. Rectal temperature during heat stress is heritable, so genetic selection for body temperature regulation could reduce effects of heat stress on production. One aim of the study was to validate the relationship between genotype and heat tolerance for single nucleotide polymorphisms (SNPs) previously associated with resistance to heat stress. A second aim was to identify new SNPs associated with heat stress resistance. Thermotolerance was assessed in lactating Holsteins during the summer by measuring rectal temperature (a direct measurement of body temperature regulation; n = 435), respiration rate (an indirect measurement of body temperature regulation, n = 450) and sweating rate (the major evaporative cooling mechanism in cattle, n = 455). The association between genotype and thermotolerance was evaluated for 19 SNPs previously associated with rectal temperature from a genomewide analysis study (GWAS), four SNPs previously associated with change in milk yield during heat stress from GWAS, 2 candidate gene SNPs previously associated with rectal temperature and respiration rate during heat stress (ATPA1A and HSP70A) and 66 SNPs in genes previously shown to be associated with reproduction, production or health traits in Holsteins. For SNPs previously associated with heat tolerance, regions of BTA4, BTA6 and BTA24 were associated with rectal temperature; regions of BTA6 and BTA24 were associated with respiration rate; and regions of BTA5, BTA26 and BTA29 were associated with sweating rate. New SNPs were identified for rectal temperature (n = 12), respiration rate (n = 8) and sweating rate (n = 3) from among those previously associated with production, reproduction or health traits. The SNP that explained the most variation were PGR and ASL for rectal temperature, ACAT2 and HSD17B7 for respiration rate, and ARL6IP1 and SERPINE2 for sweating rate. ARL6IP1 was associated with all three thermotolerance traits. In conclusion, specific genetic markers responsible for genetic variation in thermoregulation during heat stress in Holsteins were identified. These markers may prove useful in genetic selection for heat tolerance in Holstein cattle.

Short communication: Projecting milk yield using best prediction and the MilkBot lactation model

The accuracy and precision of 3 lactation models was estimated by summarizing means and variability in projection error for next-test milk and actual 305-d milk yield (M305) for 50-d intervals in a large Dairy Herd Improvement Association data set. Lactations were grouped by breed (Holstein, Jersey, and crossbred) and parity (first vs. later). A smaller, single-herd data set with both Dairy Herd Improvement Association data and daily milk weights was used to compare M305 calculated from test-day data with M305 computed by summing daily milk weights. The lactation models tested were best prediction (BP), the nonlinear MilkBot (MB) model, and a null model (NM) based on a stepwise function. The accuracy of the models was ranked (best to worst) MB, BP, and NM for later-parity cows and MB, NM, and BP for first-parity cows, with MB achieving accuracy in projecting daily milk of 0.5 kg or better in most groups. The models generally showed better accuracy after 50 d in milk. Best prediction and NM had low accuracy for crossbred cows and first-parity Holstein and Jersey cows. The MB model appears to be more precise than BP, and NM had low precision, especially for M305. Regression of model-generated M305 on summed M305 showed BP and MB to be equally efficient in ranking lactations, but MB was better at quantifying differences.

Genomic imputation and evaluation using high-density Holstein genotypes

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.

Heritability of rectal temperature and genetic correlations with production and reproduction traits in dairy cattle

Genetic selection for body temperature during heat stress might be a useful approach to reduce the magnitude of heat stress effects on production and reproduction. Objectives of the study were to estimate the genetic parameters of rectal temperature (RT) in dairy cows in freestall barns under heat stress conditions and to determine the genetic and phenotypic correlations of rectal temperature with other traits. Afternoon RT were measured in a total of 1,695 lactating Holstein cows sired by 509 bulls during the summer in North Florida. Genetic parameters were estimated with Gibbs sampling, and best linear unbiased predictions of breeding values were predicted using an animal model. The heritability of RT was estimated to be 0.17 ± 0.13. Predicted transmitting abilities for rectal temperature changed 0.0068 ± 0.0020°C/yr from (birth year) 2002 to 2008. Approximate genetic correlations between RT and 305-d milk, fat, and protein yields, productive life, and net merit were significant and positive, whereas approximate genetic correlations between RT and somatic cell count score and daughter pregnancy rate were significant and negative. Rectal temperature during heat stress has moderate heritability, but genetic correlations with economically important traits mean that selection for RT could lead to lower productivity unless methods are used to identify genes affecting RT that do not adversely affect other traits of economic importance.

Short communication: best prediction of 305-day lactation yields with regional and seasonal effects

In the United States, lactation yields are calculated using best prediction (BP), a method in which test-day (TD) data are compared with breed- and parity-specific herd lactation curves that do not account for differences among regions of the country or seasons of calving. Complete data from 538,090 lactations of 348,123 Holstein cows with lactation lengths between 250 and 500 d, records made in a single herd, at least 5 reported TD, and twice-daily milking were extracted from the national dairy database and used to construct regional and seasonal lactation curves. Herds were assigned to 1 of 7 regions of the country, individual lactations were assigned to 3-mo seasons of calving, and lactation curves for milk, fat, and protein yields were estimated by parity group for regions, seasons, and seasons within regions. Multiplicative pre-adjustment factors (MF) also were computed. The resulting lactation curves and MF were tested on a validation data set of 891,806 lactations from 400,000 Holstein cows sampled at random from the national dairy database. Mature-equivalent milk, fat, and protein yields were calculated using the standard and adjusted curves and MF, and differences between 305-d mature-equivalent yields were tested for significance. Yields calculated using 50-d intervals from 50 to 250 d in milk (DIM) and using all TD to 500 DIM allowed comparisons of predictions for records in progress (RIP). Differences in mature-equivalent milk ranged from 0 to 51 kg and were slightly larger for first-parity than for later parity cows. Milk and components yields did not differ significantly in any case. Correlations of yields for 50-d intervals with those using all TD were similar across analyses. Yields for RIP were slightly more accurate when adjusted for regional and seasonal differences.

Genetic evaluation of lactation persistency for five breeds of dairy cattle

Cows with high lactation persistency tend to produce less milk than expected at the beginning of lactation and more than expected at the end. Best prediction of lactation persistency is calculated as a function of trait-specific standard lactation curves and linear regressions of test-day deviations on days in milk. Because regression coefficients are deviations from a tipping point selected to make yield and lactation persistency phenotypically uncorrelated it should be possible to use 305-d actual yield and lactation persistency to predict yield for lactations with later endpoints. The objectives of this study were to calculate (co)variance components and breeding values for best predictions of lactation persistency of milk (PM), fat (PF), protein (PP), and somatic cell score (PSCS) in breeds other than Holstein, and to demonstrate the calculation of prediction equations for 400-d actual milk yield. Data included lactations from Ayrshire, Brown Swiss, Guernsey (GU), Jersey (JE), and Milking Shorthorn (MS) cows calving since 1997. The number of sires evaluated ranged from 86 (MS) to 3,192 (JE), and mean sire estimated breeding value for PM ranged from 0.001 (Ayrshire) to 0.10 (Brown Swiss); mean estimated breeding value for PSCS ranged from -0.01 (MS) to -0.043 (JE). Heritabilities were generally highest for PM (0.09 to 0.15) and lowest for PSCS (0.03 to 0.06), with PF and PP having intermediate values (0.07 to 0.13). Repeatabilities varied considerably between breeds, ranging from 0.08 (PSCS in GU, JE, and MS) to 0.28 (PM in GU). Genetic correlations of PM, PF, and PP with PSCS were moderate and favorable (negative), indicating that increasing lactation persistency of yield traits is associated with decreases in lactation persistency of SCS, as expected. Genetic correlations among yield and lactation persistency were low to moderate and ranged from -0.55 (PP in GU) to 0.40 (PP in MS). Prediction equations for 400-d milk yield were calculated for each breed by regression of both 305-d yield and 305-d yield and lactation persistency on 400-d yield. Goodness-of-fit was very good for both models, but the addition of lactation persistency to the model significantly improved fit in all cases. Routine genetic evaluations for lactation persistency, as well as the development of prediction equations for several lactation end-points, may provide producers with tools to better manage their herds.