Skip to main content

Mitochondrial genome evolution in species belonging to the Phialocephala fortinii s.l. - Acephala applanata species complex

Abstract

Background

Mitochondrial (mt) markers are successfully applied in evolutionary biology and systematics because mt genomes often evolve faster than the nuclear genomes. In addition, they allow robust phylogenetic analysis based on conserved proteins of the oxidative phosphorylation system. In the present study we sequenced and annotated the complete mt genome of P. subalpina, a member of the Phialocephala fortinii s.l. – Acephala applanata species complex (PAC). PAC belongs to the Helotiales, which is one of the most diverse groups of ascomycetes including more than 2,000 species. The gene order was compared to deduce the mt genome evolution in the Pezizomycotina. Genetic variation in coding and intergenic regions of the mtDNA was studied for PAC to assess the usefulness of mt DNA for species diagnosis.

Results

The mt genome of P. subalpina is 43,742 bp long and codes for 14 mt genes associated with the oxidative phosphorylation. In addition, a GIY-YIG endonuclease, the ribosomal protein S3 (Rps3) and a putative N-acetyl-transferase were recognized. A complete set of tRNA genes as well as the large and small rRNA genes but no introns were found. All protein-coding genes were confirmed by EST sequences. The gene order in P. subalpina deviated from the gene order in Sclerotinia sclerotiorum, the only other helotialean species with a fully sequenced and annotated mt genome. Gene order analysis within Pezizomycotina suggests that the evolution of gene orders is mostly driven by transpositions. Furthermore, sequence diversity in coding and non-coding mtDNA regions in seven additional PAC species was pronounced and allowed for unequivocal species diagnosis in PAC.

Conclusions

The combination of non-interrupted ORFs and EST sequences resulted in a high quality annotation of the mt genome of P. subalpina, which can be used as a reference for the annotation of other mt genomes in the Helotiales. In addition, our analyses show that mtDNA loci will be the marker of choice for future analysis of PAC communities.

Background

Phialocephala fortinii s.l., an anamorphic ascomycete [1, 2] belonging to the Helotiales [2, 3], has been identified as an ubiquitous colonizer of woody plant roots colonizing up to 90% of the roots of woody plant species [4]. The geographical distribution of PAC species ranges from polar regions [5], over temperate regions [6], to subtropical regions [7]. Phialocephala fortinii s.l. was shown to be composed of at least 21 reproductively isolated lineages, eight of which were formally described [4, 8]. In addition, a closely related but sterile species also colonizing roots of woody plant species endophytically was described as Acephala applanata[9]. These species are also known as Phialocephala fortinii s.l. – Acephala applanata species complex (PAC). PAC species form communities of up to 10 sympatrically occurring species and communities were shown to remain stable for several years [1012]. No distance-decay relationship was observed among PAC communities collected across the Northern hemisphere and therefore the Baas-Becking hypothesis that “everything is everywhere” could not be rejected for this assemblage of species [4]. Although sequencing of the internal transcribed spacer (ITS) regions of the rDNA is often regarded as a ‘gold standard’ in species diagnosis in fungi [13], the resolution of ITS sequences was not sufficient to differentiate species in this complex [10]. Instead, several classes of molecular markers were developed for PAC species assignment including PCR fingerprinting, single-copy restriction fragment length polymorphisms (RFLP), multilocus sequence typing, and microsatellites. Although each of these molecular markers supported the delineation of multiple species in this complex, with concordant cryptic species defined by all markers [14, 15], the application of these markers is laborious. In addition, introgression or incomplete lineage sorting further complicates species diagnosis in this species complex [14]. A single, short sequence that allows unequivocally diagnosing PAC species, and if possible other closely related taxa, is still missing.

Mitochondrial (mt) markers were successfully applied in evolutionary biology and systematics [1618] because mt genomes often evolve faster than the nuclear genomes especially in intergenic regions [19, 20]. The mt genomes of filamentous ascomycetes range between 24 to over 100 kb and have a circular topology. They usually contain 14 mt genes encoding proteins of the oxidative phosphorylation system (OXPHOS), the large (rnl) and small (rns) ribosomal subunits, and a varying set of tRNAs genes [21, 22]. In addition, a varying number group I and group II introns specific to fungal mt genomes often including GIY-YIG or LAGLIDAGD endonucleases were reported [2326] and several open reading frames with unknown functions were described [20, 23, 27]. In Pezizomycotina, completely sequenced and annotated mt genomes are available for members of the Eurotiomycetes and Sordariomycetes but only partial non-annotated genomes or draft annotations are available for helotialean species. The Helotiales is one of the most diverse fungal orders and is comprised of more than 350 genera including over 2,000 described species including many important plant pathogens [28].

In the course of a genome sequencing project of P. subalpina a draft sequence of the mt genome became available. We used this draft sequence and re-sequenced the complete mt genome of P. subalpina. Specifically we aimed to: (i) present a completely sequenced and annotated mt genome for the helotialean branch of the fungal tree of life, (ii) compare gene orders of P. subalpina with those found in other filamentous ascomycetes, (iii) compare the evolution of mt and nuclear genomes in PAC species and, (v) test mt loci as a tool for species diagnosis in PAC.

Results

Genome content and genome organization

The circular mt genome of P. subalpina strain UAMH 11012 is 43,742 bp long [GenBank: JN031566] with an AT-content of 72% and contains 21 open reading frames, 14 of which code for OXPHOS proteins (atp6, atp8, atp9, cox1-3, cob, nad1-6, nad4L; Figure 1). All 21 protein-coding ORFs are transcribed in the same direction and start with the canonical translation initiation codon ATG except cox1 (TTG) and cox3 (GTG). The preferred stop codon was TAA with the exception of nad3 and nad5 (both TAG). Beside the 14 OXPHOS genes, the ribosomal protein S3 (Rps3) and a GIY-YIG endonuclease were recognized. In addition, a putative N-acetyl-transferase with an Acetyltransf_1 domain at the C-terminus was predicted (Pfam Nr: PF00583, amino acid positions 361–526, E-value: 1.8E-23). blastp searches returned several hits for Acetylase_1 containing proteins in other ascomycetes. However, all of these proteins were located in the nuclear genome. Similarly, tblastn searches of the putative mt N-acetyl-transferase against a draft genome sequence of P. subalpina recognized two additional N-acetyl-transferases for P. subalpina but both nuclear proteins were substantially smaller (198 and 202 aa) than the predicted mt protein (579 aa). The functions of the remaining four ORFs are unknown because neither significant blastp hits nor conserved domains in interproscan searches were found. In addition, the large and small ribosomal RNA subunits (rns and rnl) were present in the mt genome. A number of 5,903 Roche/454 GS FLX (454) EST sequencing reads (total number of aligned bases: 1.89 Mbp) mapped to the mt genome. Reads almost exclusively mapped to coding regions of the mt genome (proteins and rRNA genes) and all ORFs except ORF_03 were partially or fully covered by ESTs. However, sequence coverage differed considerably among genes (see Additional file 1).

Figure 1
figure 1

Map of the mt genome of Phialocephala subalpina. Map displaying the circular mt genome of P. subalpina strain UAMH 11012. All open reading frames, tRNA genes and the large ribosomal RNA are transcribed clockwise.

In total, 27 tRNA genes coding for all amino acids were predicted in the mt genome. Both software tools applied resulted in the same tRNA predictions except tRNA-Cys at position 5,363, which was only predicted by rnaweasel and tRNA-Asn at position 20,122 that was exclusively predicted by trnascan-se. A putative Ochre suppressor tRNA gene (antisense codon: UUA) was predicted by trnascan-se. Protein-coding genes covered 48.1% of the mt genome, while 4.7% and 11.8% correspond to tRNA genes and rRNA genes, respectively (total coding regions: 64.6%).

Duplication of a region including the atp9 open reading frame in PAC

A duplication of an mtDNA region including the atp9 open reading frame was observed in P. subalpina (Figure 1). Beside the atp9 ORF (225 bp), the duplication included 124 bp upstream of the start codon and 37 bp downstream of the stop codon of the atp9 ORF (total length of duplication: 386 bp). The duplication was interrupted by a 488 bp insertion (Figures 1 and 2). The putative ORF_01 (123 aa) shared the first 14 bp of the duplicated atp9 ORF but the remaining 358 bp of ORF_01 including the stop codon lie in the 488 bp insertion interrupting the duplication (Figure 2, bottom).

Figure 2
figure 2

Genome organisation of the region between cox1 and nad1 for eight PAC species. Genome organisation for 8 PAC species between cox1 and nad1 covering the duplication, which includes the atp9 ORF. Boxes represent ORFs and tRNA genes. Thick lines represent the position of the duplicated region and indels (black: duplicated regions upstream and downstream the atp9 open reading frame; maroon: region of the atp9 ORF; white: deletions; grey: inserts). The duplication including the atp9 open reading frame was present in all species at the nucleotide level, although some species included inserts (e.g. Acephala applanata) or deletions (e.g. Phialocephala turicensis). ORF_01 of Phialocephala fortinii s.s. includes a premature stop codon indicated by an asterisk. An example of an alignment for the duplicated region in P. subalpina is given below the overview. The intron is removed and its length was indicated in brackets. Start and stop positions of atp9 and ORF_01 are indicated by arrows. Numbers beside the alignment represent the mt genome coordinates for P. subalpina.

Analysis of the cox1 nad1 region for seven additional PAC strains (Table 1) showed that the synteny between ORF_02 and nad1 was conserved among all PAC species and sequence identity was > 90%. In contrast, it was not possible to derive significant alignments based on nucleotide data for all 8 species between cox1 and ORF_02. Two groups were recognized that produced significant alignments over the entire cox1 nad1 sequences (Figure 2). Whereas the first group is comprised of P. turicensis P. letzii P europaea and A. applanata, the second group included P. helvetica P. uotolensis P. subalpina and P. fortinii s.s. The differences between cox1 and ORF_02 for the two groups are also reflected by the annotation of different gene features between these two groups. Species in group 1 possess an additional tRNA-Cys between cox1 and ORF_02 except for A. applanata and no insert was observed in the duplicated region including the atp9 ORF. The tRNA-Cys was not recognized by rnaweasel in A. applanata due to the accumulation of mutations in this particular region compared to the other species of group 1. In contrast, members of group 2 possess a tRNA-Ser between ORF_01 and ORF_02, which was, however, not present in P. helvetica due to a 118 bp deletion (Figure 2). In addition, P. helvetica deviated from other group 2 species by the absence of the 488 bp insertion leading to the annotation of ORF_01. Genome organisation between cox1 and nad1 reflected the known phylogenetic relatedness of PAC species with the exception of P. uotolensis. Strain 5_134_3 of P. uotolensis, which is closely related to P. turicensis based on sequences of several nuclear loci, single-copy RFLP analysis and microsatellite data [14, 15], showed a similar genome organization as P. fortinii s.s. and P. subalpina.

Table 1 PAC strains included to study the duplication of the atp9 region and the conservation of putative ORFs

The presence of a duplication of the region including the atp9 ORF could be confirmed for all eight genomes at the nucleotide level (Figure 2). Three species had two non-interrupted ORFs for the atp9 gene (P. letzii, P. europaea, P. helvetica) with identical protein sequences for P. letzii and P. europaea. The two P. helvetica atp9 protein sequences deviate at two of the 74 aa positions between the two atp9 ORFs. All other species have either an insertion leading to ORF_01 (P. uotolensis, P. fortinii s.s., P. subalpina) or a deletion (P. turicensis) in one of the atp9 ORFs (Figure 2). Similarly, A. applanata possessed an insertion of 234 bp in one of the atp9 ORFs (Figure 2). The insert shows neither direct nor inverted repeats and did not show any similarities with known nucleotide or protein sequences using blastn and blastx searches.

Conservation of putative ORFs in PAC species

In order to collect further evidence that putative ORFs code for proteins, the conservation of these ORFs was investigated in seven additional PAC species. Unfortunately, repeated attempts to amplify the region of ORF_04 failed. In contrast, we could assess the presence and conservation of the other 3 putative ORFs observed in P. subalpina in the additional PAC species. ORF_02 was observed in all eight PAC species. Despite SNPs and indels, none of the indels resulted in frameshift mutations and no premature stop codons were present (see Additional file 2). In contrast, ORF_03 was either missing completely, contained stop codons, or the start codon was missing in five of the eight species. Similarly, ORF_01 was only present in three of the eight species and P. fortinii s.s. had a premature stop codon (Figure 2).

Evolution of gene orders in the Helotiales and Pezizomycotina

Anchored genome alignments for the three helotialean species P. subalpina, S. sclerotiorum and B. cinerea show pronounced genome rearrangements between P. subalpina and the closely related species S. sclerotiorum and B. cinerea (Figure 3). In contrast, complete synteny in gene order were evident between S. sclerotiorum and B. cinerea. Phylogenetic analysis based on 12 OXPHOS protein sequences (atp6, cox1-3, cob, nad1-6, nad4L) allowed for the identification of the four classes of filamentous ascomycetes (Figure 4, Table 2). Although the ML tree shows that the Dothideomycetes and Leotiomycetes share a most recent common ancestor with the Eurotiomycetes, the BS support value (53%) and bayesian PP value (0.55) for this placement are low. Therefore, alternative hypotheses such as Dothideomycetes and Leotiomycetes sharing a most recent common ancestor with the Sordariomycetes cannot be ruled out.

Figure 3
figure 3

Genome rearrangements observed among three helotialean species. Genome rearrangements detected by Mauve genome alignments for the three helotialean species Phialocephala subalpina Sclerotinia sclerotiorum and Botrytis cinerea. Locally collinear blocks identified by mauve are given in different colours and were compared with the annotated gene features in P. subalpina. The Rps3 forms a free-standing ORF in P. subalpina and was not placed in a group-I intron located in the U11 domain of rnl[26]. The rns gene was used as anchor to linearize the three genomes.

Figure 4
figure 4

Evolution of gene orders in Pezizomycotina. A: Single ML tree (log likelihood: -50,671.86) constructed from 12 OXPHOS proteins for selected ascomycetes with complete mt genomes (see Table 2). The tree was used to map gene order rearrangements of 14 OXPHOS proteins, Rps3, and the rRNAs. Thick branches indicate bootstrap support values ≥90% in ML analysis and posterior probabilities ≥0.95 in BI. Candida albicans served as outgroup. Group numbering given on the right indicates groups of species with identical gene orders for the gene set analyzed. Ancestral gene orders are indicated next to the nodes as G followed by the group index if the gene order is identical to the gene order of a given group or as A for the gene order reconstructed by treerex (see B), except for the root node. Boxes on the braches indicate mapped rearrangements as analysed by crex and treerex. Numbers in boxes refer to the rearrangement given in C. B: Reconstructed ancestral gene orders for the Eurotiomycetes, Sordariomycetes and the Leotiomycetes. C: Evolutionary scenarios to deduce one gene order of the other for each change indicated in Figure 4A are given (trp, transpositions; tdrl, tandem-duplication-random-loss). Only the mt regions included in the operation are indicated in the Figure.

Table 2 Ascomycetous species included in gene order analysis

Based on this phylogeny, we have reconstructed putative ancestral gene orders and rearrangements leading to the contemporary gene orders of 14 OXPHOS proteins (atp6, atp8, atp9, cox1-3, cob, nad1-6, nad4L), Rps3, and rRNAs by using crex and treerex analysis (Figures 4A, B; see Additional file 3 for full input gene orders). Identical ancestral gene arrangements were observed in treerex analysis irrespective whether the Leotiomycetes and Dothideomycetes share the most recent common ancestor with the Sordariomycetes or the Eurotiomycetes (data not shown). The evolution of gene orders among species belonging to the Pezizomycotina was mainly characterized by transpositions (Figures 4A, C). In contrast, inversions and inverse-transpositions have not been found in the analysed gene orders, which excluded the two Dothideomycete species. However, the Dothideomycetes have genes that are in opposite orientation with respect to the remaining Pezizomycotina species. Moreover, both species were separated from all other Pezizomycotina species and by themselves by very long evolutionary scenarios rendering the reconstruction unreliable (Figure 4A). In our analysis, only one tandem-duplication-random-loss (trdl) event was suggested by crex analysis involving cob and cox1 to derive the gene order of S. sclerotiorum (Rearrangement 11 in Figure 4). blastx searches of cob and cox1 protein sequences against the mtDNA of S. sclerotiorum confirmed that a duplication-based rearrangement mechanism was involved because a partial cob fragment in front of cox1 was found as expected for the predicted operation (30 aa, E-value: 3.00E-09). Together with the duplication of atp9 for PAC species (see above) this indicates that duplication-based rearrangement are involved in mitogenome evolution of Pezizomycotina. The number of necessary events required to convert one gene order to the other was not always related to the phylogenetic distance between species. For example, gene orders of Podospora anserina and S. sclerotiorum could only be derived from gene orders of their closest relatives by evolutionary scenarios that include multiple transpositions and in the case of S. sclerotiorum a tdrl event.

Based on the used species and gene sampling our analysis indicates that i) the group of Beauveria bassiana (Group 1 in Figure 4) represent the ancestral gene order for the Sordariomycetes, i.e., no gene rearrangements are necessary to reach the node in the phylogenetic tree representing the Sordariomycetes. Similarly, Aspergillus niger (Group 8 in Figure 4) represents the ancestral gene order for the Eurotiomycetes. Furthermore, comparisons among the possible intermediate gene orders of the Sordariomycetes, Eurotiomycetes, and the reconstructed ancestral gene order of the Leotiomycetes indicates the intermediate position of the later (Additional file 4).

Evolution of mtDNA in PAC and searching regions for species diagnosis

A total of 7,350 bp (16.8% of the mt genome) from four regions of the mt genome were amplified in four additional PAC species to screen for polymorphisms (Table 3). Fragments included coding (cox1 rnl) and intergenic regions (Rps3- putative N-acetyl-transferase, atp9-nad1-nad4). High amounts of DNA polymorphisms were observed within coding and intergenic regions among PAC species (Additional file 5). However, DNA polymorphisms were not always evenly distributed among species within the sequenced DNA fragments. Whereas in three cases all five species contributed to a similar extent to the observed polymorphic regions (e.g. rnl), P. europaea was almost exclusively responsible for a peak in diversity observed in cox1 (data not shown). In addition, a 520 bp indel in the non-coding region spanning Rps3 and the putative N-acetyl-transferase showed a very high nucleotide diversity among the three species with this insert (P. europaea P. fortinii s.s., A. applanata) (Additional File 5). Based on these results three regions were selected and tested on a broader selection of strains. Parts of rnl (~730 bp), a fragment located between Rps3 and the putative N-acetyl-transferase (~850 bp), and a fragment including parts of nad1 and nad4 (~1,650 bp) were sequenced from 32 strains belonging to eight PAC species (Table 3). Best-fitting mutation models selected from jModeltest were F81 + I, TVM + G, TrN + G, and TPM3uf + I + G for the three individual loci and the concatenated dataset, respectively. Variation in the three fragments was considerable among and within species. Often single strains of a species were responsible for this high variation found within species (Figure 5). For example, strains 7_45_5, or T1_K92_131 were separated from the other strains of the respective species for all three loci resulting in 19 and 21 substitutions respectively, not considering indels. P. turicensis was the only species that was monomorphic for all three loci. ML and BI analysis using the concatenated dataset placed strains of each species in well-supported clades (Figure 7). In addition, the tree topology of each individual mtDNA locus and the concatenated mtDNA dataset did not deviate significantly from the topology of the concatenated sequences of four genomic regions sequenced previously [14] (Table 4). No evidence for recombination blocks among the three sequenced loci were found, although incompatible sites were regularly observed (data not shown). Incompatible sites were located within as well as between the three sequenced loci.

Figure 5
figure 5

Phylogenetic analysis for PAC species based on three mt loci. Single ML tree for eight PAC species based on DNA sequence data of three concatenated mt loci. Posterior probabilities of BI (above branches) and bootstrap values of ML analysis (below branches) are indicated. Acephala applanata was chosen as outgroup.

Table 3 PAC strains included to study the suitability of mtDNA markers for species diagnosis
Table 4 Concordance in evolution among mtDNA and nucDNA for PAC

Discussion

In the present study we sequenced and annotated the mt genome of a widely distributed group of root-inhabiting fungi belonging to the Helotiales (Leotiomycetes). Gene order analysis showed that the evolution of mt genomes in the Pezizomycotina is mainly driven by transpositions. Moreover, we show the usefulness of mtDNA loci for species diagnosis in the Phialocephala fortinii s.l. – Acephala applanata species complex (PAC).

Gene content and mt organisation in P. subalpina

The mt genome of filamentous ascomycetes normally code for 14 proteins of the oxidative phosphorylation system [22], including the sequences for the large (rnl) and small (rns) ribosomal subunits. In addition, a varying number of additional proteins with homologies to known proteins, i.e., the ribosomal protein S3 (Rps3) [23, 24], and ORFs with unknown functions may be found [20, 23, 27]. The mt genome of P. subalpina follows this rule. All 14 OXPHOS proteins, the two rRNAs and a complete set of tRNA genes were present. Nevertheless, some unique features were observed. (i) A putative N-acetyl-transferase was predicted in the mt genome. To the best of our knowledge, this is the first report of a putative N-acetyl-transferase in the mt genome of ascomycetes. Acetyl-transferases modify proteins in eukaryotes, both co- and post-translationally by transferring acetyl groups from acetyl-coenzyme A to either the a-amino group of amino-terminal residues or to the e-amino group of lysine residues at various positions [30]. In addition, N-acetyl-transferases were shown to modify several ribosomal proteins in Escherichia coli[31, 32]. Although it seems likely that the putative mt N-acetyl-transferase is involved in protein modifications, the exact function is unknown because N-acetyl-transferases can act on different groups of substrates [30]. (ii) The ribosomal protein S3 (Rps3) formed a free-standing ORF in P. subalpina (Figure 1). To the best of our knowledge this is the second report of a free-standing Rps3 gene in the Pezizomycotina. In all completely sequenced mt genomes of filamentous ascomycetes possessing Rps3, this gene is placed in a group-I intron located in the U11 domain of the rnl except for Phaeospaeria nodorum[23, 24, 26, 33]. Whereas the Rps3 of P. nodorum is large (771 aa) and includes parts of cox1[34], Rps3 found in P. subalpina is similar in length as Rps3 found in S. sclerotiorum and other ascomycetes. In addition, blastp searches revealed no segments of other mt proteins. Rps3 followed cox2 in both species and sequencing of additional mt genomes in the Dothideomycetes and Leotiomycetes will show whether cox2 Rps3 synteny may be another common position of Rps3 in mt genomes. (iii) An Ochre suppressor tRNA gene was predicted. Suppressor tRNAs were found in the genomes of many species and allow the read-through of stop codons [3537]. They provide a regulatory mechanism of gene expression allowing the production of several proteins from a single gene and were shown to be especially important for RNA viruses [38]. Moreover, suppressor tRNAs are hypothesized to play a role for the transcription of cryptic mitochondrial gene on the antisense strand [39, 40]. (iv) A duplication of a genomic region including the atp9 open reading frame interrupted by a 488 bp insert was observed in P. subalpina, leading to the annotation of ORF_0 1. Re-sequencing of the cox1-ORF_01-ORF_02-atp9-nad1 region for additional seven PAC species confirmed that three species have two intact ORFs coding for atp9. In all other species one of the two atp9 ORFs was interrupted by indels indicating that the duplication of the atp9 region followed by successive modifications predated the separation of PAC species. Interestingly, two types of sequences with no obvious similarity were found between cox1 and the duplicated atp9 regions for PAC species, which was reflected in the annotated gene features and was in accordance with the relatedness of the species. An exception was P. uotolensis strain 5_134_3, which showed a similar sequence and organisation as P. subalpina and P. fortinii s.s. In previous studies we showed for this strain that it might be the result of a hybridization event with P. subalpina. For example the same strain clustered with a P. subalpina strain for the sequence of the translation elongation factor 1-α [14] but three other nuclear loci, microsatellite analysis and single-copy RFLP analysis placed this strain with other strains of P. uotolensis[8, 14, 15].

Transcription of ORFs in the mt genome

Identification of mtDNA genes is no guarantee that they are also active, as was e.g. shown for a silent copy of the cox2 gene in the soybean mitochondrion [19, 41]. Mapping of 454 EST sequence reads showed that they were mostly restricted to ORFs and rDNA genes. All ORFs except ORF_03 were confirmed by ESTs. However, pronounced differences in the sequence coverage were observed. Several factors could explain the differences observed in the transcription profile of coding regions. Either the 454 sequencing protocol introduced non-uniform coverage in the normalized cDNA library [42] or the differences in the coverage reflect some differences in the expression of mt genes. Indeed, transcription of mt genes in the ascomycete Saccharomyces cerevisiae was shown to be far from uniform [43]. Future RT-qPCR assays for P. subalpina will show whether the observed differences in cDNA coverage reflect differences in expression profiles for mt proteins.

Do putative ORFs code for proteins?

Putative ORFs are regularly annotated in fungal genomes [20, 23, 27, 44]. However, whether they really code for proteins is rarely assessed. In the present study we used an indirect approach to gain additional evidence whether putative ORFs code for proteins. We assume that putative ORFs that code for a protein should (i) be transcribed, i.e., ESTs should match the respective regions and (ii) nucleotide sequences of the ORFs should be conserved among closely related species, in particular premature stop codons and frameshifts should not be detected. ORF_02 is an example of a transcribed and conserved ORF that most likely codes for a protein although blastp searches recognized no similar proteins in the NCBI database. In contrast, ORF_03 was either completely deleted or included premature stop codons in four species and the insert responsible for ORF_01 was only present in three species and included a frame-shift mutation in one of these species. Therefore, we hypothesize that ORF_01 and ORF_03 do not code for proteins because they are not universally conserved despite their (partial) transcription (see Additional file 1). However, definitive evidence that conserved and transcribed ORFs code for proteins is only given when the corresponding proteins are isolated [19].

Evolution of gene orders in filamentous ascomycetes

Phylogenetic analysis based on protein data of 12 OXPHOS genes resulted in a robust phylogeny with the exception of the Dothideomycetes and Leotiomycetes which clustered together, but their position was not resolved and it remained unclear whether they share the most common ancestor with the Sordariomycetes or the Eurotiomycetes. Indeed, several deep phylogenetic analysis based on nuclear genes and/or rDNA loci placed the Leotiomycetes in a cluster with the Sordariomycetes [43, 44]. To attribute for the uncertainty of the placement of Leotiomycetes, we analyzed the evolution of gene orders using both alternative hypotheses.

Gene orders were compared previously for some fungal orders [16, 45, 46] but no automated and formally well-defined approach was used to reconstruct the rearrangement history. crex is based on the notion of common intervals, which reflect genes that appear consecutively in the input gene orders. In addition, we used treerex resulting in a detailed reconstruction of the rearrangement history and our analysis revealed that the evolution of gene orders in filamentous ascomycetes is mainly driven by transposition events and that gene orders are (mostly) conserved for close relatives. The conservation of gene orders in close relatives was also shown in previous studies. For example, the transposition leading to the gene order of Verticillium dahliae (Rearrangement 3 in Figure 4) was confirmed in seven additional Verticillium species [16]. Similarly, 3 species closely related to A. niger (Aspergillus tubingensis Penicillium chrysogenum and Penicillium digitatum) were recently shown to have the same gene order as A. niger for the set of genes used in our analysis [47].

Gene order analyses were shown to be useful to confirm phylogenetic analysis in several organismal groups [4549] and in the present study we show, that the extent of gene rearrangements reflects the phylogenetic position of a species in most cases. For example, the ancestral gene order of the Leotiomycetes was most closely related to the two ancestral gene orders of the Sordariomycetes and Eurotiomycetes based on the parsimonious rearrangement scenarios reconstructed by crex reflecting an intermediate position of the Leotiomycetes (see Additional File 4). This fits well with the results of the phylogenetic analysis showing that the Leotiomycetes could share the most recent common ancestor either with the Sordariomycetes or the Eurotiomycetes. However, two exceptions were observed. First, P. anserina deviated from N. crassa by several transpositions and the loss of the atp9 gene. Similarly, the two Dothideomycete species deviated from themselves and all other Pezizomycotina by complex rearrangement scenarios. It was possible to reconstruct the possible ancestral gene order for the Sordariomycetes and the Eurotiomycetes although the presented results must be regarded as preliminary because of the limited species sampling in particular in the Leotiomycetes and the exclusion of the tRNAs due to restrictions of the software tools used. Since tRNAs are important for understanding the evolution of gene orders [49, 50] the presented results need to be verified with a data set including the tRNAs.

Evolution of mtDNA in PAC species and suitability of mt markers for species diagnosis

The amount of variation found within PAC species was high, despite the fact that strains from most of the species were derived from a single study site (Table 3). Intra-species variability in the mt genome among species belonging to ascomycetes can vary greatly. Whereas mtDNA was highly diverse in many collections of filamentous ascomycetes [20, 51, 52], low levels of diversity were observed in M. graminicola for a world-wide collection of strains [27]. The high amount of variation in PAC species could be the result of genetic exchange among mt genomes. Although mt genomes were characterized by the absence of any genetic recombination in some studies [27] an increasing number of reports show recombination in mtDNA of plants and fungi [5254]. Therefore, we searched for inconsistencies in our dataset to find evidence for recombination events. Although inconsistencies were observed among parsimony informative sites for PAC species, the inconsistencies occurred within loci as well as between loci and no recombination blocks were evident, indicating that rather parallel mutations than recombination characterize the evolution of the mt genome in PAC species.

Mt loci were often used as diagnostic markers in systematics. They received particular attention during barcoding campaigns [55] and were also successfully applied for fungal species diagnosis [18, 20, 56]. In the present study it was possible to analyze the resolution of mtDNA markers at a fine taxonomic level using a well-defined species complex. Several marker types were used to characterize cryptic species in PAC such as single-copy RFLPs, microsatellite markers and sequence markers [14, 15]. In the present study we show that two of the three examined mtDNA regions allow to distinguish all eight PAC species. Only one strain of P. uotolensis was misplaced using the sequence information of the partial rnl gene. The amount of variation found in three mt loci was high compared to nuclear loci sequenced previously [14]. For example, P. letzii and P. europaea have identical sequences for several nuclear coding and non-coding loci and were difficult to diagnose based on single-copy RFLP markers [14] but formed well supported monophyletic groups for all mt loci. So far, mating type idiomorphs were the only nuclear regions with a similar ability to distinguish PAC species [57]. However, a disadvantage of the mating type for species diagnosis is their heterothallic organisation in PAC species, which does not allow amplifying a single fragment for all strains. In this respect, mt markers developed in the present study offer advantages and mt loci can be regarded as useful diagnostic markers for PAC.

Conclusions

The availability of the complete annotated mt genome and the knowledge about the intra- and inter-species diversity in PAC provides the basis for the development of new markers to study the community ecology, population biology and evolution in this species complex. In addition, it provides a reference for the annotation of other mt genomes in the Helotiales.

Methods

Sequencing the complete mt genome of Phialocephala subalpina

In the course of a genome sequencing project of P. subalpina strain UAMH 11012, an initial Roche/454 GS FLX (454) shotgun run was performed at the Functional Genomics Centre Zurich (FGCZ, Uni/ETH Zurich) and from that a draft of the circular mt genome of P. subalpina strain UAMH 11012 became available. The draft sequence was subdivided into 12 fragments (see Additional file 6) and amplified from strain UAMH 11012 using long-range PCR in 20 μl volumes (Expand Long Range dNTPack kit, Roche, Rotkreuz, Switzerland). PCR conditions were optimized for 4 kb fragments with high AT contents by lowering the temperature during the elongation step from 68°C to 62°C and an initial elongation time of 4 min. PCR fragments were purified using the Wizard Plus SV kit (Promega, Wallisellen, Switzerland) and sequenced at Microsynth (Balgach, Switzerland).

