Expression of the urease gene of Agaricus bisporus: a tool for studying fruit body formation and post-harvest development

The First Whole Genome Sequencing of Agaricus bitorquis and Its Metabolite Profiling

Agaricus bitorquis, an emerging wild mushroom with remarkable biological activities and a distinctive oversized mushroom shape, has gained increasing attention in recent years. Despite its status as an important resource of wild edible fungi, knowledge about this mushroom is still limited. In this study, we used the Illumina NovaSeq and Nanopore PromethION platforms to sequence, de novo assemble, and annotate the whole genome and mitochondrial genome (mitogenome) of the A. bitorquis strain BH01 isolated from Bosten Lake, Xinjiang Province, China. Using the genome-based biological information, we identified candidate genes associated with mating type and carbohydrate-active enzymes in A. bitorquis. Cluster analysis based on P450 of basidiomycetes revealed the types of P450 members of A. bitorquis. Comparative genomic, mitogenomic, and phylogenetic analyses were also performed, revealing interspecific differences and evolutionary features of A. bitorquis and A. bisporus. In addition, the molecular network of metabolites was investigated, highlighting differences in the chemical composition and content of the fruiting bodies of A. bitorquis and A. bisporus. The genome sequencing provides a comprehensive understanding and knowledge of A. bitorquis and the genus Agaricus mushrooms. This work provides valuable insights into the potential for artificial cultivation and molecular breeding of A. bitorquis, which will facilitate the development of A. bitorquis in the field of edible mushrooms and functional food manufacture.

Identification of bHLH family genes in Agaricus bisporus and transcriptional regulation of arginine catabolism-related genes by AbbHLH1 after harvest

Basic helix-loop-helix (bHLH) transcription factors (TFs) are widely distributed in eukaryotes and play an important role in biological growth and development. The identification and functional analyses of bHLH genes/proteins in edible mushrooms (Agaricus bisporus) have yet to be reported. In the present study, we identified 10 putative bHLH members carrying the conserved bHLH domains. Phylogenetic analyses revealed that the 10 AbbHLHs were the closest to sequences of species belonging to 7 different fungal subgroups, which was supported by loop length, intron patterns, and key amino acid residues. The substantial increase after harvest and continuously elevated expression of AbbHLH1 during the development until the disruption of mushroom velum, and the preferential expression in cap and gill tissues suggest the important function of AbbHLH1 in postharvest development of A. bisporus. The relationship of arginine catabolism-related genes with the early stage of postharvest continuing development also was revealed by expression determination. Subcellular localization showed that AbbHLH1 could be localized in nucleus. Importantly, the electrophoretic mobility shift and dual-luciferase reporter assays showed that AbbHLH1 activated the promoters of AbOAT, AbSPDS, and AbSAMDC and suppressed the expression of AbARG, AbUREA, and AbODC, probably for the modulation of arginine catabolism and thus control of postharvest mushroom development. Taken together, the available data provide valuable functional insight into the role of AbbHLH proteins in postharvest mushrooms.

Comparative transcriptome profiling of Termitomyces sp. between monocultures in vitro and link-stipe of fungus-combs in situ

The edible mushroom Termitomyces is an agaric-type basidiomycete fungus that has a symbiotic relationship with fungus-growing termites. An understanding of the detailed development mechanisms underlying the adaptive responses of Termitomyces sp. to their growing environment is lacking. Here, we compared the transcriptome sequences of different Termitomyces sp. samples and link-stipe grown on fungus combs in situ and monocultured in vitro. The assembled reads generated 8052 unigenes. The expression profiles were highly different for 2556 differentially expressed genes (DEGs) of the treated samples, where the expression of 1312 and 1244 DEGs was upregulated in the Mycelium and link-stipe groups respectively. Functional classification of the DEGs based on both Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed an expected shift in fungal gene expression, where stress response genes whose expression was upregulated in link-stipe may adaptively be involved in cell wall hydrolysis and fusion, pathogenesis, oxidation-reduction, transporter efflux, transposon efflux and self/non-self-recognition. Urease has implications in the expression of genes involved in the nitrogen metabolism pathway, and its expression could be controlled by low-level nitrogen fixation of fungus combs. In addition, the expression patterns of eleven select genes on the basis of qRT-PCR were consistent with their changes in transcript abundance, as revealed by RNA sequencing. Taken together, these findings may be useful for enriching the knowledge concerning the Termitomyces adaptive response to in situ fungus combs compared with the response of monocultures in vitro.

Proteomic investigation of interhyphal interactions between strains of Agaricus bisporus

Hyphae of filamentous fungi undergo polar extension, bifurcation and hyphal fusion to form reticulating networks of mycelia. Hyphal fusion or anastomosis, a ubiquitous process among filamentous fungi, is a vital strategy for how fungi expand over their substrate and interact with or recognise self- and non-self hyphae of neighbouring mycelia in their environment. Morphological and genetic characterisation of anastomosis has been studied in many model fungal species, but little is known of the direct proteomic response of two interacting fungal isolates. Agaricus bisporus, the most widely cultivated edible mushroom crop worldwide, was used as an in vitro model to profile the proteomes of interacting cultures. The globally cultivated strain (A15) was paired with two distinct strains; a commercial hybrid strain and a wild isolate strain. Each co-culture presented a different interaction ranging from complete vegetative compatibility (self), lack of interactions, and antagonistic interactions. These incompatible strains are the focus of research into disease-resistance in commercial crops as the spread of intracellular pathogens, namely mycoviruses, is limited by the lack of interhyphal anastomosis. Unique proteomic responses were detected between all co-cultures. An array of cell wall modifying enzymes, plus fungal growth and morphogenesis proteins were found in significantly (P < 0.05) altered abundances. Nitrogen metabolism dominated in the intracellular proteome, with evidence of nitrogen starvation between competing, non-compatible cultures. Changes in key enzymes of A. bisporus morphogenesis were observed, particularly via increased abundance of glucanosyltransferase in competing interactions and certain chitinases in vegetative compatible interactions only. Carbohydrate-active enzyme arsenals are expanded in antagonistic interactions in A. bisporus. Pathways involved in carbohydrate metabolism and genetic information processing were higher in interacting cultures, most notably during self-recognition. New insights into the differential response of interacting strains of A. bisporus will enhance our understanding of potential barriers to viral transmission through vegetative incompatibility. Our results suggest that a differential proteomic response occurs between A. bisporus at strain-level and findings from this work may guide future proteomic investigation of fungal anastomosis.

Integrated Metabolomics and Transcriptomics Unravel the Metabolic Pathway Variations for Different Sized Beech Mushrooms

Beech mushrooms (Hypsizygus marmoreus) are largely relished for their characteristic earthy flavor, chewy-texture, and gustatory and nutritional properties in East Asian societies. Intriguingly, the aforementioned properties of beech mushroom can be subsumed under its elusive metabolome and subtle transcriptome regulating its various stages of growth and development. Herein, we carried out an integrated metabolomic and transcriptomic profiling for different sized beech mushrooms across spatial components (cap and stipe) to delineate their signature pathways. We observed that metabolite profiles and differentially expressed gene (DEGs) displayed marked synergy for specific signature pathways according to mushroom sizes. Notably, the amino acid, nucleotide, and terpenoid metabolism-related metabolites and genes were more abundant in small-sized mushrooms. On the other hand, the relative levels of carbohydrates and TCA intermediate metabolites as well as corresponding genes were linearly increased with mushroom size. However, the composition of flavor-related metabolites was varying in different sized beech mushrooms. Our study explores the signature pathways associated with growth, development, nutritional, functional and organoleptic properties of different sized beech mushrooms.

Element distribution in fruiting bodies of Lactarius pubescens with focus on rare earth elements

During growth and senescence, fungal fruiting bodies accumulate essential and non-essential elements to different extent in their compartments. This study bases on a dataset of 32 basidiocarps of the ectomycorrhizal Lactarius pubescens sampled in a former U mining area. Statistical analyses were combined with rare earth element (REE, La-Lu) patterns to study the element distribution within sporocarp compartments and between three different age classes. For this purpose, fruiting bodies were separated into stipe, pileus trama, pileipelles and lamellae, dried and digested with HNO₃. While macronutrient (e.g. K, Mg, P, S) contents resemble those of a non-mining affected site, several elements (e.g. Co, Mn) were site-specifically taken up relative to elevated soil contents. With statistics, two main element distribution groups for L. pubescens were revealed: mainly essential (Cu, Mg, Mn, P, S, Zn, Cd, Co, Ni) and mainly non-essential elements (Al, Ca, Fe, Sr, U, REE). The highest REE contents were found in pileipelles and lamellae, corresponding to relatively small cell sizes. Stipes and pileus trama had low REE contents due to their function as transport systems. During growth, light REE (La-Nd) were strongly enriched in lamellae and pileipelles. Middle REE (Sm-Dy) enrichment was found both in soil and fungal biomass. Contents of nutrients decrease with age, while non-essential elements are enriched especially in pileipelles and lamellae. A weak positive Ce anomaly appeared in the bioavailable soil fraction and in the pileipelles of younger individuals. Substrate dependent uptake thus gets reduced with sporocarp senescence, possibly due to redistribution.

Genome-wide gene expression patterns in dikaryon of the basidiomycete fungus Pleurotus ostreatus

Dikarya is a subkingdom of fungi that includes Ascomycota and Basidiomycota. The gene expression patterns of dikaryon are poorly understood. In this study, we bred a dikaryon DK13×3 by mating monokaryons MK13 and MK3, which were from the basidiospores of Pleurotus ostreatus TD300. Using RNA-Seq, we obtained the transcriptomes of the three strains. We found that the total transcript numbers in the transcriptomes of the three strains were all more than ten thousand, and the expression profile in DK13×3 was more similar to MK13 than MK3. However, the genes involved in macromolecule utilization, cellular material synthesis, stress-resistance and signal transduction were much more up-regulated in the dikaryon than its constituent monokaryons. All possible modes of differential gene expression, when compared to constituent monokaryons, including the presence/absence variation, and additivity/nonadditivity gene expression in the dikaryon may contribute to heterosis. By sequencing the urease gene poure sequences and mRNA sequences, we identified the monoallelic expression of the poure gene in the dikaryon, and its transcript was from the parental monokaryon MK13. Furthermore, we discovered RNA editing in the poure gene mRNA of the three strains. These results suggest that the gene expression patterns in dikaryons should be similar to that of diploids during vegetative growth.

Environmental regulation of reproductive phase change in Agaricus bisporus by 1-octen-3-ol, temperature and CO₂

Reproductive phase change from vegetative mycelium to the initiation of fruiting in Agaricus bisporus is regulated in large part by the sensing of environmental conditions. A model is proposed in which three separate environmental factors exert control at different stages of the reproductive developmental process change. The eight carbon volatile 1-octen-3-ol controls the early differentiation from vegetative hyphae to multicellular knots; temperature reduction is essential for the later differentiation of primodia; and carbon dioxide level exerts quantitative control on the number of fruiting bodies developed. Analysis of transcriptomic changes during the reproductive phase change was carried out with initiation-specific microarrays, and the newly published A. bisporus genome was used to analyse the promoter regions of differentially regulated genes. Our studies have shown there to be both early and late initiation responses relating to sensing of eight carbon volatiles and temperature respectively. A subset of 45 genes was transcriptionally regulated during the reproductive phase change which exhibited a range of functions including cell structure, nitrogen and carbon metabolism, and sensing and signalling. Three gene clusters linking increased transcription with developmental stage were identified. Analysis of promoter regions revealed cluster-specific conserved motifs indicative of co-ordinated regulation of transcription.

Hairpin-mediated down-regulation of the urea cycle enzyme argininosuccinate lyase in Agaricus bisporus

A double-stranded (ds) RNA hairpin-mediated down-regulation system was developed for the cultivated mushroom Agaricus bisporus, and the role of the urea cycle enzyme argininosuccinate lyase (asl) in mushroom post-harvest development was investigated. Hairpin expression vectors were constructed to initiate down-regulation of asl and introduced into A. bisporus by Agrobacterium tumefaciens-mediated transformation. Transcripts of asl were significantly reduced (93.1 and 99.9%) in two transformants and hairpin vector transgene sequences were maintained throughout sporophore development. Single and multiple hairpin integration events were observed in Southern analysis. Transformants with down-regulated asl exhibited reduced yield and cap expansion during post-harvest sporophore development. There were no detectable differences in urea levels between the hairpin-transformed and control strains. This is the first report of reduced gene expression resulting from the introduction of dsRNA hairpins in A. bisporus and the applications of this technology will facilitate functional studies in the mushroom.

A gene repertoire for nitrogen transporters in Laccaria bicolor

Ectomycorrhizal interactions established between the root systems of terrestrial plants and hyphae from soil-borne fungi are the most ecologically widespread plant symbioses. The efficient uptake of a broad range of nitrogen (N) compounds by the fungal symbiont and their further transfer to the host plant is a major feature of this symbiosis. Nevertheless, we far from understand which N form is preferentially transferred and what are the key molecular determinants required for this transfer. Exhaustive in silico analysis of N-compound transporter families were performed within the genome of the ectomycorrhizal model fungus Laccaria bicolor. A broad phylogenetic approach was undertaken for all families and gene regulation was investigated using whole-genome expression arrays. A repertoire of proteins involved in the transport of N compounds in L. bicolor was established that revealed the presence of at least 128 gene models in the genome of L. bicolor. Phylogenetic comparisons with other basidiomycete genomes highlighted the remarkable expansion of some families. Whole-genome expression arrays indicate that 92% of these gene models showed detectable transcript levels. This work represents a major advance in the establishment of a transportome blueprint at a symbiotic interface, which will guide future experiments.

GC-MS studies of the chemical composition of two inedible mushrooms of the genus Agaricus

Background

Mushrooms in the genus Agaricus have worldwide distribution and include the economically important species A. bisporus. Some Agaricus species are inedible, including A. placomyces and A. pseudopratensis, which are similar in appearance to certain edible species, yet are known to possess unpleasant odours and induce gastrointestinal problems if consumed. We have studied the chemical composition of these mushrooms using GC-MS.

Results

Our GC-MS studies on the volatile fractions and butanol extracts resulted in the identification of 44 and 34 compounds for A. placomyces and A. pseudopratensis, respectively, including fatty acids and their esters, amino acids, and sugar alcohols. The most abundant constituent in the volatiles and butanol were phenol and urea respectively. We also identified the presence of ergosterol and two Δ⁷-sterols. In addition, 5α,8α-Epidioxi-24(ξ)-methylcholesta-6,22-diene-3β-ol was isolated for the first time from both mushrooms. Our study is therefore the first report on the chemical composition of these two species.

Conclusion

The results obtained contribute to the knowledge of the chemical composition of mushrooms belonging to the Agaricus genus, and provide some explanation for the reported mild toxicity of A. placomyces and A. pseudopratensis, a phenonomenon that can be explained by a high phenol content, similar to that found in other Xanthodermatei species.

The role of ornithine aminotransferase in fruiting body formation of the mushroom Agaricus bisporus

The complete oat gene and cDNA from the commercial mushroom, Agaricus bisporus, encoding ornithine aminotransferase (OAT) was characterized. The gene encodes a 466 amino acid protein and provides the first fully reported homobasidiomycete OAT protein sequence. The gene is interrupted by ten introns, and no mitochondrial targeting motif was present pointing to a cytoplasmic localization. The function of the gene was demonstrated by complementation of a Saccharomyces cerevisiae mutant unable to utilize ornithine as a sole source of nitrogen with an A. bisporus oat cDNA construct. Northern analysis of the oat gene together with the pruA gene (encoding Delta(1)-pyrroline-5-carboxylate dehydrogenase) showed that transcripts of both genes were lower during the first stages of fruiting body development. The higher expression of the oat gene in later stages of development, suggests the importance of ornithine metabolism for the redistribution of metabolites in the developing mushroom. Hplc analysis of all amino acids revealed that ornithine levels increased during fruiting body development whereas proline levels fell.

Argininosuccinate synthetase and argininosuccinate lyase: two ornithine cycle enzymes from Agaricus bisporus

Accumulation of high quantities of urea in fruiting bodies is a known feature of larger basidiomycetes. Argininosuccinate synthetase (ASS) and argininosuccinate lyase (ASL) are two ornithine cycle enzymes catalysing the last two steps in the arginine biosynthetic pathway. Arginine is the main precursor for urea formation. In this work the nucleotide sequences of the genes and corresponding cDNAs encoding argininosuccinate synthetase (ass) and argininosuccinate lyase (asl) from Agaricus bisporus were determined. Eight and six introns were present in the ass and asl gene, respectively. The location of four introns in the asl gene were conserved among vertebrate asl genes. Deduced amino acid sequences, representing the first homobasidiomycete ASS and ASL protein sequences, were analysed and compared with their counterparts in other organisms. The ass ORF encoded for a protein of 425 amino acids with a calculated molecular mass of 47266Da. An alignment with ASS proteins from other organisms revealed high similarity with fungal and mammalian ASS proteins, 61-63% and 51-55% identity, respectively. The asl open reading frame (ORF) encoded a protein of 464 amino acids with an calculated mass of 52337Da and similar to ASS shared the highest similarity with fungal ASL proteins, 59-60% identity. Northern analyses of ass and asl during fruiting body formation and post-harvest development revealed that expression was significantly up-regulated from developmental stage 3 on for all the tissues studied. The expression reached a maximum at the later stages of fruiting body growth, stages 6 and 7. Both ass and asl genes were up-regulated within 3h after harvest showing that the induction mechanism is very sensitive to the harvest event and emphasizes the importance of the arginine biosynthetic pathway/ornithine cycle in post-harvest physiology.