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

A Novel Report of mycetoma with Spinal Spread due to Madurella fahalli from India

Eumycetoma caused by Madurella fahalii, a drug-resistant fungus, has never been reported in India. Here, we describe a fatal case of eumycetoma with spinal involvement due to M. fahalii for the first time in India.

The Mitogenomic Characterization and Phylogenetic Analysis of the Plant Pathogen Phyllosticta yuccae

Phyllosticta yuccae is an important plant pathogen causing leaf spot disease in Yucca gigantea Lem. It is imperative to note that the amount of information available about the mitogenome of this subject is severely limited. This must be addressed immediately, as it is crucial to our understanding and progress in this field. To better understand the mitogenomic characteristics of P. yuccae, we conducted its sequencing by MGISEQ. Afterwards, the mitogenome was assembled and annotated. The mitogenomic characteristics and phylogenetic placement of the P. yuccae strain KUMCC 6213 were analyzed. The study revealed that the mitogenome of P. yuccae is a circular DNA molecule, consisting of 178,540 base pairs. It contains a total of 64 genes, including 14 protein-coding genes (PCGs), 26 transfer RNA genes (tRNA), 2 ribosomal RNA genes (rRNA), and 22 open reading frame genes (ORF), accounting for 80.98% of the total size. Repetitive sequences accounted for 15.42% of the mitogenome. The analysis of codon usage indicated that the codon UUA was the most commonly utilized, whereas the amino acid Leu was the most frequently employed. A comparative analysis of mitogenomes between P. yuccae and Macrophomina phaseolina (Tassi) Goid. showed notable variations in the position and size of gene clusters, with cox1, nad4, and nad4L genes exhibiting relatively low conservation. Phylogenetic analysis based on the 14 PCGs revealed that P. yuccae has the closest genetic relationship with M. phaseolina (Botryosphaeriaceae, Botryosphaeriales). This study first reports the mitogenome of P. yuccae and validates its phylogenetic placement. The findings enhance the knowledge of mitogenomes in Botryosphaeriales, offering novel perspectives on the genetics and evolution of the plant pathogen P. yuccae. This is crucial for the accurate prevention and management of leaf spot disease in Y. gigantea.

The combination of manogepix and itraconazole is synergistic and inhibits the growth of Madurella mycetomatis in vitro but not in vivo

Mycetoma is a neglected tropical disease commonly caused by the fungus Madurella mycetomatis. Standard treatment consists of extensive treatment with itraconazole in combination with surgical excision of the infected tissue, but has a low success rate. To improve treatment outcomes, novel treatment strategies are needed. Here, we determined the potential of manogepix, a novel antifungal agent that targets the GPI-anchor biosynthesis pathway by inhibition of the GWT1 enzyme. Manogepix was evaluated by determining the minimal inhibitory concentrations (MICs) according to the CLSI-based in vitro susceptibility assay for 22 M. mycetomatis strains and by in silico protein comparison of the target protein. The synergy between manogepix and itraconazole was determined using a checkerboard assay. The efficacy of clinically relevant dosages was assessed in an in vivo grain model in Galleria mellonella larvae. MICs for manogepix ranged from <0.008 to >8 mg/l and 16/22 M. mycetomatis strains had an MIC ≥4 mg/ml. Differences in MICs were not related to differences observed in the GWT1 protein sequence. For 70% of the tested isolates, synergism was found between manogepix and itraconazole in vitro. In vivo, enhanced survival was not observed upon admission of 8.6 mg/kg manogepix, nor in combination treatment with 5.7 mg/kg itraconazole. MICs of manogepix were high, but the in vitro antifungal activity of itraconazole was enhanced in combination therapy. However, no efficacy of manogepix was found in an in vivo grain model using clinically relevant dosages. Therefore, the therapeutic potential of manogepix in mycetoma caused by M. mycetomatis seems limited.

The complete mitochondrial genome of the rice blast fungus Pyricularia oryzae Cavara 1892 strain Guy11 and phylogenetic analysis

The complete mitochondrial genome of Pyricularia oryzae Cavara 1892 strain Guy11 is 34,865 bp in length (GenBank accession number OP095391), containing 29 tRNA genes, 2 rRNA genes, and 15 protein-coding genes (PCGs). The gene order and orientation are novel compared to other Sordariomycetes species with sequenced mitogenomes in the GenBank database. Phylogenetic analysis suggests that P. oryzae Guy11 and 19 other Sordariomycetes species form a monophyletic group. The complete mitochondrial sequence of P. oryzae Guy11 will be a valuable resource for species identification, population genetics, phylogenetics, and comparative genomics studies in Sordariomycetes and Magnaporthales.

Trichoderma koningiopsis (Hypocreaceae) has the smallest mitogenome of the genus Trichoderma

Introduction

Fungal mitogenomes exhibit remarkable variation in conformation, size, gene content, arrangement and expression, including their intergenic spacers and introns.

Methods

The complete mitochondrial genome sequence of the mycoparasitic fungus Trichoderma koningiopsis was determined using the Illumina next-generation sequencing technology. We used data from our recent Illumina NGS-based project of T. koningiopsis genome sequencing to study its mitochondrial genome. The mitogenome was assembled, annotated, and compared with other fungal mitogenomes.

Results

T. koningiopsis strain POS7 mitogenome is a circular molecule of 27,560 bp long with a GC content of 27.80%. It harbors the whole complement of the 14 conserved mitochondrial protein-coding genes (PCG) such as atp6, atp8, atp9, cox1, cox2, cox3, cob, nad1, nad2, nad3, nad4, nad4L, nad5, and nad6, also found in the same gene order to other Hypocreales. The mitogenome also contains 26 transfer RNA genes (tRNAs), 5 of them with more than one copy. Other genes also present in the assembled mitochondrial genome are a small rRNA subunit and a large rRNA subunit containing ribosomal protein S3 gene. Despite the small genome size, two introns were detected in the T. koningiopsis POS7 mitogenome, one of them in cox3 gene and the other in rnl gene, accounting 7.34% of this mitogenome with a total size of 2,024 bp. A phylogenetic analysis was done using the 14 PCGs genes of T. koningiopsis strain POS7 mitogenome to compare them with those from other fungi of the Subphyla Pezizomycotina and Saccharomycotina. T. koningiopsis strain POS7 was clustered together with other representatives of Trichoderma lineage, within the Hypocreales group, which is also supported by previous phylogenetic studies based on nuclear markers.

Discussion

The mitochondrial genome of T. koningiopsis POS7 will allow further investigations into the taxonomy, phylogenetics, conservation genetics, and evolutionary biology of this important genus as well as other closely related species.

Functional annotation and comparative analysis of four Botrytis cinerea mitogenomes reported from Punjab, Pakistan

Botrytis cinerea is one of the top phytopathogenic fungus which ubiquitously cause grey mold on a variety of horticultural plants. The mechanism of respiration in the fungus occurs within the mitochondria. Mitogenomes serve as a key molecular marker for the investigation of fungal evolutionary patterns. This study aimed at the complete assembly, characterization, and comparative relationship of four mitogenomes of Botrytis cinerea strains including Kst5C, Kst14A, Kst32B, Kst33A, respectively. High throughput sequencing of four mitogenomes allowed the full assembly and annotation of these sequences. The total genome length of these 4 isolates Kst5C Kst14A, Kst32B, Kst33A was 69,986 bp, 77,303 bp, 76,204 bp and 55, 226 bp respectively. The distribution of features represented 2 ribosomal RNA genes,14 respiration encoding proteins, 1 mitochondrial ribosomal protein-encoding gene, along with varying numbers of transfer RNA genes, protein-coding genes, mobile intronic regions and homing endonuclease genes including LAGLIDADG and GIY-YIG domains were found in all four mitogenomes. The comparative analyses performed also decipher significant results for four mitogenomes among fungal isolates included in the study. This is the first report on the detailed annotation of mitogenomes as a proof for investigation of variation patterns present with in the B. cinerea causing grey mold on strawberries in Pakistan. This study will also contribute to the rapid evolutionary analysis and population patterns present among Botrytis cinerea.

Multifractal and cross-correlation analysis on mitochondrial genome sequences using chaos game representation

We characterized the multifractality and power-law cross-correlation of mitochondrial genomes of various species through the recently developed method which combines the chaos game representation theory and 2D-multifractal detrended cross-correlation analysis. In the present paper, we analyzed 32 mitochondrial genomes of different species and the obtained results show that all the analyzed data exhibit multifractal nature and power-law cross-correlation behaviour. Further, we performed a cluster analysis from the calculated scaling exponents to identify the class affiliation and its outcome is represented as a dendrogram. We suggest that this integrative approach may help the researchers to understand the phylogeny of any kingdom with their varying genome lengths and also this approach may find applications in characterizing the protein sequences, mRNA sequences, next-generation sequencing, and drug development, etc.

Revisiting the Neurospora crassa mitochondrial genome

The mitochondrial genome of Neurospora crassa has been less studied than its nuclear counterpart, yet it holds great potential for understanding the diversity and evolution of this important fungus. Here we describe a new mitochondrial DNA (mtDNA) complete sequence of a N. crassa wild type strain. The genome with 64 839 bp revealed 21 protein-coding genes and several hypothetical open reading frames with no significant homology to any described gene. Five large repetitive regions were identified across the genome, including partial or complete genes. The largest repeated region holds a partial nd2 section that was also detected in Neurospora intermedia, suggesting a rearrangement that occurred before the N. crassa speciation. Interestingly, N. crassa has a palindrome adjacent to the partial nd2 repeated region possibly related to the genomic rearrangement, which is absent in N. intermedia. Finally, we compared the sequences of the three available N. crassa complete mtDNAs and found low levels of intraspecific variability. Most differences among strains were due to small indels in noncoding regions. The revisiting of the N. crassa mtDNA forms the basis for future studies on mitochondrial genome organization and variability.

Ranked k-Spectrum Kernel for Comparative and Evolutionary Comparison of Exons, Introns, and CpG Islands

MOTIVATION

Existing k-mer based string kernel methods have been successfully used for sequence comparison. However, existing kernel methods have limitations for comparative and evolutionary comparisons of genomes due to the sensitiveness to over-represented k-mers and variable sequence lengths.

RESULTS

In this study, we propose a novel ranked k-spectrum string (RKSS) kernel. 1) RKSS kernel utilizes common k-mer sets across species, named landmarks, that can be used for comparing multiple genomes. 2) Based on the landmarks, we can use ranks of k-mers, rather than frequencies, that can produce more robust distances between genomes. To show the power of RKSS kernel, we conducted two experiments using 10 mammalian species with exon, intron, and CpG island sequences. RKSS kernel reconstructed more consistent evolutionary trees than the k-spectrum string kernel. In the subsequent experiment, for each sequence, kernel distance was calculated from 30 landmarks representing exon, intron, and CpG island sequences of 10 genomes. Based on kernel distances, concordance tests were performed and the result suggested that more information is conserved in CpG islands across species than in introns. In conclusion, our analysis suggests that the relational order, exon CpG island intron, in terms of evolutionary information contents.

Fungal Mitogenomes: Relevant Features to Planning Plant Disease Management

Mitochondrial genomes (mt-genomes) are characterized by a distinct codon usage and their autonomous replication. Mt-genomes encode highly conserved genes (mt-genes), like proteins involved in electron transport and oxidative phosphorylation but they also carry highly variable regions that are in part responsible for their high plasticity. The degree of conservation of their genes is such that they allow the establishment of phylogenetic relationships even across distantly related species. Here, we describe the mechanisms that generate changes along mt-genomes, which play key roles at enlarging the ability of fungi to adapt to changing environments. Within mt-genomes of fungal pathogens, there are dispensable as well as indispensable genes for survival, virulence and/or pathogenicity. We also describe the different complexes or mechanisms targeted by fungicides, thus addressing a relevant issue regarding disease management. Despite the controversial origin and evolution of fungal mt-genomes, the intrinsic mechanisms and molecular biology involved in their evolution will help to understand, at the molecular level, the strategies for fungal disease management.

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.

Genomic characterization of Parengyodontium americanum sp. nov

Modern genome analysis and phylogenomic methods have increased the number of fungal species, as well as enhanced appreciation of the degree of diversity within the fungal kingdom. In this context, we describe a new Parengyodontium species, P. americanum, which is phylogenetically related to the opportunistic human fungal pathogen P. album. Five unusual fungal isolates were recovered from five unique and confirmed coccidioidomycosis patients, and these isolates were subsequently submitted to detailed molecular and morphological identification procedures to determine identity. Molecular and morphological diagnostic analyses showed that the isolates belong to the Cordycipitaceae. Subsequently, three representative genomes were sequenced and annotated, and a new species, P. americanum, was identified. Using various genomic analyses, gene family expansions related to novel compounds and potential for ability to grow in diverse habitats are predicted. A general description of the genomic composition of this newly described species and comparison of genome content with Beauveria bassiana, Isaria fumosorosea and Cordyceps militaris shows a shared core genome of 6371 genes, and 148 genes that appear to be specific for P. americanum. This work provides the framework for future investigations of this interesting fungal species.

The complete mitochondrial genome of the Chan-hua fungus Isaria cicadae: a tale of intron evolution in Cordycipitaceae

Isaria cicadae is an entomogenous fungus of great medicinal value. Its nuclear genome has been reported, while its mitogenome remains unknown. Herein, we first described its mitogenome and then inferred intron evolution from both intraspecific and interspecific perspectives. The fungus represented the largest mitogenome (56.6 kb in strain CCAD02) known in Cordycipitaceae due to the presence of 25 introns interrupting nine genes. Comparison of three I. cicadae strains revealed intron presence/absence dynamics at six intron loci plus a few indels and single nucleotide polymorphisms. Phylogenetic analyses confirmed the placement of I. cicadae in Cordycipitaceae. Comparison of 10 Cordycipitaceae species revealed a high degree of synteny and conserved genetic content. They, however, varied in intron numbers (1-25 per species) with overall 34 intron loci identified, which resulted in more than twofold variations in mitogenome sizes (24.5-56.6 kb). An rnl intron encoding ribosomal protein S3 was present in all species, suggesting its early invasion in Cordycipitaceae, while further divergence occurred for this intron. The other introns identified in this study were present in some, but not all of the species and have undergone multiple gains and losses in Cordycipitaceae. This study greatly enhanced our understanding of intron evolution in Cordycipitaceae.

Mitochondrial genome of the entomophthoroid fungus Conidiobolus heterosporus provides insights into evolution of basal fungi

Entomophthoroid fungi represent an ecologically important group of fungal pathogens on insects. Here, the whole mitogenome of Conidiobolus heterosporus, one of the entomophthoroid fungi, was described and compared to those early branching fungi with available mitogenomes. The 53,364-bp circular mitogenome of C. heterosporus contained two rRNA genes, 14 standard protein-coding genes, 26 tRNA genes, and three free-standing ORFs. Thirty introns interrupted nine mitochondrial genes. Phylogenetic analysis based on mitochondrion-encoded proteins revealed that C. heterosporus was most close to Zancudomyces culisetae in the Zoopagomycota of basal fungi. Comparison on mitogenomes of 23 basal fungi revealed great variabilities in terms of mitogenome conformation (circular or linear), genetic code (codes 1, 4, or 16), AT contents (53.3-85.5%), etc. These mitogenomes varied from 12.0 to 97.3 kb in sizes, mainly due to different numbers of genes and introns. They showed frequent DNA rearrangement events and a high variability of gene order, although high synteny and conserved gene order were also present between closely related species. By reporting the first mitogenome in Entomophthoromycotina and the second in Zoopagomycota, this study greatly enhanced our understanding on evolution of basal fungi.

The mitochondrial genome of Endoconidiophora resinifera is intron rich

Endoconidiophora resinifera (=Ceratocystis resinifera) is a blue-stain fungus that occurs on conifers. The data showed that the Endoconidiophora resinifera mitochondrial genome is one of the largest mitochondrial genomes (>220 kb) so far reported among members of the Ascomycota. An exceptional large number of introns (81) were noted and differences among the four strains were restricted to minor variations in intron numbers and a few indels and single nucleotide polymorphisms. The major differences among the four strains examined are due to size polymorphisms generated by the absence or presence of mitochondrial introns. Also, these mitochondrial genomes encode the largest cytochrome oxidase subunit 1 gene (47.5 kb) reported so far among the fungi. The large size for this gene again can be attributed to the large number of intron insertions. This study reports the first mitochondrial genome for the genus Endoconidiophora, previously members of this genus were assigned to Ceratocystis. The latter genus has recently undergone extensive taxonomic revisions and the mitochondrial genome might provide loci that could be applied as molecular markers assisting in the identification of taxa within this group of economically important fungi. The large mitochondrial genome also may provide some insight on mechanisms that can lead to mitochondrial genome expansion.

Comparative analysis of the mitochondrial genome of the fungus Colletotrichum lindemuthianum, the causal agent of anthracnose in common beans

Fungi of the genus Colletotrichum are economically important and are used as models in plant-pathogen interaction studies. In this study, the complete mitochondrial genomes of two Colletotrichum lindemuthianum isolates were sequenced and compared with the mitochondrial genomes of seven species of Colletotrichum. The mitochondrial genome of C. lindemuthianum is a typical circular molecule 37,446 bp (isolate 89 A₂ 2-3) and 37,440 bp (isolate 83.501) in length. The difference of six nucleotides between the two genomes is the result of a deletion in the ribosomal protein S3 (rps3) gene in the 83.501 isolate. In addition, substitution of adenine for guanine within the rps3 gene in the mitochondrial genome of the 83.501 isolate was observed. Compared to the previously sequenced C. lindemuthianum mitochondrial genome, an exon no annotated in the cytochrome c oxidase I (cox1) gene and a non-conserved open reading frame (ncORF) were observed. The size of the mitochondrial genomes of the seven species of Colletotrichum was highly variable, being attributed mainly to the ncORF, ranging from one to 10 and also from introns ranging from one to 11 and which encode a total of up to nine homing endonucleases. This paper reports for the first time by means of transcriptome that then ncORFs are transcribed in Colletotrichum spp. Phylogeny data revealed that core mitochondrial genes could be used as an alternative in phylogenetic relationship studies in Colletotrichum spp. This work contributes to the genetic and biological knowledge of Colletotrichum spp., which is of great economic and scientific importance.

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.

Fungi that Infect Humans

Fungi must meet four criteria to infect humans: growth at human body temperatures, circumvention or penetration of surface barriers, lysis and absorption of tissue, and resistance to immune defenses, including elevated body temperatures. Morphogenesis between small round, detachable cells and long, connected cells is the mechanism by which fungi solve problems of locomotion around or through host barriers. Secretion of lytic enzymes, and uptake systems for the released nutrients, are necessary if a fungus is to nutritionally utilize human tissue. Last, the potent human immune system evolved in the interaction with potential fungal pathogens, so few fungi meet all four conditions for a healthy human host. Paradoxically, the advances of modern medicine have made millions of people newly susceptible to fungal infections by disrupting immune defenses. This article explores how different members of four fungal phyla use different strategies to fulfill the four criteria to infect humans: the Entomophthorales, the Mucorales, the Ascomycota, and the Basidiomycota. Unique traits confer human pathogenic potential on various important members of these phyla: pathogenic Onygenales comprising thermal dimorphs such as Histoplasma and Coccidioides; the Cryptococcus spp. that infect immunocompromised as well as healthy humans; and important pathogens of immunocompromised patients-Candida, Pneumocystis, and Aspergillus spp. Also discussed are agents of neglected tropical diseases important in global health such as mycetoma and paracoccidiomycosis and common pathogens rarely implicated in serious illness such as dermatophytes. Commensalism is considered, as well as parasitism, in shaping genomes and physiological systems of hosts and fungi during evolution.

Transcriptome analysis of the two unrelated fungal β-lactam producers Acremonium chrysogenum and Penicillium chrysogenum: Velvet-regulated genes are major targets during conventional strain improvement programs

BACKGROUND

Cephalosporins and penicillins are the most frequently used β-lactam antibiotics for the treatment of human infections worldwide. The main industrial producers of these antibiotics are Acremonium chrysogenum and Penicillium chrysogenum, two taxonomically unrelated fungi. Both were subjects of long-term strain development programs to reach economically relevant antibiotic titers. It is so far unknown, whether equivalent changes in gene expression lead to elevated antibiotic titers in production strains.

RESULTS

Using the sequence of PcbC, a key enzyme of β-lactam antibiotic biosynthesis, from eighteen different pro- and eukaryotic microorganisms, we have constructed a phylogenetic tree to demonstrate the distant relationship of both fungal producers. To address the question whether both fungi have undergone similar genetic adaptions, we have performed a comparative gene expression analysis of wild-type and production strains. We found that strain improvement is associated with the remodeling of the transcriptional landscape in both fungi. In P. chrysogenum, 748 genes showed differential expression, while 1572 genes from A. chrysogenum are differentially expressed in the industrial strain. Common in both fungi is the upregulation of genes belonging to primary and secondary metabolism, notably those involved in precursor supply for β-lactam production. Other genes not essential for β-lactam production are downregulated with a preference for those responsible for transport processes or biosynthesis of other secondary metabolites. Transcriptional regulation was shown to be an important parameter during strain improvement in different organisms. We therefore investigated deletion strains of the major transcriptional regulator velvet from both production strains. We identified 567 P. chrysogenum and 412 A. chrysogenum Velvet target genes. In both deletion strains, approximately 50% of all secondary metabolite cluster genes are differentially regulated, including β-lactam biosynthesis genes. Most importantly, 35-57% of Velvet target genes are among those that showed differential expression in both improved industrial strains.

CONCLUSIONS

The major finding of our comparative transcriptome analysis is that strain improvement programs in two unrelated fungal β-lactam antibiotic producers alter the expression of target genes of Velvet, a global regulator of secondary metabolism. From these results, we conclude that regulatory alterations are crucial contributing factors for improved β-lactam antibiotic titers during strain improvement in both fungi.

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.

From NGS assembly challenges to instability of fungal mitochondrial genomes: A case study in genome complexity

The presence of repetitive or non-unique DNA persisting over sizable regions of a eukaryotic genome can hinder the genome's successful de novo assembly from short reads: ambiguities in assigning genome locations to the non-unique subsequences can result in premature termination of contigs and thus overfragmented assemblies. Fungal mitochondrial (mtDNA) genomes are compact (typically less than 100 kb), yet often contain short non-unique sequences that can be shown to impede their successful de novo assembly in silico. Such repeats can also confuse processes in the cell in vivo. A well-studied example is ectopic (out-of-register, illegitimate) recombination associated with repeat pairs, which can lead to deletion of functionally important genes that are located between the repeats. Repeats that remain conserved over micro- or macroevolutionary timescales despite such risks may indicate functionally or structurally (e.g., for replication) important regions. This principle could form the basis of a mining strategy for accelerating discovery of function in genome sequences. We present here our screening of a sample of 11 fully sequenced fungal mitochondrial genomes by observing where exact k-mer repeats occurred several times; initial analyses motivated us to focus on 17-mers occurring more than three times. Based on the diverse repeats we observe, we propose that such screening may serve as an efficient expedient for gaining a rapid but representative first insight into the repeat landscapes of sparsely characterized mitochondrial chromosomes. Our matching of the flagged repeats to previously reported regions of interest supports the idea that systems of persisting, non-trivial repeats in genomes can often highlight features meriting further attention.

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

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

Characterization and phylogenetic analysis of the mitochondrial genome of Shiraia bambusicola reveals special features in the order of pleosporales

Shiraia bambusicola P. Henn. is a pathogenic fungus of bamboo, and its fruiting bodies are regarded as folk medicine. We determined and analyzed its complete mitochondrial DNA sequence (circular DNA molecule of 39,030 bp, G + C content of 25.19%). It contains the typical genes encoding proteins involved in electron transport and coupled oxidative phosphorylation (nad1-6 and nad4L, cob and cox1-3), one ATP synthase subunit (atp6), 4 hypothetical proteins, and two genes for large and small rRNAs (rnl and rns). There is a set of 32 tRNA genes comprising all 20 amino acids, and these genes are evenly distributed on the two strands. Phylogenetic analyses based on concatenated mitochondrial proteins indicated that S. bambusicola clustered with members of the order Pleosporales, which is in agreement with previous results. The gene arrangements of Dothideomycetes species contained three regions of gene orders partitioned in their mitochondrial genomes, including block 1 (nad6-atp6), block 2 (nad1-cox3) and block 3 (genes around rns). S. bambusicola displayed unique special features that differed from the other Pleosporales species, especially in the coding regions around rns (trnR-trnY). Moreover, a comparison of gene orders in mitochondrial genomes from Pezizomycotina revealed that although all encoded regions are located on the same strand in most Pezizomycotina mtDNAs, genes from Dothideomycetes species had different orientations, as well as diverse positions and colocalization of genes (such as cox3, cox1-cox2 and nad2-nad3); these distinctions were regarded as class-specific features. Interestingly, two incomplete copies of the atp6 gene were found on different strands of the mitogenomic DNA, a finding that has not been observed in the other analyzed fungal species. In our study, mitochondrial genomes from Dothideomycetes species were comprehensively analyzed for the first time, including many species that have not appeared in previous reports.

Eumycetoma and actinomycetoma--an update on causative agents, epidemiology, pathogenesis, diagnostics and therapy

Mycetoma is a chronic putrid infection of the cutaneous and subcutaneous tissue concerning predominantly the feet, and more rarely other body parts. Mycetoma can be caused by both fungi (eumycetoma) and bacteria (actinomycetoma). Mode of infection is an inoculation of the causative microorganism via small injuries of the skin. The clinical correlate of both forms of mycetoma is tumescence with abscesses, painless nodules, sinuses and discharge. The latter is commonly serous-purulent and contains grains (filamentous granules) which can be expressed for diagnostic purposes. Distinctive for both eumycetoma and actinomycetoma, are the formation of grains. Grains represent microcolonies of the microorganism in vivo in the vital tissue. The most successful treatment option for eumycetomas offers itraconazole in a dosage of 200 mg twice daily. This triazole antifungal is considered as 'gold standard' for eumycetomas. Alternatively, the cheaper ketoconazole was widely used, however, it was currently stopped by the FDA. Actinomycetomas should be treated by the combination of trimethoprim-sulphamethoxazole (co-trimoxazole 80/400 to 160/800 mg per day) and amikacin 15 mg/kg body weight per day. Mycetomas are neglected infections of the poor. They are more than a medical challenge. In rural areas of Africa, Asia and South America mycetomas lead to socio-economic consequences involving the affected patients, their families and the society in general.

Revision of agents of black-grain eumycetoma in the order Pleosporales

Eumycetoma is a chronic fungal infection characterised by large subcutaneous masses and the presence of sinuses discharging coloured grains. The causative agents of black-grain eumycetoma mostly belong to the orders Sordariales and Pleosporales. The aim of the present study was to clarify the phylogeny and taxonomy of pleosporalean agents, viz. Madurella grisea, Medicopsis romeroi (syn.: Pyrenochaeta romeroi), Nigrograna mackinnonii (syn. Pyrenochaeta mackinnonii), Leptosphaeriasenegalensis, L. tompkinsii, and Pseudochaetosphaeronema larense. A phylogenetic analysis based on five loci was performed: the Internal Transcribed Spacer (ITS), large (LSU) and small (SSU) subunit ribosomal RNA, the second largest RNA polymerase subunit (RPB2), and translation elongation factor 1-alpha (TEF1) gene. In addition, the morphological and physiological characteristics were determined. Three species were well-resolved at the family and genus level. Madurella grisea, L. senegalensis, and L. tompkinsii were found to belong to the family Trematospheriaceae and are reclassified as Trematosphaeria grisea comb. nov., Falciformispora senegalensis comb. nov., and F. tompkinsii comb. nov. Medicopsis romeroi and Pseudochaetosphaeronema larense were phylogenetically distant and both names are accepted. The genus Nigrograna is reduced to synonymy of Biatriospora and therefore N. mackinnonii is reclassified as B.mackinnonii comb. Nov. Mycetoma agents in Pleosporales were phylogenetically quite diverse despite their morphological similarity in the formation of pycnidia, except for the ascosporulating genus Falciformispora (formerly in Leptosphaeria). Most of the species diagnosed from human mycetoma were found to be related to waterborne or marine fungi, suggesting an association of the virulence factors with oligotrophism or halotolerance.

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.

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.