Annotation of the mt genome

Open reading frames (ORFs) in the mtDNA sequence of P. subalpina were searched with NCBI orf finder[58]. blastp[59] and interproscan[60] were used for homology-based function prediction of proteins. In addition, homologous regions of putative ORFs were sequenced for seven PAC species and sequences were analyzed for indels and SNPs resulting in frameshift mutations and/or pre-mature stop codons. The tRNA genes were predicted by trnascan-se v1.21[61] and rnaweasel[62] searches using default settings for mt genomes. The ribosomal RNAs were determined by comparison with sequences from other fungi using blastn. Group I and group II introns were predicted using rnaweasel[62]. Duplications in the mt genome were searched using the software dotter[63].

EST evidence for ORFs

Expression of ORFs was tested by searching for ESTs derived from a normalized EST library of strain UAMH 11012 sequenced on the 454 platform at FGCZ (Rémy Bruggmann, unpublished). Single 454 cDNA reads belonging to the mt genome were filtered using the software genomethreader[64] and the sequence coverage for each base in the mt genome was calculated (= number of reads covering a specific base in the genome).

Analyzing the duplication of the atp9 region in PAC species

A region between cox1 and nad1 included a duplication of the atp9 region in P. subalpina (see results). In order to test the presence of the duplication in other PAC species, the region between cox1 and nad1 was sequenced from additional 7 PAC species (Table 1). The entire locus was amplified using primers mtPF_F07_F27774 (TAGAGGTAATCAAACCAATG) and PF_cox1_F (AGCCCACCAAAACCTCATGC) using long-range PCR and sequenced as described above.

Evolution of gene orders in Pezizomycotina

In a first step, mt genome alignments for the helotialean species B. cinerea S. sclerotiorum and P. subalpina based on the publically available sequences were performed using the mauve 2.3.1 software [65] and locally collinear blocks identified by mauve were compared with the annotated gene features in P. subalpina. In a second step, the phylogenetic relationship among species belonging to the Pezizomycotina was analyzed using a concatenated protein dataset of 12 OXPHOS proteins (atp6 cox1-3 cob nad1-6 nad4L). Proteins coding for atp8 (48 aa) and atp9 (74 aa) are missing in one or more of the species included in the analysis (atp8: P. nodorum; atp9: P. anserina and P. nodorum) and we omitted these two genes in the phylogenetic analysis. The protein data set was aligned using mafft[66] and the mafft-homologues option was applied. The resulting alignment was trimmed in gblocks 0.91b[67] using default settings. Candida albicans served as outgroup in the phylogenetic analysis (see below).

In a third step, gene order of all 14 OXPHOS proteins, Rps3, and rRNAs (rns and rnl) was studied for 22 species in Pezizomycotina (Table 2) using crex[68] and treerex[69] analysis. A limitation that applies to most gene order analysis approaches, including crex and treerex is that they can only be applied to gene orders with an equal gene content with each gene appearing exactly once. Therefore, we had to exclude tRNA genes because the identification of homologous tRNA genes was hampered in our dataset due to gene losses and/or multiple occurrences of the same tRNA genes (see Additional file 7). Moreover, tRNAs are difficult to predict by the available tools as we have shown for P. subalpina (see Results), rendering them uncertain candidates for gene order analysis. An informal introduction to crex and treerex analysis as well as a detailed description of the performed gene order analysis is given in the Additional file 4. The Dothideomycetes were excluded from analysis because both species were separated from all other Pezizomycotina and by themselves by very long evolutionary scenarios including several inversion and tdrl events which makes the reconstruction unreliable. Similarly, the gene order of P. anserina was excluded since TreeREx cannot handle gene orders with unequal gene content (loss of atp9 in P. anserina). The scenario shown in Figure 4 is derived from the separate analysis of N. crassa and P. anserina considering the common genes only.

Searching regions for species diagnosis in PAC

The following strategy was applied to identify regions suitable for species diagnosis in PAC. In a first step, four fragments each between 1,200-2,800 bp long spanning coding and intergenic regions in the mtDNA were chosen (cox1 rnl nad1-nad4 Rps3- putative N-acetyl-transferase) and amplified in four additional PAC species (A. applanata P. fortinii s.s., P. helvetica, and P. europaea) including the phylogenetically most distantly related PAC species known to date (P. europaea and A. applanata) [14]. Coding regions were chosen because they were either shown to be suitable for species diagnosis previously [17] or have the potential to be used to diagnose also closely related species of PAC. In contrast, intergenic regions were chosen because they were shown to be highly polymorphic [20]. Nucleotide diversity for the five species was analyzed in dnasp v5.0[70] and sites[71]. Based on these results, one coding region (rnl) and two intergenic regions (nad1-nad4 Rps3- putative N-acetyl-transferase) were selected. These three loci were then tested on a broader collection of strains to study the evolution of mtDNA within PAC species and their suitability for species diagnosis. The dataset includes 32 strains belonging to 8 PAC species (Table 3). For primers used to amplify and sequence fragments see Additional file 8. Fragments were amplified in a 15 μl reaction volume using approximately 5 ng of template DNA. After an initial denaturation step for 2 min at 94°C, 31 cycles were performed each consisting of a denaturation step at 94°C for 30 s, an annealing step at 50°C for 30 s and an extension step at 60°C for 90 s [72] followed by a final extension step for 6 min at 72°C. Fragments were directly purified using an ExoSap protocol [15] and sequencing was conducted at the Genetic Diversity Centre (GDC, ETH Zurich). The nucleotide sequences of each locus were aligned with Clustal W separately and the resulting data matrices were then concatenated into one combined data matrix and subjected to phylogenetic analysis.

Phylogentic analysis

Phylogenetic trees were inferred using Bayesian inference (BI) and maximum likelihood (ML) methods for the protein and nucleotide dataset. ML analysis for the protein dataset was performed in treefinder[73] using the substitution model selected by protest[74]. ML analysis for nucleotide datasets were performed in paup[75] using the substitution model selected by jmodeltest v1.0[76]. Branch supports were provided by 1,000 bootstrap replicates for both datasets. BI trees were calculated with mrbayes 3.1[77] for protein and nucleotide datasets by running two analyses each consisting of two simultaneous runs with four heated chains per run. Each analysis was run for 5 Mio generations and trees were sampled every 100 generations. Post-burn-in was assumed when the average standard deviation of split frequencies was consistently ≤ 0.01. Post-burn-in trees were collected and the parameter and topology summarizations calculated. To ensure that the analyses reached stationarity and converged on the same topology and likelihood scores, the resulting likelihoods, tree topologies and model estimates were compared by eye.

The tree topologies derived from the nucleotide data sets (each single locus and the concatenated alignment) were compared with tree topologies from the concatenated sequences of four genomic regions (b-tubulin elongation factor 1-α pPF_018 and pPF_076) sequenced previously [14] using the approximately unbiased (AU) and the weighted Shimodaira and Hasegawa (WSH) tests [29]. Only one strain (type strain) per species was included to compare tree topologies. Both tests were calculated in treefinder using 10,000 replicates. In addition, a compatibility matrix of nucleotide substitutions was constructed using the software Sites[71] to search for possible recombination blocks in the datasets.

Abbreviations

PAC:

Phialocephala fortinii s.l. – Acephala applanata species complex

OXPHOS:

Proteins of the oxidative phosphorylation system

ORF:

Open reading frame.

References

  1. Grünig CR, Sieber TN, Rogers SO, Holdenrieder O: Genetic variability among strains of Phialocephala fortinii and phylogenetic analysis of the genus Phialocephala based on rDNA ITS sequence comparisons. Can J Bot. 2002, 80 (12): 1239-1249. 10.1139/b02-115.

    Article  Google Scholar 

  2. Wang Z, Binder M, Schoch CL, Johnston PR, Spatafora JW, Hibbett DS: Evolution of helotialean fungi (Leotiomycetes, Pezizomycotina): A nuclear rDNA phylogeny. Mol Phylogenet Evol. 2006, 41 (2): 295-312. 10.1016/j.ympev.2006.05.031.

    Article  CAS  PubMed  Google Scholar 

  3. Grünig CR, Queloz V, Duò A, Sieber TN: Phylogeny of Phaeomollisia piceae gen. sp. nov.: a dark-septate conifer-needle endophyte and its relationships to Phialocephala and Acephala. Mycol Res. 2009, 113 (2): 207-221. 10.1016/j.mycres.2008.10.005.

    Article  PubMed  Google Scholar 

  4. Queloz V, Sieber TN: Holdenrieder, McDonald BA, Grünig CR: No biogeographical pattern for a root-associated fungal species complex. Global Ecol Biogeogr. 2011, 20 (1): 160-169. 10.1111/j.1466-8238.2010.00589.x.

    Article  Google Scholar 

  5. Piercey MM, Graham SW, Currah RS: Patterns of genetic variation in Phialocephala fortinii across a broad latitudinal transect in Canada. Mycol Res. 2004, 108 (8): 955-964. 10.1017/S0953756204000528.

    Article  CAS  PubMed  Google Scholar 

  6. Grünig CR, Queloz V, Sieber TN, Holdenrieder O: Dark septate endophytes (DSE) of the Phialocephala fortinii s.l. - Acephala applanata species complex in tree roots—classification, population biology and ecology. Botany. 2008, 86 (12): 1355-1369. 10.1139/B08-108.

    Article  Google Scholar 

  7. Zhang C, Yin L, Dai S: Diversity of root-associated fungal endophytes in Rhododendron fortunei in subtropical forests of China. Mycorrhiza. 2009, 19 (6): 417-423. 10.1007/s00572-009-0246-1.

    Article  PubMed  Google Scholar 

  8. Grünig CR, Duò A, Sieber TN, Holdenrieder O: Assignment of species rank to six reproductively isolated cryptic species of the Phialocephala fortinii s.l.-Acephala applanata species complex. Mycologia. 2008, 100 (1): 47-67. 10.3852/mycologia.100.1.47.

    Article  PubMed  Google Scholar 

  9. Grünig CR, Sieber TN: Molecular and phenotypic description of the widespread root symbiont Acephala applanata gen. et sp. nov., formerly known as dark septate endophyte Type 1. Mycologia. 2005, 97 (3): 628-640. 10.3852/mycologia.97.3.628.

    Article  PubMed  Google Scholar 

  10. Grünig CR, McDonald BA, Sieber TN, Rogers SO, Holdenrieder O: Evidence for subdivision of the root-endophyte Phialocephala fortinii into cryptic species and recombination within species. Fung Genet Biol. 2004, 41 (7): 676-687. 10.1016/j.fgb.2004.03.004.

    Article  Google Scholar 

  11. Grünig CR, Duò A, Sieber TN: Population genetic analysis of Phialocephala fortinii s.l. and Acephala applanata in two undisturbed forests in Switzerland and evidence for new cryptic species. Fung Genet Biol. 2006, 43 (6): 410-421. 10.1016/j.fgb.2006.01.007.

    Article  Google Scholar 

  12. Queloz V, Grünig CR, Sieber TN: Holdenrieder O: Monitoring the spatial and temporal dynamics of a community of the tree-root endophyte Phialocephala fortinii s.l. New Phytol. 2005, 168 (3): 651-660. 10.1111/j.1469-8137.2005.01529.x.

    Article  PubMed  Google Scholar 

  13. Peay KG, Kennedy PG, Bruns TD: Fungal Community Ecology: A Hybrid Beast with a Molecular Master. Bioscience. 2008, 58 (9): 799-810. 10.1641/B580907.

    Article  Google Scholar 

  14. Grünig CR, Brunner PC, Duò A, Sieber TN: Suitability of methods for species recognition in the Phialocephala fortinii - Acephala applanata species complex using DNA analysis. Fung Genet Biol. 2007, 44 (8): 773-788. 10.1016/j.fgb.2006.12.008.

    Article  Google Scholar 

  15. Queloz V, Duo A, Sieber TN, Grünig CR: Microsatellite size homoplasies and null alleles do not affect species diagnosis and population genetic analysis in a fungal species complex. Mol Ecol Resour. 2010, 10: 348-367. 10.1111/j.1755-0998.2009.02757.x.

    Article  CAS  PubMed  Google Scholar 

  16. Pantou MP, Kouvelis VN, Typas MA: The complete mitochondrial genome of the vascular wilt fungus Verticillium dahliae: a novel gene order for Verticillium and a diagnostic tool for species identification. Curr Genet. 2006, 50 (2): 125-136. 10.1007/s00294-006-0079-9.

    Article  CAS  PubMed  Google Scholar 

  17. Seifert KA, Samson RA, Dewaard JR, Houbraken J, Lévesque CA, Moncalvo J-M, Louis-Seize G, Hebert PDN: Prospects for fungus identification using CO1 DNA barcodes, with Penicillium as a test case. Proc Natl Acad Sci USA. 2007, 104 (10): 3901-3906. 10.1073/pnas.0611691104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Kouvelis VN, Sialakouma A, Typas MA: Mitochondrial gene sequences alone or combined with ITS region sequences provide firm molecular criteria for the classification of Lecanicillium species. Mycol Res. 2008, 112: 829-844. 10.1016/j.mycres.2008.01.016.

    Article  CAS  PubMed  Google Scholar 

  19. Burger G, Gray M, Lang B: Mitochondrial genomes: anything goes. Trends Genet. 2003, 19 (12): 709-716. 10.1016/j.tig.2003.10.012.

    Article  CAS  PubMed  Google Scholar 

  20. Ghikas DV, Kouvelis VN, Typas MA: Phylogenetic and biogeographic implications inferred by mitochondrial intergenic region analyses and ITS1–5.8S-ITS2 of the entomopathogenic fungi Beauveria bassiana and B. brongniartii. BMC Microbiol. 2010, 10: 174-10.1186/1471-2180-10-174.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Gray M, Burger G, Lang B: Mitochondrial evolution. Science. 1999, 283 (5407): 1476-1481. 10.1126/science.283.5407.1476.

    Article  CAS  PubMed  Google Scholar 

  22. Lavin JL, Oguiza JA, Ramirez L, Pisabarro AG: Comparative genomics of the oxidative phosphorylation system in fungi. Fungal Genet Biol. 2008, 45 (9): 1248-1256. 10.1016/j.fgb.2008.06.005.

    Article  CAS  PubMed  Google Scholar 

  23. Wu Y, Yang J, Yang F, Liu T, Leng W, Chu Y, Jin Q: Recent dermatophyte divergence revealed by comparative and phylogenetic analysis of mitochondrial genomes. BMC Genomics. 2009, 10: 238-10.1186/1471-2164-10-238.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Sethuraman J, Majer A, Friedrich NC, Edgell DR, Hausner G: Genes within Genes: Multiple LAGLIDADG Homing Endonucleases Target the Ribosomal Protein S3 Gene Encoded within an rnl Group I Intron of Ophiostoma and Related Taxa. Mol Biol Evol. 2009, 26 (10): 2299-2315. 10.1093/molbev/msp145.

    Article  CAS  PubMed  Google Scholar 

  25. Monteiro-Vitorello CB, Hausner G, Searles DB, Gibb EA, Fulbright DW, Bertrand H: The Cryphonectria parasitica mitochondrial rns gene: Plasmid-like elements, introns and homing endonucleases. Fungal Genet Biol. 2009, 46 (11): 837-848. 10.1016/j.fgb.2009.07.005.

    Article  CAS  PubMed  Google Scholar 

  26. Gibb E, Hausner G: Optional mitochondrial introns and evidence for a homing-endonuclease gene in the mtDNA rnl gene in Ophiostoma ulmi s. lat. Mycol Res. 2005, 109: 1112-1126. 10.1017/S095375620500376X.

    Article  CAS  PubMed  Google Scholar 

  27. Torriani SFF, Goodwin SB, Kema GHJ, Pangilinan JL, McDonald BA: Intraspecific comparison and annotation of two complete mitochondrial genome sequences from the plant pathogenic fungus Mycosphaerella graminicola. Fungal Genet Biol. 2008, 45 (5): 628-637. 10.1016/j.fgb.2007.12.005.

    Article  CAS  PubMed  Google Scholar 

  28. Kirk PM, Cannon PF, David JC, Stalpers JA: (Eds): Dictionary of the fungi. 2001, CAB International, Oxon, UK, 9

    Google Scholar 

  29. Shimodaira H: An approximately unbiased test of phylogenetic tree selection. Syst Biol. 2002, 51 (3): 492-508. 10.1080/10635150290069913.

    Article  PubMed  Google Scholar 

  30. Polevoda B, Sherman F: The diversity of acetylated proteins. Genome Biol. 2002, 3 (5): reviews0006-reviews0006.6.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Yoshikawa A, Isono S, Sheback A, Isono K: Cloning and nucleotide sequencing of the genes rimI and rimJ which encode enzymes acetylating ribosomal proteins S18 and S5 of Escherichia coli K12. Mol Gen Genet. 1987, 209 (3): 481-488. 10.1007/BF00331153.

    Article  CAS  PubMed  Google Scholar 

  32. Isono K, Isono S: Ribosomal protein modification in Escherichia coli. II. Studies of a mutant lacking the N-terminal acetylation of protein S18. Mol Gen Genet. 1980, 177 (4): 645-651.

    Article  CAS  PubMed  Google Scholar 

  33. Hane JK, Lowe RGT, Solomon PS, Tan K-C, Schoch CL, Spatafora JW, Crous PW, Kodira C, Birren BW, Galagan JE: Dothideomycete-plant interactions illuminated by genome sequencing and EST analysis of the wheat pathogen Stagonospora nodorum. Plant Cell. 2007, 19 (11): 3347-3368. 10.1105/tpc.107.052829.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Sethuraman J, Majer A, Iranpour M, Hausner G: Molecular Evolution of the mtDNA Encoded rps3 Gene Among Filamentous Ascomycetes Fungi with an Emphasis on the Ophiostomatoid Fungi. J Mol Evol. 2009, 69 (4): 372-385. 10.1007/s00239-009-9291-9.

    Article  CAS  PubMed  Google Scholar 

  35. Debuchy R, Brygoo Y: Cloning of opal suppressor tRNA genes of a filamentous fungus reveals two tRNASerUGA genes with unexpected structural differences. EMBO J. 1985, 4 (13A): 3553-3556.

    PubMed Central  CAS  PubMed  Google Scholar 

  36. Grimm M, Nass A, Schüll C, Beier H: Nucleotide sequences and functional characterization of two tobacco UAG suppressor tRNA(Gln) isoacceptors and their genes. Plant Mol Biol. 1998, 38 (5): 689-697. 10.1023/A:1006068303683.

    Article  CAS  PubMed  Google Scholar 

  37. Baum M, Beier H: Wheat cytoplasmic arginine tRNA isoacceptor with a U*CG anticodon is an efficient UGA suppressor in vitro. Nucleic Acids Res. 1998, 26 (6): 1390-1395. 10.1093/nar/26.6.1390.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Beier H, Grimm M: Misreading of termination codons in eukaryotes by natural nonsense suppressor tRNAs. Nucleic Acids Res. 2001, 29 (23): 4767-4782. 10.1093/nar/29.23.4767.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Faure E, Delaye L, Tribolo S, Levasseur A, Seligmann H, Barthélémy R-M: Probable presence of an ubiquitous cryptic mitochondrial gene on the antisense strand of the cytochrome oxidase I gene. Biol Direct. 2011, 6: 56-10.1186/1745-6150-6-56.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Seligmann H: Two genetic codes, one genome: frameshifted primate mitochondrial genes code for additional proteins in presence of antisense antitermination tRNAs. Biosystems. 2011, 105 (3): 271-285. 10.1016/j.biosystems.2011.05.010.

    Article  CAS  PubMed  Google Scholar 

  41. Covello PS, Gray MW: Silent mitochondrial and active nuclear genes for subunit 2 of cytochrome c oxidase (cox2) in soybean: evidence for RNA-mediated gene transfer. EMBO J. 1992, 11 (11): 3815-3820.

    PubMed Central  CAS  PubMed  Google Scholar 

  42. Harismendy O, Ng PC, Strausberg RL, Wang XY, Stockwell TB, Beeson KY, Schork NJ, Murray SS, Topol EJ, Levy S: Evaluation of next generation sequencing platforms for population targeted sequencing studies. Genome Biol. 2009, 10 (3): R32-10.1186/gb-2009-10-3-r32.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Costanzo MC, Fox TD: Control of mitochondrial gene expression in Saccharomyces cerevisiae. Annu Rev Genet. 1990, 24: 91-113. 10.1146/annurev.ge.24.120190.000515.

    Article  CAS  PubMed  Google Scholar 

  44. Cummings DJ, McNally KL, Domenico JM, Matsuura ET: The complete DNA sequence of the mitochondrial genome of Podospora anserina. Curr Genet. 1990, 17 (5): 375-402. 10.1007/BF00334517.

    Article  CAS  PubMed  Google Scholar 

  45. Kunisawa T: Inference of the phylogenetic position of the phylum Deferribacteres from gene order comparison. Antonie Van Leeuwenhoek. 2011, 99 (2): 417-422. 10.1007/s10482-010-9492-7.

    Article  PubMed  Google Scholar 

  46. Sankoff D, Leduc G, Antoine N, Paquin B, Lang B, Cedergren R: Gene order comparisons for phylogenetic inference - Evolution of the mitochondrial genome. Proc Natl Acad Sci USA. 1992, 89 (14): 6575-6579. 10.1073/pnas.89.14.6575.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Sun X, Li H, Yu D: Complete mitochondrial genome sequence of the phytopathogenic fungus Penicillium digitatum and comparative analysis of closely related species. FEMS Microbiol Lett. 2011, 323 (1): 29-34. 10.1111/j.1574-6968.2011.02358.x.

    Article  CAS  PubMed  Google Scholar 

  48. Kunisawa T: The phylogenetic placement of the non-phototrophic, Gram-positive thermophile 'Thermobaculum terrenum' and branching orders within the phylum 'Chloroflexi' inferred from gene order comparisons. Int J Systematic Evol Microbiol. 2011, 61: 1944-1953. 10.1099/ijs.0.026088-0.

    Article  CAS  Google Scholar 

  49. Saccone C, Gissi C, Reyes A, Larizza A, Sbisa E, Pesole G: Mitochondrial DNA in metazoa: degree of freedom in a frozen event. Gene. 2002, 286 (1): 3-12. 10.1016/S0378-1119(01)00807-1.

    Article  CAS  PubMed  Google Scholar 

  50. Jühling F, Pütz J, Bernt M, Donath A, Middendorf M: Florentz C. 2011, Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements. Nucleic Acids Res, Stadler PF

    Google Scholar 

  51. Sommerhalder RJ, McDonald BA, Zhan J: Concordant evolution of mitochondrial and nuclear genomes in the wheat pathogen Phaeosphaeria nodorum. Fungal Genet Biol. 2007, 44 (8): 764-772. 10.1016/j.fgb.2007.01.003.

    Article  CAS  PubMed  Google Scholar 

  52. Anderson J, Wickens C, Khan M, Cowen L, Federspiel N, Jones T, Kohn L: Infrequent genetic exchange and recombination in the mitochondrial genome of Candida albicans. J Bacteriol. 2001, 183 (3): 865-872. 10.1128/JB.183.3.865-872.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. van Diepeningen AD, Goedbloed DJ, Slakhorst SM, Koopmanschap AB, Maas MFPM, Hoekstra RF, Debets AJM: Mitochondrial recombination increases with age in Podospora anserina. Mech Ageing Dev. 2010, 131 (5): 315-322. 10.1016/j.mad.2010.03.001.

    Article  CAS  PubMed  Google Scholar 

  54. Barr CM, Neiman M, Taylor DR: Inheritance and recombination of mitochondrial genomes in plants, fungi and animals. New Phytol. 2005, 168 (1): 39-50. 10.1111/j.1469-8137.2005.01492.x.

    Article  CAS  PubMed  Google Scholar 

  55. Seifert KA: Progress towards DNA barcoding of fungi. Mol Ecol Resour. 2009, 9: 83-89.

    Article  CAS  PubMed  Google Scholar 

  56. Ghikas DV, Kouvelis VN, Typas MA: The complete mitochondrial genome of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae: gene order and trn gene clusters reveal a common evolutionary course for all Sordariomycetes, while intergenic regions show variation. Arch Microbiol. 2006, 185 (5): 393-401. 10.1007/s00203-006-0104-x.

    Article  CAS  PubMed  Google Scholar 

  57. Zaffarano PL, Duò A, Grünig CR: Characterization of the mating type (MAT) locus in the Phialocephala fortinii s.l. -Acephala applanata species complex. Fungal Genet Biol. 2010, 47 (9): 761-772. 10.1016/j.fgb.2010.06.001.

    Article  CAS  PubMed  Google Scholar 

  58. NCBI ORF finder:.http://www.ncbi.nlm.nih.gov/gorf/gorf.html,

  59. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. EMBOSS tools:.http://www.ebi.ac.uk/Tools/pfa/iprscan/,

  61. Lowe T, Eddy S: tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25 (5): 955-964.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Lang BF, Laforest M-J, Burger G: Mitochondrial introns: a critical view. Trends Genet. 2007, 23 (3): 119-125. 10.1016/j.tig.2007.01.006.

    Article  CAS  PubMed  Google Scholar 

  63. Sonnhammer EL, Durbin R: A dot-matrix program with dynamic threshold control suited for genomic DNA and protein sequence analysis. Gene. 1995, 167 (1–2): GC1-GC10.

    CAS  PubMed  Google Scholar 

  64. Gremme G, Brendel V, Sparks M, Kurtz S: Engineering a software tool for gene structure prediction in higher organisms. Inform Software Tech. 2005, 47 (15): 965-978. 10.1016/j.infsof.2005.09.005.

    Article  Google Scholar 

  65. Darling ACE, Mau B, Blattner FR, Perna NT: Mauve: Multiple Alignment of Conserved Genomic Sequence With Rearrangements. Genome Res. 2004, 14 (7): 1394-1403. 10.1101/gr.2289704.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. MAFFT:.http://mafft.cbrc.jp/alignment/server/index.html,

  67. Gblocks 0.91b:.http://www.phylogeny.fr/version2_cgi/one_task.cgi?task_type=gblocks,

  68. Bernt M, Merkle D, Ramsch K, Fritzsch G, Perseke M, Bernhard D, Schlegel M, Stadler PF, Middendorf M: CREx: inferring genomic rearrangements based on common intervals. Bioinformatics. 2007, 23 (21): 2957-2958. 10.1093/bioinformatics/btm468.

    Article  CAS  PubMed  Google Scholar 

  69. Bernt M, Merkle D, Middendorf M: An Algorithm for Inferring Mitochondrial Genome Rearrangements in a Phylogenetic Tree. In Comparative Genomics International Workshop, RECOMB-CG: Paris, Lecture Notes in Computer Sciences (LNCS). Berlin: Springer. 2008, 2008: 143-157.

    Google Scholar 

  70. Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19 (18): 2496-2497. 10.1093/bioinformatics/btg359.

    Article  CAS  PubMed  Google Scholar 

  71. Hey J, Wakeley J: A coalescent estimator of the population recombination rate. Genetics. 1997, 145: 833-846.

    PubMed Central  CAS  PubMed  Google Scholar 

  72. Su X, Wu Y, Sifri C, Wellems T: Reduced extension temperatures required for PCR amplification of extremely A + T-rich DNA. Nucleic Acids Res. 1996, 24 (8): 1574-1575. 10.1093/nar/24.8.1574.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  73. Jobb G: TREEFINDER version of November 2010.http://www.treefinder.de/,

  74. Darriba D, Taboada GL, Doallo R, Posada D: ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011, 27 (8): 1164-1165. 10.1093/bioinformatics/btr088.

    Article  CAS  PubMed  Google Scholar 

  75. Swofford DL: PAUP: pylogenetic analysis using parsimony, version 3.1.1. 1993, Illinois Natural History Survey, Champaign, IL

    Google Scholar 

  76. Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25 (7): 1253-1256. 10.1093/molbev/msn083.

    Article  CAS  PubMed  Google Scholar 

  77. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank S. Torriani for helpful discussions on a previous draft of this manuscript. The Functional Genomics Centre Zurich (FGCZ, Uni/ETH Zurich) and the Genetic Diversity Centre (GDC, ETH Zurich) provided facilities for collecting the sequence data. This study was partially funded by Vontobel Stiftung, Zürich to CRG.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Christoph R Grünig.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AD and CRG designed the research project. AD and CRG performed all experiments and collected sequence data. RB, SZ, MB and CRG implemented analytical tools and performed analysis. CRG, RB, MB and SZ wrote the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

12864_2011_4129_MOESM1_ESM.eps

Additional file 1: Titel: Transcriptomics landscape in the mt genome of Phialocephala subalpina. Description: Sequence coverage for 5,903 454 cDNA reads (total number of aligned bases: 1.89 Mbp) that mapped to the mt genome of Phialocephala subalpina. 454 reads almost exclusively mapped to coding regions of the mt genome (proteins and rDNAs) but depths of coverage differed considerably among genes. (EPS 686 KB)

12864_2011_4129_MOESM2_ESM.doc

Additional file 2: Titel: Protein aligenemts for putative ORFs. Description: Protein alignments of two putative ORFs (ORF_02, ORF_03) for eight PAC species. (DOC 36 KB)

12864_2011_4129_MOESM3_ESM.txt

Additional file 3: Titel: Gene order data used for TreeREx and CREx analysis. Descriptions: Gene order data for 13 unique gene orders observed in species belonging to the Pezizomycotina and used for TreeREx and CREx analysis. (TXT 2 KB)

12864_2011_4129_MOESM4_ESM.pdf

Additional file 4: Titel: CREx and TreeREx analysis. Description: In the first two sections an informal introduction to crex and treerex is given. For a detailed introduction we refer to Bernt et al. (2007, 2008). Sections 3 to 6 include a detailed description of the gene order analyses performed completing the results included in the paper. (PDF 145 KB)

