Mitochondrial genome rearrangements in glomus species triggered by homologous recombination between distinct mtDNA haplotypes

General features and evolution of mitochondrial genomes in Dictyostelia (Amoebozoa)

Dictyostelia exhibits high diversity; however, mitochondrial genome data remain scarce for many genera. Although key features of some dictyostelid mitogenomes have been identified, several aspects remain unclear, including core gene content, segmental arrangement, and differences between 18S rRNA and mitochondrial gene phylogenies. To address these gaps, we reconstructed two complete mitogenomes-from Cavenderia subdiscoidea and Dictyostelium sp. (TH18CC)-and analyzed mitochondrial genes from ten genera and six additional complete mitogenomes from public databases. A comparison of eight complete mitogenomes revealed a conserved core of 39 protein-coding genes, 17 tRNA genes, and three rRNA genes. Two distinct segmental arrangements were identified: Dictyostelium (except D. purpureum) exhibits an A-C-B pattern, while other genera display an A-B-C pattern defined by the clusters nad9-atp1, trnC(GCA)-atp9, and rnl-nad3. Phylogenetic analyses based on 18S rDNA and mitochondrial rns suggest the transposition between segments B and C occurred after D. purpureum diverged from other Dictyostelium species, potentially involving tRNA gene displacement. In contrast, the mitochondrial protein-coding gene phylogeny differs from the rRNA trees, indicating that these gene sets may have evolved independently. These findings advance our understanding of dictyostelid mitogenome structure and evolution.

From Chaos Comes Order: Genetics and Genome Biology of Arbuscular Mycorrhizal Fungi

Arbuscular mycorrhizal fungi (AMF) are obligate mutualists that can enhance nutrition and growth of their plant hosts while providing protection against pathogens. AMF produce spores and hyphal networks that can carry thousands of nuclei in a continuous cytoplasm, with no evidence of sexual reproduction. This review examines the impact of genomic technologies on our view of AMF genetics and evolution. We highlight how the genetics, nuclear dynamics, and epigenetics of these prominent symbionts follow trends preserved in distant multinucleate fungal relatives. We also propose new avenues of research to improve our understanding of their nuclear biology and their intricate genetic interactions with plant hosts.

Characterization of the complete mitochondrial genome of Mucor indicus Lendn. 1930 (Mucorales: Mucoraceae), isolated from the wine fermentation system

Mucor indicus Lendn. 1930 has been widely used in food fermentation; however, its mitochondrial genome characteristics are not well understood. In this study, the complete mitochondrial genome of M. indicus was obtained, which was 61,400 bp in length with a GC content of 33%. The M. indicus mitochondrial genome was found to contain 14 core protein-coding genes, four free-standing ORFs, 18 intronic ORFs, 26 tRNAs, and two rRNA genes. Phylogenetic trees were generated for 25 early-differentiated fungi using the Bayesian inference (BI) method, which demonstrated that M. indicus is closely related to Mucor piriformis. This study provides useful information for the classification and evolution of Mucor species or other early-differentiated fungi.

Mitochondrial genomic characteristics and phylogenetic analysis of a brewing fungus, Rhizopus microsporus Tiegh. 1875 (Mucorales: Rhizopodaceae)

Rhizopus microsporus Tiegh. 1875 is widely used in a variety of industries, such as brewing, wine making, baking, and medicine production, as it has the capability to break down proteins and generate surface-active agents. To date, the mitochondrial genome features of early evolved fungi from the Rhizopus genus have not been extensively studied. Our research obtained a full mitochondrial genome of R. microsporus species, which was 43,837 bp in size and had a GC content of 24.93%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 7 intronic ORFs, 24 tRNAs, and 2 rRNA genes. Through the use of the BI phylogenetic inference method, we were able to create phylogenetic trees for 25 early differentiation fungi which strongly supported the major clades; this indicated that R. microsporus is most closely related to Rhizopus oryzae.

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 complete mitochondrial genome of Linnemannia amoeboidea (W. Gams) Vandepol & Bonito (Mortierellales: Mortierellaceae)

