| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
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 |
|---|
|
|
|---|
We investigated intraspecific diversity and genetic structures of a saprotrophic fungusThysanophora penicillioidesbased 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 |
|---|
|
|
|---|
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 2001
). Another obstacle is the difficulty in defining fungal populations because direct confirmation of their individuals is lacking in the field (Rogers and Rogers 1999
). 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 habitatdecaying 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 |
|---|
|
|
|---|
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 1960
). The distributions in the forests are discontinuous (FIG. 1
).
|
|
Ten needles per site were used for fungal detection. Brown or partly darkened brown needles were selected from each litter sample. A modified Tokumasus washing method (Tokumasu 1980
) 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 II
).
|
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 1998
). 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 manufacturers 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). Neis gene diversity formulae (Hep) were used to evaluate the genetic diversity of populations (Nei 1973
). The uncorrected distance (p-distance) between sequences of ITS or ß-tubulin were calculated using PAUP version 4.0b10 (Swofford 2002
).
Phylogenetic analysis.
Three data setsITS, ß-tubulin, and ITS and ß-tubulinwere prepared of combined data including all ITS types of T. penicillioides (TABLE II
, III
). Alignments of data sets were determined using Clustal X version 1.81 (Thompson et al 1997
) and were refined visually.
|
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 II
, III
).
| RESULTS |
|---|
|
|
|---|
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. 2
). Finally, 12 distinctive sequences (ITS types 1a5a) were detected. These ITS types were classified into five groups in which groups 1 and 2 were included more than two ITS types (FIG. 2
). Frequencies of respective ITS types and diversity indexes (Hep) within local populations are shown in TABLE IV
. 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%).
|
|
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 (2c2e); they were found only in this region. Minor ITS types of group 1 (1b1d) 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. 2
). 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 II
). 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 II
).
Genetic distances and Phylogenetic relationships.
TABLE V
shows genetic distances between ITS types. The distances within each ITS group (for groups 1, 2 and 6) are 0.0020.003. On the other hand, the genetic distance between ITS groups was an average of 0.013 for values of 0.0030.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. 2
).
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
|
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. 3
5![]()
). 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 |
|---|
| FOOTNOTES |
|---|
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 |
|---|
|
|
|---|
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:11481154.
Boucias D, Clare S, Suazo A, Funderburk J. 2000. AFLP analysis of the entomopathogen Nomuraea rileyi. Mycologia 92:638648.
Carbone I, Kohn LM. 2001. A microbial population-species interface: nested cladistic and coalescent inference with multilocus data. Mol Ecol 10:947964.
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:105113.
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:10661074.
Felsenstein J. 1985. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 39:783791.
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:7384.
Geiser DM, Frisvad JC, Taylor JW. 1998. Evolutionary relationships in Aspergillus section Fumigati inferred from partial ß-tubulin and hydrophobin DNA sequences. Mycologia 90:831845.
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:967973.
, , , 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:11771185.
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:1247312477.
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:676687.
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:595600.
Iwamoto S, Tokumasu S. 2001. Dematiaceous hyphomycetes inhabiting decaying blackish needles of Abies firma and their distribution in the Kanto district, Japan. Mycoscience 42:273279.
, , Suyama Y, Kakishima M. 2002. Molecular phylogeny of four selected species of the strictly anamorphic genus Thysanophora using nuclear ribosomal DNA sequences. Mycoscience 43:169180.
James TY, Vilgalys R. 2001. Abundance and diversity of Schizophyllum commune spore clouds in the Caribbean detected by selective sampling. Mol Eco 10:471479.
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:307312.
Kendrick WB. 1961. Hyphomycetes of conifer leaf litter. Thysanophora gen. nov. Can J Bot 39:817832.
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:915923.
Kominami K, Kobayashi Y, Tubaki K. 1952. Enumerations of the moulds of Japan. 1. species of fungi imperfecti. Nagaoa 1:725.
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:5363.
Milgroom MG. 1996. Recombination and the multilocus structure of fungal populations. Ann Rev Phytopathol 34:457477.
, Wang K, Zhou Y, Lipari SE, Kaneko S. 1996. Intercontinental population structure of the chestnut blight fungus, Cryphonectria parasitica. Mycologia 88:179190.
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:18011810.
Nei M. 1973. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci USA 70:33213323.
ODonnell 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:600623.
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:434440.
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 97121.
Scott J, Untereiner WA, Wong B, Straus NA, Malloch D. 2004. Genotypic variation in Penicillium chrysogenum from indoor environments. Mycologia 96:10951105.
Sieber-Canavesi F, Sieber TN. 1993. Successional patterns of fungal communities in needles of European silver fir (Abies alba Mill.). New Phytologist 125:149161.
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:873881.
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:145149.
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:2132.
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:48764882.
Tokumasu S. 1980. Observations on the fungal flora in pine leaf litter In: Ecology of microorganisms. Vol. 7. Tokyo: Gakkai Shuppan Center. p 130149.
. 1996. Mycofloral succession on Pinus densiflora needles on a moder site. Mycoscience 37:313321.
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:490500.
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:699710.
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:482488.
, Vilgalys R, Mitchell TG. 2000. Multiple gene genealogies reveal recent dispersion and hybridization in the human pathogenic fungus Cryptococcus neoformans. Mol Ecol 9:14711481.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |