Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout ageing

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mitochondrial genome in different tissues throughout ageing. We used ultrasensitive duplex sequencing to profile ~2.5 million mitochondrial genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloguing ~1.2 million mitochondrial somatic and ultralow-frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the light strand origin of replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared with primates with a surfeit of reactive oxygen species-associated G > T/C > A mutations, and that somatic mutations in protein-coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that 're-align' mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

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.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout aging

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here, we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mt-genome in different tissues throughout aging. We used ultra-sensitive Duplex Sequencing to profile ~2.5 million mt-genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloging ~1.2 million mitochondrial somatic and ultra low frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the Light Strand Origin of Replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared to primates with a surfeit of reactive oxygen species-associated G>T/C>A mutations, and that somatic mutations in protein coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that "re-align" mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

A Test of Haldane's Rule in Neodiprion Sawflies and Implications for the Evolution of Postzygotic Isolation in Haplodiploids

AbstractHaldane's rule-a pattern in which hybrid sterility or inviability is observed in the heterogametic sex of an interspecific cross-is one of the most widely obeyed rules in nature. Because inheritance patterns are similar for sex chromosomes and haplodiploid genomes, Haldane's rule may apply to haplodiploid taxa, predicting that haploid male hybrids will evolve sterility or inviability before diploid female hybrids. However, there are several genetic and evolutionary mechanisms that may reduce the tendency of haplodiploids to obey Haldane's rule. Currently, there are insufficient data from haplodiploids to determine how frequently they adhere to Haldane's rule. To help fill this gap, we crossed a pair of haplodiploid hymenopteran species (Neodiprion lecontei and Neodiprion pinetum) and evaluated the viability and fertility of female and male hybrids. Despite considerable divergence, we found no evidence of reduced fertility in hybrids of either sex, consistent with the hypothesis that hybrid sterility evolves slowly in haplodiploids. For viability, we found a pattern opposite to that of Haldane's rule: hybrid females, but not males, had reduced viability. This reduction was most pronounced in one direction of the cross, possibly due to a cytoplasmic-nuclear incompatibility. We also found evidence of extrinsic postzygotic isolation in hybrids of both sexes, raising the possibility that this form or reproductive isolation tends to emerge early in speciation in host-specialized insects. Our work emphasizes the need for more studies on reproductive isolation in haplodiploids, which are abundant in nature but underrepresented in the speciation literature.

Ecology and the origin of non-ephemeral species

AbstractResearch over the past three decades has shown that ecology-based extrinsic reproductive barriers can rapidly arise to generate incipient species-but such barriers can also rapidly dissolve when environments change, resulting in incipient species collapse. Understanding the evolution of unconditional, "intrinsic" reproductive barriers is therefore important for understanding the longer-term buildup of biodiversity. In this article, we consider ecology's role in the evolution of intrinsic reproductive isolation. We suggest that this topic has fallen into a gap between disciplines: while evolutionary ecologists have traditionally focused on the rapid evolution of extrinsic isolation between co-occurring ecotypes, speciation geneticists studying intrinsic isolation in other taxa have devoted little attention to the ecological context in which it evolves. We argue that for evolutionary ecology to close this gap, the field will have to expand its focus beyond rapid adaptation and its traditional model systems. Synthesizing data from several subfields, we present circumstantial evidence for and against different forms of ecological adaptation as promoters of intrinsic isolation and discuss alternative forces that may be significant. We conclude by outlining complementary approaches that can better address the role of ecology in the evolution of nonephemeral reproductive barriers and, by extension, less ephemeral species.

Patterns of Performance Variation Between Animal Hybrids and their Parents: A Meta-analysis

Hybridization is a widespread phenomenon in animals, and hybrid heterosis/breakdown could be key processes determining the evolutionary dynamics of hybrids. Indeed, hybrids are not consistently disadvantaged compared to the parental lineages, as was historically assumed. Multiple processes could lead to performance differences between parental lineages and their hybrids. Despite many studies evaluated the performance of hybrids, a quantitative synthesis is required to assess the general pattern. Here we used meta-analytic and meta-regression approaches to quantify the fitness differences between parental lineages and their hybrids, and to identify possible processes that could lead to these differences. Specifically, we tested biological and methodological parameters that could determine differences in performance between hybrids and parental lineages. Hybrid performance was extremely variable across studies, being often significantly higher or lower compared to the mean performance of their parents. Nevertheless, the averaged hybrid performance was similar to the fitness of parental lineages, with differences across studies related to how performance was assessed. Genetic divergence between parental lineages, and the approach used to identify hybrids were the parameters most strongly related to variation in hybrid performance. Performance was lower for hybrids between distantly related lineages. Furthermore, study settings and the use of imprecise approaches for hybrid identification (e.g. morphology-based) can bias assessments of performance. Studies performed on wild populations and using genetic approaches for hybrid identification detected more often a decreased hybrid performance, compared to laboratory studies. We highlight the importance of appropriate settings for a realistic understanding of the evolutionary impacts of hybridization.

Genome-wide local ancestry and evidence for mitonuclear coadaptation in African hybrid cattle populations (Bos taurus/indicus)

The phenotypic diversity of African cattle reflects adaptation to a wide range of agroecological conditions, human-mediated selection preferences, and complex patterns of admixture between the humpless Bos taurus (taurine) and humped Bos indicus (zebu) subspecies, which diverged 150-500 thousand years ago. Despite extensive admixture, all African cattle possess taurine mitochondrial haplotypes, even populations with significant zebu biparental and male uniparental nuclear ancestry. This has been interpreted as the result of human-mediated dispersal ultimately stemming from zebu bulls imported from South Asia during the last three millennia. Here, we assess whether ancestry at mitochondrially targeted nuclear genes in African admixed cattle is impacted by mitonuclear functional interactions. Using high-density SNP data, we find evidence for mitonuclear coevolution across hybrid African cattle populations with a significant increase of taurine ancestry at mitochondrially targeted nuclear genes. Our results, therefore, support the hypothesis of incompatibility between the taurine mitochondrial genome and the zebu nuclear genome.

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.

