tvor/cE, r#,: Mycologia, 95(3), 2003. pp. 480--487. 77#e 17~'Cu-PYrinl~j ,,,ay be © 200.~ by The MycologicalSociety of America, Lawrence, KS 66044-8897 q S. Co,fellah.

Taxonomy of the vinicolor species complex based on analysis of ITS sequences and microsatellite loci

Annette M. Kretzer ~ mushroom genera Suillus, Gomphidius and Chroogom- S. U.N.Y. College of Environmental Science and phus and is thought to be derived from an epigeous Forestry, Faculty of Environmental and Forest Biology, (possibly Suillu~qike) ancestor through loss of com- 350 lllick Hall, Syracuse, Ntno York 13210-2788 plex morphological characters, such as a stipe, verti- Daniel L. Luoma cally oriented tubes and forcible spore discharge .Department of Foresl Science, Oregon State University, (Bruns et al 1989). Because of its reduced morphol- 154 Peary Hall, Corvallis, Oregon 97331-7501 ogy, Rhizopogon is a taxonomically difficult genus and the current has largely been shaped by the Randy MoIina U.S.D.A. Forest Service, Pat@c Northwest Re.warch influential work of Smith and Zeller (1966). In ad- Station, 2300 SWJefferson Way, Corvallis, Oregon dition to describing numerous new species, Smith 97331 and Zeller (1966) also established a detailed subge- neric classification with two subgenera (Rhizopogonel- .Joseph W. Spatafora la 'and Rhizopogon) that comprise two sections (Rhi- De4?artment of Botany and Plant Pathology, Oregon State University, 2082 Cordley Hall, Corvallis, Oregon zopogonella, ~trulatae) and four sections ( Pd~izopogon, 97331-2902 Amylopogon, Fulviglebae, Villosuli), respectively. Spe- cies in the subgenus Rhizopogonella later have been transferred to the genus Alpova (Trappe 1975), and Abstract: We are re-addressing species concepts in what we recognize as genus llhizopogon today is es- the Rhizopogon vinicolor species complex (, sentially Smith and Zeller's subgenus 1U~izopogon plus ) using sequence data from the inter- Alpova olivaceotinclus (from section Rhizopogonella) nal-transcribed spacer (ITS) region of the nuclear ri- which has been transferred back to the genus Rhi- bosomal repeat, as well as genoLypic data from five zopogon by Grubisha et al (2001); it remains to be microsatellite loci. The R. vinicolor species complex determined in future studies, however, if other spe- by our definition includes, but is not limited to, col- cies from Alpova section Rhizopogonella will have to lections referred to as/L vinicolorSmith, R. diabolicus be transferred back to genus Rhizopogon as well. Smith, R. ochraceisporus Smith, R. parvulus Smith or Grubisha et al (2002) recently conducted an ITS tL vesiculosus Smith. Holo- and/or paratype material sequence-based phylogenetic study of the genus Rhi- for the named species is included. AJmlyses of both zopogon. They found that Smith's and Ze[ler's (1966) ITS sequences and microsatellite loci separate collec- section Amylopogon was monophyletic, section Rhizo- tions of the/7_ vinicolor species complex into two dis- pogon was polyphyletic and section Fulviglebae likely tinct clades or clusters, suggestive of two biological would have been polyphyletic, if more species had species that subsequently are referred to as R, vini- been included; those species from section Fulviglebae color sensu Kretzer et al and R. vesiculosus sensu that were included in the study formed a monophy- Kretzer et al. Choice of the latter names, as well as letic clade nested within a paraphyletic section ViUo- morphological characters, are discussed. suli. Based on these results, Grubisha et al (2002) Key words: fimgal species concepts, internal have proposed a number of changes to the subge- transcribed spacers, microsatellite markers, Rhizopo- neric classification of Rhizopogon that include reinsti- gon tution of an emended subgenus Rhizopogon, eleva- tion of sections Amylopogon and Villosuli to subgenus level and creation of two new subgenera, Roseoli and INTRODUC~rlON Versicolores (for more details see Grubisha et al 2002). Rhiz~pogon is an ectomycorrhizal genus in the Bole- Those species of section Fulviglebae that group with tales (Basidiomycota) that forms hypogeous, truffle- former section Villosuli were transferred to the new like fruit bodies. It is closely related to the epigeous subgenus ViUosuli. Subgenus Villosuli is a strongly supported, monophyletic clade that will be the focus Accepted for publication February 5, 2003. of the study presented below. Members of the sub- Corresponding autbor. E-mail: [email protected] genus Villosuli share a high level of host specificity

480 KRETZER ET AI_,: RIIIZOPOGONVINICOLOR SPECIES COMI'LEX 481

and are known to associate only with Pseudotsuga spp. They lend themselves to population genetic work be- (Massicotte et al 1994, Molina et al 1997). Another cause they are large (up to several cm across) and detailed phylogenetic study of North American col- are encountered more frequently than frtfit bodies, lections of subgenus Amylopogon recently has been making them relatively easy to sample in the field. published by Bidartondo and Bruns (2002). Spurred by our interest in population genetics of In addition to subgeneric classification schemcs, this group, we arc re-addressing taxonomic issues and there has been controversy about species delinea- species delineations within the 17~ vinicolor species tions in Rhizopogon. Most species currently recog- complex using both ITS sequences and microsatellite nized in North Aanerica date back to Zeller and markers. Characterization o[ six polymorphic loci Dodge (1918), Zeller (1941) and Smith and Zeller with trinucleotide repeat motils from 1¢. vinicolorwas (1966), with subsequent additions by Smith (1966, reported earlier (Kretzer et al 2000). Included in this 1968), Harrison and Smith (1968), Trappe and Guz- study are scveral I°&iztrpogon spp. type collections m~n (1971), Hosford (1975), Cfizares et al (1992) from which DNA has not been extracted and ana- and Allen et al (1999). However, many species mor- lyzed before. phologically are very similar and described differenc- es might represent morphological and/or ontoge- netic variation within biological species. A number of MA'FERIA[~SAND METHODS synonymies therefore have been proposed (Martin et Type collections were obtained from The University of al 1998). Grubisha et al (2002) also reported a nnm- Michigan Collection. When collections consisted of ber of irregularities with species delineations and several sporocarps, one well-preserved sporocarp was cho- identifications. Of particular interest in the context sen for DNA extraction and amplification. Tubcrculale EM of this study was the t'act that several collections of R. amdyzed in this study were collected in the spring of 1998 vinicolor were not monophyletic but appeared to in two plots, one located at Mary's Peak (44 ° 31.86' N and form a paraphyletic grade. Collections of R. diaboli- 123 ° 32.86' W) in the Oregon Coast Range and the other cus, R. ochraceisptnn~s and R. parvulus were derived at Mill Creek (44° 12.17' N and 122 ° 13.95' W) in the Oregon Cascade Range. Plots were within 40-80- and 40- from within the paraphylctic R. vinicolor grade. We 50-year-old second-growth forests, respectively, dominated collectively will refer to these taxa as the R. vinicolor by Douglas (Pseudotsuga menziesii), weslern hemlock species complex. R. vinicolor and R. ochraceisporus ( 7?~uga heterophylla) and western red cedar ( 7"huja plicata). can be differentiated primarily by subtle color differ- The plotswere 10 m x 10 m, and an II x I1 point grid ences in the reaction of the peridinm to KOH as well was established with 1 m point intervals. Soil cores (5 cm as by an olive-brown versus rusty-brown ("russet") wide and approx. 30 cm deep) were taken at each grid gleba at maturity, according to Smith and Zeller point (121 cores per plot), and soils were sifted through a (1966). In R. diabolicus, the mature gleba is "russet," W. S. Tyler Company 2 mm sieve. [n total, 30 independcnl as in the case of R. ochraceisporus, but maintains a collections of tuberculate EM were ohlatined, 15 fl'om bright rusty-cinnamon color afi.er drying. R. parvulus Mary's Peak and 15 fl'om Mill Creek. Tuberculate EM col- is distinguished by some irregularly shaped spores. lections were fi'eeze-dried be[bi'e molecular analysis and will be identified as MPI-98 (Mary's Peak) and MCI-98 (Mill Finally, we also include 1~ vesiculosus in the R. vini- Creek) ~mples, respectively. Finally, cultures of R. vinicolor color species complex because Smith calls it "scarcely T20787 and T20874 were kindly provided by Dr. Ari Jump- distinguishable" from R. vinicolor in the dried stage; ponen (now at Famsas State University). Voucher collections when fresh, R. vesiculosus is distinguished by yellow- of sporocarps from which these cultures were derived arc brown inflated cells in the epicutis (Smith arid Zellcr housed in the Fungus Collection of the Oregon State Uni- 1966). versity Herbarium (OS()#61147 and OSC#61148). Despite known difficulties with species delineations Genomic DNA was extracted from dried material, as de- in the R. vinicolor species complex, we have chosen scribed in Kretzer et al (2000). The ITS region (comprising to work with R. vinicolor as a model taxon to study ITS-I, the 5.8S rRNA gene, and ITS-If) was l'CR-amplificd the population structure of hypogeous ba~sidiomy- with primers ITSlf and ITS4. When amplification of the ceres (false-truffles). Rhizopogon vinicolor is a pre- entire region was unsuccessfid, only ITS-I was amplified with primers ITSlf and ITS2. For primer sequences, see White dominantly spring-fruiting species (Luoma et al et al (1990) and Gardes and Bruns (1993). PCR reactions 1991) that is very abundant in the Pacific Northwest. contained 50 mM KCI, 10 mM Tris-HCl (pH 9.0 at 25 C), Most important, it forms morphologically distinct ec- 2.5 mM MgC1._,, 0.1% Triton ® X-100, 0.2 mM of each dNTI', tomycorrhizae (EM) on (Pseudotsuga 0.5 IxM of each primer, 50 U/m[. Taq DNA polymerase, and menziesii) known as tuberculate EM (Zak 1971, Mas- empirical amounts of template DNA. Temperature cycling sicotte et al 1992). Tuberculate EM consists of tight conditions were: 2 rain at 94 C, fifllowed by 35--40 cycles of clusters of ectomycorrhizal roots,encased in a weft of 45 sat94C, 30 sat50C, 60s + I s/cycle at 72 Cand a hyphae that often is referred to as the "peridium". final extension of 10 rain at 72 C. PCR products were di- 482 MYCOi,O(;IA gested with AluI and IIinfl to produce restriction fragment- tance, which is [1 - (no. of shared alleles/2 X no. of loci length polymorphisms (RFI.P's) according to Gardcs and compared)] (Bowcock et al 1994). Distances based on step- Bruns (1996). For sequencing purposes, PCR products were wise mutation models were not suitable for our dataset, be- purilied by electrophoresis on 1% agarose gels followed by cause at some loci alleles sizes were not strictly spaced in extraction with the QIAquick Gel Extraction Kit (Qiagen multiples of three nucleotides as would be expected for loci Inc.). Nucleotide sequences were determined with a BigDye with trinucleotide repeats under a stepwise mutation model. Terminator sequencing kit and an ABI 373 automated se- As already was discussed, irreg-ular spacing can be caused quencer. Sequencing Analysis and SeqEd sof~vare was used by small insertions/deletions and/or substitutions in the to process raw data (PE Applied Biosystems). microsatellite flanking regions, and was particularly pro- PCR conditions for the amplification of microsatellite loci nounced when collections belonging to different species as well as primer sequences have been reported before were analyzed togethel: (Kretzer et al 2000). Sizing of PCR product-~ was performed by acrylamide gel electrophnresis on an ABI 377 auu~mated sequencer using the "GS500 Tamra" internal size standard. RESULTS Band sizes were estimated with GeneScan software (PE Ap- plied Biosystems). Because the mobility of DNA fragments A total of 14 new ITS sequences were determined in is inlluenced by base composition as well as the lluorescent this study and are listed in TABLE I with their Gen- dyes used, size units translate only roughly into numbers of Bank accession numbers. Most new sequences were hasepairs and commonly include fractions of a unit. Nev- derived from holo- and paratype material designated ertheless, they are known to be highly reproducible with by Alexander Smith (Smith and Zeller 1966) and standard deviations for fragments in the 200-300 bp range made available to us by The University of Michigan commonly -0.15 (e.g., Haberl and Tautz 1999). Alleles Fungus Collection. Collection AHS68673, the holo- therefore were scored as different whenever a break in the type for R. parvulus, was not available to us tor DNA conli~uous distrihution of allele nletusurements was detect- cal. In some cases, that meant thai. allele si~es might differ extraction because of small collection size, and DNA hy as little as 0.5 size units. Small differences of this kind extraction and amplification from collection can he explained by single basepair insertions and/or sub- AHS65267, the holotype for R. ochraceisporus, was un- stitutions that are known to occur in the microsatellite successful. All tuberculate EM collections were Ilaqking regions (e.g., Orti et al 1997). screened before sequencing by ITS-RFLP's, using the Phy]ogenefic analysis of ITS sequences was condt,cted in PCR primers ITSlf and ITS4 and the restriction en- PAUP* (Sinauer Associates Inc.). Alignment of ITS se- zymes Alul and HinfI. Two different fungal ITS- quences was perJbrmed manually with the PAUP* editor RFLP patterns were observed with the restriction en- at~d a color fi:mt. Sequences immediately adjacent to the zyme Alul, and one exemplary collection of each pat- priming sites of tile ~quencing primers were of poor qual- tern was chosen for sequencing (for details on band ity in some of tile sequences, and the respective areas fi'om the 18S, 5.8S and 28S genes were excluded fi'om the anal- sizes see discussion). ysis. Fore" additional nucleotide positions were excluded be- ITS sequences determined in this study were cause they appeared to be polymorphic, not only between aligned and analyzed together with 18 sequences taxa hut also within individual collections as indicated by from Rhizopogon subgenus Villosuli that were deter- double peaks in the sequencing chromatograms. Transi- mined in a previous study (Grubisha et al 2002). tions and transversions were weighted equally (=1 step). The complete ITS dataset included 32 taxa and 542 Alignmem gaps were treated as missing data, but parsimo- characters, of which 47 characters were parsimony ny-informative alignment gaps were recoded by changing informative. Parsimony analysis resulted in 100 one character per alignment gap to a new state 'T' (lbr most-parsimonious trees that were 81 steps long. "indel'). The new state 'T' was introduced in the same One of the most-parsimonious trees is shown in FIG. position, when alignment gaps were identical in size and positiml and in dil'l~rent positions otherwise. Change fl'om 1A, and branches snpported by all 100 trees are any nt,clcotide to an "I" was weighted as one step. Parsi- highlighted in bold lines. Numbers above branches ,nony analysis was perf~)rmed using the henristic search op- indicate support from bootstrap analysis. Collections tion with 10 random sequence additions. Bootstrap analysis of the R. vinicolorspecies complex fall into two well- was based on 10 000 replicates using the fast stepwise ad- resolved and strongly supported clades that further dition option. Tile dataset minus the exclt, ded positions was on we will refer to as R. vinicolor sensu Kretzer et al d(,posited at TreeBase. and R. vesiculosus sensu Kretzer et al (for choice of Microsatellite data were anal~ed by neighbor-joining clade names, see discussion); there was no resolu- (NJ) analysis of genotypic distances. Allele frequency-based tion between collections within either clade due to analyses were not possible because of the relatively small munber of type collections available fbr many of the species little or no sequence variation. Both R. vinicolorsen- analyzed here and because the group is taxonomically so sn Kretzer et al and R. vesiculosus sensu Kretzer et difficul! that existiqg species concepts call be tested reliably al clades include exemplary sequences from tuber- only cm type material. We used a simple allele sharing dis- culate EM, indicating that taxa in both clades form KRETZER ET AL: RItIZOPOGON V1NICOLOR SPECIES COMPLEX 483

TAm.~; I. Collections from the R_ vinicolorspecies complex studied by either ITS sequence or microsatellite analysis or both. Additional ITS sequences from Rhizopogon subgenus ViUosuli were taken from Grubisha et al. (2002) Collection No. Type Origin ITS accession IL diabolic~ AHS68404 paratype Pend Oreille Co., WA AF366386 IL diabolictL, AHS68424 paratype Pend Oreille Co., WA Grubisha et ai., 2002 IL diabolicTL~ AHS68489 paratype Bonner Co., ID AF366387 tL diabolicus AHS70305 holotype Idaho Co., ID AF366388" IL ochraceisporus AHS58928 paratype Idaho Co., ID AF366389 R. ochraceisporus AHS65963 paratype Adams Co., ID Grubisha et al., 2002 IL ochraceisporu.s AHS68395 paratype Pend Oreille Co., WA AF366390" R. ochraceisporus T17916 -- OR Grubisha et al., 2002 IL ochraceisporus T17944 -- OR Grubisha et al., 2002 R. parvulus AHS68364 paratype Bonner Co., ID Grubisha et al., 2002 IL v~ic~zlosus AHS68040 holotype Bonnet Co., ID AF366391 IL vesiculosus AHS68041 paratype Bonner Co., ID AF366392 IL vinicolor Tru2195 holotype Boise Co., ID AF263929 1L vinicolor AHS68071 paratype Bonner Co., ID AF418268 1~ vinicolor AHS68092 paratype Pend Oreille Co., WA AF263930 R. vinicolor AHS68177 paratype Bonner Co., ID -- FL vinicolor AHS68179 paratype Bonner Co., ID AF263931 tL vinicolor AHS68189 paratype Bonner Co., ID -- R. vinicolor AHS68227 paratype Pend Oreille Co., WA -- R. vinicolor AHS68576 paratype Bonner Co., ID R. vinicolor AHS68579 paratype Bonner Co., ID -- R. vinicolor AHS68595 paratype Bonner Co., ID -- R. vinicolor AHS68602 paratype Bonner Co., ID AF263932 IZ vinicolor AHS68690 paratype Bonner Co., 1D -- R. vinicolor AHS68705 paratype Bonner Co., ID -- R. vinicolor T17899 -- Benton Co., OR Grubisha et al., 2002 R. vinicolor T19383 ~ Coos Co., OR Grubisha et al., 2002 R. vinicolor T20787 ~ Adams Co., ID Grubisha et al., 2002 1.L vinicolor T20874 ~ Bonner Co., ID -- Tuberculate EM MC1-98-H9 = type I ~ Lane Co., OR AF263938 Tuberculate EM MC1-98-B0= type II -- I,ane Co., OR AF263939 Tuberculate EM MC1-98-XX ~ Lane Co., OR -- Tuberculate EM MP1-98-XX -- Benton Co., OR -- * Only the ITSI region was successfully amplified and sequenced; see methods. ec'tomycorrhizae with tuberculate morphology on were genets of R. vinicolor sensu Kretzer et al, and Douglas fir. five were genets of R. vesiculosus sensu Kretzer et al We subsequently have scored microsatellite alleles as indicated by ITS-RFLP analysis conducted earlier at five loci for 22 sporocarp voucher collections from (see first paragraph of the results section). Details on the R. vinicolor species complex (again mostly type the size and distribution of genets across both plots material, see Fie,. 1B and TAm.E I) and 30 collections will be published elsewhere. To avoid redundancy, of tuberculate EM from two plots described in the only one representative tuberculate EM sample of methods. Characterization of six polymorphic micro- each genet was included in the cluster analysis pre- satellite loci from R. vinicolor T20787 was reported sented below; in addition, one genet was omitted previously (Kretzer et al 2000). PCR primers devel- from the analysis because of missing data at one lo- oped for five of the six loci also were found to amplify cus. fragments of expected size from collections of R. ves- Calculation of multilocus pairwise distances (see iculosus sensu Kretzer et al; those were the primers methods) and neighbor-joining analysis resulted in for loci Rv02, Rv15, Rv25, Rv46 and Rv53, which con- the phenogram shown in FIG. lB. It separates collec- sequently were used in this study. Microsatellite mark- tions into two distinct clusters, which correspond to ers revealed that the 30 collections of tubet:culate EM R. vinicolorsensu Kretzer et al and R. vesiculosus sen- made in two 10 m by 10 m plots represented a total su Kretzer et al clades in Fie,. 1A. Within the two of 11 different multilocus genotypes (= genets); six clusters, however, again there is no clear evidence for 484 MYCOLOGIA

A R. vinicolor Tru2195 H R. vinicolor AHS68092 P ITS nucleotide sequences: R. vinicolor T17899 one out of 100 mp trees. R._~.(~hvinicolor T20787 R. ochraceisporus AHS65963 P 1 step raceisporus AHS58928 P clade h 89 tuberculate EM type I R. vinicolor R. diabolicus AHS68424 P sensu Kretzer et al. R. diabolicus AHS70305 H R. diabolic'us AHS68404 P R, parvulus AHS68364 P R. ochraceisporus T17916 B8 R. ochraceisporus T17944 R. ochraceisporus AHS68395 P R. vinicolor AHS68071 P R. vinicolor AHS68602 P R. vinicolor AHS68179 P clade I1: 100 R. vinicolor T19383 R. vesiculosus tuberculate EM type II sensu Kretzer et al. R. vesjculosus AHS68040 H R. vesiculosus AHS68041 P R diabolicus AHS68489 P .~R R. parksii T17679 parksii T19446 Villosulus AHS59143 R. hawkerae AHS68417 P ersii T17228 Ilescens T17681 R. colossus AHS49480 H ~R 76 R, villosulus T19466 • sp. nov. T17466

B Genotypic data from 5 microsatellite loci: NJ phenogram. o. ~L

-- o.o5 %. ,~ ~

" "" " '° Af¢

Rvi T20787

'a~8o8,o "--n:= =0 %~e.o

o FIr,. I. (A) One of 100 most-parsimonious (mp) trees obtained from nucleotide sequences of the internal-transcribed spacer region. Branches that are supported by all 100 mp trees are highlighted in bold; numbers above branches indicate support from bootstrap analysis. The tree is unrooted. (B) Neighbor-joining phenogram based on genotypic data from five microsatellite loci. Abbreviations are: Rd = R_ diabolicus, Ro = R. ochraceisporus, Rp = R. parvulus, Rve = R. vesiculosus, Rvi = R. vinicolor, H = holotype, P = paratype. MC1 and MPI samples are from tuberculate mycorrhizae; see Methods. Arrow indicates branch separating R. vinicolor sensu Kretzer et al and H. vesicula~us sensu Kretzer et al. any sort of sub-clustering that would correlate with et al cluster. This simply reflects the fact that the mi- current taxonomy or otherwise indicate the presence crosatellite markers originally were developed for R_ of multiple biological species. Internal branch vinicolor, and it is not uncommon for microsatellite lengths are longer in the R. vinicolor sensu Ki-etzer markers to be less polymorphic in nontarget species et al cluster than in the R. vesiculosus sensu Kretzer than in target species (e.g., FitzSimmons et al 1995). KRETZER ET AL: RHIZOPOGON VINICOLOR SPECIES COMPLEX 485

DISCUSSION cies names has priority over the others, and we there- fore chose R. vinicolor (sensu Kretzer et al) tor this Both the ITS sequences and the microsatellite loci group, because it currendy is the most widely used analyzed in this study separated collections of the/L name and therefore provides the greatest taxonomic vinicolor species complex into two distinct clades or continuity. clusters indicative of two biological species. These The taxonomic status of R. parvulus is more diffi- findings are in taxonomic conflict with the five or cult to interpret, largely because only a single collec- more species names applied to this group. All species tion (AHS68364 = paratype) was available to us for in question were described by Alexander Smith DNA extraction. This collection clusters tightly with (Smith and Zeller 1966), who intended to document R. vinicolor sensu Kretzer et al in the ITS tree, a find- all observed morphological diversity and to establish ing that is consistent with most morphological char- narrow species concepts, potentially including onto- acters (Smith and Zeller 1966). In the microsatellite- genetic or environmental variants as discrete species based phenogram, however, it falls just slightly out- (Smith and Thiers 1971). To address this possibility, side the R. vinicolor sensu Kretzer et al cluster. It is a re-evaluation of the main morphological characters different from R. vinicolor sensu Smith primarily in used by Smith to separate the species is warranted. A that it has slightly larger (8-10 × 4-5 p.m) and often key character in differentiating R. vinicolor sensu irregularly shaped spores. Even if R. parvulus should Smith and R_ ochraceisporus is the color reaction of constitute an independent species or hybrid, iris un- the peridium to KOH. The differences, however, not likely to confound future population genetic work only are subtle but also are influenced by the devel- because itis very rare. opmental stage. Rhizopogon vinicolorsensu Smith and In both the ITS tree and the microsatellite phen- R_ ochraceisporus also differ in the color of the gleba ogram, a large number of apparently misidentified at maturity, which is "dark olive-brown" in the for- I-L vinicolor collections group strongly and distinctly mer and "rusty cinnamon-brown" (="russet") in the with AHS68489 (discussed above) and latter. Color characteristics of immature fruit bodies R. diabolic'us with two collections of including the are given in much detail for R. vinicolor but are ab- R. vesiculosus, sent from the description of R. ochraceisporus. We holotype (AHS68040). Because again there is no in- think that R. ochraceisporus constitutes a mature color dication for sub-clustering within this group from ei- variant of R. vinicolor. Similarly, R. diabolicus is char- ther dataset, we shall refer to it collectively as R. ves- acterized by a "russet" gleba at maturity that retains iculosus sensu Kre~er et al (Fw,. 1A). As noted by a "bright rusty cinnamon" color when dried. We ex- Snfith and Zeller (1966), R. vesiculosus is "scarcely amined type material for both R. diabolicus distinguishable" from R. vinicolor in the dried stage. (AHS70305 = holo, AHS68404 = para) and R. ochra- This high degree of similarity might explain the large ceisporus (AHS65267 = holo, AHS68395 = para) and number of misidentified collections, not only among found the color of the dried glebas comparable if not Smith's own "vinicolo¢' collections but also among identical. R. diabolicus AHS68489 (=paratype) did OSC herbarium material (data not shown). When not have a "rusty cinnamon" gleba, an observation" flesh, the color of the mature gleba distinguishes that is consistent with its placement within the R. ves- both species, an observation that is consistent with iculosus group by both ITS sequence and microsat- our own (see below). Unfortunately, we have not ellite data (FIG. 1). Mature R. diabolicus sporocarps been able to observe another diagnostic character for strongly resemble R. ochraceisporus ("ochraceous" to R. vesiculosus that Smith and Zeller (1966) describe "rusty brown" peridium staining "rusty brown" in as "yellow-brown (fresh) inflated ceils" found in the KOH and dark "olive" in FeSOt, "russet" gleba), ac- epicutis that are "similar in size and shape to those cording to Smith's description, while immature spo- found in many species of section Villosuli". Although rocarps bear much resemblance to R. vinicolar, in- these cells readily can be observed in dried material cluding the white peridium with "vinaceous" flushes, from section Villosuli, Smith reports that they col- the "lilaceous" or "vinaceous" reaction to KOH and lapse in R_ vesiculosus upon drying and then are "not the weak reaction to FeSO4. Finally, spores of all so readily demonstrable" thus making use of this three species (described by predominant size and ex- character difficult to test. We have not undertaken tremes) are comparable in size (approx. 6.5-9.0 × microscopic studies of this character in fresh material 3.0--4.5 i~m) (Smith and Zeller 1966). In conclusion, but have been unable to observe the inflated cells we believe that R. vinicolor, IL ochraceisporus and R. either in Smith's collections or in our own dried col- diabolicus are synonyms, an interpretation that also is lections. In addition, R. vesiculosus sensu Kretzer et supported by neighbor-joining analysis of multilocus al differs from Smith's description of R. vesiculosus genotypic distances (FIG. 1B). None of the three spe- in having a wider range of spore lengths (approx. (5) 486 MYCOLOGIA

6--9 (10) ~m versus 6--6.5 I~m) and in growing under the R. vinicolor species complex is composed of two Douglas fir rather than lodgepole . sympatrically distributed, phylogenetic species (indi- Based on our extensive sampling of spring-fruiting cated by ITS sequence analyses), which correlate with Rhizopogon species under Douglas fir around Mary's biological species (indicated by microsatellite geno- Peak in the Oregon Coast Range, we find that the typic distances). In future population genetic work, most useful characters for differentiating sporocarps both species readily can be differentiated from either of R. vinicolor sensu Kretzer et al and R. vesiculosus reproductive or vegetative structures using ITS- sensu Kretzer et al are the colors of the fresh perid- RFLP's as described. ium and gleba. However, not all stages of develop- From a taxonomic point of view, we have shown ment are equally distinctive and more than one de- that, in Rhizopogon subgenus ViUosuli, paratype ma- velopmental stage often is needed to unambiguously terial in many cases cannot be relied on to represent assign material to either species: Both species begin the holotype. That puts us in a difficult position for with a white peridium that bruises pinkish-red making taxonomic changes with respect to R. ochra- (Smith's "vinaceous"), but only R. vinicoltrr sensu ceisporus because DNA from the holotype was not suc- Kretzer et al develops vivid yellow patches during ear- cessfully amplified. Two lines of evidence, however, ly maturity. At maturity, both species are light yellow- lead us to believe that R. vinicolor and R. ochraceispo- ish brown (Smith's "ochraceous") and turn various rus should be regarded as synonyms: (i) Morpholog- shades of brown from handling, the shades typically ical characters as discussed above do not provide reflecting the color of the mature gleba underneath strong evidence against synonymy. In particular, we (see below). In R. vinicolor sensu Kretzer et al, the believe that morphology strongly supports synonymy gleba develops from white when immature through of/L ochraceisporus and R. diabolicus (see above); the pale yellow and pale greenish yellow-brown, to dark latter in turn is supported by our molecular data to greenish brown (Smith's "olive-brown") or brown or be synonymous with R. vinicolor. (ii) The three para- more rarely reddish brown (Smith's "rusty cinna- type collections of R. ochraceisporus analyzed here ac- mon-brown" or "russet"). On the other hand, the tually did cluster within the same clade, suggesting gleba of R. ve~siculosus sensu Kretzer et al appears to that paratype material for this particular species develop from white to greenish brown (a stage that might be more consistent than other species. In the apparently is short and was not noted in the type next section, we therefore formally propose synony- description) to dark blackish-brown (dominant mization of R. vinicolor, R. ochraceisporus and/L dia- stage). Finally, we find that tuberculate EM of both bolicus. Although we believe that these changes best species also differ somewhat in morphology. The weft reflect the current state of knowledge, this study does of darkly pigmented hyphae that encases the clus- not claim to be an exhaustive treatment of the R. tered ectomycorrhizae is fluffy in R. vinicolor and at- vinicolor species complex. The study was guided pri- tached to the ectomycorrhizae in such a way that it marily by our desire to clarify species delineations in molds to the ectomycorrhizae and cannot be peeled the Pacific Northwest, which is the geographic center back readily. In R. vesicmlosus, it is appressed but de- of our ongoing population genetic work, and, tached from the ectomycorrhizae such that it can be through incorporation of type material, to provide peeled back in large patches. Although a more de- valid names for the two species identified. The selec- tailed and formalized description of morphological tion of taxa to be included in this study was based differences at the EM level would be desirable, it was largely on pre-existing molecular evidence that sug- not the goal of this study. gested particularly close affiliations of these taxa with Ultimately, sporocarps and tuberculate EM of both I~ vinicolor (Grubisha et al 2002). A filture, more species most readily and reliably can be difterentiated comprehensive study should include vinicolor-like by ITS-RFLP's with the restriction enzyme AluI. "collections from a wider geographic range, as well as When PCR primers ITSlf and ITS4 are used, a single type material from other morphologically similar spe- undigested band of size 743 bp characterizes R. ves- cies, such as R. cinnamoraeus Harrison and Smith, R. iculosus sensu Kretzer et al while three bands of sizes subcinnamorneus Smith, R. olivaceoJitscus Smith, R. 419 bp, 224 bp and 97 bp are most typical of R. vin- pachyspora Hosford and others. icolor sensu Kretzer et al (a C to T transition is oc- casionally observed in one of the restriction sites of TAXONOMY 1~ vinicolor and results in only two bands of sizes 516 bp and 224 bp). Exact band sizes have been deduced Rhizt~pogon vinicolor A. H. Smith. from nucleotide sequences. Basionym: Rhizopogon vinicolorA. H. Smith. In Smith and Our data support the conclusion that, within our Zeller, 1966: Mem. N. Y. Bot. Gard. 14 (2): 67-69. sampling range (Oregon, Washington and Idaho), Synonyms: Rhizopogon ochraceisporus A. H. Smith. In KRFTZER ET AI.: PuHIZOPO(;ONVINICOLOR SPE(:IES COMPI.EX 487

Smith and Zeller, 1966: Mere. N. Y. Bot. Gard. 14 (2): 62. Harrison KA, Smith AH. 1968. Some new species and dis- Rhizopogon diabolicu~$ A. H. Smith. In Smith and Zeller, tribution records of tOzizopol4on ill North Amcrica. Can 1966: Mere. N. E Bot. Gard. 14 (2): 64-65. J Bot 46:881-899. Hosford DR. 1975. Taxonomic studies on the genus Rhizo- pogon I. Two new species from file Pacific Northwest. ACKNOWLEDGMENTS Beih Nova Hedwigia 51:16.'3-169. The authors thank Caprice Rosato fi~r running countless Kretzer AM, Molina R, Spatafi~ra JW. 2000. Microsatellile GeneScan gels, and Nancy Adair and Gi-Ho Sung tbr help markers for the ectomycorrhizal basidiomycete Rhizo- with sequencing. Susie Dunham, Jim Trappe and two anon- pogon vinicohrr. Mol Ecol 9:1190-1191. ymous reviewers have made valuable comments at different Luoma DL, Frenkel RE, Trappe .]M 1991. Fruiting o[' hy- stages of this mannscript's development. This work would pogeous fungi in Oregon Douglas fir forests: seasonal not have been possible without the type material housed at and habitat variation. Myeologia 83:335-353. the University of Michigan Fungus Collection and provided Martin ME H6gbcrg N, Nylund J-E. 1998. Molecular anal- to us by Robert Fogel and Patricia Rogers. The work was ysis confirms morphological reclassification of Rhizo- financed byjoint venture agreement No. PNW-98-5113-1JVA pogon. Mycol Res 102:855-858. and was supported partially by co-op agreement No. 01-CA- Massicotte HB, Melville LIt, Li CY, Peterson RI.. 1992• 11261993-09~PNW, both from the U.S.D.A. Fores! Service Structural aspecls of Douglas fir [Psm~dotsul,m m~.q~ziesii Pacific Northwest Research Station. Christopher Baycura (Mirb.) Franco] tuberculate ecUmlycorrhizae. Trees 6: h~s helped with design of the cover image. 137-146• --, Molina R, Luoma DL, Smith .]E. 1994. Biolog-y of the ectomycorrhizal genus Rhizopogon: II. Patterns of LITERATURE CITED host-dingus specificity following spore inoculation of Allen MF, Trappe JM, Horton TR. 1999. NATS truffle and diverse hosts grown in monoctflture and dual culture. truffle-like fungi 8: Rhizopogon mengei sp. nov. (Bole- New Phytol 126:677-690. taceae, Basidiomycota). Mycotaxon 70:149-152. Molina R, Smith .]E, McKay D, Melville LH• 1997. Biolog), Bidartondo MI, Bruns TD. 2002. Fine-level mycorrhizal of the ectomycorrhizal genus Rhizopogon. Ill. hflluence specificity in the Monotropoideae (Ericaceae): specific- of co-cultured conifer specics on mycorrhizal specificity ity for fungal species groups. Mol Ecol 11:557-569• with the arbutoid hosts Arctostaphylos uva-ursi and A t- Bowcock AM, Ruiz-Linares A, TomfohrdeJ, Minch E, Kidd Imtus menziesii. New Phytol 137:519-528. JR, Cavalli-Sforza IX,. 1994• High resohttion of human Orti G, Pearse DE, Arise JC. 1997. l'hylogenetic assessment evolutionary trees with polymorphic microsatellites. of length variation at a microsatellite locus. Proc Natl Nature 368:455-457. Acad Sci USA 94:10745-10749. Bruns TD, Fogel R, White TJ, PalmerJD. 1989. Accelerated Smith AH. 1966. Ncw and noteworthy higher ttmgi from ew)lution of a false-truffle from a mushroom ancestor. Michigan. Mich Bot 5:18-25. Nature 339:140-142. 1968. Further studies on lOlizopog,'on. I. J Elisha Cfizares E, GarciaJ, Castillo J, Trappc J M. 1992. Hypogeous Mitchell Sci Soc 84:274-280. fungi from northern Mexico. Mycologia 84:341-359. , Thiers HD. 1971. The Boletes of Michigan. Ann FitzSimmons NN, Moritz C, Moore SS. 1995. Conservation A,'bor: University of Michigan Press. 428 p. and dynamics of microsatellite loci over 300 million ~, Zeller SM. 1966. A preliminary account of the years of marine turtle evolution. Mol Biol Evol 12:432- North American species of Rhizopogon. Mere N Y Bot 440. Gard 14:1-178. Gardes M, Bruns TD. 1993. ITS primers with enhanced Trappe JM. 1975. A revision of the genus Alpova wilh notes specificity for basidiomycetes--application to the iden- on Rhizopogon and the Melanogastraceae. Beih Nova tification of mycorrhizae and rusty. Mol Ecol 2:11.%- Hedwigia 51:279-309• 118. ~, Guzmfin G. 1971. Notes on some bypogeous fungi , 1996. ITS-RFLP matching for identitication from Mexico. Mycologia 63:317-332. of fungi. Meth Mol Biol 50:177-186. White "l[J, Bruns TD, l,ee S, Taylor JW. 1991). Amplification Grubisha LC, Trappe .]M, Molina R, Spatafora JW. 2001. and direct sequencing of fimgal ribosomal RNA genes Biology of the ectomycorrhizal genus Rhizopogon. V. tot phylogenetics. In: Innis MA, Gelfand DII, Sninsky Phylogenetic relationships in the Boletales inferred .H, White TJ, eds. PCR protocols; a guide to methods from LSU rDNA sequences. Mycologia 93:82-89. and applications. San l)iego: Academic Press. p 315- .... 2002. Biolog), of the ecto-. 322. mvcorrhizal genus Rhizopogon. VI. Re-examination of Zak B. 1971. Characterization and classification of mycor- infrageneric relationships inferred fi'om phylogenetic rhizas of Douglas fir. 11 Pseudotsuga menziesii + Rhizo- analyses of internal transcribed spacer sequences. My- pogon vinicohn: Can J Bot 4tt:1048-1079. cologia 94:607-619. Zellet" SM. 1941. Further notes on fungi. Mycologia "~'~:19fi- Haberl M, Tautz D. 1999. Comparative allele sizing can pl(~ 214.

duce inaccurate allele size diflerences lot microsatcl- • Dodge (:W. I.t)18. I¢hizo/m,,o, i,i Norl]l Amcrica. lites. Mo[ Ecol 8:1",47-1349. Ann M<~ l~,<)l (;:u'd 5:1-36.