Pleurotus ostreatus is a commercially grown mushroom species in China. However, studies on the mechanisms of the fruiting body development and stress response of P. ostreatus are still at a primary stage. In this study, we report the entire genome sequence of P. ostreatus CCMSSC03989. Then, we performed comprehensive genome-wide characterization and expression analysis of the MYB transcription factor family during a series of developmental stages and under the condition of heat stress. A 34.76 Mb genome was obtained through next-generation sequencing (NGS) and Bionano optical mapping approaches. The genome has a scaffold N50 of 1.1 Mb and contains 10.11% repeats, and 10,936 gene models were predicted. A total of 20 MYB genes (PoMYB) were identified across the genome, and the full-length open reading frames were isolated. The PoMYBs were classified into 1 repeat (1R), 2R, and 3R-MYB groups according to their MYB domain repeat numbers, and 3R-MYBs possessed relatively more introns than 1R and 2R-MYBs. Based on phylogenetic analysis, the PoMYBs were divided into four groups and showed close relationships with the MYB genes of plants and fungi. RNA-sequencing (RNA-Seq) and quantitative PCR (qPCR) analyses revealed that PoMYB expression showed stage-specific patterns in reproductive stages and could be induced by heat stress. The P. ostreatus draft genome will promote genome-wide analysis, and our study of PoMYBs will promote further functional analysis of MYB genes in mushrooms.

Catalases are ubiquitous hydrogen peroxide-detoxifying enzymes. They participate in fungal growth and development, such as mycelial growth and cellular differentiation, and in protecting fungi from oxidative damage under stressful conditions. To investigate the potential functions of catalases in Pleurotus ostreatus, we obtained two catalase genes from a draft genome sequence of P. ostreatus, and cloned and characterized them (Po-cat1 and Po-cat2). Po-cat1 (group II) and Po-cat2 (group III) encoded putative peptides of 745 and 528 amino acids, respectively. Furthermore, the gene structures were variant between Po-cat1 and Po-cat2. Further research revealed that these two catalase genes have divergent expression patterns during different developmental stages. Po-cat1/Po-cat1 was at a barely detectable level in mycelia, accumulated gradually during reproductive growth, and was maximal in separated spores. But no catalase activity of Po-cat1 was detected by native-PAGE during any part of the developmental stages. In contrast, high Po-cat2/Po-cat2 expression and Po-cat2 activity found in mycelia were gradually lost during reproductive growth, and at a minimal level in separated spores. In addition, these two genes responded differentially under 32 °C and 40 °C heat stresses. Po-cat1 was up-regulated under both temperature conditions, while Po-cat2 was up-regulated at 32 °C but down-regulated at 40 °C. The accumulation of catalase proteins correlated with gene expression. These results indicate that the two divergent catalases in P. ostreatus may play different roles during development and under heat stress.

Pleurotus tuoliensis (Bailinggu) is a commercially cultivated mushroom species with an increasing popularity in China and other Asian countries. Commercial profits are now low, mainly due to a low yield, long cultivation period and sensitivity to diseases. Breeding efforts are thus required to improve agronomical important traits. Developing saturated genetic linkage and physical maps is a start for applying genetic and molecular approaches to accelerate the precise breeding programs. Here we present a genetic linkage map for P. tuoliensis constructed by using 115 haploid monokaryons derived from a hybrid strain H6. One thousand one hundred and eighty-two SNP markers developed by 2b–RAD (type IIB restriction-site associated DNA) approach were mapped to 12 linkage groups. The map covers 1073 cM with an average marker spacing of 1.0 cM. The genome of P. tuoliensis was de novo sequenced as 40.8 Mb and consisted of 500 scaffolds (>500 bp), which showed a high level of colinearity to the genome of P. eryngii var. eryngii. A total of 97.4% SNP markers (1151) were physically localized on 78 scaffolds, and the physical length of these anchored scaffolds were 33.9 Mb representing 83.1% of the whole genome. Mating type loci A and B were mapped on separate linkage groups and identified physically on the assembled genomes. Five putative pheromone receptors and two putative pheromone precursors were identified for the mating type B locus. This study reported a first genetic linkage map integrated with physical mapping of the de novo sequenced genome and the mating type loci of an important cultivated mushroom in China, P. tuoliensis. The de novo sequenced and annotated genome, assembled using a 2b–RAD generated linkage map, provides a basis for marker-assisted breeding of this economic important mushroom species.

Oyster mushrooms are widely cultivated and consumed over the world for their easy cultivation and high biological efficiency (mushroom fresh weight/substrate dry weight × 100%). Fruiting bodies with dark caps are more and more popular according to consumer preferences, but dark varieties are rarely seen on the market.