Mapping mitonuclear epistasis using a novel recombinant yeast population

Mitochondrial genome diversity across the subphylum Saccharomycotina

Introduction

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available.

Methods

By extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades.

Results

Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements.

Discussion

As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.

Diet and mitonuclear haplotype interactions affect growth rate in a slime mould

Trait expression in metazoans is strongly influenced by the balance of macronutrients (i.e. protein, carbohydrate and fat) in the diet. At the same time, an individual's genetic background seems to regulate the magnitude of phenotypic response to a particular diet. It needs to be better understood whether interactions between diet, genetic background and trait expression are found in unicellular eukaryotes. A protist-the slime mould, Physarum polycephalum can choose diets based on protein-to-carbohydrate (P:C) content to support optimal growth rate. Yet, the role of genetic background (variation in the mitochondrial and nuclear DNAs) in mediating growth rate response to dietary P:C ratios in the slime mould is unknown. Here, we studied the effects of interactions between mitochondrial and nuclear DNA haplotypes and diet (i.e. G × G × E interactions) on the growth rate of P. polycephalum. A genetic panel of six distinct strains of P. polycephalum that differ in their mitochondrial and nuclear DNA haplotypes was used to measure growth rate across five diets that varied in their P:C ratio and total calories. We first determined the strains' growth rate (total biomass and surface area) when grown on a set menu with access to a particular diet. We then assessed whether the growth rate of strains increased on a buffet menu with access to all diets. Our findings show that the growth rate of P. polycephalum is generally higher on diets containing more carbohydrates than protein and that total calories negatively affect the growth rate. Three-way interactions between mitochondrial, nuclear haplotypes and dietary P:C ratios affected the strains' surface area of growth but not biomass. Intriguingly, strains did not increase their surface area and biomass when they had access to all diets on the buffet menu. Our findings have broad implications for our understanding of the effect of mitonuclear interactions on trait expression across diverse eukaryotic lineages.

Mitochondrial Genome Diversity across the Subphylum Saccharomycotina

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available. By extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades. Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements. As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.