Ectomycorrhizal community structure in a xeric ... - Wiley Online Library

Research

Ectomycorrhizal community structure in a xeric Quercus woodland based on rDNA sequence analysis of sporocarps and pooled roots Blackwell Publishing Ltd

Matthew E. Smith1, Greg W. Douhan1,2 and David M. Rizzo1 1

Department of Plant Pathology, University of California at Davis, Davis CA 95616, USA; 2Department of Plant Pathology, University of California at

Riverside, Riverside CA 92521, USA

Summary Author for correspondence: Matthew E. Smith Tel: +1 530 7520862 Fax: +1 530 7525809 Email: mesmith@ucdavis.edu Received: 21 November 2006 Accepted: 19 January 2007

• Quercus woodlands are key components of California’s wild landscapes, yet little is known about ectomycorrhizal (EM) fungi in these ecosystems. • We examined the EM community associated with Quercus douglasii using sporocarp surveys and by pooling EM roots and subjecting them to DNA extraction, polymerase chain reaction (PCR), cloning, restriction fragment length polymorphism (RFLP) screening and DNA sequencing. • Ectomycorrhizal root symbionts were sampled four times in 2003–04. During this time, the below-ground community structure was relatively stable; we found no evidence of taxa adapted to winter or spring conditions and only one species varied widely in occurrence between years. • The EM community from sporocarps and roots was diverse (161 species), rich in Ascomycota (46 species), and dominated by fungi with cryptic sporocarps. This included a large number of resupinate and hypogeous taxa, many of which were detected both above- and below-ground. Key words: ectomycorrhiza, fungal community, hypogeous fungi, Mediterranean climate, Pezizales, Quercus (oak), Sierra Nevada, temporal variation. New Phytologist (2007) 174: 847–863 © The Authors (2007). Journal compilation © New Phytologist (2007) doi: 10.1111/j.1469-8137.2007.02040.x

Introduction The last decade has seen explosive growth in the number of ectomycorrhizal (EM) fungal community ecology studies (Horton & Bruns, 2001). The heightened interest is partly due to the decreased cost and increased availability of powerful molecular tools first employed by Gardes & Bruns (1996) to elucidate below-ground EM dynamics. Since then studies have used similar techniques to clarify various aspects of EM ecology on roots, including effects of disturbance ( Taylor & Bruns, 1999), nitrogen deposition (Avis et al., 2003), soil type (Gehring et al., 1998), coarse woody debris (Tedersoo et al., 2003), host plant (Kennedy et al., 2003) and succession (Nara et al., 2003). Despite this increase in research, several gaps or limitations remain in our understanding of ectomycorrhizal community ecology. For example, many EM fungi on roots remain taxo-

www.newphytologist.org

nomically unknown when sporocarps cannot be matched to roots using restriction fragment length polymorphism (RFLP) of the internal transcribed spacer (ITS) region (Horton & Bruns, 2001). Although ITS sequencing is now commonly used for identification, this region provides poor resolution for fungal groups underrepresented in public databases (e.g. hypogeous fungi) (Izzo et al., 2005b). Two recent EM community studies used rDNA sequences from both ITS and the large subunit (28s) to clarify the phylogenetic position of fungi not clearly resolved by ITS alone (Tedersoo et al., 2003; Cline et al., 2005). Because of the conserved nature of the 28s rDNA, this approach provides higher-level taxonomic information for unknowns while retaining the enhanced resolution provided by ITS. Ectomycorrhizal studies also suffer from the common practice of assessing the diversity of sporocarp-producing EM fungi

847

848 Research

based primarily on macroscopic, epigeous fruiting bodies. Whether sporocarps are used to generate ITS profiles for comparison to roots (Yamada & Katsuya, 2001) or to detect ecological patterns (Richard et al., 2004), reported species assemblages may provide a skewed view of EM communities because they disregard important hypogeous and resupinate taxa as well as most Ascomycota (Köljalg et al., 2000; Weiss et al., 2004; Izzo et al., 2005b; Tedersoo et al., 2006). Comparisons between EM roots and primarily epigeous sporocarps may contribute to the widely held view that there is a disconnect between taxa on roots and those present as reproductive structures (Izzo et al., 2005b). Another common limitation of EM studies has been the use of single sampling dates to assess the diversity and dynamics of below-ground fungal communities; such studies can only provide ‘snapshots’ of communities and cannot fully account for temporal dynamics (Buée et al., 2005). Bruns (1995) hypothesized that EM communities may respond rapidly to seasonal fluctuations in temperature, moisture, and light. While we know that there can be significant species turnover at an annual scale in some ecosystems (Izzo et al., 2005a) and that EM communities change radically over longer time scales (Taylor & Bruns, 1999), few studies have tracked EM fungi over shorter times. Much of our knowledge about fungal communities on EM roots is inferred from studies of conifer-dominated ecosystems (Horton & Bruns, 2001 and references therein) but three recent studies have used molecular methods to examine EM fungal associates of the geographically widespread and ecologically diverse genus Quercus (Avis et al., 2003; Richard et al., 2005; Walker et al., 2005). In California, Quercus-dominated woodlands and savannas are key ecosystems, covering approximately one-third of the land area and hosting the greatest species diversity of any biome in the state (Standiford, 2002). However, few of California’s Quercus-dominated habitats are protected by inclusion in natural reserves, and Quercus spp. are under increasing pressure from human activities such as fire suppression, urban development and introduced pathogens (Pavlik et al., 1991; Standiford, 2002). Although Quercus are obligate ectomycorrhizal hosts (Smith & Read, 1997), little is known about their native EM mycota in California. No previous studies have examined EM communities on Quercus roots in California but data from taxonomic studies of sporocarps suggest that EM fungi associated with Quercus are diverse, with a high number of putative endemics (Parks, 1921; Nishida, 1989; Smith et al., 2006a; Smith et al., 2006b). In this study we characterized the EM fungal community associated with a drought-tolerant, deciduous, California endemic oak, Quercus douglasii, based on rDNA sequencing from sporocarps and pooled EM roots. To document the temporal dynamics of the EM fungal community, we sampled EM roots on four dates over 2 yr and collected sporocarps over 4 yr. We attempted to accurately depict the EM fungal sporocarp community by sampling equitably from different morphological types including both hypogeous and resupinate taxa.

New Phytologist (2007) 174: 847–863

Materials and Methods Study site All research was conducted in the Koch Natural Area at the University of California Sierra Foothill Research and Extension Center (UCSFREC) in the Sierra Nevada foothills of Yuba County, CA, USA (39°17′ N, 121°17′ W). The site has a Mediterranean climate characterized by cool, wet winters and hot, dry summers with terrain that consists of hills 400– 600 m above sea level. Precipitation generally occurs between October and May (annual mean 71 cm, range 23–132 cm). Temperature varies greatly between seasons (annual mean 17.8°C, range 10–43°C). Soils are derived from basic metavolcanic rocks and are classified as fine, mixed, active thermic Mollic Haploxeralfs (Argonaut series) from the Alfisol order. Soils are low in carbonates and pH ranges between 5.7 and 6.2. The site is dominated by Quercus douglasii Hook & Arn. (blue oak) but two other EM species, Quercus wislizeni A. DC. (interior live oak) and Pinus sabiniana Douglas (foothill pine), are locally abundant. Little disturbance has occurred since 1960. Before that time the area was subjected to small-scale logging, thinning, cattle grazing and frequent fires (McClaran, 1988, UCSFREC:http://danrrec.ucdavis.edu). EM root sampling A 32 × 32 m plot was established in a blue oak woodland to investigate EM community structure over time at a small spatial scale (Fig. 1). Forty-one Q. douglasii were the only trees and known EM hosts on the plot, although Q. wislizeni and P. sabiniana were present in the surrounding woodland. The nearest alternative EM host was a mature Q. wislizeni c. 15 m from the plot edge. Based on fire history and tree sizes, we estimated the stand of Q. douglasii to be ≥ 50 yr old. However, because diameter at breast height and age are weakly correlated in Q. douglasii, the exact age of the stand is not known (Standiford, 2002). Understory vegetation consisted of assorted nonectomycorrhizal herbaceous plants and scattered Toxicodendron diversilobum. A thin litter layer (mean depth 0.75 cm, range 0–2 cm) was irregularly distributed below the trees, but otherwise no soil stratification was noted in the upper 40 cm of soil. We sampled soil cores at two different times (winter and spring) in two different years (2003 and 2004) for a total of four sampling dates. The 32 m2 plot was subdivided into 16 8 × 8 m squares (Fig. 1) and on 5 March 2003 we randomly sampled one soil core from within each of the 16 sampling squares. We returned on three subsequent dates (21 May 2003, 29 February 2004, 3 May 2004) and randomly sampled a core from each 8 × 8 m square within a 25 cm radius from the first core (Fig. 1). This yielded 64 soil cores, of which 59 were successfully polymerase chain reaction (PCR)-amplified and used for both EM community and spatiotemporal analysis.

www.newphytologist.org © The Authors (2007). Journal compilation © New Phytologist (2007)

Research

Fig. 1 Diagram of sampling design for ectomycorrhizal (EM) roots of Quercus douglasii with 32 m2 plot subdivided into 8 × 8 m squares (1–16). Light grey area within square 7 was intensively sampled. Open circles represent the approximate location of Q. douglasii trees. Dark grey 25-cm radius circle depicts where EM cores were extracted at random locations within each sampling square at four times (winter and spring; 2003 and 2004).

On each of the four dates we also sampled 8–10 soil cores from a 2 × 2 m area within square 7 (Fig. 1), yielding a total of 35 soil cores. Data from these 35 soil cores were used to characterize the EM community on roots of Q. douglasii but were not used for spatiotemporal analysis. Since weather and soil conditions varied between years, sampling dates were selected based on empirical observations of soil moisture, root activity, and above ground phenology of Q. douglasii. Winter samples were collected after significant rainfall occurred, but before leaf production by Q. douglasii. This period was selected because exploratory sampling showed that few EM roots were active in early winter (M. Smith, pers. obs.). Sampling during the summer drought was impractical because of hardening of the upper soil layers. Spring samples were collected after leaf out and flowering of Q. douglasii, when the soil was rapidly drying but not yet impenetrable (see the Supplementary Material, Fig. S1). Litter was removed and soil cores of c. 900 cm3 (8.5 cm deep × 11.5 cm diameter) were collected, stored at 4°C and processed within 15 d. Shallow soil layers were examined because fine roots of Q. douglasii are most abundant in the top 20 cm (Milikin & Bledsoe, 1999). During soil collection, care was taken to minimize soil disturbance and replace leaf litter since any disruption could affect the remaining EM community. Soil samples were gently wet-sieved over a 0.12 mm sieve with deionized water and all detectable fine roots were removed with forceps under a dissecting microscope (×10–×20 magnification). Roots were trimmed into 6-cm sections and gently rinsed. Almost all live root tips were ectomycorrhizal; roots with signs of decay were discarded. The first 400–500 healthy EM tips encountered under the microscope were severed from

fine roots, cleaned, and placed in a Petri dish with water. Roots were considered healthy and colonized by EM fungi when they lacked root hairs and displayed signs of a fungal mantle (color change or emanating hyphae). To achieve randomized sampling, the Petri dish was viewed under the dissecting microscope and swirled in a circular motion to redistribute the EM tips. The first three healthy EM tips encountered were placed in a 1.5 ml microcentrifuge tube. The process was repeated until c. 100 root tips were obtained. EM tips were rinsed again with clean water then 100 root tips were lyophilized together and frozen at −80°C until DNA extraction. Sporocarp sampling Sporocarps of EM fungi were collected within a c. 200 m radius of our EM sampling plot on c. 50 expeditions between December 2000 and April 2005. Sporocarps were categorized into three morphological types: epigeous, hypogeous, and resupinate. Taxa whose sporocarps were found at the soil surface, but below the litter, were considered hypogeous regardless of morphology. To collect cryptic sporocarps, a garden cultivator was used to remove litter and soil beneath randomly selected, mature Quercus spp. Cryptic sporocarps were not collected inside the root plot because we did not want to disturb EM fungi on roots. For molecular analysis, a small piece of fresh tissue was removed from each sporocarp and stored in cetyltrimethylammonium bromide (CTAB) buffer at −20°C until DNA extraction. Voucher specimens of sporocarps have been deposited at the herbaria of UC Berkeley (UC) and Oregon State University (OSC) in the USA.

© The Authors (2007). Journal compilation © New Phytologist (2007) www.newphytologist.org


849

850 Research

Molecular protocols Sporocarp samples were ground with a micropestle and DNA was extracted with a modified cetyltrimethylammonium bromide (CTAB) method (Gardes & Bruns, 1993) or Qiagen Stool Kit (Qiagen Inc., Valencia, CA, USA). PCR was performed with forward primers ITS1F and LROR in various combinations with reverse primers ITS4, ITS4B and LR3 (White et al., 1990; Gardes & Bruns, 1993; Hopple & Vilgalys, 1994). The PCR conditions were optimized for each primer set, but the general protocol was: 94°C for 5 min followed by 25 cycles of 1 min at 94°C, 1 min at 55°C and 2 min at 72°C, followed by 7 min at 72°C. The PCR products were visualized on 1.5% agarose gels with SYBR Green I (Molecular Probes, Eugene, OR, USA). When necessary, PCR products were gel purified with BioRad Freeze ‘n’ Squeeze Tubes (Bio-Rad, Hercules, CA, USA) and reamplified with appropriate primers. The PCR products were cleaned with Microcon-PCR Spin Tubes or Montage PCR96 Filter Plates (Millipore, Billerica, MA, USA). Sequencing was performed as needed with the same primers as above with the ABI Big Dye Terminator Sequencing Kit (v3.1). Sequences were read with an ABI13730xl capillary sequencer (Applied Biosystems, Foster City, CA, USA) at the CAES Genomics Facility, UC Davis. Sequences were edited with SEQUENCHER v.4.1 (Gene Codes Inc., Ann Arbor, MI, USA). Pooled roots were treated differently from sporocarps. Before analysing roots, we tested three primer combinations – (i) ITS1F and LR3, (ii) ITS1F and LR5, and (iii) ITS1F and nLSU1221r (Hopple & Vilgalys, 1994; Schadt et al., 2003) – to determine whether different primers detected the same EM fungi. Although all primer sets recovered most of the same EM fungi in preliminary tests, we used the ITS1F/LR3 primer set because it detected several rare species not found by the others (M. Smith, unpublished). Pooled lyophilized roots were ground with a micropestle and DNA extracted using a modified CTAB method (Gardes & Bruns, 1993) and then purified with a MO-BIO Soil DNA kit (MO-BIO Laboratories, Solana Beach, CA, USA). We used 1–4 µl of DNA template and high resolution Taq polymerase (Invitrogen, Carlsbad, CA, USA) for PCR from roots. The PCR protocol was similar to that used for sporocarps except it had an initial denaturation step of 10 min, included 20 instead of 25 cycles and had 4-min extension steps instead of 2-min extension steps. Amplicons were visualized as above. Four successful replicate reactions per sample were mixed together then cloned with a TOPO-TA kit (Invitrogen). At least 48 successful clones per reaction were grown overnight in Luria–Bertani media amended with 100 µg ml−1 of ampicillin. We selected 48 clones because previous work by Landeweert et al. (2003) and Anderson et al. (2003) found that analysis of 30–50 clones was sufficient to detect most common species of fungi from diverse DNA mixtures. Cloned fragments were reamplified in PCR reactions with c. 0.5 µl of the bacterial suspension as template with the


same protocol as for sporocarps. Amplicons were digested separately with the restriction enzymes AluI and HinfI and run simultaneously on a 1.5% agarose gel. Paired digests of each clone were run in adjacent lanes and visualized together. Gel photographs were manually scored to detect all unique RFLP patterns within a core; patterns were not compared across cores. One to four representative clones of each RFLP type in each core were sequenced with ITS1F following the protocol used for sporocarps. Further sequencing was performed as needed with ITS4, LROR, and LR3. DNA sequence analysis and taxonomic determination of EM fungi The ITS sequences were examined by BLAST searches against GenBank and an in-house database using BIOEDIT (Altschul et al., 1997; Hall, 1999). Since almost half of the species detected on roots were also present as sporocarps, we used a traditional species concept to designate taxon names. The EM root sequences matching sporocarps were designated either by a genus and species epithet or, in cases where species identity was unclear, by a genus name followed by an ‘SRC’ (Sierra Research Center) or ‘SF’ (Sierra Field) voucher number. Sequences found on roots but not as sporocarps were named based on data from BLAST searches, sequence alignments and advice from taxonomic experts. We used a conservative naming approach whereby a sequence with affinities to two genera within a family was designated by the family name, followed by a number. For example, a sequence similar to Peziza and Hydnotryopsis, but not matching a sporocarp, was designated as Pezizaceae followed by a number. To confirm identifications, unknown sequences were aligned with known sequences in CLUSTALX (Chenna et al., 2003) and MACCLADE (Swofford, 2001). When ITS sequences did not provide sufficient information, we sequenced representative clones with LR3 and/or LROR and then subjected them to similar analyses. For most taxonomic groups, ITS sequences were considered representative of the same species if they differed by < 3% across the length of ITS1, 5.8 s, and ITS2. This criterion assumes that a 0.2–1.2% error rate may be generated by PCR, cloning, and unidirectional sequencing, and that ITS heterogeneity of c. 1.5% exists at small spatial scales for some species (Horton, 2002; Izzo et al., 2005a). This criterion is also similar to that of previous studies (Izzo et al., 2005a; O’Brien et al., 2005). A few phylogenetic groups had either a large number of ITS types (e.g. Thelephoraceae) or appeared to represent species complexes (e.g. G. harknessii). For these groups we designated names based on previously identified lineages (Douhan & Rizzo, 2005; Smith et al., 2006a), well supported groups, as identified by neighbor-joining analyses in PAUP* 4.0b10 (Swofford, 2001), and consultation with taxonomic experts. To detect PCR artifacts and sequences of nonEM fungi, we used a combination of sequence alignments and BLAST searches.


Research

The ITS sequences found in multiple cores or matching sporocarps were considered authentic sequences but root sequences detected only once were considered possible chimeras. These sequences were subjected to BLAST followed by alignment with taxonomically similar clones and a second alignment with all sequences from the same core. It is difficult to detect chimeras comprised of DNA from closely related taxa (Qiu et al., 2001), so it is possible that the methods described did not detect all chimeras in two speciose groups, the family Thelephoraceae and the genus Inocybe. For these groups, we assembled alignments with all sequences and searched for chimeras using the MAXCHI and CHIMAERA programs of RDP (recombination detection program) (Martin & Rybicki, 2000). Sequences were designated nonEM taxa if they: (i) showed low BLAST affinity to known EM fungi; (ii) showed high BLAST affinity to saprobic or parasitic fungi or nonfungal sequences; and (iii) were not identified as chimeras. Diversity and cluster analyses We used four diversity measures to characterize the community of EM fungi on roots. Richness (S) is the total number of species. Simpson’s Diversity Index (1–D) is a measure of the ‘likelihood that two randomly chosen individuals will be of different species’; values of 1–D range between 0 (no diversity) and 1 (infinite diversity). Pielou’s J is a measure of evenness, if all species occur at the same proportion then J = 1. The Shannon–Wiener Index (H) is a measure that accounts for evenness and richness. We did not consider EM root tips from a core as independent units and thus did not consider clones from a core as independent units (see Cline et al., 2005). Accordingly, we calculated diversity indices based on frequency rather than abundance. Frequency is defined as the number of soil cores in which a species or group was detected and relative frequency is the total number of occurrences of a species or group divided by the total number of occurrences of all taxa. Abundance data was only used as a secondary measure for calculating the importance of different EM groups. Abundance is defined as the total number of clones of a species or group and relative abundance is defined as the total number of clones of a species or group divided by the total number of clones of all taxa. Hierarchical cluster analysis is an agglomerative method that creates groups from multivariate data using a dissimilarity matrix. We used cluster analysis to categorize soil cores based on species composition of EM fungi. We limited cluster analysis to cores sampled at the largest spatial scale (n = 59, Fig. 1); intensively sampled cores from square 7 were not included. We tested whether species composition was influenced more by date or geographic position. Because cluster analyses are influenced by the selected linkage method and distance measure, we performed cluster analysis with two distance measures (Sørrenson and Euclidean distance) and four linkage methods (flexible Beta, centroid, group average, and Ward’s method).

The Sørrenson distance is influenced by rare species while Euclidean distance downweights rare species (McCune & Mefford, 1999). We selected four linkage methods that compared groups based on different criteria and produced dendrograms with varying levels of chaining. Because Ward’s method and the Sørrenson distance are incompatible, Ward’s method was only used with Euclidean distance (McCune & Mefford, 1999). Statistical analysis of clustering was performed by MRPP (Multi-Response Permutation Procedure). This is a nonparametric statistic that uses a Monte Carlo method to examine whether user defined groups are statistically different from random association. Since there are no distributional assumptions, MRPP is appropriate for community analysis (McCune & Mefford, 1999). We tested whether data grouped by each of four categories (core position, year, season, or sampling time (year × season)) were statistically different from random associations. Group weighting was calculated using the default method of PC-ORD with 1000 permutations. All diversity indices, cluster analyses, and statistical tests were performed with PCORD version 4.20 (McCune & Mefford, 1999).

Results Fungal community on EM roots During winter and spring of two consecutive years, 9400 EM root tips of Q. douglasii were sampled from 94 soil cores in a 32 m2 plot. Restriction fragment length polymorphisms were used to screen 4532 clones; 4473 of which represented EM taxa. We recovered a total of 92 EM taxa on roots (1–14 spp. per core, mean 5.56 spp. per core) (Table 1). From each of the four sampling times we analysed 20–26 soil cores and recovered 42–59 EM taxa. No arbuscular mycorrhizal (AM) fungi were detected. The EM community was dominated by several frequently occurring species with a large number of rare taxa (Table 1) and the species accumulation curve showed a characteristic exponential rate of increase similar to other EM ecosystems (data not shown). Commonly used diversity indices (Shannon–Wiener Diversity Index (H) = 1.077, Evenness ( J) = 0.63, and Simpson’s Diversity Index (1–D) = 0.525) suggest that overall diversity was high. The ITS sequences from 42 of 92 taxa (45.7%) were matched to locally collected sporocarps (sclerotia for Cenococcum), and the identities of most other taxa on roots were well-resolved to the genus or species level (Table 1). However, several taxa were only resolved at the Family or Order level. Of the 92 taxa, 36 were Ascomycota (39.1%) and 56 were Basidiomycota (60.9%). Species in the families Thelephoraceae, Sebacinales, Tuberaceae, and Cortinariaceae were the most frequent and abundant. Thelephoraceae, Pyronemataceae, and Cortinariaceae were particularly speciose with 17, 17, and 16 species, respectively (Table 2). Taxonomic groups tended to be ranked similarly whether by relative frequency or relative abundance, although



851

852 Research


Table 1 Ectomycorrhizal (EM) fungi detected as sporocarps or on the roots of Quercus douglasii


Taxon Name

Voucher1

Accession

EM cores2

Phylum3

Sporocarp type4

rDNA

Closest BLAST match5

% Similarity

Amanita sp. (715) Amanita lanei Amanita novinupta Amanita vaginata Amanita velosa Balsamia sp. Boletus amygdalinus Boletus flaviporus Boletus zellerii Cantharellus californicus nom. prov. Cenoccocum geophilum (Lineage 1)7 Cenoccocum geophilum (Lineage 2)7 Cenoccocum geophilum (Lineage 3)7 Clavariadelphus sp. (src121) Clavulina cf. cristata Clavulina cf. rugosa Clavulina sp. (SF-11) Clavulina sp. 1 Cortinarius cf. glaucopus

src715 src437 src669 src385 src431 src395 src491 src443 src444 src46

DQ974689 DQ974693 DQ974690 DQ974691 DQ974692 DQ974730 DQ974705 DQ974706 DQ974704 DQ974688

0 0 0 0 0 0 0 0 0 0

Basidio Basidio Basidio Basidio Basidio Asco Basidio Basidio Basidio Basidio

E E E E E H E E E E

ITS ITS ITS ITS ITS ITS ITS ITS ITS ITS

– – – – – – – – – –

– – – – – – – – – –

–

1

Asco

U/A

ITS, SSU intron

–

–

30

Asco

U/A

ITS, SSU intron

–

–

12

Asco

U/A

ITS, SSU intron

–

–

src121

AY818603, AY818562 AY818592, AY818560 AY818587, AY818564 DQ974709

0

Basidio

E

ITS

430/450 (95%)

src75 src661 SF-11 – –

DQ974710 DQ974712 DQ974711 EF025112 DQ974723

0 0 4 1 1

Basidio Basidio Basidio Basidio Basidio

E E E E E

ITS ITS ITS ITS ITS

Cortinarius cf. flexipes (src48)

src48

DQ974714

0

Basidio

E

ITS


src6

DQ974715

0

Basidio

E

ITS


src94

DQ974716

0

Basidio

E

ITS

Cortinarius cf. flexipes ITS type 1

–

DQ974724

6

Basidio

E

ITS


–

DQ974725

4

Basidio

E

ITS


–

DQ974722

1

Basidio

E

ITS


–

DQ974726

1

Basidio

E

ITS

Clavariadelphus truncatus (DQ097871) Clavulina sp. (DQ202266) Clavulina sp. (AY918958) Clavulina cristata (AJ889929) Clavulina sp. (AY918958) Cortinarius glaucopus (AY174785) Cortinarius cedriolens (AY083179) Cortinarius flexipes (AJ889971) Cortinarius hinnuleus (AY083184) Cortinarius paleaceus (AJ889974) Cortinarius flexipes (AJ889971) Cortinarius cf. decipiens (AJ889947) Cortinarius flexipes (AJ889971)

– –

578/604 (95%) 485/529 (91%) 538/587 (91%) 565/566 (99%) 586/604 (97%) 367/379 (96%) 476/498 (95%) 502/525 (95%) 440/475 (92%) 630/654 (96%) 349/359 (97%) 573/593 (96%)


Table 1 continued

Voucher1

Accession

EM cores2

Phylum3

Sporocarp type4

rDNA


% Similarity

Cortinarius cf. glaucopus

src262

DQ974717

0

Basidio

E

ITS

586/604 (97%)

Cortinarius sp. (src174)

src174

DQ974719

0

Basidio

E

ITS

Cortinarius sp. (src28) (subgen Phlegmacium) Cortinarius sp.1 (subgen. Phlegmacium) Cortinarius cf. trivialis Elaphomyces muricatus Entoloma sp. (src741)

src28

DQ974718

0

Basidio

E

ITS

Cortinarius glaucopus (AY174785) Cortinarius caninus (U56024.1) Cortinarius citrinus (AY174825)

–

DQ974721

1

Basidio

E

ITS

src611 src641 src741

DQ974720 DQ974740 DQ974694

0 0 0

Basidio Asco Basidio

E H E

ITS ITS ITS

Entoloma sp. (src742)

src742

DQ974695

2

Basidio

E

ITS

Ericoid ascomycete19 Gautieria cf. crispa Genabea cerebriformis Genea arenaria Genea gardneri

– src627 src637 src398 src831

1 3 4 1 4

Asco Basidio Asco Asco Asco

U/A H H H H

ITS, 28s ITS ITS, 28s ITS, 28s ITS, 28s

Genea cf. harknessii (lineage1)8 Genea cf. harknessii (lineage2)8 Genea cf. harknessii (src616)8 Genea cf. harknessii (src665)8 Genea cf. harknessii (src830)8 Genea sp. (src680)8

– – src616 src665 src830 src680

5 7 7 2 3 0

Asco Asco Asco Asco Asco Asco

H H H H H H

ITS, 28s ITS, 28s ITS, 28s ITS, 28s ITS, 28s ITS, 28s

– – – – – –

– – – – – –

Genea sp. (src680-like)8

–

3

Asco

H

ITS, 28s

–

–

Geopora cooperii var. gilkeyi Gilkeya compacta8 Gymnomyces xerophilus Hebeloma cf. oculatum

src515 src718 src433 src875

DQ974827 DQ974732 DQ206864 DQ206845 DQ206850, DQ206965 DQ218301 DQ218300 DQ206857 DQ206859 DQ206861 DQ206849, DQ206967 DQ206855, DQ206970 DQ974731 DQ206862 AY603102 DQ974696

Cortinarius citrinus (AY174825) Cortinarius favrei (AF325575) – Entoloma abortivum (AF357019) Entoloma prunuloides (DQ457633) Capronia pilosella (AF050255) – – – –

0 0 0 1

Asco Asco Basidio Basidio

H H H E

ITS ITS, 28s ITS ITS

– – – 613/630 (97%)

Hebeloma cf. velutipes

src303

DQ974697

2

Basidio

E

ITS

Helotiales19

–

DQ974826

1

Asco

U/A

ITS, 28s

Helotiales29

–

DQ974824

2

Asco

U/A

ITS

Helvella cf. crispa

src194

–

0

Asco

E

–

– – – Hebeloma oculatum (AY311525) Hebeloma velutipes (AY818351) Meliniomyces bicolor (AY394885) Hydrocina chaetocladia (AY789413) –

586/612 (95%) 629/650 (96%) 629/650 (96%) 396/401 (98%) – 556/583 (95%) 429/462 (92%) 552/599 (92%) – – – –

593/620 (95%) 587/614 (95%) 310/331 (93%) –

Research


Taxon Name

853

854 Research


Table 1 continued


Taxon Name

Voucher1

Accession

EM cores2

Phylum3

Sporocarp type4

rDNA


% Similarity

Helvella cf. lacunosa6

–

2

Asco

E

ITS, 28s

src667 – src701

1 2 0

Asco Asco Asco

H H H

ITS ITS ITS

Helvella lacunosa (U42681) – – –

367/375 (97%)

Hydnobolites californicus Hydnobolites sp. 1 Hydnoplicata purpurea nom. prov. Hygrophorus cf. gliocyclus Hygrophorus cf. roseiburneus Hygrophorus chrysodon Hygrophorus eburneus Hynotryopsis setchellii Hysterangium cf. separabile Hysterangium sp. (src642) Inocybe cf. armeniaca Inocybe cf. fraudans Inocybe cf. godeyi (src252) Inocybe cf. maculata Inocybe cf. pudica Inocybe cf. pusio (src527) Inocybe cf. rimosa (src514) Inocybe cf. sororia

DQ974831 DQ974834 DQ974733 DQ974818 DQ974734

src744 src29 src782 src13 src870 src649 src642 src225 src156 src252 src516 src554 src527 src514 src60

DQ974728 DQ974727 DQ974809 DQ974729 DQ974735 DQ974810 DQ974736 DQ974803 DQ974804 DQ974805 DQ974806 DQ974741 DQ974742 DQ974801 DQ974802

0 0 1 1 0 4 0 0 5 3 1 0 0 0 9

Basidio Basidio Basidio Basidio Asco Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio

E E E E H H H E E E E E E E E

ITS ITS ITS, 28s ITS ITS ITS ITS ITS, 28s ITS, 28s ITS, 28s ITS, 28s ITS, 28s ITS, 28s ITS, 28s ITS, 28s

– – – – – – – – – 556/572 (97%) – 594/612 (97%) – – 603/623 (96%)

Inocybe sp. 1 Inocybe sp. 2 (cf. sindonia) Inocybe sp. 3 (cf. lanatodisca) Inocybe sp. 4 (cf. flocculosa)

– – – –

DQ974813 DQ974812 DQ974811 DQ974815

11 21 18 2

Basidio Basidio Basidio Basidio

E E E E

ITS, 28s ITS, 28s ITS, 28s ITS, 28s

Inocybe sp. 5 Inocybe sp. 6 Laccaria cf. bicolor

– – src65 646/646 100% src671 src439 src743

DQ974816 DQ974743 DQ974698

1 1 2

Basidio Basidio Basidio

E E E

ITS, 28s ITS ITS

– – – – – – – – – Inocybe godeyi (AY038316) – Inocybe pudica (AY228341) – – Inocybe sp. PBM 2206 (AY732213) Inocybe godeyi (AY038316) Inocybe sindonia (AY380393) Inocybe lanatodisca (AY380382) Inocybe cf. flocculosa (AF335452) Inocybe godeyi (AY038316) Inocybe relicina (AF325664) Laccaria cf. bicolor (AY228356)

DQ974699 DQ974744 DQ974747

0 1 0

Basidio Basidio Basidio

E E E

ITS ITS ITS

Laccaria laccata (AF204814) – –

592/615 (96%) – –

src438 src299 src731 src419

DQ974746 DQ974745 DQ974825 DQ974737

1 0 8 0

Basidio Basidio Asco Basidio

E E H H

ITS ITS ITS, 28s ITS

– – – –

– – – –

Laccaria cf. laccata Lactarius alnicola Lactarius pallescens var. pallescens Lactarius substriatus Lactarius xanthogalactus Marcelleina sp. (src731) Melanogaster cf. tuberiformis

– – –

616/640 (96%) 649/684 (94%) 579/617 (93%) 1209/1235 (97%) 595/615 (96%) 201/206 (97%)


Table 1 continued

Voucher1

Accession

EM cores2

Phylum3

Sporocarp type4

rDNA


% Similarity

Octavianina sp. 1 Otidea umbrina Pachyphloeus sp. 1

– src520 –

2 1 1

Basidio Asco Asco

H E H

ITS ITS ITS

– src696 src683 – – –

1 2 12 3 1 1

Asco Asco Asco Asco Asco Asco

H H H U/A UA UA

ITS, 28s ITS, 28s ITS, 28s ITS, 28s ITS ITS, 28s

Phialophora sp. 19

–

DQ974822

1

Asco

U/A

ITS, 28s

Pyronemataceae1 Pyronemataceae2 (cf. Genea)

– –

DQ974751 DQ974832

4 1

Asco Asco

U/A U/A

ITS, 28s ITS, 28s

Octavianina sp. (AY918953) Otidea umbrina (AF086586) Pachyphloeus sp. (AY920528) Pachyphloeus sp. (AY920528) – – Scabropezia flavovirens (AY500556) Pezizales EM (AJ893245) Marcelleina pseudoanthracina (AY500538) Phialophora verrucosa (AF050283) Otidea umbrina (AY500540) Humaria hemsiphaerica (AY789389)

535/550 (97%) 515/515 100% 426/453 (94%)

Pachyphloeus sp. 2 Peziza ‘hypogeous’ sp. Peziza infossa Pezizaceae1 (cf. Scabropezia) Pezizaceae2 Pezizaceae3

DQ974748 DQ974738 DQ974830 DQ974833 DQ974819 DQ974739 DQ974817 DQ974749 DQ974687 DQ974750

Pyronemataceae3 (cf. Genea)

–

DQ974835 DQ974752

1

Asco

U/A

ITS, 28s

515/558 (92%)

Pyronemataceae4 (src840) Ramaria sp. (src827)

src840 src827

DQ974823 DQ974713

9 0

Asco Basidio

H E

ITS, 28s ITS

Russula albidula Russula basifurcata Russula brevipes var. acrior Russula cf. amoenolens Russula cf. crassotunicata Russula cremoricolor Russula cyanoxantha Russula exalbicans Russula rhodopoda Russula sp.1 (grisea group) Russula sp.2 (grisea group) Russula tenuiceps Sebacina cf. epigaea Sebacina sp. (src657)

src2 src553 src442 – src23 R-3006 R-3010 R-4006 R-4001 R-3029 src204 R-4005 src723 src657

DQ974760 DQ974829 DQ974753 DQ974763 DQ974754 DQ974755 DQ974758 DQ974759 DQ974757 DQ974762 DQ974761 DQ974756 DQ974767 DQ974768

11 0 1 1 0 0 0 0 0 0 0 0 12 3

Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio Basidio

E E E E E E E E E E E E R R

ITS ITS ITS ITS ITS ITS ITS ITS ITS ITS ITS ITS ITS, 28s ITS, 28s

Sebacina sp. (src725) Sebacina sp. (src835)

src725 src835

DQ974769 DQ974770

0 0

Basidio Basidio

R R

ITS ITS, 28s

Humaria hemsiphaerica (AY789389) Aleuria aurantia (AY544654) Ramaria flavobrunnescens (AY102864) – – – Russula pectinatoides (AY061732) Russula emetica (AY061673) – – – – Russula parazurea (AY061704) Russula aeruginea (AF418612) – Sebacina epigaea (AF490397) Tremellodendron pallidum (AF384862) Sebacina calcea (AJ427408) Sebacina epigaea (AJ966754)

580/622 (93%) – – 534/575 (92%) 593/602 (98%) 426/464 (91%) 1050/1144 (91%) 374/405 (92%) 438/481 (91%)

183/184 (99%) 293/327 (89%) – – – 638/658 (96%) 623/644 (96%) – – – – 601/635 (94%) 549/580 (94%) – 537/585 (91%) 1111/1171 (94%) 231/244 (94%) 500/522 (95%)

Research


Taxon Name

855

856 Research


Table 1 continued


Taxon Name

Voucher1

Accession

EM cores2

Phylum3

Sporocarp type4

rDNA


% Similarity

Sebacinales1

–

DQ974764

24

Basidio

R

ITS, 28s

1109/1157 (95%)

Sebacinales2

–

DQ974766

19

Basidio

R

ITS, 28s

Sebacinales3 Sebacinales4 Sistotrema sp. 1 Tarzetta sp. (src844) Thelephora cf. anthocephala Thelephoraceae1

– – – src844 src614 –

DQ974765 EF027386 DQ974821 DQ974820 DQ974771 DQ974787

4 1 1 2 5 4

Basidio Basidio Basidio Asco Basidio Basidio

R R U/A H R R

ITS, 28s ITS ITS, 28s ITS, 28s ITS ITS

Thelephoraceae10 Thelephoraceae11

– –

DQ974797 DQ974796

1 1

Basidio Basidio

R R

ITS ITS

Thelephoraceae12 (fuscocinerea group) Thelephoraceae13 Thelephoraceae2 Thelephoraceae3 Thelephoraceae4

–

DQ974795

1

Basidio

R

ITS

– – – –

DQ974828 DQ974794 DQ974793 DQ974792

1 2 2 4

Basidio Basidio Basidio Basidio

R R R R

ITS ITS ITS ITS


– –

DQ974791 DQ974790

3 10

Basidio Basidio

R R

ITS ITS


– –

DQ974789 DQ974788

3 20

Basidio Basidio

R R

ITS ITS

Tomentella cf. atramentaria Tomentella cf. coerulea1 (src839) Tomentella cf. coerulea2 (src630) Tomentella cf. ellisii Tomentella cf. fuscocinerea Tomentella cf. lateritia Tomentella cf. nitellina Tomentella sp. (src821) Tomentella sp. (src822) Tomentella sp. (src823)

src753 src839

DQ974772 DQ974773

0 0

Basidio Basidio

R R

ITS ITS

Sebacina endomycorrhiza (AF440648) Sebacina endomycorrhiza (AF440651) Sebacina incrustans (MW 573) Sebacinales EM (AY656955) Sistotrema coronilla (AF506475) – – Thelephora caryophyllea (AJ889980) Tomentella sp. (DQ068972) Tomentella ferruginea (AF272909) Tomentella stuposa (AY010277) Tomentella lapidum (AF272941) Thelephoraceae EM (AJ893307) Tomentella sp. N44 (AJ534916) Tomentella cinerascens (AF272915) Tomentella galzinii (AJ421255) Tomentella atramentaria (AF272904) Tomentella badia (AF272937) Thelephora caryophyllea (AJ889980) – –

src630

DQ974774

0

Basidio

R

ITS

–

–

src846 src813 src833 src675 src821 src822 src823

DQ974775 DQ974776 DQ974777 DQ974778 DQ974779 DQ974780 DQ974781

0 4 0 12 0 0 0

Basidio Basidio Basidio Basidio Basidio Basidio Basidio

R R R R R R R

ITS ITS ITS ITS ITS ITS ITS

– – – – Tomentella galzinii (AJ421255) Thelephoraceae EM (AJ893317) Thelephora americana (U83487)

– – – – 521/538 (96%) 641/662 (96%) 598/648 (92%)

1133/1179 (96%) 1134/1168 (97%) 528/541 (97%) 441/461 (95%) – – 490/523 (93%) 469/508 (92%) 545/555 (98%) 471/511 (92%) 387/407 (95%) 620/635 (97%) 636/690 (92%) 528/534 (98%) 433/465 (93%) 450/470 (95%) 441/465 (94%) 584/624 (93%) – –


Table 1 continued

Taxon Name

Voucher1

Accession

EM cores2

Phylum3

Sporocarp type4

rDNA


% Similarity

Tomentella sp. (src824)

src824

DQ974782

32

Basidio

R

ITS

390/418 (93%)


src834

DQ974783

0

Basidio

R

ITS


src857

DQ974784

33

Basidio

R

ITS


src859

DQ974785

0

Basidio

R

ITS


src862

DQ974786

0

Basidio

R

ITS

Tricholoma cf. flavovirens Tricholoma portentosum Tricholoma saponeaceum Tricholoma ustale Tuber californicum Tuber cf. candidum Tuber cf. whetstonense

SF-4 src522 src582 src352 src880 src625 –

DQ974703 DQ974700 DQ974701 DQ974702 DQ974799 DQ974807 DQ974800

0 0 0 0 0 43 11

Basidio Basidio Basidio Basidio Asco Asco Asco

E E E E H H H

ITS ITS ITS ITS ITS ITS, 28s ITS

Tuber sp. (src664)

src664

DQ974808

7

Asco

H

ITS, 28s

Tuber sp. (src709)

src709

DQ974798

0

Asco

H

ITS

Tomentella fusco-cinerea (AF272942) Tomentella atramentaria (AF272904) Tomentella lapidum (AF272941) Tomentella botryoides (AF272912) Tomentella subclavigera (AY010275) – – – – – – Tuber whetstonense (AY830855) Tuber oligospermum (AY515306) Tuber quercicola (AY918957)

444/470 (94%) 515/546 (94%) 392/410 (95%) 539/585 (92%) – – – – – – 499/511 (97%) 374/388 (96%) 230/246 (93%)

1

Taxa collected as sporocarps have voucher numbers: SRC (Sierra Research Center), SF (Sierra Field), and R- (Russula). –, Taxa where no sporocarp was found. The number of soil cores where each taxon was detected (out of a total of 94 cores). Taxa on roots and as sporocarps have both a voucher and EM core number. 3 Asco, Ascomycota; Basidio, Basidiomycota. 4 All taxa have been assigned to one of four sporocarp categories: epigeous (E), hypogeous (H), resupinate (R) or unknown/asexual (U/A). 5 BLAST results are not shown when sporocarps were identified to species or genus by taxonomic experts or when BLAST results were uninformative. In cases where the best BLAST matches were unidentified environmental samples, an identified sequence with a slightly lower BLAST score is shown for clarity. 6 Putatively matched to EM roots: sporocarps of H. lacunosa were found but no ITS sequence could be obtained. 7 Lineages as identified in Douhan & Rizzo (2005). 8 Lineages as identified in Smith et al. (2006a). 9 Some root associated Ascomycota were previously considered ericoid mycorrhizas or root endophytes but have recently been shown to form ectomycorrhizas (Vrålstad, 2004 and references therein), we have included several here as putative Quercus ectomycorrhizal symbionts. 2

Research


857

858 Research Table 2 Relative importance of different ectomycorrhizal fungal groups on the roots of Quercus douglasii as indicated by the number of species, the relative frequency, and the relative abundance Fungal group

Number of species (n = 92)

Relative frequency (%)

Relative abundance (%)

Taxonomic group1 Thelephoraceae Cortinariaceae Sebacinales2 Tuberaceae Pyronemataceae Cenococcum3 Pezizaceae Other Basidiomycota4 Russulaceae Other Ascomycota5

17 16 6 3 17 3 8 10 5 5

26.3 16.4 12.0 11.6 10.9 8.2 5.5 4.6 2.9 1.7

33.0 15.9 17.1 16.1 7.0 2.5 3.0 3.6 1.5 0.4

Sporocarp type Resupinate Hypogeous Epigeous Unknown/asexual

23 24 31 14

38.4 27.5 22.6 11.5

50.1 26.8 18.8 4.3

Relative frequency is defined as the total number of occurrences of the particular taxonomic or morphological group divided by the total number of occurrences of all taxa. Relative abundance is defined as the total number of clones of the particular group divided by the total number of clones of all taxa. 1 Taxa were placed in family level groups based on a combination of data from BLAST searches, sequence alignments and the advice of taxonomic experts. 2 Sebacinales is listed at the Order level because of taxonomic uncertainties within this group (Weiss et al., 2004). 3 Cenococcum is listed at the genus level because it is an asexual ascomycete of uncertain taxonomic placement (Douhan & Rizzo, 2005). 4 Other Basidiomycota includes species from several disparate genera: Clavulina, Entoloma, Gautieria, Hygrophorus, Hysterangium, Octavianina and Sistotrema. 5 Other Ascomycota includes species within the genera Helvella and Phialophora as well as three poorly resolved taxa (Helotiales1, Helotiales2 and Ericoid ascomycete1).

this pattern was not true for every group (Table 2). Table 1 shows the frequency (number of cores) of each taxon. Reporting of abundance for individual taxa was problematic because RFLP did not resolve all species within the Thelephoraceae, the Cortinarius flexipes group or between Sebacinales2 and Sebacinales4 (one core only). Thus, 666 clones were not resolved at the species level, 605 of which were Thelephoraceae. When grouped by sporocarp type, 31 taxa had epigeous or likely epigeous sporocarps (33.7%), 24 had hypogeous or likely hypogeous sporocarps (26.1%), 23 had resupinate or likely resupinate sporocarps (25.0%) and 14 had an unknown fruiting form or reproduce asexually (15.2%) (Table 2). A large number of Ascomycota and hypogeous taxa were detected relative to other studies of EM fungi associated with hosts in the Fagaceae (see the Supplementary Material, Table S1). In addition to EM fungi, we encountered 51 clones from 20 cores that were either chimeras or nonEM taxa. Chimeras were detected in nine cores and 18 nonEM taxa were identified from 13 cores (Table S2). Fungal community of EM sporocarps Sporocarps of 108 EM species were collected over a 4-yr period; 83 were Basidiomycota and 25 were Ascomycota (Table 1). Many taxa were rare and encountered only one or


a few times. When grouped by sporocarp type, 60 species produced epigeous sporocarps (55.6%), 28 produced hypogeous sporocarps (26.4%) and 20 produced resupinate sporocarps (18.5%). Most Ascomyota formed hypogeous sporocarps, and most hypogeous sporocarps were Ascomycota; only five Basidiomycota formed hypogeous sporocarps. By contrast, most epigeous species were Basidiomycota; only three Ascomycota were epigeous. All resupinate sporocarps were Basidiomycota; we detected 15 Tomentella spp. and four Sebacinales spp. The most speciose genera based on sporocarps were Tomentella (15 spp.), Russula (11 spp.), and Inocybe (8 spp.). Spatiotemporal patterns in the EM community Approximately 92% of the common EM fungal species (those found in four or more soil cores) were detected on EM roots in both seasons and in both years (Fig. 2). One exceptional species, Inocybe sororia, was not detected in 2003 but was present in nine cores during 2004. The opposite pattern was never observed; none of the common species from 2003 were absent during 2004. Among infrequently occurring species (taxa found in 2–4 cores) much wider variation was observed between seasons and years. Thirteen species were found only in 2003 or 2004 and six species were detected only in winter


Research

Fig. 2 Frequency of most commonly detected ectomycorrhizal (EM) fungi on roots by (a) year (2003 or 2004) and (b) season (winter or spring). Data from 94 soil cores.

Table 3 Evidence of persistence in space and time of nine species of ectomycorrhizal (EM) fungi on Quercus douglasii roots 20031

Taxon name

Total number of occurrences

Number of sampling grids

Winter (n = 24)

Thelephoraceae5 Pezizaceae 1 (Scabropezia clade) Helvella cf. lacunosa Entoloma sp. (src742) Laccaria cf. bicolor Pyronemataceae1 Genea gardneri3 Tuber sp. (src664)3 Genea cf. harknessii (src616)3

3 3 2 2 2 4 4 6 7

1 1 1 1 1 2 2 3 4

1 12

7 13 10

20041 Spring (n = 20)

Winter (n = 26)

Spring (n = 24)

Persistence between years2

1 12 3 4

1 12 3 4

13, 10 10 3, 10 6, 11

10 10 3 11

+ + – – – + + + +

7 10 3 6, 11

These taxa were found three to seven times over a 2-yr period but were restricted to small, 8 × 8 m plots. Numbers indicate 8 × 8 m sampling squares where the species were detected; see Fig. 1. 2 +, Species that persisted between years; –, species found only in 1 yr. 3 These taxa also found in one or two other cores on a single sampling date (these data have been omitted for clarity). 1

or spring (Fig. 2). There were no obvious groupings based on taxonomy or sporocarp type. The stability of species composition over time suggests persistence of EM taxa in small habitat patches, a pattern that was notable for several infrequently occurring species. Table 3

depicts nine species found on multiple sampling dates, but only in one or a few limited areas of the plot. Six of these species persisted in the same location(s) during both years, while the other three were found during two consecutive seasons of the same year.



859

860 Research

Fig. 3 Dendrogram depicting hierarchical cluster analysis of soil cores grouped by similarity of species composition of ectomycorrhizal fungi, generated using the Sørrenson distance measure and flexible Beta linkage method. Data are from 59 cores sampled at the largest spatial scale (8 × 8 m plots). For each number in the dendrogram, the first number designates the location within the sampling grid and the number preceded by a ‘T’ designates the sampling time (1–4). Black rectangles indicate where cores from within the same 25-cm radius, but at sampled at different times, had similar or identical species composition as indicated by adjacent positions in the dendrogram. Grey rectangles indicate cores that were close in terms of species composition, but were not directly adjacent in the dendrogram.

Hierarchical cluster analysis Regardless of whether rare species were considered equally or not, cluster analyses produced dendrograms with similar patterns. Only the dendrogram created with the Sørrenson distance measure and flexible Beta linkage method is shown (Fig. 3). Soil cores from within the same 25-cm radius tended to have similar species composition, regardless of sampling date. Consistent with this pattern, MRPP analysis indicated that species composition was correlated with core position (Euclidean P = 0.001, Sørrenson P = 0.001). By contrast, MRPP constrained by season (all samples from a season are pooled, regardless of year) or sampling date (each year × season combination treated separately) were not significantly different from random association (season: Euclidean P = 0.707, Sørrenson P = 0.515; sampling date: Euclidean P = 0.464, Sørrenson


P = 0.475). However, when MRPP was constrained by year (all samples pooled by year regardless of season) the two distance measures produced conflicting results (Euclidean P = 0.048, Sørrenson P = 0.150), suggesting a weak effect due to year.

Discussion The EM community associated with Quercus douglasii Although the diversity of EM fungi in this system was comparable to that of more mesic habitats (Avis et al., 2003; Walker et al., 2005), the Mediterranean climate of California’s Sierra Nevada foothills strongly influenced the phylogenetic affinities and sporocarp forms of the EM fungi associated with Q. douglasii. Two of the most striking features of this community were that many of the EM species on roots


Research

(c. 40%) were Ascomycota and that many (c. 26%) produced hypogeous sporocarps. We found species of Pezizales in > 80% of soil cores and encountered several pezizalean genera that have rarely or never been detected on roots (e.g. Genabea, Hydnobolites, Hydnotryopsis, Pachyphloeus and Tarzetta) (Tedersoo et al., 2006). Despite high ascomycete diversity, the EM community on roots was dominated overall by basidiomycetes, particularly species with inconspicuous sporocarps (e.g. Inocybe, Thelephoraceae). Except for Inocybe, agaricoid taxa were abundant and diverse as sporocarps but either infrequent or altogether absent on roots (e.g. Amanita, Boletus and Hygrophorus). Even Russula and Lactarius, two dominant genera on EM roots worldwide (Horton & Bruns, 2001), were poorly represented on Q. douglasii. Although the reasons for this phenomenon are unclear, preliminary data suggest that agaricoid taxa are poorly adapted for the dry microclimate beneath deciduous Q. douglasii and instead associate preferentially with two cooccurring evergreen hosts that have thicker litter (Q. wislizeni and P. sabiniana) (M. Smith and M. Morris, unpublished). Previous studies have reported a pervasive disconnect between EM fungi detected on roots and those detected as sporocarps (Gardes & Bruns, 1996; van der Heijden et al., 1999; Horton & Bruns, 2001; Yamada & Katsuya, 2001; Richard et al., 2005). Few studies, however, have intensively sampled both EM roots and sporocarps simultaneously; in particular, most studies have overlooked inconspicuous sporocarps (Izzo et al., 2005b). We understand that some taxa produce abundant sporocarps and others produce them sparingly (Luoma et al., 1991; Gardes & Bruns, 1996) but the reported disconnect between EM roots and sporocarps has often been interpreted to mean that species with ‘missing sporocarps’ are present only on EM roots or in soil (van der Heijden et al., 1999; Yamada & Katsuya, 2001). There are certainly cases where true discrepancies exist between ‘aboveground’ and ‘below-ground’ biomass; we detected at least one species that produced extremely abundant sporocarps, but was rare on EM roots (Hygrophorus eburneus; M. Smith, pers. obs.). However, in their seminal work Gardes & Bruns (1996) noted that many taxa found on roots might produce inconspicuous sporocarps below-ground or within leaf litter. In this study, we matched EM roots to sporocarps for c. 45% of the EM root community, despite the cryptic and ephemeral nature of most fungi in this arid woodland. Had we surveyed only epigeous sporocarps, the level of correspondence between fruiting bodies and roots would have been quite low (c. 20%). It is possible that unmatched taxa on roots did not produce sporocarps or produced them rarely (Gardes & Bruns, 1996; Straatsma et al., 2001), but an alternative explanation is that their sporocarps were present but difficult to detect and collect. Moreover, many ‘unmatched’ taxa were from hyperdiverse taxonomic groups (e.g. Thelephoraceae, Inocybe) or cryptic species complexes (e.g. Genea harknesii, Cortinarius flexipes group) where species are difficult or impossible to discriminate morphologically. Given these results we suggest that the

purported disconnect between EM roots and sporocarps may be heavily influenced by sampling difficulties and the inherent visual bias toward epigeous species. It is possible that the inclusion of hypogeous and resupinate taxa in future sporocarp surveys will decrease this perceived disconnect. Although the Mediterranean ecosystem and unique mycorrhizal host certainly influenced the composition of EM fungi detected in this study, it is also possible that the pooled root technique affected our perception of the community. For example, members of the C. geophilum species complex were uncommon relative to many other EM studies (Avis et al., 2003); a finding that may result from the small size of C. geophilum EM tips, which may yield less rDNA template, less PCR product and fewer clones than species with larger tips. Despite its potential pitfalls, we used a pooled root technique to avoid the low PCR success rate previously reported from individual EM roots of Fagaceae (e.g. 49–70%)(Avis et al., 2003; Walker et al., 2005). We reasoned that our pooling technique should detect the most common EM taxa but eliminate the need for morphotyping and the problem of poor DNA sequencing due to contamination or dual EM colonization. Spatiotemporal patterns of EM fungi on Quercus douglasii roots Because of the strong seasonality of the Mediterranean climate in the Sierra Nevada foothills, we expected significant seasonal variation within the EM community on roots and greater species-turnover between years. Although winter and spring sampling times were separated by only 3 months, these seasons represented biological extremes. When winter soil cores were collected, temperatures were cool, soils were saturated, and shoots and leaves were dormant. When spring cores were collected, temperatures were hot, soils were much drier, and shoots and leaves were growing (see the Supplementary Material, Fig. S1). Despite these considerable seasonal differences, several lines of evidence suggest that the EM fungal community structure on roots of Q. douglasii remained relatively stable over a 2-yr period. First, the same taxa were dominant on all four sampling dates and only one common species, I. sororia, varied substantially between years (Fig. 2). Second, hierarchical cluster analysis indicated that cores collected within a 25-cm radius contained similar species regardless of the sampling date, suggesting that most taxa either persisted in or recolonized small areas (Fig. 3). Third, persistence was further highlighted by several rare taxa that were detected only on roots within the same 25 cm radius but on successive sampling dates (Table 3). The hypothesis that the EM community would be temporally dynamic was influenced by results from two recent studies. Izzo et al. (2005a) examined the EM root community in a montane, Abies-dominated California forest and found that although species composition remained relatively stable at large spatial scales over 3 yr, common species went locally



861

862 Research

extinct within many plots. In a second study, Buée et al. (2005) intensively sampled morphotypes of Fagus in France over 1 yr and found seasonal shifts in the EM community, particularly for two dominant species. In contrast to the results of Izzo et al. (2005a), none of the common species in our plot were ephemeral and we never detected ‘local extinction’. Similarly, we did not detect the wide seasonal variation seen by Buée et al. (2005). Although sampling design probably accounts for some of the differences between our study and those of Izzo et al. (2005a) and Buée et al. (2005), differences in other factors, such as EM host species, may have affected the dynamics of the EM root communities. Two recent studies by Querejeta et al. (2003, 2006) highlight how deep-rooted, drought-adapted hosts, such as Q. douglasii, might enhance the survival and persistence of EM fungi during prolonged drought. Querejeta et al. (2003, 2007) showed that Quercus pull water from deep soil layers by hydraulic lift and translocate the water directly to EM hyphae in the rhizosphere. Consequently, despite xeric conditions above-ground, EM fungi connected to Q. douglasii roots may be buffered from the effects of drought. Sustaining EM fungi in this way presumably benefits Quercus hosts because viable symbionts can remain active longer into the summer and grow quickly at the onset of winter rains. While we did not find evidence for frequent species turnover or strong seasonal changes in the EM root community, our cluster analysis indicated that species composition was more similar within years than between years. This suggests that the EM community overall changed via ‘slow drift’ rather than ‘rapid turnover’, at least at the site and spatiotemporal scale we examined. However, despite the large number of soil cores and rDNA sequences we examined, it seems likely that sampling with increased frequency or over more years would have enhanced our view of EM temporal dynamics. When our results are viewed with the studies mentioned above, they highlight how environmental conditions, host species, experimental design, or other unknown factors might affect our understanding of temporal dynamics within EM communities. More research from diverse EM systems is needed before we can make general and definitive statements about temporal dynamics within EM fungal below-ground communities.

Acknowledgements This work was made possible by valuable contributions of many taxonomic experts including J. Trappe, K. Hansen, B. Perry, J. Ammirati, B. Matheny, U. Köljalg, S. Dunham, M. Weiss, J. Frank, S. Miller, and R. M. Davis. We thank M. Morris, K. Huryn, R. Albright, E. Smith and F. Smith for laboratory and field assistance. UCSFREC staff provided logistical support for fieldwork, A. Fremier supplied invaluable assistance during data analysis and R. Vilgalys suggested candidate 28s primers. Helpful comments on previous drafts


were provided by A. Izzo, C. Bledsoe, T. Gordon and several anonymous reviewers. Research was supported by grants to D. Rizzo by the National Science Foundation (DEB-99-81711) and to M. Smith by the Mycological Society of America and the North American Truffling Society.

References Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research 25: 3389–3402. Anderson IC, Campbell CD, Prosser JI. 2003. Potential bias of fungal 18s rDNA and internal transcribed spacer polymerase chain reaction primers for estimating fungal biodiversity in soil. Environmental Microbiology 5: 36–47. Avis PG, McLaughlin DJ, Dentinger BC, Reich PB. 2003. Long-term increase in nitrogen supply alters above- and below-ground ectomycorrhizal communities and increases dominance of Russula spp. in a temperate oak savanna. New Phytologist 160: 239–253. Bruns TB. 1995. Thoughts on the processes that maintain local species diversity of ectomycorrhizal fungi. In: Collins HP, Gilbertson GP, Klug MJ, eds. The significance and regulation of soil biodiversity. Dordrecht, the Netherlands: Kluwer Academic Publishers, 47–62. Buée M, Vairelles D, Garbaye J. 2005. Year-round monitoring of diversity and potential metabolic activity of the ectomycorrhizal community in a beech (Fagus silvatica) forest subjected to two thinning regimes. Mycorrhiza 15: 235–245. Chenna R, Sugawara H, Koide T, Lopez R, Gibson TJ, Higgins DG, Thompson JD. 2003. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Research 31: 3497–3500. Cline ET, Ammirati JE, Edmonds RL. 2005. Does proximity to mature trees influence ectomycorrhizal fungus communities of Douglas-fir seedlings? New Phytologist 166: 993–1009. Douhan GW, Rizzo DM. 2005. Phylogenetic divergence in a local population of the ectomycorrhizal fungus Cenococcum geophilum. New Phytologist 166: 263–271. Gardes M, Bruns TD. 1993. ITS primers with enhanced specificity for basidiomycetes – application to the identification of mycorrhizae and rusts. Molecular Ecology 2: 113–118. Gardes M, Bruns TD. 1996. Community structure of ectomycorrhizal fungi in a Pinus muricata forest: above-and below-ground views. Canadian Journal of Botany 74: 1572–1583. Gehring CA, Theimer TC, Whitham TG, Keim P. 1998. Ectomycorrhizal fungal community structure of pinyon pines growing in two environmental extremes. Ecology 79: 1562–1572. Hall TA. 1999. BIOEDIT: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series 41: 95–98. van der Heijden EW, de Vries FW, Kuyper TW. 1999. Mycorrhizal associations of Salix repens L. communities in succession of dune ecosystems. I. Above-ground and below-ground views of ectomycorrhizal fungi in relation to soil chemistry. Canadian Journal of Botany 77: 1821– 1832. Hopple J, Vilgalys R. 1994. Phylogenetic relationships among coprinoid taxa and allies based on data from restriction site mapping of nuclear rDNA. Mycologia 86: 96–107. Horton TR. 2002. Molecular approaches to ectomycorrhizal diversity studies: variation in ITS at a local scale. Plant and Soil 244: 29–39. Horton TR, Bruns TD. 2001. The molecular revolution in ectomycorrhizal ecology: peeking into the black box. Molecular Ecology 10: 1855–1871. Izzo AD, Agbowo J, Bruns TD. 2005a. Detection of plot-level changes in ectomycorrhizal communities across years in an old-growth mixed-conifer forest. New Phytologist 166: 619–630.


Research Izzo AD, Meyer M, Trappe JM, North M, Bruns TD. 2005b. Hypogeous ectomycorrhizal fungal species on roots and in small mammal diet in a mixed-conifer forest. Forest Science 51: 243–254. Kennedy PG, Izzo AD, Bruns TD. 2003. There is high potential for the formation of common mycorrhizal networks between understory and canopy trees in a mixed evergreen forest. Journal of Ecology 91: 1071–1080. Köljalg U, Dahlber A, Taylor AFS, Larsson E, Hallenberg N, Stenlid J, Larsson K-H, Fransson PM, Kåren O, Jonsson L. 2000. Diversity and abundance of resupinate thelephoroid fungi as ectomycorrhizal symbionts in Swedish boreal forests. Molecular Ecology 9: 1985–1996. Landeweert R, Leeflang P, Kuyper TW, Hoffland E, Rosling A, Wernars K, Smit E. 2003. Molecular identification of ectomycorrhizal mycelium in soil horizons. Applied and Environmental Microbiology 69: 327–333. Luoma DL, Frenkel RE, Trappe JM. 1991. Fruiting of hypogeous fungi in Oregon Douglas-fir forests: seasonal and habitat variation. Mycologia 83: 335–353. Martin D, Rybicki E. 2000. RDP: detection of recombination amongst aligned sequences. Bioinformatics 16: 562–563. McClaran MP. 1988. Comparison of fire history estimates between open-scarred and intact Quercus douglasii. American Midlands Naturalist 120: 432–435. McCune B, Mefford MJ. 1999. PC-ORD: Multivariate Analysis of Ecological Data v. 4.20. Glenden Beach, OR, USA; MjM Software Design. Milikin CS, Bledsoe CS. 1999. Biomass and distribution of fine and coarse roots from blue oak (Quercus douglasii) trees in northern Sierra Nevada foothills of California. Plant and Soil 214: 27–38. Nara K, Nakaya H, Wu B, Zhou Z, Hogetsu T. 2003. Underground primary succession of ectomycorrhizal fungi in a volcanic desert on Mount Fuji. New Phytologist 159: 743–756. Nishida FH. 1989. Key to the species of Inocybe in California. Mycotaxon 34: 181–196. O’Brien HE, Parrent JL, Jackson JA, Moncalvo J-M, Vilgalys R. 2005. Fungal community analysis by large-scale sequencing of environmental samples. Applied and Environmental Microbiology 71: 5544–5550. Parks HE. 1921. California hypogeous fungi – Tuberaceae. Mycologia 13: 301–314. Pavlik BM, Muick PC, Johnson SG, Popper M. 1991. Oaks of California. Los Olivos, CA, USA: Cachuma Press. Qiu X, Wu L, Huang H, McDonel PE, Palumbo AV, Tiedje JM, Zhou J. 2001. Evaluation of PCR-generated chimeras, mutations, and hetereoduplexes with 16s rRNA gene-based cloning. Applied and Environmental Microbiology 67: 880–887. Querejeta JI, Egerton-Warburton L, Allen MF. 2003. Direct nocturnal water transfer from oaks to their mycorrhizal symbionts during severe soil drying. Oecologia 134: 55–64. Querejeta JI, Egerton-Warburton L, Allen MF. 2007. Hydraulic lift may buffer rhizosphere hyphae against the negative effects of severe soil drying in a California Oak savanna. Soil Biology and Biochemistry 39: 409–417. Richard F, Moreau P-A, Selosse M-A, Gardes M. 2004. Diversity and fruiting patterns of ectomycorrhizal and litter saprobic fungi in an old-growth Mediterranean forest dominated by Quercus ilex L. Canadian Journal of Botany 82: 1711–1729. Richard F, Millot S, Gardes M, Selosse MA. 2005. Diversity and specificity of ectomycorrhizal fungi retrieved from an old-growth Mediterranean forest dominated by Quercus ilex. New Phytologist 166: 1011–1023. Schadt CW, Martin AP, Lipson DA, Schmidt SK. 2003. Seasonal dynamics of previously unknown fungal lineages in tundra soils. Science 301: 1359– 1361. Smith ME, Trappe JM, Rizzo DM. 2006a. Genea, Genabea, and Gilkeya General. nov. ascomata and ectomycorrhiza formation in a Quercus woodland. Mycologia 98: 699–716. Smith ME, Trappe JM, Rizzo DM, Miller SM. 2006b. Gymnomyces xerophilus sp. nov. (sequestrate Russulaceae), an ectomycorrhizal associate of Quercus in California. Mycological Research 110: 572–582.

Smith SE, Read DJ. 1997. Mycorrhizal symbiosis, 2nd edn. London, UK: Academic Press. Standiford RB. 2002. California’s Oak woodlands. In: McShea WJ, Healy WM, eds. Oak forest ecosystems: ecology and management for wildlife. Baltimore, MD, USA: Johns Hopkins University Press, 280–303. Straatsma G, Ayer F, Egli S. 2001. Species richness, abundance, and phenology of fungal fruit bodies over 21 years in a Swiss forest plot. Mycological Research 105: 515–523. Swofford DL. 2001. PAUP*: phylogenetic analysis using parsimony (*and other methods), version 4.0b10. Sunderland, MA, USA: Sinauer Associates. Taylor DL, Bruns TD. 1999. Community structure of ectomycorrhizal fungi in a Pinus muricata forest: minimal overlap between the mature forest and resistant propagule communites. Molecular Ecology 8: 1837–1850. Tedersoo L, Köljalg U, Hallenberg N, Larsson K. 2003. Fine scale distribution of ectomycorrhizal fungi and roots across substrate layers including course woody debris in a mixed forest. New Phytologist 159: 153–165. Tedersoo L, Hansen K, Perry BA, Kjøller R. 2006. Molecular and morphological diversity of pezizalean ectomycorrhiza. New Phytologist 170: 581–596. Vrålstad T. 2004. Are ericoid and ectomycorrhizal fungi part of a common guild? New Phytologist 164: 7–10. Walker JF, Miller OK, Horton JL. 2005. Hyperdiversity of ectomycorrhizal fungus assemblages on oak seedlings in mixed forests in the southern Appalachian Mountains. Molecular Ecology 14: 829–838. Weiss M, Selosse MA, Rexer KH, Urban A, Overwinkler F. 2004. Sebacinales: a hitherto overlooked cosm of heterobasidiomycetes with a broad mycorrhizal potential. Mycological Research 108: 1003–1010. White TJ, Bruns TD, 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, CA, USA: Academic Press. Yamada A, Katsuya K. 2001. The disparity between the number of ectomycorrhizal fungi and those producing fruiting bodies in a Pinus densiflora stand. Mycological Research 105: 957–965.

Supplementary Material The following supplementary material is available for this article online: Fig. S1 Four sampling dates for ectomycorrhizal roots of Quercus douglasii relative to temperature and precipitation during 2003 and 2004. Sampling dates are designated T1–T4. Table S1 Comparison of ectomycorrhizal fungal communities associated with Fagaceae, as detected on roots and as sporocarps, from five studies that used molecular methods. Table S2 Occurrence of nonectomycorrhizal fungi and fungal-like organisms associated with ectomycorrhizal roots of Quercus douglasii. This material is available as part of the online article from: http://www.blackwell-synergy.com/doi/abs/10.1111/ j.1469-8137.2007.02040.x (This link will take you to the article abstract). Please note: Blackwell Publishing are not responsible for the content or functionality of any supplementary materials supplied by the authors. Any queries (other than missing material) should be directed to New Phytologist Central Office.



863