Interspecific hybridization between Ganoderma lingzhi and G. applanatum through protoplast fusion.

Interspecific hybridization between Ganoderma lingzhi and G. applanatum was attempted through polyethylene glycol (PEG) induced fusion technique. The protoplast isolation procedure was simplified, and we obtained a significant number of protoplasts from both Ganoderma species. The number of protoplasts obtained was 5.27 ± 0.31 × 107/mL in G. lingzhi and 5.57 ± 0.49 × 106/mL in G. applanatum. Osmotic stabilizer NaCl (0.4 M) at pH 5.8 and enzymolysis time 3.5 h have supported high frequency of protoplast regeneration. G. lingzhi and G. applanatum regeneration frequency was 1.73 ± 0.04% and 0.23 ± 0.02%, respectively. 40% of PEG induced high number of protoplast fusion the regeneration frequency was 0.09% on a minimal medium. Two hundred fifty-two fusant colonies were isolated from the following four individual experiments. Among them, ten fusants showed the mycelial morphological difference compared to their parents and other fusant isolates. The fruiting body could be generated on oak sawdust and wheat bran substrate, and a few of them showed recombined morphology of the parental strains. The highest yield and biological efficacy (BE) were recorded in GF248, while least in GF244. The hybridity of the fusant was established based on mycelia, fruiting morphology, and PCR fingerprinting. ISSR and RAPD profile analysis of ten fusants and parents depicted that fusants contained polymorphic bands, which specified the rearrangement and deletion of DNA in the fusants. A Dendrogram was constructed based on the RAPD profile, and the clustering data exhibited two major clusters: cluster I included the G. lingzhi and Cluster II, including the G. applanatum and fusant lines. Total polysaccharide (α, ß and total glucan) content was compared with fusants and parental strains. The present study highlighted the efficient methods for protoplast isolation from Ganoderma species. PEG-induced fusants showed high polymorphic frequency index, while the phenotypic characters showed high similarity to G. applanatum. A significant difference was observed in the mushroom yield and its total polysaccharide between the fusants and parental strains.

Genomic discovery of the hypsin gene and biosynthetic pathways for terpenoids in Hypsizygus marmoreus.

BACKGROUND: Hypsizygus marmoreus (Beech mushroom) is a popular ingredient in Asian cuisine. The medicinal effects of its bioactive compounds such as hypsin and hypsiziprenol have been reported, but the genetic basis or biosynthesis of these components is unknown. RESULTS: In this study, we sequenced a reference strain of H. marmoreus (Haemi 51,987-8). We evaluated various assembly strategies, and as a result the Allpaths and PBJelly produced the best assembly. The resulting genome was 42.7 Mbp in length and annotated with 16,627 gene models. A putative gene (Hypma_04324) encoding the antifungal and antiproliferative hypsin protein with 75% sequence identity with the previously known N-terminal sequence was identified. Carbohydrate active enzyme analysis displayed the typical feature of white-rot fungi where auxiliary activity and carbohydrate-binding modules were enriched. The genome annotation revealed four terpene synthase genes responsible for terpenoid biosynthesis. From the gene tree analysis, we identified that terpene synthase genes can be classified into six clades. Four terpene synthase genes of H. marmoreus belonged to four different groups that implies they may be involved in the synthesis of different structures of terpenes. A terpene synthase gene cluster was well-conserved in Agaricomycetes genomes, which contained known biosynthesis and regulatory genes. CONCLUSIONS: Genome sequence analysis of this mushroom led to the discovery of the hypsin gene. Comparative genome analysis revealed the conserved gene cluster for terpenoid biosynthesis in the genome. These discoveries will further our understanding of the biosynthesis of medicinal bioactive molecules in this edible mushroom.

Genome Sequencing and Carbohydrate-Active Enzyme (CAZyme) Repertoire of the White Rot Fungus Flammulina elastica.

Next-generation sequencing (NGS) of the Flammulina elastica (wood-rotting basidiomycete) genome was performed to identify carbohydrate-active enzymes (CAZymes). The resulting assembly (31 kmer) revealed a total length of 35,045,521 bp (49.7% GC content). Using the AUGUSTUS tool, 12,536 total gene structures were predicted by ab initio gene prediction. An analysis of orthologs revealed that 6806 groups contained at least one F. elastica protein. Among the 12,536 predicted genes, F. elastica contained 24 species-specific genes, of which 17 genes were paralogous. CAZymes are divided into five classes: glycoside hydrolases (GHs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), glycosyltransferases (GTs), and auxiliary activities (AA). In the present study, annotation of the predicted amino acid sequences from F. elastica genes using the dbCAN CAZyme database revealed 508 CAZymes, including 82 AAs, 218 GHs, 89 GTs, 18 PLs, 59 CEs, and 42 carbohydrate binding modules in the F. elastica genome. Although the CAZyme repertoire of F. elastica was similar to those of other fungal species, the total number of GTs in F. elastica was larger than those of other basidiomycetes. This genome information elucidates newly identified wood-degrading machinery in F. elastica, offers opportunities to better understand this fungus, and presents possibilities for more detailed studies on lignocellulosic biomass degradation that may lead to future biotechnological and industrial applications.

A detailed analysis of the recombination landscape of the button mushroom Agaricus bisporus var. bisporus.

The button mushroom (Agaricus bisporus) is one of the world's most cultivated mushroom species, but in spite of its economic importance generation of new cultivars by outbreeding is exceptional. Previous genetic analyses of the white bisporus variety, including all cultivars and most wild isolates revealed that crossing over frequencies are low, which might explain the lack of introducing novel traits into existing cultivars. By generating two high quality whole genome sequence assemblies (one de novo and the other by improving the existing reference genome) of the first commercial white hybrid Horst U1, a detailed study of the crossover (CO) landscape was initiated. Using a set of 626 SNPs in a haploid offspring of 139 single spore isolates and whole genome sequencing on a limited number of homo- and heterokaryotic single spore isolates, we precisely mapped all COs showing that they are almost exclusively restricted to regions of about 100kb at the chromosome ends. Most basidia of A. bisporus var. bisporus produce two spores and pair preferentially via non-sister nuclei. Combined with the COs restricted to the chromosome ends, these spores retain most of the heterozygosity of the parent thus explaining how present-day white cultivars are genetically so close to the first hybrid marketed in 1980. To our knowledge this is the first example of an organism which displays such specific CO landscape.

Evaluation of anti-atopic dermatitis activity of Hypsizigus marmoreus extract.

Hypsizigus marmoreus is an edible mushroom that is cultivated worldwide. In this study, we investigated antiinflammatory activities of H. marmoreus extract on atopic dermatitis (AD)-like symptoms. Ethanol extract of H. marmoreus (HMEE) was administrated in powder to BALB/c mice in which AD was induced by 1-chloro 2, 4, 6-trinitrobenzene [picryl-chloride, (PCL)]. The dermatitis severity score and the thickness of the epidermis were significantly decreased following daily intake of HMEE powder (1 g/kg/day) for 5 weeks compared with a PCL-treated group. The mRNA expression of proinflammatory cytokines, such as interleukin-1 beta (IL-1ß) and interferon-gamma (IFN-γ), was significantly attenuated in the dorsal skin of the HMEE-fed mouse group compared with the PCL-treated mouse group. In addition, in concanavalin A-stimulated and lipopolysaccharide-stimulated mouse splenocytes and macrophages, levels of IL-1ß and IFN-γ production were attenuated following the addition of HMEE. Interestingly, the administration of HMEE to mouse splenocytes stimulated the production of an antiinflammatory cytokine, IL-4. However, increases in the levels of proinflammatory cytokines and nitric oxide were attenuated by treating the mouse splenocytes, mouse macrophages, and Raw 264.7 cell line with HMEE. These results strongly suggest that HMEE exhibits anti-AD activity via the regulation of inflammatory responses.

Improved accuracy of geographical origin identification of shiitake grown in sawdust medium: A compound-specific isotope model-based pilot study.

In countries like South Korea and the USA, origin labeling of shiitake grown using imported Chinese-inoculated medium is an issue. Therefore, we evaluated the use of compound-specific isotope analysis (CSIA) for the accurate identification of the geographical origin of shiitake (Korean, Chinese-inoculated medium, and Chinese); Chinese-inoculated medium shiitake were cultivated in Korea using inoculated sawdust medium from China. The CSIA-discriminant model showed an overall accuracy of 100% in the geographical classification of the original set and 96.4% for the cross-validated set. Glutamate and aspartate δ15N values were the most important variables for differentiating shiitake based on their origins. Compared to that observed upon using the bulk stable isotope analysis, the CSIA model was associated with significantly improved predictability of origin identification. Our findings elucidate the importance of isotope signatures in developing a reliable origin labeling method for shiitake cultured on the sawdust medium for the global market.

Identification of differentially expressed genes in Flammulina velutipes with anti-tyrosinase activity.

It was previously shown that fruiting bodies of the mushroom, Flammulina velutipes, exert anti-tyrosinase activity by producing a triacylglycerol characterized as 1',3'-dilinolenoyl-2'-linoleoylglycerol (LnLLn). In this study, we provide evidence that the mycelia of F. velutipes grown on glucose, but not on glycerol, exhibit anti-tyrosinase activity. To identify genes possibly involved in the process of expressing anti-tyrosinase activity in F. velutipes, a reverse transcription-polymerase chain reaction (RT-PCR) method that involves annealing control primers (ACPs) was employed. Using 120 ACPs, a total of 84 differentially expressed genes (DEGs) in F. velutipes mycelia with anti-tyrosinase activity were cloned and sequenced. Basic Local Alignment Search Tool (BLAST) searches revealed that 72 of the genes have known sequence homology. Of these, three genes involved in the synthesis and modification of fatty acids were selected and further quantified by real-time RT-PCR. These genes involved in fatty acid-biosynthetic processes can be potential candidates for further studies related to the development of a new anti-browning agent.

An origin identification model for labeling of shiitake (Lentinula edodes).

With the increasing globalization of the food trade across countries and continents, reliable identification of the geographical origin of products is critical. In this study, we describe the limitations of the current origin labeling system for non-soil-based agricultural products and suggest alternative strategies for the identification of the geographical origin of such products. An origin identification model based on stable isotope ratio analysis combined with discriminant analysis is used to evaluate the similarities and dissimilarities between domestic and foreign shiitake mushrooms, including Chinese inoculated sawdust blocks and Chinese origin. The results show a classification sensitivity of 92.0%, classification specificity of 91.5%, and overall accuracy of 93.5%. In particular, δ15N was the most important isotope marker for the identification of the origin of shiitake mushrooms. Hence, the current origin labeling system for mushroom species has to be revised to establish fair trade and avoid improper origin labeling in the global shiitake market.

Comparative Analysis of Carbohydrate Active Enzymes in the Flammulina velutipes var. lupinicola Genome.

The purpose of this study was to determine the genome sequence of Flammulina velutipes var. lupinicola based on next-generation sequencing (NGS) and to identify the genes encoding carbohydrate-active enzymes (CAZymes) in the genome. The optimal assembly (71 kmer) based on ABySS de novo assembly revealed a total length of 33,223,357 bp (49.53% GC content). A total of 15,337 gene structures were identified in the F. velutipes var. lupinicola genome using ab initio gene prediction method with Funannotate pipeline. Analysis of the orthologs revealed that 11,966 (96.6%) out of the 15,337 predicted genes belonged to the orthogroups and 170 genes were specific for F. velutipes var. lupinicola. CAZymes are divided into six classes: auxiliary activities (AAs), glycosyltransferases (GTs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), glycoside hydrolases (GHs), and carbohydrate-binding modules (CBMs). A total of 551 genes encoding CAZymes were identified in the F. velutipes var. lupinicola genome by analyzing the dbCAN meta server database (HMMER, Hotpep, and DIAMOND searches), which consisted of 54-95 AAs, 145-188 GHs, 55-73 GTs, 6-19 PLs, 13-59 CEs, and 7-67 CBMs. CAZymes can be widely used to produce bio-based products (food, paper, textiles, animal feed, and biofuels). Therefore, information about the CAZyme repertoire of the F. velutipes var. lupinicola genome will help in understanding the lignocellulosic machinery and in-depth studies will provide opportunities for using this fungus for biotechnological and industrial applications.

Evaluating Genetic Diversity of Agaricus bisporus Accessions through Phylogenetic Analysis Using Single-Nucleotide Polymorphism (SNP) Markers.

Agaricus bisporus, commonly known as the button mushroom, is widely cultivated throughout the world. To breed new strains with more desirable traits and improved adaptability, diverse germplasm, including wild accessions, is a valuable genetic resource. To better understand the genetic diversity available in A. bisporus and identify previously unknown diversity within accessions, a phylogenetic analysis of 360 Agaricus spp. accessions using single-nucleotide polymorphism genotyping was performed. Genetic relationships were compared using principal coordinate analysis (PCoA) among accessions with known origins and accessions with limited collection data. The accessions clustered into four groups based on the PCoA with regard to genetic relationships. A subset of 67 strains, which comprised a core collection where repetitive and uninformative accessions were not included, clustered into 7 groups following analysis. Two of the 170 accessions with limited collection data were identified as wild germplasm. The core collection allowed for the accurate analysis of A. bisporus genetic relationships, and accessions with an unknown pedigree were effectively grouped, allowing for origin identification, by PCoA analysis in this study.

Unusual genome expansion and transcription suppression in ectomycorrhizal Tricholoma matsutake by insertions of transposable elements.

Genome sequencing of Tricholoma matsutake revealed its unusually large size as 189.0 Mbp, which is a consequence of extraordinarily high transposable element (TE) content. We identified that 702 genes were surrounded by TEs, and 83.2% of these genes were not transcribed at any developmental stage. This observation indicated that the insertion of TEs alters the transcription of the genes neighboring these TEs. Repeat-induced point mutation, such as C to T hypermutation with a bias over "CpG" dinucleotides, was also recognized in this genome, representing a typical defense mechanism against TEs during evolution. Many transcription factor genes were activated in both the primordia and fruiting body stages, which indicates that many regulatory processes are shared during the developmental stages. Small secreted protein genes (<300 aa) were dominantly transcribed in the hyphae, where symbiotic interactions occur with the hosts. Comparative analysis with 37 Agaricomycetes genomes revealed that IstB-like domains (PF01695) were conserved across taxonomically diverse mycorrhizal genomes, where the T. matsutake genome contained four copies of this domain. Three of the IstB-like genes were overexpressed in the hyphae. Similar to other ectomycorrhizal genomes, the CAZyme gene set was reduced in T. matsutake, including losses in the glycoside hydrolase genes. The T. matsutake genome sequence provides insight into the causes and consequences of genome size inflation.

Fatty Acids and Stable Isotope Ratios in Shiitake Mushrooms (Lentinula edodes) Indicate the Origin of the Cultivation Substrate Used: A Preliminary Case Study in Korea.

Shiitake mushroom (Lentinula edodes) is commonly consumed worldwide and is cultivated in many farms in Korea using Chinese substrates owing to a lack of knowledge on how to prepare sawdust-based substrate blocks (bag cultivation). Consequently, issues related to the origin of the Korean or Chinese substrate used in shiitake mushrooms produced using bag cultivation have been reported. Here, we investigated differences in fatty acids (FAs) and stable isotope ratios (SIRs) in shiitake mushrooms cultivated using Korean and Chinese substrates under similar conditions (strain, temperature, humidity, etc.) and depending on the harvesting cycle. The total FA level decreased significantly by 5.49 mg∙g-1 as the harvesting cycle increased (p < 0.0001); however, no differences were found in FAs between shiitake mushrooms cultivated using Korean and Chinese substrates. Linoleic acid was the most abundant FA, accounting for 77-81% of the total FAs during four harvesting cycles. Moreover, the SIRs differed significantly between the Korean and Chinese substrates, and the harvesting cycles resulted in smaller maximum differences in SIR values compared to those of the cultivation substrate origins. Our findings contribute to the identification of the geographical origin of shiitake mushrooms and may have potential applications in international shiitake-mushroom markets.

Bioactive Sterol Derivatives Isolated from the Pleurotus djamor var. Roseus Induced Apoptosis in Cancer Cell Lines.

OBJECTIVE: The aim of the present study is to isolate and characterize the bioactive compounds from Pleurotus djamor against human breast cancer (MDA-MD-231) and mouse T cell lymphoma (EL4) cell lines. MATERIALS AND METHODS: Sequential fractionization and column chromatography methods were involved in compound isolation. The structures of the isolated compound were determined by NMR, GC/MS, and X-ray crystallography studies. RESULTS: The isolated compounds 1- 4 [D-mannitol (C1), ergosta-5,7,22-trien-3ß-ol (C2), 5,8- epidioxy-ergosta-6-22-dien-3ß-ol (C3), and palmitic acid (C4)] are white crystal and amorphous powder in nature. All these compounds were isolated from this mushroom for the first time. In vitro lipid peroxidation activities of isolated compounds were determined by ferric thiocyanate (FTC) and thiobarbituric acid (TBA) method. The sterol derivatives C2 and C3 compounds displayed strong antioxidant activity and were not significantly different (p<0.05) to α-tocopherol. This finding elaborates on the isolation of a cytotoxic compound C2 and C3 from P. djamor via a rapid elution method. CONCLUSION: The compound C3 has exhibited better cytotoxic activity against MDA-MD-231 and EL4 cells. The present finding and data might provide new insights into the possible therapeutic and pharmaceutical use for the design of anti-cancer drugs from this edible mushroom.

Isolation of 1',3'-dilinolenoyl'-2'-linoleoylglycerol with tyrosinase inhibitory activity from Flammulina velutipes.

This study was carried out to evaluate the inhibitory effect of Flammulina velutipes extracts on tyrosinase activity and to identify its biologically active component. The ethyl acetate and n-butanol extracts showed potent tyrosinase inhibitory activities. Subsequently, fractions of the nbutanol extract showed only a partial tyrosinase inhibitory activity. The most active compound of tyrosinase inhibitory activity was identified from the ethyl acetate extract as 1',3'-dilinolenoyl-2'-linoleoylglycerol (LnLLn) by comparing its mass, 1H-, and 13C-NMR spectral data with those previously reported in the literature. LnLLn showed tyrosinase inhibitory activity with an IC50 value of 16.1 mug/ml. These results suggest that the ethyl acetate extract of F. velutipes could be applicable for the development of a new whitening agent.

Comparative analysis of expressed sequence tags from Flammulina velutipes at different developmental stages.

Flammulina velutipes is a popular edible basidiomycete mushroom found in East Asia and is commonly known as winter mushroom. Mushroom development showing dramatic morphological changes by different environmental factors is scientifically and commercially interesting. To create a genetic database and isolate genes regulated during mushroom development, cDNA libraries were constructed from three developmental stages of mycelium, primordium, and fruit body in F. velutipes. We generated a total of 5,431 expressed sequence tags (ESTs) from randomly selected clones from the three cDNA libraries. Of these, 3,332 different unique genes (unigenes) were consistent with 2,442 (73%) singlets and 890 (27%) contigs. This corresponds to a redundancy of 39%. Using a homology search in the gene ontology database, the EST unigenes were classified into the three categories of molecular function (28%), biological process (29%), and cellular component (6%). Comparative analysis found great variations in the unigene expression pattern among the three different unigene sets generated from the cDNA libraries of mycelium, primordium, and fruit body. The 19-34% of total unigenes were unique to each unigene set and only 3% were shared among all three unigene sets. The unique and common representation in F. velutipes unigenes from the three different cDNA libraries suggests great differential gene expression profiles during the different developmental stages of F. velutipes mushroom.

Development and characterization of new microsatellite markers for the oyster mushroom (Pleurotus ostreatus).

We developed and characterized 36 polymorphic microsatellite markers for the oyster mushroom (Pleurotus ostreatus). In total, 169 alleles were identified with an average of 4.7 alleles per locus. Values for observed (HO) and expected (HE) heterozygosities ranged from 0.027 to 0.946 and from 0.027 to 0.810, respectively. Nineteen loci deviated from Hardy-Weinberg equilibrium. Significant (P<0.05) excess heterozygosity was observed at nine loci. Linkage disequilibrium (LD) was significant (P<0.05) between pairs of locus alleles. Cluster analysis revealed that five species of genus Pleurotus made a distinct group, and the individual cultivars were grouped into major five groups from G-1 to G-5. The diverse cultivars of P. ostreatus were discriminated and the other four species revealed a different section in the UPGMA tree. These microsatellite markers proved to be very useful tools for genetic studies, including assessment of the diversity and population structure of P. ostreatus.

Genomic Insights into the Fungal Lignocellulolytic Machinery of Flammulina rossica.

Next-generation sequencing (NGS) of the Flammulina rossica (wood-rotting basidiomycete) genome was performed to identify its carbohydrate-active enzymes (CAZymes). De novo genome assembly (31 kmer) revealed a total length of 35,646,506 bp (49.79% GC content). In total, 12,588 gene models of F. rossica were predicted using an ab initio gene prediction tool (AUGUSTUS). Orthologous analysis with other fungal species revealed that 7433 groups contained at least one F. rossica gene. Additionally, 12,033 (95.6%) of 12,588 genes for F. rossica proteins had orthologs among the Dikarya, and F. rossica contained 12 species-specific genes. CAZyme annotation in the F. rossica genome revealed 511 genes predicted to encode CAZymes including 102 auxiliary activities, 236 glycoside hydrolases, 94 glycosyltransferases, 19 polysaccharide lyases, 56 carbohydrate esterases, and 21 carbohydrate binding-modules. Among the 511 genes, several genes were predicted to simultaneously encode two different CAZymes such as glycoside hydrolases (GH) as well as carbohydrate-binding module (CBM). The genome information of F. rossica offers opportunities to understand the wood-degrading machinery of this fungus and will be useful for biotechnological and industrial applications.

Development of a Molecular Marker Linked to the A4 Locus and the Structure of HD Genes in Pleurotus eryngii.

Allelic differences in A and B mating-type loci are a prerequisite for the progression of mating in the genus Pleurotus eryngii; thus, the crossing is hampered by this biological barrier in inbreeding. Molecular markers linked to mating types of P. eryngii KNR2312 were investigated with randomly amplified polymorphic DNA to enhance crossing efficiency. An A4-linked sequence was identified and used to find the adjacent genomic region with the entire motif of the A locus from a contig sequenced by PacBio. The sequence-characterized amplified region marker 7-2299 distinguished A4 mating-type monokaryons from KNR2312 and other strains. A BLAST search of flanked sequences revealed that the A4 locus had a general feature consisting of the putative HD1 and HD2 genes. Both putative HD transcription factors contain a homeodomain sequence and a nuclear localization sequence; however, valid dimerization motifs were found only in the HD1 protein. The ACAAT motif, which was reported to have relevance to sex determination, was found in the intergenic region. The SCAR marker could be applicable in the classification of mating types in the P. eryngii breeding program, and the A4 locus could be the basis for a multi-allele detection marker.

Molecular Characterization of 170 New gDNA-SSR Markers for Genetic Diversity in Button Mushroom (Agaricus bisporus).

We designed 170 new simple sequence repeat (SSR) markers based on the whole-genome sequence data of button mushroom (Agaricus bisporus), and selected 121 polymorphic markers. A total of 121 polymorphic markers, the average major allele frequency (MAF) and the average number of alleles (NA) were 0.50 and 5.47, respectively. The average number of genotypes (NG), observed heterozygosity (HO), expected heterozygosity (HE), and polymorphic information content (PIC) were 6.177, 0.227, 0.619, and 0.569, respectively. Pearson's correlation coefficient showed that MAF was negatively correlated with NG (-0.683), NA (-0.600), HO (-0.584), and PIC (-0.941). NG, NA, HO, and PIC were positively correlated with other polymorphic parameters except for MAF. UPGMA clustering showed that 26 A. bisporus accessions were classified into 3 groups, and each accession was differentiated. The 121 SSR markers should facilitate the use of molecular markers in button mushroom breeding and genetic studies.

Characterization of Volatile Profiles of Six Popular Edible Mushrooms Using Headspace-Solid-Phase Microextraction Coupled with Gas Chromatography Combined with Chemometric Analysis.

The classification of six mushroom species (white beech, brown beech, button, oyster, king oyster, and enoki mushrooms) was successfully achieved using canonical discriminant analysis (CDA) on volatile metabolite data sets obtained by headspace-solid-phase microextraction gas chromatography (HS-SPME-GC). Twenty-seven major volatile compounds in six edible mushrooms were positively identified by HS-SPME-GC mass spectroscopy. The total volatile content was highest in brown beech mushroom (P < 0.05). Significant difference in volatile profile was observed between brown beach and white beech mushrooms. Button mushroom contained significantly higher contents of benzaldehyde and benzyl alcohol than the other mushrooms (P < 0.05). Oyster mushroom contained 1-octen-3-ol as the most prevalent volatile, representing 67% out of total volatiles. Hexanal (35.0%) and 1-octen-3-ol (22.5%) were the most abundant volatiles found in king oyster. Hexanal (29.1%) was the most prevalent volatile in enoki mushroom only. Several volatile pairs with very high positive correlation in their levels were identified, representing the highest correlation coefficient (r = 0.970) for the pair of t-2-octenal and 2,4-octandienal. CDA was much more efficient than principal component analysis for the differentiation of mushroom species. PRACTICAL APPLICATION: The present study provided the important information on the volatile metabolite profiles of popular six commercial mushroom species. The present data will be useful for the quality control of mushrooms cultivated in farms and mushroom products processed in food industry. The strategy of canonical discriminant analysis in combination with HS-SPME-GC could be expanded for the determining the authentication of mushroom species.