Mycologia
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

DOI: 10.3852/mycologia.97.4.812
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kim, K. M.
Right arrow Articles by Jung, H. S.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Kim, K. M.
Right arrow Articles by Jung, H. S.
Agricola
Right arrow Articles by Kim, K. M.
Right arrow Articles by Jung, H. S.
Mycologia, 97(4), 2005, pp. 812-822.
© 2005 by The Mycological Society of America

Evaluation of the monophyly of Fomitopsis using parsimony and MCMC methods


Kyung Mo Kim
Yuh-Gang Yoon
Hack Sung Jung 1

     Department of Biological Sciences, College of Natural Sciences, Seoul National University, San 56-1 Shillim-9-dong, Kwanak-gu, Seoul 151-747, Korea

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 

To evaluate the monophyly of Fomitopsis and elucidate phylogenetic relationships of its members, partial nuclear large subunit (partial 28S) ribosomal RNA genes were sequenced from 10 species of Fomitopsis and 15 related species. Phylogenetic analyses indicated that Fomitopsis was phylogenetically heterogeneous and its members were divided into three subgroups. The constrained tree excluding F. palustris (the type species of Pilatoporus) from Fomitopsis core group was rejected, thus rejecting the taxonomic concept to segregate Pilatoporus from Fomitopsis. The monophyly of taxa belonging to F. rosea complex was rejected, thus rejecting the complex definition based on morphological similarities. The exclusion of Piptoporus betulinus (the type species of Piptoporus) from Fomitopsis core group was rejected and Piptoporus proved to be heterogeneous in both best MP and MAP trees. The monophyly of F. officinalis with Fomitopsis core group also was rejected. Fomitopsis officinalis was closely related to Antrodia xantha and formed an independent lineage from Fomitopsis core group at the basal position of brown rotting fungi comprising Antrodia, Daedalea, Fomitopsis, Piptoporus and Postia. The MAP tree topology obtained from MCMC computation of Bayesian inference was similar to the one of the best MP tree based on the parsimony analysis but showed a higher likelihood score in the Kishino-Hasegawa test and reflected better evolutionary patterns for the phylogeny of Fomitopsis.

Key words: Bayesian inference, best MP tree, Kishino-Hasegawa test, MAP tree, 28S rDNA


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Fomitopsis P. Karst. is a well-known genus in the modern classification system of polypores (Pouzar 1966Go) and includes brown rotting fungi with dimitic or trimitic hyphae. Fomitopsis is classified in the Polyporaceae, a large and artificial family composed of numerous genera with poroid hymenophores (Donk 1964Go) and placed in the polyporoid clade of Homobasidiomycetes (Hibbett and Thorn 2001Go). According to Donk’s (1974)Go and Gilbertson’s and Ryvarden’s (1986Go–1987Go) generic descriptions, members of Fomitopsis have characters as follows: the whitish to bright-colored, di- to trimitic context of the cap; the generative hyphae with clamps; mostly ellipsoid to cylindric, hyaline, smooth, IKI- spores; the tendency of the hymenophore to become layered; the presence of a laccate, glabrous crust covering the cap; and the brown rot type in contrast to the white rot activity of Trametes Fr. with same types of spores and hyphal systems as Fomitopsis.

Until recently, taxonomic delimitation of Fomitopsis has been controversial and remained insufficiently resolved (Kotlaba and Pouzar 1998Go). Bondartsev and Singer (1941)Go proposed that the species currently placed in Haploporus Bondartsev & Singer, Heterobasidion Bref. and Perenniporia Murrill with dextrinoid spores or hyphae, and Rigidoporus Murrill with generative hyphae without clamps also be included in the members of Fomitopsis. Corner (1989)Go neglected the type of rot as a generic character and treated Fomitopsis as a synonym of Trametes but did not formally transfer the species of Fomitopsis to Trametes (Ryvarden 1991Go). Fomitopsis rosea (Alb. & Schwein.) P. Karst. is closely related to F. pinicola (Sw.) P. Karst., the type species of Fomitopsis, while F. cajanderi (P. Karst.) Kotl. & Pouzar is closely related to the species of Daedalea Pers. sensu stricto sharing several important morphological characters, according to Nobles’ (1971)Go study. After Nobles’ view, Donk (1974)Go suggested the removal of F. cajanderi from Fomitopsis and the inclusion into Daedalea or any other related group, although F. cajanderi and F. rosea with pinkish to pinkish-vinaceous context were thought to be closely related to each other and hardly distinguishable at species level.

Kotlaba and Pouzar (1990Go, 1998)Go suggested a narrow concept of Fomitopsis and combined F. rosea into a new genus, Rhodofomes Kotl. & Pouzar (type species: Rhodofomes roseus [Alb. & Schwein.] Kotl. & Pouzar), and F. palustris (Berk. & M.A. Curtis) Gilb. & Ryvarden into another new genus, Pilatoporus Kotl. & Pouzar (type species: Pilatoporus palustris (Berk. & M.A. Curtis) Kotl. & Pouzar). Rhodofomes is characterized by thin-walled spores, the presence of clamps on thin-walled generative hyphae, the rosy context and the absence of a resinous crust on the pileal surface. Pilatoporus shares a trimitic hyphal system and hyaline hyphae with Fomitopsis and Rhodofomes but has pseudoskeletal hyphae provided with conspicuous clamp connections (Kotlaba and Pouzar 1990Go). Carranza-Morse and Gilbertson (1986)Go defined F. rosea complex, with F. cajanderi, F. carnea (Blume & Nees) Imazeki, F. cupreorosea (Berk.) J. Carranza & Gilb., F. dochmia (Berk. & Broome) Ryvarden, F. feei (Fr.) Kreisel, F. lilacinogilva (Berk.) J.E. Wright & J.R. Deschamps and F. rosea, with these characters: the rosy context, at least when young; dimitic or trimitic hyphal systems; hyaline, smooth, cylindric to allantoid spores; and no special structures in hymenium with an exception of cystidioles present in F. cajanderi and F. dochmia. However, it has not been confirmed that F. rosea complex based upon morphological similarities is still maintained in molecular phylogenetic analyses.

To elucidate phylogenetic relationships of Fomitopsis and related members, Bayesian inference and maximum parsimony were used in this study. The Bayesian inference is useful in estimating phylogenetic trees from DNA sequence data (Yang and Rannala 1997Go, Larget and Simon 1999Go) and employs stochastic simulation of the Markov chain Monte Carlo (MCMC) method to obtain trees with highest posterior probabilities (Gelman et al 1995Go) calculated according to Bayes’ theorem (Huelsenbeck et al 2001Go). Among trees generated from MCMC computation for a given sequence dataset (Huelsenbeck and Ronquist 2001Go), the final tree can be determined with one of criteria for mean approach and maximum a priori (MAP) (Durbin et al 1998Go). Mean approach produces one consensus tree out of all generated trees. Thus, it reduces the probability of occurrence of local maxima (Durbin et al 1998Go) happening in maximum likelihood computation but increases the probability of occurrence of polytomies in tree topology. On the other hand, MAP criterion (Durbin et al 1998Go) has relatively high probability of occurrence of local maxima but confers higher resolution upon the tree than mean approach. For the resolved topology of ingroups of the phylogenetic tree, MAP approach seemed to be appropriate in phylogenetic studies at intra- and intergeneric levels and was applied here to reconstruct the Bayesian phylogenetic tree of Fomitopsis.

The purpose of this study was to evaluate the monophyly of Fomitopsis and clarify relevant taxonomic questions in Fomitopsis and related polyporoid taxa. To evaluate whether Fomitopsis is monophyletic, members of genera Antrodia P. Karst., Daedalea, Piptoporus P. Karst. and Postia Fr. were incorporated along with F. pinicola (the type species of Fomitopsis), F. officinalis (Vill.) Bondartsev & Singer, F. palustris and F. spraguei (Berk. & M.A. Curtis) Gilb. & Ryvarden, and six species of F. rosea complex TABLE IGo). Those non-Fomitopsis members are all brown rotting fungi that share many morphological characters with the members of Fomitopsis (Gilbertson and Ryvarden 1986Go). Two brown rot species, Laetiporus sulphureus (Bull.) Murrill & Phaeolus schweinitzii (Fr.) Pat., were used as outgroups to root phylogenetic trees.


View this table:
[in this window]
[in a new window]
 
TABLE I. Strains sequenced in this study

 

    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Cultures, DNA extraction, amplification and sequencing.— – Sources of fungal strains used in this study are listed (TABLE IGo). Herbarium specimens from Seoul National University Fungus Collection (SFC) were used and mycelial cultures from Centraalbureau voor Schimmelcultures (CBS) and American Type Culture Collection (ATCC) were maintained on potato-dextrose agar (PDA). Total genomic DNAs were extracted from SFC specimens or from cultured mycelia by the method of Lecellier and Silar (1994)Go with some modifications (Lee et al 2000Go). From extracted genomic DNAs, partial nuclear large subunit (partial 28S) ribosomal DNA region was amplified with primers ITS1 (White et al 1990Go) and LR3 (5'-CCGTGTTTCAAGACGGG-3', Vilgalys and Hester 1990Go) using Quick PCR Premix (GENENMED, Daejeon, Korea). Amplified products were checked by agarose gel electrophoresis and purified with an AccuPrep PCR purification kit (Bioneer Corp., Daejeon, Korea). Purified PCR products were sequenced with ABI3700 automated DNA sequencer (Applied Biosystems, Foster City, California). Primer pairs, LR3 and LR21 (5'-ACTTCAAGCGTTTCCCTTT-3', http://www.biology.duke.edu/fungi/mycolab/) were used for sequencing partial 28S ribosomal DNA.

Phylogenetic analyses.— – Using Clustal X version 1.83 (Thompson et al 1997Go) (gap opening penalty = 10.0; gap extension penalty = 0.05), 32 sequences were aligned with 12 sequences retrieved from GenBank database (http://www.ncbi.nlm.nih.gov/). Using BioEdit version 5.0.9 (Hall 1999Go), ambiguous and uninformative variable sites were excluded and an optimized sequence dataset of length 376 nucleotides was submitted to subsequent phylogenetic analyses. To determine the model of nucleotide substitution that best fits the sequence dataset, Modeltest version 3.06 (Posada and Crandall 1998Go) was used. This practice has two methods: hierarchical likelihood ratio test (hLRT) that examines models if they are accepted or rejected at each level of a hierarchical hypothesis-testing framework and Akaike information criterion (AIC) that rewards models for good fit but penalizes them for unnecessary parameters (Posada and Crandall 1998Go). The nucleotide substitution models selected by hLRT and AIC were compared and the one with a higher likelihood score was chosen in the subsequent phylogenetic analyses. Phylogenetic analyses were performed based on maximum parsimony and Bayesian inference (Huelsenbeck et al 2001Go). Parsimony analysis was conducted with PAUP 4.0b10 (Swofford 2002Go) with tree bisection reconnection (TBR) branch swapping and maxtrees unrestricted. All gaps were treated as missing data. Among most parsimonious (MP) trees generated by parsimony method, the tree with maximum likelihood score based on the selected DNA substitution model (Posada and Crandall 1998Go) was chosen for display and was indicated as best MP.

To evaluate confidence in internal nodes for the best MP tree, 1000 nonparametric bootstrap replications (resampling size = 1 000; TBR; maxtrees, unrestricted) were used (Felsenstein 1985Go). MrBayes version 3.0b4 software (Huelsenbeck and Ronquist 2001Go) was used to reconstruct the Bayesian tree of Fomitopsis and related taxa. MrBayes offers the best fitting tree for a given sequence dataset using MCMC that modifies the tree topology according to posterior probabilities through branch swapping for nearest neighbor interchange (NNI) (Huelsenbeck and Ronquist 2001Go). The Bayesian phylogenetic tree was reconstructed according to these conditions: (i) the model of nucleotide substitution resulted from Modeltest was used to reconstruct Bayesian trees; (ii) four Markov chains (number of generations = 2 000 000; frequency of sampling trees = one per 1000 generations) were set according to specific conditions for running Markov chains (Hall 2001Go); (iii) burn-in was set at 10 000 to discard trees obtained before 10 000th generation when Markov chains reached convergent and stationary likelihood values (Hall 2001Go); (iv) PAUP 4.0b10 was used to reconstruct a MAP tree with maximum posterior probability among 2000 sampled trees; (v) the posterior probability of each node was estimated based on the frequency at which the node was resolved among 2000 sampled trees with the consensus option of 50% majority-rule. The posterior probability of the node is referred to as Bayesian support value (Simmons et al 2004Go). All trees were rooted by two outgroups, L. sulphureus and P. schweinitzii, and the aligned sequences were deposited at TreeBase (SN1721). Constrained trees based on several hypotheses for the phylogenetic relationships of Fomitopsis were constructed with MacClade 4 (Maddison and Maddison 2000Go) and PAUP 4.0b10, and the Kishino-Hasegawa test (Kishino and Hasegawa 1989Go) was performed to evaluate the significance of those hypotheses. The congruence between best MP and MAP tree topologies also was examined by the Kishino-Hasegawa test.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Selection of the model of nucleotide substitution.— – We amplified partial 28S ribosomal DNA from 13 strains of 10 Fomitopsis species and 20 strains of 15 related white and brown rotting fungi and determined their nucleotide sequences. The nucleotide length determined from partial 28S ribosomal DNA was about 600 bp long. From the results of hLRT and AIC examinations in Modeltest, TrN+I+G and TVM+I+G models respectively were selected for the sequence dataset (TABLE IIGo). The two models were not significantly different in prior information of base frequency, gamma distribution shape parameter, proportion of invariable sites and rate matrix (TABLE IIGo). However, the TVM+I+G model selected by AIC method showed a higher likelihood score than that of the TrN+I+G model selected by hLRT, so the former was chosen to be the best fitting evolutionary model for the sequence dataset in this study.


View this table:
[in this window]
[in a new window]
 
TABLE II. Models, likelihood scores and prior informations resulted from MODELTEST for the sequence dataset used in this study

 
Phylogenetic analysis using parsimony method.— – In analyzing phylogenetic relationships of Fomitopsis and related taxa using parsimony method, 75 most parsimonious trees (tree length = 291 steps, CI = 0.491, RI = 0.684) were generated. Through the calculation of likelihood scores for all most parsimonious trees under the TVM+I+G model, one best MP tree with maximum likelihood score was determined and is shown (FIG. 2Go). In this tree Fomitopsis were divided into three clades. Clade A including F. cajanderi and F. rosea formed a group with Antrodia albida (Fr.) Donk, the type species of Antrodia, A. heteromorpha (Fr.) Donk, A. variiformis (Peck) Donk, Daedalea quercina (L.) Pers., Piptoporus portentosus (Berk.) G. Cunn., Postia placenta (Fr.) M.J. Larsen & Lombard and P. rennyi (Berk. & Broome) Rajchenb. Clade B containing F. pinicola, the type species of Fomitopsis, F. feei, F. palustris, F. cupreorosea, F. lilacinogilva, F. spraguei, F. dochmia and Piptoporus betulinus (Bull.) P. Karst. represented a main lineage and was indicated as Fomitopsis core group. Clade C consisted of only one species, F. officinalis. Antrodia xantha (Fr.) Ryvarden was placed at the phylogenetic base for the whole clade of the brown rot members consisting of Antrodia, Daedalea, Fomitopsis, Piptoporus and Postia.



View larger version (41K):
[in this window]
[in a new window]
 
FIG. 2. Best MP tree of partial nuclear large subunit rDNA sequences. Sequences of Laetiporus sulphureus and Phaeolus schweinitzii were used as outgroups to root tree. Fomitopsis taxa within clades are indicated in boldface and marked as black vertical bars and non-Fomitopsis as blank vertical bars. Bootstrap values are shown at nodes supported by more than 50% from 1000 replications, and bold lines are used where branches are supported by more than 90%.

 
Assessment of Fomitopsis monophyly.— – The Kishino-Hasegawa test was used to find out whether Fomitopsis species form a monophyletic group as a whole or in part through the statistical estimation of phylogenetic tree topologies. Eleven different hypotheses (TABLE IIIGo) were submitted to the Kishino-Hasegawa test. When a constrained tree excluding P. betulinus from the core group was constructed to see whether Fomitopsis core group without P. betulinus forms a natural group, the tree was rejected in favor of the unconstrained tree (FIG. 2Go) by 0.016 of P value (TABLE IIIGo). When all taxa of Fomitopsis were forced to form a monophyletic group, the resulted tree also was rejected (P = 0.003) and the unconstrained tree was proved to be best, which suggested that all Fomitopsis members are not clustered as a phylogenetically homogeneous group. The constrained tree excluding F. palustris, the type species of Pilatoporus, from Fomitopsis core group was rejected by P value of 0.000, thus rejecting the taxonomic proposal of Kotlaba and Pouzar (1990)Go to segregate a new genus Pilatoporus from Fomitopsis. For the constrained tree of Fomitopsis core group and F. cajanderi and F. rosea, the Kishino-Hasegawa test again rejected the hypothesis that F. cajanderi and F. rosea could be included in Fomitopsis core group (P = 0.004). The monophyly of six species belonging to the F. rosea complex (Carranza-Morse and Gilbertson 1986Go) evidently was rejected by the Kishino-Hasegawa test in favor of the best MP tree (P = 0.001). For the monophyletic tree of Fomitopsis core group and F. officinalis, the Kishino-Hasegawa test rejected the hypothetical tree (P = 0.006) and indicated that F. officinalis could not be included in Fomitopsis core group.


View this table:
[in this window]
[in a new window]
 
TABLE III. Results of Kishino-Hasegawa tests for the sequence dataset using the bootstrap analysis with full optimization and the two-tailed test

 
Bayesian inference using MCMC computation.— – Priors determined from the examination by the TVM+I+G model (TABLE IIGo) were applied to accomplish Bayesian phylogenetic analyses. The process of MCMC completed the assigned 2 000 000 generations and plotted likelihood scores of 2000 sampled Bayesian trees (FIG. 3Go). The 1 089 000th generation proved to be the MAP point and the 1089th tree among sampled ones stood for the MAP tree of 50% majority-rule consensus (FIG. 4Go). Bayesian support values of the 50% majority-rule consensus tree for 2000 sampled Bayesian trees were marked at nodes. When best MP and MAP trees were compared, some differences in tree topology were noted. Taxa composition of Fomitopsis core group (clade B) for two trees was basically same but F. feei and F. palustris in the best MP tree were clustered with P. betulinus, while those in the MAP tree were clustered with D. quercina and P. portentosus. For the constrained best MP tree where F. feei and F. palustris were grouped with D. quercina and P. portentosus, the Kishino-Hasegawa test (P = 0.924, TABLE IIIGo) showed that the phylogenetic position of D. quercina and P. portentosus was statistically acceptable in both trees.



View larger version (11K):
[in this window]
[in a new window]
 
FIG. 3. Plot of likelihood scores that MCMC generated using the priors determined from the values of TVM+I+G. Range of "a" is the consensus region where the log likelihood scores reached a plateau and "b" the MAP point at the 1 089 000th generation. Number of generations = 2 000 000; frequency of sampling trees = 1000; burn-in = 10 000.

 


View larger version (44K):
[in this window]
[in a new window]
 
FIG. 4. MAP tree of partial nuclear large subunit rDNA sequences constructed by MrBayes. Laetiporus sulphureus and Phaeolus schweinitzii were used as outgroups to root tree. Fomitopsis taxa within clades are indicated in boldface and marked as black vertical bars and non-Fomitopsis ones as blank vertical bars. Bayesian support values of the 50% majority-rule consensus tree for 2000 sampled Bayesian trees are marked at nodes supported by more than 50%. Bold lines show where branches are supported by more than 90%.

 
All nodes of Fomitopsis core group in the best MP tree was fully resolved, but the node of the group comprising of F. feei, F. palustris, D. quercina, P. portentosus, F. pinicola and P. betulinus in the MAP tree was not resolved and formed a trichotomy. There was a branch exchange of F. dochmia for F. spraguei between best MP and MAP trees. When their branches of the best MP tree were exchanged as in the MAP tree, the constrained tree was statistically accepted by the Kishino-Hasegawa test (P = 0.899, TABLE IIIGo). The group of F. cajanderi and F. rosea was placed at the basal position of clade A in the best MP tree but clustered with the subclade of P. placenta, P. rennyi and A. variiformis in the MAP tree. Fomitopsis officinalis formed a distinct clade of its own in both best MP and MAP trees but served as a phylogenetically basal taxon for the brown rot taxa of A and B clades and the rest. A branch exchange of F. officinalis and A. xantha between best MP and MAP trees was noted. The topological differences between both trees for F. cajanderi and F. rosea and for the exchanged positions between F. officinalis and A. xantha all were statistically acceptable (P = 0.139 and P = 0.920, TABLE IIIGo, respectively). However, when the Kishino-Hasegawa test was conducted to compare general topologies of those two trees, the likelihood score (–ln L = 1966.603 16) of the MAP tree was higher than that (–ln L = 1976.017 39) of the best MP tree by the difference of 9.414 23 in –ln L values (TABLE IVGo).


View this table:
[in this window]
[in a new window]
 
TABLE IV. Results of Kishino-Hasegawa tests on MAP and best MP tree topologies

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Fomitopsis species of typical brown rot habit formed a polyphyletic group and were divided and dispersed into three clades, suggesting that Fomitopsis is a phylogenetically heterogeneous and diverse taxonomic group (FIGS. 2Go, 4Go). In broad or narrow senses, various generic concepts have been suggested for Fomitopsis. However, our molecular analyses provided a narrow generic concept based on Fomitopsis species that were clustered as the members of Fomitopsis core group (clade B), suggesting that it is phylogenetically and/or taxonomically appropriate to differentiate F. rosea and F. cajanderi of clade A and F. officinalis of clade C from Fomitopsis at generic level. On the other hand, the segregation of F. palustris from Fomitopsis to designate as the type species of Pilatoporus was rejected by the Kishino-Hasegawa test, showing that F. palustris is a good member of Fomitopsis core group. In addition, the monophyly of six taxa belonging to F. rosea complex was rejected, suggesting that the complex definition based on morphological similarities is not phylogenetically secured. The phylogeny of Bayesian inference had an overall similar configuration of topology to that of maximum parsimony, even though there were certain differences in taxa composition of Fomitopsis core group, branch position of F. dochmia and F. spraguei and clade position of F. cajanderi and F. rosea within clade A.

Fomitopsis core group of clade B.— – In this study, the major group of Fomitopsis consisted of F. pinicola, F. feei, F. palustris, F. cupreorosea, F. lilacinogilva, F. spraguei and F. dochmia, along with P. betulinus in best MP tree and P. betulinus, D. quercina, and P. portentosus in MAP tree. The clade B, Fomitopsis core group, was supported weakly in best MP and MAP trees with less than 50% bootstrap and 50% Bayesian support values (FIGS. 2Go, 4Go). To test whether clade B is homogeneous for Fomitopsis, it needs to be examined whether F. lilacinogilva and P. betulinus belongs to Fomitopsis. Fomitopsis lilacinogilva used to be called T. lilacinogilva (Berk.) Lloyd, and the taxonomic position of F. lilacinogilva has been controversial. However, F. lilacinogilva has annual, biennial or perennial basidiocarps of brown rotting activity in comparison with most Trametes species that usually have annual basidiocarps of white rotting activity. In both best MP and MAP trees, the clade of F. cupreorosea and F. lilacinogilva was well supported by 99% bootstrap and 100% Bayesian values. These morphological and molecular data support the nomenclature of Wright and Deschamps (1975)Go, who renamed T. lilacinogilva as F. lilacinogilva, and suggest that F. lilacinogilva is a natural member of Fomitopsis. In both trees F. lilacinogilva CBS 421.84 and F. cupreorosea CBS 236.87 had no phylogenetic variation from each other, while F. lilacinogilva CBS 422.84 was differentiated phylogenetically from F. lilacinogilva CBS 421.84 and F. cupreorosea CBS 236.87, indicating a possibility that F. lilacinogilva CBS 421.84 could be a misidentified strain of F. cupreorosea.

In Fomitopsis core group, Fomitopsis is not monophyletic without the inclusion of P. betulinus, the type species of Piptoporus, which lends support to the transfer of this species to Fomitopsis. Piptoporus betulinus shares these characters with Fomitopsis: dimitic to trimitic hyphal systems with clamped generative hyphae; mostly ellipsoid to cylindric, hyaline, smooth, and IKI- basidiospores (Donk 1974Go, Gilbertson and Ryvarden 1986Go). However, the basidiocarp of P. betulinus is annual and its substrate is restricted to Betula L., which is quite different from Fomitopsis members with perennial basidiocarps living on substrates of both conifers and hardwoods. Although P. betulinus was clustered weakly with other members of Fomitopsis core group, the Kishino-Hasegawa test showed that the exclusion of P. betulinus from Fomitopsis core group was not acceptable statistically (TABLE IIIGo). The transfer of P. betulinus seemed to be plausible, but there was another Piptoporus, P. portentosus in this Fomitopsis core group and a third Piptoporus, P. soloniensis in clade A, proving that Piptoporus was also heterogeneous. The inclusion of P. betulinus into Fomitopsis certainly requires additional studies with more strains of the species and the genus to draw inclusive conclusions for Piptoporus. Daedalea quercina and P. portentosus were clustered with F. feei and F. palustris in the MAP tree, while these species were grouped with P. rennyi, A. variiformis and P. placenta in the best MP tree. The best MP tree constrained in the same way as the MAP tree was statistically accepted by the Kishino-Hasegawa test (P = 0.924, TABLE IIIGo), suggesting that two different phylogenetic positions of D. quercina and P. portentosus in both trees are dependable topologically.

Kotlaba and Pouzar (1990)Go proposed a new genus, Pilatoporus, and transferred F. palustris into it as a type species based on the presence of pseudoskeletal hyphae with conspicuous clamp connections. Although F. palustris was included in Fomitopsis core group by weak support values (FIGS. 2Go, 4Go), the Kishino-Hasegawa test (TABLE IIIGo) showed that the exclusion of F. palustris from Fomitopsis core group statistically was inappropriate, indicating that pseudoskeletal hyphae described by Kotlaba and Pouzar (1990)Go could not be a useful generic character for the segregation of F. palustris from Fomitopsis core group.

Fomitopsis cajanderi and F. rosea and Antrodia strains of clade A.— – Carranza-Morse and Gilbertson (1986)Go suggested F. rosea complex comprising F. cajanderi, F. cupreorosea, F. dochmia, F. feei, F. lilacinogilva and F. rosea. In best MP and MAP trees, F. feei, F. cupreorosea, F. lilacinogilva and F. dochmia of F. rosea complex were placed in clade B but F. cajanderi and F. rosea were positioned in clade A. In addition, the monophyly of F. rosea complex was rejected by the Kishino-Hasegawa test in the best MP tree (P = 0.001, TABLE IIIGo). These results indicate that F. rosea complex defined by Carranza-Morse and Gilbertson (1986)Go was not statistically and phylogenetically monophyletic. Fomitopsis rosea is closely related to F. pinicola, the type species of Fomitopsis, while F. cajanderi shares some morphological characters with Daedalea, according to Nobles’ (1971)Go study. However, phylogenetic positions for F. cajanderi and F. rosea in best MP and MAP trees (FIGS. 2Go, 4Go) indicate that these species are not closely related to the Fomitopsis core group including F. pinicola but distantly separated from D. quercina, the type species of Daedalea. In addition, the inclusion of F. cajanderi and F. rosea into Fomitopsis core group was rejected statistically by the Kishino-Hasegawa test (P = 0.004, TABLE IIIGo), suggesting that F. cajanderi and F. rosea should be separated from Fomitopsis. Kotlaba and Pouzar (1990)Go transferred F. rosea into their new genus Rhodofomes as a type species, but our findings showed that F. cajanderi and F. rosea were closely related to each other with little genetic variation, indicating that they should be placed in a same genus together. However, to determine whether F. cajanderi should be transferred into Rhodofomes, more taxonomic and phylogenetic studies with additional materials need to be achieved because other species are related to F. cajanderi and F. rosea in clade A.

In both trees (FIGS. 2Go, 4Go), two strains of Antrodia variiformis were clustered respectively with P. rennyi and the group of A. serialis and P. soloniensis. Gilbertson and Ryvarden (1986)Go suggested that the basidiocarp of A. variiformis is similar to that of A. serialis sharing a same type of substrate. Postia rennyi was characterized by the presence of chlamydospores along the margin of the basidiocarp. Although this fungus recently was reported in China (Ryvarden and Gilbertson 1994Go), it mostly is distributed in Europe (Gilbertson and Ryvarden 1987Go). In contrast, A. variiformis is a cosmopolitan fungus with no chlamydospores. Antrodia variiformis CBS 375.82 originally was isolated from gymnosperm wood from New York and deposited in CBS by Butin in 1982 (CBS Filamentous Fungi Database, http://www.cbs.knaw.nl/databases/index.htm). Such data indicate that A. variiformis CBS 375.82 of the USA origin clustered with P. rennyi could be a misidentified (or mishandled) strain. In addition, two strains of A. albida were placed separately as shown in best MP and MAP trees (FIGS. 2Go, 4Go). When our private ITS sequences from these two strains of A. albida (CBS 458.86, CBS 308.82) were compared with the one of A. albida AJ006680 [GenBank] from Yao et al (1999)Go and genetic distances were calculated among three sequences, sequence similarity between AJ006680 [GenBank] and CBS 458.86 was 98.6% while the one between AJ006680 [GenBank] and CBS 308.82 was 94%. The difference between two sequence similarities indicated that the CBS 308.82 strain could be another strain closely related to A. albida or a third that has been misidentified.

Fomitopsis officinalis of clade C.— – The best MP and MAP trees show that the phylogenetic position of F. officinalis is close to the basal lineage of brown rotting fungi comprising Antrodia, Daedalea, Fomitopsis, Piptoporus and Postia. In the Kishino-Hasegawa test (TABLE IIIGo), monophyly of F. officinalis with Fomitopsis core group was rejected. In fact, F. officinalis is a synonym of Laricifomes officinalis (Vill.) Kotl. & Pouzar, the type species of Laricifomes Kotl. & Pouzar. The genus Laricifomes defined by Kotlaba and Pouzar (1957)Go was described on the basis of trimitic hyphal system with hyaline hyphae, presence of clamp connections solely on thin-walled generative hyphae, white chalky context, crumbly consistency, styptic properties and pileal surface without a resinous substance. Among characters describing Laricifomes, the presence of sclerids is an especially important key to separate Laricifomes from Fomitopsis (Kotlaba and Pouzar 1998Go). Sclerids are unique hyphal bodies present in context that are contorted and irregular, often lobed and thick-walled, but their taxonomic importance is not specified fully (Gilbertson and Ryvarden 1986Go). Although F. officinalis has sclerids, there has been confusion whether the hyphal system of F. officinalis is dimitic (Gilbertson and Ryvarden 1986Go, Ryvarden and Gilbertson 1993Go) or trimitic (Kotlaba and Pouzar 1957Go, 1998Go). On the other hand, the morphological characters of F. officinalis match well with the generic concepts of Antrodia, which are characterized mainly by the dimitic hyphal system with clamped generative hyphae and cylindric-ellipsoid, hyaline, thin-walled, smooth and nonamyloid basidiospores (Gilbertson and Ryvarden 1986Go). Ryvarden and Gilbertson (1993)Go described A. macrospora Bernicchia & De Domincis as having sclerid-like elements while the others of Antrodia lack such a structure (Gilbertson and Ryvarden 1986Go, 1993). Although F. officinalis was closely related to A. xantha in best MP and MAP trees where both species are topologically exchanged, F. officinalis was not clustered with any other taxa in this study. When Antrodia including A. macrospora and Laricifomes including L. concavus (Cooke) G. Cunn. and L. maire G. Cunn. are incorporated in a future study, it will be easier to evaluate the taxonomic implication of sclerids as a key character to species and determine the exact phylogenetic position and nomenclatural status of F. officinalis.

Phylogenetic reliability of the Bayesian method.— – The difference of criteria for tree reconstruction is subject to generate different branch patterns between best MP and MAP trees. MP generates the shortest tree with minimal number of total nucleotide substitutions (Felsenstein 1988Go), but MAP tree is obtained by selecting one with maximum posterior probability among generated trees. Maximum parsimony and Bayesian inference are accomplished by different tree search algorithms, heuristic search and MCMC process, respectively. Local maxima obtained by heuristic search of parsimony are different from those obtained by MCMC computation of Bayesian inference (Durbin et al 1998Go), indicating that the difference between criteria eventually generates topological conflicts between best MP and MAP trees. If topological uncertainties at weakly supported nodes are removed through tighter taxon sampling and multi-gene analysis with large character information, topological differences between both trees could be compared more exactly.

Bayesian support values (Bayesian posterior probabilities) and bootstrap values have statistical differences. The Bayesian support value represents an estimate of the probability that one clade is correct under a given substitution model and sequence data (Miller et al 2002Go). In contrast, the bootstrap value is the frequency that the topology of clade appears when repetitive processes are conducted with sequence data resampled from original data (Miller et al 2002Go). Wilcox et al (2002)Go suggested that Bayesian support values provided much closer estimates of phylogenetic accuracy than the estimates provided by corresponding bootstrap proportions. On the other hand, Simmons et al (2004)Go recently indicated that Bayesian support values generally were overestimated in comparison to common resampling-based support measures including the bootstrap and jackknife. Thus, they advised against considering posterior probabilities as universal probabilities, although it is appropriate to consider a posterior probability the probability of truth, and recommended the term posterior probability to refer to the probability conditional on data, models and priors. However, when the Kishino-Hasegawa test determined the topologies of two phylogenetic trees, the MAP tree based on Bayesian inference showed a higher likelihood score than the best MP tree, indicating that the relationships of Fomitopsis species in the MAP tree reflects better evolutionary patterns for the present sequence dataset than those in the best MP tree. In our study, Bayesian support values in the MAP tree were found to be comparatively higher than bootstrap values for the clades in the best MP tree, suggesting that the Bayesian algorithm could be properly applied to the phylogenetic analyses of Fomitopsis groups and also might be able to provide improved inference methods for brown rotting fungi.



View larger version (13K):
[in this window]
[in a new window]
 
FIG. 1. Rate matrices of nucleotide substitutions for hLRT and AIC generated by Modeltest.

 

    ACKNOWLEDGMENTS
 
The authors are sincerely grateful to Dr Soon Gyu Hong (Department of Plant Pathology, University of Arizona, Tucson, Arizona) who carefully read the manuscript and offered valuable comments. This work was supported by the Korea Research Foundation Grant (KRF-2002-CP0377) for the biochemical and phylogenetic analyses of Fomitopsis. Kyung Mo Kim was supported by the BK21 Research Fellowship from the Ministry of Education and Human Resources Development.


    FOOTNOTES
 
Accepted for publication April 20, 2005.

1 Corresponding author. E-mail: minervas{at}snu.ac.kr


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Bondartsev AS, Singer R. 1941. Zur Systematik der Polyporaceen. Ann Mycol 39:43–65.

Carranza-Morse J, Gilbertson RL. 1986. Taxonomy of the Fomitopsis rosea complex (Aphyllophorales; Polyporaceae). Mycotaxon 25:469–486.

Corner EJH. 1989. Ad Polyporaceas VI. The genus Trametes. Nova Hedwigia 97:1–197.

Donk MA. 1964. The conspectus of the families of the Aphyllophorales. Persoonia 3:337–343.

———. 1974. Check list of European polypores. Verh Koninkl Nederl Akad Wet, Natuurk 2, 62:1–469.

Durbin R, Eddy S, Krogh A, Mitchison G. 1998. Biological sequence analysis. Cambridge, UK: Cambridge University Press. 356 p.

Felsenstein J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783–791.[CrossRef]

———. 1988. Phylogenies from molecular sequences: inference and reliability. Annu Rev Genet 22:521–565.[CrossRef][Medline]

Gelman A, Carlin J, Stern H, Rubin D. 1995. Bayesian data analysis. London, UK: Chapman and Hall/CRC. 526 p.

Gilbertson RL, Ryvarden L. 1986. North American polypores. Vol. 1. Abortiporus-Lindtneria. Oslo, Norway: Fungiflora. 433 p.

———, ———. 1987. North American polypores. Vol. 2. Megasporoporia-Wrightoporia. Oslo, Norway: Fungiflora. 451 p.

Hall BG. 2001. Phylogenetic trees made easy: a how-to manual for molecular biologists. Sunderland, Massachusetts: Sinauer Associates. 179 p.

Hall TA. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 41:95–98.

Hibbett DS, Thorn RG. 2001. Basidiomycota: Homobasidiomycetes. In: McLaughlin DJ, McLaughlin EG, Lemke PA, eds. Systematics and evolution. The Mycota VII Part B. Berlin, Germany: Springer-Verlag. p 121–168.

Huelsenbeck JP, Ronquist F. 2001. MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics 17:754–755.[Abstract/Free Full Text]

———, ———, Nielsen R, Bollback JP. 2001. Bayesian inference of phylogeny and its impact on evolutionary biology. Science 294:2310–2314.[Abstract/Free Full Text]

Kishino H, Hasegawa M. 1989. Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in hominoidea. J Mol Evol 29:170–179.[CrossRef][Medline]

Kotlaba F, Pouzar Z. 1957. Notes on classification of European pore fungi. Cesk Mykol 11:152–170.

———, ———. 1990. Type studies of polypores described by A. Pilát—III. Cesk Mykol 44:228–237.

———, ———. 1998. Notes on the division of the genus Fomitopsis (Polyporales). Folia Cryptog Estonica Fasc 33:49–52.

Larget B, Simon DL. 1999. Markov chain Monte Carlo algorithms for the Bayesian analysis of phylogenetic trees. Mol Biol Evol 16:750–759.

Lecellier G, Silar P. 1994. Rapid methods for nucleic acids extraction from Petri dish-grown mycelia. Curr Genet 25:122–123.[CrossRef][Medline]

Lee JS, Ko KS, Jung HS. 2000. Phylogenetic analysis of Xylaria based on nuclear ribosomal ITS1-5.8S-ITS2 sequences. FEMS Microbiol Lett 187:89–93.[CrossRef][Medline]

Maddison WP, Maddison DR. 2000. MacClade: analysis of phylogeny and character evolution. Version 4.0. Sunderland, Massachusetts: Sinauer Associates.

Miller RE, Buckley TR, Manos PS. 2002. An examination of the monophyly of morning glory taxa using Bayesian phylogenetic inference. Syst Biol 51:740–753.[CrossRef][Medline]

Nobles MK. 1971. Cultural characters as a guide to the taxonomy of the Polyporaceae. In: Petersen RH, ed. Evolution in the higher basidiomycetes. Knoxville: The University of Tennessee Press. p 169–196.

Posada D, Crandall KA. 1998. Modeltest: the model of DNA substitution. Bioinformatics 14:817–818.[Abstract/Free Full Text]

Pouzar Z. 1966. Studies in the taxonomy of the polypores II. Folia Geobot Phytotaxon 1:356–375.

Ryvarden L. 1991. Genera of polypores: nomenclature and taxonomy. Syn Fung 5:1–363.

———, Gilbertson RL. 1993. European polypores. Part 1. Abortiporus-Lindtneria. Syn Fung 6:1–387.

———, ———. 1994. European polypores. Part 2. Meripilus-Tyromyces. Syn Fung 7:394–743.

Simmons MP, Pickett KM, Miya M. 2004. How meaningful are Bayesian support values? Mol Biol Evol 21:188–199.[Abstract/Free Full Text]

Swofford DL. 2002. PAUP: Phylogenetic Analysis Using Parsimony. Version 4.0b10. Champaign, Illinois: Illinois Natural History Survey.

Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgin DG. 1997. The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res 24:4876–4882.

Vilgalys R, Hester M. 1990. Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species. J Bacteriol 172:4238–4246.[Abstract/Free Full Text]

White TJ, Bruns T, Lee S, Taylor J. 1990. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, eds. PCR protocols: a guide to methods and applications. San Diego, California: Academic Press. p 315–322.

Wilcox TP, Zwickl DJ, Heath TA, Hillis DM. 2002. Phylogenetic relationships of the dwarf boas and a comparison of Bayesian and bootstrap measures of phylogenetic support. Mol Phylogenet Evol 25:361–371.[CrossRef][Medline]

Wright J, Deschamps J. 1975. Revista de investigaciones agropecuarias. Ser. 5. Patologia Vegetal 12:143.

Yang Z, Rannala B. 1997. Bayesian phylogenetic inference using DNA sequences: a Markov chain Monte Carlo method. Mol Biol Evol 14:717–724.[Abstract]

Yao YJ, Pegler DN, Chase MW. 1999. Application of ITS (nrDNA) sequences in the phylogenetic study of Tyromyces s.l. Mycol Res 103:219–229.[CrossRef]




This article has been cited by other articles:


Home page
MycologiaHome page
K. M. Kim, J. S. Lee, and H. S. Jung
Fomitopsis incarnatus sp. nov. based on generic evaluation of Fomitopsis and Rhodofomes.
Mycologia, November 1, 2007; 99(6): 833 - 841.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kim, K. M.
Right arrow Articles by Jung, H. S.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Kim, K. M.
Right arrow Articles by Jung, H. S.
Agricola
Right arrow Articles by Kim, K. M.
Right arrow Articles by Jung, H. S.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS