Imputation strategies for genomic prediction using nanopore sequencing.

BACKGROUND: Genomic prediction describes the use of SNP genotypes to predict complex traits and has been widely applied in humans and agricultural species. Genotyping-by-sequencing, a method which uses low-coverage sequence data paired with genotype imputation, is becoming an increasingly popular SNP genotyping method for genomic prediction. The development of Oxford Nanopore Technologies' (ONT) MinION sequencer has now made genotyping-by-sequencing portable and rapid. Here we evaluate the speed and accuracy of genomic predictions using low-coverage ONT sequence data in a population of cattle using four imputation approaches. We also investigate the effect of SNP reference panel size on imputation performance. RESULTS: SNP array genotypes and ONT sequence data for 62 beef heifers were used to calculate genomic estimated breeding values (GEBVs) from 641 k SNP for four traits. GEBV accuracy was much higher when genome-wide flanking SNP from sequence data were used to help impute the 641 k panel used for genomic predictions. Using the imputation package QUILT, correlations between ONT and low-density SNP array genomic breeding values were greater than 0.91 and up to 0.97 for sequencing coverages as low as 0.1 × using a reference panel of 48 million SNP. Imputation time was significantly reduced by decreasing the number of flanking sequence SNP used in imputation for all methods. When compared to high-density SNP arrays, genotyping accuracy and genomic breeding value correlations at 0.5 × coverage were also found to be higher than those imputed from low-density arrays. CONCLUSIONS: Here we demonstrated accurate genomic prediction is possible with ONT sequence data from sequencing coverages as low as 0.1 × , and imputation time can be as short as 10 min per sample. We also demonstrate that in this population, genotyping-by-sequencing at 0.1 × coverage can be more accurate than imputation from low-density SNP arrays.

Developing flexible models for genetic evaluations in smallholder crossbred dairy farms.

The productivity of smallholder dairy farms is very low in developing countries. Important genetic gains could be realized using genomic selection, but genetic evaluations need to be tailored for lack of pedigree information and very small farm sizes. To accommodate this situation, we propose a flexible Bayesian model for the genetic evaluation of milk yield, which allows us to simultaneously account for nongenetic random effects for farms and varying SNP variance (BayesR model). First, we used simulations based on real genotype data from Indian crossbred dairy cattle to demonstrate that the proposed model can separate the true genetic and nongenetic parameters even for small farm sizes (2 cows on average) although with high standard errors in scenarios with low heritability. The accuracy of genomic genetic evaluation increased until farm size was approximately 5. We then applied the model to real data from 4,655 crossbred cows with 106,109 monthly test day milk records and 689,750 autosomal SNPs. We estimated a heritability of 0.16 (0.04) for milk yield and using cross-validation, a genomic estimated breeding value (GEBV) accuracy of 0.45 and bias (regression of phenotype on GEBV) of 1.04 (0.26). Estimated genetic parameters were very similar using BayesR, BayesC, and genomic BLUP approaches. Candidate genes near the top variants, IMMP2L and ARHGEF2, have been previously associated with milk protein composition, mastitis resistance, and milk cholesterol content. The estimated heritability and GEBV accuracy for milk yield are much lower than those from intensive or pasture-based systems in many countries. Further increases in the number of phenotyped and genotyped animals in farms with at least 2 cows (preferably 3-5, to allow for dropout of cows) are needed to improve the estimation of genetic effects in these smallholder dairy farms.

CELPI: trial protocol for a randomised controlled trial of a Carer End of Life Planning Intervention in people dying with dementia.

BACKGROUND: Dementia is a leading cause of death in developed nations. Despite an often distressing and symptom laden end of life, there are systematic barriers to accessing palliative care in older people dying of dementia. Evidence exists that 70% of people living with severe dementia attend an emergency department (ED) in their last year of life. The aim of this trial is to test whether a Carer End of Life Planning Intervention (CELPI), co-designed by consumers, clinicians and content specialists, improves access to end of life care for older people with severe dementia, using an ED visit as a catalyst for recognising unmet needs and specialist palliative care referral where indicated. METHODS: A randomised controlled trial (RCT) enrolling at six EDs across three states in Australia will be conducted, enrolling four hundred and forty dyads comprising a person with severe dementia aged ≥ 65 years, and their primary carer. Participants will be randomly allocated to CELPI or the control group. CELPI incorporates a structured carer needs assessment and referral to specialist palliative care services where indicated by patient symptom burden and needs assessment. The primary outcome measure is death of the person with dementia in the carer-nominated preferred location. Secondary outcomes include carer reported quality of life of the person dying of dementia, hospital bed day occupancy in the last 12 months of life, and carer stress. An economic evaluation from the perspective of a health funder will be conducted. DISCUSSION: CELPI seeks to support carers and provide optimal end of life care for the person dying of dementia. This trial will provide high level evidence as to the clinical and cost effectiveness of this intervention. TRIAL REGISTRATION: ACTRN12622000611729 registered 22/04/2022.

Overlap between eQTL and QTL associated with production traits and fertility in dairy cattle.

BACKGROUND: Identifying causative mutations or genes through which quantitative trait loci (QTL) act has proven very difficult. Using information such as gene expression may help to identify genes and mutations underlying QTL. Our objective was to identify regions associated both with production traits or fertility and with gene expression, in dairy cattle. We used three different approaches to discover QTL that are also expression QTL (eQTL): 1) estimate the correlation between local genomic estimated breeding values (GEBV) and gene expression, 2) investigate whether the 300 intervals explaining most genetic variance for a trait contain more eQTL than 300 randomly selected intervals, and 3) a colocalisation analysis. Phenotypes and genotypes up to sequence level of 35,775 dairy bulls and cows were used for QTL mapping, and gene expression and genotypes of 131 cows were used to identify eQTL. RESULTS: With all three approaches, we identified some overlap between eQTL and QTL, though the majority of QTL in our dataset did not seem to be eQTL. The most significant associations between QTL and eQTL were found for intervals on chromosome 18, where local GEBV for all traits showed a strong association with the expression of the FUK and DDX19B. Intervals whose local GEBV for a trait correlated highly significantly with the expression of a nearby gene explained only a very small part of the genetic variance for that trait. It is likely that part of these correlations were due to linkage disequilibrium (LD) in the interval. While the 300 intervals explaining most genetic variance explained most of the GEBV variance, they contained only slightly more eQTL than 300 randomly selected intervals that explained a minimal portion of the GEBV variance. Furthermore, some variants showed a high colocalisation probability, but this was only the case for few variants. CONCLUSIONS: Several reasons may have contributed to the low level of overlap between QTL and eQTL detected in our study, including a lack of power in the eQTL study and long-range LD making it difficult to separate QTL and eQTL. Furthermore, it may be that eQTL explain only a small fraction of QTL.

The cost-benefit of genomic testing of heifers and using sexed semen in pasture-based dairy herds.

Recent improvements in dairy cow fertility and female reproductive technologies offer an opportunity to apply greater selection pressure to females. This means there may be greater incentive to obtain genomic breeding values for females. We modeled the impact of changes to key parameters on the net benefit from genomic testing of heifer calves with and without usage of sexed semen. This paper builds on earlier cost-benefit studies but uses parameters relevant to pasture-based systems. A deterministic model was used to evaluate the effect on net benefit due to changes in (1) reproduction rate, (2) genomic test costs, (3) availability of parent-derived breeding values (EBVPA), and (4) replacement rate. When the use of sexed semen was included, we also considered (1) the proportion of heifers and cows mated to sexed semen, (2) decreases in conception rate in inseminations with sexed semen, and (3) the marginal return for surplus heifers. Scenarios with lower replacement rates and no availability of EBVPA had the largest net benefits. Under current Australian parameters, the net benefit of genomic testing realized over the lifetime of genotyped heifers is expected to range from A$204 to A$1,124 per 100 cows for a herd with median reproductive performance. The cost of a genomic test, a perceived barrier to many farmers, had only a small effect on net benefit. Genomic testing alone was always more profitable than using sexed semen and genomic testing together if the only benefit considered was increased genetic gain in heifer replacements. When other benefits (i.e., the higher sale price of a surplus heifer compared with a male calf) were considered, there were combinations of parameters where net benefit from using sexed semen and genomic testing was higher than the equivalent scenario with genomic testing only. Using sexed semen alongside genomic testing is most likely to be profitable when (1) used in heifers, (2) the marginal return for selling surplus heifers (sale price minus rearing costs) is greater than A$400, and (3) conception rates of no more than 10 percentage points lower than those achieved using conventional semen can be realized. Net benefit was highly dependent on the marginal return. Demonstrating that the initial investment in genomic testing can be recouped within the lifetime of the heifers tested may assist in the development of extension messages to explain the value of genomic testing females at the herd level.

Multibreed genomic prediction using multitrait genomic residual maximum likelihood and multitask Bayesian variable selection.

Genomic prediction is applicable to individuals of different breeds. Empirical results to date, however, show limited benefits in using information on multiple breeds in the context of genomic prediction. We investigated a multitask Bayesian model, presented previously by others, implemented in a Bayesian stochastic search variable selection (BSSVS) model. This model allowed for evidence of quantitative trait loci (QTL) to be accumulated across breeds or for both QTL that segregate across breeds and breed-specific QTL. In both cases, single nucleotide polymorphism effects were estimated with information from a single breed. Other models considered were a single-trait and multitrait genomic residual maximum likelihood (GREML) model, with breeds considered as different traits, and a single-trait BSSVS model. All single-trait models were applied to each of the 2 breeds separately and to the pooled data of both breeds. The data used included a training data set of 6,278 Holstein and 722 Jersey bulls, as well as 374 Jersey validation bulls. All animals had genotypes for 474,773 single nucleotide polymorphisms after editing and phenotypes for milk, fat, and protein yields. Using the same training data, BSSVS consistently outperformed GREML. The multitask BSSVS, however, did not outperform single-trait BSSVS, which used pooled Holstein and Jersey data for training. Thus, the rigorous assumption that the traits are the same in both breeds yielded a slightly better prediction than a model that had to estimate the correlation between the breeds from the data. Adding the Holstein data significantly increased the accuracy of the single-trait GREML and BSSVS in predicting the Jerseys for milk and protein, in line with estimated correlations between the breeds of 0.66 and 0.47 for milk and protein yields, whereas only the BSSVS model significantly improved the accuracy for fat yield with an estimated correlation between breeds of only 0.05. The relatively high genetic correlations for milk and protein yields, and the superiority of the pooling strategy, is likely the result of the observed admixture between both breeds in our data. The Bayesian model was able to detect several QTL in Holsteins, which likely enabled it to outperform GREML. The inability of the multitask Bayesian models to outperform a simple pooling strategy may be explained by the fact that the pooling strategy assumes equal effects in both breeds; furthermore, this assumption may be valid for moderate- to large-sized QTL, which are important for multibreed genomic prediction.

Responses of dairy cows with divergent residual feed intake as calves to metabolic challenges during midlactation and the nonlactating period.

Residual feed intake (RFI) is defined as the difference between the actual and expected feed intake required to support animal maintenance and growth. Thus, a cow with a low RFI can obtain nutrients for maintenance and growth from a reduced amount of feed compared with a cow with a high RFI. Variation in RFI is underpinned by a combination of factors, including genetics, metabolism, thermoregulation and body composition; hypothalamic-pituitary-adrenal (HPA) axis responsiveness is also a possible contributor. Responses to 3 metabolic challenges were measured in lactating and nonlactating dairy cattle. Sixteen Holstein Friesian cows with phenotypic RFI measurements that were obtained during the growth period (188-220 d old) were grouped as either low-calfhood RFI (n = 8) or high-calfhood RFI (n = 8). An ACTH (2 µg/kg of body weight), insulin (0.12 U/kg), and epinephrine (a low dose of 0.1 µg/kg and a high dose of 1.6 µg/kg of epinephrine) challenge were each conducted during both midlactation (122 ± 23.4 d in milk) and the nonlactating period (dry period; approximately 38 d after cessation of milking). Cows were housed in metabolism stalls for the challenges and were fed a diet of alfalfa cubes ad libitum for at least 10 d before the experiment (lactating cows also were offered a total of 6 kg of dry matter/d of crushed wheat grain plus minerals fed as 3 kg of dry matter at each milking) and were fasted for 12 h before the challenges. The efficiency of conversion of feed into milk (the ratio of feed consumed to milk produced over the 7 d before the experiment) during midlactation was better (lower) in low-calfhood RFI cows, although dry matter intake did not differ between RFI groups. Low-calfhood RFI cows exhibited a lower plasma cortisol response to the ACTH challenge than high-calfhood RFI cows, particularly in midlactation (-15%). The low-calfhood RFI cows had a greater plasma insulin-like growth factor-1 response to the insulin challenge and plasma fatty acid response to epinephrine compared with the high-calfhood RFI cows. These data suggest that high-calfhood RFI cows exhibit a more responsive HPA axis. As divergence in RFI measured during growth is retained (although reduced) during lactation, it is possible that energy is used to respond to HPA axis activation at the expense of production in high-calfhood RFI dairy cattle during lactation and contributes to a decrease in overall feed use efficiency.

