BACKGROUND
the globally distributed mushroom-forming genus armillaria contains species that are frequently encountered as tree pathogens in natural forests  <cit>  as well as silvicultural and agronomic systems  <cit> , and are the causal agent of the root disease armillaria root rot  <cit> .
besides being pathogens, armillaria species also play a critical role as decomposers.
extensive infection observed in stumps and roots and long possession of the substrate suggest that armillaria species contribute significantly to decomposition and mineral cycling within many forests  <cit> .
armillaria species form basidiomes—in this case, mushrooms—which serve as reproductive structures that fruit seasonally when conditions are optimal.
each of these basidiomes produces their reproductive propagules, or basidiospores, on an exposed spore-bearing surface, or hymenium.
at maturity, the basidiospores are propelled into the air column via a mechanism of forcible spore discharge known as ballistospory.
basidiospores of other fungi capable of ballistospory have been shown to disperse between continents .
based on molecular clock analyses, the origin of armillaria post-dates the gondwana break-up, leading to the hypothesis that several such long distance basidiospore dispersal events led to the global distribution of armillaria  <cit> .
one change in basidiome form that has occurred repeatedly within mushroom-forming lineages is the shift from producing basidiospores on an exposed hymenium to producing them within an enclosed hymenium.
this loss of an exposed hymenium has resulted in the loss of ballistospory  <cit> .
fungi that undergo this change will hereinafter be referred to as gasteroid fungi as they go through a process called gasteromycetation.
a suite of changes in basidiome morphology, including the development of a fully enclosed basidiospore-bearing mass  that is encased in a specialized covering , typically take place during gasteromycetation.
gasteromycetation is believed to be a unidirectional process, as gasteroid fungi have never been observed to regain the mechanism of ballistospory  <cit> .
in order to persist, gasteroid fungi have evolved a diverse array of dispersal mechanisms that engage exogenous forces.
examples include animal mycophagy for volatile-producing underground truffle-like basidiomycetes  <cit>  and the “bellows” mechanism of puffballs  <cit> .
gasteroid fungi have evolved independently many times from ballistosporic ancestors within the agaricomycetidae .
these ancestors of gasteroid fungi include agaricoid lineages  and boletoid lineages .
to date, gasteroid fungi have only been documented as ectomycorrhizal  or saprotrophic in their ecology .
the first known gasteroid fungus closely related to armillaria was described in  <dig>  <cit> .
the monotypic genus, guyanagaster t. w. henkel, aime & m. e. sm.
and the species g. necrorhizus t. w. henkel, aime & m. e. sm.
are only known from the pakaraima mountains of guyana in the guiana shield region of south america  <cit> .
this fungus produces subhypogeous basidiomes that are often attached to decaying tree roots, has a thick black peridium with pyramidal warts, a reduced stipe, and a tough, gelatinous gleba that changes from white to pink to brick red with maturation  <cit> .
this species is remarkable in the sense that it does not possess the morphological hallmarks of gasteroid fungi adapted for any of the aforementioned dispersal mechanisms .
to our knowledge, until the discovery of guyanagaster, no known gasteroid fungus has been demonstrated to have evolved from a pathogenic lineage.
guyanagaster necrorhizus basidiomes are almost always found fruiting from dead and decaying tree roots, in this case primarily of the forest-dominating trees in the genus dicymbe spruce ex benth.
.
roots to which g. necrorhizus are attached display signs of white rot, indicating a wood decay capability for this fungus.
given the association of g. necrorhizus to dead and decaying roots and its close relationship to armillaria species, it is conjectured that it is also has pathogenic capabilities, but no experimental evidence exists to confirm this.
prior multi-gene phylogenetic analyses by henkel et al.
were equivocal in that some loci suggested g. necrorhizus was sister to armillaria, while other loci showed g. necrorhizus derived from within armillaria  <cit> .
resolving the phylogenetic placement of g. necrorhizus in relation to armillaria will create a more complete picture of armillariod  evolution through time and allow us to understand traits that led to the success of this group.
a resolved phylogeny will also give us a robust framework in which to understand the evolution of guyanagaster.
here we use molecular data from g. necrorhizus and armillaria species from every major infrageneric lineage, including armillaria puiggarii speg.—the only armillaria species to be collected in the same region as g. necrorhizus—to provide a fully resolved phylogenetic hypothesis for armillarioid fungi.
we also erect a new genus, desarmillaria, composed of the exannulate mushroom-forming armillarioid species, to accommodate those species formerly placed in armillaria subgenus desarmillaria.
a second species of guyanagaster, guyanagaster lucianii, is described as new to science.
using these data we produce a time-calibrated phylogeny, which we analyze in combination with ancestral range estimations, in order to examine the morphological and biogeographic evolution of armillarioid fungi.
methods
collecting and morphological analyses
collecting expeditions to the upper potaro river basin  in the west-central pakaraima mountains of guyana were conducted during the rainy seasons of may–july of 2002– <dig>  fungi were collected within a 15-km radius of a previously established base camp  in forests dominated by dicymbe corymbosa spruce ex benth.
and dicymbe altsonii sandw.
<cit> .
in december 2013–january  <dig>  a collecting expedition to the mabura ecological reserve  in mabura hill, guyana—approximately 125 km from the potaro—was conducted.
fungi were collected within a 3-km radius of the field station  in forests with high densities of d. altsonii  <cit> .
new collections of african armillaria species were made in the dja biosphere reserve in the east province of cameroon during the rainy season of august–september  <dig> within a 2-km radius of the dja base camp .
fresh morphological characteristics and substratum relationships were described in the field.
color was described subjectively and coded according to kornerup and wanscher  <cit> , with color plates noted in parentheses.
the descriptions of the collected armillaria specimens were compared to the descriptions of validly published species.
additionally, attempts to obtain specimen cultures were made in the field by placing small pieces of unexposed basidiome tissue on potato dextrose agar   and grown for one month at room temperature.
specimens were then field-dried with silica gel.
micromorphological features on dried specimens were examined in the laboratory using an olympus bh2-rfca compound microscope.
two dried armillaria specimens from africa, th  <dig> and thdja  <dig>  and four dried guyanagaster specimens, mca  <dig>  rak  <dig>  rak  <dig> and rak  <dig>  were rehydrated in 70% ethanol and then sectioned by hand and mounted in water, melzer’s reagent and 1% aqueous congo red.
twenty randomly selected basidiospores were measured per collection under a 100× objective.
length/width q values for basidiospores are reported as qr  and qm .
specimens were deposited in the following herbaria: pul , brg , hsc  and ya .
molecular methods
dna was extracted from basidiome tissue using the wizard® genomic dna purification kit .
specimens housed in the kriebel herbarium at purdue university were used to obtain dna if sequence data for equivalent taxa were not available on the ncbi nucleotide database.
pcr reactions included  <dig>  μl of meangreen 2× taq dna polymerase pcr master mix ,  <dig>  μl of each primer  and approximately 100 ng of dna.
the final pcr reaction volume was 25 μl.
the recommended cycling conditions for each primer pair we used were followed.
to determine the identity of the recently collected armillaria specimens from guyana and africa, as well as to confirm the identity of the recently collected guyanagaster specimens, pcr was performed to acquire sequence data from the internal transcribed spacer  region , using the primer pair its1f/its4b  <cit> .
to analyze the phylogenetic relationship of the armillarioid fungi, pcr was performed on all specimens at the three following loci: nuclear ribosomal large subunit dna , elongation factor 1-α  and rna polymerase ii  genes using the following primer pairs, respectively: lror/lr <dig>  <cit> , 987 f/2218r  <cit>  and brpb2-6 f/brpb2- <dig> r  <cit> .
uncleaned pcr products were sent to beckman coulter, inc.  for sequencing.
sequences were manually edited using sequencher  <dig> .
<dig> .
the its sequences we generated were used as queries in the sequence similarity search tool, ncbi blast, to search the ncbi nucleotide database.
phylogenetic analyses
in order to resolve the phylogeny of guyanagaster and armillaria, we compiled a dataset composed of sequence data from both armillaria and guyanagaster specimens.
two closely related species in the physalacriaceae, oudemansiella mucida and strobilurus esculentus, served as outgroup taxa.
six loci were used , glyceraldehyde-3-phosphate dehydrogenase  and beta-tubulin ) to assess guyanagaster and armillaria systematic relationships through phylogenetic analysis.
not all loci were available for all specimens.
dna sequences that were not generated in this study were obtained from genbank or published whole genome sequences.
collection information for all specimens and genbank accession numbers for the included sequences are compiled in additional file  <dig>
sequences were aligned in mega  <dig>   <cit>  using the muscle algorithm  <cit>  with refinements to the alignment done manually.
the introns of ef1α were alignable across the armillaria/guyanagaster dataset and were therefore included.
phylogenies were reconstructed using maximum likelihood  and bayesian methods.
the gtr + g model of molecular evolution was selected for all data sets as determined by partitionfinder v <dig> .
<dig>  <cit> .
maximum likelihood bootstrap analysis for phylogeny and assessment of the branch support by bootstrap percentages  was performed using raxml v <dig> .
<dig>  <cit> .
one thousand bootstrap replicates were produced.
each locus was analyzed separately as well as all together in a supermatrix data set using ml.
bayesian analyses for the reporting of bayesian posterior probability  support for branches was conducted on individual loci as well as the six-gene data set using the program mr. bayes v <dig> .
<dig>  <cit> .
four simultaneous, independent runs each with four markov chain monte carlo  chains, were initiated and run at a temperature of  <dig>  for  <dig> million generations, sampling trees every  <dig> generations until the standard deviation of the split frequencies reached a final stop value of  <dig> .
we discarded the initial 10% of trees as burn in and produced a maximum clade credibility tree from the remaining trees;  <dig>  bpp represents a well-supported lineage.
divergence time estimation
divergence ages within the armillarioid clade were estimated using the fossil calibration approach described and implemented by .
molecular clock analysis was performed using beast v <dig> .
<dig>  <cit>  with xml files containing beast commands and priors assembled in beauti v <dig> .
<dig>  these files included the following analytical settings: gtr + g model of evolution, uncorrelated relaxed clock with lognormal rate distribution; tree prior was set to speciation birth-death process, running  <dig> million generations, sampling every 1000th tree.
the analysis was run three times.
the first 10% of trees were removed as burn-in after ensuring a minimum effective sample size  of  <dig> was reached for all of the parameters.
the remaining trees—representing the posterior distribution from all bayesian analyses—were combined into a single file using logcombiner v <dig> .
<dig>  this file was used to produce a summary tree using tree annotator v <dig> .
<dig>  the mean ages of supported nodes, as well as the corresponding 95% highest posterior densities , were examined from beast logfiles using tracer v <dig>   <cit> .
taxa used in this analysis are in bold in additional file  <dig>  it includes  <dig> armillaria species , two guyanagaster specimens , nine physalacriaceae species , while six species from agaricales families closely related to the physalacriaceae fide  <cit>  served as outgroup taxa.
the alignment was analyzed using five partitions: 28s was analyzed in a single partition, while two codon partitions , 3) were utilized for both ef1α and rpb <dig>  the marasmioid fungi  were calibrated based on a 90-ma fossil archaeomarasmius leggetti from mid-cretaceous amber  <cit> , following the parameter settings of  <cit> .
in beauti, this prior was set as a lognormal distribution with a mean of  <dig>  log standard deviation of  <dig> and offset of  <dig>  with mean in real space, truncated on the lower end at  <dig> and at the upper end at  <dig>
ancestral range estimation
to estimate the ancestral range of the armillarioid clade and lineages within, we used dispersal-extinction-cladogenesis  analysis in the package lagrange  <cit>  and implemented in the program rasp  <cit> .
six areas were defined: north america, europe + asia , africa, tropical south america , australasia  and temperate south america .
each taxon in our data set was assigned to areas based on its current known range.
in the scenario we tested, movement between any area at any time was unconstrained.
under this model, all areas are treated as equally probable ancestral ranges.
this model was tested under both two and three area constraints.
the consensus tree produced during the divergence time analysis was used for all lagrange calculations and the root node of the armillarioid clade was calibrated to 51 ma.
RESULTS
sequences, collections and cultures
eight its, eleven 28s, twelve ef1α and five rpb <dig> sequences were generated during this study .
the size of the sequences ranged from 514– <dig>  705– <dig>  796– <dig>  and 553–704 bp, respectively.
after the ends of the individual alignments were trimmed, the size of the aligned datasets were as follows: 28s was 906 bp; ef1α was 906 bp; rpb <dig> was 638 bp; actin- <dig> was 607 bp; gpd was 499 bp; tub was 911 bp.
the six-gene dataset was composed of a total of  <dig> armillaria specimens, representing  <dig> species, as well as four guyanagaster specimens, representing two species.
the number of taxa in each of the single-locus phylogenies is as follows: 28s had  <dig> taxa; ef1α had  <dig> taxa; rpb <dig> had  <dig> taxa; actin- <dig> had <dig> taxa; gpd had  <dig> taxa; tub had  <dig> taxa.
phylogenies for each locus, as well as a combination of loci, are available in additional file  <dig>
two armillaria specimens were collected from guyana—mca  <dig> and th  <dig>  its sequences from both shared 99% identity with armillaria puiggarii .
no significant morphological differences were found when comparing our field descriptions to the species description by  <cit> .
two armillaria specimens were collected from cameroon—th  <dig> and th dja  <dig>  its sequences from both shared 99% identity with an unidentified armillaria species collected from zimbabwe .
additionally, there was variation between the two sequences.
no significant morphological differences were found when comparing th  <dig> to the species description of armillaria camerunensis  <cit> , but there was some morphological variation when comparing th dja  <dig>  for now, we conservatively hypothesize that both of these specimens represent armillaria camerunensis.
morphological descriptions of th  <dig> and th dja  <dig> are available in additional file  <dig>
pure vegetative cultures were obtained from guyanagaster specimens mca  <dig> and rak  <dig>  as well as armillaria puiggarii specimen th  <dig>  with rhizomorphs produced in each .
after one month, differences in growth and branching pattern of the rhizomorphs can be observed between the cultures of the two guyanagaster specimens, but the rhizomorphs of neither were observed to become melanized.
in contrast, the rhizomorphs in the culture of a. puiggarii became melanized almost completely.
strain mca  <dig> is deposited in centraalbureau voor schimmelcultures with the strain accession number cbs  <dig>  while strains rak  <dig> and th  <dig> are available from the authors upon request.fig.
<dig>
a vegetative culture of guyanagaster necrorhizus mca  <dig> showing unmelanized rhizomorphs.
b vegetative culture of guyanagaster lucianii rak  <dig> showing unmelanized rhizomorphs.
c vegetative culture of armillaria puiggarii th  <dig> showing melanized rhizomorphs.
bar = one cm
resolved guyanagaster phylogenetic hypothesis
no single locus was sufficient to obtain a well-supported armillarioid phylogeny , but through the use of multiple loci, we were able to obtain a well-supported armillarioid phylogeny .
the six-gene analyses recovered a strongly supported armillarioid clade  , which is composed of all known armillaria and guyanagaster species.
the armillarioid clade is further composed of three well-supported lineages that were recovered from both analyses :  the annulate lineage , which includes annulate mushroom-forming armillarioid species,  the exannulate lineage , composed of exannulate mushroom-forming armillarioid species formerly placed in armillaria subgen.
desarmillaria and  the gasteroid lineage , composed of guyanagaster species.
the gasteroid lineage consists of two distinct species: g. necrorhizus from the potaro region and an undescribed guyanagaster species from the mabura region of guyana.
these two species are geographically separated by mountains of elevations over 1000 m, as well as kaieteur falls and the essequibo river, all within a distance of 125 km.fig.
<dig> phylogram generated from the analysis of six gene regions  from  <dig> taxa.
guyanagaster is the earliest diverging lineage and is sister to the mushroom-forming armillarioid species.
the exannulate armillarioid species are the next to diverge, and compose desarmillaria.
the annulate armillarioid species form a monophyletic lineage and compose armillaria sensu stricto.
the major lineages within armillaria are indicated in bold text at the corresponding node.
the exannulate armillarioid species are sister to the annulate armillarioid species.
strobilurus esculentus and oudemansiella mucida were selected as outgroup taxa.
black circles represent support of 90%  and  <dig>  bpp  or greater, grey circles represent support of  <dig>  bpp or greater and white circles represent support of 75% or greater.
images:  armillaria sinapina;  armillaria hinnulea;  armillaria puiggarii;  armillaria mellea;  armillaria camerunensis;  desarmillaria tabescens;  desarmillaria ectypa,  guyanagaster necrorhizus.
christian schwarz;  j. j. harrison;  todd f. elliott;  rachel a. koch;  terry w. henkel;  stephen d. russell;  tatyana svetasheva)
the annulate lineage is further composed of four well-supported lineages , with species compositions as follows:  african lineage composed of armillaria fuscipes and a. camerunensis;  melleioid lineage composed of a. mellea;  northern hemisphere lineage composed of a. altimontana, a. borealis, a. calvescens, a. cepistipes, a. gallica, a. gemina, a. sinapina and a. solidipes; and  australasian + temperate south american lineage composed of a. affinis, a. fumosa, a. hinnulea, a. limonea, a. luteobubalina, a. novae-zelandiae, a. pallidula and a. puiggarii.
the most glaring conflict between single gene and multi-gene analyses is the placement of the african lineage: in gpd, actin- <dig> and ef1α this lineage is nested within or sister to the australasian + temperate south american lineage; in tub and rpb <dig>  it is sister to a. mellea; in 28s it is sister to guyanagaster.
in the multi-gene analyses, it is consistently recovered as diverging earlier than the a. mellea lineage.
in all of the single gene analyses except for 28s, a close relationship between the exannulate and the gasteroid lineages was recovered.
age and ancestral range of armillaria and guyanagaster
table  <dig> has a complete list of the divergence dates for the well-supported nodes recovered in the phylogenetic analysis.
figure  <dig> is the time-calibrated phylogeny for the group.
the time to most recent common ancestor  of the armillarioid clade estimated in the beast analysis was 51 ma .
the tmrca of the exannulate lineage was estimated at 41 ma .
the tmrca of the annulate armillarioid lineage was estimated at 33 ma .
the date of the divergence between the two guyanagaster species was estimated at 8 ma .table  <dig> ages for nodes, confidence intervals and posterior probability on the phylogenetic tree presented in fig.
3
guyanagaster
desarmillaria + armillaria
desarmillaria
armillaria
a. borealis + a. solidipes + a. gemina
a. altimontana + a. calvescens + a. gallica
a. cepistipes + a. sinapina
a. novae-zelandiae + a. affinis + a. puiggarii
a. hinnulea + a. pallidula + a. fumosa + a. limonea + a. luteobubalina
fig.
<dig> time-calibrated phylogeny generated from bayesian analysis of three gene regions  from  <dig> armillaria species, two desarmillaria species and two guyanagaster species.
geographic origin of each specimen is indicated by the box to the left of the name: red = neotropics, blue = australasia, brown = temperate south america, lime green = eurasia, purple = north america, and turquoise = africa.
boxes at ancestral nodes correspond to the most probable ancestral range at that node as presented in additional file  <dig>  geologic epochs are noted above the time scale.
numbers at nodes correspond to lineages in table  <dig> and additional file  <dig>  dark grey bars correspond to the 95% hpd and correspond to the values in table  <dig>  the map in the lower left hand corner represents the proposed dispersal pattern
we estimated the ancestral ranges for  <dig> well-supported clades.
the most probable area was the same under range constraints of ≤  <dig> and ≤  <dig> for all nodes.
the most probable ancestral range for the most recent common ancestor  of the armillarioid clade is eurasia.
eurasia is also the most probable ancestral range for all of the annulate mushroom-forming lineages containing northern hemisphere taxa.
migration to both africa and australasia from eurasia is the most probable scenario as estimated from this analysis.
additional file  <dig> provides the probabilities for each area being the ancestral range under both area constraints; numbers in bold represent the most probable region at that node.
discussion
our phylogenetic analyses, coupled with morphological data, support a three-genus composition of the armillarioid clade: armillaria contains the annulate mushroom-forming species, desarmillaria contains the exannulate mushroom-forming species, and guyanagaster contains the known gasteroid species.
our estimates of the tmrca of the armillarioid clade suggest it arose 51 ma, which is in accord with a previous estimate  <cit> , even though a different set of taxa, loci and calibration methods were used.
during this time, earth was experiencing one of the warmest intervals of the past  <dig> million years known as the paleocene-eocene thermal maximum   <cit> .
this warm and humid climate led to the widespread expansion of area covered in tropical rainforests.
the ancestral range estimation suggests the most probable area for the origin of the mrca is the eurasian subcontinent , much of which contained tropical rainforests at that time  <cit> .
biogeography and evolution of guyanagaster and desarmillaria
our phylogenetic reconstruction suggests that guyanagaster species evolved within the earliest diverging lineage in the armillarioid clade.
the process of gasteromycetation is believed to be unidirectional; gasteroid fungi evolve from a ballistosporic ancestor, and there is no evidence that this transformation has ever been reversed  <cit> .
with this in mind, the ancestor to the armillarioid lineage was likely a gilled mushroom.
according to our estimations, between  <dig> and 8 ma the mushroom progenitor to guyanagaster underwent the gasteromycetation process, resulting in the evolution of guyanagaster.
at one point in time, this ancient lineage was composed of mushroom-forming species—all of which likely went extinct—as the only known extant members of this lineage are guyanagaster species.
the next lineage to diverge is composed of the two described desarmillaria species: d. ectypa and d. tabescens.
neither of the species in this genus have an annulus at maturity, in contrast to all other known mushroom-forming armillarioid species.
desarmillaria ectypa is different from all other armillarioid species by virtue of its ecology: whereas armillaria species are pathogens and wood decomposers, d. ectypa is restricted to peat bogs and associated with sphagnum  <cit> .
the ecology of d. tabescens is much more typical of armillaria species.
it is a primary pathogen towards introduced eucalyptus species in france, as well as a secondary pathogen to quercus species, but in many cases it is a saprotroph, colonizing stumps of quercus species  <cit> .
both guyanagaster and desarmillaria evolved within the oldest lineages in the armillarioid clade and both genera contain two known species.
armillaria represents the most recently diverged genus-level lineage within the armillarioid clade, yet has over  <dig> described species .
the older ages of the lineages from which guyanagaster and desarmillaria evolved, coupled with the fact that they contain far fewer extant species, suggest these two genera are on a different evolutionary trajectory compared to armillaria.
as stated above, our ancestral range estimation suggests the most probable ancestral range for the mrca to the armillarioid clade is eurasia, whereas species of guyanagaster have only so far been found in guyana.
one question is how the mushroom progenitor to guyanagaster dispersed from eurasia to guyana.
one possibility is that the mushroom progenitor to guyanagaster was widely dispersed within the humid, tropical rainforest habitat that predominated at the time it arose.
it is possible that guyanagaster evolved within the guiana shield and represents a highly specialized endemic to this region.
the two species of desarmillaria are thought to be more thermophilic than their sympatric armillaria species  <cit> .
this temperature adaptation could be a residual characteristic of their origin during the petm and could also explain the paucity of species in this genus.
as the climate changed, species adapted for a cool and arid climate could have competitively displaced the species adapted for a warm and humid climate, in this case, species of desarmillaria and the mushroom progenitor to guyanagaster.
the novel morphology of guyanagaster species suggests that they are highly specialized for their current habitat, and as stated above, the restriction of d. ectypa to peat bogs is a specialized habitat for armillarioid fungi.
habitat specialization has been shown to act as buffer from extinction  <cit> , and could be what has allowed species in these early diverging lineages to avoid competitive displacement by better-adapted armillarioid fungi.
whether d. tabescens occupies a specialized niche that has aided in its persistence to the present remains unknown.
biogeography and spread of armillaria
according to our analyses, armillaria sensu stricto started diversifying approximately 33 ma in eurasia, at which time the earth’s climate shifted from the relatively ice-free world to one with glacial conditions; the tropical climate during which the mrca of the armillarioid clade arose was replaced by a period of severe cooling and drying as well as an increase in seasonality  <cit> .
this led to a contraction in the area covered by tropical rainforests  <cit>  and the geographic expansion of species-poor deciduous vegetation  <cit> .
therefore, we hypothesize that armillaria species represent annulate mushrooms that survived, diversified and radiated with their hosts as they adapted to the drier, cooler and seasonal climate  of the late eocene and beyond.
armillaria species in africa form a monophyletic lineage .
single gene analyses and multi-gene analyses were ambiguous as to where this lineage diverged within the armillarioid clade .
in all multi-gene combinations, the african armillaria lineage is resolved as the earliest diverging within armillaria sensu stricto, which is similar to previous work that has also suggested that this lineage is divergent compared to the rest of the armillarioid lineage  <cit> .
our ancestral area estimation suggests that the ancestral region of the african lineage is eurasia.
during the eocene and oligocene, a shallow seaway separated africa from eurasia, acting as a barrier to dispersal.
in the early miocene , the seaway drained, effectively linking the biotas in these two regions  <cit> .
our molecular clock analysis suggests that extant species in africa started diversifying approximately 15 ma, which fits a migration from eurasia during this time.
many northern hemisphere armillaria species are found in both north america and eurasia .
in a study of a. mellea, which is a species distributed throughout the northern hemisphere, four distinct lineages were identified, corresponding to their locality: 1) western north america, 2) eastern north america, 3) western eurasia and 4) eastern eurasia  <cit> .
they suggest that this species was once widespread and is now in the process of speciation  <cit> .
within this backdrop, it is possible that armillaria species facilitated their dispersal as pathogens on the diverse flora that existed during this time.
as available pathways disappeared, ultimately ceasing migration, the species started to diverge.
from our ancestral range estimation, we can see that there were multiple armillarioid introductions into north america , suggesting a dynamic process of dispersal across the northern hemisphere.
a more comprehensive study on northern hemisphere armillaria species is necessary to determine when and how they migrated between the major landmasses.
armillaria species started diversifying in the temperate southern hemisphere  about 19 ma .
our results suggest that the most probable dispersal pathway was from eurasia to australasia.
in a plant biogeographic meta-analysis of the australasian region, it was found that after 25 ma, there was an increase of flora inputs to australia from the indomalayan region, representing the time when the continental shelf containing southeast asia collided with that containing new guinea and australia  <cit> .
the authors of this study hypothesized that habitat instability in the region at that time generated vacant niches and increased the probability of successful establishment of dispersed lineages compared with the relatively stable, saturated states existing prior to approximately 25 ma  <cit> .
in this scenario both the greater dispersal distance and relative difficulty in establishing viable populations together contributed to the infrequency of successful exchanges prior to 25 ma.
this too follows our hypothesis that armillaria species dispersed as pathogens on flora that existed during this time.
many of the same armillaria species are found in both australasia and temperate south america  <cit> .
this pattern fits with an ancient gondwanan distribution, but the lineage is far too young for this to be the case.
alternatively, the discovery of fossilized nothofagus wood in antarctic deposits from the pliocene  <cit>  suggests this could have been an overland dispersal route until 2 ma.
many of the fungi in this australasian + temperate south american lineage are thought to be close associates with nothofagus  <cit> , suggesting again that they could have dispersed with their plant associates through antarctica.
alternatively, long distance basidiospore dispersal events could account for the geographic disjunction in closely related armillaria species.
because appropriate hosts—species of nothofagus—occur in both australasia and temperate south america  <cit> , establishment after a long-distance basidiospore dispersal event could be more probable.
the rise of rhizomorphs: morphological developments in the armillarioid clade
our phylogenetic hypothesis sheds light on the evolution of rhizomorphs in the armillarioid clade.
rhizomorphs are discrete, filamentous aggregations that extend from a resource base into substrates that may not support their growth, foraging for new resource bases  <cit> .
individuals that produce rhizomorphs can occupy huge swaths of habitat and prevent other organisms from establishing, as is the case with the “humungous fungus”  <cit> .
all armillarioid species have the capacity to form rhizomorphs in culture, but only the later diverging lineages have been observed to form them in nature .table  <dig> trophic strategy, rhizomorph production in nature and known geographic range of species used in this study
armillaria affinis
armillaria altimontana
armillaria borealis
armillaria calvescens
armillaria camerunensis
armillaria cepistipes
armillaria fumosa
armillaria fuscipes
armillaria gallica
armillaria gemina
armillaria hinnulea
armillaria limonea
armillaria luteobubalina
armillaria mellea
armillaria novae-zelandiae
armillaria pallidula
armillaria puiggarii
armillaria sinapina
armillaria solidipes
desarmillaria ectypa
desarmillaria tabescens
eucalyptus, quercus
guyanagaster lucianii
guyanagaster necrorhizus
species in the earliest diverging genera in the armillarioid clade, desarmillaria and guyanagaster, have not been observed to form rhizomorphs in nature, but they do produce unmelanized rhizomorphs in culture  <cit>  .
rhizomorph production of armillaria species from africa, the next lineage to diverge, has been scantily observed in nature  <cit> , but the production of melanized rhizomorphs does occur in culture  <cit> .
the next lineage to diverge is that composed of a. mellea, which produces short-lived melanized rhizomorphs with limited growth in the soil  <cit> .
the most species rich armillarioid lineage, the northern hemisphere and australasian/temperate south america lineage, produce melanized rhizomorphs in nature .
the capacity of armillaria species to produce melanized rhizomorphs in nature is thought to be an adaptation to harsh environments  <cit> .
melanized rhizomorphs can confer advantages like protection against other microbial competitors, translocation of resources, growth from a suitable resource base into an environment that may not support growth, as well as enhancement of inoculum potential  <cit> .
additionally, melanin is thought to promote longevity and survival of rhizomorphs within the soil  <cit> , and has also been found to help other fungi survive in extreme environments  <cit> .
although all species in the armillarioid clade appear to retain the trait of rhizomorph production, only the most recently diverged appear to be adapted for melanized rhizomorph production in the current environment.
in a controlled environment, it was found that only under high oxygen availability and near saturated moisture would d. tabescens produce melanized rhizomorphs  <cit> .
it was concluded that these environmental conditions are sufficiently stringent in the climate of today, and could explain why d. tabescens has not been observed to form rhizomorphs in nature.
it is possible that the production of melanized rhizomorphs could be a driver of diversification in armillaria by conferring a competitive advantage to the species that produce them.
desarmillaria species lack an annulus at maturity, whereas armillaria species have a robust annulus at maturity.
an annulus is the remnant of a partial veil that remains after the pileus expands.
partial veils protect the immature hymenium  <cit> .
the development of a more protected hymenium could be an adaptation of armillaria to drier, more unpredictable habitats that were common when this lineage began to diversify.
CONCLUSIONS
our analyses suggest that the armillarioid clade arose in eurasia during the petm, a time marked by a warm, tropical climate.
two lineages arose during this time: the earliest diverging lineage—which eventually led to the gasteroid genus guyanagaster—and the exannulate mushroom-forming genus desarmillaria.
besides being the oldest lineages in the armillarioid clade, they are also depauperate compared to armillaria.
armillaria diverged after the shift to a cooler and more arid climate at the eocene-oligocene boundary.
the production of melanized rhizomorphs in nature and the development of a protective partial veil could be adaptations that led to the subsequent dispersal and diversification of armillaria in the much harsher, temperate climate.
the success of armillaria could have displaced now-extinct species in the exannulate and gasteroid lineages.
guyanagaster and desarmillaria species have likely persisted to the present because they are highly specialized for their habitat.
formal taxonomic descriptions
desarmillaria
r. a. koch & aime gen. et stat.
nov.
basionym—armillaria ss.
fries subgenus desarmillaria herink, sympozium o václavce obecné .
<dig> september.
lesnicka fakulta vsz brno:  <dig>   <dig>
type—armillaria socialis  herink, sympozium o václavce obecné .
<dig> september.
lesnicka fakulta vsz brno:  <dig>   <dig>
mycobank number—mb  <dig> .
description—basidiomata stipitate, usually in caespitose clusters on wood.
pileus convex to applanate, usually some shade of brown.
lamellae pallid, adnexed to adnate to subdecurrent with lamellulae.
stipe central, various colorations but often concolorous with cap, longitudinally striate, usually with fibrils on upper third.
annulus absent.
basidiospores ellipsoid to spherical, smooth, hyaline, inamyloid.
basidia clavate, four-sterigmate, hyaline.
cheilocystidia clavate, usually resembling basidioles.
pileipellis suprapellis composed of round, ellipsoid, cylindric, utriform, brown to hyaline, verrucose cells; subpellis composed of a layer of compact, shortened hyphae.
clamp connections absent.
rhizomorph production not observed in nature, while unmelanized production in culture has been observed.
saprotrophic to parasitic.
known only from the northern hemisphere.
commentary—desarmillaria includes mushroom-forming armillarioid species that lack an annulus.
this difference in morphology led singer  <cit>  to divide armillaria  into two sections based on the presence or absence of an annulus at maturity.
herink  <cit> , then, recognized armillaria as an annulate subgenus and desarmillaria as an exannulate subgenus.
additionally, members of this genus have not been observed producing rhizomorphs in the field  <cit> , which is in contrast to most armillaria species, most of which do produce rhizomorphs in the field .
rhizomorph production in nature appears to be a lost trait in this genus, as species are still able to form them in culture.
we have opted to elevate desarmillaria to genus-level for two reasons.
first, the absence of an annulus in desarmillaria species and the presence of one in armillaria species is a reliable characteristic to differentiate the two genera.
second, our phylogenetic and molecular clock analyses show desarmillaria is on a separate evolutionary trajectory compared to armillaria, meriting a separate genus.
desarmillaria tabescens
r. a. koch & aime comb.
nov.
basionym—agaricus tabescens scop., scopoli  <dig>  flora carniolica plantas corniolae indigenas no.
1537:  <dig>
mycobank number—mb  <dig> .
desarmillaria ectypa
r. a. koch & aime comb.
nov.
basionym—agaricus ectypus fr., fries  <dig>  syst.
mycol.
1:  <dig>
mycobank number—mb  <dig> .
guyanagaster lucianii
r. a. koch & aime sp.
nov
.
—holotype—brg  <dig>  isotype—pul f <dig> .
mabura ecological reserve, ,  <dig> december  <dig>  r. a. koch  <dig>
mycobank number—mb  <dig> .
representative dna barcode—rak  <dig> , its—genbank ku <dig>  28s—genbank ku <dig>  ef1α—genbank ku <dig>
etymology—lucianii = in honor of “lucian” edmund, a patamona fungal parataxonomist, who discovered the type locality of g. lucianii and has since been invaluable in collecting guyanagaster specimens in the field.
description—figures  <dig> and 4a–g.
basidiomata gasteroid, subhypogeous to hypogeous, scattered, or in linear troops, attached directly to woody roots of eperua falcata and dicymbe altsonii trees; 12–34 mm broad, 8–25 mm tall, globose to subglobose to ovoid and irregularly broadly lobate, dense, base with smooth, sterile concolorous stipe, 5–44 × 2–5 mm, attached directly to the substratum.
peridium dark brown  to black  during all stages of development, moist, tough, covered in 4– <dig> sided polygonal pyramidal warts that come to a shallow point, these  <dig> – <dig>  mm broad,  <dig>  mm tall, sharply differentiated from endoperidium; endoperidium white  to light pink , tough,  <dig> –1 mm thick, composed of matted hyphae transitioning evenly to sterile, white hyphal veins between the glebal locules; gleba composed of well-defined locules and intervening veins of hyphae; locules globose to ovate to subangular,  <dig> –2 mm broad, initially white , maturing in stages to gold , light orange-pink , to deep orange , often evenly, but rarely the locules that border the peridium and columella take on the darker shade first; hyphae separating locules initially white , hyaline and waterlogged in mature specimens; columella well-defined to not, base continuous with stipe, 6–13 × 3–6 mm, initially white , turning to brown , taking on the color of the gleba with age, cottony to gelatinous.
exoperidium individual hyphae of outer layer dark brown, individual hyphae curved to tortuously curved, thick-walled, non-gelatinous, 18–99 × 4–11 μm.
endoperidium 700–800 μm thick, well-differentiated from exoperidium and glebal locules, hyphae hyaline, 3–11 μm.
sterile hyphae separating locules similar to endoperidial hyphae, hyaline, thin-walled, branching, 3– <dig> um wide, organized in a parallel manner.
basidia not observed.
basidiospores globose, dextrinoid;  <dig> – <dig>  ×  <dig> – <dig>  μm ; immature basidiospores hyaline when immature, then darkening to light pink and then rusty brown with maturity.
columella hyphae hyaline, thin-walled, not as tightly packed, some randomly inflated at septa, 4–13 μm wide, random apical cells inflated, 13–26 × 37–57 μm.
clamp connections absent in all tissues.
taste not tested.
odor not tested.fig.
<dig>
a–g
guyanagaster lucianii .
a exterior showing dark brown peridium and irregularly lobed shape.
b subhypogeous habit.
c longitudinal section of basidioma showing immature gleba, developed columella and elongated stipe.
d maturing gleba.
e basidiospores of immature basidioma.
f basidiospores of mature basidioma.
g mature gleba.
bar = one cm in a, b, c, d and g. bar = 20 μm in e and f
habit, habitat and distribution—attached to decaying woody roots of dicymbe altsonii and eperua falcata in tropical rainforests.
known only from the type locality in the mabura ecological reserve of guyana.
specimens examined—guyana.
region  <dig> upper demerara-berbice—mabura hill, elevation 161 m; vicinity of basecamp, eduba stand  <dig>  four basidiomata,  <dig> june  <dig>  mca  <dig> ; vicinity of basecamp, patch  <dig>  one basidioma found at the base of a dead eperua falcata,  <dig> december  <dig>  rak  <dig> ; vicinity of basecamp, patch  <dig>  four basidiomata found at the base of a dead eperua falcata,  <dig> december  <dig>  rak  <dig> ;  <dig>  km north of basecamp, patch  <dig>  nine basidiomata found along the decaying roots of eperua falcata,  <dig> december  <dig>  rak  <dig> .
commentary— phylogenetic evidence suggests that the two lineages of guyanagaster represent distinct species that diverged approximately 8 ma.
the sister species share 87%  nucleotide identity at the its region and 96%  nucleotide identity at lsu.
the drastic morphological changes that occur during the maturation of guyanagaster basidiomata  as well as the morphological plasticity between individual basidiomata  make delimiting species using macromorphological characters unreliable.
however, g. lucianii can be distinguished from g. necrorhizus by its consistently larger basidiospores.
during all stages of maturation, the basidiospores of g. lucianii measure between 21–30 × 21–29 μm, compared to those of g. necrorhizus, which measure 15–19 × 15–18 μm.
additionally, the basidiospores of g. lucianii lack the pedicel that is readily apparent on the basidiospores of g. necrorhizus .
ecologically, g. lucianii appears to occupy the same niche as g. necrorhizus.
basidiomata are hypogeous to subhypogeous and growing from the roots of dead and decaying trees.
both guyanagaster species are white rotters as the decaying roots to which they are attached are all light-colored, spongy and sometimes gelatinous.
the known habitat for g. necrorhizus and g. lucianii are separated by only 125 km, so elucidation of their dispersal strategy is needed to understand the evolutionary forces that led to their speciation.
additional files
additional file 1: collection data for the specimens used ﻿in the phylogenetic analyses and genbank accession numbers for the corresponding sequence data, both generated for this study and obtained from previous studies ﻿.
additional file 2: single-locus and multi-gene phylogenies for the six loci used.
additional file 3: morphological description of armillaria camerunensis th dja  <dig> and th  <dig>
additional file 4: results from the ancestral range estimation analysis.

 the globally distributed mushroom-forming genus armillaria contains species that are frequently encountered as tree pathogens in natural forests as well as silvicultural and agronomic systems, and are the causal agent of the root disease armillaria root rot. besides being pathogens, armillaria species also play a critical role as decomposers. extensive infection observed in stumps and roots and long possession of the substrate suggest that armillaria species contribute significantly to decomposition and mineral cycling within many forests. armillaria species form basidiomes—in this case, mushrooms—which serve as reproductive structures that fruit seasonally when conditions are optimal. each of these basidiomes produces their reproductive propagules, or basidiospores, on an exposed spore-bearing surface, or hymenium. at maturity, the basidiospores are propelled into the air column via a mechanism of forcible spore discharge known as ballistospory. basidiospores of other fungi capable of ballistospory have been shown to disperse between continents. based on molecular clock analyses, the origin of armillaria post-dates the gondwana break-up, leading to the hypothesis that several such long distance basidiospore dispersal events led to the global distribution of armillaria. one change in basidiome form that has occurred repeatedly within mushroom-forming lineages is the shift from producing basidiospores on an exposed hymenium to producing them within an enclosed hymenium. this loss of an exposed hymenium has resulted in the loss of ballistospory. fungi that undergo this change will hereinafter be referred to as gasteroid fungi as they go through a process called gasteromycetation. a suite of changes in basidiome morphology, including the development of a fully enclosed basidiospore-bearing mass that is encased in a specialized covering, typically take place during gasteromycetation. gasteromycetation is believed to be a unidirectional process, as gasteroid fungi have never been observed to regain the mechanism of ballistospory. in order to persist, gasteroid fungi have evolved a diverse array of dispersal mechanisms that engage exogenous forces. examples include animal mycophagy for volatile-producing underground truffle-like basidiomycetes and the “bellows” mechanism of puffballs. gasteroid fungi have evolved independently many times from ballistosporic ancestors within the agaricomycetidae. these ancestors of gasteroid fungi include agaricoid lineages and boletoid lineages. to date, gasteroid fungi have only been documented as ectomycorrhizal or saprotrophic in their ecology. the first known gasteroid fungus closely related to armillaria was described in. the monotypic genus, guyanagaster t. w. henkel, aime & m. e. sm. and the species g. necrorhizus t. w. henkel, aime & m. e. sm. are only known from the pakaraima mountains of guyana in the guiana shield region of south america. this fungus produces subhypogeous basidiomes that are often attached to decaying tree roots, has a thick black peridium with pyramidal warts, a reduced stipe, and a tough, gelatinous gleba that changes from white to pink to brick red with maturation. this species is remarkable in the sense that it does not possess the morphological hallmarks of gasteroid fungi adapted for any of the aforementioned dispersal mechanisms. to our knowledge, until the discovery of guyanagaster, no known gasteroid fungus has been demonstrated to have evolved from a pathogenic lineage. guyanagaster necrorhizus basidiomes are almost always found fruiting from dead and decaying tree roots, in this case primarily of the forest-dominating trees in the genus dicymbe spruce ex benth.. roots to which g. necrorhizus are attached display signs of white rot, indicating a wood decay capability for this fungus. given the association of g. necrorhizus to dead and decaying roots and its close relationship to armillaria species, it is conjectured that it is also has pathogenic capabilities, but no experimental evidence exists to confirm this. prior multi-gene phylogenetic analyses by henkel et al. were equivocal in that some loci suggested g. necrorhizus was sister to armillaria, while other loci showed g. necrorhizus derived from within armillaria. resolving the phylogenetic placement of g. necrorhizus in relation to armillaria will create a more complete picture of armillariod evolution through time and allow us to understand traits that led to the success of this group. a resolved phylogeny will also give us a robust framework in which to understand the evolution of guyanagaster. here we use molecular data from g. necrorhizus and armillaria species from every major infrageneric lineage, including armillaria puiggarii speg.—the only armillaria species to be collected in the same region as g. necrorhizus—to provide a fully resolved phylogenetic hypothesis for armillarioid fungi. we also erect a new genus, desarmillaria, composed of the exannulate mushroom-forming armillarioid species, to accommodate those species formerly placed in armillaria subgenus desarmillaria. a second species of guyanagaster, guyanagaster lucianii, is described as new to science. using these data we produce a time-calibrated phylogeny, which we analyze in combination with ancestral range estimations, in order to examine the morphological and biogeographic evolution of armillarioid fungi. methods collecting and morphological analyses collecting expeditions to the upper potaro river basin in the west-central pakaraima mountains of guyana were conducted during the rainy seasons of may–july of 2002– fungi were collected within a 15-km radius of a previously established base camp in forests dominated by dicymbe corymbosa spruce ex benth. and dicymbe altsonii sandw.. in december 2013–january a collecting expedition to the mabura ecological reserve in mabura hill, guyana—approximately 125 km from the potaro—was conducted. fungi were collected within a 3-km radius of the field station in forests with high densities of d. altsonii. new collections of african armillaria species were made in the dja biosphere reserve in the east province of cameroon during the rainy season of august–september within a 2-km radius of the dja base camp. fresh morphological characteristics and substratum relationships were described in the field. color was described subjectively and coded according to kornerup and wanscher, with color plates noted in parentheses. the descriptions of the collected armillaria specimens were compared to the descriptions of validly published species. additionally, attempts to obtain specimen cultures were made in the field by placing small pieces of unexposed basidiome tissue on potato dextrose agar and grown for one month at room temperature. specimens were then field-dried with silica gel. micromorphological features on dried specimens were examined in the laboratory using an olympus bh2-rfca compound microscope. two dried armillaria specimens from africa, th and thdja and four dried guyanagaster specimens, mca rak rak and rak were rehydrated in 70% ethanol and then sectioned by hand and mounted in water, melzer’s reagent and 1% aqueous congo red. twenty randomly selected basidiospores were measured per collection under a 100× objective. length/width q values for basidiospores are reported as qr and qm. specimens were deposited in the following herbaria: pul, brg, hsc and ya. molecular methods dna was extracted from basidiome tissue using the wizard® genomic dna purification kit. specimens housed in the kriebel herbarium at purdue university were used to obtain dna if sequence data for equivalent taxa were not available on the ncbi nucleotide database. pcr reactions included  μl of meangreen 2× taq dna polymerase pcr master mix,  μl of each primer and approximately 100 ng of dna. the final pcr reaction volume was 25 μl. the recommended cycling conditions for each primer pair we used were followed. to determine the identity of the recently collected armillaria specimens from guyana and africa, as well as to confirm the identity of the recently collected guyanagaster specimens, pcr was performed to acquire sequence data from the internal transcribed spacer region, using the primer pair its1f/its4b. to analyze the phylogenetic relationship of the armillarioid fungi, pcr was performed on all specimens at the three following loci: nuclear ribosomal large subunit dna, elongation factor 1-α and rna polymerase ii genes using the following primer pairs, respectively: lror/lr, 987 f/2218r and brpb2-6 f/brpb2- r. uncleaned pcr products were sent to beckman coulter, inc. for sequencing. sequences were manually edited using sequencher.. the its sequences we generated were used as queries in the sequence similarity search tool, ncbi blast, to search the ncbi nucleotide database. phylogenetic analyses in order to resolve the phylogeny of guyanagaster and armillaria, we compiled a dataset composed of sequence data from both armillaria and guyanagaster specimens. two closely related species in the physalacriaceae, oudemansiella mucida and strobilurus esculentus, served as outgroup taxa. six loci were used, glyceraldehyde-3-phosphate dehydrogenase and beta-tubulin ) to assess guyanagaster and armillaria systematic relationships through phylogenetic analysis. not all loci were available for all specimens. dna sequences that were not generated in this study were obtained from genbank or published whole genome sequences. collection information for all specimens and genbank accession numbers for the included sequences are compiled in additional file sequences were aligned in mega using the muscle algorithm with refinements to the alignment done manually. the introns of ef1α were alignable across the armillaria/guyanagaster dataset and were therefore included. phylogenies were reconstructed using maximum likelihood and bayesian methods. the gtr + g model of molecular evolution was selected for all data sets as determined by partitionfinder v.. maximum likelihood bootstrap analysis for phylogeny and assessment of the branch support by bootstrap percentages was performed using raxml v.. one thousand bootstrap replicates were produced. each locus was analyzed separately as well as all together in a supermatrix data set using ml. bayesian analyses for the reporting of bayesian posterior probability support for branches was conducted on individual loci as well as the six-gene data set using the program mr. bayes v.. four simultaneous, independent runs each with four markov chain monte carlo chains, were initiated and run at a temperature of for million generations, sampling trees every generations until the standard deviation of the split frequencies reached a final stop value of. we discarded the initial 10% of trees as burn in and produced a maximum clade credibility tree from the remaining trees; bpp represents a well-supported lineage. divergence time estimation divergence ages within the armillarioid clade were estimated using the fossil calibration approach described and implemented by. molecular clock analysis was performed using beast v. with xml files containing beast commands and priors assembled in beauti v. these files included the following analytical settings: gtr + g model of evolution, uncorrelated relaxed clock with lognormal rate distribution; tree prior was set to speciation birth-death process, running million generations, sampling every 1000th tree. the analysis was run three times. the first 10% of trees were removed as burn-in after ensuring a minimum effective sample size of was reached for all of the parameters. the remaining trees—representing the posterior distribution from all bayesian analyses—were combined into a single file using logcombiner v. this file was used to produce a summary tree using tree annotator v. the mean ages of supported nodes, as well as the corresponding 95% highest posterior densities, were examined from beast logfiles using tracer v. taxa used in this analysis are in bold in additional file it includes armillaria species, two guyanagaster specimens, nine physalacriaceae species, while six species from agaricales families closely related to the physalacriaceae fide served as outgroup taxa. the alignment was analyzed using five partitions: 28s was analyzed in a single partition, while two codon partitions, 3) were utilized for both ef1α and rpb the marasmioid fungi were calibrated based on a 90-ma fossil archaeomarasmius leggetti from mid-cretaceous amber, following the parameter settings of. in beauti, this prior was set as a lognormal distribution with a mean of log standard deviation of and offset of with mean in real space, truncated on the lower end at and at the upper end at ancestral range estimation to estimate the ancestral range of the armillarioid clade and lineages within, we used dispersal-extinction-cladogenesis analysis in the package lagrange and implemented in the program rasp. six areas were defined: north america, europe + asia, africa, tropical south america, australasia and temperate south america. each taxon in our data set was assigned to areas based on its current known range. in the scenario we tested, movement between any area at any time was unconstrained. under this model, all areas are treated as equally probable ancestral ranges. this model was tested under both two and three area constraints. the consensus tree produced during the divergence time analysis was used for all lagrange calculations and the root node of the armillarioid clade was calibrated to 51 ma.  sequences, collections and cultures eight its, eleven 28s, twelve ef1α and five rpb sequences were generated during this study. the size of the sequences ranged from 514– 705– 796– and 553–704 bp, respectively. after the ends of the individual alignments were trimmed, the size of the aligned datasets were as follows: 28s was 906 bp; ef1α was 906 bp; rpb was 638 bp; actin- was 607 bp; gpd was 499 bp; tub was 911 bp. the six-gene dataset was composed of a total of armillaria specimens, representing species, as well as four guyanagaster specimens, representing two species. the number of taxa in each of the single-locus phylogenies is as follows: 28s had taxa; ef1α had taxa; rpb had taxa; actin- had taxa; gpd had taxa; tub had taxa. phylogenies for each locus, as well as a combination of loci, are available in additional file two armillaria specimens were collected from guyana—mca and th its sequences from both shared 99% identity with armillaria puiggarii. no significant morphological differences were found when comparing our field descriptions to the species description by. two armillaria specimens were collected from cameroon—th and th dja its sequences from both shared 99% identity with an unidentified armillaria species collected from zimbabwe. additionally, there was variation between the two sequences. no significant morphological differences were found when comparing th to the species description of armillaria camerunensis, but there was some morphological variation when comparing th dja for now, we conservatively hypothesize that both of these specimens represent armillaria camerunensis. morphological descriptions of th and th dja are available in additional file pure vegetative cultures were obtained from guyanagaster specimens mca and rak as well as armillaria puiggarii specimen th with rhizomorphs produced in each. after one month, differences in growth and branching pattern of the rhizomorphs can be observed between the cultures of the two guyanagaster specimens, but the rhizomorphs of neither were observed to become melanized. in contrast, the rhizomorphs in the culture of a. puiggarii became melanized almost completely. strain mca is deposited in centraalbureau voor schimmelcultures with the strain accession number cbs while strains rak and th are available from the authors upon request.fig. a vegetative culture of guyanagaster necrorhizus mca showing unmelanized rhizomorphs. b vegetative culture of guyanagaster lucianii rak showing unmelanized rhizomorphs. c vegetative culture of armillaria puiggarii th showing melanized rhizomorphs. bar = one cm resolved guyanagaster phylogenetic hypothesis no single locus was sufficient to obtain a well-supported armillarioid phylogeny, but through the use of multiple loci, we were able to obtain a well-supported armillarioid phylogeny. the six-gene analyses recovered a strongly supported armillarioid clade, which is composed of all known armillaria and guyanagaster species. the armillarioid clade is further composed of three well-supported lineages that were recovered from both analyses: the annulate lineage, which includes annulate mushroom-forming armillarioid species, the exannulate lineage, composed of exannulate mushroom-forming armillarioid species formerly placed in armillaria subgen. desarmillaria and the gasteroid lineage, composed of guyanagaster species. the gasteroid lineage consists of two distinct species: g. necrorhizus from the potaro region and an undescribed guyanagaster species from the mabura region of guyana. these two species are geographically separated by mountains of elevations over 1000 m, as well as kaieteur falls and the essequibo river, all within a distance of 125 km.fig. phylogram generated from the analysis of six gene regions from taxa. guyanagaster is the earliest diverging lineage and is sister to the mushroom-forming armillarioid species. the exannulate armillarioid species are the next to diverge, and compose desarmillaria. the annulate armillarioid species form a monophyletic lineage and compose armillaria sensu stricto. the major lineages within armillaria are indicated in bold text at the corresponding node. the exannulate armillarioid species are sister to the annulate armillarioid species. strobilurus esculentus and oudemansiella mucida were selected as outgroup taxa. black circles represent support of 90% and bpp or greater, grey circles represent support of bpp or greater and white circles represent support of 75% or greater. images: armillaria sinapina; armillaria hinnulea; armillaria puiggarii; armillaria mellea; armillaria camerunensis; desarmillaria tabescens; desarmillaria ectypa, guyanagaster necrorhizus. christian schwarz; j. j. harrison; todd f. elliott; rachel a. koch; terry w. henkel; stephen d. russell; tatyana svetasheva) the annulate lineage is further composed of four well-supported lineages, with species compositions as follows: african lineage composed of armillaria fuscipes and a. camerunensis; melleioid lineage composed of a. mellea; northern hemisphere lineage composed of a. altimontana, a. borealis, a. calvescens, a. cepistipes, a. gallica, a. gemina, a. sinapina and a. solidipes; and australasian + temperate south american lineage composed of a. affinis, a. fumosa, a. hinnulea, a. limonea, a. luteobubalina, a. novae-zelandiae, a. pallidula and a. puiggarii. the most glaring conflict between single gene and multi-gene analyses is the placement of the african lineage: in gpd, actin- and ef1α this lineage is nested within or sister to the australasian + temperate south american lineage; in tub and rpb it is sister to a. mellea; in 28s it is sister to guyanagaster. in the multi-gene analyses, it is consistently recovered as diverging earlier than the a. mellea lineage. in all of the single gene analyses except for 28s, a close relationship between the exannulate and the gasteroid lineages was recovered. age and ancestral range of armillaria and guyanagaster table  has a complete list of the divergence dates for the well-supported nodes recovered in the phylogenetic analysis. figure  is the time-calibrated phylogeny for the group. the time to most recent common ancestor of the armillarioid clade estimated in the beast analysis was 51 ma. the tmrca of the exannulate lineage was estimated at 41 ma. the tmrca of the annulate armillarioid lineage was estimated at 33 ma. the date of the divergence between the two guyanagaster species was estimated at 8 ma.table ages for nodes, confidence intervals and posterior probability on the phylogenetic tree presented in fig. 3 guyanagaster desarmillaria + armillaria desarmillaria armillaria a. borealis + a. solidipes + a. gemina a. altimontana + a. calvescens + a. gallica a. cepistipes + a. sinapina a. novae-zelandiae + a. affinis + a. puiggarii a. hinnulea + a. pallidula + a. fumosa + a. limonea + a. luteobubalina fig. time-calibrated phylogeny generated from bayesian analysis of three gene regions from armillaria species, two desarmillaria species and two guyanagaster species. geographic origin of each specimen is indicated by the box to the left of the name: red = neotropics, blue = australasia, brown = temperate south america, lime green = eurasia, purple = north america, and turquoise = africa. boxes at ancestral nodes correspond to the most probable ancestral range at that node as presented in additional file geologic epochs are noted above the time scale. numbers at nodes correspond to lineages in table  and additional file dark grey bars correspond to the 95% hpd and correspond to the values in table  the map in the lower left hand corner represents the proposed dispersal pattern we estimated the ancestral ranges for well-supported clades. the most probable area was the same under range constraints of ≤  and ≤  for all nodes. the most probable ancestral range for the most recent common ancestor of the armillarioid clade is eurasia. eurasia is also the most probable ancestral range for all of the annulate mushroom-forming lineages containing northern hemisphere taxa. migration to both africa and australasia from eurasia is the most probable scenario as estimated from this analysis. additional file provides the probabilities for each area being the ancestral range under both area constraints; numbers in bold represent the most probable region at that node. discussion our phylogenetic analyses, coupled with morphological data, support a three-genus composition of the armillarioid clade: armillaria contains the annulate mushroom-forming species, desarmillaria contains the exannulate mushroom-forming species, and guyanagaster contains the known gasteroid species. our estimates of the tmrca of the armillarioid clade suggest it arose 51 ma, which is in accord with a previous estimate, even though a different set of taxa, loci and calibration methods were used. during this time, earth was experiencing one of the warmest intervals of the past million years known as the paleocene-eocene thermal maximum. this warm and humid climate led to the widespread expansion of area covered in tropical rainforests. the ancestral range estimation suggests the most probable area for the origin of the mrca is the eurasian subcontinent, much of which contained tropical rainforests at that time. biogeography and evolution of guyanagaster and desarmillaria our phylogenetic reconstruction suggests that guyanagaster species evolved within the earliest diverging lineage in the armillarioid clade. the process of gasteromycetation is believed to be unidirectional; gasteroid fungi evolve from a ballistosporic ancestor, and there is no evidence that this transformation has ever been reversed. with this in mind, the ancestor to the armillarioid lineage was likely a gilled mushroom. according to our estimations, between and 8 ma the mushroom progenitor to guyanagaster underwent the gasteromycetation process, resulting in the evolution of guyanagaster. at one point in time, this ancient lineage was composed of mushroom-forming species—all of which likely went extinct—as the only known extant members of this lineage are guyanagaster species. the next lineage to diverge is composed of the two described desarmillaria species: d. ectypa and d. tabescens. neither of the species in this genus have an annulus at maturity, in contrast to all other known mushroom-forming armillarioid species. desarmillaria ectypa is different from all other armillarioid species by virtue of its ecology: whereas armillaria species are pathogens and wood decomposers, d. ectypa is restricted to peat bogs and associated with sphagnum. the ecology of d. tabescens is much more typical of armillaria species. it is a primary pathogen towards introduced eucalyptus species in france, as well as a secondary pathogen to quercus species, but in many cases it is a saprotroph, colonizing stumps of quercus species. both guyanagaster and desarmillaria evolved within the oldest lineages in the armillarioid clade and both genera contain two known species. armillaria represents the most recently diverged genus-level lineage within the armillarioid clade, yet has over described species. the older ages of the lineages from which guyanagaster and desarmillaria evolved, coupled with the fact that they contain far fewer extant species, suggest these two genera are on a different evolutionary trajectory compared to armillaria. as stated above, our ancestral range estimation suggests the most probable ancestral range for the mrca to the armillarioid clade is eurasia, whereas species of guyanagaster have only so far been found in guyana. one question is how the mushroom progenitor to guyanagaster dispersed from eurasia to guyana. one possibility is that the mushroom progenitor to guyanagaster was widely dispersed within the humid, tropical rainforest habitat that predominated at the time it arose. it is possible that guyanagaster evolved within the guiana shield and represents a highly specialized endemic to this region. the two species of desarmillaria are thought to be more thermophilic than their sympatric armillaria species. this temperature adaptation could be a residual characteristic of their origin during the petm and could also explain the paucity of species in this genus. as the climate changed, species adapted for a cool and arid climate could have competitively displaced the species adapted for a warm and humid climate, in this case, species of desarmillaria and the mushroom progenitor to guyanagaster. the novel morphology of guyanagaster species suggests that they are highly specialized for their current habitat, and as stated above, the restriction of d. ectypa to peat bogs is a specialized habitat for armillarioid fungi. habitat specialization has been shown to act as buffer from extinction, and could be what has allowed species in these early diverging lineages to avoid competitive displacement by better-adapted armillarioid fungi. whether d. tabescens occupies a specialized niche that has aided in its persistence to the present remains unknown. biogeography and spread of armillaria according to our analyses, armillaria sensu stricto started diversifying approximately 33 ma in eurasia, at which time the earth’s climate shifted from the relatively ice-free world to one with glacial conditions; the tropical climate during which the mrca of the armillarioid clade arose was replaced by a period of severe cooling and drying as well as an increase in seasonality. this led to a contraction in the area covered by tropical rainforests and the geographic expansion of species-poor deciduous vegetation. therefore, we hypothesize that armillaria species represent annulate mushrooms that survived, diversified and radiated with their hosts as they adapted to the drier, cooler and seasonal climate of the late eocene and beyond. armillaria species in africa form a monophyletic lineage. single gene analyses and multi-gene analyses were ambiguous as to where this lineage diverged within the armillarioid clade. in all multi-gene combinations, the african armillaria lineage is resolved as the earliest diverging within armillaria sensu stricto, which is similar to previous work that has also suggested that this lineage is divergent compared to the rest of the armillarioid lineage. our ancestral area estimation suggests that the ancestral region of the african lineage is eurasia. during the eocene and oligocene, a shallow seaway separated africa from eurasia, acting as a barrier to dispersal. in the early miocene, the seaway drained, effectively linking the biotas in these two regions. our molecular clock analysis suggests that extant species in africa started diversifying approximately 15 ma, which fits a migration from eurasia during this time. many northern hemisphere armillaria species are found in both north america and eurasia. in a study of a. mellea, which is a species distributed throughout the northern hemisphere, four distinct lineages were identified, corresponding to their locality: 1) western north america, 2) eastern north america, 3) western eurasia and 4) eastern eurasia. they suggest that this species was once widespread and is now in the process of speciation. within this backdrop, it is possible that armillaria species facilitated their dispersal as pathogens on the diverse flora that existed during this time. as available pathways disappeared, ultimately ceasing migration, the species started to diverge. from our ancestral range estimation, we can see that there were multiple armillarioid introductions into north america, suggesting a dynamic process of dispersal across the northern hemisphere. a more comprehensive study on northern hemisphere armillaria species is necessary to determine when and how they migrated between the major landmasses. armillaria species started diversifying in the temperate southern hemisphere about 19 ma. our results suggest that the most probable dispersal pathway was from eurasia to australasia. in a plant biogeographic meta-analysis of the australasian region, it was found that after 25 ma, there was an increase of flora inputs to australia from the indomalayan region, representing the time when the continental shelf containing southeast asia collided with that containing new guinea and australia. the authors of this study hypothesized that habitat instability in the region at that time generated vacant niches and increased the probability of successful establishment of dispersed lineages compared with the relatively stable, saturated states existing prior to approximately 25 ma. in this scenario both the greater dispersal distance and relative difficulty in establishing viable populations together contributed to the infrequency of successful exchanges prior to 25 ma. this too follows our hypothesis that armillaria species dispersed as pathogens on flora that existed during this time. many of the same armillaria species are found in both australasia and temperate south america. this pattern fits with an ancient gondwanan distribution, but the lineage is far too young for this to be the case. alternatively, the discovery of fossilized nothofagus wood in antarctic deposits from the pliocene suggests this could have been an overland dispersal route until 2 ma. many of the fungi in this australasian + temperate south american lineage are thought to be close associates with nothofagus, suggesting again that they could have dispersed with their plant associates through antarctica. alternatively, long distance basidiospore dispersal events could account for the geographic disjunction in closely related armillaria species. because appropriate hosts—species of nothofagus—occur in both australasia and temperate south america, establishment after a long-distance basidiospore dispersal event could be more probable. the rise of rhizomorphs: morphological developments in the armillarioid clade our phylogenetic hypothesis sheds light on the evolution of rhizomorphs in the armillarioid clade. rhizomorphs are discrete, filamentous aggregations that extend from a resource base into substrates that may not support their growth, foraging for new resource bases. individuals that produce rhizomorphs can occupy huge swaths of habitat and prevent other organisms from establishing, as is the case with the “humungous fungus”. all armillarioid species have the capacity to form rhizomorphs in culture, but only the later diverging lineages have been observed to form them in nature.table trophic strategy, rhizomorph production in nature and known geographic range of species used in this study armillaria affinis armillaria altimontana armillaria borealis armillaria calvescens armillaria camerunensis armillaria cepistipes armillaria fumosa armillaria fuscipes armillaria gallica armillaria gemina armillaria hinnulea armillaria limonea armillaria luteobubalina armillaria mellea armillaria novae-zelandiae armillaria pallidula armillaria puiggarii armillaria sinapina armillaria solidipes desarmillaria ectypa desarmillaria tabescens eucalyptus, quercus guyanagaster lucianii guyanagaster necrorhizus species in the earliest diverging genera in the armillarioid clade, desarmillaria and guyanagaster, have not been observed to form rhizomorphs in nature, but they do produce unmelanized rhizomorphs in culture. rhizomorph production of armillaria species from africa, the next lineage to diverge, has been scantily observed in nature, but the production of melanized rhizomorphs does occur in culture. the next lineage to diverge is that composed of a. mellea, which produces short-lived melanized rhizomorphs with limited growth in the soil. the most species rich armillarioid lineage, the northern hemisphere and australasian/temperate south america lineage, produce melanized rhizomorphs in nature. the capacity of armillaria species to produce melanized rhizomorphs in nature is thought to be an adaptation to harsh environments. melanized rhizomorphs can confer advantages like protection against other microbial competitors, translocation of resources, growth from a suitable resource base into an environment that may not support growth, as well as enhancement of inoculum potential. additionally, melanin is thought to promote longevity and survival of rhizomorphs within the soil, and has also been found to help other fungi survive in extreme environments. although all species in the armillarioid clade appear to retain the trait of rhizomorph production, only the most recently diverged appear to be adapted for melanized rhizomorph production in the current environment. in a controlled environment, it was found that only under high oxygen availability and near saturated moisture would d. tabescens produce melanized rhizomorphs. it was concluded that these environmental conditions are sufficiently stringent in the climate of today, and could explain why d. tabescens has not been observed to form rhizomorphs in nature. it is possible that the production of melanized rhizomorphs could be a driver of diversification in armillaria by conferring a competitive advantage to the species that produce them. desarmillaria species lack an annulus at maturity, whereas armillaria species have a robust annulus at maturity. an annulus is the remnant of a partial veil that remains after the pileus expands. partial veils protect the immature hymenium. the development of a more protected hymenium could be an adaptation of armillaria to drier, more unpredictable habitats that were common when this lineage began to diversify.  our analyses suggest that the armillarioid clade arose in eurasia during the petm, a time marked by a warm, tropical climate. two lineages arose during this time: the earliest diverging lineage—which eventually led to the gasteroid genus guyanagaster—and the exannulate mushroom-forming genus desarmillaria. besides being the oldest lineages in the armillarioid clade, they are also depauperate compared to armillaria. armillaria diverged after the shift to a cooler and more arid climate at the eocene-oligocene boundary. the production of melanized rhizomorphs in nature and the development of a protective partial veil could be adaptations that led to the subsequent dispersal and diversification of armillaria in the much harsher, temperate climate. the success of armillaria could have displaced now-extinct species in the exannulate and gasteroid lineages. guyanagaster and desarmillaria species have likely persisted to the present because they are highly specialized for their habitat. formal taxonomic descriptions desarmillaria r. a. koch & aime gen. et stat. nov. basionym—armillaria ss. fries subgenus desarmillaria herink, sympozium o václavce obecné. september. lesnicka fakulta vsz brno: type—armillaria socialis herink, sympozium o václavce obecné. september. lesnicka fakulta vsz brno: mycobank number—mb. description—basidiomata stipitate, usually in caespitose clusters on wood. pileus convex to applanate, usually some shade of brown. lamellae pallid, adnexed to adnate to subdecurrent with lamellulae. stipe central, various colorations but often concolorous with cap, longitudinally striate, usually with fibrils on upper third. annulus absent. basidiospores ellipsoid to spherical, smooth, hyaline, inamyloid. basidia clavate, four-sterigmate, hyaline. cheilocystidia clavate, usually resembling basidioles. pileipellis suprapellis composed of round, ellipsoid, cylindric, utriform, brown to hyaline, verrucose cells; subpellis composed of a layer of compact, shortened hyphae. clamp connections absent. rhizomorph production not observed in nature, while unmelanized production in culture has been observed. saprotrophic to parasitic. known only from the northern hemisphere. commentary—desarmillaria includes mushroom-forming armillarioid species that lack an annulus. this difference in morphology led singer to divide armillaria into two sections based on the presence or absence of an annulus at maturity. herink, then, recognized armillaria as an annulate subgenus and desarmillaria as an exannulate subgenus. additionally, members of this genus have not been observed producing rhizomorphs in the field, which is in contrast to most armillaria species, most of which do produce rhizomorphs in the field. rhizomorph production in nature appears to be a lost trait in this genus, as species are still able to form them in culture. we have opted to elevate desarmillaria to genus-level for two reasons. first, the absence of an annulus in desarmillaria species and the presence of one in armillaria species is a reliable characteristic to differentiate the two genera. second, our phylogenetic and molecular clock analyses show desarmillaria is on a separate evolutionary trajectory compared to armillaria, meriting a separate genus. desarmillaria tabescens r. a. koch & aime comb. nov. basionym—agaricus tabescens scop., scopoli flora carniolica plantas corniolae indigenas no. 1537: mycobank number—mb. desarmillaria ectypa r. a. koch & aime comb. nov. basionym—agaricus ectypus fr., fries syst. mycol. 1: mycobank number—mb. guyanagaster lucianii r. a. koch & aime sp. nov. —holotype—brg isotype—pul f. mabura ecological reserve,, december r. a. koch mycobank number—mb. representative dna barcode—rak, its—genbank ku 28s—genbank ku ef1α—genbank ku etymology—lucianii = in honor of “lucian” edmund, a patamona fungal parataxonomist, who discovered the type locality of g. lucianii and has since been invaluable in collecting guyanagaster specimens in the field. description—figures  and 4a–g. basidiomata gasteroid, subhypogeous to hypogeous, scattered, or in linear troops, attached directly to woody roots of eperua falcata and dicymbe altsonii trees; 12–34 mm broad, 8–25 mm tall, globose to subglobose to ovoid and irregularly broadly lobate, dense, base with smooth, sterile concolorous stipe, 5–44 × 2–5 mm, attached directly to the substratum. peridium dark brown to black during all stages of development, moist, tough, covered in 4– sided polygonal pyramidal warts that come to a shallow point, these –  mm broad,  mm tall, sharply differentiated from endoperidium; endoperidium white to light pink, tough, –1 mm thick, composed of matted hyphae transitioning evenly to sterile, white hyphal veins between the glebal locules; gleba composed of well-defined locules and intervening veins of hyphae; locules globose to ovate to subangular, –2 mm broad, initially white, maturing in stages to gold, light orange-pink, to deep orange, often evenly, but rarely the locules that border the peridium and columella take on the darker shade first; hyphae separating locules initially white, hyaline and waterlogged in mature specimens; columella well-defined to not, base continuous with stipe, 6–13 × 3–6 mm, initially white, turning to brown, taking on the color of the gleba with age, cottony to gelatinous. exoperidium individual hyphae of outer layer dark brown, individual hyphae curved to tortuously curved, thick-walled, non-gelatinous, 18–99 × 4–11 μm. endoperidium 700–800 μm thick, well-differentiated from exoperidium and glebal locules, hyphae hyaline, 3–11 μm. sterile hyphae separating locules similar to endoperidial hyphae, hyaline, thin-walled, branching, 3– um wide, organized in a parallel manner. basidia not observed. basidiospores globose, dextrinoid; –  ×  –  μm; immature basidiospores hyaline when immature, then darkening to light pink and then rusty brown with maturity. columella hyphae hyaline, thin-walled, not as tightly packed, some randomly inflated at septa, 4–13 μm wide, random apical cells inflated, 13–26 × 37–57 μm. clamp connections absent in all tissues. taste not tested. odor not tested.fig. a–g guyanagaster lucianii. a exterior showing dark brown peridium and irregularly lobed shape. b subhypogeous habit. c longitudinal section of basidioma showing immature gleba, developed columella and elongated stipe. d maturing gleba. e basidiospores of immature basidioma. f basidiospores of mature basidioma. g mature gleba. bar = one cm in a, b, c, d and g. bar = 20 μm in e and f habit, habitat and distribution—attached to decaying woody roots of dicymbe altsonii and eperua falcata in tropical rainforests. known only from the type locality in the mabura ecological reserve of guyana. specimens examined—guyana. region upper demerara-berbice—mabura hill, elevation 161 m; vicinity of basecamp, eduba stand four basidiomata, june mca; vicinity of basecamp, patch one basidioma found at the base of a dead eperua falcata, december rak; vicinity of basecamp, patch four basidiomata found at the base of a dead eperua falcata, december rak;  km north of basecamp, patch nine basidiomata found along the decaying roots of eperua falcata, december rak. commentary— phylogenetic evidence suggests that the two lineages of guyanagaster represent distinct species that diverged approximately 8 ma. the sister species share 87% nucleotide identity at the its region and 96% nucleotide identity at lsu. the drastic morphological changes that occur during the maturation of guyanagaster basidiomata as well as the morphological plasticity between individual basidiomata make delimiting species using macromorphological characters unreliable. however, g. lucianii can be distinguished from g. necrorhizus by its consistently larger basidiospores. during all stages of maturation, the basidiospores of g. lucianii measure between 21–30 × 21–29 μm, compared to those of g. necrorhizus, which measure 15–19 × 15–18 μm. additionally, the basidiospores of g. lucianii lack the pedicel that is readily apparent on the basidiospores of g. necrorhizus. ecologically, g. lucianii appears to occupy the same niche as g. necrorhizus. basidiomata are hypogeous to subhypogeous and growing from the roots of dead and decaying trees. both guyanagaster species are white rotters as the decaying roots to which they are attached are all light-colored, spongy and sometimes gelatinous. the known habitat for g. necrorhizus and g. lucianii are separated by only 125 km, so elucidation of their dispersal strategy is needed to understand the evolutionary forces that led to their speciation. additional files additional file 1: collection data for the specimens used ﻿in the phylogenetic analyses and genbank accession numbers for the corresponding sequence data, both generated for this study and obtained from previous studies ﻿. additional file 2: single-locus and multi-gene phylogenies for the six loci used. additional file 3: morphological description of armillaria camerunensis th dja and th additional file 4: results from the ancestral range estimation analysis. 
