Heterochiasmy and the establishment of gsdf as a novel sex determining gene in Atlantic halibut.

Atlantic Halibut (Hippoglossus hippoglossus) has a X/Y genetic sex determination system, but the sex determining factor is not known. We produced a high-quality genome assembly from a male and identified parts of chromosome 13 as the Y chromosome due to sequence divergence between sexes and segregation of sex genotypes in pedigrees. Linkage analysis revealed that all chromosomes exhibit heterochiasmy, i.e. male-only and female-only meiotic recombination regions (MRR/FRR). We show that FRR/MRR intervals differ in nucleotide diversity and repeat class content and that this is true also for other Pleuronectidae species. We further show that remnants of a Gypsy-like transposable element insertion on chr13 promotes early male specific expression of gonadal somatic cell derived factor (gsdf). Less than 4.5 MYA, this male-determining element evolved on an autosomal FRR segment featuring pre-existing male meiotic recombination barriers, thereby creating a Y chromosome. Our findings indicate that heterochiasmy may facilitate the evolution of genetic sex determination systems relying on linkage of sexually antagonistic loci to a sex-determining factor.

Comparative Population Transcriptomics Provide New Insight into the Evolutionary History and Adaptive Potential of World Ocean Krill.

Genetic variation is instrumental for adaptation to changing environments but it is unclear how it is structured and contributes to adaptation in pelagic species lacking clear barriers to gene flow. Here, we applied comparative genomics to extensive transcriptome datasets from 20 krill species collected across the Atlantic, Indian, Pacific, and Southern Oceans. We compared genetic variation both within and between species to elucidate their evolutionary history and genomic bases of adaptation. We resolved phylogenetic interrelationships and uncovered genomic evidence to elevate the cryptic Euphausia similis var. armata into species. Levels of genetic variation and rates of adaptive protein evolution vary widely. Species endemic to the cold Southern Ocean, such as the Antarctic krill Euphausia superba, showed less genetic variation and lower evolutionary rates than other species. This could suggest a low adaptive potential to rapid climate change. We uncovered hundreds of candidate genes with signatures of adaptive evolution among Antarctic Euphausia but did not observe strong evidence of adaptive convergence with the predominantly Arctic Thysanoessa. We instead identified candidates for cold-adaptation that have also been detected in Antarctic fish, including genes that govern thermal reception such as TrpA1. Our results suggest parallel genetic responses to similar selection pressures across Antarctic taxa and provide new insights into the adaptive potential of important zooplankton already affected by climate change.

Extreme Differences in Recombination Rate between the Genomes of a Solitary and a Social Bee.

Social insect genomes exhibit the highest rates of crossing over observed in plants and animals. The evolutionary causes of these extreme rates are unknown. Insight can be gained by comparing recombination rate variation across the genomes of related social and solitary insects. Here, we compare the genomic recombination landscape of the highly social honey bee, Apis mellifera, with the solitary alfalfa leafcutter bee, Megachile rotundata, by analyzing patterns of linkage disequilibrium in population-scale genome sequencing data. We infer that average recombination rates are extremely elevated in A. mellifera compared with M. rotundata. However, our results indicate that similar factors control the distribution of crossovers in the genomes of both species. Recombination rate is significantly reduced in coding regions in both species, with genes inferred to be germline methylated having particularly low rates. Genes with worker-biased patterns of expression in A. mellifera and their orthologs in M. rotundata have higher than average recombination rates in both species, suggesting that selection for higher diversity in genes involved in worker caste functions in social taxa is not the explanation for these elevated rates. Furthermore, we find no evidence that recombination has modulated the efficacy of selection among genes during bee evolution, which does not support the hypothesis that high recombination rates facilitated positive selection for new functions in social insects. Our results indicate that the evolution of sociality in insects likely entailed selection on modifiers that increased recombination rates genome wide, but that the genomic recombination landscape is determined by the same factors.

Two extended haplotype blocks are associated with adaptation to high altitude habitats in East African honey bees.

Understanding the genetic basis of adaption is a central task in biology. Populations of the honey bee Apis mellifera that inhabit the mountain forests of East Africa differ in behavior and morphology from those inhabiting the surrounding lowland savannahs, which likely reflects adaptation to these habitats. We performed whole genome sequencing on 39 samples of highland and lowland bees from two pairs of populations to determine their evolutionary affinities and identify the genetic basis of these putative adaptations. We find that in general, levels of genetic differentiation between highland and lowland populations are very low, consistent with them being a single panmictic population. However, we identify two loci on chromosomes 7 and 9, each several hundred kilobases in length, which exhibit near fixation for different haplotypes between highland and lowland populations. The highland haplotypes at these loci are extremely rare in samples from the rest of the world. Patterns of segregation of genetic variants suggest that recombination between haplotypes at each locus is suppressed, indicating that they comprise independent structural variants. The haplotype on chromosome 7 harbors nearly all octopamine receptor genes in the honey bee genome. These have a role in learning and foraging behavior in honey bees and are strong candidates for adaptation to highland habitats. Molecular analysis of a putative breakpoint indicates that it may disrupt the coding sequence of one of these genes. Divergence between the highland and lowland haplotypes at both loci is extremely high suggesting that they are ancient balanced polymorphisms that greatly predate divergence between the extant honey bee subspecies.

A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds.

BACKGROUND: The ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map. RESULTS: Each of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features. CONCLUSIONS: The improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Chromosomal inversions associated with environmental adaptation in honeybees.

Chromosomal inversions can facilitate local adaptation in the presence of gene flow by suppressing recombination between well-adapted native haplotypes and poorly adapted migrant haplotypes. East African mountain populations of the honeybee Apis mellifera are highly divergent from neighbouring lowland populations at two extended regions in the genome, despite high similarity in the rest of the genome, suggesting that these genomic regions harbour inversions governing local adaptation. Here, we utilize a new highly contiguous assembly of the honeybee genome to characterize these regions. Using whole-genome sequencing data from 55 highland and lowland bees, we find that the highland haplotypes at both regions are present at high frequencies in three independent highland populations but extremely rare elsewhere. The boundaries of both divergent regions are characterized by regions of high homology with each other positioned in opposite orientations and contain highly repetitive, long inverted repeats with homology to transposable elements. These regions are likely to represent inversion breakpoints that participate in nonallelic homologous recombination. Using long-read data, we confirm that the lowland samples are contiguous across breakpoint regions. We do not find evidence for disruption of functional sequence by these breakpoints, which suggests that the inversions are likely maintained due to their allelic content conferring local adaptation in highland environments. Finally, we identify a third divergent genomic region, which contains highly divergent segregating haplotypes that also may contain inversion variants under selection. The results add to a growing body of evidence indicating the importance of chromosomal inversions in local adaptation.

The genomic basis of adaptation to high-altitude habitats in the eastern honey bee (Apis cerana).

The eastern honey bee (Apis cerana) is of central importance for agriculture in Asia. It has adapted to a wide variety of environmental conditions across its native range in southern and eastern Asia, which includes high-altitude regions. eastern honey bees inhabiting mountains differ morphologically from neighbouring lowland populations and may also exhibit differences in physiology and behaviour. We compared the genomes of 60 eastern honey bees collected from high and low altitudes in Yunnan and Gansu provinces, China, to infer their evolutionary history and to identify candidate genes that may underlie adaptation to high altitude. Using a combination of FST -based statistics, long-range haplotype tests and population branch statistics, we identified several regions of the genome that appear to have been under positive selection. These candidate regions were strongly enriched for coding sequences and had high haplotype homozygosity and increased divergence specifically in highland bee populations, suggesting they have been subjected to recent selection in high-altitude habitats. Candidate loci in these genomic regions included genes related to reproduction and feeding behaviour in honey bees. Functional investigation of these candidate loci is necessary to fully understand the mechanisms of adaptation to high-altitude habitats in the eastern honey bee.

Identification of Multiple Loci Associated with Social Parasitism in Honeybees.