The complete mitochondrial genome of Linnemannia amoeboidea (W. Gams) Vandepol & Bonito 2020 (Strain no.: CBS 889.72) was sequenced under the next-generation sequencing platform. It was the second one in the family Mortierellaceae Luerss. 1877. The circular genome was 49,702 bp in size, with a GC content of 20.86%. Gene prediction revealed 15 PCGs, two rRNA genes, 26 tRNA genes, one rnpB gene and seven ORFs. Phylogenetic analyses showed that L. amoeboidea was closely related to Podila verticillate (Linnem.) Vandepol & Bonito 2020.

Tandem integration of circular plasmid contributes significantly to the expanded mitochondrial genomes of the green-tide forming alga Ulva meridionalis (Ulvophyceae, Chlorophyta)

Comparative mitogenomics of Ulva species have revealed remarkable variations in genome size due to the integration of exogenous DNA fragments, the proliferation of group I/II introns, and the change of repeat sequences. The genus Ulva is a species-rich taxonomic group, containing a variety of green-tide forming algae. In this study, five complete mitogenomes of the green-tide forming macroalga, Ulva meridionalis R. Horimoto and S. Shimada, were assembled and compared with the available ulvophyceae mtDNAs. The main circular mitogenomes of U. meridionalis ranged from 82.94 to 111.49 kb in size, and its 111.49-kb mitogenome was the largest Ulva mitogenome sequenced so far. The expansion of U. meridionalis mitogenomes is mainly due to the tandem integration of a 5.36-kb mitochondrial circular plasmid (pUme), as well as the proliferation of introns. An intact DNA-directed RNA polymerase gene (rpo) was present in pUme of U. meridionalis and was then detected in two putative plasmids (pUmu1 and pUmu2) found in Ulva mutabilis. The observed integration of the circular plasmid into U. meridionalis mitogenomes seems to occur via homologous recombination, and is a more recent evolutionary event. Many highly homologous sequences of these three putative plasmids can be detected in the other Ulva mtDNAs sequenced thus far, indicating the integration of different mitochondrial plasmid DNA into the mitogenomes is a common phenomenon in the evolution of Ulva mitogenomes. The random incidence of destruction of plasmid-derived rpos and open reading frames (orfs) suggests that their existence is not the original characteristic of Ulva mitogenomes and there is no selective pressure to maintain their integrity. The frequent integration and rapid divergence of plasmid-derived sequences is one of the most important evolutionary forces to shape the diversity of Ulva mitogenomes.

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

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

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

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.

Object Weighting: A New Clustering Approach to Deal with Outliers and Cluster Overlap in Computational Biology

Considerable efforts have been made over the last decades to improve the robustness of clustering algorithms against noise features and outliers, known to be important sources of error in clustering. Outliers dominate the sum-of-the-squares calculations and generate cluster overlap, thus leading to unreliable clustering results. They can be particularly detrimental in computational biology, e.g., when determining the number of clusters in gene expression data related to cancer or when inferring phylogenetic trees and networks. While the issue of feature weighting has been studied in detail, no clustering methods using object weighting have been proposed yet. Here we describe a new general data partitioning method that includes an object-weighting step to assign higher weights to outliers and objects that cause cluster overlap. Different object weighting schemes, based on the Silhouette cluster validity index, the median and two intercluster distances, are defined. We compare our novel technique to a number of popular and efficient clustering algorithms, such as K-means, X-means, DAPC and Prediction Strength. In the presence of outliers and cluster overlap, our method largely outperforms X-means, DAPC and Prediction Strength as well as the K-means algorithm based on feature weighting.

Repeats in mitochondrial and chloroplast genomes characterize the ecotypes of the Oryza

Mitochondria and chloroplast are very important organelles for organism, participating in basic life activity. Their genomes contain many repeats which can lead to a variation of genome structure. Oryza is an important genus for human beings' nutrition. Several mitochondrial and chloroplast genomes of Oryza have been sequenced, which help us to insight the distribution and evolution of the repeats in Oryza species. In this paper, we compared six mitochondrial and 13 chloroplast genomes of Oryza and found that the structures of mitochondrial genomes were more diverse than chloroplast genomes. Since repeats can change the structure of the genome, resulting in the structural diversity of the genome, we analyzed all repeats and found 31 repeats in mitochondrial and 13 repeats in chloroplast genomes. Further, we developed 21 pairs of MRS molecular markers and 12 pairs of CRS molecular markers based on mitochondrial repeats and chloroplast repeats, respectively. These molecular markers can be used to detect the repeat-mediated recombination in Oryza mitochondrial and chloroplast genomes by PCR or fluorescence quantification.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-020-01198-6.

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.

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.

Intra-specific comparison of mitochondrial genomes reveals host gene fragment exchange via intron mobility in Tremella fuciformis

Background

Mitochondrial genomic sequences are known to be variable. Comparative analyses of mitochondrial genomes can reveal the nature and extent of their variation.

Results

Draft mitochondrial genomes of 16 Tremella fuciformis isolates (TF01-TF16) were assembled from Illumina and PacBio sequencing data. Mitochondrial DNA contigs were extracted and assembled into complete circular molecules, ranging from 35,104 bp to 49,044 bp in size. All mtDNAs contained the same set of 41 conserved genes with identical gene order. Comparative analyses revealed that introns and intergenic regions were variable, whereas genic regions (including coding sequences, tRNA, and rRNA genes) were conserved. Among 24 introns detected, 11 were in protein-coding genes, 3 in tRNA genes, and the other 10 in rRNA genes. In addition, two mobile fragments were found in intergenic regions. Interestingly, six introns containing N-terminal duplication of the host genes were found in five conserved protein-coding gene sequences. Comparison of genes with and without these introns gave rise to the following proposed model: gene fragment exchange with other species can occur via gain or loss of introns with N-terminal duplication of the host genes.

Conclusions

Our findings suggest a novel mechanism of fungal mitochondrial gene evolution: partial foreign gene replacement though intron mobility.

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.

High and Variable Rates of Repeat-Mediated Mitochondrial Genome Rearrangement in a Genus of Plants

For 30 years, it has been clear that angiosperm mitochondrial genomes evolve rapidly in sequence arrangement (i.e., synteny), yet absolute rates of rearrangement have not been measured in any plant group, nor is it known how much these rates vary. To investigate these issues, we sequenced and reconstructed the rearrangement history of seven mitochondrial genomes in Monsonia (Geraniaceae). We show that rearrangements (occurring mostly as inversions) not only take place at generally high rates in these genomes but also uncover significant variation in rearrangement rates. For example, the hyperactive mitochondrial genome of Monsonia ciliata has accumulated at least 30 rearrangements over the last million years, whereas the branch leading to M. ciliata and its sister species has sustained rearrangement at a rate that is at least ten times lower. Furthermore, our analysis of published data shows that rates of mitochondrial genome rearrangement in seed plants vary by at least 600-fold. We find that sites of rearrangement are highly preferentially located in very close proximity to repeated sequences in Monsonia. This provides strong support for the hypothesis that rearrangement in angiosperm mitochondrial genomes occurs largely through repeat-mediated recombination. Because there is little variation in the amount of repeat sequence among Monsonia genomes, the variable rates of rearrangement in Monsonia probably reflect variable rates of mitochondrial recombination itself. Finally, we show that mitochondrial synonymous substitutions occur in a clock-like manner in Monsonia; rates of mitochondrial substitutions and rearrangements are therefore highly uncoupled in this group.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Comparative mitochondrial genomics of cryptophyte algae: gene shuffling and dynamic mobile genetic elements

