Selection and validation of internal control genes for quantitative real-time RT‒qPCR normalization of Phlebopus portentosus gene expression under different conditions.

Phlebopus portentosus (Berk. and Broome) Boedijn is an attractive edible mushroom and is considered the only bolete for which artificial cultivation in vitro has been achieved. Gene expression analysis has become widely used in research on edible fungi and is important for elucidating the functions of genes involved in complex biological processes. Selecting appropriate reference genes is crucial to ensuring reliable RT‒qPCR gene expression analysis results. In our study, a total of 12 candidate control genes were selected from 25 traditional housekeeping genes based on their expression stability in 9 transcriptomes of 3 developmental stages. These genes were further evaluated using geNorm, NormFinder, and RefFinder under different conditions and developmental stages. The results revealed that MSF1 domain-containing protein (MSF1), synaptobrevin (SYB), mitogen-activated protein kinase genes (MAPK), TATA-binding protein 1 (TBP1), and SPRY domain protein (SPRY) were the most stable reference genes in all sample treatments, while elongation factor 1-alpha (EF1), actin and ubiquitin-conjugating enzyme (UBCE) were the most unstably expressed. The gene SYB was selected based on the transcriptome results and was identified as a novel reference gene in P. portentosus. This is the first detailed study on the identification of reference genes in this fungus and may provide new insights into selecting genes and quantifying gene expression.

LC/MS- and GC/MS-based metabolomic profiling to determine changes in flavor quality and bioactive components of Phlebopus portentosus under low-temperature storage.

Introduction: Low temperature is the most common method used to maintain the freshness of Phlebopus portentosus during long-distance transportation. However, there is no information regarding the nutritional changes that occur in P. portentosus preserved postharvest in low temperature. Methods: In this study, the changes in flavor quality and bioactive components in fruiting bodies stored at 4 °C for different storage periods were determined through LC/MS and GC/MS analyses. Sampling was performed at 0, 3, 5, 7, and 13 days storage. Results and Discussion: Based on the results, the metabolites present in caps and stipes were different at the same period and significantly different after 7 days of storage. A total of 583 and 500 different metabolites were detected in caps and stipes, respectively, and were mainly lipids and lipid-like molecules, organic acids and derivatives, organic oxygen compounds and others. Except for prenol lipids and nucleotides, the expression levels of most metabolites increased with longer storage time. In addition, geosmin was identified as the major contributor to earthy-musty odors, and the level of geosmin was increased when the storage time was short. Conclusion: The variations in these metabolites might cause changes in flavor quality and bioactive components in P. portentosus. Variations in these metabolites were thoroughly analyzed, and the results revealed how storage processes affect the postharvest quality of P. portentosus for the first time.

Two new species of Fistulina (Agaricales, Basidiomycota) from the Northern Hemisphere.

Phylogenetic and morphological analyses on samples of Fistulina from East Asia and North America were carried out, and two new species were described, namely, Fistulina americana and Fistulina orientalis, both previously known as Fistulina hepatica. The former is characterized by lateral stipitate basidiocarps, relatively small pores (7-8 per mm), a monomitic hyphal system with both clamp connections and simple septa, and ellipsoid basidiospores of 4-4.8 × 3-3.3 µm, and the species has been found on Quercus in North-East USA. F. orientalis is characterized by lateral stipitate basidiocarps, very small pores (11-12 per mm) with pruinose dissepiments, a monomitic hyphal system with both clamp connections and simple septa, and ovoid to subglobose basidiospores of 3-4 × 2.7-3 µm, and the species has been found on Castanopsis in East Asia. Phylogenetically, samples of F. americana and F. orientalis form two new lineages nested in the Fistulina clade.

The complete mitochondrial genome of the edible mushroom Grifola frondosa.

