A complete sequence of Saccharomyces paradoxus mitochondrial genome that restores the respiration in S. cerevisiae

Dissecting the sequence and structural determinants guiding m6A deposition and evolution via inter- and intra-species hybrids

BACKGROUND

N6-methyladenosine (m6A) is the most abundant mRNA modification, and controls mRNA stability. m6A distribution varies considerably between and within species. Yet, it is unclear to what extent this variability is driven by changes in genetic sequences ('cis') or cellular environments ('trans') and via which mechanisms.

RESULTS

Here we dissect the determinants governing RNA methylation via interspecies and intraspecies hybrids in yeast and mammalian systems, coupled with massively parallel reporter assays and m6A-QTL reanalysis. We find that m6A evolution and variability is driven primarily in 'cis', via two mechanisms: (1) variations altering m6A consensus motifs, and (2) variation impacting mRNA secondary structure. We establish that mutations impacting RNA structure - even when distant from an m6A consensus motif - causally dictate methylation propensity. Finally, we demonstrate that allele-specific differences in m6A levels lead to allele-specific changes in gene expression.

CONCLUSIONS

Our findings define the determinants governing m6A evolution and diversity and characterize the consequences thereof on gene expression regulation.

A high-quality reference genome for the fission yeast Schizosaccharomyces osmophilus

Fission yeasts are an ancient group of fungal species that diverged from each other from tens to hundreds of million years ago. Among them is the preeminent model organism Schizosaccharomyces pombe, which has significantly contributed to our understandings of molecular mechanisms underlying fundamental cellular processes. The availability of the genomes of S. pombe and 3 other fission yeast species S. japonicus, S. octosporus, and S. cryophilus has enabled cross-species comparisons that provide insights into the evolution of genes, pathways, and genomes. Here, we performed genome sequencing on the type strain of the recently identified fission yeast species S. osmophilus and obtained a complete mitochondrial genome and a nuclear genome assembly with gaps only at rRNA gene arrays. A total of 5,098 protein-coding nuclear genes were annotated and orthologs for more than 95% of them were identified. Genome-based phylogenetic analysis showed that S. osmophilus is most closely related to S. octosporus and these 2 species diverged around 16 million years ago. To demonstrate the utility of this S. osmophilus reference genome, we conducted cross-species comparative analyses of centromeres, telomeres, transposons, the mating-type region, Cbp1 family proteins, and mitochondrial genomes. These analyses revealed conservation of repeat arrangements and sequence motifs in centromere cores, identified telomeric sequences composed of 2 types of repeats, delineated relationships among Tf1/sushi group retrotransposons, characterized the evolutionary origins and trajectories of Cbp1 family domesticated transposases, and discovered signs of interspecific transfer of 2 types of mitochondrial selfish elements.

Phylogenetic and ecological reevaluation of the order Onygenales

The order Onygenales is classified in the class Eurotiomycetes of the subphylum Pezizomycotina. Families in this order have classically been isolated from soil and dung, and two lineages contain causative agents of superficial, cutaneous and systemic infections in mammals. The ecology and habitat choices of the species are driven mainly by the keratin and cellulose degradation abilities. The present study aimed to investigate whether the ecological trends of the members of Onygenales can be interpreted in an evolutionary sense, linking phylogenetic parameters with habitat preferences, to achieve polyphasic definitions of the main taxonomic groups. Evolutionary processes were estimated by multiple gene genealogies and divergence time analysis. Previously described families, namely, Arthrodermataceae, Ajellomycetaceae, Ascosphaeraceae, Eremascaceae, Gymnoascaceae, Onygenaceae and Spiromastigoidaceae, were accepted in Onygenales, and two new families, Malbrancheaceae and Neogymnomycetaceae, were introduced. A number of species could not be assigned to any of the defined families. Our study provides a revised overview of the main lines of taxonomy of Onygenales, supported by multilocus analyses of ITS, LSU, TUB, TEF1, TEF3, RPB1, RPB2, and ribosomal protein 60S L10 (L1) (RP60S) sequences, combined with available data on ecology, physiology, morphology, and genomics.

Hybridization drives mitochondrial DNA degeneration and metabolic shift in a species with biparental mitochondrial inheritance

Mitochondrial DNA (mtDNA) is a cytoplasmic genome that is essential for respiratory metabolism. Although uniparental mtDNA inheritance is most common in animals and plants, distinct mtDNA haplotypes can coexist in a state of heteroplasmy, either because of paternal leakage or de novo mutations. mtDNA integrity and the resolution of heteroplasmy have important implications, notably for mitochondrial genetic disorders, speciation, and genome evolution in hybrids. However, the impact of genetic variation on the transition to homoplasmy from initially heteroplasmic backgrounds remains largely unknown. Here, we use Saccharomyces yeasts, fungi with constitutive biparental mtDNA inheritance, to investigate the resolution of mtDNA heteroplasmy in a variety of hybrid genotypes. We previously designed 11 crosses along a gradient of parental evolutionary divergence using undomesticated isolates of Saccharomyces paradoxus and Saccharomyces cerevisiae Each cross was independently replicated 48 to 96 times, and the resulting 864 hybrids were evolved under relaxed selection for mitochondrial function. Genome sequencing of 446 MA lines revealed extensive mtDNA recombination, but the recombination rate was not predicted by parental divergence level. We found a strong positive relationship between parental divergence and the rate of large-scale mtDNA deletions, which led to the loss of respiratory metabolism. We also uncovered associations between mtDNA recombination, mtDNA deletion, and genome instability that were genotype specific. Our results show that hybridization in yeast induces mtDNA degeneration through large-scale deletion and loss of function, with deep consequences for mtDNA evolution, metabolism, and the emergence of reproductive isolation.

Mitochondrial DNA duplication, recombination, and introgression during interspecific hybridization

mtDNA recombination events in yeasts are known, but altered mitochondrial genomes were not completed. Therefore, we analyzed recombined mtDNAs in six Saccharomyces cerevisiae × Saccharomyces paradoxus hybrids in detail. Assembled molecules contain mostly segments with variable length introgressed to other mtDNA. All recombination sites are in the vicinity of the mobile elements, introns in cox1, cob genes and free standing ORF1, ORF4. The transplaced regions involve co-converted proximal exon regions. Thus, these selfish elements are beneficial to the host if the mother molecule is challenged with another molecule for transmission to the progeny. They trigger mtDNA recombination ensuring the transfer of adjacent regions, into the progeny of recombinant molecules. The recombination of the large segments may result in mitotically stable duplication of several genes.

sppIDer: A Species Identification Tool to Investigate Hybrid Genomes with High-Throughput Sequencing

