Mitochondrial genome evolution in species belonging to the Phialocephala fortinii s.l. - Acephala applanata species complex
© Duò et al.; licensee BioMed Central Ltd. 2012
Received: 19 May 2011
Accepted: 4 May 2012
Published: 4 May 2012
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.
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.
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.
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 . The geographical distribution of PAC species ranges from polar regions , over temperate regions , to subtropical regions . 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. 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 [10–12]. 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 . Although sequencing of the internal transcribed spacer (ITS) regions of the rDNA is often regarded as a ‘gold standard’ in species diagnosis in fungi , the resolution of ITS sequences was not sufficient to differentiate species in this complex . 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 . 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 [16–18] 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 [23–26] 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 .
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.
Genome content and genome organization
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
PAC strains included to study the duplication of the atp9 region and the conservation of putative ORFs
Culture collection number1
GenBank accession numbers
cox1 - nad1 3
Phialocephala fortinii s.s.
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
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
PAC strains included to study the suitability of mtDNA markers for species diagnosis
Culture collection number2
GenBank accession numbers
Rps3 - N-acetyl- transferase
atp9 - nad4
Phialocephala fortinii s.s.
Concordance in evolution among mtDNA and nucDNA for PAC
Rps3- putative N-
4 nucDNA loci1
4 nucDNA loci
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 , 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 . 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 . (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, 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 [35–37]. 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 . 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-α  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  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 . 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 .
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 . 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 .
Gene order analyses were shown to be useful to confirm phylogenetic analysis in several organismal groups [45–49] 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 . 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  an increasing number of reports show recombination in mtDNA of plants and fungi [52–54]. 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  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 . 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  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 . 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.
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.
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. blastp and interproscan 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 and rnaweasel 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. Duplications in the mt genome were searched using the software dotter.
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 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  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 and the mafft-homologues option was applied. The resulting alignment was trimmed in gblocks 0.91b 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 and treerex 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) . Coding regions were chosen because they were either shown to be suitable for species diagnosis previously  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 . Nucleotide diversity for the five species was analyzed in dnasp v5.0 and sites. 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  followed by a final extension step for 6 min at 72°C. Fragments were directly purified using an ExoSap protocol  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.
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 using the substitution model selected by protest. ML analysis for nucleotide datasets were performed in paup using the substitution model selected by jmodeltest v1.0. Branch supports were provided by 1,000 bootstrap replicates for both datasets. BI trees were calculated with mrbayes 3.1 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  using the approximately unbiased (AU) and the weighted Shimodaira and Hasegawa (WSH) tests . 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 to search for possible recombination blocks in the datasets.
Phialocephala fortinii s.l. – Acephala applanata species complex
Proteins of the oxidative phosphorylation system
Open reading frame.
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.
- 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.View Article
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View Article
- 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.View ArticlePubMed
- 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.View Article
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View Article
- 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.View Article
- 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.View ArticlePubMed
- 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.View Article
- 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.View Article
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- Burger G, Gray M, Lang B: Mitochondrial genomes: anything goes. Trends Genet. 2003, 19 (12): 709-716. 10.1016/j.tig.2003.10.012.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- Gray M, Burger G, Lang B: Mitochondrial evolution. Science. 1999, 283 (5407): 1476-1481. 10.1126/science.283.5407.1476.View ArticlePubMed
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- Kirk PM, Cannon PF, David JC, Stalpers JA: (Eds): Dictionary of the fungi. 2001, CAB International, Oxon, UK, 9
- Shimodaira H: An approximately unbiased test of phylogenetic tree selection. Syst Biol. 2002, 51 (3): 492-508. 10.1080/10635150290069913.View ArticlePubMed
- Polevoda B, Sherman F: The diversity of acetylated proteins. Genome Biol. 2002, 3 (5): reviews0006-reviews0006.6.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- 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 CentralPubMed
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- 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 CentralPubMed
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- 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.View Article
- 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.View ArticlePubMed
- 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
- 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.View ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- Seifert KA: Progress towards DNA barcoding of fungi. Mol Ecol Resour. 2009, 9: 83-89.View ArticlePubMed
- 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.View ArticlePubMed
- 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.View ArticlePubMed
- NCBI ORF finder:.http://www.ncbi.nlm.nih.gov/gorf/gorf.html,
- 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.PubMed CentralView ArticlePubMed
- EMBOSS tools:.http://www.ebi.ac.uk/Tools/pfa/iprscan/,
- 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.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- 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.PubMed
- 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.View Article
- 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.PubMed CentralView ArticlePubMed
- Gblocks 0.91b:.http://www.phylogeny.fr/version2_cgi/one_task.cgi?task_type=gblocks,
- 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.View ArticlePubMed
- 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.
- 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.View ArticlePubMed
- Hey J, Wakeley J: A coalescent estimator of the population recombination rate. Genetics. 1997, 145: 833-846.PubMed CentralPubMed
- 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.PubMed CentralView ArticlePubMed
- Jobb G: TREEFINDER version of November 2010.http://www.treefinder.de/,
- 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.View ArticlePubMed
- Swofford DL: PAUP: pylogenetic analysis using parsimony, version 3.1.1. 1993, Illinois Natural History Survey, Champaign, IL
- Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25 (7): 1253-1256. 10.1093/molbev/msn083.View ArticlePubMed
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.