The genome sequence of the wood white butterfly, Leptidea sinapis (Linnaeus, 1758)

Meiotic drive against chromosome fusions in butterfly hybrids

Species frequently differ in the number and structure of chromosomes they harbor, but individuals that are heterozygous for chromosomal rearrangements may suffer from reduced fitness. Chromosomal rearrangements like fissions and fusions can hence serve as a mechanism for speciation between incipient lineages, but their evolution poses a paradox. How can rearrangements get fixed between populations if heterozygotes have reduced fitness? One solution is that this process predominantly occurs in small and isolated populations, where genetic drift can override natural selection. However, fixation is also more likely if a novel rearrangement is favored by a transmission bias, such as meiotic drive. Here, we investigate chromosomal transmission distortion in hybrids between two wood white (Leptidea sinapis) butterfly populations with extensive karyotype differences. Using data from two different crossing experiments, we uncover that there is a transmission bias favoring the ancestral chromosomal state for derived fusions, a result that shows that chromosome fusions actually can fix in populations despite being counteracted by meiotic drive. This means that meiotic drive not only can promote runaway chromosome number evolution and speciation, but also that it can be a conservative force acting against karyotypic change and the evolution of reproductive isolation. Based on our results, we suggest a mechanistic model for why chromosome fusion mutations may be opposed by meiotic drive and discuss factors contributing to karyotype evolution in Lepidoptera.

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.

High-density linkage maps and chromosome level genome assemblies unveil direction and frequency of extensive structural rearrangements in wood white butterflies (Leptidea spp.)

Karyotypes are generally conserved between closely related species and large chromosome rearrangements typically have negative fitness consequences in heterozygotes, potentially driving speciation. In the order Lepidoptera, most investigated species have the ancestral karyotype and gene synteny is often conserved across deep divergence, although examples of extensive genome reshuffling have recently been demonstrated. The genus Leptidea has an unusual level of chromosome variation and rearranged sex chromosomes, but the extent of restructuring across the rest of the genome is so far unknown. To explore the genomes of the wood white (Leptidea) species complex, we generated eight genome assemblies using a combination of 10X linked reads and HiC data, and improved them using linkage maps for two populations of the common wood white (L. sinapis) with distinct karyotypes. Synteny analysis revealed an extensive amount of rearrangements, both compared to the ancestral karyotype and between the Leptidea species, where only one of the three Z chromosomes was conserved across all comparisons. Most restructuring was explained by fissions and fusions, while translocations appear relatively rare. We further detected several examples of segregating rearrangement polymorphisms supporting a highly dynamic genome evolution in this clade. Fusion breakpoints were enriched for LINEs and LTR elements, which suggests that ectopic recombination might be an important driver in the formation of new chromosomes. Our results show that chromosome count alone may conceal the extent of genome restructuring and we propose that the amount of genome evolution in Lepidoptera might still be underestimated due to lack of taxonomic sampling.