Evolutionary relationships within the fungi: analyses of nuclear small subunit rRNA sequences.

Nucleotide sequences of the small subunit ribosomal RNA (18S) gene were used to investigate evolutionary relationships within the Fungi. The inferred tree topologies are in general agreement with traditional classifications in the following ways: (1) the Chytridiomycota and Zygomycota appear to be basal groups within the Fungi. (2) The Ascomycota and Basidiomycota are a derived monophyletic group. (3) Relationships within the Ascomycota are concordant with traditional orders and divide the hemi- and euascomycetes into distinct lineages. (4) The Basidiomycota is divided between the holobasidiomycetes and phragmobasidiomycetes. Conflicts with traditional classification were limited to weakly supported branches of the tree. Strongly supported relationships were robust to minor changes in alignment, method of analysis, and various weighting schemes. Weighting, either of transversions or by site, did not convincingly improve the status of poorly supported portions of the tree. The rate of variation at particular sites does not appear to be independent of lineage, suggesting that covariation of sites may be an important phenomenon in these genes.     

Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species.

Detailed restriction analyses of many samples often require substantial amounts of time and effort for DNA extraction, restriction digests, Southern blotting, and hybridization. We describe a novel approach that uses the polymerase chain reaction (PCR) for rapid simplified restriction typing and mapping of DNA from many different isolates. DNA fragments up to 2 kilobase pairs in length were efficiently amplified from crude DNA samples of several pathogenic Cryptococcus species, including C. neoformans, C. albidus, C. laurentii, and C. uniguttulatus. Digestion and electrophoresis of the PCR products by using frequent-cutting restriction enzymes produced complex restriction phenotypes (fingerprints) that were often unique for each strain or species. We used the PCR to amplify and analyze restriction pattern variation within three major portions of the ribosomal DNA (rDNA) repeats from these fungi. Detailed mapping of many restriction sites within the rDNA locus was determined by fingerprint analysis of progressively larger PCR fragments sharing a common primer site at one end. As judged by PCR fingerprints, the rDNA of 19 C. neoformans isolates showed no variation for four restriction enzymes that we surveyed. Other Cryptococcus spp. showed varying levels of restriction pattern variation within their rDNAs and were shown to be genetically distinct from C. neoformans. The PCR primers used in this study have also been successfully applied for amplification of rDNAs from other pathogenic and nonpathogenic fungi, including Candida spp., and ought to have wide applicability for clinical detection and other studies.     

High diversity and widespread occurrence of mitotic spore mats in ectomycorrhizal Pezizales

Fungal mitospores may function as dispersal units and/ or spermatia and thus play a role in distribution and/or mating of species that produce them. Mitospore production in ectomycorrhizal (EcM) Pezizales is rarely reported, but here we document mitospore production by a high diversity of EcM Pezizales on three continents, in both hemispheres. We sequenced the internal transcribed spacer (ITS) and partial large subunit (LSU) nuclear rDNA from 292 spore mats (visible mitospore clumps) collected in Argentina, Chile, China, Mexico and the USA between 2009 and 2012. We collated spore mat ITS sequences with 105 fruit body and 47 EcM root sequences to generate operational taxonomic units (OTUs). Phylogenetic inferences were made through analyses of both molecular data sets. A total of 48 OTUs from spore mats represented six independent EcM Pezizales lineages and included truffles and cup fungi. Three clades of seven OTUs have no known meiospore stage. Mitospores failed to germinate on sterile media, or form ectomycorrhizas on Quercus, Pinus and Populus seedlings, consistent with a hypothesized role of spermatia. The broad geographic range, high frequency and phylogenetic diversity of spore mats produced by EcM Pezizales suggests that a mitospore stage is important for many species in this group in terms of mating, reproduction and/or dispersal.     

Historical Biogeography and Diversification of Truffles in the Tuberaceae and Their Newly Identified Southern Hemisphere Sister Lineage

Truffles have evolved from epigeous (aboveground) ancestors in nearly every major lineage of fleshy fungi. Because accelerated rates of morphological evolution accompany the transition to the truffle form, closely related epigeous ancestors remain unknown for most truffle lineages. This is the case for the quintessential truffle genus Tuber, which includes species with socio-economic importance and esteemed culinary attributes. Ecologically, Tuber spp. form obligate mycorrhizal symbioses with diverse species of plant hosts including pines, oaks, poplars, orchids, and commercially important trees such as hazelnut and pecan. Unfortunately, limited geographic sampling and inconclusive phylogenetic relationships have obscured our understanding of their origin, biogeography, and diversification. To address this problem, we present a global sampling of Tuberaceae based on DNA sequence data from four loci for phylogenetic inference and molecular dating. Our well-resolved Tuberaceae phylogeny shows high levels of regional and continental endemism. We also identify a previously unknown epigeous member of the Tuberaceae – the South American cup-fungus Nothojafnea thaxteri (E.K. Cash) Gamundí. Phylogenetic resolution was further improved through the inclusion of a previously unrecognized Southern hemisphere sister group of the Tuberaceae. This morphologically diverse assemblage of species includes truffle (e.g. Gymnohydnotrya spp.) and non-truffle forms that are endemic to Australia and South America. Southern hemisphere taxa appear to have diverged more recently than the Northern hemisphere lineages. Our analysis of the Tuberaceae suggests that Tuber evolved from an epigeous ancestor. Molecular dating estimates Tuberaceae divergence in the late Jurassic (∼156 million years ago), with subsequent radiations in the Cretaceous and Paleogene. Intra-continental diversification, limited long-distance dispersal, and ecological adaptations help to explain patterns of truffle evolution and biodiversity.     

Polymorphism at the ribosomal DNA spacers and its relation to breeding structure of the widespread mushroom Schizophyllum commune

The common split-gilled mushroom Schizophyllum commune is found throughout the world on woody substrates. This study addresses the dispersal and population structure of this fungal species by studying the phylogeny and evolutionary dynamics of ribosomal DNA (rDNA) spacer regions. Extensive sampling (n = 195) of sequences of the intergenic spacer region (IGS1) revealed a large number of unique haplotypes (n = 143). The phylogeny of these IGS1 sequences revealed strong geographic patterns and supported three evolutionarily distinct lineages within the global population. The same three geographic lineages were found in phylogenetic analysis of both other rDNA spacer regions (IGS2 and ITS). However, nested clade analysis of the IGS1 phylogeny suggested the population structure of S. commune has undergone recent changes, such as a long distance colonization of western North America from Europe as well as a recent range expansion in the Caribbean. Among all spacer regions, variation in length and nucleotide sequence was observed between but not within the tandem rDNA repeats (arrays). This pattern is consistent with strong within-array and weak among-array homogenizing forces. We present evidence for the suppression of recombination between rDNA arrays on homologous chromosomes that may account for this pattern of concerted evolution.     

Development and characterization of a genetic linkage map of Cryptococcus neoformans var. neoformans using amplified fragment length polymorphisms and other markers

A segregating population of single basidiospore isolates from a sexual cross was used to generate the first moderately dense genetic linkage map of Cryptococcus neoformans var. neoformans (Serotype D). Polymorphic DNA markers were developed using amplified fragment length polymorphisms, random amplified polymorphic DNA, and gene-encoding sequences. These markers were used to analyze 100 meiotic progeny. All markers were tested for distorted segregation with a goodness of fit test. Of the total of 181 markers, 148 showed balanced (1:1) segregation ratios. Segregation distortion was observed for 33 markers. Based on all the markers, a linkage map was generated that consists of 14 major linkage groups with 127 markers, several small linkage groups, and 2 linkage groups that consist only of highly skewed markers. The genetic distance of the linkage map is 1356.3 cM. The estimated total haploid genome size for C. neoformans var. neoformans was calculated using Hulberts method and yielded a map size of 1917 cM. The number of major linkage groups correlates well with the proposed number of 13 chromosomes for C. neoformans var. neoformans. Several genes, including CAP64, CnLAC, and the mating-type locus, were mapped, and their associations were consistent with published data. To date, 6 linkage groups have been assigned to their corresponding chromosomes. This linkage map should provide a framework for the ongoing genome sequencing project and will be a useful tool for studying the genetics and pathogenicity of this important medical yeast.     

Diversity and phylogenetic affinities of foliar fungal endophytes in loblolly pine inferred by culturing and environmental PCR

We examined endophytic fungi in asymptomatic foliage of loblolly pine (Pinus taeda) in North Carolina, U.S.A., with four goals: (i) to evaluate morphotaxa, BLAST matches and groups based on sequence similarity as functional taxonomic units; (ii) to explore methods to maximize phylogenetic signal for environmental datasets, which typically contain many taxa but few characters; (iii) to compare culturing vs. culture-free methods (environmental PCR of surface sterilized foliage) for estimating endophyte diversity and species composition; and (iv) to investigate the relationships between traditional ecological indices (e.g. Shannon index) and phylogenetic diversity (PD) in estimating endophyte diversity and spatial heterogeneity. Endophytes were recovered in culture from 87 of 90 P. taeda leaves sampled, yielding 439 isolates that represented 24 morphotaxa. Sequence data from the nuclear ribosomal internal transcribed spacer (ITS) for 150 isolates revealed 59 distinct ITS genotypes that represented 24 and 37 unique groups based on 90% and 95% sequence similarity, respectively. By recoding ambiguously aligned regions to extract phylogenetic signal and implementing a conservative phylogenetic backbone constraint, we recovered well supported phylogenies based on ca. 600 bp of the nuclear ribosomal large subunit (LSUrDNA) for 72 Ascomycota and Basidiomycota, 145 cultured endophytes and 33 environmental PCR samples. Comparisons with LSUrDNA-delimited species showed that morphotaxa adequately estimated total species richness but rarely corresponded to biologically meaningful groups. ITS BLAST results were variable in their utility, but ITS genotype groups based on 90% sequence similarity were concordant with LSUrDNA-delimited species. Environmental PCR yielded more genotypes per sampling effort and recovered several distinct clades relative to culturing, but some commonly cultured clades were never found (Sordariomycetes) or were rare relative to their high frequency among cultures (Leotiomycetes). In contrast to traditional indices, PD demonstrated spatial heterogeneity in endophyte assemblages among P. taeda trees and study plots. Our results highlight the need for caution in designating taxonomic units based on gross cultural morphology or ITS BLAST matches, the utility of phylogenetic tools for extracting robust phylogenies from environmental samples, the complementarity of culturing and environmental PCR, the utility of PD relative to traditional ecological indices, and the remarkably high diversity of foliar fungal endophytes in this simplified temperate ecosystem.     

Genetic isolation between two recently diverged populations of a symbiotic fungus

Fungi are an omnipresent and highly diverse group of organisms, making up a significant part of eukaryotic diversity. Little is currently known about the drivers of fungal population differentiation and subsequent divergence of species, particularly in symbiotic, mycorrhizal fungi. Here, we investigate the population structure and environmental adaptation in Suillus brevipes (Peck) Kuntze, a wind‐dispersed soil fungus that is symbiotic with pine trees. We assembled and annotated the reference genome for Su. brevipes and resequenced the whole genomes of 28 individuals from coastal and montane sites in California. We detected two clearly delineated coast and mountain populations with very low divergence. Genomic divergence was restricted to few regions, including a region of extreme divergence containing a gene encoding for a membrane Na⁺/H⁺ exchanger known for enhancing salt tolerance in plants and yeast. Our results are consistent with a very recent split between the montane and coastal Su. brevipes populations, with few small genomic regions under positive selection and a pattern of dispersal and/or establishment limitation. Furthermore, we identify a putatively adaptive gene that motivates further functional analyses to link genotypes and phenotypes and shed light on the genetic basis of adaptive traits.     

Responses of soil cellulolytic fungal communities to elevated atmospheric CO₂ are complex and variable across five ecosystems

Elevated atmospheric CO(2) generally increases plant productivity and subsequently increases the availability of cellulose in soil to microbial decomposers. As key cellulose degraders, soil fungi are likely to be one of the most impacted and responsive microbial groups to elevated atmospheric CO(2). To investigate the impacts of ecosystem type and elevated atmospheric CO(2) on cellulolytic fungal communities, we sequenced 10,677 cbhI gene fragments encoding the catalytic subunit of cellobiohydrolase I, across five distinct terrestrial ecosystem experiments after a decade of exposure to elevated CO(2). The cbhI composition of each ecosystem was distinct, as supported by weighted Unifrac analyses (all P-values; < 0.001), with few operational taxonomic units (OTUs) being shared across ecosystems. Using a 114-member cbhI sequence database compiled from known fungi, less than 1% of the environmental sequences could be classified at the family level indicating that cellulolytic fungi in situ are likely dominated by novel fungi or known fungi that are not yet recognized as cellulose degraders. Shifts in fungal cbhI composition and richness that were correlated with elevated CO(2) exposure varied across the ecosystems. In aspen plantation and desert creosote bush soils, cbhI gene richness was significantly higher after exposure to elevated CO(2) (550 μmol mol(-1)) than under ambient CO(2) (360 μmol mol(-1) CO(2)). In contrast, while the richness was not altered, the relative abundance of dominant OTUs in desert soil crusts was significantly shifted. This suggests that responses are complex, vary across different ecosystems and, in at least one case, are OTU-specific. Collectively, our results document the complexity of cellulolytic fungal communities in multiple terrestrial ecosystems and the variability of their responses to long-term exposure to elevated atmospheric CO(2).