Single-step genomic predictions for crossbred Holstein and Jersey cattle in the United States.

In the dairy cattle sector, the number of crossbred genotypes increased in the last years, and therefore, the need for accurate genomic evaluations for crossbred animals has also increased. Thus, this study aimed to investigate the feasibility of including crossbred genotypes in multibreed, single-step genomic BLUP (ssGBLUP) evaluations. The Council of Dairy Cattle Breeding provided more than 47 million lactation records registered between 2000 and 2021 in purebred Holstein and Jersey and their crosses. A total of 27 million animals were included in the analysis, of which 1.4 million were genotyped. Milk, fat, and protein yields were analyzed in a 3-trait repeatability model using BLUP or ssGBLUP. The 2 models were validated using prediction bias and accuracy computed for genotyped cows with no records in the truncated dataset and at least one lactation in the complete dataset. Bias and accuracy were better in the genomic model than in the pedigree-based one, with accuracies for crossbred cows higher than those of purebreds, except for fat yield in Holstein. Our study shows that genotypes for crossbred animals can be included in a ssGBLUP analysis with their purebred ancestors to estimate genomic estimated breeding values in a single run.

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.

Genomic evaluation of late-term abortion in cows recorded through Dairy Herd Improvement test plans.

Late-term abortions cause significant economic loss and are of great concern for dairy herds. Late-term abortions >152 d and <251 d of gestation that terminate a lactation or initiate a new lactation have long been recorded by Dairy Herd Improvement (DHI). For 24.8 million DHI lactations, the average recorded incidence of late-term abortions across all years (2001-2018) was 1.2%. However, the 1.3% incidence of abortions reported in 2012 has declined to <1.0% incidence since 2015. Small adjustments were applied to the 82 million daughter pregnancy rate (DPR), 29 million cow conception rate (CCR), and 9 million heifer conception rate (HCR) records to account for late-term abortions more accurately. Fertility credits for CCR and HCR were changed to treat the last breeding as a failure instead of success if the next calving was coded as a late-term abortion. Similarly, when computing DPR, days open is now set to the maximum value of 250 instead of the reported days open if the next reported calving is an abortion. The test of these changes showed very small changes in standard deviation and high correlations (0.997) of adjusted predicted transmitting abilities (PTA) with official PTA from about 20,000 Holstein bulls born since 2000 with >50% reliability. For late-term abortion as a trait, estimated heritability was only 0.001 and PTA had a standard deviation of only 0.1% for recent sires with high reliability (>75%). Young animal genomic PTA have near 50% reliability but range only from -0.5 to +0.4 because of the low incidence and heritability. Genetic trend was slightly favorable and late-term abortion PTA were correlated favorably by 0.27 with net merit, 0.49 with productive life, 0.33 with livability, 0.23 with CCR, 0.20 with HCR, 0.26 with DPR, -0.31 with somatic cell score, -0.24 with daughter stillbirth, and -0.26 with daughter dystocia. Thus, PTA for late-term abortions should not be needed as a separate fertility trait and instead these minor edit changes should suffice. The PTA for late-term abortions would add little value because national evaluations for current fertility traits already account for those economic losses.

Multibreed genomic evaluation for production traits of dairy cattle in the United States using single-step genomic best linear unbiased predictor.

Official multibreed genomic evaluations for dairy cattle in the United States are based on multibreed BLUP evaluation followed by single-breed estimation of SNP effects. Single-step genomic BLUP (ssGBLUP) allows the straight computation of genomic (G)EBV in a multibreed context. This work aimed to develop ssGBLUP multibreed genomic predictions for US dairy cattle using the algorithm for proven and young (APY) to compute the inverse of the genomic relationship matrix. Only purebred Ayrshire (AY), Brown Swiss (BS), Guernsey (GU), Holstein (HO), and Jersey (JE) animals were considered. A 3-trait model with milk (MY), fat (FY), and protein (PY) yields was applied using about 45 million phenotypes recorded from January 2000 to June 2020. The whole data set included about 29.5 million animals, of which almost 4 million were genotyped. All the effects in the model were breed specific, and breed was also considered as fixed unknown parent groups. Evaluations were done for (1) each single breed separately (single); (2) HO and JE together (HO_JE); (3) AY, BS, and GU together (AY_BS_GU); (4) all the 5 breeds together (5_BREEDS). Initially, 15k core animals were used in APY for AY_BS_GU and 5_BREEDS, but larger core sets with more animals from the least represented breeds were also tested. The HO_JE evaluation had a fixed set of 30k core animals, with an equal representation of the 2 breeds, whereas HO and JE single-breed analysis involved 15k core animals. Validation for cows was based on correlations between adjusted phenotypes and (G)EBV, whereas for bulls on the regression of daughter yield deviations on (G)EBV. Because breed was correctly considered in the model, BLUP results for single and multibreed analyses were the same. Under ssGBLUP, predictability and reliability for AY, BS, and GU were on average 7% and 2% lower in 5_BREEDS compared with single-breed evaluations, respectively. However, validation parameters for these 3 breeds became better than in the single-breed evaluations when 45k animals were included in the core set for 5_BREEDS. Evaluations for Holsteins were more stable across scenarios because of the greatest number of genotyped animals and amount of data. Combining AY, BS, and GU into one evaluation resulted in predictions similar to the ones from single breed, especially when using about 30k core animals in APY. The results showed that single-step large-scale multibreed evaluations are computationally feasible, but fine tuning is needed to avoid a reduction in reliability when numerically dominant breeds are combined. Having evaluations for AY, BS, and GU separated from HO and JE may reduce inflation of GEBV for the first 3 breeds.

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.

Genomic evaluation of dairy heifer livability.

Differences in breeds and sire lines suggest the presence of a genetic component for heifer livability (HLIV). Genomic evaluation for this trait can increase profitability and improve animal health and welfare. Evaluations for HLIV were examined from 3,362,499 calf data records from heifers of all breeds born from 2009 to 2016. Data were obtained from the national cooperator database maintained by the Council on Dairy Cattle Breeding (https://www.uscdcb.com/). The total number of deaths reported was 134,753 (4.01%), which included herds with death loss between 1.5 and 25.5%. Age at death was evaluated and ranged from >2 d of age until the heifer left the herd, with a maximum of 18 mo of age. Records were not included until 3 yr after the birthdate so that live status of contemporaries could be confirmed by a calving date for those animals. Deaths observed until 2 d after birth were considered to be a stillbirth rather than a failure of HLIV. The scale used for analysis of HLIV was 0 (died) or 100 (live), and the heritability estimate was 0.7% based on sire model with restricted maximum likelihood estimation. Genomic predicted transmitting abilities for Holstein ranged from -1.6% to +1.6% with a standard deviation of 0.5%, and genomic predicted transmitting abilities for Jersey ranged from -0.5% to +0.5% with a standard deviation of 0.2%. The mean overall death loss was about 4%. Reliabilities of genomic predictions for young animals averaged 46% for Holsteins and 30% for Jerseys, and corresponding traditional parent average reliabilities averaged 16% and 12%, respectively. Correlations of HLIV were 0.44 with productive life, 0.18 to 0.22 with yield traits, and 0.29 with early first calving on proven Holstein bulls. The HLIV trait had a favorable genetic trend in recent years, likely because of the indirect selection associated with the correlated traits. The trait HLIV should receive 1% of emphasis on the Lifetime Net Merit index, resulting in economic progress worth $50,000/yr. By encouraging more comprehensive recording on calf mortality, the reliabilities of genetic predictions could increase significantly.

Genomic predictions for yield traits in US Holsteins with unknown parent groups.

The objective of this study was to assess the reliability and bias of estimated breeding values (EBV) from traditional BLUP with unknown parent groups (UPG), genomic EBV (GEBV) from single-step genomic BLUP (ssGBLUP) with UPG for the pedigree relationship matrix (A) only (SS_UPG), and GEBV from ssGBLUP with UPG for both A and the relationship matrix among genotyped animals (A22; SS_UPG2) using 6 large phenotype-pedigree truncated Holstein data sets. The complete data included 80 million records for milk, fat, and protein yields from 31 million cows recorded since 1980. Phenotype-pedigree truncation scenarios included truncation of phenotypes for cows recorded before 1990 and 2000 combined with truncation of pedigree information after 2 or 3 ancestral generations. A total of 861,525 genotyped bulls with progeny and cows with phenotypic records were used in the analyses. Reliability and bias (inflation/deflation) of GEBV were obtained for 2,710 bulls based on deregressed proofs, and on 381,779 cows born after 2014 based on predictivity (adjusted cow phenotypes). The BLUP reliabilities for young bulls varied from 0.29 to 0.30 across traits and were unaffected by data truncation and number of generations in the pedigree. Reliabilities ranged from 0.54 to 0.69 for SS_UPG and were slightly affected by phenotype-pedigree truncation. Reliabilities ranged from 0.69 to 0.73 for SS_UPG2 and were unaffected by phenotype-pedigree truncation. The regression coefficient of bull deregressed proofs on (G)EBV (i.e., GEBV and EBV) ranged from 0.86 to 0.90 for BLUP, from 0.77 to 0.94 for SS_UPG, and was 1.00 ± 0.03 for SS_UPG2. Cow predictivity ranged from 0.22 to 0.28 for BLUP, 0.48 to 0.51 for SS_UPG, and 0.51 to 0.54 for SS_UPG2. The highest cow predictivities for BLUP were obtained with the most extreme truncation, whereas for SS_UPG2, cow predictivities were also unaffected by phenotype-pedigree truncations. The regression coefficient of cow predictivities on (G)EBV was 1.02 ± 0.02 for SS_UPG2 with the most extreme truncation, which indicated the least biased predictions. Computations with the complete data set took 17 h with BLUP, 58 h with SS_UPG, and 23 h with SS_UPG2. The same computations with the most extreme phenotype-pedigree truncation took 7, 36, and 15 h, respectively. The SS_UPG2 converged in fewer rounds than BLUP, whereas SS_UPG took up to twice as many rounds. Thus, the ssGBLUP with UPG assigned to both A and A22 provided accurate and unbiased evaluations, regardless of phenotype-pedigree truncation scenario. Old phenotypes (before 2000 in this data set) did not affect the reliability of predictions for young selection candidates, especially in SS_UPG2.

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.

Investigating conception rate for beef service sires bred to dairy cows and heifers.

The widespread use of sexed semen on US dairy cows and heifers has led to an excess of replacement heifers' calves, and the sale prices for those calves are much lower than in the past. Females not selected to produce the next generation of replacement heifers are increasingly being bred to beef bulls to produce crossbred calves for beef production. The purpose of this study was to investigate the use of beef service sires bred to dairy cows and heifers and to provide a tool for dairy producers to evaluate beef service sires' conception. Sire conception rate (SCR) is a phenotypic evaluation of service sire fertility that is routinely calculated for US dairy bulls. A total of 268,174 breedings were available, which included 36 recognized beef breeds and 7 dairy breeds. Most of the beef-on-dairy inseminations (95.4%) were to Angus (AN) bulls. Because of the limited number of records among other breeds, we restricted our final evaluations to AN service sires bred to Holstein (HO) cows. Service-sire inbreeding and expected inbreeding of resulting embryo were set to zero because pedigree data for AN bulls were unavailable. There were 233,379 breedings from 1,344 AN service sire to 163,919 HO cows. A mean (SD) conception rate of 33.8% (47.3%) was observed compared with 34.3% (47.5%) for breedings with HO sires mated to HO cows. Publishable AN bulls were required to have ≥100 total matings, ≥10 matings in the most recent 12 mo, and breedings in at least 5 herds. Mean SCR reliability was 64.5% for 116 publishable bulls, with a maximum reliability of 99% based on 25,217 breedings. Average SCR was near zero (on AN base) with a range of -5.1 to 4.4. Breedings to HO heifers were also examined, which included 19,437 breedings (443 AN service sire and 15,971 HO heifers). A mean (SD) conception rate of 53.0% (49.9%) was observed, compared with 55.3% (49.7%) for breedings with a HO sire mated to a HO heifer. Beef sires were used more frequently in cows known to be problem breeders, which explains some of the difference in conception rate. Mean service number was 1.92 and 2.87 for HO heifers and 2.13 and 3.04 for HO cows mated to HO and AN sires, respectively. Mating dairy cows and heifers to beef bulls may be profitable if calf prices are higher, fertility is improved, or if practices such as sexed semen, genomic testing, and improved cow productive life allow herd owners to produce both higher quality dairy replacement and increased income from market calves.

Symposium review: How to implement genomic selection.

Genomic selection was adopted very quickly in the 10 yr after first implementation, and breeders continue to find new uses for genomic testing. Breeding values with higher reliability earlier in life are estimated by combining DNA genotypes for many thousands of loci using existing identification, pedigree, and phenotype databases for millions of animals. Quality control for both new and previous data is greatly improved by comparing genomic and pedigree relationships to correct parent-progeny conflicts and discover many additional ancestors. Many quantitative trait loci and gene tests have been added to previous assays that used only evenly spaced, highly polymorphic markers. Imputation now combines genotypes from many assays of differing marker densities. Prediction models have gradually advanced from normal or Bayesian distributions within trait and breed to single-step, multitrait, or other more complex models, such as multibreed models that may be needed for crossbred prediction. Genomic selection was initially applied to males to predict progeny performance but is now widely applied to females or even embryos to predict their own later performance. The initial focus on additive merit has expanded to include mating programs, genomic inbreeding, and recessive alleles. Many producers now use DNA testing to decide which heifers should be inseminated with elite dairy, beef, or sex-sorted semen, which should be embryo donors or recipients, or which should be sold or kept for breeding. Because some of these decisions are expensive to delay, predictions are now provided weekly instead of every few months. Predictions from international genomic databases are often more accurate and cost-effective than those from within-country databases that were previously designed for progeny testing unless local breeds, conditions, or traits differ greatly from the larger database. Selection indexes include many new traits, often with lower heritability or requiring large initial investments to obtain phenotypes, which provide further incentive to cooperate internationally. The genomic prediction methods developed for dairy cattle are now applied widely to many animal, human, and plant populations and could be applied to many more.

Symposium review: Development, implementation, and perspectives of health evaluations in the United States.

The rate at which new traits are being developed is increasing, leading to an expanding number of evaluations provided to dairy producers, especially for functional traits. This review will discuss the development and implementation of genetic evaluations for direct health traits in the United States, as well as potential future developments. Beginning in April 2018, routine official genomic evaluations for 6 direct health traits in Holsteins were made available to US producers from the Council on Dairy Cattle Breeding (Bowie, MD). Traits include resistance to milk fever, displaced abomasum, ketosis, clinical mastitis, metritis, and retained placenta. These health traits were included in net merit indices beginning in August 2018, with a total weight of approximately 2%. Previously, improvement of cow health was primarily made through changes to management practices or genetic selection on indicator traits, such as somatic cell score, productive life, or livability. Widespread genomic testing now allows for accelerated improvement of traits with low heritabilities such as health; however, phenotypes remain essential to the success of genomic evaluations. Establishment and maintenance of data pipelines is a critical component of health trait evaluations, as well as appropriate data quality control standards. Data standardization is a necessary process when multiple data sources are involved. Model refinement continues, including implementation of variance adjustments beginning with the April 2019 evaluation. Mastitis evaluations are submitted to Interbull along with somatic cell score for international validation and evaluation of udder health. Additional areas of research include evaluation of other breeds for direct health traits, use of multiple-trait models, and evaluations for additional functional traits such as calf health and feed efficiency. Future developments will require new and continued cooperation among numerous industry stakeholders. There is more information available than ever before with which to make better selection decisions; however, this also makes it increasingly important to provide accurate and unbiased information.

Genomic prediction of residual feed intake in US Holstein dairy cattle.

Genomic selection is an important tool to introduce feed efficiency into dairy cattle breeding. The goals of the current research are to estimate genomic breeding values of residual feed intake (RFI) and to assess the prediction reliability for RFI in the US Holstein population. The RFI data were collected from 4,823 lactations of 3,947 Holstein cows in 9 research herds in the United States, and were pre-adjusted to remove phenotypic correlations with milk energy, metabolic body weight, body weight change, and for several environmental effects. In the current analyses, genomic predicted transmitting abilities of milk energy and of body weight composite were included into the RFI model to further remove the genetic correlations that remained between RFI and these energy sinks. In the first part of the analyses, a national genomic evaluation for RFI was conducted for all the Holsteins in the national database using a standard multi-step genomic evaluation method and 60,671 SNP list. In the second part of the study, a single-step genomic prediction method was applied to estimate genomic breeding values of RFI for all cows with phenotypes, 5,252 elite young bulls, 4,029 young heifers, as well as their ancestors in the pedigree, using a high-density genotype chip. Theoretical prediction reliabilities were calculated for all the studied animals in the single-step genomic prediction by direct inversion of the mixed model equations. In the results, breeding values were estimated for 1.6 million genotyped Holsteins and 60 million ungenotyped Holsteins, The genomic predicted transmitting ability correlations between RFI and other traits in the index (e.g., fertility) are generally low, indicating minor correlated responses on other index traits when selecting for RFI. Genomic prediction reliabilities for RFI averaged 34% for all phenotyped animals and 13% for all 1.6 million genotyped animals. Including genomic information increased the prediction reliabilities for RFI compared with using only pedigree information. All bulls had low reliabilities, and averaged to only 16% for the top 100 net merit progeny-tested bulls. Analyses using single-step genomic prediction and high-density genotypes gave similar results to those obtained from the national evaluation. The average theoretical reliability for RFI was 18% among the elite young bulls under 5 yr old, being lower in the younger generations of elite bulls compared with older bulls. To conclude, the size of the reference population and its relationship to the predicted population remain as the limiting factors in the genomic prediction for RFI. Continued collection of feed intake data is necessary so that reliabilities can be maintained due to close relationships of phenotyped animals with breeding stock. Considering the currently low prediction reliability and high cost of data collection, focusing RFI data collection on relatives of elite bulls that will have the greatest genetic contribution to the next generation will give more gains and profit.

Discovering ancestors and connecting relatives in large genomic databases.

Genomic evaluation has improved both plant and animal breeding by allowing more accurate estimation of an individual's genetic potential. Because often only a small proportion of the population to be evaluated has been genotyped, genomic estimations rely heavily on complete pedigree information. Confirmation, discovery, and correction of parentage and connected relatives allow the creation of more complete pedigrees, which in turn increase the number of usable phenotypic records and prediction accuracy. Previous methods accounted for parent-progeny conflicts using SNP. More recently haplotype methods allowed discovery of distant relationships such as maternal grandsire (MGS) and maternal great-grandsire (MGGS) with improved accuracy. However, discovered MGS and MGGS often were not used, because no dam information was available to link them to the calf. An automated procedure to discover and fill missing maternal identification information was developed, allowing discovered MGS and MGGS to be used in imputation as well as in calculating breeding values for animals in the US dairy cattle database. An MGS was discovered for 295,136 animals with unknown dam, and the MGGS was discovered for 153,909 of these animals. A virtual maternal identification was added for animals with missing information. The effect of pedigree completion on progeny inbreeding, breeding values, and reliabilities was examined. Mean inbreeding of animals with missing maternal pedigree information was 6.69% before and 6.87% after pedigree assignment; expected future inbreeding was 7.24% before and 7.20% after assignment. Reliabilities for traditional breeding values increased from 26.6 to 32.6% for milk yield, 25.9 to 32.0% for fat yield, and 26.9 to 32.9% for protein yield; genomic reliabilities also increased slightly from 76.2 to 77.1% for milk, 76.0 to 76.9% for fat, and 76.3 to 77.3% for protein. The procedure developed for pedigree completion is a useful tool for improving accuracy of national and international evaluations and aiding producers in making better mating decisions.

Genomic predictions for crossbred dairy cattle.

Genomic evaluations are useful for crossbred as well as purebred populations when selection is applied to commercial herds. Dairy farmers had already spent more than $1 million to genotype over 32,000 crossbred animals before US genomic evaluations became available for those animals. Thus, new tools were needed to provide accurate genomic predictions for crossbreds. Genotypes for crossbreds are imputed more accurately when the imputation reference population includes purebreds. Therefore, genotypes of 6,296 crossbred animals were imputed from lower-density chips by including either 3,119 ancestors or 834,367 genotyped animals in the reference population. Crossbreds in the imputation study included 733 Jersey × Holstein F1 animals, 55 Brown Swiss × Holstein F1 animals, 2,300 Holstein backcrosses, 2,026 Jersey backcrosses, 27 Brown Swiss backcrosses, and 502 other crossbreds of various breed combinations. Another 653 animals appeared to be purebreds that owners had miscoded as a different breed. Genomic breed composition was estimated from 60,671 markers using the known breed identities for purebred, progeny-tested Holstein, Jersey, Brown Swiss, Ayrshire, and Guernsey bulls as the 5 traits (breed fractions) to be predicted. Estimates of breed composition were adjusted so that no percentages were negative or exceeded 100%, and breed percentages summed to 100%. Another adjustment set percentages above 93.5% equal to 100%, and the resulting value was termed breed base representation (BBR). Larger percentages of missing alleles were imputed by using a crossbred reference population rather than only the closest purebred reference population. Crossbred predictions were averages of genomic predictions computed using marker effects for each pure breed, which were weighted by the animal's BBR. Marker and polygenic effects were estimated separately for each breed on the all-breed scale instead of within-breed scales. For crossbreds, genomic predictions weighted by BBR were more accurate than the average of parents' breeding values and slightly more accurate than predictions using only the predominant breed. For purebreds, single-trait predictions using only within-breed data were as accurate as multi-trait predictions with allele effects in different breeds treated as correlated effects. Crossbred genomic predicted transmitting abilities were implemented by the Council on Dairy Cattle Breeding in April 2019 and will aid producers in managing their breeding programs and selecting replacement heifers.

Alternative SNP weighting for single-step genomic best linear unbiased predictor evaluation of stature in US Holsteins in the presence of selected sequence variants.

Causal variants inferred from sequence data analysis are expected to increase accuracy of genomic selection. In this work we evaluated the gain in reliability of genomic predictions, for stature in US Holsteins, when adding selected sequence variants to a pre-existent SNP chip. Two prediction methods were tested: de-regressed proofs assuming heterogeneous (genomic BLUP; GBLUP) residual variances and by single-step GBLUP (ssGBLUP) using actual phenotypes. Phenotypic data included 3,999,631 records for stature on 3,027,304 Holstein cows. Genotypes on 54,087 SNP markers (54k) were available for 26,877 bulls. Additionally, 16,648 selected sequence variants were combined with the 54k markers, for a total of 70,735 (70k) markers. In all methods, SNP in the genomic relationship matrix (G) were unweighted or weighted iteratively, with weights derived either by SNP effects squared or by a nonlinear method that resembles BayesA (nonlinear A). Reliability of genomic predictions were obtained by cross validation. With unweighted G derived from 54k markers, the reliabilities (× 100) were 72.4 for GBLUP and 75.3 for ssGBLUP. With unweighted G derived from 70k markers, the reliabilities were 73.4 and 76.0, respectively. Weighting by nonlinear A changed reliabilities to 73.3, and 75.9, respectively. Addition of selected sequence variants had a small effect on reliabilities. Weighting by quadratic functions reduced reliabilities. Weighting by nonlinear A increased reliabilities for GBLUP but had only a small effect in ssGBLUP. Reliabilities for direct genomic values extracted from ssGBLUP using unweighted G with 54k were higher than reliabilities by any GBLUP. Thus, ssGBLUP seems to capture more information than GBLUP and there is less room for extra reliability. Improvements in GBLUP may be because the weights in G change the covariance structure, which can explain a proportion of the variance that is accounted for when a heterogeneous residual variance is assumed by considering a different number of daughters per bull.

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.

Variance of gametic diversity and its application in selection programs.

The variance of gametic diversity ( σgamete2) can be used to find individuals that more likely produce progeny with extreme breeding values. The aim of this study was to obtain this variance for individuals from routine genomic evaluations, and to apply gametic variance in a selection criterion in conjunction with breeding values to improve genetic progress. An analytical approach was developed to estimate σgamete2 by the sum of binomial variances of all individual quantitative trait loci across the genome. Simulation was used to verify the predictability of this variance in a range of scenarios. The accuracy of prediction ranged from 0.49 to 0.85, depending on the scenario and model used. Compared with sequence data, SNP data are sufficient for estimating σgamete2 Results also suggested that markers with low minor allele frequency and the covariance between markers should be included in the estimation. To incorporate σgamete2 into selective breeding programs, we proposed a new index, relative predicted transmitting ability, which better utilizes the genetic potential of individuals than traditional predicted transmitting ability. Simulation with a small genome showed an additional genetic gain of up to 16% in 10 generations, depending on the number of quantitative trait loci and selection intensity. Finally, we applied σgamete2 to the US genomic evaluations for Holstein and Jersey cattle. As expected, the DGAT1 gene had a strong effect on the estimation of σgamete2 for several production traits. However, inbreeding had a small impact on gametic variability, with greater effect for more polygenic traits. In conclusion, gametic variance, a potentially important parameter for selection programs, can be easily computed and is useful for improving genetic progress and controlling genetic diversity.

Calling known variants and identifying new variants while rapidly aligning sequence data.

Whole-genome sequencing studies can identify causative mutations for subsequent use in genomic evaluations. Speed and accuracy of sequence alignment can be improved by accounting for known variant locations during alignment instead of calling the variants after alignment as in previous programs. The new programs Findmap and Findvar were compared with alignment using Burrows-Wheeler alignment (BWA) or SNAP and variant identification using Genome Analysis ToolKit (GATK) or SAMtools. Findmap stores the reference map and any known variant locations while aligning reads and counting reference and alternate alleles for each DNA source. Findmap also outputs potential new single nucleotide variant, insertion, and deletion alleles. Findvar separates likely true variants from read errors and outputs genotype probabilities. Strategies were tested using cattle, human, and a completely random reference map and simulated or actual data. Most tests simulated 10 bulls, each with 10× simulated sequence reads containing 39 million variants from the 1000 Bull Genomes Project. With 10 processors, clock times for processing 100× data were 105 h for BWA, 25 h for GATK, and 11 h for SAMtools but only about 4 h for SNAP, 3 h for Findmap, and 1 h for Findvar. Alignment programs required about the same total memory; BWA used 46 GB (4.6 GB/processor), whereas >10 processors can share the same memory in SNAP and Findmap, which used 40 and 46 GB, respectively. Findmap correctly mapped 92.9% of reads (compared with 92.6% from SNAP and 90.5% from BWA) and had high accuracy of calling alleles for known variants. For new variants, Findvar found 99.8% of single nucleotide variants, 79% of insertions, and 67% of deletions; GATK found 99.4, 95, and 90%, respectively; and SAMtools found 99.8, 12, and 16%, respectively. False positives (as percentages of true variants) were 11% of single nucleotide variants, 0.4% of insertions, and 0.3% of deletions from Findvar; 12, 8.4, and 2.9%, respectively, from GATK; and 37, 1.3, and 0.4%, respectively, from SAMtools. Advantages of Findmap and Findvar are fast processing, precise alignment, more useful data summaries, more compact output, and fewer steps. Calling known variants during alignment allows more efficient and accurate sequence-based genotyping.

Modeling missing pedigree in single-step genomic BLUP.

The objective was to compare methods of modeling missing pedigree in single-step genomic BLUP (ssGBLUP). Options for modeling missing pedigree included ignoring the missing pedigree, unknown parent groups (UPG) based on A (the numerator relationship matrix) or H (the unified pedigree and genomic relationship matrix), and metafounders. The assumptions for the distribution of estimated breeding values changed with the different models. We simulated data with heritabilities of 0.3 and 0.1 for dairy cattle populations that had more missing pedigrees for animals of lesser genetic merit. Predictions for the youngest generation and UPG solutions were compared with the true values for validation. For both traits, ssGBLUP with metafounders provided accurate and unbiased predictions for young animals while also appropriately accounting for genetic trend. Accuracy was least and bias was greatest for ssGBLUP with UPG for H for the trait with heritability of 0.3 and with UPG for A for the trait with heritability of 0.1. For the trait with heritability of 0.1 and UPG for H, the UPG accuracy (SD) was -0.49 (0.12), suggesting poor estimates of genetic trend despite having little bias for validations on young, genotyped animals. Problems with UPG estimates were likely caused by the lesser amount of information available for the lower heritability trait. Hence, UPG need to be defined differently based on the trait and amount of information. More research is needed to investigate accounting for UPG in A22 to better account for missing pedigrees for genotyped animals.