Estimated Breeding Values Research Articles

Haplotype-tagged SNPs improve genomic prediction accuracy for Fusarium head blight resistance and yield-related traits in wheat

Key messageLinkage disequilibrium (LD)-based haplotyping with subsequent SNP tagging improved the genomic prediction accuracy up to 0.07 and 0.092 for Fusarium head blight resistance and spike width, respectively, across six different models.Genomic prediction is a powerful tool to enhance genetic gain in plant breeding. However, the method is accompanied by various complications leading to low prediction accuracy. One of the major challenges arises from the complex dimensionality of marker data. To overcome this issue, we applied two pre-selection methods for SNP markers viz. LD-based haplotype-tagging and GWAS-based trait-linked marker identification. Six different models were tested with preselected SNPs to predict the genomic estimated breeding values (GEBVs) of four traits measured in 419 winter wheat genotypes. Ten different sets of haplotype-tagged SNPs were selected by adjusting the level of LD thresholds. In addition, various sets of trait-linked SNPs were identified with different scenarios from the training-test combined and only from the training populations. The BRR and RR-BLUP models developed from haplotype-tagged SNPs had a higher prediction accuracy for FHB and SPW by 0.07 and 0.092, respectively, compared to the corresponding models developed without marker pre-selection. The highest prediction accuracy for SPW and FHB was achieved with tagged SNPs pruned at weak LD thresholds (r2 < 0.5), while stringent LD was required for spike length (SPL) and flag leaf area (FLA). Trait-linked SNPs identified only from training populations failed to improve the prediction accuracy of the four studied traits. Pre-selection of SNPs via LD-based haplotype-tagging could play a vital role in optimizing genomic selection and reducing genotyping costs. Furthermore, the method could pave the way for developing low-cost genotyping methods through customized genotyping platforms targeting key SNP markers tagged to essential haplotype blocks.

Effects of genomic estimated breeding value and dietary energy to protein ratio on growth performance, carcass trait, and lipogenic gene expression in Hanwoo steer

"Genome-based precision feeding" is a concept that involves the application of customised diets to different genetic groups of cattle. We investigated the effects of the genomic estimated breeding value (gEBV) and dietary energy to protein ratio (DEP) on growth performance, carcass traits, and lipogenic gene expression in Hanwoo (Korean cattle) steers. Forty-four Hanwoo steers (BW = 636 kg, age = 26.9 months) were genotyped using the Illumina Bovine 50 K BeadChip. The gEBV was calculated using genomic best linear unbiased prediction. Animals were separated into high gEBV of marbling score or low-gMS groups based on the upper and lower 50% groupings of the reference population, respectively. Animals were assigned to one of four groups in a 2 × 2 factorial arrangement: high gMS/high DEP (0.084 MJ/g), high gMS/low DEP (0.079 MJ/g), low gMS/high DEP, and low gMS/low DEP. Steers were fed concentrate with a high or low DEP for 31 weeks. The BW tended to be higher (0.05 < P < 0.1) in the high-gMS groups compared to the low-gMS groups at 0, 4, 8, 12, and 20 weeks. The average daily gain (ADG) tended to be lower (P = 0.08) in the high-gMS group than in the low-gMS group. Final BW and measured carcass weight (CW) were positively correlated with the gEBV of carcass weight (gCW). The DEP did not affect ADG. Neither the gMS nor the DEP affected the MS and beef quality grade. The intramuscular fat (IMF) content in the longissimus thoracis (LT) tended to be higher (P = 0.08) in the high-gMS groups than in the low-gMS groups. The mRNA levels of lipogenic acetyl-CoA carboxylase and fatty acid binding protein 4 genes in the LT were higher (P < 0.05) in the high-gMS group than in the low-gMS group. Overall, the IMF content tended to be affected by the gMS, and the genetic potential (i.e., gMS) was associated with the functional activity of lipogenic gene expression. The gCW was associated with the measured BW and CW. The results demonstrated that the gMS and the gCW may be used as early prediction indexes for meat quality and growth potential of beef cattle.

Modelling growth in Suffolk and Charollais sheep populations using random regression models and validation of constrained polynomial correlation values