The genomics era has expanded our knowledge about the diversity of the living world, yet harnessing high-throughput sequencing data to investigate alternative evolutionary trajectories, such as hybridization, is still challenging. Here we present sppIDer, a pipeline for the characterization of interspecies hybrids and pure species, that illuminates the complete composition of genomes. sppIDer maps short-read sequencing data to a combination genome built from reference genomes of several species of interest and assesses the genomic contribution and relative ploidy of each parental species, producing a series of colorful graphical outputs ready for publication. As a proof-of-concept, we use the genus Saccharomyces to detect and visualize both interspecies hybrids and pure strains, even with missing parental reference genomes. Through simulation, we show that sppIDer is robust to variable reference genome qualities and performs well with low-coverage data. We further demonstrate the power of this approach in plants, animals, and other fungi. sppIDer is robust to many different inputs and provides visually intuitive insight into genome composition that enables the rapid identification of species and their interspecies hybrids. sppIDer exists as a Docker image, which is a reusable, reproducible, transparent, and simple-to-run package that automates the pipeline and installation of the required dependencies (https://github.com/GLBRC/sppIDer; last accessed September 6, 2018).

Mitochondrial Recombination and Introgression during Speciation by Hybridization

Genome recombination is a major source of genotypic diversity and contributes to adaptation and speciation following interspecies hybridization. The contribution of recombination in these processes has been thought to be largely limited to the nuclear genome because organelles are mostly uniparentally inherited in animals and plants, which prevents recombination. Unicellular eukaryotes such as budding yeasts do, however, transmit mitochondria biparentally, suggesting that during hybridization, both parents could provide alleles that contribute to mitochondrial functions such as respiration and metabolism in hybrid populations or hybrid species. We examined the dynamics of mitochondrial genome transmission and evolution during speciation by hybridization in the natural budding yeast Saccharomyces paradoxus. Using population-scale mitochondrial genome sequencing in two endemic North American incipient species SpB and SpC and their hybrid species SpC*, we found that both parental species contributed to the hybrid mitochondrial genome through recombination. We support our findings by showing that mitochondrial recombination between parental types is frequent in experimental crosses that recreate the early step of this speciation event. In these artificial hybrids, we observed that mitochondrial genome recombination enhances phenotypic variation among diploid hybrids, suggesting that it could play a role in the phenotypic differentiation of hybrid species. Like the nuclear genome, the mitochondrial genome can, therefore, also play a role in hybrid speciation.

The evolutionary history of Saccharomyces species inferred from completed mitochondrial genomes and revision in the 'yeast mitochondrial genetic code'

The yeast Saccharomyces are widely used to test ecological and evolutionary hypotheses. A large number of nuclear genomic DNA sequences are available, but mitochondrial genomic data are insufficient. We completed mitochondrial DNA (mtDNA) sequencing from Illumina MiSeq reads for all Saccharomyces species. All are circularly mapped molecules decreasing in size with phylogenetic distance from Saccharomyces cerevisiae but with similar gene content including regulatory and selfish elements like origins of replication, introns, free-standing open reading frames or GC clusters. Their most profound feature is species-specific alteration in gene order. The genetic code slightly differs from well-established yeast mitochondrial code as GUG is used rarely as the translation start and CGA and CGC code for arginine. The multilocus phylogeny, inferred from mtDNA, does not correlate with the trees derived from nuclear genes. mtDNA data demonstrate that Saccharomyces cariocanus should be assigned as a separate species and Saccharomyces bayanus CBS 380T should not be considered as a distinct species due to mtDNA nearly identical to Saccharomyces uvarum mtDNA. Apparently, comparison of mtDNAs should not be neglected in genomic studies as it is an important tool to understand the origin and evolutionary history of some yeast species.

Mitochondrial genome evolution in the Saccharomyces sensu stricto complex

Exploring the evolutionary patterns of mitochondrial genomes is important for our understanding of the Saccharomyces sensu stricto (SSS) group, which is a model system for genomic evolution and ecological analysis. In this study, we first obtained the complete mitochondrial sequences of two important species, Saccharomyces mikatae and Saccharomyces kudriavzevii. We then compared the mitochondrial genomes in the SSS group with those of close relatives, and found that the non-coding regions evolved rapidly, including dramatic expansion of intergenic regions, fast evolution of introns and almost 20-fold higher rearrangement rates than those of the nuclear genomes. However, the coding regions, and especially the protein-coding genes, are more conserved than those in the nuclear genomes of the SSS group. The different evolutionary patterns of coding and non-coding regions in the mitochondrial and nuclear genomes may be related to the origin of the aerobic fermentation lifestyle in this group. Our analysis thus provides novel insights into the evolution of mitochondrial genomes.

Mitochondrial introgression suggests extensive ancestral hybridization events among Saccharomyces species

Horizontal gene transfer (HGT) in eukaryotic plastids and mitochondrial genomes is common, and plays an important role in organism evolution. In yeasts, recent mitochondrial HGT has been suggested between S. cerevisiae and S. paradoxus. However, few strains have been explored given the lack of accurate mitochondrial genome annotations. Mitochondrial genome sequences are important to understand how frequent these introgressions occur, and their role in cytonuclear incompatibilities and fitness. Indeed, most of the Bateson-Dobzhansky-Muller genetic incompatibilities described in yeasts are driven by cytonuclear incompatibilities. We herein explored the mitochondrial inheritance of several worldwide distributed wild Saccharomyces species and their hybrids isolated from different sources and geographic origins. We demonstrated the existence of several recombination points in mitochondrial region COX2-ORF1, likely mediated by either the activity of the protein encoded by the ORF1 (F-SceIII) gene, a free-standing homing endonuclease, or mostly facilitated by A+T tandem repeats and regions of integration of GC clusters. These introgressions were shown to occur among strains of the same species and among strains of different species, which suggests a complex model of Saccharomyces evolution that involves several ancestral hybridization events in wild environments.

The Genome Sequence of Saccharomyces eubayanus and the Domestication of Lager-Brewing Yeasts

The dramatic phenotypic changes that occur in organisms during domestication leave indelible imprints on their genomes. Although many domesticated plants and animals have been systematically compared with their wild genetic stocks, the molecular and genomic processes underlying fungal domestication have received less attention. Here, we present a nearly complete genome assembly for the recently described yeast species Saccharomyces eubayanus and compare it to the genomes of multiple domesticated alloploid hybrids of S. eubayanus × S. cerevisiae (S. pastorianus syn. S. carlsbergensis), which are used to brew lager-style beers. We find that the S. eubayanus subgenomes of lager-brewing yeasts have experienced increased rates of evolution since hybridization, and that certain genes involved in metabolism may have been particularly affected. Interestingly, the S. eubayanus subgenome underwent an especially strong shift in selection regimes, consistent with more extensive domestication of the S. cerevisiae parent prior to hybridization. In contrast to recent proposals that lager-brewing yeasts were domesticated following a single hybridization event, the radically different neutral site divergences between the subgenomes of the two major lager yeast lineages strongly favor at least two independent origins for the S. cerevisiae × S. eubayanus hybrids that brew lager beers. Our findings demonstrate how this industrially important hybrid has been domesticated along similar evolutionary trajectories on multiple occasions.

Population structure of mitochondrial genomes in Saccharomyces cerevisiae

Background

Rigorous study of mitochondrial functions and cell biology in the budding yeast, Saccharomyces cerevisiae has advanced our understanding of mitochondrial genetics. This yeast is now a powerful model for population genetics, owing to large genetic diversity and highly structured populations among wild isolates. Comparative mitochondrial genomic analyses between yeast species have revealed broad evolutionary changes in genome organization and architecture. A fine-scale view of recent evolutionary changes within S. cerevisiae has not been possible due to low numbers of complete mitochondrial sequences.

Results

To address challenges of sequencing AT-rich and repetitive mitochondrial DNAs (mtDNAs), we sequenced two divergent S. cerevisiae mtDNAs using a single-molecule sequencing platform (PacBio RS). Using de novo assemblies, we generated highly accurate complete mtDNA sequences. These mtDNA sequences were compared with 98 additional mtDNA sequences gathered from various published collections. Phylogenies based on mitochondrial coding sequences and intron profiles revealed that intraspecific diversity in mitochondrial genomes generally recapitulated the population structure of nuclear genomes. Analysis of intergenic sequence indicated a recent expansion of mobile elements in certain populations. Additionally, our analyses revealed that certain populations lacked introns previously believed conserved throughout the species, as well as the presence of introns never before reported in S. cerevisiae.

Conclusions

Our results revealed that the extensive variation in S. cerevisiae mtDNAs is often population specific, thus offering a window into the recent evolutionary processes shaping these genomes. In addition, we offer an effective strategy for sequencing these challenging AT-rich mitochondrial genomes for small scale projects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1664-4) contains supplementary material, which is available to authorized users.

Post-zygotic sterility and cytonuclear compatibility limits in S. cerevisiae xenomitochondrial cybrids

Nucleo-mitochondrial interactions, particularly those determining the primary divergence of biological species, can be studied by means of xenomitochondrial cybrids, which are cells where the original mitochondria are substituted by their counterparts from related species. Saccharomyces cerevisiae cybrids are prepared simply by the mating of the ρ(0) strain with impaired karyogamy and germinating spores from other Saccharomyces species and fall into three categories. Cybrids with compatible mitochondrial DNA (mtDNA) from Saccharomyces paradoxus CBS 432 and Saccharomyces cariocanus CBS 7994 are metabolically and genetically similar to cybrids containing mtDNA from various S. cerevisiae. Cybrids with mtDNA from other S. paradoxus strains, S. cariocanus, Saccharomyces kudriavzevii, and Saccharomyces mikatae require a period of adaptation to establish efficient oxidative phosphorylation. They exhibit a temperature-sensitive phenotype, slower growth rate on a non-fermentable carbon source and a long lag phase after the shift from glucose. Their decreased respiration capacity and reduced cytochrome aa3 content is associated with the inefficient splicing of cox1I3β, the intron found in all Saccharomyces species but not in S. cerevisiae. The splicing defect is compensated in cybrids by nuclear gain-of-function and can be alternatively suppressed by overexpression of MRP13 gene for mitochondrial ribosomal protein or the MRS2, MRS3, and MRS4 genes involved in intron splicing. S. cerevisiae with Saccharomyces bayanus mtDNA is unable to respire and the growth on ethanol-glycerol can be restored only after mating to some mit (-) strains. The nucleo-mitochondrial compatibility limit of S. cerevisiae and other Saccharomyces was set between S. kudriavzevii and S. bayanus at the divergence from S. cerevisiae about 15 MYA. The MRS1-cox1 S. cerevisiae/S. paradoxus cytonuclear Dobzhansky-Muller pair has a neglible impact on the separation of species since its imperfection is compensated for by gain-of-function mutation.

