Accuracy of the unified approach in maternally influenced traits  illustrated by a simulation study in the honey bee (Apis mellifera)
 Pooja Gupta^{1, 2}Email author,
 Norbert Reinsch^{3},
 Andreas Spötter^{1},
 Tim Conrad^{2} and
 Kaspar Bienefeld^{1}
DOI: 10.1186/147121561436
© Gupta et al.; licensee BioMed Central Ltd. 2013
Received: 8 December 2012
Accepted: 30 April 2013
Published: 6 May 2013
Abstract
Background
The honey bee is an economically important species. With a rapid decline of the honey bee population, it is necessary to implement an improved genetic evaluation methodology. In this study, we investigated the applicability of the unified approach and its impact on the accuracy of estimation of breeding values for maternally influenced traits on a simulated dataset for the honey bee. Due to the limitation to the number of individuals that can be genotyped in a honey bee population, the unified approach can be an efficient strategy to increase the genetic gain and to provide a more accurate estimation of breeding values. We calculated the accuracy of estimated breeding values for two evaluation approaches, the unified approach and the traditional pedigree based approach. We analyzed the effects of different heritabilities as well as genetic correlation between direct and maternal effects on the accuracy of estimation of direct, maternal and overall breeding values (sum of maternal and direct breeding values). The genetic and reproductive biology of the honey bee was accounted for by taking into consideration characteristics such as colony structure, uncertain paternity, overlapping generations and polyandry. In addition, we used a modified numerator relationship matrix and a realistic genome for the honey bee.
Results
For all values of heritability and correlation, the accuracy of overall estimated breeding values increased significantly with the unified approach. The increase in accuracy was always higher for the case when there was no correlation as compared to the case where a negative correlation existed between maternal and direct effects.
Conclusions
Our study shows that the unified approach is a useful methodology for genetic evaluation in honey bees, and can contribute immensely to the improvement of traits of apicultural interest such as resistance to Varroa or production and behavioural traits. In particular, the study is of great interest for cases where negative correlation between maternal and direct effects and uncertain paternity exist, thus, is of relevance for other species as well. The study also provides an important framework for simulating genomic and pedigree datasets that will prove to be helpful for future studies.
Keywords
Maternal effects Singlestep approach Genetic evaluation Single nucleotide polymorphismBackground
A colony trait (e.g. honey and wax production) in the honey bee is comparable to maternally influenced traits in mammals such as birth and weaning weight; thus, it can be partitioned into the additive genetic effect of the queen (maternal genetic effects) and the additive genetic effects of the progeny workers (direct genetic effects). The queen mediates its effect through heritable characters like egg laying rate or pheromone production in the hive whereas workers affect a trait through their hoarding behaviour or production of and responsiveness to pheromones. Until now, genetic evaluation in the honey bee has been implemented using a pedigree based BLUPanimal model with maternal and direct genetic effects [1]. In the last decades, genetic evaluation strategies in agricultural animals have undergone remarkable advancement as a result of the introduction of genomic selection strategy. Genomic selection [2], which is based on highdensity molecular marker information, has now become the ‘state of the art’ method for genetic evaluation. For many livestock species, for example, in US dairy cattle [3] and pigs [4], a multistep procedure for genomic selection was proposed. This multistep procedure has certain disadvantages with respect to the honey bee. Due to economical and technical constraints, it may not be possible to genotype a large number of animals in the honey bee population. Thus, instead of a multistep procedure, we employed a singlestep unified approach in our study. The unified approach was proposed by Legarra et al. [5] and Christensen and Lund [6], and it combines full pedigree and genomic information from both genotyped and ungenotyped individuals. The advantage of this procedure over the multistep approach is that it gives a more accurate estimate of breeding values for ungenotyped animals [6, 7] and is resistant to selection bias [8]. Moreover, it is simpler to implement as compared to the multistep approach and provides an easy extension to a multitrait model [9] with maternal effects in honey bees.
We performed a simulation study to investigate the impact of the unified approach on the accuracy of estimated breeding values in honey bees. Similar to the case of the honey bee, other species also have a situation where genetic evaluation needs to account for maternal effects and uncertain paternity, e.g. weaning weight in beef cattle is a maternally influenced trait. Besides, cows can be exposed to more than one male in a herd and pasture paddock within the same breeding season, thus generating uncertainty on paternity assignments and adversely affecting the accuracy of breeding value predictions [10]. In this study, we simultaneously dealt with the effect of uncertain paternity and maternal effects on genomic predictions. Therefore, the study is of broad interest, and can be of use for other species where maternal effects or/and uncertain paternity exists. With the exception of an abstract contribution by Lourenco et al. [11] this is, to the authors’ knowledge, the first simulation study where the application of a unified approach is evaluated for a trait with direct and maternal genetic effects.
Methods
This study consisted of two main steps. In the first step, a dataset was simulated for a honey bee population, which involved the modelling and simulation of the population structure, genome, correlation between maternal and direct effects, heritability, true breeding values, genetic, phenotypic and residual variances. The simulated dataset was close to realistic scenarios and in agreement with the genetic and reproductive peculiarities of the honey bee. In the second step, genetic evaluation was performed using the unified approach and the traditional pedigree based BLUP approach.
Population structure
Base population in linkage disequilibrium
A random mating population was simulated for 1000 generations to obtain a base population in mutationdrift equilibrium with linkage disequilibrium (LD) [12]. The simulated generations were discrete and nonoverlapping. The population size was kept constant in every generation, and consisted of 500 sire queens and 50 dam queens. A sire queen represented a droneproducing queen that produced only drones whereas a dam queen represented a queen that mated with the drones to produce offspring for the next generation. The resulting base population was assumed to be noninbred and unrelated.
Mating and selection scheme
Population characteristics specific to the honey bee
To construct a population similar to that used in the genetic evaluation program of honey bees, we constructed a dummy sire and an average worker (representing direct effects) in the pedigree. Generations following the base population were overlapping and mating was polyandrous as in the normal breeding population. These characteristics are described in more detail in the following section.
Construction of a dummy sire and an average worker
A colony is formed by a queen and its progeny comprising several thousand workers. Since it is impossible to include all workers of a colony for genetic evaluation, an average worker was constructed that represented all workers of a colony. It was assumed that one average worker existed for each potentialdam queen/dam queen in the pedigree.
Modelling polyandry and overlapping generations
Pedigree, phenotypic and genomic information
A phenotypic value in the honey bee represents an observation for the whole colony and thus, cannot be decomposed into individual phenotypic values of a queen and an average worker. Therefore, both the queen and the average worker of a colony were assigned the same colony phenotypic value. It was assumed that pedigree records were available for all generations; phenotypes were available for all dam queens (and the corresponding average worker) in the base generation and all potentialdam queens (and the corresponding average worker) in all but the last generation. Genotyping information was available for all dam queens in the base generation and all potentialdam queens.
Genome
Summary of the simulated chromosome lengths and number of SNP on each chromosome of the honey bee
Chromosome  Chromosome length (in basepairs)  Number of SNP 

Chromosome 1  29 893 408  14 137 
Chromosome 2  15 549 267  6 335 
Chromosome 3  13 234 341  7 119 
Chromosome 4  12 718 334  5 589 
Chromosome 5  14 363 272  6 330 
Chromosome 6  18 472 937  7 877 
Chromosome 7  13 219 345  5 973 
Chromosome 8  13 546 544  6 235 
Chromosome 9  11 120 453  5 578 
Chromosome 10  12 965 953  5 068 
Chromosome 11  14 726 556  6 957 
Chromosome 12  11 902 654  5 812 
Chromosome 13  10 288 499  5 082 
Chromosome 14  10 253 655  4 874 
Chromosome 15  10 167 229  3 879 
Chromosome 16  7 207 165  3 155 
Total  219 629 612  100 000 
Correlation between maternal and direct effects
Studies in honey bees ([21]; Ehrhardt and Bienefeld, unpublished results) have shown that there is a strong negative correlation between maternal (queen) and direct (worker) effects. To model this, a total of 250 QTL were simulated, out of which 86 loci controlled the direct effects, 78 pleiotropic loci controlled both the direct and maternal effects and the remaining 86 loci controlled the maternal effects. To establish a negative correlation between maternal and direct effects, signs for QTL effects for maternal and direct genetic effects were chosen opposite to each other at the pleiotropic loci. The level of negative correlation was determined by the number of pleiotropic loci. No correlation between maternal and direct effects was obtained by randomly choosing signs for QTL effects for maternal and direct effects at the pleiotropic loci. The simulated value of correlation (r_{ qw }) was obtained by estimating the correlation between the maternal and direct true breeding values.
True breeding values and phenotypic values
Maternal and direct true breeding values were simulated for all dam queens of the base population and all potentialdam queens from generations 1–5. True breeding values for maternal (TBV_{ q }) and direct effects (TBV_{ w }) for a queen were calculated using the formula $\mathit{TB}{V}_{q}^{i}={\displaystyle \sum _{j}{q}_{q}^{\mathit{ij}}{a}^{j}}$ and $\mathit{TB}{V}_{w}^{i}={\displaystyle \sum _{k}{q}_{w}^{\mathit{ik}}{a}^{k}}$ where TBV_{ q }^{ i } and TBV_{ w }^{ i } are the maternal and direct true breeding values for the i^{ th } queen, respectively. q_{ q }^{ ij } and q_{ w }^{ ik } are QTL genotypes of the i^{ th } queen at the j^{ th } and k^{ th } QTL controlling the maternal and direct effects, respectively and has a value of 1 or −1 for the homozygous genotypes or 0 for the heterozygous genotype. a^{ j } and a^{ k } are allele substitution effects at the j^{ th } and k^{ th } QTL.
The overall true breeding value of a queen was the sum of its maternal and direct true breeding values. The phenotype of each queen was obtained by adding the overall true breeding value of a queen to a residual value drawn from a normal distribution N(0, σ_{ e }^{2}). The way the value for residual variance (σ_{ e }^{2}) was chosen is explained in the later section.
Genetic variance
Variance and covariance of maternal and direct effects
Variances of maternal (σ_{ q }^{2}) and direct (σ_{ w }^{2}) effects were obtained by calculating the variance of the simulated maternal and direct true breeding values, respectively. The covariance between maternal and direct effects (σ_{ qw }) was obtained by calculating the covariance between the maternal and direct true breeding values.
Total genetic variance
Usually a breeding value is defined as twice the expected deviation of an individual's progeny from the mean, or twice the ‘transmitting ability’ of an individual [22]. If we consider a complete colony as ‘offspring’ of a queen, then this colony comprises a daughter (the queen) and a family of granddaughters (the workers). These offspring express 1/2 of the mother's maternal breeding value and 1/4 of the granddam's direct breeding value. In this case, the overall true or estimated breeding value would be defined as twice the 1/2 of the maternal breeding value of a queen plus twice the 1/4 of its direct breeding value (i.e. two times the expected deviation of ‘progeny’ from the mean, provided all other relatives have average breeding values of zero). Thus, maternal and direct breeding values get a weight of 1 and 0.5, respectively. The total genetic variance (σ_{ g }^{2}) would become σ_{ q }^{2} + 0.25σ_{ w }^{2} + σ_{ qw } (the latter from 2 × 1 × 0.5 × σ_{ qw }). However, for the sake of easy comparison and interpretation, the overall breeding value was taken as a sum of the direct and maternal breeding values of a queen and the total genetic variance was taken as a sum of variance of maternal effects, direct effects and twice the covariance between them, and can be expressed as σ_{ g }^{2} = σ_{ q }^{2} + σ_{ w }^{2} + 2σ_{ qw }.
Phenotypic variance, residual variance and maternal and direct heritability
Heritability of direct effects for different values of simulated heritability of maternal effects and correlation between maternal and direct effects
Simulatedh_{ m }^{2}  Corr_{(m,d)}  Achievedh_{ d }^{2}(SE) 

