Comparison of the Mitochondrial Genome Sequences of Six Annulohypoxylon stygium Isolates Suggests Short Fragment Insertions as a Potential Factor Leading to Larger Genomic Size

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

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

Comparative analyses of Pleurotus pulmonarius mitochondrial genomes reveal two major lineages of mini oyster mushroom cultivars

Pleurotus pulmonarius, commonly known as the mini oyster mushroom, is highly esteemed for its crisp texture and umami flavor. Limited genetic diversity among P. pulmonarius cultivars raises concerns regarding its sustainable industrial production. To delve into the maternal genetic diversity of the principal P. pulmonarius cultivars, 36 cultivars and five wild isolates were subjected to de novo sequencing and assembly to generate high-quality mitogenome sequences. The P. pulmonarius mitogenomes had lengths ranging from 69,096 to 72,905 base pairs. The mitogenome sizes of P. pulmonarius and those of other mushroom species in the Pleurotus genus showed a significant positive correlation with the counts of LAGLIDAG and GIY-YIG homing endonucleases encoded by intronic open reading frames. A comparison of gene arrangements revealed an inversion of a fragment containing atp9-nad3-nad2 between P. pulmonarius and P. ostreatus. The mitogenomes of P. pulmonarius were clustered into three distinct clades, two of which were crowded with commercial cultivars. Clade I, all of which possess an inserted dpo gene, shared a maternal origin linked to an ancestral cultivar from Taiwan. Primers were designed to target the dpo gene, potentially safeguarding intellectual property rights. The wild isolates in Clade III exhibited more divergent mitogenomes, rendering them valuable for breeding.

Complete mitochondrial genome sequence of the white root rot pathogen Dematophora necatrix (Xylariaceae: Xylariales)

The mitochondrial genome of Dematophora necatrix is 121,350 base pairs in length with a G + C content of 30.19%. Phylogenetic analysis showed that D. necatrix grouped with other members of the Xylariaceae, with which its mitogenome also shares a broadly similar architecture and gene content. The D. necatrix mitogenome contains 14 protein-coding and 26 tRNA-encoding genes, as well as one copy each of the rnl, rns, rps3 and nat1 genes. However, as much as 80% of this genome is intronic or non-coding. This is likely due to expansions and rearrangements caused by the large number of group I introns and the homing endonucleases and reverse-transcriptases they encode. Our study thus provides a valuable foundation from which to explore the mitochondrion's role in the biology of D. necatrix, and also serves as a resource for investigating the pathogen's population biology and general ecology.

Insights into Fungal Mitochondrial Genomes and Inheritance Based on Current Findings from Yeast-like Fungi

The primary functions of mitochondria are to produce energy and participate in the apoptosis of cells, with them being highly conserved among eukaryotes. However, the composition of mitochondrial genomes, mitochondrial DNA (mtDNA) replication, and mitochondrial inheritance varies significantly among animals, plants, and fungi. Especially in fungi, there exists a rich diversity of mitochondrial genomes, as well as various replication and inheritance mechanisms. Therefore, a comprehensive understanding of fungal mitochondria is crucial for unraveling the evolutionary history of mitochondria in eukaryotes. In this review, we have organized existing reports to systematically describe and summarize the composition of yeast-like fungal mitochondrial genomes from three perspectives: mitochondrial genome structure, encoded genes, and mobile elements. We have also provided a systematic overview of the mechanisms in mtDNA replication and mitochondrial inheritance during bisexual mating. Additionally, we have discussed and proposed open questions that require further investigation for clarification.

Comprehensive investigations of 2-phenylethanol production by the filamentous fungus Annulohypoxylon stygium

2-Phenylethanol (2-PE) is an aromatic compound with a rose-like fragrance that is widely used in food and other industries. Yeasts have been implicated in the biosynthesis of 2-PE; however, few studies have reported the involvement of filamentous fungi. In this study, 2-PE was detected in Annulohypoxylon stygium mycelia grown in both potato dextrose broth (PDB) and sawdust medium. Among the 27 A. stygium strains investigated in this study, the strain "Jinjiling" (strain S20) showed the highest production of 2-PE. Under optimal culture conditions, the concentration of 2-PE was 2.33 g/L. Each of the key genes in Saccharomyces cerevisiae shikimate and Ehrlich pathways was found to have homologous genes in A. stygium. Upon the addition of L-phenylalanine to the medium, there was an upregulation of all key genes in the Ehrlich pathway of A. stygium, which was consistent with that of S. cerevisiae. A. stygium as an associated fungus provides nutrition for the growth of Tremella fuciformis and most spent composts of T. fuciformis contain pure A. stygium mycelium. Our study on the high-efficiency biosynthesis of 2-PE in A. stygium offers a sustainable solution by utilizing the spent compost of T. fuciformis and provides an alternative option for the production of natural 2-PE. KEY POINTS: • Annulohypoxylon stygium can produce high concentration of 2-phenylethanol. • The pathways of 2-PE biosynthesis in Annulohypoxylon stygium were analyzed. • Spent compost of Tremella fuciformis is a potential source for 2-phenylethanol.

The mitogenomes of Leptographium aureum, Leptographium sp., and Grosmannia fruticeta: expansion by introns

Introduction

Many members of the Ophiostomatales are of economic importance as they are bark-beetle associates and causative agents for blue stain on timber and in some instances contribute towards tree mortality. The taxonomy of these fungi has been challenging due to the convergent evolution of many traits associated with insect dispersal and a limited number of morphological characters that happen to be highly pleomorphic. This study examines the mitochondrial genomes for three members of Leptographium sensu lato [Leptographium aureum (also known as Grosmannia aurea), Grosmannia fruticeta (also known as Leptographium fruticetum), and Leptographium sp. WIN(M)1376)].

Methods

Illumina sequencing combined with gene and intron annotations and phylogenetic analysis were performed.

Results

Sequence analysis showed that gene content and gene synteny are conserved but mitochondrial genome sizes were variable: G. fruticeta at 63,821 bp, Leptographium sp. WIN(M)1376 at 81,823 bp and L. aureum at 104,547 bp. The variation in size is due to the number of introns and intron-associated open reading frames. Phylogenetic analysis of currently available mitochondrial genomes for members of the Ophiostomatales supports currently accepted generic arrangements within this order and specifically supports the separation of members with Leptographium-like conidiophores into two genera, with L. aureum grouping with Leptographium and G. fruticeta aligning with Grosmannia.

Discussion

Mitochondrial genomes are promising sequences for resolving evolutionary relationships within the Ophiostomatales.

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.

The first two mitochondrial genomes from Apiotrichum reveal mitochondrial evolution and different taxonomic assignment of Trichosporonales

Apiotrichum is a diverse anamorphic basidiomycetous yeast genus, and its mitogenome characterization has not been revealed. In this study, we assembled two Apiotrichum mitogenomes and compared them with mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. The mitogenomes of Apiotrichum gracile and A. gamsii comprised circular DNA molecules, with sizes of 34,648 bp and 38,096 bp, respectively. Intronic regions were found contributed the most to the size expansion of A. gamsii mitogenome. Comparative mitogenomic analysis revealed that 6.85-38.89% of nucleotides varied between tRNAs shared by the two Apiotrichum mitogenomes. The GC content of all core PCGs in A. gamsii was lower than that of A. gracile, with an average low value of 4.97%. The rps3 gene differentiated the most among Agaricomycotina, Pucciniomycotina and Ustilaginomycotina species, while nad4L gene was the most conserved in evolution. The Ka/Ks values for cob and rps3 genes were > 1, indicating the two genes may be subjected to positive selection in Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. Frequent intron loss/gain events and potential intron transfer events have been detected in evolution of Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. We further detected large-scale gene rearrangements between the 19 mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina, and fifteen of the 17 mitochondrial genes shared by Apiotrichum varied in gene arrangements. Phylogenetic analyses based on maximum likelihood and Bayesian inference methods using a combined mitochondrial gene dataset revealed different taxonomic assignment of two Apiotrichum species, wherein A. gamsii had a more closely relationship with Trichosporon asahii. This study served as the first report on mitogenomes from the genus Apiotrichum, which promotes the understanding of evolution, genomics, and phylogeny of Apiotrichum.

Comparative analyses of Flammulina filiformis mitochondrial genomes reveal high length polymorphism in intergenic regions and multiple intron gain/loss in cox1

The golden-needle mushroom Flammulina filiformis is one of the bulk mushroom products in the world. This study obtained complete mitogenomes of 44 wild isolates collected from nine provinces and two artificially bred cultivars of F. filiformis, together with three Flammulina rossica isolates and one Flammulina fennae isolate for comparison. The mitogenome of F. filiformis ranged from 83,540 bp to 90,938 bp, consisting of 14 conserved protein-coding genes (PCGs), two rRNA genes, and 25 tRNA genes. To the best of our knowledge, it contained the highest proportion of intergenic regions compared to the other known Basidiomycota mitogenomes. Introns and intergenic regions were two major contributing factors to the total size of the F. filiformis mitogenome. The conserved PCG cox3 is located in an intron of another conserved PCG, nad5. This is a unique phenomenon in all known fungal mitogenomes. Gain/loss of introns was observed in cox1, nad5, and rnl. Length polymorphism was widely observed in intergenic regions. Accordingly, primers were designed as useful markers for rapid identification of F. filiformis isolates with differentiated mitogenomes. Our findings provide a basis for further studies related to variety identification and population genetics of this economically important mushroom.

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.

Comparative Mitogenomic Analysis Reveals Intraspecific, Interspecific Variations and Genetic Diversity of Medical Fungus Ganoderma

Ganoderma species are widely distributed in the world with high diversity. Some species are considered to be pathogenic fungi while others are used as traditional medicine in Asia. In this study, we sequenced and assembled four Ganoderma complete mitogenomes, including G. subamboinense s118, G. lucidum s37, G. lingzhi s62, and G. lingzhi s74. The sizes of the four mitogenomes ranged from 50,603 to 73,416 bp. All Ganoderma specimens had a full set of core protein-coding genes (PCGs), and the rps3 gene of Ganoderma species was detected to be under positive or relaxed selection. We found that the non-conserved PCGs, which encode RNA polymerases, DNA polymerases, homing endonucleases, and unknown functional proteins, are dynamic within and between Ganoderma species. Introns were thought to be the main contributing factor in Ganoderma mitogenome size variation (p < 0.01). Frequent intron loss/gain events were detected within and between Ganoderma species. The mitogenome of G. lucidum s26 gained intron P637 in the cox3 gene compared with the other two G. lucidum mitogenomes. In addition, some rare introns in Ganoderma were detected in distinct Basidiomycetes, indicating potential gene transfer events. Comparative mitogenomic analysis revealed that gene arrangements also varied within and between Ganoderma mitogenomes. Using maximum likelihood and Bayesian inference methods with a combined mitochondrial gene dataset, phylogenetic analyses generated identical, well-supported tree topologies for 71 Agaricomycetes species. This study reveals intraspecific and interspecific variations of the Ganoderma mitogenomes, which promotes the understanding of the origin, evolution, and genetic diversity of Ganoderma species.

The mitogenome of Urnula craterium

Urnula craterium (Schwein.) Fr. (1851) has been reported from North America, Europe, and Asia, and can be a pathogen on various hardwood species. In this study we investigated the mitochondrial genome of U. craterium. The biology and taxonomy of this fungus is poorly studied and there are no mitogenomes currently available for any member of the Sarcosomataceae (Order Pezizales). The complete mitogenome of U. craterium comprises 43 967 bps and encodes 14 protein-coding genes, a complete set of tRNAs and rRNA genes. A novel feature of the mitogenome is the presence of a single subunit DNA polymerase coding region that is typically associated with linear invertron-type plasmids. The mitogenome may offer insights into the evolution of mitogenomes among members of the Pezizales with regards to gene content and order, mobile elements, and genome sizes.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

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

