Complex mitonuclear interactions and metabolic costs of mating in male seed beetles

The lack of evolutionary response to selection on mitochondrial genes through males predicts the evolution of nuclear genetic influence on male‐specific mitochondrial function, for example by gene duplication and evolution of sex‐specific expression of paralogs involved in metabolic pathways. Intergenomic epistasis may therefore be a prevalent feature of the genetic architecture of male‐specific organismal function. Here, we assess the role of mitonuclear genetic variation for male metabolic phenotypes [metabolic rate and respiratory quotient (RQ)] associated with ejaculate renewal, in the seed beetle Callosobruchus maculatus, by assaying lines with crossed combinations of distinct mitochondrial haplotypes and nuclear lineages. We found a significant increase in metabolic rate following mating relative to virgin males. Moreover, processes associated with ejaculate renewal showed variation in metabolic rate that was affected by mitonuclear interactions. Mitochondrial haplotype influenced mating‐related changes in RQ, but this pattern varied over time. Mitonuclear genotype and the energy spent during ejaculate production affected the weight of the ejaculate, but the strength of this effect varied across mitochondrial haplotypes showing that the genetic architecture of male‐specific reproductive function is complex. Our findings unveil hitherto underappreciated metabolic costs of mating and ejaculate renewal, and provide the first empirical demonstration of mitonuclear epistasis on male reproductive metabolic processes.     

Positive and purifying selection in mitochondrial genomes of a bird with mitonuclear discordance

Diversifying selection on metabolic pathways can reduce intraspecific gene flow and promote population divergence. An opportunity to explore this arises from mitonuclear discordance observed in an Australian bird Eopsaltria australis. Across >1500 km, nuclear differentiation is low and latitudinally structured by isolation by distance, whereas two highly divergent, parapatric mitochondrial lineages (>6.6% in ND2) show a discordant longitudinal geographic pattern and experience different climates. Vicariance, incomplete lineage sorting and sex‐biased dispersal were shown earlier to be unlikely drivers of the mitonuclear discordance; instead, natural selection on a female‐linked trait was the preferred hypothesis. Accordingly, here we tested for signals of positive, divergent selection on mitochondrial genes in E. australis. We used codon models and physicochemical profiles of amino acid replacements to analyse complete mitochondrial genomes of the two mitochondrial lineages in E. australis, its sister species Eopsaltria griseogularis, and outgroups. We found evidence of positive selection on at least five amino acids, encoded by genes of two oxidative phosphorylation pathway complexes NADH dehydrogenase (ND4 and ND4L) and cytochrome bc₁ (cyt‐b) against a background of widespread purifying selection on all mitochondrial genes. Three of these amino acid replacements were fixed in ND4 of the geographically most widespread E. australis lineage. The other two replacements were fixed in ND4L and cyt‐b of the geographically more restricted E. australis lineage. We discuss whether this selection may reflect local environmental adaptation, a by‐product of other selective processes, or genetic incompatibilities, and propose how these hypotheses can be tested in future.     

Direct and indirect genetic effects of sex-specific mitonuclear epistasis on reproductive ageing

Mitochondria are involved in ageing and their function requires coordinated action of both mitochondrial and nuclear genes. Epistasis between the two genomes can influence lifespan but whether this also holds for reproductive senescence is unclear. Maternal inheritance of mitochondria predicts sex differences in the efficacy of selection on mitonuclear genotypes that should result in differences between females and males in mitochondrial genetic effects. Mitonuclear genotype of a focal individual may also indirectly affect trait expression in the mating partner. We tested these predictions in the seed beetle Callosobruchus maculatus, using introgression lines harbouring distinct mitonuclear genotypes. Our results reveal both direct and indirect sex-specific effects of mitonuclear epistasis on reproductive ageing. Females harbouring coadapted mitonuclear genotypes showed higher lifetime fecundity due to slower senescence relative to novel mitonuclear combinations. We found no evidence for mitonuclear coadaptation in males. Mitonuclear epistasis not only affected age-specific ejaculate weight, but also influenced male age-dependent indirect effects on traits expressed by their female partners (fecundity, egg size, longevity). These results demonstrate important consequences of sex-specific mitonuclear epistasis for both mating partners, consistent with a role for mitonuclear genetic constraints upon sex-specific adaptive evolution.     

Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population