Mitochondrial-nuclear epistasis contributes to phenotypic variation and coadaptation in natural isolates of Saccharomyces cerevisiae

Mitochondria are essential multifunctional organelles whose metabolic functions, biogenesis, and maintenance are controlled through genetic interactions between mitochondrial and nuclear genomes. In natural populations, mitochondrial efficiencies may be impacted by epistatic interactions between naturally segregating genome variants. The extent that mitochondrial-nuclear epistasis contributes to the phenotypic variation present in nature is unknown. We have systematically replaced mitochondrial DNAs in a collection of divergent Saccharomyces cerevisiae yeast isolates and quantified the effects on growth rates in a variety of environments. We found that mitochondrial-nuclear interactions significantly affected growth rates and explained a substantial proportion of the phenotypic variances under some environmental conditions. Naturally occurring mitochondrial-nuclear genome combinations were more likely to provide growth advantages, but genetic distance could not predict the effects of epistasis. Interruption of naturally occurring mitochondrial-nuclear genome combinations increased endogenous reactive oxygen species in several strains to levels that were not always proportional to growth rate differences. Our results demonstrate that interactions between mitochondrial and nuclear genomes generate phenotypic diversity in natural populations of yeasts and that coadaptation of intergenomic interactions likely occurs quickly within the specific niches that yeast occupy. This study reveals the importance of considering allelic interactions between mitochondrial and nuclear genomes when investigating evolutionary relationships and mapping the genetic basis underlying complex traits.

A high-definition view of functional genetic variation from natural yeast genomes

The question of how genetic variation in a population influences phenotypic variation and evolution is of major importance in modern biology. Yet much is still unknown about the relative functional importance of different forms of genome variation and how they are shaped by evolutionary processes. Here we address these questions by population level sequencing of 42 strains from the budding yeast Saccharomyces cerevisiae and its closest relative S. paradoxus. We find that genome content variation, in the form of presence or absence as well as copy number of genetic material, is higher within S. cerevisiae than within S. paradoxus, despite genetic distances as measured in single-nucleotide polymorphisms being vastly smaller within the former species. This genome content variation, as well as loss-of-function variation in the form of premature stop codons and frameshifting indels, is heavily enriched in the subtelomeres, strongly reinforcing the relevance of these regions to functional evolution. Genes affected by these likely functional forms of variation are enriched for functions mediating interaction with the external environment (sugar transport and metabolism, flocculation, metal transport, and metabolism). Our results and analyses provide a comprehensive view of genomic diversity in budding yeast and expose surprising and pronounced differences between the variation within S. cerevisiae and that within S. paradoxus. We also believe that the sequence data and de novo assemblies will constitute a useful resource for further evolutionary and population genomics studies.

The mitochondrial genome impacts respiration but not fermentation in interspecific Saccharomyces hybrids

In eukaryotes, mitochondrial DNA (mtDNA) has high rate of nucleotide substitution leading to different mitochondrial haplotypes called mitotypes. However, the impact of mitochondrial genetic variant on phenotypic variation has been poorly considered in microorganisms because mtDNA encodes very few genes compared to nuclear DNA, and also because mitochondrial inheritance is not uniparental. Here we propose original material to unravel mitotype impact on phenotype: we produced interspecific hybrids between S. cerevisiae and S. uvarum species, using fully homozygous diploid parental strains. For two different interspecific crosses involving different parental strains, we recovered 10 independent hybrids per cross, and allowed mtDNA fixation after around 80 generations. We developed PCR-based markers for the rapid discrimination of S. cerevisiae and S. uvarum mitochondrial DNA. For both crosses, we were able to isolate fully isogenic hybrids at the nuclear level, yet possessing either S. cerevisiae mtDNA (Sc-mtDNA) or S. uvarum mtDNA (Su-mtDNA). Under fermentative conditions, the mitotype has no phenotypic impact on fermentation kinetics and products, which was expected since mtDNA are not necessary for fermentative metabolism. Alternatively, under respiratory conditions, hybrids with Sc-mtDNA have higher population growth performance, associated with higher respiratory rate. Indeed, far from the hypothesis that mtDNA variation is neutral, our work shows that mitochondrial polymorphism can have a strong impact on fitness components and hence on the evolutionary fate of the yeast populations. We hypothesize that under fermentative conditions, hybrids may fix stochastically one or the other mt-DNA, while respiratory environments may increase the probability to fix Sc-mtDNA.