International cooperation: The pathway for cattle genomics

Comparison of International Dairy Sire Evaluations from Meta-Analysis of National Estimated Breeding Values and Direct Analysis of Individual Animal Performance Records

Validation of simultaneous deregression of cow and bull breeding values and derivation of appropriate weights

Inner Mongolia Cashmere Goats (IMCGs) are famous for its cashmere quality and it's a unique genetic resource in China. Therefore, it is necessary to use genomic selection to improve the accuracy of selection for fleece traits in Inner Mongolia cashmere goats. The aim of this study was to determine the effect of methods (GBLUP, BayesA, BayesB, Bayesian LASSO, Bayesian Ridge Region) and the reference population size on accuracy of genomic selection in IMCGs. This study fully utilizes the pedigree and phenotype records of fleece traits in 2255 individuals, genotype of 50794 SNPs after quality control, and environmental data to perform genomic selection of fleece traits. Then GBLUP and Bayes series methods (BayesA, BayesB, Bayesian LASSO, Bayesian Ridge Region) were used to perform estimates of genetic parameter and genomic breeding value. And the accuracy of genomic estimated breeding value (GEBV) is evaluated using the five-fold cross validation method. And the analysis of variance and multiple comparison methods were used to determine the best method for genomic selection in fleece traits of IMCGs. Further the different reference population sizes (500, 1000, 1500, and 2000) was set. Then the best method was applied to estimate genome breeding values, and evaluate the impact of reference population sizes on the accuracy of genome selection for fleece traits in IMCGs. It was found that the genomic prediction accuracy for each fleece trait in IMCGs by GBLUP method is highest, and it is significantly higher than that obtained by Bayesian method. The accuracy of breeding value estimation is 58.52% -68.49%. Also, it was found that the size of the reference population has a significant impact on the accuracy of genome prediction of fleece traits. When the reference population size is 2000, the accuracy of genomic prediction for each fleece trait is significantly higher than other levels, with accuracy of 55.47% -67.87%. This provides a theoretical basis for design a reasonable genome selection plan for Inner Mongolia cashmere goats in the later stag.

In recent years, breeders have combined the use of phenotypic appraisal and the estimation of breeding values (EPD) to make genetic selection decisions in beef and dairy cattle that have resulted in a steady genetic gain of 2% per year. However, the most extensive application of genomics has occurred in the dairy industry with the estimation of molecular breeding values that has improved selection efficiency to a much higher order of magnitude. Despite a growing molecular and physiological understanding of complex traits, little is known about the genes determining the traits and their precise function, and a significant unexplained source of variation of phenotypes remains in livestock. Within this context, a more complete understanding of the genes and regulatory pathways and networks involved in economically important traits (i.e. fertility and reproduction, feed efficiency, meat quality and carcass traits) in beef cattle will provide knowledge to help improve genetic selection and reproductive management. Currently, with all the new available technologies in livestock combined with statistical methodologies, the integration of structural and functional genomics information with other –OMICS into a systems biology approach has allowed development of a better biological understanding of phenotypes complementing the traditional genetic tools and further advance identification of functional genes. As part of the genomics tool box and the HD-genotyping SNP chips, whole genome sequencing technologies are now available in cattle and extensively utilized in genetic improvement. As a part of high throughput tools available for genomic analysis, RNA-Sequencing allows measuring not only gene expression, but also examining genome structure identifying SNP and other structural variation such as insertions, deletions and splice variants. The expectation is that the integration of all these types of genomic data will accelerate the genetic improvement by improving accuracy of selection and reducing the generation interval. Combining the information from the –OMICS technologies (i.e., transcriptomics, metagenomics, metabolomics, amongst several others) together with metabolic pathways and functional/biological analysis into a systems biology approach allows the identification of functional SNP increasing the accuracy of selection. The particular benefits of new integrated high-throughput genomics technologies within a systems biology approach will most likely be used to accelerate the genetic improvement of those traits that are difficult to measure such as health, feed efficiency, methane emission and fertility and reproduction traits in beef cattle.

The main goal in animal breeding is to select individuals that have high breeding values for traits of interest as parents to produce the next generation and to do so as quickly as possible. To date, most programs rely on statistical analysis of large data bases with phenotypes on breeding populations by linear mixed model methodology to estimate breeding values on selection candidates. However, there is a long history of research on the use of genetic markers to identify quantitative trait loci and their use in marker-assisted selection but with limited implementation in practical breeding programs. The advent of high-density SNP genotyping, combined with novel statistical methods for the use of this data to estimate breeding values, has resulted in the recent extensive application of genomic or whole-genome selection in dairy cattle and research to implement genomic selection in other livestock species is underway. The high-density SNP data also provides opportunities to detect QTL and to encover the genetic architecture of quantitative traits, in terms of the distribution of the size of genetic effects that contribute to trait differences in a population. Results show that this genetic architecture differs between traits but that for most traits, over 50% of the genetic variation resides in genomic regions with small effects that are of the order of magnitude that is expected under a highly polygenic model of inheritance.

Reliabilities for genomic estimated breeding values (GEBV) were investigated by simulation for a typical dairy cattle breeding setting. Scenarios were simulated with different heritabilites (h2) and for different haplotype sizes, and seven generations with only genotypes were generated to investigate reliability of GEBV over time. A genome with 5000 single nucleotide polymorphisms (SNP) at distances of 0.1 cM and 50 quantitative trait loci (QTL) was simulated, and a Bayesian variable selection model was implemented to predict GEBV. Highest reliabilities were obtained for 10 SNP haplotypes. At optimal haplotype size, reliabilities in generation 1 without phenotypes ranged from 0.80 for h2 = 0.02 to 0.93 for h2 = 0.30, and in the seventh generation without phenotypes ranged from 0.69 for h2 = 0.02 to 0.86 for h2 = 0.30. Reliabilities of GEBV were found sufficiently high to implement dairy selection schemes without progeny testing in which case a data time-lag of two to three generations may be present. Reliabilities were also relatively high for low heritable traits, implying that genomic selection could be especially beneficial to improve the selection on, e.g. health and fertility.

BackgroundOver the last ten years, genomic selection has developed enormously. Simulations and results on real data suggest that breeding values can be predicted with high accuracy using genetic markers alone. However, to reach high accuracies, large reference populations are needed. In many livestock populations or even species, such populations cannot be established when traits are difficult or expensive to record, or when the population size is small. The value of genomic selection is then questionable.MethodsIn this study, we compare traditional breeding schemes based on own performance or progeny information to genomic selection schemes, for which the number of phenotypic records is limiting. Deterministic simulations were performed using selection index theory. Our focus was on the equilibrium response obtained after a few generations of selection. Therefore, we first investigated the magnitude of the Bulmer effect with genomic selection.ResultsResults showed that the reduction in response due to the Bulmer effect is the same for genomic selection as for selection based on traditional BLUP estimated breeding values, and is independent of the accuracy of selection. The reduction in response with genomic selection is greater than with selection based directly on phenotypes without the use of pedigree information, such as mass selection. To maximize the accuracy of genomic estimated breeding values when the number of phenotypic records is limiting, the same individuals should be phenotyped and genotyped, rather than genotyping parents and phenotyping their progeny. When the generation interval cannot be reduced with genomic selection, large reference populations are required to obtain a similar response to that with selection based on BLUP estimated breeding values based on own performance or progeny information. However, when a genomic selection scheme has a moderate decrease in generation interval, relatively small reference population sizes are needed to obtain a similar response to that with selection on traditional BLUP estimated breeding values.ConclusionsWhen the trait of interest cannot be recorded on the selection candidate, genomic selection schemes are very attractive even when the number of phenotypic records is limited, because traditional breeding requires progeny testing schemes with long generation intervals in those cases.

ObjectiveWhen evaluating individuals with the same parent and no phenotype by pedigree best linear unbiased prediction (BLUP), it is difficult to explain carcass grade difference and select individuals because they have the same value in pedigree BLUP (PBLUP). However, single step GBLUP (ssGBLUP), which can estimate the breeding value suitable for the individual by adding genotype, is more accurate than the existing method.MethodsThe breeding value and accuracy were estimated with pedigree BLUP and ssGBLUP using pedigree and genotype of 408 Hanwoo cattle from 16 families with the same parent among siblings produced by fertilized egg transplantation. A total of 14,225 Hanwoo cattle with pedigree, genotype and phenotype were used as the reference population. PBLUP obtained estimated breeding value (EBV) using the pedigree of the test and reference populations, and ssGBLUP obtained genomic EBV (GEBV) after constructing and H-matrix by integrating the pedigree and genotype of the test and reference populations.ResultsFor all traits, the accuracy of GEBV using ssGBLUP is 0.18 to 0.20 higher than the accuracy of EBV obtained with PBLUP. Comparison of EBV and GEBV of individuals without phenotype, since the value of EBV is estimated based on expected values of alleles passed down from common ancestors. It does not take Mendelian sampling into consideration, so the EBV of all individuals within the same family is estimated to be the same value. However, GEBV makes estimating true kinship coefficient based on different genotypes of individuals possible, so GEBV that corresponds to each individual is estimated rather than a uniform GEBV for each individual.ConclusionSince Hanwoo cows bred through embryo transfer have a high possibility of having the same parent, if ssGBLUP after adding genotype is used, estimating true kinship coefficient corresponding to each individual becomes possible, allowing for more accurate estimation of breeding value.

BackgroundThe markers detected by genome-wide association study (GWAS) make it possible to dissect genetic structure and diversity at many loci. This can enable a wheat breeder to reveal and used genomic loci controlling drought tolerance. This study was focused on determining the population structure of Iranian 208 wheat landraces and 90 cultivars via genotyping-by-sequencing (GBS) and also on detecting marker-trait associations (MTAs) by GWAS and genomic prediction (GS) of wheat agronomic traits for drought-tolerance breeding. GWASs were conducted using both the original phenotypes (pGWAS) and estimated breeding values (eGWAS). The bayesian ridge regression (BRR), genomic best linear unbiased prediction (gBLUP), and ridge regression-best linear unbiased prediction (rrBLUP) approaches were used to estimate breeding values and estimate prediction accuracies in genomic selection.ResultsPopulation structure analysis using 2,174,975 SNPs revealed four genetically distinct sub-populations from wheat accessions. D-Genome harbored the lowest number of significant marker pairs and the highest linkage disequilibrium (LD), reflecting different evolutionary histories of wheat genomes. From pGWAS, BRR, gBLUP, and rrBLUP, 284, 363, 359 and 295 significant MTAs were found under normal and 195, 365, 362 and 302 under stress conditions, respectively. The gBLUP with the most similarity (80.98 and 71.28% in well-watered and rain-fed environments, correspondingly) with the pGWAS method in the terms of discovered significant SNPs, suggesting the potential of gBLUP in uncovering SNPs. Results from gene ontology revealed that 29 and 30 SNPs in the imputed dataset were located in protein-coding regions for well-watered and rain-fed conditions, respectively. gBLUP model revealed genetic effects better than other models, suggesting a suitable tool for genome selection in wheat.ConclusionWe illustrate that Iranian landraces of bread wheat contain novel alleles that are adaptive to drought stress environments. gBLUP model can be helpful for fine mapping and cloning of the relevant QTLs and genes, and for carrying out trait introgression and marker-assisted selection in both normal and drought environments in wheat collections.

This study compared genomic predictions using conventional estimated breeding values (EBV) and daughter yield deviations (DYD) as response variables based on simulated data. Eight scenarios were simulated in regard to heritability (0.05 and 0.30), number of daughters per sire (30, 100, and unequal numbers with an average of 100 per sire) and numbers of genotyped sires (all or half of sires were genotyped). The simulated genome had a length of 1200 cM with 15,000 equally spaced Single-nucleotide polymorphism (SNP) markers and 500 randomly distributed Quantitative trait locus (QTL). In the simulated scenarios, the EBV approach was as effective as or slightly better than the DYD approach at predicting breeding value, dependent on simulated scenarios and statistical models. Applying a Bayesian common prior model (the same prior distribution of marker effect variance) and a linear mixed model (GBLUP), the EBV and DYD approaches provided similar genomic estimated breeding value (GEBV) reliabilities, except for scenarios with unequal numbers of daughters and half of sires without genotype, for which the EBV approach was superior to the DYD approach (by 1.2 and 2.4%). Using a Bayesian mixture prior model (mixture prior distribution of marker effect variance), the EBV approach resulted in slightly higher reliabilities of GEBV than the DYD approach (by 0.3-3.6% with an average of 1.9%), and more obvious in scenarios with low heritability, small or unequal numbers of daughters, and half of sires without genotype. Moreover, the results showed that the correlation between GEBV and conventional parent average (PA) was lower (corresponding to a relatively larger gain by including PA) when using the DYD approach than when using the EBV approach. Consequently, the two approaches led to similar reliability of an index combining GEBV and PA in most scenarios. These results indicate that EBV can be used as an alternative response variable for genomic prediction.

Estimated breeding values (EBVs) are traditionally obtained from pedigree information. However, EBVs from high-density genotypes can have higher accuracy than EBVs from pedigree information. At the same time, it has been shown that EBVs from genomic data lead to lower increases in inbreeding compared with traditional selection based on genealogies. Here we evaluate the performance with BLUP selection based on genealogical coancestry with three different genome-based coancestry estimates: (1) an estimate based on shared segments of homozygosity, (2) an approach based on SNP-by-SNP count corrected by allelic frequencies, and (3) the identity by state methodology. We evaluate the effect of different population sizes, different number of genomic markers, and several heritability values for a quantitative trait. The performance of the different measures of coancestry in BLUP is evaluated in the true breeding values after truncation selection and also in terms of coancestry and diversity maintained. Accordingly, cross-performances were also carried out, that is, how prediction based on genealogical records impacts the three other measures of coancestry and inbreeding, and viceversa. Our results show that the genetic gains are very similar for all four coancestries, but the genomic-based methods are superior to using genealogical coancestries in terms of maintaining diversity measured as observed heterozygosity. Furthermore, the measure of coancestry based on shared segments of the genome seems to provide slightly better results on some scenarios, and the increase in inbreeding and loss in diversity is only slightly larger than the other genomic selection methods in those scenarios. Our results shed light on genomic selection vs. traditional genealogical-based BLUP and make the case to manage the population variability using genomic information to preserve the future success of selection programmes.

Genomic selection (GS) uses associations between markers and phenotypes to predict the breeding values of individuals. It can be applied early in the breeding cycle to reduce the cross-to-cross generation interval and thereby increase genetic gain per unit of time. The development of cost-effective, high-throughput genotyping platforms has revolutionized plant breeding programs by enabling the implementation of GS at the scale required to achieve impact. As a result, GS is becoming routine in plant breeding, even in minor crops such as pulses. Here we examined 2,081 breeding lines from Agriculture Victoria's national lentil breeding program for a range of target traits including grain yield, ascochyta blight resistance, botrytis grey mould resistance, salinity and boron stress tolerance, 100-grain weight, seed size index and protein content. A broad range of narrow-sense heritabilities was observed across these traits (0.24-0.66). Genomic prediction models were developed based on 64,781 genome-wide SNPs using Bayesian methodology and genomic estimated breeding values (GEBVs) were calculated. Forward cross-validation was applied to examine the prediction accuracy of GS for these targeted traits. The accuracy of GEBVs was consistently higher (0.34-0.83) than BLUP estimated breeding values (EBVs) (0.22-0.54), indicating a higher expected rate of genetic gain with GS. GS-led parental selection using early generation breeding materials also resulted in higher genetic gain compared to BLUP-based selection performed using later generation breeding lines. Our results show that implementing GS in lentil breeding will fast track the development of high-yielding cultivars with increased resistance to biotic and abiotic stresses, as well as improved seed quality traits.

Complex traits in sugarcane can be accurately predicted using genome-wide DNA markers. Genomic single-step prediction is an attractive method for genomic selection in commercial breeding programs. Sugarcane breeding programs have achieved up to 1% genetic gain in key traits such as tonnes of cane per hectare (TCH), commercial cane sugar (CCS) and Fibre content over the past decades. Here, we assess the potential of genomic selection to increase the rate of genetic gain for these traits by deriving genomic estimated breeding values (GEBVs) from a reference population of 3984 clones genotyped for 26K SNP. We evaluated the three different genomic prediction approaches GBLUP, genomic single step (GenomicSS), and BayesR. GenomicSS combining pedigree and SNP information from historic and recent breeding programs achieved the most accurate predictions for most traits (0.3-0.44). This method is attractive for routine genetic evaluation because it requires relatively little modification to the existing evaluation and results in breeding value estimates for all individuals, not only those genotyped. Adding information from early-stage trials added up to 5% accuracy for CCS and Fibre, but 0% for TCH, reflecting the importance of competition effects for TCH. These GEBV accuracies are sufficiently high that, combined with the right breeding strategy, a doubling of the rate of genetic gain could be achieved. We also assessed the flowering traits days to flowering, gender and pollen viability and found high heritabilities of 0.57, 0.78 and 0.72, respectively. The GEBV accuracies indicated that genomic selection could be used to improve these traits. This could open new avenues for breeders to manage their breeding programs, for example, by synchronising flowering time and selecting males with high pollen viability.

