| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Plant Pathology and Institute of Molecular Evolutionary Genetics, Pennsylvania State University, 211 Buckhout Laboratory, University Park, Pennsylvania 16802
Donald D. Davis
Department of Plant Pathology, Pennsylvania State University, 211 Buckhout Laboratory, University Park, Pennsylvania 16802
David M. Geiser
Department of Plant Pathology and Institute of Molecular Evolutionary Genetics, Pennsylvania State University, 211 Buckhout Laboratory, University Park, Pennsylvania 16802
| ABSTRACT |
|---|
|
|
|---|
Despite mycologists interest in its unique spore-dispersal mechanism, systematic studies of the genus Sphaerobolus have received little attention. In our previous work multiple gene genealogies indicated the existence of three divergent lineages in the genus Sphaerobolus, each representing a phylogenetic species. Macro- and micromorphological analyses of colony and fruit-body characters presented here confirmed that these three phylogenetic species correspond to two known species, S. iowensis and S. stellatus, and a newly discovered species. In addition, an expanded gene genealogical analysis is presented for the three species. The new species, named Sphaerobolus ingoldii Geml, Davis et Geiser, is described based on both molecular and morphological data. In addition, while S. iowensis previously had been reported in only two localities, we found that it is as common as or more common than S. stellatus in North America. Despite the considerable amount of DNA polymorhism found in all species, nested clade analyses of S. iowensis and S. stellatus indicated little phylogeographic structure in either species, perhaps due to heavy movement mediated by human activities.
Key words: phylogenetics, Sphaerobolus ingoldii, systematics
| INTRODUCTION |
|---|
|
|
|---|
Because it produces its basidiospores in an enclosed structure, the genus Sphaerobolus has been classified as a member of the class Gasteromycetes along with fungi that undergo passive spore discharge, including birds nest fungi (Cyathus, Crucibulum etc.), puffballs (Lycoperdon, Langermannia etc.), earth balls (Scleroderma etc.), stinkhorns (Phallus, Mutinus, Pseudocolus etc.), and earth stars (Geastrum etc.). Within this form class authors have placed the artillery fungus in family Sphaerobolaceae in these orders: Sphaerobolales (Ainsworth 1971
); Nidulariales (Ulloa and Hanlin 2000
); Sclerodermatales (Hawksworth et al 1996
). However, molecular phylogenetic analyses of Hibbett et al (1997)
revealed that the Gasteromycetes represents a polyphyletic, artificial taxon that does not represent a true evolutionary group within the Basidiomycota. Based on molecular data, the genus Sphaerobolus is best placed in the gomphoid-phalloid clade with these genera as closest relatives: Geastrum, Phallus, Pseudocolus, Ramaria, Clavariadelphus, Gomphus, Gautieria etc. (Hibbett et al 1997
, Moncalvo et al 2002
). These results are incorporated in the classification of Kirk et al (2001)
, who placed Sphaerobolus in the family Geastraceae, order Phallales, class Basidiomycetes, phylum Basidiomycota.
The artillery fungus produces spherical, whitish or orange basidiocarps (12 mm diam), in which a single spore mass (gleba or peridiole) is formed. The gleba contains numerous spores of two types, uninucleate, thick-walled basidiospores and thin-walled, elongated gemmae (asexual spores) with two or more nuclei (Walker 1927
, Ingold 1972
). These two types presumably play different adaptive roles in the dual lignicolous and coprophilous survival mechanism of this fungus. Basidiospores are stimulated to germinate when exposed to proteolytic enzymes (e.g. pepsin) and relatively high temperatures as expected in the digestive system of mammal herbivores after ingestion of a discharged gleba attached to plant material (Dykstra 1982
). If no herbivore ingests the glebal mass, the gemmae can germinate directly on wood or other plant debris. Glebae have been reported to germinate after at least 11 y of storage (Walker 1927
) and thus might serve as long-term survival structures. After the germination of gemmae, colonization of the substrate begins. The characteristic bleaching associated with this white rotting fungus can be attributed to the digestion of cellulose, hemicellulose and lignin in the substrate. After colonizing the substrate for some time the fungus produces phototropic fruiting bodies. At maturity these fruiting bodies orient themselves toward the brightest light source and shoot their glebae into air (Alasoadura 1963
, Nawaz 1967
). Similar active/passive mechanisms for dispersing spore masses can be observed in other dung-inhabiting fungi as well, including Pezizomycetes (e.g. Ascobolus) and Sordariomycetes (e.g. Podospora) of the Ascomycota and Mucorales of the Zygomycota (e.g. Pilobolus). This feature is likely an adaptation to the coprophilous habit and evolved many times independently during evolution.
The Greek origin of the name Sphaerobolus ("sphere thrower") is related to the ability of the fungus to shoot its glebae a considerable distance. The genus contains two currently recognized species, S. stellatus (Tode) Pers. and S. iowensis Walker (Hawksworth et al 1996
), that are distinguished by a number of micromophological characteristics. In the developing glebae of young S. stellatus basidiocarps, some of the binucleate hyphae at the center of the globular knots have enlarged, clavate basidial primordia (Walker 1927
). New basidia proliferate in clusters, each eventually producing four (occasionally eight) basidiospores. The basidia decompose and disappear as the spores mature, making room for the enlargement and maturation of other basidia. In S. iowensis the development of primary basidia is followed by the formation of characteristic cavities or chambers in which basidiospores are produced (Walker 1927
). In both species sterile hyphae undergo disintegration in maturing basidiocarps and break down by the time the glebal mass is discharged. Almost all studies in the genus have concentrated on isolates identified as S. stellatus. However, several other names of uncertain origin appear in culture collections, including S. bombardioides, S. carpobolus, S. corii, S. crustacus, S. epigaeus, S. minimus, S. minutissimus, S. rubidus, S. sparsus and S. tubulosus. The origin of these binomials cannot be found in the literature and thus cannot be considered scientifically valid. Sphaerobolus is a cosmopolitan genus reported from Alaska, British Columbia, most parts of the continental U.S. (it is particularly abundant in the Northeast), Europe (from Greece to Iceland), Asia (including Japan), Australia, New Zealand, Africa and Latin America (Aplin 1961
, Dring 1964
, Ingold 1972
, Herrera and Perezsilva 1987
, McKenzie and Foggo 1989
, Halgrimsson et al 1992
, Zervakis et al 1998
, McGray pers comm 2002). Whether this is its natural range or whether it has been dispersed artificially is unknown.
In our studies (Geml et al 2005
) we generated nucleotide sequences of three ribosomal DNA regions and one protein coding gene region (partial mitochondrial ribosomal RNA small subunit [mtSSU], the internal transcribed spacer regions of the nuclear ribosomal gene repeat [ITS], partial nuclear large ribosomal RNA subunit [LSU] and partial translation elongation factor 1-
gene [EF 1-
]). Multiple gene genealogies inferred from maximum likelihood, Bayesian and maximum parsimony analyses of sequence data from individual loci (mtSSU, ITS, EF 1-
, and LSU) and a combined dataset (mtSSU, ITS, EF 1-
) concordantly indicated the existence of three divergent lineages in the genus Sphaerobolus, each representing a phylogenetic species. In this paper we have investigated variation in morphological characteristics corresponding to the phylogenetic species revealed, associating them with new and known taxa. We have included additional isolates of the newly discovered species for phylogenetic analyses and described it based on both molecular and morphological data. We also have sampled various artillery fungus communities and identified them based on ITS sequences to reveal more information on the geographic distribution of the three species in the U.S. In addition we analyzed the nucleotide variation in ITS sequence in S. iowensis and S. stellatus to investigate their biogeography using nested clade analyses.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
, and the entire ITS region were PCR amplified in reaction mixtures containing 37 µL PCR water, 5 µL 10x PCR buffer (0.5M KCl, 0.1M Tris HCl pH 8.3, 0.025M MgCl2), 5 µL 10x dNTPs (2mM of each dNTP), 0.1 µL AmpliTaq® DNA polymerase (Perkin-Elmer, Foster City, California), 1 µL of 10 µM forward primer and reverse primer for the region of interest, and 1 µL template DNA (100-fold dilution of original DNA solution extracted). PCR reactions were performed in a 96-well thermocycler (PTC-100 Programmable Thermal Controller, MJ Research Inc.) using this temperature program for all ribosomal DNA regions: 94 C/5 min; 34 cycles of 94 C/1 min, 53 C/1 min, 72 C/1 min; and 72 C/5 min. For the single copy gene EF 1-
a "touchdown" PCR setting was used with an annealing temperature of 65 C in the first cycle, then successively reduced by 1 C per cycle to 56 C, after which the annealing temperature was maintained at 56 C for the remaining 3036 cycles (Stephen Rehner pers comm). These primers were used for amplification: ITS5 and ITS4 for ITS (White et al 1990
(Stephen Rehner pers comm). Amplification products were electrophoresed in 3.0% agarose gels and stained with ethidium bromide for visualization of the bands. PCR products were purified directly with the QIAquick® PCR Purification Kit (QIAGEN Inc., Valencia, California) and were sequenced with the Applied Biosystems (ABI) BigDye® v. 3.0 terminator kit and an ABI 377 automated DNA sequencer (Perkin-Elmer, Foster City, California). Each sample was sequenced with the same primers used for PCR with the exception of the 1567R primer of the EF 1-
gene that was replaced by primer 1567Ra (ACH GTR CCR ATA CCA CC) in sequencing for better results (Stephen Rehner pers comm).
Phylogenetic analyses.
Sequence data obtained for both strands of each locus were edited and assembled for each isolate with Sequencher 3.1 (Gene Codes, Ann Arbor, Michigan). Sequence alignments were made by Clustal X (Thompson et al 1997
). Analyses were conducted with maximum likelihood (ML) and maximum parsimony (MP) methods in PAUP* 4b10 (Swofford 2002
) and Bayesian analysis in MrBayes 3.0 (Huelsenbeck and Ronquist 2001
). Because methods follow different theories and algorithms, congruent features found in all three types of analyses were considered meaningful. At first we excluded ambiguously aligned regions, consisting of these positions: 4763, 205255, 296308 in the mtSSU; 7984, 160221, 487499, 623666, 699703 in the ITS; and 203260, 408460, 558562 in the EF 1-
alignments. Alignments were deposited in the TreeBase database (http://www.treebase.org) under accession numbers SN2111-7316, SN2111-7317 and SN2111-7320, and sequences were deposited in GenBank (TABLES I
and II
). Different evolutionary models with varying values of base frequencies, substitution types,
-parameter of the
-distribution of variable sites and proportion of invariable sites, among other parameters, were compared via the likelihood ratio test for each locus with PAUP* and Modeltest 3.06 (Posada and Crandall 1998
) to determine the best-fit evolutionary model for both ML and Bayesian analyses. ML analyses were carried out with the heuristic search option using the tree bisection and reconnection (TBR) algorithm with 100 random sequence additions to find the overall optimum instead of local optima. To determine the statistical reliability of the trees and the stability of clades, the bootstrap test (Felsenstein 1985
) was used with full heuristic search and 100 replicates. In Bayesian phylogenetic analyses, 100 000 generations were run in four chains. The chains were sampled every 100th generation. When the likelihood scores of trees sampled approached similar values, they were considered to have converged. In each run only trees after this convergence point were included in computing the consensus tree. In the MP analyses previously excluded ambiguous regions were included after being recoded with the program INAASE 2.3b (Lutzoni et al 2000
). The code and step matrices were attached to the appropriate alignments. Gaps were treated as "new state". The bootstrap test was used with 500 replicates.
|
ITS-based identification and biogeography of species.
Eighty isolates were sampled in various parts of the U.S. or received from collections. The isolation and growth of pure cultures, DNA extraction, PCR amplification and sequencing of the ITS region were conducted in the same manner described above. Geographic distribution, frequency and GenBank accession numbers for ITS haplotype sequences are provided (TABLE II
). Species identities of sampled isolates were made based on sequence similarity to various phylogenetic species determined by our earlier work. For this purpose we added newly generated sequences to an existing alignment (Geml et al 2005
) and analyzed them with MP. Species identifications were inferred based on the phylogenetic placement of the sequences within one of the three species clades. Biogeographic patterns linked to the different phylogenetic species were investigated with nested clade analysis (NCA) (Templeton 1998). Maximum parsimony haplotype networks were generated by TCS v.1.13 (Clement et al 2000) and were used to define a series of nested clades that in turn were used to perform random, two-way contingency permutation analysis to detect any association between geographic distribution and genetic variation (Templeton 1998). The nested clade information, sample size for each haplotype and geographic location of each clade (latitude and longitude coordinates) were entered into the software package GeoDis v.2.0 (Posada et al 2000). GeoDis was used to calculate clade distance (Dc) and nested clade distance (Dn) and to test them for significance at the 5% level with a permutation technique with 1000 resampling replicates (Posada et al 2000). Dc was calculated as the average distance of all individuals in clade X from the geographic center of that clade, while Dn was the average distance of individuals in clade X from the geographic center of clades of the next highest nesting level. Where significant Dc and/or Dn values were detected, a set of criteria was used to discriminate between the effects of contemporary (e.g. gene flow) and historical (e.g. allopatric fragmentation, range expansion) processes (Posada et al 2000, Templeton 1998). Because only five S. ingoldii isolates were collected, this species was excluded from the nested clade analyses.
| RESULTS |
|---|
|
|
|---|
and the combined datasets respectively consisted of 583, 769, 628, and 1980 characters, including gaps. Of these 63, 78, 110 and 251 ambiguous positions were excluded in the intrageneric analyses. Of interest, we found that regions that could not be aligned unambiguously between the three different phylogenetic species were highly similar and sometimes identical within lineages, indicating the existence of fixed differences. Three types of long (>1 kb), independent insertions, each with a unique insertion site, length and nucleotide sequence, were detected in the mtSSU region of some of the isolates (Geml et al 2005
-parameter of
-distribution (TN + G) was selected as the best-fit evolutionary model for the EF 1-
dataset. The base frequencies were freqA = 0.2349, freqC = 0.2869, freqG = 0.2492, freqT = 0.2290. The substitution rate matrix was: R(a) [A-C] = 1.0000, R(b) [A-G] = 2.7043, R(c) [A-T] = 1.0000, R(d) [C-G] = 1.0000, R(e) [C-T] = 11.3650, R(f) [G-T] = 1.0000. The among-site rate variation was characterized by variable sites following a
-distribution with an estimated
-parameter = 0.3767. Because the HKY model was determined earlier as best-fit model for two of the three loci, that model was chosen for ML and Bayesian analyses of the combined dataset as well.
Likelihood values converged after about 16 000, 20 000, 19 000 and 16 000 generations in Bayesian analysis of the mtSSU, ITS, EF 1-
and combined datasets, respectively. The consensus trees were computed from 840, 800, 810 and 840 trees, after discarding the first 161, 201, 191 and 161 trees as burn-in. The ML phylograms and the Bayesian consensus trees generated with the models described earlier had similar topologies at each locus. After including the character matrices of the ambiguous regions recoded by INAASE, the final alignments of mtSSU, ITS, and EF 1-
for MP analyses consisted respectively of 522, 697, 520 and 1739 characters. One of the 26 most parsimonious trees of the combined dataset is shown with bootstrap values (FIG. 1
).
|
Morphological examinations.
Morphological and growth characteristics of the three phylogenetic species of Sphaerobolus are summarized (TABLE III
). Measurements were taken from 25 glebae per isolate (50 glebae per species). The only exception was S. iowensis, for which only a single fruiting culture (ATCC 52850) was obtained, producing only eight basidiocarps. All glebae recorded were produced on WSA. Although S. ingoldii and S. stellatus fruited abundantly on OA as well, S. iowensis failed to produce basidiocarps on OA, despite repeated efforts. Therefore this medium, as well as PDA that did not support any fruiting, was excluded from the gleba and basidiospore morphology evaluation. We found that significantly smaller (P < 0.001) glebae were produced by S. ingoldii (0.977 ± 0.198 mm) than produced by the other two species (S. iowensis: 1.549 ± 0.156, S. stellatus: 1.495 ± 0.289). Microscopic examinations revealed the existence of basidial chambers, a distinct character used in the species description by Walker (1927)
, in the premature glebae of S. iowensis. Neither of the other species exhibited this trait. We found significant differences in the length and width of basidiospores produced (both P < 0.001). The largest spores were produced by S. ingoldii (8.780 ± 0.761 x 5.874 ± 0.583 µm), followed by S. iowensis (7.215 ± 0.510 x 4.917 ± 0.217 µm) and S. stellatus (7.266 ± 0.498 x 4.571 ± 0.288 µm) with nonsignificant difference between the latter two. Growth rates of the three species on OA and WSA were significantly different (both P < 0.001), while we failed to reject the null hypothesis (i.e. all rates are equal) when analyzing growth on PDA (P = 0.113). Sphaerobolus ingoldii showed the fastest growth on both OA and WSA, with daily linear growth of 2.095 ± 0.036 mm and 2.089 ± 0.029 mm, respectively. The growth rate of S. iowensis was 1.310 ± 0.241 mm (OA) and 1.221 ± 0.081 mm (WSA), while S. stellatus grew at a daily rate of 1.381 ± 0.356 mm (OA) and 1.213 ± 0.288 mm (WSA). Growth rates on PDA for S. ingoldii, S. iowensis and S. stellatus were respectively 1.345 ± 0.425 mm, 1.025 ± 0.131 mm and 1.101 ± 0.484 mm, (TABLE III
). Morphological characteristics of S. ingoldii, S. iowensis and S. stellatus are shown (FIGS. 2
5, 6
9
and 1012
, respectively).
|
|
|
|
|
|
|
| TAXONOMY |
|---|
|
|
|---|
Peridium globosum,
11.5 mm diam. Basidiomata gregaria, aliquando solitaria. Superficies exoperidii saepe rimosa, subreticulata. Gleba brunnea aut nigrofusca, 0.61.3 mm diam. Basidiosporae hyalinae, glabrae, ellipsoideae, crassitunicatae, (7.2)89.5(10.2) x (5)5.46.2(7.1) µm, gemmae absentes. Coloniae in OA et WSA 2.02.2 mm, et PDA 0.91.8 mm per diem crescens. Color coloniarum in OA albus, eburneus aut alboluteus, mycelio saepe chordaceo.
Peridium globose, rarely exceeding 11.5 mm diam. Basidiocarps most often gregarious, occasionally solitary. The outer surface of the exoperidium is often rimose (cracked) in a reticulated manner. Expelled glebae dark brown to black, with a dull finish, 0.61.3 mm diam, with many basidiospores, but no gemmae. Even the ejected, hence slightly deformed (flattened) gleba rarely is larger than 1 mm. Basidiospores hyaline, smooth, elliptical, thick-walled, (7.2)89.5(10.2) x (5)5.46.2(7.1) µm, borne on basidia irregularly disposed in the gleba. Colonies on OA, WSA and PDA at 25 C exhibit respectively a daily linear growth of 2.02.2 mm, 2.02.2 mm, 0.91.8 mm. Mycelium off-white, sometimes with a yellow shade, often cord-like in appearance on OA.
Characteristic fixed DNA polymorphisms.
Characteristic fixed DNA polymorphisms used in the description of S. ingoldii sp. nov. were determined based on the sequence alignments used in the phylogenetic analyses. Only those characters present in all S. ingoldii isolates and not present in other Sphaerobolus isolates were included. These fixed DNA polymorphisms are indicated with capital and italicized letters (nucleotides) or italicized numbers (gaps) with the alignment position given. EF 1-
: gcaActc @ 161; acc/TtcA/Gac @ 304, 307; gccAtggtAcaa @ 471, 476; aca/GagTccggtgc/tG/AgtGaa @ 499, 502, 510, 513; a[-1-]gCcgtccTTcc @ 566, 572573. ITS: tgCtg @ 68; ggcCttaCacCaa @ 76, 80, 83; tcAAatttcTaa @ 104105, 111; ttAgacgTagAtTGgActcccGgT/CtTct-ATttgg/aagcCgtgagCGg @ 140, 145, 148, 150151, 153, 159, 161, 163, 166167, 175, 181182; tgActCgt @ 220, 223; acAa/gg @ 252; taAtGgg @ 277, 279; cttAta @ 295; ac/taTTaat @ 302303; gcgGttGAaCta @ 538, 541542, 544; aTCgc @ 553554; tatTaTA[-1-]a[-1-]aaCTCTAtggt @ 598, 600601, 607611; aaCacaaAttc[-8-]ttt @ 626, 631, 635642; ctAtG/AgtGAata @ 654, 656, 659660; gaTgt @ 675; gtcCAttgactTgga @ 695696, 703; taTcACTTtgtTatc @ 712, 714717, 721. mtSSU: aaTgCa @ 74 and 76; acCttGctaCGAtAag @ 98, 101, 105107, 109; caTctttGa[-1-]a[-3-]tct[-1-]Ctt[-1-]tt[-4-]gG[-3-]ga @ 200, 205, 207, 209211, 216, 219, 227230; aaCtct[-4-] gTTtGtt @ 292, 296299, 301302, 304; tacACag @ 329330; gcCga @ 399; tcgCtt @ 543. LSU: ggcCtt @ 407; tc[-1-]Ttg @ 427; tt[-1-]Ggac @ 435; tttCgaccaCtg @ 455, 461; gcAAGgg @ 473475; gccCttg[-1-]tTgta @ 511, 517; aggTctg @ 541.
Etymology. The name is a tribute to Charles T. In-gold, an outstanding British mycologist and former president of the British Mycological Society, who conducted fundamental research on the spore dispersal of Gasteromycetes, with a particular interest in Sphaerobolus.
Habitat. On rotten wood, forest litter, herbivore dung.
HOLOTY PE. USA: Michigan. Herbarium specimen (RGT 930825/04) at the USDA Forest Products Laboratory (Madison, Wisconsin) (CFMR). Originally collected and isolated by R.G. Thorn, 25 Aug 1993 as Sphaerobolus stellatus from Michigan, Kellogg Biological Station Long Term Ecological Research site (42°24'N latitude and 85°24'W longitude), near south edge of Treatment 8, on twig of Quercus alba. Living culture ex holotype: T-800 CFMR. CBS MycoBank No. 500167.
Other specimens examined. Additional isolates were received as cultures from R. Hanlin (SS19, University of Georgia), and from K. Yokoyama (IFO 9597, Institute for Fermentation, Osaka, Japan) and were collected as gleba by D.D. Davis (SS42, Hershey, Pennsylvania, and SS93, Sandusky, Ohio). Cultures are available upon request.
Known distribution. USA (Georgia, Michigan, Ohio, Pennsylvania), Japan.
| DISCUSSION |
|---|
|
|
|---|
S. ingoldii was found to produce basidiospores but no gemmae in its glebae. Gemmae are assumed to play a crucial role in the lignicolous habit of Sphaerobolus because basidiospores were found not to germinate on wood (Ingold 1972
). Because most of the isolates of S. ingoldii were found on wood, this raises the question whether the basidiospores of S. ingoldii in fact require the presence of proteolytic enzymes and high temperature to stimulate germination as had been reported in S. stellatus (Dykstra 1982
) or are able to germinate on wood and other plant debris directly after discharge. It is possible that the in vitro fruiting conditions were inappropriate for stimulating gemmae production, although gemmae were produced in vitro in S. iowensis and S. stellatus under the same conditions. Another taxon that lacks gemmae in this group of fungi is Nidulariopsis melanocarpa Greis. The genus Nidulariopsis Greis is likely a synonym for Sphaerobolus as indicated by the name Nidulariopsis iowensis (Walker) Zeller, a synonym of S. iowensis and the morphological descriptions provided by Greis (1935)
. Although N. melanocarpacharacterized by basidiocarps 2.54 mm diam, yellow exterior, and a whitish basidiocarp interiorclearly is not conspecific with S. ingoldii described herein, both species lack gemmae. Future comparative investigations should address the systematic position and taxonomy of N. melanocarpa, as well as the ecological importance of the presence or absence of gemmae.
Although S. ingoldii has a wide range of geographic distribution, it seems to be more rarely encountered than the other two species. We found only five representatives of S. ingoldii in the 80 isolates analyzed. Somewhat to our surprise, we found that S. iowensis is quite common in North America, perhaps as or more common than S. stellatus, representing more than half the isolates we encountered. There are only two previous reports of S. iowensis, both from the Midwest (Walker 1927
, Flegler 1984
). While our geographic sampling admittedly is biased toward the eastern USA, we were able to obtain both S. ingoldii and S. stellatus from other continents, but not S. iowensis. It is worth noting that all three species are sympatric in the eastern USA and Midwest, often inhabiting the same microhabitat, apparently without interbreeding.
While considerable DNA polymorhism was found in S. iowensis and S. stellatus, haplotype networks, in agreement with gene genealogies detailed earlier, little biogeographic structure was indicated in both species. In numerous instances isolates from the same region appeared in divergent portions of the networks, while isolates from distant localities often clustered together. Geographic distribution of the highest level nesting clades of S. iowensis and S. stellatus are shown (FIG. 14C, D
, respectively). The null hypothesis of no geographic association of haplotypes could not be rejected in most clades, suggesting either panmixia and/or extensive dispersal or inadequate geographic sampling, possibilities that cannot be distinguished with this method. The dual coprophilous and lignicolous ecology of Sphaerobolus species provides many possible dispersal scenarios over great geographic distances. In prehistoric times animal migration might have been a considerable way of dispersal, occasionally supplemented with the natural movements of rotten wood and other plant debris. However, recent human activities might play a role in dispersal. Such activities include the transportation of commercial wood products and livestock between distant areas. The statistically significant results obtained in one clade of each species suggest restricted gene flow due to isolation by distance for some isolates. Because restricted gene flow implies only limited movement by individuals during any given generation, the geographic centers of all the clades nested together should be close, therefore the Dc and Dn values should show similar patterns under this scenario. We were interested in the fact that, even though most clades of both S. iowensis and S. stellatus did not show statistically significant genotype-geography association, the Dc and Dn values in fact were similar in most cases, supporting the hypothesis of restricted gene flow among at least some of the isolates.
In addition, while isolates of S. stellatus came from a much wider geographic area, nucleotide variation in S. stellatus was considerably smaller than that observed among S. iowensis isolates. This observation raises the question of whether S. iowensis represents a native species with genetically diverse populations and whether S. stellatus has been introduced to North America. To answer this question is beyond the scope of this paper; however, the occurrence of isolates with identical ITS sequence (haplotype 0) in distant locations (Alaska, Colorado, Maryland, Massachusetts, Minnesota, Pennsylvania and The Netherlands) is consistent with the hypothesis of introduction and spread of S. stellatus to North America.
In this paper we provided the first systematic overview of the genus Sphaerobolus based on both molecular and morphological data that can serve as base for future research on the artillery fungi. Because our analyses included isolates predominantly from North America, future studies should include samples from other areas to further elucidate the genetic diversity and phylogenetic structure of the genus.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Corresponding author. E-mail: jgeml{at}iab.alaska.edu
| LITERATURE CITED |
|---|
|
|
|---|
Alasoadura SO. 1963. Fruiting in Sphaerobolus with special reference to light. Ann Bot N.S. 27:123145.
Aplin TEH. 1961. Sphaerobolus stellatus, a new fungus for Western Australia. West Aust Nat 8:2729.
Baayen RP, ODonnell K, Breeuwsma S, Geiser DM, Waalwijk C. 2001. Molecular relationships of fungi within the Fusarium redolens-F. hostae clade. Phytopathology 91:10371044.[Medline]
Brantley EA, Davis DD, Kuhns LJ. 2001a. Biological control of the artillery fungus, Sphaerobolus stellatus, with Trichoderma harzianum and Bacillus subtilis. J Environ Hort 19:2123.
, , . 2001b. Influence of mulch characteristics on the sporulation of the artillery fungus Sphaerobolus stellatus. J Environ Hort 19:8995.
Bruns TD, White TJ, Taylor JW. 1991. Fungal Molecular Systematics. Ann Rev Ecol Syst 22:525564.[CrossRef]
Buller AHR. 1933. Researches on fungi. Vol. 5. London: Longmans, Green & Co.
Dring DM. 1964. Gasteromycetes of west tropical Africa. Commonwealth Mycological Institute Mycological Papers. p 98109.
Dykstra MJ. 1982. A cytological examination of Sphaerobolus stellatus fruiting bodies. Mycologia 74:4453.[CrossRef]
Felsenstein J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783791.[CrossRef]
Flegler SL. 1984. An improved method for production of Sphaerobolus fruit bodies in culture. Mycologia 76:944946.[CrossRef]
Fletcher M, Cooke C. 1984. Carbohydrate changes in the developing sporophore of Sphaerobolus stellatus. Trans Br Mycol Soc 82:366369.
Geml J, Davis DD, Geiser DM. 2005. Phylogenetic analyses reveal deeply divergent species lineages in the genus Sphaerobolus (Phallales: Basidiomycota). Mol Phyl Evol 35:313322.[CrossRef][Medline]
Greis H. 1935. Nidulariopsis melanocarpa Greis nov. gen. spec. und eine neue Form von Sphaerobolus iowensis. Hedwigia 75:255266.
Halgrimsson H, Jensson E, Kristinsson H. 1992. Three new gasteromycetes discovered in Iceland. Natturufraedingurinn 61:219227.
Hawksworth DL, Kirk PM, Sutton BC, Pegler DN. 1996. Ainsworth & Bisbys dictionary of the fungi. 8th ed. New York: Oxford University Press Inc.
Herrera T, Perezsilva E. 1987. Sphaerobolusrecord of a new genus of Gasteromycetes for Mexico. Mycotaxon 28:413414.
Hibbett DS, Pine EM, Langer E, Langer G, Donoghue MJ. 1997. Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences. Proc Natl Acad Sci USA 94:1200212006.
Hillis DM, Dixon MT. 1991. Ribosomal DNAmolecular evolution and phylogenetic inference. Quart Rev Biol 66:410453.
Hopple JS, Vilgalys R. 1999. Phylogenetic relationships in the mushroom genus Coprinus and dark-spored allies based on sequence data from the nuclear gene coding for the large ribosomal subunit RNA: divergent domains, outgroups, and monophyly. Mol Phyl Evol 13: 119.[CrossRef][Medline]
Huelsenbeck JP, Ronquist F. 2001. MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics 17:754755.
Ingold CT. 1971. The glebal mass of Sphaerobolus. Trans Br Mycol Soc 56:105113.
. 1972. Sphaerobolus: the story of a fungus. Trans Br Mycol Soc 58:179195.
, Peach J. 1970. Further observations on fruiting in Sphaerobolus in relation to light. Trans Br Mycol Soc 54:211220.
Kirk PM, Cannon PF, David JC, Stalpers JA. 2001. Dictionary of the Fungi. 9th ed. Wallingford, U.K.: CAB International.
Lehman RD. 1985. Black spots on housesan insect or disease problem? Penn Dept Agric Bur Pl Ind, Reg Hort 11:1516.
Lutzoni F, Vilgalys R. 1995. Integration of morphological and molecular datasets in estimating fungal phylogenies. Can J Bot Suppl 73:S649659.[CrossRef]
, Wagner W, Reeb V, Zoller S. 2000. Integrating ambiguously aligned regions of DNA sequences in phylogenetic analyses without violating positional homology. Syst Biol 49:628651.
McKenzie EHC, Foggo MN. 1989. Fungi of New Zealand subantarctic islands. New Zeal J Bot 27:91100.
Miadlikowska J, Lutzoni F, Goward T, Zoller S, Posada D. 2003. New approach to an old problem: incorporating signal from gap-rich regions of the ITS and rDNA large subunit into phylogenetic analyses to resolve the Peltigera canina species complex. Mycologia 95:11811203.
Micheli PA. 1729. Nova Plantarum Genera. Florence, Italy.
Moncalvo JM, Vilgalys R, Redhead SA, Johnson JE, James TY, Aime MC, Hofstetter V, Verduin SJW, Larsson E, Baroni TJ, Thorn RG, Jacobsson S, Clémençon H, Miller Jr OK. 2002. One hundred and seventeen clades of Euagarics. Mol Phyl Evol 23:357400.[CrossRef][Medline]
Nawaz M. 1967. Phototropism in Sphaerobolus. Biologia 13: 514.
ODonnell K. 1996. Progress toward a phylogenetic classification of Fusarium. Sydowia 48:5770.
, Kistler HC, Cigelnik E, Ploetz RC. 1998. Multiple evolutionary origins of the fungus causing Panama disease of banana: concordant evidence from nuclear and mitochondrial gene genealogies. Proc Natl Acad Sci USA 95:20442049.
, Lutzoni FM, Ward TJ, Benny GL. 2001. Evolutionary relationships among mucoralean fungi (Zygomycota): evidence for family polyphyly on a large scale. Mycologia 93:286296.[CrossRef]
Ott RL. 1993. An introduction to statistical methods and data analysis. 4th ed. Belmont, California: Duxbury Press.
Pine EM, Hibbett DS, Donoghue MJ. 1999. Phylogenetic relationships of cantharelloid and clavarioid Homobasidiomycetes based on mitochondrial and nuclear rDNA sequences. Mycologia 91:944963.[CrossRef]
Posada D, Crandall KA. 1998. Modeltest: testing the model of DNA substitution. Bioinformatics 14:817818.
Roger AJ, Sandbloom O, Doolittle WF, Philippe H. 1999. An evaluation of elongation factor 1 alpha as a phylogenetic marker for eukaryotes. Mol Biol Evol 16:218233.[Abstract]
Schussler A, Schwarzott D, Walker C. 2001. A new fungal phylum, the Glomeromycota: phylogeny and evolution. Mycol Res 105:14131421.[CrossRef]
Swofford DL. 2002. PAUP*. Phylogenetic Analysis Using Parsimony (*and other methods), 4.0b4a. Sunderland, Massachusetts: Sinauer Associates.
Taylor JW, Jacobson DJ, Kroken S, Kasuga T, Geiser DM, Hibbett DS, Fisher MC. 2000. Phylogenetic species recognition and species concepts in fungi. Fung Gen Biol 31:2132.[CrossRef]
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. 1997. The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25:48764882.
Ulloa M, Hanlin RT. 2000. Illustrated dictionary of Mycology. St Paul, Minnesota: APS Press.
Walker LB. 1927. Development and mechanism of discharge in Sphaerobolus iowensis and S. stellatus Tode. J Eli Mitch Sci Soc 42:151178.
, Anderson EN. 1925. Relation of glycogen to spore-ejection. Mycologia 15:154159.
White TJ, Bruns T, Lee S, Taylor JW. 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. New York: Academic Press Inc. p 315322.
Zervakis G, Dimou D, Ballis C. 1998. A check-list of the Greek macrofungi including hosts and biogeographic distribution: I. Basidio Mycotax 66:273336.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |