Modeling Mito-nuclear Compatibility and Its Role in Species Identification

Mitochondrial genetic material (mtDNA) is widely used for phylogenetic reconstruction and as a barcode for species identification. The utility of mtDNA in these contexts derives from its particular molecular properties, including its high evolutionary rate, uniparental inheritance, and small size. ButmtDNA may also play a fundamental role in speciation-as suggested by recent observations of coevolution with the nuclear DNA, along with the fact that respiration depends on coordination of genes from both sources. Here, we study how mito-nuclear interactions affect the accuracy of species identification by mtDNA, as well as the speciation process itself. We simulate the evolution of a population of individuals who carry a recombining nuclear genome and a mitochondrial genome inherited maternally. We compare a null model fitness landscape that lacks any mito-nuclear interaction against a scenario in which interactions influence fitness. Fitness is assigned to individuals according to their mito-nuclear compatibility, which drives the coevolution of the nuclear and mitochondrial genomes. Depending on the model parameters, the population breaks into distinct species and the model output then allows us to analyze the accuracy of mtDNA barcode for species identification. Remarkably, we find that species identification by mtDNA is equally accurate in the presence or absence of mito-nuclear coupling and that the success of the DNA barcode derives mainly from population geographical isolation during speciation. Nevertheless, selection imposed by mito-nuclear compatibility influences the diversification process and leaves signatures in the genetic content and spatial distribution of the populations, in three ways. First, speciation is delayed and the resulting phylogenetic trees are more balanced. Second, clades in the resulting phylogenetic tree correlate more strongly with the spatial distribution of species and clusters of more similar mtDNA's. Third, there is a substantial increase in the intraspecies mtDNA similarity, decreasing the number of alleles substitutions per locus and promoting the conservation of genetic information. We compare the evolutionary patterns observed in our model to empirical data from copepods (Tigriopus californicus). We find good qualitative agreement in the geographic patterns and the topology of the phylogenetic tree, provided the model includes selection based on mito-nuclear interactions. These results highlight the role of mito-nuclear compatibility in the speciation process and its reconstruction from genetic data. (AU)