High variability of mitochondrial gene order among fungi

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

Mitochondrial genome evolution in the Saccharomyces sensu stricto complex

Exploring the evolutionary patterns of mitochondrial genomes is important for our understanding of the Saccharomyces sensu stricto (SSS) group, which is a model system for genomic evolution and ecological analysis. In this study, we first obtained the complete mitochondrial sequences of two important species, Saccharomyces mikatae and Saccharomyces kudriavzevii. We then compared the mitochondrial genomes in the SSS group with those of close relatives, and found that the non-coding regions evolved rapidly, including dramatic expansion of intergenic regions, fast evolution of introns and almost 20-fold higher rearrangement rates than those of the nuclear genomes. However, the coding regions, and especially the protein-coding genes, are more conserved than those in the nuclear genomes of the SSS group. The different evolutionary patterns of coding and non-coding regions in the mitochondrial and nuclear genomes may be related to the origin of the aerobic fermentation lifestyle in this group. Our analysis thus provides novel insights into the evolution of mitochondrial genomes.

SMRT Sequencing Revealed Mitogenome Characteristics and Mitogenome-Wide DNA Modification Pattern in Ophiocordyceps sinensis

Single molecule, real-time (SMRT) sequencing was used to characterize mitochondrial (mt) genome of Ophiocordyceps sinensis and to analyze the mt genome-wide pattern of epigenetic DNA modification. The complete mt genome of O. sinensis, with a size of 157,539 bp, is the fourth largest Ascomycota mt genome sequenced to date. It contained 14 conserved protein-coding genes (PCGs), 1 intronic protein rps3, 27 tRNAs and 2 rRNA subunits, which are common characteristics of the known mt genomes in Hypocreales. A phylogenetic tree inferred from 14 PCGs in Pezizomycotina fungi supports O. sinensis as most closely related to Hirsutella rhossiliensis in Ophiocordycipitaceae. A total of 36 sequence sites in rps3 were under positive selection, with dN/dS >1 in the 20 compared fungi. Among them, 16 sites were statistically significant. In addition, the mt genome-wide base modification pattern of O. sinensis was determined in this study, especially DNA methylation. The methylations were located in coding and uncoding regions of mt PCGs in O. sinensis, and might be closely related to the expression of PCGs or the binding affinity of transcription factor A to mtDNA. Consequently, these methylations may affect the enzymatic activity of oxidative phosphorylation and then the mt respiratory rate; or they may influence mt biogenesis. Therefore, methylations in the mitogenome of O. sinensis might be a genetic feature to adapt to the cold and low PO2 environment at high altitude, where O. sinensis is endemic. This is the first report on epigenetic modifications in a fungal mt genome.

Mitochondrial genome of the nematode endoparasitic fungus Hirsutella vermicola reveals a high level of synteny in the family Ophiocordycipitaceae

Ophiocordycipitaceae is a diverse fungal family comprising multiple ecologically, economically, medicinally, and culturally important fungal species; however, only four species of the family have available mitochondrial genomes (mitogenomes). In this study, the complete mitogenome of the nematode endoparasitic fungus Hirsutella vermicola in Ophiocordycipitaceae was sequenced, and a comparative mitogenomic analysis of Ophiocordycipitaceae was performed. We found that the 53,793-bp circular mitogenome of H. vermicola, except for standard fungal mitochondrial genes, harbors seven introns acquired possibly through lateral transfer from other fungi and three free-standing open reading frames (ORFs) coding for hypothetical proteins. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed its placement in Ophiocordycipitaceae. Comparison on five mitogenomes of Ophiocordycipitaceae revealed great variation on their sizes, from 35.2 kb in Tolypocladium ophioglossoides to 157.5 kb in Ophiocordyceps sinensis, mainly due to variable numbers of introns (from 7 to 54) as well as variable lengths of intergenic regions. The five mitogenomes, however, are highly syntenic to each other in terms of gene order, the presence of an intronic ORF encoding ribosomal protein S3 within rnl, and the nad2/nad3 joining pattern. Our study is the first report of the mitogenome of H. vermicola and has facilitated the understanding of mitogenome evolution of Ophiocordycipitaceae.

Elevated Rate of Genome Rearrangements in Radiation-Resistant Bacteria

A number of bacterial, archaeal, and eukaryotic species are known for their resistance to ionizing radiation. One of the challenges these species face is a potent environmental source of DNA double-strand breaks, potential drivers of genome structure evolution. Efficient and accurate DNA double-strand break repair systems have been demonstrated in several unrelated radiation-resistant species and are putative adaptations to the DNA damaging environment. Such adaptations are expected to compensate for the genome-destabilizing effect of environmental DNA damage and may be expected to result in a more conserved gene order in radiation-resistant species. However, here we show that rates of genome rearrangements, measured as loss of gene order conservation with time, are higher in radiation-resistant species in multiple, phylogenetically independent groups of bacteria. Comparison of indicators of selection for genome organization between radiation-resistant and phylogenetically matched, nonresistant species argues against tolerance to disruption of genome structure as a strategy for radiation resistance. Interestingly, an important mechanism affecting genome rearrangements in prokaryotes, the symmetrical inversions around the origin of DNA replication, shapes genome structure of both radiation-resistant and nonresistant species. In conclusion, the opposing effects of environmental DNA damage and DNA repair result in elevated rates of genome rearrangements in radiation-resistant bacteria.

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.

Two complete mitochondrial genomes from Praticolella mexicana Perez, 2011 (Polygyridae) and gene order evolution in Helicoidea (Mollusca, Gastropoda)

Helicoidea is a diverse group of land snails with a global distribution. While much is known regarding the relationships of helicoid taxa, comparatively little is known about the evolution of the mitochondrial genome in the superfamily. We sequenced two complete mitochondrial genomes from Praticolellamexicana Perez, 2011 representing the first such data from the helicoid family Polygyridae, and used them in an evolutionary analysis of mitogenomic gene order. We found the mitochondrial genome of Praticolellamexicana to be 14,008 bp in size, possessing the typical 37 metazoan genes. Multiple alternate stop codons are used, as are incomplete stop codons. Mitogenome size and nucleotide content is consistent with other helicoid species. Our analysis of gene order suggested that Helicoidea has undergone four mitochondrial rearrangements in the past. Two rearrangements were limited to tRNA genes only, and two involved protein coding genes.

Analysis of mitochondrial genetic diversity of Ustilago maydis in Mexico

The current understanding of the genetic diversity of the phytopathogenic fungus Ustilago maydis is limited. To determine the genetic diversity and structure of U. maydis, 48 fungal isolates were analyzed using mitochondrial simple sequence repeats (SSRs). Tumours (corn smut or 'huitlacoche') were collected from different Mexican states with diverse environmental conditions. Using bioinformatic tools, five microsatellites were identified within intergenic regions of the U. maydis mitochondrial genome. SSRMUM4 was the most polymorphic marker. The most common repeats were hexanucleotides. A total of 12 allelic variants were identified, with a mean of 2.4 alleles per locus. An estimate of the genetic diversity using analysis of molecular variance (AMOVA) revealed that the highest variance component is within states (84%), with moderate genetic differentiation between states (16%) (F_ST= 0.158). A dendrogram generated using the unweighted paired-grouping method with arithmetic averages (UPGMA) and the Bayesian analysis of population structure grouped the U. maydis isolates into two subgroups (K = 2) based on their shared SSRs.

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.

The diversity of mtDNA rns introns among strains of Ophiostoma piliferum, Ophiostoma pluriannulatum and related species

Background

Based on previous studies, it was suspected that the mitochondrial rns gene within the Ophiostomatales is rich in introns. This study focused on a collection of strains representing Ophiostoma piliferum, Ophiostoma pluriannulatum and related species that cause blue-stain; these fungi colonize the sapwood of trees and impart a dark stain. This reduces the value of the lumber. The goal was to examine the mtDNA rns intron landscape for these important blue stain fungi in order to facilitate future annotation of mitochondrial genomes (mtDNA) and to potentially identify mtDNA introns that can encode homing endonucleases which may have applications in biotechnology.

Results

Comparative sequence analysis identified five intron insertion sites among the ophiostomatoid fungi examined. Positions mS379 and mS952 harbor group II introns, the mS379 intron encodes a reverse transcriptase, and the mS952 intron encodes a potential homing endonuclease. Positions mS569, mS1224, and mS1247 have group I introns inserted and these encode intact or eroded homing endonuclease open reading frames (ORF). Phylogenetic analysis of the intron ORFs showed that they can be found in the same insertion site in closely and distantly related species.

Conclusions

Based on the molecular markers examined (rDNA internal transcribed spacers and rns introns), strains representing O. pilifera, O. pluriannulatum and Ophiostoma novae-zelandiae could not be resolved. Phylogenetic studies suggest that introns are gained and lost and that horizontal transfer could explain the presence of related intron in distantly related fungi. With regard to the mS379 group II intron, this study shows that mitochondrial group II introns and their reverse transcriptases may also follow the life cycle previously proposed for group I introns and their homing endonucleases. This consists of intron invasion, decay of intron ORF, loss of intron, and possible reinvasion.

Electronic supplementary material

The online version of this article (doi:10.1186/s40064-016-3076-6) contains supplementary material, which is available to authorized users.

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.

Intron Derived Size Polymorphism in the Mitochondrial Genomes of Closely Related Chrysoporthe Species

In this study, the complete mitochondrial (mt) genomes of Chrysoporthe austroafricana (190,834 bp), C. cubensis (89,084 bp) and C. deuterocubensis (124,412 bp) were determined. Additionally, the mitochondrial genome of another member of the Cryphonectriaceae, namely Cryphonectria parasitica (158,902 bp), was retrieved and annotated for comparative purposes. These genomes showed high levels of synteny, especially in regions including genes involved in oxidative phosphorylation and electron transfer, unique open reading frames (uORFs), ribosomal RNAs (rRNAs) and transfer RNAs (tRNAs), as well as intron positions. Comparative analyses revealed signatures of duplication events, intron number and length variation, and varying intronic ORFs which highlighted the genetic diversity of mt genomes among the Cryphonectriaceae. These mt genomes showed remarkable size polymorphism. The size polymorphism in the mt genomes of these closely related Chrysoporthe species was attributed to the varying number and length of introns, coding sequences and to a lesser extent, intergenic sequences. Compared to publicly available fungal mt genomes, the C. austroafricana mt genome is the second largest in the Ascomycetes thus far.

The mitochondrial genome of the ethanol-metabolizing, wine cellar mold Zasmidium cellare is the smallest for a filamentous ascomycete

Fungi in the class Dothideomycetes often live in extreme environments or have unusual physiology. One of these, the wine cellar mold Zasmidium cellare, produces thick curtains of mycelia in cellars with high humidity, and its ability to metabolize volatile organic compounds is thought to improve air quality. Whether these abilities have affected its mitochondrial genome is not known. To fill this gap, the circular-mapping mitochondrial genome of Z. cellare was sequenced and, at only 23 743 bp, is the smallest reported for a filamentous fungus. Genes were encoded on both strands with a single change of direction, different from most other fungi but consistent with the Dothideomycetes. Other than its small size, the only unusual feature of the Z. cellare mitochondrial genome was two copies of a 110-bp sequence that were duplicated, inverted and separated by approximately 1 kb. This inverted-repeat sequence confused the assembly program but appears to have no functional significance. The small size of the Z. cellare mitochondrial genome was due to slightly smaller genes, lack of introns and non-essential genes, reduced intergenic spacers and very few ORFs relative to other fungi rather than a loss of essential genes. Whether this reduction facilitates its unusual biology remains unknown.

Overview of genomic and bioinformatic resources for Zymoseptoria tritici

Zymoseptoria tritici (syn. Mycosphaerella graminicola, Septoria tritici) is a haploid fungus belonging to the class Dothideomycetes. It is the causal agent of septoria leaf blotch - one of the world's most significant diseases of wheat. Here we review the genomic and bioinformatic resources that have been generated for Z. tritici. These include the whole-genome reference assembly for isolate IPO323, genome resequencing of alternate isolates, mitochondrial genome sequences, transcriptome sequences and expression data, and annotations of gene structure and function. We also highlight important advances in our fundamental knowledge of genome evolution and its effects on adaptation and pathogenicity in Z. tritici that have been facilitated by these resources.

Complete mitochondrial genome of the Verticillium-wilt causing plant pathogen Verticillium nonalfalfae

Verticillium nonalfalfae is a fungal plant pathogen that causes wilt disease by colonizing the vascular tissues of host plants. The disease induced by hop isolates of V. nonalfalfae manifests in two different forms, ranging from mild symptoms to complete plant dieback, caused by mild and lethal pathotypes, respectively. Pathogenicity variations between the causal strains have been attributed to differences in genomic sequences and perhaps also to differences in their mitochondrial genomes. We used data from our recent Illumina NGS-based project of genome sequencing V. nonalfalfae to study the mitochondrial genomes of its different strains. The aim of the research was to prepare a V. nonalfalfae reference mitochondrial genome and to determine its phylogenetic placement in the fungal kingdom. The resulting 26,139 bp circular DNA molecule contains a full complement of the 14 "standard" fungal mitochondrial protein-coding genes of the electron transport chain and ATP synthase subunits, together with a small rRNA subunit, a large rRNA subunit, which contains ribosomal protein S3 encoded within a type IA-intron and 26 tRNAs. Phylogenetic analysis of this mitochondrial genome placed it in the Verticillium spp. lineage in the Glomerellales group, which is also supported by previous phylogenetic studies based on nuclear markers. The clustering with the closely related Verticillium dahliae mitochondrial genome showed a very conserved synteny and a high sequence similarity. Two distinguishing mitochondrial genome features were also found-a potential long non-coding RNA (orf414) contained only in the Verticillium spp. of the fungal kingdom, and a specific fragment length polymorphism observed only in V. dahliae and V. nubilum of all the Verticillium spp., thus showing potential as a species specific biomarker.

Complete mitochondrial genome of the medicinal fungus Ophiocordyceps sinensis

As part of a genome sequencing project for Ophiocordyceps sinensis, strain 1229, a complete mitochondrial (mt) genome was assembled as a single circular dsDNA of 157,510 bp, one of the largest reported for fungi. Conserved genes including the large and small rRNA subunits, 27 tRNA and 15 protein-coding genes, were identified. In addition, 58 non-conserved open reading frames (ncORFs) in the intergenic and intronic regions were also identified. Transcription analyses using RNA-Seq validated the expression of most conserved genes and ncORFs. Fifty-two introns (groups I and II) were found within conserved genes, accounting for 68.5% of the genome. Thirty-two homing endonucleases (HEs) with motif patterns LAGLIDADG (21) and GIY-YIG (11) were identified in group I introns. The ncORFs found in group II introns mostly encoded reverse transcriptases (RTs). As in other hypocrealean fungi, gene contents and order were found to be conserved in the mt genome of O. sinensis, but the genome size was enlarged by longer intergenic regions and numerous introns. Intergenic and intronic regions were composed of abundant repetitive sequences usually associated with mobile elements. It is likely that intronic ncORFs, which encode RTs and HEs, may have contributed to the enlarged mt genome of O. sinensis.

Evolutionary defined role of the mitochondrial DNA in fertility, disease and ageing

BACKGROUND

The endosymbiosis of an alpha-proteobacterium and a eubacterium a billion years ago paved the way for multicellularity and enabled eukaryotes to flourish. The selective advantage for the host was the acquired ability to generate large amounts of intracellular hydrogen-dependent adenosine triphosphate. The price was increased reactive oxygen species (ROS) inside the eukaryotic cell, causing high mutation rates of the mitochondrial DNA (mtDNA). According to the Muller's ratchet theory, this accumulation of mutations in asexually transmitted mtDNA would ultimately lead to reduced reproductive fitness and eventually extinction. However, mitochondria have persisted over the course of evolution, initially due to a rapid, extreme evolutionary reduction of the mtDNA content. After the phylogenetic divergence of eukaryotes into animals, fungi and plants, differences in evolution of the mtDNA occurred with different adaptations for coping with the mutation burden within these clades. As a result, mitochondrial evolutionary mechanisms have had a profound effect on human adaptation, fertility, healthy reproduction, mtDNA disease manifestation and transmission and ageing. An understanding of these mechanisms might elucidate novel approaches for treatment and prevention of mtDNA disease.

METHODS

The scientific literature was investigated to determine how mtDNA evolved in animals, plants and fungi. Furthermore, the different mechanisms of mtDNA inheritance and of balancing Muller's ratchet in these species were summarized together with the consequences of these mechanisms for human health and reproduction.

RESULTS

Animal, plant and fungal mtDNA have evolved differently. Animals have compact genomes, little recombination, a stable number of genes and a high mtDNA copy number, whereas plants have larger genomes with variable gene counts, a low mtDNA copy number and many recombination events. Fungal mtDNA is somewhere in between. In plants, the mtDNA mutation rate is kept low by effective ROS defence and efficient recombination-mediated mtDNA repair. In animal mtDNA, these mechanisms are not or less well-developed and the detrimental mutagenesis events are controlled by a high mtDNA copy number in combination with a genetic bottleneck and purifying selection during transmission. The mtDNA mutation rates in animals are higher than in plants, which allow mobile animals to adapt more rapidly to various environmental conditions in terms of energy production, whereas static plants do not have this need. Although at the level of the species, these mechanisms have been extremely successful, they can have adverse effects for the individual, resulting, in humans, in severe or unpredictably segregating mtDNA diseases, as well as fertility problems and unhealthy ageing.

CONCLUSIONS

Understanding the forces and processes that underlie mtDNA evolution among different species increases our knowledge on the detrimental consequences that individuals can have from these evolutionary end-points. Alternative outcomes in animals, fungi and plants will lead to a better understanding of the inheritance of mtDNA disorders and mtDNA-related fertility problems. These will allow the development of options to ameliorate, cure and/or prevent mtDNA diseases and mtDNA-related fertility problems.

Characterization and phylogenetic analysis of the mitochondrial genome of Shiraia bambusicola reveals special features in the order of pleosporales

Shiraia bambusicola P. Henn. is a pathogenic fungus of bamboo, and its fruiting bodies are regarded as folk medicine. We determined and analyzed its complete mitochondrial DNA sequence (circular DNA molecule of 39,030 bp, G + C content of 25.19%). It contains the typical genes encoding proteins involved in electron transport and coupled oxidative phosphorylation (nad1-6 and nad4L, cob and cox1-3), one ATP synthase subunit (atp6), 4 hypothetical proteins, and two genes for large and small rRNAs (rnl and rns). There is a set of 32 tRNA genes comprising all 20 amino acids, and these genes are evenly distributed on the two strands. Phylogenetic analyses based on concatenated mitochondrial proteins indicated that S. bambusicola clustered with members of the order Pleosporales, which is in agreement with previous results. The gene arrangements of Dothideomycetes species contained three regions of gene orders partitioned in their mitochondrial genomes, including block 1 (nad6-atp6), block 2 (nad1-cox3) and block 3 (genes around rns). S. bambusicola displayed unique special features that differed from the other Pleosporales species, especially in the coding regions around rns (trnR-trnY). Moreover, a comparison of gene orders in mitochondrial genomes from Pezizomycotina revealed that although all encoded regions are located on the same strand in most Pezizomycotina mtDNAs, genes from Dothideomycetes species had different orientations, as well as diverse positions and colocalization of genes (such as cox3, cox1-cox2 and nad2-nad3); these distinctions were regarded as class-specific features. Interestingly, two incomplete copies of the atp6 gene were found on different strands of the mitogenomic DNA, a finding that has not been observed in the other analyzed fungal species. In our study, mitochondrial genomes from Dothideomycetes species were comprehensively analyzed for the first time, including many species that have not appeared in previous reports.

Analysis of the complete mitochondrial genome of Pochonia chlamydosporia suggests a close relationship to the invertebrate-pathogenic fungi in Hypocreales

Background

The fungus Pochonia chlamydosporia parasitizes nematode eggs and has become one of the most promising biological control agents (BCAs) for plant-parasitic nematodes, which are major agricultural pests that cause tremendous economic losses worldwide. The complete mitochondrial (mt) genome is expected to open new avenues for understanding the phylogenetic relationships and evolution of the invertebrate-pathogenic fungi in Hypocreales.

Results

The complete mitogenome sequence of P. chlamydosporia is 25,615 bp in size, containing the 14 typical protein-coding genes, two ribosomal RNA genes, an intronic ORF coding for a putative ribosomal protein (rps3) and a set of 23 transfer RNA genes (trn) which recognize codons for all amino acids. Sequence similarity studies and syntenic gene analyses show that 87.02% and 58.72% of P. chlamydosporia mitogenome sequences match 90.50% of Metarhizium anisopliae sequences and 61.33% of Lecanicillium muscarium sequences with 92.38% and 86.04% identities, respectively. A phylogenetic tree inferred from 14 mt proteins in Pezizomycotina fungi supports that P. chlamydosporia is most closely related to the entomopathogenic fungus M. anisopliae. The invertebrate-pathogenic fungi in Hypocreales cluster together and clearly separate from a cluster comprising plant-pathogenic fungi (Fusarium spp.) and Hypocrea jecorina. A comparison of mitogenome sizes shows that the length of the intergenic regions or the intronic regions is the major size contributor in most of mitogenomes in Sordariomycetes. Evolutionary analysis shows that rps3 is under positive selection, leading to the display of unique evolutionary characteristics in Hypocreales. Moreover, the variability of trn distribution has a clear impact on gene order in mitogenomes. Gene rearrangement analysis shows that operation of transposition drives the rearrangement events in Pezizomycotina, and most events involve in trn position changes, but no rearrangement was found in Clavicipitaceae.

Conclusions

We present the complete annotated mitogenome sequence of P. chlamydosporia. Based on evolutionary and phylogenetic analyses, we have determined the relationships between the invertebrate-pathogenic fungi in Hypocreales. The invertebrate-pathogenic fungi in Hypocreales referred to in this paper form a monophyletic group sharing a most recent common ancestor. Our rps3 and trn gene order results also establish a foundation for further exploration of the evolutionary trajectory of the fungi in Hypocreales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0341-8) contains supplementary material, which is available to authorized users.