Accelerating wheat breeding for end-use quality with multi-trait genomic predictions incorporating near infrared and nuclear magnetic resonance-derived phenotypes.

KEY MESSAGE: Using NIR and NMR predictions of quality traits overcomes a major barrier for the application of genomic selection to accelerate improvement in grain end-use quality traits of wheat. Grain end-use quality traits are among the most important in wheat breeding. These traits are difficult to breed for, as their assays require flour quantities only obtainable late in the breeding cycle, and are expensive. These traits are therefore an ideal target for genomic selection. However, large reference populations are required for accurate genomic predictions, which are challenging to assemble for these traits for the same reasons they are challenging to breed for. Here, we use predictions of end-use quality derived from near infrared (NIR) or nuclear magnetic resonance (NMR), that require very small amounts of flour, as well as end-use quality measured by industry standard assay in a subset of accessions, in a multi-trait approach for genomic prediction. The NIR and NMR predictions were derived for 19 end-use quality traits in 398 accessions, and were then assayed in 2420 diverse wheat accessions. The accessions were grown out in multiple locations and multiple years, and were genotyped for 51208 SNP. Incorporating NIR and NMR phenotypes in the multi-trait approach increased the accuracy of genomic prediction for most quality traits. The accuracy ranged from 0 to 0.47 before the addition of the NIR/NMR data, while after these data were added, it ranged from 0 to 0.69. Genomic predictions were reasonably robust across locations and years for most traits. Using NIR and NMR predictions of quality traits overcomes a major barrier for the application of genomic selection for grain end-use quality traits in wheat breeding.

Genomic prediction of reproduction traits for Merino sheep.

Economically important reproduction traits in sheep, such as number of lambs weaned and litter size, are expressed only in females and later in life after most selection decisions are made, which makes them ideal candidates for genomic selection. Accurate genomic predictions would lead to greater genetic gain for these traits by enabling accurate selection of young rams with high genetic merit. The aim of this study was to design and evaluate the accuracy of a genomic prediction method for female reproduction in sheep using daughter trait deviations (DTD) for sires and ewe phenotypes (when individual ewes were genotyped) for three reproduction traits: number of lambs born (NLB), litter size (LSIZE) and number of lambs weaned. Genomic best linear unbiased prediction (GBLUP), BayesR and pedigree BLUP analyses of the three reproduction traits measured on 5340 sheep (4503 ewes and 837 sires) with real and imputed genotypes for 510 174 SNPs were performed. The prediction of breeding values using both sire and ewe trait records was validated in Merino sheep. Prediction accuracy was evaluated by across sire family and random cross-validations. Accuracies of genomic estimated breeding values (GEBVs) were assessed as the mean Pearson correlation adjusted by the accuracy of the input phenotypes. The addition of sire DTD into the prediction analysis resulted in higher accuracies compared with using only ewe records in genomic predictions or pedigree BLUP. Using GBLUP, the average accuracy based on the combined records (ewes and sire DTD) was 0.43 across traits, but the accuracies varied by trait and type of cross-validations. The accuracies of GEBVs from random cross-validations (range 0.17-0.61) were higher than were those from sire family cross-validations (range 0.00-0.51). The GEBV accuracies of 0.41-0.54 for NLB and LSIZE based on the combined records were amongst the highest in the study. Although BayesR was not significantly different from GBLUP in prediction accuracy, it identified several candidate genes which are known to be associated with NLB and LSIZE. The approach provides a way to make use of all data available in genomic prediction for traits that have limited recording.

Including nonadditive genetic effects in mating programs to maximize dairy farm profitability.

We compared the outcome of mating programs based on different evaluation models that included nonadditive genetic effects (dominance and heterozygosity) in addition to additive effects. The additive and dominance marker effects and the values of regression on average heterozygosity were estimated using 632,003 single nucleotide polymorphisms from 7,902 and 7,510 Holstein cows with calving interval and production (milk, fat, and protein yields) records, respectively. Expected progeny values were computed based on the estimated genetic effects and genotype probabilities of hypothetical progeny from matings between the available genotyped cows and the top 50 young genomic bulls. An index combining the traits based on their economic values was developed and used to evaluate the performance of different mating scenarios in terms of dollar profit. We observed that mating programs with nonadditive genetic effects performed better than a model with only additive effects. Mating programs with dominance and heterozygosity effects increased milk, fat, and protein yields by up to 38, 1.57, and 1.21 kg, respectively. The inclusion of dominance and heterozygosity effects decreased calving interval by up to 0.70 d compared with random mating. The average reduction in progeny inbreeding by the inclusion of nonadditive genetic effects in matings compared with random mating was between 0.25 to 1.57 and 0.64 to 1.57 percentage points for calving interval and production traits, respectively. The reduction in inbreeding was accompanied by an average of A$8.42 (Australian dollars) more profit per mating for a model with additive, dominance, and heterozygosity effects compared with random mating. Mate allocations that benefit from nonadditive genetic effects can improve progeny performance only in the generation where it is being implemented, and the gain from specific combining abilities cannot be accumulated over generations. Continuous updating of genomic predictions and mate allocation programs are required to benefit from nonadditive genetic effects in the long term.

Exploiting biological priors and sequence variants enhances QTL discovery and genomic prediction of complex traits.

BACKGROUND: Dense SNP genotypes are often combined with complex trait phenotypes to map causal variants, study genetic architecture and provide genomic predictions for individuals with genotypes but no phenotype. A single method of analysis that jointly fits all genotypes in a Bayesian mixture model (BayesR) has been shown to competitively address all 3 purposes simultaneously. However, BayesR and other similar methods ignore prior biological knowledge and assume all genotypes are equally likely to affect the trait. While this assumption is reasonable for SNP array genotypes, it is less sensible if genotypes are whole-genome sequence variants which should include causal variants. RESULTS: We introduce a new method (BayesRC) based on BayesR that incorporates prior biological information in the analysis by defining classes of variants likely to be enriched for causal mutations. The information can be derived from a range of sources, including variant annotation, candidate gene lists and known causal variants. This information is then incorporated objectively in the analysis based on evidence of enrichment in the data. We demonstrate the increased power of BayesRC compared to BayesR using real dairy cattle genotypes with simulated phenotypes. The genotypes were imputed whole-genome sequence variants in coding regions combined with dense SNP markers. BayesRC increased the power to detect causal variants and increased the accuracy of genomic prediction. The relative improvement for genomic prediction was most apparent in validation populations that were not closely related to the reference population. We also applied BayesRC to real milk production phenotypes in dairy cattle using independent biological priors from gene expression analyses. Although current biological knowledge of which genes and variants affect milk production is still very incomplete, our results suggest that the new BayesRC method was equal to or more powerful than BayesR for detecting candidate causal variants and for genomic prediction of milk traits. CONCLUSIONS: BayesRC provides a novel and flexible approach to simultaneously improving the accuracy of QTL discovery and genomic prediction by taking advantage of prior biological knowledge. Approaches such as BayesRC will become increasing useful as biological knowledge accumulates regarding functional regions of the genome for a range of traits and species.

Leveraging genetically simple traits to identify small-effect variants for complex phenotypes.

BACKGROUND: Polymorphisms underlying complex traits often explain a small part (less than 1 %) of the phenotypic variance (σ2P). This makes identification of mutations underling complex traits difficult and usually only a subset of large-effect loci are identified. One approach to identify more loci is to increase sample size of experiments but here we propose an alternative. The aim of this paper is to use secondary phenotypes for genetically simple traits during the QTL discovery phase for complex traits. We demonstrate this approach in a dairy cattle data set where the complex traits were milk production phenotypes (fat, milk and protein yield; fat and protein percentage in milk) measured on thousands of individuals while secondary (potentially genetically simpler) traits are detailed milk composition traits (measurements of individual protein abundance, mineral and sugar concentrations; and gene expression). RESULTS: Quantitative trait loci (QTL) were identified using 11,527 Holstein cattle with milk production records and up to 444 cows with milk composition traits. There were eight regions that contained QTL for both milk production and a composition trait, including four novel regions. One region on BTAU1 affected both milk yield and phosphorous concentration in milk. The QTL interval included the gene SLC37A1, a phosphorous antiporter. The most significant imputed sequence variants in this region explained 0.001 σ2P for milk yield, and 0.11 σ2P for phosphorus concentration. Since the polymorphisms were non-coding, association mapping for SLC37A1 gene expression was performed using high depth mammary RNAseq data from a separate group of 371 lactating cows. This confirmed a strong eQTL for SLC37A1, with peak association at the same imputed sequence variants that were most significant for phosphorus concentration. Fitting any of these variants as covariables in the association analysis removed the QTL signal for milk production traits. Plausible causative mutations in the casein complex region were also identified using a similar strategy. CONCLUSIONS: Milk production traits in dairy cows are typical complex traits where polymorphisms explain only a small portion of the phenotypic variance. However, here we show that these mutations can have larger effects on secondary traits, such as concentrations of minerals, proteins and sugars in the milk, and expression levels of genes in mammary tissue. These larger effects were used to successfully map variants for milk production traits. Genetically simple traits also provide a direct biological link between possible causal mutations and the effect of these mutations on milk production.

Genetics of complex traits: prediction of phenotype, identification of causal polymorphisms and genetic architecture.

Complex or quantitative traits are important in medicine, agriculture and evolution, yet, until recently, few of the polymorphisms that cause variation in these traits were known. Genome-wide association studies (GWAS), based on the ability to assay thousands of single nucleotide polymorphisms (SNPs), have revolutionized our understanding of the genetics of complex traits. We advocate the analysis of GWAS data by a statistical method that fits all SNP effects simultaneously, assuming that these effects are drawn from a prior distribution. We illustrate how this method can be used to predict future phenotypes, to map and identify the causal mutations, and to study the genetic architecture of complex traits. The genetic architecture of complex traits is even more complex than previously thought: in almost every trait studied there are thousands of polymorphisms that explain genetic variation. Methods of predicting future phenotypes, collectively known as genomic selection or genomic prediction, have been widely adopted in livestock and crop breeding, leading to increased rates of genetic improvement.

Polyceraty (multi-horns) in Damara sheep maps to ovine chromosome 2.

Polyceraty (presence of multiple horns) is rare in modern day ungulates. Although not found in wild sheep, polyceraty does occur in a small number of domestic sheep breeds covering a wide geographical region. Damara are fat-tailed hair sheep, from the south-western region of Africa, which display polyceraty, with horn number ranging from zero to four. We conducted a genome-wide association study for horn number with 43 Damara genotyped with 606 006 SNP markers. The analysis revealed a region with multiple significant SNPs on ovine chromosome 2, in a location different from the mutation for polled in sheep on chromosome 10. The causal mutation for polyceraty was not identified; however, the region associated with polyceraty spans nine HOXD genes, which are critical in embryonic development of appendages. Mutations in HOXD genes are implicated in polydactly phenotypes in mice and humans. There was no evidence for epistatic interactions contributing to polyceraty. This is the first report on the genetic mechanisms underlying polyceraty in the under-studied Damara.

Technical note: Equivalent genomic models with a residual polygenic effect.

Routine genomic evaluations in animal breeding are usually based on either a BLUP with genomic relationship matrix (GBLUP) or single nucleotide polymorphism (SNP) BLUP model. For a multi-step genomic evaluation, these 2 alternative genomic models were proven to give equivalent predictions for genomic reference animals. The model equivalence was verified also for young genotyped animals without phenotypes. Due to incomplete linkage disequilibrium of SNP markers to genes or causal mutations responsible for genetic inheritance of quantitative traits, SNP markers cannot explain all the genetic variance. A residual polygenic effect is normally fitted in the genomic model to account for the incomplete linkage disequilibrium. In this study, we start by showing the proof that the multi-step GBLUP and SNP BLUP models are equivalent for the reference animals, when they have a residual polygenic effect included. Second, the equivalence of both multi-step genomic models with a residual polygenic effect was also verified for young genotyped animals without phenotypes. Additionally, we derived formulas to convert genomic estimated breeding values of the GBLUP model to its components, direct genomic values and residual polygenic effect. Third, we made a proof that the equivalence of these 2 genomic models with a residual polygenic effect holds also for single-step genomic evaluation. Both the single-step GBLUP and SNP BLUP models lead to equal prediction for genotyped animals with phenotypes (e.g., reference animals), as well as for (young) genotyped animals without phenotypes. Finally, these 2 single-step genomic models with a residual polygenic effect were proven to be equivalent for estimation of SNP effects, too.

Design of a low-density SNP chip for the main Australian sheep breeds and its effect on imputation and genomic prediction accuracy.

Genotyping sheep for genome-wide SNPs at lower density and imputing to a higher density would enable cost-effective implementation of genomic selection, provided imputation was accurate enough. Here, we describe the design of a low-density (12k) SNP chip and evaluate the accuracy of imputation from the 12k SNP genotypes to 50k SNP genotypes in the major Australian sheep breeds. In addition, the impact of imperfect imputation on genomic predictions was evaluated by comparing the accuracy of genomic predictions for 15 novel meat traits including carcass and meat quality and omega fatty acid traits in sheep, from 12k SNP genotypes, imputed 50k SNP genotypes and real 50k SNP genotypes. The 12k chip design included 12 223 SNPs with a high minor allele frequency that were selected with intermarker spacing of 50-475 kb. SNPs for parentage and horned or polled tests also were represented. Chromosome ends were enriched with SNPs to reduce edge effects on imputation. The imputation performance of the 12k SNP chip was evaluated using 50k SNP genotypes of 4642 animals from six breeds in three different scenarios: (1) within breed, (2) single breed from multibreed reference and (3) multibreed from a single-breed reference. The highest imputation accuracies were found with scenario 2, whereas scenario 3 was the worst, as expected. Using scenario 2, the average imputation accuracy in Border Leicester, Polled Dorset, Merino, White Suffolk and crosses was 0.95, 0.95, 0.92, 0.91 and 0.93 respectively. Imputation scenario 2 was used to impute 50k genotypes for 10 396 animals with novel meat trait phenotypes to compare genomic prediction accuracy using genomic best linear unbiased prediction (GBLUP) with real and imputed 50k genotypes. The weighted mean imputation accuracy achieved was 0.92. The average accuracy of genomic estimated breeding values (GEBVs) based on only 12k data was 0.08 across traits and breeds, but accuracies varied widely. The mean GBLUP accuracies with imputed 50k data more than doubled to 0.21. Accuracies of genomic prediction were very similar for imputed and real 50k genotypes. There was no apparent impact on accuracy of GEBVs as a result of using imputed rather than real 50k genotypes, provided imputation accuracy was >90%.

Including overseas performance information in genomic evaluations of Australian dairy cattle.

In dairy cattle, the rate of genetic gain from genomic selection depends on reliability of direct genomic values (DGV). One option to increase reliabilities could be to increase the size of the reference set used for prediction, by using genotyped bulls with daughter information in countries other than the evaluating country. The increase in reliabilities of DGV from using this information will depend on the extent of genotype by environment interaction between the evaluating country and countries contributing information, and whether this is correctly accounted for in the prediction method. As the genotype by environment interaction between Australia and Europe or North America is greater than between Europe and North America for most dairy traits, ways of including information from other countries in Australian genomic evaluations were examined. Thus, alternative approaches for including information from other countries and their effect on the reliability and bias of DGV of selection candidates were assessed. We also investigated the effect of including overseas (OS) information on reliabilities of DGV for selection candidates that had weaker relationships to the current Australian reference set. The DGV were predicted either using daughter trait deviations (DTD) for the bulls with daughters in Australia, or using this information as well as OS information by including deregressed proofs (DRP) from Interbull for bulls with only OS daughters in either single trait or bivariate models. In the bivariate models, DTD and DRP were considered as different traits. Analyses were performed for Holstein and Jersey bulls for milk yield traits, fertility, cell count, survival, and some type traits. For Holsteins, the data used included up to 3,580 bulls with DTD and up to 5,720 bulls with only DRP. For Jersey, about 900 bulls with DTD and 1,820 bulls with DRP were used. Bulls born after 2003 and genotyped cows that were not dams of genotyped bulls were used for validation. The results showed that the combined use of DRP on bulls with OS daughters only and DTD for Australian bulls in either the single trait or bivariate model increased the coefficient of determination [(R(2)) (DGV,DTD)] in the validation set, averaged across 6 main traits, by 3% in Holstein and by 5% in Jersey validation bulls relative to the use of DTD only. Gains in reliability and unbiasedness of DGV were similar for the single trait and bivariate models for production traits, whereas the bivariate model performed slightly better for somatic cell count in Holstein. The increase in R(2) (DGV,DTD) as a result of using bulls with OS daughters was relatively higher for those bulls and cows in the validation sets that were less related to the current reference set. For example, in Holstein, the average increase in R(2) for milk yield traits when DTD and DRP were used in a single trait model was 23% in the least-related cow group, but only 3% in the most-related cow group. In general, for both breeds the use of DTD from domestic sources and DRP from Interbull in a single trait or bivariate model can increase reliability of DGV for selection candidates.

Hot topic: Definition and implementation of a breeding value for feed efficiency in dairy cows.

