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

Mitonuclear interactions impact aerobic metabolism in hybrids and may explain mitonuclear discordance in young, naturally hybridizing bird lineages

Understanding genetic incompatibilities and genetic introgression between incipient species are major goals in evolutionary biology. Mitochondrial genes evolve rapidly and exist in dense gene networks with coevolved nuclear genes, suggesting that mitochondrial respiration may be particularly susceptible to disruption in hybrid organisms. Mitonuclear interactions have been demonstrated to contribute to hybrid dysfunction between deeply divergent taxa crossed in the laboratory, but there are few empirical examples of mitonuclear interactions between younger lineages that naturally hybridize. Here, we use controlled hybrid crosses and high-resolution respirometry to provide the first experimental evidence in a bird that inter-lineage mitonuclear interactions impact mitochondrial aerobic metabolism. Specifically, respiration capacity of the two mitodiscordant backcrosses (with mismatched mitonuclear combinations) differs from one another, although they do not differ significantly from the parental groups or mitoconcordant backcrosses as we would expect of mitonuclear disruptions. In the wild hybrid zone between these subspecies, the mitochondrial cline centre is shifted west of the nuclear cline centre, which is consistent with the direction of our experimental results. Our results therefore demonstrate asymmetric mitonuclear interactions that impact the capacity of cellular mitochondrial respiration and may help to explain the geographic discordance between mitochondrial and nuclear genomes observed in the wild.

A lethal mitonuclear incompatibility in complex I of natural hybrids

The evolution of reproductive barriers is the first step in the formation of new species and can help us understand the diversification of life on Earth. These reproductive barriers often take the form of hybrid incompatibilities, in which alleles derived from two different species no longer interact properly in hybrids1-3. Theory predicts that hybrid incompatibilities may be more likely to arise at rapidly evolving genes4-6 and that incompatibilities involving multiple genes should be common7,8, but there has been sparse empirical data to evaluate these predictions. Here we describe a mitonuclear incompatibility involving three genes whose protein products are in physical contact within respiratory complex I of naturally hybridizing swordtail fish species. Individuals homozygous for mismatched protein combinations do not complete embryonic development or die as juveniles, whereas those heterozygous for the incompatibility have reduced complex I function and unbalanced representation of parental alleles in the mitochondrial proteome. We find that the effects of different genetic interactions on survival are non-additive, highlighting subtle complexity in the genetic architecture of hybrid incompatibilities. Finally, we document the evolutionary history of the genes involved, showing signals of accelerated evolution and evidence that an incompatibility has been transferred between species via hybridization.

Minuca panema (Coelho, 1972): Resurrection of a Fiddler Crab Species from Brazil Closely Related to Minuca burgersi (Holthuis, 1967) (Crustacea, Decapoda, Brachyura, Ocypodidae)

We redescribe a species of fiddler crab, Minuca panema (Coelho, 1972), from the Atlantic coast of South America. It is closely related to M. mordax (Smith, 1870), and until now, the taxon has been considered to be synonymous with another closely related species Minuca burgersi (Holthuis, 1967). However, we found that two clades of M. burgersi sensu lato were restricted to the Caribbean Basin. This distribution differs from than that of M. panema, which occurs primarily along the eastern coast of South America, ranging from the island of Trinidad to Praia da Armação, Santa Catarina, Brazil. Based on our field studies, the geographical boundary between M. burgersi sensu stricto and M. panema is the Tobago Basin, north of Trinidad. Since the two species diverged only 3 to 4 million years ago, as dated from the phylogeny of the genus Minuca Bott 1954, there are few reliable morphological features that can be used to distinguish them clearly. In live crabs, there is a striking difference in coloration between the cherryred South American M. panema and the rusty-red Caribbean M. burgersi sensu lato. In males, the pattern of tubercles on the inner surface of the major cheliped varies between the two species. In females, the vulva is slightly larger in M. burgersi sensu stricto. Ocean tides and currents together with siltation owing to freshwater outflow from the Amazon and Orinoco rivers most likely have driven the divergence of these species. In the Caribbean, small tidal amplitudes have minimized long-distance gene flow in M. burgersi sensu stricto from isolated insular lagoons. In contrast, large tidal amplitudes and exposed habitats on riverbanks along the eastern Atlantic coast of South America have enabled long-distance dispersal in M. panema. DNA analysis reveals that haplotypes of cytochrome c oxidase subunit 1 are not shared between the species. Since natural selection and/or genetic drift have yet to produce extensive morphological divergences between M. panema and M. burgersi sensu stricto, we speculate that changes in the genes regulating mitochondrial DNA functions have led to speciation at the molecular level. According to the mitonuclear compatibility concept, we propose that mitochondrial DNA may be at the forefront of speciation events and that co-evolved mitonuclear interactions are responsible for some of the earliest genetic incompatibilities arising among isolated populations.

High quality de novo genome assembly of the non-conventional yeast Kazachstania bulderi describes a potential low pH production host for biorefineries

Kazachstania bulderi is a non-conventional yeast species able to grow efficiently on glucose and δ-gluconolactone at low pH. These unique traits make K. bulderi an ideal candidate for use in sustainable biotechnology processes including low pH fermentations and the production of green chemicals including organic acids. To accelerate strain development with this species, detailed information of its genetics is needed. Here, by employing long read sequencing we report a high-quality phased genome assembly for three strains of K. bulderi species, including the type strain. The sequences were assembled into 12 chromosomes with a total length of 14 Mb, and the genome was fully annotated at structural and functional levels, including allelic and structural variants, ribosomal array and mating type locus. This high-quality reference genome provides a resource to advance our fundamental knowledge of biotechnologically relevant non-conventional yeasts and to support the development of genetic tools for manipulating such strains towards their use as production hosts in biotechnological processes.

Genetic incompatibilities in reciprocal hybrids between populations of Tigriopus californicus with low to moderate mitochondrial sequence divergence

All mitochondrial-encoded proteins and RNAs function through interactions with nuclear-encoded proteins, which are critical for mitochondrial performance and eukaryotic fitness. Coevolution maintains inter-genomic (i.e., mitonuclear) compatibility within a taxon, but hybridization can disrupt coevolved interactions, resulting in hybrid breakdown. Thus, mitonuclear incompatibilities may be important mechanisms underlying reproductive isolation and, potentially, speciation. Here we utilize Pool-seq to assess the effects of mitochondrial genotype on nuclear allele frequencies in fast- and slow-developing reciprocal inter-population F2 hybrids between relatively low-divergence populations of the intertidal copepod Tigriopus californicus. We show that mitonuclear interactions lead to elevated frequencies of coevolved (i.e., maternal) nuclear alleles on two chromosomes in crosses between populations with 1.5% or 9.6% fixed differences in mitochondrial DNA nucleotide sequence. However, we also find evidence of excess mismatched (i.e., noncoevolved) alleles on three or four chromosomes per cross, respectively, and of allele frequency differences consistent with effects involving only nuclear loci (i.e., unaffected by mitochondrial genotype). Thus, our results for low-divergence crosses suggest an underlying role for mitonuclear interactions in variation in hybrid developmental rate, but despite substantial effects of mitonuclear coevolution on individual chromosomes, no clear bias favoring coevolved interactions overall.

The Rhyacophila fasciata Species Complex (Trichoptera: Rhyacophilidae) in Central Europe with description of a new species, Rhyacophila loeffleri Valladolid & Waringer, n. sp., based on morphological, genetic and ecological evidence

The morphology of the larva, pupa, male, and female of Rhyacophila loeffleri Valladolid & Waringer, n. sp., from Austria, Czech Republic, France, and Italy is described. Data from molecular analysis are also provided, based on mitochondrial cytochrome c oxidase subunit I (mtCOI), including sequences from German and Swiss specimens, and compared with those of samples of the nominate species R. fasciata Hagen 1859, as well as with other species and subspecies in the Rhyacophila fasciata Species Complex: R. septentrionis McLachlan 1865, R. denticulata McLachlan 1879, R. sociata Navs 1916, R. kykladica Malicky & Sipahiler 1993, R. delici Kuini & Valladolid 2020, R. fasciata viteceki Valladolid & Kuini 2020 and R. macedonica Karaouzas, Valladolid & Ibrahimi 2022. Our study revealed morphological differences between R. fasciata and R. loeffleri, as well as genetic differences among the other taxa of the R. fasciata Complex, supporting the recognition of R. loeffleri as a new species. Also, Rhyacophila fasciata viteceki is recognized as a distinct species, R. viteceki (status promotus), based on results of the molecular analysis.

A chromosome-scale genome assembly and evaluation of mtDNA variation in the willow leaf beetle Chrysomela aeneicollis

The leaf beetle Chrysomela aeneicollis has a broad geographic range across Western North America but is restricted to cool habitats at high elevations along the west coast. Central California populations occur only at high altitudes (2,700-3,500 m) where they are limited by reduced oxygen supply and recent drought conditions that are associated with climate change. Here, we report a chromosome-scale genome assembly alongside a complete mitochondrial genome and characterize differences among mitochondrial genomes along a latitudinal gradient over which beetles show substantial population structure and adaptation to fluctuating temperatures. Our scaffolded genome assembly consists of 21 linkage groups; one of which we identified as the X chromosome based on female/male whole genome sequencing coverage and orthology with Tribolium castaneum. We identified repetitive sequences in the genome and found them to be broadly distributed across all linkage groups. Using a reference transcriptome, we annotated a total of 12,586 protein-coding genes. We also describe differences in putative secondary structures of mitochondrial RNA molecules, which may generate functional differences important in adaptation to harsh abiotic conditions. We document substitutions at mitochondrial tRNA molecules and substitutions and insertions in the 16S rRNA region that could affect intermolecular interactions with products from the nuclear genome. This first chromosome-level reference genome will enable genomic research in this important model organism for understanding the biological impacts of climate change on montane insects.

Mitogenome-based phylogenomics provides insights into the positions of the enigmatic sinensis group and the sanguinolenta group in Selaginellaceae (Lycophyte)

Spikemoss (Selaginellaceae) is one of the basal lineages of vascular plants. This family has a single genus Selaginella which consists of about 750 extant species. The phylogeny of Selaginellaceae has been extensively studied mainly based on plastid DNA and a few nuclear sequences. However, the placement of the enigmatic sinensis group is a long-term controversy because of the long branch in the plastid DNA phylogeny. The sanguinolenta group is also a phylogenetically problematic clade owing to two alternative positions resulted from different datasets. Here, we newly sequenced 34 mitochondrial genomes (mitogenomes) of individuals representing all seven subgenera and major clades in Selaginellaceae. We assembled the draft mitogenomes and annotated the genes and performed phylogenetic analyses based on the shared 17 mitochondrial genes. Our major results include: (1) all the assembled mitogenomes have complicated structures, unparalleled high GC content and a small gene content set, and the positive correlations among GC content, substitution rates and the number of RNA editing sites hold; (2) the sinensis group was well supported as a member of subg. Stachygynandrum; (3) the sanguinolenta group was strongly resolved as sister to all other Selaginella species except for subg. Selaginella. This study demonstrates the potential of mitogenome data in providing novel insights into phylogenetically recalcitrant problems.

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.

Mito-nuclear selection induces a trade-off between species ecological dominance and evolutionary lifespan

Mitochondrial and nuclear genomes must be co-adapted to ensure proper cellular respiration and energy production. Mito-nuclear incompatibility reduces individual fitness and induces hybrid infertility, which can drive reproductive barriers and speciation. Here, we develop a birth-death model for evolution in spatially extended populations under selection for mito-nuclear co-adaptation. Mating is constrained by physical and genetic proximity, and offspring inherit nuclear genomes from both parents, with recombination. The model predicts macroscopic patterns including a community's species diversity, species abundance distribution, speciation and extinction rates, as well as intraspecific and interspecific genetic variation. We explore how these long-term outcomes depend upon the parameters of reproduction: individual fitness governed by mito-nuclear compatibility, constraints on mating compatibility and ecological carrying capacity. We find that strong selection for mito-nuclear compatibility reduces the equilibrium number of species after a radiation, increasing species' abundances and simultaneously increasing both speciation and extinction rates. The negative correlation between species diversity and diversification rates in our model agrees with the broad empirical pattern of lower diversity and higher speciation/extinction rates in temperate regions, compared to the tropics. We conclude that these empirical patterns may be caused in part by latitudinal variation in metabolic demands and corresponding variation in selection for mito-nuclear function.

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.

Widespread genomic signatures of reproductive isolation and sex-specific selection in the Eastern Yellow Robin, Eopsaltria australis

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane’s Rule.

Life barcoded by DNA barcodes

The modern concept of DNA-based barcoding for cataloguing biodiversity was proposed in 2003 by first adopting an approximately 600 bp fragment of the mitochondrial COI gene to compare via nucleotide alignments with known sequences from specimens previously identified by taxonomists. Other standardized regions meeting barcoding criteria then are also evolving as DNA barcodes for fast, reliable and inexpensive assessment of species composition across all forms of life, including animals, plants, fungi, bacteria and other microorganisms. Consequently, global DNA barcoding campaigns have resulted in the formation of many online workbenches and databases, such as BOLD system, as barcode references, and facilitated the development of mini-barcodes and metabarcoding strategies as important extensions of barcode techniques. Here we intend to give an overview of the characteristics and features of these barcode markers and major reference libraries existing for barcoding the planet’s life, as well as to address the limitations and opportunities of DNA barcodes to an increasingly broader community of science and society.

Mitonuclear interplay in yeast: from speciation to phenotypic adaptation

Saccharomyces yeasts have evolved into an important model system to study mitonuclear incompatibilities, thanks to recent advances in the field of sequencing, yeast hybridisation and multigenerational breeding. Yeast hybrids contain two homologous proteomes but retain only one type of mitochondria allowing studies on the effect of mitochondria on phenotype and gene expression. Here, we discuss the recent developments in the growing field of yeast mitogenomics spanning from the impact that this organelle has in shaping yeast fitness and genome evolution to the dissection of molecular determinants of mitonuclear incompatibilities. Applying the state-of-the-art genetic tools to a broader range of natural yeast species from different environments will help progress the field and untap the mitochondrial potential in strain development.

A practical approach to revise the Oxynoemacheilus bergianus species group (Teleostei: Nemacheilidae)

The Oxynoemacheilus bergianus species group is revised based on tree topology (ML, NJ, MP), distance (K2P and ASAP) and Poisson tree process analyses of DNA barcode data tested against morphometric and morphological characters including colour patterns. The O. bergianus species group is distinguished from other Oxynoemacheilus groups based on morphological characters: its constituent species have a slender caudal peduncle, a suborbital flap in the male, a mottled or blotched colour pattern, and lack bold, black spots on the caudal-fin base. It is also supported as a monophyletic unit in our molecular analysis.                 The O. bergianus group includes 10 molecular clades following congruently well-supported NJ, MP and ML based entities. Species described as O. bergianus, O. banarescui, O. erdali, O. fatsaensis, O. samanticus, and O. simavicus from Turkey, O. lenkoranensis from Azerbaijan, and O. longipinnis and O. parvinae from Iran belong to this species group. The group includes also four unnamed molecular clades. We were unable to detect external differences between any of the molecular clades in colour pattern or any morphometric or morphological characters examined. In the 10 molecular clades in the O. bergianus species group, the intraclade K2P distance ranges from 0.0–1.8% while the distances between molecular clades ranges from 0.6–5.9%. To resolve the species diversity of this group, we also analysed the intraspecific and interspecific variability in the K2P distance of DNA barcode data from 53 other Oxynoemacheilus species. Here, the intraspecific variability ranges from 0.0–2.4% while the interspecific K2P distance ranges from 1.2–20.8%. In the O. bergianus species group, only four groups are detected by the mPTP species delimitation approach distinguished by a K2P distance of 2.9% or more. We treat these four groups as valid species, corresponding to O. banarescui, O. bergianus, O. fatsaensis, and O. simavicus. Oxynoemacheilus samanticus from the Kızılırmak and Seyhan drainages, O. lenkoranensis from the Caspian basin, O. erdali from the Euphrates, and O. longipinnis and O. parvinae from the Tigris drainage are treated as synonyms of O. bergianus. Fishes from an unnamed molecular clade from the upper Tigris, and from a second unnamed clade from the upper Euphrates, are both identified as O. bergianus. Oxynoemacheilus bergianus might be a junior synonym of O. bergi from the Kura. The distribution range of O. simavicus, described from the Simav drainage in the Marmara basin, is expanded to the east and two molecularly differentiated population groups occur in the Sakarya drainage, the Büyük Melen River and potentially in other adjacent coastal streams. Oxynoemacheilus fatsaensis, described from the coastal stream Elekçi in northern Anatolia, is also widespread in the Yeşilırmak drainage. Morphological characters proposed to distinguish O. fatsaensis from the other species of the O. bergianus group could not be confirmed by our data on fishes from the Yeşilırmak. This study also discusses the theoretical background, our reasons for conducting this revision in the way we did, and what the alternatives would be.  

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.

Genome and cuticular hydrocarbon-based species delimitation shed light on potential drivers of speciation in a Neotropical ant species complex

Geographic separation that leads to the evolution of reproductive isolation between populations generally is considered the most common form of speciation. However, speciation may also occur in the absence of geographic barriers due to phenotypic and genotypic factors such as chemical cue divergence, mating signal divergence, and mitonuclear conflict. Here, we performed an integrative study based on two genome-wide techniques (3RAD and ultraconserved elements) coupled with cuticular hydrocarbon (CHC) and mitochondrial (mt) DNA sequence data, to assess the species limits within the Ectatomma ruidum species complex, a widespread and conspicuous group of Neotropical ants for which heteroplasmy (i.e., presence of multiple mtDNA variants in an individual) has been recently discovered in some populations from southeast Mexico. Our analyses indicate the existence of at least five distinct species in this complex: two widely distributed across the Neotropics, and three that are restricted to southeast Mexico and that apparently have high levels of heteroplasmy. We found that species boundaries in the complex did not coincide with geographic barriers. We therefore consider possible roles of alternative drivers that may have promoted the observed patterns of speciation, including mitonuclear incompatibility, CHC differentiation, and colony structure. Our study highlights the importance of simultaneously assessing different sources of evidence to disentangle the species limits of taxa with complicated evolutionary histories.

Predictors of genomic differentiation within a hybrid taxon

Hybridization is increasingly recognized as an important evolutionary force. Novel genetic methods now enable us to address how the genomes of parental species are combined in hybrid lineages. However, we still do not know the relative importance of admixed proportions, genome architecture and local selection in shaping hybrid genomes. Here, we take advantage of the genetically divergent island populations of Italian sparrow on Crete, Corsica and Sicily to investigate the predictors of genomic variation within a hybrid taxon. We test if differentiation is affected by recombination rate, selection, or variation in ancestry proportions. We find that the relationship between recombination rate and differentiation is less pronounced within hybrid lineages than between the parent species, as expected if purging of minor parent ancestry in low recombination regions reduces the variation available for differentiation. In addition, we find that differentiation between islands is correlated with differences in signatures of selection in two out of three comparisons. Signatures of selection within islands are correlated across all islands, suggesting that shared selection may mould genomic differentiation. The best predictor of strong differentiation within islands is the degree of differentiation from house sparrow, and hence loci with Spanish sparrow ancestry may vary more freely. Jointly, this suggests that constraints and selection interact in shaping the genomic landscape of differentiation in this hybrid species.

Genomes of hybrid lineages are mosaics of those of their parent species and harbour variation that has the potential to facilitate adaptation when hybrids encounter diverse environments. However, genetic incompatibilities between parental species can also act to limit possible combinations of parental alleles, constraining hybrid genome formation. What is the relative importance of selection and constraints in form of admixture proportions and genomic architecture in this process? We investigated this in the Italian sparrow, a hybrid species resulting from past hybridization between the house and Spanish sparrow. Using three independent hybrid lineages, we addressed how their genomes, harbouring different parental combinations, have evolved. We examined the roles of selection due to divergent local adaptation, recombination and purging of genetic incompatibilities in predicting differentiation. We found that selection against incompatibilities may constrain hybrid genome composition. In addition, signals of local selection as well as estimates of differentiation were correlated across populations, and outliers were shared among the hybrid lineages more often than expected by chance. Overall, our results suggest that in the Italian sparrow selection interacts with constraints linked to genetic incompatibilities affecting which sections of the genome can readily diverge among hybrid lineages.

Migration restores hybrid incompatibility driven by mitochondrial-nuclear sexual conflict

In the mitochondrial genome, sexual asymmetry in transmission allows the accumulation of male-harming mutations since selection acts only on the effect of the mutation in females. Called the 'Mother's Curse', this phenomenon induces a selective pressure for nuclear variants that compensate for this reduction in male fitness. Previous work has demonstrated the existence of these interactions and their potential to act as Dobzhansky-Muller incompatibilities, contributing to reproductive isolation between populations. However, it is not clear how readily they would give rise to and sustain hybrid incompatibilities. Here, we use computer simulations in SLiM 3 to investigate the consequences of sexually antagonistic mitochondrial-nuclear interactions in a subdivided population. We consider distinct migration schemes and vary the chromosomal location, and consequently the transmission pattern, of nuclear restorers. Disrupting these co-evolved interactions results in less-fit males, skewing the sex ratio toward females. Restoration of male fitness depends on both the chromosomal location of nuclear restorer loci and the migration scheme. Our results show that these interactions may act as Dobzhansky-Muller incompatibilities, but their strength is not enough to drive population isolation. Overall, this model shows the varied ways in which populations can respond to migration's disruption of co-evolved mitochondrial-nuclear interactions.

Mitonuclear incompatibility as a hidden driver behind the genome ancestry of African admixed cattle

Background

Africa is an important watershed in the genetic history of domestic cattle, as two lineages of modern cattle, Bos taurus and B. indicus, form distinct admixed cattle populations. Despite the predominant B. indicus nuclear ancestry of African admixed cattle, B. indicus mitochondria have not been found on the continent. This discrepancy between the mitochondrial and nuclear genomes has been previously hypothesized to be driven by male-biased introgression of Asian B. indicus into ancestral African B. taurus. Given that this hypothesis mandates extreme demographic assumptions relying on random genetic drift, we propose a novel hypothesis of selection induced by mitonuclear incompatibility and assess these hypotheses with regard to the current genomic status of African admixed cattle.

Results

By analyzing 494 mitochondrial and 235 nuclear genome sequences, we first confirmed the genotype discrepancy between mitochondrial and nuclear genome in African admixed cattle: the absence of B. indicus mitochondria and the predominant B. indicus autosomal ancestry. We applied approximate Bayesian computation (ABC) to assess the posterior probabilities of two selection hypotheses given this observation. The results of ABC indicated that the model assuming both male-biased B. indicus introgression and selection induced by mitonuclear incompatibility explains the current genomic discrepancy most accurately. Subsequently, we identified selection signatures at autosomal loci interacting with mitochondria that are responsible for integrity of the cellular respiration system. By contrast with B. indicus-enriched genome ancestry of African admixed cattle, local ancestries at these selection signatures were enriched with B. taurus alleles, concurring with the key expectation of selection induced by mitonuclear incompatibility.

Conclusions

Our findings support the current genome status of African admixed cattle as a potential outcome of male-biased B. indicus introgression, where mitonuclear incompatibility exerted selection pressure against B. indicus mitochondria. This study provides a novel perspective on African cattle demography and supports the role of mitonuclear incompatibility in the hybridization of mammalian species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01206-x.

Sea surface temperature, rather than land mass or geographic distance, may drive genetic differentiation in a species complex of highly dispersive seabirds

Seabirds, particularly Procellariiformes, are highly mobile organisms with a great capacity for long dispersal, though simultaneously showing high philopatry, two conflicting life-history traits that may lead to contrasted patterns of genetic population structure. Landmasses were suggested to explain differentiation patterns observed in seabirds, but philopatry, isolation by distance, segregation between breeding and nonbreeding zones, and oceanographic conditions (sea surface temperatures) may also contribute to differentiation patterns. To our knowledge, no study has simultaneously contrasted the multiple factors contributing to the diversification of seabird species, especially in the gray zone of speciation. We conducted a multilocus phylogeographic study on a widespread seabird species complex, the little shearwater complex, showing highly homogeneous morphology, which led to considerable taxonomic debate. We sequenced three mitochondrial and six nuclear markers on all extant populations from the Atlantic (lherminieri) and Indian Oceans (bailloni), that is, five nominal lineages from 13 populations, along with one population from the eastern Pacific Ocean (representing the dichrous lineage). We found sharp differentiation among populations separated by the African continent with both mitochondrial and nuclear markers, while only mitochondrial markers allowed characterizing the five nominal lineages. No differentiation could be detected within these five lineages, questioning the strong level of philopatry showed by these shearwaters. Finally, we propose that Atlantic populations likely originated from the Indian Ocean. Within the Atlantic, a stepping-stone process accounts for the current distribution. Based on our divergence time estimates, we suggest that the observed pattern of differentiation mostly resulted from historical and current variation in sea surface temperatures.

Estimating the genetic diversity of Pacific salmon and trout using multigene eDNA metabarcoding

Genetic diversity underpins species conservation and management goals, and ultimately determines a species' ability to adapt. Using freshwater environmental DNA (eDNA) samples, we examined mitochondrial genetic diversity using multigene metabarcode sequence data from four Oncorhynchus species across 16 sites in Oregon and northern California. Our multigene metabarcode panel included targets commonly used in population genetic NADH dehydrogenase 2 (ND2), phylogenetic cytochrome c oxidase subunit 1 (COI) and eDNA (12S ribosomal DNA) screening. The ND2 locus showed the greatest within-species haplotype diversity for all species, followed by COI and then 12S rDNA for all species except Oncorhynchus kisutch. Sequences recovered for O. clarkii clarkii were either identical to, or one mutation different from, previously characterized haplotypes (95.3% and 4.5% of reads, respectively). The greatest diversity in O. c. clarkii was among coastal watersheds, and subsets of this diversity were shared with fish in inland watersheds. However, coastal streams and the Umpqua River watershed appear to harbour unique haplotypes. Sequences from O. mykiss revealed a disjunction between the Willamette watershed and southern watersheds suggesting divergent histories. We also identified similarities between populations in the northern Deschutes and southern Klamath watersheds, consistent with previously hypothesized connections between the two via inland basins. Oncorhynchus kisutch was only identified in coastal streams and the Klamath River watershed, with most diversity concentrated in the coastal Coquille watershed. Oncorhynchus tshawytscha was only observed at one site, but contained multiple haplotypes at each locus. The characterization of genetic diversity at multiple loci expands the knowledge gained from eDNA sampling and provides crucial information for conservation actions and genetic management.

Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats

Background

Mitochondrial function involves the interplay between mitochondrial and nuclear genomes. Such mitonuclear interactions can be disrupted by the introgression of mitochondrial DNA between taxa or divergent populations. Previous studies of several model systems (e.g. Drosophila) indicate that the disruption of mitonuclear interactions, termed mitonuclear mismatch, can alter nuclear gene expression, yet few studies have focused on natural populations.

Results

Here we study a naturally introgressed population in the secondary contact zone of two subspecies of the intermediate horseshoe bat (Rhinolophus affinis), in which individuals possess either mitonuclear matched or mismatched genotypes. We generated transcriptome data for six tissue types from five mitonuclear matched and five mismatched individuals. Our results revealed strong tissue-specific effects of mitonuclear mismatch on nuclear gene expression with the largest effect seen in pectoral muscle. Moreover, consistent with the hypothesis that genes associated with the response to oxidative stress may be upregulated in mitonuclear mismatched individuals, we identified several such gene candidates, including DNASE1L3, GPx3 and HSPB6 in muscle, and ISG15 and IFI6 in heart.

Conclusion

Our study reveals how mitonuclear mismatch arising from introgression in natural populations is likely to have fitness consequences. Underlying the processes that maintain mitonuclear discordance is a step forward to understand the role of mitonuclear interactions in population divergence and speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00424-x.

A response to Justen et al. 2020: Estimating hybridization rates in the wild: Easier said than done?

We consider four key challenges related to estimating per-individual rates of hybridization in wild birds: (1) what is the meaning of the term "hybrid"?, (2) the importance of distinguishing between shared DNA sequences and on-going hybridization between populations, (3) the perils of focusing exclusively on known hybrid zones, and (4) the implications of very low rates of per individual hybridization. Because our focus is on using phenotype to recognize hybrids, we define a hybrid as an individual with a phenotype that is intermediate between two parental species. We emphasize the value of quantifying the rate of between-species mating that is ongoing in current populations and distinguish this endeavor from estimates of gene flow between populations based on genomic analysis, which can reflect both current and ancient hybridization. We restate the importance of quantifying per individual rates of hybridization among all birds without prejudging which birds are presumed to engage in hybridization. And finally, we emphasize that evidence for strong prezygotic sorting is not necessarily evidence that mate choice is a driver of speciation.

Recovery from hybrid breakdown reveals a complex genetic architecture of mitonuclear incompatibilities

Reproductive isolation is often achieved when genes that are neutral or beneficial in their genomic background become functionally incompatible in a foreign genomic background, causing inviability, sterility or other forms of low fitness in hybrids. Recent studies suggest that mitonuclear interactions are among the initial incompatibilities to evolve at early stages of population divergence across taxa. Yet, the genomic architecture of mitonuclear incompatibilities has rarely been elucidated. We employ an experimental evolution approach starting with low-fitness F₂ interpopulation hybrids of the copepod Tigriopus californicus, in which frequencies of compatible and incompatible nuclear alleles change in response to an alternative mitochondrial background. After about nine generations, we observe a generalized increase in population size and in survivorship, suggesting efficiency of selection against maladaptive phenotypes. Whole genome sequencing of evolved populations showed some consistent allele frequency changes across three replicates of each reciprocal cross, but markedly different patterns between mitochondrial backgrounds. In only a few regions (~6.5% of the genome), the same parental allele was overrepresented irrespective of the mitochondrial background. About 33% of the genome showed allele frequency changes consistent with divergent selection, with the location of these genomic regions strongly differing between mitochondrial backgrounds. In 87% and 89% of these genomic regions, the dominant nuclear allele matched the associated mitochondrial background, consistent with mitonuclear co-adaptation. These results suggest that mitonuclear incompatibilities have a complex polygenic architecture that differs between populations, potentially generating genome-wide barriers to gene flow between closely related taxa.

Evaluating evidence of mitonuclear incompatibilities with the sex chromosomes in an avian hybrid zone

The exploration of hybrid zones and the intergenomic conflicts exposed through hybridization provide windows into the processes of divergence and speciation. Sex chromosomes and mitonuclear incompatibilities have strong associations with the genetics of hybrid dysfunction. In ZW sex-determining systems, maternal co-inheritance of the mitochondrial and W chromosomes immediately exposes incompatibilities between these maternal contributions of one species and the Z chromosome of another. We analyze mitochondrial and Z chromosome admixture in the long-tailed finch (Poephila acuticauda) of Australia, where hybridizing subspecies differ prominently in Z chromosome genotype and in bill color, yet the respective centers of geographic admixture for these two traits are offset by 350 km. We report two well-defined mitochondrial clades that diverged ∼0.5 million years ago. Mitochondrial contact is geographically co-located within a hybrid zone of Z chromosome admixture and is displaced from bill color admixture by nearly 400 km. Consistent with Haldane's rule expectations, hybrid zone females are significantly less likely than males to carry an admixed Z chromosome or have mismatched Z-mitochondrial genotypes. Furthermore, there are significantly fewer than expected mitonuclear mismatches in hybrid zone females and paternal backcross males. Results suggest a potential for mitonuclear/sex chromosome incompatibilities in the emergence of reproductive isolation in this system.

The Role of Mitonuclear Incompatibility in Bipolar Disorder Susceptibility and Resilience Against Environmental Stressors

It has been postulated that mitochondrial dysfunction has a significant role in the underlying pathophysiology of bipolar disorder (BD). Mitochondrial functioning plays an important role in regulating synaptic transmission, brain function, and cognition. Neuronal activity is energy dependent and neurons are particularly sensitive to changes in bioenergetic fluctuations, suggesting that mitochondria regulate fundamental aspects of brain function. Vigorous evidence supports the role of mitochondrial dysfunction in the etiology of BD, including dysregulated oxidative phosphorylation, general decrease of energy, altered brain bioenergetics, co-morbidity with mitochondrial disorders, and association with genetic variants in mitochondrial DNA (mtDNA) or nuclear-encoded mitochondrial genes. Despite these advances, the underlying etiology of mitochondrial dysfunction in BD is unclear. A plausible evolutionary explanation is that mitochondrial-nuclear (mitonuclear) incompatibility leads to a desynchronization of machinery required for efficient electron transport and cellular energy production. Approximately 1,200 genes, encoded from both nuclear and mitochondrial genomes, are essential for mitochondrial function. Studies suggest that mitochondrial and nuclear genomes co-evolve, and the coordinated expression of these interacting gene products are essential for optimal organism function. Incompatibilities between mtDNA and nuclear-encoded mitochondrial genes results in inefficiency in electron flow down the respiratory chain, differential oxidative phosphorylation efficiency, increased release of free radicals, altered intracellular Ca²⁺ signaling, and reduction of catalytic sites and ATP production. This review explores the role of mitonuclear incompatibility in BD susceptibility and resilience against environmental stressors.

A DNA barcode-based survey of wild urban bees in the Loire Valley, France

The current decline of wild bees puts important ecosystem services such as pollination at risk. Both inventory and monitoring programs are needed to understand the causes of wild bee decline. Effective insect monitoring relies on both mass-trapping methods coupled with rapid and accurate identifications. Identifying wild bees using only morphology can be challenging, in particular, specimens from mass-trapped samples which are often in poor condition. We generated DNA barcodes for 2931 specimens representing 157 species (156 named and one unnamed species) and 28 genera. Automated cluster delineation reveals 172 BINs (Barcodes Index Numbers). A total of 36 species (22.93%) were found in highly urbanized areas. The majority of specimens, representing 96.17% of the species barcoded form reciprocally exclusive groups, allowing their unambiguous identification. This includes several closely related species notoriously difficult to identify. A total of 137 species (87.26%) show a "one-to-one" match between a named species and the BIN assignment. Fourteen species (8.92%) show deep conspecific lineages with no apparent morphological differentiation. Only two species pairs shared the same BIN making their identification with DNA barcodes alone uncertain. Therefore, our DNA barcoding reference library allows reliable identification by non-experts for the vast majority of wild bee species in the Loire Valley.

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.

Integrating Mitochondrial Aerobic Metabolism into Ecology and Evolution

Biologists have long appreciated the critical role that energy turnover plays in understanding variation in performance and fitness among individuals. Whole-organism metabolic studies have provided key insights into fundamental ecological and evolutionary processes. However, constraints operating at subcellular levels, such as those operating within the mitochondria, can also play important roles in optimizing metabolism over different energetic demands and time scales. Herein, we explore how mitochondrial aerobic metabolism influences different aspects of organismal performance, such as through changing adenosine triphosphate (ATP) and reactive oxygen species (ROS) production. We consider how such insights have advanced our understanding of the mechanisms underpinning key ecological and evolutionary processes, from variation in life-history traits to adaptation to changing thermal conditions, and we highlight key areas for future research.

