Long-Term Evaluation of Breeding Scheme Alternatives for Endangered Honeybee Subspecies

Comparison of pooled semen insemination and single colony insemination as sustainable honeybee breeding strategies

Instrumental insemination of honeybees allows for two opposing breeding strategies. In single colony insemination (SCI), all drones to inseminate a queen are taken from one colony. In pooled semen insemination (PSI), sperm of many genetically diverse drones is mixed and queens are fertilized from the resulting drone pool. While SCI allows for maximum pedigree control, proponents of PSI claim to reduce inbreeding and maintain genetic variance. Using stochastic simulation studies, we compared genetic progress and inbreeding rates in small honeybee populations under SCI and PSI. Four different selection criteria were covered: estimated breeding values (EBV), phenotypes, true breeding values (TBV) and random selection. Under EBV-based truncation selection, SCI yielded 9.0% to 44.4% higher genetic gain than PSI, but had vastly increased inbreeding rates. Under phenotypical or TBV selection, the gap between SCI and PSI in terms of genetic progress narrowed. Throughout, PSI yielded lower inbreeding rates than SCI, but the differences were only substantial under EBV truncation selection. As a result, PSI did not appear as a viable breeding strategy owing to its incompatibility with modern methods of genetic evaluation. Instead, SCI is to be preferred but instead of strict truncation selection, strategies to avoid inbreeding need to be installed.

The Potential of Instrumental Insemination for Sustainable Honeybee Breeding

Mating control is crucial in honeybee breeding and commonly guaranteed by bringing virgin queens to isolated mating stations (IMS) for their nuptial flights. However, most breeding programs struggle to provide sufficiently many IMS. Research institutions routinely perform instrumental insemination of honeybees, but its potential to substitute IMS in breeding programs has not been sufficiently studied. We performed stochastic simulations to compare instrumental insemination strategies and mating on IMS in terms of genetic progress and inbreeding development. We focused on the role of paternal generation intervals, which can be shortened to two years with instrumental insemination in comparison to three years when using IMS. After 70 years, instrumental insemination yielded up to 42% higher genetic gain than IMS strategies-particularly with few available mating sites. Inbreeding rates with instrumental insemination and IMS were comparable. When the paternal generation interval in instrumental insemination was stretched to three years, the number of drone producers required for sustainable breeding was reduced substantially. In contrast, when shortening the interval to two years, it yielded the highest generational inbreeding rates (up to 2.28%). Overall, instrumental insemination with drones from a single colony appears as a viable strategy for honeybee breeding and a promising alternative to IMS.

New insights into the criteria of functional heterozygosity of the Apis mellifera complementary sex determining gene–Discovery of a functional allele pair differing by a single amino acid

The complementary sex determiner (csd) gene is responsible for controlling the sex-determination molecular switch in western honey bees (Apis mellifera): bees that are heterozygous for csd develop into females, whereas bees that are hemizygous or homozygous develop into males. The homozygous diploid males are destroyed at an early stage of their development. It has been proposed that the minimal number of amino acid differences between two csd alleles needed to fully determine femaleness is five and it has also been shown that smaller differences may result in forming an evolutionary intermediate that is not fully capable of female determination, but has increased fitness compared to the homozygous genotype. In this study, we have implemented a terminal restriction length polymorphism-based method of identifying and distinguishing paternal alleles in a given bee colony and assigning them to a particular maternal allele in order to gather information on large number of functional csd pairs and also to identify, to some extent, genotypes that are underrepresented or absent in bee colonies. The main finding of this study is the identification of a fully functional genotype consisting of csd alleles that differed from each other by a one amino acid position. The individuals carrying this genotype expressed only female-specific transcripts of feminizer and double-sex genes. By comparing the sequences differences between the csd pair identified in our study with those described earlier, we conclude that functional heterozygosity of the csd gene is dependent not only on the number of the amino acid differences but also on the sequence context and position of the change. The discovery of a functional allele pair differing by a single amino acid also implies that the generation of a new csd specificity may also occur during a single mutation step with no need for evolutionary intermediates accumulating further mutations.

Consequences of incorrect genetic parameter estimates for single-trait and multi-trait genetic evaluations in honeybees

Genetic and residual variances of traits are important input parameters for best linear unbiased prediction (BLUP) breeding value estimation. In honeybees, estimates of these variances are often associated with large standard errors, entailing a risk to perform genetic evaluations under wrong premises. The consequences hereof have not been sufficiently studied. In particular, there are no adequate investigations on this topic accounting for multi-trait selection or genetic peculiarities of the honeybee. We performed simulation studies and explored the consequences of selection for honeybee populations with a broad range of true and assumed genetic parameters. We found that in single-trait evaluations, the response to selection was barely compromised by assuming erroneous parameters, so that reductions in genetic progress after 20 years never exceeded 21%. Phenotypic selection appeared inferior to BLUP selection, particularly under low heritabilities. Parameter choices for genetic evaluation had great effects on inbreeding development. By wrongly assuming high heritabilities, inbreeding rates were reduced by up to 74%. When parallel selection was performed for two traits, the right choice of genetic parameters appeared considerably more crucial as several incorrect premises yielded inadvertent negative selection for one of the traits. This phenomenon occurred in multiple constellations in which the selection traits expressed a negative genetic correlation. It was not reflected in the estimated breeding values. Our results indicate that breeding efforts heavily rely on detailed knowledge on genetic parameters, particularly when multi-trait selection is performed. Thus, considerable effort should be invested into precise parameter estimations.

Influence of model selection and data structure on the estimation of genetic parameters in honeybee populations

Estimating genetic parameters of quantitative traits is a prerequisite for animal breeding. In honeybees, the genetic variance separates into queen and worker effects. However, under data paucity, parameter estimations that account for this peculiarity often yield implausible results. Consequently, simplified models that attribute all genetic contributions to either the queen (queen model) or the workers (worker model) are often used to estimate variance components in honeybees. However, the causes for estimations with the complete model (colony model) to fail and the consequences of simplified models for variance estimates are little understood. We newly developed the necessary theory to compare parameter estimates that were achieved by the colony model with those of the queen and worker models. Furthermore, we performed computer simulations to quantify the influence of model choice, estimation algorithm, true genetic parameters, rates of controlled mating, apiary sizes, and phenotype data completeness on the success of genetic parameter estimations. We found that successful estimations with the colony model were only possible if at least some of the queens mated controlled on mating stations. In that case, estimates were largely unbiased if more than 20% of the colonies had phenotype records. The simplified queen and worker models proved more stable and yielded plausible parameter estimates for almost all settings. Results obtained from these models were unbiased when mating was uncontrolled, but with controlled mating, the simplified models consistently overestimated heritabilities. This study elucidates the requirements for variance component estimation in honeybees and provides the theoretical groundwork for simplified honeybee models.

A simulation study of a honeybee breeding scheme accounting for polyandry, direct and maternal effects on colony performance

BACKGROUND

Efficient breeding programs are difficult to implement in honeybees due to their biological specificities (polyandry and haplo-diploidy) and complexity of the traits of interest, with performances being measured at the colony scale and resulting from the joint effects of tens of thousands of workers (called direct effects) and of the queen (called maternal effects). We implemented a Monte Carlo simulation program of a breeding plan designed specifically for Apis mellifera's populations to assess the impact of polyandry versus monoandry on colony performance, inbreeding level and genetic gain depending on the individual selection strategy considered, i.e. complete mass selection or within-family (maternal lines) selection. We simulated several scenarios with different parameter setups by varying initial genetic variances and correlations between direct and maternal effects, the selection strategy and the polyandry level. Selection was performed on colony phenotypes.

RESULTS

All scenarios showed strong increases in direct breeding values of queens after 20 years of selection. Monoandry led to significantly higher direct than maternal genetic gains, especially when a negative correlation between direct and maternal effects was simulated. However, the relative increase in these genetic gains depended also on their initial genetic variability and on the selection strategy. When polyandry was simulated, the results were very similar with either 8 or 16 drones mated to each queen. Across scenarios, polyandrous mating resulted in equivalent or higher gains in performance than monoandrous mating, but with considerably lower inbreeding rates. Mass selection conferred a ~ 20% increase in performance compared to within-family selection, but was also accompanied by a strong increase in inbreeding levels (25 to 50% higher).

CONCLUSIONS

Our study is the first to compare the long-term effects of polyandrous versus monoandrous mating in honeybee breeding. The latter is an emergent strategy to improve specific traits, such as resistance to varroa, which can be difficult or expensive to phenotype. However, if used during several generations in a closed population, monoandrous mating increases the inbreeding level of queens much more than polyandrous mating, which is a strong limitation of this strategy.

Short-term effects of controlled mating and selection on the genetic variance of honeybee populations

Directional selection in a population yields reduced genetic variance due to the Bulmer effect. While this effect has been thoroughly investigated in mammals, it is poorly studied in social insects with biological peculiarities such as haplo-diploidy or the collective expression of traits. In addition to the natural adaptation to climate change, parasites, and pesticides, honeybees increasingly experience artificial selection pressure through modern breeding programs. Besides selection, many honeybee breeding schemes introduce controlled mating. We investigated which individual effects selection and controlled mating have on genetic variance. We derived formulas to describe short-term changes of genetic variance in honeybee populations and conducted computer simulations to confirm them. Thereby, we found that the changes in genetic variance depend on whether the variance is measured between queens (inheritance criterion), worker groups (selection criterion), or both (performance criterion). All three criteria showed reduced genetic variance under selection. In the selection and performance criteria, our formulas and simulations showed an increased genetic variance through controlled mating. This newly described effect counterbalanced and occasionally outweighed the Bulmer effect. It could not be observed in the inheritance criterion. A good understanding of the different notions of genetic variance in honeybees, therefore, appears crucial to interpreting population parameters correctly.

A theoretical derivation of response to selection with and without controlled mating in honeybees

BACKGROUND

In recent years, the breeding of honeybees has gained significant scientific interest, and numerous theoretical and practical improvements have been made regarding the collection and processing of their performance data. It is now known that the selection of high-quality drone material is crucial for mid to long-term breeding success. However, there has been no conclusive mathematical theory to explain these findings.

METHODS

We derived mathematical formulas to describe the response to selection of a breeding population and an unselected passive population of honeybees that benefits indirectly from genetic improvement in the breeding population via migration. This was done under the assumption of either controlled or uncontrolled mating of queens in the breeding population.

RESULTS

Our model equations confirm what has been observed in simulation studies. In particular, we have proven that the breeding population and the passive population will show parallel genetic gain after some years and we were able to assess the responses to selection for different breeding strategies. Thus, we confirmed the crucial importance of controlled mating for successful honeybee breeding. When compared with data from simulation studies, the derived formulas showed high coefficients of determination [Formula: see text] in cases where many passive queens had dams from the breeding population. For self-sufficient passive populations, the coefficients of determination were lower ([Formula: see text]) if the breeding population was under controlled mating. This can be explained by the limited simulated time-frame and lower convergence rates.

CONCLUSION

The presented theoretical derivations allow extrapolation of honeybee-specific simulation results for breeding programs to a wide range of population parameters. Furthermore, they provide general insights into the genetic dynamics of interdependent populations, not only for honeybees but also in a broader context.

Substantial Genetic Progress in the International Apis mellifera carnica Population Since the Implementation of Genetic Evaluation

Simple Summary

The Apis mellifera carnica subspecies of the honeybee is known for its gentleness and good honey yield. In the early 20th century, systematic breeding efforts began. Breeding progress was slow before the introduction of modern techniques of genetic evaluation in the mid 1990s. Here, the results of the official breeding value estimation in BeeBreed.eu are analyzed. From about 2000 onward, breeding progress accelerated. The result is a considerable gain in honey yield and desirable properties without increased inbreeding coefficients. The success of A. m. carnica breeding shows the potential of genetic evaluation.

Abstract

The Apis mellifera carnica subspecies of the honeybee has long been praised for its gentleness and good honey yield before systematic breeding efforts began in the early 20th century. However, before the introduction of modern techniques of genetic evaluation (best linear unbiased prediction, BLUP) and a computerized data management in the mid 1990s, genetic progress was slow. Here, the results of the official breeding value estimation in BeeBreed.eu are analyzed to characterize breeding progress and inbreeding. From about the year 2000 onward, the genetic progression accelerated and resulted in a considerable gain in honey yield and desirable properties without increased inbreeding coefficients. The prognostic quality of breeding values is demonstrated by a retrospective analysis. The success of A. m. carnica breeding shows the potential of BLUP-based breeding values and serves as an example for a large-scale breeding program.