Phylogenomics reveals extensive misidentification of fungal strains from the genus Aspergillus.

Modern taxonomic classification is often based on phylogenetic analyses of a few molecular markers, although single-gene studies are still common. Here, we leverage genome-scale molecular phylogenetics (phylogenomics) of species and populations to reconstruct evolutionary relationships in a dense data set of 710 fungal genomes from the biomedically and technologically important genus Aspergillus. To do so, we generated a novel set of 1,362 high-quality molecular markers specific for Aspergillus and provided profile Hidden Markov Models for each, facilitating their use by others. Examining the resulting phylogeny helped resolve ongoing taxonomic controversies, identified new ones, and revealed extensive strain misidentification (7.59% of strains were previously misidentified), underscoring the importance of population-level sampling in species classification. These findings were corroborated using the current standard, taxonomically informative loci. These findings suggest that phylogenomics of species and populations can facilitate accurate taxonomic classifications and reconstructions of the Tree of Life.IMPORTANCEIdentification of fungal species relies on the use of molecular markers. Advances in genomic technologies have made it possible to sequence the genome of any fungal strain, making it possible to use genomic data for the accurate assignment of strains to fungal species (and for the discovery of new ones). We examined the usefulness and current limitations of genomic data using a large data set of 710 publicly available genomes from multiple strains and species of the biomedically, agriculturally, and industrially important genus Aspergillus. Our evolutionary genomic analyses revealed that nearly 8% of publicly available Aspergillus genomes are misidentified. Our work highlights the usefulness of genomic data for fungal systematic biology and suggests that systematic genome sequencing of multiple strains, including reference strains (e.g., type strains), of fungal species will be required to reduce misidentification errors in public databases.

Sulfur assimilation using gaseous carbonyl sulfide by the soil fungus Trichoderma harzianum.

Fungi have the capacity to assimilate a diverse range of both inorganic and organic sulfur compounds. It has been recognized that all sulfur sources taken up by fungi are in soluble forms. In this study, we present evidence that fungi can utilize gaseous carbonyl sulfide (COS) for the assimilation of a sulfur compound. We found that the filamentous fungus Trichoderma harzianum strain THIF08, which has constitutively high COS-degrading activity, was able to grow with COS as the sole sulfur source. Cultivation with 34S-labeled COS revealed that sulfur atom from COS was incorporated into intracellular metabolites such as glutathione and ergothioneine. COS degradation by strain THIF08, in which as much of the moisture derived from the agar medium as possible was removed, indicated that gaseous COS was taken up directly into the cell. Escherichia coli transformed with a COS hydrolase (COSase) gene, which is clade D of the ß-class carbonic anhydrase subfamily enzyme with high specificity for COS but low activity for CO2 hydration, showed that the COSase is involved in COS assimilation. Comparison of sulfur metabolites of strain THIF08 revealed a higher relative abundance of reduced sulfur compounds under the COS-supplemented condition than the sulfate-supplemented condition, suggesting that sulfur assimilation is more energetically efficient with COS than with sulfate because there is no redox change of sulfur. Phylogenetic analysis of the genes encoding COSase, which are distributed in a wide range of fungal taxa, suggests that the common ancestor of Ascomycota, Basidiomycota, and Mucoromycota acquired COSase at about 790-670 Ma.IMPORTANCEThe biological assimilation of gaseous CO2 and N2 involves essential processes known as carbon fixation and nitrogen fixation, respectively. In this study, we found that the fungus Trichoderma harzianum strain THIF08 can grow with gaseous carbonyl sulfide (COS), the most abundant and ubiquitous gaseous sulfur compound, as a sulfur source. When the fungus grew in these conditions, COS was assimilated into sulfur metabolites, and the key enzyme of this assimilation process is COS hydrolase (COSase), which specifically degrades COS. Moreover, the pathway was more energy efficient than the typical sulfate assimilation pathway. COSase genes are widely distributed in Ascomycota, Basidiomycota, and Mucoromycota and also occur in some Chytridiomycota, indicating that COS assimilation is widespread in fungi. Phylogenetic analysis of these genes revealed that the acquisition of COSase in filamentous fungi was estimated to have occurred at about 790-670 Ma, around the time that filamentous fungi transitioned to a terrestrial environment.

Transcriptomics of Temporal- versus Substrate-Specific Wood Decay in the Brown-Rot Fungus Fibroporia radiculosa.

Brown-rot fungi lack many enzymes associated with complete wood degradation, such as lignin-attacking peroxidases, and have developed alternative mechanisms for rapid wood breakdown. To identify the effects of culture conditions and wood substrates on gene expression, we grew Fibroporia radiculosa in submerged cultures containing Wiley milled wood (5 days) and solid wood wafers (30 days), using aspen, pine, and spruce as a substrate. The comparative analysis revealed that wood species had a limited effect on the transcriptome: <3% of genes were differentially expressed between different wood species substrates. The comparison between gene expression during growth on milled wood and wood wafer conditions, however, indicated that the genes encoding plant cell wall-degrading enzymes, such as glycoside hydrolases and peptidases, were activated during growth on wood wafers, confirming previous reports. On the other hand, it was shown for the first time that the genes encoding Fenton chemistry enzymes, such as hydroquinone biosynthesis enzymes and oxidoreductases, were activated during submerged growth on ground wood. This illustrates the diversity of wood-decay reactions encoded in fungi and activated at different stages of this process.

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.

Proteome of the Wood Decay Fungus Fomitopsis pinicola Is Altered by Substrate.

The brown rot fungus Fomitopsis pinicola efficiently depolymerizes wood cellulose via the combined activities of oxidative and hydrolytic enzymes. Mass spectrometric analyses of culture filtrates identified specific proteins, many of which were differentially regulated in response to substrate composition.

Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes.

Introduction: Mushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi. Objectives: This study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution. Methods: We analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution. Results: We detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation. Conclusion: Adaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

Global phylogeny of the Shiitake mushroom and related Lentinula species uncovers novel diversity and suggests an origin in the Neotropics.

Lentinula (Basidiomycota, Agaricales) includes the most widely cultivated mushroom in the world, Lentinula edodes, also known as shiitake (Japanese) or xiang-gu (Chinese). At present, nine species are recognized in the genus, based on morphology, mating criteria, and geographic distribution. However, analyses of internal transcribed spacers (ITS) of ribosomal RNA genes have suggested that there are cryptic lineages. We analyzed a global-scale phylogenetic dataset from 325 Lentinula individuals from 24 countries in Asia-Australasia and the Americas plus Madagascar, with 325 sequences of ITS, 80 LSU sequences, and 111 sequences of translation elongation factor (tef1-α) genes. We recovered 15 independent lineages (Groups 1-15) that may correspond to species. Lineages in Asia-Australasia (Groups 1-5) and the Americas plus Madagascar (Groups 6-15) formed sister clades. Four lineages are represented only by sequences from single individuals and require further molecular sampling, including L. aff. raphanica (Group 7), L. ixodes (Group 8), L. boryana (Group 12), and L. aff. aciculospora (Group 14). Groups 1 and 5 are here referred to L. edodes and L. aff. edodes, respectively. However, these groups most likely represent the same species and are only recognized as (unsupported) monophyletic lineages by maximum likelihood analyses of ITS alone. Other putative species resolved here include L. lateritia (Group 2), L. novae-zelandieae (Group 3), L. aff. lateritia (Group 4), L. raphanica (Group 6), L. aff. detonsa (Group 9), L. detonsa (Group 10), L. guzmanii sp. nov. (Group 11), L. aciculospora (Group 13), and L. madagasikarensis (Group 15). Groups 9-12 represent the "L. boryana complex". Molecular clock and historical biogeographic analyses suggest that the most recent common ancestor (MRCA) of Lentinula can be placed in the middle Oligocene, ca. 30 million years ago (Ma), and had a likely presence in neotropical America. The MRCA of Lentinula in the Americas and Madagascar lived ca. 22 Ma in the Neotropics and the MRCA of Lentinula in Asia-Australasia lived ca. 6 Ma in Oceania. Given the current knowledge about plate tectonics and paleoclimatic models of the last 30 Myr, our phylogenetic hypothesis suggests that the extant distribution of Lentinula is likely to have arisen, in large part, due to long-distance dispersal. Lentinula collections include at least four dubious taxa that need further taxonomic studies: L. reticeps from the USA (Ohio); L. guarapiensis from Paraguay; Lentinus puiggarii from Brazil (São Paulo); and "L. platinedodes" from Vietnam. Approximately ten of the fifteen Groups are reported on Fagaceae, which appears to be the ancestral substrate of Lentinula.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi.

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Evolutionary Morphogenesis of Sexual Fruiting Bodies in Basidiomycota: Toward a New Evo-Devo Synthesis.

