Linkage map of the honey bee, Apis mellifera, based on RAPD markers

Causes and consequences of a complex recombinational landscape in the ant Cardiocondyla obscurior

Eusocial Hymenoptera have the highest recombination rates among all multicellular animals studied so far, but it is unclear why this is and how this affects the biology of individual species. A high-resolution linkage map for the ant Cardiocondyla obscurior corroborates genome-wide high recombination rates reported for ants (8.1 cM/Mb). However, recombination is locally suppressed in regions that are enriched with TEs, that have strong haplotype divergence, or that show signatures of epistatic selection in C. obscurior The results do not support the hypotheses that high recombination rates are linked to phenotypic plasticity or to modulating selection efficiency. Instead, genetic diversity and the frequency of structural variants correlate positively with local recombination rates, potentially compensating for the low levels of genetic variation expected in haplodiploid social Hymenoptera with low effective population size. Ultimately, the data show that recombination contributes to within-population polymorphism and to the divergence of the lineages within C. obscurior.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes

BACKGROUND

The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses.

RESULTS

Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants.

CONCLUSIONS

Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

Tyramine and its Amtyr1 receptor modulate attention in honey bees (Apis mellifera)

Animals must learn to ignore stimuli that are irrelevant to survival and attend to ones that enhance survival. When a stimulus regularly fails to be associated with an important consequence, subsequent excitatory learning about that stimulus can be delayed, which is a form of nonassociative conditioning called 'latent inhibition'. Honey bees show latent inhibition toward an odor they have experienced without association with food reinforcement. Moreover, individual honey bees from the same colony differ in the degree to which they show latent inhibition, and these individual differences have a genetic basis. To investigate the mechanisms that underly individual differences in latent inhibition, we selected two honey bee lines for high and low latent inhibition, respectively. We crossed those lines and mapped a Quantitative Trait Locus for latent inhibition to a region of the genome that contains the tyramine receptor gene Amtyr1 [We use Amtyr1 to denote the gene and AmTYR1 the receptor throughout the text.]. We then show that disruption of Amtyr1 signaling either pharmacologically or through RNAi qualitatively changes the expression of latent inhibition but has little or slight effects on appetitive conditioning, and these results suggest that AmTYR1 modulates inhibitory processing in the CNS. Electrophysiological recordings from the brain during pharmacological blockade are consistent with a model that AmTYR1 indirectly regulates at inhibitory synapses in the CNS. Our results therefore identify a distinct Amtyr1-based modulatory pathway for this type of nonassociative learning, and we propose a model for how Amtyr1 acts as a gain control to modulate hebbian plasticity at defined synapses in the CNS. We have shown elsewhere how this modulation also underlies potentially adaptive intracolonial learning differences among individuals that benefit colony survival. Finally, our neural model suggests a mechanism for the broad pleiotropy this gene has on several different behaviors.

A Global Survey of the Full-Length Transcriptome of Apis mellifera by Single-Molecule Long-Read Sequencing

As important pollinators, honey bees play a crucial role in both maintaining the ecological balance and providing products for humans. Although several versions of the western honey bee genome have already been published, its transcriptome information still needs to be refined. In this study, PacBio single-molecule sequencing technology was used to sequence the full-length transcriptome of mixed samples from many developmental time points and tissues of A. mellifera queens, workers and drones. A total of 116,535 transcripts corresponding to 30,045 genes were obtained. Of these, 92,477 transcripts were annotated. Compared to the annotated genes and transcripts on the reference genome, 18,915 gene loci and 96,176 transcripts were newly identified. From these transcripts, 136,554 alternative splicing (AS) events, 23,376 alternative polyadenylation (APA) sites and 21,813 lncRNAs were detected. In addition, based on the full-length transcripts, we identified many differentially expressed transcripts (DETs) between queen, worker and drone. Our results provide a complete set of reference transcripts for A. mellifera that dramatically expand our understanding of the complexity and diversity of the honey bee transcriptome.

Assemblies of the genomes of parasitic wasps using meta-assembly and scaffolding with genetic linkage

Safe, effective biological-control introductions against invasive pests depend on narrowly host-specific natural enemies with the ability to adapt to a changing environment. As part of a project on the genetic architectures of these traits, we assembled and annotated the genomes of two aphid parasitoids, Aphelinus atriplicis and Aphelinus certus. We report here several assemblies of A. atriplicis made with Illumina and PacBio data, which we combined into a meta-assembly. We scaffolded the meta-assembly with markers from a genetic map of hybrids between A. atriplicis and A. certus. We used this genetic-linkage scaffolded (GLS) assembly of A. atriplicis to scaffold a de novo assembly of A. certus. The de novo assemblies of A. atriplicis differed in contiguity, and the meta-assembly of these assemblies was more contiguous than the best de novo assembly. Scaffolding with genetic-linkage data allowed chromosomal-level assembly of the A. atriplicis genome and scaffolding a de novo assembly of A. certus with this GLS assembly, greatly increased the contiguity of the A. certus assembly to the point where it was also at the chromosomal-level. However, completeness of the A. atriplicis assembly, as measured by percent complete, single-copy BUSCO hymenopteran genes, varied little among de novo assemblies and was not increased by meta-assembly or genetic scaffolding. Furthermore, the greater contiguity of the meta-assembly and GLS assembly had little or no effect on the numbers of genes identified, the proportions with homologs or functional annotations. Increased contiguity of the A. certus assembly provided modest improvement in assembly completeness, as measured by percent complete, single-copy BUSCO hymenopteran genes. The total genic sequence increased, and while the number of genes declined, gene length increased, which together suggest greater accuracy of gene models. More contiguous assemblies provide uses other than gene annotation, for example, identifying the genes associated with quantitative trait loci and understanding of chromosomal rearrangements associated with speciation.

Societies to genes: can we get there from here?

Understanding the organization and evolution of social complexity is a major task because it requires building an understanding of mechanisms operating at different levels of biological organization from genes to social interactions. I discuss here, a unique forward genetic approach spanning more than 30 years beginning with human-assisted colony-level selection for a single social trait, the amount of pollen honey bees (Apis mellifera L.) store. The goal was to understand a complex social trait from the social phenotype to genes responsible for observed trait variation. The approach combined the results of colony-level selection with detailed studies of individual behavior and physiology resulting in a mapped, integrated phenotypic architecture composed of correlative relationships between traits spanning anatomy, physiology, sensory response systems, and individual behavior that affect individual foraging decisions. Colony-level selection reverse engineered the architecture of an integrated phenotype of individuals resulting in changes in the social trait. Quantitative trait locus (QTL) studies combined with an exceptionally high recombination rate (60 kb/cM), and a phenotypic map, provided a genotype-phenotype map of high complexity demonstrating broad QTL pleiotropy, epistasis, and epistatic pleiotropy suggesting that gene pleiotropy or tight linkage of genes within QTL integrated the phenotype. Gene expression and knockdown of identified positional candidates revealed genes affecting foraging behavior and confirmed one pleiotropic gene, a tyramine receptor, as a target for colony-level selection that was under selection in two different tissues in two different life stages. The approach presented here has resulted in a comprehensive understanding of the structure and evolution of honey bee social organization.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Possible Epigenetic Origin of a Recurrent Gynandromorph Pattern in Megachile Wild Bees

Simple Summary

Gynandromorphs, i.e., individuals with a mix of male and female body parts, are known for many species of insects and other animals with separate sexes. This anomaly is generally regarded as the result of localized genetic mutations in sex-determining genes. We analyzed the specific mix of male and female characters in naturally occurring gynandromorphs of 21 species of the wild bee genus Megachile and found a recurrent pattern. Based on the regularity of this pattern, and the current knowledge on sex determination and sex differentiation in the relatively closely-related honey bee, we argue that the origin of these composite phenotypes is possibly epigenetic, rather than genetic, i.e., produced by some defects in the maintenance of the regulatory signals that control sex differentiation at the level of single cell lineages, rather than triggered by genetic mutations.

Abstract

Gynandromorphs, i.e., individuals with a mix of male and female traits, are common in the wild bees of the genus Megachile (Hymenoptera, Apoidea). We described new transverse gynandromorphs in Megachile pilidens Alfkeen, 1924 and analyze the spatial distribution of body parts with male vs. female phenotype hitherto recorded in the transverse gynandromorphs of the genus Megachile. We identified 10 different arrangements, nine of which are minor variants of a very general pattern, with a combination of male and female traits largely shared by the gynandromorphs recorded in 20 out of 21 Megachile species in our dataset. Based on the recurrence of the same gynandromorph pattern, the current knowledge on sex determination and sex differentiation in the honey bee, and the results of recent gene-knockdown experiments in these insects, we suggest that these composite phenotypes are possibly epigenetic, rather than genetic, mosaics, with individual body parts of either male or female phenotype according to the locally expressed product of the alternative splicing of sex-determining gene transcripts.

Tyramine and its receptor TYR1 linked behavior QTL to reproductive physiology in honey bee workers (Apis mellifera)

Honey bees (Apis mellifera) provide an excellent model for studying how complex social behavior evolves and is regulated. Social behavioral traits such as the division of labor have been mapped to specific genomic regions in quantitative trait locus (QTL) studies. However, relating genomic mapping to gene function and regulatory mechanism remains a big challenge for geneticists. In honey bee workers, division of labor is known to be regulated by reproductive physiology, but the genetic basis of this regulation remains unknown. In this case, QTL studies have identified tyramine receptor 1 (TYR1) as a candidate gene in region pln2, which is associated with multiple worker social traits and reproductive anatomy. Tyramine (TA), a neurotransmitter, regulates physiology and behavior in diverse insect species including honey bees. Here, we examine directly the effects of TYR1 and TA on worker reproductive physiology, including ovariole number, ovary function and the production of vitellogenin (VG, an egg yolk precursor). First, we used a pharmacology approach to demonstrate that TA affects ovariole number during worker larval development and increases ovary maturation during the adult stage. Second, we used a gene knockdown approach to show that TYR1 regulates vg transcription in adult workers. Finally, we estimated correlations in gene expression and propose that TYR1 may regulate vg transcription by coordinating hormonal and nutritional signals. Taken together, our results suggest TYR1 and TA play important roles in regulating worker reproductive physiology, which in turn regulates social behavior. Our study exemplifies a successful forward-genetic strategy going from QTL mapping to gene function.

Experimental psychology meets behavioral ecology: what laboratory studies of learning polymorphisms mean for learning under natural conditions, and vice versa

Behavior genetics, and specifically the study of learning and memory, has benefitted immensely from the development of powerful forward- and reverse-genetic methods for investigating the relationships between genes and behavior. Application of these methods in controlled laboratory settings has led to insights into gene-behavior relationships. In this perspective article, we argue that the field is now poised to make significant inroads into understanding the adaptive value of heritable variation in behavior in natural populations. Studies of natural variation with several species, in particular, are now in a position to complement laboratory studies of mechanisms, and sometimes this work can lead to counterintuitive insights into the mechanism of gene action on behavior. We make this case using a recent example from work with the honey bee, Apis mellifera.

Local Variation in Recombination Rates of the Honey Bee (Apis mellifera) Genome among Samples from Six Disparate Populations

Meiotic recombination is an essential component of eukaryotic sexual reproduction but its frequency varies within and between genomes. Although it is well-established that honey bees have a high recombination rate with about 20 cM/Mbp, the proximate and ultimate causes of this exceptional rate are poorly understood. Here, we describe six linkage maps of the Western Honey Bee Apis mellifera that were produced with consistent methodology from samples from distinct parts of the species' near global distribution. We compared the genome-wide rates and distribution of meiotic crossovers among the six maps and found considerable differences. Overall similarity of local recombination rates among our samples was unrelated to geographic or phylogenetic distance of the populations that our samples were derived from. However, the limited sampling constrains the interpretation of our results because it is unclear how representative these samples are. In contrast to previous studies, we found only in two datasets a significant relation between local recombination rate and GC content. Focusing on regions of particularly increased or decreased recombination in specific maps, we identified several enriched gene ontologies in these regions and speculate about their local adaptive relevance. These data are contributing to an increasing comparative effort to gain an understanding of the intra-specific variability of recombination rates and their evolutionary role in honey bees and other social insects.

A genetic switch for worker nutrition-mediated traits in honeybees

Highly social insects are characterized by caste dimorphism, with distinct size differences of reproductive organs between fertile queens and the more or less sterile workers. An abundance of nutrition or instruction via diet-specific compounds has been proposed as explanations for the nutrition-driven queen and worker polyphenism. Here, we further explored these models in the honeybee (Apis mellifera) using worker nutrition rearing and a novel mutational screening approach using the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) method. The worker nutrition-driven size reduction of reproductive organs was restricted to the female sex, suggesting input from the sex determination pathway. Genetic screens on the sex determination genes in genetic females for size polyphenism revealed that doublesex (dsx) mutants display size-reduced reproductive organs irrespective of the sexual morphology of the organ tissue. In contrast, feminizer (fem) mutants lost the response to worker nutrition-driven size control. The first morphological worker mutants in honeybees demonstrate that the response to nutrition relies on a genetic program that is switched “ON” by the fem gene. Thus, the genetic instruction provided by the fem gene provides an entry point to genetically dissect the underlying processes that implement the size polyphenism.

