Mitochondrial genome evolution in a single protoploid yeast species

Multiple rearrangements and low inter- and intra-species mitogenome sequence variation in the Heterobasidion annosum s.l. species complex

Introduction

Mitochondria are essential organelles in the eukaryotic cells and responsible for the energy production but are also involved in many other functions including virulence of some fungal species. Although the evolution of fungal mitogenomes have been studied at some taxonomic levels there are still many things to be learned from studies of closely related species.

Methods

In this study, we have analyzed 60 mitogenomes in the five species of the Heterobasidion annosum sensu lato complex that all are necrotrophic pathogens on conifers.

Results and Discussion

Compared to other fungal genera the genomic and genetic variation between and within species in the complex was low except for multiple rearrangements. Several translocations of large blocks with core genes have occurred between the five species and rearrangements were frequent in intergenic areas. Mitogenome lengths ranged between 108 878 to 116 176 bp, mostly as a result of intron variation. There was a high degree of homology of introns, homing endonuclease genes, and intergenic ORFs among the five Heterobasidion species. Three intergenic ORFs with unknown function (uORF6, uORF8 and uORF9) were found in all five species and was located in conserved synteny blocks. A 13 bp long GC-containing self-complementary palindrome was discovered in many places in the five species that were optional in presence/absence. The within species variation is very low, among 48 H. parviporum mitogenomes, there was only one single intron exchange, and SNP frequency was 0.28% and indel frequency 0.043%. The overall low variation in the Heterobasidion annosum sensu lato complex suggests a slow evolution of the mitogenome.

The complete mitochondrial genome of Ophiocordyceps gracilis and its comparison with related species

In this study, the complete mitochondrial genome of O. gracilis was sequenced and assembled before being compared with related species. As the second largest mitogenome reported in the family Ophiocordycipitaceae, the mitogenome of O. gracilis (voucher OG201301) is a circular DNA molecule of 134,288 bp that contains numerous introns and longer intergenomic regions. UCA was detected as anticodon in tRNA-Sec of O. gracilis, while comparative mitogenome analysis of nine Ophiocordycipitaceae fungi indicated that the order and contents of PCGs and rRNA genes were considerably conserved and could descend from a common ancestor in Ophiocordycipitaceae. In addition, the expansion of mitochondrial organization, introns, gene length, and order of O. gracilis were determined to be similar to those of O. sinensis, which indicated common mechanisms underlying adaptive evolution in O. gracilis and O. sinensis. Based on the mitochondrial gene dataset (15 PCGs and 2 RNA genes), a close genetic relationship between O. gracilis and O. sinensis was revealed through phylogenetic analysis. This study is the first to investigate the molecular evolution, phylogenetic pattern, and genetic structure characteristics of mitogenome in O. gracilis. Based on the obtained results, the mitogenome of O. gracilis can increase understanding of the genetic diversity and evolution of cordycipitoid fungi.

Comparison of mitochondrial genomes provides insights into intron dynamics and evolution in Botryosphaeria dothidea and B. kuwatsukai

Botryosphaeria dothidea is one of the most common fungal pathogens on a large number of hosts worldwide. Botryosphaeria dothidea and B. kuwatsukai are also the main causal agents of apple ring rot. In this study, we sequenced, assembled and annotated the circular mitogenomes of 12 diverse B. dothidea isolates (105.7-114.8 kb) infecting various plants including apple, and five diverse B. kuwatsukai isolates (118.0-124.6 kb) from apple. B. dothidea mitogenomes harboured a set of 29-31 introns and 48-52 ORFs. In contrast, B. kuwatsukai mitogenomes harboured more introns (32-34) and ORFs (51-54). The variation in mitogenome sizes was associated mainly with different numbers of introns and insertions of mobile genetic elements. Interestingly, B. dothidea and B. kuwatsukai displayed distinct intron distribution patterns, with three intron loci showing presence/absence dynamics in each species. Large numbers of introns (57% in B. dothidea and 49% in B. kuwatsukai) were most likely obtained through horizontal transfer from non-Dothideomycetes. The mitochondrial gene phylogeny supported the differentiation of the two species. Overall, this study sheds light into the mitochondrial evolution of the plant pathogens B. dothidea and B. kuwatsukai, and intron distribution patterns could be useful markers for studies on population diversity.

Population Analysis and Evolution of Saccharomyces cerevisiae Mitogenomes

The study of mitogenomes allows the unraveling of some paths of yeast evolution that are often not exposed when analyzing the nuclear genome. Although both nuclear and mitochondrial genomes are known to determine phenotypic diversity and fitness, no concordance has yet established between the two, mainly regarding strains' technological uses and/or geographical distribution. In the current work, we proposed a new method to align and analyze yeast mitogenomes, overcoming current difficulties that make it impossible to obtain comparable mitogenomes for a large number of isolates. To this end, 12,016 mitogenomes were considered, and we developed a novel approach consisting of the design of a reference sequence intended to be comparable between all mitogenomes. Subsequently, the population structure of 6646 Saccharomyces cerevisiae mitogenomes was assessed. Results revealed the existence of particular clusters associated with the technological use of the strains, in particular regarding clinical isolates, laboratory strains, and yeasts used for wine-associated activities. As far as we know, this is the first time that a positive concordance between nuclear and mitogenomes has been reported for S. cerevisiae, in terms of strains' technological applications. The results obtained highlighted the importance of including the mtDNA genome in evolutionary analysis, in order to clarify the origin and history of yeast species.

Intra-specific comparison of mitochondrial genomes reveals host gene fragment exchange via intron mobility in Tremella fuciformis

Background

Mitochondrial genomic sequences are known to be variable. Comparative analyses of mitochondrial genomes can reveal the nature and extent of their variation.

Results

Draft mitochondrial genomes of 16 Tremella fuciformis isolates (TF01-TF16) were assembled from Illumina and PacBio sequencing data. Mitochondrial DNA contigs were extracted and assembled into complete circular molecules, ranging from 35,104 bp to 49,044 bp in size. All mtDNAs contained the same set of 41 conserved genes with identical gene order. Comparative analyses revealed that introns and intergenic regions were variable, whereas genic regions (including coding sequences, tRNA, and rRNA genes) were conserved. Among 24 introns detected, 11 were in protein-coding genes, 3 in tRNA genes, and the other 10 in rRNA genes. In addition, two mobile fragments were found in intergenic regions. Interestingly, six introns containing N-terminal duplication of the host genes were found in five conserved protein-coding gene sequences. Comparison of genes with and without these introns gave rise to the following proposed model: gene fragment exchange with other species can occur via gain or loss of introns with N-terminal duplication of the host genes.

Conclusions

Our findings suggest a novel mechanism of fungal mitochondrial gene evolution: partial foreign gene replacement though intron mobility.

Comparative Analyses of Mitochondrial Genomes Provide Evolutionary Insights Into Nematode-Trapping Fungi

Predatory fungi in Orbiliaceae (Ascomycota) have evolved a diversity of trapping devices that enable them to trap and kill nematodes, other small animals, and protozoans. These trapping devices include adhesive hyphae, adhesive knobs, adhesive networks, constricting rings, and non-constricting rings. Their diversity and practical importance have attracted significant attention from biologists, making them excellent model organisms for studying adaptative evolution and as biological control agents against parasitic nematodes. The putative origins and evolutionary relationships among these carnivorous fungi have been investigated using nuclear protein-encoding genes, but their patterns of mitogenome relationships and divergences remain unknown. Here we analyze and compare the mitogenomes of 12 fungal strains belonging to eight species, including six species representing all four types of nematode trapping devices and two from related but non-predatory fungi. All 12 analyzed mitogenomes were of circular DNA molecules, with lengths ranging from 146,101 bp to 280,699 bp. Gene synteny analysis revealed several gene rearrangements and intron transfers among the mitogenomes. In addition, the number of protein coding genes (PCGs), GC content, AT skew, and GC skew varied among these mitogenomes. The increased number and total size of introns were the main contributors to the length differences among the mitogenomes. Phylogenetic analyses of the protein-coding genes indicated that mitochondrial and nuclear genomes evolved at different rates, and signals of positive selection were found in several genes involved in energy metabolism. Our study provides novel insights into the evolution of nematode-trapping fungi and shall facilitate further investigations of this ecologically and agriculturally important group of fungi.

Evidence of Extensive Intraspecific Noncoding Reshuffling in a 169-kb Mitochondrial Genome of a Basidiomycetous Fungus

Comparative genomics of fungal mitochondrial genomes (mitogenomes) have revealed a remarkable pattern of rearrangement between and within major phyla owing to horizontal gene transfer and recombination. The role of recombination was exemplified at a finer evolutionary time scale in basidiomycetes group of fungi as they display a diversity of mitochondrial DNA inheritance patterns. Here, we assembled mitogenomes of six species from the Hymenochaetales order of basidiomycetes and examined 59 mitogenomes from 2 genetic lineages of Phellinus noxius. Gene order is largely collinear, while intergene regions are major determinants of mitogenome size variation. Substantial sequence divergence was found in shared introns consistent with high horizontal gene transfer frequency observed in yeasts, but we also identified a rare case where an intron was retained in five species since speciation. In contrast to the hyperdiversity observed in nuclear genomes of Phellinus noxius, mitogenomes’ intraspecific polymorphisms at protein-coding sequences are extremely low. Phylogeny network based on introns revealed turnover as well as exchange of introns between two lineages. Strikingly, some strains harbor a mosaic origin of introns from both lineages. Analysis of intergenic sequence indicated substantial differences between and within lineages, and an expansion may be ongoing as a result of exchange between distal intergenes. These findings suggest that the evolution in mitochondrial DNAs is usually lineage specific but chimeric mitotypes are frequently observed, thus capturing the possible evolutionary processes shaping mitogenomes in a basidiomycete. The large mitogenome sizes reported in various basidiomycetes appear to be a result of interspecific reshuffling of intergenes.

Intraspecific Diversity of Fission Yeast Mitochondrial Genomes

The fission yeast Schizosaccharomyces pombe is an important model organism, but its natural diversity and evolutionary history remain under-studied. In particular, the population genomics of the S. pombe mitochondrial genome (mitogenome) has not been thoroughly investigated. Here, we assembled the complete circular-mapping mitogenomes of 192 S. pombe isolates de novo, and found that these mitogenomes belong to 69 nonidentical sequence types ranging from 17,618 to 26,910 bp in length. Using the assembled mitogenomes, we identified 20 errors in the reference mitogenome and discovered two previously unknown mitochondrial introns. Analyzing sequence diversity of these 69 types of mitogenomes revealed two highly distinct clades, with only three mitogenomes exhibiting signs of inter-clade recombination. This diversity pattern suggests that currently available S. pombe isolates descend from two long-separated ancestral lineages. This conclusion is corroborated by the diversity pattern of the recombination-repressed K-region located between donor mating-type loci mat2 and mat3 in the nuclear genome. We estimated that the two ancestral S. pombe lineages diverged about 31 million generations ago. These findings shed new light on the evolution of S. pombe and the data sets generated in this study will facilitate future research on genome evolution.

Mitochondrial genome in Hypsizygus marmoreus and its evolution in Dikarya

BACKGROUND

Hypsizygus marmoreus, a high value commercialized edible mushroom is widely cultivated in East Asia, and has become one of the most popular edible mushrooms because of its rich nutritional and medicinal value. Mitochondria are vital organelles, and play various essential roles in eukaryotic cells.

RESULTS

In this study, we provide the Hypsizygus marmoreus mitochondrial (mt) genome assembly: the circular sequence is 102,752 bp in size and contains 15 putative protein-coding genes, 2 ribosomal RNAs subunits and 28 tRNAs. We compared the mt genomes of the 27 fungal species in the Pezizomycotina and Basidiomycotina subphyla, with the results revealing that H. marmoreus is a sister to Tricholoma matsutake and the phylogenetic distribution of this fungus based on the mt genome. Phylogenetic analysis shows that Ascomycetes mitochondria started to diverge earlier than that of Basidiomycetes and supported the robustness of the hyper metric tree. The fungal sequences are highly polymorphic and gene order varies significantly in the dikarya data set, suggesting a correlation between the gene order and divergence time in the fungi mt genome. To detect the mt genome variations in H. marmoreus, we analyzed the mtDNA sequences of 48 strains. The phylogeny and variation sited type statistics of H. marmoreus provide clear-cut evidence for the existence of four well-defined cultivations isolated lineages, suggesting female ancestor origin of H. marmoreus. Furthermore, variations on two loci were further identified to be molecular markers for distinguishing the subgroup containing 32 strains of other strains. Fifteen conserved protein-coding genes of mtDNAs were analyzed, with fourteen revealed to be under purifying selection in the examined fungal species, suggesting the rapid evolution was caused by positive selection of this gene.

CONCLUSIONS

Our studies have provided new reference mt genomes and comparisons between species and intraspecies with other strains, and provided future perspectives for assessing diversity and origin of H. marmoreus.

Genetic Drift and Indel Mutation in the Evolution of Yeast Mitochondrial Genome Size

Mitochondrial genomes (mitogenomes) are remarkably diverse in genome size and organization, but the origins of dynamic mitogenome architectures are still poorly understood. For instance, the mutational burden hypothesis postulates that the drastic difference between large plant mitogenomes and streamlined animal mitogenomes can be driven by their different mutation rates. However, inconsistent trends between mitogenome sizes and mutation rates have been documented in several lineages. These conflicting results highlight the need of systematic and sophisticated investigations on the evolution and diversity of mitogenome architecture. This study took advantage of the strikingly variable mitogenome size among different yeast species and also among intraspecific strains, examined sequence dynamics of introns, GC-clusters, tandem repeats, mononucleotide repeats (homopolymers) and evaluated their contributions to genome size variation. The contributions of these sequence features to mitogenomic variation are dependent on the timescale, over which extant genomes evolved from their last common ancestor, perhaps due to a combination of different turnover rates of mobile sequences, variable insertion spaces, and functional constraints. We observed a positive correlation between mitogenome size and the level of genetic drift, suggesting that mitogenome expansion in yeast is likely driven by multiple types of sequence insertions in a primarily nonadaptive manner. Although these cannot be explained directly by the mutational burden hypothesis, our results support an important role of genetic drift in the evolution of yeast mitogenomes.

Fungal mitochondrial genomes and genetic polymorphisms

Mitochondria are the powerhouses of eukaryotic cells, responsible for ATP generation and playing a role in a diversity of cellular and organismal functions. Different from the majority of other intracellular membrane structures, mitochondria contain their own genetic materials that are capable of independent replication and inheritance. In this mini-review, we provide brief summaries of fungal mitochondrial genome structure, size, gene content, inheritance, and genetic variation. We pay special attention to the relative genetic polymorphisms of the mitochondrial vs nuclear genomes at the population level within individual fungal species. Among the 20 species/groups of species reviewed here, there is a range of variation among genes and species in the relative nuclear and mitochondrial genetic polymorphisms. Interestingly, most (15/20) showed a greater genetic diversity for nuclear genes and genomes than for mitochondrial genes and genomes, with the remaining five showing similar or slower nuclear genome genetic variations. This fungal pattern is different from the dominant pattern in animals, but more similar to that in plants. At present, the mechanisms for the variations among fungal species and the overall low level of mitochondrial sequence polymorphisms are not known. The increasing availability of population genomic data should help us reveal the potential genetic and ecological factors responsible for the observed variations.

Mt-rps3 is an ancient gene which provides insight into the evolution of fungal mitochondrial genomes

The nuclear ribosomal protein S3 (Rps3) is implicated in the assembly of the ribosomal small subunit. Fungi and plants present a gene copy in their mitochondrial (mt) genomes. An analysis of 303 complete fungal mt genomes showed that, when rps3 is found, it is either a free-standing gene or an anchored gene within the omega intron of the rnl gene. Early divergent fungi, Basidiomycota and all yeasts but the CTG group belong to the first case, and Pezizomycotina to the second. Its position, size and genetic code employed are conserved within species of the same Order. Size variability is attributed to different number of repeats. These repeats consist of AT-rich sequences. MtRps3 proteins lack the KH domain, necessary for binding to rRNA, in their N-terminal region. Their C-terminal region is conserved in all Domains of life. Phylogenetic analysis showed that nuclear and mtRps3 proteins are descendants of archaeal and a-proteobacterial homologues, respectively. Thus, fungal mt-rps3 gene is an ancient gene which evolved within the endosymbiotic model and presents different evolutionary routes: (a) coming from a-proteobacteria, it was relocated to another region of the mt genome, (b) via its insertion to the omega intron, it was transferred to the nucleus and/or got lost, and (c) it was re-routed to the mt genome again. Today, Basidiomycota and Saccharomycetales seem to follow the first evolutionary route and almost all Pezizomycotina support the second scenario with their exceptions being the result of the third scenario, i.e., the gene's re-entry to the mt genome.

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.

Mobile Introns Shape the Genetic Diversity of Their Host Genes

Self-splicing introns populate several highly conserved protein-coding genes in fungal and plant mitochondria. In fungi, many of these introns have retained their ability to spread to intron-free target sites, often assisted by intron-encoded endonucleases that initiate the homing process. Here, leveraging population genomic data from Saccharomyces cerevisiae, Schizosaccharomyces pombe, and Lachancea kluyveri, we expose nonrandom patterns of genetic diversity in exons that border self-splicing introns. In particular, we show that, in all three species, the density of single nucleotide polymorphisms increases as one approaches a mobile intron. Through multiple lines of evidence, we rule out relaxed purifying selection as the cause of uneven nucleotide diversity. Instead, our findings implicate intron mobility as a direct driver of host gene diversity. We discuss two mechanistic scenarios that are consistent with the data: either endonuclease activity and subsequent error-prone repair have left a mutational footprint on the insertion environment of mobile introns or nonrandom patterns of genetic diversity are caused by exonic coconversion, which occurs when introns spread to empty target sites via homologous recombination. Importantly, however, we show that exonic coconversion can only explain diversity gradients near intron-exon boundaries if the conversion template comes from outside the population. In other words, there must be pervasive and ongoing horizontal gene transfer of self-splicing introns into extant fungal populations.

Comparative mitochondrial genomics toward exploring molecular markers in the medicinal fungus Cordyceps militaris

Cordyceps militaris is a fungus used for developing health food, but knowledge about its intraspecific differentiation is limited due to lack of efficient markers. Herein, we assembled the mitochondrial genomes of eight C. militaris strains and performed a comparative mitochondrial genomic analysis together with three previously reported mitochondrial genomes of the fungus. Sizes of the 11 mitochondrial genomes varied from 26.5 to 33.9 kb mainly due to variable intron contents (from two to eight introns per strain). Nucleotide variability varied according to different regions with non-coding regions showing higher variation frequency than coding regions. Recombination events were identified between some locus pairs but seemed not to contribute greatly to genetic variations of the fungus. Based on nucleotide diversity fluctuations across the alignment of all mitochondrial genomes, molecular markers with the potential to be used for future typing studies were determined.

Competition assays and physiological experiments of soil and phyllosphere yeasts identify Candida subhashii as a novel antagonist of filamentous fungi

BACKGROUND

While recent advances in next generation sequencing technologies have enabled researchers to readily identify countless microbial species in soil, rhizosphere, and phyllosphere microbiomes, the biological functions of the majority of these species are unknown. Functional studies are therefore urgently needed in order to characterize the plethora of microorganisms that are being identified and to point out species that may be used for biotechnology or plant protection. Here, we used a dual culture assay and growth analyses to characterise yeasts (40 different isolates) and their antagonistic effect on 16 filamentous fungi; comprising plant pathogens, antagonists, and saprophytes.

RESULTS

Overall, this competition screen of 640 pairwise combinations revealed a broad range of outcomes, ranging from small stimulatory effects of some yeasts up to a growth inhibition of more than 80% by individual species. On average, yeasts isolated from soil suppressed filamentous fungi more strongly than phyllosphere yeasts and the antagonistic activity was a species-/isolate-specific property and not dependent on the filamentous fungus a yeast was interacting with. The isolates with the strongest antagonistic activity were Metschnikowia pulcherrima, Hanseniaspora sp., Cyberlindnera sargentensis, Aureobasidium pullulans, Candida subhashii, and Pichia kluyveri. Among these, the soil yeasts (C. sargentensis, A. pullulans, C. subhashii) assimilated and/or oxidized more di-, tri- and tetrasaccharides and organic acids than yeasts from the phyllosphere. Only the two yeasts C. subhashii and M. pulcherrima were able to grow with N-acetyl-glucosamine as carbon source.

CONCLUSIONS

The competition assays and physiological experiments described here identified known antagonists that have been implicated in the biological control of plant pathogenic fungi in the past, but also little characterised species such as C. subhashii. Overall, soil yeasts were more antagonistic and metabolically versatile than yeasts from the phyllosphere. Noteworthy was the strong antagonistic activity of the soil yeast C. subhashii, which had so far only been described from a clinical sample and not been studied with respect to biocontrol. Based on binary competition assays and growth analyses (e.g., on different carbon sources, growth in root exudates), C. subhashii was identified as a competitive and antagonistic soil yeast with potential as a novel biocontrol agent against plant pathogenic fungi.

Dissection of quantitative traits by bulk segregant mapping in a protoploid yeast species

Since more than a decade ago, Saccharomyces cerevisiae has been used as a model to dissect complex traits, revealing the genetic basis of a large number of traits in fine detail. However, to have a more global view of the genetic architecture of traits across species, the examination of the molecular basis of phenotypes within non-conventional species would undoubtedly be valuable. In this respect, the Saccharomycotina yeasts represent ideal and potential non-model organisms. Here we sought to assess the feasibility of genetic mapping by bulk segregant analysis in the protoploid Lachancea kluyveri (formerly S. kluyveri) yeast species, a distantly related species to S. cerevisiae For this purpose, we designed a fluorescent mating-type marker, compatible with any mating-competent strains representative of this species, to rapidly create a large population of haploid segregants (>10(5) cells). Quantitative trait loci can be mapped by selecting and sequencing an enriched pool of progeny with extreme phenotypic values. As a test bed, we applied this strategy and mapped the causal loci underlying halotolerance phenotypes in L. kluyveri Overall, this study demonstrates that bulk segregant mapping is a powerful way for investigating the genetic basis of natural variations in non-model yeast organisms and more precisely in L. kluyveri.

Analysis of the Mitochondrial Genome in Hypomyces aurantius Reveals a Novel Twintron Complex in Fungi

Hypomyces aurantius is a mycoparasite that causes cobweb disease, a most serious disease of cultivated mushrooms. Intra-species identification is vital for disease control, however the lack of genomic data makes development of molecular markers challenging. Small size, high copy number, and high mutation rate of fungal mitochondrial genome makes it a good candidate for intra and inter species differentiation. In this study, the mitochondrial genome of H. H.a0001 was determined from genomic DNA using Illumina sequencing. The roughly 72 kb genome shows all major features found in other Hypocreales: 14 common protein genes, large and small subunit rRNAs genes and 27 tRNAs genes. Gene arrangement comparison showed conserved gene orders in Hypocreales mitochondria are relatively conserved, with the exception of Acremonium chrysogenum and Acremonium implicatum. Mitochondrial genome comparison also revealed that intron length primarily contributes to mitogenome size variation. Seventeen introns were detected in six conserved genes: five in cox1, four in rnl, three in cob, two each in atp6 and cox3, and one in cox2. Four introns were found to contain two introns or open reading frames: cox3-i2 is a twintron containing two group IA type introns; cox2-i1 is a group IB intron encoding two homing endonucleases; and cox1-i4 and cox1-i3 both contain two open reading frame (ORFs). Analyses combining secondary intronic structures, insertion sites, and similarities of homing endonuclease genes reveal two group IA introns arranged side by side within cox3-i2. Mitochondrial data for H. aurantius provides the basis for further studies relating to population genetics and species identification.

Lachancea quebecensis sp. nov., a yeast species consistently isolated from tree bark in the Canadian province of Québec

A thorough sampling of maple, oak, birch, and apple tree bark in North America yielded a set of isolates that represent a yeast species not yet formally described. The strains obtained were all isolated from the Canadian province of Québec. These four isolates have identical electrophoretic karyotypes, distinct from other species of the genus Lachancea, and are most closely related to the formally recognized species Lachancea thermotolerans according to the D1/D2 domain of the LSU rDNA gene and 5.8S–ITS region. Previous studies revealed the existence of a population of strains closely related to L. thermotolerans, with unique D1/D2 sequences and the ability to grow on melibiose, which is also true for these isolates. The sequences obtained here (for the D1/D2, and 5.8S–ITS region) are identical among the four strains, and in a phylogenetic analysis of the D1/D2 region, the strains form a distinct clade with the previously described population closely related to L. thermotolerans, composed of isolates from Japan, as well as from the provinces of Ontario and Québec in Canada. On the basis of select physiological and phylogenetic characteristics, a novel ascosporogenous yeast species, Lachancea quebecensis sp. nov., is proposed. The type strain LL11_022T ( = CBS 14138T = CLIB 1763T = UCDFST 15-106T) was isolated from maple tree bark in the Station Duchesnay, QC region of Québec, Canada. The MycoBank number is MB811749.

Large-Scale Survey of Intraspecific Fitness and Cell Morphology Variation in a Protoploid Yeast Species

It is now clear that the exploration of the genetic and phenotypic diversity of nonmodel species greatly improves our knowledge in biology. In this context, we recently launched a population genomic analysis of the protoploid yeast Lachancea kluyveri (formerly Saccharomyces kluyveri), highlighting a broad genetic diversity (π = 17 × 10⁻³) compared to the yeast model organism, S. cerevisiae (π = 4 × 10⁻³). Here, we sought to generate a comprehensive view of the phenotypic diversity in this species. In total, 27 natural L. kluyveri isolates were subjected to trait profiling using the following independent approaches: (i) analyzing growth in 55 growth conditions and (ii) investigating 501 morphological changes at the cellular level. Despite higher genetic diversity, the fitness variance observed in L. kluyveri is lower than that in S. cerevisiae. However, morphological features show an opposite trend. In addition, there is no correlation between the origins (ecological or geographical) of the isolate and the phenotypic patterns, demonstrating that trait variation follows neither population history nor source environment in L. kluyveri. Finally, pairwise comparisons between growth rate correlation and genetic diversity show a clear decrease in phenotypic variability linked to genome variation increase, whereas no such a trend was identified for morphological changes. Overall, this study reveals for the first time the phenotypic diversity of a distantly related species to S. cerevisiae. Given its genetic properties, L. kluyveri might be useful in further linkage mapping analyses of complex traits, and could ultimately provide a better insight into the evolution of the genotype–phenotype relationship across yeast species.

Differences in environmental stress response among yeasts is consistent with species-specific lifestyles

Defining how organisms respond to environmental change has always been an important step toward understanding their adaptive capacity and physiology. Variation in transcription during stress has been widely described in model species, especially in the yeast Saccharomyces cerevisiae, which helped to shape general rules regarding how cells cope with environmental constraints, as well as to decipher the functions of many genes. Comparison of the environmental stress response (ESR) across species is essential to obtaining better insight into the common and species-specific features of stress defense. In this context, we explored the transcriptional landscape of the yeast Lachancea kluyveri (formerly Saccharomyces kluyveri) in response to diverse stresses, using RNA sequencing. We investigated variation in gene expression and observed a link between genetic plasticity and environmental sensitivity. We identified the ESR genes in this species and compared them to those already found in S. cerevisiae We observed common features between the two species, as well as divergence in the regulatory networks involved. Of interest, some changes were related to differences in species lifestyle. Thus we were able to decipher how adaptation to stress has evolved among different yeast species. Finally, by analyzing patterns of coexpression, we were able to propose potential biological functions for 42% of genes and also annotate 301 genes for which no function could be assigned by homology. This large data set allowed for the characterization of the evolution of gene regulation and provides an efficient tool for assessing gene function.

Complete mitochondrial genome of the Verticillium-wilt causing plant pathogen Verticillium nonalfalfae

Verticillium nonalfalfae is a fungal plant pathogen that causes wilt disease by colonizing the vascular tissues of host plants. The disease induced by hop isolates of V. nonalfalfae manifests in two different forms, ranging from mild symptoms to complete plant dieback, caused by mild and lethal pathotypes, respectively. Pathogenicity variations between the causal strains have been attributed to differences in genomic sequences and perhaps also to differences in their mitochondrial genomes. We used data from our recent Illumina NGS-based project of genome sequencing V. nonalfalfae to study the mitochondrial genomes of its different strains. The aim of the research was to prepare a V. nonalfalfae reference mitochondrial genome and to determine its phylogenetic placement in the fungal kingdom. The resulting 26,139 bp circular DNA molecule contains a full complement of the 14 "standard" fungal mitochondrial protein-coding genes of the electron transport chain and ATP synthase subunits, together with a small rRNA subunit, a large rRNA subunit, which contains ribosomal protein S3 encoded within a type IA-intron and 26 tRNAs. Phylogenetic analysis of this mitochondrial genome placed it in the Verticillium spp. lineage in the Glomerellales group, which is also supported by previous phylogenetic studies based on nuclear markers. The clustering with the closely related Verticillium dahliae mitochondrial genome showed a very conserved synteny and a high sequence similarity. Two distinguishing mitochondrial genome features were also found-a potential long non-coding RNA (orf414) contained only in the Verticillium spp. of the fungal kingdom, and a specific fragment length polymorphism observed only in V. dahliae and V. nubilum of all the Verticillium spp., thus showing potential as a species specific biomarker.

Population genomics of yeasts: towards a comprehensive view across a broad evolutionary scale

With the advent of high-throughput technologies for sequencing, the complete description of the genetic variation that occurs in populations, also known as population genomics, is foreseeable but far from being reached. Explaining the forces that govern patterns of genetic variation is essential to elucidate the evolutionary history of species. Genetic variation results from a wide assortment of evolutionary forces, among which mutation, selection, recombination and drift play major roles in shaping genomes. In addition, exploring the genetic variation within a population also corresponds to the first step towards dissecting the genotype-phenotype relationship. In this context, yeast species are of particular interest because they represent a unique resource for studying the evolution of intraspecific genetic diversity in a phylum spanning a broad evolutionary scale. Here, we briefly review recent progress in yeast population genomics and provide some perspective on this rapidly evolving field. In fact, we truly believe that it is of interest to supplement comparative and early population genomic studies with the deep sequencing of more extensive sets of individuals from the same species. In parallel, it would be more than valuable to uncover the intraspecific variation of a large number of unexplored species, including those that are closely and more distantly related. Altogether, these data would enable substantially more powerful genomic scans for functional dissection.

The Genetics of Non-conventional Wine Yeasts: Current Knowledge and Future Challenges

Saccharomyces cerevisiae is by far the most widely used yeast in oenology. However, during the last decade, several other yeasts species has been purposed for winemaking as they could positively impact wine quality. Some of these non-conventional yeasts (Torulaspora delbrueckii, Metschnikowia pulcherrima, Pichia kluyveri, Lachancea thermotolerans, etc.) are now proposed as starters culture for winemakers in mixed fermentation with S. cerevisiae, and several others are the subject of various studies (Hanseniaspora uvarum, Starmerella bacillaris, etc.). Along with their biotechnological use, the knowledge of these non-conventional yeasts greatly increased these last 10 years. The aim of this review is to describe the last updates and the current state-of-art of the genetics of non-conventional yeasts (including S. uvarum, T. delbrueckii, S. bacillaris, etc.). We describe how genomics and genetics tools provide new data into the population structure and biodiversity of non-conventional yeasts in winemaking environments. Future challenges will lie on the development of selection programs and/or genetic improvement of these non-conventional species. We discuss how genetics, genomics and the advances in next-generation sequencing will help the wine industry to develop the biotechnological use of non-conventional yeasts to improve the quality and differentiation of wines.

Whole-Genome Sequencing and Intraspecific Analysis of the Yeast Species Lachancea quebecensis

The gold standard in yeast population genomics has been the model organism Saccharomyces cerevisiae. However, the exploration of yeast species outside the Saccharomyces genus is essential to broaden the understanding of genome evolution. Here, we report the analyses of whole-genome sequences of nineisolates from the recently described yeast species Lachancea quebecensis. The genome of one isolate was assembled and annotated, and the intraspecific variability within L. quebecensis was surveyed by comparing the sequences from the eight other isolates to this reference sequence. Our study revealed that these strains harbor genomes with an average nucleotide diversity of π = 2 × 10⁻³ which is slightly lower, although on the same order of magnitude, as that previously determined for S. cerevisiae (π = 4 × 10⁻³). Our results show that even though these isolates were all obtained from a relatively isolated geographic location, the same ecological source, and represent a smaller sample size than is available for S. cerevisiae, the levels of divergence are similar to those observed in this model species. This divergence is essentially linked to the presence of two distinct clusters delineated according to geographic location. However, even with relatively similar ranges of genome divergence, L. quebecensis has an extremely low global phenotypic variance of 0.062 compared with 0.59 previously determined in S. cerevisiae.

The complete nucleotide sequence of the mitochondrial genome of Drosophila formosana (Diptera: Drosophilidae)

Drosophila formosana (Diptera: Drosophilidae) belongs to the Drosophilidae group of Drosophila. The mitochondrial genome sequence of Drosophila formosana is determined in this study. Mitochondrion of D. formosana is a circular DNA molecule of the 16 100 nucleotide pairs (bp) that contains one encoding region including 37 genes and 1 non-coding A + T-rich region. The similarity and typicality have been showed by the structure and organization analysis. All genes are arranged in the circular DNA molecule. In addition to DN5 that use GTG start codon, all other protein-coding genes (PCGs) start with an ATN start codon. Ten protein-coding genes stop with the termination codon TAN, while other protein-coding genes (PCGs) used incomplete termination codon TA- (cox2, nad5, nad1). The A + T-rich region with a length of 1088 bp is located between rrnS and trnI. The mitochondrial genome of D. formosana has been completely sequenced for the first time in this study.

The evolution of sex: A new hypothesis based on mitochondrial mutational erosion: Mitochondrial mutational erosion in ancestral eukaryotes would favor the evolution of sex, harnessing nuclear recombination to optimize compensatory nuclear coadaptation

The evolution of sex in eukaryotes represents a paradox, given the "twofold" fitness cost it incurs. We hypothesize that the mutational dynamics of the mitochondrial genome would have favored the evolution of sexual reproduction. Mitochondrial DNA (mtDNA) exhibits a high-mutation rate across most eukaryote taxa, and several lines of evidence suggest that this high rate is an ancestral character. This seems inexplicable given that mtDNA-encoded genes underlie the expression of life's most salient functions, including energy conversion. We propose that negative metabolic effects linked to mitochondrial mutation accumulation would have invoked selection for sexual recombination between divergent host nuclear genomes in early eukaryote lineages. This would provide a mechanism by which recombinant host genotypes could be rapidly shuffled and screened for the presence of compensatory modifiers that offset mtDNA-induced harm. Under this hypothesis, recombination provides the genetic variation necessary for compensatory nuclear coadaptation to keep pace with mitochondrial mutation accumulation.

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.

Extensive Horizontal Transfer and Homologous Recombination Generate Highly Chimeric Mitochondrial Genomes in Yeast

The frequency of horizontal gene transfer (HGT) in mitochondrial DNA varies substantially. In plants, HGT is relatively common, whereas in animals it appears to be quite rare. It is of considerable importance to understand mitochondrial HGT across the major groups of eukaryotes at a genome-wide level, but so far this has been well studied only in plants. In this study, we generated ten new mitochondrial genome sequences and analyzed 40 mitochondrial genomes from the Saccharomycetaceae to assess the magnitude and nature of mitochondrial HGT in yeasts. We provide evidence for extensive, homologous-recombination-mediated, mitochondrial-to-mitochondrial HGT occurring throughout yeast mitochondrial genomes, leading to genomes that are highly chimeric evolutionarily. This HGT has led to substantial intraspecific polymorphism in both sequence content and sequence divergence, which to our knowledge has not been previously documented in any mitochondrial genome. The unexpectedly high frequency of mitochondrial HGT in yeast may be driven by frequent mitochondrial fusion, relatively low mitochondrial substitution rates and pseudohyphal fusion to produce heterokaryons. These findings suggest that mitochondrial HGT may play an important role in genome evolution of a much broader spectrum of eukaryotes than previously appreciated and that there is a critical need to systematically study the frequency, extent, and importance of mitochondrial HGT across eukaryotes.

Mitochondrial genome evolution in yeasts: an all-encompassing view

Mitochondria are important organelles that harbor their own genomes encoding a key set of proteins that ensure respiration and provide the eukaryotic cell with energy. Recent advances in high-throughput sequencing technologies present a unique opportunity to explore mitochondrial (mt) genome evolution. The Saccharomycotina yeasts have proven to be the leading organisms for mt comparative and population genomics. In fact, the explosion of complete yeast mt genome sequences has allowed for a broader view of the mt diversity across this incredibly diverse subphylum, both within and between closely related species. Here, we summarize the present state of yeast mitogenomics, including the currently available data and what it reveals concerning the diversity of content, organization, structure and evolution of mt genomes.

Comparison of mitochondrial genomes provides insights into intron dynamics and evolution in the caterpillar fungus Cordyceps militaris

Intra-specific comparison of mitochondrial genomes can help elucidate the evolution of a species, however it has not been performed for hypocrealean fungi that form diverse symbiotic associations with other organisms. In this study, comparative analyses of three completely sequenced mitochondrial genomes of a hypocrealean fungus, Cordyceps militaris, the type species of Cordyceps genus, revealed that the introns were the main contributors to mitochondrial genome size variations among strains. Mitochondrial genes in C. militaris have been invaded by group I introns in at least eight positions. PCR assays of various C. militaris isolates showed abundant variations of intron presence/absence among strains at seven of the eight intronic loci. Although the ancestral intron pattern was inferred to contain all eight introns, loss and/or gain events occurred for seven of the eight introns. These introns invaded the C. militaris mitochondrial genome probably by horizontal transfer from other fungi, and intron insertions into intronless genes in C. militaris were accompanied by co-conversions of upstream exon sequences especially for those introns targeting protein-coding genes. We also detected phylogenetic congruence between the intron and exon trees at each individual locus, consistent with the ancestral mitochondria of C. militaris as having all eight introns. This study helps to explain the evolution of C. militaris mitochondrial genomes and will facilitate population genetic studies of this medicinally important fungus.

Protocols and programs for high-throughput growth and aging phenotyping in yeast

In microorganisms, and more particularly in yeasts, a standard phenotyping approach consists in the analysis of fitness by growth rate determination in different conditions. One growth assay that combines high throughput with high resolution involves the generation of growth curves from 96-well plate microcultivations in thermostated and shaking plate readers. To push the throughput of this method to the next level, we have adapted it in this study to the use of 384-well plates. The values of the extracted growth parameters (lag time, doubling time and yield of biomass) correlated well between experiments carried out in 384-well plates as compared to 96-well plates or batch cultures, validating the higher-throughput approach for phenotypic screens. The method is not restricted to the use of the budding yeast Saccharomyces cerevisiae, as shown by consistent results for other species selected from the Hemiascomycete class. Furthermore, we used the 384-well plate microcultivations to develop and validate a higher-throughput assay for yeast Chronological Life Span (CLS), a parameter that is still commonly determined by a cumbersome method based on counting "Colony Forming Units". To accelerate analysis of the large datasets generated by the described growth and aging assays, we developed the freely available software tools GATHODE and CATHODE. These tools allow for semi-automatic determination of growth parameters and CLS behavior from typical plate reader output files. The described protocols and programs will increase the time- and cost-efficiency of a number of yeast-based systems genetics experiments as well as various types of screens.

Population genomics reveals chromosome-scale heterogeneous evolution in a protoploid yeast

Yeast species represent an ideal model system for population genomic studies but large-scale polymorphism surveys have only been reported for species of the Saccharomyces genus so far. Hence, little is known about intraspecific diversity and evolution in yeast. To obtain a new insight into the evolutionary forces shaping natural populations, we sequenced the genomes of an expansive worldwide collection of isolates from a species distantly related to Saccharomyces cerevisiae: Lachancea kluyveri (formerly S. kluyveri). We identified 6.5 million single nucleotide polymorphisms and showed that a large introgression event of 1 Mb of GC-rich sequence in the chromosomal arm probably occurred in the last common ancestor of all L. kluyveri strains. Our population genomic data clearly revealed that this 1-Mb region underwent a molecular evolution pattern very different from the rest of the genome. It is characterized by a higher recombination rate, with a dramatically elevated A:T → G:C substitution rate, which is the signature of an increased GC-biased gene conversion. In addition, the predicted base composition at equilibrium demonstrates that the chromosome-scale compositional heterogeneity will persist after the genome has reached mutational equilibrium. Altogether, the data presented herein clearly show that distinct recombination and substitution regimes can coexist and lead to different evolutionary patterns within a single genome.

Population genomic analysis reveals highly conserved mitochondrial genomes in the yeast species Lachancea thermotolerans

The increasing availability of mitochondrial (mt) sequence data from various yeasts provides a tool to study genomic evolution within and between different species. While the genomes from a range of lineages are available, there is a lack of information concerning intraspecific mtDNA diversity. Here, we analyzed the mt genomes of 50 strains from Lachancea thermotolerans, a protoploid yeast species that has been isolated from several locations (Europe, Asia, Australia, South Africa, and North / South America) and ecological sources (fruit, tree exudate, plant material, and grape and agave fermentations). Protein-coding genes from the mtDNA were used to construct a phylogeny, which reflected a similar, yet less resolved topology than the phylogenetic tree of 50 nuclear genes. In comparison to its sister species Lachancea kluyveri, L. thermotolerans has a smaller mt genome. This is due to shorter intergenic regions and fewer introns, of which the latter are only found in COX1. We revealed that L. kluyveri and L. thermotolerans share similar levels of intraspecific divergence concerning the nuclear genomes. However, L. thermotolerans has a more highly conserved mt genome with the coding regions characterized by low rates of nonsynonymous substitution. Thus, in the mt genomes of L. thermotolerans, stronger purifying selection and lower mutation rates potentially shape genome diversity in contract to what was found for L. kluyveri, demonstrating that the factors driving mt genome evolution are different even between closely related species.

A genome-wide map of mitochondrial DNA recombination in yeast

In eukaryotic cells, the production of cellular energy requires close interplay between nuclear and mitochondrial genomes. The mitochondrial genome is essential in that it encodes several genes involved in oxidative phosphorylation. Each cell contains several mitochondrial genome copies and mitochondrial DNA recombination is a widespread process occurring in plants, fungi, protists, and invertebrates. Saccharomyces cerevisiae has proved to be an excellent model to dissect mitochondrial biology. Several studies have focused on DNA recombination in this organelle, yet mostly relied on reporter genes or artificial systems. However, no complete mitochondrial recombination map has been released for any eukaryote so far. In the present work, we sequenced pools of diploids originating from a cross between two different S. cerevisiae strains to detect recombination events. This strategy allowed us to generate the first genome-wide map of recombination for yeast mitochondrial DNA. We demonstrated that recombination events are enriched in specific hotspots preferentially localized in non-protein-coding regions. Additionally, comparison of the recombination profiles of two different crosses showed that the genetic background affects hotspot localization and recombination rates. Finally, to gain insights into the mechanisms involved in mitochondrial recombination, we assessed the impact of individual depletion of four genes previously associated with this process. Deletion of NTG1 and MGT1 did not substantially influence the recombination landscape, alluding to the potential presence of additional regulatory factors. Our findings also revealed the loss of large mitochondrial DNA regions in the absence of MHR1, suggesting a pivotal role for Mhr1 in mitochondrial genome maintenance during mating. This study provides a comprehensive overview of mitochondrial DNA recombination in yeast and thus paves the way for future mechanistic studies of mitochondrial recombination and genome maintenance.

Comparative mitochondrial genomics within and among yeast species of the Lachancea genus

Yeasts are leading model organisms for mitochondrial genome studies. The explosion of complete sequence of yeast mitochondrial (mt) genomes revealed a wide diversity of organization and structure between species. Recently, genome-wide polymorphism survey on the mt genome of isolates of a single species, Lachancea kluyveri, was also performed. To compare the mitochondrial genome evolution at two hierarchical levels: within and among closely related species, we focused on five species of the Lachancea genus, which are close relatives of L. kluyveri. Hence, we sequenced the complete mt genome of L. dasiensis, L. nothofagi, L. mirantina, L. fantastica and L. meyersii. The phylogeny of the Lachancea genus was explored using these data. Analysis of intra- and interspecific variability across the whole Lachancea genus led to the same conclusions regarding the mitochondrial genome evolution. These genomes exhibit a similar architecture and are completely syntenic. Nevertheless, genome sizes vary considerably because of the variations of the intergenic regions and the intron content, contributing to mitochondrial genome plasticity. The high variability of the intergenic regions stands in contrast to the high level of similarity of protein sequences. Quantification of the selective constraints clearly revealed that most of the mitochondrial genes are under purifying selection in the whole genus.