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.     

Mitonuclear coevolution as the genesis of speciation and the mitochondrial DNA barcode gap

Mitochondrial genes are widely used in taxonomy and systematics because high mutation rates lead to rapid sequence divergence and because such changes have long been assumed to be neutral with respect to function. In particular, the nucleotide sequence of the mitochondrial gene cytochrome c oxidase subunit 1 has been established as a highly effective DNA barcode for diagnosing the species boundaries of animals. Rarely considered in discussions of mitochondrial evolution in the context of systematics, speciation, or DNA barcodes, however, is the genomic architecture of the eukaryotes: Mitochondrial and nuclear genes must function in tight coordination to produce the complexes of the electron transport chain and enable cellular respiration. Coadaptation of these interacting gene products is essential for organism function. I extend the hypothesis that mitonuclear interactions are integral to the process of speciation. To maintain mitonuclear coadaptation, nuclear genes, which code for proteins in mitochondria that cofunction with the products of mitochondrial genes, must coevolve with rapidly changing mitochondrial genes. Mitonuclear coevolution in isolated populations leads to speciation because population‐specific mitonuclear coadaptations create between‐population mitonuclear incompatibilities and hence barriers to gene flow between populations. In addition, selection for adaptive divergence of products of mitochondrial genes, particularly in response to climate or altitude, can lead to rapid fixation of novel mitochondrial genotypes between populations and consequently to disruption in gene flow between populations as the initiating step in animal speciation. By this model, the defining characteristic of a metazoan species is a coadapted mitonuclear genotype that is incompatible with the coadapted mitochondrial and nuclear genotype of any other population.     

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.     

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.     

SEXUAL SELECTION AFFECTS THE EVOLUTION OF LIFESPAN AND AGEING IN THE DECORATED CRICKET GRYLLODES SIGILLATUS

Recent work suggests that sexual selection can influence the evolution of ageing and lifespan by shaping the optimal timing and relative costliness of reproductive effort in the sexes. We used inbred lines of the decorated cricket, Gryllodes sigillatus, to estimate the genetic (co)variance between age‐dependent reproductive effort, lifespan, and ageing within and between the sexes. Sexual selection theory predicts that males should die sooner and age more rapidly than females. However, a reversal of this pattern may be favored if reproductive effort increases with age in males but not in females. We found that male calling effort increased with age, whereas female fecundity decreased, and that males lived longer and aged more slowly than females. These divergent life‐history strategies were underpinned by a positive genetic correlation between early‐life reproductive effort and ageing rate in both sexes, although this relationship was stronger in females. Despite these sex differences in life‐history schedules, age‐dependent reproductive effort, lifespan, and ageing exhibited strong positive intersexual genetic correlations. This should, in theory, constrain the independent evolution of these traits in the sexes and may promote intralocus sexual conflict. Our study highlights the importance of sexual selection to the evolution of sex differences in ageing and lifespan in G. sigillatus.     

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.     

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.     

IGF‐I regulates the age‐dependent signaling peptide humanin

Aging is influenced by endocrine pathways including the growth hormone/insulin‐like growth factor‐1 (GH/IGF) axis. Mitochondrial function has also been linked to the aging process, but the relevant mitochondrial signals mediating the effects of mitochondria are poorly understood. Humanin is a novel signaling peptide that acts as a potent regulator of cellular stress responses and protects from a variety of in vitro and in vivo toxic and metabolic insults. The circulating levels of humanin decline with age in mice and humans. Here, we demonstrate a negative correlation between the activity of the GH‐IGF axis and the levels of humanin, as well as a positive correlation between humanin and lifespan in mouse models with altered GH/IGF‐I axis. Long‐lived, GH‐deficient Ames mice displayed elevated humanin levels, while short‐lived GH‐transgenic mice have reduced humanin levels. Furthermore, treatment with GH or IGF‐I reduced circulating humanin levels in both mice and human subjects. Our results indicate that GH and IGF are potent regulators of humanin levels and that humanin levels correlate with lifespan in mice. This suggests that humanin represents a circulating mitochondrial signal that participates in modulating the aging process, adding a coordinated mitochondrial element to the endocrine regulation of aging.     

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.     

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.     

