Whole genome de novo sequencing and genome annotation of the world popular cultivated edible mushroom, Lentinula edodes

Chromosome-level assembly of Dictyophora rubrovolvata genome using third-generation DNA sequencing and Hi-C analysis

Dictyophora rubrovolvata, a rare edible mushroom with both nutritional and medicinal values, was regarded as the "queen of the mushroom" for its attractive appearance. Dictyophora rubrovolvata has been widely cultivated in China in recent years, and many researchers were focusing on its nutrition, culture condition, and artificial cultivation. Due to a lack of genomic information, research on bioactive substances, cross breeding, lignocellulose degradation, and molecular biology is limited. In this study, we report a chromosome-level reference genome of D. rubrovolvata using the PacBio single-molecule real-time-sequencing technique and high-throughput chromosome conformation capture (Hi-C) technologies. A total of 1.83 Gb circular consensus sequencing reads representing ∼983.34 coverage of the D. rubrovolvata genome were generated. The final genome was assembled into 136 contigs with a total length of 32.89 Mb. The scaffold and contig N50 length were 2.71 and 2.48 Mb, respectively. After chromosome-level scaffolding, 11 chromosomes with a total length of 28.24 Mb were constructed. Genome annotation further revealed that 9.86% of the genome was composed of repetitive sequences, and a total of 508 noncoding RNA (rRNA: 329, tRNA: 150, ncRNA: 29) were annotated. In addition, 9,725 protein-coding genes were predicted, among which 8,830 (90.79%) genes were predicted using homology or RNA-seq. Benchmarking Universal Single-Copy Orthologs results further revealed that there were 80.34% complete single-copy fungal orthologs. In this study, a total of 360 genes were annotated as belonging to the carbohydrate-active enzymes family. Further analysis also predicted 425 cytochromes P450 genes, which can be classified into 41 families. This highly accurate, chromosome-level reference genome of D. rubrovolvata will provide essential genomic information for understanding the molecular mechanism in its fruiting body formation during morphological development and facilitate the exploitation of medicinal compounds produced by this mushroom.

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.

Proteomic characterization of Shiitake (Lentinula edodes) post-harvest fruit bodies grown on hardwood logs and isolation of an antibacterial serine protease inhibitor

Lentinula edodes (Shiitake) is one of the most heavily cultivated mushrooms in the world with proven antioxidant and antibacterial properties, among others. Evidence indicates that the choice of mushroom cultivation technique strongly influences the production of bioactive compounds, but to date the nature of many of these compounds has not been fully established. This work focuses on the proteomic characterization of L. edodes to highlight the main active processes two days after harvest and elucidates the proteins involved in the known antioxidant and antibacterial proprieties of Shiitake fruit bodies cultivated on oak logs. A label-free approach allowed us to identify a total of 2702 proteins which were mainly involved in carbohydrate and protein metabolism, cell growth and replication, indicating that several developmental processes remain active in fruit bodies post-harvest. Proteins with antioxidant activities were identified, indicating the contribution of proteins to the antioxidant properties of L. edodes extracts. Antibacterial assays also reveal the activity of a serine protease inhibitor that strongly accumulates in the post-harvest fruit body grown on oak logs. Overall, this study contributes to the understanding of the impact of the log cultivation method on the production of Shiitake mushrooms richest in high-value bioactive compounds.

Comparative structural analysis on the mitochondrial DNAs from various strains of Lentinula edodes

The evolution of mitochondria through variations in mitochondrial DNA (mtDNA) is one of the intriguing questions in eukaryotic cells. In order to assess the causes of the variations in mitochondria, the mtDNAs of the 21 strains of Lentinula edodes were assembled for this study, and analyzed together with four published mtDNA sequences. The mtDNAs were within the sizes of 117 kb ~ 122 kb. The gene number was observed consistent except for two mtDNAs, which carry a duplicated trnG1-trnG2 unit or a putative gene deletion. The size variation was largely attributed to the number of introns, repeated sequences, transposable elements (TEs), and plasmid-related sequences. Intron loss and gain were found from cox1, rnl, and rns of three mtDNAs. Loss of two introns in cox1 of KY217797.1 reduced its size by 2.7 kb, making it the smallest cox1 gene (8.4 kb) among the cox1s of the 25 mtDNAs, whereas gain of a Group II intron (2.65 kb) and loss of a Group I intron (1.7 kb) in cox1 of MF774813.1 resulted in the longest cox1 (12 kb). In rnl of L. edodes, we discovered four intron insertion consensus sequences which were unique to basidiomycetes but not ascomycetes. Differential incorporation of introns was the primary cause of the rnl size polymorphism. Homing endonucleases (HEGs) were suggestively involved in the mobilization of the introns because all of the introns have HEG genes of the LAGRIDADG or GIY-YIG families with the conserved HEG cleavage sites. TEs contributed to 11.04% of the mtDNA size in average, of which 7.08% was LTR-retrotransposon and 3.96% was DNA transposon, whereas the repeated sequences covered 4.6% of the mtDNA. The repeat numbers were variable in a strain-dependent manner. Both the TEs and repeated sequences were mostly found in the intronic and intergenic regions. Lastly, two major deletions were found in the plasmid-related sequence regions (pol2-pol3 and pol1-atp8) in the five mtDNAs. Particularly, the 6.8 kb-long deletion at pol2-pol3 region made MF774813.1 the shortest mtDNA of all. Our results demonstrate that mtDNA is a dynamic molecule that persistently evolves over a short period of time by insertion/deletion and repetition of DNA segments at the strain level.