BACKGROUND

Cryptophytes are an ecologically important group of algae comprised of phototrophic, heterotrophic and osmotrophic species. This lineage is of great interest to evolutionary biologists because their plastids are of red algal secondary endosymbiotic origin. Cryptophytes have a clear phylogenetic affinity to heterotrophic eukaryotes and possess four genomes: host-derived nuclear and mitochondrial genomes, and plastid and nucleomorph genomes of endosymbiotic origin.

RESULTS

To gain insight into cryptophyte mitochondrial genome evolution, we sequenced the mitochondrial DNAs of five species and performed a comparative analysis of seven genomes from the following cryptophyte genera: Chroomonas, Cryptomonas, Hemiselmis, Proteomonas, Rhodomonas, Storeatula and Teleaulax. The mitochondrial genomes were similar in terms of their general architecture, gene content and presence of a large repeat region. However, gene order was poorly conserved. Characteristic features of cryptophyte mtDNAs included large syntenic clusters resembling α-proteobacterial operons that encode bacteria-like rRNAs, tRNAs, and ribosomal protein genes. The cryptophyte mitochondrial genomes retain almost all genes found in many other eukaryotes including the nad, sdh, cox, cob, and atp genes, with the exception of sdh2 and atp3. In addition, gene cluster analysis showed that cryptophytes possess a gene order closely resembling the jakobid flagellates Jakoba and Reclinomonas. Interestingly, the cox1 gene of R. salina, T. amphioxeia, and Storeatula species was found to contain group II introns encoding a reverse transcriptase protein, as did the cob gene of Storeatula species CCMP1868.

CONCLUSIONS

These newly sequenced genomes increase the breadth of data available from algae and will aid in the identification of general trends in mitochondrial genome evolution. While most of the genomes were highly conserved, extensive gene arrangements have shuffled gene order, perhaps due to genome rearrangements associated with hairpin-containing mobile genetic elements, tRNAs with palindromic sequences, and tandem repeat sequences. The cox1 and cob gene sequences suggest that introns have recently been acquired during cryptophyte evolution. Comparison of phylogenetic trees based on plastid and mitochondrial genome data sets underscore the different evolutionary histories of the host and endosymbiont components of present-day cryptophytes.

Mitochondrial Recombination and Introgression during Speciation by Hybridization

Genome recombination is a major source of genotypic diversity and contributes to adaptation and speciation following interspecies hybridization. The contribution of recombination in these processes has been thought to be largely limited to the nuclear genome because organelles are mostly uniparentally inherited in animals and plants, which prevents recombination. Unicellular eukaryotes such as budding yeasts do, however, transmit mitochondria biparentally, suggesting that during hybridization, both parents could provide alleles that contribute to mitochondrial functions such as respiration and metabolism in hybrid populations or hybrid species. We examined the dynamics of mitochondrial genome transmission and evolution during speciation by hybridization in the natural budding yeast Saccharomyces paradoxus. Using population-scale mitochondrial genome sequencing in two endemic North American incipient species SpB and SpC and their hybrid species SpC*, we found that both parental species contributed to the hybrid mitochondrial genome through recombination. We support our findings by showing that mitochondrial recombination between parental types is frequent in experimental crosses that recreate the early step of this speciation event. In these artificial hybrids, we observed that mitochondrial genome recombination enhances phenotypic variation among diploid hybrids, suggesting that it could play a role in the phenotypic differentiation of hybrid species. Like the nuclear genome, the mitochondrial genome can, therefore, also play a role in hybrid speciation.

Real-time PCR quantification of arbuscular mycorrhizal fungi: does the use of nuclear or mitochondrial markers make a difference?

Root colonization by arbuscular mycorrhizal fungi (AMF) can be quantified by different approaches. We compared two approaches that enable discrimination of specific AMF taxa and are therefore emerging as alternative to most commonly performed microscopic quantification of AMF in roots: quantitative real-time PCR (qPCR) using markers in nuclear ribosomal DNA (nrDNA) and mitochondrial ribosomal DNA (mtDNA). In a greenhouse experiment, Medicago truncatula was inoculated with four isolates belonging to different AMF species (Rhizophagus irregularis, Claroideoglomus claroideum, Gigaspora margarita and Funneliformis mosseae). The AMF were quantified in the root samples by qPCR targeted to both markers, microscopy and contents of AMF-specific phospholipid fatty acids (PLFA). Copy numbers of nrDNA and mtDNA were closely related within all isolates; however, the slopes and intercepts of the linear relationships significantly differed among the isolates. Across all isolates, a large proportion of variance in nrDNA copy numbers was explained by root colonization intensity or contents of AMF-specific PLFA, while variance in mtDNA copy numbers was mainly explained by differences among AMF isolates. We propose that the encountered inter-isolate differences in the ratios of mtDNA and nrDNA copy numbers reflect different physiological states of the isolates. Our results suggest that nrDNA is a more suitable marker region than mtDNA for the quantification of multiple AMF taxa as its copy numbers are better related to fungal biomass across taxa than are copy numbers of mtDNA.

Frequent heteroplasmy and recombination in the mitochondrial genomes of the basidiomycete mushroom Thelephora ganbajun

In the majority of sexual eukaryotes, the mitochondrial genomes are inherited uniparentally. As a result, individual organisms are homoplasmic, containing mitochondrial DNA (mtDNA) from a single parent. Here we analyzed the mitochondrial genotypes in Clade I of the gourmet mushroom Thelephora ganbajun from its broad geographic distribution range. A total of 299 isolates from 28 geographic locations were sequenced at three mitochondrial loci: the mitochondrial small ribosomal RNA gene, and the cytochrome c oxidase subunits I (COX1) and III (COX3) genes. Quantitative PCR analyses showed that the strains had about 60–160 copies of mitochondrial genomes per cell. Interestingly, while no evidence of heteroplasmy was found at the 12S rRNA gene, 262 of the 299 isolates had clear evidence of heterogeneity at either the COX1 (261 isolates) or COX3 (12 isolates) gene fragments. The COX1 heteroplasmy was characterized by two types of introns residing at different sites of the same region and at different frequencies among the isolates. Allelic association analyses of the observed mitochondrial polymorphic nucleotide sites suggest that mtDNA recombination is common in natural populations of this fungus. Our results contrast the prevailing view that heteroplasmy, if exists, is only transient in basidiomycete fungi.

The fate of mitochondria after infection of the Mucoralean fungus Absidia glauca by the fusion parasite Parasitella parasitica: comparison of mitochondrial genomes in zygomycetes

Absidia glauca and Parasitella parasitica constitute a versatile experimental system for studying horizontal gene transfer between a mucoralean host and its fusion parasite. The A. glauca chondriome has a length of approximately 63 kb and a GC content of 28%. The chondriome of P. parasitica is larger, 83 kb, and contains 31% GC base pairs. These mtDNAs contain the standard fungal mitochondrial gene set, small and large subunit rRNAs, plus ribonuclease P RNA. Comparing zygomycete chondriomes reveals an unusually high number of homing endonuclease genes in P. parasitica, substantiating the mobility of intron elements independent of host-parasite interactions.

Molecular diagnostic toolkit for Rhizophagus irregularis isolate DAOM-197198 using quantitative PCR assay targeting the mitochondrial genome

Rhizophagus irregularis (previously named Glomus irregulare) is one of the most widespread and common arbuscular mycorrhizal fungal (AMF) species. It has been recovered worldwide in agricultural and natural soils, and the isolate DAOM-197198 has been utilized as a commercial inoculant for two decades. Despite the ecological and economical importance of this taxon, specific markers for quantification of propagules by quantitative real-time PCR (qPCR) are extremely limited and none have been rigorously validated for quality control of manufactured products such as biofertilizers. From the sequencing of 14 complete AMF mitochondrial (mt) genomes, a qPCR assay using a hydrolysis probe designed in the single copy cox3-rnl intergenic region was tested and validated to specifically and accurately quantify the spores of R. irregularis isolate DAOM-197198. Specificity tests were performed using standard PCR and qPCR, and results clearly showed that the primers specifically amplified the isolate DAOM-197198, yielding a PCR product of 106 bp. According to the qPCR analyses on spores produced in vitro, the average copy number of mt genomes per spore was 3172 ± 304 SE (n = 6). Quantification assays were successfully undertaken on known and unknown samples in liquid suspensions and commercial dry formulations to show the accuracy, precision, robustness, and reproducibility of the qPCR assay. This study provides a powerful molecular toolkit specifically designed to quantify spores of the model AMF isolate DAOM-197198. The approach of molecular toolkit used in our study could be applied to other AMF taxa and will be useful to research institutions and governmental and industrial laboratories running routine quality control of AMF-based products.

The large (134.9 kb) mitochondrial genome of the glomeromycete Funneliformis mosseae

Funneliformis mosseae is among the most ecologically and economically important glomeromycete species and occurs both in natural and disturbed areas in a wide range of habitats and climates. In this study, we report the sequencing of the complete mitochondrial (mt) genome of F. mosseae isolate FL299 using 454 pyrosequencing and Illumina HiSeq technologies. This mt genome is a full-length circular chromosome of 134,925 bp, placing it among the largest mitochondrial DNAs (mtDNAs) in the fungal kingdom. A comparative analysis with publically available arbuscular mycorrhizal fungal mtDNAs revealed that the mtDNA of F. mosseae FL299 contained a very large number of insertions contributing to its expansion. The gene synteny was completely reshuffled compared to previously published glomeromycotan mtDNAs and several genes were oriented in an anti-sense direction. Furthermore, the presence of different types of introns and insertions in rnl (14 introns) made this gene very distinctive in Glomeromycota. The presence of alternative genetic codes in both initiation (GUG) and termination (UGA) codons was another new feature in this mtDNA compared to previously published glomeromycotan mt genomes. The phylogenetic analysis inferred from the analysis of 14 protein mt genes confirmed the position of the Glomeromycota clade as a sister group of Mortierellomycotina. This mt genome is the largest observed so far in Glomeromycota and the first mt genome within the Funneliformis clade, providing new opportunities to better understand their evolution and to develop molecular markers.

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

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

Group II Intron-Mediated Trans-Splicing in the Gene-Rich Mitochondrial Genome of an Enigmatic Eukaryote, Diphylleia rotans

Although mitochondria have evolved from a single endosymbiotic event, present day mitochondria of diverse eukaryotes display a great range of genome structures, content and features. Group I and group II introns are two features that are distributed broadly but patchily in mitochondrial genomes across branches of the tree of eukaryotes. While group I intron-mediated trans-splicing has been reported from some lineages distantly related to each other, findings of group II intron-mediated trans-splicing has been restricted to members of the Chloroplastida. In this study, we found the mitochondrial genome of the unicellular eukaryote Diphylleia rotans possesses currently the second largest gene repertoire. On the basis of a probable phylogenetic position of Diphylleia, which is located within Amorphea, current mosaic gene distribution in Amorphea must invoke parallel gene losses from mitochondrial genomes during evolution. Most notably, although the cytochrome c oxidase subunit (cox) 1 gene was split into four pieces which located at a distance to each other, we confirmed that a single mature mRNA that covered the entire coding region could be generated by group II intron-mediated trans-splicing. This is the first example of group II intron-mediated trans-splicing outside Chloroplastida. Similar trans-splicing mechanisms likely work for bipartitely split cox2 and nad3 genes to generate single mature mRNAs. We finally discuss origin and evolution of this type of trans-splicing in D. rotans as well as in eukaryotes.

