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

Comparative analysis of mitochondrial genomes in Ceratocystis fimbriata complex across diverse hosts

The decline ofAcacia mangiumWilld. in Malaysia, especially in Sabah since 2010, is primarily due to Ceratocystiswilt and canker disease (CWCD) caused by theCeratocystis fimbriataEllis & Halst. complex. This study was aimed to investigate the mitochondrial genome architecture of two differentC. fimbriatacomplex isolates from Malaysia: one fromA. mangiumin Pahang (FRIM1162) and another fromEucalyptus pellitain Sarawak (FRIM1441). This research employed Next-Generation Sequencing (NGS) to contrast genomes from diverse hosts with nine additional mitochondrial sequences, identifying significant genetic diversity and mutational hotspots in the mitochondrial genome alignment. The mitochondrial genome-based phylogenetic analysis revealed a significant genetic relationship between the studied isolates and theC. fimbriatacomplex in the South American Subclade, indicating that theC. fimbriatacomplex discovered in Malaysia isC. manginecans. The comparative mitochondrial genome demonstrates the adaptability of the complex due to mobile genetic components and genomic rearrangements in the studiedfungal isolates. This research enhances our knowledge of the genetic diversity and evolutionary patterns within theC. fimbriatacomplex, aiding in a deeper understanding of fungal disease development and host adaption processes. The acquired insights are crucial for creating specific management strategies for CWCD, improving the overall understanding of fungal disease evolution and control.

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

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

Highly Reactive Group I Introns Ubiquitous in Pathogenic Fungi

Systemic fungal infections are a growing public health threat, and yet viable antifungal drug targets are limited as fungi share a similar proteome with humans. However, features of RNA metabolism and the noncoding transcriptomes in fungi are distinctive. For example, fungi harbor highly structured RNA elements that humans lack, such as self-splicing introns within key housekeeping genes in the mitochondria. However, the location and function of these mitochondrial riboregulatory elements has largely eluded characterization. Here we used an RNA-structure-based bioinformatics pipeline to identify the group I introns interrupting key mitochondrial genes in medically relevant fungi, revealing their fixation within a handful of genetic hotspots and their ubiquitous presence across divergent phylogenies of fungi, including all highest priority pathogens such as Candida albicans, Candida auris, Aspergillus fumigatus and Cryptococcus neoformans. We then biochemically characterized two representative introns from C. albicans and C. auris, demonstrating their exceptionally efficient splicing catalysis relative to previously-characterized group I introns. Indeed, the C. albicans mitochondrial intron displays extremely rapid catalytic turnover, even at ambient temperatures and physiological magnesium ion concentrations. Our results unmask a significant new set of players in the RNA metabolism of pathogenic fungi, suggesting a promising new type of antifungal drug target.

Reexamination of Aspergillus cristatus phylogeny in dark tea: Characteristics of the mitochondrial genome

To enhance our understanding of Aspergillus cristatus, an important functional microorganism, the characteristics of its mitochondrial genome were analyzed and compared with related species. The mitochondrial genome of A. cristatus was determined to be 77,649 bp in length, with 15 protein-coding regions. Notably, its length surpassed that of the other species, primarily attributable to the intron length. Gene order exhibited significant variations, with greater conservation observed in the genus Penicillium compared to Aspergillus. Phylogenetic tree analyses indicated that the genera Aspergillus and Penicillium are closely related but monophyletic. Furthermore, the phylogenetic tree constructed based on protein-coding genes effectively distinguished all strains with high branching confidence. This approach provides a robust reflection of the evolutionary relationship between A. cristatus and its related species, offering potential for the development of molecular markers suitable for Aspergillus and Penicillium.

Mitogenome Variations in a Global Population of Aspergillus fumigatus

Aspergillus fumigatus is a ubiquitous, critical priority human fungal pathogen. Despite its clinical importance, there is limited knowledge regarding the variations of the genome within mitochondria, the powerhouse organelle within eukaryotic cells. In this study, we leveraged publicly available, raw, whole genome sequence data isolates from 1939 to investigate the variations in the mitochondrial genomes of A. fumigatus. These isolates were isolated from 22 countries on six continents, as well as from outer space and from within the International Space Station. In total, our analysis revealed 39 mitochondrial single nucleotide polymorphisms (mtSNPs) within this global sample, and, together, these 39 mtSNPs grouped the 1939 isolates into 79 mitochondrial multilocus genotypes (MLGs). Among the 79 MLGs, 39 were each distributed in at least two countries and 30 were each shared by at least two continents. The two most frequent MLGs were also broadly distributed: MLG11 represented 420 isolates from 11 countries and four continents and while MLG79 represented 418 isolates from 18 countries and five continents, consistent with long-distance dispersals of mitogenomes. Our population genetic analyses of the mtSNPs revealed limited differentiation among continental populations, but highly variable genetic differences among national populations, largely due to localized clonal expansions of different MLGs. Phylogenetic analysis and Discriminant Analysis of Principal Components of mtSNPs suggested the presence of at least three mitogenome clusters. Linkage disequilibrium, Index of Association, and phylogenetic incompatibility analyses collectively suggested evidence for mitogenome recombination in natural populations of A. fumigatus. In addition, sequence read depth analyses revealed an average ratio of ~20 mitogenomes per nuclear genome in this global population, but the ratios varied among strains within and between certain geographic populations. Together, our results suggest evidence for organelle dynamics, genetic differentiation, recombination, and both widespread and localized clonal expansion of the mitogenomes in the global A. fumigatus population.

The First Five Mitochondrial Genomes for the Family Nidulariaceae Reveal Novel Gene Rearrangements, Intron Dynamics, and Phylogeny of Agaricales

The family Nidulariaceae, consisting of five genera including Cyathus, is a unique group of mushrooms commonly referred to as bird's nest fungi due to their striking resemblance to bird's nests. These mushrooms are considered medicinal mushrooms in Chinese medicine and have received attention in recent years for their anti-neurodegenerative properties. However, despite the interest in these mushrooms, very little is known about their mitochondrial genomes (mitogenomes). This study is the first comprehensive investigation of the mitogenomes of five Nidulariaceae species with circular genome structures ranging in size from 114,236 bp to 129,263 bp. Comparative analyses based on gene content, gene length, tRNA, and codon usage indicate convergence within the family Nidulariaceae and heterogeneity within the order Agaricales. Phylogenetic analysis based on a combined mitochondrial conserved protein dataset provides a well-supported phylogenetic tree for the Basidiomycetes, which clearly demonstrates the evolutionary relationships between Nidulariaceae and other members of Agaricales. Furthermore, phylogenetic inferences based on four different gene sets reveal the stability and proximity of evolutionary relationships within Agaricales. These results reveal the uniqueness of the family Nidulariaceae and its similarity to other members of Agaricales; provide valuable insights into the origin, evolution, and genetics of Nidulariaceae species; and enrich the fungal mitogenome resource. This study will help to expand the knowledge and understanding of the mitogenomes in mushrooms.

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

Introduction

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

Methods

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

Results and Discussion

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

A highly quality genome sequence of Penicillium oxalicum species isolated from the root of Ixora chinensis in Vietnam

Penicillium oxalicum has been reported as a multienzyme-producing fungus and is widely used in industry due to great potential for cellulase release. Until now, there are only 10 available genome assemblies of P. oxalicum species deposited in the GenBank database. In this study, the genome of the I1R1 strain isolated from the root of Ixora chinensis was completely sequenced by Pacbio Sequel sequencing technology, assembled into 8 chromosomes with the genome size of 30.8 Mb, as well as a mitogenome of 26 kb. The structural and functional analyses of the I1R1 genome revealed gene model annotations encoding an enzyme set involved in significant metabolic processes, along with cytochrome P450s and secondary metabolite biosynthesis. The comparative analysis of the P. oxalicum species based on orthology and gene family duplications indicated their large and closed pan-genome of 9,500 orthologous groups. This is valuable data for future phylogenetic and population genomics studies.

Evolution and codon usage bias of mitochondrial and nuclear genomes in Aspergillus section Flavi

The fungal genus Aspergillus contains a diversity of species divided into taxonomic sections of closely related species. Section Flavi contains 33 species, many of industrial, agricultural, or medical relevance. Here, we analyze the mitochondrial genomes (mitogenomes) of 20 Flavi species-including 18 newly assembled mitogenomes-and compare their evolutionary history and codon usage bias patterns to their nuclear counterparts. Codon usage bias refers to variable frequencies of synonymous codons in coding DNA and is shaped by a balance of neutral processes and natural selection. All mitogenomes were circular DNA molecules with highly conserved gene content and order. As expected, genomic content, including GC content, and genome size differed greatly between mitochondrial and nuclear genomes. Phylogenetic analysis based on 14 concatenated mitochondrial genes predicted evolutionary relationships largely consistent with those predicted by a phylogeny constructed from 2,422 nuclear genes. Comparing similarities in interspecies patterns of codon usage bias between mitochondrial and nuclear genomes showed that species grouped differently by patterns of codon usage bias depending on whether analyses were performed using mitochondrial or nuclear relative synonymous usage values. We found that patterns of codon usage bias at gene level are more similar between mitogenomes of different species than the mitogenome and nuclear genome of the same species. Finally, we inferred that, although most genes-both nuclear and mitochondrial-deviated from the neutral expectation for codon usage, mitogenomes were not under translational selection while nuclear genomes were under moderate translational selection. These results contribute to the study of mitochondrial genome evolution in filamentous fungi.

Telomere-to-telomere genome sequence of the model mould pathogen Aspergillus fumigatus

The pathogenic fungus Aspergillus fumigatus is a major etiological agent of fungal invasive and chronic diseases affecting tens of millions of individuals worldwide. Draft genome sequences of two clinical isolates (Af293 and A1163) are commonly used as reference genomes for analyses of clinical and environmental strains. However, the reference sequences lack coverage of centromeres, an accurate sequence for ribosomal repeats, and a comprehensive annotation of chromosomal rearrangements such as translocations and inversions. Here, we used PacBio Single Molecule Real-Time (SMRT), Oxford Nanopore and Illumina HiSeq sequencing for de novo genome assembly and polishing of two laboratory reference strains of A. fumigatus, CEA10 (parental isolate of A1163) and its descendant A1160. We generated full length chromosome assemblies and a comprehensive telomere-to-telomere coverage for CEA10 and near complete assembly of A1160 including ribosomal repeats and the sequences of centromeres, which we discovered to be composed of long transposon elements. We envision these high-quality reference genomes will become fundamental resources to study A. fumigatus biology, pathogenicity and virulence, and to discover more effective treatments against diseases caused by this fungus.

The fungus Aspergillus fumigatus causes invasive and chronic diseases worldwide. Here, Bowyer et al. use long-read and short-read sequencing to generate complete chromosome assemblies and telomere-to-telomere coverage for two isolates, thus providing high-quality reference genomes as fundamental resources to study this pathogen.

Towards Understanding the Function of Aegerolysins

Aegerolysins are remarkable proteins. They are distributed over the tree of life, being relatively widespread in bacteria and fungi, but also present in some insects, plants, protozoa, and viruses. Despite their abundance in cells of certain developmental stages and their presence in secretomes, only a few aegerolysins have been studied in detail. Their function, in particular, is intriguing. Here, we summarize previously published findings on the distribution, molecular interactions, and function of these versatile aegerolysins. They have very diverse protein sequences but a common fold. The machine learning approach of the AlphaFold2 algorithm, which incorporates physical and biological knowledge of protein structures and multisequence alignments, provides us new insights into the aegerolysins and their pore-forming partners, complemented by additional genomic support. We hypothesize that aegerolysins are involved in the mechanisms of competitive exclusion in the niche.

Intraspecific comparison of mitochondrial genomes reveals the evolution in medicinal fungus Ganoderma lingzhi

Several mitogenomes of the genus Ganoderma have been assembled, but intraspecific comparisons of mitogenomes in Ganoderma lingzhi have not been reported. In this study, 19 G. lingzhi mitogenomes were assembled and analyzed combined with three mitogenomes of G. lingzhi from GenBank in term of the characteristics, evolution, and phylogeny. The results showed that the mitogenomes of the G. lingzhi strains are closed circular ranging from 49.23 kb to 68.37 kb. The genetic distance, selective pressure, and base variation indicate that the 14 common protein coding genes were highly conserved. The differences in introns, open reading frames, and repetitive sequences in the mitogenome were the main factors leaded to the variations in mitogenome. The introns were horizontally transferred in mitogenomes, and the differences between introns in the same insertion, which were primarily caused by the repetitive sequence, showed that the introns may be under degeneration. Besides, the frequent insertion and deletion of introns showed an evolutionary rate faster than protein coding genes. Phylogenetic analysis showed that the G. lingzhi strains gathered with high support, and those with the same intron distribution law had closer clustering relationships.

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

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

The importance of complete and high-quality genome sequences in Aspergillus niger research

The possibility to sequence the entire genome of an organism revolutionized the fields of biology and biotechnology. The first genome sequence of the important filamentous fungus Aspergillus niger was obtained in 2007, 11 years after the release of the first eukaryotic genome sequence. From that moment, genomics of A. niger has seen major progresses, facilitated by the advances in the sequencing technologies and in the methodologies for gene function prediction. However, there are still challenges to face when trying to obtain complete genomes, equipped with all the repetitive sequences that they contain and without omitting the mitochondrial sequences. The aim of this perspective article is to discuss the current status of A. niger genomics and draw attention to the open challenges that the fungal community should address to move research of this important fungus forward.

Phylogenetic and ecological reevaluation of the order Onygenales

The order Onygenales is classified in the class Eurotiomycetes of the subphylum Pezizomycotina. Families in this order have classically been isolated from soil and dung, and two lineages contain causative agents of superficial, cutaneous and systemic infections in mammals. The ecology and habitat choices of the species are driven mainly by the keratin and cellulose degradation abilities. The present study aimed to investigate whether the ecological trends of the members of Onygenales can be interpreted in an evolutionary sense, linking phylogenetic parameters with habitat preferences, to achieve polyphasic definitions of the main taxonomic groups. Evolutionary processes were estimated by multiple gene genealogies and divergence time analysis. Previously described families, namely, Arthrodermataceae, Ajellomycetaceae, Ascosphaeraceae, Eremascaceae, Gymnoascaceae, Onygenaceae and Spiromastigoidaceae, were accepted in Onygenales, and two new families, Malbrancheaceae and Neogymnomycetaceae, were introduced. A number of species could not be assigned to any of the defined families. Our study provides a revised overview of the main lines of taxonomy of Onygenales, supported by multilocus analyses of ITS, LSU, TUB, TEF1, TEF3, RPB1, RPB2, and ribosomal protein 60S L10 (L1) (RP60S) sequences, combined with available data on ecology, physiology, morphology, and genomics.

Quantitative proteomic analysis of Neosartorya pseudofischeri ascospores subjected to heat treatment

Neosartorya pseudofischeri, a heat- resistant fungus, was isolated from contaminated blueberry juice in our laboratory. To obtain a better understanding of the cellular process of heat- resistant fungus, a quantitative proteomic investigation employing stable isotope labeling by peptide demethylation was conducted on changes in intracellular proteins of N. pseudofischeri ascospores subjected to heat treatment at 93 °C for 0, 1 or 8 min. In total, 811 proteins were identified and quantified. Using the normalized ratio of protein abundance, proteins that changed more than two- fold after heat treatment were identified as significantly increased or decreased proteins and grouped into four clusters based on their quantitative changes. Decreased proteins were found mainly involved in the central carbon metabolism, heat stress responses, reactive oxygen intermediates elimination and translation events. A group of proteins in relation to toxicant degradation and antibiotic neutralization linking to environmental adaptability and tolerance of heat- resistant, was also identified. These findings provide insights into protein changes of N. pseudofischeri ascospores and lay foundations for further investigations on heat- resistant molds using targeted quantitative approaches. SIGNIFICANCE OF THE RESEARCH: Heat- resistant fungi can survive pasteurization processes, and subsequently germinate and grow to cause spoilage of food products, leading to significant economic losses for the food industry and potential health risk for the consumers. To obtain a better understanding of the cellular process of heat- resistant fungi, a quantitative proteomic approach employing stable isotope labeling by peptide demethylation has been used to investigate the intracellular protein changes of N. pseudofischeri ascospores isolated from blueberry juice and subjected to heat treatment. 150 significantly changed proteins were grouped into four clusters based on their quantitative changes. The significant decrease in protein abundance in response to heat treatment revealed possible mechanism that N. pseudofischeri ascospores could survive the heat treatment. This is the first proteomic profile report for N. pseudofischeri. These findings provide insights into protein changes of N. pseudofischeri and lay foundations for further investigations on heat- resistant fungi using targeted quantitative approaches to evaluate the efficiency of thermal treatment for processed food products.

Global Characterization of Fungal Mitogenomes: New Insights on Genomic Diversity and Dynamism of Coding Genes and Accessory Elements

Fungi comprise a great diversity of species with distinct ecological functions and lifestyles. Similar to other eukaryotes, fungi rely on interactions with prokaryotes and one of the most important symbiotic events was the acquisition of mitochondria. Mitochondria are organelles found in eukaryotic cells whose main function is to generate energy through aerobic respiration. Mitogenomes (mtDNAs) are double-stranded circular or linear DNA from mitochondria that may contain core genes and accessory elements that can be replicated, transcribed, and independently translated from the nuclear genome. Despite their importance, investigative studies on the diversity of fungal mitogenomes are scarce. Herein, we have evaluated 788 curated fungal mitogenomes available at NCBI database to assess discrepancies and similarities among them and to better understand the mechanisms involved in fungal mtDNAs variability. From a total of 12 fungal phyla, four do not have any representative with available mitogenomes, which highlights the underrepresentation of some groups in the current available data. We selected representative and non-redundant mitogenomes based on the threshold of 90% similarity, eliminating 81 mtDNAs. Comparative analyses revealed considerable size variability of mtDNAs with a difference of up to 260 kb in length. Furthermore, variation in mitogenome length and genomic composition are generally related to the number and length of accessory elements (introns, HEGs, and uORFs). We identified an overall average of 8.0 (0–39) introns, 8.0 (0–100) HEGs, and 8.2 (0–102) uORFs per genome, with high variation among phyla. Even though the length of the core protein-coding genes is considerably conserved, approximately 36.3% of the mitogenomes evaluated have at least one of the 14 core coding genes absent. Also, our results revealed that there is not even a single gene shared among all mitogenomes. Other unusual genes in mitogenomes were also detected in many mitogenomes, such as dpo and rpo, and displayed diverse evolutionary histories. Altogether, the results presented in this study suggest that fungal mitogenomes are diverse, contain accessory elements and are absent of a conserved gene that can be used for the taxonomic classification of the Kingdom Fungi.

Complete Mitochondrial Genome Sequences of Nine Aspergillus flavus Strains

Nuclear genome sequences incompletely characterize the genomic content and thus the genetic diversity of fungal species. Here, we present the complete mitochondrial genome sequences of nine Aspergillus flavus strains, providing useful information for inter- and intraspecific analyses.

Aspergillus fumigatus versus Genus Aspergillus: Conservation, Adaptive Evolution and Specific Virulence Genes

Aspergillus is an important fungal genus containing economically important species, as well as pathogenic species of animals and plants. Using eighteen fungal species of the genus Aspergillus, we conducted a comprehensive investigation of conserved genes and their evolution. This also allows us to investigate the selection pressure driving the adaptive evolution in the pathogenic species A. fumigatus. Among single-copy orthologs (SCOs) for A. fumigatus and the closely related species A. fischeri, we identified 122 versus 50 positively selected genes (PSGs), respectively. Moreover, twenty conserved genes of unknown function were established to be positively selected and thus important for adaption. A. fumigatus PSGs interacting with human host proteins show over-representation of adaptive, symbiosis-related, immunomodulatory and virulence-related pathways, such as the TGF-β pathway, insulin receptor signaling, IL1 pathway and interfering with phagosomal GTPase signaling. Additionally, among the virulence factor coding genes, secretory and membrane protein-coding genes in multi-copy gene families, 212 genes underwent positive selection and also suggest increased adaptation, such as fungal immune evasion mechanisms (aspf2), siderophore biosynthesis (sidD), fumarylalanine production (sidE), stress tolerance (atfA) and thermotolerance (sodA). These genes presumably contribute to host adaptation strategies. Genes for the biosynthesis of gliotoxin are shared among all the close relatives of A. fumigatus as an ancient defense mechanism. Positive selection plays a crucial role in the adaptive evolution of A. fumigatus. The genome-wide profile of PSGs provides valuable targets for further research on the mechanisms of immune evasion, antimycotic targeting and understanding fundamental virulence processes.

Genome sequencing of the neotype strain CBS 554.65 reveals the MAT1-2 locus of Aspergillus niger

Background

Aspergillus niger is a ubiquitous filamentous fungus widely employed as a cell factory thanks to its abilities to produce a wide range of organic acids and enzymes. Its genome was one of the first Aspergillus genomes to be sequenced in 2007, due to its economic importance and its role as model organism to study fungal fermentation. Nowadays, the genome sequences of more than 20 A. niger strains are available. These, however, do not include the neotype strain CBS 554.65.

Results

The genome of CBS 554.65 was sequenced with PacBio. A high-quality nuclear genome sequence consisting of 17 contigs with a N50 value of 4.07 Mbp was obtained. The assembly covered all the 8 centromeric regions of the chromosomes. In addition, a complete circular mitochondrial DNA assembly was obtained. Bioinformatic analyses revealed the presence of a MAT1-2-1 gene in this genome, contrary to the most commonly used A. niger strains, such as ATCC 1015 and CBS 513.88, which contain a MAT1-1-1 gene. A nucleotide alignment showed a different orientation of the MAT1–1 locus of ATCC 1015 compared to the MAT1–2 locus of CBS 554.65, relative to conserved genes flanking the MAT locus. Within 24 newly sequenced isolates of A. niger half of them had a MAT1–1 locus and the other half a MAT1–2 locus. The genomic organization of the MAT1–2 locus in CBS 554.65 is similar to other Aspergillus species. In contrast, the region comprising the MAT1–1 locus is flipped in all sequenced strains of A. niger.

Conclusions

This study, besides providing a high-quality genome sequence of an important A. niger strain, suggests the occurrence of genetic flipping or switching events at the MAT1–1 locus of A. niger. These results provide new insights in the mating system of A. niger and could contribute to the investigation and potential discovery of sexuality in this species long thought to be asexual.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07990-8.

Uncovering Diagnostic Value of Mitogenome for Identification of Cryptic Species Fusarium graminearum Sensu Stricto

Fungal complexes are often composed of morphologically nearly indistinguishable species with high genetic similarity. However, despite their close relationship, they can exhibit distinct phenotypic differences in pathogenicity and production of mycotoxins. Many plant pathogenic and toxigenic fungi have been shown to consist of such cryptic species. Identification of cryptic species in economically important pathogens has added value in epidemiologic studies and provides opportunities for better control. Analysis of mitochondrial genomes or mitogenomics opens up dimensions for improved diagnostics of fungi, especially when efficient recovery of DNA is problematic. In comparison to nuclear DNA, mitochondrial DNA (mtDNA) can be amplified with improved efficacy due to its multi-copy nature. However, to date, only a few studies have demonstrated the usefulness of mtDNA for identification of cryptic species within fungal complexes. In this study, we explored the value of mtDNA for identification of one of the most important cereal pathogens Fusarium graminearum sensu stricto (F.g.). We found that homing endonucleases (HEGs), which are widely distributed in mitogenomes of fungi, display small indel polymorphism, proven to be potentially species specific. The resulting small differences in their lengths may facilitate further differentiation of F.g. from the other cryptic species belonging to F. graminearum species complex. We also explored the value of SNP analysis of the mitogenome for typing F.g. The success in identifying F.g. strains was estimated at 96%, making this tool an attractive complement to other techniques for identification of F.g.

Chromosome assembled and annotated genome sequence of Aspergillus flavus NRRL 3357

Aspergillus flavus is an opportunistic pathogen of crops, including peanuts and maize, and is the second leading cause of aspergillosis in immunocompromised patients. A. flavus is also a major producer of the mycotoxin, aflatoxin, a potent carcinogen, which results in significant crop losses annually. The A. flavus isolate NRRL 3357 was originally isolated from peanut and has been used as a model organism for understanding the regulation and production of secondary metabolites, such as aflatoxin. A draft genome of NRRL 3357 was previously constructed, enabling the development of molecular tools and for understanding population biology of this particular species. Here, we describe an updated, near complete, telomere-to-telomere assembly and re-annotation of the eight chromosomes of A. flavus NRRL 3357 genome, accomplished via long-read PacBio and Oxford Nanopore technologies combined with Illumina short-read sequencing. A total of 13,715 protein-coding genes were predicted. Using RNA-seq data, a significant improvement was achieved in predicted 5’ and 3’ untranslated regions, which were incorporated into the new gene models.

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 Comparative Analyses of the Complete Mitochondrial Genomes of Fungal Endosymbionts in Sogatella furcifera, White-Backed Planthoppers

Sogatella furcifera Horvath, commonly known as the white-backed planthoppers (WBPH), is an important pest in East Asian rice fields. Fungal endosymbiosis is widespread among planthoppers in the infraorder Fulgoromorpha and suborder Auchenorrhyncha. We successfully obtained complete mitogenome of five WBPH fungal endosymbionts, belonging to the Ophiocordycipitaceae family, from next-generation sequencing (NGS) reads obtained from S. furcifera samples. These five mitogenomes range in length from 55,390 bp to 55,406 bp, which is shorter than the mitogenome of the fungal endosymbiont found in Ricania speculum, black planthoppers. Twenty-eight protein-coding genes (PCGs), 12 tRNAs, and 2 rRNAs were found in the mitogenomes. Two single-nucleotide polymorphisms, two insertions, and three deletions were identified among the five mitogenomes, which were fewer in number than those of four species of Ophiocordycipitaceae, Ophiocordyceps sinensis, Hirsutella thompsonii, Hirsutella rhossiliensis, and Tolypocladium inflatum. Noticeably short lengths (up to 18 bp) of simple sequence repeats were identified in the five WBPH fungal endosymbiont mitogenomes. Phylogenetic analysis based on conserved PCGs across 25 Ophiocordycipitaceae mitogenomes revealed that the five mitogenomes were clustered with that of R. speculum, forming an independent clade. In addition to providing the full mitogenome sequences, obtaining complete mitogenomes of WBPH endosymbionts can provide insights into their phylogenetic positions without needing to isolate the mtDNA from the host. This advantage is of value to future studies involving fungal endosymbiont mitogenomes.

Recent and Ongoing Horizontal Transfer of Mitochondrial Introns Between Two Fungal Tree Pathogens

Two recently introduced fungal plant pathogens (Ceratocystis lukuohia and Ceratocystis huliohia) are responsible for Rapid ‘ōhi‘a Death (ROD) in Hawai‘i. Despite being sexually incompatible, the two pathogens often co-occur in diseased ‘ōhi‘a sapwood, where genetic interaction is possible. We sequenced and annotated 33 mitochondrial genomes of the two pathogens and related species, and investigated 35 total Ceratocystis mitogenomes. Ten mtDNA regions [one group I intron, seven group II introns, and two autonomous homing endonuclease (HE) genes] were heterogeneously present in C. lukuohia mitogenomes, which were otherwise identical. Molecular surveys with specific primers showed that the 10 regions had uneven geographic distribution amongst populations of C. lukuohia. Conversely, identical orthologs of each region were present in every studied isolate of C. huliohia regardless of geographical origin. Close relatives of C. lukuohia lacked or, rarely, had few and dissimilar orthologs of the 10 regions, whereas most relatives of C. huliohia had identical or nearly identical orthologs. Each region included or worked in tandem with HE genes or reverse transcriptase/maturases that could facilitate interspecific horizontal transfers from intron-minus to intron-plus alleles. These results suggest that the 10 regions originated in C. huliohia and are actively moving to populations of C. lukuohia, perhaps through transient cytoplasmic contact of hyphal tips (anastomosis) in the wound surface of ‘ōhi‘a trees. Such contact would allow for the transfer of mitochondria followed by mitochondrial fusion or cytoplasmic exchange of intron intermediaries, which suggests that further genomic interaction may also exist between the two pathogens.

Pan-Mitogenomics Approach Discovers Diversity and Dynamism in the Prominent Brown Rot Fungal Pathogens

Monilinia fructicola and Monilinia laxa species are the most destructive and economically devastating fungal plant pathogens causing brown rot disease on stone and pome fruits worldwide. Mitochondrial genomes (mitogenomes) play critical roles influencing the mechanisms and directions of the evolution of fungal pathogens. The pan-mitogenomics approach predicts core and accessory regions of the mitochondrial genomes and explains the gain or loss of variation within and between species. The present study is a fungal pan-mitogenome of M. fructicola (N = 8) and M. laxa (N = 8) species. The completely sequenced and annotated mitogenomes showed high variability in size within and between the species. The mitogenomes of M. laxa were larger, ranging from 178,351 to 179,780bp, than the mitogenomes of M. fructicola, ranging from 158,607 to 167,838bp. However, size variation within the species showed that M. fructicola isolates were more variable in the size range than M. laxa isolates. All the mitogenomes included conserved mitochondrial genes, as well as variable regions including different mobile introns encoding homing endonucleases or maturase, non-coding introns, and repetitive elements. The linear model analysis supported the hypothesis that the mitogenome size expansion is due to presence of variable (accessory) regions. Gene synteny was mostly conserved among all samples, with the exception for order of the rps3 in the mitogenome of one isolate. The mitogenomes presented AT richness; however, A/T and G/C skew varied among the mitochondrial genes. The purifying selection was detected in almost all the protein-coding genes (PCGs) between the species. However, cytochrome b was the only gene showing a positive selection signal among the total samples. Combined datasets of amino acid sequences of 14 core mitochondrial PCGs and rps3 obtained from this study together with published mitochondrial genome sequences from some other species from Heliotales were used to infer a maximum likelihood (ML) phylogenetic tree. ML tree indicated that both Monilinia species highly diverged from each other as well as some other fungal species from the same order. Mitogenomes harbor much information about the evolution of fungal plant pathogens, which could be useful to predict pathogenic life strategies.

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 of a Plant Fungal Pathogen Pseudocercospora fijiensis (Mycosphaerellaceae), Comparative Analysis and Diversification Times of the Sigatoka Disease Complex Using Fossil Calibrated Phylogenies

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

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

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

Mitogenomes of Two Phallus Mushroom Species Reveal Gene Rearrangement, Intron Dynamics, and Basidiomycete Phylogeny

