Functional coadaptation between cytochrome c and cytochrome c oxidase within allopatric populations of a marine copepod

Mechanisms Maintaining Mitochondrial DNA Polymorphisms: The Role of Mito-Nuclear Interactions, Sex-Specific Selection, and Genotype-by-Environment Interactions in Drosophila subobscura

Experimental mito-nuclear introgression lines (MNILs) were established by backcrossing isofemale lines of D. subobscura originating from the same populations. MNILs were subjected to a series of life-history experiments designed to test the fitness of the bearers of different combinations of two main mtDNA haplotypes on their own nuclear background, as well as on the background of the opposite haplotype. By having 11 replicas of the four mito-nuclear combinations, we could test not only the adaptive significance of the differences between the two main haplotypes but also the influence of additional variation present within each of the 11 combinations on fitness. Testing the fitness of individuals of both sexes enabled us to examine if sex-specific selection has a role in maintaining the frequencies of the two mtDNA haplotypes in nature. Conducting the fitness assays on two different temperatures enabled us to test whether different temperatures favor specific mtDNA haplotypes or mito-nuclear genotypes and consequently promote stable sympatric mtDNA variation. The results show weak signature of genotype-by-environment interactions, and no sex-specific selection regarding differences between the two main haplotypes. However, individual models across different life-history components showed these two mechanisms at play in promoting mtDNA variability present in specific mito-nuclear crosses. Our models show that mito-nuclear interactions are, in fact, more important as units of selection.

Is There Hybrid Vigor in Dogs? Oxidative Stress and Cytokine Concentrations in Low- To Mid-Content Wolf-Dog Hybrids

AbstractChanges of deleterious traits in mitochondria within hybrids of two different species are likely to be reflected across whole-animal phenotypes. Two processes linked to mitochondrial fitness are oxidative stress and inflammation. Here, plasma oxidative stress (lipid oxidative damage, total antioxidant capacity, and activities of catalase, glutathione peroxidase, and superoxide dismutase concentration) and cytokine concentrations (IL-1β, IL-6, and TNF-α) from wolf-dog hybrids of low to mid wolf content are presented and then compared with previously published values for similar-sized dogs and gray wolves. Results indicate that lipid oxidative damage and catalase activity were not significantly different across species and hybrids. Total antioxidant capacity and glutathione peroxidase were significantly lower in wolf-dog hybrids than in domestic dogs and wolves (although total antioxidant capacity significantly increased with wolf content), but superoxide dismutase concentration was higher in hybrids. Thus, it seems that a low percentage of wolf content decreases several aspects of antioxidants but without any accumulating lipid oxidative damage. Additionally, wolf-dog hybrids had higher IL-1β and IL-6 concentrations but lower TNF-α concentrations than domestic dogs. And there was a significantly positive correlation between percentage of wolf content and IL-1β. These data imply either a mitonuclear incompatibility or a nuclear-nuclear incompatibility within wolf-dog hybrids, a case that does not support heterosis.

Mitonuclear effects on sex ratio persist across generations in interpopulation hybrids

Eukaryotic energy production requires tight coordination between nuclear and mitochondrial gene products. Because males and females often have different energetic strategies, optimal mitonuclear coordination may be sex-specific. Previous work found evidence for sex-specific mitonuclear effects in the copepod Tigriopus californicus by comparing two parental lines and their reciprocal F1 crosses. However, an alternative hypothesis is that the patterns were driven by the parental source of nuclear alleles. Here, we test this alternative hypothesis by extending the same cross to F2 hybrids, which receive both maternal and paternal nuclear alleles from F1 hybrids. Results confirm mitonuclear effects on sex ratio, with distorted ratios persisting from the F1 to F2 generations, despite reduced fitness in F2 hybrids. No sex-by-cross interactions were found for other phenotypic traits measured. Mitochondrial DNA content was higher in females. Both routine metabolic rate and oxidative DNA damage were lower in F2 hybrids than in parentals. The persistence of sex-specific mitonuclear effects, even in the face of F2 hybrid breakdown, attests to the magnitude of these effects, which contribute to the maintenance of within-population mitochondrial DNA polymorphisms.

Evolutionary codependency: insights into the mitonuclear interaction landscape from experimental and wild Caenorhabditis nematodes

Aided by new technologies, the upsurgence of research into mitochondrial genome biology during the past 15 years suggests that we have misunderstood, and perhaps dramatically underestimated, the ongoing biological and evolutionary significance of our long-time symbiotic partner. While we have begun to scratch the surface of several topics, many questions regarding the nature of mutation and selection in the mitochondrial genome, and the nature of its relationship to the nuclear genome, remain unanswered. Although best known for their contributions to studies of developmental and aging biology, Caenorhabditis nematodes are increasingly recognized as excellent model systems to advance understanding in these areas. We review recent discoveries with relevance to mitonuclear coevolution and conflict and offer several fertile areas for future work.

Metformin disrupts Danio rerio metabolism at environmentally relevant concentrations: A full life-cycle study

Metformin (MET), an anti-diabetic pharmaceutical of large-scale consumption, is increasingly detected in surface waters. However, current knowledge on the long-term effects of MET on non-target organisms is limited. The present study aimed to investigate the effects of MET in the model freshwater teleost Danio rerio, following a full life-cycle exposure to environmentally relevant concentrations (390 to 14 423 ng/L). Considering that the mode of action (MoA) of MET on non-target organisms remains underexplored and that MET may act through similar human pathways, i.e., lipid and energy metabolisms, biochemical markers were used to determine cholesterol and triglycerides levels, as well as mitochondrial complex I activity in zebrafish liver. Also, the hepatosomatic index as an indication of metabolic disruption, and the expression levels of genes involved in MET's putative MoA, i.e. acaca, acadm, cox5aa, idh3a, hmgcra, prkaa1, were determined, the last by qRT-PCR. A screening of mRNA transcripts, associated with lipid and energy metabolisms, and other signaling pathways potentially involved in MET-induced toxicity were also assessed using an exploratory RNA-seq analysis. The findings here reported indicate that MET significantly disrupted critical biochemical and molecular processes involved in zebrafish metabolism, such as cholesterol and fatty acid biosynthesis, mitochondrial electron transport chain and tricarboxylic acid cycle, concomitantly to changes on the hepatosomatic index. Likewise, MET impacted other relevant pathways mainly associated with cell cycle, DNA repair and steroid hormone biosynthesis, here reported for the first time in a non-target aquatic organism. Non-monotonic dose response curves were frequently detected in biochemical and qRT-PCR data, with higher effects observed at 390 and 2 929 ng/L MET treatments. Collectively, the results suggest that environmentally relevant concentrations of MET severely disrupt D. rerio metabolism and other important biological processes, supporting the need to revise the proposed environmental quality standard (EQS) and predicted no-effect concentration (PNEC) for MET.

