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

Clonal ants reveal a potentially hidden meiotic feature

Meiosis is essential for eukaryotic reproduction and provides the basis for Mendel's segregation laws. A recent study by Lacy et al. identified a significant deviation from these laws in a clonal ant, hinting at a potentially overlooked meiotic feature. This discovery may have broader implications for recombination in nonclonal eukaryotes.

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.

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.

Recombination mapping of the Brazilian stingless bee Frieseomelitta varia confirms high recombination rates in social hymenoptera

Background

Meiotic recombination is a fundamental genetic process that shuffles allele combinations and promotes accurate segregation of chromosomes. Analyses of the ubiquitous variation of recombination rates within and across species suggest that recombination is evolving adaptively. All studied insects with advanced eusociality have shown exceptionally high recombination rates, which may represent a prominent case of adaptive evolution of recombination. However, our understanding of the relationship between social evolution and recombination rates is incomplete, partly due to lacking empirical data. Here, we present a linkage map of the monandrous, advanced eusocial Brazilian stingless bee, Frieseomelitta varia, providing the first recombination analysis in the diverse Meliponini (Hymenoptera, Apidae).

Results

Our linkage map includes 1417 markers in 19 linkage groups. This map spans approximately 2580 centimorgans, and comparisons to the physical genome assembly indicate that it covers more than 75 % of the 275 Megabasepairs (Mbp) F. varia genome. Thus, our study results in a genome-wide recombination rate estimate of 9.3–12.5 centimorgan per Mbp. This value is higher than estimates from nonsocial insects and comparable to other highly social species, although it does not support our prediction that monandry and strong queen-worker caste divergence of F. varia lead to even higher recombination rates than other advanced eusocial species.

Conclusions

Our study expands the association between elevated recombination and sociality in the order Hymenoptera and strengthens the support for the hypothesis that advanced social evolution in hymenopteran insects invariably selects for high genomic recombination rates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07987-3.

Eusociality in snapping shrimps is associated with larger genomes and an accumulation of transposable elements

Despite progress uncovering the genomic underpinnings of sociality, much less is known about how social living affects the genome. In different insect lineages, for example, eusocial species show both positive and negative associations between genome size and structure, highlighting the dynamic nature of the genome. Here, we explore the relationship between sociality and genome architecture in Synalpheus snapping shrimps that exhibit multiple origins of eusociality and extreme interspecific variation in genome size. Our goal is to determine whether eusociality leads to an accumulation of repetitive elements and an increase in genome size, presumably due to reduced effective population sizes resulting from a reproductive division of labor, or whether an initial accumulation of repetitive elements leads to larger genomes and independently promotes the evolution of eusociality through adaptive evolution. Using phylogenetically informed analyses, we find that eusocial species have larger genomes with more transposable elements (TEs) and microsatellite repeats than noneusocial species. Interestingly, different TE subclasses contribute to the accumulation in different species. Phylogenetic path analysis testing alternative causal relationships between sociality and genome architecture is most consistent with the hypothesis that TEs modulate the relationship between sociality and genome architecture. Although eusociality appears to influence TE accumulation, ancestral state reconstruction suggests moderate TE abundances in ancestral species could have fueled the initial transitions to eusociality. Ultimately, we highlight a complex and dynamic relationship between genome and social evolution, demonstrating that sociality can influence the evolution of the genome, likely through changes in demography related to patterns of reproductive skew.

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.

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.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.

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.

The impact of GC bias on phylogenetic accuracy using targeted enrichment phylogenomic data

The field of sequence based phylogenetic analyses is currently being transformed by novel hybrid-based targeted enrichment methods, such as the use of ultraconserved elements (UCEs). Rather than analyzing relationships among organisms using a small number of genes, these methods now allow us to evaluate relationships with many hundreds to thousands of individual gene loci. However, the inclusion of thousands of loci does not necessarily overcome the long-standing challenge of incongruence among phylogenetic trees derived from different genes or gene regions. One factor that impacts the level of incongruence in phylogenomic data sets is the level of GC bias. GC rich gene regions are prone to higher recombination rates than AT rich regions, driven by a process referred to as "GC biased gene conversion". As a result, high GC content can be negatively associated with phylogenetic accuracy, but the extent to which this impacts incongruence among UCEs is currently unstudied. We investigated the impact of GC content on phylogeny reconstruction using in silico captured UCE data for the corbiculate bees (Hymenoptera: Apidae). The phylogeny of this group has been the subject of extensive study, and incongruence among gene trees is thought to be a source of phylogenetic error. We conducted coalescent- and concatenation-based analyses of 810 individual gene loci from all 13 currently available bee genomes, including 8 corbiculate taxa. Both coalescent- and concatenation-based methods converged on a single topology for the corbiculate tribes. In contrast to concatenation, the coalescent-based methods revealed significant topological conflict at nodes involving the orchid bees (Euglossini) and honeybees (Apini). Partitioning the loci by GC content reveals decreasing support for the inferred topology with increasing GC bias. Based on the results of this study, we report the first evidence that GC biased gene conversion may contribute to topological incongruence in studies based on ultraconserved elements.

The Genome and Methylome of a Subsocial Small Carpenter Bee, Ceratina calcarata

Understanding the evolution of animal societies, considered to be a major transition in evolution, is a key topic in evolutionary biology. Recently, new gateways for understanding social evolution have opened up due to advances in genomics, allowing for unprecedented opportunities in studying social behavior on a molecular level. In particular, highly eusocial insect species (caste-containing societies with nonreproductives that care for siblings) have taken center stage in studies of the molecular evolution of sociality. Despite advances in genomic studies of both solitary and eusocial insects, we still lack genomic resources for early insect societies. To study the genetic basis of social traits requires comparison of genomes from a diversity of organisms ranging from solitary to complex social forms. Here we present the genome of a subsocial bee, Ceratina calcarata This study begins to address the types of genomic changes associated with the earliest origins of simple sociality using the small carpenter bee. Genes associated with lipid transport and DNA recombination have undergone positive selection in C. calcarata relative to other bee lineages. Furthermore, we provide the first methylome of a noneusocial bee. Ceratina calcarata contains the complete enzymatic toolkit for DNA methylation. As in the honey bee and many other holometabolous insects, DNA methylation is targeted to exons. The addition of this genome allows for new lines of research into the genetic and epigenetic precursors to complex social behaviors.

Genomic sources of phenotypic novelty in the evolution of eusociality in insects

Genomic resources are now available for closely related species that vary in social behavior, providing insight on the genomics of social evolution. Changes in the architecture of gene regulatory networks likely influence the evolutionary trajectory of social traits. Evolutionarily novel genes are likely important in the evolution of social diversity among insects, but it is unclear whether new genes played a driving role in the advent or elaboration of eusociality or if they were instead a result of other genomic features of eusociality. The worker phenotype appears to be the center of genetic novelty, but the mechanisms for this remain unresolved. Future studies are needed to understand how genetic novelty arises, becomes incorporated into existing gene regulatory networks, and the effects this has on the evolution of social traits in closely related social and solitary species.

Climbing the social ladder: the molecular evolution of sociality

Genomic tools are allowing us to dissect the roles of genes and genetic architecture in social evolution, and eusocial insects are excellent models. Numerous hypotheses for molecular evolution of eusociality have been proposed, ranging from regulatory shifts in 'old' genes to rapid evolution of 'new' genes. A broad model to explain this major transition in evolution has been lacking. We provide a synthetic framework centered on the idea that different evolutionary processes dominate during different transitional stages, beginning with changes in gene regulation and culminating in novel genes later on. By considering multiple mechanisms as we 'climb the social ladder', we can test whether the transitions from solitary to simple sociality to complex sociality represent incremental changes or genetic revolutions.

Complementary sex determination in the parasitic wasp Diachasmimorpha longicaudata

We studied the sex determination in Diachasmimorpha longicaudata, a parasitoid braconid wasp widely used as biological control agent of fruit pest tephritid flies. We tested the complementary sex determination hypothesis (CSD) known in at least 60 species of Hymenoptera. According to CSD, male or female development depends on the allelic composition of one sex locus (single-locus CSD) or multiple sex loci (multiple-locus CSD). Hemizygote individuals are normal haploid males, and heterozygotes for at least one sex locus are normal diploid females, but homozygotes for all the sex loci are diploid males. In order to force the occurrence of diploid males in D. longicaudata, we established highly inbred lines and examined their offspring using chromosome counting, flow cytometry, and sex ratio analysis. We found that when mother-son crosses were studied, this wasp produced about 20% of diploid males out of the total male progeny. Our results suggest that this parasitoid may represent the second genus with multiple-locus CSD in Hymenoptera. Knowledge about the sex determination system in D. longicaudata is relevant for the improvement of mass rearing protocols of this species. This information also provides the necessary background for further investigations on the underlying molecular mechanisms of sex determination in this species, and a better insight into the evolution of this pathway in Hymenoptera in particular and insects in general.

Genomic correlates of recombination rate and its variability across eight recombination maps in the western honey bee (Apis mellifera L.)

Background

Meiotic recombination has traditionally been explained based on the structural requirement to stabilize homologous chromosome pairs to ensure their proper meiotic segregation. Competing hypotheses seek to explain the emerging findings of significant heterogeneity in recombination rates within and between genomes, but intraspecific comparisons of genome-wide recombination patterns are rare. The honey bee (Apis mellifera) exhibits the highest rate of genomic recombination among multicellular animals with about five cross-over events per chromatid.

Results

Here, we present a comparative analysis of recombination rates across eight genetic linkage maps of the honey bee genome to investigate which genomic sequence features are correlated with recombination rate and with its variation across the eight data sets, ranging in average marker spacing ranging from 1 Mbp to 120 kbp. Overall, we found that GC content explained best the variation in local recombination rate along chromosomes at the analyzed 100 kbp scale. In contrast, variation among the different maps was correlated to the abundance of microsatellites and several specific tri- and tetra-nucleotides.

Conclusions

The combined evidence from eight medium-scale recombination maps of the honey bee genome suggests that recombination rate variation in this highly recombining genome might be due to the DNA configuration instead of distinct sequence motifs. However, more fine-scale analyses are needed. The empirical basis of eight differing genetic maps allowed for robust conclusions about the correlates of the local recombination rates and enabled the study of the relation between DNA features and variability in local recombination rates, which is particularly relevant in the honey bee genome with its exceptionally high recombination rate.

Electronic supplementary material

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

Recombination, chromosome number and eusociality in the Hymenoptera

Extraordinarily high rates of recombination have been observed in some eusocial species. The most popular explanation is that increased recombination increases genetic variation among workers, which in turn increases colony performance, for example by increasing parasite resistance. However, support for the generality of higher recombination rates among eusocial organisms remains weak, due to low sample size and a lack of phylogenetic independence of observations. Recombination rate, although difficult to measure directly, is correlated with chromosome number. As predicted, several authors have noted that chromosome numbers are higher among the eusocial species of Hymenoptera (ants, bees and wasps). Here, we present a formal comparative analysis of karyotype data from 1567 species of Hymenoptera. Contrary to earlier studies, we find no evidence for an absolute difference between chromosome number in eusocial and solitary species of Hymenoptera. However, we find support for an increased rate of chromosome number change in eusocial taxa. We show that among eusocial taxa colony size is able to explain some of the variation in chromosome number: intermediate-sized colonies have more chromosomes than those that are either very small or very large. However, we were unable to detect effects of a number of other colony characteristics predicted to affect recombination rate - including colony relatedness and caste number. Taken together, our results support the view that a eusocial lifestyle has led to variable selection pressure for increased recombination rates, but that identifying the factors contributing to this variable selection will require further theoretical and empirical effort.

Causes and consequences of crossing-over evidenced via a high-resolution recombinational landscape of the honey bee

Background

Social hymenoptera, the honey bee (Apis mellifera) in particular, have ultra-high crossover rates and a large degree of intra-genomic variation in crossover rates. Aligned with haploid genomics of males, this makes them a potential model for examining the causes and consequences of crossing over. To address why social insects have such high crossing-over rates and the consequences of this, we constructed a high-resolution recombination atlas by sequencing 55 individuals from three colonies with an average marker density of 314 bp/marker.

Results

We find crossing over to be especially high in proximity to genes upregulated in worker brains, but see no evidence for a coupling with immune-related functioning. We detect only a low rate of non-crossover gene conversion, contrary to current evidence. This is in striking contrast to the ultrahigh crossing-over rate, almost double that previously estimated from lower resolution data. We robustly recover the predicted intragenomic correlations between crossing over and both population level diversity and GC content, which could be best explained as indirect and direct consequences of crossing over, respectively.

Conclusions

Our data are consistent with the view that diversification of worker behavior, but not immune function, is a driver of the high crossing-over rate in bees. While we see both high diversity and high GC content associated with high crossing-over rates, our estimate of the low non-crossover rate demonstrates that high non-crossover rates are not a necessary consequence of high recombination rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0566-0) contains supplementary material, which is available to authorized users.

Interindividual variability in social insects - proximate causes and ultimate consequences

Individuals within social groups often show consistent differences in behaviour across time and context. Such interindividual differences and the evolutionary challenge they present have recently generated considerable interest. Social insects provide some of the most familiar and spectacular examples of social groups with large interindividual differences. Investigating these within-group differences has a long research tradition, and behavioural variability among the workers of a colony is increasingly regarded as fundamental for a key feature of social insects: division of labour. The goal of this review is to illustrate what we know about both the proximate mechanisms underlying behavioural variability among the workers of a colony and its ultimate consequences; and to highlight the many open questions in this research field. We begin by reviewing the literature on mechanisms that potentially introduce, maintain, and adjust the behavioural differentiation among workers. We highlight the fact that so far, most studies have focused on behavioural variability based on genetic variability, provided by e.g. multiple mating of the queen, while other mechanisms that may be responsible for the behavioural differentiation among workers have been largely neglected. These include maturational, nutritional and environmental influences. We further discuss how feedback provided by the social environment and learning and experience of adult workers provides potent and little-explored sources of differentiation. In a second part, we address what is known about the potential benefits and costs of increased behavioural variability within the workers of a colony. We argue that all studies documenting a benefit of variability so far have done so by manipulating genetic variability, and that a direct test of the effect of behavioural variability on colony productivity has yet to be provided. We emphasize that the costs associated with interindividual variability have been largely overlooked, and that a better knowledge of the cost/benefit balance of behavioural variability is crucial for our understanding of the evolution of the mechanisms underlying the social organization of insect societies. We conclude by highlighting what we believe to be promising but little-explored avenues for future research on how within-colony variability has evolved and is maintained. We emphasize the need for comparative studies and point out that, so far, most studies on interindividual variability have focused on variability in individual response thresholds, while the significance of variability in other parameters of individual response, such as probability and intensity of the response, has been largely overlooked. We propose that these parameters have important consequences for the colony response. Much more research is needed to understand if and how interindividual variability is modulated in order to benefit division of labour, homeostasis and ultimately colony fitness in social insects.

The ant genomes have been invaded by several types of mariner transposable elements

To date, only three types of full-length mariner elements have been described in ants, each one in a different genus of the Myrmicinae subfamily: Sinvmar was isolated from various Solenopsis species, Myrmar from Myrmica ruginodis, and Mboumar from Messor bouvieri. In this study, we report the coexistence of three mariner elements (Tnigmar-Si, Tnigmar-Mr, and Tnigmar-Mb) in the genome of a single species, Tapinoma nigerrimum (subfamily Dolichoderinae). Molecular evolutionary analyses of the nucleotide sequence data revealed a general agreement between the evolutionary history of most the elements and the ant species that harbour them, and suggest that they are at the vertical inactivation stage of the so-called Mariner Life Cycle. In contrast, significantly reduced levels of synonymous divergence between Mboumar and Tnigmar-Mb and between Myrmar and Botmar (a mariner element isolated from Bombus terrestris), relative to those observed between their hosts, suggest that these elements arrived to the species that host them by horizontal transfer, long after the species' split. The horizontal transfer events for the two pairs of elements could be roughly dated within the last 2 million years and about 14 million years, respectively. As would be expected under this scenario, the coding sequences of the youngest elements, Tnigmar-Mb and Mboumar, are intact and, thus, potentially functional. Each mariner element has a different chromosomal distribution pattern according to their stage within the Mariner Life Cycle. Finally, a new defective transposable element (Azteca) has also been found inserted into the Tnigmar-Mr sequences showing that the ant genomes have been invaded by at least four different types of mariner elements.

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.

Recombination is associated with the evolution of genome structure and worker behavior in honey bees

The rise of insect societies, marked by the formation of reproductive and sterile castes, represents a major unsolved mystery in evolution. Across several independent origins of sociality, the genomes of social hymenopterans share two peculiar attributes: high recombination and low but heterogeneous GC content. For example, the genome of the honey bee, Apis mellifera, represents a mosaic of GC-poor and GC-rich regions with rates of recombination an order of magnitude higher than in humans. However, it is unclear how heterogeneity in GC content arises, and how it relates to the expression and evolution of worker traits. Using population genetic analyses, we demonstrate a bias in the allele frequency and fixation rate of derived C or G mutations in high-recombination regions, consistent with recombination's causal influence on GC-content evolution via biased gene conversion. We also show that recombination and biased gene conversion actively maintain the heterogeneous GC content of the honey bee genome despite an overall A/T mutation bias. Further, we found that GC-rich genes and intergenic regions have higher levels of genetic diversity and divergence relative to GC-poor regions, also consistent with recombination's causal influence on the rate of molecular evolution. Finally, we found that genes associated with behavior and those with worker-biased expression are found in GC-rich regions of the bee genome and also experience high rates of molecular evolution. Taken together, these findings suggest that recombination acts to maintain a genetically diverse and dynamic part of the genome where genes underlying worker behavior evolve more quickly.