First two mitochondrial genomes for the order Filobasidiales reveal novel gene rearrangements and intron dynamics of Tremellomycetes

In the present study, two mitogenomes from the Filobasidium genus were assembled and compared with other Tremellomycetes mitogenomes. The mitogenomes of F. wieringae and F. globisporum both comprised circular DNA molecules, with sizes of 27,861 bp and 71,783 bp, respectively. Comparative mitogenomic analysis revealed that the genetic contents, tRNAs, and codon usages of the two Filobasidium species differed greatly. The sizes of the two Filobasidium mitogenomes varied greatly with the introns being the main factor contributing to mitogenome expansion in F. globisporum. Positive selection was observed in several protein-coding genes (PCGs) in the Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina species, including cob, cox2, nad2, and rps3 genes. Frequent intron loss/gain events were detected to have occurred during the evolution of the Tremellomycetes mitogenomes, and the mitogenomes of 17 species from Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina have undergone large-scale gene rearrangements. Phylogenetic analyses based on Bayesian inference and the maximum likelihood methods using a combined mitochondrial gene set generated identical and well-supported phylogenetic trees, wherein Filobasidium species had close relationships with Trichosporonales species. This study, which is the first report on mitogenomes from the order Filobasidiales, provides a basis for understanding the genomics, evolution, and taxonomy of this important fungal group.

Mitogenome of a stink worm (Annelida: Travisiidae) includes degenerate group II intron that is also found in five congeneric species

Mitogenomes are useful for inferring phylogenetic relationships between organisms. Although the mitogenomes of Annelida, one of the most morphologically and ecologically diverse metazoan groups have been well sequenced, those of several families remain unexamined. This study determined the first mitogenome from the family Travisiidae (Travisia sanrikuensis), analyzed its mitogenomic features, and reconstructed a phylogeny of Sedentaria. The monophyly of the Terebellida + Arenicolida + Travisiidae clade is supported by molecular phylogenetic analysis. The placement of Travisiidae is unclear because of the lack of mitogenomes from closely related lineages. An unexpected intron appeared within the cox1 gene of T. sanrikuensis and in the same positions of five undescribed Travisia spp. Although the introns are shorter (790–1386 bp) than other group II introns, they can be considered degenerate group II introns due to type II intron maturase open reading frames, found in two of the examined species, and motifs characteristic of group II introns. This is likely the first known case in metazoans where mitochondrial group II introns obtained by a common ancestor are conserved in several descendants. Insufficient evolutionary time for intron loss in Travisiidae, or undetermined mechanisms may have helped maintain the degenerate introns.

The First Mitochondrial Genome of Ciborinia camelliae and Its Position in the Sclerotiniaceae Family

Ciborinia camelliae is the causal agent of camellia flower blight (CFB). It is a hemibiotrophic pathogen, inoperculate Discomycete of the family Sclerotiniaceae. It shows host and organ specificity infecting only flowers of species belonging to the genus Camellia, causing serious damage to the ornamental component of the plant. In this work, the first mitochondrial genome of Ciborinia camellia is reported. The mitogenome was obtained by combining Illumina short read and Nanopore long read technology. To resolve repetitive elements, specific primers were designed and used for Sanger sequencing. The manually curated mitochondrial DNA (mtDNA) of the Italian strain DSM 112729 is a circular sequence of 114,660 bp, with 29.6% of GC content. It contains two ribosomal RNA genes, 33 transfer RNAs, one RNase P gene, and 62 protein-coding genes. The latter include one gene coding for a ribosomal protein (rps3) and the 14 typical proteins involved in the oxidative metabolism. Moreover, a partial mtDNA assembled from a contig list was obtained from the deposited genome assembly of a New Zealand strain of C. camelliae. The present study contributes to understanding the mitogenome arrangement and the evolution of this phytopathogenic fungus in comparison to other Sclerotiniaceae species and confirms the usefulness of mitochondrial analysis to define phylogenetic positioning of this newly sequenced species.

Characterization of the mitochondrial genomes of three powdery mildew pathogens reveals remarkable variation in size and nucleotide composition

Powdery mildews comprise a large group of economically important phytopathogenic fungi. However, limited information exists on their mitochondrial genomes. Here, we assembled and compared the mitochondrial genomes of the powdery mildew pathogens Blumeria graminis f. sp. tritici, Erysiphe pisi, and Golovinomyces cichoracearum. Included in the comparative analysis was also the mitochondrial genome of Erysiphe necator that was previously analysed. The mitochondrial genomes of the four Erysiphales exhibit a similar gene content and organization but a large variation in size, with sizes ranging from 109800 bp in B. graminis f. sp. tritici to 332165 bp in G. cichoracearum, which is the largest mitochondrial genome of a fungal pathogen reported to date. Further comparative analysis revealed an unusual bimodal GC distribution in the mitochondrial genomes of B. graminis f. sp. tritici and G. cichoracearum that was not previously observed in fungi. The cytochrome b (cob) genes of E. necator, E. pisi, and G. cichoracearum were also exceptionally rich in introns, which in turn harboured rare open reading frames encoding reverse transcriptases that were likely acquired horizontally. Golovinomyces cichoracearum had also the longest cob gene (45 kb) among 703 fungal cob genes analysed. Collectively, these results provide novel insights into the organization of mitochondrial genomes of powdery mildew pathogens and represent valuable resources for population genetics and evolutionary studies.

Comparative analysis of Malassezia furfur mitogenomes and the development of a mitochondria-based typing approach

