The mitochondrial genome of Endoconidiophora resinifera is intron rich

Assemble and comparative analysis of the mitochondrial genome of Rhododendron delavayi: Insights into phylogenetic relationships and genomic variations

Rhododendron delavayi, a notable ornamental plant primarily found in regions of China like Yunnan and Guizhou provinces, holds substantial horticultural value. To elucidate the systematic phylogenetic relationships and organelle genomic differences within R. delavayi and related Rhododendron species, we conducted sequencing and assembly of the complete mitochondrial genome of R. delavayi. The full-length mitochondrial genome of it was a singular circular molecule spanning 1,009,263 bp, comprising 53 protein-coding genes, including 18 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 32 protein-coding genes. A total of 1,182 simple sequence repeats (SSRs) loci were identified in the R. delavayi mitochondrial genome, primarily consisting of single nucleotide, dinucleotide, and trinucleotide repeats. Nucleotide diversity analysis highlighted five genes (atp6, atp9, cox2, nad1, and rpl10) with the highest diversity within the mitochondrial genomes of Rhododendron genus. Comparative analysis of the mitochondrial genome of R. delavayi with those of four other Rhododendron species indicated complex rearrangements in 21 genes, including rps4, nad6, rps3, atp6, cob, atp9, nad7, among others. The mitochondrial phylogenetic tree revealed a close relationship between R. delavayi and R. decorum, forming a sister clade to R. × pulchrum and R. simsii. Furthermore, 126 plastid-to-mitochondrial gene transfers in R. delavayi were identified, ranging from 30 bp to 19,385 bp. These fragments collectively constituted 47.54 % and 9.52 % of the chloroplast and mitochondrial genomes (202,169 bp), respectively. Complex mitochondrial-to-mitochondrial transfers were also observed, with 843 identified fragments totaling 312,036 bp (30.92 % of the mitochondrial genome). Segments exceeding 10 kb may mediate homologous recombination within the mitochondrial molecules. Remarkably, our study underscores that the mitochondrial genome of R. delavayi was the largest reported within the Rhododendron genus to date. The intricate rearrangements observed in the mitochondrial genomes of Rhododendron species, alone with the identification of five potential molecular marker sites, provided valuable insights for species classification and parentage identification within the Rhododendron genus.

The complete mitochondrial genomes of five critical phytopathogenic Bipolaris species: features, evolution, and phylogeny

In the present study, three mitogenomes from the Bipolaris genus (Bipolaris maydis, B. zeicola, and B. oryzae) were assembled and compared with the other two reported Bipolaris mitogenomes (B. oryzae and B. sorokiniana). The five mitogenomes were all circular DNA molecules, with lengths ranging from 106,403 bp to 135,790 bp. The mitogenomes of the five Bipolaris species mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). The PCG length, AT skew and GC skew showed large variability among the 13 PCGs in the five mitogenomes. Across the 13 core PCGs tested, nad6 had the least genetic distance among the 16 Pleosporales species we investigated, indicating that this gene was highly conserved. In addition, the Ka/Ks values for all 12 core PCGs (excluding rps3) were < 1, suggesting that these genes were subject to purifying selection. Comparative mitogenomic analyses indicate that introns were the main factor contributing to the size variation of Bipolaris mitogenomes. The introns of the cox1 gene experienced frequent gain/loss events in Pleosporales species. The gene arrangement and collinearity in the mitogenomes of the five Bipolaris species were almost highly conserved within the genus. Phylogenetic analysis based on combined mitochondrial gene datasets showed that the five Bipolaris species formed well-supported topologies. This study is the first report on the mitogenomes of B. maydis and B. zeicola, as well as the first comparison of mitogenomes among Bipolaris species. The findings of this study will further advance investigations into the population genetics, evolution, and genomics of Bipolaris species.

In silico prediction of polyketide biosynthetic gene clusters in the genomes of Hypericum-borne endophytic fungi

BACKGROUND

The search for new bioactive natural compounds with anticancer activity is still of great importance. Even though their potential for diagnostics and treatment of cancer has already been proved, the availability is still limited. Hypericin, a naphthodianthrone isolated essentially from plant source Hypericum perforatum L. along with other related anthraquinones and bisanthraquinones belongs to this group of compounds. Although it has been proven that hypericin is synthesized by the polyketide pathway in plants, none of the candidate genes coding for key enzymes has been experimentally validated yet. Despite the rare occurrence of anthraquinones in plants, their presence in microorganisms, including endophytic fungi, is quite common. Unlike plants, several biosynthetic genes grouped into clusters (BGCs) in fungal endophytes have already been characterized.

RESULTS

The aim of this work was to predict, identify and characterize the anthraquinone BGCs in de novo assembled and functionally annotated genomes of selected endophytic fungal isolates (Fusarium oxysporum, Plectosphaerella cucumerina, Scedosporium apiospermum, Diaporthe eres, Canariomyces subthermophilus) obtained from different tissues of Hypericum spp. The number of predicted type I polyketide synthase (PKS) BGCs in the studied genomes varied. The non-reducing type I PKS lacking thioesterase domain and adjacent discrete gene encoding protein with product release function were identified only in the genomes of C. subthermophilus and D. eres. A candidate bisanthraquinone BGC was predicted in C. subthermophilus genome and comprised genes coding the enzymes that catalyze formation of the basic anthraquinone skeleton (PKS, metallo-beta-lactamase, decarboxylase, anthrone oxygenase), putative dimerization enzyme (cytochrome P450 monooxygenase), other tailoring enzymes (oxidoreductase, dehydrogenase/reductase), and non-catalytic proteins (fungal transcription factor, transporter protein).

CONCLUSIONS

The results provide an insight into genetic background of anthraquinone biosynthesis in Hypericum-borne endophytes. The predicted bisanthraquinone gene cluster represents a basis for functional validation of the candidate biosynthetic genes in a simple eukaryotic system as a prospective biotechnological alternative for production of hypericin and related bioactive anthraquinones.

Comparative genomic and phylogenetic analyses of mitochondrial genomes of hawthorn (Crataegus spp.) in Northeast China

Hawthorn (Crataegus spp.) plants are major sources of health food and medicines. Twenty species and seven variations of Crataegus are present in China. A variety of unique Crataegus species was found in their natural distribution in northeast China. In the present study, we assembled and annotated the mitochondrial genomes of five Crataegus species from northeastern China. The sizes of the newly sequenced mitochondrial genomes ranged from 245,907 bp to 410,837 bp. A total of 45-55 genes, including 12-19 transfer RNA genes, three ribosomal RNA genes, and 29-33 protein-coding genes (PCGs) were encoded by these mitochondrial genomes. Seven divergent hotspot regions were identified by comparative analyses: atp6, nad3, ccmFN, matR, nad1, nad5, and rps1. The most conserved genes among the Crataegus species, according to the whole-genome correlation analysis, were nad1, matR, nad5, ccmFN, cox1, nad4, trnQ-TTG, trnK-TTT, trnE-TTC, and trnM-CAT. Horizontal gene transfer between organellar genomes was common in Crataegus plants. Based on the phylogenetic trees of mitochondrial PCGs, C. maximowiczii, C. maximowiczii var. ninganensis, and C. bretschneideri shared similar maternal relationships. This study improves Crataegus mitochondrial genome resources and offers important insights into the taxonomy and species identification of this genus.

Morphological, phylogenetic, and genomic evidence reveals the causal agent of thread blight disease of cacao in Peru is a new species of Marasmius in the section Neosessiles, Marasmius infestans sp. nov

The thread blight disease (TBD) of cacao ( Theobroma cacao) in the department of Amazonas, Peru was recently reported to be caused by Marasmius tenuissimus (sect. Neosessiles). This same species is known to be the main causal agent of TBD in West Africa. However, some morphological characteristics, such as the presence of rhizomorphs, the almost exclusively white color, and pileus sizes less than 5 mm, among others, differ to the description of M. tenuissimus. Therefore, we aimed to conduct a taxonomic revision of the cacao-TBD causal agent in Peru, by using thorough micro and macro morphological, phylogenetic, and nuclear and mitochondrial genomic approaches. We showed that the causal agent of TBD of cacao in Amazonas, Peru, belongs to a new species, Marasmius infestans sp. nov. This study enriches our knowledge of species in the sect. Neosessiles, and strongly suggests that the M. tenuissimus species complex is highly diverse.

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.

Mitochondrial genome annotation with MFannot: a critical analysis of gene identification and gene model prediction

Compared to nuclear genomes, mitochondrial genomes (mitogenomes) are small and usually code for only a few dozen genes. Still, identifying genes and their structure can be challenging and time-consuming. Even automated tools for mitochondrial genome annotation often require manual analysis and curation by skilled experts. The most difficult steps are (i) the structural modelling of intron-containing genes; (ii) the identification and delineation of Group I and II introns; and (iii) the identification of moderately conserved, non-coding RNA (ncRNA) genes specifying 5S rRNAs, tmRNAs and RNase P RNAs. Additional challenges arise through genetic code evolution which can redefine the translational identity of both start and stop codons, thus obscuring protein-coding genes. Further, RNA editing can render gene identification difficult, if not impossible, without additional RNA sequence data. Current automated mito- and plastid-genome annotators are limited as they are typically tailored to specific eukaryotic groups. The MFannot annotator we developed is unique in its applicability to a broad taxonomic scope, its accuracy in gene model inference, and its capabilities in intron identification and classification. The pipeline leverages curated profile Hidden Markov Models (HMMs), covariance (CMs) and ERPIN models to better capture evolutionarily conserved signatures in the primary sequence (HMMs and CMs) as well as secondary structure (CMs and ERPIN). Here we formally describe MFannot, which has been available as a web-accessible service (https://megasun.bch.umontreal.ca/apps/mfannot/) to the research community for nearly 16 years. Further, we report its performance on particularly intron-rich mitogenomes and describe ongoing and future developments.

Complete mitochondrial genome assembly and comparison of Camellia sinensis var. Assamica cv. Duntsa

Camellia sinensis var. Assamica cv. Duntsa (C.duntsa), a valuable Theaceae from Hunan Province, has been looked at as a precious tea resource by local farmers because of its economic and ecological value. Genomics study on C.duntsa is essential for the domestication and enhancement of tea tree varieties. In the present study, we used a hybrid approach based on Illumina and PacBio data to sequence and assemble the mitochondrial genome of C.duntsa. The mitochondrial genome of C.duntsa was estimated to be 1,081,996 base pairs (bp) and eighty-one genes consisting of one pseudogene, three ribosomal RNA (rRNA) genes, thirty transfer RNA (tRNA) genes, and forty-seven protein-coding genes (PCGs). Tetramer repetitions made up 43.90% of simple sequence repeats (SSRs). The codon usage bias of the Theaceae mitochondrial gene atp9 was altered by mutation, but the codon usage of other genes was shaped by natural selection. Besides, there are eighteen gene-containing homologous regions between the chloroplast and mitochondrial genomes of C. duntsa.Some genomes including atp8, cox1, cox3, nad7, nad9, rpl16, rpl2, rps19, rps4, and sdh4 are absent in the mitochondrial genome of several Theaceae plant. However, C. duntsa maintains these genes integrity and functionality. Another gene, rps16, is either lacking from the mitochondrial genome of C. duntsa or is present as a pseudogene. C. duntsa and C. sinensis (OM809792) are very similar, as shown by a collinear match across four species of Theaceae; the most conservative genes are nad5, atp9, cox2, rps3, trnA-TGC, trnI-GAT, rrn18, trnV-GAC, and ccmFN. Similarly, the genome's phylogenetic trees revealed that C. duntsa was the sister species to C. sinensis. The results confirmed that the C. duntsa and C. sinensis (OM809792) mitochondrial genome underwent gene rearrangement.In general, our results shows that genomic information from organelles can help us understand plant phylogeny and can also be used to make molecular markers and study how genetic traits change over time. Our research will contribute to the population genetics and evolution of tea plant.

Characterization, comparative phylogenetic, and gene transfer analyses of organelle genomes of Rhododendron × pulchrum

Rhododendron × pulchrum, an important horticultural species, is widely distributed in Europe, Asia, and North America. To analyze the phylogenetic and organelle genome information of R. × pulchrum and its related species, the organelle genome of R. × pulchrum was sequenced and assembled. The complete mitochondrial genome showed lineage DNA molecules, which were 816,410 bp long and contained 64 genes, namely 24 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 37 protein-coding genes. The chloroplast genome of R. × pulchrum was reassembled and re-annotated; the results were different from those of previous studies. There were 42 and 46 simple sequence repeats (SSR) identified from the mitochondrial and chloroplast genomes of R. × pulchrum, respectively. Five genes (nad1, nad2, nad4, nad7, and rps3) were potentially useful molecular markers. The R. × pulchrum mitochondrial genome collinear alignment among five species of the Ericaceae showed that the mitochondrial genomes of these related species have a high degree of homology with R. × pulchrum in this gene region, and the most conservative genes were trnC-GCA, trnD-GUC, trnM-CAU, trnN-GUU, trnY-GUA, atp4, nad4, nad2, nad5, ccmC, and rrn26. The phylogenetic trees of mitochondrial genome showed that R. simsii was a sister to R. × pulchrum. The results verified that there was gene rearrangement between R. × pulchrum and R. simsii mitochondrial genomes. The codon usage bias of 10 Ericaceae mitochondrial genes and 7 Rhododendron chloroplast genes were influenced by mutation, while other genes codon usages had undergone selection. The study identified 13 homologous fragments containing gene sequences between the chloroplast and mitochondrial genomes of R. × pulchrum. Overall, our results illustrate the organelle genome information could explain the phylogenetics of plants and could be used to develop molecular markers and genetic evolution. Our study will facilitate the study of population genetics and evolution in Rhododendron and other genera in Ericaceae.

The first two mitochondrial genomes for the genus Ramaria reveal mitochondrial genome evolution of Ramaria and phylogeny of Basidiomycota

In the present study, we assembled and analyzed the mitogenomes of two Ramaria species. The assembled mitogenomes of Ramaria cfr. rubripermanens and R. rubella were circularized, with sizes of 126,497 bp and 143,271 bp, respectively. Comparative mitogenome analysis showed that intron region contributed the most (contribution rate, 43.74%) to the size variations of Ramaria mitogenomes. The genetic contents, gene length, tRNAs, and codon usages of the two Ramaria mitogenomes varied greatly. In addition, the evolutionary rates of different core protein coding genes (PCGs) in Phallomycetidae mitogenomes varied. We detected large-scale gene rearrangements between Phallomycetidae mitogenomes, including gene displacement and tRNA doubling. A total of 4499 bp and 7746 bp aligned fragments were detected between the mitochondrial and nuclear genomes of R. cfr. rubripermanens and R. rubella, respectively, indicating possible gene transferring events. We further found frequent intron loss/gain and potential intron transfer events in Phallomycetidae mitogenomes during the evolution, and the mitogenomes of R. rubella contained a novel intron P44. Phylogenetic analyses using both Bayesian inference (BI) and Maximum Likelihood (ML) methods based on a combined mitochondrial gene dataset obtained an identical and well-supported phylogenetic tree for Basidiomycota, wherein R. cfr. rubripermanens and Turbinellus floccosus are sister species. This study served as the first report on mitogenomes from the genus Ramaria, which provides a basis for understanding the evolution, genetics, and taxonomy of this important fungal group.

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.

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

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

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.

The strange mitochondrial genomes of Metschnikowia yeasts

While sequencing and characterizing the mitochondrial genomes of 71 strains from the yeast genus Metschnikowia [1] (close cousin to the model species Candida albicans), we uncovered one of the most extreme examples of mitochondrial genome architectural diversity observed to date. These Metschnikowia mitochondrial DNAs (mtDNAs) capture nearly the entire known gene-size and intron-content range for cox1 and cob across all eukaryotic life and show remarkable differences in structure and noncoding content. This genomic variation can be seen both among species and between strains of the same species, raising the question: why are Metschnikowia mitogenomes so malleable?

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.

Mitochondrial genomes of the human pathogens Coccidioides immitis and Coccidioides posadasii

Fungal mitochondrial genomes encode genes involved in crucial cellular processes, such as oxidative phosphorylation and mitochondrial translation, and the molecule has been used as a molecular marker for population genetics studies. Coccidioides immitis and C. posadasii are endemic fungal pathogens that cause coccidioidomycosis in arid regions across both American continents. To date, approximately 150 Coccidioides isolates have been sequenced to infer patterns of variation in nuclear genomes. However, less attention has been given to the mitochondrial genomes of Coccidioides. In this report, we describe the assembly and annotation of mitochondrial reference genomes for two representative strains of C. posadasii and C. immitis, as well as assess population variation among 77 selected genomes. The sizes of the circular-mapping molecules are 68.2 Kb in C. immitis and 75.1 Kb in C. posadasii. We identify 14 mitochondrial protein-coding genes common to most fungal mitochondria, which are largely syntenic across different populations and species of Coccidioides. Both Coccidioides species are characterized by a large number of group I and II introns, harboring twice the number of elements as compared to closely related Onygenales. The introns contain complete or truncated ORFs with high similarity to homing endonucleases of the LAGLIDADG and GIY-YIG families. Phylogenetic comparisons of mitochondrial and nuclear genomes show extensive phylogenetic discordance suggesting that the evolution of the two types of genetic material is not identical. This work represents the first assessment of mitochondrial genomes among isolates of both species of Coccidioides, and provides a foundation for future functional work.

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.

The mitochondrial genome of the grape powdery mildew pathogen Erysiphe necator is intron rich and exhibits a distinct gene organization

Powdery mildews are notorious fungal plant pathogens but only limited information exists on their genomes. Here we present the mitochondrial genome of the grape powdery mildew fungus Erysiphe necator and a high-quality mitochondrial gene annotation generated through cloning and Sanger sequencing of full-length cDNA clones. The E. necator mitochondrial genome consists of a circular DNA sequence of 188,577 bp that harbors a core set of 14 protein-coding genes that are typically present in fungal mitochondrial genomes, along with genes encoding the small and large ribosomal subunits, a ribosomal protein S3, and 25 mitochondrial-encoded transfer RNAs (mt-tRNAs). Interestingly, it also exhibits a distinct gene organization with atypical bicistronic-like expression of the nad4L/nad5 and atp6/nad3 gene pairs, and contains a large number of 70 introns, making it one of the richest in introns mitochondrial genomes among fungi. Sixty-four intronic ORFs were also found, most of which encoded homing endonucleases of the LAGLIDADG or GIY-YIG families. Further comparative analysis of five E. necator isolates revealed 203 polymorphic sites, but only five were located within exons of the core mitochondrial genes. These results provide insights into the organization of mitochondrial genomes of powdery mildews and represent valuable resources for population genetic and evolutionary studies.

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.

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

MOTIVATION

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

RESULTS

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

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.

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.

Panorama of intron dynamics and gene rearrangements in the phylum Basidiomycota as revealed by the complete mitochondrial genome of Turbinellus floccosus

In the present study, the complete mitogenome of Turbinellus floccosus was sequenced, assembled, and compared with other basidiomycete mitogenomes. The mitogenome of T. floccosus consists of a circular DNA molecule, with a size of 62,846 bp. Gene arrangement analysis indicated that large-scale gene rearrangements occurred in the levels of family and genus of basidiomycete species, and the mitogenome of T. floccosus contained a unique gene order. A significant correlation between the number of introns and the mitochondrial genome size of Basidiomycota were detected (P < 0.01). A total of 896 introns were detected in the core protein-coding genes (PCGs) of 74 basidiomycete species, and the cox1 gene was the largest host gene of basidiomycete introns. Intron position class (Pcls) P383 in the cox1 gene was the most common intron in Basidiomycota, which distributed in 40 of 74 basidiomycete species. In addition, frequent intron loss/gain events were detected in basidiomycete species. More than 50% of bases around insertion sites (- 15 bp to 15 bp) of Pcls from different species were conservative, indicating site preferences of intron insertions in Basidiomycota. Further analysis showed that 76.09% of introns tended to insert downstream to a T base in Basidiomycota. Phylogenetic analysis for 74 basidiomycetes indicated mitochondrial genes are effective molecular markers for phylogeny of basidiomycetes. The study served as the first report on the mitogenome from the family Gomphaceae, which will help to understand the intron origin and evolution in Basidiomycota. KEY POINTS: • The mitogenome of Turbinellus floccosus had a unique gene arrangement. • Intron loss/gain events were detected in the 74 basidiomycete species. • Introns tend to insert downstream of a T base in basidiomycete mitogenomes.

Updates and Comparative Analysis of the Mitochondrial Genomes of Paracoccidioides spp. Using Oxford Nanopore MinION Sequencing

The mitochondrial genome of the Paracoccidioides brasiliensis reference isolate Pb18 was first sequenced and described by Cardoso et al. (2007), as a circular genome with a size of 71.3 kb and containing 14 protein coding genes, 25 tRNAs, and the large and small subunits of ribosomal RNA. Later in 2011, Desjardins et al. (2011) obtained partial assemblies of mitochondrial genomes of P. lutzii (Pb01), P. americana (Pb03), and P. brasiliensis sensu stricto (Pb18), although with a size of only 43.1 kb for Pb18. Sequencing errors or other limitations resulting from earlier technologies, and the advantages of NGS (short and long reads), prompted us to improve and update the mtDNA sequences and annotations of two Paracoccidioides species. Using Oxford Nanopore and Illumina read sequencing, we generated high-quality complete de novo mitochondrial genome assemblies and annotations for P. brasiliensis (Pb18) and P. americana (Pb03). Both assemblies were characterized by an unusually long spacer or intron region (>50 kb) between exons 2 and 3 of the nad5 gene, which was moderately conserved between Pb03 and Pb18 but not similar to other reported sequences, except for an unassigned contig in the 2011 assembly of Pb03. The reliability of the insert missing from previous mtDNA genome assemblies was confirmed by inspection of the individual Nanopore read sequences containing nad5 coding DNA, and experimentally by PCR for Pb18. We propose that the insert may aid replication initiation and may be excised to produce a smaller structural variant. The updated mtDNA genomes should enable more accurate SNP and other comparative or evolutionary analyses and primer/probe designs. A comparative analysis of the mtDNA from 32 isolates of Paracoccidioides spp., using the SNPs of the aligned mitochondrial genomes, showed groupings within the brasiliensis species complex that were largely consistent with previous findings from only five mitochondrial loci.

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.

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.

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.

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 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.

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

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