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 δ¹⁵N 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.

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, δ¹⁵N 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.

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.

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⁻¹ 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.

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.

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.

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 (M_AF) and the average number of alleles (N_A) were 0.50 and 5.47, respectively. The average number of genotypes (N_G), observed heterozygosity (H_O), expected heterozygosity (H_E), and polymorphic information content (PIC) were 6.177, 0.227, 0.619, and 0.569, respectively. Pearson's correlation coefficient showed that M_AF was negatively correlated with N_G (-0.683), N_A (-0.600), H_O (-0.584), and PIC (-0.941). N_G, N_A, H_O, and PIC were positively correlated with other polymorphic parameters except for M_AF. 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.

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-2₂₉₉ 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.

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.

Development of Polymorphic Simple Sequence Repeat Markers using High-Throughput Sequencing in Button Mushroom (Agaricus bisporus)

The white button mushroom (Agaricus bisporus) is one of the most widely cultivated species of edible mushroom. Despite its economic importance, relatively little is known about the genetic diversity of this species. Illumina paired-end sequencing produced 43,871,558 clean reads and 69,174 contigs were generated from five offspring. These contigs were subsequently assembled into 57,594 unigenes. The unigenes were annotated with reference genome in which 6,559 unigenes were associated with clusters, indicating orthologous genes. Gene ontology classification assigned many unigenes. Based on genome data of the five offspring, 44 polymorphic simple sequence repeat (SSR) markers were developed. The major allele frequency ranged from 0.42 to 0.92. The number of genotypes and the number of alleles ranged from 1 to 4, and from 2 to 4, respectively. The observed heterozygosity and the expected heterozygosity ranged from 0.00 to 1.00, and from 0.15 to 0.64, respectively. The polymorphic information content value ranged from 0.14 to 0.57. The genetic distances and UPGMA clustering discriminated offspring strains. The SSR markers developed in this study can be applied in polymorphism analyses of button mushroom and for cultivar discrimination.

Genome-Wide Comparison of Carbohydrate-Active Enzymes (CAZymes) Repertoire of Flammulina ononidis

Whole-genome sequencing of Flammulina ononidis, a wood-rotting basidiomycete, was performed to identify genes associated with carbohydrate-active enzymes (CAZymes). A total of 12,586 gene structures with an average length of 2009 bp were predicted by the AUGUSTUS tool from a total 35,524,258 bp length of de novo genome assembly (49.76% GC). Orthologous analysis with other fungal species revealed that 7051 groups contained at least one F. ononidis gene. In addition, 11,252 (89.5%) of 12,586 genes for F. ononidis proteins had orthologs among the Dikarya, and F. ononidis contained 8 species-specific genes, of which 5 genes were paralogous. CAZyme prediction revealed 524 CAZyme genes, including 228 for glycoside hydrolases, 21 for polysaccharide lyases, 87 for glycosyltransferases, 61 for carbohydrate esterases, 87 with auxiliary activities, and 40 for carbohydrate-binding modules in the F. ononidis genome. This genome information including CAZyme repertoire will be useful to understand lignocellulolytic machinery of this white rot fungus F. ononidis.

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5159-y) contains supplementary material, which is available to authorized users.

[Clinical analysis of chronic sternal osteomyelitis with sinus tract after cardiovascular surgery]

Objective: To explore the diagnosis, therapy and prevention method of chronic sternal osteomyelitis with sinus tract after cardiovascular surgery. Methods: A total of 53 patients with chronic sternal sinus tract after cardiovascular surgery between January 2000 and January 2016. After definite diagnosis by contrast fistulography and CT scanning, all the patients received combined modality therapy including debridement, musculocutaneous flap transplantation and intermediate thickness free skin graft transplantation if necessary. Results: One patient died of false aneurysm due to the sternal sinus tract infection, there were no peri-operative death for all the left 52 patients. Forty-five patients had primary healing and 7 patients had secondary healing. All the patients became total recovery within 3-12 weeks after operation and maintained well during the 5-18 months' follow-up. Conclusions: For the patients with chronic sternal osteomyelitis, operative therapy should be performed as soon as possible once the diagnosis is confirmed. Combined modality therapy including debridement, musculocutaneous flap transplantation and intermediate thickness free skin graft transplantation is confirmed to be effective and secure.

Anti-inflammatory effect of lovastatin is mediated via the modulation of NF-κB and inhibition of HDAC1 and the PI3K/Akt/mTOR pathway in RAW264.7 macrophages

Lovastatin is a 3-hydroxy-3-methylglutaryl-CoA reductase inhibitor that is clinically used for the prevention of cardiovascular diseases. Although it has been reported that lovastatin has anti-inflammatory properties in several studies, how lovastatin regulates the inflammation is still unclear. To evaluate the effect of lovastatin on nitric oxide production (NO) in RAW264.7 macrophages, NO production assay was performed. Also, cell viability was measured to confirm cytotoxicity. Level of tumor necrosis factor-α (TNF-α) transcription was measured by reverse transcription polymerase chain reaction (RT-PCR) from total RNA in RAW264.7 cells. Western blot analysis and immunofluorescence staining were used to investigate the regulation of lovastatin on the expression, phosphorylation, and nuclear translocation of cellular proteins. The results of the present study revealed that lovastatin reduced nitric oxide production via the reduction of inducible nitric oxide synthase (iNOS) expression in lipopolysaccharide (LPS)-stimulated RAW264.7 macrophage cells. The mRNA level of TNF-α was reduced in presence of lovastatin. In addition, lovastatin downregulated histone deacetylase 1 (HDAC1), resulting in the accumulation of acetylated histone H3 and heat shock protein 70. Furthermore, the expression of phosphoinositide 3-kinase catalytic subunits α and β was reduced under lovastatin treatment, and the phosphorylation of Akt and mammalian target of rapamycin was consequently inhibited. Lovastatin also inhibited the phosphorylation of inhibitor of nuclear factor (NF)-κBα and the translocation of NF-κB into the nucleus. Therefore, the present study demonstrates that lovastatin inhibits the expression of pro-inflammatory mediators, including iNOS and TNF-α, through the suppression of HDAC1 expression, PI3K/Akt phosphorylation and NF-κB translocation in LPS-stimulated RAW264.7 macrophage cells.

Ribosomal Intergenic Spacer 1 Based Characterization of Button Mushroom (Agaricus bisporus) Strains

Breeding the button mushroom requires genetic information about its strains. This study was undertaken to genetically characterize four domestically bred button mushroom strains (Saea, Saejung, Saedo, Saeyeon cultivars) and to assess the possibility of using the intergenic spacer 1 (IGS1) region of rDNA as a genetically variable region in the genetic characterization. For the experiment, 34 strains of Agaricus bisporus, two strains of A. bitorquis, and one strain of A. silvaticus, from 17 countries were used. Nucleotide sequence analysis of IGS1 rDNA in these 37 Agaricus strains confirmed that genetic variations exist, not only among the four domestic strains, but also between the four domestic strains and foreign strains. Crossing two different haploid strains of A. bisporus seems to generate genetic variation in the IGS1 region in their off-spring haploid strains. Phylogenetic analysis based on the IGS1 sequence revealed all A. bisporus strains could be differentiated from A. silvaticus and A. bitorquis strains. Five genetic groups were resolved among A. bisporus strains. Saejung and Saeyeon cultivars formed a separate genetic group. Our results suggest that IGS1 could be complementarily applied in the polymorphism analysis of button mushroom.

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.

Taxonomic Position and Species Identity of the Cultivated Yeongji 'Ganoderma lucidum' in Korea

Ganoderma lucidum has a long history of use as a traditional medicine in Asian countries. However, the taxonomy of Ganoderma species remains controversial, since they were initially classified on the basis of their morphological characteristics. Recently, it was proposed that G. lucidum from China be renamed as G. sichuanense or G. lingzhi. In the present study, phylogenetic analysis using the internal transcribed spacer region rDNA sequences of the Ganoderma species indicated that all strains of the Korean 'G. lucidum' clustered into one group together with G. sichuanense and G. lingzhi from China. However, strains from Europe and North American, which were regarded as true G. lucidum, were positioned in a clearly different group. In addition, the average size of the basidiospores from the Korean cultivated Yeongji strains was similar to that of G. lingzhi. Based on these results, we propose that the Korean cultivated Yeongji strains of 'G. lucidum' should be renamed as G. lingzhi.

Genetic and Biochemical Characterization of Monokaryotic Progeny Strains of Button Mushroom (Agaricus bisporus)

To promote the selection of promising monokaryotic strains of button mushroom (Agaricus bisporus) during breeding, 61 progeny strains derived from basidiospores of two different lines of dikaryotic parental strains, ASI1038 and ASI1346, were analyzed by nucleotide sequencing of the intergenic spacer I (IGS I) region in their rDNA and by extracellular enzyme assays. Nineteen different sizes of IGS I, which ranged from 1,301 to 1,348 bp, were present among twenty ASI1346-derived progeny strains, while 15 different sizes of IGS I, which ranged from 700 to 1,347 bp, were present among twenty ASI1038-derived progeny strains. Phylogenetic analysis of the IGS sequences revealed that different clades were present in both the ASI10388- and ASI1346-derived progeny strains. Plating assays of seven kinds of extracellular enzymes (β-glucosidase, avicelase, CM-cellulase, amylase, pectinase, xylanase, and protease) also revealed apparent variation in the ability to produce extracellular enzymes among the 40 tested progeny strains from both parental A. bisporus strains. Overall, this study demonstrates that characterization of IGS I regions and extracellular enzymes is useful for the assessment of the substrate-degrading ability and heterogenicity of A. bisporus monokaryotic strains.

Genome-Wide Identification and Characterization of Novel Laccase Genes in the White-Rot Fungus Flammulina velutipes

The aim of this study was to identify and characterize new Flammulina velutipes laccases from its whole-genome sequence. Of the 15 putative laccase genes detected in the F. velutipes genome, four new laccase genes (fvLac-1, fvLac-2, fvLac3, and fvLac-4) were found to contain four complete copper-binding regions (ten histidine residues and one cysteine residue) and four cysteine residues involved in forming disulfide bridges, fvLac-1, fvLac-2, fvLac3, and fvLac-4, encoding proteins consisting of 516, 518, 515, and 533 amino acid residues, respectively. Potential N-glycosylation sites (Asn-Xaa-Ser/Thr) were identified in the cDNA sequence of fvLac-1 (Asn-454), fvLac-2 (Asn-437 and Asn-455), fvLac-3 (Asn-111 and Asn-237), and fvLac4 (Asn-402 and Asn-457). In addition, the first 19~20 amino acid residues of these proteins were predicted to comprise signal peptides. Laccase activity assays and reverse transcription polymerase chain reaction analyses clearly reveal that CuSO4 affects the induction and the transcription level of these laccase genes.

Identification of degenerate nuclei and development of a SCAR marker for Flammulina velutipes

Flammulina velutipes is one of the major edible mushrooms in the world. Recently, abnormalities that have a negative impact on crop production have been reported in this mushroom. These symptoms include slow vegetative growth, a compact mycelial mat, and few or even no fruiting bodies. The morphologies and fruiting capabilities of monokaryons of wild-type and degenerate strains that arose through arthrospore formation were investigated through test crossing. Only one monokaryotic group of the degenerate strains and its hybrid strains showed abnormal phenotypes. Because the monokaryotic arthrospore has the same nucleus as the parent strain, these results indicated that only one aberrant nucleus of the two nuclei in the degenerate strain was responsible for the degeneracy. A sequence-characterized amplified region marker that is linked to the degenerate monokaryon was identified based on a polymorphic sequence that was generated using random primers. Comparative analyses revealed the presence of a degenerate-specific genomic region in a telomere, which arose via the transfer of a genomic fragment harboring a putative helicase gene. Our findings have narrowed down the potential molecular targets responsible for this phenotype for future studies and have provided a marker for the detection of degenerate strains.

Complete mitochondrial genome of the ectomycorrhizal fungus Tricholoma matsutake

Tricholoma matsutake is an ectomycorrhizal fungus belonging to the class Agaricomycetes and the phylum Basidiomycota. Here, we report the complete mitochondrial genome sequence of T. matsutake. The mitochondrial DNA (mtDNA) of T. matsutake seems to be concatenated, and harbors 76,037 nucleotides, with an overall GC content of 20.69%. The mtDNA comprises 15 protein-coding genes, 28 tRNA genes, and 2 rRNA genes. A comparative mitogenomic analysis between Agaricomycetes and other classes (Exobasidiomycetes, Pucciniomycetes, Tremellomycetes, and Ustilaginomycetes) revealed that species belonging to the class Agaricomycetes preferably harbor the codon UUA (coding for leucine) in their mtDNAs, whereas other classes prefer the codon CUA. This bias in the usage of the leucine codon could be the result of evolutionary divergence between the Agaricomycetes and other classes. The T. matsutake mtDNA sequence thus provides insight into the evolutionary characteristics of the Agaricomycetes mtDNAs.

Identification and functional analysis of pheromone and receptor genes in the B3 mating locus of Pleurotus eryngii

Pleurotus eryngii has recently become a major cultivated mushroom; it uses tetrapolar heterothallism as a part of its reproductive process. Sexual development progresses only when the A and B mating types are compatible. Such mating incompatibility occasionally limits the efficiency of breeding programs in which crossing within loci-shared strains or backcrossing strategies are employed. Therefore, understanding the mating system in edible mushroom fungi will help provide a short cut in the development of new strains. We isolated and identified pheromone and receptor genes in the B3 locus of P. eryngii and performed a functional analysis of the genes in the mating process by transformation. A genomic DNA library was constructed to map the entire mating-type locus. The B3 locus was found to contain four pheromone precursor genes and four receptor genes. Remarkably, receptor PESTE3.3.1 has just 34 amino acid residues in its C-terminal cytoplasmic region; therefore, it seems likely to be a receptor-like gene. Real-time quantitative RT-PCR (real-time qRT-PCR) revealed that most pheromone and receptor genes showed significantly higher expression in monokaryotic cells than dikaryotic cells. The pheromone genes PEphb3.1 and PEphb3.3 and the receptor gene PESTE3.3.1 were transformed into P5 (A3B4). The transformants were mated with a tester strain (A4B4), and the progeny showed clamp connections and a normal fruiting body, which indicates the proposed role of these genes in mating and fruiting processes. This result also confirms that PESTE3.3.1 is a receptor gene. In this study, we identified pheromone and receptor genes in the B3 locus of P. eryngii and found that some of those genes appear to play a role in the mating and fruiting processes. These results might help elucidate the mechanism of fruiting differentiation and improve breeding efficiency.

Overproduction of Laccase by the White-Rot Fungus Pleurotus ostreatus Using Apple Pomace as Inducer

Laccase activity of Pleurotus ostreatus is significantly increased by the addition of apple pomace. Among various conditions, the best concentration of apple pomace and cultivation time for the production of laccase by P. ostreatus was 2.5% and 9 days, respectively. Reverse transcription polymerase chain reaction analyses of laccase isoenzyme genes, including pox1, pox3, pox4, poxc, poxa3, and poxa1b, revealed a clear effect of apple pomace on transcription induction. Our findings reveal that the use of apple pomace can be a model for the valuable addition of similar wastes and for the development of a solid-state fermenter and commercial production of oyster mushroom P. ostreatus.

Cloning and mRNA expression analysis of the gene encoding phenylalanine ammonia-lyase of the ectomycorrhizal fungus Tricholoma matsutake

The ectomycorrhizal fungus Tricholoma matsutake grows symbiotically with Pinus densiflora. Phenylalanine ammonia-lyase (E.C. 4.3.1.24) catalyzes the conversion of L-phenylalanine to trans-cinnamic acid. The role of fungal phenylalanine ammonia-lyase, however, has not been clear until now. In this study, the gene encoding phenylalanine ammonia-lyase (PAL), which was isolated from T. matsutake, was cloned and characterized. The PAL gene (tmpal) consists of 2,160 nucleotides, coding for a polypeptide containing 719 amino acid residues. The deduced amino acid sequence of tmpal from T. matsutake shows high identity (70%) with that from Laccaria bicolor. Comparative analysis of the PAL genes among T. matsutake and other species of the class Agaricomycetes showed that both active sites and binding sites were significantly conserved among these genes. The transcriptional analysis of the PAL gene revealed a differential gene expression pattern depending on the developmental stages (mycelium, primordium, stipe, pileus, and gills) of T. matsutake. These results suggest that the PAL gene in T. matsutake plays an important role in multiple physiological functions.

Complete mitochondrial genome of the entomopathogenic fungus Beauveria pseudobassiana (Ascomycota, Cordycipitaceae)

The complete mitochondrial genome of the entomopathogenic fungus Beauveria pseudobassiana (Ascomycota, Cordycipitaceae) was determined and found to be 28,006 bp in length. It encodes genes for 14 proteins, 2 ribosomal RNA subunits, 25 transfer RNAs with the synteny identical to those of B. bassiana and B. brongniartii. The overall base composition is 36.7% A, 35.8% T, 15.3% C and 12.2% G with an AT content of 72.5%. Group-І introns were found in large rRNA gene (rnl) and cox2 gene that include rps3 gene and putative GIY-YIG homing endonuclease, respectively. The mitochondrial genome of B. pseudobassiana will contribute to the investigation of the phylogenetic relationship, taxonomic resolution and biogeography of Beauveria.

Identification of mating type loci and development of SCAR marker genetically linked to the B3 locus in Pleurotus eryngii

In order to estimate how diverse the mating types in Pleurotus eryngii from different regions are, pairings between monokaryons derived from inter- and intragroups were done. Sixteen and 15 alleles were identified at loci A and B from the 12 strains. In the P. eryngii KNR2312, widely used for commercial production, four mating loci, A3, A4, B3, and B4, were determined. Those loci, except A3, were found in 4 strains out of 12 strains. To improve breeding efficiency, especially in mating type determination, RAPD and BSA were performed to screen for a mating type specific marker. The SCAR marker 13- 2(2100) was developed based on the RAPD-derived sequence typing B3 locus. The sequence analysis of 13-2(2100) revealed that it contained a conserved domain, the STE3 superfamily, and consensus sequences like the TATA box and GC box. It seems likely that the SCAR marker region is a part of the pheromone receptor gene.

Comparative genomics of the mating-type loci of the mushroom Flammulina velutipes reveals widespread synteny and recent inversions

BACKGROUND

Mating-type loci of mushroom fungi contain master regulatory genes that control recognition between compatible nuclei, maintenance of compatible nuclei as heterokaryons, and fruiting body development. Regions near mating-type loci in fungi often show adapted recombination, facilitating the generation of novel mating types and reducing the production of self-compatible mating types. Compared to other fungi, mushroom fungi have complex mating-type systems, showing both loci with redundant function (subloci) and subloci with many alleles. The genomic organization of mating-type loci has been solved in very few mushroom species, which complicates proper interpretation of mating-type evolution and use of those genes in breeding programs.

METHODOLOGY/PRINCIPAL FINDINGS

We report a complete genetic structure of the mating-type loci from the tetrapolar, edible mushroom Flammulina velutipes mating type A3B3. Two matB3 subloci, matB3a that contains a unique pheromone and matB3b, were mapped 177 Kb apart on scaffold 1. The matA locus of F. velutipes contains three homeodomain genes distributed over 73 Kb distant matA3a and matA3b subloci. The conserved matA region in Agaricales approaches 350 Kb and contains conserved recombination hotspots showing major rearrangements in F. velutipes and Schizophyllum commune. Important evolutionary differences were indicated; separation of the matA subloci in F. velutipes was diverged from the Coprinopsis cinerea arrangement via two large inversions whereas separation in S. commune emerged through transposition of gene clusters.

CONCLUSIONS/SIGNIFICANCE

In our study we determined that the Agaricales have very large scale synteny at matA (∼350 Kb) and that this synteny is maintained even when parts of this region are separated through chromosomal rearrangements. Four conserved recombination hotspots allow reshuffling of large fragments of this region. Next to this, it was revealed that large distance subloci can exist in matB as well. Finally, the genes that were linked to specific mating types will serve as molecular markers in breeding.

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.

Electrophoretic karyotyping and construction of a bacterial artificial chromosome library of the winter mushroom Flammulina velutipes

An electrophoretic karyotype of Korean Flammulina velutipes was obtained using contour-clamped homogeneous electric field (CHEF) gel electrophoresis. The monokaryotic progenies (4019-20 and 4019-18) and a dikaryotic strain (4019-2018) had 6-8 chromosomes, with sizes ranging from 1.6 to 5.8 megabase (Mb) pairs. Among the 3 strains that were examined, strain 4019-20 resolved into at least 7 chromosomes, with a total genome size of approximately 26.7Mb. We selected several chromosome-specific genes from the cDNA library of F. velutipes using Southern hybridization analysis. In order to determine the whole genome sequence, we constructed a bacterial artificial chromosome (BAC) library of the 4019-20 strain. The BAC library comprised 3840 clones, and the insert size ranged from 60 to 228kb, with an average size of 136kb.

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.

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.

Screening of Biodegradable Function of Indigenous Ligno-degrading Mushroom Using Dyes

The process of biodegradation in lingo-cellulosic materials is critically relevant to biospheric carbon. The study of this natural process has largely involved laboratory investigations, focused primarily on the biodegradation and recycling of agricultural by-products, generally using basidiomycetes species. In order to collect super white rot fungi and evaluate its ability to degrade lingo-cellulosic material, 35 fungal strains, collected from forests, humus soil, livestock manure, and dead trees, were screened for enzyme activities and their potential to decolorize the commercially used Poly-R 478 dye. In the laccase enzymatic analysis chemical test, 33 white rot fungi and 2 brown rot fungi were identified. The degradation ability of polycyclic aromatic hydrocarbons (PAHs) according to the utilized environmental conditions was higher in the mushrooms grown in dead trees and fallen leaves than in the mushrooms grown in humus soil and livestock manure. Using Poly-R 478 dye to assess the PAH-degradation activity of the identified strains, four strains, including Agrocybe pediades, were selected. The activities of laccase, MnP, and Lip of the four strains with PAH-degrading ability were highest in Pleurotus incarnates. 87 fungal strains, collected from forests, humus soil, livestock manure, and dead trees, were screened for enzyme activities and their potential to decolorize the commercially used Poly-R 478 dye on solid media. Using Poly-R 478 dye to assess the PAHdegrading activity of the identified strains, it was determined that MKACC 51632 and 52492 strains evidenced superior activity in static and shaken liquid cultures. Subsequent screening on plates containing the polymeric dye poly R-478, the decolorization of which is correlated with lignin degradation, resulted in the selection of a strain of Coriolus versicolor, MKACC52492, for further study, primarily due to its rapid growth rate and profound ability to decolorize poly R-478 on solid media. Considering our findings using Poly-R 478 dye to evaluate the PAH-degrading activity of the identified strains, Coriolus versicolor, MKACC 52492 was selected as a favorable strain. Coriolus versicolor, which was collected from Mt. Yeogi in Suwon, was studied for the production of the lignin-modifying enzymes laccase, manganese-dependent peroxidase (MnP), and lignin peroxidase (LiP).

Plant growth promotion and Penicillium citrinum

BACKGROUND

Endophytic fungi are known plant symbionts. They produce a variety of beneficial metabolites for plant growth and survival, as well as defend their hosts from attack of certain pathogens. Coastal dunes are nutrient deficient and offer harsh, saline environment for the existing flora and fauna. Endophytic fungi may play an important role in plant survival by enhancing nutrient uptake and producing growth-promoting metabolites such as gibberellins and auxins. We screened roots of Ixeris repenes (L.) A. Gray, a common dune plant, for the isolation of gibberellin secreting endophytic fungi.

RESULTS

We isolated 15 endophytic fungi from the roots of Ixeris repenes and screened them for growth promoting secondary metabolites. The fungal isolate IR-3-3 gave maximum plant growth when applied to waito-c rice and Atriplex gemelinii seedlings. Analysis of the culture filtrate of IR-3-3 showed the presence of physiologically active gibberellins, GA1, GA3, GA4 and GA7 (1.95 ng/ml, 3.83 ng/ml, 6.03 ng/ml and 2.35 ng/ml, respectively) along with other physiologically inactive GA5, GA9, GA12, GA15, GA19, GA20 and, GA24. The plant growth promotion and gibberellin producing capacity of IR-3-3 was much higher than the wild type Gibberella fujikuroi, which was taken as control during present study. GA5, a precursor of bioactive GA3 was reported for the first time in fungi. The fungal isolate IR-3-3 was identified as a new strain of Penicillium citrinum (named as P. citrinum KACC43900) through phylogenetic analysis of 18S rDNA sequence.

CONCLUSION

Isolation of new strain of Penicillium citrinum from the sand dune flora is interesting as information on the presence of Pencillium species in coastal sand dunes is limited. The plant growth promoting ability of this fungal strain may help in conservation and revegetation of the rapidly eroding sand dune flora. Penicillium citrinum is already known for producing mycotoxin citrinin and cellulose digesting enzymes like cellulase and endoglucanase, as well as xylulase. Gibberellins producing ability of this fungus and the discovery about the presence of GA5 will open new aspects of research and investigations.

Cloning and developmental expression of a metzincin family metalloprotease cDNA from oyster mushroom Pleurotus ostreatus

A cDNA clone, PoMTP, encoding a putative metzincin family metalloprotease was isolated from the expressed sequence tags of a basidiomycete Pleurotus ostreatus. The 5'-end sequence of PoMTP was determined by the 5'-RACE method. Full-length cDNA sequence (1140 bp) of PoMTP contained a 870 bp open reading frame encoding a protein product of 290 amino acids in addition to a 99 bp of 5'-untranslated sequence and a 171 bp of 3'-untranslated sequence with a poly(A) tail. The deduced amino-acid sequences of PoMTP contained an extensive zinc-binding consensus sequence and a so-called Met-turn sequence which are typical for the metzincin family of metalloproteases, indicating that the PoMTP protein belongs to the metzincin metalloproteases. Four cysteine residues were also observed in the zinc-binding region of PoMTP amino-acid sequence, which are known to be important for the structure and the function of some subfamilies of the metzincins. Comparison of the PoMTP in sequence database showed no significant homology with functionally known metalloproteases of Armillaria mellea, Grifola frondosa, Lentinula edodes, Pleurotus ostreatus, Schizophyllum commune and Tricholoma saponaceum in mushroom. Northern blot and qunatitative RT-PCR analyses indicated the PoMTP mRNA to be abundant at primordial and fruit body stages, but scarce at the mycelial stage, suggesting that the PoMTP metalloprotease plays an important role in mushroom fruiting.