Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U.

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

Cytonuclear interplay in auto- and allopolyploids: a multifaceted perspective from the Festuca-Lolium complex.

Restoring cytonuclear stoichiometry is necessary after whole-genome duplication (WGD) and interspecific/intergeneric hybridization in plants. We investigated this phenomenon in auto- and allopolyploids of the Festuca-Lolium complex providing insights into the mechanisms governing cytonuclear interactions in early polyploid and hybrid generations. Our study examined the main processes potentially involved in restoring the cytonuclear balance after WGD comparing diploids and new and well-established autopolyploids. We uncovered that both the number of chloroplasts and the number of chloroplast genome copies were significantly higher in the newly established autopolyploids and grew further in more established autopolyploids. The increase in the copy number of the chloroplast genome exceeded the rise in the number of chloroplasts and fully compensated for the doubling of the nuclear genome. In addition, changes in nuclear and organelle gene expression were insignificant. Allopolyploid Festuca × Lolium hybrids displayed potential structural conflicts in parental protein variants within the cytonuclear complexes. While biased maternal allele expression has been observed in numerous hybrids, our results suggest that its role in cytonuclear stabilization in the Festuca × Lolium hybrids is limited. This study provides insights into the restoration of the cytonuclear stoichiometry, yet it emphasizes the need for future research to explore post-transcriptional regulation and its impact on cytonuclear gene expression stoichiometry. Our findings may enhance the understanding of polyploid plant evolution, with broader implications for the study of cytonuclear interactions in diverse biological contexts.

Relaxed selection on male mitochondrial genes in DUI bivalves eases the need for mitonuclear coevolution.

Mitonuclear coevolution is an important prerequisite for efficient energy production in eukaryotes. However, many bivalve taxa experience doubly uniparental inheritance (DUI) and have sex-specific mitochondrial (mt) genomes, providing a challenge for mitonuclear coevolution. We examined possible mechanisms to reconcile mitonuclear coevolution with DUI. No nuclear-encoded, sex-specific OXPHOS paralogs were found in the DUI clam Ruditapes philippinarum, refuting OXPHOS paralogy as a solution in this species. It is also unlikely that mt changes causing disruption of nuclear interactions are strongly selected against because sex-specific mt-residues or those under positive selection in M mt genes were not depleted for contacting nuclear-encoded residues. However, M genomes showed consistently higher dN /dS ratios compared to putatively ancestral F genomes in all mt OXPHOS genes and across all DUI species. Further analyses indicated that this was consistently due to relaxed, not positive selection on M vs. F mt OXPHOS genes. Similarly, selection was relaxed on the F genome of DUI species compared to species with strict maternal inheritance. Coupled with recent physiological and molecular evolution studies, we suggest that relaxed selection on M mt function limits the need to maintain mitonuclear interactions in M genomes compared to F genomes. We discuss our findings with regard to OXPHOS function and the origin of DUI.

A wAlbB Wolbachia Transinfection Displays Stable Phenotypic Effects across Divergent Aedes aegypti Mosquito Backgrounds.

Aedes mosquitoes harboring intracellular Wolbachia bacteria are being released in arbovirus and mosquito control programs. With releases taking place around the world, understanding the contribution of host variation to Wolbachia phenotype is crucial. We generated a Wolbachia transinfection (wAlbBQ) in Aedes aegypti and performed backcrossing to introduce the infection into Australian or Malaysian nuclear backgrounds. Whole Wolbachia genome sequencing shows that the wAlbBQ transinfection is nearly identical to the reference wAlbB genome, suggesting few changes since the infection was first introduced to A. aegypti over 15 years ago. However, these sequences were distinct from other available wAlbB genome sequences, highlighting the potential diversity of wAlbB in natural Aedes albopictus populations. Phenotypic comparisons demonstrate the effects of wAlbB infection on egg hatching and nuclear background on fecundity and body size but no interactions between wAlbB infection and nuclear background for any trait. The wAlbB infection was stable at high temperatures and showed perfect maternal transmission and cytoplasmic incompatibility regardless of the host background. Our results demonstrate the stability of wAlbB across host backgrounds and point to its long-term effectiveness for controlling arbovirus transmission and mosquito populations. IMPORTANCEWolbachia bacteria are being used to control the transmission of dengue virus and other arboviruses by mosquitoes. For Wolbachia release programs to be effective globally, Wolbachia infections must be stable across mosquito populations from different locations. In this study, we transferred Wolbachia (strain wAlbB) to Aedes aegypti mosquitoes with an Australian genotype and introduced the infection to Malaysian mosquitoes through backcrossing. We found that the phenotypic effects of Wolbachia are stable across both mosquito backgrounds. We sequenced the genome of wAlbB and found very few genetic changes despite spending over 15 years in a novel mosquito host. Our results suggest that the effects of Wolbachia infections are likely to remain stable across time and host genotype.

Nuclear-cytoplasmic balance: whole genome duplications induce elevated organellar genome copy number.

The plant genome is partitioned across three distinct subcellular compartments: the nucleus, mitochondria, and plastids. Successful coordination of gene expression among these organellar genomes and the nuclear genome is critical for plant function and fitness. Whole genome duplication (WGD) events in the nucleus have played a major role in the diversification of land plants and are expected to perturb the relative copy number (stoichiometry) of nuclear, mitochondrial, and plastid genomes. Thus, elucidating the mechanisms whereby plant cells respond to the cytonuclear stoichiometric imbalance that follows WGDs represents an important yet underexplored question in understanding the evolutionary consequences of genome doubling. We used droplet digital PCR to investigate the relationship between nuclear and organellar genome copy numbers in allopolyploids and their diploid progenitors in both wheat and Arabidopsis. Polyploids exhibit elevated organellar genome copy numbers per cell, largely preserving the cytonuclear stoichiometry observed in diploids despite the change in nuclear genome copy number. To investigate the timescale over which cytonuclear stoichiometry may respond to WGD, we also estimated the organellar genome copy number in Arabidopsis synthetic autopolyploids and in a haploid-induced diploid line. We observed corresponding changes in organellar genome copy number in these laboratory-generated lines, indicating that at least some of the cellular response to cytonuclear stoichiometric imbalance is immediate following WGD. We conclude that increases in organellar genome copy numbers represent a common response to polyploidization, suggesting that maintenance of cytonuclear stoichiometry is an important component in establishing polyploid lineages.

Elucidating Mitochondrial DNA Markers of Ogura-Based CMS Lines in Indian Cauliflowers (Brassica oleracea var. botrytis L.) and Their Floral Abnormalities Due to Diversity in Cytonuclear Interactions.

Mitochondrial markers can be used to differentiate diverse mitotypes as well as cytoplasms in angiosperms. In cauliflower, cultivation of hybrids is pivotal in remunerative agriculture and cytoplasmic male sterile lines constitute an important component of the hybrid breeding. In diversifying the source of male sterility, it is essential to appropriately differentiate among the available male sterile cytoplasms in cauliflower. PCR polymorphism at the key mitochondrial genes associated with male sterility will be instrumental in analyzing, molecular characterization, and development of mitotype-specific markers for differentiation of different cytoplasmic sources. Presence of auto- and alloplasmic cytonuclear combinations result in complex floral abnormalities. In this context, the present investigation highlighted the utility of organelle genome-based markers in distinguishing cytoplasm types in Indian cauliflowers and unveils the epistatic effects of the cytonuclear interactions influencing floral phenotypes. In PCR-based analysis using a set of primers targeted to orf-138, 76 Indian cauliflower lines depicted the presence of Ogura cytoplasm albeit the amplicons generated exhibited polymorphism within the ofr-138 sequence. The polymorphic fragments were found to be spanning over 200-280 bp and 410-470 bp genomic regions of BnTR4 and orf125, respectively. Sequence analysis revealed that such cytoplasmic genetic variations could be attributed to single nucleotide polymorphisms and insertion or deletions of 31/51 nucleotides. The cytoplasmic effects on varying nuclear-genetic backgrounds rendered an array of floral abnormalities like reduction in flower size, fused flowers, splitted style with the exposed ovule, absence of nonfunctional stamens, and petaloid stamens. These floral malformations caused dysplasia of flower structure affecting female fertility with inefficient nectar production. The finding provides an important reference to ameliorate understanding of mechanism of cytonuclear interactions in floral organ development in Brassicas. The study paves the way for unraveling developmental biology of CMS phenotypes in eukaryotic organisms and intergenomic conflict in plant speciation.

Biased Gene Retention in the Face of Introgression Obscures Species Relationships.

Phylogenomic analyses are recovering previously hidden histories of hybridization, revealing the genomic consequences of these events on the architecture of extant genomes. We applied phylogenomic techniques and several complementary statistical tests to show that introgressive hybridization appears to have occurred between close relatives of Arabidopsis, resulting in cytonuclear discordance and impacting our understanding of species relationships in the group. The composition of introgressed and retained genes indicates that selection against incompatible cytonuclear and nuclear-nuclear interactions likely acted during introgression, whereas linkage also contributed to genome composition through the retention of ancient haplotype blocks. We also applied divergence-based tests to determine the species branching order and distinguish donor from recipient lineages. Surprisingly, these analyses suggest that cytonuclear discordance arose via extensive nuclear, rather than cytoplasmic, introgression. If true, this would mean that most of the nuclear genome was displaced during introgression whereas only a small proportion of native alleles were retained.

The toad's warts: Discordance creates bumpy expectations of mitochondrial-nuclear evolution between species.

Discordance between the mitochondrial and nuclear genomes is a prevalent phenomenon in nature, in which the underlying processes responsible are considered to be important in shaping genetic variation in natural populations. Among the evolutionary processes that best explain such genomic mismatches incomplete lineage sorting and introgression are commonly identified, however, many studies are unable to distinguish between these hypotheses, which has become a major challenge in the field. In this issue of Molecular Ecology, Firneno et al. (2020) present an elegant exploration of mitochondrial-nuclear discordance in Mesoamerican toads. Integrating genome-scale and spatial data to test between these hypotheses within an empirical model testing framework, they find strong support that incomplete lineage sorting explains the observed discordance. Their work, along with many previous articles in Molecular Ecology, highlights the commonality of mito-nuclear discordance among species despite the expectations of tightly concerted mitochondrial and nuclear genome evolution. It is increasingly clear that the nuclear genomes of many species are (at least for short periods of evolutionary time) functionally compatible with multiple, divergent mitochondrial haplotypes. As such, we suggest future research not only seeks to understand the processes causing spatial mito-nuclear discordance (e.g. incomplete lineage sorting, introgression), but also explores those that maintain discordance through time and space (e.g. relaxed selection on mito-nuclear interactions, heterozygosity, population demographics). We also discuss the vital role that taxonomy plays in interpreting patterns of mito-nuclear discordance when data-consistent yet differing taxonomies are used, such as treating allopatrically distributed taxa as multiple isolated populations versus multiple micro-endemic species.

Cytonuclear Genetic Incompatibilities in Plant Speciation.

Due to the endosymbiotic origin of organelles, a pattern of coevolution and coadaptation between organellar and nuclear genomes is required for proper cell function. In this review, we focus on the impact of cytonuclear interaction on the reproductive isolation of plant species. We give examples of cases where species exhibit barriers to reproduction which involve plastid-nuclear or mito-nuclear genetic incompatibilities, and describe the evolutionary processes at play. We also discuss potential mechanisms of hybrid fitness recovery such as paternal leakage. Finally, we point out the possible interplay between plant mating systems and cytonuclear coevolution, and its consequence on plant speciation.

Robust Cytonuclear Coordination of Transcription in Nascent Arabidopsis thaliana Autopolyploids.

Polyploidy is hypothesized to cause dosage imbalances between the nucleus and the other genome-containing organelles (mitochondria and plastids), but the evidence for this is limited. We performed RNA-seq on Arabidopsis thaliana diploids and their derived autopolyploids to quantify the degree of inter-genome coordination of transcriptional responses to nuclear whole genome duplication in two different organs (sepals and rosette leaves). We show that nuclear and organellar genomes exhibit highly coordinated responses in both organs. First, organelle genome copy number increased in response to nuclear whole genome duplication (WGD), at least partially compensating for altered nuclear genome dosage. Second, transcriptional output of the different cellular compartments is tuned to maintain diploid-like levels of relative expression among interacting genes. In particular, plastid genes and nuclear genes whose products are plastid-targeted show coordinated down-regulation, such that their expression levels relative to each other remain constant across ploidy levels. Conversely, mitochondrial genes and nuclear genes with mitochondrial targeting show either constant or coordinated up-regulation of expression relative to other nuclear genes. Thus, cytonuclear coordination is robust to changes in nuclear ploidy level, with diploid-like balance in transcript abundances achieved within three generations after nuclear whole genome duplication.

Do angiosperms with highly divergent mitochondrial genomes have altered mitochondrial function?

Angiosperm mitochondrial (mt) genes are generally slow-evolving, but multiple lineages have undergone dramatic accelerations in rates of nucleotide substitution and extreme changes in mt genome structure. While molecular evolution in these lineages has been investigated, very little is known about their mt function. Some studies have suggested altered respiration in individual taxa, although there are several reasons why mt variation might be neutral in others. Here, we develop a new protocol to characterize respiration in isolated plant mitochondria and apply it to species of Silene with mt genomes that are rapidly evolving, highly fragmented, and exceptionally large (~11 Mbp). This protocol, complemented with traditional measures of plant fitness, cytochrome c oxidase activity assays, and fluorescence microscopy, was also used to characterize inter- and intraspecific variation in mt function. Contributions of the individual "classic" OXPHOS complexes, the alternative oxidase, and external NADH dehydrogenases to overall mt respiratory flux were found to be similar to previously studied angiosperms with more typical mt genomes. Some differences in mt function could be explained by inter- and intraspecific variation. This study suggests that Silene species with peculiar mt genomes still show relatively normal mt respiration. This may be due to strong purifying selection on mt variants, coevolutionary responses in the nucleus, or a combination of both. Future experiments should explore such questions using a comparative framework and investigating other lineages with unusual mitogenomes.

Potential limits to the benefits of admixture during biological invasion.

Species introductions often bring together genetically divergent source populations, resulting in genetic admixture. This geographic reshuffling of diversity has the potential to generate favourable new genetic combinations, facilitating the establishment and invasive spread of introduced populations. Observational support for the superior performance of admixed introductions has been mixed, however, and the broad importance of admixture to invasion questioned. Under most underlying mechanisms, admixture's benefits should be expected to increase with greater divergence among and lower genetic diversity within source populations, though these effects have not been quantified in invaders. We experimentally crossed source populations differing in divergence in the invasive plant Centaurea solstitialis. Crosses resulted in many positive (heterotic) interactions, but fitness benefits declined and were ultimately negative at high source divergence, with patterns suggesting cytonuclear epistasis. We explored the literature to assess whether such negative epistatic interactions might be impeding admixture at high source population divergence. Admixed introductions reported for plants came from sources with a wide range of genetic variation, but were disproportionately absent where there was high genetic divergence among native populations. We conclude that while admixture is common in species introductions and often happens under conditions expected to be beneficial to invaders, these conditions may be constrained by predictable negative genetic interactions, potentially explaining conflicting evidence for admixture's benefits to invasion.

The Mitonuclear Dimension of Neanderthal and Denisovan Ancestry in Modern Human Genomes.

Some human populations interbred with Neanderthals and Denisovans, resulting in substantial contributions to modern-human genomes. Therefore, it is now possible to use genomic data to investigate mechanisms that shaped historical gene flow between humans and our closest hominin relatives. More generally, in eukaryotes, mitonuclear interactions have been argued to play a disproportionate role in generating reproductive isolation. There is no evidence of mtDNA introgression into modern human populations, which means that all introgressed nuclear alleles from archaic hominins must function on a modern-human mitochondrial background. Therefore, mitonuclear interactions are also potentially relevant to hominin evolution. We performed a detailed accounting of mtDNA divergence among hominin lineages and used population-genomic data to test the hypothesis that mitonuclear incompatibilities have preferentially restricted the introgression of nuclear genes with mitochondrial functions. We found a small but significant underrepresentation of introgressed Neanderthal alleles at such nuclear loci. Structural analyses of mitochondrial enzyme complexes revealed that these effects are unlikely to be mediated by physically interacting sites in mitochondrial and nuclear gene products. We did not detect any underrepresentation of introgressed Denisovan alleles at mitochondrial-targeted loci, but this may reflect reduced power because locus-specific estimates of Denisovan introgression are more conservative. Overall, we conclude that genes involved in mitochondrial function may have been subject to distinct selection pressures during the history of introgression from archaic hominins but that mitonuclear incompatibilities have had, at most, a small role in shaping genome-wide introgression patterns, perhaps because of limited functional divergence in mtDNA and interacting nuclear genes.

Causes and Consequences of Rapidly Evolving mtDNA in a Plant Lineage.

Understanding mechanisms of coevolution between nuclear and mitochondrial (mt) genomes is a defining challenge in eukaryotic genetics. The angiosperm genus Silene is a natural system to investigate the causes and consequences of mt mutation rate variation because closely related species have highly divergent rates. In Silene species with fast-evolving mtDNA, nuclear genes that encode mitochondrially targeted proteins (N-mt genes) are also fast-evolving. This correlation could indicate positive selection to compensate for mt mutations, but might also result from a recent relaxation of selection. To differentiate between these interpretations, we used phylogenetic and population-genetic methods to test for positive and relaxed selection in three classes of N-mt genes (oxidative phosphorylation genes, ribosomal genes, and "RRR" genes involved in mtDNA recombination, replication, and repair). In all three classes, we found that species with fast-evolving mtDNA had: 1) elevated dN/dS, 2) an excess of nonsynonymous divergence relative to levels of intraspecific polymorphism, which is a signature of positive selection, and 3) no clear signals of relaxed selection. "Control" genes exhibited comparatively few signs of positive selection. These results suggest that high mt mutation rates can create selection on N-mt genes and that relaxed selection is an unlikely cause of recent accelerations in the evolution of N-mt genes. Because mt-RRR genes were found to be under positive selection, it is unlikely that elevated mt mutation rates in Silene were caused by inactivation of these mt-RRR genes. Therefore, the causes of extreme increases in angiosperm mt mutation rates remain uncertain.

The on-again, off-again relationship between mitochondrial genomes and species boundaries.

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co-evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.

Intraspecific genetic variation in hosts affects regulation of obligate heritable symbionts.

Symbiotic relationships promote biological diversification by unlocking new ecological niches. Over evolutionary time, hosts and symbionts often enter intimate and permanent relationships, which must be maintained and regulated for both lineages to persist. Many insect species harbor obligate, heritable symbiotic bacteria that provision essential nutrients and enable hosts to exploit niches that would otherwise be unavailable. Hosts must regulate symbiont population sizes, but optimal regulation may be affected by the need to respond to the ongoing evolution of symbionts, which experience high levels of genetic drift and potential selection for selfish traits. We address the extent of intraspecific variation in the regulation of a mutually obligate symbiosis, between the pea aphid (Acyrthosiphon pisum) and its maternally transmitted symbiont, Buchnera aphidicola Using experimental crosses to identify effects of host genotypes, we measured symbiont titer, as the ratio of genomic copy numbers of symbiont and host, as well as developmental time and fecundity of hosts. We find a large (>10-fold) range in symbiont titer among genetically distinct aphid lines harboring the same Buchnera haplotype. Aphid clones also vary in fitness, measured as developmental time and fecundity, and genetically based variation in titer is correlated with host fitness, with higher titers corresponding to lower reproductive rates of hosts. Our work shows that obligate symbiosis is not static but instead is subject to short-term evolutionary dynamics, potentially reflecting coevolutionary interactions between host and symbiont.

Positive Selection in Rapidly Evolving Plastid-Nuclear Enzyme Complexes.

Rates of sequence evolution in plastid genomes are generally low, but numerous angiosperm lineages exhibit accelerated evolutionary rates in similar subsets of plastid genes. These genes include clpP1 and accD, which encode components of the caseinolytic protease (CLP) and acetyl-coA carboxylase (ACCase) complexes, respectively. Whether these extreme and repeated accelerations in rates of plastid genome evolution result from adaptive change in proteins (i.e., positive selection) or simply a loss of functional constraint (i.e., relaxed purifying selection) is a source of ongoing controversy. To address this, we have taken advantage of the multiple independent accelerations that have occurred within the genus Silene (Caryophyllaceae) by examining phylogenetic and population genetic variation in the nuclear genes that encode subunits of the CLP and ACCase complexes. We found that, in species with accelerated plastid genome evolution, the nuclear-encoded subunits in the CLP and ACCase complexes are also evolving rapidly, especially those involved in direct physical interactions with plastid-encoded proteins. A massive excess of nonsynonymous substitutions between species relative to levels of intraspecific polymorphism indicated a history of strong positive selection (particularly in CLP genes). Interestingly, however, some species are likely undergoing loss of the native (heteromeric) plastid ACCase and putative functional replacement by a duplicated cytosolic (homomeric) ACCase. Overall, the patterns of molecular evolution in these plastid-nuclear complexes are unusual for anciently conserved enzymes. They instead resemble cases of antagonistic coevolution between pathogens and host immune genes. We discuss a possible role of plastid-nuclear conflict as a novel cause of accelerated evolution.

"Darwin's corollary" and cytoplasmic incompatibility induced by Cardinium may contribute to speciation in Encarsia wasps (Hymenoptera: Aphelinidae).

The potential importance of cytoplasmic incompatibility (CI)-inducing bacterial symbionts in speciation of their arthropod hosts has been debated. Theoretical advances have led to a consensus that a role is plausible when CI is combined with other isolating barriers. However, the insect model systems Nasonia and Drosophila are the only two experimental examples documented. Here, we analyzed the components of reproductive isolation between the parasitoid wasp Encarsia suzannae, which is infected by the CI-inducing symbiont Cardinium, and its uninfected sibling species Encarsia gennaroi. Laboratory crosses demonstrated that: (1) sexual isolation is incomplete; (2) hybrid offspring production is greatly reduced in the interspecific CI cross; (3) viable hybrids may be produced by curing E. suzannae males of Cardinium with antibiotics; (4) hybrid offspring production in the reciprocal cross is greatly reduced by hybrid inviability due to genetic incompatibilities; (5) hybrid sterility is nearly complete in both directions at the F1 stage. Thus, asymmetrical hybrid incompatibilities and CI act as complementary isolating mechanisms. We propose a new model for contributions of CI symbionts to speciation, with CI reducing gene flow between species in one direction, and in the other, a symbiont sweep resulting in accelerated mtDNA evolution, negative cytonuclear interactions, and hybrid incompatibilities.

The Roles of Mutation, Selection, and Expression in Determining Relative Rates of Evolution in Mitochondrial versus Nuclear Genomes.

Eukaryotes rely on proteins encoded by the nuclear and mitochondrial (mt) genomes, which interact within multisubunit complexes such as oxidative-phosphorylation enzymes. Although selection is thought to be less efficient on the asexual mt genome, in bilaterian animals the ratio of nonsynonymous to synonymous substitutions (ω) is lower in mt- compared with nuclear-encoded OXPHOS subunits, suggesting stronger effects of purifying selection in the mt genome. Because high levels of gene expression constrain protein sequence evolution, one proposed resolution to this paradox is that mt genes are expressed more highly than nuclear genes. To test this hypothesis, we investigated expression and sequence evolution of mt and nuclear genes from 84 diverse eukaryotes that vary in mt gene content and mutation rate. We found that the relationship between mt and nuclear ω values varied dramatically across eukaryotes. In contrast, transcript abundance is consistently higher for mt genes than nuclear genes, regardless of which genes happen to be in the mt genome. Consequently, expression levels cannot be responsible for the differences in ω Rather, 84% of the variance in the ratio of ω values between mt and nuclear genes could be explained by differences in mutation rate between the two genomes. We relate these findings to the hypothesis that high rates of mt mutation select for compensatory changes in the nuclear genome. We also propose an explanation for why mt transcripts consistently outnumber their nuclear counterparts, with implications for mitonuclear protein imbalance and aging.