Random regression modelling has been used across multiple animal species to model longitudinal data. The random regression model for growth accounts for the genetic correlation between measures of the same trait over time and the wide environmental variability in growth, but this requires adequate weight records across the age range. However, contemporary management practices in sheep in the United Kingdom generally focus on growing lambs and neglect mature weight recordings. This study examined modelling strategies for growth data in Suffolk and Charollais sheep, provided by the Agriculture and Horticulture Development Board, with polynomial random regression modelling with many early life weight recordings but limited weight recordings in mature animals. Two methods were employed to model the data. In Method A, missing mature weight records were predicted for those animals that did not have a recorded mature weight. The animals were sorted into groups based on the identity of their sires and the year in which the animal was born. Mature weight values were predicted within each group with a multiple regression model. The dataset, including predicted values, was analysed with random regression models using polynomials and simple linear regression for animal and permanent environmental (PE) effects. In Method B, the dataset with missing mature weight records was analysed using a random linear regression animal model with random animal and PE effects. Due to problems of convergence because the parameters were close to the boundary space, fixing the correlation between the intercept and slope of the Legendre polynomial at different levels was investigated. The heritability values resulting from the model with a fixed correlation between intercept and slope parameters at 0.5 for the Suffolk dataset resulted in heritability values ranging from 0.2 to 0.5 from 1 to 619 days of age. Corresponding estimates for the Charollais dataset ranged from 0.18 to 0.49 from 1 to 640 days of age. For the Suffolk data, the genetic correlations ranged from 1.00 to 0.08 between weight at day 1 to weight at day 619, while for the Charollais, the correlations ranged from 1.00 to 0.05 from 1 to 640 days of age. Validation procedures were undertaken using a multitrait approach to examine the estimated breeding values when the correlation between the intercept and slope are fixed at different levels. The results indicated that fixing the correlation at 0.5 gave the most appropriate estimates for the Suffolk and Charollais datasets.

Genomics combined with UAS data enhances prediction of grain yield in winter wheat.

With the human population continuing to increase worldwide, there is pressure to employ novel technologies to increase genetic gain in plant breeding programs that contribute to nutrition and food security. Genomic selection (GS) has the potential to increase genetic gain because it can accelerate the breeding cycle, increase the accuracy of estimated breeding values, and improve selection accuracy. However, with recent advances in high throughput phenotyping in plant breeding programs, the opportunity to integrate genomic and phenotypic data to increase prediction accuracy is present. In this paper, we applied GS to winter wheat data integrating two types of inputs: genomic and phenotypic. We observed the best accuracy of grain yield when combining both genomic and phenotypic inputs, while only using genomic information fared poorly. In general, the predictions with only phenotypic information were very competitive to using both sources of information, and in many cases using only phenotypic information provided the best accuracy. Our results are encouraging because it is clear we can enhance the prediction accuracy of GS by integrating high quality phenotypic inputs in the models.

Genetic Parameters for a Weighted Analysis of Survivability in Dairy Cattle.

The genetic parameters for the survival of Holstein cows, analysed in nine consecutive time periods during the first three calving intervals, were estimated. The earlier the animals are culled, the more they are informationally underestimated. This undervaluing can be remedied by using a weighted analysis that balances the amount of information. If the method of estimating breeding values changes, the genetic parameters will also change. The Holstein cattle dataset from 2005 to 2017 used in this study included 1,813,636 survival records from 298,290 cows. The pedigree with three generations of ancestors included 660,476 individuals. Linear repeatability models estimated genetic parameters for overall and functional survivability. Due to weights, heritability increased from 0.013 to 0.057. Repeatability with weights was 0.505. The standard deviations of breeding values were 1.75 and 2.18 without weights and 6.04 and 6.20 with weights. Including weights in the calculation increased the additive variance proportion and the breeding values' reliabilities. We conclude that the main contribution of the weighted method we have presented is to compensate for the lack of records in culled individuals with a positive impact on the reliability of the breeding value.

Genomic selection strategies for clonally propagated crops

