Mating-Type Locus Organization and Mating-Type Chromosome Differentiation in the Bipolar Edible Button Mushroom Agaricus bisporus

Multiple independent losses of crossover interference during yeast evolutionary history

Meiotic recombination is essential for the accurate chromosome segregation and the generation of genetic diversity through crossover and gene conversion events. Although this process has been studied extensively in a few selected model species, understanding how its properties vary across species remains limited. For instance, the ancestral ZMM pathway that generates interference-dependent crossovers has undergone multiple losses throughout evolution, suggesting variations in the regulation of crossover formation. In this context, we first characterized the meiotic recombination landscape and properties of the Kluyveromyces lactis budding yeast. We then conducted a comprehensive analysis of 29,151 recombination events (19, 212 COs and 9, 939 NCOs) spanning 577 meioses in the five budding yeast species Saccharomyces cerevisiae, Saccharomyces paradoxus, Lachancea kluyveri, Lachancea waltii and K. lactis. Eventually, we found that the Saccharomyces yeasts displayed higher recombination rates compared to the non-Saccharomyces yeasts. In addition, bona fide crossover interference and associated crossover homeostasis were detected in the Saccharomyces species only, adding L. kluyveri and K. lactis to the list of budding yeast species that lost crossover interference. Finally, recombination hotspots, although highly conserved within the Saccharomyces yeasts are not conserved beyond the Saccharomyces genus. Overall, these results highlight great variability in the recombination landscape and properties through budding yeasts evolution.

Why do sex chromosomes progressively lose recombination?

Progressive recombination loss is a common feature of sex chromosomes. Yet, the evolutionary drivers of this phenomenon remain a mystery. For decades, differences in trait optima between sexes (sexual antagonism) have been the favoured hypothesis, but convincing evidence is lacking. Recent years have seen a surge of alternative hypotheses to explain progressive extensions and maintenance of recombination suppression: neutral accumulation of sequence divergence, selection of nonrecombining fragments with fewer deleterious mutations than average, sheltering of recessive deleterious mutations by linkage to heterozygous alleles, early evolution of dosage compensation, and constraints on recombination restoration. Here, we explain these recent hypotheses and dissect their assumptions, mechanisms, and predictions. We also review empirical studies that have brought support to the various hypotheses.

Decoding the chromosome-scale genome of the nutrient-rich Agaricus subrufescens: a resource for fungal biology and biotechnology

Agaricus subrufescens, also known as the "sun mushroom," has significant nutritional and medicinal value. However, its short shelf life due to the browning process results in post-harvest losses unless it's quickly dehydrated. This restricts its availability to consumers in the form of capsules. A genome sequence of A. subrufescens may lead to new cultivation alternatives or the application of gene editing strategies to delay the browning process. We assembled a chromosome-scale genome using a hybrid approach combining Illumina and Nanopore sequencing. The genome was assembled into 13 chromosomes and 31 unplaced scaffolds, totaling 44.5 Mb with 96.5% completeness and 47.24% GC content. 14,332 protein-coding genes were identified, with 64.6% of the genome covered by genes and 23.41% transposable elements. The mitogenome was circularized and encoded fourteen typical mitochondrial genes. Four polyphenol oxidase (PPO) genes and the Mating-type locus were identified. Phylogenomic analysis supports the placement of A. subrufescens in the Agaricomycetes clade. This is the first available genome sequence of a strain of the "sun mushroom." Results are available through a Genome Browser (https://plantgenomics.ncc.unesp.br/gen.php?id=Asub) and can support further fungal biological and genomic studies.

Transcription factors: switches for regulating growth and development in macrofungi

Macrofungi (or mushrooms) act as an extraordinarily important part to human health due to their nutritional and/or medicinal value, but the detailed researches in growth and development mechanisms have yet to be explored further. Transcription factors (TFs) play indispensable roles in signal transduction and affect growth, development, and metabolism of macrofungi. In recent years, increasing research effort has been employed to probe the relationship between the development of macrofungi and TFs. Herein, the present review comprehensively summarized the functional TFs researched in macrofungi, including modulating mycelial growth, fructification, sclerotial formation, sexual reproduction, spore formation, and secondary metabolism. Meanwhile, the possible effect mechanisms of TFs on the growth and development of some macrofungi were also revealed. Specific examples of functional characterizations of TFs in macrofungi (such as Schizophyllum commune and Coprinopsis cinerea) were described to a better comprehension of regulatory effect. Future research prospects in the field of TFs of macrofungi are discussed. We illustrated the functional versatility of the TFs in macrofungi based on specific examples. A systematical realization of the interaction and possible mechanisms between TFs and macrofungi can supply possible solutions to regulate genetic characteristics, which supply novel insights into the regulation of growth, development and metabolism of macrofungi. KEY POINTS: • The functional TFs researched in macrofungi were summarized. • The possible effect mechanisms of TFs in macrofungal were described. • The multiple physiological functions of TFs in macrofungi were discussed.

Genetic and Molecular Evidence of a Tetrapolar Mating System in the Edible Mushroom Grifola frondosa

Grifola frondosa is a valuable edible fungus with high nutritional and medicinal values. The mating systems of fungi not only offer practical strategies for breeding, but also have far-reaching effects on genetic variability. Grifola frondosa has been considered as a sexual species with a tetrapolar mating system based on little experimental data. In the present study, one group of test crosses and six groups of three-round mating experiments from two parental strains were conducted to determine the mating system in G. frondosa. A chi-squared test of the results of the test-cross mating experiments indicated that they satisfied Mendelian segregation, while a series of three-round mating experiments showed that Mendelian segregation was not satisfied, implying a segregation distortion phenomenon in G. frondosa. A genomic map of the G. frondosa strain, y59, grown from an LMCZ basidiospore, with 40.54 Mb and 12 chromosomes, was generated using genome, transcriptome and Hi-C sequencing technology. Based on the genomic annotation of G. frondosa, the mating-type loci A and B were located on chromosomes 1 and 11, respectively. The mating-type locus A coded for the β-fg protein, HD1, HD2 and MIP, in that order. The mating-type locus B consisted of six pheromone receptors (PRs) and five pheromone precursors (PPs) in a crossed order. Moreover, both HD and PR loci may have only one sublocus that determines the mating type in G. frondosa. The nonsynonymous SNP and indel mutations between the A1B1 and A2B2 mating-type strains and the reference genome of y59 only occurred on genes HD2 and PR1/2, preliminarily confirming that the mating type of the y59 strain was A1B2 and not A1B1. Based on the genetic evidence and the more reliable molecular evidence, the results reveal that the mating system of G. frondosa is tetrapolar. This study has important implications for the genetics and hybrid breeding of G. frondosa.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

Structural Analysis of the A Mating Type Locus and Development of the Mating Type Marker of Agaricus bisporus var. bisporus

Karyotyping in Agaricus bisporus is crucial for both the isolation of homokaryotic strains and the confirmation of dikaryon establishment. For the verification of the karyotype, the A mating type loci of two homokaryotic strains, H39 and H97, were analyzed through comparative sequence analysis. The two loci showed major differences in two sequence regions designated as Region 1 and Region 2. H97 had a putative DNA transposon in Region 1 that had target site duplications (TSDs), terminal inverted repeats (TIRs), and a loop sequence, in contrast to H39, which only had the insertional target sequence. Homologous sequences of the transposon were discovered in the two different chromosomes of H97 and in one of H39, all of which have different TSDs but share high sequence homology in TIR. Region 2 shared three consensus sequences between H97 and H39. However, it was only from H97 that a large insertional sequence of unknown origin was discovered between the first and second consensus sequences. The difference in length in Region 1, employed for the verification of the A mating type, resulted in the successful verification of mating types in the heterokaryotic and homokaryotic strains. This length difference enables the discrimination between homo- and heterokaryotic spores by PCR. The present study suggests that the A mating type locus in A. bisporus H97 has evolved through transposon insertion, allowing the discrimination of the mating type, and thus the nuclear type, between A. bisporus H97 and H39.

Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

Experimental Outcrossing in Agaricus bisporus Revealed a Major and Unexpected Involvement of Airborne Mycelium Fragments

Agaricus bisporus var. bisporus, the button mushroom, has a predominantly pseudohomothallic life cycle. Most of its spores are heterokaryotic and give rise to fertile heterokaryons. However, previous studies have suggested that outcrossing should not be rare in wild populations. In order to discover how outcrossing occurs, we experimentally favored it between aerial propagules of a fruiting donor mycelium and a delayed receiver mycelium that only invaded culture trays. Two donor/receiver pairs were studied, and potentially hybrid basidiomata collected on the receiver trays were analyzed with a mitochondrial marker, two unlinked nuclear CAPS markers, then haplotype markers based on DNA sequences obtained after PCR cloning of the rDNA ITS region and the fruk gene. For one of the two pairs, most basidiomata were hybrids between the donor and the receiver. Genotyping of the hybrids revealed only two genotypes consistent with outcrossing involving airborne mycelium fragments rather than basidiospores. The resident receiver heterokaryon that provided its mitochondria to the hybrid basidiomata is suspected to have had a trophic contribution to their growth and successful fruiting. The high level of heterozygosity and the cultivar introgression previously revealed in wild populations of this pseudohomothallic species may result from outcrossing involving airborne pieces of mycelium.

Onset and stepwise extensions of recombination suppression are common in mating-type chromosomes of Microbotryum anther-smut fungi

Sex chromosomes and mating-type chromosomes can display large genomic regions without recombination. Recombination suppression often extended stepwise with time away from the sex- or mating-type-determining genes, generating evolutionary strata of differentiation between alternative sex or mating-type chromosomes. In anther-smut fungi of the Microbotryum genus, recombination suppression evolved repeatedly, linking the two mating-type loci and extended multiple times in regions distal to the mating-type genes. Here, we obtained high-quality genome assemblies of alternative mating types for four Microbotryum fungi. We found an additional event of independent chromosomal rearrangements bringing the two mating-type loci on the same chromosome followed by recombination suppression linking them. We also found, in a new clade analysed here, that recombination suppression between the two mating-type loci occurred in several steps, with first an ancestral recombination suppression between one of the mating-type locus and its centromere; later, completion of recombination suppression up to the second mating-type locus occurred independently in three species. The estimated dates of recombination suppression between the mating-type loci ranged from 0.15 to 3.58 million years ago. In total, this makes at least nine independent events of linkage between the mating-type loci across the Microbotryum genus. Several mating-type locus linkage events occurred through the same types of chromosomal rearrangements, where similar chromosome fissions at centromeres represent convergence in the genomic changes leading to the phenotypic convergence. These findings further highlight Microbotryum fungi as excellent models to study the evolution of recombination suppression.

Whole-Genome Sequencing and Comparative Genomics Analysis of the Wild Edible Mushroom (Gomphus purpuraceus) Provide Insights into Its Potential Food Application and Artificial Domestication

