Mycologia
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

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

Thysanophora penicillioides includes multiple genetically diverged groups that coexist respectively in Abies mariesii forests in Japan


Susumu Iwamoto 1
Seiji Tokumasu

     Sugadaira Montane Research Center, University of Tsukuba, Sugadaira, Sanada, Nagano 386-2201, Japan

Yoshihisa Suyama

     Graduate School of Agricultural Science, Tohoku University, Kawatabi, Naruko, Miyagi 989-6711, Japan

Makoto Kakishima

     Institute of Agriculture and Forestry, University of Tsukuba, Tsukuba, Ibaraki 305-8572, Japan

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 

We investigated intraspecific diversity and genetic structures of a saprotrophic fungus—Thysanophora penicillioides—based on sequences of nuclear ribosomal internal transcribed spacer (ITS) in 15 discontinuous Abies mariesii forests of Japan. In such a well-defined morphological species, numerous unexpected ITS variations were revealed: 12 ITS sequence types detected in 254 isolates collected from 15 local populations were classified into five ITS sequence groups. Maximally, four ITS groups consisted of seven ITS types coexisting in one population. However, group 1 was dominant with approximately 65%; in particular, one haplotype, 1a, was most dominant with approximately 60% in respective populations. Therefore, few differences were recognized in genetic structure among local populations, implying that the gene flow of each lineage of the fungus occurs among local populations without geographic limitations. However, minor haplotypes in some ITS groups were found only in restricted areas, suggesting that they might expand steadily from their places of origin to neighboring A. mariesii forests. Aggregating sequence data of seven European strains and four North American strains from various substrates to those of Japanese strains, 18 ITS sequence types and 28 variable sites were recognized. They were clustered into nine lineages by phylogenetic analyses of the ß-tubulin and combined ITS and ß-tubulin datasets. According to phylogenetic species recognition by the concordance of genealogies, respective lineages correspond to phylogenetic species. Plural phylogenetic species coexist in a local population in an A. mariesii forest in Japan.

Key words: Biogeography, ß-tubulin, Cryptic species, Fungal population, Gene flow, rDNA ITS


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Analyses of population genetic structures based on DNA polymorphisms have been applied for various molecular ecological studies of fungi. Those analyses have examined gene flow, genetic differentiation and isolation, phylogeography, and dispersion ability (Milgroom 1996Go, Rogers and Rogers 1999Go, Taylor et al. 2000Go). Recently, studies of intraspecific diversity of many economically important fungal taxa have been undertaken energetically to examine mushrooms (e.g. Xu et al 1997Go, Hughes et al 1998Go, Coetzee et al 2000Go, James and Vilgalys 2001Go), plant pathogenic fungi (e.g. Boeger et al 1993Go, Trigiano et al 1995Go, Milgroom et al 1996Go, Koenig et al 1997Go, Carbone and Kohn 2001Go), human pathogenic fungi (e.g. Gräser et al 1996Go, Xu et al 2000Go), and entomopathogenic fungi (e.g. Boucias et al 2000Go, Jensen and Eilenberg 2001Go, Fargues et al 2002Go). Less information exists regarding saprotrophic fungi that are not related directly to economic activity. Notwithstanding, recent environmental changes derived from human activity have required such information related to saprotrophic fungi because their distributions are not influenced greatly by artificial effects.

Some obstacles remain for biogeography and population studies of fungi. One is that collecting large samples of fungi from a wide area is difficult. In particular, habitats and distribution patterns of most saprotrophic fungi have not been observed clearly. Therefore, efficient sample collection is important for the study of population genetic structures (James and Vilgalys 2001Go). Another obstacle is the difficulty in defining fungal populations because direct confirmation of their individuals is lacking in the field (Rogers and Rogers 1999Go). Still, the distribution of a given fungus is probably discontinuous because the host or suitable substrate for the fungus is distributed discontinuously in most cases. Therefore, it is reasonable that a litter fungal species inhabiting leaf litter under a given host community is defined as a local population of the fungal species.

We selected a saprotrophic microfungus, Thysanophora penicillioides (Roum.) W. B. Kendr. as a target fungal species for fungal population study. Advantages of using the fungus include the following. 1) It is possible to collect many isolates because ecology of the fungus is well known. It appears in a specific habitat—decaying fir needles. 2) Its population is definable. One host, Abies mariesii Mast. (o-shirabiso fir) is dominant in discontinuous subalpine forests throughout Japan. For that reason, we can readily recognize species assemblages in respective forests. 3) It is a morphologically well-defined species. Therefore, we can isolate the fungus from substrates without confusion with other species in the same niche.

We conducted a study to reveal intraspecific diversity and population genetic structures of Thysanophora penicillioides in Abies mariesii forests in the upland mountains interspersed with central to eastern parts of Japan. This study revealed that the fungus had unexpectedly great genetic variation and probably contains plural cryptic species. This study elucidates the population genetic structures of A. mariesii forests and the genetic diversity of the fungus through analyses of the ITS sequences and the ß-tubulin gene sequences based on Japanese isolates and several European and North American isolates from public culture collections.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Target fungus.— – As a target fungus, we selected a dematiaceous hyphomycete: Thysanophora penicillioides. It produces minute dry phialoconidia and no ascomata have been found. It mainly inhabits decaying needles of pinaceous trees, especially Abies and Picea (Kendrick 1961Go, Sieber-canavesi and Sieber 1993Go, van Maanen and Gourbière 1997Go). Its natural distribution roughly overlaps with that of Pinaceae in the northern hemisphere. We infer that the fungus is commonly distributed over wide areas of Japan because investigators have reported that the fungus was isolated from fallen conifer needles and soils in several regions of Japan (Kominami et al 1952Go, Aoki et al 1990Go, Tokumasu 1996Go, Iwamoto and Tokumasu 2001Go), and fir species that appear to be main hosts of the fungus, distribute all over Japan except in subtropical regions such as Okinawa’s islands.

Fungal isolation.— – Forests of Abies mariesii were selected as sampling sites for this study. The fir is distributed only in subalpine zones in high mountains located in central and eastern Honshu (Hayashi 1960Go). The distributions in the forests are discontinuous (FIG. 1Go).


Figure 1
View larger version (24K):
[in this window]
[in a new window]
 
FIG. 1. Natural distribution area of Abies mariesii and local populations of Thysanophora penicillioides in Japan. Gray shadows show the distribution area of A. mariesii after Hayashi (1960)Go. Black circles show local populations of T. penicillioides under A. mariesii forests in the study.

 
Fifteen forests of A. mariesii (Hakkoda = HD, Hachimantai = HT, Kayachine = HY, Kurikoma = KK, Zao = ZO, Higashi-Azuma = HA, Nishi-Azuma = NA, Myoko = MY, Naeba = NB, Tateyama = TY, Nishi-Hotaka = NH, Ontake = OT, Yatsugatake = YG, Asashi = AS, Senjyo = SJ) were selected from the entire distributional area of A. mariesii. Individuals of T. penicillioides inhabiting respective forests were inferred to be respective local populations of the species (TABLE IGo and FIG. 1Go).


View this table:
[in this window]
[in a new window]
 
TABLE I. Local populations of Thysanophora penicillioides in Abies mariesii forests

 
Decaying fir needles were collected at 12–35 sites (average 21 sites) in each forest during Aug–Oct 2000 (TABLE IGo). Distances between the sites within populations were at least 20 m. Samples for each site were collected by cutting a 10 x 10 cm section of litter from the forest floor. The samples were then brought back to the laboratory in a sterile paper bag.

Ten needles per site were used for fungal detection. Brown or partly darkened brown needles were selected from each litter sample. A modified Tokumasu’s washing method (Tokumasu 1980Go) was used for washing soil particles and other microorganisms out of the needles: ten needles placed in each test tube were washed with detergent water (0.005% aerosol-OT) and sterilized water for 3 min using an automatic mixer (S-100; Taitec Co.). After washing, the needles were dried for 12 h on sterilized filter paper on a clean bench. Subsequently, each of the needles was incubated at 20 C on a corn meal agar (Nissui Pharmaceutical Co. Ltd.) plate. After incubation, they were observed photomicroscopically for a month. One isolate per plate was reserved as the representative site when T. penicillioides appeared on plural needles in a plate. Isolates were incubated at 20 C on corn meal agar plates (Nissui Pharmaceutical Co. Ltd. or Sigma-Aldrich Corp.) for DNA techniques.

European and North American isolates.— – Seven European and four North American strains were obtained from public culture collections from Centraalbureau voor Schimmelcultures (CBS), Institute for fermentation, Osaka (IFO), and the University of Alberta Microfungus Collection & Herbarium (UAMH). Their strains were isolated in eight countries over two continents and from various substrates (TABLE IIGo).


View this table:
[in this window]
[in a new window]
 
TABLE II. Strains originated from Europe and North America

 
PCR amplification.— – DNA fragments of ITS1, 5.8 S and ITS2 rDNA with 3' end of 18 S and 5' end of 28 S rDNA partial sequences (about 600 bp) were amplified with newly designed primers—ITP1 (5'-GAGCCGGAAAGTTGGTCAAA-3') and ITP2 (5'-CGCTTCACTCGCCGTTACTG-3')—with a HotStarTaq Master Mix (Qiagen Inc.) by the direct PCR method without DNA extraction (Suyama et al 1996Go, Iwamoto et al 2002Go) under a GeneAmp PCR System 2400 (Applied Biosystems). Two-step PCR method was used with these conditions: initial denaturation at 95 C for 15 min, 40 cycles of denaturation at 94 C for 20 s, and annealing at 60 C for 1 min, a final extension period at 72 C for 10 min, and a 4 C soak. A QiAquick PCR Purification Kit (Qiagen Inc.) purified the PCR products.

The ß-tubulin gene partial fragments including three non-coding regions (about 360 bp) were amplified with benA1 and benA2 designed for Aspergillus species. (Geiser et al 1998Go). The amplification method was the same as that for ITS except that the annealing temperature was 58 C.

Sequencing.— – Sequencing reactions were performed with a BigDyeTM Terminator Cycle Sequencing FS Ready Reaction Kit (Applied Biosystems) according to the manufacturer’s instructions under a GeneAmp PCR System 2400 (Applied Biosystems). The primer ITP1 was used for reading upper sequences of ITS in all isolates. Representative isolates for phylogenetic analyses and isolates with minor changes were determined in both directions using the ITP2 primer and the mutation sites were confirmed. For isolates used in phylogenetic analyses, the ß-tubulin sequence was determined in both directions using benA1 and benA2 primers. Sequencing reaction products were purified with a DyeEx Spin Kit (Qiagen Inc.). Direct sequencing was performed using a sequencing system (ABI PRISM 377 DNA sequencing system; Applied Biosystems).

Data analyses.— – Each distinctive sequence (allele) was defined as each ITS type. Among ITS types, sequences with variation of less than two sites were classified into the same group (ITS group). Nei’s gene diversity formulae (Hep) were used to evaluate the genetic diversity of populations (Nei 1973Go). The uncorrected distance (p-distance) between sequences of ITS or ß-tubulin were calculated using PAUP version 4.0b10 (Swofford 2002Go).

Phylogenetic analysis.— – Three data sets—ITS, ß-tubulin, and ITS and ß-tubulin—were prepared of combined data including all ITS types of T. penicillioides (TABLE IIGo, IIIGo). Alignments of data sets were determined using Clustal X version 1.81 (Thompson et al 1997Go) and were refined visually.


View this table:
[in this window]
[in a new window]
 
TABLE III. Strains used in phylogenetic analyses and accession numbers of their DNA sequences

 
Phylogenetic trees were constructed by the maximum parsimony (MP) method using PAUP. MP analyses were performed with the heuristic search option and 1000 replications of bootstrap analysis (Felsenstein 1985Go). Consistency index (CI), retention index (RI), rescaled consistency index (RC), and homoplasy index (HI) were measured using PAUP. A partition homogeneity test was performed using PAUP to investigate the statistical congruence between the ITS and ß-tubulin data sets.

Sequences used for phylogenetic analyses were deposited in DDBJ. Accession numbers for ITS are AB175229 [GenBank] –AB175252 [GenBank] , AB213264 [GenBank] –AB213271 [GenBank] and numbers for ß-tubulin are AB175253 [GenBank] –AB175276 [GenBank] , AB213272 [GenBank] –AB213279 [GenBank] (TABLE IIGo, IIIGo).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Appearance of T. penicillioides from decaying A. mariesii needles.— Thysanophora penicillioides appeared in all forests with high frequencies of 55.6–100% (TABLE IGo), thereby confirming that it abundantly inhabited decaying A. mariesii needles in Japan. Ultimately, 254 isolates were collected from 324 examined sites of 15 local populations (TABLE IGo).

Genetic diversity and genetic structures based on ITS in Abies mariesii forests in Japan.— ITS of all Japanese isolates were determined and 21 variable sites were recognized in the aligned sequences of ITS regions (FIG. 2Go). Finally, 12 distinctive sequences (ITS types 1a–5a) were detected. These ITS types were classified into five groups in which groups 1 and 2 were included more than two ITS types (FIG. 2Go). Frequencies of respective ITS types and diversity indexes (Hep) within local populations are shown in TABLE IVGo. Overall, type 1a was dominant (59.8%), followed by type 2a (19.3%), whereas other ITS types had low frequencies (less than 6.7%). Among ITS groups, group 1 was dominant (65.4%), followed by group 2 (25.2%).


Figure 2
View larger version (24K):
[in this window]
[in a new window]
 
FIG. 2. Variable sites between sequence types in the ITS regions of Thysanophora penicillioides. ITS type 1a is used as the reference sequence. In the variable sites, the same nucleotides with the reference sequence are noted by dots (·). Gap sites are denoted by hyphens (-). The numbers above sequences respectively indicate the site numbers from the 5' terminal of ITS1 and ITS2.

 

View this table:
[in this window]
[in a new window]
 
TABLE IV. Allele frequency and gene diversity (Hep) within each local population of Thysanophora penicillioides

 
In a single population, 2–7 ITS types and 1–4 ITS groups were observed. The HY population had the highest diversity (Hep = 0.704): seven ITS types belonging to four ITS groups were detected, whereas MY and NH populations had the lowest diversity (Hep = 0.272) because they were almost entirely occupied by one ITS type (1a). The SJ population was the only local population that consisted of an ITS group.

Type 1a was the only ITS type that was detected in all local populations; its allele frequency exceeded 0.500 in ten local populations. Type 2a occurred at relatively high frequency in three northern and three southern local populations. Remarkably, the three northern populations (HD, HT and HY) had minor haplotypes belonging to group 2 (2c–2e); they were found only in this region. Minor ITS types of group 1 (1b–1d) were restricted to western A. mariesii forests (MY, NB, TY, NH, YG, AS and SJ). Type 3a was detected from two adjacent local populations (KK and ZO) and type 5a in two geographically separate northern populations (HY and HA). A member of minor ITS group Type 4a occurred in a wide area.

Genetic diversity in European and North American isolates based on ITS.— – When European and North American isolates were added, 28 variable sites were found in ITS of the fungus (FIG. 2Go). Four ITS types and two ITS groups were detected from 11 isolates collected in Europe and North America. It is noteworthy that none had been detected in Japan; new ITS types 6a, 6b, 6c and 7a were classified into new ITS groups 6 and 7 (TABLE IIGo). Type 6a was dominant in the examined isolates and detected from both Europe (four countries) and North America (two countries). Types 6b and 6c were detected respectively in one European isolate, whereas type 7a was found in three North American isolates (TABLE IIGo).

Genetic distances and Phylogenetic relationships.— – TABLE VGo shows genetic distances between ITS types. The distances within each ITS group (for groups 1, 2 and 6) are 0.002–0.003. On the other hand, the genetic distance between ITS groups was an average of 0.013 for values of 0.003–0.029. The distance of group 3 vs. group 7 was not very distant (0.003) because the variation between groups was mainly contributed by some indels (FIG. 2Go).


View this table:
[in this window]
[in a new window]
 
TABLE V. Uncorrected distances between alleles in the ITS regions of Thysanophora penicillioides

 
We also determined the ß-tubulin sequences of representative isolates selected based on the ITS type (TABLES IIGo, IIIGo). The analyzed ß-tubulin region was more variable than ITS and many alleles were detected (TABLES IIGo, IIIGo). The genetic distance matrix between each ß-tubulin allele is shown in TABLE VIGo. Considered based on the classification of its groups, grouping for groups 1–5 of ß-tubulin well corresponded to that of ITS; the distances within their groups and among groups were 0–0.028 and 0.031–0.082 respectively (TABLE VIGo). As for groups 6 and 7, however, the grouping of ITS was inconsistent to that of ß-tubulin: the ß-tubulin distances among isolates indicating ITS types 6a or 6b vs. 6c (CBS 345.33, 348.64, 553.68, 576.68. UAMH 3391, 4691, 9768 vs. CBS 292.60) were 0.036–0.042 and 0.041 between two groups of three isolates, indicating ITS type 7a (CBS 314.56 vs. UAMH 9248, NBRC 9011). These values corresponded with that between group vs. group for groups 1–5 (TABLES VIGo). In other words, these isolates accumulated genetic variations in the ß-tubulin locus but few or no genetic changes in ITS. Therefore, the strain CBS 292.60 indicating ITS type 6c is treated as a member of ß-tubulin group 8 (TABLE IIGo). Similarly, UAMH 9248 and NBRC 9011, indicating ITS type 7a, are treated as ß-tubulin group 9 (TABLE IIGo). Consequently, the genetic distance within a group ranged 0.003–0.028 and that between groups was 0.031–0.094 overall. This grouping was supported by phylogenetic analyses.


View this table:
[in this window]
[in a new window]
 
TABLE VI. Uncorrected distances between alleles in the ß-tubulin gene of Thysanophora penicillioides

 
A maximum parsimony analysis was undertaken to reveal phylogenetic relationships among alleles found in the study. The ITS data set including gap sites as a fifth character consisted of 609 sites with 21 parsimony-informative sites and eight generated phylogenetic trees (length = 34 steps, CI = 0.824, RI = 0.961, RC = 0.791, HI = 0.176). Six of eight trees showed that six ITS groups (Groups 1–6) classified based on ITS sequence variations were clustered into independent clades. The remaining trees showed that no type 2c was included in the group 2 clade. All trees indicated that type 7a was the ancestral type of group 6 on MP analysis. Bootstrap analyses supported the six clades with 58–100% probability (FIG. 3Go). On the other hand, 18 MP trees were constructed using the ITS data set without the gap sites (603 sites including 15 parsimony-informative sites, length = 26 steps, CI = 0.808, RI = 0.957, RC = 0.773, HI = 0.192). Although the bootstrap values of groups 1–3 and 5 were not different from the data set including gaps, that of group 4 and group 6 clades were very low in comparison to it (FIG. 3Go).


Figure 3
View larger version (13K):
[in this window]
[in a new window]
 
FIG. 3. The midpoint rooting tree of Thysanophora penicillioides estimated from the ITS data set by the MP analysis. Strains used in the analysis and their ITS types are figured in the terminal branches. The numbers above branches are bootstrap values from 1000 replicates and the numbers below branches are branch lengths. The branch lengths are shown below branches. Each number of ITS group of T. penicillioides is shown in the right of tree.

 
The ß-tubulin data set with gap sites excluded had fewer sites (346 sites, 67 parsimony-informative sites), but showed more variance than the ITS data sets. Four trees were constructed using MP analysis (length = 103 steps; CI = 0.845, RI = 0.950, RC = 0.802, HI = 0.155). Those trees indicated five ß-tubulin groups for Japanese isolates with more than 80% probability. Members of their ß-tubulin groups corresponded to the members of respective ITS groups (FIG. 4Go). In the group 1 clade, three subclades were generated along with two subclades in the group 2 clade. Of the isolates indicating the ITS group 6, four isolates of ITS types 6a or 6b were clustered to one clade comprising two subclades. However, one isolate (CBS 292.60) of ITS type 6c was not clustered with it; its lineage was treated as 8th group on the ß-tubulin tree (FIG. 4Go). Isolates of ITS group 7 were divisible into two lineages that were genetically distant from each other. These two lineages were not supported with a significant bootstrap value and one lineage was treated as 9th group (FIG. 4Go). Relationships among these four lineages were unknown by analyses of the ß-tubulin data set.


Figure 4
View larger version (15K):
[in this window]
[in a new window]
 
FIG. 4. The midpoint rooting tree of Thysanophora penicillioides estimated from the ß-tubulin data set by MP analysis. Strains used in this analysis and their ITS types are figured in the terminal branches. The numbers above branches are bootstrap values from 1000 replicates and the numbers below branches are branch lengths. The branch lengths are shown below branches. Each number of the ß-tubulin group of T. penicillioides is shown in the right of tree.

 
A partition homogeneity test indicated the datasets could include ITS and ß-tubulin data sets (P value when ITS gap sites were included/deleted = 0.136/ 0.202). The combined data set, which was deleted in gap sites (949 sites, 82 parsimony-informative sites) generated two MP trees (length = 131 steps; CI = 0.824, RI = 0.947, RC = 0.781, HI = 0.176). Although they were similar to results of the ß-tubulin analysis, the bootstrap value of each cluster increased: clades 1–6 and 9 were supported by high bootstrap values (more than 93%). In addition, Japanese isolates and Europe/North American isolates were clustered separately with 84% probability.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Intraspecific diversity of Thysanophora penicillioides.— That study evaluated the intraspecific diversity of the fungus based on Japanese, European, and North American isolates. The fungus, a well defined morphological species (Kendrick 1961Go) was found to have uncommon genetic variations within a species. Respective genetic variations of regions among alleles were 8.57%, 5.56% and 6.29% in ITS1, ITS2, and total. These values apparently differ from those of related taxa belonging to Penicillium subgenus Penicillium, a close relativity of the genus Thysanophora (Iwamoto et al 2002Go), for which very low variations at ITS regions were indicated (Skouboe et al 1999Go, Scott et al 2004Go, Peterson 2004Go). The phylogenetic analyses showed that multiple lineages were recognizable in the fungus. Especially, the combined data of ITS and ß-tubulin strongly support nine lineages with high bootstrap values (93–100%, FIG. 5Go). According to phylogenetic species recognition by the concordance of genealogies (GCPSR, Taylor et al 2000Go), each lineage corresponds to a phylogenetic species. However, we should treat the fungus as a species complex that includes several cryptic species because we have found no clear phenotypic differentiations among lineages. Examples for which morphological characteristics are consistent with species boundaries by GCPSR have been reported for many fungi. In rare cases, combinations of GCPSR and other phenotypic characteristics or geographical information have been used to discern morphologically indistinguishable species complexes (Fisher et al 2002Go, O’Donnell et al 2004Go). Whether this species complex deserves to be separated depends on the result of further morphological, physiological and phylogenetic studies.


Figure 5
View larger version (15K):
[in this window]
[in a new window]
 
FIG. 5. The midpoint rooting tree of Thysanophora penicillioides estimated from the ITS and ß-tubulin combined data set by MP analysis. Strains used in the analysis and their ITS types are figured in the terminal branches. The numbers above branches are bootstrap values from 1000 replicates and the numbers below branches are branch lengths. The branch lengths are shown below branches. The numbers of intraspecific groups recognized by the analysis are shown in the right of tree.

 
Genetic structures of Thysanophora penicillioides species complex.— The 254 isolates used in the study were collected in the natural distributional area of A. mariesii in Japan. It is a very small part in the entire distribution range. Nevertheless, genetically diverse isolates of 12 kinds of ITS sequences categorizable into five ITS groups were found in such a small area. The local population-based survey revealed that, at most, four phylogenetically divergent groups coexisted in a forest. These results indicate that the diversity of the fungi in Japan is very high and phylogenetically differentiated individuals occurring in a morphological species have shared the same niche: decaying fir needles. Cases in which morphologically similar but genetically divergent groups inhabit the same niche have been reported for Botrytis cinerea (Giraud et al 1997Go, 1999Go), Colletotrichum species (Lu et al 2004Go), Lophodermium pinastrii (Müller et al 2001Go), and Phialocephara fortinii (Grünig et al 2004Go). Grünig et al. (2004)Go reported that different cryptic species in the P. fortinii species complex occurred in the same root segments of Norway spruce. The authors isolated two ITS types from a single decaying needle several times (data not shown).

Each lineage of the species complex has undergone gene dispersion without geographic limitation among local populations because the distribution of each lineage shows a wide expansion through the studied areas in A. mariesii forests. However, ITS type 1a was detected frequently in all local populations: it was dominant in most local populations even though multiple ITS types coexisted in all populations. The fact that group 1, including ITS type 1a, is overwhelmingly dominant and other groups are minor in A. mariesii forests of Japan is interesting, but its reason remains unclear. Perhaps members of group 1 acquired more advantageous characteristics for fitness than other groups in the ecosystem.

In marked contrast, the minor ITS types within groups 1 and 2 were recorded only in restricted areas. It is assumed that types 1b, 1c, 1d were derived from type 1a, and types 2c, 2d, 2e from type 2a based on results of parsimony analysis (FIGS. 3Go–5GoGo). Probably, they were generated by a recent point mutation at ITS of types 1a or 2a. Therefore, their distributional areas might remain small even though they began to disperse from their places of origin.

We initially surmised that any clonal lineage might be distributed around the world and that various lineages found in Japanese population were also detected in overseas samples. Thereby, the fungus might be capable of dispersing its propagules throughout the world for a relatively short period because of abundant production of small dry conidia; it might be colonizable on fallen needles of various pinaceous trees that are distributed over a wide area. For sequencing the strains out of Japan, we obtained strains of Europe and North America via public culture collections. Results showed that ITS types of European and North American isolates differed from Japanese isolates and that they formed separate clades for Japanese isolates in phylogenetic trees. Against our expectations, this fact suggests that genetic diversification occurred independently in two separate regions and that numerous differentiated populations of the fungus might exist on a global scale. Although genetic structures among three separate regions might be different, we cannot draw a conclusion about the relationships between genetic structures of the three regions because the isolates obtained from European and North American populations are too few.

This study detected numerous genetic variations using ITS sequence polymorphism in an anamorphic fungus inhabiting leaf litter of Pinaceaeous forests. Results revealed genetically diverged groups, suggesting a species complex including multiple cryptic species. These facts foreshadow that more genetic variations for the species might be detectable in other areas.


    ACKNOWLEDGMENTS
 
We thank Dr. S. Inaba and Dr. H. Saito for assisting collection of litter samples in the Tohoku district.


    FOOTNOTES
 
Accepted for publication November 11, 2005.

1 Corresponding author, Present address: BioFrontier Laboratories, Kyowa Hakko Kogyo Co. Ltd., Asahi-machi, Machida, Tokyo 194-8533, Japan, E-mail: susumu.iwamoto{at}kyowa.co.jp


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Aoki T, Tokumasu S, Tubaki K. 1990. Fungal succession on momi fir needles. Trans Mycol Soc Japan 31:355–374.

Boeger JM, Chen RS, McDonald BA. 1993. Gene flow between geographic populations of Mycosphaerella graminicola (anamorph Septria tritici) detected with restriction fragment length polymorphism markers. Phytopathol 83:1148–1154.

Boucias D, Clare S, Suazo A, Funderburk J. 2000. AFLP analysis of the entomopathogen Nomuraea rileyi. Mycologia 92:638–648.

Carbone I, Kohn LM. 2001. A microbial population-species interface: nested cladistic and coalescent inference with multilocus data. Mol Ecol 10:947–964.

Coetzee MPA, Wingfield BD, Harrington TC, Dalevi D, Coutinho TA, Wingfield MJ. 2000. Geographical diversity of Armillaria mellea s. s. based on phylogenetic analysis. Mycologia 92:105–113.

Fargues J, Bon M-J, Manguin S, Couteaudier Y. 2002. Genetic variability among Paecilomyces fumosoroseus isolates from various geographical and host origins based on the rDNA-ITS regions. Mycol Res 106:1066–1074.

Felsenstein J. 1985. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 39:783–791.

Fisher MC, Koenig GL, White TJ, Taylor JW. 2002. Molecular and phenotypic description of Coccidioides posadasii sp. nov., previously recognized as the non-California population of Coccidioides immitis. Mycologia 94:73–84.

Geiser DM, Frisvad JC, Taylor JW. 1998. Evolutionary relationships in Aspergillus section Fumigati inferred from partial ß-tubulin and hydrophobin DNA sequences. Mycologia 90:831–845.

Giraud T, Fortini D, Levis C, Lamarque C, Leroux P, LoBuglio K, Brygoo Y. 1999. Two sibling species of the Botrytis cinerea complex, transpose and vacuma, are found in sympatry on numerous host plants. Phytopathol 89:967–973.

———, ———, ———, Leroux P, Brygoo Y. 1997. RFLP markers show genetic recombination in Botryotinia fuckeliana (Botrytis cinerea) and transposable elements reveal two sympatric species. Mol Biol Evol 14:1177–1185.

Gräser Y, Volovsek M, Arrington J, Schönian G, Presber W, Mitchell TG, Vilgalys R. 1996. Molecular markers reveal that populations structure of the human pathogen Candida albicans exhibits both clonality and recombination. Proc Natl Acad Sci USA 93:12473–12477.

Grünig CR, McDonald MA, Sieber TN, Rogers SO, Holdenrieder O. 2004. Evidence for subdivision of the root-endophyte Phialocephara fortinii into cryptic species and recombination within species. Fungal Gent Biol 41:676–687.

Hayashi Y. 1960. Taxonomical and phytogeographical study of Japanese conifer. Tokyo: Norin Syuppan.

Hughes KW, Toyohara TL, Peterson RH. 1998. DNA sequence and RFLP analysis of Pleurotopsis longinqua from three disjunct populations. Mycologia 90:595–600.

Iwamoto S, Tokumasu S. 2001. Dematiaceous hyphomycetes inhabiting decaying blackish needles of Abies firma and their distribution in the Kanto district, Japan. Mycoscience 42:273–279.

———, ———, Suyama Y, Kakishima M. 2002. Molecular phylogeny of four selected species of the strictly anamorphic genus Thysanophora using nuclear ribosomal DNA sequences. Mycoscience 43:169–180.

James TY, Vilgalys R. 2001. Abundance and diversity of Schizophyllum commune spore clouds in the Caribbean detected by selective sampling. Mol Eco 10:471–479.