The emergence of ecotypes in a parasitoid wasp: a case of incipient sympatric speciation in Hymenoptera?

Background

To understand which reproductive barriers initiate speciation is a major question in evolutionary research. Despite their high species numbers and specific biology, there are only few studies on speciation in Hymenoptera. This study aims to identify very early reproductive barriers in a local, sympatric population of Nasonia vitripennis (Walker 1836), a hymenopterous parasitoid of fly pupae. We studied ecological barriers, sexual barriers, and the reduction in F1-female offspring as a postmating barrier, as well as the population structure using microsatellites.

Results

We found considerable inbreeding within female strains and a population structure with either three or five subpopulation clusters defined by microsatellites. In addition, there are two ecotypes, one parasitizing fly pupae in bird nests and the other on carrion. The nest ecotype is mainly formed from one of the microsatellite clusters, the two or four remaining microsatellite clusters form the carrion ecotype. There was slight sexual isolation and a reduction in F1-female offspring between inbreeding strains from the same microsatellite clusters and the same ecotypes. Strains from different microsatellite clusters are separated by a reduction in F1-female offspring. Ecotypes are separated only by ecological barriers.

Conclusions

This is the first demonstration of very early reproductive barriers within a sympatric population of Hymenoptera. It demonstrates that sexual and premating barriers can precede ecological separation. This indicates the complexity of ecotype formation and highlights the general need for more studies within homogenous populations for the identification of the earliest barriers in the speciation process.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01938-y.

Mitonuclear interactions and introgression genomics of macaque monkeys (Macaca) highlight the influence of behaviour on genome evolution

In most macaques, females are philopatric and males migrate from their natal ranges, which results in pronounced divergence of mitochondrial genomes within and among species. We therefore predicted that some nuclear genes would have to acquire compensatory mutations to preserve compatibility with diverged interaction partners from the mitochondria. We additionally expected that these sex-differences would have distinctive effects on gene flow in the X and autosomes. Using new genomic data from 29 individuals from eight species of Southeast Asian macaque, we identified evidence of natural selection associated with mitonuclear interactions, including extreme outliers of interspecies differentiation and metrics of positive selection, low intraspecies polymorphism and atypically long runs of homozygosity associated with nuclear-encoded genes that interact with mitochondria-encoded genes. In one individual with introgressed mitochondria, we detected a small but significant enrichment of autosomal introgression blocks from the source species of her mitochondria that contained genes which interact with mitochondria-encoded loci. Our analyses also demonstrate that sex-specific demography sculpts genetic exchange across multiple species boundaries. These findings show that behaviour can have profound but indirect effects on genome evolution by influencing how interacting components of different genomic compartments (mitochondria, the autosomes and the sex chromosomes) move through time and space.

Restoring fertility in yeast hybrids: Breeding and quantitative genetics of beneficial traits

Hybrids between species can harbor a combination of beneficial traits from each parent and may exhibit hybrid vigor, more readily adapting to new harsher environments. Interspecies hybrids are also sterile and therefore an evolutionary dead end unless fertility is restored, usually via auto-polyploidisation events. In the Saccharomyces genus, hybrids are readily found in nature and in industrial settings, where they have adapted to severe fermentative conditions. Due to their hybrid sterility, the development of new commercial yeast strains has so far been primarily conducted via selection methods rather than via further breeding. In this study, we overcame infertility by creating tetraploid intermediates of Saccharomyces interspecies hybrids to allow continuous multigenerational breeding. We incorporated nuclear and mitochondrial genetic diversity within each parental species, allowing for quantitative genetic analysis of traits exhibited by the hybrids and for nuclear-mitochondrial interactions to be assessed. Using pooled F12 generation segregants of different hybrids with extreme phenotype distributions, we identified quantitative trait loci (QTLs) for tolerance to high and low temperatures, high sugar concentration, high ethanol concentration, and acetic acid levels. We identified QTLs that are species specific, that are shared between species, as well as hybrid specific, in which the variants do not exhibit phenotypic differences in the original parental species. Moreover, we could distinguish between mitochondria-type-dependent and -independent traits. This study tackles the complexity of the genetic interactions and traits in hybrid species, bringing hybrids into the realm of full genetic analysis of diploid species, and paves the road for the biotechnological exploitation of yeast biodiversity.

Improving mosquito control strategies with population genomics

Mosquito control strategies increasingly apply knowledge from population genomics research. This review highlights recent applications to three research domains: mosquito invasions, insecticide resistance evolution, and rear and release programs. Current research trends follow developments in reference assemblies, either as improvements to existing assemblies (particularly Aedes) or assemblies for new taxa (particularly Anopheles). With improved assemblies, studies of invasive and rear and release target populations are better able to incorporate adaptive as well as demographic hypotheses. New reference assemblies are aiding comparisons of insecticide resistance across sister taxa while helping resolve taxon boundaries amidst frequent introgression. Anopheles gene drive deployments and improved Aedes genome assemblies should lead to a convergence in research aims for Anopheles and Aedes in the coming years.

Maternal reproductive output and F1 hybrid fitness may influence contact zone dynamics