this subgenus is elevated to genus level to accommodate the exannulate mushroom-forming armillarioid species. the final lineage to diverge is that composed of annulate mushroom-forming armillarioid species, in what is now armillaria sensu stricto. the six-locus phylogenetic analysis resolves guyanagaster as the earliest diverging lineage in the armillarioid clade. three loci—28s, ef1α and rpb2—were analyzed in a partitioned nucleotide data set to infer divergence dates and ancestral range estimations for well-supported, monophyletic lineages. in this study, we conducted phylogenetic analyses to fully resolve the relationship of g. necrorhizus with armillaria. the molecular clock analysis and ancestral range estimation suggest the most recent common ancestor to the armillarioid lineage arose million years ago in eurasia. the armillarioid lineage evolved in eurasia during the height of tropical rainforest expansion about million years ago, a time marked by a warm and wet global climate. desarmillaria represents an armillarioid lineage that was likely much more widespread in the past. the production of rhizomorphs is an adaptation to harsh environments, and could be a driver of diversification in armillaria by conferring a competitive advantage to the species that produce them. the first gasteroid fungus with close affinities to armillaria—guyanagaster necrorhizus—was described from the neotropical rainforests of guyana. guyanagaster likely evolved from a gilled mushroom ancestor and could represent a highly specialized endemic in the guiana shield. armillaria root rotcamerooneocenefungal taxonomygasteromycetationguiana shieldmelaninmushroom evolutionphysalacriaceaesystematicshttp://dx.doi.org/ /100000001national science foundationdeb-0918591aime m. catherine issue-copyright-statement© the author 2017 armillaria species represent those that evolved after the shift in climate from warm and tropical to cool and arid during the late eocene. data sets containing guyanagaster from two collecting localities, along with a global sampling of armillaria species—including newly collected specimens from guyana and africa—at six loci were used. no species in either desarmillaria or guyanagaster are known to produce melanized rhizomorphs in nature, whereas almost all armillaria species are known to produce them. the next lineage to diverge is that composed of species in armillaria subgenus desarmillaria.
 armillaria is a globally distributed mushroom-forming genus composed primarily of plant pathogens. species in this genus are prolific producers of rhizomorphs, or vegetative structures, which, when found, are often associated with infection. because of their importance as plant pathogens, understanding the evolutionary origins of this genus and how it gained a worldwide distribution is of interest. the first gasteroid fungus with close affinities to armillaria—guyanagaster necrorhizus—was described from the neotropical rainforests of guyana. in this study, we conducted phylogenetic analyses to fully resolve the relationship of g. necrorhizus with armillaria. data sets containing guyanagaster from two collecting localities, along with a global sampling of armillaria species—including newly collected specimens from guyana and africa—at six loci were used. three loci—28s, ef1α and rpb2—were analyzed in a partitioned nucleotide data set to infer divergence dates and ancestral range estimations for well-supported, monophyletic lineages.  the six-locus phylogenetic analysis resolves guyanagaster as the earliest diverging lineage in the armillarioid clade. the next lineage to diverge is that composed of species in armillaria subgenus desarmillaria. this subgenus is elevated to genus level to accommodate the exannulate mushroom-forming armillarioid species. the final lineage to diverge is that composed of annulate mushroom-forming armillarioid species, in what is now armillaria sensu stricto. the molecular clock analysis and ancestral range estimation suggest the most recent common ancestor to the armillarioid lineage arose million years ago in eurasia. a new species, guyanagaster lucianii sp. nov. from guyana, is described.  the armillarioid lineage evolved in eurasia during the height of tropical rainforest expansion about million years ago, a time marked by a warm and wet global climate. species of guyanagaster and desarmillaria represent extant taxa of these early diverging lineages. desarmillaria represents an armillarioid lineage that was likely much more widespread in the past. guyanagaster likely evolved from a gilled mushroom ancestor and could represent a highly specialized endemic in the guiana shield. armillaria species represent those that evolved after the shift in climate from warm and tropical to cool and arid during the late eocene. no species in either desarmillaria or guyanagaster are known to produce melanized rhizomorphs in nature, whereas almost all armillaria species are known to produce them. the production of rhizomorphs is an adaptation to harsh environments, and could be a driver of diversification in armillaria by conferring a competitive advantage to the species that produce them. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords armillaria root rotcamerooneocenefungal taxonomygasteromycetationguiana shieldmelaninmushroom evolutionphysalacriaceaesystematicshttp://dx.doi.org/ /100000001national science foundationdeb-0918591aime m. catherine issue-copyright-statement© the author 2017 