The complete mitochondrial genome of Fusarium oxysporum: Insights into fungal mitochondrial evolution

An Alternative Self-Splicing Intron Lifecycle Revealed by Dynamic Intron Turnover in Epichloë Endophyte Mitochondrial Genomes

Self-splicing group I and II introns are selfish genetic elements that are widely yet patchily distributed across the tree of life. Their selfish behavior comes from super-Mendelian inheritance behaviors, collectively called "homing", which allow them to rapidly spread within populations to the specific genomic sites they home into. Observations of self-splicing intron evolutionary dynamics have led to the formulation of an intron "lifecycle" model where, once fixed in a population, the introns lose selection for homing and undergo an extensive period of degradation until their eventual loss. Here, we find that self-splicing introns are common in the mitochondrial genomes of Epichloë species, endophytic fungi that live in symbioses with grasses. However, these introns show substantial intron presence-absence polymorphism, with our analyses suggesting that these result from a combination of vertical intron inheritance coupled with multiple invasion and loss events over the course of Epichloë evolution. Surprisingly, we find little evidence for the extensive intron degradation expected under the existing intron lifecycle model. Instead, these introns in Epichloë appear to be lost soon after fixation, suggesting that Epichloë self-splicing introns have a different lifecycle. However, rapid intron loss alone cannot explain our results, indicating that additional factors, such as the evolution of homing suppressors, also contribute to Epichloë self-splicing intron dynamics. This work shows that self-splicing introns have more diverse evolutionary dynamics than previously appreciated.

Characterisation and comparative analysis of mitochondrial genomes of false, yellow, black and blushing morels provide insights on their structure and evolution

Morchella species have considerable significance in terrestrial ecosystems, exhibiting a range of ecological lifestyles along the saprotrophism-to-symbiosis continuum. However, the mitochondrial genomes of these ascomycetous fungi have not been thoroughly studied, thereby impeding a comprehensive understanding of their genetic makeup and ecological role. In this study, we analysed the mitogenomes of 30 Morchellaceae species, including yellow, black, blushing and false morels. These mitogenomes are either circular or linear DNA molecules with lengths ranging from 217 to 565 kbp and GC content ranging from 38% to 48%. Fifteen core protein-coding genes, 28-37 tRNA genes and 3-8 rRNA genes were identified in these Morchellaceae mitogenomes. The gene order demonstrated a high level of conservation, with the cox1 gene consistently positioned adjacent to the rnS gene and cob gene flanked by apt genes. Some exceptions were observed, such as the rearrangement of atp6 and rps3 in Morchellaimportuna and the reversed order of atp6 and atp8 in certain morel mitogenomes. However, the arrangement of the tRNA genes remains conserved. We additionally investigated the distribution and phylogeny of homing endonuclease genes (HEGs) of the LAGLIDADG (LAGs) and GIY-YIG (GIYs) families. A total of 925 LAG and GIY sequences were detected, with individual species containing 19-48HEGs. These HEGs were primarily located in the cox1, cob, cox2 and nad5 introns and their presence and distribution displayed significant diversity amongst morel species. These elements significantly contribute to shaping their mitogenome diversity. Overall, this study provides novel insights into the phylogeny and evolution of the Morchellaceae.

Co-evolution of large inverted repeats and G-quadruplex DNA in fungal mitochondria may facilitate mitogenome stability: the case of Malassezia

Mitogenomes are essential due to their contribution to cell respiration. Recently they have also been implicated in fungal pathogenicity mechanisms. Members of the basidiomycetous yeast genus Malassezia are an important fungal component of the human skin microbiome, linked to various skin diseases, bloodstream infections, and they are increasingly implicated in gut diseases and certain cancers. In this study, the comparative analysis of Malassezia mitogenomes contributed to phylogenetic tree construction for all species. The mitogenomes presented significant size and gene order diversity which correlates to their phylogeny. Most importantly, they showed the inclusion of large inverted repeats (LIRs) and G-quadruplex (G4) DNA elements, rendering Malassezia mitogenomes a valuable test case for elucidating the evolutionary mechanisms responsible for this genome diversity. Both LIRs and G4s coexist and convergently evolved to provide genome stability through recombination. This mechanism is common in chloroplasts but, hitherto, rarely found in mitogenomes.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

Saccharomycotina yeasts belong to diverse clades within the kingdom of fungi and are important to human everyday life. This work investigates the evolutionary relationships among these yeasts from a mitochondrial (mt) genomic perspective. A comparative study of 155 yeast mt genomes representing all major phylogenetic lineages of Saccharomycotina was performed, including genome size and content variability, intron and intergenic regions’ diversity, genetic code alterations, and syntenic variation. Findings from this study suggest that mt genome size diversity is the result of a ceaseless random process, mainly based on genetic recombination and intron mobility. Gene order analysis revealed conserved syntenic units and many occurring rearrangements, which can be correlated with major evolutionary events as shown by the phylogenetic analysis of the concatenated mt protein matrix. For the first time, molecular dating indicated a slower mt genome divergence rate in the early stages of yeast evolution, in contrast with a faster rate in the late evolutionary stages, compared to their nuclear time divergence. Genetic code reassignments of mt genomes are a perpetual process happening in many different parallel evolutionary steps throughout the evolution of Saccharomycotina. Overall, this work shows that phylogenetic studies based on the mt genome of yeasts highlight major evolutionary events.

Mitochondrial Transcription of Entomopathogenic Fungi Reveals Evolutionary Aspects of Mitogenomes

Entomopathogenic fungi and more specifically genera Beauveria and Metarhizium have been exploited for the biological control of pests. Genome analyses are important to understand better their mode of action and thus, improve their efficacy against their hosts. Until now, the sequences of their mitochondrial genomes were studied, but not at the level of transcription. Except of yeasts and Neurospora crassa, whose mt gene transcription is well described, in all other Ascomycota, i.e., Pezizomycotina, related information is extremely scarce. In this work, mt transcription and key enzymes of this function were studied. RT-PCR experiments and Northern hybridizations reveal the transcriptional map of the mt genomes of B. bassiana and M. brunneum species. The mt genes are transcribed in six main transcripts and undergo post-transcriptional modifications to create single gene transcripts. Promoters were determined in both mt genomes with a comparative in silico analysis, including all known information from other fungal mt genomes. The promoter consensus sequence is 5'-ATAGTTATTAT-3' which is in accordance with the definition of the polycistronic transcripts determined with the experiments described above. Moreover, 5'-RACE experiments in the case of premature polycistronic transcript nad1-nad4-atp8-atp6 revealed the 5' end of the RNA transcript immediately after the in silico determined promoter, as also found in other fungal species. Since several conserved elements were retrieved from these analyses compared to the already known data from yeasts and N. crassa, the phylogenetic analyses of mt RNA polymerase (Rpo41) and its transcriptional factor (Mtf1) were performed in order to define their evolution. As expected, it was found that fungal Rpo41 originate from the respective polymerase of T7/T3 phages, while the ancestor of Mtf1 is of alpha-proteobacterial origin. Therefore, this study presents insights about the fidelity of the mt single-subunit phage-like RNA polymerase during transcription, since the correct identification of mt promoters from Rpo41 requires an ortholog to bacterial sigma factor, i.e., Mtf1. Thus, a previously proposed hypothesis of a phage infected alpha-proteobacterium as the endosymbiotic progenitor of mitochondrion is confirmed in this study and further upgraded by the co-evolution of the bacterial (Mtf1) and viral (Rpo41) originated components in one functional unit.

LAMP Assay for Distinguishing Fusarium oxysporum and Fusarium commune in Lotus (Nelumbo nucifera) Rhizomes

Yields of edible rhizome from cultivation of the perennial hydrophyte lotus (Nelumbo nucifera) can be severely reduced by rhizome rot disease caused by Fusarium species. There is a lack of rapid field-applicable methods for detection of these pathogens on lotus plants displaying symptoms of rhizome rot. Fusarium commune (91%) and Fusarium oxysporum (9%) were identified at different frequencies from lotus samples showing symptoms of rhizome rot. Because these two species can cause different severity of disease and their morphology is similar, molecular diagnostic-based methods to detect these two species were developed. Based on the comparison of the mitochondrial genome of the two species, three specific DNA loci targets were found. The designed primer sets for conventional PCR, quantitative PCR, and loop-mediated isothermal amplification (LAMP) precisely distinguished the above two species when isolated from lotus and other plants. The LAMP detection limits were 10 pg/μl and 1 pg/μl of total DNA for F. commune and F. oxysporum, respectively. We also carried out field-mimicked experiments on lotus seedlings and rhizomes (including inoculated samples and field-diseased samples), and the results indicated that the LAMP primer sets and the supporting portable methods are suitable for rapid diagnosis of the lotus disease in the field. The LAMP-based detection method will aid in the rapid identification of whether F. oxysporum or F. commune is infecting lotus plants with symptoms of rhizome rot and can facilitate efficient pesticide use and prevent disease spread through vegetative propagation of Fusarium-infected lotus rhizomes.

The Coevolution of Fungal Mitochondrial Introns and Their Homing Endonucleases (GIY-YIG and LAGLIDADG)

Fungal mitochondrial (mt) genomes exhibit great diversity in size which is partially attributed to their variable intergenic regions and most importantly to the inclusion of introns within their genes. These introns belong to group I or II, and both of them are self-splicing. The majority of them carry genes encoding homing endonucleases, either LAGLIDADG or GIY-YIG. In this study, it was found that these intronic homing endonucleases genes (HEGs) may originate from mt free-standing open reading frames which can be found nowadays in species belonging to Early Diverging Fungi as “living fossils.” A total of 487 introns carrying HEGs which were located in the publicly available mt genomes of representative species belonging to orders from all fungal phyla was analyzed. Their distribution in the mt genes, their insertion target sequence, and the phylogenetic analyses of the HEGs showed that these introns along with their HEGs form a composite structure in which both selfish elements coevolved. The invasion of the ancestral free-standing HEGs in the introns occurred through a perpetual mechanism, called in this study as “aenaon” hypothesis. It is based on recombination, transpositions, and horizontal gene transfer events throughout evolution. HEGs phylogenetically clustered primarily according to their intron hosts and secondarily to the mt genes carrying the introns and their HEGs. The evolutionary models created revealed an “intron-early” evolution which was enriched by “intron-late” events through many different independent recombinational events which resulted from both vertical and horizontal gene transfers.

Comparative genome analyses suggest a hemibiotrophic lifestyle and virulence differences for the beech bark disease fungal pathogens Neonectria faginata and Neonectria coccinea

Neonectria faginata and Neonectria coccinea are the causal agents of the insect-fungus disease complex known as beech bark disease (BBD), known to cause mortality in beech forest stands in North America and Europe. These fungal species have been the focus of extensive ecological and disease management studies, yet less progress has been made toward generating genomic resources for both micro- and macro-evolutionary studies. Here, we report a 42.1 and 42.7 mb highly contiguous genome assemblies of N. faginata and N. coccinea, respectively, obtained using Illumina technology. These species share similar gene number counts (12,941 and 12,991) and percentages of predicted genes with assigned functional categories (64 and 65%). Approximately 32% of the predicted proteomes of both species are homologous to proteins involved in pathogenicity, yet N. coccinea shows a higher number of predicted mitogen-activated protein kinase genes, virulence determinants possibly contributing to differences in disease severity between N. faginata and N. coccinea. A wide range of genes encoding for carbohydrate-active enzymes capable of degradation of complex plant polysaccharides and a small number of predicted secretory effector proteins, secondary metabolite biosynthesis clusters and cytochrome oxidase P450 genes were also found. This arsenal of enzymes and effectors correlates with, and reflects, the hemibiotrophic lifestyle of these two fungal pathogens. Phylogenomic analysis and timetree estimations indicated that the N. faginata and N. coccinea species divergence may have occurred at ∼4.1 million years ago. Differences were also observed in the annotated mitochondrial genomes as they were found to be 81.7 kb (N. faginata) and 43.2 kb (N. coccinea) in size. The mitochondrial DNA expansion observed in N. faginata is attributed to the invasion of introns into diverse intra- and intergenic locations. These first draft genomes of N. faginata and N. coccinea serve as valuable tools to increase our understanding of basic genetics, evolutionary mechanisms and molecular physiology of these two nectriaceous plant pathogenic species.

A highly conserved mechanism for the detoxification and assimilation of the toxic phytoproduct L-azetidine-2-carboxylic acid in Aspergillus nidulans

Plants produce toxic secondary metabolites as defense mechanisms against phytopathogenic microorganisms and predators. L-azetidine-2-carboxylic acid (AZC), a toxic proline analogue produced by members of the Liliaceae and Agavaciae families, is part of such a mechanism. AZC causes a broad range of toxic, inflammatory and degenerative abnormalities in human and animal cells, while it is known that some microorganisms have evolved specialized strategies for AZC resistance. However, the mechanisms underlying these processes are poorly understood. Here, we identify a widespread mechanism for AZC resistance in fungi. We show that the filamentous ascomycete Aspergillus nidulans is able to not only resist AZC toxicity but also utilize it as a nitrogen source via GABA catabolism and the action of the AzhA hydrolase, a member of a large superfamily of detoxifying enzymes, the haloacid dehalogenase-like hydrolase (HAD) superfamily. This detoxification process is further assisted by the NgnA acetyltransferase, orthologue of Mpr1 of Saccharomyces cerevisiae. We additionally show that heterologous expression of AzhA protein can complement the AZC sensitivity of S. cerevisiae. Furthermore, a detailed phylogenetic analysis of AzhA homologues in Fungi, Archaea and Bacteria is provided. Overall, our results unravel a widespread mechanism for AZC resistance among microorganisms, including important human and plant pathogens.

Revealing the high variability on nonconserved core and mobile elements of Austropuccinia psidii and other rust mitochondrial genomes

Mitochondrial genomes are highly conserved in many fungal groups, and they can help characterize the phylogenetic relationships and evolutionary biology of plant pathogenic fungi. Rust fungi are among the most devastating diseases for economically important crops around the world. Here, we report the complete sequence and annotation of the mitochondrial genome of Austropuccinia psidii (syn. Puccinia psidii), the causal agent of myrtle rust. We performed a phylogenomic analysis including the complete mitochondrial sequences from other rust fungi. The genome composed of 93.299 bp has 73 predicted genes, 33 of which encoded nonconserved proteins (ncORFs), representing almost 45% of all predicted genes. A. psidii mtDNA is one of the largest rust mtDNA sequenced to date, most likely due to the abundance of ncORFs. Among them, 33% were within intronic regions of diverse intron groups. Mobile genetic elements invading intron sequences may have played significant roles in size but not shaping of the rust mitochondrial genome structure. The mtDNAs from rust fungi are highly syntenic. Phylogenetic inferences with 14 concatenated mitochondrial proteins encoded by the core genes placed A. psidii according to phylogenetic analysis based on 18S rDNA. Interestingly, cox1, the gene with the greatest number of introns, provided phylogenies not congruent with the core set. For the first time, we identified the proteins encoded by three A. psidii ncORFs using proteomics analyses. Also, the orf208 encoded a transmembrane protein repressed during in vitro morphogenesis. To the best of our knowledge, we presented the first report of a complete mtDNA sequence of a member of the family Sphaerophragmiacea.

The mitochondrial genome contribution to the phylogeny and identification of Metarhizium species and strains

The genus Metarhizium is composed of entomopathogenic fungal biological control agents (BCAs) used for invertebrate pest control. The phylogenetic relationships of species within this genus are still under scrutiny as several cryptic species can be found. In this work, the mitochondrial (mt) genome of Metarhizium brunneum ARSEF 4556 was fully sequenced and a comparative genome analysis was conducted with 7 other available mt genomes, belonging to 5 Metarhizium species: M. anisopliae, M. brunneum, M. robertsii, M. guizhouense and M. majus. Results showed that Metarhizium demonstrates greater conserved stability than other fungal mt genomes. Furthermore, this analysis located 7 diverse regions in both intergenic domains and gene fragments which were ideal for species/strain discrimination. The sequencing of these regions revealed several SNPs among 38 strains tested, 11 of which were uncharacterized. Single gene phylogenies presented variable results which may be used further for intra-species discrimination. Phylogenetic trees based on the concatenation of mt domains and the nuclear ITS1-5.8S-ITS2 region showed discrimination of the species studied and allowed the identification of uncharacterized strains. These were mostly placed within species M. anisopliae and M. brunneum. Five strains clustered together in a clade related to M. brunneum, suggesting that they comprise a cryptic species.

Diversity of Mobile Genetic Elements in the Mitogenomes of Closely Related Fusarium culmorum and F. graminearum sensu stricto Strains and Its Implication for Diagnostic Purposes

Much of the mitogenome variation observed in fungal lineages seems driven by mobile genetic elements (MGEs), which have invaded their genomes throughout evolution. The variation in the distribution and nucleotide diversity of these elements appears to be the main distinction between different fungal taxa, making them promising candidates for diagnostic purposes. Fungi of the genus Fusarium display a high variation in MGE content, from MGE-poor (Fusarium oxysporum and Fusarium fujikuroi species complex) to MGE-rich mitogenomes found in the important cereal pathogens F. culmorum and F. graminearum sensu stricto. In this study, we investigated the MGE variation in these latter two species by mitogenome analysis of geographically diverse strains. In addition, a smaller set of F. cerealis and F. pseudograminearum strains was included for comparison. Forty-seven introns harboring from 0 to 3 endonucleases (HEGs) were identified in the standard set of mitochondrial protein-coding genes. Most of them belonged to the group I intron family and harbored either LAGLIDADG or GIY-YIG HEGs. Among a total of 53 HEGs, 27 were shared by all fungal strains. Most of the optional HEGs were irregularly distributed among fungal strains/species indicating ancestral mosaicism in MGEs. However, among optional MGEs, one exhibited species-specific conservation in F. culmorum. While in F. graminearum s.s. MGE patterns in cox3 and in the intergenic spacer between cox2 and nad4L may facilitate the identification of this species. Thus, our results demonstrate distinctive traits of mitogenomes for diagnostic purposes of Fusaria.

The Mitochondrial Genome of the Phytopathogenic Fungus Bipolaris sorokiniana and the Utility of Mitochondrial Genome to Infer Phylogeny of Dothideomycetes

A number of species in Bipolaris are important plant pathogens. Due to a limited number of synapomorphic characters, it is difficult to perform species identification and to estimate phylogeny of Bipolaris based solely on morphology. In this study, we sequenced the complete mitochondrial genome of Bipolaris sorokiniana, and presented the detailed annotation of the genome. The B. sorokiniana mitochondrial genome is 137,775 bp long, and contains two ribosomal RNA genes, 12 core protein-coding genes, 38 tRNA genes. In addition, two ribosomal protein genes (rps3 gene and rps5 gene) and the fungal mitochondrial RNase P gene (rnpB) are identified. The large genome size is mostly determined by the presence of numerous intronic and intergenic regions. A total of 28 introns are inserted in eight core protein-coding genes. Together with the published mitochondrial genome sequences, we conducted a preliminary phylogenetic inference of Dothideomycetes under various datasets and substitution models. The monophyly of Capnodiales, Botryosphaeriales and Pleosporales are consistently supported in all analyses. The Venturiaceae forms an independent lineage, with a distant phylogenetic relationship to Pleosporales. At the family level, the Mycosphaerellaceae, Botryosphaeriaceae. Phaeosphaeriaceae, and Pleosporaceae are recognized in the majority of trees.

Population Genomic Analysis Reveals a Highly Conserved Mitochondrial Genome in Fusarium asiaticum

Fusarium asiaticum is one of the pivotal members of the Fusarium graminearum species complex (FGSC) causing Fusarium head blight (FHB) on wheat, barley and rice in large parts of Asia. Besides resulting in yield losses, FHB also causes the accumulation of mycotoxins such as nivalenol (NIV) and deoxynivalenol (DON). The aim of this study was to conduct population studies on F. asiaticum from Southern China through mitochondrial genome analyses. All strains were isolated from wheat or rice from several geographic areas in seven provinces in Southern China. Based on geographic location and host, 210 isolates were selected for next generation sequencing, and their mitogenomes were assembled by GRAbB and annotated to explore the mitochondrial genome variability of F. asiaticum. The F. asiaticum mitogenome proves extremely conserved and variation is mainly caused by absence/presence of introns harboring homing endonuclease genes. These variations could be utilized to develop molecular markers for track and trace of migrations within and between populations. This study illustrates how mitochondrial introns can be used as markers for population genetic analysis. SNP analysis demonstrate the occurrence of mitochondrial recombination in F. asiaticum as was previously found for F. oxysporum and implied for F. graminearum. Furthermore, varying degrees of genetic diversity and recombination showed a high association with different geographic regions as well as with cropping systems. The mitogenome of F. graminearum showed a much higher SNP diversity while the interspecies intron variation showed no evidence of gene flow between the two closely related and sexual compatible species.

Characterization of Nuclear and Mitochondrial Genomes of Two Tobacco Endophytic Fungi Leptosphaerulina chartarum and Curvularia trifolii and Their Contributions to Phylogenetic Implications in the Pleosporales

The symbiont endophytic fungi in tobacco are highly diverse and difficult to classify. Here, we sequenced the genomes of Curvularia trifolii and Leptosphaerulina chartarum isolated from tobacco plants. Finally, 41.68 Mb and 37.95 Mb nuclear genomes were sequenced for C. trifolii and L. chartarum with the scaffold N50, accounting for 638.94 Kb and 284.12 Kb, respectively. Meanwhile, we obtained 68,926 bp and 59,100 bp for their mitochondrial genomes. To more accurately classify C. trifolii and L. chartarum, we extracted seven nuclear genes and 12 mitochondrial genes from these two genomes and their closely related species. The genes were then used for calculation of evolutionary rates and for phylogenetic analysis. Results showed that it was difficult to achieve consistent results using a single gene due to their different evolutionary rates, while the phylogenetic trees obtained by combining datasets showed stable topologies. It is, therefore, more accurate to construct phylogenetic relationships for endophytic fungi based on multi-gene datasets. This study provides new insights into the distribution and characteristics of endophytic fungi in tobacco.

Complete mitochondrial genome sequence of lettuce pathogenic fungus, Fusarium oxysporum f. sp. lactucae 09-002

Fusarium oxysporum is a famous plant pathogenic filamentous fungus. Here, we report the complete mitochondrial genome sequence of F. oxysporum f. sp. lactucae isolated from the lettuce field in Suwon area, Korea. Total length of the mitochondrial genome is 45,020 bp and it encodes 42 genes (15 protein-coding genes, two rRNAs, and 25 tRNAs). Nucleotide sequence of coding region takes over 32.7%, and overall GC content is 32.4%. Phylogenetic tree presented that F. oxysporum f. sp. lactucae 09-002 was clustered with Fusarium commune not like another F. oxysporum mitochondrial genomes, requiring further analyses.

A Sequence Resource of Autosomes and Additional Chromosomes in the Peach Pathotype of Alternaria alternata

Alternaria alternata is a generally saprophytic fungus. Its genome consists of 10 autosomes, while some strains have one or two additional chromosomes, called a conditionally dispensable chromosome (CDC). A CDC is not required for reproduction but confers host-specific pathogenicity. We sequenced the genome of the peach pathotype of A. alternata using Nanopore and assembled it into 20 sequences. The 10 largest sequences corresponded to 10 gapless sequences of A. solani autosomes, and 1 sequence was a mitochondrial genome. Nine other sequences may be derived from CDCs because of lack of similarity with autosomes of other Alternaria spp. The sequence information could provide novel insights into genomes of Alternaria spp. and on the biosynthesis of a novel host-specific toxin in the peach pathotype of A. alternata.

Complete mitochondrial genome sequence of lettuce pathogenic fungus, Fusarium oxysporum f. sp. lactucae 16-086

Fusarium oxysporum f. sp. lactucae 016-086 is a plant pathogenic filamentous fungus isolated from wilted lettuce in Korea. We reported complete mitochondrial genome sequence of F. oxysporum f. sp. lactucae 016-086. Total length of this mitogenome is 41,826 bp and it encoded 42 genes (14 protein-coding genes, 2 rRNAs, and 26 tRNAs). Nucleotide sequence of coding region takes over 30.6%, and overall GC content is 32.5%. Phylogenetic tree of Fusarium mitochondrial genomes presented distinct clades along with nine formae speciales. This mitogenome will contribute distinguishing formae speciales of F. oxysoporum claearly with additional mitogenomes sequenced in the near future.

Maternal mitochondrial inheritance in two Fusarium pathogens of prickly ash (Zanthoxylum bungeanum) in northern China

Mitochondrial inheritance in Fusarium zanthoxyli and F. continuum, two canker-inducing pathogens of prickly ash (Zanthoxylum bungeanum) in northern China, was investigated by genotyping ascospore progeny obtained from laboratory crosses. Polymorphic regions of the mitochondrial genomes (mitogenomes) that contained indels and single-nucleotide polymorphisms (SNPs) were identified via comparative analyses of the complete mitogenomes of the parents used in the intraspecific crosses. A reciprocal genetic cross of F. zanthoxyli NRRL 66714 × NRRL 66285, and a separate cross of F. continuum ♀ NRRL 66286 × ♂ NRRL 66218, revealed that mitochondria were only inherited from the maternal parent. In addition, the reciprocal cross demonstrated that mitochondrial inheritance is not linked to mating type. Gene order in the circular mitogenomes of the prickly ash pathogens was identical to that previously reported for other fusaria and members of the Hypocreales, except that the TRNL tRNAs were duplicated in F. zanthoxyli NRRL 66714. The genomes contained 14 polypeptide-encoding genes involved in oxidative respiration, one intron-encoded ribosomal protein (rps3) gene, two ribosomal RNA (rRNA) genes, and 26-28 tRNA genes. The F. zanthoxyli mitogenomes were 80.9 and 98.7 kb in length, whereas those of F. continuum were considerably shorter and nearly identical in length at 63.4 kb. The significant differences in mitogenome length were primarily due to variable numbers of introns and open reading frames (ORFs) encoding hypothetical proteins.

Genome analysis of the ubiquitous boxwood pathogen Pseudonectria foliicola

Boxwood (Buxus spp.) are broad-leaved, evergreen landscape plants valued for their longevity and ornamental qualities. Volutella leaf and stem blight, caused by the ascomycete fungi Pseudonectria foliicola and P. buxi, is one of the major diseases affecting the health and ornamental qualities of boxwood. Although this disease is less severe than boxwood blight caused by Calonectria pseudonaviculata and C. henricotiae, its widespread occurrence and disfiguring symptoms have caused substantial economic losses to the ornamental industry. In this study, we sequenced the genome of P. foliicola isolate ATCC13545 using Illumina technology and compared it to other publicly available fungal pathogen genomes to better understand the biology of this organism. A de novo assembly estimated the genome size of P. foliicola at 28.7 Mb (425 contigs; N50 = 184,987 bp; avg. coverage 188×), with just 9,272 protein-coding genes. To our knowledge, P. foliicola has the smallest known genome within the Nectriaceae. Consistent with the small size of the genome, the secretome, CAzyme and secondary metabolite profiles of this fungus are reduced relative to two other surveyed Nectriaceae fungal genomes: Dactylonectria macrodidyma JAC15-245 and Fusarium graminearum Ph-1. Interestingly, a large cohort of genes associated with reduced virulence and loss of pathogenicity was identified from the P. foliicola dataset. These data are consistent with the latest observations by plant pathologists that P. buxi and most likely P. foliicola, are opportunistic, latent pathogens that prey upon weak and stressed boxwood plants.

Mitochondrial introgression and interspecies recombination in the Fusarium fujikuroi species complex

The Fusarium fujikuroi species complex (FFSC) is an economically important monophyletic lineage in the genus Fusarium. Incongruence observed among mitochondrial gene trees, as well as the multiple non-orthologous copies of the internal transcribed spacer region of the ribosomal RNA genes, suggests that the origin and history of this complex likely involved interspecies gene flow. Based on this hypothesis, the mitochondrial genomes of non-conspecific species should harbour signatures of introgression or introgressive hybridization. The aim of this study was therefore to search for recombination between the mitochondrial genomes of different species in the FFSC. Using methods based on mt genome sequence similarity, five significant recombinant regions in both gene and intergenic regions were detected. Using coalescent-based methods and the sequences for individual mt genes, various ancestral recombination events between different lineages of the FFSC were also detected. These findings suggest that interspecies gene flow and introgression are likely to have played key roles in the evolution of the FFSC at both ancient and more recent time scales.

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

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

Complete mitogenome of the high ethanol production fungus Fusarium oxysporum Mh2-2

Fusarium spp. are significantly important plant pathogens, and some of them are ethanol-producing strains. During infection and/or ethanol production, Fusarium requires a plenty of energy that is mainly provided by mitochondria. Here we report the first mitogenome from a selected Fusarium oxysporum strain mh2-2 that produces ethanol from glucose and xylose. The size of this mitogenome, 46 kb, is different from the size of any reported Fusarium mitogenome. Our results provide insight into the functions and evolution of mitochondrial genes and genomes.

van der Lee TA, Waalwijk C, van Diepeningen AD. Mitochondrial genomes reveal recombination in the presumed asexual Fusarium oxysporum species complex

BACKGROUND

The Fusarium oxysporum species complex (FOSC) contains several phylogenetic lineages. Phylogenetic studies identified two to three major clades within the FOSC. The mitochondrial sequences are highly informative phylogenetic markers, but have been mostly neglected due to technical difficulties.

RESULTS

A total of 61 complete mitogenomes of FOSC strains were de novo assembled and annotated. Length variations and intron patterns support the separation of three phylogenetic species. The variable region of the mitogenome that is typical for the genus Fusarium shows two new variants in the FOSC. The variant typical for Fusarium is found in members of all three clades, while variant 2 is found in clades 2 and 3 and variant 3 only in clade 2. The extended set of loci analyzed using a new implementation of the genealogical concordance species recognition method support the identification of three phylogenetic species within the FOSC. Comparative analysis of the mitogenomes in the FOSC revealed ongoing mitochondrial recombination within, but not between phylogenetic species.

CONCLUSIONS

The recombination indicates the presence of a parasexual cycle in F. oxysporum. The obstacles hindering the usage of the mitogenomes are resolved by using next generation sequencing and selective genome assemblers, such as GRAbB. Complete mitogenome sequences offer a stable basis and reference point for phylogenetic and population genetic studies.

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

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

Intron Derived Size Polymorphism in the Mitochondrial Genomes of Closely Related Chrysoporthe Species

In this study, the complete mitochondrial (mt) genomes of Chrysoporthe austroafricana (190,834 bp), C. cubensis (89,084 bp) and C. deuterocubensis (124,412 bp) were determined. Additionally, the mitochondrial genome of another member of the Cryphonectriaceae, namely Cryphonectria parasitica (158,902 bp), was retrieved and annotated for comparative purposes. These genomes showed high levels of synteny, especially in regions including genes involved in oxidative phosphorylation and electron transfer, unique open reading frames (uORFs), ribosomal RNAs (rRNAs) and transfer RNAs (tRNAs), as well as intron positions. Comparative analyses revealed signatures of duplication events, intron number and length variation, and varying intronic ORFs which highlighted the genetic diversity of mt genomes among the Cryphonectriaceae. These mt genomes showed remarkable size polymorphism. The size polymorphism in the mt genomes of these closely related Chrysoporthe species was attributed to the varying number and length of introns, coding sequences and to a lesser extent, intergenic sequences. Compared to publicly available fungal mt genomes, the C. austroafricana mt genome is the second largest in the Ascomycetes thus far.

Mitochondrial genome evolution in yeasts: an all-encompassing view

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

Analysis of the complete mitochondrial genome of Pochonia chlamydosporia suggests a close relationship to the invertebrate-pathogenic fungi in Hypocreales

Background

The fungus Pochonia chlamydosporia parasitizes nematode eggs and has become one of the most promising biological control agents (BCAs) for plant-parasitic nematodes, which are major agricultural pests that cause tremendous economic losses worldwide. The complete mitochondrial (mt) genome is expected to open new avenues for understanding the phylogenetic relationships and evolution of the invertebrate-pathogenic fungi in Hypocreales.

Results

The complete mitogenome sequence of P. chlamydosporia is 25,615 bp in size, containing the 14 typical protein-coding genes, two ribosomal RNA genes, an intronic ORF coding for a putative ribosomal protein (rps3) and a set of 23 transfer RNA genes (trn) which recognize codons for all amino acids. Sequence similarity studies and syntenic gene analyses show that 87.02% and 58.72% of P. chlamydosporia mitogenome sequences match 90.50% of Metarhizium anisopliae sequences and 61.33% of Lecanicillium muscarium sequences with 92.38% and 86.04% identities, respectively. A phylogenetic tree inferred from 14 mt proteins in Pezizomycotina fungi supports that P. chlamydosporia is most closely related to the entomopathogenic fungus M. anisopliae. The invertebrate-pathogenic fungi in Hypocreales cluster together and clearly separate from a cluster comprising plant-pathogenic fungi (Fusarium spp.) and Hypocrea jecorina. A comparison of mitogenome sizes shows that the length of the intergenic regions or the intronic regions is the major size contributor in most of mitogenomes in Sordariomycetes. Evolutionary analysis shows that rps3 is under positive selection, leading to the display of unique evolutionary characteristics in Hypocreales. Moreover, the variability of trn distribution has a clear impact on gene order in mitogenomes. Gene rearrangement analysis shows that operation of transposition drives the rearrangement events in Pezizomycotina, and most events involve in trn position changes, but no rearrangement was found in Clavicipitaceae.

Conclusions

We present the complete annotated mitogenome sequence of P. chlamydosporia. Based on evolutionary and phylogenetic analyses, we have determined the relationships between the invertebrate-pathogenic fungi in Hypocreales. The invertebrate-pathogenic fungi in Hypocreales referred to in this paper form a monophyletic group sharing a most recent common ancestor. Our rps3 and trn gene order results also establish a foundation for further exploration of the evolutionary trajectory of the fungi in Hypocreales.

Electronic supplementary material

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

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Complete mitochondrial genome of the cephalosporin-producing fungus Acremonium chrysogenum

Complete nucleotide sequence of the 27,266 bp mitochondrial genome of the сephalosporin C producing fungus Acremonium chrysogenum have been determined using whole genome shotgun sequencing approach. The circular mapping molecule encodes a usual set of mitochondrial proteins and RNA genes, including large and small ribosomal RNAs, 19 proteins and 26 tRNA genes and contains 2 introns. All structural genes are located on one strand and are apparently transcribed in one direction. Comparative analysis of this and previously sequenced Pezizomycotina mtDNAs revealed more extensive similarity between A. chrysogenum genome and those of Fusarium clade and specific synthenic patterns characteristic of Hypocrealean mitogenomes. Phylogenetic analysis based on catenated mitochondrial protein sequences confirmed current taxonomic position of A. chrysogenum within Hyprocreales and related taxa.

Evidence for inter-specific recombination among the mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex

Background

The availability of mitochondrial genomes has allowed for the resolution of numerous questions regarding the evolutionary history of fungi and other eukaryotes. In the Gibberella fujikuroi species complex, the exact relationships among the so-called “African”, “Asian” and “American” Clades remain largely unresolved, irrespective of the markers employed. In this study, we considered the feasibility of using mitochondrial genes to infer the phylogenetic relationships among Fusarium species in this complex. The mitochondrial genomes of representatives of the three Clades (Fusarium circinatum, F. verticillioides and F. fujikuroi) were characterized and we determined whether or not the mitochondrial genomes of these fungi have value in resolving the higher level evolutionary relationships in the complex.

Results

Overall, the mitochondrial genomes of the three species displayed a high degree of synteny, with all the genes (protein coding genes, unique ORFs, ribosomal RNA and tRNA genes) in identical order and orientation, as well as introns that share similar positions within genes. The intergenic regions and introns generally contributed significantly to the size differences and diversity observed among these genomes. Phylogenetic analysis of the concatenated protein-coding dataset separated members of the Gibberella fujikuroi complex from other Fusarium species and suggested that F. fujikuroi (“Asian” Clade) is basal in the complex. However, individual mitochondrial gene trees were largely incongruent with one another and with the concatenated gene tree, because six distinct phylogenetic trees were recovered from the various single gene datasets.

Conclusion

The mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex are remarkably similar to those of the previously characterized Fusarium species and Sordariomycetes. Despite apparently representing a single replicative unit, all of the genes encoded on the mitochondrial genomes of these fungi do not share the same evolutionary history. This incongruence could be due to biased selection on some genes or recombination among mitochondrial genomes. The results thus suggest that the use of individual mitochondrial genes for phylogenetic inference could mask the true relationships between species in this complex.

Multiple copies of cytochrome oxidase 1 in species of the fungal genus Fusarium

Using data from published mitochondrial or complete genomes, we developed and tested primers for amplification and sequencing of the barcode region of cytochrome oxidase 1 (COX1) of the fungal genus Fusarium, related genera of the order Hypocreales, and degenerate primers for fungi in the subdivision Pezizomycotina. The primers were successful for amplifying and sequencing COX1 barcodes from 13 genera of Hypocreales (Acremonium, Beauveria, Clonostachys, Emericellopsis, Fusarium, Gliocladium, Hypocrea, Lanatonectria, Lecanicillium, Metarhizium, Monocillium, Neonectria and Stilbella), 22 taxa of Fusarium, and two genera in other orders (Arthrosporium, Monilochaetes). Parologous copies of COX1 occurred in several strains of Fusarium. In some, copies of the same length were detected either by heterozygous bases in otherwise clean sequences or in different replicates of amplification and sequencing events; this may indicate multiple transcribed copies. Other strains included one or two introns. Two intron insertion sites had at least two nonhomologous intron sequences among Fusarium species. Irrespective of whether the multiple copy issue could be resolved by sequencing RNA transcripts, developing a precise COX1-based barcoding system for Fusarium may not be feasible. The overall divergence among homologous COX1 sequences obtained so far is rather low, with many species sharing identical sequences.

Global analyses of Ceratocystis cacaofunesta mitochondria: from genome to proteome

Background

The ascomycete fungus Ceratocystis cacaofunesta is the causal agent of wilt disease in cacao, which results in significant economic losses in the affected producing areas. Despite the economic importance of the Ceratocystis complex of species, no genomic data are available for any of its members. Given that mitochondria play important roles in fungal virulence and the susceptibility/resistance of fungi to fungicides, we performed the first functional analysis of this organelle in Ceratocystis using integrated “omics” approaches.

Results

The C. cacaofunesta mitochondrial genome (mtDNA) consists of a single, 103,147-bp circular molecule, making this the second largest mtDNA among the Sordariomycetes. Bioinformatics analysis revealed the presence of 15 conserved genes and 37 intronic open reading frames in C. cacaofunesta mtDNA. Here, we predicted the mitochondrial proteome (mtProt) of C. cacaofunesta, which is comprised of 1,124 polypeptides - 52 proteins that are mitochondrially encoded and 1,072 that are nuclearly encoded. Transcriptome analysis revealed 33 probable novel genes. Comparisons among the Gene Ontology results of the predicted mtProt of C. cacaofunesta, Neurospora crassa and Saccharomyces cerevisiae revealed no significant differences. Moreover, C. cacaofunesta mitochondria were isolated, and the mtProt was subjected to mass spectrometric analysis. The experimental proteome validated 27% of the predicted mtProt. Our results confirmed the existence of 110 hypothetical proteins and 7 novel proteins of which 83 and 1, respectively, had putative mitochondrial localization.

Conclusions

The present study provides the first partial genomic analysis of a species of the Ceratocystis genus and the first predicted mitochondrial protein inventory of a phytopathogenic fungus. In addition to the known mitochondrial role in pathogenicity, our results demonstrated that the global function analysis of this organelle is similar in pathogenic and non-pathogenic fungi, suggesting that its relevance in the lifestyle of these organisms should be based on a small number of specific proteins and/or with respect to differential gene regulation. In this regard, particular interest should be directed towards mitochondrial proteins with unknown function and the novel protein that might be specific to this species. Further functional characterization of these proteins could enhance our understanding of the role of mitochondria in phytopathogenicity.

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.

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.

Mitochondrial genomes of yeasts of the Yarrowia clade

Candida alimentaria, Candida deformans, Candida galli, and Candida phangngensis have been recently reported to be the close relatives of Yarrowia lipolytica. To explore this clade of yeasts, we sequenced the mitochondrial genome (mtDNA) of these four species and compared it with the mtDNA of Y. lipolytica. The five mtDNAs exhibit a similar architecture and a high level of similarity of protein coding sequences. Genome sizes are variable, ranging from 28 017 bp in C. phangngensis to 48 508 bp in C. galli, mainly because of the variations in intron size and number. All introns are of group I, except for a group II intron inserted in the cob gene of a single species, C. galli. Putative endonuclease coding sequences were present in most group I introns, but also twice as free-standing ORFs in C. galli. Phylogenetic relationships of the five species were explored using protein alignments. No close relative of the Yarrowia clade could be identified, but protein and rRNA gene orders were partially conserved in the mtDNA of Candida salmanticensis.

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.

LTR retrotransposons in fungi

Transposable elements with long terminal direct repeats (LTR TEs) are one of the best studied groups of mobile elements. They are ubiquitous elements present in almost all eukaryotic genomes. Their number and state of conservation can be a highlight of genome dynamics. We searched all published fungal genomes for LTR-containing retrotransposons, including both complete, functional elements and remnant copies. We identified a total of over 66,000 elements, all of which belong to the Ty1/Copia or Ty3/Gypsy superfamilies. Most of the detected Gypsy elements represent Chromoviridae, i.e. they carry a chromodomain in the pol ORF. We analyzed our data from a genome-ecology perspective, looking at the abundance of various types of LTR TEs in individual genomes and at the highest-copy element from each genome. The TE content is very variable among the analyzed genomes. Some genomes are very scarce in LTR TEs (<50 elements), others demonstrate huge expansions (>8000 elements). The data shows that transposon expansions in fungi usually involve an increase both in the copy number of individual elements and in the number of element types. The majority of the highest-copy TEs from all genomes are Ty3/Gypsy transposons. Phylogenetic analysis of these elements suggests that TE expansions have appeared independently of each other, in distant genomes and at different taxonomical levels. We also analyzed the evolutionary relationships between protein domains encoded by the transposon pol ORF and we found that the protease is the fastest evolving domain whereas reverse transcriptase and RNase H evolve much slower and in correlation with each other.

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.

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.