The complete mitochondrial genome of the cavity-nesting honeybee, Apis cerana (Insecta: Hymenoptera: Apidae) from Borneo.

The complete mitochondrial genome of the cavity-nesting honeybee Apis cerana from Sabah on Borneo Island was analyzed using next-generation sequencing. The mitochondrial genome of A. cerana was a circular molecule of 15,884 bp and was similar to that of the other cavity-nesting honeybee species. The average AT content in the A. cerana mitochondrial genome was 84.4%. It was predicted to contain 13 protein-coding, 22 tRNA, and two rRNA genes, along with one A + T-rich control region.

The complete mitochondrial DNA sequence of endemic honeybee Apis nuluensis (Insecta: Hymenoptera: Apidae) inhabiting Mount Kinabalu in Sabah Province, Borneo Island.

The cavity-nesting honeybee Apis nuluensis inhabits only the highlands of Mount Kinabalu of Sabah, Borneo Island. The mitochondrial genome is a circular molecule of approximately 1.6 kb that includes 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one AT-rich control region. The average AT content was 84.5%. The start codons ATC, ATG, and ATT were found in one, three, and nine genes, respectively, whereas the stop codon TAA was observed in all genes. The phylogenetic relationship, inferred using 13 PCGs, was consistent with that reported in previous studies that predicted a sister taxon relationship between A. nuluensis and A. cerana.

The complete mitochondrial genome of the cavity-nesting honeybee, Apis koschevnikovi (Insecta: Hymenoptera: Apidae).

In this study, we analyzed the complete mitochondrial genome of the cavity-nesting honeybee, A. koschevnikovi. The mitochondrial genome of A. koschevnikovi was observed to be a circular molecule of 15,278 bp and was similar to that of the other cavity-nesting honeybee species. The average AT content in the A. koschevnikovi mitochondrial genome was 84%. It was predicted to contain 13 protein-coding, 24 tRNA and two rRNA genes, along with one A + T-rich control region, besides three tRNA-Met repeats.

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.

Cross-species correlation between queen mating numbers and worker ovary sizes suggests kin conflict may influence ovary size evolution in honeybees.

During social evolution, the ovary size of reproductively specialized honey bee queens has dramatically increased while their workers have evolved much smaller ovaries. However, worker division of labor and reproductive competition under queenless conditions are influenced by worker ovary size. Little comparative information on ovary size exists in the different honey bee species. Here, we report ovariole numbers of freshly dissected workers from six Apis species from two locations in Southeast Asia. The average number of worker ovarioles differs significantly among species. It is strongly correlated with the average mating number of queens, irrespective of body size. Apis dorsata, in particular, is characterized by numerous matings and very large worker ovaries. The relation between queen mating number and ovary size across the six species suggests that individual selection via reproductive competition plays a role in worker ovary size evolution. This indicates that genetic diversity, generated by multiple mating, may bear a fitness cost at the colony level.

A simple and distinctive microbiota associated with honey bees and bumble bees.

Specialized relationships with bacteria often allow animals to exploit a new diet by providing a novel set of metabolic capabilities. Bees are a monophyletic group of Hymenoptera that transitioned to a completely herbivorous diet from the carnivorous diet of their wasp ancestors. Recent culture-independent studies suggest that a set of distinctive bacterial species inhabits the gut of the honey bee, Apis mellifera. Here we survey the gut microbiotae of diverse bee and wasp species to test whether acquisition of these bacteria was associated with the transition to herbivory in bees generally. We found that most bee species lack phylotypes that are the same or similar to those typical of A. mellifera, rejecting the hypothesis that this dietary transition was symbiont-dependent. The most common bacteria in solitary bee species are a widespread phylotype of Burkholderia and the pervasive insect associate, Wolbachia. In contrast, several social representatives of corbiculate bees do possess distinctive bacterial phylotypes. Samples of A. mellifera harboured the same microbiota as in previous surveys, and closely related bacterial phylotypes were identified in two Asian honey bees (Apis andreniformis and Apis dorsata) and several bumble bee (Bombus) species. Potentially, the sociality of Apis and Bombus species facilitates symbiont transmission and thus is key to the maintenance of a more consistent gut microbiota. Phylogenetic analyses provide a more refined taxonomic placement of the A. mellifera symbionts.

Rapid evolution of immune proteins in social insects.

The existence of behavioral traits connected to defense against pathogens manifests the importance of pathogens in the evolution of social insects. However, very little is known about how pathogen pressure has affected the molecular evolution of genes involved in their innate immune system. We have studied the sequence evolution of several immune genes in ants and honeybees. The results show high rates of evolution in both ants and honeybees as measured by the ratio of amino acid changes to silent nucleotide changes, the ratio being clearly higher than in Drosophila immune genes or in nonimmunity genes of bees. This conforms to our expectations based on high pathogen pressure in social insects. The codon-based likelihood method found clear evidence of positive selection only in one ant gene, even though positive selection has earlier been found in both ant and termite immune genes. There is now indication that selection on the amino acid composition of the immune-related genes has been an important part in the fight against pathogens by social insects. However, we cannot distinguish in all the cases whether the high observed d(N)/d(S) ratio results from positive selection within a restricted part of the studied genes or from relaxation of purifying selection associated with effective measures of behaviorally based colony-level defenses.

Evidence for convergent nucleotide evolution and high allelic turnover rates at the complementary sex determiner gene of Western and Asian honeybees.

Our understanding of the impact of recombination, mutation, genetic drift, and selection on the evolution of a single gene is still limited. Here we investigate the impact of all these evolutionary forces at the complementary sex determiner (csd) gene that evolves under a balancing mode of selection. Females are heterozygous at the csd gene and males are hemizygous; diploid males are lethal and occur when csd is homozygous. Rare alleles thus have a selective advantage, are seldom lost by the effect of genetic drift, and are maintained over extended periods of time when compared with neutral polymorphisms. Here, we report on the analysis of 17, 19, and 15 csd alleles of Apis cerana, Apis dorsata, and Apis mellifera honeybees, respectively. We observed great heterogeneity of synonymous (piS) and nonsynonymous (piN) polymorphisms across the gene, with a consistent peak in exons 6 and 7. We propose that exons 6 and 7 encode the potential specifying domain (csd-PSD) that has accumulated elevated nucleotide polymorphisms over time by balancing selection. We observed no direct evidence that balancing selection favors the accumulation of nonsynonymous changes at csd-PSD (piN/piS ratios are all <1, ranging from 0.6 to 0.95). We observed an excess of shared nonsynonymous changes, which suggest that strong evolutionary constraints are operating at csd-PSD resulting in the independent accumulation of the same nonsynonymous changes in different alleles across species (convergent evolution). Analysis of csd-PSD genealogy revealed relatively short average coalescence times ( approximately 6 Myr), low average synonymous nucleotide diversity (piS < 0.09), and a lack of trans-specific alleles that substantially contrasts with previously analyzed loci under strong balancing selection. We excluded the possibility of a burst of diversification after population bottlenecking and intragenic recombination as explanatory factors, leaving high turnover rates as the explanation for this observation. By comparing observed allele richness and average coalescence times with a simplified model of csd-coalescence, we found that small long-term population sizes (i.e., N(e) < 10(4)), but not high mutation rates, can explain short maintenance times, implicating a strong historical impact of genetic drift on the molecular evolution of highly social honeybees.