- Research article
- Open Access
Evolutionary history of Serpulaceae (Basidiomycota): molecular phylogeny, historical biogeography and evidence for a single transition of nutritional mode
BMC Evolutionary Biology volume 11, Article number: 230 (2011)
The fungal genus Serpula (Serpulaceae, Boletales) comprises several saprotrophic (brown rot) taxa, including the aggressive house-infecting dry rot fungus Serpula lacrymans. Recent phylogenetic analyses have indicated that the ectomycorrhiza forming genera Austropaxillus and Gymnopaxillus cluster within Serpula. In this study we use DNA sequence data to investigate phylogenetic relationships, historical biogeography of, and nutritional mode transitions in Serpulaceae.
Our results corroborate that the two ectomycorrhiza-forming genera, Austropaxillus and Gymnopaxillus, form a monophyletic group nested within the saprotrophic genus Serpula, and that the Serpula species S. lacrymans and S. himantioides constitute the sister group to the Austropaxillus-Gymnopaxillus clade. We found that both vicariance (Beringian) and long distance dispersal events are needed to explain the phylogeny and current distributions of taxa within Serpulaceae. Our results also show that the transition from brown rot to mycorrhiza has happened only once in a monophyletic Serpulaceae, probably between 50 and 22 million years before present.
This study supports the growing understanding that the same geographical barriers that limit plant- and animal dispersal also limit the spread of fungi, as a combination of vicariance and long distance dispersal events are needed to explain the present patterns of distribution in Serpulaceae. Our results verify the transition from brown rot to ECM within Serpulaceae between 50 and 22 MyBP.
Studying fungal distributions in an evolutionary context is relatively new, mainly due to the long-standing perception that fungi are more or less free from dispersal barriers, and that fungal distributions are primarily controlled by the distribution of hosts and substrata (see ). However, this notion of "cosmopolitan species" has recently been challenged by molecularly based studies [2–4]. Fungal distributional patterns can be complicated by the fact that cryptic species and polyphyly often are revealed in molecular based fungal studies. Thus, molecular phylogenetic approaches have proven important for the identification of unique lineages with restricted distributions . Although it is obvious that hosts impose some restriction on the distribution of their parasite or symbiont (e.g. ), Taylor et al.  argue that very few (if any) organisms can be considered free of entanglement with other organisms. In fact, Parrent et al.  suggested that high levels of host specificity increase the effect of isolation-by-distance.
Vicariance versus long distance dispersal are two alternative explanations to widely disjunct distributions, and have been identified to operate in most groups of organisms [7, 8]. It is becoming increasingly accepted that the same barriers to dispersal operate in fungi as in other groups of organisms [1, 2]. Vicariance has given explanatory power for disjunct distributions in some extomycorrhizal fungi such as pacific boletes  and truffle-formin fungi in the Hysterangiales , while dispersal best explains the distribution of the wood-rot fungi Ganoderma .
Molecular dating makes it possible to relate organismal evolution to major ecological and geological events. The combination of fossils that correlate to divergence events in organismal phylogenies with increasingly realistic models of molecular evolution has increased the accuracy of molecular dating . Reliable fossils are preferred for calibration of molecular clocks. However, dated geographic and ecological events as well as age estimates of symbiotic or parasitic partners have also been used [13–15], either in combination with fossil calibration or as the best available alternative when no reliably determined and correctly assigned fossils exist. Unfortunately, the fossil record of fungi is sparse , and several studies have attempted to date evolutionary splits of fungi using various calibration strategies (e.g. [15–20]). An additional complicating factor in molecular dating is between-lineage rate heterogeneity of molecular evolution, which has been extensively documented in fungi (see  and references therein). However, methods have been developed that account for rate heterogeneity by various sophisticated relaxed clock models (e.g. BEAST ).
In recent times, the fungal tree of life has been largely redrawn by molecular phylogenetics. This has in turn improved our understanding of how different nutritional modes have evolved [22–24]. The Agaricomycetes is a class in the fungal kingdom that contains most of the forest floor mushroom-forming fungi . Currently saprotrophy (possibly white rot) is suggested to be the plesiomorphic state of the Agaricomycetes , while brown rot and ectomycorrhiza have evolved several times independently from the ancient white rot state. Ectomycorrhizal (ECM) fungi form mutualistic interactions between plant roots and fungal mycelia. The symbiosis provides the plant with important nutrients such as phosphorous and nitrogen, and the fungi receive photosynthetic products in return [27, 28]. Mycorrhizal symbioses are ubiquitous in natural plant environments [29–31]. Saprotrophy by white rot decomposes lignocellulose and utilize all parts of the wood. Brown rot, on the other hand, only deforms the lignin to get access to the cellulose and hemicelluose, which are then efficiently depolymerized and absorbed by the fungus [32, 33].
Recent work shows that ECM is a nutritional mode that has evolved several times independently in both plants and fungi (both Basidiomycota and Ascomycota) [28, 34]. Hibbett and Matheny  concluded that there have been many independent origins of ECM associations of Agaricomycetes with both Angiosperms (at least 8) and Gymnosperms (at least 6-8). The monophyletic basidiomycete order Boletales comprises various nutritional modes, including ECM, saprotrophy and parasitism [35, 36]. Brown rot has been suggested as the plesiomorphic state for the Boletales .
The genus Serpula (Pers.) Gray in the Boletales consists of saprotrophic taxa that mainly degrade conifer substrates . Serpula species produce annual, brownish, and resupinate fruitbodies (flat structures adhered to a woody or soil substrate with a morphological cariable sexual spore-producing surface; in the case of Serupla this spore producing surface is wrinkled or merulioid) [38, 39]. Many Serpula species appear frequently on structural timbers in buildings, including S. himantioides (Fr.) P. Karst., S. incrassata (Berk. & M.A. Curtis) Donk and S. lacrymans (Wulfen) J. Schröt . Serpula lacrymans is very aggressive in buildings and produces what is commonly known as dry rot (ref). The taxonomic circumscriptions of Serpula himantioides and S. lacrymans have varied through time: they are lumped as one based on morphology  or recognized as two separate species based on mating experiments . The morphotaxon S. lacrymans includes two varieties. One of them (var. shastensis Harmsen) seems to have a natural distribution in the Cascade Range in western North America, while the other (var. lacrymans [Wulfen] J. Schröt.) has a natural distribution in northeast Asia from where it has spread to most continents and colonised buildings (see  and references therein). The morphotaxon S. himantioides includes multiple (at least five) cryptic species [44, 45], most of them with a primary affinity to southern South America and North America. Serpula incrassata seems to have a strictly North American distribution . A less studied species, Serpula puverulenta (Fr.) Bondartsev, which also appears in houses, has at least a European and probably also a North American distribution [46, 47]. Serpula similis (Berk. & Broome) Ginns is known from tropical regions in Africa and southeast Asia, where it is sampled on bamboo and hardwoods [48, 49].
Molecular phylogenetic analyses have indicated that the ECM-forming genera Austropaxillus Bresinsky & Jarosch and Gymnopaxillus E. Horak are nested within Serpula, and the family Serpulaceae Jarosch & Bresinsky was described comprising these three genera . Austropaxillus species have fruiting bodies with a stipe (stem) and pileus (cap), and a bear a gilled hymenophore (spore-producing surface). In contrast, Gymnopaxillus encompasses truffle-like hypogeous (under ground) fruit bodies . Due to morphological and ecological similarities, these genera were earlier included in Paxillus Fr. . Austropaxillus and Gymnopaxillus form ECM with roots of Nothofagus and Eucalyptus trees and are known only from the temperate Southern Hemisphere.
The aim of the present study was threefold. First, we wanted to further investigate phylogenetic relationships within Serpulaceae (previously described in [37, 52]), in particular the relationship between the saprotrophic Serpula and the mycorrhiza-forming Austropaxillus and Gymnopaxillus. Secondly, we wanted to relate major Serpulaceae divergence events and ecological transitions to geologic and climatic conditions of the past. Finally, we wanted to reveal how many times the transition from brown rot to mycorrhizal symbiosis has occurred in Serpulaceae including reconstructing its ancestral biographical range within a historical geological context. We addressed these issues using a molecular phylogenetic approach, including Bayesian age estimation and various ancestral area analyses, on a multi-locus dataset.
We generated DNA sequence data from 30 accessions of five Serpula species and six accessions representing six species of Austropaxillus (Additional file 1: Specimens included in this study). In the phylogenetic analysis, we included an extensive set of outgroups from both the Boletales (41 accessions) and its sister group the Atheliales, the latter used for rooting purposes (Dataset1). For the purpose of dating, we reduced the boletalean outgroup to 24 accessions, expanded the athelialean outgroup to seven accessions, and included 1-24 accessions from each of nine additional fungal orders (Dataset2). Sequences not generated for the present study were taken from Binder et al.  and Jarosch and Bresinsky . We used a total of 79 accessions for the phylogenetic analyses (Dataset1) and 109 accessions in the dating analysis (Dataset2). A separate phylogenetic analysis was performed on a reduced dataset of 39 accessions (Serpulaceae and a single outgroup taxon; Dataset3) to depict the phylogenetic position of Gymnopaxillus, for which only sequences of a single genetic region are available. All included accessions are listed in Additional file 1.
DNA extraction, PCR amplification and DNA sequencing
A small amount of fungal tissue was homogenized on a Mixer Mill (MM301, Retsch GmbH & Co., Haan, Germany) before we extracted total DNA following the 2% CTAB miniprep method described by Murray and Thompson  with minor modifications from Gardes and Bruns . We dissolved the dried DNA pellet in 100 μL milli-Q H2O, and used further dilutions for the molecular work. We have deposited the DNA aliquots in the fungal DNA collection at the Institute of Biology, UoO, Norway.
Five nuclear DNA regions were amplified: the nuclear ribosomal large subunit (LSU) using the primers LR0R/LR5 , the nuclear ribosomal small subunit (SSU) using the primers PNS1/NS41 , the nuclear ribosomal 5.8S using the primers ITS5 and ITS4 , a fragment of the gene that encodes the second largest subunit of the RNA polymerase II (rpb2) using the primers b3.1F and b6R2 , and the translation elongation factor 1α (tef 1) using the primers EF595F and EF1160R .
We generated amplicons using the following PCR protocol: 4 min at 94°C; followed by 35 cycles of 25 s at 94°C, 30 s at 56°C (rpb2 53°C; ITS and tef 154°C), and 60 s (30 s ITS and tef 1) at 72°C; followed by a 10 min extension at 72°C and an indefinite hold at 4°C. For all PCR reactions we used the PuReTaq Ready-To-Go™ PCR Beads (GEhealthcare, Waukesha, WI) in a 25 μl reaction on an Eppendorf thermocycler (Mastercycler, Hamburg, Germany). We purified the PCR-products using 2 μl 10 times diluted ExoSAP-IT (GEhealthcare) to 8 μl PCR product, incubated at 37°C for 45 minutes followed by 15 minutes at 80°C. Cycle sequencing was performed by the CEES ABI-laboratory http://www.mn.uio.no/bio/english/research/about/infrastructure/abi-lab/index.html using the ABI BigDye Terminator sequencing buffer and v3.1 Cycle Sequencing kit (Life Technologies, Carlsbad, CA). Sequences were processed on an ABI 3730 DNA analyser (Life Technologies). Sequences were assembled and edited using BioEdit 7 . All sequences have been deposited in GenBank, and accession numbers are given in Additional file 1.
We aligned the sequences manually using BioEdit 7. We assembled three datasets: Dataset1, consisting of a 2856 basepair (bp) long alignment that includes five genetic regions (nuclear ribosomal LSU, SSU and 5.8S, rpb2 and tef 1), available as Additional file 2; Dataset2, consisting of a 6172 bp long alignment including longer stretches of the above-mentioned markers, available as Additional file 3; Dataset3, consisting of a 918 bp long alignment including LSU only and a total of 38 Serpulaceae accessions (incl. Gymnopaxillus) and Bondarcevomyces for rooting, available as Additional file 4. Due to restricted sequence data, we did not include Gymnopaxillus in the final analyses, but Dataset3 was analyzed to confirm the phylogenetic position of Gymnopaxillus within Serpulaceae and together with Austropaxillus.
Dataset1 and 2 were analyzed using both Maximum Parsimony (MP) and Bayesian Inference (BI) phylogentic methods. Dataset3 was analyzed using MP only. We performed preliminary MP phylogenetic analyses on each DNA region separately to control visually for between-region incongruence prior to concatenation and final phylogenetic analyses. For the MP analyses of the concatenated datasets and Dataset3 we used TNT . We performed heuristic searches with 2000 random addition sequences and TBR branch swapping, saving ten trees per replication. The resulting trees were swapped with TBR saving up to 10 000 trees. We set the collapsing rule to minimum length = 0 and random seed was set to "time". Jackknife  resampling was performed with 2000 replicates (10 random entry orders and 10 trees saved each repetition), 36% deletion, absolute frequencies as output, and a cut-off value of 50%.
The best-fit evolutionary models were estimated using the Akaike information criterion (AIC) and the software MrModeltest 2.3  at the Bioportal http://www.bioportal.uio.no. We divided the datasets (Dataset1 and Dataset2) into three partitions: (1) the ribosomal regions (SSU, LSU and 5.8S; collectively nrDNA), (2) rpb2, and (3) tef 1. The best-fit models for these three parts were: GTR+I+G for the ribosomal regions, HKY+I+G for rpb2, and SYM+I+G for tef 1.
For the BI analyses, we analyzed a partitioned (according to the above established parts) concatenated alignment of all five nuclear regions using MrBayes 3.1.2 [65, 66] at the Bioportal. Bayesian posterior probabilities were determined twice by running one cold and three heated chains for six million generations, saving trees every 1000th generation. To test whether the Markov Chain converged, we monitored the standard deviation of split frequencies (SDSF), and discarded as burn-in the generations prior to the point when the SDSF fell below 0.01 when comparing two independent runs. A 50% majority rule consensus tree was used to calculate posterior probabilities.
Divergence time estimation
We used Dataset2 (five DNA regions and 109 accessions) and the computer program Bayesian Evolutionary Analysis Sampling Trees (BEAST) 1.6.1  for estimating the posterior probability distribution of divergence times in Serpulaceae. BEAST co-estimates phylogeny and divergence times . Details about the Bayesian approach used in BEAST have been described by others .
Xml-files for the BEAST analyses were constructed using BEAUti 1.6.1 (BEAST package). We analyzed the data under four different treatments; different clock models (strict versus uncorrelated lognormal relaxed [UCLD]), and different prior distributions for the calibrations (exponential versus lognormal). The xml-files of uncorrelated lognormal relaxed clock with exponential and lognormal prior distributions are available as Additional file 5 and Additional file 6, respectively. Because among branch rate heterogeneity within the data was high, we also performed estimations without long-branch taxa (indicated by arrows in Figure 1). However, as divergence time estimates from the analyses with and without long-branch taxa were highly similar, only results from the analysis including the long-branch taxa is presented (see Results). We partitioned the data according to marker and codon position (only the two coding genes; [1+2],3). Partitioning of data into markers and appropriate codon usage has been shown to produce shorter node heights than unpartitioned analyses . A Yule tree prior was employed in all runs, which assumes a constant speciation rate for each branch in the tree. A lognormal distribution was used for the priors for all calibrated nodes because it fixes the minimum age and allows the maximum age to be sampled following a lognormal distribution with no hard upper limit.
Since there are no known fossils of Serpulaceae we used an extensive outgroup for which fossils are available. We included three primary calibration points (corresponding to the numbered nodes in Figure 2): (1) Quatsinoporites cranhamii Smith, Currah and Stockey from the Early Creatceous (Barremian) ; (2) Archeomarasmius leggetti Hibbett, Grimaldi and Donoghue from the early Late Cretaceous [70, 71]; and, (3) a permineralized "sulloid" ECM associated with Pinus roots of at least 50 million years . We used the following lognormal settings (Mean in Real Space): (1) logmean = 2.0, logstdev = 0.5, and offset = 125.0; (2) logmean = 2.5, logstdev = 0.5, and offset = 90.0; and, (3) logmean = 1.5, logstdev = 0.5, and offset = 50.0. The exponential settings were: (1) mean = 3.0 and offset = 130.0; (2) mean = 3.2 and offset = 90.0; and, (3) mean = 2.0 and offset = 50.0. These lognormal/exponential settings correspond approximately to the following median and 95% HPD (square-bracketed; numbers rounded off to nearest million): (1) 127 [126-129]/132 [130-139]; (2) 92 (2) 92 [91-95]/92 [90-99]; and (3) 51 [50-53]/51 [50-56]. BEAST may get stuck on local optima because it does not employ a coupled MCMC. Thus, ten additional monophyletic sets of taxa (indicated with black dots on Figure 2) were predefined in order to obtain a topology that corresponds to the phylogeny obtained from a previous study  that used a more robust phylogenetic method (e.g. MrBayes) and more taxa. The Markov chains were run for ten million generations, sampling and saving every 1000th tree. Convergence of the chain to stationary distributions was controlled by inspection of the MCMC samples (Trace) in each analysis using the program Tracer 1.5 (BEAST package), and 1000 trees (i.e. 10%) were discarded as burn-in prior to summarizing the posterior distribution of trees and indentifying the maximum clade credibility (MCC) tree using TreeAnnotator 1.6.1 (BEAST package). Each treatment was analyzed five times to determine if the independent runs converged on the same posterior distribution. In order to obtain effective sample sizes (ESS) of more than 100 for all parameters, we combined the log output files from separate runs with the same treatment using LogCombiner 1.6.1 (BEAST package). To compare the different clock treatments, we monitored the "CoefficientOfVariation" for each partition in the relaxed clock analyses, and discarded a strict clock evolution if the distribution did not bump up against zero.
Ancestral area analyses
We predefined unit areas based on present-day distributions of Serpulaceae taxa. The following seven unit areas were used: (A) western North America; (B) eastern North America; (C) western Eurasia; (D) eastern Eurasia; (E) southern South America; (F) Australia + New Zealand; and, (G) Old World Tropical. For each terminal taxon we coded the unit area(s) in which it presently occurs. The variety Serpula lacrymans var. shastensis has been recorded once in the Caucasus , outside its range as defined in this study. However, this observation has not been confirmed by other studies.
We used two alternative event-based ancestral area methods: (1) a Bayesian approach to dispersal-vicariance analysis (DIVA ) as implemented in the computer software Reconstruct Ancestral States in Phylogenies (RASP ) and, (2) the dispersal-extinction-cladogenesis analysis (DEC ) implemented in the computer program LAGRANGE. The Bayesian approach in RASP is built on the source code of MrBayes 3.1.2, and is a generalization of the methods used in Olsson et al.  and Sanmartin et al. . We chose the F81 model for the Bayesian MCMC analyses, allowing for different rates of change among ancestral areas. The Bayesian posterior probabilities were determined twice by running 10 chains fifty thousand generations, saving reconstructions every 100th generations.
The ancestral area analyses were conducted on the posterior distribution of dated Serpulaceae trees estimated from BEAST. RASP analyses were carried out using 103 randomly selected trees from BEAST.
The DEC analyses in LAGRANGE enables maximum likelihood estimation of range inheritance scenarios at cladogenesis events by modelling transitions between discrete states (ranges) along phylogenetic branches as a function of time. The DEC script was produced using a web-based analysis configuration tool available from http://www.reelab.net/lagrange/configurator/index.
In the MP analysis of Dataset1, we obtained 469 most parsimonious trees (MPTs) of 4485 steps and with a reconciled consistency index (RC) of 0.2 and a homoplasy index (HI) of 0.69. Our results from the BI and MP analyses were largely congruent, but variously resolved. The Bayesian 50% majority rule consensus tree is presented with branch support (posterior probabilities [PP] and parsimony jackknifing [JK]) superimposed in Figure 1. Results render Serpulaceae monophyletic with Tapinellineae as phylogenetic sister (Figure 1). The Serpulaceae/Tapinellineae group consists of brown rot species except for the Austropaxillus clade. The brown rot species belonging to Coniophora and Coniophoraceae appear as a monophyletic group separated from the Serpulaceae/Tapinellineae clade (Figure 1).
The saprotrophic Serpula himantioides and S. lacrymans appear as sister taxa in the phylogeny and together form a monophyletic group (Figure 1). Serpula lacrymans is further subdivided in two lineages representing the two described varieties, var. lacrymans and var. shastensis (see ). Austropaxillus forms a monophyletic sister to the S. lacrymans/S. himantioides group that makes Serpula paraphyletic as currently circumscribed (Figure 1). A clade comprising the saprotrophic species Serpula incrassata, S. pulverulenta and S. similis is recovered as sister to all other Serpulaceae (Figure 1). The branch leading to S. similis is very long, indicating an increased rate of molecular evolution in this lineage.
The MP analysis of Dataset3 produced 16 MTPs of 311 steps, an RC of 0.90, and an HI of 0.06. The strict consensus tree with JK values superimposed is presented in the Additional files. A monophyletic Gymnopaxillus group renders Austropaxillus paraphyletic (Additional file 7, Figure. S1).
Divergence time estimation
The posterior distribution of the "CoefficientOfVariation" suggests rejection of a strict clock for all gene partitions. Relaxed (UCLD) clock analyses with lognormal versus exponential priors for the calibrations estimated rather similar age estimates of Serpulaceae divergences (Table 1). The maximum clade credibility (MCC) tree obtained from the relaxed clock/lognormal treatment is presented in Figure 2. A subsection of this chronogram, including Serpulaceae only, is presented in Figure 3. Mean nodal age estimates for numbered nodes in Figure 3 and their respective 95% highest posterior density (HPD) intervals are reported in Table 1.
According to the estimated divergence times, Serpulaceae diverged from the Tapinellineae during the Late Cretaceous, and the mean crown-group age of Serpulaceae was estimated to about 64 (80-58) MyBP (Figure 2; Table 1; Figure 3, clade 1). The three Serpula species, S. incrassata, S. similis and S. pulverulenta (Figure 3, clade 2) appear to be the oldest species within Serpulaceae, with divergences among sampled species evident from the Late Cretaceous to Early Eocene (71-48 MyBP).
The divergence between the ECM group and the brown rot species (Austropaxillus vs. S. lacrymans and S. himatioides) is estimated to have happened between the Late Cretaceous to Late Eocene (66-35 MyBP; Figure 3, node 5). The posterior mean crown-group age of Austropaxillus is approximately 22 (32-12) My (Figure 3, node 6). Subsequent speciation events of sampled Austropaxillus species seem to suggest a fairly constant rate of speciation. The divergence of the two widely distributed Serpula lacrymans and S. himantioides is estimated to have occurred in Miocene; the mean age estimate being 16 MyBP (Figure 3, node 9). Intraspecific diversification leading to varieties and cryptic species within S. lacrymans and S. himantioides is estimated to have occurred during the Miocene and Pliocene (Figure 3).
Ancestral area analyses
The inferred historical biogeographic scenarios from analyses using DEC (from LAGRANGE) and RASP are summarized in Table 1 and Figure 3. The inferred ancestral areas at internal nodes estimated using the Bayesian RASP correspond largely to the results obtained from the maximum likelihood estimates in DEC. RASP estimated more combined ancestral areas than DEC, and the RASP estimates have generally a higher relative probability for each ancestral area (Table 1).
The maximum likelihood reconstruction of ancestral areas for the internal nodes of Serpulaceae show that the most recent common ancestor (tmrca) of clades 2-4 (Figure 3, node 5) most likely appeared in western North America during the Paleocene-Eocene boundary. Range expansion of one daughter lineage, Austropaxillus (Figure 3, clade 2), into South America, Australia and/or New Zealand some time during the Eocene to Mid Miocene followed this. The divergence into one mainly southern South American clade (node 8) and one Australian/New Zealand clade (node 7) is estimated to have happened about 15 (24-7) MyBP. The second daughter lineage, tmrca of S. himantioides and S. lacrymans (Figure 3, clades 3-4, node 9), most likely occurred in western North America during the Early Eocene to Mid Miocene.
Molecular phylogenetics of Serpulaceae
Our results based on five nuclear DNA regions corroborate the monophyly of Serpulaceae sensu Jarosch . Moreover, a close relationship between the saprotrophic Serpula and the ectomycorrhizal Austropaxillus is supported, as already indicated in previous studies based on more restricted taxon samplings and fewer genetic markers [37, 50]. Gymnopaxillus nested within the Australian Austropaxillus clade (Additional file 7, Figure. S1) in the nrLSU analysis. Traditionally, Serpula has been recognised as a member of Coniophoraceae [39, 78]. Our molecular results clearly show that Serpula is only distantly related to Coniophora, supporting the treatment of Serpulaceae  as a separate family that comprises the three genera Austropaxillus, Gymnopaxillus and Serpula. That Austropaxillus and Gymnopaxillus are nested within Serpula suggests taxonomic changes are needed to accommodate this monophyly. Unfortunately we were unable to include samples from South American Gymnopaxillus species, so our results call for further analyses including more markers and samples from this group.
Historical biogeography of Serpulaceae
Our results indicate a Late Cretaceous origin of extant Serpulaceae, between 94 and 74 MyBP (Figure 2) in temperate North America (Figure 3). The Gymnosperm genus Pinus is the main substrate of Serpulaceae in the Northern Hemisphere. Pinus, is estimated to have evolved during the Cretaceous, between 155 and 87 MyBP . Thus the age of Pinus corresponds well with the timing of the evolution of Serpulaceae on Pinus wood. Moreover, the fossil record supports the presence of Pinaceae members in the high-latitude and high-altitude regions of North America during the early Cenozoic . This supports our result, as the main host of Serpulaceae was available at the time of origin.
The initial divergence of extant Austropaxillus taxa into one mainly southern South American clade and one Australian/New Zealand clade about 22 (32-12) MyBP (Figure 3, node 6) must represent a long distance dispersal event, and is consistent with a biogeographical tract referred to by Crisp et al as a pan-temperate tract . Nor does the relatively recent split between northern and southern hemisphere taxa in Serpulaceae support an ancient Gondwanan distribution of Austropaxillus. However, there have been fierce discussions of the time and climatic conditions of the Gondwanan break-up . Land bridges have been suggested between Antarctica and Australia up until 28 MyBP [82, 83]. During this period Antarctica was probably glaciated  and may have been too hostile for temperate fungi. Thus, even if we cannot rule out vicariance, we consider long distance dispersal a more plausible explanation. Moreover, the divergence times estimated within Serpulaceae can be expected to be somewhat younger due to the node density effect , supporting a long-distance dispersal event.
The divergence between Austropaxillus macnabbi from New Zealand and A. muelleri from Tasmania (Figure 3, node 7) must have resulted from a more recent dispersal event. New Zealand was separated from South America and Australia about 80 MyBP , i.e. about 60 My years before the split at node 7 (24-7 MyBP). Interestingly, the same recent divergence between species from New Zealand and Australia have also been found for Nothofagus , the host of Austropaxillus, and its fungal parasite Cyttaria . Long distance dispersal is also suggested in order to explain the divergence between Austropaxillus sp. from Tasmania and A. statuum from southern South America (Figure 3). However, as the sister-relationship between the two species is weakly supported, more data is needed to confirm or reject this long distance dispersal hypothesis. Additional fungal disjunctions that are too young to be caused by Gondwanan breakup, and thus are explained by long-distance dispersal have been reported elsewhere (e.g. [10, 11, 18]).
Assuming western North America as the true ancestral area at node 5 (Figure 3), the presence of southern South American species in both clade 2 (Austropaxillus) and clade 4 (S. himantioides) must be explained by independent long-distance dispersal events. This can be concluded, because Laurasia (incl. North America) and Gondwana (incl. southern South America) were separated by the east-west Tethyan seaway from the latest Jurassic (c 150 MyBP) to the end of Pliocene (c 3 MyBP; ). This seaway represented a major dispersal barrier to most terrestrial organisms. Not until the establishment of the Great American Interchange (GAI; ca. 3 MyBP) was there land connection between North and South America (Figure 3). Although the two dispersals happened at different times, none correspond to time periods of land-connection between North and South America. Moreover, the large distributional gap between North American and southern South American sister-taxa of Serpula (see inset map in Figure 3) favours long distance dispersals to historically continuous distributions with subsequent vicariance. Still, the distinctness of the disjunct lineages indicates that such inter-area dispersals are not frequent. Nevertheless, historical long distance dispersal has been reported in several mycogeographic studies (e.g. [11, 18, 88]), and dispersal between northern and southern temperate regions of the world seems to have occurred in other fungal groups also (e.g. [10, 20]) and plants .
All taxa represented in clade 3 (Figure 3) belong to the same morphospecies, Serpula lacrymans. The chronogram shows that the divergence of the two varieties occurred about 9 (13-4) MyBP (Figure 3, node 10). The presence of western North American and eastern Eurasian phylogenetic sister-taxa are indicative of a trans-Beringian distribution of their most recent common ancestor with subsequent vicariance (Figure 3). About 14-3.5 MyBP, there was continuous boreal forests across Beringia (Figure 3, BLB2; ). It seems likely therefore, that S. lacrymans had a continuous distribution from northern North America and into Eurasia and that var. lacrymans became differentiated from var. shastensis due to Beringan vicariance (Figure 3). Comparison of ancient permafrost-preserved (plant-parasitic) fungal DNA from north-eastern Siberia to modern tundra macrofungi indicate a shift in fungal diversity following the last ice-age , a scenario probably reiterated many times during the Pleistocene producing allopatric populations. Migration across a North Atlantic landbridge (NALB) has sometimes been used as an alternative explanation to Nearctic-Palearctic disjunctions, but this bridge disappeared by Mid Eocene (Figure 3; ) -- i.e. about 30 My before the estimated divergence of the two varieties. Historical Beringian distributions with subsequent vicariance and range expansion southwards have been shown also in Amanita  and Tricholoma , and we suspect that additional Beringian disjunctions would have been identified if the Holarctic was not treated as a single area in mycogeography (see e.g. [11, 20]).
Evolution of nutritional modes in Serpulaceae
According to our results, the transition to ECM has happened only once in Serpulaceae, in the lineage leading to Austropaxillus (Figure 1). Transition from a saprotrophic to ectomycorrhizal nutritional mode is a common ecological transition . The common ancestor of Austropaxillus and Gymnopaxillus diverged from the saprotrophic clade about 50 MyBP. The temperature decline and a drier climate in the end of the Eocene and the onset of Oligocene may have promoted this transition from saprotrophy to mycorrhiza. When growing in nature, S. lacrymans mainly occurs on substrates in close contact with the ground (H. Kauserud, personal observations), and it seems to be a geophilous fungus that may transport water from the ground and into the substrate. Whether such an adaptation could be a first step towards mycorrhizal symbiosis is a subject for future investigation.
Our results indicate that the mycorrhizal Austropaxillus/Gymnopaxillus lineage diverged from the saprotrophic S. lacrymans/S. himantioides group about 50 (66-35) MyBP and that the radiation of extant Austropaxillus/Gymnopaxillus taxa commenced about 22 (32-12) MyBP (Figure 3). A more recent split between Austropaxillus and S. lacrymans (53.1-15 MyBp) was suggested in a recent paper analyzing a smaller dataset of genome-sequenced fungi . The differences between these studies can probably largely be explained by taxon sampling. Because of the focus of this study, we have included more samples from the Serpulaceae compared to other fungal families; this could give a somewhat older age according to the node density effect  as argued above. Fossilized ECM demonstrates that ECM associations had evolved at least 50 MyBP, probably with a Pinus-host . It has been suggested that the radiation of ECM fungi happened as the obligate ECM hosts (Pinaceae and Fagales) became dominant in temperate forests as a consequence of drying and cooling from the Late Eocene [16, 20].
Our analyses support the monophyly of Serpulaceae sensu Jarosch and a Late Cretaceous origin of the family. A single transition from brown rot to ECM occurred within the Serpulaceae. The dating analysis indicates that this transition happened between 50 and 22 MyBP. Most Serpulaceae divergences seem to have happened within continents (Australia, South American, and North America). A combination of long distance dispersals for the southern hemisphere taxa and vicariance for the northern hemisphere taxa is used to explain the distribution of extant Serpulaceae taxa. The dating analysis indicates that the divergence of the two varieties of Serpula lacrymans occurred between 13 and 4 MyBP due to Beringian vicariance.
Peay KG, Bidartondo MI, Arnold AE: Not every fungus is everywhere: scaling to the biogeography of fungal-plant interactions across roots, shoots and ecosystems. New Phytologist. 2010, 185: 878-882. 10.1111/j.1469-8137.2009.03158.x.
Lumbsch TH, Buchanan PK, May TW, Mueller GM: Phylogeography and biogeography of fungi. Mycological Research. 2008, 112: 423-424. 10.1016/j.mycres.2008.02.002.
Taylor JW, Turner E, Townsend JP, Dettman JR, Jacobson D: Eukaryotic microbes, species recognition and the geographic limitation of species: examples from the kingdom fungi. Proceedings of the Royal Society B: Biological Sciences. 2006, 361: 1947-1963. 10.1098/rstb.2006.1923.
O'Donnell K, Rooney AP, Mills GL, Kuo M, Weber NS: Phylogeny and historical biogeography of true morels (Morchella) reveals an early Cretaceous origin and high continental endemism and provincalism in the Holarctic. Fungal Genetics and Biology. 2011, 48: 252-265. 10.1016/j.fgb.2010.09.006.
Martin F, Díez J, Dell B, Delaruelle C: Phylogeography of the ectomycorrhizal Pisolithus species as inferred from nuclear ribosomal DNA ITS sequences. New Phytologist. 2002, 153: 345-357. 10.1046/j.0028-646X.2001.00313.x.
Parrent JL, Garbelotto M, Gilbert GS: Population genetic structure of the polypore Datronia caperata in fragmented mangrove forests. Mycological Research. 2004, 108: 403-410. 10.1017/S0953756204009773.
Wiley EO: Vicariance biogeography. Annual Review of Ecology and Systematics. 1988, 19: 513-542. 10.1146/annurev.es.19.110188.002501.
de Queiroz A: The resurrection of oceanic dispersal in historical biogeography. Trends in ecology & evolution. 2005, 20: 68-73. 10.1016/j.tree.2004.11.006.
Halling RE, Osmundson TW, Neves M-A: Pacific boletes: Implications for biogeographic relationships. Mycological Research. 2008, 112: 437-447. 10.1016/j.mycres.2007.11.021.
Hosaka K, Castellano MA, Spatafora JW: Biogeography of Hysterangiales (Phallomycetidae, Basidiomycota). Mycological Research. 2008, 112: 448-462. 10.1016/j.mycres.2007.06.004.
Moncalvo JM, Buchanan PK: Molecular evidence for long distance dispersal across the Southern Hemisphere in the Ganoderma applanatum-australe species complex (Basidiomycota). Mycological Research. 2008, 112: 425-436. 10.1016/j.mycres.2007.12.001.
Berbee ML, Taylor JW: Dating the molecular clock in fungi-how close are we?. Fungal Biology Reviews. 2010, 24: 1-16. 10.1016/j.fbr.2010.03.001.
Allwood J, Gleeson D, Mayer G, Daniels S, Beggs JR, Buckley TR: Support for vicariant origins of the New Zealand Onychophora. Journal of Biogeography. 2010, 37: 669-681. 10.1111/j.1365-2699.2009.02233.x.
Benavides E, Baum R, Snell HM, Snell HL, Sites JW: Island biogeography of Galapagos lava lizards (Tropiduridae: Microlophus): species diversity and colonization of the archipelago. Evolution. 2009, 63: 1606-1626. 10.1111/j.1558-5646.2009.00617.x.
Berbee ML, Taylor JW: Dating the evolutionary radiations of the true fungi. Canadian Journal of Botany. 1993, 71: 1114-1127. 10.1139/b93-131.
Bruns TD, Szaro TM, Gardes M, Cullings KW, Pan JJ, Taylor DL, Horton TR, Kretzer A, Garbelotto M, Li Y: A sequence database for the identification of ectomycorrhizal basidiomycetes by phylogenetic analysis. Molecular Ecology. 1998, 7: 257-272. 10.1046/j.1365-294X.1998.00337.x.
Geml J, Geiser DM, Royse DJ: Molecular evolution of Agaricus species based on ITS and LSU rDNA sequences. Mycological Progress. 2004, 3: 157-176. 10.1007/s11557-006-0086-8.
Hibbett D: Shiitake mushrooms and molecular clocks: historical biogeography of Lentinula. Journal of Biogeography. 2001, 28: 231-241.
Jeandroz S, Murat C, Wang Y, Bonfante P, Le Tacon F: Molecular phylogeny and historical biogeography of the genus Tuber, the "true truffles". Journal of Biogeography. 2008, 35: 815-829. 10.1111/j.1365-2699.2007.01851.x.
Matheny PB, Aime MC, Bougher NL, Buyck B, Desjardin DE, Horak E, Kropp BR, Lodge DJ, Soytong K, Trappe JM, Hibbett DS: Out of the Palaeotropics? Historical biogeography and diversification of the cosmopolitan ectomycorrhizal mushroom family Inocybaceae. Journal of Biogeography. 2009, 36: 577-592. 10.1111/j.1365-2699.2008.02055.x.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 2007, 7: 214-10.1186/1471-2148-7-214.
Hibbett D: A Phylogenetic overview of the Agaricomycotina. Mycologia. 2006, 98: 917-925. 10.3852/mycologia.98.6.917.
James TY, Kauff F, Schoch CL, Matheny PB, Hofstetter V, Cox CJ, Celio G, Gueidan C, Fraker E, Miadlikowska J, Lumbsch TH, Rauhut A, Reeb V, Arnold EA, Amtoft A, Stajich JE, Hosaka K, Sung G-H, Johnson D, O'Rourke B, Crockett M, Binder M, Curtis JM, Slot JC, Wang Z, Wilson AW, Schüler A, Longcore JE, O'Donnell K, Mozley-Standridge S, et al: Reconstructing the early evolution of Fungi using a six-gene phylogeny. Nature. 2006, 443: 818-822. 10.1038/nature05110.
Binder M, Hibbett DS: Higher-level phylogenetic relationships of Homobasidiomycetes (mushroom-forming fungi) inferred from four rDNA regions. Molecular Phylogenetics and Evolution. 2002, 22: 76-90. 10.1006/mpev.2001.1043.
Hibbett DS: After the gold rush, or before the flood? Evolutionary morphology of mushroom-forming fungi (Agaricomycetes) in the early 21st century. Mycological Research. 2007, 111: 1001-1018. 10.1016/j.mycres.2007.01.012.
Hibbett D, Matheny PB: The relative ages of ectomycorrhizal mushrooms and their plant hosts estimated using Bayesian relaxed molecular clock analyses. BMC Biology. 2009, 7: 13-10.1186/1741-7007-7-13.
Smith SE, Read DJ: Mycorrhizal symbiosis. 2008, Oxford: Elsevier
Tedersoo L, May TW, Smith ME: Ectomycorrhizal lifestyle in fungi: global diversity, distribution, and evolution of phylogenetic lineages. Mycorrhiza. 2010, 20: 217-263. 10.1007/s00572-009-0274-x.
Read DJ, Perez-Moreno J: Mycorrhizas and nutrient cycling in ecosystems-a journey towards relevance?. New Phytologist. 2003, 157: 475-492. 10.1046/j.1469-8137.2003.00704.x.
Selosse MA, Le Tacon F: The land flora: a phototroph-fungus partnership?. Trends in Ecology & Evolution. 1998, 13: 15-20. 10.1016/S0169-5347(97)01230-5.
Brundrett M: Coevolution of roots and mycorrhizas of land plants. New Phytol. 2002, 154: 275-304. 10.1046/j.1469-8137.2002.00397.x.
Boddy L, Frankland JC, van West P, eds: Ecology of saprotrophic Basidiomycetes. 2008, Oxford: Elsevier
Hibbett DS, Donoghue MJ: Analysis of character correlations among wood decay mechanisms, mating systems, and substrate ranges in Homobasidiomycetes. Systematic Biology. 2001, 50: 215-242.
Hibbett D, Gilbert L-B, Donoghue M: Evolutionary instability of ectomycorrhizal symbioses in basidiomycetes. Nature. 2000, 407: 506-508. 10.1038/35035065.
Binder M, Hibbett DS, Larsson K-H, Langer E, Langer G: The phylogenetic distribution of resupinate forms across the major clades of mushroom-forming fungi (Homobasiodiomycetes). Systematic Biodiversity. 2005, 3: 113-157. 10.1017/S1477200005001623.
Hibbett D, Binder M: Evolution of complex fruiting-body morphologies in homobasidiomycetes. Proceedings of the Royal Society B: Biological Sciences. 2002, 269: 1963-1969. 10.1098/rspb.2002.2123.
Binder M, Hibbett DS: Molecular systematics and biological diversification of Boletales. Mycologia. 2006, 98: 971-981. 10.3852/mycologia.98.6.971.
Falck R: Die Meruliusfäule des Bauholzes. Neue Untersuchungen über Unterscheidung, Verbreitung, Entstehung und Bekämpfung des echten Hausscwammes. 1912, Jena: G. Fischer, 6:
Hallenberg N, Eriksson J: The Lachnocladiaceae and Coniophoraceae of North Europe. 1985, Oslo: Fungiflora A/S
Bech-Andersen J: The dry rot fungus and other fungi in houses. 1995, Holte, Denmark: Hussvamp Laboratoriet
Cooke WB: The genera Serpula and Meruliporia. Mycologia. 1957, 49: 197-225. 10.2307/3755630.
Harmsen L, Bakshi BK, Choudhury TG: Relationship between Merulius lacrymans and M. himantioides. Nature. 1958, 181: 1011-
Kauserud H, Svegården IB, Saetre G-P, Knudsen H, Stensrud Ø, Schmidt O, Doi S, Sugiyama T, Högberg N: Asian origin and rapid global spread of the destructive dry rot fungus Serpula lacrymans. Molecular Ecology. 2007, 16: 3350-3360. 10.1111/j.1365-294X.2007.03387.x.
Carlsen T, Engh IB, Decock C, Rajchenberg M, Kauserud H: Multiple cryptic species with divergent substrate affinities in the Serpula himantioides species complex. Fungal Biology. 2011, 115: 54-61. 10.1016/j.funbio.2010.10.004.
Kauserud H, Stensrud Ø, Decock C, Shalchian-Tabrizi K, Schumacher T: Multiple gene genealogies and AFLPs suggest cryptic speciation and long-distance dispersal in the basidiomycete Serpula himantioides (Boletales). Molecular Ecology. 2006, 15: 421-431.
Harmsen L: De danske Merulius-arter. Botanisk Tidsskrift. 1954, 50: 146-162.
Christiansen MP: Danish resupinate fungi. Part II: Homobasidiomycetes. Dansk Botanisk Arkiv. 1960, 19: 57-388.
Carlier F-X, Bitew A, Castillo G, Decock C: Some Coniophoraceae (Basidiomycetes, Boletales) from the Ethiopian highlands: Coniophora bimacrospora, sp. nov. and a note on the phylogenetic relationships of Serpula similis and Gyrodontium. Cryptogamie, Mycologie. 2004, 25: 261-275.
Ginns J: The genus Merulius IV. Species proposed by Berkeley and Curtis, and by Berkeley and Broome. Mycologia. 1971, 63: 219-236. 10.2307/3757757.
Jarosch M: Zur Molekularen Systematik der Boletales: Coniophorineae, Paxillineae und Suillineae. 2001, Berlin-Stuttgart: J. Cramer, 191:
Claridge A, Trappe J, Castellano M: Australasian truffle-like fungi. X. Gymnopaxillus (Basidiomycota, Austropaxillaceae). Australian Systematic Botany. 2001, 14: 273-281. 10.1071/SB00012.
Bresinsky A, Jarosch M, Fischer M, Schönberger I, Wittmann-Bresinsky B: Phylogenetic relationships within Paxillus s.l. (Basidiomycetes, Boletales): Separation of a Southern Hemisphere Genus. Plant Biology. 1999, 1: 327-333. 10.1111/j.1438-8677.1999.tb00260.x.
Binder M, Larsson K-H, Matheny B, Hibbett DS: New orders of Agaricomycetidae Amylocorticiales ord. nov. and Jaapiales ord. nov.: Early diverging clades of Agaricomycetidae were dominated by corticoid forms. Mycologia. 2010, 102: 865-880. 10.3852/09-288.
Murray MG, Thompson WF: Rapid isolation of high molecular weight plant DNA. Nucleic Acids Research. 1980, 8: 4321-4325. 10.1093/nar/8.19.4321.
Gardes M, Bruns TD: ITS primers with enhanced specificity for basidiomycetes-application to the identification of mycorrhizae and rusts. Molecular Ecology. 1993, 2: 113-118. 10.1111/j.1365-294X.1993.tb00005.x.
Vilgalys R, Hester M: Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species. Journal of Bacteriology. 1990, 172: 4238-4246.
Hibbett DS: Phylogenetic evidence for horizontal transmission of group I introns in the nuclear ribosomal DNA of mushroom-forming fungi. Molecular Biology and Evolution. 1996, 13: 903-917.
White TJ, Bruns T, Taylor JW: Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR Protocols: a guide to methods and applications. Edited by: Innis MA, Gelfand DH, Sninsky JJ, White TJ. 1990, San Diego, California: Academic Press, 315-321.
Matheny P, Wang Z, Binder M, Curtis J, Lim Y, Nilsson R, Hughes K, Hofstetter V, Ammirati J, Schoch C, Langer E, Langer G, McLaughlin DJ, Wilson AW, Frøslev T, Ge Z-W, Kerrigan RW, Slot JC, Yang Z-L, Baroni TJ, Fischer M, Hosaka K, Matsuura K, Seidl T, Vauras J, Hibbett DS: Contributions of rpb2 and tef1 to the phylogeny of mushrooms and allies (Basidiomycota, Fungi). Molecular Phylogenetics and Evolution. 2007, 43: 430-451. 10.1016/j.ympev.2006.08.024.
Kauserud H, Schumacher T: Outcrossing or inbreeding: DNA markers provide evidence for type of reproductive mode in Phellinus nigrolimitatus (Basidiomycota). Mycological Research. 2001, 105: 676-683. 10.1017/S0953756201004191.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic acids symposium series. 1999, 41: 95-98.
Goloboff PA, Farris JS, Nixon KC: TNT, a free program for phylogenetic analysis. Cladistics-the International Journal of the Willi Hennig Society. 2008, 24: 774-786.
Farris JS: The future of phylogeny reconstruction. Zoologica Scripta. 1997, 26: 303-311. 10.1111/j.1463-6409.1997.tb00420.x.
Nylander JAA: MrModeltest 2.3. 2004, Program distributed by the author. Evolutionary Biology Centre, Uppsala University
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.
Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Drummond A, Ho S, Phillips M, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biology. 2006, 4: e88-10.1371/journal.pbio.0040088.
Ryberg M, Matheny PB: Dealing with incomplete taxon sampling and diversification of a large clade of mushroom-forming fungi. Evolution. 2011
Smith SY, Currah RS, Stockey RA: Cretaceous and Eocene poroid hymenophores from Vancouver Island, British Columbia. Mycologia. 2004, 96: 180-186. 10.2307/3762001.
Hibbett DS, Grimaldi D, Donoghue MJ: Cretaceous mushrooms in amber. Nature. 1995, 377-487.
Hibbett DS, Grimaldi D, Donoghue MJ: Fossil mushrooms from Miocene and Creaceous amber and the evolution of homobasidiomycetes. Amercian Journal of Botany. 1997, 84: 981-991. 10.2307/2446289.
LePage BA: The evolution, biogeography and palaeoecology of the Pinaceae based on fossil and extant representatives. ISHS Acta Horticulturae. 2003, 615: 29-52.
Ronquist F: Dispersal-vicariance analysis: A new approach to the quantification of historical biogeography. Systematic Biology. 1997, 46: 195-203. 10.1093/sysbio/46.1.195.
Yu Y, Harris AJ, He X: RASP (Reconstruct Ancestal State in Phylogenies).
Ree RH, Moore BR, Webb CO, Donoghue MJ: A likelihood framework for inferring the evolution of geographic range on phylogenetic trees. Evolution. 2005, 59: 2299-2311.
Olsson U, Alström P, Gelang M, Ericson PGP, Sundberg P: Phylogeography of Indonesian and Sino-Himalayan region bush warblers (Cettia, Aves). Molecular Phylogenetics and Evolution. 2006, 41: 556-565. 10.1016/j.ympev.2006.05.009.
Sanmartin I, Van Der Mark P, Ronquist F: Inferring dispersal: a bayesian approach to phyleogeny-based island biogeography, with special reference to the Canary Islands. Journal of Biogeography. 2008, 35: 428-449. 10.1111/j.1365-2699.2008.01885.x.
Donk MA: Notes on Malesian Fungi 1. Bulletin of the Botanic Gardens Buitenzorg. 1948, 3: 473-482.
Won H, Renner SS: Dating dispersal and radiation in the Gymnosperm Gnetum (Gnetales)-clock calibration when outgroup relationships are uncertain. Systematic Biology. 2006, 55: 610-622. 10.1080/10635150600812619.
Crisp MD, West JG, Linder HP: Biogeography of the terrestrial flora. 1999, 1: 2
Upchurch P: Gondwanan break-up: legacies of a lost world?. Trends in Ecology & Evolution. 2008, 23: 229-236. 10.1016/j.tree.2007.11.006.
McLoughlin S: The breakup history of Gondwana and its impact on pre-Cenozoic floristic provincialism. Australian Journal of Botany. 2001, 49: 271-300. 10.1071/BT00023.
Sanmartín I, Ronquist F: Southern Hemisphere biogeography inferred by event-based models: Plant versus animal patterns. Systematic Biology. 2004, 53: 216-243. 10.1080/10635150490423430.
Hugall AF, Lee MSY: The likelihood node densitiy effect and consequences for evolutionary studies of molecular rates. Evolution. 2007, 61: 2293-2307. 10.1111/j.1558-5646.2007.00188.x.
Cook LG, Crisp MD: Not so ancient: the extant crown group of Nothofagus represents a post-Gondwanan radiation. Proceedings of the Royal Society B: Biological Sciences. 2005, 272: 2535-2544. 10.1098/rspb.2005.3219.
Peterson KR, Pfister DH, Bell CD: Cophylogeny and biogeography of the fungal parasite Cyttaria and its host Nothofagus, southern beech. Mycologia. 2010, 102: 1417-1425. 10.3852/10-048.
Lomolino MV, Riddle BR, Brown JH: Biogeography. 2005, Sunderland, Massachusetts: Sinauer Associates, 3
Takamatsu S, Matsuda S, Niinomi S, Havrylenko M: Molecular phylogeny supports a northern hemisphere origin of Golouinomyces (Ascomycota: Erysiphales). Mycological Research. 2006, 110: 1093-1101. 10.1016/j.mycres.2006.07.005.
Sanmartín I, Enghoff H, Ronquist F: Patterns of animal dispersal, vicariance and diversification in the Holartic. Biological Journal of the Linnean Society. 2001, 73: 345-390.
Lydolph MC, Jacobsen J, Arctander P, Gilbert MTP, Gilichinsky DA, Hansen AJ, Willerslev E, Lange L: Beringian paleoecology inferred from permafrost-preserved fungal DNA. Applied and Environmental Microbiology. 2005, 71: 1012-1017. 10.1128/AEM.71.2.1012-1017.2005.
Geml J, Laursen GA, O'Neill K, Nusbaum HC, Taylor DL: Beringian origins and cryptic speciation events in the fly agaric. Molecular Ecology. 2006, 15: 225-239.
Chapela IH, Garbelotto M: Phylogeography and evolution in matsutake and close allies inferred by analyses of ITS sequences and AFLPs. Mycologia. 2004, 96: 730-741. 10.2307/3762107.
Eastwood D, Floudas D, Binder M, Majcherczyk A, Schneider P, Aerts A, Asiegbu FO, Baker SE, Barry K, Bendiksby M, Blumentritt M, Coutinho PM, Cullen D, de Vries RP, Gathman A, Goodell B, Henrissat B, Ihmark K, Kauserud H, Kohler A, LaButti K, Lapidus JL, Lee Y-H, Lindquist E, Lilly W, Lucas S, Morin E, Murat C, Oguiza JA, Park J, et al: The plant cell wall decomposing machinery underlies the functional diversity of forest fungi. Science. 2011
LePage B, Currah R, Stockey R, Rothwell G: Fossil ectomycorrhizae from the Middle Eocene. American Journal of Botany. 1997, 84: 410-412. 10.2307/2446014.
We acknowledge Margit Jarosch and Fungal Herbarium of New Zealand (PDD) for providing study material of Austropaxillus, Andreas Bresinsky for helpful information about Austropaxillus. Teresa Lebel, Marie L. Davey and two referees are acknowledged for helpful comments that improved the manuscript, and the Research Council of Norway and the University of Oslo for financial support. MaB was supported by AFTOL2 (DEB-0732968).
IS and MiB participated in analyzing the data and wrote most of the manuscript. IBE generated data, established alignments, participated in analyzing the data, and wrote parts of the manuscript. MaB generated data, established alignments and commented upon the manuscript. TC participated in analyzing the data and commented upon the manuscript. HK conceived of the study, and participated in its coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Specimens included in this study. A table including origin, Isolate ID and GenBank accession numbers for each specimen in the study. (DOC 236 KB)
Additional file 2: . A 2856 bp allignment of the five nuclear regions 5.8S, nrLSU, nrSSU, rpb2, tef1 for 79 accessions of Serpulaceae taxa including an extensive boletalean outgroup (41 accessions) and two athelialean accessions for rooting the tree. (TXT 223 KB)
Additional file 3: . A 6172 bp allignment of the five nuclear regions 5.8S, nrLSU, nrSSU, rpb2, tef1 for 109 accessions of Boletales and nine additional fungal orders. (TXT 685 KB)
Additional file 4: Allignment of Dataset3. Allignment of the nuclear ribosomal LSU region for 39 accessions of the Serpulaceae and one outgroup. (TXT 39 KB)
Additional file 5: Xml file using an exponential prior distribution for Dataset2. An .xml file created by BEAUti 1.6.1 for analyses of the data using an uncorrelated exponential relaxed molecular clock in BEAST. (XML 717 KB)
Additional file 6: Xml file using a lognormal prior distribution for Dataset2. An .xml file created by BEAUti 1.6.1 for analyses of the data using an uncorrelated lognormal relaxed molecular clock in BEAST. (XML 717 KB)
Additional file 7: Gymnopaxillus. Strict consensus tree from a maximum parsimony analysis of the nuclear ribosomal LSU region of 38 Serpulaceae taxa. Bondarcevomyces taxi was used as outgroup. The 16 most parsimonious trees (MPTs) were 311 steps long and had a rescaled consistency index of 0.90 and homoplasy index of 0.06. Jackknife support values (2000 replicates) are superimposed on the branches, showing that the included Gymnopaxillus spp. (grey shading) are nested within Austropaxillus with high support. (PDF 156 KB)
About this article
Cite this article
Skrede, I., Engh, I.B., Binder, M. et al. Evolutionary history of Serpulaceae (Basidiomycota): molecular phylogeny, historical biogeography and evidence for a single transition of nutritional mode. BMC Evol Biol 11, 230 (2011). https://doi.org/10.1186/1471-2148-11-230
- Maximum Parsimony Analysis
- Nutritional Mode
- High Posterior Density
- Ancestral Area
- Molecular Phylogenetic Approach