Phallus indusiatus and Phallus echinovolvatus are edible bamboo mushrooms with pharmacological properties. We sequenced, assembled, annotated, and compared the mitogenomes of these species. Both mitogenomes were composed of circular DNA molecules, with sizes of 89,139 and 50,098 bp, respectively. Introns were the most important factor in mitogenome size variation within the genus Phallus. Phallus indusiatus, P. echinovolvatus, and Turbinellus floccosus in the subclass Phallomycetidae have conservative gene arrangements. Large-scale gene rearrangements were observed in species representing 42 different genera of Basidiomycetes. A variety of intron position classes were found in the 44 Basidiomycete species analyzed. A novel group II intron from the P. indusiatus mitogenome was compared with other fungus species containing the same intron, and we demonstrated that the insertion sites of the intron had a base preference. Phylogenetic analyses based on combined gene datasets yielded well-supported Bayesian posterior probability (BPP = 1) topologies. This indicated that mitochondrial genes are reliable molecular markers for analyzing the phylogenetic relationships of the Basidiomycetes. This is the first study of the mitogenome of the genus Phallus, and it increases our understanding of the population genetics and evolution of bamboo mushrooms and related species.

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

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

Complete mitochondrial genome sequence of Aspergillus flavus SRRC1009: insight of intraspecific variations on A. flavus mitochondrial genomes

The mitogenome of Aspergillus flavus SRRC1009 was sequenced to investigate intraspecific variations on mitochondrial genomes of A. flavus. It shows 29,202 bp with a typical configuration of Aspergillus mitogenome. Sixteen SNPs and 22 INDELs and 17 SNPs and 27 INDELs were identified against AflaGuard® and JQ355000, respectively. Phylogenetic trees present in the three A. flavus mitochondrial genomes were clustered with A. oryzae mitochondrial genome in one clade.

First Record of the Complete Mitochondrial Genome of a Saprotrophic and Opportunistic Human Pathogenic Fungus, Scopulariopsis brevicaulis

Scopulariopsis brevicaulis is a widely distributed soil fungus known as a common saprotroph of biodegradation. It is also an opportunistic human pathogen that can produce various secondary metabolites. Here, we report the first complete mitochondrial genome sequence of S. brevicaulis isolated from air in South Korea. Total length of the mitochondrial genome is 28,829 bp and encoded 42 genes (15 protein-coding genes, 2 rRNAs, and 25 tRNAs). Nucleotide sequence of coding region takes over 26.2%, and overall GC content is 27.6%. Phylogenetic trees present that S. brevicaulis is clustered with Lomentospora prolificans with presenting various mitochondrial genome length.

First characterization of the complete mitochondrial genome of fungal plant-pathogen Monilinia laxa which represents the mobile intron rich structure

Monilinia laxa is an important fungal plant pathogen causing brown rot on many stone and pome fruits worldwide. Mitochondrial genome (mitogenome) plays a critical role in evolutionary biology of the organisms. This study aimed to characterize the complete mitogenome of M. laxa by using next-generation sequencing and approaches of de novo assembly and annotation. The total length of the mitogenome of M. laxa was 178,357 bp, and its structure was circular. GC content of the mitogenome was 30.1%. Annotation of the mitogenome presented 2 ribosomal RNA (rRNA) genes, 32 transfer RNA genes (tRNA), 1 gene encoding mitochondrial ribosomal protein S3, 14 protein-coding genes and 15 open reading frame encoding hypothetical proteins. Moreover, the group I mobile introns encoding homing endonucleases including LAGLIDADG and GIY-YIG families were found both within coding regions (genic) and intergenic regions of the mitogenome, indicating an enlarged size and a dynamic structure of the mitogenome. Furthermore, a comparative mitogenomic analysis was performed between M. laxa and the three closely related fungal phytopathogen species (Botryotinia fuckeliana, Sclerotinia sclerotiorum and, S. borealis). Due to the number and distribution of introns, the large extent of structural rearrangements and diverse mitogenome sizes were detected among the species investigated. Monilinia laxa presented the highest number of homing endonucleases among the fungal species considered in the analyses. This study is the first to report a detailed annotation of the mitogenome of an isolate of M. laxa, providing a solid basis for further investigations of mitogenome variations for the other Monilinia pathogens causing brown rot disease.

Genetic Underpinnings of Host Manipulation by Ophiocordyceps as Revealed by Comparative Transcriptomics

Ant-infecting Ophiocordyceps fungi are globally distributed, host manipulating, specialist parasites that drive aberrant behaviors in infected ants, at a lethal cost to the host. An apparent increase in activity and wandering behaviors precedes a final summiting and biting behavior onto vegetation, which positions the manipulated ant in a site beneficial for fungal growth and transmission. We investigated the genetic underpinnings of host manipulation by: (i) producing a high-quality hybrid assembly and annotation of the Ophiocordyceps camponoti-floridani genome, (ii) conducting laboratory infections coupled with RNAseq of O. camponoti-floridani and its host, Camponotus floridanus, and (iii) comparing these data to RNAseq data of Ophiocordyceps kimflemingiae and Camponotus castaneus as a powerful method to identify gene expression patterns that suggest shared behavioral manipulation mechanisms across Ophiocordyceps-ant species interactions. We propose differentially expressed genes tied to ant neurobiology, odor response, circadian rhythms, and foraging behavior may result by activity of putative fungal effectors such as enterotoxins, aflatrem, and mechanisms disrupting feeding behaviors in the ant.

Sequencing and phylogenetic analysis of mitochondrial genome of Aspergillus cristatus

