Genome Sequence of the Edible Cultivated Mushroom Lentinula edodes (Shiitake) Reveals Insights into Lignocellulose Degradation

Complete mitochondrial genomes of edible mushrooms: features, evolution, and phylogeny

Edible mushrooms are an important food source with high nutritional and medicinal value. They are a useful source for studying phylogenetic evolution and species divergence. The exploration of the evolutionary relationships among these species conventionally involves analyzing sequence variations within their complete mitochondrial genomes, which range from 31,854 bp (Cordyceps militaris) to 197,486 bp (Grifolia frondosa). The study of the complete mitochondrial genomes of edible mushrooms has emerged as a critical field of research, providing important insights into fungal genetic makeup, evolution, and phylogenetic relationships. This review explores the mitochondrial genome structures of various edible mushroom species, highlighting their unique features and evolutionary adaptations. By analyzing these genomes, robust phylogenetic frameworks are constructed to elucidate mushrooms lineage relationships. Furthermore, the exploration of different variations of mitochondrial DNA presents novel opportunities for enhancing mushroom cultivation biotechnology and medicinal applications. The mitochondrial genomic features are essential for improving agricultural practices and ensuring food security through improved crop productivity, disease resistance, and nutritional qualities. The current knowledge about the mitochondrial genomes of edible mushrooms is summarized in this review, emphasising their significance in both scientific research and practical applications in bioinformatics and medicine.

Near-gapless genome and transcriptome analyses provide insights into fruiting body development in Lentinula edodes

Fruiting body development in macrofungi is an intensive research subject. In this study, high-quality genomes were assembled for two sexually compatible monokaryons from a heterokaryotic Lentinula edodes strain WX1, and variations in L. edodes genomes were analyzed. Specifically, differential gene expression and allele-specific expression (ASE) were analyzed using the two monokaryotic genomes and transcriptome data from four different stages of fruiting body development in WX1. Results revealed that after aeration, mycelia sensed cell wall stress, pheromones, and a decrease in CO2 concentration, leading to up-regulated expression in genes related to cell adhesion, cell wall remodeling, proteolysis, and lipid metabolism, which may promote primordium differentiation. Aquaporin genes and those related to proteolysis, mitosis, lipid, and carbohydrate metabolism may play important roles in primordium development, while genes related to tissue differentiation and sexual reproduction were active in fruiting body. Several essential genes for fruiting body development were allele-specifically expressed and the two nuclear types could synergistically regulate fruiting body development by dominantly expressing genes with different functions. ASE was probably induced by long terminal repeat-retrotransposons. Findings here contribute to the further understanding of the mechanism of fruiting body development in macrofungi.

Pleurotus ostreatus as a model mushroom in genetics, cell biology, and material sciences

Pleurotus ostreatus, also known as the oyster mushroom, is a popular edible mushroom cultivated worldwide. This review aims to survey recent progress in the molecular genetics of this fungus and demonstrate its potential as a model mushroom for future research. The development of modern molecular genetic techniques and genome sequencing technologies has resulted in breakthroughs in mushroom science. With efficient transformation protocols and multiple selection markers, a powerful toolbox, including techniques such as gene knockout and genome editing, has been developed, and numerous new findings are accumulating in P. ostreatus. These include molecular mechanisms of wood component degradation, sexual development, protein secretion systems, and cell wall structure. Furthermore, these techniques enable the identification of new horizons in enzymology, biochemistry, cell biology, and material science through protein engineering, fluorescence microscopy, and molecular breeding. KEY POINTS: • Various genetic techniques are available in Pleurotus ostreatus. • P. ostreatus can be used as an alternative model mushroom in genetic analyses. • New frontiers in mushroom science are being developed using the fungus.

Comparative Transcriptome Analysis of Candidate Genes Associated with Mycelia Growth from a He-Ne Laser with Pulsed Light Mutant of Phellinus igniarius (Agaricomycetes)

A mutant Phellinus igniarius JQ9 with higher mycelial production was screened out by He-Ne laser with pulsed light irradiation, the mechanism underlying the higher mycelial production is still unknown. This study aims to obtain a comprehensive transcriptome assembly during the Ph. igniarius liquid fermentation and characterize the key genes associated with the mycelial growth and metabolism in Ph. igniarius JQ9. Our transcriptome data of Ph. iniarius JQ9 and the wild strain were obtained with the Illumina platform comparative transcriptome sequencing technology. The results showed that among all the 346 differentially expressed genes (DEGs), 245 were upregulated and 101 were downregulated. Candidate genes encoding endoglucanase, beta-glucosidase, cellulose 1,4-beta-cellobiosidase, glycoside hydrolase family 61 protein, were proposed to participate in the carbohydrate utilization from KEGG enrichment of the starch and sucrose metabolism pathways were upregulated in Ph. igniarius JQ9. In addition, three candidate genes encoding the laccase and another two candidate genes related with the cell growth were higher expressed in Ph. igniarius JQ9 than in the wild type of strain (CK). Analysis of these data revealed that increased these related carbohydrate metabolism candidate genes underlying one crucial way may cause the higher mycelia production.

Hydrated lime promoted the polysaccharide content and affected the transcriptomes of Lentinula edodes during brown film formation

Brown film formation, a unique developmental stage in the life cycle of Lentinula edodes, is essential for the subsequent development of fruiting bodies in L. edodes cultivation. The pH of mushroom growth substrates are usually adjusted with hydrated lime, yet the effects of hydrated lime on cultivating L. edodes and the molecular mechanisms associated with the effects have not been studied systemically. We cultivated L. edodes on substrates supplemented with 0% (CK), 1% (T1), 3% (T2), and 5% (T3) hydrated lime (Ca (OH)2), and applied transcriptomics and qRT-PCR to study gene expression on the brown film formation stage. Hydrated lime increased polysaccharide contents in L. edodes, especially in T2, where the 5.3% polysaccharide content was approximately 1.5 times higher than in the CK. The addition of hydrated lime in the substrate promoted laccase, lignin peroxidase and manganese peroxidase activities, implying that hydrated lime improved the ability of L. edodes to decompose lignin and provide nutrition for its growth and development. Among the annotated 9,913 genes, compared to the control, 47 genes were up-regulated and 52 genes down-regulated in T1; 73 genes were up-regulated and 44 were down-regulated in T2; and 125 genes were up-regulated and 65 genes were down-regulated in T3. Differentially expressed genes (DEGs) were enriched in the amino acid metabolism, lipid metabolism and carbohydrate metabolism related pathways. The carbohydrate-active enzyme genes up-regulated in the hydrated lime treatments were mostly glycosyl hydrolase genes. The results will facilitate future optimization of L. edodes cultivation techniques and possibly shortening the production cycle.

Decoding the chromosome-scale genome of the nutrient-rich Agaricus subrufescens: a resource for fungal biology and biotechnology

Agaricus subrufescens, also known as the "sun mushroom," has significant nutritional and medicinal value. However, its short shelf life due to the browning process results in post-harvest losses unless it's quickly dehydrated. This restricts its availability to consumers in the form of capsules. A genome sequence of A. subrufescens may lead to new cultivation alternatives or the application of gene editing strategies to delay the browning process. We assembled a chromosome-scale genome using a hybrid approach combining Illumina and Nanopore sequencing. The genome was assembled into 13 chromosomes and 31 unplaced scaffolds, totaling 44.5 Mb with 96.5% completeness and 47.24% GC content. 14,332 protein-coding genes were identified, with 64.6% of the genome covered by genes and 23.41% transposable elements. The mitogenome was circularized and encoded fourteen typical mitochondrial genes. Four polyphenol oxidase (PPO) genes and the Mating-type locus were identified. Phylogenomic analysis supports the placement of A. subrufescens in the Agaricomycetes clade. This is the first available genome sequence of a strain of the "sun mushroom." Results are available through a Genome Browser (https://plantgenomics.ncc.unesp.br/gen.php?id=Asub) and can support further fungal biological and genomic studies.

Multi-stage nuclear transcriptomic insights of morphogenesis and biparental role changes in Lentinula edodes