The outcome of secondary contact between divergent lineages or species may be influenced by both the reproductive traits of parental species and the fitness of offspring; however, their relative contributions have rarely been evaluated, particularly in longer-lived vertebrate species. We performed pure and reciprocal laboratory crosses between Ctenophorus decresii (tawny dragon) and C. modestus (swift dragon) to examine how parental reproductive traits and ecologically relevant offspring fitness traits may explain contact zone dynamics in the wild. The two species meet in a contact zone of post-F1 hybrids with asymmetric backcrossing and predominantly C. modestus mtDNA haplotypes. We found no evidence for reduced parental fecundity or offspring fitness for F1 hybrid crosses. However, maternal reproductive strategy differed between species, irrespective of the species of their mate. Ctenophorus modestus females had higher fecundity and produced more and larger clutches with lower embryonic mortality. Parental species also influenced sex ratios and offspring traits, with C. modestus ♀ × C. decresii ♂ hybrids exhibiting higher trait values for more fitness measures (growth rate, sprint speed, bite force) than offspring from all other pairings. Together, these patterns are consistent with the prevalence of C. modestus mtDNA in the contact zone, and asymmetric backcrossing likely reflects fitness effects that manifest in the F2 generation. Our results highlight how parental species can influence multiple offspring traits in different ways, which together may combine to influence offspring fitness and shape contact zone dynamics.

Genomic divergence landscape in recurrently hybridizing Chironomus sister taxa suggests stable steady state between mutual gene flow and isolation

Divergence is mostly viewed as a progressive process often initiated by selection targeting individual loci, ultimately resulting in ever increasing genomic isolation due to linkage. However, recent studies show that this process may stall at intermediate stable equilibrium states without achieving complete genomic isolation. We tested the extent of genomic isolation between two recurrently hybridizing nonbiting midge sister taxa, Chironomus riparius and Chironomus piger, by analyzing the divergence landscape. Using a principal component-based method, we estimated that only about 28.44% of the genomes were mutually isolated, whereas the rest was still exchanged. The divergence landscape was fragmented into isolated regions of on average 30 kb, distributed throughout the genome. Selection and divergence time strongly influenced lengths of isolated regions, whereas local recombination rate only had minor impact. Comparison of divergence time distributions obtained from several coalescence-simulated divergence scenarios with the observed divergence time estimates in an approximate Bayesian computation framework favored a short and concluded divergence event in the past. Most divergence happened during a short time span about 4.5 million generations ago, followed by a stable equilibrium between mutual gene flow through ongoing hybridization for the larger part of the genome and isolation in some regions due to rapid purifying selection of introgression, supported by high effective population sizes and recombination rates.

Pervasive Mitonuclear Coadaptation Underlies Fast Development in Interpopulation Hybrids of a Marine Crustacean

Cellular energy production requires coordinated interactions between genetic components from the nuclear and mitochondrial genomes. This coordination results in coadaptation of interacting elements within populations. Interbreeding between divergent gene pools can disrupt coadapted loci and result in hybrid fitness breakdown. While specific incompatible loci have been detected in multiple eukaryotic taxa, the extent of the nuclear genome that is influenced by mitonuclear coadaptation is not clear in any species. Here, we used F₂ hybrids between two divergent populations of the copepod Tigriopus californicus to examine mitonuclear coadaptation across the nuclear genome. Using developmental rate as a measure of fitness, we found that fast-developing copepods had higher ATP synthesis capacity than slow developers, suggesting variation in developmental rates is at least partly associated with mitochondrial dysfunction. Using Pool-seq, we detected strong biases for maternal alleles across 7 (of 12) chromosomes in both reciprocal crosses in high-fitness hybrids, whereas low-fitness hybrids showed shifts toward the paternal population. Comparison with previous results on a different hybrid cross revealed largely different patterns of strong mitonuclear coadaptation associated with developmental rate. Our findings suggest that functional coadaptation between interacting nuclear and mitochondrial components is reflected in strong polygenic effects on this life-history phenotype, and reveal that molecular coadaptation follows independent evolutionary trajectories among isolated populations.

Germline mutagenesis of Nasonia vitripennis through ovarian delivery of CRISPR-Cas9 ribonucleoprotein

CRISPR/Cas9 gene editing is a powerful technology to study the genetics of rising model organisms, such as the jewel wasp Nasonia vitripennis. However, current methods involving embryonic microinjection of CRISPR reagents are challenging. Delivery of Cas9 ribonucleoprotein into female ovaries is an alternative that has only been explored in a small handful of insects, such as mosquitoes, whiteflies and beetles. Here, we developed a simple protocol for germline gene editing by injecting Cas9 ribonucleoprotein in adult N. vitripennis females using either ReMOT control (Receptor-Mediated Ovary Transduction of Cargo) or BAPC (Branched Amphiphilic Peptide Capsules) as ovary delivery methods. For ReMOT Control we used the Drosophila melanogaster-derived peptide 'P2C' fused to EGFP to visualize the ovary delivery, and fused to Cas9 protein for gene editing of the cinnabar gene using saponin as an endosomal escape reagent. For BAPC we optimized the concentrations of protein, sgRNA and the transfection reagent. We demonstrate delivery of protein cargo such as EGFP and Cas9 into developing oocytes via P2C peptide and BAPC. Additionally, somatic and germline gene editing were demonstrated. This approach will greatly facilitate CRISPR-applied genetic manipulation in this and other rising model organisms.

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.

Extreme heterogeneity of human mitochondrial DNA from organelles to populations

Contrary to the long-held view that most humans harbour only identical mitochondrial genomes, deep resequencing has uncovered unanticipated extreme genetic variation within mitochondrial DNA (mtDNA). Most, if not all, humans contain multiple mtDNA genotypes (heteroplasmy); specific patterns of variants accumulate in different tissues, including cancers, over time; and some variants are preferentially passed down or suppressed in the maternal germ line. These findings cast light on the origin and spread of mtDNA mutations at multiple scales, from the organelle to the human population, and challenge the conventional view that high percentages of a mutation are required before a new variant has functional consequences.

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.

The importance of intrinsic postzygotic barriers throughout the speciation process

Intrinsic postzygotic barriers can play an important and multifaceted role in speciation, but their contribution is often thought to be reserved to the final stages of the speciation process. Here, we review how intrinsic postzygotic barriers can contribute to speciation, and how this role may change through time. We outline three major contributions of intrinsic postzygotic barriers to speciation. (i) reduction of gene flow: intrinsic postzygotic barriers can effectively reduce gene exchange between sympatric species pairs. We discuss the factors that influence how effective incompatibilities are in limiting gene flow. (ii) early onset of species boundaries via rapid evolution: intrinsic postzygotic barriers can evolve between recently diverged populations or incipient species, thereby influencing speciation relatively early in the process. We discuss why the early origination of incompatibilities is expected under some biological models, and detail how other (and often less obvious) incompatibilities may also serve as important barriers early on in speciation. (iii) reinforcement: intrinsic postzygotic barriers can promote the evolution of subsequent reproductive isolation through processes such as reinforcement, even between relatively recently diverged species pairs. We incorporate classic and recent empirical and theoretical work to explore these three facets of intrinsic postzygotic barriers, and provide our thoughts on recent challenges and areas in the field in which progress can be made. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

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.

Mitochondria and the Origin of Species: Bridging Genetic and Ecological Perspectives on Speciation Processes

Mitochondria have been known to be involved in speciation through the generation of Dobzhansky-Muller incompatibilities, where functionally neutral co-evolution between mitochondrial and nuclear genomes can cause dysfunction when alleles are recombined in hybrids. We propose that adaptive mitochondrial divergence between populations can not only produce intrinsic (Dobzhansky-Muller) incompatibilities, but could also contribute to reproductive isolation through natural and sexual selection against migrants, post-mating prezygotic isolation, as well as by causing extrinsic reductions in hybrid fitness. We describe how these reproductive isolating barriers can potentially arise through adaptive divergence of mitochondrial function in the absence of mito-nuclear coevolution, a departure from more established views. While a role for mitochondria in the speciation process appears promising, we also highlight critical gaps of knowledge: (1) many systems with a potential for mitochondrially-mediated reproductive isolation lack crucial evidence directly linking reproductive isolation and mitochondrial function; (2) it often remains to be seen if mitochondrial barriers are a driver or a consequence of reproductive isolation; (3) the presence of substantial gene flow in the presence of mito-nuclear incompatibilities raises questions whether such incompatibilities are strong enough to drive speciation to completion; and (4) it remains to be tested how mitochondrial effects on reproductive isolation compare when multiple mechanisms of reproductive isolation coincide. We hope this perspective and the proposed research plans help to inform future studies of mitochondrial adaptation in a manner that links genotypic changes to phenotypic adaptations, fitness, and reproductive isolation in natural systems, helping to clarify the importance of mitochondria in the formation and maintenance of biological diversity.

Thermal tolerance and fish heart integrity: fatty acids profiles as predictors of species resilience

The cardiovascular system is a major limiting system in thermal adaptation, but the exact physiological mechanisms underlying responses to thermal stress are still not completely understood. Recent studies have uncovered the possible role of reactive oxygen species production rates of heart mitochondria in determining species' upper thermal limits. The present study examines the relationship between individual response to a thermal challenge test (CT_max), susceptibility to peroxidation of membrane lipids, heart fatty acid profiles and cardiac antioxidant enzyme activities in two salmonid species from different thermal habitats (Salvelinus alpinus, Salvelinus fontinalis) and their hybrids. The susceptibility to peroxidation of membranes in the heart was negatively correlated with individual thermal tolerance. The same relationship was found for arachidonic and eicosapentaenoic acid. Total H₂O₂ buffering activity of the heart muscle was higher for the group with high thermal resistance. These findings underline a potential general causative relationship between sensitivity to oxidative stress, specific fatty acids, antioxidant activity in the cardiac muscle and thermal tolerance in fish and likely other ectotherms. Heart fatty acid profile could be indicative of species resilience to global change, and more importantly the plasticity of this trait could predict the adaptability of fish species or populations to changes in environmental temperature.

Mitonuclear conflict and cooperation govern the integration of genotypes, phenotypes and environments

The mitonuclear genome is the most successful co-evolved mutualism in the history of life on Earth. The cross-talk between the mitochondrial and nuclear genomes has been shaped by conflict and cooperation for more than 1.5 billion years, yet this system has adapted to countless genomic reorganizations by each partner, and done so under changing environments that have placed dramatic biochemical and physiological pressures on evolving lineages. From putative anaerobic origins, mitochondria emerged as the defining aerobic organelle. During this transition, the two genomes resolved rules for sex determination and transmission that made uniparental inheritance the dominant, but not a universal pattern. Mitochondria are much more than energy-producing organelles and play crucial roles in nutrient and stress signalling that can alter how nuclear genes are expressed as phenotypes. All of these interactions are examples of genotype-by-environment (GxE) interactions, gene-by-gene (GxG) interactions (epistasis) or more generally context-dependent effects on the link between genotype and phenotype. We provide evidence from our own studies in Drosophila, and from those of other systems, that mitonuclear interactions—either conflicting or cooperative—are common features of GxE and GxG. We argue that mitonuclear interactions are an important model for how to better understand the pervasive context-dependent effects underlying the architecture of complex phenotypes. Future research in this area should focus on the quantitative genetic concept of effect size to place mitochondrial links to phenotype in a proper context.

This article is part of the theme issue ‘Linking the mitochondrial genotype to phenotype: a complex endeavour’.

Investigating mitonuclear interactions in human admixed populations

