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.

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.

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.

Simulation studies to optimize genomic selection in honey bees.

BACKGROUND: With the completion of a single nucleotide polymorphism (SNP) chip for honey bees, the technical basis of genomic selection is laid. However, for its application in practice, methods to estimate genomic breeding values need to be adapted to the specificities of the genetics and breeding infrastructure of this species. Drone-producing queens (DPQ) are used for mating control, and usually, they head non-phenotyped colonies that will be placed on mating stations. Breeding queens (BQ) head colonies that are intended to be phenotyped and used to produce new queens. Our aim was to evaluate different breeding program designs for the initiation of genomic selection in honey bees. METHODS: Stochastic simulations were conducted to evaluate the quality of the estimated breeding values. We developed a variation of the genomic relationship matrix to include genotypes of DPQ and tested different sizes of the reference population. The results were used to estimate genetic gain in the initial selection cycle of a genomic breeding program. This program was run over six years, and different numbers of genotyped queens per year were considered. Resources could be allocated to increase the reference population, or to perform genomic preselection of BQ and/or DPQ. RESULTS: Including the genotypes of 5000 phenotyped BQ increased the accuracy of predictions of breeding values by up to 173%, depending on the size of the reference population and the trait considered. To initiate a breeding program, genotyping a minimum number of 1000 queens per year is required. In this case, genetic gain was highest when genomic preselection of DPQ was coupled with the genotyping of 10-20% of the phenotyped BQ. For maximum genetic gain per used genotype, more than 2500 genotyped queens per year and preselection of all BQ and DPQ are required. CONCLUSIONS: This study shows that the first priority in a breeding program is to genotype phenotyped BQ to obtain a sufficiently large reference population, which allows successful genomic preselection of queens. To maximize genetic gain, DPQ should be preselected, and their genotypes included in the genomic relationship matrix. We suggest, that the developed methods for genomic prediction are suitable for implementation in genomic honey bee breeding programs.

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.

The importance of controlled mating in honeybee breeding.

BACKGROUND: Controlled mating procedures are widely accepted as a key aspect of successful breeding in almost all animal species. In honeybees, however, controlled mating is hard to achieve. Therefore, there have been several attempts to breed honeybees using free-mated queens. In such breeding schemes, selection occurs only on the maternal path since the drone sires are random samples of the population. The success rates of breeding approaches without controlled mating have so far not been investigated on a theoretical or simulation-based level. METHODS: Stochastic simulation studies were carried out to examine the chances of success in honeybee breeding with and without controlled mating. We investigated the influence of different sizes of breeding populations (500, 1000, 2000 colonies per year) and unselected passive populations (0, 500, 1000, 2000, infinitely many colonies per year) on selection for a maternally (queen) and directly (worker group) influenced trait with moderate ([Formula: see text]) or strong ([Formula: see text]) negative correlation between the two effects. The simulations described 20 years of selection. RESULTS: Our simulations showed a reduction of breeding success between 47 and 99% if mating was not controlled. In the most drastic cases, practically no genetic gain could be generated without controlled mating. We observed that in the trade-off between selection for direct or maternal effects, the absence of mating control leads to a shift in favor of maternal effects. Moreover, we describe the implications of different breeding strategies on the unselected passive population that benefits only indirectly via the transfer of queens or drones from the breeding population. We show that genetic gain in the passive population develops parallel to that of the breeding population. However, we found a genetic lag that became significantly smaller as more breeding queens served as dams of queens in the passive population. CONCLUSIONS: We conclude that even when unwanted admixture of subspecies can be excluded in natural matings, controlled mating is imperative for successful breeding efforts. This is especially highlighted by the strong positive impact that controlled mating in the breeding population has on the unselected passive population.

Differential proliferation rates generate patterns of mechanical tension that orient tissue growth.

Orientation of cell divisions is a key mechanism of tissue morphogenesis. In the growing Drosophila wing imaginal disc epithelium, most of the cell divisions in the central wing pouch are oriented along the proximal-distal (P-D) axis by the Dachsous-Fat-Dachs planar polarity pathway. However, cells at the periphery of the wing pouch instead tend to orient their divisions perpendicular to the P-D axis despite strong Dachs polarization. Here, we show that these circumferential divisions are oriented by circumferential mechanical forces that influence cell shapes and thus orient the mitotic spindle. We propose that this circumferential pattern of force is not generated locally by polarized constriction of individual epithelial cells. Instead, these forces emerge as a global tension pattern that appears to originate from differential rates of cell proliferation within the wing pouch. Accordingly, we show that localized overgrowth is sufficient to induce neighbouring cell stretching and reorientation of cell division. Our results suggest that patterned rates of cell proliferation can influence tissue mechanics and thus determine the orientation of cell divisions and tissue shape.

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.

Novel automated tracking analysis of particles subjected to shear flow: kindlin-3 role in B cells.