Key messageFor genomic selection in clonally propagated crops with diploid (-like) meiotic behavior to be effective, crossing parents should be selected based on genomic predicted cross-performance unless dominance is negligible.For genomic selection (GS) in clonal breeding programs to be effective, parents should be selected based on genomic predicted cross-performance unless dominance is negligible. Genomic prediction of cross-performance enables efficient exploitation of the additive and dominance value simultaneously. Here, we compared different GS strategies for clonally propagated crops with diploid (-like) meiotic behavior, using strawberry as an example. We used stochastic simulation to evaluate six combinations of three breeding programs and two parent selection methods. The three breeding programs included (1) a breeding program that introduced GS in the first clonal stage, and (2) two variations of a two-part breeding program with one and three crossing cycles per year, respectively. The two parent selection methods were (1) parent selection based on genomic estimated breeding values (GEBVs) and (2) parent selection based on genomic predicted cross-performance (GPCP). Selection of parents based on GPCP produced faster genetic gain than selection of parents based on GEBVs because it reduced inbreeding when the dominance degree increased. The two-part breeding programs with one and three crossing cycles per year using GPCP always produced the most genetic gain unless dominance was negligible. We conclude that (1) in clonal breeding programs with GS, parents should be selected based on GPCP, and (2) a two-part breeding program with parent selection based on GPCP to rapidly drive population improvement has great potential to improve breeding clonally propagated crops.

Increasing genomic prediction accuracy for unphenotyped full-sib families by modeling additive and dominance effects with large datasets in white spruce.

Genomic selection is becoming a standard technique in plant breeding and is now being introduced into forest tree breeding. Despite promising results to predict the genetic merit of superior material based on their additive breeding values, many studies and operational programs still neglect non-additive effects and their potential for enhancing genetic gains. Using two large comprehensive datasets totaling 4,066 trees from 146 full-sib families of white spruce (Picea glauca (Moench) Voss), we evaluated the effect of the inclusion of dominance on the precision of genetic parameter estimates and on the accuracy of conventional pedigree-based (ABLUP-AD) and genomic-based (GBLUP-AD) models. While wood quality traits were mostly additively inherited, considerable non-additive effects and lower heritabilities were detected for growth traits. For growth, GBLUP-AD better partitioned the additive and dominance effects into roughly equal variances, while ABLUP-AD strongly overestimated dominance. The predictive abilities of breeding and total genetic value estimates were similar between ABLUP-AD and GBLUP-AD when predicting individuals from the same families as those included in the training dataset. However, GBLUP-AD outperformed ABLUP-AD when predicting for new unphenotyped families that were not represented in the training dataset, with, on average, 22% and 53% higher predictive ability of breeding and genetic values, respectively. Resampling simulations showed that GBLUP-AD required smaller sample sizes than ABLUP-AD to produce precise estimates of genetic variances and accurate predictions of genetic values. Still, regardless of the method used, large training datasets were needed to estimate additive and non-additive genetic variances precisely. This study highlights the different quantitative genetic architectures between growth and wood traits. Furthermore, the usefulness of genomic additive-dominance models for predicting new families should allow practicing mating allocation to maximize the total genetic values for the propagation of elite material.

Validation with single-step SNPBLUP shows that evaluations can continue using a single mean of genotyped individuals, even with multiple breeds

BackgroundIn genomic prediction, it is common to centre the genotypes of single nucleotide polymorphisms based on the allele frequencies in the current population, rather than those in the base generation. The mean breeding value of non-genotyped animals is conditional on the mean performance of genotyped relatives, but can be corrected by fitting the mean performance of genotyped individuals as a fixed regression. The associated covariate vector has been referred to as a ‘J-factor’, which if fitted as a fixed effect can improve the accuracy and dispersion bias of sire genomic estimated breeding values (GEBV). To date, this has only been performed on populations with a single breed. Here, we investigated whether there was any benefit in fitting a separate J-factor for each breed in a three-way crossbred population, and in using pedigree-based expected or genome-based estimated breed fractions to define the J-factors.ResultsFor body weight at 7 days, dispersion bias decreased when fitting multiple J-factors, but only with a low proportion of genotyped individuals with selective genotyping. On average, the mean regression coefficients of validation records on those of GEBV increased with one J-factor compared to none, and further increased with multiple J-factors. However, for body weight at 35 days this was not observed. The accuracy of GEBV remained unchanged regardless of the J-factor method used. Differences between the J-factor methods were limited with correlations approaching 1 for the estimated covariate vector, the estimated coefficients of the regression on the J-factors, and the GEBV.ConclusionsBased on our results and in the particular design analysed here, i.e. all the animals with phenotype are of the same type of crossbreds, fitting a single J-factor should be sufficient, to reduce dispersion bias. Fitting multiple J-factors may reduce dispersion bias further but this depends on the trait and genotyping rate. For the crossbred population analysed, fitting multiple J-factors has no adverse consequences and if this is done, it does not matter if the breed fractions used are based on the pedigree-expectation or the genomic estimates. Finally, when GEBV are estimated from crossbred data, any observed bias can potentially be reduced by including a straightforward regression on actual breed proportions.

