Evolutionary Rate Correlation between Mitochondrial-Encoded and Mitochondria-Associated Nuclear-Encoded Proteins in Insects

ERCnet: Phylogenomic Prediction of Interaction Networks in the Presence of Gene Duplication

Assigning gene function from genome sequences is a rate-limiting step in molecular biology research. A protein's position within an interaction network can potentially provide insights into its molecular mechanisms. Phylogenetic analysis of evolutionary rate covariation (ERC) in protein sequence has been shown to be effective for large-scale prediction of functional relationships and interactions. However, gene duplication, gene loss, and other sources of phylogenetic incongruence are barriers for analyzing ERC on a genome-wide basis. Here, we developed ERCnet, a bioinformatic program designed to overcome these challenges, facilitating efficient all-versus-all ERC analyses for large protein sequence datasets. We simulated proteome datasets and found that ERCnet achieves combined false positive and negative error rates well below 10% and that our novel "branch-by-branch" length measurements outperforms "root-to-tip" approaches in most cases, offering a valuable new strategy for performing ERC. We also compiled a sample set of 35 angiosperm genomes to test the performance of ERCnet on empirical data, including its sensitivity to user-defined analysis parameters such as input dataset size and branch-length measurement strategy. We investigated the overlap between ERCnet runs with different species samples to understand how species number and composition affect predicted interactions and to identify the protein sets that consistently exhibit ERC across angiosperms. Our systematic exploration of the performance of ERCnet provides a roadmap for design of future ERC analyses to predict functional interactions in a wide array of genomic datasets. ERCnet code is freely available at https://github.com/EvanForsythe/ERCnet.

Kingdom-specific lipid unsaturation calibrates sequence evolution in membrane arm subunits of eukaryotic respiratory complexes

Sequence evolution of protein complexes (PCs) is constrained by protein-protein interactions (PPIs). PPI-interfaces are predominantly conserved and hotspots for disease-related mutations. How do lipid-protein interactions (LPIs) constrain sequence evolution of membrane-PCs? We explore Respiratory Complexes (RCs) as a case study as these allow to compare sequence evolution in subunits exposed to both lipids in inner-mitochondrial membrane (IMM) and lipid-free aqueous matrix. We find that lipid-exposed surfaces of the IMM-subunits but not of the matrix subunits are populated with non-PPI disease-causing mutations signifying LPIs in stabilizing RCs. Further, IMM-subunits including their exposed surfaces show high intra-kingdom sequence conservation but remarkably diverge beyond. Molecular Dynamics simulation suggests contrasting LPIs of structurally superimposable but sequence-wise diverged IMM-exposed helices of Complex I (CI) subunit Ndufa1 from human and Arabidopsis depending on kingdom-specific unsaturation of cardiolipin fatty acyl chains. in cellulo assays consolidate inter-kingdom incompatibility of Ndufa1-helices due to the lipid-exposed amino acids. Plant-specific unsaturated fatty acids in human cells also trigger CI-instability. Taken together, we posit that altered LPIs calibrate sequence evolution at the IMM-arms of eukaryotic RCs.

Insect Mitochondrial Genomics: A Decade of Progress

The past decade has seen the availability of insect genomic data explode, with mitochondrial (mt) genome data seeing the greatest growth. The widespread adoption of next-generation sequencing has solved many earlier methodological limitations, allowing the routine sequencing of whole mt genomes, including from degraded or museum specimens and in parallel to nuclear genomic projects. The diversity of available taxa now allows finer-scale comparisons between mt and nuclear phylogenomic analyses; high levels of congruence have been found for most orders, with some significant exceptions (e.g., Odonata, Mantodea, Diptera). The evolution of mt gene rearrangements and their association with haplodiploidy have been tested with expanded taxonomic sampling, and earlier proposed trends have been largely supported. Multiple model systems have been developed based on findings unique to insects, including mt genome fragmentation (lice and relatives) and control region duplication (thrips), allowing testing of hypothesized evolutionary drivers of these aberrant genomic phenomena. Finally, emerging research topics consider the contributions of mt genomes to insect speciation and habitat adaption, with very broad potential impacts. Integration between insect mt genomic research and other fields within entomology continues to be our field's greatest opportunity and challenge.

Mitonuclear compatibility is maintained despite relaxed selection on male mitochondrial DNA in bivalves with doubly uniparental inheritance

Mitonuclear coevolution is common in eukaryotes, but bivalve lineages that have doubly uniparental inheritance (DUI) of mitochondria may be an interesting example. In this system, females transmit mtDNA (F mtDNA) to all offspring, while males transmit a different mtDNA (M mtDNA) solely to their sons. Molecular evolution and functional data suggest oxidative phosphorylation (OXPHOS) genes encoded in M mtDNA evolve under relaxed selection due to their function being limited to sperm only (vs. all other tissues for F mtDNA). This has led to the hypothesis that mitonuclear coevolution is less important for M mtDNA. Here, we use comparative phylogenetics, transcriptomics, and proteomics to understand mitonuclear interactions in DUI bivalves. We found nuclear OXPHOS proteins coevolve and maintain compatibility similarly with both F and M mtDNA OXPHOS proteins. Mitochondrial recombination did not influence mitonuclear compatibility and nuclear-encoded OXPHOS genes were not upregulated in tissues with M mtDNA to offset dysfunction. Our results support that selection maintains mitonuclear compatibility with F and M mtDNA despite relaxed selection on M mtDNA. Strict sperm transmission, lower effective population size, and higher mutation rates may explain the evolution of M mtDNA. Our study highlights that mitonuclear coevolution and compatibility may be broad features of eukaryotes.

Rapid evolution of mitochondrion-related genes in haplodiploid arthropods

BACKGROUND

Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes.

RESULTS

We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants.

CONCLUSIONS

Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans.

Novel Mito-Nuclear Combinations Facilitate the Global Invasion of a Major Agricultural Crop Pest

A fundamental understanding of the underlying mechanisms involved in biological invasions is crucial to developing effective risk assessment and control measures against invasive species. The fall armyworm (FAW), Spodoptera frugiperda, is a highly invasive pest that has rapidly spread from its native Americas into much of the Eastern Hemisphere, with a highly homogeneous nuclear genetic background. However, the exact mechanism behind its rapid introduction and propagation remains unclear. Here, a systematic investigation is conducted into the population dynamics of FAW in China from 2019 to 2021 and found that FAW individuals carrying "rice" mitochondria (FAW-mR) are more prevalent (>98%) than that with "corn" mitochondria (FAW-mC) at the initial stage of the invasion and in newly-occupied non-overwintering areas. Further fitness experiments show that the two hybrid-strains of FAW exhibit different adaptions in the new environment in China, and this may have been facilitated by amino acid changes in mitochondrial-encoded proteins. FAW-mR used increases energy metabolism, faster wing-beat frequencies, and lower wing loadings to drive greater flight performance and subsequent rapid colonization of new habitats. In contrast, FAW-mC individuals adapt with more relaxed mitochondria and shuttle energetics into maternal investment, observed as faster development rate and higher fecundity. The presence of two different mitochondria types within FAW has the potential to significantly expand the range of damage and enhance competitive advantage. Overall, the study describes a novel invasion mechanism displayed by the FAW population that facilitates its expansion and establishment in new environments.

Identification of a longevity gene through evolutionary rate covariation of insect mito-nuclear genomes

Oxidative phosphorylation, essential for energy metabolism and linked to the regulation of longevity, involves mitochondrial and nuclear genes. The functions of these genes and their evolutionary rate covariation (ERC) have been extensively studied, but little is known about whether other nuclear genes not targeted to mitochondria evolutionarily and functionally interact with mitochondrial genes. Here we systematically examined the ERC of mitochondrial and nuclear benchmarking universal single-copy ortholog (BUSCO) genes from 472 insects, identifying 75 non-mitochondria-targeted nuclear genes. We found that the uncharacterized gene CG11837-a putative ortholog of human DIMT1-regulates insect lifespan, as its knockdown reduces median lifespan in five diverse insect species and Caenorhabditis elegans, whereas its overexpression extends median lifespans in fruit flies and C. elegans and enhances oxidative phosphorylation gene activity. Additionally, DIMT1 overexpression protects human cells from cellular senescence. Together, these data provide insights into the ERC of mito-nuclear genes and suggest that CG11837 may regulate longevity across animals.

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Polyploid plants take cytonuclear perturbations in stride

Hybridization in plants is often accompanied by nuclear genome doubling (allopolyploidy), which has been hypothesized to perturb interactions between nuclear and organellar (mitochondrial and plastid) genomes by creating imbalances in the relative copy number of these genomes and producing genetic incompatibilities between maternally derived organellar genomes and the half of the allopolyploid nuclear genome from the paternal progenitor. Several evolutionary responses have been predicted to ameliorate these effects, including selection for changes in protein sequences that restore cytonuclear interactions; biased gene retention/expression/conversion favoring maternal nuclear gene copies; and fine-tuning of relative cytonuclear genome copy numbers and expression levels. Numerous recent studies, however, have found that evolutionary responses are inconsistent and rarely scale to genome-wide generalities. The apparent robustness of plant cytonuclear interactions to allopolyploidy may reflect features that are general to allopolyploids such as the lack of F2 hybrid breakdown under disomic inheritance, and others that are more plant-specific, including slow sequence divergence in organellar genomes and preexisting regulatory responses to changes in cell size and endopolyploidy during development. Thus, cytonuclear interactions may only rarely act as the main barrier to establishment of allopolyploid lineages, perhaps helping to explain why allopolyploidy is so pervasive in plant evolution.

Association Analysis Provides Insights into Plant Mitonuclear Interactions

Cytonuclear interaction refers to the complex and ongoing process of coevolution between nuclear and organelle genomes, which are responsible for cellular respiration, photosynthesis, lipid metabolism, etc. and play a significant role in adaptation and speciation. There have been a large number of studies to detect signatures of cytonuclear interactions. However, identification of the specific nuclear and organelle genetic polymorphisms that are involved in these interactions within a species remains relatively rare. The recent surge in whole genome sequencing has provided us an opportunity to explore cytonuclear interaction from a population perspective. In this study, we analyzed a total of 3,439 genomes from 7 species to identify signals of cytonuclear interactions by association (linkage disequilibrium) analysis of variants in both the mitochondrial and nuclear genomes across flowering plants. We also investigated examples of nuclear loci identified based on these association signals using subcellular localization assays, gene editing, and transcriptome sequencing. Our study provides a novel perspective on the investigation of cytonuclear coevolution, thereby enriching our understanding of plant fitness and offspring sterility.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Genomic Signatures of Mitonuclear Coevolution in Mammals

Mitochondrial (mt) and nuclear-encoded proteins are integrated in aerobic respiration, requiring co-functionality among gene products from fundamentally different genomes. Different evolutionary rates, inheritance mechanisms, and selection pressures set the stage for incompatibilities between interacting products of the two genomes. The mitonuclear coevolution hypothesis posits that incompatibilities may be avoided if evolution in one genome selects for complementary changes in interacting genes encoded by the other genome. Nuclear compensation, in which deleterious mtDNA changes are offset by compensatory nuclear changes, is often invoked as the primary mechanism for mitonuclear coevolution. Yet, direct evidence supporting nuclear compensation is rare. Here, we used data from 58 mammalian species representing eight orders to show strong correlations between evolutionary rates of mt and nuclear-encoded mt-targeted (N-mt) proteins, but not between mt and non-mt-targeted nuclear proteins, providing strong support for mitonuclear coevolution across mammals. N-mt genes with direct mt interactions also showed the strongest correlations. Although most N-mt genes had elevated dN/dS ratios compared to mt genes (as predicted under nuclear compensation), N-mt sites in close contact with mt proteins were not overrepresented for signs of positive selection compared to noncontact N-mt sites (contrary to predictions of nuclear compensation). Furthermore, temporal patterns of N-mt and mt amino acid substitutions did not support predictions of nuclear compensation, even in positively selected, functionally important residues with direct mitonuclear contacts. Overall, our results strongly support mitonuclear coevolution across ∼170 million years of mammalian evolution but fail to support nuclear compensation as the major mode of mitonuclear coevolution.

Mito-nuclear coevolution and phylogenetic artifacts: the case of bivalve mollusks

Mito-nuclear phylogenetic discordance in Bivalvia is well known. In particular, the monophyly of Amarsipobranchia (Heterodonta + Pteriomorphia), retrieved from mitochondrial markers, contrasts with the monophyly of Heteroconchia (Heterodonta + Palaeoheterodonta), retrieved from nuclear markers. However, since oxidative phosphorylation nuclear markers support the Amarsipobranchia hypothesis instead of the Heteroconchia one, interacting subunits of the mitochondrial complexes ought to share the same phylogenetic signal notwithstanding the genomic source, which is different from the signal obtained from other nuclear markers. This may be a clue of coevolution between nuclear and mitochondrial genes. In this work we inferred the phylogenetic signal from mitochondrial and nuclear oxidative phosphorylation markers exploiting different phylogenetic approaches and added two more datasets for comparison: genes of the glycolytic pathway and genes related to the biogenesis of regulative small noncoding RNAs. All trees inferred from mitochondrial and nuclear subunits of the mitochondrial complexes support the monophyly of Amarsipobranchia, regardless of the phylogenetic pipeline. However, not every single marker agrees with this topology: this is clearly visible in nuclear subunits that do not directly interact with the mitochondrial counterparts. Overall, our data support the hypothesis of a coevolution between nuclear and mitochondrial genes for the oxidative phosphorylation. Moreover, we suggest a relationship between mitochondrial topology and different nucleotide composition between clades, which could be associated to the highly variable gene arrangement in Bivalvia.

Mitochondrial Genomes of Two Asexual Trichogramma (Hymenoptera: Trichogrammatidae) Strains and Comparison with Their Sexual Relatives

Simple Summary

Sexual reproduction is dominant in animals, while asexual lineages are rare and evolutionarily short-lived. However, sexual reproduction has substantial costs, such as male production, inputs to courtship and mating, increased risk of predator exposure, and sexually transmitted diseases. A large body of theories has been proposed to explain the paradox of sex. One favored explanation is that asexuals are more likely to accumulate a greater number of deleterious mutations, known as Muller’s ratchet. Trichogramma is a genus of egg parasitoid wasps and is widely used as a biological control agent for agricultural and forest pests. With asexual lineages in at least 16 species, Trichogramma provides an excellent model to investigate the causes and consequences of asexual reproduction. In this study, we sequenced and assembled two asexual Trichogramma mitogenomes, representing two divergent origins of asexual reproduction. The asexual T. pretiosum is induced by the endosymbiont Wolbachia, while T. cacoeciae presumably originates from interspecific hybridization. To test Muller’s ratchet hypothesis, we compared these two asexual mitogenomes with their sexual relatives and found no association between asexual reproduction and mutation accumulation. This study provides a basis for further investigation into mitochondrial evolution and asexual reproduction in Trichogramma.

Abstract