Independent mitochondrial and nuclear exchanges arising in Rhizophagus irregularis crossed-isolates support the presence of a mitochondrial segregation mechanism

BACKGROUND

Arbuscular mycorrhizal fungi (AMF) are members of the phylum Glomeromycota, an early divergent fungal lineage that forms symbiotic associations with the large majority of land plants. These organisms are asexual obligate biotrophs, meaning that they cannot complete their life cycle in the absence of a suitable host. These fungi can exchange genetic information through hyphal fusions (i.e. anastomosis) with genetically compatible isolates belonging to the same species. The occurrence of transient mitochondrial length-heteroplasmy through anastomosis between geographically distant Rhizophagus irregularis isolates was previously demonstrated in single spores resulting from crossing experiments. However, (1) the persistence of this phenomenon in monosporal culture lines from crossed parental isolates, (2) its correlation with nuclear exchanges and (3) the potential mechanisms responsible for mitochondrial inheritance are still unknown. Using the AMF model organism R. irregularis, we tested whether the presence of a heteroplasmic state in progeny spores was linked to the occurrence of nuclear exchanges and whether the previously observed heteroplasmic state persisted in monosporal in vitro crossed-culture lines. We also investigated the presence of a putative mitochondrial segregation apparatus in Glomeromycota by identifying proteins similar to those found in other fungal groups.

RESULTS

We observed the occurrence of biparental inheritance both for mitochondrial and nuclear markers tested in single spores obtained from crossed-isolates. However, only one parental mitochondrial DNA and nuclear genotype were recovered in each monosporal crossed-cultures, with an overrepresentation of certain mitochondrial haplotypes. These results strongly support the presence of a nuclear-independent mitochondrial segregation mechanism in R. irregularis. Furthermore, a nearly complete set of genes was identified with putative orthology to those found in other fungi and known to be associated with the mitochondrial segregation in Saccharomyces cerevisiae and filamentous fungi.

CONCLUSIONS

Our findings suggest that mitochondrial segregation might take place either during spore formation or colony development and that it might be independent of the nuclear segregation machinery. We present the basic building blocks for a better understanding of the mitochondrial inheritance process and segregation in these important symbiotic fungi. The comprehension of these processes is of great importance since it has been shown that different segregated lines of the same isolate can have variable effects on the host plant.

Genome diversity in arbuscular mycorrhizal fungi

Arbuscular mycorrhizal fungi (Glomeromycota) are the most widespread and important symbionts of plants. They cannot be cultured without plants, are apparently asexual, and have multiple nuclei in a common cytoplasm. There is evidence for genetic variation among nuclei, and for segregation of this variation during growth, but these findings remain contentious. Recently, two papers have reported whole genome sequences for a strain of Rhizophagus irregularis; both suggest that genetic variation among nuclei is low. Genome assembly is very incomplete, though, so significant nuclear diversity cannot be excluded. While the diversity of nuclear genomes remains unresolved, multiple complete mitochondrial genomes are now available; there is virtually no variation within isolates, but significant variation between them.

The mitochondrial genome of the glomeromycete Rhizophagus sp. DAOM 213198 reveals an unusual organization consisting of two circular chromosomes

