The complete DNA sequence of the mitochondrial genome of Podospora anserina

The complete mitochondrial genome of Penicillium expansum: Insights into the fungal evolution and phylogeny

Penicillium expansum is an important phytopathogenic fungus that causes blue mold disease. In this study, the novel mitochondrial genome of P. expansum was sequenced, assembled, annotated, and compared with the previously published Penicillium mitogenomes. P. expansum mitogenome is composed of circular DNA molecules with a genome size of 25,496 bp. It encodes 16 protein-encoding genes (PCGs), two rRNA genes, and 25 tRNA genes. Comparative analysis with six other Penicillium species revealed that gene length, GC content, AT skew, and GC skew were variable among the core protein-coding genes. The Penicillium species' gene synteny analysis identified several gene rearrangements. Among the core 15 PCGs, atp8 had the lowest K2P genetic distance, which shows that this gene is highly conserved. The Ka/Ks value of most PCGs was less than 1, which shows that these genes have undergone purifying selection. Phylogenetic analysis based on 14 concatenated core mitochondrial genes revealed that P. expansum shares a close relationship with P. solitum. This study served as a first report on the complete mitochondrial genome of P. expansum and its comparative analysis that will contribute to population genetics and rapid evolutionary studies among Penicillium species.

Impact of Mitochondrial Architecture, Function, Redox Homeostasis, and Quality Control on Organismic Aging: Lessons from a Fungal Model System

Significance: Mitochondria are eukaryotic organelles with various essential functions. They are both the source and the targets of reactive oxygen species (ROS). Different branches of a mitochondrial quality control system (mQCS), such as ROS balancing, degradation of damaged proteins, or whole mitochondria, can mitigate the adverse effects of ROS stress. However, the capacity of mQCS is limited. Overwhelming this capacity leads to dysfunctions and aging. Strategies to interfere into mitochondria-dependent human aging with the aim to increase the healthy period of life, the health span, rely on the precise knowledge of mitochondrial functions. Experimental models such as Podospora anserina, a filamentous fungus with a clear mitochondrial aging etiology, proved to be instrumental to reach this goal. Recent Advances: Investigations of the P. anserina mQCS revealed that it is constituted by a complex network of different branches. Moreover, mitochondrial architecture and lipid homeostasis emerged to affect aging. Critical Issues: The regulation of the mQCS is only incompletely understood. Details about the involved signaling molecules and interacting pathways remain to be elucidated. Moreover, most of the currently generated experimental data were generated in well-controlled experiments that do not reflect the constantly changing natural life conditions and bear the danger to miss relevant aspects leading to incorrect conclusions. Future Directions: In P. anserina, the precise impact of redox signaling as well as of molecular damaging for aging remains to be defined. Moreover, natural fluctuation of environmental conditions needs to be considered to generate a realistic picture of aging mechanisms as they developed during evolution.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.

Refining Mitochondrial Intron Classification With ERPIN: Identification Based on Conservation of Sequence Plus Secondary Structure Motifs

Mitochondrial genomes—in particular those of fungi—often encode genes with a large number of Group I and Group II introns that are conserved at both the sequence and the RNA structure level. They provide a rich resource for the investigation of intron and gene structure, self- and protein-guided splicing mechanisms, and intron evolution. Yet, the degree of sequence conservation of introns is limited, and the primary sequence differs considerably among the distinct intron sub-groups. It makes intron identification, classification, structural modeling, and the inference of gene models a most challenging and error-prone task—frequently passed on to an “expert” for manual intervention. To reduce the need for manual curation of intron structures and mitochondrial gene models, computational methods using ERPIN sequence profiles were initially developed in 2007. Here we present a refinement of search models and alignments using the now abundant publicly available fungal mtDNA sequences. In addition, we have tested in how far members of the originally proposed sub-groups are clearly distinguished and validated by our computational approach. We confirm clearly distinct mitochondrial Group I sub-groups IA1, IA3, IB3, IC1, IC2, and ID. Yet, IB1, IB2, and IB4 ERPIN models are overlapping substantially in predictions, and are therefore combined and reported as IB. We have further explored the conversion of our ERPIN profiles into covariance models (CM). Current limitations and prospects of the CM approach will be discussed.

Comparative Analysis of Mitochondrial Genome Features among Four Clonostachys Species and Insight into Their Systematic Positions in the Order Hypocreales

The mycoparasite fungi of Clonostachys have contributed to the biological control of plant fungal disease and nematodes. The Clonostachys fungi strains were isolated from Ophiocordyceps highlandensis, Ophiocordycepsnigrolla and soil, which identified as Clonostachyscompactiuscula, Clonostachysrogersoniana, Clonostachyssolani and Clonostachys sp. To explore the evolutionary relationship between the mentioned species, the mitochondrial genomes of four Clonostachys species were sequenced and assembled. The four mitogenomes consisted of complete circular DNA molecules, with the total sizes ranging from 27,410 bp to 42,075 bp. The GC contents, GC skews and AT skews of the mitogenomes varied considerably. Mitogenomic synteny analysis indicated that these mitogenomes underwent gene rearrangements. Among the 15 protein-coding genes within the mitogenomes, the nad4L gene exhibited the least genetic distance, demonstrating a high degree of conservation. The selection pressure analysis of these 15 PCGs were all below 1, indicating that PCGs were subject to purifying selection. Based on protein-coding gene calculation of the significantly supported topologies, the four Clonostachys species were divided into a group in the phylogenetic tree. The results supplemented the database of mitogenomes in Hypocreales order, which might be a useful research tool to conduct a phylogenetic analysis of Clonostachys. Additionally, the suitable molecular marker was significant to study phylogenetic relationships in the Bionectriaceae family.

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.

Comparative mitogenomics of Agaricomycetes: Diversity, abundance, impact and coding potential of putative open-reading frames

The mitochondrion is an organelle found in eukaryote organisms, and it is vital for different cellular pathways. The mitochondrion has its own DNA molecule and, because its genetic content is relatively conserved, despite the variation of size and structure, mitogenome sequences have been widely used as a promising molecular biomarker for taxonomy and evolution in fungi. In this study, the mitogenomes of two fungal species of Agaricomycetes class, Phellinotus piptadeniae and Trametes villosa, were assembled and annotated for the first time. We used these newly sequenced mitogenomes for comparative analyses with other 55 mitogenomes of Agaricomycetes available in public databases. Mitochondrial DNA (mtDNA) size and content are highly variable and non-coding and intronic regions, homing endonucleases (HEGs), and unidentified ORFs (uORFs) significantly contribute to the total size of the mitogenome. Furthermore, accessory genes (most of them as HEGs) are shared between distantly related species, most likely as a consequence of horizontal gene transfer events. Conversely, uORFs are only shared between taxonomically related species, most probably as a result of vertical evolutionary inheritance. Additionally, codon usage varies among mitogenomes and the GC content of mitochondrial features may be used to distinguish coding from non-coding sequences. Our results also indicated that transposition events of mitochondrial genes to the nuclear genome are not common. Despite the variation of size and content of the mitogenomes, mitochondrial genes seemed to be reliable molecular markers in our time-divergence analysis, even though the nucleotide substitution rates of mitochondrial and nuclear genomes of fungi are quite different. We also showed that many events of mitochondrial gene shuffling probably happened amongst the Agaricomycetes during evolution, which created differences in the gene order among species, even those of the same genus. Altogether, our study revealed new information regarding evolutionary dynamics in Agaricomycetes.

Complete mitochondrial genome sequence of Afla-Guard®, commercially available non-toxigenic Aspergillus flavus

Afla-Guard^® is a commercial non-toxigenic Aspergillus flavus strain used to decrease aflatoxin contamination level in field. Its mitochondrial genome was sequenced, showing that its length is 29,208 bp with typical configuration of Aspergillus mitochondrial genome. 17 SNPs and 27 INDELs were identified by comparing with previous A. flavus mitochondrial genome. Phylogenetic trees present that A. flavus of Afla-Guard^® was clustered with the previous A. flavus mitochondrial genome.

Exploring the Relationship Among Divergence Time and Coding and Non-coding Elements in the Shaping of Fungal Mitochondrial Genomes

The order Hypocreales (Ascomycota) is composed of ubiquitous and ecologically diverse fungi such as saprobes, biotrophs, and pathogens. Despite their phylogenetic relationship, these species exhibit high variability in biomolecules production, lifestyle, and fitness. The mitochondria play an important role in the fungal biology, providing energy to the cells and regulating diverse processes, such as immune response. In spite of its importance, the mechanisms that shape fungal mitogenomes are still poorly understood. Herein, we investigated the variability and evolution of mitogenomes and its relationship with the divergence time using the order Hypocreales as a study model. We sequenced and annotated for the first time Trichoderma harzianum mitochondrial genome (mtDNA), which was compared to other 34 mtDNAs species that were publicly available. Comparative analysis revealed a substantial structural and size variation on non-coding mtDNA regions, despite the conservation of copy number, length, and structure of protein-coding elements. Interestingly, we observed a highly significant correlation between mitogenome length, and the number and size of non-coding sequences in mitochondrial genome. Among the non-coding elements, group I and II introns and homing endonucleases genes (HEGs) were the main contributors to discrepancies in mitogenomes structure and length. Several intronic sequences displayed sequence similarity among species, and some of them are conserved even at gene position, and were present in the majority of mitogenomes, indicating its origin in a common ancestor. On the other hand, we also identified species-specific introns that advocate for the origin by different mechanisms. Investigation of mitochondrial gene transfer to the nuclear genome revealed that nuclear copies of the nad5 are the most frequent while atp8, atp9, and cox3 could not be identified in any of the nuclear genomes analyzed. Moreover, we also estimated the divergence time of each species and investigated its relationship with coding and non-coding elements as well as with the length of mitogenomes. Altogether, our results demonstrated that introns and HEGs are key elements on mitogenome shaping and its presence on fast-evolving mtDNAs could be mostly explained by its divergence time, although the intron sharing profile suggests the involvement of other mechanisms on the mitochondrial genome evolution, such as horizontal transference.

The mitochondrial genome of Morchella importuna (272.2 kb) is the largest among fungi and contains numerous introns, mitochondrial non-conserved open reading frames and repetitive sequences

The complete mitochondrial genome of Morchella importuna, the famous edible and medicinal mushroom, was assembled as a 272,238 bp single circular dsDNA. As the largest mitogenome among fungi, it exhibits several distinct characteristics. The mitogenome of M. importuna encoded 14 core conserved mitochondrial protein-coding genes and 151 mitochondrial non-conserved open reading frames (ncORFs) were predicted, of which 61 were annotated as homing endonuclease genes, and 108 were confirmed to be expressed during the vegetative growth stages of M. importuna. In addition, 34 introns were identified in seven core genes (cob, cox1, cox2, cox3, nad1, nad4 and nad5) and two rRNA genes (rrnS and rrnL) with a length from 383 bp to 7453 bp, and eight large introns with a length range of 2340 bp to 7453 bp contained multiple intronic mtORFs. Moreover, 34 group I (IA, IB, IC1, IC2, ID and derived group I introns) and four group II intron domains were identified for the 34 introns, including five hybrid ones. Furthermore, the M. importuna mitogenome showed the presence of about 18.7% mitogenomic interspersed repeats. These and the aforementioned ncORFs and introns, contributed to the enlarged size of the mitogenome.

The fungal mitochondrial Nad5 pan-genic intron landscape

An intron landscape was prepared for the fungal mitochondrial nad5 gene. A hundred and eighty-eight fungal species were examined and a total of 265 introns were noted to be located in 29 intron insertion sites within the examined nad5 genes. Two hundred and sixty-three introns could be classified as group I types and two group II introns were noted. One additional group II intron module was identified nested within a composite group I intron. Based on features related to RNA secondary structures, introns can be classified into different subtypes and it was observed that intron insertion-sites are biased towards phase 0 and they appear to be specific to an intron type. Intron landscapes could be used as a guide map to predict the location of fungal mtDNA mobile introns, which are composite elements that include a ribozyme component and in some instances open reading frames encoding homing endonucleases or reverse transcriptases and all of these have applications in biotechnology.

The mitochondrial translocase of the inner membrane PaTim54 is involved in defense response and longevity in Podospora anserina

Fungi are very successful microorganisms capable of colonizing virtually any ecological niche where they must constantly cope with competitors including fungi, bacteria and nematodes. We have shown previously that the ascomycete Podopora anserina exhibits Hyphal Interference (HI), an antagonistic response triggered by direct contact of competing fungal hyphae. When challenged with Penicillium chrysogenum, P. anserina produces hydrogen peroxide at the confrontation and kills the hyphae of P. chrysogenum. Here, we report the characterization of the PDC²²¹⁸ mutant affected in HI. When challenged with P. chrysogenum, the PDC²²¹⁸ mutant produces a massive oxidative burst at the confrontation. However, this increased production of hydrogen peroxide is not correlated to increased cell death in P. chrysogenum. Hence, the oxidative burst and cell death in the challenger are uncoupled in PDC²²¹⁸. The gene affected in PDC²²¹⁸ is PaTim54, encoding the homologue of the budding yeast mitochondrial inner membrane import machinery component Tim54p. We show that PaTim54 is essential in P. anserina and that the phenotypes displayed by the PDC²²¹⁸ mutant, renamed PaTim54²²¹⁸, are the consequence of a drastic reduction in the expression of PaTim54. Among these pleiotropic phenotypes, PDC²²¹⁸-PaTim54²²¹⁸- displays increased lifespan, a phenotype in line with the observed mitochondrial defects in the mutant.

Annotation and analysis of the mitochondrial genome of Coniothyrium glycines, causal agent of red leaf blotch of soybean, reveals an abundance of homing endonucleases

Coniothyrium glycines, the causal agent of soybean red leaf blotch, is a USDA APHIS-listed Plant Pathogen Select Agent and potential threat to US agriculture. Sequencing of the C. glycines mt genome revealed a circular 98,533-bp molecule with a mean GC content of 29.01%. It contains twelve of the mitochondrial genes typically involved in oxidative phosphorylation (atp6, cob, cox1-3, nad1-6, and nad4L), one for a ribosomal protein (rps3), four for hypothetical proteins, one for each of the small and large subunit ribosomal RNAs (rns and rnl) and a set of 30 tRNAs. Genes were encoded on both DNA strands with cox1 and cox2 occurring as adjacent genes having no intergenic spacers. Likewise, nad2 and nad3 are adjacent with no intergenic spacers and nad5 is immediately followed by nad4L with an overlap of one base. Thirty-two introns, comprising 54.1% of the total mt genome, were identified within eight protein-coding genes and the rnl. Eighteen of the introns contained putative intronic ORFs with either LAGLIDADG or GIY-YIG homing endonuclease motifs, and an additional eleven introns showed evidence of truncated or degenerate endonuclease motifs. One intron possessed a degenerate N-acetyl-transferase domain. C. glycines shares some conservation of gene order with other members of the Pleosporales, most notably nad6-rnl-atp6 and associated conserved tRNA clusters. Phylogenetic analysis of the twelve shared protein coding genes agrees with commonly accepted fungal taxonomy. C. glycines represents the second largest mt genome from a member of the Pleosporales sequenced to date. This research provides the first genomic information on C. glycines, which may provide targets for rapid diagnostic assays and population studies.

Genomic analyses reveal low mitochondrial and high nuclear diversity in the cyclosporin-producing fungus Tolypocladium inflatum

Mitochondrial DNA is generally regarded to evolve faster than nuclear DNA in animals, whereas if this is also true in fungi remains unclear. Herein, we annotate the first complete mitochondrial genome (mitogenome) of the cyclosporin-producing fungus Tolypocladium inflatum and report the genome-wide sequence variations among five isolates originating from distantly separated localities. We found that T. inflatum has among the most compact of fungal mitogenomes; its 25 kb DNA molecule encodes all standard fungal mitochondrial genes and harbors only one intron. Transcriptional analyses validated the expression of most conserved genes. We found several uncommon repetitive elements and evidence of gene transfer from the mitochondrion to the nucleus. Phylogenetic analyses confirmed the placement of T. inflatum in the fungal order Hypocreales although there was uncertainty on its family-level affiliation. Comparative genomic analyses among the five isolates identified an overall lower level of intraspecific variation in mitogenomes than in nuclear genomes; however, both the nuclear and mitochondrial genomes revealed similar isolate relationships, not correlating with geographic sources of these isolates. Our study shed new insights into the evolution of the medicinally important ascomycete T. inflatum.

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

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.

Characterization of the aodA, dnmA, mnSOD and pimA genes in Aspergillus nidulans

Mitochondria play key roles in cellular energy generation and lifespan of most eukaryotes. To understand the functions of four nuclear-encoded genes predicted to be related to the maintenance of mitochondrial morphology and function in Aspergillus nidulans, systematic characterization was carried out. The deletion and overexpression mutants of aodA, dnmA, mnSOD and pimA encoding alternative oxidase, dynamin related protein, manganese superoxide dismutase and Lon protease, respectively, were generated and examined for their growth, stress tolerances, respiration, autolysis, cell death, sterigmatocystin production, hyphal morphology and size, and mitochondrial superoxide production as well as development. Overall, genetic manipulation of these genes had less effect on cellular physiology and ageing in A. nidulans than that of their homologs in another fungus Podospora anserina with a well-characterized senescence. The observed interspecial phenotypic differences can be explained by the dissimilar intrinsic stabilities of the mitochondrial genomes in A. nidulans and P. anserina. Furthermore, the marginally altered phenotypes observed in A. nidulans mutants indicate the presence of effective compensatory mechanisms for the complex networks of mitochondrial defense and quality control. Importantly, these findings can be useful for developing novel platforms for heterologous protein production, or on new biocontrol and bioremediation technologies based on Aspergillus species.

Complete mitochondrial genome of the medicinal fungus Ophiocordyceps sinensis

As part of a genome sequencing project for Ophiocordyceps sinensis, strain 1229, a complete mitochondrial (mt) genome was assembled as a single circular dsDNA of 157,510 bp, one of the largest reported for fungi. Conserved genes including the large and small rRNA subunits, 27 tRNA and 15 protein-coding genes, were identified. In addition, 58 non-conserved open reading frames (ncORFs) in the intergenic and intronic regions were also identified. Transcription analyses using RNA-Seq validated the expression of most conserved genes and ncORFs. Fifty-two introns (groups I and II) were found within conserved genes, accounting for 68.5% of the genome. Thirty-two homing endonucleases (HEs) with motif patterns LAGLIDADG (21) and GIY-YIG (11) were identified in group I introns. The ncORFs found in group II introns mostly encoded reverse transcriptases (RTs). As in other hypocrealean fungi, gene contents and order were found to be conserved in the mt genome of O. sinensis, but the genome size was enlarged by longer intergenic regions and numerous introns. Intergenic and intronic regions were composed of abundant repetitive sequences usually associated with mobile elements. It is likely that intronic ncORFs, which encode RTs and HEs, may have contributed to the enlarged mt genome of O. sinensis.

The past, present and future of mitochondrial genomics: have we sequenced enough mtDNAs?

The year 2014 saw more than a thousand new mitochondrial genome sequences deposited in GenBank—an almost 15% increase from the previous year. Hundreds of peer-reviewed articles accompanied these genomes, making mitochondrial DNAs (mtDNAs) the most sequenced and reported type of eukaryotic chromosome. These mtDNA data have advanced a wide range of scientific fields, from forensics to anthropology to medicine to molecular evolution. But for many biological lineages, mtDNAs are so well sampled that newly published genomes are arguably no longer contributing significantly to the progression of science, and in some cases they are tying up valuable resources, particularly journal editors and referees. Is it time to acknowledge that as a research community we have published enough mitochondrial genome papers? Here, I address this question, exploring the history, milestones and impacts of mitochondrial genomics, the benefits and drawbacks of continuing to publish mtDNAs at a high rate and what the future may hold for such an important and popular genetic marker. I highlight groups for which mtDNAs are still poorly sampled, thus meriting further investigation, and recommend that more energy be spent characterizing aspects of mitochondrial genomes apart from the DNA sequence, such as their chromosomal and transcriptional architectures. Ultimately, one should be mindful before writing a mitochondrial genome paper. Consider perhaps sending the sequence directly to GenBank instead, and be sure to annotate it correctly before submission.

Evolutionary defined role of the mitochondrial DNA in fertility, disease and ageing

BACKGROUND

The endosymbiosis of an alpha-proteobacterium and a eubacterium a billion years ago paved the way for multicellularity and enabled eukaryotes to flourish. The selective advantage for the host was the acquired ability to generate large amounts of intracellular hydrogen-dependent adenosine triphosphate. The price was increased reactive oxygen species (ROS) inside the eukaryotic cell, causing high mutation rates of the mitochondrial DNA (mtDNA). According to the Muller's ratchet theory, this accumulation of mutations in asexually transmitted mtDNA would ultimately lead to reduced reproductive fitness and eventually extinction. However, mitochondria have persisted over the course of evolution, initially due to a rapid, extreme evolutionary reduction of the mtDNA content. After the phylogenetic divergence of eukaryotes into animals, fungi and plants, differences in evolution of the mtDNA occurred with different adaptations for coping with the mutation burden within these clades. As a result, mitochondrial evolutionary mechanisms have had a profound effect on human adaptation, fertility, healthy reproduction, mtDNA disease manifestation and transmission and ageing. An understanding of these mechanisms might elucidate novel approaches for treatment and prevention of mtDNA disease.

METHODS

The scientific literature was investigated to determine how mtDNA evolved in animals, plants and fungi. Furthermore, the different mechanisms of mtDNA inheritance and of balancing Muller's ratchet in these species were summarized together with the consequences of these mechanisms for human health and reproduction.

RESULTS

Animal, plant and fungal mtDNA have evolved differently. Animals have compact genomes, little recombination, a stable number of genes and a high mtDNA copy number, whereas plants have larger genomes with variable gene counts, a low mtDNA copy number and many recombination events. Fungal mtDNA is somewhere in between. In plants, the mtDNA mutation rate is kept low by effective ROS defence and efficient recombination-mediated mtDNA repair. In animal mtDNA, these mechanisms are not or less well-developed and the detrimental mutagenesis events are controlled by a high mtDNA copy number in combination with a genetic bottleneck and purifying selection during transmission. The mtDNA mutation rates in animals are higher than in plants, which allow mobile animals to adapt more rapidly to various environmental conditions in terms of energy production, whereas static plants do not have this need. Although at the level of the species, these mechanisms have been extremely successful, they can have adverse effects for the individual, resulting, in humans, in severe or unpredictably segregating mtDNA diseases, as well as fertility problems and unhealthy ageing.

CONCLUSIONS

Understanding the forces and processes that underlie mtDNA evolution among different species increases our knowledge on the detrimental consequences that individuals can have from these evolutionary end-points. Alternative outcomes in animals, fungi and plants will lead to a better understanding of the inheritance of mtDNA disorders and mtDNA-related fertility problems. These will allow the development of options to ameliorate, cure and/or prevent mtDNA diseases and mtDNA-related fertility problems.

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.

Comparative genomics of the major fungal agents of human and animal Sporotrichosis: Sporothrix schenckii and Sporothrix brasiliensis

BACKGROUND

The fungal genus Sporothrix includes at least four human pathogenic species. One of these species, S. brasiliensis, is the causal agent of a major ongoing zoonotic outbreak of sporotrichosis in Brazil. Elsewhere, sapronoses are caused by S. schenckii and S. globosa. The major aims on this comparative genomic study are: 1) to explore the presence of virulence factors in S. schenckii and S. brasiliensis; 2) to compare S. brasiliensis, which is cat-transmitted and infects both humans and cats with S. schenckii, mainly a human pathogen; 3) to compare these two species to other human pathogens (Onygenales) with similar thermo-dimorphic behavior and to other plant-associated Sordariomycetes.

RESULTS

The genomes of S. schenckii and S. brasiliensis were pyrosequenced to 17x and 20x coverage comprising a total of 32.3 Mb and 33.2 Mb, respectively. Pair-wise genome alignments revealed that the two species are highly syntenic showing 97.5% average sequence identity. Phylogenomic analysis reveals that both species diverged about 3.8-4.9 MYA suggesting a recent event of speciation. Transposable elements comprise respectively 0.34% and 0.62% of the S. schenckii and S. brasiliensis genomes and expansions of Gypsy-like elements was observed reflecting the accumulation of repetitive elements in the S. brasiliensis genome. Mitochondrial genomic comparisons showed the presence of group-I intron encoding homing endonucleases (HE's) exclusively in S. brasiliensis. Analysis of protein family expansions and contractions in the Sporothrix lineage revealed expansion of LysM domain-containing proteins, small GTPases, PKS type1 and leucin-rich proteins. In contrast, a lack of polysaccharide lyase genes that are associated with decay of plants was observed when compared to other Sordariomycetes and dimorphic fungal pathogens, suggesting evolutionary adaptations from a plant pathogenic or saprobic to an animal pathogenic life style.

CONCLUSIONS

Comparative genomic data suggest a unique ecological shift in the Sporothrix lineage from plant-association to mammalian parasitism, which contributes to the understanding of how environmental interactions may shape fungal virulence. . Moreover, the striking differences found in comparison with other dimorphic fungi revealed that dimorphism in these close relatives of plant-associated Sordariomycetes is a case of convergent evolution, stressing the importance of this morphogenetic change in fungal pathogenesis.

The 203 kbp mitochondrial genome of the phytopathogenic fungus Sclerotinia borealis reveals multiple invasions of introns and genomic duplications

Here we report the complete sequence of the mitochondrial (mt) genome of the necrotrophic phytopathogenic fungus Sclerotinia borealis, a member of the order Helotiales of Ascomycetes. The 203,051 bp long mtDNA of S. borealis represents one of the largest sequenced fungal mt genomes. The large size is mostly determined by the presence of mobile genetic elements, which include 61 introns. Introns contain a total of 125,394 bp, are scattered throughout the genome, and are found in 12 protein-coding genes and in the ribosomal RNA genes. Most introns contain complete or truncated ORFs that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Integrations of mobile elements are also evidenced by the presence of two regions similar to fragments of inverton-like plasmids. Although duplications of some short genome regions, resulting in the appearance of truncated extra copies of genes, did occur, we found no evidences of extensive accumulation of repeat sequences accounting for mitochondrial genome size expansion in some other fungi. Comparisons of mtDNA of S. borealis with other members of the order Helotiales reveal considerable gene order conservation and a dynamic pattern of intron acquisition and loss during evolution. Our data are consistent with the hypothesis that horizontal DNA transfer has played a significant role in the evolution and size expansion of the S. borealis mt genome.

Overexpression of Pa_1_10620 encoding a mitochondrial Podospora anserina protein with homology to superoxide dismutases and ribosomal proteins leads to lifespan extension

In biological systems, reactive oxygen species (ROS) represent 'double edged swords': as signaling molecules they are essential for proper development, as reactive agents they cause molecular damage and adverse effects like degeneration and aging. A well-coordinated control of ROS is therefore of key importance. Superoxide dismutases (SODs) are enzymes active in the detoxification of superoxide. The number of isoforms of these proteins varies among species. Here we report the characterization of the putative protein encoded by Pa_1_10620 that has been previously annotated to code for a mitochondrial ribosomal protein but shares also sequence domains with SODs. We report that the gene is transcribed in P. anserina cultures of all ages and that the encoded protein localizes to mitochondria. In strains overexpressing Pa_1_10620 in a genetic background in which PaSod3, the mitochondrial MnSOD of P. anserina, is deleted, no SOD activity could be identified in isolated mitochondria. However, overexpression of the gene leads to lifespan extension suggesting a pro-survival function of the protein in P. anserina.

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.

Establishment and interpretation of the genome sequence of the phytopathogenic fungus Rhizoctonia solani AG1-IB isolate 7/3/14

Anastomosis group AG1-IB isolates of the anamorphic basidiomycetous fungus Rhizoctonia solani Kühn affect various agricultural and horticultural important crops including bean, rice, soybean, figs, hortensia, cabbage and lettuce. To gain insights into the genome structure and content, the first draft genome sequence of R. solani AG1-IB isolate 7/3/14 was established. Four complete runs on the Genome Sequencer (GS) FLX platform (Roche Applied Science) yielding approx. a 25-fold coverage of the R. solani genome were accomplished. Assembly of the sequence reads by means of the gsAssembler software version 2.6 applying the heterozygotic mode resulted in numerous contigs and scaffolds and a predicted size of 87.1 Mb for the diploid status of the genome. 'Contig-length vs. read-count' analysis revealed that the assembled contigs can be classified into five different groups. Detailed BLAST-analysis revealed that most contigs of group II feature high-scoring matches to other contigs of the same group suggesting that distinguishable allelic variants exist for many genes. Due to the supposed diploid and heterokaryotic nature of R. solani AG1-IB 7/3/14, this result has been anticipated. However, the heterokaryotic character of the isolate is not really supported by sequencing data obtained for the isolate R. solani AG1-IB 7/3/14. Coverage of group III contigs is twice as high as for group II contigs which can also be explained by the diploid status of the genome and indistinguishable alleles on homologous chromosomes. Assembly of sequence data led to the identification of the rRNA unit (group V contigs) and the mitochondrial (mt) genome (group IV contigs) which is a circular molecule of 162,751 bp in size featuring a GC-content of 36.4%. The R. solani 7/3/14 mt-genome is one of the largest fungal mitochondrial genomes known to date. Its large size essentially is due to the presence of numerous non-conserved hypothetical ORFs and introns. Gene prediction for the R. solani AG1-IB 7/3/14 genome was conducted by the Augustus Gene Prediction Software for Eukaryotes (version 2.6.) applying the parameter set for the fungus Coprinopsis cinerea okayama 7#130. Gene prediction and annotation provided first insights into the R. solani AG1-IB 7/3/14 gene structure and content. In total, 12,422 genes were predicted. The average number of exons per gene is five. Exons have a mean length of 214 bp, whereas introns on average are 66 bp in length. Annotation of the genome revealed that 4169 of 12,422 genes could be assigned to KOG functional categories.

Sequencing of mitochondrial genomes of nine Aspergillus and Penicillium species identifies mobile introns and accessory genes as main sources of genome size variability

BACKGROUND

The genera Aspergillus and Penicillium include some of the most beneficial as well as the most harmful fungal species such as the penicillin-producer Penicillium chrysogenum and the human pathogen Aspergillus fumigatus, respectively. Their mitochondrial genomic sequences may hold vital clues into the mechanisms of their evolution, population genetics, and biology, yet only a handful of these genomes have been fully sequenced and annotated.

RESULTS

Here we report the complete sequence and annotation of the mitochondrial genomes of six Aspergillus and three Penicillium species: A. fumigatus, A. clavatus, A. oryzae, A. flavus, Neosartorya fischeri (A. fischerianus), A. terreus, P. chrysogenum, P. marneffei, and Talaromyces stipitatus (P. stipitatum). The accompanying comparative analysis of these and related publicly available mitochondrial genomes reveals wide variation in size (25-36 Kb) among these closely related fungi. The sources of genome expansion include group I introns and accessory genes encoding putative homing endonucleases, DNA and RNA polymerases (presumed to be of plasmid origin) and hypothetical proteins. The two smallest sequenced genomes (A. terreus and P. chrysogenum) do not contain introns in protein-coding genes, whereas the largest genome (T. stipitatus), contains a total of eleven introns. All of the sequenced genomes have a group I intron in the large ribosomal subunit RNA gene, suggesting that this intron is fixed in these species. Subsequent analysis of several A. fumigatus strains showed low intraspecies variation. This study also includes a phylogenetic analysis based on 14 concatenated core mitochondrial proteins. The phylogenetic tree has a different topology from published multilocus trees, highlighting the challenges still facing the Aspergillus systematics.

CONCLUSIONS

The study expands the genomic resources available to fungal biologists by providing mitochondrial genomes with consistent annotations for future genetic, evolutionary and population studies. Despite the conservation of the core genes, the mitochondrial genomes of Aspergillus and Penicillium species examined here exhibit significant amount of interspecies variation. Most of this variation can be attributed to accessory genes and mobile introns, presumably acquired by horizontal gene transfer of mitochondrial plasmids and intron homing.

Phylogenetic analysis of the complete mitochondrial genome of Madurella mycetomatis confirms its taxonomic position within the order Sordariales

BACKGROUND

Madurella mycetomatis is the most common cause of human eumycetoma. The genus Madurella has been characterized by overall sterility on mycological media. Due to this sterility and the absence of other reliable morphological and ultrastructural characters, the taxonomic classification of Madurella has long been a challenge. Mitochondria are of monophyletic origin and mitochondrial genomes have been proven to be useful in phylogenetic analyses.

RESULTS

The first complete mitochondrial DNA genome of a mycetoma-causative agent was sequenced using 454 sequencing. The mitochondrial genome of M. mycetomatis is a circular DNA molecule with a size of 45,590 bp, encoding for the small and the large subunit rRNAs, 27 tRNAs, 11 genes encoding subunits of respiratory chain complexes, 2 ATP synthase subunits, 5 hypothetical proteins, 6 intronic proteins including the ribosomal protein rps3. In phylogenetic analyses using amino acid sequences of the proteins involved in respiratory chain complexes and the 2 ATP synthases it appeared that M. mycetomatis clustered together with members of the order Sordariales and that it was most closely related to Chaetomium thermophilum. Analyses of the gene order showed that within the order Sordariales a similar gene order is found. Furthermore also the tRNA order seemed mostly conserved.

CONCLUSION

Phylogenetic analyses of fungal mitochondrial genomes confirmed that M. mycetomatis belongs to the order of Sordariales and that it was most closely related to Chaetomium thermophilum, with which it also shared a comparable gene and tRNA order.

Mitochondrial genomes from the lichenized fungi Peltigera membranacea and Peltigera malacea: features and phylogeny

Mitochondrial genomes from the fungal partners of two terricolous foliose lichen symbioses, Peltigera membranacea and Peltigera malacea, have been determined using metagenomic approaches, including RNA-seq. The roughly 63 kb genomes show all the major features found in other Pezizomycotina, such as unidirectional transcription, 14 conserved protein genes, genes for the two subunit rRNAs and for a set of 26 tRNAs used in translating the 62 amino acid codons. In one of the tRNAs a CAU anticodon is proposed to be modified, via the action of the nuclear-encoded enzyme, tRNA Ile lysidine synthase, so that it recognizes the codon AUA (Ile) instead of AUG (Met). The overall arrangements and sequences of the two circular genomes are similar, the major difference being the inversion and deterioration of a gene encoding a type B DNA polymerase. Both genomes encode the RNA component of RNAse P, a feature seldom found in ascomycetes. The difference in genome size from the minimal ascomycete mitochondrial genomes is largely due to 17 and 20 group I introns, respectively, most associated with homing endonucleases and all found within protein-coding genes and the gene encoding the large subunit rRNA. One new intron insertion point was found, and an unusually small exon of seven nucleotides (nt) was identified and verified by RNA sequencing. Comparative analysis of mitochondrion-encoded proteins places the Peltigera spp., representatives of the class Lecanoromycetes, close to Leotiomycetes, Dothidiomycetes, and Sordariomycetes, in contrast to phylogenies found using nuclear genes.

Mitochondrial genome evolution in species belonging to the Phialocephala fortinii s.l. - Acephala applanata species complex

Background

Mitochondrial (mt) markers are successfully applied in evolutionary biology and systematics because mt genomes often evolve faster than the nuclear genomes. In addition, they allow robust phylogenetic analysis based on conserved proteins of the oxidative phosphorylation system. In the present study we sequenced and annotated the complete mt genome of P. subalpina, a member of the Phialocephala fortinii s.l. – Acephala applanata species complex (PAC). PAC belongs to the Helotiales, which is one of the most diverse groups of ascomycetes including more than 2,000 species. The gene order was compared to deduce the mt genome evolution in the Pezizomycotina. Genetic variation in coding and intergenic regions of the mtDNA was studied for PAC to assess the usefulness of mt DNA for species diagnosis.

Results

The mt genome of P. subalpina is 43,742 bp long and codes for 14 mt genes associated with the oxidative phosphorylation. In addition, a GIY-YIG endonuclease, the ribosomal protein S3 (Rps3) and a putative N-acetyl-transferase were recognized. A complete set of tRNA genes as well as the large and small rRNA genes but no introns were found. All protein-coding genes were confirmed by EST sequences. The gene order in P. subalpina deviated from the gene order in Sclerotinia sclerotiorum, the only other helotialean species with a fully sequenced and annotated mt genome. Gene order analysis within Pezizomycotina suggests that the evolution of gene orders is mostly driven by transpositions. Furthermore, sequence diversity in coding and non-coding mtDNA regions in seven additional PAC species was pronounced and allowed for unequivocal species diagnosis in PAC.

Conclusions

The combination of non-interrupted ORFs and EST sequences resulted in a high quality annotation of the mt genome of P. subalpina, which can be used as a reference for the annotation of other mt genomes in the Helotiales. In addition, our analyses show that mtDNA loci will be the marker of choice for future analysis of PAC communities.

Complete mitochondrial genome of compactin-producing fungus Penicillium solitum and comparative analysis of Trichocomaceae mitochondrial genomes

We determined the complete mitochondrial genome sequence of the compactin-producing fungus Penicillium solitum strain 20-01. The 28 601-base pair circular-mapping DNA molecule encodes a characteristic set of mitochondrial proteins and RNA genes and is intron-free. All 46 protein- and RNA-encoding genes are located on one strand and apparently transcribed in one direction. Comparative analysis of this mtDNA and previously sequenced but unannotated mitochondrial genomes of several medically and industrially important species of the Aspergillus/Penicillium group revealed their extensive similarity in terms of size, gene content and sequence, which is also reflected in the almost perfect conservation of mitochondrial gene order in Penicillium and Aspergillus. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed the monophyletic origin of Eurotiomycetes.

Approaches to Fungal Genome Annotation

Fungal genome annotation is the starting point for analysis of genome content. This generally involves the application of diverse methods to identify features on a genome assembly such as protein-coding and non-coding genes, repeats and transposable elements, and pseudogenes. Here we describe tools and methods leveraged for eukaryotic genome annotation with a focus on the annotation of fungal nuclear and mitochondrial genomes. We highlight the application of the latest technologies and tools to improve the quality of predicted gene sets. The Broad Institute eukaryotic genome annotation pipeline is described as one example of how such methods and tools are integrated into a sequencing center's production genome annotation environment.

Genome-wide gene expression profiling of fertilization competent mycelium in opposite mating types in the heterothallic fungus Podospora anserina

Background

Mating-type loci in yeasts and ascomycotan filamentous fungi (Pezizomycotina) encode master transcriptional factors that play a critical role in sexual development. Genome-wide analyses of mating-type-specification circuits and mating-type target genes are available in Saccharomyces cerevisiae and Schizosaccharomyces pombe; however, no such analyses have been performed in heterothallic (self-incompatible) Pezizomycotina. The heterothallic fungus Podospora anserina serves as a model for understanding the basic features of mating-type control. Its mat+ and mat− mating types are determined by dissimilar allelic sequences. The mat− sequence contains three genes, designated FMR1, SMR1 and SMR2, while the mat+ sequence contains one gene, FPR1. FMR1 and FPR1 are the major regulators of fertilization, and this study presents a genome-wide view of their target genes and analyzes their target gene regulation.

Methodology/Principal Findings

The transcriptomic profiles of the mat+ and mat− strains revealed 157 differentially transcribed genes, and transcriptomic analysis of fmr1⁻ and fpr1⁻ mutant strains was used to determine the regulatory actions exerted by FMR1 and FPR1 on these differentially transcribed genes. All possible combinations of transcription repression and/or activation by FMR1 and/or FPR1 were observed. Furthermore, 10 additional mating-type target genes were identified that were up- or down-regulated to the same level in mat+ and mat− strains. Of the 167 genes identified, 32 genes were selected for deletion, which resulted in the identification of two genes essential for the sexual cycle. Interspecies comparisons of mating-type target genes revealed significant numbers of orthologous pairs, although transcriptional profiles were not conserved between species.

Conclusions/Significance

This study represents the first comprehensive genome-wide analysis of mating-type direct and indirect target genes in a heterothallic filamentous fungus. Mating-type transcription factors have many more target genes than are found in yeasts and exert a much greater diversity of regulatory actions on target genes, most of which are not directly related to mating.

The subunit composition of mitochondrial NADH:ubiquinone oxidoreductase (complex I) from Pichia pastoris

Respiratory complex I (NADH:quinone oxidoreductase) is an entry point to the electron transport chain in the mitochondria of many eukaryotes. It is a large, multisubunit enzyme with a hydrophilic domain in the matrix and a hydrophobic domain in the mitochondrial inner membrane. Here we present a comprehensive analysis of the protein composition and post-translational modifications of complex I from Pichia pastoris, using a combination of proteomic and bioinformatic approaches. Forty-one subunits were identified in P. pastoris complex I, comprising the 14 core (conserved) subunits and 27 supernumerary subunits; seven of the core subunits are mitochondrial encoded. Three of the supernumerary subunits (named NUSM, NUTM, and NUUM) have not been observed previously in any species of complex I. However, homologues to all three of them are present in either Yarrowia lipolytica or Pichia angusta complex I. P. pastoris complex I has 39 subunits in common with Y. lipolytica complex I, 37 in common with N. crassa complex I, and 35 in common with the bovine enzyme. The mitochondrial encoded subunits (translated by the mold mitochondrial genetic code) retain their N-α-formyl methionine residues. At least eight subunits are N-α-acetylated, but the N-terminal modifications of the nuclear encoded subunits are not well-conserved. A combination of two methods of protein separation (SDS-PAGE and HPLC) and three different mass spectrometry techniques (peptide mass fingerprinting, tandem MS and molecular mass measurements) were required to define the protein complement of P. pastoris complex I. This requirement highlights the need for inclusive and comprehensive strategies for the characterization of challenging membrane-bound protein complexes containing both hydrophilic and hydrophobic components.

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.

Mitochondrial recombination increases with age in Podospora anserina

With uniparental inheritance of mitochondria, there seems little reason for homologous recombination in mitochondria, but the machinery for mitochondrial recombination is quite well-conserved in many eukaryote species. In fungi and yeasts heteroplasmons may be formed when strains fuse and transfer of organelles takes place, making it possible to study mitochondrial recombination when introduced mitochondria contain different markers. A survey of wild-type isolates from a local population of the filamentous fungus Podospora anserina for the presence of seven optional mitochondrial introns indicated that mitochondrial recombination does take place in nature. Moreover the recombination frequency appeared to be correlated with age: the more rapidly ageing fraction of the population had a significantly lower linkage disequilibrium indicating more recombination. Direct confrontation experiments with heterokaryon incompatible strains with different mitochondrial markers at different (relative) age confirmed that mitochondrial recombination increases with age. We propose that with increasing mitochondrial damage over time, mitochondrial recombination - even within a homoplasmic population of mitochondria - is a mechanism that may restore mitochondrial function.

Heteroplasmy and stoichiometric complexity of plant mitochondrial genomes--though this be madness, yet there's method in't

Mitochondrial heteroplasmy is defined as the coexistence of divergent mitochondrial genotypes in a cell. The ratio of the alternative genomes may be variable, but in plants, the usually prevalent main genome is accompanied by sublimons--substoichiometric mitochondrial DNA (mtDNA) molecules. Plant mitochondrial heteroplasmy was originally viewed as being associated with pathological mutations or was found in non-natural plant populations. Currently, it is considered to be a common situation in plants. Recent years have changed the previous view on the role of homologous recombination, small-scale mutations, and paternal leakage of mtDNA in the generation of heteroplasmy. Newly developed sensitive techniques have allowed the precise estimation of mtDNA stoichiometry. Mechanisms of maintenance and transmission of heteroplasmic genomes, including DNA recombination and replication, as well as mitochondrial fusion and fission, have been studied. This review describes the high level of plant mitochondrial genome complication--the 'madness' resulting from the heteroplasmic state and explains the method hidden in this madness. Heteroplasmy is described as the evolutionary strategy of uniparentally inherited plant mitochondrial genomes which do not undergo sexual recombination. In order to compensate for this deficiency, alternative types of mtDNA are substoichiometrically accumulated as a reservoir of genetic variability and may undergo accelerated evolution. Occasionally, sublimons are selected and amplified in the process called substoichiometric shifting, to take over the role of the main genome. Alternative mitochondrial genomes may recombine, yielding new mtDNA variants, or segregate during plant growth resulting in plants with mosaic phenotypes. Two opposite roles of mitochondrial heteroplasmy with respect to acceleration or counteracting of mutation accumulation are also discussed. Finally, nuclear control of heteroplasmy and substoichiometric shifting is described.

Insights into the evolution of mitochondrial genome size from complete sequences of Citrullus lanatus and Cucurbita pepo (Cucurbitaceae)

The mitochondrial genomes of seed plants are unusually large and vary in size by at least an order of magnitude. Much of this variation occurs within a single family, the Cucurbitaceae, whose genomes range from an estimated 390 to 2,900 kb in size. We sequenced the mitochondrial genomes of Citrullus lanatus (watermelon: 379,236 nt) and Cucurbita pepo (zucchini: 982,833 nt)--the two smallest characterized cucurbit mitochondrial genomes--and determined their RNA editing content. The relatively compact Citrullus mitochondrial genome actually contains more and longer genes and introns, longer segmental duplications, and more discernibly nuclear-derived DNA. The large size of the Cucurbita mitochondrial genome reflects the accumulation of unprecedented amounts of both chloroplast sequences (>113 kb) and short repeated sequences (>370 kb). A low mutation rate has been hypothesized to underlie increases in both genome size and RNA editing frequency in plant mitochondria. However, despite its much larger genome, Cucurbita has a significantly higher synonymous substitution rate (and presumably mutation rate) than Citrullus but comparable levels of RNA editing. The evolution of mutation rate, genome size, and RNA editing are apparently decoupled in Cucurbitaceae, reflecting either simple stochastic variation or governance by different factors.

Suppression of mitochondrial DNA instability of autosomal dominant forms of progressive external ophthalmoplegia-associated ANT1 mutations in Podospora anserina

Maintenance and expression of mitochondrial DNA (mtDNA) are essential for the cell and the organism. In humans, several mutations in the adenine nucleotide translocase gene ANT1 are associated with multiple mtDNA deletions and autosomal dominant forms of progressive external ophthalmoplegia (adPEO). The mechanisms underlying the mtDNA instability are still obscure. A current hypothesis proposes that these pathogenic mutations primarily uncouple the mitochondrial inner membrane, which secondarily causes mtDNA instability. Here we show that the three adPEO-associated mutations equivalent to A114P, L98P, and V289M introduced into the Podospora anserina ANT1 ortholog dominantly cause severe growth defects, decreased reactive oxygen species production (ROS), decreased mitochondrial inner membrane potential (Deltapsi), and accumulation of large-scale mtDNA deletions leading to premature death. Interestingly, we show that, at least for the adPEO-type M106P and A121P mutant alleles, the associated mtDNA instability cannot be attributed only to a reduced membrane potential or to an increased ROS level since it can be suppressed without restoration of the Deltapsi or modification of the ROS production. Suppression of mtDNA instability due to the M106P and A121P mutations was obtained by an allele of the rmp1 gene involved in nucleo-mitochondrial cross- talk and also by an allele of the AS1 gene encoding a cytosolic ribosomal protein. In contrast, the mtDNA instability caused by the S296M mutation was not suppressed by these alleles.