Alternative double strand break repair pathways shape the evolution of high recombination in the honey bee, <scp> <i>Apis mellifera</i> </scp>

Social insects, particularly honey bees, have exceptionally high genomic frequencies of genetic recombination. This phenomenon and underlying mechanisms are poorly understood. To characterise the patterns of crossovers and gene conversion in the honey bee genome, a recombination map of 187 honey bee brothers was generated by whole-genome resequencing. Recombination events were heterogeneously distributed without many true hotspots. The tract lengths between phase shifts were bimodally distributed, indicating distinct crossover and gene conversion events. While crossovers predominantly occurred in G/C-rich regions and seemed to cause G/C enrichment, the gene conversions were found predominantly in A/T-rich regions. The nucleotide composition of sequences involved in gene conversions that were associated with or distant from crossovers corresponded to the differences between crossovers and gene conversions. These combined results suggest two types of DNA double-strand break repair during honey bee meiosis: non-canonical homologous recombination, leading to gene conversion and A/T enrichment of the genome, and the canonical homologous recombination based on completed double Holliday Junctions, which can result in gene conversion or crossover and is associated with G/C bias. This G/C bias may be selected for to balance the A/T-rich base composition of eusocial hymenopteran genomes. The lack of evidence for a preference of the canonical homologous recombination for double-strand break repair suggests that the high genomic recombination rate of honey bees is mainly the consequence of a high rate of double-strand breaks, which could in turn result from the life history of honey bees and their A/T-rich genome.