Complete mitochondrial genome of the medicinal mushroom Ganoderma lucidum

Ganoderma lucidum is one of the well-known medicinal basidiomycetes worldwide. The mitochondrion, referred to as the second genome, is an organelle found in most eukaryotic cells and participates in critical cellular functions. Elucidating the structure and function of this genome is important to understand completely the genetic contents of G. lucidum. In this study, we assembled the mitochondrial genome of G. lucidum and analyzed the differential expressions of its encoded genes across three developmental stages. The mitochondrial genome is a typical circular DNA molecule of 60,630 bp with a GC content of 26.67%. Genome annotation identified genes that encode 15 conserved proteins, 27 tRNAs, small and large rRNAs, four homing endonucleases, and two hypothetical proteins. Except for genes encoding trnW and two hypothetical proteins, all genes were located on the positive strand. For the repeat structure analysis, eight forward, two inverted, and three tandem repeats were detected. A pair of fragments with a total length around 5.5 kb was found in both the nuclear and mitochondrial genomes, which suggests the possible transfer of DNA sequences between two genomes. RNA-Seq data for samples derived from three stages, namely, mycelia, primordia, and fruiting bodies, were mapped to the mitochondrial genome and qualified. The protein-coding genes were expressed higher in mycelia or primordial stages compared with those in the fruiting bodies. The rRNA abundances were significantly higher in all three stages. Two regions were transcribed but did not contain any identified protein or tRNA genes. Furthermore, three RNA-editing sites were detected. Genome synteny analysis showed that significant genome rearrangements occurred in the mitochondrial genomes. This study provides valuable information on the gene contents of the mitochondrial genome and their differential expressions at various developmental stages of G. lucidum. The results contribute to the understanding of the functions and evolution of fungal mitochondrial DNA.

Phylogenetic findings suggest possible new habitat and routes of infection of human eumyctoma

Eumycetoma is a traumatic fungal infection in tropical and subtropical areas that may lead to severe disability. Madurella mycetomatis is one of the prevalent etiologic agents in arid Northeastern Africa. The source of infection has not been clarified. Subcutaneous inoculation from plant thorns has been hypothesized, but attempts to detect the fungus in relevant material have remained unsuccessful. The present study aims to find clues to reveal the natural habitat of Madurella species using a phylogenetic approach, i.e. by comparison of neighboring taxa with known ecology. Four species of Madurella were included in a large data set of species of Chaetomium, Chaetomidium, Thielavia, and Papulaspora (n = 128) using sequences of the universal fungal barcode gene rDNA ITS and the partial LSU gene sequence. Our study demonstrates that Madurella species are nested within the Chaetomiaceae, a family of fungi that mainly inhabit animal dung, enriched soil, and indoor environments. We hypothesize that cattle dung, ubiquitously present in rural East Africa, plays a significant role in the ecology of Madurella. If cow dung is an essential factor in inoculation by Madurella, preventative measures may involve the use of appropriate footwear in addition to restructuring of villages to reduce the frequency of contact with etiologic agents of mycetoma. On the other hand, the Chaetomiaceae possess a hidden clinical potential which needs to be explored.

Eumycetoma caused by Madurella mycetomatis is a common subcutaneous, mutilating fungal infection endemic in arid climate zones. Still there are many controversies on the route of infection, but traumatic inoculation of the subcutaneous tissue with the thorn or soil causative organism through minor skin trauma is a popular theory. This is due to the fact that, the origin and natural habitat of Madurella species, the prevalent mycetoma agents are still unknown. In order to predict the natural habitat of M. mycetomatis we investigated its phylogenetic relationships to species with known ecology. Two genes phylogeny based on LSU and ITS was performed for the species of the genus Madurella and representative genera from the family of Chaetomiaceae. Our findings confirmed that Madurella species are phylogenetically member of the family Chaetomiaceae. Members of this family are often found in dung and manure-enriched soil. We therefore suggest that animal dung, abundantly present in endemic villages, could be a possible niche for Madurella and plays an essential role in the onset of eumycetoma. This will help in understanding the origin of the disease and could be a base for future in depth study to investigate the presence of Madurella in dung from endemic areas.