Diversification, Introgression, and Rampant Cytonuclear Discordance in Rocky Mountains Chipmunks (Sciuridae: Tamias)

Evidence from natural systems suggests that hybridization between animal species is more common than traditionally thought, but the overall contribution of introgression to standing genetic variation within species remains unclear for most animal systems. Here, we use targeted exon-capture to sequence thousands of nuclear loci and complete mitochondrial genomes from closely related chipmunk species in the Tamias quadrivittatus group that are distributed across the Great Basin and the central and southern Rocky Mountains of North America. This recent radiation includes six overlapping, ecologically distinct species (T. canipes, T. cinereicollis, T. dorsalis, T. quadrivittatus, T. rufus, and T. umbrinus) that show evidence for widespread introgression across species boundaries. Such evidence has historically been derived from a handful of markers, typically focused on mitochondrial loci, to describe patterns of introgression; consequently, the extent of introgression of nuclear genes is less well characterized. We conducted a series of phylogenomic and species-tree analyses to resolve the phylogeny of six species in this group. In addition, we performed several population genomic analyses to characterize nuclear genomes and infer coancestry among individuals. Furthermore, we used emerging quartets-based approaches to simultaneously infer the species tree (SVDquartets) and identify introgression (HyDe). We found that, in spite of rampant introgression of mitochondrial genomes between some species pairs (and sometimes involving up to three species), there appears to be little to no evidence for nuclear introgression. These findings mirror other genomic results where complete mitochondrial capture has occurred between chipmunk species in the absence of appreciable nuclear gene flow. The underlying causes of recurrent massive cytonuclear discordance remain unresolved in this group but mitochondrial DNA appears highly misleading of population histories as a whole. Collectively, it appears that chipmunk species boundaries are largely impermeable to nuclear gene flow and that hybridization, while pervasive with respect to mtDNA, has likely played a relatively minor role in the evolutionary history of this group.

Mitogenomics of the endemic Ethiopian rats: looking for footprints of adaptive evolution in sky islands

Organisms living in high altitude must adapt to environmental conditions with hypoxia and low temperature, e.g. by changes in the structure and function of proteins associated with oxidative phosphorylation in mitochondria. Here we analysed the signs of adaptive evolution in 27 mitogenomes of endemic Ethiopian rats (Stenocephalemys), where individual species adapted to different elevation. Significant signals of positive selection were detected in 10 of the 13 mitochondrial protein-coding genes, with a majority of functional substitutions in the NADH dehydrogenase complex. Higher frequency of positively selected sites was found in phylogenetic lineages corresponding to Afroalpine specialists.

Towards understanding the evolutionary dynamics of mtDNA

Historically, mtDNA was considered a selectively neutral marker that was useful for estimating the population genetic history of the maternal lineage. Over time there has been an increasing appreciation of mtDNA and mitochondria in maintaining cellular and organismal health. Beyond energy production, mtDNA and mitochondria have critical cellular roles in signalling. Here we briefly review the structure of mtDNA and the role of the mitochondrion in energy production. We then discuss the predictions that can be obtained from quaternary structure modelling and focus on mitochondrial complex I. Complex I is the primary entry point for electrons into the electron transport system is the largest respiratory complex of the chain and produces about 40% of the proton flux used to synthesize ATP. A focus of the review is Drosophila's utility as a model organism to study the selective advantage of specific mutations. However, we note that the incorporation of insights from a multitude of systems is necessary to fully understand the range of roles that mtDNA has in organismal fitness. We speculate that dietary changes can illicit stress responses that influence the selective advantage of specific mtDNA mutations and cause spatial and temporal fluctuations in the frequencies of mutations. We conclude that developing our understanding of the roles mtDNA has in determining organismal fitness will enable increased evolutionary insight and propose we can no longer assume it is evolving as a strictly neutral marker without testing this hypothesis.

Mitochondrial-nuclear coadaptation revealed through mtDNA replacements in Saccharomyces cerevisiae

BACKGROUND

Mitochondrial function requires numerous genetic interactions between mitochondrial- and nuclear- encoded genes. While selection for optimal mitonuclear interactions should result in coevolution between both genomes, evidence for mitonuclear coadaptation is challenging to document. Genetic models where mitonuclear interactions can be explored are needed.

RESULTS

We systematically exchanged mtDNAs between 15 Saccharomyces cerevisiae isolates from a variety of ecological niches to create 225 unique mitochondrial-nuclear genotypes. Analysis of phenotypic profiles confirmed that environmentally-sensitive interactions between mitochondrial and nuclear genotype contributed to growth differences. Exchanges of mtDNAs between strains of the same or different clades were just as likely to demonstrate mitonuclear epistasis although epistatic effect sizes increased with genetic distances. Strains with their original mtDNAs were more fit than strains with synthetic mitonuclear combinations when grown in media that resembled isolation habitats.

CONCLUSIONS

This study shows that natural variation in mitonuclear interactions contributes to fitness landscapes. Multiple examples of coadapted mitochondrial-nuclear genotypes suggest that selection for mitonuclear interactions may play a role in helping yeasts adapt to novel environments and promote coevolution.

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.

Genetic hitchhiking, mitonuclear coadaptation, and the origins of mt DNA barcode gaps

DNA barcoding based on mitochondrial (mt) nucleotide sequences is an enigma. Neutral models of mt evolution predict DNA barcoding cannot work for recently diverged taxa, and yet, mt DNA barcoding accurately delimits species for many bilaterian animals. Meanwhile, mt DNA barcoding often fails for plants and fungi. I propose that because mt gene products must cofunction with nuclear gene products, the evolution of mt genomes is best understood with full consideration of the two environments that impose selective pressure on mt genes: the external environment and the internal genomic environment. Moreover, it is critical to fully consider the potential for adaptive evolution of not just protein products of mt genes but also of mt transfer RNAs and mt ribosomal RNAs. The tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. Accordingly, the most important factor determining whether or not mt DNA sequences diagnose species boundaries may be the extent to which the mt chromosomes engage in recombination.

A mitochondrial genetic divergence proxy predicts the reproductive compatibility of mammalian hybrids

Numerous pairs of evolutionarily divergent mammalian species have been shown to produce hybrid offspring. In some cases, F1 hybrids are able to produce F2s through matings with F1s. In other instances, the hybrids are only able to produce offspring themselves through backcrosses with a parent species owing to unisexual sterility (Haldane's Rule). Here, we explicitly tested whether genetic distance, computed from mitochondrial and nuclear genes, can be used as a proxy to predict the relative fertility of the hybrid offspring resulting from matings between species of terrestrial mammals. We assessed the proxy's predictive power using a well-characterized felid hybrid system, and applied it to modern and ancient hominins. Our results revealed a small overlap in mitochondrial genetic distance values that distinguish species pairs whose calculated distances fall within two categories: those whose hybrid offspring follow Haldane's Rule, and those whose hybrid F1 offspring can produce F2s. The strong correlation between genetic distance and hybrid fertility demonstrated here suggests that this proxy can be employed to predict whether the hybrid offspring of two mammalian species will follow Haldane's Rule.

Modeling Mito-nuclear Compatibility and Its Role in Species Identification

Mitochondrial genetic material (mtDNA) is widely used for phylogenetic reconstruction and as a barcode for species identification. The utility of mtDNA in these contexts derives from its particular molecular properties, including its high evolutionary rate, uniparental inheritance, and small size. But mtDNA may also play a fundamental role in speciation-as suggested by recent observations of coevolution with the nuclear DNA, along with the fact that respiration depends on coordination of genes from both sources. Here, we study how mito-nuclear interactions affect the accuracy of species identification by mtDNA, as well as the speciation process itself. We simulate the evolution of a population of individuals who carry a recombining nuclear genome and a mitochondrial genome inherited maternally. We compare a null model fitness landscape that lacks any mito-nuclear interaction against a scenario in which interactions influence fitness. Fitness is assigned to individuals according to their mito-nuclear compatibility, which drives the coevolution of the nuclear and mitochondrial genomes. Depending on the model parameters, the population breaks into distinct species and the model output then allows us to analyze the accuracy of mtDNA barcode for species identification. Remarkably, we find that species identification by mtDNA is equally accurate in the presence or absence of mito-nuclear coupling and that the success of the DNA barcode derives mainly from population geographical isolation during speciation. Nevertheless, selection imposed by mito-nuclear compatibility influences the diversification process and leaves signatures in the genetic content and spatial distribution of the populations, in three ways. First, speciation is delayed and the resulting phylogenetic trees are more balanced. Second, clades in the resulting phylogenetic tree correlate more strongly with the spatial distribution of species and clusters of more similar mtDNA's. Third, there is a substantial increase in the intraspecies mtDNA similarity, decreasing the number of alleles substitutions per locus and promoting the conservation of genetic information. We compare the evolutionary patterns observed in our model to empirical data from copepods (Tigriopus californicus). We find good qualitative agreement in the geographic patterns and the topology of the phylogenetic tree, provided the model includes selection based on mito-nuclear interactions. These results highlight the role of mito-nuclear compatibility in the speciation process and its reconstruction from genetic data.[Mito-nuclear coevolution; mtDNA barcode; parapatry; phylogeny.].

Phylogenetic relationships between genera Dorcadion, Lamia, Morimus, Herophila and some other Lamiinae (Coleoptera: Cerambycidae) based on chromosome and CO1 gene sequence comparison

A dual molecular and cytogenetic study was performed with the aim to improve the controversial systematic classification of some species of Lamiinae (Coleoptera: Cerambycidae). The karyotypes of species belonging to genera Morimus, Herophila, Dorcadion, Neodorcadion and Lamia share a number of characters, which differentiate them from other species, belonging to genera Phytoecia, Parmena and Monochamus. The karyotypes of the last three species comprise 20 chromosomes, mostly metacentric or sub-metacentric, as in the presumed Cerambycidae ancestors. The karyotypes of the former species share many characters derived from the Lamiinae ancestors by a number of chromosome fissions and inversions indicating their monophyly. Comparisons of the CO1 gene sequence also show the monophyly of Morimus, Lamia, Herophila and Dorcadion and their distant relationship with others. These convergent results allow us to propose a phylogenetic classification of these genera, which places the monospecific genus Lamia close to Dorcadion, clearly separates Dorcadion and Neodorcadion and places Herophila closer to Morimus than to Dorcadion/Lamia. The genus Morimus is the most derived. CO1 mutations loosely separate the forms M. asper and M. funereus, which have similar karyotypes and behaviour and copulate in captivity. The form M. ganglebaueri may have a funereus X asper hybrid origin.

Mitonuclear mismatch alters performance and reproductive success in naturally introgressed populations of a montane leaf beetle

Coordination between nuclear and mitochondrial genomes is critical to metabolic processes underlying animals' ability to adapt to local environments, yet consequences of mitonuclear interactions have rarely been investigated in populations where individuals with divergent mitochondrial and nuclear genomes naturally interbreed. Genetic variation in the leaf beetle Chrysomela aeneicollis was assessed along a latitudinal thermal gradient in California's Sierra Nevada. Variation at mitochondrial cytochrome oxidase II (COII) and the nuclear gene phosphoglucose isomerase (PGI) shows concordance and was significantly greater along a 65 km transect than 10 other loci. STRUCTURE analyses using neutral loci identified a southern and northern subpopulation, which interbreed in the central drainage Bishop Creek. COII and PGI were used as indicators of mitochondrial and nuclear genetic variation in field and laboratory experiments conducted on beetles from this admixed population. Fecundity, larval development rate, running speed and male mating frequency were higher for beetles with geographically "matched" than "mismatched" mitonuclear genotypes. Effects of mitonuclear mismatch were largest for individuals with northern nuclear genotypes possessing southern mitochondria and were most pronounced after heat treatment or at high elevation. These findings suggest that mitonuclear incompatibility diminishes performance and reproductive success in nature, effects that could intensify at environmental extremes.

Mitonuclear Compensatory Coevolution

In bilaterian animals, the mitochondrial genome is small, haploid, does not typically recombine, and is subject to accumulation of deleterious alleles via Muller's ratchet. These basic features of the genomic architecture present a paradox: mutational erosion of these genomes should lead to decline in mitochondrial function over time, yet no such decline is observed. Compensatory coevolution, whereby the nuclear genome evolves to compensate for the deleterious alleles in the mitochondrial genome, presents a potential solution to the paradox of Muller's ratchet without loss of function. Here, I review different proposed forms of mitonuclear compensatory coevolution. Empirical evidence from diverse eukaryotic taxa supports the mitonuclear compensatory coevolution hypothesis, but the ubiquity and importance of such compensatory coevolution remains a topic of debate.

Cytonuclear Genetic Incompatibilities in Plant Speciation

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

Strong selective effects of mitochondrial DNA on the nuclear genome

Oxidative phosphorylation, the primary source of cellular energy in eukaryotes, requires gene products encoded in both the nuclear and mitochondrial genomes. As a result, functional integration between the genomes is essential for efficient adenosine triphosphate (ATP) generation. Although within populations this integration is presumably maintained by coevolution, the importance of mitonuclear coevolution in key biological processes such as speciation and mitochondrial disease has been questioned. In this study, we crossed populations of the intertidal copepod Tigriopus californicus to disrupt putatively coevolved mitonuclear genotypes in reciprocal F₂ hybrids. We utilized interindividual variation in developmental rate among these hybrids as a proxy for fitness to assess the strength of selection imposed on the nuclear genome by alternate mitochondrial genotypes. Developmental rate varied among hybrid individuals, and in vitro ATP synthesis rates of mitochondria isolated from high-fitness hybrids were approximately two-fold greater than those of mitochondria isolated from low-fitness individuals. We then used Pool-seq to compare nuclear allele frequencies for high- or low-fitness hybrids. Significant biases for maternal alleles were detected on 5 (of 12) chromosomes in high-fitness individuals of both reciprocal crosses, whereas maternal biases were largely absent in low-fitness individuals. Therefore, the most fit hybrids were those with nuclear alleles that matched their mitochondrial genotype on these chromosomes, suggesting that mitonuclear effects underlie individual-level variation in developmental rate and that intergenomic compatibility is critical for high fitness. We conclude that mitonuclear interactions can have profound impacts on both physiological performance and the evolutionary trajectory of the nuclear genome.

Interactions between cytoplasmic and nuclear genomes confer sex-specific effects on lifespan in Drosophila melanogaster

Genetic variation outside of the cell nucleus can affect the phenotype. The cytoplasm is home to the mitochondria, and in arthropods often hosts intracellular bacteria such as Wolbachia. Although numerous studies have implicated epistatic interactions between cytoplasmic and nuclear genetic variation as mediators of phenotypic expression, two questions remain. Firstly, it remains unclear whether outcomes of cyto-nuclear interactions will manifest differently across the sexes, as might be predicted given that cytoplasmic genomes are screened by natural selection only through females as a consequence of their maternal inheritance. Secondly, the relative contribution of mitochondrial genetic variation to other cytoplasmic sources of variation, such as Wolbachia infection, in shaping phenotypic outcomes of cyto-nuclear interactions remains unknown. Here, we address these questions, creating a fully crossed set of replicated cyto-nuclear populations derived from three geographically distinct populations of Drosophila melanogaster, measuring the lifespan of males and females from each population. We observed that cyto-nuclear interactions shape lifespan and that the outcomes of these interactions differ across the sexes. Yet, we found no evidence that placing the cytoplasms from one population alongside the nuclear background of others (generating putative cyto-nuclear mismatches) leads to decreased lifespan in either sex. Although it was difficult to partition mitochondrial from Wolbachia effects, our results suggest at least some of the cytoplasmic genotypic contribution to lifespan was directly mediated by an effect of sequence variation in the mtDNA. Future work should explore the degree to which cyto-nuclear interactions result in sex differences in the expression of other components of organismal life history.

Distinct patterns of hybridization across a suture zone in a coral reef fish (Dascyllus trimaculatus)

Hybrid zones are natural laboratories for investigating the dynamics of gene flow, reproductive isolation, and speciation. A predominant marine hybrid (or suture) zone encompasses Christmas Island (CHR) and Cocos (Keeling) Islands (CKE), where 15 different instances of interbreeding between closely related species from Indian and Pacific Oceans have been documented. Here, we report a case of hybridization between genetically differentiated Pacific and Indian Ocean lineages of the three-spot dascyllus, Dascyllus trimaculatus (Rüppell, 1829). Field observations indicate there are subtle color differences between Pacific and Indian Ocean lineages. Most importantly, population densities of color morphs and genetic analyses (mitochondrial DNA and SNPs obtained via RADSeq) suggest that the pattern of hybridization within the suture zone is not homogeneous. At CHR, both color morphs were present, mitochondrial haplotypes of both lineages were observed, and SNP analyses revealed both pure and hybrid genotypes. Meanwhile, in CKE, the Indian Ocean color morphs were prevalent, only Indian Ocean mitochondrial haplotypes were observed, and SNP analysis showed hybrid individuals with a large proportion (~80%) of their genotypes assigning to the Indian Ocean lineage. We conclude that CHR populations are currently receiving an influx of individuals from both ocean basins, with a greater influence from the Pacific Ocean. In contrast, geographically isolated CKE populations appear to be self-recruiting and with more influx of individuals from the Indian Ocean. Our research highlights how patterns of hybridization can be different at scales of hundreds of kilometers, due to geographic isolation and the history of interbreeding between lineages.

First genetic evidence that invasive bullhead (Cottus L. 1758) in Scotland is of English origin and the difficulty of resolving the European Cottus species taxonomy

The European bullhead (Cottus gobio) is widely distributed across Europe, and within the UK is native to England and Wales, where it is protected under the Habitats Directive. In Scotland, however, the species is considered invasive and thriving populations are recorded in the Forth and Clyde river catchments, and the Ale Water in the Scottish Borders. The genetic identity of the Scottish populations has not been established. There is also debate about the status of the European bullhead and its validity as single species, a species complex with several unresolved species, or distinct different species in its European distribution range. There is therefore a need to determine the taxonomy and likely source of the novel Scottish populations. Genetic analyses using cytochrome oxidase 1 (COI) mitochondrial DNA sequences were undertaken on specimens from the Forth and Clyde catchments, and combined with the results of morphological characteristics to provide a comprehensive assessment of the taxonomic classification for Scottish bullheads. There was considerable variation in morphological characteristics between populations within Scotland and a wider range of variability than previously recorded for English populations. Genetically the Scottish populations were very closely related to English specimens, supporting the hypothesis of introduction directly from England to Scotland. In terms of broader relationships, Scottish specimens are genetically more closely related to the ostensible species Chabot fluviatile Cottus perifretum, which has been suggested as one of a complex of species across Europe. Morphologically they exhibit characteristics on the spectrum between C. perifretum and C. gobio. There is an urgent need for the clarification of the taxonomy of Cottus sp(p). to avoid confusion in future publications, legislation and management practices relating to bullheads throughout the UK and Europe.

Beyond the Powerhouse: Integrating Mitonuclear Evolution, Physiology, and Theory in Comparative Biology

Eukaryotes are the outcome of an ancient symbiosis and as such, eukaryotic cells fundamentally possess two genomes. As a consequence, gene products encoded by both nuclear and mitochondrial genomes must interact in an intimate and precise fashion to enable aerobic respiration in eukaryotes. This genomic architecture of eukaryotes is proposed to necessitate perpetual coevolution between the nuclear and mitochondrial genomes to maintain coadaptation, but the presence of two genomes also creates the opportunity for intracellular conflict. In the collection of papers that constitute this symposium volume, scientists working in diverse organismal systems spanning vast biological scales address emerging topics in integrative, comparative biology in light of mitonuclear interactions.

Sexually Antagonistic Mitonuclear Coevolution in Duplicate Oxidative Phosphorylation Genes

Mitochondrial function is critical in eukaryotes. To maintain an adequate supply of energy, precise interactions must be maintained between nuclear- and mitochondrial-encoded gene products. Such interactions are paramount in chimeric enzymes such as the oxidative phosphorylation (OXPHOS) complexes. Mutualistic coevolution between the two genomes has therefore been suggested to be a critical, ubiquitous feature of eukaryotes that acts to maintain cellular function. However, mitochondrial genomes can also act selfishly and increase their own transmission at the expense of organismal function. For example, male-harming mutations are predisposed to accumulate in mitochondrial genomes due to their maternal inheritance ("mother's curse"). Here, we investigate sexually antagonistic mitonuclear coevolution in nuclear-encoded OXPHOS paralogs from mammals and Drosophila. These duplicate genes are highly divergent but must interact with the same set of mitochondrial-encoded genes. Many such paralogs show testis-specific expression, prompting previous hypotheses suggesting they may have evolved under selection to counteract male-harming mitochondrial mutations. We found increased rates of evolution in OXPHOS paralogs with testis-specific expression in mammals and Drosophila, supporting this hypothesis. However, further analyses suggested such patterns may be due to relaxed, not positive selection, especially in Drosophila. Structural data also suggest that mitonuclear interactions do not play a major role in the evolution of many OXPHOS paralogs in a consistent way. In conclusion, no single OXPHOS paralog met all our criteria for being under selection to counteract male-harming mitochondrial mutations. We discuss alternative explanations for the drastic patterns of evolution in these genes, including mutualistic mitonuclear coevolution, adaptive subfunctionalization after gene duplication, and relaxed selection on OXPHOS in male tissues.