Evolution of mitochondrial genomes in fungi

The Mitochondrial Genome of a Plant Fungal Pathogen Pseudocercospora fijiensis (Mycosphaerellaceae), Comparative Analysis and Diversification Times of the Sigatoka Disease Complex Using Fossil Calibrated Phylogenies

Mycosphaerellaceae is a highly diverse fungal family containing a variety of pathogens affecting many economically important crops. Mitochondria play a crucial role in fungal metabolism and in the study of fungal evolution. This study aims to: (i) describe the mitochondrial genome of Pseudocercospora fijiensis, and (ii) compare it with closely related species (Sphaerulina musiva, S. populicola, P. musae and P. eumusae) available online, paying particular attention to the Sigatoka disease’s complex causal agents. The mitochondrial genome of P. fijiensis is a circular molecule of 74,089 bp containing typical genes coding for the 14 proteins related to oxidative phosphorylation, 2 rRNA genes and a set of 38 tRNAs. P. fijiensis mitogenome has two truncated cox1 copies, and bicistronic transcription of nad2-nad3 and atp6-atp8 confirmed experimentally. Comparative analysis revealed high variability in size and gene order among selected Mycosphaerellaceae mitogenomes likely to be due to rearrangements caused by mobile intron invasion. Using fossil calibrated Bayesian phylogenies, we found later diversification times for Mycosphaerellaceae (66.6 MYA) and the Sigatoka disease complex causal agents, compared to previous strict molecular clock studies. An early divergent Pseudocercospora fijiensis split from the sister species P. musae + P. eumusae 13.31 MYA while their sister group, the sister species P. eumusae and P. musae, split from their shared common ancestor in the late Miocene 8.22 MYA. This newly dated phylogeny suggests that species belonging to the Sigatoka disease complex originated after wild relatives of domesticated bananas (section Eumusae; 27.9 MYA). During this time frame, mitochondrial genomes expanded significantly, possibly due to invasions of introns into different electron transport chain genes.

Mitochondrial evolution in the entomopathogenic fungal genus Beauveria

Species in the fungal genus Beauveria are pathogens of invertebrates and have been commonly used as the active agent in biopesticides. After many decades with few species described, recent molecular approaches to classification have led to over 25 species now delimited. Little attention has been given to the mitochondrial genomes of Beauveria but better understanding may led to insights into the nature of species and evolution in this important genus. In this study, we sequenced the mitochondrial genomes of four new strains belonging to Beauveria bassiana, Beauveria caledonica and Beauveria malawiensis, and compared them to existing mitochondrial sequences of related fungi. The mitochondrial genomes of Beauveria ranged widely from 28,806 to 44,135 base pairs, with intron insertions accounting for most size variation and up to 39% (B. malawiensis) of the mitochondrial length due to introns in genes. Gene order of the common mitochondrial genes did not vary among the Beauveria sequences, but variation was observed in the number of transfer ribonucleic acid genes. Although phylogenetic analysis using whole mitochondrial genomes showed, unsurprisingly, that B. bassiana isolates were the most closely related to each other, mitochondrial codon usage suggested that some B. bassiana isolates were more similar to B. malawiensis and B. caledonica than the other B. bassiana isolates analyzed.

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.

Unveiling the mystery of mitochondrial DNA replication in yeasts

Conventional DNA replication is initiated from specific origins and requires the synthesis of RNA primers for both the leading and lagging strands. In contrast, the replication of yeast mitochondrial DNA is origin-independent. The replication of the leading strand is likely primed by recombinational structures and proceeded by a rolling circle mechanism. The coexistent linear and circular DNA conformers facilitate the recombination-based initiation. The replication of the lagging strand is poorly understood. Re-evaluation of published data suggests that the rolling circle may also provide structures for the synthesis of the lagging-strand by mechanisms such as template switching. Thus, the coupling of recombination with rolling circle replication and possibly, template switching, may have been selected as an economic replication mode to accommodate the reductive evolution of mitochondria. Such a replication mode spares the need for conventional replicative components, including those required for origin recognition/remodelling, RNA primer synthesis and lagging-strand processing.

Characterization and phylogenetic analysis of the mitochondrial genome of Glarea lozoyensis indicates high diversity within the order Helotiales

BACKGROUND

Glarea lozoyensis is a filamentous fungus used for the industrial production of non-ribosomal peptide pneumocandin B0. In the scope of a whole genome sequencing the complete mitochondrial genome of the fungus has been assembled and annotated. It is the first one of the large polyphyletic Helotiaceae family. A phylogenetic analysis was performed based on conserved proteins of the oxidative phosphorylation system in mitochondrial genomes.

RESULTS

The total size of the mitochondrial genome is 45,038 bp. It contains the expected 14 genes coding for proteins related to oxidative phosphorylation,two rRNA genes, six hypothetical proteins, three intronic genes of which two are homing endonucleases and a ribosomal protein rps3. Additionally there is a set of 33 tRNA genes. All genes are located on the same strand. Phylogenetic analyses based on concatenated mitochondrial protein sequences confirmed that G. lozoyensis belongs to the order of Helotiales and that it is most closely related to Phialocephala subalpina. However, a comparison with the three other mitochondrial genomes known from Helotialean species revealed remarkable differences in size, gene content and sequence. Moreover, it was found that the gene order found in P. subalpina and Sclerotinia sclerotiorum is not conserved in G. lozoyensis.

CONCLUSION

The arrangement of genes and other differences found between the mitochondrial genome of G. lozoyensis and those of other Helotiales indicates a broad genetic diversity within this large order. Further mitochondrial genomes are required in order to determine whether there is a continuous transition between the different forms of mitochondrial genomes or G. lozoyensis belongs to a distinct subgroup within Helotiales.

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

We determined the complete sequence of 71 355-bp-long mitochondrial genome from Saccharomyces paradoxus entirely by direct sequencing of purified mitochondrial DNA (mtDNA). This mtDNA possesses the same features as its close relative Saccharomyces cerevisiae - A + T content 85.9%, set of genes coding for the three components of cytochrome oxidase, cytochrome b, three subunits of ATPase, both ribosomal subunits, gene for ribosomal protein, rnpB gene, tRNA package (24) and yeast genetic code. Genes are interrupted by nine group I and group II introns, two of which are in positions unknown in S. cerevisiae, but recognized in Saccharomyces pastorianus. The gene products are related to S. cerevisiae, and the identity of amino acid residues varies from 100% for cox2 to 83% for rps3. The remarkable differences from S. cerevisiae are (1) different gene order (translocation of trnF-trnT1-trnV-cox3-trnfM-rnpb-trnP and transposition of trnW-rns), (2) occurrence of two unusual GI introns, (3) eight active ori elements, and (4) reduced number of GC clusters and divergent intergenic spacers. Despite these facts, the sequenced S. paradoxus mtDNA introduced to S. cerevisiae was able to support the respiratory function to the same extent as the original mtDNAs.

The mitochondrial genome of Moniliophthora roreri, the frosty pod rot pathogen of cacao

In this study, we report the sequence of the mitochondrial (mt) genome of the Basidiomycete fungus Moniliophthora roreri, which is the etiologic agent of frosty pod rot of cacao (Theobroma cacao L.). We also compare it to the mtDNA from the closely-related species Moniliophthora perniciosa, which causes witches' broom disease of cacao. The 94 Kb mtDNA genome of M. roreri has a circular topology and codes for the typical 14 mt genes involved in oxidative phosphorylation. It also codes for both rRNA genes, a ribosomal protein subunit, 13 intronic open reading frames (ORFs), and a full complement of 27 tRNA genes. The conserved genes of M. roreri mtDNA are completely syntenic with homologous genes of the 109 Kb mtDNA of M. perniciosa. As in M. perniciosa, M. roreri mtDNA contains a high number of hypothetical ORFs (28), a remarkable feature that make Moniliophthoras the largest reservoir of hypothetical ORFs among sequenced fungal mtDNA. Additionally, the mt genome of M. roreri has three free invertron-like linear mt plasmids, one of which is very similar to that previously described as integrated into the main M. perniciosa mtDNA molecule. Moniliophthora roreri mtDNA also has a region of suspected plasmid origin containing 15 hypothetical ORFs distributed in both strands. One of these ORFs is similar to an ORF in the mtDNA gene encoding DNA polymerase in Pleurotus ostreatus. The comparison to M. perniciosa showed that the 15 Kb difference in mtDNA sizes is mainly attributed to a lower abundance of repetitive regions in M. roreri (5.8 Kb vs 20.7 Kb). The most notable differences between M. roreri and M. perniciosa mtDNA are attributed to repeats and regions of plasmid origin. These elements might have contributed to the rapid evolution of mtDNA. Since M. roreri is the second species of the genus Moniliophthora whose mtDNA genome has been sequenced, the data presented here contribute valuable information for understanding the evolution of fungal mt genomes among closely-related species.

Phylogenetic and biogeographic implications inferred by mitochondrial intergenic region analyses and ITS1-5.8S-ITS2 of the entomopathogenic fungi Beauveria bassiana and B. brongniartii

BACKGROUND

The entomopathogenic fungi of the genus Beauveria are cosmopolitan with a variety of different insect hosts. The two most important species, B. bassiana and B. brongniartii, have already been used as biological control agents of pests in agriculture and as models for the study of insect host - pathogen interactions. Mitochondrial (mt) genomes, due to their properties to evolve faster than the nuclear DNA, to contain introns and mobile elements and to exhibit extended polymorphisms, are ideal tools to examine genetic diversity within fungal populations and genetically identify a species or a particular isolate. Moreover, mt intergenic region can provide valuable phylogenetic information to study the biogeography of the fungus.

RESULTS

The complete mt genomes of B. bassiana (32,263 bp) and B. brongniartii (33,920 bp) were fully analysed. Apart from a typical gene content and organization, the Beauveria mt genomes contained several introns and had longer intergenic regions when compared with their close relatives. The phylogenetic diversity of a population of 84 Beauveria strains -mainly B. bassiana (n = 76) - isolated from temperate, sub-tropical and tropical habitats was examined by analyzing the nucleotide sequences of two mt intergenic regions (atp6-rns and nad3-atp9) and the nuclear ITS1-5.8S-ITS2 domain. Mt sequences allowed better differentiation of strains than the ITS region. Based on mt and the concatenated dataset of all genes, the B. bassiana strains were placed into two main clades: (a) the B. bassiana s. l. and (b) the "pseudobassiana". The combination of molecular phylogeny with criteria of geographic and climatic origin showed for the first time in entomopathogenic fungi, that the B. bassiana s. l. can be subdivided into seven clusters with common climate characteristics.

CONCLUSIONS

This study indicates that mt genomes and in particular intergenic regions provide molecular phylogeny tools that combined with criteria of geographic and climatic origin can subdivide the B. bassiana s.l. entomopathogenic fungi into seven clusters with common climate characteristics.

The complete mitochondrial genome of yarrowia lipolytica