Ecological speciation by sympatric host shifts in a clade of herbivorous sea slugs, with introgression and localized mitochondrial capture between species

Host shifting in insect-plant systems was historically important to the development of ecological speciation theory, yet surprisingly few studies have examined whether host shifting drives diversification of marine herbivores. When small-bodied consumers feed and also mate on a preferred host, disruptive selection can split a population into host races despite gene flow. Support for host shifts is notably lacking for invertebrates associated with macroalgae, where the scale of dispersal by planktonic larvae often far exceeds the grain of host patchiness, and adults are typically less specialized than terrestrial herbivores. Here, we present a candidate example of ecological speciation in a clade of sea slugs that primarily consume green algae in the genus Caulerpa, including highly invasive species. Ancestral character state reconstructions supported 'sea grapes' (C. racemosa, C. lentillifera) as the ancestral host for a tropical radiation of 12 Elysia spp., with one shift onto alternative Caulerpa spp. in the Indo-Pacific. A Caribbean radiation of three species included symaptric host shifts to Rhipocephalus brevicaulis in the ancestor of E. pratensis Ortea & Espinosa, 1996, and to C. prolifera in E. hamanni Krug, Vendetti & Valdes 2016, plus a niche expansion to a range of Caulerpa spp. in E. subornata Verrill, 1901. All three species are broadly sympatric across the Caribbean but are host-partitioned at a fine grain, and distinct by morphology and at nuclear loci. However, non-recombining mtDNA revealed a history of gene flow between E. pratensis and E. subornata: COI haplotypes from E. subornata were 10.4% divergent from E. pratensis haplotypes from four sites, but closely related to all E. pratensis haplotypes sampled from six Bahamian islands, indicating historical introgression and localized "mitochondrial capture." Disruptive selective likely fueled divergence and adaptation to distinct host environments, indicating ecological speciation may be an under-appreciated driver of diversification for marine herbivores as well as epibionts and other resource specialists.

Frequent Paternal Mitochondrial Inheritance and Rapid Haplotype Frequency Shifts in Copepod Hybrids

Mitochondria are assumed to be maternally inherited in most animal species, and this foundational concept has fostered advances in phylogenetics, conservation, and population genetics. Like other animals, mitochondria were thought to be solely maternally inherited in the marine copepod Tigriopus californicus, which has served as a useful model for studying mitonuclear interactions, hybrid breakdown, and environmental tolerance. However, we present PCR, Sanger sequencing, and Illumina Nextera sequencing evidence that extensive paternal mitochondrial DNA (mtDNA) transmission is occurring in inter-population hybrids of T. californicus. PCR on four types of crosses between three populations (total sample size of 376 F1 individuals) with 20% genome-wide mitochondrial divergence showed 2% to 59% of F1 hybrids with both paternal and maternal mtDNA, where low and high paternal leakage values were found in different cross directions of the same population pairs. Sequencing methods further verified nucleotide similarities between F1 mtDNA and paternal mtDNA sequences. Interestingly, the paternal mtDNA in F1s from some crosses inherited haplotypes that were uncommon in the paternal population. Compared to some previous research on paternal leakage, we employed more rigorous methods to rule out contamination and false detection of paternal mtDNA due to non-functional nuclear mitochondrial DNA fragments. Our results raise the potential that other animal systems thought to only inherit maternal mitochondria may also have paternal leakage, which would then affect the interpretation of past and future population genetics or phylogenetic studies that rely on mitochondria as uniparental markers.

Phenotypic disruption of cuticular hydrocarbon production in hybrids between sympatric species of Hawaiian picture-wing Drosophila

Interspecies hybrids can express phenotypic traits far outside the range of parental species. The atypical traits of hybrids provide insight into differences in the factors that regulate the expression of these traits in the parental species. In some cases, the unusual phenotypic traits of hybrids can lead to phenotypic dysfunction with hybrids experiencing reduced survival or reproduction. Cuticular hydrocarbons (CHCs) in insects are important phenotypic traits that serve several functions, including desiccation resistance and pheromones for mating. We used gas chromatography mass spectrometry to investigate the differences in CHC production between two closely related sympatric Hawaiian picture-wing Drosophila species, Drosophila heteroneura and D. silvestris, and their F1 and backcross hybrid offspring. CHC profiles differed between males of the two species, with substantial sexual dimorphism in D. silvestris but limited sexual dimorphism in D. heteroneura. Surprisingly, F1 hybrids did not produce three CHCs, and the abundances of several other CHCs occurred outside the ranges present in the two parental species. Backcross hybrids produced all CHCs with greater variation than observed in F1 or parental species. Overall, these results suggest that the production of CHCs was disrupted in F1 and backcross hybrids, which may have important consequences for their survival or reproduction.

Mitonuclear Mismatch is Associated With Increased Male Frequency, Outcrossing, and Male Sperm Size in Experimentally-Evolved C. elegans

