Mycologia
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

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 Shen, Q.
Right arrow Articles by Royse, D. J.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Shen, Q.
Right arrow Articles by Royse, D. J.
Agricola
Right arrow Articles by Shen, Q.
Right arrow Articles by Royse, D. J.
Mycologia, 94(3), 2002, pp. 472-482.
© 2002 by The Mycological Society of America

Molecular phylogenetic analysis of Grifola frondosa (maitake) reveals a species partition separating eastern North American and Asian isolates


Qing Shen
David M. Geiser
Daniel J. Royse 1

     Department of Plant Pathology, The Pennsylvania State University, University Park, Pennsylvania 16802

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 

A phylogenetic analysis was performed on 51 isolates of the commercially valuable basidiomycete, Grifola frondosa (maitake), using sequences from the Internal Transcribed Spacers and 5.8S region of the nuclear ribosomal DNA (rDNA) and a portion of the ß-tubulin gene. The ß-tubulin gene provided more than twice as many variable characters as the ITS/5.8S regions. The isolates analyzed comprised 21 from eastern North America, 27 from Asia, one from Europe, and two of unknown geographic origin, one of which was the major US commercial production strain in use. Grifola sordulenta was used as an outgroup. Combined and separate analysis of both genes showed a partition separating Asian versus eastern North American isolates. Bootstrap analysis showed strong support for these clades in the ß-tubulin data alone and in the combined data. The major commercial isolate of unknown geographic origin is apparently of Asian descent based on its grouping within the Asian clade. The single European isolate analyzed was distinct from both the eastern North American and Asian clades. These results indicate strong support for a species partition separating eastern North American and Asian isolates of G. frondosa, despite previous studies indicating no morphological distinction between them.

Key words: ß-tubulin, biogeography, genetics, ITS, mushroom cultivation, nucleotide variation


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Grifola frondosa (Dickson : Fr.) S.F.Gray is a white rot fungus widely distributed in Asia (Zhao and Zhang 1992Citation ), North America (Gilbertson and Ryvarden 1986Citation ) and Europe (Breitenbach and Kräntzlin 1986Citation ). In North America, it is found throughout regions of the east, Midwest, southeast and the Pacific Northwest. It occurs on a variety of hardwoods, particularly oak and chestnut (Breitenbach and Kräntzlin 1986Citation ). Commonly called maitake or hen-of-the-woods, it is considered a choice edible mushroom with unique culinary and medicinal qualities. Maitake is marketed throughout Asia and, because of increased consumer demand, its commercial production has grown dramatically in Asia and the United States (Chang 1999Citation , Royse 1997Citation ).

Traditional classification of G. frondosa was based solely on morphological characters. Grifola frondosa was first named Boletus frondosus by Dickson (1785)Citation and later, Fries (1821)Citation changed the name to Polyporus frondosus Dicks. : Fr. Even today, P. frondosus, the synonym of G. frondosa, is widely used. The genus Grifola S.F.Gray was first applied by Gray (1821)Citation and described as a polypore with large compound basidiomes. Previous taxonomic investigations by Gilbertson and Ryvarden (1986)Citation and Zhao and Zhang (1992)Citation described similar morphological characters shared between North American and Asian isolates, although these studies separately analyzed North American and mainland China isolates, respectively. Both studies recognized G. frondosa (Dicks. : Fr.) S.F.Gray as the only species in the genus Grifola. However, another Grifola species, G. sordulenta, was identified by Singer in 1969. Recent studies by Hibbett et al (2000)Citation based on mitochondrial and nuclear small and large subunit rRNA gene sequences showed that G. frondosa is related to other polypores including Laetiporus and Ganoderma.

Biogeographic phylogenetic structure is known to exist in various polypores and agarics, such as Pleurotus (Vilgalys and Sun 1994Citation ), Lentinula (Hibbett et al 1998Citation , Thon and Royse 1999Citation ), Panellus stypticus (Jin et al 2001Citation ) and Flammulina (Methven et al 2000Citation ). Genetic selection and improvement of cultivated isolates in commercial mushroom production may be facilitated by a better understanding of the biogeographic phylogenetic structure of G. frondosa germplasm. Molecular approaches have been used extensively for examining phylogenetic relationships in other edible fungi (for example see Hibbett et al 1995Citation , Thon and Royse 1999Citation , and Vilgalys and Sun 1994Citation ).

