Evolution and genetic architecture of disassortative mating at a locus under heterozygote advantage

The interplay of local adaptation and gene flow may lead to the formation of supergenes

Supergenes are genetic architectures resulting in the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of alleles at linked loci is often facilitated by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, colour or behavioural polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the formation of supergenes remains poorly understood. Here, using a two-island population genetic model, we explore how gene flow and the evolution of overdominant chromosomal inversions may jointly lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations, both under disruptive selection, leads to an increase in frequency of poorly adapted, immigrant haplotypes. Indeed, rare allelic combinations, such as immigrant haplotypes, are more frequently reshuffled by recombination than common allelic combinations, and therefore benefit from the recombination suppression generated by inversions. When an inversion capturing a locally adapted haplotype spreads but is associated with a fitness cost hampering its fixation (e.g. a recessive mutation load), the maintenance of a non-inverted haplotype in the population is enhanced; under certain conditions, the immigrant haplotype persists alongside the inverted local haplotype, while the standard local haplotype disappears. This establishes a stable, local polymorphism with two non-recombining haplotypes encoding alternative adaptive strategies, that is, a supergene. These results bring new light to the importance of local adaptation, overdominance, and gene flow in the formation of supergenes and inversion polymorphisms in general.

How chromosomal inversions reorient the evolutionary process

Inversions are structural mutations that reverse the sequence of a chromosome segment and reduce the effective rate of recombination in the heterozygous state. They play a major role in adaptation, as well as in other evolutionary processes such as speciation. Although inversions have been studied since the 1920s, they remain difficult to investigate because the reduced recombination conferred by them strengthens the effects of drift and hitchhiking, which in turn can obscure signatures of selection. Nonetheless, numerous inversions have been found to be under selection. Given recent advances in population genetic theory and empirical study, here we review how different mechanisms of selection affect the evolution of inversions. A key difference between inversions and other mutations, such as single nucleotide variants, is that the fitness of an inversion may be affected by a larger number of frequently interacting processes. This considerably complicates the analysis of the causes underlying the evolution of inversions. We discuss the extent to which these mechanisms can be disentangled, and by which approach.

A supergene in seaweed flies modulates male traits and female perception

Supergenes, tightly linked sets of alleles, offer some of the most spectacular examples of polymorphism persisting under long-term balancing selection. However, we still do not understand their evolution and persistence, especially in the face of accumulation of deleterious elements. Here, we show that an overdominant supergene in seaweed flies, Coelopa frigida, modulates male traits, potentially facilitating disassortative mating and promoting intraspecific polymorphism. Across two continents, the Cf-Inv(1) supergene strongly affected the composition of male cuticular hydrocarbons (CHCs) but only weakly affected CHC composition in females. Using gas chromatography-electroantennographic detection, we show that females can sense male CHCs and that there may be differential perception between genotypes. Combining our phenotypic results with RNA-seq data, we show that candidate genes for CHC biosynthesis primarily show differential expression for Cf-Inv(1) in males but not females. Conversely, candidate genes for odorant detection were differentially expressed in both sexes but showed high levels of divergence between supergene haplotypes. We suggest that the reduced recombination between supergene haplotypes may have led to rapid divergence in mate preferences as well as increasing linkage between male traits, and overdominant loci. Together this probably helped to maintain the polymorphism despite deleterious effects in homozygotes.

Genotype and growth rate influence female mate preference in Xiphophorus multilineatus: Potential selection to optimize mortality-growth rate tradeoff

The extent to which mate preferences are adaptive requires a better understanding of the factors that influence variation in mate preferences. Xiphophorus multilineatus is a live-bearing fish with males that exhibit alternative reproductive tactics (courter/sneaker). We examined the influence of a female's genotype (courter vs sneaker lineage), growth rate, and social experience on mate preference for courter as compared to sneaker males. We found that females with a sneaker genotype and slower growth rates had stronger mate preferences for the faster growing courter males than females with a courter genotype, regardless of mating experience with one or both types of males. In addition, the relationship between strength of preference and growth rate depended on a females' genotype; females with sneaker genotypes decreased their preference as their growth rates increased, a pattern that trended in the opposite direction for females from the courter genotypes. Disassortative mating preferences are predicted to evolve when heterozygous offspring benefit from increased fitness. Given male tactical dimorphism in growth rates and a mortality-growth rate tradeoff previously detected in this species, the variation in mating preferences for the male tactics we detected may be under selection to optimize the mortality-growth rate tradeoff for offspring.

Cryptic recessive lethality of a supergene controlling social organization in ants

Supergenes are clusters of linked loci that control complex phenotypes, such as alternative forms of social organization in ants. Explaining the long-term maintenance of supergenes is challenging, particularly when the derived haplotype lacks homozygous lethality and causes gene drive. In the Alpine silver ant, Formica selysi, a large and ancient social supergene with two haplotypes, M and P, controls colony social organization. Single-queen colonies only contain MM females, while multiqueen colonies contain MP and PP females. The derived P haplotype, found only in multiqueen colonies, selfishly enhances its transmission through maternal effect killing, which could have led to its fixation. A population genetic model showed that a stable social polymorphism can only be maintained under a narrow set of conditions, which includes partial assortative mating by social form (which is known to occur in the wild), and low fitness of PP queens. With a combination of field and laboratory experiments, we show that the P haplotype has deleterious effects on female fitness. The survival rate of PP queens and workers was around half that of other genotypes. Moreover, P-carrying queens had lower fertility and fecundity compared to other queens. We discuss how cryptic lethal effects of the P haplotype help stabilize this ancient polymorphism.

Association mapping of colour variation in a butterfly provides evidence that a supergene locks together a cluster of adaptive loci

Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between sets of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining a landmark-free phenotyping algorithm with multivariate genome-wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We show that the supergene controlling the striking wing pattern polymorphism displayed by this species contains several independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are, however, associated with colour variations only in a subset of morphs where the phenotype is controlled by derived inversion forms, indicating that they were recruited after the formation of the inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Mutation accumulation opposes polymorphism: supergenes and the curious case of balanced lethals

Supergenes offer spectacular examples of long-term balancing selection in nature, but their origin and maintenance remain a mystery. Reduced recombination between arrangements, a critical aspect of many supergenes, protects adaptive multi-trait phenotypes but can lead to mutation accumulation. Mutation accumulation can stabilize the system through the emergence of associative overdominance (AOD), destabilize the system, or lead to new evolutionary outcomes. One outcome is the formation of maladaptive balanced lethal systems, where only heterozygotes remain viable and reproduce. We investigated the conditions under which these different outcomes occur, assuming a scenario of introgression after divergence. We found that AOD aided the invasion of a new supergene arrangement and the establishment of a polymorphism. However, this polymorphism was easily destabilized by further mutation accumulation, which was often asymmetric, disrupting the quasi-equilibrium state. Mechanisms that accelerated degeneration tended to amplify asymmetric mutation accumulation between the supergene arrangements and vice-versa. As the evolution of balanced lethal systems requires symmetric degeneration of both arrangements, this leaves only restricted conditions for their evolution, namely small population sizes and low rates of gene conversion. The dichotomy between the persistence of polymorphism and degeneration of supergene arrangements likely underlies the rarity of balanced lethal systems in nature.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

A large deletion at the cortex locus eliminates butterfly wing patterning