Observation of Genetic Gain with Instrumental Insemination of Honeybee Queens

Controlling mating in the honeybee (Apis mellifera) is part of one of the greatest challenges for the beekeeping industry's genetic selection programs due to specific characteristics of their reproduction. Several techniques for supervising honeybee mating with relative effective control have been developed over the years to allow honeybee selection. As part of this project, we compared the genetic gains for several colony performance traits, obtained using the BLUP-animal method, according to the selection pressure applied in controlled reproduction (directed fertilization versus instrumental insemination). Our results show similar genetic gains for hygienic behavior and honey production between colonies whether queens were fertilized naturally or via instrumental insemination, as well as similar or lower genetic gains for colonies with queens inseminated for spring development. In addition, we noticed greater fragility in queens following insemination. These findings show that instrumental insemination is an effective tool for reproductive control in genetic selection and for estimating breeding values more precisely. However, this technique does not result in queens of superior genetic quality for commercial purposes.

Pedigree and genome-based patterns of homozygosity in the South African Ayrshire, Holstein, and Jersey breeds.

The erosion of genetic diversity limits long-term genetic gain and impedes the sustainability of livestock production. In the South African (SA) dairy industry, the major commercial dairy breeds have been applying estimated breeding values (EBVs) and/or have been participating in Multiple Across Country Evaluations (MACE). The transition to genomic estimated breeding values (GEBVs) in selection strategies requires monitoring of the genetic diversity and inbreeding of current genotyped animals, especially considering the comparatively small population sizes of global dairy breeds in SA. This study aimed to perform a homozygosity-based evaluation of the SA Ayrshire (AYR), Holstein (HST), and Jersey (JER) dairy cattle breeds. Three sources of information, namely 1) single nucleotide polymorphism (SNP) genotypes (3,199 animals genotyped for 35,572 SNPs) 2) pedigree records (7,885 AYR; 28,391 HST; 18,755 JER), and 3) identified runs of homozygosity (ROH) segments were used to quantify inbreeding related parameters. The lowest pedigree completeness was for the HST population reducing from a value of 0.990 to 0.186 for generation depths of one to six. Across all breeds, 46.7% of the detected ROH were between 4 megabase pairs (Mb) and 8Mb in length. Two conserved homozygous haplotypes were identified in more than 70% of the JER population on Bos taurus autosome (BTA) 7. The JER breed displayed the highest level of inbreeding across all inbreeding coefficients. The mean (± standard deviation) pedigree-based inbreeding coefficient (FPED) ranged from 0.051 (±0.020) for AYR to 0.062 (±0.027) for JER, whereas SNP-based inbreeding coefficients (FSNP) ranged from 0.020 (HST) to 0.190 (JER) and ROH-based inbreeding coefficients, considering all ROH segment coverage (FROH), ranged from 0.053 (AYR) to 0.085 (JER). Within-breed Spearman correlations between pedigree-based and genome-based estimates ranged from weak (AYR: 0.132 between FPED and FROH calculated for ROH <4Mb in size) to moderate (HST: 0.584 between FPED and FSNP). Correlations strengthened between FPED and FROH as the ROH length category was considered lengthened, suggesting a dependency on breed-specific pedigree depth. The genomic homozygosity-based parameters studied proved useful in investigating the current inbreeding status of reference populations genotyped to implement genomic selection in the three most prominent South African dairy cattle breeds.

Exploring the statistical nature of independent chromosome segments

Previous research in genomics developed the concept of independent chromosome segments and suggested that genomic prediction works on capturing the segment effects rather than on LD to individual QTL. The number of independent chromosome segments (Me) was posited as 4NeL, a function of effective population size (Ne) and genome length (L) in Morgans (M). The objective of this study was to determine whether the physical segments are approximately consecutive haplotype blocks of length ¼ M, with the number of haplotype blocks for each physical location equal to Ne. In a simulated population, the number of animals randomly selected as reference animals is represented as Na. For all animals, the genome was split into equal-sized segments represented as Ns. For each specific location, segments of non-reference animals were assigned the most similar segment of one reference animal. Genomic analyses estimated the value of the segment effects, and breeding values were a sum of all segment effects for a specific animal. Accuracies of segment effects were calculated by correlating the true breeding values (TBV) and the breeding values based on segments. Segment effect accuracies were compared with the true accuracy calculated by the correlation of TBV and genomic estimated breeding values (GEBV) computed using GBLUP. Accuracies were maximized at Na=Ne, Ns=4 L, but they were not as high as in GBLUP. Accuracies may be smaller using the statistical concept of segments due to approximations based on computing limitations, as the origin of each segment and the recombination sites were unknown in the simulation. Therefore, random animals served as reference animals, and each non-reference animal received the most similar segment of a reference animal instead of a linear combination of such segments. Genomic selection acts partially on ¼ M long chromosome segments, and using the statistical definition of segments moderately explains the accuracy. More complex simulations are needed to investigate the issue thoroughly.

Single-step genome-wide association analyses of claw horn lesions in Holstein cattle using linear and threshold models

BackgroundLameness in dairy cattle is primarily caused by foot lesions including the claw horn lesions (CHL) of sole haemorrhage (SH), sole ulcers (SU), and white line disease (WL). This study investigated the genetic architecture of the three CHL based on detailed animal phenotypes of CHL susceptibility and severity. Estimation of genetic parameters and breeding values, single-step genome-wide association analyses, and functional enrichment analyses were performed.ResultsThe studied traits were under genetic control with a low to moderate heritability. Heritability estimates of SH and SU susceptibility on the liability scale were 0.29 and 0.35, respectively. Heritability of SH and SU severity were 0.12 and 0.07, respectively. Heritability of WL was relatively lower, indicating stronger environmental influence on the presence and development of WL than the other two CHL. Genetic correlations between SH and SU were high (0.98 for lesion susceptibility and 0.59 for lesion severity), whereas genetic correlations of SH and SU with WL also tended to be positive. Candidate quantitative trait loci (QTL) were identified for all CHL, including some on Bos taurus chromosome (BTA) 3 and 18 with potential pleiotropic effects associated with multiple foot lesion traits. A genomic window of 0.65 Mb on BTA3 explained 0.41, 0.50, 0.38, and 0.49% of the genetic variance for SH susceptibility, SH severity, WL susceptibility, and WL severity, respectively. Another window on BTA18 explained 0.66, 0.41, and 0.70% of the genetic variance for SH susceptibility, SU susceptibility, and SU severity, respectively. The candidate genomic regions associated with CHL harbour annotated genes that are linked to immune system function and inflammation responses, lipid metabolism, calcium ion activities, and neuronal excitability.ConclusionsThe studied CHL are complex traits with a polygenic mode of inheritance. Most traits exhibited genetic variation suggesting that animal resistance to CHL can be improved with breeding. The CHL traits were positively correlated, which will facilitate genetic improvement for resistance to CHL as a whole. Candidate genomic regions associated with lesion susceptibility and severity of SH, SU, and WL provide insights into a global profile of the genetic background underlying CHL and inform genetic improvement programmes aiming at enhancing foot health in dairy cattle.

Validating the female fertility estimated breeding value in Australian commercial dairy herds

We conducted a retrospective cohort study to validate the efficacy of the Australian multitrait fertility estimated breeding value (EBV). We did this by determining its associations with phenotypic measures of reproductive performance (i.e., submission rate, first service conception rate, and early calving). Our secondary aim was to report the associations between these reproductive outcomes and management and climate-related factors hypothesized to affect fertility. Our study population included 38 pasture-based dairy herds from the northern Victorian irrigation region in Australia. We collected records for 86,974 cows with 219,156 lactations and 438,578 mating events from the date on which managers started herd recording until December 2016, comprising both fertility-related data such as insemination records, calving dates, and pregnancy test results, and systems-related data such as production, herd size, and calving pattern. We also collected hourly data from 2004 to 2017 from the closest available weather station to account for climate-related factors (i.e., temperature humidity index; THI). Multilevel Cox proportional hazard models were used to analyze time-to-event outcomes (days to first service, days to cow calving following the planned herd calving start date), and multilevel logistic regression models for binomial outcomes (conception to first service) in the Holstein-Friesian and Jersey breeds. A 1-unit increase in daughter fertility EBV was associated with a 5.4 and 8.2% increase in the daily hazard of calving in the Holstein-Friesian and Jersey breeds respectively. These are relative increases (i.e., a Holstein-Friesian herd with a 60% 6-wk in-calf rate would see an improvement to 63.2% with a 1-unit increase in herd fertility EBV). Similar results were obtained for submission and conception rate. Associations between 120-d milk yield and reproductive outcome were complicated by interactions with 120-d protein percentage and calving age, depending on the breed and outcome. In general, we found that the reproductive performance of high milk-yielding animals deteriorated faster with age than low milk-yielding animals, and high protein percentage exacerbated the differences between low and high milk-yielding animals. Climate-related factors were also associated with fertility, with a 1-unit increase in maximum THI decreasing first service conception rate by 1.2% for Holstein-Friesians but having no statistically significant association in the Jersey breed. However, THI had a negative association in both breeds on the daily hazard of calving. Our study validates the efficacy of the daughter fertility EBV for improving herd reproductive performance and identifies significant associations between 120-d milk and protein yields and THI on the fertility of Australian dairy cows.

Selective genotyping strategies for a sib test scheme of a broiler breeder program

BackgroundIn broiler breeding, genotype-by-environment interaction is known to result in a genetic correlation between body weight measured in bio-secure and commercial environments that is substantially less than 1. Thus, measuring body weights on sibs of selection candidates in a commercial environment and genotyping them could increase genetic progress. Using real data, the aim of this study was to evaluate which genotyping strategy and which proportion of sibs placed in the commercial environment should be genotyped to optimize a sib-testing breeding program in broilers. Phenotypic body weight and genomic information were collected on all sibs raised in a commercial environment, which allowed to retrospectively analyze different sampling strategies and genotyping proportions.ResultsAccuracies of genomic estimated breeding values (GEBV) obtained with the different genotyping strategies were assessed by computing their correlation with GEBV obtained when all sibs in the commercial environment were genotyped. Results showed that, compared to random sampling (RND), genotyping sibs with extreme phenotypes (EXT) resulted in higher GEBV accuracy across all genotyping proportions, especially for genotyping proportions of 12.5% or 25%, which resulted in correlations of 0.91 vs 0.88 for 12.5% and 0.94 vs 0.91 for 25% genotyped. Including pedigree on birds with phenotype in the commercial environment that were not genotyped increased accuracy at lower genotyping proportions, especially for the RND strategy (correlations of 0.88 vs 0.65 at 12.5% and 0.91 vs 0.80 at 25%), and a smaller but still substantial increase in accuracy for the EXT strategy (0.91 vs 0.79 for 12.5% and 0.94 vs 0.88 for 25% genotyped). Dispersion bias was virtually absent for RND if 25% or more birds were genotyped. However, GEBV were considerably inflated for EXT, especially when the proportion genotyped was low, which was further exacerbated if the pedigree of non-genotyped sibs was excluded.ConclusionsWhen less than 75% of all animals placed in a commercial environment are genotyped, it is recommended to use the EXT strategy, because it yields the highest accuracy. However, caution should be taken when interpreting the resulting GEBV because they will be over-dispersed. When 75% or more of the animals are genotyped, random sampling is recommended because it yields virtually no bias of GEBV and results in similar accuracies as the EXT strategy.

First large-scale genomic prediction in the honey bee

Genomic selection has increased genetic gain in several livestock species, but due to the complicated genetics and reproduction biology not yet in honey bees. Recently, 2970 queens were genotyped to gather a reference population. For the application of genomic selection in honey bees, this study analyzes the accuracy and bias of pedigree-based and genomic breeding values for honey yield, three workability traits, and two traits for resistance against the parasite Varroa destructor. For breeding value estimation, we use a honey bee-specific model with maternal and direct effects, to account for the contributions of the workers and the queen of a colony to the phenotypes. We conducted a validation for the last generation and a five-fold cross-validation. In the validation for the last generation, the accuracy of pedigree-based estimated breeding values was 0.12 for honey yield, and ranged from 0.42 to 0.61 for the workability traits. The inclusion of genomic marker data improved these accuracies to 0.23 for honey yield, and a range from 0.44 to 0.65 for the workability traits. The inclusion of genomic data did not improve the accuracy of the disease-related traits. Traits with high heritability for maternal effects compared to the heritability for direct effects showed the most promising results. For all traits except the Varroa resistance traits, the bias with genomic methods was on a similar level compared to the bias with pedigree-based BLUP. The results show that genomic selection can successfully be applied to honey bees.

A transgenerational study on the effect of great-granddam birth month on granddaughter EBV for production traits in Italian Simmental cattle

Heat tolerance is a key feature of resilient animals. Offspring of animals that suffer environmental stress during pregnancy could show physiological, morphological, and metabolic modifications. This is due to a dynamic reprogramming of the epigenetics of the mammalian genome that occurs in the early life cycle. Thus, the aim of this study was to investigate the extent of the transgenerational effect of heat stress during the pregnancy of Italian Simmental cows. The effects of dam and granddam birth months (as indicator of pregnancy period) on their daughter and granddaughter estimated breeding values (EBV) for some dairy traits as well as of the temperature-humidity index (THI) during the pregnancy were tested. A total of 128,437 EBV (milk, fat, and protein yields, and somatic cell score) were provided by the Italian Association of Simmental Breeders. The best birth months (of both dam and granddam) for milk yield and protein yield were May and June, whereas the worst were January and March. Great-granddam pregnancies developed during the winter and spring seasons positively affected the EBV for milk and protein yields of their great-granddaughters; in contrast, pregnancies during summer and autumn had negative effects. These findings were confirmed by the effects of maximum and minimum THI in different parts of the great-granddam pregnancy on the performances of their great-granddaughters. Thus, a negative effect of high temperatures during the pregnancy of female ancestors was observed. Results of the present study suggest a transgenerational epigenetic inheritance in Italian Simmental cattle due to environmental stressors.

Evaluation of the Crossbreeding Scheme and Farmers’ Perception of Awassi and Dorper Crossbred Sheep

The objectives of this study were to characterize the performance of crossbred sheep and crossbreeding scheme and to solicit farmers’ opinions about crossbred sheep. A total of 238 respondents were selected from five districts and data were collected through a personal interview, focus group discussion, measurement of animals, and field observation. Continuous type of data were analyzed using different procedures of SAS, while the estimated breeding value (EBV) for distributed rams was estimated by the WOMBAT software fitting animal model. The findings showed that sheep were the most significant species and that they were primarily raised for meat consumption and income generation. The proportions of crossbred sheep were 79.6% in Angot, 61.3% in Dessie zuria, 64.2% in Gazo, 80.6% in Legambo, and 27.5% in the Kobo area. The sheep crossbreeding program was constrained by feed shortage, lack of improved genotype, and diseases. The ongoing sheep crossbreeding program lacks a proper recording scheme, rams were not selected based on genetic merit, not exchanged in time, lack periodic importation of genetically unrelated rams, and lacks fixing the maximum exotic gene level suitable for the production system. The EBVs of Dorper crossbred rams for birth weight and weaning weight were 0.007 kg and 0.273 kg, respectively. Tikur sheep had an earlier age at first lambing (AFL) than Awassi x Tikur sheep. Likewise, the indigenous Wollo sheep had earlier AFL, short lambing interval, and produced a large number of lambs per lifetime than their Awassi crossbreds. Dorper crossbreds were preferred (9.30 - 72.2 times greater and P &lt; 0.001 ) for their growth rate, physical appearance, preference in the market or price, and milk production of ewes compared with indigenous Tumele sheep. The Awassi crossbred sheep were preferred (odds ratio = 12.7– 90.0 and P &lt; 0.001 ) due to their good physical appearance, fast growth rate, wool yield, milk yield, and better preference in the market compared with Tikur and Wollo sheep breeds. The implementation of the crossbreeding program needs some sort of revision, monitoring, and periodic evaluation. Besides, it must be accompanied by improved management to exploit the expected benefits from the crossbreeding program in the low-input production system.

Sample size determination for training set optimization in genomic prediction

Key messageA practical approach is developed to determine a cost-effective optimal training set for selective phenotyping in a genomic prediction study. An R function is provided to facilitate the application of the approach.Genomic prediction (GP) is a statistical method used to select quantitative traits in animal or plant breeding. For this purpose, a statistical prediction model is first built that uses phenotypic and genotypic data in a training set. The trained model is then used to predict genomic estimated breeding values (GEBVs) for individuals within a breeding population. Setting the sample size of the training set usually takes into account time and space constraints that are inevitable in an agricultural experiment. However, the determination of the sample size remains an unresolved issue for a GP study. By applying the logistic growth curve to identify prediction accuracy for the GEBVs and the training set size, a practical approach was developed to determine a cost-effective optimal training set for a given genome dataset with known genotypic data. Three real genome datasets were used to illustrate the proposed approach. An R function is provided to facilitate widespread application of this approach to sample size determination, which can help breeders to identify a set of genotypes with an economical sample size for selective phenotyping.

Single-step genomic evaluation for growth traits in a Mexican Braunvieh cattle population.

The objective was to compare (pedigree-based) BLUP, genomic BLUP (GBLUP), and single-step GBLUP (ssGBLUP) methods for genomic evaluation of growth traits in a Mexican Braunvieh cattle population. Birth (BW), weaning (WW), and yearling weight (YW) data of a Mexican Braunvieh cattle population were analyzed with BLUP, GBLUP, and ssGBLUP methods. These methods are differentiated by the additive genetic relationship matrix included in the model and the animals under evaluation. The predictive ability of the model was evaluated using random partitions of the data in training and testing sets, consistently predicting about 20% of genotyped animals on all occasions. For each partition, the Pearson correlation coefficient between adjusted phenotypes for fixed effects and non-genetic random effects and the estimated breeding values (EBV) were computed. The random contemporary group (CG) effect explained about 50, 45, and 35% of the phenotypic variance in BW, WW, and YW, respectively. For the three methods, the CG effect explained the highest proportion of the phenotypic variances (except for YW-GBLUP). The heritability estimate obtained with GBLUP was the lowest for BW, while the highest heritability was obtained with BLUP. For WW, the highest heritability estimate was obtained with BLUP, the estimates obtained with GBLUP and ssGBLUP were similar. For YW, the heritability estimates obtained with GBLUP and BLUP were similar, and the lowest heritability was obtained with ssGBLUP. Pearson correlation coefficients between adjusted phenotypes for non-genetic effects and EBVs were the highest for BLUP, followed by ssBLUP and GBLUP. The successful implementation of genetic evaluations that include genotyped and non-genotyped animals in our study indicate a promising method for use in genetic improvement programs of Braunvieh cattle. Our findings showed that simultaneous evaluation of genotyped and non-genotyped animals improved prediction accuracy for growth traits even with a limited number of genotyped animals.

Adjusted fence height: an improved phenotype for the genetic evaluation of show jumping performance in Warmblood horses

BackgroundShow jumping is one of the most popular disciplines in the horse sector, which makes success in show jumping competitions an important breeding goal for many studbooks. Therefore, the genetic evaluation of show jumping performance is of major interest and this is the case for two Belgian Warmblood studbooks: the Belgian Warmblood horse and Zangersheide. In this study, first an improved phenotype for show jumping performance was developed, i.e. adjusted fence height based on a new non-arbitrary method to scale ranking and competition level, which are two major components of success in competitions. Second, we assessed the importance of including a rider effect in genetic models for show jumping performance, this effect being under debate in sport horse breeding. Third, genetic models based on elementary performances and one model based on a summarized performance were compared in terms of model fit, heritabilities and the stability of estimated breeding values to define the most suitable one for the genetic evaluation of show jumping performance.ResultsIn this study, more than 600,000 Belgian competition records and almost 81,000 horses were used. Genetic evaluations were developed based on elementary performances (Blom-transformed ranking and adjusted fence height) and on a summarized performance (highest level achieved). Estimated heritabilities of Blom-transformed ranking, adjusted fence height and highest level achieved were 0.09, 0.12 and 0.39, respectively. Including a rider effect improved the models for genetic evaluations. Estimated genetic correlations between the studied models were moderate to high (rg = 0.60–0.99). With the best fit model, the accuracy of the estimated breeding value (EBV) for adjusted fence height reached 0.70 for a larger number of stallions and for stallions that tended to be younger.ConclusionsWe recommend breeders to implement this new phenotype ‘adjusted fence height’ in breeding programs. It is moderately to highly correlated with Blom-transformed ranking and highest level achieved, a proxy for lifetime success, and is available for selection candidates from an early age onwards.