Lethal epistasis maintains strong linkage disequilibrium between unlinked supergenes

Linkage disequilibrium (LD) between adaptive gene combinations is typically maintained through physical linkage and suppressed recombination. Although epistatic interactions can maintain LD between unlinked loci, this mechanism is rarely documented in nature. Supergenes, regions of suppressed recombination containing tightly linked loci, typically exemplify the first mechanism, in which genomic rearrangements lock together coadapted alleles. Here, we demonstrate a rare case of epistasis-driven LD between two supergenes on different chromosomes in the European ant Formica cinerea : one controlling colony queen number (chromosome 3) and another determining sexual body size (chromosome 9). We show that these supergenes assort independently according to Mendelian expectations during meiosis, yet exhibit high LD between the multi-queen haplotype P 2 and the miniaturizing haplotype 9r, with small queens and males occurring exclusively in multi-queen colonies. Mismatched genotype combinations (P 2 without 9r and vice versa) are severely underrepresented among all adult castes, with mismatched males experiencing complete mortality. We show that this pattern cannot be explained by meiotic drive or maternal-effect killing, indicating strong postzygotic epistatic selection maintaining the observed LD. The fitness costs of these epistatic interactions are substantial but critically depend on mating combinations: heterozygous small queens achieve 1.5-2.5 times higher fitness (based on offspring genotype viability) when mated to small males (P 2 -9r) compared to large males. Our findings provide empirical evidence for epistatic interactions between unlinked supergenes maintaining LD in the absence of physical linkage.