Comparative analysis of multiple genomes in a phylogenetic framework dramatically improves the precision and sensitivity of evolutionary inference, producing more robust results than single-genome analyses can provide. The genomes of 12 Drosophila species, ten of which are presented here for the first time (sechellia, simulans, yakuba, erecta, ananassae, persimilis, willistoni, mojavensis, virilis and grimshawi), illustrate how rates and patterns of sequence divergence across taxa can illuminate evolutionary processes on a genomic scale. These genome sequences augment the formidable genetic tools that have made Drosophila melanogaster a pre-eminent model for animal genetics, and will further catalyse fundamental research on mechanisms of development, cell biology, genetics, disease, neurobiology, behaviour, physiology and evolution. Despite remarkable similarities among these Drosophila species, we identified many putatively non-neutral changes in protein-coding genes, non-coding RNA genes, and cis-regulatory regions. These may prove to underlie differences in the ecology and behaviour of these diverse species.

Mitochondrial transcription, translation, and respiration require interactions between genes encoded in two distinct genomes, generating the potential for mutations in nuclear and mitochondrial genomes to interact epistatically and cause incompatibilities that decrease fitness. Mitochondrial-nuclear epistasis for fitness has been documented within and between populations and species of diverse taxa, but rarely has the genetic or mechanistic basis of these mitochondrial–nuclear interactions been elucidated, limiting our understanding of which genes harbor variants causing mitochondrial–nuclear disruption and of the pathways and processes that are impacted by mitochondrial–nuclear coevolution. Here we identify an amino acid polymorphism in the Drosophila melanogaster nuclear-encoded mitochondrial tyrosyl–tRNA synthetase that interacts epistatically with a polymorphism in the D. simulans mitochondrial-encoded tRNATyr to significantly delay development, compromise bristle formation, and decrease fecundity. The incompatible genotype specifically decreases the activities of oxidative phosphorylation complexes I, III, and IV that contain mitochondrial-encoded subunits. Combined with the identity of the interacting alleles, this pattern indicates that mitochondrial protein translation is affected by this interaction. Our findings suggest that interactions between mitochondrial tRNAs and their nuclear-encoded tRNA synthetases may be targets of compensatory molecular evolution. Human mitochondrial diseases are often genetically complex and variable in penetrance, and the mitochondrial–nuclear interaction we document provides a plausible mechanism to explain this complexity.

Efficient mitochondrial function requires physical interactions between the proteins encoded by the mitochondrial and nuclear genomes. Co-evolution between these genomes may result in the accumulation of incompatibilities between divergent lineages. We test whether mitochondrialnuclear incompatibilities have accumulated within the Drosophila melanogaster species subgroup by combining divergent mitochondrial and nuclear lineages and quantifying the effects on relative fitness. Precise placement of nine mtDNAs from D. melanogaster, D. simulans and D. mauritiana into two D. melanogaster nuclear genetic backgrounds reveals significant mitochondrial-nuclear epistasis affecting fitness in females. Combining the mitochondrial genomes with three different D. melanogaster X chromosomes reveals significant epistasis for male fitness between X-linked and mitochondrial variation. However, we find no evidence that the more than 500 fixed differences between the mitochondrial genomes of D. melanogaster and the D. simulans species complex are incompatible with the D. melanogaster nuclear genome. Rather, the interactions of largest effect occur between mitochondrial and nuclear polymorphisms that segregate within species of the D. melanogaster species subgroup. We propose that a low mitochondrial substitution rate, resulting from a low mutation rate and/or efficient purifying selection, precludes the accumulation of mitochondrial-nuclear incompatibilities among these Drosophila species.