Mitochondrial (mt) genomes are intensively studied in Ascomycota and Basidiomycota, but they are poorly documented in basal fungal lineages. In this study, we sequenced the complete mtDNA of Rhizophagus sp. DAOM 213198, a close relative to Rhizophagus irregularis, a widespread, ecologically and economical relevant species belonging to Glomeromycota. Unlike all other known taxonomically close relatives harboring a full-length circular chromosome, mtDNA of Rhizophagus sp. reveals an unusual organization with two circular chromosomes of 61,964 and 29,078 bp. The large chromosome contained nine protein-coding genes (atp9, nad5, cob, nad4, nad1, nad4L, cox1, cox2, and atp8), small subunit rRNA gene (rns), and harbored 20 tRNA-coding genes and 10 orfs, while the small chromosome contained five protein-coding genes (atp6, nad2, nad3, nad6, and cox3), large subunit rRNA gene (rnl) in addition to 5 tRNA-coding genes, and 8 plasmid-related DNA polymerases (dpo). Although structural variation of plant mt genomes is well documented, this study is the first report of the presence of two circular mt genomes in arbuscular mycorrhizal fungi. Interestingly, the presence of dpo at the breakage point in intergenes cox1-cox2 and rnl-atp6 for large and small mtDNAs, respectively, could be responsible for the conversion of Rhizophagus sp. mtDNA into two chromosomes. Using quantitative real-time polymerase chain reaction, we found that both mtDNAs have an equal abundance. This study reports a novel mtDNA organization in Glomeromycota and highlights the importance of studying early divergent fungal lineages to describe novel evolutionary pathways in the fungal kingdom.

Intraisolate mitochondrial genetic polymorphism and gene variants coexpression in arbuscular mycorrhizal fungi

Arbuscular mycorrhizal fungi (AMF) are multinucleated and coenocytic organisms, in which the extent of the intraisolate nuclear genetic variation has been a source of debate. Conversely, their mitochondrial genomes (mtDNAs) have appeared to be homogeneous within isolates in all next generation sequencing (NGS)-based studies. Although several lines of evidence have challenged mtDNA homogeneity in AMF, extensive survey to investigate intraisolate allelic diversity has not previously been undertaken. In this study, we used a conventional polymerase chain reaction -based approach on selected mitochondrial regions with a high-fidelity DNA polymerase, followed by cloning and Sanger sequencing. Two isolates of Rhizophagus irregularis were used, one cultivated in vitro for several generations (DAOM-197198) and the other recently isolated from the field (DAOM-242422). At different loci in both isolates, we found intraisolate allelic variation within the mtDNA and in a single copy nuclear marker, which highlighted the presence of several nonsynonymous mutations in protein coding genes. We confirmed that some of this variation persisted in the transcriptome, giving rise to at least four distinct nad4 transcripts in DAOM-197198. We also detected the presence of numerous mitochondrial DNA copies within nuclear genomes (numts), providing insights to understand this important evolutionary process in AMF. Our study reveals that genetic variation in Glomeromycota is higher than what had been previously assumed and also suggests that it could have been grossly underestimated in most NGS-based AMF studies, both in mitochondrial and nuclear genomes, due to the presence of low-level mutations.

Sequencing and comparison of the mitochondrial COI gene from isolates of Arbuscular Mycorrhizal Fungi belonging to Gigasporaceae and Glomeraceae families

Arbuscular Mycorrhizal Fungi (AMF) are well known for their ecological importance and their positive influence on plants. The genetics and phylogeny of this group of fungi have long been debated. Nuclear markers are the main tools used for phylogenetic analyses, but they have sometimes proved difficult to use because of their extreme variability. Therefore, the attention of researchers has been moving towards other genomic markers, in particular those from the mitochondrial DNA. In this study, 46 sequences of different AMF isolates belonging to two main clades Gigasporaceae and Glomeraceae have been obtained from the mitochondrial gene coding for the Cytochrome c Oxidase I (COI), representing the largest dataset to date of AMF COI sequences. A very low level of divergence was recorded in the COI sequences from the Gigasporaceae, which could reflect either a slow rate of evolution or a more recent evolutionary divergence of this group. On the other hand, the COI sequence divergence between Gigasporaceae and Glomeraceae was high, with synonymous divergence reaching saturated levels. This work also showed the difficulty in developing valuable mitochondrial markers able to effectively distinguish all Glomeromycota species, especially those belonging to Gigasporaceae, yet it represents a first step towards the development of a full mtDNA-based dataset which can be used for further phylogenetic investigations of this fungal phylum.