Roel F. Veerkamp

Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle

The 1000 bull genomes project supports the goal of accelerating the rates of genetic gain in domestic cattle while at the same time considering animal health and welfare by providing the annotated sequence variants and genotypes of key ancestor bulls. In the first phase of the 1000 bull genomes project, we sequenced the whole genomes of 234 cattle to an average of 8.3-fold coverage. This sequencing includes data for 129 individuals from the global Holstein-Friesian population, 43 individuals from the Fleckvieh breed and 15 individuals from the Jersey breed. We identified a total of 28.3 million variants, with an average of 1.44 heterozygous sites per kilobase for each individual. We demonstrate the use of this database in identifying a recessive mutation underlying embryonic death and a dominant mutation underlying lethal chrondrodysplasia. We also performed genome-wide association studies for milk production and curly coat, using imputed sequence variants, and identified variants associated with these traits in cattle.

Dairy cattle breeding objectives combining yield, survival and calving interval for pasture-based systems in Ireland under different milk quota scenarios

Economic values in Irish pounds (IR£) were calculated for milk, fat and protein yields (IR£ cow−1 year−1 kg−1), survival (IR£ cow−1 year−1 % survival−1), and calving interval (IR£ cow−1 year−1 day−1) using a farm model for pasture-based systems in Ireland. Economic values for yield were calculated by changing one of these traits whilst keeping the other traits at the default level, i.e. the economic value for calving interval does not include the costs of higher culling due to a longer calving interval. Herd parameters (e.g. number of milking cows, replacements, youngstock and calving pattern), milk production, energy requirements, feeding ration, land use and labour requirements were re-adjusted to calculate economic performance. Three scenarios were simulated as follows: (S1) milk and fat% quota with a fixed number of cows per farm and quota leasing; (S2) non-quota scenario with a fixed number of cows per farm; and (S3) milk and fat% quota with a fixed output per farm. Sensitivity analysis showed that, with high quota costs and/or lower fat prices, the weighting of fat yield was close to becoming negative in the quota leasing scenario (S1). There was major re-rankings among the top 1000 bulls due to a negative weighting on fat. The results for S1 suggested that the economic values in profit per cow per year are −IR£0.06 per kg milk, IR£0.68 per kg butterfat, IR£4.49 per kg protein, IR£8.98 per % survival, and −IR£1.63 per day calving interval, and the corresponding values in additive genetic standard deviation units were −0.47, 0.19, 1.00, 0.56 and −0.21, respectively.

Accuracy of multi-trait genomic selection using different methods

BackgroundGenomic selection has become a very important tool in animal genetics and is rapidly emerging in plant genetics. It holds the promise to be particularly beneficial to select for traits that are difficult or expensive to measure, such as traits that are measured in one environment and selected for in another environment. The objective of this paper was to develop three models that would permit multi-trait genomic selection by combining scarcely recorded traits with genetically correlated indicator traits, and to compare their performance to single-trait models, using simulated datasets.MethodsThree (SNP) Single Nucleotide Polymorphism based models were used. Model G and BCπ0 assumed that contributed (co)variances of all SNP are equal. Model BSSVS sampled SNP effects from a distribution with large (or small) effects to model SNP that are (or not) associated with a quantitative trait locus. For reasons of comparison, model A including pedigree but not SNP information was fitted as well.ResultsIn terms of accuracies for animals without phenotypes, the models generally ranked as follows: BSSVS > BCπ0 > G > > A. Using multi-trait SNP-based models, the accuracy for juvenile animals without any phenotypes increased up to 0.10. For animals with phenotypes on an indicator trait only, accuracy increased up to 0.03 and 0.14, for genetic correlations with the evaluated trait of 0.25 and 0.75, respectively.ConclusionsWhen the indicator trait had a genetic correlation lower than 0.5 with the trait of interest in our simulated data, the accuracy was higher if genotypes rather than phenotypes were obtained for the indicator trait. However, when genetic correlations were higher than 0.5, using an indicator trait led to higher accuracies for selection candidates. For different combinations of traits, the level of genetic correlation below which genotyping selection candidates is more effective than obtaining phenotypes for an indicator trait, needs to be derived considering at least the heritabilities and the numbers of animals recorded for the traits involved.

The Effect of Linkage Disequilibrium and Family Relationships on the Reliability of Genomic Prediction

Although the concept of genomic selection relies on linkage disequilibrium (LD) between quantitative trait loci and markers, reliability of genomic predictions is strongly influenced by family relationships. In this study, we investigated the effects of LD and family relationships on reliability of genomic predictions and the potential of deterministic formulas to predict reliability using population parameters in populations with complex family structures. Five groups of selection candidates were simulated by taking different information sources from the reference population into account: (1) allele frequencies, (2) LD pattern, (3) haplotypes, (4) haploid chromosomes, and (5) individuals from the reference population, thereby having real family relationships with reference individuals. Reliabilities were predicted using genomic relationships among 529 reference individuals and their relationships with selection candidates and with a deterministic formula where the number of effective chromosome segments (Me) was estimated based on genomic and additive relationship matrices for each scenario. At a heritability of 0.6, reliabilities based on genomic relationships were 0.002 ± 0.0001 (allele frequencies), 0.022 ± 0.001 (LD pattern), 0.018 ± 0.001 (haplotypes), 0.100 ± 0.008 (haploid chromosomes), and 0.318 ± 0.077 (family relationships). At a heritability of 0.1, relative differences among groups were similar. For all scenarios, reliabilities were similar to predictions with a deterministic formula using estimated Me. So, reliabilities can be predicted accurately using empirically estimated Me and level of relationship with reference individuals has a much higher effect on the reliability than linkage disequilibrium per se. Furthermore, accumulated length of shared haplotypes is more important in determining the reliability of genomic prediction than the individual shared haplotype length.

Accuracy of imputation to whole-genome sequence data in Holstein Friesian cattle

BackgroundThe use of whole-genome sequence data can lead to higher accuracy in genome-wide association studies and genomic predictions. However, to benefit from whole-genome sequence data, a large dataset of sequenced individuals is needed. Imputation from SNP panels, such as the Illumina BovineSNP50 BeadChip and Illumina BovineHD BeadChip, to whole-genome sequence data is an attractive and less expensive approach to obtain whole-genome sequence genotypes for a large number of individuals than sequencing all individuals. Our objective was to investigate accuracy of imputation from lower density SNP panels to whole-genome sequence data in a typical dataset for cattle.MethodsWhole-genome sequence data of chromosome 1 (1737 471 SNPs) for 114 Holstein Friesian bulls were used. Beagle software was used for imputation from the BovineSNP50 (3132 SNPs) and BovineHD (40 492 SNPs) beadchips. Accuracy was calculated as the correlation between observed and imputed genotypes and assessed by five-fold cross-validation. Three scenarios S40, S60 and S80 with respectively 40%, 60%, and 80% of the individuals as reference individuals were investigated.ResultsMean accuracies of imputation per SNP from the BovineHD panel to sequence data and from the BovineSNP50 panel to sequence data for scenarios S40 and S80 ranged from 0.77 to 0.83 and from 0.37 to 0.46, respectively. Stepwise imputation from the BovineSNP50 to BovineHD panel and then to sequence data for scenario S40 improved accuracy per SNP to 0.65 but it varied considerably between SNPs.ConclusionsAccuracy of imputation to whole-genome sequence data was generally high for imputation from the BovineHD beadchip, but was low from the BovineSNP50 beadchip. Stepwise imputation from the BovineSNP50 to the BovineHD beadchip and then to sequence data substantially improved accuracy of imputation. SNPs with a low minor allele frequency were more difficult to impute correctly and the reliability of imputation varied more. Linkage disequilibrium between an imputed SNP and the SNP on the lower density panel, minor allele frequency of the imputed SNP and size of the reference group affected imputation reliability.

Association of leptin gene polymorphisms with serum leptin concentration in dairy cows

Leptin is a hormone produced by adipocytes, and its expression is regulated by body fatness and energy balance. This study describes the association of four leptin gene polymorphisms in dairy cows (R4C, A59V, RFLP1, and BM1500) with circulating leptin concentrations during the periparturient period. A59V is located at a between-species conserved region of leptin, and R4C might have effect on the tertiary structure of the leptin protein because of the presence of an extra cystein. RFLP1 is an intronic SNP and BM1500 is a microsatellite located 3.6 kb downstream of the leptin locus. The four polymorphisms were genotyped in 323 HF heifers with known pedigree. Leptin concentrations were determined biweekly from 30 days before until 80 days after parturition. The effect of genotype on leptin concentrations was modeled by fitting a spline in ASREML describing leptin concentrations as a function of days relative to parturition for each genotype/allele. Surprisingly, associations were found during pregnancy, but not during lactation. This indicates that the polymorphism could be more effective during pregnancy. If further studies demonstrate that more leptin-binding protein (Ob-Re) is present in this stage, it is hypothesized that a structural difference in the leptin protein could cause a sub-optimal binding stringency to Ob-Re. Free leptin could be cleared faster than bound leptin, and this could result in lower leptin concentrations during pregnancy for the polymorphism. The effects found might be ascribed to R4C. However, more study on the Ob-Re receptor, like binding stringencies between R4C and wild-type leptin and glycosylation during pregnancy, would provide more insight in the results found.

The use of mid-infrared spectrometry to predict body energy status of Holstein cows

Energy balance, especially in early lactation, is known to be associated with subsequent health and fertility in dairy cows. However, its inclusion in routine management decisions or breeding programs is hindered by the lack of quick, easy, and inexpensive measures of energy balance. The objective of this study was to evaluate the potential of mid-infrared (MIR) analysis of milk, routinely available from all milk samples taken as part of large-scale milk recording and milk payment operations, to predict body energy status and related traits in lactating dairy cows. The body energy status traits investigated included energy balance and body energy content. The related traits of body condition score and energy intake were also considered. Measurements on these traits along with milk MIR spectral data were available on 17 different test days from 268 cows (418 lactations) and were used to develop the prediction equations using partial least squares regression. Predictions were externally validated on different independent subsets of the data and the results averaged. The average accuracy of predicting body energy status from MIR spectral data was as high as 75% when energy balance was measured across lactation. These predictions of body energy status were considerably more accurate than predictions obtained from the sometimes proposed fat-to-protein ratio in milk. It is not known whether the prediction generated from MIR data are a better reflection of the true (unknown) energy status than the actual energy status measures used in this study. However, results indicate that the approach described may be a viable method of predicting individual cow energy status for a large scale of application.

Improved accuracy of genomic prediction for dry matter intake of dairy cattle from combined European and Australian data sets

With the aim of increasing the accuracy of genomic estimated breeding values for dry matter intake (DMI) in dairy cattle, data from Australia (AU), the United Kingdom (UK), and the Netherlands (NL) were combined using both single-trait and multi-trait models. In total, DMI records were available on 1,801 animals, including 843 AU growing heifers with records on DMI measured over 60 to 70 d at approximately 200 d of age, and 359 UK and 599 NL lactating heifers with records on DMI during the first 100 d in milk. The genotypes used in this study were obtained from the Illumina Bovine 50K chip (Illumina Inc., San Diego, CA). The AU, UK, and NL genomic data were matched using the single nucleotide polymorphism (SNP) name. Quality controls were applied by carefully comparing the genotypes of 40 bulls that were available in each data set. This resulted in 30,949 SNP being used in the analyses. Genomic predictions were estimated with genomic REML, using ASReml software. The accuracy of genomic prediction was evaluated in 11 validation sets; that is, at least 3 validation sets per country were defined. The reference set (in which animals had both DMI phenotypes and genotypes) was either AU or Europe (UK and NL) or a multi-country reference set consisting of all data except the validation set. When DMI for each country was treated as the same trait, use of a multi-country reference set increased the accuracy of genomic prediction for DMI in UK, but not in AU and NL. Extending the model to a bivariate (AU-EU) or trivariate (AU-UK-NL) model increased the accuracy of genomic prediction for DMI in all countries. The highest accuracies were estimated for all countries when data were analyzed with a trivariate model, with increases of up to 5.5% compared with univariate models within countries.

Genomic selection for feed efficiency in dairy cattle.

Feed is a major component of variable costs associated with dairy systems and is therefore an important consideration for breeding objectives. As a result, measures of feed efficiency are becoming popular traits for genetic analyses. Already, several countries account for feed efficiency in their breeding objectives by approximating the amount of energy required for milk production, maintenance, etc. However, variation in actual feed intake is currently not captured in dairy selection objectives, although this could be possible by evaluating traits such as residual feed intake (RFI), defined as the difference between actual and predicted feed (or energy) intake. As feed intake is expensive to accurately measure on large numbers of cows, phenotypes derived from it are obvious candidates for genomic selection provided that: (1) the trait is heritable; (2) the reliability of genomic predictions are acceptable to those using the breeding values; and (3) if breeding values are estimated for heifers, rather than cows then the heifer and cow traits need to be correlated. The accuracy of genomic prediction of dry matter intake (DMI) and RFI has been estimated to be around 0.4 in beef and dairy cattle studies. There are opportunities to increase the accuracy of prediction, for example, pooling data from three research herds (in Australia and Europe) has been shown to increase the accuracy of genomic prediction of DMI from 0.33 within country to 0.35 using a three-country reference population. Before including RFI as a selection objective, genetic correlations with other traits need to be estimated. Weak unfavourable genetic correlations between RFI and fertility have been published. This could be because RFI is mathematically similar to the calculation of energy balance and failure to account for mobilisation of body reserves correctly may result in selection for a trait that is similar to selecting for reduced (or negative) energy balance. So, if RFI is to become a selection objective, then including it in an overall multi-trait selection index where the breeding objective is net profit is sensible, as this would allow genetic correlations with other traits to be properly accounted for. If genetic parameters are accurately estimated then RFI is a logical breeding objective. If there is uncertainty in these, then DMI may be preferable.

International genetic evaluations for feed intake in dairy cattle through the collation of data from multiple sources

Feed represents a large proportion of the variable costs in dairy production systems. The omission of feed intake measures explicitly from national dairy cow breeding objectives is predominantly due to a lack of information from which to make selection decisions. However, individual cow feed intake data are available in different countries, mostly from research or nucleus herds. None of these data sets are sufficiently large enough on their own to generate accurate genetic evaluations. In the current study, we collate data from 10 populations in 9 countries and estimate genetic parameters for dry matter intake (DMI). A total of 224,174 test-day records from 10,068 parity 1 to 5 records of 6,957 cows were available, as well as records from 1,784 growing heifers. Random regression models were fit to the lactating cow test-day records and predicted feed intake at 70 d postcalving was extracted from these fitted profiles. The random regression model included a fixed polynomial regression for each lactation separately, as well as herd-year-season of calving and experimental treatment as fixed effects; random effects fit in the model included individual animal deviation from the fixed regression for each parity as well as mean herd-specific deviations from the fixed regression. Predicted DMI at 70 d postcalving was used as the phenotype for the subsequent genetic analyses undertaken using an animal repeatability model. Heritability estimates of predicted cow feed intake 70 d postcalving was 0.34 across the entire data set and varied, within population, from 0.08 to 0.52. Repeatability of feed intake across lactations was 0.66. Heritability of feed intake in the growing heifers was 0.20 to 0.34 in the 2 populations with heifer data. The genetic correlation between feed intake in lactating cows and growing heifers was 0.67. A combined pedigree and genomic relationship matrix was used to improve linkages between populations for the estimation of genetic correlations of DMI in lactating cows; genotype information was available on 5,429 of the animals. Populations were categorized as North America, grazing, other low input, and high input European Union. Albeit associated with large standard errors, genetic correlation estimates for DMI between populations varied from 0.14 to 0.84 but were stronger (0.76 to 0.84) between the populations representative of high-input production systems. Genetic correlations with the grazing populations were weak to moderate, varying from 0.14 to 0.57. Genetic evaluations for DMI can be undertaken using data collated from international populations; however, genotype-by-environment interactions with grazing production systems need to be considered.