Hybridization between divergent lineages generates new allelic combinations. One mechanism that can hinder the formation of hybrid populations is mitonuclear incompatibility, that is, dysfunctional interactions between proteins encoded in the nuclear and mitochondrial genomes (mitogenomes) of diverged lineages. Theoretically, selective pressure due to mitonuclear incompatibility can affect genotypes in a hybrid population in which nuclear genomes and mitogenomes from divergent lineages admix. To directly and thoroughly observe this key process, we de novo sequenced the 747‐Mb genome of the coastal goby, Chaenogobius annularis, and investigated its integrative genomic phylogeographics using RNA‐sequencing, RAD‐sequencing, genome resequencing, whole mitogenome sequencing, amplicon sequencing, and small RNA‐sequencing. Chaenogobius annularis populations have been geographically separated into Pacific Ocean (PO) and Sea of Japan (SJ) lineages by past isolation events around the Japanese archipelago. Despite the divergence history and potential mitonuclear incompatibility between these lineages, the mitogenomes of the PO and SJ lineages have coexisted for generations in a hybrid population on the Sanriku Coast. Our analyses revealed accumulation of nonsynonymous substitutions in the PO‐lineage mitogenomes, including two convergent substitutions, as well as signals of mitochondrial lineage‐specific selection on mitochondria‐related nuclear genes. Finally, our data implied that a microRNA gene was involved in resolving mitonuclear incompatibility. Our integrative genomic phylogeographic approach revealed that mitonuclear incompatibility can affect genome evolution in a natural hybrid population.     

Conservative and compensatory evolution in oxidative phosphorylation complexes of angiosperms with highly divergent rates of mitochondrial genome evolution

Interactions between nuclear and mitochondrial gene products are critical for eukaryotic cell function. Nuclear genes encoding mitochondrial‐targeted proteins (N‐mt genes) experience elevated rates of evolution, which has often been interpreted as evidence of nuclear compensation in response to elevated mitochondrial mutation rates. However, N‐mt genes may be under relaxed functional constraints, which could also explain observed increases in their evolutionary rate. To disentangle these hypotheses, we examined patterns of sequence and structural evolution in nuclear‐ and mitochondrial‐encoded oxidative phosphorylation proteins from species in the angiosperm genus Silene with vastly different mitochondrial mutation rates. We found correlated increases in N‐mt gene evolution in species with fast‐evolving mitochondrial DNA. Structural modeling revealed an overrepresentation of N‐mt substitutions at positions that directly contact mutated residues in mitochondrial‐encoded proteins, despite overall patterns of conservative structural evolution. These findings support the hypothesis that selection for compensatory changes in response to mitochondrial mutations contributes to the elevated rate of evolution in N‐mt genes. We discuss these results in light of theories implicating mitochondrial mutation rates and mitonuclear coevolution as drivers of speciation and suggest comparative and experimental approaches that could take advantage of heterogeneity in rates of mtDNA evolution across eukaryotes to evaluate such theories.     

Comparative Genomics of Mitochondrial DNA in Drosophila simulans

The current study compares the nucleotide variation among 22 complete mitochondrial genomes of the three distinct Drosophila simulans haplotypes with intron 1 of the alcohol dehydrogenase-related locus. This is the first study to investigate the sequence variation of multiple complete mitochondrial genomes within distinct mitochondrial haplotypes of a single species. Patterns of variation suggest distinct forces are influencing the evolution of mitochondrial DNA (mtDNA) and autosomal DNA in D. simulans. First, there is little variation within each mtDNA haplotype but strong differentiation among them. In contrast, there is no support for differentiation of the mitochondrial haplotypes at the autosomal locus. Second, there is a significant deficiency of mitochondrial variation in each haplotype relative to the autosomal locus. Third, the ratio of nonsynonymous to synonymous substitutions is not equal in all branches of the well-resolved phylogeny. There is an excess of nonsynonymous substitutions relative to synonymous substitutions within each D. simulans haplotype. This result is similar to that previously observed within the mtDNA of distinct species. A single evolutionary force may be causally linked to the observed patterns of mtDNA variation—a rickettsia-like microorganism, Wolbachia pipientis, which is known to directly influence mitochondrial evolution but have a less direct influence on autosomal loci. Received: 16 September 1999 / Accepted: 14 March 2000     

Evolutionary neutrality of mtDNA introgression: evidence from complete mitogenome analysis in roe deer

Introgressive hybridization offers a unique platform for studying the molecular basis of natural selection acting on mitogenomes. Most of the mtDNA protein‐coding genes are extremely conserved; however, some of the observed variations have potentially adaptive significance. Here, we evaluated whether the evolution of mtDNA in closely related roe deer species affected by widespread mtDNA introgression is neutral or adaptive. We characterized and compared 16 complete mitogenomes of European (Capreolus capreolus) and Siberian (C. pygargus) roe deer, including four of Siberian origin introgressed into European species. The average sequence divergence of species‐specific lineages was estimated at 2.8% and varied across gene classes. Only 21 of 315 fixed differences identified in protein‐coding genes represented nonsynonymous changes. Only three of them were determined to have arisen in the C. pygargus lineage since the time to the most recent common ancestor (TMRCA) of both Capreolus species, reflecting a decelerated evolutionary ratio. The almost four‐fold higher d_N/d_S ratio described for the European roe deer lineage is constrained by overall purifying selection, especially pronounced in the ND4 and ND5 genes. We suggest that the highly divergent C. capreolus lineage could have maintained a capability for genomic incorporation of the well‐preserved and almost ancestral type of mtDNA present in C. pygargus. Our analyses did not indicate any signs of positive selection for Siberian roe deer mtDNA, suggesting that the present widespread introgression is evolutionarily neutral.     

Adaptive evolution of the <Emphasis Type="Italic">Homo</Emphasis> mitochondrial genome

Adaptive evolution of 12 protein-coding mitochondrial genes in members of genus Homo (Denisova hominin (H. sp. Altai), Neanderthals (H. neanderthalensis) and modern humans (H. sapiens)) has been evaluated by assessing the pattern of changes in the physicochemical properties of amino acid replacements during primate evolution. It has been found that molecular adaptation (positive destabilizing selection) in Homo becomes apparent in the form of 12 radical amino acid replacements accompanied with statistically significant (P < 0.001) changes of physicochemical properties that probably had functional consequences. These replacements occurred at the stage of a common ancestor of Homo (in CO2 and CytB genes) as well as with the appearance of the common ancestor of Neanderthals and modern humans (in CO1 and ND5 genes). Radical amino acid replacements were mainly revealed in the cytochrome c oxidase complex IV and cytochrome bc1 complex III, thus coinciding with the general trend of increasing nonsynonymous changes in mtDNA genes coding subunits of complexes’ III and IV proteins in anthropoid primates.     

Differential fitness of mitochondrial DNA in perturbation cage studies correlates with global abundance and population history in Drosophila simulans