In colonies of the honeybee Apis mellifera, the queen is usually the only reproductive female, which produces new females (queens and workers) by laying fertilized eggs. However, in one subspecies of A. mellifera, known as the Cape bee (A. m. capensis), worker bees reproduce asexually by thelytoky, an abnormal form of meiosis where two daughter nucleii fuse to form single diploid eggs, which develop into females without being fertilized. The Cape bee also exhibits a suite of phenotypes that facilitate social parasitism whereby workers lay such eggs in foreign colonies so their offspring can exploit their resources. The genetic basis of this switch to social parasitism in the Cape bee is unknown. To address this, we compared genome variation in a sample of Cape bees with other African populations. We find genetic divergence between these populations to be very low on average but identify several regions of the genome with extreme differentiation. The regions are strongly enriched for signals of selection in Cape bees, indicating that increased levels of positive selection have produced the unique set of derived phenotypic traits in this subspecies. Genetic variation within these regions allows unambiguous genetic identification of Cape bees and likely underlies the genetic basis of social parasitism. The candidate loci include genes involved in ecdysteroid signaling and juvenile hormone and dopamine biosynthesis, which may regulate worker ovary activation and others whose products localize at the centrosome and are implicated in chromosomal segregation during meiosis. Functional analysis of these loci will yield insights into the processes of reproduction and chemical signaling in both parasitic and non-parasitic populations and advance understanding of the process of normal and atypical meiosis.

Extreme recombination frequencies shape genome variation and evolution in the honeybee, Apis mellifera.

Meiotic recombination is a fundamental cellular process, with important consequences for evolution and genome integrity. However, we know little about how recombination rates vary across the genomes of most species and the molecular and evolutionary determinants of this variation. The honeybee, Apis mellifera, has extremely high rates of meiotic recombination, although the evolutionary causes and consequences of this are unclear. Here we use patterns of linkage disequilibrium in whole genome resequencing data from 30 diploid honeybees to construct a fine-scale map of rates of crossing over in the genome. We find that, in contrast to vertebrate genomes, the recombination landscape is not strongly punctate. Crossover rates strongly correlate with levels of genetic variation, but not divergence, which indicates a pervasive impact of selection on the genome. Germ-line methylated genes have reduced crossover rate, which could indicate a role of methylation in suppressing recombination. Controlling for the effects of methylation, we do not infer a strong association between gene expression patterns and recombination. The site frequency spectrum is strongly skewed from neutral expectations in honeybees: rare variants are dominated by AT-biased mutations, whereas GC-biased mutations are found at higher frequencies, indicative of a major influence of GC-biased gene conversion (gBGC), which we infer to generate an allele fixation bias 5 - 50 times the genomic average estimated in humans. We uncover further evidence that this repair bias specifically affects transitions and favours fixation of CpG sites. Recombination, via gBGC, therefore appears to have profound consequences on genome evolution in honeybees and interferes with the process of natural selection. These findings have important implications for our understanding of the forces driving molecular evolution.

Genomewide analysis of admixture and adaptation in the Africanized honeybee.

Genetic exchange by hybridization or admixture can make an important contribution to evolution, and introgression of favourable alleles can facilitate adaptation to new environments. A small number of honeybees (Apis mellifera) with African ancestry were introduced to Brazil ~60 years ago, which dispersed and hybridized with existing managed populations of European origin, quickly spreading across much of the Americas in an example of a massive biological invasion. Here, we analyse whole-genome sequences of 32 Africanized honeybees sampled from throughout Brazil to study the effect of this process on genome diversity. By comparison with ancestral populations from Europe and Africa, we infer that these samples have 84% African ancestry, with the remainder from western European populations. However, this proportion varies across the genome and we identify signals of positive selection in regions with high European ancestry proportions. These observations are largely driven by one large gene-rich 1.4-Mbp segment on chromosome 11 where European haplotypes are present at a significantly elevated frequency and likely confer an adaptive advantage in the Africanized honeybee population. This region has previously been implicated in reproductive traits and foraging behaviour in worker bees. Finally, by analysing the distribution of ancestry tract lengths in the context of the known time of the admixture event, we are able to infer an average generation time of 2.0 years. Our analysis highlights the processes by which populations of mixed genetic ancestry form and adapt to new environments.

Acoelomorph flatworms are deuterostomes related to Xenoturbella.

Xenoturbellida and Acoelomorpha are marine worms with contentious ancestry. Both were originally associated with the flatworms (Platyhelminthes), but molecular data have revised their phylogenetic positions, generally linking Xenoturbellida to the deuterostomes and positioning the Acoelomorpha as the most basally branching bilaterian group(s). Recent phylogenomic data suggested that Xenoturbellida and Acoelomorpha are sister taxa and together constitute an early branch of Bilateria. Here we assemble three independent data sets-mitochondrial genes, a phylogenomic data set of 38,330 amino-acid positions and new microRNA (miRNA) complements-and show that the position of Acoelomorpha is strongly affected by a long-branch attraction (LBA) artefact. When we minimize LBA we find consistent support for a position of both acoelomorphs and Xenoturbella within the deuterostomes. The most likely phylogeny links Xenoturbella and Acoelomorpha in a clade we call Xenacoelomorpha. The Xenacoelomorpha is the sister group of the Ambulacraria (hemichordates and echinoderms). We show that analyses of miRNA complements have been affected by character loss in the acoels and that both groups possess one miRNA and the gene Rsb66 otherwise specific to deuterostomes. In addition, Xenoturbella shares one miRNA with the ambulacrarians, and two with the acoels. This phylogeny makes sense of the shared characteristics of Xenoturbellida and Acoelomorpha, such as ciliary ultrastructure and diffuse nervous system, and implies the loss of various deuterostome characters in the Xenacoelomorpha including coelomic cavities, through gut and gill slits.

Limited parallelism in genetic adaptation to brackish water bodies in European sprat and Atlantic herring.

The European sprat is a small plankton-feeding clupeid present in the northeastern Atlantic Ocean, the Mediterranean Sea as well as in the brackish Baltic Sea and Black Sea. This species is the target of a major fishery and therefore an accurate characterization of its genetic population structure is crucial to delineate proper stock assessments that aid ensuring the fishery's sustainability. Here we present (i) a draft genome assembly, (ii) pooled whole genome sequencing of 19 population samples covering most of the species' distribution range, and (iii) the design and test of a SNP-chip resource and use this to validate the population structure inferred from pooled sequencing. These approaches revealed, using the populations sampled here, three major groups of European sprat: Oceanic, Coastal, and Brackish with limited differentiation within groups even over wide geographical stretches. Genetic structure is largely driven by six large putative inversions that differentiate Oceanic and Brackish sprats, while Coastal populations display intermediate frequencies of haplotypes at each locus. Interestingly, populations from the Baltic and the Black Seas share similar frequencies of haplotypes at these putative inversions despite their distant geographic location. The closely related clupeids European sprat and Atlantic herring both show genetic adaptation to the brackish Baltic Sea, providing an opportunity to explore the extent of genetic parallelism. This analysis revealed limited parallelism because out of 125 independent loci detected in the Atlantic herring, three showed sharp signals of selection that overlapped between the two species and contained single genes such as PRLRA, which encodes the receptor for prolactin, a freshwater-adapting hormone in euryhaline species, and THRB, a receptor for thyroid hormones, important both for metabolic regulation and the development of red cone photoreceptors.

Discovery of N-[[1-[2-(tert-butylcarbamoylamino)ethyl]-4-(hydroxymethyl)-4-piperidyl]methyl]-3,5-dichloro-benzamide as a selective T-type calcium channel (Cav3.2) inhibitor.

The T-type calcium channel inhibitor Mibefradil was reported to protect the heart from atrial remodeling, a key process involved in the development of atrial fibrillation and arrhythmias. Mibefradil is not a selective T-type calcium channel inhibitor and also affects the function of different ion channels. Our aim was to develop a selective T-type calcium channel inhibitor to validate the importance of T-type-related pharmacology in atrial fibrillation. Structural optimisation of a previously disclosed hit series focussed on minimising exposure to the central nervous system and improving pharmacokinetic properties, while maintain adequate potency and selectivity. This resulted in the design of N-[[1-[2-(tert-butylcarbamoylamino)ethyl]-4-(hydroxymethyl)-4-piperidyl]methyl]-3,5-dichloro-benzamide, a novel, selective, peripherally restricted chemical probe to verify the role of T-type calcium channel inhibition on atrial fibrillation protection.

The Genomic Basis of Adaptation to High Elevations in Africanized Honey Bees.