We here report the complete nucleotide sequence of the 47.9 kb mitochondrial (mt) genome from the obligate aerobic yeast Yarrowia lipolytica. It encodes, all on the same strand, seven subunits of NADH: ubiquinone oxidoreductase (ND1-6, ND4L), apocytochrome b (COB), three subunits of cytochrome oxidase (COX1, 2, 3), three subunits of ATP synthetase (ATP6, 8 and 9), small and large ribosomal RNAs and an incomplete set of tRNAs. The Y. lipolytica mt genome is very similar to the Hansenula wingei mt genome, as judged from blocks of conserved gene order and from sequence homology. The extra DNA in the Y. lipolytica mt genome consists of 17 group 1 introns and stretches of A+Trich sequence, interspersed with potentially transposable GC clusters. The usual mould mt genetic code is used. Interestingly, there is no tRNA able to read CGN (arginine) codons. CGN codons could not be found in exonic open reading frames, whereas they do occur in intronic open reading frames. However, several of the intronic open reading frames have accumulated mutations and must be regarded as pseudogenes. We propose that this may have been triggered by the presence of untranslatable CGN codons. This sequence is available under EMBL Accession No. AJ307410.

Recent dermatophyte divergence revealed by comparative and phylogenetic analysis of mitochondrial genomes

Background

Dermatophytes are fungi that cause superficial infections of the skin, hair, and nails. They are the most common agents of fungal infections worldwide. Dermatophytic fungi constitute three genera, Trichophyton, Epidermophyton, and Microsporum, and the evolutionary relationships between these genera are epidemiologically important. Mitochondria are considered to be of monophyletic origin and mitochondrial sequences offer many advantages for phylogenetic studies. However, only one complete dermatophyte mitochondrial genome (E. floccosum) has previously been determined.

Results

The complete mitochondrial DNA sequences of five dermatophyte species, T. rubrum (26,985 bp), T. mentagrophytes (24,297 bp), T. ajelloi (28,530 bp), M. canis (23,943 bp) and M. nanum (24,105 bp) were determined. These were compared to the E. floccosum sequence. Mitochondrial genomes of all 6 species were found to harbor the same set of genes arranged identical order indicating that these dermatophytes are closely related. Genome size differences were largely due to variable lengths of non-coding intergenic regions and the presence/absence of introns. Phylogenetic analyses based on complete mitochondrial genomes reveals that the divergence of the dermatophyte clade was later than of other groups of pathogenic fungi.

Conclusion

This is the first systematic comparative genomic study on dermatophytes, a highly conserved and recently-diverged lineage of ascomycota fungi. The data reported here provide a basis for further exploration of interrelationships between dermatophytes and will contribute to the study of mitochondrial evolution in higher fungi.

Complete mitochondrial genome sequence of the wine yeast Candida zemplinina: intraspecies distribution of a novel group-IIB1 intron with eubacterial affiliations

The mtDNA of the ascomycetous wine yeast Candida zemplinina is a circularly mapping genome of 23,114 bp. It contains 35 genes coding for the seven basic subunits of oxidative phosporylation found in yeasts (the genes encoding for NADH oxidoreductase subunits are absent), the ribosomal protein Var1, two rRNAs and 25 tRNA genes. Although protein phylogenetic analysis showed a divergent mitochondrial genome, several traits appeared preserved. The conserved gene blocks between the mtDNAs of C. zemplinina and Candida glabrata were maintained and changes in gene order and putative promoters were due to restricted genome reshuffling. New heterogeneous hairpin elements were identified scattered throughout cox1 introns. The large subunit rRNA gene harboured the first group-IIB1 intron containing a putative active reverse transcriptase (RT) in mitochondrial genomes of fungi. Phylogenetic analysis of the RT protein confirmed its closer relationship to eubacterial intronic RTs, while being only distantly related to all other fungal mitochondrial group-II introns and RTs. The findings point towards an early migration event of a eubacterial group-II intron to the mitochondrial genome of C. zemplinina.

A functional core of the mitochondrial genome maintenance protein Mgm101p in Saccharomyces cerevisiae determined with a temperature-conditional allele

Analysis of Mgm101p isolated from mitochondria shows that the mature protein of 27.6 kDa lacks 22 amino acids from the N-terminus. This mitochondrial targeting sequence has been incorporated in the design of oligonucleotides used to determine a functional core of Mgm101p. Progressive deletions, although retaining the targeting sequence, reveal that 76 N-terminal and six C-terminal amino acids of Mgm101p can be removed without altering the ability to complement an mgm101-1(ts) temperature-sensitive mutant. However, this active core is unable to complement mgm101 null mutants, suggesting that the Mgm101p might need to form a dimer or multimer to be functional in vivo. The active core, enriched in basic residues, contains 165 amino acids with a pI of 9.2. Alignment with 22 Mgm101p sequences from other lower eukaryotes shows that a number of amino acids are highly conserved in this region. Random mutagenesis confirms that certain critical amino acids required for function are invariant across the 23 proteins. Searches in the PFAM database revealed a low level of structural similarity between the active core and the Rad52 protein family.

The mitochondrial genome of the wine yeast Hanseniaspora uvarum: a unique genome organization among yeast/fungal counterparts

The complete sequence of the apiculate wine yeast Hanseniaspora uvarum mtDNA has been determined and analysed. It is an extremely compact linear molecule containing the shortest functional region ever found in fungi (11 094 bp long), flanked by Type 2 telomeric inverted repeats. The latter contained a 2704-bp-long subterminal region and tandem repeats of 839-bp units. In consequence, a population of mtDNA molecules that differed at the number of their telomeric reiterations was detected. The functional region of the mitochondrial genome coded for 32 genes, which included seven subunits of respiratory complexes and ATP synthase (the genes encoding for NADH oxidoreductase subunits were absent), two rRNAs and 23 tRNA genes which recognized codons for all amino acids. A single intron interrupted the cytochrome oxidase subunit 1 gene. A number of reasons contributed towards its strikingly small size, namely: (1) the remarkable size reduction (by >40%) of the rns and rnl genes; (2) that most tRNA genes and five of the seven protein-coding genes were the shortest among known yeast homologs; and (3) that the noncoding regions were restricted to 5.1% of the genome. In addition, the genome showed multiple changes in the orientation of transcription and the gene order differed drastically from other yeasts. When all protein coding gene sequences were considered as one unit and were compared with the corresponding molecules from all other complete mtDNAs of yeasts, the phylogenetic trees constructed robustly supported its placement basal to the yeast species of the 'Saccharomyces complex', demonstrating the advantage of this approach over single-gene or multigene approaches of unlinked genes.

Complete DNA sequences of the mitochondrial genomes of the pathogenic yeasts Candida orthopsilosis and Candida metapsilosis: insight into the evolution of linear DNA genomes from mitochondrial telomere mutants

We determined complete mitochondrial DNA sequences of the two yeast species, Candida orthopsilosis and Candida metapsilosis, and compared them with the linear mitochondrial genome of their close relative, C.parapsilosis. Mitochondria of all the three species harbor compact genomes encoding the same set of genes arranged in the identical order. Differences in the length of these genomes result mainly from the presence/absence of introns. Multiple alterations were identified also in the sequences of the ribosomal and transfer RNAs, and proteins. However, the most striking feature of C.orthopsilosis and C.metapsilosis is the existence of strains differing in the molecular form of the mitochondrial genome (circular-mapping versus linear). Their analysis opens a unique window for understanding the role of mitochondrial telomeres in the stability and evolution of molecular architecture of the genome. Our results indicate that the circular-mapping mitochondrial genome derived from the linear form by intramolecular end-to-end fusions. Moreover, we suggest that the linear mitochondrial genome evolved from a circular-mapping form present in a common ancestor of the three species and, at the same time, the emergence of mitochondrial telomeres enabled the formation of linear monomeric DNA forms. In addition, comparison of isogenic C.metapsilosis strains differing in the form of the organellar genome suggests a possibility that, under some circumstances, the linearity and/or the presence of telomeres provide a competitive advantage over a circular-mapping mitochondrial genome.

The complete mitochondrial genome of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae: gene order and trn gene clusters reveal a common evolutionary course for all Sordariomycetes, while intergenic regions show variation

The mitochondrial genome (mtDNA) of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae, with a total size of 24,673 bp, was one of the smallest known mtDNAs of Pezizomycotina. It contained the 14 typical genes coding for proteins related to oxidative phosphorylation, the two rRNA genes, a single intron that harbored an intronic ORF coding for a putative ribosomal protein (rps) within the large rRNA gene (rnl), and a set of 24 tRNA genes which recognized codons for all amino acids, except proline and valine. Gene order comparison with all known mtDNAs of Sordariomycetes illustrated a highly conserved genome organization for all the protein- and rRNA-coding genes, as well as three clusters of tRNA genes. By considering all mitochondrial essential protein-coding genes as one unit a phylogenetic study of these small genomes strongly supported the common evolutionary course of Sordariomycetes (100% bootstrap support) and highlighted the advantages of analyzing small genomes (mtDNA) over single genes. In addition, comparative analysis of three intergenic regions demonstrated sequence variability that can be exploited for intra- and inter-specific identification of Metarhizium.

Sequence analysis of three mitochondrial DNA molecules reveals interesting differences among Saccharomyces yeasts

The complete sequences of mitochondrial DNA (mtDNA) from the two budding yeasts Saccharomyces castellii and Saccharomyces servazzii, consisting of 25 753 and 30 782 bp, respectively, were analysed and compared to Saccharomyces cerevisiae mtDNA. While some of the traits are very similar among Saccharomyces yeasts, others have highly diverged. The two mtDNAs are much more compact than that of S.cerevisiae and contain fewer introns and intergenic sequences, although they have almost the same coding potential. A few genes contain group I introns, but group II introns, otherwise found in S.cerevisiae mtDNA, are not present. Surprisingly, four genes (ATP6, COX2, COX3 and COB) in the mtDNA of S.servazzii contain, in total, five +1 frameshifts. mtDNAs of S.castellii, S.servazzii and S.cerevisiae contain all genes on the same strand, except for one tRNA gene. On the other hand, the gene order is very different. Several gene rearrangements have taken place upon separation of the Saccharomyces lineages, and even a part of the transcription units have not been preserved. It seems that the mechanism(s) involved in the generation of the rearrangements has had to ensure that all genes stayed encoded by the same DNA strand.

The global genetic structure of the wheat pathogen Mycosphaerella graminicola is characterized by high nuclear diversity, low mitochondrial diversity, regular recombination, and gene flow

A total of 1673 Mycosphaerella graminicola strains were assayed for DNA fingerprints and restriction fragment length polymorphism (RFLP) markers in the nuclear and mitochondrial genomes. The isolates were collected from 17 wheat fields located in 11 countries on five continents over a six year period (1989-1995). Our results indicate that genetic diversity in the nuclear genome of this fungus was high for all but three of the field populations surveyed and that populations sampled from different continents had similar frequencies for the most common RFLP alleles. Hierarchical analysis revealed that more than 90% of global gene diversity was distributed within a wheat field, while approximately 5% of gene diversity was distributed among fields within regions and approximately 3% was distributed among regions on different continents. These findings suggest that gene flow has occurred on a global scale. On average, each leaf was colonized by a different nuclear genotype. In contrast, only seven mtDNA haplotypes were detected among the 1673 isolates and the two most common mtDNA haplotypes represented approximately 93% of the world population, consistent with a selective sweep. Analysis of multilocus associations indicated that all field populations were in gametic equilibrium, suggesting that sexual recombination is a regular occurrence globally.

A new PCR-based method for monitoring inoculated wine fermentations

A new PCR-based method has been developed to monitor inoculated wine fermentations. The method is based on the variation in the number and position of introns in the mitochondrial gene COX1. Oligonucleotide primers homologous to the regions flanking the Saccharomyces cerevisiae COX1 introns have been designed and tested for S. cerevisiae wine yeast strain differentiation. Four primers were selected for their subsequent use in a multiplex PCR reaction and have proved to be very effective in uncovering polymorphism in natural and commercial yeast strains. An important point is that the speed and simplicity of the technique, which does not require the isolation of DNA, allows early detection of the starter yeast strain throughout the fermentation process. The main advantage for the wineries is that the must sample can be used directly for the PCR reaction obtaining very fast results (in approximately 8 h). This allows the wine industries to intervene quickly if necessary.

Mechanisms influencing the evolution of resistance to Qo inhibitor fungicides

Fungicides inhibiting the mitochondrial respiration of plant pathogens by binding to the cytochrome bc1 enzyme complex (complex III) at the Qo site (Qo inhibitors, QoIs) were first introduced to the market in 1996. After a short time period, isolates resistant to QoIs were detected in field populations of a range of important plant pathogens including Blumeria graminis Speer f sp tritici, Sphaerotheca fuliginea (Schlecht ex Fr) Poll, Plasmopara viticola (Berk & MA Curtis ex de Bary) Berl & de Toni, Pseudoperonospora cubensis (Berk & MA Curtis) Rost, Mycosphaerella fijiensis Morelet and Venturia inaequalis (Cooke) Wint. In most cases, resistance was conferred by a point mutation in the mitochondrial cytochrome b (cyt b) gene leading to an amino-acid change from glycine to alanine at position 143 (G143A), although additional mutations and mechanisms have been claimed in a number of organisms. Transformation of sensitive protoplasts of M fijiensis with a DNA fragment of a resistant M fijiensis isolate containing the mutation yielded fully resistant transformants, demonstrating that the G143A substitution may be the most powerful transversion in the cyt b gene conferring resistance. The G143A substitution is claimed not to affect the activity of the enzyme, suggesting that resistant individuals may not suffer from a significant fitness penalty, as was demonstrated in B graminis f sp tritici. It is not known whether this observation applies also for other pathogen species expressing the G143A substitution. Since fungal cells contain a large number of mitochondria, early mitotic events in the evolution of resistance to QoIs have to be considered, such as mutation frequency (claimed to be higher in mitochondrial than nuclear DNA), intracellular proliferation of mitochondria in the heteroplasmatic cell stage, and cell to cell donation of mutated mitochondria. Since the cyt b gene is located in the mitochondrial genome, inheritance of resistance in filamentous fungi is expected to be non-Mendelian and, therefore, in most species uniparental. In the isogamous fungus B graminis f sp tritici, crosses of sensitive and resistant parents yielded cleistothecia containing either sensitive or resistant ascospores and the segregation pattern for resistance in the F1 progeny population was 1:1. In the anisogamous fungus V inaequalis, donation of resistance was maternal and the segregation ratio 1:0. In random mating populations, the sex ratio (mating type distribution) is generally assumed to be 1:1. Therefore, the overall proportion of sensitive and resistant individuals in unselected populations is expected to be 1:1. Evolution of resistance to QoIs will depend mainly on early mitotic events; the selection process for resistant mutants in populations exposed to QoI treatments may follow mechanisms similar to those described for resistance controlled by single nuclear genes in other fungicide classes. It will remain important to understand how the mitochondrial nature of QoI resistance and factors such as mutation, recombination, selection and migration might influence the evolution of QoI resistance in different plant pathogens.

Analysis of restriction profiles of mitochondrial DNA from Sporothrix schenckii

Mitochondrial DNA (mtDNA) diversity was analyzed in 42 clinical isolates of Sporothrix schenckii from Mexico (n = 29), Guatemala (n = 4) and Colombia (n = 9). Based on HaeIII restriction digestion profiles, the isolates were classified into eight types. In addition to 24 mtDNA types previously reported in another study, 6 new types were found in this study. Most of the strains belong to type 14 and type 30, the former restricted to Mexico, whereas the latter was distributed in Mexico, Guatemala and Colombia. The new types (25-30) were identified in Mexico, Guatemala and Colombia. Restriction-fragment length polymorphism in mtDNA of S. schenckii revealed high levels of genetic variation attributable to differences in restriction sites as well as in mtDNA size. Based on genetic distances S. schenckii types were clustered into two main groups.

Diversity in organization and the origin of gene orders in the mitochondrial DNA molecules of the genus Saccharomyces

Sequencing of the Saccharomyces cerevisiae nuclear and mitochondrial genomes provided a new background for studies on the evolution of the genomes. In this study, mitochondrial genomes of a number of Saccharomyces yeasts were mapped by restriction enzyme analysis, the orders of the genes were determined, and two of the genes were sequenced. The genome organization, i.e., the size, presence of intergenic sequences, and gene order, as well as polymorphism within the coding regions, indicate that Saccharomyces mtDNA molecules are dynamic structures and have undergone numerous changes during their evolution. Since the separation and sexual isolation of different yeast lineages, the coding parts have been accumulating point mutations, presumably in a linear manner with the passage of time. However, the accumulation of other changes may not have been a simple function of time. Larger mtDNA molecules belonging to Saccharomyces sensu stricto yeasts have acquired extensive intergenic sequences, including guanosine-cytosine-rich clusters, and apparently have rearranged the gene order at higher rates than smaller mtDNAs belonging to the Saccharomyces sensu lato yeasts. While within the sensu stricto group transposition has been a predominant mechanism for the creation of novel gene orders, the sensu lato yeasts could have used both transposition- and inversion-based mechanisms.