Genomic selection refers to the use of dense markers covering the whole genome to estimate the breeding value of selection candidates for a quantitative trait. This paper considers prediction of breeding value based on a linear combination of the markers. In this case the best estimate of each marker's effect is the expectation of the effect conditional on the data. To calculate this requires a prior distribution of marker effects. If the marker effects are normally distributed with constant variance, BLUP can be used to calculate the estimated effects of the markers and hence the estimated breeding value (EBV). In this case the model is equivalent to a conventional animal model in which the relationship matrix among the animals is estimated from the markers instead of the pedigree. The accuracy of the EBV can approach 1.0 but a very large amount of data is required. An alternative model was investigated in which only some markers have non-zero effects and these effects follow a reflected exponential distribution. In this case the expected effect of a marker is a non-linear function of the data such that apparently small effects are regressed back almost to zero and consequently these markers can be deleted from the model. The accuracy in this case is considerably higher than when marker effects are normally distributed. If genomic selection is practiced for several generations the response declines in a manner that can be predicted from the marker allele frequencies. Genomic selection is likely to lead to a more rapid decline in the selection response than phenotypic selection unless new markers are continually added to the prediction of breeding value. A method to find the optimum index to maximise long term selection response is derived. This index varies the weight given to a marker according to its frequency such that markers where the favourable allele has low frequency receive more weight in the index.

Genomic evaluation by single-step procedure allows efficient implanting of genetic SNP markers into complicated procedure of random regression test-day-model of milk traits. Test-day records and pedigree files on nation-wide scale are combined with genomic relationship into common evaluation of both genotyped and ungenotyped animals. Due to strong import of foreign sperm into small national populations is the reliability of evaluation of young animals low. This is particularly seen in evaluation of young bulls, which frequently have both parents foreign. Genomic evaluation helps and notably improves reliability of evaluation. ssGBLUP procedure is advantageous especially for small populations. Domestic genomic evaluation of young animals has medium to high correlation with foreign Interbull values. Interbull conversion of values of bulls according MACE, which works with progeny tested bulls, is more reliable than conversion according GMACE procedure, which works with genomic evaluation of young animals. Keywords: genomic breeding value, ssGBLUP, test day model, MACE, GMACE References Bauer, J. et al. (2014) Approximation of the reliability of single-step genomic breeding values for dairy cattle in the Czech Republic. Anim. Sci. Papers and Reports, 32, pp. 301-306. Bauer, J., Přibyl, J. and Vostrý, L. (2015) Contribution of domestic production records and Interbull EBV on approximate reliabilities of single-step genomic breeding values in dairy cattle. Czech J. Anim. Sci., 60, 263-267. Christensen, O.F. and Lund, M.S. (2010) Genomic prediction when some animals are not genotyped. Genet.Sel.Evol. 42, pp. 2. Forni, S., Aguilar I. and Misztal, I. (2011) Different genomic relationship matrices for single step analysis using phenotypic, pedigree and genomic information. Genet. Sel. Evol., 43, pp. 1. Legarra A., Aguilar I. and Misztal, I. (2009) A relationship matrix including full pedigree and genomic information. J. Dairy Sci., 92, pp. 4656-4663. Meuwissen, T.H.E., Hayes, B.J. and Goddard, M.E. (2001) Prediction of total genetic value using genome-wide dense marker maps. Genetics, 157, pp. 1819–1829. Misztal, I., Legarra A. and Aguilar, I. (2009) Computing procedures for genetic evaluation including phenotypic, full pedigree, and genomic information. J. Dairy Sci., 92, pp. 4648–4655. Pesek, P., Přibyl, J. and Vostrý, L . (2015) Genetic variances of SNP loci for milk yield in dairy cattle. J. Appl. Genet., 56, pp. 339-347. Plemdat, (2015) Descriptions of Breeding values Evaluation. Retrieved on 10 th June 2015 From www.plemdat.cz. Přibyl, J. et al. (2014) Domestic and Interbull information in the single step genomic evaluation of Holstein milk production. Czech J. Anim. Sci., 59, pp. 409-415. Přibyl, J. et al. (2015) Domestic estimated breeding values and genomic enhanced breeding values of bulls in comparison with their foreign genomic enhanced breeding values. Animal, 9, pp. 1635-1642. Vitezica, Z.G. et al. (2011) Bias in genomic predictions for populations under selection. Genet. Res. (Camb), 93, pp. 357-366. Zavadilova, L., Jamrozik, J. and Schaeffer, L.R. (2005a) Genetic parameters for test- day model with random regressions for production traits of Czech Holstein cattle. Czech J. Anim. Sci., 50, pp. 142-154. Zavadilova, L., Němcova, E. and Wolf, J. (2005b) Definition of subgroups for fixed regression in the test-day animal model for milk production of Holstein cattle in the Czech Republic. Czech J. Anim. Sci., 50, pp. 7-13. Zavadilova, L. et al. (2014) Single-step genomic evaluation for linear type traits of Holstein cows in Czech Republic. Anim. Sci. Papers and Reports vol. 32, pp. 201-208.