Population genomic analysis uncovers environmental stress-driven selection and adaptation of Lentinula edodes population in China

Population genetic structure of Hymenopellis radicata germplasm resources based on genome re-sequencing

Through whole-genome re-sequencing of 18 Hymenopellis radicata germplasm resources collected from diverse regions in China, we identified significant variations in the form of Single Nucleotide Polymorphisms (SNPs) and Insertions and Deletions (InDels). These variations were comprehensively annotated, shedding light on the mutation types present in the entire genome of the H. radicata germplasm. This analysis revealed the number and position information of each mutation and provided insights into the overall genomic landscape of H. radicata germplasm. Utilizing SNP data, we delved into the population structure of the 18 H. radicata germplasm resources. The results indicated the presence of 2,335,179 Indel sites and 12,050,448 SNP sites. The population structure analysis unveiled two distinct subgroups among the H. radicata germplasm resources. Phenotypic statistics, principal component analysis, and phylogenetic tree results echoed the findings of the population structure analysis. Different strains of H. radicata from various regions in China exhibited notable differences in genetic diversity, mycelial growth rate, yield, and fruiting body characteristics. Significant disparities were observed between the two subgroups, while strains within each subgroup shared common characteristics. This research establishes a solid foundation for integrating H. radicata into diverse breeding programs. The data underscore the potential of H. radicata for genetic improvement and exploitation in breeding initiatives, paving the way for future advancements in this field.

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.

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.

Prevalence and diversity of mycoviruses occurring in Chinese Lentinula edodes germplasm resource

Incidence and banding patterns of virus-like dsRNA elements in 215 Chinese genetically diverse Lentinula edodes strains collected from wide geographic distribution (or producing areas) were first investigated, and 17 viruses were identified including eight novel viruses. The results revealed a 63.3% incidence of dsRNA elements in the cultivated strains and a 67.2% incidence in the wild strains. A total of 10 distinguishable dsRNAs ranging from 0.6 to 12 kbp and 12 different dsRNA patterns were detected in the positive strains. The molecular information of these dsRNA elements was characterized, and the molecular information of the other 12 different viral sequences with (+) ssRNA genome was revealed in four L. edodes strains with complex dsRNA banding patterns. RT-PCR was also done to verify the five dsRNA viruses and 12 (+) ssRNA ones. The results presented may enrich our understanding of L. edodes virus diversity, and will promote further research on virus-host interactions. IMPORTANCE: Viral infections involve complicated interactions including benign, harmful or possibly beneficial to hosts. Sometimes environment could lead to a transition in lifestyles from persistent to acute, resulting in a disease phenotype. The quality of spawn, such as the vulnerability to infection of viruses, is therefore important for mushroom production. Lentinula edodes, a wood rot basidiomycete fungus, was widely cultivated in the world for its edible and medicinal properties. In this study, the profile of dsRNA elements from Chinese genetically diverse L. edodes strains collected from wide geographic distribution or producing areas was first investigated. The molecular information of the dsRNA elements was characterized. Additionally, 12 different viral sequences with (+) ssRNA genome from four L. edodes strains with complex dsRNA banding patterns were identified. The results presented here will broaden our knowledge about mushroom viruses, and promote further studies of L. edodes production and the interaction between viruses and L. edodes.

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.

Growth modulation and metabolic responses of Ganoderma boninense to salicylic acid stress

Various phenolic compounds have been screened against Ganoderma boninense, the fungal pathogen causing basal stem rot in oil palms. In this study, we focused on the effects of salicylic acid (SA) on the growth of three G. boninense isolates with different levels of aggressiveness. In addition, study on untargeted metabolite profiling was conducted to investigate the metabolomic responses of G. boninense towards salicylic acid. The inhibitory effects of salicylic acid were both concentration- (P < 0.001) and isolate-dependent (P < 0.001). Also, growth-promoting effect was observed in one of the isolates at low concentrations of salicylic acid where it could have been utilized by G. boninense as a source of carbon and energy. Besides, adaptation towards salicylic acid treatment was evident in this study for all isolates, particularly at high concentrations. In other words, inhibitory effect of salicylic acid treatment on the fungal growth declined over time. In terms of metabolomics response to salicylic acid treatment, G. boninense produced several metabolites such as coumarin and azatyrosine, which suggests that salicylic acid modulates the developmental switch in G. boninense towards the defense mode for its survival. Furthermore, the liquid chromatography time-of-flight mass spectrometry (LC-TOF-MS) analysis showed that the growth of G. boninense on potato dextrose agar involved at least four metabolic pathways: amino acid metabolism, lipid pathway, tryptophan pathway and phenylalanine pathway. Overall, there were 17 metabolites that contributed to treatment separation, each with P<0.005. The release of several antimicrobial metabolites such as eudistomin I may enhance G. boninense's competitiveness against other microorganisms during colonisation. Our findings demonstrated the metabolic versatility of G. boninense towards changes in carbon sources and stress factors. G. boninense was shown to be capable of responding to salicylic acid treatment by switching its developmental stage.

Construction of a CRISPR/Cas9-Mediated Genome Editing System in Lentinula edodes

CRISPR/Cas9 genome editing systems have been established in a broad range of eukaryotic species. Herein, we report the first method for genetic engineering in pyogo (shiitake) mushrooms (Lentinula edodes) using CRISPR/Cas9. For in vivo expression of guide RNAs (gRNAs) targeting the mating-type gene HD1 (LeA1), we identified an endogenous LeU6 promoter in the L. edodes genome. We constructed a plasmid containing the LeU6 and glyceraldehyde-3-phosphate dehydrogenase (LeGPD) promoters to express the Cas9 protein. Among the eight gRNAs we tested, three successfully disrupted the LeA1 locus. Although the CRISPR-Cas9-induced alleles did not affect mating with compatible monokaryotic strains, disruption of the transcription levels of the downstream genes of LeHD1 and LeHD2 was detected. Based on this result, we present the first report of a simple and powerful genetic manipulation tool using the CRISPR/Cas9 toolbox for the scientifically and industrially important edible mushroom, L. edodes.

Molecular identification of Shiitake [Lentinula edodes Berk (Pegler)] and production of secondary metabolites with biotechnological potential

The Shitake mushroom (Lentinula edodes) is the second most-consumed mushroom in the world; in Colombia, it is cultivated and commercialized on a small scale in some supermarkets. Little is known about the precedence, nutritional and medicinal properties of Shiitake produced in Colombia. In this study, four shiitake isolates were grown in Colombia (LEUCO1, LEUCO2, LEUCO3, and LEUCO4) were sequenced in their ITS genes and evaluated for the production of three medicinal metabolites, eritadenine, ergotioneine and β-glucans (1,3-1, 6), using submerged culture. Genetic analysis revealed that all the isolates were close and related to the Japanese strain Cr62. LEUCO1 and LEUCO2 showed a distance of 0.000, as well as LEUCO3 and LEUCO4. All four isolates produced erythadenin in a range of 26.3-8.6 mg / L, with the best performance of LEUCO1 at 26.3 mg / L (p <0.05). Ergotioneine was produced with similar statistical yields in all the isolates with an average of 0.7 mg / g of dry weight biomass (DW). Β-glucans (1.3-1.6) were produced with yields of 5.6 - 3.8% of DW biomass, with the best values ​​for LEUCO2 and the lowest for LEUCO4 (p <0.05). In conclusion, we identified low genetic diversity in the four isolates, corresponding to two haplotypes with minimal genetic difference between them, related to the Japanese strain Cr62, indicating that Colombian farmers cultivate almost the same strains of shitake. Secondary metabolites, eritadenine, β-glucans and ergotioneine were found in promising yields useful for the pharmaceutical and food industries. More studies should be conducted to improve the yield of shitake metabolites through new growing conditions for industrial production and to find metabolic pathways and related genes.

Analysis of Genetic Diversity and Population Structure of Wild Strains and Cultivars Using Genomic SSR Markers in Lentinula edodes

In this study, the genetic diversity and the population structure of 77 wild strains and 23 cultivars of Lentinula edodes from Korea were analyzed using 20 genomic SSRs, and their genetic relationship was investigated. The tested strains of L. edodes were divided into three sub-groups consisting of only wild strains, mainly wild strains and several cultivars, and mainly cultivars and several wild strains by distance-based analysis. Using model-based analysis, L. edodes strains were divided into two subpopulations; the first one consisting of only wild strains and the second one with mainly cultivars and several wild strains. Moreover, AMOVA analysis revealed that the genetic variation in the cultivars was higher than that in the wild strains. The expected and observed heterozygosity and values indicating the polymorphic information content of L. edodes cultivars from Korea were also higher than that of the wild strains. Based on these results, we presume that the cultivars in Korea have developed by using numerous strains from other countries. In conclusion, the usage of wild strains for the development of new cultivars could improve the adaptability of L. edodes to biotic and abiotic stress.

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

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

Mitogenome types of two Lentinula edodes sensu lato populations in China

China has two populations of Lentinula edodes sensu lato as follows: L. edodes sensu stricto and an unexcavated morphological species respectively designated as A and B. In a previous study, we found that the nuclear types of the two populations are distinct and that both have two branches (A1, A2, B1 and B2) based on the internal transcribed spacer 2 (ITS2) sequence. In this paper, their mitogenome types were studied by resequencing 20 of the strains. The results show that the mitogenome type (mt) of ITS2-A1 was mt-A1, that of ITS2-A2 was mt-A2, and those of ITS2-B1 and ITS2-B2 were mt-B. The strains with heterozygous ITS2 types had one mitogenome type, and some strains possessed a recombinant mitogenome. This indicated that there may be frequent genetic exchanges between the two populations and both nuclear and mitochondrial markers were necessary to identify the strains of L. edodes sensu lato. In addition, by screening SNP diversity and comparing four complete mitogenomes among mt-A1, mt-A2 and mt-B, the cob, cox3, nad2, nad3, nad4, nad5, rps3 and rrnS genes could be used to identify mt-A and mt-B and that the cox1, nad1 and rrnL genes could be used to identify mt-A1, mt-A2 and mt-B.

Effects of Medium Composition and Genetic Background on Agrobacterium-Mediated Transformation Efficiency of Lentinula edodes

The establishment of genetic transformation method is crucial for the functional genomics research in filamentous fungi. Although the transformation method has been developed in several types of fungi, a highly efficient and convenient transformation system is desperately needed in Lentinula edodes. Present work established the Agrobacterium-mediated transformation (ATMT) of basidiomycete L. edodes in both monokaryon and dikaryon mycelia by using constructed binary plasmid pCAMBIA-1300-GFP. Then, the transformation efficiency of ATMT was evaluated by using different mediums for recipient incubation and different varieties of L. edodes. The results showed that in dikaryon strain W1, the positive hygromycin-resistant transformants was observed in all medium with the positive frequency of selected transformants that ranged from 0 to 30%. While in the monokaryon strain W1-26, only the millet medium group obtained positive transformants with a positive frequency of 75.48%. Moreover, three dikaryotic wild strains (YS55, YS3334, and YS3357) and two dikaryotic cultivated strains (W1 and S606) showed the highest transformation efficiency, with 32.96% of the germination frequency, and 85.12% of positive frequency for hygromycin-resistant transformants. This work demonstrated that Agrobacterium-mediated transformation was successfully performed in L. edodes, and the genotype of recipients as well as the medium for mycelial incubation were suggested to play key roles in determining the transformation efficiency. These findings may provide new avenues for the genetic modification of edible mushroom and may extend the cognition of DNA-mediated transformation in filamentous fungi.

Transcriptome analysis reveals the role of nitric oxide in Pleurotus eryngii responses to Cd2+ stress

Pleurotus eryngii is widely cultivated in China. However, our understanding of its transcriptional response to heavy metal stress and the underlying mechanism of nitric oxide (NO) in enhancing its tolerance to heavy metals is limited. In the present study, RNA-seq was used to generate large transcript sequences from P. eryngii exposed to cadmium chloride (CdCl₂) and exogenous NO. A total of 45,833 unigenes were assembled from the P. eryngii transcriptome, of which 32,333 (70.54%) unigenes matched known proteins in the nr database. Transcriptional analysis revealed that putative genes encoding heat shock proteins (HSPs) and genes participating in glycerolipid metabolism and steroid biosynthesis were significantly up-regulated in P. eryngii exposed to 50 μM Cd (P < 0.05). P. eryngii mycelia exposed to extremely high levels of heavy metals showed an increase in biomass when exogenous NO was added to the culture. The collaboration of putative oxidoreductase, dehydrogenase, reductase, transferase genes and transcription factors such as "GTPase activator activity", "transcription factor complex", "ATP binding", "GTP binding", and "enzyme activator activity", which were significantly up-regulated in samples induced by exogenous NO, contributed to the enhancement of P. eryngii tolerance to extremely high levels of heavy metals. The study provides a new insight into the transcriptional response of P. eryngii to extremely high levels of heavy metals and the mechanism of NO in enhancing heavy metal tolerance.