Mitochondria are often referred to as the powerhouse of the cell. However, research linking intraspecific differences in organismal fitness with genotypic mitochondrial DNA (mtDNA) variation has been hampered by the lack of variation in experimentally tractable species. This study examines whether fly lines harbouring three distinct Drosophila simulans mtDNA types (siI, -II and -III) exhibit differential fitness in laboratory perturbation cages. Comparison of the pre-perturbation and post-perturbation data shows that both the mtDNA and mitonuclear interactions have a significant and repeatable effect on the frequency of flies with specific genotypes in population cages (siII > -III > -I) and that coadapted mitonuclear interactions are greatest in the siI type. The rank order of mtDNA frequency correlates with the observed worldwide distribution of the haplogroups while mitonuclear interactions are most significant in the siI haplogroup that is likely to have been subject to repeated population bottlenecks. One possible explanation for the maintenance of the least fit siI haplogroup on Pacific islands is that it is protected from extinction by Wolbachia infection.     

Mitonuclear Coevolution,but not Nuclear Compensation,Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.     

Comparative genomics and the evolution of human mitochondrial DNA: assessing the effects of selection

This article provides evidence that selection has been a significant force during the evolution of the human mitochondrial genome. Both gene-by-gene and whole-genome approaches were used here to assess selection in the 560 mitochondrial DNA (mtDNA) coding-region sequences that were used previously for reduced-median-network analysis. The results of the present analyses were complex, in that the action of selection was not indicated by all tests, but this is not surprising, in view of the characteristics and limitations of the different analytical methods. Despite these limitations, there is evidence for both gene-specific and lineage-specific variation in selection. Whole-genome sliding-window approaches indicated a lack of selection in large-scale segments of the coding region. In other tests, we analyzed the ratio of nonsynonymous-to-synonymous substitutions in the 13 protein-encoding mtDNA genes. The most straightforward interpretation of those results is that negative selection has acted on the mtDNA during evolution. Single-gene analyses indicated significant departures from neutrality in the CO1, ND4, and ND6 genes, although the data also suggested the possible operation of positive selection on the AT6 gene. Finally, our results and those of other investigators do not support a simple model in which climatic adaptation has been a major force during human mtDNA evolution.     

Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans,but not in chimpanzees.

To test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution, nucleotide sequences were determined for the 1041 bp of the NADH dehydrogenase subunit 2 (ND2) gene in 20 geographically diverse humans and 20 common chimpanzees. Contingency tests of neutrality were performed using four mutational categories for the ND2 molecule: synonymous and nonsynonymous mutations in the transmembrane regions, and synonymous and nonsynonymous mutations in the surface regions. The following three topological mutational categories were also used: intraspecific tips, intraspecific interiors, and interspecific fixed differences. The analyses reveal a significantly greater number of nonsynonymous polymorphisms within human transmembrane regions than expected based on interspecific comparisons, and they are inconsistent with a neutral equilibrium model. This pattern of excess nonsynonymous polymorphism is not seen within chimpanzees. Statistical tests of neutrality, such as TAJIMA''s D test, and the D and F tests proposed by FU and LI, indicate an excess of low frequency polymorphisms in the human data, but not in the chimpanzee data. This is consistent with recent directional selection, a population bottleneck or background selection of slightly deleterious mutations in human mtDNA samples. The analyses further support the idea that mitochondrial genome evolution is governed by selective forces that have the potential to affect its use as a "neutral" marker in evolutionary and population genetic studies.     

Within‐population genetic effects of mtDNA on metabolic rate in Drosophila subobscura

A growing body of research supports the view that within‐species sequence variation in the mitochondrial genome (mtDNA) is functional, in the sense that it has important phenotypic effects. However, most of this empirical foundation is based on comparisons across populations, and few studies have addressed the functional significance of mtDNA polymorphism within populations. Here, using mitonuclear introgression lines, we assess differences in whole‐organism metabolic rate of adult Drosophila subobscura fruit flies carrying either of three different sympatric mtDNA haplotypes. We document sizeable, up to 20%, differences in metabolic rate across these mtDNA haplotypes. Further, these mtDNA effects are to some extent sex specific. We found no significant nuclear or mitonuclear genetic effects on metabolic rate, consistent with a low degree of linkage disequilibrium between mitochondrial and nuclear genes within populations. The fact that mtDNA haplotype variation within a natural population affects metabolic rate, which is a key physiological trait with important effects on life‐history traits, adds weight to the emergent view that mtDNA haplotype variation is under natural selection and it revitalizes the question as to what processes act to maintain functional mtDNA polymorphism within populations.     

Mitochondrial rate variation among lineages of passerine birds

The order Passeriformes comprises the majority of extant avian species. Analyses of molecular data have provided important insights into the evolution of this diverse order. However, molecular estimates of the evolutionary and demographic timescales of passerine species have been hindered by a lack of reliable calibrations. This has led to a reliance on the application of standard substitution rates to mitochondrial DNA data, particularly rates estimated from analyses of the gene encoding cytochrome b (CYTB). To investigate patterns of rate variation across passerine lineages, we used a Bayesian phylogenetic approach to analyse the protein‐coding genes of 183 mitochondrial genomes. We found that the most commonly used mitochondrial marker, CYTB, has low variation in rates across passerine lineages. This lends support to its widespread use as a molecular clock in birds. However, we also found that the patterns of among‐lineage rate variation in CYTB are only weakly related to the evolutionary rate of the mitochondrial genome as a whole. Our analyses confirmed the presence of mutational saturation at third codon positions across the protein‐coding genes of the mitochondrial genome, reinforcing the view that these sites should be excluded in studies of deep passerine relationships. The results of our analyses have provided information that will be useful for molecular‐clock studies of passerine evolution.     

Sex‐specific mitonuclear epistasis and the evolution of mitochondrial bioenergetics,ageing, and life history in seed beetles

The role of mitochondrial DNA for the evolution of life‐history traits remains debated. We examined mitonuclear effects on the activity of the multisubunit complex of the electron transport chain (ETC) involved in oxidative phosphorylation (OXPHOS) across lines of the seed beetle Acanthoscelides obtectus selected for a short (E) or a long (L) life for more than >160 generations. We constructed and phenotyped mitonuclear introgression lines, which allowed us to assess the independent effects of the evolutionary history of the nuclear and the mitochondrial genome. The nuclear genome was responsible for the largest share of divergence seen in ageing. However, the mitochondrial genome also had sizeable effects, which were sex‐specific and expressed primarily as epistatic interactions with the nuclear genome. The effects of mitonuclear disruption were largely consistent with mitonuclear coadaptation. Variation in ETC activity explained a large proportion of variance in ageing and life‐history traits and this multivariate relationship differed somewhat between the sexes. In conclusion, mitonuclear epistasis has played an important role in the laboratory evolution of ETC complex activity, ageing, and life histories and these are closely associated. The mitonuclear architecture of evolved differences in life‐history traits and mitochondrial bioenergetics was sex‐specific.     

NATURALLY OCCURRING MITOCHONDRIAL DNA HAPLOTYPES EXHIBIT METABOLIC DIFFERENCES: INSIGHT INTO FUNCTIONAL PROPERTIES OF MITOCHONDRIA

Linking the mitochondrial genotype and the organismal phenotype is of paramount importance in evolution of mitochondria. In this study, we determined the differences in catalytic properties of mitochondria dictated by divergences in the siII and siIII haplogroups of Drosophila simulans using introgressions of siII mtDNA type into the siIII nuclear background. We used a novel in situ method (permeabilized fibers) that allowed us to accurately measure the consumption of oxygen by mitochondria in constructed siII‐introgressed flies and in siIII‐control flies. Our results showed that the catalytic capacity of the electron transport system is not impaired by introgressions, suggesting that the functional properties of mitochondria are tightly related to the mtDNA haplogroup and not to the nuclear DNA or to the mito‐nuclear interactions. This is the first study, to our knowledge, that demonstrates a naturally occurring haplogroup can confer specific functional differences in aspects of mitochondrial metabolism. This study illustrates the importance of mtDNA changes on organelle evolution and highlights the potential bioenergetic and metabolic impacts that divergent mitochondrial haplogroups may have upon a wide variety of species including humans.