A new breeding value that combines the amount of feed saved through improved metabolic efficiency with predicted maintenance requirements is described. The breeding value includes a genomic component for residual feed intake (RFI) combined with maintenance requirements calculated from either a genomic or pedigree estimated breeding value (EBV) for body weight (BW) predicted using conformation traits. Residual feed intake is only available for genotyped Holsteins; however, BW is available for all breeds. The RFI component of the "feed saved" EBV has 2 parts: Australian calf RFI and Australian lactating cow RFI. Genomic breeding values for RFI were estimated from a reference population of 2,036 individuals in a multi-trait analysis including Australian calf RFI (n=843), Australian lactating cow RFI (n=234), and UK and Dutch lactating cow RFI (n=958). In all cases, the RFI phenotypes were deviations from a mean of 0, calculated by correcting dry matter intake for BW, growth, and milk yield (in the case of lactating cows). Single nucleotide polymorphism effects were calculated from the output of genomic BLUP and used to predict breeding values of 4,106 Holstein sires that were genotyped but did not have RFI phenotypes themselves. These bulls already had BW breeding values calculated from type traits, from which maintenance requirements in kilograms of feed per year were inferred. Finally, RFI and the feed required for maintenance (through BW) were used to calculate a feed saved breeding value and expressed as the predicted amount of feed saved per year. Animals that were 1 standard deviation above the mean were predicted to eat 66 kg dry matter less per year at the same level of milk production. In a data set of genotyped Holstein sires, the mean reliability of the feed saved breeding value was 0.37. For Holsteins that are not genotyped and for breeds other than Holsteins, feed saved is calculated using BW only. From April 2015, feed saved has been included as part of the Australian national selection index, the Balanced Performance Index (BPI). Selection on the BPI is expected to lead to modest gains in feed efficiency.

How old are quantitative trait loci and how widely do they segregate?

The mutations that cause genetic variation in quantitative traits could be old and segregate across many breeds or they could be young and segregate only within one breed. This has implications for our understanding of the evolution of quantitative traits and for genomic prediction to improve livestock. We investigated the age of quantitative trait loci (QTL) for milk production traits identified as segregating in Holstein dairy cattle. We use a multitrait method and found that six of 11 QTL also segregate in Jerseys. Variants identified as Holstein-only QTL were fixed or rare [minor allele frequency (MAF) < 0.05] in Jersey. The age of the QTL mutations appears to vary from perhaps 2000 to 50,000 generations old. The older QTL tend to have high derived allele frequencies and often segregate across both breeds. Holstein-only QTL were often embedded within longer haplotypes, supporting the conclusion that they are typically younger mutations that have occurred more recently than QTL that segregate in both breeds. A reference population for genomic prediction using both Holsteins and Jersey cattle incorrectly predicted a QTL in Jersey cattle when the QTL only segregates in Holsteins. Overcoming this error should help to make genomic prediction more accurate in smaller breeds.

Toward genomic prediction from whole-genome sequence data: impact of sequencing design on genotype imputation and accuracy of predictions.

Genomic prediction from whole-genome sequence data is attractive, as the accuracy of genomic prediction is no longer bounded by extent of linkage disequilibrium between DNA markers and causal mutations affecting the trait, given the causal mutations are in the data set. A cost-effective strategy could be to sequence a small proportion of the population, and impute sequence data to the rest of the reference population. Here, we describe strategies for selecting individuals for sequencing, based on either pedigree relationships or haplotype diversity. Performance of these strategies (number of variants detected and accuracy of imputation) were evaluated in sequence data simulated through a real Belgian Blue cattle pedigree. A strategy (AHAP), which selected a subset of individuals for sequencing that maximized the number of unique haplotypes (from single-nucleotide polymorphism panel data) sequenced gave good performance across a range of variant minor allele frequencies. We then investigated the optimum number of individuals to sequence by fold coverage given a maximum total sequencing effort. At 600 total fold coverage (x 600), the optimum strategy was to sequence 75 individuals at eightfold coverage. Finally, we investigated the accuracy of genomic predictions that could be achieved. The advantage of using imputed sequence data compared with dense SNP array genotypes was highly dependent on the allele frequency spectrum of the causative mutations affecting the trait. When this followed a neutral distribution, the advantage of the imputed sequence data was small; however, when the causal mutations all had low minor allele frequencies, using the sequence data improved the accuracy of genomic prediction by up to 30%.