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

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

Structural insights into latency of the metallopeptidase ulilysin (lysargiNase) and its unexpected inhibition by a sulfonyl-fluoride inhibitor of serine peptidases

Peptidases are regulated by latency and inhibitors, as well as compatibilization and cofactors. Ulilysin from Methanosarcina acetivorans, also called lysargiNase, is an archaeal metallopeptidase (MP) that is biosynthesized as a zymogen with a 60-residue N-terminal prosegment (PS). In the presence of calcium, it self-activates to yield the mature enzyme, which specifically cleaves before basic residues and thus complements trypsin in proteomics workflows. Here, we obtained a low-resolution crystal structure of proulilysin, in which 28 protomers arranged as 14 dimers form a continuous double helix of 544 Å pitch that parallels cell axis b of the crystal. The PS includes two α-helices and obstructs the active-site cleft of the catalytic domain (CD) by traversing it in the opposite orientation of a substrate, and a cysteine blocks the catalytic zinc according to a "cysteine-switch mechanism". Moreover, the PS interacts through its first helix with an "S-loop" of the CD, which acts as an "activation segment" that lacks one of two essential calcium cations. Upon PS removal during maturation, the S-loop adopts its competent conformation and binds the second calcium ion. Next, we found that in addition to general MP inhibitors, ulilysin was competitively and reversibly inhibited by 4-(2-aminoethyl)benzenesulfonyl fluoride (AEBSF; K_i = 4 μM). This is a compound that normally forms an irreversible covalent complex with serine peptidases but does not inhibit MPs. A high-resolution crystal structure of the complex revealed that the inhibitor penetrates the specificity pocket of ulilysin. A primary amine of the inhibitor salt-bridges an aspartate at the pocket bottom, thus mimicking the basic side chain of substrates. In contrast, the sulfonyl fluoride warhead is not involved and the catalytic zinc ion is freely accessible. Thus, the usage of inhibitor cocktails of peptidases, which typically contain AEBSF at ∼25-fold higher concentrations than the determined K_i, should be avoided when working with ulilysin. Finally, the structure of the complex, which occurred as a crystallographic dimer recurring in previous mature ulilysin structures, unveiled an N-terminal product fragment that delineated the non-primed side of the cleft. These results complement prior structures of ulilysin with primed-side product fragments and inhibitors.

Comparative transcriptome analysis revealed candidate genes involved in fruiting body development and sporulation in Ganoderma lucidum

Ganoderma lucidum is an edible mushroom highly regarded in the traditional Chinese medicine. To better understand the molecular mechanisms underlying fruiting body development in G. lucidum, transcriptome analysis based on RNA sequencing was carried out on different developmental stages: mycelium (G1); primordium (G2); young fruiting body (G3); mature fruiting body (G4); fruiting body in post-sporulation stage (G5). In total, 26,137 unigenes with an average length of 1078 bp were de novo assembled. Functional annotation of transcriptomes matched 72.49% of the unigenes to known proteins available in at least one database. Differentially expressed genes (DEGs) were identified between the evaluated stages: 3135 DEGs in G1 versus G2; 120 in G2 versus G3; 3919 in G3 versus G4; and 1012 in G4 versus G5. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs identified in G1 versus G2 revealed that, in addition to global and overview maps, enriched pathways were related to amino acid metabolism and carbohydrate metabolism. In contrast, DEGs identified in G2 versus G3 were mainly assigned to the category of metabolism of amino acids and their derivatives, comprising mostly upregulated unigenes. In addition, highly expressed unigenes associated with the transition between different developmental stages were identified, including those encoding hydrophobins, cytochrome P450s, extracellular proteases, and several transcription factors. Meanwhile, highly expressed unigenes related to meiosis such as DMC1, MSH4, HOP1, and Mek1 were also analyzed. Our study provides important insights into the molecular mechanisms underlying fruiting body development and sporulation in G. lucidum.

Extracellular proteases and laccases produced by Pleurotus ostreatus PoB: the effects of proteases on laccase activity

Laccases are enzymes produced by plants and white rot fungi, such as Pleurotus ostreatus, with industrial applications. Fungal laccases have been widely studied, and investigations, such as those involving recombinant DNA technology or adding inducers, have been made to increase laccase production. On the other hand, it has been proposed that extracellular proteases could decrease laccase activity when both types of enzymes are produced by P. ostreatus. The aim of this work was to evaluate the effects of proteases on the activity of extracellular laccases produced by P. ostreatus PoB in submerged culture. Results showed that P. ostreatus PoB produced alkaline, acidic, and neutral proteases. Protease activity was quantified, and the highest activity at alkaline pH (9.0) was 5.63 IU/L (192 h), that at acidic pH (2.0) was 3.38 IU/L (192 h), and that at neutral pH (7.0) was 6.20 IU/L (312 h). The protease activity decreased in the presence of different protease inhibitors, as phenylmethylsulfonyl fluoride (PMSF), EDTA, pepstatin A, and a cocktail of protease inhibitors. Laccase activity was determined in cultures with and without protease inhibitors. In the control culture (without inhibitor), the highest laccase specific activity was 99.88 IU/mg protein. In cultures with PMSF, pepstatin A, or a cocktail of protease inhibitors, laccase activity increased by approximately 1.35-fold (138 IU/mg protein) with respect to the control culture. The inhibitor EDTA did not produce a positive effect on extracellular laccase activity. These results suggest that laccase activity is affected by the actions of acidic and neutral extracellular proteases.

Peptidyl-Lys metalloendopeptidase (Lys-N) purified from dry fruit of Grifola frondosa demonstrates "mirror"digestion property with lysyl endopeptidase (Lys-C)

RATIONALE

Lys-N, also known as lysine-specific metalloendopeptidase, functions as the "sister" enzyme of lysyl endopeptidase (Lys-C) in proteomic research. Its digestion specificity at the N-terminal lysine residue makes it a very useful tool in proteomics analysis, especially in mass spectrometry (MS)-based de novo sequencing of proteins.

METHODS

Here we present a complete production process of highly purified Lys-N from dry fruit of Grifola frondosa (maitake mushroom). The purification process includes one step of microfiltration plus one step of UF/DF (ultrafiltration used in tandem with a diafiltration method) recovery and four steps of chromatographic purification.

RESULTS

The overall yield of the process was approximately 6.7 mg Lys-N protein/kg dry fruit of G. frondosa. The assay data demonstrated that the purified Lys-N exhibited high enzymatic activity and specificity.

CONCLUSIONS

The novel production process provides for the first time the extraction of Lys-N from dry fruit of G. frondosa. The process is also stable and scalable, and provides an economic way of producing the enzyme in large quantities for MS-based proteomics and other biological studies.

Genome Survey and Transcriptome Analysis on Mycelia and Primordia of Agaricus blazei

Agaricus blazei, a type of edible straw-rotting mushroom with somewhat sweet taste and fragrance of almonds, has attracted considerable scientific and practical attention. High-throughput Illumina PE150 and PacBio RSII platform were employed to generate a genomic sequence. De novo assembly generated 36 contigs with 38,686,133 bp in size, containing 10,119 putative predicted genes. Additionally, we also studied transcriptional regulation of the mycelia and the primordia for exploration of genes involved in fruiting body formation. Expression profiling analysis revealed that 2,164 genes were upregulated in mycelia and 1,557 in primordia. Functional enrichment showed that differentially expressed genes associated with response to stress, ribosome biogenesis, arginine biosynthesis, and steroid biosynthesis pathway were more active in fruiting body. The genome and transcriptome analysis of A. blazei provide valuable sequence resources and contribute to our understanding of genes related to the biosynthesis pathway of polysaccharide and benzaldehyde, as well as the fruiting body formation.

Proteases of Wood Rot Fungi with Emphasis on the Genus Pleurotus

Proteases are present in all living organisms and they play an important role in physiological conditions. Cell growth and death, blood clotting, and immune defense are all examples of the importance of proteases in maintaining homeostasis. There is growing interest in proteases due to their use for industrial purposes. The search for proteases with specific characteristics is designed to reduce production costs and to find suitable properties for certain industrial sectors, as well as good producing organisms. Ninety percent of commercialized proteases are obtained from microbial sources and proteases from macromycetes have recently gained prominence in the search for new enzymes with specific characteristics. The production of proteases from saprophytic basidiomycetes has led to the identification of various classes of proteases. The genus Pleurotus has been extensively studied because of its ligninolytic enzymes. The characteristics of this genus are easy cultivation techniques, high yield, low nutrient requirements, and excellent adaptation. There are few studies in the literature about proteases of Pleurotus spp. This review gathers together information about proteases, especially those derived from basidiomycetes, and aims at stimulating further research about fungal proteases because of their physiological importance and their application in various industries such as biotechnology and medicine.

Reprolysin metalloproteases from Ixodes persulcatus, Rhipicephalus sanguineus and Rhipicephalus microplus ticks

Metalloproteases (MPs) have been considered essential for blood feeding and other physiological functions in several hematophagous animals, including ticks. We report the characterization of MP sequences of three important ticks from Asia, Africa and America: Ixodes persulcatus (Ip-MPs), Rhipicephalus sanguineus (Rs-MPs) and R. microplus (BrRm-MPs). Amino acid sequence identity between R. microplus and R. sanguineus MPs ranged from 76 to 100 %, and identities among I. persulcatus, I. ricinus and I. scapularis MP sequences ranged from 88 to 97 %. This high sequence identity and typical functional motifs show that all sequences are MPs. The presence of a zinc binding site, a Met-turn and cysteine rich domain at the C-terminal region indicates that these proteins belong to the reproplysin family of MPs. Differences in amino acid sequences of BrRm-MP1, BrRm-MP2, BrRm-MP4 and BrRm-MP5 (from Porto Alegre strain ticks) were 6, 2, 7 and 5 %, respectively, when compared with sequences deposited in GenBank for the same genes from other R. microplus isolates. Analyses of MPs predicted that they have various highly antigenic regions. Semi-quantitative RT-PCR analysis revealed the presence of transcripts in salivary glands of partially and fully fed female ticks. None of these transcripts were observed in males (except BrRm-MP4) and eggs. These enzymes may be functional components required during tick feeding to manipulate host defenses and support tick hematophagy.

Sequencing and comparative analysis of the straw mushroom (Volvariella volvacea) genome

Volvariella volvacea, the edible straw mushroom, is a highly nutritious food source that is widely cultivated on a commercial scale in many parts of Asia using agricultural wastes (rice straw, cotton wastes) as growth substrates. However, developments in V. volvacea cultivation have been limited due to a low biological efficiency (i.e. conversion of growth substrate to mushroom fruit bodies), sensitivity to low temperatures, and an unclear sexuality pattern that has restricted the breeding of improved strains. We have now sequenced the genome of V. volvacea and assembled it into 62 scaffolds with a total genome size of 35.7 megabases (Mb), containing 11,084 predicted gene models. Comparative analyses were performed with the model species in basidiomycete on mating type system, carbohydrate active enzymes, and fungal oxidative lignin enzymes. We also studied transcriptional regulation of the response to low temperature (4°C). We found that the genome of V. volvacea has many genes that code for enzymes, which are involved in the degradation of cellulose, hemicellulose, and pectin. The molecular genetics of the mating type system in V. volvacea was also found to be similar to the bipolar system in basidiomycetes, suggesting that it is secondary homothallism. Sensitivity to low temperatures could be due to the lack of the initiation of the biosynthesis of unsaturated fatty acids, trehalose and glycogen biosyntheses in this mushroom. Genome sequencing of V. volvacea has improved our understanding of the biological characteristics related to the degradation of the cultivating compost consisting of agricultural waste, the sexual reproduction mechanism, and the sensitivity to low temperatures at the molecular level which in turn will enable us to increase the industrial production of this mushroom.

Transcriptome and proteome exploration to provide a resource for the study of Agrocybe aegerita

BACKGROUND

Agrocybe aegerita, the black poplar mushroom, has been highly valued as a functional food for its medicinal and nutritional benefits. Several bioactive extracts from A. aegerita have been found to exhibit antitumor and antioxidant activities. However, limited genetic resources for A. aegerita have hindered exploration of this species.

METHODOLOGY/PRINCIPAL FINDINGS

To facilitate the research on A. aegerita, we established a deep survey of the transcriptome and proteome of this mushroom. We applied high-throughput sequencing technology (Illumina) to sequence A. aegerita transcriptomes from mycelium and fruiting body. The raw clean reads were de novo assembled into a total of 36,134 expressed sequences tags (ESTs) with an average length of 663 bp. These ESTs were annotated and classified according to Gene Ontology (GO), Clusters of Orthologous Groups (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways. Gene expression profile analysis showed that 18,474 ESTs were differentially expressed, with 10,131 up-regulated in mycelium and 8,343 up-regulated in fruiting body. Putative genes involved in polysaccharide and steroid biosynthesis were identified from A. aegerita transcriptome, and these genes were differentially expressed at the two stages of A. aegerita. Based on one-dimensional gel electrophoresis (1-DGE) coupled with electrospray ionization liquid chromatography tandem MS (LC-ESI-MS/MS), we identified a total of 309 non-redundant proteins. And many metabolic enzymes involved in glycolysis were identified in the protein database.

CONCLUSIONS/SIGNIFICANCE

This is the first study on transcriptome and proteome analyses of A. aegerita. The data in this study serve as a resource of A. aegerita transcripts and proteins, and offer clues to the applications of this mushroom in nutrition, pharmacy and industry.

Opposing developmental functions of Agrocybe aegerita galectin (AAL) during mycelia differentiation

Mycelia of basidiomycetes differentiating into fruiting body is a controlled developmental process, however the underlying molecular mechanism remains unknown. In previous work, a novel fungal Agrocybe aegerita galectin (AAL) was isolated from A. aegerita in our laboratory. AAL was shown to promote mycelial differentiation in A. aegerita and Auricularia polytricha, indicating that AAL might function as a conserved fruiting initiator during basidiomycete mycelia development. In the current work, we investigate the role of AAL in mycelia differentiation and fruiting body formation. First, the expression and localization of AAL in mycelia, primordium and fruiting body were assessed by Western blotting and immunohistochemistry. AAL was found to be ubiquitously expressed in the primordium and fruiting body but not in the mycelia. AAL facilitated mycelia congregation and promoted fruiting body production when AAL was applied on mycelia. At the same time, when AAL was spread on potato dextrose agar (PDA) medium prior to mycelia inoculation, mycelia exhibited slowed growth rates, resulting in mycelia cords formation and inhibition of fruiting body formation. The 5' regulatory sequence of aal was cloned by 'genome walking'. Here, we show that aal lack introns in the coding region and the upstream 740 bp sequence was characterized by the existence of core promoter elements, which included: two CCAAT boxes (-535/-280), a GC box (-145), a TATA box (-30) and a fungal leader intron within the 5' UTR. The identification of regulatory expression elements may provide an explanation to the stage-specific and high-level expression of aal during fruiting development.

Substrate specificity of a metalloprotease of the pappalysin family revealed by an inhibitor and a product complex

Human pappalysin-1 is a multi-domain metalloprotease engaged in the homeostasis of insulin-like growth factors and the founding member of the pappalysin family within the metzincin clan of metalloproteases. We have recently identified an archaeal relative, ulilysin, encompassing only the protease domain. It is a 262-residue active protease with a novel 3D structure with two subdomains separated by an active-site cleft. Despite negligible overall sequence similarity, noticeable similarity is found with other metzincin prototypes, adamalysins/ADAMs and matrix metalloproteinases. Ulilysin has been crystallised in a product complex with an arginine-valine dipeptide occupying the active-site S(1') and S(2') positions and in a complex with the broad-spectrum hydroxamic acid-based metalloprotease inhibitor, batimastat. This molecule inhibits mature ulilysin with an IC(50) value of 61 microM under the conditions assayed. The binding of batimastat to ulilysin evokes binding to vertebrate matrix metalloproteases but is much weaker. These data give insight into substrate specificity and mechanism of action and inhibition of the novel pappalysin family.

Molecular Analysis of Ulilysin, the Structural Prototype of a New Family of Metzincin Metalloproteases

The metzincin clan encompasses several families of zinc-dependent metalloproteases with proven function both in physiology and pathology. They act either as broad spectrum protein degraders or as sheddases, operating through limited proteolysis. Among the structurally uncharacterized metzincin families are the pappalysins, of which the most thoroughly studied member is human pregnancy-associated plasma protein A (PAPP-A), a heavily glycosylated 170-kDa multidomain protein specifically cleaving insulin-like growth factor (IGF)-binding proteins (IGFBPs). Proulilysin is a 38-kDa archaeal protein that shares sequence similarity with PAPP-A but encompasses only the pro-domain and the catalytic domain. It undergoes calcium-mediated autolytic activation, and the mature protein adopts a three-dimensional structure with two subdomains separated by an active site cleft containing the catalytic zinc ion. This structure is reminiscent of human members of the adamalysin/ADAMs (a disintegrin and a metalloprotease) family of metzincins. A bound dipeptide yields information on the substrate specificity of ulilysin, which specifically hydrolyzes IGFBP-2 to -6, insulin, and extracellular matrix proteins but not IGFBP-1 or IGF-II. Accordingly, ulilysin has higher proteolytic efficiency and a broader substrate specificity than human PAPP-A. The structure of ulilysin represents a prototype for the catalytic domain of pappalysins.

Genetic linkage map and expression analysis of genes expressed in the lamellae of the edible basidiomycete Pleurotus ostreatus

Pleurotus ostreatus is an industrially cultivated basidiomycete with nutritional and environmental applications. Its genome contains 35 Mbp organized in 11 chromosomes. There is currently available a genetic linkage map based predominantly on anonymous molecular markers complemented with the mapping of QTLs controlling growth rate and industrial productivity. To increase the saturation of the existing linkage maps, we have identified and mapped 82 genes expressed in the lamellae. Their manual annotation revealed that 34.1% of the lamellae-expressed and 71.5% of the lamellae-specific genes correspond to previously unknown sequences or to hypothetical proteins without a clearly established function. Furthermore, the expression pattern of some genes provides an experimental basis for studying gene regulation during the change from vegetative to reproductive growth. Finally, the identification of various differentially regulated genes involved in protein metabolism suggests the relevance of these processes in fruit body formation and maturation.

Purification and characterization of fibrinolytic enzyme from cultured mycelia of Armillaria mellea

A fibrinolytic enzyme was purified from the cultured mycelia of Armillaria mellea by ion-exchange chromatography followed by gel filtration, and was designated A. mellea metalloprotease (AMMP). The purification protocol resulted in a 627-fold purification of the enzyme, with a final yield of 6.05%. The apparent molecular mass of the purified enzyme was estimated to be 21kDa by SDS-PAGE, fibrin-zymography and gel filtration chromatography, which revealed a monomeric form of the enzyme. The optimal reaction pH value and temperature were, pH 6.0, and 33 degrees C, respectively. This protease effectively hydrolyzed fibrinogen, preferentially digesting the Aalpha-chain over the Bbeta- and r-chains. Enzyme activity was inhibited by Cu(2+) and Co(2+), but enhanced by the addition of Ca(2+) and Mg(2+) ions. Furthermore, AMMP activity was potently inhibited by EDTA, and was found to exhibit a higher specificity for the substrate S-2586 for chymotrypsin, indicating that the enzyme is a chymotrypsin-like metalloprotease. The first 24 amino acid residues of the N-terminal sequence were MFSLSSRFFLYTLCL SAVAVSAAP, which is extremely similar to the 24 amino acid residues of the N-terminal sequence of the fruiting body of A. mellea. These data suggest that the fibrinolytic enzyme AMMP, obtained from the A. mellea exhibits a profound fibrinolytic activity. The mycelia of A. mellea may thus represent a potential source of new therapeutic agents to treat thrombosis.