Within‐population Y‐linked genetic variation for lifespan in Drosophila melanogaster

The view that the Y chromosome is of little importance for phenotypic evolution stems from early studies of Drosophila melanogaster. This species’ Y chromosome contains only 13 protein‐coding genes, is almost entirely heterochromatic and is not necessary for male viability. Population genetic theory further suggests that non‐neutral variation can only be maintained at the Y chromosome under special circumstances. Yet, recent studies suggest that the D. melanogaster Y chromosome trans‐regulates hundreds to thousands of X and autosomal genes. This finding suggests that the Y chromosome may play a far more active role in adaptive evolution than has previously been assumed. To evaluate the potential for the Y chromosome to contribute to phenotypic evolution from standing genetic variation, we test for Y‐linked variation in lifespan within a population of D. melanogaster. Assessing variation for lifespan provides a powerful test because lifespan (i) shows sexual dimorphism, which the Y is primarily predicted to contribute to, (ii) is influenced by many genes, which provides the Y with many potential regulatory targets and (iii) is sensitive to heterochromatin remodelling, a mechanism through which the Y chromosome is believed to regulate gene expression. Our results show a small but significant effect of the Y chromosome and thus suggest that the Y chromosome has the potential to respond to selection from standing genetic variation. Despite its small effect size, Y‐linked variation may still be important, in particular when evolution of sexual dimorphism is genetically constrained elsewhere in the genome.     

Stress profiling of longevity mutants identifies Afg3 as a mitochondrial determinant of cytoplasmic mRNA translation and aging

Although environmental stress likely plays a significant role in promoting aging, the relationship remains poorly understood. To characterize this interaction in a more comprehensive manner, we examined the stress response profiles for 46 long‐lived yeast mutant strains across four different stress conditions (oxidative, ER, DNA damage, and thermal), grouping genes based on their associated stress response profiles. Unexpectedly, cells lacking the mitochondrial AAA protease gene AFG3 clustered strongly with long‐lived strains lacking cytosolic ribosomal proteins of the large subunit. Similar to these ribosomal protein mutants, afg3Δ cells show reduced cytoplasmic mRNA translation, enhanced resistance to tunicamycin that is independent of the ER unfolded protein response, and Sir2‐independent but Gcn4‐dependent lifespan extension. These data demonstrate an unexpected link between a mitochondrial protease, cytoplasmic mRNA translation, and aging.     

Hybrid crosses and the genetic basis of interspecific divergence in lifespan in Pristionchus nematodes

Characterizing the genetic basis of among‐species variation in lifespan is a major goal of evolutionary gerontology research, but the very feature that defines separate species – the inability to interbreed – makes achieving this goal impractical, if not impossible, for most taxa. Pristionchus nematodes provide an intriguing system for tackling this problem, as female lifespan varies among species that can be crossed to form viable (although infertile) hybrids. By conducting reciprocal crosses among three species – two dioecious (long‐lived Pristionchus exspectatus and short‐lived Pristionchus arcanus) and one androdioecious (short‐lived Pristionchus pacificus) – we found that female lifespan was long for all hybrids, consistent with the hypothesis that the relatively short lifespans seen for P. pacificus hermaphrodites and P. arcanus females are caused by independent, recessive alleles that are masked in hybrid genomes. Cross‐direction had a small effect on survivorship for crosses involving P. exspectatus, indicating that nuclear–mitochondrial interactions may also influence Pristionchus longevity. Our findings suggest that long lifespan in P. exspectatus reflects the realization of an ancestral potential for extended longevity in the P. pacificus species complex. This work demonstrates the utility of interspecific hybrids for ageing research and provides a foundation for future work on the genetic architecture of interspecific lifespan variation.     

Differential response to larval crowding of a long‐ and a short‐lived medfly biotype

Response of endophytic fruit fly species (Tephritidae) to larval crowding is a form of scramble competition that may affect important life history traits of adults, such as survival and reproduction. Recent empirical evidence demonstrates large differences in adult life history traits, especially longevity, among Mediterranean fruit fly (Ceratitis capitata; "medfly") biotypes obtained from different regions of the world. However, whether the evolution of long lifespan is associated with response to stress induced by larval crowding has not been fully elucidated. We investigated, under constant laboratory conditions, the response of a short‐ and a long‐lived medfly biotypes to stress induced by larval crowding. Survival and development of larvae and pupae and the size of resulting pupae were recorded. The lifespan and age‐specific egg production patterns of the obtained adults were recorded. Our findings reveal that increased larval density reduced immature survival (larvae and pupae) in the short‐lived biotype but had rather neutral effects on the longed‐lived one. Only larvae of the long‐lived biotype were capable of prolonging their developmental duration under the highest crowding regime to successfully pupate and emerge as adults. Response of emerging adults to larvae crowding conditions was similar in the two medfly biotypes. Those individuals emerging from high larval density regimes had reduced longevity and fecundity. Long‐lived biotype individuals, however, appeared to suffer a higher cost in longevity compared with the short‐lived one. The importance of our findings to understand the evolution of long lifespan is discussed.     

A metabolic signature for long life in the Caenorhabditis elegans Mit mutants

Mit mutations that disrupt function of the mitochondrial electron transport chain can, inexplicably, prolong Caenorhabditis elegans lifespan. In this study we use a metabolomics approach to identify an ensemble of mitochondrial‐derived α‐ketoacids and α‐hydroxyacids that are produced by long‐lived Mit mutants but not by other long‐lived mutants or by short‐lived mitochondrial mutants. We show that accumulation of these compounds is dependent on concerted inhibition of three α‐ketoacid dehydrogenases that share dihydrolipoamide dehydrogenase (DLD) as a common subunit, a protein previously linked in humans with increased risk of Alzheimer's disease. When the expression of DLD in wild‐type animals was reduced using RNA interference we observed an unprecedented effect on lifespan – as RNAi dosage was increased lifespan was significantly shortened, but, at higher doses, it was significantly lengthened, suggesting that DLD plays a unique role in modulating length of life. Our findings provide novel insight into the origin of the Mit phenotype.     

Mitochondrial variation in small brown planthoppers linked to multiple traits and probably reflecting a complex evolutionary trajectory

While it has been proposed in several taxa that the mitochondrial genome is associated with adaptive evolution to different climatic conditions, making links between mitochondrial haplotypes and organismal phenotypes remains a challenge. Mitonuclear discordance occurs in the small brown planthopper (SBPH), Laodelphax striatellus, with one mitochondrial haplogroup (HGI) more common in the cold climate region of China relative to another form (HGII) despite strong nuclear gene flow, providing a promising model to investigate climatic adaptation of mitochondrial genomes. We hypothesized that cold adaptation through HGI may be involved, and considered mitogenome evolution, population genetic analyses, and bioassays to test this hypothesis. In contrast to our hypothesis, chill‐coma recovery tests and population genetic tests of selection both pointed to HGII being involved in cold adaptation. Phylogenetic analyses revealed that HGII is nested within HGI, and has three nonsynonymous changes in ND2, ND5 and CYTB in comparison to HGI. These molecular changes likely increased mtDNA copy number, cold tolerance and fecundity of SBPH, particularly through a function‐altering amino acid change involving M114T in ND2. Nuclear background also influenced fecundity and chill recovery (i.e., mitonuclear epistasis) and protein modelling indicates possible nuclear interactions for the two nonsynonymous changes in ND2 and CYTB. The high occurrence frequency of HGI in the cold climate region of China remains unexplained, but several possible reasons are discussed. Overall, our study points to a link between mtDNA variation and organismal‐level evolution and suggests a possible role of mitonuclear interactions in maintaining mtDNA diversity.