Malassezia furfur is a yeast species belonging to Malasseziomycetes, Ustilaginomycotina and Basidiomycota that is found on healthy warm-blooded animal skin, but also involved in various skin disorders like seborrheic dermatitis/dandruff and pityriasis versicolor. Moreover, Malassezia are associated with bloodstream infections, Crohn's disease and pancreatic carcinoma. Recent advances in Malassezia genomics and genetics have focused on the nuclear genome. In this work, we present the M. furfur mitochondrial (mt) genetic heterogenicity with full analysis of 14 novel and six available M. furfur mt genomes. The mitogenome analysis reveals a mt gene content typical for fungi, including identification of variable mt regions suitable for intra-species discrimination. Three of them, namely the trnK–atp6 and cox3–nad3 intergenic regions and intron 2 of the cob gene, were selected for primer design to identify strain differences. Malassezia furfur strains belonging to known genetic variable clusters, based on AFLP and nuclear loci, were assessed for their mt variation using PCR amplification and sequencing. The results suggest that these mt regions are excellent molecular markers for the typing of M. furfur strains and may provide added value to nuclear regions when assessing evolutionary relationships at the intraspecies level.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

An Update on Trichoderma Mitogenomes: Complete De Novo Mitochondrial Genome of the Fungal Biocontrol Agent Trichoderma harzianum (Hypocreales, Sordariomycetes), an Ex-Neotype Strain CBS 226.95, and Tracing the Evolutionary Divergences of Mitogenomes in Trichoderma

Members of the genus Trichoderma (Hypocreales), widely used as biofungicides, biofertilizers, and as model fungi for the industrial production of CAZymes, have actively been studied for the applications of their biological functions. Recently, the study of the nuclear genomes of Trichoderma has expanded in the directions of adaptation and evolution to gain a better understanding of their ecological traits. However, Trichoderma's mitochondria have received much less attention despite mitochondria being the most necessary element for sustaining cell life. In this study, a mitogenome of the fungus Trichoderma harzianum CBS 226.95 was assembled de novo. A 27,632 bp circular DNA molecule was revealed with specific features, such as the intronless of all core PCGs, one homing endonuclease, and a putative overlapping tRNA, on a closer phylogenetic relationship with T. reesei among hypocrealean fungi. Interestingly, the mitogenome of T. harzianum CBS 226.95 was predicted to have evolved earlier than those of other Trichoderma species and also assumed with a selection pressure in the cox3. Considering the bioavailability, both for the ex-neotype strain of the T. harzianum species complex and the most globally representative commercial fungal biocontrol agent, our results on the T. harzianum CBS 226.95 mitogenome provide crucial information which will be helpful criteria in future studies on Trichoderma.

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.

Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes

The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.

Intraspecific Mitochondrial DNA Comparison of Mycopathogen Mycogone perniciosa Provides Insight Into Mitochondrial Transfer RNA Introns

Mycogone perniciosa is the main causative agent of wet bubble disease, which causes severe damage to the production of the cultivated mushroom Agaricus bisporus around the world. Whole-genome sequencing of 12 isolates of M. perniciosa was performed using the Illumina sequencing platform, and the obtained paired-end reads were used to assemble complete mitochondrial genomes. Intraspecific comparisons of conserved protein-coding genes, transfer RNA (tRNA) and ribosomal RNA (rRNA) genes, introns, and intergenic regions were conducted. Five different mitochondrial DNA (mtDNA) haplotypes were detected among the tested isolates, ranging from 89,080 to 93,199 bp in length. All of the mtDNAs contained the same set of 14 protein-coding genes and 2 rRNA and 27 tRNA genes, which shared high sequence similarity. In contrast, the number, insertion sites, and sequences of introns varied greatly among the mtDNAs. Eighteen of 43 intergenic regions differed among the isolates, reflecting 65 single nucleotide polymorphisms, 76 indels, and the gain/loss of nine long fragments. Intraspecific comparison revealed that two introns were located within tRNA genes, which is the first detailed description of mitochondrial tRNA introns. Intronic sequence comparison within the same insertion sites revealed the formation process of two introns, which also illustrated a fast evolutionary rate of introns among M. perniciosa isolates. Based on the intron distribution pattern, a pair of universal primers and four pairs of isolate-specific primers were designed and were used to identify the five mtDNA types. In summary, the rapid gain or loss of mitochondrial introns could be an ideal marker for population genetics analysis.

Comparative Mitogenomic Analysis Reveals Dynamics of Intron Within and Between Tricholoma Species and Phylogeny of Basidiomycota

The genus of Tricholoma is a group of important ectomycorrhizal fungi. The overlapping of morphological characteristics often leads to the confusion of Tricholoma species classification. In this study, the mitogenomes of five Tricholoma species were sequenced based on the next-generation sequencing technology, including T. matsutake SCYJ1, T. bakamatsutake, T. terreum, T. flavovirens, and T. saponaceum. These five mitogenomes were all composed of circular DNA molecules, with sizes ranging from 49,480 to 103,090 bp. Intergenic sequences were considered to be the main factor contributing to size variations of Tricholoma mitogenomes. Comparative mitogenomic analysis showed that the introns of the Agaricales mitogenome experienced frequent loss/gain events. In addition, potential gene transfer was detected between the mitochondrial and nuclear genomes of the five species of Tricholoma. Evolutionary analysis showed that the rps3 gene of the Tricholoma species was under positive selection or relaxed selection in the evolutionary process. In addition, large-scale gene rearrangements were detected between some Tricholoma species. Phylogenetic analysis using the Bayesian inference and maximum likelihood methods based on a combined mitochondrial gene set yielded identical and well-supported tree topologies. This study promoted the understanding of the genetics, evolution, and phylogeny of the Tricholoma genus and related species.

The mitochondrial genome of Ophiostoma himal-ulmi and comparison with other fungi causing Dutch elm disease

The mitochondrial genome of Ophiostoma himal-ulmi, a species endemic to the Western Himalayas and one of the fungi that cause Dutch elm disease, has been sequenced and characterized. The mitochondrial genome was compared with other available genomes for members of the Ophiostomatales, including other agents of Dutch elm disease (Ophiostoma ulmi, Ophiostoma novo-ulmi subspecies novo-ulmi, and Ophiostoma novo-ulmi subspecies americana), and it was observed that gene synteny is highly conserved, and variability among members of the fungi that cause Dutch-elm disease is primarily due to the number of intron insertions. Among the fungi that cause Dutch elm disease that we examined, O. himal-ulmi has the largest mitochondrial genomes (ranging from 94 934 to 111 712 bp), owing to the expansion of the number of introns.

The Mitogenomes of Ophiostoma minus and Ophiostoma piliferum and Comparisons With Other Members of the Ophiostomatales

Fungi assigned to the Ophiostomatales are of economic concern as many are blue-stain fungi and some are plant pathogens. The mitogenomes of two blue-stain fungi, Ophiostoma minus and Ophiostoma piliferum, were sequenced and compared with currently available mitogenomes for other members of the Ophiostomatales. Species representing various genera within the Ophiostomatales have been examined for gene content, gene order, phylogenetic relationships, and the distribution of mobile elements. Gene synteny is conserved among the Ophiostomatales but some members were missing the atp9 gene. A genome wide intron landscape has been prepared to demonstrate the distribution of the mobile genetic elements (group I and II introns and homing endonucleases) and to provide insight into the evolutionary dynamics of introns among members of this group of fungi. Examples of complex introns or nested introns composed of two or three intron modules have been observed in some species. The size variation among the mitogenomes (from 23.7 kb to about 150 kb) is mostly due to the presence and absence of introns. Members of the genus Sporothrix sensu stricto appear to have the smallest mitogenomes due to loss of introns. The taxonomy of the Ophiostomatales has recently undergone considerable revisions; however, some lineages remain unresolved. The data showed that genera such as Raffaelea appear to be polyphyletic and the separation of Sporothrix sensu stricto from Ophiostoma is justified.

Comparative mitochondrial genome analysis reveals intron dynamics and gene rearrangements in two Trametes species

Trametes species are efficient wood decomposers that are widespread throughout the world. Mitogenomes have been widely used to understand the phylogeny and evolution of fungi. Up to now, two mitogenomes from the Trametes genus have been revealed. In the present study, the complete mitogenomes of two novel Trametes species, Trametes versicolor and T. coccinea, were assembled and compared with other Polyporales mitogenomes. Both species contained circular DNA molecules, with sizes of 67,318 bp and 99,976 bp, respectively. Comparative mitogenomic analysis indicated that the gene number, length and base composition varied between the four Trametes mitogenomes we tested. In addition, all of the core protein coding genes in Trametes species were identified and subjected to purifying selection. The mitogenome of T. coccinea contained the largest number of introns among the four Trametes species tested, and introns were considered the main factors contributing to size variations of Polyporales. Several novel introns were detected in the Trametes species we assembled, and introns identified in Polyporales were found to undergo frequent loss/gain events. Large-scale gene rearrangements were detected between closely related Trametes species, including gene inversions, insertions, and migrations. A well-supported phylogenetic tree for 77 Basidiomycetes was obtained based on the combined mitochondrial gene set using 2 phylogenetic inference methods. The results showed that mitochondrial genes are effective molecular markers for understanding the phylogeny of Basidiomycetes. This study is the first to report the mitogenome rearrangement and intron dynamics of Trametes species, which shed light on the evolution of Trametes and other related species.

The mitogenomes of two saprophytic Boletales species (Coniophora) reveals intron dynamics and accumulation of plasmid-derived and non-conserved genes

The order Boletales is a group of fungi with complex life styles, which include saprophytic and ectomycorrhizal mushroom-forming fungi. In the present study, the complete mitogenomes of two saprophytic Boletales species, Coniophora olivacea, and C. puteana, were assembled and compared with mitogenomes of ectomycorrhizal Boletales. Both mitogenomes comprised circular DNA molecules with sizes of 78,350 bp and 79,655 bp, respectively. Comparative mitogenomic analysis indicated that the two saprophytic Boletales species contained more plasmid-derived (7 on average) and unknown functional genes (12 on average) than the four ectomycorrhizal Boletales species previously reported. In addition, the core protein coding genes, nad2 and rps3, were found to be subjected to positive selection pressure between some Boletales species. Frequent intron gain/loss events were detected in Boletales and Basidiomycetes, and several novel intron classes were found in two Coniophora species. A total of 33 introns were detected in C. olivacea, and most were found to have undergone contraction in the C. olivacea mitogenome. Mitochondrial genes of Coniophora species were found to have undergone large-scale gene rearrangements, and the accumulation of intra-genomic repeats in the mitogenome was considered as one of the main contributing factors. Based on combined mitochondrial gene sets, we obtained a well-supported phylogenetic tree for 76 Basidiomycetes, demonstrating the utility of mitochondrial gene analysis for inferring Basidiomycetes phylogeny. The study served as the first report on the mitogenomes of the family Coniophorineae, which will help to understand the origin and evolution patterns of Boletales species with complex lifestyles.

Mitochondrial evolution in the entomopathogenic fungal genus Beauveria

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

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.

The 256 kb mitochondrial genome of Clavaria fumosa is the largest among phylum Basidiomycota and is rich in introns and intronic ORFs

In the present study, the complete mitogenome of Clavaria fumosa, was sequenced, assembled, and compared. The complete mitogenome of C. fumosa is 256,807 bp in length and is the largest mitogenomes among all Basidiomycota mitogenomes reported. Comparative mitogenomic analysis indicated that the C. fumosa mitogenome contained the most introns and intronic ORFs among all fungal mitogenomes. Large intergenic regions, intronic regions, accumulation of repeat sequences and plasmid-derived genes together promoted the size expansion of the C. fumosa mitogenome. In addition, the rps3 gene was found subjected to positive selection between some Agaricales species. We found frequent intron gain/loss events in Agaricales mitogenomes, and four novel intron classes were detected in the C. fumosa mitogenome. Large-scale gene rearrangements were found occurred in Agaricales species and the C. fumosa mitogenome had a unique gene arrangement which differed from other Agaricales species. Phylogenetic analysis for 76 Basidiomycetes based on combined mitochondrial gene sets indicated that mitochondrial genes could be used as effective molecular markers for reconstructing evolution of Basidiomycota. The study served as the first report on the mitogenomes of the family Clavariaceae, which will promote the understanding of the genetics, evolution and taxonomy of C. fumosa and related species.

The complete mitochondrial genome of medicinal fungus Taiwanofungus camphoratus reveals gene rearrangements and intron dynamics of Polyporales

Taiwanofungus camphoratus is a highly valued medicinal mushroom that is endemic to Taiwan, China. In the present study, the mitogenome of T. camphoratus was assembled and compared with other published Polyporales mitogenomes. The T. camphoratus mitogenome was composed of circular DNA molecules, with a total size of 114,922 bp. Genome collinearity analysis revealed large-scale gene rearrangements between the mitogenomes of Polyporales, and T. camphoratus contained a unique gene order. The number and classes of introns were highly variable in 12 Polyporales species we examined, which proved that numerous intron loss or gain events occurred in the evolution of Polyporales. The Ka/Ks values for most core protein coding genes in Polyporales species were less than 1, indicating that these genes were subject to purifying selection. However, the rps3 gene was found under positive or relaxed selection between some Polyporales species. Phylogenetic analysis based on the combined mitochondrial gene set obtained a well-supported topology, and T. camphoratus was identified as a sister species to Laetiporus sulphureus. This study served as the first report on the mitogenome in the Taiwanofungus genus, which will provide a basis for understanding the phylogeny and evolution of this important fungus.

Characterization of Two Mitochondrial Genomes and Gene Expression Analysis Reveal Clues for Variations, Evolution, and Large-Sclerotium Formation in Medical Fungus Wolfiporia cocos

Wolfiporia cocos, a precious mushroom with a long history as an edible food and Asian traditional medicine, remains unclear in the genetic mechanism underlying the formation of large sclerotia. Here, two complete circular mitogenomes (BL16, 135,686 bp and MD-104 SS10, 124,842 bp, respectively) were presented in detail first. The salient features in the mitogenomes of W. cocos include an intron in the tRNA (trnQ-UUG²), and an obvious gene rearrangement identified between the two mitogenomes from the widely geographically separated W. cocos strains. Genome comparison and phylogenetic analyses reveal some variations and evolutional characteristics in W. cocos. Whether the mitochondrion is functional in W. cocos sclerotium development was investigated by analyzing the mitogenome synteny of 10 sclerotium-forming fungi and mitochondrial gene expression patterns in different W. cocos sclerotium-developmental stages. Three common homologous genes identified across ten sclerotium-forming fungi were also found to exhibit significant differential expression levels during W. cocos sclerotium development. Most of the mitogenomic genes are not expressed in the mycelial stage but highly expressed in the sclerotium initial or developmental stage. These results indicate that some of mitochondrial genes may play a role in the development of sclerotium in W. cocos, which needs to be further elucidated in future studies. This study will stimulate new ideas on cytoplasmic inheritance of W. cocos and facilitate the research on the role of mitochondria in large sclerotium formation.

The First Mitochondrial Genome for Geastrales (Sphaerobolus stellatus) Reveals Intron Dynamics and Large-Scale Gene Rearrangements of Basidiomycota

In this study, the mitogenome of artillery fungus, Sphaerobolus stellatus, was assembled and compared with other Basidiomycota mitogenomes. The Sphaerobolus stellatus mitogenome was composed of circular DNA molecules, with a total size of 152,722 bp. Accumulation of intergenic and intronic sequences contributed to the Sphaerobolus stellatus mitogenome becoming the fourth largest mitogenome among Basidiomycota. We detected large-scale gene rearrangements in Basidiomycota mitogenomes, and the Sphaerobolus stellatus mitogenome contains a unique gene order. The quantity and position classes of intron varied between 75 Basidiomycota species we tested, indicating frequent intron loss/gain events occurred in the evolution of Basidiomycota. A novel intron position classes (P1281) was detected in the Sphaerobolus stellatus mitogenome, without any homologous introns from other Basidiomycota species. A pair of fragments with a total length of 9.12 kb in both the nuclear and mitochondrial genomes of Sphaerobolus stellatus was detected, indicating possible gene transferring events. Phylogenetic analysis based on the combined mitochondrial gene set obtained well-supported tree topologies (Bayesian posterior probabilities ≥ 0.99; bootstrap values ≥98). This study served as the first report on the mitogenome from the order Geastrales, which will promote the understanding of the phylogeny, population genetics, and evolution of the artillery fungus, Sphaerobolus stellatus.

Comparative mitogenome analysis of two ectomycorrhizal fungi (Paxillus) reveals gene rearrangement, intron dynamics, and phylogeny of basidiomycetes

In this study, the mitogenomes of two Paxillus species were assembled, annotated and compared. The two mitogenomes of Paxillus involutus and P. rubicundulus comprised circular DNA molecules, with the size of 39,109 bp and 41,061 bp, respectively. Evolutionary analysis revealed that the nad4L gene had undergone strong positive selection in the two Paxillus species. In addition, 10.64 and 36.50% of the repetitive sequences were detected in the mitogenomes of P. involutus and P. rubicundulus, respectively, which might transfer between mitochondrial and nuclear genomes. Large-scale gene rearrangements and frequent intron gain/loss events were detected in 61 basidiomycete species, which revealed large variations in mitochondrial organization and size in Basidiomycota. In addition, the insertion sites of the basidiomycete introns were found to have a base preference. Phylogenetic analysis of the combined mitochondrial gene set gave identical and well-supported tree topologies, indicating that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships of Basidiomycota. This study is the first report on the mitogenomes of Paxillus, which will promote a better understanding of their contrasted ecological strategies, molecular evolution and phylogeny of these important ectomycorrhizal fungi and related basidiomycete species.

Comparative Mitogenome Analysis Reveals Mitochondrial Genome Differentiation in Ectomycorrhizal and Asymbiotic Amanita Species

In this present study, we assembled and analyzed the mitogenomes of two asymbiotic and six ectomycorrhizal Amanita species based on next-generation sequencing data. The size of the eight Amanita mitogenomes ranged from 37,341 to 137,428 bp, and we considered introns to be one of the main factors contributing to the size variation of Amanita. The introns of the cox1 gene experienced frequent gain/loss events in Amanita; and the intron position class cox1P386 was lost in the six ectomycorrhizal Amanita species. In addition, ectomycorrhizal Amanita species had more repetitive sequences and fewer intergenic sequences than asymbiotic Amanita species in their mitogenomes. Large-scale gene rearrangements were detected in the Amanita species we tested, including gene displacements and inversions. On the basis of the combined mitochondrial gene set, we reconstructed the phylogenetic relationships of 66 Basidiomycetes. The six ectomycorrhizal Amanita species were of single origin, and the two saprophytic Amanita species formed two distinct clades. This study is the first to elucidate the functions of the mitogenome in the evolution and ecological adaptation of Amanita species.

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.

A functional twintron, 'zombie' twintrons and a hypermobile group II intron invading itself in plant mitochondria

The occurrence of group II introns in plant mitochondrial genomes is strikingly different between the six major land plant clades, contrasting their highly conserved counterparts in chloroplast DNA. Their present distribution likely reflects numerous ancient intron gains and losses during early plant evolution before the emergence of seed plants. As a novelty for plant organelles, we here report on five cases of twintrons, introns-within-introns, in the mitogenomes of lycophytes and hornworts. An internal group II intron interrupts an intron-borne maturase of an atp9 intron in Lycopodiaceae, whose splicing precedes splicing of the external intron. An invasive, hypermobile group II intron in cox1, has conquered nine further locations including a previously overlooked sdh3 intron and, most surprisingly, also itself. In those cases, splicing of the external introns does not depend on splicing of the internal introns. Similar cases are identified in the mtDNAs of hornworts. Although disrupting a group I intron-encoded protein in one case, we could not detect splicing of the internal group II intron in this ‘mixed’ group I/II twintron. We suggest the name ‘zombie’ twintrons (half-dead, half-alive) for such cases where splicing of external introns does not depend any more on prior splicing of fossilized internal introns.

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

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

The complete mitochondrial genome sequence of Scolopendra mutilans L. Koch, 1878 (Scolopendromorpha, Scolopendridae), with a comparative analysis of other centipede genomes

Scolopendramutilans L. Koch, 1878 is an important Chinese animal with thousands of years of medicinal history. However, the genomic information of this species is limited, which hinders its further application. Here, the complete mitochondrial genome (mitogenome) of S.mutilans was sequenced and assembled by next-generation sequencing. The genome is 15,011 bp in length, consisting of 13 protein-coding genes (PCGs), 14 tRNA genes, and two rRNA genes. Most PCGs start with the ATN initiation codon, and all PCGs have the conventional stop codons TAA and TAG. The S.mutilans mitogenome revealed nine simple sequence repeats (SSRs), and an obviously lower GC content compared with other seven centipede mitogenomes previously sequenced. After analysis of homologous regions between the eight centipede mitogenomes, the S.mutilans mitogenome further showed clear genomic rearrangements. The phylogenetic analysis of eight centipedes using 13 conserved PCG genes was finally performed. The phylogenetic reconstructions showed Scutigeromorpha as a separate group, and Scolopendromorpha in a sister-group relationship with Lithobiomorpha and Geophilomorpha. Collectively, the S.mutilans mitogenome provided new genomic resources, which will improve its medicinal research and applications in the future.

Evidence of Extensive Intraspecific Noncoding Reshuffling in a 169-kb Mitochondrial Genome of a Basidiomycetous Fungus

Comparative genomics of fungal mitochondrial genomes (mitogenomes) have revealed a remarkable pattern of rearrangement between and within major phyla owing to horizontal gene transfer and recombination. The role of recombination was exemplified at a finer evolutionary time scale in basidiomycetes group of fungi as they display a diversity of mitochondrial DNA inheritance patterns. Here, we assembled mitogenomes of six species from the Hymenochaetales order of basidiomycetes and examined 59 mitogenomes from 2 genetic lineages of Phellinus noxius. Gene order is largely collinear, while intergene regions are major determinants of mitogenome size variation. Substantial sequence divergence was found in shared introns consistent with high horizontal gene transfer frequency observed in yeasts, but we also identified a rare case where an intron was retained in five species since speciation. In contrast to the hyperdiversity observed in nuclear genomes of Phellinus noxius, mitogenomes’ intraspecific polymorphisms at protein-coding sequences are extremely low. Phylogeny network based on introns revealed turnover as well as exchange of introns between two lineages. Strikingly, some strains harbor a mosaic origin of introns from both lineages. Analysis of intergenic sequence indicated substantial differences between and within lineages, and an expansion may be ongoing as a result of exchange between distal intergenes. These findings suggest that the evolution in mitochondrial DNAs is usually lineage specific but chimeric mitotypes are frequently observed, thus capturing the possible evolutionary processes shaping mitogenomes in a basidiomycete. The large mitogenome sizes reported in various basidiomycetes appear to be a result of interspecific reshuffling of intergenes.

Natural diversity in the predatory behavior facilitates the establishment of a robust model strain for nematode-trapping fungi