12864_2011_4129_MOESM5_ESM.eps

Additional file 5: Titel: Polymorphism screening in PAC. Description: Amount and distribution of polymorphisms found at four mt loci for five species of the Phialocephala fortinii s.l. – Acephala applanata species complex (PAC). The four loci included approx. 7,800 bp representing 17.8% of the mt genome. The position of protein coding regions and rDNAs are indicated. Thick bars represent the fragment chosen for testing mt loci for species diagnosis in PAC (8 species, 32 strains). (EPS 984 KB)

12864_2011_4129_MOESM6_ESM.xls

Additional file 6: Titel: List of primers used to sequence the mt genome of Phialocephala subalpina. Description: List of primers used to amplify and sequence the complete mt genome of Phialocephala subalpina. (XLS 26 KB)

12864_2011_4129_MOESM7_ESM.xlsx

Additional file 7: Titel: Presence/absence of tRNA genes in species belonging to the Pezizomycotina included in crex and treerex analysis. (XLSX 57 KB)

12864_2011_4129_MOESM8_ESM.xls

Additional file 8: Titel: Primers used for phylogenetic analyis in the Phialocephala fortinii s.l. – Acephala applanata species complex (PAC). Description: Primers used to amplify three mtDNA loci in eight species belonging to PAC. (XLS 18 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Duò, A., Bruggmann, R., Zoller, S. et al. Mitochondrial genome evolution in species belonging to the Phialocephala fortinii s.l. - Acephala applanata species complex. BMC Genomics 13, 166 (2012). https://doi.org/10.1186/1471-2164-13-166

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-13-166

Keywords