UVB Irradiation-Induced Transcriptional Changes in Lignin- and Flavonoid Biosynthesis and Indole/Tryptophan-Auxin-Responsive Genes in Rice Seedlings

Global warming accelerates the destruction of the ozone layer, increasing the amount of UVB reaching the Earth’s surface, which in turn alters plant growth and development. The effects of UVB-induced alterations of plant secondary and cell wall metabolism were previously documented; however, there is little knowledge of its effects on rice seedlings during the developmental phase of leaves. In this study, we examined secondary metabolic responses to UVB stress using a transcriptomic approach, focusing on the biosynthetic pathways for lignin, flavonoid, and indole/tryptophan-auxin responses. As new leaves emerged, they were irradiated with UVB for 5 days (for 3 h/day⁻¹). The genes encoding the enzymes related to lignin (4CL, CAD, and POD) and flavonoid biosynthesis (CHS, CHI, and FLS) were highly expressed on day 1 (younger leaves) and day 5 (older leaves) after UVB irradiation. The expression of the genes encoding the enzymes related to tryptophan biosynthesis (AS, PRT, PRAI, IGPS, and TS) increased on day 3 of UVB irradiation, and the level of tryptophan increased and showed the same temporal pattern of occurrence as the expression of the cognate gene. Interestingly, the genes encoding BBX4 and BBX11, negative regulators of UVB signaling, and SAUR27 and SAUR55, auxin response enzymes, were downregulated on day 3 of UVB irradiation. When these results are taken together, they suggest that secondary metabolic pathways in rice seedlings are influenced by the interaction between UVB irradiation and the leaf developmental stage. Thus, the strategies of protection against, adaptation to, and mitigation of UVB might be delicately regulated, and, in this context, our data provide valuable information to understand UVB-induced secondary metabolism in rice seedlings.

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

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

Improved Foods Using Enzymes from Basidiomycetes

Within the kingdom of fungi, the division Basidiomycota represents more than 30,000 species, some with huge genomes indicating great metabolic potential. The fruiting bodies of many basidiomycetes are appreciated as food (“mushrooms”). Solid-state and submerged cultivation processes have been established for many species. Specifically, xylophilic fungi secrete numerous enzymes but also form smaller metabolites along unique pathways; both groups of compounds may be of interest to the food processing industry. To stimulate further research and not aim at comprehensiveness in the broad field, this review describes some recent progress in fermentation processes and the knowledge of fungal genetics. Processes with potential for food applications based on lipases, esterases, glycosidases, peptidases and oxidoreductases are presented. The formation and degradation of colourants, the degradation of harmful food components, the formation of food ingredients and particularly of volatile and non-volatile flavours serve as examples. In summary, edible basidiomycetes are foods—and catalysts—for food applications and rich donors of genes to construct heterologous cell factories for fermentation processes. Options arise to support the worldwide trend toward greener, more eco-friendly and sustainable processes.

Oxidative Stress and Autophagy Are Important Processes in Post Ripeness and Brown Film Formation in Mycelium of Lentinula edodes

Lentinula edodes (Berk.) Pegler, the shiitake mushroom, is one of the most important mushrooms in the global mushroom industry. Although mycelium post ripeness and brown film (BF) formation are crucial for fruiting body initiation, the underlying molecular mechanisms of BF formation are largely unknown. In this study, proteomic quantification (relative and absolute) and metabolomic profiling of L. edodes were performed using isobaric tags and gas chromatography-mass spectroscopy, respectively. A total of 2,474 proteins were identified, which included 239 differentially expressed proteins. Notably, several proteins associated with autophagy were upregulated, including RPD3, TOR1, VAC8, VPS1, and VPS27. Transmission electron microscopy also indicated that autophagy occurred in post ripeness and BF formation. In time-dependent analysis of the metabolome, metabolites associated with oxidative stress and autophagy changed significantly, including mannitol, trehalose, myo-inositol, glucose, leucine, valine, glutamine, and 4-aminobutyric acid. Thus, oxidative stress and autophagy were important processes in post ripeness and BF formation in L. edodes, and new insights were gained into molecular mechanisms at proteome and metabolome levels.

Chromosomal genome and population genetic analyses to reveal genetic architecture, breeding history and genes related to cadmium accumulation in Lentinula edodes

BACKGROUND

Lentinula edodes (Berk.) is the second most productive mushroom in the world. It contains compounds effective for antiviral, antitumor, antioxidant and immune regulation. Although genomes have previously been reported for this species, a high-quality chromosome-level reference for L. edodes is unavailable. This hinders detailed investigation of population genetics, breeding history of strains and genes related to environmental stress responses.

RESULTS

A high-quality chromosome-level genome was constructed. We separated a monokaryon from protoplasts of the commercial L. edodes strain L808 and assembled the genome of L. edodes using PacBio long-read and Illumina short-read sequencing, along with the high-throughput chromatin conformation capture (Hi-C) technique. We assembled a 45.87 Mb genome, and 99% of the sequences were anchored onto 10 chromosomes. The contig and scaffold N50 length were 2.17 and 4.94 Mb, respectively. Over 96% of the complete Benchmarking Universal Single-Copy Orthologs (BUSCO) were identified, and 9853 protein-coding genes were predicted. We performed population genome resequencing using 34 wild strains and 65 commercial cultivars of L. edodes originating from China, Japan, the United States and Australia. Based on whole-genome variants, we showed substantial differences in the Chinese wild population, which divided into different branches according to the main areas of their geographical distribution. We also determined the breeding history of L. edodes at the molecular level, and demonstrated that the cultivated strains in China mainly originated from wild strains from China and Northeast Asia. Phenotypic analysis showed that 99 strains exhibited differences on the Cd accumulation. Three significant loci in the of L. edodes genome were identified using the genome-wide association study (GWAS) of Cd accumulation traits. Functional genes associated with Cd accumulation traits were related to DNA ligase and aminoacyl tRNA synthetase, indicating that DNA damage repair and in vivo protein translation may be responses to Cd stress.

CONCLUSIONS

A high-quality chromosome-level genome and population genetic data of L. edodes provide genetic resources for functional genomic, evolutionary and artificial breeding studies for L. edodes.

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.

Whole Genome Sequencing and Annotation of Naematelia aurantialba (Basidiomycota, Edible-Medicinal Fungi)

Naematelia aurantialba is a rare edible fungus with both nutritional and medicinal values and especially rich in bioactive polysaccharides. However, due to the lack of genomic information, researches on the mining of active compounds, artificial breeding and cultivation, genetics, and molecular biology are limited. To facilitate the medicinal and food applications of N. aurantialba, we sequenced and analyzed the whole genome of N. aurantialba for the first time. The 21-Mb genome contained 15 contigs, and a total of 5860 protein-coding genes were predicted. The genome sequence shows that 296 genes are related to polysaccharide synthesis, including 15 genes related to nucleoside-activated sugar synthesis and 11 genes related to glucan synthesis. The genome also contains genes and gene clusters for the synthesis of other active substances, including terpenoids, unsaturated fatty acids, and bioactive proteins. In addition, it was also found that N. aurantialba was more closely related to Naematelia encephala than to Tremella fuciformis. In short, this study provides a reference for molecular cognition of N. aurantialba and related researches.

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.

Construction of a CRISPR/Cas9-Mediated Genome Editing System in Lentinula edodes

CRISPR/Cas9 genome editing systems have been established in a broad range of eukaryotic species. Herein, we report the first method for genetic engineering in pyogo (shiitake) mushrooms (Lentinula edodes) using CRISPR/Cas9. For in vivo expression of guide RNAs (gRNAs) targeting the mating-type gene HD1 (LeA1), we identified an endogenous LeU6 promoter in the L. edodes genome. We constructed a plasmid containing the LeU6 and glyceraldehyde-3-phosphate dehydrogenase (LeGPD) promoters to express the Cas9 protein. Among the eight gRNAs we tested, three successfully disrupted the LeA1 locus. Although the CRISPR-Cas9-induced alleles did not affect mating with compatible monokaryotic strains, disruption of the transcription levels of the downstream genes of LeHD1 and LeHD2 was detected. Based on this result, we present the first report of a simple and powerful genetic manipulation tool using the CRISPR/Cas9 toolbox for the scientifically and industrially important edible mushroom, L. edodes.

The Genomic and Transcriptomic Analyses of Floccularia luteovirens, a Rare Edible Fungus in the Qinghai-Tibet Plateau, Provide Insights into the Taxonomy Placement and Fruiting Body Formation

Floccularia luteovirens is a famous and precious edible mushroom (Huang Mogu) on the Qinghai–Tibet plateau that has a unique flavor and remarkable medical functions. Herein, we report a reference-grade 27 Mb genome of F. luteovirens containing 7068 protein-coding genes. The genome component and gene functions were predicted. Genome ontology enrichment and pathway analyses indicated the potential production capacity for terpenoids, polyketides and polysaccharides. Moreover, 16 putative gene clusters and 145 genes coding for secondary metabolites were obtained, including guadinomine and melleolides. In addition, phylogenetic and comparative genomic analyses shed light on the precise classification of F. luteovirens suggesting that it belongs to the genus Floccularia instead of Armillaria. RNA-sequencing and comparative transcriptomic analysis revealed differentially expressed genes during four developmental stages of F. luteovirens, that of which helps to identify important genes regulating fruiting body formation for strain modification. This study will provide insight into artificial cultivation and increase the production of useful metabolites.

Whole-genome assembly of Ganoderma leucocontextum (Ganodermataceae, Fungi) discovered from the Tibetan Plateau of China

Ganoderma leucocontextum, a newly discovered species of Ganodermataceae in China, has diverse pharmacological activities. G. leucocontextum was widely cultivated in southwest China, but the systematic genetic study has been impeded by the lack of a reference genome. Herein, we present the first whole-genome assembly of G. leucocontextum based on the Illumina and Nanopore platform from high-quality DNA extracted from a monokaryon strain (DH-8). The generated genome was 50.05 Mb in size with a N50 scaffold size of 3.06 Mb, 78,206 coding sequences and 13,390 putative genes. Genome completeness was assessed using the Benchmarking Universal Single-Copy Orthologs (BUSCO) tool, which identified 96.55% of the 280 Fungi BUSCO genes. Furthermore, differences in functional genes of secondary metabolites (terpenoids) were analyzed between G. leucocontextum and G. lucidum. G. leucocontextum has more genes related to terpenoids synthesis compared to G. lucidum, which may be one of the reasons why they exhibit different biological activities. This is the first genome assembly and annotation for G. leucocontextum, which would enrich the toolbox for biological and genetic studies in G. leucocontextum.

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.

Transcriptome of different fruiting stages in the cultivated mushroom Cyclocybe aegerita suggests a complex regulation of fruiting and reveals enzymes putatively involved in fungal oxylipin biosynthesis

BACKGROUND

Cyclocybe aegerita (syn. Agrocybe aegerita) is a commercially cultivated mushroom. Its archetypal agaric morphology and its ability to undergo its whole life cycle under laboratory conditions makes this fungus a well-suited model for studying fruiting body (basidiome, basidiocarp) development. To elucidate the so far barely understood biosynthesis of fungal volatiles, alterations in the transcriptome during different developmental stages of C. aegerita were analyzed and combined with changes in the volatile profile during its different fruiting stages.

RESULTS

A transcriptomic study at seven points in time during fruiting body development of C. aegerita with seven mycelial and five fruiting body stages was conducted. Differential gene expression was observed for genes involved in fungal fruiting body formation showing interesting transcriptional patterns and correlations of these fruiting-related genes with the developmental stages. Combining transcriptome and volatilome data, enzymes putatively involved in the biosynthesis of C8 oxylipins in C. aegerita including lipoxygenases (LOXs), dioxygenases (DOXs), hydroperoxide lyases (HPLs), alcohol dehydrogenases (ADHs) and ene-reductases could be identified. Furthermore, we were able to localize the mycelium as the main source for sesquiterpenes predominant during sporulation in the headspace of C. aegerita cultures. In contrast, changes in the C8 profile detected in late stages of development are probably due to the activity of enzymes located in the fruiting bodies.

CONCLUSIONS

In this study, the combination of volatilome and transcriptome data of C. aegerita revealed interesting candidates both for functional genetics-based analysis of fruiting-related genes and for prospective enzyme characterization studies to further elucidate the so far barely understood biosynthesis of fungal C8 oxylipins.

Transcriptional Changes of Cell Wall Organization Genes and Soluble Carbohydrate Alteration during Leaf Blade Development of Rice Seedlings

Plant cell walls have two constituent parts with different components and developmental stages. Much of the mystery concerning the mechanisms of synthesis, decomposition, modification, and so forth, has been resolved using omics and microscopic techniques. However, it still remains to be determined how cell wall development progresses over time after leaf emergence. Our focus in the present study was to expand our knowledge of the molecular mechanisms associated with cell wall synthesis in rice leaf blade during three distinct stages (sink, sink-to-source transition, and source). The RNA-seq, quantitative reverse transcription PCR (qRT-PCR) and carbohydrate concentrations were evaluated using developing fifth leaf blades harvested at different time points. The results revealed that some of the essential genes for the primary cell wall (PCW) were highly upregulated in the sink-to-source transition compared to the sink stage, whereas those essential to the secondary cell wall (SCW) displayed relatively higher levels (p < 0.05) during the source stage. The concentrations of soluble carbohydrates differed via type rather than stage; we observed higher monosaccharides during the sink stage and higher di- and oligo-saccharides during the sink-to-source transition and source stages. In conclusion, our findings suggest that the transcriptional regulation of plant cell wall biosynthesis genes are both synchronistic with and independent of, and directly and indirectly governed by, the abundance of soluble carbohydrates in the developing leaf blade, and, finally, raffinose is likely to play a transport role comparable to sucrose.

RNA-Seq-based high-resolution linkage map reveals the genetic architecture of fruiting body development in shiitake mushroom, Lentinula edodes

•

We constructed a reference genetic map of Lentinula edodes.

•

We re-assembled a chromosome-level genome of L. edodes.

•

We disclosed three hotspots regions for fruiting body-related traits in shiitake.

•

We scanned candidate genes for fruiting body-related traits.

Fruiting body development (FBD) of mushroom-forming fungi has attracted tremendous interest. However, the genetic and molecular basis of FBD is poorly known. Here, using Lentinula edodes (shiitake) as a model, we deciphered the genetic architecture underlying fruiting body-related traits (FBRTs) by combined genomic, genetic and phenotypic data. Using RNA-Seq of fruiting bodies from 110 dikaryons in a bi-parental mapping population, we constructed an ultra-high-density genetic map of L. edodes (Lemap2.0) with a total length of 810.14 cM, which covered 81.7% of the shiitake genome. A total of 94 scaffolds of the shiitake genome were aligned to Lemap2.0 and re-anchored into nine pseudo-chromosomes. Then via quantitative trait locus (QTL) analysis, we disclosed an outline of the genetic architecture of FBD in shiitake. Twenty-nine QTLs and three main genomic regions associated with FBD of shiitake were identified. Using meta-QTL analysis, seven pleiotropic QTLs for multiple traits were detected, which contributed to the correlations of FBRTs. In the mapped QTLs, the expression of 246 genes were found to significantly correlate with the phenotypic traits. Thirty-three of them were involved in FBD and could represent candidate genes controlling the shape and size of fruiting bodies. Collectively, our findings have advanced our understanding of the genetic regulation of FBD in shiitake and mushroom-forming fungi at large.

Functional Analysis of a Novel ABL (Abnormal Browning Related to Light) Gene in Mycelial Brown Film Formation of Lentinula edodes

Lentinula edodes is a globally important edible mushroom species that is appreciated for its medicinal properties as well as its nutritional value. During commercial cultivation, a mycelial brown film forms on the surface of the sawdust growth medium at the late vegetative stage. Mycelial film formation is a critical developmental process that contributes to the quantity and quality of the mushroom yield. However, little is known regarding the genetic underpinnings of brown film formation on the surface of mycelial tissue. A novel causal gene associated with the formation of the mycelial brown film, named ABL (Abnormal browning related to light), was identified in this study. The comparative genetic analysis by dihybrid crosses between normal and abnormal browning film cultivars demonstrated that a single dominant allele was responsible for the abnormal mycelium browning phenotype. Whole-genome sequencing analysis of hybrid isolates revealed five missense single-nucleotide polymorphisms (SNPs) in the ABL locus of individuals forming abnormal partial brown films. Additional whole-genome resequencing of a further 16 cultivars showed that three of the five missense SNPs were strongly associated with the abnormal browning phenotype. Overexpression of the dominant abl-D allele in a wild-type background conferred the abnormal mycelial browning phenotype upon transformants, with slender hyphae observed as a general defective mycelial growth phenotype. Our methodology will aid the future discovery of candidate genes associated with favorable traits in edible mushrooms. The discovery of a novel gene, ABL, associated with mycelial film formation will facilitate marker-associated breeding in L. edodes.

Genetically Independent Tetranucleotide to Hexanucleotide Core Motif SSR Markers for Identifying Lentinula edodes Cultivars

For the purpose of protecting the rights of Lentinula edodes breeders, we developed a new simple sequence repeat (SSR) marker set consisting only of genetically independent tetranucleotide or longer core motifs. Using available genome sequences for five L. edodes strains, we designed primers for 13 SSR markers that amplified polymorphic sequences in 20 L. edodes cultivars. We evaluated the independence of every possible marker pair based on genotype data. Consequently, eight genetically independent markers were selected. The polymorphic information content values of the markers ranged from 0.269 to 0.764, with an average of 0.409. The markers could distinguish among 20 L. edodes cultivars and produced highly repeatable and reproducible results. The markers developed in this study will enable the precise identification of L. edodes cultivars, and may be useful for protecting breeders' rights.

Genomic Analyses Provide Insights Into the Evolutionary History and Genetic Diversity of Auricularia Species

Species in the genus Auricularia play important roles for people’s food and nutrition especially Auricularia cornea and A. heimuer. To understand their evolutionary history, genome structure, and population-level genetic variation, we performed a high-quality genome sequencing of Auricularia cornea and the corresponding comparative genomic analysis. The genome size of A. cornea was similar to Auricularia subglabra, but 1.5 times larger than that of A. heimuer. Several factors were responsible for genome size variation including gene numbers, repetitive elements, and gene lengths. Phylogenomic analysis revealed that the estimated divergence time between A. heimuer and other Auricularia is ∼79.1 million years ago (Mya), while the divergence between A. cornea and A. subglabra occurred in ∼54.8 Mya. Population genomic analysis also provided insight into the demographic history of A. cornea and A. heimuer, indicating that their populations fluctuated over time with global climate change during Marine Isotope Stage 5-2. Moreover, despite the highly similar external morphologies of A. cornea and A. heimuer, their genomic properties were remarkably different. The A. cornea genome only shared 14% homologous syntenic blocks with A. heimuer and possessed more genes encoding carbohydrate-active enzymes and secondary metabolite biosynthesis proteins. The cross-taxa transferability rates of simple sequence repeat (SSR) and insertion or deletion (InDel) markers within the genus Auricularia were also lower than that previously observed for species within the same genus. Taken together, these results indicate a high level of genetic differentiation between these two Auricularia species. Consequently, our study provides new insights into the genomic evolution and genetic differentiation of Auricularia species that will facilitate future genetic breeding.

Chitinases Play a Key Role in Stipe Cell Wall Extension in the Mushroom Coprinopsis cinerea

The elongation growth of the mushroom stipe is a characteristic but not well-understood morphogenetic event of basidiomycetes. We found that extending native stipe cell walls of Coprinopsis cinerea were associated with the release of N-acetylglucosamine and chitinbiose and with chitinase activity. Two chitinases among all detected chitinases from C. cinerea, ChiE1 and ChiIII, reconstituted heat-inactivated stipe wall extension and released N-acetylglucosamine and chitinbiose. Interestingly, both ChiE1 and ChiIII hydrolyze insoluble crystalline chitin powder, while other C. cinerea chitinases do not, suggesting that crystalline chitin components of the stipe cell wall are the target of action for ChiE1 and ChiIII. ChiE1- or ChiIII-reconstituted heat-inactivated stipe walls showed maximal extension activity at pH 4.5, consistent with the optimal pH for native stipe wall extension in vitro; ChiE1- or ChiIII-reconstituted heat-inactivated stipe wall extension activities were associated with stipe elongation growth regions; and the combination of ChiE1 and ChiIII showed a synergism to reconstitute heat-inactivated stipe wall extension at a low action concentration. Field emission scanning electron microscopy (FESEM) images showed that the inner surface of acid-induced extended native stipe cell walls and ChiE1- or ChiIII-reconstituted extended heat-inactivated stipe cell walls exhibited a partially broken parallel microfibril architecture; however, these broken transversely arranged microfibrils were not observed in the unextended stipe cell walls that were induced by neutral pH buffer or heat inactivation. Double knockdown of ChiE1 and ChiIII resulted in the reduction of stipe elongation, mycelium growth, and heat-sensitive cell wall extension of native stipes. These results indicate a chitinase-hydrolyzing mechanism for stipe cell wall extension.IMPORTANCE A remarkable feature in the development of basidiomycete fruiting bodies is stipe elongation growth that results primarily from manifold cell elongation. Some scientists have suggested that stipe elongation is the result of enzymatic hydrolysis of cell wall polysaccharides, while other scientists have proposed the possibility that stipe elongation results from nonhydrolytic disruption of the hydrogen bonds between cell wall polysaccharides. Here, we show direct evidence for a chitinase-hydrolyzing mechanism of stipe cell wall elongation in the model mushroom Coprinopsis cinerea that is different from the expansin nonhydrolysis mechanism of plant cell wall extension. We presumed that in the growing stipe cell walls, parallel chitin microfibrils are tethered by β-1,6-branched β-1,3-glucans, and that the breaking of the tether by chitinases leads to separation of these microfibrils to increase their spacing for insertion of new synthesized chitin and β-1,3-glucans under turgor pressure in vivo.

Variable Number Tandem Repeats in the Mitochondrial DNA of Lentinula edodes

Variable number tandem repeats (VNTRs) in mitochondrial DNA (mtDNA) of Lentinula edodes are of interest for their role in mtDNA variation and their application as genetic marker. Sequence analysis of three L. edodes mtDNAs revealed the presence of VNTRs of two categories. Type I VNTRs consist of two types of repeat units in a symmetric distribution, whereas Type II VNTRs contain tandemly arrayed repeats of 7- or 17-bp DNA sequences. The number of repeat units was variable depending on the mtDNA of different strains. Using the variations in VNTRs as a mitochondrial marker and the A mating type as a nuclear type marker, we demonstrated that one of the two nuclei in the donor dikaryon preferentially enters into the monokaryotic cytoplasm to establish a new dikaryon which still retains the mitochondria of the monokaryon in the individual mating. Interestingly, we found 6 VNTRs with newly added repeat units from the 22 mates, indicating that elongation of VNTRs occurs during replication of mtDNA. This, together with comparative analysis of the repeating pattern, enables us to propose a mechanistic model that explains the elongation of Type I VNTRs through reciprocal incorporation of basic repeat units, 5’-TCCCTTTAGGG-3’ and its complementary sequence (5’-CCCTAAAGGGA-3’).

Whole genome sequencing and genome annotation of the wild edible mushroom, Russula griseocarnosa

Russula griseocarnosa is a species of edible ectomycorrhizal fungi with medicinal properties that grows in southern China. Total DNA was isolated from a fresh fruiting body of R. griseocarnosa and subjected to sequencing using Illumina Hiseq with the PacBio RS sequencing platform. Here, we present the 64.81 Mb draft genome map of R. griseocarnosa based on 471 scaffolds and 16,128 coding protein genes. The gene annotation of protein coding genes was used to obtain corresponding annotations by blastp. Phylogenetic analysis revealed a close evolutionary relationship of R. griseocarnosa to Heterobasidion irregulare and Stereum hirsutum in the core Russulales clade. The R. griseocarnosa genome encodes a repertoire of enzymes engaged in carbohydrate and polysaccharide metabolism, along with cytochrome P450s and secondary metabolite biosynthesis. The genome content of R. griseocarnosa provides insights into the genetic basis of its reported medicinal properties and serves as a reference for comparative genomics of fungi.

Genome Sequencing Illustrates the Genetic Basis of the Pharmacological Properties of Gloeostereum incarnatum

Gloeostereum incarnatum is a precious edible mushroom that is widely grown in Asia and known for its useful medicinal properties. Here, we present a high-quality genome of G. incarnatum using the single-molecule real-time (SMRT) sequencing platform. The G. incarnatum genome, which is the first complete genome to be sequenced in the family Cyphellaceae, was 38.67 Mbp, with an N50 of 3.5 Mbp, encoding 15,251 proteins. Based on our phylogenetic analysis, the Cyphellaceae diverged ~174 million years ago. Several genes and gene clusters associated with lignocellulose degradation, secondary metabolites, and polysaccharide biosynthesis were identified in G. incarnatum, and compared with other medicinal mushrooms. In particular, we identified two terpenoid-associated gene clusters, each containing a gene encoding a sesterterpenoid synthase adjacent to a gene encoding a cytochrome P450 enzyme. These clusters might participate in the biosynthesis of incarnal, a known bioactive sesterterpenoid produced by G. incarnatum. Through a transcriptomic analysis comparing the G. incarnatum mycelium and fruiting body, we also demonstrated that the genes associated with terpenoid biosynthesis were generally upregulated in the mycelium, while those associated with polysaccharide biosynthesis were generally upregulated in the fruiting body. This study provides insights into the genetic basis of the medicinal properties of G. incarnatum, laying a framework for future characterization of bioactive proteins and pharmaceutical uses of this fungus.

Comparative transcriptome analysis identified candidate genes involved in mycelium browning in Lentinula edodes

Background

Lentinula edodes is one of the most popular edible mushroom species in the world and contains useful medicinal components, such as lentinan. The light-induced formation of brown film on the vegetative mycelial tissues of L. edodes is an important process for ensuring the quantity and quality of this edible mushroom. To understand the molecular mechanisms underlying this critical developmental process in L. edodes, we characterized the morphological phenotypic changes in a strain, Chamaram, associated with abnormal brown film formation and compared its genome-wide transcriptional features.

Results

In the present study, we performed genome-wide transcriptome analyses of different vegetative mycelium growth phenotypes, namely, early white, normal brown, and defective dark yellow partial brown films phenotypes which were exposed to different light conditions. The analysis revealed the identification of clusters of genes specific to the light-induced brown film phenotypes. These genes were significantly associated with light sensing via photoreceptors such as FMN- and FAD-bindings, signal transduction by kinases and GPCRs, melanogenesis via activation of tyrosinases, and cell wall degradation by glucanases, chitinases, and laccases, which suggests these processes are involved in the formation of mycelial browning in L. edodes. Interestingly, hydrophobin genes such as SC1 and SC3 exhibited divergent expression levels in the normal and abnormal brown mycelial films, indicating the ability of these genes to act in fruiting body initiation and formation of dikaryotic mycelia. Furthermore, we identified the up-regulation of glycoside hydrolase domain-containing genes in the normal brown film but not in the abnormal film phenotype, suggesting that cell wall degradation in the normal brown film phenotype is crucial in the developmental processes related to the initiation and formation of fruiting bodies.

Conclusions

This study systematically analysed the expression patterns of light-induced browning-related genes in L. edodes. Our findings provide information for further investigations of browning formation mechanisms in L. edodes and a foundation for future L. edodes breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5509-4) contains supplementary material, which is available to authorized users.

TaF: a web platform for taxonomic profile-based fungal gene prediction

INTRODUCTION

The accurate prediction and annotation of gene structures from the genome sequence of an organism enable genome-wide functional analyses to obtain insight into the biological properties of an organism.

OBJECTIVES

We recently developed a highly accurate filamentous fungal gene prediction pipeline and web platform called TaF. TaF is a homology-based gene predictor employing large-scale taxonomic profiling to search for close relatives in genome queries.

METHODS

TaF pipeline consists of four processing steps; (1) taxonomic profiling to search for close relatives to query, (2) generation of hints for determining exon-intron boundaries from orthologous protein sequence data of the profiled species, (3) gene prediction by combination of ab inito and evidence-based prediction methods, and (4) homology search for gene models.

RESULTS

TaF generates extrinsic evidence that suggests possible exon-intron boundaries based on orthologous protein sequence data, thus reducing false-positive predictions of gene structure based on distantly related orthologs data. In particular, the gene prediction method using taxonomic profiling shows very high accuracy, including high sensitivity and specificity for gene models, suggesting a new approach for homology-based gene prediction from newly sequenced or uncharacterized fungal genomes, with the potential to improve the quality of gene prediction.

CONCLUSION

TaF will be a useful tool for fungal genome-wide analyses, including the identification of targeted genes associated with a trait, transcriptome profiling, comparative genomics, and evolutionary analysis.

Comparative transcriptome analysis of dikaryotic mycelia and mature fruiting bodies in the edible mushroom Lentinula edodes

Lentinula edodes is a popular cultivated edible mushroom with high nutritional and medicinal value. To understand the regulation of gene expression in the dikaryotic mycelium and mature fruiting body in the commercially important Korean L. edodes strain, we first performed comparative transcriptomic analysis, using Illumina HiSeq platform. De novo assembly of these sequences revealed 11,675 representative transcripts in two different stages of L. edodes. A total of 9,092 unigenes were annotated and subjected to Gene Ontology, EuKaryotic Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Gene expression analysis revealed that 2,080 genes were differentially expressed, with 1,503 and 577 upregulated in the mycelium and a mature fruiting body, respectively. Analysis of 18 KEGG categories indicated that fruiting body-specific transcripts were significantly enriched in ‘replication and repair’ and ‘transcription’ pathways, which are important for premeiotic replication, karyogamy, and meiosis during maturation. We also searched for fruiting body-specific proteins such as aspartic protease, gamma-glutamyl transpeptidase, and cyclohexanone monooxygenase, which are involved in fruiting body maturation and isolation of functional substances. These transcriptomes will be useful in elucidating the molecular mechanisms of mature fruiting body development and beneficial properties, and contribute to the characterization of novel genes in L. edodes.

Nucleus-Selective Expression of Laccase Genes in the Dikaryotic Strain of Lentinula edodes

In mating of Lentinula edodes, dikaryotic strains generated from certain monokaryotic strains such as the B2 used in this study tend to show better quality of fruiting bodies regardless of the mated monokaryotic strains. Unlike B2, dikaryotic strains generated from B16 generally show low yields, with deformed or underdeveloped fruiting bodies. This indicates that the two nuclei in the cytoplasm do not contribute equally to the physiology of dikaryotic L. edodes, suggesting an expression bias in the allelic genes of the two nuclei. To understand the role of each nucleus in dikaryotic strains, we investigated single nucleotide polymorphisms (SNPs) in laccase genes of monokaryotic strains to reveal nuclear origin of the expressed mRNAs in dikaryotic strain. We performed reverse transcription PCR (RT-PCR) analysis using total RNAs extracted from dikaryotic strains (A5B2, A18B2, and A2B16) as well as from compatible monokaryotic strains (A5, A18, and B2 for A5B2 and A18B2; A2 and B16 for A2B16). RT-PCR results revealed that Lcc1, Lcc2, Lcc4, Lcc7, and Lcc10 were the mainly expressed laccase genes in the L. edodes genome. To determine the nuclear origin of these laccase genes, the genomic DNA sequences in monokaryotic strains were analyzed, thereby revealing five SNPs in Lcc4 and two in Lcc7. Subsequent sequence analysis of laccase mRNAs expressed in dikaryotic strains revealed that these were almost exclusively expressed from B2-originated nuclei in A5B2 and A18B2 whereas B16 nucleus did not contribute to laccase expression in A2B16 strain. This suggests that B2 nucleus dominates the expression of allelic genes, thereby governing the physiology of dikaryons.

Long-read transcriptome data for improved gene prediction in Lentinula edodes

Lentinula edodes is one of the most popular edible mushrooms in the world and contains useful medicinal components such as lentinan. The whole-genome sequence of L. edodes has been determined with the objective of discovering candidate genes associated with agronomic traits, but experimental verification of gene models with correction of gene prediction errors is lacking. To improve the accuracy of gene prediction, we produced 12.6 Gb of long-read transcriptome data of variable lengths using PacBio single-molecule real-time (SMRT) sequencing and generated 36,946 transcript clusters with an average length of 2.2 kb. Evidence-driven gene prediction on the basis of long- and short-read RNA sequencing data was performed; a total of 16,610 protein-coding genes were predicted with error correction. Of the predicted genes, 42.2% were verified to be covered by full-length transcript clusters. The raw reads have been deposited in the NCBI SRA database under accession number PRJNA396788.

Lentinula edodes Genome Survey and Postharvest Transcriptome Analysis

Lentinula edodes is a popular, cultivated edible and medicinal mushroom. Lentinula edodes is susceptible to postharvest problems, such as gill browning, fruiting body softening, and lentinan degradation. We constructed a de novo assembly draft genome sequence and performed gene prediction for Lentinula edodesDe novo assembly was carried out using short reads from paired-end and mate-paired libraries and by using long reads by PacBio, resulting in a contig number of 1,951 and an N₅₀ of 1 Mb. Furthermore, we predicted genes by Augustus using transcriptome sequencing (RNA-seq) data from the whole life cycle of Lentinula edodes, resulting in 12,959 predicted genes. This analysis revealed that Lentinula edodes lacks lignin peroxidase. To reveal genes involved in the loss of quality of Lentinula edodes postharvest fruiting bodies, transcriptome analysis was carried out using serial analysis of gene expression (SuperSAGE). This analysis revealed that many cell wall-related enzymes are upregulated after harvest, such as β-1,3-1,6-glucan-degrading enzymes in glycoside hydrolase (GH) families GH5, GH16, GH30, GH55, and GH128, and thaumatin-like proteins. In addition, we found that several chitin-related genes are upregulated, such as putative chitinases in GH family 18, exochitinases in GH20, and a putative chitosanase in GH family 75. The results suggest that cell wall-degrading enzymes synergistically cooperate for rapid fruiting body autolysis. Many putative transcription factor genes were upregulated postharvest, such as genes containing high-mobility-group (HMG) domains and zinc finger domains. Several cell death-related proteins were also upregulated postharvest.IMPORTANCE Our data collectively suggest that there is a rapid fruiting body autolysis system in Lentinula edodes The genes for the loss of postharvest quality newly found in this research will be targets for the future breeding of strains that keep fresh longer than present strains. De novoLentinula edodes genome assembly data will be used for the construction of a complete Lentinula edodes chromosome map for future breeding.

Development of 44 Novel Polymorphic SSR Markers for Determination of Shiitake Mushroom (Lentinula edodes) Cultivars

The shiitake mushroom (Lentinulaedodes) is one of the most popular edible mushrooms in the world and has attracted attention for its value in medicinal and pharmacological uses. With recent advanced research and techniques, the agricultural cultivation of the shiitake mushroom has been greatly increased, especially in East Asia. Additionally, demand for the development of new cultivars with good agricultural traits has been greatly enhanced, but the development processes are complicated and more challenging than for other edible mushrooms. In this study, we developed 44 novel polymorphic simple sequence repeat (SSR) markers for the determination of shiitake mushroom cultivars based on a whole genome sequencing database of L. edodes. These markers were found to be polymorphic and reliable when screened in 23 shiitake mushroom cultivars. For the 44 SSR markers developed in this study, the major allele frequency ranged from 0.13 to 0.94; the number of genotypes and number of alleles were each 2-11; the observed and expected heterozygosity were 0.00-1.00 and 0.10-0.90, respectively; and the polymorphic information content value ranged from 0.10 to 0.89. These new markers can be used for molecular breeding, the determination of cultivars, and other applications.

Genome Sequence of the Edible Cultivated Mushroom Lentinula edodes (Shiitake) Reveals Insights into Lignocellulose Degradation

Lentinula edodes, one of the most popular, edible mushroom species with a high content of proteins and polysaccharides as well as unique aroma, is widely cultivated in many Asian countries, especially in China, Japan and Korea. As a white rot fungus with lignocellulose degradation ability, L. edodes has the potential for application in the utilization of agriculture straw resources. Here, we report its 41.8-Mb genome, encoding 14,889 predicted genes. Through a phylogenetic analysis with model species of fungi, the evolutionary divergence time of L. edodes and Gymnopus luxurians was estimated to be 39 MYA. The carbohydrate-active enzyme genes in L. edodes were compared with those of the other 25 fungal species, and 101 lignocellulolytic enzymes were identified in L. edodes, similar to other white rot fungi. Transcriptome analysis showed that the expression of genes encoding two cellulases and 16 transcription factor was up-regulated when mycelia were cultivated for 120 minutes in cellulose medium versus glucose medium. Our results will foster a better understanding of the molecular mechanism of lignocellulose degradation and provide the basis for partial replacement of wood sawdust with agricultural wastes in L. edodes cultivation.