The culinary-medicinal mushroom Grifola frondosa is widely cultivated in East Asia. In this study, the complete mitochondrial genome of G. frondosa was determined using Illumina sequencing. The circular molecule was 197,486 bp in length with a content of 25.01% GC, which was one of the largest mitochondrial genomes in the order Polyporales. A total of 39 known genes encoding 13 common mitochondrial genes, 24 tRNA genes, 1 ribosomal protein s3 gene (rps3), and 1 DNA polymerase gene (dpo) were predicted in this genome. The phylogenetic analysis showed that G. frondosa clustered together with Sparassis crispa, Laetiporus sulphureus, Wolfiporia cocos, and Taiwanofungus camphoratus. The complete mitochondrial genome reported here may provide new insight into genetic information and evolution for further studies.

Whole-Genome and Transcriptome Sequencing of Phlebopus portentosus Reveals Its Associated Ectomycorrhizal Niche and Conserved Pathways Involved in Fruiting Body Development.

Phlebopus portentosus (Berk. and Broome) Boedijin, a widely consumed mushroom in China and Thailand, is the first species in the order Boletaceae to have been industrially cultivated on a large scale. However, to date, the lignocellulose degradation system and molecular basis of fruiting body development in P. portentosus have remained cryptic. In the present study, genome and transcriptome sequencing of P. portentosus was performed during the mycelium (S), primordium (P), and fruiting body (F) stages. A genome of 32.74 Mb with a 48.92% GC content across 62 scaffolds was obtained. A total of 9,464 putative genes were predicted from the genome, of which the number of genes related to plant cell wall-degrading enzymes was much lower than that of some saprophytic mushrooms with specific ectomycorrhizal niches. Principal component analysis of RNA-Seq data revealed that the gene expression profiles at all three stages were different. The low expression of plant cell wall-degrading genes also confirmed the limited ability to degrade lignocellulose. The expression profiles also revealed that some conserved and specific pathways were enriched in the different developmental stages of P. portentosus. Starch and sucrose metabolic pathways were enriched in the mycelium stage, while DNA replication, the proteasome and MAPK signaling pathways may be associated with maturation. These results provide a new perspective for understanding the key pathways and hub genes involved in P. portentosus development.

Deep Sequencing, Nested PCR, and Denaturing Gradient Gel Electrophoresis Reveal a Wider Distribution of Chinese Caterpillar Mushroom, Ophiocordyceps sinensis (Ascomycetes), in Native Soil Types.

Ophiocordyceps sinensis appears as stroma emerging from underground sclerotium enclosed by the skeleton of Thitarodes moth larvae. However, the actual distribution of the fungus in soil still remains unclarified. In this study, 40 soil samples were used for detection of O. sinensis to confirm its distribution in native habitats using denaturing gradient gel electrophoresis, nested internal transcribed spacer (ITS) PCR, and 454 pyrosequencing methods. The soil samples included six types: Os, where both stromata and host moth larvae were found; NL, representing no signs of stromata, but where moth larvae were found; NOs, where neither stroma nor moth larvae were found; BS, with bare soil without the presence of stroma of O. sinensis or moth larvae; AF, from soil surrounding the stroma; and MP, soil particles firmly wrapping the sclerotium of O. sinensis. Of 40 samples tested, 36 showed positive detection of O. sinensis by at least one of the three detection methods, with positive detection in all six sample types at all five sites. The results showed that traces of O. sinensis can be detected in locations with no macroscopically visible evidence of the fungus or its host and at least 100 m away from such locations.

Diverse function and regulation of CmSnf1 in entomopathogenic fungus Cordyceps militaris.

SNF1/AMPK protein kinases play important roles in fungal development and activation of catabolite-repressed genes. In this study, we characterized the role of SNF1 ortholog in Cordyceps militaris (CmSnf1). The vegetative growth of a CmSnf1 deletion mutant was (ΔCmSnf1) reduced by 42.2% with arabinose as a sole carbon source. Most strikingly, the ΔCmSnf1 produced only a few conidia and exhibited delayed conidial germination. We found that CmSnf1 was necessary for mycelium to penetrate the insect cuticle to form the fruiting body on silkworm pupae, consistent with the down-regulation of chitinase- and protease-encoding genes in ΔCmSnf1. However, cordycepin content increased by more than 7 times in culture supernatants. Correspondingly, the relative expression levels of cordycepin gene cluster members were also elevated. In particular, the expression of cns4 associated with cordycepin transfer was up-regulated >10-fold. Furthermore, transcriptional analysis showed that CmSnf1 regulated the expression of genes involved in cell autophagy and oxidative stress tolerance. We speculated that under environmental stress, both the ATG and SNF1 pathways might collaborate to sustain adverse environments. Our study provides an initial framework to probe the diverse function and regulation of CmSnf1 in C. militaris, which will shed more light on the direction of molecular improvement of medicinal fungi.

Bacterial Profiling and Dynamic Succession Analysis of Phlebopus portentosus Casing Soil Using MiSeq Sequencing.

Phlebopus portentosus (Berk. and Broome) Boedijin is a popular edible mushroom found in China and Thailand. To date, P. portentosus is the only species in the order Boletales that can be successfully cultivated worldwide. The use of a casing layer or casing soil overlaying the substrate is a crucial step in the production of this mushroom. In this study, bacterial profiling and dynamic succession analyses of casing soil during the cultivation of P. portentosus were performed. One hundred and fifty samples were collected, and MiSeq sequencing of the V3-V4 region of the 16S rRNA gene was conducted. After performing a decontamination procedure, only 38 samples were retained, including 6 casing soil-originated samples (OS), 6 casing soil samples (FHCS) and 5 upper substrate samples (FHCU) from the period of complete colonization by mycelia; 6 casing soil samples (PCS) and 5 upper substrate samples (PCU) from the primordium period; and 6 casing soil samples (FCS) and 4 upper substrate samples (FCU) from fruit body period. The results revealed that bacterial diversity increased sharply from the hyphal to the primordium stage and then decreased during harvesting. The non-metric multidimensional scaling (NMDS) ordination and analysis of similarities (ANOSIM) analysis suggested that the community composition during different stages was significantly different in casing soil. The most abundant phyla in all of the samples were Proteobacteria, Chloroflexi, Acidobacteria, Actinobacteria, Saccharibacteria, and Bacteroidetes. Burkholderia was the most abundant genus in all the samples except the OS samples. The relative abundance of Burkholderia in the FHCS samples (55.79%) decreased to 35.14% in the PCS samples and then increased to 45.60% in the FCS samples. The abundances of Acidobacterium, Rhizobium, Acidisphaera, Bradyrhizobium, and Bacillus increased from the FHCS to PCS samples. The linear discriminant analysis (LDA) effect size (LEfSe) suggested that Acidobacterium and Acidisphaera are micromarkers for PCS, whereas Bradyrhizobium, Roseiarcus, and Pseudolabrys were associated with fruit body stages. The network analyses resulted in 23 edges, including 4 negative and 19 positive edges. Extensive mutualistic interactions may occur among casing soil bacteria. Furthermore, these bacteria play important roles in mycelial elongation, primordium formations, and the production of increased yields.

Identification of a Novel Anti-cancer Protein, FIP-bbo, from Botryobasidium botryosum and Protein Structure Analysis using Molecular Dynamic Simulation.

Fungal immunoregulatory proteins (FIP) are effective small molecule proteins with broad-spectrum immunomodulatory and anti-cancer activities and can be potential agents for the development of clinical drugs and health food additives. In this study, a new member of FIP named FIP-bbo was obtained through Botryobasidium botryosum genome mining. FIP-bbo has the typical characteristics of FIP but is genetically distant from other FIPs. Recombinant FIP-bbo (rFIP-bbo) was produced in an optimized E. coli expression system, and the pure protein was isolated using a Ni-NTA column. Antineoplastic experiments suggested that FIP-bbo is similar to LZ-8 in inhibiting various cancer cells (Hela, Spac-1, and A549) at lower concentrations, but it is not as potent as LZ-8. The molecular mechanism by which FIP-bbo, FIP-fve, and LZ-8 are cytotoxic to cancer cells has been discussed based on molecular dynamics simulation. Point mutations that may improve the thermal stability of FIP-fve and FIP-bbo were predicted. These results not only present a new candidate protein for the development of anticancer adjuvants, but also provide an approach for designing FIPs with high anticancer activity.

The Reliability of DNA Sequences in Public Databases Belonging to the Most Economically Important Shiitake Culinary-Medicinal Mushroom Lentinus edodes (Agaricomycetes) in Asia.

Large numbers of DNA sequences deposited in the International Nucleotide Sequence Databases (INSD) are erroneously annotated. The erroneous information may lead to misleading conclusions or cause great economic losses to farmers. Lentinus edodes (= Lentinula edodes (Berk.) Pegler) is one of the most important and popular culinary-medicinal mushrooms with a high nutritional value. In this study, experimental and in silico methods were used to correct the sequences annotated as L. edodes in the INSD. A total of 3,426 nucleotide entries were retrieved from public databases, including 140 different types of genetic sequences. Excluding 1,893 genome sequences, the most abundant signatures represented by ITS (258) and IGS1 (259) sequences accounted for 33.23% of the total entries. A total of 3,058 sequences were annotated correctly, 350 were indeterminate, and 18 were annotated erroneously based on the two methods. Correction of sequences will be beneficial for species identification and annotation. Phylogenic analysis based on ITS sequences suggested that L. edodes segregate in four clades in the tree based on ITS sequences. The isolates from China were distributed into two clades. In L. edodes, the intraspecific variation of the ITS2 sequences was much higher than that of the ITS1 sequences. In addition, the genetic diversity of the L. edodes sequences from China was much higher than that of any other regions included in this study. The northwest and southwest regions of China were L. edodes diversity centers.

Evaluation of the ribosomal DNA internal transcribed spacer (ITS), specifically ITS1 and ITS2, for the analysis of fungal diversity by deep sequencing.

The nuclear ribosomal DNA internal transcribed spacer (ITS) has been widely used to assess the fungal composition in different environments by deep sequencing. To evaluate the ITS in the analysis of fungal diversity, comparisons of the clustering and taxonomy generated by sequencing with different portions of the whole fragment were conducted in this study. For a total of 83,120 full-length ITS sequences obtained from the UNITE database, it was found that, on average, ITS1 varied more than ITS2 within the kingdom Fungi; this variation included length and GC content variations and polymorphisms, with some polymorphisms specific to particular fungal groups. The taxonomic accuracy for ITS was higher than that for ITS1 or ITS2. The commonly used operational taxonomic unit (OTU) for evaluating fungal diversity and richness assigned several species to a single OTU even with clustering at 99.00% sequence similarity. The clustering and taxonomic capacities did not differ between ITS1 and ITS2. However, the OTU commonality between ITS1 and ITS2 was very low. To test this observation further, 219,741 pyrosequencing reads, including 39,840 full-length ITS sequences, were obtained from 10 soil samples and were clustered into OTUs. The pyrosequencing results agreed with the results of the in silico analysis. ITS1 might overestimate the fungal diversity and richness. Analyses using ITS, ITS1 and ITS2 yielded several different taxa, and the taxonomic preferences for ITS and ITS2 were similar. The results demonstrated that ITS2 alone might be a more suitable marker for revealing the operational taxonomic richness and taxonomy specifics of fungal communities when the full-length ITS is not available.

Computational Insights into the Molecular Mechanism of the High Immunomodulatory Activity of LZ-8 Protein Isolated from the Lingzhi or Reishi Medicinal Mushroom Ganoderma lucidum (Agaricomycetes).

Ling zhi-8 (LZ-8) is the first fungal immunomodulatory protein (FIP) isolated from the lingzhi or reishi medicinal mushroom, Ganoderma lucidum. LZ-8 effectively induces interleukin 2 expression and secretion by forming a stable homodimer, and it is regarded as a good candidate to become a new therapeutic agent and/or functional food supplement. However, the molecular mechanism by which LZ-8 dimerization influences the regulation of interleukin 2 is not clear. In this study we performed structure-based multiple alignment of LZ-8 and an FIP from Volvariella volvacea, compared the electrostatic potential of their protein surfaces, and developed a model summarizing the unique electrostatic interaction in LZ-8 dimerization. In addition, further electrostatic potential and virtual amino acid mutation analyses suggested that L10, W12, and D45 are the key amino acid residues responsible for the protein's high immunomodulatory activity. These findings may provide useful insights into the design and construction of a new FIP mutant for use in treating and preventing autoimmune diseases.

ITS1/5.8S/ITS2, a Good Marker for Initial Classification of Shiitake Culinary-Medicinal Lentinus edodes (Agaricomycetes) Strains in China.

China is home to rich wild and cultivated strains of Lentinus edodes, an important edible and medicinal mushroom. Artificial selection of L. edodes has a long history, and the widely cultivated strains belong to populations different from those of most wild strains. Internal transcribed spacer (ITS) regions have been used as good markers to identify L. edodes populations. But because ITS regions exhibit incomplete concerted evolution, the use of an ITS to identify L. edodes populations has been questioned. The objective of this study was to determine whether the ITS region is suitable for identifying L. edodes populations and which populations the widely cultivated strains and the most wild strains belong to by investigating intraindividual and differential ITS polymorphisms between 44 cultivars and 44 wild strains of L. edodes in China. Intraindividual ITS polymorphism is common in L. edodes strains, and most strains possessed 2 different ITS sequences, which came from their heterokaryons. The genetic polymorphisms of ITS1, 5.8S, and ITS2 in L. edodes strains are distinct. All strains were divided into one 5.8S type (5.8S-A), 2 ITS1 types (ITS1-A and ITS1-B), and 2 ITS2 types (ITS2-A and ITS2-B), which were subdivided into 2 branches (ITS2-A1 and ITS2-A2; ITS2-B1 and ITS2-B2). ITS1/5.8S/ITS2 could be used as a good marker in preliminary classification of L. edodes strains in China. It not only exhibited classified information of ITS1, 5.8S, and ITS2 for each strain at the same time, it also indicated whether the strain was heterozygous. The 44 cultivated strains were mainly the A/A/A1 type, and the 44 wild strains were mainly the A/A/A2 and other mixed types.

The complete mitochondrial genome of the Basidiomycete edible fungus Hypsizygus marmoreus.

The complete mitochondrial genome of the edible fungus Hypsizygus marmoreus was published in this paper. It was determined using Pacbio and Illumina sequencing. The complete mitochondrial DNA (mtDNA) is 106,417 bp in length with a GC content of 31.74%, which was the fourth large mitogenome in Agaricales. The circular mitogenome encoded 67 protein-coding genes and one ribosomal RNAs (rns). Among these genes, 13 conserved protein-coding genes were determined in the genome, including 6 subunits of NAD dehydrogenase (nad1-4, 4L and 6), three cytochrome oxidases (cox1-3), one apocytochrome b (cob) and three ATP synthases (atp6, apt 8 and apt 9). The phylogenic analysis confirmed that H. marmoreus (Lyophyllaceae) clustered together with Tricholoma matsutake (Tricholomataceae).

Comparison and Validation of Putative Pathogenicity-Related Genes Identified by T-DNA Insertional Mutagenesis and Microarray Expression Profiling in Magnaporthe oryzae.

High-throughput technologies of functional genomics such as T-DNA insertional mutagenesis and microarray expression profiling have been employed to identify genes related to pathogenicity in Magnaporthe oryzae. However, validation of the functions of individual genes identified by these high-throughput approaches is laborious. In this study, we compared two published lists of genes putatively related to pathogenicity in M. oryzae identified by T-DNA insertional mutagenesis (comprising 1024 genes) and microarray expression profiling (comprising 236 genes), respectively, and then validated the functions of some overlapped genes between the two lists by knocking them out using the method of target gene replacement. Surprisingly, only 13 genes were overlapped between the two lists, and none of the four genes selected from the overlapped genes exhibited visible phenotypic changes on vegetative growth, asexual reproduction, and infection ability in their knockout mutants. Our results suggest that both of the lists might contain large proportions of unrelated genes to pathogenicity and therefore comparing the two gene lists is hardly helpful for the identification of genes that are more likely to be involved in pathogenicity as we initially expected.

Comparative Proteomic Analysis of Light-Induced Mycelial Brown Film Formation in Lentinula edodes.

Light-induced brown film (BF) formation by the vegetative mycelium of Lentinula edodes is important for ensuring the quantity and quality of this edible mushroom. Nevertheless, the molecular mechanism underlying this phenotype is still unclear. In this study, a comparative proteomic analysis of mycelial BF formation in L. edodes was performed. Seventy-three protein spots with at least a twofold difference in abundance on two-dimensional electrophoresis (2DE) maps were observed, and 52 of them were successfully identified by matrix-assisted laser desorption/ionization tandem time-of-flight mass spectrometry (MALDI-TOF/TOF/MS). These proteins were classified into the following functional categories: small molecule metabolic processes (39%), response to oxidative stress (5%), and organic substance catabolic processes (5%), followed by oxidation-reduction processes (3%), single-organism catabolic processes (3%), positive regulation of protein complex assembly (3%), and protein metabolic processes (3%). Interestingly, four of the proteins that were upregulated in response to light exposure were nucleoside diphosphate kinases. To our knowledge, this is the first proteomic analysis of the mechanism of BF formation in L. edodes. Our data will provide a foundation for future detailed investigations of the proteins linked to BF formation.

The genome of Pleurotus eryngii provides insights into the mechanisms of wood decay.

Pleurotus eryngii (DC.) Quél. is widely used for bioconverting lignocellulosic byproducts into biofuel and value added products. Sequencing and annotating the genome of a monokaryon strain P. eryngii 183 allows us to gain a better understanding of carbohydrate-active enzymes (CAZymes) and oxidoreductases for degradation of lignocellulose in white-rot fungi. The genomic data provides insights into genomic basis of degradation mechanisms of lignin and cellulose and may pave new avenues for lignocellulose bioconversion.

Discovery and Characterization of the Highly Active Fungal Immunomodulatory Protein Fip-vvo82.

Volvaria volvacea (Bull. ex Fr.) Sing, an important edible and medicinal macro-fungus, has been used to remedy various diseases for hundreds of years in East Asia. To identify key proteins with the unique therapeutic activity in V. volvacea, we conducted a genomewide comparison of V. volvacea protein families and those of other edible fungi that lack therapeutic functions and identified seven fungal immunomodulatory proteins (FIPs) in V. volvacea. On the basis of the predicted physiological and biochemical characteristics of the seven FIPs, the novel Fip-vvo82 was inferred to have high immunomodulatory activity; this was confirmed by molecular and immunological experiments and further characterized by modeling the three-dimensional structure and protein-protein docking. This is the first study to show that V. volvacea has more than one FIP.

Anthropometric Characteristics of Chinese Professional Female Marathoners and Predicted Variables for Their Personal Bests.

To investigate anthropometric characteristics of Chinese professional female marathoners and suitable predicted variables correlated with their personal bests (PB), 96 Chinese female long-distance runners were divided into international (< 2 h 34 min), national (2 h 34 min~2 h 45 min) and average (2 h 45 min~3 h 19 min) levels according to their PB in marathon during the process of talent identification for London Olympic Games. Selected anthropometric variables, including height, body mass, percentages of body fat, girths, breadths, lengths and skin-folds were measured. Only iliac crest skin-fold of international athletes was significantly lower than it is in national group. Girth of forearm and lower limbs, length of lower limbs, and all skin-folds of national athletes were significantly lower than those from average level group. Percentages of body fat, girth of forearm and calf, length of lower limbs, and skin-folds at sites of subscapular, abdominal and iliac crest of athletes from average level group were significantly higher than those in international athletes. Positive correlation was found between forearm girth and PB, and between the subscapular, abdominal, iliac crest and triceps surae skin-folds and PB for total athletes. Negative correlation between biiliac breadth and PB in international athletes, and positive correlations between abdominal and triceps surae skin-folds and PB in national athletes were found. For average runners, high positive correlation was found between upper arm girth and PB, and between subscapular, abdominal, iliac crest and triceps surae skin-folds and PB. The findings suggested that compared to stride length, stride frequency and efficiency were more important factors influencing running performance, which were in accordance with running technique in Chinese female marathoners.