Based on six offspring with different mitochondrial (M) and parental nuclear (N) genotypes, the multi-stage morphological characteristics and nuclear transcriptomes of Lentinula edodes were compared to investigate morphogenesis mechanisms during cultivation, the key reason for cultivar resistance to genotype changes, and regulation related to biparental role changes. Six offspring had specific transcriptomic data and morphological characteristics that were mainly regulated by the two parental nuclei, followed by the cytoplasm, at different growth stages. Importing a wild N genotype easily leads to failure or instability of fruiting; however, importing wild M genotypes may improve cultivars. Major facilitator superfamily (MFS) transporter genes encoding specific metabolites in spawns may play crucial roles in fruiting body formation. Pellets from submerged cultivation and spawns from sawdust substrate cultivation showed different carbon metabolic pathways, especially in secondary metabolism, degradation of lignin, cellulose and hemicellulose, and plasma membrane transport (mainly MFS). When the stage of small young pileus (SYP) was formed on the surface of the bag, the spawns inside were mainly involved in nutrient accumulation. Just broken pileus (JBP) showed a different expression of plasma membrane transporter genes related to intracellular material transport compared to SYP and showed different ribosomal proteins and cytochrome P450 functioning in protein biosynthesis and metabolism than near spreading pileus (NSP). Biparental roles mainly regulate offspring metabolism, growth, and morphogenesis by differentially expressing specific genes during different vegetative growth stages. Additionally, some genes encoding glycine-rich RNA-binding proteins, F-box, and folliculin-interacting protein repeat-containing proteins may be related to multi-stage morphogenesis. KEY POINTS: • Replacement of nuclear genotype is not suitable for cultivar breeding of L. edodes. • Some genes show a biparental role-divergent expression at mycelial growth stage. • Transcriptomic changes of some sawdust substrate cultivation stages have been elucidated.

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.

Lignin-degrading enzymes from a pathogenic canker-rot fungus Inonotus obliquus strain IO-B2

Inonotus obliquus is a pathogenic fungus found in living trees and has been widely used as a traditional medicine for cancer therapy. Although lignocellulose-degrading enzymes are involved in the early stages of host infection, the parasitic life cycle of this fungus has not been fully understood. In this study, we aimed to investigate the activities of laccase (Lac), manganese peroxidase (MnP), and lignin peroxidase (LiP) from I. obliquus cultivated in Kirk's medium. The fungus was subjected to genome sequencing, and genes related to wood degradation were identified. The draft genome sequence of this fungus comprised 21,203 predicted protein-coding genes, of which 134 were estimated to be related to wood degradation. Among these, 47 genes associated with lignin degradation were found to have the highest number of mnp genes. Furthermore, we cloned the cDNA encoding a putative MnP, referred to as IoMnP1, and characterized its molecular structure. The results show that IoMnP1 has catalytic properties analogous to MnP. Phylogenetic analysis also confirmed that IoMnP1 was closely related to the MnPs from Pyrrhoderma noxium, Fomitiporia mediterranea, and Sanghuangporus baumii, which belong to the same family of Hymenochaetaceae. From the above results, we suggest that IoMnP1 is a member of MnPs.

Functional Clustering of Metabolically Related Genes Is Conserved across Dikarya

Transcriptional regulation is vital for organismal survival, with many layers and mechanisms collaborating to balance gene expression. One layer of this regulation is genome organization, specifically the clustering of functionally related, co-expressed genes along the chromosomes. Spatial organization allows for position effects to stabilize RNA expression and balance transcription, which can be advantageous for a number of reasons, including reductions in stochastic influences between the gene products. The organization of co-regulated gene families into functional clusters occurs extensively in Ascomycota fungi. However, this is less characterized within the related Basidiomycota fungi despite the many uses and applications for the species within this clade. This review will provide insight into the prevalence, purpose, and significance of the clustering of functionally related genes across Dikarya, including foundational studies from Ascomycetes and the current state of our understanding throughout representative Basidiomycete species.

A global phylogenomic analysis of the shiitake genus Lentinula

Lentinula is a broadly distributed group of fungi that contains the cultivated shiitake mushroom, L. edodes. We sequenced 24 genomes representing eight described species and several unnamed lineages of Lentinula from 15 countries on four continents. Lentinula comprises four major clades that arose in the Oligocene, three in the Americas and one in Asia-Australasia. To expand sampling of shiitake mushrooms, we assembled 60 genomes of L. edodes from China that were previously published as raw Illumina reads and added them to our dataset. Lentinula edodes sensu lato (s. lat.) contains three lineages that may warrant recognition as species, one including a single isolate from Nepal that is the sister group to the rest of L. edodes s. lat., a second with 20 cultivars and 12 wild isolates from China, Japan, Korea, and the Russian Far East, and a third with 28 wild isolates from China, Thailand, and Vietnam. Two additional lineages in China have arisen by hybridization among the second and third groups. Genes encoding cysteine sulfoxide lyase (lecsl) and γ-glutamyl transpeptidase (leggt), which are implicated in biosynthesis of the organosulfur flavor compound lenthionine, have diversified in Lentinula. Paralogs of both genes that are unique to Lentinula (lecsl 3 and leggt 5b) are coordinately up-regulated in fruiting bodies of L. edodes. The pangenome of L. edodes s. lat. contains 20,308 groups of orthologous genes, but only 6,438 orthogroups (32%) are shared among all strains, whereas 3,444 orthogroups (17%) are found only in wild populations, which should be targeted for conservation.

Laccase mediated delignification of wasted and non-food agricultural biomass: Recent developments and challenges

Utilization of microbial laccases is considered as the cleaner and target specific biocatalytic mechanism for the recovery of cellulose and hemicelluloses from nonfood and wasted agricultural, lignocellulosic biomass (LCB). The extent of lignin removal by laccase depends on the biochemical composition of biomass and the redox potential (E⁰) of the biocatalyst. Intensive research efforts are going on all over the world for the recognition of appropriate and easily available agricultural lignocellulosic feedstocks to exploit maximally for the production of value-added bioproducts and biofuels. In such circumstances, laccase can play a major role as a leading biocatalyst and potent substitute for chemical based deconstruction of the lignocellulosic materials. The limited commercialization of laccase at an industrial scale has been feasible due to its full working efficiency mostly expressed in the presence of cost intensive redox mediators only. Although, recently there are some reports that came on the mediator free biocatalysis of enzyme but still not considerably explored and neither understood in depth. The present review will address the various research gaps and shortcomings that acted as the big hurdles before the complete exploitation of laccases at an industrial scale. Further, this article also reveals insights on different microbial laccases and their diverse functional environmental conditions that affect the deconstruction process of LCB.

Transcriptome Analysis of Antrodia cinnamomea Mycelia from Different Wood Substrates

Antrodia cinnamomea, an edible and medicinal fungus with significant economic value and application prospects, is rich in terpenoids, benzenoids, lignans, polysaccharides, and benzoquinone, succinic and maleic derivatives. In this study, the transcriptome of A. cinnamomea cultured on the wood substrates of Cinnamomum glanduliferum (YZM), C. camphora (XZM), and C. kanehirae (NZM) was sequenced using the high-throughput sequencing technology Illumina HiSeq 2000, and the data were assembled by de novo strategy to obtain 78,729 Unigenes with an N50 of 4,463 bp. Compared with public databases, about 11,435, 6,947, and 5,994 Unigenes were annotated to the Non-Redundant (NR), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genome (KEGG), respectively. The comprehensive analysis of the mycelium terpene biosynthesis-related genes in A. cinnamomea revealed that the expression of acetyl-CoA acetyltransferase (AACT), acyl-CoA dehydrogenase (MCAD), 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA), mevalonate pyrophosphate decarboxylase (MVD), and isopentenyl diphosphate isomerase (IDI) was significantly higher on NZM compared to the other two wood substrates. Similarly, the expression of geranylgeranyltransferase (GGT) was significantly higher on YZM compared to NZM and XZM, and the expression of farnesyl transferase (FTase) was significantly higher on XZM. Furthermore, the expressions of 2,3-oxidized squalene cyclase (OCS), squalene synthase (SQS), and squalene epoxidase (SE) were significantly higher on NZM. Overall, this study provides a potential approach to explore the molecular regulation mechanism of terpenoid biosynthesis in A. cinnamomea.

Lignin Promotes Mycelial Growth and Accumulation of Polyphenols and Ergosterol in Lentinula edodes

It has been demonstrated that lignin was efficiently degraded by Lentinula edodes (L. edodes). However, the process of lignin degradation and utilization by L. edodes has not been discussed in detail. Therefore, the effects of lignin on L. edodes mycelium growth, chemical compositions, and phenolic profiles were investigated herein. It has been revealed that 0.10% lignin acted as the most effective concentration to accelerate mycelia growth, which yielded the highest biomass of 5.32 ± 0.07 g/L. Furthermore, a 0.10% concentration of lignin promoted the accumulation of phenolic compounds, especially protocatechuic acid, with peak value of 48.5 ± 1.2 μg/g. In contrast, the higher concentration of lignin (0.20%) exerted an inhibitory effect on the growth of L. edodes. Overall, the application of lignin at the optimal concentration of 0.10% could not only enhance the mycelial growth but also accumulate the phenolic acids and raise the nutritional and medical values of L. edodes.

Increasing the production of the bioactive compounds in medicinal mushrooms: an omics perspective

Macroscopic fungi, mainly higher basidiomycetes and some ascomycetes, are considered medicinal mushrooms and have long been used in different areas due to their pharmaceutically/nutritionally valuable bioactive compounds. However, the low production of these bioactive metabolites considerably limits the utilization of medicinal mushrooms both in commerce and clinical trials. As a result, many attempts, ranging from conventional methods to novel approaches, have been made to improve their production. The novel strategies include conducting omics investigations, constructing genome-scale metabolic models, and metabolic engineering. So far, genomics and the combined use of different omics studies are the most utilized omics analyses in medicinal mushroom research (both with 31% contribution), while metabolomics (with 4% contribution) is the least. This article is the first attempt for reviewing omics investigations in medicinal mushrooms with the ultimate aim of bioactive compound overproduction. In this regard, the role of these studies and systems biology in elucidating biosynthetic pathways of bioactive compounds and their contribution to metabolic engineering will be highlighted. Also, limitations of omics investigations and strategies for overcoming them will be provided in order to facilitate the overproduction of valuable bioactive metabolites in these valuable organisms.

The genome of Lyophyllum shimeji provides insight into the initial evolution of ectomycorrhizal fungal genomes

Mycorrhizae are one of the most fundamental symbioses between plants and fungi, with ectomycorrhizae being the most widespread in boreal forest ecosystems. Ectomycorrhizal fungi are hypothesized to have evolved convergently from saprotrophic ancestors in several fungal clades, especially members of the subdivision Agaricomycotina. Studies on fungal genomes have identified several typical characteristics of mycorrhizal fungi, such as genome size expansion and decreases in plant cell-wall degrading enzymes (PCWDEs). However, genomic changes concerning the evolutionary transition to the ectomycorrhizal lifestyle are largely unknown. In this study, we sequenced the genome of Lyophyllum shimeji, an ectomycorrhizal fungus that is phylogenetically related to saprotrophic species and retains some saprotroph-like traits. We found that the genome of Ly. shimeji strain AT787 lacks both incremental increases in genome size and reduced numbers of PCWDEs. Our findings suggest that the previously reported common genomic traits of mycorrhizal fungi are not essential for the ectomycorrhizal lifestyle, but are a result of abolishing saprotrophic activity. Since Ly. shimeji is commercially consumed as an edible mushroom, the newly available genomic information may also impact research designed to enhance the cultivation of this mushroom.

Complete genome sequences and comparative secretomic analysis for the industrially cultivated edible mushroom Lyophyllum decastes reveals insights on evolution and lignocellulose degradation potential

Lyophyllum decastes, also known as Luronggu in China, is a culinary edible and medicinal mushroom that was widely cultivated in China in recent years. In the present study, the complete high-quality genome of two mating compatible L. decastes strain was sequenced. The L. decastes LRG-d1-1 genome consists of 47.7 Mb in 15 contigs with a contig N90 of 2.08 Mb and 14,499 predicted gene models. Phylogenetic analysis revealed that L. decastes exhibits a close evolutionary relationship to the Termitomyces and Hypsizygus genus and was diverged from H. marmoreus  ~ 45.53 Mya ago. Mating A loci of L. decastes compose of five and four HD genes in two monokaryotic strains, respectively. Mating B loci compose of five STE genes in both two monokaryotic strains. To accelerate the cross-breeding process, we designed four pairs of specific primers and successfully detected both mating types in L. decastes. As a wood-rotting mushroom, a total of 541 genes accounting for 577 CAZymes were identified in the genome of L. decastes. Proteomic analysis revealed that 1,071 proteins including 182 CAZymes and 258 secreted enzymes were identified from four groups (PDB, PDB + bran, PDB + cotton hull, and PDB + sawdust). Two laccases and a quinone reductase were strongly overproduced in lignin-rich cultures, and the laccases were among the top-3 secreted proteins, suggesting an important role in the synergistic decomposition of lignin. These results revealed the robustness of the lignocellulose degradation capacity of L. decastes. This is the first study to provide insights into the evolution and lignocellulose degradation of L. decastes.

Proteomic characterization of Shiitake (Lentinula edodes) post-harvest fruit bodies grown on hardwood logs and isolation of an antibacterial serine protease inhibitor

Lentinula edodes (Shiitake) is one of the most heavily cultivated mushrooms in the world with proven antioxidant and antibacterial properties, among others. Evidence indicates that the choice of mushroom cultivation technique strongly influences the production of bioactive compounds, but to date the nature of many of these compounds has not been fully established. This work focuses on the proteomic characterization of L. edodes to highlight the main active processes two days after harvest and elucidates the proteins involved in the known antioxidant and antibacterial proprieties of Shiitake fruit bodies cultivated on oak logs. A label-free approach allowed us to identify a total of 2702 proteins which were mainly involved in carbohydrate and protein metabolism, cell growth and replication, indicating that several developmental processes remain active in fruit bodies post-harvest. Proteins with antioxidant activities were identified, indicating the contribution of proteins to the antioxidant properties of L. edodes extracts. Antibacterial assays also reveal the activity of a serine protease inhibitor that strongly accumulates in the post-harvest fruit body grown on oak logs. Overall, this study contributes to the understanding of the impact of the log cultivation method on the production of Shiitake mushrooms richest in high-value bioactive compounds.

Characteristics of the Genome, Transcriptome and Ganoderic Acid of the Medicinal Fungus Ganoderma lingzhi

Ganoderma (Ganodermaceae) is a genus of edible and medicinal mushrooms that create a diverse set of bioactive compounds. Ganoderma lingzhi has been famous in China for more than 2000 years for its medicinal properties. However, the genome information of G. lingzhi has not been characterized. Here, we characterized its 49.15-Mb genome, encoding 13,125 predicted genes which were sequenced by the Illumina and PacBio platform. A wide spectrum of carbohydrate-active enzymes, with a total number of 519 CAZymes were identified in G. lingzhi. Then, the genes involved in sexual recognition and ganoderic acid (GA, key bioactive metabolite) biosynthesis were characterized. In addition, we identified and deduced the possible structures of 20 main GA constituents by UPLC-ESI-MS/MS, including a new special ganochlearic acid A. Furthermore, 3996 novel transcripts were discovered, and 9276 genes were predicted to have the possibility of alternative splicing from RNA-Seq data. The alternative splicing genes were enriched for functional categories involved in protein processing, endocytosis, and metabolic activities by KEGG. These genomic, transcriptomic, and GA constituents’ resources would enrich the toolbox for biological, genetic, and secondary metabolic pathways studies in G. lingzhi.

Comparative structural analysis on the mitochondrial DNAs from various strains of Lentinula edodes

The evolution of mitochondria through variations in mitochondrial DNA (mtDNA) is one of the intriguing questions in eukaryotic cells. In order to assess the causes of the variations in mitochondria, the mtDNAs of the 21 strains of Lentinula edodes were assembled for this study, and analyzed together with four published mtDNA sequences. The mtDNAs were within the sizes of 117 kb ~ 122 kb. The gene number was observed consistent except for two mtDNAs, which carry a duplicated trnG1-trnG2 unit or a putative gene deletion. The size variation was largely attributed to the number of introns, repeated sequences, transposable elements (TEs), and plasmid-related sequences. Intron loss and gain were found from cox1, rnl, and rns of three mtDNAs. Loss of two introns in cox1 of KY217797.1 reduced its size by 2.7 kb, making it the smallest cox1 gene (8.4 kb) among the cox1s of the 25 mtDNAs, whereas gain of a Group II intron (2.65 kb) and loss of a Group I intron (1.7 kb) in cox1 of MF774813.1 resulted in the longest cox1 (12 kb). In rnl of L. edodes, we discovered four intron insertion consensus sequences which were unique to basidiomycetes but not ascomycetes. Differential incorporation of introns was the primary cause of the rnl size polymorphism. Homing endonucleases (HEGs) were suggestively involved in the mobilization of the introns because all of the introns have HEG genes of the LAGRIDADG or GIY-YIG families with the conserved HEG cleavage sites. TEs contributed to 11.04% of the mtDNA size in average, of which 7.08% was LTR-retrotransposon and 3.96% was DNA transposon, whereas the repeated sequences covered 4.6% of the mtDNA. The repeat numbers were variable in a strain-dependent manner. Both the TEs and repeated sequences were mostly found in the intronic and intergenic regions. Lastly, two major deletions were found in the plasmid-related sequence regions (pol2-pol3 and pol1-atp8) in the five mtDNAs. Particularly, the 6.8 kb-long deletion at pol2-pol3 region made MF774813.1 the shortest mtDNA of all. Our results demonstrate that mtDNA is a dynamic molecule that persistently evolves over a short period of time by insertion/deletion and repetition of DNA segments at the strain level.

Comparative Transcriptome Analysis Provides Insights Into the Mechanism by Which 2,4-Dichlorophenoxyacetic Acid Improves Thermotolerance in Lentinula edodes

As the widest cultivated edible mushroom worldwide, Lentinula edodes suffers serious yield and quality losses from heat stress during growth and development, and in our previous study, exogenous 2,4-Dichlorophenoxyacetic acid (2,4-D) was found to improve the thermotolerance of L. edodes strain YS3357, but the molecular mechanism remains unclear. Here, we explored the potential protective mechanism of exogenous 2,4-D against heat stress by transcriptome analysis. 2,4-D possible improve the thermotolerance of L. edodes through regulating antioxidant genes, transcription factors, energy-provision system, membrane fluidity, and cell wall remodeling. Furthermore, 2,4-D was also found to regulate the saturation levels of fatty acids and ATP content in L. edodes mycelium under heat stress. This study proposed a regulatory network of 2,4-D in regulating L. edodes response to heat stress, providing a theoretical basis for improving L. edodes thermotolerance, and facilitating the understanding of the molecular mechanism of exogenous hormones in alleviating abiotic stress damage to macrofungi.

Transcriptome and Differentially Expressed Gene Profiles in Mycelium, Primordium and Fruiting Body Development in Stropharia rugosoannulata

Stropharia rugosoannulata uses straw as a growth substrate during artificial cultivation and has been widely promoted in China. However, its fruiting body formation and development processes have not been elucidated. In this study, the developmental transcriptomes were analyzed at three stages: the mycelium (G-S), primordium (P-S) and fruiting body (M-F) stages. A total of 9690 differentially expressed genes (DEGs) were identified in the different developmental stages. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that these DEGs were involved mainly in hydrolase activity, structural molecule activity and oxidoreductase activity as well as xenobiotic biodegradation and metabolism and energy metabolism pathways. We further found that the higher expression of most carbohydrate enzyme (i.e., GH, CE, CBM, AA and PL) genes in the hyphal (i.e., G-S) stage was related mainly to substrate degradation, while the upregulation of glycosyltransferase (GT) gene expression in the P-S and M-F stages may be related to cell wall synthesis. In addition, we found that CO₂-sensing-related genes (i.e., CA-2, CA-3, PKA-1 and PKA-2) were upregulated in the P-S and M-F stages, heat shock protein genes (HSP60 and HSP90) were significantly downregulated in the P-S stage and upregulated in the M-F stage and the transcription factors (i.e., steA, MYB, nosA, HAP1, and GATA-4/5/6) involved in growth and development were significantly upregulated in the P-S stage. These results suggest that environmental factors (i.e., CO₂ and temperature) and transcription factors may play a key role in primordium formation. In short, this study provides new insights into the study of stimulating primordia formation affecting the development of fruiting bodies of S. rugosoannulata.

Systems biology-guided understanding of white-rot fungi for biotechnological applications: a review

Plant-derived biomass is the most abundant biogenic carbon source on Earth. Despite this, only a small clade of organisms known as white-rot fungi (WRF) can efficiently break down both the polysaccharide and lignin components of plant cell walls. This unique ability imparts a key role for WRF in global carbon cycling and highlights their potential utilization in diverse biotechnological applications. To date, research on WRF has primarily focused on their extracellular ‘digestive enzymes’ whereas knowledge of their intracellular metabolism remains underexplored. Systems biology is a powerful approach to elucidate biological processes in numerous organisms, including WRF. Thus, here we review systems biology methods applied to WRF to date, highlight observations related to their intracellular metabolism, and conduct comparative extracellular proteomic analyses to establish further correlations between WRF species, enzymes, and cultivation conditions. Lastly, we discuss biotechnological opportunities of WRF as well as challenges and future research directions.

Global phylogeny of the Shiitake mushroom and related Lentinula species uncovers novel diversity and suggests an origin in the Neotropics

Lentinula (Basidiomycota, Agaricales) includes the most widely cultivated mushroom in the world, Lentinula edodes, also known as shiitake (Japanese) or xiang-gu (Chinese). At present, nine species are recognized in the genus, based on morphology, mating criteria, and geographic distributions. However, analyses of internal transcribed spacers (ITS) of ribosomal RNA genes have suggested that there are cryptic lineages. We analyzed a global-scale phylogenetic dataset from 325 Lentinula individuals from 24 countries in Asia-Australasia and the Americas plus Madagascar, with 325 sequences of ITS, 80 LSU sequences, and 111 sequences of translation elongation factor (tef1-α) genes. We recovered 15 independent lineages (Groups 1-15) that may correspond to species. Lineages in Asia-Australasia (Groups 1-5) and the Americas plus Madagascar (Groups 6-15) formed sister clades. Four lineages are represented only by sequences from single individuals and require further molecular sampling, including L. aff. raphanica (Group 7), L. ixodes (Group 8), L. boryana (Group 12), and L. aff. aciculospora (Group 14). Groups 1 and 5 are here referred to L. edodes and L. aff. edodes, respectively. However, these groups most likely represent the same species and are only recognized as (unsupported) monophyletic lineages by maximum likelihood analyses of ITS alone. Other putative species resolved here include L. lateritia (Group 2), L. novae-zelandieae (Group 3), L. aff. lateritia (Group 4), L. raphanica (Group 6), L. aff. detonsa (Group 9), L. detonsa (Group 10), L. guzmanii sp. nov. (Group 11), L. aciculospora (Group 13), and L. madagasikarensis (Group 15). Groups 9-12 represent the "L. boryana complex". Molecular clock and historical biogeographic analyses suggest that the most recent common ancestor (MRCA) of Lentinula can be placed in the middle Oligocene, ca. 30 million years ago (Ma), and had a likely presence in neotropical America. The MRCA of Lentinula in the Americas and Madagascar lived ca. 22 Ma in the Neotropics and the MRCA of Lentinula in Asia-Australasia lived ca. 6 Ma in Oceania. Given the current knowledge about plate tectonics and paleoclimatic models of the last 30 Myr, our phylogenetic hypothesis suggests that the extant distribution of Lentinula is likely to have arisen, in large part, due to long-distance dispersal. Lentinula collections include at least four dubious taxa that need further taxonomic studies: L. reticeps from the USA (Ohio); L. guarapiensis from Paraguay; Lentinus puiggarii from Brazil (São Paulo); and "L. platinedodes" from Vietnam. Approximately ten of the fifteen Groups are reported on Fagaceae, which appears to be the ancestral substrate of Lentinula.

Physiological mechanisms by which gypsum increases the growth and yield of Lentinula edodes

Lentinula edodes is one of the most important commercially cultivated edible mushrooms. It is well known that gypsum (CaSO₄·2H₂O) supplementation in sawdust medium increases the yield of L. edodes, while the physiological mechanisms remain unclear. Our previous study showed that the acidification of the medium to pH 3.5-4.0 was essential for the growth of L. edodes. In this study, it was found that the oxalic acid excreted by L. edodes was responsible for the acidification of the medium. The biosynthesis of oxalic acid was regulated by the ambient pH and buffer capacity of the medium. To acidify the sawdust medium, the concentrations of total and soluble oxalate were 51.1 mmol/kg and 10.8 mmol/kg, respectively. However, when the concentration of soluble oxalate was 8.0 mmol/kg, the mycelial growth rate decreased by 29% compared with the control. Soluble oxalate was toxic to L. edodes, while soluble sulfate was nontoxic. CaSO₄ reacted with soluble oxalate to form nontoxic insoluble CaC₂O₄ and the strong acid H₂SO₄. When the CaSO₄ supplemented in sawdust medium was more than 25 mmol/kg, the soluble oxalate decreased to less than 1 mmol/kg, and the mycelial growth rate increased by 32% compared with the control. In conclusion, gypsum improved the growth and yield by relieving the toxicity of oxalate and facilitating the acidification of sawdust medium. KEY POINTS: • L. edodes excretes oxalic acid to acidify the ambient environment for growth. • Soluble oxalate is toxic to L. edodes. • Gypsum increases growth by reacting with oxalate to relieve its toxicity.

Oxidative Stress and Autophagy Are Important Processes in Post Ripeness and Brown Film Formation in Mycelium of Lentinula edodes

Lentinula edodes (Berk.) Pegler, the shiitake mushroom, is one of the most important mushrooms in the global mushroom industry. Although mycelium post ripeness and brown film (BF) formation are crucial for fruiting body initiation, the underlying molecular mechanisms of BF formation are largely unknown. In this study, proteomic quantification (relative and absolute) and metabolomic profiling of L. edodes were performed using isobaric tags and gas chromatography-mass spectroscopy, respectively. A total of 2,474 proteins were identified, which included 239 differentially expressed proteins. Notably, several proteins associated with autophagy were upregulated, including RPD3, TOR1, VAC8, VPS1, and VPS27. Transmission electron microscopy also indicated that autophagy occurred in post ripeness and BF formation. In time-dependent analysis of the metabolome, metabolites associated with oxidative stress and autophagy changed significantly, including mannitol, trehalose, myo-inositol, glucose, leucine, valine, glutamine, and 4-aminobutyric acid. Thus, oxidative stress and autophagy were important processes in post ripeness and BF formation in L. edodes, and new insights were gained into molecular mechanisms at proteome and metabolome levels.

Chromosomal genome and population genetic analyses to reveal genetic architecture, breeding history and genes related to cadmium accumulation in Lentinula edodes

BACKGROUND

Lentinula edodes (Berk.) is the second most productive mushroom in the world. It contains compounds effective for antiviral, antitumor, antioxidant and immune regulation. Although genomes have previously been reported for this species, a high-quality chromosome-level reference for L. edodes is unavailable. This hinders detailed investigation of population genetics, breeding history of strains and genes related to environmental stress responses.

RESULTS

A high-quality chromosome-level genome was constructed. We separated a monokaryon from protoplasts of the commercial L. edodes strain L808 and assembled the genome of L. edodes using PacBio long-read and Illumina short-read sequencing, along with the high-throughput chromatin conformation capture (Hi-C) technique. We assembled a 45.87 Mb genome, and 99% of the sequences were anchored onto 10 chromosomes. The contig and scaffold N50 length were 2.17 and 4.94 Mb, respectively. Over 96% of the complete Benchmarking Universal Single-Copy Orthologs (BUSCO) were identified, and 9853 protein-coding genes were predicted. We performed population genome resequencing using 34 wild strains and 65 commercial cultivars of L. edodes originating from China, Japan, the United States and Australia. Based on whole-genome variants, we showed substantial differences in the Chinese wild population, which divided into different branches according to the main areas of their geographical distribution. We also determined the breeding history of L. edodes at the molecular level, and demonstrated that the cultivated strains in China mainly originated from wild strains from China and Northeast Asia. Phenotypic analysis showed that 99 strains exhibited differences on the Cd accumulation. Three significant loci in the of L. edodes genome were identified using the genome-wide association study (GWAS) of Cd accumulation traits. Functional genes associated with Cd accumulation traits were related to DNA ligase and aminoacyl tRNA synthetase, indicating that DNA damage repair and in vivo protein translation may be responses to Cd stress.

CONCLUSIONS

A high-quality chromosome-level genome and population genetic data of L. edodes provide genetic resources for functional genomic, evolutionary and artificial breeding studies for L. edodes.

Haplotype-Resolved Genome Analyses Reveal Genetically Distinct Nuclei within a Commercial Cultivar of Lentinula edodes

Lentinula edodes is a tetrapolar basidiomycete with two haploid nuclei in each cell during most of their life cycle. Understanding the two haploid nuclei genome structures and their interactions on growth and fruiting body development has significant practical implications, especially for commercial cultivars. In this study, we isolated and assembled the two haploid genomes from a commercial strain of L. edodes using Illumina, HiFi, and Hi-C technologies. The total genome lengths were 50.93 Mb and 49.80 Mb for the two monokaryons SP3 and SP30, respectively, with each assembled into 10 chromosomes with 99.63% and 98.91% anchoring rates, respectively, for contigs more than 100 Kb. Genome comparisons suggest that two haploid nuclei likely derived from distinct genetic ancestries, with ~30% of their genomes being unique or non-syntenic. Consistent with a tetrapolar mating system, the two mating-type loci A (matA) and B (matB) of L. edodes were found located on two different chromosomes. However, we identified a new but incomplete homeodomain (HD) sublocus at ~2.8 Mb from matA in both monokaryons. Our study provides a solid foundation for investigating the relationships among cultivars and between cultivars and wild strains and for studying how two genetically divergent nuclei coordinate to regulate fruiting body formation in L. edodes.

Genetic structure and evolutionary diversity of mating-type (MAT) loci in Hypsizygus marmoreus

The mating compatibility in fungi is generally governed by genes located within a single or two unlinked mating type (MAT) loci. Hypsizygus marmoreus is an edible mushroom in the order Agaricales with a tetrapolar system, which contains two unlinked MAT loci-homeodomain (HD) transcription factor genes and pheromone/pheromone receptor genes (P/R). In this study, we analyzed the genetic structure and diversity of MAT loci in tetrapolar system of H. marmoreus through sequencing of 54 heterokaryon and 8 homokaryon strains. Although within the HD loci, the gene order was conserved, the gene contents were variable, and the HD loci haplotypes were further classified into four types. By analyzing the structure, phylogeny, and the HD transmissibility based on the progeny of these four HD mating-type loci types, we found that they were heritable and tightly linked at the HD loci. The P/R loci genes were found to comprise three pheromone receptors, three pheromones, and two pheromone receptor-like genes. Intra- and inter-specific phylogenetic analyses of pheromone receptors revealed that the STE3 genes were divided into three groups, and we thus theorize that they diverged before speciation. Comparative analysis of the MAT regions among 73 Basidiomycete species indicated that the diversity of HD and P/R loci in Agaricales and Boletales may contribute to mating compatibility. The number of HD genes were not correlated with the tetrapolar or bipolar systems. In H. marmoreus, the expression levels of these genes at HD and P/R loci of compatible strains were found higher than in those of homonuclear/homokaryotic strains, indicating that these mating genes acted as switches for mating processes. Further collinear analysis of HD loci in interspecific species found that HD loci contains conserved recombination hotspots showing major rearrangements in Coprinopsis cinerea and Schizophyllum commune, suggesting different mechanisms for evolution of physically linked MAT loci in these groups. It seems likely that gene rearrangements are common in Agaricales fungi around HD loci. Together, our study provides insights into the genomic basis of mating compatibility in H. marmoreus.

Metabolic Profiles, Bioactive Compounds, and Antioxidant Capacity in Lentinula edodes Cultivated on Log versus Sawdust Substrates

Lentinula edodes (shiitake) is a popular nutritious edible mushroom with a desirable aroma and flavor. Traditional cultivation of L. edodes on beds of logs has been replaced by cultivation on sawdust, but the effects of cultivation changes on L. edodes mushrooms have not been well characterized. We determined the metabolic profile, bioactive compounds, and antioxidant capacity in L. edodes grown on log or sawdust substrates. Metabolic profiles of L. edodes extracts were determined by ¹H nuclear magnetic resonance (NMR) and ultra-performance liquid chromatography to quadrupole time-of-flight mass spectrometry. Principal component analysis score plots from ¹H NMR analysis showed clear differences between samples. Concentrations of primary metabolites, especially amino acids, generally decreased in L. edodes grown on logs compared to sawdust. Phenolic compounds showed variations in concentration depending on the cultivation method. Bioactive compounds and their antioxidant capacity were analyzed spectrophotometrically. L. edodes cultivated on logs had high concentrations of bioactive compounds with strong antioxidant capacity compared to L. edodes cultivated on sawdust. Thus, the concentration of primary metabolites was high in L. edodes grown on sawdust, which produces a high growth rate. In contrast, log-cultivated L. edodes, which were similar to wild mushrooms, had high levels of bioactive compounds and high antioxidant capacity. This information is useful for determining optimal cultivation conditions for nutritional and medicinal uses of L. edodes mushrooms.

The Genomic and Transcriptomic Analyses of Floccularia luteovirens, a Rare Edible Fungus in the Qinghai-Tibet Plateau, Provide Insights into the Taxonomy Placement and Fruiting Body Formation

Floccularia luteovirens is a famous and precious edible mushroom (Huang Mogu) on the Qinghai–Tibet plateau that has a unique flavor and remarkable medical functions. Herein, we report a reference-grade 27 Mb genome of F. luteovirens containing 7068 protein-coding genes. The genome component and gene functions were predicted. Genome ontology enrichment and pathway analyses indicated the potential production capacity for terpenoids, polyketides and polysaccharides. Moreover, 16 putative gene clusters and 145 genes coding for secondary metabolites were obtained, including guadinomine and melleolides. In addition, phylogenetic and comparative genomic analyses shed light on the precise classification of F. luteovirens suggesting that it belongs to the genus Floccularia instead of Armillaria. RNA-sequencing and comparative transcriptomic analysis revealed differentially expressed genes during four developmental stages of F. luteovirens, that of which helps to identify important genes regulating fruiting body formation for strain modification. This study will provide insight into artificial cultivation and increase the production of useful metabolites.

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

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

Screening candidate genes related to volatile synthesis in shiitake mushrooms and construction of regulatory networks to effectively improve mushroom aroma

BACKGROUND

Metabolite formation is a physiological stress response during the growth and development of shiitake mushrooms (Lentinula edodes). The characteristic flavor metabolites are important quality components in shiitake mushrooms. To investigate the formation mechanisms of characteristic flavor metabolites, transcriptome analyses were performed on shiitake mushrooms harvested at different growth stages.

RESULTS

In total, 30 genes related to the synthesis of characteristic volatiles of mushrooms were identified via screening. Through KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of the selected genes and correlation analyses of gene expressions, the main volatile synthesis pathways were determined as histidine metabolism, glutathione metabolism and biosynthesis of unsaturated fatty acids. Gene cluster and correlation analyses were performed to clarify the combined effects of different genes in the enzymatic reactions. Further, a correlation network of candidate genes was built based on the gene expression levels.

CONCLUSION

The activities of flavor synthases and the content of characteristic flavor metabolites were analyzed; the enzyme activity changes and metabolic product distribution sites were clarified. A synthesis and regulation network was constructed for the candidate genes and characteristic volatiles, and information was obtained for 16 hub genes. Moreover, it was essential to identify and characterize the key genes and synthases involved in the synthesis of the characteristic volatiles of mushrooms. This information provides us with a better understanding of the biosynthesis and regulation of the volatiles, which will lay the foundation for improving the quality of shiitake mushrooms. © 2021 Society of Chemical Industry.

Haploid Genome Analysis Reveals a Tandem Cluster of Four HSP20 Genes Involved in the High-Temperature Adaptation of Coriolopsis trogii

Coriolopsis trogii is a typical thermotolerant basidiomycete fungus, but its thermotolerance mechanisms are currently unknown. In this study, two monokaryons of C. trogii strain Ct001 were assembled: Ct001_29 had a genome assembly size of 38.85 Mb and encoded 13,113 genes, while Ct001_31 was 40.19 Mb in length and encoded 13,309 genes. Comparative intra- and interstrain genomic analysis revealed the rich genetic diversity of C. trogii, which included more than 315,194 single-nucleotide polymorphisms (SNPs), 30,387 insertion/deletions (indels), and 1,460 structural variations. Gene family analysis showed that the expanded families of C. trogii were functionally enriched in lignocellulose degradation activities. Furthermore, a total of 14 allelic pairs of heat shock protein 20 (HSP20) genes were identified in the C. trogii genome. The expression profile obtained from RNA sequencing (RNA-Seq) showed that four tandem-duplicated allelic pairs, HSP20.5 to HSP20.8, had more than 5-fold higher expression at 35°C than at 25°C. In particular, HSP20.5 and HSP20.8 were the most highly expressed HSP20 genes. Allelic expression bias was found for HSP20.5 and HSP20.8; the expression of Ct29HSP20.8 was at least 1.34-fold higher than that of Ct31HSP20.8, and that of Ct31HSP20.5 was at least 1.5-fold higher than that of Ct29HSP20.5. The unique structural and expression profiles of the HSP20 genes revealed by these haplotype-resolved genomes provide insight into the molecular mechanisms of high-temperature adaptation in C. trogii. IMPORTANCE Heat stress is one of the most frequently encountered environmental stresses for most mushroom-forming fungi. Currently available fungal genomes are mostly haploid because high heterozygosity hinders diploid genome assembly. Here, two haplotype genomes of C. trogii, a thermotolerant basidiomycete, were assembled separately. A conserved tandem cluster of four HSP20 genes showing allele-specific expression was found to be closely related to high-temperature adaptation in C. trogii. The obtained haploid genomes and their comparison offer a more thorough understanding of the genetic background of C. trogii. In addition, the responses of HSP20 genes at 35°C, which may contribute to the growth and survival of C. trogii at high temperatures, could inform the selection and breeding of elite strains in the future.

Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes

•

We detected three subgroups of L. edodes with robust phenotypic differentiation.

•

The three subgroups were diverged 36,871 generations ago.

•

L. edodes cultivars in China might originate from the vicinity of Northeast China.

•

We dissected the genetic basis of adaptive evolution in L. edodes.

•

Genes related to fruiting body development are involved in adaptive evolution.

Introduction

Mushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi.

Objectives

This study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution.

Methods

We analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution.

Results

We detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation.

Conclusion

Adaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

Complete Genomic Characterization and Identification of Saccharomycopsisphalluae sp. nov., a Novel Pathogen Causes Yellow Rot Disease on Phallus rubrovolvatus

“Hongtuozhusun” (Phallus rubrovolvatus) is an important edible and medicinal mushroom endemic to Southwest China. However, yellow rot disease is a severe disease of P. rubrovolvatus that occurs extensively in Guizhou Province. It has caused major economic losses and hinders the development of the P. rubrovolvatus industry. In this study, 28 microorganism strains were isolated from diseased fruiting bodies of P. rubrovolvatus at various stages, two of which were confirmed to be pathogenic based on Koch’s postulates. These two strains are introduced herein as Saccharomycopsisphalluae sp. nov. based on morphological, physiological, and molecular analysis. We reported a high-quality de novo sequencing and assembly of the S. phalluae genome using single-molecule real-time sequencing technology. The whole genome was approximately 14.148 Mb with a G+C content of 43.55%. Genome assembly generated 8 contigs with an N50 length of 1,822,654 bp. The genome comprised 5966 annotated protein-coding genes. This is the first report of mushroom disease caused by Saccharomycopsis species. We expect that the information on genome properties, particularly in pathogenicity-related genes, assist in developing effective control measures in order to prevent severe losses and make amendments in management strategies.

Genomic and Experimental Investigations of Auriscalpium and Strobilurus Fungi Reveal New Insights into Pinecone Decomposition

Saprophytic fungi (SPF) play vital roles in ecosystem dynamics and decomposition. However, because of the complexity of living systems, our understanding of how SPF interact with each other to decompose organic matter is very limited. Here we studied their roles and interactions in the decomposition of highly specialized substrates between the two genera Auriscalpium and Strobilurus fungi-colonized fallen pinecones of the same plant sequentially. We obtained the genome sequences from seven fungal species with three pairs: A. orientale-S. luchuensis, A. vulgare-S. stephanocystis and A. microsporum-S. pachcystidiatus/S. orientalis on cones of Pinus yunnanensis, P. sylvestris and P. armandii, respectively, and the organic profiles of substrate during decomposition. Our analyses revealed evidence for both competition and cooperation between the two groups of fungi during decomposition, enabling efficient utilization of substrates with complementary profiles of carbohydrate active enzymes (CAZymes). The Auriscalpium fungi are highly effective at utilizing the primary organic carbon, such as lignin, and hemicellulose in freshly fallen cones, facilitated the invasion and colonization by Strobilurus fungi. The Strobilurus fungi have genes coding for abundant CAZymes to utilize the remaining organic compounds and for producing an arsenal of secondary metabolites such as strobilurins that can inhibit other fungi from colonizing the pinecones.

Effect of Pretreated Colza Straw on the Growth and Extracellular Ligninolytic Enzymes Production by Lentinula edodes and Ganoderma lucidum

Lentinula edodes 3565 and Ganoderma lucidum 9621 were compared for their ability to produce lignocellulolytic enzymes in submerged (SM) and surface liquid (SL) fermentation of hydrolysed colza straw lignin waste that remained after the production of furfural and bioethanol (CS lignin). Application of cultivated mushrooms to dispose of pretreated colza straw agricultural waste is an approach to decrease the quantity of residual lignin while simultaneously obtaining active substances, e.g., the ligninolytic enzyme complex from mycelium. The effect of adding CS lignin to culture media on the yield of L. edodes and G. lucidum mycelium and extracellular laccase activity was studied. It was revealed that the mycelial growth of G. lucidum on solid media was significantly improved by adding CS lignin. Laccase activity during SL cultivation of L. edodes on medium with CS lignin gradually increased over the experiment starting on day 21 and peaked at 520 U/mL on day 28. G. lucidum expressed the maximum laccase activity, 540 U/mL, during the first 14 days of mycelium SM cultivation. Extracellular laccase activity was enhanced about 35- to 40-fold at cultivation of L. edodes and about 10- to 15-fold in the case of G. lucidum by supplementing liquid culture media with CS lignin.

Transcriptional response of mushrooms to artificial sun exposure

Climate change causes increased tree mortality leading to canopy loss and thus sun-exposed forest floors. Sun exposure creates extreme temperatures and radiation, with potentially more drastic effects on forest organisms than the current increase in mean temperature. Such conditions might potentially negatively affect the maturation of mushrooms of forest fungi. A failure of reaching maturation would mean no sexual spore release and, thus, entail a loss of genetic diversity. However, we currently have a limited understanding of the quality and quantity of mushroom-specific molecular responses caused by sun exposure. Thus, to understand the short-term responses toward enhanced sun exposure, we exposed mushrooms of the wood-inhabiting forest species Lentinula edodes, while still attached to their mycelium and substrate, to artificial solar light (ca. 30°C and 100,000 lux) for 5, 30, and 60 min. We found significant differentially expressed genes at 30 and 60 min. Eukaryotic Orthologous Groups (KOG) class enrichment pointed to defense mechanisms. The 20 most significant differentially expressed genes showed the expression of heat-shock proteins, an important family of proteins under heat stress. Although preliminary, our results suggest mushroom-specific molecular responses to tolerate enhanced sun exposure as expected under climate change. Whether mushroom-specific molecular responses are able to maintain fungal fitness under opening forest canopies remains to be tested.

De novo Sequencing and Comparative Transcriptome Analyses Provide First Insights Into Polysaccharide Biosynthesis During Fruiting Body Development of Lentinula edodes

Polysaccharides separated from Lentinula edodes are well known for their medicinal properties. However, the precise molecular mechanisms of polysaccharide biosynthesis in L. edodes remain unclear. In this study, the fruiting bodies of L. edodes in four developmental stages with significant differences in polysaccharide yield were collected, and the characteristics of polysaccharides were studied. De novo sequencing and comparative transcriptomic analysis were performed by using high-throughput Illumina RNA-sequencing. KS1P30, KS2P30, KS3P30, and KS4P30 were obtained from the four developmental stages, respectively, by hot water extraction and 30% ethanol precipitation. These four polysaccharides had good immune activity in vitro; all of them were β-glucopyranose with a high molecular weight. Glucose was the main monosaccharide component of these polysaccharides. High-quality clean reads (57.88, 53.17, 53.28, and 47.56 million for different growth stages) and mapping ratios ranging from 84.75 to 90.11% were obtained. In total, 11,493 (96.56%) unigenes and 18,924 (97.46%) transcripts were successfully annotated in five public databases. The biosynthetic pathway and related genes of LEFP30 were mined. The molecular mechanism of LEFP30 yield change in the different developmental stages was predicted. The results provide some insights into the possible mechanisms involved in the biosynthetic pathway of this kind of polysaccharide in L. edodes fruiting bodies. They also indicate that candidate genes can be used as important resources for biotechnology and molecular breeding to regulate L. edodes fruiting body polysaccharide biosynthesis.

Genomic Analysis of Sarcomyxa edulis Reveals the Basis of Its Medicinal Properties and Evolutionary Relationships

Yuanmo [Sarcomyxa edulis (Y.C. Dai, Niemelä & G.F. Qin) T. Saito, Tonouchi & T. Harada] is an important edible and medicinal mushroom endemic to Northeastern China. Here we report the de novo sequencing and assembly of the S. edulis genome using single-molecule real-time sequencing technology. The whole genome was approximately 35.65 Mb, with a G + C content of 48.31%. Genome assembly generated 41 contigs with an N50 length of 1,772,559 bp. The genome comprised 9,364 annotated protein-coding genes, many of which encoded enzymes involved in the modification, biosynthesis, and degradation of glycoconjugates and carbohydrates or enzymes predicted to be involved in the biosynthesis of secondary metabolites such as terpene, type I polyketide, siderophore, and fatty acids, which are responsible for the pharmacodynamic activities of S. edulis. We also identified genes encoding 1,3-β-glucan synthase and endo-1,3(4)-β-glucanase, which are involved in polysaccharide and uridine diphosphate glucose biosynthesis. Phylogenetic and comparative analyses of Basidiomycota fungi based on a single-copy orthologous protein indicated that the Sarcomyxa genus is an independent group that evolved from the Pleurotaceae family. The annotated whole-genome sequence of S. edulis can serve as a reference for investigations of bioactive compounds with medicinal value and the development and commercial production of superior S. edulis varieties.

Sequencing and Comparative Genomic Analysis of a Highly Metal-Tolerant Penicillium janthinellum P1 Provide Insights Into Its Metal Tolerance

Heavy metal pollution is a global knotty problem and fungi hold promising potential for the remediation of wastewater containing heavy metals. Here, a new highly chromium-tolerance species, Penicillium janthinellum P1, is investigated. The genome of P1 was sequenced and assembled into 30 Mb genome size containing 10,955 predicted protein-coding genes with a GC content of 46.16% through an integrated method of Illumina short-read sequencing and single-molecule real-time Pacific Biosciences sequencing platforms. Through a phylogenetic analysis with model species of fungi, the evolutionary divergence time of Penicillium janthinellum P1 and Penicillium oxalicum 114-2 was estimated to be 74 MYA. 33 secondary metabolism gene clusters were identified via antiSMASH software, mainly including non-ribosomal peptide synthase genes and T1 polyketide synthase genes. 525 genes were annotated to encode enzymes that act on carbohydrates, involving 101 glucose-degrading enzymes and 24 polysaccharide synthase. By whole-genome sequence analysis, large numbers of metal resistance genes were found in strain P1. Especially ABC transporter and Superoxide dismutase ensure that the P1 fungus can survive in a chromium-polluted environment. ChrA and ChrR were also identified as key genes for chromium resistance. Analysis of their genetic loci revealed that the specific coding-gene arrangement may account for the fungus’s chromium resistance. Genetic information and comparative analysis of Penicillium janthinellum are valuable for further understanding the mechanism of high resistance to heavy metal chromium, and gene loci analysis provides a new perspective for identifying chromium-resistant strains.

RNA-Seq-based high-resolution linkage map reveals the genetic architecture of fruiting body development in shiitake mushroom, Lentinula edodes

•

We constructed a reference genetic map of Lentinula edodes.

•

We re-assembled a chromosome-level genome of L. edodes.

•

We disclosed three hotspots regions for fruiting body-related traits in shiitake.

•

We scanned candidate genes for fruiting body-related traits.

Fruiting body development (FBD) of mushroom-forming fungi has attracted tremendous interest. However, the genetic and molecular basis of FBD is poorly known. Here, using Lentinula edodes (shiitake) as a model, we deciphered the genetic architecture underlying fruiting body-related traits (FBRTs) by combined genomic, genetic and phenotypic data. Using RNA-Seq of fruiting bodies from 110 dikaryons in a bi-parental mapping population, we constructed an ultra-high-density genetic map of L. edodes (Lemap2.0) with a total length of 810.14 cM, which covered 81.7% of the shiitake genome. A total of 94 scaffolds of the shiitake genome were aligned to Lemap2.0 and re-anchored into nine pseudo-chromosomes. Then via quantitative trait locus (QTL) analysis, we disclosed an outline of the genetic architecture of FBD in shiitake. Twenty-nine QTLs and three main genomic regions associated with FBD of shiitake were identified. Using meta-QTL analysis, seven pleiotropic QTLs for multiple traits were detected, which contributed to the correlations of FBRTs. In the mapped QTLs, the expression of 246 genes were found to significantly correlate with the phenotypic traits. Thirty-three of them were involved in FBD and could represent candidate genes controlling the shape and size of fruiting bodies. Collectively, our findings have advanced our understanding of the genetic regulation of FBD in shiitake and mushroom-forming fungi at large.

Comparative Analysis of Carbohydrate Active Enzymes in the Flammulina velutipes var. lupinicola Genome

The purpose of this study was to determine the genome sequence of Flammulina velutipes var. lupinicola based on next-generation sequencing (NGS) and to identify the genes encoding carbohydrate-active enzymes (CAZymes) in the genome. The optimal assembly (71 kmer) based on ABySS de novo assembly revealed a total length of 33,223,357 bp (49.53% GC content). A total of 15,337 gene structures were identified in the F. velutipes var. lupinicola genome using ab initio gene prediction method with Funannotate pipeline. Analysis of the orthologs revealed that 11,966 (96.6%) out of the 15,337 predicted genes belonged to the orthogroups and 170 genes were specific for F. velutipes var. lupinicola. CAZymes are divided into six classes: auxiliary activities (AAs), glycosyltransferases (GTs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), glycoside hydrolases (GHs), and carbohydrate-binding modules (CBMs). A total of 551 genes encoding CAZymes were identified in the F. velutipes var. lupinicola genome by analyzing the dbCAN meta server database (HMMER, Hotpep, and DIAMOND searches), which consisted of 54-95 AAs, 145-188 GHs, 55-73 GTs, 6-19 PLs, 13-59 CEs, and 7-67 CBMs. CAZymes can be widely used to produce bio-based products (food, paper, textiles, animal feed, and biofuels). Therefore, information about the CAZyme repertoire of the F. velutipes var. lupinicola genome will help in understanding the lignocellulosic machinery and in-depth studies will provide opportunities for using this fungus for biotechnological and industrial applications.

Identification and Functional Characterization of a Novel Immunomodulatory Protein From Morchella conica SH

A novel fungal immunomodulatory protein (FIP) was found in the precious medical and edible mushroom Morchella conica SH, defined as FIP-mco, which belongs to the FIP family. Phylogenetic analyses of FIPs from different origins were performed using Neighbor-Joining method. It was found that FIP-mco belonged to a new branch of the FIP family and may evolved from a different ancestor compared with most other FIPs. The cDNA sequence of FIP-mco was cloned and expressed in the yeast Pichia Pastoris X33. The recombinant protein of FIP-mco (rFIP-mco) was purified by agarose Ni chromatography and determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and Western blot analysis. The protein rFIP-mco could significantly suppress the proliferation of A549 and HepG2 cells at the concentration of 15 and 5 μg/ml, respectively, and inhibited the migration and invasion of human A549 and HepG2 cells at the concentration of 15 and 30 μg/ml respectively in vitro. Further, rFIP-mco can significantly reduce the expression levels of TNF-α, IL-1β, and IL-6 in the THP1 cells (human myeloid leukemia mononuclear cells). In order to explore the potential mechanism of the cytotoxicity effect of rFIP-mco on A549 and HepG2 cells, cell cycle and apoptosis assay in the two cancer cells were conducted. The results demonstrated that G0/G1 to S-phase arrest and increased apoptosis may contribute to the proliferation inhibition by rFIP-mco in the two cancer cells. Molecular mechanism of rFIP-mco's reduction effect on the inflammatory cytokines was also studied by suppression of the NF-κB signaling pathway. It showed that suppression of NF-κB signaling is responsible for the reduction of inflammatory cytokines by rFIP-mco. The results indicated the prospect of FIP-mco from M. conica SH as an effective and feasible source for cancer therapeutic studies and medical applications.

Functional Analysis of a Novel ABL (Abnormal Browning Related to Light) Gene in Mycelial Brown Film Formation of Lentinula edodes

Lentinula edodes is a globally important edible mushroom species that is appreciated for its medicinal properties as well as its nutritional value. During commercial cultivation, a mycelial brown film forms on the surface of the sawdust growth medium at the late vegetative stage. Mycelial film formation is a critical developmental process that contributes to the quantity and quality of the mushroom yield. However, little is known regarding the genetic underpinnings of brown film formation on the surface of mycelial tissue. A novel causal gene associated with the formation of the mycelial brown film, named ABL (Abnormal browning related to light), was identified in this study. The comparative genetic analysis by dihybrid crosses between normal and abnormal browning film cultivars demonstrated that a single dominant allele was responsible for the abnormal mycelium browning phenotype. Whole-genome sequencing analysis of hybrid isolates revealed five missense single-nucleotide polymorphisms (SNPs) in the ABL locus of individuals forming abnormal partial brown films. Additional whole-genome resequencing of a further 16 cultivars showed that three of the five missense SNPs were strongly associated with the abnormal browning phenotype. Overexpression of the dominant abl-D allele in a wild-type background conferred the abnormal mycelial browning phenotype upon transformants, with slender hyphae observed as a general defective mycelial growth phenotype. Our methodology will aid the future discovery of candidate genes associated with favorable traits in edible mushrooms. The discovery of a novel gene, ABL, associated with mycelial film formation will facilitate marker-associated breeding in L. edodes.

Transcriptomic Analysis of Two Lentinula edodes Genotypes With Different Cadmium Accumulation Ability

Lentinula edodes, also known as Xiang'gu, is commonly eaten in cultures around the world. However, L. edodes is particularly susceptible to enrichment with heavy metals, particularly cadmium (Cd), which is toxic to human health. Understanding the molecular mechanism and mining key genes involved in Cd enrichment will facilitate genetic modification of L. edodes strains. Two L. edodes genotypes, Le4625 (with higher Cd enrichment capability) and Le4606 (with lower Cd enrichment capability) were used in this study. The Cd concentrations in the mycelia of the tested genotypes differed significantly after Cd (0.1 mg/L) exposure; and the Cd content in Le4625 (1.390 ± 0.098 mg/kg) was approximately three-fold that in Le4606 (0.440 ± 0.038 mg/kg) after 7 h of Cd exposure. A total of 24,592 transcripts were assessed by RNA-Seq to explore variance in Cd accumulation. Firstly, differentially expressed genes (DEGs) were analyzed separately following Cd exposure. In comparison with Ld4625, Ld4606 showed a greater number of Cd-induced changes in transcription. In Ld4606, DEGs following Cd exposure were associated with transmembrane transport, glutathione transfer and cytochrome P450, indicating that these genes could be involved in Cd resistance in L. edodes. Next, Le4606 and Le4625 were exposed to Cd, after which DEGs were identified to explore genetic factors affecting Cd accumulation. After Cd exposure, DEGs between Le4606 and Le4625 encoded proteins involved in multiple biological pathways, including transporters on the membrane, cell wall modification, oxidative stress response, translation, degradation, and signaling pathways. Cadmium enrichment in cells may activate MAPK signaling and the anti-oxidative stress response, which can subsequently alter signal transduction and the intracellular oxidation/reduction balance. Furthermore, several possible candidate genes involved in the Cd accumulation were identified, including the major facilitator superfamily genes, heat shock proteins, and laccase 11, a multicopper oxidase. This comparison of the transcriptomes of two L. edodes strains with different capacities for Cd accumulation provides valuable insight into the cultivation of mushrooms with less Cd enrichment and also serves as a reference for the construction of engineered strains for environmental pollution control.

Multilocus phylogeny- and fruiting feature-assisted delimitation of European Cyclocybe aegerita from a new Asian species complex and related species

Cyclocybe aegerita (synonym: Agrocybe aegerita) is a widely cultivated edible and reportedly almost cosmopolitan mushroom species that serves as a model fungus for basidiome formation and as producer of useful natural products and enzymes. Focusing on strains from different continents, here, we present a phylogenetic analysis of this species and some adjacent taxa that employs four phylogenetic markers. In addition, we tested the strains’ capability to fructify on agar media. Our analysis reveals that “C. aegerita sensu lato” splits up into the following two well-supported monophyletic geographic lineages: a European clade and an Asian clade. The European one is closely associated with the Chinese species Cyclocybe salicaceicola. In contrast, the Asian lineage, which we preliminarily designate as Cyclocybe chaxingu agg., may comprise several species (species complex) and clusters with the Pacific species Cyclocybe parasitica (New Zealand). In addition, fruiting properties differ across C. aegerita and its Asian and Pacific relatives; however, strains from the Asian clade and C. parasitica tend to form larger basidiomes with relatively big caps and long stipes and strains from the European clade exhibit a more variable fruiting productivity with the tendency to form more basidiomes, with smaller caps and shorter stipes. Moreover, some strains showed individual fruiting patterns, such as the preference to fruit where they were exposed to injuring stimuli. In conclusion, the delimitation of the newly delimited Asian species complex from our multilocus phylogeny of “C. aegerita sensu lato”, which is supported by phenotypic data, depicts an exemplary case of biogeographic diversity within a previously thought homogeneous species of near worldwide distribution.

Draft genomic sequence of Armillaria gallica 012m: insights into its symbiotic relationship with Gastrodia elata

Armillaria species (Basidiomycota, Physalacriaceae) are well known as plant pathogens related to serious root rot disease on various trees in forests and plantations. Interestingly, some Armillaria species are essential symbionts of the rare Chinese medicinal herb Gastrodia elata, a rootless and leafless orchid used for over 2000 years. In this work, an 87.3-M draft genome of Armillaria gallica 012m strain, which was symbiotic with G. elata, was assembled. The genome includes approximately 23.6% repetitive sequences and encodes 26,261 predicted genes. In comparison with other four genomes of Armillaria, the following gene families related to pathogenicity/saprophytic phase, including cytochrome P450 monooxygenases, carbohydrate-active enzyme AA3, and hydrophobins, were significantly contracted in A. gallica 012m. These characteristics may be beneficial for G. elata to get less injuries. The genome-guided analysis of differential expression between rhizomorph (RH) and vegetative mycelium (VM) showed that a total of 2549 genes were differentially expressed, including 632 downregulated genes and 1917 upregulated genes. In the RH, most differentially expressed genes (DEGs) related to pathogenicity were significantly upregulated. To further elucidate gene function, Gene Ontology enrichment analysis showed that the upregulated DEGs significantly grouped into monooxygenase activity, hydrolase activity, glucosidase activity, extracellular region, fungal cell wall, response to xenobiotic stimulus, response to toxic substance, etc. These phenomena indicate that RH had better infection ability than VM. The infection ability of RH may be beneficial for G. elata to obtain nutrition, because the rhizomorph constantly infected the nutritional stems of G. elata and formed the hyphae that can be digested by G. elata. These results clarified the characteristics of A. gallica 012m and the reason why the strain 012m can establish a symbiotic relationship with G. elata in some extent from the perspective of genomics.