Aspergillus cristatus are the dominantly present microorganisms in dark tea. The whole mitochondrial genome sequence of A. cristatus was sequenced and reported in this study. The mitochondrial genome in A. cristatushas a full length of 77,649 bp, which is reported to be the longest among the mitochondrial genomes of Aspergillus species. The basesincluding A (34.14％), T (37.64％), C (15.61％) and G (12.61％) are found in their genome. A total of 42 genes (15 protein-coding genes, lrRNA/srRNA and 25 tRNAs) are encoded by the mitochondrial genome of this fungus. Phylogenetic analysis showed a closest relationship betweenA. pseudoglaucusand the taxonomic status of A. cristatus.

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.

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.

Aspergillus niger citrate exporter revealed by comparison of two alternative citrate producing conditions

Currently, there is no consensus regarding the mechanism underlying Aspergillus niger citrate biosynthesis and secretion. We hypothesise that depending on the experimental setup, extracellular citrate accumulation can have fundamentally different underlying transcriptomic landscapes. We show that varying the amount and type of supplement of an arginine auxotrophic A. niger strain results in transcriptional down-regulation of citrate metabolising enzymes in the condition in which more citrate is accumulated extracellularly. This contrasts with the transcriptional adaptations when increased citrate production is triggered by iron limitation. By combining gene expression data obtained from these two very distinct experimental setups with hidden Markov models and transporter homology approaches, we were able to compile a shortlist of the most likely citrate transporter candidates. Two candidates (An17g01710 and An09g06720m.01) were heterologously expressed in the yeast Saccharomyces cerevisiae, and one of the resultant mutants showed the ability to secrete citrate. Our findings provide steps in untangling the complex interplay of different mechanisms underlying A. niger citrate accumulation, and we demonstrate how a comparative transcriptomics approach complemented with further bioinformatics analyses can be used to pinpoint a fungal citrate exporter.

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

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

Promising Perspectives for Detection, Identification, and Quantification of Plant Pathogenic Fungi and Oomycetes through Targeting Mitochondrial DNA

Fungi and oomycetes encompass many pathogens affecting crops worldwide. Their effective control requires screening pathogens across the local and international trade networks along with the monitoring of pathogen inocula in the field. Fundamentals to all of these concerns are their efficient detection, identification, and quantification. The use of molecular markers showed the best promise in the field of plant pathogen diagnostics. However, despite the unquestionable benefits of DNA-based methods, two significant limitations are associated with their use. The first limitation concerns the insufficient level of sensitivity due to the very low and uneven distribution of pathogens in plant material. The second limitation pertains to the inability of widely used diagnostic assays to detect cryptic species. Targeting mtDNA appears to provide a solution to these challenges. Its high copy number in microbial cells makes mtDNA an attractive target for developing highly sensitive assays. In addition, previous studies on different pathogen taxa indicated that mitogenome sequence variation could improve cryptic species delimitation accuracy. This review sheds light on the potential application of mtDNA for pathogen diagnostics. This paper covers a brief description of qPCR and DNA barcoding as two major strategies enabling the diagnostics of plant pathogenic fungi and oomycetes. Both strategies are discussed along with the potential use of mtDNA, including their strengths and weaknesses.

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

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

Genome-scale data resolve ancestral rock-inhabiting lifestyle in Dothideomycetes (Ascomycota)

Dothideomycetes is the most diverse fungal class in Ascomycota and includes species with a wide range of lifestyles. Previous multilocus studies have investigated the taxonomic and evolutionary relationships of these taxa but often failed to resolve early diverging nodes and frequently generated inconsistent placements of some clades. Here, we use a phylogenomic approach to resolve relationships in Dothideomycetes, focusing on two genera of melanized, extremotolerant rock-inhabiting fungi, Lichenothelia and Saxomyces, that have been suggested to be early diverging lineages. We assembled phylogenomic datasets from newly sequenced (4) and previously available genomes (238) of 242 taxa. We explored the influence of tree inference methods, supermatrix vs. coalescent-based species tree, and the impact of varying amounts of genomic data. Overall, our phylogenetic reconstructions provide consistent and well-supported topologies for Dothideomycetes, recovering Lichenothelia and Saxomyces among the earliest diverging lineages in the class. In addition, many of the major lineages within Dothideomycetes are recovered as monophyletic, and the phylogenomic approach implemented strongly supports their relationships. Ancestral character state reconstruction suggest that the rock-inhabiting lifestyle is ancestral within the class.

Mitochondrial genome in Hypsizygus marmoreus and its evolution in Dikarya

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

Complete mitochondrial genome sequence of a xerophilic fungus, Aspergillus pseudoglaucus

Aspergillus pseudoglaucus is a xerophilic filamentous fungus which can produce various secondary metabolites. Here, we reported the complete mitochondrial genome sequence of A. pseudoglaucus isolated from Meju, a soybean brick in Korea. Its mitochondrial genome was successfully assembled from raw reads sequenced using MiSeq by Velvet and SOAPGapCloser. Total length of the mitochondrial genome is 53,882 bp, which is third longest among known Aspergillus mitochondrial genomes and encoded 58 genes (30 protein-coding genes including hypothetical ORFs, two rRNAs, and 26 tRNAs). Nucleotide sequence of coding regions takes over 66.6% and overall GC content is 27.8%. Phylogenetic trees present that A. pseudoglaucus is located outside of section Nidulantes. Additional researches will be required for clarifying phylogenetic position of section Aspergillus.

Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species

Background

Species in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time.

Results

Here we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes.

Conclusions

Our study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5732-z) contains supplementary material, which is available to authorized users.