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

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

Modeling Recombination Rate as a Quantitative Trait Reveals New Insight into Selection in Humans

Meiotic recombination is both a fundamental biological process required for proper chromosomal segregation during meiosis and an important genomic parameter that shapes major features of the genomic landscape. However, despite the central importance of this phenotype, we lack a clear understanding of the selective pressures that shape its variation in natural populations, including humans. While there is strong evidence of fitness costs of low rates of recombination, the possible fitness costs of high rates of recombination are less defined. To determine whether a single lower fitness bound can explain the variation in recombination rates observed in human populations, we simulated the evolution of recombination rates as a sexually dimorphic quantitative trait. Under each scenario, we statistically compared the resulting trait distribution with the observed distribution of recombination rates from a published study of the Icelandic population. To capture the genetic architecture of recombination rates in humans, we modeled it as a moderately complex trait with modest heritability. For our fitness function, we implemented a hyperbolic tangent curve with several flexible parameters to capture a wide range of existing hypotheses. We found that costs of low rates of recombination alone are likely insufficient to explain the current variation in recombination rates in both males and females, supporting the existence of fitness costs of high rates of recombination in humans. With simulations using both upper and lower fitness boundaries, we describe a parameter space for the costs of high recombination rates that produces results consistent with empirical observations.

Base Composition, Codon Usage, and Patterns of Gene Sequence Evolution in Butterflies

Coding sequence evolution is influenced by both natural selection and neutral evolutionary forces. In many species, the effects of mutation bias, codon usage, and GC-biased gene conversion (gBGC) on gene sequence evolution have not been detailed. Quantification of how these forces shape substitution patterns is therefore necessary to understand the strength and direction of natural selection. Here, we used comparative genomics to investigate the association between base composition and codon usage bias on gene sequence evolution in butterflies and moths (Lepidoptera), including an in-depth analysis of underlying patterns and processes in one species, Leptidea sinapis. The data revealed significant G/C to A/T substitution bias at third codon position with some variation in the strength among different butterfly lineages. However, the substitution bias was lower than expected from previously estimated mutation rate ratios, partly due to the influence of gBGC. We found that A/T-ending codons were overrepresented in most species, but there was a positive association between the magnitude of codon usage bias and GC-content in third codon positions. In addition, the tRNA-gene population in L. sinapis showed higher GC-content at third codon positions compared to coding sequences in general and less overrepresentation of A/T-ending codons. There was an inverse relationship between synonymous substitutions and codon usage bias indicating selection on synonymous sites. We conclude that the evolutionary rate in Lepidoptera is affected by a complex interaction between underlying G/C -> A/T mutation bias and partly counteracting fixation biases, predominantly conferred by overall purifying selection, gBGC, and selection on codon usage.

AmelHap: Leveraging drone whole-genome sequence data to create a honey bee HapMap

Honey bee, Apis mellifera, drones are typically haploid, developing from an unfertilized egg, inheriting only their queen's alleles and none from the many drones she mated with. Thus the ordered combination or 'phase' of alleles is known, making drones a valuable haplotype resource. We collated whole-genome sequence data for 1,407 drones, including 45 newly sequenced Scottish drones, collectively representing 19 countries, 8 subspecies and various hybrids. Following alignment to Amel_HAv3.1, variant calling and quality filtering, we retained 17.4 M high quality variants across 1,328 samples with a genotyping rate of 98.7%. We demonstrate the utility of this haplotype resource, AmelHap, for genotype imputation, returning >95% concordance when up to 61% of data is missing in haploids and up to 12% of data is missing in diploids. AmelHap will serve as a useful resource for the community for imputation from low-depth sequencing or SNP chip data, accurate phasing of diploids for association studies, and as a comprehensive reference panel for population genetic and evolutionary analyses.

Complex population structure and haplotype patterns in the Western European honey bee from sequencing a large panel of haploid drones

Honey bee subspecies originate from specific geographical areas in Africa, Europe and the Middle East, and beekeepers interested in specific phenotypes have imported genetic material to regions outside of the bees' original range for use either in pure lines or controlled crosses. Moreover, imported drones are present in the environment and mate naturally with queens from the local subspecies. The resulting admixture complicates population genetics analyses, and population stratification can be a major problem for association studies. To better understand Western European honey bee populations, we produced a whole genome sequence and single nucleotide polymorphism (SNP) genotype data set from 870 haploid drones and demonstrate its utility for the identification of nine genetic backgrounds and various degrees of admixture in a subset of 629 samples. Five backgrounds identified correspond to subspecies, two to isolated populations on islands and two to managed populations. We also highlight several large haplotype blocks, some of which coincide with the position of centromeres. The largest is 3.6 Mb long and represents 21% of chromosome 11, with two major haplotypes corresponding to the two dominant genetic backgrounds identified. This large naturally phased data set is available as a single vcf file that can now serve as a reference for subsequent populations genomics studies in the honey bee, such as (i) selecting individuals of verified homogeneous genetic backgrounds as references, (ii) imputing genotypes from a lower-density data set generated by an SNP-chip or by low-pass sequencing, or (iii) selecting SNPs compatible with the requirements of genotyping chips.

Fragile, unfaithful and persistent Ys-on how meiosis can shape sex chromosome evolution

Sex-linked inheritance is a stark exception to Mendel's Laws of Heredity. Here we discuss how the evolution of heteromorphic sex chromosomes (mainly the Y) has been shaped by the intricacies of the meiotic programme. We propose that persistence of Y chromosomes in distantly related mammalian phylogroups can be explained in the context of pseudoautosomal region (PAR) size, meiotic pairing strategies, and the presence of Y-borne executioner genes that regulate meiotic sex chromosome inactivation. We hypothesise that variation in PAR size can be an important driver for the evolution of recombination frequencies genome wide, imposing constraints on Y fate. If small PAR size compromises XY segregation during male meiosis, the stress of producing aneuploid gametes could drive function away from the Y (i.e., a fragile Y). The Y chromosome can avoid fragility either by acquiring an achiasmatic meiotic XY pairing strategy to reduce aneuploid gamete production, or gain meiotic executioner protection (a persistent Y). Persistent Ys will then be under strong pressure to maintain high recombination rates in the PAR (and subsequently genome wide), as improper segregation has fatal consequences for germ cells. In the event that executioner protection is lost, the Y chromosome can be maintained in the population by either PAR rejuvenation (extension by addition of autosome material) or gaining achiasmatic meiotic pairing, the alternative is Y loss. Under this dynamic cyclic evolutionary scenario, understanding the meiotic programme in vertebrate and invertebrate species will be crucial to further understand the plasticity of the rise and fall of heteromorphic sex chromosomes.

A Novel Strategy to Reveal the Landscape of Crossovers in an F1 Hybrid Population of Populus deltoides and Populus simonii

Although the crossover (CO) patterns of different species have been extensively investigated, little is known about the landscape of CO patterns in Populus because of its high heterozygosity and long-time generation. A novel strategy was proposed to reveal the difference of CO rate and interference between Populus deltoides and Populus simonii using their F1 hybrid population. We chose restriction site-associated DNA (RAD) tags that contained two SNPs, one only receiving the CO information from the female P. deltoides and the other from the male P. simonii. These RAD tags allowed us to investigate the CO patterns between the two outbred species, instead of using the traditional backcross populations in inbred lines. We found that the CO rate in P. deltoides was generally greater than that in P. simonii, and that the CO interference was a common phenomenon across the two genomes. The COs landscape of the different Populus species facilitates not only to understand the evolutionary mechanism for adaptability but also to rebuild the statistical model for precisely constructing genetic linkage maps that are critical in genome assembly in Populus. Additionally, the novel strategy could be applied in other outbred species for investigating the CO patterns.

A genomic and morphometric analysis of alpine bumblebees: Ongoing reductions in tongue length but no clear genetic component

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyse population structure in these populations. We generated a genome assembly of B. balteatus. We then analysed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait probably has a complex genetic and environmental basis.

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.

The Effects of GC-Biased Gene Conversion on Patterns of Genetic Diversity among and across Butterfly Genomes

Recombination reshuffles the alleles of a population through crossover and gene conversion. These mechanisms have considerable consequences on the evolution and maintenance of genetic diversity. Crossover, for example, can increase genetic diversity by breaking the linkage between selected and nearby neutral variants. Bias in favor of G or C alleles during gene conversion may instead promote the fixation of one allele over the other, thus decreasing diversity. Mutation bias from G or C to A and T opposes GC-biased gene conversion (gBGC). Less recognized is that these two processes may-when balanced-promote genetic diversity. Here, we investigate how gBGC and mutation bias shape genetic diversity patterns in wood white butterflies (Leptidea sp.). This constitutes the first in-depth investigation of gBGC in butterflies. Using 60 resequenced genomes from six populations of three species, we find substantial variation in the strength of gBGC across lineages. When modeling the balance of gBGC and mutation bias and comparing analytical results with empirical data, we reject gBGC as the main determinant of genetic diversity in these butterfly species. As alternatives, we consider linked selection and GC content. We find evidence that high values of both reduce diversity. We also show that the joint effects of gBGC and mutation bias can give rise to a diversity pattern which resembles the signature of linked selection. Consequently, gBGC should be considered when interpreting the effects of linked selection on levels of genetic diversity.

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (d_XY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.

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.