Despite its substantial costs, sexual reproduction dominates in animals. One popular explanation for the paradox of sex is that asexual reproduction is more likely to accumulate deleterious mutations than sexual reproduction. To test this hypothesis, we compared the mitogenomes of two asexual wasp strains, Trichogramma cacoeciae and T. pretiosum, to their sexual relatives. These two asexual strains represent two different transition mechanisms in Trichogramma from sexual to asexual reproduction. Asexual T. pretiosum is induced by Wolbachia, while T. cacoeciae presumably originated from interspecific hybridization. We sequenced and assembled complete mitochondrial genomes of asexual T. cacoeciae and T. pretiosum. Compared to four sexual relatives, we found no evidence of higher mutation accumulation in asexual Trichogramma mitogenomes than in their sexual relatives. We also did not detect any relaxed selection in asexual Trichogramma mitogenomes. In contrast, the intensified selection was detected in Nad1 and Nad4 of the asexual T. pretiosum mitogenome, suggesting more purifying selection. In summary, no higher mitochondrial mutation accumulation was detected in these two asexual Trichogramma strains. This study provides a basis for further investigating mitochondrial evolution and asexual reproduction in Trichogramma.

Generation of a mitochondrial protein compendium in Dictyostelium discoideum

The social ameba Dictyostelium discoideum has emerged as a powerful model to study mitochondrial genetics and bioenergetics. However, a comprehensive inventory of mitochondrial proteins that is critical to understanding mitochondrial processes has yet to be curated. Here, we utilized high-throughput multiplexed protein quantitation and homology analyses to generate a high-confidence mitochondrial protein compendium consisting of 936 proteins. Our proteomic approach, which utilizes mass spectrometry in combination with mathematical modeling, was validated through mitochondrial targeting sequence prediction and live-cell imaging. Our final compendium consists of 936 proteins. Nearly, a third of D. discoideum mitochondrial proteins do not have homologs in humans, budding yeasts, or an ancestral alphaproteobacteria. Additionally, we leverage our compendium to highlight the complexity of metabolic reprogramming during starvation-induced development. Our compendium lays a foundation to investigate mitochondrial processes that are unique in ameba and to understand the functions of conserved mitochondrial proteins in D. discoideum.

Global Patterns of Subgenome Evolution in Organelle-Targeted Genes of Six Allotetraploid Angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary "mismatches" between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Rapid sequence evolution is associated with genetic incompatibilities in the plastid Clp complex

KEY MESSAGE

Replacing the native clpP1 gene in the Nicotiana plastid genome with homologs from different donor species showed that the extent of genetic incompatibilities depended on the rate of sequence evolution. The plastid caseinolytic protease (Clp) complex plays essential roles in maintaining protein homeostasis and comprises both plastid-encoded and nuclear-encoded subunits. Despite the Clp complex being retained across green plants with highly conserved protein sequences in most species, examples of extremely accelerated amino acid substitution rates have been identified in numerous angiosperms. The causes of these accelerations have been the subject of extensive speculation but still remain unclear. To distinguish among prevailing hypotheses and begin to understand the functional consequences of rapid sequence divergence in Clp subunits, we used plastome transformation to replace the native clpP1 gene in tobacco (Nicotiana tabacum) with counterparts from another angiosperm genus (Silene) that exhibits a wide range in rates of Clp protein sequence evolution. We found that antibiotic-mediated selection could drive a transgenic clpP1 replacement from a slowly evolving donor species (S. latifolia) to homoplasmy but that clpP1 copies from Silene species with accelerated evolutionary rates remained heteroplasmic, meaning that they could not functionally replace the essential tobacco clpP1 gene. These results suggest that observed cases of rapid Clp sequence evolution are a source of epistatic incompatibilities that must be ameliorated by coevolutionary responses between plastid and nuclear subunits.

The role of mitonuclear incompatibilities in allopatric speciation

Aerobic metabolism in eukaryotic cells requires extensive interactions between products of the nuclear and mitochondrial genomes. Rapid evolution of the mitochondrial genome, including fixation of both adaptive and deleterious mutations, creates intrinsic selection pressures favoring nuclear gene mutations that maintain mitochondrial function. As this process occurs independently in allopatry, the resulting divergence between conspecific populations can subsequently be manifest in mitonuclear incompatibilities in inter-population hybrids. Such incompatibilities, mitonuclear versions of Bateson-Dobzhansky-Muller incompatibilities that form the standard model for allopatric speciation, can potentially restrict gene flow between populations, ultimately resulting in varying degrees of reproductive isolation. The potential role of mitonuclear incompatibilities in speciation is further enhanced where mtDNA substitution rates are elevated compared to the nuclear genome and where population structure maintains allopatry for adequate time to evolve multiple mitonuclear incompatibilities. However, the fact that mitochondrial introgression occurs across species boundaries has raised questions regarding the efficacy of mitonuclear incompatibilities in reducing gene flow. Several scenarios now appear to satisfactorily explain this phenomenon, including cases where differences in mtDNA genetic load may drive introgression or where co-introgression of coadapted nuclear genes may support the function of introgressed mtDNA. Although asymmetries in reproductive isolation between taxa are consistent with mitonuclear incompatibilities, interactions between autosomes and sex chromosomes yield similar predictions that are difficult to disentangle. With regard to establishing reproductive isolation while in allopatry, existing studies clearly suggest that mitonuclear incompatibilities can contribute to the evolution of barriers to gene flow. However, there is to date relatively little definitive evidence supporting a primary role for mitonuclear incompatibilities in the speciation process.

Novel ACE2 protein interactions relevant to COVID-19 predicted by evolutionary rate correlations

Angiotensin-converting enzyme 2 (ACE2) is the cell receptor that the coronavirus SARS-CoV-2 binds to and uses to enter and infect human cells. COVID-19, the pandemic disease caused by the coronavirus, involves diverse pathologies beyond those of a respiratory disease, including micro-thrombosis (micro-clotting), cytokine storms, and inflammatory responses affecting many organ systems. Longer-term chronic illness can persist for many months, often well after the pathogen is no longer detected. A better understanding of the proteins that ACE2 interacts with can reveal information relevant to these disease manifestations and possible avenues for treatment. We have undertaken an approach to predict candidate ACE2 interacting proteins which uses evolutionary inference to identify a set of mammalian proteins that “coevolve” with ACE2. The approach, called evolutionary rate correlation (ERC), detects proteins that show highly correlated evolutionary rates during mammalian evolution. Such proteins are candidates for biological interactions with the ACE2 receptor. The approach has uncovered a number of key ACE2 protein interactions of potential relevance to COVID-19 pathologies. Some proteins have previously been reported to be associated with severe COVID-19, but are not currently known to interact with ACE2, while additional predicted novel ACE2 interactors are of potential relevance to the disease. Using reciprocal rankings of protein ERCs, we have identified strongly interconnected ACE2 associated protein networks relevant to COVID-19 pathologies. ACE2 has clear connections to coagulation pathway proteins, such as Coagulation Factor V and fibrinogen components FGA, FGB, and FGG, the latter possibly mediated through ACE2 connections to Clusterin (which clears misfolded extracellular proteins) and GPR141 (whose functions are relatively unknown). ACE2 also connects to proteins involved in cytokine signaling and immune response (e.g. XCR1, IFNAR2 and TLR8), and to Androgen Receptor (AR). The ERC prescreening approach has elucidated possible functions for relatively uncharacterized proteins and possible new functions for well-characterized ones. Suggestions are made for the validation of ERC-predicted ACE2 protein interactions. We propose that ACE2 has novel protein interactions that are disrupted during SARS-CoV-2 infection, contributing to the spectrum of COVID-19 pathologies.

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 d_N /d_S 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.

Evolutionary rate covariation identifies SLC30A9 (ZnT9) as a mitochondrial zinc transporter

Recent advances in genome sequencing have led to the identification of new ion and metabolite transporters, many of which have not been characterized. Due to the variety of subcellular localizations, cargo and transport mechanisms, such characterization is a daunting task, and predictive approaches focused on the functional context of transporters are very much needed. Here we present a case for identifying a transporter localization using evolutionary rate covariation (ERC), a computational approach based on pairwise correlations of amino acid sequence evolutionary rates across the mammalian phylogeny. As a case study, we find that poorly characterized transporter SLC30A9 (ZnT9) coevolves with several components of the mitochondrial oxidative phosphorylation chain, suggesting mitochondrial localization. We confirmed this computational finding experimentally using recombinant human SLC30A9. SLC30A9 loss caused zinc mishandling in the mitochondria, suggesting that under normal conditions it acts as a zinc exporter. We therefore propose that ERC can be used to predict the functional context of novel transporters and other poorly characterized proteins.

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.

Phenotypic Variation in Mitochondria-Related Performance Traits Across New Zealand Snail Populations

Mitochondrial function is critical for energy homeostasis and should shape how genetic variation in metabolism is transmitted through levels of biological organization to generate stability in organismal performance. Mitochondrial function is encoded by genes in two distinct and separately inherited genomes-the mitochondrial genome and the nuclear genome-and selection is expected to maintain functional mito-nuclear interactions. The documented high levels of polymorphism in genes involved in these mito-nuclear interactions and wide variation for mitochondrial function demands an explanation for how and why variability in such a fundamental trait is maintained. Potamopyrgus antipodarum is a New Zealand freshwater snail with coexisting sexual and asexual individuals and, accordingly, contrasting systems of separate vs. co-inheritance of nuclear and mitochondrial genomes. As such, this snail provides a powerful means to dissect the evolutionary and functional consequences of mito-nuclear variation. The lakes inhabited by P. antipodarum span wide environmental gradients, with substantial across-lake genetic structure and mito-nuclear discordance. This situation allows us to use comparisons across reproductive modes and lakes to partition variation in cellular respiration across genetic and environmental axes. Here, we integrated cellular, physiological, and behavioral approaches to quantify variation in mitochondrial function across a diverse set of wild P. antipodarum lineages. We found extensive across-lake variation in organismal oxygen consumption and behavioral response to heat stress and differences across sexes in mitochondrial membrane potential but few global effects of reproductive mode. Taken together, our data set the stage for applying this important model system for sexual reproduction and polyploidy to dissecting the complex relationships between mito-nuclear variation, performance, plasticity, and fitness in natural populations.

Genome-wide signatures of plastid-nuclear coevolution point to repeated perturbations of plastid proteostasis systems across angiosperms

Nuclear and plastid (chloroplast) genomes experience different mutation rates, levels of selection, and transmission modes, yet key cellular functions depend on their coordinated interactions. Functionally related proteins often show correlated changes in rates of sequence evolution across a phylogeny [evolutionary rate covariation (ERC)], offering a means to detect previously unidentified suites of coevolving and cofunctional genes. We performed phylogenomic analyses across angiosperm diversity, scanning the nuclear genome for genes that exhibit ERC with plastid genes. As expected, the strongest hits were highly enriched for genes encoding plastid-targeted proteins, providing evidence that cytonuclear interactions affect rates of molecular evolution at genome-wide scales. Many identified nuclear genes functioned in post-transcriptional regulation and the maintenance of protein homeostasis (proteostasis), including protein translation (in both the plastid and cytosol), import, quality control, and turnover. We also identified nuclear genes that exhibit strong signatures of coevolution with the plastid genome, but their encoded proteins lack organellar-targeting annotations, making them candidates for having previously undescribed roles in plastids. In sum, our genome-wide analyses reveal that plastid-nuclear coevolution extends beyond the intimate molecular interactions within chloroplast enzyme complexes and may be driven by frequent rewiring of the machinery responsible for maintenance of plastid proteostasis in angiosperms.

Evolutionary Rates are Correlated Between Buchnera Endosymbionts and the Mitochondrial Genomes of Their Aphid Hosts

The evolution of bacterial endosymbiont genomes is strongly influenced by host-driven selection. Factors affecting host genome evolution will potentially affect endosymbiont genomes in similar ways. One potential outcome is correlations in molecular rates between the genomes of the symbiotic partners. Recently, we presented the first evidence of such correlations between the mitochondrial genomes of cockroaches and the genomes of their endosymbiont (Blattabacterium cuenoti). Here we investigate whether similar patterns are found in additional host-symbiont partners. We use partial genome data from multiple strains of the bacterial endosymbionts Buchnera aphidicola and Sulcia muelleri, and the mitochondrial genomes of their sap-feeding insect hosts. Both endosymbionts show phylogenetic congruence with the mitochondria of their hosts, a result that is expected due to their identical mode of inheritance. We compared root-to-tip distances and branch lengths of phylogenetically independent species pairs. Both analyses showed a highly significant correlation of molecular rates between the genomes of Buchnera and the mitochondrial genomes of their hosts. A similar correlation was detected between Sulcia and their hosts, but was not statistically significant. Our results indicate that evolutionary rate correlations between hosts and long-term symbionts may be a widespread phenomenon.

Patterns of variation of mutation rates of mitochondrial and nuclear genes of gastropods

BACKGROUND

Although mitochondrial DNA (mtDNA) of many animals tends to mutate at higher rates than nuclear DNA (nuDNA), a recent survey of mutation rates of various animal groups found that the gastropod family Bradybaenidae (suborder Helicina) shows a nearly 40-fold difference in mutation rates of mtDNA ([Formula: see text]_m) and nuDNA ([Formula: see text]_n), while other gastropod taxa exhibit only two to five-fold differences. To determine if Bradybaenidae represents an outlier within Gastropoda, I compared estimated values of [Formula: see text]_m/[Formula: see text]_n of additional gastropod groups. In particular, I reconstructed mtDNA and nuDNA gene trees of 121 datasets that include members of various clades contained within the gastropod subclasses Caenogastropoda, Heterobranchia, Patellogastropoda, and Vetigastropoda and then used total branch length estimates of these gene trees to infer [Formula: see text]_m/[Formula: see text]_n.

RESULTS

Estimated values of [Formula: see text]_m/[Formula: see text]_n range from 1.4 to 91.9. Datasets that exhibit relatively large values of [Formula: see text]_m/[Formula: see text]_n (i.e., > 20), however, show relatively lower estimates of [Formula: see text]_n (and not elevated [Formula: see text]_m) in comparison to groups with lower values. These datasets also tend to contain sequences of recently diverged species. In addition, datasets with low levels of phylogenetic breadth (i.e., contain members of single genera or families) exhibit higher values of [Formula: see text]_m/[Formula: see text]_n than those with high levels (i.e., those that contain representatives of single superfamilies or higher taxonomic ranks).

CONCLUSIONS

Gastropods exhibit considerable variation in estimates of [Formula: see text]_m/[Formula: see text]_n. Large values of [Formula: see text]_m/[Formula: see text]_n that have been calculated for Bradybaenidae and other gastropod taxa may be overestimated due to possible sampling artifacts or processes that depress estimates of total molecular divergence of nuDNA in groups that recently diversified.

Episodic evolution of coadapted sets of amino acid sites in mitochondrial proteins

The rate of evolution differs between protein sites and changes with time. However, the link between these two phenomena remains poorly understood. Here, we design a phylogenetic approach for distinguishing pairs of amino acid sites that evolve concordantly, i.e., such that substitutions at one site trigger subsequent substitutions at the other; and also pairs of sites that evolve discordantly, so that substitutions at one site impede subsequent substitutions at the other. We distinguish groups of amino acid sites that undergo coordinated evolution and evolve discordantly from other such groups. In mitochondrion-encoded proteins of metazoans and fungi, we show that concordantly evolving sites are clustered in protein structures. By analysing the phylogenetic patterns of substitutions at concordantly and discordantly evolving site pairs, we find that concordant evolution has two distinct causes: epistatic interactions between amino acid substitutions and episodes of selection independently affecting substitutions at different sites. The rate of substitutions at concordantly evolving groups of protein sites changes in the course of evolution, indicating episodes of selection limited to some of the lineages. The phylogenetic positions of these changes are consistent between proteins, suggesting common selective forces underlying them.

The mode and rate of evolution of a protein site depends on the effect of its mutations on protein fitness. The fitness effect of a mutation itself can change in the course of evolution for at least two reasons. First, it can be modulated by substitutions occurring at other sites, a phenomenon called epistasis. Second, changes in selection can be non-epistatic, affecting sites independently of one another. Here, we analyse substitutions accumulated by the evolving lineages of the five proteins encoded by the mitochondrial genomes of thousands of species of metazoans and fungi. We show that substitutions at different amino acid sites occur in a coordinated fashion, and this coordination is caused both by epistasis and by episodes of selection affecting groups of sites. We partition each protein into several groups of concordantly evolving sites such that evolution of sites from different groups is discordant, and show that the proteins encoded by the mitochondrial genome consist of coevolving structural blocks. Some of these blocks have a clear functional specialization, e.g. are associated with interfaces between proteins composing respiratory complexes. Together, our results reveal a previously unrecognized complexity in the causes of variation in evolutionary rates between protein sites.

Comparative analysis reveals the expansion of mitochondrial DNA control region containing unusually high G-C tandem repeat arrays in Nasonia vitripennis

Insect mitochondrial DNA (mtDNA) ranges from 14 to 19 kbp, and the size difference is attributed to the AT-rich control region. Jewel wasps have a parasitoid lifestyle, which may affect mitochondria function and evolution. We sequenced, assembled, and annotated mitochondrial genomes in Nasonia and outgroup species. Gene composition and order are conserved within Nasonia, but they differ from other parasitoids by two large inversion events that were not reported before. We observed a much higher substitution rate relative to the nuclear genome and mitochondrial introgression between N. giraulti and N. oneida, which is consistent with previous studies. Most strikingly, N. vitripennis mtDNA has an extremely long control region (7665 bp), containing twenty-nine 217 bp tandem repeats and can fold into a super-cruciform structure. In contrast to tandem repeats commonly found in other mitochondria, these high-copy repeats are highly conserved (98.7% sequence identity), much longer in length (approximately 8 Kb), extremely GC-rich (50.7%), and CpG-rich (percent CpG 19.4% vs. 1.1% in coding region), resulting in a 23 kbp mtDNA beyond the typical size range in insects. These N. vitripennis-specific mitochondrial repeats are not related to any known sequences in insect mitochondria. Their evolutionary origin and functional consequences warrant further investigations.

Disrupted architecture and fast evolution of the mitochondrial genome of Argeia pugettensis (Isopoda): implications for speciation and fitness

BACKGROUND

Argeia pugettensis is an isopod species that parasitizes other crustaceans. Its huge native geographic range spans the Pacific from China to California, but molecular data are available only for a handful of specimens from North-American populations. We sequenced and characterised the complete mitogenome of a specimen collected in the Yellow Sea.

RESULTS

It exhibited a barcode (cox1) similarity level of only 87-89% with North-American populations, which is unusually low for conspecifics. Its mitogenome is among the largest in isopods (≈16.5 Kbp), mostly due to a large duplicated palindromic genomic segment (2 Kbp) comprising three genes. However, it lost a segment comprising three genes, nad4L-trnP-nad6, and many genes exhibited highly divergent sequences in comparison to isopod orthologues, including numerous mutations, deletions and insertions. Phylogenetic and selection analyses corroborated that this is one of the handful of most rapidly evolving available isopod mitogenomes, and that it evolves under highly relaxed selection constraints (as opposed to positive selection). However, its nuclear 18S gene is highly conserved, which suggests that rapid evolution is limited to its mitochondrial genome. The cox1 sequence analysis indicates that elevated mitogenomic evolutionary rates are not shared by North-American conspecifics, which suggests a breakdown of cox1 barcoding in this species.

CONCLUSIONS

A highly architecturally disrupted mitogenome and decoupling of mitochondrial and nuclear rates would normally be expected to have strong negative impacts on the fitness of the organism, so the existence of this lineage is a puzzling evolutionary question. Additional studies are needed to assess the phylogenetic breadth of this disrupted mitochondrial architecture and its impact on fitness.

An evolutionary, or "Mitocentric" perspective on cellular function and disease

The incidence of common, metabolic diseases (e.g. obesity, cardiovascular disease, diabetes) with complex genetic etiology has been steadily increasing nationally and globally. While identification of a genetic model that explains susceptibility and risk for these diseases has been pursued over several decades, no clear paradigm has yet been found to disentangle the genetic basis of polygenic/complex disease development. Since the evolution of the eukaryotic cell involved a symbiotic interaction between the antecedents of the mitochondrion and nucleus (which itself is a genetic hybrid), we suggest that this history provides a rational basis for investigating whether genetic interaction and co-evolution of these genomes still exists. We propose that both mitochondrial and Mendelian, or "mito-Mendelian" genetics play a significant role in cell function, and thus disease risk. This paradigm contemplates the natural variation and co-evolution of both mitochondrial and nuclear DNA backgrounds on multiple mitochondrial functions that are discussed herein, including energy production, cell signaling and immune response, which collectively can influence disease development. At the nexus of these processes is the economy of mitochondrial metabolism, programmed by both mitochondrial and nuclear genomes.

Convergent patterns of evolution of mitochondrial oxidative phosphorylation (OXPHOS) genes in electric fishes

The ability to generate and detect electric fields has evolved in several groups of fishes as a means of communication, navigation and, occasionally, predation. The energetic burden required can account for up to 20% of electric fishes' daily energy expenditure. Despite this, molecular adaptations that enable electric fishes to meet the metabolic demands of bioelectrogenesis remain unknown. Here, we investigate the molecular evolution of the mitochondrial oxidative phosphorylation (OXPHOS) complexes in the two most diverse clades of weakly electric fishes-South American Gymnotiformes and African Mormyroidea, using codon-based likelihood approaches. Our analyses reveal that although mitochondrial OXPHOS genes are generally subject to strong purifying selection, this constraint is significantly reduced in electric compared to non-electric fishes, particularly for complexes IV and V. Moreover, analyses of concatenated mitochondrial genes show strong evidence for positive selection in complex I genes on the two branches associated with the independent evolutionary origins of electrogenesis. These results suggest that adaptive evolution of proton translocation in the OXPHOS cellular machinery may be associated with the evolution of bioelectrogenesis. Overall, we find striking evidence for remarkably similar effects of electrogenesis on the molecular evolution of mitochondrial OXPHOS genes in two independently derived clades of electrogenic fishes. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Disentangling the intertwined roles of mutation, selection and drift in the mitochondrial genome

Understanding and quantifying the rates of change in the mitochondrial genome is a major component of many areas of biological inquiry, from phylogenetics to human health. A critical parameter in understanding rates of change is estimating the mitochondrial mutation rate (mtDNA MR). Although the first direct estimates of mtDNA MRs were reported almost 20 years ago, the number of estimates has not grown markedly since that time. This is largely owing to the challenges associated with time- and labour-intensive mutation accumulation (MA) experiments. But even MA experiments do not solve a major problem with estimating mtDNA MRs-the challenge of disentangling the role of mutation from other evolutionary forces acting within the cell. Now that it is widely understood that any newly generated mutant allele in the mitochondria will initially be at very low frequency (1/N, where N is the number of mtDNA molecules in the cell), the importance of understanding the effective population size (Ne) of the mtDNA and the size of genetic bottlenecks during gametogenesis and development has come into the spotlight. In addition to these factors regulating the role of genetic drift, advances in our understanding of mitochondrial replication and turnover allow us to more easily envision how natural selection within the cell might favour or purge mutations in multi-copy organellar genomes. Here, we review the unique features of the mitochondrial genome that pose a challenge for accurate MR estimation and discuss ways to overcome those challenges. Estimates of mtDNA MRs remain one of the most widely used parameters in biology, thus accurate quantification and a deeper understanding of how and why they may vary within and between individuals, populations and species is an important goal. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Evolutionary rates are correlated between cockroach symbionts and mitochondrial genomes

Bacterial endosymbionts evolve under strong host-driven selection. Factors influencing host evolution might affect symbionts in similar ways, potentially leading to correlations between the molecular evolutionary rates of hosts and symbionts. Although there is evidence of rate correlations between mitochondrial and nuclear genes, similar investigations of hosts and symbionts are lacking. Here, we demonstrate a correlation in molecular rates between the genomes of an endosymbiont (Blattabacterium cuenoti) and the mitochondrial genomes of their hosts (cockroaches). We used partial genome data for multiple strains of B. cuenoti to compare phylogenetic relationships and evolutionary rates for 55 cockroach/symbiont pairs. The phylogenies inferred for B. cuenoti and the mitochondrial genomes of their hosts were largely congruent, as expected from their identical maternal and cytoplasmic mode of inheritance. We found a correlation between evolutionary rates of the two genomes, based on comparisons of root-to-tip distances and on comparisons of the branch lengths of phylogenetically independent species pairs. Our results underscore the profound effects that long-term symbiosis can have on the biology of each symbiotic partner.

CyMIRA: The Cytonuclear Molecular Interactions Reference for Arabidopsis

The function and evolution of eukaryotic cells depend upon direct molecular interactions between gene products encoded in nuclear and cytoplasmic genomes. Understanding how these cytonuclear interactions drive molecular evolution and generate genetic incompatibilities between isolated populations and species is of central importance to eukaryotic biology. Plants are an outstanding system to investigate such effects because of their two different genomic compartments present in the cytoplasm (mitochondria and plastids) and the extensive resources detailing subcellular targeting of nuclear-encoded proteins. However, the field lacks a consistent classification scheme for mitochondrial- and plastid-targeted proteins based on their molecular interactions with cytoplasmic genomes and gene products, which hinders efforts to standardize and compare results across studies. Here, we take advantage of detailed knowledge about the model angiosperm Arabidopsis thaliana to provide a curated database of plant cytonuclear interactions at the molecular level. CyMIRA (Cytonuclear Molecular Interactions Reference for Arabidopsis) is available at http://cymira.colostate.edu/ and https://github.com/dbsloan/cymira and will serve as a resource to aid researchers in partitioning evolutionary genomic data into functional gene classes based on organelle targeting and direct molecular interaction with cytoplasmic genomes and gene products. It includes 11 categories (and 27 subcategories) of different cytonuclear complexes and types of molecular interactions, and it reports residue-level information for cytonuclear contact sites. We hope that this framework will make it easier to standardize, interpret, and compare studies testing the functional and evolutionary consequences of cytonuclear interactions.