The development of sexual fruiting bodies is one of the most complex morphogenetic processes in fungi. Mycologists have long been fascinated by the morphological and developmental diversity of fruiting bodies; however, evolutionary developmental biology of fungi still lags significantly behind that of animals or plants. Here, we summarize the current state of knowledge on fruiting bodies of mushroom-forming Basidiomycota, focusing on phylogenetic and developmental biology. Phylogenetic approaches have revealed a complex history of morphological transformations and convergence in fruiting body morphologies. Frequent transformations and convergence is characteristic of fruiting bodies in contrast to animals or plants, where main body plans are highly conserved. At the same time, insights into the genetic bases of fruiting body development have been achieved using forward and reverse genetic approaches in selected model systems. Phylogenetic and developmental studies of fruiting bodies have each yielded major advances, but they have produced largely disjunct bodies of knowledge. An integrative approach, combining phylogenetic, developmental, and functional biology, is needed to achieve a true fungal evolutionary developmental biology (evo-devo) synthesis for fungal fruiting bodies.

Retracted and Republished from: "Substrate-Specific Differential Gene Expression and RNA Editing in the Brown Rot Fungus Fomitopsis pinicola"

Wood-decaying fungi tend to have characteristic substrate ranges that partly define their ecological niche. Fomitopsis pinicola is a brown rot species of Polyporales that is reported on 82 species of softwoods and 42 species of hardwoods. We analyzed gene expression levels of F. pinicola from submerged cultures with ground wood powder (sampled at 5 days) or solid wood wafers (sampled at 10 and 30 days), using aspen, pine, and spruce substrates (aspen was used only in submerged cultures). Fomitopsis pinicola expressed similar sets of wood-degrading enzymes typical of brown rot fungi across all culture conditions and time points. Nevertheless, differential gene expression was observed across all pairwise comparisons of substrates and time points. Genes exhibiting differential expression encode diverse enzymes with known or potential function in brown rot decay, including laccase, benzoquinone reductase, aryl alcohol oxidase, cytochrome P450s, and various glycoside hydrolases. Comparing transcriptomes from submerged cultures and wood wafers, we found that culture conditions had a greater impact on global expression profiles than substrate wood species. These findings highlight the need for standardization of culture conditions in studies of gene expression in wood-decaying fungi. IMPORTANCE All species of wood-decaying fungi occur on a characteristic range of substrates (host plants), which may be broad or narrow. Understanding the mechanisms that allow fungi to grow on particular substrates is important for both fungal ecology and applied uses of different feedstocks in industrial processes. We grew the wood-decaying polypore Fomitopsis pinicola on three different wood species­aspen, pine, and spruce­under various culture conditions. We found that F. pinicola is able to modify gene expression (transcription levels) across different substrate species and culture conditions. Many of the genes involved encode enzymes with known or predicted functions in wood decay. This study provides clues to how wood-decaying fungi may adjust their arsenal of decay enzymes to accommodate different host substrates.

Genomic Analysis Enlightens Agaricales Lifestyle Evolution and Increasing Peroxidase Diversity.

As actors of global carbon cycle, Agaricomycetes (Basidiomycota) have developed complex enzymatic machineries that allow them to decompose all plant polymers, including lignin. Among them, saprotrophic Agaricales are characterized by an unparalleled diversity of habitats and lifestyles. Comparative analysis of 52 Agaricomycetes genomes (14 of them sequenced de novo) reveals that Agaricales possess a large diversity of hydrolytic and oxidative enzymes for lignocellulose decay. Based on the gene families with the predicted highest evolutionary rates-namely cellulose-binding CBM1, glycoside hydrolase GH43, lytic polysaccharide monooxygenase AA9, class-II peroxidases, glucose-methanol-choline oxidase/dehydrogenases, laccases, and unspecific peroxygenases-we reconstructed the lifestyles of the ancestors that led to the extant lignocellulose-decomposing Agaricomycetes. The changes in the enzymatic toolkit of ancestral Agaricales are correlated with the evolution of their ability to grow not only on wood but also on leaf litter and decayed wood, with grass-litter decomposers as the most recent eco-physiological group. In this context, the above families were analyzed in detail in connection with lifestyle diversity. Peroxidases appear as a central component of the enzymatic toolkit of saprotrophic Agaricomycetes, consistent with their essential role in lignin degradation and high evolutionary rates. This includes not only expansions/losses in peroxidase genes common to other basidiomycetes but also the widespread presence in Agaricales (and Russulales) of new peroxidases types not found in wood-rotting Polyporales, and other Agaricomycetes orders. Therefore, we analyzed the peroxidase evolution in Agaricomycetes by ancestral-sequence reconstruction revealing several major evolutionary pathways and mapped the appearance of the different enzyme types in a time-calibrated species tree.

Fruiting body form, not nutritional mode, is the major driver of diversification in mushroom-forming fungi.

With ∼36,000 described species, Agaricomycetes are among the most successful groups of Fungi. Agaricomycetes display great diversity in fruiting body forms and nutritional modes. Most have pileate-stipitate fruiting bodies (with a cap and stalk), but the group also contains crust-like resupinate fungi, polypores, coral fungi, and gasteroid forms (e.g., puffballs and stinkhorns). Some Agaricomycetes enter into ectomycorrhizal symbioses with plants, while others are decayers (saprotrophs) or pathogens. We constructed a megaphylogeny of 8,400 species and used it to test the following five hypotheses regarding the evolution of morphological and ecological traits in Agaricomycetes and their impact on diversification: 1) resupinate forms are plesiomorphic, 2) pileate-stipitate forms promote diversification, 3) the evolution of gasteroid forms is irreversible, 4) the ectomycorrhizal (ECM) symbiosis promotes diversification, and 5) the evolution of ECM symbiosis is irreversible. The ancestor of Agaricomycetes was a saprotroph with a resupinate fruiting body. There have been 462 transitions in the examined morphologies, including 123 origins of gasteroid forms. Reversals of gasteroid forms are highly unlikely but cannot be rejected. Pileate-stipitate forms are correlated with elevated diversification rates, suggesting that this morphological trait is a key to the success of Agaricomycetes. ECM symbioses have evolved 36 times in Agaricomycetes, with several transformations to parasitism. Across the entire 8,400-species phylogeny, diversification rates of ectomycorrhizal lineages are no greater than those of saprotrophic lineages. However, some ECM lineages have elevated diversification rates compared to their non-ECM sister clades, suggesting that the evolution of symbioses may act as a key innovation at local phylogenetic scales.

Fungal functional ecology: bringing a trait-based approach to plant-associated fungi.

Fungi play many essential roles in ecosystems. They facilitate plant access to nutrients and water, serve as decay agents that cycle carbon and nutrients through the soil, water and atmosphere, and are major regulators of macro-organismal populations. Although technological advances are improving the detection and identification of fungi, there still exist key gaps in our ecological knowledge of this kingdom, especially related to function. Trait-based approaches have been instrumental in strengthening our understanding of plant functional ecology and, as such, provide excellent models for deepening our understanding of fungal functional ecology in ways that complement insights gained from traditional and -omics-based techniques. In this review, we synthesize current knowledge of fungal functional ecology, taxonomy and systematics and introduce a novel database of fungal functional traits (FunFun ). FunFun is built to interface with other databases to explore and predict how fungal functional diversity varies by taxonomy, guild, and other evolutionary or ecological grouping variables. To highlight how a quantitative trait-based approach can provide new insights, we describe multiple targeted examples and end by suggesting next steps in the rapidly growing field of fungal functional ecology.

Transcriptomic atlas of mushroom development reveals conserved genes behind complex multicellularity in fungi.

The evolution of complex multicellularity has been one of the major transitions in the history of life. In contrast to simple multicellular aggregates of cells, it has evolved only in a handful of lineages, including animals, embryophytes, red and brown algae, and fungi. Despite being a key step toward the evolution of complex organisms, the evolutionary origins and the genetic underpinnings of complex multicellularity are incompletely known. The development of fungal fruiting bodies from a hyphal thallus represents a transition from simple to complex multicellularity that is inducible under laboratory conditions. We constructed a reference atlas of mushroom formation based on developmental transcriptome data of six species and comparisons of >200 whole genomes, to elucidate the core genetic program of complex multicellularity and fruiting body development in mushroom-forming fungi (Agaricomycetes). Nearly 300 conserved gene families and >70 functional groups contained developmentally regulated genes from five to six species, covering functions related to fungal cell wall remodeling, targeted protein degradation, signal transduction, adhesion, and small secreted proteins (including effector-like orphan genes). Several of these families, including F-box proteins, expansin-like proteins, protein kinases, and transcription factors, showed expansions in Agaricomycetes, many of which convergently expanded in multicellular plants and/or animals too, reflecting convergent solutions to genetic hurdles imposed by complex multicellularity among independently evolved lineages. This study provides an entry point to studying mushroom development and complex multicellularity in one of the largest clades of complex eukaryotic organisms.

Evolution of substrate-specific gene expression and RNA editing in brown rot wood-decaying fungi.

Fungi that decay wood have characteristic associations with certain tree species, but the mechanistic bases for these associations are poorly understood. We studied substrate-specific gene expression and RNA editing in six species of wood-decaying fungi from the 'Antrodia clade' (Polyporales, Agaricomycetes) on three different wood substrates (pine, spruce, and aspen) in submerged cultures. We identified dozens to hundreds of substrate-biased genes (i.e., genes that are significantly upregulated in one substrate relative to the other two substrates) in each species, and these biased genes are correlated with their host ranges. Evolution of substrate-biased genes is associated with gene family expansion, gain and loss of genes, and variation in cis- and trans- regulatory elements, rather than changes in protein coding sequences. We also demonstrated widespread RNA editing events in the Antrodia clade, which differ from those observed in the Ascomycota in their distribution, substitution types, and the genomic environment. Moreover, we found that substrates could affect editing positions and frequency, including editing events occurring in mRNA transcribed from wood-decay-related genes. This work shows the extent to which gene expression and RNA editing differ among species and substrates, and provides clues into mechanisms by which wood-decaying fungi may adapt to different hosts.

Contemporaneous radiations of fungi and plants linked to symbiosis.

Interactions between fungi and plants, including parasitism, mutualism, and saprotrophy, have been invoked as key to their respective macroevolutionary success. Here we evaluate the origins of plant-fungal symbioses and saprotrophy using a time-calibrated phylogenetic framework that reveals linked and drastic shifts in diversification rates of each kingdom. Fungal colonization of land was associated with at least two origins of terrestrial green algae and preceded embryophytes (as evidenced by losses of fungal flagellum, ca. 720 Ma), likely facilitating terrestriality through endomycorrhizal and possibly endophytic symbioses. The largest radiation of fungi (Leotiomyceta), the origin of arbuscular mycorrhizae, and the diversification of extant embryophytes occurred ca. 480 Ma. This was followed by the origin of extant lichens. Saprotrophic mushrooms diversified in the Late Paleozoic as forests of seed plants started to dominate the landscape. The subsequent diversification and explosive radiation of Agaricomycetes, and eventually of ectomycorrhizal mushrooms, were associated with the evolution of Pinaceae in the Mesozoic, and establishment of angiosperm-dominated biomes in the Cretaceous.

Genomics and Development of Lentinus tigrinus: A White-Rot Wood-Decaying Mushroom with Dimorphic Fruiting Bodies.

Lentinus tigrinus is a species of wood-decaying fungi (Polyporales) that has an agaricoid form (a gilled mushroom) and a secotioid form (puffball-like, with enclosed spore-bearing structures). Previous studies suggested that the secotioid form is conferred by a recessive allele of a single locus. We sequenced the genomes of one agaricoid (Aga) strain and one secotioid (Sec) strain (39.53-39.88 Mb, with 15,581-15,380 genes, respectively). We mated the Sec and Aga monokaryons, genotyped the progeny, and performed bulked segregant analysis (BSA). We also fruited three Sec/Sec and three Aga/Aga dikaryons, and sampled transcriptomes at four developmental stages. Using BSA, we identified 105 top candidate genes with nonsynonymous SNPs that cosegregate with fruiting body phenotype. Transcriptome analyses of Sec/Sec versus Aga/Aga dikaryons identified 907 differentially expressed genes (DEGs) along four developmental stages. On the basis of BSA and DEGs, the top 25 candidate genes related to fruiting body development span 1.5 Mb (4% of the genome), possibly on a single chromosome, although the precise locus that controls the secotioid phenotype is unresolved. The top candidates include genes encoding a cytochrome P450 and an ATP-dependent RNA helicase, which may play a role in development, based on studies in other fungi.