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.

Genomic determination of breeding systems and trans-specific evolution of HD MAT genes in suilloid fungi

Studying the signatures of evolution can help to understand genetic processes. Here, we demonstrate how the existence of balancing selection can be used to identify the breeding systems of fungi from genomic data. The breeding systems of fungi are controlled by self-incompatibility loci that determine mating types between potential mating partners, resulting in strong balancing selection at the loci. Within the fungal phylum Basidiomycota, two such self-incompatibility loci, namely HD MAT locus and P/R MAT locus, control mating types of gametes. Loss of function at one or both MAT loci results in different breeding systems and relaxes the MAT locus from balancing selection. By investigating the signatures of balancing selection at MAT loci, one can infer a species' breeding system without culture-based studies. Nevertheless, the extreme sequence divergence among MAT alleles imposes challenges for retrieving full variants from both alleles when using the conventional read-mapping method. Therefore, we employed a combination of read-mapping and local de novo assembly to construct haplotypes of HD MAT alleles from genomes in suilloid fungi (genera Suillus and Rhizopogon). Genealogy and pairwise divergence of HD MAT alleles showed that the origins of mating types predate the split between these two closely related genera. High sequence divergence, trans-specific polymorphism, and the deeply diverging genealogy confirm the long-term functionality and multiallelic status of HD MAT locus in suilloid fungi. This work highlights a genomics approach to studying breeding systems regardless of the culturability of organisms based on the interplay between evolution and genetics.

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

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

Homokaryotic High-Quality Genome Assembly of Medicinal Fungi Wolfiporia hoelen Reveals Auto-Regulation and High-Temperature Adaption of Probable Two-Speed Genome

Sclerotia of Wolfiporia hoelen are one of the most important traditional Chinese medicines and are commonly used in China, Japan, Korea, and other Asian countries. In the present study, we presented the first high-quality homokaryotic genome of W. hoelen with 14 chromosomes which was evaluated with assembly index, telomere position detection, and whole-genome collinearity. A 64.44 Mb genome was assembled with a Contig N50 length of 3.76 Mb. The imbalanced distribution of transposons and chromosome characters revealed the probable two-speed genome of W. hoelen. High consistency between methylation and transposon conserved the genome stability. The expansion of the gene family about signal transduction and nutritional transport has intimate relationships with sclerotial formation. Up-regulation of expression for distinctive decomposition enzymes, ROS clearance genes, biosynthesis of unsaturated fatty acids, and change of the cell wall components maintained high-speed growth of mycelia that may be the high-temperature adaption strategy of W. hoelen. Further, the analysis of mating-control genes demonstrated that HD3 probably had no function on mating recognition, with the HD protein in a distant genetic with known species. Overall, the high-quality genome of W. hoelen provided crucial information for genome structure and stability, high-temperature adaption, and sexual and asexual process.

Genome sequencing of Inonotus obliquus reveals insights into candidate genes involved in secondary metabolite biosynthesis

Background

Inonotus obliquus is an important edible and medicinal mushroom that was shown to have many pharmacological activities in preclinical trials, including anti-inflammatory, antitumor, immunomodulatory, and antioxidant effects. However, the biosynthesis of these pharmacological components has rarely been reported. The lack of genomic information has hindered further molecular characterization of this mushroom.

Results

In this study, we report the genome of I. obliquus using a combined high-throughput Illumina NovaSeq with Oxford Nanopore PromethION sequencing platform. The de novo assembled 38.18 Mb I. obliquus genome was determined to harbor 12,525 predicted protein-coding genes, with 81.83% of them having detectable sequence similarities to others available in public databases. Phylogenetic analysis revealed the close evolutionary relationship of I. obliquus with Fomitiporia mediterranea and Sanghuangporus baumii in the Hymenochaetales clade. According to the distribution of reproduction-related genes, we predict that this mushroom possesses a tetrapolar heterothallic reproductive system. The I. obliquus genome was found to encode a repertoire of enzymes involved in carbohydrate metabolism, along with 135 cytochrome P450 proteins. The genome annotation revealed genes encoding key enzymes responsible for secondary metabolite biosynthesis, such as polysaccharides, polyketides, and terpenoids. Among them, we found four polyketide synthases and 20 sesquiterpenoid synthases belonging to four more types of cyclization mechanism, as well as 13 putative biosynthesis gene clusters involved in terpenoid synthesis in I. obliquus.

Conclusions

To the best of our knowledge, this is the first reported genome of I. obliquus; we discussed its genome characteristics and functional annotations in detail and predicted secondary metabolic biosynthesis-related genes, which provides genomic information for future studies on its associated molecular mechanism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08511-x.

Fruit body formation and intra-species DNA polymorphism in Japanese Wolfiporia cocos strains

Poria, the dried sclerotium of Wolfiporia cocos, is a medicinal mushroom that is widely used in traditional Japanese medicine. The fruit body of W. cocos is rarely found in the natural environment in Japan, therefore an optimized technique for fruit body formation is essential for producing new strains through crossbreeding and for biological research. Here, we developed a cultivation technique for fruit body formation of W. cocos using three strains collected from different areas of Japan. When mycelia were cultured on sawdust-based medium after liquid medium culture, all strains successfully formed fruit bodies as a brown honeycomb-like structure. Furthermore, we analyzed single nucleotide polymorphisms of the three strains using the STE3-like pheromone receptor protein gene, STE3.2, and found a genetic marker for discriminating one strain from the others. The results are expected to promote extensive studies on crossbreeding and domestic production of W. cocos.

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.

Sexual Crossing, Chromosome-Level Genome Sequences, and Comparative Genomic Analyses for the Medicinal Mushroom Taiwanofungus Camphoratus (Syn. Antrodia Cinnamomea, Antrodia Camphorata)

Taiwanofungus camphoratus mushrooms are a complementary and alternative medicine for hangovers, cancer, hypertension, obesity, diabetes, and inflammation. Though Taiwanofungus camphoratus has attracted considerable biotechnological and pharmacological attention, neither classical genetic nor genomic approaches have been properly established for it. We isolated four sexually competent monokaryons from two T. camphoratus dikaryons used for the commercial cultivation of orange-red (HC1) and milky-white (SN1) mushrooms, respectively. We also sequenced, annotated, and comparatively analyzed high-quality and chromosome-level genome sequences of these four monokaryons. These genomic resources represent a valuable basis for understanding the biology, evolution, and secondary metabolite biosynthesis of this economically important mushrooms. We demonstrate that T. camphoratus has a tetrapolar mating system and that HC1 and SN1 represent two intraspecies isolates displaying karyotypic variation. Compared with several edible mushroom model organisms, T. camphoratus underwent a significant contraction in the gene family and individual gene numbers, most notably for plant, fungal, and bacterial cell-wall-degrading enzymes, explaining why T. camphoratus mushrooms are rare in natural environments, are difficult and time-consuming to artificially cultivate, and are susceptible to fungal and bacterial infections. Our results lay the foundation for an in-depth T. camphoratus study, including precise genetic manipulation, improvements to mushroom fruiting, and synthetic biology applications for producing natural medicinal products. IMPORTANCETaiwanofungus camphoratus (Tc) is a basidiomycete fungus that causes brown heart rot of the aromatic tree Cinnamomum kanehirae. The Tc fruiting bodies have been used to treat hangovers, abdominal pain, diarrhea, hypertension, and other diseases first by aboriginal Taiwanese and later by people in many countries. To establish classical genetic and genomic approaches for this economically important medicinal mushroom, we first isolated and characterized four sexually competent monokaryons from two dikaryons wildly used for commercial production of Tc mushrooms. We applied PacBio single molecule, real-time sequencing technology to determine the near-completed genome sequences of four monokaryons. These telomere-to-telomere and gapless haploid genome sequences reveal all genomic variants needed to be studied and discovered, including centromeres, telomeres, retrotransposons, mating type loci, biosynthetic, and metabolic gene clusters. Substantial interspecies diversities are also discovered between Tc and several other mushroom model organisms, including Agrocybe aegerita, Coprinopsis cinerea, and Schizophyllum commune, and Ganoderma lucidum.

Genomic and Transcriptomic Insight of Giant Sclerotium Formation of Wood-Decay Fungi

Many fungi form persistent and dormant sclerotia with compact hardened mycelia during unfavorable circumstances. While most of these sclerotia are small in size, Wolfiporia cocos, a wood-decay fungus, grows into giant sclerotia, which are mainly composed of polysaccharides of linear (1→3)-β-D-glucans. To explore the underlying mechanism of converting sophisticated wood polysaccharides for biosynthesis of highly homogenized glucans in W. cocos, we sequenced and assembled the genome of a cultivated W. cocos strain (WCLT) in China. The 62-Mb haploid genome contains 44.2% repeat sequences, of which, 48.0% are transposable elements (TEs). Contrary to the genome of W. cocos from North America, WCLT has independently undergone a partial genome duplication (PGD) event. The large-scale TE insertion and PGD occurrence overlapped with an archeological Pleistocene stage of low oxygen and high temperature, and these stresses might have induced the differences in sclerotium due to geographical distribution. The wood decomposition enzymes, as well as sclerotium-regulator kinases, aquaporins, and highly expanded gene families such as NAD-related families, together with actively expressed 1,3-β-glucan synthase for sclerotium polysaccharides, all have contributed to the sclerotium formation and expansion. This study shall inspire further exploration on how fungi convert wood into simple glucans in the sclerotium of W. cocos.

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

In heterothallic basidiomycete fungi, sexual compatibility is restricted by mating types, typically controlled by two loci: PR, encoding pheromone precursors and pheromone receptors, and HD, encoding two types of homeodomain transcription factors. We analysed the single mating-type locus of the commercial button mushroom variety, Agaricus bisporus var. bisporus, and of the related variety burnettii. We identified the location of the mating-type locus using genetic map and genome information, corresponding to the HD locus, the PR locus having lost its mating-type role. We found the mip1 and β-fg genes flanking the HD genes as in several Agaricomycetes, two copies of the β-fg gene, an additional HD2 copy in the reference genome of A. bisporus var. bisporus and an additional HD1 copy in the reference genome of A. bisporus var. burnettii. We detected a 140 kb-long inversion between mating types in an A. bisporus var. burnettii heterokaryon, trapping the HD genes, the mip1 gene and fragments of additional genes. The two varieties had islands of transposable elements at the mating-type locus, spanning 35 kb in the A. bisporus var. burnettii reference genome. Linkage analyses showed a region with low recombination in the mating-type locus region in the A. bisporus var. burnettii variety. We found high differentiation between β-fg alleles in both varieties, indicating an ancient event of recombination suppression, followed more recently by a suppression of recombination at the mip1 gene through the inversion in A. bisporus var. burnettii and a suppression of recombination across whole chromosomes in A. bisporus var. bisporus, constituting stepwise recombination suppression as in many other mating-type chromosomes and sex chromosomes.

Macrosynteny analysis between Lentinula edodes and Lentinula novae-zelandiae reveals signals of domestication in Lentinula edodes

The basidiomycete fungus Lentinula novae-zelandiae is endemic to New Zealand and is a sister taxon to Lentinula edodes, the second most cultivated mushroom in the world. To explore the biology of this organism, a high-quality chromosome level reference genome of L. novae-zelandiae was produced. Macrosyntenic comparisons between the genome assembly of L. novae-zelandiae, L. edodes and a set of three genome assemblies of diverse species from the Agaricomycota reveal a high degree of macrosyntenic restructuring within L. edodes consistent with signal of domestication. These results show L. edodes has undergone significant genomic change during the course of its evolutionary history, likely a result of its cultivation and domestication over the last 1000 years.

Distinguishing Homokaryons and Heterokaryons in Medicinal Polypore Mushroom Wolfiporia cocos (Agaricomycetes) Based on Cultural and Genetic Characteristics

The sclerotia of Wolfiporia cocos are a kind of traditional medicine and food with excellent benefits and are widely used in China, Japan, and other Asian countries. The mating system of fungi is not only of practical importance for breeding but also has profound effects on genetic variability and molecular evolution. However, the lack of clamp connections in W. cocos increases the difficulty of research on mating systems. In this study, homokaryons and heterokaryons were distinguished by comparing the characteristics of culture, fruiting tests, and molecular markers, which was further demonstrated by k-mer analysis based on Illumina sequencing. Uninucleate, binucleate, and nuclei-free condition basidiospores of W. cocos were observed, and binucleate basidiospores were the most predominant. Brown-type colonies, slow growth rates in both PDA medium and sawdust substrate, and neutral pH after the growth of mycelia and unfruiting were found to be the morphological and growth characteristics of homokaryotic strains. Primers SSR37 and 38 were screened to identify homokaryons. K-mer analysis based on Illumina sequencing exhibited different heterozygous ratios for homokaryons and heterokaryons. The results revealed that pseudo-homothallism was the predominant mode of reproduction in the Chinese population of W. cocos, and heterothallism also existed in all probability. This study will be helpful for the cross-breeding of this precious medicinal mushroom and for understanding its evolution and population structure.

The Wolfiporia cocos Genome and Transcriptome Shed Light on the Formation of Its Edible and Medicinal Sclerotium

Wolfiporia cocos (F. A. Wolf) has been praised as a food delicacy and medicine for centuries in China. Here, we present the genome and transcriptome of the Chinese strain CGMCC5.78 of W. cocos. High-confidence functional prediction was made for 9277 genes among the 10,908 total predicted gene models in the W. cocos genome. Up to 2838 differentially expressed genes (DEGs) were identified to be related to sclerotial development by comparing the transcriptomes of mycelial and sclerotial tissues. These DEGs are involved in mating processes, differentiation of fruiting body tissues, and metabolic pathways. A number of genes encoding enzymes and regulatory factors related to polysaccharide and triterpenoid production were strikingly regulated. A potential triterpenoid gene cluster including the signature lanosterol synthase (LSS) gene and its modified components were annotated. In addition, five nonribosomal peptide synthase (NRPS)-like gene clusters, eight polyketide synthase (PKS) gene clusters, and 15 terpene gene clusters were discovered in the genome. The differential expression of the velevt family proteins, transcription factors, carbohydrate-active enzymes, and signaling components indicated their essential roles in the regulation of fungal development and secondary metabolism in W. cocos. These genomic and transcriptomic resources will be valuable for further investigations of the molecular mechanisms controlling sclerotial formation and for its improved medicinal applications.

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.

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.

Nine draft genome sequences of Claviceps purpurea s.lat., including C. arundinis, C. humidiphila, and C. cf. spartinae, pseudomolecules for the pitch canker pathogen Fusarium circinatum, draft genome of Davidsoniella eucalypti, Grosmannia galeiformis, Quambalaria eucalypti, and Teratosphaeria destructans

This genome announcement includes draft genomes from Claviceps purpurea s.lat., including C. arundinis, C. humidiphila and C. cf. spartinae. The draft genomes of Davidsoniella eucalypti, Quambalaria eucalypti and Teratosphaeria destructans, all three important eucalyptus pathogens, are presented. The insect associate Grosmannia galeiformis is also described. The pine pathogen genome of Fusarium circinatum has been assembled into pseudomolecules, based on additional sequence data and by harnessing the known synteny within the Fusarium fujikuroi species complex. This new assembly of the F. circinatum genome provides 12 pseudomolecules that correspond to the haploid chromosome number of F. circinatum. These are comparable to other chromosomal assemblies within the FFSC and will enable more robust genomic comparisons within this species complex.

Population structure of an orchid mycorrhizal fungus with genus-wide specificity

Fundamental life history processes of mycorrhizal fungi with inconspicuous fruiting bodies can be difficult to elucidate. In this study we investigated the species identities and life history of the orchid mycorrhizal Tulasnella fungi, which associate with the south eastern Australia orchid genus Chiloglottis. Tulasnella prima was the primary partner and was found to be associated with all 17 Chiloglottis species across a range of >1000 km, and to occur in the two edaphic conditions investigated (soil and sphagnum hammocks). Another Tulasnella species (T. sphagneti) appears to be restricted to moist conditions of alpine sphagnum hammocks. The population genetic structure of the widespread species T. prima, was investigated at 10 simple sequence repeat (SSR) markers and at four cross-amplified SSR loci for T. sphagneti. For both taxa, no sharing of multilocus genotypes was found between sites, but clones were found within sites. Evidence for inbreeding within T. prima was found at 3 of 5 sites. Significant genetic differentiation was found within and between taxa. Significant local positive spatial genetic autocorrelation was detected among non-clonal isolates at the scale of two metres. Overall, the population genetic patterns indicated that in Tulasnella mating occurs by inbreeding and dispersal is typically restricted to short-distances.

A revised family-level classification of the Polyporales (Basidiomycota)

Polyporales is strongly supported as a clade of Agaricomycetes, but the lack of a consensus higher-level classification within the group is a barrier to further taxonomic revision. We amplified nrLSU, nrITS, and rpb1 genes across the Polyporales, with a special focus on the latter. We combined the new sequences with molecular data generated during the PolyPEET project and performed Maximum Likelihood and Bayesian phylogenetic analyses. Analyses of our final 3-gene dataset (292 Polyporales taxa) provide a phylogenetic overview of the order that we translate here into a formal family-level classification. Eighteen clades are assigned a family name, including three families described as new (Cerrenaceae fam. nov., Gelatoporiaceae fam. nov., Panaceae fam. nov.) and fifteen others (Dacryobolaceae, Fomitopsidaceae, Grifolaceae, Hyphodermataceae, Incrustoporiaceae, Irpicaceae, Ischnodermataceae, Laetiporaceae, Meripilaceae, Meruliaceae, Phanerochaetaceae, Podoscyphaceae, Polyporaceae, Sparassidaceae, Steccherinaceae). Three clades are given informal names (/hypochnicium,/climacocystis and/fibroporia + amyloporia). Four taxa (Candelabrochete africana, Mycoleptodonoides vassiljevae, Auriporia aurea, and Tyromyces merulinus) cannot be assigned to a family within the Polyporales. The classification proposed here provides a framework for further taxonomic revision and will facilitate communication among applied and basic scientists. A survey of morphological, anatomical, physiological, and genetic traits confirms the plasticity of characters previously emphasized in taxonomy of Polyporales.

Comparative Genomics of the Ectomycorrhizal Sister Species Rhizopogon vinicolor and Rhizopogon vesiculosus (Basidiomycota: Boletales) Reveals a Divergence of the Mating Type B Locus

Divergence of breeding system plays an important role in fungal speciation. Ectomycorrhizal fungi, however, pose a challenge for the study of reproductive biology because most cannot be mated under laboratory conditions. To overcome this barrier, we sequenced the draft genomes of the ectomycorrhizal sister species Rhizopogon vinicolor Smith and Zeller and R. vesiculosus Smith and Zeller (Basidiomycota, Boletales)-the first genomes available for Basidiomycota truffles-and characterized gene content and organization surrounding their mating type loci. Both species possess a pair of homeodomain transcription factor homologs at the mating type A-locus as well as pheromone receptor and pheromone precursor homologs at the mating type B-locus. Comparison of Rhizopogon genomes with genomes from Boletales, Agaricales, and Polyporales revealed synteny of the A-locus region within Boletales, but several genomic rearrangements across orders. Our findings suggest correlation between gene content at the B-locus region and breeding system in Boletales with tetrapolar species possessing more diverse gene content than bipolar species. Rhizopogon vinicolor possesses a greater number of B-locus pheromone receptor and precursor genes than R. vesiculosus, as well as a pair of isoprenyl cysteine methyltransferase genes flanking the B-locus compared to a single copy in R. vesiculosus Examination of dikaryotic single nucleotide polymorphisms within genomes revealed greater heterozygosity in R. vinicolor, consistent with increased rates of outcrossing. Both species possess the components of a heterothallic breeding system with R. vinicolor possessing a B-locus region structure consistent with tetrapolar Boletales and R. vesiculosus possessing a B-locus region structure intermediate between bipolar and tetrapolar Boletales.

Fungal Sex: The Basidiomycota

Fungi of the Basidiomycota, representing major pathogen lineages and mushroom-forming species, exhibit diverse means to achieve sexual reproduction, with particularly varied mechanisms to determine compatibilities of haploid mating partners. For species that require mating between distinct genotypes, discrimination is usually based on both the reciprocal exchange of diffusible mating pheromones, rather than sexes, and the interactions of homeodomain protein signals after cell fusion. Both compatibility factors must be heterozygous in the product of mating, and genetic linkage relationships of the mating pheromone/receptor and homeodomain genes largely determine the complex patterns of mating-type variation. Independent segregation of the two compatibility factors can create four haploid mating genotypes from meiosis, referred to as tetrapolarity. This condition is thought to be ancestral to the basidiomycetes. Alternatively, cosegregation by linkage of the two mating factors, or in some cases the absence of the pheromone-based discrimination, yields only two mating types from meiosis, referred to as bipolarity. Several species are now known to have large and highly rearranged chromosomal regions linked to mating-type genes. At the population level, polymorphism of the mating-type genes is an exceptional aspect of some basidiomycete fungi, where selection under outcrossing for rare, intercompatible allelic variants is thought to be responsible for numbers of mating types that may reach several thousand. Advances in genome sequencing and assembly are yielding new insights by comparative approaches among and within basidiomycete species, with the promise to resolve the evolutionary origins and dynamics of mating compatibility genetics in this major eukaryotic lineage.

The cacao pathogen Moniliophthora roreri (Marasmiaceae) possesses biallelic A and B mating loci but reproduces clonally

The cacao pathogen Moniliophthora roreri belongs to the mushroom-forming family Marasmiaceae, but it has never been observed to produce a fruiting body, which calls to question its capacity for sexual reproduction. In this study, we identified potential A (HD1 and HD2) and B (pheromone precursors and pheromone receptors) mating genes in M. roreri. A PCR-based method was subsequently devised to determine the mating type for a set of 47 isolates from across the geographic range of the fungus. We developed and generated an 11-marker microsatellite set and conducted association and linkage disequilibrium (standardized index of association, IA(s)) analyses. We also performed an ancestral reconstruction analysis to show that the ancestor of M. roreri is predicted to be heterothallic and tetrapolar, which together with sliding window analyses support that the A and B mating loci are likely unlinked and follow a tetrapolar organization within the genome. The A locus is composed of a pair of HD1 and HD2 genes, whereas the B locus consists of a paired pheromone precursor, Mr_Ph4, and receptor, STE3_Mr4. Two A and B alleles but only two mating types were identified. Association analyses divided isolates into two well-defined genetically distinct groups that correlate with their mating type; IA(s) values show high linkage disequilibrium as is expected in clonal reproduction. Interestingly, both mating types were found in South American isolates but only one mating type was found in Central American isolates, supporting a prior hypothesis of clonal dissemination throughout Central America after a single or very few introductions of the fungus from South America.

Development of crossbreeding high-yield-potential strains for commercial cultivation in the medicinal mushroom Wolfiporia cocos (Higher Basidiomycetes)

Wolfiporia cocos is a well-known medicinal mushroom, and its dried sclerotia has been widely used as a traditional medicine in China, Japan, and other Asian countries for centuries. However, long-term asexual reproduction of the breeding system in W. cocos results in a current universal degeneration of cultivated strains. To develop a W. cocos breeding program that will benefit commercial cultivation, we previously developed an optimum method for indoor induction of W. cocos fruiting bodies and clarified the nature of preponderant binuclear sexual basidiospores. In this paper, we first show that the majority of W. cocos single-spore isolates cannot form sclerotium in field cultivation. We then investigated the possibility of breeding new strains by crossbreeding. Three types of mating reactions were observed in both intra-strain pairings and inter-strain pairings, and a total of fifty-five hybrids were selected by antagonistic testing and allele-specific polymerase chain reaction (PCR). Field cultivation of hybrids demonstrated that some hybrids can form sclerotium via two cultivated methods. Two new high-yield strains were identified. This report will stimulate new thinking on W. cocos and promote further extensive studies on crossbreeding in W. cocos, a new topic related to the development of more efficient protocols for the discrimination of hybrids in W. cocos.

Cloning of the A Mating-Type Locus from Lepista nuda and Characterization of Its Genetic Structure

Lepista nuda (Bull. ex Fr.) Cooke (Tricholomataceae) is an edible fungus with both economic and medical value. Identification of its mating-type loci is important for promoting breeding programs in L. nuda. The A mating-type locus of L. nuda and its flanking region were cloned and characterized in the present study. It contained two homeodomain transcription factor genes (called lna1 and lna2). Lna1 and Lna2 protein harbored conserved motif of homeodomain transcription factor protein. The novel finding of this study was that the gene order around the A locus was mip, lna2, lna1, and β-fg in L. nuda, which was differed from other edible fungi. In addition, lna1 and lna2 showed divergent, inward transcriptional direction. The phylogenetic tree of HD proteins showed that L. nuda Lna1 and Lna2 were phylogenetically related with Laccaria bicolor. Our results revealed that the A mating-type locus had been subjected to gene rearrangements relative to all other basidiomycetes.

Evolution of Mating Systems in Basidiomycetes and the Genetic Architecture Underlying Mating-Type Determination in the Yeast Leucosporidium scottii

In most fungi, sexual reproduction is bipolar; that is, two alternate sets of genes at a single mating-type (MAT) locus determine two mating types. However, in the Basidiomycota, a unique (tetrapolar) reproductive system emerged in which sexual identity is governed by two unlinked MAT loci, each of which controls independent mechanisms of self/nonself recognition. Tetrapolar-to-bipolar transitions have occurred on multiple occasions in the Basidiomycota, resulting, for example, from linkage of the two MAT loci into a single inheritable unit. Nevertheless, owing to the scarcity of molecular data regarding tetrapolar systems in the earliest-branching lineage of the Basidiomycota (subphylum Pucciniomycotina), it is presently unclear if the last common ancestor was tetrapolar or bipolar. Here, we address this question, by investigating the mating system of the Pucciniomycotina yeast Leucosporidium scottii. Using whole-genome sequencing and chromoblot analysis, we discovered that sexual reproduction is governed by two physically unlinked gene clusters: a multiallelic homeodomain (HD) locus and a pheromone/receptor (P/R) locus that is biallelic, thereby dismissing the existence of a third P/R allele as proposed earlier. Allele distribution of both MAT genes in natural populations showed that the two loci were in strong linkage disequilibrium, but independent assortment of MAT alleles was observed in the meiotic progeny of a test cross. The sexual cycle produces fertile progeny with similar proportions of the four mating types, but approximately 2/3 of the progeny was found to be nonhaploid. Our study adds to others in reinforcing tetrapolarity as the ancestral state of all basidiomycetes.

Basidioascus undulatus: genome, origins, and sexuality

Basidioascus undulatus is a soil basidiomycete belonging to the order Geminibasidiales. The taxonomic status of the order was unclear as originally it was only tentatively classified in the class Wallemiomycetes. The fungi in Geminibasidiales have an ambiguously defined sexual cycle. In this study, we sequenced the genome of B. undulatus to gain insights into its sexuality and evolutionary origins. The assembled genome draft was approximately 32 Mb in size, had a median nucleotide coverage of 24X, and contained 6123 predicted genes. Previous morphological descriptions of B. undulatus relied on interpretation of putative sexual structures. In this study, nuclear staining and confocal microscopy showed meiosis occurring in basidia and genome analysis confirmed the existence of genes involved in meiosis and mating. Using 35 protein-coding genes extracted from genomic information, phylogenomic and molecular dating analyses confirmed that B. undulatus indeed belongs to a lineage distantly related to Wallemia while retaining a basal position in Agaricomycotina. These results, combined with differences in septal pore morphology, led us to move the order Geminibasidiales out of the Wallemiomycetes and into the new class Geminibasidiomycetes cl. nov. Finally, the concept of Agaricomycotina is emended to include both Wallemiomycetes and Geminibasidiomycetes.

Genome analysis of medicinal Ganoderma spp. with plant-pathogenic and saprotrophic life-styles

Ganoderma is a fungal genus belonging to the Ganodermataceae family and Polyporales order. Plant-pathogenic species in this genus can cause severe diseases (stem, butt, and root rot) in economically important trees and perennial crops, especially in tropical countries. Ganoderma species are white rot fungi and have ecological importance in the breakdown of woody plants for nutrient mobilization. They possess effective machineries of lignocellulose-decomposing enzymes useful for bioenergy production and bioremediation. In addition, the genus contains many important species that produce pharmacologically active compounds used in health food and medicine. With the rapid adoption of next-generation DNA sequencing technologies, whole genome sequencing and systematic transcriptome analyses become affordable approaches to identify an organism's genes. In the last few years, numerous projects have been initiated to identify the genetic contents of several Ganoderma species, particularly in different strains of Ganoderma lucidum. In November 2013, eleven whole genome sequencing projects for Ganoderma species were registered in international databases, three of which were already completed with genomes being assembled to high quality. In addition to the nuclear genome, two mitochondrial genomes for Ganoderma species have also been reported. Complementing genome analysis, four transcriptome studies on various developmental stages of Ganoderma species have been performed. Information obtained from these studies has laid the foundation for the identification of genes involved in biological pathways that are critical for understanding the biology of Ganoderma, such as the mechanism of pathogenesis, the biosynthesis of active components, life cycle and cellular development, etc. With abundant genetic information becoming available, a few centralized resources have been established to disseminate the knowledge and integrate relevant data to support comparative genomic analyses of Ganoderma species. The current review carries out a detailed comparison of the nuclear genomes, mitochondrial genomes and transcriptomes from several Ganoderma species. Genes involved in biosynthetic pathways such as CYP450 genes and in cellular development such as matA and matB genes are characterized and compared in detail, as examples to demonstrate the usefulness of comparative genomic analyses for the identification of critical genes. Resources needed for future data integration and exploitation are also discussed.

Genomic and transcriptomic analyses of the medicinal fungus Antrodia cinnamomea for its metabolite biosynthesis and sexual development

Antrodia cinnamomea, a polyporus mushroom of Taiwan, has long been used as a remedy for cancer, hypertension, and hangover, with an annual market of over $100 million (US) in Taiwan. We obtained a 32.15-Mb genome draft containing 9,254 genes. Genome ontology enrichment and pathway analyses shed light on sexual development and the biosynthesis of sesquiterpenoids, triterpenoids, ergostanes, antroquinonol, and antrocamphin. We identified genes differentially expressed between mycelium and fruiting body and 242 proteins in the mevalonate pathway, terpenoid pathways, cytochrome P450s, and polyketide synthases, which may contribute to the production of medicinal secondary metabolites. Genes of secondary metabolite biosynthetic pathways showed expression enrichment for tissue-specific compounds, including 14-α-demethylase (CYP51F1) in fruiting body for converting lanostane to ergostane triterpenoids, coenzymes Q (COQ) for antroquinonol biosynthesis in mycelium, and polyketide synthase for antrocamphin biosynthesis in fruiting body. Our data will be useful for developing a strategy to increase the production of useful metabolites.

The genome of the white-rot fungus Pycnoporus cinnabarinus: a basidiomycete model with a versatile arsenal for lignocellulosic biomass breakdown

BACKGROUND

Saprophytic filamentous fungi are ubiquitous micro-organisms that play an essential role in photosynthetic carbon recycling. The wood-decayer Pycnoporus cinnabarinus is a model fungus for the study of plant cell wall decomposition and is used for a number of applications in green and white biotechnology.

RESULTS

The 33.6 megabase genome of P. cinnabarinus was sequenced and assembled, and the 10,442 predicted genes were functionally annotated using a phylogenomic procedure. In-depth analyses were carried out for the numerous enzyme families involved in lignocellulosic biomass breakdown, for protein secretion and glycosylation pathways, and for mating type. The P. cinnabarinus genome sequence revealed a consistent repertoire of genes shared with wood-decaying basidiomycetes. P. cinnabarinus is thus fully equipped with the classical families involved in cellulose and hemicellulose degradation, whereas its pectinolytic repertoire appears relatively limited. In addition, P. cinnabarinus possesses a complete versatile enzymatic arsenal for lignin breakdown. We identified several genes encoding members of the three ligninolytic peroxidase types, namely lignin peroxidase, manganese peroxidase and versatile peroxidase. Comparative genome analyses were performed in fungi displaying different nutritional strategies (white-rot and brown-rot modes of decay). P. cinnabarinus presents a typical distribution of all the specific families found in the white-rot life style. Growth profiling of P. cinnabarinus was performed on 35 carbon sources including simple and complex substrates to study substrate utilization and preferences. P. cinnabarinus grew faster on crude plant substrates than on pure, mono- or polysaccharide substrates. Finally, proteomic analyses were conducted from liquid and solid-state fermentation to analyze the composition of the secretomes corresponding to growth on different substrates. The distribution of lignocellulolytic enzymes in the secretomes was strongly dependent on growth conditions, especially for lytic polysaccharide mono-oxygenases.

CONCLUSIONS

With its available genome sequence, P. cinnabarinus is now an outstanding model system for the study of the enzyme machinery involved in the degradation or transformation of lignocellulosic biomass.

Toward genome-enabled mycology

Genome-enabled mycology is a rapidly expanding field that is characterized by the pervasive use of genome-scale data and associated computational tools in all aspects of fungal biology. Genome-enabled mycology is integrative and often requires teams of researchers with diverse skills in organismal mycology, bioinformatics and molecular biology. This issue of Mycologia presents the first complete fungal genomes in the history of the journal, reflecting the ongoing transformation of mycology into a genome-enabled science. Here, we consider the prospects for genome-enabled mycology and the technical and social challenges that will need to be overcome to grow the database of complete fungal genomes and enable all fungal biologists to make use of the new data.