The decoupled evolution of the organellar genomes of Silene nutans leads to distinct roles in the speciation process

There is growing evidence that cytonuclear incompatibilities (i.e. disruption of cytonuclear coadaptation) might contribute to the speciation process. In a former study, we described the possible involvement of plastid-nuclear incompatibilities in the reproductive isolation between four lineages of Silene nutans (Caryophyllaceae). Because organellar genomes are usually cotransmitted, we assessed whether the mitochondrial genome could also be involved in the speciation process, knowing that the gynodioecious breeding system of S. nutans is expected to impact the evolutionary dynamics of this genome. Using hybrid capture and high-throughput DNA sequencing, we analyzed diversity patterns in the genic content of the organellar genomes in the four S. nutans lineages. Contrary to the plastid genome, which exhibited a large number of fixed substitutions between lineages, extensive sharing of polymorphisms between lineages was found in the mitochondrial genome. In addition, numerous recombination-like events were detected in the mitochondrial genome, loosening the linkage disequilibrium between the organellar genomes and leading to decoupled evolution. These results suggest that gynodioecy shaped mitochondrial diversity through balancing selection, maintaining ancestral polymorphism and, thus, limiting the involvement of the mitochondrial genome in evolution of hybrid inviability between S. nutans lineages.