Genome assembly of M. spongiola and comparative genomics of the genus Morchella provide initial insights into taxonomy and adaptive evolution.

Morchella spongiola is a highly prized mushroom for its delicious flavor and medical value and is one of the most flourishing, representative, and dominant macrofungi in the Qilian Mountains of the Qinghai-Tibet Plateau subkingdoms (QTPs). However, the understanding of M. spongiola remains largely unknown, and its taxonomy is ambiguous. In this study, we redescribed a unique species of M. spongiola, i.e., micromorphology, molecular data, genomics, and comparative genomics, and the historical biogeography of M. spongiola were estimated for 182 single-copy homologous genes. A high-quality chromosome-level reference genome of M. spongiola M12-10 was obtained by combining PacBio HiFi data and Illumina sequencing technologies; it was approximately 57.1 Mb (contig N50 of 18.14 Mb) and contained 9775 protein-coding genes. Comparative genome analysis revealed considerable conservation and unique characteristics between M. spongiola M12-10 and 32 other Morchella species. Molecular phylogenetic analysis indicated that M. spongiola M12-10 is similar to the M. prava/Mes-7 present in sandy soil near rivers, differentiating from black morels ~ 43.06 Mya (million years ago), and diverged from M. parva/Mes-7 at approximately 12.85 Mya (in the Miocene epoch), which is closely related to the geological activities in the QTPs (in the Neogene). Therefore, M. spongiola is a unique species rather than a synonym of M. vulgaris/Mes-5, which has a distinctive grey-brown sponge-like ascomata. This genome of M. spongiola M12-10 is the first published genome sequence of the species in the genus Morchella from the QTPs, which could aid future studies on functional gene identification, germplasm resource management, and molecular breeding efforts, as well as evolutionary studies on the Morchella taxon in the QTPs.

Distribution of rhizosphere fungi of Kobresia humilis on the Qinghai-Tibet Plateau.

Kobresia humilis is a major species in the alpine meadow communities of the Qinghai-Tibet Plateau (QTP); it plays a crucial role in maintaining the ecological balance of these meadows. Nevertheless, little is known about the rhizosphere fungi associated with K. humilis on the Qinghai Tibet Plateau. In this study, we used Illumina Miseq to investigate the fungal diversity, community structure, and ecological types in the root and rhizosphere soil of K. humilis across eight areas on the QTP and analyzed the correlation between rhizosphere fungi of K. humilis and environmental factors. A total of 19,423 and 25,101 operational taxonomic units (OTUs) were obtained from the roots and rhizosphere soil of K. humilis. These were classified into seven phyla, 25 classes, 68 orders, 138 families, and 316 genera in the roots, and nine phyla, 31 classes, 76 orders, 152 families, and 407 genera in the rhizosphere soil. There were 435 and 415 core OTUs identified in root and rhizosphere soil, respectively, which were categorized into 68 and 59 genera, respectively, with 25 shared genera. Among them, the genera with a relative abundance >1% included Mortierella, Microscypha, Floccularia, Cistella, Gibberella, and Pilidium. Compared with the rhizosphere soil, the roots showed five differing fungal community characteristics, as well as differences in ecological type, and in the main influencing environmental factors. First, the diversity, abundance, and total number of OTUs in the rhizosphere soil of K. humilis were higher than for the endophytic fungi in the roots by 11.85%, 9.85%, and 22.62%, respectively. The composition and diversity of fungal communities also differed between the eight areas. Second, although saprotroph-symbiotrophs were the main ecological types in both roots and rhizosphere soil; there were 62.62% fewer pathotrophs in roots compared to the rhizosphere soil. Thirdly, at the higher altitude sites (3,900-4,410 m), the proportion of pathotroph fungi in K. humilis was found to be lower than at the lower altitude sites (3,200-3,690 m). Fourthly, metacommunity-scale network analysis showed that during the long-term evolutionary process, ZK (EICZK = 1) and HY (EICHY = 1) were critical sites for development of the fungal community structure in the roots and rhizosphere soil of K. humilis, respectively. Fifthly, canonical correspondence analysis (CCA) showed that key driving factors in relation to the fungal community were longitude (R2 = 0.5410) for the root community and pH (R2 = 0.5226) for the rhizosphere soil community. In summary, these results show that K. humilis fungal communities are significantly different in the root and rhizosphere soil and at the eight areas investigated, indicating that roots select for specific microorganisms in the soil. This is the first time that the fungal distribution of K. humilis on the QTP in relation to long-term evolutionary processes has been investigated. These findings are critical for determining the effects of environmental variables on K. humilis fungal communities and could be valuable when developing guidance for ecological restoration and sustainable utilization of the biological resources of the QTP.

Corrigendum: Out of the Qinghai-Tibetan plateau: origin, evolution and historical biogeography of Morchella (both Elata and Esculenta clades).

[This corrects the article DOI: 10.3389/fmicb.2022.1078663.].

The Molecular Mechanism of Yellow Mushroom (Floccularia luteovirens) Response to Strong Ultraviolet Radiation on the Qinghai-Tibet Plateau.

The Qinghai-Tibet Plateau (QTP) is the highest plateau in the world, and its ultraviolet (UV) radiation is much greater than that of other regions in the world. Yellow mushroom (Floccularia luteovirens) is a unique and widely distributed edible fungus on the QTP. However, the molecular mechanism of F. luteovirens's response to strong UV radiation remains unclear. Herein, we reported the 205 environmental adaptation and information processing genes from genome of F. luteovirens. In addition, we assembled the RNA sequence of UV-affected F. luteovirens at different growth stages. The results showed that in response to strong UV radiation, a total of 11,871 significantly different genes were identified, of which 4,444 genes in the vegetative mycelium (VM) stage were significantly different from the young fruiting bodies (YFB) stage, and only 2,431 genes in the YFB stage were significantly different from fruiting bodies (FB) stage. A total of 225 differentially expressed genes (DEGs) were found to be involved in environmental signal transduction, biochemical reaction preparation and stress response pathway, pigment metabolism pathway, and growth cycle regulation, so as to sense UV radiation, promote repair damage, regulate intracellular homeostasis, and reduce oxidative damage of UV radiation. On the basis of these results, a molecular regulation model was proposed for the response of F. luteovirens to strong UV radiation. These results revealed the molecular mechanism of adaptation of F. luteovirens adapting to strong UV radiation, and provided novel insights into mechanisms of fungi adapting to extreme environmental conditions on the QTP; the production the riboflavin pigment of the endemic fungi (Yellow mushroom) in the QTP was one of the response to extreme environment of the strong UV radiation.

Microbial response to multiple-level addition of grass organic matter in lake sediments with different salinity.

Water surface expansion of saline lakes usually causes the inundation of surrounding grassland, leading to the increase of terrestrial grass organic matter (OM) input to the lakes and the decrease of lake salinity. However, the influence of terrestrial grass OM input increase and salinity decrease on organic carbon mineralization and microbial community composition remains unknown in saline lakes. Here, microbial mineralization of terrestrial grass (Achnatherum splendens) OM at different quantity levels in lake sediments with different salinity was investigated by performing microcosm experiments. The results showed that the CO2 production rates increased with the increase of grass OM supply in the studied sediments with different salinity, which may be driven by certain microbial groups (e.g. Bacteroidota, Firmicutes, and Ascomycota). The increase of grass OM supply reduced the richness of prokaryotic community, which will decrease the size and complexity of the studied microbial networks, but increase the interaction between prokaryotic and fungal taxa. Taken together, our results suggest that the increase of terrestrial grass OM input caused by lake expansion would enhance the mineralization of organic carbon and affect the community composition and interactions of related microorganisms in lake sediments with different salinity.

Soil metabolomics and bacterial functional traits revealed the responses of rhizosphere soil bacterial community to long-term continuous cropping of Tibetan barley.

Continuous cropping often leads to an unbalanced soil microbial community, which in turn negatively affects soil functions. However, systematic research of how these effects impact the bacterial composition, microbial functional traits, and soil metabolites is lacking. In the present study, the rhizosphere soil samples of Tibetan barley continuously monocropped for 2 (CCY02), 5 (CCY05), and 10 (CCY10) years were collected. By utilizing 16S high-throughput sequencing, untargeted metabolomes, and quantitative microbial element cycling smart chips, we examined the bacterial community structure, soil metabolites, and bacterial functional gene abundances, respectively. We found that bacterial richness (based on Chao1 and Phylogenetic Diversity [PD] indices) was significantly higher in CCY02 and CCY10 than in CCY05. As per principal component analysis (PCA), samples from the continuous monocropping year tended to share more similar species compositions and soil metabolites, and exhibited distinct patterns over time. The results of the Procrustes analysis indicated that alterations in the soil metabolic profiles and bacterial functional genes after long-term continuous cropping were mainly mediated by soil microbial communities (P < 0.05). Moreover, 14 genera mainly contributed to the sample dissimilarities. Of these, five genera were identified as the dominant shared taxa, including Blastococcus, Nocardioides, Sphingomonas, Bacillus, and Solirubrobacter. The continuous cropping of Tibetan barley significantly increased the abundances of genes related to C-degradation (F = 9.25, P = 0.01) and P-cycling (F = 5.35, P = 0.03). N-cycling significantly negatively correlated with bacterial diversity (r =  - 0.71, P = 0.01). The co-occurrence network analysis revealed that nine hub genera correlated with most of the functional genes and a hub taxon, Desulfuromonadales, mainly co-occurred with the metabolites via both negative and positive correlations. Collectively, our findings indicated that continuous cropping significantly altered the bacterial community structure, functioning of rhizosphere soils, and soil metabolites, thereby providing a comprehensive understanding of the effects of the long-term continuous cropping of Tibetan barley.

Nitrite- and N2O-reducing bacteria respond differently to ecological factors in saline lakes.

The distribution of nitrite- and N2O-reducing bacteria is key to potential N2O emission from lakes. However, such information in highland saline lakes remains unknown. Here, we investigated the abundance and community composition of nitrite- and N2O-reducing bacteria in the sediments of six saline lakes on the Qing-Tibetan Plateau. The studied lakes covered a wide range of salinity (1.0-340.0 g/L). Results showed that in the studied saline lake sediments, nitrite-reducing bacteria were significantly more abundant than N2O-reducing bacteria, and their abundances ranged 7.14 × 103-8.26 × 108 and 1.18 × 106-6.51 × 107 copies per gram sediment (dry weight), respectively. Nitrite-reducing bacteria were mainly affiliated with α-, ß- and γ-Proteobacteria, with ß- and α-Proteobacteria being dominant in low- and high-salinity lakes, respectively; N2O-reducing bacterial communities mainly consisted of Proteobacteria (α-, ß-, γ- and δ-subgroups), Bacteroidetes, Verrucomicrobia, Actinobacteria, Chloroflexi, Gemmatimonadetes and Balneolaeota, with Proteobacteria and Bacteroidetes/Verrucomicrobia dominating in low- and high-salinity lakes, respectively. The nitrite- and N2O-reducing bacterial communities showed distinct responses to ecological factors, and they were mainly regulated by mineralogical and physicochemical factors, respectively. In response to salinity change, the community composition of nitrite-reducing bacteria was more stable than that of N2O-reducing bacteria. These findings suggest that nitrite- and N2O-reducing bacteria may prefer niches with different salinity.

Out of the Qinghai-Tibetan plateau: Origin, evolution and historical biogeography of Morchella (both Elata and Esculenta clades).

Introduction: Morchella has become a research hotspot because of its wide distribution, delicious taste, and phenotypic plasticity. The Qinghai-Tibet Plateau subkingdoms (QTPs) are known as the cradle of Ice age biodiversity. However, the diversity of Morchella in the QTPs has been poorly investigated, especially in phylogenetic diversity, origin, and biogeography. Methods: The genealogical concordance phylogenetic species recognition (GCPSR, based on Bayesian evolutionary analysis using sequences from the internal transcribed spacer (ITS), nuclear large subunit rDNA (nrLSU), translation elongation factor 1-α (EF1-α), and the largest and second largest subunits of RNA polymerase II (RPB1 and RPB2)), differentiation time estimation, and ancestral region reconstruction were used to infer Morchella's phylogenetic relationships and historical biogeography in the QTPs. Results: Firstly, a total of 18 Morchella phylogenetic species are recognized in the QTPs, including 10 Elata clades and 8 Esculenta clades of 216 individuals Secondly, the divergences of the 18 phylogenetic species were 50.24-4.20 Mya (Eocene-Pliocene), which was closely related to the geological activities in the QTPs. Furthermore, the ancestor of Morchella probably originated in the Northern regions (Qilian Shan, Elata cade) and southwestern regions (Shangri-La, Esculenta clade) of QTPs and might have migrated from North America (Rufobrunnea clade) via Beringian Land Bridge (BLB) and Long-Distance Dispersal (LDD) expansions during the Late Cretaceous. Moreover, as the cradle of species origin and diversity, the fungi species in the QTPs have spread out and diffused to Eurasia and South Africa starting in the Paleogene Period. Conclusion: This is the first report that Esculenta and Elata clade of Morchella originated from the QTPs because of orogenic, and rapid differentiation of fungi is strongly linked to geological uplift movement and refuge in marginal areas of the QTPs. Our findings contribute to increasing the diversity of Morchella and offer more evidence for the origin theory of the QTPs.

Morphological changes and bioaccumulation in response to cadmium exposure in Morchella spongiola, a fungus with potential for detoxification.

Morchella is a genus of edible fungi with strong resistance to cadmium (Cd) and the ability to accumulate it in its mycelia. However, the mechanisms underlying Cd resistance in Morchella remain unknown. In the present study, morphological and physiological responses to Cd were evaluated in the mycelia of Morchella spongiola. Variations in hyphal micromorphology, including twisting, folding, and kinking, in mycelia exposed to different Cd concentrations (0.15, 0.9, 1.5, 2.4, 5.0 mg/L) were observed using scanning electron microscopy. Deposition of Cd precipitates on cell surfaces (at Cd concentrations >2.4 mg/L) was shown by SEM-EDS. Transmission electron microscopy analysis of cells exposed to different concentrations of Cd revealed the loss of intracellular structures and the localization of Cd depositions inside and outside the cell. FTIR analysis showed that functional groups such as C=O, -OH, -NH, and -CH could be responsible for Cd binding on the cell surface of M. spongiola. In addition, intracellular accumulation was observed in cultures at low Cd concentrations (<0.9 mg/L), while extracellular adsorption occurred at higher concentrations. These results provide valuable information on the Cd tolerance mechanism in M. spongiola and provide a robust foundation for further studies on fungal bioremediation strategies.

Transcriptome Analysis and Expression Profiling of Molecular Responses to Cd Toxicity in Morchella spongiola

Morchella is a genus of fungi with the ability to concentrate Cd both in the fruit-body and mycelium. However, the molecular mechanisms conferring resistance to Cd stress in Morchella are unknown. Here, RNA-based transcriptomic sequencing was used to identify the genes and pathways involved in Cd tolerance in Morchella spongiola. 7444 differentially expressed genes (DEGs) were identified by cultivating M. spongiola in media containing 0.15, 0.90, or 1.50 mg/L Cd2+ . The DEGs were divided into six sub-clusters based on their global expression profiles. GO enrichment analysis indicated that numerous DEGs were associated with catalytic activity, cell cycle control, and the ribosome. KEGG enrichment analysis showed that the main pathways under Cd stress were MAPK signaling, oxidative phosphorylation, pyruvate metabolism, and propanoate metabolism. In addition, several DEGs encoding ion transporters, enzymaticon-enzymatic antioxidants, and transcription factors were identified. Based on these results, a preliminary gene regulatory network was firstly proposed to illustrate the molecular mechanisms of Cd detoxification in M. spongiola. These results provide valuable insights into the Cd tolerance mechanism of M. spongiola and constitute a robust foundation for further studies on detoxification mechanisms in macrofungi that could potentially lead to the development of new and improved fungal bioremediation strategies.

Transcriptome Analysis and Expression Profiling of Molecular Responses to Cd Toxicity in Morchella spongiola

Morchella is a genus of fungi with the ability to concentrate Cd both in the fruit-body and mycelium. However, the molecular mechanisms conferring resistance to Cd stress in Morchella are unknown. Here, RNA-based transcriptomic sequencing was used to identify the genes and pathways involved in Cd tolerance in Morchella spongiola. 7444 differentially expressed genes (DEGs) were identified by cultivating M. spongiola in media containing 0.15, 0.90, or 1.50 mg/L Cd2+ . The DEGs were divided into six sub-clusters based on their global expression profiles. GO enrichment analysis indicated that numerous DEGs were associated with catalytic activity, cell cycle control, and the ribosome. KEGG enrichment analysis showed that the main pathways under Cd stress were MAPK signaling, oxidative phosphorylation, pyruvate metabolism, and propanoate metabolism. In addition, several DEGs encoding ion transporters, enzymaticon-enzymatic antioxidants, and transcription factors were identified. Based on these results, a preliminary gene regulatory network was firstly proposed to illustrate the molecular mechanisms of Cd detoxification in M. spongiola. These results provide valuable insights into the Cd tolerance mechanism of M. spongiola and constitute a robust foundation for further studies on detoxification mechanisms in macrofungi that could potentially lead to the development of new and improved fungal bioremediation strategies.

Minerals play key roles in driving prokaryotic and fungal communities in the surface sediments of the Qinghai-Tibetan lakes.

There is limited knowledge of the relative influences of deterministic and stochastic processes on prokaryotic and fungal communities in lake sediments. In this study, we surveyed the prokaryotic and fungal community compositions and their influencing factors in 23 surface sediments from six lakes on the Qinghai-Tibetan Plateau (QTP) with the use of Illumina sequencing. The results showed the distribution of prokaryotic and fungal communities in the studied QTP lake sediments was shaped by different assembly processes, with prokaryotes primarily governed by variable selection and homogenizing dispersal (accounting for 57.9% and 37.3% of the observed variations) and fungi being mainly regulated by variable selection, non-dominant processes and homogenizing dispersal (38.3%, 43.7% and 13.7%, respectively). Regarding the variable selection, mineralogical variables played key roles in shaping prokaryotic and fungal community structures. Collectively, these findings expand current knowledge concerning the influences of deterministic (e.g. variable selection) and stochastic processes (e.g. homogenizing dispersal and non-dominant processes) on the prokaryotic and fungal distribution in the QTP lakes.

Surviving onshore soil microbial communities differ among the Qing-Tibetan lakes with different salinity.

Little is known about the onshore microbial contribution to the microbial communities in nearby lakes and its response to salinity. In this study, transplanting experiments were established by caging onshore soils with dialysis bags followed by in situ 50-day incubation in nearby lakes with different salinity on the Qinghai-Tibetan Plateau. At the end of the experiment, geochemical and microbial analyses were performed on the original soils, caged soils and lake waters and sediments at the incubation sites. The results showed that the salinity increased significantly (P < 0.05) in the caged soils and such salinity increases showed significant (P < 0.05) positive correlation with the salinity of the studied lakes. The microbial community composition and predicted functions in the caged soils were significantly (P < 0.05) changed in comparison with their corresponding original soils, and such variation could be mainly explained by the succession of members of the Proteobacteria, Bacteroidetes and Actinobacteria from the original soils to their corresponding caged soils. The onshore microbial contribution appeared to be limited (up to 11.2% for sediment and negligible for water, respectively) to nearby lake microbial communities. Nevertheless, the survival of onshore soil microbial communities was mainly limited by the salinity of the receiving lakes.

[Purification and characterization of endoglucanase Egn21 from Fusarium sp. Q7-31T].

Objective: The objective of this research was to study plant cell wall degradation enzymes from Fusarium sp. Q7-31T. Methods: Strain was cultured in liquid medium with 1% (W/V) peptone as nitrogen source, 0.5% (W/V) oat straw as carbon source, 120 r/min shaking at 20 °C for 3 days. The endoglucanase Egn21 was purified by using Sephacry S-100 chromatography and DEAE-sepharose ion-exchange column chromatography. Then the enzymatic properties and MADIL-TOF-TOF identification were analyzed. Results: The molecular weight and isoelectric point (pI) of Egn21 was 44.25 kDa and 4.91, respectively. Egn21 had optimal activity with carboxymethyl cellulose at 40 °C and pH 6.0, stable at 45 °C and pH between 5.0 and 8.0, inhibited by Fe2+, Ca2+, K+, Na+, Mn2+ and inactivated by Hg2+, whereas Co2+, Zn2+ and Mg2+ had no effect. Conclusion: The enzymatic properties and MADIL-TOF-TOF results suggested that Egn21 belongs to GH5 family.

[Purification, identification and characterization of an endoglucanase Egn20 from Fusarium sp. Q7-31T].

OBJECTIVE: The endoglucanase from Fusarium sp. Q7-31T was isolated, purified, identified and characterized to provide data for enzyme system of Fusarium sp. . [Methods] Strain was cultured in liquid fermentation with oat straw as carbon source, the endoglucanase was purified by using Sephacry S-100 chromatography and DEAE-sepharose ion-exchange column chromatography and the enzymatic properties were studied. The protein was identified using MADIL-TOF-TOF. RESULTS: An endoglucanase was purified and named Egn20. The molecular weight was 55.37 kDa and isoelectric point (pI) was 7.44. Egn20 had optimal activity with carboxymethyl cellulose at 40 degrees C and pH 6.0, stabilized at 45 degrees C and pH 5.0 - 7.0, activated by Fe2+, inhibited by Na+, Ca2+, Mg2+, Zn2+, K+ and inactivated by Hg2+. The enzymatic properties and MADIL-TOF-TOF results suggested that Egn20 belongs to GH7 family. CONCLUSION: Our results may provide important data for the study of Fusarium sp. enzyme system.

Cloning of a novel xylanase gene from a newly isolated Fusarium sp. Q7-31 and its expression in Escherichia coli

A strain of Q7-31 was isolated from Qinghai-Tibet Plateau and was identified as Fusarium sp. based on its morphological characteristics and ITS rDNA gene sequence analysis. It has the highest capacity of degrading cell wall activity compared with other 11 strains. To do research on its xylanase activity of Fusarium sp. Q7-31 while the degrading the rice cell walls, the complete gene xyn8 that encodes endo-1, 4-¥â-xylanase secreted by Fusarium sp. Q7-31 was cloned and sequenced. The coding region of the gene is separated by two introns of 56bp and 55bp. It encodes 230 amino acid residues of a protein with a calculated molecular weight of 25.7 kDa. The animo acids sequence of xyn8 gene has higher similarity with those of family 11 of glycosyl hydrolases reported from other microorganisms. The nature peptide encodeing cDNA was subcloned into pGEX5x-1 expression vector. The recombinant plasmid was expressed in Escherichia coli BL21-CodonPlus (DE3)-RIL, and xylanase activity was measured. The expression fusion protein was identified by SDS-PAGE and Western blotting, a new specific band of about 52kDa was identified when induced by IPTG. Enzyme activity assay verified the recombinants proteins as a xylanase. A maxium activity of 2.34U/ mg, the xylanase had optimal activity at pH 6.0 and temperature 40¨¬C .

Cloning of a novel xylanase gene from a newly isolated Fusarium sp. Q7-31 and its expression in Escherichia coli.

A strain of Q7-31 was isolated from Qinghai-Tibet Plateau and was identified as Fusarium sp. based on its morphological characteristics and ITS rDNA gene sequence analysis. It has the highest capacity of degrading cell wall activity compared with other 11 strains. To do research on its xylanase activity of Fusarium sp. Q7-31 while the degrading the rice cell walls, the complete gene xyn8 that encodes endo-1, 4-ß-xylanase secreted by Fusarium sp. Q7-31 was cloned and sequenced. The coding region of the gene is separated by two introns of 56bp and 55bp. It encodes 230 amino acid residues of a protein with a calculated molecular weight of 25.7 kDa. The animo acids sequence of xyn8 gene has higher similarity with those of family 11 of glycosyl hydrolases reported from other microorganisms. The nature peptide encodeing cDNA was subcloned into pGEX5x-1 expression vector. The recombinant plasmid was expressed in Escherichia coli BL21-CodonPlus (DE3)-RIL, and xylanase activity was measured. The expression fusion protein was identified by SDS-PAGE and Western blotting, a new specific band of about 52kDa was identified when induced by IPTG. Enzyme activity assay verified the recombinants proteins as a xylanase. A maxium activity of 2.34U/ mg, the xylanase had optimal activity at pH 6.0 and temperature 40℃.