Gomphus purpuraceus (Iwade) Yokoyama is a species of wild fungi that grows in southwest China, considered an edible and medicinal fungus with potential commercial prospects. However, the detailed mechanisms related to the development of mycelium and the formation of the fruiting body are unclear. To obtain a comprehensive overview of genetic features, whole-genome and comparative genomics analyses of G. purpuraceus were performed. High-quality DNA was extracted from the mycelium, which was isolated from a fresh fruiting body of G. purpuraceus. The DNA sample was subjected to sequencing using Illumina and Oxford Nanopore sequencing platforms. A genome assembly totaling 40.15 Mb in 50 contigs with an N50 length of 2.06 Mb was generated, and 8705 putative predicted genes were found. Subsequently, phylogenetic analysis revealed a close evolutionary relationship between G. purpuraceus and Gomphus bonarii. Moreover, a total of 403 carbohydrate-active enzymes (CAZymes) were identified in G. purpuraceus, which included 147 glycoside hydrolases (GHs), 85 glycosyl transferases (GTs), 8 polysaccharide lyases (PLs), 76 carbohydrate esterases (CEs), 57 auxiliary activities (AAs) and 30 carbohydrate-binding modules (CBMs). Compared with the other 13 fungi (Laccaria bicolor, Russula virescens, Boletus edulis, etc.), the number and distribution of CAZymes in G. purpuraceus were similar to other mycorrhizal fungi. Furthermore, the optimization of culture medium for G. purpuraceus showed the efficient utilization of disaccharides such as sucrose and maltose. The genome of G. purpuraceus provides new insights into its niche, food applications and potential artificial domestication.

Genetic Analyses of Discrete Geographic Samples of a Golden Chanterelle in Canada Reveal Evidence for Recent Regional Differentiation

The wild edible mushroom Cantharellus enelensis is a recently described species of the golden chanterelles found in eastern North America. At present, the genetic diversity and population structure of C. enelensis are not known. In this study, we analyzed a total of 230 fruiting bodies of C. enelensis that were collected from three regions of Canada: near the east and west coasts of Newfoundland (NFLD), with 110 fruiting bodies each, and around Hamilton, Ontario (10 fruiting bodies). Among the 110 fruiting bodies from each coast in NFLD, 10 from 2009 were without specific site information, while 100 sampled in 2010 were from each of five patches separated by at least 100 m from each other. Each fruiting body was genotyped at three microsatellite loci. Among the total 28 multilocus genotypes (MLGs) identified, 2 were shared among all three regions, 4 were shared between 2 of the t3hree regions, and the remaining 22 were each found in only 1 region. Minimal spanning network analyses revealed several region-specific MLG clusters, consistent with geographic specific mutation and expansion. Though the most frequently observed MLGs were shared among local (patch) and regional populations, population genetic analyses revealed that both local and regional geographic separations contributed significantly to the observed genetic variation in the total sample. All three regional populations showed excess heterozygosity; for the eastern NFLD population, we reject the null hypothesis of Hardy–Weinberg equilibrium (HWE) at all three loci. However, the analyses of clone-corrected samples revealed that most loci were in HWE. Together, our results suggest that the three discrete regional populations of C. enelensis were likely colonized from a common refugium since the last ice age. However, the local and regional populations are diverging from each other through mutation, drift, and selection at least partly due to heterozygous advantage.

Chromosome-Wide Characterization of Intragenic Crossover in Shiitake Mushroom, Lentinula edodes

Meiotic crossover plays a critical role in generating genetic variations and is a central component of breeding. However, our understanding of crossover in mushroom-forming fungi is limited. Here, in Lentinula edodes, we characterized the chromosome-wide intragenic crossovers, by utilizing the single-nucleotide polymorphisms (SNPs) datasets of an F₁ haploid progeny. A total of 884 intragenic crossovers were identified in 110 single-spore isolates, the majority of which were closer to transcript start sites. About 71.5% of the intragenic crossovers were clustered into 65 crossover hotspots. A 10 bp motif (GCTCTCGAAA) was significantly enriched in the hotspot regions. Crossover frequencies around mating-type A (MAT-A) loci were enhanced and formed a hotspot in L. edodes. Genome-wide quantitative trait loci (QTLs) mapping identified sixteen crossover-QTLs, contributing 8.5–29.1% of variations. Most of the detected crossover-QTLs were co-located with crossover hotspots. Both cis- and trans-QTLs contributed to the nonuniformity of crossover along chromosomes. On chr2, we identified a QTL hotspot that regulated local, global crossover variation and crossover hotspot in L. edodes. These findings and observations provide a comprehensive view of the crossover landscape in L. edodes, and advance our understandings of conservation and diversity of meiotic recombination in mushroom-forming fungi.