In honeybees, nutrition drives dimorphic size development of reproductive organs in fertile queens and sterile workers. A study using the first induced morphological mutants in honeybees demonstrates that this developmental plasticity requires a genetic program that is switched on by the “feminizer” gene.

In honeybees, nutrition drives dimorphic size development of reproductive organs in fertile queens and sterile workers. The first induced morphological mutants in honeybees demonstrate that this developmental plasticity requires a genetic program that is switched “ON” by the feminizer (fem) gene.

A New Metazoan Recombination Rate Record and Consistently High Recombination Rates in the Honey Bee Genus Apis Accompanied by Frequent Inversions but Not Translocations

Western honey bees (Apis mellifera) far exceed the commonly observed 1–2 meiotic recombination events per chromosome and exhibit the highest Metazoan recombination rate (20 cM/Mb) described thus far. However, the reasons for this exceptional rate of recombination are not sufficiently understood. In a comparative study, we report on the newly constructed genomic linkage maps of Apis florea and Apis dorsata that represent the two honey bee lineages without recombination rate estimates so far. Each linkage map was generated de novo, based on SNP genotypes of haploid male offspring of a single female. The A. florea map spans 4,782 cM with 1,279 markers in 16 linkage groups. The A. dorsata map is 5,762 cM long and contains 1,189 markers in 16 linkage groups. Respectively, these map sizes result in average recombination rate estimates of 20.8 and 25.1 cM/Mb. Synteny analyses indicate that frequent intra-chromosomal rearrangements but no translocations among chromosomes accompany the high rates of recombination during the independent evolution of the three major honey bee lineages. Our results imply a common cause for the evolution of very high recombination rates in Apis. Our findings also suggest that frequent homologous recombination during meiosis might increase ectopic recombination and rearrangements within but not between chromosomes. It remains to be investigated whether the resulting inversions may have been important in the evolutionary differentiation between honey bee species.

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M-lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C-lineage). Effective conservation actions require reliable molecular tools to identify pure-bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q-values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.

The complete mitochondrial genome of the Cape honey bee, Apis mellifera capensis Esch. (Insecta: hymenoptera: apidae)

We characterized the complete mitogenome sequence of the Cape honey bee, Apis mellifera capensis, from South Africa. The circle genome is 16,470 bp in length, with the base composition of 43.2% A, 9.6% C, 5.6% G, and 41.5% T. The assembled mitogenome has 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNA genes, and one control region. All protein-coding genes are initiated by ATT, ATC, ATG or ATA codons and are terminated by the typical stop codon TAA. The heavy strand encodes four protein-coding genes, eight tRNAs, and two rRNAs. The light strand encodes nine protein-coding genes and 14 tRNAs. The complete mitogenome sequence of A.m. capensis is identical to the gene arrangement found in other A. mellifera mitogenomes and it provides essential and important DNA molecular data for further phylogenetic and evolutionary analysis of members of the genus Apis.

Kimball R, L. Braun E, Ellis JD. The complete mitochondrial genome of the hybrid honey bee, Apis mellifera capensis × Apis mellifera scutellata, from South Africa

We characterized the complete mitogenome sequence of the South African hybrid honey bee Apis mellifera capensis × Apis mellifera scutellata using genome skimming. The mitochondrial genome was a circular molecule 16,340 bp in length with a gene organization identical to that of the other A. mellifera mitogenomes. The base composition is 43.2% A, 9.7% C, 5.6% G, and 41.5% T, with an A + T content of 84.7%. The mitogenome had 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs genes, and one control region. All PCGs were initiated by ATT, ATG, ATA, and ATC codons and were terminated by a TAA stop codon. The heavy strand encodes four PCGs, eight tRNAs, and two rRNAs. The light strand encodes nine PCGS and 14 tRNAs. A phylogenetic analysis of the PCGs reveals a close relationship between this hybrid honey bee and other Apis spp.

Transcriptome profiling of Diachasmimorpha longicaudata towards useful molecular tools for population management

Background

Diachasmimorpha longicaudata (Hymenoptera: Braconidae) is a solitary parasitoid of Tephritidae (Diptera) fruit flies of economic importance currently being mass-reared in bio-factories and successfully used worldwide. A peculiar biological aspect of Hymenoptera is its haplo-diploid life cycle, where females (diploid) develop from fertilized eggs and males (haploid) from unfertilized eggs. Diploid males were described in many species and recently evidenced in D. longicaudata by mean of inbreeding studies. Sex determination in this parasitoid is based on the Complementary Sex Determination (CSD) system, with alleles from at least one locus involved in early steps of this pathway. Since limited information is available about genetics of this parasitoid species, a deeper analysis on D. longicaudata’s genomics is required to provide molecular tools for achieving a more cost effective production under artificial rearing conditions.

Results

We report here the first transcriptome analysis of male-larvae, adult females and adult males of D. longicaudata using 454-pyrosequencing. A total of 469766 reads were analyzed and 8483 high-quality isotigs were assembled. After functional annotation, a total of 51686 unigenes were produced, from which, 7021 isotigs and 20227 singletons had at least one BLAST hit against the NCBI non-redundant protein database. A preliminary comparison of adult female and male evidenced that 98 transcripts showed differential expression profiles, with at least a 10-fold difference. Among the functionally annotated transcripts we detected four sequences potentially involved in sex determination and three homologues to two known genes involved in the sex determination cascade. Finally, a total of 4674SimpleSequence Repeats (SSRs) were in silico identified and characterized.

Conclusion

The information obtained here will significantly contribute to the development of D. longicaudata functional genomics, genetics and population-based genome studies. Thousands of new microsatellite markers were identified as toolkits for population genetics analysis. The transcriptome characterized here is the starting point to elucidate the molecular bases of the sex determination mechanism in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2759-2) contains supplementary material, which is available to authorized users.

Genome-Wide Association Study of a Varroa-Specific Defense Behavior in Honeybees (Apis mellifera)

Honey bees are exposed to many damaging pathogens and parasites. The most devastating is Varroa destructor, which mainly affects the brood. A promising approach for preventing its spread is to breed Varroa-resistant honey bees. One trait that has been shown to provide significant resistance against the Varroa mite is hygienic behavior, which is a behavioral response of honeybee workers to brood diseases in general. Here, we report the use of an Affymetrix 44K SNP array to analyze SNPs associated with detection and uncapping of Varroa-parasitized brood by individual worker bees (Apis mellifera). For this study, 22 000 individually labeled bees were video-monitored and a sample of 122 cases and 122 controls was collected and analyzed to determine the dependence/independence of SNP genotypes from hygienic and nonhygienic behavior on a genome-wide scale. After false-discovery rate correction of the P values, 6 SNP markers had highly significant associations with the trait investigated (α < 0.01). Inspection of the genomic regions around these SNPs led to the discovery of putative candidate genes.

Population-level consequences of complementary sex determination in a solitary parasitoid

Background

Sex determination mechanisms are known to be evolutionarily labile but the factors driving transitions in sex determination mechanisms are poorly understood. All insects of the Hymenoptera are haplodiploid, with males normally developing from unfertilized haploid eggs. Under complementary sex determination (CSD), diploid males can be produced from fertilized eggs that are homozygous at the sex locus. Diploid males have near-zero fitness and thus represent a genetic load, which is especially severe under inbreeding. Here, we study mating structure and sex determination in the parasitoid Cotesia vestalis to investigate what may have driven the evolution of two complementary sex determination loci in this species.

Results

We genotyped Cotesia vestalis females collected from eight fields in four townships in Western Taiwan. 98 SNP markers were developed by aligning Illumina sequence reads of pooled DNA of eight different females against a de novo assembled genome of C. vestalis. This proved to be an efficient method for this non-model species and provides a resource for future use in related species. We found significant genetic differentiation within the sampled population but variation could not be attributed to sampling locations by AMOVA. Non-random mating was detected, with 8.1% of matings between siblings. Diploid males, detected by flow cytometry, were produced at a rate of 1.4% among diploids.

Conclusions

We think that the low rate of diploid male production is best explained by a CSD system with two independent sex loci, supporting laboratory findings on the same species. Fitness costs of diploid males in C. vestalis are high because diploid males can mate with females and produce infertile triploid offspring. This severe fitness cost of diploid males combined with non-random mating may have resulted in evolution from single locus CSD to CSD with two independent loci.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0340-2) contains supplementary material, which is available to authorized users.

A selective sweep in a Varroa destructor resistant honeybee (Apis mellifera) population

The mite Varroa destructor is one of the most dangerous parasites of the Western honeybee (Apis mellifera) causing enormous colony losses worldwide. Various chemical treatments for the control of the Varroa mite are currently in use, which, however, lead to residues in bee products and often to resistance in mites. This facilitated the exploration of alternative treatment methods and breeding for mite resistant honeybees has been in focus for breeders in many parts of the world with variable results. Another approach has been applied to a honeybee population on Gotland (Sweden) that was exposed to natural selection and survived Varroa-infestation for more than 10years without treatment. Eventually this population became resistant to the parasite by suppressing the reproduction of the mite. A previous QTL mapping study had identified a region on chromosome 7 with major loci contributing to the mite resistance. Here, a microsatellite scan of the significant candidate QTL regions was used to investigate potential footprints of selection in the original population by comparing the study population on Gotland before (2000) and after selection (2007). Genetic drift had caused an extreme loss of genetic diversity in the 2007 population for all genetic markers tested. In addition to this overall reduction of heterozygosity, two loci on chromosome 7 showed an even stronger and significant reduction in diversity than expected from genetic drift alone. Within the selective sweep eleven genes are annotated, one of them being a putative candidate to interfere with reduced mite reproduction. A glucose-methanol-choline oxidoreductase (GMCOX18) might be involved in changing volatiles emitted by bee larvae that might be essential to trigger oogenesis in Varroa.

Highly efficient integration and expression of piggyBac-derived cassettes in the honeybee (Apis mellifera)

Honeybees (Apis mellifera), which are important pollinators of plants, display remarkable individual behaviors that collectively contribute to the organization of a complex society. Advances in dissecting the complex processes of honeybee behavior have been limited in the recent past due to a lack of genetic manipulation tools. These tools are difficult to apply in honeybees because the unit of reproduction is the colony, and many interesting phenotypes are developmentally specified at later stages. Here, we report highly efficient integration and expression of piggyBac-derived cassettes in the honeybee. We demonstrate that 27 and 20% of queens stably transmitted two different expression cassettes to their offspring, which is a 6- to 30-fold increase in efficiency compared with those generally reported in other insect species. This high efficiency implies that an average beekeeping facility with a limited number of colonies can apply this tool. We demonstrated that the cassette stably and efficiently expressed marker genes in progeny under either an artificial or an endogenous promoter. This evidence of efficient expression encourages the use of this system to inhibit gene functions through RNAi in specific tissues and developmental stages by using various promoters. We also showed that the transgenic marker could be used to select transgenic offspring to be employed to facilitate the building of transgenic colonies via the haploid males. We present here the first to our knowledge genetic engineering tool that will efficiently allow for the systematic detection and better understanding of processes underlying the biology of honeybees.

High-resolution linkage map for two honeybee chromosomes: the hotspot quest

Meiotic recombination is a fundamental process ensuring proper disjunction of homologous chromosomes and allele shuffling in successive generations. In many species, this cellular mechanism occurs heterogeneously along chromosomes and mostly concentrates in tiny fragments called recombination hotspots. Specific DNA motifs have been shown to initiate recombination in these hotspots in mammals, fission yeast and drosophila. The aim of this study was to check whether recombination also occurs in a heterogeneous fashion in the highly recombinogenic honeybee genome and whether this heterogeneity can be connected with specific DNA motifs. We completed a previous picture drawn from a routine genetic map built with an average resolution of 93 kb. We focused on the two smallest honeybee chromosomes to increase the resolution and even zoomed at very high resolution (3.6 kb) on a fragment of 300 kb. Recombination rates measured in these fragments were placed in relation with occurrence of 30 previously described motifs through a Poisson regression model. A selection procedure suitable for correlated variables was applied to keep significant motifs. These fine and ultra-fine mappings show that recombination rate is significantly heterogeneous although poorly contrasted between high and low recombination rate, contrarily to most model species. We show that recombination rate is probably associated with the DNA methylation state. Moreover, three motifs (CGCA, GCCGC and CCAAT) are good candidates of signals promoting recombination. Their influence is however moderate, doubling at most the recombination rate. This discovery extends the way to recombination dissection in insects.

A SNP based high-density linkage map of Apis cerana reveals a high recombination rate similar to Apis mellifera

Background

The Eastern honey bee, Apis cerana Fabricius, is distributed in southern and eastern Asia, from India and China to Korea and Japan and southeast to the Moluccas. This species is also widely kept for honey production besides Apis mellifera. Apis cerana is also a model organism for studying social behavior, caste determination, mating biology, sexual selection, and host-parasite interactions. Few resources are available for molecular research in this species, and a linkage map was never constructed. A linkage map is a prerequisite for quantitative trait loci mapping and for analyzing genome structure. We used the Chinese honey bee, Apis cerana cerana to construct the first linkage map in the Eastern honey bee.

Results

F2 workers (N = 103) were genotyped for 126,990 single nucleotide polymorphisms (SNPs). After filtering low quality and those not passing the Mendel test, we obtained 3,000 SNPs, 1,535 of these were informative and used to construct a linkage map. The preliminary map contains 19 linkage groups, we then mapped the 19 linkage groups to 16 chromosomes by comparing the markers to the genome of A. mellfiera. The final map contains 16 linkage groups with a total of 1,535 markers. The total genetic distance is 3,942.7 centimorgans (cM) with the largest linkage group (180 loci) measuring 574.5 cM. Average marker interval for all markers across the 16 linkage groups is 2.6 cM.

Conclusion

We constructed a high density linkage map for A. c. cerana with 1,535 markers. Because the map is based on SNP markers, it will enable easier and faster genotyping assays than randomly amplified polymorphic DNA or microsatellite based maps used in A. mellifera.

RESTseq--efficient benchtop population genomics with RESTriction Fragment SEQuencing

We present RESTseq, an improved approach for a cost efficient, highly flexible and repeatable enrichment of DNA fragments from digested genomic DNA using Next Generation Sequencing platforms including small scale Personal Genome sequencers. Easy adjustments make it suitable for a wide range of studies requiring SNP detection or SNP genotyping from fine-scale linkage mapping to population genomics and population genetics also in non-model organisms. We demonstrate the validity of our approach by comparing two honeybee and several stingless bee samples.

High-resolution linkage analyses to identify genes that influence Varroa sensitive hygiene behavior in honey bees

Varroa mites (V. destructor) are a major threat to honey bees (Apis melilfera) and beekeeping worldwide and likely lead to colony decline if colonies are not treated. Most treatments involve chemical control of the mites; however, Varroa has evolved resistance to many of these miticides, leaving beekeepers with a limited number of alternatives. A non-chemical control method is highly desirable for numerous reasons including lack of chemical residues and decreased likelihood of resistance. Varroa sensitive hygiene behavior is one of two behaviors identified that are most important for controlling the growth of Varroa populations in bee hives. To identify genes influencing this trait, a study was conducted to map quantitative trait loci (QTL). Individual workers of a backcross family were observed and evaluated for their VSH behavior in a mite-infested observation hive. Bees that uncapped or removed pupae were identified. The genotypes for 1,340 informative single nucleotide polymorphisms were used to construct a high-resolution genetic map and interval mapping was used to analyze the association of the genotypes with the performance of Varroa sensitive hygiene. We identified one major QTL on chromosome 9 (LOD score = 3.21) and a suggestive QTL on chromosome 1 (LOD = 1.95). The QTL confidence interval on chromosome 9 contains the gene ‘no receptor potential A’ and a dopamine receptor. ‘No receptor potential A’ is involved in vision and olfaction in Drosophila, and dopamine signaling has been previously shown to be required for aversive olfactory learning in honey bees, which is probably necessary for identifying mites within brood cells. Further studies on these candidate genes may allow for breeding bees with this trait using marker-assisted selection.

Fine-scale linkage mapping reveals a small set of candidate genes influencing honey bee grooming behavior in response to Varroa mites

Populations of honey bees in North America have been experiencing high annual colony mortality for 15-20 years. Many apicultural researchers believe that introduced parasites called Varroa mites (V. destructor) are the most important factor in colony deaths. One important resistance mechanism that limits mite population growth in colonies is the ability of some lines of honey bees to groom mites from their bodies. To search for genes influencing this trait, we used an Illumina Bead Station genotyping array to determine the genotypes of several hundred worker bees at over a thousand single-nucleotide polymorphisms in a family that was apparently segregating for alleles influencing this behavior. Linkage analyses provided a genetic map with 1,313 markers anchored to genome sequence. Genotypes were analyzed for association with grooming behavior, measured as the time that individual bees took to initiate grooming after mites were placed on their thoraces. Quantitative-trait-locus interval mapping identified a single chromosomal region that was significant at the chromosome-wide level (p<0.05) on chromosome 5 with a LOD score of 2.72. The 95% confidence interval for quantitative trait locus location contained only 27 genes (honey bee official gene annotation set 2) including Atlastin, Ataxin and Neurexin-1 (AmNrx1), which have potential neurodevelopmental and behavioral effects. Atlastin and Ataxin homologs are associated with neurological diseases in humans. AmNrx1 codes for a presynaptic protein with many alternatively spliced isoforms. Neurexin-1 influences the growth, maintenance and maturation of synapses in the brain, as well as the type of receptors most prominent within synapses. Neurexin-1 has also been associated with autism spectrum disorder and schizophrenia in humans, and self-grooming behavior in mice.

Multiple mating but not recombination causes quantitative increase in offspring genetic diversity for varying genetic architectures

Explaining the evolution of sex and recombination is particularly intriguing for some species of eusocial insects because they display exceptionally high mating frequencies and genomic recombination rates. Explanations for both phenomena are based on the notion that both increase colony genetic diversity, with demonstrated benefits for colony disease resistance and division of labor. However, the relative contributions of mating number and recombination rate to colony genetic diversity have never been simultaneously assessed. Our study simulates colonies, assuming different mating numbers, recombination rates, and genetic architectures, to assess their worker genotypic diversity. The number of loci has a strong negative effect on genotypic diversity when the allelic effects are inversely scaled to locus number. In contrast, dominance, epistasis, lethal effects, or limiting the allelic diversity at each locus does not significantly affect the model outcomes. Mating number increases colony genotypic variance and lowers variation among colonies with quickly diminishing returns. Genomic recombination rate does not affect intra- and inter-colonial genotypic variance, regardless of mating frequency and genetic architecture. Recombination slightly increases the genotypic range of colonies and more strongly the number of workers with unique allele combinations across all loci. Overall, our study contradicts the argument that the exceptionally high recombination rates cause a quantitative increase in offspring genotypic diversity across one generation. Alternative explanations for the evolution of high recombination rates in social insects are therefore needed. Short-term benefits are central to most explanations of the evolution of multiple mating and high recombination rates in social insects but our results also apply to other species.

Genetic variability of arrhenotokous and thelytokous Venturia canescens (Hymenoptera)

The ichneumonid wasp Venturia canescens (Hymenoptera) has been studied extensively for foraging behaviour and population dynamics of sexually (arrhenotokous) and parthenogenetically (thelytokous) reproducing individuals. Here we report the development of a set of microsatellite markers for V.canescens and use them to show that arrhenotokous individuals have more genetic variability than thelytokous ones, which are even homozygous for all tested loci. Crosses between arrhenotokous individuals suggested one marker, Vcan071, to be linked with the Complementary Sex Determiner (CSD) locus and one, Vcan109, with the Virus Like Protein (vlp-p40) locus. The genome size of V. canescens was estimated to be 274–279 Mb. We discuss how both reproductive modes can give rise to the observed genetic variability and how the new markers can be used for future genetic studies of V. canescens.

Pleiotropy of segregating genetic variants that affect honey bee worker life expectancy

In contrast to many other complex traits, the natural genetic architecture of life expectancy has not been intensely studied, particularly in non-model organisms, such as the honey bee (Apis mellifera L.). Multiple factors that determine honey bee worker lifespan have been identified and genetic analyses have been performed on some of those traits. Several of the traits are included in a suite of correlated traits that form the pollen hoarding syndrome, which was named after the behavior to store surplus pollen in the nest and is tied to social evolution. Here, seven quantitative trait loci that had previously been identified for their effects on different aspects of the pollen hoarding syndrome were studied for their genetic influence on the survival of adult honey bee workers. To gain a more comprehensive understanding of the genetic architecture of worker longevity, a panel of 280 additional SNP markers distributed across the genome was also tested. Allelic distributions were compared between young and old bees in two backcross populations of the bi-directionally selected high- and low-pollen hoarding strain. Our results suggest a pleiotropic effect of at least one of the behavioral quantitative trait loci on worker longevity and one significant and several other putative genetic effects in other genomic regions. At least one locus showed evidence for strong antagonistic pleiotropy and several others suggested genetic factors that influence pre-emergence survival of worker honey bees. Thus, the predicted association between worker lifespan and the pollen hoarding syndrome was supported at the genetic level and the magnitude of the identified effects also strengthened the view that naturally segregating genetic variation can have major effects on age-specific survival probability in the wild.

Fine scale analysis of crossover and non-crossover and detection of recombination sequence motifs in the honeybee (Apis mellifera)

BACKGROUND

Meiotic exchanges are non-uniformly distributed across the genome of most studied organisms. This uneven distribution suggests that recombination is initiated by specific signals and/or regulations. Some of these signals were recently identified in humans and mice. However, it is unclear whether or not sequence signals are also involved in chromosomal recombination of insects.

METHODOLOGY

We analyzed recombination frequencies in the honeybee, in which genome sequencing provided a large amount of SNPs spread over the entire set of chromosomes. As the genome sequences were obtained from a pool of haploid males, which were the progeny of a single queen, an oocyte method (study of recombination on haploid males that develop from unfertilized eggs and hence are the direct reflect of female gametes haplotypes) was developed to detect recombined pairs of SNP sites. Sequences were further compared between recombinant and non-recombinant fragments to detect recombination-specific motifs.

CONCLUSIONS

Recombination events between adjacent SNP sites were detected at an average distance of 92 bp and revealed the existence of high rates of recombination events. This study also shows the presence of conversion without crossover (i. e. non-crossover) events, the number of which largely outnumbers that of crossover events. Furthermore the comparison of sequences that have undergone recombination with sequences that have not, led to the discovery of sequence motifs (CGCA, GCCGC, CCGCA), which may correspond to recombination signals.

Efficacy of random primer-pair arrays in plant genome analysis: a case study of Cucumis (Cucurbitaceae) for identification of wild and cultivated species

The efficacy of random primer-pair arrays compared to conventional RAPD method with a single decamer primer was evaluated using DNA from two species of Cucumis. The banding patterns of amplicons revealed enhanced utility of primer-pair arrays over conventional RAPDs, producing more bands and a higher degree of polymorphism, both at intra- and inter-specific levels. Amplification produced by both methods clearly distinguished a wild from a cultivated species of the genus Cucumis. The main advantage of the primer-pair RAPD over single-primer-based RAPD is the increase in the number of reactions and amplification products in the form of novel/unique bands with a limited number of primers. It also enables the generation of reliable amplicons with a large number of polymorphic bands, which can be linked to gene-governing traits, allowing sequence-characterized partial genome analysis.

Regulation of behaviorally associated gene networks in worker honey bee ovaries

Several lines of evidence support genetic links between ovary size and division of labor in worker honey bees. However, it is largely unknown how ovaries influence behavior. To address this question, we first performed transcriptional profiling on worker ovaries from two genotypes that differ in social behavior and ovary size. Then, we contrasted the differentially expressed ovarian genes with six sets of available brain transcriptomes. Finally, we probed behavior-related candidate gene networks in wild-type ovaries of different sizes. We found differential expression in 2151 ovarian transcripts in these artificially selected honey bee strains, corresponding to approximately 20.3% of the predicted gene set of honey bees. Differences in gene expression overlapped significantly with changes in the brain transcriptomes. Differentially expressed genes were associated with neural signal transmission (tyramine receptor, TYR) and ecdysteroid signaling; two independently tested nuclear hormone receptors (HR46 and ftz-f1) were also significantly correlated with ovary size in wild-type bees. We suggest that the correspondence between ovary and brain transcriptomes identified here indicates systemic regulatory networks among hormones (juvenile hormone and ecdysteroids), pheromones (queen mandibular pheromone), reproductive organs and nervous tissues in worker honey bees. Furthermore, robust correlations between ovary size and neuraland endocrine response genes are consistent with the hypothesized roles of the ovaries in honey bee behavioral regulation.

Recombination rate variation in closely related species

Despite their importance to successful meiosis and various evolutionary processes, meiotic recombination rates sometimes vary within species or between closely related species. For example, humans and chimpanzees share virtually no recombination hotspot locations in the surveyed portion of the genomes. However, conservation of recombination rates between closely related species has also been documented, raising an apparent contradiction. Here, we evaluate how and why conflicting patterns of recombination rate conservation and divergence may be observed, with particular emphasis on features that affect recombination, and the scale and method with which recombination is surveyed. Additionally, we review recent studies identifying features influencing fine-scale and broad-scale recombination patterns and informing how quickly recombination rates evolve, how changes in recombination impact selection and evolution in natural populations, and more broadly, which forces influence genome evolution.

A high recombination rate in eusocial Hymenoptera: evidence from the common wasp Vespula vulgaris

Background

High recombination rates have previously been detected in two groups of eusocial insects; honeybees and ants. In this study we estimate recombination rate in a eusocial wasp Vespula vulgaris that represents a third phylogenetic lineage within eusocial hymenopterans.

Results

A genetic linkage map of V. vulgaris based on 210 markers shows that the total map length is 2129 cM and the recombination rate is 9.7 cM/Mb (or 103 kb/cM). The present estimate in V. vulgaris is somewhat smaller than in the honeybee Apis mellifera and intermediate between the estimates from two ant species (Acromyrmex echinatior, Pogonomyrmex rugosus). Altogether, the estimates from these eusocial species are higher than in any other insect reported so far.

Conlusions

The four species (V. vulgaris, A. mellifera, A. echinatior, P. rugosus) are characterized by advanced eusociality with large colonies, clear queen-worker dimorphism and well developed task specialization. They also have colonies with a single, normally multiply inseminated (polyandrous) queen. Benefits of genotypic diversity within colonies (e.g. through improved task specialization or pathogen and parasite resistance) may have selected for both polyandry and high recombination rate in such advanced eusocial insects.

Hymenoptera Genome Database: integrated community resources for insect species of the order Hymenoptera

The Hymenoptera Genome Database (HGD) is a comprehensive model organism database that caters to the needs of scientists working on insect species of the order Hymenoptera. This system implements open-source software and relational databases providing access to curated data contributed by an extensive, active research community. HGD contains data from 9 different species across ∼200 million years in the phylogeny of Hymenoptera, allowing researchers to leverage genetic, genome sequence and gene expression data, as well as the biological knowledge of related model organisms. The availability of resources across an order greatly facilitates comparative genomics and enhances our understanding of the biology of agriculturally important Hymenoptera species through genomics. Curated data at HGD includes predicted and annotated gene sets supported with evidence tracks such as ESTs/cDNAs, small RNA sequences and GC composition domains. Data at HGD can be queried using genome browsers and/or BLAST/PSI-BLAST servers, and it may also be downloaded to perform local searches. We encourage the public to access and contribute data to HGD at: http://HymenopteraGenome.org.

Comparative linkage mapping suggests a high recombination rate in all honeybees

Meiotic recombination is required for proper chromosome assortment, and accordingly, 1-2 chiasmata per chromosome are found in most species. However, observed recombination rates deviate in some cases from neutral expectations between and within genomes and may play an important role in adaptive evolution. One potentially important argument for an adaptive evolution of recombination rates is the exceptionally high genome-wide recombination rates of social Hymenoptera, in particular the Western honeybee, Apis mellifera. It has the highest metazoan recombination rate reported so far. Proximate or ultimate causes for this elevated recombination rate have not yet been resolved. In a comparative study, we investigated meiotic recombination in the red dwarf honeybee Apis florea. Microsatellite markers developed for A. mellifera were genotyped in a natural mapping population of A. florea. From these genotypes, we calculated local recombination rates, using the physical distances from A. mellifera. In addition to a few comparisons of intervals across the genome, we particularly focused on chromosomes 3 and 12. Confirming marker synteny, we found that recombination rates in A. florea are as high as or higher than those in A. mellifera. Our results are limited to select genomic regions but suggest that A. florea also exhibits an exceptionally high genome-wide recombination rate. This trait may thus occur genus wide. Although our study cannot identify a single explanation for the high rates of recombination in Apis, it favors hypotheses that apply to the entire genus. Furthermore, we conclude that the genome structure of the 2 species has been largely conserved, at least in the parts we investigated.

Recombination and its impact on the genome of the haplodiploid parasitoid wasp Nasonia

Homologous meiotic recombination occurs in most sexually reproducing organisms, yet its evolutionary advantages are elusive. Previous research explored recombination in the honeybee, a eusocial hymenopteran with an exceptionally high genome-wide recombination rate. A comparable study in a non-social member of the Hymenoptera that would disentangle the impact of sociality from Hymenoptera-specific features such as haplodiploidy on the evolution of the high genome-wide recombination rate in social Hymenoptera is missing. Utilizing single-nucleotide polymorphisms (SNPs) between two Nasonia parasitoid wasp genomes, we developed a SNP genotyping microarray to infer a high-density linkage map for Nasonia. The map comprises 1,255 markers with an average distance of 0.3 cM. The mapped markers enabled us to arrange 265 scaffolds of the Nasonia genome assembly 1.0 on the linkage map, representing 63.6% of the assembled N. vitripennis genome. We estimated a genome-wide recombination rate of 1.4-1.5 cM/Mb for Nasonia, which is less than one tenth of the rate reported for the honeybee. The local recombination rate in Nasonia is positively correlated with the distance to the center of the linkage groups, GC content, and the proportion of simple repeats. In contrast to the honeybee genome, gene density in the parasitoid wasp genome is positively associated with the recombination rate; regions of low recombination are characterized by fewer genes with larger introns and by a greater distance between genes. Finally, we found that genes in regions of the genome with a low recombination frequency tend to have a higher ratio of non-synonymous to synonymous substitutions, likely due to the accumulation of slightly deleterious non-synonymous substitutions. These findings are consistent with the hypothesis that recombination reduces interference between linked sites and thereby facilitates adaptive evolution and the purging of deleterious mutations. Our results imply that the genomes of haplodiploid and of diploid higher eukaryotes do not differ systematically in their recombination rates and associated parameters.

Sex determination in honeybees: two separate mechanisms induce and maintain the female pathway

Sex determination in honeybees is realized by the csd and the fem gene that establish and maintain, throughout development, sexual fates via the control of alternative splicing.

Organisms have evolved a bewildering diversity of mechanisms to generate the two sexes. The honeybee (Apis mellifera) employs an interesting system in which sex is determined by heterozygosity at a single locus (the Sex Determination Locus) harbouring the complementary sex determiner (csd) gene. Bees heterozygous at Sex Determination Locus are females, whereas bees homozygous or hemizygous are males. Little is known, however, about the regulation that links sex determination to sexual differentiation. To investigate the control of sexual development in honeybees, we analyzed the functions and the regulatory interactions of genes involved in the sex determination pathway. We show that heterozygous csd is only required to induce the female pathway, while the feminizer (fem) gene maintains this decision throughout development. By RNAi induced knockdown we show that the fem gene is essential for entire female development and that the csd gene exclusively processes the heterozygous state. Fem activity is also required to maintain the female determined pathway throughout development, which we show by mosaic structures in fem-repressed intersexuals. We use expression of Fem protein in males to demonstrate that the female maintenance mechanism is controlled by a positive feedback splicing loop in which Fem proteins mediate their own synthesis by directing female fem mRNA splicing. The csd gene is only necessary to induce this positive feedback loop in early embryogenesis by directing splicing of fem mRNAs. Finally, fem also controls the splicing of Am-doublesex transcripts encoding conserved male- and female-specific transcription factors involved in sexual differentiation. Our findings reveal how the sex determination process is realized in honeybees differing from Drosophila melanogaster.

Sexual differentiation is a fundamental process in the animal kingdom, and different species have evolved a bewildering diversity of mechanisms to generate the two sexes in the proper proportions. Sex determination in honeybees (Apis mellifera) provides an interesting and unusual system to study, as it is governed by heterozygosity of a single locus harbouring the complementary sex determiner gene (csd), in contrast to the well-studied sex chromosome system of Drosophila melanogaster. We show that the female sex determination pathway is exclusively induced by the csd gene in early embryogenesis. Later on and throughout development this inductive signal is maintained via a positive feedback loop of the feminizer (fem) gene, in which the Fem protein mediates its own synthesis. The findings reveal how the sex determination process in honeybees is realized by the regulation and function of two genes differing from Drosophila.

The genetic basis of transgressive ovary size in honeybee workers

Ovarioles are the functional unit of the female insect reproductive organs and the number of ovarioles per ovary strongly influences egg-laying rate and fecundity. Social evolution in the honeybee (Apis mellifera) has resulted in queens with 200-360 total ovarioles and workers with usually 20 or less. In addition, variation in ovariole number among workers relates to worker sensory tuning, foraging behavior, and the ability to lay unfertilized male-destined eggs. To study the genetic architecture of worker ovariole number, we performed a series of crosses between Africanized and European bees that differ in worker ovariole number. Unexpectedly, these crosses produced transgressive worker phenotypes with extreme ovariole numbers that were sensitive to the social environment. We used a new selective pooled DNA interval mapping approach with two Africanized backcrosses to identify quantitative trait loci (QTL) underlying the transgressive ovary phenotype. We identified one QTL on chromosome 11 and found some evidence for another QTL on chromosome 2. Both QTL regions contain plausible functional candidate genes. The ovariole number of foragers was correlated with the sugar concentration of collected nectar, supporting previous studies showing a link between worker physiology and foraging behavior. We discuss how the phenotype of extreme worker ovariole numbers and the underlying genetic factors we identified could be linked to the development of queen traits.

Novel microsatellite DNA loci for Bombus terrestris (Linnaeus, 1758)

We present details and characteristics of 123 novel polymorphic microsatellite DNA loci for Bombus terrestris. Thirty-four of these loci have been tested in nine other Bombus species and 25 of them showed polymorphisms in at least one species. These microsatellite DNA loci together with the already established 60 loci will be useful for characterizing wild and managed populations of B. terrestris and other Bombus species as well as for detailed genetic studies in including mapping studies and genome annotations.

An olfactory shift is associated with male perfume differentiation and species divergence in orchid bees

Saltational changes may underlie the diversification of pheromone communication systems in insects, which are normally under stabilizing selection favoring high specificity in signals and signal perception. In orchid bees (Euglossini), the production of male signals depends on the sense of smell: males collect complex blends of volatiles (perfumes) from their environment, which are later emitted as pheromone analogs at mating sites. We analyzed the behavioral and antennal response to perfume components in two male morphotypes of Euglossa cf. viridissima from Mexico, which differ in the number of mandibular teeth. Tridentate males collected 2-hydroxy-6-nona-1,3-dienyl-benzaldehyde (HNDB) as the dominant component of their perfume. In bidentate males, blends were broadly similar but lacked HNDB. Population genetic analysis revealed that tri- and bidentate males belong to two reproductively isolated lineages. Electroantennogram tests (EAG and GC-EAD) showed substantially lower antennal responses to HNDB in bidentate versus tridentate males, revealing for the first time a mechanism by which closely related species acquire different chemical compounds from their habitat. The component-specific differences in perfume perception and collection in males of two sibling species are in agreement with a saltational, olfaction-driven mode of signal perfume evolution. However, the response of females to the diverged signals remains unknown.

Identification of Quantitative Trait Loci and candidate genes influencing ethanol sensitivity in honey bees

Invertebrate models have greatly furthered our understanding of ethanol sensitivity and alcohol addiction. The honey bee (Apis mellifera), a widely used behavioral model, is valuable for comparative studies. A quantitative trait locus (QTL) mapping experiment was designed to identify QTL and genes influencing ethanol vapor sensitivity. A backcross mating between ethanol-sensitive and resistant lines resulted in worker offspring that were tested for sensitivity to the sedative effects of alcohol. A linkage map was constructed with over 500 amplified fragment length polymorphism (AFLP) and sequence-tagged site (STS) markers. Four QTL were identified from three linkage groups with log of odds ratio (LOD) scores of 2.28, 2.26, 2.23, and 2.02. DNA from markers within and near QTL were cloned and sequenced, and this data was utilized to integrate our map with the physical honey bee genome. Many candidate genes were identified that influence synaptic transmission, neuronal growth, and detoxification. Others affect lipid synthesis, apoptosis, alcohol metabolism, cAMP signaling, and electron transport. These results are relevant because they present the first search for QTL that affect resistance to acute ethanol exposure in an invertebrate, could be useful for comparative genomic purposes, and lend credence to the use of honey bees as biomedical models of alcohol metabolism and sensitivity.

Sex determination in the hymenoptera

The dominant and ancestral mode of sex determination in the Hymenoptera is arrhenotokous parthenogenesis, in which diploid females develop from fertilized eggs and haploid males develop from unfertilized eggs. We discuss recent progress in the understanding of the genetic and cytoplasmic mechanisms that make arrhenotoky possible. The best-understood mode of sex determination in the Hymenoptera is complementary sex determination (CSD), in which diploid males are produced under conditions of inbreeding. The gene mediating CSD has recently been cloned in the honey bee and has been named the complementary sex determiner. However, CSD is only known from 4 of 21 hymenopteran superfamilies, with some taxa showing clear evidence of the absence of CSD. Sex determination in the model hymenopteran Nasonia vitripennis does not involve CSD, but it is consistent with a form of genomic imprinting in which activation of the female developmental pathway requires paternally derived genes. Some other hymenopterans are not arrhenotokous but instead exhibit thelytoky or paternal genome elimination.

A third-generation microsatellite-based linkage map of the honey bee, Apis mellifera, and its comparison with the sequence-based physical map

The meiotic map of the honey bee is presented, including the main features that emerged from comparisons with the sequence-based physical map. The map is based on 2,008 markers and is about 40 M long, corresponding to a recombination rate of 22 cM/Mb.

Background:

The honey bee is a key model for social behavior and this feature led to the selection of the species for genome sequencing. A genetic map is a necessary companion to the sequence. In addition, because there was originally no physical map for the honey bee genome project, a meiotic map was the only resource for organizing the sequence assembly on the chromosomes.

Results:

We present the genetic (meiotic) map here and describe the main features that emerged from comparison with the sequence-based physical map. The genetic map of the honey bee is saturated and the chromosomes are oriented from the centromeric to the telomeric regions. The map is based on 2,008 markers and is about 40 Morgans (M) long, resulting in a marker density of one every 2.05 centiMorgans (cM). For the 186 megabases (Mb) of the genome mapped and assembled, this corresponds to a very high average recombination rate of 22.04 cM/Mb. Honey bee meiosis shows a relatively homogeneous recombination rate along and across chromosomes, as well as within and between individuals. Interference is higher than inferred from the Kosambi function of distance. In addition, numerous recombination hotspots are dispersed over the genome.

Conclusion:

The very large genetic length of the honey bee genome, its small physical size and an almost complete genome sequence with a relatively low number of genes suggest a very promising future for association mapping in the honey bee, particularly as the existence of haploid males allows easy bulk segregant analysis.