As the genetic basis of natural and domesticated variation has been described in recent years, a number of hotspot genes have been repeatedly identified as the targets of selection, Heliconius butterflies display a spectacular diversity of pattern variants in the wild and the genetic basis of these patterns has been well-described. Here, we sought to identify the mechanism behind an unusual pattern variant that is instead found in captivity, the ivory mutant, in which all scales on both the wings and body become white or yellow. Using a combination of autozygosity mapping and coverage analysis from 37 captive individuals, we identify a 78-kb deletion at the cortex wing patterning locus, a gene which has been associated with wing pattern evolution in H. melpomene and 10 divergent lepidopteran species. This deletion is undetected among 458 wild Heliconius genomes samples, and its dosage explains both homozygous and heterozygous ivory phenotypes found in captivity. The deletion spans a large 5′ region of the cortex gene that includes a facultative 5′UTR exon detected in larval wing disk transcriptomes. CRISPR mutagenesis of this exon replicates the wing phenotypes from coding knock-outs of cortex, consistent with a functional role of ivory-deleted elements in establishing scale color fate. Population demographics reveal that the stock giving rise to the ivory mutant has a mixed origin from across the wild range of H. melpomene, and supports a scenario where the ivory mutation occurred after the introduction of cortex haplotypes from Ecuador. Homozygotes for the ivory deletion are inviable while heterozygotes are the targets of artificial selection, joining 40 other examples of allelic variants that provide heterozygous advantage in animal populations under artificial selection by fanciers and breeders. Finally, our results highlight the promise of autozygosity and association mapping for identifying the genetic basis of aberrant mutations in captive insect populations.

Balanced polymorphisms and their divergence in a Heliconius butterfly

The evolution of mimicry in similarly defended prey is well described by the Müllerian mimicry theory, which predicts the convergence of warning patterns in order to gain the most protection from predators. However, despite this prediction, we can find great diversity of color patterns among Müllerian mimics such as Heliconius butterflies in the neotropics. Furthermore, some species have evolved the ability to maintain multiple distinct warning patterns in single populations, a phenomenon known as polymorphic mimicry. The adaptive benefit of these polymorphisms is questionable since variation from the most common warning patterns is expected to be disadvantageous as novel signals are punished by predators naive to them. In this study, we use artificial butterfly models throughout Central and South America to characterize the selective pressures maintaining polymorphic mimicry in Heliconius doris. Our results highlight the complexity of positive frequency-dependent selection, the principal selective pressure driving convergence among Müllerian mimics, and its impacts on interspecific variation of mimetic warning coloration. We further show how this selection regime can both limit and facilitate the diversification of mimetic traits.

When Do Opposites Attract? A Model Uncovering the Evolution of Disassortative Mating

AbstractDisassortative mating is a rare form of mate preference that promotes the persistence of polymorphism. While the evolution of assortative mating and its consequences for trait variation and speciation have been extensively studied, the conditions enabling the evolution of disassortative mating are still poorly understood. Mate preferences increase the risk of missing mating opportunities, a cost that can be compensated by a greater fitness of offspring. Heterozygote advantage should therefore promote the evolution of disassortative mating, which maximizes the number of heterozygous offspring. From the analysis of a two-locus diploid model with one locus controlling the mating cue under viability selection and the other locus coding for the level of disassortative preference, we show that heterozygote advantage and negative frequency-dependent viability selection acting at the cue locus promote the evolution of disassortative preferences. We predict conditions of evolution of disassortative mating coherent with selection regimes acting on traits observed in the wild. We also show that disassortative mating generates sexual selection, which disadvantages heterozygotes at the cue locus, limiting the evolution of disassortative preferences. Altogether, our results partially explain why this behavior is rare in natural populations.

Mutation load at a mimicry supergene sheds new light on the evolution of inversion polymorphisms

Chromosomal inversions are ubiquitous in genomes and often coordinate complex phenotypes, such as the covariation of behavior and morphology in many birds, fishes, insects or mammals^1-11. However, why and how inversions become associated with polymorphic traits remains obscure. Here we show that despite a strong selective advantage when they form, inversions accumulate recessive deleterious mutations that generate frequency-dependent selection and promote their maintenance at intermediate frequency. Combining genomics and in vivo fitness analyses in a model butterfly for wing-pattern polymorphism, Heliconius numata, we reveal that three ecologically advantageous inversions have built up a heavy mutational load from the sequential accumulation of deleterious mutations and transposable elements. Inversions associate with sharply reduced viability when homozygous, which prevents them from replacing ancestral chromosome arrangements. Our results suggest that other complex polymorphisms, rather than representing adaptations to competing ecological optima, could evolve because chromosomal rearrangements are intrinsically prone to carrying recessive harmful mutations.