A range of different genetic architectures underpin local adaptation in nature. Honey bees (Apis mellifera) in the Eastern African Mountains harbor high frequencies of two chromosomal inversions that likely govern adaptation to this high-elevation habitat. In the Americas, honey bees are hybrids of European and African ancestries and adaptation to latitudinal variation in climate correlates with the proportion of these ancestries across the genome. It is unknown which, if either, of these forms of genetic variation governs adaptation in honey bees living at high elevations in the Americas. Here, we performed whole-genome sequencing of 29 honey bees from both high- and low-elevation populations in Colombia. Analysis of genetic ancestry indicated that both populations were predominantly of African ancestry, but the East African inversions were not detected. However, individuals in the higher elevation population had significantly higher proportions of European ancestry, likely reflecting local adaptation. Several genomic regions exhibited particularly high differentiation between highland and lowland bees, containing candidate loci for local adaptation. Genes that were highly differentiated between highland and lowland populations were enriched for functions related to reproduction and sperm competition. Furthermore, variation in levels of European ancestry across the genome was correlated between populations of honey bees in the highland population and populations at higher latitudes in South America. The results are consistent with the hypothesis that adaptation to both latitude and elevation in these hybrid honey bees are mediated by variation in ancestry at many loci across the genome.

Analyses of expressed sequence tags in Neurospora reveal rapid evolution of genes associated with the early stages of sexual reproduction in fungi.

BACKGROUND: The broadly accepted pattern of rapid evolution of reproductive genes is primarily based on studies of animal systems, although several examples of rapidly evolving genes involved in reproduction are found in diverse additional taxa. In fungi, genes involved in mate recognition have been found to evolve rapidly. However, the examples are too few to draw conclusions on a genome scale. RESULTS: In this study, we performed microarray hybridizations between RNA from sexual and vegetative tissues of two strains of the heterothallic (self-sterile) filamentous ascomycete Neurospora intermedia, to identify a set of sex-associated genes in this species. We aligned Expressed Sequence Tags (ESTs) from sexual and vegetative tissue of N. intermedia to orthologs from three closely related species: N. crassa, N. discreta and N. tetrasperma. The resulting four-species alignments provided a dataset for molecular evolutionary analyses. Our results confirm a general pattern of rapid evolution of fungal sex-associated genes, compared to control genes with constitutive expression or a high relative expression during vegetative growth. Among the rapidly evolving sex-associated genes, we identified candidates that could be of importance for mating or fruiting-body development. Analyses of five of these candidate genes from additional species of heterothallic Neurospora revealed that three of them evolve under positive selection. CONCLUSIONS: Taken together, our study represents a novel finding of a genome-wide pattern of rapid evolution of sex-associated genes in the fungal kingdom, and provides a list of candidate genes important for reproductive isolation in Neurospora.

How the worm got its pharynx: phylogeny, classification and Bayesian assessment of character evolution in Acoela.

Acoela are marine microscopic worms currently thought to be the sister taxon of all other bilaterians. Acoels have long been used as models in evolutionary scenarios, and generalized conclusions about acoel and bilaterian ancestral features are frequently drawn from studies of single acoel species. There is no extensive phylogenetic study of Acoela and the taxonomy of the 380 species is chaotic. Here we use two nuclear ribosomal genes and one mitochondrial gene in combination with 37 morphological characters in an analysis of 126 acoel terminals (about one-third of the described species) to estimate the phylogeny and character evolution of Acoela. We present an estimate of posterior probabilities for ancestral character states at 31 control nodes in the phylogeny. The overall reconstruction signal based on the shape of the posterior distribution of character states was computed for all morphological characters and control nodes to assess how well these were reconstructed. The body-wall musculature appears more clearly reconstructed than the reproductive organs. Posterior similarity to the root was calculated by averaging the divergence between the posterior distributions at the nodes and the root over all morphological characters. Diopisthoporidae is the sister group to all other acoels and has the highest posterior similarity to the root. Convolutidae, including several "model" acoels, is most divergent. Finally, we present a phylogenetic classification of Acoela down to the family level where six previous family level taxa are synonymized.

Discovery of 4-morpholino-pyrimidin-6-one and 4-morpholino-pyrimidin-2-one-containing Phosphoinositide 3-kinase (PI3K) p110ß isoform inhibitors through structure-based fragment optimisation.

The discovery of 4-morpholino-pyrimidin-6-one and 4-morpholino-pyrimidin-2-one-containing inhibitors of Phosphoinositide 3-kinases (PI3K) p110ß isoform is reported. Structure-based optimisation of the original fragment hit resulted in lead compounds with improvements in ligand efficiency, lipophilicity efficiency, p110ß potency and selectivity over p110α.

Discovery of phosphoinositide 3-kinases (PI3K) p110ß isoform inhibitor 4-[2-hydroxyethyl(1-naphthylmethyl)amino]-6-[(2S)-2-methylmorpholin-4-yl]-1H-pyrimidin-2-one, an effective antithrombotic agent without associated bleeding and insulin resistance.

Structure-based evolution of the original fragment leads resulted in the identification of 4-[2-hydroxyethyl(1-naphthylmethyl)amino]-6-[(2S)-2-methylmorpholin-4-yl]-1H-pyrimidin-2-one, (S)-21, a potent, selective phosphoinositide 3-kinases (PI3K) p110ß isoform inhibitor with favourable in vivo antiplatelet effect. Despite its antiplatelet action, (S)-21 did not significantly increase bleeding time in dogs. Additionally, due to its enhanced selectivity over p110α, (S)-21 did not induce any insulin resistance in rats.

A comprehensive phylogeny of Neurospora reveals a link between reproductive mode and molecular evolution in fungi.

The filamentous ascomycete genus Neurospora encompasses taxa with a wide range of reproductive modes. Sexual reproduction in this genus can be divided into three major modes; heterothallism (self-incompatibility), homothallism (self-compatibility) and pseudohomothallism (partial self-compatibility). In addition to the sexual pathway, most of the heterothallic taxa propagate with morphologically distinct, vegetative dissemination propagules (macroconidia), while this feature is undetected in the majority of the homothallic taxa. In this study, we used sequence information of seven nuclear gene loci from 43 taxa (295 of the possible 301 locus-by-taxon combinations) to create a phylogeny of Neurospora. The results suggest that transitions in reproductive mode have occurred at multiple times within this group of fungi. Although a homothallic ancestor would imply fewer switches in reproductive mode, we argue that the ancestor of Neurospora was likely heterothallic and that homothallism has evolved independently at least six times in the evolutionary history of the genus. Furthermore, the two pseudohomothallic taxa of Neurospora (N. tetrasperma and N. tetraspora) represent two independent origins of pseudohomothallism. Likelihood ratio tests of substitution rates among branches in the phylogeny indicate that reproductive mode is an important factor driving genome evolution in Neurospora. First, an increased level of non-synonymous/synonymous substitutions in branches delineating homothallic taxa was found, suggesting a reduced efficiency of purifying selection in these taxa. Furthermore, elevated nucleotide substitution rates were found in heterothallic, conidia-producing, lineages as compared to the homothallic non-conidiating lineages. The latter finding is likely due to the presence of conidia, i.e., a higher rate of mitotic divisions inducing mutations, and/or that the homothallic taxa have evolved a lower mutation rate to avoid genomic degeneration.

Substantial Heritable Variation in Recombination Rate on Multiple Scales in Honeybees and Bumblebees.

Meiotic recombination shuffles genetic variation and promotes correct segregation of chromosomes. Rates of recombination vary on several scales, both within genomes and between individuals, and this variation is affected by both genetic and environmental factors. Social insects have extremely high rates of recombination, although the evolutionary causes of this are not known. Here, we estimate rates of crossovers and gene conversions in 22 colonies of the honeybee, Apis mellifera, and 9 colonies of the bumblebee, Bombus terrestris, using direct sequencing of 299 haploid drone offspring. We confirm that both species have extremely elevated crossover rates, with higher rates measured in the highly eusocial honeybee than the primitively social bumblebee. There are also significant differences in recombination rate between subspecies of honeybee. There is substantial variation in genome-wide recombination rate between individuals of both A. mellifera and B. terrestris and the distribution of these rates overlap between species. A large proportion of interindividual variation in recombination rate is heritable, which indicates the presence of variation in trans-acting factors that influence recombination genome-wide. We infer that levels of crossover interference are significantly lower in honeybees compared to bumblebees, which may be one mechanism that contributes to higher recombination rates in honeybees. We also find a significant increase in recombination rate with distance from the centromere, mirrored by methylation differences. We detect a strong transmission bias due to GC-biased gene conversion associated with noncrossover gene conversions. Our results shed light on the mechanistic causes of extreme rates of recombination in social insects and the genetic architecture of recombination rate variation.