Nematode-trapping fungi (NTF) are a group of specialized microbial predators that consume nematodes when food sources are limited. Predation is initiated when conserved nematode ascaroside pheromones are sensed, followed by the development of complex trapping devices. To gain insights into the coevolution of this interkingdom predator-prey relationship, we investigated natural populations of nematodes and NTF that we found to be ubiquitous in soils. Arthrobotrys species were sympatric with various nematode species and behaved as generalist predators. The ability to sense prey among wild isolates of Arthrobotrys oligospora varied greatly, as determined by the number of traps after exposure to Caenorhabditis elegans While some strains were highly sensitive to C. elegans and the nematode pheromone ascarosides, others responded only weakly. Furthermore, strains that were highly sensitive to the nematode prey also developed traps faster. The polymorphic nature of trap formation correlated with competency in prey killing, as well as with the phylogeny of A. oligospora natural strains, calculated after assembly and annotation of the genomes of 20 isolates. A chromosome-level genome assembly and annotation were established for one of the most sensitive wild isolates, and deletion of the only G-protein β-subunit-encoding gene of A. oligospora nearly abolished trap formation. In summary, our study establishes a highly responsive A. oligospora wild isolate as a model strain for the study of fungus-nematode interactions and demonstrates that trap formation is a fitness character in generalist predators of the nematode-trapping fungus family.

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

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

The fungal mitochondrial Nad5 pan-genic intron landscape

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

The complete mitochondrial genomes of two model ectomycorrhizal fungi (Laccaria): features, intron dynamics and phylogenetic implications

Laccaria amethystine and L. bicolor have served as model species for studying the life history and genetics of ectomycorrhizal fungi. However, the characterizations and variations of their mitogenomes are still unknown. In the present study, the mitogenomes of the two Laccaria species were assembled, annotated, and compared. The two mitogenomes of L. amethystine and L. bicolor comprised circular DNA molecules, with the sizes of 65,156 bp and 95,304 bp, respectively. Genome collinearity analysis revealed large-scale gene rearrangements between the two Laccaria species. Comparative mitogenome analysis indicated the introns of cox1 genes in Agaricales experienced frequent lost/gain eveants, which promoted the organization and size variations in Agaricales mitogenomes. Evolutionary analysis indicated the core protein-coding genes in the two mitogenomes were subject to strong pressure of purifying selection. Phylogenetic analysis using the Bayesian inference (BI) and Maximum likelihood (ML) methods based on a combined mitochondrial gene set resulted in identical and well-supported tree topologies, wherein the two Laccaria species were most closely related to Coprinopsis cinerea. This study severed as the first study on the mitogenomes of Laccaria species, which promoted a comprehensive understanding of the genetics and evolution of the model ectomycorrhizal fungi.

Comparative Mitochondrial Genome Analysis of Two Ectomycorrhizal Fungi (Rhizopogon) Reveals Dynamic Changes of Intron and Phylogenetic Relationships of the Subphylum Agaricomycotina

In the present study, we assembled and compared two mitogenomes from the Rhizopogon genus. The two mitogenomes of R. salebrosus and R. vinicolor comprised circular DNA molecules, with the sizes of 66,704 bp and 77,109 bp, respectively. Comparative mitogenome analysis indicated that the length and base composition of protein coding genes (PCGs), rRNA genes and tRNA genes varied between the two species. Large fragments aligned between the mitochondrial and nuclear genomes of both R. salebrosus (43.41 kb) and R. vinicolor (12.83 kb) indicated that genetic transfer between mitochondrial and nuclear genomes has occurred over evolutionary time of Rhizopogon species. Intronic regions were found to be the main factors contributing to mitogenome expansion in R. vinicolor. Variations in the number and type of introns in the two mitogenomes indicated that frequent intron loss/gain events occurred during the evolution of Rhizopogon species. Phylogenetic analyses based on Bayesian inference (BI) and Maximum likelihood (ML) methods using a combined mitochondrial gene set yielded identical and well-supported tree topologies, wherein Rhizopogon species showed close relationships with Agaricales species. This is the first study of mitogenomes within the genus Rhizopogon, and it provides a basis for understanding the evolution and differentiation of mitogenomes from the ectomycorrhizal fungal genus.

Comparative transcriptomic analysis identified differentially expressed genes and pathways involved in the interaction between Tremella fuciformis and Annulohypoxylon stygium

Tremella fuciformis is an edible and medicinal white jelly mushroom. It has a life cycle that interacts with its companion fungus Annulohypoxylon stygium, both in natural conditions and artificial cultivation. RNA sequencing (RNA-Seq) was used to study the interaction between T. fuciformis and A. stygium by constructing 5 libraries, including the individual T. fuciformis mycelium (1), the T. fuciformis mycelium after interaction with A. stygium (2), the dual mycelia after interaction (3), the A. stygium mycelium after interaction with T. fuciformis (4), and the individual A. stygium mycelium (5). 33.4 G data and 46,871 Unigenes were generated from de novo splicing. For identification of differentially expressed genes (DEGs) related to interaction, we analyzed the expression data of DEGs1-vs-2 ∩ DEGs1-vs-3, and DEGs5-vs-4 ∩ DEGs5-vs-3. DEGs1-vs-2 ∩ DEGs1-vs-3, and DEGs5-vs-4 ∩ DEGs5-vs-3 data showed 614 DEGs and 1537 DEGs, respectively. The 614 DEGs for T. fuciformis and 1537 DEGs for A. stygium were analyzed by GO annotation and were assigned to biology process, cell composition, and molecular functions. The DEGs were used to match the KEGG database. In T. fuciformis, the pathways are primarily enriched various amino acids metabolism, pentose and glucuronate interconversions. In A. stygium, the pathways are primarily enriched in the biosynthesis of secondary metabolites, biosynthesis of antibiotics, starch and sucrose metabolism. The expression patterns of DEGs determined by qRT-PCR were consistent with those obtained by RNA-Seq, thus validating the reliability of our RNA-Seq data. Future studies of the functions of these interesting genes will be helpful to understand the mechanisms by which T. fuciformis interacts with A. stygium. This will also provide a reference for other research on interacting microorganisms.