To function properly, mitochondria utilize products of 37 mitochondrial and >1,000 nuclear genes, which should be compatible with each other. Discordance between mitochondrial and nuclear genetic ancestry could contribute to phenotypic variation in admixed populations. Here, we explored potential mitonuclear incompatibility in six admixed human populations from the Americas: African Americans, African Caribbeans, Colombians, Mexicans, Peruvians and Puerto Ricans. By comparing nuclear versus mitochondrial ancestry in these populations, we first show that mitochondrial DNA (mtDNA) copy number decreases with increasing discordance between nuclear and mtDNA ancestry. The direction of this effect is consistent across mtDNA haplogroups of different geographic origins. This observation indicates suboptimal regulation of mtDNA replication when its components are encoded by nuclear and mtDNA genes with different ancestry. Second, while most populations analysed exhibit no such trend, in African Americans and Puerto Ricans, we find a significant enrichment of ancestry at nuclear-encoded mitochondrial genes towards the source populations contributing the most prevalent mtDNA haplogroups (African and Native American, respectively). This possibly reflects compensatory effects of selection in recovering mitonuclear interactions optimized in the source populations. Our results provide evidence of mitonuclear interactions in human admixed populations and we discuss their implications for human health and disease.

Admixture Effects on Coevolved Metabolic Systems

Oxidative phosphorylation (OXPHOS) is the primary energy generating system in eukaryotic organisms. The complexes within the OXPHOS pathway are of mixed genomic origin. Although most subunit-coding genes are located within the nuclear genome, several genes are coded for in the mitochondrial genome. There is strong evidence to support coadaptation between the two genomes in these OXPHOS gene regions in order to create tight protein interactions necessary for a functional energetics system. In this study, we begin to assess the physiological impact of separating coevolved protein motifs that make up the highly conserved energy production pathway, as we hypothesize that divergent matings will significantly diminish the protein interactions and therefore hinder efficient OXPHOS activity We measured mitochondrial activity in high energy-demanding tissues from six strains of Mus musculus with varying degrees of mixed ancestral background. Mice with divergent mitochondrial and nuclear backgrounds consistently yielded lower mitochondrial activity. Bioinformatic analysis of common single nucleotide variants across the nuclear and mitochondrial genomes failed to identify any non-synonymous variants that could account for the energetic differences, suggesting that interpopulational mating between ancestrally distinct groups influences energy production efficiency.

The Genetic Architecture of Intra-Species Hybrid Mito-Nuclear Epistasis

Genetic variants that are neutral within, but deleterious between, populations (Dobzhansky-Muller Incompatibilities) are thought to initiate hybrid dysfunction and then to accumulate and complete the speciation process. To identify the types of genetic differences that might initiate speciation, it is useful to study inter-population (intra-species) hybrids that exhibit reduced fitness. In Caenorhabditis briggsae, a close relative of the nematode C. elegans, such minor genetic incompatibilities have been identified. One incompatibility between the mitochondrial and nuclear genomes reduces the fitness of some hybrids. To understand the nuclear genetic architecture of this epistatic interaction, we constructed two sets of recombinant inbred lines by hybridizing two genetically diverse wild populations. In such lines, selection is able to eliminate deleterious combinations of alleles derived from the two parental populations. The genotypes of surviving hybrid lines thus reveal favorable allele combinations at loci experiencing selection. Our genotype data from the resulting lines are consistent with the interpretation that the X alleles participate in epistatic interactions with autosomes and the mitochondrial genome. We evaluate this possibility given predictions that mitochondria-X epistasis should be more prevalent than between mitochondria and autosomes. Our empirical identification of inter-genomic linkage disequilibrium supports the body of literature indicating that the accumulation of mito-nuclear genetic incompatibilities might initiate the speciation process through the generation of less-fit inter-population hybrids.

Comparative genomics of the miniature wasp and pest control agent Trichogramma pretiosum

Background

Trichogrammatids are minute parasitoid wasps that develop within other insect eggs. They are less than half a millimeter long, smaller than some protozoans. The Trichogrammatidae are one of the earliest branching families of Chalcidoidea: a diverse superfamily of approximately half a million species of parasitoid wasps, proposed to have evolved from a miniaturized ancestor. Trichogramma are frequently used in agriculture, released as biological control agents against major moth and butterfly pests. Additionally, Trichogramma are well known for their symbiotic bacteria that induce asexual reproduction in infected females. Knowledge of the genome sequence of Trichogramma is a major step towards further understanding its biology and potential applications in pest control.

Results

We report the 195-Mb genome sequence of Trichogramma pretiosum and uncover signatures of miniaturization and adaptation in Trichogramma and related parasitoids. Comparative analyses reveal relatively rapid evolution of proteins involved in ribosome biogenesis and function, transcriptional regulation, and ploidy regulation. Chalcids also show loss or especially rapid evolution of 285 gene clusters conserved in other Hymenoptera, including many that are involved in signal transduction and embryonic development. Comparisons between sexual and asexual lineages of Trichogramma pretiosum reveal that there is no strong evidence for genome degradation (e.g., gene loss) in the asexual lineage, although it does contain a lower repeat content than the sexual lineage. Trichogramma shows particularly rapid genome evolution compared to other hymenopterans. We speculate these changes reflect adaptations to miniaturization, and to life as a specialized egg parasitoid.

Conclusions

The genomes of Trichogramma and related parasitoids are a valuable resource for future studies of these diverse and economically important insects, including explorations of parasitoid biology, symbiosis, asexuality, biological control, and the evolution of miniaturization. Understanding the molecular determinants of parasitism can also inform mass rearing of Trichogramma and other parasitoids for biological control.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0520-9) contains supplementary material, which is available to authorized users.

The molecular evolutionary dynamics of oxidative phosphorylation (OXPHOS) genes in Hymenoptera

BACKGROUND

The primary energy-producing pathway in eukaryotic cells, the oxidative phosphorylation (OXPHOS) system, comprises proteins encoded by both mitochondrial and nuclear genes. To maintain the function of the OXPHOS system, the pattern of substitutions in mitochondrial and nuclear genes may not be completely independent. It has been suggested that slightly deleterious substitutions in mitochondrial genes are compensated by substitutions in the interacting nuclear genes due to positive selection. Among the four largest insect orders, Coleoptera (beetles), Hymenoptera (sawflies, wasps, ants, and bees), Diptera (midges, mosquitoes, and flies) and Lepidoptera (moths and butterflies), the mitochondrial genes of Hymenoptera exhibit an exceptionally high amino acid substitution rate while the evolution of nuclear OXPHOS genes is largely unknown. Therefore, Hymenoptera is an excellent model group for testing the hypothesis of positive selection driving the substitution rate of nuclear OXPHOS genes. In this study, we report the evolutionary rate of OXPHOS genes in Hymenoptera and test for evidence of positive selection in nuclear OXPHOS genes of Hymenoptera.

RESULTS

Our analyses revealed that the amino acid substitution rate of mitochondrial and nuclear OXPHOS genes in Hymenoptera is higher than that in other studied insect orders. In contrast, the amino acid substitution rate of non-OXPHOS genes in Hymenoptera is lower than the rate in other insect orders. Overall, we found the dN/dS ratio of the nuclear OXPHOS genes to be higher in Hymenoptera than in other insect orders. However, nuclear OXPHOS genes with high dN/dS ratio did not always exhibit a high amino acid substitution rate. Using branch-site and site model tests, we identified various codon sites that evolved under positive selection in nuclear OXPHOS genes.

CONCLUSIONS

Our results showed that nuclear OXPHOS genes in Hymenoptera are evolving faster than the genes in other three insect orders. The branch test suggested that while some nuclear OXPHOS genes in Hymenoptera show a signature of positive selection, the pattern is not consistent across all nuclear OXPHOS genes. As only few codon sites were under positive selection, we suggested that positive selection might not be the only factor contributing to the rapid evolution of nuclear OXPHOS genes in Hymenoptera.

Hybridization increases mitochondrial production of reactive oxygen species in sunfish

Mitochondrial dysfunction and oxidative stress have been suggested to be possible mechanisms underlying hybrid breakdown, as a result of mito-nuclear incompatibilities in respiratory complexes of the electron transport system. However, it remains unclear whether hybridization increases the production of reactive oxygen species (ROS) by mitochondria. We used high-resolution respirometry and fluorometry on isolated liver mitochondria to examine mitochondrial physiology and ROS emission in naturally occurring hybrids of pumpkinseed (Lepomis gibbosus) and bluegill (L. macrochirus). ROS emission was greater in hybrids than in both parent species when respiration was supported by complex I (but not complex II) substrates, and was associated with increases in lipid peroxidation. However, respiratory capacities for oxidative phosphorylation, phosphorylation efficiency, and O₂  kinetics in hybrids were intermediate between those in parental species. Flux control ratios of capacities for electron transport (measured in uncoupled mitochondria) relative to oxidative phosphorylation suggested that the limiting influence of the phosphorylation system is reduced in hybrids. This likely helped offset impairments in electron transport capacity and complex III activity, but contributed to augmenting ROS production. Therefore, hybridization can increase mitochondrial ROS production, in support of previous suggestions that mitochondrial dysfunction can induce oxidative stress and thus contribute to hybrid breakdown.

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.

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

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

Cytoplasmic-Nuclear Incompatibility Between Wild Isolates of Caenorhabditis nouraguensis

How species arise is a fundamental question in biology. Species can be defined as populations of interbreeding individuals that are reproductively isolated from other such populations. Therefore, understanding how reproductive barriers evolve between populations is essential for understanding the process of speciation. Hybrid incompatibility (for example, hybrid sterility or lethality) is a common and strong reproductive barrier in nature. Here we report a lethal incompatibility between two wild isolates of the nematode Caenorhabditis nouraguensis Hybrid inviability results from the incompatibility between a maternally inherited cytoplasmic factor from each strain and a recessive nuclear locus from the other. We have excluded the possibility that maternally inherited endosymbiotic bacteria cause the incompatibility by treating both strains with tetracycline and show that hybrid death is unaffected. Furthermore, cytoplasmic-nuclear incompatibility commonly occurs between other wild isolates, indicating that this is a significant reproductive barrier within C. nouraguensis We hypothesize that the maternally inherited cytoplasmic factor is the mitochondrial genome and that mitochondrial dysfunction underlies hybrid death. This system has the potential to shed light on the dynamics of divergent mitochondrial-nuclear coevolution and its role in promoting speciation.

Asymmetric energetic costs in reciprocal-cross hybrids between carnivorous mice (Onychomys)

Aerobic respiration is a fundamental physiological trait dependent on coordinated interactions between gene products of the mitochondrial and nuclear genomes. Mitonuclear mismatch in interspecific hybrids may contribute to reproductive isolation by inducing reduced viability (or even complete inviability) due to increased metabolic costs. However, few studies have tested for effects of mitonuclear mismatch on respiration at the whole-organism level. We explored how hybridization affects metabolic rate in closely related species of grasshopper mice (genus Onychomys) to better understand the role of metabolic costs in reproductive isolation. We measured metabolic rate across a range of temperatures to calculate basal metabolic rate (BMR) and cold-induced metabolic rate (MR_c) in O. leucogaster, O. torridus and O. arenicola, and in reciprocal F₁ hybrids between the latter two species. Within the genus, we found a negative correlation between mass-specific BMR and body mass. Although O. arenicola was smaller than O. torridus, hybrids from both directions of the cross resembled O. arenicola in body mass. In contrast, hybrid BMR was strongly influenced by the direction of the cross: reciprocal F₁ hybrids were different from each other but indistinguishable from the maternal species. In addition, MR_c was not significantly different between hybrids and either parental species. These patterns indicate that metabolic costs are not increased in Onychomys F₁ hybrids and, while exposure of incompatibilities in F₂ hybrids cannot be ruled out, suggest that mitonuclear mismatch does not act as a primary barrier to gene flow. Maternal matching of BMR is suggestive of a strong effect of mitochondrial genotype on metabolism in hybrids. Together, our findings provide insight into the metabolic consequences of hybridization, a topic that is understudied in mammals.

Disentangling a Holobiont - Recent Advances and Perspectives in Nasonia Wasps

The parasitoid wasp genus Nasonia (Hymenoptera: Chalcidoidea) is a well-established model organism for insect development, evolutionary genetics, speciation, and symbiosis. The host-microbiota assemblage which constitutes the Nasonia holobiont (a host together with all of its associated microbes) consists of viruses, two heritable bacterial symbionts and a bacterial community dominated in abundance by a few taxa in the gut. In the wild, all four Nasonia species are systematically infected with the obligate intracellular bacterium Wolbachia and can additionally be co-infected with Arsenophonus nasoniae. These two reproductive parasites have different transmission modes and host manipulations (cytoplasmic incompatibility vs. male-killing, respectively). Pioneering studies on Wolbachia in Nasonia demonstrated that closely related Nasonia species harbor multiple and mutually incompatible Wolbachia strains, resulting in strong symbiont-mediated reproductive barriers that evolved early in the speciation process. Moreover, research on host-symbiont interactions and speciation has recently broadened from its historical focus on heritable symbionts to the entire microbial community. In this context, each Nasonia species hosts a distinguishable community of gut bacteria that experiences a temporal succession during host development and members of this bacterial community cause strong hybrid lethality during larval development. In this review, we present the Nasonia species complex as a model system to experimentally investigate questions regarding: (i) the impact of different microbes, including (but not limited to) heritable endosymbionts, on the extended phenotype of the holobiont, (ii) the establishment and regulation of a species-specific microbiota, (iii) the role of the microbiota in speciation, and (iv) the resilience and adaptability of the microbiota in wild populations subjected to different environmental pressures. We discuss the potential for easy microbiota manipulations in Nasonia as a promising experimental approach to address these fundamental aspects.

Hybrid Dysfunction Expressed as Elevated Metabolic Rate in Male Ficedula Flycatchers

Studies of ecological speciation are often biased towards extrinsic sources of selection against hybrids, resulting from intermediate hybrid morphology, but the knowledge of how genetic incompatibilities accumulate over time under natural conditions is limited. Here we focus on a physiological trait, metabolic rate, which is central to life history strategies and thermoregulation but is also likely to be sensitive to mismatched mitonuclear interactions. We measured the resting metabolic rate of male collared, and pied flycatchers as well as of naturally occurring F1 hybrid males, in a recent hybrid zone. We found that hybrid males had a higher rather than intermediate metabolic rate, which is indicative of hybrid physiological dysfunction. Fitness costs associated with elevated metabolic rate are typically environmentally dependent and exaggerated under harsh conditions. By focusing on male hybrid dysfunction in an eco-physiological trait, our results contribute to the general understanding of how combined extrinsic and intrinsic sources of hybrid dysfunction build up under natural conditions.

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

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

Gene Regulatory Evolution During Speciation in a Songbird

Over the last decade, tremendous progress has been made toward a comparative understanding of gene regulatory evolution. However, we know little about how gene regulation evolves in birds, and how divergent genomes interact in their hybrids. Because of the unique features of birds – female heterogamety, a highly conserved karyotype, and the slow evolution of reproductive incompatibilities – an understanding of regulatory evolution in birds is critical to a comprehensive understanding of regulatory evolution and its implications for speciation. Using a novel complement of analyses of replicated RNA-seq libraries, we demonstrate abundant divergence in brain gene expression between zebra finch (Taeniopygia guttata) subspecies. By comparing parental populations and their F1 hybrids, we also show that gene misexpression is relatively rare among brain-expressed transcripts in male birds. If this pattern is consistent across tissues and sexes, it may partially explain the slow buildup of postzygotic reproductive isolation observed in birds relative to other taxa. Although we expected that the action of genetic drift on the island-dwelling zebra finch subspecies would be manifested in a higher rate of trans regulatory divergence, we found that most divergence was in cis regulation, following a pattern commonly observed in other taxa. Thus, our study highlights both unique and shared features of avian regulatory evolution.

Cytonuclear interactions affect adaptive traits of the annual plant Arabidopsis thaliana in the field

Although the contribution of cytonuclear interactions to plant fitness variation is relatively well documented at the interspecific level, the prevalence of cytonuclear interactions at the intraspecific level remains poorly investigated. In this study, we set up a field experiment to explore the range of effects that cytonuclear interactions have on fitness-related traits in Arabidopsis thaliana To do so, we created a unique series of 56 cytolines resulting from cytoplasmic substitutions among eight natural accessions reflecting within-species genetic diversity. An assessment of these cytolines and their parental lines scored for 28 adaptive whole-organism phenotypes showed that a large proportion of phenotypic traits (23 of 28) were affected by cytonuclear interactions. The effects of these interactions varied from slight but frequent across cytolines to strong in some specific parental pairs. Two parental pairs accounted for half of the significant pairwise interactions. In one parental pair, Ct-1/Sha, we observed symmetrical phenotypic responses between the two nuclear backgrounds when combined with specific cytoplasms, suggesting nuclear differentiation at loci involved in cytonuclear epistasis. In contrast, asymmetrical phenotypic responses were observed in another parental pair, Cvi-0/Sha. In the Cvi-0 nuclear background, fecundity and phenology-related traits were strongly affected by the Sha cytoplasm, leading to a modified reproductive strategy without penalizing total seed production. These results indicate that natural variation in cytoplasmic and nuclear genomes interact to shape integrative traits that contribute to adaptation, thereby suggesting that cytonuclear interactions can play a major role in the evolutionary dynamics ofA. thaliana.

Mitochondrial-Nuclear Epistasis Impacts Fitness and Mitochondrial Physiology of Interpopulation Caenorhabditis briggsae Hybrids

