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.

Gene expression of the heat stress response in bovine peripheral white blood cells and milk somatic cells in vivo.

Heat stress in dairy cattle leads to reduction in feed intake and milk production as well as the induction of many physiological stress responses. The genes implicated in the response to heat stress in vivo are not well characterised. With the aim of identifying such genes, an experiment was conducted to perform differential gene expression in peripheral white blood cells and milk somatic cells in vivo in 6 Holstein Friesian cows in thermoneutral conditions and in 6 Holstein Friesian cows exposed to a short-term moderate heat challenge. RNA sequences from peripheral white blood cells and milk somatic cells were used to quantify full transcriptome gene expression. Genes commonly differentially expressed (DE) in both the peripheral white blood cells and in milk somatic cells were associated with the cellular stress response, apoptosis, oxidative stress and glucose metabolism. Genes DE in peripheral white blood cells of cows exposed to the heat challenge compared to the thermoneutral control were related to inflammation, lipid metabolism, carbohydrate metabolism and the cardiovascular system. Genes DE in milk somatic cells compared to the thermoneutral control were involved in the response to stress, thermoregulation and vasodilation. These findings provide new insights into the cellular adaptations induced during the response to short term moderate heat stress in dairy cattle and identify potential candidate genes (BDKRB1 and SNORA19) for future research.

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.

Sequencing the mosaic genome of Brahman cattle identifies historic and recent introgression including polled.

Brahman cattle have a Bos indicus and Bos taurus mosaic genome, as a result of the process used to create the breed (repeat backcrossing of Bos taurus females to Bos indicus bulls). With the aim of identifying Bos taurus segments in the Brahman genome at sequence level resolution, we sequenced the genomes of 46 influential Brahman bulls. Using 36 million variants identified in the sequences, we searched for regions close to fixation for Bos indicus or Bos taurus segments that were longer than expected by chance (from simulation of the breed formation history of Brahman cattle). Regions close to fixation for Bos indicus content were enriched for protein synthesis genes, while regions of higher Bos taurus content included genes of the G-protein coupled receptor family (including genes implicated in puberty, such as THRS). The region with the most extreme Bos taurus enrichment was on chromosome 14 surrounding PLAG1. The introgressed Bos taurus allele at PLAG1 increases stature and the high frequency of the allele likely reflects strong selection for the trait. Finally, we provide evidence that the polled mutation in Brahmans, a desirable trait under very strong recent selection, is of Celtic origin and is introgressed from Bos taurus.

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.

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.

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.

Rapid Discovery of De Novo Deleterious Mutations in Cattle Enhances the Value of Livestock as Model Species.

In humans, the clinical and molecular characterization of sporadic syndromes is often hindered by the small number of patients and the difficulty in developing animal models for severe dominant conditions. Here we show that the availability of large data sets of whole-genome sequences, high-density SNP chip genotypes and extensive recording of phenotype offers an unprecedented opportunity to quickly dissect the genetic architecture of severe dominant conditions in livestock. We report on the identification of seven dominant de novo mutations in CHD7, COL1A1, COL2A1, COPA, and MITF and exploit the structure of cattle populations to describe their clinical consequences and map modifier loci. Moreover, we demonstrate that the emergence of recessive genetic defects can be monitored by detecting de novo deleterious mutations in the genome of bulls used for artificial insemination. These results demonstrate the attractiveness of cattle as a model species in the post genomic era, particularly to confirm the genetic aetiology of isolated clinical case reports in humans.

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.

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.

Genomewide association study of methane emissions in Angus beef cattle with validation in dairy cattle.

Methane (CH) is a product of enteric fermentation in ruminants, and it represents around 17% of global CH emissions. There has been substantial effort from the livestock scientific community toward tools that can help reduce this percentage. One approach is to select for lower emitting animals. To achieve this, accurate genetic parameters and identification of the genomic basis of CH traits are required. Therefore, the objectives of this study were 1) to perform a genomewide association study to identify SNP associated with several CH traits in Angus beef cattle (1,020 animals) and validate them in a lactating Holstein population (population 1 [POP1]; 205 animals); 2) to validate significant SNP for DMI and weight at test (WT) from a second Holstein population, from a previous study (population 2 [POP2]; 903 animals), in an Angus population; and 3) to evaluate 2 different residual CH traits and determine if the genes associated with CH also control residual CH traits. Phenotypes calculated for the genotyped Angus population included CH production (MeP), CH yield (MeY), CH intensity (MI), DMI, and WT. The Holstein population (POP1) was multiparous, with phenotypes on CH traits (MeP, MeY, and MI) plus genotypes. Additionally, 2 CH traits, residual genetic CH (RGM) and residual phenotypic CH (RPM), were calculated by adjusting MeP for DMI and WT. Estimated heritabilities in the Angus population were 0.30, 0.19, and 0.15 for MeP, RGM, and RPM, respectively, and genetic correlations of MeP with DMI and WT were 0.83 and 0.80, respectively. Estimated heritabilities in Holstein POP1 were 0.23, 0.30, and 0.42 for MeP, MeY, and MI, respectively. Strong associations with MeP were found on chromosomes 4, 12, 14, 20, and 30 at < 0.001, and those chromosomes also had significant SNP for DMI in Holstein POP1. In the Angus population, the number of significant SNP for MeP at < 0.005 was 3,304, and approximately 630 of those SNP also were important for DMI and WT. When a set (approximately 3,300) of significant SNP for DMI and WT in the Angus population was used to estimate genetic parameters for MeP and MeY in Holstein POP1, the genetic variance and, consequently, the heritability slightly increased, meaning that most of the genetic variation is largely captured by these SNP. Residual traits could be a good option to include in the breeding goal, as this would facilitate selection for lower emitting animals without compromising DMI and WT.

Genomic Selection Improves Heat Tolerance in Dairy Cattle.

Dairy products are a key source of valuable proteins and fats for many millions of people worldwide. Dairy cattle are highly susceptible to heat-stress induced decline in milk production, and as the frequency and duration of heat-stress events increases, the long term security of nutrition from dairy products is threatened. Identification of dairy cattle more tolerant of heat stress conditions would be an important progression towards breeding better adapted dairy herds to future climates. Breeding for heat tolerance could be accelerated with genomic selection, using genome wide DNA markers that predict tolerance to heat stress. Here we demonstrate the value of genomic predictions for heat tolerance in cohorts of Holstein cows predicted to be heat tolerant and heat susceptible using controlled-climate chambers simulating a moderate heatwave event. Not only was the heat challenge stimulated decline in milk production less in cows genomically predicted to be heat-tolerant, physiological indicators such as rectal and intra-vaginal temperatures had reduced increases over the 4 day heat challenge. This demonstrates that genomic selection for heat tolerance in dairy cattle is a step towards securing a valuable source of nutrition and improving animal welfare facing a future with predicted increases in heat stress events.

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.

Genomic heritabilities and genomic estimated breeding values for methane traits in Angus cattle.

Enteric methane emissions from beef cattle are a significant component of total greenhouse gas emissions from agriculture. The variation between beef cattle in methane emissions is partly genetic, whether measured as methane production, methane yield (methane production/DMI), or residual methane production (observed methane production - expected methane production), with heritabilities ranging from 0.19 to 0.29. This suggests methane emissions could be reduced by selection. Given the high cost of measuring methane production from individual beef cattle, genomic selection is the most feasible approach to achieve this reduction in emissions. We derived genomic EBV (GEBV) for methane traits from a reference set of 747 Angus animals phenotyped for methane traits and genotyped for 630,000 SNP. The accuracy of GEBV was tested in a validation set of 273 Angus animals phenotyped for the same traits. Accuracies of GEBV ranged from 0.29 ± 0.06 for methane yield and 0.35 ± 0.06 for residual methane production. Selection on GEBV using the genomic prediction equations derived here could reduce emissions for Angus cattle by roughly 5% over 10 yr.

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.