Molecular identification of wine yeasts at species or strain level: a case study with strains from two vine-growing areas of Greece

The composition of wine yeast populations, present during spontaneous fermentation of musts from two wine-producing areas of Greece (Amyndeon and Santorini) and followed for two consecutive years, were studied using a range of molecular techniques. Internal Transcribed Spacer (ITS) ribotyping was convincingly applied for yeast species identification, proving its usefulness as a reliable tool for the rapid characterization of species composition in yeast population studies. Restriction Fragment Length Polymorphism (RFLP) of mitochondrial DNA (mtDNA) was shown to be a convenient criterion for the detection of intraspecies genetic diversity of both Saccharomyces and non-Saccharomyces isolate populations. Similarly, polymorphism of amplified delta interspersed element sequences provided an additional criterion for S. cerevisiae strain differentiation. Comparative analysis of S. cerevisiae genetic diversity, using mtDNA restriction patterns and delta-amplification profiles, showed a similar discriminative power of the two techniques. However, by combining these approaches it was possible to distinguish/characterize strains of the same species and draw useful conclusions about yeast diversity during alcoholic fermentation. The most significant findings in population dynamics of yeasts in the spontaneous fermentations were (i) almost complete absence of non-S.cerevisiae species from fermentations of must originating from the island Santorini, (ii) a well recorded strain polymorphism in populations of non-Saccharomyces species originating from Amyndeon and (iii) an unexpected polymorphism concerning S. cerevisiae populations, much greater than ever reported before in similar studies with wine yeasts of other geographical regions.

The petite mutation in yeasts: 50 years on

Fifty years ago it was reported that baker's yeast, Saccharomyces cerevisiae, can form "petite colonie" mutants when treated with the DNA-targeting drug acriflavin. To mark the jubilee of studies on cytoplasmic inheritance, a review of the early work will be presented together with some observations on current developments. The primary emphasis is to address the questions of how loss of mtDNA leads to lethality (rho 0-lethality) in petite-negative yeasts and how S. cerevisiae tolerates elimination of mtDNA. Recent investigation have revealed that rho 0-lethality can be suppressed by specific mutations in the alpha, beta, and gamma subunits of the mitochondrial F1-ATPase of the petite-negative yeast Kluyveromyces lactis and by the nuclear ptp alleles in Schizosaccharomyces pombe. In contrast, inactivation of genes coding for F1-ATPase alpha and beta subunits and disruption of AAC2, PGS1/PEL1, and YME1 genes in S. cerevisiae convert this petite-positive yeast into a petite-negative form. Studies on nuclear genes affecting dependence on mtDNA have provided important insight into the functions provided by the mitochondrial genome and the maintenance of structural and functional integrity of the mitochondrial inner membrane.

Extrachromosomal genetic elements in fungi

All fungi like eukaryotes possess mitochondria, which are the sites of the oxydative phosphorylation. As eukaryote evolution depends on oxygenic atmosphere, these organisms are primarily aerobic. Except a small group of strict anaerobes (those which lost the capacity of oxydative pathways living in special milieu in association with rumen of grass-eating animals) all fungi can utilize various compounds as carbon sources via oxidative phosphorylation pathways resulting in high energy yield. Certain groups of fungi – i.e. most of the yeasts – under anaerobe conditions, are able to supply themselves with lower levels of fermentation energy, too exhibiting a slow growing capacity utilizing the same amount of carbon source. The mutation of mitochondrial genome or mitochondrial functions encoded by nuclear genes of these fungi might result in a so-calledpetitephenotype producing small colonies on solid media due to their slow growing capacity. These mutants can utilize only fermentable carbon sources. Filamentous fungi have only limited possibilities to produce such phenotypes. ExceptZygomycetes(where the shortage of oxygen induces dimorphic transitions) filamentous fungi can grow and develop their vegetative and sexual reproductive structures only in aerobe milieu. However amongNeurosporaspecies there are several mitochondrial mutations resulting in morphological phenotypes. These are due to the lower energy level provided by the reduced capacity of cytochrome-oxidase enzymes. These mutants (e.g.pokystopper) can be considered aspetiteanalogues. The complete loss of mitochondrial functions – such asrho zerocharacter in yeast – cannot be survived by filamentous fungi.Podospora anserinaand some of its close relatives exhibit a so-calledsenescencephenotype, which means that the growing hyphae in the youngest part of the colonies stop growing and start to die within a short period of time. This phenomenon – discussed below – is also connected to reduced function of mitochondria.The first part of this paper gives a short overview of the genetic organization of mitochondria of fungi, based on the most recent data of three filamentous fungi:Aspergillus nidulansNeurospora crassaandPodospora anserina. Their data are compared to those of the well-characterizedSaccharomyces cerevisiae. In the second part we summarize what we know about other extrachrosomal elements, such as DNA plasmids of various origins and structures, and dsRNAs or virus like particles (VLP). Also discussed are their roles and/or putative functions in the life of the fungi.

Sex and conflict

Evolutionary conflict occurs when the deterministic spread of an allele lowers the fitness either of its bearer or of other individuals in the population, leading to selection for suppressors. Sex promotes conflict because associations between alleles are temporary. Differing selection on males and females, sexual selection, and differences in transmission patterns between classes of nuclear and cytoplasmic genes can all give rise to conflict. Inert Y chromosomes, uniparental inheritance of cytoplasmic genes, mating strains and sexes, and many features of sexual behavior may have evolved in part as a result of evolutionary conflict. Estimates of its quantitative importance, however, are still needed.

Genome structure and gene content in protist mitochondrial DNAs

Although the collection of completely sequenced mitochondrial genomes is expanding rapidly, only recently has a phylogenetically broad representation of mtDNA sequences from protists (mostly unicellular eukaryotes) become available. This review surveys the 23 complete protist mtDNA sequences that have been determined to date, commenting on such aspects as mitochondrial genome structure, gene content, ribosomal RNA, introns, transfer RNAs and the genetic code and phylogenetic implications. We also illustrate the utility of a comparative genomics approach to gene identification by providing evidence that orfB in plant and protist mtDNAs is the homolog of atp8 , the gene in animal and fungal mtDNA that encodes subunit 8 of the F0portion of mitochondrial ATP synthase. Although several protist mtDNAs, like those of animals and most fungi, are seen to be highly derived, others appear to be have retained a number of features of the ancestral, proto-mitochondrial genome. Some of these ancestral features are also shared with plant mtDNA, although the latter have evidently expanded considerably in size, if not in gene content, in the course of evolution. Comparative analysis of protist mtDNAs is providing a new perspective on mtDNA evolution: how the original mitochondrial genome was organized, what genes it contained, and in what ways it must have changed in different eukaryotic phyla.

Transcript mapping and processing of mitochondrial RNA in the chlorophyte alga Prototheca wickerhamii

The detailed transcript map of the circular 55328 bp mitochondrial (mt) genome from the colourless chlorophycean alga Prototheca wickerhamii has been determined. On each half of this genome the genes are encoded only on one DNA strand, forming transcriptional units comprising variable numbers of genes. With the exception of four genes coding for ribosomal proteins, transcripts of the three rRNA genes and all protein-coding genes have been detected by both northern analysis and primer extension experiments. Polycistronic transcripts of protein coding and tRNA genes were verified by northern analyses, primer extension and RNAse mapping experiments. The 5' and 3' ends of different RNA species are often located in close proximity to putative stem-loop structures and some 5' termini of mRNAs coincide with the 3' end of tRNAs located immediately upstream. Transcript mapping in a putative promoter region revealed two different possible transcription initiation sites; no significant sequence homology to putative mt promoters from higher plants could be found. In addition, two out of three group I introns residing in the cox1 gene were found to be self-splicing in vitro under reaction conditions developed for related mt introns from a filamentous fungus. Mitochondrial gene expression of P. wickerhamii and of filamentous fungi has several features in common, such as intron splicing and the processing of longer polycistronic transcripts. The similarities in RNA maturation between higher-plant and P. wickerhamii mitochondria are less pronounced, since plants rarely use tRNAs as processing signals for their relatively short mitochondrial co-transcripts.

Fungal origin by horizontal transfer of a plant mitochondrial group I intron in the chimeric CoxI gene of Peperomia

We present phylogenetic evidence that a group I intron in an angiosperm mitochondrial gene arose recently by horizontal transfer from a fungal donor species. A 1,716-bp fragment of the mitochondrial coxI gene from the angiosperm Peperomia polybotrya was amplified via the polymerase chain reaction and sequenced. Comparison to other coxI genes revealed a 966-bp group I intron, which, based on homology with the related yeast coxI intron aI4, potentially encodes a 279-amino-acid site-specific DNA endonuclease. This intron, which is believed to function as a ribozyme during its own splicing, is not present in any of 19 coxI genes examined from other diverse vascular plant species. Phylogenetic analysis of intron origin was carried out using three different tree-generating algorithms, and on a variety of nucleotide and amino acid data sets from the intron and its flanking exon sequences. These analyses show that the Peperomia coxI gene intron and exon sequences are of fundamentally different evolutionary origin. The Peperomia intron is more closely related to several fungal mitochondrial introns, two of which are located at identical positions in coxI, than to identically located coxI introns from the land plant Marchantia and the green alga Prototheca. Conversely, the exon sequence of this gene is, as expected, most closely related to other angiosperm coxI genes. These results, together with evidence suggestive of co-conversion of exonic markers immediately flanking the intron insertion site, lead us to conclude that the Peperomia coxI intron probably arose by horizontal transfer from a fungal donor, using the double-strand-break repair pathway. The donor species may have been one of the symbiotic mycorrhizal fungi that live in close obligate association with most plants.

Linear mitochondrial genome organization in vivo in the genus Pythium

Pulsed-field gel electrophoresis (PFGE) of isolates of Pythium oligandrum with linear mitochondrial genomes revealed a distinct band in ethidium bromide-stained gels similar in size to values estimated by restriction mapping of mitochondrial DNA (mtDNA). Southern analysis confirmed that these bands were mtDNA and indicated that linear genomes were present in unit-length size as well as multimers. Isolates of this species with circular mtDNA restriction maps also had low levels of linear mono- and multimers visualized by Southern analysis of PFGE gels. Examination of 17 additional species revealed similar results; three species had distinct linear mtDNA bands in ethidium bromide-stained gels while the remainder had linear mono- and multi-mers in lower amounts detected only by Southern analysis. Sequence analysis of an isolate of P. oligandrum with a primarily circular mitochondrial genomic map and a low amount of linear molecules revealed that the small unique region of the circular map (which corresponded to the terminal region of linear genomes) was flanked by palindromic intrastrand complementary sequences separated by a unique 194-bp sequence. Sequences with similarity to ATPase9 coding regions from other organisms were located adjacent to this region. Sequences with similarity to mitochondrial origins of replication and autonomously replicating sequences were also located in this region: their potential involvement in the generation of linear molecules is discussed.

The complete mitochondrial DNA sequence of Hansenula wingei reveals new characteristics of yeast mitochondria

The complete 27,694-bp mitochondrial (mt) DNA sequence of Hansenula wingei, which is a typical budding yeast and contains circular mitochondrial DNA, has been determined. The mt sequence contains genes encoding large and small ribosomal RNAs, 25 tRNAs, three subunits of cytochrome c oxidase (subunits 1, 2 and 3), three subunits of ATPase (subunits 6, 8 and 9), apocytochrome b, seven subunits of NADH dehydrogenase (subunits 1, 2, 3, 4, 4L, 5 and 6), and a ribosomal protein, VAR1. The VAR1 gene is considered to be a typical yeast type. This is consistent with data on DNA and the deduced amino-acid sequence homology comparisons of genes ubiquitous in yeast and fungi. However, we have identified seven genes encoding NADH dehydrogenase subunits, which are not found in other yeast mitochondrial genomes, thus placing the H. wingei mitochondrial genome in a unique position. In addition the H. wingei mitochondrial genome also encodes one tRNA pseudogene and one short unidentified ORF. The genome is compact with only two introns both of which contain an ORF. One intron lies in the large rRNA gene while the other is situated in the cytochrome c oxidase subunit-1 gene. The conserved nonanucleotide motif (A/T)TATAAG (T/A)(A/T), which is a transcription start signal in Saccharomyces cerevisiae mitochondria, has also been found in the H. wingei mitochondrial genome. The codon assignments for ATA and CTN in H. wingei mitochondria are different from those in S. cerevisiae mitochondria. These results indicate a unique and novel structure for the H. wingei mitochondrial genome in terms of characteristics which are typical for both yeast and for filamentous fungi. This is the first complete mt DNA sequence report in yeast.

Sexuality of mitochondria: fusion, recombination, and plasmids

Mitochondrial fusion, recombination, and mobile genetic elements, which are essential for mitochondrial sexuality, are well established in various organisms. The recombination of mitochondrial DNA (mtDNA) depends upon fusion between parental mitochondria, and between their mtDNA-containing areas (mt-nuclei), to allow pairing between the parental mtDNAs. Such mitochondrial fusion followed by recombination may be called "mitochondrial sex." We have identified a novel mitochondrial plasmid named mF. This plasmid is apparently responsible for promoting mitochondrial fusion and crosses over with mtDNA in successive sexual crosses with mF- strains. Only in mF+ strains carrying the mF plasmid did small spherical mitochondria fuse which subsequently underwent fusion between the mt-nuclei that contained the mtDNA derived from individual mitochondria. Several successive mitochondrial divisions followed, accompanied by mt-nuclear divisions. The resulting mitochondria contained recombinant mtDNA with the mF plasmid. Such features remind us also of the bacterial conjugative plasmids such as F plasmid. Therefore, in the final part of this chapter, we discuss the origin of sex and its relationship to the sexuality of mitochondria.

Skewed base compositions, asymmetric transition matrices, and phylogenetic invariants

Evolutionary inference methods that assume equal DNA base compositions and symmetric nucleotide substitution matrices, where these assumptions do not hold, are likely to group species on the basis of similar base compositions rather than true phylogenetic relationships. We propose an invariants-based method for dealing with this problem. An invariant QT of a tree T under a k-state Markov model, where a generalized time parameter is identified with the E edges of T, allows us to recognize whether data on N observed species can be associated with the N terminal vertices of T in the sense of having been generated on T rather than on any other tree with N terminals. The form of the generalized time parameter is a positive determinant matrix in some semigroup S of stochastic matrices. The invariance is with respect to the choice of the set of E matrices in S, one associated with each of the E edges of T. We apply a general "empirical" method of finding invariants of a parametrized functional form. It involves calculating the probability f of all KN data possibilities for each of m sets of E matrices in S to associate with the edges of T, then solving for the parameters using the m equations of form Q(f) = 0. We discuss the problems of finding asymmetric models satisfying the property of semigroup closure, of finding asymmetric models that admit invariants at all, and of the computational complexity of the method. We propose a class of semigroups Sc containing matrices of form [formula: see text] to account for A+T versus G+C asymmetries in DNA base composition. Quadratic invariants are obtained for rooted trees with three and with four terminals. In the latter case the smallest set of algebraically independent invariants is sought. These invariants are applied to data pertaining the fungal evolution and to the origin of mitochondria as bacterial endosymbionts.

Larger rearranged mitochondrial genomes in Dekkera/Brettanomyces yeasts are more closely related than smaller genomes with a conserved gene order

Mitochondrial genomes from yeasts in the Dekkera/Brettanomyces/Eeniella group vary in size from 28 to 101 kb. Mapping of genes has shown that the three smallest genomes, of 28-42 kb, have the same gene order, whereas the three larger mitochondrial DNAs of 57-101 kb are rearranged relative to the smaller molecules and between themselves. To examine the relationships between these genomes, a phylogenetic tree has been constructed by sequence comparison of the mitochondrial-encoded cytochrome oxidase subunit gene (COX2) from the six species. Contrary to expectation, the tree shows that the larger rearranged genomes are more closely related than the smaller mtDNAs. This result indicates that the gene order of the smaller mtDNAs (28-42 kb) is ancestral and that larger mtDNA molecules (57-101 kb) are more prone to rearrangement than smaller forms.