0.150  0  0.162 (0.005) 
0.150  −0.46  0.155 (0.005) 
0.250  0  0.270 (0.008) 
0.250  −0.46  0.259 (0.008) 
0.350  0  0.377 (0.011) 
0.350  −0.46  0.362 (0.011) 
Estimation of breeding values
where y is a vector of records of the colonies, b is a vector of fixed effects, u_{ 1 } is a vector of random direct effects, u_{ 2 } is a vector of random maternal effects, e is a vector of random residual effects, X is an incidence matrix relating observations to the corresponding environment, Z_{ 1 } and Z_{ 2 } are the incidence matrices relating observations to the corresponding direct effects and maternal effects, respectively.
Estimation of breeding values was done using the following two approaches: (1) the traditional BLUP approach (PED_BLUP) based on a numerator relationship matrix (A) constructed from pedigree information and (2) the unified approach (UNI_BLUP) based on a combined relationship matrix (H) constructed from pedigree and genomic information.
Relationship matrix constructed from pedigree data
Relationship matrix constructed from pedigree and genomic data
In the honey bee pedigree, a dummy sire and an average worker represent a group of individuals and thus, it is not possible to get individual genotyping data. Moreover, it is not possible to obtain genotyping information from all queens in the population. Using the unified approach is advantageous for honey bees as genomic information for genotyped queens can be integrated with pedigree information from genotyped as well as nongenotyped individuals resulting in a combined relationship matrix H. A genomic matrix (G) was constructed for the 2550 queens with genotyping data. Different methods have been developed to derive the G matrix [24, 25]. We chose a methodology proposed by VanRaden [24]. The G matrix was obtained from ZZ’/2 ∑ p_{ i }(1 − p_{ i }), where Z is equal to M − P, M is the matrix specifying marker alleles inherited by each individual and P is equal to 2(p_{ i } − 0.5) with p_{ i } being the frequency of second allele at locus i in the base population. In order to avoid a singular G matrix [5, 7, 24], a modified matrix (G_{ w }) was constructed using a weighing factor (w), given as G_{ w } = w G + (1 − w)A_{ 11 }. Christensen and Lund [6] suggested that (1 − w) could be interpreted as the relative weight on the polygenic effect. Aguilar et al. [7] reported that the weights were not critical, and using a value of 0.95 or 0.98 caused negligible difference in the results. For this study, the value of w was taken as 0.99 [6].
Simulated values for the genetic and residual variance were used for estimating the breeding values. For both approaches, statistics for the achieved heritability of direct effects and accuracies for the overall, maternal and direct estimated breeding values were based on 20 replicated simulations. The accuracy was reported as a correlation between the estimated and true breeding values [26] for 500 ‘juvenile queens’ constituted by potentialdam queens in the last generation and 2550 ‘all queens’ constituted by dam queens in the base population and potentialdam queens in all generations. All calculations were performed in MATLAB.
Results
Accuracy of the overall estimated breeding values
The accuracy of overall EBV (the sum of maternal and direct breeding values) in case of 500 juvenile queens and all 2550 queens in the pedigree
h _{ m } ^{2}  Method  Corr_{(m,d)}  Accuracy for JQ (SE)  Accuracy for AQ (SE) 

0.15  UNI  0  0.468 ^{a,b,c,d} (0.010)  0.661 ^{a,b,c,d} (0.005) 
PED  0  0.363 (0.017)  0.603 (0.007)  
UNI  −0.46  0.381 ^{a,b,c,d} (0.021)  0.555 ^{a,b,c,d} (0.010)  
PED  −0.46  0.295 (0.023)  0.489 (0.009)  
0.25  UNI  0  0.542 ^{a,b,c,e} (0.009)  0.756 ^{a,b,c,e} (0.006) 
PED  0  0.420 (0.015)  0.710 (0.008)  
UNI  −0.46  0.449 ^{a,b,c} (0.018)  0.640 ^{a,b,c,e} (0.009)  
PED  −0.46  0.348 (0.021)  0.577 (0.008)  
0.35  UNI  0  0.604 ^{a,b,d,e} (0.009)  0.832 ^{a,b,d,e} (0.008) 
PED  0  0.467 (0.012)  0.800 (0.010)  
UNI  −0.46  0.498 ^{a,b,d} (0.017)  0.700 ^{a,b,d,e} (0.008)  
PED  −0.46  0.388 (0.019)  0.642 (0.008) 
For juvenile queens, the accuracy of overall estimated breeding values was significantly higher with the UNI_BLUP approach (P < 0.05) as compared to the PED_BLUP approach for all values of heritability and correlation between maternal and direct effects. The increase in accuracy by UNI_BLUP was approximately 0.1 (or 29%) for most of the cases.
Similar to juvenile queens, the accuracy of overall estimated breeding values for all queens was higher with the UNI_BLUP approach (P < 0.05) than the PED_BLUP for all values of heritability and correlation between maternal and direct effects. The percentage increase in accuracy for the case of no correlation between maternal and direct effects at maternal heritabilities of 0.15, 0.25 and 0.35 was approximately 9.6%, 6.5% and 4.0%, respectively. In case of a negative correlation of −0.46, the percentage increase in accuracy was approximately 13.5%, 10.9% and 9.0% at maternal heritabilities of 0.15, 0.25 and 0.35, respectively.
From these results we can conclude that the UNI_BLUP approach performed better than the PED_BLUP and the accuracy of overall estimated breeding values increased considerably with the UNI_BLUP approach.
Accuracy of the maternal and direct estimated breeding values
The accuracy of direct and maternal EBV in case of 500 juvenile queens and 2550 all queens in the pedigree
h _{ m } ^{2}  Method  Corr_{(m,d)}  Accuracy of direct EBV for JQ (SE)  Accuracy of maternal EBV for JQ (SE)  Accuracy of direct EBV for AQ (SE)  Accuracy of maternal EBV for AQ (SE) 

0.15  UNI  0  0.323 ^{a,b,c,d} (0.015)  0.279 ^{a,b,c,d} (0.015)  0.446 ^{a,b,c,d} (0.011)  0.420 ^{a,b,c,d} (0.009) 
PED  0  0.227 (0.023)  0.225 (0.019)  0.406 (0.012)  0.381 (0.010)  
UNI  −0.46  0.115 ^{b,d} (0.023)  0.127 ^{b} (0.031)  0.225 ^{b,d} (0.018)  0.223 ^{b,d} (0.013)  
PED  −0.46  0.059 (0.025)  0.103 (0.030)  0.186 (0.016)  0.208 (0.012)  
0.25  UNI  0  0.373 ^{a,b,c} (0.016)  0.330 ^{a,b,c,e} (0.014)  0.510 ^{b,c,e} (0.013)  0.482 ^{a,b,c,e} (0.010) 
PED  0  0.268 (0.021)  0.260 (0.018)  0.474 (0.013)  0.447 (0.011)  
UNI  −0.46  0.154 ^{b} (0.023)  0.154 ^{b} (0.031)  0.272 ^{b} (0.017)  0.257 ^{b} (0.014)  
PED  −0.46  0.085 (0.025)  0.125 (0.030)  0.231 (0.016)  0.240 (0.013)  
0.35  UNI  0  0.418 ^{a,b,d} (0.017)  0.371 ^{a,b,d,e} (0.014)  0.566 ^{b,d,e} (0.015)  0.527 ^{b,d,e} (0.011) 
PED  0  0.307 (0.018)  0.287 (0.017)  0.538 (0.015)  0.496 (0.013)  
UNI  −0.46  0.186 ^{a,b,d} (0.024)  0.173 ^{b} (0.031)  0.308 ^{b,d} (0.017)  0.280 ^{b,d} (0.015)  
PED  −0.46  0.110 (0.025)  0.138 (0.030)  0.268 (0.016)  0.258 (0.014) 
Effect of correlation and heritability
Both low heritability and negative correlation contribute to a lower genetic variance which leads to a decrease in the accuracy. The accuracy of overall estimated breeding values was reduced as a result of negative correlation in comparison to the case where maternal and genetic effects had no correlation (Table 3; P < 0.05). Similarly, the accuracy of overall estimated breeding values increased as the heritability increased (Table 3; P < 0.05). The only exception, where no significant difference was observed, was between maternal heritabilities of 0.25 and 0.35 at a negative correlation of −0.46 for juvenile queens, although the accuracy was higher for high heritability. This can be explained by the fact that the impact of using genomic information is smaller for traits with high heritablities.
The accuracy of maternal and direct estimated breeding values (Table 4) was higher for high values of heritability indicating a similar trend as the overall estimated breeding values. The difference was significant in most cases for any two compared values of heritability with no correlation and between heritability of 0.15 and 0.35 with negative correlation.
Discussion
The study provided comparative insight into genetic evaluation performed using: (1) the traditional BLUP approach based on pedigree data and (2) the unified approach based on both pedigree and marker data. In this study, we investigated the accuracy of overall, direct and maternal estimated breeding values as well as the influence of heritability of the trait and the genetic correlation between maternal and direct effects on the accuracy.
It has been reported in honey bees that most economically important traits have low to medium heritability [27–29]. Therefore, we also simulated heritabilities in the same range. The extremely negative estimates of genetic correlation between maternal and direct effects have often been questioned ([30]; Ehrhardt and Bienefeld, unpublished results), therefore, we simulated a general value of correlation of −0.46 which exists in other species as well [31–33] and compared it to a case with no correlation between maternal and direct effects.
Unlike previous studies [6, 7, 34], our study takes into account the influence of maternal and direct effects. We observed that the accuracy of overall estimated breeding values (Table 3) increased considerably with the unified approach for all scenarios of heritability and correlation with significant Pvalues (< 0.05). A higher gain in the accuracy of overall estimated breeding values was observed for juvenile animals. It is desired that the gain in accuracy is higher for juvenile animals as they are the subsequent candidates for selection. This may consequently help speeding up the selection procedure as a result of reduction in the generation interval. Similar gain in accuracy was reported in other studies. For example, in another study [34], the accuracy of estimated breeding values for genotyped female pigs was reported to be 0.22 with the pedigree based approach whereas it ranged from 0.28 to 0.49 with the unified approach depending on the G matrix. Likewise, Christensen and Lund [6] reported an accuracy of 0.66 with the onestep unified approach and 0.35 with the pedigree based approach.
In our study, the accuracies of maternal and direct estimated breeding values for the pedigree based approach (PED_BLUP) with maternal and direct heritability of 0.15 were 0.38 and 0.41 at no correlation and 0.21 and 0.19 at a correlation of −0.46, respectively. In an earlier pedigree based study by Roehe and Kennedy [35], the accuracy of maternal and direct estimated breeding values was reported to be 0.21 (0.21) and 0.38 (0.28) for the case of no correlation and 0.19 (0.18) and 0.31 (0.23) for a negative correlation of −0.5 in female (male) pigs for maternal and direct heritability of 0.05 and 0.1, respectively. These estimates were based on a pedigree based complete animal model with maternal effects. The difference in the accuracies between our study and to that reported by Roehe and Kennedy [35] can be a result of dissimilarities between the two studies such as the construction of the numerator relationship matrix, value of simulated maternal and direct heritability, random selection of the individuals, number of generations simulated, population structure and size. Nonetheless, the comparison of results of the pedigree based approach with the study from Roehe and Kennedy [35] helps to assess and validate the values of accuracy of maternal and direct estimated breeding values obtained in our study. In our study, the accuracies of maternal and direct estimated breeding values were higher for the UNI_BLUP approach as compared to the PED_BLUP approach, but the difference between UNI_BLUP and PED_BLUP approaches were not significant for some cases (Table 4). Thus, in order to achieve maximum gain from implementing the unified approach, a proper investigation into the cost benefits and the relative improvement in genetic gain is required for traits selected solely on the basis of maternal or direct breeding values. Nevertheless, the sum of maternal and direct effects is still the most important criterion for selection and the use of only direct or maternal effects is not helpful for the honey bee.
A complexity associated with the estimation of breeding values for maternally influenced traits is that negative correlation between maternal and direct effects severely impedes the response to selection [35, 36]. Additionally, it leads to a decrease in the total genetic variance resulting in lowered accuracies. As shown in the results, the accuracy of estimated breeding values improved significantly in case of negative correlation with the unified approach as compared to the pedigree based approach. We propose that the extra gain from genomic selection (versus pedigreeBLUP) is larger, when the correlations between direct and maternal effects are negative, compared to scenarios with positive correlations. This is because some of the markers capture maternal (direct) genetic differences, induced by loci without pleiotropic direct (maternal) effects. These parts of the genetic variation are of special interest in case of negative correlations between direct and maternal effects, since they provide, at least partially, opportunity for achieving both maternal and direct genetic gain in the desired direction. In contrast to this, each positive maternal (direct) gain stemming from pleiotropic loci is counterbalanced by a negative direct (maternal) one. Obviously, this kind of distinction does not matter at all in case of positive correlations and is even impossible to be exploited with only pedigree information at hand. So, genomic selection offers special advantages in cases, where maternal and direct effects are negatively correlated through the pleiotropic action of common loci, a phenomenon probably deserving more attention and research also in other species. Additionally, the increase in accuracies can be attributed to the genomic matrix which is able to provide a more precise measure of genetic relatedness. The numerator relationship matrix uses pedigree information to derive the probability of genes to be identical by descent that gives an estimate of the relatedness of individuals. The genomic matrix, in contrast, uses highdensity marker information and thus, can identify genes that are identical by state and may be shared through common ancestor not recorded in the pedigree [34]. Hence, it provides a more accurate measure for the relationship between individuals. It also enables better differentiation among closely related individuals since it captures Mendelian sampling with greater precision. Thus, the use of a marker based relationship matrix in the unified approach greatly helps to improve the accuracy of estimated breeding values for low heritability traits and/or negatively correlated traits, e.g. traits with negatively correlated maternal and direct effects.
Conclusions
To provide a comparison between genetic evaluation methods based on the unified approach and the pedigree based approach, we modelled a complex scenario by taking into consideration characteristics such as varying heritability and correlation between maternal and direct genetic effects, uncertain paternity and other genetic and reproductive peculiarities of the honey bee. To the best of our knowledge, this is the first study that describes the use of molecular marker data for genetic evaluation in honey bees by employing the unified approach. The study provides background knowledge about the simulation of a genomic and a pedigree dataset in honey bees for genetic evaluation, therefore, it can serve as an important framework for future studies. Studies in other species [7, 34] have already optimized the approach with respect to the construction of genomic matrix and computational solving procedures. Thus, additional investigation is needed in future to improve the methodology in the honey bee. The unified approach is a progressive step in the genetic evaluation program of honey bees that will facilitate to reduce the rates of inbreeding, improve the genetic gain and response to selection.
Abbreviations
 BLUP:

Best Linear Unbiased Prediction
 LD:

Linkage Disequilibrium
 MAF:

Minor Allele Frequency
 QTL:

Quantitative Trait Loci
 SNP:

Single Nucleotide Polymorphism
Declarations
Acknowledgements
The work was a part of the project "Marker assisted selection of Varroatolerant honey bees" financially supported by the German Federal Ministry of Food, Agriculture and Consumer Protection (BMELV) [2808HS009].
Authors’ Affiliations
References
 Bienefeld K, Ehrhardt K, Reinhardt F: Genetic evaluation in the honey bee considering queen and worker effects – A BLUPanimal model approach. Apidologie. 2007, 38: 7785. 10.1051/apido:2006050.View Article
 Meuwissen THE, Hayes BJ, Goddard ME: Prediction of total genetic value using genomewide dense marker maps. Genetics. 2001, 157: 18191829.PubMed CentralPubMed
 VanRaden PM, Van Tassell CP, Wiggans GR, Sonstegard TS, Schnabel RD, Taylor JF, Schenkel FS: Invited review: Reliability of genomic predictions for North American Holstein bulls. J Dairy Sci. 2009, 92: 1624. 10.3168/jds.20081514.View ArticlePubMed
 Ostersen T, Christensen OF, Henryon M, Nielsen B, Su G, Madsen P: Deregressed EBV as the response variable yield more reliable genomic predictions than traditional EBV in purebred pigs. Genet Sel Evol. 2011, 43: 3810.1186/129796864338.PubMed CentralView ArticlePubMed
 Legarra A, Aguilar I, Misztal I: A relationship matrix including full pedigree and genomic information. J Dairy Sci. 2009, 92: 46564663. 10.3168/jds.20092061.View ArticlePubMed
 Christensen OF, Lund MS: Genomic prediction when some animals are not genotyped. Genet Sel Evol. 2010, 42: 210.1186/12979686422.PubMed CentralView ArticlePubMed
 Aguilar I, Misztal I, Johnson DL, Legarra A, Tsuruta S, Lawlor TJ: Hot topic: A unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J Dairy Sci. 2010, 93: 743752. 10.3168/jds.20092730.View ArticlePubMed
 Vitezica ZG, Aguilar I, Misztal I, Legarra A: Bias in genomic predictions for populations under selection. Genet Res (Camb). 2011, 93: 357366. 10.1017/S001667231100022X.View Article
 Misztal I, Legarra A, Aguilar I: Computing procedures for genetic evaluation including phenotypic, full pedigree, and genomic information. J Dairy Sci. 2009, 92: 46484655. 10.3168/jds.20092064.View ArticlePubMed
 Cardoso FF, Tempelman RJ: Genetic evaluation of beef cattle accounting for uncertain paternity. Livest Prod Sci. 2004, 89: 109120. 10.1016/j.livprodsci.2004.02.006.View Article
 Lourenco DAL, Misztal I, Wang H, Aguilar I, Tsuruta S: Accuracies with different genomic models for traits with maternal effects [abstract]. J Anim Sci. 2012, 90: s3J Dairy Sci 2012, 95:s2
 Gupta P, Conrad T, Spötter A, Reinsch N, Bienefeld K: Simulating a base population in honey bee for molecular genetic studies. Genet Sel Evol. 2012, 44: 1410.1186/129796864414.PubMed CentralView ArticlePubMed
 Nolan WJ: Bee breeding. United States Department of Agriculture Yearbook. 1937, Washington, DC: United States Department of Agriculture, 13961418.
 Rothenbuhler WC: Genetics and breeding of the honey bee. Annu Rev Entomol. 1958, 3: 161180. 10.1146/annurev.en.03.010158.001113.View Article
 The Honeybee Genome Sequencing Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443: 931949. 10.1038/nature05260.PubMed CentralView Article
 Honey Bee Genome Project. http://www.hgsc.bcm.tmc.edu/content/honeybeegenomeproject,
 Honey Bee Genome Resources. http://www.ncbi.nlm.nih.gov/genome/?term=apis%20mellifera,
 Haldane JBS: The combination of linkage values and the calculation of distances between the loci of linked factors. J Genet. 1919, 8: 299309.View Article
 Beye M, Gattermeier I, Hasselmann M, Gempe T, Schioett M, Baines JF, Schlipalius D, Mougel F, Emore C, Rueppell O, Sirviö A, GuzmánNovoa E, Hunt G, Solignac M, Page RE: Exceptionally high levels of recombination across the honey bee genome. Genome Res. 2006, 16: 13391344. 10.1101/gr.5680406.PubMed CentralView ArticlePubMed
 Spötter A, Gupta P, Nürnberg G, Reinsch N, Bienefeld K: Development of a 44K SNP assay focussing on the analysis of a varroaspecific defense behaviour in honey bees (Apis mellifera carnica). Mol Ecol Resour. 2012, 12: 323332. 10.1111/j.17550998.2011.03106.x.View ArticlePubMed
 Bienefeld K, Pirchner F: Genetic correlations among several colony characters in the honey bee (Hymenoptera: Apidae) taking queen and worker effects into account. Ann Entomol Soc Am. 1991, 84: 324331.View Article
 Falconer DS, Mackay TFC: Introduction to quantitative genetics. 1996, Essex: Longman
 Quaas RL, Pollak EJ: Mixed model methodology for farm and ranch beef cattle testing programs. J Anim Sci. 1980, 51: 12771287.
 VanRaden PM: Efficient Methods to Compute Genomic Predictions. J Dairy Sci. 2008, 91: 44144423. 10.3168/jds.20070980.View ArticlePubMed
 Meuwissen THE, Luan T, Woolliams JA: The unified approach to the use of genomic and pedigree information in genomic evaluations revisited. J Anim Breed Genet. 2011, 128: 429439. 10.1111/j.14390388.2011.00966.x.View ArticlePubMed
 Mrode RA: Linear models for the prediction of animal breeding values. 2005, Wallingford: CABI PublishingView Article
 Bienefeld K, Pirchner F: Heritabilities for several colony traits in the honeybee (Apis mellifera carnica). Apidologie. 1990, 21: 175183. 10.1051/apido:19900302.View Article
 Boecking O, Bienefeld K, Drescher W: Heritability of the Varroaspecific hygienic behaviour in honey bees (Hymenoptera: Apidae). J Anim Breed Genet. 2000, 117: 417424. 10.1046/j.14390388.2000.00271.x.View Article
 CostaMaia FM, de Toledo VAA, Martins EN, LinoLourenço DA, Sereia MJ, de Oliveira CAL, Faquinello P, Halak AL: Estimates of covariance components for hygienic behavior in Africanized honeybees (Apis mellifera). R Bras Zootec. 2011, 40: 19091916. 10.1590/S151635982011000900010.View Article
 Bijma P: Estimating maternal genetic effects in livestock. J Anim Sci. 2006, 84: 800806.PubMed
 Larsgard AG, Olesen I: Genetic parameters for direct and maternal effects on weights and ultrasonic muscle and fat depth of lambs. Livest Prod Sci. 1998, 55: 273278. 10.1016/S03016226(98)001341.View Article
 Splan RK, Cundiff LV, Dikeman ME, Van Vleck LD: Estimates of parameters between direct and maternal genetic effects for weaning weight and direct genetic effects for carcass traits in crossbred cattle. J Anim Sci. 2002, 80: 31073111.PubMed
 Safari E, Fogarty NM, Gilmour AR: A review of genetic parameter estimates for wool, growth, meat and reproduction traits in sheep. Livest Prod Sci. 2005, 92: 271289. 10.1016/j.livprodsci.2004.09.003.View Article
 Forni S, Aguilar I, Misztal I: Different genomic relationship matrices for singlestep analysis using phenotypic, pedigree and genomic information. Genet Sel Evol. 2011, 43: 110.1186/12979686431.PubMed CentralView ArticlePubMed
 Roehe R, Kennedy BW: The influence of maternal effects on accuracy of evaluation of litter size in swine. J Anim Sci. 1993, 71: 23532364.PubMed
 Willham RL: The role of maternal effects in animal breeding: III. Biometrical aspects of maternal effects in animals. J Anim Sci. 1972, 35: 12881293.PubMed
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.