We provide a partial test of the mitonuclear sex hypothesis with the first controlled study of how male frequencies and rates of outcrossing evolve in response to mitonuclear mismatch by allowing replicate lineages of C. elegans nematodes containing either mitochondrial or nuclear mutations of electron transport chain (ETC) genes to evolve under three sexual systems: facultatively outcrossing (wildtype), obligately selfing, and obligately outcrossing. Among facultatively outcrossing lines, we found evolution of increased male frequency in at least one replicate line of all four ETC mutant backgrounds tested—nuclear isp-1, mitochondrial cox-1 and ctb-1, and an isp-1 IV; ctb-1M mitonuclear double mutant—and confirmed for a single line set (cox-1) that increased male frequency also resulted in successful outcrossing. We previously found the same result for lines evolved from another nuclear ETC mutant, gas-1. For several lines in the current experiment, however, male frequency declined to wildtype levels (near 0%) in later generations. Male frequency did not change in lines evolved from a wildtype control strain. Additional phenotypic assays of lines evolved from the mitochondrial cox-1 mutant indicated that evolution of high male frequency was accompanied by evolution of increased male sperm size and mating success with tester females, but that it did not translate into increased mating success with coevolved hermaphrodites. Rather, hermaphrodites’ self-crossed reproductive fitness increased, consistent with sexually antagonistic coevolution. In accordance with evolutionary theory, males and sexual outcrossing may be most beneficial to populations evolving from a state of low ancestral fitness (gas-1, as previously reported) and less beneficial or deleterious to those evolving from a state of higher ancestral fitness (cox-1). In support of this idea, the obligately outcrossing fog-2 V; cox-1 M lines exhibited no fitness evolution compared to their ancestor, while facultatively outcrossing lines showed slight upward evolution of fitness, and all but one of the obligately selfing xol-1 X; cox-1 M lines evolved substantially increased fitness—even beyond wildtype levels. This work provides a foundation to directly test the effect of reproductive mode on the evolutionary dynamics of mitonuclear genomes, as well as whether compensatory mutations (nuclear or mitochondrial) can rescue populations from mitochondrial dysfunction.

The transcriptomic signature of obligate parthenogenesis

Investigating the origin of parthenogenesis through interspecific hybridization can provide insight into how meiosis may be altered by genetic incompatibilities, which is fundamental for our understanding of the formation of reproductive barriers. Yet the genetic mechanisms giving rise to obligate parthenogenesis in eukaryotes remain understudied. In the microcrustacean Daphnia pulex species complex, obligately parthenogenetic (OP) isolates emerged as backcrosses of two cyclically parthenogenetic (CP) parental species, D. pulex and D. pulicaria, two closely related but ecologically distinct species. We examine the genome-wide expression in OP females at the early resting egg production stage, a life-history stage distinguishing OP and CP reproductive strategies, in comparison to CP females of the same stage from the two parental species. Our analyses of the expression data reveal that underdominant and overdominant genes are abundant in OP isolates, suggesting widespread regulatory incompatibilities between the parental species. More importantly, underdominant genes (i.e., genes with expression lower than both parentals) in the OP isolates are enriched in meiosis and cell-cycle pathways, indicating an important role of underdominance in the origin of obligate parthenogenesis. Furthermore, metabolic and biosynthesis pathways enriched with overdominant genes (i.e., expression higher than both parentals) are another genomic signature of OP isolates.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

Evaluating Assisted Gene Flow in Marginal Populations of a High Mountain Species

Many species cannot either migrate or adapt at the rate of temperature increases due to climate warming. Therefore, they need active conservation strategies to avoid extinction. Facilitated adaptation actions, such assisted gene flow, aim at the increase of the evolutionary resilience of species affected by global change. In elevational gradients, marginal populations at the lower elevation edges are experiencing earlier snowmelt and higher temperatures, which force them to adapt to the new conditions by modifying their phenology. In this context, advancing the onset of flowering and seed germination times are crucial to ensure reproductive success and increase seedling survival prior to summer drought. Assisted gene flow may bring adaptive alleles and increase genetic diversity that can help throughout ontogeny. The main aim of this work is to assess the effects that different gene flow treatments could have on the desired trait changes in marginal populations. Accordingly, we established a common garden experiment in which we assayed four different gene flow treatments between Silene ciliata Pourr. (Caryophyllaceae) populations located in similar and different elevation edges, belonging to the same and different mountains. As a control treatment, within-population crosses of low elevation edge populations were performed. The resulting seeds were sown and the germination and flowering onset dates of the resulting plants recorded, as well as the seedling survival. Gene flow between populations falling on the same mountain and same elevation and gene flow from high-elevation populations from a different mountain to low-elevation populations advanced seed germination time with respect to control crosses. No significant effects of gene flow on seedling survival were found. All the gene flow treatments delayed the onset of flowering with respect to control crosses and this effect was more pronounced in among-mountain gene flows. The results of this study highlight two important issues that should be thoroughly studied before attempting to apply assisted gene flow in practical conservation situations. Firstly, among-populations gene flow can trigger different responses in crucial traits throughout the ontogeny of plant species. Secondly, the population provenance of gene flow is determinant and plays a significant role on the effects of gene flow.

The life-history fitness of F1 hybrids of the microcrustacean Daphnia pulex and D. pulicaria (Crustacea, Anomopoda)

Negative interaction between alleles that arise independently in diverging populations (i.e., Dobzhansky-Muller incompatibilities) can cause reduction of fitness in their hybrids. However, heterosis in hybrids can emerge if hybridization breaks down detrimental epistatic interaction within parental lineages. In this study, we examined the life-history fitness of the inter-specific F₁s of two recently diverged microcrustacean species Daphnia pulex and D. pulicaria as well as intra-specific crosses of D. pulex. We identified heterosis in two out of five life-history traits in the inter-specific F₁s. According to theories that heterosis can transiently emerge in early speciation, the observation of heterosis in these life-history traits suggests that there are no major genetic incompatibilities between these two species affecting these traits and that D. pulex and D. pulicaria are at an early stage of speciation.

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.

Embryo, larval, and juvenile staging of Lytechinus pictus from fertilization through sexual maturation

BACKGROUND

Sea urchin embryos have been used for more than a century in the study of fertilization and early development. However, several of the species used, such as Strongylocentrotus purpuratus, have long generation times making them suboptimal for transgenerational studies.

RESULTS

