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. complex multicellularity | evolution | fungi | comparative genomics | fruiting body development F ungi represent a diverse lineage of complex multicellular organisms with a unique evolutionary history compared with complex multicellular animals, embryophytes, florideophytes, and laminarean brown algae (1-4). Within the fungal kingdom, complex multicellularity is discussed mostly in the context of fruiting bodies, which are found in at least eight independent lineages (2), of which the Pezizomycotina (Ascomycota) and the Agaricomycetes (Basidiomycota) contain the vast majority of species. The mushroom-forming fungi (Agaricomycetes) comprise >21,000 species and originated 350 million years ago (5), approximately coinciding with the origin of tetrapods. Fruiting bodies of mushroom-forming fungi have immense importance in agriculture, ecology, and medicine; they represent an important and sustainable food source, with favorable medicinal properties (e.g., antitumor, immunomodulatory) (6). Mushroom-forming fungi share a single origin of fruiting body formation that probably dates to the most recent common ancestor of the Agaricomycetes, Dacrymycetes, and Tremellomycetes (2).Fruiting body development in mushroom-forming fungi has been subject to surprisingly few studies (see, e.g., refs. 7-10), result...

T he genus Armillaria causes root rot disease in both gymnoand angiosperms, in forests, parks, and even vineyards in more than 500 host plant species 1 across the world. Most Armillaria species are facultative necrotrophs, which, after colonizing and killing the root cambium, transition to a saprobic phase, decomposing dead woody tissues of the host. As saprotrophs, Armillaria spp. are white rot (WR) fungi, which can efficiently decompose all components of plant cell walls, including lignin, (hemi-)cellulose and pectin 2 . They produce fleshy fruiting bodies (honey mushrooms) that appear in large clumps around infected plants and produce sexual spores. The vegetative phase of Armillaria is predominantly diploid rather than dikaryotic like most basidiomycetes.Individuals of Armillaria can reach immense sizes and include the 'humongous fungus' , one of the largest terrestrial organisms on Earth 3 , measuring up to 965 hectares and 600 tons 4 , and can display a mutation rate ≅ 3 orders of magnitude lower than most filamentous fungi 5 . Individuals reach this immense size via growing rhizomorphs, dark mycelial strings 1-4 mm wide that allow the fungus to bridge gaps between food sources or host plants 1,6 (hence the name shoestring root rot). Rhizomorphs develop through the aggregation and coordinated parallel growth of hyphae, similar to some fruiting body tissues 7,8 . As migratory and exploratory organs, rhizomorphs can grow approximately 1 m yr −1 and cross several metres underground in search for new hosts, although roles in uptake and longrange translocation of nutrients have also been proposed 1,9,10 . Root contact by rhizomorphs is the main mode of infection by the fungus, which makes the prevention of recurrent infection in Armillariacontaminated areas particularly difficult 1 . Despite their huge impact on forestry, horticulture and agriculture, the genetics of the pathogenicity of Armillaria species is poorly understood. The only -omics data published so far have highlighted a substantial repertoire of plant cell wall degrading enzymes (PCWDE) and secreted proteins, among others, in A. mellea and A. solidipes 11,12 , while analyses of the genomes of other pathogenic basidiomycetes (such as Moniliophthora 13,14 , Heterobasidion 15 and Rhizoctonia 16 ) identified genes coding for PCWDEs, secreted and effector proteins or secondary metabolism (SM) as putative pathogenicity factors. However, the lifecycle and unique dispersal strategy of Armillaria prefigure other evolutionary routes to pathogenicity, which, along with other potential genomic factors (such as transposable elements 17 ) are not yet known.Here, we investigate genome evolution and the origin of pathogenicity in Armillaria using comparative genomics, transcriptomics