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.

Mitochondrial comparative genomics and phylogenetic signal assessment of mtDNA among arbuscular mycorrhizal fungi.

Mitochondrial (mt) genes, such as cytochrome C oxidase genes (cox), have been widely used for barcoding in many groups of organisms, although this approach has been less powerful in the fungal kingdom due to the rapid evolution of their mt genomes. The use of mt genes in phylogenetic studies of Dikarya has been met with success, while early diverging fungal lineages remain less studied, particularly the arbuscular mycorrhizal fungi (AMF). Advances in next-generation sequencing have substantially increased the number of publically available mtDNA sequences for the Glomeromycota. As a result, comparison of mtDNA across key AMF taxa can now be applied to assess the phylogenetic signal of individual mt coding genes, as well as concatenated subsets of coding genes. Here we show comparative analyses of publically available mt genomes of Glomeromycota, augmented with two mtDNA genomes that were newly sequenced for this study (Rhizophagus irregularis DAOM240159 and Glomus aggregatum DAOM240163), resulting in 16 complete mtDNA datasets. R. irregularis isolate DAOM240159 and G. aggregatum isolate DAOM240163 showed mt genomes measuring 72,293bp and 69,505bp with G+C contents of 37.1% and 37.3%, respectively. We assessed the phylogenies inferred from single mt genes and complete sets of coding genes, which are referred to as "supergenes" (16 concatenated coding genes), using Shimodaira-Hasegawa tests, in order to identify genes that best described AMF phylogeny. We found that rnl, nad5, cox1, and nad2 genes, as well as concatenated subset of these genes, provided phylogenies that were similar to the supergene set. This mitochondrial genomic analysis was also combined with principal coordinate and partitioning analyses, which helped to unravel certain evolutionary relationships in the Rhizophagus genus and for G. aggregatum within the Glomeromycota. We showed evidence to support the position of G. aggregatum within the R. irregularis 'species complex'.

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.

Detection of a transient mitochondrial DNA heteroplasmy in the progeny of crossed genetically divergent isolates of arbuscular mycorrhizal fungi.

Nonself fusion and nuclear genetic exchange have been documented in arbuscular mycorrhizal fungi (AMF), particularly in Rhizophagus irregularis. However, mitochondrial transmission accompanying nonself fusion of genetically divergent isolates remains unknown. Here, we tested the hypothesis that mitochondrial DNA (mtDNA) heteroplasmy occurs in the progeny of spores, obtained by crossing genetically divergent mtDNAs in R. irregularis isolates. Three isolates of geographically distant locations were used to investigate nonself fusions and mtDNA transmission to the progeny. We sequenced two additional mtDNAs of two R. irregularis isolates and developed isolate-specific size-variable markers in intergenic regions of these isolates and those of DAOM-197198. We achieved three crossing combinations in pre-symbiotic and symbiotic phases. Progeny spores per crossing combination were genotyped using isolate-specific markers. We found evidence that nonself recognition occurs between isolates originating from different continents both in pre-symbiotic and symbiotic phases. Genotyping patterns of individual spores from the progeny clearly showed the presence of markers of the two parental mtDNA haplotypes. Our results demonstrate that mtDNA heteroplasmy occurs in the progeny of the crossed isolates. However, this heteroplasmy appears to be a transient stage because all the live progeny spores that were able to germinate showed only one mtDNA haplotype.

Rapid mitochondrial genome evolution through invasion of mobile elements in two closely related species of arbuscular mycorrhizal fungi.

Arbuscular mycorrhizal fungi (AMF) are common and important plant symbionts. They have coenocytic hyphae and form multinucleated spores. The nuclear genome of AMF is polymorphic and its organization is not well understood, which makes the development of reliable molecular markers challenging. In stark contrast, their mitochondrial genome (mtDNA) is homogeneous. To assess the intra- and inter-specific mitochondrial variability in closely related Glomus species, we performed 454 sequencing on total genomic DNA of Glomus sp. isolate DAOM-229456 and we compared its mtDNA with two G. irregulare isolates. We found that the mtDNA of Glomus sp. is homogeneous, identical in gene order and, with respect to the sequences of coding regions, almost identical to G. irregulare. However, certain genomic regions vary substantially, due to insertions/deletions of elements such as introns, mitochondrial plasmid-like DNA polymerase genes and mobile open reading frames. We found no evidence of mitochondrial or cytoplasmic plasmids in Glomus species, and mobile ORFs in Glomus are responsible for the formation of four gene hybrids in atp6, atp9, cox2, and nad3, which are most probably the result of horizontal gene transfer and are expressed at the mRNA level. We found evidence for substantial sequence variation in defined regions of mtDNA, even among closely related isolates with otherwise identical coding gene sequences. This variation makes it possible to design reliable intra- and inter-specific markers.

Group I intron-mediated trans-splicing in mitochondria of Gigaspora rosea and a robust phylogenetic affiliation of arbuscular mycorrhizal fungi with Mortierellales.

Gigaspora rosea is a member of the arbuscular mycorrhizal fungi (AMF; Glomeromycota) and a distant relative of Glomus species that are beneficial to plant growth. To allow for a better understanding of Glomeromycota, we have sequenced the mitochondrial DNA of G. rosea. A comparison with Glomus mitochondrial genomes reveals that Glomeromycota undergo insertion and loss of mitochondrial plasmid-related sequences and exhibit considerable variation in introns. The gene order between the two species is almost completely reshuffled. Furthermore, Gigaspora has fragmented cox1 and rns genes, and an unorthodox initiator tRNA that is tailored to decoding frequent UUG initiation codons. For the fragmented cox1 gene, we provide evidence that its RNA is joined via group I-mediated trans-splicing, whereas rns RNA remains in pieces. According to our model, the two cox1 precursor RNA pieces are brought together by flanking cox1 exon sequences that form a group I intron structure, potentially in conjunction with the nad5 intron 3 sequence. Finally, we present analyses that address the controversial phylogenetic association of Glomeromycota within fungi. According to our results, Glomeromycota are not a separate group of paraphyletic zygomycetes but branch together with Mortierellales, potentially also Harpellales.