Domestication through clandestine cultivation constrained genetic diversity in magic mushrooms relative to naturalized populations

Fungi that are edible or fermentative were domesticated through selective cultivation of their desired traits. Domestication is often associated with inbreeding or selfing, which may fix traits other than those under selection, and causes an overall decrease in heterozygosity. A hallucinogenic mushroom, Psilocybe cubensis, was domesticated from its niche in livestock dung for production of psilocybin. It has caused accidental poisonings since the 1940s in Australia, which is a population hypothesized to be introduced from an unknown center of origin. We sequenced genomes of 38 isolates from Australia and compared them with 86 genomes of commercially available cultivars to determine (1) whether P. cubensis was introduced to Australia, and (2) how domestication has impacted commercial cultivars. Our analyses of genome-wide SNPs and single-copy orthologs showed that the Australian population is naturalized, having recovered its effective population size after a bottleneck when it was introduced, and it has maintained relatively high genetic diversity based on measures of nucleotide and allelic diversity. In contrast, domesticated cultivars generally have low effective population sizes and hallmarks of selfing and clonal propagation, including low genetic diversity, low heterozygosity, high linkage disequilibrium, and low allelic diversity of mating-compatibility genes. Analyses of kinship show that most cultivars are founded from related populations. Alleles in the psilocybin gene cluster are identical across most cultivars of P. cubensis with low diversity across coding sequence; however, unique allelic diversity in Australia and some cultivars may translate to differences in biosynthesis of psilocybin and its analogs.

Cultivation, chemistry, and genome of Psilocybe zapotecorum

Psilocybe zapotecorum is a strongly blue-bruising psilocybin mushroom used by indigenous groups in southeastern Mexico and beyond. While this species has a rich history of ceremonial use, research into its chemistry and genetics have been limited. Herein, we detail mushroom morphology and report on cultivation parameters, chemical profile, and the full genome sequence of P. zapotecorum . First, growth and cloning methods are detailed that are simple, and reproducible. In combination with high resolution microscopic analysis, the strain was barcoded, confirming species-level identification. Full genome sequencing reveals the architecture of the psilocybin gene cluster in P. zapotecorum, and can serve as a reference genome for Psilocybe Clade I. Characterization of the tryptamine profile revealed a psilocybin concentration of 17.9±1.7 mg/g, with a range of 10.6-25.7 mg/g (n=7), and similar tryptamines (psilocin, baeocystin, norbaeocystin, norpsilocin, aeruginascin, 4-HO-tryptamine, and tryptamine) in lesser concentrations for a combined tryptamine concentration of 22.5±3.2 mg/g. These results show P. zapotecorum to be a potent - and variable - Psilocybe mushroom. Chemical profiling, genetic analysis, and cultivation assist in demystifying these mushrooms. As clinical studies with psilocybin gain traction, understanding the diversity of psilocybin mushrooms will assure that psilocybin therapy does not become synonymous with psilocybin mushrooms.

A global phylogenomic analysis of the shiitake genus Lentinula

Lentinula is a broadly distributed group of fungi that contains the cultivated shiitake mushroom, L. edodes. We sequenced 24 genomes representing eight described species and several unnamed lineages of Lentinula from 15 countries on four continents. Lentinula comprises four major clades that arose in the Oligocene, three in the Americas and one in Asia-Australasia. To expand sampling of shiitake mushrooms, we assembled 60 genomes of L. edodes from China that were previously published as raw Illumina reads and added them to our dataset. Lentinula edodes sensu lato (s. lat.) contains three lineages that may warrant recognition as species, one including a single isolate from Nepal that is the sister group to the rest of L. edodes s. lat., a second with 20 cultivars and 12 wild isolates from China, Japan, Korea, and the Russian Far East, and a third with 28 wild isolates from China, Thailand, and Vietnam. Two additional lineages in China have arisen by hybridization among the second and third groups. Genes encoding cysteine sulfoxide lyase (lecsl) and γ-glutamyl transpeptidase (leggt), which are implicated in biosynthesis of the organosulfur flavor compound lenthionine, have diversified in Lentinula. Paralogs of both genes that are unique to Lentinula (lecsl 3 and leggt 5b) are coordinately up-regulated in fruiting bodies of L. edodes. The pangenome of L. edodes s. lat. contains 20,308 groups of orthologous genes, but only 6,438 orthogroups (32%) are shared among all strains, whereas 3,444 orthogroups (17%) are found only in wild populations, which should be targeted for conservation.

Genome sequencing progenies of magic mushrooms (Psilocybe subaeruginosa) identifies tetrapolar mating and gene duplications in the psilocybin pathway

Knowledge of breeding systems and genetic diversity is critical to select and combine desired traits that advance new cultivars in agriculture and horticulture. Mushrooms that produce psilocybin, magic mushrooms, may potentially be used in therapeutic and wellness industries, and stand to benefit from genetic improvement. We studied haploid siblings of Psilocybe subaeruginosa to resolve the genetics behind mating compatibility and advance knowledge of breeding. Our results show that mating in P. subaeruginosa is tetrapolar, with compatibility controlled at a homeodomain locus with one copy each of HD1 and HD2, and a pheromone/receptor locus with four homologs of the receptor gene STE3. An additional two pheromone/receptor loci homologous to STE3 do not appear to regulate mating compatibility. Alleles in the psilocybin gene cluster did not vary among the five siblings and were likely homozygous in the parent. Psilocybe subaeruginosa and its relatives have three copies of PsiH genes but their impact on production of psilocybin and its analogues is unknown. Genetic improvement in Psilocybe will require access to genetic diversity from the centre of origin of different species, identification of genes behind traits, and strategies to avoid inbreeding depression.

A metabolomic platform to identify and quantify polyphenols in coffee and related species using liquid chromatography mass spectrometry

INTRODUCTION

Products of plant secondary metabolism, such as phenolic compounds, flavonoids, alkaloids, and hormones, play an important role in plant growth, development, stress resistance. The plant family Rubiaceae is extremely diverse and abundant in Central America and contains several economically important genera, e.g. Coffea and other medicinal plants. These are known for the production of bioactive polyphenols (e.g. caffeine and quinine), which have had major impacts on human society. The overall goal of this study was to develop a high-throughput workflow to identify and quantify plant polyphenols.

METHODS

First, a method was optimized to extract over 40 families of phytochemicals. Then, a high-throughput metabolomic platform has been developed to identify and quantify 184 polyphenols in 15 min.

RESULTS

The current metabolomics study of secondary metabolites was conducted on leaves from one commercial coffee variety and two wild species that also belong to the Rubiaceae family. Global profiling was performed using liquid chromatography high-resolution time-of-flight mass spectrometry. Features whose abundance was significantly different between coffee species were discriminated using statistical analysis and annotated using spectral databases. The identified features were validated by commercially available standards using our newly developed liquid chromatography tandem mass spectrometry method.

DISCUSSION

Caffeine, trigonelline and theobromine were highly abundant in coffee leaves, as expected. Interestingly, wild Rubiaceae leaves had a higher diversity of phytochemicals in comparison to commercial coffee: defense-related molecules, such as phenylpropanoids (e.g., cinnamic acid), the terpenoid gibberellic acid, and the monolignol sinapaldehyde were found more abundantly in wild Rubiaceae leaves.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

Giant Starship Elements Mobilize Accessory Genes in Fungal Genomes

Accessory genes are variably present among members of a species and are a reservoir of adaptive functions. In bacteria, differences in gene distributions among individuals largely result from mobile elements that acquire and disperse accessory genes as cargo. In contrast, the impact of cargo-carrying elements on eukaryotic evolution remains largely unknown. Here, we show that variation in genome content within multiple fungal species is facilitated by Starships, a newly discovered group of massive mobile elements that are 110 kb long on average, share conserved components, and carry diverse arrays of accessory genes. We identified hundreds of Starship-like regions across every major class of filamentous Ascomycetes, including 28 distinct Starships that range from 27 to 393 kb and last shared a common ancestor ca. 400 Ma. Using new long-read assemblies of the plant pathogen Macrophomina phaseolina, we characterize four additional Starships whose activities contribute to standing variation in genome structure and content. One of these elements, Voyager, inserts into 5S rDNA and contains a candidate virulence factor whose increasing copy number has contrasting associations with pathogenic and saprophytic growth, suggesting Voyager's activity underlies an ecological trade-off. We propose that Starships are eukaryotic analogs of bacterial integrative and conjugative elements based on parallels between their conserved components and may therefore represent the first dedicated agents of active gene transfer in eukaryotes. Our results suggest that Starships have shaped the content and structure of fungal genomes for millions of years and reveal a new concerted route for evolution throughout an entire eukaryotic phylum.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Ecological generalism drives hyperdiversity of secondary metabolite gene clusters in xylarialean endophytes

Although secondary metabolites are typically associated with competitive or pathogenic interactions, the high bioactivity of endophytic fungi in the Xylariales, coupled with their abundance and broad host ranges spanning all lineages of land plants and lichens, suggests that enhanced secondary metabolism might facilitate symbioses with phylogenetically diverse hosts. Here, we examined secondary metabolite gene clusters (SMGCs) across 96 Xylariales genomes in two clades (Xylariaceae s.l. and Hypoxylaceae), including 88 newly sequenced genomes of endophytes and closely related saprotrophs and pathogens. We paired genomic data with extensive metadata on endophyte hosts and substrates, enabling us to examine genomic factors related to the breadth of symbiotic interactions and ecological roles. All genomes contain hyperabundant SMGCs; however, Xylariaceae have increased numbers of gene duplications, horizontal gene transfers (HGTs) and SMGCs. Enhanced metabolic diversity of endophytes is associated with a greater diversity of hosts and increased capacity for lignocellulose decomposition. Our results suggest that, as host and substrate generalists, Xylariaceae endophytes experience greater selection to diversify SMGCs compared with more ecologically specialised Hypoxylaceae species. Overall, our results provide new evidence that SMGCs may facilitate symbiosis with phylogenetically diverse hosts, highlighting the importance of microbial symbioses to drive fungal metabolic diversity.

De Novo Gene Birth, Horizontal Gene Transfer, and Gene Duplication as Sources of New Gene Families Associated with the Origin of Symbiosis in Amanita

By introducing novel capacities and functions, new genes and gene families may play a crucial role in ecological transitions. Mechanisms generating new gene families include de novo gene birth, horizontal gene transfer, and neofunctionalization following a duplication event. The ectomycorrhizal (ECM) symbiosis is a ubiquitous mutualism and the association has evolved repeatedly and independently many times among the fungi, but the evolutionary dynamics enabling its emergence remain elusive. We developed a phylogenetic workflow to first understand if gene families unique to ECM Amanita fungi and absent from closely related asymbiotic species are functionally relevant to the symbiosis, and then to systematically infer their origins. We identified 109 gene families unique to ECM Amanita species. Genes belonging to unique gene families are under strong purifying selection and are upregulated during symbiosis, compared with genes of conserved or orphan gene families. The origins of seven of the unique gene families are strongly supported as either de novo gene birth (two gene families), horizontal gene transfer (four), or gene duplication (one). An additional 34 families appear new because of their selective retention within symbiotic species. Among the 109 unique gene families, the most upregulated gene in symbiotic cultures encodes a 1-aminocyclopropane-1-carboxylate deaminase, an enzyme capable of downregulating the synthesis of the plant hormone ethylene, a common negative regulator of plant-microbial mutualisms.

Diversification of Fungal Chitinases and Their Functional Differentiation in Histoplasma capsulatum

Chitinases enzymatically hydrolyze chitin, a highly abundant and utilized polymer of N-acetyl-glucosamine. Fungi are a rich source of chitinases; however, the phylogenetic and functional diversity of fungal chitinases are not well understood. We surveyed fungal chitinases from 373 publicly available genomes, characterized domain architecture, and conducted phylogenetic analyses of the glycoside hydrolase (GH18) domain. This large-scale analysis does not support the previous division of fungal chitinases into three major clades (A, B, C) as chitinases previously assigned to the “C” clade are not resolved as distinct from the “A” clade. Fungal chitinase diversity was partly shaped by horizontal gene transfer, and at least one clade of bacterial origin occurs among chitinases previously assigned to the “B” clade. Furthermore, chitin-binding domains (including the LysM domain) do not define specific clades, but instead are found more broadly across clades of chitinases. To gain insight into biological function diversity, we characterized all eight chitinases (Cts) from the thermally dimorphic fungus, Histoplasma capsulatum: six A clade, one B clade, and one formerly classified C clade chitinases. Expression analyses showed variable induction of chitinase genes in the presence of chitin but preferential expression of CTS3 in the mycelial stage. Activity assays demonstrated that Cts1 (B-I), Cts2 (A-V), Cts3 (A-V), Cts4 (A-V) have endochitinase activities with varying degrees of chitobiosidase function. Cts6 (C-I) has activity consistent with N-acetyl-glucosaminidase exochitinase function and Cts8 (A-II) has chitobiase activity. These results suggest chitinase activity is variable even within subclades and that predictions of functionality require more sophisticated models.

The Architecture of Metabolism Maximizes Biosynthetic Diversity in the Largest Class of Fungi

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or “specialized” metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

Ecology: Fungal Mimics Dupe Animals by Transforming Plants

Many fungi transform host tissues to benefit their own reproduction. A recent study investigates a fungus that converts its plant host's reproductive tissues into ornate flower mimics. These 'pseudoflowers' present complex cues that may enlist insects to facilitate fungal dispersal.

Repeated gain and loss of a single gene modulates the evolution of vascular plant pathogen lifestyles

Vascular plant pathogens travel long distances through host veins, leading to life-threatening, systemic infections. In contrast, nonvascular pathogens remain restricted to infection sites, triggering localized symptom development. The contrasting features of vascular and nonvascular diseases suggest distinct etiologies, but the basis for each remains unclear. Here, we show that the hydrolase CbsA acts as a phenotypic switch between vascular and nonvascular plant pathogenesis. cbsA was enriched in genomes of vascular phytopathogenic bacteria in the family Xanthomonadaceae and absent in most nonvascular species. CbsA expression allowed nonvascular Xanthomonas to cause vascular blight, while cbsA mutagenesis resulted in reduction of vascular or enhanced nonvascular symptom development. Phylogenetic hypothesis testing further revealed that cbsA was lost in multiple nonvascular lineages and more recently gained by some vascular subgroups, suggesting that vascular pathogenesis is ancestral. Our results overall demonstrate how the gain and loss of single loci can facilitate the evolution of complex ecological traits.

Dynamics in Secondary Metabolite Gene Clusters in Otherwise Highly Syntenic and Stable Genomes in the Fungal Genus Botrytis

Fungi of the genus Botrytis infect >1,400 plant species and cause losses in many crops. Besides the broad host range pathogen Botrytis cinerea, most other species are restricted to a single host. Long-read technology was used to sequence genomes of eight Botrytis species, mostly pathogenic on Allium species, and the related onion white rot fungus, Sclerotium cepivorum. Most assemblies contained <100 contigs, with the Botrytis aclada genome assembled in 16 gapless chromosomes. The core genome and pan-genome of 16 Botrytis species were defined and the secretome, effector, and secondary metabolite repertoires analyzed. Among those genes, none is shared among all Allium pathogens and absent from non-Allium pathogens. The genome of each of the Allium pathogens contains 8-39 predicted effector genes that are unique for that single species, none stood out as potential determinant for host specificity. Chromosome configurations of common ancestors of the genus Botrytis and family Sclerotiniaceae were reconstructed. The genomes of B. cinerea and B. aclada were highly syntenic with only 19 rearrangements between them. Genomes of Allium pathogens were compared with ten other Botrytis species (nonpathogenic on Allium) and with 25 Leotiomycetes for their repertoire of secondary metabolite gene clusters. The pattern was complex, with several clusters displaying patchy distribution. Two clusters involved in the synthesis of phytotoxic metabolites are at distinct genomic locations in different Botrytis species. We provide evidence that the clusters for botcinic acid production in B. cinerea and Botrytis sinoallii were acquired by horizontal transfer from taxa within the same genus.

Phlox Species Show Quantitative and Qualitative Resistance to a Population of Powdery Mildew Isolates from the Eastern United States

Ornamental plants in the genus Phlox are extensively planted in landscapes and home gardens around the world. A major limitation to a more widespread use of these plants is their susceptibility to powdery mildew (PM). In this study, we used multilocus sequence typing (MLST) analysis to gain insights into the population diversity of 32 Phlox PM pathogen (Golovinomyces magnicellulatus and Podosphaera sp.) isolates collected from the eastern United States and relate it to the ability to overcome host resistance. Low genetic diversity and a lack of structure were found within our population. Whole genome comparison of two isolates was used to support low genetic diversity evidence found with the MLST analysis. Recombination was suggested by the incongruences observed in the six phylogenetic trees generated from the housekeeping genes TEF-1α, CSI, ITS, IGS, H3, and TUB. Contrasting with low genetic diversity, we found high phenotypic diversity when using 10 of the 32 isolates to evaluate host resistance in four different Phlox species (P. paniculata 'Dunbar Creek', P. amoena OPGC 3598, P. glaberrima OPGC 3594, and P. subulata OPGC 4185) using in vitro bioassays. We observed quantitative and qualitative resistance in all Phlox species and a consistent low disease severity in our control, P. paniculata 'Dunbar Creek'. Taken together, the results generated in this study constitute a robust screening of popular Phlox germplasm that can be incorporated into breeding programs for PM resistance and provides significant information on the evolution of PM pathogens.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata

BACKGROUND

Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution.

RESULTS

A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting.

CONCLUSIONS

The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

Psychoactive plant- and mushroom-associated alkaloids from two behavior modifying cicada pathogens

Entomopathogenic fungi routinely kill their hosts before releasing infectious spores, but a few species keep insects alive while sporulating, which enhances dispersal. Transcriptomics- and metabolomics-based studies of entomopathogens with post-mortem dissemination from their parasitized hosts have unraveled infection processes and host responses. However, the mechanisms underlying active spore transmission by Entomophthoralean fungi in living insects remain elusive. Here we report the discovery, through metabolomics, of the plant-associated amphetamine, cathinone, in four Massospora cicadina-infected periodical cicada populations, and the mushroom-associated tryptamine, psilocybin, in annual cicadas infected with Massospora platypediae or Massospora levispora, which likely represent a single fungal species. The absence of some fungal enzymes necessary for cathinone and psilocybin biosynthesis along with the inability to detect intermediate metabolites or gene orthologs are consistent with possibly novel biosynthesis pathways in Massospora. The neurogenic activities of these compounds suggest the extended phenotype of Massospora that modifies cicada behavior to maximize dissemination is chemically-induced.

Chromosome rearrangements shape the diversification of secondary metabolism in the cyclosporin producing fungus Tolypocladium inflatum

BACKGROUND

Genes involved in production of secondary metabolites (SMs) in fungi are exceptionally diverse. Even strains of the same species may exhibit differences in metabolite production, a finding that has important implications for drug discovery. Unlike in other eukaryotes, genes producing SMs are often clustered and co-expressed in fungal genomes, but the genetic mechanisms involved in the creation and maintenance of these secondary metabolite biosynthetic gene clusters (SMBGCs) remains poorly understood.

RESULTS

In order to address the role of genome architecture and chromosome scale structural variation in generating diversity of SMBGCs, we generated chromosome scale assemblies of six geographically diverse isolates of the insect pathogenic fungus Tolypocladium inflatum, producer of the multi-billion dollar lifesaving immunosuppressant drug cyclosporin, and utilized a Hi-C chromosome conformation capture approach to address the role of genome architecture and structural variation in generating intraspecific diversity in SMBGCs. Our results demonstrate that the exchange of DNA between heterologous chromosomes plays an important role in generating novelty in SMBGCs in fungi. In particular, we demonstrate movement of a polyketide synthase (PKS) and several adjacent genes by translocation to a new chromosome and genomic context, potentially generating a novel PKS cluster. We also provide evidence for inter-chromosomal recombination between nonribosomal peptide synthetases located within subtelomeres and uncover a polymorphic cluster present in only two strains that is closely related to the cluster responsible for biosynthesis of the mycotoxin aflatoxin (AF), a highly carcinogenic compound that is a major public health concern worldwide. In contrast, the cyclosporin cluster, located internally on chromosomes, was conserved across strains, suggesting selective maintenance of this important virulence factor for infection of insects.

CONCLUSIONS

This research places the evolution of SMBGCs within the context of whole genome evolution and suggests a role for recombination between chromosomes in generating novel SMBGCs in the medicinal fungus Tolypocladium inflatum.

Horizontal gene cluster transfer increased hallucinogenic mushroom diversity

Secondary metabolites are a heterogeneous class of chemicals that often mediate interactions between species. The tryptophan‐derived secondary metabolite, psilocin, is a serotonin receptor agonist that induces altered states of consciousness. A phylogenetically disjunct group of mushroom‐forming fungi in the Agaricales produce the psilocin prodrug, psilocybin. Spotty phylogenetic distributions of fungal compounds are sometimes explained by horizontal transfer of metabolic gene clusters among unrelated fungi with overlapping niches. We report the discovery of a psilocybin gene cluster in three hallucinogenic mushroom genomes, and evidence for its horizontal transfer between fungal lineages. Patterns of gene distribution and transmission suggest that synthesis of psilocybin may have provided a fitness advantage in the dung and late wood‐decay fungal niches, which may serve as reservoirs of fungal indole‐based metabolites that alter behavior of mycophagous and wood‐eating invertebrates. These hallucinogenic mushroom genomes will serve as models in neurochemical ecology, advancing the (bio)prospecting and synthetic biology of novel neuropharmaceuticals.

Fungal Gene Cluster Diversity and Evolution

Metabolic gene clusters (MGCs) have provided some of the earliest glimpses at the biochemical machinery of yeast and filamentous fungi. MGCs encode diverse genetic mechanisms for nutrient acquisition and the synthesis/degradation of essential and adaptive metabolites. Beyond encoding the enzymes performing these discrete anabolic or catabolic processes, MGCs may encode a range of mechanisms that enable their persistence as genetic consortia; these include enzymatic mechanisms to protect their host fungi from their inherent toxicities, and integrated regulatory machinery. This modular, self-contained nature of MGCs contributes to the metabolic and ecological adaptability of fungi. The phylogenetic and ecological patterns of MGC distribution reflect the broad diversity of fungal life cycles and nutritional modes. While the origins of most gene clusters are enigmatic, MGCs are thought to be born into a genome through gene duplication, relocation, or horizontal transfer, and analyzing the death and decay of gene clusters provides clues about the mechanisms selecting for their assembly. Gene clustering may provide inherent fitness advantages through metabolic efficiency and specialization, but experimental evidence for this is currently limited. The identification and characterization of gene clusters will continue to be powerful tools for elucidating fungal metabolism as well as understanding the physiology and ecology of fungi.

Evolution of a Heavy Metal Homeostasis/Resistance Island Reflects Increasing Copper Stress in Enterobacteria

Copper homeostasis in bacteria is challenged by periodic elevation of copper levels in the environment, arising from both natural sources and human inputs. Several mechanisms have evolved to efflux copper from bacterial cells, including thecus(copper sensing copper efflux system), andpco(plasmid-borne copper resistance system) systems. The genes belonging to these two systems can be physically clustered in a Copper Homeostasis and Silver Resistance Island (CHASRI) on both plasmids and chromosomes in Enterobacteria. Increasing use of copper in agricultural and industrial applications raises questions about the role of human activity in the evolution of novel copper resistance mechanisms. Here we present evidence that CHASRI emerged and diversified in response to copper deposition across aerobic and anaerobic environments. An analysis of diversification rates and a molecular clock model suggest that CHASRI experienced repeated episodes of elevated diversification that could correspond to peaks in human copper production. Phylogenetic analyses suggest that CHASRI originated in a relative ofEnterobacter cloacaeas the ultimate product of sequential assembly of several pre-existing two-gene modules. Once assembled, CHASRI dispersed via horizontal gene transfer within Enterobacteriaceae and also to certain members of Shewanellaceae, where the originalpcomodule was replaced by a divergentpcohomolog. Analyses of copper stress mitigation suggest that CHASRI confers increased resistance aerobically, anaerobically, and during shifts between aerobic and anaerobic environments, which could explain its persistence in facultative anaerobes and emergent enteric pathogens.

Clustering of two genes putatively involved in cyanate detoxification evolved recently and independently in multiple fungal lineages

Fungi that have the enzymes cyanase and carbonic anhydrase show a limited capacity to detoxify cyanate, a fungicide employed by both plants and humans. Here, we describe a novel two-gene cluster that comprises duplicated cyanase and carbonic anhydrase copies, which we name the CCA gene cluster, trace its evolution across Ascomycetes, and examine the evolutionary dynamics of its spread among lineages of the Fusarium oxysporum species complex (hereafter referred to as the FOSC), a cosmopolitan clade of purportedly clonal vascular wilt plant pathogens. Phylogenetic analysis of fungal cyanase and carbonic anhydrase genes reveals that the CCA gene cluster arose independently at least twice and is now present in three lineages, namely Cochliobolus lunatus, Oidiodendron maius, and the FOSC. Genome-wide surveys within the FOSC indicate that the CCA gene cluster varies in copy number across isolates, is always located on accessory chromosomes, and is absent in FOSC’s closest relatives. Phylogenetic reconstruction of the CCA gene cluster in 163 FOSC strains from a wide variety of hosts suggests a recent history of rampant transfers between isolates. We hypothesize that the independent formation of the CCA gene cluster in different fungal lineages and its spread across FOSC strains may be associated with resistance to plant-produced cyanates or to use of cyanate fungicides in agriculture.

The evolution of fungal metabolic pathways

Fungi contain a remarkable range of metabolic pathways, sometimes encoded by gene clusters, enabling them to digest most organic matter and synthesize an array of potent small molecules. Although metabolism is fundamental to the fungal lifestyle, we still know little about how major evolutionary processes, such as gene duplication (GD) and horizontal gene transfer (HGT), have interacted with clustered and non-clustered fungal metabolic pathways to give rise to this metabolic versatility. We examined the synteny and evolutionary history of 247,202 fungal genes encoding enzymes that catalyze 875 distinct metabolic reactions from 130 pathways in 208 diverse genomes. We found that gene clustering varied greatly with respect to metabolic category and lineage; for example, clustered genes in Saccharomycotina yeasts were overrepresented in nucleotide metabolism, whereas clustered genes in Pezizomycotina were more common in lipid and amino acid metabolism. The effects of both GD and HGT were more pronounced in clustered genes than in their non-clustered counterparts and were differentially distributed across fungal lineages; specifically, GD, which was an order of magnitude more abundant than HGT, was most frequently observed in Agaricomycetes, whereas HGT was much more prevalent in Pezizomycotina. The effect of HGT in some Pezizomycotina was particularly strong; for example, we identified 111 HGT events associated with the 15 Aspergillus genomes, which sharply contrasts with the 60 HGT events detected for the 48 genomes from the entire Saccharomycotina subphylum. Finally, the impact of GD within a metabolic category was typically consistent across all fungal lineages, whereas the impact of HGT was variable. These results indicate that GD is the dominant process underlying fungal metabolic diversity, whereas HGT is episodic and acts in a category- or lineage-specific manner. Both processes have a greater impact on clustered genes, suggesting that metabolic gene clusters represent hotspots for the generation of fungal metabolic diversity.

Ecology drives the distribution of specialized tyrosine metabolism modules in fungi

Gene clusters encoding accessory or environmentally specialized metabolic pathways likely play a significant role in the evolution of fungal genomes. Two such gene clusters encoding enzymes associated with the tyrosine metabolism pathway (KEGG #00350) have been identified in the filamentous fungus Aspergillus fumigatus. The l-tyrosine degradation (TD) gene cluster encodes a functional module that facilitates breakdown of the phenolic amino acid, l-tyrosine through a homogentisate intermediate, but is also involved in the production of pyomelanin, a fungal pathogenicity factor. The gentisate catabolism (GC) gene cluster encodes a functional module likely involved in phenolic compound degradation, which may enable metabolism of biphenolic stilbenes in multiple lineages. Our investigation of the evolution of the TD and GC gene clusters in 214 fungal genomes revealed spotty distributions partially shaped by gene cluster loss and horizontal gene transfer (HGT). Specifically, a TD gene cluster shows evidence of HGT between the extremophilic, melanized fungi Exophiala dermatitidis and Baudoinia compniacensis, and a GC gene cluster shows evidence of HGT between Sordariomycete and Dothideomycete grass pathogens. These results suggest that the distribution of specialized tyrosine metabolism modules is influenced by both the ecology and phylogeny of fungal species.

Repeated loss of an anciently horizontally transferred gene cluster in Botrytis

At least five of the six genes of the bikaverin secondary metabolic gene cluster were shown to have undergone horizontal transfer (HGT) from a Fusarium donor to the Botrytis lineage. Of these five, two enzyme-encoding genes are found as pseudogenes in B. cinerea whereas two regulatory genes and the transporter remain intact. To reconstruct the evolutionary events leading to decay of this gene cluster and infer a more precise timing of its transfer, we examined the genomes of nine additional broadly sampled Botrytis species. We found evidence that a Botrytis ancestor acquired the entire gene cluster through an ancient HGT that occurred before the diversification of the genus. During the subsequent evolution and diversification of the genus, four of the 10 genomes appear to have lost the gene cluster, while in the other six the cluster is in various stages of degeneration. Across the Botrytis genomes, the modes of gene decay in the cluster differed between enzyme-encoding genes, which had higher rates of transition to or retention of pseudogenes and were universally inactivated, and regulatory genes (particularly the non-pathway-specific regulator bik4), which more frequently appeared intact. Consistent with these results, the regulatory genes bik4 and bik5 showed stronger evidence of transcriptional expression than other bikaverin genes under multiple conditions in B. cinerea. These results could be explained by pleiotropy in the bikaverin regulatory genes either through rewiring or their interaction with more central pathways or by constraints on the order of gene loss driven by the intrinsic toxicity of the pathway. Our finding that most of the bikaverin pathway genes have been lost or pseudogenized in these Botrytis genomes suggests that the incidence of HGT of gene cluster-encoded metabolic pathways might be higher than what is possible to be inferred from isolated genome analyses.

Fragmentation of an aflatoxin-like gene cluster in a forest pathogen

Plant pathogens use a complex arsenal of weapons, such as toxic secondary metabolites, to invade and destroy their hosts. Knowledge of how secondary metabolite pathways evolved is central to understanding the evolution of host specificity. The secondary metabolite dothistromin is structurally similar to aflatoxins and is produced by the fungal pine pathogen Dothistroma septosporum. Our study focused on dothistromin genes, which are widely dispersed across one chromosome, to determine whether this unusual distributed arrangement evolved from an ancestral cluster. We combined comparative genomics and population genetics approaches to elucidate the origins of the dispersed arrangement of dothistromin genes over a broad evolutionary time-scale at the phylum, class and species levels. Orthologs of dothistromin genes were found in two major classes of fungi. Their organization is consistent with clustering of core pathway genes in a common ancestor, but with intermediate cluster fragmentation states in the Dothideomycetes fungi. Recombination hotspots in a D. septosporum population matched sites of gene acquisition and cluster fragmentation at higher evolutionary levels. The results suggest that fragmentation of a larger ancestral cluster gave rise to the arrangement seen in D. septosporum. We propose that cluster fragmentation may facilitate metabolic retooling and subsequent host adaptation of plant pathogens.

The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes

Wood is a major pool of organic carbon that is highly resistant to decay, owing largely to the presence of lignin. The only organisms capable of substantial lignin decay are white rot fungi in the Agaricomycetes, which also contains non-lignin-degrading brown rot and ectomycorrhizal species. Comparative analyses of 31 fungal genomes (12 generated for this study) suggest that lignin-degrading peroxidases expanded in the lineage leading to the ancestor of the Agaricomycetes, which is reconstructed as a white rot species, and then contracted in parallel lineages leading to brown rot and mycorrhizal species. Molecular clock analyses suggest that the origin of lignin degradation might have coincided with the sharp decrease in the rate of organic carbon burial around the end of the Carboniferous period.

The ASP3 locus in Saccharomyces cerevisiae originated by horizontal gene transfer from Wickerhamomyces

The asparagine degradation pathway in the S288c laboratory strain of Saccharomyces cerevisiae is comprised of genes located at two separate loci. ASP1 is located on chromosome IV and encodes for cytosolic l-asparaginase I, whereas ASP3 contains a gene cluster located on chromosome XII comprised of four identical genes, ASP3-1, ASP3-2, ASP3-3, and ASP3-4, which encode for cell wall-associated l-asparaginase II. Interestingly, the ASP3 locus appears to be only present, in variable copy number, in S. cerevisiae strains isolated from laboratory or industrial environments and is completely absent from the genomes of 128 diverse fungal species. Investigation of the evolutionary history of ASP3 across these 128 genomes as well as across the genomes of 43 S. cerevisiae strains shows that ASP3 likely arose in a S. cerevisiae strain via horizontal gene transfer (HGT) from, or a close relative of, the wine yeast Wickerhamomyces anomalus, which co-occurs with S. cerevisiae in several biotechnological processes. Thus, because the ASP3 present in the S288c laboratory strain of S. cerevisiae is induced in response to nitrogen starvation, its acquisition may have aided yeast adaptation to artificial environments. Our finding that the ASP3 locus in S. cerevisiae originated via HGT further highlights the importance of gene sharing between yeasts in the evolution of their remarkable metabolic diversity.

The evolutionary imprint of domestication on genome variation and function of the filamentous fungus Aspergillus oryzae

The domestication of animals, plants, and microbes fundamentally transformed the lifestyle and demography of the human species [1]. Although the genetic and functional underpinnings of animal and plant domestication are well understood, little is known about microbe domestication [2-6]. Here, we systematically examined genome-wide sequence and functional variation between the domesticated fungus Aspergillus oryzae, whose saccharification abilities humans have harnessed for thousands of years to produce sake, soy sauce, and miso from starch-rich grains, and its wild relative A. flavus, a potentially toxigenic plant and animal pathogen [7]. We discovered dramatic changes in the sequence variation and abundance profiles of genes and wholesale primary and secondary metabolic pathways between domesticated and wild relative isolates during growth on rice. Our data suggest that, through selection by humans, an atoxigenic lineage of A. flavus gradually evolved into a "cell factory" for enzymes and metabolites involved in the saccharification process. These results suggest that whereas animal and plant domestication was largely driven by Neolithic "genetic tinkering" of developmental pathways, microbe domestication was driven by extensive remodeling of metabolism.

In silico characterization and molecular evolutionary analysis of a novel superfamily of fungal effector proteins

Most fungal plant pathogens secrete effector proteins during pathogenesis to manipulate their host's defense and promote disease. These are so highly diverse in sequence and distribution, they are essentially considered as species-specific. However, we have recently shown the presence of homologous effectors in fungal species of the Dothideomycetes class. One such example is Ecp2, an effector originally described in the tomato pathogen Cladosporium fulvum but later detected in the plant pathogenic fungi Mycosphaerella fijiensis and Mycosphaerella graminicola as well. Here, using in silico sequence-similarity searches against a database of 135 fungal genomes and GenBank, we extend our queries for homologs of Ecp2 to the fungal kingdom and beyond, and further study their history of diversification. Our analyses show that Ecp2 homologs are members of an ancient and widely distributed superfamily of putative fungal effectors, which we term Hce2 for Homologs of C. fulvum Ecp2. Molecular evolutionary analyses show that the superfamily originated and diversified within the fungal kingdom, experiencing multiple lineage-specific expansions and losses that are consistent with the birth-and-death model of gene family evolution. Newly formed paralogs appear to be subject to diversification early after gene duplication events, whereas at later stages purifying selection acts to preserve diversity and the newly evolved putative functions. Some members of the Hce2 superfamily are fused to fungal Glycoside Hydrolase family 18 chitinases that show high similarity to the Zymocin killer toxin from the dairy yeast Kluyveromyces lactis, suggesting an analogous role in antagonistic interactions. The observed high rates of gene duplication and loss in the Hce2 superfamily, combined with diversification in both sequence and possibly functions within and between species, suggest that Hce2s are involved in adaptation to stresses and new ecological niches. Such findings address the need to rationalize effector biology and evolution beyond the perspective of solely host-microbe interactions.

Horizontal transfer and death of a fungal secondary metabolic gene cluster

A cluster composed of four structural and two regulatory genes found in several species of the fungal genus Fusarium (class Sordariomycetes) is responsible for the production of the red pigment bikaverin. We discovered that the unrelated fungus Botrytis cinerea (class Leotiomycetes) contains a cluster of five genes that is highly similar in sequence and gene order to the Fusarium bikaverin cluster. Synteny conservation, nucleotide composition, and phylogenetic analyses of the cluster genes indicate that the B. cinerea cluster was acquired via horizontal transfer from a Fusarium donor. Upon or subsequent to the transfer, the B. cinerea gene cluster became inactivated; one of the four structural genes is missing, two others are pseudogenes, and the fourth structural gene shows an accelerated rate of nonsynonymous substitutions along the B. cinerea lineage, consistent with relaxation of selective constraints. Interestingly, the bik4 regulatory gene is still intact and presumably functional, whereas bik5, which is a pathway-specific regulator, also shows a mild but significant acceleration of evolutionary rate along the B. cinerea lineage. This selective preservation of the bik4 regulator suggests that its conservation is due to its likely involvement in other non-bikaverin-related biological processes in B. cinerea. Thus, in addition to novel metabolism, horizontal transfer of wholesale metabolic gene clusters might also be contributing novel regulation.

Horizontal transfer of a large and highly toxic secondary metabolic gene cluster between fungi

Genes involved in intermediary and secondary metabolism in fungi are frequently physically linked or clustered. For example, in Aspergillus nidulans the entire pathway for the production of sterigmatocystin (ST), a highly toxic secondary metabolite and a precursor to the aflatoxins (AF), is located in a ∼54 kb, 23 gene cluster. We discovered that a complete ST gene cluster in Podospora anserina was horizontally transferred from Aspergillus. Phylogenetic analysis shows that most Podospora cluster genes are adjacent to or nested within Aspergillus cluster genes, although the two genera belong to different taxonomic classes. Furthermore, the Podospora cluster is highly conserved in content, sequence, and microsynteny with the Aspergillus ST/AF clusters and its intergenic regions contain 14 putative binding sites for AflR, the transcription factor required for activation of the ST/AF biosynthetic genes. Examination of ∼52,000 Podospora expressed sequence tags identified transcripts for 14 genes in the cluster, with several expressed at multiple life cycle stages. The presence of putative AflR-binding sites and the expression evidence for several cluster genes, coupled with the recent independent discovery of ST production in Podospora [1], suggest that this HGT event probably resulted in a functional cluster. Given the abundance of metabolic gene clusters in fungi, our finding that one of the largest known metabolic gene clusters moved intact between species suggests that such transfers might have significantly contributed to fungal metabolic diversity. PAPERFLICK:

Horizontal transfer of a nitrate assimilation gene cluster and ecological transitions in fungi: a phylogenetic study

High affinity nitrate assimilation genes in fungi occur in a cluster (fHANT-AC) that can be coordinately regulated. The clustered genes include nrt2, which codes for a high affinity nitrate transporter; euknr, which codes for nitrate reductase; and NAD(P)H-nir, which codes for nitrite reductase. Homologs of genes in the fHANT-AC occur in other eukaryotes and prokaryotes, but they have only been found clustered in the oomycete Phytophthora (heterokonts). We performed independent and concatenated phylogenetic analyses of homologs of all three genes in the fHANT-AC. Phylogenetic analyses limited to fungal sequences suggest that the fHANT-AC has been transferred horizontally from a basidiomycete (mushrooms and smuts) to an ancestor of the ascomycetous mold Trichoderma reesei. Phylogenetic analyses of sequences from diverse eukaryotes and eubacteria, and cluster structure, are consistent with a hypothesis that the fHANT-AC was assembled in a lineage leading to the oomycetes and was subsequently transferred to the Dikarya (Ascomycota+Basidiomycota), which is a derived fungal clade that includes the vast majority of terrestrial fungi. We propose that the acquisition of high affinity nitrate assimilation contributed to the success of Dikarya on land by allowing exploitation of nitrate in aerobic soils, and the subsequent transfer of a complete assimilation cluster improved the fitness of T. reesei in a new niche. Horizontal transmission of this cluster of functionally integrated genes supports the "selfish operon" hypothesis for maintenance of gene clusters.

Diversification of NRT2 and the Origin of Its Fungal Homolog

We investigated the origin and diversification of the high-affinity nitrate transporter NRT2 in fungi and other eukaryotes using Bayesian and maximum parsimony methods. To assess the higher-level relationships and origins of NRT2 in eukaryotes, we analyzed 200 amino acid sequences from the Nitrate/Nitrite Porter (NNP) Family (to which NRT2 belongs), including 55 fungal, 41 viridiplantae (green plants), 11 heterokonts (stramenopiles), and 87 bacterial sequences. To assess evolution of NRT2 within fungi and other eukaryotes, we analyzed 116 amino acid sequences of NRT2 from 58 fungi, 40 viridiplantae (green plants), 1 rhodophyte, and 5 heterokonts, rooted with 12 bacterial sequences. Our results support a single origin of eukaryotic NRT2 from 1 of several clades of mostly proteobacterial NNP transporters. The phylogeny of bacterial NNP transporters does not directly correspond with bacterial taxonomy, apparently due to ancient duplications and/or horizontal gene transfer events. The distribution of NRT2 in the eukaryotes is patchy, but the NRT2 phylogeny nonetheless supports the monophyly of major groups such as viridiplantae, flowering plants, monocots, and eudicots, as well as fungi, ascomycetes, basidiomycetes, and agaric mushrooms. At least 1 secondary origin of eukaryotic NRT2 via horizontal transfer to the fungi is suggested, possibly from a heterokont donor. Our analyses also suggest that there has been a horizontal transfer of nrt2 from a basidiomycete fungus to an ascomycete fungus and reveal a duplication of nrt2 in the ectomycorrhizal mushroom genus, Hebeloma.

Major clades of Agaricales: a multilocus phylogenetic overview

An overview of the phylogeny of the Agaricales is presented based on a multilocus analysis of a six-gene region supermatrix. Bayesian analyses of 5611 nucleotide characters of rpb1, rpb1-intron 2, rpb2 and 18S, 25S, and 5.8S ribosomal RNA genes recovered six major clades, which are recognized informally and labeled the Agaricoid, Tricholomatoid, Marasmioid, Pluteoid, Hygrophoroid and Plicaturopsidoid clades. Each clade is discussed in terms of key morphological and ecological traits. At least 11 origins of the ectomycorrhizal habit appear to have evolved in the Agaricales, with possibly as many as nine origins in the Agaricoid plus Tricholomatoid clade alone. A family-based phylogenetic classification is sketched for the Agaricales, in which 30 families, four unplaced tribes and two informally named clades are recognized.

Reconstructing the early evolution of Fungi using a six-gene phylogeny

The ancestors of fungi are believed to be simple aquatic forms with flagellated spores, similar to members of the extant phylum Chytridiomycota (chytrids). Current classifications assume that chytrids form an early-diverging clade within the kingdom Fungi and imply a single loss of the spore flagellum, leading to the diversification of terrestrial fungi. Here we develop phylogenetic hypotheses for Fungi using data from six gene regions and nearly 200 species. Our results indicate that there may have been at least four independent losses of the flagellum in the kingdom Fungi. These losses of swimming spores coincided with the evolution of new mechanisms of spore dispersal, such as aerial dispersal in mycelial groups and polar tube eversion in the microsporidia (unicellular forms that lack mitochondria). The enigmatic microsporidia seem to be derived from an endoparasitic chytrid ancestor similar to Rozella allomycis, on the earliest diverging branch of the fungal phylogenetic tree.

Contributions of rpb2 and tef1 to the phylogeny of mushrooms and allies (Basidiomycota, Fungi)

A phylogeny of the fungal phylum Basidiomycota is presented based on a survey of 160 taxa and five nuclear genes. Two genes, rpb2, and tef1, are presented in detail. The rpb2 gene is more variable than tef1 and recovers well-supported clades at shallow and deep taxonomic levels. The tef1 gene recovers some deep and ordinal-level relationships but with greater branch support from nucleotides compared to amino acids. Intron placement is dynamic in tef1, often lineage-specific, and diagnostic for many clades. Introns are fewer in rpb2 and tend to be highly conserved by position. When both protein-coding loci are combined with sequences of nuclear ribosomal RNA genes, 18 inclusive clades of Basidiomycota are strongly supported by Bayesian posterior probabilities and 16 by parsimony bootstrapping. These numbers are greater than produced by single genes and combined ribosomal RNA gene regions. Combination of nrDNA with amino acid sequences, or exons with third codon positions removed, produces strong measures of support, particularly for deep internodes of Basidiomycota, which have been difficult to resolve with confidence using nrDNA data alone. This study produces strong boostrap support and significant posterior probabilities for the first time for the following monophyletic groups: (1) Ustilaginomycetes plus Hymenomycetes, (2) an inclusive cluster of hymenochaetoid, corticioid, polyporoid, Thelephorales, russuloid, athelioid, Boletales, and euagarics clades, (3) Thelephorales plus the polyporoid clade, (4) the polyporoid clade, and (5) the cantharelloid clade. Strong support is also recovered for the basal position of the Dacrymycetales in the Hymenomycetidae and paraphyly of the Exobasidiomycetidae.

New Asian species of the genus Anamika (euagarics, hebelomatoid clade) based on morphology and ribosomal DNA sequences

Two dark-spored agaric species from Asia are placed in the genus Anamika (Agaricales or euagarics clade). This result is supported by ITS and nLSU-rDNA sequences with strong measures of branch support, in addition to several morphological and ecological similarities. An inclusive ITS study was performed using a mixed model Bayesian analysis that suggests the derived status of Anamika within Hebeloma, thereby rendering Hebeloma a paraphyletic genus. However, the monophyly of Hebeloma cannot be rejected outright given ITS and nLSU-rDNA data. Thus, we propose two new Asian species in Anamika: A. angustilamellata sp. nov. from dipterocarp and fagaceous forests of southwestern China and northern Thailand; and A. lactariolens comb. nov., a Japanese species originally described in the genus Alnicola. A complete description of A. angustilamellata, including illustrations, is provided.