Comparative and transcriptional analysis of the predicted secretome in the lignocellulose-degrading basidiomycete fungusPleurotus ostreatus

Transcription factor Pofst3 regulates Pleurotus ostreatus development by targeting multiple biological pathways

Pleurotus ostreatus is a popular edible mushroom cultivated worldwide. However, the mechanism of P. ostreatus primordia formation is unclear. Pofst3 is a MHR superfamily transcription factor, which has the function of regulating primordia formation. In this study, the target genes of Pofst3 in P. ostreatus were identified by DAP-Seq approach at the genome level, 1481 peaks were obtained and the Pofst3 binding motif sequence was GARGRVGARGAR. The interaction between transcription factor Pofst3 and this motif was confirmed in vitro and in vivo through electrophoretic mobility shift (EMSA) and yeast one-hybrid screening (Y1H) assays. Among the top 20 GO enrichment results, most were related to transcriptional regulation, and some transcription factor encoding genes, such as HMG-box (gene_5346), MADS-box (gene_86), FOG (gene_6211) and RFX (gene_3496) were obtained. Besides basic metabolism, MAPK signaling pathway, Inositol phosphate metabolism, Glycosylphosphatidylinositol (GPI)-anchor biosynthesis and Pentose phosphate pathway were significantly enriched in the KEGG pathway analysis. The expression levels of randomly selected 11 genes, some transcription factor genes, and genes involved in metabolic pathways in wild and Pofst3 transgenic P. ostreatus strains indicated that target genes likely involved in the development of the P. ostreatus primordia. These results indicated that transcription factor Pofst3 ultimately negatively regulated the development of P. ostreatus primordia very likely through regulating a series of biological pathways.

Comparative analyses of Pleurotus pulmonarius mitochondrial genomes reveal two major lineages of mini oyster mushroom cultivars

Pleurotus pulmonarius, commonly known as the mini oyster mushroom, is highly esteemed for its crisp texture and umami flavor. Limited genetic diversity among P. pulmonarius cultivars raises concerns regarding its sustainable industrial production. To delve into the maternal genetic diversity of the principal P. pulmonarius cultivars, 36 cultivars and five wild isolates were subjected to de novo sequencing and assembly to generate high-quality mitogenome sequences. The P. pulmonarius mitogenomes had lengths ranging from 69,096 to 72,905 base pairs. The mitogenome sizes of P. pulmonarius and those of other mushroom species in the Pleurotus genus showed a significant positive correlation with the counts of LAGLIDAG and GIY-YIG homing endonucleases encoded by intronic open reading frames. A comparison of gene arrangements revealed an inversion of a fragment containing atp9-nad3-nad2 between P. pulmonarius and P. ostreatus. The mitogenomes of P. pulmonarius were clustered into three distinct clades, two of which were crowded with commercial cultivars. Clade I, all of which possess an inserted dpo gene, shared a maternal origin linked to an ancestral cultivar from Taiwan. Primers were designed to target the dpo gene, potentially safeguarding intellectual property rights. The wild isolates in Clade III exhibited more divergent mitogenomes, rendering them valuable for breeding.

Effects of hap2 deletion on mnp/vp transcription in Pleurotus ostreatus grown on lignocellulosic substrates

The regulatory mechanisms governing expression of genes encoding lignin-modifying enzymes (LME) in white-rot fungi remain largely unexplored. Although molecular cloning has identified CCAAT-boxes frequently located 5'-upstream of these genes, their role in transcriptional regulation is not well understood. This study examines the function of hap2, a gene encoding a hypothetical protein homologous to a component of the CCAAT-binding Hap complex, in the white-rot fungus Pleurotus ostreatus. Deletion of hap2 resulted in significantly reduced Mn2+-dependent peroxidase activity and lignin-degrading capacity compared to the parental strain 20b grown on beech wood sawdust (BWS) medium. Real-time PCR revealed that vp2 transcript levels were significantly lower in hap2 deletants than in 20b grown when cultured on the three solid media consisting of BWS, holocellulose, or Avicel, but not on yeast-malt-glucose (YMG) agar plates. Additionally, glutathione S-transferase (GST) pulldown and electrophoretic mobility shift assays demonstrated that recombinant P. ostreatus Hap2, Hap3, and Hap5 expressed in Escherichia coli form a complex capable of binding to the CCAAT sequence 5'-upstream of vp2 in vitro. These results suggest that Hap2, as part of the CCAAT-binding complex, is essential for transcriptional upregulation of vp2 in P. ostreatus growing on lignocellulosic substrates. KEY POINTS: • P. ostreatus hap2 deletants were generated. • Lignin-degrading capacity was significantly reduced in the hap2 deletants. • vp2 was significantly downregulated upon hap2 deletion.

Parascedosporium putredinis NO1 tailors its secretome for different lignocellulosic substrates

Parascedosporium putredinis NO1 is a plant biomass-degrading ascomycete with a propensity to target the most recalcitrant components of lignocellulose. Here we applied proteomics and activity-based protein profiling (ABPP) to investigate the ability of P. putredinis NO1 to tailor its secretome for growth on different lignocellulosic substrates. Proteomic analysis of soluble and insoluble culture fractions following the growth of P. putredinis NO1 on six lignocellulosic substrates highlights the adaptability of the response of the P. putredinis NO1 secretome to different substrates. Differences in protein abundance profiles were maintained and observed across substrates after bioinformatic filtering of the data to remove intracellular protein contamination to identify the components of the secretome more accurately. These differences across substrates extended to carbohydrate-active enzymes (CAZymes) at both class and family levels. Investigation of abundant activities in the secretomes for each substrate revealed similar variation but also a high abundance of "unknown" proteins in all conditions investigated. Fluorescence-based and chemical proteomic ABPP of secreted cellulases, xylanases, and β-glucosidases applied to secretomes from multiple growth substrates for the first time confirmed highly adaptive time- and substrate-dependent glycoside hydrolase production by this fungus. P. putredinis NO1 is a promising new candidate for the identification of enzymes suited to the degradation of recalcitrant lignocellulosic feedstocks. The investigation of proteomes from the biomass bound and culture supernatant fractions provides a more complete picture of a fungal lignocellulose-degrading response. An in-depth understanding of this varied response will enhance efforts toward the development of tailored enzyme systems for use in biorefining.IMPORTANCEThe ability of the lignocellulose-degrading fungus Parascedosporium putredinis NO1 to tailor its secreted enzymes to different sources of plant biomass was revealed here. Through a combination of proteomic, bioinformatic, and fluorescent labeling techniques, remarkable variation was demonstrated in the secreted enzyme response for this ascomycete when grown on multiple lignocellulosic substrates. The maintenance of this variation over time when exploring hydrolytic polysaccharide-active enzymes through fluorescent labeling, suggests that this variation results from an actively tailored secretome response based on substrate. Understanding the tailored secretomes of wood-degrading fungi, especially from underexplored and poorly represented families, will be important for the development of effective substrate-tailored treatments for the conversion and valorization of lignocellulose.

Altered Expression of Two Small Secreted Proteins (ssp4 and ssp6) Affects the Degradation of a Natural Lignocellulosic Substrate by Pleurotus ostreatus

Pleurotus ostreatus is a white-rot fungus that can degrade lignin in a preferential manner using a variety of extracellular enzymes, including manganese and versatile peroxidases (encoded by the vp1-3 and mnp1-6 genes, respectively). This fungus also secretes a family of structurally related small secreted proteins (SSPs) encoded by the ssp1-6 genes. Using RNA sequencing (RNA-seq), we determined that ssp4 and ssp6 are the predominant members of this gene family that were expressed by P. ostreatus during the first three weeks of growth on wheat straw. Downregulation of ssp4 in a strain harboring an ssp RNAi construct (KDssp1) was then confirmed, which, along with an increase in ssp6 transcript levels, coincided with reduced lignin degradation and the downregulation of vp2 and mnp1. In contrast, we observed an increase in the expression of genes related to pectin and side-chain hemicellulose degradation, which was accompanied by an increase in extracellular pectin-degrading capacity. Genome-wide comparisons between the KDssp1 and the wild-type strains demonstrated that ssp silencing conferred accumulated changes in gene expression at the advanced cultivation stages in an adaptive rather than an inductive mode of transcriptional response. Based on co-expression networking, crucial gene modules were identified and linked to the ssp knockdown genotype at different cultivation times. Based on these data, as well as previous studies, we propose that P. ostreatus SSPs have potential roles in modulating the lignocellulolytic and pectinolytic systems, as well as a variety of fundamental biological processes related to fungal growth and development.

A bioinformatic workflow for in silico secretome prediction with the lignocellulose degrading ascomycete fungus Parascedosporium putredinis NO1

The increasing availability of microbial genome sequences provides a reservoir of information for the identification of new microbial enzymes. Genes encoding proteins engaged in extracellular processes are of particular interest as these mediate the interactions microbes have with their environments. However, proteomic analysis of secretomes is challenging and often captures intracellular proteins released through cell death and lysis. Secretome prediction workflows from sequence data are commonly used to filter proteins identified through proteomics but are often simplified to a single step and are not evaluated bioinformatically for their effectiveness. Here, a workflow to predict a fungal secretome was designed and applied to the coding regions of the Parascedosporium putredinis NO1 genome. This ascomycete fungus is an exceptional lignocellulose degrader from which a new lignin-degrading enzyme has previously been identified. The 'secretome isolation' workflow is based on two strategies of localisation prediction and secretion prediction each utilising multiple available tools. The workflow produced three final secretomes with increasing levels of stringency. All three secretomes showed increases in functional annotations for extracellular processes and reductions in annotations for intracellular processes. Multiple sequences isolated as part of the secretome lacked any functional annotation and made exciting candidates for novel enzyme discovery.

Experimental evidence that lignin-modifying enzymes are essential for degrading plant cell wall lignin by Pleurotus ostreatus using CRISPR/Cas9

Lignin-modifying enzymes (LMEs), which include laccases (Lacs), manganese peroxidases (MnPs), versatile peroxidases (VPs), and lignin peroxidases (LiPs), have been considered key factors in lignin degradation by white-rot fungi because they oxidize lignin model compounds and depolymerize synthetic lignin in vitro. However, it remains unclear whether these enzymes are essential/important in the actual degradation of natural lignin in plant cell walls. To address this long-standing issue, we examined the lignin-degrading abilities of multiple mnp/vp/lac mutants of Pleurotus ostreatus. One vp2/vp3/mnp3/mnp6 quadruple-gene mutant was generated from a monokaryotic wild-type strain PC9 using plasmid-based CRISPR/Cas9. Also, two vp2/vp3/mnp2/mnp3/mnp6, two vp2/vp3/mnp3/mnp6/lac2 quintuple-gene mutants, and two vp2/vp3/mnp2/mnp3/mnp6/lac2 sextuple-gene mutants were generated. The lignin-degrading abilities of the sextuple and vp2/vp3/mnp2/mnp3/mnp6 quintuple-gene mutants on the Beech wood sawdust medium reduced drastically, but not so much for those of the vp2/vp3/mnp3/mnp6/lac2 mutants and the quadruple mutant strain. The sextuple-gene mutants also barely degraded lignin in Japanese Cedar wood sawdust and milled rice straw. Thus, this study presented evidence that the LMEs, especially MnPs and VPs, play a crucial role in the degradation of natural lignin by P. ostreatus for the first time.

Transcriptome Metabolic Characterization of Tuber borchii SP1-A New Spanish Strain for In Vitro Studies of the Bianchetto Truffle

Truffles are ascomycete hypogeous fungi belonging to the Tuberaceae family of the Pezizales order that grow in ectomycorrhizal symbiosis with tree roots, and they are known for their peculiar aromas and flavors. The axenic culture of truffle mycelium is problematic because it is not possible in many cases, and the growth rate is meager when it is possible. This limitation has prompted searching and characterizing new strains that can be handled in laboratory conditions for basic and applied studies. In this work, a new strain of Tuber borchii (strain SP1) was isolated and cultured, and its transcriptome was analyzed under different in vitro culture conditions. The results showed that the highest growth of T. borchii SP1 was obtained using maltose-enriched cultures made with soft-agar and in static submerged cultures made at 22 °C. We analyzed the transcriptome of this strain cultured in different media to establish a framework for future comparative studies, paying particular attention to the central metabolic pathways, principal secondary metabolite gene clusters, and the genes involved in producing volatile aromatic compounds (VOCs). The results showed a transcription signal for around 80% of the annotated genes. In contrast, most of the transcription effort was concentrated on a limited number of genes (20% of genes account for 80% of the transcription), and the transcription profile of the central metabolism genes was similar in the different conditions analyzed. The gene expression profile suggests that T. borchii uses fermentative rather than respiratory metabolism in these cultures, even in aerobic conditions. Finally, there was a reduced expression of genes belonging to secondary metabolite clusters, whereas there was a significative transcription of those involved in producing volatile aromatic compounds.

Silencing of the Laccase (lacc2) Gene from Pleurotus ostreatus Causes Important Effects on the Formation of Toxocyst-like Structures and Fruiting Body

A wide variety of biological functions, including those involved in the morphogenesis process of basidiomycete fungi, have been attributed to laccase enzymes. In this work, RNA interference (RNAi) was used to evaluate the role of the laccase (lacc2) gene of Pleurotus ostreatus PoB. Previously, transformant strains of P. ostreatus were obtained and according to their level of silencing they were classified as light (T7), medium (T21) or severe (T26 and T27). The attenuation of the lacc2 gene in these transformants was determined by RT-PCR. Silencing of lacc2 resulted in a decrease in laccase activity between 30 and 55%, which depended on the level of laccase expression achieved. The silenced strains (T21, T26, and T27) displayed a delay in the development of mycelium on potato dextrose agar (PDA) medium, whereas in the cultures grown on wheat straw, we found that these strains were incapable of producing aerial mycelium, primordia, and fruiting bodies. Scanning electron microscopy (SEM) showed the presence of toxocyst-like structures. The highest abundance of these structures was observed in the wild-type (PoB) and T7 strains. However, the abundance of toxocysts decreased in the T21 and T26 strains, and in T27 they were not detected. These results suggest that the presence and abundance of toxocyst-like structures are directly related to the development of fruiting bodies. Furthermore, our data confirm that lacc2 is involved in the morphogenesis process of P. ostreatus.

Pleurotus pulmonarius Strain: Arsenic(III)/Cadmium(II) Accumulation, Tolerance, and Simulation Application in Environmental Remediation

The arsenic (As, III) and cadmium (Cd, II) accumulation and tolerance traits of a new strain Pleurotus pulmonarius MT were evaluated, and the utilization of the strain for repairing contaminated liquid and soil was explored. The hypha cultivated in potato dextrose agar (PDA) exhibited medium or high Cd accumulation (0 to 320 mg/L), medium Cd tolerance (maximum tolerated concentration, MTC ≥ 640 mg/L), medium As accumulation (0 to 80 mg/L), and high As tolerance (MTC > 1280 mg/L). The hypha has application potential in processes related to the removal of Cd and As in aqueous pollutants at concentrations of 80 mg/L Cd and 20 mg/L As. The trends obtained for the fruiting bodies of P. pulmonarius MT seemed to deviate from those of the hypha of this strain. The results show that the fruiting bodies featured medium As accumulation (0 to 40 mg/kg), medium As tolerance (MTC > 160 mg/kg), medium Cd accumulation (0 to 10 mg/kg), and high Cd tolerance (MTC > 1280 mg/kg). The fruiting bodies of P. pulmonarius MT were utilized in processes related to the recovery of Cd and As in substrates, that is, 12% contaminated soil mixed with 50 mg/kg Cd and 200 mg/kg As; thus, the hypha and fruiting bodies of P. pulmonarius MT can be used for the decontamination of water and soil containing As(III) and Cd(II).

Biocatalyst collection and heterologous expression of sesquiterpene synthases from basidiomycetous fungi: Discovery of a novel sesquiterpene hydrocarbon

Basidiomycetes produce a wide variety of sesquiterpenoids, which attract significant interest in pharmaceutical and industrial applications. Structural diversification of sesquiterpenoids is performed by sesquiterpene synthases (STSs), which produce a wide array of backbone structures; therefore, functional characterization and increased biocatalyst collection of STSs are important for expanding scientific knowledge and meeting the needs of advanced biotechnology. Gene identification and functional annotation of STSs from the basidiomycetous fungi Agaricus bisporus, Auriscalpium vulgare, Lepista nuda, Pleurotus ostreatus and Trametes versicolor were conducted. Through these investigations, the catalytic functions of 30 STSs were revealed using recombinant enzymes heterologously expressed in Saccharomyces cerevisiae. Furthermore, the unique function of an STS from P. ostreatus, PoSTS-06, was revealed to be the production of a novel sesquiterpene hydrocarbon that we named pleostene. The absolute structure of pleostene was determined by NMR spectroscopy and X-ray crystallography using the crystalline sponge method.

Purification and Characterization of Class III Lipase from a White-Rot Fungus Pleurotus ostreatus

Pleurotus ostreatus is an edible white-rot fungus with lignocellulosic biomass degrading enzymes that have been studied extensively. However, until now, lipolytic enzymes from P. ostreatus, which degrade extractives in lignocellulosic biomass, have not been purified and characterized. In this study, P. ostreatus was inoculated into the rapeseed oil containing culture to induce lipase. The lipase in the culture broth was successfully purified to homogeneity by chromatographic methods. The molecular weight of the purified lipase was 27 kDa, and its optimal pH and temperature were 5.0 and 30 °C, respectively. The purified lipase showed high activity with the substrates 4-methylumbelliferyl (4-MU) decanoate (C10:0) and 4-MU oleate (C18:1), and no activity with 4-MU acetate (C2:0) and 4-MU butyrate (C4:0). The amino acid sequences and substrate specificities of the purified lipase suggested that it belonged to class III. Kinetic parameters measurements (Km and Vmax) showed that 4-MU palmitate had a high affinity for the purified lipase, and it was the substrate most efficiently hydrolyzed by the purified lipase.

Application of omics technology in the research on edible fungi

Edible fungus is a large fungus distributed all over the world and used as food and medicine. But people's understanding of edible fungi is not as much as that of ordinary crops, so people have started a number of research on edible fungi in recent years. With the development of science and technology, omics technology has gradually walked into people's vision. Omics technology has high sensitivity and wide application range, which is favored by researchers. The application of omics technology to edible fungus research is a major breakthrough, which has transferred edible fungus research from artificial cultivation to basic research. Now omics technology in edible fungi has been flexibly combined with other research methods, involving multiple studies of edible fungus, such as genetic breeding, growth and development, stress resistance, and the use of special components in edible fungus as pharmaceutical additives. It is believed that in the future, the research of edible fungi will also be brought to a deeper level with the help of omics technology. This paper introduces the application progress of modern omics technology to the study on edible fungi and mentions the application prospect of edible fungi research with the constant development of omics technology, thereby providing ideas for the follow-up in-depth research on edible fungi.

Towards Understanding the Function of Aegerolysins

Aegerolysins are remarkable proteins. They are distributed over the tree of life, being relatively widespread in bacteria and fungi, but also present in some insects, plants, protozoa, and viruses. Despite their abundance in cells of certain developmental stages and their presence in secretomes, only a few aegerolysins have been studied in detail. Their function, in particular, is intriguing. Here, we summarize previously published findings on the distribution, molecular interactions, and function of these versatile aegerolysins. They have very diverse protein sequences but a common fold. The machine learning approach of the AlphaFold2 algorithm, which incorporates physical and biological knowledge of protein structures and multisequence alignments, provides us new insights into the aegerolysins and their pore-forming partners, complemented by additional genomic support. We hypothesize that aegerolysins are involved in the mechanisms of competitive exclusion in the niche.

Bacterial Infection Induces Ultrastructural and Transcriptional Changes in the King Oyster Mushroom (Pleurotus eryngii)

Pleurotus eryngii (king oyster mushroom) is a commercially important mushroom with high nutritional and economic value. However, soft rot disease, caused by the pathogenic bacterium Erwinia beijingensis, poses a threat to its quality and production. Morphological and ultrastructural observations of P. eryngii were conducted at early, middle, and late stages of infection. At 2 days postinoculation (dpi), small yellow spots on the fruiting body were observed. The infected tissue displayed hyphal deformations and breaks at 5 dpi. At 9 dpi, damage to cell wall integrity and absence of intact cellular organelles were observed and the diseased fruiting bodies were unable to grow normally. Transcriptome analysis identified 4,296 differentially expressed genes in the fruiting body following infection. In fact, broad transcriptional reprogramming was observed in infected fruiting bodies compared to controls. The affected pathways included antioxidant systems, peroxisome biogenesis, autophagy, and oxidation-reduction. More specifically, pex genes were downregulated during infection, indicating impaired peroxisome homeostasis and redox balance. Additionally, genes encoding chitinase, β-1,3-glucanase, and proteases associated with cell wall degradation were upregulated in infected P. eryngii. This study provides insights into the responses of P. eryngii during soft rot disease and facilitates the understanding of the pathogenic process of bacteriosis in mushrooms.

IMPORTANCEPleurotus eryngii (king oyster mushroom) is a popular and economically valuable edible mushroom; however, it suffers from various bacterial diseases, including soft rot disease caused by the bacterium Erwinia beijingensis. Here, we examined bacterial infection of the mushroom through morphological and ultrastructural observations as well as transcriptome analysis. Pathogen attack damaged the cell structure of P. eryngii, including the cell wall, and also induced high levels of reactive oxygen species. These results were reflected in differential gene expression in P. eryngii as a response to the pathogenic bacteria, including genes involved in antioxidant systems, peroxisome biogenesis, autophagy, oxidation-reduction, ribosome biogenesis, and cell-wall degradation, among others. This study provides insights into the structural and molecular responses of P. eryngii during soft rot disease, improving our understanding and the potential control of the pathogenic process of bacteriosis in mushrooms.

Influence of Different Vegetation Types on Soil Physicochemical Parameters and Fungal Communities

This study assessed the effects of Betula dahurica (BD), Betula platyphylla (BP), Larix gmelinii (LG), Quercus mongolica (QM), and a mixed conifer–broadleaf forest composed of LG and QM (LGQM) on the soil physicochemical parameters and community structure of fungi in the Zhongyangzhan Black-billed Capercaillie Nature Reserve. Fungal community structures were characterized via ITS rRNA sequencing. The effects of soil parameters on the community structure of soil fungi were assessed by Pearson correlation analysis and redundancy analysis (RDA). LGQM exhibited lower C/N, available nitrogen (AN), total phosphorus (TP), and available phosphorus (AP) compared with the QM broadleaf forest. The fungal Shannon and Simpson diversity indices were highest in BP, whereas LG exhibited the highest ACE index. The Basidiomycota, Ascomycota, Mortierellomycota, and Mucoromycota fungal phyla were dominant across all vegetation types. Each of the different vegetation types studied herein exhibited a unique fungal community structure. The RDA results indicated that fungal community structures were primarily shaped by the total N, available N, and available P of soil. Our findings thus indicated that forests restored with different species of trees may exhibit variations in soil quality and characteristics despite sharing the same climate. Furthermore, broadleaved and coniferous forests exhibited a unique fungal community diversity and composition.

Basidiomycota Fungi and ROS: Genomic Perspective on Key Enzymes Involved in Generation and Mitigation of Reactive Oxygen Species

Our review includes a genomic survey of a multitude of reactive oxygen species (ROS) related intra- and extracellular enzymes and proteins among fungi of Basidiomycota, following their taxonomic classification within the systematic classes and orders, and focusing on different fungal lifestyles (saprobic, symbiotic, pathogenic). Intra- and extracellular ROS metabolism-involved enzymes (49 different protein families, summing 4170 protein models) were searched as protein encoding genes among 63 genomes selected according to current taxonomy. Extracellular and intracellular ROS metabolism and mechanisms in Basidiomycota are illustrated in detail. In brief, it may be concluded that differences between the set of extracellular enzymes activated by ROS, especially by H2O2, and involved in generation of H2O2, follow the differences in fungal lifestyles. The wood and plant biomass degrading white-rot fungi and the litter-decomposing species of Agaricomycetes contain the highest counts for genes encoding various extracellular peroxidases, mono- and peroxygenases, and oxidases. These findings further confirm the necessity of the multigene families of various extracellular oxidoreductases for efficient and complete degradation of wood lignocelluloses by fungi. High variations in the sizes of the extracellular ROS-involved gene families were found, however, among species with mycorrhizal symbiotic lifestyle. In addition, there are some differences among the sets of intracellular thiol-mediation involving proteins, and existence of enzyme mechanisms for quenching of intracellular H2O2 and ROS. In animal- and plant-pathogenic species, extracellular ROS enzymes are absent or rare. In these fungi, intracellular peroxidases are seemingly in minor role than in the independent saprobic, filamentous species of Basidiomycota. Noteworthy is that our genomic survey and review of the literature point to that there are differences both in generation of extracellular ROS as well as in mechanisms of response to oxidative stress and mitigation of ROS between fungi of Basidiomycota and Ascomycota.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Lipid-Binding Aegerolysin from Biocontrol Fungus Beauveria bassiana

Fungi are the most common pathogens of insects and thus important regulators of their populations. Lipid-binding aegerolysin proteins, which are commonly found in the fungal kingdom, may be involved in several biologically relevant processes including attack and defense against other organisms. Aegerolysins act alone or together with membrane-attack-complex/perforin (MACPF)-like proteins to form transmembrane pores that lead to cell lysis. We performed an in-depth bioinformatics analysis of aegerolysins in entomopathogenic fungi and selected a candidate aegerolysin, beauveriolysin A (BlyA) from Beauveria bassiana. BlyA was expressed as a recombinant protein in Escherichia coli, and purified to further determine its functional and structural properties, including lipid-binding ability. Aegerolysins were found to be encoded in genomes of entomopathogenic fungi, such as Beauveria, Cordyceps, Metarhizium and Ophiocordyceps. Detailed bioinformatics analysis revealed that they are linked to MACPF-like genes in most genomes. We also show that BlyA interacts with an insect-specific membrane lipid. These results were placed in the context of other fungal and bacterial aegerolysins and their partner proteins. We believe that aegerolysins play a role in promoting the entomopathogenic and antagonistic activity of B. bassiana, which is an active ingredient of bioinsecticides.

Assessing the Biodegradation of Vulcanised Rubber Particles by Fungi Using Genetic, Molecular and Surface Analysis

Millions of tonnes of tyre waste are discarded annually and are considered one of the most difficult solid wastes to recycle. A sustainable alternative for the treatment of vulcanised rubber is the use of microorganisms that can biotransform polymers and aromatic compounds and then assimilate and mineralise some of the degradation products. However, vulcanised rubber materials present great resistance to biodegradation due to the presence of highly hydrophobic cross-linked structures that are provided by the additives they contain and the vulcanisation process itself. In this work, the biodegradation capabilities of 10 fungal strains cultivated in PDA and EM solid medium were studied over a period of 4 weeks. The growth of the strains, the mass loss of the vulcanised rubber particles and the surface structure were analysed after the incubation period. With the white rot fungi Trametes versicolor and Pleurotus ostreatus, biodegradation percentages of 7.5 and 6.1%, respectively, were achieved. The FTIR and SEM-EDS analyses confirmed a modification of the abundance of functional groups and elements arranged on the rubber surface, such as C, O, S, Si, and Zn, due to the biological treatment employed. The availability of genomic sequences of P. ostreatus and T. versicolor in public repositories allowed the analysis of the genetic content, genomic characteristics and specific components of both fungal species, determining some similarities between both species and their relationship with rubber biodegradation. Both fungi presented a higher number of sequences for laccases and manganese peroxidases, two extracellular enzymes responsible for many of the oxidative reactions reported in the literature. This was confirmed by measuring the laccase and peroxidase activity in cultures of T. versicolor and P. ostreatus with rubber particles, reaching between 2.8 and 3.3-times higher enzyme activity than in the absence of rubber. The integrative analysis of the results, supported by genetic and bioinformatics tools, allowed a deeper analysis of the biodegradation processes of vulcanised rubber. It is expected that this type of analysis can be used to find more efficient biotechnological solutions in the future.

Strain Degeneration in Pleurotus ostreatus: A Genotype Dependent Oxidative Stress Process Which Triggers Oxidative Stress, Cellular Detoxifying and Cell Wall Reshaping Genes

Strain degeneration has been defined as a decrease or loss in the yield of important commercial traits resulting from subsequent culture, which ultimately leads to Reactive Oxygen Species (ROS) production. Pleurotus ostreatus is a lignin-producing nematophagous edible mushroom. Mycelia for mushroom production are usually maintained in subsequent culture in solid media and frequently show symptoms of strain degeneration. The dikaryotic strain P. ostreatus (DkN001) has been used in our lab as a model organism for different purposes. Hence, different tools have been developed to uncover genetic and molecular aspects of this fungus. In this work, strain degeneration was studied in a full-sib monokaryotic progeny of the DkN001 strain with fast (F) and slow (S) growth rates by using different experimental approaches (light microscopy, malondialdehyde levels, whole-genome transcriptome analysis, and chitosan effect on monokaryotic mycelia). The results obtained showed that: (i) strain degeneration in P. ostreatus is linked to oxidative stress, (ii) the oxidative stress response in monokaryons is genotype dependent, (iii) stress and detoxifying genes are highly expressed in S monokaryons with symptoms of strain degeneration, (iv) chitosan addition to F and S monokaryons uncovered the constitutive expression of both oxidative stress and cellular detoxifying genes in S monokaryon strains which suggest their adaptation to oxidative stress, and (v) the overexpression of the cell wall genes, Uap1 and Cda1, in S monokaryons with strain degeneration phenotype indicates cell wall reshaping and the activation of High Osmolarity Glycerol (HOG) and Cell Wall Integrity (CWI) pathways. These results could constitute a hallmark for mushroom producers to distinguish strain degeneration in commercial mushrooms.

Overexpressing Pleurotus ostreatus rho1b results in transcriptional upregulation of the putative cellulolytic enzyme-encoding genes observed in ccl1 disruptants

The transcriptional expression pattern of lignocellulolytic enzyme-encoding genes in white-rot fungi differs depending on the culture conditions. Recently, it was shown that 13 putative cellulolytic enzyme-encoding genes were significantly upregulated in most Pleurotus ostreatus ligninolysis-deficient mutant strains on beech wood sawdust medium. However, the mechanisms by which this transcriptional shift is triggered remain unknown. In this study, we identified one mechanism. Our previous study implied that histone H3 N-dimethylation at lysine 4 level possibly affects the shift; therefore, we analysed the expression pattern in the disruptants of P. ostreatus ccl1, which encodes a putative component of the COMPASS complex mediating the methylation. The results showed upregulation of 5 of the 13 cellulolytic enzyme-encoding genes. We also found that rho1b, encoding a putative GTPase regulating signal transduction pathways, was upregulated in the ccl1 disruptants and ligninolysis-deficient strains. Upregulation of at least three of the five cellulolytic enzyme-encoding genes was observed in rho1b-overexpressing strains but not in ccl1/rho1b double-gene disruptants, during the 20-day culture period. These results suggest that Rho1b may be involved in the upregulation of cellulolytic enzyme-encoding genes observed in the ccl1 disruptants. Furthermore, we suggest that Mpk1b, a putative Agaricomycetes-specific mitogen-activated protein kinase, functions downstream of Rho1b.

The growth of marine fungi on seaweed polysaccharides produces cerato-platanin and hydrophobin self-assembling proteins

The marine fungi Paradendryphiela salina and Talaromyces pinophilus degrade and assimilate complex substrates from plants and seaweed. Additionally, these fungi secrete surface-active proteins, identified as cerato-platanins and hydrophobins. These hydrophobic proteins have the ability to self-assemble forming amyloid-like aggregates and play an essential role in the growth and development of the filamentous fungi. It is the first time that one cerato-platanin (CP) is identified and isolated from P. salina (PsCP) and two Class I hydrophobins (HFBs) from T. pinophilus (TpHYD1 and TpHYD2). Furthermore, it is possible to extract cerato-platanins and hydrophobins using marine fungi that can feed on seaweed biomass, and through a submerged liquid fermentation process. The propensity to aggregate of these proteins has been analyzed using different techniques such as Thioflavin T fluorescence assay, Fourier-transform Infrared Spectroscopy, and Atomic Force Microscopy. Here, we show that the formation of aggregates of PsCP and TpHYD, was influenced by the carbon source from seaweed. This study highlighted the potential of these self-assembling proteins generated from a fermentation process with marine fungi and with promising properties such as conformational plasticity with extensive applications in biotechnology, pharmacy, nanotechnology, and biomedicine.

Genome sequence of the oyster mushroom Pleurotus ostreatus strain PC9

The oyster mushroom Pleurotus ostreatus is a basidiomycete commonly found in the rotten wood and it is one of the most cultivated edible mushrooms globally. Pleurotus ostreatus is also a carnivorous fungus, which can paralyze and kill nematodes within minutes. However, the molecular mechanisms of the predator-prey interactions between P. ostreatus and nematodes remain unclear. PC9 and PC15 are two model strains of P. ostreatus and the genomes of both strains have been sequenced and deposited at the Joint Genome Institute (JGI). These two monokaryotic strains exhibit dramatic differences in growth, but because PC9 grows more robustly in laboratory conditions, it has become the strain of choice for many studies. Despite the fact that PC9 is the common strain for investigation, its genome is fragmentary and incomplete relative to that of PC15. To overcome this problem, we used PacBio long reads and Illumina sequencing to assemble and polish a more integrated genome for PC9. Our PC9 genome assembly, distributed across 17 scaffolds, is highly contiguous and includes five telomere-to-telomere scaffolds, dramatically improving the genome quality. We believe that our PC9 genome resource will be useful to the fungal research community investigating various aspects of P. ostreatus biology.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

The synergistic actions of hydrolytic genes reveal the mechanism of Trichoderma harzianum for cellulose degradation

Bioprospecting genes and proteins related to plant biomass degradation is an attractive approach for the identification of target genes for biotechnological purposes, especially those with potential applications in the biorefinery industry that can enhance second-generation ethanol production technology. Trichoderma harzianum is a potential candidate for cellulolytic enzyme prospection and production. Herein, the enzymatic activities, transcriptome, exoproteome, and coexpression networks of the T. harzianum strain CBMAI-0179 were examined under biomass degradation conditions. We identified differentially expressed genes (DEGs) and carbohydrate-active enzyme (CAZyme) genes related to plant biomass degradation and compared them with those of strains from congeneric species (T. harzianum IOC-3844 and T. atroviride CBMAI-0020). T. harzianum CBMAI-0179 harbors strain- and treatment-specific CAZyme genes and transcription factors. We detected important proteins related to biomass degradation, including β-glucosidases, endoglucanases, cellobiohydrolases, lytic polysaccharide monooxygenases, endo-1,4-β-xylanases and β-mannanases. Based on coexpression networks, an enriched cluster with degradative enzymes was described, and the subnetwork of CAZymes revealed strong correlations among important secreted proteins and differentially expressed CAZyme genes. Our results provide valuable information for future studies on the genetic regulation of plant cell wall-degrading enzymes. This knowledge can be exploited for the improvement of enzymatic reactions in biomass degradation for bioethanol production.

Genomic analyses of two Italian oyster mushroom Pleurotus pulmonarius strains

Pleurotus mushrooms are among the most cultivated fungi in the world and are highly valuable for food, medicine, and biotechnology industries. Furthermore, Pleurotus species are carnivorous fungi; they can rapidly paralyze and kill nematodes when nutrient-deprived. The predator-prey interactions between Pleurotus and nematodes are still widely unexplored. Moreover, the molecular mechanisms and the genes involved in the carnivorous behavior of Pleurotus mushrooms remain a mystery. We are attempting to understand the interactions between Pleurotus mushrooms and their nematode prey through genetic and genomic analyses. Two single spores (ss2 and ss5) isolated from a fruiting body of Pleurotus pulmonarius exhibited significant differences in growth and toxicity against nematodes. Thus, using PacBio long reads, we assembled and annotated two high-quality genomes for these two isolates of P. pulmonarius. Each of these assemblies contains 23 scaffolds, including 6 (ss2) and 8 (ss5) telomere-to-telomere scaffolds, and they are among the most complete assembled genomes of the Pleurotus species. Comparative analyses identified the genomic differences between the two P. pulmonarius strains. In sum, this work provides a genomic resource that will be invaluable for better understanding the Italian oyster mushroom P. pulmonarius.

Secretomic insight into the biomass hydrolysis potential of the phytopathogenic fungus Chrysoporthe cubensis

The phytopathogenic fungus Chrysoporthe cubensis has a great capacity to produce highly efficient enzymes for the hydrolysis of lignocellulosic biomass. The bioinfosecretome of C. cubensis was identified by computational predictions of secreted proteins combined with protein analysis using 1D-LC-MS/MS. The in silico secretome predicted 562 putative genes capable of encoding secreted proteins, including 273 CAZymes. Proteomics analysis confirmed the existence of 313 proteins, including 137 CAZymes classified as Glycosyl Hydrolases (GH), Polysaccharide Lyases (PL), Carbohydrate Esterases (CE) and Auxiliary Activities enzymes (AA), which indicates the presence of classical and oxidative cellulolytic mechanisms. The enzymes diversity in the extract shows fungal versatility to act in complex biomasses. This study provides an insight into the lignocellulose-degradation mechanisms by C. cubensis and allows the identification of the enzymes that are potentially useful in improving industrial process of bioconversion of lignocellulose. SIGNIFICANCE: Chrysoporthe cubensis is an important deadly canker pathogen of commercially cultivated Eucalyptus species. The effective depolymerisation of the recalcitrant plant cell wall performed by this fungus is closely related to its high potential of lignocellulolytic enzymes secretion. Since the degradation of biomass occurs in nature almost exclusively by enzyme secretion systems, it is reasonable to suggest that the identification of C. cubensis lignocellulolytic enzymes is relevant in contributing to new sustainable alternatives for industrial solutions. As far as we know, this work is the first accurate proteomic evaluation of the enzymes secreted by this species of fungus. The integration of the gel-based proteomic approach, the bioinformatic prediction of the secretome and the analyses of enzymatic activity are powerful tools in the evaluation of biotechnological potential of C. cubensis in producing carbohydrate-active enzymes. In addition, analysis of the C. cubensis secretome grown in wheat bran draws attention to this plant pathogen and its extracellular enzymatic machinery, especially regarding the identification of promising new enzymes for industrial applications. The results from this work allowed for explanation and reinforce previous research that revealed C. cubensis as a strong candidate to produce enzymes to hydrolyse sugarcane bagasse and similar substrates.

Comparative transcriptional analyses of Pleurotus ostreatus mutants on beech wood and rice straw shed light on substrate-biased gene regulation

Distinct wood degraders occupying their preferred habitats have biased enzyme repertoires that are well fitted to their colonized substrates. Pleurotus ostreatus, commonly found on wood, has evolved its own enzyme-producing traits. In our previous study, transcriptional shifts in several P. ostreatus delignification-defective mutants, including Δhir1 and Δgat1 strains, were analyzed, which revealed the downregulation of ligninolytic genes and the upregulation of cellulolytic and xylanolytic genes when compared to their parental strain 20b on beech wood sawdust medium (BWS). In this study, rice straw (RS) was used as an alternative substrate to examine the transcriptional responses of P. ostreatus to distinct substrates. The vp1 gene and a cupredoxin-encoding gene were significantly upregulated in the 20b strain on RS compared with that on BWS, reflecting their distinct regulation patterns. The overall expression level of genes encoding glucuronidases was also higher on RS than on BWS, showing a good correlation with the substrate composition. Transcriptional alterations in the mutants (Δhir1 or Δgat1 versus 20b strain) on RS were similar to those on BWS, and the extracellular lignocellulose-degrading enzyme activities and lignin-degrading ability of the mutants on RS were consistent with the transcriptional alterations of the corresponding enzyme-encoding genes. However, transcripts of specific genes encoding enzymes belonging to the same CAZyme family exhibited distinct alteration patterns in the mutant strains grown on RS compared to those grown on BWS. These findings provide new insights into the molecular mechanisms underlying the transcriptional regulation of lignocellulolytic genes in P. ostreatus.Key Points• P. ostreatus expressed variable enzymatic repertoire-related genes in response to distinct substrates.• A demand to upregulate the cellulolytic genes seems to be present in ligninolysis-deficient mutants.• The regulation of some specific genes probably driven by the demand is dependent on the substrate.

Targeted disruption of hir1 alters the transcriptional expression pattern of putative lignocellulolytic genes in the white-rot fungus Pleurotus ostreatus

Pleurotus ostreatus is frequently used in molecular genetics and genomic studies on white-rot fungi because various molecular genetic tools and relatively well-annotated genome databases are available. To explore the molecular mechanisms underlying wood lignin degradation by P. ostreatus, we performed mutational analysis of a newly isolated mutant UVRM28 that exhibits decreased lignin-degrading ability on the beech wood sawdust medium. We identified that a mutation in the hir1 gene encoding a putative histone chaperone, which probably plays an important role in DNA replication-independent nucleosome assembly, is responsible for the mutant phenotype. The expression pattern of ligninolytic genes was altered in hir1 disruptants. The most highly expressed gene vp2 was significantly inactivated, whereas the expression of vp1 was remarkably upregulated (300-400 fold) at the transcription level. Conversely, many cellulolytic and xylanolytic genes were upregulated in hir1 disruptants. Chromatin immunoprecipitation analysis suggested that the histone modification status was altered in the 5'-upstream regions of some of the up- and down-regulated lignocellulolytic genes in hir1 disruptants compared with that in the 20b strain. Hence, our data provide new insights into the regulatory mechanisms of lignocellulolytic genes in P. ostreatus.

Interruption of an MSH4 homolog blocks meiosis in metaphase I and eliminates spore formation in Pleurotus ostreatus

Pleurotus ostreatus, one of the most widely cultivated edible mushrooms, produces high numbers of spores causing severe respiratory health problems for people, clogging of filters and spoilage of produce. A non-sporulating commercial variety (SPOPPO) has been successfully introduced into the market in 2006. This variety was generated by introgression breeding of a natural mutation into a commercial variety. Our cytological studies revealed that meiosis in the natural and derived sporeless strains was blocked in metaphase I, apparently resulting in a loss of spore formation. The gene(s) underlying this phenotype were mapped to an 80 kb region strongly linked to sporelessness and identified by transformation of wild type genes of this region into a sporeless strain. Sporulation was restored by re-introduction of the DNA sequence encoding the P. ostreatus meiotic recombination gene MSH4 homolog (poMSH4). Subsequent molecular analysis showed that poMSH4 in the sporeless P. ostreatus was interrupted by a DNA fragment containing a region encoding a CxC5/CxC6 cysteine cluster associated with Copia-type retrotransposons. The block of meiosis in metaphase I by a poMSH4 null mutant suggests that this protein plays an essential role in both Class I and II crossovers in mushrooms, similar to animals (mice), but unlike in plants. MSH4 was previously shown to be a target for breeding of sporeless varieties in P. pulmonarius, and the null mutant of the MSH4 homolog of S. commune (scMSH4) confers an extremely low level of spore formation. We propose that MSH4 homologs are likely to be a breeding target for sporeless strains both within Pleurotus sp. and in other Agaricales.

De novo assembly of Auricularia polytricha transcriptome and discovery of genes involved in the degradation of lignocellulose

Auricularia polytricha belonging to Basidiomycota has the ability to degrade lignocellulose. However, there has been no resource in public databases examining the transcriptome of A. polytricha. In this study, high-throughput sequencing platform BGISEQ-500 was used to generate large amount of transcript sequences from A. polytricha for gene discovery and molecular marker development. A total of 28,102 unigenes were discovered from the assembly of clean reads. In addition, functional categorization of the gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) metabolic pathways revealed several important biological processes. GO annotation analysis presented 47 categories, with the major subcategories being catalytic activity, binding, cellular process, metabolic process, and cell. Among the five functional categories and 21 subcategories of processes discovered from KEGG, global and overview maps, carbohydrate metabolism, transport, and catabolism are the main subcategories. Furthermore, among the unigenes related to lignocellulosic degradation discovered by KEGG pathway enrichment analysis, 2, 5, and 16 unigenes in de novo assembly of A. polytricha transcriptome were found to relate to cellulose, hemicellulose, and lignin degradation, respectively. The study provided valuable information on the degradation of lignocellulose to facilitate research on the degradation mechanism, molecular marker, functional research, gene mapping, and other multigenomic studies of species containing lignocellulose.

Secretomes of medically important fungi reflect morphological and phylogenetic diversity

Secretome represents a main target for understanding the mechanisms of fungal adaptation. In the present study, we focus on the secretomes of fungi associated with infections in humans and other mammals in order to explore relationships between the diverse morphological and phylogenetic groups. Almost all the mammalian pathogenic fungi analyzed have secretome sizes smaller than 1000 proteins and, secreted proteins comprise between 5% and 10% of the total proteome. As expected, the correlation pattern between the secretome size and the total proteome was similar to that described in previous secretome studies of fungi. With regard to the morphological groups, minimum secretome sizes of less than 250 secreted proteins and low values for the fraction of secreted proteins are shown in mammalian pathogenic fungi with reduced proteomes such as microsporidia, atypical fungi and some species of yeasts and yeast-like fungi (Malassezia). On the other hand, filamentous fungi have significantly more secreted proteins and the highest numbers are present in species of filamentous fungi that also are plant or insect pathogens (Fusarium verticilloides, Fusarium oxysporum and Basidiobolus meristosporus). With respect to phylogeny, there are also variations in secretome size across fungal subphyla: Microsporidia, Taphrinomycotina, Ustilagomycotina and Saccharomycotina contain small secretomes; whereas larger secretomes are found in Agaricomycotina, Pezizomycotina, Mucoromycotina and Entomophthoromycotina. Finally, principal component analysis (PCA) was conducted on the complete secretomes. The PCA results revealed that, in general, secretomes of fungi belonging to the same morphological group or subphyla cluster together. In conclusion, our results point out that in medically important fungi there is a relationship between the secretome and the morphological group or phylogenetic classification.

Transcriptional shifts in delignification-defective mutants of the white-rot fungus Pleurotus ostreatus

White-rot fungi efficiently degrade lignin and, thus, play a pivotal role in the global carbon cycle. However, the mechanisms of lignin degradation are largely unknown. Recently, mutations in four genes, namely wtr1, chd1, pex1, and gat1, were shown to abrogate the wood lignin-degrading ability of Pleurotus ostreatus. In this study, we conducted a comparative transcriptome analysis to identify genes that are differentially expressed in ligninolysis-deficient mutant strains. Putative ligninolytic genes that are highly expressed in parental strains are significantly downregulated in the mutant strains. On the contrary, many putative cellulolytic and xylanolytic genes are upregulated in the chd1-1, Δpex1, and Δgat1 strains. Identifying transcriptional alterations in mutant strains could provide new insights into the regulatory mechanisms of lignocellulolytic genes in P. ostreatus.

Glucose counteracts wood-dependent induction of lignocellulolytic enzyme secretion in monokaryon and dikaryon submerged cultures of the white-rot basidiomycete Pleurotus ostreatus

The secretome complexity and lignocellulose degrading capacity of Pleurotus ostreatus monokaryons mkPC9 and mkPC15 and mated dikaryon dkN001 were studied in submerged liquid cultures containing wood, glucose, and wood plus glucose as carbon sources. The study revealed that this white-rot basidiomycete attacks all the components of the plant cell wall. P. ostreatus secretes a variety of glycoside hydrolases, carbohydrate esterases, and polysaccharide lyases, especially when wood is the only carbon source. The presence of wood increased the secretome complexity, whereas glucose diminished the secretion of enzymes involved in cellulose, hemicellulose and pectin degradation. In contrast, the presence of glucose did not influence the secretion of redox enzymes or proteases, which shows the specificity of glucose on the secretion of cellulolytic enzymes. The comparison of the secretomes of monokaryons and dikaryons reveals that secretome complexity is unrelated to the nuclear composition of the strain.

Whole-genome and time-course dual RNA-Seq analyses reveal chronic pathogenicity-related gene dynamics in the ginseng rusty root rot pathogen Ilyonectria robusta

Ilyonectria robusta causes rusty root rot, the most devastating chronic disease of ginseng. Here, we for the first time report the high-quality genome of the I. robusta strain CD-56. Time-course (36 h, 72 h, and 144 h) dual RNA-Seq analysis of the infection process was performed, and many genes, including candidate effectors, were found to be associated with the progression and success of infection. The gene expression profile of CD-56 showed a trend of initial inhibition and then gradually returned to a profile similar to that of the control. Analyses of the gene expression patterns and functions of pathogenicity-related genes, especially candidate effector genes, indicated that the stress response changed to an adaptive response during the infection process. For ginseng, gene expression patterns were highly related to physiological conditions. Specifically, the results showed that ginseng defenses were activated by CD-56 infection and persisted for at least 144 h thereafter but that the mechanisms invoked were not effective in preventing CD-56 growth. Moreover, CD-56 did not appear to fully suppress plant defenses, even in late stages after infection. Our results provide new insight into the chronic pathogenesis of CD-56 and the comprehensive and complex inducible defense responses of ginseng root to I. robusta infection.

Comparative transcriptomics and transcriptional regulation analysis of enhanced laccase production induced by co-culture of Pleurotus eryngii var. ferulae with Rhodotorula mucilaginosa

The co-culturing of Pleurotus eryngii var. ferulae and Rhodotorula mucilaginosa was confirmed in our previous studies to be an efficient strategy to improve laccase production by submerged fermentation. To determine the possible regulation principles underlying this behaviour, comparative transcriptomic analysis was performed on P. eryngii var. ferulae to investigate the differential expression of genes in co-culture. RNA-seq analysis showed that genes concerning xenobiotic biodegradation and expenditure of energy were upregulated. However, genes related to oxidative stress were downregulated. In addition, the transcription levels of laccase isoenzymes were not consistent in the co-culture system: 3 laccase genes (lacc1, lacc2, lacc12) were upregulated, and 3 laccase genes (lacc4, lacc6, lacc9) were downregulated. The enhancement in laccase activity can be due to upregulation of a laccase heterodimer encoded by the genes lacc2 and ssPOXA3a (or ssPOXA3b), whose expression levels were increased by 459% and 769% (or 585% for ssPOXA3b) compared with those of a control, respectively. β-Carotene produced by R. mucilaginosa upregulated the transcription of lacc2 only. Combining these results with an analysis of cis-acting responsive elements indicated that four transcription factors (TFs) had potential regulatory effects on the transcription of laccase genes. It was supposed that TFa regulated lacc transcription by binding with methyl jasmonate and heat shock response elements. The expression of TFb, TFc, and TFd was regulated by β-carotene. However, β-carotene had no effect on TFa expression. These results provide a possible mechanism for the regulation of laccase gene transcription in the co-culture system and are also beneficial for the future intensification of fungal laccase production.

Pleurotus eryngii Genomes Reveal Evolution and Adaptation to the Gobi Desert Environment

Pleurotus eryngii (King Oyster) is one of the most highly prized edible mushrooms. Among the diverse varieties within P. eryngii, P. eryngii var. eryngii is the commonest one, with a worldwide distribution, while P. eryngii var. ferulae is only distributed in Europe and China, and is especially adapted to the Gobi Desert in Xinjiang Autonomous Region of China. However, little is known about the genome-wide pattern of evolution and adaptation to the divergent environments of P. eryngii. Here, we present the high-quality genome sequences of P. eryngii var. eryngii strain PEE81 originating from Europe and P. eryngii var. ferulae strain PEF12 originating from the Gobi Desert of China. The assembled genome sizes of PEE81 and PEF12 were 53.6 and 48.0 Mbp, respectively, which are larger than other reported genomes in the genus Pleurotus. We propose that the selective amplification of long terminal repeat (LTR) retrotransposons increases the genome size of the genus Pleurotus, and may play a key role in driving their rapid species diversification. Molecular clock analyses of five Pleurotus species, namely PEE81, PEF12, P. tuoliensis, P. ostreatus and P. cf. floridanus suggest that the divergence estimates of the genus Pleurotus over time scales ranged from ∼4 to ∼38 million years ago (Mya), and PEE81 and PEF12 diverged at ∼13 Mya. The whole genome resequencing of 33 geographically diverse strains of P. eryngii var. eryngii and var. ferulae was then performed and the genome variation among and within these two populations were investigated. Comparative analyses of these two populations detected several candidate genes related to stress responses and DNA repair that are putatively involved in adaptation to the Gobi Desert environment. These findings offer insights into genome evolution of the genus Pleurotus and provide valuable genomic resources for King Oyster mushroom breeding.

Manipulating the Expression of Small Secreted Protein 1 (Ssp1) Alters Patterns of Development and Metabolism in the White-Rot Fungus Pleurotus ostreatus

The function of small secreted proteins (SSPs) in saprotrophic fungi is, for the most part, unknown. The white-rot mushroom Pleurotus ostreatus produces considerable amounts of SSPs at the onset of secondary metabolism, during colony development, and in response to chemical compounds such as 5-hydroxymethylfurfural and aryl alcohols. Genetic manipulation of Ssp1, by knockdown (KDssp1) or overexpression (OEssp1), indicated that they are, in fact, involved in the regulation of the ligninolytic system. To elucidate their potential involvement in fungal development, quantitative secretome analysis was performed during the trophophase and the idiophase and at a transition point between the two growth phases. The mutations conferred a time shift in the secretion and expression patterns: OEssp1 preceded the entrance to idiophase and secondary metabolism, while KDssp1 was delayed. This was also correlated with expression patterns of selected genes. The KDssp1 colony aged at a slower pace, accompanied by a slower decline in biomass over time. In contrast, the OEssp1 strain exhibited severe lysis and aging of the colony at the same time point. These phenomena were accompanied by variations in yellow pigment production, characteristic of entrance of the wild type into idiophase. The pigment was produced earlier and in a larger amount in the OEssp1 strain and was absent from the KDssp1 strain. Furthermore, the dikaryon harboring OEssp1 exhibited a delay in the initiation of fruiting body formation as well as earlier aging. We propose that Ssp1 might function as a part of the fungal communication network and regulate the pattern of fungal development and metabolism in P. ostreatus IMPORTANCE Small secreted proteins (SSPs) are common in fungal saprotrophs, but their roles remain elusive. As such, they comprise part of a gene pool which may be involved in governing fungal lifestyles not limited to symbiosis and pathogenicity, in which they are commonly referred to as "effectors." We propose that Ssp1 in the white-rot fungus Pleurotus ostreatus regulates the transition from primary to secondary metabolism, development, aging, and fruiting body initiation. Our observations uncover a novel regulatory role of effector-like SSPs in a saprotroph, suggesting that they may act in fungal communication as well as in response to environmental cues. The presence of Ssp1 homologues in other fungal species supports a common potential role in environmental sensing and fungal development.

De novo transcriptome analysis of Pleurotus djamor to identify genes encoding CAZymes related to the decomposition of corn stalk lignocellulose

CAZymes play a very important role in the biotransformation of corn stalk biomass, which is an important resource for sustainable development. Pleurotus djamor can produce CAZymes related to the decomposition of corn stalk lignocellulose biomass in sole corn stalk substrate; however, little is known about their encoding genes. In order to identify CAZymes encoding genes, RNA high-throughput sequencing of P. djamor was performed in this study. The results showed that a core set of 70 upregulated genes encoding putative CAZymes were revealed. They encode 19 kinds of CAZymes in total, of which there are 4 EGLs, 8 CBHs, 5 BGLs, and 12 LPMOs related to cellulose degradation, 8 XYNs, 1 XYL, 2 AGUs, 3 ABFs, 2 AGLs, and 2 AXEs related to hemicellulose degradation, and 5 LACCs, 2 MnPs, 5 VPs, 3 CDHs, 1 AAO, 1 GOX, 1 AOX, 2 GAOXs, and 3 GLOXs related to lignin degradation. This variety suggests that CAZymes may play a very important role in decomposing the lignocellulose biomass of corn stalk. This is the first study to report the de novo transcriptome sequencing of P. djamor, which will produce a dataset of genes encoding CAZymes, thereby laying the foundation to elucidate the degradation mechanism of corn stalk biomass and boost the biotransformation of corn stalk biomass resources.

Mapping the Secretome and Its N-Linked Glycosylation of Pleurotus eryngii and Pleurotus ostreatus Grown on Hemp Stalks

Our previous research showed that Pleurotus eryngii and Pleurotus ostreatus were effective fungi for pretreatment of industrial hemp stalks to improve enzymatic saccharification. The secretomes of these two fungi were analyzed to search for the effective enzyme cocktails degrading hemp lignin during the pretreatment process. In total, 169 and 155 proteins were identified in Pleurotus eryngii and Pleurotus ostreatus, respectively, and 50% of the proteins involved in lignocellulose degradation were CAZymes. Because most of the extracellular proteins secreted by fungi are glycosylated proteins, the N-linked glycosylation of enzymes could be mapped. In total, 27 and 24 N-glycosylated peptides were detected in Pleurotus eryngii and Pleurotus ostreatus secretomes, respectively. N-Glycosylated peptides of laccase, GH92, exoglucanase, phenol oxidase, α-galactosidase, carboxylic ester hydrolase, and pectin lyase were identified. Deglycosylation could decrease enzymatic saccharification of hemp stalks. The activities of laccase, α-galactosidase, and phenol oxidase and the thermal stability of laccase were reduced after deglycosylation.

Biochemical characterization of a key laccase-like multicopper oxidase of artificially cultivable Morchella importuna provides insights into plant-litter decomposition

The aim of this study is to determine the key laccase-encoding gene in the life cycle of Morchella importuna SCYDJ1-A1, and to characterize the biochemical properties of the laccase. Two laccase-like multicopper oxidase (LMCO) genes were identified in the genome of M. importuna SCYDJ1-A1 as putative laccase-encoding genes. The two genes, belonging to Auxiliary Activity family 1 subfamily 3, were named as MiLacA and MiLacB. Phylogenetic analysis of deduced amino acid sequences showed that MiLacA is closest to a LMCO of M. importuna 22J1, while MiLacB had low similarity with known Morchella LMCOs. Real-time quantitative PCR results showed that MiLacA was expressed at much higher levels than MiLacB throughout the entire course of artificial cultivation. MiLacA was overexpressed in Pichia pastoris as a recombinant protein. Biochemical characterization of the purified enzyme showed that MiLacA simultaneously possessed laccase and polyphenol-oxidase activities. MiLacA could be strongly inhibited by Fe²⁺, which is unusual. The optimum pH was four and optimum temperature was 60 °C. The enzyme retained over 74% of the laccase activity after 16-h incubation at 60 °C, which means that its thermostability is at the forefront among the currently known laccases. Our findings may help to elucidate how the laccase of M. importuna is involved in decaying lignin in plant litter, and could also provide a candidate thermostable laccase for potential industrial application.

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

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

Comparative genomics of 40 edible and medicinal mushrooms provide an insight into the evolution of lignocellulose decomposition mechanisms

Diversity comparison and phylogenetic analyses of carbohydrate-active enzymes (CAZymes), auxiliary activities (AAs) and cytochromes P450 among 40 fungi, which are based on different nutritional pathways, help clarify and explain their divergence and improvement of various life-styles. Molecular clock analyses allow us to understand the evolutionary and developmental rules in decomposition gene families. Our results suggested that fungi in different ecological types acquired an obvious preference on specific decomposing gene families during evolutionary selection. White rot and litter saprotrophic fungi possessed more complete types of varied degradation gene families and were superior in quantities. With evolution and development of lignocellulose decomposition mechanism, certain families (like CBM1, GH6, GH7, GH10, and CYP53) disappeared in brown rot fungi and symbiotic fungi. In addition, the earlier time of phylogenetic divergence determined the more integrated and larger decomposition families. And various gains and losses in gene quantity of varied decomposition families led in particularly phylogenetic clades or nodes, then accelerated in forming varied ecotypes of species.

Molecular characterization and overexpression of mnp6 and vp3 from Pleurotus ostreatus revealed their involvement in biodegradation of cotton stalk lignin

Fungal secretory heme peroxidase (Class II POD) plays a significant role in biomass conversion due to its lignin-degrading activity. In this study, genome-wide identification and bioinformatics were performed to analyze P leurotus ostreatus peroxidases (PoPODs). A total of six manganese peroxidases (MnPs) and three versatile peroxidases (VPs) were obtained. Bioinformatics analysis and qRT-PCR showed that P. ostreatus mnp6 (Pomnp6) and P. ostreatus vp3 (Povp3) could be involved in lignin degradation. Both Pomnp6 and Povp3 transgenetic fungi showed significantly increased lignin degradation of cotton stalks. 1H-NMR revealed that Pomnp6 and Povp3 may preferentially degrade S-lignin in cotton stalks and mainly break β-O-4' bond linkages and hydroxyl. These results support the possible utility of Pomnp6 and Povp3 in natural straw resources and development of sustainable energy.

Spatiotemporal Dynamics of Molecular Messaging in Bacterial Co-Cultures Studied by Multimodal Chemical Imaging

Microbial community behavior is coupled to a set of genetically-regulated chemical signals that correlate with cell density - the quorum sensing (QS) system - and there is growing appreciation that the QS-regulated behavior of bacteria is chemically, spatially, and temporally complex. In addition, while it has been known for some time that different species use different QS networks, we are beginning to appreciate that different strains of the same bacterial species also differ in their QS networks. Here we combine mass spectrometric imaging (MSI) and confocal Raman microscopy (CRM) approaches to investigate co-cultures involving different strains (FRD1 and PAO1C) of the same species (Pseudomonas aeruginosa) as well as those involving different species (P. aeruginosa and E. coli). Combining MSI and CRM makes it possible to supersede the limits imposed by individual imaging approaches and enables the spatial mapping of individual bacterial species and their microbial products within a mixed bacterial community growing in situ on surfaces. MSI is used to delineate the secretion of a specific rhamnolipid surfactant as well as alkyl quinolone (AQ) messengers between FRD1 and PAO1C strains of P. aeruginosa, showing that the spatial distribution and production rate of AQ messengers in PAO1C far outstrips that of FRD1. In the case of multiple species, CRM is used to show that the prolific secretion of AQs by the PAO1C strain of P. aeruginosa is used to mediate its interaction with co-cultured E. coli.

Carbon sources and XlnR-dependent transcriptional landscape of CAZymes in the industrial fungus Talaromyces versatilis: when exception seems to be the rule

Background

Research on filamentous fungi emphasized the remarkable redundancy in genes encoding hydrolytic enzymes, the similarities but also the large differences in their expression, especially through the role of the XlnR/XYR1 transcriptional activator. The purpose of this study was to evaluate the specificities of the industrial fungus Talaromyces versatilis, getting clues into the role of XlnR and the importance of glucose repression at the transcriptional level, to provide further levers for cocktail production.

Results

By studying a set of 62 redundant genes representative of several categories of enzymes, our results underlined the huge plasticity of transcriptional responses when changing nutritional status. As a general trend, the more heterogeneous the substrate, the more efficient to trigger activation. Genetic modifications of xlnR led to significant reorganisation of transcriptional patterns. Just a minimal set of genes actually fitted in a simplistic model of regulation by a transcriptional activator, and this under specific substrates. On the contrary, the diversity of xlnR⁺ versus ΔxlnR responses illustrated the existence of complex and unpredicted patterns of co-regulated genes that were highly dependent on the culture condition, even between genes that encode members of a functional category of enzymes. They notably revealed a dual, substrate-dependant repressor-activator role of XlnR, with counter-intuitive transcripts regulations that targeted specific genes. About glucose, it appeared as a formal repressive sugar as we observed a massive repression of most genes upon glucose addition to the mycelium grown on wheat straw. However, we also noticed a positive role of this sugar on the basal expression of a few genes, (notably those encoding cellulases), showing again the strong dependence of these regulatory mechanisms upon promoter and nutritional contexts.

Conclusions

The diversity of transcriptional patterns appeared to be the rule, while common and stable behaviour, both within gene families and with fungal literature, the exception. The setup of a new biotechnological process to reach optimized, if not customized expression patterns of enzymes, hence appeared tricky just relying on published data that can lead, in the best scenario, to approximate trends. We instead encourage preliminary experimental assays, carried out in the context of interest to reassess gene responses, as a mandatory step before thinking in (genetic) strategies for the improvement of enzyme production in fungi.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1062-8) contains supplementary material, which is available to authorized users.

Integrative visual omics of the white-rot fungus Polyporus brumalis exposes the biotechnological potential of its oxidative enzymes for delignifying raw plant biomass

BACKGROUND

Plant biomass conversion for green chemistry and bio-energy is a current challenge for a modern sustainable bioeconomy. The complex polyaromatic lignin polymers in raw biomass feedstocks (i.e., agriculture and forestry by-products) are major obstacles for biomass conversions. White-rot fungi are wood decayers able to degrade all polymers from lignocellulosic biomass including cellulose, hemicelluloses, and lignin. The white-rot fungus Polyporus brumalis efficiently breaks down lignin and is regarded as having a high potential for the initial treatment of plant biomass in its conversion to bio-energy. Here, we describe the extraordinary ability of P. brumalis for lignin degradation using its enzymatic arsenal to break down wheat straw, a lignocellulosic substrate that is considered as a biomass feedstock worldwide.

RESULTS

We performed integrative multi-omics analyses by combining data from the fungal genome, transcriptomes, and secretomes. We found that the fungus possessed an unexpectedly large set of genes coding for Class II peroxidases involved in lignin degradation (19 genes) and GMC oxidoreductases/dehydrogenases involved in generating the hydrogen peroxide required for lignin peroxidase activity and promoting redox cycling of the fungal enzymes involved in oxidative cleavage of lignocellulose polymers (36 genes). The examination of interrelated multi-omics patterns revealed that eleven Class II Peroxidases were secreted by the fungus during fermentation and eight of them where tightly co-regulated with redox cycling enzymatic partners.

CONCLUSION

As a peculiar feature of P. brumalis, we observed gene family extension, up-regulation and secretion of an abundant set of versatile peroxidases and manganese peroxidases, compared with other Polyporales species. The orchestrated secretion of an abundant set of these delignifying enzymes and redox cycling enzymatic partners could contribute to the delignification capabilities of the fungus. Our findings highlight the diversity of wood decay mechanisms present in Polyporales and the potentiality of further exploring this taxonomic order for enzymatic functions of biotechnological interest.

Genome-Wide Characterization and Expression Analyses of Pleurotus ostreatus MYB Transcription Factors during Developmental Stages and under Heat Stress Based on de novo Sequenced Genome

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.