In order to identify the earliest genetic changes that precipitate species formation, it is useful to study genetic incompatibilities that cause only mild dysfunction when incompatible alleles are combined in an interpopulation hybrid. Such hybridization within the nematode species Caenorhabditis briggsae has been suggested to result in selection against certain combinations of nuclear and mitochondrial alleles, raising the possibility that mitochondrial-nuclear (mitonuclear) epistasis reduces hybrid fitness. To test this hypothesis, cytoplasmic-nuclear hybrids (cybrids) were created to purposefully disrupt any epistatic interactions. Experimental analysis of the cybrids suggests that mitonuclear discord can result in decreased fecundity, increased lipid content, and increased mitochondrial reactive oxygen species levels. Many of these effects were asymmetric with respect to cross direction, as expected if cytoplasmic-nuclear Dobzhansky-Muller incompatibilities exist. One such effect is consistent with the interpretation that disrupting coevolved mitochondrial and nuclear loci impacts mitochondrial function and organismal fitness. These findings enhance efforts to study the genesis, identity, and maintenance of genetic incompatibilities that precipitate the speciation process.

Gene flow between Drosophila yakuba and Drosophila santomea in subunit V of cytochrome c oxidase: A potential case of cytonuclear cointrogression

Introgression is the effective exchange of genetic information between species through natural hybridization. Previous genetic analyses of the Drosophila yakuba—D. santomea hybrid zone showed that the mitochondrial genome of D. yakuba had introgressed into D. santomea and completely replaced its native form. Since mitochondrial proteins work intimately with nuclear‐encoded proteins in the oxidative phosphorylation (OXPHOS) pathway, we hypothesized that some nuclear genes in OXPHOS cointrogressed along with the mitochondrial genome. We analyzed nucleotide variation in the 12 nuclear genes that form cytochrome c oxidase (COX) in 33 Drosophila lines. COX is an OXPHOS enzyme composed of both nuclear‐ and mitochondrial‐encoded proteins and shows evidence of cytonuclear coadaptation in some species. Using maximum‐likelihood methods, we detected significant gene flow from D. yakuba to D. santomea for the entire COX complex. Interestingly, the signal of introgression is concentrated in the three nuclear genes composing subunit V, which shows population migration rates significantly greater than the background level of introgression in these species. The detection of introgression in three proteins that work together, interact directly with the mitochondrial‐encoded core, and are critical for early COX assembly suggests this could be a case of cytonuclear cointrogression.

Gametic selection, developmental trajectories, and extrinsic heterogeneity in Haldane's rule

Deciphering the genetic and developmental causes of the disproportionate rarity, inviability, and sterility of hybrid males, Haldane's rule, is important for understanding the evolution of reproductive isolation between species. Moreover, extrinsic and prezygotic factors can contribute to the magnitude of intrinsic isolation experienced between species with partial reproductive compatibility. Here, we use the nematodes Caenorhabditis briggsae and C. nigoni to quantify the sensitivity of hybrid male viability to extrinsic temperature and developmental timing, and test for a role of mito-nuclear incompatibility as a genetic cause. We demonstrate that hybrid male inviability manifests almost entirely as embryonic, not larval, arrest and is maximal at the lowest rearing temperatures, indicating an intrinsic-by-extrinsic interaction to hybrid inviability. Crosses using mitochondrial substitution strains that have reciprocally introgressed mitochondrial and nuclear genomes show that mito-nuclear incompatibility is not a dominant contributor to postzygotic isolation and does not drive Haldane's rule in this system. Crosses also reveal that competitive superiority of X-bearing sperm provides a novel means by which postmating prezygotic factors exacerbate the rarity of hybrid males. These findings highlight the important roles of gametic, developmental, and extrinsic factors in modulating the manifestation of Haldane's rule.

Mitonuclear Ecology

Eukaryotes were born of a chimeric union between two prokaryotes--the progenitors of the mitochondrial and nuclear genomes. Early in eukaryote evolution, most mitochondrial genes were lost or transferred to the nucleus, but a core set of genes that code exclusively for products associated with the electron transport system remained in the mitochondrion. The products of these mitochondrial genes work in intimate association with the products of nuclear genes to enable oxidative phosphorylation and core energy production. The need for coadaptation, the challenge of cotransmission, and the possibility of genomic conflict between mitochondrial and nuclear genes have profound consequences for the ecology and evolution of eukaryotic life. An emerging interdisciplinary field that I call "mitonuclear ecology" is reassessing core concepts in evolutionary ecology including sexual reproduction, two sexes, sexual selection, adaptation, and speciation in light of the interactions of mitochondrial and nuclear genomes.

Hybrid incompatibilities are affected by dominance and dosage in the haplodiploid wasp Nasonia

Study of genome incompatibilities in species hybrids is important for understanding the genetic basis of reproductive isolation and speciation. According to Haldane's rule hybridization affects the heterogametic sex more than the homogametic sex. Several theories have been proposed that attribute asymmetry in hybridization effects to either phenotype (sex) or genotype (heterogamety). Here we investigate the genetic basis of hybrid genome incompatibility in the haplodiploid wasp Nasonia using the powerful features of haploid males and sex reversal. We separately investigate the effects of heterozygosity (ploidy level) and sex by generating sex reversed diploid hybrid males and comparing them to genotypically similar haploid hybrid males and diploid hybrid females. Hybrid effects of sterility were more pronounced than of inviability, and were particularly strong in haploid males, but weak to absent in diploid males and females, indicating a strong ploidy level but no sex specific effect. Molecular markers identified a number of genomic regions associated with hybrid inviability in haploid males that disappeared under diploidy in both hybrid males and females. Hybrid inviability was rescued by dominance effects at some genomic regions, but aggravated or alleviated by dosage effects at other regions, consistent with cytonuclear incompatibilities. Dosage effects underlying Bateson–Dobzhansky–Muller (BDM) incompatibilities need more consideration in explaining Haldane's rule in diploid systems.