Shear flow assays are used to mimic the influence of physiological shear force in diverse situations such as leukocyte rolling and arrest on the vasculature, capture of nanoparticles, and bacterial adhesion. Analysis of such assays usually involves manual counting, is labor-intensive, and is subject to bias. We have developed the Leukotrack program that incorporates a novel (to our knowledge) segmentation routine capable of reliable detection of cells in phase contrast images. The program also automatically tracks rolling cells in addition to those that are more firmly attached and migrating in random directions. We demonstrate its use in the analysis of lymphocyte arrest mediated by one or more active conformations of the integrin LFA-1. Activation of LFA-1 is a multistep process that depends on several proteins including kindlin-3, the protein that is mutated in leukocyte adhesion deficiency-III patients. We find that the very first stage of LFA-1-mediated attaching is unable to proceed in the absence of kindlin-3. Our evidence indicates that kindlin-3-mediated high-affinity LFA-1 controls both the early transient integrin-dependent adhesions in addition to the final stable adhesions made under flow conditions.

Mucosal vaccination against tuberculosis using inert bioparticles.

Needle-free, mucosal immunization is a highly desirable strategy for vaccination against many pathogens, especially those entering through the respiratory mucosa, such as Mycobacterium tuberculosis. Unfortunately, mucosal vaccination against tuberculosis (TB) is impeded by a lack of suitable adjuvants and/or delivery platforms that could induce a protective immune response in humans. Here, we report on a novel biotechnological approach for mucosal vaccination against TB that overcomes some of the current limitations. This is achieved by coating protective TB antigens onto the surface of inert bacterial spores, which are then delivered to the respiratory tract. Our data showed that mice immunized nasally with coated spores developed humoral and cellular immune responses and multifunctional T cells and, most importantly, presented significantly reduced bacterial loads in their lungs and spleens following pathogenic challenge. We conclude that this new vaccine delivery platform merits further development as a mucosal vaccine for TB and possibly also other respiratory pathogens.

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.

Heritability of Apis mellifera recapping behavior and suppressed mite reproduction as resistance traits towards Varroa destructor.

The selection of honeybee strains resistant to the ectoparasitic mite Varroa destructor is generally considered as one of the most sustainable ways of coping with this major bee parasite. Thus, breeding efforts increasingly focus on resistance parameters in addition to common beekeeping traits like honey yield and gentleness. In every breeding effort, the success strongly depends on the quantifiability and heritability of the traits accounted. To find the most suitable traits among the manifold variants to assess Varroa resistance, it is necessary to evaluate how easily a trait can be measured (i.e., testing effort) in relation to the underlying heritability (i.e., expected transfer to the following generation). Various possible selection traits are described as beneficial for colony survival in the presence of Varroa destructor and therefore are measured in breeding stocks around the globe. Two of them in particular, suppressed mite reproduction (SMR, sensu lato any reproductive failure of mother mites) and recapping of already sealed brood cells have recently gained increasing attention among the breeders because they closely resemble resistance mechanisms of some Varroa-surviving honeybee populations. However, it was still unknown whether the genetic background of the trait is sufficient for targeted selection. We therefore investigated the heritabilities and genetic correlations for SMR and REC, distinguishing between recapping of infested cells (RECinf) and all cells (RECall), on an extensive dataset of Buckfast and Carniolan stock in Germany. With an accessible h² of 0.18 and 0.44 for SMR and an accessible h² of 0.44 and 0.40 for RECinf, both traits turned out to be very promising for further selection in the Buckfast and Carnica breeding population, respectively.

Metannogen: annotation of biological reaction networks.

MOTIVATION: Semantic annotations of the biochemical entities constituting a biological reaction network are indispensable to create biologically meaningful networks. They further heighten efficient exchange, reuse and merging of existing models which concern present-day systems biology research more often. Two types of tools for the reconstruction of biological networks currently exist: (i) several sophisticated programs support graphical network editing and visualization. (ii) Data management systems permit reconstruction and curation of huge networks in a team of scientists including data integration, annotation and cross-referencing. We seeked ways to combine the advantages of both approaches. RESULTS: Metannogen, which was previously developed for network reconstruction, has been considerably improved. From now on, Metannogen provides sbml import and annotation of networks created elsewhere. This permits users of other network reconstruction platforms or modeling software to annotate their networks using Metannogen's advanced information management. We implemented word-autocompletion, multipattern highlighting, spell check, brace-expansion and publication management, and improved annotation, cross-referencing and team work requirements. Unspecific enzymes and transporters acting on a spectrum of different substrates are efficiently handled. The network can be exported in sbml format where the annotations are embedded in line with the miriam standard. For more comfort, Metannogen may be tightly coupled with the network editor such that Metannogen becomes an additional view for the focused reaction in the network editor. Finally, Metannogen provides local single user, shared password protected multiuser or public access to the annotation data. AVAILABILITY: Metannogen is available free of charge at: http://www.bioinformatics.org/strap/metannogen/ or http://3d-alignment.eu/metannogen/. CONTACT: christoph.gille@charite.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

ECM degradation in the Drosophila abdominal epidermis initiates tissue growth that ceases with rapid cell-cycle exit.

During development, multicellular organisms undergo stereotypical patterns of tissue growth in space and time. How developmental growth is orchestrated remains unclear, largely due to the difficulty of observing and quantitating this process in a living organism. Drosophila histoblast nests are small clusters of progenitor epithelial cells that undergo extensive growth to give rise to the adult abdominal epidermis and are amenable to live imaging. Our quantitative analysis of histoblast proliferation and tissue mechanics reveals that tissue growth is driven by cell divisions initiated through basal extracellular matrix degradation by matrix metalloproteases secreted by the neighboring larval epidermal cells. Laser ablations and computational simulations show that tissue mechanical tension does not decrease as the histoblasts fill the abdominal epidermal surface. During tissue growth, the histoblasts display oscillatory cell division rates until growth termination occurs through the rapid emergence of G0/G1 arrested cells, rather than a gradual increase in cell-cycle time as observed in other systems such as the Drosophila wing and mouse postnatal epidermis. Different developing tissues can therefore achieve their final size using distinct growth termination strategies. Thus, adult abdominal epidermal development is characterized by changes in the tissue microenvironment and a rapid exit from the cell cycle.

FASIMU: flexible software for flux-balance computation series in large metabolic networks.

BACKGROUND: Flux-balance analysis based on linear optimization is widely used to compute metabolic fluxes in large metabolic networks and gains increasingly importance in network curation and structural analysis. Thus, a computational tool flexible enough to realize a wide variety of FBA algorithms and able to handle batch series of flux-balance optimizations is of great benefit. RESULTS: We present FASIMU, a command line oriented software for the computation of flux distributions using a variety of the most common FBA algorithms, including the first available implementation of (i) weighted flux minimization, (ii) fitness maximization for partially inhibited enzymes, and (iii) of the concentration-based thermodynamic feasibility constraint. It allows batch computation with varying objectives and constraints suited for network pruning, leak analysis, flux-variability analysis, and systematic probing of metabolic objectives for network curation. Input and output supports SBML. FASIMU can work with free (lp_solve and GLPK) or commercial solvers (CPLEX, LINDO). A new plugin (faBiNA) for BiNA allows to conveniently visualize calculated flux distributions. The platform-independent program is an open-source project, freely available under GNU public license at http://www.bioinformatics.org/fasimu including manual, tutorial, and plugins. CONCLUSIONS: We present a flux-balance optimization program whose main merits are the implementation of thermodynamics as a constraint, batch series of computations, free availability of sources, choice on various external solvers, and the flexibility on metabolic objectives and constraints.

HepatoNet1: a comprehensive metabolic reconstruction of the human hepatocyte for the analysis of liver physiology.

We present HepatoNet1, the first reconstruction of a comprehensive metabolic network of the human hepatocyte that is shown to accomplish a large canon of known metabolic liver functions. The network comprises 777 metabolites in six intracellular and two extracellular compartments and 2539 reactions, including 1466 transport reactions. It is based on the manual evaluation of >1500 original scientific research publications to warrant a high-quality evidence-based model. The final network is the result of an iterative process of data compilation and rigorous computational testing of network functionality by means of constraint-based modeling techniques. Taking the hepatic detoxification of ammonia as an example, we show how the availability of nutrients and oxygen may modulate the interplay of various metabolic pathways to allow an efficient response of the liver to perturbations of the homeostasis of blood compounds.

IGERS: inferring Gibbs energy changes of biochemical reactions from reaction similarities.

Mathematical analysis and modeling of biochemical reaction networks requires knowledge of the permitted directionality of reactions and membrane transport processes. This information can be gathered from the standard Gibbs energy changes (DeltaG(0)) of reactions and the concentration ranges of their reactants. Currently, experimental DeltaG(0) values are not available for the vast majority of cellular biochemical processes. We propose what we believe to be a novel computational method to infer the unknown DeltaG(0) value of a reaction from the known DeltaG(0) value of the chemically most similar reaction. The chemical similarity of two arbitrary reactions is measured by the relative number (T) of co-occurring changes in the chemical attributes of their reactants. Testing our method across a validated reference set of 173 biochemical reactions with experimentally determined DeltaG(0) values, we found that a minimum reaction similarity of T = 0.6 is required to infer DeltaG(0) values with an error of <10 kJ/mol. Applying this criterion, our method allows us to assign DeltaG(0) values to 458 additional reactions of the BioPath database. We believe our approach permits us to minimize the number of DeltaG(0) measurements required for a full coverage of a given reaction network with reliable DeltaG(0) values.

Superimpose: a 3D structural superposition server.

The Superimposé webserver performs structural similarity searches with a preference towards 3D structure-based methods. Similarities can be detected between small molecules (e.g. drugs), parts of large structures (e.g. binding sites of proteins) and entire proteins. For this purpose, a number of algorithms were implemented and various databases are provided. Superimposé assists the user regarding the selection of a suitable combination of algorithm and database. After the computation on our server infrastructure, a visual assessment of the results is provided. The structure-based in silico screening for similar drug-like compounds enables the detection of scaffold-hoppers with putatively similar effects. The possibility to find similar binding sites can be of special interest in the functional analysis of proteins. The search for structurally similar proteins allows the detection of similar folds with different backbone topology. The Superimposé server is available at: http://bioinformatics.charite.de/superimpose.