Here, we present an overview of the development of a rapidly developing echinoderm species, Lytechinus pictus, from fertilization through sexual maturation. When grown at room temperature (20°C) embryos complete the first cell cycle in 90 minutes, followed by subsequent cleavages every 45 minutes, leading to hatching at 9 hours postfertilization (hpf). The swimming embryos gastrulate from 12 to 36 hpf and produce the cells which subsequently give rise to the larval skeleton and immunocytes. Larvae begin to feed at 2 days and metamorphose by 3 weeks. Juveniles reach sexual maturity at 4 to 6 months of age, depending on individual growth rate.

CONCLUSIONS

This staging scheme lays a foundation for future studies in L. pictus, which share many of the attractive features of other urchins but have the key advantage of rapid development to sexual maturation. This is significant for multigenerational and genetic studies newly enabled by CRISPR-CAS mediated gene editing.

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.

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’.

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

The mitochondrion is a pivotal organelle for energy production, and includes components encoded by both the mitochondrial and nuclear genomes. Functional and evolutionary interactions are expected between the nuclear- and mitochondrial-encoded components. The topic is of broad interest in biology, with implications to genetics, evolution, and medicine. Here, we compare the evolutionary rates of mitochondrial proteins and ribosomal RNAs to rates of mitochondria-associated nuclear-encoded proteins, across the major orders of holometabolous insects. There are significant evolutionary rate correlations (ERCs) between mitochondrial-encoded and mitochondria-associated nuclear-encoded proteins, which are likely driven by different rates of mitochondrial sequence evolution and correlated changes in the interacting nuclear-encoded proteins. The pattern holds after correction for phylogenetic relationships and considering protein conservation levels. Correlations are stronger for both nuclear-encoded OXPHOS proteins that are in contact with mitochondrial OXPHOS proteins and for nuclear-encoded mitochondrial ribosomal amino acids directly contacting the mitochondrial rRNAs. We find that ERC between mitochondrial- and nuclear-encoded proteins is a strong predictor of nuclear-encoded proteins known to interact with mitochondria, and ERC shows promise for identifying new candidate proteins with mitochondrial function. Twenty-three additional candidate nuclear-encoded proteins warrant further study for mitochondrial function based on this approach, including proteins in the minichromosome maintenance helicase complex.

Assessing the fitness consequences of mitonuclear interactions in natural populations

Metazoans exist only with a continuous and rich supply of chemical energy from oxidative phosphorylation in mitochondria. The oxidative phosphorylation machinery that mediates energy conservation is encoded by both mitochondrial and nuclear genes, and hence the products of these two genomes must interact closely to achieve coordinated function of core respiratory processes. It follows that selection for efficient respiration will lead to selection for compatible combinations of mitochondrial and nuclear genotypes, and this should facilitate coadaptation between mitochondrial and nuclear genomes (mitonuclear coadaptation). Herein, we outline the modes by which mitochondrial and nuclear genomes may coevolve within natural populations, and we discuss the implications of mitonuclear coadaptation for diverse fields of study in the biological sciences. We identify five themes in the study of mitonuclear interactions that provide a roadmap for both ecological and biomedical studies seeking to measure the contribution of intergenomic coadaptation to the evolution of natural populations. We also explore the wider implications of the fitness consequences of mitonuclear interactions, focusing on central debates within the fields of ecology and biomedicine.

Rapid divergence, molecular evolution, and morphological diversification of coastal host-parasite systems from southern Brazil

This study assessed the role of historical processes on the geographic isolation, molecular evolution, and morphological diversification of host-parasite populations from the southern Brazilian coast. Adult specimens of Scleromystax barbatus and Scleromystax macropterus were collected from the sub-basin of the Nhundiaquara River and the sub-basin of the Paranaguá Bay, state of Paraná, Brazil. Four species of Gyrodactylus were recovered from the body surface of both host species. Morphometric analysis of Gyrodactylus spp. and Scleromystax spp. indicated that subpopulations of parasites and hosts could be distinguished from different sub-basins and locations, but the degree of morphological differentiation seems to be little related to geographic distance between subpopulations. Phylogenetic relationships based on DNA sequences of Gyrodactylus spp. and Scleromystax spp. allowed distinguishing lineages of parasites and hosts from different sub-basins. However, the level of genetic structuring of parasites was higher in comparison to host species. Evidence of positive selection in mtDNA sequences is likely associated with local adaptation of lineages of parasites and hosts. A historical demographic analysis revealed that populations of Gyrodactylus and Scleromystax have expanded in the last 250 000 years. The genetic variation of parasites and hosts is consistent with population-specific selection, population expansions, and recent evolutionary co-divergence.

The evolution of hybrid fitness during speciation

The evolution of postzygotic reproductive isolation is an important component of speciation. But before isolation is complete there is sometimes a phase of heterosis in which hybrid fitness exceeds that of the two parental species. The genetics and evolution of heterosis and postzygotic isolation have typically been studied in isolation, precluding the development of a unified theory of speciation. Here, we develop a model that incorporates both positive and negative gene interactions, and accounts for the evolution of both heterosis and postzygotic isolation. We parameterize the model with recent data on the fitness effects of 10,000 mutations in yeast, singly and in pairwise epistatic combinations. The model makes novel predictions about the types of interactions that contribute to declining hybrid fitness. We reproduce patterns familiar from earlier models of speciation (e.g. Haldane's Rule and Darwin's Corollary) and identify new mechanisms that may underlie these patterns. Our approach provides a general framework for integrating experimental data from gene interaction networks into speciation theory and makes new predictions about the genetic mechanisms of speciation.

Sex-specific stress tolerance, proteolysis, and lifespan in the invertebrate Tigriopus californicus

Because stress tolerance and longevity are mechanistically and phenotypically linked, the sex with higher acute stress tolerance might be expected to also live longer. On the other hand, the association between stress tolerance and lifespan may be complicated by tradeoffs between acute tolerance and long-term survival. Here we use the copepod Tigriopus californicus to test for sex differences in stress resistance, proteolytic activity and longevity. Unlike many model organisms, this species does not have sex chromosomes. However, substantial sex differences were still observed. Females were found to have superior tolerance to a range of acute stressors (high temperature, high salinity, low salinity, copper and bisphenol A (BPA)) across a variety of treatments including different populations, pure vs. hybrid crosses, and different shading environments. Upregulation of proteolytic capacity - one molecular mechanism for responding to acute stress - was also found to be sexually dimorphic. In the combined stress treatment of chronic copper exposure followed by acute heat exposure, proteolytic capacity was suppressed for males. Females, however, maintained a robust proteolytic stress response. While females consistently showed greater tolerance to short-term stress, lifespan was largely equivalent between the two sexes under both benign conditions and mild thermal stress. Our findings indicate that short-term stress tolerance does not predict long-term survival under relatively mild conditions.

Sex and Mitonuclear Adaptation in Experimental Caenorhabditis elegans Populations

To reveal phenotypic and functional genomic patterns of mitonuclear adaptation, a laboratory adaptation study with Caenorhabditis elegans nematodes was conducted in which independently evolving lines were initiated from a low-fitness mitochondrial electron transport chain (ETC) mutant, gas-1 Following 60 generations of evolution in large population sizes with competition for food resources, two distinct classes of lines representing different degrees of adaptive response emerged: a low-fitness class that exhibited minimal or no improvement compared to the gas-1 mutant ancestor, and a high-fitness class containing lines that exhibited partial recovery of wild-type fitness. Many lines that achieved higher reproductive and competitive fitness levels were also noted to evolve high frequencies of males during the experiment, consistent with adaptation in these lines having been facilitated by outcrossing. Whole-genome sequencing and analysis revealed an enrichment of mutations in loci that occur in a gas-1-centric region of the C. elegans interactome and could be classified into a small number of functional genomic categories. A highly nonrandom pattern of mitochondrial DNA mutation was observed within high-fitness gas-1 lines, with parallel fixations of nonsynonymous base substitutions within genes encoding NADH dehydrogenase subunits I and VI. These mitochondrial gene products reside within ETC complex I alongside the nuclear-encoded GAS-1 protein, suggesting that rapid adaptation of select gas-1 recovery lines was driven by fixation of compensatory mitochondrial mutations.

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.

Genetic Variation for Mitochondrial Function in the New Zealand Freshwater Snail Potamopyrgus antipodarum

The proteins responsible for mitochondrial function are encoded by 2 different genomes with distinct inheritance regimes, rendering rigorous inference of genotype-phenotype connections intractable for all but a few model systems. Asexual organisms provide a powerful means to address these challenges because offspring produced without recombination inherit both nuclear and mitochondrial genomes from a single parent. As such, these offspring inherit mitonuclear genotypes that are identical to the mitonuclear genotypes of their parents and siblings but different from those of other asexual lineages. Here, we compared mitochondrial function across distinct asexual lineages of Potamopyrgus antipodarum, a New Zealand freshwater snail model for understanding the evolutionary consequences of asexuality. Our analyses revealed substantial phenotypic variation across asexual lineages at 3 levels of biological organization: mitogenomic, organellar, and organismal. These data demonstrate that different asexual lineages have different mitochondrial function phenotypes, likely reflecting heritable variation (i.e., the raw material for evolution) for mitochondrial function in P. antipodarum. The discovery of this variation combined with the methods developed here sets the stage to use P. antipodarum to study central evolutionary questions involving mitochondrial function, including whether mitochondrial mutation accumulation influences the maintenance of sexual reproduction in natural populations.

Comparative biochemistry of cytochrome c oxidase in animals

Cytochrome c oxidase (COX), the terminal enzyme of the electron transport system, is central to aerobic metabolism of animals. Many aspects of its structure and function are highly conserved, yet, paradoxically, it is also an important model for studying the evolution of the metabolic phenotype. In this review, part of a special issue honouring Peter Hochachka, we consider the biology of COX from the perspective of comparative and evolutionary biochemistry. The approach is to consider what is known about the enzyme in the context of conventional biochemistry, but focus on how evolutionary researchers have used this background to explore the role of the enzyme in biochemical adaptation of animals. In synthesizing the conventional and evolutionary biochemistry, we hope to identify synergies and future research opportunities. COX represents a rare opportunity for researchers to design studies that span the breadth of biology: molecular genetics, protein biochemistry, enzymology, metabolic physiology, organismal performance, evolutionary biology, and phylogeography.

Cellular and population level processes influence the rate, accumulation and observed frequency of inherited and somatic mtDNA mutations

Mitochondria are found in all animals and have the unique feature of containing multiple copies of their own small, circular DNA genome (mtDNA). The rate and pattern of mutation accumulation in the mtDNA are influenced by molecular, cellular and population level processes. We distinguish between inherited and somatic mtDNA mutations and review evidence for the often-made assumption that mutations accumulate at a higher rate in mtDNA than in nuclear DNA (nDNA). We conclude that the whole genome mutation accumulation rate is higher for mtDNA than for nDNA but include the caveat that rates overlap considerably between the individual mtDNA- and nDNA-encoded genes. Next, we discuss the postulated causal mechanisms for the high rate of mtDNA mutation accumulation in both inheritance and in somatic cells. Perhaps unexpectedly, mtDNA is resilient to many mutagens of nDNA but is prone to errors of replication. We then consider the influence of maternal inheritance, recombination and selection on the observed accumulation pattern of inherited mtDNA mutations. Finally, we discuss environmental influences of temperature and diet on the observed frequency of inherited and somatic mtDNA mutations. We conclude that it is necessary to understand the cellular processes to fully interpret the pattern of mutations and how they influence our interpretations of evolution and disease.

Drosophila mitotypes determine developmental time in a diet and temperature dependent manner

It is well known that specific mitochondrial (mt) DNA mutations can reduce organismal fitness and influence mitochondrial-nuclear interactions. However, determining specific mtDNA mutations that are beneficial has been elusive. In this study, we vary the diet and environmental temperature to study larval development time of two Drosophila melanogaster mitotypes (Alstonville and Dahomey), in two nuclear genetic backgrounds, and investigate developmental differences through weight, feeding rate, and movement. To manipulate the diet, we utilize the nutritional geometric framework to manipulate isocaloric diets of differing macronutrient ratios (1:2 and 1:16 protein: carbohydrate (P:C) ratios) and raise flies at three temperatures (19°C, 23°C and 27°C). Larvae with Dahomey mtDNA develop more slowly than Alstonville when fed the 1:2 P:C diet at all temperatures and developed more quickly when fed the 1:16 P:C diet at 23°C and 27°C. We determined that Dahomey larvae eat more, move less, and weigh more than Alstonville larvae when raised on the 1:16 P:C diet and that these physiological responses are modified by temperature. We suggest that 1 (or more) of 4 mtDNA changes is likely responsible for the observed effects and posit the mtDNA changes moderate a physiological trade-off between consumption and foraging.

Evolved genetic and phenotypic differences due to mitochondrial-nuclear interactions

The oxidative phosphorylation (OxPhos) pathway is responsible for most aerobic ATP production and is the only pathway with both nuclear and mitochondrial encoded proteins. The importance of the interactions between these two genomes has recently received more attention because of their potential evolutionary effects and how they may affect human health and disease. In many different organisms, healthy nuclear and mitochondrial genome hybrids between species or among distant populations within a species affect fitness and OxPhos functions. However, what is less understood is whether these interactions impact individuals within a single natural population. The significance of this impact depends on the strength of selection for mito-nuclear interactions. We examined whether mito-nuclear interactions alter allele frequencies for ~11,000 nuclear SNPs within a single, natural Fundulus heteroclitus population containing two divergent mitochondrial haplotypes (mt-haplotypes). Between the two mt-haplotypes, there are significant nuclear allele frequency differences for 349 SNPs with a p-value of 1% (236 with 10% FDR). Unlike the rest of the genome, these 349 outlier SNPs form two groups associated with each mt-haplotype, with a minority of individuals having mixed ancestry. We use this mixed ancestry in combination with mt-haplotype as a polygenic factor to explain a significant fraction of the individual OxPhos variation. These data suggest that mito-nuclear interactions affect cardiac OxPhos function. The 349 outlier SNPs occur in genes involved in regulating metabolic processes but are not directly associated with the 79 nuclear OxPhos proteins. Therefore, we postulate that the evolution of mito-nuclear interactions affects OxPhos function by acting upstream of OxPhos.

Integrative Approaches for Studying Mitochondrial and Nuclear Genome Co-evolution in Oxidative Phosphorylation

In animals, interactions among gene products of mitochondrial and nuclear genomes (mitonuclear interactions) are of profound fitness, evolutionary, and ecological significance. Most fundamentally, the oxidative phosphorylation (OXPHOS) complexes responsible for cellular bioenergetics are formed by the direct interactions of 13 mitochondrial-encoded and ∼80 nuclear-encoded protein subunits in most animals. It is expected that organisms will develop genomic architecture that facilitates co-adaptation of these mitonuclear interactions and enhances biochemical efficiency of OXPHOS complexes. In this perspective, we present principles and approaches to understanding the co-evolution of these interactions, with a novel focus on how genomic architecture might facilitate it. We advocate that recent interdisciplinary advances assist in the consolidation of links between genotype and phenotype. For example, advances in genomics allow us to unravel signatures of selection in mitochondrial and nuclear OXPHOS genes at population-relevant scales, while newly published complete atomic-resolution structures of the OXPHOS machinery enable more robust predictions of how these genes interact epistatically and co-evolutionarily. We use three case studies to show how integrative approaches have improved the understanding of mitonuclear interactions in OXPHOS, namely those driving high-altitude adaptation in bar-headed geese, allopatric population divergence in Tigriopus californicus copepods, and the genome architecture of nuclear genes coding for mitochondrial functions in the eastern yellow robin.

Mitochondrial-nuclear co-evolution leads to hybrid incompatibility through pentatricopeptide repeat proteins

Mitochondrial-nuclear incompatibility has a major role in reproductive isolation between species. However, the underlying mechanism and driving force of mitochondrial-nuclear incompatibility remain elusive. Here, we report a pentatricopeptide repeat-containing (PPR) protein, Ccm1, and its interacting partner, 15S rRNA, to be involved in hybrid incompatibility between two yeast species, Saccharomyces cerevisiae and Saccharomyces bayanus S. bayanus-Ccm1 has reduced binding affinity for S. cerevisiae-15S rRNA, leading to respiratory defects in hybrid cells. This incompatibility can be rescued by single mutations on several individual PPR motifs, demonstrating the highly evolvable nature of PPR proteins. When we examined other PPR proteins in the closely related Saccharomyces sensu stricto yeasts, about two-thirds of them showed detectable incompatibility. Our results suggest that fast co-evolution between flexible PPR proteins and their mitochondrial RNA substrates may be a common driving force in the development of mitochondrial-nuclear hybrid incompatibility.

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

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

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.

Chromosome-Wide Impacts on the Expression of Incompatibilities in Hybrids of Tigriopus californicus

Chromosome rearrangements such as inversions have been recognized previously as contributing to reproductive isolation by maintaining alleles together that jointly contribute to deleterious genetic interactions and postzygotic reproductive isolation. In this study, an impact of potential incompatibilities merely residing on the same chromosome was found in crosses of populations of the copepod Tigriopus californicus When genetically divergent populations of this copepod are crossed, hybrids show reduced fitness, and deviations from expected genotypic ratios can be used to determine regions of the genome involved in deleterious interactions. In this study, a set of markers was genotyped for a cross of two populations of T. californicus, and these markers show widespread deviations from Mendelian expectations, with entire chromosomes showing marked skew. Despite the importance of mtDNA/nuclear interactions in incompatibilities in this system in previous studies, in these crosses the expected patterns stemming from these interactions are not widely apparent. Females lack recombination in this species, and a striking difference is observed between male and female backcrosses. This suggests that the maintenance of multiple loci on individual chromosomes can enable some incompatibilities, perhaps playing a similar role in the initial rounds of hybridization to chromosomal rearrangements in preserving sets of alleles together that contribute to incompatibilities. Finally, it was observed that candidate pairs of incompatibility regions are not consistently interacting across replicates or subsets of these crosses, despite the repeatability of the deviations at many of the single loci themselves, suggesting that more complicated models of Dobzhansky-Muller incompatibilities may need to be considered.

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.

Fifty years of co-evolution and beyond: integrating co-evolution from molecules to species

Fifty years after Ehrlich and Raven's seminal paper, the idea of co-evolution continues to grow as a key concept in our understanding of organic evolution. This concept has not only provided a compelling synthesis between evolutionary biology and community ecology, but has also inspired research that extends beyond its original scope. In this article, we identify unresolved questions about the co-evolutionary process and advocate for the integration of co-evolutionary research from molecular to interspecific interactions. We address two basic questions: (i) What is co-evolution and how common is it? (ii) What is the unit of co-evolution? Both questions aim to explore the heart of the co-evolutionary process. Despite the claim that co-evolution is ubiquitous, we argue that there is in fact little evidence to support the view that reciprocal natural selection and coadaptation are common in nature. We also challenge the traditional view that co-evolution only occurs between traits of interacting species. Co-evolution has the potential to explain evolutionary processes and patterns that result from intra- and intermolecular biochemical interactions within cells, intergenomic interactions (e.g. nuclear-cytoplasmic) within species, as well as intergenomic interactions mediated by phenotypic traits between species. Research that bridges across these levels of organization will help to advance our understanding of the importance of the co-evolutionary processes in shaping the diversity of life on Earth.

Reproduction and development in Halocaridina rubra Holthuis, 1963 (Crustacea: Atyidae) clarifies larval ecology in the Hawaiian anchialine ecosystem

Larvae in aquatic habitats often develop in environments different from those they inhabit as adults. Shrimp in the Atyidae exemplify this trend, as larvae of many species require salt or brackish water for development, while adults are freshwater-adapted. An exception within the Atyidae family is the "anchialine clade," which are euryhaline as adults and endemic to habitats with subterranean fresh and marine water influences. Although the Hawaiian anchialine atyid Halocaridina rubra is a strong osmoregulator, its larvae have never been observed in nature. Moreover, larval development in anchialine species is poorly studied. Here, reproductive trends in laboratory colonies over a 5-y period are presented from seven genetic lineages and one mixed population of H. rubra; larval survivorship under varying salinities is also discussed. The presence and number of larvae differed significantly among lineages, with the mixed population being the most prolific. Statistical differences in reproduction attributable to seasonality also were identified. Larval survivorship was lowest (12% settlement rate) at a salinity approaching fresh water and significantly higher in brackish and seawater (88% and 72%, respectively). Correlated with this finding, identifiable gills capable of ion transport did not develop until metamorphosis into juveniles. Thus, early life stages of H. rubra are apparently excluded from surface waters, which are characterized by lower and fluctuating salinities. Instead, these stages are restricted to the subterranean (where there is higher and more stable salinity) portion of Hawaii's anchialine habitats due to their inability to tolerate low salinities. Taken together, these data contribute to the understudied area of larval ecology in the anchialine ecosystem.

Speciation and reduced hybrid female fertility in house mice

In mammals, intrinsic postzygotic isolation has been well studied in males but has been less studied in females, despite the fact that female gametogenesis and pregnancy provide arenas for hybrid sterility or inviability that are absent in males. Here, we asked whether inviability or sterility is observed in female hybrids of Mus musculus domesticus and M. m. musculus, taxa which hybridize in nature and for which male sterility has been well characterized. We looked for parent-of-origin growth phenotypes by measuring adult body weights in F1 hybrids. We evaluated hybrid female fertility by crossing F1 females to a tester male and comparing multiple reproductive parameters between intrasubspecific controls and intersubspecific hybrids. Hybrid females showed no evidence of parent-of-origin overgrowth or undergrowth, providing no evidence for reduced viability. However, hybrid females had smaller litter sizes, reduced embryo survival, fewer ovulations, and fewer small follicles relative to controls. Significant variation in reproductive parameters was seen among different hybrid genotypes, suggesting that hybrid incompatibilities are polymorphic within subspecies. Differences in reproductive phenotypes in reciprocal genotypes were observed and are consistent with cyto-nuclear incompatibilities or incompatibilities involving genomic imprinting. These findings highlight the potential importance of reduced hybrid female fertility in the early stages of speciation.

A Genetic Incompatibility Accelerates Adaptation in Yeast

During mismatch repair (MMR) MSH proteins bind to mismatches that form as the result of DNA replication errors and recruit MLH factors such as Mlh1-Pms1 to initiate excision and repair steps. Previously, we identified a negative epistatic interaction involving naturally occurring polymorphisms in the MLH1 and PMS1 genes of baker’s yeast. Here we hypothesize that a mutagenic state resulting from this negative epistatic interaction increases the likelihood of obtaining beneficial mutations that can promote adaptation to stress conditions. We tested this by stressing yeast strains bearing mutagenic (incompatible) and non-mutagenic (compatible) mismatch repair genotypes. Our data show that incompatible populations adapted more rapidly and without an apparent fitness cost to high salt stress. The fitness advantage of incompatible populations was rapid but disappeared over time. The fitness gains in both compatible and incompatible strains were due primarily to mutations in PMR1 that appeared earlier in incompatible evolving populations. These data demonstrate a rapid and reversible role (by mating) for genetic incompatibilities in accelerating adaptation in eukaryotes. They also provide an approach to link experimental studies to observational population genomics.

In nature, bacterial populations with high mutation rates can adapt faster to new environments by acquiring beneficial mutations. However, such populations also accumulate harmful mutations that reduce their fitness. We show that the model eukaryote baker’s yeast can use a similar mutator strategy to adapt to new environments. The mutator state that we observed resulted from an incompatibility involving two genes, MLH1 and PMS1, that work together to remove DNA replication errors through a spellchecking mismatch repair mechanism. This incompatibility can occur through mating between baker’s yeast from different genetic backgrounds, yielding mutator offspring containing an MLH1-PMS1 combination not present in either parent. Interestingly, these offspring adapted more rapidly to stress, compared to the parental strains, and did so without an overall loss in fitness. DNA sequencing analyses of baker’s yeast strains from across the globe support the presence of incompatible hybrid yeast strains in nature. These observations provide a powerful model to understand how the segregation of defects in DNA mismatch repair can serve as an effective strategy to enable eukaryotes to adapt to changing environments.

Reproductive isolation of hybrid populations driven by genetic incompatibilities

Despite its role in homogenizing populations, hybridization has also been proposed as a means to generate new species. The conceptual basis for this idea is that hybridization can result in novel phenotypes through recombination between the parental genomes, allowing a hybrid population to occupy ecological niches unavailable to parental species. Here we present an alternative model of the evolution of reproductive isolation in hybrid populations that occurs as a simple consequence of selection against genetic incompatibilities. Unlike previous models of hybrid speciation, our model does not incorporate inbreeding, or assume that hybrids have an ecological or reproductive fitness advantage relative to parental populations. We show that reproductive isolation between hybrids and parental species can evolve frequently and rapidly under this model, even in the presence of substantial ongoing immigration from parental species and strong selection against hybrids. An interesting prediction of our model is that replicate hybrid populations formed from the same pair of parental species can evolve reproductive isolation from each other. This non-adaptive process can therefore generate patterns of species diversity and relatedness that resemble an adaptive radiation. Intriguingly, several known hybrid species exhibit patterns of reproductive isolation consistent with the predictions of our model.

Mitochondrial genomes of the Baltic clam Macoma balthica (Bivalvia: Tellinidae): setting the stage for studying mito-nuclear incompatibilities

BACKGROUND

Allopatric divergence across lineages can lead to post-zygotic reproductive isolation upon secondary contact and disrupt coevolution between mitochondrial and nuclear genomes, promoting emergence of genetic incompatibilities. A previous F ST scan on the transcriptome of the Baltic clam Macoma balthica highlighted several genes potentially involved in mito-nuclear incompatibilities (MNIs). As proteins involved in the mitochondrial oxidative phosphorylation (OXPHO) chain are prone to MNIs and can contribute to the maintenance of genetic barriers, the mitochondrial genomes of six Ma. balthica individuals spanning two secondary contact zones were sequenced using the Illumina MiSeq plateform.

RESULTS

The mitogenome has an approximate length of 16,806 bp and encodes 13 protein-coding genes, 2 rRNAs and 22 tRNAs, all located on the same strand. atp8, a gene long reported as rare in bivalves, was detected. It encodes 42 amino acids and is putatively expressed and functional. A large unassigned region was identified between rrnS and tRNA (Met) and could likely correspond to the Control Region. Replacement and synonymous mutations were mapped on the inferred secondary structure of all protein-coding genes of the OXPHO chain. The atp6 and atp8 genes were characterized by background levels of replacement mutations, relative to synonymous mutations. However, most nad genes (notably nad2 and nad5) were characterized by an elevated proportion of replacement mutations.

CONCLUSIONS

Six nearly complete mitochondrial genomes were successfully assembled and annotated, providing the necessary roadmap to study MNIs at OXPHO loci. Few replacement mutations were mapped on mitochondrial-encoded ATP synthase subunits, which is in contrast with previous data on nuclear-encoded subunits. Conversely, the high population divergence and the prevalence of non-synonymous mutations at nad genes are congruent with previous observations from the nuclear transcriptome. This further suggest that MNIs between subunits of Complex I of the OXPHO chain, coding for NADH dehydrogenase, may play a role in maintaining barriers to gene flow in Ma. balthica.

Trace elements in major marketed marine bivalves from six northern coastal cities of China: concentrations and risk assessment for human health

One hundred and fifty nine samples of nine edible bivalve species (Argopecten irradians, Chlamys farreri, Crassostrea virginica, Lasaea nipponica, Meretrix meretrix, Mytilus edulis, Ruditapes philippinarum, Scapharca subcrenata and Sinonovacula constricta) were randomly collected from eight local seafood markets in six big cities (Dalian, Qingdao, Rizhao, Weifang, Weihai and Yantai) in the northern coastal areas of China for the investigation of trace element contamination. As, Cd, Cr, Cu, Hg, Pb and Zn were quantified. The risk of these trace elements to humans through bivalve consumption was then assessed. Results indicated that the concentrations of most of the studied trace element varied significantly with species: the average concentration of Cu in C. virginica was an order of magnitude higher than that in the remaining species; the average concentration of Zn was also highest in C. virginica; the average concentration of As, Cd and Pb was highest in R. philippinarum, C. farreri and A. irradians, respectively. Spatial differences in the concentrations of elements were generally less than those of interspecies, yet some elements such as Cr and Hg in the samples from different cities showed a significant difference in concentrations for some bivalve species. Trace element concentrations in edible tissues followed the order of Zn>Cu>As>Cd>Cr>Pb>Hg generally. Statistical analysis (one-way ANOVA) indicated that different species examined showed different bioaccumulation of trace elements. There were significant correlations between the concentrations of some elements. The calculated hazard quotients indicated in general that there was no obvious health risk from the intake of trace elements through bivalve consumption. But care must be taken considering the increasing amount of seafood consumption.

Mitonuclear interactions: evolutionary consequences over multiple biological scales

Fundamental biological processes hinge on coordinated interactions between genes spanning two obligate genomes--mitochondrial and nuclear. These interactions are key to complex life, and allelic variation that accumulates and persists at the loci embroiled in such intergenomic interactions should therefore be subjected to intense selection to maintain integrity of the mitochondrial electron transport system. Here, we compile evidence that suggests that mitochondrial-nuclear (mitonuclear) allelic interactions are evolutionarily significant modulators of the expression of key health-related and life-history phenotypes, across several biological scales--within species (intra- and interpopulational) and between species. We then introduce a new frontier for the study of mitonuclear interactions--those that occur within individuals, and are fuelled by the mtDNA heteroplasmy and the existence of nuclear-encoded mitochondrial gene duplicates and isoforms. Empirical evidence supports the idea of high-resolution tissue- and environment-specific modulation of intraindividual mitonuclear interactions. Predicting the penetrance, severity and expression patterns of mtDNA-induced mitochondrial diseases remains a conundrum. We contend that a deeper understanding of the dynamics and ramifications of mitonuclear interactions, across all biological levels, will provide key insights that tangibly advance our understanding, not only of core evolutionary processes, but also of the complex genetics underlying human mitochondrial disease.