To analyze biogeographic structure within G. frondosa, partial regions of rDNA and ß-tubulin genes were analyzed and phylogenetic relationships were inferred among isolates from North America and Asia. Results showed a strongly supported partition separating isolates from North America and Asia, suggesting that these represent separate phylogenetic species.


    MATERIAL AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Cultures – Fifty one isolates of G. frondosa and one isolate of G. sordulenta were examined (Table I ), including all available isolates from the Pennsylvania State University Mushroom Culture Collection (PSUMCC) and the American Type Culture Collection (ATCC), as well as an isolate of G. sordulenta (Argentina), which was used as an outgroup. Isolates represented various geographic origins including Asia (27), United States (21), Europe (1), and unknown origin (2). All cultures were maintained on potato dextrose agar supplemented with 1.5g/L of yeast extract (PDYA).


View this table:
[in this window]
[in a new window]
 
TABLE I. List of species, isolate code, source, geographic origin, substrate and locality of Grifola frondosa used for this study

 

View this table:
[in this window]
[in a new window]
 
TABLE I. Continued

 
DNA extraction – Fresh mycelium (100 mg) was used to isolate DNA following the LETS extraction procedure (Chen et al 1999Citation ). Cultures were grown in 50 mL of potato dextrose yeast broth (PDYB) for 20 to 30 d at room temperature. Mycelium was harvested by vacuum filtration on Whatman (grade #1) filter paper, and washed once with distilled water. One-hundred milligrams of fresh mycelium was homogenized with a micropestle in a 1.5-mL microfuge tube containing 400 µL of LETS buffer [20 mM Tris-HCL (pH 7.8), 100 mM LiCl, 10 mM EDTA, and 0.5% SDS, pH 7.8]. Five-hundred mL of a 25:24:1 (vol/vol/vol) mixture of phenol-chloroform-pentanol were added, and the tube was vortexed for 30 s and centrifuged at 14 000 x g for 10 min. The upper portion of the aqueous phase (~250 µL) was recovered, and DNA was precipitated by the addition of 0.5 mL of cold absolute ethanol with incubation at -20 C for 10 min. DNA was collected by centrifugation at 14 000 x g for 10 min, dried, and resuspended in 20 µL of sterile water. DNA preparations were diluted with sterile water and used as template for PCR amplification.

PCR amplification and sequencing – PCR was performed in 25 µL reactions with a 96-well PCR cycler (PTC-100 Programmable Thermal Controller, MJ Research, Inc.), using 10 ng DNA template, one U of Taq DNA polymerase (Promega, Madison, Wisconsin), 0.2 mM of each dNTP, 2 mM MgCl2, 0.1% Triton, and 0.5 µM of each primer. Amplification of ITS-1, ITS-2, and 5.8 S rDNA was performed by utilizing primers ITS1 (White et al 1990Citation ) and ALR0 (5'-CATATGCTTAAGTTCAGCGGG-3') (R. Vilgalys pers comm). PCR reactions for ITS regions were performed using the following parameters: 94°C/1 min; 35 cycles of 94 C/15 s, 60 C/30 s, 72 C/1 min; and 72 C/5 min. PCR Reactions for ß-tubulin regions were performed with primers (designed by authors) BTG5F (5'-CGTTGTGCCCAGTCCTAAGGTG-3') and BTG8R (5'-GTTCTTGCTCTGCACGTTCTG-3') (Fig. 1 ) with the following parameters: 94 C/2 min; 35 cycles of 94 C/15 s, 57 C/30 s, 72 C/1 min; and 72 C/7 min. Amplification products were electrophoresed on a 1.0% agarose gel and checked to ensure that a single DNA band was produced of the expected size (~600bp for ITS PCR products and ~680bp for ß-tubulin PCR products). For sequencing, the PCR products were purified directly from reactions using the Wizard PCR Preps System (Promega Corp., Madison, Wisconsin) and the concentration adjusted to 20 ng/µL. Sequencing reactions were performed using an ABI dye-terminator kit (ABI/Perkin-Elmer) and analyzed using an ABI Prism® Model 377 automated sequencing system (Applied Biosystems, Foster City, California). Sequences have been deposited in GenBank (accession numbers AY049091AY049142 for rDNA ITS sequences, AY049143AY049194 for ß-tubulin sequences).



View larger version (6K):
[in this window]
[in a new window]
 
 FIG. 1. Locations of primers used for PCR-amplification of the ß-tubulin gene in Grifola frondosa and G. sordulenta (G01). Numbers indicate exons

 
Sequence data analysis – Sequence ends were trimmed using the SeqMan II module in the Lasergene package (DNAStar, Inc. Madison, Wisconsin) and adjusted manually. All sequences then were edited and initially aligned using the clustal W algorithm (Higgins et al 1991Citation ) in the Lasergene package (DNAStar, Inc. Madison, Wisconsin). Multiple alignment parameters used were gap penalty = 10 and gap length penalty = 10. Final alignments then were optimized visually. Intron/exon junctions in the ß-tubulin sequences were inferred based on comparisons with the known Schizophyllum commune sequence (Russo et al 1992Citation ). Levels of molecular sequence divergence in each of the data sets were compared by calculating pairwise estimates of nucleotide substitution rates by the two-parameter method of Kimura. Phylogenetic analyses were completed using PAUP Version 4.0b4a (Swofford 2000Citation ). A neighbor-joining (NJ) tree was constructed using the Kimura 2-parameter model. The stability of clades was evaluated by bootstrap tests with 1000 replications (Felsenstein 1985Citation , Hills and Bull 1993Citation ). A maximum parsimony (MP) analysis was performed using 1000 heuristic searches with random taxon addition searches and TBR branch swapping with MAXTREES unrestricted. Other indices for the generated topology, including tree length, consistency index (CI), and retention index (RI) were calculated. A strict consensus of the minimum length MP trees was calculated. Gaps were considered missing data for all analyses. The sequence alignment has been deposited in TreeBASE (accession number M1020).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Analysis of rDNA ITS sequences – Amplification of the ITS-1, ITS-2 and 5.8S ribosomal DNA repeat yielded fragments of approximately 600 bp. Characteristics of nucleotide variation present in these regions of Grifola frondosa and its allies are summarized in Table II . Nucleotide variation among isolates of Grifola frondosa was 5.4% in the rDNA region, and 14.3% between G. frondosa and the G. sordulenta outgroup. Total variable nucleotide sites were 82 with 25 phylogenetically informative sites.


View this table:
[in this window]
[in a new window]
 
TABLE II. Site variation within the ITS-1, 5.8 S, and ITS-2 gene region of Grifola frondosa and Grifola sordulenta

 
A Neighbor-Joining analysis (Fig. 2 ) based on rDNA ITS sequences identified two clades within G. frondosa. The eastern North American clade included all of the U.S. isolates, while the Asian clade consisted of Asian isolates. The eastern North American clade received 83% bootstrap support, but the Asian clade did not receive >50% support. The single European isolate (WC493) fell on a branch basal to the eastern North American clade, but its connection to this clade was not strongly supported (67%). Two isolates of unknown geographic origin (WC685 and WC828, the major US commercial isolate) were placed within the Asian clade. The maximum parsimony (MP) analysis produced 424 equally parsimonious trees (length = 113 steps, consistency index = 0.796, retention index = 0.911), which were similar in topology to the neighbor-joining tree.



View larger version (26K):
[in this window]
[in a new window]
 
 FIG. 2. Phylogenetic analysis of 51 Grifola frondosa isolates based on rDNA ITS sequences using the neighbor-joining method with distance analysis calculated by the Kimura 2-parameter model. Geographic origin is listed beside isolate codes. Numbers on branches represent bootstrap values obtained from 1000 replications (values greater than 60% were shown). Sidebars represent inferred clades based on geographic origin

 
Analysis of ß-tubulin sequences – The characteristics of nucleotide variation present in partial ß-tubulin sequence regions of Grifola frondosa and G. sordulenta are summarized in Table III . Among isolates of Grifola frondosa, nucleotide variation was 12.2% for all sites. Most of the variation occurred in introns with 14.8% for intron 5, 23.7% for intron 6, and 32.3% for intron 7. Less variation was observed in exons with 8.3% for exon 6 and 5.5% for exon 7. The sequences of G. frondosa and G. sordulenta showed much more variation over the total alignment (30.2%) and across the introns and exons. The most variable region was intron 5 (56.9%) and exon 7 was the most conserved region (18.3%) within the ß-tubulin gene sequence regions analyzed. No amino acid changes were inferred among isolates of G. frondosa, but two amino acid changes (glutamate and isoleucine in G. frondosa, but glutamine and valine in G. sordulenta, both in exon 6) were observed between isolates of G. frondosa and G. sordulenta. Most nucleotide substitutions in exons were observed in the third codon positions. Variable nucleotide sites totaled 177. The phylogenetically informative sites totaled 62 with 37 sites from introns.


View this table:
[in this window]
[in a new window]
 
TABLE III. Site variation within the partial ß-tubulin gene region of Grifola frondosa and Grifola sordulenta

 
A Neighbor-Joining analysis (Fig. 3 ) based on partial ß-tubulin gene sequences gave similar results to those based on the rDNA dataset, showing two distinct Eastern North American and Asian clades within G. frondosa. The Eastern North American clade received 99% bootstrap support, whereas the Asian clade received 74% support. In contrast to the rDNA tree, where the European isolate grouped basal to the Eastern North American clade, the ß-tubulin tree placed the European isolate (WC493) on a branch basal to the American clade, with 95% bootstrap support. As in the rDNA analysis, the two isolates of unknown geographic origin grouped with the Asian isolates. The maximum parsimony (MP) analysis produced 48 100 trees (length = 251 steps, consistency index = 0.793, retention index = 0.932) until PAUP was aborted because of lack of memory. The strict consensus of these trees was similar in topology to the neighbor-joining tree. A phylogenetic analysis without intron 5 was also performed because of the high variation in this region. No appreciable difference was found (data not shown).



View larger version (30K):
[in this window]
[in a new window]
 
 FIG. 3. Phylogenetic analysis of 51 Grifola frondosa isolates based on partial ß-tubulin gene sequences using the neighbor-joining method with distance analysis calculated by the Kimura 2-parameter model. Geographic origin is shown beside isolate codes. Numbers on branches represent bootstrap values obtained from 1000 replications (values greater than 70% were shown). Sidebars represent inferred clades based on geographic origin

 
Analysis of combined rDNA ITS and partial ß-tubulin gene sequence data – A Neighbor-Joining analysis (Fig. 4 ) based on combined rDNA ITS and partial ß-tubulin gene sequence data supported most of the results produced by rDNA and ß-tubulin separately, with much higher bootstrap support. The Eastern North American and Asian clades were strongly supported by high bootstrap values (100% and 98%, respectively). The European isolate (WC493) grouped basal to the Asian clade with 89% bootstrap support, and agrees with the results from the ß-tubulin data alone. The two isolates of unknown geographic origin (WC828 and WC685) grouped within the Asian clade. The maximum parsimony (MP) analysis of the combined dataset produced 405 equally parsimonious trees (length = 382 steps, consistency index = 0.741, retention index = 0.901). The strict consensus tree (not shown) retained a similar topology to the NJ trees.



View larger version (29K):
[in this window]
[in a new window]
 
 FIG. 4. Phylogenetic analysis of 51 Grifola frondosa isolates based on a combination of rDNA and partial ß-tubulin gene sequences using the neighbor-joining method with distance analysis calculated by the Kimura 2-parameter model. Geographic origin is shown beside isolate codes. Numbers on branches represent bootstrap values obtained from 1000 replications (values greater than 80% were shown). Sidebars represent inferred clades based on geographic origin

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The Asian and eastern North American isolates appear to form discrete groups, but it is not known whether these groups overlap in Europe, which has yet to be sampled intensively. The 51 isolates of G. frondosa analyzed in this study clearly clustered into two clades based on eastern North American and Asian origins in all analyses of rDNA, ß-tubulin, and combined sequences. Molecular data presented here resolve groups that have not been distinguished based on morphology (Gilbertson and Ryvarden 1986Citation , Zhao and Zhang 1992Citation ). These studies included descriptions of basidiomes, basidiospores, habitats and hyphal context systems. However, neither study compared North American and Asian isolates side-by-side. No mating tests between the North American and Asian isolates have been conducted to determine if they could be different biological species.

The name Grifola frondosa (Dicks.) Gray was applied to the basionym Boletus frondosus J. Dicks., which was described based on a specimen collected in Yorkshire, England (Dickson 1785Citation , Gray 1821Citation ). Although only a single, non-type European isolate was analyzed in the current study, the phylogenetic data suggest the possibility that a distinct European lineage may exist in G. frondosa. If this is the case, then phylogenetic taxonomic revisions would require that the name G. frondosa be applied to the European lineage, with different names for the eastern North American and Asian clades. Taxonomic changes must await study of the type material and an expanded sampling of isolates from Europe and the British Isles. Grifola frondosa is common in Europe (Breitenbach and Kränzlin 1986Citation ), but WC493 was the only culture available for this study. The ß-tubulin analysis and combined analysis both showed WC493 to form a distinct basal branch related to the Asian clade. However, because most of the phylogenetically informative sites came from the ß-tubulin sequence, this inferred relationship between WC493 and the Asian clade cannot be considered strong. Isolates from the Pacific Northwest of North America and the Southern Hemisphere were not sampled in this study. Further study of Grifola spp. and relatives from these areas may uncover additional biogeographic patterns.

Continental phylogeographic structure is common in basidiomycete macrofungi. Intersterility groups in the Pleurotus ostreatus species complex correlate well with continental biogeography, and also with an ITS rDNA phylogeny (Vilgalys and Sun 1994Citation ). Phylogenetic partitions may exist that do not correlate with intersterility barriers. In Lentinula (shiitake), strong continental phylogeographic structure is evidenced based on rDNA phylogenies, including fairly distinct Asian and Eastern North American clades (Hibbett et al 1995Citation , Hibbett et al 1998Citation ). Our results are consistent with the proposal that a biogeographic connection exists between basidiomycete macrofungi from eastern Asia and temperate North America (Wu and Mueller 1997Citation ). In the genus Suillus, this observation is supported by the inference that North American species tend to have closely related Asian sister taxa, based on ITS rDNA data (Wu et al 2000Citation ). The split-gill fungus Schizophyllum commune shows worldwide interfertility and a limited degree of phylogeographic structure in its cosmopolitan range (Raper et al 1958Citation , James et al 1999Citation , James et al 2001Citation ).

Multilocus phylogenetics provides a powerful tool for the recognition of species, as phylogenetic partitions shared among different loci indicate a historical reproductive barrier between clades (= genealogical concordance phylogenetic species recognition or GCPSR, Taylor et al 2000Citation ). The shared partition between Asian and Eastern North American isolates indicated by both genes suggests that there is a reproductive barrier between these groups, meeting the criteria of GCPSR. However, we cannot say whether the reproductive barrier is intrinsic (i.e., reflecting intersterility) or extrinsic (i.e., reflecting geographic separation). Because phylogenetic partitioning can precede the evolution of intersterility, it cannot be assumed that Asian and Eastern North American isolates of G. frondosa are intersterile. Indeed, strong geographic and phylogenetic partitions, albeit inferred from single genes, can be observed within intersterility groups in the genus Pleurotus (Isikhuemhen et al 2000Citation , R. Vilgalys, pers comm), and despite worldwide interfertility (Raper et al 1958Citation ), some degree of continental biogeographic structure can be found in the split-gill fungus Schizophyllum commune (James et al 1999Citation , James et al 2001Citation ). Recombination may occur among isolates within the two geographic lineages of G. frondosa, but evidence for that cannot be extrapolated from these data. Neither the ß-tubulin nor the ITS rDNA trees show much strongly supported structure within lineages, and the partition homogeneity test (PHT) or incongruence length difference (ILD) test cannot be applied to these datasets to test for historical recombination because of relatively high levels of homoplasy (Farris et al 1994Citation , Koufopanou et al 1997Citation , Barker and Lutzoni 2000Citation ).

Both neighbor joining (NJ) and most parsimonious (MP) trees derived from all DNA datasets revealed a consistent grouping of U.S. commercial cultivar WC828 in the Asian clade. This suggested that WC828 has an Asian origin and is closely related to Asian commercial cultivars. It is known that molecular data can be effectively used to select unique shiitake genotypes for evaluation of biological efficiency, quality and average weight of mushrooms (Diehle and Royse 1986Citation , Levanon et al 1993Citation ). So a better understanding of the phylogenetic relationships of Grifola frondosa may help the selection and breeding of commercial lines and help to improve commercial cultivation of these mushrooms.

The ß-tubulin gene region provided more than twice as many phylogenetically informative nucleotide sites as did the ITS rDNA region, in a similar-sized amplicon. This result is consistent with findings in other fungi (e.g., O'Donnell et al 2000Citation ) and shows that ß-tubulin may be a more powerful tool for analyzing the intraspecific phylogeography of basidiomycete macrofungi than ITS rDNA. Other partial protein-coding genes such as translation elongation factor 1-{alpha} may also be extremely useful.


    FOOTNOTES
 
1 Corresponding author, djr4{at}psu.edu Back

Accepted for publication October 1, 2001.


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Barker FK, Lutzoni FM., 2000 Spurious rejection of partition homogeneity by the ILD test: a simulation study Inoculum 51:16

Breitenbach J, Kräntzlin F., 1986 Fungi of Switzerland, Vol. 2. Non-gilled fungi Lucerne: Verlag Mykologia. 412 p

Chang ST., 1999 World production of cultivated edible and medicinal mushrooms in 1997 with emphasis on Lentinus edodes (Berk.) Sing. in China Int J Med Mushrooms 1:291-300

Chen X, Romaine CP, Ospina-Giraldo MD, Royse DJ., 1999 A polymerase chain reaction-based test for the identification of Trichoderma harzianum biotypes 2 and 4, responsible for the worldwide green mold epidemic in cultivated Agaricus bisporus Appl Microbio Biotech 52:246-250

Dickson J., 1785 Fasciculus Plantarum Crytogamicarum Britanniae Londini: Nicol

Diehle DA, Royse DJ., 1986 Shiitake cultivation on sawdust: evaluation of selected genotypes for biological efficiency and mushroom size Mycologia 78:929-933

Farris JS, Kallersjo M, Kluge AG, Bult C., 1994 Testing significance of congruence Cladistics 10:315-319

Felsenstein J., 1985 Confidence limits on phylogenies: an approach using the bootstrap Evolution 39:783-791

Fries EM., 1821 Systema Mycologicum: Sistens Fungorum Ordines, Genera et Species, Huc U.S.que Cognitas, Quas ad Normam Methodi Naturalis Determinavit, Disposuit Atque Descripsit Lundae: ex officina Berlingiana. 355 p

Gilbertson RL, Ryvarden L., 1986 North American Polypores Norway: Fungiflora, Oslo. p 332–333

Gray SF., 1821 A natural arrangement of British plants, according to their relations to each other as pointed out by Jussieu, De Candolle, Brown, C including those cultivated for use, with an introduction to botany, in which the terms newly introduced are explained London: Baldwin, Cradock and Joy. 643 p

Hibbett DS, Gilbert LB, Donoghue MJ., 2000 Evolutionary instability of ectomycorrhizal symbioses in basidiomycetes Nature 407:506-508[Medline]

Hibbett DS, Hansen K, Donoghue MJ., 1998 Phylogeny and biogeography of Lentinula inferred from an expanded rDNA dataset Mycol Res 102:1041-1049

Hibbett DS, Fukumasa-Nakai Y, Tsuneda A, Donoghue MJ., 1995 Phylogenetic diversity in shiitake inferred from nuclear ribosomal DNA sequences Mycologia 87:618-638

Higgins DG, Bleasby AJ, Fuchs R., 1991 CLUSTAL W: improved software for multiple sequence alignment CABIOS 8:189-191[Abstract/Free Full Text]

Hills DM, Bull JJ., 1993 An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis Syst Biol 42:182-192

Isikhuemhen OS, Moncalvo JM, Nerud F, Vilgalys R., 2000 Mating compatibility and phylogeography in Pleurotus tuberregium Mycol Res 104:732-737

James TY, Moncalvo JM, Li S, Vilgalys R., 2001 Polymorphism at the ribosomal DNA spacers and its relation to breeding structure of the widespread mushroom Schizophyllum commune Genetics 157:149-161[Abstract/Free Full Text]

James TY, Porter D, Hamrick JL, Vilgalys R., 1999 Evidence for limited intercontinental gene flow in the cosmopolitan mushroom Schizophyllum commune Evolution 53:1665-1677

Jin J, Hughes KW, Petersen RH., 2001 Biogeographical patterns in Panellus stypticus Mycologia 93:309-316

Koufopanou V, Burt A, Taylor JW., 1997 Concordance of gene genealogies reveals reproductive isolation in the pathogenic fungus Coccidioides immitis Proc Natl Acad Sci USA 94:5478-5482[Abstract/Free Full Text]

Levanon D, Rothschild N, Danai O, Masaphy S., 1993 Strain selection for cultivation of shiitake mushrooms (Lentinus edodes) on straw Bioresource Tech 45:9-12

Methven A, Hughes KW, Petersen RH., 2000 Flammulina RFLP patterns identify species and show biogeographical patterns within species Mycologia 92:1064-1070

O'Donnell K, Nirenberg HI, Aoki T, Cigelnik E., 2000 A multigene phylogeny of the Gibberella fujikuroi species complex: detection of additional phylogenetically distinct species Mycoscience 41:61-78

Raper JR, Krongelb JS, Baxter MG., 1958 The number and distribution of incompatibility factors in Schizophyllum commune Am Nat 92:221-232

Royse DJ., 1997 Specialty mushrooms and their cultivation Hort Rev 19:59-97

Russo P, Juuti JT, Raudaskoski M., 1992 Cloning, sequence and expression of a ß-tubulin-encoding gene in the homobasidiomycete Schizophyllum commune Gene 119:175-182[Medline]

Singer R., 1969 Mycoflora Australis Lehre: Cramer. p 382–383

Swofford DL., 2000 PAUP*. Phylogenetic analysis using parsimony (*and other methods), 4.0b4a Sunderland, Massachusetts: Sinauer Associates

Taylor JW, Jacobson DJ, Kroken S, Kasuga T, Geiser DM, Hibbett DS, Fisher MC., 2000 Phylogenetic species recognition and species concepts in fungi Fungal Gen Bio 31:21-32

Thon MR, Royse DJ., 1999 Evidence for two independent lineages of shiitake of the Americas (Lentinula boryana) based on rDNA and ß-tubulin genes sequences Mol Phylogen Evol 13:520-524[Medline]

Vilgalys R, Sun BL., 1994 Ancient and recent patterns of geographic speciation in the oyster mushroom Pleurotus revealed by phylogenetic analysis of ribosomal DNA sequences Proc Natl Acad Sci USA 91:4599-4603[Abstract/Free Full Text]

White T, Bruns T, Lee S, Taylor J., 1990 Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics In Innis et al, eds. PCR protocols: a guide to methods and applications. New York: Academic Press. p 315–322

Wu QX, Mueller GM., 1997 Biogeographic relationships between the macrofungi of temperate eastern Asia and eastern North America Can J Bot 75:2108-2116

Wu QX, Mueller GM, Lutzoni FM, Huang YQ, Guo SY., 2000 Phylogenetic and biogeographic relationships of eastern Asian and eastern North American disjunct Suillus species (Fungi) as inferred from nuclear ribosomal RNA ITS sequences Mol Phyl Evol 17:37-47[Medline]

Zhao JD, Zhang X., 1992 The polypores of china Berlin: Bibliotheca Mycologica. J. Cramer. p 209–210




This article has been cited by other articles:


Home page
MicrobiologyHome page
L. Yu, B. Zhang, G. R. Szilvay, R. Sun, J. Janis, Z. Wang, S. Feng, H. Xu, M. B. Linder, and M. Qiao
Protein HGFI from the edible mushroom Grifola frondosa is a novel 8 kDa class I hydrophobin that forms rodlets in compressed monolayers
Microbiology, June 1, 2008; 154(6): 1677 - 1685.
[Abstract] [Full Text] [PDF]


Home page
MycologiaHome page
J. Nuytinck, A. Verbeken, and S. L. Miller
Worldwide phylogeny of Lactarius section Deliciosi inferred from ITS and glyceraldehyde-3-phosphate dehydrogenase gene sequences.
Mycologia, November 1, 2007; 99(6): 820 - 832.
[Abstract] [Full Text] [PDF]


Home page
MicrobiologyHome page
G. I. Zervakis, J.-M. Moncalvo, and R. Vilgalys
Molecular phylogeny, biogeography and speciation of the mushroom species Pleurotus cystidiosus and allied taxa
Microbiology, March 1, 2004; 150(3): 715 - 726.
[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 Shen, Q.
Right arrow Articles by Royse, D. J.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Shen, Q.
Right arrow Articles by Royse, D. J.
Agricola
Right arrow Articles by Shen, Q.
Right arrow Articles by Royse, D. J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS