Dissecting contributions of directional and balancing selection to trajectories of mitochondrial haplotype evolution in Drosophila melanogaster.

Emerging evidence suggests mtDNA haplotypes contribute to fitness variation and local adaptation, with directional thermal selection and negative frequency-dependent selection shaping haplotype diversity. However, their interplay remains unexplored. We conducted experimental evolution using Drosophila melanogaster populations from opposite ends of an Australian latitudinal cline (Melbourne and Townsville), exposing them to contrasting temperatures (17°C versus 27°C) and varying starting frequencies of two mtDNA haplotypes (A1 and B1) that occur at appreciable frequencies in these populations. We paired this with population genetic simulations to estimate selection and its influence on haplotype trajectories. Haplotype frequencies were influenced by interactions involving temperature, starting frequency, and nuclear genomic background (Melbourne, Townsville, or admixed). Although prior work predicted A1 should be favoured at the warmer temperature and B1 at the cooler temperature, A1 was generally favoured across both temperatures. Simulations supported directional selection in populations evolving at 17°C in the Melbourne background; otherwise dynamics were best explained by balancing selection shaped by negative frequency-dependent fitness effects. Patterns also varied across nuclear backgrounds, suggestive of mito-nuclear epistasis. These findings challenge a simple thermal adaptation model of mtDNA dynamics, suggesting that mtDNA evolution is shaped by interacting effects of temperature, frequency-dependence, nuclear background and experimental environment.