Intraspecific Diversity of Fission Yeast Mitochondrial Genomes

Schizosaccharomyces pombe as a fundamental model for research on mitochondrial gene expression: Progress, achievements and outlooks

Schizosaccharomyces pombe (fission yeast) is an attractive model for mitochondrial research. The organism resembles human cells in terms of mitochondrial inheritance, mitochondrial transport, sugar metabolism, mitogenome structure and dependence of viability on the mitogenome (the petite-negative phenotype). Transcriptions of these genomes produce only a few polycistronic transcripts, which then undergo processing as per the tRNA punctuation model. In general, the machinery for mitochondrial gene expression is structurally and functionally conserved between fission yeast and humans. Furthermore, molecular research on S. pombe is supported by a considerable number of experimental techniques and database resources. Owing to these advantages, fission yeast has significantly contributed to biomedical and fundamental research. Here, we review the current state of knowledge regarding S. pombe mitochondrial gene expression, and emphasise the pertinence of fission yeast as both a model and tool, especially for studies on mitochondrial translation.

Reinstatement of the fission yeast species Schizosaccharomyces versatilis Wickerham et Duprat, a sibling species of Schizosaccharomyces japonicus

Schizosaccharomyces japonicus Yukawa et Maki (1931) and Schizosaccharomyces versatilis Wickerham et Duprat (1945) have been treated as varieties of S. japonicus or as conspecific, based on various approaches including mating trials and nDNA/nDNA optical reassociation studies. However, the type strains of S. japonicus and S. versatilis differ by five substitutions (99.15% identity) and one 1-bp indel in the sequences of the D1/D2 domain of the 26S rRNA gene, and 23 substitutions (96.3% identity) and 31-bp indels in the sequences of internal transcribed spacer (ITS) of rRNA, suggesting that they may not be conspecific. To reassess their taxonomic status, we conducted mating trials and whole-genome analyses. Mating trials using the type strains showed a strong but incomplete prezygotic sterility barrier, yielding interspecies mating products at two orders of magnitude lower efficiency than intraspecies matings. These mating products, which were exclusively allodiploid hybrids, were unable to undergo the haplontic life cycle of the parents. We generated chromosome-level gap-less genome assemblies for both type strains. Whole genome sequences yielded an average nucleotide identity (ANI) of 86.4%, indicating clear separation of S. japonicus and S. versatilis. Based on these findings, we propose the reinstatement of S. versatilis as a distinct species (holotype strain: CBS 103T and ex-types: NRRL Y-1026, NBRC 1607, ATCC 9987, PYCC 7100; Mycobank no.: 847838).

Diverse modes of chromosome terminal deletion in spontaneous canavanine-resistant Schizosaccharomyces pombe mutants

Canavanine resistance has been used to analyze mutation rates in the fission yeast Schizosaccharomyces pombe . However, the genetic basis of canavanine resistance in this organism remains incompletely understood. Here, we performed whole genome sequencing on five spontaneously arising canavanine-resistant S. pombe mutants, including the can2-1 mutant isolated in the 1970s. This analysis revealed that three mutants, including can2-1 , experienced terminal deletions of the left arm of chromosome II, leading to the loss of multiple amino acid transporter genes. Interestingly, these three mutants underwent chromosome terminal deletion through distinct mechanisms, including homology-driven translocation, homology-independent chromosome fusion, and de novo telomere addition. Our findings shed new light on the genetic basis of canavanine resistance and mechanisms underlying chromosome terminal deletions in fission yeast.

Schizosaccharomyces lindneri sp. nov., a fission yeast occurring in honey

Two strains of fission yeast were isolated from honey. They differ from the type strain of Schizosaccharomyces octosporus by three substitutions in the D1/D2 domain of the nuclear 26S large subunit ribosomal RNA (rRNA) gene sequence, resulting in a 99.5% identity. In the internal transcribed spacer (ITS) region (consisting of ITS1, 5.8S rDNA, and ITS2), the strains differ from S. octosporus by 16 gaps and 91 substitutions, which is equivalent to an identity of 88.1%. Genome sequencing on one of the new strains revealed that the average nucleotide identity (ANI) between its genome and the reference genome of S. octosporus is 90.43% and there exist major genome rearrangements between the two genomes. Mating analysis revealed that S. octosporus and one of the new strains are completely reproductively separated. A strong prezygotic barrier exists and the few mating products consist of diploid hybrids that do not form recombinant ascospores. In the new strains, asci are either zygotic, arising from conjugation, or they develop without conjugation from asexual cells (azygotic). Compared to the currently recognized Schizosaccharomyces species, the spectrum of nutrients that are assimilated by the new strains is restricted. Of the 43 carbohydrates that were included in the physiological standard tests, only 7 were assimilated. According to the results of the genome sequence analysis, the mating trials, and the phenotypic characterization, the new species Schizosaccharomyces lindneri is described to accommodate the two strains (holotype: CBS 18203^T  and ex-type: MUCL 58363; MycoBank no.: MB 847838).

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.

Insight into the evolutionary and domesticated history of the most widely cultivated mushroom Agaricus bisporus via mitogenome sequences of 361 global strains

Agaricus bisporus is the most widely cultivated edible mushroom in the world with a only around three hundred years known history of cultivation. Therefore, it represents an ideal organism not only to investigate the natural evolutionary history but also the understanding on the evolution going back to the early era of domestication. In this study, we generated the mitochondrial genome sequences of 352 A. bisporus strains and 9 strains from 4 closely related species around the world. The population mitogenomic study revealed all A. bisporus strains can be divided into seven clades, and all domesticated cultivars present only in two of those clades. The molecular dating analysis showed this species origin in Europe on 4.6 Ma and we proposed the main dispersal routes. The detailed mitogenome structure studies showed that the insertion of the plasmid-derived dpo gene caused a long fragment (MIR) inversion, and the distributions of the fragments of dpo gene were strictly in correspondence with these seven clades. Our studies also showed A. bisporus population contains 30 intron distribution patterns (IDPs), while all cultivars contain only two IDPs, which clearly exhibit intron loss compared to the others. Either the loss occurred before or after domestication, that could suggest that the change facilitates their adaptation to the cultivated environment.

S. pombe wtf drivers use dual transcriptional regulation and selective protein exclusion from spores to cause meiotic drive

Meiotic drivers bias gametogenesis to ensure their transmission into more than half the offspring of a heterozygote. In Schizosaccharomyces pombe, wtf meiotic drivers destroy the meiotic products (spores) that do not inherit the driver from a heterozygote, thereby reducing fertility. wtf drivers encode both a Wtfpoison protein and a Wtfantidote protein using alternative transcriptional start sites. Here, we analyze how the expression and localization of the Wtf proteins are regulated to achieve drive. We show that transcriptional timing and selective protein exclusion from developing spores ensure that all spores are exposed to Wtf4poison, but only the spores that inherit wtf4 receive a dose of Wtf4antidote sufficient for survival. In addition, we show that the Mei4 transcription factor, a master regulator of meiosis, controls the expression of the wtf4poison transcript. This transcriptional regulation, which includes the use of a critical meiotic transcription factor, likely complicates the universal suppression of wtf genes without concomitantly disrupting spore viability. We propose that these features contribute to the evolutionary success of the wtf drivers.

Insights into the ecology of Schizosaccharomyces species in natural and artificial habitats

The fission yeast genus Schizosaccharomyces contains important model organisms for biological research. In particular, S. pombe is a widely used model eukaryote. So far little is known about the natural and artificial habitats of species in this genus. Finding out where S. pombe and other fission yeast species occur and how they live in their habitats can promote better understanding of their biology. Here we investigate in which substrates S. pombe, S. octosporus, S. osmophilus and S. japonicus are present. To this end about 2100 samples consisting of soil, tree sap fluxes, fresh fruit, dried fruit, honey, cacao beans, molasses and other substrates were analyzed. Effective isolation methods that allow efficient isolation of the above mentioned species were developed. Based on the frequency of isolating different fission yeast species in various substrates and on extensive literature survey, conclusions are drawn on their ecology. The results suggest that the primary habitat of S. pombe and S. octosporus is honeybee honey. Both species were also frequently detected on certain dried fruit like raisins, mango or pineapple to which they could be brought by the honey bees during ripening or during drying. While S. pombe was regularly isolated from grape mash and from fermented raw cacao beans S. octosporus was never isolated from fresh fruit. The main habitat of S. osmophilus seems to be solitary bee beebread. It was rarely isolated from raisins. S. japonicus was mainly found in forest substrates although it occurs on fruit and in fruit fermentations, too.

Canavanine resistance mutation can1-1 in Schizosaccharomyces pombe is a missense mutation in the ubiquitin ligase adaptor gene any1

In Schizosaccharomyces pombe, the can1-1 mutation confers resistance to the toxic arginine analog canavanine. This mutation has been assumed to disrupt a gene encoding an arginine transporter. In PomBase, the gene SPBC18H10.16 is currently designated can1. Here, we sequenced the genomes of three can1-1 strains. No mutations were found in SPBC18H10.16. Instead, these strains harbor an R175C mutation in the gene any1 (SPBC18H10.20c). any1 encodes an α-arrestin that acts as a ubiquitin ligase adaptor to downregulate plasma membrane amino acid transporters. Our findings indicate that can1-1 is not a loss-of-function mutation in an amino acid transporter gene, but a possible gain-of-function mutation in a gene encoding a negative regulator of amino acid transporters.

The search for Schizosaccharomyces fission yeasts in environmental meta-transcriptomes

Fission yeast is an important model organism in evolutionary genetics and cell biology research. Nevertheless, most research is limited to a single laboratory strain and knowledge of its natural occurrence is limited, which reduces our understanding of its life history and hinders isolation of new strains from nature. Understanding the natural diversity of fission yeast can provide insight into its genetic and phenotypic diversity and the evolutionary processes that shaped these. Here we aimed to identify candidate natural habitats of fission yeasts by searching through a large collection of publicly available environmental metatranscriptomic datasets. Using a custom pipeline, we processed over 13,000 NCBI SRA accessions, from a wide range of 34 different environmental categories. Overall, we found a very low abundance of putative yeast transcripts, with most fission yeast signatures coming from the categories of 'food' and 'terrestrial arthropods'. Additionally, a signal could be found in a variety of marine and fresh aquatic habitats. Our results do not provide a conclusive answer on the natural habitat of fission yeasts, but our analysis further narrows the range of locations where fission yeasts naturally occur.

Reactivation of transposable elements following hybridization in fission yeast

Hybridization is thought to reactivate transposable elements (TEs) that were efficiently suppressed in the genomes of the parental hosts. Here, we provide evidence for this “genomic shock hypothesis” in the fission yeast Schizosaccharomyces pombe. In this species, two divergent lineages (Sp and Sk) have experienced recent, likely human-induced, hybridization. We used long-read sequencing data to assemble genomes of 37 samples derived from 31 S. pombe strains spanning a wide range of ancestral admixture proportions. A comprehensive TE inventory revealed exclusive presence of long terminal repeat (LTR) retrotransposons. Sequence analysis of active full-length elements, as well as solo LTRs, revealed a complex history of homologous recombination. Population genetic analyses of syntenic sequences placed insertion of many solo LTRs before the split of the Sp and Sk lineages. Most full-length elements were inserted more recently, after hybridization. With the exception of a single full-length element with signs of positive selection, both solo LTRs and, in particular, full-length elements carry signatures of purifying selection indicating effective removal by the host. Consistent with reactivation upon hybridization, the number of full-length LTR retrotransposons, varying extensively from zero to 87 among strains, significantly increases with the degree of genomic admixture. This study gives a detailed account of global TE diversity in S. pombe, documents complex recombination histories within TE elements, and provides evidence for the “genomic shock hypothesis.”

Diverse mating phenotypes impact the spread of wtf meiotic drivers in Schizosaccharomyces pombe

Meiotic drivers are genetic elements that break Mendel's law of segregation to be transmitted into more than half of the offspring produced by a heterozygote. The success of a driver relies on outcrossing (mating between individuals from distinct lineages) because drivers gain their advantage in heterozygotes. It is, therefore, curious that Schizosaccharomyces pombe, a species reported to rarely outcross, harbors many meiotic drivers. To address this paradox, we measured mating phenotypes in S. pombe natural isolates. We found that the propensity for cells from distinct clonal lineages to mate varies between natural isolates and can be affected both by cell density and by the available sexual partners. Additionally, we found that the observed levels of preferential mating between cells from the same clonal lineage can slow, but not prevent, the spread of a wtf meiotic driver in the absence of additional fitness costs linked to the driver. These analyses reveal parameters critical to understanding the evolution of S. pombe and help explain the success of meiotic drivers in this species.

Genome Comparisons of the Fission Yeasts Reveal Ancient Collinear Loci Maintained by Natural Selection

Fission yeasts have a unique life history and exhibit distinct evolutionary patterns from other yeasts. Besides, the species demonstrate stable genome structures despite the relatively fast evolution of their genomic sequences. To reveal what could be the reason for that, comparative genomic analyses were carried out. Our results provided evidence that the structural and sequence evolution of the fission yeasts were correlated. Moreover, we revealed ancestral locally collinear blocks (aLCBs), which could have been inherited from their last common ancestor. These aLCBs proved to be the most conserved regions of the genomes as the aLCBs contain almost eight genes/blocks on average in the same orientation and order across the species. Gene order of the aLCBs is mainly fission-yeast-specific but supports the idea of filamentous ancestors. Nevertheless, the sequences and gene structures within the aLCBs are as mutable as any sequences in other parts of the genomes. Although genes of certain Gene Ontology (GO) categories tend to cluster at the aLCBs, those GO enrichments are not related to biological functions or high co-expression rates, they are, rather, determined by the density of essential genes and Rec12 cleavage sites. These data and our simulations indicated that aLCBs might not only be remnants of ancestral gene order but are also maintained by natural selection.

Intraspecific Mitochondrial DNA Comparison of Mycopathogen Mycogone perniciosa Provides Insight Into Mitochondrial Transfer RNA Introns

Mycogone perniciosa is the main causative agent of wet bubble disease, which causes severe damage to the production of the cultivated mushroom Agaricus bisporus around the world. Whole-genome sequencing of 12 isolates of M. perniciosa was performed using the Illumina sequencing platform, and the obtained paired-end reads were used to assemble complete mitochondrial genomes. Intraspecific comparisons of conserved protein-coding genes, transfer RNA (tRNA) and ribosomal RNA (rRNA) genes, introns, and intergenic regions were conducted. Five different mitochondrial DNA (mtDNA) haplotypes were detected among the tested isolates, ranging from 89,080 to 93,199 bp in length. All of the mtDNAs contained the same set of 14 protein-coding genes and 2 rRNA and 27 tRNA genes, which shared high sequence similarity. In contrast, the number, insertion sites, and sequences of introns varied greatly among the mtDNAs. Eighteen of 43 intergenic regions differed among the isolates, reflecting 65 single nucleotide polymorphisms, 76 indels, and the gain/loss of nine long fragments. Intraspecific comparison revealed that two introns were located within tRNA genes, which is the first detailed description of mitochondrial tRNA introns. Intronic sequence comparison within the same insertion sites revealed the formation process of two introns, which also illustrated a fast evolutionary rate of introns among M. perniciosa isolates. Based on the intron distribution pattern, a pair of universal primers and four pairs of isolate-specific primers were designed and were used to identify the five mtDNA types. In summary, the rapid gain or loss of mitochondrial introns could be an ideal marker for population genetics analysis.

Schizosaccharomyces pombe Ppr10 is required for mitochondrial translation

The mitochondrial genome encodes key components of the oxidative phosphorylation (OXPHOS) system, whose expression is essential for mitochondrial functions. We have previously shown that deletion of the Schizosaccharomyces pombe ppr10 encoding a pentatricopeptide repeat protein severely reduces the mature levels of intron-containing mitochondrial transcripts cox1 and cob1, and severely impairs mitochondrial translation. In this study, we examined the possibility that the reduced levels of Cox1 and Cob1 proteins in cells were due to lowered levels of cox1 and cob1 mRNAs. We found that deletion of ppr10 did not affect the levels of mature cox1 and cob1 mRNAs in a mitochondrial intronless background. However, synthesis of Cox1 and Cob1 proteins were still severely affected by deletion of ppr10 in a mitochondrial intronless background. Consistent with this, we found that deletion of mitochondrial introns could not rescue the respiratory growth defect of Δppr10 cells. Our results reveal that Ppr10 is not required for the stability of cox1 and cob1 mRNAs, and provide further support for the idea that Ppr10 plays a critical role in mitochondrial translation.

Large inverted repeats identified by intra-specific comparison of mitochondrial genomes provide insights into the evolution of Agrocybe aegerita

Genomic structure and content of Agrocybe aegerita mitochondrial DNA contain essential information regarding the evolution of this gourmet mushroom. In this study, eight isolates of A. aegerita were sequenced and assembled into complete mitochondrial genomes. The mtDNA of the isolate Ag0067 contained two genotypes, both of which were quadripartite architecture consisting of two identical inverted repeats, separated by a small single-copy region and a large single-copy region. The only difference was opposite directions of the small single-copy region. The mtDNAs ranged from 116,329 bp to 134,035 bp, harboring two large identical inverted repeats. Genes of plasmid-origin were present in regions flanked by inverted repeat ID2. Most of the core genes evolved at a relatively low rate, whereas five tRNA genes located in corresponding regions of Ag0002:1-14000 and Ag0002:50001-61000 showed higher diversity. A long fragment inversion (10 Kb) was suggested to have occurred during the differentiation of two main clades, leading to two different gene orders. The number and distribution of the introns varied greatly among the A. aegerita mtDNAs. Fast invasion of short insertions likely resulted in the diversity of introns as well as other non-coding regions, increasing the variation of the mtDNAs. We raised a model about the evolution of the large repeats to explain the unusual features of A. aegerita mtDNAs. This study constructed quadripartite architecture of A. aegerita mtDNAs analogous to chloroplast DNA, proposed an interconversion model of the divergent mitochondrial genotypes with large inverted repeats. The findings could increase our knowledge of fungal evolution.

Assembly and Analysis of the Complete Mitochondrial Genome of Capsella bursa-pastoris

Shepherd's purse (Capsella bursa-pastoris) is a cosmopolitan annual weed and a promising model plant for studying allopolyploidization in the evolution of angiosperms. Though plant mitochondrial genomes are a valuable source of genetic information, they are hard to assemble. At present, only the complete mitogenome of C. rubella is available out of all species of the genus Capsella. In this work, we have assembled the complete mitogenome of C. bursa-pastoris using high-precision PacBio SMRT third-generation sequencing technology. It is 287,799 bp long and contains 32 protein-coding genes, 3 rRNAs, 25 tRNAs corresponding to 15 amino acids, and 8 open reading frames (ORFs) supported by RNAseq data. Though many repeat regions have been found, none of them is longer than 1 kbp, and the most frequent structural variant originated from these repeats is present in only 4% of the mitogenome copies. The mitochondrial DNA sequence of C. bursa-pastoris differs from C. rubella, but not from C. orientalis, by two long inversions, suggesting that C. orientalis could be its maternal progenitor species. In total, 377 C to U RNA editing sites have been detected. All genes except cox1 and atp8 contain RNA editing sites, and most of them lead to non-synonymous changes of amino acids. Most of the identified RNA editing sites are identical to corresponding RNA editing sites in A. thaliana.