Advances in Understanding Mating Type Gene Organization in the Mushroom-Forming Fungus Flammulina velutipes

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

Homeodomain 1 Genes of the Different HD Subloci of Flammulina velutipes Can Activate the HD Pathway and Are Involved in Mating, Clamp Cell Formation, and Upregulation of FvClp1

Flammulina velutipes has two independent and functional mating type factors, HD and PR. The HD locus contains two separate subloci: HD-a and HD-b. In this study, we investigated the roles of Hd1 genes of the HD-a and HD-b subloci in the process of mating, clamp cell formation, and regulation of FvClp1 (F. velutipes clampless1 gene) gene expression in F. velutipes. To this end, we introduced Hd1 genes from mating compatible strains into F. velutipes monokaryon L11. Overexpression of Hd1 gene FvHd-a1-1 of the HD-a sublocus resulted in the formation of pseudoclamps in L11 monokaryons. L11 mutants overexpressing the Hd1 gene FvHd-b1-2 of the HD-b sublocus also similarly developed pseudoclamps in the L11 monokaryons. Moreover, these mutant L11 monokaryons produced complete clamps when crossed with monokaryotic strains that differed at the PR loci, i.e., when selective activation of the PR pathway was obtained through crossing. Thus, Hd1 genes of the two different HD subloci in F. velutipes can activate the HD mating type pathway and induce clamp cell formation. In addition, activation of the HD pathway resulted in upregulation of the FvClp1 gene. Finally, to complete clamp cell formation, activation of the PR pathway appears to be essential. Overall, these findings were beneficial for deepening our understanding of sexual reproduction and fruiting body development of edible fungi.

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.

Structure Analysis of the MatA Locus of Sexual Compatibility in the Edible Mushroom Pleurotus ostreatus

The edible oyster mushroom Pleurotus ostreatus is one of the most cultivated species worldwide. Morphogenesis associated with the maturation of fruit bodies is controlled by two unlinked loci of sexual compatibility matA and matB with multiple alleles (tetrapolar system of sexual compatibility). Quantitative analysis of the alleles of mating compatibility loci in 17 natural isolates collected in the Moscow region was performed in mon-mon (monokaryons-monokaryon) and di-mon (dikaryon-monokaryon) crossings. Four monokaryotic testers strains which were heteroallelic at both mating type loci were obtained for each of the five natural mushroom isolates by using original technique of sterile spore prints on Petri dishes and mon-mon crossing. Twelve natural isolates were crossed via di-mon mating with the four monokaryotic testers M-38. Genetic analysis of the alleles of sexual compatibility loci in 17 natural isolates revealed multiple alleles at both loci: at least ten alleles at matA locus and eight alleles at matB locus. Structural organization analysis of the matA locus was performed in silico for homokaryotic strains PC9 and PC15 based on the whole-genome sequencing data available at DOE Joint Genome Institute. The matA locus has an extremely divergent structure: there are one copy of the homeodomain gene hd1 and one copy of the hd2 gene in the PC9 strain, whereas the matA locus of the PC15 strain is composed by two copies of hd1.1 and hd1.2 genes (class HD1 homeodomain proteins) and one copy of hd2 gene (class HD2 proteins). Comprehensive analysis of amino acid sequences of HD1 and HD2 homeodomain proteins demonstrated that the proteins have a globular structure with the nuclear localization and contain a variable N-terminus and a more conserved DNA-binding domain with a specific conserved motif  WFXNXR in the third ɑ-helix. The results suggest that multiple alleles of the matA locus of sexual compatibility in basidiomycete fungi is achieved due to both different copy number of the coding hd genes within the locus and the variability of the coding gene sequences.

Whole Genome Sequence of an Edible Mushroom Oudemansiella raphanipes (Changgengu)

Oudemansiella raphanipes, considered as a well-known culinary edible mushroom with a high content of natural bioactive substances, is widely cultivated in China with the commercial name Changgengu. However, due to the lack of genomic data, molecular and genetic study on O. raphanipes is rare. To obtain a comprehensive overview of genetic characteristics and enhance the value of O. raphanipes, two mating-compatible monokaryons isolated from the dikaryon were applied for de novo genome sequencing and assembly using Nanopore and /or Illumina sequencing platforms. One of the monokaryons, O. raphanipes CGG-A-s1, was annotated with 21,308 protein-coding genes, of which 56 were predicted to be involved in the biosynthesis of secondary metabolites such as terpene, type I PKS, NRPS, and siderophore. Phylogenetic and comparative analysis of multiple fungi genomes revealed a close evolutionary relationship between O. raphanipes and Mucidula mucid based on single-copy orthologous protein genes. Significant collinearity was detected between O. raphanipes and Flammulina velutipes on the synteny of inter-species genomes. 664 CAZyme genes in CGG-A-s1 were identified with GHs and AAs families significantly elevated when compared with the other 25 sequenced fungi, indicating a strong wood degradation ability. Furthermore, the mating type locus analysis revealed that CGG-A-s1 and CGG-A-s2 were conserved in the gene organization of the mating A locus but various in that of the mating B locus. The genome resource of O. raphanipes will provide new insights into its development of genetic studies and commercial production of high-quality varieties.

Determination and Analysis of Hyper-Variable A Mating Types in Wild Strains of Lentinula edodes in Korea

The diversity of A mating type in wild strains of Lentinula edodes was extensively analyzed to characterize and utilize them for developing new cultivars. One hundred twenty-three A mating type alleles, including 67 newly discovered alleles, were identified from 106 wild strains collected for the past four decades in Korea. Based on previous studies and current findings, a total of 130 A mating type alleles have been found, 124 of which were discovered from wild strains, indicating the hyper-variability of A mating type alleles of L. edodes. About half of the A mating type alleles in wild strains were found in more than two strains, whereas the other half of the alleles were found in only one strain. About 90% of A mating type combinations in dikaryotic wild strains showed a single occurrence. Geographically, diverse A mating type alleles were intensively located in the central region of the Korean peninsula, whereas only allele A17 was observed throughout Korea. We also found the conservation of the TCCCAC motif in addition to the previously reported motifs, including ATTGT, ACAAT, and GCGGAG, in the intergenic regions of A mating loci. Sequence comparison among some alleles indicated that accumulated mutation and recombination would contribute to the diversification of A mating type alleles in L. edodes. Our data support the rapid evolution of A mating locus in L. edodes, and would help to understand the characteristics of A mating loci of wild strains in Korea and help to utilize them for developing new cultivars.

Whole-genome sequence of a high-temperature edible mushroom Pleurotus giganteus (zhudugu)

Most of the sequenced wood-rotting edible mushroom produce fruiting body at relatively low temperatures. Little information has been known about the high-temperature wood-rotting mushroom. Here, we performed de novo sequencing and assembly of the genome of a high-temperature edible mushroom Pleurotus giganteus from a monokaryotic strain zhudugu2 using the Illumina and Pac-Bio CLR sequencing technologies. P. giganteus, also known as Zhudugu in China, is a well-known culinary edible mushroom that has been widely distributed and cultivated in China, Southeast Asia, and South Asia. The genome consists of 40.00 Mb in 27 contigs with a contig N50 of 4.384 Mb. Phylogenetic analysis reveals that P. giganteus and other strains in Pleurotus clustered in one clade. Phylogenetic analysis and average nucleotide identity analysis indicated that the P. giganteus genome showed a closer relationship with other Pleurotus species. Chromosome collinearity analysis revealed a high level of collinearity between P. ostreatus and P. giganteus. There are 12,628 protein-coding genes annotated in this monoploid genome. A total of 481 enzymes accounting for 514 carbohydrate-active enzymes (CAZymes) terms were identified in the P. giganteus genome, including 15 laccases and 10 class II peroxidases predicted in the genome, which revealed the robustness of lignocellulose degradation capacity of P. giganteus. The mating-A type locus of P. giganteus consisted of a pair of homeodomain mating-type genes HD1 and HD2. The mating-B type locus of P. giganteus consisted of at least four pheromone receptor genes and three pheromone genes. The genome is not only beneficial for the genome-assisted breeding of this mushroom but also helps us to understand the high-temperature tolerance of the edible mushroom.

Large-scale fungal strain sequencing unravels the molecular diversity in mating loci maintained by long-term balancing selection

Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of strains. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci define mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty strains of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. The sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. 17,550 mating types were predicted based on the allelic classes. In vitro crosses allowed us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. We conclude that the genetic diversity of mating loci in Trichaptum is due to long-term balancing selection, with limited recombination and duplication activity. The large number of sequenced strains highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Haplotype-Resolved Genome Analyses Reveal Genetically Distinct Nuclei within a Commercial Cultivar of Lentinula edodes

Lentinula edodes is a tetrapolar basidiomycete with two haploid nuclei in each cell during most of their life cycle. Understanding the two haploid nuclei genome structures and their interactions on growth and fruiting body development has significant practical implications, especially for commercial cultivars. In this study, we isolated and assembled the two haploid genomes from a commercial strain of L. edodes using Illumina, HiFi, and Hi-C technologies. The total genome lengths were 50.93 Mb and 49.80 Mb for the two monokaryons SP3 and SP30, respectively, with each assembled into 10 chromosomes with 99.63% and 98.91% anchoring rates, respectively, for contigs more than 100 Kb. Genome comparisons suggest that two haploid nuclei likely derived from distinct genetic ancestries, with ~30% of their genomes being unique or non-syntenic. Consistent with a tetrapolar mating system, the two mating-type loci A (matA) and B (matB) of L. edodes were found located on two different chromosomes. However, we identified a new but incomplete homeodomain (HD) sublocus at ~2.8 Mb from matA in both monokaryons. Our study provides a solid foundation for investigating the relationships among cultivars and between cultivars and wild strains and for studying how two genetically divergent nuclei coordinate to regulate fruiting body formation in L. edodes.

Genetic structure and evolutionary diversity of mating-type (MAT) loci in Hypsizygus marmoreus

The mating compatibility in fungi is generally governed by genes located within a single or two unlinked mating type (MAT) loci. Hypsizygus marmoreus is an edible mushroom in the order Agaricales with a tetrapolar system, which contains two unlinked MAT loci-homeodomain (HD) transcription factor genes and pheromone/pheromone receptor genes (P/R). In this study, we analyzed the genetic structure and diversity of MAT loci in tetrapolar system of H. marmoreus through sequencing of 54 heterokaryon and 8 homokaryon strains. Although within the HD loci, the gene order was conserved, the gene contents were variable, and the HD loci haplotypes were further classified into four types. By analyzing the structure, phylogeny, and the HD transmissibility based on the progeny of these four HD mating-type loci types, we found that they were heritable and tightly linked at the HD loci. The P/R loci genes were found to comprise three pheromone receptors, three pheromones, and two pheromone receptor-like genes. Intra- and inter-specific phylogenetic analyses of pheromone receptors revealed that the STE3 genes were divided into three groups, and we thus theorize that they diverged before speciation. Comparative analysis of the MAT regions among 73 Basidiomycete species indicated that the diversity of HD and P/R loci in Agaricales and Boletales may contribute to mating compatibility. The number of HD genes were not correlated with the tetrapolar or bipolar systems. In H. marmoreus, the expression levels of these genes at HD and P/R loci of compatible strains were found higher than in those of homonuclear/homokaryotic strains, indicating that these mating genes acted as switches for mating processes. Further collinear analysis of HD loci in interspecific species found that HD loci contains conserved recombination hotspots showing major rearrangements in Coprinopsis cinerea and Schizophyllum commune, suggesting different mechanisms for evolution of physically linked MAT loci in these groups. It seems likely that gene rearrangements are common in Agaricales fungi around HD loci. Together, our study provides insights into the genomic basis of mating compatibility in H. marmoreus.

Discovery and Functional Study of a Novel Genomic Locus Homologous to Bα-Mating-Type Sublocus of Lentinula edodes

The interaction of mating pheromone and pheromone receptor from the B mating-type locus is the first step in the activation of the mushroom mating signal transduction pathway. The B mating-type locus of Lentinula edodes is composed of Bα and Bβ subloci, each of which contains genes for mating pheromone and pheromone receptor. Allelic variations in both subloci generate multiple B mating-types through which L. edodes maintains genetic diversity. In addition to the B mating-type locus, our genomic sequence analysis revealed the presence of a novel chromosomal locus 43.3 kb away from the B mating-type locus, containing genes for a pair of mating pheromones (PHBN1 and PHBN2) and a pheromone receptor (RCBN). The new locus (Bα-N) was homologous to the Bα sublocus, but unlike the multiallelic Bα sublocus, it was highly conserved across the wild and cultivated strains. The interactions of RcbN with various mating pheromones from the B and Bα-N mating-type loci were investigated using yeast model that replaced endogenous yeast mating pheromone receptor STE2 with RCBN. The yeast mating signal transduction pathway was only activated in the presence of PHBN1 or PHBN2 in the RcbN producing yeast, indicating that RcbN interacts with self-pheromones (PHBN1 and PHBN2), not with pheromones from the B mating-type locus. The biological function of the Bα-N locus was suggested to control the expression of A mating-type genes, as evidenced by the increased expression of two A-genes HD1 and HD2 upon the treatment of synthetic PHBN1 and PHBN2 peptides to the monokaryotic strain of L. edodes.

The Fvclp1 gene regulates mycelial growth and fruiting body development in edible mushroom Flammulina velutipes

Fruiting body development in Agaricomycetes represents the most complex and unclear process in the fungi. Mating type pathways (A and B) and transcription factors are important regulators in the sexual development of mushrooms. It is known that clampless1 (clp1) is an additional gene that participate under the homeodomain (HD) genes in the matA pathway and clp1 inactivation blocks clamps formation in Coprinopsis cinerea. In this study we identified and analyzed a homologous Fvclp1 gene in the edible mushroom Flammulina velutipes. The coding sequence of the Fvclp1 was 1011 bp without intron interruption, encoding a protein of 336 amino acids. To exhibit the role of Fvclp1 in clamp development and fruiting body formation, knockdown and overexpression mutants were prepared. No significant difference was observed in the monokaryotic hyphal morphology of overexpression and knockdown transformants. In the dikaryotic hyphae from the compatible crossings between the wild-type L22 strain and Fvclp1 knockdown or overexpression mutants, clamp connections developed. However, knockdown mutants could generate fewer fruiting bodies than the wild-type strain. On the contrary, reduced mycelial growth rate but improved fruiting ability was observed in the dikaryotic Fvclp1 overexpression mutants as compared to the wild-type strain. These results indicate that Fvclp1 is necessary and actively involved in fruiting body development in F. velutipes. Overall, these findings suggest that further studies on the function of Fvclp1 would advance our understanding of sexual reproduction and fruiting body development in edible mushrooms.

Molecular and genetic evidence for a tetrapolar mating system in Sparassis latifolia

Sparassis latifolia is a valuable edible fungus cultivated in East Asia that is rich in β-glucans. Understanding the mating system and sexual life cycle is important not only for breeding programs to improve strains but also for studies on speciation and population structures. In the present study, mating experiments using monokaryons derived from two different parental strains were performed. Chi-squared test indicated satisfied Mendel segregation, which supported a tetrapolar mating system. A search in the genome for homologs to the well-defined homeodomain and pheromone/receptors, as well as frequently found flanking genes, resulted in the identification of known mating-type loci previously identified in tetrapolar basidiomycetes, each represented by two idiomorphic alleles on separate contigs. Deficiency of the β-flanking protein in S. latifolia and S. crispa around the MAT-A locus may be explained by the locus being rich in transposable elements adjacent to HD genes. Monokaryotic mycelia are characterized by a slower growth rate and a relative lack of aerial mycelia compared with the parental strain. Chlamydospores can be produced in both monokaryotic and dikaryotic mycelial stages. We provide genetic and molecular evidence for the mating system of S. latifolia, a finding that will be helpful for the cross-breeding of this mushroom.

Transcriptome data reveal conserved patterns of fruiting body development and response to heat stress in the mushroom-forming fungus Flammulina filiformis

Mushroom-forming fungi are complex multicellular organisms that form the basis of a large industry, yet, our understanding of the mechanisms of mushroom development and its responses to various stresses remains limited. The winter mushroom (Flammulina filiformis) is cultivated at a large commercial scale in East Asia and is a species with a preference for low temperatures. This study investigated fruiting body development in F. filiformis by comparing transcriptomes of 4 developmental stages, and compared the developmental genes to a 200-genome dataset to identify conserved genes involved in fruiting body development, and examined the response of heat sensitive and -resistant strains to heat stress. Our data revealed widely conserved genes involved in primordium development of F. filiformis, many of which originated before the emergence of the Agaricomycetes, indicating co-option for complex multicellularity during evolution. We also revealed several notable fruiting-specific genes, including the genes with conserved stipe-specific expression patterns and the others which related to sexual development, water absorption, basidium formation and sporulation, among others. Comparative analysis revealed that heat stress induced more genes in the heat resistant strain (M1) than in the heat sensitive one (XR). Of particular importance are the hsp70, hsp90 and fes1 genes, which may facilitate the adjustment to heat stress in the early stages of fruiting body development. These data highlighted novel genes involved in complex multicellular development in fungi and aid further studies on gene function and efforts to improve the productivity and heat tolerance in mushroom-forming fungi.

Investigation of Mating Pheromone-Pheromone Receptor Specificity in Lentinula edodes

The B mating-type locus of Lentinula edodes, a representative edible mushroom, is highly complex because of allelic variations in the mating pheromone receptors (RCBs) and the mating pheromones (PHBs) in both the Bα and Bβ subloci. The complexity of the B mating-type locus, five Bα subloci with five alleles of RCB1 and nine PHBs and three Bβ subloci with 3 alleles of RCB2 and five PHBs, has led us to investigate the specificity of the PHB-RCB interaction because the interaction plays a key role in non-self-recognition. In this study, the specificities of PHBs to RCB1-2 and RCB1-4 from the Bα sublocus and RCB2-1 from the Bb sublocus were investigated using recombinant yeast strains generated by replacing STE2, an endogenous yeast mating pheromone receptor, with the L. edodes RCBs. Fourteen synthetic PHBs with C-terminal carboxymethylation but without farnesylation were added to the recombinant yeast cells and the PHB-RCB interaction was monitored by the expression of the FUS1 gene-a downstream gene of the yeast mating signal pathway. RCB1-2 (Bα2) was activated by PHB1 (4.3-fold) and PHB2 (2.1-fold) from the Bα1 sublocus and RCB1-4 (Bα4) was activated by PHB5 (3.0-fold) and PHB6 (2.7-fold) from the Bα2 sublocus and PHB13 (3.0-fold) from the Bα5 sublocus. In particular, PHB3 from Bβ2 and PHB9 from Bβ3 showed strong activation of RCB2-1 of the Bβ1 sublocus by 59-fold. The RCB-PHB interactions were confirmed in the monokaryotic S1-10 strain of L. edodes by showing increased expression of clp1, a downstream gene of the mating signal pathway and the occurrence of clamp connections after the treatment of PHBs. These results indicate that a single PHB can interact with a non-self RCB in a sublocus-specific manner for the activation of the mating pheromone signal pathways in L. edodes.

A Single Transcription Factor (PDD1) Determines Development and Yield of Winter Mushroom (Flammulina velutipes)

Most of the edible mushrooms cannot be cultivated or have low yield under industrial conditions, partially due to the lack of knowledge on how basidioma (fruiting body) development is regulated. From winter mushroom (Flammulina velutipes), one of the most popular industrially cultivated mushrooms, a transcription factor, PDD1, with a high-mobility group (HMG)-box domain was identified based on its increased transcription during basidioma development. pdd1 knockdown by RNA interference affected vegetative growth and dramatically impaired basidioma development. A strain with an 89.9% reduction in the level of pdd1 transcription failed to produce primordia, while overexpression of pdd1 promoted basidioma development. When the transcriptional level of pdd1 was increased to 5 times the base level, the mushroom cultivation time was shortened by 9.8% and the yield was increased by at least 33%. RNA sequencing (RNA-seq) analysis revealed that pdd1 knockdown downregulated 331 genes and upregulated 463 genes. PDD1 positively regulated several genes related to fruiting, including 6 pheromone receptor-encoding genes, 3 jacalin-related lectin-encoding genes, FVFD16, and 2 FVFD16 homolog-encoding genes. PDD1 is a novel transcription factor with regulatory function in basidioma development found in industrially cultivated mushrooms. Since its orthologs are widely present in fungal species of the Basidiomycota phylum, PDD1 might have important application prospects in mushroom breeding.IMPORTANCE Mushrooms are sources of food and medicine and provide abundant nutrients and bioactive compounds. However, most of the edible mushrooms cannot be cultivated commercially due to the limited understanding of basidioma development. From winter mushroom (Flammulina velutipes; also known as Enokitake), one of the most commonly cultivated mushrooms, we identified a novel transcription factor, PDD1, positively regulating basidioma development. PDD1 increases expression during basidioma development. Artificially increasing its expression promoted basidioma formation and dramatically increased mushroom yield, while reducing its expression dramatically impaired its development. In its PDD1 overexpression mutants, mushroom number, height, yield, and biological efficiency were significantly increased. PDD1 regulates the expression of some genes that are important in or related to basidioma development. PDD1 is the first identified transcription factor with defined functions in mushroom development among commercially cultivated mushroom species, and it might be useful in mushroom breeding.

Comparative Transcriptomics of Flammulina filiformis Suggests a High CO2 Concentration Inhibits Early Pileus Expansion by Decreasing Cell Division Control Pathways

Carbon dioxide is commonly used as one of the significant environmental factors to control pileus expansion during mushroom cultivation. However, the pileus expansion mechanism related to CO₂ is still unknown. In this study, the young fruiting bodies of a popular commercial mushroom Flammulina filiformis were cultivated under different CO₂ concentrations. In comparison to the low CO₂ concentration (0.05%), the pileus expansion rates were significantly lower under a high CO₂ concentration (5%). Transcriptome data showed that the up-regulated genes enriched in high CO₂ concentration treatments mainly associated with metabolism processes indicated that the cell metabolism processes were active under high CO₂ conditions. However, the gene ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with cell division processes contained down-regulated genes at both 12 h and 36 h under a high concentration of CO₂. Transcriptome and qRT-PCR analyses demonstrated that a high CO₂ concentration had an adverse effect on gene expression of the ubiquitin–proteasome system and cell cycle–yeast pathway, which may decrease the cell division ability and exhibit an inhibitory effect on early pileus expansion. Our research reveals the molecular mechanism of inhibition effects on early pileus expansion by elevated CO₂, which could provide a theoretical basis for a CO₂ management strategy in mushroom cultivation.

A second HD mating type sublocus of Flammulina velutipes is at least di-allelic and active: new primers for identification of HD-a and HD-b subloci

Background

Sexual development in Flammulina velutipes is controlled by two different mating type loci (HD and PR). The HD locus contains homeodomain (Hd) genes on two separate HD subloci: HD-a and HD-b. While the functionality of the HD-b sublocus has been largely confirmed, the status and content of the HD-a sublocus has remained unclear.

Methods

To examine the function of the HD-a sublocus, genome sequences of a series of F. velutipes strains were analyzed and tested through series of amplification by specific primer sets. Furthermore, activity of di-allelic HD-a locus was confirmed by crossing strains with different combinations of HD-a and HD-b subloci.

Results

Sublocus HD-b contained a large variety of fixed Hd1/Hd2 gene pairs, while the HD-a sublocus either contained a conserved Hd2 gene or, a newly discovered Hd1 gene that was also conserved. Identification of whole HD loci, that is, the contents of HD-a and HD-b subloci in a strain, revealed that strains with similar HD-b subloci could still form normal dikaryons if the two genes at the HD-a sublocus differed. At least di-allelic HD-a sublocus, is thus indicated to be actively involved in mating type compatibility.

Conclusions

HD-a sublocus is active and di-allelic. Using the new information on the HD subloci, primers sets were developed that specifically amplify HD-a or HD-b subloci in the majority of F. velutipes strains. In this way, unknown HD mating types of F. velutipes can now be quickly identified, and HD mating type compatibility conferred by HD-a or HD-b can be confirmed by PCR.

Diversity of A mating type in Lentinula edodes and mating type preference in the cultivated strains

Diversity of A mating type in Lentinula edodes has been assessed by analysis of A mating loci in 127 strains collected from East Asia. It was discovered that hypervariable sequence region with an approximate length of 1 kb in the A mating locus, spanning 5' region of HD2-intergenic region-5' region of HD1, could represent individual A mating type as evidenced by comprehensive mating analysis. The sequence analysis revealed 27 A mating type alleles from 96 cultivated strains and 48 alleles from 31 wild strains. Twelve of them commonly appeared, leaving 63 unique A mating type alleles. It was also revealed that only A few A mating type alleles such as A1, A4, A5, and A7 were prevalent in the cultivated strains, accounting for 62.5% of all A mating types. This implies preferred selection of certain A mating types in the process of strain development and suggests potential role of A mating genes in the expression of genes governing mushroom quality. Dominant expression of an A mating gene HD1 was observed from A1 mating locus, the most prevalent A allele, in A1-containing dikaryons. However, connections between HD1 expression and A1 preference in the cultivated strains remain to be verified. The A mating type was highly diverse in the wild strains. Thirty-six unique A alleles were discovered from relatively small and confined area of mountainous region in Korean peninsula. The number will further increase because no A allele has been recurrently observed in the wild strains and thus newly discovered strain will have good chances to contain new A allele. The high diversity in small area also suggests that the A mating locus has evolved rapidly and thus its diversity will further increase.

A genetic linkage map of Pleurotus tuoliensis integrated with physical mapping of the de novo sequenced genome and the mating type loci

Background

Pleurotus tuoliensis (Bailinggu) is a commercially cultivated mushroom species with an increasing popularity in China and other Asian countries. Commercial profits are now low, mainly due to a low yield, long cultivation period and sensitivity to diseases. Breeding efforts are thus required to improve agronomical important traits. Developing saturated genetic linkage and physical maps is a start for applying genetic and molecular approaches to accelerate the precise breeding programs.

Results

Here we present a genetic linkage map for P. tuoliensis constructed by using 115 haploid monokaryons derived from a hybrid strain H6. One thousand one hundred and eighty-two SNP markers developed by 2b–RAD (type IIB restriction-site associated DNA) approach were mapped to 12 linkage groups. The map covers 1073 cM with an average marker spacing of 1.0 cM. The genome of P. tuoliensis was de novo sequenced as 40.8 Mb and consisted of 500 scaffolds (>500 bp), which showed a high level of colinearity to the genome of P. eryngii var. eryngii. A total of 97.4% SNP markers (1151) were physically localized on 78 scaffolds, and the physical length of these anchored scaffolds were 33.9 Mb representing 83.1% of the whole genome. Mating type loci A and B were mapped on separate linkage groups and identified physically on the assembled genomes. Five putative pheromone receptors and two putative pheromone precursors were identified for the mating type B locus.

Conclusions

This study reported a first genetic linkage map integrated with physical mapping of the de novo sequenced genome and the mating type loci of an important cultivated mushroom in China, P. tuoliensis. The de novo sequenced and annotated genome, assembled using a 2b–RAD generated linkage map, provides a basis for marker-assisted breeding of this economic important mushroom species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4421-z) contains supplementary material, which is available to authorized users.