Jensen AB, Eilenberg J. 2001. Genetic variation within the insect-pathogenic genus Entomophthora, focusing on the E. muscae complex, using PCR-RFLP of the ITS and the LSU rDNA. Mycol Res 105:307–312.

Kendrick WB. 1961. Hyphomycetes of conifer leaf litter. Thysanophora gen. nov. Can J Bot 39:817–832.

Koenig RL, Ploetz RC, Kistler HC. 1997. Fusarium oxysporum f. sp. cubense consists of a small number of divergent and globally distributed clonal lineages. Phytopathol 87:915–923.

Kominami K, Kobayashi Y, Tubaki K. 1952. Enumerations of the moulds of Japan. 1. species of fungi imperfecti. Nagaoa 1:7–25.

Lu G, Cannon PF, Reid A, Simmons CM. 2004. Diversity and molecular relationships of endophytic Colletotrichum isolates from the Iwokrama Forest Reserve, Guyana. Mycol Res 108:53–63.

Milgroom MG. 1996. Recombination and the multilocus structure of fungal populations. Ann Rev Phytopathol 34:457–477.

———, Wang K, Zhou Y, Lipari SE, Kaneko S. 1996. Intercontinental population structure of the chestnut blight fungus, Cryphonectria parasitica. Mycologia 88:179–190.

Müller MM, Valjakka R, Suokko A, Hantula J. 2001. Diversity of endophytic fungi of single Norway spruce needles and their role as pioneer decomposers. Mol Ecol 10:1801–1810.

Nei M. 1973. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci USA 70:3321–3323.

O’Donnell K, Ward TJ, Geiser DM, Kistler HC, Aoki T. 2004. Genealogical concordance between the mating type locus and seven other nuclear genes supports formal recognition of nine phylogenetically distinct species within the Fusarium graminearum clade. Fungal Genet Biol 41:600–623.

Peterson SW. 2004. Multilocus DNA sequence analysis shows that Penicillium biourgeianum is a distinct species closely related to P. brevicompactum and P. olsonii. Mycol Res 108:434–440.

Rogers SO, Rogers MAM. 1999. Gene flow in fungi. In: Worrall JJ. ed. Structure and dynamics of fungal populations. Dordrecht, UK: Kluwer Academic Publishers. p 97–121.

Scott J, Untereiner WA, Wong B, Straus NA, Malloch D. 2004. Genotypic variation in Penicillium chrysogenum from indoor environments. Mycologia 96:1095–1105.

Sieber-Canavesi F, Sieber TN. 1993. Successional patterns of fungal communities in needles of European silver fir (Abies alba Mill.). New Phytologist 125:149–161.

Skouboe P, Frisvad JC, Taylor JW, Lauritsen D, Boysen M, Rossen L. 1999. Phylogenetic analysis of nucleotide sequences from the ITS region of terverticillate Penicillium species. Mycol Res 103:873–881.

Suyama Y, Kawamuro K, Kinoshita I, Yoshimura K, Tsumura Y, Takahara H. 1996. DNA sequence from a fossil pollen of Abies spp. from Pleistocene peat. Genes Genet Syst 71:145–149.

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

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

Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. 1997. The ClustalX Windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nuc Acids Res 24:4876–4882.

Tokumasu S. 1980. Observations on the fungal flora in pine leaf litter In: Ecology of microorganisms. Vol. 7. Tokyo: Gakkai Shuppan Center. p 130–149.

———. 1996. Mycofloral succession on Pinus densiflora needles on a moder site. Mycoscience 37:313–321.

Trigiano RN, Caetano-Anolles GC, Bassam BJ, Windham MT. 1995. DNA amplification fingerprinting provides evidence that Discula destructive, the cause of dogwood anthracnose in North America, is an introduced pathogen. Mycologia 87:490–500.

van Maanen A, Gourbière F. 1997. Host and geographical distribution of Verticicladium trifidum, Thysanophora penicillioides, and similar fungi on decaying coniferous needles. Can J Bot 75:699–710.

Xu J, Kerrigan RW, Callac P, Horgen PA, Anderson JB. 1997. Genetic structure of natural populations of Agaricus bisporus, the commercial button mushroom. Journal of Heredity 88:482–488.

———, Vilgalys R, Mitchell TG. 2000. Multiple gene genealogies reveal recent dispersion and hybridization in the human pathogenic fungus Cryptococcus neoformans. Mol Ecol 9:1471–1481.





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


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS